Evolutionary history and global spread of the Mycobacterium tuberculosis Beijing lineage

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.

MTBseq: a comprehensive pipeline for whole genome sequence analysis of Mycobacterium tuberculosis complex isolates

Analyzing whole-genome sequencing data of Mycobacterium tuberculosis complex (MTBC) isolates in a standardized workflow enables both comprehensive antibiotic resistance profiling and outbreak surveillance with highest resolution up to the identification of recent transmission chains. Here, we present MTBseq, a bioinformatics pipeline for next-generation genome sequence data analysis of MTBC isolates. Employing a reference mapping based workflow, MTBseq reports detected variant positions annotated with known association to antibiotic resistance and performs a lineage classification based on phylogenetic single nucleotide polymorphisms (SNPs). When comparing multiple datasets, MTBseq provides a joint list of variants and a FASTA alignment of SNP positions for use in phylogenomic analysis, and identifies groups of related isolates. The pipeline is customizable, expandable and can be used on a desktop computer or laptop without any internet connection, ensuring mobile usage and data security. MTBseq and accompanying documentation is available from https://github.com/ngs-fzb/MTBseq_source.

Genetic diversity and distribution dynamics of multidrug-resistant Mycobacterium tuberculosis isolates in Nepal

Multidrug-resistant tuberculosis (MDR-TB) is an emerging public health problem in Nepal. Despite the implementation of a successful TB control program in Nepal, notifications of MDR-TB are increasing, yet the reasons are unknown. The objective of this study was to understand the genetic diversity and epidemiological characteristics of MDR-Mycobacterium tuberculosis (MTB) isolates in Nepal. We isolated and genotyped 498 MDR-MTB isolates collected from April 2009 to March 2013 and analyzed the patients' background information. Our results showed that the lineage 2 (Beijing family) was the most predominant lineage (n = 241; 48.4%), followed by lineage 3 (n = 153, 30.7%). Lineage 4 was the third most prevalent (n = 73, 14.5%) followed by lineage 1 (n = 32, 6.4%). The lineages were significantly associated with geographic region, ethnic group, age and sex of patients. The Beijing genotype was found to have an important role in transmitting MDR-TB in Nepal and was significantly associated with the eastern region, mongoloid ethnic group and younger age group. We conclude that early diagnosis and treatment including molecular-epidemiological surveillance of MDR-TB cases will help to control transmission of MDR-TB in Nepal.

The relationship between transmission time and clustering methods in Mycobacterium tuberculosis epidemiology

BACKGROUND

Tracking recent transmission is a vital part of controlling widespread pathogens such as Mycobacterium tuberculosis. Multiple methods with specific performance characteristics exist for detecting recent transmission chains, usually by clustering strains based on genotype similarities. With such a large variety of methods available, informed selection of an appropriate approach for determining transmissions within a given setting/time period is difficult.

METHODS

This study combines whole genome sequence (WGS) data derived from 324 isolates collected 2005-2010 in Kinshasa, Democratic Republic of Congo (DRC), a high endemic setting, with phylodynamics to unveil the timing of transmission events posited by a variety of standard genotyping methods. Clustering data based on Spoligotyping, 24-loci MIRU-VNTR typing, WGS based SNP (Single Nucleotide Polymorphism) and core genome multi locus sequence typing (cgMLST) typing were evaluated.

FINDINGS

Our results suggest that clusters based on Spoligotyping could encompass transmission events that occurred almost 200 years prior to sampling while 24-loci-MIRU-VNTR often represented three decades of transmission. Instead, WGS based genotyping applying low SNP or cgMLST allele thresholds allows for determination of recent transmission events, e.g. in timespans of up to 10 years for a 5 SNP/allele cut-off.

INTERPRETATION

With the rapid uptake of WGS methods in surveillance and outbreak tracking, the findings obtained in this study can guide the selection of appropriate clustering methods for uncovering relevant transmission chains within a given time-period. For high resolution cluster analyses, WGS-SNP and cgMLST based analyses have similar clustering/timing characteristics even for data obtained from a high incidence setting.

Global variation in bacterial strains that cause tuberculosis disease: a systematic review and meta-analysis

BACKGROUND

The host, microbial, and environmental factors that contribute to variation in tuberculosis (TB) disease are incompletely understood. Accumulating evidence suggests that one driver of geographic variation in TB disease is the local ecology of mycobacterial genotypes or strains, and there is a need for a comprehensive and systematic synthesis of these data. The objectives of this study were to (1) map the global distribution of genotypes that cause TB disease and (2) examine whether any epidemiologically relevant clinical characteristics were associated with those genotypes.

METHODS

We performed a systematic review of PubMed and Scopus to create a comprehensive dataset of human TB molecular epidemiology studies that used representative sampling techniques. The methods were developed according to the Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA). We extracted and synthesized data from studies that reported prevalence of bacterial genotypes and from studies that reported clinical characteristics associated with those genotypes.

RESULTS

The results of this study are twofold. First, we identified 206 studies for inclusion in the study, representing over 200,000 bacterial isolates collected over 27 years in 85 countries. We mapped the genotypes and found that, consistent with previously published maps, Euro-American lineage 4 and East Asian lineage 2 strains are widespread, and West African lineages 5 and 6 strains are geographically restricted. Second, 30 studies also reported transmission chains and 4 reported treatment failure associated with genotypes. We performed a meta-analysis and found substantial heterogeneity across studies. However, based on the data available, we found that lineage 2 strains may be associated with increased risk of transmission chains, while lineages 5 and 6 strains may be associated with reduced risk, compared with lineage 4 strains.

CONCLUSIONS

This study provides the most comprehensive systematic analysis of the evidence for diversity in bacterial strains that cause TB disease. The results show both geographic and epidemiological differences between strains, which could inform our understanding of the global burden of TB. Our findings also highlight the challenges of collecting the clinical data required to inform TB diagnosis and treatment. We urge future national TB programs and research efforts to prioritize and reinforce clinical data collection in study designs and results dissemination.

Compensatory evolution drives multidrug-resistant tuberculosis in Central Asia

Bacterial factors favoring the unprecedented multidrug-resistant tuberculosis (MDR-TB) epidemic in the former Soviet Union remain unclear. We utilized whole genome sequencing and Bayesian statistics to analyze the evolutionary history, temporal emergence of resistance and transmission networks of MDR Mycobacterium tuberculosis complex isolates from Karakalpakstan, Uzbekistan (2001-2006). One clade (termed Central Asian outbreak, CAO) dating back to 1974 (95% HPD 1969-1982) subsequently acquired resistance mediating mutations to eight anti-TB drugs. Introduction of standardized WHO-endorsed directly observed treatment, short-course in Karakalpakstan in 1998 likely selected for CAO-strains, comprising 75% of sampled MDR-TB isolates in 2005/2006. CAO-isolates were also identified in a published cohort from Russia (2008-2010). Similarly, the presence of mutations supposed to compensate bacterial fitness deficits was associated with transmission success and higher drug resistance rates. The genetic make-up of these MDR-strains threatens the success of both empirical and standardized MDR-TB therapies, including the newly WHO-endorsed short MDR-TB regimen in Uzbekistan.

Insights into the processes that drive the evolution of drug resistance in Mycobacterium tuberculosis

At present, the successful transmission of drug-resistant Mycobacterium tuberculosis, including multidrug-resistant (MDR) and extensively drug-resistant (XDR) strains, in human populations, threatens tuberculosis control worldwide. Differently from many other bacteria, M. tuberculosis drug resistance is acquired mainly through mutations in specific drug resistance-associated genes. The panel of mutations is highly diverse, but depends on the affected gene and M. tuberculosis genetic background. The variety of genetic profiles observed in drug-resistant clinical isolates underlines different evolutionary trajectories towards multiple drug resistance, although some mutation patterns are prominent. This review discusses the intrinsic processes that may influence drug resistance evolution in M. tuberculosis, such as mutation rate, drug resistance-associated mutations, fitness cost, compensatory mutations and epistasis. This knowledge should help to better predict the risk of emergence of highly resistant M. tuberculosis strains and to develop new tools and strategies to limit the development and spread of MDR and XDR strains.

