Evidence for Host-Bacterial Co-evolution via Genome Sequence Analysis of 480 Thai Mycobacterium tuberculosis Lineage 1 Isolates

The first insight into Mycobacterium tuberculosis complex isolates in the lower northern region in Thailand

BACKGROUND

Tuberculosis (TB) remains an important infectious disease and different genotypes have been reported. This study aimed to investigate the genetic diversity and molecular epidemiology of TB in the lower northern region of Thailand, where genotyping data are limited.

METHODS

A total of 159 Mycobacterium tuberculosis complex (MTBC) isolates from this region were genotyped by spoligotyping and the major spoligotypes were further subdivided by the mycobacterial interspersed repetitive unit-variable number tandem repeat (MIRU-VNTR) method.

RESULTS

Spoligotyping identified 34 types and classified them into 14 clusters. East African-Indian (EAI) groups were the most frequent (44.7%), followed by Beijing (36.5%), with a higher prevalence of drug resistance. By 15-loci MIRU-VNTR typing, the major groups of the Beijing and EAI2_NTB were further differentiated into 44 and 21 subtypes forming 9 and 5 subclusters with cluster rates of 0.26 and 0.44, respectively. The Hunter-Gaston Discriminatory Index among the Beijing and EAI2_NTB groups were 0.987 and 0.931, respectively, indicating high diversity.

CONCLUSIONS

This is the first look at the MTBC genotypes in the lower northern region of Thailand, which could aid in understanding the distribution and potential spread of MTBC and Mycobacterium bovis in the target region to support TB control in Thailand.

Towards the reconstruction of a global TB history using a new pipeline "TB-Annotator"

Mycobacterium tuberculosis complex (MTBC) has a population structure consisting of 9 human and animal lineages. The genomic diversity within these lineages is a pathogenesis factor that affects virulence, transmissibility, host response, and antibiotic resistance. Hence it is important to develop improved information systems for tracking and understanding the spreading and evolution of genomes. We present results obtained thanks to a new informatics platform for computational biology of MTBC, that uses a convenience sample from public/private SRAs, designated as TB-Annotator. Version 1 was a first interactive graphic-based web tool based on 15,901 representative genomes. Version 2, still interactive, is a more sophisticated database, developed using the Snakemake Workflow Management System (WMS) that allows an unsupervised global and scalable analysis of the content of the USA National Center for Biotechnology Information Short Read Archives database. This platform analyzes nucleotide variants, the presence/absence of genes, known regions of difference and detect new deletions, the insertion sites of mobile genetic elements, and allows phylogenetic trees to be built, imported in a graphical interface and interactively analyzed between the data and the tree. The objective of TB-Annotator is triple: detect recent epidemiological links, reconstruct distant phylogeographical histories as well as perform more complex phenotypic/genotypic Genome-Wide Association Studies (GWAS). In this paper, we compare the various taxonomic SNPs-based labels and hierarchies previously described in recent reference papers for L1, and present a comparative analysis that allows identification of alias and thus provides the basis of a future unifying naming scheme for L1 sublineages. We present a global phylogenetic tree built with RAxML-NG, and one on L2; at the time of writing, we characterized about 200 sublineages, with many new ones; a detail tree for Modern L2 and a hierarchical scheme allowing to facilitate L2 lineage assignment are also presented.

A revised SNP-based barcoding scheme for typing Mycobacterium tuberculosis complex isolates

The development of whole-genome sequencing technologies is gradually leading to a more detailed description of the population structure of the Mycobacterium tuberculosis complex (MTBC). In this study, we correlated previously published classifications on a collection of more than 10,000 genomes and proposed a new, comprehensive nomenclature that unifies the existing ones. In total, we identified 169 lineages and sublineages of M. tuberculosis/M. africanum and 9 animal-adapted species. For the purpose of organizing these genotypes in a more streamlined manner, we stratified them into five hierarchical levels. To represent the classification and compare it with the reference, we compiled a confirmatory data set of 670 high-quality isolates, which includes all genotypes and species of MTBC, and this confirmatory data set can serve as a basis for further studies. We proposed a set of 213 robust barcoding single-nucleotide polymorphisms and a suitable workflow for reliable differentiation of genotypes and species within the complex. This work integrates the results of all the major systematized studies to date to provide an understanding of the global diversity of the MTBC population structure. The results of this work may ultimately help to reliably determine the pathogen genotype and associate it with traits that reflect its prevalence, virulence, vaccination, and treatment efficiency, as well as to reliably find natural features revealed during its spread. IMPORTANCE Through years of research into the Mycobacterium tuberculosis complex (MTBC), a number of ambiguous phylogenetic classifications have emerged, which often overlap with one another. In the present study, we have combined all major studies on MTBC classification and inferred a unified, most complete to date classification and accompanying SNP barcodes.

Analysis of whiB7 in Mycobacterium tuberculosis reveals novel AT-hook deletion mutations

Mutations in whiB7 have been associated with both hypersusceptibility and resistance to various antibiotics in Mycobacterium tuberculosis (Mtb). Unlocking the secrets of antibiotic resistance in the bacterium, we examined mutations in the coding sequences of whiB7 of over 40,000 diverse Mtb isolates. Our results unveil the dominant c.191delG (Gly64delG) mutation, present in all members of the lineage L1.2.2 and its impact on WhiB7's conserved GVWGG-motif, causing conformational changes and deletion of the C-terminal AT-hook. Excitingly, we discovered six unique mutations associated with partial or total deletion of the AT-hook, specific to certain sublineages. Our findings suggest the selective pressures driving these mutations, underlining the potential of genomics to advance our understanding of Mtb's antibiotic resistance. As tuberculosis remains a global health threat, our study offers valuable insights into the diverse nature and functional consequences of whiB7 mutations, paving the way for the development of novel therapeutic interventions.

Comparison of in silico predicted Mycobacterium tuberculosis spoligotypes and lineages from whole genome sequencing data

