Evaluation of spoligotyping, SNPs and customised MIRU-VNTR combination for genotyping Mycobacterium tuberculosis clinical isolates in Madagascar

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.

Development, Evaluation, and Implementation of a House-Made Targeted Next-Generation Sequencing Spoligotyping in a French Laboratory

Epidemiological studies investigating transmission chains of tuberculosis are undertaken worldwide to tackle its spread. CRISPR locus diversity, called spoligotyping, is a widely used genotyping assay for Mycobacterium tuberculosis complex (MTBC) characterization. Herein, we developed a house-made targeted next-generation sequencing (tNGS) spoligotyping, and compared its outputs with those of membrane-based spoligotyping. A total of 144 clinical MTBC strains were retrospectively selected to be representative of the local epidemiology. Data analysis of a training set allowed for the setting of “presence”/“absence” thresholds for each spacer to maximize the sensibility and specificity related to the membrane-based spoligotyping. The thresholds above, in which the spacer was considered present, were 50 read per millions for spacers 10 and 14, 20,000 for spacers 20, 21, and 31, and 1000 for the other spacers. The confirmation of these thresholds was performed using a validation set. The overall agreement on the training and validation sets was 97.5% and 93.8%, respectively. The discrepancies concerned six strains: Two for spacer 14, two for spacer 31, and two for spacer 32. The tNGS spoligotyping, whose thresholds were finely-tuned during a careful bioinformatics pipeline development process, appears be a technique that is reliable, inexpensive, free of handling errors, and automatable through automatic transfer into the laboratory computer system.

The shortcut of mycobacterial interspersed repetitive unit-variable number tandem repeat typing for Mycobacterium tuberculosis differentiation

The 24-loci mycobacterial interspersed repetitive unit-variable number tandem repeat (MIRU-VNTR) genotyping has been used as an international standard method for Mycobacterium tuberculosis (Mtb) genotyping. However, different optimized VNTR loci sets for improving the discrimination of specific Mtb genotypes have been proposed. In this regard, we investigated the efficacy of accumulation of the percentage differences (APDs) compared with the least absolute shrinkage and selection operator (LASSO) regression strategy to identify a customized genotype-specific VNTR loci set which provides a resolution comparable to 24-loci MIRU-VNTR in divergent Mtb populations. We utilized Spoligotyping and 24-loci MIRU-VNTR typing for genotyping 306 Mtb isolates. The APD and LASSO regression approaches were used to identify a customized VNTR set in our studied isolates. Besides, the Hunter-Gaston discriminatory index (HGDI), sensitivity, and specificity of each selected loci set were calculated based on both strategies. The selected loci based on LASSO regression compared with APD-based loci showed a better discriminatory power for identifying all studied genotypes except for T genotype, which APD-based loci showed promising discriminative power. Our findings suggested the LASSO regression rather than the APD approach is more effective in the determination of possible discriminative VNTR loci set to precise discrimination of our studied Mtb population and may be beneficial to be used in finding reduced number loci sets in other Mtb genotypes or sublineages. Moreover, we proposed customized genotype-specific MIRU-VNTR loci sets based on the LASSO regression and APD approaches for precise Mtb strains identification. As the proposed VNTR sets offered a comparable discriminatory power to the standard 24 MIRU-VNTR loci set could be promising alternatives to the standard genotyping for using in resource-limited settings.

Consistency of Mycobacterium tuberculosis Complex Spoligotyping between the Membrane-Based Method and In Silico Approach

To tackle the spread of tuberculosis (TB), epidemiological studies are undertaken worldwide to investigate TB transmission chains. Clustered regulatory interspaced short palindromic repeats (CRISPR) locus diversity, also called spoligotyping, is a widely used genotyping assay for the characterization of Mycobacterium tuberculosis complex (MTBC). We compared herein the spoligotyping of MTBC clinical isolates using a membrane-based method (following an initial PCR step) and whole-genome sequencing (WGS)-based method (i.e., in silico spoligotyping). All MTBC strains isolated at the Lyon University Hospital, France, between November 2016 and December 2020 were included (n = 597). Spoligotyping profiles were also used for species identification among the MTBC. Outputs of both methods were analyzed, and discrepant results were investigated thanks to CRISPRbuilder-TB. The overall agreement was 85.7%. Spacer discrepancies observed between the methods were due to the insertion of IS6110 within the direct repeat (DR) sequence upstream or downstream of spacers, mutated DR sequences, or truncated spacers. Discrepancies did not impact species identification. Although spoligotyping-based species identification was inconclusive for 29 isolates, SNP-based phylogeny conducted after WGS allowed the identification of 23 M. tuberculosis (Mtb), 2 M. canettii, and 4 mixed MTBC infections. WGS yielded very few discrepancies compared to membrane-based spoligotyping. Overall agreement was significantly improved (92.4%) by the CRISPR locus reconstruction using CRISPRbuilder-TB for the MTBC isolates with the shared international type 53 in silico spoligotyping. A smooth transition from the membrane-based to the in silico-based genotyping of M. tuberculosis isolates is, therefore, possible for TB diagnosis and epidemiologic survey. IMPORTANCE Whole-genome sequencing (WGS) has profoundly transformed the perspectives of tuberculosis (TB) diagnosis, providing a better discriminatory power to determine relatedness between Mycobacterium tuberculosis complex (MTBC) isolates. Previous genotyping approaches, such as spoligotyping consisting of an initial PCR step followed by reverse dot hybridization, are currently being replaced by WGS. Several pipelines have been developed to extract a spoligotype from WGS data (in silico spoligotyping) allowing for the continuity of MTBC molecular surveys before and after WGS implementation. The present study found very good overall agreement between hybridization to membrane-based spoligotyping and in silico spoligotyping, indicating the possibility of a smooth transition from the traditional to the in silico-based genotyping of MTBC isolates for TB diagnosis and epidemiological survey.

