Genetic Characterization and Comparative Genome Analysis of Brucella melitensis Isolates from India

Phylogenetic Analysis and Comparative Genomics of Brucella abortus and Brucella melitensis Strains in Egypt

Brucellosis is a notifiable disease induced by a facultative intracellular Brucella pathogen. In this study, eight Brucella abortus and eighteen Brucella melitensis strains from Egypt were annotated and compared with RB51 and REV1 vaccines respectively. RAST toolkit in the BV-BRC server was used for annotation, revealing genome length of 3,250,377 bp and 3,285,803 bp, 3289 and 3323 CDS, 48 and 49 tRNA genes, the same number of rRNA (3) genes, 583 and 586 hypothetical proteins, 2697 and 2726 functional proteins for B. abortus and B. melitensis respectively. B. abortus strains exhibit a similar number of candidate genes, while B. melitensis strains showed some differences, especially in the SRR19520422 Faiyum strain. Also, B. melitensis clarified differences in antimicrobial resistance genes (KatG, FabL, MtrA, MtrB, OxyR, and VanO-type) in SRR19520319 Faiyum and (Erm C and Tet K) in SRR19520422 Faiyum strain. Additionally, the whole genome phylogeny analysis proved that all B. abortus strains were related to vaccinated animals and all B. melitensis strains of Menoufia clustered together and closely related to Gharbia, Dameitta, and Kafr Elshiek. The Bowtie2 tool identified 338 (eight B. abortus) and 4271 (eighteen B. melitensis) single nucleotide polymorphisms (SNPs) along the genomes. These variants had been annotated according to type and impact. Moreover, thirty candidate genes were predicted and submitted at GenBank (24 in B. abortus) and (6 in B. melitensis). This study contributes significant insights into genetic variation, virulence factors, and vaccine-related associations of Brucella pathogens, enhancing our knowledge of brucellosis epidemiology and evolution in Egypt.

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.

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.

A novel multiplex qPCR‑HRM assay for the simultaneous detection of four abortive zoonotic agents in cattle, sheep, and goats

Abortifacient pathogens induce substantial economic losses in the livestock industry worldwide, and many of these pathogens are zoonotic, impacting human health. As Brucella spp., Coxiella burnetii, Leptospira spp., and Listeria monocytogenes cause abortion, rapid differential molecular diagnostic tests are needed to facilitate early and accurate detection of abortion to establish effective control measures. However, the available molecular methods are laborious, time-consuming, or costly. Therefore, we developed and validated a novel multiplex real-time polymerase chain reaction (qPCR) method based on high-resolution melting (HRM) curve analysis to simultaneously detect and differentiate four zoonotic abortifacient agents in cattle, goats, and sheep. Our HRM assay generated four well-separated melting peaks allowing the differentiation between the four zoonotic abortifacients. Out of 216 DNA samples tested, Brucella spp. was detected in 45 samples, Coxiella burnetii in 57 samples, Leptospira spp. in 12 samples, and Listeria monocytogenes in 19 samples, co-infection with Brucella spp. and Coxiella burnetii in 41 samples, and 42 samples were negative. This assay demonstrated good analytical sensitivity, specificity, and reproducibility. This is a valuable rapid, cost-saving, and reliable diagnostic tool for detecting individual and co-infections for zoonotic abortifacient agents in ruminants.

Genome Sequencing and Comparative Genomics of Indian Isolates of Brucella melitensis

Brucella melitensis causes small ruminant brucellosis and a zoonotic pathogen prevalent worldwide. Whole genome phylogeny of all available B. melitensis genomes (n = 355) revealed that all Indian isolates (n = 16) clustered in the East Mediterranean lineage except the ADMAS-GI strain. Pangenome analysis indicated the presence of limited accessory genomes with few clades showing specific gene presence/absence pattern. A total of 43 virulence genes were predicted in all the Indian strains of B. melitensis except 2007BM-1 (ricA and wbkA are absent). Multilocus sequence typing (MLST) analysis indicated all except one Indian strain (ADMAS-GI) falling into sequence type (ST 8). In comparison with MLST, core genome phylogeny indicated two major clusters (>70% bootstrap support values) among Indian strains. Clusters with <70% bootstrap support values represent strains with diverse evolutionary origins present among animal and human hosts. Genetic relatedness among animal (sheep and goats) and human strains with 100% bootstrap values shows its zoonotic transfer potentiality. SNP-based analysis indicated similar clustering to that of core genome phylogeny. Among the Indian strains, the highest number of unique SNPs (112 SNPs) were shared by a node that involved three strains from Tamil Nadu. The node SNPs involved several peptidase genes like U32, M16 inactive domain protein, clp protease family protein, and M23 family protein and mostly represented non-synonymous (NS) substitutions. Vaccination has been followed in several parts of the world to prevent small ruminant brucellosis but not in India. Comparison of Indian strains with vaccine strains showed that M5 is genetically closer to most of the Indian strains than Rev.1 strain. The presence of most of the virulence genes among all Indian strains and conserved core genome compositions suggest the use of any circulating strain/genotypes for the development of a vaccine candidate for small ruminant brucellosis in India.