Bacterial strain-types in the Mycobacterium tuberculosis complex underlie tuberculosis disease, and have been associated with drug resistance, transmissibility, virulence, and host-pathogen interactions. Spoligotyping was developed as a molecular genotyping technique used to determine strain-types, though recent advances in whole genome sequencing (WGS) technology have led to their characterization using SNP-based sub-lineage nomenclature. Notwithstanding, spoligotyping remains an important tool and there is a need to study the congruence between spoligotyping-based and SNP-based sub-lineage assignation. To achieve this, an in silico spoligotype prediction method ("Spolpred2") was developed and integrated into TB-Profiler. Lineage and spoligotype predictions were generated for > 28 k isolates and the overlap between strain-types was characterized. Major spoligotype families detected were Beijing (25.6%), T (18.6%), LAM (13.1%), CAS (9.4%), and EAI (8.3%), and these broadly followed known geographic distributions. Most spoligotypes were perfectly correlated with the main MTBC lineages (L1-L7, plus animal). Conversely, at lower levels of the sub-lineage system, the relationship breaks down, with only 65% of spoligotypes being perfectly associated with a sub-lineage at the second or subsequent levels of the hierarchy. Our work supports the use of spoligotyping (membrane or WGS-based) for low-resolution surveillance, and WGS or SNP-based systems for higher-resolution studies.

Nanoscale clustering of mycobacterial ligands and DC-SIGN host receptors are key determinants for pathogen recognition

The bacterial pathogen Mycobacterium tuberculosis binds to the C-type lectin DC-SIGN (dendritic cell-specific intercellular adhesion molecule 3-grabbing nonintegrin) on dendritic cells to evade the immune system. While DC-SIGN glycoconjugate ligands are ubiquitous among mycobacterial species, the receptor selectively binds pathogenic species from the M. tuberculosis complex (MTBC). Here, we unravel the molecular mechanism behind this intriguing selective recognition by means of a multidisciplinary approach combining single-molecule atomic force microscopy with Förster resonance energy transfer and bioassays. Molecular recognition imaging of mycobacteria demonstrates that the distribution of DC-SIGN ligands markedly differs between Mycobacterium bovis Bacille Calmette-Guérin (BCG) (model MTBC species) and Mycobacterium smegmatis (non-MTBC species), the ligands being concentrated into dense nanodomains on M. bovis BCG. Upon bacteria-host cell adhesion, ligand nanodomains induce the recruitment and clustering of DC-SIGN. Our study highlights the key role of clustering of both ligands on MTBC species and DC-SIGN host receptors in pathogen recognition, a mechanism that might be widespread in host-pathogen interactions.

Genomic Sequencing Profiles of Mycobacterium tuberculosis in Mandalay Region, Myanmar

This study aimed to characterize whole-genome sequencing (WGS) information of Mycobacterium tuberculosis (Mtb) in the Mandalay region of Myanmar. It was a cross-sectional study conducted with 151 Mtb isolates obtained from the fourth nationwide anti-tuberculosis (TB) drug-resistance survey. Frequency of lineages 1, 2, 3, and 4 were 55, 65, 9, and 22, respectively. The most common sublineage was L1.1.3.1 (n = 31). Respective multi-drug resistant tuberculosis (MDR-TB) frequencies were 1, 1, 0, and 0. Four clusters of 3 (L2), 2 (L4), 2 (L1), and 2 (L2) isolates defined by a 20-single-nucleotide variant (SNV) cutoff were detected. Simpson's index for sublineages was 0.0709. Such high diversity suggests that the area probably had imported Mtb from many geographical sources. Relatively few genetic clusters and MDR-TB suggest there is a chance the future control will succeed if it is carried out properly.

Distribution and origins of Mycobacterium tuberculosis L4 in Southeast Asia

Molecular and genomic studies have revealed that Mycobacterium tuberculosis Lineage 4 (L4, Euro-American lineage) emerged in Europe before becoming distributed around the globe by trade routes, colonial migration and other historical connections. Although L4 accounts for tens or hundreds of thousands of tuberculosis (TB) cases in multiple Southeast Asian countries, phylogeographical studies have either focused on a single country or just included Southeast Asia as part of a global analysis. Therefore, we interrogated public genomic data to investigate the historical patterns underlying the distribution of L4 in Southeast Asia and surrounding countries. We downloaded 6037 genomes associated with 29 published studies, focusing on global analyses of L4 and Asian studies of M. tuberculosis. We identified 2256 L4 genomes including 968 from Asia. We show that 81 % of L4 in Thailand, 51 % of L4 in Vietnam and 9 % of L4 in Indonesia belong to sub-lineages of L4 that are rarely seen outside East and Southeast Asia (L4.2.2, L4.4.2 and L4.5). These sub-lineages have spread between East and Southeast Asian countries, with no recent European ancestor. Although there is considerable uncertainty about the exact direction and order of intra-Asian M. tuberculosis dispersal, due to differing sampling frames between countries, our analysis suggests that China may be the intermediate location between Europe and Southeast Asia for two of the three predominantly East and Southeast Asian L4 sub-lineages (L4.2.2 and L4.5). This new perspective on L4 in Southeast Asia raises the possibility of investigating host population-specific evolution and highlights the need for more structured sampling from Southeast Asian countries to provide more certainty of the historical and current routes of dispersal.

Investigating the Diversity of Tuberculosis Spoligotypes with Dimensionality Reduction and Graph Theory

The spoligotype is a graphical description of the CRISPR locus present in Mycobacterium tuberculosis, which has the particularity of having only 68 possible spacers. This spoligotype, which can be easily obtained either in vitro or in silico, allows to have a summary information of lineage or even antibiotic resistance (when known to be associated to a particular cluster) at a lower cost. The objective of this article is to show that this representation is richer than it seems, and that it is under-exploited until now. We first recall an original way to represent these spoligotypes as points in the plane, allowing to highlight possible sub-lineages, particularities in the animal strains, etc. This graphical representation shows clusters and a skeleton in the form of a graph, which led us to see these spoligotypes as vertices of an unconnected directed graph. In this paper, we therefore propose to exploit in detail the description of the variety of spoligotypes using a graph, and we show to what extent such a description can be informative.

Identification of region of difference and H37Rv-related deletion in Mycobacterium tuberculosis complex by structural variant detection and genome assembly

