The within-host population dynamics of Mycobacterium tuberculosis vary with treatment efficacy

Advantages of updated WHO mutation catalog combined with existing whole-genome sequencing-based approaches for Mycobacterium tuberculosis resistance prediction

BACKGROUND

The WHO recently released a second edition of the mutation catalog for predicting drug resistance in Mycobacterium tuberculosis (MTB). This study evaluated its effectiveness compared to existing whole-genome sequencing (WGS)-based prediction methods and proposes a novel approach for its optimization.

METHODS

We tested the accuracy of five tools-the WHO catalog, TB Profiler, SAM-TB, GenTB, and MD-CNN-for predicting drug susceptibility on a global dataset of 36,385 MTB isolates with high-quality phenotypic drug susceptibility testing (DST) and WGS data. By integrating the genotypic DST predictions of these five tools in an ensemble machine learning framework, we developed an improved computational model for MTB drug susceptibility prediction. We then validated the ensemble model on 860 MTB isolates with phenotypic and WGS data collected in Shenzhen, China (2013-2019) and Valencia, Spain (2014-2016).

RESULTS

Among the five genotypic DST tools for predicting susceptibility to ten drugs, MD-CNN exhibited the highest overall performance (AUC 92.1%; 95% CI 89.8-94.4%). The WHO catalog demonstrated the highest specificity of 97.3% (95% CI 95.8-98.4%), while TB Profiler had the best sensitivity at 79.5% (95% CI 71.8-86.2%). The ensemble machine learning model (AUC 93.4%; 95% CI 91.4-95.4%) outperformed all of the five individual tools, with a specificity of 95.4% (95% CI 93.0-97.6%) and a sensitivity of 84.1% (95% CI 78.8-88.8%), principally due to considerable improvements in second-line drug resistance predictions (AUC 91.8%; 95% CI 89.6-94.0%).

CONCLUSIONS

The second edition of the WHO MTB mutation catalog does not, by itself, perform better than existing tools for predicting MTB drug resistance. An integrative approach combining the WHO catalog with other genotypic DST methods significantly enhances prediction accuracy.

Ecology, global diversity and evolutionary mechanisms in the Mycobacterium tuberculosis complex

With the COVID-19 pandemic receding, tuberculosis (TB) is again the number one cause of human death to a single infectious agent. TB is caused by bacteria that belong to the Mycobacterium tuberculosis complex (MTBC). Recent advances in genome sequencing have provided new insights into the ecology and evolution of the MTBC. This includes the discovery of new phylogenetic lineages within the MTBC, a deeper understanding of the host tropism among the various animal-adapted lineages, enhanced knowledge on the evolutionary dynamics of antimicrobial resistance and transmission, as well as a better grasp of the within-host MTBC diversity. Moreover, advances in long-read sequencing are increasingly highlighting the relevance of structural genomic variation in the MTBC. These findings not only shed new light on the biology and epidemiology of TB, but also give rise to new questions and research avenues. The purpose of this Review is to summarize these new insights and discuss their implications for global TB control.

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.

Molecular characteristics of drug-susceptible Mycobacterium tuberculosis clinical isolates based on treatment duration

BACKGROUND

In this study, we performed comparative genomic and transcriptomic analysis of clinical isolates of Mycobacterium tuberculosis collected from patients with drug-susceptible tuberculosis (DS-TB). The clinical isolates were categorized based on treatment duration: standard 6 months or >6 months.

METHODS

Study participants were recruited from a 2016 to 2018 tuberculosis cohort, and clinical M. tuberculosis isolates were collected from the sputum of patients with tuberculosis. We analyzed the genome and transcriptome of the isolated M. tuberculosis.

RESULTS

Genomic analysis revealed a specific non-synonymous single-nucleotide polymorphism in pe_pgrs9 and ppe34, exclusive to the group treated for >6 months. Transcriptomic analysis revealed increased expression of various virulence-associated protein family genes and decreased expression of ribosomal protein genes and ppe38 genes in the group treated for >6 months.

CONCLUSION

The identified genetic variation and gene expression patterns may influence treatment outcomes by modulating host immune responses, increasing virulence, and potentially contributing to persister cell formation in M. tuberculosis. This study provides insights into the genetic and transcriptomic factors associated with prolonged DS-TB treatment. However, our study identified molecular characteristics using a small sample size, and further detailed studies are warranted.

Rifampicin tolerance and growth fitness among isoniazid-resistant clinical Mycobacterium tuberculosis isolates from a longitudinal study

Antibiotic tolerance in Mycobacterium tuberculosis reduces bacterial killing, worsens treatment outcomes, and contributes to resistance. We studied rifampicin tolerance in isolates with or without isoniazid resistance (IR). Using a minimum duration of killing assay, we measured rifampicin survival in isoniazid-susceptible (IS, n=119) and resistant (IR, n=84) isolates, correlating tolerance with bacterial growth, rifampicin minimum inhibitory concentrations (MICs), and isoniazid-resistant mutations. Longitudinal IR isolates were analyzed for changes in rifampicin tolerance and genetic variant emergence. The median time for rifampicin to reduce the bacterial population by 90% (MDK90) increased from 1.23 days (IS) and 1.31 days (IR) to 2.55 days (IS) and 1.98 days (IR) over 15-60 days of incubation, indicating fast and slow-growing tolerant sub-populations. A 6 log10-fold survival fraction classified tolerance as low, medium, or high, showing that IR is linked to increased tolerance and faster growth (OR = 2.68 for low vs. medium, OR = 4.42 for low vs. high, p-trend = 0.0003). High tolerance in IR isolates was associated with rifampicin treatment in patients and genetic microvariants. These findings suggest that IR tuberculosis should be assessed for high rifampicin tolerance to optimize treatment and prevent the development of multi-drug-resistant tuberculosis.

An additional proofreader contributes to DNA replication fidelity in mycobacteria

