Microbiology
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


Microbiology 153 (2007), 160-168; DOI  10.1099/mic.0.2006/001149-0
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Supplementary Table
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via CrossRef
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Vitorino, L.
Right arrow Articles by Zé-Zé, L.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Vitorino, L.
Right arrow Articles by Zé-Zé, L.
Agricola
Right arrow Articles by Vitorino, L.
Right arrow Articles by Zé-Zé, L.
Microbiology 153 (2007), 160-168; DOI  10.1099/mic.0.2006/001149-0
© 2007 Society for General Microbiology

Rickettsiae phylogeny: a multigenic approach

Liliana Vitorino1, Ivo M. Chelo1, Fatima Bacellar2 and Líbia Zé-Zé1,{dagger}

1 Universidade de Lisboa, Faculdade de Ciências. Centro de Genética e Biologia Molecular and Instituto de Ciência Aplicada e Tecnologia. Edifício ICAT, Campus FCUL, Campo Grande, 1749-016 Lisboa, Portugal
2 CEVDI, Instituto Nacional de Saúde Dr Ricardo Jorge, INSA, Av. Padre Cruz, 1649-016 Lisboa, Portugal

Correspondence
Líbia Zé-Zé
libia.zeze{at}insa.min-saude.pt


    ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS AND DISCUSSION
 REFERENCES
 
The development of molecular taxonomic methods has provided a large amount of data in the reorganization of Rickettsiae taxonomy. Nevertheless, phylogenetic relationships among some groups and species delimitation remain unclear. To clarify rickettsial phylogeny, a multigenic approach was used for the first time for the genus Rickettsia, based on simultaneous analyses of eight loci: atpA, recA, virB4, dnaA, dnaK, rrl-rrf internal transcribed spacer, ompA and gltA. Concatenation of different nucleotide sequences resulted in an improvement in phylogenetic resolution when compared to single gene data. This multigenic approach has enabled the differentiation of many groups, including the spotted fever group which includes a great number of closely related species. The reliability of some previously recognized groups was evaluated.


Abbreviations: dN/dS, non-synonymous/synonymous substitutions distances ratio; ITS, internal transcribed spacer; ML, maximum-likelihood; MLST, multilocus sequence typing; MP, maximum-parsimony; NJ, neighbour-joining; SFG, spotted fever group; TG, typhus group

A table of GenBank/EMBL/DDBJ accession numbers for the sequences reported in this paper is available as supplementary data with the online version of this paper.

{dagger}Present address: CEVDI, Instituto Nacional de Saúde Dr Ricardo Jorge, INSA, Av. Padre Cruz, 1649-016 Lisboa, Portugal.


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS AND DISCUSSION
 REFERENCES
 
The genus Rickettsia belongs to the Alphaproteobacteria and represents a group of obligate intracellular endosymbionts of eukaryotic cells (Weiss & Dasch, 1991Down). The genus subdivides into the conventionally well defined typhus group (TG) and spotted fever group (SFG), based mainly on phenotypic and serological features. The genus Rickettsia also includes several marginal species, such as Rickettsia canadensis and Rickettsia bellii, the most divergent species within the Rickettsiae (Stothard et al., 1994Down). The recently sequenced R. bellii genome was shown to be larger than the other rickettsial genomes sequenced so far and exhibits little colinearity with any of them, suggesting that this species may have retained several ancestral features lost in other lineages in the course of evolution (Ogata et al., 2006Down).

The TG includes Rickettsia prowazekii and Rickettsia typhi, which are the aetiological agents of insect-borne epidemic and endemic typhus, respectively. The SFG consists of several species, some of which are the causative agents of well known tick-borne diseases, such as Rocky Mountain spotted fever (Rickettsia rickettsii) and Mediterranean spotted fever (Rickettsia conorii) (Raoult & Roux, 1997Down).

Rickettsiae phylogeny has been addressed by sequence analyses of certain genes, mostly those belonging to the sca (surface cell antigen) family (Roux & Raoult, 1995Down, 2000Down; Roux et al., 1997Down; Fournier et al., 1998Down; Sekeyova et al., 2001Down; Ngwamidiba et al., 2006Down), enabling the division of this genus into six groups represented by the species: Rickettsia rickettsii, Rickettsia massiliae, Rickettsia helvetica, Rickettsia akari, Rickettsia prowazekii and Rickettsia canadensis. Nevertheless, a robust and accurate phylogeny of this genus has not yet been accomplished and ambiguities remain in the phylogenetic position of some species. Moreover, for the last two decades, more novel Rickettsia species have been described, contributing to an increase in the complexity of rickettsial taxonomy.

