# Genotypic Characterization ofRickettsia bellii Reveals Distinct Lineages in the United States and South America

Authors: Felipe S. Krawczak; Marcelo B. Labruna; Joy A. Hecht; Christopher D. Paddock; Sandor E. Karpathy
Journal: BioMed Research International
 (2018)

---

## Abstract

The bacteriumRickettsia bellii belongs to a basal group of rickettsiae that diverged prior to the pathogenic spotted fever group and typhus groupRickettsiaspecies. Despite a diverse representation ofR. bellii across more than 25 species of hard and soft ticks in the American continent, phylogeographical relationships among strains of this basal group-Rickettsia species are unknown; the work described here explores these relationships. DNA was extracted from 30R. belliitick isolates: 15 from the United States, 14 from Brazil, and 1 from Argentina. A total of 2,269 aligned nucleotide sites of 3 protein coding genes (gltA,atpA, andcoxA) and 2 intergenic regions (rpmE-tRNAfmet andRC1027-xthA2) were concatenated and subjected to phylogenetic analysis by Bayesian methods. Results showed a separation of almost all isolates between North and South Americas, suggesting that they have radiated within their respective continents. Phylogenetic positions of the 30 isolates could be a result of not only their geographical origin but also the tick hosts they have coevolved with. WhetherR. belliioriginated with ticks in North or South America remains obscure, as our analyses did not show evidence for greater genetic divergence ofR. belliiin either continent.


---

## Body

## 1. Introduction

Members of the genusRickettsia (Rickettsiales: Rickettsiaceae) are Gram-negative obligate intracellular bacteria, usually in association with arthropods, although a number of distinct genotypes have also been described from leeches, amoeba, ciliate, and hydra [1]. Phylogenetic analysis-based studies have classified theRickettsiaspecies into five major groups, namely, the spotted fever group (SFG), the typhus group (TG), the transitional group (TRG), theRickettsia canadensis group, and theRickettsia bellii group [1, 2]. This later group occupies a basal position in all major phylogenetic studies, which indicates early divergence within the genus [1, 3, 4].Rickettsia bellii was formally described in 1983 based on isolates obtained from multiple species of ixodids (hard ticks) and argasids (soft ticks) in the United States [5]. In the United States,R. belliihas been reported to be infecting the following 8 tick species:Dermacentor variabilis, Dermacentor occidentalis, Dermacentor andersoni, Dermacentor albipictus, Dermacentor parumapertus, Haemaphysalis leporispalustris, Ornithodoros concanensis,andArgas cooleyi [5]. In recent years, investigators in Latin America have identifiedR. belliiin ixodid ticks throughout Latin America, including El Salvador, Costa Rica, Panama, Colombia, Brazil, Peru, and Argentina [6–13]. In Central and South America,R. belliihas been detected in at least 19 tick species, namely,Ixodes loricatus, Haemaphysalis juxtakochi, Amblyomma aureolatum, Amblyomma dubitatum, Amblyomma humerale, Amblyomma incisum, Amblyomma neumanni, Amblyomma longirostre, Amblyomma naponense, Amblyomma nodosum, Amblyomma ovale, Amblyomma oblongoguttatum, Amblyomma parvum, Amblyomma pseudoconcolor, Amblyomma rotundatum, Amblyomma sabanerae, Amblyomma scalpturatum, Amblyomma tigrinum, andAmblyomma varium [6, 7, 12, 14–16].The wide host range and extraordinarily broad distribution of thisRickettsia species are intriguing and althoughR. bellii is generally considered as nonpathogenic for animals and humans [2, 17], it could play an important role in the ecology and epidemiology of other pathogenic tick-borne SFG rickettsiae in the Americas [18]. Despite the diverse representation ofR. bellii across more than 25 species of hard and soft ticks in the American continent, the phylogeographical relationships of this basal group-Rickettsia species has not been examined; the work described here explores these relationships, based on a hypothesis thatR. belliicould have evolved firstly in one continent or a particular tick group and then radiated to other continents or other tick groups.


## 2. Materials and Methods

### 2.1. Rickettsial Isolates

