Genotyping of human Brucella melitensis biovar 3 isolated from Shanxi Province in China by MLVA16 and HOOF

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.

Molecular typing methods to characterize Brucella spp. from animals: A review

Brucellosis is an infectious disease of animals that can infect humans. The disease causes significant economic losses and threatens human health. A timely and accurate disease diagnosis plays a vital role in the identification of brucellosis. In addition to traditional diagnostic methods, molecular methods allow diagnosis and typing of the causative agent of brucellosis. This review will discuss various methods, such as Bruce-ladder, Suiladder, high-resolution melt analysis, restriction fragment length polymorphism, multilocus sequence typing, multilocus variable-number tandem repeat analysis, and whole-genome sequencing single-nucleotide polymorphism, for the molecular typing of Brucella and discuss their advantages and disadvantages.

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.

An Epidemiological Study of Brucellosis in Different Animal Species from the Al-Qassim Region, Saudi Arabia

Brucellosis is a zoonotic bacterial illness that affects humans and a variety of domestic animals, especially ruminants. It is mostly transmitted through the consumption of contaminated drinks, foods, undercooked meat, or unpasteurized milk or contact with infected animals. Therefore, the present study aimed to investigate the seroprevalence of brucellosis in camels, sheep, and goat herds in the Qassim region, Saudi Arabia, using commonly used diagnostic serological procedures such as the Rose Bengal test (RBT), complement fixation test (CFT), and enzyme-linked immunosorbent assay (ELISA). The seroprevalence of brucellosis in camels, sheep, and goats was determined in the selected areas using a cross-sectional study design and a total of 690 farm animals of both sexes of different ages from the three animal species (274 camels, 227 sheep, and 189 goats). According to RBT results, 65 sera were positive for brucellosis, including 15 (5.47%) for camels, 32 (14.09%) for sheep, and 18 (9.50%) for goats. CFT and c-ELISA were performed as confirmatory tests on positive samples resulting from RBT. With c-ELISA, 60 serum samples were confirmed positive, in 14 (5.10%), 30 (13.21%), and 16 (8.46%) camels, sheep, and goats, respectively. There were 59 serum samples confirmed as positive for CFT, including 14 (5.11%), 29 (12.77%), and 16 (8.46%) for camels, sheep, and goats, respectively. Overall, the highest seroprevalence of brucellosis was found in sheep while the least was found in camels from the three tests (RBT, c-ELISA, and CFT). The highest seroprevalence of brucellosis was found in sheep while the least seroprevalence was found in camels. There was also a higher seroprevalence of brucellosis among female animals than males as well as among old animals than young animals. The study, thus, demonstrates brucellosis seroprevalence among farm animals (camels, sheep, and goats) and the significance of intervention measures against brucellosis incidence in both humans and animals through the creation of public awareness and other relevant policy measures such as livestock vaccination, effective hygiene management, and adequate quarantine or serological analysis for newly introduced animals.

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.

Molecular investigation of infection sources and transmission chains of brucellosis in Zhejiang, China

In the present study, a total of 7793 samples from 5 different types of hosts were collected and tested, with a seroprevalence of 2.4% (184/7793). Although the seroprevalence of human and animal brucellosis is relatively low, numbers of human brucellosis cases reported have increased continuously from 2004 to 2018. A total of 118 Brucella strains containing 4 biotypes were obtained, including Brucella melitensis bv.1 (n = 8) and bv.3 (n = 106), Brucella abortus bv.3 (n = 3) and bv.7 (n = 1). Twenty-one shared MLVA-16 genotypes, each composed of 2 to 19 strains obtained from different hosts, suggest the occurrence of a brucellosis outbreak epidemic with multiple source points and laboratory infection events. Moreover, 30 shared MLVA-16 genotypes were observed among 59.6% (68/114) B. melitensis isolates from Zhejiang and strains from other 21 different provinces, especially northern provinces, China. The analysis highlighted the imported nature of the strains from all over the northern provinces with a dominant part from the developed areas of animal husbandry. These data revealed a potential transmission pattern of brucellosis in this region, due to introduced infected sheep leading to a brucellosis outbreak epidemic, and eventually causing multiple laboratory infection events. It is urgent to strengthen the inspection and quarantine of the introduced animals.

Molecular identification of Brucella species and biovars associated with animal and human infection in Iran

