Search

Objectives AGILE is a Phase Ib/IIa platform for rapidly evaluating COVID-19 treatments. In this trial (NCT04746183) we evaluated the safety and optimal dose of molnupiravir in participants with early symptomatic infection. Methods We undertook a dose-escalating, open-label, randomized-controlled (standard-of-care) Bayesian adaptive Phase I trial at the Royal Liverpool and Broadgreen Clinical Research Facility. Participants (adult outpatients with PCR-confirmed SARS-CoV-2 infection within 5 days of symptom onset) were randomized 2:1 in groups of 6 participants to 300, 600 and 800 mg doses of molnupiravir orally, twice daily for 5 days or control. A dose was judged unsafe if the probability of 30% or greater dose-limiting toxicity (the primary outcome) over controls was 25% or greater. Secondary outcomes included safety, clinical progression, pharmacokinetics and virological responses. Results Of 103 participants screened, 18 participants were enrolled between 17 July and 30 October 2020. Molnupiravir was well tolerated at 300, 600 and 800 mg doses with no serious or severe adverse events. Overall, 4 of 4 (100%), 4 of 4 (100%) and 1 of 4 (25%) of the participants receiving 300, 600 and 800 mg molnupiravir, respectively, and 5 of 6 (83%) controls, had at least one adverse event, all of which were mild (≤grade 2). The probability of ≥30% excess toxicity over controls at 800 mg was estimated at 0.9%. Conclusions Molnupiravir was safe and well tolerated; a dose of 800 mg twice daily for 5 days was recommended for Phase II evaluation.

BackgroundAGILE is a phase Ib/IIa platform for rapidly evaluating COVID-19 treatments. In this trial (NCT04746183) we evaluated the safety and optimal dose of molnupiravir in participants with early symptomatic infection.MethodsWe undertook a dose-escalating, open-label, randomised-controlled (standard-of-care) Bayesian adaptive phase I trial at the Royal Liverpool and Broadgreen Clinical Research Facility. Participants (adult outpatients with PCR-confirmed SARS-CoV-2 infection within 5 days of symptom onset) were randomised 2:1 in groups of 6 participants to 300mg, 600mg and 800mg doses of molnupiravir orally, twice daily for 5 days or control. A dose was judged unsafe if the probability of 30% or greater dose-limiting toxicity (the primary outcome) over controls was higher than 25%. Secondary outcomes included safety, clinical progression, pharmacokinetics and virologic responses.ResultsOf 103 volunteers screened, 18 participants were enrolled between 17 July and 30 October 2020. Molnupiravir was well tolerated at 400, 600 or 800mg doses with no serious or severe adverse events. Overall, 4 of 4 (100%), 4 of 4 (100%) and 1 of 4 (25%) of the participants receiving 300, 600 and 800mg molnupiravir respectively, and 5 of 6 (83%) controls, had at least one adverse event, all of which were mild (≤grade 2). The probability of ≥30% excess toxicity over controls at 800mg was estimated at 0.9%.ConclusionMolnupiravir was safe and well tolerated; a dose of 800mg twice-daily for 5 days was recommended for Phase II evaluation.