THE DETECTION OF EPIDEMIC SUBTYPES OF BEIJING GENOTYPE OF MYCOBACTERIUM TUBERCULOSIS CIRCULATED IN THE PRIMORSKY KRAI

Background. The Far East is the territory with high rate of incidence and prevalence of tuberculosis. Cases of tuberculosis caused by epidemic strains have high frequency of MDR and XDR. It is important to study the prevalence of TB in areas with a high burden of infection, to which the Far East belongs. The aim of the research is to carry out genotyping of strains and assess the prevalence of CC1 and CC2 subtypes in the territory of Primorsky Krai. Materials and methods. The DNAs of 99 clinical isolates of MBT from Primorsky Krai have been genotyped by the 24-locus MIRU-VNTR and RD105/RD207. Results. The dominant number of strains pertained to Beijing genotype (59.6 %). The express method revealed 22 isolates of the CC2/W148 subtype, which had 6 different MIRU-VNTR-24 profile. According to MLVA classification MtbC 15-9, the most common among the isolates of CC2/W148 profile is 100-32 (59.1 %). Among these profiles the highest frequency of MDR/XDR was recorded – 69,2 %. According to the results of the express analysis, 39 isolates with 26 different MIRU-VNTR-24 profiles belonged to the CC1 subtype, of which the dominant number belonged to 99-32 and 94-32. Conclusions. The methods of express genotyping of epidemic subtypes of the Beijing genotype are very important for epidemiological surveillance and clinical practice. The developed methods allow to define a wider range of strains than previously used methods.

In vivo evolution of drug-resistant Mycobacterium tuberculosis in patients during long-term treatment

Background

In the current scenario, the drug-resistant tuberculosis is a significant challenge in the control of tuberculosis worldwide. In order to investigate the in vivo evolution of drug-resistant M. tuberculosis, the present study envisaged sequencing of the draft genomes of 18 serial isolates from four pre-extensively drug-resistant (pre-XDR) tuberculosis patients for continuous genetic alterations.

Results

All of the isolates harbored single nucleotide polymorphisms (SNPs) ranging from 1303 to 1309 with M. tuberculosis H37Rv as the reference. SNPs ranged from 0 to 12 within patients. The evolution rates were higher than the reported SNPs of 0.5 in the four patients. All the isolates exhibited mutations at sites of known drug targets, while some contained mutations in uncertain drug targets including folC, proZ, and pyrG. The compensatory substitutions for rescuing these deleterious mutations during evolution were only found in RpoC I491T in one patient. Many loci with microheterogeneity showed transient mutations in different isolates. Ninety three SNPs exhibited significant association with refractory pre-XDR TB isolates.

Conclusions

Our results showed evolutionary changes in the serial genetic characteristics of the pre-XDR TB patients due to accumulation of the fixed drug-resistant related mutations, and the transient mutations under continuous antibiotics pressure over several years.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5010-5) contains supplementary material, which is available to authorized users.

Harmonized Genome Wide Typing of Tubercle Bacilli Using a Web-Based Gene-By-Gene Nomenclature System

Background

Global tuberculosis (TB) control is challenged by uncontrolled transmission of Mycobacterium tuberculosis complex (Mtbc) strains, esp. of multidrug (MDR) or extensively resistant (XDR) variants. Precise analysis of transmission networks is the basis to trace outbreak M/XDR clones and improve TB control. However, classical genotyping tools lack discriminatory power due to the high similarity of strains of particular successful lineages, e.g. Beijing or outbreak strains. This can be overcome by whole genome sequencing (WGS) approaches, but these are not yet standardized to facilitate larger investigations encompassing different laboratories or outbreak tracing across borders.

Methods

We established and improved a whole genome gene-by-gene multi locus sequence typing approach encompassing a stable set of core genome genes (cgMLST) and linked it to a web-based nomenclature server (cgMLST.org) facilitating assignment and storage of allele numbers.

Findings

We evaluated and refined a previously suggested cgMLST schema by using a reference strain set (n = 251) reflecting the global diversity of the Mtbc. A set of 2891 genes showed excellent performance with at least 97% of the genes reliably identified in strains of all Mtbc lineages and in discriminating outbreak strains. cgMLST allele numbers were automatically retrieved from and stored at cgMLST.org.

Interpretation

The refined cgMLST schema provides high resolution genome-based typing of clinical strains of all Mtbc lineages. Combined with a web-based nomenclature server, it facilitates rapid, high-resolution, and harmonized tracing of clinical Mtbc strains needed for prospective local and global surveillance.

Genetic diversity of Mycobacterium tuberculosis complex isolates circulating in an area with high tuberculosis incidence: Using 24-locus MIRU-VNTR method

We aimed to determine the genetic diversity, phylogenetic relationship and transmission dynamics of Mycobacterium tuberculosis complex (MTBC) genotypes in an area with high tuberculosis (TB) incidence. A set of 164 MTBC isolates from new TB patients of Golestan province, Iran, were subjected to genotyping using the standard 24-locus MIRU-VNTR method. Recent TB transmission was evaluated and phylogenetic relationships were analysed by minimum spanning tree and cluster-graph methods. Among the 164 isolates, 132 distinct patterns were detected. The 48 clustered isolates (29.3%) were distributed into 16 clusters ranging in size from 2 to 12 isolates. The most frequent genotype was Central Asian Strain/Delhi (CAS/Delhi) (n = 67, 40.8%), followed by NEW-1 (n = 53, 32.3%) and Beijing (n = 19, 11.6%) genotypes. Thirty five (72.9%) of NEW-1 isolates were recovered from immigrant patients and 84.2% (n = 16) of Beijing genotypes recovered from native cases. Statistically significant association was found between clustering and smoking (p = 0.047), drug addiction (p = 0.01) and prison history (p = 0.003). The estimated proportion of recent transmission was 19.5%. Presence of highly diverse MTBC isolates circulating in this province without a dominant genotype might be a consequence of importation of various genotypes in this area.

A Quantitative Evaluation of MIRU-VNTR Typing Against Whole-Genome Sequencing for Identifying Mycobacterium tuberculosis Transmission: A Prospective Observational Cohort Study

Background

Mycobacterial Interspersed Repetitive Unit-Variable Number Tandem Repeat (MIRU-VNTR) typing is widely used in high-income countries to determine Mycobacterium tuberculosis relatedness. Whole-genome sequencing (WGS) is known to deliver greater specificity, but no quantitative prospective comparison has yet been undertaken.

Methods

We studied isolates from the English Midlands, sampled consecutively between 1 January 2012 and 31 December 2015. In addition to routinely performed MIRU-VNTR typing, DNA was extracted from liquid cultures and sequenced using Illumina technology. Demographic and epidemiological data for the relevant patients were extracted from the Enhanced Tuberculosis Surveillance system run by Public Health England. Closely related samples, defined using a threshold of five single nucleotide variants (SNVs), were compared to samples with identical MIRU-VNTR profiles, to samples from individuals with shared epidemiological risk factors, and to those with both characteristics.

Findings

1999 patients were identified for whom at least one M. tuberculosis isolate had been MIRU-VNTR typed and sequenced. Comparing epidemiological risk factors with close genetic relatedness, only co-residence had a positive predictive value of over 5%. Excluding co-resident individuals, 18.6% of patients with identical MIRU-VNTR profiles were within 5 SNVs. Where patients also shared social risk factors and ethnic group, this rose to 48%. Only 8% of MIRU-VNTR linked pairs in lineage 1 were within 5 SNV, compared to 31% in lineage 4.

