Genetic diversity and population structure of Mycobacterium tuberculosis in Casablanca, a Moroccan city with high incidence of tuberculosis

Mutations Associated with Rifampicin Resistance in Mycobacterium tuberculosis Isolates from Moroccan Patients: Systematic Review

Background

In recent years, the treatment of tuberculosis has been threatened by the increasing number of patients with drug resistance, especially rifampicin resistance, which is the most effective first-line antibiotic against Mycobacterium tuberculosis.

Methods

We performed a systematic review of the literature by searching the PubMed database for studies of rifampicin-resistant Mycobacterium tuberculosis (MTB) isolates from Moroccan patients, published between 2010 and 2020. The aim of this review was to quantify the frequency of the most common mutations associated with rifampicin resistance, to describe the frequency at which these mutations co-occur. Identified studies were critically appraised according to the Quality Assessment of Diagnostic Accuracy Studies-2 (QUADAS-2) tool.

Results

6 studies met our inclusion criteria. Results show that 99.36% of MTB isolates had a single-point mutation, and the most commonly mutated codon of rpoB gene is 531 with 70.33% of phenotypically resistant strains. However, 10.38% of MTB strains phenotypically resistant to RIF did not exhibit any mutation in the rpoB gene.

Conclusion

Identification of a resistance-associated mutation to rifampicin can be a good marker of drug-resistant TB, but lack of a mutation in the target sequence must be interpreted with caution.

Molecular Typing of Mycobacterium Tuberculosis Complex by 24-Locus Based MIRU-VNTR Typing in Conjunction with Spoligotyping to Assess Genetic Diversity of Strains Circulating in Morocco

Background

Standard 24-locus Mycobacterial Interspersed Repetitive Unit Variable Number Tandem Repeat (MIRU-VNTR) typing allows to get an improved resolution power for tracing TB transmission and predicting different strain (sub) lineages in a community.

Methodology

During 2010–2012, a total of 168 Mycobacterium tuberculosis Complex (MTBC) isolates were collected by cluster sampling from 10 different Moroccan cities, and centralized by the National Reference Laboratory of Tuberculosis over the study period. All isolates were genotyped using spoligotyping, and a subset of 75 was genotyped using 24-locus based MIRU-VNTR typing, followed by first line drug susceptibility testing. Corresponding strain lineages were predicted using MIRU-VNTRplus database.

Principal Findings

Spoligotyping resulted in 137 isolates in 18 clusters (2–50 isolates per cluster: clustering rate of 81.54%) corresponding to a SIT number in the SITVIT database, while 31(18.45%) patterns were unique of which 10 were labelled as “unknown” according to the same database. The most prevalent spoligotype family was LAM; (n = 81 or 48.24% of isolates, dominated by SIT42, n = 49), followed by Haarlem (23.80%), T superfamily (15.47%), >Beijing (2.97%), > U clade (2.38%) and S clade (1.19%). Subsequent 24-Locus MIRU-VNTR typing identified 64 unique types and 11 isolates in 5 clusters (2 to 3isolates per cluster), substantially reducing clusters defined by spoligotyping only. The single cluster of three isolates corresponded to two previously treated MDR-TB cases and one new MDR-TB case known to be contact a same index case and belonging to a same family, albeit residing in 3 different administrative regions. MIRU-VNTR loci 4052, 802, 2996, 2163b, 3690, 1955, 424, 2531, 2401 and 960 were highly discriminative in our setting (HGDI >0.6).

Conclusions

24-locus MIRU-VNTR typing can substantially improve the resolution of large clusters initially defined by spoligotyping alone and predominating in Morocco, and could therefore be used to better study tuberculosis transmission in a population-based, multi-year sample context.

Prevalence of Extensively Drug Resistant Tuberculosis among Archived Multidrug Resistant Tuberculosis Isolates in Zimbabwe