Rickettsia bellii was isolated from various species of hard ticks collected in South America and North America, as described in the original publication of each of the isolates listed in Table 1. A total of 30 isolates ofR. bellii from 13 species of hard ticks from 3 countries were used for the analysis (Figure 1). These included 1 from Argentina, 14 from Brazil, and 15 from USA (Table 1). The Brazilian and Argentinean isolates are available at the Rickettsial Collection of the Laboratory of Parasitic Diseases of Faculdade de Medicina Veterinária e Zootecnia (FMVZ) of the University of São Paulo (USP), and the US isolates are available at the Rickettsial Zoonoses Branch at the Centers for Disease Control and Prevention (CDC), Atlanta, Georgia, USA.Table 1
Rickettsia bellii isolates from South and North America used in the present study.
Number∗ Isolate name Tick host Geographic origin Reference (1) Mogi Amblyomma aureolatum Mogi das Cruzes-SP, Brazil [19] (2) Ad-MG Amblyomma dubitatum Guarda-Mor-MG, Brazil [16] (3) Ad-CORD-SP A. dubitatum Cordeirópolis -SP, Brazil [20] (4) Ad-25-INT A. dubitatum Ribeirão Grande-SP, Brazil [20] (5) PNSM Amblyomma incisum Cubatão-SP, Brazil [21] (6) AO Amblyomma ovale Ribeirão Grande-SP, Brazil [22] (7) HJ#4 Haemaphysalis juxtakochi Ribeirão Grande-SP, Brazil [23] (8) IL-Mogi Ixodes loricatus Mogi das Cruzes-SP, Brazil [24] (9) Ap GSV 136 Amblyomma parvum Chapada Gaúcha-MG, Brazil [25] (10) IL-RS1 I. loricatus Derrubadas-RS, Brazil [26] (11) RB-CL I. loricatus Cerro Largo-RS, Brazil [13] (12) HJ#1 H. juxtakochi São Paulo-SP, Brazil [23] (13) Ap-MS Amblyomma pseudoconcolor Corumbá-MS, Brazil [16] (14) P-10 A. ovale Peruíbe-SP, Brazil [27] (15) An4 Amblyomma neumanni Deán Funes-Córdoba, Argentina [28] (16) 369-C Dermacentor variabilis Washington Co., Arkansas, USA [5] (17) CA13-1 Dermacentor variabilis Yolo Co., California, USA [29] (18) CA13-9 D. variabilis Yolo Co., California, USA [29] (19) CA13-17 D. variabilis Yolo Co., California, USA [29] (20) Putah Creek D. variabilis Solano Co., California, USA [29] (21) Yolo D. variabilis Yolo Co., California, USA [29] (22) Stevenson Bridge D. variabilis Yolo Co., California, USA [29] (23) OSU 83-1223 D. variabilis Knox County, Ohio, USA [30] (24) TX15-1 Dermacentor parumapertus Brewster County, Texas, USA [31] (25) OSU 83-117 D. variabilis Licking County, Ohio, USA [30] (26) OSU 85-1299 D. variabilis Morrow County, Ohio, USA [30] (27) OSU 83-452 D. variabilis Coshocton County, Ohio, USA [30] (28) Skull Valley D. parumapertus Tooele County, Utah, USA [31] (29) OSU 85-389 D. variabilis Ohio, USA [30] (30) CA-459 Haemaphysalis leporispalustris Mendocino Co., California, USA [32] ∗These numbers are represented in Figure 1 as the geographical location of each isolate.Figure 1
Molecular phylogenetic analysis of 30 isolates ofRickettsia bellii from North and South America.


### 2.2. DNA Extraction and PCR

