Brucella abortus Strain 2308 Wisconsin Genome: Importance of the Definition of Reference Strains

Brucellosis is a bacterial infectious disease affecting a wide range of mammals and a neglected zoonosis caused by species of the genetically homogenous genus Brucella. As in most studies on bacterial diseases, research in brucellosis is carried out by using reference strains as canonical models to understand the mechanisms underlying host pathogen interactions. We performed whole genome sequencing analysis of the reference strain B. abortus 2308 routinely used in our laboratory, including manual curated annotation accessible as an editable version through a link at https://en.wikipedia.org/wiki/Brucella#Genomics. Comparison of this genome with two publically available 2308 genomes showed significant differences, particularly indels related to insertional elements, suggesting variability related to the transposition of these elements within the same strain. Considering the outcome of high resolution genomic techniques in the bacteriology field, the conventional concept of strain definition needs to be revised.

Introduction

Brucella is a bacterial genus responsible for brucellosis, a disease in animals causing infertility, pre-term birth, or abortion (Moreno and Moriyón, 2006). It is also one of the most worldwide spread bacterial zoonosis, not only causing human suffering but also representing a significant economic burden on animal industries. Because this severe and debilitating disease has not been adequately addressed at some national and international level, it is considered by WHO (2014) as one of the “forgotten neglected zoonosis,” constituting a major burden for poor rural communities.

The main etiological agent of cattle brucellosis is B. abortus associated with abortion, infertility, reproductive failure, and decreased milk production (Harmon et al., 1989; Xavier et al., 2009; Neta et al., 2010). Humans usually become infected when in contact with infected animals or their derived products, particularly non-pasteurized dairy products (Spink, 1956). Several B. abortus reference strains have been described and used as models to understand Brucella pathogenesis or as challenge strains for vaccine testing (Meyer and Morgan, 1973; Chain et al., 2005; Halling et al., 2005; Singh et al., 2015).

Genetic drift causing loss of virulence or antigenic properties has been reported in several Brucella vaccines strains (Bosseray, 1991; Mukherjee et al., 2005; Jacob et al., 2006; Miranda et al., 2013), but very little is known on the genetic stability of reference strains.

Brucella abortus strain 2308 was originally described as a highly virulent strain recovered in 1940 from an aborted fetus of a cow which had been in contact with cattle experimentally infected with a mixture of B. abortus cultures (Jones et al., 1965). Since then it has been widely used as a reference and challenge strain within the brucellosis research community (Trant et al., 2010; Dabral et al., 2014; Miranda et al., 2015; Olsen et al., 2015; Truong et al., 2016; Zhu et al., 2016). Whole Genome Sequencing (WGS) of this strain was first carried out by Chain et al. (2005), who found conservation in chromosome synteny with other sequenced B. abortus genomes and also some differences that were suggested to be strain specific. In 2009, the results of a WGS analysis from a strain named B. abortus 2308A were publically available under accession number GCA_000182625.1. No additional information regarding the isolates used in these two WGS projects was provided.

We performed WGS, functional annotation and manual curation of reference strain B. abortus 2308 kept at the Tropical Disease Research Program, Veterinary School, National University in Costa Rica and compared it with these two published B. abortus 2308 genomes. Significant differences were found among the genomes, challenging the idea of reference strains as a non-changing entities in time and among laboratories. The use and communication of standardized quality control and experimental design protocols could help interpretation and follow up of reported results.

Materials and Methods

Strain Description

