Use of cluster-graphs from spoligotyping data to study genotype similarities and a comparison of three indices to quantify recent tuberculosis transmission among culture positive cases in French Guiana during a eight year period
© Guernier et al; licensee BioMed Central Ltd. 2008
Received: 06 November 2007
Accepted: 14 April 2008
Published: 14 April 2008
French Guiana has the highest tuberculosis (TB) burden among all French departments, with a strong increase in the TB incidence over the last few years. It is now uncertain how best to explain this incidence. The objective of this study was to compare three different methods evaluating the extent of recent TB transmission in French Guiana.
We conducted a population-based molecular epidemiology study of tuberculosis in French Guiana based on culture-positive TB strains (1996 to 2003, n = 344) to define molecular relatedness between isolates, i.e. potential transmission events. Phylogenetic relationships were inferred by comparing two methods: a "cluster-graph" method based on spoligotyping results, and a minimum spanning tree method based on both spoligotyping and variable number of tandem DNA repeats (VNTR). Furthermore, three indices attempting to reflect the extent of recent TB transmission (RTIn, RTIn-1 and TMI) were compared.
Molecular analyses showed a total amount of 120 different spoligotyping patterns and 273 clinical isolates (79.4%) that were grouped in 49 clusters. The comparison of spoligotypes from French Guiana with an international spoligotype database (SpolDB4) showed that the majority of isolates belonged to major clades of M. tuberculosis (Haarlem, 22.6%; Latin American-Mediterranean, 23.3%; and T, 32.6%). Indices designed to quantify transmission of tuberculosis gave the following values: RTIn = 0.794, RTIn-1 = 0.651, and TMI = 0.146.
Our data showed a high number of Mycobacterium tuberculosis clusters, suggesting a high level of recent TB transmission, nonetheless an estimation of transmission rate taking into account cluster size and mutation rate of genetic markers showed a low ongoing transmission rate (14.6%). Our results indicate an endemic mode of TB transmission in French Guiana, with both resurgence of old spatially restricted genotypes, and a significant importation of new TB genotypes by migration of TB infected persons from neighgouring high-incidence countries.
The resurgence of human tuberculosis (TB) in many industrialized countries in the late 1980s and early 1990s  has focused critical attention on a disease that was thought to be under control, and as such, has been neglected for too long. Effective prevention and control of Mycobacterium tuberculosis, the causative agent of TB, must be based on a clear understanding of how the disease is transmitted, how infection becomes established, and how infection progresses to clinical disease. This has motivated investigations on TB transmission among persons at higher risk, such as the homeless [2, 3] or HIV-infected persons [4, 5]. Furthermore, the improvement of molecular techniques to characterize Mycobacterium strains, especially M. tuberculosis, and the concomitant accumulation of epidemiological data over large spatial scales has increased the possibility to trace transmission routes of Mycobacterium outbreaks (e.g. M. leprae in ). The identification of M. tuberculosis strains can be used to identify patients which have some apparently identical genotypes, thus potentially associated to a same source of contamination  and designated as «clustered cases» [8, 9]. The existence of genetic clusters is thought to be indicative of recent TB transmission [10, 11], in opposition to cases arising either from imports or from reactivation of latent infection. Featuring epidemiological links between patients allow inferences about the level of ongoing TB transmission, thus the level of epidemicity/endemicity of the infection in a defined geographical area, and underline some potential risk factors associated to the disease . For example, when an unusual number of TB cases occurs over time, it is possible to do early molecular analyses to determine whether this cluster of cases represents temporal coincidence, or a genetic cluster due to a local chain of transmission. If fingerprinting demonstrates different strains, the cases are not due to transmission, and there is no need for further epidemiological evaluation.
The underlined hypothesis of the use of molecular information for the evaluation of ongoing transmission is that mutations are slower than transmission rate , which is not true for all genetic markers. A large number of molecular methods are available for the characterization of infra-specific M. tuberculosis strains. Nonetheless the choice of a method strongly affects the results, as low mutation rates tend to less discriminate the genetic differences, whereas high mutation rates tend to increase the number of clusters considered . This article investigates the genotypic diversity of M. tuberculosis strains isolated in French Guiana over the 1996–2003 period, characterized by two molecular methods, and attempts to determine potential transmission between infected persons. We used a recent method proposed by Tanaka & Francis  to visualise molecular epidemiological data, i.e. cluster-graphs, which also provides a framework for evaluating ongoing transmission. We also compared the method used to previous indices attempting to reflect the extent of recent TB transmission.
French Guiana is a French overseas department of 32,432 sq mi. It is located on the north-east of South America, between Suriname and Brazil. The region has the highest TB rates among all French departments, with an average annual incidence of 63.1 TB cases per 100,000 population estimated during the period 1996–2003 . Among all TB cases identified from January 1, 1996 to December 31, 2003, 345 were confirmed by isolation and cell culture of M. tuberculosis by the Mycobacterium Laboratory of the Institut Pasteur de la Guyane (IPG). A total of 344 clinical isolates of M. tuberculosis were isolated from patients with culture-positive TB, representing 342 patients (two patients only were contaminated twice, by two different strains). The three missing cases (0,86%) correspond to contaminated cultures. After a prior explanation of the recommended investigations for the diagnostic purposes by a physician, patients were asked to provide oral consent and the pathological specimens were obtained from suspected tuberculosis patients residing in French Guiana. The cultures were performed using Löwenstein-Jensen slants at 37°C at the Pasteur Institute of French Guiana. Positive cultures (where no identification was attached) were sent to the Pasteur Institute of Guadeloupe which has been serving as a regional reference laboratory for tuberculosis and mycobacterial for the Caribbean since 1994. The cultures were duly identified as M. tuberculosis complex using classical biochemical tests and the AccuProbe test (GenProbe Inc., San Diego, CA) and subjected to drug-susceptibility testing using the proportional method. The bacteriological results were communicated to the physicians and/or hospital services responsible for patient care and treatment, through the Pasteur Institute of French Guiana. DNA-based molecular typing was also performed (see below).
Molecular typing was performed using spoligotyping and VNTRs analyses. It was previously suggested that a combined spoligotyping-5 VNTR loci genotyping scheme may have a discriminatory power close to that of IS6110 RFLP, the historical "golden standard" in tuberculosis genotyping . Hence, for ease and resources reasons, we choosed this scheme, which was less expensive in time and money than the combined spoligotyping-12 VNTR loci scheme, i.e. MIRUs . Recent improvement in VNTR typing now recommended a new "universal" 15 or even 24 VNTR loci genotyping scheme for M. tuberculosis .
The bacterial DNA was prepared by the classical cetyl-trimethyl ammonium bromide method (CTAB)  and used for spoligotyping (spacer oligonucleotide typing) and Variable Number of DNA Tandem Repeats (VNTRs) analyses. Spoligotyping was performed using a home-made membrane with 43 covalently bound oligonucleotides and PCR was achieved using primers designated DRa and DRb, with DRa biotinylated 5' to amplify the whole DR region as described previously . As spoligotyping used alone is known to overestimate the number of epidemiological links (depending on the settings, but by around 30%), it was suggested to be used in association with another rapid fingerprinting technique . Hence, the isolates clustered by spoligotyping (and only those first-clustered isolates) were further subtyped using VNTRs as described previously  in order to provide a second independent indicator of clonality. The number of copies for each exact tandem repeat (ETR) was documented as a five-digit numbers representing allele profiles ETR-A to ETR-E .
Spoligotype patterns are designated as 43-character-long strings consisting of black and white squares, indicating respectively the presence or the absence of an individual spacer. The totality of the 344 isolates was entered in this binary format as Excel (Microsoft, Cupertino, CA) spreadsheets. The spoligotype designations (shared types or ST, defined as an identical spoligotype found in ≥ 2 individual patient isolates) were attributed by comparing the patterns obtained with those included in an international spoligotyping database held at the Pasteur Institute of Guadeloupe, designated as spolDB4 . At the time of the matching analysis, the updated spolDB4 version contained 31,642 patterns distributed into 2,393 shared types in 114 countries (an online version of this database is now available at http://www.pasteur-guadeloupe.fr:8081/SITVITDemo/). Patterns referenced only once are designated as true orphans (in opposition to pseudo-orphans, i.e. isolates found to be unique in this study but for whom a counterpart exists in the database) and are not labelled with an ST number. Major phylogenetic clades were assigned to STs according to signatures provided in SpolDB4 (the reader is referred to the original paper  for a detailed description). These included specific signatures for various M. tuberculosis sub-species (M. bovis, M. microti, M. caprae, M. pinipedii, M. africanum), as well as rules defining major lineages/sub-lineages for M. tuberculosis stricto sensu. The latter included the Central-Asian (CAS) clade and 2 sub-lineages; the East-African-Indian (EAI) clade and 9 sub-lineages; the Haarlem (H) clade and 3 sub-lineages; the Latin-American-Mediterranean (LAM) clade and 12 presumed sub-lineages; the "Manu" family and 3 sub-lineages; the S clade; the IS6110-low banding X clade and 3 sub-lineages; and an ill-defined T clade with 5 sub-lineages. The T clade corresponds to Principal Genetic Groups (PGG) 2 and 3, and remains defined by default (no specific signature sequences). M. tuberculosis families are phylogeographically-specific, as generally indicated by their names.
Phylogeography and cluster-graphs
Data were analysed by grouping isolates into clusters of identical genotypes and by organizing all genotypes present in the overall sample according to genetic relationships, i.e. by constructing a phylogeny. We used and compared two methods for inferring potential phylogenetic relationships: (1) a classification based on spoligotyping results, which we call a "cluster-graph" and which is based on parsimony principles, and (2) a minimum spanning tree (MST) method which is based on a spoligotyping and VNTR combined distance matrix, computed using the BioNumerics software (Applied Maths, Sint Marten-Latems, Belgium).
The cluster-graph method could not be applied to the combined spoligotyping-VNTR data, as the manual inference of tree edges between genotypes needs to be supported by the number of mutations events between different genotypes. The VNTR analyses giving occurrence of ETR-A to E do not give any information about the genetic distance between different VNTR patterns. VNTR analyses were included on the construction of the minimum spanning tree method, in order to provide further discrimination between strains for isolates with identical spoligotyping results, Distance matrices are calculated and the obtained groupings can be used to build a dendrogram.
Estimation of transmission
where is an independent estimate of the mutation rate of the genetic marker, and v 1 is the number of single-step mutation events inferred from the data . We used a mutation rate for spoligotypes of 0.039 events per year computed from previous estimates in the literature [14, 29]. All three indices are expressed as a ratio from 0 to 1, which can also be expressed in percentage. It must be interpreted in terms of "percent of isolates in cluster", which is indicative of the rate of recent TB transmission. Thus, a result of 0.5 indicates that 50% of the isolates of a data set were probably related to a recent transmission event.
Molecular typing results
The spoligotyping was performed as a first-line screening method on 344 M. tuberculosis clinical isolates from patients living in French Guiana, followed by VNTR performed on 273 clustered isolates. Among the 344 typed isolates, spoligotyping generated 49 clustered patterns totalling 273 isolates (with 2 to 48 patients per cluster), 25 true orphan patterns, and 46 pseudo-orphans, i.e. patterns already found in SpolDB4 but present as single in this study, with a total amount of 120 different patterns. Analysis of the frequency of major spoligotypes with SpolDB4 allowed a differentiation between ubiquitous types present in all the continents and types more specific to French Guiana. Several types were endemic in French Guiana (i.e. ST 66, 76, 94, 1084, 1340, 1526 and 25 orphan isolates), representing 25.8% of the shared-types, and 12.2% of the total isolates. Ubiquitous types represent 40.8% of the 120 individual patterns, i.e. 70.6% of the isolates.
Number and frequencies of isolates clustered by spoligotyping alone and further subtyped by VNTRs.
Isolates clustered by spoligotyping
Isolates clustered by spoligotyping and VNTR
Total number of isolates studied
Number of clusters
Mean number of isolates per cluster
Number of clustered isolates (and %)
Number of unclustered isolates (and %)
Detailed results of clustering for the five largest spoligotype clusters with combined spoligotyping – 5 VNTR loci genotyping.
Number of different VNTR
Number of isolates/cluster
6, 5, 5, 4, 4, 3, 2, 2, 2, 2
24, 7, 4, 2, 2
Estimation of transmission
Applying the indices previously described to our French Guiana data set gives the following values: RTIn = 0.794, RTIn-1 = 0.651, and TMI = 0.146, with n = 344 isolates, n c = 273 isolates in cluster, c = 49 clusters, g = 120 spoligotypes and v 1 = 82 edges involving a 1-step change. The number v 1 is equal to the number of genotypes g in the data set minus the number of connected components, which is calculated using the cluster-graph (Figure 3). The RTI values here are very high compared to the same statistic computed on spoligotypes data (TMI), which means a lower recent transmission rate evaluated with TMI method. This difference is further explained in the discussion.
Molecular typing results
Our genotypic results seems to show a low level of TB transmission as evidence through: (1) a high diversity of spoligotypes (120 individual patterns for 344 isolates), (2) a low mean cluster size (Table 1; Figure 2A) in spite of the high level of clustering (79.4% of isolates included in clusters; Table 1), and (3) a high degree of unique STs (59.2%), i.e. ST that cannot be related to a recent event of transmission. Moreover, the biggest clusters (48 and 39 isolates respectively; see Figure 2A) concern ubiquitous spoligotypes 50 and 53, widely spread on all the five continents (spreading index SI equal to 20.95 and 28.09 respectively ), and for which interpretation in terms of recent transmission remains difficult. However, we should keep in mind the suboptimal choice of the spoligotyping-VNTR genotyping scheme, which provided only a slight or no improvement compared to spoligotyping only. However, for practical reasons, and given the time frame covered (1996–2003), a full MIRU12 retrospective characterization of our sample could not be performed.
Cluster-graph and minimum spanning tree
Phylogeography from spoligotyping
Our results are coherent with migration and colonization history of human populations in French Guiana. The Haarlem clade (23% of our strains) is considered to be a clade of European descent [33–35] and related spoligotypes might have been introduced during previous colonization of the three Guiana (i.e. French Guiana, Suriname and Guyana) by migrants from France, Holland and Great Britain, respectively. Occurrence of 23% of spoligotypes belonging to the LAM clade (Latin-American and Mediterranean) might be best explained by population flows with neighbouring Latin-American countries, in particular Brazil, or by historical Lusitanian or Spanish settling. Concerning the EAI clade (East-African Indian; 12% of STs), the absence of many likely-ancestral STs was noticeable (e.g. ST 236, the putative EAI ancestor, and 7 out of 10 absent STs on the one-mutation circle illustrated in Figure 3), whereas distant genotypes were censused, presumably indicating a more recent expansion. The genotypes present in French Guiana might have been introduced by successive arrivals of Hakka settlers from Taiwan, Singapour, Hong-Kong, Vietnam and continental China since 1820 http://www.tlfq.ulaval.ca/axl/amsudant/guyanefr2.htm, and more recently, by settlers from the Lao-Hmong community which arrived in French Guyana in 1977 . But there are obvious limitations to what can be achieved by spoligotyping, and the dominance of T strains in the current study (33%) may (at best) demonstrate one of these. The T clade indeed remains an ill-defined ubiquist family of M. tuberculosis found all over the world [33, 34], thus providing few phylogeographical information.
Recent TB transmission in French Guiana
The calculated TMI indicates a recent transmission rate of 14.6% on the 1996–2003 period, instead of 79.4% and 65.1% of recent TB transmission calculated with RTIn and RTIn-1 methods respectively. In a previous study , the rate of recent transmission was calculated to be 49.3% in French Guiana, using the RTIn-1 method based on the combined spoligotyping and VNTRs genotyping. The difference between our RTIn-1 result and the previous study is probably due to a different sampling and a smaller time frame of the present study (396 strains from 1994 to 2003 in  as compared to 344 isolates from 1996 to 2003 in the present study). Comparing RTI and TMI indices, we obtained a lower ongoing transmission rate with TMI method. The three methods are based on the same hypothesis about the potential epidemiological link between patients sharing the same M. tuberculosis strains. But the major weakness with RTI methods is their failure to account for strains diversity and mutation. The relationship between sample size n and the number of observed genotypes g is complex and depends on sampling, mutation, and population history, which were not integrated in previous RTI indices. In addition to being useful in visualizing data, cluster-graphs provide information about the abundance of particular genotypes in the sample, as well as partial information about possible evolutionary relationships, and can then be used as a way of understanding indices measuring the severity of outbreaks. The TMI result is the only one to be coherent with genotyping results described above (high genotypic diversity, small clusters, many orphans) which supposed a low level of recent TB transmission. The concomitant hypothesis is an endemic mode of TB transmission in French Guiana, with little local contamination occurring from person to person, but an overall persistence of the disease due to the large number of imported cases from countries with high TB incidence rates.
New "transmission mutation index" method
There are at least two biases that can occur when calculating TMI, and that may conduct to strictly opposite conclusions. The first possible bias is the emphasis of a single epidemic picture, whereas there are several parallel epidemics. For example, when considering two genotypes g1 and g2 differing by one single mutation step, the two isolates will be associated in the cluster-graph and thus related to a same epidemic, unless they potentially acquired the related genotypes independently by convergence, due to homoplasia. The infinite allele model, which underlies the TMI method, does not take into account such possibility, which can definitely occur.
The second possible bias is exactly the opposite, i.e. the emphasis of several epidemics, whereas there is only one. A sample is, by nature, incomplete. Thus, when isolates included in a cluster-graph are sampled from a population, some cases might be missing, which will induce a possible under-estimation of a cluster size, and/or over-estimation of the number of different clusters. If one intermediate genotype is missing between two isolates, we will detect two epidemics unless there is a single transmission event.
A third bias can also be discussed here, considering the hypothesis we made that, when person A contaminates person B, their TB strains may differ at maximum by one mutation event. In fact, this hypothesis is likely to be true for spoligotypes, which have a low mutation rate. So, the possible bias due to points 1 and 2 are believed to be greater than this last one.
We must also notice that, because genotyping requires the availability of a viable isolate of M. tuberculosis, the population for which genotyping is performed must be a subset of all cases with positive cultures. In this study, 342 of 345 confirmed cases were spoligotyped. In turn, cases with positive cultures must be a subset of all cases (culture positive, negative, or not done). A previous study on the Ile-de-Cayenne, French Guiana, showed that 61.7% of identified cases were confirmed by isolation and cellule culture of M. tuberculosis . This must be taken into account when interpreting our results.
Last but not least, one should also pinpoint that no true analysis of patient demographical and medical data could be performed within the frame of this study. This is no doubt a strong limitation, since it is well known that molecular epidemiology performs well when collaboration between conventional and molecular epidemiology is achieved. The main reason relates to the difficulty to provide adequate human resources that would work for such a long period on the epidemiology of tuberculosis; however, even if a careful case-by-case analysis of the patient file data could not be assessed within the frame of the study, we fully assume that our results are a mirror of the epidemic situation of tuberculosis in French Guiana.
Our results confirmed that the use of cluster-graphs from spoligotyping data to derive a TMI index provides more sensitive results than previous RTI indices in order to measure the severity of outbreaks. It allowed us to show an endemic mode of TB transmission in French Guiana associated to low local transmission of the disease, with both resurgence of old spatially restricted genotypes (endemic to French Guiana), and an importation of new TB genotypes by migration of TB infected people from high-incidence countries. Endemic genotypes might also be due, to a lower extent, to local mutation events. Those results highlight the need for specific strategies of TB control in this region, considering those new epidemiological hypotheses.
This study also suggested that our methods choice in order to get a clear picture of the TB snapshot in a high prevalence setting such as French Guiana was indeed suboptimal. Genotyping of pathogens is a fast moving field that also requires long-term experience and stability of methods to provide golden standards in molecular characterization. A decade of IS6110-RFLP predominance, has now been efficiently challenged by the power of the Multi Locus VNTR analysis (MLVA) approach, which has proven to be useful for many pathogens. Meanwhile, the spoligotyping approach, once thought to be not discriminative enough, has revealed to be an excellent long-term choice in TB genotyping, not only because it is cheap and robust, but also because the underlying role of the DR (a member of the CRISPR family of repeats) is now better understood , and because updated databases are regularly maintained. An ultimate player in this game, i.e. SNP genotyping, could become the most useful tool during the next decade [38, 39]. However, for the time-being, there is no emerging consensus on what is the best genotyping scheme to understand TB transmission based on molecular characterization of clinical isolates.
Availability and requirements
At the time of the matching analysis, the updated spolDB4 version contained 31,642 patterns distributed into 2,393 shared types in 114 countries (an online version of this database is now available at http://www.pasteur-guadeloupe.fr:8081/SITVITDemo/).
The genotypes present in French Guiana might have been introduced by successive arrivals of Hakka settlers from Taiwan, Singapour, Hong-Kong, Vietnam and continental China since 1820 http://www.tlfq.ulaval.ca/axl/amsudant/guyanefr2.htm, and more recently, by settlers from the Lao-Hmong community which arrived in French Guyana in 1977.
VG acknowledges the French Ministry of Higher Education and Research for providing a Ph.D. fellowship. The research done at the Pasteur Institute of Guadeloupe benefited through a Ph.D. fellowship awarded to KB (partially financed by the European Social Funds provided through the Regional Council of Guadeloupe), and research grants awarded to NR by the International Network of the Pasteur Institutes, Direction of International Affairs, Institut Pasteur, Paris. JFG acknowledges the Institut de Recherche pour le Développement and the Centre National de la Recherche Scientifique for their financial, administrative and technical support. Thanks are due to Dr. Andrew Francis for its help on applying the method of cluster-graph to our data set and for very useful comments, and to Drs. P. Duval and B. Maubert (Pasteur Institute of French Guiana) for providing samples. We are also grateful to Drs. N. Quintard, B. Cottrelle, D. Louvel, B. Carme, J. Thonnon, R. Pradineau, G. Guillot, and A.-M. Bourbigot in French Guiana for their help during epidemiological investigations around clustered cases.
- Hopewell PC: Tuberculosis control: how the world has changed since 1990. Bulletin of the World Health Organization. 2002, 80 (6): 427-PubMedPubMed CentralGoogle Scholar
- Barnes PF, El-Hajj H, Preston-Martin S, Cave MD, Jones BE, Otaya M, Pogoda J, Eisenach KD: Transmission of tuberculosis among the urban homeless. JAMA. 1996, 275: 305-307. 10.1001/jama.275.4.305.View ArticlePubMedGoogle Scholar
- Lukacs J, Tubak V, Mester J, David S, Bartfai Z, Kubica T, Niemann S, Somoskovi A: Conventional and Molecular Epidemiology of Tuberculosis in Homeless Patients in Budapest, Hungary. J Clin Microbiol. 2004, 42 (12): 5931-5934. 10.1128/JCM.42.12.5931-5934.2004.View ArticlePubMedPubMed CentralGoogle Scholar
- Reichler MR, Bur S, Reves R, Mangura B, Thompson V, Ford J, Castro KG: Results of testing for human immunodeficiency virus infection among recent contacts of infectious tuberculosis cases in the United States. Int J Tuberc Lung Dis. 2003, 7 (12 Suppl 3): S471-8.PubMedGoogle Scholar
- Schaaf HS, Krook S, Hollemans DW, Warren RM, Donald PR, Hesseling AC: Recurrent culture-confirmed tuberculosis in human immunodeficiency virus-infected children. Pediatr Infect Dis J. 2005, 24 (8): 685-691. 10.1097/01.inf.0000172933.22481.36.View ArticlePubMedGoogle Scholar
- Monot M, Honore N, Garnier T, Araoz R, Coppée JY, Lacroix C, Sow S, Spencer JS, Truman RW, Williams DL, Gelber R, Virmond M, Flageul B, Cho S-N, Ji B, Paniz-Mondolfi A, Convit J, Young S, Fine PE, Rasolofo V, Brennan PJ, Cole ST: On the origin of leprosy. Science. 2005, 308 (5724): 1040-1042. 10.1126/science/1109759.View ArticlePubMedGoogle Scholar
- Small P, Behr M: Molecular fingerprinting of Mycobacterium tuberculosis: how can it help the clinician?. Clinical Infectious Diseases. 1997, 25: 806-810. 10.1086/515550.View ArticlePubMedGoogle Scholar
- Crawford JT, Braden CR, Schable BA, Onorato IM: National Tuberculosis Genotyping and Surveillance Network: design and methods. Emerg Infect Dis. 2002, 8 (11): 1192-1196.View ArticlePubMedPubMed CentralGoogle Scholar
- Moström P, Gordon M, Sola C, Ridell M, Rastogi N: Methods used in the molecular epidemiology of tuberculosis. Clin Microbiol Infect. 2002, 8 (11): 694-704. 10.1046/j.1469-0691.2002.00460.x.View ArticlePubMedGoogle Scholar
- Alland D, Kalkut G, Moss A, McAdam R, Hahn JA, Bosworth W, Drucker E, Bloom BR: Transmission of tuberculosis in New York City. An analysis by DNA fingerprinting and conventional epidemiologic methods. N Engl J Med. 1994, 330 (24): 1710-1716. 10.1056/NEJM199406163302403.View ArticlePubMedGoogle Scholar
- Small PM, Hopewell PC, Singh SP, Paz A, Parsonnet J, Ruston DC, Schecter GF, Daley CL, Schoolnik GK: The epidemiology of tuberculosis in San Francisco: A population-based study using conventional and molecular methods. N Engl J Med. 1994, 330: 1703-1709. 10.1056/NEJM199406163302402.View ArticlePubMedGoogle Scholar
- Seidler A, Nienhaus A, Diel R: The transmission of tuberculosis in the light of new molecular biological approaches. Occup Environ Med. 2004, 61 (2): 96-102. 10.1136/oem.2003.008573.View ArticlePubMedPubMed CentralGoogle Scholar
- Kulaga S, Behr MA, Schwartzman K: Genetic fingerprinting in the study of tuberculosis transmission. CMAJ. 1999, 161 (9): 1165-1169.PubMedPubMed CentralGoogle Scholar
- Tanaka MM, Francis AR: Methods of quantifying and visualising outbreaks of tuberculosis using genotypic information. Infect Genet Evol. 2005, 5 (1): 35-43. 10.1016/j.meegid.2004.06.001.View ArticlePubMedGoogle Scholar
- Guernier V, Guegan JF, Deparis X: An evaluation of the actual incidence of tuberculosis in French Guiana using a capture-recapture model. Microbes Infect. 2006, 8 (3): 721-727. 10.1016/j.micinf.2005.09.009.View ArticlePubMedGoogle Scholar
- Filliol I, Ferdinand S, Negroni L, Sola C, Rastogi N: Molecular typing of Mycobacterium tuberculosis based on variable number of tandem DNA repeats used alone and in association with spoligotyping. J Clin Microbiol. 2000, 38 (7): 2520-2524.PubMedPubMed CentralGoogle Scholar
- Sola C, Filliol I, Legrand E, Lesjean S, Locht C, Supply P, Rastogi N: Genotyping of the Mycobacterium tuberculosis complex using MIRUs: association with VNTR and spoligotyping for molecular epidemiology and evolutionary genetics. Infect Genet Evol. 2003, 3 (2): 125-133. 10.1016/S1567-1348(03)00011-X.View ArticlePubMedGoogle Scholar
- Supply P, Allix C, Lesjean S, Cardoso-Oelemann M, Rüsch-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: Proposal for standardization of optimized mycobacterial interspersed repetitive unit-variable-number tandem repeat typing of Mycobacterium tuberculosis. J Clin Microbiol. 2006, 44 (12): 4498-4510. 10.1128/JCM.01392-06.View ArticlePubMedPubMed CentralGoogle Scholar
- van Soolingen D, de Haas PE, Hermans PWM, van Embden JDA: DNA fingerprinting of Mycobacterium tuberculosis. Methods Enzymol. 1994, 235: 196-205.View ArticlePubMedGoogle Scholar
- Kamerbeek J, Schouls L, Kolk A, van Agterveld M, van Soolingen D, Kuijper S, Bunschoten A, Molhuizen H, Shaw R, Goyal M, van Embden J: Simultaneous detection and strain differentiation of Mycobacterium tuberculosis for diagnosis and epidemiology. J Clin Microbiol. 1997, 35: 907-914.PubMedPubMed CentralGoogle Scholar
- Goguet de la Salmoniere YO, Li HM, Torrea G, Bunschoten A, van Embden J, Gicquel B: Evaluation of spoligotyping in a study of the transmission of Mycobacterium tuberculosis. J Clin Microbiol. 1997, 35: 2210-2214.PubMedPubMed CentralGoogle Scholar
- Brudey K, Filliol I, Ferdinand S, Guernier V, Duval P, Maubert P, Sola S, Rastogi N: Long-term population-based genotyping study of Mycobacterium tuberculosis complex isolates in the French Departments of the Americas. J Clin Microbiol. 2006, 44 (1): 183-191. 10.1128/JCM.44.1.183-191.2006.View ArticlePubMedPubMed CentralGoogle Scholar
- Frothingham R, Meeker-O'Connell WA: Genetic diversity in the Mycobacterium tuberculosis complex based on variable numbers of tandem DNA repeats. Microbiology. 1998, 144: 1189-1196.View ArticlePubMedGoogle Scholar
- Brudey K, Driscoll JR, Rigouts L, Prodinger WM, Gori A, Al-Hajoj SA, Allix C, Aristimuño 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, Narvskaïa O, Ngeow YF, Naumann L, Niemann S, Parwati I, Rahim Z, Rasolofo-Razanamparany V, Rasolonavalona T, Rossetti ML, Rüsch-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: Mycobacterium tuberculosis complex genetic diversity: mining the fourth international spoligotyping database (SpolDB4) for classification, population genetics, and epidemiology. BMC Microbiol. 2006, 6: 23-10.1186/1471-2180-6-23.View ArticlePubMedPubMed CentralGoogle Scholar
- Enright AJ, Ouzounis CA: BioLayout – an automatic graph layout algorithm for similarity visualization. Bioinformatics. 2001, 17 (9): 853-854. 10.1093/bioinformatics/17.9.853.View ArticlePubMedGoogle Scholar
- Goldovsky L, Cases I, Enright AJ, Ouzounis CA: BioLayoutJava: versatile network visualisation of structural and functional relationships. Applied Bioinformatics. 2005, 4 (1): 71-74. 10.2165/00822942-200504010-00009.View ArticlePubMedGoogle Scholar
- Glynn JR, Vynnycky E, Fine PE: Influence of sampling on estimates of clustering and recent transmission of Mycobacterium tuberculosis derived from DNA fingerprinting techniques. Am J Epidemiol. 1999, 149: 366-371.View ArticlePubMedGoogle Scholar
- Murray M: Sampling bias in the molecular epidemiology of tuberculosis. Emerg Infect Dis. 2002, 8 (4): 363-369.View ArticlePubMedPubMed CentralGoogle Scholar
- Rosenberg NA, Tsolaki AG, Tanaka MM: Estimating change rates of genetic markers using serial samples: applications to the transposon IS 611O in Mycobacterium tuberculosis. Theor Popul Biol. 2003, 63 (4): 347-363. 10.1016/S0040-5809(03)00010-8.View ArticlePubMedGoogle Scholar
- Filliol I, Driscoll JR, van Soolingen D, Kreiswirth BN, Kremer K, Valétudie G, Dang DA, Barlow R, Banerjee D, Bifani PJ, Brudey K, Cataldi A, Cooksey RC, Cousins DV, Dale JW, Dellagostin OA, Drobniewski F, Engelmann G, Ferdinand S, Gascoyne-Binzi D, Gordon M, Gutierrez MC, Haas WH, Heersma H, Kassa-Kelembho E, Ho ML, Makristathis A, Mammina C, Martin G, Moström P, Mokrousov I, Narbonne V, Narvskaya O, Nastasi A, Niobe-Eyangoh SN, Pape JW, Rasolofo-Razanamparany V, Ridell M, Rossetti ML, Stauffer F, Suffys PN, Takiff H, Texier-Maugein J, Vincent V, de Waard JH, Sola C, Rastogi N: Snapshot of moving and expanding clones of Mycobacterium tuberculosis and their global distribution assessed by spoligotyping in an international study. J Clin Microbiol. 2003, 41 (5): 1963-1970. 10.1128/JCM.41.5.1963-1970.2003.View ArticlePubMedPubMed CentralGoogle Scholar
- Warren RM, Streicher EM, Sampson SL, van der Spuy GD, Richardson M, Nguyen D, Behr MA, Victor TC, van Helden PD: Microevolution of the Direct Repeat region of Mycobacterium tuberculosis: implications for interpretation of spoligotyping data. J Clin Microbiol. 2002, 40 (12): 4457-4465. 10.1128/JCM.40.12.4457-4465.2002.View ArticlePubMedPubMed CentralGoogle Scholar
- Filliol I, Driscoll JR, Van Soolingen D, Kreiswirth BN, Kremer K, Valétudie G, Anh DD, Barlow R, Banerjee D, Bifani PJ, Brudey K, Cataldi A, Cooksey RC, Cousins DV, Dale JW, Dellagostin OA, Drobniewski F, Engelmann G, Ferdinand S, Gascoyne-Binzi D, Gordon M, Gutierrez MC, Haas WH, Heersma H, Källenius G, Kassa-Kelembho E, Koivula T, Ly HM, Makristathis A, Mammina C, Martin G, Moström P, Mokrousov I, Narbonne V, Narvskaya O, Nastasi A, Niobe-Eyangoh SN, Pape JW, Rasolofo-Razanamparany V, Ridell M, Rossetti ML, Stauffer F, Suffys PN, Takiff H, Texier-Maugein J, Vincent V, De Waard JH, Sola C, Rastogi N: Global distribution of Mycobacterium tuberculosis spoligotypes. Emerg Infect Dis. 2002, 8: 1347-1349.View ArticlePubMedPubMed CentralGoogle Scholar
- Singh UB, Suresh N, Bhanu NV, Arora J, Pant H, Sinha S, Aggarwal RC, Singh S, Pande JN, Sola C, Rastogi N, Seth P: Predominant tuberculosis spoligotypes, Delhi, India. Emerg Infect Dis. 2004, 10: 1138-1142.View ArticlePubMedPubMed CentralGoogle Scholar
- Sola C, Filliol I, Guttierez C, Mokrousov I, Vincent V, Rastogi N: Spoligotype database of Mycobacterium tuberculosis: biogeographical distribution of shared types and epidemiological and phylogenetic perspectives. Emerg Infect Dis. 2001, 7: 390-396.View ArticlePubMedPubMed CentralGoogle Scholar
- Kremer K, van Soolingen D, Frothingham R, Haas WH, Hermans PW, Martín C, Palittapongarnpim P, Plikaytis BB, Riley LW, Yakrus MA, Musser JM, van Embden JD: Comparison of methods based on different molecular epidemiological markers for typing of Mycobacterium tuberculosis complex strains: interlaboratory study of discriminatory power and reproducibility. J Clin Microbiol. 1999, 37: 2607-2618.PubMedPubMed CentralGoogle Scholar
- Barret J, (Dir): Atlas illustré de la Guyane. 2001, Limoges: IRD, 2Google Scholar
- Barrangou R, Fremaux C, Deveau H, Richards M, Boyaval P, Moineau S, Romero DA, Horvath P: CRISPR provides acquired resistance against viruses in prokaryotes. Science. 2007, 315 (5819): 1709-1712. 10.1126/science.1138140.View ArticlePubMedGoogle Scholar
- Filliol I, Motiwala AS, Cavatore M, Qi W, Hazbón 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: 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. 2006, 188 (2): 759-772. 10.1128/JB.188.2.759-772.2006.View ArticlePubMedPubMed CentralGoogle Scholar
- Gutacker MM, Mathema B, Soini H, Shashkina E, Kreiswirth BN, Graviss EA, Musser JM: Single-nucleotide polymorphism-based population genetic analysis of Mycobacterium tuberculosis strains from 4 geographic sites. J Infect Dis. 2006, 193 (1): 121-128. 10.1086/498574.View ArticlePubMedGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2334/8/46/prepub
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.