DNA from theR. bellii isolates from South America was extracted using the DNeasy Blood and Tissue Kit (Qiagen, Valencia, CA) and the QIAamp DNA Mini Kit (Qiagen) was used for the 15 isolates from North America, all in accordance with the manufacturer’s recommendations. Amplification of fragments of five rickettsial genes and fifteen intergenic regions was attempted with the primer pairs listed in Table 2. Each PCR reaction consisted of 2 μl template DNA, 20 picomoles of each primer, and 10 μl Taq PCR Master Mix (Qiagen), while cycling conditions included a 1-minute incubation at 95°C followed by 35 cycles of a 30-second denaturation at 95°C, a 30-second annealing incubation (Table 2), and a 1-minute extension at 72°C. This was followed by a final 10-minute extension at 72°C. Gradient PCR was used to optimize the annealing temperatures inR. bellii for each primer pair.Table 2
Primer pairs used for amplification of rickettsial genes or intergenic regions in the present study.
Primer pair Target Forward primer (5′ to 3′) Reverse primer (5′ to 3′) Annealing temperature (°C) Amplicon size (nt) Reference (1) mppA-purC GCAATTATCGGTCCGAATG TTTCATTTATTTGTCTCAAAATTCA 45 to 551 160 [33] (2) nusG-rplK CAGTTGCAATATTGGTAAAGCA CAGCAGCTGGAATTATCAAGTT 45 to 551 270 [33] (3) murG-RC0563 GAAGAAAAGAAGGGCATAAGCTA CAAGCTGAAAGTAAAAACATTCC 40 to 521 293 [33] (4) ntrY-rpsU AGCTGCTGTTGCTAAAGTAAAAA CAAGAAGCAGCAAGAAGACAGA 52 to 641 363 [33] (5) dksA-xerC TCCCATAGGTAATTTAGGTGTTTC TACTACCGCATATCCAATTAAAAA 45 to 551 416 [33] (6) spo0J-abcT1 AAAGATTTGGAAGAATTAGACTTGAT TTTGCTTAAACCAACCATTTCA 45 to 551 259 [33] (7) fabZ-lpxD TGTTAGGATCGATTTTAAGTACTCTATCT TGGATTGGCATAGACAATCTATTA 45 to 551 195 [33] (8) RC1137-tlc-5 CGGGATAACGCCGAGTAATA ATGCCGCTCTGAATTTGTTT 45 to 551 264 [33] (9) RC0230-RC0231 TGCACCCGCCTAAAACTAAC ATGGTCGGCCGTAGAAAAA 45 to 551 232 [33] (10) groES-RC0970 CTTGCATCGGCTTTTCTTTT AGCTTTGAGCTGATGGGCTA 45 to 551 215 [33] (11) rrf-pyrH GAGCTTTCTCCATCTTTTCTTG AAAGGGGAATATACGACAATTGAG 45 to 551 238 [33] (12) t R N A G l y-tRNATyr AGCTTGGAAGGCTGGAACTC ATCCTTCTCCCTCCACCACT 45 to 551 148 [33] (13) pal-RC1201 TGCAAGCACACATAATGCAA TCAAAATCGATTCCTCTTTTCC 45 to 551 216 [33] (14) atpA ATCAAGCGTTGCACAGATAG GGAAGTGCCGTAAGTGAACC 58 Unknown∗ [1] (15) rpmE-tRNAfMet TTCCGGAAATGTAGTAAATCAATC TCAGGTTATGAGCCTGACGA 54 144 [33] (16) RC1027-xthA2 GGTATGTAAATGAGCCTTATCAATACT TCAGTAGTATAAGTAGCTCCTGCTGTC 54 351 [33] (17) gltA GTCTACTGCTTCGTGTAGATCAAC GGCTGACCTATAGAATATTTATAAGAC 54 408 [25] (18) coxA ACAGCCGTTGATATGGCTA CATATTCCAACCGGCAAAAG 58 Unknown∗ [1] (19) coxA GGTGCTCCTGATATGGCATT CATATTCCAGCCGGCAAAAG 58 Unknown∗ [1] (20) atpA ATCAAGCGTTGCACAGATAG CGACTTACCGAAATACCGAC 56 449 [1, 34] ∗Data not given in the original publication; 1tested in gradient PCR.


### 2.3. DNA Purification and Sequencing

DNA fragments amplified by PCR were visualized using a UV lamp in a 1.5% agarose gel containing 0.1μg/ml ethidium bromide. PCR products of the appropriate size were cut from the gel and then purified using the Wizard SV gel and PCR clean-up system (Promega, Madison, WI). Sequencing reactions were prepared using one microliter of purified PCR product and the BigDye Terminator v3.1 Cycle Sequencing Kit (Applied Biosystems, Foster City, CA) according to the manufacturer’s instructions and then sequenced using an ABI 3100 genetic analyzer (Applied Biosystems). Each PCR amplicon was sequenced at least once in both directions.


### 2.4. Sequence Alignment and Phylogenetic Inferences