All procedures involving live B. abortus were carried out according to the “Reglamento de Bioseguridad de la CCSS 39975-0,” year 2012, after the “Decreto Ejecutivo #30965-S,” year 2002 and research protocol approved by SIA 0434-14 from the National University, Costa Rica. A vial of B. abortus strain 2308 was obtained at the Tropical Disease Research Program, Veterinary School, National University, Costa Rica, from Dr. Ignacio Moriyón at University of Navarra, Pamplona, Spain. He received it in 1983 from Dr. Lois M. Jones from the laboratory of Prof. David T. Berman at University of Wisconsin Madison as a lyophilized vial coming originally from the National Animal Disease Center at Ames, IA, USA. From this, bacteria were expanded in trypticase soy agar, assayed biochemically to assure properties (Alton et al., 1988) and stored at -70°C in 20% glycerol. For master seed preparation, recommendations for identification, maintenance and to rule out attenuation were followed as suggested (Alton et al., 1988). Therefore, before master seed storage, a passage in mice was performed using 10⁵CFU from an overnight culture in 100 μL PBS for intraperitoneal inoculation. Bacteria were recovered from the spleen after 3 weeks after platting in trypticase soy agar and a single colony expanded in trypticase soy broth overnight (Bosseray, 1991; Grilló et al., 2012). Trypticase soy broth aliquots were prepared with 20% glycerol and stored at -70°C as a master seed. Protocols for experimentation with mice were revised and approved by the Comité Institucional para el Cuido y Uso de los Animales of the Universidad de Costa Rica (CICUA-47-12) and were in agreement with the corresponding law, Ley de Bienestar de los Animales, of Costa Rica (law 7451 on animal welfare). Mice were housed in the animal building of the Veterinary School, Universidad Nacional, Costa Rica. Animals were kept in cages with water and food ad libitum under biosafety containment conditions previous to and during the experiment.

Finally, a master seed aliquot was grown in tryptic soy agar for phenotypical and biochemical characterization (Alton et al., 1988), MLVA-16 (Le Flèche et al., 2006), Bruce-ladder PCR (López-Goñi et al., 2011) and infection of Raw macrophages as described (Palacios-Chaves et al., 2011). Aliquots were used only once and discharged.

Due to its origin, this strain is from herein referred as B. abortus strain 2308 Wisconsin (B. abortus 2308W).

WGS, Assembly and Annotation

A bacteriological loop sample from a frozen master seed aliquot of B. abortus 2308W was inoculated in trypticase soy broth and grown overnight. DNA was extracted using the Promega Wizard Genomic DNA Purification kit, and stored at -70°C until used.

WGS was performed at the Wellcome Trust Sanger Institute on Illumina platforms according to in house protocols (Quail et al., 2009, 2012), resulting in 4 396 650 reads, depth of coverage: 128.76 ± 34.84, an error rate of 0.005 and a duplication rate of 0.0017.

The SNPs discovery performed in this study used raw reads with a base quality score of Q33. Sequence variation was called by BCFtools¹ if the depth of coverage is greater than 5, the variant is present in at least 75% of reads at that position, the variant is present on both strands, and the mapping quality is greater than 30. Each SNP was manually investigated using BamView (Carver et al., 2010).

For WGS assembly and alignment, sequencing reads were de novo assembled using Velvet Optimiser (Zerbino and Birney, 2008; Page et al., 2016) and 18 contigs (N50 = 294779) were ordered using Abacas (Assefa et al., 2009) against B. abortus 9–941 (accession numbers NC_006932 and NC_006933). To detect mis-assemblies, raw data were mapped back against the 2308W genome assembly using SMALT 0.5.8² which resulted in 99.89% mapping. Two more assembly controls were performed by mapping the 2308W raw reads against 2308 reference genome (NC_007618 and NC_007624) or 9–941 reference genome. This resulted in 98.88 and 98.77% mapping, respectively.

Annotation was automatically transferred from a previously WGS obtained from B. abortus 2308 (NC_007618 and NC_007624). In order to facilitate the 2308W annotation review, a BLAST comparison (Altschul et al., 1997) between 2308W chromosomes, B. abortus 9–941, B. melitensis 16M (accession NC_003317.1 and NC_003318.1) and B. suis (NC_004310.3 and NC_004311.2) was performed. Visualizations were done with Artemis and BLAST comparisons with the Artemis Comparison Tool (Rutherford et al., 2000). Each coding region (CDS) was checked manually and curation was performed according to available experimental information and literature search on B. abortus 2308. Major findings were summarized in Supplementary Table S1 and a link to an editable version is available at https://en.wikipedia.org/wiki/Brucella#Genomics.

Sequencing info has been deposited at the European Nucleotide Archive (ENA)³ under the accession code ERS568782.

Genome Comparison among B. abortus 2308 Strains

A BLAST comparison using the genome here described, B. abortus 2308W (assembly accession ERS568782) and two previously published genomes: B. abortus 2308 (Chain et al., 2005; accession number NC_007618 and NC_007624); and B. abortus str. 2308A (assembly accession GCA_000182625.1; contigs accession numbers ACOR01000001-ACOR01000009) was performed. Visualizations were done with Artemis and ACT (Rutherford et al., 2000).

