Whole genome sequencing versus traditional genotyping for investigation of a Mycobacterium tuberculosis outbreak: a longitudinal molecular epidemiological study

Integration of Whole-Genome Sequencing into Infection Control Practices: the Potential and the Hurdles

Microbial whole-genome sequencing (WGS) is poised to transform many of the currently used approaches in medical microbiology. Recent reports on the application of WGS to understand genetic evolution and reconstruct transmission pathways have provided valuable information that will influence infection control practices. While this technology holds great promise, obstacles to full implementation remain. Two articles in this issue of the Journal of Clinical Microbiology (S. Octavia, Q. Wang, M. M. Tanaka, S. Kaur, V. Sintchenko, and R. Lan, J Clin Microbiol 53:1063–1071, 2015, doi:10.1128/JCM.03235-14, andS. J. Salipante, D. J. SenGupta, L. A. Cummings, T. A. Land, D. R. Hoogestraat, and B. T. Cookson, J Clin Microbiol 53:1072–1079, 2015, doi:10.1128/JCM.03385-14) describe the breadth of application of WGS to the field of clinical epidemiology.

Mycobacterial DNA extraction for whole-genome sequencing from early positive liquid (MGIT) cultures

We developed a low-cost and reliable method of DNA extraction from as little as 1 ml of early positive mycobacterial growth indicator tube (MGIT) cultures that is suitable for whole-genome sequencing to identify mycobacterial species and predict antibiotic resistance in clinical samples. The DNA extraction method is based on ethanol precipitation supplemented by pretreatment steps with a MolYsis kit or saline wash for the removal of human DNA and a final DNA cleanup step with solid-phase reversible immobilization beads. The protocol yielded ≥0.2 ng/μl of DNA for 90% (MolYsis kit) and 83% (saline wash) of positive MGIT cultures. A total of 144 (94%) of the 154 samples sequenced on the MiSeq platform (Illumina) achieved the target of 1 million reads, with <5% of reads derived from human or nasopharyngeal flora for 88% and 91% of samples, respectively. A total of 59 (98%) of 60 samples that were identified by the national mycobacterial reference laboratory (NMRL) as Mycobacterium tuberculosis were successfully mapped to the H37Rv reference, with >90% coverage achieved. The DNA extraction protocol, therefore, will facilitate fast and accurate identification of mycobacterial species and resistance using a range of bioinformatics tools.

Pathogen genome bioinformatics

Recent advances in DNA sequencing technology have made the whole-genome sequencing of pathogens in a clinically relevant turn-around time both technically and economically feasible. The DNA sequencing of pathogens with epidemic potential offers new and exciting opportunities for high-resolution public health surveillance. This chapter outlines major methods and bioinformatics tools for pathogen genome characterization, the identification of infectious disease clusters, as well as for genomics-guided biosurveillance. Existing challenges are also considered.

Reemergence and amplification of tuberculosis in the Canadian arctic

BACKGROUND

Between November 2011 and November 2012, a Canadian village of 933 persons had 50 culture-positive cases of tuberculosis, with 49 sharing the same genotype.

METHODS

We performed Illumina-based whole-genome sequencing on Mycobacterium tuberculosis isolates from this village, during and before the outbreak. Phylogenetic trees were generated using the maximum likelihood method.

RESULTS

Three distinct genotypes were identified. Strain I (n = 7) was isolated in 1991-1996. Strain II (n = 8) was isolated in 1996-2004. Strain III (n = 62) first appeared in 2007 and did not arise from strain I or II. Within strain III, there were 3 related but distinct clusters: IIIA, IIIB, and IIIC. Between 2007 and 2010, cluster IIIA predominated (11 of 22 vs 2 of 40; P < .001), whereas in 2011-2012 clusters IIIB (n = 18) and IIIC (n = 20) predominated over cluster IIIA (n = 11). Combined evolutionary and epidemiologic analysis of strain III cases revealed that the outbreak in 2011-2012 was the result of ≥6 temporally staggered events, spanning from 1 reactivation case to a point-source outbreak of 20 cases.

CONCLUSIONS

After the disappearance of 2 strains of M. tuberculosis in this village, its reemergence in 2007 was followed by an epidemiologic amplification, affecting >5% of the population.

Mathematical Modelling and Tuberculosis: Advances in Diagnostics and Novel Therapies

As novel diagnostics, therapies, and algorithms are developed to improve case finding, diagnosis, and clinical management of patients with TB, policymakers must make difficult decisions and choose among multiple new technologies while operating under heavy resource constrained settings. Mathematical modelling can provide helpful insight by describing the types of interventions likely to maximize impact on the population level and highlighting those gaps in our current knowledge that are most important for making such assessments. This review discusses the major contributions of TB transmission models in general, namely, the ability to improve our understanding of the epidemiology of TB. We focus particularly on those elements that are important to appropriately understand the role of TB diagnosis and treatment (i.e., what elements of better diagnosis or treatment are likely to have greatest population-level impact) and yet remain poorly understood at present. It is essential for modellers, decision-makers, and epidemiologists alike to recognize these outstanding gaps in knowledge and understand their potential influence on model projections that may guide critical policy choices (e.g., investment and scale-up decisions).

Current and emerging Legionella diagnostics for laboratory and outbreak investigations

Legionnaires' disease (LD) is an often severe and potentially fatal form of bacterial pneumonia caused by an extensive list of Legionella species. These ubiquitous freshwater and soil inhabitants cause human respiratory disease when amplified in man-made water or cooling systems and their aerosols expose a susceptible population. Treatment of sporadic cases and rapid control of LD outbreaks benefit from swift diagnosis in concert with discriminatory bacterial typing for immediate epidemiological responses. Traditional culture and serology were instrumental in describing disease incidence early in its history; currently, diagnosis of LD relies almost solely on the urinary antigen test, which captures only the dominant species and serogroup, Legionella pneumophila serogroup 1 (Lp1). This has created a diagnostic "blind spot" for LD caused by non-Lp1 strains. This review focuses on historic, current, and emerging technologies that hold promise for increasing LD diagnostic efficiency and detection rates as part of a coherent testing regimen. The importance of cooperation between epidemiologists and laboratorians for a rapid outbreak response is also illustrated in field investigations conducted by the CDC with state and local authorities. Finally, challenges facing health care professionals, building managers, and the public health community in combating LD are highlighted, and potential solutions are discussed.

Whole genome sequencing: a new paradigm in the surveillance and control of human tuberculosis

Whole Genome Sequencing (WGS) is emerging as a very powerful tool for the management, outbreak analyses, surveillance and determining drug resistance of human infectious pathogens including Mycobacterium tuberculosis and MRSA. WGS can also discriminate relapse TB from re-infection and the resolution provided by WGS has no comparison to conventional technologies. With current cost coming down to <£70 per bacterial genome, WGS has emerged as an alternative to all the existing technologies put together. We discuss the advantage and disadvantages of WGS and whether it can become a point of care tool in not just developed countries but also in developing countries which have a huge TB burden. The likely utility of WGS for other pathogens and also in characterizing holobionts is also discussed.