The resulting sequences were assembled using Geneious® 9.1.4 software (Biomatters Ltd., Auckland, New Zealand) and aligned using ClustalW [35] and MEGA 6.0.6 software [36]. The resulting alignment was examined by eye to ensure proper alignment of the sequences and the simple indel-coding method [37] was used to remove insertions/deletions. MrBayes 3.2.6 [38] was used in Geneious® to perform a Bayesian phylogenetic analysis. The Jukes-Cantor model was used in an analysis consisting of 10,000,000 generations (1,000,000 burn-ins). The analysis included 3 heated chains, and an effective sample size (ESS) of 901 was achieved for all parameters.


### 2.5. Accession Numbers

The GenBank accession numbers for the DNA sequences generated in this study for the 30R. belliiisolates shown in Table 1 are the following:gltA gene (MF154866–MF154895),atpA gene (MF154926–MF154955),coxA gene (MF154896–MF154925),rpmE-tRNAfmet intergenic region (MF154956–MF154985), andRC1027-xthA2intergenic region (MF154986–MF155015).


## 2.1. Rickettsial Isolates

Rickettsia bellii was isolated from various species of hard ticks collected in South America and North America, as described in the original publication of each of the isolates listed in Table 1. A total of 30 isolates ofR. bellii from 13 species of hard ticks from 3 countries were used for the analysis (Figure 1). These included 1 from Argentina, 14 from Brazil, and 15 from USA (Table 1). The Brazilian and Argentinean isolates are available at the Rickettsial Collection of the Laboratory of Parasitic Diseases of Faculdade de Medicina Veterinária e Zootecnia (FMVZ) of the University of São Paulo (USP), and the US isolates are available at the Rickettsial Zoonoses Branch at the Centers for Disease Control and Prevention (CDC), Atlanta, Georgia, USA.Table 1
Rickettsia bellii isolates from South and North America used in the present study.
Number∗ Isolate name Tick host Geographic origin Reference (1) Mogi Amblyomma aureolatum Mogi das Cruzes-SP, Brazil [19] (2) Ad-MG Amblyomma dubitatum Guarda-Mor-MG, Brazil [16] (3) Ad-CORD-SP A. dubitatum Cordeirópolis -SP, Brazil [20] (4) Ad-25-INT A. dubitatum Ribeirão Grande-SP, Brazil [20] (5) PNSM Amblyomma incisum Cubatão-SP, Brazil [21] (6) AO Amblyomma ovale Ribeirão Grande-SP, Brazil [22] (7) HJ#4 Haemaphysalis juxtakochi Ribeirão Grande-SP, Brazil [23] (8) IL-Mogi Ixodes loricatus Mogi das Cruzes-SP, Brazil [24] (9) Ap GSV 136 Amblyomma parvum Chapada Gaúcha-MG, Brazil [25] (10) IL-RS1 I. loricatus Derrubadas-RS, Brazil [26] (11) RB-CL I. loricatus Cerro Largo-RS, Brazil [13] (12) HJ#1 H. juxtakochi São Paulo-SP, Brazil [23] (13) Ap-MS Amblyomma pseudoconcolor Corumbá-MS, Brazil [16] (14) P-10 A. ovale Peruíbe-SP, Brazil [27] (15) An4 Amblyomma neumanni Deán Funes-Córdoba, Argentina [28] (16) 369-C Dermacentor variabilis Washington Co., Arkansas, USA [5] (17) CA13-1 Dermacentor variabilis Yolo Co., California, USA [29] (18) CA13-9 D. variabilis Yolo Co., California, USA [29] (19) CA13-17 D. variabilis Yolo Co., California, USA [29] (20) Putah Creek D. variabilis Solano Co., California, USA [29] (21) Yolo D. variabilis Yolo Co., California, USA [29] (22) Stevenson Bridge D. variabilis Yolo Co., California, USA [29] (23) OSU 83-1223 D. variabilis Knox County, Ohio, USA [30] (24) TX15-1 Dermacentor parumapertus Brewster County, Texas, USA [31] (25) OSU 83-117 D. variabilis Licking County, Ohio, USA [30] (26) OSU 85-1299 D. variabilis Morrow County, Ohio, USA [30] (27) OSU 83-452 D. variabilis Coshocton County, Ohio, USA [30] (28) Skull Valley D. parumapertus Tooele County, Utah, USA [31] (29) OSU 85-389 D. variabilis Ohio, USA [30] (30) CA-459 Haemaphysalis leporispalustris Mendocino Co., California, USA [32] ∗These numbers are represented in Figure 1 as the geographical location of each isolate.Figure 1
Molecular phylogenetic analysis of 30 isolates ofRickettsia bellii from North and South America.


