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

Deciphering the Tangible Spatio-Temporal Spread of a 25-Year Tuberculosis Outbreak Boosted by Social Determinants

Outbreak strains of Mycobacterium tuberculosis are promising candidates as targets in the search for intrinsic determinants of transmissibility, as they are responsible for many cases with sustained transmission; however, the use of low-resolution typing methods and restricted geographical investigations represent flaws in assessing the success of long-lived outbreak strains. We can now address the nature of outbreak strains by combining large genomic data sets and phylodynamic approaches. We retrospectively sequenced the whole genome of representative samples assigned to an outbreak circulating in the Canary Islands (the GC strain) since 1993, which accounts for ~20% of local tuberculosis cases. We selected a panel of specific single nucleotide polymorphism (SNP) markers for an in-silico search for additional outbreak-related sequences within publicly available tuberculosis genomic data. Using this information, we inferred the origin, spread, and epidemiological parameters of the GC strain. Our approach allowed us to accurately trace the historical and more recent dispersion of the GC strain. We provide evidence of a highly successful nature within the Canarian archipelago but limited expansion abroad. Estimation of epidemiological parameters from genomic data disagree with a distinctive biology of the GC strain. With the increasing availability of genomic data allowing for the accurate inference of strain spread and critical epidemiological parameters, we can now revisit the link between Mycobacterium tuberculosis genotypes and transmission, as is routinely carried out for SARS-CoV-2 variants of concern. We demonstrate that social determinants rather than intrinsically higher bacterial transmissibility better explain the success of the GC strain. Importantly, our approach can be used to trace and characterize strains of interest worldwide. IMPORTANCE Infectious disease outbreaks represent a significant problem for public health. Tracing outbreak expansion and understanding the main factors behind emergence and persistence remain critical to effective disease control. Our study allows researchers and public health authorities to use Whole-Genome Sequencing-based methods to trace outbreaks, and shows how available epidemiological information helps to evaluate the factors underpinning outbreak persistence. Taking advantage of all the freely available information placed in public repositories, researchers can accurately establish the expansion of an outbreak beyond original boundaries, and determine the potential risk of a strain to inform health authorities which, in turn, can define target strategies to mitigate expansion and persistence. Finally, we show the need to evaluate strain transmissibility in different geographic contexts to unequivocally associate spread to local or pathogenic factors, an important lesson taken from genomic surveillance of SARS-CoV-2.

Methods Combining Genomic and Epidemiological Data in the Reconstruction of Transmission Trees: A Systematic Review

In order to better understand transmission dynamics and appropriately target control and preventive measures, studies have aimed to identify who-infected-whom in actual outbreaks. Numerous reconstruction methods exist, each with their own assumptions, types of data, and inference strategy. Thus, selecting a method can be difficult. Following PRISMA guidelines, we systematically reviewed the literature for methods combing epidemiological and genomic data in transmission tree reconstruction. We identified 22 methods from the 41 selected articles. We defined three families according to how genomic data was handled: a non-phylogenetic family, a sequential phylogenetic family, and a simultaneous phylogenetic family. We discussed methods according to the data needed as well as the underlying sequence mutation, within-host evolution, transmission, and case observation. In the non-phylogenetic family consisting of eight methods, pairwise genetic distances were estimated. In the phylogenetic families, transmission trees were inferred from phylogenetic trees either simultaneously (nine methods) or sequentially (five methods). While a majority of methods (17/22) modeled the transmission process, few (8/22) took into account imperfect case detection. Within-host evolution was generally (7/8) modeled as a coalescent process. These practical and theoretical considerations were highlighted in order to help select the appropriate method for an outbreak.

Reporting practices for genomic epidemiology of tuberculosis: a systematic review of the literature using STROME-ID guidelines as a benchmark

BACKGROUND

Pathogen genomics have become increasingly important in infectious disease epidemiology and public health. The Strengthening the Reporting of Molecular Epidemiology for Infectious Diseases (STROME-ID) guidelines were developed to outline a minimum set of criteria that should be reported in genomic epidemiology studies to facilitate assessment of study quality. We evaluate such reporting practices, using tuberculosis as an example.

METHODS

For this systematic review, we initially searched MEDLINE, Embase Classic, and Embase on May 3, 2017, using the search terms "tuberculosis" and "genom* sequencing". We updated this initial search on April 23, 2019, and also included a search of bioRxiv at this time. We included studies in English, French, or Spanish that recruited patients with microbiologically confirmed tuberculosis and used whole genome sequencing for typing of strains. Non-human studies, conference abstracts, and literature reviews were excluded. For each included study, the number and proportion of fulfilled STROME-ID criteria were recorded by two reviewers. A comparison of the mean proportion of fulfilled STROME-ID criteria before and after publication of the STROME-ID guidelines (in 2014) was done using a two-tailed t test. Quasi-Poisson regression and tobit regression were used to examine associations between study characteristics and the number and proportion of fulfilled STROME-ID criteria. This study was registered with PROSPERO, CRD42017064395.