The removal of mis-incorporated nucleotides by proofreading activity ensures DNA replication fidelity. Whereas the ε-exonuclease DnaQ is a well-established proofreader in the model organism Escherichia coli, it has been shown that proofreading in a majority of bacteria relies on the polymerase and histidinol phosphatase (PHP) domain of replicative polymerase, despite the presence of a DnaQ homolog that is structurally and functionally distinct from E. coli DnaQ. However, the biological functions of this type of noncanonical DnaQ remain unclear. Here, we provide independent evidence that noncanonical DnaQ functions as an additional proofreader for mycobacteria. Using the mutation accumulation assay in combination with whole-genome sequencing, we showed that depletion of DnaQ in Mycolicibacterium smegmatis leads to an increased mutation rate, resulting in AT-biased mutagenesis and increased insertions/deletions in the homopolymer tract. Our results showed that mycobacterial DnaQ binds to the β clamp and functions synergistically with the PHP domain proofreader to correct replication errors. Furthermore, the loss of dnaQ results in replication fork dysfunction, leading to attenuated growth and increased mutagenesis on subinhibitory fluoroquinolones potentially due to increased vulnerability to fork collapse. By analyzing the sequence polymorphism of dnaQ in clinical isolates of Mycobacterium tuberculosis (Mtb), we demonstrated that a naturally evolved DnaQ variant prevalent in Mtb lineage 4.3 may enable hypermutability and is associated with drug resistance. These results establish a coproofreading model and suggest a division of labor between DnaQ and PHP domain proofreader. This study also provides real-world evidence that a mutator-driven evolutionary pathway may exist during the adaptation of Mtb.

Mycobacterium tuberculosis is less likely to acquire pathogenic mutations during latent infection than during active disease

Most people infected with Mycobacterium tuberculosis (Mtb) are believed to be in a state of latent tuberculosis (TB) infection (LTBI). Although LTBI is asymptomatic and not infectious, there is a risk of developing active disease even decades after infection. Here, to characterize mutations acquired during LTBI, we collected and analyzed Mtb genomes from seven Japanese patient pairs, each pair consisting of two active TB patients whose starting dates of developing active disease were >3 years apart; one had a high suspicion of LTBI before developing active disease, whereas the other did not. Thereafter, we compared these genomes with those of longitudinal sample pairs within a host of chronic active TB infections combined with public data. The bacterial populations in patients with LTBI were genetically more homogeneous and accumulated single nucleotide polymorphisms (SNPs) slower than those from active disease. Moreover, the lower proportion of nonsynonymous SNPs indicated weaker selective pressures during LTBI than active disease. Finally, the different mutation spectrums indicated different mutators between LTBI and active disease. These results suggest that the likelihood of the acquisition of mutations responsible for antibiotic resistance and increased virulence was lower in the Mtb population from LTBI than active disease.IMPORTANCEControlling latent tuberculosis (TB) infection (LTBI) activation is an effective strategy for TB elimination, where understanding Mycobacterium tuberculosis (Mtb) dynamics within the host plays an important role. Previous studies on chronic active disease reported that Mtb accumulated genomic mutations within the host, possibly resulting in acquired drug resistance and increased virulence. However, several reports suggest that fewer mutations accumulate during LTBI than during the active disease, but the associated risk is largely unknown. Here, we analyzed the genomic dynamics of Mtb within the host during LTBI. Our results statistically suggest that Mtb accumulates mutations during LTBI, but most mutations are under low selective pressures, which induce mutations responsible for drug resistance and virulence. Thus, we propose that LTBI acts as a source for new TB disease rather than as a period for in-host genome evolution.

Rifampicin tolerance and growth fitness among isoniazid-resistant clinical Mycobacterium tuberculosis isolates: an in-vitro longitudinal study

Antibiotic tolerance in Mycobacterium tuberculosis leads to less effective bacterial killing, poor treatment responses and resistant emergence. Therefore, we investigated the rifampicin tolerance of M. tuberculosis isolates, with or without pre-existing isoniazid-resistance. We determined the in-vitro rifampicin survival fraction by minimum duration of killing assay in isoniazid susceptible (IS, n=119) and resistant (IR, n=84) M. tuberculosis isolates. Then we correlated the rifampicin tolerance with bacterial growth, rifampicin minimum inhibitory concentrations (MICs) and isoniazid-resistant mutations. The longitudinal IR isolates collected from patients were analyzed for changes in rifampicin tolerance and associated emergence of genetic variants. The median duration of rifampicin exposure reducing the M. tuberculosis surviving fraction by 90% (minimum duration of killing-MDK90) increased from 1.23 (95%CI 1.11; 1.37) and 1.31 (95%CI 1.14; 1.48) to 2.55 (95%CI 2.04; 2.97) and 1.98 (95%CI 1.69; 2.56) days, for IS and IR respectively, during 15 to 60 days of incubation. This indicated the presence of fast and slow growing tolerant sub-populations. A range of 6 log 10 -fold survival fraction enabled classification of tolerance as low, medium or high and revealed IR association with increased tolerance with faster growth (OR=2.68 for low vs. medium, OR=4.42 for low vs. high, P -trend=0.0003). The high tolerance in IR isolates was specific to those collected during rifampicin treatment in patients and associated with bacterial genetic microvariants. Furthermore, the high rifampicin tolerant IR isolates have survival potential similar to multi-drug resistant isolates. These findings suggest that IR tuberculosis needs to be evaluated for high rifampicin tolerance to improve treatment regimen and prevent the risk of MDR-TB emergence.

Genetic factors associated with acquired phenotypic drug resistance and its compensatory evolution during tuberculosis treatment

OBJECTIVES

We elucidated the factors, evolution, and compensation of antimicrobial resistance (AMR) in Mycobacterium tuberculosis (MTB) isolates under dual pressure from the intra-host environment and anti-tuberculosis (anti-TB) drugs.

METHODS

This retrospective case-control study included 337 patients with pulmonary tuberculosis from 15 clinics in Tianjin, China, with phenotypic drug susceptibility testing results available for at least two time points between January 1, 2009 and December 31, 2016. Patients in the case group exhibited acquired AMR to isoniazid (INH) or rifampicin (RIF), while those in the control group lacked acquired AMR. The whole-genome sequencing (WGS) was conducted on 149 serial longitudinal MTB isolates from 46 patients who acquired or reversed phenotypic INH/RIF-resistance during treatment. The genetic basis, associated factors, and intra-host evolution of acquired phenotypic INH/RIF-resistance were elucidated using a combined analysis.