[High-throughput sequencing: towards a genome-based diagnosis in infectious diseases]

During a pathogen outbreak, the emergency resides in the identification and characterization of the infectious agent. In addition to the traditional phenotypic methods which are still widely used, the molecular biology is nowadays a common approach of clinical microbiology labs and the pathogen can be identified by comparing its molecular fingerprint to a data-bank. High-throughput sequencing should allow overcoming this single identification to exploit the whole information encoded in the pathogen genome. This evolution, supported by an increasing number of proof-of-concept studies, should result in moving from detection through fingerprints to the use of the pathogen whole genome; this forensic profile should allow the adaptation of the treatment to the pathogen specificities. From concept to routine use, many parameters need to be considered to promote high-throughput sequencing as a powerful tool to help physicians and clinicians in microbiological investigations.

Dissemination of cephalosporin resistance genes between Escherichia coli strains from farm animals and humans by specific plasmid lineages

Third-generation cephalosporins are a class of β-lactam antibiotics that are often used for the treatment of human infections caused by Gram-negative bacteria, especially Escherichia coli. Worryingly, the incidence of human infections caused by third-generation cephalosporin-resistant E. coli is increasing worldwide. Recent studies have suggested that these E. coli strains, and their antibiotic resistance genes, can spread from food-producing animals, via the food-chain, to humans. However, these studies used traditional typing methods, which may not have provided sufficient resolution to reliably assess the relatedness of these strains. We therefore used whole-genome sequencing (WGS) to study the relatedness of cephalosporin-resistant E. coli from humans, chicken meat, poultry and pigs. One strain collection included pairs of human and poultry-associated strains that had previously been considered to be identical based on Multi-Locus Sequence Typing, plasmid typing and antibiotic resistance gene sequencing. The second collection included isolates from farmers and their pigs. WGS analysis revealed considerable heterogeneity between human and poultry-associated isolates. The most closely related pairs of strains from both sources carried 1263 Single-Nucleotide Polymorphisms (SNPs) per Mbp core genome. In contrast, epidemiologically linked strains from humans and pigs differed by only 1.8 SNPs per Mbp core genome. WGS-based plasmid reconstructions revealed three distinct plasmid lineages (IncI1- and IncK-type) that carried cephalosporin resistance genes of the Extended-Spectrum Beta-Lactamase (ESBL)- and AmpC-types. The plasmid backbones within each lineage were virtually identical and were shared by genetically unrelated human and animal isolates. Plasmid reconstructions from short-read sequencing data were validated by long-read DNA sequencing for two strains. Our findings failed to demonstrate evidence for recent clonal transmission of cephalosporin-resistant E. coli strains from poultry to humans, as has been suggested based on traditional, low-resolution typing methods. Instead, our data suggest that cephalosporin resistance genes are mainly disseminated in animals and humans via distinct plasmids.

The rapid global rise of infections caused by Escherichia coli that are resistant to clinically relevant antimicrobials, including third-generation cephalosporins, is cause for concern. The intestinal tract of livestock, in particular poultry, is an important reservoir for drug resistant E. coli, but it is unknown to what extent these bacteria can spread to humans. Food is thought to be an important source because drug-resistant E. coli have been detected in animals raised for meat consumption and in meat products. Previous studies that used traditional, low-resolution, genetic typing methods found that drug resistant E. coli present in humans and poultry were indistinguishable from each other, suggesting dissemination of these bacteria through the food-chain to humans. However, by applying high-resolution, whole-genome sequencing methods, we did not find evidence for such transmission of bacteria through the food-chain. Instead, by employing a novel approach for the reconstruction of mobile genetic elements from whole-genome sequence data, we discovered that genetically unrelated E. coli isolates from both humans and animal sources carried nearly identical plasmids that encode third-generation cephalosporin resistance determinants. Our data suggest that cephalosporin resistance is mainly disseminated via the transfer of mobile genetic elements between animals and humans.

Added value of whole-genome sequencing for management of highly drug-resistant TB

Objectives

Phenotypic drug susceptibility testing (DST) for Mycobacterium tuberculosis takes several weeks to complete and second-line DST is often poorly reproducible, potentially leading to compromised clinical decisions. Following a fatal case of XDR TB, we investigated the potential benefit of using whole-genome sequencing to generate an in silico drug susceptibility profile.

Methods

The clinical course of the patient was reviewed, assessing the times at which phenotypic DST data became available and changes made to the therapeutic regimen. Whole-genome sequencing was performed on the earliest available isolate and variants associated with drug resistance were identified.

Results

The final DST report, including second-line drugs, was issued 10 weeks after patient presentation and 8 weeks after initial growth of M. tuberculosis. In the interim, the patient may have received a compromised regimen that had the potential to select for further drug resistance. The in silico susceptibility profile, extrapolated from evolving evidence in the literature, provided comparable or superior data to the DST results for second-line drugs and could be generated in a much shorter timeframe.

Conclusions

We propose routine whole-genome sequencing of all MDR M. tuberculosis isolates in adequately resourced settings. This will improve individual patient care, monitor for transmission events and advance our understanding of resistance-associated mutations.

Clinical detection and characterization of bacterial pathogens in the genomics era

The availability of genome sequences obtained using next-generation sequencing (NGS) has revolutionized the field of infectious diseases. Indeed, more than 38,000 bacterial and 5,000 viral genomes have been sequenced to date, including representatives of all significant human pathogens. These tremendous amounts of data have not only enabled advances in fundamental biology, helping to understand the pathogenesis of microorganisms and their genomic evolution, but have also had implications for clinical microbiology. Here, we first review the current achievements of genomics in the development of improved diagnostic tools, including those that are now available in the clinic, such as the design of PCR assays for the detection of microbial pathogens, virulence factors or antibiotic-resistance determinants, or the design of optimized culture media for 'unculturable' pathogens. We then review the applications of genomics to the investigation of outbreaks, either through the design of genotyping assays or the direct sequencing of the causative strains. Finally, we discuss how genomics might change clinical microbiology in the future.

Marked microevolution of a unique Mycobacterium tuberculosis strain in 17 years of ongoing transmission in a high risk population

The transmission and persistence of Mycobacterium tuberculosis within high risk populations is a threat to tuberculosis (TB) control. In the current study, we used whole genome sequencing (WGS) to decipher the transmission dynamics and microevolution of M. tuberculosis ON-A, an endemic strain responsible for an ongoing outbreak of TB in an urban homeless/under-housed population. Sixty-one M. tuberculosis isolates representing 57 TB cases from 1997 to 2013 were subjected to WGS. Sequencing data was integrated with available epidemiological information and analyzed to determine how the M. tuberculosis ON-A strain has evolved during almost two decades of active transmission. WGS offers higher discriminatory power than traditional genotyping techniques, dividing the M. tuberculosis ON-A strain into 6 sub-clusters, each defined by unique single nucleotide polymorphism profiles. One sub-cluster, designated ON-A^NM (Natural Mutant; 26 isolates from 24 cases) was also defined by a large, 15 kb genomic deletion. WGS analysis reveals the existence of multiple transmission chains within the same population/setting. Our results help validate the utility of WGS as a powerful tool for identifying genomic changes and adaptation of M. tuberculosis.