Interpretation

In the setting studied, this molecular epidemiological study shows MIRU-VNTR typing and epidemiological risk factors are poorly predictive of close genomic relatedness, assessed by SNV. MIRU-VNTR performance varies markedly by lineage.

Funding

Public Health England, Health Innovation Challenge Fund, NIHR Health Protection Research Unit Oxford, NIHR Oxford Biomedical Research Centre.

•

M. tuberculosis relatedness was measured by MIRU-VNTR typing and single nucleotide variants (SNV).

•

In consecutive isolates in England, only 20% of samples with identical MIRU-VNTR profiles were within 5 SNVs.

•

MIRU-VNTR is more predictive of close relatedness in lineage 4 (Euro-American) lineages than in other lineages.

Recombination Signal in Mycobacterium tuberculosis Stems from Reference-guided Assemblies and Alignment Artefacts

DNA acquisition via genetic recombination is considered advantageous as it has the potential to bring together beneficial mutations that emerge independently within a population. Furthermore, recombination is considered to contribute to the maintenance of genome stability by purging slightly deleterious mutations. The prevalence of recombination differs among prokaryotic species and depends on the accessibility of DNA transfer mechanisms. An exceptional example is the human pathogen Mycobacterium tuberculosis (MTB) where no clear transfer mechanisms have been so far characterized and the presence of recombination is questioned. Here, we analyze completely assembled MTB genomes in search for evidence of recombination. We find that putative recombination events are enriched in strains reconstructed by reference-guided assembly and in regions with unreliable alignments. In addition, assembly and alignment artefacts introduce phylogenetic signals that are conflicting the established MTB phylogeny. Our results reveal that the so far reported recombination events in MTB are likely to stem from methodological artefacts. We conclude that no reliable signal of recombination is observed in the currently available MTB genomes. Moreover, our study demonstrates the limitations of reference-guided genome assembly for phylogenetic reconstructions. Rigorously de novo assembled genomes of high quality are mandatory in order to distinguish true evolutionary signal from noise, in particular for low diversity species such as MTB.

Geospatial distribution of Mycobacterium tuberculosis genotypes in Africa

Objective

To investigate the distribution of Mycobacterium tuberculosis genotypes across Africa.

Methods

The SITVIT2 global repository and PUBMED were searched for spoligotype and published genotype data respectively, of M. tuberculosis from Africa. M. tuberculosis lineages in Africa were described and compared across regions and with those from 7 European and 6 South-Asian countries. Further analysis of the major lineages and sub-lineages using Principal Component analysis (PCA) and hierarchical cluster analysis were done to describe clustering by geographical regions. Evolutionary relationships were assessed using phylogenetic tree analysis.

Results

A total of 14727 isolates from 35 African countries were included in the analysis and of these 13607 were assigned to one of 10 major lineages, whilst 1120 were unknown. There were differences in geographical distribution of major lineages and their sub-lineages with regional clustering. Southern African countries were grouped based on high prevalence of LAM11-ZWE strains; strains which have an origin in Portugal. The grouping of North African countries was due to the high percentage of LAM9 strains, which have an origin in the Eastern Mediterranean region. East African countries were grouped based on Central Asian (CAS) and East-African Indian (EAI) strain lineage possibly reflecting historic sea trade with Asia, while West African Countries were grouped based on Cameroon lineage of unknown origin. A high percentage of the Haarlem lineage isolates were observed in the Central African Republic, Guinea, Gambia and Tunisia, however, a mixed distribution prevented close clustering.

Conclusions

This study highlighted that the TB epidemic in Africa is driven by regional epidemics characterized by genetically distinct lineages of M. tuberculosis. M. tuberculosis in these regions may have been introduced from either Europe or Asia and has spread through pastoralism, mining and war. The vast array of genotypes and their associated phenotypes should be considered when designing future vaccines, diagnostics and anti-TB drugs.

Whole genome SNP analysis suggests unique virulence factor differences of the Beijing and Manila families of Mycobacterium tuberculosis found in Hawaii

While tuberculosis (TB) remains a global disease, the WHO estimates that 62% of the incident TB cases in 2016 occurred in the WHO South-East Asia and Western Pacific regions. TB in the Pacific is composed predominantly of two genetic families of Mycobacterium tuberculosis (Mtb): Beijing and Manila. The Manila family is historically under-studied relative to the families that comprise the majority of TB in Europe and North America (e.g. lineage 4), and it remains unclear why this lineage has persisted in Filipino populations despite the predominance of more globally successful Mtb lineages in most of the world. The Beijing family is of particular interest as it is increasingly associated with drug resistance throughout the world. Both of these lineages are important to the State of Hawaii, where they comprise over two-thirds of TB cases. Here, we performed whole genome sequencing on 82 Beijing family, Manila family, and outgroup clinical Mtb isolates from Hawaii to identify lineage-specific SNPs (SNPs found in all isolates from their respective families, and exclusively in those families) in established virulence factor genes. Six non-silent lineage-specific virulence factor SNPs were found in the Beijing family, including mutations in alternative sigma factor sigG and polyketide synthases pks5 and pks7. The Manila family displayed more than eleven non-silent lineage-specific and characteristic virulence factor mutations, including in genes coding for MCE-family protein Mce1B, two mutations in fatty-acid-AMP ligase FadD26, and virulence-regulating transcriptional regulator VirS. This study further identified an ancient clade that shared some virulence factor mutations with the Manila family, and investigated the relationship of those and other “Manila-like” spoligotypes to the Manila family with this SNP dataset. This work identified a set of virulence genes that are worth pursuing to determine potential differences in transmission or virulence displayed by these two Mtb families.

MIRU-profiler: a rapid tool for determination of 24-loci MIRU-VNTR profiles from assembled genomes of Mycobacterium tuberculosis

Background

Tuberculosis (TB) resulted in an estimated 1.7 million deaths in the year 2016. The disease is caused by the members of Mycobacterium tuberculosis complex, which includes Mycobacterium tuberculosis, Mycobacterium bovis and other closely related TB causing organisms. In order to understand the epidemiological dynamics of TB, national TB control programs often conduct standardized genotyping at 24 Mycobacterial-Interspersed-Repetitive-Units (MIRU)-Variable-Number-of-Tandem-Repeats (VNTR) loci. With the advent of next generation sequencing technology, whole-genome sequencing (WGS) has been widely used for studying TB transmission. However, an open-source software that can connect WGS and MIRU-VNTR typing is currently unavailable, which hinders interlaboratory communication. In this manuscript, we introduce the MIRU-profiler program which could be used for prediction of MIRU-VNTR profile from WGS of M. tuberculosis.

Implementation

The MIRU-profiler is implemented in shell scripting language and depends on EMBOSS software. The in-silico workflow of MIRU-profiler is similar to those described in the laboratory manuals for genotyping M. tuberculosis. Given an input genome sequence, the MIRU-profiler computes alleles at the standard 24-loci based on in-silico PCR amplicon lengths. The final output is a tab-delimited text file detailing the 24-loci MIRU-VNTR pattern of the input sequence.

Validation

The MIRU-profiler was validated on four datasets: complete genomes from NCBI-GenBank (n = 11), complete genomes for locally isolated strains sequenced using PacBio (n = 4), complete genomes for BCG vaccine strains (n = 2) and draft genomes based on 250 bp paired-end Illumina reads (n = 106).

Results