## 2.2. DNA Extraction and PCR

DNA from theR. bellii isolates from South America was extracted using the DNeasy Blood and Tissue Kit (Qiagen, Valencia, CA) and the QIAamp DNA Mini Kit (Qiagen) was used for the 15 isolates from North America, all in accordance with the manufacturer’s recommendations. Amplification of fragments of five rickettsial genes and fifteen intergenic regions was attempted with the primer pairs listed in Table 2. Each PCR reaction consisted of 2 μl template DNA, 20 picomoles of each primer, and 10 μl Taq PCR Master Mix (Qiagen), while cycling conditions included a 1-minute incubation at 95°C followed by 35 cycles of a 30-second denaturation at 95°C, a 30-second annealing incubation (Table 2), and a 1-minute extension at 72°C. This was followed by a final 10-minute extension at 72°C. Gradient PCR was used to optimize the annealing temperatures inR. bellii for each primer pair.Table 2
Primer pairs used for amplification of rickettsial genes or intergenic regions in the present study.
Primer pair Target Forward primer (5′ to 3′) Reverse primer (5′ to 3′) Annealing temperature (°C) Amplicon size (nt) Reference (1) mppA-purC GCAATTATCGGTCCGAATG TTTCATTTATTTGTCTCAAAATTCA 45 to 551 160 [33] (2) nusG-rplK CAGTTGCAATATTGGTAAAGCA CAGCAGCTGGAATTATCAAGTT 45 to 551 270 [33] (3) murG-RC0563 GAAGAAAAGAAGGGCATAAGCTA CAAGCTGAAAGTAAAAACATTCC 40 to 521 293 [33] (4) ntrY-rpsU AGCTGCTGTTGCTAAAGTAAAAA CAAGAAGCAGCAAGAAGACAGA 52 to 641 363 [33] (5) dksA-xerC TCCCATAGGTAATTTAGGTGTTTC TACTACCGCATATCCAATTAAAAA 45 to 551 416 [33] (6) spo0J-abcT1 AAAGATTTGGAAGAATTAGACTTGAT TTTGCTTAAACCAACCATTTCA 45 to 551 259 [33] (7) fabZ-lpxD TGTTAGGATCGATTTTAAGTACTCTATCT TGGATTGGCATAGACAATCTATTA 45 to 551 195 [33] (8) RC1137-tlc-5 CGGGATAACGCCGAGTAATA ATGCCGCTCTGAATTTGTTT 45 to 551 264 [33] (9) RC0230-RC0231 TGCACCCGCCTAAAACTAAC ATGGTCGGCCGTAGAAAAA 45 to 551 232 [33] (10) groES-RC0970 CTTGCATCGGCTTTTCTTTT AGCTTTGAGCTGATGGGCTA 45 to 551 215 [33] (11) rrf-pyrH GAGCTTTCTCCATCTTTTCTTG AAAGGGGAATATACGACAATTGAG 45 to 551 238 [33] (12) t R N A G l y-tRNATyr AGCTTGGAAGGCTGGAACTC ATCCTTCTCCCTCCACCACT 45 to 551 148 [33] (13) pal-RC1201 TGCAAGCACACATAATGCAA TCAAAATCGATTCCTCTTTTCC 45 to 551 216 [33] (14) atpA ATCAAGCGTTGCACAGATAG GGAAGTGCCGTAAGTGAACC 58 Unknown∗ [1] (15) rpmE-tRNAfMet TTCCGGAAATGTAGTAAATCAATC TCAGGTTATGAGCCTGACGA 54 144 [33] (16) RC1027-xthA2 GGTATGTAAATGAGCCTTATCAATACT TCAGTAGTATAAGTAGCTCCTGCTGTC 54 351 [33] (17) gltA GTCTACTGCTTCGTGTAGATCAAC GGCTGACCTATAGAATATTTATAAGAC 54 408 [25] (18) coxA ACAGCCGTTGATATGGCTA CATATTCCAACCGGCAAAAG 58 Unknown∗ [1] (19) coxA GGTGCTCCTGATATGGCATT CATATTCCAGCCGGCAAAAG 58 Unknown∗ [1] (20) atpA ATCAAGCGTTGCACAGATAG CGACTTACCGAAATACCGAC 56 449 [1, 34] ∗Data not given in the original publication; 1tested in gradient PCR.


