Whole-genome-based phylogeny and divergence of the genus Brucella

Genomic insights into Brucella melitensis in India: stability of ST8 and the role of virulence genes in regional adaptations

Brucella melitensis is a highly infectious zoonotic pathogen responsible for brucellosis, which significantly affects both human and livestock health worldwide. This study employed whole-genome sequencing (WGS) to analyze the genetic diversity of 24 B. melitensis isolates from India. Pangenome analysis revealed a highly conserved nature with the involved strains having very limited accessory genes. Multilocus Sequence Typing (MLST) identified sequence type ST8 as predominant among Indian strains. Analysis of virulence genes revealed a total of 43 virulence-related genes in all strains, emphasizing their critical role in the pathogenicity of B. melitensis. Unique gene profiles and distinct phylogenetic clusters suggest regional adaptations and evolutionary pressures. The comprehensive genomic insights from this study help to elucidate the geographic distribution and interspecies transmission of Indian strains, highlighting the importance of targeted brucellosis control measures in India. Additionally, the identification of conserved virulence genes involved in immune evasion and intracellular survival highlights their importance in the bacterium's pathogenicity. This research contributes to the global understanding of B. melitensis genomic diversity, providing valuable insights for broader epidemiological studies and brucellosis management strategies worldwide.IMPORTANCEB. melitensis is a significant cause of illness in both humans and animals, particularly in India, where the disease remains a major concern. This study highlights that only a few genetic types of the bacteria are circulating in the region, which means control efforts can be better focused on these specific types. By understanding the unique characteristics of Indian strains, and how these strains spread and adapt, this research offers valuable guidance for improving brucellosis prevention strategies. These insights can help in developing more effective diagnostic tools, enhancing vaccination efforts, and strengthening disease control programs to reduce the impact of brucellosis on public health and livestock industries.

Molecular epidemiology of brucellosis in Asia: insights from genotyping analyses

Brucellosis infects humans and animals worldwide but is particularly prevalent in Asia. In many Asian countries, molecular diagnostic tools for accurate molecular diagnostics and molecular epidemiology are lacking. Nonetheless, some countries have conducted in-depth molecular epidemiological studies. The objective of this study was to reveal the genetic relationships, geographic origins, and distributions of Brucella strains across Asia for two primary species, B. abortus and B. melitensis. For this, we systematically searched genotyping data from published studies on the molecular epidemiology of Brucella species for both humans and livestock in Asia. We used data from multilocus sequence typing (MLST), multiple-locus variable-number tandem repeat analysis (MLVA), and whole genome sequencing analysis of Brucella strains. We also analyzed the MLVA genotypes of 129 B. abortus isolates and 242 B. melitensis isolates with known origins in Asia from an online MLVA database using MLVA-11 data in minimum spanning trees and MLVA-16 data in neighbor-joining trees. We found that the B. melitensis East Mediterranean lineage is predominant across the continent, with only a small number of samples from the Africa and Americas lineages, and none from the West Mediterranean lineage. The "abortus C" genotype was the most common group of B. abortus in Asia, with limited genetic variation for this species. Several studies also reported that Near Eastern countries frequently encounter human brucellosis cases of B. abortus from genotypes 42 and 43. Our study highlights the inconsistent collection of genetic data for Brucella species across Asia and a need for more extensive sampling in most countries. Finally, a consistent nomenclature is necessary to define various groupings of strains within a lineage (i.e., clade) so uniform terminology should denote particular genetic groups that are understood by all researchers.

Closely Related Brucella Species Widely Differ in their Vegetative and Intracellular Growth

Growth rate is a key prokaryotic trait that allows for estimating fitness and understanding cell metabolism. While it has been well studied in model organisms, there is limited data on slow-growing bacteria. In particular, there is a lack of quantitative studies on Brucella species. This genus includes important microorganisms that are causative agents of brucellosis, one of the most widespread bacterial zoonoses, affecting several species of animals, including humans. Brucella species exhibit approximately 97% genomic similarity. Despite this, Brucella species show different host preferences, zoonotic risks, and pathogenicity. After more than one hundred years of research, numerous aspects of Brucella biology, such as in vivo and in vitro growth, remain poorly characterized. In this work, we analyzed vegetative and intracellular growth of the classical Brucella species (B. abortus, B. melitensis, B. suis, B. ovis, and B. canis). Strikingly, each species displayed distinct growth parameters in culture. Doubling time (DT) ranged from 2.7 hs-1 in B. suis to 18 h-1 for B. ovis. In the context of intracellular infection of J774A.1 phagocytic cells, DT was longer, but it widely varied across species, closely correlating with the growth observed in vitro. Overall, and despite high genome similarity, we also found species-specific growth parameters in the intracellular cell cycle.

Whole-genome sequencing-based analysis of Brucella species isolated from ruminants in various regions of Türki̇ye

BACKGROUND

Brucellosis, a zoonotic disease in Türkiye, which has significant direct and indirect impacts on the healthcare system and livestock. This study, which aimed to investigate the differences among Brucella spp. isolates originating from different regions of Türkiye, for implications for public health and veterinary medicine.

METHOD

Twenty-one isolates from ruminants and two isolates from humans obtained from various regions of Türkiye were utilized in the study. The isolates were identified and biotyped using traditional microbiological procedures, and whole-genome sequencing (WGS) was performed. This was followed by single nucleotide polymorphism (SNP)--based phylogenetic analysis and WGS-based analysis of virulence and resistance genes. Additionally, phenotypic antimicrobial resistance and phage susceptibilities were determined. The obtained data were then compared for concordance, ensuring the validity and reliability of the results.

RESULTS

Our study, employing culture methods, polymerase chain reaction (PCR), and WGS analyses, identified 11 Brucella melitensis (bv 3 (n = 9), one each bv 1 and bv 2) and 12 B. abortus (bv 3 (n = 11), bv 9 (n = 1)) isolates All B. abortus isolates were of bovine origin, while the B. melitensis isolates were from sheep (n = 7), goat (n = 1), ram (n = 1), and humans (n = 2). In the whole-genome SNP-based phylogenetic tree, all B. melitensis strains were found to be of the IIb subtype of genotype II associated with the Eastern Mediterranean lineage. Ten different genotypes were identified in the SNP analysis of the isolates, with a maximum SNP difference of 278 and a minimum SNP difference of 4 among these genotypes. According to the WGS-SNP-based phylogenetic tree of B. abortus isolates, they were grouped in clade C1. In the SNP analysis, where ten different genotypes were identified, the SNP difference among these genotypes was a maximum of 316 and a minimum of 6. In the in silico MLST analysis performed with WGS data, B. melitensis isolates were identified as ST8 and ST102 genotypes, while B. abortus isolates were identified as ST2 and ST3 genotypes. The dominant genotypes were ST8 for B. melitensis and ST2 for B. abortus, respectively. Virulence gene analysis conducted based on WGS data of the 23 B. abortus and B. melitensis isolates revealed 43 virulence gene-associated regions in all strains, irrespective of species, host, or isolation year. Although classical resistance-related genes were not detected by WGS-based antimicrobial resistance gene analysis, phenotypic resistance analysis revealed resistance to azithromycin, rifampin, and trimethoprim/sulfamethoxazole in B. abortus and B. melitensis isolates.

CONCLUSION

Both B. melitensis and B. abortus were circulating species in animals and human. The dominant genotypes were ST8 for B. melitensis and ST2 for B. abortus, respectively. All B. melitensis strains were found to be of the IIb subtype of genotype II associated with the Eastern Mediterranean lineage, while B. abortus isolates, they were grouped in clade C1. Further, a comprehensive study with a sufficient number of isolates covering all regions of Türkiye would provide more accurate information about the current epidemiological situation in the country.

Whole-genome analysis reveals phylogenetic and demographic history of Eurasian perch

