Comparative genomic analysis of Brucella abortus vaccine strain 104M reveals a set of candidate genes associated with its virulence attenuation

Rapid Identification of Brucella Genus and Species In Silico and On-Site Using Novel Probes with CRISPR/Cas12a

Human brucellosis caused by Brucella is a widespread zoonosis that is prevalent in many countries globally. The high homology between members of the Brucella genus and Ochrobactrum spp. often complicates the determination of disease etiology in patients. The efficient and reliable identification and distinction of Brucella are of primary interest for both medical surveillance and outbreak purposes. A large amount of genomic data for the Brucella genus was analyzed to uncover novel probes containing single-nucleotide polymorphisms (SNPs). GAMOSCE v1.0 software was developed based on the above novel eProbes. In conjunction with clinical requirements, an RPA-Cas12a detection method was developed for the on-site determination of B. abortus and B. melitensis by fluorescence and lateral flow dipsticks (LFDs). We demonstrated the potential of these probes for rapid and accurate detection of the Brucella genus and five significant Brucella species in silico using GAMOSCE. GAMOSCE was validated on different Brucella datasets and correctly identified all Brucella strains, demonstrating a strong discrimination ability. The RPA-Cas12a detection method showed good performance in detection in clinical blood samples and veterinary isolates. We provide both in silico and on-site methods that are convenient and reliable for use in local hospitals and public health programs for the detection of brucellosis.

A novel Galleria mellonella experimental model for zoonotic pathogen Brucella

Brucellosis is a major threat to public health and animal husbandry. Several in vivo vertebrate models, such as mice, guinea pigs, and nonhuman primates, have been used to study Brucella pathogenesis, bacteria-host interactions, and vaccine efficacy. However, these models have limitations whereas the invertebrate Galleria mellonella model is a cost-effective and ethical alternative. The aim of the present study was to examine the invertebrate G. mellonella as an in vivo infection model for Brucella. Infection assays were employed to validate the fitness of the larval model for Brucella infection and virulence evaluation. The protective efficacy of immune sera was evaluated by pre-incubated with a lethal dose of bacteria before infection. The consistency between the mouse model and the larval model was confirmed by assessing the protective efficacy of two Brucella vaccine strains. The results show that G. mellonella could be infected by Brucella strains, in a dose- and temperature-dependent way. Moreover, this larval model can effectively evaluate the virulence of Brucella strains in a manner consistent with that of mammalian infection models. Importantly, this model can assess the protective efficacy of vaccine immune sera within a day. Further investigation implied that haemolymph played a crucial role in the protective efficacy of immune sera. In conclusion, G. mellonella could serve as a quick, efficient, and reliable model for evaluating the virulence of Brucella strains and efficacy of immune sera in an ethical manner.

Russian collection of Brucella abortus vaccine strains: annotation, implementation and genomic analysis

Over the past 10 years, immunization of cattle in Russia has been performed using vaccines from Brucella abortus strains 82, 19 and 75/79. To prevent brucellosis in small ruminants, two vaccines have been used, from the Brucella melitensis strain REV-1 and the B. abortus strain 19; note that twice as many animals have been immunized with the former vaccine than with the latter vaccine. The disadvantage of using these preparations is the formation of prolonged post-vaccination seropositivity, which is especially pronounced in animals after immunization with vaccines from B. abortus strain 19 and B. melitensis strain REV-1. This study aims to perform the whole genome sequencing of Brucella vaccine strains from the Russian collection. A bioinformatics analysis of the genomic data proved that the vaccine strains 75/79AB, 82, R-1096, and the KV 17/100 belong to ST-2, 104 M to ST-1, KV 13/100 to ST-5. This analysis allowed us to characterize vaccine strains's phylogenetic relationships and to prove the close relation of vaccine strains 75/79AB, 82, R-1096. Also, we defined candidate mutations in genes pmm, wbdA, wbkA, wboA, and eryB, which could be responsible for the attenuated virulence of vaccine strains. The complete genomic sequences of B. abortus strains make further studies of bacterial pathogenicity determinants and virulence phenotype feasible, as well as their use in quality control of animal medicines.

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.