FINDINGS

976 titles and abstracts were identified by our primary search, with an additional 16 studies identified in bioRxiv. 114 full texts (published between 2009 and 2019) were eligible for inclusion. The mean proportion of STROME-ID criteria fulfilled was 50% (SD 12; range 16-75). The proportion of criteria fulfilled was similar before and after STROME-ID publication (51% [SD 11] vs 46% [14], p=0·26). The number of criteria reported (among those applicable to all studies) was not associated with impact factor, h-index, country of affiliation of senior author, or sample size of isolates. Similarly, the proportion of criteria fulfilled was not associated with these characteristics, with the exception of a sample size of isolates of 277 or more (the highest quartile). In terms of reproducibility, 100 (88%) studies reported which bioinformatic tools were used, but only 33 (33%) reported corresponding version numbers. Sequencing data were available for 86 (75%) studies.

INTERPRETATION

The reporting of STROME-ID criteria in genomic epidemiology studies of tuberculosis between 2009 and 2019 was low, with implications for assessment of study quality. The considerable proportion of studies without bioinformatics version numbers or sequencing data available highlights a key concern for reproducibility.

Computational pan-genome mapping and pairwise SNP-distance improve detection of Mycobacterium tuberculosis transmission clusters

Next-generation sequencing based base-by-base distance measures have become an integral complement to epidemiological investigation of infectious disease outbreaks. This study introduces PANPASCO, a computational pan-genome mapping based, pairwise distance method that is highly sensitive to differences between cases, even when located in regions of lineage specific reference genomes. We show that our approach is superior to previously published methods in several datasets and across different Mycobacterium tuberculosis lineages, as its characteristics allow the comparison of a high number of diverse samples in one analysis—a scenario that becomes more and more likely with the increased usage of whole-genome sequencing in transmission surveillance.

Tuberculosis still is a threat to global health. It is essential to detect and interrupt transmissions to stop the spread of this infectious disease. With the rising use of next-generation sequencing methods, its application in the surveillance of Mycobacterium tuberculosis has become increasingly important in the last years. The main goal of molecular surveillance is the identification of patient-patient transmission and cluster detection. The mutation rate of M. tuberculosis is very low and stable. Therefore, many existing methods for comparative analysis of isolates provide inadequate results since their resolution is too limited. There is a need for a method that takes every detectable difference into account. We developed PANPASCO, a novel approach for comparing pairs of isolates using all genomic information available for each pair. We combine improved SNP-distance calculation with the use of a pan-genome incorporating more than 100 M. tuberculosis reference genomes representing lineages 1-4 for read mapping prior to variant detection. We thereby enable the collective analysis and comparison of similar and diverse isolates associated with different M. tuberculosis strains.

Molecular Epidemiology of Tuberculosis in British Columbia, Canada: A 10-Year Retrospective Study

Background

Understanding regional molecular epidemiology allows for the development of more efficient tuberculosis prevention strategies in low-incidence settings.

Methods

We analyzed 24-locus mycobacterial interspersed repetitive-unit–variable-number tandem repeat (MIRU-VNTR) genotyping for 2290 Mycobacterium tuberculosis clinical isolates collected in the province of British Columbia (BC), Canada, in 2005–2014. Laboratory data for each isolate were linked to case-level clinical and demographic data. These data were used to describe the molecular epidemiology of tuberculosis across the province.

Results

We detected >1500 distinct genotypes across the 4 major M. tuberculosis lineages, reflecting BC’s diverse population. Disease site and clustering rates varied across lineages, and MIRU-VNTR was used to group the 2290 isolates into 189 clusters (2–70 isolates per cluster), with an overall clustering rate of 42.4% and an estimated local transmission rate of 34.1%. Risk factors for clustering varied between Canadian-born and foreign-born individuals; the former had increased odds (odds ratio, 7.8; 95% confidence interval [CI], 6.2–9.6) of belonging to a genotypic cluster, although nearly one-quarter of clusters included both Canadian- and foreign-born persons. Large clusters (≥10 cases) occurred more frequently within the M. tuberculosis Euro-American lineage, and individual-level risk factors associated with belonging to a large cluster included being Canadian born (adjusted odds ratio, 3.3; 95% CI, 2.3–4.8), residing in a rural area (2.3; 1.2–4.5), and illicit drug use (2.0; 1.2–3.4).

