Novel methods included in SpolLineages tool for fast and precise prediction of Mycobacterium tuberculosis complex spoligotype families

Whole-genome sequencing-based genetic diversity, transmission dynamics, and drug-resistant mutations in Mycobacterium tuberculosis isolated from extrapulmonary tuberculosis patients in western Ethiopia

Background

Extrapulmonary tuberculosis (EPTB) refers to a form of Tuberculosis (TB) where the infection occurs outside the lungs. Despite EPTB being a devastating disease of public health concern, it is frequently overlooked as a public health problem. This study aimed to investigate genetic diversity, identify drug-resistance mutations, and trace ongoing transmission chains.

Methods

A cross-sectional study was undertaken on individuals with EPTB in western Ethiopia. In this study, whole-genome sequencing (WGS) was employed to analyze Mycobacterium tuberculosis (MTB) samples obtained from EPTB patients. Out of the 96 genomes initially sequenced, 89 met the required quality standards for genetic diversity, and drug-resistant mutations analysis. The data were processed using robust bioinformatics tools.

Results

Our analysis reveals that the majority (87.64%) of the isolates can be attributed to Lineage-4 (L4), with L4.6.3 and L4.2.2.2 emerging as the predominant sub-lineages, constituting 34.62% and 26.92%, respectively. The overall clustering rate and recent transmission index (RTI) were 30 and 17.24%, respectively. Notably, 7.87% of the isolates demonstrated resistance to at least one anti-TB drug, although multi-drug resistance (MDR) was observed in only 1.12% of the isolates.

Conclusions

The genetic diversity of MTBC strains in western Ethiopia was found to have low inter-lineage diversity, with L4 predominating and exhibiting high intra-lineage diversity. The notably high clustering rate in the region implies a pressing need for enhanced TB infection control measures to effectively disrupt the transmission chain. It's noteworthy that 68.75% of resistance-conferring mutations went undetected by both GeneXpert MTB/RIF and the line probe assay (LPA) in western Ethiopia. The identification of resistance mutations undetected by both GeneXpert and LPA, along with the detection of mixed infections through WGS, emphasizes the value of adopting WGS as a high-resolution approach for TB diagnosis and molecular epidemiological surveillance.

An overview of next generation sequencing strategies and genomics tools used for tuberculosis research

Tuberculosis (TB) is a grave public health concern and is considered the foremost contributor to human mortality resulting from infectious disease. Due to the stringent clonality and extremely restricted genomic diversity, conventional methods prove inefficient for in-depth exploration of minor genomic variations and the evolutionary dynamics operating in Mycobacterium tuberculosis (M.tb) populations. Until now, the majority of reviews have primarily focused on delineating the application of whole-genome sequencing (WGS) in predicting antibiotic resistant genes, surveillance of drug resistance strains, and M.tb lineage classifications. Despite the growing use of next generation sequencing (NGS) and WGS analysis in TB research, there are limited studies that provide a comprehensive summary of there role in studying macroevolution, minor genetic variations, assessing mixed TB infections, and tracking transmission networks at an individual level. This highlights the need for systematic effort to fully explore the potential of WGS and its associated tools in advancing our understanding of TB epidemiology and disease transmission. We delve into the recent bioinformatics pipelines and NGS strategies that leverage various genetic features and simultaneous exploration of host-pathogen protein expression profile to decipher the genetic heterogeneity and host-pathogen interaction dynamics of the M.tb infections. This review highlights the potential benefits and limitations of NGS and bioinformatics tools and discusses their role in TB detection and epidemiology. Overall, this review could be a valuable resource for researchers and clinicians interested in NGS-based approaches in TB research.

Genetic diversity of Mycobacterium tuberculosis strains isolated from spiritual holy water site attendees in Northwest Ethiopia. A cross-sectional study

Background

The genetic diversity of Mycobacterium tuberculosis complex (MTBC) strains was characterized among isolates from individuals with pulmonary tuberculosis (PTB) symptoms attended holy water sites (HWSs) in the Amhara region, Ethiopia.

Methods

A cross-sectional study was done from June 2019 to March 2020 to describe the genetic diversity and drug-resistance profiles of MTBC isolates. Sputum specimens were collected and cultured in the Löwenstein-Jensen culture medium. Line Probe Assay, MTBDRplus VER 2.0, and MTBDRsl VER 2.0 were used to detect first-and second-line anti-TB drug-resistance patterns. A spoligotyping technique was utilized to characterize the genetic diversity. Statistical analysis was performed using STATA 15.

Results