The contemporary diversity and distribution of species are shaped by their evolutionary and ecological history. This can be deciphered with the help of phylogenetic and demographic analysis methods, ideally combining and supplementing information from mitochondrial and nuclear genomes. In this study, we investigated the demographic history of Eurasian perch (Perca fluviatilis), a highly adaptable teleost with a distribution range across Eurasia. We combined whole-genome resequencing data with available genomic resources to analyse the phylogeny, phylogeography, and demographic history of P. fluviatilis populations from Europe and Siberia. We identified five highly diverged evolutionary mtDNA lineages, three of which show a strong signal of admixture in the Baltic Sea region. The estimated mean divergence time between these lineages ranged from 0.24 to 1.42 million years. Based on nuclear genomes, two distinct demographic trajectories were observed in European and Siberian samples reflecting contrasting demographic histories ca. 30,000-100,000 years before the present. A comparison of mtDNA and nuclear DNA evolutionary trees and AMOVA revealed concordances, as well as incongruences, between the two types of data, most likely reflecting recent postglacial colonization and hybridization events. Overall, our findings demonstrate the power and usefulness of genome-wide information for delineating historical processes that have shaped the genome of P. fluviatilis. We also highlight the added value of data-mining existing transcriptomic resources to complement novel sequence data, helping to shed light on putative glacial refugia and postglacial recolonization routes.

Comparative genomics of Brucella abortus and Brucella melitensis unravels the gene sharing, virulence factors and SNP diversity among the standard, vaccine and field strains

Brucella abortus and Brucella melitensis are the primary etiological agents of brucellosis in large and small ruminants, respectively. There are limited comparative genomic studies involving Brucella strains that explore the relatedness among both species. In this study, we involved strains (n=44) representing standard, vaccine and Indian field origin for pangenome, single nucleotide polymorphism (SNP) and phylogenetic analysis. Both species shared a common gene pool representing 2884 genes out of a total 3244 genes. SNP-based phylogenetic analysis indicated higher SNP diversity among B. melitensis (3824) strains in comparison to B. abortus (540) strains, and a clear demarcation was identified between standard/vaccine and field strains. The analysis for virulence genes revealed that virB3, virB7, ricA, virB5, ipx5, wbkC, wbkB, and acpXL genes were highly conserved in most of the Brucella strains. Interestingly, virB10 gene was found to have high variability among the B. abortus strains. The cgMLST analysis revealed distinct sequence types for the standard/vaccine and field strains. B. abortus strains from north-eastern India fall within similar sequence type differing from other strains. In conclusion, the analysis revealed a highly shared core genome among two Brucella species. SNP analysis revealed B. melitensis strains exhibit high diversity as compared to B. abortus strains. Strains with absence or high polymorphism of virulence genes can be exploited for the development of novel vaccine candidates effective against both B. abortus and B. melitensis.

Vibrio vulnificus mutation rate: an in vitro approach

Vibrio vulnificus is a multi-host pathogenic species currently subdivided into five phylogenetic lineages (L) plus one pathovar with the ability to infect fish due to a transmissible virulence plasmid. This plasmid (or a fragment of it) has been transmitted between lineages within the species, contributing to the evolution of V. vulnificus. This study aimed to provide an experimental approximation to the V. vulnificus mutation rate by determining spontaneous mutation rates from bacterial cultures of representants of the different lineages by whole-genome sequencing. To this purpose, synonymous SNP differences, i.e., spontaneous mutation not subjected to the evolutive forces, between initial and final culture after serial growth were evaluated and used for mutation rate calculation.

Whole-genome sequencing for genetic diversity analysis of Iranian Brucella spp. isolated from humans and livestock

Brucellosis is one of the most common zoonoses in the Middle East. It is causing economic losses to the livestock industry and has a great public health concern. Little is known about the genetic diversity and distribution of brucellae in Iran. Therefore, forty Brucella spp. strains (B. abortus and B. melitensis) isolated from animals and humans were analyzed by whole genome sequencing (WGS) technology using single nucleotide polymorphism (SNP) analysis and core genome multilocus sequence typing (cgMLST). Brucella isolates were obtained from lymph nodes (cows and camels), milk (cows, camels and sheep), and aborted foetus samples (sheep and goats), as well as cerebrospinal fluid and blood of humans. The isolates were originating from thirteen provinces of Iran and isolated between 2015 and 2020. According to in-silico MLST, ST8 and ST2 were the most frequent sequence types in B. melitensis and B. abortus, respectively. Based on phylogeographic reconstruction using cgSNP analysis, the investigated Iranian B. melitensis strains belonged to the American and Mediterranean lineages of the B. melitensis phylogeny. Furthermore, cgSNP analysis revealed a similarity between Iranian B. abortus isolates and strains from Iraq and Egypt. Therefore, the origin of the Iranian strains can be suggested to be strains from neighboring and Middle East countries. Moreover, cgMLST analysis showed that the Iranian B. melitensis strains were closely relative to strains recovered from sheep and humans in Iraq, Afghanistan, Syria, Turkmenistan, and Pakistan. In the current panel of strains, cgMLST and cgSNP analysis provided an appropriate and accurate tool for effective traceback analyses for Brucella spp. from Iran. The results of cgSNP and cgMLST helped to understand the geographic distribution and interspecies transmission of Iranian strains and highlight the importance of specific brucellosis control measures in Iran with regard to the One-Health approach.

Genomic distribution of the insertion sequence IS711 reveal a potential role in Brucella genome plasticity and host preference

The Insertion Sequence 711 (IS711) is linked to the Brucella genus. Mapping the genomic distribution of IS711 can help understand this insertion element's biological and evolutionary role. This work aimed to delineate the genomic distribution of the IS711 element and to study its association with Brucella evolution. A total of 124 genomes representing 9 Brucella species were searched using BLASTn sequence alignment tool to identify complete and truncated copies of IS711. Based on the genomic context, each IS711 locus was assigned a code using the initial letters of its neighboring genes. Various tools were used to annotate the neighboring genes and determine the shared synteny around orthologous IS711 loci. The tool Islandviewer 4 was used to scan for genomic islands. The Codon Tree method was used to build phylogenetic trees of B. melitensis, B. abortus, and B. suis genomes. The phylogenetic trees of the three species were analyzed, taking into account the genomic distribution patterns of IS711. The result of IS711 frequency analysis showed a relatively conserved number of copies/genome for the different species and for some biovars. The analysis showed that Brucella species with a relatively low IS711 copy number (4-8 copies/genome) are linked to domestic animals as primary hosts and have potential for zoonotic transmission. However, species with a relatively higher copy number (12-30 copies/genome) are less zoonotic and tend to be linked with wild animals as primary hosts. Analyzing the genomic distribution map of IS711 loci showed several unique patterns of IS711 distribution that are correlated with the evolution of Brucella species and biovars. The results also showed that 46.2% of the conserved IS711 elements are located within genomic islands. Based on our results and previous data, we postulate a model explaining the IS711 role in Brucella evolution. We assume that during the transition from a free-living to an intracellular lifestyle, a descendant of the Brucella genus had acquired a progenitor sequence of the IS711. Subsequently, a burst in IS711 transposition occurred. This parasitic expansion can be deleterious and has to be counteracted by evolutionary forces to prevent lineage extension and to promote adaptation to host. Similar to other plasmid-free pathogenic α-Proteobacteria bacteria, the balance of expansion and reduction of insertion elements could be one of the mechanisms to control genome reduction and streamlining. We hypothesize that the IS711-mediated genomic changes and other small sequence nucleotide changes in specific orthologous genes could significantly contribute to Brucella's evolution and adaptation to different animal hosts.

Seroprevalence of Brucella Infection in Wild Boars (Sus scrofa) of Bavaria, Germany, 2019 to 2021 and Associated Genome Analysis of Five B. suis Biovar 2 Isolates

Brucella species are highly pathogenic zoonotic agents and are found in vertebrates all over the world. To date, Germany is officially declared free from brucellosis and continuous surveillance is currently limited to farm ruminants. However, porcine brucellosis, mostly caused by B. suis biovar 2, is still found in wild boars and hares. In the present study, a three-year monitoring program was conducted focusing on the wild boar population in the state of Bavaria. Serologic screening of 11,956 animals and a direct pathogen detection approach, including a subset of 681 tissue samples, was carried out. The serologic incidence was 17.9%, which is in approximate accordance with previously published results from various European countries. Applying comparative whole genome analysis, five isolated B. suis biovar 2 strains from Bavaria could be assigned to three known European genetic lineages. One isolate was closely related to another strain recovered in Germany in 2006. Concluding, porcine brucellosis is endemic in Bavaria and the wild boar population represents a reservoir for genetically distinct B. suis biovar 2 strains. However, the transmission risk of swine brucellosis to humans and farm animals is still regarded as minor due to low zoonotic potential, awareness, and biosafety measures.

