Automated high-throughput genotyping for study of global epidemiology of Mycobacterium tuberculosis based on mycobacterial interspersed repetitive units

Molecular epidemiology of Mycobacterium tuberculosis in Ecuador: Recent advances and future challenges

Tuberculosis (TB) is one of the three leading causes of death from a single infectious agent, Mycobacterium tuberculosis (MTB), together with COVID-19 and HIV/AIDS. This disease places a heavy burden on countries with low socio-economic development and aggravates existing inequalities. For the year 2021, estimations for Ecuador were 8500 TB cases, of which 370 were associated to multiple drug resistance (TB-MDR), and 1160 deaths. In the same year, Ecuador notified 5973 total cases, 401 of them were TB-MDR, pointing out an under diagnosis problem. The few molecular epidemiology studies available conclude that L4 is the most prevalent MTB lineage in Ecuador (with LAM as the main L4 sublineage), but L2-Beijing family is also present at low prevalence. Nevertheless, with less than 1 % MTB isolates genetically characterized by either MIRU-VNTR, spolygotyping or WGS to date, molecular epidemiology research must me improved to assist the TB surveillance and control program in Ecuador.

Bridging the gap between molecular and genomic epidemiology in tuberculosis: inferring MIRU-VNTR patterns from genomic data

The transition from MIRU-VNTR-based epidemiology studies in tuberculosis (TB) to genomic epidemiology has transformed how we track transmission. However, short-read sequencing is poor at analyzing repetitive regions such as the MIRU-VNTR loci. This causes a gap between the new genomic data and the large amount of information stored in historical databases. Long-read sequencing could bridge this knowledge gap by allowing analysis of repetitive regions. However, the feasibility of extracting MIRU-VNTRs from long reads and linking them to historical data has not been evaluated. In our study, an in silico arm, consisting of inference of MIRU patterns from long-read sequences (using MIRUReader program), was compared with an experimental arm, involving standard amplification and fragment sizing. We analyzed overall performance on 39 isolates from South Africa and confirmed reproducibility in a sample enriched with 62 clustered cases from Spain. Finally, we ran 25 consecutive incident cases, demonstrating the feasibility of correctly assigning new clustered/orphan cases by linking data inferred from genomic analysis to MIRU-VNTR databases. Of the 3,024 loci analyzed, only 11 discrepancies (0.36%) were found between the two arms: three attributed to experimental error and eight to misassigned alleles from long-read sequencing. A second round of analysis of these discrepancies resulted in agreement between the experimental and in silico arms in all but one locus. Adjusting the MIRUReader program code allowed us to flag potential in silico misassignments due to suboptimal coverage or unfixed double alleles. Our study indicates that long-read sequencing could help address potential chronological and geographical gaps arising from the transition from molecular to genomic epidemiology of tuberculosis.

IMPORTANCE

The transition from molecular epidemiology in tuberculosis (TB), based on the analysis of repetitive regions (VNTR-based genotyping), to genomic epidemiology transforms in the precision with which we track transmission. However, short-read sequencing, the most common method for performing genomic analysis, is poor at analyzing repetitive regions. This means that we face a gap between the new genomic data and the large amount of information stored in historical databases, which is also an obstacle to cross-national surveillance involving settings where only molecular data are available. Long-read sequencing could help bridge this knowledge gap by allowing analysis of repetitive regions. Our study demonstrates that MIRU-VNTR patterns can be successfully inferred from long-read sequences, allowing the correct assignment of new cases as clustered/orphan by linking new data extracted from genomic analysis to historical MIRU-VNTR databases. Our data may provide a starting point for bridging the knowledge gap between the molecular and genomic eras in tuberculosis epidemiology.

SAGAS: Simulated annealing and greedy algorithm scheduler for laboratory automation

During laboratory automation of life science experiments, coordinating specialized instruments and human experimenters for various experimental procedures is important to minimize the execution time. In particular, the scheduling of life science experiments requires the consideration of time constraints by mutual boundaries (TCMB) and can be formulated as the "scheduling for laboratory automation in biology" (S-LAB) problem. However, existing scheduling methods for the S-LAB problems have difficulties in obtaining a feasible solution for large-size scheduling problems at a time sufficient for real-time use. In this study, we proposed a fast schedule-finding method for S-LAB problems, SAGAS (Simulated annealing and greedy algorithm scheduler). SAGAS combines simulated annealing and the greedy algorithm to find a scheduling solution with the shortest possible execution time. We have performed scheduling on real experimental protocols and shown that SAGAS can search for feasible or optimal solutions in practicable computation time for various S-LAB problems. Furthermore, the reduced computation time by SAGAS enables us to systematically search for laboratory automation with minimum execution time by simulating scheduling for various laboratory configurations. This study provides a convenient scheduling method for life science automation laboratories and presents a new possibility for designing laboratory configurations.

Hi-plex deep amplicon sequencing for identification, high-resolution genotyping and multidrug resistance prediction of Mycobacterium leprae directly from patient biopsies by using Deeplex Myc-Lep

BACKGROUND

Expansion of antimicrobial resistance monitoring and epidemiological surveillance are key components of the WHO strategy towards zero leprosy. The inability to grow Mycobacterium leprae in vitro precludes routine phenotypic drug susceptibility testing, and only limited molecular tests are available. We evaluated a culture-free targeted deep sequencing assay, for mycobacterial identification, genotyping based on 18 canonical SNPs and 11 core variable-number tandem-repeat (VNTR) markers, and detection of rifampicin, dapsone and fluoroquinolone resistance-associated mutations in rpoB/ctpC/ctpI, folP1, gyrA/gyrB, respectively, and hypermutation-associated mutations in nth.

METHODS

The limit of detection (LOD) was determined using DNA of M. leprae reference strains and from 246 skin biopsies and 74 slit skin smears of leprosy patients, with genome copies quantified by RLEP qPCR. Sequencing results were evaluated versus whole genome sequencing (WGS) data of 14 strains, and versus VNTR-fragment length analysis (FLA) results of 89 clinical specimens.

FINDINGS

The LOD for sequencing success ranged between 80 and 3000 genome copies, depending on the sample type. The LOD for minority variants was 10%. All SNPs detected in targets by WGS were identified except in a clinical sample where WGS revealed two dapsone resistance-conferring mutations instead of one by Deeplex Myc-Lep, due to partial duplication of the sulfamide-binding domain in folP1. SNPs detected uniquely by Deeplex Myc-Lep were missed by WGS due to insufficient coverage. Concordance with VNTR-FLA results was 99.4% (926/932 alleles).

INTERPRETATION

Deeplex Myc-Lep may help improve the diagnosis and surveillance of leprosy. Gene domain duplication is an original putative drug resistance-related genetic adaptation in M. leprae.

FUNDING

EDCTP2 programme supported by the European Union (grant number RIA2017NIM-1847 -PEOPLE). EDCTP, R2Stop: Effect:Hope, The Mission To End Leprosy, the Flemish Fonds Wetenschappelijk Onderzoek.

Genetic Diversity and Primary Drug Resistance of Mycobacterium tuberculosis Beijing Genotype Strains in Northwestern Russia

The Beijing genotype is the main family of Mycobacterium tuberculosis in Russia. We analyzed its diversity and drug resistance in provinces across Northwestern Russia to identify the epidemiologically relevant Beijing strains. The study collection included 497 isolates from newly-diagnosed tuberculosis (TB) patients. Bacterial isolates were subjected to drug-susceptibility testing and genotyping. The Beijing genotype was detected in 57.5% (286/497); 50% of the Beijing strains were multidrug-resistant (MDR). Central Asian/Russian and B0/W148 groups included 176 and 77 isolates, respectively. MDR was more frequent among B0/W148 strains compared to Central Asian/Russian strains (85.7% vs. 40.3%, p < 0.0001). Typing of 24 minisatellite loci of Beijing strains revealed 82 profiles; 230 isolates were in 23 clusters. The largest Central Asian/Russian types were 94-32 (n = 75), 1065-32 (n = 17), and 95-32 (n = 12). B0/W148 types were 100-32 (n = 59) and 4737-32 (n = 5). MDR was more frequent in types 1065-32 (88.2%), 100-32 (83.1%), and 4737-32 (100%). In contrast, type 9391-32 (n = 9) included only drug-susceptible strains. To conclude, M. tuberculosis Beijing genotype is dominant in Northwestern Russia, and an active transmission of overwhelmingly MDR B0/W148 types explains the reported increase of MDR-TB. The presence of MDR-associated minor variants (type 1071-32/ancient Beijing and Central Asia Outbreak strain) in some of the studied provinces also requires attention.

Whole genome sequencing for tuberculosis in Victoria, Australia: A genomic implementation study from 2017 to 2020

Background