The digital MIRU-VNTR results were identical to the experimental genotyping results for complete genomes of locally isolated strains, BCG vaccine strains and five out of 11 genomes from the NCBI-GenBank. For draft genomes based on short Illumina reads, 21 out of 24 loci were inferred with a high accuracy, while a number of inaccuracies were recorded for three specific loci (ETRA, QUB11b and QUB26). One of the unique features of the MIRU-profiler was its ability to process multiple genomes in a batch. This feature was tested on all complete M. tuberculosis genome (n = 157), for which results were successfully obtained in approximately 14 min.

Conclusion

The MIRU-profiler is a rapid tool for inference of digital MIRU-VNTR profile from the assembled genome sequences. The tool can accurately infer repeat numbers at the standard 24 or 21/24 MIRU-VNTR loci from the complete or draft genomes respectively. Thus, the tool is expected to bridge the communication gap between the laboratories using WGS and those using the conventional MIRU-VNTR typing.

Insight into multidrug-resistant Beijing genotype Mycobacterium tuberculosis isolates in Myanmar

OBJECTIVES

Myanmar is a World Health Organization high tuberculosis (TB) burden country with a high multidrug-resistant (MDR)-TB burden. Of significance, a high prevalence of the Beijing genotype of Mycobacterium tuberculosis (MTB) among MDR-MTB has been reported previously. A detailed genetic characterization of TB clinical isolates was performed in order to explore whether there is an association between the prevalence of the Beijing MTB genotype and MDR-TB in Myanmar.

METHODS

A total of 265 MDR-MTB clinical isolates collected in 2010 and 2012 were subjected to spoligotyping, mycobacterial interspersed repetitive unit-variable number tandem repeat (MIRU-VNTR) analysis, single nucleotide polymorphism (SNP) typing, and drug resistance-associated gene sequencing, including rpoC to detect potential compensatory evolution.

RESULTS

Of the total MDR-MTB isolates, 79.2% (210/265) were of the Beijing genotype, the majority of which were the 'modern' subtype. Beijing genotype isolates were differentiated by 15-locus MIRU-VNTR and a high clustering rate (53.0%) was observed in the modern subtype. These MIRU-VNTR patterns were similar to Beijing genotype clones spreading across Russia and Central Asia. A high prevalence of katG Ser315Thr, and genetic evidence of extensive drug resistance (XDR) and pre-XDR and compensatory mutations in rpoC were observed among clustered isolates.

CONCLUSIONS

MDR-MTB strains of the Beijing genotype might be spreading in Myanmar and present a major challenge to TB control in this country.

Inferring demographic parameters in bacterial genomic data using Bayesian and hybrid phylogenetic methods

Background

Recent developments in sequencing technologies make it possible to obtain genome sequences from a large number of isolates in a very short time. Bayesian phylogenetic approaches can take advantage of these data by simultaneously inferring the phylogenetic tree, evolutionary timescale, and demographic parameters (such as population growth rates), while naturally integrating uncertainty in all parameters. Despite their desirable properties, Bayesian approaches can be computationally intensive, hindering their use for outbreak investigations involving genome data for a large numbers of pathogen isolates. An alternative to using full Bayesian inference is to use a hybrid approach, where the phylogenetic tree and evolutionary timescale are estimated first using maximum likelihood. Under this hybrid approach, demographic parameters are inferred from estimated trees instead of the sequence data, using maximum likelihood, Bayesian inference, or approximate Bayesian computation. This can vastly reduce the computational burden, but has the disadvantage of ignoring the uncertainty in the phylogenetic tree and evolutionary timescale.

Results

We compared the performance of a fully Bayesian and a hybrid method by analysing six whole-genome SNP data sets from a range of bacteria and simulations. The estimates from the two methods were very similar, suggesting that the hybrid method is a valid alternative for very large datasets. However, we also found that congruence between these methods is contingent on the presence of strong temporal structure in the data (i.e. clocklike behaviour), which is typically verified using a date-randomisation test in a Bayesian framework. To reduce the computational burden of this Bayesian test we implemented a date-randomisation test using a rapid maximum likelihood method, which has similar performance to its Bayesian counterpart.

Conclusions

Hybrid approaches can produce reliable inferences of evolutionary timescales and phylodynamic parameters in a fraction of the time required for fully Bayesian analyses. As such, they are a valuable alternative in outbreak studies involving a large number of isolates.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1210-5) contains supplementary material, which is available to authorized users.

RD5-mediated lack of PE_PGRS and PPE-MPTR export in BCG vaccine strains results in strong reduction of antigenic repertoire but little impact on protection

Tuberculosis is the deadliest infectious disease worldwide. Although the BCG vaccine is widely used, it does not efficiently protect against pulmonary tuberculosis and an improved tuberculosis vaccine is therefore urgently needed. Mycobacterium tuberculosis uses different ESX/Type VII secretion (T7S) systems to transport proteins important for virulence and host immune responses. We recently reported that secretion of T7S substrates belonging to the mycobacteria-specific Pro-Glu (PE) and Pro-Pro-Glu (PPE) proteins of the PGRS (polymorphic GC-rich sequences) and MPTR (major polymorphic tandem repeat) subfamilies required both a functional ESX-5 system and a functional PPE38/71 protein for secretion. Inactivation of ppe38/71 and the resulting loss of PE_PGRS/PPE-MPTR secretion were linked to increased virulence of M. tuberculosis strains. Here, we show that a predicted total of 89 PE_PGRS/PPE-MPTR surface proteins are not exported by certain animal-adapted strains of the M. tuberculosis complex including M. bovis. This Δppe38/71-associated secretion defect therefore also occurs in the M. bovis-derived tuberculosis vaccine BCG and could be partially restored by introduction of the M. tuberculosis ppe38-locus. Epitope mapping of the PPE-MPTR protein PPE10, further allowed us to monitor T-cell responses in splenocytes from BCG/M. tuberculosis immunized mice, confirming the dependence of PPE10-specific immune-induction on ESX-5/PPE38-mediated secretion. Restoration of PE_PGRS/PPE-MPTR secretion in recombinant BCG neither altered global antigenic presentation or activation of innate immune cells, nor protective efficacy in two different mouse vaccination-infection models. This unexpected finding stimulates a reassessment of the immunomodulatory properties of PE_PGRS/PPE-MPTR proteins, some of which are contained in vaccine formulations currently in clinical evaluation.

One of the major findings of the pioneering Mycobacterium tuberculosis H37Rv genome sequencing project was the identification of the highly abundant PE and PPE proteins, named after their N-terminal motifs Pro–Glu (PE) or Pro–Pro–Glu (PPE). Within the 20 years of research since then, many claims were made that PE/PPE proteins, including the two large subgroups encoded by repetitive sequences with very high GC content (PE_PGRS and PPE-MPTR families), are exported to the bacterial surface or beyond, and show broad immunomodulatory impact on host-pathogen interaction. We thus screened strains from different branches of the M. tuberculosis complex, including the attenuated Mycobacterium bovis BCG vaccine strains, for their capacity to export PE_PGRS/PPE-MPTR proteins. Strikingly, we found that BCG strains were unable to export the plethora of PE_PGRS/PPE-MPTR proteins due to the absence of the region of difference RD5, which in M. tuberculosis encodes PPE38, required for PE_PGRS/PPE-MPTR export. Surprisingly, the restoration of PE_PGRS/PPE-MPTR export by genetic complementation in recombinant BCG did not result in immunomodulatory changes or altered protection in mouse models. Our results thus put into perspective the numerous reports on virulence-associated immunomodulatory impact of individual PE_PGRS and PPE-MPTR proteins and open novel questions on their biological function(s).

A complete high-quality MinION nanopore assembly of an extensively drug-resistant Mycobacterium tuberculosis Beijing lineage strain identifies novel variation in repetitive PE/PPE gene regions

