Although eradicated throughout much of the developed world, bovine brucellosis continues to be common in Southern Italy [1]. Brucella abortus, the etiological agent, is a Gram-negative non-motile and non-sporulating bacterium that forms coccobacilli with oxidative metabolism [2]. At the beginning of the 20th century, B. abortus was recognized by Benhard Bang as the causative agent of epizootic abortion in cattle. Only later was this disease connected to Malta Fever, a disease caused by a closely related species, B. melitensis, which humans typically contract from infected sheep and goats [3]. Meyer classified these and other species in the genus Brucella [4]. In cattle, the disease presents primarily as reproductive disorders such as abortion, resulting in substantial economic losses to agriculture [5]. For example, the economic burden of bovine brucellosis in Latin America has been estimated at $600 million annually [6]. Brucella abortus is a significant zoonotic agent, and humans are typically infected by consumption of raw dairy products, exposure to aerosolized bacteria at slaughter, or during veterinary care and animal husbandry [7].

The Italian government has conducted a country-wide eradication program for cattle brucellosis since 1994 [8]. Animal cases are currently limited to seven regions of southern Italy, with the highest prevalence of infection in areas of Sicily, Calabria, and Apulia [9]. This eradication program, along with strict regulations on cattle movements, has reduced the prevalence and geographic distribution of brucellosis in Italy [10]. In order to achieve the eradication of brucellosis, the identification of the sources of infection is essential, as is determining whether the same strains are circulating within a region or continuously reintroduced from other areas or countries. Additionally, determining the host species and source regions are important components of epidemiological tracking of transmission to humans.

Molecular epidemiology has been challenging in Brucella due to relatively high levels of nucleotide similarity throughout the genus [11]. Traditional molecular approaches involving biochemical testing are only able to categorize samples into broad phenotypic groups of biovars; eight different biovars have been described for B. abortus (biovars 1–7 and 9) [12]. More recent genetic analyses using variable number tandem repeat (VNTR) markers have given much higher resolution and can be highly informative in some cases [13–16], but homoplasy in VNTR markers and a lack of sampling for some regions make it challenging to place B. abortus into a global phylogeny [17].

We used whole genome sequencing (WGS) to place 11 Italian B. abortus isolates onto a global context. We also identified specific single nucleotide polymorphisms (SNPs) that define the three major Italian clades in our study, and then developed clade-specific assays to genotype our Italian collection of 261 isolates. This genotyping allowed us to rapidly and inexpensively identify the lineages of B. abortus circulating in Southern Italy and place Italian strains in a global phylogeny. In Italy, where brucellosis remains a major public health concern, the information generated by this genotyping will also support epidemiological investigations of human infections.

Whole genome comparisons generated approximately 9,000 putative SNPs that were used to construct a phylogenetic tree using maximum parsimony. Italian isolates of B. abortus were part of three clades that are basal to the most common worldwide lineage, which we refer to as the biovar 1, 2, 4 lineage (hereafter biovar 1/2/4); a lineage that is ubiquitous in North America and parts of Europe (Fig. 1). Isolates from western Italy and a closely related sample from France form a “West Italia” subclade defined by 83 SNPs, the top most clade in Fig. 1 containing Italian isolates in red. A broader clade that also includes isolates from Poland, Spain, and China is defined by 91 SNPs. Members of the West Italia subclade and broader clade are most closely related to the biovar 1/2/4 lineage, which is consistent with the biochemical designation of most of these isolates as biovar 1.

The other Italian isolates are part of an even more basal and considerably more diverse group of B. abortus genomes. Unlike the biovar 1/2/4 lineage, this basal group is comprised largely of isolates from biovars 3, 5, 6, 7, 9—isolates that are broadly distributed across Europe, Asia, and Africa. Substantial diversity among African isolates suggests this continent may be the region of origin of B. abortus as a species (JTF, unpubl. data) but additional sampling is necessary to test this hypothesis. Inconsistencies in biovar assignment compared to their placement in the phylogenetic tree suggest that biovar classification in B. abortus does not consistently reflect genetic relationships or that biochemical determination of biovars is unreliable.

Six of our isolates from across Italy form a “Trans Italia” subclade along with isolates from Mongolia, southern Europe, and Egypt, which are defined by 44 SNPs. The broader clade also includes an isolate from Turkey. The 169 SNPs exclusive to this clade suggests a fairly long time period separates these isolates from the common ancestor they share with genomes from Uganda and Iraq. Isolates from eastern Italy form an “East Italia” subclade, with 157 SNPs exclusive to the branch. The broader clade contains an isolate from Greece, which shares 145 exclusive SNPs with the East Italia cluster. Members of the East Italia subclade are closely related, although these three isolates came from three different regions.

We developed 13 Melt-MAMA assays to genotype a larger collection of 261 samples and to better define the population substructure of B. abortus in Southern Italy. Six of the Melt-MAMA assays produced consistent results and corresponded to two assays each for the West Italia, Trans Italia, and East Italia subclades and the broader clades to which they belong (Table 1). Our genotyping indicated that 13 Italian isolates came from the West Italia subclade, 214 isolates were part of the Trans Italia subclade, and 34 isolates were part of the East Italia subclade (Additional file 1: Table S1). We then mapped the geographic distributions of the clades to determine their phylogeographic patterns. The West Italia genotypes consisted entirely of isolates from biovar 1 detected in the three western Italian districts of Campania, Basilicata, and Calabria (Fig. 2). The East Italia group was found mostly in the eastern regions of southern Italy: Apulia, Molise and Abruzzo (Fig. 2). Finally, the Trans Italia group was distributed throughout endemic regions of Southern Italy, suggesting its ubiquity and dominance in the region (Fig. 2).