We conducted a cross-sectional study of second line drug resistance patterns and genetic diversity of MDR-TB isolates archived at the BRTI-TB Laboratory, Harare, between January 2007 and December 2011. DSTs were performed for second line antituberculosis drugs. XDR-TB strains were defined as MDR-TB strains with resistance to either kanamycin and ofloxacin or capreomycin and ofloxacin. Strain types were identified by spoligotyping. No resistance to any second line drugs was shown in 73% of the isolates, with 23% resistant to one or two drugs but not meeting the definition of XDR-TB. A total of 26 shared types were identified, and 18 (69%) matched preexisting shared types in the current published spoligotype databases. Of the 11 out of 18 clustered SITs, 4 predominant (>6 isolates per shared type) were identified. The most and least abundant types were SIT 1468 (LAM 11-ZWE) with 12 (18%) isolates and SIT 53 (T1) with 6 (9%) isolates, respectively. XDR-TB strains are rare in Zimbabwe, but the high proportion of “pre-XDR-TB” strains and treatment failure cases is of concern. The genetic diversity of the MDR-TB strains showed no significant association between SITs and drug resistance.

Contribution of spoligotyping and MIRU-VNTRs to characterize prevalent Mycobacterium tuberculosis genotypes infecting tuberculosis patients in Morocco

In the present study, Mycobacterium tuberculosis complex (MTBC) clinical isolates from culture-positive TB patients in Morocco were studied by spoligotyping and 12-loci MIRU-VNTR typing methods to characterize prevalent genotypes (n = 219 isolates from 208 patients). Spoligotyping resulted in 39 unique patterns and 167 strains in 30 clusters (2-50 strains per cluster). Comparison with international database showed that 29 of 39 unique patterns matched existing shared spoligotype international types (SITs). Nine shared types containing 10 strains were newly created (SIT 2891 to SIT 2899); this led to the description of 69 SITs with 206 strains and two orphan patterns. The most prevalent spoligotype was SIT42 (LAM; n = 50 or 24% of isolates). The repartition of strains according to major MTBC clades was as follows LAM (46.1%)> Haarlem (26%) >ill-defined T superfamily (22.6%) and S clade (0.96%). On the other hand, Beijing, CAS (Central Asian) and EAI (East-African Indian) strains were absent in this setting. Subsequent 12-Loci MIRU typing resulted in a total of 25 SIT/MIT clusters (n = 66 isolates, 2-6 isolates per cluster), with a resulting recent transmission rate of 22.3%. The MIRU-VNTR patterns corresponded to 69 MITs for 138 strains and 46 orphan patterns. The most frequent patterns were MIT43 (n = 8), MIT9 (n = 7) and MIT42 (n = 7). HGDI analysis of the 12 MIRU loci showed that loci 10, 23 and 40 were highly discriminative in our setting. The results also underlined the usefulness of spoligotyping and MIRU-VNTR to detect mixed infections among certain of our TB patients. Globally, the results obtained showed that TB is almost exclusively transmitted in Morocco through evolutionary-modern MTBC lineages belonging to principal genetic groups 2/3 strains (Haarlem, LAM, T), with a high level of biodiversity seen by MIRU typing. This study provides with a 1st global snapshot of MTBC population structure in Morocco, and validates the potential use of spoligotyping in conjunction with minisatellites for future investigations in Morocco that should in future ideally include optimized 15- or 24-loci MIRU-VNTRs.

A model-based Bayesian estimation of the rate of evolution of VNTR loci in Mycobacterium tuberculosis

Variable numbers of tandem repeats (VNTR) typing is widely used for studying the bacterial cause of tuberculosis. Knowledge of the rate of mutation of VNTR loci facilitates the study of the evolution and epidemiology of Mycobacterium tuberculosis. Previous studies have applied population genetic models to estimate the mutation rate, leading to estimates varying widely from around to per locus per year. Resolving this issue using more detailed models and statistical methods would lead to improved inference in the molecular epidemiology of tuberculosis. Here, we use a model-based approach that incorporates two alternative forms of a stepwise mutation process for VNTR evolution within an epidemiological model of disease transmission. Using this model in a Bayesian framework we estimate the mutation rate of VNTR in M. tuberculosis from four published data sets of VNTR profiles from Albania, Iran, Morocco and Venezuela. In the first variant, the mutation rate increases linearly with respect to repeat numbers (linear model); in the second, the mutation rate is constant across repeat numbers (constant model). We find that under the constant model, the mean mutation rate per locus is (95% CI: ,)and under the linear model, the mean mutation rate per locus per repeat unit is (95% CI: ,). These new estimates represent a high rate of mutation at VNTR loci compared to previous estimates. To compare the two models we use posterior predictive checks to ascertain which of the two models is better able to reproduce the observed data. From this procedure we find that the linear model performs better than the constant model. The general framework we use allows the possibility of extending the analysis to more complex models in the future.