As the 2308A genome is reported as contigs, raw sequencing data were downloaded and de novo assembly using Velvet Optimizer (Zerbino and Birney, 2008) was performed followed by ordering with Abacas (Assefa et al., 2009), using B. abortus 9–941 genome as reference as per the analysis of our strain.

Results and Discussion

For comparative analysis, we performed WGS analysis of the 2308 strain that has been used in our laboratory and referred here as B. abortus 2308W, to distinguish its genome from previously reported ones (2308 and 2308A). The genome consists of two chromosomes, chromosome I is 2.10 kb in size and chromosome II is 1.16 kb. Automated annotation and manual curation were performed and summarized in Supplementary Table S1. An editable spreadsheet was also created in order to facilitate updates in the annotation or additional relevant comments (e.g., virulence, function, mutants availability) and is available through a link at https://en.wikipedia.org/wiki/Brucella#Genomics.

A BLAST comparison of this genome with those two previously reported showed major differences summarized in Figure 1 and detailed in Supplementary Table S2. Three major deletions in chromosome I in B. abortus 2308W were detected relative to B. abortus 2308 and 2308A of 1.13, 3.4, and 5.7 kb. It is important to point out that these three deletions are located at the end of contigs and surrounded by repetitive sequences so mayrepresent incompletely assembled contiguous sequences, rather than genuine deletions. Regions containing repetitive sequences are proven difficult to assemble and sequence, regardless the sequencing technology used (Bidmos and Bayliss, 2014). The 1.13 kb deletion contains the 2308 loci BAB1_0934 to BAB1_0937 associated to IS711 elements and transposases. The 3.4 kb deletion includes loci BAB1_1102–BAB1_1104; the first two CDSs are predicted proteins of unknown function and the third one is predicted to be a site-specific recombinase, DNA invertase. These loci fall within a larger 8.1 kb genomic island named GI-1, encoding mainly predicted proteins and phage-related proteins (Rajashekara et al., 2004). Five CDSs are part of the 5.7 kb region absent in chromosome I of 2308W (BAB1_2221–BAB1_2225). They encode tRNAs and rRNAs genes with copies elsewhere in the genome (Supplementary Table S2). These copies are represented by a higher number of reads as compared with average for the rest of the CDS (9700 versus 3200 reads), suggesting that contig breaks are caused by mis-assembly of the repeated copies of the tRNA and rRNA genes during the assembly.

FIGURE 1

FIGURE 1. Graphic representation of BLAST comparison of Brucella abortus 2308 (light gray bar), 2308W (light green bar) and 2308A (dark gray bar), visualized in ACT. Base pairs in kb coordinates are shown in the upper line for each chromosome. The red squares over the bars highlight an approximate region where the main differences among the genomes are found, and a zoom in is represented below. The red bars in the middle of the genomes’ bars indicate the regions are identical. White segments indicate deletions and the blue “hourglass-like” figure means an inversion. Chromosome I of 2308W shows the absence of three regions, which are present in the other two genomes. These regions coincide with contig breaks. At chromosome II, a 5.2 kb insertion is observed in 2308W and 2308A, and an inversion is present in 2308W. All deletions and inversions are surrounded by insertion elements.

An insertion of a 5.2 kb segment in chromosome II was observed in 2308W and 2308A, relative to 2308 (Chain et al., 2005). This region displayed even read coverage as compared to the rest of the 2308W genome. It includes four genes encoding sugar binding proteins, proteins involved in nitrogen metabolism as well as the transcriptional activator FtrB involved in regulation of carbon and amino acid metabolism (BAW_20862–BAW_20865). The adjacent loci (BAW_20861 and BAW_20866) contain sequence partial deletions in 2308 (corresponding to BAB2_0903) or 2308W (corresponding to BAB2_0904). Also evident is an inversion of the 47 kb GI-5 (Rajashekara et al., 2004) in 2308W as compared to 2308 and 2308A that results in deletion of BAB2_1075, encoding an IS3 transposase. Notably, insertion elements were detected alongside all of the indel regions. These observations show not only how unstable the B. abortus 2308 genome is but also that chromosomal rearrangements are a source of genetic variability, probably as consequence of IS mediated transposition events (Tsoktouridis et al., 2003; Ocampo-Sosa and García-Lobo, 2008; Mancilla et al., 2012).

