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.

Correction to: Whole genome sequencing Mycobacterium tuberculosis directly from sputum identifies more genetic diversity than sequencing from culture.

He authors reported that one of the authors' names was typeset incorrectly in the authorship list.

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.

Population-level emergence of bedaquiline and clofazimine resistance-associated variants among patients with drug-resistant tuberculosis in southern Africa: a phenotypic and phylogenetic analysis.

BACKGROUND: Bedaquiline and clofazimine are important drugs in the treatment of drug-resistant tuberculosis and are commonly used across southern Africa, although drug susceptibility testing is not routinely performed. In this study, we did a genotypic and phenotypic analysis of drug-resistant Mycobacterium tuberculosis isolates from cohort studies in hospitals in KwaZulu-Natal, South Africa, to identify resistance-associated variants (RAVs) and assess the extent of clofazimine and bedaquiline cross-resistance. We also used a comprehensive dataset of whole-genome sequences to investigate the phylogenetic and geographical distribution of bedaquiline and clofazimine RAVs in southern Africa. METHODS: In this study, we included M tuberculosis isolates reported from the PRAXIS study of patients with drug-resistant tuberculosis treated with bedaquiline (King Dinuzulu Hospital, Durban) and three other cohort studies of drug-resistant tuberculosis in other KwaZulu-Natal hospitals, and sequential isolates from six persistently culture-positive patients with extensively drug-resistant tuberculosis at the KwaZulu-Natal provincial referral laboratory. Samples were collected between 2013 and 2019. Microbiological cultures were done as part of all parent studies. We sequenced whole genomes of included isolates and measured bedaquiline and clofazimine minimum inhibitory concentrations (MICs) for isolates identified as carrying any Rv0678 variant or previously published atpE, pepQ, and Rv1979c RAVs, which were the subject of the phenotypic study. We combined all whole-genome sequences of M tuberculosis obtained in this study with publicly available sequence data from other tuberculosis studies in southern Africa (defined as the countries of the Southern African Development Community), including isolates with Rv0678 variants identified by screening public genomic databases. We used this extended dataset to reconstruct phylogenetic relationships across lineage 2 and 4 M tuberculosis isolates. FINDINGS: We sequenced the whole genome of 648 isolates from 385 patients with drug-resistant tuberculosis recruited into cohort studies in KwaZulu-Natal, and 28 isolates from six patients from the KwaZulu-Natal referral laboratory. We identified 30 isolates with Rv0678 RAVs from 16 (4%) of 391 patients. We did not identify any atpE, pepQ, or Rv1979c RAVs. MICs were measured for 21 isolates with Rv0678 RAVs. MICs were above the critical concentration for bedaquiline resistance in nine (43%) of 21 isolates, in the intermediate category in nine (43%) isolates, and within the wild-type range in three (14%) isolates. Clofazimine MICs in genetically wild-type isolates ranged from 0·12-0·5 µg/mL, and in isolates with RAVs from 0·25-4·0 µg/mL. Phylogenetic analysis of the extended dataset including M tuberculosis isolates from southern Africa resolved multiple emergences of Rv0678 variants in lineages 2 and 4, documented two likely nosocomial transmission events, and identified the spread of a possibly bedaquiline and clofazimine cross-resistant clone in eSwatini. We also identified four patients with pepQ frameshift mutations that may confer resistance. INTERPRETATION: Bedaquiline and clofazimine cross-resistance in southern Africa is emerging repeatedly, with evidence of onward transmission largely due to Rv0678 mutations in M tuberculosis. Roll-out of bedaquiline and clofazimine treatment in the setting of limited drug susceptibility testing could allow further spread of resistance. Designing strong regimens would help reduce the emergence of resistance. Drug susceptibility testing is required to identify where resistance does emerge. FUNDING: Wellcome Trust, National Institute of Allergy and Infectious Diseases and National Center for Advancing Translational Sciences of the National Institutes of Health.