Previous Article | Next Article ![]()
Journal of Clinical Microbiology, April 2005, p. 1587-1593, Vol. 43, No. 4
0095-1137/05/$08.00+0 doi:10.1128/JCM.43.4.1587-1593.2005
Copyright © 2005, American Society for Microbiology. All Rights Reserved.
M'Hammed Sarih,2,
Fatima Jouda,3
Edmond Godfroid,4
Lise Gern,3
Ali Bouattour,1
Guy Baranton,5 and
Daniele Postic5*
Institut Pasteur, Tunis, Tunisia,1 Institut Pasteur du Maroc, Casablanca, Morocco,2 Institut de Zoologie, Université de Neuchâtel, Neuchâtel, Switzerland,3 Université Libre de Bruxelles, Gosselies, Belgium,4 Institut Pasteur, Paris, France5
Received 29 March 2004/ Returned for modification 17 October 2004/ Accepted 7 November 2004
|
|
|---|
|
|
|---|
In contrast to these species, little data are available concerning Borrelia lusitaniae. Only sporadic isolates have been characterized from ticks collected in the Czech Republic, Moldavia, Ukraine, and Belarus (20), in Slovakia (11), in Spain (3), in Poland (25), in France (30), in Switzerland (15-17), or in Turkey (13), and these isolates did not show any significant diversity. However, a study conducted in Portugal revealed a large diversity among B. lusitaniae strains isolated from a local tick population (6). Moreover, it is noteworthy that, among the three first Portuguese isolates identified as B. lusitaniae, two different types were demonstrated (20, 26). Hence, these studies suggested a higher diversity among B. lusitaniae around the Mediterranean basin. However, all these studies only focused on the small ribosomal rrf-rrl spacer and, therefore, the global diversity of this species remained unknown.
Recently, we reported that B. lusitaniae was, by far, the most prevalent species, ranging from 96.6 to 100% of the B. burgdorferi sensu lato species identified in Tunisia and Morocco (35, 44, 46). The report of some Lyme borreliosis cases in North Africa (1, 27) and the recent isolation in Portugal of B. lusitaniae from the skin of a human patient with chronic lesions (5) prompted us to characterize this species as a potential human pathogen. Therefore, we investigated the molecular characterization and the genetic diversity of B. lusitaniae isolated from free-living Ixodes ricinus ticks in Tunisia and Morocco. For these purposes we studied two single loci (ospC and ospA) and the noncoding intergenic spacer (rrf-rrl), as well as the whole genome.
|
|
|---|
Samplings at these sites were conducted between 2000 and 2003. Questing ticks were collected by blanket dragging the vegetation, as previously reported (35, 44). Only adult ticks were used in this study.
Borrelia culture and DNA extraction. Live field-collected I. ricinus ticks were dipped in 70% ethanol, rinsed in distilled water, and then ground up and inoculated into 2 ml of BSK-H medium (Sigma) supplemented with 7% gelatin and 1% antibiotic mixture for Borrelia (Sigma). Tightly capped tubes were incubated at 34°C and examined once a week by dark-field microscopy. Subcultures were performed without antibiotics and into a larger volume. For pulsed-field gel electrophoresis (PFGE) only isolates maintained in pure cultures were used, whereas for PCR experiments contaminated cultures and/or cultures containing few spirochetes were used to extract DNA by boiling at 100°C for 10 min. Strains and DNA used in this study are shown in Table 1.
|
View this table: [in a new window] |
TABLE 1. B. lusitaniae from Maghrebin countries used in this work and the targets studied for each isolatea
|
PFGE. Previously described procedures were used for the preparation of high-molecular-weight genomic DNAs and for PFGE (4, 9). Rare-cutting endonucleases SacII, BssHII, and MluI were used to differentiate isolates. Lambda concatemers (monomer size, 48.5 kbp) purchased from New England Biolabs, Beverly, Mass., were used as a size control. All fragments on the gels resulting from restrictions were transformed as binary data, where each fragment was scored as 1 (present) or 0 (absent). These data yielded a binary matrix that was used to draw a phylogenetic tree.
Analysis of ospC diversity by SSCP-PCR. To screen for the genetic diversity of the ospC gene, we used a single-strand conformation polymorphism (SSCP) analysis. SSCP is a screening method based on the secondary structure of a single-stranded DNA fragment. A 277-bp fragment in the ospC gene was amplified by PCR using the primer set F-SC3 and R-OspC9 under the conditions previously described (19). Briefly, after denaturation, PCR samples were loaded on a polyacrylamide gel. Electrophoresis was conducted in a temperature-controlled electrophoresis system (GenePhor; Amersham Pharmacia Biotech) at 6°C with a first run at 600 V, 25 mA, and 15 W for 10 min and then at 600 V, 37 mA, and 21 W for 2.5 h. Gels were revealed by silver staining according to the manufacturer's instructions (Plus One DNA silver staining kit; Amersham Pharmacia Biotech).
Sequencing and phylogenetic analysis of sequences. Intergenic spacer PCR products were sequenced by Genome Express (Meylan, France). Amplification products from the ospA gene were sequenced using the BigDye Terminator cycle sequencing kit (Applied Biosystems) as described in reference 12. Primers used in this study were direct and reverse primers 1 and 2, as previously described (12).
Sequences were aligned by using Clustal V software and manually by using VSM software and then analyzed by either the unweighted pair group with mathematical average (UPGMA) (36) or the neighbor-joining (NJ) method (34). Phylogenetic trees were drawn with Mega software (18). Sequences of North African B. lusitaniae isolates were compared to additional sequences from B. lusitaniae and from other species available in data banks.
|
|
|---|
![]() View larger version (52K): [in a new window] |
FIG. 1. UPGMA-rooted tree obtained with MEGA software from the ribosomal rrf-rrl spacer sequences. Sequences issued from data banks are followed by the sequence accession number. AF, Africa; E, Europe; AS, Asia; Mo, Morocco; Tu, Tunisia; Po, Portugal; Tk, Turkey; Sp, Spain; Mol, Moldavia.
|
Six different restriction patterns were observed, as shown in Table 2. A large majority of MseI patterns (n = 456) were similar to that of the type strain PotiB2T. Sequences of strains exhibiting a similar restriction pattern differed from the PotiB2T sequence by only one or two nucleotide deletions. Interestingly, the MseI and DraI patterns of strain M3-22 were indistinguishable from those of strain PotiB2T on the gels. However, sequencing revealed an insertion of 15 nucleotides containing an MseI and DraI restriction site leading to an additional 16-bp-long fragment not visible on the gels. The remaining sequence was identical to that of strain PotiB2T. Consequently, this pattern may be underscored when no sequencing is performed. Four Tunisian strains were identical to strain TD12 (44). Two previously undescribed patterns were identified and named MT26 and MA3-13 (Table 2). Five Moroccan strains exhibited the MT26 pattern, and two strains exhibited the MA3-13 pattern. These patterns have been confirmed by sequencing as belonging to B. lusitaniae, close to PotiB2T. Only 11 strains showed the typical pattern of strain PotiB3, one of the first genotypes of B. lusitaniae described in Portugal.
|
View this table: [in a new window] |
TABLE 2. Restriction fragment analysis of rrl-rrf amplicons
|
Whole genome analyzed by PFGE. Thirty-three North African B. lusitaniae isolates were analyzed by PFGE after DNA restriction by three enzymes, SacII, BssHII (Fig. 2), and MluI, and compared to the Portuguese strains PotiB2T and PotiB3. Only discrete restriction fragments of >50 kb were included in the analysis.
![]() View larger version (101K): [in a new window] |
FIG. 2. BssHII macrorestriction patterns of B. lusitaniae strains obtained with the Bio-Rad apparatus with the pulse time ramped from 2 to 70 s for 30 h. (A) Moroccan strains. The lanes contained lambda concatemers (lanes 1 and 11), PotiB2 (lane 2), MT1 (lane 3), MT7 (lane 4), MT18 (lane 5), MT20 (lane 6), MT24 (lane 7), MT25 (lane 8), MT27 (lane 9), and PotiB3 (lane 10). (B) Tunisian strains. The lanes contained lambda concatemers (lanes 1 and 13), PotiB2 (lane 2), PotiB3 (lane 3), TT840 (lane 4), TT867 (lane 5), TT894 (lane 6), TT899 (lane 7), TT908 (lane 8), TT916 (lane 9), TT918 (lane 10), TT919 (lane 11), and TT920 (lane 12).
|
Although restriction by SacII led to the higher polymorphism, it did not allow recognition of the PotiB3 genotype, which can be easily distinguished from PotiB2 after restriction by BssHII and MluI.
All PFGE data yielded a matrix of 35 strains by 25 characters that was subjected to a phylogenetic analysis by NJ and UPGMA methods. The pulsotype PotiB3 clustered separately from other B. lusitaniae pulsotypes resolved in two large divisions. The largest division contained isolates that segregated with PotiB2T (Fig. 3). As shown in Fig. 3, 33 Moroccan and Tunisian strains were equally distributed into 19 alleles.
![]() View larger version (22K): [in a new window] |
FIG. 3. Phylogenetic analysis of polymorphism after DNA restriction with MluI, BssHII, and SacII and resolution by PFGE. All fragments were transformed into a binary matrix used to generate a tree by the UPGMA method.
|
ospC gene mobility classes. Thirty-three DNAs, 16 from ticks collected in Morocco and 17 from ticks collected in Tunisia, were analyzed by SSCP-PCR. A consistent diversity was found. With the exception of two couples of Tunisian strains that each shared a similar pattern (data not shown), each DNA studied was characterized by a specific SSCP pattern (Fig. 4).
![]() View larger version (125K): [in a new window] |
FIG. 4. ospC SSCP patterns of Moroccan B. lusitaniae strains. Lane 1, MT2; lane 2, MT4; lane 3, MT8; lane 4, MT9; lane 5, MT10; lane 6, MT11; lane 7, MT12; lane 8, MT14; lane 9, MT16; lane 10, B. lusitaniae PotiB1; lane 11, B. lusitaniae PotiB2T; lane 12, B. lusitaniae PotiB3; lane 13, MT17; lane 14, MT20; lane 15, MT22; lane 16, MT24; lane 17, MT27; lane 18, MT13, lane 19, MT26.
|
|
|
|---|
Most authors have analyzed the global genetic variability of B. burgdorferi sensu lato by the analysis of the macrorestriction polymorphism with one enzyme (MluI). If we considered the MluI restriction polymorphism, only three patterns were recorded among 35 strains. By using three enzymes, we observed a global diversity similar to that reported for B. burgdorferi sensu stricto or B. garinii after restriction by only one enzyme, confirming a low diversity (4, 9, 10, 23). Our results did not evidence any clear clustering among strains isolated in Tunisia or Morocco (Fig. 3).
In order to obtain a higher level of discrimination among B. lusitaniae strains, we also focused on two loci, ospA and ospC. The former is subjected to an adaptive pressure, whereas the latter is subjected to a diversifying frequency-dependent selection (2, 8), with these two different mechanisms leading to distinct ranges of variability. The ospC diversity should be due to an immune pressure from locally prominent hosts, preventing a reinfection by an antigenically related strain. Therefore, whatever the host or vector spectrum, the diversity is expected to be high. This was confirmed for B. lusitaniae in North Africa. The genetic variability of the ospC gene of pathogenic B. burgdorferi sensu lato strains has been extensively studied (14, 21, 37, 39). The diversity of the ospC gene is mostly acquired by lateral transfer between species or within species. However, given both the restricted geographic distribution of B. lusitaniae and the scarcity of isolates from other species in North Africa, incorporation of sequences from other B. burgdorferi sensu lato species is unlikely. To clarify this point, it would be worthwhile to extend the study to B. lusitaniae strains from other geographic areas. ospA gene sequence analyses allow identification of species within B. burgdorferi sensu lato (23, 42, 43), except for B. valaisiana. It is noteworthy that the ospA sequences from Portuguese Poti strains clustered separately, although close to North African sequences (data not shown). No precise correlation could be established between segregation of ospA and macrorestriction patterns. The relative level of conservation of ospA sequences from B. lusitaniae could be attributed to the restricted geographic extension and to the interaction with only one vector species and possibly one predominant reservoir host. This fact contrasts with the diversity reported for B. valaisiana, which is more widely distributed and is associated with several tick species, requiring adaptation to a broad range of vectors (12, 22, 24, 32, 40).
Our analysis of the restriction polymorphism of the small rrf-rrl spacer emphasized a high monomorphism that was confirmed by the conservation of sequences within B. lusitaniae isolates from North Africa. Fourteen alleles were recorded for the whole B. lusitaniae sequences. Among these 14 alleles, 11 contained 14 available sequences from Portuguese isolates, whereas only 7 alleles contained 79 sequences from North African isolates. Only three alleles present among North African isolates were specific to this continent and had no counterpart in Portugal. Therefore, a considerable rrf-rrl sequence heterogeneity seems to characterize Portuguese strains, in contrast to North African strains. However, De Michelis et al. (6) reported that different alleles were not equally represented, with 50% of Portuguese strains belonging to allele I (strain GT058). It is noteworthy that 80.7% of North African sequences also belonged to this allele (Fig. 1), as did the Moldavian strain Ir345 (20, 29). The description of three new alleles from North African strains does not significantly increase the global genetic diversity previously known for the species B. lusitaniae on this locus. Given the large number of sequences determined, it is likely that the global diversity of the species has been assessed. It is noteworthy that no sequence from strains isolated in North Africa fell in the same allele as strain PoHL isolated from a Portuguese patient.
It must be emphasized that the genetic heterogeneity involving the ribosomal spacer was essentially acquired by insertion or deletion of short DNA fragments (data not shown). Very few mutations were observed along the sequences. If we consider that this locus is not subjected to any selection pressure and mostly reflects the molecular evolutionary clock, we could hypothesize that either North African B. lusitaniae strains recently evolved from a common ancestor or only a few prevalent genotypes have successfully settled in North Africa. Considering that the whole diversity of the species is present in Portugal, it is likely that B. lusitaniae has evolved for a considerable time in Portugal and has been recently introduced in North Africa. Therefore, the analysis of a chromosomal locus makes obvious a founding event leading to the recent emergence of some North African clones sampled from the Portuguese diversity.
Our findings with a large population of B. lusitaniae strains from North Africa support the conclusion that a high proportion of this population belongs to a limited number of closely related genotypes. This is in accordance with the clonal structure of the Borrelia population as previously demonstrated (7). Some cases of Lyme borreliosis have been reported in Tunisia, based on clinical and serological features (1), and two cases of facial palsy have been reported in Morocco (27). Moreover, a recent study indicated that the inoculation of B. lusitaniae into susceptible mice induced pathological features (45). Finally, B. lusitaniae was recently isolated, for the first time, from a chronic skin lesion in a Portuguese patient (5). If this species is actually pathogenic for humans, the risk of Lyme borreliosis must be underlined in North Africa. However, the pathogenicity could be restricted to distinct clones which have not yet settled in North Africa. Additional studies are needed to clarify the evolutionary events leading to the emergence of a restricted population of bacteria in a given geographical area.
We thank K. Dellagi and M. Hassar, Directors of the Pasteur Institute of Tunisia and Morocco, respectively, for their constant support. The language of the manuscript was revised by Iain Old.
H.Y. and M'H.S. contributed equally to this work. ![]()
|
|
|---|
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»