According to the National Collection of Cultures (NCTC) in England, reference strains are “stipulated in internationally recognized standard methods as definitive control strains for various microbiological testing procedures⁴.” The fact that Brucella is regarded as a genetically homogeneous genus, but significant differences are present in B. abortus 2308 stored in different laboratories raises a broader question and challenges the concept of “reference strain.” In the case of B. abortus 2308 this is even more relevant, since this strain has been regarded as the canonical challenge organism in vaccine trials where different results about its survival in mice models are reported (Montaraz and Winter, 1986; Ko et al., 2002; Yang et al., 2006). There are probably several reasons for these contrasting results. Some, such as differences due to strain handling can be accounted by using standardized quality control and experimental design protocols. Probably more difficult to control are genome changes. Reference strains are usually isolates obtained from clinical cases used in research labs as infection models, and as such, are subject of genetic modification according to the environment (Hopkins et al., 1981; Turse et al., 2011). A detailed description of the strain used, including a publically available WGS, as well as growth, propagation and maintenance conditions included with reported results could help readers to better assess the reported observations.

Author Contributions

MS-E performed bioinformatics analysis, manual curation, and prepared tables and figures. NR-V performed molecular typing and manual curation. AC-Z, CJ-R, RR, and DC performed manual curation. CC provided data for manual curation. EB-C and CC-D performed biochemical typing and provided data for manual curation. KB performed bioinformatics analysis and edited the manuscript. EC-O and EM analyzed data. NT performed bioinformatics analysis and manual curation. XD and JL performed manual curation, organized, and analyzed data. CG-V designed the study, analyzed data, performed manual curation, and wrote the manuscript. All authors edited the final version of the manuscript.

Funding

This study was funded by Fondos del Sistema del Consejo Nacional de Rectores (FEES-CONARE) Costa Rica, and Wellcome Trust. NR-V, AC-Z, and CJ-R were partially sponsored by University of Costa Rica scholarships. KB is in receipt of a WT clinical research fellowship (106690/A/14/Z). KB and NT were funded by the Wellcome Trust, grant number 098051

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

We are grateful to Daphnne Garita for technical assistance. We thank Dr. Andrew Page and Dr. Gordon Dougan for helpful discussions.

Supplementary Material

The Supplementary Material for this article can be found online at: http://journal.frontiersin.org/article/10.3389/fmicb.2016.01557

TABLE S1 | 2308W annotation.

TABLE S2 | Genomic differences among 2308W, 2308 and 2308A.

Footnotes

1. ^https://samtools.github.io/bcftools/bcftools.html
2. ^http://www.sanger.ac.uk/resources/software/smalt/
3. ^http://www.ebi.ac.uk/ena/
4. ^https://www.phe-culturecollections.org.uk/collections/nctc.aspx

References

Alton, G. G., Jones, L. M., Angus, R. D., and Verger, J. M. (1988). Techniques for the Brucellosis Laboratory. Paris: Institut National de la Recherche Agronomique.

Google Scholar

Altschul, S. F., Madden, T. L., Schäffer, A. A., Zhang, J., Zhang, Z., Miller, W., et al. (1997). Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 25, 3389–3402. doi: 10.1093/nar/25.17.3389

CrossRef Full Text | Google Scholar

Assefa, S., Keane, T. M., Otto, T. D., Newbold, C., and Berriman, M. (2009). ABACAS: algorithm-based automatic contiguation of assembled sequences. Bioinformatics 25, 1968–1969. doi: 10.1093/bioinformatics/btp347

CrossRef Full Text | Google Scholar

Bidmos, F. A., and Bayliss, C. D. (2014). Genomic and global approaches to unravelling how hypermutable sequences influence bacterial pathogenesis. Pathogens 3, 164–184. doi: 10.3390/pathogens3010164

CrossRef Full Text | Google Scholar

Bosseray, N. (1991). Brucella melitensis Rev. 1 living attenuated vaccine: stability of markers, residual virulence and immunogenicity in mice. Biologicals 19, 355–363. doi: 10.1016/S1045-1056(05)80025-9

CrossRef Full Text | Google Scholar

Carver, T., Bohme, U., Otto, T. D., Parkhill, J., and Berriman, M. (2010). BamView: viewing mapped read alignment data in the context of the reference sequence. Bioinformatics 26, 676–677. doi: 10.1093/bioinformatics/btq010