As recommend by Stackebrandt et al. (2002)Down, the use of several genes scattered across the bacterial genome boosts the informative level of phylogenetic data, improving the discriminatory power and the robustness of the phylogenetic relationships. Ideally, one should use data from housekeeping genes, which are ubiquitous within the population and tend to evolve at a moderate rate, and from hypervariable genes which might encode outer-surface proteins and typically evolve at an accelerated rate (Cooper & Feil, 2004Down). Hence, we have analysed eight target nucleotide sequences, atpA, recA, virB4, dnaA, dnaK, rrl-rrf internal transcribed spacer (ITS), ompA and gltA, selected as representative of different genomic functions believed to be active in most rickettsial species. In the present study, we compared the phylogenetic trees derived from different datasets and analysed for the first time the discriminatory capacity of concatenated sequences to infer stable phylogenetic relationships within rickettsial species. With this multigenic approach we have been able to achieve a good overall phylogenetic resolution at different levels of genetic resemblance.


    METHODS
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS AND DISCUSSION
 REFERENCES
 
Bacterial strains, culture conditions and DNA extraction.
Strains used are listed in Table 1Down. The rickettsiae were propagated on Vero cell (E-6 clone) monolayers in Eagle's minimal medium (MEM; Gibco-BRL) supplemented with 1 % (v/v) glutamine and 8 % (v/v) fetal bovine serum at 32 °C for 5–6 days. Infected cell cultures were monitored for the presence of rickettsiae. Harvesting was done when the degree of infection was optimal, as determined by Gimenez staining (Gimenez, 1964Down). The cultures were frozen and unfrozen three times, and then centrifuged at 121 g for 30 min. The supernatants were recovered, then the rickettsiae were harvested by 30 min centrifugation at 7740 g and resuspended in ultra-pure water. DNA was extracted using Qiagen columns according to the DNAeasy protocol for the Animal Tissue Kit.


View this table:
[in this window]
[in a new window]

 
Table 1. Rickettsial strains used in this study

 
PCR amplification and gene sequencing.
Primers used to amplify and sequence the several gene targets are listed in Table 2Down. PCR was performed in a 50 µl reaction mixture containing 1 pmol each primer, 200 µM each dATP, dGTP, dCTP and dTTP (Invitrogen), 1.75 U Taq polymerase (Invitrogen), 2 mM MgCl2, 0.5x BSA, 1x Taq buffer and 50–100 ng genomic DNA. Amplification was carried out in a DNA thermocycler (TGradient; Biometra) using the following conditions: 5 min of initial denaturation at 96 °C, then 35 cycles at 94 °C for 1 min, 52 °C for 1 min and 72 °C for 1 min. The amplification was completed by holding for 5 min at 72 °C to allow complete extension of the PCR product. When no amplification product was obtained, the annealing temperature was lowered to 50 °C and the MgCl2 concentration was raised to 2.5 mM. PCR products were visualized by ethidium bromide staining after electrophoresis in a 1 % (w/v) agarose gel and their sizes were estimated by comparison with a molecular mass standard (1 kb plus DNA ladder; Gibco-BRL). The PCR products were purified using the Jet Quick-PCR Purification Kit (Genomed) as described by the manufacturer. The purified PCR products were sequenced in both directions using an automated DNA capillary sequencer CEQ 2000-XL (Beckman Coulter) and a dye-labelled dideoxy termination method (DTCS, Dye Terminator Cycle Sequencer Start Kit; Beckman Coulter). Sequences were proofread, edited and merged into composite sequences using CEQ 8000 software (Beckman Coulter). Accession numbers of nucleotide sequences are listed in supplementary Table S1 (available with the online version of this paper). All sequence data were posted at the multilocus sequence typing (MLST) website, http://rickettsia.mlst.net (Aanensen & Spratt, 2005Down).


