Genomic comparisons of Brucella spp. and closely related bacteria using base compositional and proteome based methods

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.

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.

Selection of reference genes for measuring the expression of aiiO in Ochrobactrum quorumnocens A44 using RT-qPCR

Reverse transcription quantitative PCR (RT-qPCR), a method of choice for quantification of gene expression changes, requires stably expressed reference genes for normalization of data. So far, no reference genes were established for the Alphaproteobacteria of the genus Ochrobactrum. Here, we determined reference genes for gene expression studies in O. quorumnocens A44. Strain A44 was cultured under 10 different conditions and the stability of expression of 11 candidate genes was evaluated using geNorm, NormFinder and BestKeeper. Most stably expressed genes were found to be rho, gyrB and rpoD. Our results can facilitate the choice of reference genes in the related Ochrobactrum strains. O. quorumnocens A44 is able to inactivate a broad spectrum of N-acyl homoserine lactones (AHLs) - the quorum sensing molecules of many Gram-negative bacteria. This activity is attributed to AiiO hydrolase, yet it remains unclear whether AHLs are the primary substrate of this enzyme. Using the established RT-qPCR setup, we found that the expression of the aiiO gene upon exposure to two AHLs, C6-HLS and 3OC12-HSL, does not change above the 1-fold significance threshold. The implications of this finding are discussed in the light of the role of quorum sensing-interfering enzymes in the host strains.

Intracellular invasion and survival of Brucella neotomae, another possible zoonotic Brucella species

In 1967, Brucella neotomae was first isolated from Neotoma lepida, the dessert wood rat, in Utah. With little infection data since its discovery, the zoonotic potential of this Brucella species is largely unknown. Recent reports of isolation from human cerebrospinal fluid, along with current literature suggest that B. neotomae has the ability to infect various hosts and cell types. In this report we extend the knowledge of B. neotomae ATCC 23459’s intracellular invasion and survival abilities to a variety of cell lines through gentamicin protection assays. Some of the phagocytic and epithelial cell lines from various mammalian species represent characteristics of some cell types that could be encountered by Brucella in potential hosts. It was found that B. neotomae ATCC 23459 exhibits generally lower intracellular bacterial CFUs compared to the mouse-passaged strain of B. neotomae ATCC 23459, B. suis 1330, and B. abortus 2308. Ultimately, these observations provide a small piece of the puzzle in the investigation of the breadth of B. neotomae’s pathogenic potential.

Ochrobactrum quorumnocens sp. nov., a quorum quenching bacterium from the potato rhizosphere, and comparative genome analysis with related type strains

Ochrobactrum spp. are ubiquitous bacteria attracting growing attention as important members of microbiomes of plants and nematodes and as a source of enzymes for biotechnology. Strain Ochrobactrum sp. A44^T was isolated from the rhizosphere of a field-grown potato in Gelderland, the Netherlands. The strain can interfere with quorum sensing (QS) of Gram-negative bacteria through inactivation of N-acyl homoserine lactones (AHLs) and protect plant tissue against soft rot pathogens, the virulence of which is governed by QS. Phylogenetic analysis based on 16S rRNA gene alone and concatenation of 16S rRNA gene and MLSA genes (groEL and gyrB) revealed that the closest relatives of A44^T are O. grignonense OgA9a^T, O. thiophenivorans DSM 7216^T, O. pseudogrignonense CCUG 30717^T, O. pituitosum CCUG 50899^T, and O. rhizosphaerae PR17^T. Genomes of all six type strains were sequenced, significantly expanding the possibility of genome-based analyses in Ochrobactrum spp. Average nucleotide identity (ANIb) and genome-to-genome distance (GGDC) values for A44^T and the related strains were below the single species thresholds (95% and 70%, respectively), with the highest scores obtained for O. pituitosum CCUG 50899^T (87.31%; 35.6%), O. rhizosphaerae PR17^T (86.80%; 34.3%), and O. grignonense OgA9a^T (86.30%; 33.6%). Distinction of A44^T from the related type strains was supported by chemotaxonomic and biochemical analyses. Comparative genomics revealed that the core genome for the newly sequenced strains comprises 2731 genes, constituting 50–66% of each individual genome. Through phenotype-to-genotype study, we found that the non-motile strain O. thiophenivorans DSM 7216^T lacks a cluster of genes related to flagella formation. Moreover, we explored the genetic background of distinct urease activity among the strains. Here, we propose to establish a novel species Ochrobactrum quorumnocens, with A44^T as the type strain (= LMG 30544^T = PCM 2957^T).