Whole genome sequencing (WGS) is increasingly used by tuberculosis (TB) programs to monitor Mycobacterium tuberculosis (Mtb) transmission. We aimed to characterise the molecular epidemiology of TB and Mtb transmission in the low-incidence setting of Victoria, Australia, and assess the utility of WGS.

Methods

WGS was performed on all first Mtb isolates from TB cases from 2017 to 2020. Potential clusters (≤12 single nucleotide polymorphisms [SNPs]) were investigated for epidemiological links. Transmission events in highly-related (≤5 SNPs) clusters were classified as likely or possible, based on the presence or absence of an epidemiological link, respectively. Case characteristics and transmission settings (as defined by case relationship) were summarised. Poisson regression was used to examine associations with secondary case number.

Findings

Of 1844 TB cases, 1276 (69.2%) had sequenced isolates, with 182 (14.2%) in 54 highly-related clusters, 2-40 cases in size. Following investigation, 140 cases (11.0% of sequenced) were classified as resulting from likely/possible local-transmission, including 82 (6.4%) for which transmission was likely. Common identified transmission settings were social/religious (26.4%), household (22.9%) and family living in different households (7.1%), but many were uncertain (41.4%). While household transmission featured in many clusters (n = 24), clusters were generally smaller (median = 3 cases) than the fewer that included transmission in social/religious settings (n = 12, median = 7.5 cases). Sputum-smear-positivity was associated with higher secondary case numbers.

Interpretation

WGS results suggest Mtb transmission commonly occurs outside the household in our low-incidence setting. Further work is required to optimise the use of WGS in public health management of TB.

Funding

The Victorian Tuberculosis Program receives block funding for activities including case management and contact tracing from the Victorian Department of Health. No specific funding for this report was received by manuscript authors or the Victorian Tuberculosis Program, and the funders had no role in the study design, data collection, data analysis, interpretation or report writing.

Compilation of 10 Years of MIRU-VNTR Data: Canadian National Tuberculosis Laboratory's Experience

Tuberculosis is a significant cause of morbidity worldwide and is a priority at the provincial and federal levels in Canada. It is known that tuberculosis transmission networks are complex and span many years as well as different jurisdictions and countries. MIRU-VNTR is a universal tuberculosis genotyping method that utilizes a 24-loci pattern and it has shown promise in identifying inter and intrajurisdictional clusters within Canada. MIRU-VNTR data collected over 10 years from the National Reference Centre for Mycobacteriology (NRCM) were analyzed in this study. Some clusters were unique to a single province/territory, while others spanned multiple provinces and/or territories in Canada. The use of a universal laboratory test can enhance contact tracing, provide geographical information on circulating genotypes, and hence, aid in tuberculosis investigation by public health. The housing of all data on one platform, technical ease of the method, easy exchange of data between jurisdictions, and strong collaboration with laboratories and surveillance units at the provincial and federal levels have the potential to identify possible outbreaks in real time.

Genetic diversity of drug-resistant Mycobacterium tuberculosis clinical isolates from Khuzestan province, Iran

The emergence of drug-resistant strains of the Mycobacterium tuberculosis (MTB) has challenged tuberculosis control programs. So far, few studies using the 24-locus mycobacterial interspersed repetitive unit variable number tandem repeats (MIRU-VNTR) have investigated the genetic diversity of MTB in Iran. This study aimed to determine the genetic diversity of MTB isolates resistant to first-line anti-tuberculosis drugs using 24-locus MIRU-VNTR in southwestern Iran. Out of 6620 MTB clinical isolates, 29 resistant isolates to one or more isoniazid, rifampin, and ethambutol were detected using drug susceptibility testing by the proportional method. The manual 24-locus MIRU-VNTR was used to determine the MTB resistant isolates’ phylogenetic relationship. MIRU-VNTRplus web application tools were applied to analyze the associated data. Using 24-locus MIRU-VNTR, 13.8% of isolates (n = 4) were distributed in two clusters, and the remaining 86.2% (n = 25) showed a unique pattern. Four clonal complexes were observed in the minimum spanning tree based on the double-locus variant. Most isolates belonged to Delhi/CAS (34.5%, 10/29) and NEW-1 (24.1%, 7/29) sub-lineages, followed by EAI and LAM with a frequency of 6.9% (2/29) and 3.5% (1/29), respectively. Eight isolates (27.6%) did not match any genotype in the database. The 24-locus MIRU-VNTR showed a high discriminatory power; however, the 15-locus and 12-locus set analyses were more discriminative. Our study revealed a high degree of genetic diversity among drug-resistant MTB isolates, which could be interpreted as the low rate of person-to-person transmission in this region. The 15-locus MIRU-VNTR would be recommended for preliminary genotyping of drug-resistant MTB.

Detection of Mycobacterium tuberculosis multiple strains in sputum samples from patients with pulmonary tuberculosis in south western Uganda using MIRU-VNTR

Infections with multiple strains of Mycobacterium tuberculosis are now widely recognized as a common occurrence. Identification of patients infected with multiple strains provides both insight into the disease dynamics and the epidemiology of tuberculosis. Analysis of Mycobacterial Interspersed Repetitive Unit-Variable-Number Tandem Repeats (MIRU-VNTR) has been shown to be highly sensitive in detecting multiple M. tuberculosis strains even in sputum. The goal of this study was to identify cases of multiple M. tuberculosis strain infections among patients diagnosed with pulmonary tuberculosis in Southwestern Uganda and assessment of factors associated with multiple strain infections. DNA extracted directly from 78 sputum samples, each from an individual patient, was analyzed using the standard 24 loci MIRU-VNTR typing. Five (6.4%) of the 78 patients were infected with multiple strains of M. tuberculosis with all of them being the newly diagnosed cases while two-thirds of them were co-infected with HIV. Exact regression analysis projected that the natives were more likely to harbor multiple strains (OR; 0.981, 95% CI 0–7.926) as well as those with a high microbial load (OR; 0.390, 95% CI 0–3.8167). Despite these findings being not statistically significant due to the small sample size, this points to a critical component of disease dynamics that has clinical implications and emphasizes a need for a study using a larger cohort. It is also essential to study the potential factors associated with higher risk of exposure to newly diagnosed and HIV positive patients at the community level. In addition, our ability to detect multiple M. tuberculosis strains using the standard 24 loci MIRU-VNTR typing especially with allelic diversity in loci 2059 and 3171, which are excluded from the 15-locus MIRU-VNTR, lead us to recommend the use of this genotyping technique, especially in areas with tuberculosis endemicity similar to this study.

Optimal Scheduling for Laboratory Automation of Life Science Experiments with Time Constraints

In automated laboratories consisting of multiple different types of instruments, scheduling algorithms are useful for determining the optimal allocations of instruments to minimize the time required to complete experimental procedures. However, previous studies on scheduling algorithms for laboratory automation have not emphasized the time constraints by mutual boundaries (TCMBs) among operations, which is important in procedures involving live cells or unstable biomolecules. Here, we define the "scheduling for laboratory automation in biology" (S-LAB) problem as a scheduling problem for automated laboratories in which operations with TCMBs are performed by multiple different instruments. We formulate an S-LAB problem as a mixed-integer programming (MIP) problem and propose a scheduling method using the branch-and-bound algorithm. Simulations show that our method can find the optimal schedules of S-LAB problems that minimize overall execution time while satisfying the TCMBs. Furthermore, we propose the use of our scheduling method for the simulation-based design of job definitions and laboratory configurations.

A Prospective Cohort Study on the Prevalent and Recurrent Tuberculosis Isolates Using the MIRU-VNTR Typing

The study aims to describe the clustering characteristics of Mycobacterium tuberculosis (M.tb) strains circulating in eastern China and determine the ratio of relapse and reinfection in recurrent patients. We recruited sputum smear-positive pulmonary tuberculosis cases from five cities of Jiangsu Province, China, during August 2013 and December 2015. Patients were followed for the treatment outcomes and recurrence based on a cohort design. M.tb strains were isolated and genotyped using the 12-locus MIRU-VNTR. The Beijing family was identified by the extended Region of Difference (RD) analysis. The Hunter-Gaston Discriminatory Index (HGDI) was used to judge the resolution ability of MIRU-VNTR. The odds ratio (OR) together with 95% confidence interval (CI) were used to estimate the strength of association. We performed a cluster analysis on 2098 M.tb isolates and classified them into 545 genotypes and five categories (I, 0.19%; II, 0.43%; III, 3.34%; IV, 77.46%; V, 18.59%). After adjusting for potential confounders, the Beijing family genotype (OR = 118.63, 95% CI: 79.61–176.79, P = 0.001) was significantly related to the dominant strain infections. Patients infected with non-dominant strains had a higher risk of the pulmonary cavity (OR = 1.39, 95% CI: 1.01–1.91, P = 0.046). Among 37 paired recurrent cases, 22 (59.46%) were determined as endogenous reactivation, and 15 (40.54%) were exogenous reinfection. The type of M.tb strains prevalent in Jiangsu Province is relatively single. Beijing family strains infection is dominant in local tuberculosis cases. Endogenous reactivation appears to be a major cause of recurrent tuberculosis in Eastern China. This finding emphasizes the importance of case follow-up and monitoring after the completion of antituberculosis treatment.

