Previous Article | Next Article ![]()
Journal of Clinical Microbiology, April 2005, p. 1546-1551, Vol. 43, No. 4
0095-1137/05/$08.00+0 doi:10.1128/JCM.43.4.1546-1551.2005
Copyright © 2005, American Society for Microbiology. All Rights Reserved.
Mycobacteriology Unit, Institute of Tropical Medicine, Antwerp,1 Laboratory of Microbiology,2 BCCM/LMG Culture Collection, University Ghent, Ghent, Belgium,3 Laboratoire des Biopuces,4 Laboratoire des Mécanismes Moléculaires de la Pathogenèse Bactérienne, INSERM U629, Institut Pasteur de Lille, Lille Cedex, France5
Received 5 October 2004/ Returned for modification 19 November 2004/ Accepted 6 December 2004
|
|
|---|
|
|
|---|
Other molecular techniques have been applied with relatively modest gains in further unraveling aspects of the epidemiology of the disease. The genetic relationships among M. ulcerans strains from different geographical regions have been inferred from comparisons of insertion sequences and 16S rRNA, housekeeping, and structural gene sequences. These markers have indicated restricted genetic variation in M. ulcerans, especially among strains from the same geographic origin (5, 16, 21, 22), strongly hampering the determination of the sources of infection at a local level.
Tandem repeat (TR) DNA sequences are important sources of polymorphism in the genome of many eukaryotes and prokaryotes. Genomic regions showing polymorphism due to different numbers of TR motifs in different strains or individuals are described as variable number tandem repeat (VNTR) loci. Allele-length polymorphism associated with VNTR loci can be indexed by PCR amplification. Especially in highly monomorphic and clonal species such as Mycobacterium tuberculosis (8, 13, 17, 19, 20, 25), Bacillus anthracis (9), and Yersinia pestis (10), VNTR arrays represent polymorphism islets in the otherwise highly conserved genomic backbones and have therefore been very useful in differentiating and studying the genetic relationships of strains of these species.
The main objective of this study was to investigate the existence of VNTR loci in M. ulcerans and to assess their utility for the discrimination of strains of this organism.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Strain information
|
Location of tandem repeats in M. marinum and M ulcerans. The unassembled M. marinum genome sequence available as of August 2003 at the Sanger Centre website (http://www.sanger.ac.uk/Projects/M_marinum/) was screened for the presence of TR loci by using the Tandem Repeat Finder program of the Department of Biomathematical Science of Mount Sinai School of Medicine (http://c3.biomath.mssm.edu/trf.html).
BLAST searches were performed on the M. ulcerans genome sequences available as of August 2003 on the BURULIST Web server of the Pasteur Institute (http://genopole.pasteur.fr/Mulc/BuruList.html) to establish the presence of homologous TR loci in the M. ulcerans genome sequence. The sequences were also aligned with the M. tuberculosis genome sequence available on the TUBERCULIST Web server (http://genolist.pasteur.fr/TubercuList/).
VNTR PCR. The primers for PCR amplification (Table 2) were designed on the basis of the flanking sequences of each TR locus of M. marinum by the use of Oligo 5.0 software (National Biosciences). The PCR was carried out using a Hotstar Taq DNA polymerase kit (QIAGEN). Sample DNA (purified or crude extract) (3 µl) was added to 27 µl of a PCR mix containing 15.5 µl of water, 3 µl of 10x PCR buffer (containing 1.5 mM MgCl2 at the final concentration), 0.4 µM of each primer, 0.2 mM of each deoxynucleotide triphosphate (Roche), 1x Q-solution, and 0.1 µl of HotstarTaq DNA polymerase (0.5 U).
|
View this table: [in a new window] |
TABLE 2. Characteristics of VNTR loci in M. marinum and M. ulcerans
|
Genetic distance analysis. The analysis of the genetic relationships on the basis of the VNTR amplicon size profiles was performed using the neighbor-joining algorithm included in the PAUP software package (D. Swofford, phylogenetic analysis using parsimony, 4.0 beta version; Sinauer Associates, Inc., Sunderland, Mass.).
|
|
|---|
The loci were given numerical labels 1 to 19. A total of 13 of them were found to be polymorphic in a small collection of M. marinum strains. Table 2 shows their locus information. Most of them, with the exception of locus 18, contained repeats with complete consensus sizes plus one partial repeat, with lengths ranging from half to nearly complete units. An alignment of the sequences of each locus with the homologous sequences in the M. ulcerans genome revealed complete identity of TR unit sequences among both mycobacterial species at loci 8, 9, and 14. In the other loci, there were nucleotide deletions of from 2 to 17 bp and/or from 2 to 14 nucleotide substitutions in the M. ulcerans units compared to those of M. marinum. The sequences flanking the repeat sequences were found to be highly conserved in both species and homologous to M. tuberculosis open reading frames encoding products of diverse functional categories (Table 2). Most of the corresponding TR sequences were found in intergenic positions in these regions. At least one of these homologous regions (locus 13) corresponds to a known VNTR locus in M. tuberculosis, namely, the senX3-regX3 locus (25).
Typing of strains and cluster analysis. The size polymorphism of the 19 TR loci was analyzed by PCR using a collection of 23 M. ulcerans isolates, 5 M. marinum isolates, and an intermediate isolate. This latter strain contains a low copy number of the M. ulcerans-specific insertion sequence IS2404, while many other of its phenotypic and genetic characteristics are similar to those of M. marinum (4). The amplification was robust and yielded reproducible results for all loci and all isolates. A total of 13 of these loci were polymorphic among the M. marinum control strains; nine were also polymorphic in M. ulcerans (Fig. 1). (Only the results obtained for the nine polymorphic loci in both species are shown).
![]() View larger version (35K): [in a new window] |
FIG. 1. Dendrogram showing the genetic relationships of the M. marinum and M. ulcerans isolates determined using nine VNTR loci. The dendrogram was built using the neighbor-joining algorithm and ITM 1717 and 7732 as outgroups, as described in Material and Methods. The VNTR amplicon sizes and the corresponding genotype designations are indicated at the right. The linkage distance scale and the VNTR locus numbers are indicated at the bottom. MM, M. marinum; MU, M. ulcerans; MI, intermediate isolate.
|
Eight different M. ulcerans genotypes, five M. marinum genotypes, and a unique intermediate strain genotype were found by VNTR typing (Fig. 1). Although the test panel is of limited size, the results suggest differential variability among the nine loci. Loci 18 and 19 were the most variable loci, with a total of seven and six alleles overall, respectively. Locus 19 was the most discriminatory at the intraspecies level, showing four different alleles in M. ulcerans and four in M. marinum (Fig. 2A and B, respectively). At the other extremity, locus 6 was among the least discriminatory, with only three alleles in total, of which two and three were found in M. ulcerans and M. marinum, respectively. Nevertheless, only loci 6 and 14 could be used to differentiate between the Surinam and French Guyana M. ulcerans isolates (Fig. 1).
![]() View larger version (11K): [in a new window] |
FIG. 2. PCR analysis of VNTR locus 19. (A) M. ulcerans and the intermediate strain. Lanes: L, 100-bp DNA ladder; 1, ITM 5142; 2, ITM 5147; 3, ITM 94-1324; 4, ITM 9537; 5, ITM 94-1331; 6, ITM 94-1328; 7, ITM 98-912; 8, ITM 8756; 9, ITM 842; 10, ITM 7922; 11, ITM 5114; 12, ITM 5143; 13, ITM 00-1026; 14, 50-bp ladder (locus 19); 15, ITM 96-658; 16, ITM 97-111; 17, ITM 97-483; 18, ITM 97-104; 19, ITM 5152; 20, ITM 5155; 21, ITM 97-610; 22, ITM 97-648; 23, ITM 97-680; 24, ITM 94-662; 25, ITM 94-511; 26, negative control. (B) M. marinum. Lanes: L, 100-bp DNA ladder; 27, ITM 8022; 28, ITM 1717; 29, ITM 94-56; 30, ITM 1320; 31, ITM 7732; 32, negative control.
|
A first group in the dendrogram included one Australian genotype (alias Victoria) (22), two genotypes designated Southeast Asian and Papua New Guinea of M. ulcerans, and one African genotype. The Australian genotype included two isolates from Australia (Victoria), while the Southeast Asian genotype comprised the third isolate from Australia (Queensland), the Malaysian isolate, and one Papua New Guinea isolate (ITM 94-1331). The genotype of the other Papua New Guinea isolate (ITM 9537) was distinct from the Southeast Asian genotype by only one locus. Strikingly, the African genotype included the 11 M. ulcerans isolates from the seven different African countries tested; these isolates did not display any differences among the nine VNTR loci investigated. Similarly, the isolates from Japan and China had identical profiles at all the nine loci and constituted the Asian genotype. The two Mexican isolates had also identical profiles at all nine loci and were designated the Mexican genotype. This Mexican type was found to be closely related to the Asian type. The two South American isolates from Surinam and French Guyana both gave unique but closely related genotypes, referred to as the Surinam and Guyana genotypes, respectively. These strains appeared to be more distant from the other M. ulcerans strains.
The five M. marinum isolates had five unique genotypes but could be classified into two subgroups, showing some correlation with their source of isolation. These two groups included the two isolates from Belgium and the one from Nicaragua as one group and the two isolates from the United States as the other group. The isolate of the intermediate species (ITM 00-1026) was found to be distantly related to the M. ulcerans Asian-Mexican genotype.
|
|
|---|
The set of nine VNTR loci revealed eight genotypes among the 23 M. ulcerans reference isolates tested. Using largely overlapping sets of isolates, six genotypes have previously been found by multilocus sequence typing (MLST) analysis of eight structural and house-keeping genes (21) and by restriction fragment length polymorphism (RFLP) of IS2404 (6), whereas nine genotypes were found with 2426-PCR, a repetitive-element PCR targeting IS2404 and IS2606 (22) (Table 1). In contrast, amplified fragment length polymorphism (AFLP) revealed only three M. ulcerans subtypes (6). Although the different typing methods were not applied to exactly the same strain collections, these results suggests that the typing method used on the basis of the set of nine VNTR loci described here is at least as discriminatory as the previous typing methods.
VNTR typing groups M. ulcerans isolates according to their geographic origin similarly to MLST results. Both methods cluster the Malaysian and Papua New Guinean isolates ITM 94-1331 as the Southeast Asian type. VNTR analysis provided additional resolution by discriminating one Papua New Guinean (strain 11878/70 or ITM 9537) genotype, but this genotype was consistently closely related by a single-locus difference to the South East Asian genotype. Numerical analyses of the profiles generated by both methods indicate the closest evolutionary link between strains from Southeast Asia, Australia, and Africa and a comparatively more distant relation among genotypes from other regions. However, at variance with 2426-PCR results (22) but consistent with MLST results (21), multilocus VNTR analysis shows the Asian genotype (Japan-China) to be much more closely related to the Mexican genotype than to the Surinam genotype. Considering that the most reliable method for determining genetic relatedness is sequence comparison, the congruence between VNTR-based and MLST-based (21) groupings indicates that reliable phylogenetic inferences can be made from this set of markers over the scale of evolutionary divergence considered. Furthermore, such congruence is consistent with the clonal population structure of M. ulcerans apparent from MLST analysis (21).
It is particularly noteworthy that despite the resolution power of these markers, we did not find any size differences in any of nine VNTR loci tested among the 11 African isolates coming from the seven different countries. Consistent with the absence of any sequence variation among African strains identified so far (21, 22), these findings highlight the extremely high clonal homogeneity of M. ulcerans in Africa and lend support to the hypothesis of a recent distribution of the organism across this continent (21). The level of genetic conservation is not so extreme within other geographic regions. In addition to discriminating the isolates from Papua New Guinea by one locus (see above), VNTR analysis could discriminate between the isolates from French Guyana and Surinam by two loci, in contrast to IS2404 RFLP results (5) and AFLP results (6). This finding identifies for the first time two closely related but different M. ulcerans genotypes in South America. Furthermore, the difference between these two isolates and the Mexican isolates, also apparent from MLST analysis (21), was much greater (at seven loci out of nine) than that between the two South American isolates. This is the most extensive difference among M. ulcerans strains originating from a same continent observed so far by VNTR analysis and suggests that the corresponding genotypes may have arisen independently in the Americas.
Beyond this comparison with the Mexican isolates, the South American isolates appeared distantly related to the M. ulcerans strains from the other geographic regions in general. Consistently, the isolate from Surinam appeared as the most distal geographical genotype among the M. ulcerans isolates tested by MLST (21) and was distinct from M. ulcerans isolates from the other geographic regions by the sequence of the 3'-terminal end of its 16S rRNA gene, which is identical to that of M. marinum strains (16). Similarly, the isolate of the intermediate species (ITM 00-1026) was found to be distantly related to the Asian-Mexican M. ulcerans genotypes and one of the M. marinum subgroups. However, the analysis of more deep-rooted relationships should be considered with some caution because of the limited degree of confidence inherent in investigations conducted using a set presently restricted to nine VNTR loci.
In conclusion, this study supported the potential of a VNTR-based genotyping method for M. ulcerans and M. marinum. Since VNTR genotyping has also been successfully applied to the M. tuberculosis (8, 17, 19, 20, 25) and M. avium (2, 3) complexes, it is likely that they can be used for the typing of other species of the mycobacterial genus as well. This method offers several practical advantages over AFLP- or RFLP-based typing methods, such as easier DNA extraction (from boiled colonies) suitable for typing and easier interpretation of the banding patterns. The panel of nine markers identified here permitted the discrimination of M. ulcerans isolates within some but not all geographic regions and offers a resolution power comparable to or slightly higher than that of these methods and of MLST. The identification of these nine loci can therefore be considered as a first promising step towards the potential use of VNTRs as markers for molecular epidemiological studies of BU within given geographic regions. This set of loci most likely represents only a fraction of the VNTR patrimony of M. ulcerans, as numerous and diverse arrays of TRs identified in M. marinum have homologs in M. ulcerans which have not been explored yet. Conversely, it is possible that in M. ulcerans, some VNTR loci, such as in the recently identified giant plasmid encoding its macrolide toxin (24), exist that have no counterparts in the M. marinum genome and thus could not be detected by our homology-based strategy. Hence, the publication of the complete M. ulcerans genome sequence will provide an opportunity for the discovery of additional tandem repeat loci, which should enhance strain resolution. Higher resolution will be especially useful to further assess the genetic homogeneity of African M. ulcerans isolates.
This work was partly supported by a grant from the Fund for Scientific Research, Flanders (Flanders, Belgium) (F. W. O. Vlaanderen; contracts G. 0301.01 and G. 0471.03N), by the Institut Pasteur de Lille, and by INSERM. A.A. was supported by a grant from the Damien Foundation (Brussels, Belgium). P.S. is a Researcher of the Centre National de la Recherche Scientifique (CNRS).
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»