Identification of Brucella abortus biovar 4 of bovine origin in Colombia

The objective of this study was to identify twelve Brucella abortus isolates of bovine origin from the department of Nariño in Colombia up to the biovar level. These isolates are included in the collection of the Germplasm Bank of Microorganisms of Animal Health Interest - Bacteria and Virus (BGSA-BV). The identification was carried out through conventional methods such as macro and microscopic morphological descriptions, enzymatic activity, biochemical profile, substrate use and sensitivity to dyes. Complementary genotypic characterization was carried out using multiplex PCR for B. abortus, Brucella melitensis, Brucella ovis, and Brucella suis-Erytritol (AMOS-ERY-PCR), RFLP-IS711, by southern blot hybridization, as well as by the multiple locus variable number of tandem repeat analysis (MLVA) using the ery gene and the insertion sequence IS711 and variable number of tandem repeats (VNTR) as molecular markers. The results of the phenotypic and molecular characterization allowed to identify twelve isolates as B. abortus biovar 4 as well as to differentiate field from vaccine strains. This is the first study on the phenotypic and molecular identification of B. abortus isolates in Colombia. It was concluded that the phenotypic and molecular identification of twelve isolates as B. abortus biovar 4 could be achieved using conventional and molecular techniques with enough resolution power. The identification of these isolates to the biovar level in taxonomic and epidemiological terms will allow the use of this genetic resource as reference strains in future research. This finding constitutes the basis for identifying biotypes not previously reported in the country that might be useful to support brucellosis survey programs in Colombia.

Whole-genome sequencing and antimicrobial resistance in Brucella melitensis from a Norwegian perspective

Brucellosis is a rarely encountered infection in Norway. The aim of this study was to explore all Brucella melitensis isolates collected in Norway from 1999 to 2016 in relation to origin of infection and antimicrobial resistance patterns. A total of 23 isolates were analysed by whole-genome sequencing and compared with selected sequences of B. melitensis available from NCBI. Additionally, SNP analysis in antibiotic resistance determining genes was performed. The majority belonged to the East Mediterranean clade (genotype II), while the remaining isolates belonged to the African clade (genotype III). These results indicate that human brucellosis in Norway is related to travels or migration from the Middle East, Asia or Africa, in accordance with results from Germany, Denmark and Sweden. Antibiotic susceptibility patterns were determined by broth microdilution method and/or gradient strip method. All isolates were susceptible for all tested antibiotics, except for rifampicin where phenotypical results indicated resistance or intermediate resistance in all isolates based on broth microdilution method, and in four isolates based on gradient strip testing. In contrast, screening of the rpoB gene did not reveal any mutations in the previously described rpoB "hot spot" regions related to rifampicin resistance, indicating overestimation of resistance based on phenotypical results.

Application of real-time quantitative PCR assays for detecting marine Brucella spp. in fish

Brucella ceti and Brucella pinnipedialis have been documented as occurring in marine mammals, and B. ceti has been identified in 3 naturally acquired human cases. Seroconversion and infection patterns in Pacific Northwest harbor seals ( Phoca vitulina richardii) and North Atlantic hooded seals ( Cystophora cristata) indicate post-weaning exposure through prey consumption or lungworm infection, suggesting fish and possibly invertebrates play an epizootiologic role in marine Brucella transmission and possible foodborne risk to humans. We determined if real-time quantitative PCR (qPCR) assays can detect marine Brucella DNA in fish DNA. Insertion sequence (IS) 711 gene and sequence type (ST)27 primer-probe sets were used to detect Brucella associated with marine mammals and human zoonotic infections, respectively. First, DNA extracts from paired-species fish (containing 2 species) samples were tested and determined to be Brucella DNA negative using both IS 711 and ST27 primer-probe sets. A representative paired-species fish DNA sample was spiked with decreasing concentrations of B. pinnipedialis DNA to verify Brucella detection by the IS 711 primer-probe within fish DNA. A standard curve, developed using isolated DNA from B. pinnipedialis, determined the limit of detection. Finally, the IS 711 primer-probe was used to test Atlantic cod ( Gadus morhua) DNA extracts experimentally infected with the B. pinnipedialis hooded seal strain. In culture-positive cod tissue, the IS 711 limit of detection was ~1 genome copy of Brucella. Agreement between culture and PCR results for the 9 positive and 9 negative cod tissues was 100%. Although a larger sample set is required for validation, our study shows that qPCR can detect marine Brucella in fish.

