Molecular epidemiological investigation of Brucella melitensis circulating in Mongolia by MLVA16

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.

Emerging diversity and ongoing expansion of the genus Brucella

Remarkable genetic diversity and breadth of host species has been uncovered in the Brucella genus over the past decade, fundamentally changing our concept of what it means to be a Brucella. From ocean fishes and marine mammals, to pond dwelling amphibians, forest foxes, desert rodents, and cave-dwelling bats, Brucella have revealed a variety of previously unknown niches. Classical microbiological techniques have been able to help us classify many of these new strains but at times have limited our ability to see the true relationships among or within species. The closest relatives of Brucella are soil bacteria and the adaptations of Brucella spp. to live intracellularly suggest that the genus has evolved to live in vertebrate hosts. Several recently discovered species appear to have phenotypes that are intermediate between soil bacteria and core Brucella, suggesting that they may represent ancestral traits that were subsequently lost in the traditional species. Remarkably, the broad relationships among Brucella species using a variety of sequence and fragment-based approaches have been upheld when using comparative genomics with whole genomes. Nonetheless, genomes are required for fine-scale resolution of many of the relationships and for understanding the evolutionary history of the genus. We expect that the coming decades will reveal many more hosts and previously unknown diversity in a wide range of environments.

Human brucellosis epidemiology in the pastoral area of Hulun Buir city, Inner Mongolia autonomous region, China, between 2003 and 2018

Human brucellosis represents a serious public health concern in Hulun Buir and requires a comprehensive epidemiologic analysis to define adapted control measures. The present study describes the case numbers, constituent ratios and incidence rate of human brucellosis. Conventional biotyping, that is abortus, melitensis, ovis and suis (AMOS)-PCR and multi-locus variable-number tandem repeat analysis (MLVA) were used to characterize the Brucella strains. Between 2003 and 2018, a total of 23,897 human brucellosis cases were reported, with an incidence rate of 56.03/100,000, which is 20 times higher than the country's average incidence. This incidence rate increased year after year, culminating in 2005 and decreased between 2011 and 2018. Because Hulun Buir relies on a nomadic livestock rearing system, brucellosis spreads easily among different animal species and humans. In Xin Barag Left Banner and Xin Barag Right Banner, the incidence rates were, respectively, 226.54/100,000 and 199.10/100,000, exceeding those observed in other areas. Most of the cases occurred in the 25- to 45-year-old group, accounting for 65.74% of the cases (15,709/23,897), and among farmers, accounting for 66.71% (15,942/23,897). The male to female incidence ratio was 2.67:1. The higher incidence in younger people and the large gender ratio reflected the unique traditional production and lifestyle of nomads. Most reported cases were observed from April to June, indicating that more than 40% of the cases were related to the delivery of domestic livestock. The biotyping showed that the 44 isolated strains were all B. melitensis, including 12 Brucella melitensis biovar (bv) 1 and 32 B. melitensis bv. 3. The strains displayed a genetic similarity of 80%-100%. Our hypothesis is that human brucellosis outbreak in this region may be originating from a limited source of infection, so further investigation is necessary. The epidemic situation of human brucellosis in Hulun Buir is extremely serious, strengthened surveillance and control in animals' brucellosis should be priority.

Investigation of Genetic Relatedness of Brucella Strains in Countries Along the Silk Road

In this study, MLVA (multiple-locus variable-number tandem repeat analysis) genotype data of Brucella strains from 11 countries along the Silk Road were downloaded from the MLVAbank. MLVA data of strains were applied to the constructed Minimum Spanning Tree to explore the species/biovars distribution, geographic origins, and genetic relationships of the strains analyzed. Moreover, whole-genome sequencing–single-nucleotide polymorphism (WGS-SNP) phylogenetic analysis of the genome of Brucella melitensis strains from GenBank was performed to discriminate the relatedness of strains further and investigate the transmission pattern of B. melitensis brucellosis. A total of 1,503 Brucella strains were analyzed in this study: 431 Brucella abortus strains (29.8%), 1,009 B. melitensis strains (65.7%), and 63 Brucella suis strains (4.5%). B. melitensis biovar 3 was the dominant species and was shown to be widespread in all of the examined regions, suggesting that the prevention and surveillance of the B. melitensis population are a main challenge in these countries. A wide host spectrum was observed for this Brucella population; many animal reservoirs are a potential reason for the continuous brucellosis circulation in these countries. Although the B. abortus strains from the examined regions had common geographic origins, only a few shared genotypes were observed in different countries. These data revealed that the majority of B. abortus strains were spreading within the national borders. However, the B. melitensis strains from Italy originated from a Western Mediterranean lineage; strains from the other 10 countries originated from Eastern Mediterranean lineage, and this lineage was shared by strains from three to nine different countries, suggesting that the introduction and reintroduction of the disease in the 10 countries might have occurred in the past. Furthermore, the most shared MLVA-16 genotypes were formed in the B. melitensis strains from China, Kazakhstan, and Turkey, suggesting that the introduction and trade in sheep and goats have occurred frequently in these countries. WGS-SNP analysis showed that the B. melitensis in this study originated from the Malta (Italy) region. According to their territorial affiliation between four clade strains from these countries in genotype B, the absence of a clear differentiation suggests that strains continuously expand and spread in countries along with Silk Road. Active exchange and trade of animals (sheep and goats) among these countries are reasonable explanations. B. suis strains from different nations showed unique geographic origins and epidemiological characteristics. Therefore, there is an urgent need for the control of transfer and trade of infected sheep (goats) in countries along the Silk Road, namely, the strengthening of the entry–exit quarantine of sheep and goats and improvements in the diagnosis of animal brucellosis.