Evolution of extensively drug-resistant Mycobacterium tuberculosis from a susceptible ancestor in a single patient

Background

Mycobacterium tuberculosis is characterized by a low mutation rate and a lack of genetic recombination. Yet, the rise of extensively resistant strains paints a picture of a microbe with an impressive adaptive potential. Here we describe the first documented case of extensively drug-resistant tuberculosis evolved from a susceptible ancestor within a single patient.

Results

Genome sequences of nine serial M. tuberculosis isolates from the same patient uncovered a dramatic turnover of competing lineages driven by the emergence, and subsequent fixation or loss of single nucleotide polymorphisms. For most drugs, resistance arose through independent emergence of mutations in more than one clone, of which only one ultimately prevailed as the clone carrying it expanded, displacing the other clones in the process. The vast majority of mutations identified over 3.5 years were either involved in drug resistance or hitchhiking in the genetic background of these. Additionally, RNA-sequencing of isolates grown in the absence of drug challenge revealed that the efflux-associated iniBAC operon was up-regulated over time, whereas down-regulated genes include those involved in mycolic acid synthesis.

Conclusions

We observed both rapid acquisitions of resistance to antimicrobial compounds mediated by individual mutations as well as a gradual increase in fitness in the presence of antibiotics, likely driven by stable gene expression reprogramming. The rapid turnover of resistance mutations and hitchhiking neutral mutations has major implications for inferring tuberculosis transmission events in situations where drug resistance evolves within transmission chains.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-014-0490-3) contains supplementary material, which is available to authorized users.

Evaluation of large genotypic Mycobacterium tuberculosis clusters: contributions from remote and recent transmission

Tuberculosis genotypic clustering is used as a proxy for recent transmission. The association between clustering and recent transmission becomes problematic when the genotyping method lacks specificity in defining a cluster, as well as for clusters with extensive jurisdictional histories and/or common genotypes. We investigated the four largest spoligotype/12 loci MIRU-VNTR-defined clusters in Harris County, Texas from 2006-2012 to determine their historical contribution to tuberculosis morbidity, estimate the contributions from recent and remote transmission, and determine the impact of secondary genotyping on cluster definition. The clusters contained 189, 64, 51 and 38 cases. Each cluster was linked to cluster(s) previously identified by Houston Tuberculosis Initiative; 3 since 1995 and the fourth in 2002. Among cases for which timing of Mycobacterium tuberculosis transmission relative to tuberculosis disease could be ascertained, nearly equal proportions were associated with recent and remote transmission. The extent to which genotyping with an additional 12 MIRU-VNTR loci modified the cluster definition varied from little or no impact for the two smaller clusters to moderate impact for the larger clusters. Tuberculosis control measures to reduce morbidity associated with large clusters must involve strategies to identify and treat individuals who recently acquired infection, as well as persons infected for years.

Tracking a tuberculosis outbreak over 21 years: strain-specific single-nucleotide polymorphism typing combined with targeted whole-genome sequencing

BACKGROUND

Whole-genome sequencing (WGS) is increasingly used in molecular-epidemiological investigations of bacterial pathogens, despite cost- and time-intensive analyses. We combined strain-specific single-nucleotide polymorphism (SNP) typing and targeted WGS to investigate a tuberculosis cluster spanning 21 years in Bern, Switzerland.

METHODS

On the basis of genome sequences of 3 historical outbreak Mycobacterium tuberculosis isolates, we developed a strain-specific SNP-typing assay to identify further cases. We screened 1642 patient isolates and performed WGS on all identified cluster isolates. We extracted SNPs to construct genomic networks. Clinical and social data were retrospectively collected.

RESULTS

We identified 68 patients associated with the outbreak strain. Most received a tuberculosis diagnosis in 1991-1995, but cases were observed until 2011. Two thirds were homeless and/or substance abusers. Targeted WGS revealed 133 variable SNP positions among outbreak isolates. Genomic network analyses suggested a single origin of the outbreak, with subsequent division into 3 subclusters. Isolates from patients with confirmed epidemiological links differed by 0-11 SNPs.

CONCLUSIONS

Strain-specific SNP genotyping allowed rapid and inexpensive identification of M. tuberculosis outbreak isolates in a population-based strain collection. Subsequent targeted WGS provided detailed insights into transmission dynamics. This combined approach could be applied to track bacterial pathogens in real time and at high resolution.

Age-specific risks of tuberculosis infection from household and community exposures and opportunities for interventions in a high-burden setting

We analyzed data from a large population-based prospective cohort study of household contacts of tuberculosis patients in Lima, Peru, to estimate the importance of within-household transmission relative to community-based transmission. We identified all adults (older than 15 years of age) who had incident pulmonary tuberculosis diagnosed at any of 106 public health centers in Lima from September 2009 to August 2012. A total of 14,041 household contacts of 3,446 index patients were assessed for tuberculosis infection and disease. We compared the prevalence of latent tuberculosis infection (LTBI) among persons who had received the Bacillus Calmette-Guérin vaccine in households with and without a microbiologically confirmed index case to estimate the age-specific risk of infection and excess risk of LTBI from household and community exposures. We found that the risk of infection from household and community sources increased from birth until 20 years of age. However, a large proportion of infections among child and young-adult household contacts could have been the result of household exposure. Excess infection risk associated with household exposure accounted for 58% (95% confidence interval: 47, 66) of LTBI prevalence among exposed children younger than 1 year of age, 48% (95% confidence interval: 39, 57) among 10-year-old children, and 44% (95% confidence interval: 34, 51) among 15-year-old adolescents. These findings suggest that expanded access to preventive therapy for older children and young-adult household contacts of known tuberculosis cases may be beneficial.

Standing of nucleic acid testing strategies in veterinary diagnosis laboratories to uncover Mycobacterium tuberculosis complex members