RESULTS

Anti-TB interruption duration of ≥30 days showed association with acquired phenotypic INH/RIF resistance (aOR = 2·2, 95% CI, 1·0-5·1) and new rpoB mutations (p = 0·024). The MTB evolution was 1·2 (95% CI, 1·02-1·38) single nucleotide polymorphisms per genome per year under dual pressure from the intra-host environment and anti-TB drugs. AMR-associated mutations occurred before phenotypic AMR appearance in cases with acquired phenotypic INH (10 of 16) and RIF (9 of 22) resistances.

DISCUSSION

Compensatory evolution may promote the fixation of INH/RIF-resistance mutations and affect phenotypic AMR. The TB treatment should be adjusted based on gene sequencing results, especially in persistent culture positivity during treatment, which highlights the clinical importance of WGS in identifying reinfection and AMR acquisition before phenotypic drug susceptibility testing.

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.

Mycobacterial DnaQ is an Alternative Proofreader Ensuring DNA Replication Fidelity

Remove of mis-incorporated nucleotides ensures replicative fidelity. Although the ε-exonuclease DnaQ is a well-established proofreader in the model organism Escherichia coli, proofreading in mycobacteria relies on the polymerase and histidinol phosphatase (PHP) domain of replicative polymerase despite the presence of an alternative DnaQ homolog. Here, we show that depletion of DnaQ in Mycolicibacterium smegmatis results in increased mutation rate, leading to AT-biased mutagenesis and elevated insertions/deletions in homopolymer tract. We demonstrated that mycobacterial DnaQ binds to the b-clamp and functions synergistically with the PHP domain to correct replication errors. Further, we found that the mycobacterial DnaQ sustains replicative fidelity upon chromosome topological stress. Intriguingly, we showed that a naturally evolved DnaQ variant prevalent in clinical Mycobacterium tuberculosis isolates enables hypermutability and is associated with extensive drug resistance. These results collectively establish that the alternative DnaQ functions in proofreading, and thus reveal that mycobacteria deploy two proofreaders to maintain replicative fidelity.

Genomic markers of drug resistance in Mycobacterium tuberculosis populations with minority variants

Minority variants of Mycobacterium tuberculosis harboring mutations conferring resistance can become dominant populations during tuberculosis (TB) treatment, leading to treatment failure. Our understanding of drug-resistant within-host subpopulations and the frequency of resistance-conferring mutations in minority variants remains limited. M. tuberculosis sequences recovered from liquid cultures of culture-confirmed TB cases notified between January 2017 and December 2021 in New South Wales, Australia were examined. Potential drug resistance-conferring minority variants were identified using LoFreq, and mixed populations of different M. tuberculosis strains (≥100 SNPs apart) were examined using QuantTB. A total of 1831 routinely sequenced M. tuberculosis strains were included in the analysis. Drug resistance-conferring minority variants were detected in 3.5% (65/1831) of sequenced cultures; 84.6% (55/65) had majority strains that were drug susceptible and 15.4% (10/65) had majority strains that were drug resistant. Minority variants with high-confidence drug resistance-conferring mutations were 1.5 times more common when the majority strains were drug resistant. Mixed M. tuberculosis strain populations were documented in 10.0% (183/1831) of specimens. Minority variants with high-confidence drug resistance-conferring mutations were more frequently detected in mixed M. tuberculosis strain populations (2.7%, 5/183) than in single strain populations (0.6%, 10/1648; P = 0.01). Drug-resistant minority variants require monitoring in settings that implement routine M. tuberculosis sequencing. The frequency with which drug-resistant minority variants are detected is likely influenced by pre-culture requirement. Culture-independent sequencing methods should provide a more accurate reflection of drug-resistant subpopulations.

The Mycobacterium tuberculosis genome at 25 years: lessons and lingering questions

First achieved in 1998 by Cole et al., the complete genome sequence of Mycobacterium tuberculosis continues to provide an invaluable resource to understand tuberculosis (TB), the leading cause of global infectious disease mortality. At the 25-year anniversary of this accomplishment, we describe how insights gleaned from the M. tuberculosis genome have led to vital tools for TB research, epidemiology, and clinical practice. The increasing accessibility of whole-genome sequencing across research and clinical settings has improved our ability to predict antibacterial susceptibility, to track epidemics at the level of individual outbreaks and wider historical trends, to query the efficacy of the bacille Calmette-Guérin (BCG) vaccine, and to uncover targets for novel antitubercular therapeutics. Likewise, we discuss several recent efforts to extract further discoveries from this powerful resource.

The Relative Positioning of Genotyping and Phenotyping for Tuberculosis Resistance Screening in Two EU National Reference Laboratories in 2023

The routine use of whole genome sequencing (WGS) as a reference typing technique for Mycobacterium tuberculosis epidemiology combined with the catalogued and extensive knowledge base of resistance-associated mutations means an initial susceptibility prediction can be derived from all cultured isolates in our laboratories based on WGS data alone. Preliminary work has confirmed, in our low-burden settings, these predictions are for first-line drugs, reproducible, robust with an accuracy similar to phenotypic drug susceptibility testing (pDST) and in many cases able to also predict the level of resistance (MIC). Routine screening for drug resistance by WGS results in approximately 80% of the isolates received being predicted as fully susceptible to the first-line drugs. Parallel testing with both WGS and pDST has demonstrated that routine pDST of genotypically fully susceptible isolates yields minimal additional information. Thus, rather than re-confirming all fully sensitive WGS-based predictions, we suggest that a more efficient use of available mycobacterial culture capacity in our setting is the development of a more extensive and detailed pDST targeted at any mono or multi-drug-resistant isolates identified by WGS screening. Phenotypic susceptibility retains a key role in the determination of an extended susceptibility profile for mono/multi-drugresistant isolates identified by WGS screening. The pDST information collected is also needed to support the development of future catalogues of resistance-associated mutations.

Rapid Detection of Extensive Drug Resistance by Xpert MTB/XDR Optimizes Therapeutic Decision-Making in Rifampin-Resistant Tuberculosis Patients

