Brucella neotomae infection in humans, Costa Rica.
B. neotomae, isolated in 1957 from wood rats (Neotoma lepida) in North America (2), has been considered nonzoonotic (3). It has been isolated from target organs of experimentally infected mice and guinea pigs (4,5). We report the isolation of B. neotomae from cerebrospinal fluid samples from 2 human patients with neurobrucellosis.
In 2008, a Brucella sp. isolate was submitted to the Tropical Diseases Research Center at the Universidad de Costa Rica. This isolate (denoted strain bneohCR2) was cultured from a cerebrospinal fluid sample obtained from a 64-year-old male patient at one of the main hospitals in San Jose, Costa Rica. In 2011, another isolate (denoted strain bneohCR1) was recovered from a cerebrospinal fluid sample from a 51-year-old male patient at a regional hospital in Costa Rica. As is common for other patients with brucellosis, the blood leukocyte count for each patient was almost within the reference range, and C-reactive protein level was within reference range. Both patients showed clinical signs compatible with neurobrucellosis (6), had positive Rose Bengal test results, and recovered after receiving 1 month of streptomycin (750 mg/d) and 3 months of doxycycline (100 mg/12 h).
Further bacteriologic analysis (7,8) confirmed that the isolates were a Brucella sp. (Table). Real-time PCR high-resolution melting analysis (9) confirmed genus designation but was inconclusive regarding species designation. Bruce-ladder multiplex PCR (10) and multiple loci variable number of tandem repeats-16-loci panel analysis (http:// mlva.u-psud.fr/brucella/;_Figure 1) indicated that the profiles for both DNA isolates corresponded to profiles for B. neotomae. Analysis of bneohCR2 by multiplex single-nucleotide polymorphism (SNP) primer extension assay (11) and by matrix-assisted laser desorption/ionization time-of-flight mass spectrometry of protein extracts (12) (online Technical Appendix Table 1, https://wwwnc.cdc.gov/EID/ article/23/6/16-2018-Techapp1.pdf) confirmed that the isolate was B. neotomae.
We performed whole-genome sequencing of bneohCR1 and bneohCR2 and resequencing of reference strain B. neotomae 5K33. Sequencing data were deposited at the European Nucleotide Archive (http://www.ebi.ac.uk/ena/) under accession codes ERS1563929 (bneohCR1), ERS1563928 (bneohCR2), and ERS1624467 (SK33). To place the bneohCR1 and bneohCR2 in a phylogenetic context, publicly available reads from 51 Brucella whole-genome sequences (online Technical Appendix Table 2) were aligned and then mapped to B. suis 1330 by using SMALT version 0.5.8 (ftp:// ftp.sanger.ac.uk/pub/resources/software/smalt/). Reads from bneohCR1 and bneohCR2 genomes mapped to 98.6% of the B. suis 1330 genome. SNPs were called from the alignment by use of Samtools (http://samtools.sourceforge.net/), and 34,307 variable sites across all isolates were extracted by using SNP sites (13). The resulting alignment of SNPs (https://github.com/sanger-pathogens/Farm_blast) comparison of this pseudo-molecule and the bneohCR draft genomes, generated by assembly with Velvet (13), showed that the genomic loci known as 26.5 kb, 12 kb, and GI-6 that are absent in B. neotomae (14) are also absent in the queried genomes.
This report of B. neotomae as a cause of zoonotic disease raises questions about possible underrepresentation of reported cases. This study also has implications for brucellosis diagnosis. Specifically, the differences among B. neotomae and the other Brucella species at the biochemical level are subtle. The major difference between B. neotomae and B. was used for maximum-likelihood phylogenetic reconstruction by use of RAxML version 7.0.4 (https://github.com/ stamatak/standard-RAxML). The generated phylogenetic tree confirmed that the bneohCR isolates clustered together with B. neotomae reference strain 5K33 (ENA accession no. JMSC01, assembly accession no. GCA_00742255.1) (Figure 2). Isolates bneohCR1 and bneohCR2 differed from the reference genome by 174 and 160 SNPs, respectively. This number of SNPs is smaller than that between B. abortus 9-941 and B. abortus 2308 (214 SNPs), which are 2 well-recognized strains of the same Brucella species.
Analysis of 23 previously reported genomic islands or anomalous genomic loci (14) was performed for both bneohCR genomes. For this analysis, a "genomic-island pseudo-molecule" was constructed by concatenation of 23 genomic regions obtained from different Brucella genomes. BLAST abortus, the main cause of human brucellosis in Costa Rica, is the presence of oxidase activity in B. abortus, which is assessed subjectively (7,8). Because B. neotomae has not, until now, been considered zoonotic, some cases of brucellosis reported as being caused by atypical zoonotic classical Brucella might have been misdiagnosed cases of B. neotomae infection. The introduction of whole-genome sequencing into the clinical field will thus improve diagnosis.
A lack of epidemiologic information with regard to the 2 patients reported here precluded the investigation of exposure or contact with hosts known to harbor B. neotomae. Further studies are needed to establish which animals may act as reservoirs for B. neotomae in Costa Rica.
In summary, B. neotomae should be considered a zoonosis risk for infection in humans. Incorporation of molecular techniques for diagnosis will help resolve the Brucella genus homogeneity obtained when only biochemical assays are used.
We are grateful to Heilyn Estrada-Murillo for isolate information, Daphnne Garita for technical assistance, Gabriela Hernandez-Mora for provision of serotype-specific antiserum, and Gordon Dougan for helpful discussions.
This work was funded by Fondos del Sistema del Consejo Nacional de Rectores, Costa Rica, and Wellcome Trust. N.R.-V. and C.J.-R. were partially sponsored by University of Costa Rica scholarships. K.S.B. is funded by Wellcome Trust Fellowship 106690/A/14/Z. K.S.B. and N.R.T. were funded by Wellcome Trust grant no. 098051.
All procedures involving live Brucella were carried out according to the "Reglamento de Bioseguridad de la CCSS39975-0," year 2012, after "Decreto Ejecutivo no. 0965-S," year 2002 and protocol approved by SIA0434-14 Universidad Nacional, Costa Rica. These genetic resources were accessed in Costa Rica according to the Biodiversity Law no. 7788 and the Convention on Biological Diversity, under the terms of respect to equal and fair distribution of benefits among those who provided such resources.
Miss Suarez-Esquivel is a PhD student at Programa de Doctorado en Ciencias as Universidad de Costa Rica, working at Programa de Investigacion en Enfermedades Tropicales at Universidad Nacional. Her primary research interest is the evolution of Brucella species and their host specificity.
(1.) World Health Organization. The control of neglected zoonotic diseases. In: NZD4 Organising Committee, editor. WHO conference report. Geneva: WHO Press; 2014. p. 23-35.
(2.) Stoenner HG, Lackman DB. A new species of Brucella isolated from the desert wood rat, Neotoma lepida Thomas. Am J Vet Res. 1957; 18:947-51.
(3.) Moreno E. Retrospective and prospective perspectives on zoonotic brucellosis. Front Microbiol. 2014; 5:213. http://dx.doi.org/10.3389/ fmicb.2014.00213
(4.) Stoenner HG. The behavior of Brucella neotomae and Brucella suis in reciprocal superinfection experiments in mice and guinea pigs. Am J Vet Res. 1963; 24:376-80.
(5.) Gibby IW, Gibby AM. Host-parasite relationships with Brucella neotomae. J Bacteriol. 1965; 89:9-16.
(6.) Ceran N, Turkoglu R, Erdem I, Inan A, Engin D, Tireli H, et al. Neurobrucellosis: clinical, diagnostic, therapeutic features and outcome. Unusual clinical presentations in an endemic region. Braz J Infect Dis. 2011; 15:52-9. http://dx.doi.org/10.1016/ S1413-8670(11)70140-4
(7.) Alton GG, Jones LM, Angus RD, Verger JM. Techniques for the brucellosis laboratory. Paris: Institut National de la Recherche Agronomique; 1988. p. 39-41.
(8.) Weyant R, Popovic T, Bragg S. Basic laboratory protocols for the presumptive identification of Brucella species. Atlanta: Centers for Disease Control and Prevention; 2001. p. 1-16.
(9.) Winchell JM, Wolff BJ, Tiller R, Bowen MD, Hoffmaster AR. Rapid identification and discrimination of Brucella isolates by use of real-time PCR and high-resolution melt analysis. J Clin Microbiol. 2010; 48:697-702. http://dx.doi.org/10.1128/ JCM.02021-09
(10.) Garci'a-Yoldi D, Marin CM, de Miguel M-J, Munoz PM, Vizmanos JL, Lopez-Goni I. Multiplex PCR assay for the identification and differentiation of all Brucella species and the vaccine strains Brucella abortus S19 and RB51 and Brucella melitensis Revl. Clin Chem. 2006; 52:779-81. http://dx.doi.org/ 10.1373/clinchem.2005.062596
(11.) Scott JC, Koylass MS, Stubberfield MR, Whatmore AM. Multiplex assay based on single-nucleotide polymorphisms for rapid identification of Brucella isolates at the species level. Appl Environ Microbiol. 2007; 73:7331-7. http://dx.doi.org/10.1128/ AEM.00976-07
(12.) Isidoro-Ayza M, Ruiz-Villalobos N, Perez L, Guzman-Verri C, Munoz PM, Alegre F, et al. Brucella ceti infection in dolphins from the western Mediterranean Sea. BMC Vet Res. 2014; 10:206. http://dx.doi.org/10.1186/s12917-014-0206-7
(13.) Page AJ, Taylor B, Delaney AJ, Soares J, Seemann T, Keane JA, et al. SNP-sites: rapid efficient extraction of SNPs from multi-FASTA alignments. Microb Genomics. 2016 Apr 29. http://dx.doi.org/10.1099/mgen.0.000056
(14.) Mancilla M. The Brucella genomic islands. In: Lopez-Goni I, O'Callaghan D, editors. Brucella: Molecular Microbiology and Genomics. Poole (UK): Caister Academic Press; 2012. p. 36-57.
(15.) Page AJ, De Silva N, Hunt M, Quail MA, Parkhill J, Harris SR, et al. Robust high-throughput prokaryote de novo assembly and improvement pipeline for Illumina data. Microb Genomics. Microbiology Society; 2016 Aug 25. http://dx.doi.org/10.1099/ mgen.0.000083
Marcela Suarez-Esquivel, Nazareth Ruiz-Villalobos, Cesar Jimenez-Rojas, Elfas Barquero-Calvo, Carlos Chacon-Dfaz, Eunice Vfquez-Ruiz, Norman Rojas-Campos, Kate S. Baker,  Gerardo Oviedo-Sanchez, Ernesto Amuy, Esteban Chaves-Olarte, Nicholas R. Thomson, Edgardo Moreno, Caterina Guzman-Verri
Author affiliations: Universidad Nacional, Heredia, Costa Rica (M. Suarez-Esquivel, N. Ruiz-Villalobos, C. Jimenez-Rojas, E. Barquero-Calvo, E. Vlquez-Ruiz, E. Moreno, C. Guzman-Verri); Universidad de Costa Rica, San Jose, Costa Rica (E. Barquero-Calvo, C. Chacon-DIaz, N. Rojas-Campos, G. Oviedo-Sanchez, E. Chaves-Olarte, E. Moreno, C. Guzman-Verri); Wellcome Trust Sanger Institute, Hinxton, UK (K.S. Baker, N.R. Thomson); Caja Costarricense del Seguro Social, Puntarenas, Costa Rica (E. Amuy)
 Current affiliation: University of Liverpool, Liverpool, UK.
Address for correspondence: Caterina Guzman-Verri, Programa de Investigation en Enfermedades Tropicales, Escuela de Medicina Veterinaria, Campus Benjamin Nunez, Universidad Nacional, Heredia Aptdo 304, 3000 Costa Rica; email email@example.com
Caption: Figure 1. Dendrogram based on multiple locus variable number of tandem repeats-16-loci panel analysis of Brucella spp. (http://mlva.u-psud.fr/brucella/) and clinical isolates from human cerebrospinal fluid samples from 2 patients with brucellosis. The isolates bneohCRI and bneohCR2 (red branches) showed a pattern consistent with previously reported profiles for Brucella neotomae (blue shading). Black, gray, and tic marks are used to differentiate between adjacent species. Arrows separate small neighboring clusters and indicate the B. neotomae cluster.
Caption: Figure 2. Phylogenetic tree based on 34,307 single-nucleotide polymorphisms (SNPs) found among 51 Brucella genome sequences. The clinical isolates bneohCRi and bneohCR2 cluster with B. neotomae 5K33 and differ by 164 SNPs. A different color is used to represent each Brucella species. Dotted red lines denote the 3 B. neotomae isolates, which overlap at the tip of the branch because of the high identity among them.
Table. Differential biochemical profile of Brucella isolates from 2 men with neurobrucellosis, Costa Rica, 2008 and 2011 Analysis bneohCRI bneohCR2 Biochemical tests Oxidase - - Citrate utilization - - Nitrate reduction + + C[O.sub.2] required - - [H.sub.2]S production + + Urease activity, h 0-0.5 + 0-0.5+ Growth in presence of dyes Thionin 20 [micro]g/mL - - Basic fuchsin 20 [micro]g/mL - - A + + M - -
|Printer friendly Cite/link Email Feedback|
|Author:||Suarez-Esquivel, Marcela; Ruiz-Villalobos, Nazareth; Jimenez-Rojas, Cesar; Barquero-Calvo, Elfas; Ch|
|Publication:||Emerging Infectious Diseases|
|Date:||Jun 1, 2017|
|Previous Article:||Congenital malformations of calves infected with Shamonda virus, Southern Japan.|
|Next Article:||Isolated case of Marburg virus disease, Kampala, Uganda, 2014.|