View this table:
[in this window]
[in a new window]

 
Table 2. Primers used for PCR amplification and sequencing

 
Data analysis.
For phylogenetic inference, nine different alignments were created, eight corresponding to the alignment of sequences of each gene: atpA, recA, dnaA, rrl-rrf ITS, virB4, gltA, dnaK and ompA and one alignment corresponding to the concatenation of all gene sequences. All alignments were made using the CLUSTAL W program (Thompson et al., 1997Down). For protein-encoding genes, DNA alignments were made based on the corresponding amino acid alignments by using CodonAlign 2.0 (Hall, 2004Down). Aligned regions with gaps were removed from the alignment before phylogenetic analysis.

Phylogenetic analysis including neighbour-joining (NJ), maximum-parsimony (MP) and maximum-likelihood (ML) trees of DNA sequence alignments were conducted using PAUP* 4.0b10 software (Swofford, 2000Down). Genetic distances used in NJ trees are Kimura two-parameter model distances with a transition : transversion ratio of 2 : 1. Non-synonymous/synonymous substitutions distances (dN/dS) ratios, as described by Nei & Gojobori (1986)Down, were calculated using MEGA3 software (Kumar et al., 2004Down). Analysis of protein sequences were done using several programs present in the PHYLIP software package, version 3.6 (Felsenstein, 1989Down). Bootstrap analysis was made with 1000 replicates except in ML where only 100 replicates were generated. For the Bayesian approach, MrBayes software version 3.0b4 (Huelsenbeck & Ronquist, 2001Down) was used. Each analysis consisted of 2.0x107 generations from a random starting tree and four Markov chains (with default heating values) sampled every 100 generations. The first 8000 sampled trees were discarded resulting in a set of 12 000 analysed trees sampled after stationarity. To prevent reaching only apparent stationarity, two separate runs were made for each analysis. Hierarchical likelihood ratio tests were conducted using a batch file supplied with MODELTEST 3.06 (Posada & Crandall, 1998Down) to provide the evolutionary models used in ML and Bayesian analysis. Based on the work by Stothard et al. (1994)Down, R. bellii was used as outgroup.


    RESULTS AND DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS AND DISCUSSION
 REFERENCES
 
Phylogenetic relationships inferred from throughout the genome are viewed as accurate estimates of those that would be obtained from an analysis of the entire genome, since the effects of location-dependent processes in sequence evolution can be reduced by sampling sites from different regions (Cummings et al., 1995Down). Hence, several target nucleotide sequences were selected as representative of different genomic functions (Table 2Up) believed to have an important role in the rickettsial life cycle, and with regard to their scattered distribution in the rickettsial chromosome (Andersson et al., 1998Down; Ogata et al., 2001Down). With this approach we expected to use sequences that were not contiguous in the rickettsial genome to minimize the effect of location-dependent nucleotide evolution (Cummings et al., 1995Down). Moreover, we used both housekeeping and hypervariable genes to provide a framework for studying closely related strains as well as different species (Cooper & Feil, 2004Down). Although our main goal was not rickettsial strain characterization, considering that this data can be used in an MLST scheme, all sequence data were posted to an MLST database as a focal point for future storage and retrieval of rickettsial sequence, strain and host data.

Thirty-four rickettsial strains were analysed. For R. bellii and R. canadensis the PCR amplification of some of the genes was difficult or even impossible (supplementary Table S1, available with the online version of this paper). This was not unexpected, since a greater genomic divergence is recognized for these two species, preventing the amplification of some genes even when carried out with consensus primers (designed from known Rickettsia spp. sequences), as described by Sekeyova et al. (2001)Down.

Due to the genetic closeness of rickettsial species, nucleotide sequences were used for phylogenetic analysis instead of protein sequences. Phylogenetic trees based on amino acid sequences were drawn at a preliminary stage (results not shown), but their lack of resolution was obvious even for the most divergent gene (ompA).

Given that the aim of this study was to conduct a wide-ranging clarification of the phylogenetic relationships within the genus Rickettsia by the use of different genes, we decided to make a joint analysis based on the assembly of all the sequences, to guarantee improvement of phylogenetic tree robustness. Incongruence length difference tests have been used for data concatenation; these tests can be good indicators of different evolutionary histories and are used to prevent the combination of incongruent data (Cunningham, 1997Down; Bull et al., 2003Down). To assess the adequacy of the all-gene tree to explain the evolution of each locus, we also performed Shimodaira–Hasegawa tests with re-sampling estimated log likelihood optimization (Goldman et al., 2000Down). In this test, different tree topologies were compared and their adequacy was tested based upon the comparison of their log likelihood values. Usually a P value <=0.05 is used for the rejection of a given tree topology representing a hypothesis that might explain the data. Beyond the topology given by the all-data tree, the set of tested trees included the ML tree of each locus. Shimodaira–Hasegawa tests indicated that the tree constructed from concatenated data provides a plausible explanation for the evolutionary history of every gene (data not shown). Furthermore, single-gene trees were also obtained, enabling comparison of both approaches.