The Xpert MTB/XDR assay met the critical need for etiologic diagnosis of tuberculosis and rifampin resistance in previous studies. However, its benefits in tailoring the treatment regimen and improving the outcome for patients with rifampin-resistant tuberculosis (RR-TB) require further investigation. In this study, the Xpert MTB/XDR assay was used to determine the resistance profile of second-line drugs for RR-TB patients in two registered multicenter clinical trials, TB-TRUST (NCT03867136) and TB-TRUST-plus (NCT04717908), with the aim of testing the efficacy of all-oral shorter regimens in RR-TB patients in China. Patients would receive the fluoroquinolone-based all-oral shorter regimen, the injectable-containing regimen, or the bedaquiline-based regimen depending on fluoroquinolone susceptibility by using Xpert MTB/XDR. Among the 497 patients performed with Xpert MTB/XDR, 128 (25.8%) had infections resistant to fluoroquinolones and/or second-line injectable drugs (SLIDs). A total of 371 participants were recruited for the trials, and whole-genome sequencing (WGS) was performed on all corresponding culture-positive baseline strains. Taking the WGS results as the standard, the accuracy of the Xpert MTB/XDR assay in terms of resistance detection was 95.2% to 99.0% for all drugs. A total of 33 cases had inconsistent results, 9 of which were due to resistance heterogeneity. Most of the patients (241/281, 85.8%) had sputum culture conversion at 2 months. In conclusion, the Xpert MTB/XDR assay has the potential to serve as a quick reflex test in patients with RR-TB, as detected via Xpert MTB/RIF, to provide a reliable drug susceptibility profile of the infecting Mycobacterium tuberculosis strain and to initiate optimized treatment promptly.

Quantification of Isoniazid-Heteroresistant Mycobacterium tuberculosis Using Droplet Digital PCR

The quantitative detection of drug-resistance mutations in Mycobacterium tuberculosis (MTB) is critical for determining the drug resistance status of a sample. We developed a drop-off droplet digital PCR (ddPCR) assay targeting all major isoniazid (INH)-resistant mutations. The ddPCR assay consisted of three reactions: reaction A detects mutations at katG S315; reaction B detects inhA promoter mutations; and reaction C detects ahpC promoter mutations. All reactions could quantify 1%-50% of mutants in the presence of the wild-type, ranging from 100 to 50,000 copies/reaction. Clinical evaluation with 338 clinical isolates yielded clinical sensitivity of 94.5% (95% confidence interval [CI] = 89.1%-97.3%) and clinical specificity of 97.6% (95% CI = 94.6%-99.0%) compared with the traditional drug susceptibility testing (DST). Further clinical evaluation using 194 nucleic acid-positive MTB sputum samples revealed clinical sensitivity of 87.8% (95% CI = 75.8%-94.3%) and clinical specificity of 96.5% (95% CI = 92.2%-98.5%) in comparison with DST. All the mutant and heteroresistant samples detected by the ddPCR assay but susceptible by DST were confirmed by combined molecular assays, including Sanger sequencing, mutant-enriched Sanger sequencing and a commercial melting curve analysis-based assay. Finally, the ddPCR assay was used to monitor longitudinally the INH-resistance status and the bacterial load in nine patients undergoing treatment. Overall, the developed ddPCR assay could be an indispensable tool for quantification of INH-resistant mutations in MTB and bacterial loads in patients.

Mutation bias and the predictability of evolution

Predicting evolutionary outcomes is an important research goal in a diversity of contexts. The focus of evolutionary forecasting is usually on adaptive processes, and efforts to improve prediction typically focus on selection. However, adaptive processes often rely on new mutations, which can be strongly influenced by predictable biases in mutation. Here, we provide an overview of existing theory and evidence for such mutation-biased adaptation and consider the implications of these results for the problem of prediction, in regard to topics such as the evolution of infectious diseases, resistance to biochemical agents, as well as cancer and other kinds of somatic evolution. We argue that empirical knowledge of mutational biases is likely to improve in the near future, and that this knowledge is readily applicable to the challenges of short-term prediction. This article is part of the theme issue 'Interdisciplinary approaches to predicting evolutionary biology'.

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.

Tuberculosis treatment failure associated with evolution of antibiotic resilience

The widespread use of antibiotics has placed bacterial pathogens under intense pressure to evolve new survival mechanisms. Genomic analysis of 51,229 Mycobacterium tuberculosis (Mtb)clinical isolates has identified an essential transcriptional regulator, Rv1830, herein called resR for resilience regulator, as a frequent target of positive (adaptive) selection. resR mutants do not show canonical drug resistance or drug tolerance but instead shorten the post-antibiotic effect, meaning that they enable Mtb to resume growth after drug exposure substantially faster than wild-type strains. We refer to this phenotype as antibiotic resilience. ResR acts in a regulatory cascade with other transcription factors controlling cell growth and division, which are also under positive selection in clinical isolates of Mtb. Mutations of these genes are associated with treatment failure and the acquisition of canonical drug resistance.

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.

Rapid adaptation of a complex trait during experimental evolution of Mycobacterium tuberculosis

Tuberculosis (TB), caused by Mycobacterium tuberculosis (M. tb), is a leading cause of death due to infectious disease. TB is not traditionally associated with biofilms, but M. tb biofilms are linked with drug and immune tolerance and there is increasing recognition of their contribution to the recalcitrance of TB infections. Here, we used M. tb experimental evolution to investigate this complex phenotype and identify candidate loci controlling biofilm formation. We identified novel candidate loci, adding to our understanding of the genetic architecture underlying M. tb biofilm development. Under selective pressure to grow as a biofilm, regulatory mutations rapidly swept to fixation and were associated with changes in multiple traits, including extracellular matrix production, cell size, and growth rate. Genetic and phenotypic paths to enhanced biofilm growth varied according to the genetic background of the parent strain, suggesting that epistatic interactions are important in M. tb adaptation to changing environments.

Whole Genome Sequencing Identifies Novel Mutations Associated With Bedaquiline Resistance in Mycobacterium tuberculosis

