Mycobacterium tuberculosis belonging to family LAM and sublineage RD Rio: common strains in Southern Brazil for over 10 years

RESEARCH PAPER

Mycobacterium tuberculosis belonging to family LAM and sublineage RD^Rio: common strains in Southern Brazil for over 10 years

Renata Oliveira Soares; Maíra Bidart de Macedo; Andrea von Groll; Pedro Eduardo Almeida da Silva

Universidade Federal do Rio Grande, Rio Grande, RS, Brazil

^{Correspondence}

ABSTRACT

A sublineage of Mycobacterium tuberculosis called RD^Rio was described in 2007. Although only recently described, this strain may have been present previously in the population, and its identification in clinical isolates will elucidate bacterial transmission dynamics and host-pathogen interactions. This study evaluated the clonal diversity of the RD^Rio sublineage in clinical isolates from Rio Grande-RS obtained between 1998 and 2001. Among the 45 samples analyzed by the MIRU-VNTR method, there were six clusters with two samples each and 33 orphan strains with unique pattern. The strains were distributed across several different lineages including LAM (34.04%), X (14.89%), Haarlem (12.77%), UgandaI (10.64%), S (4.26%), NEW-1 (2.13%) and Cameroon (2.13%); 14.89% of the strains matched to multiple lineages. RD^Rio strains were present in 28.9% of the samples and 81.25% of the identified strains belonged to the LAM family. The high clonal diversity observed in this study is a constant feature in this region. The RD^Rio sublineage has been in Rio Grande-RS since 1998. The continued monitoring of RD^Rio in clinical isolates will enhance the understanding of its epidemiological significance.

Key words:Mycobacterium tuberculosis, MIRU-VNTR, RDRio.

Introduction

Tuberculosis (TB) is a global public health problem. There were 9.27 million new TB cases in 2007 with an incidence of 139 cases per 100,000 people (World Health Organization, 2009). Brazil, with 72,000 new cases reported annually, is among the 22 countries that account for 80% of the TB cases worldwide (World Health Organization, 2009). Rio Grande city, located in southern Rio Grande do Sul, has an incidence of approximately 70/100,000 within the municipalities prioritized for TB control (Ministério Da Saúde Do Brasil, 2010).

Several molecular biology tools have become available to support TB management. For instance, Mycobacterium tuberculosis genotyping determines the geographical and spatial distribution of strains, the predominance of different genotypes and potential bacterial adaptations to specific human populations (Brudey et al., 2006; Filliol et al., 2006; Gagneux et al., 2006). Genotyping assays using MIRU-VNTR (Mycobacterial Interspersed Repetitive Units Variable Number Tandem Repeat) amplify independent loci and quantify the repeat regions at each locus (Viedma et al., 2005). Repeat pattern similarities are used to classify new isolates into clades, genotype families and/or strains (Supply et al., 2001, 2003).

Lineage identification with MIRU-VNTR is facilitated by MIRU-VNTRplus database, which contains genotyping data from a reference collection of 186 isolates representing the main M. tuberculosis lineages. This collection includes M. tuberculosis strains of the following lineages: W/Beijing, Cameroon, Delhi/Central Asian, East African-Indian, Ghana, Haarlem, Latin American-Mediterranean (LAM), Turkish, S, Uganda I and II, Ural, and X (Allix-Béguec et al., 2008; Supply et al., 2006).

Recently, the predominant genotype of M. tuberculosis in Rio de Janeiro, Brazil, was identified (Lazzarini et al., 2007). This sublineage, RD^Rio, is a derivative of the LAM lineage and has a deletion of ten genes (26.317 kb), or approximately 0.6% of the genome (Cole 2002; Lazzarini et al., 2007). RD^Rio strains may have higher transmission rates and could cause more severe TB (Lazzarini et al., 2008). In this study, the RD^Rio strains in isolates from Rio Grande between 1998 and 2001 were characterized, and the major lineages present in the population were studied. Despite being only recently described, RD^Rio strains could possibly have been circulating for much longer, and its identification in older clinical isolates may lead to a better understanding of its transmission dynamics and interactions with its host.

Materials and Methods

Samples

Forty-five M. tuberculosis strains were isolated from patients admitted to the University Hospital of Rio Grande, RS, Brazil, between 1998 and 2001. These strains are archived in the Mycobacteriology Laboratory, Faculty of Medicine at Federal University of Rio Grande.