Objectives To estimate UK prevalence and incidence of clinically significant carbapenemase-producing Enterobacteriaceae (CPE), and to determine epidemiological characteristics, laboratory methods and infection prevention and control (IPC) measures in acute care facilities. Methods A 6 month survey was undertaken in November 2013-April 2014 in 21 sentinel UK laboratories as part of the European Survey on Carbapenemase-Producing Enterobacteriaceae (EuSCAPE) project. Up to 10 consecutive, non-duplicate, clinically significant and carbapenem-non-susceptible isolates of Escherichia coli or Klebsiella pneumoniae were submitted to a reference laboratory. Participants answered a questionnaire on relevant laboratory methods and IPC measures. Results Of 102 isolates submitted, 89 (87%) were non-susceptible to ≥1 carbapenem, and 32 (36%) were confirmed as CPE. CPE were resistant to most antibiotics, except colistin (94% susceptible), gentamicin (63%), tigecycline (56%) and amikacin (53%). The prevalence of CPE was 0.02% (95% CI = 0.01%-0.03%). The incidence of CPE was 0.007 per 1000 patient-days (95% CI = 0.005-0.010), with north-west England the most affected region at 0.033 per 1000 patient-days (95% CI = 0.012-0.072). Recommended IPC measures were not universally followed, notably screening high-risk patients on admission (applied by 86%), using a CPE 'flag' on patients' records (70%) and alerting neighbouring hospitals when transferring affected patients (only 30%). Most sites (86%) had a laboratory protocol for CPE screening, most frequently using chromogenic agar (52%) or MacConkey/CLED agars with carbapenem discs (38%). Conclusions The UK prevalence and incidence of clinically significant CPE is currently low, but these MDR bacteria affect most UK regions. Improved IPC measures, vigilance and monitoring are required.

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. Electronic supplementary material The online version of this article (10.1186/s12864-018-4988-z) contains supplementary material, which is available to authorized users.