Drug resistant tuberculosis cases from the Copperbelt province and Northern regions of Zambia: Genetic diversity, demographic and clinical characteristics

Tuberculosis (TB) caused by Mycobacterium tuberculosis remains a major cause of death worldwide. Diverse genotypes have been demonstrated to drive the epidemiology of drug resistant (DR-) TB globally. Currently, there is limited knowledge on the genotypes and transmission dynamics of M. tuberculosis in Zambia. This study aimed to describe the genotypes of DR-TB from the Copperbelt and Northern regions of Zambia. Molecular typing tools of insertion sequence 6110-restriction fragment length polymorphism (IS6110-RFLP) and spacer oligonucleotide typing (spoligotyping) were applied. We demonstrate that diverse genotypes are associated with DR-TB in Zambia. The predominant genotype was lineage 4; other strains belonged to lineage 2 and 3. Genotypes previously identified as driving the epidemiology of drug susceptible TB have been identified as drivers of DR-TB. Genotyping analysis showed clustering of strains among patients from different regions of the country; suggesting that DR-TB is widespread. Molecular findings combined with phenotypic and epidemiologic findings play a critical role in identifying circulating genotypes and possible transmission chains. Clustering of drug resistant strains was demonstrated to be 48% and 86% according to IS6110-RFLP and spoligotyping, respectively. However, gaps in clinical and demographic data skew the interpretation, and call for data collection policy improvements.

The Relevance of Genomic Epidemiology for Control of Tuberculosis in West Africa

Tuberculosis (TB), an airborne infectious disease caused by Mycobacterium tuberculosis complex (MTBC), remains a global health problem. West Africa has a unique epidemiology of TB that is characterized by medium- to high-prevalence. Moreover, the geographical restriction of M. africanum to the sub-region makes West Africa have an extra burden to deal with a two-in-one pathogen. The region is also burdened with low case detection, late reporting, poor treatment adherence leading to development of drug resistance and relapse. Sporadic studies conducted within the subregion report higher burden of drug resistant TB (DRTB) than previously thought. The need for more sensitive and robust tools for routine surveillance as well as to understand the mechanisms of DRTB and transmission dynamics for the design of effective control tools, cannot be overemphasized. The advancement in molecular biology tools including traditional fingerprinting and next generation sequencing (NGS) technologies offer reliable tools for genomic epidemiology. Genomic epidemiology provides in-depth insight of the nature of pathogens, circulating strains and their spread as well as prompt detection of the emergence of new strains. It also offers the opportunity to monitor treatment and evaluate interventions. Furthermore, genomic epidemiology can be used to understand potential emergence and spread of drug resistant strains and resistance mechanisms allowing the design of simple but rapid tools. In this review, we will describe the local epidemiology of MTBC, highlight past and current investigations toward understanding their biology and spread as well as discuss the relevance of genomic epidemiology studies to TB control in West Africa.

Genotypic diversity of multi- and pre-extremely drug-resistant Mycobacterium tuberculosis isolates from Morocco

In Morocco, the prevalence of multidrug resistant tuberculosis (MDR-TB) continues to increase especially within previously treated cases; these MDR cases may evolve to extensively drug resistant tuberculosis (XDR-TB) raising major concern to TB control programs. From an epidemiological window, scarce informations are available about the genetic diversity of Mycobacterium tuberculosis (MTB) strains fueling these forms of resistance. The aim of this study was to assess to genetic diversity of MDR-MTB strains. Hence, this prospective study was conducted on patients diagnosed with MDR-TB at Pasteur Institute of Casablanca from 2010 to 2013. A total of 70 MDR-MTB isolates were genotyped by spoligotyping and 15-loci MIRU-VNTR methods. Spoligotyping generated four orphan patterns, five unique profiles whereas 61 strains were grouped in nine clusters (2 to 25 strains per cluster), the clustering rates being 87.1%. Subtyping by 15 loci MIRU-VNTR splitted all clusters already established by spoligotyping and generated 70 unique profiles not recognized in SITVIT2 database; clustering rate was equal to zero. HGDI analysis of 15 loci MIRU demonstrated that eight out of 15 loci were highly discriminant. Of note, all pre-XDR strains belongs to many clades, meaning that there no association between gyrA mutants and particular clade. Overall, the data generated by this study (i) describe the population structure of MDR MTBC in Morocco which is highly homogenous, (ii) confirm that TB in Morocco is almost exclusively transmitted by modern and evolutionary lineages with high level of biodiversity seen by MIRU, and (iii) validate the use of optimized 15-loci MIRU-VNTR format for future investigations in Morocco.

Variable number tandem repeats mediate the expression of proximal genes

Variable number tandem repeats (VNTRs) account for significant genetic variation in many organisms. In humans, VNTRs have been implicated in both Mendelian and complex disorders, but are largely ignored by genomic pipelines due to the complexity of genotyping and the computational expense. We describe adVNTR-NN, a method that uses shallow neural networks to genotype a VNTR in 18 seconds on 55X whole genome data, while maintaining high accuracy. We use adVNTR-NN to genotype 10,264 VNTRs in 652 GTEx individuals. Associating VNTR length with gene expression in 46 tissues, we identify 163 "eVNTRs". Of the 22 eVNTRs in blood where independent data is available, 21 (95%) are replicated in terms of significance and direction of association. 49% of the eVNTR loci show a strong and likely causal impact on the expression of genes and 80% have maximum effect size at least 0.3. The impacted genes are involved in diseases including Alzheimer's, obesity and familial cancers, highlighting the importance of VNTRs for understanding the genetic basis of complex diseases.

Mycobacterium tuberculosis Beijing Lineage and Risk for Tuberculosis in Child Household Contacts, Peru

Few studies have prospectively compared the relative transmissibility and propensity to cause disease of Mycobacterium tuberculosis Beijing strains with other human-adapted strains of the M. tuberculosis complex. We assessed the effect of Beijing strains on the risk for M. tuberculosis infection and disease progression in 9,151 household contacts of 2,223 culture-positive pulmonary tuberculosis (TB) patients in Lima, Peru. Child contacts exposed to Beijing strains were more likely than child contacts exposed to non-Beijing strains to be infected at baseline, by 12 months of follow-up, and during follow-up. We noted an increased but nonsignificant tendency for child contacts to develop TB. Beijing strains were not associated with TB in adult contacts. These findings suggest that Beijing strains are more transmissible in children than are non-Beijing strains.

Genotypic diversity of Mycobacterium tuberculosis isolates and its association with drug-resistance status in Xinjiang, China

Xinjiang is a unique region of Central Asian part of China. It is widely noted for high tuberculosis burden and particularly for growing prevalence of drug resistance. Understanding genotypic distribution of Mycobacterium tuberculosis could help clarify unknown causes for the spread of drug-resistant strains. We analyzed 986 M. tuberculosis isolates collected from Xinjiang. Two genotyping schemes, i.e., spoligotyping and multiple-locus variable number tandem repeats (VNTR), were used to determine the phylogenetic lineages and their association with drug-resistances. The M. tuberculosis isolates studied displayed wide distribution of spoligotypic lineages, including Beijing, T, CAS, Ural, LAM, MANU, H, X, EAI, S, Microti, and BOV. The dominant Beijing lineage showed statistical difference from non-Beijing lineages in patients ages (P < 0.001), ethnic groups (P < 0.001) and resistance of three or more drugs (P = 0.008). Further analysis of the year of 2017 subset (n = 257) using VNTR scheme revealed an extremely high discrimination power (Hunter-Gaston discriminatory index = 0.9994). Cluster analysis showed a much lower recent transmission index (7.93%), indicating that the high drug-resistant tuberculosis in this region was mainly caused by reactivation or inappropriate therapy rather than by recent transmission. These data would be valuable for making and implementing policies for improving tuberculosis treatment and care in Xinjiang.

Human tuberculosis and Mycobacterium tuberculosis complex: A review on genetic diversity, pathogenesis and omics approaches in host biomarkers discovery