Genome expansion in bacteria: the curios case of Chlamydia trachomatis

Background

Recent findings indicated that a correlation between genomic % AT and genome size within strains of microbial species was predominantly associated with the uptake of foreign DNA. One species however, Chlamydia trachomatis, defied any explanation. In the present study 79 fully sequenced C. trachomatis genomes, representing ocular- (nine strains), urogenital- (36 strains) and lymphogranuloma venereum strains (LGV, 22 strains), in three pathogroups, in addition to 12 laboratory isolates, were scrutinized with the intent of elucidating the positive correlation between genomic AT content and genome size.

Results

The average size difference between the strains of each pathogroup was largely explained by the incorporation of genetic fragments. These fragments were slightly more AT rich than their corresponding host genomes, but not enough to justify the difference in AT content between the strains of the smaller genomes lacking the fragments. In addition, a genetic region predominantly found in the ocular strains, which had the largest genomes, was on average more GC rich than the host genomes of the urogenital strains (58.64 % AT vs. 58.69 % AT), which had the second largest genomes, implying that the foreign genetic regions cannot alone explain the association between genome size and AT content in C. trachomatis. 23,492 SNPs were identified for all 79 genomes, and although the SNPs were on average slightly GC rich (~47 % AT), a significant association was found between genome-wide SNP AT content, for each pathogroup, and genome size (p < 0.001, R² = 0.86) in the C. trachomatis strains.

Conclusions

The correlation between genome size and AT content, with respect to the C. trachomatis pathogroups, was explained by the incorporation of genetic fragments unique to the ocular and/or urogenital strains into the LGV- and urogential strains in addition to the genome-wide SNP AT content differences between the three pathogroups.

Electronic supplementary material

The online version of this article (doi:10.1186/s13104-015-1464-6) contains supplementary material, which is available to authorized users.

Entrance and survival of Brucella pinnipedialis hooded seal strain in human macrophages and epithelial cells

Marine mammal Brucella spp. have been isolated from pinnipeds (B. pinnipedialis) and cetaceans (B. ceti) from around the world. Although the zoonotic potential of marine mammal brucellae is largely unknown, reports of human disease exist. There are few studies of the mechanisms of bacterial intracellular invasion and multiplication involving the marine mammal Brucella spp. We examined the infective capacity of two genetically different B. pinnipedialis strains (reference strain; NTCT 12890 and a hooded seal isolate; B17) by measuring the ability of the bacteria to enter and replicate in cultured phagocytes and epithelial cells. Human macrophage-like cells (THP-1), two murine macrophage cell lines (RAW264.7 and J774A.1), and a human malignant epithelial cell line (HeLa S3) were challenged with bacteria in a gentamicin protection assay. Our results show that B. pinnipedialis is internalized, but is then gradually eliminated during the next 72 – 96 hours. Confocal microscopy revealed that intracellular B. pinnipedialis hooded seal strain colocalized with lysosomal compartments at 1.5 and 24 hours after infection. Intracellular presence of B. pinnipedialis hooded seal strain was verified by transmission electron microscopy. By using a cholesterol-scavenging lipid inhibitor, entrance of B. pinnipedialis hooded seal strain in human macrophages was significantly reduced by 65.8 % (± 17.3), suggesting involvement of lipid-rafts in intracellular entry. Murine macrophages invaded by B. pinnipedialis do not release nitric oxide (NO) and intracellular bacterial presence does not induce cell death. In summary, B. pinnipedialis hooded seal strain can enter human and murine macrophages, as well as human epithelial cells. Intracellular entry of B. pinnipedialis hooded seal strain involves, but seems not to be limited to, lipid-rafts in human macrophages. Brucella pinnipedialis does not multiply or survive for prolonged periods intracellulary.

POGO-DB--a database of pairwise-comparisons of genomes and conserved orthologous genes