Conclusions

Although tuberculosis in BC largely arises through reactivation of latent tuberculosis in foreign-born persons, locally transmitted infections occur in discrete populations with distinct disease and risk factor profiles, representing groups for targeted interventions.

Mycobacterium avium in Community and Household Water, Suburban Philadelphia, Pennsylvania, USA, 2010-2012

Attention to environmental sources of Mycobacterium avium complex (MAC) infection is a vital component of disease prevention and control. We investigated MAC colonization of household plumbing in suburban Philadelphia, Pennsylvania, USA. We used variable-number tandem-repeat genotyping and whole-genome sequencing with core genome single-nucleotide variant analysis to compare M. avium from household plumbing biofilms with M. avium isolates from patient respiratory specimens. M. avium was recovered from 30 (81.1%) of 37 households, including 19 (90.5%) of 21 M. avium patient households. For 11 (52.4%) of 21 patients with M. avium disease, isolates recovered from their respiratory and household samples were of the same genotype. Within the same community, 18 (85.7%) of 21 M. avium respiratory isolates genotypically matched household plumbing isolates. Six predominant genotypes were recovered across multiple households and respiratory specimens. M. avium colonizing municipal water and household plumbing may be a substantial source of MAC pulmonary infection.

Role and value of whole genome sequencing in studying tuberculosis transmission

BACKGROUND

Tuberculosis (TB) remains a serious public health threat worldwide. Theoretically ultimate resolution of whole genome sequencing (WGS) for Mycobacterium tuberculosis complex (MTBC) strain classification makes this technology very attractive for epidemiological investigations.

OBJECTIVES

To summarize the evidence available in peer-reviewed publications on the role and place of WGS in detection of TB transmission.

SOURCES

A total of 69 peer-reviewed publications identified in Pubmed database.

CONTENT

Evidence from >30 publications suggests that a cut-off value of fewer than six single nucleotide polymorphisms between strains efficiently excludes cases that are not the result of recent transmission and could be used for the identification of drug-sensitive isolates involved in direct human-to-human TB transmission. Sensitivity of WGS to identify epidemiologically linked isolates is high, reaching 100% in eight studies with specificity (17%-95%) highly dependent on the settings. Drug resistance and specific phylogenetic lineages may be associated with accelerated mutation rates affecting genetic distances. WGS can be potentially used to distinguish between true relapses and re-infections but in high-incidence low-diversity settings this would require consideration of epidemiological links and minority alleles. Data from four studies looking into within-host diversity highlight a need for developing criteria for acceptance or rejection of WGS relatedness results depending on the proportion of minority alleles.

IMPLICATIONS

WGS will potentially allow for more targeted public health actions preventing unnecessary investigations of false clusters. Consensus on standardization of raw data quality control processing criteria, analytical pipelines and reporting language is yet to be reached.

Universal genotyping reveals province-level differences in the molecular epidemiology of tuberculosis

OBJECTIVES

Compare the molecular epidemiology of tuberculosis (TB) between two large Canadian provinces-Ontario and British Columbia (BC)-to identify genotypic clusters within and across both provinces, allowing for an improved understanding of genotype data and providing context to more accurately identify clusters representing local transmission.

DESIGN

We compared 24-locus Mycobacterial Interspersed Repetitive Units-Variable Number of Tandem Repeats (MIRU-VNTR) genotyping for 3,314 Ontario and 1,602 BC clinical Mycobacterium tuberculosis isolates collected from 2008 through 2014. Laboratory data for each isolate was linked to case-level records to obtain clinical and demographic data.

RESULTS

The demographic characteristics of persons with TB varied between provinces, most notably in the proportion of persons born outside Canada, which was reflected in the large number of unique genotypes (n = 3,461). The proportion of clustered isolates was significantly higher in BC. Substantial clustering amongst non-Lineage 4 TB strains was observed within and across the provinces. Only two large clusters (≥10 cases/cluster) representing within province transmission had interprovincial genotype matches.

CONCLUSION

We recommend expanding analysis of shared genotypes to include neighbouring jurisdictions, and implementing whole genome sequencing to improve identification of TB transmission, recognize outbreaks, and monitor changing trends in TB epidemiology.

Whole genome sequencing for improved understanding of Mycobacterium tuberculosis transmission in a remote circumpolar region

