Genome-wide unique insertion sequences among five Brucella species and demonstration of differential identification of Brucella by multiplex PCR assay

Brucella melitensis: Divergence Among Indian Strains and Genetic Characterization of a Strain Isolated from Cattle

Brucella melitensis primarily affects sheep, goats and is associated with brucellosis in humans, which is one of the world's most widespread neglected zoonotic disease. The current study attempted the determination of genetic diversity through comparative genome analysis of B. melitensis strains reported from India with other countries. The study also reports the isolation and identification of B. melitensis BMNDDB8664 from a cow with a history of abortion, whole-genome sequencing (WGS), determination of virulence factors, genotyping, and comparative genome analysis. Multilocus sequence typing, Multiple locus variable number of tandem repeats analysis (MLVA), and WGS based phylogeny revealed the predominance of ST-8 and genotypes (116 and II respectively) that clustered to the East Mediterranean lineage. Identification of hitherto unreported genotypes by MLVA also indicated the existence and circulation of West Mediterranean and American lineages in India. Though the AMOS-PCR results suggest the BMNDDB8664 isolate as Brucella abortus, the outcomes from multiplex PCR, ribosomal multilocus sequence typing, and WGS analysis confirmed it as B. melitensis. The analysis revealed the presence of adeF gene (aids conferring resistance to fluoro-quinolone and tetracyclines). The isolate lacked two important T4SS genes virB2 and virB7 genes (roles in infection and rifampicin resistance respectively) and also lacked the Brucella suis mprF gene that aids intracellular survival. Further, BMNDDB8664 lacked some of the genes associated with LPS synthesis (wbkB, wbkC) and transport (wzm, wzt) and hence, is most likely a rough strain. WGS-based phylogenetic analysis revealed close genetic relatedness of this BMNDDB8664 with a sheep isolate and two human isolates. The results prompt systematic, broad-based epidemiological studies on brucella infection at the species level. For effective control of human brucellosis, a concerted One Health approach with studies encircling the identification of aetiology at species, strain level to find their prevalence, spread, and inter-host transmission patterns need to be understood, for better design and implementation of effective control strategies in India and other endemic regions.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12088-023-01081-w.

Design of multi-epitope vaccine candidate against Brucella type IV secretion system (T4SS)

Brucellosis is a common zoonosis, which is caused by Brucella infection, and Brucella often infects livestock, leading to abortion and infertility. At present, human brucellosis remains one of the major public health problems in China. According to previous research, most areas in northwest China, including Xinjiang, Tibet, and other regions, are severely affected by Brucella. Although there are vaccines against animal Brucellosis, the effect is often poor. In addition, there is no corresponding vaccine for human Brucellosis infection. Therefore, a new strategy for early prevention and treatment of Brucella is needed. A multi-epitope vaccine should be developed. In this study, we identified the antigenic epitopes of the Brucella type IV secretion system VirB8 and Virb10 using an immunoinformatics approach, and screened out 2 cytotoxic T lymphocyte (CTL) epitopes, 9 helper T lymphocyte (HTL) epitopes, 6 linear B cell epitopes, and 6 conformational B cell epitopes. These advantageous epitopes are spliced together through different linkers to construct a multi-epitope vaccine. The silico tests showed that the multi-epitope vaccine was non-allergenic and had a strong interaction with TLR4 molecular docking. In immune simulation results, the vaccine construct may be useful in helping brucellosis patients to initiate cellular and humoral immunity. Overall, our findings indicated that the multi-epitope vaccine construct has a high-quality structure and suitable characteristics, which may provide a theoretical basis for the development of a Brucella vaccine.

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.

Rapid, ultrasensitive, and highly specific identification of Brucella abortus utilizing multiple cross displacement amplification combined with a gold nanoparticles-based lateral flow biosensor