Mycobacterium tuberculosis complex (MTBC) refers to a group of mycobacteria encompassing nine members of closely related species that causes tuberculosis in animals and humans. Among the nine members, Mycobacterium tuberculosis (M. tuberculosis) remains the main causative agent for human tuberculosis that results in high mortality and morbidity globally. In general, MTBC species are low in diversity but exhibit distinctive biological differences and phenotypes among different MTBC lineages. MTBC species are likely to have evolved from a common ancestor through insertions/deletions processes resulting in species speciation with different degrees of pathogenicity. The pathogenesis of human tuberculosis is complex and remains poorly understood. It involves multi-interactions or evolutionary co-options between host factors and bacterial determinants for survival of the MTBC. Granuloma formation as a protection or survival mechanism in hosts by MTBC remains controversial. Additionally, MTBC species are capable of modulating host immune response and have adopted several mechanisms to evade from host immune attack in order to survive in humans. On the other hand, current diagnostic tools for human tuberculosis are inadequate and have several shortcomings. Numerous studies have suggested the potential of host biomarkers in early diagnosis of tuberculosis, in disease differentiation and in treatment monitoring. "Multi-omics" approaches provide holistic views to dissect the association of MTBC species with humans and offer great advantages in host biomarkers discovery. Thus, in this review, we seek to understand how the genetic variations in MTBC lead to species speciation with different pathogenicity. Furthermore, we also discuss how the host and bacterial players contribute to the pathogenesis of human tuberculosis. Lastly, we provide an overview of the journey of "omics" approaches in host biomarkers discovery in human tuberculosis and provide some interesting insights on the challenges and directions of "omics" approaches in host biomarkers innovation and clinical implementation.

First approach to the population structure of Mycobacterium tuberculosis complex in the indigenous population in Puerto Nariño-Amazonas, Colombia

Introduction

Tuberculosis affects vulnerable groups to a greater degree, indigenous population among them.

Objective

To determine molecular epidemiology of clinical isolates of Mycobacterium tuberculosis circulating in an indigenous population through Spoligotyping and 24-loci MIRU-VNTR.

Methodology

A descriptive cross-sectional study was conducted in 23 indigenous communities of Puerto Nariño-Amazonas, Colombia. Recovered clinical isolates were genotyped. For genotyping analyzes global SITVIT2 database and the MIRU-VNTRplus web portal were used.

Results

74 clinical isolates were recovered. Genotyping of clinical isolates by spoligotyping determined 5 different genotypes, all of them belonged to Euro-American lineage. By MIRU-VNTR typing, a total of 14 different genotypes were recorded. Furthermore, polyclonal infection was found in two patients from the same community. The combination of the two methodologies determined the presence of 19 genotypes, 8 formed clusters with 63 clinical isolates in total. Based on epidemiological information, it was possible to establish a potential chain of active transmission in 10/63 (15.9%) patients.

Conclusions

High genomic homogeneity was determined in the indigenous population suggesting possible chains of active transmission. The results obtained showed that specific genotypes circulating among the indigenous population of Colombia are significantly different from those found in the general population.

Methods for Detecting Mycobacterial Mixed Strain Infections-A Systematic Review