A better understanding of the genomic changes that facilitate the emergence and spread of drug-resistant Mycobacterium tuberculosis strains is currently required. Here, we report the use of the MinION nanopore sequencer (Oxford Nanopore Technologies) to sequence and assemble an extensively drug-resistant (XDR) isolate, which is part of a modern Beijing sub-lineage strain, prevalent in Western Province, Papua New Guinea. Using 238-fold coverage obtained from a single flow-cell, de novo assembly of nanopore reads resulted into one contiguous assembly with 99.92 % assembly accuracy. Incorporation of complementary short read sequences (Illumina) as part of consensus error correction resulted in a 4 404 064 bp genome with 99.98 % assembly accuracy. This assembly had an average nucleotide identity of 99.7 % relative to the reference genome, H37Rv. We assembled nearly all GC-rich repetitive PE/PPE family genes (166/168) and identified variants within these genes. With an estimated genotypic error rate of 5.3 % from MinION data, we demonstrated identification of variants to include the conventional drug resistance mutations, and those that contribute to the resistance phenotype (efflux pumps/transporter) and virulence. Reference-based alignment of the assembly allowed detection of deletions and insertions. MinION sequencing provided a fully annotated assembly of a transmissible XDR strain from an endemic setting and showed its utility to provide further understanding of genomic processes within Mycobacterium tuberculosis.

Single nucleotide polymorphisms of Beijing lineage Mycobacterium tuberculosis toxin-antitoxin system genes: Their role in the changes of protein activity and evolution

The article investigates SNP in genes of toxin-antitoxin systems type II in Mycobacterium tuberculosis Beijing lineage strains and their possible role in the development and formation of new sublineages. We established the catalog of SNPs in 142 TA systems genes in 1349 sequenced genomes of the M. tuberculosis Beijing lineage. Based on the catalog, 15 new sublineages were identified as part of Beijing lineages by non-synonymous SNP in 21 genes of TA systems. We discovered three toxin genes with mutations specific for epidemiologically dangerous sublineages Beijing-modern (vapC37 A46G, vapC38 T143C) and Beijing-B0/W148 (vapC12 A95G). We proved the functional significance of these polymorphisms by cloning these genes wild-type and with marker mutations for the Beijing lineage vapC12 (A95G), vapC37 (A46G), vapC38 (T143C). In vitro study of their activities revealed effect of mutations on the RNase activity of toxin proteins. Mutations in vapC37 and vapC38 decreased toxin activity, and mutation in the vapC12 increased it. We cloned the toxin vapC37 gene of Mycobacterium smegmatis mc² 155 in both allelic variants: without mutation and with A46G mutation, specific for the Beijing-modern lineage. It was shown that this mutation leads to a loss of toxicity.

Frequent transmission of the Mycobacterium tuberculosis Beijing lineage and positive selection for the EsxW Beijing variant in Vietnam

To examine the transmission dynamics of Mycobacterium tuberculosis (Mtb) isolated from tuberculosis patients in Ho Chi Minh City, Vietnam, we sequenced the whole genomes of 1,635 isolates and compared these with 3,144 isolates from elsewhere. The data identify an underlying burden of disease caused by the endemic Mtb lineage 1 associated with the activation of long-term latent infection, and a threefold higher burden associated with the more recently introduced Beijing lineage and lineage 4 Mtb strains. We find that Beijing lineage Mtb is frequently transferred between Vietnam and other countries, and detect higher levels of transmission of Beijing lineage strains within this host population than the endemic lineage 1 Mtb. Screening for parallel evolution of Beijing lineage-associated SNPs in other Mtb lineages as a signal of positive selection, we identify an alteration in the ESX-5 type VII-secreted protein EsxW, which could potentially contribute to the enhanced transmission of Beijing lineage Mtb in Vietnamese and other host populations.

Frequent transmission of the Mycobacterium tuberculosis Beijing lineage and positive selection for the EsxW Beijing variant in Vietnam

To examine the transmission dynamics of Mycobacterium tuberculosis (Mtb) isolated from tuberculosis patients in Ho Chi Minh City, Vietnam, we sequenced the whole genomes of 1,635 isolates and compared these with 3,144 isolates from elsewhere. The data identify an underlying burden of disease caused by the endemic Mtb lineage 1 associated with the activation of long-term latent infection, and a threefold higher burden associated with the more recently introduced Beijing lineage and lineage 4 Mtb strains. We find that Beijing lineage Mtb is frequently transferred between Vietnam and other countries, and detect higher levels of transmission of Beijing lineage strains within this host population than the endemic lineage 1 Mtb. Screening for parallel evolution of Beijing lineage-associated SNPs in other Mtb lineages as a signal of positive selection, we identify an alteration in the ESX-5 type VII-secreted protein EsxW, which could potentially contribute to the enhanced transmission of Beijing lineage Mtb in Vietnamese and other host populations.

Importance of differential identification of Mycobacterium tuberculosis strains for understanding differences in their prevalence, treatment efficacy, and vaccine development

Tuberculosis (TB), caused by Mycobacterium tuberculosis (Mtb), remains a serious global health problem in the 21^st century because of its high mortality. Mtb is an extremely successful human-adapted pathogen that displays a multifactorial ability to control the host immune response and to evade killing by drugs, resulting in the breakdown of BCG vaccine-conferred anti-TB immunity and development of multidrug-resistant (MDR) and extensively drug-resistant (XDR) Mtb. Although genetic components of the genomes of the Mtb complex strains are highly conserved, showing over 99% similarity to other bacterial genera, recently accumulated evidence suggests that the genetic diversity of the Mtb complex strains has implications for treatment outcomes, development of MDR/XDR Mtb, BCG vaccine efficacy, transmissibility, and epidemiological outbreaks. Thus, new insights into the pathophysiological features of the Mtb complex strains are required for development of novel vaccines and for control of MDR/XDR Mtb infection, eventually leading to refinement of treatment regimens and the health care system. Many studies have focused on the differential identification of Mtb complex strains belonging to different lineages because of differences in their virulence and geographical dominance. In this review, we discuss the impact of differing genetic characteristics among Mtb complex strains on vaccine efficacy, treatment outcome, development of MDR/XDR Mtb strains, and epidemiological outbreaks by focusing on the best-adapted human Mtb lineages. We further explore the rationale for differential identification of Mtb strains for more effective control of TB in clinical and laboratory settings by scrutinizing current diagnostic methods.

Genetic profiling of Mycobacterium bovis strains from slaughtered cattle in Eritrea

Mycobacterium bovis (M.bovis) is the main causative agent for bovine tuberculosis (BTB) and can also be the cause of zoonotic tuberculosis in humans. In view of its zoonotic nature, slaughterhouse surveillance, potentially resulting in total or partial condemnation of the carcasses and organs, is conducted routinely. Spoligotyping, VNTR profiling, and whole genome sequencing (WGS) of M. bovis isolated from tissues with tuberculosis-like lesions collected from 14 cattle at Eritrea’s largest slaughterhouse in the capital Asmara, were conducted.The 14 M. bovis isolates were classified into three different spoligotype patterns (SB0120, SB0134 and SB0948) and six VNTR profiles. WGS results matched those of the conventional genotyping methods and further discriminated the six VNTR profiles into 14 strains. Furthermore, phylogenetic analysis of the M. bovis isolates suggests two independent introductions of BTB into Eritrea possibly evolving from a common ancestral strain in Europe.This molecular study revealed the most important strains of M. bovis in Eritrea and their (dis)similarities with the strains generally present in East Africa and Europe, as well as potential routes of introduction of M. bovis. Though the sample size is small, the current study provides important information as well as platform for future in-depth molecular studies on isolates from both the dairy and the traditional livestock sectors in Eritrea and the region. This study provides information onthe origin of some of the M. bovis strains in Eritrea, its genetic diversity, evolution and patterns of spread between dairy herds. Such information is essential in the development and implementation of future BTB control strategy for Eritrea.