The analysis of the concatenated gene set (Fig. 1Down) was made using a Bayesian approach where cluster confidence is given by posterior probabilities (Huelsenbeck et al., 2001Down).


Figure 1
View larger version (20K):
[in this window]
[in a new window]

 
Fig. 1. MrBayes tree using concatenated sequences: atpA, recA, dnaA, rrl-rrf ITS, virB4, gltA, dnaK and ompA. The radial tree layout is shown to give a better visualization of the overall Rickettsiae phylogeny, and the phylogram is shown to highlight the lineage diversity within the SFG (a, in grey). In the phylogram posterior probabilities and ML bootstrap values are indicated to provide branch support. Bar, 10 % sequence divergence.

 
To account for different gene-specific evolutionary rates and rate patterns, the data were partitioned for the several genes, and an appropriate evolutionary model was used for each gene (Castoe et al., 2004Down). In addition, the confidence of the groups obtained by this approach was also evaluated by bootstrap values taken from ML trees. However, in the latter case no partition of the data was made and a single model of DNA evolution was computed for the entire set. While the chosen model for the ITS region and dnaK was HKY (Hasegawa et al., 1985Down) with gamma variation, for all the others the GTR model (Lanave et al., 1984Down), also with gamma variation, was selected. In atpA, a proportion of invariable sites were considered. In the ML analysis, the GTR+G+I model was used for the entire set.

The analysis of the all-gene phylogenetic tree (Fig. 1Up) allows the recognition of several rickettsial groups whose phylogenetic clustering seems to be robust, with ML bootstrap values above 70 % and posterior probability values close to 1.00. Most of these groups were also recovered in MP and NJ trees (results not shown), indicating that the increase in tree resolution can be attributed to an increase in phylogenetic signal and not to the reconstruction method used. Moreover, these results are compatible with single-gene phylogenies (Fig. 2Down), since there are no mutually exclusive groups with significant bootstrap support. In Fig. 2Down, NJ trees of single-gene phylogenies are shown. Bootstrap proportions are given only for the NJ method, but branches leading to groups also supported by MP and ML trees, with bootstrap values >70 %, are highlighted.


Figure 2
Figure 2
View larger version (46K):
[in this window]
[in a new window]

 
Fig. 2. NJ tree obtained from (a) atpA, (b) recA, (c) dnaA, (d) rrl-rrf ITS, (e) virB4, (f) gltA, (g) dnaK and (h) ompA sequences, and Kimura two-parameter model distances. Bootstrap proportions are given only for the NJ method. Only bootstrap values >70 % are shown. Branches leading to groups also supported by MP and ML trees are highlighted.

 
The phylogeny obtained with all the data also has the advantage of containing information from genes with different rates of evolution that are informative at different levels. This is particularly important since the phylogenetic inferences to be made are aimed both at the comparison of genetically close taxa of the traditional SFG and of the Rickettsia genus as a whole. ompA and the ITS region, which in the SFG have a mean genetic distance (Kimura two-parameter distance with a transition : transversion ratio 2 : 1) of 0.057 nucleotide substitutions per site, are suitable for the clarification of phylogenetic relationships within the SFG. In contrast, phylogenies obtained from atpA, recA, dnaA, virB, gltA and dnaK, in which the mean intra-SFG group distances vary between 0.018 and 0.027, are more suited to the establishment of more basal relationships, providing good delineation of the main groups of the Rickettsiae. As an example, one can see that phylogenies obtained from recA, virB and dnaK (Fig. 2Up) give strong support to most of the groups suggested by Sekeyova et al. (2001)Down and reviewed by Parola et al. (2005)Down.