Risk factors for Brucella spp. and Coxiella burnetii infection among small ruminants in Eastern India

Small ruminants are the main reservoirs for brucellosis and coxiellosis, two zoonotic diseases affecting livestock production, and posing a public health threat in India. Understanding disease prevalence and risk factors associated with small ruminant infection can help mitigate disease transmission.

We report a cross-sectional survey in the states of Assam and Odisha in Eastern India. We interviewed 244 farmers to assess knowledge, attitude and practices relevant to brucellosis and coxiellosis infection. Serum samples from 411 goats and 21 sheep were analysed using enzyme-linked immunosorbent assay and Rose-Bengal Brucella agglutination plate test. Higher Brucella and Coxiella burnetii seroprevalence were found in Odisha (22% and 11.5%, respectively) than Assam (9.8% and 1.6%, respectively), and certain districts in Odisha were at higher risk. No association was found between seropositive animals and clinical signs, a challenge when attempting to identify seropositive animals in the herd. None of the farmers interviewed were aware of brucellosis, its aetiology, clinical form, or zoonotic risk. This study acts as a first indication of the extent of these diseases among small ruminants in these Indian states, highlighting how farming practices are associated with increased risk of infection. More research is urgently needed to mitigate zoonoses transmission in this region.

Genome-scale approach to study the genetic relatedness among Brucella melitensis strains

Brucellosis is an important zoonotic disease that affects both humans and animals. To date, laboratory surveillance is still essentially based on the traditional MLVA-16 methodology and the associated epidemiological information is frequently scarce. Our goal was to contribute to the improvement of Brucella spp. surveillance through the implementation of a whole genome sequencing (WGS) approach. We created a curated ready-to-use species-specific wgMLST scheme enrolling a panel of 2656 targets (http://doi.org/10.5281/zenodo.3575026) and used this schema to perform a retrospective analysis of the genetic relatedness among B. melitensis strains causing human infection in Portugal (a country where brucellosis is an endemic disease) from 2010 to 2018. The strains showed a phylogenetic clustering within genotype II (25 out of 36) and IV (4 out of 36), and shared clades with strains isolated from countries with which Portugal has intense food trading, tourism and similar eating habits, such as Spain, Italy and Greece. In addition, our results point to the identification of strong associations between B. melitensis strains, likely underlying missed "outbreaks" as 22 out of the 36 strains showed genetic linkage with others. In fact, the applied gene-by-gene approach grouped these strains into six genetic clusters each one containing putative epidemiological links. Nevertheless, more studies will be needed in order to define the appropriate range of cut-offs (probable non-static cut-offs) that best illustrate the association between genetic linkage and epidemiological information and may serve as alerts for the health authorities. The release of this freely available and scalable schema contributes to the required technological transition for laboratorial surveillance of brucellosis and will facilitate the assessment of ongoing and future outbreaks in order to prevent the transmission spread.

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.

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.

Whole genome sequencing of Brucella melitensis isolated from 57 patients in Germany reveals high diversity in strains from Middle East

Brucellosis, a worldwide common bacterial zoonotic disease, has become quite rare in Northern and Western Europe. However, since 2014 a significant increase of imported infections caused by Brucella (B.) melitensis has been noticed in Germany. Patients predominantly originated from Middle East including Turkey and Syria. These circumstances afforded an opportunity to gain insights into the population structure of Brucella strains. Brucella-isolates from 57 patients were recovered between January 2014 and June 2016 with culture confirmed brucellosis by the National Consultant Laboratory for Brucella. Their whole genome sequences were generated using the Illumina MiSeq platform. A whole genome-based SNP typing assay was developed in order to resolve geographically attributed genetic clusters. Results were compared to MLVA typing results, the current gold-standard of Brucella typing. In addition, sequences were examined for possible genetic variation within target regions of molecular diagnostic assays. Phylogenetic analyses revealed spatial clustering and distinguished strains from different patients in either case, whereas multiple isolates from a single patient or technical replicates showed identical SNP and MLVA profiles. By including WGS data from the NCBI database, five major genotypes were identified. Notably, strains originating from Turkey showed a high diversity and grouped into seven subclusters of genotype II. MLVA analysis congruently clustered all isolates and predominantly matched the East Mediterranean genetic clade. This study confirms whole-genome based SNP-analysis as a powerful tool for accurate typing of B. melitensis. Furthermore it allows special allocation and therefore provides useful information on the geographic origin for trace-back analysis. However, the lack of reliable metadata in public databases often prevents a resolution below geographic regions or country levels and corresponding precise trace-back analysis. Once this obstacle is resolved, WGS-derived bacterial typing adds an important method to complement epidemiological surveys during outbreak investigations. This is the first report of a detailed genetic investigation of an extensive collection of B. melitensis strains isolated from human cases in Germany.

Draft Genome Sequences of Five Clinical Strains of Brucella melitensis Isolated from Patients Residing in Kuwait

Human brucellosis is a neglected and underrecognized infection of widespread geographic distribution. Brucellosis is present on all inhabited continents and endemic in many areas of the world, including Kuwait and the Middle East. Here, we present draft genome assemblies of five Brucella melitensis strains isolated from brucellosis patients in Kuwait.