<figure><img src="image_1.png" /></figure>

Check for
updates

# Genetic Diversity of *Brucella melitensis* in Kazakhstan in Relation to World-Wide Diversity

Elena Shevtsova<fnref n="1" />, Gilles Vergnaud², Alexandr Shevtsov<fnref n="1" />*, Alexandr Shustov<fnref n="1" />, Kalysh Berdimuratova<fnref n="1" />, Kasim Mukanov<fnref n="1" />, Marat Syzdykov³, Andrey Kuznetsov³, Larissa Lukhnova³, Uinkul Izbanova³, Maxim Filipenko<fnref n="1" />,⁴,⁵ and Yerlan Ramankulov<fnref n="1" />,⁶

## OPEN ACCESS

**Edited by:**
Axel Cloeckaert,
Institut National de la Recherche
Agronomique (INRA), France

**Reviewed by:**
Herbert Tomaso,
Friedrich-Loeffler-Institut, Germany
Jacques Xavier Godfroid,
UiT The Arctic University of Norway,
Norway

**\*Correspondence:**
Alexandr Shevtsov
ncbshevtsov@gmail.com

**Specialty section:**
This article was submitted to
Infectious Diseases,
a section of the journal
Frontiers in Microbiology

**Received:** 03 March 2019
**Accepted:** 31 July 2019
**Published:** 13 August 2019

**Citation:**
Shevtsova E, Vergnaud G, Shevtsov A, Shustov A, Berdimuratova K, Mukanov K, Syzdykov M, Kuznetsov A, Lukhnova L, Izbanova U, Filipenko M and Ramankulov Y (2019) Genetic Diversity of Brucella melitensis in Kazakhstan in Relation to World-Wide Diversity. *Front. Microbiol.* 10:1897. doi: 10.3389/fmicb.2019.01897

We describe the genetic diversity of 1327 *Brucella* strains from human patients in Kazakhstan using multiple-locus variable-number tandem repeat (VNTR) analysis (MLVA). All strains were assigned to the *Brucella melitensis* East Mediterranean group and clustered into 16 MLVA11 genotypes, nine of which are reported for the first time. MLVA11 genotype 116 predominates (86.8%) and is present all over Kazakhstan indicating existence and temporary preservation of a “founder effect” among *B. melitensis* strains circulating in Central Eurasia. The diversity pattern observed in humans is highly similar to the pattern previously reported in animals. The diversity observed by MLVA suggested that the epidemiological status of brucellosis in Kazakhstan is the result of the introduction of a few lineages, which have subsequently diversified at the most unstable tandem repeat loci. This investigation will allow to select the most relevant strains for testing these hypotheses via whole genome sequencing and to subsequently adjust the genotyping scheme to the Kazakhstan epidemiological situation.

**Keywords:** *Brucella melitensis*, Kazakhstan, genotyping, genetic diversity, multiple-locus variable-number tandem repeat (VNTR) analysis (MLVA)

## INTRODUCTION

Brucellosis is a zoonotic infection affecting many mammals including marine species as well as humans. Infection is endemic in many countries of the world, the etiological agent being a Gram-negative bacteria of the genus *Brucella*. Among the 12 species currently proposed in this genus, *Brucella melitensis*, *B. abortus*, and *B. suis* are highly dangerous for humans and cause disease with severe complications and chronic process (Meyer, 1990; Godfroid et al., 2011). Rare cases of infection of people with other *Brucella* species are also recorded (Ficht, 2010). Despite its low mortality rates, brucellosis is a very important public health problem in Kazakhstan. Husbandry suffer direct economic losses from brucellosis due to reduced productivity, culling of livestock and costs of associated measures. In spite of the control strategies, brucellosis remains a major economic problem for agriculture, with an annual direct cost of more than $45 million and a social problem in health care with an average disability-adjusted life year of 0.5 per case (Charypkhan et al., 2019).

<footnote n="1">¹ National Center for Biotechnology, Nur-Sultan, Kazakhstan, ² Institute for Integrative Biology of the Cell (I2BC), CEA, CNRS, Univ. Paris-Sud, Université Paris-Saclay, Gif-sur-Yvette, France, ³ Kazakh Scientific Center for Quarantine and Zoonotic Diseases, Almaty, Kazakhstan, ⁴ Institute of Chemical Biology and Fundamental Medicine, Novosibirsk, Russia, ⁵ Synthetic Biology Department, Novosibirsk State University, Novosibirsk, Russia, ⁶ School of Science and Technology, Nazarbayev University, Nur-Sultan, Kazakhstan</footnote>
In addition, aerosol transmission and low infectious doses (10–100 bacteria) potentially allow using *Brucella* in acts of bioterrorism (Pappas et al., 2006b; Neubauer, 2010; Doganay and Doganay, 2013).

Control of animal infections is essential to prevent human infections. In developed countries, epidemiological measures were successfully undertaken to eliminate brucellosis in domestic animals, resulting in that all human cases are recorded among migrants, or associated with traveling to endemic areas, consumption of imported products, and in rare cases contacts with wild animals (Gwida et al., 2012). On the contrary, in many countries of the Mediterranean region, Southern and Central America, Africa, Asia, Arabian Peninsula, Indian subcontinent, Eastern Europe, and Middle East, the brucellosis epidemiological situation remains complex allowing to describe it as a continuing epidemic (Nicoletti, 2010).