Bedaquiline (BDQ), a new antitubercular agent, has been used to treat drug-resistant tuberculosis (TB). Although mutations in atpE, rv0678, and pepQ confer major resistance to BDQ, the mechanisms of resistance to BDQ in vitro and in clinical settings have not been fully elucidated. We selected BDQ-resistant mutants from 7H10 agar plates containing 0.5 mg/L BDQ (the critical concentration) and identified mutations associated with BDQ resistance through whole genome sequencing and Sanger sequencing. A total of 1,025 mutants were resistant to BDQ. We randomly selected 168 mutants for further analysis and discovered that 157/168 BDQ-resistant mutants harbored mutations in rv0678, which encodes a transcriptional regulator that represses the expression of the efflux pump, MmpS5–MmpL5. Moreover, we found two mutations with high frequency in rv0678 at nucleotide positions 286–287 (CG286–287 insertion; accounting for 26.8% [45/168]) and 198–199 (G198, G199 insertion, and G198 deletion; accounting for 14.3% [24/168]). The other mutations were dispersed covering the entire rv0678 gene. Moreover, we found that one new gene, glpK, harbors a G572 insertion; this mutation has a high prevalence (85.7%; 144/168) in the isolated mutants, and the minimum inhibitory concentration (MIC) assay demonstrated that it is closely associated with BDQ resistance. In summary, we characterized 168/1,025 mutants resistant to BDQ and found that mutations in rv0678 confer the primary mechanism of BDQ resistance. Moreover, we identified a new gene (glpK) involved in BDQ resistance. Our study offers new insights and valuable information that will contribute to rapid identification of BDQ-resistant isolates in clinical settings.

The microbial ecology of Escherichia coli in the vertebrate gut

Escherichia coli has a rich history as biology's 'rock star', driving advances across many fields. In the wild, E. coli resides innocuously in the gut of humans and animals but is also a versatile pathogen commonly associated with intestinal and extraintestinal infections and antimicrobial resistance-including large foodborne outbreaks such as the one that swept across Europe in 2011, killing 54 individuals and causing approximately 4000 infections and 900 cases of haemolytic uraemic syndrome. Given that most E. coli are harmless gut colonizers, an important ecological question plaguing microbiologists is what makes E. coli an occasionally devastating pathogen? To address this question requires an enhanced understanding of the ecology of the organism as a commensal. Here, we review how our knowledge of the ecology and within-host diversity of this organism in the vertebrate gut has progressed in the 137 years since E. coli was first described. We also review current approaches to the study of within-host bacterial diversity. In closing, we discuss some of the outstanding questions yet to be addressed and prospects for future research.

The impact of frequently neglected model violations on bacterial recombination rate estimation: a case study in Mycobacterium canettii and Mycobacterium tuberculosis

Mycobacterium canettii is a causative agent of tuberculosis in humans, along with the members of the Mycobacterium tuberculosis complex. Frequently used as an outgroup to the M. tuberculosis complex in phylogenetic analyses, M. canettii is thought to offer the best proxy for the progenitor species that gave rise to the complex. Here, we leverage whole-genome sequencing data and biologically relevant population genomic models to compare the evolutionary dynamics driving variation in the recombining M. canettii with that in the nonrecombining M. tuberculosis complex, and discuss differences in observed genomic diversity in the light of expected levels of Hill-Robertson interference. In doing so, we highlight the methodological challenges of estimating recombination rates through traditional population genetic approaches using sequences called from populations of microorganisms and evaluate the likely mis-inference that arises owing to a neglect of common model violations including purifying selection, background selection, progeny skew, and population size change. In addition, we compare performance when full within-host polymorphism data are utilized, versus the more common approach of basing analyses on within-host consensus sequences.

Exploiting genomics to mitigate the public health impact of antimicrobial resistance

Antimicrobial resistance (AMR) is a major global public health threat, which has been largely driven by the excessive use of antimicrobials. Control measures are urgently needed to slow the trajectory of AMR but are hampered by an incomplete understanding of the interplay between pathogens, AMR encoding genes, and mobile genetic elements at a microbial level. These factors, combined with the human, animal, and environmental interactions that underlie AMR dissemination at a population level, make for a highly complex landscape. Whole-genome sequencing (WGS) and, more recently, metagenomic analyses have greatly enhanced our understanding of these processes, and these approaches are informing mitigation strategies for how we better understand and control AMR. This review explores how WGS techniques have advanced global, national, and local AMR surveillance, and how this improved understanding is being applied to inform solutions, such as novel diagnostic methods that allow antimicrobial use to be optimised and vaccination strategies for better controlling AMR. We highlight some future opportunities for AMR control informed by genomic sequencing, along with the remaining challenges that must be overcome to fully realise the potential of WGS approaches for international AMR control.

Eosinophils are part of the granulocyte response in tuberculosis and promote host resistance in mice

Host resistance to Mycobacterium tuberculosis (Mtb) infection requires the activities of multiple leukocyte subsets, yet the roles of the different innate effector cells during tuberculosis are incompletely understood. Here we uncover an unexpected association between eosinophils and Mtb infection. In humans, eosinophils are decreased in the blood but enriched in resected human tuberculosis lung lesions and autopsy granulomas. An influx of eosinophils is also evident in infected zebrafish, mice, and nonhuman primate granulomas, where they are functionally activated and degranulate. Importantly, using complementary genetic models of eosinophil deficiency, we demonstrate that in mice, eosinophils are required for optimal pulmonary bacterial control and host survival after Mtb infection. Collectively, our findings uncover an unexpected recruitment of eosinophils to the infected lung tissue and a protective role for these cells in the control of Mtb infection in mice.

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.

Rifampicin exposure reveals within-host Mycobacterium tuberculosis diversity in patients with delayed culture conversion