It is generally recognized that the phylogenetic relationships of the Rickettsiae are characterized by the presence of several groups with one or a few identifiable individuals, and that one group, the SFG, contains most of the recognized species. It is interesting to see that all of the groups recognized so far (Sekeyova et al., 2001Down) join each other in the middle of the phylogenetic tree, being connected by relatively short internal branches. Not surprisingly, the confidence of the clustering of any of these groups is hampered by low support in either bootstrap or posterior probability values. It should also be noted that these results do not support a monophyletic nature for the SFG. Both concatenated and single-gene sequence data indicate instead that SFG is divided into two groups: a monophyletic cluster, which includes the R. rickettsii group, the R. massiliae group (as defined by Sekeyova et al., 2001Down) and Rickettsia montana, and a second group, consisting of Rickettsia australis, Rickettsia akari and Rickettsia felis. The only exception to this evolutionary history is presented by the virB4 phylogeny, in which the SFG is seen as monophyletic, but without support. Regarding the R. akari group, only R. australis and R. akari form a well supported monophyletic cluster. This is to be expected, taking into consideration the hypothetical early separation of R. felis (as suggested in virB4, ITS and gltA phylogenies; Fig. 2d, e and gUp). As can be clearly seen from the radial tree in Fig. 1Up, R. felis branches from a central position in Rickettsia phylogeny. This could explain the variable phylogenetic positions or support values for R. felis in our study and also in other studies using different genes (Fournier et al., 2003Down; Zhu et al., 2005Down; Ngwamidiba et al., 2006Down). This relatedness between both the SFG and TG rickettsiae could account for the presence in the R. felis genome of genes exclusively shared with R. typhi; this may be the reason for their serological cross-reactivity, as pointed out by Ogata et al. (2005)Down.

The previously mentioned monophyletic subgroup of the SFG rickettsiae is distinguished by the presence of many genetically close isolates. An overall good phylogenetic resolution is nevertheless possible, with many increasingly specific subgroups being recognized. In the R. rickettsii group, as mentioned by Sekeyova et al. (2001)Down the R. conorii complex is well established and can be readily subdivided into a group containing R. conorii Malish-like isolates, a group of Israeli-tick-typhus-like isolates and Rickettsia sp. A-167, the latter two being most closely related. This is in agreement with the resolution of this group as presented by Zhu et al. (2005)Down. Although in this study, only Portuguese Israeli tick typhus isolates were used, ompA phylogenies including Israeli tick typhus strain ISTT (accession no. U83441) resulted in the same clear separation within the R. conorii complex (results not shown).

Another well defined monophyletic group includes not only Rickettsia sibirica strains, but also Rickettsia africae. Although the separation of these two can be made based on the all-gene phylogeny (most probably as a result of the information present in the ompA gene; see Fig. 2fUp), the proximity of R. africae to the R. sibirica group is also clear and becomes almost indistinguishable in recA phylogeny (Fig. 2bUp). The monophyly of this clade has already been recognized and it also contains Rickettsia parkeri and Rickettsia sp. strain S (Roux et al., 1997Down; Fournier et al., 1998Down; Roux & Raoult, 2000Down; Sekeyova et al., 2001Down), which were not included in this study. Since evidence points towards it being a sister taxon to R. conorii, the attribution of a taxonomically defined rank such as species should be subject to further discussion.

In the remainder of the R. rickettsii group there is a distinct cluster with several Rickettsia slovaca strains (including two Portuguese isolates) and, of special notice, is the consistent position of Rickettsia japonica as the most divergent line in the group.

As found in the analysis of Sekeyova et al. (2001)Down, the R. massiliae group, consisting of R. massiliae, Rickettsia sp. Bar29, Rickettsia rhipicephali and Rickettsia aeschlimannii is well supported, although only in the ompA single-gene tree (Fig. 2fUp) is it recovered with high bootstrap values by all methods. As shown in Fig. 1Up the phylogenetic relationship between the R. massiliae group, the R. rickettsii group and Rickettsia montanensis is uncertain in this analysis. This fact, which can be seen as a difficulty in the resolution of the basal lines with regard to this diversification, gives credit to the previously recognized (Sekeyova et al., 2001Down) variable position of R. montanensis, thus care should be taken when considering its inclusion in the R. massiliae group, as suggested recently (Parola et al., 2005Down).

We should also draw attention to the phylogenetic position of R. canadensis, which is one of the most variable. In the dnaA gene phylogeny R. canadensis is even hypothesized as a member of the SFG, with high bootstrap proportions and posterior probability values. Nevertheless, this phylogeny is unique because of the long phylogenetic line leading to R. canadensis, most probably a result of an unusually high genetic divergence of this species. Indeed, R. canadensis was initially considered as a member of the TG on the basis of antigenic similarities; however, phylogenetic analysis (Roux et al., 1997Down; Fournier et al., 1998Down) has confirmed that this species constitutes a separate group from the SFG and TG rickettsiae. Remarkably, in spite of its unclear phylogenetic position, this species seems to be closer to R. helvetica than to members of the TG in the concatenated data analysis (Fig. 1Up).

