Kim Mallard

Rapid determination of anti-tuberculosis drug resistance from whole-genome sequences

Mycobacterium tuberculosis drug resistance (DR) challenges effective tuberculosis disease control. Current molecular tests examine limited numbers of mutations, and although whole genome sequencing approaches could fully characterise DR, data complexity has restricted their clinical application. A library (1,325 mutations) predictive of DR for 15 anti-tuberculosis drugs was compiled and validated for 11 of them using genomic-phenotypic data from 792 strains. A rapid online ‘TB-Profiler’ tool was developed to report DR and strain-type profiles directly from raw sequences. Using our DR mutation library, in silico diagnostic accuracy was superior to some commercial diagnostics and alternative databases. The library will facilitate sequence-based drug-susceptibility testing.

Large-scale whole genome sequencing of M. tuberculosis provides insights into transmission in a high prevalence area

To improve understanding of the factors influencing tuberculosis transmission and the role of pathogen variation, we sequenced all available specimens from patients diagnosed over 15 years in a whole district in Malawi. Mycobacterium tuberculosis lineages were assigned and transmission networks constructed, allowing ≤10 single nucleotide polymorphisms (SNPs) difference. We defined disease as due to recent infection if the network-determined source was within 5 years, and assessed transmissibility from forward transmissions resulting in disease. High-quality sequences were available for 1687 disease episodes (72% of all culture-positive episodes): 66% of patients linked to at least one other patient. The between-patient mutation rate was 0.26 SNPs/year (95% CI 0.21–0.31). We showed striking differences by lineage in the proportion of disease due to recent transmission and in transmissibility (highest for lineage-2 and lowest for lineage-1) that were not confounded by immigration, HIV status or drug resistance. Transmissions resulting in disease decreased markedly over time. DOI: http://dx.doi.org/10.7554/eLife.05166.001

SpolPred: rapid and accurate prediction of Mycobacterium tuberculosis spoligotypes from short genomic sequences

Summary: Spoligotyping is a well-established genotyping technique based on the presence of unique DNA sequences in Mycobacterium tuberculosis (Mtb), the causal agent of tuberculosis disease (TB). Although advances in sequencing technologies are leading to whole-genome bacterial characterization, tens of thousands of isolates have been spoligotyped, giving a global view of Mtb strain diversity. To bridge the gap, we have developed SpolPred, a software to predict the spoligotype from raw sequence reads. Our approach is compared with experimentally and de novo assembly determined strain types in a set of 44 Mtb isolates. In silico and experimental results are identical for almost all isolates (39/44). However, SpolPred detected five experimentally false spoligotypes and was more accurate and faster than the assembling strategy. Application of SpolPred to an additional seven isolates with no laboratory data led to types that clustered with identical experimental types in a phylogenetic analysis using single-nucleotide polymorphisms. Our results demonstrate the usefulness of the tool and its role in revealing experimental limitations. Availability and implementation: SpolPred is written in C and is available from www.pathogenseq.org/spolpred. Contact: francesc.coll@lshtm.ac.uk Supplementary information: Supplementary data are available at Bioinformatics Online.

Recurrence due to Relapse or Reinfection With Mycobacterium tuberculosis: A Whole-Genome Sequencing Approach in a Large, Population-Based Cohort With a High HIV Infection Prevalence and Active Follow-up

Background. Recurrent tuberculosis is a major health burden and may be due to relapse with the original strain or reinfection with a new strain. Methods. In a population-based study in northern Malawi, patients with tuberculosis diagnosed from 1996 to 2010 were actively followed after the end of treatment. Whole-genome sequencing with approximately 100-fold coverage was performed on all available cultures. Results of IS6110 restriction fragment-length polymorphism analyses were available for cultures performed up to 2008. Results. Based on our data, a difference of ≤10 single-nucleotide polymorphisms (SNPs) was used to define relapse, and a difference of >100 SNPs was used to define reinfection. There was no evidence of mixed infections among those classified as reinfections. Of 1471 patients, 139 had laboratory-confirmed recurrences: 55 had relapse, and 20 had reinfection; for 64 type of recurrence was unclassified. Almost all relapses occurred in the first 2 years. Human immunodeficiency virus infection was associated with reinfection but not relapse. Relapses were associated with isoniazid resistance, treatment before 2007, and lineage-3 strains. We identified several gene variants associated with relapse. Lineage-2 (Beijing) was overrepresented and lineage-1 underrepresented among the reinfecting strains (P = .004). Conclusions. While some of the factors determining recurrence depend on the patient and their treatment, differences in the Mycobacterium tuberculosis genome appear to have a role in both relapse and reinfection.

Whole Genome Sequencing Based Characterization of Extensively Drug- Resistant Mycobacterium tuberculosis Isolates from Pakistan