Mycobacterium tuberculosis (Mtb) genetic micro-diversity in clinical isolates may underline mycobacterial adaptation to tuberculosis (TB) infection and provide insights to anti-TB treatment response and emergence of resistance. Herein we followed within-host evolution of Mtb clinical isolates in two cohorts of TB patients, either with delayed Mtb culture conversion (> 2 months), or with fast culture conversion (< 2 months). We captured the genetic diversity of Mtb isolates obtained in each patient, by focusing on minor variants detected as unfixed single nucleotide polymorphisms (SNPs). To unmask antibiotic tolerant sub-populations, we exposed these isolates to rifampicin (RIF) prior to whole genome sequencing (WGS) analysis. Thanks to WGS, we detected at least 1 unfixed SNP within the Mtb isolates for 9/15 patients with delayed culture conversion, and non-synonymous (ns) SNPs for 8/15 patients. Furthermore, RIF exposure revealed 9 additional unfixed nsSNP from 6/15 isolates unlinked to drug resistance. By contrast, in the fast culture conversion cohort, RIF exposure only revealed 2 unfixed nsSNP from 2/20 patients. To better understand the dynamics of Mtb micro-diversity, we investigated the variant composition of a persistent Mtb clinical isolate before and after controlled stress experiments mimicking the course of TB disease. A minor variant, featuring a particular mycocerosates profile, became enriched during both RIF exposure and macrophage infection. The variant was associated with drug tolerance and intracellular persistence, consistent with the pharmacological modeling predicting increased risk of treatment failure. A thorough study of such variants not necessarily linked to canonical drug-resistance, but which are prone to promote anti-TB drug tolerance, may be crucial to prevent the subsequent emergence of resistance. Taken together, the present findings support the further exploration of Mtb micro-diversity as a promising tool to detect patients at risk of poorly responding to anti-TB treatment, ultimately allowing improved and personalized TB management.

Tuberculosis (TB) is caused by Mycobacterium tuberculosis (Mtb), bacteria that are able to persist inside the patient for many months or years, thus requiring long antibiotic treatments. Here we focused on TB patients with delayed response to treatment and we performed genetic characterization of Mtb isolates to search for sub-populations that may tolerate anti-TB drugs. We found that Mtb cultured from 9/15 patients contained different sub-populations, and in vitro drug exposure revealed Mtb sub-populations in 6/15 isolates, none related to known drug-resistance mechanisms. By contrast, drug exposure revealed Mtb sup-populations in 2/20 isolates in the control cohort of patients with fast culture conversion. Furthermore, we characterized a Mtb variant isolated from a sub-population growing in the presence of rifampicin (RIF), a major anti-TB drug. We found that this variant featured a modified lipidic envelope, and that it was able to develop in the presence of RIF and inside human macrophage cells. We performed pharmacological modelling and found that this kind of variant may be related to a poor response to treatment. In conclusion, searching for particular Mtb sub-populations may help to detect patients at risk of treatment failure and provide additional guidance for TB management.

Genomic analyses of Mycobacterium tuberculosis from human lung resections reveal a high frequency of polyclonal infections

Polyclonal infections occur when at least two unrelated strains of the same pathogen are detected in an individual. This has been linked to worse clinical outcomes in tuberculosis, as undetected strains with different antibiotic resistance profiles can lead to treatment failure. Here, we examine the amount of polyclonal infections in sputum and surgical resections from patients with tuberculosis in the country of Georgia. For this purpose, we sequence and analyse the genomes of Mycobacterium tuberculosis isolated from the samples, acquired through an observational clinical study (NCT02715271). Access to the lung enhanced the detection of multiple strains (40% of surgery cases) as opposed to just using a sputum sample (0-5% in the general population). We show that polyclonal infections often involve genetically distant strains and can be associated with reversion of the patient's drug susceptibility profile over time. In addition, we find different patterns of genetic diversity within lesions and across patients, including mutational signatures known to be associated with oxidative damage; this suggests that reactive oxygen species may be acting as a selective pressure in the granuloma environment. Our results support the idea that the magnitude of polyclonal infections in high-burden tuberculosis settings is underestimated when only testing sputum samples.

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.

Emergence of additional drug resistance during treatment of multidrug-resistant tuberculosis in China: a prospective cohort study

OBJECTIVES

Little is known about how additional second-line drug resistance emerges during multidrug-resistant tuberculosis (MDR-TB) treatment. The present study aimed to investigate the influence of microevolution, exogenous reinfection and mixed infection on second-line drug resistance during the recommended 2-year MDR-TB treatment.

METHODS

Individuals with MDR-TB were enrolled between 2013 and 2016 in a multicentre prospective observational cohort study and were followed up for 2 years until treatment completion. Whole-genome sequencing (WGS) was applied for serial Mycobacterium tuberculosis isolates from study participants throughout the treatment, to study the role of microevolution, exogenous reinfection and mixed infection in the development of second-line drug resistance.

RESULTS

Of the 286 enrolled patients with MDR-TB, 63 (22.0%) M. tuberculosis isolates developed additional drug resistance during the MDR-TB treatment, including 5 that fulfilled the criteria of extensively drug-resistant TB. By comparing WGS data of serial isolates retrieved from the patients throughout treatment, 41 (65.1%) of the cases of additional second-line drug resistance were the result of exogenous reinfection, 18 (28.6%) were caused by acquired drug resistance, i.e. microevolution, while the remaining 4 (6.3%) were caused by mixed infections with drug-resistant and drug-susceptible strains. In multivariate analysis, previous TB treatment (adjusted hazard ratio (aHR) 2.51, 95% CI 1.51-4.18), extensive disease on chest X-ray (aHR 3.39, 95% CI 2.03-5.66) and type 2 diabetes mellitus (aHR 4.00, 95% CI 2.22-7.21) were independent risk factors associated with the development of additional second-line drug resistance.

CONCLUSIONS

A large proportion of additional second-line drug resistance emerging during MDR-TB treatment was attributed to exogenous reinfection, indicating the urgency of infection control in health facilities as well as the need for repeated drug susceptibility testing throughout MDR-TB treatment.

Mycobacterial fatty acid catabolism is repressed by FdmR to sustain lipogenesis and virulence