Mixed strain infection (MSI) refers to the concurrent infection of a susceptible host with multiple strains of a single pathogenic species. Known to occur in humans and animals, MSIs deserve special consideration when studying transmission dynamics, evolution, and treatment of mycobacterial diseases, notably tuberculosis in humans and paratuberculosis (or Johne's disease) in ruminants. Therefore, a systematic review was conducted to examine how MSIs are defined in the literature, how widespread the phenomenon is across the host species spectrum, and to document common methods used to detect such infections. Our search strategy identified 121 articles reporting MSIs in both humans and animals, the majority (78.5%) of which involved members of the Mycobacterium tuberculosis complex, while only a few (21.5%) examined non-tuberculous mycobacteria (NTM). In addition, MSIs exist across various host species, but most reports focused on humans due to the extensive amount of work done on tuberculosis. We reviewed the strain typing methods that allowed for MSI detection and found a few that were commonly employed but were associated with specific challenges. Our review notes the need for standardization, as some highly discriminatory methods are not adapted to distinguish between microevolution of one strain and concurrent infection with multiple strains. Further research is also warranted to examine the prevalence of NTM MSIs in both humans and animals. In addition, it is envisioned that the accurate identification and a better understanding of the distribution of MSIs in the future will lead to important information on the epidemiology and pathophysiology of mycobacterial diseases.

Long-term molecular surveillance provides clues on a cattle origin for Mycobacterium bovis in Portugal

Animal tuberculosis (TB), caused by Mycobacterium bovis, is maintained in Portugal in a multi-host system, with cattle, red deer and wild boar, playing a central role. However, the ecological processes driving transmission are not understood. The main aim of this study was thus to contribute to the reconstruction of the spatiotemporal history of animal TB and to refine knowledge on M. bovis population structure in order to inform novel intervention strategies. A collection of 948 M. bovis isolates obtained during long-term surveillance (2002–2016, 15 years) of cattle (n = 384), red deer (n = 303) and wild boar (n = 261), from the main TB hotspot areas, was characterized by spoligotyping and 8 to 12-loci MIRU-VNTR. Spoligotyping identified 64 profiles and MIRU-VNTR distinguished 2 to 36 subtypes within each spoligotype, enabling differentiation of mixed or clonal populations. Common genotypic profiles within and among livestock and wildlife in the same spatiotemporal context highlighted epidemiological links across hosts and regions, as for example the SB0119-M205 genotype shared by cattle in Beja district or SB0121-M34 shared by the three hosts in Castelo Branco and Beja districts. These genomic data, together with metadata, were integrated in a Bayesian inference framework, identifying five ancestral M. bovis populations. The phylogeographic segregation of M. bovis in specific areas of Portugal where the disease persists locally is postulated. Concurrently, robust statistics indicates an association of the most probable ancient population with cattle and Beja, providing a clue on the origin of animal TB epidemics. This relationship was further confirmed through a multinomial probability model that assessed the influence of host species on spatiotemporal clustering. Two significant clusters were identified, one that persisted between 2004 and 2010, in Beja district, with Barrancos county at the centre, overlapping the central TB core area of the Iberian Peninsula, and highlighting a significant higher risk associated to cattle. The second cluster was predominant in the 2012–2016 period, holding the county Rosmaninhal at the centre, in Castelo Branco district, for which wild boar contributed the most in relative risk. These results provide novel quantitative insights beyond empirical perceptions, that may inform adaptive TB control choices in different regions.

Circulating strains of Mycobacterium tuberculosis: 24 loci MIRU-VNTR analysis in Bangladesh

Bangladesh is among the high burden countries for tuberculosis (TB) and multidrug resistant TB (MDR-TB). As the genetic diversity and distinct phylogeographic distribution of Mycobacterium tuberculosis are responsible for regional differences in drug resistance, this cross sectional study was conducted to identify the circulating M. tuberculosis strains belonging to different lineages among pulmonary tuberculosis and, to investigate the contribution of distinct M. tuberculosis lineages to rifampicin resistant (RR) and rifampicin sensitive (RS) TB. A total of 40 RR and 20 RS isolates were enrolled in this study, all of which confirmed as M. tuberculosis by MPT 64 antigen detection. Furthermore, all isolates were genotyped by 24 loci Mycobacterial Interspersed Repetitive Units-Variable Number of Tandem Repeats (MIRU-VNTR), thus comprising the first study to employ this approach in Bangladesh. Beijing was the predominant lineage (26.8%) followed by EAI (23.2%), Delhi/CAS (16.1%), H37Rv (8.9%), Haarlem (7.1%), LAM (5.4%), Cameroon (3.6%) and a NEW-1 (1.8%). Four (7.1%) isolates remained as unidentified. Beijing strains were the significantly predominant (36.8%; p = 0.0135) among the RR isolates in comparison with other strains whereas EAI was the predominant (38.8%) lineage among RS isolates. Also, approximately 13% RR isolates showed genotypic resistance against fluoroquinolones by LPA and, hence, classed as pre-XDR TB albeit no specific lineage was found associated with these latter strains. A low transmission rate (10.5%) and high genetic diversity was detected in this setting with all the clustered strains herein identified belonging to the Beijing lineage. This study highlights 24 loci MIRU-VNTR analysis as a powerful tool for genotyping of Mycobacterium tuberculosis in this setting as it shows a high discriminatory index (0.81).

Molecular epidemiology of Mycobacterium tuberculosis isolates in Iran using spoligotyping

Spoligotyping can help assess the transmission of Mycobacterium tuberculosis strains. We aimed to study the genotyping of M. tuberculosis isolated from patients with tuberculosis from the west of Iran by spoligotyping. Forty-seven M. tuberculosis isolates were collected from the west of Iran. All samples were cultured on Löwenstein-Jensen medium incubated at 37°C for 8 weeks. Bacterial isolates were identified as M. tuberculosis using standard biochemical tests. Drug resistance patterns of M. tuberculosis to rifampicin and isoniazid were determined, and multidrug-resistant (MDR) strains were isolated. After DNA extraction, spoligotyping was performed. We found new spoligotypes 4162 and 4163, which correlated with atypical lineage. Atypical and unknown lineages also had correlations with the MDR tuberculosis rate (4%). The most prevalent spoligointernational types were orphan (34%), 2669 (23.4%) and 127 (14.8%) types. The most prevalent clades were Ural-2 (NEW-1) (25.53%) and atypical (23.40%) lineages. The predominant clade was Ural-2 (NEW-1) and an atypical lineage restricted to Iran. The rate of MDR was low. Knowledge of the circulating isolates in the west of Iran will help implement control programmes, so knowledge of the dynamic transmission of local isolates is crucial.

Mycobacterium tuberculosis complex in wildlife: Review of current applications of antemortem and postmortem diagnosis

Tuberculosis (TB) is a chronic inflammatory and zoonotic disease caused by Mycobacterium tuberculosis complex (MTBC) members, which affects various domestic animals, wildlife, and humans. Some wild animals serve as reservoir hosts in the transmission and epidemiology of the disease. Therefore, the monitoring and surveillance of both wild and domestic hosts are critical for prevention and control strategies. For TB diagnosis, the single intradermal tuberculin test or the single comparative intradermal tuberculin test, and the gamma-interferon test, which is regarded as an ancillary test, are used. Postmortem examination can identify granulomatous lesions compatible with a diagnosis of TB. In contrast, smears of the lesions can be stained for acid-fast bacilli, and samples of the affected organs can be subjected to histopathological analyses. Culture is the gold standard test for isolating mycobacterial bacilli because it has high sensitivity and specificity compared with other methods. Serology for antibody detection allows the testing of many samples simply, rapidly, and inexpensively, and the protocol can be standardized in different laboratories. Molecular biological analyses are also applicable to trace the epidemiology of the disease. In conclusion, reviewing the various techniques used in MTBC diagnosis can help establish guidelines for researchers when choosing a particular diagnostic method depending on the situation at hand, be it disease outbreaks in wildlife or for epidemiological studies. This is because a good understanding of various diagnostic techniques will aid in monitoring and managing emerging pandemic threats of infectious diseases from wildlife and also preventing the potential spread of zoonotic TB to livestock and humans. This review aimed to provide up-to-date information on different techniques used for diagnosing TB at the interfaces between wildlife, livestock, and humans.

Extensive Homoplasy but No Evidence of Convergent Evolution of Repeat Numbers at MIRU Loci in Modern Mycobacterium tuberculosis Lineages

More human deaths have been attributable to Mycobacterium tuberculosis than any other pathogen, and the epidemic is sustained by ongoing transmission. Various typing schemes have been developed to identify strain-specific differences and track transmission dynamics in affected communities, with recent introduction of whole genome sequencing providing the most accurate assessment. Mycobacterial interspersed repetitive unit (MIRU) typing is a family of variable number tandem repeat schemes that have been widely used to study the molecular epidemiology of M. tuberculosis. MIRU typing was used in most well-resourced settings to perform routine molecular epidemiology. Instances of MIRU homoplasy have been observed in comparison with sequence-based phylogenies, limiting its discriminatory value. A fundamental question is whether the observed homoplasy arises purely through stochastic processes, or whether there is evidence of natural selection. We compared repeat numbers at 24 MIRU loci with a whole genome sequence-based phylogeny of 245 isolates representing three modern M. tuberculosis lineages. This analysis demonstrated extensive homoplasy of repeat numbers, but did not detect any evidence of natural selection of repeat numbers, at least since the ancestral branching of the three modern lineages of M. tuberculosis. In addition, we observed good sensitivity but poor specificity and positive predictive values of MIRU-24 to detect clusters of recent transmission, as defined by whole-genome single nucleotide polymorphism analysis. These findings provide mechanistic insight, and support a transition away from VNTR-based typing toward sequence-based typing schemes for both research and public health purposes.

Phylogenetic diversity of Mycobacterium tuberculosis in two geographically distinct locations in Botswana - The Kopanyo Study

Mycobacterium tuberculosis complex (MTBC) is divided into several major phylogenetic lineages, with differential distribution globally. Using population-based data collected over a three year period, we performed 24-locus Mycobacterial Interspersed Repeat Unit - Variable Number Tandem Repeat (MIRU-VNTR) genotyping on all culture isolates from two districts of the country that differ in tuberculosis (TB) incidence (Gaborone, the capital, and Ghanzi in the Western Kalahari). The study objective was to characterize the molecular epidemiology of TB in these districts. Overall phylogenetic diversity mirrored that reported from neighboring Republic of South Africa, but differences in the two districts were marked. All four major lineages of M. tuberculosis were found in Gaborone, but only three of the four major lineages were found in Ghanzi. Strain diversity was lower in Ghanzi, with a large proportion (38%) of all isolates having an identical MIRU-VNTR result, compared to 6% of all isolates in Gaborone with the same MIRU-VNTR result. This study demonstrates localized differences in strain diversity by two districts in Botswana, and contributes to a growing characterization of MTBC diversity globally.

Exploring the usefulness of molecular epidemiology of tuberculosis in Africa: a systematic review

BACKGROUND

Tuberculosis (TB) is caused by Mycobacterium tuberculosis complex (MTBC) and remains a serious global public health threat, especially in resource-limited settings such as the African region. Recent developments in molecular epidemiology tools have significantly improved our understanding of TB transmission patterns and revealed the high genetic diversity of TB isolates across geographical entities in Africa. This study reports the results of a systematic review of current knowledge about MTBC strain diversity and geographical distribution in African regions.

METHODS

Search tools (PubMed, Embase, Popline, OVID and Africa Wide Information) were employed to identify the relevant literature about prevalence, strain diversity, and geographic distribution of MTBC infection in Africa.

RESULTS

A total of 59 articles from 739 citations met our inclusion criteria. Most articles reported about patients with presumptive pulmonary TB (73%), fewer reports were on retreatment and treatment failure cases (12%), and presumptive drug resistance cases (3%). Spoligotyping was the most used, alone in 21 studies and in parallel with either the Mycobacterial Interspersed Repetitive Units Variable Number of Tandem Repeats or the Restriction Fragment Length Polymorphism. Various TB lineages were observed across the African continent, with the originally European lineage 4 spotted in all countries studied.

CONCLUSION

TB molecular epidemiology tools have substantially improved our understanding of the MTBC circulating isolates, their evolution, and diversity in this highly endemic region of Africa. We found that only TB lineage 4 is present throughout all the continent and the clusters identified provides an extended insight into the disease transmission dynamics.

Most LAM Mycobacterium tuberculosis strains in Venezuela, but not SIT605, belong to the RDRio subfamily

Tuberculosis is a global public health problem that is resurgent in Venezuela, with 13 thousand estimated new cases in 2018. Strains of the Mycobacterium tuberculosis RD^Rio, subfamily belong to the Latín American Mediterranean (LAM) family and are a major cause of TB in Rio de Janeiro, Brazil. LAM strains predominate in Venezuela, where spoligotype SIT605 is common, but surprisingly not found elsewhere. We sought to assess the presence of RD^Rio strains in tuberculosis patients in different regions of Venezuela and determine whether SIT605 also belongs to the RD^Rio subfamily. Using spoligotyping and MIRU-VNTR 24 loci, we identified 86 clinical LAM and SIT605 isolates from the Venezuelan capital Caracas and several Venezuelan states. Region of difference deletion loci RD174 and RD^Rio, and also IS1561 were used to identify strains of the RD^Rio subfamily, while IS6110 at position 932,204 and the Ag85C¹⁰³ polymorphism were used to validate SIT 605 as a LAM family strain. We found that 69.8% of the isolates were RD^Río, including 94.3% of strains isolated in Caracas, 17.9% isolated in the state of Carabobo, the two strains analyzed from Delta Amacuro, and one each from Sucre, Apure and Aragua states. RD^Rio was in 100% of: SIT17 (LAM 2); SIT20 (LAM 1); SITs 93, 1694, 1696, 960, 1367 (LAM 5); and SITs 216 (LAM 9); but only 75% of SIT42 (LAM 9) strains. Thus, most of the LAM strains in Venezuela belong to the RD^Río subfamily. SIT 605 strains, although LAM, are not in the RD^Río subfamily.

Molecular epidemiology of drug resistant Mycobacterium tuberculosis in Africa: a systematic review

Background

The burden of drug resistant tuberculosis in Africa is largely driven by the emergence and spread of multidrug resistant (MDR) and extensively drug resistant (XDR) Mycobacterium tuberculosis strains. MDR-TB is defined as resistance to isoniazid and rifampicin, while XDR-TB is defined as MDR-TB with added resistance to any of the second line injectable drugs and any fluoroquinolone.

The highest burden of drug resistant TB is seen in countries further experiencing an HIV epidemic. The molecular mechanisms of drug resistance as well as the evolution of drug resistant TB strains have been widely studied using various genotyping tools. The study aimed to analyse the drug resistant lineages in circulation and transmission dynamics of these lineages in Africa by describing outbreaks, nosocomial transmission and migration. Viewed as a whole, this can give a better insight into the transmission dynamics of drug resistant TB in Africa.

Methods

A systematic review was performed on peer reviewed original research extracted from PubMed reporting on the lineages associated with drug resistant TB from African countries, and their association with outbreaks, nosocomial transmission and migration. The search terms “Tuberculosis AND drug resistance AND Africa AND (spoligotyping OR molecular epidemiology OR IS6110 OR MIRU OR DNA fingerprinting OR RFLP OR VNTR OR WGS)” were used to identify relevant articles reporting the molecular epidemiology of drug resistant TB in Africa.

Results

Diverse genotypes are associated with drug resistant TB in Africa, with variations in strain predominance within the continent. Lineage 4 predominates across Africa demonstrating the ability of “modern strains” to adapt and spread easily. Most studies under review reported primary drug resistance as the predominant type of transmission. Drug resistant TB strains are associated with community and nosocomial outbreaks involving MDR- and XDR-TB strains. The under-use of molecular epidemiological tools is of concern, resulting in gaps in knowledge of the transmission dynamics of drug resistant TB on the continent.

Conclusions

Genetic diversity of M. tuberculosis strains has been demonstrated across Africa implying that diverse genotypes are driving the epidemiology of drug resistant TB across the continent.

Polyclonal infection as a new scenario in Mycobacterium caprae epidemiology

Portugal is one of the European Union countries with an ongoing eradication program for bovine tuberculosis (TB), which does not include systematic goat testing. However, surveillance in small ruminants is increasingly important, since goat and sheep can harbour Mycobacterium caprae and be an infection source to cattle with impact in the success of bovine TB control. Furthermore, the information regarding the epidemiology and biology of M. caprae is quite limited comparing to the cognate bovine-adapted ecotype, M. bovis. In this work, we applied spoligotyping and MIRU-VNTR (Mycobacterial Interspersed Repetitive Units-Variable Number of Tandem Repeats) to M. caprae obtained between 2003 and 2014 from 55 animal hosts, including goat (n = 29), cattle (n = 21), sheep (n = 1) and wild boar (n = 4) from Portugal. The molecular analysis revealed a unique spoligotyping pattern (SB0157) and 24 MIRU types. Genotyping of serial M. caprae from herds with recurrent outbreaks enabled further discrimination of epidemiologically related isolates, supporting a clonal structure in Portugal and denoting the emergence of clonal diversity at the herd level, more apparent for MIRU4. Results suggest a founder effect and adaptive genotypic divergence, paving the way for sympatric speciation. Double allele findings at MIRU4 in over 20 % of infected animals indicates that co-infection and in vivo microevolution may be frequent in the goat-adapted ecotype. While polyclonal infection appears common in M. caprae epidemiology, the functional significance of subtle genotypic variations remains to be disclosed, namely at the interface with the host, to expand knowledge on the epidemiology and biology of this neglected ecotype.

Analysis of Mycobacterium tuberculosis genetic lineages circulating in Riga and Riga region, Latvia, isolated between 2008 and 2012

Although the number of new tuberculosis (TB) cases registered per year has decreased by 3-fold between 2001 and 2017 in Latvia, the TB incidence and rates of multidrug resistant TB in this Baltic country remain substantially higher than in most other European countries. Molecular typing methods of Mycobacterium tuberculosis (MTB) play an important role both in clinical studies of the disease and the epidemiological investigations, allowing to describe and characterize the pathogen's population structure and spread of particular genotypes. Aim of this study was to examine the prevalence of MTB lineages in Riga and Riga region of Latvia within a five-year period (2008-2012), and to evaluate the discriminatory power (DP) of spoligotyping, standard 24-locus MIRU-VNTR and IS6110-RFLP methods in this setting. The results showed that the main MTB spoligotype families were Beijing (25.3%) and LAM (24.3%), followed by T (22.1%), Ural (11.2%), Haarlem (6.6%) and X superfamily (3.4%). This distribution remained stable over the five consecutive years. 67.6% of MTB isolates were pan-susceptible, and 32.4% were resistant to any drug; multi-drug resistance was found in 5.8% of MTB strains, and 7.6% of MTB isolates were extensively drug-resistant. Drug resistance was associated with SIT1, SIT283 and SIT42 genotypes, while SIT1 and SIT42 were overrepresented among multi drug-resistant MTB strains. Overall, DP of spoligotyping method alone was 0.8953, while DP of both 24-locus MIRU-VNTR analysis and IS6110 RFLP was higher (DP = 0.9846 and 0.9927, respectively), mainly due to the improvement of the resolution for the Beijing strains. In conclusion, this work represents the first comprehensive molecular epidemiological description of TB in Latvia, highlighting the high genetic diversity of MTB strains circulating in Riga and Riga region. In combination with detailed epidemiological data this approach was helpful for the in-depth understanding of epidemiological processes in settings where the Next-Gen sequencing is not available as a routine method.

Accuracy of whole-genome sequencing to determine recent tuberculosis transmission: an 11-year population-based study in Hamburg, Germany

Controlling human-to-human tuberculosis (TB) transmission is key for achieving the targets of the End TB Strategy set by the World Health Organization (WHO) [1, 2]. Stopping TB transmission, in large cities especially, is a challenging top priority worldwide [3]. Metropolitan areas have higher TB case notification rates than the rest of a country, as they concentrate high-risk groups, such as homeless people, drug users and migrants often from (other) high TB incidence settings. Opportunities for transmission are amplified by population density and complex social interactions, regularly leading to large, temporally extended transmission networks [3]. Targeted interventions to interrupt transmission require the combination of effective genotyping of TB strains with enhanced epidemiological investigation. While classic IS6110 DNA fingerprinting and 24-locus MIRU–VNTR (mycobacterial interspersed repetitive units–variable number of tandem repeats) typing provide standardised and easily computable typing results with an online nomenclature system, several studies have now demonstrated that whole-genome sequencing (WGS) has a superior discriminatory power, allowing for an unparalleled resolution of outbreak strains [4–10]. However, predictivity of WGS for detecting transmission in metropolitan areas has not yet been quantified versus most deterministic references, i.e. tangible epidemiological links identified by ad hoc investigation, at extended time and population scales.

WGS typing with a five-SNP c:d:ut-off delineates recent transmission chains with highest accuracy and also provides high-resolution resistance patterns, thus enabling direct clinical benefits http://bit.ly/2Pk37Wo

Applicability of using 15 MIRU-VNTR loci for genotyping of Mycobacterium avium subsp. paratuberculosis from two cattle farms in Egypt

Mycobacterium avium subspecies paratuberculosis (MAP) is a notorious infectious agent that causes Johne's disease which leads to serious economic losses in cattle farms all over the world. The Lack of accurate epidemiological and molecular data is a major barrier to the implementation of disease control strategies. Basically, the tracing of infections requires rapid detection of the widely spreading genotypes with the ability to determine isolates from common and different sources. This study aimed to evaluate the applicability and discriminatory power of 15 mycobacterial interspersed repetitive unit (MIRU)-variable number tandem repeat (VNTR) loci of M. tuberculosis for MAP genotyping. Additionally, detection of the most efficient loci combinations for molecular epidemiological investigations of MAP isolates. The discriminatory capacity and applicability of 15 known loci [2 exact tandem repeat (ETR) loci, 6 MIRU loci, 4 Mtub loci, and 3 Queen's University of Belfast (QUB) group loci] were assessed using 26 isolates from two cattle herds (Holstein Frisian) in El buhaira and Giza Governorates at north of Egypt. The results proved the presence of 12 different genotypes. All the used loci gave Hunter-Gaston discrimination index of DI = 0.963 while the ten loci (Mtub04, MIRU10, QUB11b, MIRU26, QUB26, QUB4156, MIRU04, ETRC, Mtub30, and Mtub39) were highly discriminating with DI = 0.956. Moreover, the five loci (Mtub21, MIRU31, MIRU16, MIRU40, and ETRA) gave moderate discriminatory power with DI = 0.839. The MIRU31 locus expressed no polymorphism among strains. MIRU-VNTR typing generally proved applicability and high discriminatory power with DI = 0.963. The ten highly discriminating DI = 0.956 proved to be the most suitable for the first-line genotyping of MAP from cattle, with nearly similar resolving ability as all the 15 loci. MIRU-VNTR proved fastness, efficiency, and feasibility in genotyping of MAP from cattle in Egypt.