Nucleic acid testing (NAT) designate any molecular approach used for the detection, identification, and characterization of pathogenic microorganisms, enabling the rapid, specific, and sensitive diagnostic of infectious diseases, such as tuberculosis. These assays have been widely used since the 90s of the last century in human clinical laboratories and, subsequently, also in veterinary diagnostics. Most NAT strategies are based in the polymerase chain reaction (PCR) and its several enhancements and variations. From the conventional PCR, real-time PCR and its combinations, isothermal DNA amplification, to the nanotechnologies, here we review how the NAT assays have been applied to decipher if and which member of the Mycobacterium tuberculosis complex is present in a clinical sample. Recent advances in DNA sequencing also brought new challenges and have made possible to generate rapidly and at a low cost, large amounts of sequence data. This revolution with the high-throughput sequencing (HTS) technologies makes whole genome sequencing (WGS) and metagenomics the trendiest NAT strategies, today. The ranking of NAT techniques in the field of clinical diagnostics is rising, and we provide a SWOT (Strengths, Weaknesses, Opportunities, and Threats) analysis with our view of the use of molecular diagnostics for detecting tuberculosis in veterinary laboratories, notwithstanding the gold standard being still the classical culture of the agent. The complementary use of both classical and molecular diagnostics approaches is recommended to speed the diagnostic, enabling a fast decision by competent authorities and rapid tackling of the disease.

Prospective genotyping of Mycobacterium tuberculosis from fresh clinical samples

Shorter time-to-result is key for improving molecular-guided epidemiological investigation of tuberculosis (TB) cases. We performed a prospective study to evaluate the use of standardized MIRU-VNTR (mycobacterial interspersed repetitive-unit-variable-number tandem-repeat) typing of Mycobacterium tuberculosis directly on 79 fresh clinical samples from 26 TB patients consecutively enrolled over a 17-month period. Overall, complete 24-locus types were obtained for 18 out of the 26 (69.2%) patients and 14 of the 16 grade 3+ and grade 2+ samples (87.5%). The degree of completion of the genotypes obtained significantly correlated with smear microscopy grade both for 26 first samples (p = 0.0003) and for 53 follow-up samples (p = 0.002). For 20 of the 26 patients for whom complete or even incomplete M. tuberculosis isolate genotypes were obtained, typing applied to the clinical samples allowed the same unambiguous conclusions regarding case clustering or uniqueness as those that could have been drawn based on the corresponding cultured isolates. Standard 24 locus MIRU-VNTR typing of M. tuberculosis can be applied directly to fresh clinical samples, with typeability depending on the bacterial load in the sample.

Strain classification of Mycobacterium tuberculosis isolates in Brazil based on genotypes obtained by spoligotyping, mycobacterial interspersed repetitive unit typing and the presence of large sequence and single nucleotide polymorphism

Rio de Janeiro is endemic for tuberculosis (TB) and presents the second largest prevalence of the disease in Brazil. Here, we present the bacterial population structure of 218 isolates of Mycobacterium tuberculosis, derived from 186 patients that were diagnosed between January 2008 and December 2009. Genotypes were generated by means of spoligotyping, 24 MIRU-VNTR typing and presence of fbpC¹⁰³, RD^Rio and RD174. The results confirmed earlier data that predominant genotypes in Rio de Janeiro are those of the Euro American Lineages (99%). However, we observed differences between the classification by spoligotyping when comparing to that of 24 MIRU-VNTR typing, being respectively 43.6% vs. 62.4% of LAM, 34.9% vs. 9.6% of T and 18.3% vs. 21.5% of Haarlem. Among isolates classified as LAM by MIRU typing, 28.0% did not present the characteristic spoligotype profile with absence of spacers 21 to 24 and 32 to 36 and we designated these conveniently as “LAM-like”, 79.3% of these presenting the LAM-specific SNP fbpC¹⁰³. The frequency of RD^Rio and RD174 in the LAM strains, as defined both by spoligotyping and 24 MIRU-VNTR loci, were respectively 11% and 15.4%, demonstrating that RD174 is not always a marker for LAM/RD^Rio strains. We conclude that, although spoligotyping alone is a tool for classification of strains of the Euro-American lineage, when combined with MIRU-VNTRs, SNPs and RD typing, it leads to a much better understanding of the bacterial population structure and phylogenetic relationships among strains of M. tuberculosis in regions with high incidence of TB.

An evaluation of alternative methods for constructing phylogenies from whole genome sequence data: a case study with Salmonella

Comparative genomics based on whole genome sequencing (WGS) is increasingly being applied to investigate questions within evolutionary and molecular biology, as well as questions concerning public health (e.g., pathogen outbreaks). Given the impact that conclusions derived from such analyses may have, we have evaluated the robustness of clustering individuals based on WGS data to three key factors: (1) next-generation sequencing (NGS) platform (HiSeq, MiSeq, IonTorrent, 454, and SOLiD), (2) algorithms used to construct a SNP (single nucleotide polymorphism) matrix (reference-based and reference-free), and (3) phylogenetic inference method (FastTreeMP, GARLI, and RAxML). We carried out these analyses on 194 whole genome sequences representing 107 unique Salmonella enterica subsp. enterica ser. Montevideo strains. Reference-based approaches for identifying SNPs produced trees that were significantly more similar to one another than those produced under the reference-free approach. Topologies inferred using a core matrix (i.e., no missing data) were significantly more discordant than those inferred using a non-core matrix that allows for some missing data. However, allowing for too much missing data likely results in a high false discovery rate of SNPs. When analyzing the same SNP matrix, we observed that the more thorough inference methods implemented in GARLI and RAxML produced more similar topologies than FastTreeMP. Our results also confirm that reproducibility varies among NGS platforms where the MiSeq had the lowest number of pairwise differences among replicate runs. Our investigation into the robustness of clustering patterns illustrates the importance of carefully considering how data from different platforms are combined and analyzed. We found clear differences in the topologies inferred, and certain methods performed significantly better than others for discriminating between the highly clonal organisms investigated here. The methods supported by our results represent a preliminary set of guidelines and a step towards developing validated standards for clustering based on whole genome sequence data.

KvarQ: targeted and direct variant calling from fastq reads of bacterial genomes

BACKGROUND

High-throughput DNA sequencing produces vast amounts of data, with millions of short reads that usually have to be mapped to a reference genome or newly assembled. Both reference-based mapping and de novo assembly are computationally intensive, generating large intermediary data files, and thus require bioinformatics skills that are often lacking in the laboratories producing the data. Moreover, many research and practical applications in microbiology require only a small fraction of the whole genome data.

RESULTS

We developed KvarQ, a new tool that directly scans fastq files of bacterial genome sequences for known variants, such as single nucleotide polymorphisms (SNP), bypassing the need of mapping all sequencing reads to a reference genome and de novo assembly. Instead, KvarQ loads "testsuites" that define specific SNPs or short regions of interest in a reference genome, and directly synthesizes the relevant results based on the occurrence of these markers in the fastq files. KvarQ has a versatile command line interface and a graphical user interface. KvarQ currently ships with two "testsuites" for Mycobacterium tuberculosis, but new "testsuites" for other organisms can easily be created and distributed. In this article, we demonstrate how KvarQ can be used to successfully detect all main drug resistance mutations and phylogenetic markers in 880 bacterial whole genome sequences. The average scanning time per genome sequence was two minutes. The variant calls of a subset of these genomes were validated with a standard bioinformatics pipeline and revealed >99% congruency.