Host-derived fatty acids are an important carbon source for pathogenic mycobacteria during infection. How mycobacterial cells regulate the catabolism of fatty acids to serve the pathogenicity, however, remains unknown. Here, we identified a TetR-family transcriptional factor, FdmR, as the key regulator of fatty acid catabolism in the pathogen Mycobacterium marinum by combining use of transcriptomics, chromatin immunoprecipitation followed by sequencing, dynamic 13C-based flux analysis, metabolomics, and lipidomics. An M. marinum mutant deficient in FdmR was severely attenuated in zebrafish larvae and adult zebrafish. The mutant showed defective growth but high substrate consumption on fatty acids. FdmR was identified as a long-chain acyl-coenzyme A (acyl-CoA)-responsive repressor of genes involved in fatty acid degradation and modification. We demonstrated that FdmR functions as a valve to direct the flux of exogenously derived fatty acids away from β-oxidation toward lipid biosynthesis, thereby avoiding the overactive catabolism and accumulation of biologically toxic intermediates. Moreover, we found that FdmR suppresses degradation of long-chain acyl-CoAs endogenously synthesized through the type I fatty acid synthase. By modulating the supply of long-chain acyl-CoAs for lipogenesis, FdmR controls the abundance and chain length of virulence-associated lipids and mycolates and plays an important role in the impermeability of the cell envelope. These results reveal that despite the fact that host-derived fatty acids are used as an important carbon source, overactive catabolism of fatty acids is detrimental to mycobacterial cell growth and pathogenicity. This study thus presents FdmR as a potentially attractive target for chemotherapy.

Local adaptation of Mycobacterium tuberculosis on the Tibetan Plateau

During its global dispersal, Mycobacterium tuberculosis (Mtb) has encountered varied geographic environments and host populations. Although local adaptation seems to be a plausible model for describing long-term host-pathogen interactions, genetic evidence for this model is lacking. Here, we analyzed 576 whole-genome sequences of Mtb strains sampled from different regions of high-altitude Tibet. Our results show that, after sequential introduction of a few ancestral strains, the Tibetan Mtb population diversified locally while maintaining strict separation from the Mtb populations on the lower altitude plain regions of China. The current population structure and estimated past population dynamics suggest that the modern Beijing sublineage strains, which expanded over most of China and other global regions, did not show an expansion advantage in Tibet. The mutations in the Tibetan strains showed a higher proportion of A > G/T > C transitions than strains from the plain regions, and genes encoding DNA repair enzymes showed evidence of positive selection. Moreover, the long-term Tibetan exclusive selection for truncating mutations in the thiol-oxidoreductase encoding sseA gene suggests that Mtb was subjected to local selective pressures associated with oxidative stress. Collectively, the population genomics of Mtb strains in the relatively isolated population of Tibet provides genetic evidence that Mtb has adapted to local environments.

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.

In-host population dynamics of Mycobacterium tuberculosis complex during active disease

Tuberculosis (TB) is a leading cause of death globally. Understanding the population dynamics of TB’s causative agent Mycobacterium tuberculosis complex (Mtbc) in-host is vital for understanding the efficacy of antibiotic treatment. We use longitudinally collected clinical Mtbc isolates that underwent Whole-Genome Sequencing from the sputa of 200 patients to investigate Mtbc diversity during the course of active TB disease after excluding 107 cases suspected of reinfection, mixed infection or contamination. Of the 178/200 patients with persistent clonal infection >2 months, 27 developed new resistance mutations between sampling with 20/27 occurring in patients with pre-existing resistance. Low abundance resistance variants at a purity of ≥19% in the first isolate predict fixation in the subsequent sample. We identify significant in-host variation in 27 genes, including antibiotic resistance genes, metabolic genes and genes known to modulate host innate immunity and confirm several to be under positive selection by assessing phylogenetic convergence across a genetically diverse sample of 20,352 isolates.

Analysis of Serial Multidrug-Resistant Tuberculosis Strains Causing Treatment Failure and Within-Host Evolution by Whole-Genome Sequencing

The cure rate of multidrug-resistant tuberculosis (MDR-TB) is relatively low in China. The reasons for the treatment failure and within-host evolution during treatment have not been sufficiently studied. All MDR-TB patients receiving standard treatment from January 2014 to September 2016 at a designated TB Hospital in Zhejiang Province were retrospectively included and grouped according to their known treatment outcome. Clinical information was collected. Baseline strains of all patients and serial strains of treatment-failure patients were revived. Drug susceptibility tests (DSTs) of 14 drugs and single nucleotide polymorphism (SNP) analysis based on whole-genome sequencing (WGS) were performed. The genetic distance and within-host evolution were investigated based on SNPs. In total, 20 treatment failure patients and 74 patients who succeeded in treatment were included. The number of effective drugs for patients who failed treatment was no more than three. Eighteen (90.0%) treatment-failure patients were characterized by a continuous infection of the primary strain, of which 14 patients (77.8%) developed phenotypic or genotypic acquired drug resistance under ineffective treatment. Acquired resistance to amikacin and moxifloxacin (2.0 mg/ml) was detected most frequently, in 5 and 4 patients, respectively. The insufficient number of effective drugs in the combined treatment regimen was the main reason for MDR-TB treatment failure. The study emphasizes the importance of DST for second-line drugs when implementing the second-line drug regimen in MDR-TB patients. For patients with risk factors for MDR-TB, DST of second-line antituberculosis drugs should be performed at initiation of treatment. Second-line drugs should be selected based on the results of DST to avoid acquired resistance. WGS detects low-frequency resistance mutations and heterogeneous resistance with high sensitivity, which is of great significance for guiding clinical treatment and preventing acquired resistance.IMPORTANCE Few studies have focused on the reasons for the low cure rate of multidrug-resistant tuberculosis in China and within-host evolution during treatment, which is of great significance for improving clinical treatment regimens. Acquired resistance events were common during the ineffective treatment, among which resistance to amikacin and high-level moxifloxacin were the most common. The main reason for the treatment failure of MDR-TB patients was insufficient effective drugs, which may lead to higher levels of drug resistance in MDR-TB strains. Therefore, the study emphasizes the importance of DST in the development of second-line treatment regimen when there is a risk of MDR. By performing whole-genome sequencing of serial strains from patients with treatment failure, we found that WGS can detect low-frequency resistance mutations and heterogeneous resistance with high sensitivity. It is thus recommended to conduct drug susceptibility tests at the beginning of treatment and repeat the DST when the sputum bacteria remain positive.