The livestock sector plays a major role in poverty and hunger reduction in the vast majority of Africa, as a source of food, cash income, manure, draught power, transportation, savings, insurance and social status. However, for livestock to play this vital role, the impact of diseases of economic and zoonotic importance need to be reduced. Bovine tuberculosis, mainly caused by Mycobacterium bovis, is such an infectious disease. Slaughterhouse (gross pathology) surveillance, followed by bacterial culture and genotyping, are options to identify the disease-causing agents, their distribution, and enabling trace back of the sources of infections, in order to prevent their re-introduction and spread. Unfortunately, genotyping is by far not generally introduced in the continent. In the present study, tissues with tuberculosis-like lesions were collected from the Asmara municipal slaughterhouse, the largest slaughterhouse in Eritrea, and bacterial culture, classical Mycobacterium tuberculosis complex typing (Spoligotyping and VNTR profiling), as well as whole genome sequencing (WGS) were used to gain insight into the spatial and temporal distribution, genetic diversity and evolution of M. bovis strains circulating in Eritrean dairy cattle. The results revealed (dis)similarities of the Eritrean M. bovis strains with the strains generally present in Africa and Europe, potential routes of introduction to Eritrea and genetic diversity of the M. bovis strains. Future in-depth molecular studies including more samples from dairy cattle as well as cattle and goats from the traditional livestock sector are recommended.

New insights into the transposition mechanisms of IS6110 and its dynamic distribution between Mycobacterium tuberculosis Complex lineages

The insertion Sequence IS6110, only present in the pathogens of the Mycobacterium tuberculosis Complex (MTBC), has been the gold-standard epidemiological marker for TB for more than 25 years, but biological implications of IS6110 transposition during MTBC adaptation to humans remain elusive. By studying 2,236 clinical isolates typed by IS6110-RFLP and covering the MTBC, we remarked a lineage-specific content of IS6110 being higher in modern globally distributed strains. Once observed the IS6110 distribution in the MTBC, we selected representative isolates and found a correlation between the normalized expression of IS6110 and its abundance in MTBC chromosomes. We also studied the molecular regulation of IS6110 transposition and we found a synergistic action of two post-transcriptional mechanisms: a -1 ribosomal frameshift and a RNA pseudoknot which interferes translation. The construction of a transcriptionally active transposase resulted in 20-fold increase of the transposition frequency. Finally, we examined transposition in M. bovis and M. tuberculosis during laboratory starvation and in a mouse infection model of TB. Our results shown a higher transposition in M. tuberculosis, that preferably happens during TB infection in mice and after one year of laboratory culture, suggesting that IS6110 transposition is dynamically adapted to the host and to adverse growth conditions.

Since the pioneering discovery of transposition by Barbara McClintock in eukaryotes and later in prokaryotes by Robert W. Hedges and Alan E. Jacob, it has become clear the key role of mobile genetics elements in chromosome remodelling, microbial evolution and host adaptation. The insertion sequence IS6110 is widely recognized for its utility in TB diagnosis and epidemiology because it is only present in the M. tuberculosis Complex (MTBC) and its transposition provides an excellent chromosomal polymorphic variability allowing the study of recent TB transmission. This inherent feature of IS6110 leads us to hypothesize that IS6110 plays a crucial role during the TB infectious cycle. However, the biological significance of IS6110 has been hindered by its almost exclusive use as an epidemiological marker. Here, we study the regulatory mechanisms and the distribution of IS6110 in the different MTBC lineages. We discuss the potential biological implications of IS6110, that is much more than an excellent TB epidemiological tool. Since IS6110 could play an important role in the adaptation of MTBC to the host, this study opens new avenues to decipher the biological roles of IS6110 in TB pathogenesis.

Changing patterns of human migrations shaped the global population structure of Mycobacterium tuberculosis in France

Mycobacterium tuberculosis (Mtb) exhibits a structured phylogeographic distribution worldwide linked with human migrations. We sought to infer how the interactions between distinct human populations shape the global population structure of Mtb on a regional scale. We applied the recently described timescaled haplotypic density (THD) technique on 638 minisatellite-based Mtb genotypes from French tuberculosis patients. THD with a long-term (200 y) timescale indicated that Mtb population in France had been mostly influenced by interactions with Eastern and Southern Europe and, to a lesser extent, Northern and Middle Africa, consistent with historical migrations favored by geographic proximity or commercial exchanges with former French colonies. Restricting the timescale to 20 y, THD identified a sustained influence of Northern Africa, but not Europe where tuberculosis incidence decreased sharply. Evolving interactions between human populations, thus, measurably influence the local population structure of Mtb. Relevant information on such interactions can be inferred using THD from Mtb genotypes.

Clonal expansion across the seas as seen through CPLP-TB database: A joint effort in cataloguing Mycobacterium tuberculosis genetic diversity in Portuguese-speaking countries

Tuberculosis (TB) remains a major health problem within the Community of Portuguese Language Speaking Countries (CPLP). Despite the marked variation in TB incidence across its member-states and continued human migratory flux between countries, a considerable gap in the knowledge on the Mycobacterium tuberculosis population structure and strain circulation between the countries still exists. To address this, we have assembled and analysed the largest CPLP M. tuberculosis molecular and drug susceptibility dataset, comprised by a total of 1447 clinical isolates, including 423 multidrug-resistant isolates, from five CPLP countries. The data herein presented reinforces Latin American and Mediterranean (LAM) strains as the hallmark of M. tuberculosis populational structure in the CPLP coupled with country-specific differential prevalence of minor clades. Moreover, using high-resolution typing by 24-loci MIRU-VNTR, six cross-border genetic clusters were detected, thus supporting recent clonal expansion across the Lusophone space.