Of 560 PTB-symptomatic participants, 122 (21.8%) were culture-positive cases. Spoligotyping of 116 isolates revealed diverse MTBC sublineages, with four major lineages: Euro-American (EA) (Lineage 4), East-African-Indian (EAI) (Lineage 3), Ethiopian (ETH) (Lineage 7), East Asian (EA) (Lineage 2). The majority (96.6%) of the isolates were EA (lineage 4) and EAI, with proportions of 54.3% and 42.2%, respectively. A total of 31 spoligotype patterns were identified, 26 of which were documented in the SITVIT2 database. Of these, there were 15 unique spoligotypes, while eleven were grouped with 2-17 isolates. SIT149/T3-ETH (n = 17), SIT26/CAS1-DELHI (n = 16), SIT25/CAS1-DELHI (n = 12), and SIT52/T2 (n = 11) spoligotypes were predominant. A rare spoligotype pattern: SIT41/Turkey and SIT1/Beijing, has also been identified in North Shewa. The overall clustering rate of sub-lineages with known SIT was 76.4%.Of the 122 culture-positive isolates tested, 16.4% were resistant to rifampicin (RIF) and/or isoniazid (INH). Multidrug-resistant TB (MDR-TB) was detected in 12.3% of isolates, five of which were fluoroquinolones (FLQs) resistant. SIT149/T3-ETH and SIT21/CAS1-KILI sublineages showed a higher proportion of drug resistance.

Conclusions

Diverse MTBC spoligotypes were identified, with the T and CAS families and EA (lineage 4) predominating. A high prevalence of drug-resistant TB, with SIT149/T3-ETH and CAS1-KILI sublineages comprising a greater share, was observed. A study with large sample size and a sequencing method with stronger discriminatory power is warranted to understand better the genetic diversity of circulating MTBC in this cohort of study, which would help to adopt targeted interventions.

Population structure and spatial distribution of Mycobacterium tuberculosis in Ethiopia

Ethiopia is one of the countries with a high tuberculosis (TB) burden, yet little is known about the spatial distribution of Mycobacterium tuberculosis (Mtb) lineages. This study identifies the spoligotyping of 1735 archived Mtb isolates from the National Drug Resistance Survey, collected between November 2011 and June 2013, to investigate Mtb population structure and spatial distribution. Spoligotype International Types (SITs) and lineages were retrieved from online databases. The distribution of lineages was evaluated using Fisher's exact test and logistic regression models. The Global Moran's Index and Getis-Ord Gi statistic were utilized to identify hotspot areas. Our results showed that spoligotypes could be interpreted and led to 4 lineages and 283 spoligotype patterns in 91% of the isolates, including 4% of those with multidrug/rifampicin resistance (MDR/RR) TB. The identified Mtb lineages were lineage 1 (1.8%), lineage 3 (25.9%), lineage 4 (70.6%) and lineage 7 (1.6%). The proportion of lineages 3 and 4 varied by regions, with lineage 3 being significantly greater than lineage 4 in reports from Gambella (AOR = 4.37, P < 0.001) and Tigray (AOR = 3.44, P = 0.001) and lineage 4 being significantly higher in Southern Nations Nationalities and Peoples Region (AOR = 1.97, P = 0.026) than lineage 3. Hotspots for lineage 1 were located in eastern Ethiopia, while a lineage 7 hotspot was identified in northern and western Ethiopia. The five prevalent spoligotypes, which were SIT149, SIT53, SIT25, SIT37 and SIT26 account for 42.8% of all isolates under investigation, while SIT149, SIT53 and SIT21 account for 52-57.8% of drug-resistant TB cases. TB and drug resistant TB are mainly caused by lineages 3 and 4, and significant proportions of the prevalent spoligotypes also influence drug-resistant TB and the total TB burden. Regional variations in lineages may result from both local and cross-border spread.

Genetic diversity of Mycobacterium tuberculosis isolates from the central, eastern and southeastern Ethiopia

Introduction

The population structure of Mycobacterium tuberculosis complex (MTBC) in Ethiopia is diverse but dominated by Euro-American (Lineage 4) and East-African-Indian (Lineage 3) lineages. The objective of this study was to describe the genetic diversity of MTBC isolates in Central, Eastern and Southeastern Ethiopia.

Methods

A total of 223 MTBC culture isolates obtained from patients referred to Adama and Harar TB reference laboratories were spoligotyped. Demographic and clinical characteristics were collected.

Results

Six major lineages: Euro-American (Lineage 4), East-African-Indian (Lineage 3), East Asian (Lineage 2), Indo-Oceanic (Lineage 1), Mycobacterium africanum (Lineage 5 and Lineage 6) and Ethiopian (Lineage 7) were identified. The majority (94.6 %) of the isolates were Euro-American and East-African-Indian, with proportions of 75.3 % and 19.3 %, respectively. Overall, 77 different spoligotype patterns were identified of which 42 were registered in the SITVIT2 database. Of these, 27 spoligotypes were unique, while 15 were clustered with 2-49 isolates. SIT149/T3_ETH (n = 49), SIT53/T1 (n = 33), SIT21/CAS1_Kili (n = 24) and SIT41/Turkey (n = 11) were the dominant spoligotypes. A rare Beijing spoligotype pattern, SIT541, has also been identified in Eastern Ethiopia. The overall clustering rate of sub-lineages with known SIT was 71.3 %. Age group (25-34) was significantly associated with clustering.

Conclusion

