Genetic characterization of the Guinea-Bissau family of Mycobacterium tuberculosis complex strains

Prediction of drug resistance by Sanger sequencing of Mycobacterium tuberculosis complex strains isolated from multidrug resistant tuberculosis suspect patients in Ethiopia

Background

Ethiopia is one of the high multidrug-resistant tuberculosis (MDR-TB) burden countries. However, phenotypic drug susceptibility testing can take several weeks due to the slow growth of Mycobacterium tuberculosis complex (MTBC) strains. In this study, we assessed the performance of a Sanger sequencing approach to predict resistance against five anti-tuberculosis drugs and the pattern of resistance mediating mutations.

Methods

We enrolled 226 MTBC culture-positive MDR-TB suspects and collected sputum specimens and socio-demographic and TB related data from each suspect between June 2015 and December 2016 in Addis Ababa, Ethiopia. Phenotypic drug susceptibility testing (pDST) for rifampicin, isoniazid, pyrazinamide, ethambutol, and streptomycin using BACTEC MGIT 960 was compared with the results of a Sanger sequencing analysis of seven resistance determining regions in the genes rpoB, katG, fabG-inhA, pncA, embB, rpsL, and rrs.

Result

DNA isolation for Sanger sequencing was successfully extracted from 92.5% (209/226) of the MTBC positive cultures, and the remaining 7.5% (17/226) strains were excluded from the final analysis. Based on pDST results, drug resistance proportions were as follows: isoniazid: 109/209 (52.2%), streptomycin: 93/209 (44.5%), rifampicin: 88/209 (42.1%), ethambutol: 74/209 (35.4%), and pyrazinamide: 69/209 (33.0%). Resistance against isoniazid was mainly mediated by the mutation katG S315T (97/209, 46.4%) and resistance against rifampicin by rpoB S531L (58/209, 27.8%). The dominating resistance-conferring mutations for ethambutol, streptomycin, and pyrazinamide affected codon 306 in embB (48/209, 21.1%), codon 88 in rpsL (43/209, 20.6%), and codon 65 in pncA (19/209, 9.1%), respectively. We observed a high agreement between phenotypic and genotypic DST, such as 89.9% (at 95% confidence interval [CI], 84.2%–95.8%) for isoniazid, 95.5% (95% CI, 91.2%–99.8%) for rifampicin, 98.6% (95% CI, 95.9–100%) for ethambutol, 91.3% (95% CI, 84.6–98.1%) for pyrazinamide and 57.0% (95% CI, 46.9%–67.1%) for streptomycin.

Conclusion

We detected canonical mutations implicated in resistance to rifampicin, isoniazid, pyrazinamide, ethambutol, and streptomycin. High agreement with phenotypic DST results for all drugs renders Sanger sequencing promising to be performed as a complementary measure to routine phenotypic DST in Ethiopia. Sanger sequencing directly from sputum may accelerate accurate clinical decision-making in the future.

Genetic diversity of Mycobacterium tuberculosis strains circulating in Botswana

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

High prevalence of shared international type 53 among Mycobacterium tuberculosis complex strains in retreated patients from Côte d'Ivoire

BACKGROUND

Genotyping methods are useful tools to provide information on tuberculosis epidemic. They can allow a better response from health authorities and the implementation of measures for tuberculosis control. This study aimed to identify the main lineages and clades of Mycobacterium tuberculosis complex strains circulating in Côte d'Ivoire.

METHODS/MAIN FINDINGS

Strains isolated from sputum samples of patients ongoing retreatment from all the country were characterized by spoligotyping and by MIRU-VNTR. Profiles obtained by spoligotyping were first compared to the SITVIT/SpolDB4 database for family assignment. Of 194 strains analysed, 146 (75.3%) belonged to the T lineage. The most predominant spoligotype was the shared international type 53 with 135 strains (69.6%). In contrast with neighbouring countries, LAM (11 strains, 5.7%) and H (9 strains 4.6%) lineages were slightly represented. Only 3 Beijing strains (1.5%) and 4 strains of Mycobacterium africanum (2%) were found. Analysis of the results obtained with MIRU-VNTR revealed also a high level of clustering.

CONCLUSION/SIGNIFICANCE

The population of Mycobacterium tuberculosis complex strains among retreatment cases in Côte d'Ivoire exhibits a low diversity, allowing to assume recent transmission and locally based infection.

The genome of Mycobacterium africanum West African 2 reveals a lineage-specific locus and genome erosion common to the M. tuberculosis complex

Background

M. africanum West African 2 constitutes an ancient lineage of the M. tuberculosis complex that commonly causes human tuberculosis in West Africa and has an attenuated phenotype relative to M. tuberculosis.

Methodology/Principal Findings

