Identifying mixed Mycobacterium tuberculosis infections from whole genome sequence data.

BACKGROUND: Mixed, polyclonal Mycobacterium tuberculosis infection occurs in natural populations. Developing an effective method for detecting such cases is important in measuring the success of treatment and reconstruction of transmission between patients. Using whole genome sequence (WGS) data, we assess two methods for detecting mixed infection: (i) a combination of the number of heterozygous sites and the proportion of heterozygous sites to total SNPs, and (ii) Bayesian model-based clustering of allele frequencies from sequencing reads at heterozygous sites. RESULTS: In silico and in vitro artificially mixed and known pure M. tuberculosis samples were analysed to determine the specificity and sensitivity of each method. We found that both approaches were effective in distinguishing between pure strains and mixed infection where there was relatively high (> 10%) proportion of a minor strain in the mixture. A large dataset of clinical isolates (n = 1963) from the Karonga Prevention Study in Northern Malawi was tested to examine correlations with patient characteristics and outcomes with mixed infection. The frequency of mixed infection in the population was found to be around 10%, with an association with year of diagnosis, but no association with age, sex, HIV status or previous tuberculosis. CONCLUSIONS: Mixed Mycobacterium tuberculosis infection was identified in silico using whole genome sequence data. The methods presented here can be applied to population-wide analyses of tuberculosis to estimate the frequency of mixed infection, and to identify individual cases of mixed infections. These cases are important when considering the evolution and transmission of the disease, and in patient treatment

The infectiousness of tuberculosis patients coinfected with HIV

The current understanding of airborne tuberculosis (TB) transmission is based on classic 1950s studies in which guinea pigs were exposed to air from a tuberculosis ward. Recently we recreated this model in Lima, Peru, and in this paper we report the use of molecular fingerprinting to investigate patient infectiousness in the current era of HIV infection and multidrug-resistant (MDR) TB

A Geographically-Restricted but Prevalent Mycobacterium tuberculosis Strain Identified in the West Midlands Region of the UK between 1995 and 2008