POGO-DB (http://pogo.ece.drexel.edu/) provides an easy platform for comparative microbial genomics. POGO-DB allows users to compare genomes using pre-computed metrics that were derived from extensive computationally intensive BLAST comparisons of >2000 microbes. These metrics include (i) average protein sequence identity across all orthologs shared by two genomes, (ii) genomic fluidity (a measure of gene content dissimilarity), (iii) number of ‘orthologs’ shared between two genomes, (iv) pairwise identity of the 16S ribosomal RNA genes and (v) pairwise identity of an additional 73 marker genes present in >90% prokaryotes. Users can visualize these metrics against each other in a 2D plot for exploratory analysis of genome similarity and of how different aspects of genome similarity relate to each other. The results of these comparisons are fully downloadable. In addition, users can download raw BLAST results for all or user-selected comparisons. Therefore, we provide users with full flexibility to carry out their own downstream analyses, by creating easy access to data that would normally require heavy computational resources to generate. POGO-DB should prove highly useful for researchers interested in comparative microbiology and benefit the microbiome/metagenomic communities by providing the information needed to select suitable phylogenetic marker genes within particular lineages.

Entry and elimination of marine mammal Brucella spp. by hooded seal (Cystophora cristata) alveolar macrophages in vitro

A high prevalence of Brucellapinnipedialis serology and bacteriology positive animals has been found in the Northeast Atlantic stock of hooded seal (Cystophoracristata); however no associated gross pathological changes have been identified. Marine mammal brucellae have previously displayed different infection patterns in human and murine macrophages. To investigate if marine mammal Brucella spp. are able to invade and multiply in cells originating from a presumed host species, we infected alveolar macrophages from hooded seal with a B. pinnipedialis hooded seal isolate. Hooded seal alveolar macrophages were also challenged with B. pinnipedialis reference strain (NCTC 12890) from harbor seal (Phocavitulina), B. ceti reference strain (NCTC 12891) from harbor porpoise (Phocoenaphocoena) and a B. ceti Atlantic white-sided dolphin (Lagenorhynchusacutus) isolate (M83/07/1), to evaluate possible species-specific differences. Brucella suis 1330 was included as a positive control. Alveolar macrophages were obtained by post mortem bronchoalveolar lavage of euthanized hooded seals. Phenotyping of cells in the lavage fluid was executed by flow cytometry using the surface markers CD14 and CD18. Cultured lavage cells were identified as alveolar macrophages based on morphology, expression of surface markers and phagocytic ability. Alveolar macrophages were challenged with Brucella spp. in a gentamicin protection assay. Following infection, cell lysates from different time points were plated and evaluated quantitatively for colony forming units. Intracellular presence of B. pinnipedialis hooded seal isolate was verified by immunocytochemistry. Our results show that the marine mammal brucellae were able to enter hooded seal alveolar macrophages; however, they did not multiply intracellularly and were eliminated within 48 hours, to the contrary of B. suis that showed the classical pattern of a pathogenic strain. In conclusion, none of the four marine mammal strains tested were able to establish a persistent infection in primary alveolar macrophages from hooded seal.

Yeast proteome variations reveal different adaptive responses to grape must fermentation

Saccharomyces cerevisiae and S. uvarum are two domesticated species of the Saccharomyces sensu stricto clade that diverged around 100 Ma after whole-genome duplication. Both have retained many duplicated genes associated with glucose fermentation and are characterized by the ability to achieve grape must fermentation. Nevertheless, these two species differ for many other traits, indicating that they underwent different evolutionary histories. To determine how the evolutionary histories of S. cerevisiae and S. uvarum are mirrored on the proteome, we analyzed the genetic variability of the proteomes of domesticated strains of these two species by quantitative mass spectrometry. Overall, 445 proteins were quantified. Massive variations of protein abundances were found, that clearly differentiated the two species. Abundance variations in specific metabolic pathways could be related to phenotypic traits known to discriminate the two species. In addition, proteins encoded by duplicated genes were shown to be differently recruited in each species. Comparing the strain differentiation based on the proteome variability to those based on the phenotypic and genetic variations further revealed that the strains of S. uvarum and some strains of S. cerevisiae displayed similar fermentative performances despite strong proteomic and genomic differences. Altogether, these results indicate that the ability of S. cerevisae and S. uvarum to complete grape must fermentation arose through different evolutionary roads, involving different metabolic pathways and duplicated genes.

Metagenomic microbial community profiling using unique clade-specific marker genes

Metagenomic shotgun sequencing data can identify microbes populating a microbial community and their proportions, but existing taxonomic profiling methods are inefficient for increasingly large data sets. We present an approach that uses clade-specific marker genes to unambiguously assign reads to microbial clades more accurately and >50× faster than current approaches. We validated our metagenomic phylogenetic analysis tool, MetaPhlAn, on terabases of short reads and provide the largest metagenomic profiling to date of the human gut. It can be accessed at http://huttenhower.sph.harvard.edu/metaphlan/.

Brucella genomics as we enter the multi-genome era

The genus Brucella includes species considered among the worlds most important zoonotic pathogens, with brucellosis remaining a significant problem in large parts of the world. Over the last decade a number of Brucella genomes have been fully sequenced providing new insights into this relatively poorly understood group of organisms. In the forthcoming months and years, the availability of many additional genomes should help in further understanding of the evolution, host specificity and pathogenicity of this group as well as providing a resource to further improve epidemiological typing tools. This article describes progress to date and looks forward to the opportunities that should be afforded as we enter an era of multiple, freely available, Brucella genome sequences.

A review of Brucella infection in marine mammals, with special emphasis on Brucella pinnipedialis in the hooded seal (Cystophora cristata)

Brucella spp. were isolated from marine mammals for the first time in 1994. Two novel species were later included in the genus; Brucella ceti and Brucella pinnipedialis, with cetaceans and seals as their preferred hosts, respectively. Brucella spp. have since been isolated from a variety of marine mammals. Pathological changes, including lesions of the reproductive organs and associated abortions, have only been registered in cetaceans. The zoonotic potential differs among the marine mammal Brucella strains. Many techniques, both classical typing and molecular microbiology, have been utilised for characterisation of the marine mammal Brucella spp. and the change from the band-based approaches to the sequence-based approaches has greatly increased our knowledge about these strains. Several clusters have been identified within the B. ceti and B. pinnipedialis species, and multiple studies have shown that the hooded seal isolates differ from other pinniped isolates. We describe how different molecular methods have contributed to species identification and differentiation of B. ceti and B. pinnipedialis, with special emphasis on the hooded seal isolates. We further discuss the potential role of B. pinnipedialis for the declining Northwest Atlantic hooded seal population.

The genome sequence of Brucella pinnipedialis B2/94 sheds light on the evolutionary history of the genus Brucella

Background

Since the discovery of the Malta fever agent, Brucella melitensis, in the 19th century, six terrestrial mammal-associated Brucella species were recognized over the next century. More recently the number of novel Brucella species has increased and among them, isolation of species B. pinnipedialis and B. ceti from marine mammals raised many questions about their origin as well as on the evolutionary history of the whole genus.

Results

We report here on the first complete genome sequence of a Brucella strain isolated from marine mammals, Brucella pinnipedialis strain B2/94. A whole gene-based phylogenetic analysis shows that five main groups of host-associated Brucella species rapidly diverged from a likely free-living ancestor close to the recently isolated B. microti. However, this tree lacks the resolution required to resolve the order of divergence of those groups. Comparative analyses focusing on a) genome segments unshared between B. microti and B. pinnipedialis, b) gene deletion/fusion events and c) positions and numbers of Brucella specific IS711 elements in the available Brucella genomes provided enough information to propose a branching order for those five groups.