Mutation and Transmission Profiles of Second-Line Drug Resistance in Clinical Isolates of Drug-Resistant Mycobacterium tuberculosis From Hebei Province, China

The emergence of drug-resistant tuberculosis (TB) is involved in ineffective treatment of TB, especially multidrug resistant/extensively resistant TB (MDR/XDR-TB), leading to acquired resistance and transmission of drug-resistant strains. Second-line drugs (SLD), including both fluoroquinolones and injectable drugs, were commonly proved to be the effective drugs for treatment of drug-resistant TB. The purpose of this study was to investigate the prevalence of SLD-resistant strains and its specific mutations in drug-resistant Mycobacterium tuberculosis clinical isolates, and to acknowledge the transmission pattern of SLD resistance strains in Hebei. The genes gyrA, gyrB, rrs, eis promoter and tlyA of 257 drug-resistant clinical isolates were sequenced to identify mutations that could be responsible for resistance against fluoroquinolones and second-line injectable drugs. Each isolate was genotyped by Spoligotyping and 15-loci MIRU-VNTR. Our results indicated that 48.2% isolates were resistant to at least one of five SLD. Of them, 37.7% isolates were resistant to fluoroquinolones and 24.5% isolates were resistant to second-line injectable drugs. Mutations in genes gyrA, gyrB, rrs, eis promoter and tlyA were detected in 73 (75.3%), 7 (7.2%), 24 (38.1%), 5 (7.9%), and 3 (4.8%) isolates, respectively. The most prevalent mutations were the D94G (23.7%) in gyrA gene and the A1401G (33.3%) in rrs gene. A combination of gyrA, rrs and eis promoter can act as a valuable predicator for predicting XDR phenotype. These results highlight the development of rapid diagnosis are the effective manners for the control of SLD-TB or XDR-TB.