We found a heterogeneous population structure of MTBC dominated by T and CAS families, and the Euro-American lineage. The identification of the Beijing strain, particularly the rare SIT541 spoligotype in Eastern Ethiopia, warrants a heightened surveillance plan, as little is known about this genotype. A large-scale investigation utilizing a tool with superior discriminatory power, such as whole genome sequencing, is necessary to gain a thorough understanding of the genetic diversity of MTBC in the nation, which would help direct the overall control efforts.

simpiTB - a pipeline designed to extract meaningful information from whole genome sequencing data of Mycobacterium tuberculosis complex, allows to combine genomic, phylogenetic and clustering analyses in existing SITVIT databases

Data obtained from new sequencing technologies are evolving rapidly, leading to the development of specific bioinformatic tools, pipelines and softwares. Several algorithms and tools are today available allowing a better identification and description of Mycobacterium tuberculosis complex (MTBC) isolates worldwide. Our approach consists in applying existing methods to analyze DNA sequencing data (from FASTA or FASTQ files), and tentatively extract meaningful information that would facilitate identification as well as a better understanding and management of MTBC isolates (taking into account whole genome sequencing and classical genotyping data). The aim of this study is to propose a pipeline analysis allowing to potentially simplify MTBC data analysis by providing different ways to interpret genomic or genotyping information based on existing tools. Furthermore, we propose a "reconciledTB" list making a link with results directly obtained from whole genome sequencing (WGS) data and results obtained from classical genotyping analysis (data inferred from SpoTyping and MIRUReader). Data visualization graphics and trees generated provide additional elements to better understand and confer associations among information overlap analyses. Additionally, comparison between data entered in an international genotyping database (SITVITEXTEND) and ensuing data obtained from the pipeline not only provide meaningful information, but further suggest that simpiTB could also be suitable for new data integration in specific TB genotyping databases.

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.

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.

KaruBioNet: a network and discussion group for a better collaboration and structuring of bioinformatics in Guadeloupe (French West Indies)

Summary

Sequencing and other biological data are now more frequently available and at a lower price. Mutual tools and strategies are needed to analyze the huge amount of heterogeneous data generated by several research teams and devices. Bioinformatics represents a growing field in the scientific community globally. This multidisciplinary field provides a great amount of tools and methods that can be used to conduct scientific studies in a more strategic way. Coordinated actions and collaborations are needed to find more innovative and accurate methods for a better understanding of real-life data. A wide variety of organizations are contributing to KaruBioNet in Guadeloupe (French West Indies), a Caribbean archipelago. The purpose of this group is to foster collaboration and mutual aid among people from different disciplines using a 'one health' approach, for a better comprehension and surveillance of humans, plants or animals' health and diseases. The KaruBioNet network particularly aims to help researchers in their studies related to 'omics' data, but also more general aspects concerning biological data analysis. This transdisciplinary network is a platform for discussion, sharing, training and support between scientists interested in bioinformatics and related fields. Starting from a little archipelago in the Caribbean, we envision to facilitate exchange between other Caribbean partners in the future, knowing that the Caribbean is a region with non-negligible biodiversity which should be preserved and protected. Joining forces with other Caribbean countries or territories would strengthen scientific collaborative impact in the region. Information related to this network can be found at: http://www.pasteur-guadeloupe.fr/karubionet.html. Furthermore, a dedicated 'Galaxy KaruBioNet' platform is available at: http://calamar.univ-ag.fr/c3i/galaxy_karubionet.html. Availability and implementation Information about KaruBioNet is availabe at: http://www.pasteur-guadeloupe.fr/karubionet.html.

Contact

dcouvin@pasteur-guadeloupe.fr.

Supplementary information

Supplementary data are available at Bioinformatics Advances online.

SITVITBovis—a publicly available database and mapping tool to get an improved overview of animal and human cases caused by Mycobacterium bovis

Limited data are available for bovine tuberculosis and the infections it can cause in humans and other mammals. We therefore constructed a publicly accessible SITVITBovis database that incorporates genotyping and epidemiological data on Mycobacterium bovis. It also includes limited data on Mycobacterium caprae (previously synonymous with the name M. bovis subsp. Caprae) that can infect both animals and humans. SITVITBovis incorporates data on 25,741 isolates corresponding to 60 countries of origin (75 countries of isolation). It reports a total of 1000 spoligotype patterns: 537 spoligotype international types (SITs, containing 25 278 clinical isolates) and 463 orphan patterns, allowing a wide overview of the geographic distribution of various phylogenetical sublineages (BOV_1, BOV_2, BOV_3 and BOV_4-CAPRAE). The SIT identifiers of the SITVITBovis were compared to the SB numbers of the Mbovis.org database to facilitate crosscheck among databases. Note that SITVITBovis also contains limited information on mycobacterial interspersed repetitive units-variable number of tandem repeats when available. Significant differences were observed when comparing age/gender of human isolates as well as various hosts. The database includes information on the regions where a strain was isolated as well as hosts involved, making it possible to see geographic trends. SITVITBovis is publicly accessible at: http://www.pasteur-guadeloupe.fr:8081/SITVIT_Bovis. Finally, a future second version is currently in progress to allow query of associated whole-genome sequencing data.

Database URLhttp://www.pasteur-guadeloupe.fr:8081/SITVIT_Bovis