Developing customized stepwise MIRU-VNTR typing for tuberculosis surveillance in Georgia

Introduction

Mycobacterial Interspersed Repetitive Units–Variable Tandem Repeats (MIRU-VNTR) typing has been widely used for molecular epidemiological studies of tuberculosis (TB). However, genotyping tools for Mycobacterium tuberculosis (Mtb) may be limiting in some settings due to high cost and workload. In this study developed a customized stepwise MIRU-VNTR typing that prioritizes high discriminatory loci and validated this method using penitentiary system cohort in the country of Georgia.

Methods

We used a previously generated MIRU-VNTR dataset from recurrent TB cases (32 cases) in Georgia and a new dataset of TB cases from the penitentiary system (102 cases) recruited from 2014 to 2015. A Hunter-Gaston Discriminatory Index (HGDI) was calculated utilizing a 24 standard loci panel, to select high discriminatory power loci, subsequently defined as the customized Georgia-specific set of loci for initial typing. The remaining loci were scored and hierarchically grouped for second and third step typing of the cohort. We then compared the processing time and costs of the customized stepwise method to the standard 24-loci method.

Results

For the customized Georgia-specific set that was used for initial typing, 10 loci were selected with a minimum value of 0.32 to the highest HGDI score locus. Customized 10 loci (step 1) typing of 102 Mtb patient isolates revealed 35.7% clustered cases. This proportion was reduced to 19.5% after hierarchical application of 2nd and 3rd step typing with the corresponding groups of loci. Our customized stepwise MIRU-VNTR genotyping approach reduced the quantity of samples to be typed and therefore overall processing time and costs by 42.6% each.

Conclusion

Our study shows that our customized stepwise MIRU-VNTR typing approach is a valid alternative of standard MIRI-VNTR typing panels for molecular epidemiological investigation in Georgia that saves time, workload and costs. Similar approaches could be developed for other settings.

A Comprehensive Evaluation of GeneLEAD VIII DNA Platform Combined to Deeplex Myc-TB® Assay to Detect in 8 Days Drug Resistance to 13 Antituberculous Drugs and Transmission of Mycobacterium tuberculosis Complex Directly From Clinical Samples

The GeneLEAD VIII (Diagenode, Belgium) is a new, fully automated, sample-to-result precision instrument for the extraction of DNA and PCR detection of Mycobacterium tuberculosis complex (MTBC) directly from clinical samples. The Deeplex Myc-TB^® assay (Genoscreen, France) is a diagnostic kit based on the deep sequencing of a 24-plexed amplicon mix allowing simultaneously the detection of resistance to 13 antituberculous (antiTB) drugs and the determination of spoligotype. We evaluated the performance of a strategy combining the both mentioned tools to detect directly from clinical samples, in 8 days, MTBC and its resistance to 13 antiTB drugs, and identify potential transmission of strains from patient-to-patient. Using this approach, we screened 112 clinical samples (65 smear-negative) and 94 MTBC cultured strains. The sensitivity and the specificity of the GeneLEAD/Deeplex Myc-TB approach for MTBC detection were 79.3% and 100%, respectively. One hundred forty successful Deeplex Myc-TB results were obtained for 46 clinical samples and 94 strains, a total of 85.4% of which had a Deeplex Myc-TB susceptibility and resistance prediction consistent with phenotypic drug susceptibility testing (DST). Importantly, the Deeplex Myc-TB assay was able to detect 100% of the multidrug-resistant (MDR) MTBC tested. The lowest concordance rates were for pyrazinamide, ethambutol, streptomycin, and ethionamide (84.5%, 81.5%, 73%, and 55%, respectively) for which the determination of susceptibility or resistance is generally difficult with current tools. One of the main difficulties of Deeplex Myc-TB is to interpret the non-synonymous uncharacterized variants that can represent up to 30% of the detected single nucleotide variants. We observed a good level of concordance between Deeplex Myc-TB-spoligotyping and MIRU-VNTR despite a lower discriminatory power for spoligotyping. The median time to obtain complete results from clinical samples was 8 days (IQR 7–13) provided a high-throughput NGS sequencing platform was available. Our results highlight that the GeneLEAD/Deeplex Myc-TB approach could be a breakthrough in rapid diagnosis of MDR TB in routine practice.