CrossRef Full Text | Google Scholar

Chain, P., Comerci, D. J., Tolmasky, M. E., Larimer, F. W., Malfatti, S. A., Vergez, L. M., et al. (2005). Whole-genome analyses of speciation events in pathogenic Brucellae. Infect. Immun. 73, 8353–8361. doi: 10.1128/IAI.73.12.8353-8361.2005

CrossRef Full Text | Google Scholar

Dabral, N., Martha-Moreno-Lafont, Sriranganathan, N., and Vemulapalli, R. (2014). Oral immunization of mice with gamma-irradiated Brucella neotomae induces protection against intraperitoneal and intranasal challenge with virulent B. abortus 2308. PLoS ONE 9:e107180. doi: 10.1371/journal.pone.0107180

CrossRef Full Text | Google Scholar

Grilló, M.-J., Blasco, J. M., Gorvel, J. P., Moriyón, I., and Moreno, E. (2012). What have we learned from brucellosis in the mouse model? Vet. Res. 43:29. doi: 10.1186/1297-9716-43-29

CrossRef Full Text | Google Scholar

Halling, S. M., Peterson-Burch, B. D., Bricker, B. J., Zuerner, R. L., Qing, Z., Li, L.-L., et al. (2005). Completion of the genome sequence of Brucella abortus and comparison to the highly similar genomes of Brucella melitensis and Brucella suis. J. Bacteriol. 187, 2715–2726. doi: 10.1128/JB.187.8.2715-2726.2005

CrossRef Full Text | Google Scholar

Harmon, B. G., Adams, L. G., Templeton, J. W., and Smith, R. (1989). Macrophage function in mammary glands of Brucella abortus-infected cows and cows that resisted infection after inoculation of Brucella abortus. Am. J. Vet. Res. 50, 459–465.

Google Scholar

Hopkins, I. G., Thornton, D. H., and Hebert, C. N. (1981). Comparison of Brucella abortus (strain 544) variants. J. Biol. Stand. 9, 421–429. doi: 10.1016/S0092-1157(81)80033-9

CrossRef Full Text | Google Scholar

Jacob, J., Hort, G. M., Overhoff, P., and Mielke, M. E. A. (2006). In vitro and in vivo characterization of smooth small colony variants of Brucella abortus S19. Microbes Infect. 8, 363–371. doi: 10.1016/j.micinf.2005.07.003

CrossRef Full Text | Google Scholar

Jones, L. M., Montgomery, V., and Wilson, J. (1965). Characteristics of carbon dioxide-independent cultures of Brucella abortus isolated from cattle vaccinated with strain 19. J. Infect. Dis. 115, 312–320. doi: 10.1093/infdis/115.3.312

CrossRef Full Text | Google Scholar

Ko, J., Gendron-Fitzpatrick, A., Ficht, T. A., and Splitter, G. A. (2002). Virulence criteria for Brucella abortus strains as determined by interferon regulatory factor 1-deficient mice. Infect. Immun. 70, 7004–7012. doi: 10.1128/IAI.70.12.7004-7012.2002

CrossRef Full Text | Google Scholar

Le Flèche, P., Jacques, I., Grayon, M., Al Dahouk, S., Bouchon, P., Denoeud, F., et al. (2006). Evaluation and selection of tandem repeat loci for a Brucella MLVA typing assay. BMC Microbiol. 6:9. doi: 10.1186/1471-2180-6-9

CrossRef Full Text | Google Scholar

López-Goñi, I., García-Yoldi, D., Marín, C. M., de Miguel, M. J., Barquero-Calvo, E., Guzmán-Verri, C., et al. (2011). New Bruce-ladder multiplex PCR assay for the biovar typing of Brucella suis and the discrimination of Brucella suis and Brucella canis. Vet. Microbiol. 154, 152–155. doi: 10.1016/j.vetmic.2011.06.035

CrossRef Full Text | Google Scholar

Mancilla, M., Marín, C. M., Blasco, J. M., Zárraga, A. M., López-Goñi, I., and Moriyón, I. (2012). Spontaneous excision of the O-polysaccharide wbkA glycosyltranferase gene is a cause of dissociation of smooth to rough Brucella colonies. J. Bacteriol. 194, 1860–1867. doi: 10.1128/JB.06561-11

CrossRef Full Text | Google Scholar