Determination of Virulence-Associated Genes and Antimicrobial Resistance Profiles in Brucella Isolates Recovered from Humans and Animals in Iran Using NGS Technology

Brucellosis is a common zoonotic disease in Iran. Antimicrobial-resistant (AMR) Brucella isolates have been reported from different developing countries, posing an imminent health hazard. The objective of this study was to evaluate AMR and virulence-associated factors in Brucella isolates recovered from humans and animals in different regions of Iran using classical phenotyping and next generation sequencing (NGS) technology. Our findings revealed that B. melitensis is the most common species in bovines, small ruminants and camels. B. abortus was isolated only from one human case. Probable intermediate or resistant phenotype patterns for rifampicin, trimethoprim-sulfamethoxazole, ampicillin-sulbactam and colistin were found. Whole genome sequencing (WGS) identified mprF, bepG, bepF, bepC, bepE, and bepD in all isolates but failed to determine other classical AMR genes. Forty-three genes associated with five virulence factors were identified in the genomes of all Brucella isolates, and no difference in the distribution of virulence-associated genes was found. Of them, 27 genes were associated with lipopolysaccharide (LPS), 12 genes were related to a type IV secretion system (virB1-B12), two were associated with the toll-interleukin-1 receptor (TIR) domain-containing proteins (btpA, btpB), one gene encoded the Rab2 interacting conserved protein A (ricA) and one was associated with the production of cyclic β-1,2 glucans (cgs). This is the first investigation reporting the molecular-based AMR and virulence factors in brucellae isolated from different animal hosts and humans in Iran. Iranian B. abortus and B. melitensis isolates are still in vitro susceptible to the majority of antibiotics used for the treatment of human brucellosis. WGS failed to determine classical AMR genes and no difference was found in the distribution of virulence-associated genes in all isolates. Still, the absence of classical AMR genes in genomes of resistant strains is puzzling, and investigation of phenotypic resistance mechanisms at the proteomic and transcriptomic levels is needed.

Evidente-a visual analytics tool for data enrichment in SNP-based phylogenetic trees

Motivation

A common practice in the analysis of pathogens and their strains is using single-nucleotide polymorphisms (SNPs) to reconstruct their evolutionary history. However, genome-wide SNP-based phylogenetic trees are rarely analyzed without any further information. Including the underlying SNP data together with further metadata on the respective samples in the exploration process can facilitate linking the genomic and phenotypic properties of the samples.

Results

We introduce Efficient VIsual analytics tool for Data ENrichment in phylogenetic TreEs (Evidente), a web-application that provides an interactive visual analysis interface for the simultaneous interrogation of phylogenetic relationships, genome-wide SNP data and metadata for samples of an organism. Besides visualizing the phylogenetic tree, Evidente classifies SNPs as supporting or non-supporting of the tree structures and shows the distribution of both types of SNPs among samples and clades of interest. Furthermore, additional metadata can be included in the visualization. Lastly, Evidente includes an enrichment analysis to identify over-represented genomic features encoded by GO-terms within the clades of the tree. We demonstrate the usability of Evidente with the data of the pathogens Treponema pallidum and Mycobacterium leprae.

Availability and implementation

Evidente is available at the TueVis visualization web server at https://evidente-tuevis.cs.uni-tuebingen.de/, it can also be run locally.

Supplementary information

Supplementary data are available at Bioinformatics Advances online.

Assays for Identification and Differentiation of Brucella Species: A Review

Brucellosis is one of the most important and widespread bacterial zoonoses worldwide. Cases are reported annually across the range of known infectious species of the genus Brucella. Globally, Brucella melitensis, primarily hosted by domestic sheep and goats, affects large proportions of livestock herds, and frequently spills over into humans. While some species, such as Brucella abortus, are well controlled in livestock in areas of North America, the Greater Yellowstone Ecosystem supports the species in native wild ungulates with occasional spillover to livestock. Elsewhere in North America, other Brucella species still infect domestic dogs and feral swine, with some associated human cases. Brucella spp. patterns vary across space globally with B. abortus and B. melitensis the most important for livestock control. A myriad of other species within the genus infect a wide range of marine mammals, wildlife, rodents, and even frogs. Infection in humans from these others varies with geography and bacterial species. Control in humans is primarily achieved through livestock vaccination and culling and requires accurate and rapid species confirmation; vaccination is Brucella spp.-specific and typically targets single livestock species for distribution. Traditional bacteriology methods are slow (some media can take up to 21 days for bacterial growth) and often lack the specificity of molecular techniques. Here, we summarize the molecular techniques for confirming and identifying specific Brucella species and provide recommendations for selecting the appropriate methods based on need, sensitivity, and laboratory capabilities/technology. As vaccination/culling approaches are costly and logistically challenging, proper diagnostics and species identification are critical tools for targeting surveillance and control.

Molecular characterization of zoonotic Brucella species isolated from animal and human samples in Iran

Brucellosis is an endemic infection in Iran and represents a serious health problem in humans and livestock causing important economic losses. The objective of this study was to undertake molecular characterization of Brucella spp. isolated from humans and livestock in several provinces of Iran including by multi-locus sequence typing (MLST), in order to understand the genotypes circulating in Iran and their relationship to genotypes globally. A total of 23 Brucella isolates were isolated from eight milk samples (seven cows, and one camel), human blood samples (seven), bovine lymph nodes (two), and samples from aborted fetuses (three sheep, two cows, and one goat). Phenotypic and molecular identification of Brucella isolates was performed on all isolated bacteria and showed that all were either Brucella melitensis or Brucella abortus. B. melitensis was associated with ovine/caprine and camel samples, most human isolates, and a significant minority of cattle isolates. In contrast B. abortus from livestock was associated only with isolations from bovine samples, as well as a single human sample. These results indicate that both B. melitensis and B. abortus contribute to the human brucellosis burden in Iran. B. melitensis isolates comprised three MLST-9 genotypes, the common and globally distributed ST8, a single representative of ST7, and several additional examples of ST102, a genotype previously only reported in a single isolate from a human brucellosis case believed to be acquired through travel to Iran. B. abortus isolates represented two globally common MLST-9 genotypes (ST1 and ST2), with relationships to biotype and other PCR-based typing methods consistent with previous observations. The results provide the basis for further studies examining the molecular epidemiology of Brucella circulating in Iran and the relationships of local isolates to those present globally.

High-Resolution Melting PCR as Rapid Genotyping Tool for Brucella Species

Brucella sp. are the causative agents of brucellosis. One of the main characteristics of the Brucella genus concerns its very high genetic homogeneity. To date, classical bacteriology typing is still considered as the gold standard assay for direct diagnosis of Brucella. Molecular approaches are routinely used for the identification of Brucella at the genus level. However, genotyping is more complex, and to date, no method exists to quickly assign a strain into species and biovar levels, and new approaches are required. Next generation sequencing (NGS) opened a new era into the diagnosis of bacterial diseases. In this study, we designed a high-resolution melting (HRM) method for the rapid screening of DNA and direct assignment into one of the 12 species of the Brucella genus. This method is based on 17 relevant single nucleotide polymorphisms (SNPs), identified and selected from a whole genome SNP (wgSNP) analysis based on 988 genomes (complete and drafts). These markers were tested against the collection of the European Reference Laboratory (EU-RL) for brucellosis (1440 DNAs extracted from Brucella strains). The results confirmed the reliability of the panel of 17 SNP markers, allowing the differentiation of each species of Brucella together with biovars 1, 2, and 3 of B. suis and vaccine strain Rev1 (B. melitensis) within 3 h, which is a considerable gain of time for brucellosis diagnosis. Therefore, this genotyping tool provides a new and quick alternative for Brucella identification based on SNPs with the HRM-PCR assay.

Brucella and Its Hidden Flagellar System