Rapid typing of carbapenem-resistant Acinetobacter baumannii and Acinetobacter nosocomialis by multiplex Pan- and OXA-PCR assays

Introduction. Outbreaks of carbapenem-resistant A. baumannii and A. nosocomialis have occurred worldwide in healthcare settings. Rapid and reliable molecular typing of bacterial isolates is vital for the effective surveillance of institutional outbreaks. The Pan-PCR and OXA-PCR assays are two multiplex PCR-based assays for the molecular typing of Acinetobacter species.Gap statement. However, few studies have investigated the discriminatory power of two multiplex PCR assays in in the genotyping of Acinetobacter species.Aim. We aimed to evaluate the efficacies of the Pan-PCR and OXA-PCR assays for molecular typing of A. baumannii and A. nosocomialis.Methodology. A total of 105 carbapenem-resistant A. baumannii isolates (CRABs) and 93 carbapenem-resistant A. nosocomialis isolates (CRANs) obtained from blood cultures were used for molecular typing by the Pan-PCR and OXA-PCR assays and two multilocus sequence typing (MLST) schemes.Results. The isolates were individually divided into 12 and 21 different sequence types via the Pasteur and Oxford MLST schemes, respectively. Additionally, these isolates were distinguished into 18 different types by the Pan-PCR and OXA-PCR assays. The results of the Pan-PCR and OXA-PCR assays distinguished CRABs and CRANs with a sensitivity of 98.13 % and a specificity of 100 %.Conclusion. The Pan-PCR and OXA-PCR assays are promising alternative methods for rapid molecular typing of CRABs and CRANs in a routine laboratory setting.

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.

Enhancing the differentiation of specific genotypes in Mycobacterium tuberculosis population

Today, significant attention is directed towards the global lineages and sublineages of Mycobacterium tuberculosis (Mtb). NEW-1 (SIT 127) and CAS1-Delhi (SIT 26) strains are recognized as growing and circulating Mtb genotypes, especially in Asian countries. It is crucial to develop or enhance Mtb genotyping methods for a more accurate and simple differentiation of these families. We used 24-loci mycobacterial interspersed repetitive unit-variable number tandem repeat (MIRU-VNTR) typing for genotyping 217 Mtb isolates. To select the optimal MIRU-VNTR loci, we calculated the Hunter-Gaston discriminatory index (HGDI), allelic diversity, and accumulation of percentage differences (APDs) between the strains among different groups of genotypes (NEW-1 and non-NEW-1; CAS1-Delhi and non-CAS). Finally, the minimum spanning tree was constructed for clustering analysis. In the NEW-1 population, loci with APD > 60% were found to have a high discriminatory power. VNTR loci with APD > 50% showed high discrimination power for the CAS population. Our findings suggest that APDs, which are valuable for the selection of VNTR loci sets, may improve the discriminatory power of MIRU-VNTR typing for identification of Mtb genotypes in specific regions.

Introducing the Best Six Loci in Mycobacterial Interspersed Repetitive Unit-Variable-Number Tandem Repeat (MIRU-VNTR) Typing for Mycobacterium Tuberculosis Genotyping

BACKGROUND

Tuberculosis (TB) still remains endemic worldwide making epidemiological studies essential to mitigating efforts implicated in identifying its source, controlling, and preventing the spread of dangerous strains amongst humans such as Mycobacterium tuberculosis (Mtb).

METHODS

In this study, we sought to determine the 6 Mycobacterial Interspersed Repetitive Unit-Variable-Number Tandem Repeat (MIRU-VNTR) loci with high discriminatory powers for Mtb genotyping as well as the loci with the highest and the lowest discriminatory powers for MIRU-VNTR. To conduct our search, we used several databases such as science direct, Embase (Elsevier), Web of Science, Scopus and Medline via PubMed. Searches were performed using key words including: Mycobacterium tuberculosis, MIRU-VNTR, Allele diversity, Genetic diversity and human patient. Finally, 56 articles were selected after filtering out titles, abstracts and full texts.