Meyer, M. E., and Morgan, W. J. B. (1973). Designation of neotype strains and of biotype reference strains for species of the genus Brucella meyer and shaw. Int. J. Syst. Bacteriol. 23, 135–141. doi: 10.1099/00207713-23-2-135

CrossRef Full Text | Google Scholar

Miranda, K. L., Dorneles, E. M. S., Pauletti, R. B., Poester, F. P., and Lage, A. P. (2015). Brucella abortus S19 and RB51 vaccine immunogenicity test: evaluation of three mice (BALB/c, Swiss and CD-1) and two challenge strains (544 and 2308). Vaccine 33, 507–511. doi: 10.1016/j.vaccine.2014.11.056

CrossRef Full Text | Google Scholar

Miranda, K. L., Poester, F. P., Minharro, S., Dorneles, E. M. S., Stynen, A. P. R., and Lage, A. P. (2013). Evaluation of Brucella abortus S19 vaccines commercialized in Brazil: immunogenicity, residual virulence and MLVA15 genotyping. Vaccine 31, 3014–3018. doi: 10.1016/j.vaccine.2013.04.054

CrossRef Full Text | Google Scholar

Montaraz, J. A., and Winter, A. J. (1986). Comparison of living and nonliving vaccines for Brucella abortus in BALB/c mice. Infect. Immun. 53, 245–251.

Google Scholar

Moreno, E., and Moriyón, I. (2006). “The genus Brucella,” in The Prokaryotes, eds M. Dworkin, S. Falkow, E. Rosenberg, K.-H. Schleifer, and E. Stackebrandt (New York, NY: Springer), 315–456. doi: 10.1007/0-387-30745-1

CrossRef Full Text | Google Scholar

Mukherjee, F., Jain, J., Grilló, M. J., Blasco, J. M., and Nair, M. (2005). Evaluation of Brucella abortus S19 vaccine strains by bacteriological tests, molecular analysis of ery loci and virulence in BALB/c mice. Biologicals 33, 153–160. doi: 10.1016/j.biologicals.2005.04.003

CrossRef Full Text | Google Scholar

Neta, A. V. C., Mol, J. P. S., Xavier, M. N., Paixão, T. A., Lage, A. P., and Santos, R. L. (2010). Pathogenesis of bovine brucellosis. Vet. J. 184, 146–155. doi: 10.1016/j.tvjl.2009.04.010

CrossRef Full Text | Google Scholar

Ocampo-Sosa, A. A., and García-Lobo, J. M. (2008). Demonstration of IS711 transposition in Brucella ovis and Brucella pinnipedialis. BMC Microbiol. 8:17. doi: 10.1186/1471-2180-8-17

CrossRef Full Text | Google Scholar

Olsen, S. C., McGill, J. L., Sacco, R. E., and Hennager, S. G. (2015). Immune responses of bison and efficacy after booster vaccination with Brucella abortus strain RB51. Clin. Vaccine Immunol. 22, 440–447. doi: 10.1128/CVI.00746-14

CrossRef Full Text | Google Scholar

Page, A. J., De Silva, N., Hunt, M., Quail, M. A., Parkhill, J., Harris, S. R., et al. (2016). Robust high throughput prokaryote de novo assembly and improvement pipeline for Illumina data. Microbiol. Soc. 2, 1–7. doi: 10.1101/052688

CrossRef Full Text | Google Scholar

Palacios-Chaves, L., Conde-Álvarez, R., Gil-Ramírez, Y., Zúñiga-Ripa, A., Barquero-Calvo, E., Chacón-Díaz, C., et al. (2011). Brucella abortus ornithine lipids are dispensable outer membrane components devoid of a marked pathogen-associated molecular pattern. PLoS ONE 6:e16030. doi: 10.1371/journal.pone.0016030

CrossRef Full Text | Google Scholar

Quail, M. A., Kozarewa, I., Smith, F., Scally, A., Stephens, P. J., Durbin, R., et al. (2009). A large genome centre’s improvements to the Illumina sequencing system. Nat. Methods 5, 1005–1010. doi: 10.1038/nmeth.1270.A

CrossRef Full Text

Quail, M. A., Otto, T. D., Gu, Y., Harris, S. R., Skelly, T. F., McQuillan, J. A., et al. (2012). Optimal enzymes for amplifying sequencing libraries. Nat. Methods 9, 10–11. doi: 10.1038/nmeth.1814