## 2.3. DNA Purification and Sequencing

DNA fragments amplified by PCR were visualized using a UV lamp in a 1.5% agarose gel containing 0.1μg/ml ethidium bromide. PCR products of the appropriate size were cut from the gel and then purified using the Wizard SV gel and PCR clean-up system (Promega, Madison, WI). Sequencing reactions were prepared using one microliter of purified PCR product and the BigDye Terminator v3.1 Cycle Sequencing Kit (Applied Biosystems, Foster City, CA) according to the manufacturer’s instructions and then sequenced using an ABI 3100 genetic analyzer (Applied Biosystems). Each PCR amplicon was sequenced at least once in both directions.


## 2.4. Sequence Alignment and Phylogenetic Inferences

The resulting sequences were assembled using Geneious® 9.1.4 software (Biomatters Ltd., Auckland, New Zealand) and aligned using ClustalW [35] and MEGA 6.0.6 software [36]. The resulting alignment was examined by eye to ensure proper alignment of the sequences and the simple indel-coding method [37] was used to remove insertions/deletions. MrBayes 3.2.6 [38] was used in Geneious® to perform a Bayesian phylogenetic analysis. The Jukes-Cantor model was used in an analysis consisting of 10,000,000 generations (1,000,000 burn-ins). The analysis included 3 heated chains, and an effective sample size (ESS) of 901 was achieved for all parameters.


## 2.5. Accession Numbers

The GenBank accession numbers for the DNA sequences generated in this study for the 30R. belliiisolates shown in Table 1 are the following:gltA gene (MF154866–MF154895),atpA gene (MF154926–MF154955),coxA gene (MF154896–MF154925),rpmE-tRNAfmet intergenic region (MF154956–MF154985), andRC1027-xthA2intergenic region (MF154986–MF155015).


## 3. Results

Among the 20 primer pairs used for amplification of rickettsial DNA fragments, only the primer pairs # 15, 16, 17, 19, and 20 (Table2) were successful in amplifying DNA from theR. belliiisolates. These primers corresponded to 3 protein coding genes (gltA,atpA,andcoxA) and 2 intergenic regions (rpmE-tRNAfmet and RC1027-xthA2). The nucleotide sequences from the five loci for each of the 30 isolates were concatenated and used to perform a Bayesian analysis of their phylogeny. There was a separation of the isolates into at least three clades: one representing all South American isolates and the United States isolate CA-459 cultivated from aH. leporispalustris tick collected in northern California [39], one comprising 6 isolates fromD. variabilis collected in northern California [29], and a third clade comprising 5 isolates fromD. variabilis collected in Ohio [30] (Figure 1). North American isolates fromD. parumapertusand the 369-C strain fromD. variabilis did not group within any of the above clades. Within the South American clade, there were two smaller well-supported clades (each with a posterior probability of 1.0), one consisting of three isolates associated withI. loricatusticks and one consisting of twoA. ovale-associated isolates.Overall, there were very few polymorphisms (<0.5%) between the 30R. belliiisolates. In the 2,269-nucleotide alignment, identity values between the isolates varied from 99.56 to 100% (Table S1). DNA sequences of the isolates composing the main South American clade were 99.69–100% identical and at the same time 99.56–99.74% identical to the clade containingD. variabilisisolates from California and 99.74–99.91% identical to the clade containingD. variabilisisolates from Ohio. Within-clade identities were 99.74–100% for theD. variabilisisolates from Ohio and 100% for theD. variabilisisolates from California; the Ohio clade was 99.74% identical to the California clade. Finally, sequences of the three North American isolates (two fromD. parumapertusand one fromD. variabilis) that did not group within any main clade were 99.87–99.96% identical to each other and 99.60–99.91% identical to the remaining isolates.


## 4. Discussion