Brucellosis is a costly contagious disease of human, domestic and wild animals. It is a serious health problem in Iran causing significant economic losses therefore, control approaches to prevent its spread are of great importance. In Iran, the species and biovars of virulent Brucella species are still under-reported due to the inadequate diagnostic protocols and insufficient laboratory facilities. The objective of this study was to characterize Brucella isolates obtained from passive animal and human surveillance in Iran from 2011 to 2018 in order to understand the current epidemiological situation of the disease. A total of 419 samples (milk, blood, cerebrospinal fluid, abomasum content and aborted fetus tissues) were collected from 65 cases/case series (human and animals) and examined bacteriologically. The initially identified Brucella isolates were further characterized using phenotypic and molecular approaches. All recovered isolates were either B. abortus or B. melitensis. The infection in sheep appeared to be exclusively associated with B. melitensis, but both B. abortus and B. melitensis were common in bovine samples. Samples from one sheep and one goat were confirmed to be infected by the B. melitensis vaccine strain Rev1. In spite of B. abortus burden in animals (14 cases in cattle and camel), brucellosis in human was predominantly associated with B. melitensis (15 cases). The results confirmed that B. melitensis biovar 1 and B. abortus biovar 3 remain the most prevalent biovars in Iran. This report builds a picture of the significance of different Brucella species in different hosts in Iran and provides applicable information for the healthcare professionals about the public health risks of brucellosis and relevant preventive strategies.

Multiple Locus Variable-Number Tandem-Repeat and Single-Nucleotide Polymorphism-Based Brucella Typing Reveals Multiple Lineages in Brucella melitensis Currently Endemic in China

Brucellosis is a worldwide zoonotic disease caused by Brucella spp. In China, brucellosis is recognized as a reemerging disease mainly caused by Brucella melitensis specie. To better understand the currently endemic B. melitensis strains in China, three Brucella genotyping methods were applied to 110 B. melitensis strains obtained in past several years. By MLVA genotyping, five MLVA-8 genotypes were identified, among which genotypes 42 (1-5-3-13-2-2-3-2) was recognized as the predominant genotype, while genotype 63 (1-5-3-13-2-3-3-2) and a novel genotype of 1-5-3-13-2-4-3-2 were second frequently observed. MLVA-16 discerned a total of 57 MLVA-16 genotypes among these Brucella strains, with 41 genotypes being firstly detected and the other 16 genotypes being previously reported. By BruMLSA21 typing, six sequence types (STs) were identified, among them ST8 is the most frequently seen in China while the other five STs were firstly detected and designated as ST137, ST138, ST139, ST140, and ST141 by international multilocus sequence typing database. Whole-genome sequence (WGS)-single-nucleotide polymorphism (SNP)-based typing and phylogenetic analysis resolved Chinese B. melitensis strains into five clusters, reflecting the existence of multiple lineages among these Chinese B. melitensis strains. In phylogeny, Chinese lineages are more closely related to strains collected from East Mediterranean and Middle East countries, such as Turkey, Kuwait, and Iraq. In the next few years, MLVA typing will certainly remain an important epidemiological tool for Brucella infection analysis, as it displays a high discriminatory ability and achieves result largely in agreement with WGS-SNP-based typing. However, WGS-SNP-based typing is found to be the most powerful and reliable method in discerning Brucella strains and will be popular used in the future.

Molecular epidemiological characterization of Brucella isolates from sheep and yaks in northwest China

Animal brucellosis is a re-emerging disease in China with high prevalence in the northwest region. A total of 66 isolates of Brucella were recovered from sheep and yaks in the Inner Mongolia, Xinjiang, Qinghai and Gansu provinces of northwest China in 2015 and 2016. Using classical biotyping and the Brucella AMOS PCR assay, all isolates were identified as Brucella melitensis biovar 3 (n = 58), B. melitensis biovar 1 (n = 1), Brucella abortus (n = 5), or Brucella suis biovar 3 (n = 2), and B. melitensis biovar 3 was found to be mainly responsible for sheep brucellosis in northwest China. Multilocus variable-number tandem-repeat analysis (MLVA) was used to identify the epidemiological relationships among the isolates and to assess their genetic diversity. Multilocus variable-number tandem-repeat analysis-16 identified 46 genotypes in these populations, including 37 unique and nine shared genotypes. Multilocus variable-number tandem-repeat analysis-11 showed that 71% of the isolates (47 of 66) were genotype 116 (1-5-3-13-2-2-3-2-4-41-8), a characteristic subgroup of the East Mediterranean group, showing that isolates from different geographical areas exhibit similar epidemiological characteristics in different regions and may be epidemiologically linked. Multilocus variable-number tandem-repeat analysis-11 also revealed that an isolate from Inner Mongolia had a novel genotype, 369 (1-5-3-13-2-2-3-2-7-41-8). Multilocus variable-number tandem-repeat analysis-16 genotyping of northwest China Brucella isolates allows a better understanding of the epidemiology of animal brucellosis in this region. This study is the first analysis of B. melitensis in Gansu province, and the results confirmed that in this province, isolates of this species are disorderly and unsystematic.

Species identification and molecular typing of human Brucella isolates from Kuwait