CrossRef Full Text | Google Scholar

Rajashekara, G., Glasner, J. D., Glover, D. A., and Splitter, G. A. (2004). Comparative whole-genome hybridization reveals genomic islands in Brucella species. J. Bacteriol. 186, 5040–5051. doi: 10.1128/JB.186.15.5040-5051.2004

CrossRef Full Text | Google Scholar

Rutherford, K., Parkhill, J., Crook, J., Horsnell, T., Rice, P., Rajandream, M.-A., et al. (2000). Artemis: sequence visualization and annotation. Bioinformatics 16, 944–945. doi: 10.1093/bioinformatics/16.10.944

CrossRef Full Text | Google Scholar

Singh, D. K., Kumar, A., Tiwari, A. K., Sankarasubramanian, J., Vishnu, U. S., Sridhar, J., et al. (2015). Draft genome sequence of Brucella abortus virulent strain 544. Genome Announc. 3:e00419-15. doi: 10.1128/genomeA.00419-15

CrossRef Full Text | Google Scholar

Spink, W. (1956). “The evolution of the concept that brucellosis is a disease of animals and man,” in The Nature of Brucellosis (Minneapolis, MN: University of Minnesota Press), 464.

Trant, C. G. M. C., Lacerda, T. L. S., Carvalho, N. B., Azevedo, V., Rosinha, G. M. S., Salcedo, S. P., et al. (2010). The Brucella abortus phosphoglycerate kinase mutant is highly attenuated and induces protection superior to that of vaccine strain 19 in immunocompromised and immunocompetent mice. Infect. Immun. 78, 2283–2291. doi: 10.1128/IAI.01433-09

CrossRef Full Text | Google Scholar

Truong, Q. L., Cho, Y., Park, S., Park, B.-K., and Hahn, T.-W. (2016). Brucella abortus mutants lacking ATP-binding cassette transporter proteins are highly attenuated in virulence and confer protective immunity against virulent B. abortus challenge in BALB/c mice. Microb. Pathog. 95, 175–185. doi: 10.1016/j.micpath.2016.04.009

CrossRef Full Text | Google Scholar

Tsoktouridis, G., Merz, C. A., Manning, S. P., Giovagnoli-Kurtz, R., Williams, L. E., Mujer, C. V., et al. (2003). Molecular characterization of Brucella abortus chromosome II recombination. J. Bacteriol. 185, 6130–6136. doi: 10.1128/JB.185.20.6130-6136.2003

CrossRef Full Text | Google Scholar

Turse, J. E., Pei, J., and Ficht, T. A. (2011). Lipopolysaccharide-deficient Brucella variants arise spontaneously during infection. Front. Microbiol. 2:54. doi: 10.3389/fmicb.2011.00054

CrossRef Full Text | Google Scholar

WHO (2014). “The control of neglected zoonotic diseases,” in WHO Conference Report, ed. NZD4 organising committee, (Geneva: WHO Press), 48.

Google Scholar

Xavier, M. N., Paixão, T. A., Poester, F. P., Lage, A. P., and Santos, R. L. (2009). Pathological, immunohistochemical and bacteriological study of tissues and milk of cows and fetuses experimentally infected with Brucella abortus. J. Comp. Pathol. 140, 149–157. doi: 10.1016/j.jcpa.2008.10.004

CrossRef Full Text | Google Scholar

Yang, X., Becker, T., Walters, N., and Pascual, D. W. (2006). Deletion of znuA virulence factor attenuates Brucella abortus and confers protection against wild-type challenge. Infect. Immun. 74, 3874–3879. doi: 10.1128/IAI.01957-05

CrossRef Full Text | Google Scholar

Zerbino, D. R., and Birney, E. (2008). Velvet: algorithms for de novo short read assembly using de Bruijn graphs. Genome Res. 18, 821–829. doi: 10.1101/gr.074492.107

CrossRef Full Text | Google Scholar

Zhu, L., Feng, Y., Zhang, G., Jiang, H., Zhang, Z., Wang, N., et al. (2016). Brucella suis strain 2 vaccine is safe and protective against heterologous Brucella spp. infections. Vaccine 34, 395–400. doi: 10.1016/j.vaccine.2015.09.116

CrossRef Full Text | Google Scholar