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.

The recent introduction of Angiostrongylus cantonensis and its intermediate host Achatina fulica into Guadeloupe detected by phylogenetic analyses

BACKGROUND

Angiostrongylus cantonensis (rat lungworm) is the main pathogen responsible for eosinophilic meningitis in humans. One of its intermediate snail hosts, Achatina fulica, was already present in many countries around the world before it appeared in the West Indies in the late 1980s. In the French territories in the Caribbean and northern South America, the first cases of human neuroangiostrongyliasis were reported in Martinique, Guadeloupe and French Guiana in 2002, 2013 and 2017, respectively. In order to better characterize angiostrongyliasis in Guadeloupe, particularly its geographical origin and route of introduction, we undertook molecular characterization of adult worms of Angiostrongylus cantonensis and its intermediate host Achatina fulica.

METHODS

Genomic DNA of adult Angiostrongylus cantonensis and Achatina fulica was extracted and amplified by polymerase chain reaction (PCR) targeting the mitochondrial genes cytochrome B and C for A. cantonensis and 16S ribosomal RNA for A. fulica. The PCR products were sequenced and studied by phylogenetic analysis.

RESULTS

Cytochrome B and cytochrome C molecular markers indicate a monophyletic lineage of A. cantonensis adult worms in Guadeloupe. Two sequences of A. fulica were identified.

CONCLUSIONS

These results confirm the recent introduction of both Angiostrongylus cantonensis and Achatina fulica into Guadeloupe. Achatina fulica in Guadeloupe shares a common origin with those in Barbados and New Caledonia, while Angiostrongylus cantonensis in Guadeloupe shares a common origin with those in Brazil, Hawaii and Japan.

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

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

Draft genome sequence of Enterobacter chengduensis ECC445, isolated from fresh water in the West Indies

OBJECTIVES

The Enterobacter cloacae complex is considered an important opportunistic pathogen. It comprises many members that remain difficult to delineate by phenotypic approaches. Despite its importance in human infection, there is a lack of information on associated members in other compartments. Here we report the first de novo assembled and annotated whole-genome sequence of a E. chengduensis strain isolated from the environment.

DATA DESCRIPTION

ECC445 specimen was isolated in 2018 from a drinking water catchment point in Guadeloupe. It was clearly related to E. chengduensis species according to hsp60 typing and genomic comparison. Its whole-genome sequence is 5,211,280-bp long divided into 68 contigs, and presents a G + C content of 55.78%. This genome and associated datasets provided here will serve as a useful resource for further analyses of this rarely reported Enterobacter species.

Limited Transmission of Klebsiella pneumoniae among Humans, Animals, and the Environment in a Caribbean Island, Guadeloupe (French West Indies)

Guadeloupe (French West Indies), a Caribbean island, is an ideal place to study the reservoirs of the Klebsiella pneumoniae species complex (KpSC) and identify the routes of transmission between human and nonhuman sources due to its insularity, small population size, and small area. Here, we report an analysis of 590 biological samples, 546 KpSC isolates, and 331 genome sequences collected between January 2018 and May 2019. The KpSC appears to be common whatever the source. Extended-spectrum-β-lactamase (ESBL)-producing isolates (21.4%) belonged to K. pneumoniae sensu stricto (phylogroup Kp1), and all but one were recovered from the hospital setting. The distribution of species and phylogroups across the different niches was clearly nonrandom, with a distinct separation of Kp1 and Klebsiella variicola (Kp3). The most frequent sequence types (STs) (≥5 isolates) were previously recognized as high-risk multidrug-resistant (MDR) clones, namely, ST17, ST307, ST11, ST147, ST152, and ST45. Only 8 out of the 63 STs (12.7%) associated with human isolates were also found in nonhuman sources. A total of 22 KpSC isolates were defined as hypervirulent: 15 associated with human infections (9.8% of all human isolates), 4 (8.9%) associated with dogs, and 3 (15%) associated with pigs. Most of the human isolates (33.3%) belonged to the globally successful sublineage CG23-I. ST86 was the only clone shared by a human and a nonhuman (dog) source. Our work shows the limited transmission of KpSC isolates between human and nonhuman sources and points to the hospital setting as a cornerstone of the spread of MDR clones and antibiotic resistance genes. IMPORTANCE In this study, we characterized the presence and genomic features of isolates of the Klebsiella pneumoniae species complex (KpSC) from human and nonhuman sources in Guadeloupe (French West Indies) in order to identify the reservoirs and routes of transmission. This is the first study in an island environment, an ideal setting that limits the contribution of external imports. Our data showed the limited transmission of KpSC isolates between the different compartments. In contrast, we identified the hospital setting as the epicenter of antibiotic resistance due to the nosocomial spread of successful multidrug-resistant (MDR) K. pneumoniae clones and antibiotic resistance genes. Ecological barriers and/or limited exposure may restrict spread from the hospital setting to other reservoirs and vice versa. These results highlight the need for control strategies focused on health care centers, using genomic surveillance to limit the spread, particularly of high-risk clones, of this important group of MDR pathogens.

Emergence of a Novel Lineage and Wide Spread of a blaCTX-M-15/IncHI2/ST1 Plasmid among Nosocomial Enterobacter in Guadeloupe

Between April 2018 and August 2019, a total of 135 strains of Enterobacter cloacae complex (ECC) were randomly collected at the University Hospital Center of Guadeloupe to investigate the structure and diversity of the local bacterial population. These nosocomial isolates were initially identified genetically by the hsp60 typing method, which revealed the clinical relevance of E. xiangfangensis (n = 69). Overall, 57/94 of the third cephalosporin-resistant strains were characterized as extended-spectrum-β-lactamase (ESBL) producers, and their whole-genome was sequenced using Illumina technology to determine the clonal relatedness and diffusion of resistance genes. We found limited genetic diversity among sequence types (STs). ST114 (n = 13), ST1503 (n = 9), ST53 (n = 5) and ST113 (n = 4), which belong to three different Enterobacter species, were the most prevalent among the 57 ESBL producers. The blaCTXM-15 gene was the most prevalent ESBL determinant (56/57) and was in most cases associated with IncHI2/ST1 plasmid replicon carriage (36/57). To fully characterize this predominant blaCTXM-15/IncHI2/ST1 plasmid, four isolates from different lineages were also sequenced using Oxford Nanopore sequencing technology to generate long-reads. Hybrid sequence analyses confirmed the circulation of a well-conserved plasmid among ECC members. In addition, the novel ST1503 and its associated species (ECC taxon 4) were analyzed, in view of its high prevalence in nosocomial infections. These genetic observations confirmed the overall incidence of nosocomial ESBL Enterobacteriaceae infections acquired in this hospital during the study period, which was clearly higher in Guadeloupe (1.59/1000 hospitalization days) than in mainland France (0.52/1,000 hospitalization days). This project revealed issues and future challenges for the management and surveillance of nosocomial and multidrug-resistant Enterobacter in the Caribbean.