Branch	Location in B. abortus 2308	SNP allele	Melt-MAMA primer	Melt-MAMA Primer sequences (5′ to 3′)	PCR mix Concentration	Annealing Temp.	Melting Temp.
Branch 1 West Italia	951620 NC_007618.1	C/T	A	GCGACCGCATCCTCTAaTC	0.4 μM		80.8 ° C
			D	cggggcggggcggggGCGACCGCATCCTCTAcTT	0.1 μM	56 ° C	83.4 ° C
			C	CATGGTAGGCAGGCTGGTT	0.5 μM
Branch 2	1336380	A/G	A	CGAAGTTTCCCGACACCcA	0.4 μM		81.2 ° C
West Italia	NC_007618.1		D	cggggcggggcggggCGAAGTTTCCCGACACCtG	0.1 μM	56 ° C	83.9 ° C
			C	GATGAAGAAATCGATGAGCG	0.5 μM
Branch 1	104052	G/T	A	TGGAGGACGAAGTCTTCGaTG	0.25 μM		78 ° C
East Italia	NC_007618.1		D	cggggcggggcggggGGAGGACGAAGTCTTCGcTT	0.25 μM	56 ° C	82.9 ° C
			C	ATCCTCGACAGTTGCGGC	0.5 μM
Branch 2	360724	C/A	A	TACCGCGGCCCGTCTAtC	0.25 μM		81.6 ° C
East Italia	NC_007618.1		D	cggggcggggcggggTACCGCGGCCCGTCTAaA	0.25 μM	56 ° C	83.7 ° C
			C	AATCAAAGCGGCTTGCCTC	0.5 μM
Branch 1	740090	T/C	A	CAAAAGCGGACGGTCAgAT	0.4 μM		78.9 ° C
Trans Italia	NC_007618.1		D	cggggcggggcggggTCAAAAGCGGACGGTCAtAC	0.1 μM	56 ° C	80.9 ° C
			C	CCCTCAACCCAGACGTTCG	0.5 μM
Branch 2	864308	C/T	A	AGGCATTACCGTTTCAGATTtC	0.4 μM		78.5 ° C
Trans Italia	NC_007618.1		D	cggggcggggcggggAGGCATTACCGTTTCAGATTaT	0.1 μM	56 ° C	81.9 ° C
			C	TCTTTCCCGGCTCCGAAT	0.5 μM

Melt-MAMA primer abbreviations: A ancestral, D derived, C common

Details of primers concentration used in each reaction with proper annealing temperature is given. Melting temperatures detected for both the amplicons derived by the specific primer binding is given according to the experiments run on an Applied Biosystems 7500 Real-Time PCR System

Whole genome sequencing and comparison to a large global collection of genomes allowed us to determine the evolutionary history of B. abortus in Italy. Despite the strains being isolated exclusively in 2011, the use of multi-locus variable number tandem repeat (VNTR) analysis (MLVA) in a previous study allowed us to select and then sequence a diverse collection of isolates from throughout endemic districts in Southern Italy [17]. Challenges with VNTRs limited our phylogenetic inference for B. abortus for the country. However, genome-wide phylogenetics using SNPs identified that Italian strains come from three major clades, all of which are substantially different than the common biovar 1/2/4 lineage found in much of Europe and North America. Genotyping assays supported these findings in a larger collection of 200 isolates from Italy.

We believe that the relatedness of Italian strains to lineages from other parts of Europe, Asia, and Africa suggests that B. abortus in Italy has a very different evolutionary history than the dominant lineage of biovar 1/2/4 that was brought to the New World by various Europeans in infected cattle. Two or three different domestication events of wild cattle in Eurasia have occurred [18], and then cattle, with their associated diseases, were introduced by humans throughout the world [19]. Italian farms contain several European cattle breeds, but also contain several primitive cattle breeds such as Podolica, Marchigiana, Chianina, and Romagnola, which are all directly linked to the ancient cattle Bos primigenius from Eurasia [20]. Current brucellosis programs strictly control cattle movements, but these programs have been only recently established, and B. abortus has likely been influenced by several centuries of introductions and disseminations of infected cattle throughout Italy. Thus, contemporary B. abortus diversity and distributions may reflect past mixing of cattle of different origins.

The newly developed Melt-MAMA assays specific to the three Italian clades and subclades allowed us to screen a large collection of isolates and quickly and inexpensively assign them to one of three genetic groups. Future testing of additional isolates using these six SNP assays during eradication programs will allow them to be typed more rapidly and confidently than biotyping, and could serve Italian diagnostic laboratories as a first line assay. Genotyping coupled with knowledge of herd movement through use of risk maps, spatial modeling, and Social Network Analysis (SNA) techniques should enhance the reliability of the Italian brucellosis eradication plan. Although a multiplex PCR approach has been developed recently [14], and the MLVA-16 has been shown to be robust and informative [21], these techniques remain cumbersome for data production and standardization. Perhaps a more reasonable approach would be to pair the SNP assays with a limited number of VNTR loci to achieve a higher resolution for molecular epidemiology purposes.

We recognize the capabilities of whole genome sequencing, but we are also aware that this technology is not yet accessible to all labs and requires skills that are not present in many diagnostic laboratories. Nonetheless, for some outbreak investigations, MLVA-16 or whole genome sequencing may be essential when complete epidemiological traceback is required. Ultimately, our study demonstrates the utility of WGS SNPs paired with extensive epidemiologic data for analyzing the distribution of B. abortus isolates throughout endemic regions.