CONCLUSION

KvarQ is a user-friendly tool that directly extracts relevant information from fastq files. This enables researchers and laboratory technicians with limited bioinformatics expertise to scan and analyze raw sequencing data in a matter of minutes. KvarQ is open-source, and pre-compiled packages with a graphical user interface are available at http://www.swisstph.ch/kvarq.

Culture-independent detection and characterisation of Mycobacterium tuberculosis and M. africanum in sputum samples using shotgun metagenomics on a benchtop sequencer

Tuberculosis remains a major global health problem. Laboratory diagnostic methods that allow effective, early detection of cases are central to management of tuberculosis in the individual patient and in the community. Since the 1880s, laboratory diagnosis of tuberculosis has relied primarily on microscopy and culture. However, microscopy fails to provide species- or lineage-level identification and culture-based workflows for diagnosis of tuberculosis remain complex, expensive, slow, technically demanding and poorly able to handle mixed infections. We therefore explored the potential of shotgun metagenomics, sequencing of DNA from samples without culture or target-specific amplification or capture, to detect and characterise strains from the Mycobacterium tuberculosis complex in smear-positive sputum samples obtained from The Gambia in West Africa. Eight smear- and culture-positive sputum samples were investigated using a differential-lysis protocol followed by a kit-based DNA extraction method, with sequencing performed on a benchtop sequencing instrument, the Illumina MiSeq. The number of sequence reads in each sputum-derived metagenome ranged from 989,442 to 2,818,238. The proportion of reads in each metagenome mapping against the human genome ranged from 20% to 99%. We were able to detect sequences from the M. tuberculosis complex in all eight samples, with coverage of the H37Rv reference genome ranging from 0.002X to 0.7X. By analysing the distribution of large sequence polymorphisms (deletions and the locations of the insertion element IS6110) and single nucleotide polymorphisms (SNPs), we were able to assign seven of eight metagenome-derived genomes to a species and lineage within the M. tuberculosis complex. Two metagenome-derived mycobacterial genomes were assigned to M. africanum, a species largely confined to West Africa; the others that could be assigned belonged to lineages T, H or LAM within the clade of "modern" M. tuberculosis strains. We have provided proof of principle that shotgun metagenomics can be used to detect and characterise M. tuberculosis sequences from sputum samples without culture or target-specific amplification or capture, using an accessible benchtop-sequencing platform, the Illumina MiSeq, and relatively simple DNA extraction, sequencing and bioinformatics protocols. In our hands, sputum metagenomics does not yet deliver sufficient depth of coverage to allow sequence-based sensitivity testing; it remains to be determined whether improvements in DNA extraction protocols alone can deliver this or whether culture, capture or amplification steps will be required. Nonetheless, we can foresee a tipping point when a unified automated metagenomics-based workflow might start to compete with the plethora of methods currently in use in the diagnostic microbiology laboratory.

Diversity and evolution of Mycobacterium tuberculosis: moving to whole-genome-based approaches

Genotyping of clinical Mycobacterium tuberculosis complex (MTBC) strains has become a standard tool for epidemiological tracing and for the investigation of the local and global strain population structure. Of special importance is the analysis of the expansion of multidrug (MDR) and extensively drug-resistant (XDR) strains. Classical genotyping and, more recently, whole-genome sequencing have revealed that the strains of the MTBC are more diverse than previously anticipated. Globally, several phylogenetic lineages can be distinguished whose geographical distribution is markedly variable. Strains of particular (sub)lineages, such as Beijing, seem to be more virulent and associated with enhanced resistance levels and fitness, likely fueling their spread in certain world regions. The upcoming generalization of whole-genome sequencing approaches will expectedly provide more comprehensive insights into the molecular and epidemiological mechanisms involved and lead to better diagnostic and therapeutic tools.

Temporal dynamics of Mycobacterium tuberculosis genotypes in New South Wales, Australia

BACKGROUND

Molecular epidemiology of Mycobacterium tuberculosis, its transmission dynamics and population structure have become important determinants of targeted tuberculosis control programs. Here we describe recent changes in the distribution of M. tuberculosis genotypes in New South Wales (NSW), Australia and compared strain types with drug resistance, site of disease and demographic data.

METHODS

We evaluated all culture-confirmed newly identified tuberculosis cases in NSW, Australia, from 2010-2012. M. tuberculosis population structure and clustering rates were assessed using 24-loci Mycobacterial interspersed repetitive unit (MIRU) analysis and compared to MIRU data from 2006-2008.

RESULTS

Of 1177 tuberculosis cases, 1128 (95.8%) were successfully typed. Beijing and East African Indian (EAI) lineage strains were most common (27.6% and 28.5%, respectively) with EAI strains increasing in relative abundance from 11.8% in 2006-2008 to 28.5% in 2010-2012. Few cases of multi-drug resistant tuberculosis were identified (18; 1.7%). Compared to 12-loci, 24-loci MIRU provided improved cluster resolution with 695 (61.6%) and 227 (20.1%) clustered cases identified, respectively. Detailed analysis of the largest cluster identified (an 11 member Beijing cluster) revealed wide geographic diversity in the absence of documented social contact.

CONCLUSIONS

EAI strains of M. tuberculosis recently overtook Beijing family as a prevalent cause of tuberculosis in NSW, Australia. This lineage appeared to be less commonly related to multi-drug resistant tuberculosis as compared to Beijing strain lineage. The resolution provided by 24-loci MIRU typing was insufficient for reliable assessment of transmissions, especially of Beijing family strains.

Progenitor “Mycobacterium canettii” clone responsible for lymph node tuberculosis epidemic, Djibouti

Two outbreaks among expatriate children were caused by an epidemic clone from the Horn of Africa.

“Mycobacterium canettii,” an opportunistic human pathogen living in an unknown environmental reservoir, is the progenitor species from which Mycobacterium tuberculosis emerged. Since its discovery in 1969, most of the ≈70 known M. canettii strains were isolated in the Republic of Djibouti, frequently from expatriate children and adults. We show here, by whole-genome sequencing, that most strains collected from February 2010 through March 2013, and associated with 2 outbreaks of lymph node tuberculosis in children, belong to a unique epidemic clone within M. canettii. Evolution of this clone, which has been recovered regularly since 1983, may mimic the birth of M. tuberculosis. Thus, recognizing this organism and identifying its reservoir are clinically important.

Whole-genome sequencing to control antimicrobial resistance

Following recent improvements in sequencing technologies, whole-genome sequencing (WGS) is positioned to become an essential tool in the control of antibiotic resistance, a major threat in modern healthcare. WGS has already found numerous applications in this area, ranging from the development of novel antibiotics and diagnostic tests through to antibiotic stewardship of currently available drugs via surveillance and the elucidation of the factors that allow the emergence and persistence of resistance. Numerous proof-of-principle studies have also highlighted the value of WGS as a tool for day-to-day infection control and, for some pathogens, as a primary diagnostic tool to detect antibiotic resistance. However, appropriate data analysis platforms will need to be developed before routine WGS can be introduced on a large scale.

