Previous Article | Next Article ![]()
Journal of Clinical Microbiology, September 2005, p. 4708-4712, Vol. 43, No. 9
0095-1137/05/$08.00+0 doi:10.1128/JCM.43.9.4708-4712.2005
Copyright © 2005, American Society for Microbiology. All Rights Reserved.
Pierre-Edouard Fournier,1,
Hiroyuki Ogata,2 and
Didier Raoult1*
Unité des Rickettsies, CNRS UMR 6020, IFR 48, Faculté de Médecine, Université de la Méditerranée, 27 Blvd. Jean Moulin, 13385 Marseille Cedex 5, France,1 Information Génomique et Structurale, CNRS UPR2589, 31 Chemin Joseph Aiguier, 13402 Marseille Cedex 20, France2
Received 1 February 2005/ Returned for modification 9 May 2005/ Accepted 10 June 2005
|
|
|---|
|
|
|---|
Present-day threats posed by R. prowazekii include natural exposure and, possibly, the use of laboratory-manipulated strains as agents of biological warfare (27). R. prowazekii has been classified on the B list of potential bioterrorism agents by the Centers for Diseases Control and Prevention (Atlanta, GA). A method of differentiating strains of R. prowazekii would be very useful in determining the cause of outbreaks that might occur.
Rickettsial species were first characterized by serotyping with antisera raised in mice and subsequently by protein analysis with sodium dodecyl sulfate-polyacrylamide gel electrophoresis (6) and pulsed-field gel electrophoresis (31). The latter method could differentiate R. prowazekii strains Breinl and Evir (14). However, since rickettsiae are strict intracellular bacteria, serotyping, sodium dodecyl sulfate-polyacrylamide gel electrophoresis, and pulsed-field gel electrophoresis are often difficult, time consuming, and expensive and lack interlaboratory reproducibility. Because rickettsiae express few remarkable phenotypic characteristics, their precise identification and phylogenic classification have mainly become dependent on the study of sequences of various genes they contain. The 16S rRNA and gltA genes, which were studied initially, have been demonstrated to lack intraspecies variability (32, 34). Other studies have focused on a gene family encoding cell surface-exposed proteins (sca genes). Among these, ompA, ompB and sca4, also named gene D, have proven useful to infer reliable phylogenetic relationships among Rickettsia spp. (18, 33, 35). However, although sca genes were more variable than the 16S rRNA and gltA genes at the interspecies level and had a certain degree of intraspecies variability, they had only limited interstrain variability (18, 33, 35). Similarly, ompA, which is one of the most discriminatory genes among spotted fever group rickettsiae, cannot be used for typhus group rickettsiae because it only occurs in the group as a remnant gene (25). Using DNA microarray, Ge et al. found that the genetic variation between the Breinl and Madrid E strains was only 3% (20), thus confirming the unsuitability of coding DNA for typing rickettsiae at the strain level. Recently, we developed a genotyping method for Yersinia pestis named multispacer typing (MST) (12). This method, based on the comparison of several intergenic spacer sequences, has proven useful among strains of Y. pestis (12), R. conorii (19), Bartonella quintana (16), and Coxiella burnetii (22).
In this report, we applied MST based on 25 intergenic spacers to a total of 15 R. prowazekii strains or DNA amplicons of the same organism made from human body lice. By comparing these sequences, we were able to estimate the usefulness of MST for genotyping R. prowazekii at the strain level.
|
|
|---|
Rickettsia prowazekii strains and amplicons. In our study, we used cultures of four R. prowazekii strains (Table 1) and compared their sequences to those obtained from the DNA, extracted as described below, from the Madrid E strain and provided by Patrick Rozmajzl (Naval Medical Research Center, Silver Spring, MD). Rickettsial strains were propagated onto L929 cell monolayers (ATCC CCL NCTC clone 929) at 35°C in Eagle's minimal essential medium (Seromed, Berlin, Germany) supplemented with 4% fetal bovine serum (Seromed) and 2 mM glutamine. When Gimenez-stained cells (21) were heavily infected (after 3 to 5 days of culture), they were harvested, centrifuged (12,000 x g for 10 min), resuspended in minimal essential medium, and stored at 70°C until processed further.
|
View this table: [in a new window] |
TABLE 1. Rickettsia prowazekii strains and louse amplicons used in our study
|
Selection of target sequences. We aligned the genome sequences of R. conorii (GenBank accession number NC_003103) and R. prowazekii (accession number NC_000963) using BLAST software (1) and identified conserved or degraded fragments within intergenic sequences. The 25 most variable intergenic spacers were selected from 100- to 500-bp sequences separating two consecutive genes in both genomes which had a BLASTN value of <75 between both genomes.
DNA extraction and PCR-based sequencing method. Genomic DNA was extracted from rickettsial cultures using the QIAamp Tissue kit (QIAGEN, Hilden, Germany) according to the manufacturer's instructions.
The PCR primers used in this study were obtained from Eurogentec (Seraing, Belgium) and are described in Table 2. Their specificity was verified using BLAST software (1). PCRs were carried out in a PTC-200 automated thermal cycler (MJ Research, Waltham, MA). Two microliters of the DNA preparation was amplified in a 50-µl reaction mixture containing 50 pM of each primer; 200 µM (each) dATP, dCTP, dGTP, and dTTP (Invitrogen); 1 U eLONGase polymerase (Invitrogen, Gaithersburg, MD), 2 µl of eLONGase buffer A, and 8 µl of eLONGase buffer B. The following conditions were used for amplification: an initial 3 min of denaturation at 94°C was followed by 40 cycles of denaturation for 30 s at 94°C, annealing for 30 s at various temperatures indicated in Table 1, and extension for 1 min at 68°C. Amplification was completed by holding the reaction mixture for 3 min at 68°C to allow complete extension of the PCR products. PCR products were purified using a QIAquick Spin PCR purification kit (QIAGEN) as described by the manufacturer. Sequencing reactions were carried out using the d-Rhodamine Terminator cycle sequencing ready reaction kit with Amplitaq Polymerase FS (Perkin-Elmer, Coignieres, France) as described by the manufacturer. For all PCR products, sequences from both DNA strands were determined twice. Sequencing products were resolved using an ABI 3100 automated sequencer (Perkin-Elmer). Sequence analysis was performed using the software package ABI Prism DNA Sequencing Analysis Software version 3.0 (Perkin-Elmer). Sterile water was used as a negative control in each assay.
|
View this table: [in a new window] |
TABLE 2. Sequences of primers, amplicon sizes, and annealing temperatures used in this studyb
|
Nucleotide sequence accession numbers. Sequences determined in this study were deposited in GenBank under the following accession numbers: AY695448, AY695447, and AY695449 for the rpmE/tRNAfMet genotypes A, B, and C, respectively; and AY695454 and AY695452 for the serS/virB4 genotypes A and B, respectively.
|
|
|---|
|
View this table: [in a new window] |
TABLE 3. Classification of R. prowazekii and louse amplicons with MST genotyping
|
When the sequences from both spacers from R. prowazekii were compared to those of other rickettsial genomes, the degree of nucleotide sequence similarity ranged from 89.6% by comparison with R. akari to 91.4% with R. typhi for the rpmE/tRNAfMet spacer, compared with 99.6 to 100% among R. prowazekii strains. For the serS/virB4 spacer, the degree of nucleotide sequence similarity ranged from 88.1% by comparison with R. typhi to 92.4% with R. akari, compared with 99.6 to 100% among R. prowazekii strains.
Multispacer typing. Combining the results obtained from the analysis of the two variable spacers enabled us to identify four genotypes among the R. prowazekii strains and louse amplicons we studied (Table 3). Genotype 1 contained the Breinl strain and the Rw26862, Rw26877, Bur12727, Bur12728, and Bur12729 louse amplicons; genotype 2 was made of strains Madrid E, Evir, BatnaRp22, and Kuzina and of the Rw26875, Rw26879, and Bur12726 louse amplicons; genotype 3 contained only the Rw26860 louse amplicon; and genotype 4 included only the Bur12749 louse amplicon.
When sequences from the rpmE/tRNAfMet and serS/virB4 spacers were concatenated, the dendrograms obtained with the three different tree-building analysis methods used showed similar organization. All R. prowazekii strains and amplicons were classified within two clusters. The first cluster contained members of genotypes 1, 3, and 4, with representatives from genotypes 3 and 4 being issued from genotype 1. Representatives of genotype 2 clustered into a second group. R. typhi was positioned as the outermost taxon of studied species and strains (Fig. 1).
![]() View larger version (12K): [in a new window] |
FIG. 1. Phylogenetic relationships between 15 R. prowazekii strains and louse amplicons inferred by comparison of the sequences of the rpmE/tRNAfMet and serS/virB4 intergenic spacers using the neighbor-joining method. The scale bar represents a 2.5% nucleotide sequence divergence. Bootstrap values are indicated at the nodes.
|
|
|
|---|
The bioterrorism risk has recently reemerged, which highlighted the need for unambiguous, discriminatory strain characterization schemes. Prior to our study, there was no genotyping method described for R. prowazekii at the strain level. The development of a typing method for rickettsial strains has become crucial with the classification of R. prowazekii as a potential agent of bioterrorism. The recent finding by Ge et al. that limited variations occurred between coding sequences of two R. prowazekii strains (20) confirmed that there was little intraspecies variability in various gene sequences (18, 32, 33, 35). However, about 24% of the R. prowazekii genome consists of noncoding DNA, a very high percentage compared to that in other microbes sequenced so far (25). It has been suggested that intergenic spacer sequences are an important source of bacterial genome variability because they do not undergo selection pressure (11). For rickettsiae, it has been suggested that most of the intergenic sequences of R. prowazekii and R. conorii consist of decayed genes that are no longer active but have not yet been totally eliminated from the genome (2, 25). Recently, we demonstrated that the MST method, based on the comparison of intergenic spacer sequences, was a rational method for genotyping Y. pestis (12), R. conorii (19), B. quintana (16), and C. burnetii (22) at the strain level. In the present study, the comparison of the R. conorii and R. prowazekii genomes, which exhibit a high degree of colinearity, enabled us to select the 25 intergenic sequences exhibiting the highest interspecies variability. Two of these also had interstrain variability in the R. prowazekii strains and louse amplicons we studied. Sequences from the two variable spacers from R. prowazekii were clearly different from, and thus could not be confounded with, those of other Rickettsia species. When variations in these two variable spacers were considered, we could identify four genotypes among 15 R. prowazekii strains and louse amplicons. Although the number of strains or louse amplicons we studied may seem small, these include strains from three of the four current endemic foci of louse-borne typhus. R. prowazekii strains were classified within two genotypes, with genotype 1 incorporating the Breinl strain isolated during a large typhus outbreak in Poland following World War I (38) and genotype 2, including three European strains and one Algerian strain, isolated during World War II and later (4, 7, 9, 13). The body louse amplicons we studied were classified within four genotypes. The amplicons of the lice from Burundi could be separated into two genotypes, as could those made from the lice from Rwanda. As these lice were collected in a refugee camp in Burundi and in a jail in Rwanda during the same period in 2001, our data show that more than one strain can circulate in lice, and thus, louse-infested populations may be threatened by several strains of R. prowazekii at one time. These findings support our strategy of selecting intergenic sequences with the greatest interspecies variability as targets for strain typing and emphasize the usefulness of a typing technique as sensitive as ours to trace epidemic strains and differentiate natural from intentional outbreaks. The phylogenic organization of studied strains matched their genotypic classification (Fig. 1). Genotypes 3 and 4 were likely to be issued from the same lineage as members of genotype 1.
Initially, rather than amplifying and sequencing intergenic spacers from strain Madrid E, we used spacer sequences available in the genome sequence from R. prowazekii strain Madrid E (GenBank accession number NC_000963). We were surprised to find nucleotide differences between the genome sequence and strain Evir at seven positions within three spacers (see Results). These differences classified Madrid E within a unique MST type (data not shown). Strains Madrid E and Evir were described as having a common origin (4). Madrid E, an attenuated mutant of R. prowazekii obtained by passage through yolk sacs, was first described in 1943 (10) and proposed as a vaccine (26). In 1972, Balayeva and Nikolskaya reported that the virulence of the Madrid E strain was enhanced following passage through mice or guinea pigs (4, 5), and the revertant, virulent strain was named Evir. However, Wisseman suspected that a laboratory contamination of the Madrid E strain by a virulent R. prowazekii strain was the source of the emergence of the Evir strain rather than a spontaneous reversion to a virulent state (37). As the instability of attenuation of Madrid E raised questions about its suitability as a live vaccine antigen against epidemic typhus (30), it is an important issue to determine whether Evir and Madrid E have a common origin. Because of this polemic, the sequence differences we observed, and the fact that we have previously demonstrated that cell culture passages do not alter MST genotypes in R. conorii (19), we considered it possible that Madrid E and Evir were issued from different lineages. However, an anonymous reviewer suggested that these sequence differences could be explained by errors in the genome sequence. In consequence, we asked Patrick Rozmajzl to provide us with DNA from R. prowazekii strain Madrid E and determined the sequences of the three spacers discordant between Evir and the Madrid E genome sequence. The spacer sequences we obtained from strains Madrid E and Evir were carefully checked nucleotide by nucleotide in both directions and were found to be identical between both strains. Therefore, we are confident that Madrid E and Evir belong to the same MST genotype and speculate that the mutations observed in the R. prowazekii genome are the result of sequencing errors.
The combined use of variable spacer sequences, which we have named multispacer typing (MST), could identify four genotypes among 15 R. prowazekii strains and body louse amplicons. This technique, which is easy to use, may be applied for tracking strains and may even be applied directly to clinical specimens or body lice, which have been demonstrated to be useful field specimens for diagnosis and surveillance of epidemic typhus outbreaks (17). Our use of MST enabled us to show that several strains of R. prowazekii were involved in the outbreak of epidemic typhus in Africa in 2001.
No exterior funding was obtained for this research.
The authors of this research do not have any conflict of interest.
In memory of Natasha Balayeva, a famous rickettsiologist and friend. ![]()
Y.Z. and P.-E.F. 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»