getSequenceInfo: a suite of tools allowing to get genome sequence information from public repositories

Background

Biological sequences are increasing rapidly and exponentially worldwide. Nucleotide sequence databases play an important role in providing meaningful genomic information on a variety of biological organisms.

Results

The getSequenceInfo software tool allows to access sequence information from various public repositories (GenBank, RefSeq, and the European Nucleotide Archive), and is compatible with different operating systems (Linux, MacOS, and Microsoft Windows) in a programmatic way (command line) or as a graphical user interface. getSequenceInfo or gSeqI v1.0 should help users to get some information on queried sequences that could be useful for specific studies (e.g. the country of origin/isolation or the release date of queried sequences). Queries can be made to retrieve sequence data based on a given kingdom and species, or from a given date. This program allows the separation between chromosomes and plasmids (or other genetic elements/components) by arranging each component in a given folder. Some basic statistics are also performed by the program (such as the calculation of GC content for queried assemblies). An empirically designed nucleotide ratio is calculated using nucleotide information in order to tentatively provide a “NucleScore” for studied genome assemblies. Besides the main gSeqI tool, other additional tools have been developed to perform various tasks related to sequence analysis.

Conclusion

The aim of this study is to democratize the use of public repositories in programmatic ways, and to facilitate sequence data analysis in a pedagogical perspective. Output results are available in FASTA, FASTQ, Excel/TSV or HTML formats. The program is freely available at: https://github.com/karubiotools/getSequenceInfo. getSequenceInfo and supplementary tools are partly available through the recently released Galaxy KaruBioNet platform (http://calamar.univ-ag.fr/c3i/galaxy_karubionet.html).

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-022-04809-5.

High Prevalence of bla CTXM-1/IncI1-Iγ/ST3 Plasmids in Extended-Spectrum β-Lactamase-Producing Escherichia coli Isolates Collected From Domestic Animals in Guadeloupe (French West Indies)

Extended-spectrum β-lactamase-producing Enterobacteriaceae (ESBL-E) have been classified in the group of resistant bacteria of highest priority. We determined the prevalence of ESBL-E collected in feces from household and shelter pets in Guadeloupe (French West Indies). A single rectal swab was taken from 125 dogs and 60 cats between June and September 2019. The prevalence of fecal carriage of ESBL-E was 7.6% (14/185, 95% CI: 4.2-12.4), within the range observed worldwide. The only risk factor associated with a higher prevalence of ESBL-E rectal carriage was a stay in a shelter, suggesting that refuges could be hotspots for their acquisition. All but one (Klebsiella pneumoniae from a cat) were Escherichia coli. We noted the presence of a bla _CTX-M-1/IncI1-Iγ/sequence type (ST3) plasmid in 11 ESBL-producing E. coli isolates belonging to ST328 (n = 6), ST155 (n = 4) and ST953 (n = 1). A bla _CTX-M-15 gene was identified in the three remaining ESBL-E isolates. The bla _CTX-M-1 and most of the antimicrobial resistance genes were present in a well-conserved large conjugative IncI1-Iγ/ST3 plasmid characterized by two accessory regions containing antibiotic resistance genes. The plasmid has been detected worldwide in E. coli isolates from humans and several animal species, such as food-producing animals, wild birds and pets, and from the environment. This study shows the potential role of pets as a reservoir of antimicrobial-resistant bacteria or genes for humans and underlines the importance of basic hygiene measures by owners of companion animals.

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

Wide Distribution and Specific Resistance Pattern to Third-Generation Cephalosporins of Enterobacter cloacae Complex Members in Humans and in the Environment in Guadeloupe (French West Indies)

Species belonging to Enterobacter cloacae complex have been isolated in numerous environments and samples of various origins. They are also involved in opportunistic infections in plants, animals, and humans. Previous prospection in Guadeloupe (French West Indies) indicated a high frequency of E. cloacae complex strains resistant to third-generation cephalosporins (3GCs) in a local lizard population (Anolis marmoratus), but knowledge of the distribution and resistance of these strains in humans and the environment is limited. The aim of this study was to compare the distribution and antibiotic susceptibility pattern of E. cloacae complex members from different sources in a “one health” approach and to find possible explanations for the high level of resistance in non-human samples. E. cloacae complex strains were collected between January 2017 and the end of 2018 from anoles, farm animals, local fresh produce, water, and clinical human samples. Isolates were characterized by the heat-shock protein 60 gene-fragment typing method, and whole-genome sequencing was conducted on the most frequent clusters (i.e., C-VI and C-VIII). The prevalence of resistance to 3GCs was relatively high (56/346, 16.2%) in non-human samples. The associated resistance mechanism was related to an AmpC overproduction; however, in human samples, most of the resistant strains (40/62) produced an extended-spectrum beta-lactamase. No relation was found between resistance in isolates from wild anoles (35/168) and human activities. Specific core-genome phylogenetic analysis highlighted an important diversity in this bacterial population and no wide circulation among the different compartments. In our setting, the mutations responsible for resistance to 3GCs, especially in ampD, were diverse and not compartment specific. In conclusion, high levels of resistance in non-human E. cloacae complex isolates are probably due to environmental factors that favor the selection of these resistant strains, and this will be explored further.

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.

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

