|Year : 2015 | Volume
| Issue : 3 | Page : 191-195
Prevalence of Latin-American-Mediterranean genetic family in population structure of Mycobacterium tuberculosis in Bulgaria
Violeta Valcheva1, Nalin Rastogi2, Igor Mokrousov3
1 Department of Infectious Microbiology, The Stephan Angeloff Institute of Microbiology, n Academy of Sciences, Sofia 1113, Bulgaria
2 WHO Supranational Tuberculosis Reference Laboratory, Tuberculosis and Mycobacteria Unit, Institut Pasteur de Guadeloupe, Abymes 97183, Guadeloupe, France
3 Laboratory of Molecular Microbiology, St. Petersburg Pasteur Institute, 197101 St. Petersburg, Russia
|Date of Web Publication||23-Feb-2017|
St. Petersburg Pasteur Institute, 14 Mira Street, 197101 St. Petersburg
The Stephan Angeloff Institute of Microbiology, 26 Acad. G. Bonchev Street, Sofia 1113
Source of Support: None, Conflict of Interest: None
Tuberculosis (TB) control remains an important public health priority for Bulgaria. The population structure of Mycobacterium tuberculosis is clonal and certain genetic families of this species (e.g., Latin-American-Mediterranean [LAM]) have attracted more attention due to their global dissemination and/or particular pathogenic properties, e.g., association with multidrug resistance (MDR). The aim of this study was to evaluate the prevalence of the M. tuberculosis LAM family in Bulgaria based on the use of different molecular markers. A total of 101 previously spoligotyped M. tuberculosis strains were studied by LAM-specific PCR assay to detect an insertion of IS6110 in the specific genome region. On the whole, clear-cut results were obtained for most strains; spoligotype-based family was reassigned in some of them. At the same time, double bands were amplified in some cases and warrant further validation studies of this method. The higher MDR rate among LAM versus other genotype isolates was observed (P = 0.04). In conclusion, these results suggest a low (<4%) prevalence rate of LAM in Bulgaria (that is similar to its Balkan neighbors) and highlight the importance of using robust markers for correct detection of the LAM family.
Keywords: Mycobacterium tuberculosis, LAM (Latin-American-Mediterranean) family, Genotyping, Spoligotyping, IS 6110, Bulgaria
|How to cite this article:|
Valcheva V, Rastogi N, Mokrousov I. Prevalence of Latin-American-Mediterranean genetic family in population structure of Mycobacterium tuberculosis in Bulgaria. Int J Mycobacteriol 2015;4:191-5
|How to cite this URL:|
Valcheva V, Rastogi N, Mokrousov I. Prevalence of Latin-American-Mediterranean genetic family in population structure of Mycobacterium tuberculosis in Bulgaria. Int J Mycobacteriol [serial online] 2015 [cited 2020 Jan 22];4:191-5. Available from: http://www.ijmyco.org/text.asp?2015/4/3/191/200821
| Introduction|| |
Tuberculosis (TB) remains an important public health issue in Bulgaria. Although the number of new cases has declined since 2006 (41/100,000), the TB incidence rate in Bulgaria is still sufficiently high (26.7/100,000 in 2013), as well as multidrug-resistant TB (MDR-TB) rates among new TB cases (5.8%). Clinical and epidemiological characterization of the major Mycobacterium tuberculosis (M. tuberculosis) subpopulations circulating in different geographic regions is an essential step for developing better diagnostics, therapeutics, and vaccines. The integration of molecular analysis of M. tuberculosis with an epidemiological analysis of clinical information provides a new tool to assess possible associations between M. tuberculosis strain types and the clinical and epidemiological characteristics of the disease .
The population structure of M. tuberculosis is clonal and certain genetic families of this species have attracted more attention due to their global dissemination and/or pathogenic properties. A number of studies demonstrated an association with drug resistance of the Latin-American-Mediterranean (LAM) family ,,. On the other hand, other studies have shown high transmissibility of this strain, but not a significant association with increased virulence ,. LAM was initially suggested by analysis of the global spoligotyping dataset  and its spoligotyping signature is absence of signals 21–24 and 33–36 ([Figure 1]). However, abridged spoligoprofiles with long blocks of deleted spacers (e.g., SIT125) or homoplasy spoligoprofiles (e.g., SIT803) make identification of such isolates uncertain ,,. To solve a phylogenetic ambiguity, the use of other molecular markers/methods may be helpful, in particular, MIRU-VNTR cluster analysis using reference LAM isolates (e.g., using MIRU-VNTRplus.org) or use of phylogenetically robust SNPs ,. An IS 6110 insertion in the specific genome region between genes lpqQ and Rv0836c (coding for hypothetical proteins) has also been suggested for detection of LAM , ([Figure 2]) and evaluated in the present study.
|Figure 1: Visual presentation of the decision rules for definition of the LAM and S spoligotype families and ambiguous spoligotypes.|
Click here to view
Correct detection of the LAM family of M. tuberculosis has both phylogenetic and practical interest. The aim of this study was to evaluate the prevalence of LAM among M. tuberculosis strains circulating in Bulgaria based on the use of different molecular markers.
| Materials and methods|| |
A study sample included 101 M. tuberculosis isolates obtained from different regions of the country and previously characterized by spoligotyping and other methods . The DNA of the studied strains was extracted from 4 to 6 weeks Löwenstein–Jensen medium culture using the recommended method . Spoligotyping was performed as described previously . The individual spoligotyping patterns were compared with the international database SITVIT2 (Institut Pasteur de Guadeloupe) that is an updated version of the published SITVIT_WEB database .
Analysis of the IS 6110 element specific for the LAM genetic family was performed as described previously . In brief, a 205-bp band indicates a LAM strain due to the presence of an IS 6110 element in a specific site in the genome, whereas a 141-bp band indicates a non-LAM strain lacking the IS 6110 element in this site ([Figure 3]).
|Figure 3: LAM-specific PCR: (a) schematic view of the genome region harboring LAM-specific IS6110 insertion in complete genome of strain F11. Three primers (arrows) are used in one reaction: LAM-F, LAM-R, XhoI. Non-LAM strain (without this IS6110-insertion): primers LAM-R and LAM-F amplify 141 bp fragment. LAM strain (with this IS6110 insertion): primers LAMR and XhoI amplify 205 bp fragment (b) gel-electrophoresis. Lane 4 – LAM strain; lanes 8–9: mixed profile. M, molecular weights marker “100 bp DNA ladder” (GE Healthcare).|
Click here to view
To minimize the risk of laboratory cross-contamination during PCR amplification, each procedure (preparation of the PCR mixes, the addition of the DNA, the PCR amplification, and the electrophoretic fractionation) was conducted in physically separated rooms. Negative controls (water) were included to control for reagent contamination.
A 2 × 2 χ2 test was used to detect any significant difference between the two groups. Yates corrected χ2 and P -values were calculated with 95% confidence interval at http://www.medcalc.org/calc/odds_ratio.php online resource.
| Results and discussion|| |
Application of the published rules for the definition of the major spoligotype clades and comparison with SITVIT2 global database permitted the assignment of most of the strains to the known spoligotype families. At the same time, some types were assigned to the convergent family (e.g., LAM/S) and use of further methods would be required. The PCR approach was used to define the LAM family  and all strains were tested for the presence or absence of the LAM-specific IS 6110 insertion ([Figure 3] ; [Table S1]). Comparison of results obtained by different methods is shown in Table S1 where the spoligoprofiles are sorted for easier visualization of data. This permitted the identification of 3 differently sized groups of strains. The first group included 11 strains with 2 amplified bands. The multiplex PCR targeted an insertion of IS 6110 in the specific genome region and was designed in such a way that either one of the other bands are amplified in case of the presence or absence of such insertion. Hence, these 11 strains present an apparent discrepancy. A double band is apriori impossible since this genome target appears to be unique. Only one of these strains was assigned to LAM in SITVIT_WEB (strain 54, SIT1280, although it was labeled as LAM/S in the previous version SpolDB4). Other strains with a double band belonged to T, X, Haarlem and unknown families; neither had a LAM-specific deletion block in the spoligoprofile (deleted signals 21–24). Furthermore, in 4 strains both bands were equally strong, thus excluding a possibility of contamination. 24-loci VNTR typing presented an additional proof of the lack of technical error (contamination): no strains had multiple alleles (not shown). An amplification of the additional band may be specific but not related to the targeted genome region. Further study of these cases is warranted.
The second group included 86 strains with a single amplified band, specific for the absence of the IS6110 insertion in the target locus, hence indicative of other than the LAM family (Table S1). Few observations can be made regarding the discrepancies with spoligotype-defined families. Four of the studied strains belonged to SIT41 defined in SITVIT_WEB as LAM7-TUR. However, no LAM-specific band was amplified in these strains; hence they cannot be defined as LAM. Indeed, another online resource for M. tuberculosis molecular typing (MIRU-VNTRplus.org) defines the family of these strains as TUR, based on the analysis of different molecular markers, including VNTRs and SNPs. It is also interesting to note that all studied strains of the SIT125, which had an uncertain family definition (LAM/S in earlier SpolDB4 and T2/LAM3 in SITVIT_WEB), did not have a LAM-specific band; hence their family status cannot be LAM. Finally, two strains of SIT1588 had a single non-LAM band which is intriguing, since these strains had a prototype LAM signature (see SIT42 in [Figure 1]), in addition to the deleted first signals. This profile may be a result of the convergent evolution of the DR locus, and more strains of this spoligotype SIT1588 should be studied in order to verify its family status.
Regarding the third and smallest group of strains, it included 4 strains with only LAM-specific band amplified. One of them–strain 52–was assigned to LAM in SITVIT_WEB; strain 8 was assigned to T5-RUS1, but it was recently shown to belong to LAM . In its turn, strain 126 (‘new’ spoligotype, ‘unknown’ family in SITVIT_WEB) may indeed belong to LAM as confirmed herein by LAM-PCR. At the same time, strain 81 (SIT453/T) had LAM-specific band amplified, but intact spacers 21–24; this discrepancy is unexpected and questions the value of this IS 6110 -PCR as LAM-specific.
The importance of the correct detection of LAM is explained by the fact that the heterogeneous genetic family of M. tuberculosis , LAM, has been shown in geographically distant settings to demonstrate remarkable pathogenic features. First, in Brazil, RD-Rio sublineage of LAM accounts for 37% of the total TB burden and was shown to be associated with pulmonary cavitation. Since cavitary TB is associated with a higher sputum-bacillary load, this finding supports the hypothesis that RD-Rio M. tuberculosis is associated with a more “severe” disease as a strategy to increase transmission, at least in some ethnic groups . Secondly, LAM-RUS sublineage in central Russia (along with the Beijing genotype) was shown to be associated with MDR and clustering: the level of drug resistance in new cases was almost twice as high as the estimated average national level . In the present study, 2 of the 4 isolates with a single LAM-specific PCR fragment were MDR, while two others were susceptible. Regarding 97 non-LAM isolates, 71 of them were susceptible and 10 were MDR. However, the higher MDR rate among LAM versus other strains (2/4 versus 10/97; P = 0.04) should be regarded with caution due to the very small sample size of the LAM group, and further study is warranted.
Geographically, the 4 confirmed LAM isolates were not evenly dispersed across Bulgaria, but rather represented northern (Vratsa, Veliko Tarnovo and Pleven) and western (Sofia) parts of the country. Further placing these results in the wider regional context permitted this study to reveal an intriguing gradient of the LAM prevalence. Data on some countries are absent or scarce and estimations based only on spoligotypes are uncertain. These estimations were obtained keeping in mind that: (i) SIT4 cannot be assigned to any family based on spoligotyping alone; and (ii) LAM-TUR and LAM-CAM do not belong to LAM, but rather represent separate families within the Euro-American lineage. It appears that similar to Bulgaria (4%), the LAM rate is low/negligible in other Balkan countries: Turkey – 3.4%; Greece – 1–4%; Romania – 0% (SITVIT_WEB); and Albania – 0% . In contrast, the LAM rate is higher in the northeastern neighbors of Bulgaria beyond the Balkans, reaching 17% in Southern Ukraine  , 25%-39% in European Russia, and 42% in Belarus . In its turn, this situation correlates with the prevalence of the rare Bulgarian LAM isolates mainly in the northern part of the country.
In conclusion, application of the LAM-specific PCR revealed a low (3%-4% depending on the marker used) prevalence rate of LAM family in the M. tuberculosis population in this country. Compared with spoligotyping, a phylogenetic family of 38 strains was revised or questioned: 27 strains were shown not to belong to LAM, while 11 more strains showed apparently discrepant results. This situation questions the global utility of such PCR and highlights the importance of using multiple markers for the molecular detection of the LAM family.
| Ethical approval|| |
| Funding|| |
| Conflict of interest|| |
N.R. and I.M. are Editorial Board members of the International Journal of Mycobacteriology.
| Acknowledgement|| |
Dr. Nadya Markova is gratefully acknowledged for her kind support.
| Appendix A. Supplementary data|| |
Supplementary data associated with this article can be found, in the online version, at http://dx.doi.org/10.1016/j.ijmyco.2015.04.003.
| References|| |
S. Gagneux, K. DeRiemer, T. Van, M. Kato-Maeda, B.C. de Jong, S. Narayanan, et al, Variable host-pathogen compatibility inMycobacterium tuberculosis
, Proc. Natl. Acad. Sci. U.S.A. 103 (2006) 2869–2873.
S. Dubiley, A. Ignatova, T. Mukhina, A. Nizova, S.Blagodatskikh, et al, Molecular epidemiology of tuberculosisin the Tula area, Central Russia, before the introduction of thedirectly observed therapy strategy, Clin. Microbiol. Infect. 16 (2010) 1421–1426.
E.R. Dalla Costa, C.L. Lazzarini, P.F. Perizzolo, C.A. Díaz, F.S.Spies, L.L. Costa, et al, Mycobacterium tuberculosis
of the RDRiogenotype is the predominant cause of tuberculosis andassociated with multidrug resistance in Porto Alegre City,South Brazil, J. Clin. Microbiol. 51 (2013) 1071–1077.
S. David, E.L. Duarte, C.Q. Leite, J.N. Ribeiro, J.N. Maio, E.Paixão, et al, Implication of the RD(Rio) Mycobacteriumtuberculosis sublineage in multidrug resistant tuberculosis inPortugal, Infect. Genet. Evol. 12 (2012) 1362–1367.
L.C. Lazzarini, S.M. Spindola, H. Bang, A.L. Gibson, S.Weisenberg, W. da Silva Carvalho, et al, RDRio Mycobacteriumtuberculosis infection is associated with a higher frequency ofcavitary pulmonary disease, J. Clin. Microbiol. 46 (2008) 2175–2183.
B.C. de Barbosa, L.C. Lazzarini, A.R. Elias, J.A. Leung, S.B.Ribeiro, M.G. da Silva, et al, Tuberculosis caused by RDRioMycobacterium tuberculosis
is not associated with differentialclinical features, Int. J. Tuberc. Lung Dis. 16 (2012) 1377–1382.
C. Sola, I. Filliol, E. Legrand, I. Mokrousov, N. Rastogi,Mycobacterium tuberculosis
phylogeny reconstruction based oncombined numerical analysis with IS1081, IS6110, VNTR, andDR-based spoligotyping suggests the existence of two newphylogeographical clades, J. Mol. Evol. 53 (2001) 680–689.
V. Valcheva, I. Mokrousov, S. Panaiotov, E. Bachiiska, T. Zozio,C. Sola, et al, Bulgarian specificity and controversialphylogeography of Mycobacterium tuberculosis
spoligotype ST125_BGR, FEMS Immunol. Med. Microbiol. 59 (2010) 90–99.
I. Mokrousov, A. Vyazovaya, O. Narvskaya, Mycobacteriumtuberculosis Latin-American Mediterranean family and itssublineages in the light of robust evolutionary markers, J.Bacteriol. 196 (2014) 1833–1841.
I. Mokrousov, W.W. Jiao, K. Wan, A. Shen, Stranger in astrange land: Ibero-American strain of Mycobacteriumtuberculosis in Tibet, China, Infect. Genet. Evol. 26 (2014) 323–326.
E. Abadia, J. Zhang, T. dos Vultos, V. Ritacco, K. Kremer, E.Aktas, et al, Resolving lineage assignation on Mycobacteriumtuberculosis clinical isolates classified by spoligotyping with anew high-throughput 3R SNPs based method, Infect. Genet.Evol. 10 (2010) 1066–1074.
S. Homolka, M. Projahn, S. Feuerriegel, T. Ubben, R. Diel, et al,High resolution discrimination of clinical Mycobacteriumtuberculosis complex strains based on single nucleotidepolymorphisms, PLoS One 7 (2012) e39855.
L.C. Lazzarini, R.C. Huard, N.L. Boechat, H.M. Gomes, M.C.Oelemann, N. Kurepina, et al, Discovery of a novelMycobacterium tuberculosis
lineage that is a major cause oftuberculosis in Rio de Janeiro, Brazil, J. Clin. Microbiol. 45 (2007) 3891–3902.
B.J. Marais, T.C. Victor, A.C. Hesseling, M. Barnard, A. Jordaan,W. Brittle, et al, Beijing and Haarlem genotypes areoverrepresented among children with drug-resistanttuberculosis in the Western Cape Province of South Africa, J.Clin. Microbiol. 44 (2006) 3539–3543.
V. Valcheva, I. Mokrousov, O. Narvskaya, N. Rastogi, N.Markova, Utility of new 24-locus variable-number tandemrepeattyping for discriminating Mycobacterium tuberculosis
clinical isolates collected in Bulgaria, J. Clin. Microbiol. 46 (2008) 3005–3011.
J. van Embden, M. Cave, J. Crawford, J. Dale, K. Eisenach, B.Gicquel, et al, Strain identification of Mycobacteriumtuberculosis by DNA fingerprinting: recommendations for astandardized methodology, J. Clin. Microbiol. 31 (1993) 406–409.
J. Kamerbeek, L. Schouls, A. Kolk, M. van Agterveld, D. vanSoolingen, S. Kuijper, et al, Simultaneous detection andstrain differentiation of Mycobacterium tuberculosis
fordiagnosis and epidemiology, J. Clin. Microbiol. 35 (1997) 907–914.
C. Demay, B. Liens, T. Burguie`re, V. Hill, D. Couvin, J. Millet,et al, SITVITWEB – a publicly available internationalmultimarker database for studying Mycobacterium tuberculosis
genetic diversity and molecular epidemiology, Infect. Genet.Evol. 12 (2012) 755–766.
S. Tafaj, J. Zhang, Y. Hauck, C. Pourcel, H. Hafizi, et al, Firstinsight into genetic diversity of the Mycobacterium tuberculosis
complex in Albania obtained by multilocus variable-numbertandem-repeat analysis and spoligotyping reveals thepresence of Beijing multidrug-resistant isolates, J. Clin.Microbiol. 47 (2009) 1581–1584.
V.V. Nikolayevskyy, T.J. Brown, Y.I. Bazhora, A.A. Asmolov,Y.M. Balabanova, F.A. Drobniewski, Molecular epidemiologyand prevalence of mutations conferring rifampicin andisoniazid resistance in Mycobacterium tuberculosis
strains fromthe southern Ukraine, Clin. Microbiol. Infect. 13 (2007) 129–138.
[Figure 1], [Figure 2], [Figure 3]
|This article has been cited by|
||Peculiar features of the Mycobacterium tuberculosis population structure in Albania
| ||Silva Tafaj,Igor Mokrousov,Emanuele Borroni,Alberto Trovato,Perlat Kapisyzi,Donika Bardhi,Hasan Hafizi,Silvana Bala,Anyla Bulo,Silvia Bino,Nalin Rastogi,Daniela Cirillo |
| ||Infection, Genetics and Evolution. 2019; : 104136 |
|[Pubmed] | [DOI]|
||Latin-American-Mediterranean lineage of Mycobacterium tuberculosis: Human traces across pathogenís phylogeography
| ||Igor Mokrousov,Anna Vyazovaya,Tomotada Iwamoto,Yuriy Skiba,Ilva Pole,Svetlana Zhdanova,Kentaro Arikawa,Viacheslav Sinkov,Tatiana Umpeleva,Violeta Valcheva,Maria Alvarez Figueroa,Renate Ranka,Inta Jansone,Oleg Ogarkov,Viacheslav Zhuravlev,Olga Narvskaya |
| ||Molecular Phylogenetics and Evolution. 2016; 99: 133 |
|[Pubmed] | [DOI]|