Conclusions

In this study, it appears that the closest relatives of marine mammal Brucella sp. are B. ovis and Brucella sp. NVSL 07-0026 isolated from a baboon, followed by B. melitensis and B. abortus strains, and finally the group consisting of B. suis strains, including B. canis and the group consisting of the single B. neotomae species. We were not able, however, to resolve the order of divergence of the two latter groups.

Identification of a Brucella spp. secreted effector specifically interacting with human small GTPase Rab2

Bacteria of the Brucella genus are facultative intracellular class III pathogens. These bacteria are able to control the intracellular trafficking of their vacuole, presumably by the use of yet unknown translocated effectors. To identify such effectors, we used a high-throughput yeast two-hybrid screen to identify interactions between putative human phagosomal proteins and predicted Brucella spp. proteins. We identified a specific interaction between the human small GTPase Rab2 and a Brucella spp. protein named RicA. This interaction was confirmed by GST-pull-down with the GDP-bound form of Rab2. A TEM-β-lactamase-RicA fusion was translocated from Brucella abortus to RAW264.7 macrophages during infection. This translocation was not detectable in a strain deleted for the virB operon, coding for the type IV secretion system. However, RicA secretion in a bacteriological culture was still observed in a ΔvirB mutant. In HeLa cells, a ΔricA mutant recruits less GTP-locked myc-Rab2 on its Brucella-containing vacuoles, compared with the wild-type strain. We observed altered kinetics of intracellular trafficking and faster proliferation of the B. abortusΔricA mutant in HeLa cells, compared with the wild-type control. Altogether, the data reported here suggest RicA as the first reported effector with a proposed function for B. abortus.