Brucella abortus (B. abortus) as an important infectious agent of bovine brucellosis cannot be ignored, especially in countries/regions dominated by animal husbandry. Thus, the development of an ultrasensitive and highly specific identification technique is an ideal strategy to control the transmission of bovine brucellosis. In this report, a novel detection protocol, which utilizes multiple cross displacement amplification (MCDA) combined with a gold nanoparticles-based lateral flow biosensor (AuNPs-LFB) targeting the BruAb2_0168 gene was successfully devised and established for the identification of B. abortus (termed B. abortus-MCDA-LFB). Ten specific primers containing engineered C1-FAM (carboxyfluorescein) and D1-biotin primers were designed according to the MCDA reaction mechanism. These genomic DNA extracted from various bacterial strains and whole blood samples were used to optimize and evaluate the B. abortus-MCDA-LFB assay. As a result, the optimal reaction conditions for the B. abortus-MCDA-LFB assay were 66°C for 40 min. The limit of detection of the B. abortus-MCDA-LFB was 10 fg/μl (~3 copies/μl) for genomic DNA extracted from pure cultures of B. abortus isolate. Meanwhile, the B. abortus-MCDA-LFB assay accurately identified all tested B. abortus strains, and there was no cross-reaction with non-B. abortus pathogens. Moreover, the detection workflow of the B. abortus-MCDA-LFB assay for whole blood samples can be completed within 70 min, and the cost of a single test is approximately 5.0 USD. Taken together, the B. abortus-MCDA-LFB assay is a visual, fast, ultrasensitive, low-cost, easy-to-operate, and highly specific detection method, which can be used as a rapid identification tool for B. abortus infections.

Design of a multi-epitope vaccine candidate against Brucella melitensis

Brucella is a typical facultative intracellular bacterium that can cause zoonotic infections. For Brucella, it is difficult to eliminate with current medical treatment. Therefore, a multi-epitope vaccine (MEV) should be designed to prevent Brucella infection. For this purpose, we applied the reverse vaccinology approach from Omp10, Omp25, Omp31 and BtpB. Finally, we obtained 13 cytotoxic T lymphocyte (CTL) epitopes, 17 helper T lymphocyte (HTL) epitopes, 9 linear B cell epitopes, and 2 conformational B cell epitopes for further study. To keep the protein folded normally, we linked AAY, GPGPG, and KK to CTL epitopes, HTL epitopes, and B cell epitopes, respectively. The N-terminal of the vaccine peptide is supplemented with appropriate adjuvants to enhance immunogenicity. To evaluate its immunogenicity, stability, safety, and feasibility, a final MEV containing 806 amino acids was constructed by linking linkers and adjuvants. In addition, molecular docking and molecular dynamics simulations were performed to verify the affinity and stability of the MEV-TLR4. Then, codon adaptation and in silico cloning studies were carried out to identify the possible codons for expressing the MEV. In animal experiments, the results demonstrated that the MEV had high immunogenicity. Collectively, this study provided a theoretical basis for the development of a Brucella vaccine.

Genetic and molecular Omp25 analyses from worldwide Brucella canis strains: Possible mutational influences in protein function

Brucella canis is responsible for canine brucellosis, a neglected zoonotic disease. The omp25 gene has been described as an important marker for Brucella intra-species differentiation, in addition to the ability to interact with the host immune system. Therefore, this study investigated the omp25 sequence from B. canis strains associated to a phylogenetic characterization and the unveiling of the molecular structure. In vitro analyses comprised DNA extraction, PCR, and sequencing of omp25 from 19 B. canis strains. Moreover, in silico analyses were performed at nucleotide level for phylogenetic characterization and evolutionary history of B. canis omp25 gene; and in amino acid level including modeling, dynamics, and epitope prediction of B. canis Omp25 protein. Here, we identified a new mutation, L109P, which diverges the worldwide omp25 sequences in two large branches. Interestingly, this mutation appears to have epidemiology importance, based on a geographical distribution of B. canis strains. Structural and molecular dynamics analyses of Omp25 revealed that Omp25L109P does not sustain its native β-barrel. Likewise, the conformation of B-cell epitope on the mutated region was changed in Omp25L109P protein. Even without an evolutive marker, the new identified mutation appears to affect the basic function of B. canis Omp25 protein, which could indicate virulence adaptation for some B. canis strains in a context of geographical disposition.