# New insights into phylogeography of worldwide Brucella canis isolates by comparative genomics-based approaches: focus on Brazil

Acacia Ferreira Vicente, Guillaume Girault, Yannick Corde, Mateus Souza Ribeiro Mioni, Lara Borges Keid, Maryne Jaÿ, Jane Megid, Virginie Mick

<figure><img src="image_1.png" /></figure>

## ► To cite this version:

Acacia Ferreira Vicente, Guillaume Girault, Yannick Corde, Mateus Souza Ribeiro Mioni, Lara Borges Keid, et al.. New insights into phylogeography of worldwide Brucella canis isolates by comparative genomics-based approaches: focus on Brazil. BMC Genomics, 2018, 19, 10.1186/s12864-018-5001-6. hal-02623500

HAL Id: hal-02623500

https://hal.inrae.fr/hal-02623500v1

Submitted on 26 May 2020

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>
RESEARCH ARTICLE

Open Access

<figure><img src="image_3.png" /></figure>

CrossMark

# New insights into phylogeography of worldwide *Brucella canis* isolates by comparative genomics-based approaches: focus on Brazil

Acácia Ferreira Vicente<fnref n="1" />,²†, Guillaume Girault<fnref n="1" />*†, Yannick Corde<fnref n="1" />,⁴, Mateus Souza Ribeiro Mioni², Lara Borges Keid³, Maryne Jay<fnref n="1" />,⁵, Jane Megid² and Virginie Mick<fnref n="1" />*

**Abstract**

**Background:** Canine brucellosis, due to *Brucella canis*, is a worldwide zoonosis that remains endemic in South America, including Brazil. Implementation of powerful whole-genome sequencing approaches allowed exploring the *Brucella* genus considered as monomorphic, with, to date, more than 500 genomes available in public databases. Nevertheless, with under-representation of *B. canis* genomes –only twenty complete or draft genomes–, lack of knowledge about this species is still considerable. This report describes a comparative genomics-based phylogeographic investigation of 53 *B. canis* strains, including 28 isolates paired-end sequenced in this work.

**Results:** Obtained results allow identifying a SNP panel species-specific to *B. canis* of 1086 nucleotides. In addition, high-resolution analyses assess the epidemiological relationship between worldwide isolates. Our findings show worldwide strains are distributed among 2 distinct lineages. One of them seems to be specific to South American strains, including Brazil. *B. canis* South American strains may be identified by a SNP panel of 15 nucleotides, whereas a 22 SNP panel is sufficient to define contamination origin from Brazil. These results lead to the proposal of a possible spread route for dog brucellosis through South America. Additionally, whole-genome analyses highlight the remarkable genomic stability of *B. canis* strains over time and the sustainability of the infection in São Paulo over 12 year-period.

**Conclusions:** Significant increase of *B. canis* genomes available in public databases provides new insights into *B. canis* infection in South America, including Brazil, as well as in the world, and also offers new perspectives for the *Brucella* genus *largo* sensu.

**Keywords:** *Brucella canis*, Dog brucellosis, Brazil, São Paulo, Europe, Whole genome sequencing, SNP, Phylogeography

**Background**

Canine brucellosis is a worldwide zoonosis caused by *Brucella canis* [1–3]. This bacterium is usually associated with dogs and occasionally causes brucellosis in humans [1]. It is considered as the main cause of reproductive failure in

dogs, and responsible for important economic losses to kennels [3, 4].

*B. canis* infections, mainly diagnosed by serological methods and bacteriological evidence [5], remain endemic in South America, including Brazil with high dog population [3], raising public and animal health concerns. Although some studies have emphasized clinical and epidemiological importance of *B. canis* in the canine population, especially in confinement areas, such as kennels that increase the probability of disease transmission [6], canine brucellosis is underestimated in animals and

<figure><img src="image_4.png" /></figure>

<footnote>* Correspondence: guillaume.girault@anses.fr; virginie.mick@anses.fr</footnote>

<footnote>†Acácia Ferreira Vicente and Guillaume Girault contributed equally to this work.</footnote>