Brucella, a Gram-negative bacterium with a high infective capacity and a wide spectrum of hosts in the animal world, is found in terrestrial and marine mammals, as well as amphibians. This broad spectrum of hosts is closely related to the non-classical virulence factors that allow this pathogen to establish its replicative niche, colonizing epithelial and immune system cells, evading the host's defenses and defensive response. While motility is the primary role of the flagellum in most bacteria, in Brucella, the flagellum is involved in virulence, infectivity, cell growth, and biofilm formation, all of which are very important facts in a bacterium that to date has been described as a non-motile organism. Evidence of the expression of these flagellar proteins that are present in Brucella makes it possible to hypothesize certain evolutionary aspects as to where a free-living bacterium eventually acquired genetic material from environmental microorganisms, including flagellar genes, conferring on it the ability to reach other hosts (mammals), and, under selective pressure from the environment, can express these genes, helping it to evade the immune response. This review summarizes relevant aspects of the presence of flagellar proteins and puts into context their relevance in certain functions associated with the infective process. The study of these flagellar genes gives the genus Brucella a very high infectious versatility, placing it among the main organisms in urgent need of study, as it is linked to human health by direct contact with farm animals and by eventual transmission to the general population, where flagellar genes and proteins are of great relevance.

The one hundred year journey of the genus Brucella (Meyer and Shaw 1920)

The genus Brucella, described by Meyer and Shaw in 1920, comprises bacterial pathogens of veterinary and public health relevance. For 36 years, the genus came to include three species that caused brucellosis in livestock and humans. In the second half of the 20th century, bacteriologists discovered five new species and several 'atypical' strains in domestic animals and wildlife. In 1990, the Brucella species were recognized as part of the Class Alphaproteobacteria, clustering with pathogens and endosymbionts of animals and plants such as Bartonella, Agrobacterium and Ochrobactrum; all bacteria that live in close association with eukaryotic cells. Comparisons with Alphaproteobacteria contributed to identify virulence factors and to establish evolutionary relationships. Brucella members have two circular chromosomes, are devoid of plasmids, and display close genetic relatedness. A proposal, asserting that all brucellae belong to a single species with several subspecies debated for over 70 years, was ultimately rejected in 2006 by the subcommittee of taxonomy, based on scientific, practical, and biosafety considerations. Following this, the nomenclature of having multiples Brucella species prevailed and defined according to their molecular characteristics, host preference, and virulence. The 100-year history of the genus corresponds to the chronicle of scientific efforts and the struggle for understanding brucellosis.

Whole-Genome Sequencing for Tracing the Genetic Diversity of Brucella abortus and Brucella melitensis Isolated from Livestock in Egypt

Brucellosis is a highly contagious zoonosis that occurs worldwide. Whole-genome sequencing (WGS) has become a widely accepted molecular typing method for outbreak tracing and genomic epidemiology of brucellosis. Twenty-nine Brucella spp. (eight B. abortus biovar 1 and 21 B. melitensis biovar 3) were isolated from lymph nodes, milk, and fetal abomasal contents of infected cattle, buffaloes, sheep, and goats originating from nine districts in Egypt. The isolates were identified by microbiological methods and matrix-assisted laser desorption ionization time-of-flight mass spectrometry (MALDI-TOF MS). Differentiation and genotyping were confirmed using multiplex PCR. Illumina MiSeq^® was used to sequence the 29 Brucella isolates. Using MLST typing, ST11 and ST1 were identified among B. melitensis and B. abortus, respectively. Brucella abortus and B. melitensis isolates were divided into two main clusters (clusters 1 and 2) containing two and nine distinct genotypes by core-genome SNP analysis, respectively. The genotypes were irregularly distributed over time and space in the study area. Both Egyptian B. abortus and B. melitensis isolates proved to be genomically unique upon comparison with publicly available sequencing from strains of neighboring Mediterranean, African, and Asian countries. The antimicrobial resistance mechanism caused by mutations in rpoB, gyrA, and gyrB genes associated with rifampicin and ciprofloxacin resistance were identified. To the best of our knowledge, this is the first study investigating the epidemiology of Brucella isolates from livestock belonging to different localities in Egypt based on whole genome analysis.

Adhesive Functions or Pseudogenization of Type Va Autotransporters in Brucella Species

Adhesion to host cells is a key step for successful infection of many bacterial pathogens and may define tropism to different host tissues. To do so, bacteria display adhesins on their surfaces. Brucella is an intracellular pathogen capable of proliferating in a wide variety of cell types. It has been described that BmaC, a large protein that belongs to the classical (type Va) autotransporter family, is required for efficient adhesion of Brucella suis strain 1330 to epithelial cells and fibronectin. Here we show that B. suis 1330 harbors two other type Va autotransporters (BmaA and BmaB), which, although much smaller, share significant sequence similarities with BmaC and contain the essential domains to mediate proper protein translocation to the bacterial surface. Gain and loss of function studies indicated that BmaA, BmaB, and BmaC contribute, to a greater or lesser degree, to adhesion of B. suis 1330 to different cells such as synovial fibroblasts, osteoblasts, trophoblasts, and polarized epithelial cells as well as to extracellular matrix components. It was previously shown that BmaC localizes to a single bacterial pole. Interestingly, we observed here that, similar to BmaC, the BmaB adhesin is localized mostly at a single cell pole, reinforcing the hypothesis that Brucella displays an adhesive pole. Although Brucella species have strikingly similar genomes, they clearly differ in their host preferences. Mainly, the differences identified between species appear to be at loci encoding surface proteins. A careful in silico analysis of the putative type Va autotransporter orthologues from several Brucella strains showed that the bmaB locus from Brucella abortus and both, the bmaA and bmaC loci from Brucella melitensis are pseudogenes in all strains analyzed. Results reported here evidence that all three autotransporters play a role in the adhesion properties of B. suis 1330. However, Brucella spp. exhibit extensive variations in the repertoire of functional adhesins of the classical autotransporter family that can be displayed on the bacterial surface, making them an interesting target for future studies on host preference and tropism.

Emerging diversity and ongoing expansion of the genus Brucella

Remarkable genetic diversity and breadth of host species has been uncovered in the Brucella genus over the past decade, fundamentally changing our concept of what it means to be a Brucella. From ocean fishes and marine mammals, to pond dwelling amphibians, forest foxes, desert rodents, and cave-dwelling bats, Brucella have revealed a variety of previously unknown niches. Classical microbiological techniques have been able to help us classify many of these new strains but at times have limited our ability to see the true relationships among or within species. The closest relatives of Brucella are soil bacteria and the adaptations of Brucella spp. to live intracellularly suggest that the genus has evolved to live in vertebrate hosts. Several recently discovered species appear to have phenotypes that are intermediate between soil bacteria and core Brucella, suggesting that they may represent ancestral traits that were subsequently lost in the traditional species. Remarkably, the broad relationships among Brucella species using a variety of sequence and fragment-based approaches have been upheld when using comparative genomics with whole genomes. Nonetheless, genomes are required for fine-scale resolution of many of the relationships and for understanding the evolutionary history of the genus. We expect that the coming decades will reveal many more hosts and previously unknown diversity in a wide range of environments.

Canine brucellosis in Costa Rica reveals widespread Brucella canis infection and the recent introduction of foreign strains

Brucellosis is a prevalent disease in Costa Rica (CR), with an increasing number of human infections. Close to half of homes in CR have one or more dogs, corresponding to ∼1.4 million canines, most of them in the Central Valley within or near the cities of San José, Heredia, and Alajuela. From 302 dog sera collected from this region, 19 were positive for Brucella canis antigens, and five had antibodies against smooth lipopolysaccharide, suggesting infections by both B. canis and other Brucella species. B. canis strains were isolated in the Central Valley from 26 kennel dogs and three pet dogs, all displaying clinical signs of canine brucellosis. We detected three recent introductions of different B. canis strains in kennels: two traced from Mexico and one from Panama. Multiple locus-variable number tandem repeats (MLVA-16) and whole-genome sequencing (WGSA) analyses showed that B. canis CR strains comprise three main lineages. The tree topologies obtained by WGSA and MLVA-16 just partially agreed, indicating that the latter analysis is not suitable for phylogenetic studies. The fatty acid methyl ester analysis resolved five different B. canis groups, showing less resolution power than the MLVA-16 and WGSA. Lactobacillic acid was absent in linages I and II but present in linage III, supporting the recent introductions of B. canis strains from Mexico. B. canis displaying putative functional cyclopropane synthase for the synthesis of lactobacillic acid are phylogenetically intertwined with B. canis with non-functional protein, indicating that mutations have occurred independently in the various lineages.