Bioinformatic tools are currently being developed to better understand the Mycobacterium tuberculosis complex (MTBC). Several approaches already exist for the identification of MTBC lineages using classical genotyping methods such as mycobacterial interspersed repetitive units-variable number of tandem DNA repeats and spoligotyping-based families. In the recently released SITVIT2 proprietary database of the Institut Pasteur de la Guadeloupe, a large number of spoligotype families were assigned by either manual curation/expertise or using an in-house algorithm. In this study, we present two complementary data-driven approaches allowing fast and precise family prediction from spoligotyping patterns. The first one is based on data transformation and the use of decision tree classifiers. In contrast, the second one searches for a set of simple rules using binary masks through a specifically designed evolutionary algorithm. The comparison with the three main approaches in the field highlighted the good performances of our contributions and the significant runtime gain. Finally, we propose the 'SpolLineages' software tool (https://github.com/dcouvin/SpolLineages), which implements these approaches for MTBC spoligotype families' identification.

Population structure of multidrug-resistant Mycobacterium tuberculosis clinical isolates in Colombia

Emergence of multidrug-resistant (MDR) Mycobacterium tuberculosis complex (MTBC) isolates is a major public health problem that threatens progress made in tuberculosis (TB) care and control worldwide. In Colombia, the prevalence of MDR tuberculosis (MDR-TB) has increased slowly but steadily since 2001. However, the population structure of the MDR-TB strains circulating in Colombia is sparsely known. In this work, 203 MDR isolates isolated in 2012-2013 were collected, and characterized by spoligotyping, followed by 24-loci MIRU-VNTR (data available for 190 isolates). The most prevalent genotypes corresponded to SIT42/LAM9 (12.81%), SIT62/H1 (10.34%), and SIT190/Beijing (10.34%). A fine analysis showed that although the MDR strains came from 29 of the 33 departments of Colombia, the distribution of these main lineages was not at random and depended on the city of isolation (p-value <0.000001). Both LAM and Beijing lineage strains were significantly associated with MDR-TB (p-value <0.0001): LAM lineage was associated with 2 patterns of MDR, namely combined resistance to INH + Rifampin (HR), and to SHRE (Streptomycin + INH + Rifampin + Ethambutol), while the Beijing lineage strains were essentially associated with MDR (SHRE). Interestingly, distribution of genotypic lineages in function of drug resistance information (e.g. pansusceptible vs. MDR) was different in our setting as compared to other countries in Latin America. However, MIRU-VNTR patterns were unique for all strains, an observation that did not support active transmission of circulating MDR clones.

Antimicrobial Resistance in Wildlife in Guadeloupe (French West Indies): Distribution of a Single bla CTX-M-1/IncI1/ST3 Plasmid Among Humans and Wild Animals

Limited data are available on the contribution of wildlife to the spread of antibacterial resistance. We determined the prevalence of resistance to antibiotics in Escherichia coli isolates collected from wild animals in 2013 and 2014 and the genetic basis for resistance to third-generation cephalosporin in Guadeloupe. We recovered 52 antibiotic-resistant (AR) E. coli strains from 48 of the 884 (5.4%) wild animals tested (46 iguanas, 181 birds, 289 anoles, and 368 rodents at 163 sampling sites). Rodents had higher rates of carriage (n = 38, 10.3%) than reptiles and birds (2.4% and 1.1%, respectively, p < 0.001). A significant association (p < 0.001) was found between the degree of anthropization and the frequency of AR E. coli carriage for all species. The carriage rate of ciprofloxacin- and cefotaxime-resistant isolates was 0.7% (6/884) and 1.5% (13/884), respectively. Most (65.4%) AR E. coli were multi-drug resistant, and the prevalence of extended-spectrum beta-lactamase (ESBL)-producing E. coli was low (n = 7, 0.8%) in all species. Eight ESBL-producing E. coli were recovered, two genetically unrelated isolates being found in one bird. These isolates and 20 human invasive ESBL E. coli isolates collected in Guadeloupe during the same period were investigated by whole genome sequencing. bla_CTX–M–1 was the only ESBL gene shared by three animal classes (humans, n = 2; birds, n = 2; rodents, n = 2). The bla_CTX–M–1 gene and most of the antimicrobial resistance genes were present in a large conjugative IncI1 plasmid that was highly similar (>99% nucleotide identity) to ESBL-carrying plasmids found in several countries in Europe and in Australia. Although the prevalence of ESBL-producing E. coli isolates was very low in wild animals, it is of concern that the well-conserved IncI1 plasmid-carrying bla_CTX–M–1 is widespread and occurs in various E. coli strains from animals and humans.

Local adaptive evolution of two distinct clades of Beijing and T families of Mycobacterium tuberculosis in Chongqing: a Bayesian population structure and phylogenetic study

Background

Beijing sub-pedigree 2 (BSP2) and T sub-lineage 6 (TSL6) are two clades belonging to Beijing and T family of Mycobacterium tuberculosis (MTB), respectively, defined by Bayesian population structure analysis based on 24-loci mycobacterial interspersed repetitive unit-variable number of tandem repeats (MIRU-VNTR). Globally, over 99% of BSP2 and 89% of TSL6 isolates were distributed in Chongqing, suggesting their possible local adaptive evolution. The objective of this paper is to explore whether BSP2 and TSL6 originated by their local adaptive evolution from the specific isolates of Beijing and T families in Chongqing.

Methods

The genotyping data of 16 090 MTB isolates were collected from laboratory collection, published literatures and SITVIT database before subjected to Bayesian population structure analysis based on 24-loci MIRU-VNTR. Spacer Oligonucleotide Forest (Spoligoforest) and 24-loci MIRU-VNTR-based minimum spanning tree (MST) were used to explore their phylogenetic pathways, with Bayesian demographic analysis for exploring the recent demographic change of TSL6.

Results

Phylogenetic analysis suggested that BSP2 and TSL6 in Chongqing may evolve from BSP4 and TSL5, respectively, which were locally predominant in Tibet and Jiangsu, respectively. Spoligoforest showed that Beijing and T families were genetically distant, while the convergence of the MIRU-VNTR pattern of BSP2 and TSL6 was revealed by WebLogo. The demographic analysis concluded that the recent demographic change of TSL6 might take 111.25 years.

Conclusions

BSP2 and TSL6 clades might originate from BSP4 and TSL5, respectively, by their local adaptive evolution in Chongqing. Our study suggests MIRU-VNTR be combined with other robust markers for a more comprehensive genotyping approach, especially for families of clades with the same MIRU-VNTR pattern.

CRISPRCasdb a successor of CRISPRdb containing CRISPR arrays and cas genes from complete genome sequences, and tools to download and query lists of repeats and spacers

In Archaea and Bacteria, the arrays called CRISPRs for 'clustered regularly interspaced short palindromic repeats' and the CRISPR associated genes or cas provide adaptive immunity against viruses, plasmids and transposable elements. Short sequences called spacers, corresponding to fragments of invading DNA, are stored in-between repeated sequences. The CRISPR-Cas systems target sequences homologous to spacers leading to their degradation. To facilitate investigations of CRISPRs, we developed 12 years ago a website holding the CRISPRdb. We now propose CRISPRCasdb, a completely new version giving access to both CRISPRs and cas genes. We used CRISPRCasFinder, a program that identifies CRISPR arrays and cas genes and determine the system's type and subtype, to process public whole genome assemblies. Strains are displayed either in an alphabetic list or in taxonomic order. The database is part of the CRISPR-Cas++ website which also offers the possibility to analyse submitted sequences and to download programs. A BLAST search against lists of repeats and spacers extracted from the database is proposed. To date, 16 990 complete prokaryote genomes (16 650 bacteria from 2973 species and 340 archaea from 300 species) are included. CRISPR-Cas systems were found in 36% of Bacteria and 75% of Archaea strains. CRISPRCasdb is freely accessible at https://crisprcas.i2bc.paris-saclay.fr/.

CRISPRCasdb a successor of CRISPRdb containing CRISPR arrays and cas genes from complete genome sequences, and tools to download and query lists of repeats and spacers

In Archaea and Bacteria, the arrays called CRISPRs for 'clustered regularly interspaced short palindromic repeats' and the CRISPR associated genes or cas provide adaptive immunity against viruses, plasmids and transposable elements. Short sequences called spacers, corresponding to fragments of invading DNA, are stored in-between repeated sequences. The CRISPR-Cas systems target sequences homologous to spacers leading to their degradation. To facilitate investigations of CRISPRs, we developed 12 years ago a website holding the CRISPRdb. We now propose CRISPRCasdb, a completely new version giving access to both CRISPRs and cas genes. We used CRISPRCasFinder, a program that identifies CRISPR arrays and cas genes and determine the system's type and subtype, to process public whole genome assemblies. Strains are displayed either in an alphabetic list or in taxonomic order. The database is part of the CRISPR-Cas++ website which also offers the possibility to analyse submitted sequences and to download programs. A BLAST search against lists of repeats and spacers extracted from the database is proposed. To date, 16 990 complete prokaryote genomes (16 650 bacteria from 2973 species and 340 archaea from 300 species) are included. CRISPR-Cas systems were found in 36% of Bacteria and 75% of Archaea strains. CRISPRCasdb is freely accessible at https://crisprcas.i2bc.paris-saclay.fr/.

Two tales: Worldwide distribution of Central Asian (CAS) versus ancestral East-African Indian (EAI) lineages of Mycobacterium tuberculosis underlines a remarkable cleavage for phylogeographical, epidemiological and demographical characteristics

The East African Indian (EAI) and Central Asian (CAS) lineages of Mycobacterium tuberculosis complex (MTBC) mainly infect tuberculosis (TB) patients in the eastern hemisphere which contains many of the 22 high TB burden countries including China and India. We investigated if phylogeographical, epidemiological and demographical characteristics for these 2 lineages differed in SITVIT2 database. Genotyping results and associated data (age, sex, HIV serology, drug resistance) on EAI and CAS lineages (n = 10,974 strains) were extracted. Phylogenetic and Bayesian, and other statistical analyses were used to compare isolates. The male/female sex ratio was 907/433 (2.09) for the EAI group vs. 881/544 (1.62) for CAS (p-value<0.002). The proportion of younger patients aged 0-20 yrs. with CAS lineage was significantly higher than for EAI lineage (18.07% vs. 10.85%, p-value<0.0001). The proportion of multidrug resistant and extensively drug resistant TB among CAS group (30.63% and 1.03%, respectively) was significantly higher than in the EAI group (12.14% and 0.29%, respectively; p-value<0.0001). Lastly, the proportion of HIV+ patients was 20.34% among the EAI group vs. 3.46% in the CAS group (p-value<0.0001). This remarkable split observed between various parameters for these 2 lineages was further corroborated by their geographic distribution profile (EAI being predominantly found in Eastern-Coast of Africa, South-India and Southeast Asia, while CAS was predominantly found in Afghanistan, Pakistan, North India, Nepal, Middle-east, Libya, Sudan, Ethiopia, Kenya and Tanzania). Some geo-specificities were highlighted. This study demonstrated a remarkable cleavage for aforementioned characteristics of EAI and CAS lineages, showing a North-South divide along the tropic of cancer in Eastern hemisphere-mainly in Asia, and partly prolonged along the horn of Africa. Such studies would be helpful to better comprehend prevailing TB epidemic in context of its historical spread and evolutionary features, and provide clues to better treatment and patient-care in countries and regions concerned by these lineages.

Macro-geographical specificities of the prevailing tuberculosis epidemic as seen through SITVIT2, an updated version of the Mycobacterium tuberculosis genotyping database

In order to provide a global overview of genotypic, epidemiologic, demographic, phylogeographical, and drug resistance characteristics related to the prevailing tuberculosis (TB) epidemic, we hereby report an update of the 6th version of the international genotyping database SITVIT2. We also make all the available information accessible through a dedicated website (available at http://www.pasteur-guadeloupe.fr:8081/SITVIT2). Thanks to the public release of SITVIT2 which is currently the largest international multimarker genotyping database with a compilation of 111,635 clinical isolates from 169 countries of patient origin (131 countries of isolation, representing 1032 cities), our major aim is to highlight macro- and micro-geographical cleavages and phylogeographical specificities of circulating Mycobacterium tuberculosis complex (MTBC) clones worldwide. For this purpose, we retained strains typed by the most commonly used PCR-based methodology for TB genotyping, i.e., spoligotyping based on the polymorphism of the direct repeat (DR) locus, 5-loci Exact Tandem Repeats (ETRs), and MIRU-VNTR minisatellites used in 12-, 15-, or 24-loci formats. We describe the SITVIT2 database and integrated online applications that permit to interrogate the database using easy drop-down menus to draw maps, graphics and tables versus a long list of parameters and variables available for individual clinical isolates (year and place of isolation, origin, sex, and age of patient, drug-resistance, etc.). Available tools further allow to generate phylogenetical snapshot of circulating strains as Lineage-specific WebLogos, as well as minimum spanning trees of their genotypes in conjunction with their geographical distribution, drug-resistance, demographic, and epidemiologic characteristics instantaneously; whereas online statistical analyses let a user to pinpoint phylogeographical specificities of circulating MTBC lineages and conclude on actual demographic trends. Available associated information on gender (n = 18,944), age (n = 16,968), drug resistance (n = 19,606), and HIV serology (n = 2673), allowed to draw some important conclusions on TB geo-epidemiology; e.g. a positive correlation exists between certain Mycobacterium tuberculosis lineages (such as CAS and Beijing) and drug resistance (p-value<.001), while other lineages (such as LAM, X, and BOV) are more frequently associated with HIV-positive serology (p-value<.001). Besides, availability of information on the year of isolation of strains (range 1759-2012), also allowed to make tentative correlations between drug resistance information and lineages - portraying probable evolution trends over time and space. To conclude, the present approach of geographical mapping of predominant clinical isolates of tubercle bacilli causing the bulk of the disease both at country and regional level in conjunction with epidemiologic and demographic characteristics allows to shed new light on TB geo-epidemiology in relation with the continued waves of peopling and human migration.

Mycobacterial biomaterials and resources for researchers

There are many resources available to mycobacterial researchers, including culture collections around the world that distribute biomaterials to the general scientific community, genomic and clinical databases, and powerful bioinformatics tools. However, many of these resources may be unknown to the research community. This review article aims to summarize and publicize many of these resources, thus strengthening the quality and reproducibility of mycobacterial research by providing the scientific community access to authenticated and quality-controlled biomaterials and a wealth of information, analytical tools and research opportunities.

Geospatial distribution of Mycobacterium tuberculosis genotypes in Africa

Objective

To investigate the distribution of Mycobacterium tuberculosis genotypes across Africa.

Methods

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

Results

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

Conclusions

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

CRISPRCasFinder, an update of CRISRFinder, includes a portable version, enhanced performance and integrates search for Cas proteins

CRISPR (clustered regularly interspaced short palindromic repeats) arrays and their associated (Cas) proteins confer bacteria and archaea adaptive immunity against exogenous mobile genetic elements, such as phages or plasmids. CRISPRCasFinder allows the identification of both CRISPR arrays and Cas proteins. The program includes: (i) an improved CRISPR array detection tool facilitating expert validation based on a rating system, (ii) prediction of CRISPR orientation and (iii) a Cas protein detection and typing tool updated to match the latest classification scheme of these systems. CRISPRCasFinder can either be used online or as a standalone tool compatible with Linux operating system. All third-party software packages employed by the program are freely available. CRISPRCasFinder is available at https://crisprcas.i2bc.paris-saclay.fr.

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

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

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

•

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

•

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

•

Common Mycobacterium tuberculosis population structure denote historical strain flow.

•

Cross-border clusters suggest recent intercontinental tuberculosis transmission.

•

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

Analysis of Mycobacterium tuberculosis Genotypic Lineage Distribution in Chile and Neighboring Countries

Tuberculosis (TB), caused by the pathogen Mycobacterium tuberculosis (MTB), remains a disease of high importance to global public health. Studies into the population structure of MTB have become vital to monitoring possible outbreaks and also to develop strategies regarding disease control. Although Chile has a low incidence of MTB, the current rates of migration have the potential to change this scenario. We collected and analyzed a total of 458 M. tuberculosis isolates (1 isolate per patient) originating from all 15 regions of Chile. The isolates were genotyped using the spoligotyping method and the data obtained were analyzed and compared with the SITVIT2 database. A total of 169 different patterns were identified, of which, 119 patterns (408 strains) corresponded to Spoligotype International Types (SITs) and 50 patterns corresponded to orphan strains. The most abundantly represented SITs/lineages were: SIT53/T1 (11.57%), SIT33/LAM3 (9.6%), SIT42/LAM9 (9.39%), SIT50/H3 (5.9%), SIT37/T3 (5%); analysis of the spoligotyping minimum spanning tree as well as spoligoforest were suggestive of a recent expansion of SIT42, SIT50 and SIT37; all of which potentially evolved from SIT53. The most abundantly represented lineages were LAM (40.6%), T (34.1%) and Haarlem (13.5%). LAM was more prevalent in the Santiago (43.6%) and Concepción (44.1%) isolates, rather than the Iquique (29.4%) strains. The proportion of X lineage was appreciably higher in Iquique and Concepción (11.7% in both) as compared to Santiago (1.6%). Global analysis of MTB lineage distribution in Chile versus neighboring countries showed that evolutionary recent lineages (LAM, T and Haarlem) accounted together for 88.2% of isolates in Chile, a pattern which mirrored MTB lineage distribution in neighboring countries (n = 7378 isolates recorded in SITVIT2 database for Peru, Brazil, Paraguay, and Argentina; and published studies), highlighting epidemiological advantage of Euro-American lineages in this region. Finally, we also observed exclusive emergence of patterns SIT4014/X1 and SIT4015 (unknown lineage signature) that have hitherto been found exclusively in Chile, indicating that conditions specific to Chile, along with the unique genetic makeup of the Chilean population, might have allowed for a possible co-evolution leading to the success of these emerging genotypes.

Mapping of Mycobacterium tuberculosis Complex Genetic Diversity Profiles in Tanzania and Other African Countries

The aim of this study was to assess and characterize Mycobacterium tuberculosis complex (MTBC) genotypic diversity in Tanzania, as well as in neighbouring East and other several African countries. We used spoligotyping to identify a total of 293 M. tuberculosis clinical isolates (one isolate per patient) collected in the Bunda, Dar es Salaam, Ngorongoro and Serengeti areas in Tanzania. The results were compared with results in the SITVIT2 international database of the Pasteur Institute of Guadeloupe. Genotyping and phylogeographical analyses highlighted the predominance of the CAS, T, EAI, and LAM MTBC lineages in Tanzania. The three most frequent Spoligotype International Types (SITs) were: SIT21/CAS1-Kili (n = 76; 25.94%), SIT59/LAM11-ZWE (n = 22; 7.51%), and SIT126/EAI5 tentatively reclassified as EAI3-TZA (n = 18; 6.14%). Furthermore, three SITs were newly created in this study (SIT4056/EAI5 n = 2, SIT4057/T1 n = 1, and SIT4058/EAI5 n = 1). We noted that the East-African-Indian (EAI) lineage was more predominant in Bunda, the Manu lineage was more common among strains isolated in Ngorongoro, and the Central-Asian (CAS) lineage was more predominant in Dar es Salaam (p-value<0.0001). No statistically significant differences were noted when comparing HIV status of patients vs. major lineages (p-value = 0.103). However, when grouping lineages as Principal Genetic Groups (PGG), we noticed that PGG2/3 group (Haarlem, LAM, S, T, and X) was more associated with HIV-positive patients as compared to PGG1 group (Beijing, CAS, EAI, and Manu) (p-value = 0.03). This study provided mapping of MTBC genetic diversity in Tanzania (containing information on isolates from different cities) and neighbouring East African and other several African countries highlighting differences as regards to MTBC genotypic distribution between Tanzania and other African countries. This work also allowed underlining of spoligotyping patterns tentatively grouped within the newly designated EAI3-TZA lineage (remarkable by absence of spacers 2 and 3, and represented by SIT126) which seems to be specific to Tanzania. However, further genotyping information would be needed to confirm this specificity.

First insight into the molecular epidemiology of Mycobacterium tuberculosis in Santa Catarina, southern Brazil

Molecular epidemiology of Mycobacterium tuberculosis is useful for understanding disease transmission dynamics, and to establish strategic measures for TB control and prevention. The aim of this study was to analyze clinical, epidemiological and molecular characteristics of MTBC clinical isolates from Santa Catarina state, southern Brazil. During one-year period, 406 clinical isolates of MTBC were collected from Central Laboratory of Public Health and typed by spoligotyping. Demographic and clinical data were collected from the Brazilian National Mandatory Disease Reporting System. The majority of cases occurred in highest population densities regions and about 50% had some condition associated with TB. Among all isolates, 5.7% were MDR, which showed association with drug addiction. LAM was the most predominant lineage with 47.5%, followed by the T superfamily with 25.9% and Haarlem with 12.3%. The MST showed two major groups: the first was formed mainly by the LAM lineage and the second was mainly formed by the T and Haarlem lineages. Others lineages were distributed in peripheral positions. This study provides the first insight into the population structure of M. tuberculosis in SC State. Spoligotyping and other genotyping analyses are important to establish strategic measures for TB control and prevention.

Mycobacterium tuberculosis causing tuberculous lymphadenitis in Maputo, Mozambique

Background

The zoonosis bovine tuberculosis (TB) is known to be responsible for a considerable proportion of extrapulmonary TB. In Mozambique, bovine TB is a recognised problem in cattle, but little has been done to evaluate how Mycobacterium bovis has contributed to human TB. We here explore the public health risk for bovine TB in Maputo, by characterizing the isolates from tuberculous lymphadenitis (TBLN) cases, a common manifestation of bovine TB in humans, in the Pathology Service of Maputo Central Hospital, in Mozambique, during one year.

Results

Among 110 patients suspected of having TBLN, 49 had a positive culture result. Of those, 48 (98 %) were positive for Mycobacterium tuberculosis complex and one for nontuberculous mycobacteria. Of the 45 isolates analysed by spoligotyping and Mycobacterial Interspersed Repetitive Unit – Variable Number Tandem Repeat (MIRU-VNTR), all were M. tuberculosis. No M. bovis was found.

Cervical TBLN, corresponding to 39 (86.7 %) cases, was the main cause of TBLN and 66.7 % of those where from HIV positive patients.

We found that TBLN in Maputo was caused by a variety of M. tuberculosis strains. The most prevalent lineage was the EAI (n = 19; 43.2 %). Particular common spoligotypes were SIT 48 (EAI1_SOM sublineage), SIT 42 (LAM 9), SIT 1 (Beijing) and SIT53 (T1), similar to findings among pulmonary cases.

Conclusions

M. tuberculosis was the main etiological agent of TBLN in Maputo. M. tuberculosis genotypes were similar to the ones causing pulmonary TB, suggesting that in Maputo, cases of TBLN arise from the same source as pulmonary TB, rather than from an external zoonotic source. Further research is needed on other forms of extrapulmonary TB and in rural areas where there is high prevalence of bovine TB in cattle, to evaluate the risk of transmission of M. bovis from cattle to humans.

Electronic supplementary material

The online version of this article (doi:10.1186/s12866-015-0603-5) contains supplementary material, which is available to authorized users.

Predominance of Uganda genotype of Mycobacterium tuberculosis isolated from Ugandan patients with tuberculous lymphadenitis

Background

In Uganda, the emerging Uganda genotype of Mycobacterium tuberculosis is the most common cause of pulmonary tuberculosis (PTB), and accounts for up to 70 % of isolates. Extrapulmonary TB (EPTB) is less studied in Uganda.

Methods

Molecular characterization using deletion analysis and spoligotyping was performed on 121 M. tuberculosis isolates from lymph node fine needle biopsy aspirates of Ugandan patients with tuberculous lymphadenitis. The evolutionary relationships and worldwide distribution of the spoligotypes were analyzed.

Results

Mycobacterium tuberculosis was the only cause of EPTB in this study. The T2 sublineage was the most predominant lineage and the Uganda genotype was the dominant genotype. There were 54 spoligotype patterns among the 121 study isolates. The dominant spoligotypes were shared international types (SIT) SIT420, SIT53, SIT 135, SIT 128 and SIT590 in descending order. All but SIT420 were previously reported in pulmonary TB in this setting. The phylogenetic analysis showed a long descendant branch of spoligotypes belonging to the T2-Uganda sublineage containing specifically SITs 135, 128 and 420.

Conclusion

In most cases, the spoligotypes were similar to those causing PTB, but the Uganda genotype was found to be less common in EPTB than previously reported for PTB in Uganda. The phylogenetic analysis and the study of the worldwide distribution of clustered spoligotypes indicate an ongoing evolution of the Uganda genotype, with the country of Uganda at the center of this evolution.

Electronic supplementary material

The online version of this article (doi:10.1186/s13104-015-1362-y) contains supplementary material, which is available to authorized users.

Genotypic characterization and historical perspective of Mycobacterium tuberculosis among older and younger Finns, 2008-2011

The Mycobacterium tuberculosis genotypes obtained from elderly Finns were assessed and compared with those obtained from younger Finns to comprehend the epidemiology of tuberculosis (TB) in Finland. From 2008 to 2011, a total of 1021 M. tuberculosis isolates were characterized by spoligotyping and 15-locus mycobacterial interspersed repetitive units-variable number tandem repeat typing. In total, 733 Finnish-born cases were included in the study, of which 466 (64%) were born before 1945 (older Finns). Of these, 63 (14%) shared an M. tuberculosis genotype with foreign-born or younger Finnish cases (born after 1945), and 59 (13%) shared a genotype with older Finnish cases. Eighty-five per cent had a unique genotypic profile while 70% belonged to T or Haarlem families, suggesting that ongoing transmission is infrequent among young and elderly Finns. Simultaneous reactivation of TB among older Finns was the most likely cause for clustering. As most isolates belonged to Haarlem or T, Finland was most likely affected by a similar TB epidemic at the beginning of the twentieth century as that seen in Sweden and Norway. Younger Finns were significantly more likely to be clustered (56% versus 27%, p<0.001), have pulmonary TB (87% versus 71%, p<0.001) and to be sputum smear positive (57% versus 48%, p<0.05) indicating that the risk of TB transmission from younger Finns is likely to be larger than from older Finns. The M. tuberculosis isolates from elderly Finns were associated with dominant lineages of the early twentieth century and differed from the heterogeneous lineages found among younger TB patients. Additionally, younger TB patients were more likely to transmit TB than elderly Finns.

Genetic diversity and drug susceptibility profile of Mycobacterium tuberculosis isolated from different regions of India

OBJECTIVES

Molecular genotyping profiles of Mycobacterium tuberculosis (MTB) provide a valuable insight into the evolution and transmission of the bacilli. Due to the lack of comprehensive national level data from India on this subject, we performed this study to determine the recent trends and distribution of various MTB lineages circulating in India.

METHODS

A total of 628 MTB isolates were obtained from North, West, South, Central and Eastern India. Spoligotyping and drug susceptibility testing was performed by using manufacturer's instructions.

RESULTS

Spoligotyping detected 102 distinct spoligo-patterns. A total of 536 (85.3%) isolates were distributed into 85 SITs which matched the pre-existing database, whereas 17 SITs were newly created for 34 (5.4%) isolates. Overall, CAS family genotype was predominant, comprising 222 (35.4%) isolates, followed by EAI in 152 (24.2%), Beijing in 108 (17.2%), Manu in 41 (6.5%), T in 30 (4.8%), H in 6 (0.9%), X in 3 (0.5%) and one (0.2%) each in Ural and AFRI. Drug susceptibility testing identified 134 (21.3%) isolates as multi drug resistant (MDR).

CONCLUSIONS

The CAS lineage had a pan India presence but EAI lineage was confined to southern parts of India. Beijing genotype of MTB was significantly associated (p-value <0.0001) with MDR.

Phylogenetic associations with drug-resistant Mycobacterium tuberculosis isolates in a paediatric population

SETTING

Department of Microbiology, National Institute of Tuberculosis and Respiratory Diseases, New Delhi, India.

OBJECTIVES

As paediatric tuberculosis (TB) is a surrogate marker for actively transmitted disease in a community, we investigated drug resistance patterns of 97 Mycobacterium tuberculosis complex strains isolated from children and explored their phylogenetic associations.

DESIGN

A total of 111 paediatric patients who attended the out-patient department during the study period 2009-2011 and whose sputum samples were sent to the Microbiology Department for liquid culture and drug susceptibility testing (DST) were included in this study. DST and spoligotyping were performed on cultures positive for M. tuberculosis complex.

RESULTS

DST against four first-line drugs showed that 31 of 97 (32%) strains were pan-susceptible, while 66/97 (68%) were resistant to at least one drug, including 55/97 (56.7%) that were resistant to at least isoniazid and rifampicin (i.e., multidrug-resistant). The majority of the isolates (n = 81/90, 90%) belonged to the principal genetic group 1 strains, the most predominant spoligotyping clusters being spoligotyping international type (SIT)1/Beijing (n = 28), SIT26/CAS1-Delhi (n = 27) and SIT53/T1 (n = 6).

CONCLUSION

The involvement of Beijing and CAS1-Delhi clades in paediatric TB patients suggests that these two lineages play a major role in ongoing active transmission.

The Association between Mycobacterium Tuberculosis Genotype and Drug Resistance in Peru

Background

The comparison of Mycobacterium tuberculosis bacterial genotypes with phenotypic, demographic, geospatial and clinical data improves our understanding of how strain lineage influences the development of drug-resistance and the spread of tuberculosis.

Methods

To investigate the association of Mycobacterium tuberculosis bacterial genotype with drug-resistance. Drug susceptibility testing together with genotyping using both 15-loci MIRU-typing and spoligotyping, was performed on 2,139 culture positive isolates, each from a different patient in Lima, Peru. Demographic, geospatial and socio-economic data were collected using questionnaires, global positioning equipment and the latest national census.

Results

The Latin American Mediterranean (LAM) clade (OR 2.4, p<0.001) was significantly associated with drug-resistance and alone accounted for more than half of all drug resistance in the region. Previously treated patients, prisoners and genetically clustered cases were also significantly associated with drug-resistance (OR's 2.5, 2.4 and 1.8, p<0.001, p<0.05, p<0.001 respectively).

Conclusions

Tuberculosis disease caused by the LAM clade was more likely to be drug resistant independent of important clinical, genetic and socio-economic confounding factors. Explanations for this include; the preferential co-evolution of LAM strains in a Latin American population, a LAM strain bacterial genetic background that favors drug-resistance or the "founder effect" from pre-existing LAM strains disproportionately exposed to drugs.

A first insight on the population structure of Mycobacterium tuberculosis complex as studied by spoligotyping and MIRU-VNTRs in Santiago, Chile

Tuberculosis (TB) remains a significant public health problem worldwide, but the ecology of the prevalent mycobacterial strains, and their transmission, can vary depending on country and region. Chile is a country with low incidence of TB, that has a geographically isolated location in relation to the rest of South American countries due to the Andes Mountains, but recent migration from neighboring countries has changed this situation. We aimed to assess the genotypic diversity of Mycobacterium tuberculosis complex (MTBC) strains in Santiago, Chile, and compare with reports from other Latin-American countries. We analyzed MTBC isolates from pulmonary tuberculosis cases collected between years 2008 and 2013 in Central Santiago, using two genotyping methods: spoligotyping and 12-loci mycobacterial interspersed repetitive unit-variable number of tandem repeats (MIRU-VNTRs). Data obtained were analyzed and compared to the SITVIT2 database. Mean age of the patients was 47.5 years and 61% were male; 11.6% were migrants. Of 103 strains (1 isolate/patient) included, there were 56 distinct spoligotype patterns. Of these, 16 strains (15.5%) corresponded to orphan strains in the SITVIT2 database, not previously reported. Latin American and Mediterranean (LAM) (34%) and T (33%) lineages were the most prevalent strains, followed by Haarlem lineage (16.5%). Beijing family was scarcely represented with only two cases (1.9%), one of them isolated from a Peruvian migrant. The most frequent clustered spoligotypes were SIT33/LAM3 (10.7%), SIT53/T1 (8.7%), SIT50/H3 (7.8%), and SIT37/T3 (6.8%). We conclude that LAM and T genotypes are the most prevalent genotypes of MTBC in Santiago, Chile, and together correspond to almost two thirds of analyzed strains, which is similar to strain distribution reported from other countries of Latin America. Nevertheless, the high proportion of SIT37/T3, which was rarely found in other Latin American countries, may underline a specific history or demographics of Chile related to probable human migrations and evolutions.

Study of Mycobacterium tuberculosis complex genotypic diversity in Malaysia reveals a predominance of ancestral East-African-Indian lineage with a Malaysia-specific signature

Background

Tuberculosis (TB) still constitutes a major public health problem in Malaysia. The identification and genotyping based characterization of Mycobacterium tuberculosis complex (MTBC) isolates causing the disease is important to determine the effectiveness of the control and surveillance programs.

Objectives

This study intended a first assessment of spoligotyping-based MTBC genotypic diversity in Malaysia followed by a comparison of strains with those prevailing in neighboring countries by comparison with an international MTBC genotyping database.

Methods

Spoligotyping was performed on a total of 220 M. tuberculosis clinical isolates collected in Kelantan and Kuala Lumpur. The results were compared with the SITVIT2 international database of the Pasteur Institute of Guadeloupe.

Results

Spoligotyping revealed 77 different patterns: 22 corresponded to orphan patterns while 55 patterns containing 198 isolates were assigned a Spoligo International Type (SIT) designation in the database (the latter included 6 newly created SITs). The eight most common SITs grouped 141 isolates (5 to 56 strains per cluster) as follows: SIT1/Beijing, n = 56, 25.5%; SIT745/EAI1-SOM, n = 33, 15.0%; SIT591/EAI6-BGD1, n = 13, 5.9%; SIT256/EAI5, n = 12, 5.5%; SIT236/EAI5, n = 10, 4.6%; SIT19/EAI2-Manila, n = 9, 4.1%; SIT89/EAI2-Nonthaburi, n = 5, 2.3%; and SIT50/H3, n = 3, 1.4%. The association between city of isolation and lineages was statistically significant; Haarlem and T lineages being higher in Kuala Lumpur (p<0.01). However, no statistically significant differences were noted when comparing drug resistance vs. major lineages, nor between gender and clades.

Conclusions

The ancestral East-African-Indian (EAI) lineage was most predominant followed by the Beijing lineage. A comparison of strains with those prevailing in neighboring countries in South Asia, East Asia and South East Asia underlined the phylogeographical specificity of SIT745 for Malaysia, and its probable ongoing evolution with locally evolved strains sharing a specific signature characterized by absence of spacers 37, 38, and 40. Pending complementary genotyping confirmation, we propose that SIT745/EAI-SOM is tentatively reclassified as SIT745/EAI-MYS.

Characterization of the genetic diversity of extensively-drug resistant Mycobacterium tuberculosis clinical isolates from pulmonary tuberculosis patients in Peru

Background

Peru holds the fourth highest burden of tuberculosis in the Americas. Despite an apparently well-functioning DOTS control program, the prevalence of multidrug resistant tuberculosis (MDR-TB) continues to increase. To worsen this situation, cases of extensively drug resistance tuberculosis (XDR-TB) have been detected. Little information exists about the genetic diversity of drug-susceptible vs. MDR-TB and XDR-TB.

Methods

Cryopreserved samples of XDR strains from 2007 to 2009 (second semester), were identified and collected. Starting from 227 frozen samples, a total of 142 XDR-TB strains of Mycobacterium tuberculosis complex (MTBC; 1 isolate per patient) were retained for this study. Each strain DNA was analyzed by spoligotyping and the 15-loci Mycobacterial Interspersed Repetitive Unit (MIRU-15).

Results

Among the 142 isolates analyzed, only 2 samples (1.41%) could not be matched to any lineage. The most prevalent sublineage was Haarlem (43.66%), followed by T (27.46%), LAM (16.2%), Beijing (9.15%), and X clade (1.41%). Spoligotype analysis identified clustering for 128/142 (90.1%) isolates vs. 49/142 (34.5%) with MIRUs. Of the samples, 90.85% belonged to retreated patients. The drug resistant profile demonstrated that 62.67% showed resistance to injectable drugs capreomycin (CAP) and kanamycin (KAN) vs. 15.5% to CAP alone and 21.8% to KAN alone. The SIT219/T1 and SIT50/H3 were the most prevalent patterns in our study. The spoligoforest analysis showed that SIT53/T1 was at the origin of many of the T lineage strains as well as a big proportion of Haarlem lineage strains (SIT50/H3, followed by SIT47/H1, SIT49/H3, and SIT2375/H1), as opposed to the SIT1/Beijing strains that did not appear to evolve into minor Beijing sublineages among the XDR-TB strains.

Conclusion

In contrast with other Latin-American countries where LAM sublineage is the most predominant, we found the Haarlem to be the most common followed by T sublineage among the XDR-TB strains.

Strain classification of Mycobacterium tuberculosis isolates in Brazil based on genotypes obtained by spoligotyping, mycobacterial interspersed repetitive unit typing and the presence of large sequence and single nucleotide polymorphism

Rio de Janeiro is endemic for tuberculosis (TB) and presents the second largest prevalence of the disease in Brazil. Here, we present the bacterial population structure of 218 isolates of Mycobacterium tuberculosis, derived from 186 patients that were diagnosed between January 2008 and December 2009. Genotypes were generated by means of spoligotyping, 24 MIRU-VNTR typing and presence of fbpC¹⁰³, RD^Rio and RD174. The results confirmed earlier data that predominant genotypes in Rio de Janeiro are those of the Euro American Lineages (99%). However, we observed differences between the classification by spoligotyping when comparing to that of 24 MIRU-VNTR typing, being respectively 43.6% vs. 62.4% of LAM, 34.9% vs. 9.6% of T and 18.3% vs. 21.5% of Haarlem. Among isolates classified as LAM by MIRU typing, 28.0% did not present the characteristic spoligotype profile with absence of spacers 21 to 24 and 32 to 36 and we designated these conveniently as “LAM-like”, 79.3% of these presenting the LAM-specific SNP fbpC¹⁰³. The frequency of RD^Rio and RD174 in the LAM strains, as defined both by spoligotyping and 24 MIRU-VNTR loci, were respectively 11% and 15.4%, demonstrating that RD174 is not always a marker for LAM/RD^Rio strains. We conclude that, although spoligotyping alone is a tool for classification of strains of the Euro-American lineage, when combined with MIRU-VNTRs, SNPs and RD typing, it leads to a much better understanding of the bacterial population structure and phylogenetic relationships among strains of M. tuberculosis in regions with high incidence of TB.

Population structure among mycobacterium tuberculosis isolates from pulmonary tuberculosis patients in Colombia

Background

Phylogeographic composition of M. tuberculosis populations reveals associations between lineages and human populations that might have implications for the development of strategies to control the disease. In Latin America, lineage 4 or the Euro-American, is predominant with considerable variations among and within countries. In Colombia, although few studies from specific localities have revealed differences in M. tuberculosis populations, there are still areas of the country where this information is lacking, as is a comparison of Colombian isolates with those from the rest of the world.

Principal Findings

A total of 414 M. tuberculosis isolates from adult pulmonary tuberculosis cases from three Colombian states were studied. Isolates were genotyped using IS6110-restriction fragment length polymorphism (RFLP), spoligotyping, and 24-locus Mycobacterial interspersed repetitive units variable number tandem repeats (MIRU-VNTRs). SIT42 (LAM9) and SIT62 (H1) represented 53.3% of isolates, followed by 8.21% SIT50 (H3), 5.07% SIT53 (T1), and 3.14% SIT727 (H1). Composite spoligotyping and 24-locus MIRU- VNTR minimum spanning tree analysis suggest a recent expansion of SIT42 and SIT62 evolved originally from SIT53 (T1). The proportion of Haarlem sublineage (44.3%) was significantly higher than that in neighboring countries. Associations were found between M. tuberculosis MDR and SIT45 (H1), as well as HIV-positive serology with SIT727 (H1) and SIT53 (T1).

Conclusions

This study showed the population structure of M. tuberculosis in several regions from Colombia with a dominance of the LAM and Haarlem sublineages, particularly in two major urban settings (Medellín and Cali). Dominant spoligotypes were LAM9 (SIT 42) and Haarlem (SIT62). The proportion of the Haarlem sublineage was higher in Colombia compared to that in neighboring countries, suggesting particular conditions of co-evolution with the corresponding human population that favor the success of this sublineage.

Molecular epidemiology of tuberculosis in Finland, 2008-2011

In industrialized countries the majority of tuberculosis (TB) cases are linked to immigration. In Finland, most cases are still Finnish born but the number of foreign born cases is steadily increasing. In this 4-year population based study, the TB situation in Finland was characterized by a genotypic analysis of Mycobacterium tuberculosis isolates. A total of 1048 M. tuberculosis isolates (representing 99.4% of all culture positive cases) were analyzed by spoligotyping and MIRU. Spoligotype lineages belonging to the Euro-American family were predominant among the Finnish isolates, particularly T (n=346, 33.0%) and Haarlem (n=237, 22.6%) strains. The lineage signature was unknown for 130 (12.4%) isolates. Out of the 17 multi-drug resistant TB strains, 10 (58.8%) belonged to the Beijing lineage. In total, 23 new SIT designations were given and 51 orphan strains were found, of which 58 patterns were unique to Finland. Phylogeographical TB mapping as compared to neighboring countries showed that the population structure in Finland most closely resembled that observed in Sweden. By combining spoligotyping and MIRU results, 98 clusters comprising 355 isolates (33.9%) were found. Only 10 clusters contained both Finnish and foreign born cases. In conclusion, a large proportion of the M. tuberculosis isolates were from Finnish born elderly patients. Moreover, many previously unidentified spoligotype profiles and isolates belonging to unknown lineages were encountered.