Validation of high throughput sequencing and microbial forensics applications

High throughput sequencing (HTS) generates large amounts of high quality sequence data for microbial genomics. The value of HTS for microbial forensics is the speed at which evidence can be collected and the power to characterize microbial-related evidence to solve biocrimes and bioterrorist events. As HTS technologies continue to improve, they provide increasingly powerful sets of tools to support the entire field of microbial forensics. Accurate, credible results allow analysis and interpretation, significantly influencing the course and/or focus of an investigation, and can impact the response of the government to an attack having individual, political, economic or military consequences. Interpretation of the results of microbial forensic analyses relies on understanding the performance and limitations of HTS methods, including analytical processes, assays and data interpretation. The utility of HTS must be defined carefully within established operating conditions and tolerances. Validation is essential in the development and implementation of microbial forensics methods used for formulating investigative leads attribution. HTS strategies vary, requiring guiding principles for HTS system validation. Three initial aspects of HTS, irrespective of chemistry, instrumentation or software are: 1) sample preparation, 2) sequencing, and 3) data analysis. Criteria that should be considered for HTS validation for microbial forensics are presented here. Validation should be defined in terms of specific application and the criteria described here comprise a foundation for investigators to establish, validate and implement HTS as a tool in microbial forensics, enhancing public safety and national security.

Whole-genome sequencing to detect recent transmission of Mycobacterium tuberculosis in settings with a high burden of tuberculosis

Whole genome sequencing (WGS) of Mycobacterium tuberculosis has been used to trace the transmission of M. tuberculosis, the causative agent of tuberculosis (TB). Previously published studies using WGS were conducted in developed countries with a low TB burden. We sought to evaluate the relative usefulness of traditional VNTR and SNP typing methods, WGS and epidemiological investigations to study the recent transmission of M. tuberculosis in a high TB burden country. We conducted epidemiological investigations of 42 TB patients whose M. tuberculosis isolates were classified into three clusters based on variable-number tandem repeat (VNTR) typing. We applied WGS to 32 (76.2%) of the 42 strains and calculated the pairwise genomic distances between strains within each cluster. Eighteen (56.3%) of the 32 strains had genomic differences ≥100 SNPs with every other strain, suggesting that direct transmission did not likely occurred. Ten strains were grouped into four WGS-based clusters with genomic distances ≤5 SNPs within each cluster, and confirmed epidemiological links were identified in two of these clusters. Our results indicate that WGS provides reliable resolution for tracing the transmission of M. tuberculosis in high TB burden settings. The high resolution of WGS is particularly useful to confirm or exclude the possibility of direct transmission events defined by traditional typing methods.

PolyTB: a genomic variation map for Mycobacterium tuberculosis