<footnote n="1">¹EU/OIE/FAO & National Reference Laboratory for animal Brucellosis, Animal Health Laboratory, Paris-Est University/Anses, Maisons-Alfort, France
Full list of author information is available at the end of the article</footnote>
human and its epidemiological aspects are still poorly understood.

*Brucella* is a well-known genus, whose taxonomy is in constant evolution, with the ongoing description of new species [7–9]. To date, *Brucella* genomes are composed of 2 circular chromosomes, except for *B. suis* biovar 3 (only 1 chromosome), with a global genome size of approximately 3.3 Mb [10]. All *Brucella* species are highly related to each other genetically, with sequence similarity values of 98% to 100% in the core genome [11]. Despite its genetic homogeneity, genotyping and phylogenetic approaches based on multiple genomic markers, like MLST –Multi Locus Sequence Typing –[12] or MLVA –Multiple-Locus Variable number tandem repeat Analysis–[13], are robust tools for wide- and fine-scale epidemiological/taxonomical investigations. Moreover, high-resolution studies based on *Brucella* whole genome strategies have been reported most recently, providing new insights into the genus [14–19]. Nevertheless, lack of knowledge about *B. canis* is still considerable. Indeed, although more than 500 *Brucella* genomes are available on public databases, *B. canis* is poorly represented with only 20 complete and draft genomes available [20–25].

The purpose of this work was a comparative genomics-based phylogeographic investigation of 27 worldwide *B. canis* field strains, with a focus on Brazil, in order to improve knowledge about *B. canis* epidemiology in the world, and especially in South America.

# Methods

## *Brucella* strains

Twenty-eight dog and human *B. canis* strains –including the *B. canis* reference strain Rm6–66 (ATCC 23365)–were isolated from routine epidemiological veterinary investigations conducted in São Paulo, Brazil, in 2005 and 2015 (Keid, L.B., personal data) or were obtained from the ANSES collection, especially from Europe (Additional file 1: Table S1).

## Phenotypic characterization

All isolates were confirmed as *B. canis* using conventional *Brucella* typing methods, based on CO₂ requirement, H₂S production, oxidase test, urea hydrolysis, agglutination with monospecific sera, fuchsin, thionin and safranin dye sensitivity and phage typing [26].

## Molecular studies

Genomic DNA was extracted from *Brucella* cultures using the High Pure PCR Template Preparation Kit (Roche Diagnostics, France), according to the manufacturer's instructions. Molecular species confirmation was performed by Suis-Ladder multiplex PCR assay as previously described [27].

In this study, 28 *B. canis* strains (Additional file 1: Table S1), including the *B. canis* reference strain Rm6–66 (ATCC 23365), were paired-end sequenced at the Genoscreen institute (Lille, France) on Illumina platforms –HiSeq2000 and MiSeq (2 × 250 bp)–. An average 50-fold sequencing depth was obtained (100-fold coverage for Rm6–66). In addition, complete and draft genomes, as well as Sequence Read Archives (SRA, raw reads), belonging to the *B. canis* species (n = 25) were retrieved from Pathosystems Resource Integration Center (PATRIC) and NCBI center (Additional file 1: Table S1).

To discard misassemblies, raw sequencing reads were trimmed using Trimmomatic-0.36 (phred33, minimum 50 bp length). Chimeric genomes from databases have been generated by merging both chromosomes 1 and 2 to compare complete and draft genomes. To homogenize data, sets of sequencing pseudo-reads were created in silico on contigs from public databases using ART program with a 250 bp length and a 50-fold depth [28]. The mapping step of short read datasets on FASTQ format from this study as well as from pseudo-reads generated with ART was realized using BWA algorithm implemented in BioNumerics 7.6.1 (Applied Maths, Belgium) with 90% sequence similarity against the Rm6–66 *B. canis* reference genome (CP007758.1, CP007759.1). A set of SNPs was determined for each genome sequencing data respect to the Rm6–66 genome, using the BioNumerics whole genome SNP (wgSNP) module. Several position filters were applied on the SNP matrix: (i) contiguous SNPs were removed (if found in a 10 bp-window), (ii) position mask for repeated elements, including tandem repeats, (iii) a minimum 20-fold coverage for each SNP was required, and (iv) ambiguous (i.e. non-ACGT bases) and unreliable bases (i.e. Ns) were discarded. The refined SNP matrix was used to generate a minimum spanning tree using maximum bootstrap maximum likelihood approach, allowing phylogenetic analyses.