Immunological pathways of macrophage response to Brucella ovis infection

As the molecular mechanisms of Brucella ovis pathogenicity are not completely clear, we have applied a transcriptome approach to identify the differentially expressed genes (DEGs) in RAW264.7 macrophage infected with B. ovis. The DEGs related to immune pathway were identified by Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) functional enrichment analysis. Quantitative real-time PCR (qRT-PCR) was performed to validate the transcriptome sequencing data. In total, we identified 337 up-regulated and 264 down-regulated DEGs in B. ovis-infected group versus mock group. Top 20 pathways were enriched by KEGG analysis and 20 GO by functional enrichment analysis in DEGs involved in the molecular function, cellular component, and biological process and so on, which revealed multiple immunological pathways in RAW264.7 macrophage cells in response to B. ovis infection, including inflammatory response, immune system process, immune response, cytokine activity, chemotaxis, chemokine-mediated signaling pathway, chemokine activity, and CCR chemokine receptor binding. qRT-PCR results showed Ccl2 (ENSMUST00000000193), Ccl2 (ENSMUST00000124479), Ccl3 (ENSMUST00000001008), Hmox1 (ENSMUST00000005548), Hmox1 (ENSMUST00000159631), Cxcl2 (ENSMUST00000075433), Cxcl2 (ENSMUST00000200681), Cxcl2 (ENSMUST00000200919), and Cxcl2 (ENSMUST00000202317). Our findings firstly elucidate the pathways involved in B. ovis-induced host immune response, which may lay the foundation for revealing the bacteria–host interaction and demonstrating the pathogenic mechanism of B. ovis.

Molecular Characteristics of Brucella Isolates Collected From Humans in Hainan Province, China

Brucellosis has been reported in several regions of Hainan Province, but the extent of the disease has not been fully elucidated. Conventional biotyping methods, multiple locus variable number tandem repeats analysis (MLVA), and single-nucleotide polymorphisms (SNPs) from draft genome sequencing were employed to characterize the strains. There were four biovars (Brucella melitensis bv. 1, 2, and 3 and Brucella suis bv. 3) detected, which showed that the biovar diversity of Brucella in Hainan is higher than in other areas of China. Both B. melitensis bv. 3 and B. suis bv. 3 were dominant species and showed epidemiology patterns that were compatible with both southern and northern China. Eight of MLVA-11 genotypes were known (31, 111, 116, 120, 136, 291, 297, and 345), and the remaining seven were novel (HN11-1 to HN11-7); these data showed that Brucella strains in this study had multiple geographic origins and exhibited characteristics of origin and evolution of co-existing imported and Hainan specific lineage. A total of 41 strains were found, belonging to 37 unique genotypes that each represented a single strain, which suggests that these strains were not directly related epidemiologically and indicates that the epidemic characteristics of human brucellosis in Hainan was dominated by sporadic strains. The high HGDI values were observed in MLVA-8, MLVA-11, and MLVA-16 among two species, suggesting considerable genetic diversity among these species. MST is characterized based on MLVA-16 that was found both throughout China and on a global level and showed that strains of this study had significant genetic differences with strains from many parts of the globe and seemingly represent a unique genetic lineage. Whole-genome SNP analysis showed that four B. melitensis were closely related to strains from China's northern provinces, and the source of infection was partly of human brucellosis in this province that may have been from these regions. The B. suis were closely related to strains from the United States, and further investigation of the transportation of animals, such as pigs, is needed to elucidate the origins of these strains.

Reasons behind the epidemiological situation of brucellosis in the Republic of Kazakhstan

Brucellosis is endemic in the Republic of Kazakhstan, particularly in agro-pastoral areas. The purpose of this research is to study the epidemiological situation of brucellosis in livestock recorded in the Republic of Kazakhstan, and to identify the reasons why anti-brucellosis measures were not effective. The research was performed on statistical data provided by the Republican Veterinary Laboratory (RVL), the Ministry of Agriculture of the Republic of Kazakhstan. The analysis touched upon the prevalence of Brucellosis in cattle and small ruminants (sheep and goats) in 13 regions of the Republic of Kazakhstan in 2012-2016. Aside from that, Research Institute for Biological Safety Problems conducted screening assays that involved 11,889 samples of blood and tissues from said animals. The risks of developing Brucellosis were assessed for each particular region. The comparison of studies conducted in 2012-2016 reveals an increase in the prevalence of Brucellosis in cattle in the following regions: West Kazakhstan, Karaganda Region and Pavlodar Region. For small ruminants, growing prevalence was observed in the Kostanay Region, Jambyl Region and Almaty Region. Between 2014 and 2016, the incidence rate had a growing trend, with a high in 2014 and 2015. The lowest prevalence rate during the following years (2012-2016) was in the Mangystau Region. The Enzyme-Linked Immunosorbent Assay (ELISE) test applied in the research among other tests provided the best results. The main risk factors involve epidemiology and sanitary measures, which are undertaken in the Republic of Kazakhstan, geography of the region with focuses of infection, and randomness of spread.

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.