Genomic DNA extraction

Each bacterial pellet was transferred to a microtube and resuspended in 300 μL of Tris-EDTA (TE). The bacteria were inactivated by heating at 80°C for 30 min. The microtube was centrifuged at 5000 g for 5 min, and the supernatant containing the DNA was removed. The extraction products were stored at -20°C.

MIRU-VNTR genotyping

The MIRU-VNTR 12 loci method was performed according to Supply et al. (2000). The results from PCR of each locus (2, 4, 40,10, 16, 20, 23, 24, 26, 27, 31 and 40) were combined to form a 12-digit allelic profile that was used to analyze the genetic relationships.

Strain identification was performed with tools on the MIRU-VNTRplus website (http://www.miru-vntrplus.org). The distance between each locus was compared with the data from 186 strains in the website reference collection. While this distance is adjustable, we initially selected 0.17, which corresponds to four loci of difference tolerance when MIRU-VNTR is the only method used (Allix-Béguec et al., 2008). Similar strains were grouped, and based on the database reference sequences, strains in the samples of interest were identified.

Analysis of genetic relationships

Dendrograms were constructed with the tools on the MIRU-VNTRplus website (Allix-Béguec et al., 2008), and clusters were grouped by UPGMA (Unweighted Pair Group Method with Arithmetic Mean). Clusters were defined with at least two strains of M. tuberculosis isolates with identical patterns in different patients. The discriminatory power of MIRU-VNTR was calculated using the discriminatory index of Hunter-Gaston (HGDI). The HGDI was calculated by the following equation: HGDI = 1 - [1/N(N-1) ∑ n_j (n_j-1)] with j = 1, where N is the total number of strains, is the total number of different patterns and n_j is the number of strains belonging to each pattern (Hunter and Gaston, 1988). The allelic diversity (h) of MIRU-VNTR in each of the 12 loci was calculated by the equation h = 1 - ∑ xi² [n/(n - 1)], where xi is the allele frequency at each locus and n is the number isolates (Selander et al., 2006).

Identification of the RD^Rio lineage strains

The identification of RDRio lineage strains was performed with a PCR protocol modified as previously described (Gibson et al., 2008). PCRs were performed in 30 μL containing 1.5 mM MgCl₂, 10 mM Tris-HCl (pH 8.3), 400 mM of each primer (RDRioBrg and IS1561), 200 mM deoxyribonucleoside triphosphates (Invitrogen^TM), 6% glycerol, 1 U Taq DNA polymerase (Invitrogen^TM) and 1.5 μL of isolate DNA.

DNA amplifications were performed under the following conditions: 1 cycle at 95 °C for 5 min; 35 cycles at 95 °C for 1 min, 60 °C for 1 min, and 72 °C for 2 min; and 1 cycle at 72 °C for 10 min. PCR products were electrophoresed on a 1.5% agarose gel for 1 hour at 90 V, and the presence of either a 1175 bp band, corresponding to an RDRio strain, or a 530 bp band, corresponding to non-RDRio strains, was observed.

Results and Discussion

MIRU-VNTR

Thirty-nine different MIRU-VNTR patterns were detected in the 45 strains analyzed. Twelve strains were grouped into six clusters of two samples each, and 33 strains had unique patterns (Figure 1).

The analysis of allelic diversity at each locus (Table 1) showed that loci 10, 16 and 40 were highly discriminatory (h ≥ 0.6), while loci 23, 26, 27 and 31 were moderately discriminatory (h ≥ 0.3). Loci 2, 4, 20, 24 and 39 had low discriminatory power (h < 0.3). These results are similar to those in multiple previous studies (1, 9, 17, 19, 21). As noted by Sola et al. (2003), there is a hierarchy of polymorphism, such that loci 10, 23, 26, 31 and 40 have greater discriminatory power than the others.

The MIRU-VNTR 12 loci method is quick and easy, and it has been evaluated in several studies in different locations (Kang et al., 2009). Because it has a discriminatory power equivalent to RFLP (Restriction Fragment Length Polymorphism), this method is an alternative to IS6110-RFLP for epidemiological studies (Mazars et al., 2001; Supply et al., 2000). As in other studies, we found that the HGDI of MIRU-VNTR was 0.994, indicating a high discriminatory power (Sun et al., 2004; Silva et al., 2009; Von Groll et al., 2010).

This technique has limitations, and it could be combined with an additional genotyping method for greater accuracy (Cowan et al., 2005; Van Deutekom et al., 2005; Kang et al., 2009). Different combinations of MIRU and other loci have been evaluated. In future studies, MIRU-VNTR could be standardized to include 15 loci for epidemiological studies and 24 loci for phylogenetic studies of M. tuberculosis (Supply et al., 2006).

The strains were distributed among the following lineages: LAM (16 strains, 34.04% of the total); X (7, 14.89%); Haarlem (6, 12.77%); Uganda I (5, 10.64%); S (2, 4.26%); NEW-1 (1, 2.13%); and Cameroon (1, 2.13%). Seven strains (14.89%) were matches to multiple lineages (i.e., more than one family was within the same phylogenetic distance).

Identification of RD^Rio strains

Of the 45 strains analyzed, 13 (28.9%) belonged to the RD^Rio lineage, and the other 32 strains were grouped as non-RD^Rio. In previous studies, the RD^Rio lineage comprised 30% of all isolates in Rio de Janeiro (Lazzarini et al., 2007), 37% in Belo Horizonte (Lazzarini et al., 2008) and 38% in Rio Grande (Von Groll et al., 2010). Because this lineage was also detected in other parts of the world (Gibson et al., 2008), it is not restricted to Brazil. Of the six clusters grouped by MIRU-VNTR analysis, four were formed by RD^Rio strains, and the other two clusters were composed of non-RD^Rio strains.

Similar to other studies, all RD^Rio lineage strains identified here belonged to the LAM family, which is responsible for 15% of TB cases worldwide (Brudey et al., 2006; Lazzarini et al., 2007; Lazzarini et al., 2008; Von Groll et al., 2010). We found that 81.25% of the LAM strains found in this study were from the RDRio sublineage, but more studies are needed to evaluate any possible selective advantages of LAM RD^Rio strains.

Conclusion

High genetic diversity within circulating M. tuberculosis was observed in this study and others; this feature may contribute to high TB prevalence (Silva et al., 2009; Von Groll et al., 2010). However, the absence of clonal dispersion may also indicate that effective medical treatment would result in reduced transmission. The RD^Rio sublineage has been present in Rio Grande-RS since 1998. The continued monitoring of RD^Rio in clinical isolates will elucidate its epidemiological significance.

Submitted: July 29, 2011

Approved: April 04, 2013

All the content of the journal, except where otherwise noted, is licensed under a Creative Commons License CC BY-NC.

• Ali A, Hasan Z, Tanveer M, Siddiqui AR, Ghebremichael S, Kallenius G, Hasan R (2007) Characterization of Mycobacterium tuberculosis central asian strain I using mycobacterial interspersed repetitive unit genotyping. BMC Microbiol 7:76.
• Allix-Béguec C, Harmsen D, Weniger T, Supply P, Niemann S (2008) Evaluation and user-strategy of MIRU-VNTRplus, a multifunctional database for online analysis of genotyping data and phylogenetic identification of Mycobacterium tuberculosis complex isolates. J Clin Microbiol 46:2692-2699.
• Brudey K, Driscoll JR, Rigouts L, Prodinger WM, Gori A, Al-Hajoj SA, Allix C, Aristimuno L, Arora J, Baumanis V, Binder L, Cafrune P, Cataldi A, Cheong S, Diel R, Ellermeier C, Evans JT, Fauville-Dufaux M, Ferdinand S, Garcia De Viedma D, Garzelli C, Gazzola L, Gomes HM, Guttierez MC, Hawkey PM, Van Helden PD, Kadival GV, Kreiswirth BN, Kremer K, Kubin M, Kulkarni SP, Liens B, Lillebaek T, Ho ML, Martin C, Martin C, Mokrousov I, Narvskaia O, Ngeow YF, Naumann L, Niemann S, Parwati I, Rahim Z, Rasolofo-Razanamparany V, Rasolonavalona T, Rossetti ML, Rusch-Gerdes S, Sajduda A, Samper S, Shemyakin IG, Singh UB, Somoskovi A, Skuce RA, Van Soolingen D, Streicher EM, Suffys PN, Tortoli E, Tracevska T, Vincent V, Victor TC, Warren RM, Yap SF, Zaman K, Portaels F, Rastogi N, Sola C (2006) Mycobacterium tuberculosis complex genetic diversity: mining the fourth international spoligotyping database (SpolDB4) for classification, population genetics and epidemiology. BMC Microbiol 6:23.
• Cole ST (2002) Comparative and functional genomics of the Mycobacterium tuberculosis complex. Microbiol 148:2919-2928.
• Cowan LS, Diem L, Monson T, Wand P, Temporado D, Oemig TV, Crawford JT (2005) Evaluation of a two-step approach for large-scale, prospective genotyping of Mycobacterium tuberculosis isolates in the United States. J Clin Microbiol 43:688-695.
• Filliol I, Motiwala AS, Cavatore M, Qi W, Hazbo MH, Bobadilla Del Valle M, Fyfe J, García-García L, Rastogi N, Sola C, Zozio T, Guerrero MI, León CI, Crabtree J, Angiuoli S, Eisenach KD, Durmaz R, Joloba ML, Rendón A, Sifuentes-Osornio J, Ponce De León A, Cave MD, Fleischmann R, Whittam TS, Alland D (2006) Global phylogeny of Mycobacterium tuberculosis based on single nucleotide polymorphism (SNP) analysis: insights into tuberculosis evolution, phylogenetic accuracy of other DNA fingerprinting systems, and recommendations for a minimal standard SNP set. J Bacteriol 188:759-772.
• Gagneux S, Deriemer K, Van T, Kato-Maeda M, De Jong BC, Narayanan S, Nicol M, Niemann S, Kremer K, Gutierrez MC, Hilty M, Hopewell PC, Small PM (2006) Variable host-pathogen compatibility in Mycobacterium tuberculosis. Proc Natl Acad Sci U S A 103:2869-2873.
• Gibson AL, Huard RC, Gey Van Pittius NC, Lazzarini LC, Driscoll J, Kurepina N, Zozio T, Sola C, Spindola SM, Kritski AL, Fitzgerald D, Kremer K, Mardassi H, Chitale P, Brinkworth J, Garcia De Viedma D, Gicquel B, Pape JW, Van Soolingen D, Kreiswirth BN, Warren RM, Van Helden PD, Rastogi N, Suffys PN, Lapa E Silva J, Ho JL (2008) Application of sensitive and specific molecular methods to uncover global dissemination of the major RDRio Sublineage of the Latin American-Mediterranean Mycobacterium tuberculosis spoligotype family. J Clin Microbiol 46:1259-1267.
• Han H, Wang F, Xiao Y, Ren Y, Chao Y, Guo A, Ye L (2007) Utility of mycobacterial interspersed repetitive unit typing for differentiating Mycobacterium tuberculosis isolates in Wuhan, China. J Med Microbiol 56:1219-1223.
• Hunter PR, Gaston MA (1988) Numerical index of discriminatory ability of typing system: an application of Simpson's index of diversity. J Clin Microbiol 26:2465-2466.
• Kang H, Ryoo S, Park Y, Lew W (2009) Evaluation of the selected 12-locus MIRU for genotyping beijing family Mycobacterium Tuberculosis in Korea. Tuberc Respir Dis 67:499-505.
• Lazzarini LCO, Huard RC, Boechat NL, Gomes HM, Oelemann MC, Kurepina N, Shashkina E, Mello FC, Gibson AL, Virginio MJ, Marsico AG, Butler WR, Kreiswirth BN, Suffys PN, Lapa E Silva JR, Ho JL (2007) Discovery of a novel Mycobacterium Tuberculosis lineage that Is a major cause of tuberculosis in Rio de Janeiro, Brazil. J Clin Microbiol 45:3891 - 3902.
• Lazzarini LCO, Spindola SM, Bang H, Gibson AL, Weisenberg S, Da Silva Carvalho W, Augusto CJ, Huard RC, Kritski AL, Ho JL (2008) RDRio Mycobacterium tuberculosis infection is associated with a higher frequency of Cavitary pulmonary disease. J Clin Microbiol 46
• Mazars E, Lesjean S, Banuls AL, Gilbert M, Vincent V, Gicquel B, Tibayrenc M, Locht C, Supply P (2001) High-resolution minisatellite based typing as a portable approach to global analysis of Mycobacterium tuberculosis molecular epidemiology. Proc Natl Acad Sci USA 98:1901-1906.
• Ministério Da Saúde Do Brasil (2010) Sistema de informação de agravos de notificação http://dtr2004.saude.gov.br/sinanweb/
• Selander RK, Caugant DA, Ochman H, Musser JM, Gilmour MN, Whittam TS (1986) Methods of multilocus enzyme electrophoresis for bacterial population genetics and systematics. Appl Environ Microbiol 51:873-884.
• Silva ABSS, Von Groll A, Félix C, Conceição FR, Spies FS, Scaini CJ, Rossetti ML, Borsuk S, Dellagostin OA, Almeida Da Silva PE (2009) Clonal diversity of Mycobacterium tuberculosis isolated in a sea port city in Brazil. Tuberculosis 89:433-437.
• Sola C, Filliol I, Legrand E, Lesjean S, Locht C, Supply P, Rastogi N (2003) Genotyping of the Mycobacterium tuberculosis complex using MIRUs: association whit VNTR and spolygotyping for molecular epidemiology and evolutionary genetics. Infect Genet Evo 3:125-133.
• Sun Y, Bellamy R, Lee ASG, Tang S, Ravindran S Et (2004) Use of mycobacterial interspersed repetitive unit-variable-number tandem repeat typing to examine genetic diversity of Mycobacterium tuberculosis in Singapore. J Clin Microbiol 42:1986-1993.
• Supply P, Allix C, Lesjean S, Cardoso-Oelemann M, Rusch-Gerdes S, Willery E, Savine E, De Haas P, Van Deutekom H, Roring S, Bifani P, Kurepina N, Kreiswirth B, Sola C, Rastogi N, Vatin V, Gutierrez MC, Fauville M, Niemann S, Skuce R, Kremer K, Locht C, Van Soolingen D (2006) Proposal for standardization of optimized mycobacterial interspersed repetitive unit-variable-number tandem repeat typing of Mycobacterium tuberculosis. J Clin Microbiol 44:4498-4510.
• Supply P, Lesjean S, Savine E, Kremer K, Van Soolingen D, Locht C (2001) Automated high-throughput genotyping for study of global epidemiology of Mycobacterium tuberculosis based on mycobacterial interspersed repetitive units. J Clin Microbiol 39:3563-3571.
• Supply P, Mazars E, Lesjean S, Vincent V, Gicquel B, Locht C (2000) Variable human minisatellite-like regions in the Mycobacterium tuberculosis genome. Mol Microbiol 36:762-771.
• Supply P, Warren RM, Bañuls AL, Lesjean S, Van Der Spuy GD, Lewis LA, Tibayrenc M, Van Helden PD, Locht C (2003) Linkage disequilibrium between minisatellite loci supports clonal evolution of Micobacterium tuberculosis in a high tuberculosis incidence area. Mol Microbiol 47:529-538.
• Van Deutekom H, Supply P, Haas PE, Willery E, Hoijng SP, Locht C, Coutinho RA, Van Soolingen D (2005) Molecular typing of Mycobacterium tuberculosis by mycobacterial interspersed repetitive unit-variable-number tandem repeat analysis, a more accurate method for identifying epidemiological links between patients with tuberculosis. J Clin Microbiol 43:4473-4479.
• Viedma DG, Rodriguez NA, Andres S, Serrano MJR, Bouza E (2005) Characterization of clonal complexity in tuberculosis by mycobacterial interspersed repetitive unit-variable-number tandem repeat typing. J Clin Microbiol 43:5660-5664.
• Von Groll A, Martin A, Felix C, Prata PF, Honscha G, Portaels F, Vandame P, Da Silva PE, Palomino JC (2010) Fitness study of the RDRio lineage and Latin American-Mediterranean family of Mycobacterium tuberculosis in the city of Rio Grande, Brazil. FEMS Immunol Med Microbiol 58:119-127.
• World Health Organization (2009) Global Tuberculosis Control 2009 - Surveillance, Planning, Financing www://who.int/tb/publications/global_report/es

Publication Dates

• Publication in this collection
  27 Mar 2014
• Date of issue
  Dec 2013

History

• Received
  29 July 2011
• Accepted
  04 Apr 2013