Exploring the "Latin American Mediterranean" family and the RDRio lineage in Mycobacterium tuberculosis isolates from Paraguay, Argentina and Venezuela

Background

The Latin American & Mediterranean (LAM) spoligotype family is one of the most successful genotype of Mycobacterium tuberculosis worldwide and particularly prevalent in South-America. Within this family, a sublineage named Region of Difference Rio (RD^Rio) was reported initially in Brazil and is characterized by a genomic deletion of about 26.3 kb. This lineage seems to show a specific adaptation to the Euro-Latin American population. In this context, we sought to evaluate the LAM family and the presence of the RD^Rio genotype in samples from three Latin American countries including Paraguay, Venezuela and Argentina. To detect LAM strains reliably we applied a typing scheme using spoligotyping, 12 loci MIRU-VNTR, the Ag85C¹⁰³ SNP and the regions of difference RD^Rio and RD174. IS6110-RFLP results were also used when available.

Results

Genotyping of 413 M. tuberculosis isolates from three Latin-American countries detected LAM (46%) and the ill-defined T clade (16%) as the most frequent families. The highest clustering rate was detected in the sample population from the city of Caracas in Venezuela. We observed considerable differences in the presence of the RD^Rio lineage, with high frequency in Caracas-Venezuela (55%) and low frequency in Buenos Aires-Argentina (11%) and Paraguay (10%). The molecular markers (RD174, Ag85C¹⁰³, MIRU02-MIRU40 signature) of the RD^Rio lineage were essentially confirmed. For the LAM family, the most polymorphic loci were MIRU40, MIRU31, MIRU10, MIRU26, MIRU16 and the least polymorphic MIRU24, MIRU20, MIRU04, MIRU23.

Conclusions

Our results suggest a differential adaptation of LAM-sublineages in neighboring populations and that RD^Rio strains spread regionally with different rates of distribution. The Ag85C SNP and RDs (RD174, RD^Rio) tested in this study can in fact facilitate molecular epidemiological studies of LAM strains in endemic settings and low-income countries.

Electronic supplementary material

The online version of this article (10.1186/s12866-019-1479-6) contains supplementary material, which is available to authorized users.

Highly successful treatment outcome of multidrug-resistant and genetic diversity of multidrug-resistant Mycobacterium tuberculosis strains in Rwanda

OBJECTIVE

To determine prevalent MDR-TB genotypes and describe treatment outcome and bacteriology conversion in MDR-TB patients.

METHODS

Review of laboratory records of 173 MDR-TB patients from all over Rwanda who initiated treatment under programmatic management of MDR-TB (PMDT) between 2014 and 2015. Fifty available archived isolates were genotyped by mycobacterial interspersed repetitive units - variable number of tandem repeats (MIRU-VNTR) genotyping.

RESULT

Of the 170 patients whose outcome was known, 114 (66.3%) were cured and 36 (21%) completed the treatment, giving a successful outcome (cured and completed) of 150 (87.3%) patients. Of 20 MDR-TB patients with unfavourable treatment outcome, 18 died, one failed and one defaulted/stopped treatment. Of the 18 patients who died, 11 (61%) were HIV-coinfected. The treatment outcome was successful for 93.9% among HIV negative and 81.8% among HIV-coinfected patients (P = 0.02). Sputum smear conversion occurred in 3, 46, 57 and 78 patients before 2, 3, 4 and 6 months, respectively, with median time of sputum smear and culture conversion at 3 months. The 44 MDR-TB isolates with MIRU-VNTR result, showed high genetic diversity with low clustering rate (9.09%) and Uganda II being the most prevalent sub-family lineage detected in 68.2% of isolates. Beijing family was the least common genotype detected (2.3%, 1 isolate).

CONCLUSION

The high success rates for MDR-TB treatment achieved in Rwanda were comparable to outcomes observed in resource-rich settings with HIV being an independent risk factor for poor treatment outcome. High genetic diversity and low clustering rate reported here suggest that reactivation of previous disease plays an important role in the transmission of MDR-TB in Rwanda.

Genetic diversity of Mycobacterium tuberculosis strains circulating in Botswana

BACKGROUND

Molecular typing of Mycobacterium tuberculosis (M.tb) isolates can inform Tuberculosis (TB) control programs on the relative proportion of transmission driving the TB epidemic. There is limited data on the M. tb genotypes that are circulating in Botswana. The aim of this study was to generate baseline data on the genetic diversity of M.tb isolates circulating in the country.

METHODS

A total of 461 M.tb isolates received at the Botswana National Tuberculosis Reference Laboratory between March 2012 and October 2013 were included in this study. Drug susceptibility testing was conducted using the BD BACTEC MGIT 960 System. M.tb strains were genotyped using spoligotyping and spoligotype patterns were compared with existing patterns in the SITVIT Web database. A subset of drug resistant isolates which formed spoligo clusters (n = 65) was additionally genotyped with 12-loci MIRU. Factors associated with drug resistance and clustering were evaluated using logistic regression.

RESULTS

Of the 461 isolates genotyped, 458 showed 108 distinct spoligotype patterns. The predominant M.tb lineages were Lineage 4 (81.9%), Lineage 2 (9%) and Lineage 1 (7.2%). The predominant spoligotype families within Lineage 4 were LAM (33%), S (14%), T (16%), X (16%). Three hundred and ninety-two (86%) isolates could be grouped into 44 clusters (2-46 isolates per cluster); giving a clustering rate of 76%. We identified 173 (37.8%) drug resistant isolates, 48 (10.5%) of these were multi-drug resistant. MIRU typing of the drug resistant isolates allowed grouping of 46 isolates into 14 clusters, giving a clustering rate of 49.2%. There was no association between age, sex, treatment category, region and clustering.

CONCLUSIONS

This study highlights the complexity of the TB epidemic in Botswana with multiple strains contributing to disease and provides baseline data on the population structure of M.tb strains in Botswana.

A first insight into the application of high discriminatory MIRU-VNTR typing using QIAxcel technology for genotyping Mycobacterium bovis isolated from the Delta area in Egypt