Few studies have used genomic epidemiology to understand tuberculosis (TB) transmission in rural and remote settings - regions often unique in history, geography and demographics. To improve our understanding of TB transmission dynamics in Yukon Territory (YT), a circumpolar Canadian territory, we conducted a retrospective analysis in which we combined epidemiological data collected through routine contact investigations with clinical and laboratory results. Mycobacterium tuberculosis isolates from all culture-confirmed TB cases in YT (2005-2014) were genotyped using 24-locus Mycobacterial Interspersed Repetitive Units-Variable Number of Tandem Repeats (MIRU-VNTR) and compared to each other and to those from the neighbouring province of British Columbia (BC). Whole genome sequencing (WGS) of genotypically clustered isolates revealed three sustained transmission networks within YT, two of which also involved BC isolates. While each network had distinct characteristics, all had at least one individual acting as the probable source of three or more culture-positive cases. Overall, WGS revealed that TB transmission dynamics in YT are distinct from patterns of spread in other, more remote Northern Canadian regions, and that the combination of WGS and epidemiological data can provide actionable information to local public health teams.

Highlights on the Application of Genomics and Bioinformatics in the Fight Against Infectious Diseases: Challenges and Opportunities in Africa

Genomics and bioinformatics are increasingly contributing to our understanding of infectious diseases caused by bacterial pathogens such as Mycobacterium tuberculosis and parasites such as Plasmodium falciparum. This ranges from investigations of disease outbreaks and pathogenesis, host and pathogen genomic variation, and host immune evasion mechanisms to identification of potential diagnostic markers and vaccine targets. High throughput genomics data generated from pathogens and animal models can be combined with host genomics and patients’ health records to give advice on treatment options as well as potential drug and vaccine interactions. However, despite accounting for the highest burden of infectious diseases, Africa has the lowest research output on infectious disease genomics. Here we review the contributions of genomics and bioinformatics to the management of infectious diseases of serious public health concern in Africa including tuberculosis (TB), dengue fever, malaria and filariasis. Furthermore, we discuss how genomics and bioinformatics can be applied to identify drug and vaccine targets. We conclude by identifying challenges to genomics research in Africa and highlighting how these can be overcome where possible.

Technology and tuberculosis control: the OUT-TB Web experience

Objective

Develop a tool to disseminate integrated laboratory, clinical, and demographic case data necessary for improved contact tracing and outbreak detection of tuberculosis (TB).

Methods

In 2007, the Public Health Ontario Laboratories implemented a universal genotyping program to monitor the spread of TB strains within Ontario. Ontario Universal Typing of TB (OUT-TB) Web utilizes geographic information system (GIS) technology with a relational database platform, allowing TB control staff to visualize genotyping matches and microbiological data within the context of relevant epidemiological and demographic data.

Results

OUT-TB Web is currently available to the 8 health units responsible for >85% of Ontario's TB cases and is a valuable tool for TB case investigation. Users identified key features to implement for application enhancements, including an e-mail alert function, customizable heat maps for visualizing TB and drug-resistant cases, socioeconomic map layers, a dashboard providing TB surveillance metrics, and a feature for animating the geographic spread of strains over time.

Conclusion

OUT-TB Web has proven to be an award-winning application and a useful tool. Developed and enhanced using regular user feedback, future versions will include additional data sources, enhanced map and line-list filter capabilities, and development of a mobile app.

Bacterial genome sequencing in clinical microbiology: a pathogen-oriented review

In recent years, whole-genome sequencing (WGS) has been perceived as a technology with the potential to revolutionise clinical microbiology. Herein, we reviewed the literature on the use of WGS for the most commonly encountered pathogens in clinical microbiology laboratories: Escherichia coli and other Enterobacteriaceae, Staphylococcus aureus and coagulase-negative staphylococci, streptococci and enterococci, mycobacteria and Chlamydia trachomatis. For each pathogen group, we focused on five different aspects: the genome characteristics, the most common genomic approaches and the clinical uses of WGS for (i) typing and outbreak analysis, (ii) virulence investigation and (iii) in silico antimicrobial susceptibility testing. Of all the clinical usages, the most frequent and straightforward usage was to type bacteria and to trace outbreaks back. A next step toward standardisation was made thanks to the development of several new genome-wide multi-locus sequence typing systems based on WGS data. Although virulence characterisation could help in various particular clinical settings, it was done mainly to describe outbreak strains. An increasing number of studies compared genotypic to phenotypic antibiotic susceptibility testing, with mostly promising results. However, routine implementation will preferentially be done in the workflow of particular pathogens, such as mycobacteria, rather than as a broadly applicable generic tool. Overall, concrete uses of WGS in routine clinical microbiology or infection control laboratories were done, but the next big challenges will be the standardisation and validation of the procedures and bioinformatics pipelines in order to reach clinical standards.