Whole Genome Sequence Analysis of Brucella abortus Isolates from Various Regions of South Africa

The availability of whole genome sequences in public databases permits genome-wide comparative studies of various bacterial species. Whole genome sequence-single nucleotide polymorphisms (WGS-SNP) analysis has been used in recent studies and allows the discrimination of various Brucella species and strains. In the present study, 13 Brucella spp. strains from cattle of various locations in provinces of South Africa were typed and discriminated. WGS-SNP analysis indicated a maximum pairwise distance ranging from 4 to 77 single nucleotide polymorphisms (SNPs) between the South African Brucella abortus virulent field strains. Moreover, it was shown that the South African B. abortus strains grouped closely to B. abortus strains from Mozambique and Zimbabwe, as well as other Eurasian countries, such as Portugal and India. WGS-SNP analysis of South African B. abortus strains demonstrated that the same genotype circulated in one farm (Farm 1), whereas another farm (Farm 2) in the same province had two different genotypes. This indicated that brucellosis in South Africa spreads within the herd on some farms, whereas the introduction of infected animals is the mode of transmission on other farms. Three B. abortus vaccine S19 strains isolated from tissue and aborted material were identical, even though they originated from different herds and regions of South Africa. This might be due to the incorrect vaccination of animals older than the recommended age of 4-8 months or might be a problem associated with vaccine production.

Genomic insights into Brucella

Brucellosis is a zoonotic disease caused by certain species of Brucella. Each species has its preferred host animal, though it can infect other animals too. For a longer period, only six classical species were recognized in the genus Brucella. No vaccine is available for human brucellosis. Therefore, human brucellosis can be controlled only by controlling brucellosis in animals. The genus is now expanding with the newly isolated atypical strains from various animals, including marine mammals. Presently, 12 species of Brucella have been recognized. The first genome of Brucella was released in 2002, and today, we have more than 1500 genomes of Brucella spp. isolated worldwide. Multiple genome sequences are available for the major zoonotic species, B. abortus, B. melitensis, and B. suis. The Brucella genome has two chromosomes with the approximate sizes of 2.1 and 1.2 Mbp. The genome of Brucella is highly conserved across all the species at the nucleotide level. One of the unanswered questions is what makes host preference in different species of Brucella. Here, I summarize the recent advancements in the Brucella genomics research.

Non-adaptive Evolution of Trimeric Autotransporters in Brucellaceae

Brucella species are Gram-negative, facultative intracellular pathogens. They are the main cause of brucellosis, which has led to a global health burden. Adherence of the pathogen to the host cells is the first step in the infection process. The bacteria can adhere to various biotic and abiotic surfaces using their outer membrane proteins. Trimeric autotransporter adhesins (TAAs) are modular homotrimers of various length and domain complexity. They are a diverse, and widespread gene family constituting the type Vc secretion pathway. These adhesins have been established as virulence factors in Brucellaceae. To date, no comprehensive and exhaustive study has been performed on the trimeric autotransporter family in the genus. In the present study, various bioinformatics tools were used to provide a novel evolutionary insight into the sequence and structure of this protein family in Brucellaceae. To this end, a dataset of all trimeric autotransporters from the Brucella genomes was built. Analyses included but were not limited to sequence alignment, phylogenetic tree constructions, codon-based test for selection, clustering of the sequences, and structure (primary to quaternary) predictions. Batch analyzes of the dataset suggested the existence of a few structural domains within the whole population. BatA from the B. abortus 2308 genome was selected as a reference to describe the features of these structural domains. Furthermore, we examined the structural basis for the observed rigidity and resiliency of the protein structure through a molecular dynamics evaluation, which led us to deduce that the random drift results in the non-adaptive evolution of the trimeric autotransporter genes in the Brucella genus. Notably, the modifications have occurred across the genus without interference of gene transmission.

Brucella Genomics: Macro and Micro Evolution

Brucella organisms are responsible for one of the most widespread bacterial zoonoses, named brucellosis. The disease affects several species of animals, including humans. One of the most intriguing aspects of the brucellae is that the various species show a ~97% similarity at the genome level. Still, the distinct Brucella species display different host preferences, zoonotic risk, and virulence. After 133 years of research, there are many aspects of the Brucella biology that remain poorly understood, such as host adaptation and virulence mechanisms. A strategy to understand these characteristics focuses on the relationship between the genomic diversity and host preference of the various Brucella species. Pseudogenization, genome reduction, single nucleotide polymorphism variation, number of tandem repeats, and mobile genetic elements are unveiled markers for host adaptation and virulence. Understanding the mechanisms of genome variability in the Brucella genus is relevant to comprehend the emergence of pathogens.

First identification of a Brucella abortus biovar 4 strain from yak in Tibet, China

Brucellosis is one of the major zoonotic diseases in the world. In China, understanding on its causative agent Brucella is still limited. Recently, we isolated a Brucella strain XZ19-1 from yak in Lhasa, Tibet. Phenotypical characterization proved that it belongs to B. abortus biovar 4, a biotype that has never been reported in China. MLVA-16 genotyping revealed a novel profile (4-5-3-12-2-2-3-3-8-32-8-5-4-3-3-3) in this strain, while MLST sequence typing demonstrated that it belongs to ST 71. Furthermore, the whole genome of XZ19-1 strain was sequenced. Subsequent phylogenetic analysis demonstrated that XZ19-1was genetically more closely related to B. abortus strains originated from European countries rather than to those collected from China previously. Isolation and identification of XZ19-1 strain may thus indicate a unique Brucella lineage existing in Qing-Tibet plateau. These findings will help to improve the diagnosis and epidemiological studies of brucellosis in animals and human in this part of China.

ORF-based binarized structure network analysis of plasmids (OSNAp), a novel approach to core gene-independent plasmid phylogeny

OBJECTIVES

Systematic comparison of multiple plasmids remains challenging. We aimed to develop a new method for phylogenetic analysis of plasmids, open reading frame (ORF)-based binarized structure network analysis of plasmids (OSNAp).

METHODS

With the OSNAp, the genetic structures of plasmids in a given plasmid group are expressed as binary sequences based on the presence or absence of ORFs regardless of their positions or directions. As a proof-of-concept, ORFs were collected from 101 complete I1 plasmid sequences, and their corresponding binary sequences were generated. A tree was generated using the neighbor-net, an algorithm for constructing phylogenetic networks based on distance between taxa, to visualize the plasmid phylogeny drawn from binary sequences. The results were compared with those of plasmid sequence types (pSTs) defined by plasmid multilocus sequence typing (pMLST).

RESULTS

All I1 plasmids were placed on the phylogenetic tree constructed from the binary sequences. Most plasmids belonging to the same pSTs had Dice indices of ≥0.95 and were placed in the same OSNAp split. On the other hand, pST12 plasmids were distributed on separate splits due to differences in ORFs not used in pMLST, suggesting improved differentiation of the plasmids with OSNAp compared with pMLST.

CONCLUSION

OSNAp is a novel holistic approach to assess relatedness of a population of plasmids in a given plasmid group based on nucleotide sequence data. It provides higher discrimination than pMLST, which may prove useful in tracing bacteria that harbor plasmids of shared origins.

Taxonomic Organization of the Family Brucellaceae Based on a Phylogenomic Approach

Deciphering the evolutionary history of pathogenic bacteria and their near neighbors may help to understand the genetic or ecological bases which led to their pathogenic behavior. The Brucellaceae family comprises zoonotic pathogenic species belonging to the genus Brucella as well as the environmental genus Ochrobactrum for which some species are considered as opportunistic pathogens. Here, we used a phylogenomic approach including a set of 145 Brucellaceae genomes representative of the family diversity and more than 40 genomes of the order Rhizobiales to infer the taxonomic relationships between the family’s species. Our results clarified some unresolved phylogenetic ambiguities, conducting to the exclusion of Mycoplana spp. out of the family Brucellaceae and the positioning of all Brucella spp. as a single genomic species within the current Ochrobactrum species diversity. Additional analyses also revealed that Ochrobactrum spp. separate into two clades, one comprising mostly environmental species while the other one includes the species considered as pathogens (Brucella spp.) or opportunistic pathogens (mainly O. anthropi, O. intermedium, and O. pseudintermedium). Finally, we show that O. intermedium is undergoing a beginning of genome reduction suggestive of an ongoing ecological niche specialization, and that some lineages of O. intermedium and O. anthropi may shift toward an adaption to the human host.