Given the zoonotic nature of brucellosis, identification of sources of the infection and tracking of transmission paths are of high importance for epidemiologic surveillance. Reliable discrimination among *Brucella* species and biovars are important for treatment, due to differences in pathogenicity between various *Brucella*. DNA-based genotyping methods are preferred because molecular methods minimize the risks of infection of laboratory personnel and phenotypic methods have a much lower discriminatory power. For *Brucella*, a number of genotyping methods have been proposed including AMOS-PCR, Brucel ladder, tandem-repeats polymorphism [multiple-locus VNTR analysis (MLVA)] and whole genome SNP analysis (wgSNP). wgSNP analysis is by far the most powerful and robust approach, but its current cost does not allow yet to use it as a routine first line assay (Ledwaba et al., 2019). MLVA typing is a compromise in terms of discriminatory power, cost and speed to provide an overview of genetic diversity and assist the selection of strains for whole genome sequencing (Scholz and Vergnaud, 2013). Availability of freely accessible databases with MLVA genotypes allow to track the spreading of specific *Brucella* biovars globally (Kattar et al., 2008; Kilic et al., 2011; Ferreira et al., 2012; Vergnaud et al., 2018). Fundamental for this work is genotyping and consolidation of data from all endemic regions.

Kazakhstan is located in the center of Eurasia. It covers a territory of 2.725 million km² for 18 million inhabitants. The country shares land borders with Russia, China, Kyrgyzstan, Uzbekistan, Turkmenistan. With regard to brucellosis, Kazakhstan is a hyperendemic area considering the very high incidence rates in the human population and farm animals (Pappas et al., 2006a). According to evaluations in the Former Soviet Union (FSU), more than 40% of newly diagnosed brucellosis cases are from Kazakhstan (Studentsov, 1975). From 1999 to 2016, 38,557 new cases of human brucellosis were recorded, with annual counts of cases ranging from 1443 (2014) to 3596 (2004) (Syzdykov et al., 2013; Shevtsova et al., 2016; Daugaliyeva et al., 2018). Within these 18 years, the majority of cases (89.5%) occurred in southern and eastern Kazakhstan (Almaty, Zhambyl, Kyzylorda, East-Kazakhstan, and South-Kazakhstan regions). 75% of the cases of newly diagnosed brucellosis were in the age group 10–39 years, and 85.8% were rural residents. Most strains recovered from patients

were identified as *B. melitensis* (Mizanbayeva et al., 2009). The epidemiological situation correlates with occurrence of brucellosis in small ruminants, e.g., the mentioned geographic regions hold more than 70% of all identified seropositive animals among small ruminants (data collected during 2008–2015) (Espembetov et al., 2017). From 2006 to 2015, the average level of seropositive small ruminants ranged from 0.15 to 0.3% with an increase in 2008 (0.65%) and 2013 (0.5%). The fluctuation from 0.01 to 2.5% of the percentage of seropositive animals at the district level is uncorrelated with the level of infection in humans (Daugaliyeva et al., 2018; Charypkhan et al., 2019). In cattle, the incidence of seropositive animals varied from 0.3 in 2006 to 0.6% in 2013–2015. The high incidence of brucellosis reported in cattle during 2009–2011 (1.85%) is most probably an artifact originating from a change in diagnostic procedures used in Kazakhstan. During this period, all cattle was screened using only an indirect enzyme-linked immunosorbent assay (iELISA) and classical methods such as the rose Bengal test (RBT), serum agglutination test (SAT), and complement fixation test (CFT) were not used for confirmation of the brucellosis diagnosis, until the low specificity of the iELISA became evident. Starting from 2012, the classical methods were again made compulsory for confirmation of brucellosis diagnosis which resulted in reduction of the reported incidence (Espembetov et al., 2017). In addition to the main reservoir hosts, brucellosis caused by *B. melitensis* or *B. abortus* is registered in camels, horses, and dogs (Zheludkov and Tsirelson, 2010). *Brucella* antibodies have been detected in more than 60 wild mammals (Shevtsova et al., 2016).

The high morbidity and large territory of Kazakhstan make it essential to describe the genetic diversity and distribution of genotypes among circulating strains of *Brucella* in order to support efforts for disease control. Previous reports have described the genetic diversity of *Brucella* isolated from animals in Kazakhstan (Shevtsov et al., 2015; Shevtsova et al., 2016; Daugaliyeva et al., 2018; Yespembetov et al., 2019). The purpose of the present study was to determine the genetic diversity of the strains infecting humans at the whole country scale.

# MATERIALS AND METHODS

## Ethics Statement

A formal institutional ethical review process and approval was not required for this study as all strains investigated here were collected as a part of standard clinical investigation of patients with suspected brucellosis and the strains were anonymized. All *Brucella* strains used in this study were obtained from the collection maintained by the Masgut Aykimbayev Kazakh Scientific Center of Quarantine and Zoonotic Diseases (KSCQZD) in accordance to Kazakhstan regulations.

## Clinical Strains Characterization

A total of 1383 *Brucella* strains were isolated from patients seropositive to brucellosis antigen (agglutination test 1:200 and higher) during 2015–2016 in 13 regional infectious diseases hospitals. The hemocultures were produced by the two-phase method proposed by Castaneda (Castaneda, 1947;
Yagupsky, 1999). All strains were Gram negative, agglutinated with polyvalent brucellosis serum, had oxidase and catalase activity, did not produce H₂S, synthesized urease, and were capable of growing in plain atmospheric conditions. Biotyping was performed by KSCQZD by routine test (Al Dahouk et al., 2003).

## DNA Preparation, Genotyping, and Data Analysis

Bacterial culture was scrapped from the surface of the solid agar medium, resuspended in 500 µl *Tris* 10 mM EDTA 1 mM pH 8 (TE) buffer and inactivated by adding an equal volume of chloroform. DNA was isolated using the QIAamp DNA Mini Kit (Qiagen, United States). Species identification was confirmed

by the multiplex PCR “Bruce-ladder” assay (Lopez-Goni et al., 2008). The 16S rDNA fragment was sequenced to confirm taxonomic identification of the genus and to exclude presence of extraneous microflora.