Mycobacterium tuberculosis complex (MTBC), the main cause of TB in humans and animals, is an extreme example of genetic homogeneity, whereas it is still nevertheless separated into various lineages by numerous typing methods, which differ in phenotype, virulence, geographic distribution, and host preference. The large sequence polymorphism (LSP), incorporating region of difference (RD) and H37Rv-related deletion (RvD), is considered to be a powerful means of constructing phylogenetic relationships within MTBC. Although there have been many studies on LSP already, focusing on the distribution of RDs in MTBC and their impact on MTB phenotypes, a crumb of new lineages or sub-lineages have been excluded and RvDs have received less attention. We, therefore, sampled a dataset of 1,495 strains, containing 113 lineages from the laboratory collection, to screen for RDs and RvDs by structural variant detection and genome assembly, and examined the distribution of RvDs in MTBC, including RvD2, RvD5, and cobF region. Consistent with genealogical delineation by single nucleotide polymorphism (SNP), we identified 125 RDs and 5 RvDs at the species, lineage, or sub-lineage levels. The specificities of RDs and RvDs were further investigated in the remaining 10,218 strains, suggesting that most of them were highly specific to distinct phylogenetic groups, could be used as stable genetic markers in genotyping. More importantly, we identified 34 new lineage or evolutionary branch specific RDs and 2 RvDs, also demonstrated the distribution of known RDs and RvDs in MTBC. This study provides novel details about deletion events that have occurred in distinct phylogenetic groups and may help to understand the genealogical differentiation.

Distribution and identification of Mycobacterium tuberculosis lineage in Kashgar prefecture

Objectives

Kashgar prefecture is an important transportation and trade hub with a high incidence of tuberculosis. The following study analyzed the composition and differences in Mycobacterium tuberculosis (M.tb) lineage and specific tags to distinguish the lineage of the M.tb in Kashgar prefecture, thus providing a basis for the classification and diagnosis of tuberculosis in this area.

Methods

Whole-genome sequencing (WGS) of 161 M.tb clinical strains was performed. The phylogenetic tree was constructed using Maximum Likelihood (ML) based on single nucleotide polymorphisms (SNPs) and verified through principal component analysis (PCA). The composition structure of M.tb in different regions was analyzed by combining geographic information.

Results

M.tb clinical strains were composed of lineage 2 (73/161, 45.34%), lineage 3 (52/161, 32.30%) and lineage 4 (36/161, 22.36%). Moreover, the 3 lineages were subdivided into 11 sublineages, among which lineage 2 included lineage 2.2.2/Asia Ancestral 1 (9/73, 12.33%), lineage 2.2.1-Asia Ancestral 2 (9/73, 12.33%), lineage 2.2.1-Asia Ancestral 3 (18/73, 24.66%), and lineage 2.2.1-Modern Beijing (39/73, 53.42%). Lineage 3 included lineage 3.2 (14/52, 26.92%) and lineage 3.3 (38/52, 73.08%), while lineage 4 included lineage 4.1 (3/36, 8.33%), lineage 4.2 (2/36, 5.66%), lineage 4.4.2 (1/36, 2.78%), lineage 4.5 (28/36, 77.78%) and lineage 4.8 (2/36, 5.66%), all of which were consistent with the PCA results. One hundred thirty-six markers were proposed for discriminating known circulating strains. Reconstruction of a phylogenetic tree using the 136 SNPs resulted in a tree with the same number of delineated clades. Based on geographical location analysis, the composition of Lineage 2 in Kashgar prefecture (45.34%) was lower compared to other regions in China (54.35%-90.27%), while the composition of Lineage 3 (32.30%) was much higher than in other regions of China (0.92%-2.01%), but lower compared to the bordering Pakistan (70.40%).

Conclusion

Three lineages were identified in M.tb clinical strains from Kashgar prefecture, with 136 branch-specific SNP. Kashgar borders with countries that have a high incidence of tuberculosis, such as Pakistan and India, which results in a large difference between the M.tb lineage and sublineage distribution in this region and other provinces of China.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12879-022-07307-4.

A review on the application of bioinformatics tools in food microbiome studies

There is currently a transformed interest toward understanding the impact of fermentation on functional food development due to growing consumer interest on modified health benefits of sustainable foods. In this review, we attempt to summarize recent findings regarding the impact of Next-generation sequencing and other bioinformatics methods in the food microbiome and use prediction software to understand the critical role of microbes in producing fermented foods. Traditionally, fermentation methods and starter culture development were considered conventional methods needing optimization to eliminate errors in technique and were influenced by technical knowledge of fermentation. Recent advances in high-output omics innovations permit the implementation of additional logical tactics for developing fermentation methods. Further, the review describes the multiple functions of the predictions based on docking studies and the correlation of genomic and metabolomic analysis to develop trends to understand the potential food microbiome interactions and associated products to become a part of a healthy diet.

Whole-genome single nucleotide variant phylogenetic analysis of Mycobacterium tuberculosis Lineage 1 in endemic regions of Asia and Africa

Mycobacterium tuberculosis (Mtb) lineage 1 (L1) contributes considerably to the disease morbidity. While whole genome sequencing (WGS) is increasingly used for studying Mtb, our understanding of genetic diversity of L1 remains limited. Using phylogenetic analysis of WGS data from endemic range in Asia and Africa, we provide an improved genotyping scheme for L1. Mapping deletion patterns of the 68 direct variable repeats (DVRs) in the CRISPR region of the genome onto the phylogeny provided supporting evidence that the CRISPR region evolves primarily by deletion, and hinted at a possible Southeast Asian origin of L1. Both phylogeny and DVR patterns clarified some relationships between different spoligotypes, and highlighted the limited resolution of spoligotyping. We identified a diverse repertoire of drug resistance mutations. Altogether, this study demonstrates the usefulness of WGS data for understanding the genetic diversity of L1, with implications for public health surveillance and TB control. It also highlights the need for more WGS studies in high-burden but underexplored regions.

Connection between two historical tuberculosis outbreak sites in Japan, Honshu, by a new ancestral Mycobacterium tuberculosis L2 sublineage

By gathering 680 publicly available Sequence Read Archives from isolates of Mycobacterium tuberculosis complex (MTBC) including 190 belonging to the lineage 2 Beijing , and using an in-house bioinformatical pipeline, the TB-Annotator , that analyses more than 50 000 characters, we describe herein a new L2 sublineage from 20 isolates found in the Tochigi province, (Japan), that we designate as asia ancestral 5 (AAnc5). These isolates harbour a number of specific criteria (42 SNPs) and their intra-cluster pairwise distance suggests historical and not epidemiological transmission. These isolates harbour a mutation in rpoC , and do not fulfil, any of the modern Beijing lineage criteria, nor any of the other ancestral Beijing lineages described so far. Asia ancestral 5 isolates do not possess mutT2 58 and ogt 12 characteristics of modern Beijing , but possess ancestral Beijing SNPs characteristics. By looking into the literature, we found a reference isolate ID381, described in Kobe and Osaka belonging to the ‘G3’ group, sharing 36 out of the 42 specific SNPs found in AAnc5. We also assessed the intermediate position of the asia ancestral 4 (AAnc4) sublineage recently described in Thailand and propose an improved classification of the L2 that now includes AAnc4 and AAnc5. By increasing the recruitment into TB-Annotator to around 3000 genomes (including 642 belonging to L2), we confirmed our results and discovered additional historical ancestral L2 branches that remain to be investigated in more detail. We also present, in addition, some anthropological and historical data from Chinese and Japan history of tuberculosis, as well as from Korea, that could support our results on L2 evolution. This study shows that the reconstruction of the early history of tuberculosis in Asia is likely to reveal complex patterns since its emergence.