Polymorphisms in Brucella Carbonic Anhydrase II Mediate CO2 Dependence and Fitness in vivo

Some Brucella isolates are known to require an increased concentration of CO₂ for growth, especially in the case of primary cultures obtained directly from infected animals. Moreover, the different Brucella species and biovars show a characteristic pattern of CO₂ requirement, and this trait has been included among the routine typing tests used for species and biovar differentiation. By comparing the differences in gene content among different CO₂-dependent and CO₂-independent Brucella strains, we have confirmed that carbonic anhydrase (CA) II is the enzyme responsible for this phenotype in all the Brucella strains tested. Brucella species contain two CAs of the β family, CA I and CA II; genetic polymorphisms exist for both of them in different isolates, but only those putatively affecting the activity of CA II correlate with the CO₂ requirement of the corresponding isolate. Analysis of these polymorphisms does not allow the determination of CA I functionality, while the polymorphisms in CA II consist of small deletions that cause a frameshift that changes the C-terminus of the protein, probably affecting its dimerization status, essential for the activity. CO₂-independent mutants arise easily in vitro, although with a low frequency ranging from 10^–6 to 10^–10 depending on the strain. These mutants carry compensatory mutations that produce a full-length CA II. At the same time, no change was observed in the sequence coding for CA I. A competitive index assay designed to evaluate the fitness of a CO₂-dependent strain compared to its corresponding CO₂-independent strain revealed that while there is no significant difference when the bacteria are grown in culture plates, growth in vivo in a mouse model of infection provides a significant advantage to the CO₂-dependent strain. This could explain why some Brucella isolates are CO₂ dependent in primary isolation. The polymorphism described here also allows the in silico determination of the CO₂ requirement status of any Brucella strain.

Brucella: not your ‘typical’ intracellular pathogen

Currently the genus Brucella consists of a group of bacteria that are genetically monospecific yet phenotypically diverse, and a recent genetic and phenotypic divergent group known as ‘atypical' Brucellae. The host range is extremely varied and includes mammals, including humans, terrestrial animals and marine mammals, but now extends to reptiles and amphibians. Almost all Brucella species are zoonotic. The disease collectively termed Brucellosis leads to abortion and reproductive disease in animals, whereas human infection presents as a non-specific undulating fever accompanied by general malaise, chills, joint pain, muscle aches, genitourinary disease and adverse pregnancy outcomes. These Gram-negative coccobacilli invade and replicate in the host macrophages where they can limit the effects of the host immune system and antibiotic treatment. Due to the phenotypic and genotypic diversity and close relationship with Ochrobactrum species, the genus Brucella presents challenges for accurate identification and recognition of new species.

Laboratory Diagnosis of Human Brucellosis

The clinical presentation of brucellosis in humans is variable and unspecific, and thus, laboratory corroboration of the diagnosis is essential for the patient's proper treatment. The diagnosis of brucellar infections can be made by culture, serological tests, and nucleic acid amplification assays. Modern automated blood culture systems enable detection of acute cases of brucellosis within the routine 5- to 7-day incubation protocol employed in clinical microbiology laboratories, although a longer incubation and performance of blind subcultures may be needed for protracted cases. Serological tests, though they lack specificity and provide results that may be difficult to interpret in individuals repeatedly exposed to Brucella organisms, nevertheless remain a diagnostic cornerstone in resource-poor countries. Nucleic acid amplification assays combine exquisite sensitivity, specificity, and safety and enable rapid diagnosis of the disease. However, long-term persistence of positive molecular test results in patients that have apparently fully recovered is common and has unclear clinical significance and therapeutic implications. Therefore, as long as there are no sufficiently validated commercial tests or studies that demonstrate an adequate interlaboratory reproducibility of the different homemade PCR assays, cultures and serological methods will remain the primary tools for the diagnosis and posttherapeutic follow-up of human brucellosis.

Molecular characterization of Brucella species from Zimbabwe

Brucella abortus and B. melitensis have been reported in several studies in animals in Zimbabwe but the extent of the disease remains poorly known. Thus, characterizing the circulating strains is a critical first step in understanding brucellosis in the country. In this study we used an array of molecular assays including AMOS-PCR, Bruce-ladder, multiple locus variable number tandem repeats analysis (MLVA) and single nucleotide polymorphisms from whole genome sequencing (WGS-SNP) to characterize Brucella isolates to the species, biovar, and individual strain level. Sixteen Brucella strains isolated in Zimbabwe at the Central Veterinary laboratory from various hosts were characterized using all or some of these assays. The strains were identified as B. ovis, B. abortus, B. canis and B. suis, with B. canis being the first report of this species in Zimbabwe. Zimbabwean strains identified as B. suis and B. abortus were further characterized with whole genome sequencing and were closely related to reference strains 1330 and 86/8/59, respectively. We demonstrate the range of different tests that can be performed from simple assays that can be run in laboratories lacking sophisticated instrumentation to whole genome analyses that currently require substantial expertise and infrastructure often not available in the developing world.

Comparative Ecology of Bartonella and Brucella Infections in Wild Carnivores

Phylogenetic sister clades Bartonella and Brucella within the order Rhizobiales present some common biological characteristics as well as evident differences in adaptations to their mammalian reservoirs. We reviewed published data on Bartonella and Brucella infections in wild carnivores to compare the ecology of these bacteria in relatively similar host environments. Arthropod vectors are the main mechanism for Bartonella species transmission between mammalian hosts. The role of arthropods in transmission of Brucella remains disputed, however experimental studies and reported detection of Brucella in arthropods indicate potential vector transmission. More commonly, transmission of Brucella occurs via contact exposure to infected animals or the environment contaminated with their discharges. Of 26 species of carnivores tested for both Bartonella and Brucella, 58% harbored either. Among them were bobcats, African lions, golden jackals, coyotes, wolves, foxes, striped skunks, sea otters, raccoons, and harbor seals. The most common species of Bartonella in wild carnivores was B. henselae, found in 23 species, followed by B. rochalimae in 12, B. clarridgeiae in ten, and B. vinsonii subsp. berkhoffii in seven. Among Brucella species, Br. abortus was reported in over 30 terrestrial carnivore species, followed by Br. canis in seven. Marine carnivores, such as seals and sea lions, can host Br. pinnipedialis. In contrast, there is no evidence of a Bartonella strain specific for marine mammals. Bartonella species are present practically in every sampled species of wild felids, but of 14 Brucella studies of felids, only five reported Brucella and those were limited to detection of antibodies. We found no reports of Bartonella in bears while Brucella was detected in these animals. There is evident host-specificity of Bartonella species in wild carnivores (e.g., B. henselae in felids and B. vinsonii subsp. berkhoffii in canids). A co-adaptation of Brucella with terrestrial wild carnivore hosts is not as straightforward as in domestic animals. Wild carnivores often carry the same pathogens as their domesticated relatives (cats and dogs), but the risk of exposure varies widely because of differences in biology, distribution, and historical interactions.

Learning from the master: targets and functions of the CtrA response regulator in Brucella abortus and other alpha-proteobacteria

The α-proteobacteria are a fascinating group of free-living, symbiotic and pathogenic organisms, including the Brucella genus, which is responsible for a worldwide zoonosis. One common feature of α-proteobacteria is the presence of a conserved response regulator called CtrA, first described in the model bacterium Caulobacter crescentus, where it controls gene expression at different stages of the cell cycle. Here, we focus on Brucella abortus and other intracellular α-proteobacteria in order to better assess the potential role of CtrA in the infectious context. Comparative genomic analyses of the CtrA control pathway revealed the conservation of specific modules, as well as the acquisition of new factors during evolution. The comparison of CtrA regulons also suggests that specific clades of α-proteobacteria acquired distinct functions under its control, depending on the essentiality of the transcription factor. Other CtrA-controlled functions, for instance motility and DNA repair, are proposed to be more ancestral. Altogether, these analyses provide an interesting example of the plasticity of a regulation network, subject to the constraints of inherent imperatives such as cell division and the adaptations to diversified environmental niches.