As in the TG group, R. bellii is phylogenetically distant from all the other members of the Rickettsiae, which is in agreement with its distinct genomic (Ogata et al., 2006Down) and phenotypic (Philip et al., 1983Down) characteristics. Indeed, R. bellii is believed to be the only species of the first divergent line of the genus Rickettsia (Stothard et al., 1994Down). Despite its unquestionable use as an outgroup for other rickettsial species, it is impossible to recognize the order of divergence of the major groups due to the short length of internal branches, as already mentioned. This issue might be overemphasized by the missing data not only for R. bellii but also for R. canadensis. In fact, even if four gene sequences were obtained for these two species, only two of them would be common to both. This could seriously undermine inferences resulting from the presence of data from both species that would contribute to the clarification of the relatedness of the major rickettsial groups.

In conclusion, the combined use of the sequences of several housekeeping genes has led to the construction of a multi-gene phylogeny, which helps in clarifying the taxonomic structure of some groups of interest within the genus Rickettsia, particularly within the SFG. This group is not monophyletic, highlighting the limitations of serotyping in the delineation of taxonomic groups. Our study supports the distinctiveness of R. montanensis, which should thus be considered separate from the R. massiliae and R. rickettsii groups. In this analysis we have used several Portuguese isolates, which enabled us to take intraspecies variability into account. However, little resolution could be achieved at this level of relatedness. Obviously the multigenic approach to Rickettsiae phylogeny would be reinforced with the addition of more rickettsial species and new isolates that are constantly being described.


    ACKNOWLEDGEMENTS
 
The ICAT/FCUL research team involved in this work strongly acknowledges its scientific leader, Professor Rogério Tenreiro, for general scientific supervision and engagement in providing the necessary resources and facilities. We also thank Edward J. Feil for helpful discussion of the final version of this manuscript. This publication made use of the Multi Locus Sequence Typing website at Imperial College London developed by David Aanensen and Man-Suen Chan, and funded by the Wellcome Trust. L. Vitorino and I. M. Chelo are recipients of FCT research grants SFRH/BD/10676/2002 and SFRH/BD/10675/2002, respectively.

Edited by: J. Parkhill


    REFERENCES
 TOP
 ABSTRACT
 INTRODUCTION
 METHODS
 RESULTS AND DISCUSSION
 REFERENCES
 
Aanensen, D. M. & Spratt, B. G. (2005). The multilocus sequence typing network: mlst.net. Nucleic Acids Res 33 (Web Server issue), W728–W733.[Abstract/Free Full Text]

Andersson, S. G. E., Zomorodipour, A., Andersson, J. O., Sicheritz-Pontén, T., Alsmark, U. C. M., Podowski, R. M., Näslund, A. K., Eriksson, A.-S., Winkler, H. H. & Kurland, C. G. (1998). The genome sequence of Rickettsia prowazekii and the origin of mitochondria. Nature 396, 133–140.[CrossRef][Medline]

Bull, J. J., Huelsenbeck, J. P., Cunningham, C. W., Swofford, D. L. & Waddell, P. J. (2003). Partitioning and combining data in phylogenetic analysis. Syst Biol 42, 384–397.[CrossRef]

Castoe, T. A., Doan, T. M. & Parkinson, C. L. (2004). Data partitions and complex models in Bayesian analysis: the phylogeny of gymnophthalmid lizards. Syst Biol 53, 448–469.[CrossRef][Medline]

Cooper, J. E. & Feil, E. J. (2004). Multilocus sequence typing – what is resolved? Trends Microbiol 12, 373–377.[CrossRef][Medline]

Cummings, M. P., Otto, S. P. & Wakeley, J. (1995). Sampling properties of DNA sequence data in phylogenetic analysis. Mol Biol Evol 12, 814–822.[Abstract]

Cunningham, C. W. (1997). Can three incongruence tests predict when data should be combined? Mol Biol Evol 14, 733–740.[Abstract]

Felsenstein, J. (1989). PHYLIP – Phylogeny Inference Package (Version 3.2). Cladistics 5, 164–166.