Evaluation and Differential Diagnosis of a Genetic Marked Brucella Vaccine A19ΔvirB12 for Cattle

Brucella abortus is an important zoonotic pathogen that causes severe economic loss to husbandry and poses a threat to human health. The B. abortus A19 live vaccine has been extensively used to prevent bovine brucellosis in China. However, it is difficult to distinguish the serological response induced by A19 from that induced by natural infection. In this study, a novel genetically marked vaccine, A19ΔvirB12, was generated and evaluated. The results indicated that A19ΔvirB12 was able to provide effective protection against B. abortus 2308 (S2308) challenge in mice. Furthermore, the safety and protective efficacy of A19ΔvirB12 have been confirmed in natural host cattle. Additionally, the VirB12 protein allowed for serological differentiation between the S2308 challenge/natural infection and A19ΔvirB12 vaccination. However, previous studies have found that the accuracy of the serological detection based on VirB12 needs to be improved. Therefore, we attempted to identify potential supplementary antigens with differential diagnostic functions by combining label-free quantitative proteomics and protein chip technology. Twenty-six proteins identified only in S2308 were screened; among them, five proteins were considered as potential supplementary antigens. Thus, the accuracy of the differential diagnosis between A19ΔvirB12 immunization and field infection may be improved through multi-antigen detection. In addition, we explored the possible attenuation factors of Brucella vaccine strain. Nine virulence factors were downregulated in A19ΔvirB12. The downregulation pathways of A19ΔvirB12 were significantly enriched in quorum sensing, ATP-binding cassette transporter, and metabolism. Several proteins related to cell division were significantly downregulated, while some proteins involved in transcription were upregulated in S2308. In conclusion, our results contribute to the control and eradication of brucellosis and provide insights into the mechanisms underlying the attenuation of A19ΔvirB12.

Screening of potential vaccine candidates against pathogenic Brucella spp. using compositive reverse vaccinology

Brucella spp. are Gram-negative, facultative intracellular bacteria that cause brucellosis in humans and various animals. The threat of brucellosis has increased, yet currently available live attenuated vaccines still have drawbacks. Therefore, subunit vaccines, produced using protein antigens and having the advantage of being safe, cost-effective and efficacious, are urgently needed. In this study, we used core proteome analysis and a compositive RV methodology to screen potential broad-spectrum antigens against 213 pathogenic strains of Brucella spp. with worldwide geographic distribution. Candidate proteins were scored according to six biological features: subcellular localization, antigen similarity, antigenicity, mature epitope density, virulence, and adhesion probability. In the RV analysis, a total 32 candidate antigens were picked out. Of these, three proteins were selected for assessment of immunogenicity and preliminary protection in a mouse model: outer membrane protein Omp19 (used as a positive control), type IV secretion system (T4SS) protein VirB8, and type I secretion system (T1SS) protein HlyD. These three antigens with a high degree of conservation could induce specific humoral and cellular immune responses. Omp19, VirB8 and HlyD could substantially reduce the organ bacterial load of B. abortus S19 in mice and provide varying degrees of protection. In this study, we demonstrated the effectiveness of this unique strategy for the screening of potential broad-spectrum antigens against Brucella. Further evaluation is needed to identify the levels of protection conferred by the vaccine antigens against wild-type pathogenic Brucella species challenge.

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.

Different Impacts of MucR Binding to the babR and virB Promoters on Gene Expression in Brucella abortus 2308

The protein MucR from Brucella abortus has been described as a transcriptional regulator of many virulence genes. It is a member of the Ros/MucR family comprising proteins that control the expression of genes important for the successful interaction of α-proteobacteria with their eukaryotic hosts. Despite clear evidence of the role of MucR in repressing virulence genes, no study has been carried out so far demonstrating the direct interaction of this protein with the promoter of its target gene babR encoding a LuxR-like regulator repressing virB genes. In this study, we show for the first time the ability of MucR to bind the promoter of babR in electrophoretic mobility shift assays demonstrating a direct role of MucR in repressing this gene. Furthermore, we demonstrate that MucR can bind the virB gene promoter. Analyses by RT-qPCR showed no significant differences in the expression level of virB genes in Brucella abortus CC092 lacking MucR compared to the wild-type Brucella abortus strain, indicating that MucR binding to the virB promoter has little impact on virB gene expression in B. abortus 2308. The MucR modality to bind the two promoters analyzed supports our previous hypothesis that this is a histone-like protein never found before in Brucella.