Markers proposed by Le Flèche et al. (2006) and Al Dahouk et al. (2007) were used for MLVA. The 16 loci are commonly divided into three panels. The first panel (also called MLVA8) consists of the eight tandem repeat loci with the largest repeat units, Bruce06, Bruce08, Bruce12, Bruce42, Bruce43, Bruce45, and Bruce55. The eight other tandem repeats have smaller repeat units up to eight base-pairs. The second panel (panel 2A) includes Bruce18, Bruce19, and Bruce21. Panel 2B includes the last five loci: Bruce04, Bruce07, Bruce09, Bruce16, and Bruce30. The combination of panels MLVA8 and 2A is called MLVA11, while the combination of all three panels (16 loci) is designated MLVA16. The MLVA11 panel allows for species identification in the *Brucella* genus and grouping of strains in accordance to the global geographic distribution. The panel 2B loci are highly discriminatory and their combination with MLVA11 is used in tracking local outbreaks. New variable-number tandem repeat (VNTR) alleles were confirmed by sequencing the PCR product.

Primers for multiplex PCR and their combinations were used as previously described (Garofolo et al., 2013; Shevtsov et al., 2015). PCR amplification products were diluted 1/70, and aliquots (3 µl) of the diluted samples were run on a capillary electrophoresis machine (DNA Analyzer 3730xl, Applied Biosystems, Japan), using the POP7 polymer and LIZ 1200 size standard. Size analysis of VNTR repeats was performed using GeneMapper 4.1 (Applied Biosystems).

To visualize clustering relations, maximum parsimony analysis dendrograms were constructed using BioNumerics 7.6 (Applied Maths, Sint-Martens-Latem, Belgium). The Hunter and Gaston diversity index (HGDI) was used to describe discriminatory capacity of each locus, as well as of the panels MLVA8, MLVA11, and MLVA16 (Hunter and Gaston, 1988). Computations were performed using the Internet resources at www.hpa-bioinfotools.org.uk/cgi-bin/DICI/DICI.pl. Geographic coordinates (longitude, latitude) were deduced from available geographic origin data using the BioNumerics geographical plugin. The geographic location map was drawn using QGIS 3.4.8-Madeira (QGIS Geographic Information System. Open Source Geospatial Foundation)<fnref n="1" />.

<table><caption>**TABLE 1** | Hunter and Gaston diversity index (HGDI) for loci and MLVA panels for 1327 *B. melitensis* strains from human patients in Kazakhstan.</caption><thead><tr><th>Locus</th><th>Diversity index</th><th>Confidence interval</th><th>Number of genotypes observed</th><th>Max(pI)*</th></tr></thead><tbody><tr><td>Bruce06</td><td>0.000</td><td>0.000-0.006</td><td>1</td><td>1.000</td></tr><tr><td>Bruce08</td><td>0.025</td><td>0.020-0.031</td><td>3</td><td>0.987</td></tr><tr><td>Bruce11</td><td>0.000</td><td>0.000-0.006</td><td>1</td><td>1.000</td></tr><tr><td>Bruce12</td><td>0.000</td><td>0.000-0.006</td><td>1</td><td>1.000</td></tr><tr><td>Bruce42</td><td>0.000</td><td>0.000-0.006</td><td>1</td><td>1.000</td></tr><tr><td>Bruce43</td><td>0.142</td><td>0.130-0.153</td><td>4</td><td>0.924</td></tr><tr><td>Bruce45</td><td>0.000</td><td>0.000-0.006</td><td>1</td><td>1.000</td></tr><tr><td>Bruce55</td><td>0.000</td><td>0.000-0.006</td><td>1</td><td>1.000</td></tr><tr><td>MLVA8</td><td>0.163</td><td>0.150-0.175</td><td>7</td><td>0.913</td></tr><tr><td>Bruce18</td><td>0.049</td><td>0.041-0.056</td><td>3</td><td>0.975</td></tr><tr><td>Bruce19</td><td>0.052</td><td>0.044-0.059</td><td>6</td><td>0.974</td></tr><tr><td>Bruce21</td><td>0.000</td><td>0.000-0.006</td><td>1</td><td>1.000</td></tr><tr><td>MLVA11</td><td>0.242</td><td>0.228-0.256</td><td>16</td><td>0.868</td></tr><tr><td>Bruce04</td><td>0.770</td><td>0.764-0.775</td><td>11</td><td>0.349</td></tr><tr><td>Bruce07</td><td>0.292</td><td>0.277-0.306</td><td>6</td><td>0.835</td></tr><tr><td>Bruce09</td><td>0.016</td><td>0.012-0.021</td><td>4</td><td>0.992</td></tr><tr><td>Bruce16</td><td>0.786</td><td>0.782-0.791</td><td>12</td><td>0.320</td></tr><tr><td>Bruce30</td><td>0.632</td><td>0.624-0.640</td><td>6</td><td>0.503</td></tr><tr><td>MLVA16</td><td>0.989</td><td>0.988-0.989</td><td>310</td><td>0.050</td></tr></tbody></table>

<table><caption>**TABLE 2** | HGDI for the panel 2B for East Mediterranean *B. melitensis* from various countries.</caption><thead><tr><th>Country</th><th>Bruce04</th><th>Bruce07</th><th>Bruce09</th><th>Bruce16</th><th>Bruce30</th><th>Number of strains</th></tr></thead><tbody><tr><td>Kazakhstan</td><td>0.770</td><td>0.292</td><td>0.016</td><td>0.786</td><td>0.632</td><td>1327</td></tr><tr><td>France</td><td>0.851</td><td>0.720</td><td>0.858</td><td>0.880</td><td>0.540</td><td>479</td></tr><tr><td>Spain</td><td>0.807</td><td>0.778</td><td>0.371</td><td>0.814</td><td>0.680</td><td>403</td></tr><tr><td>China</td><td>0.750</td><td>0.114</td><td>0.109</td><td>0.834</td><td>0.754</td><td>303</td></tr><tr><td>Turkey</td><td>0.762</td><td>0.487</td><td>0.195</td><td>0.782</td><td>0.737</td><td>233</td></tr><tr><td>Italy</td><td>0.798</td><td>0.784</td><td>0.894</td><td>0.817</td><td>0.070</td><td>225</td></tr><tr><td>Portugal</td><td>0.760</td><td>0.737</td><td>0.321</td><td>0.777</td><td>0.451</td><td>149</td></tr><tr><td>Peru</td><td>0</td><td>0.164</td><td>0.689</td><td>0.18</td><td>0</td><td>114</td></tr><tr><td>Greece</td><td>0.774</td><td>0.279</td><td>0.57</td><td>0.778</td><td>0.75</td><td>78</td></tr></tbody></table>