Application of whole genome sequence analysis to the study of Mycobacterium tuberculosis in Nunavut, Canada

Canada has one of the lowest rates of tuberculosis (TB) in the world, however, among certain sub-populations, disease incidence rates approach those observed in sub-Saharan Africa, and other high incidence regions. In this study, we applied mycobacterial interspersed repetitive unit (MIRU) variable number of tandem repeat (VNTR) and whole genome sequencing (WGS) to the analysis of Mycobacterium tuberculosis isolates obtained from Northern communities in the territory of Nunavut. WGS was carried out using the Illumina MiSeq, with identified variants used to infer phylogenetic relationships and annotated to infer functional implications. Additionally, the sequencing data from these isolates were augmented with publically available WGS to evaluate data from the Nunavut outbreak in the broader Canadian context. In this study, isolates could be classified into four major clusters by MIRU-VNTR analysis. These could be further resolved into sub-clusters using WGS. No evidence for antimicrobial resistance, either genetic or phenotypic, was observed in this cohort. Among most subjects with multiple samples, reactivation/incomplete treatment likely contributed to recurrence. However, isolates from two subjects appeared more likely to have occurred via reinfection, based on the large number of genomic single nucleotide variants detected. Finally, although quite distinct from previously reported Canadian MTB strains, isolates obtained from Nunavut clustered most closely with a cohort of samples originating in the Nunavik region of Northern Quebec. This study demonstrates the benefit of using WGS for discriminatory analysis of MTB in Canada, especially in high incidence regions. It further emphasizes the importance of focusing epidemiological intervention efforts on interrupting transmission chains of endemic TB throughout Northern communities, rather than relying on strategies applied in regions where the majority of TB cases result from importation of foreign strains.

Genomic Epidemiology of a Major Mycobacterium tuberculosis Outbreak: Retrospective Cohort Study in a Low-Incidence Setting Using Sparse Time-Series Sampling

Since 1992, Denmark has documented the largest outbreak of tuberculosis in Scandinavia ascribed to a single genotype, termed C2/1112-15. As of spring 2017, the International Reference Laboratory of Mycobacteriology in Copenhagen has collected and identified isolates from more than a thousand cases belonging to this outbreak via routine mycobacterial interspersed repetitive units-variable number of tandem repeats typing. Here, we present a retrospective analysis of the C2/1112-15 dataset, based on whole-genome data from a sparse time series consisting of 5 randomly selected isolates from 23 years of sampling. Even if these data are derived from only 12% of the collected isolates, we have been able to extract important key information, such as mutation rate and conserved single-nucleotide polymorphisms to identify discrete transmission chains, as well as the possible historical origins of the outbreak.

Whole transcriptomic and proteomic analyses of an isogenic M. tuberculosis clinical strain with a naturally occurring 15 Kb genomic deletion

Tuberculosis remains one of the most difficult to control infectious diseases in the world. Many different factors contribute to the complexity of this disease. These include the ability of the host to control the infection which may directly relate to nutritional status, presence of co-morbidities and genetic predisposition. Pathogen factors, in particular the ability of different Mycobacterium tuberculosis strains to respond to the harsh environment of the host granuloma, which includes low oxygen and nutrient availability and the presence of damaging radical oxygen and nitrogen species, also play an important role in the success of different strains to cause disease. In this study we evaluated the impact of a naturally occurring 12 gene 15 Kb genomic deletion on the physiology and virulence of M. tuberculosis. The strains denominated ON-A WT (wild type) and ON-A NM (natural mutant) were isolated from a previously reported TB outbreak in an inner city under-housed population in Toronto, Canada. Here we subjected these isogenic strains to transcriptomic (via RNA-seq) and proteomic analyses and identified several gene clusters with differential expression in the natural mutant, including the DosR regulon and the molybdenum cofactor biosynthesis genes, both of which were found in lower abundance in the natural mutant. We also demonstrated lesser virulence of the natural mutant in the guinea pig animal model. Overall, our findings suggest that the ON-A natural mutant is less fit to cause disease, but nevertheless has the potential to cause extended transmission in at-risk populations.

Outbreak of tuberculosis among substance users and homeless people in Greater Montréal, Canada, 2003-2016

Background

In Canada, active tuberculosis (TB) is found mainly among migrants from endemic countries and Indigenous populations. However, cases of active tuberculosis in substance users and homeless persons have been reported in Greater Montréal since 2003.

Objective