In search of candidate genes underlying these differences, the genome of M. africanum West African 2 was sequenced using classical capillary sequencing techniques. Our findings reveal a unique sequence, RD900, that was independently lost during the evolution of two important lineages within the complex: the “modern” M. tuberculosis group and the lineage leading to M. bovis. Closely related to M. bovis and other animal strains within the M. tuberculosis complex, M. africanum West African 2 shares an abundance of pseudogenes with M. bovis but also with M. africanum West African clade 1. Comparison with other strains of the M. tuberculosis complex revealed pseudogenes events in all the known lineages pointing toward ongoing genome erosion likely due to increased genetic drift and relaxed selection linked to serial transmission-bottlenecks and an intracellular lifestyle.

Conclusions/Significance

The genomic differences identified between M. africanum West African 2 and the other strains of the Mycobacterium tuberculosis complex may explain its attenuated phenotype, and pave the way for targeted experiments to elucidate the phenotypic characteristic of M. africanum. Moreover, availability of the whole genome data allows for verification of conservation of targets used for the next generation of diagnostics and vaccines, in order to ensure similar efficacy in West Africa.

Mycobacterium africanum, a close relative of M. tuberculosis, is studied for the following reasons: M. africanum is commonly isolated from West African patients with tuberculosis yet has not spread beyond this region, it is more common in HIV infected patients, and it is less likely to lead to tuberculosis after one is exposed to an infectious case. Understanding this organism's unique biology gets a boost from the decoding of its genome, reported in this issue. For example, genome analysis reveals that M. africanum contains a region shared with “ancient” lineages in the M. tuberculosis complex and other mycobacterial species, which was lost independently from both M. tuberculosis and M. bovis. This region encodes a protein involved in transmembrane transport. Furthermore, M. africanum has lost genes, including a known virulence gene and genes for vitamin synthesis, in addition to an intact copy of a gene that may increase its susceptibility to antibiotics that are insufficiently active against M. tuberculosis. Finally, the genome sequence and analysis reported here will aid in the development of new diagnostics and vaccines against tuberculosis, which need to take into account the differences between M. africanum and other species in order to be effective worldwide.

The Guinea-Bissau family of Mycobacterium tuberculosis complex revisited

The Guinea-Bissau family of strains is a unique group of the Mycobacterium tuberculosis complex that, although genotypically closely related, phenotypically demonstrates considerable heterogeneity. We have investigated 414 M. tuberculosis complex strains collected in Guinea-Bissau between 1989 and 2008 in order to further characterize the Guinea-Bissau family of strains. To determine the strain lineages present in the study sample, binary outcomes of spoligotyping were compared with spoligotypes existing in the international database SITVIT2. The major circulating M. tuberculosis clades ranked in the following order: AFRI (n = 195, 47.10%), Latin-American-Mediterranean (LAM) (n = 75, 18.12%), ill-defined T clade (n = 53, 12.8%), Haarlem (n = 37, 8.85%), East-African-Indian (EAI) (n = 25, 6.04%), Unknown (n = 12, 2.87%), Beijing (n = 7, 1.68%), X clade (n = 4, 0.96%), Manu (n = 4, 0.97%), CAS (n = 2, 0.48%). Two strains of the LAM clade isolated in 2007 belonged to the Cameroon family (SIT61). All AFRI isolates except one belonged to the Guinea-Bissau family, i.e. they have an AFRI_1 spoligotype pattern, they have a distinct RFLP pattern with low numbers of IS6110 insertions, and they lack the regions of difference RD7, RD8, RD9 and RD10, RD701 and RD702. This profile classifies the Guinea-Bissau family, irrespective of phenotypic biovar, as part of the M. africanum West African 2 lineage, or the AFRI_1 sublineage according to the spoligtyping nomenclature. Guinea-Bissau family strains display a variation of biochemical traits classically used to differentiate M. tuberculosis from M. bovis. Yet, the differential expression of these biochemical traits was not related to any genes so far investigated (narGHJI and pncA). Guinea-Bissau has the highest prevalence of M. africanum recorded in the African continent, and the Guinea-Bissau family shows a high phylogeographical specificity for Western Africa, with Guinea-Bissau being the epicenter. Trends over time however indicate that this family of strains is waning in most parts of Western Africa, including Guinea-Bissau (p = 0.048).

Molecular diversity of Mycobacterium tuberculosis isolates from patients with pulmonary tuberculosis in Mozambique

BACKGROUND

Mozambique is one of the countries with the highest burden of tuberculosis (TB) in Sub-Saharan Africa, and information on the predominant genotypes of Mycobacterium tuberculosis circulating in the country are important to better understand the epidemic. This study determined the predominant strain lineages that cause TB in Mozambique.

RESULTS

A total of 445 M. tuberculosis isolates from seven different provinces of Mozambique were characterized by spoligotyping and resulting profiles were compared with the international spoligotyping database SITVIT2.The four most predominant lineages observed were: the Latin-American Mediterranean (LAM, n = 165 or 37%); the East African-Indian (EAI, n = 132 or 29.7%); an evolutionary recent but yet ill-defined T clade, (n = 52 or 11.6%); and the globally-emerging Beijing clone, (n = 31 or 7%). A high spoligotype diversity was found for the EAI, LAM and T lineages.

CONCLUSIONS