Fournier, P. E., Roux, V. & Raoult, D. (1998). Phylogenetic analysis of the spotted fever Rickettsiae by study of the outer surface protein rOmpA. Int J Syst Bacteriol 48, 839–849.[Abstract/Free Full Text]

Fournier, P. E., Dumler, J. S., Greub, G., Zhang, J., Wu, Y. & Raoult, D. (2003). Gene sequence-based criteria for identification of new Rickettsia isolates and description of Rickettsia heilongjiangensis sp. nov. J Clin Microb 41, 5456–5465.[Abstract/Free Full Text]

Gimenez, D. F. (1964). Staining Rickettsiae in yolk-sac cultures. Stain Technol 39, 135–140.[Medline]

Goldman, N., Anderson, J. P. & Rodrigo, A. G. (2000). Likelihood-based tests of topologies in phylogenetics. Syst Biol 49, 652–670.[CrossRef][Medline]

Hall, B. G. (2004). Phylogenetic Trees Made Easy: a How-To Manual for Molecular Biologists, 2nd edn. Sunderland, MA: Sinauer Associates.

Hasegawa, M., Kishino, H. & Yano, T. (1985). Dating the human–ape splitting by a molecular clock of mitochondrial DNA. J Mol Evol 22, 160–174.[CrossRef][Medline]

Huelsenbeck, J. P. & Ronquist, F. R. (2001). MrBayes: Bayesian inference of phylogenetic trees. Bioinformatics 17, 754–755.[Abstract/Free Full Text]

Huelsenbeck, J. P., Ronquist, F., Nielsen, R. & Bollback, J. P. (2001). Bayesian inference of phylogeny and its impact on evolutionary biology. Science 294, 2310–2314.[CrossRef][Medline]

Kumar, S., Tamura, K. & Nei, M. (2004). MEGA3: integrated software for Molecular Evolutionary Genetic Analysis and sequence alignment. Brief Bioinform 5, 150–163.[Abstract/Free Full Text]

Lanave, C., Preparata, G., Saccone, C. & Serio, G. (1984). A new method for calculating evolutionary substitution rates. J Mol Evol 20, 86–93.[CrossRef][Medline]

Nei, M. & Gojobori, T. (1986). Simple methods for estimating the number of synonymous and nonsynonymous nucleotide substitutions. Mol Biol Evol 3, 418–426.[Abstract]

Ngwamidiba, M., Blanc, G., Raoult, D. & Fournier, P.-E. (2006). ScaI, a previously undescribed paralog from autotransporter protein-encoding genes in Rickettsia species. BMC Microbiol 20, 6–12.[CrossRef]

Ogata, H., Audic, S., Renesto-Audiffren, P., Fournier, P.-E., Barbe, V., Samson, D., Roux, V., Cossart, P., Weissenbach, J. & other authors (2001). Mechanisms of evolution in Rickettsia conorii and R. prowazekii. Science 293, 2093–2098.[Abstract/Free Full Text]

Ogata, H., Renesto, P., Audic, S., Robert, C., Blanc, G., Fournier, P. E., Parinello, H., Claverie, J. M. & Raoult, D. (2005). The genome sequence of Rickettsia felis identifies the first putative conjugative plasmid in an obligate intracellular parasite. PLoS Biol 3, e248.[CrossRef][Medline]

Ogata, H., La Scola, B., Audic, S., Renesto, P., Blanc, G., Robert, C., Fournier, P. E., Claverie, J. M. & Raoult, D. (2006). Genome sequence of Rickettsia bellii illuminates the role of Amoebae in gene exchanges between intracellular pathogens. PLoS Genetics 2, 733–744.

Philip, R. N., Casper, E. A., Anacker, R. L., Cory, J., Hayes, S. F., Burgdorfer, W. & Yunker, C. E. (1983). Rickettsia bellii sp. nov., a tick-borne rickettsia, widely distributed in the United States, that is distinct from the spotted fever and typhus biogroups. Int J Syst Bacteriol 33, 94–106.[Abstract/Free Full Text]

Parola, P., Paddock, C. D. & Raoult, D. (2005). Tick-borne rickettsioses around the world: emerging diseases challenging old concepts. Clin Microbiol Rev 18, 719–756.[Abstract/Free Full Text]

Posada, D. & Crandall, K. A. (1998). MODELTEST: testing the model of DNA substitution. Bioinformatics 14, 817–818.[Abstract/Free Full Text]