To describe the Montréal TB outbreak in terms of the sociodemographic characteristics, risk factors and clinical characteristics of cases, as well as the intensity of public health interventions, the follow-up and identification of locations of potential transmission.

Methods

All cases of active tuberculosis with the same genotype of interest residing in Quebec and epidemiologically linked cases were included in the analysis. Data were retrospectively extracted from routine public health investigations. Characteristics of cases were summarized using Excel. Spatial analysis of locations frequented during cases' infectiousness periods was performed.

Results

Between January 2003 and February 2016 a total of 35 cases were identified. Most (86%) were non-Indigenous people born in Canada. Of these, 28 had several risk factors, including substance use (93%), alcohol abuse (64%), homelessness (46%), comorbidities such as HIV coinfection (36%) and advanced stage of the disease. Seven cases without risk factors were all close contacts of cases. Intensity of case management by public health authorities was high. Locations frequented by cases with risk factors included crack houses, shelters and rehabilitation centers in Montréal's downtown core and a residential setting in a suburban area.

Conclusion

TB outbreaks can occur in marginalized Canadian-born urban populations, especially those with substance use. Tailored interventions in this population may be needed for screening, and earlier identification of both latent and active TB and better linkage to care.

A brief primer on genomic epidemiology: lessons learned from Mycobacterium tuberculosis

Genomics is now firmly established as a technique for the investigation and reconstruction of communicable disease outbreaks, with many genomic epidemiology studies focusing on revealing transmission routes of Mycobacterium tuberculosis. In this primer, we introduce the basic techniques underlying transmission inference from genomic data, using illustrative examples from M. tuberculosis and other pathogens routinely sequenced by public health agencies. We describe the laboratory and epidemiological scenarios under which genomics may or may not be used, provide an introduction to sequencing technologies and bioinformatics approaches to identifying transmission-informative variation and resistance-associated mutations, and discuss how variation must be considered in the light of available clinical and epidemiological information to infer transmission.

Whole Genome Sequence Analysis of a Large Isoniazid-Resistant Tuberculosis Outbreak in London: A Retrospective Observational Study

BACKGROUND

A large isoniazid-resistant tuberculosis outbreak centred on London, United Kingdom, has been ongoing since 1995. The aim of this study was to investigate the power and value of whole genome sequencing (WGS) to resolve the transmission network compared to current molecular strain typing approaches, including analysis of intra-host diversity within a specimen, across body sites, and over time, with identification of genetic factors underlying the epidemiological success of this cluster.

METHODS AND FINDINGS

We sequenced 344 outbreak isolates from individual patients collected over 14 y (2 February 1998-22 June 2012). This demonstrated that 96 (27.9%) were indistinguishable, and only one differed from this major clone by more than five single nucleotide polymorphisms (SNPs). The maximum number of SNPs between any pair of isolates was nine SNPs, and the modal distance between isolates was two SNPs. WGS was able to reveal the direction of transmission of tuberculosis in 16 cases within the outbreak (4.7%), including within a multidrug-resistant cluster that carried a rare rpoB mutation associated with rifampicin resistance. Eleven longitudinal pairs of patient pulmonary isolates collected up to 48 mo apart differed from each other by between zero and four SNPs. Extrapulmonary dissemination resulted in acquisition of a SNP in two of five cases. WGS analysis of 27 individual colonies cultured from a single patient specimen revealed ten loci differed amongst them, with a maximum distance between any pair of six SNPs. A limitation of this study, as in previous studies, is that indels and SNPs in repetitive regions were not assessed due to the difficulty in reliably determining this variation.

CONCLUSIONS

Our study suggests that (1) certain paradigms need to be revised, such as the 12 SNP distance as the gold standard upper threshold to identify plausible transmissions; (2) WGS technology is helpful to rule out the possibility of direct transmission when isolates are separated by a substantial number of SNPs; (3) the concept of a transmission chain or network may not be useful in institutional or household settings; (4) the practice of isolating single colonies prior to sequencing is likely to lead to an overestimation of the number of SNPs between cases resulting from direct transmission; and (5) despite appreciable genomic diversity within a host, transmission of tuberculosis rarely results in minority variants becoming dominant. Thus, whilst WGS provided some increased resolution over variable number tandem repeat (VNTR)-based clustering, it was insufficient for inferring transmission in the majority of cases.

Tracing Mycobacterium tuberculosis transmission by whole genome sequencing in a high incidence setting: a retrospective population-based study in East Greenland