Mycobacterium bovis is a notorious infectious agent leading to serious economic losses for cattle farms worldwide. Analysis of the widely spreading genotypes is vital for tracing infections, understanding transmission dynamics, and controlling the cluster growth. This study aimed to evaluate the discrimination ability of 15 mycobacterial interspersed repetitive unit-variable number tandem repeats (MIRU-VNTR) loci and to assess the extremely efficient loci subset for molecular epidemiological investigations of M.bovis from farms in the Delta area of Egypt. The discriminating ability of MIRU-VNTR genotyping using 15 loci {2 exact tandem repeat (ETR) loci, 6 MIRU loci, 4 Mtub loci, and 3 Queen's University of Belfast (QUB) group loci)} were evaluated on 61 M.bovis isolates from cattle (Holstein Frisian) and buffaloes. The results indicate that there are 48 genotypes: 3 unique genotypes and 45 genotypes with shared similarities. Using the MIRU-VNTRplus database, M.bovis ID 7540/01 and ID 5346/02 were the nearest lineages to both groups. Six loci (MIRU10, QUB11b, QUB26, ETRA, Mtub30, and Mtub39) were highly discriminating, seven other loci (Mtub21, MIRU26, QUB4156, MIRU04 (ETRD), MIRU16, MIRU 40, and ETRC) gave moderate discriminatory power, and the last two loci (Mtub04 and MIRU31) were poorly discriminative. MIRU-VNTR typing generally proved efficacy and high discriminatory power, with a collective allele's diversification of 0.9641. Both the six highly discriminating (DI = 0.9492) and the seven moderately discriminating loci (DI = 0.9269) evidenced to be suitable for M.bovis first-step initial genotyping from cattle herds in Egypt. MIRU-VNTR is rapid and effective in the genotyping of M.bovis from cattle and buffaloes in Egypt.

Retrospective Analysis of False-Positive and Disputed Rifampin Resistance Xpert MTB/RIF Assay Results in Clinical Samples from a Referral Hospital in Hunan, China

Concerns about the specificity of the Xpert MTB/RIF (Xpert) assay have arisen, as false-positive errors in the determination of Mycobacterium tuberculosis complex (MTBC) infection and rifampin (RIF) resistance in clinical practice have been reported. Here, we investigated 33 cases where patients were determined to be RIF susceptible using the Bactec MGIT 960 (MGIT) culture system but RIF resistant using the Xpert assay. Isolates from two of these patients were found not to have any mutations in the rifampin resistance determining region (RRDR) region of rpoB and had good treatment outcomes with first-line antituberculosis (anti-TB) drugs. The remaining 31 patients included 5 new cases and 26 previously treated patients. A large number of well-documented disputed mutations, including Leu511Pro, Asp516Tyr, His526Asn, His526Leu, His526Cys, and Leu533Pro, were detected, and mutations, including a 508 to 509 deletion and His526Gly, were described here as disputed mutations for the first time. Twenty-one (81%) of the 26 previously treated patients had poor treatment outcomes, and isolates from 19 (90%) of these 21 patients were resistant to isoniazid (INH) as determined using the MGIT culture system. Twenty-seven of the 31 isolates with disputed rpoB mutations were phenotypically resistant to INH, 21 (78%) being predicted by GenoType MTBDRplus to have a high level of INH resistance. Most (77.4%) of the isolates with disputed mutations were of the Beijing lineage. These findings have implications for the interpretation of false-positive and disputed rifampin resistance Xpert MTB/RIF results in clinical samples and provide guidance on how clinicians should manage patients carrying isolates with disputed rpoB mutations.

A Clinical TB Detection Method Based on Molecular Typing Technique with Quality Control

The gold standard for diagnosing pulmonary Mycobacterium tuberculosis (TB) is the detection of tubercle bacillus in patient sputum samples. However, current methods either require long waiting times to culture the bacteria or have a risk of getting false-positive results due to cross-contamination. In this study, a method to detect tubercle bacillus based on the molecular typing technique is presented. This method can detect genetic units, variable number of tandem repeat (VNTR), which are the characteristic of tuberculosis (TB), and performs quality control using a mathematical model, ensuring the reliability of the results. Compared to other methods, the proposed method was able to process and diagnose a large volume of samples in a run time of six hours, with high sensitivity and specificity. Our method is also in the pipeline for implementation in clinical testing. Reliable and confirmed results are stored into a database, and these data are used to further refine the model. As the volume of data processed from reliable samples increases, the diagnostic power of the model improves. In addition to improving the quality control scheme, the collected data can be also used to support other TB research, such as that regarding the evolution of the tubercle bacillus.

Genotypic diversity of Mycobacterium tuberculosis isolates from North-Central Indian population

BACKGROUND

Different strains of Mycobacterium tuberculosis (MTB) are known to have different epidemiological and clinical characteristics. Some of them are widely distributed and associated with drug resistance, whereas others are locally predominated. Molecular epidemiological investigations have always been beneficial in identifying new strains and studying their transmission dynamics. Sahariya a primitive tribe of North Madhya Pradesh, India, has already been reported to have high prevalence of tuberculosis (TB) than their non-tribal neighbours. However, the information about MTB genotypes prevalent in Sahariya tribe and their non-tribal neighbours is not available.

METHODS

A total of 214 clinical isolates representing Sahariya tribe and non-tribes were analyzed by spoligotyping and MIRU-VNTR typing.

RESULTS

The EAI3_IND/SIT11 genotype was observed as major genotype in Sahariya tribe followed by CAS1_Delhi/SIT26 genotype. A 3.04 fold higher risk of getting TB with EAI3_IND/SIT11 genotype was observed in Sahariya as compared to the non-tribal population. The EAI_IND/SIT11 genotype also found to have more number of MDR-TB cases in Sahariya as well as true and possible transmission links. In Sahariya tribe, 3 clusters (6 isolates) reflected true transmission links, whereas 8 clusters consisted of 26 isolates revealed possible transmission links within the same geographical location or nearby houses.

CONCLUSION

The present study highlighted the predominance of EAI3_IND/SIT11 genotype in Sahariya tribe followed by CAS1_Delhi/SIT26 genotype. Combined approach of MIRU-VNTR typing and spoligotyping was observed more favourable in discrimination of MTB genotypes. Further, longitudinal studies using whole genome sequencing can provide more insights into genetic diversity, drug resistance and transmission dynamics of these prevalent genotypes.

Evaluation of the 15 and 24-loci MIRU-VNTR genotyping tools with spoligotyping in the identification of Mycobacterium tuberculosis strains and their genetic diversity in molecular epidemiology studies

BACKGROUND

The transmission dynamics of Mycobacterium tuberculosis (Mtb) using various genotyping tools has been studied globally and a particular tool for genotyping Mtb is the mycobacterial interspersed repetitive units-variable number tandem repeats (MIRU-VNTR). Tuberculosis (TB) remains an important public health problem worldwide and Sri Lanka being a country of tourist destination; because of major development projects undergoing, it has a high proportion of tourists and immigrants from Asia and Europe that are characterized with highest TB incidences and drug-resistant clinical isolates. Hence, in order to address the question of Mtb genetic diversity, we investigated the discriminatory power of both MIRU-VNTR typing of 15 and 24 loci with spoligotyping to differentiate Mtb isolates.

METHOD

Acid-fast bacilli positive sputum samples (n = 150) from first visit patients were collected. Decontamination of sputum and extraction of genomic DNA were carried out using standard techniques. The isolates were characterized by MIRU-VNTR for both the 15 and 24 loci and spoligotyping.

RESULTS

In our study population, MIRU-VNTR 15 and 24 loci did not show a significant difference among the identified M. tuberculosis strains. However, MIRU 24 loci yielded an additional strain LAM, which is of T1 origin. 15 loci strain grouping had more clusters of strains grouped together while 24 loci differentiated the same cluster of strains into distinct strain types.

CONCLUSION

We conclude that the use of 15-locus MIRU-VNTR typing is sufficient for a first-line epidemiological study to genotype M. tuberculosis, but the additional discriminatory power of 24 loci MIRU-VNTR has been able to differentiate samples within highly homologous groups.

Modulatory effects of verapamil in rifampicin activity against Mycobacterium tuberculosis

AIM

To evaluate modulatory effect of verapamil (VP) in rifampicin (RIF) activity and its effect in efflux pumps (EPs) transcript levels in Mycobacterium tuberculosis.

MATERIALS & METHODS

RIF and VP minimal inhibitory concentration, combinatory effect and detection of mutations were determined in 16 isolates. EPs transcript levels were determined in four isolates by real-time PCR after exposure to drugs.

RESULTS

VP showed good combinatory effect among RIF-resistant isolates. This effect was also observed in the relative transcript levels of EPs, mainly after 72 h of exposure, depending on the EP gene, genotype and the resistance profile of the isolate.

CONCLUSION

Additional regulatory mechanisms in the EP activities, as well as, interactions with other drug-specific resistance mechanisms need further investigation in M. tuberculosis.

Integrating standardized whole genome sequence analysis with a global Mycobacterium tuberculosis antibiotic resistance knowledgebase

Drug-resistant tuberculosis poses a persistent public health threat. The ReSeqTB platform is a collaborative, curated knowledgebase, designed to standardize and aggregate global Mycobacterium tuberculosis complex (MTBC) variant data from whole genome sequencing (WGS) with phenotypic drug susceptibility testing (DST) and clinical data. We developed a unified analysis variant pipeline (UVP) ( https://github.com/CPTR-ReSeqTB/UVP ) to identify variants and assign lineage from MTBC sequence data. Stringent thresholds and quality control measures were incorporated in this open source tool. The pipeline was validated using a well-characterized dataset of 90 diverse MTBC isolates with conventional DST and DNA Sanger sequencing data. The UVP exhibited 98.9% agreement with the variants identified using Sanger sequencing and was 100% concordant with conventional methods of assigning lineage. We analyzed 4636 publicly available MTBC isolates in the ReSeqTB platform representing all seven major MTBC lineages. The variants detected have an above 94% accuracy of predicting drug based on the accompanying DST results in the platform. The aggregation of variants over time in the platform will establish confidence-graded mutations statistically associated with phenotypic drug resistance. These tools serve as critical reference standards for future molecular diagnostic assay developers, researchers, public health agencies and clinicians working towards the control of drug-resistant tuberculosis.