Sources of Multidrug Resistance in Patients With Previous Isoniazid-Resistant Tuberculosis Identified Using Whole Genome Sequencing: A Longitudinal Cohort Study

BACKGROUND

Meta-analysis of patients with isoniazid-resistant tuberculosis (TB) given standard first-line anti-TB treatment indicated an increased risk of multidrug-resistant TB (MDR-TB) emerging (8%), compared to drug-sensitive TB (0.3%). Here we use whole genome sequencing (WGS) to investigate whether treatment of patients with preexisting isoniazid-resistant disease with first-line anti-TB therapy risks selecting for rifampicin resistance, and hence MDR-TB.

METHODS

Patients with isoniazid-resistant pulmonary TB were recruited and followed up for 24 months. Drug susceptibility testing was performed by microscopic observation drug susceptibility assay, mycobacterial growth indicator tube, and by WGS on isolates at first presentation and in the case of re-presentation. Where MDR-TB was diagnosed, WGS was used to determine the genomic relatedness between initial and subsequent isolates. De novo emergence of MDR-TB was assumed where the genomic distance was 5 or fewer single-nucleotide polymorphisms (SNPs), whereas reinfection with a different MDR-TB strain was assumed where the distance was 10 or more SNPs.

RESULTS

Two hundred thirty-nine patients with isoniazid-resistant pulmonary TB were recruited. Fourteen (14/239 [5.9%]) patients were diagnosed with a second episode of TB that was multidrug resistant. Six (6/239 [2.5%]) were identified as having evolved MDR-TB de novo and 6 as having been reinfected with a different strain. In 2 cases, the genomic distance was between 5 and 10 SNPs and therefore indeterminate.

CONCLUSIONS

In isoniazid-resistant TB, de novo emergence and reinfection of MDR-TB strains equally contributed to MDR development. Early diagnosis and optimal treatment of isoniazid-resistant TB are urgently needed to avert the de novo emergence of MDR-TB during treatment.

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.

The population genomics of within-host Mycobacterium tuberculosis

Recent progress in genomic sequencing from patient samples has allowed for the first detailed insight into the within-host genetic diversity of Mycobacterium tuberculosis (M.TB), revealing remarkably low levels of variation. While this has often been attributed to low mutation rates, other factors have been described, including resistance evolution (i.e., selective sweeps), widespread purifying and background selection, and, more recently, progeny skew. Here we review recent findings pertaining to the processes governing the evolutionary dynamics of M.TB, discuss their implications for improving our understanding of this important human pathogen, and make recommendations for future work. Significantly, this emerging evolutionary framework involving the joint estimation of demographic, selective, and reproductive processes is forming a new paradigm for the study of within-host pathogen evolution that will be widely applicable across organisms.

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.

Whole-genome sequencing and Mycobacterium tuberculosis: Challenges in sample preparation and sequencing data analysis

The numbers of patients with tuberculosis (TB) caused by resistant strains are still alarming. Therefore, it is necessary to determine resistance more quickly and precisely, than it is with the currently used phenotypic and genotypic methods. In recent years, technological advances have been made and the whole-genome sequencing (WGS) method has been introduced as a part of routine diagnostics in clinical laboratories. Comparing a wide range of mycobacterial genomic variations with a reference genome leads to a consistent evaluation of molecular-epidemiology and resistance of Mycobacterium tuberculosis (M. tuberculosis) to a wide range of anti-TB drugs. The quality of the obtained sequencing data is closely related to the type of sample and the method used for DNA extraction and sequencing library preparation. Moreover, the correct interpretation of results is also influenced by a bioinformatic data processing. A large number of bioinformatics pipelines are currently available, the sensitivity of which varies due to the different sizes of databases containing relevant mutations. This review focuses on the individual steps included in the sequencing workflow and factors that may affect the interpretation of final results.

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.

Evaluating the contributions of purifying selection and progeny-skew in dictating within-host Mycobacterium tuberculosis evolution

The within-host evolutionary dynamics of tuberculosis (TB) remain unclear, and underlying biological characteristics render standard population genetic approaches based upon the Wright-Fisher model largely inappropriate. In addition, the compact genome combined with an absence of recombination is expected to result in strong purifying selection effects. Thus, it is imperative to establish a biologically relevant evolutionary framework incorporating these factors in order to enable an accurate study of this important human pathogen. Further, such a model is critical for inferring fundamental evolutionary parameters related to patient treatment, including mutation rates and the severity of infection bottlenecks. We here implement such a model and infer the underlying evolutionary parameters governing within-patient evolutionary dynamics. Results demonstrate that the progeny skew associated with the clonal nature of TB severely reduces genetic diversity and that the neglect of this parameter in previous studies has led to significant mis-inference of mutation rates. As such, our results suggest an underlying de novo mutation rate that is considerably faster than previously inferred, and a progeny distribution differing significantly from Wright-Fisher assumptions. This inference represents a more appropriate evolutionary null model, against which the periodic effects of positive selection, associated with drug-resistance for example, may be better assessed.

Contaminant DNA in bacterial sequencing experiments is a major source of false genetic variability

Background

Contaminant DNA is a well-known confounding factor in molecular biology and in genomic repositories. Strikingly, analysis workflows for whole-genome sequencing (WGS) data commonly do not account for errors potentially introduced by contamination, which could lead to the wrong assessment of allele frequency both in basic and clinical research.

Results

We used a taxonomic filter to remove contaminant reads from more than 4000 bacterial samples from 20 different studies and performed a comprehensive evaluation of the extent and impact of contaminant DNA in WGS. We found that contamination is pervasive and can introduce large biases in variant analysis. We showed that these biases can result in hundreds of false positive and negative SNPs, even for samples with slight contamination. Studies investigating complex biological traits from sequencing data can be completely biased if contamination is neglected during the bioinformatic analysis, and we demonstrate that removing contaminant reads with a taxonomic classifier permits more accurate variant calling. We used both real and simulated data to evaluate and implement reliable, contamination-aware analysis pipelines.

Conclusion

As sequencing technologies consolidate as precision tools that are increasingly adopted in the research and clinical context, our results urge for the implementation of contamination-aware analysis pipelines. Taxonomic classifiers are a powerful tool to implement such pipelines.