RESULTS

Loci with high discriminatory powers included MIRU10 and MIRU26, while MIRU2, MIRU20, MIRU24 and ETRD had poor discriminatory powers. According to previous data in the literature, the loci MIRU10, MIRU26, MIRU40, QUB 26, QUB 11b and Mtub21 have high discriminatory powers.

CONCLUSION

Therefore, these loci recommended for genotyping Mtb to save time and cost and to ensure the production of reliable results.

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.

Evaluation of a gene-by-gene approach for prospective whole-genome sequencing-based surveillance of multidrug resistant Mycobacterium tuberculosis

Whole-genome sequencing (WGS) offers unprecedented resolution for tracking Mycobacterium tuberculosis transmission and antibiotic-resistance spread. Still, the establishment of standardized WGS-based pipelines and the definition of epidemiological clusters based on genetic relatedness are under discussion. We aimed to implement a dynamic gene-by-gene approach, fully relying on freely available software, for prospective WGS-based tuberculosis surveillance, demonstrating its application for detecting transmission chains by retrospectively analysing all M/XDR strains isolated in 2013-2017 in Portugal. We observed a good correlation between genetic relatedness and epidemiological links, with strongly epilinked clusters displaying mean pairwise allele differences (AD) always below 0.3% (ratio of mean AD over the total number of shared loci between same-cluster strains). This data parallels the genetic distances acquired by the core-SNV analysis, while providing higher resolution and epidemiological concordance than MIRU-VNTR genotyping. The dynamic analysis of strain sub-sets (i.e., increasing the number of shared loci within each sub-set) also strengthens the confidence in detecting epilinked clusters. This gene-by-gene strategy also offers several practical benefits (e.g., reliance on freely-available software, scalability and low computational requirements) that further consolidated its suitability for a timely and robust prospective WGS-based laboratory surveillance of M/XDR-TB cases.

Tuberculosis in the 21th century: Current status of diagnostic methods

Tuberculosis is an infectious bacterial disease with a high mortality rate worldwide constituting a serious public health problem. The diagnostic methods commonly used by health professionals are slow and expensive and the results may take about sixty days which will cause a delay in administrating the most proper treatment to the patient, as well as increase health care costs and infection transmission possibility. Patients infected simultaneously with human immunodeficiency virus and Mycobacterium tuberculosis are a constant and worrying challenge for the scientific community which will research and develop new methods of diagnosis, new drugs and new therapies. Nowadays there are new tuberculosis diagnosis methods and some of which are already in clinical trial phases. These methods have high sensitivity, but do not replace the microbiological examination for isolation and culture of Mycobacterium spp. However, in clinical practice, microbiological, imaging, clinical and epidemiological data integration provide the best diagnosis and treatment possible. Consequently, throughout this paper, the different methods of diagnosis of human tuberculosis with its advantages and disadvantages will be covered, describing new omics and ultra-fast methods to increase knowledge and obtain a rapid diagnosis of tuberculosis.

Genotypic diversity of Mycobacterium tuberculosis in Buenos Aires, Argentina

Buenos Aires is an overpopulated port city historically inhabited by people of European descent. Together with its broader metropolitan area, the city exhibits medium tuberculosis rates, and receives migrants, mainly from tuberculosis highly endemic areas of Argentina and neighboring countries. This work was aimed to gain insight into the Mycobacterium tuberculosis population structure in two suburban districts of Buenos Aires which are illustrative of the overall situation of tuberculosis in Argentina. The Lineage 4 Euro-American accounted for >99% of the 816 isolates analyzed (one per patient). Frequencies of spoligotype families were T 35.9%, LAM 33.2%, Haarlem 19.5%, S 3.2%, X 1.5%, Ural 0.7%, BOV 0.2%, Beijing 0.2%, and Cameroon 0.2%. Unknown signatures accounted for 5.3% isolates. Of 55 spoligotypes not matching any extant shared international type (SIT) in SITVIT database, 22 fitted into 15 newly-issued SITs. Certain autochthonous South American genotypes were found to be actively evolving. LAM3, which is wild type for RD^rio, was the predominant LAM subfamily in both districts and the RD^rio signature was rare among autochthonous, newly created, SITs and orphan patterns. Two genotypes that are rarely observed in neighboring countries ̶ SIT2/H2 and SIT159/T1 Tuscany ̶ were conspicuously represented in Argentina. The infrequent Beijing patterns belonged to Peruvian patients. We conclude that the genotype diversity observed reflects the influence of the Hispanic colonization and more recent immigration waves from Mediterranean and neighboring countries. Unlike in Brazil, the RD^rio type does not play a major role in the tuberculosis epidemic in Buenos Aires.