A Comprehensive Map of Mycobacterium tuberculosis Complex Regions of Difference

Mycobacterium tuberculosis complex (MTBC) species are classic examples of genetically monomorphic microorganisms due to their low genetic variability. Whole-genome sequencing made it possible to describe both the main species within the complex and M. tuberculosis lineages and sublineages. This differentiation is based on single nucleotide polymorphisms (SNPs) and large sequence polymorphisms in the so-called regions of difference (RDs). Although a number of studies have been performed to elucidate RD localizations, their distribution among MTBC species, and their role in the bacterial life cycle, there are some inconsistencies and ambiguities in the localization of RDs in different members of the complex. To address this issue, we conducted a thorough search for all possible deletions in the WGS data collection comprising 721 samples representing the full MTBC diversity. Discovered deletions were compared with a list of all previously described RDs. As with the SNP-based analysis, we confirmed the specificities of 79 regions at the species, lineage, or sublineage level, 17 of which are described for the first time. We also present RDscan (https://github.com/dbespiatykh/RDscan), an open-source workflow, which detects deletions from short-read sequencing data and correlates the results with high-specificity RDs, curated in this study. Testing of the workflow on a collection comprising ∼7,000 samples showed a high specificity of the found RDs. This study provides novel details that can contribute to a better understanding of the species differentiation within the MTBC and can help to determine how individual clusters evolve within various MTBC species.

IMPORTANCE Reductive genome evolution is one of the most important and intriguing adaptation strategies of different living organisms to their environment. Mycobacterium offers several notorious examples of either naturally reduced (Mycobacterium leprae) or laboratory-reduced (Mycobacterium bovis BCG) genomes. Mycobacterium tuberculosis complex has its phylogeny unambiguously framed by large sequence polymorphisms that present unidirectional unique event changes. In the present study, we curated all known regions of difference and analyzed both Mycobacterium tuberculosis and animal-adapted MTBC species. For 79 loci, we have shown a relationship with phylogenetic units, which can serve as a marker for diagnosing or studying biological effects. Moreover, intersections were found for some loci, which may indicate the nonrandomness of these processes and the involvement of these regions in the adaptation of bacteria to external conditions.

Clusters of Drug-Resistant Mycobacterium tuberculosis Detected by Whole-Genome Sequence Analysis of Nationwide Sample, Thailand, 2014-2017

Multidrug-resistant tuberculosis (MDR TB), pre-extensively drug-resistant tuberculosis (pre-XDR TB), and extensively drug-resistant tuberculosis (XDR TB) complicate disease control. We analyzed whole-genome sequence data for 579 phenotypically drug-resistant M. tuberculosis isolates (28% of available MDR/pre-XDR and all culturable XDR TB isolates collected in Thailand during 2014–2017). Most isolates were from lineage 2 (n = 482; 83.2%). Cluster analysis revealed that 281/579 isolates (48.5%) formed 89 clusters, including 205 MDR TB, 46 pre-XDR TB, 19 XDR TB, and 11 poly–drug-resistant TB isolates based on genotypic drug resistance. Members of most clusters had the same subset of drug resistance-associated mutations, supporting potential primary resistance in MDR TB (n = 176/205; 85.9%), pre-XDR TB (n = 29/46; 63.0%), and XDR TB (n = 14/19; 73.7%). Thirteen major clades were significantly associated with geography (p<0.001). Clusters of clonal origin contribute greatly to the high prevalence of drug-resistant TB in Thailand.

Phenotypic and genotypic features of the Mycobacterium tuberculosis lineage 1 subgroup in central Vietnam

Mycobacterium tuberculosis (Mtb) has different features depending on different geographic areas. We collected Mtb strains from patients with smear-positive pulmonary tuberculosis in Da Nang, central Vietnam. Using a whole genome sequencing platform, including genome assembly complemented by long-read-sequencing data, genomic characteristics were studied. Of 181 Mtb isolates, predominant Vietnamese EAI4_VNM and EAI4-like spoligotypes (31.5%), ZERO strains (5.0%), and part of EAI5 (11.1%) were included in a lineage-1 (L1) sublineage, i.e., L1.1.1.1. These strains were found less often in younger people, and they genetically clustered less frequently than other modern strains. Patients infected with ZERO strains demonstrated less lung infiltration. A region in RD2bcg spanning six loci, i.e., PE_PGRS35, cfp21, Rv1985c, Rv1986, Rv1987, and erm(37), was deleted in EAI4_VNM, EAI4-like, and ZERO strains, whereas another 118 bp deletion in furA was specific only to ZERO strains. L1.1.1.1-sublineage-specific deletions in PE_PGRS4 and PE_PGRS22 were also identified. RD900, seen in ancestral lineages, was present in majority of the L1 members. All strains without IS6110 (5.0%) had the ZERO spoligo-pattern. Distinctive features of the ancestral L1 strains provide a basis for investigation of the modern versus ancestral Mtb lineages and allow consideration of countermeasures against this heterogeneous pathogen.

Genomic epidemiology of tuberculosis in eastern Malaysia: insights for strengthening public health responses

Tuberculosis is a leading public health priority in eastern Malaysia. Knowledge of the genomic epidemiology of tuberculosis can help tailor public health interventions. Our aims were to determine tuberculosis genomic epidemiology and characterize resistance mutations in the ethnically diverse city of Kota Kinabalu, Sabah, located at the nexus of Malaysia, Indonesia, Philippines and Brunei. We used an archive of prospectively collected Mycobacterium tuberculosis samples paired with epidemiological data. We collected sputum and demographic data from consecutive consenting outpatients with pulmonary tuberculosis at the largest tuberculosis clinic from 2012 to 2014, and selected samples from tuberculosis inpatients from the tertiary referral centre during 2012-2014 and 2016-2017. Two hundred and eight M. tuberculosis sequences were available for analysis, representing 8 % of cases notified during the study periods. Whole-genome phylogenetic analysis demonstrated that most strains were lineage 1 (195/208, 93.8 %), with the remainder being lineages 2 (8/208, 3.8 %) or 4 (5/208, 2.4 %). Lineages or sub-lineages were not associated with patient ethnicity. The lineage 1 strains were diverse, with sub-lineage 1.2.1 being dominant (192, 98 %). Lineage 1.2.1.3 isolates were geographically most widely distributed. The greatest diversity occurred in a border town sub-district. The time to the most recent common ancestor for the three major lineage 1.2.1 clades was estimated to be the year 1966 (95 % HPD 1948-1976). An association was found between failure of culture conversion by week 8 of treatment and infection with lineage 2 (4/6, 67 %) compared with lineage 1 strains (4/83, 5 %) (P<0.001), supporting evidence of greater virulence of lineage 2 strains. Eleven potential transmission clusters (SNP difference ≤12) were identified; at least five included people living in different sub-districts. Some linked cases spanned the whole 4-year study period. One cluster involved a multidrug-resistant tuberculosis strain matching a drug-susceptible strain from 3 years earlier. Drug resistance mutations were uncommon, but revealed one phenotype-genotype mismatch in a genotypically multidrug-resistant isolate, and rare nonsense mutations within the katG gene in two isolates. Consistent with the regionally mobile population, M. tuberculosis strains in Kota Kinabalu were diverse, although several lineage 1 strains dominated and were locally well established. Transmission clusters - uncommonly identified, likely attributable to incomplete sampling - showed clustering occurring across the community, not confined to households or sub-districts. The findings indicate that public health priorities should include active case finding and early institution of tuberculosis management in mobile populations, while there is a need to upscale effective contact investigation beyond households to include other contacts within social networks.

Estimation of the global burden of Mycobacterium tuberculosis lineage 1

Tuberculosis is still problematic as it affects large numbers of people globally. Mycobacterium tuberculosis Lineage 1 (L1) or Indo Oceanic Lineage, one of widespread major lineages, has a specific geographic distribution and high mortality. It is highly diverse and endemic in several high burden countries. However, studies on the global burden of L1 and its sublineages remain limited. This may lead to the underestimation of the importance of its variance in developing and applying tuberculosis control measures. This study aimed to estimate the number of patients infected with M. tuberculosis L1 and its sublineages worldwide. The proportion of L1 among tuberculosis patients was searched in published reports from countries around the world and the number of patients was calculated based on a WHO report on country incidences and populations. The numbers of patients infected with the five major sublineages, namely L1.1.1, L1.1.2, L1.1.3, L1.2.1, and L1.2.2 were estimated where information was available. It was found that L1 accounted for 28% of global tuberculosis cases in 2012 and 2018. Over 80% of the L1 global burden was in India, the Philippines, Indonesia and Bangladesh, which are also among the countries with highest absolute numbers of tuberculosis patients in the world. Globally, the estimated number of patients infected with M. tuberculosis L1.2.1 and L1.1.2 was over 1.1 million and of patients infected with L1.1.1 was about 200,000. This study demonstrated that L1 contributes significantly to the global burden of tuberculosis. To achieve the End TB Strategy, more attention needs to be paid to the responses of M. tuberculosis L1 to various control measures.

CRISPRbuilder-TB: "CRISPR-builder for tuberculosis". Exhaustive reconstruction of the CRISPR locus in mycobacterium tuberculosis complex using SRA

Mycobacterium tuberculosis complex (MTC) CRISPR locus diversity has long been studied solely investigating the presence/absence of a known set of spacers. Unveiling the genetic mechanisms of its evolution requires a more exhaustive reconstruction in a large amount of representative strains.

In this article, we point out and resolve, with a new pipeline, the problem of CRISPR reconstruction based directly on short read sequences in M. tuberculosis. We first show that the process we set up, that we coin as “CRISPRbuilder-TB” (https://github.com/cguyeux/CRISPRbuilder-TB), allows an efficient reconstruction of simulated or real CRISPRs, even when including complex evolutionary steps like the insertions of mobile elements. Compared to more generalist tools, the whole process is much more precise and robust, and requires only minimal manual investigation. Second, we show that more than 1/3 of the currently complete genomes available for this complex in the public databases contain largely erroneous CRISPR loci. Third, we highlight how both the classical experimental in vitro approach and the basic in silico spoligotyping provided by existing analytic tools miss a whole diversity of this locus in MTC, by not capturing duplications, spacer and direct repeats variants, and IS6110 insertion locations. This description is extended in a second article that describes MTC-CRISPR diversity and suggests general rules for its evolution. This work opens perspectives for an in-depth exploration of M. tuberculosis CRISPR loci diversity and of mechanisms involved in its evolution and its functionality, as well as its adaptation to other CRISPR locus-harboring bacterial species.

In this article, we tackle the bioinformatical issue of the reconstruction of the Mycobacterium tuberculosis complex CRISPR locus using short read sequences without requiring genome assembly. We first show that many complete genomes, as found in public databases and often reconstructed by de novo assemblies, often contain errors on this locus as well as on other repeated sequences. We provide an in-depth description of our new method, designated as ‘CRISPRbuilder-TB’, and we show that our method provides much more exhaustive and reliable information (on DR variants, spacer diversity, global structure) than Crass and CRISPR_detector. The new and unsuspected genomic diversity we detected is described in a companion paper. Scripts are available to adapt the tool to other species.

Evolutionary Genetics of Mycobacterium tuberculosis and HIV-1: "The Tortoise and the Hare"

The already enormous burden caused by Mycobacterium tuberculosis and Human Immunodeficiency Virus type 1 (HIV-1) alone is aggravated by co-infection. Despite obvious differences in the rate of evolution comparing these two human pathogens, genetic diversity plays an important role in the success of both. The extreme evolutionary dynamics of HIV-1 is in the basis of a robust capacity to evade immune responses, to generate drug-resistance and to diversify the population-level reservoir of M group viral subtypes. Compared to HIV-1 and other retroviruses, M. tuberculosis generates minute levels of genetic diversity within the host. However, emerging whole-genome sequencing data show that the M. tuberculosis complex contains at least nine human-adapted phylogenetic lineages. This level of genetic diversity results in differences in M. tuberculosis interactions with the host immune system, virulence and drug resistance propensity. In co-infected individuals, HIV-1 and M. tuberculosis are likely to co-colonize host cells. However, the evolutionary impact of the interaction between the host, the slowly evolving M. tuberculosis bacteria and the HIV-1 viral "mutant cloud" is poorly understood. These evolutionary dynamics, at the cellular niche of monocytes/macrophages, are also discussed and proposed as a relevant future research topic in the context of single-cell sequencing.

Identification and in silico functional prediction of lineage-specific SNPs distributed in DosR-related proteins and resuscitation-promoting factor proteins of Mycobacterium tuberculosis

One-third of the world population is infected by Mycobacterium tuberculosis, which may persist in the latent or dormant state. Bacteria can shift to dormancy when encountering harsh conditions such as low oxygen, nutrient starvation, high acidity and host immune defenses. Genes related to the dormancy survival regulator (DosR) regulon are responsible for the inhibition of aerobic respiration and replication, which is required to enter dormancy. Conversely, resuscitation-promoting factor (rpf) proteins participate in reactivation from dormancy and the development of active tuberculosis (TB). Many DosR regulon and rpf proteins are immunodominant T cell antigens that are highly expressed in latent TB infection. They could serve as TB vaccine candidates and be used for diagnostic development. We explored the genetic polymorphisms of 50 DosR-related genes and 5 rpf genes among 1,170 previously sequenced clinical M. tuberculosis genomes. Forty-three lineage- or sublineage-specific nonsynonymous single nucleotide polymorphisms (nsSNPs) were identified. Ten nsSNPs were specific to all Mtb isolates belonging to lineage 1 (L1). Two common sublineages, the Beijing family (L2.2) and EAI2 (L1.2.1), differed at as many as 26 lineage- or sublineage-specific SNPs. DosR regulon genes related to membrane proteins and the rpf family possessed mean dN/dS ratios greater than one, suggesting that they are under positive selection. Although the T cell epitope regions of DosR-related and rpf antigens were quite conserved, we found that the epitopes in L1 had higher rates of genetic polymorphisms than the other lineages. Some mutations in immunogenic epitopes of the antigens were specific to particular M. tuberculosis lineages. Therefore, the genetic diversity of the DosR regulon and rpf proteins might impact the adaptation of M. tuberculosis to the dormant state and the immunogenicity of latency antigens, which warrants further investigation.

Genomic evidence supporting the clonal expansion of extensively drug-resistant tuberculosis bacteria belonging to a rare protoBeijing genotype

Tuberculosis disease (TB), caused by Mycobacterium tuberculosis, is a major public health issue in Thailand. The high prevalence of modern Beijing (Lineage 2.2.1) strains has been associated with multi- and extensively drug-resistant infections (MDR-, XDR-TB), complicating disease control. The impact of rarer proto-Beijing (L2.1) strains is less clear. In our study of thirty-seven L2.1 clinical isolates spanning thirteen years, we found a high prevalence of XDR-TB cases (32.4%). With ≤ 12 pairwise SNP distances, 43.2% of L2.1 patients belong to MDR-TB or XDR-TB transmission clusters suggesting a high level of clonal expansion across four Thai provinces. All XDR-TB (100%) were likely due to transmission rather than inadequate treatment. We found a 47 mutation signature and a partial deletion of the fadD14 gene in the circulating XDR-TB cluster, which can be used for surveillance of this rare and resilient M. tuberculosis strain-type that is causing increasing health burden. We also detected three novel deletion positions, a deletion of 1285 bp within desA3 (Rv3230c), large deletions in the plcB, plcA, and ppe38 gene which may play a role in the virulence, pathogenesis or evolution of the L2.1 strain-type.

Unexpected diversity of CRISPR unveils some evolutionary patterns of repeated sequences in Mycobacterium tuberculosis

BACKGROUND

Diversity of the CRISPR locus of Mycobacterium tuberculosis complex has been studied since 1997 for molecular epidemiology purposes. By targeting solely the 43 spacers present in the two first sequenced genomes (H37Rv and BCG), it gave a biased idea of CRISPR diversity and ignored diversity in the neighbouring cas-genes.

RESULTS

We set up tailored pipelines to explore the diversity of CRISPR-cas locus in Short Reads. We analyzed data from a representative set of 198 clinical isolates as evidenced by well-characterized SNPs. We found a relatively low diversity in terms of spacers: we recovered only the 68 spacers that had been described in 2000. We found no partial or global inversions in the sequences, letting always the Direct Variant Repeats (DVR) in the same order. In contrast, we found an unexpected diversity in the form of: SNPs in spacers and in Direct Repeats, duplications of various length, and insertions at various locations of the IS6110 insertion sequence, as well as blocks of DVR deletions. The diversity was in part specific to lineages. When reconstructing evolutionary steps of the locus, we found no evidence for SNP reversal. DVR deletions were linked to recombination between IS6110 insertions or between Direct Repeats.

CONCLUSION

This work definitively shows that CRISPR locus of M. tuberculosis did not evolve by classical CRISPR adaptation (incorporation of new spacers) since the last most recent common ancestor of virulent lineages. The evolutionary mechanisms that we discovered could be involved in bacterial adaptation but in a way that remains to be identified.

Geno-Spatial Distribution of Mycobacterium Tuberculosis and Drug Resistance Profiles in Myanmar-Thai Border Area

Worldwide, studies investigating the relationship between the lineage of Mycobacterium tuberculosis (MTB) across geographic areas has empowered the “End TB” program and understand transmission across national boundaries. Genomic diversity of MTB varies with geographical locations and ethnicity. Genomic diversity can also affect the emergence of drug resistance. In Myanmar, we still have limited genetic information about geographical, ethnicity, and drug resistance linkage to MTB genetic information. This study aimed to describe the geno-spatial distribution of MTB and drug resistance profiles in Myanmar–Thailand border areas. A cross-sectional study was conducted with a total of 109 sequenced isolates. The lineages of MTB and the potential associated socio-demographic, geographic and clinical factors were analyzed using Fisher’s exact tests. p value of statistically significance was set at < 0.05. We found that 67% of the isolates were lineage 1 (L1)/East-African-Indian (EAI) (n = 73), followed by lineage 2 (L2)/Beijing (n = 26), lineage 4 (L4)/European American (n = 6) and lineage 3 (L3)/Delhi/Central Asian (n = 4). “Gender”, “type of TB patient”, “sputum smear grading” and “streptomycin resistance” were significantly different with the lineages of MTB. Sublineages of L1, which had never been reported elsewhere in Myanmar, were detected in this study area. Moreover, both ethnicity and lineage of MTB significantly differed in distribution by patient location. Diversity of the lineage of MTB and detection of new sublineages suggested that this small area had been resided by a heterogeneous population group who actively transmitted the disease. This information on distribution of lineage of MTB can be linked in the future with those on the other side of the border to evaluate cross-border transmission.

24-locus MIRU-VNTR and Spoligotyping analysis of drug-resistant Mycobacterium tuberculosis strains isolated from Northeastern Thailand

Tuberculosis, caused by Mycobacterium tuberculosis (MTB) infection, remains a global health problem with increased concerns due to drug-resistant tuberculosis. However, molecular genotyping profiles may give insight of the transmission of TB in a particular region. The present study aimed to characterize the genetic diversity of drug-resistant MTB and evaluate primer sets applied for the epidemiological study of circulating MTB in Northeastern Thailand. A total of 92 MTB isolates, resistant to rifampicin and/or isoniazid, were collected from the Office of Disease Prevention and Control between 2013 and 2016. All isolates were genotyped by 24-locus MIRU-VNTR typing combined with spoligotyping. We also analyzed the distributions of drug susceptibility pattern and demographic data among different genotypes. In comparison with different loci sets, discriminatory power based on 12, 15, 24 standard primers were investigated. Eighty-six particular profiles were found; among the patterns, two clusters were produced in 8 strains. East African Indians (EAI) were the most prevalent strains (33 isolates, 35.87%) followed by Beijing (30 isolates, 32.61%), with 23 unknown isolates strains also found. The HGDI based on combination of 24 loci analysis and spoligotyping was 0.9962. The number of tandem repeat generated was highly discriminant (HGDI>0.6) at locus 580 (0.66), 960 (0.67), 2163b (0.73), 2165 (0.62), 2461 (0.68) 3690 (0.73) and 4052 (0.79), respectively. In contrast, the diversity at locus 154 and 2059 was not revealed. The results emphasized that 24-locus MIRU-VNTR and spoligotyping could be useful for epidemiological surveillance of drug-resistant MTB in this region. At a given allelic diversity, 7 primer sets containing MIRU04, MIRU10, QUB2163b, ETRA, ETRB, Mtub39 and QUB26 may be considered for screening the VNTR patterns. In addition, this study gathered both demographics and genotypic data within the same investigation for further tuberculosis prevention and control.

Homoplastic single nucleotide polymorphisms contributed to phenotypic diversity in Mycobacterium tuberculosis

Homoplastic mutations are mutations independently occurring in different clades of an organism. The homoplastic changes may be a result of convergence evolution due to selective pressures. Reports on the analysis of homoplastic mutations in Mycobacterium tuberculosis have been limited. Here we characterized the distribution of homoplastic single nucleotide polymorphisms (SNPs) among genomes of 1,170 clinical M. tuberculosis isolates. They were present in all functional categories of genes, with pe/ppe gene family having the highest ratio of homoplastic SNPs compared to the total SNPs identified in the same functional category. Among the pe/ppe genes, the homoplastic SNPs were common in a relatively small number of homologous genes, including ppe18, the protein of which is a component of a promising candidate vaccine, M72/AS01E. The homoplastic SNPs in ppe18 were particularly common among M. tuberculosis Lineage 1 isolates, suggesting the need for caution in extrapolating the results of the vaccine trial to the population where L1 is endemic in Asia. As expected, homoplastic SNPs strongly associated with drug resistance. Most of these mutations are already well known. However, a number of novel mutations associated with streptomycin resistance were identified, which warrants further investigation. A SNP in the intergenic region upstream of Rv0079 (DATIN) was experimentally shown to increase transcriptional activity of the downstream gene, suggesting that intergenic homoplastic SNPs should have effects on the physiology of the bacterial cells. Our study highlights the potential of homoplastic mutations to produce phenotypic changes. Under selective pressure and during interaction with the host, homoplastic mutations may confer advantages to M. tuberculosis and deserve further characterization.

Bayesian reconstruction of Mycobacterium tuberculosis transmission networks in a high incidence area over two decades in Malawi reveals associated risk factors and genomic variants

Understanding host and pathogen factors that influence tuberculosis (TB) transmission can inform strategies to eliminate the spread of Mycobacterium tuberculosis (Mtb). Determining transmission links between cases of TB is complicated by a long and variable latency period and undiagnosed cases, although methods are improving through the application of probabilistic modelling and whole-genome sequence analysis. Using a large dataset of 1857 whole-genome sequences and comprehensive metadata from Karonga District, Malawi, over 19 years, we reconstructed Mtb transmission networks using a two-step Bayesian approach that identified likely infector and recipient cases, whilst robustly allowing for incomplete case sampling. We investigated demographic and pathogen genomic variation associated with transmission and clustering in our networks. We found that whilst there was a significant decrease in the proportion of infectors over time, we found higher transmissibility and large transmission clusters for lineage 2 (Beijing) strains. By performing evolutionary convergence testing (phyC) and genome-wide association analysis (GWAS) on transmitting versus non-transmitting cases, we identified six loci, PPE54, accD2, PE_PGRS62, rplI, Rv3751 and Rv2077c, that were associated with transmission. This study provides a framework for reconstructing large-scale Mtb transmission networks. We have highlighted potential host and pathogen characteristics that were linked to increased transmission in a high-burden setting and identified genomic variants that, with validation, could inform further studies into transmissibility and TB eradication.

TbD1 deletion as a driver of the evolutionary success of modern epidemic Mycobacterium tuberculosis lineages

Mycobacterium tuberculosis (Mtb) strains are classified into different phylogenetic lineages (L), three of which (L2/L3/L4) emerged from a common progenitor after the loss of the MmpS6/MmpL6-encoding Mtb-specific deletion 1 region (TbD1). These TbD1-deleted “modern” lineages are responsible for globally-spread tuberculosis epidemics, whereas TbD1-intact “ancestral” lineages tend to be restricted to specific geographical areas, such as South India and South East Asia (L1) or East Africa (L7). By constructing and characterizing a panel of recombinant TbD1-knock-in and knock-out strains and comparison with clinical isolates, here we show that deletion of TbD1 confers to Mtb a significant increase in resistance to oxidative stress and hypoxia, which correlates with enhanced virulence in selected cellular, guinea pig and C3HeB/FeJ mouse infection models, the latter two mirroring in part the development of hypoxic granulomas in human disease progression. Our results suggest that loss of TbD1 at the origin of the L2/L3/L4 Mtb lineages was a key driver for their global epidemic spread and outstanding evolutionary success.

Mycobacterium tuberculosis (Mtb) modern strains emerged from a common progenitor after the loss of Mtb-specific deletion 1 region (TbD1). Here, the authors show that deletion of TbD1 correlates with enhanced Mtb virulence in animal models, mirroring the development of hypoxic granulomas in human disease progression.

Genotyping of Mycobacterium tuberculosis spreading in Hanoi, Vietnam using conventional and whole genome sequencing methods

Hanoi is the capital of Vietnam, one of the 30 countries with a high tuberculosis (TB) burden. Fundamental data on the molecular epidemiology of the disease is required for future TB management. To identify lineages and genotypes of Mycobacterium tuberculosis (Mtb), conventional genotyping data from clinical isolates of the Hanoi area was compared with whole genome sequencing (WGS) analysis from 332 of 470 samples. It was obtained from lineage-specific single nucleotide variants (SNVs), large sequence polymorphisms, spoligotyping, and variable number of tandem repeats (VNTR) analysis using mycobacterial interspersed repetitive unit (MIRU) and Japan anti-tuberculosis association (JATA) locus sets. This information was directly compared with results obtained from WGS. Mini-satellite repeat unit variants were identified using BLAST search against concatenated short read sequences, the RepUnitTyping tool. WGS analysis revealed that the Mtb strains tested are diverse and classified into lineage (L) 1, 2 and 4 (24.7, 57.2 and 18.1% respectively). The majority of the L2 strains were further divided into ancient and modern Beijing genotypes, and most of the L1 group were EAI4_VNM strains. Although conventional PCR-based genotyping results were mostly consistent with information obtained through WGS analysis, in-depth analysis identified aberrant deletions and spacers that may cause discordance. JATA-VNTR sets, including hypervariable loci, separated large Beijing genotypic clusters generated by MIRU15 into smaller groups. The distribution of repeat unit variants observed within 33 VNTR loci showed clear variation depending on the three lineages. WGS-based pairwise-SNV differences within VNTR-defined genotypic clusters were greater in L1 than in L2 and L4 (P = .001). Direct comparisons between results of PCR-based genotyping and in silico analysis of WGS data would bridge a gap between classical and modern technologies during this transition period, and provide further information on Mtb genotypes in specific geographical areas.

A novel Ancestral Beijing sublineage of Mycobacterium tuberculosis suggests the transition site to Modern Beijing sublineages

Global Mycobacterium tuberculosis population comprises 7 major lineages. The Beijing strains, particularly the ones classified as Modern groups, have been found worldwide, frequently associated with drug resistance, younger ages, outbreaks and appear to be expanding. Here, we report analysis of whole genome sequences of 1170 M. tuberculosis isolates together with their patient profiles. Our samples belonged to Lineage 1-4 (L1-L4) with those of L1 and L2 being equally dominant. Phylogenetic analysis revealed several new or rare sublineages. Differential associations between sublineages of M. tuberculosis and patient profiles, including ages, ethnicity, HIV (human immunodeficiency virus) infection and drug resistance were demonstrated. The Ancestral Beijing strains and some sublineages of L4 were associated with ethnic minorities while L1 was more common in Thais. L2.2.1.Ancestral 4 surprisingly had a mutation that is typical of the Modern Beijing sublineages and was common in Akha and Lahu tribes who have migrated from Southern China in the last century. This may indicate that the evolutionary transition from the Ancestral to Modern Beijing sublineages might be gradual and occur in Southern China, where the presence of multiple ethnic groups might have allowed for the circulations of various co-evolving sublineages which ultimately lead to the emergence of the Modern Beijing strains.

Genome-wide mutational biases fuel transcriptional diversity in the Mycobacterium tuberculosis complex

The Mycobacterium tuberculosis complex (MTBC) members display different host-specificities and virulence phenotypes. Here, we have performed a comprehensive RNAseq and methylome analysis of the main clades of the MTBC and discovered unique transcriptional profiles. The majority of genes differentially expressed between the clades encode proteins involved in host interaction and metabolic functions. A significant fraction of changes in gene expression can be explained by positive selection on single mutations that either create or disrupt transcriptional start sites (TSS). Furthermore, we show that clinical strains have different methyltransferases inactivated and thus different methylation patterns. Under the tested conditions, differential methylation has a minor direct role on transcriptomic differences between strains. However, disruption of a methyltransferase in one clinical strain revealed important expression differences suggesting indirect mechanisms of expression regulation. Our study demonstrates that variation in transcriptional profiles are mainly due to TSS mutations and have likely evolved due to differences in host characteristics.

Mycobacterium tuberculosis lineage 1 genetic diversity in Pará, Brazil, suggests common ancestry with east-African isolates potentially linked to historical slave trade

Lineage 1 (L1) is one of seven Mycobacterium tuberculosis complex (MTBC) lineages. The objective of this study was to improve the complex taxonomy of L1 using phylogenetic SNPs, and to look for the origin of the main L1 sublineage prevalent in Para, Brazil. We developed a high-throughput SNPs-typing assay based on 12-L1-specific SNPs. This assay allowed us to experimentally retrieve SNP patterns on nine of these twelve SNPs in 277 isolates previously tentatively assigned to L1 spoligotyping-based sublineages. Three collections were used: Pará-Brazil (71); RIVM, the Netherlands (102), Madagascar (104). One-hundred more results were generated in Silico using the PolyTB database. Based on the final SNPs combination, the samples were classified into 11 clusters (C1-C11). Most isolates within a SNP-based cluster shared a mutual spoligotyping-defined lineage. However, L1/EAI1-SOM (SIT48) and L1/EAI6-BGD1 (SIT591) showed a poor correlation with SNP data and are not monophyletic. L1/EAI8-MDG and L1/EAI3-IND belonged to C5; this result suggests that they share a common ancestor. L1.1.3/SIT129, a spoligotype pattern found in SNPs-cluster C6, was found to be shared between Pará/Brazil and Malawi. SIT129 was independently found to be highly prevalent in Mozambique, which suggests a migration history from East-Africa to Brazil during the 16th-18th slave trade period to Northern Brazil.

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.