On the reconstruction of the ancestral bacterial genomes in genus Mycobacterium and Brucella

BACKGROUND

To reconstruct the evolution history of DNA sequences, novel models of increasing complexity regarding the number of free parameters taken into account in the sequence evolution, as well as faster and more accurate algorithms, and statistical and computational methods, are needed. More particularly, as the principal forces that have led to major structural changes are genome rearrangements (such as translocations, fusions, and so on), understanding their underlying mechanisms, among other things via the ancestral genome reconstruction, are essential. In this problem, since finding the ancestral genomes that minimize the number of rearrangements in a phylogenetic tree is known to be NP-hard for three or more genomes, heuristics are commonly chosen to obtain approximations of the exact solution. The aim of this work is to show that another path is possible.

RESULTS

Various algorithms and software already deal with the difficult nature of the problem of reconstruction of the ancestral genome, but they do not function with precision, in particular when indels or single nucleotide polymorphisms fall into repeated sequences. In this article, and despite the theoretical NP-hardness of the ancestral reconstruction problem, we show that an exact solution can be found in practice in various cases, encompassing organelles and some bacteria. A practical example proves that an accurate reconstruction, which also allows to highlight homoplasic events, can be obtained. This is illustrated by the reconstruction of ancestral genomes of two bacterial pathogens, belonging in Mycobacterium and Brucella genera.

CONCLUSIONS

By putting together automatically reconstructed ancestral regions with handmade ones for problematic cases, we show that an accurate reconstruction of ancestors of the Brucella genus and of the Mycobacterium tuberculosis complex is possible. By doing so, we are able to investigate the evolutionary history of each pathogen by computing their common ancestors. They can be investigated extensively, by studying the gene content evolution over time, the resistance acquisition, and the impacts of mobile elements on genome plasticity.

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

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.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5001-6) contains supplementary material, which is available to authorized users.

Core Genome Multilocus Sequence Typing and Single Nucleotide Polymorphism Analysis in the Epidemiology of Brucella melitensis Infections

The use of whole-genome sequencing (WGS) using next-generation sequencing (NGS) technology has become a widely accepted method for microbiology laboratories in the application of molecular typing for outbreak tracing and genomic epidemiology. Several studies demonstrated the usefulness of WGS data analysis through single-nucleotide polymorphism (SNP) calling from a reference sequence analysis for Brucella melitensis, whereas gene-by-gene comparison through core-genome multilocus sequence typing (cgMLST) has not been explored so far. The current study developed an allele-based cgMLST method and compared its performance to that of the genome-wide SNP approach and the traditional multilocus variable-number tandem repeat analysis (MLVA) on a defined sample collection. The data set was comprised of 37 epidemiologically linked animal cases of brucellosis as well as 71 isolates with unknown epidemiological status, composed of human and animal samples collected in Italy. The cgMLST scheme generated in this study contained 2,704 targets of the B. melitensis 16M reference genome. We established the potential criteria necessary for inclusion of an isolate into a brucellosis outbreak cluster to be ≤6 loci in the cgMLST and ≤7 in WGS SNP analysis. Higher phylogenetic distance resolution was achieved with cgMLST and SNP analysis than with MLVA, particularly for strains belonging to the same lineage, thereby allowing diverse and unrelated genotypes to be identified with greater confidence. The application of a cgMLST scheme to the characterization of B. melitensis strains provided insights into the epidemiology of this pathogen, and it is a candidate to be a benchmark tool for outbreak investigations in human and animal brucellosis.

Correct Identification of Ochrobactrum anthropi From Blood Culture Using 16rRNA Sequencing: A First Case Report in an Immunocompromised Patient in Mexico

The present report describes the misidentification of Brucella spp. from a positive blood culture using traditional microbiology tests. A molecular test identified the bacterium as Ochrobactrum anthropi. According to the information available, this report is the first to include this type of case in Mexico.

Global evolution and phylogeography of Brucella melitensis strains

Background

Brucellosis is a bacterial zoonotic disease. Annually in the world more than 500,000 new cases of brucellosis in humans are registered. In this study, we propose an evolutionary model of the historical distribution of B. melitensis based on the full-genomic SNP analysis of 98 strains.

Results

We performed an analysis of the SNP of the complete genomes of 98 B. melitensis strains isolated in different geographical regions of the world to obtain relevant information on the population structure, genetic diversity and the evolution history of the species. Using genomic sequences of 21 strains of B. melitensis isolated in Russia and WGS data from the NCBI database, it was possible to identify five main genotypes and 13 species genotypes for analysis. Data analysis based on the Bayesian Phylogenetics and Phylogeography method allowed to determine the regions of geographical origin and the expected pathways of distribution of the main lines (genotypes and subgenotypes) of the pathogen.

Conclusions

Within the framework of our study, the model of global evolution and phylogeography of B. melitensis strains isolated in various regions of the planet was proposed for the first time. The sets of unique specific SNPs described in our study, for all identified genotypes and subgenotypes, can be used to develop new bacterial typing and identification systems for B. melitensis.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4762-2) contains supplementary material, which is available to authorized users.

Typing and comparative genome analysis of Brucella melitensis isolated from Lebanon

Brucella melitensis is the main causative agent of the zoonotic disease brucellosis. This study aimed at typing and characterizing genetic variation in 33 Brucella isolates recovered from patients in Lebanon. Bruce-ladder multiplex PCR and PCR-RFLP of omp31, omp2a and omp2b were performed. Sixteen representative isolates were chosen for draft-genome sequencing and analyzed to determine variations in virulence, resistance, genomic islands, prophages and insertion sequences. Comparative whole-genome single nucleotide polymorphism analysis was also performed. The isolates were confirmed to be B. melitensis. Genome analysis revealed multiple virulence determinants and efflux pumps. Genome comparisons and single nucleotide polymorphisms divided the isolates based on geographical distribution but revealed high levels of similarity between the strains. Sequence divergence in B. melitensis was mainly due to lateral gene transfer of mobile elements. This is the first report of an in-depth genomic characterization of B. melitensis in Lebanon.

Brucellosis: Evolution and expected comeback

Brucellosis is a serious infectious disease which causes great direct and indirect economic loses for animal holders worldwide such as the reduction of milk and meat production through abortions/culling of positive reactors, the expense of disease control/eradication and farmers compensation. Although the disease was eradicated from most of the industrial countries, it remains one of the most common zoonotic diseases in developing countries being responsible for more than 500,000 new cases yearly. Brucella is considered to be a bioterrorism organism due to its low infectious doses (10-100 bacteria), capability of persistence in the environment, rapid transmission via different routes including aerosols, and finally due to its difficult treatment by antibiotics.There are many reasons to believe that a new comeback of brucellosis may occur in near future. This expectation is supported by the recent discovery of new atypical Brucella species with new genetic properties and the recent reports of (man to man) disease transmission as will be discussed later. The development of new concepts and measurements for disease control is urgently required. In the present review, the evolution of Brucella and the different factors favoring its comeback are discussed.

[Characterization of the genetic variability of field strains of Brucella canis isolated in Antioquia]

Brucella canis is a facultative intracellular pathogen responsible for canine brucellosis, a zoonotic disease that affects canines, causing abortions and reproductive failure; and the production of non-specific symptoms in humans. In 2005 the presence of B. canis in Antioquia was demonstrated and the strains were identified as type 2. The sequencing of the genome of a field strain denoted Brucella canis str. Oliveri, showed species-specific indel events, which led us to investigate the genomic characteristics of the B. canis strain isolated and to establish the phylogenetic relationships and the divergence time of B. canis str. Oliveri. Conventional PCR sequencing was performed in 30 field strains identifying 5 indel events recognized in B. canis str. Oliveri. ADN from Brucella suis, Brucella melitensis and vaccine strains from Brucella abortus were used as control, and it was determined that all of the studied field strains shared 4 out of the 5 indels of the sequenced Oliveri strain, indicating the presence of more than one strain circulating in the region. Phylogenetic analysis was performed with 24 strains of Brucella using concatenated sequences of genetic markers for species differentiation. The molecular clock hypothesis and Tajima's relative rate test were tested, showing that the Oliveri strain, similarly to other canis species, diverged from B. suis. The molecular clock hypothesis between Brucella species was rejected and an evolution rate and a similar genetic distance between the B. canis were demonstrated.