## RESULTS

A total of 1383 strains *Brucella* collected from 13 regions of Kazakhstan were studied. All strains could be assigned to the *Brucella* genus based on analysis of a fragment of the nucleotide sequence of 16S rDNA and to *B. melitensis* according to the Bruce-ladder profile. The biovar was established for 500 strains. Ninety-five (19%) were assigned to biovar 1, six (1.2%) to biovar 2 and 399 (79.8%) were assigned to biovar 3 (Supplementary Table S1). The 16 VNTR loci constituting the MLVA assay were successfully amplified in all DNAs. Five samples showed a mixed

<footnote>*Max(pI) = fraction of samples that have the most frequent repeat number in this locus (range 0.0-1.0).</footnote>

<footnote n="1">¹http://qgis.org</footnote>
<figure><img src="image_2.png" /></figure>

**FIGURE 1** | MLVA11 in 2900 East Mediterranean *B. melitensis* entries including 1580 from *Brucella* MLVAbank version 4.3. The 1327 strains from the present study are colored. The biovar assignment is indicated when known. The numbers are the MLVA11 genotypes. The genotypes are connected by branches with a length of one or two corresponding to the categorical distance between them.

pattern for the moderately variable Bruce11 locus. Twenty-seven samples demonstrated a mixed pattern in one of the highly variable loci constituting panel 2A and 2B: Bruce16 (15 isolates), Bruce18 (four isolates), Bruce30 (three isolates), Bruce19 and Bruce07 (two isolates each), Bruce04 (one isolate). Twenty-four additional isolates showed mixed patterns at more than two highly variable loci, in agreement with previous observations (Whatmore et al., 2006; Maquart et al., 2009). The 56 samples showing mixed patterns were excluded from further analyses. A new Bruce43 allele (206 bp, called 4) was observed in strains NCB#h-2315, NCB#h-3375 and NCB#h-3362. Two new Bruce19