Background: We describe the identification of, and risk factors for, the single most prevalent Mycobacterium tuberculosis strain in the West Midlands region of the UK.Methodology/Principal Findings: Prospective 15-locus MIRU-VNTR genotyping of all M. tuberculosis isolates in the West Midlands between 2004 and 2008 was undertaken. Two retrospective epidemiological investigations were also undertaken using univariable and multivariable logistic regression analysis. The first study of all TB patients in the West Midlands between 2004 and 2008 identified a single prevalent strain in each of the study years (total 155/3,056 (5%) isolates). This prevalent MIRU-VNTR profile (32333 2432515314 434443183) remained clustered after typing with an additional 9-loci MIRU-VNTR and spoligotyping. The majority of these patients (122/155, 79%) resided in three major cities located within a 40 km radius. From the apparent geographical restriction, we have named this the "Mercian" strain. A multivariate analysis of all TB patients in the West Midlands identified that infection with a Mercian strain was significantly associated with being UK-born (OR = 9.03, 95% CI = 4.56-17.87, p 65 years old (OR = 0.25, 95% CI = 0.09-0.67, p < 0.01). A second more detailed investigation analyzed a cohort of 82 patients resident in Wolverhampton between 2003 and 2006. A significant association with being born in the UK remained after a multivariate analysis (OR = 9.68, 95% CI = 2.00-46.78, p < 0.01) and excess alcohol intake and cannabis use (OR = 6.26, 95% CI = 1.45-27.02, p = .01) were observed as social risk factors for infection.Conclusions/Significance: The continued consistent presence of the Mercian strain suggests ongoing community transmission. Whilst significant associations have been found, there may be other common risk factors yet to be identified. Future investigations should focus on targeting the relevant risk groups and elucidating the biological factors that mediate continued transmission of this strain

Detecting Mutations in the Mycobacterium tuberculosis Pyrazinamidase Gene pncA to Improve Infection Control and Decrease Drug Resistance Rates in Human Immunodeficiency Virus Coinfection.

Hospital infection control measures are crucial to tuberculosis (TB) control strategies within settings caring for human immunodeficiency virus (HIV)-positive patients, as these patients are at heightened risk of developing TB. Pyrazinamide (PZA) is a potent drug that effectively sterilizes persistent Mycobacterium tuberculosis bacilli. However, PZA resistance associated with mutations in the nicotinamidase/pyrazinamidase coding gene, pncA, is increasing. A total of 794 patient isolates obtained from four sites in Lima, Peru, underwent spoligotyping and drug resistance testing. In one of these sites, the HIV unit of Hospital Dos de Mayo (HDM), an isolation ward for HIV/TB coinfected patients opened during the study as an infection control intervention: circulating genotypes and drug resistance pre- and postintervention were compared. All other sites cared for HIV-negative outpatients: genotypes and drug resistance rates from these sites were compared with those from HDM. HDM patients showed high concordance between multidrug resistance, PZA resistance according to the Wayne method, the two most common genotypes (spoligotype international type [SIT] 42 of the Latino American-Mediterranean (LAM)-9 clade and SIT 53 of the T1 clade), and the two most common pncA mutations (G145A and A403C). These associations were absent among community isolates. The infection control intervention was associated with 58-92% reductions in TB caused by SIT 42 or SIT 53 genotypes (odds ratio [OR] = 0.420, P = 0.003); multidrug-resistant TB (OR = 0.349, P < 0.001); and PZA-resistant TB (OR = 0.076, P < 0.001). In conclusion, pncA mutation typing, with resistance testing and spoligotyping, was useful in identifying a nosocomial TB outbreak and demonstrating its resolution after implementation of infection control measures

Mapping of Mycobacterium tuberculosis Complex Genetic Diversity Profiles in Tanzania and Other African Countries

The aim of this study was to assess and characterize Mycobacterium tuberculosis complex (MTBC) genotypic diversity in Tanzania, as well as in neighbouring East and other several African countries. We used spoligotyping to identify a total of 293 M. tuberculosis clinical isolates (one isolate per patient) collected in the Bunda, Dar es Salaam, Ngorongoro and Serengeti areas in Tanzania. The results were compared with results in the SITVIT2 international database of the Pasteur Institute of Guadeloupe. Genotyping and phylogeographical analyses highlighted the predominance of the CAS, T, EAI, and LAM MTBC lineages in Tanzania. The three most frequent Spoligotype International Types (SITs) were: SIT21/CAS1-Kili (n = 76; 25.94%), SIT59/LAM11-ZWE (n = 22; 7.51%), and SIT126/EAI5 tentatively reclassified as EAI3-TZA (n = 18; 6.14%). Furthermore, three SITs were newly created in this study (SIT4056/EAI5 n = 2, SIT4057/T1 n = 1, and SIT4058/EAI5 n = 1). We noted that the East-African-Indian (EAI) lineage was more predominant in Bunda, the Manu lineage was more common among strains isolated in Ngorongoro, and the Central-Asian (CAS) lineage was more predominant in Dar es Salaam (p-value<0.0001). No statistically significant differences were noted when comparing HIV status of patients vs. major lineages (p-value = 0.103). However, when grouping lineages as Principal Genetic Groups (PGG), we noticed that PGG2/3 group (Haarlem, LAM, S, T, and X) was more associated with HIV-positive patients as compared to PGG1 group (Beijing, CAS, EAI, and Manu) (p-value = 0.03). This study provided mapping of MTBC genetic diversity in Tanzania (containing information on isolates from different cities) and neighbouring East African and other several African countries highlighting differences as regards to MTBC genotypic distribution between Tanzania and other African countries. This work also allowed underlining of spoligotyping patterns tentatively grouped within the newly designated EAI3-TZA lineage (remarkable by absence of spacers 2 and 3, and represented by SIT126) which seems to be specific to Tanzania. However, further genotyping information would be needed to confirm this specificity

ClassTR: Classifying Within-Host Heterogeneity Based on Tandem Repeats with Application to Mycobacterium tuberculosis Infections.

Genomic tools have revealed genetically diverse pathogens within some hosts. Within-host pathogen diversity, which we refer to as "complex infection", is increasingly recognized as a determinant of treatment outcome for infections like tuberculosis. Complex infection arises through two mechanisms: within-host mutation (which results in clonal heterogeneity) and reinfection (which results in mixed infections). Estimates of the frequency of within-host mutation and reinfection in populations are critical for understanding the natural history of disease. These estimates influence projections of disease trends and effects of interventions. The genotyping technique MLVA (multiple loci variable-number tandem repeats analysis) can identify complex infections, but the current method to distinguish clonal heterogeneity from mixed infections is based on a rather simple rule. Here we describe ClassTR, a method which leverages MLVA information from isolates collected in a population to distinguish mixed infections from clonal heterogeneity. We formulate the resolution of complex infections into their constituent strains as an optimization problem, and show its NP-completeness. We solve it efficiently by using mixed integer linear programming and graph decomposition. Once the complex infections are resolved into their constituent strains, ClassTR probabilistically classifies isolates as clonally heterogeneous or mixed by using a model of tandem repeat evolution. We first compare ClassTR with the standard rule-based classification on 100 simulated datasets. ClassTR outperforms the standard method, improving classification accuracy from 48% to 80%. We then apply ClassTR to a sample of 436 strains collected from tuberculosis patients in a South African community, of which 92 had complex infections. We find that ClassTR assigns an alternate classification to 18 of the 92 complex infections, suggesting important differences in practice. By explicitly modeling tandem repeat evolution, ClassTR helps to improve our understanding of the mechanisms driving within-host diversity of pathogens like Mycobacterium tuberculosis

Transmission of multidrug-resistant tuberculosis in the UK: a cross-sectional molecular and epidemiological study of clustering and contact tracing

BACKGROUND: Between 2000 and 2012 the number of multidrug-resistant (MDR) tuberculosis cases in the UK increased from 28 per year to 81 per year. We investigated the proportion of MDR tuberculosis cases arising from transmission in the UK and associated risk factors. METHOD: We identified patients with MDR tuberculosis notified in England, Wales, and Northern Ireland between Jan 1, 2004, and Dec 31, 2007, by linking national laboratory and surveillance data. Data for laboratory isolates, including drug sensitivities and 24-mycobacterial interspersed repetitive-unit-variable-number tandem repeat (MIRU-VNTR) typing were obtained routinely from the National Tuberculosis Reference laboratories as part of national tuberculosis surveillance. We investigated clusters of cases with indistinguishable MIRU-VNTR profiles to identify epidemiological links. We calculated transmission using the n-1 method and established associated risk factors by logistic regression. We also assessed the likelihood of transmission to additional secondary active tuberculosis cases, identified through conventional contact tracing. FINDINGS: 204 patients were diagnosed with MDR tuberculosis in the study period; 189 (92·6%) had an MIRU-VNTR profile. We identified 12 clusters containing 40 individuals and 149 unique strains. The proportion of cases attributable to recent transmission, on the basis of molecular data, was 15% (40 cases clustered-12 clusters/189 with a strain type). The proportion of cases attributable to recent transmission (ie, transmission within the UK) after adjustment for epidemiological links was 8·5% (22 cases with epidemiological links-six clusters/189 cases with a strain type). Being UK born (odds ratio 4·81; 95% CI 2·03-11·36, p=0·0005) and illicit drug use (4·75; 1·19-18·96, p=0·026) were significantly associated with clustering. The most common transmission setting was the household but 21 of 22 of epidemiological links were missed by conventional contact tracing. 13 secondary active tuberculosis cases identified by conventional contact tracing were mostly contacts of patients with MDR tuberculosis from countries of high tuberculosis burden. 11 (85%) of 13 shared the same country of birth as the index case, of whom ten did not share a strain type or drug resistance pattern. INTERPRETATION: Transmission of MDR tuberculosis in the UK is low and associated with being UK born or illicit drug use. MIRU-VNTR typing with cluster investigation was more successful at identifying transmission events than conventional contact tracing. Individuals with tuberculosis who have had contact with a known MDR tuberculosis source case from a country of high tuberculosis burden should have drug-sensitivity testing on isolates to ensure appropriate treatment is given. FUNDING: Public Health England

The relative fitness of drug-resistant Mycobacterium tuberculosis: a modelling study of household transmission in Peru.

The relative fitness of drug-resistant versus susceptible bacteria in an environment dictates resistance prevalence. Estimates for the relative fitness of resistant Mycobacterium tuberculosis (Mtb) strains are highly heterogeneous and mostly derived from in vitro experiments. Measuring fitness in the field allows us to determine how the environment influences the spread of resistance. We designed a household structured, stochastic mathematical model to estimate the fitness costs associated with multidrug resistance (MDR) carriage in Mtb in Lima, Peru during 2010-2013. By fitting the model to data from a large prospective cohort study of TB disease in household contacts, we estimated the fitness, relative to susceptible strains with a fitness of 1, of MDR-Mtb to be 0.32 (95% credible interval: 0.15-0.62) or 0.38 (0.24-0.61), if only transmission or progression to disease, respectively, was affected. The relative fitness of MDR-Mtb increased to 0.56 (0.42-0.72) when the fitness cost influenced both transmission and progression to disease equally. We found the average relative fitness of MDR-Mtb circulating within households in Lima, Peru during 2010-2013 to be significantly lower than concurrent susceptible Mtb If these fitness levels do not change, then existing TB control programmes are likely to keep MDR-TB prevalence at current levels in Lima, Peru