Brucellosis is a zoonotic disease of major concern in Kuwait and the Middle East. Human brucellosis can be caused by several Brucella species with varying degree of pathogenesis, and relapses are common after apparently successful therapy. The classical biochemical methods for identification of Brucella are time-consuming, cumbersome, and provide information limited to the species level only. In contrast, molecular methods are rapid and provide differentiation at intra-species level. In this study, four molecular methods [16S rRNA gene sequencing, real-time PCR, enterobacterial repetitive intergenic consensus (ERIC)-PCR and multilocus variable-number tandem-repeat analysis (MLVA)-8, MLVA-11 and MLVA-16 were evaluated for the identification and typing of 75 strains of Brucella isolated in Kuwait. 16S rRNA gene sequencing of all isolates showed 90-99% sequence identity with B. melitensis and real-time PCR with genus- and species- specific primers identified all isolates as B. melitensis. The results of ERIC-PCR suggested the existence of 75 ERIC genotypes of B. melitensis with a discriminatory index of 0.997. Cluster classification of these genotypes divided them into two clusters, A and B, diverging at ~25%. The maximum number of genotypes (n = 51) were found in cluster B5. MLVA-8 analysis identified all isolates as B. melitensis, and MLVA-8, MLVA-11 and MLVA-16 typing divided the isolates into 10, 32 and 71 MLVA types, respectively. Furthermore, the combined minimum spanning tree analysis demonstrated that, compared to MLVA types discovered all over the world, the Kuwaiti isolates were a distinct group of MLVA-11 and MLVA-16 types in the East Mediterranean Region.

Molecular epidemiology of Brucella abortus isolated from cattle in Brazil, 2009-2013

The aims of the present study were to genotype Brucella abortus strains isolated from cattle in Brazil between 2009 and 2013, and to analyze their distribution to support the Programa Nacional de Controle e Erradicação de Brucelose e Tuberculose (PNCEBT) (National Brucellosis and Tuberculosis Control and Eradication Program). One hundred forty B. abortus strains isolated from cattle in Brazil between 2009 and 2013 were genotyped using a set of 18 variable number of tandem repeats (VNTR) (MLVA16+HOOF-Print 3 and 4). The multiple locus VNTR analysis (MLVA) composed by eight markers (MLVA8) revealed eight different genotypes among B. abortus strains, including five previously described and three new ones. Analysis of the MLVA16 loci revealed fifty-eight distinct genotypes, from which three were identical, thirty-eight were considered very close, and seventeen were considered distant compared to those previously described and deposited in MLVAbank. Analysis of the HOOF-Prints 3 and 4 revealed the larger number of different alleles among all VNTR assessed, exhibiting maximum resolution when associated with MLVA16 markers. This study also provides insights on the genotypes of B. abortus circulating in Brazil, which certainly contribute for the better understanding of the epidemiology and control of bovine brucellosis in the country. Moreover, our data showed a high genetic diversity among the B. abortus strains isolated between 2009 and 2013, and a close relationship among these strains and Brazilian B. abortus deposited by MLVAbank.

MLVA Genotyping Characteristics of Human Brucella melitensis Isolated from Ulanqab of Inner Mongolia, China

Brucellosis is a serious public health problem in Ulanqab, which is a region located in the middle of the Inner Mongolia Autonomous Region adjacent to Shanxi and Hebei provinces. The disease is prevalent in both the latter provinces and Ulanqab with the highest prevalence of brucellosis occurring in Inner Mongolia. The MLVA-16 scheme is a genotyping tool for assessing genetic diversity and relationships among isolates. Moreover, this genotyping tool can also be applied to epidemiological trace-back investigations. This study reports the occurrence of at least two B. melitensis biovars (1 and 3) in Ulanqab, encompassing 22 and 94 isolates, respectively. B. melitensis biovar 3 was the predominant biovar in the area examined. Panel 1 (MLVA-8) identified three genotypes (42, 63, and 114), with genotype 42 (n = 101) representing 87% of the tested strains. MLVA-11 identified eight genotypes (116, 111, 297, 291, and 342-345) from 116 of the analyzed isolates. All of these isolates were identified as belonging to the East Mediterranean group. Genotype 116 (n = 94) was the predominant genotype and represented 81% of the isolates. The isolates pertaining to this genotype were distributed throughout most of Ulanqab and neighboring regions. The MLVA-16 scheme showed the presence of 69 genotypes, with 46 genotypes being represented by single isolates. This analysis revealed that Ulanqab brucellosis cases had epidemiologically unrelated and sporadic characteristics. The remaining 23 genotypes were shared (between a total of 70 isolates) with each genotype being represented by two to eight isolates. These data indicate that these cases were epidemiologically related. MLVA genotyping confirmed the occurrence of a multipoint outbreak epidemic and intrafamilial brucellosis. Extensive genotype-sharing events were observed among isolates from different regions of Ulanqab and from other provinces of China. These findings suggest either a lack of control of animal movement between different regions or the circulation of contaminated animal products in the market. Our study is the first comprehensive genotyping and genetic analysis of B. melitensis in Ulanqab. We believe that this study will help to improve the effectiveness of brucellosis control programs.