This deposit is composed by a publication in which the IGC's authors have had the role of collaboration (it's a collaboration publication). This type of deposit in ARCA is in restrictedAccess (it can't be in open access to the public), and can only be accessed by two ways: either by requesting a legal copy from the author (the email contact present in this deposit) or by visiting the following link: https://www.nature.com/articles/s41588-018-0074-3 This publication hasn't any creative commons license associated. This deposit is composed simultaneously by the original published article and also by the "correction" for the published article (erratum). The link for the original article: https://www.nature.com/articles/s41588-017-0029-0#Ack1 This deposit is composed by the main article plus the supplementary materials of the publication. In the version of this article initially published, the URL listed for TubercuList was incorrect. The correct URL is https://mycobrowser.epfl.ch/. The error has been corrected in the HTML and PDF versions of the article. To characterize the genetic determinants of resistance to antituberculosis drugs, we performed a genome-wide association study (GWAS) of 6,465 Mycobacterium tuberculosis clinical isolates from more than 30 countries. A GWAS approach within a mixed-regression framework was followed by a phylogenetics-based test for independent mutations. In addition to mutations in established and recently described resistance-associated genes, novel mutations were discovered for resistance to cycloserine, ethionamide and para-aminosalicylic acid. The capacity to detect mutations associated with resistance to ethionamide, pyrazinamide, capreomycin, cycloserine and para-aminosalicylic acid was enhanced by inclusion of insertions and deletions. Odds ratios for mutations within candidate genes were found to reflect levels of resistance. New epistatic relationships between candidate drug-resistance-associated genes were identified. Findings also suggest the involvement of efflux pumps (drrA and Rv2688c) in the emergence of resistance. This study will inform the design of new diagnostic tests and expedite the investigation of resistance and compensatory epistatic mechanisms. The project was supported by the KAUST faculty baseline research fund (BAS/1/1020-01-01) to A.P. The authors wish to thank members of the KAUST Bioscience Core laboratory who sequenced samples. We thank the Wellcome Trust Sanger Institute core and pathogen sequencing and informatics teams who were involved in the Malawi and Uganda studies. The work was funded in part by the Wellcome Trust (grant numbers WT096249/Z/11/B, WT088559MA, WT081814/Z/06/Z and WT098051) and the Wellcome Trust–Burroughs Wellcome Fund Infectious Diseases Initiative grant (number 063410/ABC/00/Z). F.C. was the recipient of a Bloomsbury College PhD Studentship and was supported by the Wellcome Trust (201344/Z/16/Z); J. Perdigão received a Fundação para a Ciência e a Tecnologia (Portugal) postdoctoral fellowship fund (SFRH/BPD/95406/2013). The Calouste Gulbenkian Foundation, the Institute Gulbenkian in Lisbon and the European Society of Clinical Microbiology and Infectious Diseases supported the research of C.P., J. Perdigão, I.P. and M.V. J. Phelan is funded by a BBSRC PhD studentship. T.G.C. is funded by the Medical Research Council UK (grant numbers MR/K000551/1, MR/M01360X/1, MR/N010469/1 and MC_PC_15103). N.F. is funded by the Medical Research Council UK (grant number MR/K020420/1). T.M. is supported by the Ministry of Health, Labor and Welfare of Japan (H21-Shinkou-Ippan-008 and H24-Shinkou-Ippan-010). info:eu-repo/semantics/publishedVersion

International audience; To characterize the genetic determinants of resistance to antituberculosis drugs, we performed a genome-wide association study (GWAS) of 6,465 $Mycobacterium\ tuberculosis$ clinical isolates from more than 30 countries. A GWAS approach within a mixed-regression framework was followed by a phylogenetics-based test for independent mutations. In addition to mutations in established and recently described resistance-associated genes, novel mutations were discovered for resistance to cycloserine, ethionamide and para-aminosalicylic acid. The capacity to detect mutations associated with resistance to ethionamide, pyrazinamide, capreomycin, cycloserine and para-aminosalicylic acid was enhanced by inclusion of insertions and deletions. Odds ratios for mutations within candidate genes were found to reflect levels of resistance. New epistatic relationships between candidate drug-resistance-associated genes were identified. Findings also suggest the involvement of efflux pumps ($drrA$ and $Rv2688c$) in the emergence of resistance. This study will inform the design of new diagnostic tests and expedite the investigation of resistance and compensatory epistatic mechanisms.

Background Mycobacterium tuberculosis (MTB) PE_PGRS genes belong to the PE multigene family. Although the function of PE_PGRS genes is unknown, it is hypothesized that the PE_PGRS genes may be associated with antigenic variability in MTB. Material and methods Whole genome sequencing analysis was performed on (n = 37) extensively drug-resistant (XDR) MTB strains from Pakistan, which included Lineage 1 (East African Indian, n = 2); Other lineage 1 (n = 3); Lineage 3 (Central Asian, n = 24); Other lineage 3 (n = 4); Lineage 4 (X3, n = 1) and T group (n = 3) MTB strains. Results There were 107 SNPs identified from the analysis of 42 PE_PGRS genes; of these, 13 were non-synonymous SNPs (nsSNPs). The nsSNPs identified in PE_PGRS genes – 6, 9 and 10 – were common in all EAI, CAS, Other lineages (1 and 3), T1 and X3. Deletions (DELs) in PE_PGRS genes – 3 and 19 – were observed in 17 (80.9%) CAS1 and 6 (85.7%) in Other lineages (1 and 3) XDR MTB strains, while DELs in the PE_PGRS49 were observed in all CAS1, CAS, CAS2 and Other lineages (1 and 3) XDR MTB strains. All CAS, EAI and Other lineages (1 and 3) strains showed insertions (INS) in PE_PGRS6 gene, while INS in the PE_PGRS genes 19 and 33 were observed in 20 (95.2%) CAS1, all CAS, CAS2, EAI and Other lineages (1 and 3) XDR MTB strains. Conclusion Genetic diversity in PE_PGRS genes contributes to antigenic variability and may result in increased immunogenicity of strains. This is the first study identifying variations in nsSNPs and INDELs in the PE_PGRS genes of XDR-TB strains from Pakistan. It highlights common genetic variations which may contribute to persistence.

Background Combating 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. Methods To 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. Results The 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. Conclusions Using 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. Electronic supplementary material The online version of this article (doi:10.1186/s12916-016-0575-9) contains supplementary material, which is available to authorized users.

The need for improved rapid diagnostic tests for tuberculosis disease has prompted interest in the volatile organic compounds (VOCs) emitted by Mycobacterium tuberculosis complex bacteria. We have investigated VOCs emitted by Mycobacterium bovis BCG grown on Lowenstein–Jensen media using selected ion flow tube mass spectrometry and thermal desorption-gas chromatography-mass spectrometry. Compounds observed included dimethyl sulphide, 3-methyl-1-butanol, 2-methyl-1-propanol, butanone, 2-methyl-1-butanol, methyl 2-methylbutanoate, 2-phenylethanol and hydrogen sulphide. Changes in levels of acetaldehyde, methanol and ammonia were also observed. The compounds identified are not unique to M. bovis BCG, and further studies are needed to validate their diagnostic value. Investigations using an ultra-rapid gas chromatograph with a surface acoustic wave sensor (zNose) demonstrated the presence of 2-phenylethanol (PEA) in the headspace of cultures of M. bovis BCG and Mycobacterium smegmatis, when grown on Lowenstein–Jensen supplemented with glycerol. PEA is a reversible inhibitor of DNA synthesis. It is used during selective isolation of gram-positive bacteria and may also be used to inhibit mycobacterial growth. PEA production was observed to be dependent on growth of mycobacteria. Further study is required to elucidate the metabolic pathways involved and assess whether this compound is produced during in vivo growth of mycobacteria.

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 IS 6110 -based restriction fragment length polymorphism (IS 6110 -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.

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.

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, eLife digest Tuberculosis is an important public health threat around the globe and is particularly common in developing countries. It is difficult to control the spread of the disease because the bacteria that cause it can spread when an infected individual coughs or sneezes. It may take years for an infected individual to develop symptoms of tuberculosis so it can be hard to trace the source of an outbreak, and people infected with HIV are particularly susceptible to the disease. The bacterium that causes the majority of cases of tuberculosis is called Mycobacterium tuberculosis. There are several different varieties or ‘lineages’ of M. tuberculosis, and it is thought that they may vary in their ability to spread and cause disease. However, the results of previous studies have been inconsistent and there also seems to be a lot of variation between strains within the same lineage. In this study, Guerra-Assunção et al. used an approach called whole genome sequencing alongside more traditional methods to study the spread of tuberculosis in Malawi. They sequenced the genomes of every available sample of M. tuberculosis collected from patients in the Karonga district of Malawi over a 15-year period. This produced high-quality DNA sequence data about the bacteria responsible for almost 1700 cases of disease. Using this massive amount of data, Guerra-Assunção et al. constructed networks that showed how the bacteria had spread in the community. This revealed that there were differences between the ability of the various M. tuberculosis lineages to cause disease and to spread in communities. For example, lineage 1 was less likely than the other lineages to cause disease soon after infecting an individual and was less able to spread. The data also show that the proportion of cases of disease due to recent infection declined substantially during the 15-year period. This indicates that the tuberculosis and HIV control programmes in the area have been successful. Guerra-Assunção et al.'s findings show that it is possible to understand how tuberculosis is transmitted on a large scale. The next challenge is to understand why the lineages differ in their ability to cause disease and spread between individuals. DOI: http://dx.doi.org/10.7554/eLife.05166.002

We studied genomic variation in a previously selected collection of isogenic Mycobacterium tuberculosis laboratory strains subjected to one or two rounds of antibiotic selection. Whole genome sequencing analysis identified eleven single, unique mutations (four synonymous, six non-synonymous, one intergenic), in addition to drug resistance-conferring mutations, that were fixed in the genomes of six monoresistant strains. Eight loci, present as minority variants (five non-synonymous, three synonymous) in the genome of the susceptible parent strain, became fixed in the genomes of multiple daughter strains. None of these mutations are known to be involved with drug resistance. Our results confirm previously observed genomic stability for M. tuberculosis, although the parent strain had accumulated allelic variants at multiple locations in an antibiotic-free in vitro environment. It is therefore likely to assume that these so-called hitchhiking mutations were co-selected and fixed in multiple daughter strains during antibiotic selection. The presence of multiple allelic variations, accumulated under non-selective conditions, which become fixed during subsequent selective steps, deserves attention. The wider availability of 'deep' sequencing methods could help to detect multiple bacterial (sub)populations within patients with high resolution and would therefore be useful in assisting in the detailed investigation of transmission chains.

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.

The emergence of difficult to treat multi drug-resistant tuberculosis currently threatens control efforts in all regions of the world. Genotyping of M. tuberculosis strains using molecular markers has the potential to identify outbreaks and track transmission of drug resistant strains. To be of value to the public health practitioner genotyping must be discriminatory, rapid and readily available. The preferred method of typing M. tuberculosis is MIRU-VNTR where 15 loci are examined by PCR. However, implementation requires the ability to accurately size the resulting amplicons. Analysis by traditional gel electrophoresis is tedious and slow and MIRU typing has so far been largely restricted to laboratories with sophisticated high throughput analytical facilities. ScreenTape (LAB901, Loanhead, Scotland) is a fully automated electophoretic device for rapid analysis of nucleic acids. It is simple to use and provides accurate sizing of amplicons in under 10 min. When applied to MIRU-VNTR analysis of 20 drug resistant M. tuberculosis isolates, the assignment of alleles was in agreement with results obtained from a commercial genotyping company on 97.1% occasions, compared to 93.4% for traditional gel electrophoresis. Correlation coefficients for predicted fragment length were 0.9911 and 0.9919 for ScreenTape and gel electrophoresis, respectively. Whereas the 20 isolates had previously been identified as having two spoligotype patterns, 14 MIRU patterns were observed demonstrating the superior discriminatory power of MIRU-VNTR typing for M. tuberculosis. ScreenTape appears to offer rapid discriminatory genotyping and may prove an attractive option for public health laboratories charged with controlling infectious tuberculosis disease.

Tests based on bacteriophage replication enable rapid screening of Mycobacterium tuberculosis for drug resistance. We describe a novel broth-based colorimetric method for detecting phage replication. When clinical isolates were tested by this novel method, high concordance was observed with both the traditional phage assay and gene mutation analysis for detection of resistance to rifampin.

Multidrug- (MDR) and extensively drug-resistant (XDR) tuberculosis (TB) present a challenge to disease control and elimination goals. Lisbon, Portugal, has a high TB incidence rate and unusual and successful XDR-TB strains that have been found in circulation for almost two decades. For the last 20 years, a continued circulation of two phylogenetic clades, Lisboa3 and Q1, which are highly associated with MDR and XDR, have been observed. In recent years, these strains have been well characterized regarding the molecular basis of drug resistance and have been inclusively subjected to whole genome sequencing (WGS). Researchers have been studying the genomic diversity of strains circulating in Lisbon and its genomic determinants through cutting-edge next generation sequencing. An enormous amount of whole genome sequence data are now available for the most prevalent and clinically relevant strains circulating in Lisbon. It is the persistence, prevalence and rapid evolution towards drug resistance that has prompted researchers to investigate the properties of these strains at the genomic level and in the future at a global transcriptomic level. Seventy Mycobacterium tuberculosis (MTB) isolates, mostly recovered in Lisbon, were genotyped by 24-loci Mycobacterial Interspersed Repetitive Unit – Variable Number of Tandem Repeats (MIRU-VNTR) and the genomes sequenced using a next generation sequencing platform – Illumina HiSeq 2000. The genotyping data revealed three major clusters associated with MDR-TB (Lisboa3-A, Lisboa3-B and Q1), two of which are associated with XDR-TB (Lisboa3-B and Q1), whilst the genomic data contributed to elucidating the phylogenetic positioning of circulating MDR-TB strains, showing a high predominance of a single SNP cluster group 5. Furthermore, a genome-wide phylogeny analysis from these strains, together with 19 publicly available genomes of MTB clinical isolates, revealed two major clades responsible for MDR/XDR-TB in the region: Lisboa3 and Q1. The data presented by this study contributes to the expanding knowledge of MTB genomic diversity yielding insights on microevolution and identification of novel compensatory mutations. Additionally, the analysis of non-synonymous/synonymous ratios revealed heterogeneities across the chromosome, genotype and Clusters of Orthologous Groups, highlighting possible and different evolution strategies. Overall, the results that are presented support the notion of an increasing genomic diversity that may support both setting and host adaptation.

Objectives The global increase in drug resistance in Mycobacterium tuberculosis (MTB) strains increases the focus on improved molecular diagnostics for MTB. Extensively drug-resistant (XDR) – TB is caused by MTB strains resistant to rifampicin, isoniazid, fluoroquinolone and aminoglycoside antibiotics. Resistance to anti-tuberculous drugs has been associated with single nucleotide polymorphisms (SNPs), in particular MTB genes. However, there is regional variation between MTB lineages and the SNPs associated with resistance. Therefore, there is a need to identify common resistance conferring SNPs so that effective molecular-based diagnostic tests for MTB can be developed. This study investigated used whole genome sequencing (WGS) to characterize 37 XDR MTB isolates from Pakistan and investigated SNPs related to drug resistance. Methods XDR-TB strains were selected. DNA was extracted from MTB strains, and samples underwent WGS with 76-base-paired end fragment sizes using Illumina paired end HiSeq2000 technology. Raw sequence data were mapped uniquely to H37Rv reference genome. The mappings allowed SNPs and small indels to be called using SAMtools/BCFtools. Results This study found that in all XDR strains, rifampicin resistance was attributable to SNPs in the rpo B RDR region. Isoniazid resistance-associated mutations were primarily related to kat G codon 315 followed by inh A S94A. Fluoroquinolone resistance was attributable to gyr A 91–94 codons in most strains, while one did not have SNPs in either gyr A or gyr B. Aminoglycoside resistance was mostly associated with SNPs in rrs , except in 6 strains. Ethambutol resistant strains had emb B codon 306 mutations, but many strains did not have this present. The SNPs were compared with those present in commercial assays such as LiPA Hain MDRTBsl, and the sensitivity of the assays for these strains was evaluated. Conclusions If common drug resistance associated with SNPs evaluated the concordance between phenotypic and genotypic testing, the results would be rifampicin (100%), isoniazid (89%), fluoroquinolones (95%), aminoglycoside (81%) and ethambutol (61%). This work highlights the importance of expanded targets for drug resistance detection in MTB isolates.

BACKGROUND: Multidrug- (MDR) and extensively drug resistant (XDR) tuberculosis (TB) presents a challenge to disease control and elimination goals. In Lisbon, Portugal, specific and successful XDR-TB strains have been found in circulation for almost two decades. RESULTS: In the present study we have genotyped and sequenced the genomes of 56 Mycobacterium tuberculosis isolates recovered mostly from Lisbon. The genotyping data revealed three major clusters associated with MDR-TB, two of which are associated with XDR-TB. Whilst the genomic data contributed to elucidate the phylogenetic positioning of circulating MDR-TB strains, showing a high predominance of a single SNP cluster group 5. Furthermore, a genome-wide phylogeny analysis from these strains, together with 19 publicly available genomes of Mycobacterium tuberculosis clinical isolates, revealed two major clades responsible for M/XDR-TB in the region: Lisboa3 and Q1 (LAM).The data presented by this study yielded insights on microevolution and identification of novel compensatory mutations associated with rifampicin resistance in rpoB and rpoC. The screening for other structural variations revealed putative clade-defining variants. One deletion in PPE41, found among Lisboa3 isolates, is proposed to contribute to immune evasion and as a selective advantage. Insertion sequence (IS) mapping has also demonstrated the role of IS6110 as a major driver in mycobacterial evolution by affecting gene integrity and regulation. CONCLUSIONS: Globally, this study contributes with novel genome-wide phylogenetic data and has led to the identification of new genomic variants that support the notion of a growing genomic diversity facing both setting and host adaptation. This work was partially supported by Project Ref. SDH49: “Early Molecular Detection of M/XDRTB in the Great Lisbon Healthcare Region” from Fundação Calouste Gulbenkian (FCG, Portugal) and PTDC/SAU-EPI/122400/ 2010 from Fundação para a Ciência e Tecnologia (FCT). The sequencing was funded by the KAUST Research Fund. J. Perdigão, D. Machado and C. Silva were supported by FCT grants SFRH/BPD/95406/2013, SFRH/BD/65060/2009 and SFRH/BD/73579/2010, respectively. TGC is funded by the Medical Research Council (UK) and Wellcome Trust.

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.

Background:\ud \ud Despite increasing interest in possible differences in virulence and transmissibility between different genotypes of M. tuberculosis, very little is known about how genotypes within a population change over decades, or about\ud relationships to HIV infection.\ud \ud Methods and Principal Findings:\ud \ud In a population-based study in rural Malawi we have examined smears and cultures from tuberculosis patients over a 20-year period using spoligotyping. Isolates were grouped into spoligotype families and lineages following previously published criteria. Time trends, HIV status, drug resistance and outcome were examined by spoligotype family and lineage. In addition, transmissibility was examined among pairs of cases with known epidemiological contact by assessing the proportion of transmissions confirmed for each lineage, on the basis of IS6110 RFLP similarity of the M tuberculosis strains. 760 spoligotypes were obtained from smears from 518 patients from 1986–2002, and 377 spoligotypes from cultures from 347 patients from 2005–2008. There was good consistency in patients with multiple specimens. Among 781 patients with first episode tuberculosis, the majority (76%) had Lineage 4 (‘‘European/American’’) strains; 9% had Lineage 3 (‘‘East-African/Indian’’); 8% Lineage 1 (‘‘Indo-Oceanic’’); and 2% Lineage 2 (‘‘East-Asian’’); others unclassifiable. Over time the proportion of Lineage 4 decreased from .90% to 60%, with an increase in the other 3 lineages (p,0.001). Lineage 1 strains were more common in those with HIV infection, even after adjusting for age, sex and year.\ud There were no associations with drug resistance or outcome, and no differences by lineage in the proportion of pairs in which transmission was confirmed.\ud \ud Conclusions: \ud \ud This is the first study to describe long term trends in the four M. tuberculosis lineages in a population. Lineage 4 has probably been longstanding in this population, with relatively recent introductions and spread of Lineages1–3, perhaps influenced by the HIV epidemic.

Summary Tuberculosis (TB) patients with strains common to other recent cases (‘clustering’) suggest recent transmission. HIV status and age may affect proportions clustered. We investigated TB clustering by HIV and age in a population-based study in Malawi. Among 746 patients, HIV infection increased the proportion clustered. Sex-period-adjusted odds ratios for the association of HIV and clustering were 1.26 (95% CI 0.4–4.1) for ages 15–25 years, 1.40 (0.9–2.3) for 25–50 years and 10.44 (2.3–47.9) for >50 years and remained stable over two periods examined. These results suggest that HIV increases the proportion of TB due to recent transmission in the elderly.

Background Resistance 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. Methods In 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. Results Initial 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$ when testing a batch of 20 isolates on a single 96 well plate. Results were obtained in 48 hours. Conclusion The 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.

Background Sudan 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. Methods In 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). Results Multi 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. Conclusions We 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.

According to a recent paper by Blackwood and colleagues (1), the Bacillus cereus and B. anthracis ATCC type strains share identical 16S rRNA gene sequences. From the study of a number of additional phylogenetic markers in 27 strains of the B. cereus group, they concluded that the rpoB gene sequence was the most discriminatory marker, as the B. cereus type strain revealed four mismatches at gene positions considered key positions for B. anthracis (5). We fully agree with the notion of the importance of a polyphasic approach to species identification among the Bacillus cereus group. However, in our experience rpoB may not be specific enough for the identification of B. anthracis and B. cereus to the species level. (This research was presented as an abstract at the 103rd General Meeting of the American Society for Microbiology in Washington, D.C., 2003.) We isolated a Ba813-positive member of the B. cereus group with the phenotypic markers of B. cereus. A 19-year-old male developed wound infection and septic shock after surgery because of femur osteosarcoma. A Bacillus sp. was isolated from blood cultures and wound swabs in pure culture. Grey, dry, beta-hemolytic colonies grew on 5% sheep blood agar. The gram-positive, rod-shaped bacteria showed ellipsoidal central and subterminal spores. The isolate ifik1001803 was motile, resistant to penicillin G, and catalase positive. Analysis of the cellular fatty acids (MIDI system; Microbial ID, Inc., Newark, Del.) was consistent with B. cereus (without 17:1 anteiso A; with ISO 17:1 w10c) (personal communication of R. Paisley, MIDI, Inc., Newark, Del.). Sequence analysis of the 16S rRNA full gene showed a 100% match with B. anthracis type 6 (6). PCR amplifications of two virulence factors of B. anthracis (2), the protective antigen gene pag and the capsule gene fragment capC, were negative. PCR amplification of the Ba813 chromosomal fragment (4) gave a band of the expected size; sequencing of the PCR product revealed a 5-bp mismatch to B. anthracis ({"type":"entrez-nucleotide","attrs":{"text":"U46157","term_id":"1389560","term_text":"U46157"}}U46157). In addition, the sequence of the rpoB gene of ifik1001803 was identical to that of B. anthracis at three of the four key positions, but it differed from it by 2 bp (Table ​(Table1).1). The variability of the rpoB gene sequence among B. cereus strains and ifik1001803 was in the range of 3 to 5 bp. Interestingly, B. cereus ATCC 49064, the B. cereus strain most closely related to ifik1001803, had a match with B. anthracis at one of the four key positions. From the perspective of sequence information, ifik1001803 might be viewed as a transitional composite of B. cereus and B. anthracis, supporting the recent concept of B. anthracis being a lineage of B. cereus (3). In this respect, B. cereus ATCC 49064 might represent another transitional member of the B. cereus group. TABLE 1. Alignment of B. cereus and B. anthracis sequence information of the rpoB gene Our results suggest that position 84 of the rpoB gene may represent the only key position for the differentiation of B. cereus from B. anthracis. However, the number of strains examined as well the choice of strains studied may influence the outcome of phylogenetic studies. In addition, a true polyphasic approach would require the consideration of independent information, such as morphological and physiological characteristics, including pathogenicity.

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. Bloomsbury Research Fund; Fundação para a Ciência e Tecnologia (grant number PTDC/SAU-ESA/71208/2006); King Abdullah University of Science and Technology (KAUST); Medical Research Council UK (grant number MR/K000551/1); Wellcome Trust (grant number 096249/Z/11/B)

Background: Mycobacterium tuberculosis (MTB) PE_PGRS genes belong to the PE multi-gene family. Although the function of the members of the PE_PGRS multi-gene family is not yet known, it is hypothesized that the PE_PGRS genes may be associated with genetic variability. Material and Methods: Whole genome sequencing analysis was performed on (n = 37) extensively drug resistant (XDR) MTB strains from Pakistan which included Central Asian (n = 23), East African Indian (n = 2), X3 (n = 1), T group (n = 3) and Orphan (n = 8) MTB strains. Results: By analyzing 42 PE_PGRS genes, 111 SNPs were identified, of which 13 were non-synonymous SNPs (nsSNPs). The nsSNPs identified in the PE_PGRS genes were as follows: 6, 9, 10 and 55 present in each of the CAS, EAI, Orphan, T1 and X3 XDR MTB strains studied. Deletions in PE_PGRS genes: 19, 21 and 23 were observed in 7 (35.0%) CAS1 and 3 (37.5%) in Orphan XDR MTB strains, while deletions in the PE_PGRS genes: 49 and 50 were observed in 36 (95.0%) CAS1 and all CAS, CAS2 and Orphan XDR MTB strains. An insertion in PE_PGRS6 gene was observed in all CAS, EAI3 and Orphan, while insertions in the PE_PGRS genes 19 and 33 were observed in 19 (95%) CAS1 and all CAS, CAS2, EAI3 and Orphan XDR MTB strains. Conclusion: Genetic diversity in PE_PGRS genes contributes to antigenic variability and may result in increased immunogenicity of strains. This is the first study identifying variations in nsSNPs, Insertions and Deletions in the PE_PGRS genes of XDR-TB strains from Pakistan. It highlights common genetic variations which may contribute to persistence.