In East Greenland, a dramatic increase of tuberculosis (TB) incidence has been observed in recent years. Classical genotyping suggests a genetically similar Mycobacterium tuberculosis (Mtb) strain population as cause, however, precise transmission patterns are unclear. We performed whole genome sequencing (WGS) of Mtb isolates from 98% of culture-positive TB cases through 21 years (n = 182) which revealed four genomic clusters of the Euro-American lineage (mainly sub-lineage 4.8 (n = 134)). The time to the most recent common ancestor of lineage 4.8 strains was found to be 100 years. This sub-lineage further diversified in the 1970s, and massively expanded in the 1990s, a period of lowered TB awareness in Greenland. Despite the low genetic strain diversity, WGS data revealed several recent short-term transmission events in line with the increasing incidence in the region. Thus, the isolated setting and the uniformity of circulating Mtb strains indicated that the majority of East Greenlandic TB cases originated from one or few strains introduced within the last century. Thereby, the study shows the consequences of even short interruptions in TB control efforts in previously TB high incidence areas and demonstrates the potential role of WGS in detecting ongoing micro epidemics, thus guiding public health efforts in the future.

Interpreting whole genome sequencing for investigating tuberculosis transmission: a systematic review

BACKGROUND

Whole genome sequencing (WGS) is becoming an important part of epidemiological investigations of infectious diseases due to greater resolution and cost reductions compared to traditional typing approaches. Many public health and clinical teams will increasingly use WGS to investigate clusters of potential pathogen transmission, making it crucial to understand the benefits and assumptions of the analytical methods for investigating the data. We aimed to understand how different approaches affect inferences of transmission dynamics and outline limitations of the methods.

METHODS

We comprehensively searched electronic databases for studies that presented methods used to interpret WGS data for investigating tuberculosis (TB) transmission. Two authors independently selected studies for inclusion and extracted data. Due to considerable methodological heterogeneity between studies, we present summary data with accompanying narrative synthesis rather than pooled analyses.

RESULTS

Twenty-five studies met our inclusion criteria. Despite the range of interpretation tools, the usefulness of WGS data in understanding TB transmission often depends on the amount of genetic diversity in the setting. Where diversity is small, distinguishing re-infections from relapses may be impossible; interpretation may be aided by the use of epidemiological data, examining minor variants and deep sequencing. Conversely, when within-host diversity is large, due to genetic hitchhiking or co-infection of two dissimilar strains, it is critical to understand how it arose. Greater understanding of microevolution and mixed infection will enhance interpretation of WGS data.

CONCLUSIONS

As sequencing studies have sampled more intensely and integrated multiple sources of information, the understanding of TB transmission and diversity has grown, but there is still much to be learnt about the origins of diversity that will affect inferences from these data. Public health teams and researchers should combine epidemiological, clinical and WGS data to strengthen investigations of transmission.

Recombination in pe/ppe genes contributes to genetic variation in Mycobacterium tuberculosis lineages

Background

Approximately 10 % of the Mycobacterium tuberculosis genome is made up of two families of genes that are poorly characterized due to their high GC content and highly repetitive nature. The PE and PPE families are typified by their highly conserved N-terminal domains that incorporate proline-glutamate (PE) and proline-proline-glutamate (PPE) signature motifs. They are hypothesised to be important virulence factors involved with host-pathogen interactions, but their high genetic variability and complexity of analysis means they are typically disregarded in genome studies.

Results

To elucidate the structure of these genes, 518 genomes from a diverse international collection of clinical isolates were de novo assembled. A further 21 reference M. tuberculosis complex genomes and long read sequence data were used to validate the approach. SNP analysis revealed that variation in the majority of the 168 pe/ppe genes studied was consistent with lineage. Several recombination hotspots were identified, notably pe_pgrs3 and pe_pgrs17. Evidence of positive selection was revealed in 65 pe/ppe genes, including epitopes potentially binding to major histocompatibility complex molecules.

Conclusions

This, the first comprehensive study of the pe and ppe genes, provides important insight into M. tuberculosis diversity and has significant implications for vaccine development.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2467-y) contains supplementary material, which is available to authorized users.

Preliminary evaluation of exome sequencing to identify genetic markers of susceptibility to tuberculosis disease

Background

Recent studies have shown that certain human genetic polymorphisms could be associated with susceptibility to tuberculosis (TB) infection and disease. Advances in next generation sequencing include the ability to rapidly sequence the entire human exome. These new technologies can be exploited to identify new associations of human genetic polymorphisms and TB infection and disease. In this preliminary study we compared two different strategies for sequencing of the human exome in a small sample set consisting of three individuals with a history of TB disease and two individuals with latent TB infection.

Findings