Genetically typing the bacterium responsible for tuberculosis is useful for understanding the evolutionary and epidemiological characteristics of the disease. Typing methods based on variable number tandem repeat (VNTR) loci are increasingly being used. These loci, which are composed of repeated units, mutate by increasing or decreasing in the number of these repeats. Knowledge of the mutation rate of molecular markers facilitates the epidemiological interpretation of the observed genetic variation in a sample of bacterial isolates. Few studies have examined the rate of mutation at these markers and estimates to date have varied considerably. To address this problem we develop a stochastic model of evolution of these markers and then estimate their mutation rate using approximate Bayesian computation. We examine two alternative forms of the mutation process. The observed data are from four published data sets of tuberculosis bacterial isolates sampled in Albania, Iran, Morocco and Venezuela. We find that these markers have fairly high rates of mutation compared with estimates from previous studies.

Highly structured genetic diversity of the Mycobacterium tuberculosis population in Djibouti

Djibouti is an East African country with a high tuberculosis incidence. This study was conducted over a 2-month period in Djibouti, during which 62 consecutive patients with pulmonary tuberculosis (TB) were included. Genetic characterization of Mycobacterium tuberculosis, using mycobacterial interspersed repetitive-unit variable-number tandem-repeat typing and spoligotyping, was performed. The genetic and phylogenetic analysis revealed only three major families (Central Asian, East African Indian and T). The high diversity and linkage disequilibrium within each family suggest a long period of clonal evolution. A Bayesian approach shows that the phylogenetic structure observed in our sample of 62 isolates is very likely to be representative of the phylogenetic structure of the M. tuberculosis population in the total number of TB cases.

Genetic diversity among Mycobacterium tuberculosis isolates from Mexican patients

Forty-six isolates of the Mycobacterium tuberculosis complex were typified by PCR of the IS6110 region and by Mycobacterium bovis specific primers JB21/JB22. Isolate MVG01 was typified as M. bovis, being the first record of a case of human tuberculosis caused by this species in Mexico. RAPD-PCR was used to describe the genetic diversity of the remaining 45 M. tuberculosis complex isolates. The corrected genotypic diversity value calculated for the analyzed population was 0.96, the estimated mean gene diversity was 0.235, and the corrected Shannon-Weiner index was 2.15. All allele-loci combinations generated showed significant linkage disequilibria. The distribution of genetic variation was analyzed both by the unweighted pair group method with arithmetic averages clustering and by principal coordinates analysis. Unweighted pair group method with arithmetic averages clustering resulted in a tree with four main clusters and one unclustered strain (MVG20), the principal coordinates analysis strain distribution pattern being consistent with this grouping. The obtained results suggest that the studied isolates belong to a clonal population having significant genetic diversity. Our genetic diversity results are comparable with those reported for other populations of M. tuberculosis, although only three RAPD primers were used.

Characterization of Mycobacterium tuberculosis Central Asian Strain 1 using mycobacterial interspersed repetitive unit genotyping

BACKGROUND

The Central Asian Strain 1 (CAS1) genogroup of Mycobacterium tuberculosis (MTB) is the most prevalent in Pakistan, India and Bangladesh. Mycobacterial interspersed repetitive units variable number tandem repeat (MIRU-VNTR) typing is a reliable and reproducible method for differentiation of MTB isolates. However, information of its utility in determining the diversity of CAS1 strain is limited. We performed standard 12 loci based MIRU-VNTR typing on previously spoligotyped CAS1 strains and 'unique' strains in order to evaluate its discriminatory power for these isolates.