Genotyping of Brucella melitensis and Brucella abortus strains in Kazakhstan using MLVA-15

Currently, although the prevalence of brucellosis in Kazakhstan remains high, there are limited data available on the genetic diversity of circulating Brucella strains. Here, MLVA was employed to genotype a panel of 102 Brucella isolates collected from eight Kazakh regions and neighboring countries (Russia, Kyrgyzstan) during the period 1935-2017. MLVA-11 analysis classified 64 B. abortus strains into genotypes 72, 82, 331, 71, 341 and 69, while one genotype was novel, having no correspondence within the MLVA international database. MLVA-11 analysis of 37 B. melitensis strains showed 100% identity with genotypes 116, 114 and 11. One B. suis strain was classified into genotype 33. Phylogeography based on MLVA-15 demonstrated that all B. abortus and B. melitensis strains belonged to "Abortus C" and "East Mediterranean" lineages, respectively. B. abortus strains from Kazakhstan and Russia resulted genetically related to Portuguese, Brazilian and US isolates, suggesting ancient spread of these lineages from Europe westwards to South America and eastwards to Turkey, Russia and Asia. Most of Kazakh B. melitensis isolates were related to strains circulating in China, likely due to long-term trading partnerships between the two countries. In fine-scale MLVA-15 analysis, 17 B. abortus and 12 B. melitensis genotypes were identified; among them 12 are novel. Interestingly, epidemiological information supporting molecular data were retrieved for two clusters within the B. abortus group, thus proving that MLVA is an appropriate tool for effective traceback analyses. Our findings suggest that molecular genotyping should be applied systematically to support control plans for eradication of brucellosis in Kazakhstan.

Identification of genetic variants of Brucella spp. through genome-wide association studies

Brucellosis is an important zoonotic disease caused by Brucella spp. We present a phylogeny of 552 strains based on genome-wide single nucleotide polymorphisms (SNPs) determined by an alignment-free k-mer approach. A total of 138,029 SNPs were identified from 552 Brucella genomes. Of these, 31,152 and 106,877 were core and non-core SNPs, respectively. Based on pan-genome analysis 11,937 and 972 genes were identified as pan and core genome, respectively. The pan-genome-wide analysis studies (Pan-GWAS) could not identify the group-specific variants in Brucella spp. Therefore, we focused on SNP based genome-wide association studies (SNP-GWAS) to identify the species-specific genetic determinants in Brucella spp. Phylogenetic tree representing eleven recognized Brucella spp. showed 16 major lineages. We identified 143 species-specific SNPs in Brucella abortus that are conserved in 311 B. abortus genomes. Of these, 141 species-specific SNPs were confined in the positively significant SNPs of B. abortus using SNP-GWAS. Since conserved in all the B. abortus genomes studied, these SNPs might have originated very early during the evolution of B. abortus and might be responsible for the evolution of B. abortus with cattle as the preferred host. Similarly, we identified 383 species-specific SNPs conserved in 132 Brucella melitensis genomes. Of these 379 species-specific SNPs were identified as positively associated using GWAS. Interestingly, >98% of the SNPs that are significantly, positively associated with the traits showed 100% sensitivity and 100% specificity. These identified species-specific core-SNPs identified in Brucella genomes could be responsible for the speciation and their respective host adaptation.

Evolution and genome specialization of Brucella suis biovar 2 Iberian lineages

Background

Swine brucellosis caused by B. suis biovar 2 is an emergent disease in domestic pigs in Europe. The emergence of this pathogen has been linked to the increase of extensive pig farms and the high density of infected wild boars (Sus scrofa). In Portugal and Spain, the majority of strains share specific molecular characteristics, which allowed establishing an Iberian clonal lineage. However, several strains isolated from wild boars in the North-East region of Spain are similar to strains isolated in different Central European countries.

Results

Comparative analysis of five newly fully sequenced B. suis biovar 2 strains belonging to the main circulating clones in Iberian Peninsula, with publicly available Brucella spp. genomes, revealed that strains from Iberian clonal lineage share 74% similarity with those reference genomes. Besides the 210 kb translocation event present in all biovar 2 strains, an inversion with 944 kb was presented in chromosome I of strains from the Iberian clone. At left and right crossover points, the inversion disrupted a TRAP dicarboxylate transporter, DctM subunit, and an integral membrane protein TerC. The gene dctM is well conserved in Brucella spp. except in strains from the Iberian clonal lineage. Intraspecies comparative analysis also exposed a number of biovar-, haplotype- and strain-specific insertion-deletion (INDELs) events and single nucleotide polymorphisms (SNPs) that could explain differences in virulence and host specificities. Most discriminative mutations were associated to membrane related molecules (29%) and enzymes involved in catabolism processes (20%). Molecular identification of both B. suis biovar 2 clonal lineages could be easily achieved using the target-PCR procedures established in this work for the evaluated INDELs.

Conclusion

Whole-genome analyses supports that the B. suis biovar 2 Iberian clonal lineage evolved from the Central-European lineage and suggests that the genomic specialization of this pathogen in the Iberian Peninsula is independent of a specific genomic event(s), but instead driven by allopatric speciation, resulting in the establishment of a new ecovar.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4113-8) contains supplementary material, which is available to authorized users.

Caprine brucellosis: A historically neglected disease with significant impact on public health

Caprine brucellosis is a chronic infectious disease caused by the gram-negative cocci-bacillus Brucella melitensis. Middle- to late-term abortion, stillbirths, and the delivery of weak offspring are the characteristic clinical signs of the disease that is associated with an extensive negative impact in a flock's productivity. B. melitensis is also the most virulent Brucella species for humans, responsible for a severely debilitating and disabling illness that results in high morbidity with intermittent fever, chills, sweats, weakness, myalgia, abortion, osteoarticular complications, endocarditis, depression, anorexia, and low mortality. Historical observations indicate that goats have been the hosts of B. melitensis for centuries; but around 1905, the Greek physician Themistokles Zammit was able to build the epidemiological link between "Malta fever" and the consumption of goat milk. While the disease has been successfully managed in most industrialized countries, it remains a significant burden on goat and human health in the Mediterranean region, the Middle East, Central and Southeast Asia (including India and China), sub-Saharan Africa, and certain areas in Latin America, where approximately 3.5 billion people live at risk. In this review, we describe a historical evolution of the disease, highlight the current worldwide distribution, and estimate (by simple formula) the approximate costs of brucellosis outbreaks to meat- and milk-producing farms and the economic losses associated with the disease in humans. Successful control leading to eradication of caprine brucellosis in the developing world will require a coordinated Global One Health approach involving active involvement of human and animal health efforts to enhance public health and improve livestock productivity.

Fast and accurate phylogeny reconstruction using filtered spaced-word matches

Motivation

Word-based or ‘alignment-free’ algorithms are increasingly used for phylogeny reconstruction and genome comparison, since they are much faster than traditional approaches that are based on full sequence alignments. Existing alignment-free programs, however, are less accurate than alignment-based methods.

Results

We propose Filtered Spaced Word Matches (FSWM), a fast alignment-free approach to estimate phylogenetic distances between large genomic sequences. For a pre-defined binary pattern of match and don’t-care positions, FSWM rapidly identifies spaced word-matches between input sequences, i.e. gap-free local alignments with matching nucleotides at the match positions and with mismatches allowed at the don’t-care positions. We then estimate the number of nucleotide substitutions per site by considering the nucleotides aligned at the don’t-care positions of the identified spaced-word matches. To reduce the noise from spurious random matches, we use a filtering procedure where we discard all spaced-word matches for which the overall similarity between the aligned segments is below a threshold. We show that our approach can accurately estimate substitution frequencies even for distantly related sequences that cannot be analyzed with existing alignment-free methods; phylogenetic trees constructed with FSWM distances are of high quality. A program run on a pair of eukaryotic genomes of a few hundred Mb each takes a few minutes.

Availability and Implementation

The program source code for FSWM including a documentation, as well as the software that we used to generate artificial genome sequences are freely available at http://fswm.gobics.de/

Supplementary information

Supplementary data are available at Bioinformatics online.