| Literature DB >> 27031729 |
Volker Fingerle1, Michael Pritsch2, Martin Wächtler3, Gabriele Margos1, Sabine Ruske2, Jette Jung4, Thomas Löscher2, Clemens Wendtner3, Andreas Wieser2,4,5.
Abstract
A 26 year-old female patient presented to the Tropical Medicine outpatient unit of the Ludwig Maximilians-University in Munich with febrile illness after returning from Southern Africa, where she contracted a bite by a large mite-like arthropod, most likely a soft-tick. Spirochetes were detected in Giemsa stained blood smears and treatment was started with doxycycline for suspected tick-borne relapsing fever. The patient eventually recovered after developing a slight Jarisch-Herxheimer reaction during therapy. PCR reactions performed from EDTA-blood revealed a 16S rRNA sequence with 99.4% similarity to both, Borrelia duttonii, and B. parkeri. Further sequences obtained from the flagellin gene (flaB) demonstrated genetic distances of 0.066 and 0.097 to B. parkeri and B. duttonii, respectively. Fragments of the uvrA gene revealed genetic distance of 0.086 to B. hermsii in genetic analysis and only distant relations with classic Old World relapsing fever species. This revealed the presence of a novel species of tick-borne relapsing fever spirochetes that we propose to name "Candidatus Borrelia kalaharica", as it was contracted from an arthropod bite in the Kalahari Desert belonging to both, Botswana and Namibia, a region where to our knowledge no relapsing fever has been described so far. Interestingly, the novel species shows more homology to New World relapsing fever Borrelia such as B. parkeri or B. hermsii than to known Old World species such as B. duttonii or B. crocidurae.Entities:
Mesh:
Substances:
Year: 2016 PMID: 27031729 PMCID: PMC4816561 DOI: 10.1371/journal.pntd.0004559
Source DB: PubMed Journal: PLoS Negl Trop Dis ISSN: 1935-2727
Estimates of evolutionary divergence between 16S rRNA sequences.
The number of base substitutions per site from between sequences are shown. Analyses were conducted using the Kimura 2-parameter model. There were a total of 473 positions in the final dataset.
| 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | 11 | 12 | 13 | 14 | 15 | 16 | 17 | 18 | 19 | 20 | ||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 1 | |||||||||||||||||||||
| 2 | gi|49387237| | 0,004 | |||||||||||||||||||
| 3 | gi|576098524| | 0,006 | 0,006 | ||||||||||||||||||
| 4 | gi|507148154| | 0,006 | 0,006 | 0,004 | |||||||||||||||||
| 5 | gi|201083369| | 0,006 | 0,006 | 0,004 | 0,000 | ||||||||||||||||
| 6 | gi|119861577| | 0,006 | 0,006 | 0,000 | 0,004 | 0,004 | |||||||||||||||
| 7 | gi|329130340| | 0,009 | 0,009 | 0,006 | 0,002 | 0,002 | 0,006 | ||||||||||||||
| 8 | gi|1814213| | 0,009 | 0,009 | 0,006 | 0,002 | 0,002 | 0,006 | 0,004 | |||||||||||||
| 9 | gi|1814216| | 0,011 | 0,011 | 0,009 | 0,004 | 0,004 | 0,009 | 0,006 | 0,006 | ||||||||||||
| 10 | gi|770587662| | 0,013 | 0,013 | 0,006 | 0,011 | 0,011 | 0,006 | 0,013 | 0,013 | 0,015 | |||||||||||
| 11 | gi|659364680| | 0,017 | 0,017 | 0,011 | 0,011 | 0,011 | 0,011 | 0,013 | 0,013 | 0,015 | 0,017 | ||||||||||
| 12 | gi|576099683| | 0,019 | 0,019 | 0,017 | 0,022 | 0,022 | 0,017 | 0,024 | 0,024 | 0,026 | 0,015 | 0,028 | |||||||||
| 13 | gi|607196840| | 0,019 | 0,019 | 0,013 | 0,017 | 0,017 | 0,013 | 0,019 | 0,019 | 0,022 | 0,015 | 0,006 | 0,026 | ||||||||
| 14 | gi|700323429| | 0,021 | 0,021 | 0,024 | 0,024 | 0,024 | 0,024 | 0,026 | 0,024 | 0,028 | 0,021 | 0,035 | 0,026 | 0,032 | |||||||
| 15 | gi|53801415| | 0,026 | 0,026 | 0,019 | 0,024 | 0,024 | 0,019 | 0,026 | 0,024 | 0,028 | 0,022 | 0,017 | 0,033 | 0,011 | 0,035 | ||||||
| 16 | gi|890832120 | 0,028 | 0,028 | 0,022 | 0,026 | 0,026 | 0,022 | 0,028 | 0,026 | 0,030 | 0,024 | 0,019 | 0,026 | 0,013 | 0,039 | 0,017 | |||||
| 17 | gi|793775| | 0,028 | 0,028 | 0,026 | 0,030 | 0,030 | 0,026 | 0,033 | 0,030 | 0,035 | 0,024 | 0,028 | 0,028 | 0,026 | 0,009 | 0,028 | 0,032 | ||||
| 18 | gi|672590398| | 0,030 | 0,035 | 0,033 | 0,037 | 0,037 | 0,033 | 0,039 | 0,037 | 0,042 | 0,030 | 0,044 | 0,030 | 0,042 | 0,015 | 0,044 | 0,044 | 0,015 | |||
| 19 | gi|342856386| | 0,030 | 0,035 | 0,033 | 0,037 | 0,037 | 0,033 | 0,039 | 0,037 | 0,042 | 0,030 | 0,044 | 0,026 | 0,042 | 0,015 | 0,044 | 0,039 | 0,015 | 0,009 | ||
| 20 | gi|793785| | 0,035 | 0,035 | 0,033 | 0,037 | 0,037 | 0,033 | 0,039 | 0,037 | 0,042 | 0,030 | 0,035 | 0,035 | 0,032 | 0,015 | 0,035 | 0,039 | 0,006 | 0,017 | 0,017 |
Estimates of evolutionary divergence between flaB sequences.
The number of base substitutions per site from between sequences are shown. Analyses were conducted using the Kimura 2-parameter model. The analysis involved 18 nucleotide sequences. There were a total of 640 positions in the final dataset.
| 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | 11 | 12 | 13 | 14 | 15 | 16 | 17 | 18 | ||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 1 | |||||||||||||||||||
| 2 | gi|576100539| | 0,059 | |||||||||||||||||
| 3 | gi|119861577| | 0,064 | 0,064 | ||||||||||||||||
| 4 | gi|770587662| | 0,071 | 0,069 | 0,071 | |||||||||||||||
| 5 | gi|576098524| | 0,066 | 0,069 | 0,011 | 0,071 | ||||||||||||||
| 6 | gi|576093260| | 0,080 | 0,079 | 0,074 | 0,085 | 0,069 | |||||||||||||
| 7 | gi|300086788| | 0,087 | 0,093 | 0,101 | 0,111 | 0,099 | 0,101 | ||||||||||||
| 8 | gi|201084691| | 0,098 | 0,104 | 0,112 | 0,118 | 0,110 | 0,116 | 0,022 | |||||||||||
| 9 | gi|201083369| | 0,097 | 0,104 | 0,112 | 0,118 | 0,110 | 0,112 | 0,019 | 0,003 | ||||||||||
| 10 | gi|384934107| | 0,097 | 0,104 | 0,109 | 0,114 | 0,110 | 0,119 | 0,024 | 0,014 | 0,014 | |||||||||
| 11 | gi|576102789| | 0,092 | 0,083 | 0,086 | 0,092 | 0,092 | 0,100 | 0,106 | 0,117 | 0,117 | 0,117 | ||||||||
| 12 | gi|333463402| | 0,099 | 0,102 | 0,114 | 0,120 | 0,112 | 0,114 | 0,021 | 0,005 | 0,002 | 0,016 | 0,119 | |||||||
| 13 | gi|57169094| | 0,112 | 0,116 | 0,119 | 0,127 | 0,123 | 0,117 | 0,130 | 0,136 | 0,136 | 0,136 | 0,087 | 0,134 | ||||||
| 14 | gi|342856386| | 0,141 | 0,141 | 0,150 | 0,158 | 0,156 | 0,168 | 0,162 | 0,176 | 0,176 | 0,168 | 0,176 | 0,178 | 0,194 | |||||
| 15 | gi|672590398| | 0,148 | 0,154 | 0,148 | 0,164 | 0,154 | 0,168 | 0,160 | 0,174 | 0,174 | 0,170 | 0,174 | 0,176 | 0,196 | 0,044 | ||||
| 16 | gi|858732| | 0,150 | 0,147 | 0,152 | 0,164 | 0,158 | 0,172 | 0,170 | 0,180 | 0,180 | 0,172 | 0,188 | 0,182 | 0,206 | 0,049 | 0,049 | |||
| 17 | gi|51572834| | 0,160 | 0,162 | 0,162 | 0,184 | 0,164 | 0,178 | 0,172 | 0,182 | 0,182 | 0,178 | 0,194 | 0,184 | 0,209 | 0,062 | 0,059 | 0,016 | ||
| 18 | gi|6626249| | 0,154 | 0,156 | 0,154 | 0,172 | 0,160 | 0,180 | 0,164 | 0,176 | 0,174 | 0,172 | 0,180 | 0,176 | 0,194 | 0,061 | 0,061 | 0,066 | 0,068 |
Fig 2Molecular phylogenetic analysis using the maximum likelihood method based on partial sequences of (A) 16S rRNA (473 bp), (B) . The taxa labels show the NCBI accession number, Borrelia species, and strain information (if available). The strain investigated in this study (indicated by a black dot)–although acquired in Africa–clusters more closely to New World RF species in flaB (B) and uvrA (C) phylogenies. The trees with the highest log likelihood are shown. The percentage of trees in which the associated taxa clustered together is shown next to the branches. Scale bar: substitutions per site. RF = relapsing fever; LB = Lyme borreliosis
Estimates of evolutionary divergence between uvrA sequences.
The number of base substitutions per site from between sequences are shown. There were a total of 900 positions in the final dataset.
| 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 | 9 | 10 | 11 | 12 | 13 | 14 | 15 | ||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| 1 | ||||||||||||||||
| 2 | gi|770587662| | 0,086 | ||||||||||||||
| 3 | gi|576099683| | 0,108 | 0,100 | |||||||||||||
| 4 | gi|576098524| | 0,110 | 0,092 | 0,128 | ||||||||||||
| 5 | gi|119861577| | 0,117 | 0,087 | 0,125 | 0,021 | |||||||||||
| 6 | gi|752506972 | 0,139 | 0,106 | 0,149 | 0,093 | 0,098 | ||||||||||
| 7 | 0,166 | 0,154 | 0,166 | 0,161 | 0,161 | 0,169 | ||||||||||
| 8 | gb|CP000976.1 | 0,168 | 0,155 | 0,160 | 0,157 | 0,156 | 0,167 | 0,046 | ||||||||
| 9 | gi|384934107| | 0,170 | 0,156 | 0,161 | 0,156 | 0,154 | 0,169 | 0,050 | 0,008 | |||||||
| 10 | gi|201084691| | 0,167 | 0,153 | 0,158 | 0,156 | 0,154 | 0,163 | 0,047 | 0,006 | 0,009 | ||||||
| 11 | gi|700323429| | 0,266 | 0,247 | 0,269 | 0,255 | 0,248 | 0,261 | 0,281 | 0,256 | 0,254 | 0,256 | |||||
| 12 | gi|672590398| | 0,270 | 0,247 | 0,263 | 0,250 | 0,240 | 0,245 | 0,260 | 0,253 | 0,250 | 0,253 | 0,075 | ||||
| 13 | gi|342856386| | 0,258 | 0,247 | 0,267 | 0,248 | 0,243 | 0,253 | 0,266 | 0,259 | 0,256 | 0,259 | 0,081 | 0,070 | |||
| 14 | gi|635290980| | 0,275 | 0,255 | 0,271 | 0,257 | 0,252 | 0,265 | 0,275 | 0,258 | 0,256 | 0,256 | 0,076 | 0,061 | 0,066 | ||
| 15 | gi|51572834| | 0,269 | 0,255 | 0,266 | 0,253 | 0,243 | 0,263 | 0,268 | 0,253 | 0,251 | 0,251 | 0,080 | 0,059 | 0,067 | 0,017 |
Clinical laboratory values during the course of hospitalization.
| Value (normal range) | Day 1 | Day 2 | Day 3 | Day 4 | |
|---|---|---|---|---|---|
| 11am | 7pm | 8am | 10am | 11am | |
| Leucocytes/μl (3500–9800) | 7000 | 5800 | 8800 | 7100 | 6500 |
| Erythrocytes/μl (4.1–5.1) | 5.2 | 4.7 | 4.1 | 4.0 | 4.3 |
| Haematocrit % (36.0–48.0) | 44.9 | 41.4 | 36.5 | 36.6 | 39.3 |
| Hemoglobin g/dl (12.0–16.0) | 15.4 | 14.3 | 12.1 | 12.1 | 13.2 |
| Thrombocytes /μl (140 000–360 000) | 188 000 | 157 000 | 137 000 | 125 000 | 192 000 |
| CRP mg/dl (<5.0) | 4.7 | 6.5 | not done | 8.0 | 4.1 |