METHODS

Twelve loci based MIRU- VNTR typing was used to type 178 CAS1 and 189 'unique' MTB strains. The discriminatory index for each of the loci was calculated using the Hunter Gaston Discriminatory Index (HGDI). A subset of these strains (n = 78) were typed using IS6110 restriction fragment length polymorphism (RFLP). MIRU-VNTR profiles were studied together with their drug susceptibility patterns.

RESULTS

A total of 349 MIRU patterns were obtained for the 367 strains tested. The CAS1 strains were subdivided into 160 distinct patterns; 15 clusters of 2 strains each, 1 cluster of four strains and 144 unique patterns. Using HGDI, seven MIRU loci, (numbers 26, 31, 27, 16, 10, 39, and 40) were found to be "highly discriminatory" (DI: >or=0.6), four MIRU loci (numbers 20, 24, 23, and 4) were "moderately discriminatory" (DI: 0.3-0.59), and one locus (number 2) was "poorly discriminatory" (DI< 0.3). Loci 26 and 31 were the most discriminatory for the CAS1 isolates. Amongst 'unique' strains in addition to loci 26, 31, 27, 16, 10, 39, and 40, locus 23 was highly discriminatory, while no locus was poorly discriminating. DI values for loci 4, 10 and 26 were significantly lower (P-value < .01) in CAS1 strains than in 'unique' strains. The association between CAS1 strains and MDR was not found to be significant (p value = 0.21).

CONCLUSION

We propose that MIRU typing could be used to estimate the phylogenetic relatedness amongst prevalent CAS1 strains, for which MIRU loci 26, 31, 16, 10, 27, 39 and 40 were found to be the most discriminatory.

Tuberculosis transmission in a high incidence area: a retrospective molecular epidemiological study of Mycobacterium tuberculosis in Casablanca, Morocco

Like in most developing countries, tuberculosis represents a major public health problem in Morocco. This paper describes the first study combining molecular and conventional epidemiology of tuberculosis in Casablanca, the economic capital of this country. Molecular fingerprinting of the genomic DNA recovered from cultures of sputum of 150 patients was performed by MIRU-VNTR. This molecular marker revealed that 53.1% of the total cases were clustered. These cases were classified into 23 clusters ranging in size from 2 to 13 patients, suggesting a rate of 37% of recent transmission in the sample under study. In a multivariate analysis, there were no independent predictors of clustering. However, the clinical form was associated with drug resistance (odds ratio=9.9; P value=0.0006). The phylogenetic analysis showed that the heterogeneity found in this population includes also the members from a same patient family, and that the 2 major families distributed in Casablanca were the Latin-American-Mediterranean (LAM) and Haarlem families. All the results of this work allow to understand better the tuberculosis transmission in Casablanca, and suggest that different clones of M. tuberculosis seem to circulate in this city, and that the reactivation of latent infections would be mainly responsible for the endemic situation of this disease. These findings indicate also that the transmission of TB in Morocco is not optimally controlled, and that efforts for control strategies should be sustained in all developing countries where the incidence of TB is high and still raising.

Random amplified polymorphic DNA (RAPD) analysis of Mycobacterium tuberculosis strains in India

The usefulness of random amplification of polymorphic DNA (RAPD) analysis for typing Indian strains of M. tuberculosis was investigated. M. tuberculosis H37Rv, M. tuberculosis DT and 42 clinical isolates of M. tuberculosis were subjected to RAPD-PCR using 7 random decamer primers. All 7 primers were found to be differentiated and produced specific RAPD profiles. The polymorphic amplicons served as RAPD markers for M. tuberculosis. The dendrograms, obtained by different primers, showed the discriminatory ability of the primers. RAPD analysis provided a rapid and easy means of identifying polymorphism in M. tuberculosis isolates, and it was found to be a valuable alternative epidemiological tool. In addition, the results of the present study showed heterogeneity in the M. tuberculosis strains in the population studied.