alleles (106 bp, called 17 and 208 bp called 51) were observed in three strains (NCB#h-3378, NCB#h-3379, and NCB#h-3409) and one strain (NCB#h-2223), respectively. These new alleles were confirmed by sequencing of the PCR products.

The HGDI index of the 16 VNTR loci and of panels MLVA8, MLVA11, MLVA16 are indicated in **Table 1**. In panel MLVA11, only 5 among the 11 loci were polymorphic in the present collection. In panel 2B, all five loci were variable, with loci Bruce07 and Bruce09 being the least variable. The HGDI for MLVA8 and MLVA11 were 0.163 and 0.242, respectively. The HGDI for MLVA16 was 0.989. In countries with more than
<table><caption>**TABLE 3** | Distribution of genotypes by MLVA8 and MLVA11 in Kazakhstan.</caption><thead><tr><th>MLVA11 profile</th><th>MLVA8 genotype (number of strains/%)</th><th>MLVA11 genotype (number of isolates/%)</th></tr></thead><tbody><tr><td>1-5-3-13-2-2-3-2-4-41-8*</td><td>42 (1211/91.25)</td><td>116(1152/86.8)</td></tr><tr><td>1-5-3-13-2-2-3-2-5-41-8</td><td></td><td>180(29/2.18)</td></tr><tr><td>1-5-3-13-2-2-3-2-4-46-8</td><td></td><td>297(20/1.51)</td></tr><tr><td>1-5-3-13-2-2-3-2-3-41-8</td><td></td><td>342(2/0.15)</td></tr><tr><td>1-5-3-13-2-2-3-2-4-17-8</td><td></td><td>388(4/0.3)†</td></tr><tr><td>1-5-3-13-2-2-3-2-4-36-8</td><td></td><td>389(3/0.23)†</td></tr><tr><td>1-5-3-13-2-2-3-2-4-51-8</td><td></td><td>393(1/0.15)†</td></tr><tr><td>1-5-3-13-2-3-3-2-4-41-8</td><td>63 (89/6.71)</td><td>111(80/6.03)</td></tr><tr><td>1-5-3-13-2-3-3-2-4-46-8</td><td></td><td>386(6/0.45)†</td></tr><tr><td>1-5-3-13-2-3-3-2-4-36-8</td><td></td><td>392(2/0.15)†</td></tr><tr><td>1-5-3-13-2-3-3-2-4-43-8</td><td></td><td>394(1/0.15)†</td></tr><tr><td>1-3-3-13-2-2-3-2-4-41-8</td><td>62 (11/0.83)</td><td>114(11/0.83)</td></tr><tr><td>1-5-3-13-2-1-3-2-4-41-8</td><td>114 (7/0.53)</td><td>291(7/0.53)</td></tr><tr><td>1-4-3-13-2-2-3-2-4-41-8</td><td>115 (4/0.3)</td><td>387(4/0.3)†</td></tr><tr><td>1-5-3-13-2-4-3-2-4-41-8</td><td>196 (3/0.23) †</td><td>390(3/0.23)†</td></tr><tr><td>1-4-3-13-2-3-3-2-4-41-8</td><td>197(2/0.15) †</td><td>391(2/0.15)†</td></tr></tbody></table>

50 entries in MLVA bank, variable HGDI levels could be observed (Table 2). In Asian countries including Kazakhstan, panel2B loci Bruce07 and Bruce09 consistently showed low HGDI values. Locus Bruce03 is invariable in strains sampled from Africa and South America, and moderately variable in

Italy. In most European countries, all loci of the panel 2B are highly discriminatory.

All *B. melitensis* isolates isolated in Kazakhstan cluster within the East Mediterranean group congruently defined by MLVA, MLST and wgSNP (Whatmore et al., 2016; Vergnaud et al., 2018). In agreement with previous observations, all three *B. melitensis* biovars are observed within the East Mediterranean group. In the whole *Brucella* genus, 383 MLVA11 genotypes are defined, 100 of which belong to the East Mediterranean group (Figure 1). The 1327 *B. melitensis* strains from human patients isolated in Kazakhstan were clustered into 16 genotypes, nine of which numbered 386–394 have not been reported previously (Table 3). The most prevalent, genotype 116, was found in 86.8% of samples, genotype 111 in 6.03%, genotype 180 in 2.18%, and genotype 297 in 1.51%. The 12 other genotypes are rare (<1%). Genotype 116 is widespread in Kazakhstan whereas genotype 111 is detected in the hyperendemic areas in the southern part of Kazakhstan. The other genotypes have small prevalence in Kazakhstan and tend to be present in Southern Kazakhstan although genotypes 180 and 297 were also recovered in Northern Kazakhstan (Figure 2 and Supplementary Table S2).

The MLVA genotypes deposited in *Brucella* MLVAbank version 4.3 can be used to compare the incidence of these genotypes in different regions of the world. Genotypes 116 and 111 represent 50 and 3.2%, respectively of the more than 1500 East Mediterranean strains present in MLVAbank (Supplementary Table S1). In Portugal and China, genotype 116 accounts for more than 77% of cases, 37% in Spain and the

<figure><img src="image_3.png" /><figcaption>**FIGURE 2** | Map of Kazakhstan showing the distribution of *B. melitensis* MLVA11 genotypes of the 1327 strains from human patients. The color code reflects the MLVA11 genotype. The number of strains assigned to a given genotype is indicated in detail in **Supplementary Table S2**.</figcaption></figure>

<footnote>*MLVA8 in bold.†New genotype.</footnote>
<figure><img src="image_4.png" /></figure>

prevalence descends further in Turkey (16%) and France (10%) (Figure 3). Genotype 116 is predominant in Asian countries, as is evidenced by the identification of only this genotype in *B. melitensis* isolates from India and Mongolia, although the sampled populations were small. Genotype 111 is also relatively frequent in Asia. The rare genotypes 291 and 387 have been found only in Kazakhstan and China, whereas genotype 114 was found only in Kazakhstan and Turkey. Genotype 125, which is the second most frequent (18.5%) in East Mediterranean strains represented in MLVAbank and is the most prevalent genotype in Greece and Turkey, was not detected in Kazakhstan.

The MLVA16 clustered the 1327 strains into 310 genotypes, 131 of which were represented by a single isolate, 53 by two isolates, and 33 by three isolates. The most prevalent genotype was found in 66 strains (Figure 4). Among the 29 most frequent genotypes (found in 10 or more isolates) one genotype with 17 isolates was restricted to one geographic region (West Kazakhstan), whereas other most prevalent genotypes were present in a wider area. Among 61 medium-prevalence genotypes (e.g., found in four to nine samples), geographic linkage was established for ten genotypes: four genotypes are confined to Almaty region, five genotypes to East Kazakhstan, and one to West Kazakhstan. Figure 5 shows the diversity of genotypes observed in humans, together with diversity previously reported in animals (Shevtsov et al., 2015).

Clustering analysis using all 16 VNTR loci (MLVA16) was conducted including *B. melitensis* East Mediterranean data

from MLVAbank. Shared genotypes were found for 83 (among 313) genotypes representing more than 57% of the isolates (Figure 6). The highest number of shared genotypes was found for strains circulating in China and Kazakhstan. Fifty-four MLVA16 genotypes are common between 530 isolates from Kazakhstan and 172 strains from China. Shared genotypes were also observed with *Brucella* strains isolated in other Asian and European countries. Eleven MLVA16 genotypes were present in a large geographic space, including Kazakhstan (242 isolates), China and other Asian and Middle-East countries (Turkey, Saudi Arabia, India, Mongolia, and Kuwait). Identical MLVA16 genotypes were observed among strains from Kazakhstan and Syria (one genotype, one strain for each country), Kazakhstan and Afghanistan (1 genotype, 10 and 1 strain(s), respectively), Kazakhstan and India (one genotype, one strain from each country), Kazakhstan and Saudi Arabia (1 genotype, 10 and 1 strain(s), respectively), Kazakhstan and Turkey (5 genotypes, 26 and 7 strains, respectively). Ten MLVA16 genotypes from Kazakhstan (comprising 37 strains) were found in European countries: Spain (6 genotypes, 11 strains), Portugal (4 genotypes, 6 strains), and France and Belgium (2 genotypes, 2 strains).

DISCUSSION

We report here the largest investigation of the genetic diversity of *Brucella* in a single country. More than 1300 strains were
<figure><img src="image_5.png" /></figure>

**FIGURE 4** | Maximum parsimony analysis using MLVA16 for 1152 genotype 116 *B. melitensis* strains isolated from humans in Kazakhstan. The color code reflect the geographic origin of the strains (South: Almaty, Qyzylorda, Zhambyl, South-Kazakhstan; West: Aqtobe, Atyrau, West Kazakhstan; East: East Kazakhstan; North: Akmola, Kostanay, Karaganda, Pavlodar, North Kazakhstan). Inset: color code by region for a subset of the most frequent MLVA16 genotypes showing their wide geographic spreading. Branch lengths correspond to one up to two differences among the 16 values constituting the MLVA16 genotype.
<figure><img src="image_6.png" /><figcaption>**FIGURE 5** | Maximum parsimony analysis using MLVA16 data from 1455 *B. melitensis* strains isolated from humans and animals in Kazakhstan. The color code reflects the host origin (cattle: 1 strain, genotype 116; dog: 5 strains, 2 belong to genotype 111, 3 to genotype 116; human: 1327 strains; sheep: 122 strains).</figcaption></figure>

collected from human patients and genotyped using tandem repeats polymorphism. Despite the fact that both *B. melitensis* and *B. abortus* are circulating among farm animals, and that the live vaccines from strains Rev1 (*B. melitensis* "Americas"), S19, SR82, and RB51 (*B. abortus*) are used in Kazakhstan, all strains from humans in this study were identified as *B. melitensis* "East-Mediterranean" (Shevtsova et al., 2016; Daugaliyeva et al., 2018). *B. melitensis* biovar 3 is predominant in this collection (79.8%), whereas biovars 1 and 2 account for 19 and 1.2%, respectively. Predominance of biovar 3 is typical in neighboring countries as well as in the whole Eurasian region (FAO, 2010; Liu et al., 2017). The results obtained show a very good correlation between the genotypes present in livestock in Kazakhstan and the genotypes observed in humans, as expected if infections occurred locally. Brucellosis is a typical zoonotic disease and better monitoring of *Brucella* circulation in animal populations will hopefully contribute to the successful control of human brucellosis.

Multiple-locus VNTR analysis genotyping assigned all strains to the *B. melitensis* "East-Mediterranean" group with predominance of the MLVA11 genotypes 116 and 111 (86.8 and 6.03%, respectively). The recovery of genotype 116 across most of Kazakhstan may be explained by extensive movement of livestock between regions by nomadic herders. The nomadic style is a historic way of livestock husbandry in Kazakhstan and was preserved until the beginning of the 20th century. Nomadism

implies regular movement of herds to new pastures many times per year (Mukhatova, 2014). Similarly, uncontrolled movement of livestock is described as the main source for geographic propagation of the infection (Syzdykov et al., 2014).

In contrast, genotype 111 is almost exclusively observed in South Kazakhstan, Almaty, Zhambyl, East-Kazakhstan, and Qyzylorda regions. These regions suffer from the most aggravated epidemiological situation. They account for 70% of infected sheep in Kazakhstan (during 5-year observation, 2011–2015), and hold 67% of small ruminants population in the country (Syzdykov et al., 2014; Espembetov et al., 2017). Two of the mentioned Kazakhstan’s regions share long borders with China. Information is lacking on genetic diversity in the other countries bordering Kazakhstan along the southwestern border.

The distribution of MLVA genotypes in Kazakhstan circulating in animals and humans closely replicated the genotypic distribution in the brucellosis-hyperendemic regions of Northern China, and more than 50% of the isolates from Kazakhstan have analogs among strains from China. This similarity supposedly reflects a common circulation in the neighboring countries, although it is not clear if the cross-border circulation still exist. MLVA typing does not allow establishing at what historical epoch the genotypic structure shaped into its current form. Trade ties and livestock exchange between the countries have long history, archeological finds
<figure><img src="image_7.png" /><figcaption>**FIGURE 6** | Geographic origin of strains present in MLVAbank *Brucella* 4_3 and showing an MLVA16 genotype present in Kazakhstan.</figcaption></figure>

indicate trade relations between the nomadic peoples who inhabited Central Asia and China long before the Common Era (CE) (Dicosmo, 1996). In addition, 4500 years ago, seasonal nomadic pastoralist routes were formed from modern southern Kazakhstan to Xinjiang Uygur Autonomous Region, covering the territory of more than 70% of the high-mountain route of the silk road (Frachetti et al., 2017). Cross-border nomadic pastoralism played an important role in the functioning of the Silk Road until the 15th century CE (Hermes et al., 2018). Trade relations between China and nomadic pastoralists of Kazakhstan continued until the formation of the Soviet Union. During the Soviet period, cattle was run from Mongolia to the Western parts of the Soviet Union passing through Kazakhstan (Kuznetsov, 1962). In modern history there is no trade in live

cattle between Kazakhstan and China, but actively developing trade in livestock products.

## CONCLUSION

Multiple-locus VNTR analysis genotyping constitute a convenient medium-resolution classification assay for large-scale investigations. Follow-up studies based on whole genome sequencing of selected strains will be necessary to more precisely decipher the dynamics of strain circulation both within Kazakhstan and with Kazakhstan's neighboring countries. Such studies will also permit to design new low-cost genotyping assays tailored for Central Asia. These assays may be based on
a combination of selected tandem repeats including currently
unused loci and of key SNPs.

DATA AVAILABILITY

The datasets generated for this study can be found in http://
microbesgenotyping.i2bc.paris-saclay.fr.

AUTHOR CONTRIBUTIONS

GV edited the original manuscript. AShe conceptualized and
designed the experiments. AShe, ES, and GV analyzed the data.
ES, AShe, AShu, KB, KM, MF, and YR were involved in the
genotyping and wrote the manuscript. MS, AK, LL, and UI did

bacteriological researches and wrote the manuscript. All authors
have read and approved the final version of the manuscript.

FUNDING

This study was supported by the Ministry of Education and
Science of the Kazakhstan (Grant No. AP05133053).

SUPPLEMENTARY MATERIAL

The Supplementary Material for this article can be found
online at: https://www.frontiersin.org/articles/10.3389/fmicb.
2019.01897/full#supplementary-material

REFERENCES

Al Dahouk, S., Le Flèche, P., Nockler, K., Jacques, I., Grayon, M., Scholz, H. C., et al. (2007). Evaluation of *Brucella* MLVA typing for human brucellosis. *J. Microbiol. Methods* 69, 137–145. doi: 10.1016/j.mimet.2006.12.015

Al Dahouk, S., Tomaso, H., Nockler, K., Neubauer, H., and Frangoulidis, D. (2003). Laboratory-based diagnosis of brucellosis—a review of the literature. Part I: techniques for direct detection and identification of *Brucella* spp. *Clin. Lab.* 49, 487–505.

Castaneda, M. R. (1947). A practical method for routine blood cultures in brucellosis. *Proc. Soc. Exp. Biol. Med.* 64:114. doi: 10.3181/00379727-64-15717

Charypkhan, D., Sultanov, A. A., Ivanov, N. P., Baramova, S. A., Taitubayev, M. K., and Torgerson, P. R. (2019). Economic and health burden of brucellosis in Kazakhstan. *Zoonoses Public Health* 66, 487–494. doi: 10.1111/zph.12582

Daugaliyeva, A., Sultanov, A., Usserbayev, B., Baramova, S., Modesto, P., Adambayeva, A., et al. (2018). Genotyping of *Brucella melitensis* and *Brucella abortus* strains in Kazakhstan using MLVA-15. *Infect. Genet. Evol.* 58, 135–144. doi: 10.1016/j.meegid.2017.12.022

Dicosmo, N. (1996). Ancient xinjiang between central asia and china - the nomadic factor. *Anthropol. Archeol. Eurasia* 34, 87–101

Doganay, G. D., and Doganay, M. (2013). *Brucella* as a potential agent of bioterrorism. *Recent Pat. Antiinfect. Drug Discov.* 8, 27–33.

Espembetov, B. A., Syrvm, N. S., and Zinina, N. N. (2017). Monitoring and analysis of the epizootic situation of animal brucellosis in Kazakhstan in 2011-2015. *Bull. Ulyanovsk State Agric. Acad.* 37, 92-96. doi: 10.18286/1816-4501-2017-1-92-96

FAO (2010). *Brucella Melitensis in Eurasia and the Middle East*. Rome: FAO Animal Production and Health Proceedings.

Ferreira, A. C., Chambel, L., Tenreiro, T., Cardoso, R., Flor, L., Dias, I. T., et al. (2012). MLVA16 typing of portuguese human and animal *Brucella melitensis* and *Brucella abortus* isolates. *PLoS One* 7:e42514. doi: 10.1371/journal.pone.0042514

Ficht, T. (2010). *Brucella taxonomy* and evolution. *Future Microbiol.* 5, 859–866. doi: 10.2217/fmb.10.52

Frachetti, M. D., Smith, C. E., Traub, C. M., and Williams, T. (2017). Nomadic ecology shaped the highland geography of Asia's silk roads. *Nature* 543, 193–198. doi: 10.1038/nature21696

Garofolo, G., Ancora, M., and Di Giannatale, E. (2013). MLVA-16 loci panel on *Brucella* spp. using multiplex PCR and multicolor capillary electrophoresis. *J. Microbiol. Methods* 92, 103–107. doi: 10.1016/j.mimet.2012.11.007

Godfroid, J., Scholz, H. C., Barbier, T., Nicolas, C., Wattiau, P., Fretin, D., et al. (2011). Brucellosis at the animal/ecosystem/human interface at the beginning of the 21st century. *Prev. Vet. Med.* 102, 118–131. doi: 10.1016/j.prevetmed.2011.04.007

Gwida, M., Neubauer, H., Ilhan, Z., Schmoock, G., Melzer, F., Nockler, K., et al. (2012). Cross-border molecular tracing of brucellosis in Europe. *Comp. Immunol. Microbiol. Infect. Dis.* 35, 181–185. doi: 10.1016/j.cimid.2011.12.012

Hermes, T. R., Frachetti, M. D., Bullion, E. A., Maksudov, F., Mustafokulov, S., and Makarewicz, C. A. (2018). Urban and nomadic isotopic niches reveal dietary connectivities along central Asia's silk roads. *Sci. Rep.* 8:5177. doi: 10.1038/s41598-018-22995-22992

Hunter, P. R., and Gaston, M. A. (1988). Numerical index of the discriminatory ability of typing systems: an application of Simpson's index of diversity. *J. Clin. Microbiol.* 26, 2465-2466.

Kattar, M. M., Jaafar, R. F., Araj, G. F., Le Flèche, P., Matar, G. M., Abi Rached, R., et al. (2008). Evaluation of a multilocus variable-number tandem-repeat analysis scheme for typing human *Brucella* isolates in a region of brucellosis endemicity. *J. Clin. Microbiol.* 46, 3935-3940. doi: 10.1128/JCM.00464-468

Kilic, S., Ivanov, I. N., Durmaz, R., Bayraktar, M. R., Ayaslioglu, E., Uyanik, M. H., et al. (2011). Multiple-locus variable-number tandem-repeat analysis genotyping of human *Brucella* isolates from Turkey. *J. Clin. Microbiol.* 49, 3276-3283. doi: 10.1128/JCM.02538-2510

Kuznetsov, V. S. (1962). *Kazakhsko-Kitayskiye Torgovyye Otnosheniya v Kontse XVIII veke [Kazakh-Chinese trade relations at the end of the XVIII century]*. Almaty: Publishing house of the Academy of Sciences of the Kazakh SSR.

Le Flèche, P., Jacques, I., Grayon, M., Al Dahouk, S., Bouchon, P., Denoeud, F., et al. (2006). Evaluation and selection of tandem repeat loci for a *Brucella* MLVA typing assay. *BMC Microbiol.* 6:9. doi: 10.1186/1471-2180-6-9

Ledwaba, M. B., Gomo, C., Lekota, K. E., Le Flèche, P., Hassim, A., Vergnaud, G., et al. (2019). Molecular characterization of *Brucella* species from Zimbabwe. *PLoS Negl. Trop Dis.* 13:e0007311. doi: 10.1371/journal.pntd.0007311

Liu, Z. G., Di, D. D., Wang, M., Liu, R. H., Zhao, H. Y., Piao, D. R., et al. (2017). MLVA genotyping characteristics of human *Brucella melitensis* isolated from ulanqab of inner mongolia. China. *Front. Microbiol.* 8:6. doi: 10.3389/fmicb.2017.00006

Lopez-Goni, I., Garcia-Yoldi, D., Marin, C. M., de Miguel, M. J., Munoz, P. M., Blasco, J. M., et al. (2008). Evaluation of a multiplex PCR assay (Bruce-ladder) for molecular typing of all *Brucella* species, including the vaccine strains. *J. Clin. Microbiol.* 46, 3484–3487. doi: 10.1128/JCM.00837-838

Maquart, M., Le Flèche, P., Foster, G., Tryland, M., Ramisse, F., Djonne, B., et al. (2009). MLVA-16 typing of 295 marine mammal *Brucella* isolates from different animal and geographic origins identifies 7 major groups within *Brucella ceti* and *Brucella pinnipedialis*. *BMC Microbiol.* 9:145. doi: 10.1186/1471-2180-9-145

Meyer, M. E. (1990). "Evolutionary development and taxonomy of the genus *Brucella*," in *Advances in Brucellosis Research*, ed. L. G. Adams (Galveston, TX: Texas A&M University Press), 12-35.

Mizanbayeva, S., Smits, H. L., Zhalilova, K., Abdoel, T. H., Kozakov, S., Ospanov, K. S., et al. (2009). The evaluation of a user-friendly lateral flow assay for the serodiagnosis of human brucellosis in Kazakhstan. *Diagn. Microbiol. Infect. Dis.* 65, 14–20. doi: 10.1016/j.diagmicrobio.2009.05.002

Mukhatova, O. (2014). Kazakh economy In XVii-Xviii centuries. 3rd world conference on educational technology researches 2013. *Wcetr* 2013, 218–223. doi: 10.1016/j.sbspro.2014.04.107

Neubauer, H. (2010). Brucellosis: new demands in a changing world. *Prilozi* 31, 209–217.
Nicoletti, P. (2010). Brucellosis: past, present and future. *Prilozi* 31, 21–32.

Pappas, G., Papadimitriou, P., Akritidis, N., Christou, L., and Tsianos, E. V. (2006a). The new global map of human brucellosis. *Lancet Infect. Dis.* 6, 91–99. doi: 10.1016/S1473-3099(06)70382-70386

Pappas, G., Panagopoulou, P., Christou, L., and Akritidis, N. (2006b). *Brucella* as a biological weapon. *Cell. Mol. Life Sci.* 63, 2229–2236. doi: 10.1007/s00018-006-6311-6314

Scholz, H. C., and Vergnaud, G. (2013). Molecular characterisation of *Brucella* species. *Rev. Sci. Tech.* 32, 149–162. doi: 10.20506/rst.32.1.2189

Shevtsov, A., Ramanculov, E., Shevtsova, E., Kairzhanova, A., Tarlykov, P., Filipenko, M., et al. (2015). Genetic diversity of *Brucella abortus* and *Brucella melitensis* in Kazakhstan using MLVA-16. *Infect. Genet. Evol.* 34, 173–180. doi: 10.1016/j.meegid.2015.07.008

Shevtsova, E., Shevtsov, A., Mukanov, K., Filipenko, M., Kamalova, D., Sytnik, I., et al. (2016). Epidemiology of brucellosis and genetic diversity of *Brucella abortus* in Kazakhstan. *PLoS One* 11:e0167496. doi: 10.1371/journal.pone.0167496

Studentsov, Ë. P. (1975). *Brutsellez zhivotnykh [Brucellosis in animals]*. Almaty: Kaynar.

Syzdykov, M. S., Kuznetsov, A. N., Kazakov, S. V., Daulbayeva, S. F., Duysenova, A. K., Berezovskiy, D. V., et al. (2014). Analiz prostranstvenno-vremennogo raspredeleniya brutselloza cheloveka i zhivotnykh s ispol'zovaniyem Geograficheskikh informatsionnykh tekhnologiy—[Analysis of the spatial and temporal distribution of human and animal brucellosis using Geographic information technology]. *Gigiyena, epidemiologiya i immunobiologiya - [Hygiene, Epidemiology and Immunobiology]* 72, 24–26.

Syzdykov, M. S., Kuznetsov, A. N., Kazakov, S. V., Espembetov, B. A., and Duysenova, A. K. (2013). Spatio-temporal analysis of the incidence of brucellosis using GIS technology [in Rus.]. *Med. Kyrgyzstan* 4, 148–151.

Vergnaud, G., Hauck, Y., Christiany, D., Daoud, B., Pourcel, C., Jacques, I., et al. (2018). Genotypic expansion within the population structure of classical *Brucella* species revealed by MLVA16 Typing of 1404 *Brucella* isolates from

different animal and geographic origins, 1974-2006. *Front. Microbiol.* 9:1545. doi: 10.3389/fmicb.2018.01545

Whatmore, A. M., Koylass, M. S., Muchowski, J., Edwards-Smallbone, J., Gopaul, K. K., and Perrett, L. L. (2016). Extended multilocus sequence analysis to describe the global population structure of the genus *Brucella*: phylogeography and relationship to biovars. *Front. Microbiol.* 7:2049. doi: 10.3389/fmicb.2016.02049

Whatmore, A. M., Shankster, S. J., Perrett, L. L., Murphy, T. J., Brew, S. D., Thirlwall, R. E., et al. (2006). Identification and characterization of variable-number tandem-repeat markers for typing of *Brucella* spp. *J. Clin. Microbiol.* 44, 1982–1993. doi: 10.1128/JCM.02039-2035

Yagupsky, P. (1999). Detection of *Brucellae* in blood cultures. *J. Clin. Microbiol.* 37, 3437–3442.

Yespembetov, B. A., Syrym, N. S., Zinina, N. N., Sarmykova, M. K., Konbayeva, G. M., Basybekov, S. Z., et al. (2019). Phenotypic and genotypic characteristics of *Brucella* isolates from the republic of kazakhstan. *Trop Anim. Health Prod.* doi: 10.1007/s11250-019-01941-y [Epub ahead of print].

Zheludkov, M. M., and Tsirelson, L. E. (2010). Reservoirs of *Brucella* infection in nature. *Biol. Bull.* 37, 709–715. doi: 10.1134/s106235901007006x

**Conflict of Interest Statement:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Copyright © 2019 Shevtsova, Vergnaud, Shevtsov, Shustov, Berdimuratova, Mukanov, Syzdykov, Kuznetsov, Lukhnova, Izbanova, Filipenko and Ramankulov. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.