First isolation and characterization of Brucella suis from yak

Brucella spp., facultative intracellular pathogens that can persistently colonize animal host cells and cause zoonosis, affect public health and safety. A Brucella strain was isolated from yak in Qinghai Province. To detect whether this isolate could cause an outbreak of brucellosis and to reveal its genetic characteristics, several typing and whole-genome sequencing methods were applied to identify its species and genetic characteristics. Phylogenetic analysis based on MLVA and whole-genome sequencing revealed the genetic characteristics of the isolated strain. The results showed that the isolated strain is a B. suis biovar 1 smooth strain, and this isolate was named B. suis QH05. The results of comparative genomics and MLVA showed that B. suis QH05 is not a vaccine strain. Comparison with other B. suis strains isolated from humans and animals indicated that B. suis QH05 may be linked to specific animal and human sources. In conclusion, B. suis QH05 does not belong to the Brucella epidemic species in China, and as the first isolation of B. suis from yak, this strain expands the host range of B. suis.

A rapid minor groove binder PCR method for distinguishing the vaccine strain Brucella abortus 104M

Background

Brucellosis is a widespread zoonotic disease caused by Gram-negative Brucella bacteria. Immunisation with attenuated vaccine is an effective method of prevention, but it can interfere with diagnosis. Live, attenuated Brucella abortus strain 104M has been used for the prevention of human brucellosis in China since 1965. However, at present, no fast and reliable method exists that can distinguish this strain from field strains. Single nucleotide polymorphism (SNP)-based assays offer a new approach for such discrimination. SNP-based minor groove binder (MGB) and Cycleave assays have been used for rapid identification of four Brucella vaccine strains (B. abortus strains S19, A19 and RB51, and B. melitensis Rev1). The main objective of this study was to develop a PCR assay for rapid and specific detection of strain 104M.

Results

We developed a SNP-based MGB PCR assay that could successfully distinguish strain 104M from 18 representative strains of Brucella (B. abortus biovars 1, 2, 3, 4, 5, 6, 7 and 9, B. melitensis biovars 1, 2 and 3, B. suis biovars 1, 2, 3 and 4, B. canis, B. neotomae, and B. ovis), four Brucella vaccine strains (A19, S19, S2, M5), and 55 Brucella clinical field strains. The assay gave a negative reaction with four non-Brucella species (Escherichia coli, Pasteurella multocida, Streptococcus suis and Pseudomonas aeruginosa). The minimum sensitivity of the assay, evaluated using 10-fold dilutions of chromosomal DNA, was 220 fg for the 104M strain and 76 fg for the single non-104M Brucella strain tested (B. abortus A19). The assay was also reproducible (intra- and inter-assay coefficients of variation = 0.006–0.022 and 0.012–0.044, respectively).

Conclusions

A SNP-based MGB PCR assay was developed that could straightforwardly and unambiguously distinguish B. abortus vaccine strain 104M from non-104M Brucella strains. Compared to the classical isolation and identification approaches of bacteriology, this real-time PCR assay has substantial advantages in terms of simplicity and speed, and also reduces potential exposure to live Brucella. The assay developed is therefore a simple, rapid, sensitive, and specific tool for brucellosis diagnosis and control.

Genomic analysis of the original Elberg Brucella melitensis Rev.1 vaccine strain reveals insights into virulence attenuation

The live attenuated Brucella melitensis Rev.1 Elberg-originated vaccine strain has been widely used to control brucellosis in small ruminants. However, despite extensive research, the molecular mechanisms underlying the attenuation of this strain are still unknown. In the current study, we conducted a comprehensive comparative analysis of the whole-genome sequence of Rev.1 against that of the virulent reference strain, B. melitensis 16M. This analysis revealed five regions of insertion and three regions of deletion within the Rev.1 genome, among which, one large region of insertion, comprising 3,951 bp, was detected in the Rev.1 genome. In addition, we found several missense mutations within important virulence-related genes, which may be used to determine the mechanism underlying virulence attenuation. Collectively, our findings provide new insights into the Brucella virulence mechanisms and, therefore, may serve as a basis for the rational design of new Brucella vaccines.