Raoult, D. & Roux, V. (1997). Rickettsioses as paradigms of new or emerging infectious diseases. Clin Microbiol Rev 10, 694–719.[Abstract]

Regnery, R. L., Spruill, C. L. & Plikaytis, B. D. (1991). Genotypic identification of Rickettsiae and estimation of interspecies sequence divergence for portions of two rickettsial genes. J Bacteriol 173, 1576–1589.[Abstract/Free Full Text]

Roux, V. & Raoult, D. (1995). Phylogenetic analysis of the genus Rickettsia by 16S rDNA sequencing. Res Microbiol 146, 385–396.[Medline]

Roux, V. & Raoult, D. (2000). Phylogenetic analysis of members of the genus Rickettsia using the gene encoding the outer-membrane protein rOmpB (ompB). Int J Syst Evol Microbiol 50, 1449–1455.[Abstract]

Roux, V., Rydkina, E., Eremeeva, M. & Raoult, D. (1997). Citrate synthase gene comparison, a new tool for phylogenetic analysis and its application for the Rickettsiae. Int J Syst Bacteriol 47, 252–262.[Abstract/Free Full Text]

Sekeyova, Z., Roux, V. & Raoult, D. (2001). Phylogeny of Rickettsia spp. inferred by comparing sequences of ‘gene D’. which encodes an intracytoplasmic protein. Int J Syst Evol Microbiol 51, 1353–1360.[Abstract]

Stackebrandt, E., Frederiksen, W., Garrity, G. M., Grimont, P. A. D., Kämpfer, P., Maiden, M. C. J., Nesme, X., Rosselló-Mora, R., Swings, J. & other authors (2002). Report of the ad hoc committee for the re-evaluation of the species definition in bacteriology. Int J Syst Evol Microbiol 52, 1043–1047.[Abstract]

Stothard, D. R., Clark, J. B. & Fuerst, P. A. (1994). Ancestral divergence of Rickettsia bellii from the spotted fever and typhus groups of Rickettsia and antiquity of the genus Rickettsia. Int J Syst Bacteriol 44, 798–804.[Abstract/Free Full Text]

Swofford, D. L. (2000). PAUP* – Phylogenetic Analysis Using Parsimony and other methods (software). Sunderland, MA: Sinauer Associates.

Thompson, J. D., Higgins, D. G. & Gibson, T. J. (1997). The CLUSTALX–Windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tool. Nucleic Acids Res 25, 4876–4882.[Abstract/Free Full Text]

Vitorino, L., Zé-Zé, L., Sousa, A., Bacellar, F. & Tenreiro, R. (2003). rRNA intergenic spacer regions for phylogenetic analysis of Rickettsia species. Ann N Y Acad Sci 990, 726–733.[Medline]

Weiss, E. & Dasch, G. A. (1991). Introduction to the Rickettsiales and other parasitic or mutualistic prokaryotes. In The Prokaryotes, 2nd edn, pp. 2402–2406. Edited by A. Balows, H. G. Trüper, M. Dworking, W. Harder & K.-H. Schleifer. New York: Springer.

Zhu, Y., Fournier, P.-E., Eremeeva, M. & Raoult, D. (2005). Proposal to create subspecies of Rickettsia conorii based on multi-locus sequence typing and an emended description of Rickettsia conorii. BMC Microbiol 5, 11.[CrossRef][Medline]

Received 11 August 2006; revised 3 October 2006; accepted 12 October 2006.


This article has been cited by other articles:


Home page
Int. J. Syst. Evol. Microbiol.Home page
A. Alexandre, M. Laranjo, J. P. W. Young, and S. Oliveira
dnaJ is a useful phylogenetic marker for alphaproteobacteria
Int J Syst Evol Microbiol, December 1, 2008; 58(12): 2839 - 2849.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Supplementary Table
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Right arrow Citation Map
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via CrossRef
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Vitorino, L.
Right arrow Articles by Zé-Zé, L.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Vitorino, L.
Right arrow Articles by Zé-Zé, L.
Agricola
Right arrow Articles by Vitorino, L.
Right arrow Articles by Zé-Zé, L.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
INT J SYST EVOL MICROBIOL MICROBIOLOGY J GEN VIROL
J MED MICROBIOL ALL SGM JOURNALS
Copyright © 2007 Society for General Microbiology.