The TB epidemic in Mozambique is caused by a wide diversity of spoligotypes with predominance of LAM, EAI, T and Beijing lineages.

Distinct genotypic profiles of the two major clades of Mycobacterium africanum

Background

Mycobacterium tuberculosis is the principal etiologic agent of human tuberculosis (TB) and a member of the M. tuberculosis complex (MTC). Additional MTC species that cause TB in humans and other mammals include Mycobacterium africanum and Mycobacterium bovis. One result of studies interrogating recently identified MTC phylogenetic markers has been the recognition of at least two distinct lineages of M. africanum, known as West African-1 and West African-2.

Methods

We screened a blinded non-random set of MTC strains isolated from TB patients in Ghana (n = 47) for known chromosomal region-of-difference (RD) loci and single nucleotide polymorphisms (SNPs). A MTC PCR-typing panel, single-target standard PCR, multi-primer PCR, PCR-restriction fragment analysis, and sequence analysis of amplified products were among the methods utilized for the comparative evaluation of targets and identification systems. The MTC distributions of novel SNPs were characterized in the both the Ghana collection and two other diverse collections of MTC strains (n = 175 in total).

Results

The utility of various polymorphisms as species-, lineage-, and sublineage-defining phylogenetic markers for M. africanum was determined. Novel SNPs were also identified and found to be specific to either M. africanum West African-1 (Rv1332⁵²³; n = 32) or M. africanum West African-2 (nat⁷⁵¹; n = 27). In the final analysis, a strain identification approach that combined multi-primer PCR targeting of the RD loci RD9, RD10, and RD702 was the most simple, straight-forward, and definitive means of distinguishing the two clades of M. africanum from one another and from other MTC species.

Conclusion

With this study, we have organized a series of consistent phylogenetically-relevant markers for each of the distinct MTC lineages that share the M. africanum designation. A differential distribution of each M. africanum clade in Western Africa is described.

Novel genetic polymorphisms that further delineate the phylogeny of the Mycobacterium tuberculosis complex

In a previous report, we described a PCR protocol for the differentiation of the various species of the Mycobacterium tuberculosis complex (MTC) on the basis of genomic deletions (R. C. Huard, L. C. de Oliveira Lazzarini, W. R. Butler, D. van Soolingen, and J. L. Ho, J. Clin. Microbiol. 41:1637-1650, 2003). That report also provided a broad cross-comparison of several previously identified, phylogenetically relevant, long-sequence and single-nucleotide polymorphisms (LSPs and SNPs, respectively). In the present companion report, we expand upon the previous work (i) by continuing the evaluation of known MTC phylogenetic markers in a larger collection of tubercle bacilli (n = 125), (ii) by evaluating additional recently reported MTC species-specific and interspecific polymorphisms, and (iii) by describing the identification and distribution of a number of novel LSPs and SNPs. Notably, new genomic deletions were found in various Mycobacterium tuberculosis strains, new species-specific SNPs were identified for "Mycobacterium canettii," Mycobacterium microti, and Mycobacterium pinnipedii, and, for the first time, intraspecific single-nucleotide DNA differences were discovered for the dassie bacillus, the oryx bacillus, and the two Mycobacterium africanum subtype I variants. Surprisingly, coincident polymorphisms linked one M. africanum subtype I genotype with the dassie bacillus and M. microti with M. pinnipedii, thereby suggesting closer evolutionary ties within each pair of species than had been previously thought. Overall, the presented data add to the genetic definitions of several MTC organisms as well as fine-tune current models for the evolutionary history of the MTC.

Genomic analysis distinguishes Mycobacterium africanum

Mycobacterium africanum is thought to comprise a unique species within the Mycobacterium tuberculosis complex. M. africanum has traditionally been identified by phenotypic criteria, occupying an intermediate position between M. tuberculosis and M. bovis according to biochemical characteristics. Although M. africanum isolates present near-identical sequence homology to other species of the M. tuberculosis complex, several studies have uncovered large genomic regions variably deleted from certain M. africanum isolates. To further investigate the genomic characteristics of organisms characterized as M. africanum, the DNA content of 12 isolates was interrogated by using Affymetrix GeneChip. Analysis revealed genomic regions of M. tuberculosis deleted from all isolates of putative diagnostic and biological consequence. The distribution of deleted sequences suggests that M. africanum subtype II isolates are situated among strains of "modern" M. tuberculosis. In contrast, other M. africanum isolates (subtype I) constitute two distinct evolutionary branches within the M. tuberculosis complex. To test for an association between deleted sequences and biochemical attributes used for speciation, a phenotypically diverse panel of "M. africanum-like" isolates from Guinea-Bissau was tested for these deletions. These isolates clustered together within one of the M. africanum subtype I branches, irrespective of phenotype. These results indicate that convergent biochemical profiles can be independently obtained for M. tuberculosis complex members, challenging the traditional approach to M. tuberculosis complex speciation. Furthermore, the genomic results suggest a rational framework for defining M. africanum and provide tools to accurately assess its prevalence in clinical specimens.