Relative Quantitative Proteomic Analysis of Brucella abortus Reveals Metabolic Adaptation to Multiple Environmental Stresses

Brucella spp. are facultative intracellular pathogens that cause chronic brucellosis in humans and animals. The virulence of Brucella primarily depends on its successful survival and replication in host cells. During invasion of the host tissue, Brucella is simultaneously subjected to a variety of harsh conditions, including nutrient limitation, low pH, antimicrobial defenses, and extreme levels of reactive oxygen species (ROS) via the host immune response. This suggests that Brucella may be able to regulate its metabolic adaptation in response to the distinct stresses encountered during its intracellular infection of the host. An investigation into the differential proteome expression patterns of Brucella grown under the relevant stress conditions may contribute toward a better understanding of its pathogenesis and adaptive response. Here, we utilized a mass spectrometry-based label-free relative quantitative proteomics approach to investigate and compare global proteomic changes in B. abortus in response to eight different stress treatments. The 3 h short-term in vitro single-stress and multi-stress conditions mimicked the in vivo conditions of B. abortus under intracellular infection, with survival rates ranging from 3.17 to 73.17%. The proteomic analysis identified and quantified a total of 2,272 proteins and 74% of the theoretical proteome, thereby providing wide coverage of the B. abortus proteome. By including eight distinct growth conditions and comparing these with a control condition, we identified a total of 1,221 differentially expressed proteins (DEPs) that were significantly changed under the stress treatments. Pathway analysis revealed that most of the proteins were involved in oxidative phosphorylation, ABC transporters, two-component systems, biosynthesis of secondary metabolites, the citrate cycle, thiamine metabolism, and nitrogen metabolism; constituting major response mechanisms toward the reconstruction of cellular homeostasis and metabolic balance under stress. In conclusion, our results provide a better understanding of the global metabolic adaptations of B. abortus associated with distinct environmental stresses. The identification of proteins necessary for stress resistance is crucial toward elucidating the infectious process in order to control brucellosis, and may facilitate the discovery of novel therapeutic targets and effective vaccines.

A comprehensive proteogenomic study of the human Brucella vaccine strain 104 M

BACKGROUND

Brucella spp. are Gram-negative, facultative intracellular pathogens that cause brucellosis in both humans and animals. The B. abortus vaccine strain 104 M is the only vaccine available in China for the prevention of brucellosis in humans. Although the B. abortus 104 M genome has been fully sequenced, the current genome annotations are not yet complete. In addition, the main mechanisms underpinning its residual toxicity and vaccine-induced immune protection have yet to be elucidated. Mapping the proteome of B. abortus 104 M will help to improve genome annotation quality, thereby facilitating a greater understanding of its biology.

RESULTS

In this study, we utilized a proteogenomic approach that combined subcellular fractionation and peptide fractionation to perform a whole-proteome analysis and genome reannotation of B. abortus 104 M using high-resolution mass spectrometry. In total, 1,729 proteins (56.3% of 3,072) including 218 hypothetical proteins were identified using the culture conditions that were employed this study. The annotations of the B. abortus 104 M genome were also refined following identification and validation by reverse transcription-PCR. In addition, 14 pivotal virulence factors and 17 known protective antigens known to be involved in residual toxicity and immune protection were confirmed at the protein level following induction by the 104 M vaccine. Moreover, a further insight into the cell biology of multichromosomal bacteria was obtained following the elucidation of differences in protein expression levels between the small and large chromosomes.

CONCLUSIONS

The work presented in this report used a proteogenomic approach to perform whole-proteome analysis and genome reannotation in B. abortus 104 M; this work helped to improve genome annotation quality. Our analysis of virulence factors, protective antigens and other protein effectors provided the basis for further research to elucidate the mechanisms of residual toxicity and immune protection induced by the 104 M vaccine. Finally, the potential link between replication dynamics, gene function, and protein expression levels in this multichromosomal bacterium was detailed.