Here we performed for the first time a phylogenetic analysis of multiple isolates ofR. bellii.Because these isolates represented a number of different locations in North and South America, our intention was to infer phylogeographical relationships. Indeed, there was a separation of almost all isolates between the two continents, suggesting that they have radiated within their respective continents. Although the posterior probability of the separation of the South American clade from the North American isolates was not strong (PP = 0.59), this topology was also achieved using both the GTR and HKY85 evolutionary models (supplemental material), giving additional evidence for this separation. At the same time, it is noteworthy that while the North American clades represented isolates exclusively fromDermacentorticks, the South American clade was composed of isolates from the generaAmblyomma, Haemaphysalis,andIxodes.Under such circumstances, the phylogenetic positions of the 30 isolates of the present study could be a result of not only their geographical origin but also the tick hosts in which they were isolated from. This later assumption is corroborated by the position of the isolate fromH. leporispalustris, which, despite being from North America, grouped within the South American clade, where at least two other isolates fromHaemaphysalisticks were present.To our knowledge, the number and diversity of tick species infected withR. belliiare the largest and broadest described among species in the genusRickettsia. The magnitude of different host species indicates horizontal transmission among tick populations, especially in South America, where isolates from different tick species and genera grouped together in the same clade. It is generally accepted thatR. belliiis not pathogenic for vertebrates; in fact,R. belliihas never been isolated from a vertebrate host [2]. On the other hand, there has been serological evidence of animal natural infection or exposure byR. bellii [40, 41]. WhileR. belliihorizontal transmission among ticks via vertebrate host cannot be discarded, another likely mechanism could be via tick parasitoids, since the parasitism of the hymenopteranIxodiphagusspp. is relatively common among different tick genera in the Americas [42, 43]. In fact, a recent study provided molecular detection of tick-borne rickettsiae inIxodiphaguswasps that had emerged from ticks [44], highlighting the possibility that rickettsial organisms could be shared by ticks and their parasitoids. Regardless of the mechanisms of horizontal transmission, the tendency of eachR. belliigenotype to be associated with a different tick species (Figure 1) suggests that horizontal transmission was more efficient at earlier times; thereafter, most of theR. belliiisolates are likely to have coevolved specifically with their specific tick species host, possibly towards a symbiotic association. In fact, analysis of the genome of the type strain ofR. belliihas shown features compatible with various symbiotic bacteria, such as a large genome size with high coding capacity, in contrast to the reduced genome size with low coding capacity of the pathogenic rickettsiae [4, 45].It has been proposed that theR. bellii group diverged prior to the division between the SFG and the TG, forming a basal group that also includes various herbivorous arthropod symbionts [1, 46]. While the genusRickettsiawas considered to be approximately 150 million years old, the divergence of theR. bellii group was estimated to have occurred much more recently, around 50 million years ago [1]. Analyses of the oldest tick fossils, from the Cretaceous, indicate that extant tick genera(Amblyomma, Ornithodoros) were soundly established around 100 million years ago [47–49]. From these data, it can be inferred thatR. belliilikely radiated with its principal hosts approximately 50 million years ago. While ticks and rickettsiae are distributed in all continents, it is noteworthy that, under natural conditions and without human interference, tick species from the New World do not occur in the Old World and vice versa; the only exceptions are a few tick species that are associated with transcontinental marine birds [50]. Indeed, this scenario has accounted for the presence ofR. bellii-infected ticks restricted to the New World, regardless of the ability of this bacterium to infect a vast array of tick genera and species. On the other hand, whetherR. belliioriginated with ticks in North or South America remains obscure, since our phylogenetic analyses did not show any evidence for greater genetic divergence ofR. belliiin any of the two continents. Further analyses encompassing moreR. belliiisolates from different tick genera and species, encompassing more molecular markers, should provide cues for the origin and radiation ofR. belliiin the New World.


## 5. Conclusion

Phylogeographical analysis of 30 strains (15 from North America and 15 from South America) isolated from 13 species of 4 genera shows a clear differentiation between most of the North and South American isolates, indicating geographic isolation between isolates of these two continents. Additionally, this analysis separated isolates ofR. bellii by the species of tick from which these were isolated, indicating that the isolates could have coevolved with their tick vectors over time.



---

*Source: 8505483-2018-04-08.xml*