Improved molecular diagnostic methods for detection drug resistance in Mycobacterium tuberculosis (MTB) strains are required. Resistance to first- and second- line anti-tuberculous drugs has been associated with single nucleotide polymorphisms (SNPs) in particular genes. However, these SNPs can vary between MTB lineages therefore local data is required to describe different strain populations. We used whole genome sequencing (WGS) to characterize 37 extensively drug-resistant (XDR) MTB isolates from Pakistan and investigated 40 genes associated with drug resistance. Rifampicin resistance was attributable to SNPs in the rpoB hot-spot region. Isoniazid resistance was most commonly associated with the katG codon 315 (92%) mutation followed by inhA S94A (8%) however, one strain did not have SNPs in katG, inhA or oxyR-ahpC. All strains were pyrazimamide resistant but only 43% had pncA SNPs. Ethambutol resistant strains predominantly had embB codon 306 (62%) mutations, but additional SNPs at embB codons 406, 378 and 328 were also present. Fluoroquinolone resistance was associated with gyrA 91–94 codons in 81% of strains; four strains had only gyrB mutations, while others did not have SNPs in either gyrA or gyrB. Streptomycin resistant strains had mutations in ribosomal RNA genes; rpsL codon 43 (42%); rrs 500 region (16%), and gidB (34%) while six strains did not have mutations in any of these genes. Amikacin/kanamycin/capreomycin resistance was associated with SNPs in rrs at nt1401 (78%) and nt1484 (3%), except in seven (19%) strains. We estimate that if only the common hot-spot region targets of current commercial assays were used, the concordance between phenotypic and genotypic testing for these XDR strains would vary between rifampicin (100%), isoniazid (92%), flouroquinolones (81%), aminoglycoside (78%) and ethambutol (62%); while pncA sequencing would provide genotypic resistance in less than half the isolates. This work highlights the importance of expanded targets for drug resistance detection in MTB isolates.

Low-cost rapid detection of rifampicin resistant tuberculosis using bacteriophage in Kampala, Uganda.

BackgroundResistance to anti-tuberculosis drugs is a serious public health problem. Multi-drug resistant tuberculosis (MDR-TB), defined as resistance to at least rifampicin and isoniazid, has been reported in all regions of the world. Current phenotypic methods of assessing drug susceptibility of M. tuberculosis are slow. Rapid molecular methods to detect resistance to rifampicin have been developed but they are not affordable in some high prevalence countries such as those in sub Saharan Africa. A simple multi-well plate assay using mycobacteriophage D29 has been developed to test M. tuberculosis isolates for resistance to rifampicin. The purpose of this study was to investigate the performance of this technology in Kampala, Uganda.MethodsIn a blinded study 149 M. tuberculosis isolates were tested for resistance to rifampicin by the phage assay and results compared to those from routine phenotypic testing in BACTEC 460. Three concentrations of drug were used 2, 4 and 10 μg/ml. Isolates found resistant by either assay were subjected to sequence analysis of a 81 bp fragment of the rpoB gene to identify mutations predictive of resistance. Four isolates with discrepant phage and BACTEC results were tested in a second phenotypic assay to determine minimal inhibitory concentrations.ResultsInitial analysis suggested a sensitivity and specificity of 100% and 96.5% respectively for the phage assay used at 4 and 10 μg/ml when compared to the BACTEC 460. However, further analysis revealed 4 false negative results from the BACTEC 460 and the phage assay proved the more sensitive and specific of the two tests. Of the 39 isolates found resistant by the phage assay 38 (97.4%) were found to have mutations predictive of resistance in the 81 bp region of the rpoB gene. When used at 2 μg/ml false resistant results were observed from the phage assay. The cost of reagents for testing each isolate was estimated to be 1.3US

Molecular Detection of Mixed Infections of Mycobacterium tuberculosis Strains in Sputum Samples from Patients in Karonga District, Malawi

when testing a batch of 20 isolates on a single 96 well plate. Results were obtained in 48 hours.ConclusionThe phage assay can be used for screening of isolates for resistance to rifampicin, with high sensitivity and specificity in Uganda. The test may be useful in poorly resourced laboratories as a rapid screen to differentiate between rifampicin susceptible and potential MDR-TB cases.

Mycobacterium tuberculosis whole genome sequencing and protein structure modelling provides insights into anti-tuberculosis drug resistance