To make this data available to the scientific community and public health authorities we developed CPLP-TB (available at http://cplp-tb.ff.ulisboa.pt), an online database coupled with web-based tools for exploratory data analysis. As a public health tool, it is expected to contribute to improved knowledge on the M. tuberculosis population structure and strain circulation within the CPLP, thus supporting the risk assessment of strain-specific trends.

•

The Community of Portuguese Speaking Countries (CPLP) occupies a vast geographical area.

•

Three CPLP countries are shortlisted in the WHO's list of Top 30 high-burden countries.

•

Common Mycobacterium tuberculosis population structure denote historical strain flow.

•

Cross-border clusters suggest recent intercontinental tuberculosis transmission.

•

CPLP-TB: a novel strain database and framework for collaborative studies and strain tracing.

TnSeq of Mycobacterium tuberculosis clinical isolates reveals strain-specific antibiotic liabilities

Once considered a phenotypically monomorphic bacterium, there is a growing body of work demonstrating heterogeneity among Mycobacterium tuberculosis (Mtb) strains in clinically relevant characteristics, including virulence and response to antibiotics. However, the genetic and molecular basis for most phenotypic differences among Mtb strains remains unknown. To investigate the basis of strain variation in Mtb, we performed genome-wide transposon mutagenesis coupled with next-generation sequencing (TnSeq) for a panel of Mtb clinical isolates and the reference strain H37Rv to compare genetic requirements for in vitro growth across these strains. We developed an analytic approach to identify quantitative differences in genetic requirements between these genetically diverse strains, which vary in genomic structure and gene content. Using this methodology, we found differences between strains in their requirements for genes involved in fundamental cellular processes, including redox homeostasis and central carbon metabolism. Among the genes with differential requirements were katG, which encodes the activator of the first-line antitubercular agent isoniazid, and glcB, which encodes malate synthase, the target of a novel small-molecule inhibitor. Differences among strains in their requirement for katG and glcB predicted differences in their response to these antimicrobial agents. Importantly, these strain-specific differences in antibiotic response could not be predicted by genetic variants identified through whole genome sequencing or by gene expression analysis. Our results provide novel insight into the basis of variation among Mtb strains and demonstrate that TnSeq is a scalable method to predict clinically important phenotypic differences among Mtb strains.

Biology of clinical strains of Mycobacterium tuberculosis with varying levels of transmission

Transmission of Mycobacterium tuberculosis bacilli from one individual to another is the basis of the disease process. While considerable emphasis has been placed on the role of host mechanisms of resistance in establishing or preventing new infection, far less has been expended on understanding possible factors operative at the bacterial level. In this study we established a panel of clinical isolates of M. tuberculosis strains obtained from the Western Cape region of South Africa, each of which had been carefully tracked in terms of their degree of transmission in the community. Each of the panel were used to infect guinea pigs with 15-20 bacilli by aerosol exposure and the course of the infection then determined. Strains with different degrees of transmission could not be distinguished in terms of their capacity to grow in the main target organs of infected animals. However, rather surprisingly, while strains with no evidence of transmission [NOT] in general caused moderate to severe lung damage, this parameter in animals infected with highly transmitted [HT] strains was mostly mild. In terms of TH1 immunity these signals were strongest in these latter animals, as was IL-17 gene expression, whereas minimal signals for regulatory molecules including IL-10 and FoxP3 were seen across the entire panel. In terms of T cell numbers, responses of both CD4 and CD8 were both far faster and far higher in animals infected with the HT strains. At the gene expression level we observed a major three-fold difference [both up and down] between NOT and HT strains, but in terms of proteins of key interest only a few [including PD-L1 and HIF-3] showed major differences between the two groups. Overall, it was apparent that NOT strains were far more inflammatory that HT strains, and appeared to trigger a much larger number of genes, possibly explaining the observed damage to the lungs and progressive pathology. In contrast, the HT strains, while equally virulent, were more immunogenic and developed much stronger T cell responses, while keeping lung damage to a minimum. Hence, in terms of trying to explain the capacity of these strains to cause transmission, these results are clearly paradoxical.

Mycobacterium tuberculosis complex genotypes circulating in Nigeria based on spoligotyping obtained from Ziehl-Neelsen stained slides extracted DNA

METHODS

All State TB control programmes in Nigeria were requested to submit 25-50 smear-positive Ziehl-Neelsen (ZN) stained slides for screening during 2013-2014. DNA was extracted from 929 slides for spoligotyping and drug-resistance analysis using microbead-based flow-cytometry suspension arrays.

RESULTS

Spoligotyping results were obtained for 549 (59.1%) of 929 samples. Lineage 4 Cameroon sublineage (L4.6.2) represented half of the patterns, Mycobacterium africanum (L5 and L6) represented one fifth of the patterns, and all other lineages, including other L4 sublineages, represented one third of the patterns. Sublineage L4.6.2 was mostly identified in the north of the country whereas L5 was mostly observed in the south and L6 was scattered. The spatial distribution of genotypes had genetic geographic gradients. We did not obtain results enabling the detection of drug-resistance mutations.

CONCLUSION/SIGNIFICANCE

We present the first national snapshot of the M. tuberculosis spoligotypes circulating in Nigeria based on ZN slides. Spoligotyping data can be obtained in a rapid and high-throughput manner with DNA extracted from ZN-stained slides, which may potentially improve our understanding of the genetic epidemiology of TB.

Mycobacterium tuberculosis Whole Genome Sequences From Southern India Suggest Novel Resistance Mechanisms and the Need for Region-Specific Diagnostics

Background.

India is home to 25% of all tuberculosis cases and the second highest number of multidrug resistant cases worldwide. However, little is known about the genetic diversity and resistance determinants of Indian Mycobacterium tuberculosis, particularly for the primary lineages found in India, lineages 1 and 3.

Methods.

We whole genome sequenced 223 randomly selected M. tuberculosis strains from 196 patients within the Tiruvallur and Madurai districts of Tamil Nadu in Southern India. Using comparative genomics, we examined genetic diversity, transmission patterns, and evolution of resistance.

Results.

Genomic analyses revealed (11) prevalence of strains from lineages 1 and 3, (11) recent transmission of strains among patients from the same treatment centers, (11) emergence of drug resistance within patients over time, (11) resistance gained in an order typical of strains from different lineages and geographies, (11) underperformance of known resistance-conferring mutations to explain phenotypic resistance in Indian strains relative to studies focused on other geographies, and (11) the possibility that resistance arose through mutations not previously implicated in resistance, or through infections with multiple strains that confound genotype-based prediction of resistance.

Conclusions.

In addition to substantially expanding the genomic perspectives of lineages 1 and 3, sequencing and analysis of M. tuberculosis whole genomes from Southern India highlight challenges of infection control and rapid diagnosis of resistant tuberculosis using current technologies. Further studies are needed to fully explore the complement of diversity and resistance determinants within endemic M. tuberculosis populations.

Mycobacterium tuberculosis Transmission among Elderly Persons, Yamagata Prefecture, Japan, 2009-2015

In many countries with low to moderate tuberculosis (TB) incidence, cases have shifted to elderly persons. It is unclear, however, whether these cases are associated with recent Mycobacterium tuberculosis transmission or represent reactivation of past disease. During 2009–2015, we performed a population-based TB investigation in Yamagata Prefecture, Japan, using in-depth contact tracing and 24-loci variable-number tandem-repeat typing optimized for Beijing family M. tuberculosis strains. We analyzed 494 strains, of which 387 (78.3%) were derived from elderly patients. Recent transmission with an epidemiologic link was confirmed in 22 clusters (70 cases). In 17 (77.3%) clusters, the source patient was elderly; 11 (64.7%) of the 17 clusters occurred in a hospital or nursing home. In this setting, the increase in TB cases was associated with M. tuberculosis transmissions from elderly persons. Prevention of transmission in places where elderly persons gather will be an effective strategy for decreasing TB incidence among predominantly elderly populations.

New Mycobacterium tuberculosis Complex Sublineage, Brazzaville, Congo

Tuberculosis is a leading cause of illness and death in Congo. No data are available about the population structure and transmission dynamics of the Mycobacterium tuberculosis complex strains prevalent in this central Africa country. On the basis of single-nucleotide polymorphisms detected by whole-genome sequencing, we phylogenetically characterized 74 MTBC isolates from Brazzaville, the capital of Congo. The diversity of the study population was high; most strains belonged to the Euro-American lineage, which split into Latin American Mediterranean, Uganda I, Uganda II, Haarlem, X type, and a new dominant sublineage named Congo type (n = 26). Thirty strains were grouped in 5 clusters (each within 12 single-nucleotide polymorphisms), from which 23 belonged to the Congo type. High cluster rates and low genomic diversity indicate recent emergence and transmission of the Congo type, a new Euro-American sublineage of MTBC.

What Is Resistance? Impact of Phenotypic versus Molecular Drug Resistance Testing on Therapy for Multi- and Extensively Drug-Resistant Tuberculosis

Rapid and accurate drug susceptibility testing (DST) is essential for the treatment of multi- and extensively drug-resistant tuberculosis (M/XDR-TB). We compared the utility of genotypic DST assays with phenotypic DST (pDST) using Bactec 960 MGIT or Löwenstein-Jensen to construct M/XDR-TB treatment regimens for a cohort of 25 consecutive M/XDR-TB patients and 15 possible anti-TB drugs. Genotypic DST results from Cepheid GeneXpert MTB/RIF (Xpert) and line probe assays (LPAs; Hain GenoType MTBDRplus 2.0 and MTBDRsl 2.0) and whole-genome sequencing (WGS) were translated into individual algorithm-derived treatment regimens for each patient. We further analyzed if discrepancies between the various methods were due to flaws in the genotypic or phenotypic test using MIC results. Compared with pDST, the average agreement in the number of drugs prescribed in genotypic regimens ranged from just 49% (95% confidence interval [CI], 39 to 59%) for Xpert and 63% (95% CI, 56 to 70%) for LPAs to 93% (95% CI, 88 to 98%) for WGS. Only the WGS regimens did not contain any drugs to which pDST showed resistance. Importantly, MIC testing revealed that pDST likely underestimated the true rate of resistance for key drugs (rifampin, levofloxacin, moxifloxacin, and kanamycin) because critical concentrations (CCs) were too high. WGS can be used to rule in resistance even in M/XDR strains with complex resistance patterns, but pDST for some drugs is still needed to confirm susceptibility and construct the final regimens. Some CCs for pDST need to be reexamined to avoid systematic false-susceptible results in low-level resistant isolates.

Drug-Resistant tuberculosis - primary transmission and management

The DOTS strategy assisted global tuberculosis (TB) control, but was unable to prevent the emergence and spread of drug-resistant strains. Genomic evidence confirms the transmission of drug-resistant Mycobacterium tuberculosis strains in many different settings, indicative of epidemic spread. These findings emphasise the need for enhanced infection control measures in health care and congregate settings. Young children in TB endemic areas are particularly vulnerable. Although advances in TB drug and vaccine development are urgently needed, improved access to currently available preventive therapy and treatment for drug resistant TB could reduce the disease burden and adverse outcomes experienced by children. We review new insights into the transmission dynamics of drug resistant TB, the estimated disease burden in children and optimal management strategies to consider.

Antimicrobial resistance in Mycobacterium tuberculosis: mechanistic and evolutionary perspectives

Antibiotic-resistant Mycobacterium tuberculosis strains are threatening progress in containing the global tuberculosis epidemic. Mycobacterium tuberculosis is intrinsically resistant to many antibiotics, limiting the number of compounds available for treatment. This intrinsic resistance is due to a number of mechanisms including a thick, waxy, hydrophobic cell envelope and the presence of drug degrading and modifying enzymes. Resistance to the drugs which are active against M. tuberculosis is, in the absence of horizontally transferred resistance determinants, conferred by chromosomal mutations. These chromosomal mutations may confer drug resistance via modification or overexpression of the drug target, as well as by prevention of prodrug activation. Drug resistance mutations may have pleiotropic effects leading to a reduction in the bacterium's fitness, quantifiable e.g. by a reduction in the in vitro growth rate. Secondary so-called compensatory mutations, not involved in conferring resistance, can ameliorate the fitness cost by interacting epistatically with the resistance mutation. Although the genetic diversity of M. tuberculosis is low compared to other pathogenic bacteria, the strain genetic background has been demonstrated to influence multiple aspects in the evolution of drug resistance. The rate of resistance evolution and the fitness costs of drug resistance mutations may vary as a function of the genetic background.

Comparative Whole-Genomic Analysis of an Ancient L2 Lineage Mycobacterium tuberculosis Reveals a Novel Phylogenetic Clade and Common Genetic Determinants of Hypervirulent Strains

Background: Development of improved therapeutics against tuberculosis (TB) is hindered by an inadequate understanding of the relationship between disease severity and genetic diversity of its causative agent, Mycobacterium tuberculosis. We previously isolated a hypervirulent M. tuberculosis strain H112 from an HIV-negative patient with an aggressive disease progression from pulmonary TB to tuberculous meningitis—the most severe manifestation of tuberculosis. Human macrophage challenge experiment demonstrated that the strain H112 exhibited significantly better intracellular survivability and induced lower level of TNF-α than the reference virulent strain H37Rv and other 123 clinical isolates.

Aim: The present study aimed to identify the potential genetic determinants of mycobacterial virulence that were common to strain H112 and hypervirulent M. tuberculosis strains of the same phylogenetic clade isolated in other global regions.

Methods: A low-virulent M. tuberculosis strain H54 which belonged to the same phylogenetic lineage (L2) as strain H112 was selected from a collection of 115 clinical isolates. Both H112 and H54 were whole-genome-sequenced using PacBio sequencing technology. A comparative genomics approach was adopted to identify mutations present in strain H112 but absent in strain H54. Subsequently, an extensive phylogenetic analysis was conducted by including all publically available M. tuberculosis genomes. Single-nucleotide-polymorphisms (SNPs) and structural variations (SVs) common to hypervirulent strains in the global collection of genomes were considered as potential genetic determinants of hypervirulence.

Results:Sequencing data revealed that both H112 and H54 were identified as members of the same sub-lineage L2.2.1. After excluding the lineage-related mutations shared between H112 and H54, we analyzed the phylogenetic relatedness of H112 with global collection of M. tuberculosis genomes (n = 4,338), and identified a novel phylogenetic clade in which four hypervirulent strains isolated from geographically diverse regions were clustered together. All hypervirulent strains in the clade shared 12 SNPs and 5 SVs with H112, including those affecting key virulence-associated loci, notably, a deleterious SNP (rv0178 p. D150E) within mce1 operon and an intergenic deletion (854259_ 854261delCC) in close-proximity to phoP.

Conclusion: The present study identified common genetic factors in a novel phylogenetic clade of hypervirulent M. tuberculosis. The causative role of these mutations in mycobacterial virulence should be validated in future study.

Tracing the origin of Treponema pallidum in China using next-generation sequencing

Syphilis is a systemic sexually transmitted disease caused by Treponema pallidum ssp. pallidum (TPA). The origin and genetic background of Chinese TPA strains remain unclear. We identified a total of 329 single-nucleotide variants (SNVs) in eight Chinese TPA strains using next-generation sequencing. All of the TPA strains were clustered into three lineages, and Chinese TPA strains were grouped in Lineage 2 based on phylogenetic analysis. The phylogeographical data showed that TPA strains originated earlier than did T. pallidum ssp. pertenue (TPE) and T. pallidum ssp. endemicum (TPN) strains and that Chinese TPA strains might be derived from recombination between Lineage 1 and Lineage 3. Moreover, we found through a homology modeling analysis that a nonsynonymous substitution (I415F) in the PBP3 protein might affect the structural flexibility of PBP3 and the binding constant for substrates based on its possible association with penicillin resistance in T. pallidum. Our findings provide new insight into the molecular foundation of the evolutionary origin of TPA and support the development of novel diagnostic/therapeutic technology for syphilis.

Multi-clonal evolution of multi-drug-resistant/extensively drug-resistant Mycobacterium tuberculosis in a high-prevalence setting of Papua New Guinea for over three decades

An outbreak of multi-drug resistant (MDR) tuberculosis (TB) has been reported on Daru Island, Papua New Guinea. Mycobacterium tuberculosis strains driving this outbreak and the temporal accrual of drug resistance mutations have not been described. Whole genome sequencing of 100 of 165 clinical isolates referred from Daru General Hospital to the Supranational reference laboratory, Brisbane, during 2012–2015 revealed that 95 belonged to a single modern Beijing sub-lineage strain. Molecular dating suggested acquisition of streptomycin and isoniazid resistance in the 1960s, with potentially enhanced virulence mediated by an mycP1 mutation. The Beijing sub-lineage strain demonstrated a high degree of co-resistance between isoniazid and ethionamide (80/95; 84.2 %) attributed to an inhA promoter mutation combined with inhA and ndh coding mutations. Multi-drug resistance, observed in 78/95 samples, emerged with the acquisition of a typical rpoB mutation together with a compensatory rpoC mutation in the 1980s. There was independent acquisition of fluoroquinolone and aminoglycoside resistance, and evidence of local transmission of extensively drug resistant (XDR) strains from 2009. These findings underline the importance of whole genome sequencing in informing an effective public health response to MDR/XDR TB.