## Nucleotide sequence accession numbers

Raw sequencing data of five representative strains of this whole-genome sequencing project (PRJEB22763) have been deposited in the European Nucleotide Archive (ENA) (http://www.ebi.ac.uk/ena/data/view/PRJEB22763) under accession numbers ERR2136545, ERR2136546, ERR2136547, ERR2136548 and ERR2136549 (Additional file 1: Table S1).

# Results

All isolates investigated in this study were confirmed as *B. canis* by phenotypic and molecular (Suis-Ladder) approaches (data not shown). A total of 27 *B. canis* field strains, as well as the *B. canis* reference strain Rm6–66, were characterized by paired-end WGS. After applying different filters, a maximum likelihood tree (Fig. 1) was
generated from 53 *B. canis* complete genomes and rooted with the reference strain *B. melitensis* biovar 1 16 M. A 0.25% homoplasia was considered. Bootstrap values supported strong confidence for each branch (mean: 97.43%; intervals 63–100%). A total of 7327 chromosomal SNPs were identified, including 1086 nucleotides specific to *B. canis*.

The *B. canis* investigated strains were distributed into two distinct clades (Fig. 1). The first clade, named here lineage 1 (n = 26), was mainly composed of strains isolated in Europe, Asia and USA, including the reference strain. Lineage 1 diverged from the most recent common ancestor –MRCA–of all *B. canis* isolates with only 3 SNPs. Two branches were then identified and diverged from the node with similar SNP numbers (17 and 21 SNPs). The second clade, named here lineage 2 (n = 27),

contained all the strains isolated from South America (Argentina, Brazil, Chile, Colombia), as well as 2 North and Central American strains (USA, Mexico) and only one from Europe (Sweden).

Asia

Four Asian strains were represented in this study, isolated from China (n = 2), Japan (n = 1) and South Korea (n = 1). They clustered together into a distinct sub-branch in lineage 1 (Fig. 1), among a polytomy of 3 branches, including too a sub-branch with European isolates (Finland, Sweden, Serbia, and Germany) and one American isolate (USA), and a sub-branch in singleton (South Africa). These four strains analyzed in this study seemed to be defined by a specific SNP panel of 12 markers.

<figure><img src="image_5.png" /></figure>
## Europe

All strains isolated from Europe included in this study (n = 11), except for one, were clustered in lineage 1 (Fig. 1). Curiously, only one strain from Sweden -SVA13- was identified in lineage 2, 87 SNP distant from MRCA. This strain has been isolated from a dog imported from Spain for breeding [23]. Interestingly, one other Spanish strain -A592-was included in this comparative analysis. It clustered in lineage 1, together with French isolates, reference strains and USA strains, respectively with a distance of 69, 70 and 73 SNPs. Similarly, the *B. canis* German strain E291 was close to a Serbian strain -E286-with 23 SNP difference, whereas a 105 SNP difference was observed with another Serbian strain 04-2330-1 (strain from ANSES collection previously sequenced by Broad Institute Center). Finnish strains exhibited too some differences, with a 78 SNP distance from each other. These results might suggest the existence of regional polymorphism or the circulation of different *B. canis* clones among a same country.

## North America

Seven strains from North America (USA) were included in this study. Surprisingly, the great majority (n = 6) clustered into lineage 1 (Fig. 1), and only one strain -2,009,013,648-[25] isolated from Arizona in 2009 belonged to lineage 2, including all South and Central American strains (Fig. 1). Interestingly, Arizona shares common boundaries with Mexico. Regarding closeness of both strains, results suggest a same cross-border contamination origin.

Within lineage 1, the distribution of USA strains was heterogeneous since strains belonged to distinct sub-clusters -e.g. strains 2,010,009,751 and 2,009,004,498 respectively isolated from Massachusetts in 2010 and from Louisiana in 2009 [25] were 140 SNP distant between each other-. Although almost 40 *B. canis* raw sequencing data have recently been added to NCBI SRA (PRJNA369091), faced with lack of information and metadata poorly filled, only 4 of these strains have been included in this study. They showed weak polymorphism (distance between furthest strains = 32 SNPs) and clustered together with *B. canis* reference genomes.

Despite possible sampling bias, WgSNP results seemed to indicate an important genetic diversity of *B. canis* strains circulating in USA.

## Central and South America

This study describes a comparative genomics-based investigation of Brazilian *B. canis* strains. To date, lack of

knowledge regarding *B. canis* epidemiology in South America, and especially in Brazil, is reported.

All *B. canis* strains isolated from South and Central America were located among lineage 2 (Fig. 1). Lineage 2 split into two distinct branches. The minor one contained the Mexican isolate Mex51, clustered together with an USA isolate and a Swedish strain imported from Spain. The major branch allowed grouping all South American strains into a unique, exclusive and homogeneous cluster, with a node that could represent the MRCA of South American *B. canis* strains -SA-MRCA -. The two most distant strains among South America isolates harbored a 87 SNP difference (Fig. 2).

This South American branch split in two sub-branches. A first sub-branch contained strains isolated from Colombia -Oliveri and CNGB 1172-and radiated respectively 37 and 38 SNPs away from SA-MRCA. Both Colombian isolates were very similar, with only 3 SNP difference despite isolation from different hosts (human and dog, respectively for CNGB 1172 and for Oliveri), suggesting contamination by the same strain. A second sub-branch radiated 7 SNPs away from SA-MRCA and split in 3 sub-clusters. Except for two singletons, isolated from Chile -SCL and CNGB 513-, only Brazilian strains, as well as one Argentinian strain (CNGB 1324) constituted this sub-branch. Interestingly, the strain B003 isolated from South Brazil, adjacent to Argentina, seemed closer to the Argentina strain than other Brazilian strains, with a minimal distance of 50 SNPs, suggesting an infection contamination from Argentina to Brazil. All other Brazilian strains formed a single complex, split into 2 branches that have diverged very early (1 and 2 SNPs) and harbored very weak polymorphism (maximal distance of 25 SNPs). In addition, length of different branches identified in this study varied very little, from 3 to 14 SNPs, even over a 17 year-- period (1998-2015).

Thus, a 12 SNP panel seemed to be sufficient to characterize the *B. canis* South American strains, although *B. canis* contamination origin from Brazil might be identified using a 22 SNP panel (Additional file 2: Table S2).

## Focus on São Paulo (Brazil)

Seventeen *B. canis* isolates collected from São Paulo, Brazil, over 12-year-period (2003-2015) were investigated in this study. Fourteen were isolated from three well-identified kennels, respectively both in 2005 (k1, k2) and one in 2015 (k3), and three strains were isolated in 2003 from not reported kennels (Additional file 1: Table S1). Interestingly, no trade/historical exchange between kennels was reported.

As expected, wgSNP analyses showed a minor SNP difference, independent of isolation year, between São
<figure><img src="image_6.png" /><figcaption>**Fig. 2** SNP distance between different *B. canis* strains isolated from South America. This figure indicates the distance observed in SNPs between the strains isolated from South America in 4 different countries. A: Argentina (dashed lines); B: Brazil (full lines); C1: Colombia (dotted lines); C2: Chile (crossed lines). [Blank map available on Wikimedia commons]</figcaption></figure>

Paulo strains, ranging from 3 to 25 nucleotides (Fig. 1). Most *B. canis* strains isolated among a same kennel were strictly identical or harbored no significant difference, as strains collected in the same isolation year (e.g. E087, E122 from k1 and E143, E243, E246, E257 from k2).

## Discussion

To improve understanding of the genetic diversity of *B. canis*, we genotyped 27 worldwide *B. canis* field isolates by performing whole genome SNP analysis. Indeed, this approach is able to infer the relationships among worldwide *B. canis* isolates.

Comparative genomic analysis of 53 dog and human *B. canis* strains was performed. Robustness of the wgSNP method applied in this study was assessed by sequencing data comparison of 2 isolates from ANSES collection (Additional file 1: Table S1), independently sequenced in this study and previously by the Broad Institute center. Thus, C280 and SVA10 strains showed a strictly identical genome, as well as A590 and 96–7258 genomes harbored only one SNP difference. These results highlighted in vitro genomic stability of *B. canis* genomes despite successive cultural steps susceptible to induce several mutational events. Similarly and

unsurprisingly, wgSNP results showed too extreme in vivo genomic stability over time, e.g. two samples isolated in China (BCB018 and 118) [20, 21] with only 12 SNP difference in a 20 year-interval.

No host specificity was observed. Nevertheless, SNP analysis indicated a spatial distribution of the isolates that could not be correlated with a time-span scale. A slight, almost inexistent genetic diversity was observed among Brazilian isolates. According to the more parsimonious hypothesis, our results suggested one single introduction could have led to the divergence of two Brazilian sub-branches observed here. Moreover, the short branch length, even over a long time-period, suggested circulation of a dominant clone in South America and emphasizing genetic stability of *B. canis* genomes over time. It is interesting to note *B. canis* South American strains could be identified by a SNP panel of 12 nucleotides, whereas a 22 SNP panel is sufficient to refine Brazilian contamination origin of *B. canis* strains, assuming contamination origin of Brazilian B003 strain is from Argentina.

Here, we investigated potential genetic diversity over 12-year-period in a restricted geographical area: from different São Paulo kennels. Absence of significant
difference whatever the isolation year between *B. canis* studied isolates suggested a same contamination origin and/or the circulation of a dominant clone. In addition, the distance observed over a 10-year period between strains isolated from a same location was not greater (e.g. E267 vs E258, or E267 vs B009), suggesting the sustainability of *B. canis* infection over time in São Paulo.

The potential of whole-genome sequencing leads to propose a possible spread route for dog brucellosis through South America. Indeed, our findings suggested the presence of *B. canis* in South America probably might be resulted from an introduction from USA and/or from Mexico (Fig. 1) to Colombia (maximum distance of 84, 85 and 87 SNPs respectively with Chile, Brazil and Argentina), followed by spread to Brazil, Argentina and Chile (Fig. 2). On the current state of available genome dataset –Argentina, Chile and Brazil are part of a polytomic branch (Fig. 1)–, it is difficult to hypothesize into the details this second evolving step of the *B. canis* infection through South America. Indeed, a largest dataset of representative strains of these countries is required to propose a strongest hypothesis.

## Conclusions

This report describes a comparative genomics-based phylogeographic investigation of *B. canis* field strains isolated from Brazil. Results obtained allow assessing the epidemiological relationship between worldwide strains and hypothesizing a possible spread route for dog brucellosis through South America. In addition, whole genome analyses highlight the remarkable genomic stability of *B. canis* strains over time and the sustainability of the infection in São Paulo.

Significant increase of *B. canis* genomes available in public databases, resulting from this work, provides new insights not only into *B. canis* infection in South America, including Brazil, as well in the world, but also offers new perspectives for the *Brucella* genus *largo* sensu.

## Additional files

**Additional file 1: Table S1.** *B. canis* strains investigated in this study and from public databases and whole sequencing data. (DOCX 19 kb)

**Additional file 2 Table S2.** Determination of 22 SNP panel specific to Brazilian *B. canis* strains. (DOCX 37 kb)

### Abbreviations

MLST: Multi Locus Sequence Typing; MLVA: Multiple-Locus Variable number tandem repeat Analysis; MRCA: Most Recent Common Ancestor; SNP: Single Nucleotide Polymorphism; WGS: Whole Genome Sequencing; wgSNP: Whole genome SNP

### Acknowledgements

The authors wish to express their gratitude to Marion Holzapfel, Bruna Letícia Devidé Ribeiro and Fernando José Paganini Listoni for the excellent technical assistance.

### Funding

The PhD researcher was supported by FAPESP (2014/00889-0) and CAPES (88887.117958/2016-00) programs. Research work was supported by the UE project Bru-EPIDIA (291815-FP7/ERANET/ANIHWA). The funding agencies were not involved in the design of the study, collection, analysis, and interpretation of data and in writing the manuscript.

### Availability of data and materials

Most of the data generated during this study are included in this article and supplementary information files and available in the NCBI Genome database (https://www.ncbi.nlm.nih.gov/genome/?term=Brucella+canis) or PATRIC database (https://www.patricbrc.org/view/Taxonomy/234#view_tab=genomes).

Datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.

### Authors' contributions

Conceived the study: JM, VM; Isolation: AFV, MSRM, LBK; Shipment support: MJ; Phenotypic and molecular identification: AFV, YC, MSRM, LBK; performed DNA extractions: AFV, GG, YC; performed whole genome sequencing analysis: AFV, GG; supervised the laboratory work JM, VM; wrote the paper, AFV, GG, V.M. All authors participated in the presentation and discussion of results. All authors read and approved the final manuscript.

### Ethics approval and consent to participate

Ethical approval (CEUA 0033/2017) was required to conduct epidemiological veterinary investigations in São Paulo, Brazil, and to acquire the samples analyzed in this study.

### Consent for publication

Not applicable.

### Competing interests

The authors declare that they have no competing interests.

## Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

### Author details

¹EU/OIE/FAO & National Reference Laboratory for animal Brucellosis, Animal Health Laboratory, Paris-Est University/Anses, Maisons-Alfort, France.

²Molecular Biology Laboratory, Dept. Veterinary Hygiene and Public Health, FMVZ, UNESP, Botucatu, Brazil. ³Dept. Veterinary Medicine, University of Animal Science and Food Engineering, USP, Pirassununga, Brazil. ⁴Present Address: Physiology of reproduction and behaviour joint research unit, INRA Val de Loire Centre, Nouzilly, France. ⁵Present Address: Anses, Laboratoire de Lyon, UMR Mycoplasmoses des Ruminants, Lyon, France Université de Lyon, VetAgro Sup, UMR Mycoplasmoses des Ruminants, Marcy L'Etoile, France.

Received: 24 January 2018 Accepted: 8 August 2018

Published online: 28 August 2018

### References

1. Brucellosis in humans and animals. WHO guidance. Geneva: World Health Organization (WHO); 2006.

2. Hofer E, Bagó Z, Revilla-Fernández S, Melzer F, Tomaso H, López-Goñi I, Fasching G, Schmoll F. First detection of *Brucella canis* infections in a breeding kennel in Austria. New Microbiol. 2012;35:507-10.

3. Keid LB, Chiebao DP, Batinga MC, Faita T, Diniz JA, Oliveira TM, Ferreira HL, Soares RM. *Brucella canis* infection in dogs from commercial breeding kennels in Brazil. Transbound Emerg Dis. 2017;64:691-7. https://doi.org/10.1111/tbed.12632.

4. Gyuranecz M, Szeredi L, Ronai Z, Dénes B, Dencso L, Dán A, Pálmai N, Hauser Z, Lami E, Makrai L, Erdélyi K, Jánosi S. Detection of *Brucella canis*-induced reproductive diseases in a kennel. J Vet Diagn Investig. 2011;23:143-7. https://doi.org/10.1177/104063871102300127.

5. Manual of diagnostic tests and vaccines for terrestrial animals. World Organization for Animal Health (OIE). 2016. http://www.oie.int/en/standard-setting/terrestrial-manual/access-online/. Accessed 2016.
6. Hollett RB. Canine brucellosis: outbreaks and compliance. Theriogenology. 2006;66:575-87. https://doi.org/10.1016/j.theriogenology.2006.04.011.

7. Scholz HC, Revilla-Fernandez S, Al Dahouk S, Hammerl JA, Zygmunt MS, Cloeckaert A, Koylass M, Whatmore AM, Blom J, Vergnaud G, Witte A, Aistleitner K, Hofer E. *Brucella vulpis* sp. nov. isolated from mandibular lymph nodes of red foxes (*Vulpes vulpes*). Int J Syst Evol Microbiol. 2016;66:2090-8. https://doi.org/10.1099/ijsem.0.000998.

8. Soler-Llorens PF, Quance CR, Lawhon SD, Stuber TP, Edwards JF, Ficht TA, Robbe-Austerman S, O'Callaghan D, Keriel A. *A Brucella* spp. isolate from a Pac-Man Frog (*Ceratophrys omata*) reveals characteristics departing from classical *Brucella*. Front Cell Infect Microbiol. 2016;6:116. https://doi.org/10.3389/fcimb.2016.00116.

9. Whatmore AM, Davison N, Cloeckaert A, Al Dahouk S, Zygmunt MS, Brew SD, Perrett LL, Koylass MS, Vergnaud G, Quance C, Scholz HC, Dick EJ Jr, Hubbard G, Schlabritz-Loutsevitch NE. *Brucella papionis* sp. nov., isolated from baboons (*Papio spp*.). Int J Syst Evol Microbiol. 2014;64:4120-8. https://doi.org/10.1099/ijs.0.065482-0.

10. Halling SM, Peterson-Burch BD, Bricker BJ, Zuerner RL, Qing Z, Li LL, Kapur V, Alt DP, Olsen SC. Completion of the genome sequence of *Brucella abortus* and comparison to the highly similar genomes of *Brucella melitensis* and *Brucella suis*. J Bacteriol. 2005;187:2715-26. https://doi.org/10.1128/JB.187.8.2715-2726.2005.

11. Scholz HC, Vergnaud G. Molecular characterisation of *Brucella* species. Rev Sci Tech. 2013;32:149-62.

12. Whatmore AM, Perrett LL, MacMillan AP. Characterisation of the genetic diversity of *Brucella* by multilocus sequencing. BMC Microbiol. 2007;32:149-62. https://doi.org/10.1186/1471-2180-7-34.

13. Le Fleche P, Jacques I, Grayon M, Al Dahouk S, Bouchon P, Denoeud F, Nöckler K, Neubauer H, Guilloteau LA, Vergnaud G. Evaluation and selection of tandem repeat loci for a Brucella MLVA typing assay. BMC Microbiol. 2006;6:9. https://doi.org/10.1186/1471-2180-6-9.

14. Foster JT, Beckstrom-Sternberg SM, Pearson T, Beckstrom-Sternberg JS, Chain PS, Roberto FF, Hnath J, Brettin T, Keim P. Whole-genome-based phylogeny and divergence of the genus *Brucella*. J Bacteriol. 2009;191:2864-70. https://doi.org/10.1128/JB.01581-08.

15. Koylass MS, King AC, Edwards-Smallbone J, Gopaul KK, Perrett LL, Whatmore AM. Comparative performance of SNP typing and 'Bruce-ladder' in the discrimination of *Brucella suis* and *Brucella canis*. Vet Microbiol. 2010;142:450-4. https://doi.org/10.1016/j.vetmic.2009.10.013.

16. Foster JT, Price LB, Beckstrom-Sternberg SM, Pearson T, Brown WD, Kiesling DM, Allen CA, Liu CM, Beckstrom-Sternberg J, Roberto FF, Keim P. Genotyping of *Brucella* species using clade specific SNPs. BMC Microbiol. 2012;12:110. https://doi.org/10.1186/1471-2180-12-110.

17. Sanchez-Jimenez MM, Isaza JP, Alzate JF, Olivera-Angel M. Comparison of *Brucella canis* genomes isolated from different countries shows multiple variable regions. Genomics. 2015;106:43-51. https://doi.org/10.1016/j.ygeno.2015.03.004.

18. Tan KK, Tan YC, Chang LY, Lee KW, Nore SS, Yee WY, Mat Isa MN, Jafar FL, Hoh CC, AbuBakar S. Full genome SNP-based phylogenetic analysis reveals the origin and global spread of *Brucella melitensis*. BMC Genomics. 2015;16:93. https://doi.org/10.1186/s12864-015-1294-x.

19. Sankarasubramanian J, Vishnu US, Gunasekaran P, Rajendhran J. A genome-wide SNP-based phylogenetic analysis distinguishes different biovars of *Brucella suis*. Infect Genet Evol. 2016;41:213-7. https://doi.org/10.1016/j.meegid.2016.04.012.

20. Gao G, Li J, Li T, Zhang Z, Wang L, Yuan X, Wang Y, Xu J, Ke Y, Huang L, Wang D, Chen Z, Xu X. Complete genome sequence of *Brucella canis* strain 118, a strain isolated from canine. J Bacteriol. 2012;194:6680. https://doi.org/10.1128/JB.01853-12.

21. Wang Y, Ke Y, Zhen Q, Yuan X, Xu J, Qiu Y, Wang Z, Li T, Wang D, Huang L, Chen Z. Complete genome sequence of *Brucella canis* BCB018, a strain isolated from a human patient. J Bacteriol. 2012;194:6697-8. https://doi.org/10.1128/JB.01811-12.

22. Kim JS, Jeong W, Jeoung HY, Song JY, Kim H, Beak JH, Parisutham V, Lee SK, Kim JW, Kim JY, Jung SC, Her M, An DJ. Complete genome sequence of *Brucella canis* strain HSK A52141, isolated from the blood of an infected dog. J Bacteriol. 2012;194:5134. https://doi.org/10.1128/JB.01125-12.

23. Kaden R, Agren J, Ferrari S, Lindberg M, Backman S, Wahab T. Whole-genome sequence of *Brucella canis* strain SVA13, isolated from an infected dog. Genome Announc. 2014;2 https://doi.org/10.1128/genomeA.00700-14.

24. Minogue TD, Daligault HA, Davenport KW, Bishop-Lilly KA, Broomall SM, Bruce DC, Chain PS, Chertkov O, Coyne SR, Frey KG, Gibbons HS, Jaissle J, Koroleva GI, Ladner JT, Lo CC, Palacios GF, Redden CL, Rosenzweig CN, Scholz MB, Xu Y, Johnson SL. Whole-genome sequences of 24 *Brucella strains*. Genome Announc. 2014;2 https://doi.org/10.1128/genomeA.00915-14.

25. Viana MVC, Wattam AR, Batra DG, Boisvert S, Brettin TS, Frace M, Xia F, Azevedo V, Tiller R, Hoffmaster AR. Genome sequences of three *Brucella canis* strains isolated from humans and a dog. Genome Announc. 2017;5 https://doi.org/10.1128/genomeA.01688-16.

26. Alton GG, Jones LM, Angus RD, Verger JM. Techniques for the brucellosis laboratory. Paris: Institut National de la Recherche Agronomique (INRA); 1988.

27. Lopez-Goñi I, Garcia-Yoldi D, Marin CM, de Miguel MJ, Barquero-Calvo E, Guzman-Verri C, Albert D, Garin-Bastuji B. New Bruce-ladder multiplex PCR assay for the biovar typing of *Brucella suis* and the discrimination of *Brucella suis* and *Brucella canis*. Vet Microbiol. 2011;154:152-5. https://doi.org/10.1016/j.vetmic.2011.06.035.

28. Huang W, Li L, Myers JR, Marth GT. ART: a next-generation sequencing read simulator. Bioinformatics. 2012;28:593-4. https://doi.org/10.1093/bioinformatics/btr708.

Ready to submit your research? Choose BMC and benefit from:

* fast, convenient online submission

* thorough peer review by experienced researchers in your field

* rapid publication on acceptance

* support for research data, including large and complex data types

* gold Open Access which fosters wider collaboration and increased citations

* maximum visibility for your research: over 100M website views per year

At BMC, research is always in progress.

Learn more biomedcentral.com/submissions

<figure><img src="image_7.png" /></figure>