Sequencing of the entire exome of the five participants using Agilent SureSelect kit resulted in the identification of 1611 single nucleotide polymorphisms (SNPs) that were only present in the individuals with a history of active TB but not in the latent TB cases. Alternatively, sequencing of 4000 target genes available in the TruSight kit resulted in identification of 182 SNPs only present in the active TB cases and not in the latent TB participants. The overlap of the two kits was 112 SNPs.

Conclusions

Even though this pilot study was restricted to a small number of participants, we demonstrated the feasibility of using exome sequencing technologies to mine potential genetic associations of susceptibility to TB disease and presented a number of potential targets that can be further explore in larger research trials.

Electronic supplementary material

The online version of this article (doi:10.1186/s13104-015-1740-5) contains supplementary material, which is available to authorized users.

Whole genome sequencing reveals genomic heterogeneity and antibiotic purification in Mycobacterium tuberculosis isolates

BACKGROUND

Whole genome sequencing has revolutionised the interrogation of mycobacterial genomes. Recent studies have reported conflicting findings on the genomic stability of Mycobacterium tuberculosis during the evolution of drug resistance. In an age where whole genome sequencing is increasingly relied upon for defining the structure of bacterial genomes, it is important to investigate the reliability of next generation sequencing to identify clonal variants present in a minor percentage of the population. This study aimed to define a reliable cut-off for identification of low frequency sequence variants and to subsequently investigate genetic heterogeneity and the evolution of drug resistance in M. tuberculosis.

METHODS

Genomic DNA was isolated from single colonies from 14 rifampicin mono-resistant M. tuberculosis isolates, as well as the primary cultures and follow up MDR cultures from two of these patients. The whole genomes of the M. tuberculosis isolates were sequenced using either the Illumina MiSeq or Illumina HiSeq platforms. Sequences were analysed with an in-house pipeline.

RESULTS

Using next-generation sequencing in combination with Sanger sequencing and statistical analysis we defined a read frequency cut-off of 30% to identify low frequency M. tuberculosis variants with high confidence. Using this cut-off we demonstrated a high rate of genetic diversity between single colonies isolated from one population, showing that by using the current sequencing technology, single colonies are not a true reflection of the genetic diversity within a whole population and vice versa. We further showed that numerous heterogeneous variants emerge and then disappear during the evolution of isoniazid resistance within individual patients. Our findings allowed us to formulate a model for the selective bottleneck which occurs during the course of infection, acting as a genomic purification event.

CONCLUSIONS

Our study demonstrated true levels of genetic diversity within an M. tuberculosis population and showed that genetic diversity may be re-defined when a selective pressure, such as drug exposure, is imposed on M. tuberculosis populations during the course of infection. This suggests that the genome of M. tuberculosis is more dynamic than previously thought, suggesting preparedness to respond to a changing environment.

Clinical value of whole-genome sequencing of Mycobacterium tuberculosis

Whole-genome sequencing (WGS) is now common as a result of new technologies that can rapidly sequence a complete bacterial genome for US$500 or less. Many studies have addressed questions about tuberculosis with WGS, and knowing the sequence of the entire genome, rather than only a few fragments, has greatly increased the precision of molecular epidemiology and contact tracing. Additionally, topics such as the mutation rate, drug resistance, the target of new drugs, and the phylogeny and evolution of the Mycobacterium tuberculosis complex bacteria have been elucidated by WGS. Nonetheless, WGS has not explained differences in transmissibility between strains, or why some strains are more virulent than others or more prone to development of multidrug resistance. With advances in technology, WGS of clinical specimens could become routine in high-income countries; however, its relevance will probably depend on easy to use software to efficiently process the sequences produced and accessible genomic databases that can be mined in future studies.

A microbiological revolution meets an ancient disease: improving the management of tuberculosis with genomics

Tuberculosis (TB) is an ancient disease with an enormous global impact. Despite declining global incidence, the diagnosis, phenotyping, and epidemiological investigation of TB require significant clinical microbiology laboratory resources. Current methods for the detection and characterization of Mycobacterium tuberculosis consist of a series of laboratory tests varying in speed and performance, each of which yields incremental information about the disease. Since the sequencing of the first M. tuberculosis genome in 1998, genomic tools have aided in the diagnosis, treatment, and control of TB. Here we summarize genomics-based methods that are positioned to be introduced in the modern clinical TB laboratory, and we highlight how recent advances in genomics will improve the detection of antibiotic resistance-conferring mutations and the understanding of M. tuberculosis transmission dynamics and epidemiology. We imagine the future TB clinic as one that relies heavily on genomic interrogation of the M. tuberculosis isolate, allowing for more rapid diagnosis of TB and real-time monitoring of outbreak emergence.

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.