Tuberculosis (TB) caused by Mycobacterium tuberculosis (Mtb) is the second major cause of death from an infectious disease worldwide. Recent advances in DNA sequencing are leading to the ability to generate whole genome information in clinical isolates of M. tuberculosis complex (MTBC). The identification of informative genetic variants such as phylogenetic markers and those associated with drug resistance or virulence will help barcode Mtb in the context of epidemiological, diagnostic and clinical studies. Mtb genomic datasets are increasingly available as raw sequences, which are potentially difficult and computer intensive to process, and compare across studies. Here we have processed the raw sequence data (>1500 isolates, eight studies) to compile a catalogue of SNPs (n = 74,039, 63% non-synonymous, 51.1% in more than one isolate, i.e. non-private), small indels (n = 4810) and larger structural variants (n = 800). We have developed the PolyTB web-based tool (http://pathogenseq.lshtm.ac.uk/polytb) to visualise the resulting variation and important meta-data (e.g. in silico inferred strain-types, location) within geographical map and phylogenetic views. This resource will allow researchers to identify polymorphisms within candidate genes of interest, as well as examine the genomic diversity and distribution of strains. PolyTB source code is freely available to researchers wishing to develop similar tools for their pathogen of interest.

Whole-genome-based Mycobacterium tuberculosis surveillance: a standardized, portable, and expandable approach

Whole-genome sequencing (WGS) allows for effective tracing of Mycobacterium tuberculosis complex (MTBC) (tuberculosis pathogens) transmission. However, it is difficult to standardize and, therefore, is not yet employed for interlaboratory prospective surveillance. To allow its widespread application, solutions for data standardization and storage in an easily expandable database are urgently needed. To address this question, we developed a core genome multilocus sequence typing (cgMLST) scheme for clinical MTBC isolates using the Ridom SeqSphere(+) software, which transfers the genome-wide single nucleotide polymorphism (SNP) diversity into an allele numbering system that is standardized, portable, and not computationally intensive. To test its performance, we performed WGS analysis of 26 isolates with identical IS6110 DNA fingerprints and spoligotyping patterns from a longitudinal outbreak in the federal state of Hamburg, Germany (notified between 2001 and 2010). The cgMLST approach (3,041 genes) discriminated the 26 strains with a resolution comparable to that of SNP-based WGS typing (one major cluster of 22 identical or closely related and four outlier isolates with at least 97 distinct SNPs or 63 allelic variants). Resulting tree topologies are highly congruent and grouped the isolates in both cases analogously. Our data show that SNP- and cgMLST-based WGS analyses facilitate high-resolution discrimination of longitudinal MTBC outbreaks. cgMLST allows for a meaningful epidemiological interpretation of the WGS genotyping data. It enables standardized WGS genotyping for epidemiological investigations, e.g., on the regional public health office level, and the creation of web-accessible databases for global TB surveillance with an integrated early warning system.

Genome-wide Mycobacterium tuberculosis variation (GMTV) database: a new tool for integrating sequence variations and epidemiology

Background

Tuberculosis (TB) poses a worldwide threat due to advancing multidrug-resistant strains and deadly co-infections with Human immunodeficiency virus. Today large amounts of Mycobacterium tuberculosis whole genome sequencing data are being assessed broadly and yet there exists no comprehensive online resource that connects M. tuberculosis genome variants with geographic origin, with drug resistance or with clinical outcome.

Description

Here we describe a broadly inclusive unifying Genome-wide Mycobacterium tuberculosis Variation (GMTV) database, (http://mtb.dobzhanskycenter.org) that catalogues genome variations of M. tuberculosis strains collected across Russia. GMTV contains a broad spectrum of data derived from different sources and related to M. tuberculosis molecular biology, epidemiology, TB clinical outcome, year and place of isolation, drug resistance profiles and displays the variants across the genome using a dedicated genome browser. GMTV database, which includes 1084 genomes and over 69,000 SNP or Indel variants, can be queried about M. tuberculosis genome variation and putative associations with drug resistance, geographical origin, and clinical stages and outcomes.

Conclusions

Implementation of GMTV tracks the pattern of changes of M. tuberculosis strains in different geographical areas, facilitates disease gene discoveries associated with drug resistance or different clinical sequelae, and automates comparative genomic analyses among M. tuberculosis strains.

Selected insights from application of whole-genome sequencing for outbreak investigations

PURPOSE OF REVIEW

The advent of high-throughput whole-genome sequencing has the potential to revolutionize the conduct of outbreak investigation. Because of its ultimate resolution power for differentiating between closely related pathogen strains, whole-genome sequencing could augment the traditional epidemiologic investigations of infectious disease outbreaks.

RECENT FINDINGS

The combination of whole-genome sequencing and intensive epidemiologic analysis provided new insights on the sources and transmission dynamics of large-scale epidemics caused by Escherichia coli and Vibrio cholerae, nosocomial outbreaks caused by methicillin-resistant Staphylococcus aureus, Klebsiella pneumoniae, Mycobacterium abscessus, community-centered outbreaks caused by Mycobacterium tuberculosis, and natural disaster-associated outbreaks caused by environmentally acquired molds.

SUMMARY

When combined with traditional epidemiologic investigation, whole-genome sequencing has proven useful for elucidating the sources and transmission dynamics of disease outbreaks. Development of a fully automated bioinformatics pipeline for the analysis of whole-genome sequence data is much needed to make this powerful tool more widely accessible.

Profiling of rpoB mutations and MICs for rifampin and rifabutin in Mycobacterium tuberculosis

Resistance to rifampin (RIF) and rifabutin (RFB) in Mycobacterium tuberculosis is associated with mutations within an 81-bp region of the rpoB gene (RIF resistance-determining region [RRDR]). Previous studies have shown that certain mutations in this region are more likely to confer high levels of RIF resistance, while others may be found in phenotypically susceptible isolates. In this study, we sought to determine the relationship between the MICs of RIF and RFB and rpoB RRDR mutations in 32 multidrug-resistant (MDR), 4 RIF-monoresistant, and 5 susceptible M. tuberculosis clinical isolates. The MICs were determined using the MGIT 960 system. Mutations in the rpoB RRDR were determined by Sanger sequencing. RpoB proteins with mutations S531L (a change of S to L at position 531), S531W, H526Y, and H526D and the double mutation D516A-R529Q were associated with high MICs for RIF and RFB. Five isolates carrying the mutations L511P, H526L, H526N, and D516G-S522L were found to be susceptible to RIF. Several mutations were associated with resistance to RIF and susceptibility to RFB (F514FF, D516V, and S522L). Whole-genome sequencing of two MDR isolates without rpoB RRDR mutations revealed a mutation outside the RRDR (V146F; RIF MIC of 50 μg/ml). The implications of the polymorphisms identified in the second of these isolates in RIF resistance need to be further explored. Our study further establishes a correlation between the mutations and the MICs of RIF and, also, RFB in M. tuberculosis. Several rpoB mutations were identified in RIF- and RFB-susceptible isolates. The clinical significance of these findings requires further exploration. Until then, a combination of phenotypic and molecular testing is advisable for drug susceptibility testing.

High-Throughput Sequencing, a VersatileWeapon to Support Genome-Based Diagnosis in Infectious Diseases: Applications to Clinical Bacteriology

The recent progresses of high-throughput sequencing (HTS) technologies enable easy and cost-reduced access to whole genome sequencing (WGS) or re-sequencing. HTS associated with adapted, automatic and fast bioinformatics solutions for sequencing applications promises an accurate and timely identification and characterization of pathogenic agents. Many studies have demonstrated that data obtained from HTS analysis have allowed genome-based diagnosis, which has been consistent with phenotypic observations. These proofs of concept are probably the first steps toward the future of clinical microbiology. From concept to routine use, many parameters need to be considered to promote HTS as a powerful tool to help physicians and clinicians in microbiological investigations. This review highlights the milestones to be completed toward this purpose.

Supersize me: how whole-genome sequencing and big data are transforming epidemiology

•

Whole-genome sequencing is used for forensic epidemiology.

•

Big data can transform forensic epidemiology.

•

Clustering, biases, wildlife reservoirs, and emerging infections can all be addressed.

•

Phylodynamics approaches to integrate epidemiological and evolutionary data have been highly successful but still face scientific challenges.

In epidemiology, the identification of ‘who infected whom’ allows us to quantify key characteristics such as incubation periods, heterogeneity in transmission rates, duration of infectiousness, and the existence of high-risk groups. Although invaluable, the existence of many plausible infection pathways makes this difficult, and epidemiological contact tracing either uncertain, logistically prohibitive, or both. The recent advent of next-generation sequencing technology allows the identification of traceable differences in the pathogen genome that are transforming our ability to understand high-resolution disease transmission, sometimes even down to the host-to-host scale. We review recent examples of the use of pathogen whole-genome sequencing for the purpose of forensic tracing of transmission pathways, focusing on the particular problems where evolutionary dynamics must be supplemented by epidemiological information on the most likely timing of events as well as possible transmission pathways. We also discuss potential pitfalls in the over-interpretation of these data, and highlight the manner in which a confluence of this technology with sophisticated mathematical and statistical approaches has the potential to produce a paradigm shift in our understanding of infectious disease transmission and control.

Integrating social contact and environmental data in evaluating tuberculosis transmission in a South African township

BACKGROUND

Population models of tuberculosis transmission have not accounted for social contact structure and the role of the environment in which tuberculosis is transmitted.

METHODS

We utilized extensions to the Wells-Riley model of tuberculosis transmission, using exhaled carbon dioxide as a tracer gas, to describe transmission patterns in an endemic community. Drawing upon social interaction data and carbon dioxide measurements from a South African township, we created an age-structured model of tuberculosis transmission in households, public transit, schools, and workplaces. We fit the model to local data on latent tuberculosis prevalence by age.

RESULTS

Most tuberculosis infections (84%) were estimated to occur outside of one's own household. Fifty percent of infections among young adults (ages 15-19) occurred in schools, due to high contact rates and poor ventilation. Despite lower numbers of contacts in workplaces, assortative mixing among adults with high rates of smear-positive tuberculosis contributed to transmission in this environment. Households and public transit were important sites of transmission between age groups.

CONCLUSIONS

Consistent with molecular epidemiologic estimates, a minority of tuberculosis transmission was estimated to occur within households, which may limit the impact of contact investigations. Further work is needed to investigate the role of schools in tuberculosis transmission.

eCAMBer: efficient support for large-scale comparative analysis of multiple bacterial strains

BACKGROUND

Inconsistencies are often observed in the genome annotations of bacterial strains. Moreover, these inconsistencies are often not reflected by sequence discrepancies, but are caused by wrongly annotated gene starts as well as mis-identified gene presence. Thus, tools are needed for improving annotation consistency and accuracy among sets of bacterial strain genomes.

RESULTS

We have developed eCAMBer, a tool for efficiently supporting comparative analysis of multiple bacterial strains within the same species. eCAMBer is a highly optimized revision of our earlier tool, CAMBer, scaling it up for significantly larger datasets comprising hundreds of bacterial strains. eCAMBer works in two phases. First, it transfers gene annotations among all considered bacterial strains. In this phase, it also identifies homologous gene families and annotation inconsistencies. Second, eCAMBer, tries to improve the quality of annotations by resolving the gene start inconsistencies and filtering out gene families arising from annotation errors propagated in the previous phase.

CONCLUSIONS

[corrected] eCAMBer efficiently identifies and resolves annotation inconsistencies among closely related bacterial genomes. It outperforms other competing tools both in terms of running time and accuracy of produced annotations. Software, user manual, and case study results are available at the project website: http://bioputer.mimuw.edu.pl/ecamber.

Assessment of Mycobacterium tuberculosis transmission in Oxfordshire, UK, 2007-12, with whole pathogen genome sequences: an observational study

Background

Patients born outside the UK have contributed to a 20% rise in the UK’s tuberculosis incidence since 2000, but their effect on domestic transmission is not known. Here we use whole-genome sequencing to investigate the epidemiology of tuberculosis transmission in an unselected population over 6 years.

Methods

We identified all residents with Oxfordshire postcodes with a Mycobacterium tuberculosis culture or a clinical diagnosis of tuberculosis between Jan 1, 2007, and Dec 31, 2012, using local databases and checking against the national Enhanced Tuberculosis Surveillance database. We used Illumina technology to sequence all available M tuberculosis cultures from identified cases. Sequences were clustered by genetic relatedness and compared retrospectively with contact investigations. The first patient diagnosed in each cluster was defined as the index case, with links to subsequent cases assigned first by use of any epidemiological linkage, then by genetic distance, and then by timing of diagnosis.

Findings

Although we identified 384 patients with a diagnosis of tuberculosis, country of birth was known for 380 and we sequenced isolates from 247 of 269 cases with culture-confirmed disease. 39 cases were genomically linked within 13 clusters, implying 26 local transmission events. Only 11 of 26 possible transmissions had been previously identified through contact tracing. Of seven genomically confirmed household clusters, five contained additional genomic links to epidemiologically unidentified non-household members. 255 (67%) patients were born in a country with high tuberculosis incidence, conferring a local incidence of 109 cases per 100 000 population per year in Oxfordshire, compared with 3·5 cases per 100 000 per year for those born in low-incidence countries. However, patients born in the low-incidence countries, predominantly UK, were more likely to have pulmonary disease (adjusted odds ratio 1·8 [95% CI 1·2–2·9]; p=0·009), social risk factors (4·4 [2·0–9·4]; p<0·0001), and be part of a local transmission cluster (4·8 [1·6–14·8]; p=0·006).

Interpretation

Although inward migration has contributed to the overall tuberculosis incidence, our findings suggest that most patients born in high-incidence countries reactivate latent infection acquired abroad and are not involved in local onward transmission. Systematic screening of new entrants could further improve tuberculosis control, but it is important that health care remains accessible to all individuals, especially high-risk groups, if tuberculosis control is not to be jeopardised.

Microevolution during an Anthrax outbreak leading to clonal heterogeneity and penicillin resistance

Anthrax is a bacterial disease primarily affecting grazing animals but it can also cause severe disease in humans. We have used genomic epidemiology to study microevolution of the bacterium in a confined outbreak in cattle which involved emergence of an antibiotic-resistant phenotype. At the time of death, the animals contained a heterogeneous population of Single Nucleotide Variants (SNVs), some being clonal but most being subclonal. We found that independent isolates from the same carcass had similar levels of SNV differences as isolates from different animals. Furthermore the relative levels of subclonal populations were different in different locations in the same carcass. The heterogeneity appeared to be derived in part from heterogeneity in the infectious dose. The resistance phenotype was linked to clonal mutations in an anti-sigma factor gene and in one case was preceded by an acquisition of a hypermutator phenotype. In another animal, small subclonal populations were observed with counteracting mutations that had turned off the resistance genes. In summary, this study shows the importance of accounting for both acquired and inherited heterogeneity when doing high-resolution infection tracing and when estimating the risks associated with penicillin treatment.

Whole genome sequencing reveals complex evolution patterns of multidrug-resistant Mycobacterium tuberculosis Beijing strains in patients

Multidrug-resistant (MDR) Mycobacterium tuberculosis complex (MTBC) strains represent a major threat for tuberculosis (TB) control. Treatment of MDR-TB patients is long and less effective, resulting in a significant number of treatment failures. The development of further resistances leads to extensively drug-resistant (XDR) variants. However, data on the individual reasons for treatment failure, e.g. an induced mutational burst, and on the evolution of bacteria in the patient are only sparsely available. To address this question, we investigated the intra-patient evolution of serial MTBC isolates obtained from three MDR-TB patients undergoing longitudinal treatment, finally leading to XDR-TB. Sequential isolates displayed identical IS6110 fingerprint patterns, suggesting the absence of exogenous re-infection. We utilized whole genome sequencing (WGS) to screen for variations in three isolates from Patient A and four isolates from Patient B and C, respectively. Acquired polymorphisms were subsequently validated in up to 15 serial isolates by Sanger sequencing. We determined eight (Patient A) and nine (Patient B) polymorphisms, which occurred in a stepwise manner during the course of the therapy and were linked to resistance or a potential compensatory mechanism. For both patients, our analysis revealed the long-term co-existence of clonal subpopulations that displayed different drug resistance allele combinations. Out of these, the most resistant clone was fixed in the population. In contrast, baseline and follow-up isolates of Patient C were distinguished each by eleven unique polymorphisms, indicating an exogenous re-infection with an XDR strain not detected by IS6110 RFLP typing. Our study demonstrates that intra-patient microevolution of MDR-MTBC strains under longitudinal treatment is more complex than previously anticipated. However, a mutator phenotype was not detected. The presence of different subpopulations might confound phenotypic and molecular drug resistance tests. Furthermore, high resolution WGS analysis is necessary to accurately detect exogenous re-infection as classical genotyping lacks discriminatory power in high incidence settings.

Genomics and metagenomics in medical microbiology

Over the last two decades, sequencing tools have evolved from laborious time-consuming methodologies to real-time detection and deciphering of genomic DNA. Genome sequencing, especially using next generation sequencing (NGS) has revolutionized the landscape of microbiology and infectious disease. This deluge of sequencing data has not only enabled advances in fundamental biology but also helped improve diagnosis, typing of pathogen, virulence and antibiotic resistance detection, and development of new vaccines and culture media. In addition, NGS also enabled efficient analysis of complex human micro-floras, both commensal, and pathological, through metagenomic methods, thus helping the comprehension and management of human diseases such as obesity. This review summarizes technological advances in genomics and metagenomics relevant to the field of medical microbiology.