ABSTRACT The occurrence of mixed infections of Mycobacterium tuberculosis is no longer disputed. However, their frequency, and the impact they may have on our understanding of tuberculosis (TB) pathogenesis and epidemiology, remains undetermined. Most previous studies of frequency applied genotyping techniques to cultured M. tuberculosis isolates and found mixed infections to be rare. PCR-based techniques may be more sensitive for detecting multiple M. tuberculosis strains and can be applied to sputum. To date, one study in South Africa has used a PCR approach and suggested that mixed infection could be common. We investigated mixed infections in northern Malawi using two lineage-specific PCR assays targeting the Latin American-Mediterranean (LAM) and Beijing lineages. Compared with spoligotyping, the specificity and sensitivity of both assays was 100%. From 160 culture-positive sputa, mixed LAM and non-LAM strains were detected in 4 sputa belonging to 2 (2.8%) patients. Both patients were HIV positive, with no history of TB. Cultured isolates from both patients showed only LAM by PCR and spoligotyping. In a set of 377 cultured isolates, 4 were mixed LAM and non-LAM. Only one showed evidence of more than one M. tuberculosis strain using IS6110-based restriction fragment length polymorphism (IS6110-RFLP) and spoligotyping analyses. Corresponding sputa for the 4 isolates were unavailable. Mixed Beijing and non-Beijing strains were not detected in this study. Mixed infections appear to be rare in our setting and are unlikely to affect findings based on DNA fingerprinting data. Molecular methods, which avoid the selective nature of culture and target distinct strains, are well suited to detection of mixed infections.

Tuberculosis in Sudan: a study of Mycobacterium tuberculosis strain genotype and susceptibility to anti-tuberculosis drugs.

BackgroundCombating the spread of drug resistant tuberculosis is a global health priority. Whole genome association studies are being applied to identify genetic determinants of resistance to anti-tuberculosis drugs. Protein structure and interaction modelling are used to understand the functional effects of putative mutations and provide insight into the molecular mechanisms leading to resistance.MethodsTo investigate the potential utility of these approaches, we analysed the genomes of 144 Mycobacterium tuberculosis clinical isolates from The Special Programme for Research and Training in Tropical Diseases (TDR) collection sourced from 20 countries in four continents. A genome-wide approach was applied to 127 isolates to identify polymorphisms associated with minimum inhibitory concentrations for first-line anti-tuberculosis drugs. In addition, the effect of identified candidate mutations on protein stability and interactions was assessed quantitatively with well-established computational methods.ResultsThe analysis revealed that mutations in the genes rpoB (rifampicin), katG (isoniazid), inhA-promoter (isoniazid), rpsL (streptomycin) and embB (ethambutol) were responsible for the majority of resistance observed. A subset of the mutations identified in rpoB and katG were predicted to affect protein stability. Further, a strong direct correlation was observed between the minimum inhibitory concentration values and the distance of the mutated residues in the three-dimensional structures of rpoB and katG to their respective drugs binding sites.ConclusionsUsing the TDR resource, we demonstrate the usefulness of whole genome association and convergent evolution approaches to detect known and potentially novel mutations associated with drug resistance. Further, protein structural modelling could provide a means of predicting the impact of polymorphisms on drug efficacy in the absence of phenotypic data. These approaches could ultimately lead to novel resistance mutations to improve the design of tuberculosis control measures, such as diagnostics, and inform patient management.

Elucidating emergence and transmission of multidrug-resistant tuberculosis in treatment experienced patients by whole genome sequencing

BackgroundSudan is a large country with a diverse population and history of civil conflict. Poverty levels are high with a gross national income per capita of less than two thousand dollars. The country has a high burden of tuberculosis (TB) with an estimated 50,000 incident cases during 2009, when the estimated prevalence was 209 cases per 100,000 of the population. Few studies have been undertaken on TB in Sudan and the prevalence of drug resistant disease is not known.MethodsIn this study Mycobacterium tuberculosis isolates from 235 patients attending three treatment centers in Sudan were screened for susceptibility to isoniazid, rifampicin, ethambutol and streptomycin by the proportion method on Lowenstein Jensen media. 232 isolates were also genotyped by spoligotyping. Demographic details of patients were recorded using a structured questionnaire. Statistical analyses were conducted to examine the associations between drug resistance with risk ratios computed for a set of risk factors (gender, age, case status - new or relapse, geographic origin of the patient, spoligotype, number of people per room, marital status and type of housing).ResultsMulti drug-resistant tuberculosis (MDR-TB), being resistance to at least rifampicin and isoniazid, was found in 5% (95% CI: 2,8) of new cases and 24% (95% CI: 14,34) of previously treated patients. Drug resistance was associated with previous treatment with risk ratios of 3.51 (95% CI: 2.69-4.60; p < 0.001) for resistance to any drug and 5.23 (95% CI: 2.30-11.90; p < 0.001) for MDR-TB. Resistance was also associated with the geographic region of origin of the patient, being most frequently observed in patients from the Northern region and least in the Eastern region with risk ratios of 7.43 (95%CI:3.42,16.18; p: < 0.001) and 14.09 (95%CI:1.80,110.53; p:0.026) for resistance to any drug and MDR-TB. The major genotype observed was of the Central Asia spoligotype family (CAS1_Delhi), representing 49% of the 232 isolates examined.ConclusionsWe conclude that emergence of drug resistant tuberculosis has the potential to be a serious public health problem in Sudan and that strengthened tuberculosis control and improved monitoring of therapy is needed. Further surveillance is required to fully ascertain the extent of the problem.