<?xml version="1.0" encoding="utf-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" article-type="research-article" dtd-version="2.3" xml:lang="EN">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Microbiol.</journal-id>
<journal-title>Frontiers in Microbiology</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Microbiol.</abbrev-journal-title>
<issn pub-type="epub">1664-302X</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.3389/fmicb.2024.1320845</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Microbiology</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Epidemiological changes and molecular characteristics of <italic>Brucella</italic> strains in Ningxia, China</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author" equal-contrib="yes">
<name>
<surname>Liu</surname>
<given-names>Guangtian</given-names>
</name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="author-notes" rid="fn0001"><sup>&#x2020;</sup></xref>
<uri xlink:href="https://loop.frontiersin.org/people/1329450/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/funding-acquisition/"/>
<role content-type="https://credit.niso.org/contributor-roles/investigation/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
<contrib contrib-type="author" equal-contrib="yes">
<name>
<surname>Ma</surname>
<given-names>Xueping</given-names>
</name>
<xref ref-type="aff" rid="aff2"><sup>2</sup></xref>
<xref ref-type="author-notes" rid="fn0001"><sup>&#x2020;</sup></xref>
<uri xlink:href="https://loop.frontiersin.org/people/2581279/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/investigation/"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
<role content-type="https://credit.niso.org/contributor-roles/resources/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Zhang</surname>
<given-names>Ruiqing</given-names>
</name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
<role content-type="https://credit.niso.org/contributor-roles/formal-analysis/"/>
<role content-type="https://credit.niso.org/contributor-roles/project-administration/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>L&#x00FC;</surname>
<given-names>Jufen</given-names>
</name>
<xref ref-type="aff" rid="aff3"><sup>3</sup></xref>
<uri xlink:href="https://loop.frontiersin.org/people/2545330/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/supervision/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Zhou</surname>
<given-names>Pan</given-names>
</name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
<role content-type="https://credit.niso.org/contributor-roles/data-curation/"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Liu</surname>
<given-names>Bofei</given-names>
</name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<role content-type="https://credit.niso.org/contributor-roles/project-administration/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
<role content-type="https://credit.niso.org/contributor-roles/resources/"/>
<role content-type="https://credit.niso.org/contributor-roles/supervision/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Liu</surname>
<given-names>Tao</given-names>
</name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<role content-type="https://credit.niso.org/contributor-roles/data-curation/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
<role content-type="https://credit.niso.org/contributor-roles/software/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Ren</surname>
<given-names>Hui</given-names>
</name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<role content-type="https://credit.niso.org/contributor-roles/supervision/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
<role content-type="https://credit.niso.org/contributor-roles/data-curation/"/>
<role content-type="https://credit.niso.org/contributor-roles/formal-analysis/"/>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Liu</surname>
<given-names>Zhiguo</given-names>
</name>
<xref ref-type="aff" rid="aff4"><sup>4</sup></xref>
<xref ref-type="corresp" rid="c001"><sup>&#x002A;</sup></xref>
<uri xlink:href="https://loop.frontiersin.org/people/369611/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/data-curation/"/>
<role content-type="https://credit.niso.org/contributor-roles/software/"/>
<role content-type="https://credit.niso.org/contributor-roles/conceptualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Li</surname>
<given-names>Zhenjun</given-names>
</name>
<xref ref-type="aff" rid="aff4"><sup>4</sup></xref>
<xref ref-type="corresp" rid="c002"><sup>&#x002A;</sup></xref>
<uri xlink:href="https://loop.frontiersin.org/people/671873/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/supervision/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
<role content-type="https://credit.niso.org/contributor-roles/formal-analysis/"/>
<role content-type="https://credit.niso.org/contributor-roles/project-administration/"/>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Jiang</surname>
<given-names>Xuefeng</given-names>
</name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="corresp" rid="c003"><sup>&#x002A;</sup></xref>
<role content-type="https://credit.niso.org/contributor-roles/project-administration/"/>
<role content-type="https://credit.niso.org/contributor-roles/supervision/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-review-editing/"/>
</contrib>
</contrib-group>
<aff id="aff1"><sup>1</sup><institution>The Fourth People&#x2019;s Hospital of Ningxia Hui Autonomous Region</institution>, <addr-line>Yinchuan</addr-line>, <country>China</country></aff>
<aff id="aff2"><sup>2</sup><institution>Ningxia Hui Autonomous Region Center for Disease Control and Prevention</institution>, <addr-line>Yinchuan</addr-line>, <country>China</country></aff>
<aff id="aff3"><sup>3</sup><institution>The College of Public Health, Ningxia Medical University</institution>, <addr-line>Yinchuan</addr-line>, <country>China</country></aff>
<aff id="aff4"><sup>4</sup><institution>National Key Laboratory of Intelligent Tracking and Forecasting for Infectious Diseases, National Institute for Communicable Disease Control and Prevention, Chinese Center for Disease Control and Prevention</institution>, <addr-line>Beijing</addr-line>, <country>China</country></aff>
<author-notes>
<fn id="fn0002" fn-type="edited-by"><p>Edited by: Georgia Vrioni, National and Kapodistrian University of Athens, Greece</p></fn>
<fn id="fn0003" fn-type="edited-by"><p>Reviewed by: Maryam Dadar, Razi Vaccine and Serum Research Institute, Iran</p>
<p>Sanja Duvnjak, Croatian Veterinary Institute, Croatia</p></fn>
<corresp id="c001">&#x002A;Correspondence: Zhiguo Liu, <email>wlcblzg@126.com</email></corresp>
<corresp id="c002">Zhenjun Li, <email>lizhenjun@icdc.cn</email></corresp>
<corresp id="c003">Xuefeng Jiang, <email>jiang_xf2008@126.com</email></corresp>
<fn id="fn0001" fn-type="equal"><p><sup>&#x2020;</sup>These authors have contributed equally to this work</p></fn>
</author-notes>
<pub-date pub-type="epub">
<day>19</day>
<month>01</month>
<year>2024</year>
</pub-date>
<pub-date pub-type="collection">
<year>2024</year>
</pub-date>
<volume>15</volume>
<elocation-id>1320845</elocation-id>
<history>
<date date-type="received">
<day>13</day>
<month>10</month>
<year>2023</year>
</date>
<date date-type="accepted">
<day>03</day>
<month>01</month>
<year>2024</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#x00A9; 2024 Liu, Ma, Zhang, L&#x00FC;, Zhou, Liu, Liu, Ren, Liu, Li and Jiang.</copyright-statement>
<copyright-year>2024</copyright-year>
<copyright-holder>Liu, Ma, Zhang, L&#x00FC;, Zhou, Liu, Liu, Ren, Liu, Li and Jiang</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/">
<p>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.</p>
</license>
</permissions>
<abstract>
<sec id="sec1">
<title>Objective</title>
<p>Human brucellosis causes serious public health concerns in Ningxia, China.</p>
</sec>
<sec id="sec2">
<title>Methods</title>
<p>This study employed epidemiological, bacteriological, and multiple-locus variable-number tandem repeat analysis (MLVA) methods to conduct an epidemiological investigation, which is necessary for devising tailored control strategies.</p>
</sec>
<sec id="sec3">
<title>Results</title>
<p>Between 1958 and 2022, 29,892 cases were reported, with an average annual number of cases and incidence of 467 and 7.1/100,000, respectively. The epidemic situation gradually worsened, with cases escalating from 26 cases in 2005 to 6,292 in 2022, with the incidence rate rising from 0.441 in 2005 to 86.83 in 2022. Geographically, the disease spread from a single affected county in 2004 to encompass all 22 counties in 2022. Yanchi County had the highest incidence, followed by the Hongsibao and Tongxin counties. These data suggest that <italic>Brucella</italic> infection has become a rampant regional concern in human brucellosis. Between 1958 and 2019, a total of 230 <italic>Brucella</italic> strains were identified across four studied hosts. These strains comprised four species with 12 biovars, including <italic>B. melitensis</italic> bv. 1, bv. 2, bv. 3, <italic>B. abortus</italic> bv. 1, bv. 3, bv. 4, bv. 5, bv. 6, bv. 7, <italic>B. suis</italic> bv. 1 and bv. 3, and <italic>B. canis</italic>. These data highlight the high species/biovars and host diversity of the <italic>Brucella</italic> population, posing a substantial challenge to brucellosis surveillance. There was an apparent transition from multiple species/biovars historically to the current dominance of a single species, <italic>B. melitensis</italic>, emphasizing the requirement for strengthening surveillance of <italic>B. melitensis</italic>. Genotypes 42 and 116, constituting 96.2% of the total number of genotypes, predominated in panel 1 and <italic>MLVA-11</italic>, indicating that all strains belong to the East Mediterranean lineage. <italic>MLVA</italic> cluster analysis revealed persistent transmission of dominant circulating genotypes, presenting an epidemic pattern characterized primarily by epidemiologically related cases with a few sporadic cases. Strains in this study exhibited high genetic homogeneity with strains from the Northwest, and those from Kazakhstan and Mongolia.</p>
</sec>
<sec id="sec4">
<title>Conclusion</title>
<p>The epidemic situation of human brucellosis has gradually worsened; the rampant epidemic of the disease has become a regional concern. The present study highlights that implementing the of targeted surveillance and intervention strategies is urge.</p>
</sec>
</abstract>
<kwd-group>
<kwd>human brucellosis</kwd>
<kwd>epidemiology</kwd>
<kwd>bacteriology</kwd>
<kwd>MLVA</kwd>
<kwd>molecular</kwd>
</kwd-group>
<counts>
<fig-count count="4"/>
<table-count count="2"/>
<equation-count count="0"/>
<ref-count count="36"/>
<page-count count="9"/>
<word-count count="6346"/>
</counts>
<custom-meta-wrap>
<custom-meta>
<meta-name>section-at-acceptance</meta-name>
<meta-value>Infectious Agents and Disease</meta-value>
</custom-meta>
</custom-meta-wrap>
</article-meta>
</front>
<body>
<sec sec-type="intro" id="sec5">
<title>Introduction</title>
<p>Human brucellosis, caused by <italic>Brucella</italic> spp., is a global public health concern (<xref ref-type="bibr" rid="ref24">Pappas et al., 2006</xref>; <xref ref-type="bibr" rid="ref11">Di Bonaventura et al., 2021</xref>). This disease exerts negative effects on human health and the animal farming industry. <italic>Brucella</italic> strains can result in severe clinical manifestations in the population, including fever, sweat, and fatigue, along with substantial economic losses in animal farming due to factors such as abortion and decreased milk productivity (<xref ref-type="bibr" rid="ref6">Bagheri Nejad et al., 2020</xref>; <xref ref-type="bibr" rid="ref16">Jamil et al., 2021</xref>). Although some developed countries have nearly eliminated this disease, the majority of low-income countries remain vulnerable to this zoonotic threat, including Africa, Latin America, Asia, and China (<xref ref-type="bibr" rid="ref14">Gwida et al., 2010</xref>; <xref ref-type="bibr" rid="ref31">Zhu et al., 2020</xref>). Since the 2000s, the incidence of human brucellosis has been continuously increasing and is distributed across all mainland regions in China, including the Ningxia Hui Autonomous Region (Ningxia) (<xref ref-type="bibr" rid="ref19">Lai et al., 2017</xref>). Ningxia has a history of being an endemic area for brucellosis. The first survey of human brucellosis was conducted from 1958 to 1962, involving the sampling of 18,698 people, among whom 3,185 tested positive, resulting in a positivity rate of 17.03%. At that time, immunization had not yet been widely promoted, leading to the infection rate gradually increasing. Subsequently, a second brucellosis survey conducted from 1974 to 1980 revealed a positivity rate of 4.34% (68,609/1580,082). Following the implementation of mass ruminant immunization, the human brucellosis epidemic has been relatively stable, with 76.19% of the counties in the region reaching the national standard for stable brucellosis control (<xref ref-type="bibr" rid="ref10">Deqiu et al., 2002</xref>).</p>
<p>Although vaccination and test-slaughter program animals, and surveillance and intervention plan in human population were exerted continuously. Since 2008, human brucellosis has steadily increased, with the current incidence surpassing historical levels. Conducting timely epidemiological analysis of resurging brucellosis areas is vital to enhance disease surveillance and improve intervention and treatment outcomes (<xref ref-type="bibr" rid="ref12">Franco et al., 2007</xref>; <xref ref-type="bibr" rid="ref27">Yagupsky et al., 2019</xref>; <xref ref-type="bibr" rid="ref22">Matle et al., 2021</xref>). Additionally, pathogen surveillance is necessary for the formulation of targeted control measures and to provide a source for molecular epidemiological investigations (<xref ref-type="bibr" rid="ref20">Laine et al., 2022</xref>). Therefore, the aim of this study was to illustrate the epidemiology of human brucellosis and the molecular association of <italic>Brucella</italic> strains in Ningxia to provide valuable insights for the development of targeted measures to halt its spread. To achieve this, we employed epidemiological analysis, pathogen surveillance, and molecular epidemiological approaches, including multiple-locus variable-number tandem repeat analysis (MLVA), which has become a routine epidemiological tool for distinguishing isolates from specific geographic origins and conducting trace-back analyses (<xref ref-type="bibr" rid="ref17">Kili&#x00E7; et al., 2011</xref>; <xref ref-type="bibr" rid="ref9">De Massis et al., 2015</xref>).</p>
</sec>
<sec sec-type="methods" id="sec6">
<title>Methods</title>
<sec id="sec7">
<title>Epidemiological data collection and processing, spatial cluster analysis</title>
<p>Data including the number of human brucellosis cases and incidence rates in Ningxia from 2003 to 2022 were obtained from the China Information System for Disease Control and Prevention, while epidemic data from 1958 to 2002 were acquired from the Annals of Health in China. Epidemic indices, such as cases and incidence rates, were utilized to depict brucellosis epidemic characteristics. The distribution maps of the annual incidence rate of human brucellosis from 1958 to 2022 were visualized using ArcGIS version 10.7 software (Esri; Redlands, California, United States), based on city boundaries. Subsequently, global spatial autocorrelation analysis was performed using the spatial autocorrelation tool in ArcGIS software (Version 10.6, Environmental Systems Research Institute, Redlands, United States). The global Moran&#x2019;s <italic>I</italic> value served as a measure, where Moran&#x2019;s <italic>I</italic>&#x2009;=&#x2009;0 indicates that the regionalized variables are randomly distributed in space, Moran&#x2019;s <italic>I</italic>&#x2009;&#x003E;&#x2009;0 suggests a positive correlation, and Moran&#x2019;s <italic>I</italic>&#x2009;&#x003C;&#x2009;0 implies a negative correlation. Local autocorrelation is present only when |<italic>z</italic>|&#x2009;&#x2265;&#x2009;1.96 and <italic>p</italic>&#x2009;&#x2264;&#x2009;0.05 (<italic>&#x03B1;</italic>&#x2009;=&#x2009;0.05). Local Moran&#x2019;s <italic>I</italic> significance maps and cluster maps were utilized to illustrate the four different types of local correlation and the significance of the corresponding Moran&#x2019;s <italic>I</italic> index. High&#x2013;high regions denote geographical areas with a high incidence of brucellosis, encircled by other areas also exhibiting high incidences. Conversely, low&#x2013;low regions represent geographical locations with low incidences of brucellosis, surrounded by areas with similarly low incidences. High&#x2013;high and low&#x2013;low regions signify clusters of geographic regions with analogous brucellosis incidence values, whereas high&#x2013;low and low&#x2013;high regions indicate spatial outliers (<xref ref-type="bibr" rid="ref9005">Luc, 1995</xref>).</p>
</sec>
<sec id="sec8">
<title>Sample collection, bacteriological isolation, and bio-typing of <italic>Brucella</italic> strains</title>
<p>The sampling data of strains isolated from the 1950s&#x2013;1980s historically collected, sampling information is missing. After 2009, strains were isolated from human blood samples, 570 blood samples of patients were collected, and <italic>Brucella</italic> strains were isolated using a bacteriological approach (<xref ref-type="bibr" rid="ref11">Di Bonaventura et al., 2021</xref>). Briefly, 5&#x2013;10&#x2009;mL of fresh blood was sampled and injected into a biphasic culture flask within a biosafety cabinet, followed by incubation at 37&#x00B0;C for a minimum of 30 continuous days. The flask was gently shaken after observation every other day. The Rose Bengal plate test and serum agglutination test were used for brucellosis screening and diagnosis (<xref ref-type="bibr" rid="ref5">Araj, 2010</xref>). The suspected clones were further incubated, and bio-typing was identified using phenotypic methods previously described (<xref ref-type="bibr" rid="ref27">Yagupsky et al., 2019</xref>), including CO<sub>2</sub> requirement for growth, H<sub>2</sub>S production, growth inhibition by basic fuchsin and thionin, agglutination with monospecific antisera (A and M), and the phage lysis test.</p>
</sec>
<sec id="sec9">
<title>DNA extraction and MLVA genotyping of <italic>Brucella</italic> strains in this study</title>
<p>DNA from 105 (out of 230) available strains was prepared according to the manufacturer&#x2019;s protocol (Qiagen, Heidelberg, Germany). All examined strains were confirmed to be <italic>B. melitensis</italic> through the detection of a 731&#x2009;bp band specific to this species by AMOS-PCR (<xref ref-type="bibr" rid="ref9001">Bricker and Halling, 1994</xref>). Subsequently, these 105 <italic>B. melitensis</italic> strains, comprising 48 strains from 2019, 21 from 2013, 10 from 2012, 3 from 2011, 4 from 2010, 6 from 2009, 12 from 1979, and 1 from 1973 (<xref rid="SM1" ref-type="supplementary-material">Supplementary Table S1</xref>), were subjected to the MLVA approach (<xref ref-type="bibr" rid="ref9003">Le Fl&#x00E8;che et al., 2006</xref>; <xref ref-type="bibr" rid="ref21">Liu et al., 2019</xref>) for cluster analysis. Briefly, the genotyping approach, as previously described (<xref ref-type="bibr" rid="ref9004">Li et al., 2020</xref>), involved using panel 1 (MLVA-8), MLVA-11, and MLVA-16 to characterize the strains.</p>
<p>Polymorphisms at each locus and panels were analyzed using the Hunter-Gaston Diversity index (<xref ref-type="bibr" rid="ref9002">Hunter, 1990</xref>). Phylogenetic analysis was performed using unweighted pair group method using arithmetic averages in Bionumerics version 5.1 (Applied Maths, St-Martens-Latem, Belgium). To further explore the molecular epidemiology of <italic>B. melitensis</italic> strains at national (<italic>n</italic>&#x2009;=&#x2009;595) (<xref rid="SM2" ref-type="supplementary-material">Supplementary Table S2</xref>) from previous study (<xref ref-type="bibr" rid="ref31">Zhu et al., 2020</xref>) and international scales (<italic>n</italic>&#x2009;=&#x2009;750, <xref rid="SM3" ref-type="supplementary-material">Supplementary Table S3</xref>) (MLVA genotyping database V1.4.0, <ext-link xlink:href="https://microbesgenotyping.i2bc.paris-saclay.fr/databases" ext-link-type="uri">https://microbesgenotyping.i2bc.paris-saclay.fr/databases</ext-link>), and genetic relationships between the 124 isolates were evaluated using a minimum spanning tree (MST) constructed by PHYLOViZ 2.0 software (<xref ref-type="bibr" rid="ref23">Nascimento et al., 2016</xref>).</p>
</sec>
</sec>
<sec sec-type="results" id="sec10">
<title>Results</title>
<sec id="sec11">
<title>Epidemiological characteristics of human brucellosis from 1958 to 2022</title>
<p>A total of 29,892 cases were reported between 1958 and 2022, with an average annual number of cases and incidence rates of 467 and 7.1/100,000, respectively. Most cases (6,295) were reported in 2022, with an incidence rate of 86.83/100,000 (<xref ref-type="fig" rid="fig1">Figure 1</xref>). The epidemic phases of human brucellosis in Ningxia can be divided into three phases. The first phase, spanning from 1958 to 1967, was characterized as a low epidemic period, with an incidence rate increasing from 0.05/100,000 in 1959 to 10.87/100,000 in 1965. This was followed by a controlled phase from 1968 to 2004, with an incidence rate consistently below 0.45/100,000. Since 2005, it has entered a re-emerging epidemic phase. After 2004, the number of cases gradually increased annually, rising from 26 cases in 2005 to 6,292 in 2022. The incidence per 100,000 increased from 0.441 in 2005 to 86.83 in 2022.</p>
<fig position="float" id="fig1">
<label>Figure 1</label>
<caption><p>Human brucellosis epidemic profile between 1958 and 2022 in Ningxia, China.</p></caption>
<graphic xlink:href="fmicb-15-1320845-g001.tif"/>
</fig>
</sec>
<sec id="sec12">
<title>Geographic distribution and spatial autocorrelation analysis of human brucellosis</title>
<p>Geographically, the affected area continuously expanded from a single county in 2004 to encompass all 22 counties in 2022 (<xref ref-type="fig" rid="fig2">Figure 2</xref>). The cumulative number of cases across geographical regions ranged from 220 to 3,679, and the incidence rate ranged from 4.22/100,000 to 127.99/100,000 (<xref ref-type="fig" rid="fig2">Figure 2</xref>). Yanchi county reported the highest number of cases (3,679), followed by Tongxi county (2,702) and Hongsibao district (2,679), with average annual incidence rates (/100,000) of 127.99, 43.72, and 77.83, respectively (<xref ref-type="fig" rid="fig2">Figure 2</xref>). Conversely, Longde county had the lowest number of cases (220) and an annual incidence rate of 9.64/100,000. Notably, the greatest disease burden was from Wuzhong city, encompassing Yanchi county, Tongxi county, and Hongsibao district, which share borders with Inner Mongolia, Gansu, Shaanxi, and Shanxi provinces. Additionally, a spatial autocorrelation analysis based on the annual average incidence revealed a significant spatial cluster at the county level, with a Moran Index of 0.534, a <italic>Z</italic> score of 2.56, and a <italic>p</italic>-value &#x2264;0.05 (0.01). Subsequently, a local autocorrelation analysis indicated a high&#x2013;high cluster area in Hongsibao district and found five low&#x2013;low cluster areas, including Pingluo county, Helanshan city, Xiaxi district, Jinfeng district, and Yongning county (<xref rid="SM5" ref-type="supplementary-material">Supplementary Figure S1</xref>).</p>
<fig position="float" id="fig2">
<label>Figure 2</label>
<caption><p>The average annual incidence of brucellosis across 22 counties in Ningxia from 2004 to 2022. Numbers in the figure showed the territory (<italic>n</italic>&#x2009;=&#x2009;230) and MLVA-16 genotypes (<italic>n</italic>&#x2009;=&#x2009;100) distribute profile of <italic>Brucella</italic> strains.</p></caption>
<graphic xlink:href="fmicb-15-1320845-g002.tif"/>
</fig>
</sec>
<sec id="sec13">
<title>Geographic distribution of species/biovars and host profiles of <italic>Brucella</italic> strains</title>
<p>A total of 230 <italic>Brucella</italic> strains were collected, and isolated between 1950 and 2019 (<xref ref-type="table" rid="tab1">Table 1</xref>), out of which 138 strains were collected from 1950s to 1980s, and identified them based on conventional bio-typing approaches. The remainder 92 strains were isolated from 2009 to 2019, and confirmed as <italic>B. melitensis</italic> through detection of a 731&#x2009;bp <italic>B. melitensis</italic> specific band, conventional bio-typing tests used to determinate the biovar of 92 <italic>B. melitensis</italic> strains. The 230 strains were collected from 14 regions, with the majority originating from Lingwu City (<italic>n</italic>&#x2009;=&#x2009;58), followed by 48 from Yanchi, 32 from Wuzhong county, and 24 from Yinchuan. The range of strains in the other 11 areas was 1&#x2013;24 (<xref ref-type="fig" rid="fig2">Figure 2</xref> and <xref ref-type="table" rid="tab1">Table 1</xref>). The 230 <italic>Brucella</italic> strains in the present study exhibited high species/biovar diversity, encompassing four species and 12 biovars (<xref ref-type="table" rid="tab2">Table 2</xref>). Among these, <italic>B. melitensis</italic> strains were the predominant population in this region. Furthermore, <italic>B. melitensis</italic> strains exhibited the highest host diversity, which included humans, cattle, sheep, and dogs. <italic>B. abortus</italic> was found in cattle and sheep, while <italic>B. suis</italic> and <italic>B. canis</italic> were isolated from humans and dogs (<xref ref-type="table" rid="tab2">Table 2</xref>).</p>
<table-wrap position="float" id="tab1">
<label>Table 1</label>
<caption><p>The regional distribution of <italic>Brucella</italic> strains from 1950 to 2019.</p></caption>
<table frame="hsides" rules="groups">
<thead>
<tr>
<th align="left" valign="top">Counties</th>
<th align="center" valign="top">50s</th>
<th align="center" valign="top">60s</th>
<th align="center" valign="top">70s</th>
<th align="center" valign="top">80s</th>
<th align="center" valign="top">2009&#x2013;2013</th>
<th align="center" valign="top">2019</th>
<th align="center" valign="top">Total</th>
</tr>
</thead>
<tbody>
<tr>
<td align="left" valign="top">Xixia</td>
<td align="center" valign="top">1</td>
<td align="center" valign="top">5</td>
<td align="center" valign="top">13</td>
<td align="center" valign="top">4</td>
<td align="center" valign="top">1</td>
<td align="center" valign="top">/</td>
<td align="center" valign="top">24</td>
</tr>
<tr>
<td align="left" valign="top">Helan</td>
<td align="center" valign="top">/</td>
<td align="center" valign="top">9</td>
<td align="center" valign="top">8</td>
<td align="center" valign="top">1</td>
<td align="center" valign="top">/</td>
<td align="center" valign="top">/</td>
<td align="center" valign="top">18</td>
</tr>
<tr>
<td align="left" valign="top">Pingluo</td>
<td align="center" valign="top">/</td>
<td align="center" valign="top">7</td>
<td align="center" valign="top">1</td>
<td align="center" valign="top">4</td>
<td align="center" valign="top">2</td>
<td align="center" valign="top">/</td>
<td align="center" valign="top">14</td>
</tr>
<tr>
<td align="left" valign="top">Dawukou</td>
<td align="center" valign="top">/</td>
<td align="center" valign="top">/</td>
<td align="center" valign="top">4</td>
<td align="center" valign="top">2</td>
<td align="center" valign="top">/</td>
<td align="center" valign="top">/</td>
<td align="center" valign="top">6</td>
</tr>
<tr>
<td align="left" valign="top">Huinong district</td>
<td align="center" valign="top">/</td>
<td align="center" valign="top">5</td>
<td align="center" valign="top">/</td>
<td align="center" valign="top">/</td>
<td align="center" valign="top">/</td>
<td align="center" valign="top">/</td>
<td align="center" valign="top">5</td>
</tr>
<tr>
<td align="left" valign="top">Lingwu</td>
<td align="center" valign="top">/</td>
<td align="center" valign="top">54</td>
<td align="center" valign="top">/</td>
<td align="center" valign="top">/</td>
<td align="center" valign="top">/</td>
<td align="center" valign="top">4</td>
<td align="center" valign="top">58</td>
</tr>
<tr>
<td align="left" valign="top">Zhongning</td>
<td align="center" valign="top">/</td>
<td align="center" valign="top">/</td>
<td align="center" valign="top">/</td>
<td align="center" valign="top">1</td>
<td align="center" valign="top">/</td>
<td align="center" valign="top">/</td>
<td align="center" valign="top">1</td>
</tr>
<tr>
<td align="left" valign="top">Tongxin</td>
<td align="center" valign="top">/</td>
<td align="center" valign="top">/</td>
<td align="center" valign="top">1</td>
<td align="center" valign="top">/</td>
<td align="center" valign="top">/</td>
<td align="center" valign="top">/</td>
<td align="center" valign="top">1</td>
</tr>
<tr>
<td align="left" valign="top">Shapotou district</td>
<td align="center" valign="top">/</td>
<td align="center" valign="top">/</td>
<td align="center" valign="top">/</td>
<td align="center" valign="top">6</td>
<td align="center" valign="top">/</td>
<td align="center" valign="top">/</td>
<td align="center" valign="top">6</td>
</tr>
<tr>
<td align="left" valign="top">Yanchi</td>
<td align="center" valign="top">/</td>
<td align="center" valign="top">/</td>
<td align="center" valign="top">1</td>
<td align="center" valign="top">1</td>
<td align="center" valign="top">12</td>
<td align="center" valign="top">34</td>
<td align="center" valign="top">48</td>
</tr>
<tr>
<td align="left" valign="top">Pengyang</td>
<td align="center" valign="top">/</td>
<td align="center" valign="top">/</td>
<td align="center" valign="top">/</td>
<td align="center" valign="top">2</td>
<td align="center" valign="top">/</td>
<td align="center" valign="top">/</td>
<td align="center" valign="top">2</td>
</tr>
<tr>
<td align="left" valign="top">Jingyuan</td>
<td align="center" valign="top">/</td>
<td align="center" valign="top">/</td>
<td align="center" valign="top">/</td>
<td align="center" valign="top">1</td>
<td align="center" valign="top">/</td>
<td align="center" valign="top">/</td>
<td align="center" valign="top">1</td>
</tr>
<tr>
<td align="left" valign="top">Haiyuan</td>
<td align="center" valign="top">/</td>
<td align="center" valign="top">/</td>
<td align="center" valign="top">/</td>
<td align="center" valign="top">7</td>
<td align="center" valign="top">/</td>
<td align="center" valign="top">/</td>
<td align="center" valign="top">7</td>
</tr>
<tr>
<td align="left" valign="top">Litong district</td>
<td align="center" valign="top">/</td>
<td align="center" valign="top">/</td>
<td align="center" valign="top">/</td>
<td align="center" valign="top">/</td>
<td align="center" valign="top">17</td>
<td align="center" valign="top">3</td>
<td align="center" valign="top">20</td>
</tr>
<tr>
<td align="left" valign="top">Qingtongxia district</td>
<td align="center" valign="top">/</td>
<td align="center" valign="top">/</td>
<td align="center" valign="top">/</td>
<td align="center" valign="top">/</td>
<td align="center" valign="top">10</td>
<td align="center" valign="top">2</td>
<td align="center" valign="top">12</td>
</tr>
<tr>
<td align="left" valign="top">Hongshibao</td>
<td align="center" valign="top">/</td>
<td align="center" valign="top">/</td>
<td align="center" valign="top">/</td>
<td align="center" valign="top">/</td>
<td align="center" valign="top">2</td>
<td align="center" valign="top">/</td>
<td align="center" valign="top">2</td>
</tr>
<tr>
<td align="left" valign="top">Total</td>
<td align="center" valign="top">1</td>
<td align="center" valign="top">80</td>
<td align="center" valign="top">28</td>
<td align="center" valign="top">29</td>
<td align="center" valign="top">44</td>
<td align="center" valign="top">43 (5<sup>#</sup>)</td>
<td align="center" valign="top">230</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<p>&#x201C;/&#x201D;: no strain. &#x201C;#&#x201D;: the location of one strain was unknown, and four strains were from border regions, including three in Inner Mongolia, and one in Shaanxi Province.</p>
</table-wrap-foot>
</table-wrap>
<table-wrap position="float" id="tab2">
<label>Table 2</label>
<caption><p>Host species, species/biovars, and proportions of isolated <italic>Brucella</italic> strains from 1950 to 2019.</p></caption>
<table frame="hsides" rules="groups">
<thead>
<tr>
<th align="left" valign="top" rowspan="2">Host</th>
<th align="center" valign="top" colspan="3"><italic>B. melitensis</italic></th>
<th align="center" valign="top" colspan="6"><italic>B. abortus</italic></th>
<th align="center" valign="top" colspan="2"><italic>B. suis</italic></th>
<th align="center" valign="top" rowspan="2"><italic>B. canis</italic></th>
<th align="center" valign="top" rowspan="2">Total</th>
<th align="center" valign="top" rowspan="2">Proportion (%)</th>
</tr>
<tr>
<th align="center" valign="top">bv. 1</th>
<th align="center" valign="top">bv. 2</th>
<th align="center" valign="top">bv. 3</th>
<th align="center" valign="top">bv. 1</th>
<th align="center" valign="top">bv. 3</th>
<th align="center" valign="top">bv. 4</th>
<th align="center" valign="top">bv. 5</th>
<th align="center" valign="top">bv. 6</th>
<th align="center" valign="top">bv. 7</th>
<th align="center" valign="top">bv. 1</th>
<th align="center" valign="top">bv. 3</th>
</tr>
</thead>
<tbody>
<tr>
<td align="left" valign="middle">Patient</td>
<td align="center" valign="middle">21</td>
<td align="center" valign="middle">8</td>
<td align="center" valign="middle">130</td>
<td align="center" valign="middle">/</td>
<td align="center" valign="middle">/</td>
<td align="center" valign="middle">/</td>
<td align="center" valign="middle">/</td>
<td align="center" valign="middle">/</td>
<td align="center" valign="middle">/</td>
<td align="center" valign="middle">4</td>
<td align="center" valign="middle">2</td>
<td align="center" valign="middle">/</td>
<td align="center" valign="top">165</td>
<td align="char" valign="middle" char=".">71.74%</td>
</tr>
<tr>
<td align="left" valign="middle">Cattle</td>
<td align="center" valign="middle">/</td>
<td align="center" valign="middle">1</td>
<td align="center" valign="middle">/</td>
<td align="center" valign="middle">3</td>
<td align="center" valign="middle">7</td>
<td align="center" valign="middle">1</td>
<td align="center" valign="middle">/</td>
<td align="center" valign="middle">3</td>
<td align="center" valign="middle">/</td>
<td align="center" valign="middle">/</td>
<td align="center" valign="middle">/</td>
<td align="center" valign="middle">/</td>
<td align="center" valign="top">15</td>
<td align="char" valign="middle" char=".">6.52%</td>
</tr>
<tr>
<td align="left" valign="middle">Sheep</td>
<td align="center" valign="middle">18</td>
<td align="center" valign="middle">18</td>
<td align="center" valign="middle">6</td>
<td align="center" valign="middle">/</td>
<td align="center" valign="middle">/</td>
<td align="center" valign="middle">/</td>
<td align="center" valign="middle">1</td>
<td align="center" valign="middle">1</td>
<td align="center" valign="middle">2</td>
<td align="center" valign="middle">/</td>
<td align="center" valign="middle">1</td>
<td align="center" valign="middle">/</td>
<td align="center" valign="top">47</td>
<td align="char" valign="middle" char=".">20.43%</td>
</tr>
<tr>
<td align="left" valign="middle">Dogs</td>
<td align="center" valign="middle">/</td>
<td align="center" valign="middle">1</td>
<td align="center" valign="middle">/</td>
<td align="center" valign="middle">/</td>
<td align="center" valign="middle">/</td>
<td align="center" valign="middle">/</td>
<td align="center" valign="middle">/</td>
<td align="center" valign="middle">/</td>
<td align="center" valign="middle">/</td>
<td align="center" valign="middle">/</td>
<td align="center" valign="middle">/</td>
<td align="center" valign="middle">2</td>
<td align="center" valign="top">3</td>
<td align="char" valign="middle" char=".">1.30%</td>
</tr>
<tr>
<td align="left" valign="middle">Total</td>
<td align="center" valign="middle">39</td>
<td align="center" valign="middle">28</td>
<td align="center" valign="middle">136</td>
<td align="center" valign="middle">3</td>
<td align="center" valign="middle">7</td>
<td align="center" valign="middle">1</td>
<td align="center" valign="middle">1</td>
<td align="center" valign="middle">4</td>
<td align="center" valign="middle">2</td>
<td align="center" valign="middle">4</td>
<td align="center" valign="middle">3</td>
<td align="center" valign="middle">2</td>
<td align="center" valign="top">230</td>
<td align="char" valign="middle" char=".">100%</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<p>&#x201C;/&#x201D;: no strain.</p>
</table-wrap-foot>
</table-wrap>
</sec>
<sec id="sec14">
<title>MLVA genotyping characteristics of 105 <italic>B. melitensis</italic> strains</title>
<p>Based on the <italic>HGDI</italic>, the <italic>MLVA-16</italic> panels exhibited a high discriminatory power for <italic>B. melitensis</italic>, with polymorphism levels of 0.9815 (<xref ref-type="table" rid="tab1">Table 1</xref>). The eight loci of panel 1, and <italic>bruce21</italic> of panel 2A, showed only one allele (<italic>HGDI</italic>&#x2009;=&#x2009;0). In contrast, the greatest variability was detected in panel 2B (diversity index &#x2267; 0.6548) within loci <italic>bruce04</italic> (<italic>HGDI</italic>&#x2009;=&#x2009;0.8190), <italic>bruce16</italic> (<italic>HGDI</italic>&#x2009;=&#x2009;0.8176), and <italic>bruce30</italic> (<italic>HGDI</italic>&#x2009;=&#x2009;0.6548). In Panel 1 (<italic>MLVA-8</italic>), all strains exhibited genotype 42, while two <italic>MLVA-11</italic> genotypes were identified: 116 (<italic>n</italic>&#x2009;=&#x2009;102) and 108 (<italic>n</italic>&#x2009;=&#x2009;3), with 116 being the dominant circulating genotype and widely distributed throughout the area. This data suggests that all strains belong to the East Mediterranean lineage with a common geographic origin. Furthermore, 105 strains clustered into 11 clusters (<italic>a</italic>&#x2009;~&#x2009;<italic>k</italic>) harboring the 61 <italic>MLVA-16</italic> genotypes (<italic>GTs</italic> 1&#x2013;61) (<xref ref-type="fig" rid="fig2">Figures 2</xref>, <xref ref-type="fig" rid="fig3">3</xref>), 39 of which were represented by unique strains. The other 22 were shared genotypes, with associated genotypes shared between two and ten isolates. The 22 shared genotypes comprised 66 clustered strains, with a clustering rate of 62.9% (66/105; <xref ref-type="fig" rid="fig3">Figure 3</xref>). These data revealed a co-existing transmission pattern, predominantly of epidemiologically related cases and limited sporadic cases. The most frequently observed genotype, genotype 27, was present in ten isolates obtained from three regions over a decade. The five strains containing the second most frequently observed genotypes, genotypes 18 and 25, were isolated from three regions. Seven shared GTs (GTs 1, 15, 19, 25, 29, 35, and 47) were each found in two to four isolates from different counties. These data suggest that multiple cross-regional transmission events were caused by a common source of infection. Another seven shared GTs (GTs 2, 11, 22, 23, 42, 53, and 54) consisted of strains from the same counties; additionally, strains with defined epidemiological links isolated from two families (<italic>NX201201</italic> and <italic>NX201206</italic>; and <italic>NX201303</italic> and <italic>NX201305</italic>) clustered into identical genotypes (<xref ref-type="fig" rid="fig3">Figure 3</xref>), demonstrating that multiple outbreak events within the same county were caused by a common source, respectively. Furthermore, 13 shared GTs containing strains isolated in different years (<xref ref-type="fig" rid="fig3">Figure 3</xref>) indicate the ongoing endemic presence of predominant circulating genotypes in this region.</p>
<fig position="float" id="fig3">
<label>Figure 3</label>
<caption><p>MLVA dendrogram of 105 strains from this study. The columns show the identification numbers (key), MLVA-16 genotypes (GT), and MLVA-11 genotypes (panels 1 and 2A), species/biovar., host species, geographic location, and the year of strain isolation. Regions of isolated strains are indicated by color, and shared MLVA-16 genotypes are indicated by black, bold, italic font.</p></caption>
<graphic xlink:href="fmicb-15-1320845-g003.tif"/>
</fig>
</sec>
<sec id="sec15">
<title>Investigation of molecular epidemiological relationships at the two scales</title>
<p>To investigate the molecular epidemiological relationships among strains, an <italic>MLVA-16</italic> comparison analysis was conducted on both national and international scales. A total of 1,330 <italic>B. melitensis</italic> strains from 30 provinces/regions, previously reported in our study, were utilized for this purpose. The <italic>MLVA</italic> comparison revealed that 21 shared genotypes (1&#x2013;21) existed between strains from Ningxia and those from 31 other provinces (<xref rid="SM4" ref-type="supplementary-material">Supplementary Table S4</xref>). Strains from these shared genotypes almost spanned all provinces, including Inner Mongolia, Xinjiang, Qinghai, Shanxi, Hebei, Liaoning, and Guangdong, as well as Guangxi, Jiangsu, Jiangxi, Fujian, Zhejiang, and Hainan provinces. Internationally, strains from this study exhibited identical <italic>MLVA-16</italic> genotypes with strains mostly from Kazakhstan (A, E, G, and H), and Mongolia (B,H), but also from Turkey (C), Syria (D), India (E), and Saudi Arabia (F) (<xref ref-type="fig" rid="fig4">Figure 4</xref>).</p>
<fig position="float" id="fig4">
<label>Figure 4</label>
<caption><p>The genetic comparison of strains from the present study to strains from global scale. MST for <italic>Brucella melitensis</italic> generated using the MLVA-16 data from 750 strains from East Mediterranean lineage from 32 countries (<xref rid="SM6" ref-type="supplementary-material">Supplementary Figure S2</xref>). Finally, 19 strains from seven countries (<xref rid="SM1" ref-type="supplementary-material">Supplementary Table S5</xref>) that have closely relationship to strains from the present study were selected to further compared the genetic relationship. The MST of strains were constructed by the PHYLOZIV 2.0 software, and color-coding according to country where strains isolated (pink in this study, blue in Kazakhstan, brown in Saudi Arabia, green in Mongolia, purple in Syria, orange in India, blackish green in Turkey), and circle size indicate the numbers of strains. The strains from this study and other countries formed eight <bold>(A&#x2013;H)</bold> complete match MLVA-16 genotypes, the number on the branch line represents the number of allele differences.</p></caption>
<graphic xlink:href="fmicb-15-1320845-g004.tif"/>
</fig>
</sec>
</sec>
<sec sec-type="discussion" id="sec16">
<title>Discussion</title>
<p>In this study, we performed a preliminary analysis of the epidemiological and clinical characteristics of human brucellosis in Ningxia. We observed three phases of human brucellosis prevalence: the epidemic phase, control phase, and re-emergence phase. The incidence of human brucellosis was high during 1955&#x2013;1978, then decreased dramatically during 1979&#x2013;1994, after which brucellosis reemerged in 1995. The observed re-emergence phase is consistent with that reported in Shanxi, in which the incidence increased dramatically from 7.0/100,000 in 2005 to 23.5/100,000 in 2014, with an average annual increase of 14.5% (<xref ref-type="bibr" rid="ref7">Chen et al., 2016</xref>). However, some regions have not experienced a re-emergence phase. Notably, no positive cases have been reported after 15&#x2009;years of continuous serological surveillance in Zhongwei county. Consequently, the brucellosis surveillance point in Zhongwei county was discontinued in 2005, resulting in a lack of prevention and control measures in this district. The re-emergence of human brucellosis is associated with many factors. Key contributing factors include the long-term persistence of infection in livestock (sheep and goats) and their movement. Additionally, with the continuous expansion of animal husbandry and the frequent introduction of livestock, resources for quarantine, immunization, and compensation are inadequate, resulting in the long-term persistence of animal brucellosis epidemics. Previous efforts to address this contributing factor have included newborn immunization of female lambs for three consecutive years between 2015 and 2017 to curb the rapid spread (<xref ref-type="bibr" rid="ref15">Duan et al., 2023</xref>). This resulted in the formation of an effective immunization barrier, and consequently, the incidence of human brucellosis decreased rapidly during this period. Therefore, the implementation of ongoing animal immunization is recommended in Ningxia.</p>
<p>Our results revealed that the incidence of brucellosis was concentrated in Yanchi county, Hongsibao district, and Tongxin county. These data raise questions about the effectiveness of public health measures. This is consistent with a previous study that indicated that high-incidence areas were mostly concentrated in central Ningxia and displayed a certain degree of correlation with the number of sheep in stock (<xref ref-type="bibr" rid="ref30">Zhao et al., 2020</xref>). These counties are renowned for sheep breeding; for instance, Yanchi County was designated as the National Animal and Poultry Genetic Resources Reserve in 2000. Furthermore, this county also successfully registered the Yanchi Tan Sheep breed as a National Brand within the national germplasm resources core protection area in 2021.</p>
<p>Moreover, Ningxia is home to the largest concentration of Hui people in China, whose diet and lifestyle practices align closely with Islamic traditions. The primary sources of income and animal protein for the local population are sheep and goat farming. While this industry forms the backbone of the county&#x2019;s poverty reduction efforts, the processing, sale, and consumption of mutton may pose infection risks to humans. This offers a plausible explanation for why the strains from Yanchi county led to multiple-point outbreaks of human brucellosis. A similar effect was observed in Shaanxi, where the increased brucellosis incidence was closely related to the development of large-scale small ruminant (goat and sheep) farms in Guanzhong and some southern Shaanxi regions (<xref ref-type="bibr" rid="ref3">An et al., 2021</xref>). Notably, Yanchi County is located in the easternmost region of Ningxia, adjacent to the Inner Mongolia Autonomous Region, Gansu (<xref ref-type="bibr" rid="ref26">Wei et al., 2019</xref>), and Shaanxi (<xref ref-type="bibr" rid="ref4">An et al., 2022</xref>), which are cluster areas of human brucellosis (<xref ref-type="bibr" rid="ref28">Yang et al., 2020</xref>). Previous studies showed that the highest seropositivity was observed in Yulin City, Shaanxi, which borders regions that have a high incidence of human brucellosis, including Inner Mongolia, Ningxia, and Gansu (<xref ref-type="bibr" rid="ref4">An et al., 2022</xref>). The high disease burden in this region raises questions about the effectiveness of public health measures. Therefore, further genome molecular association surveys of strains from the border regions are necessary to elucidate the transmission chain.</p>
<p>In this study, we observed high species/biovars and host diversity in Ningxia, constituting a relatively complete niche for the spread of <italic>Brucella</italic> strains. The potential for cross-host transmission of <italic>Brucella</italic> strains may create an optimal environment for their persistent circulation in this region. The results from MLVA genotyping suggested that the dominant circulating genotypes persistently expanded over an extended time span, contributing to both intra-county and cross-county human brucellosis epidemics. Moreover, the 39 singletons each represented one strain, indicating sporadic cases, including four strains from the border regions of Inner Mongolia and Shaanxi Province. This is consistent with the results of MLVA-16 genotyping of 70 <italic>B. melitensis</italic> isolates from Egypt, which revealed that 51 genotypes were represented by single isolates, while the remaining 19 genotypes were shared among 67 isolates, suggesting a combination of sporadic and epidemiologically related characteristics in <italic>B. melitensis</italic> infection (<xref ref-type="bibr" rid="ref25">Sayour et al., 2020</xref>). Further genomic epidemiological studies are necessary to support the present results and illustrate the phylogenetic relationships of <italic>B. melitensis</italic> strains. <italic>B. melitensis</italic> has been the dominant circulating population for nearly eight decades, responsible for the majority of human brucellosis cases. A similar high prevalence of <italic>B. melitensis</italic> has been observed in small ruminants in Turkey, as evidenced by isolations from sheep, goat, and cattle abortion samples from various farms across seven geographical regions (<xref ref-type="bibr" rid="ref2">Akar and Erganis, 2022</xref>). Furthermore, <italic>Brucella</italic> strains isolated from all regions of Kyrgyzstan were identified as <italic>B. melitensis</italic> (<xref ref-type="bibr" rid="ref18">Kydyshov et al., 2022</xref>). Notably, Guangxi province has witnessed considerable shifts in the species/biovar composition of <italic>Brucella</italic> strains, with previously predominant populations, <italic>B. suis</italic> and <italic>B. canis</italic>, being replaced by recently emerging <italic>B. melitensis</italic> strains (<xref ref-type="bibr" rid="ref21">Liu et al., 2019</xref>). Furthermore, the predominant circulating genotypes in this region were similar to strains from the Northwest, including Inner Mongolia, Xinjiang, and Qinghai. These data indicate that brucellosis has become a significant regional issue, with a high risk of spillover to southern areas due to animal trade and their products. Additionally, strains from this study showed complete consistency in <italic>MLVA-16</italic> genotypes with strains from other Asian countries, such as Kazakhstan and Mongolia. The breeding and trade of small ruminants have contributed to the strain&#x2019;s spread due to similar lifestyle and domestic farming practices in these agricultural and pastoral areas in China, Kazakhstan, and Mongolia. WGS-SNP analysis showed that strains from these countries lack clear territorial differentiation, demonstrating that strains continuously expand and spread in countries along the Silk Road (<xref ref-type="bibr" rid="ref21">Liu et al., 2019</xref>). The active exchange and trade of sheep and goats among these countries are one of the main driving factors. Further genomic epidemiology investigation of strains from multiple hosts will provide a new genetic map to better depict the transmission route and evolution profile of strains.</p>
<p>Based on the interpretations of our results, it is recommended to strengthen publicity and education and mobilize farmers and herders to consciously participate in brucellosis control work. First, enhancing communication and education collaboration with agricultural and animal husbandry departments to raise awareness should be a focus, particularly among high-risk groups. This should include promoting healthy lifestyles to reduce the epidemic of brucellosis (<xref ref-type="bibr" rid="ref8">Dadar et al., 2020</xref>; <xref ref-type="bibr" rid="ref13">Gong et al., 2023</xref>). Second, high-risk populations should be educated about the risks associated with privately slaughtering, trading, transferring, or consuming animals infected with brucellosis, those that died of disease, or those that died of unknown causes, especially during the spring calving and lambing season. They should be encouraged to implement rigorous personal protection and field disinfection measures (<xref ref-type="bibr" rid="ref1">Abu Shaqra, 2000</xref>). Additionally, the migration of sheep from high-to low-infection areas should be strictly limited to prevent further spread (<xref ref-type="bibr" rid="ref29">Zhang et al., 2014</xref>). Moreover, targeted prevention and control measures should be implemented, including collaborative publicity campaigns with local Islamic church organizations, to ensure that local people can participate in monitoring and control while respecting different cultural differences.</p>
<p>Although there were some new insights in this study, there may be potential limitations to be considered. Firstly, whole genome sequencing of some representative strains is needed to better illustrate the epidemiological transmission pattern and evolution characteristics of strains. Secondly, the limited number of strains from animals suggests that future studies should include larger sample sizes to verify the results of this study.</p>
</sec>
<sec sec-type="conclusions" id="sec17">
<title>Conclusion</title>
<p>Our analysis indicates an ongoing increase in the incidence rate of human brucellosis in Ningxia, from 0.441 (per 100,000) in 2005 to 86.83 (per 100,000) in 2022. Geographically, the affected area has continuously expanded, progressing from one affected county in 2004 to all 22 counties by 2022. <italic>B. melitensis</italic> emerged as the dominant circulating strain, responsible for the majority of cases. Molecular epidemiological investigations conducted using MLVA-16 genotyping have revealed the persistence of dominant genotypes in both intra-county and cross-county epidemics of human brucellosis. Brucellosis has become a serious regional issue, and there is a high risk of spillover to southern areas due to the trade in animals and their products.</p>
<p>The results of this study have provided a foundation for recommendations on the prevention and control of brucellosis epidemics. Further genomic epidemiology studies of circulating <italic>Brucella</italic> strains will provide new insights into the evolving epidemiological patterns of human brucellosis in Ningxia.</p>
</sec>
<sec sec-type="data-availability" id="sec18">
<title>Data availability statement</title>
<p>The original contributions presented in the study are included in the article/<xref ref-type="sec" rid="sec23">Supplementary material</xref>, further inquiries can be directed to the corresponding authors.</p>
</sec>
<sec sec-type="ethics-statement" id="sec19">
<title>Ethics statement</title>
<p>The studies involving humans were approved by Ethics Committee of the Ningxia Fourth People&#x2019;s Hospital (no. NXFPH-0624). The studies were conducted in accordance with the local legislation and institutional requirements. Written informed consent for participation in this study was provided by the participants&#x2019; legal guardians/next of kin.</p>
</sec>
<sec sec-type="author-contributions" id="sec20">
<title>Author contributions</title>
<p>GL: Methodology, Funding acquisition, Investigation, Writing &#x2013; review &#x0026; editing. XM: Investigation, Methodology, Writing &#x2013; review &#x0026; editing, Resources. RZ: Writing &#x2013; review &#x0026; editing, Formal analysis, Project administration. JL: Supervision, Writing &#x2013; review &#x0026; editing. PZ: Writing &#x2013; review &#x0026; editing, Data curation, Methodology. BL: Project administration, Writing &#x2013; review &#x0026; editing, Resources, Supervision. TL: Data curation, Writing &#x2013; review &#x0026; editing, Software. HR: Supervision, Writing &#x2013; review &#x0026; editing, Data curation, Formal analysis. ZgL: Data curation, Software, Conceptualization, Methodology, Writing &#x2013; original draft. ZjL: Supervision, Writing &#x2013; review &#x0026; editing, Formal analysis, Project administration. XJ: Project administration, Supervision, Writing &#x2013; review &#x0026; editing.</p>
</sec>
</body>
<back>
<sec sec-type="funding-information" id="sec21">
<title>Funding</title>
<p>The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This study was supported by the Open Project in State Key Laboratory of Infectious Disease Control and Prevention (no. 2023SKLD302), and Key Research and Development Project (no. 2021YFC2301001; Project Name: Control and Prevention of Emerging Infectious Diseases, Institute No. 20072002), and Key R&#x0026;D projects of Ningxia Hui Autonomous Region (2022BEG03161). The funders played no role in the study design, data collection and analysis, decision to publish, or preparation of the manuscript.</p>
</sec>
<sec sec-type="COI-statement" id="sec22">
<title>Conflict of interest</title>
<p>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.</p>
<p>The author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision.</p>
</sec>
<sec id="sec100" sec-type="disclaimer">
<title>Publisher&#x2019;s note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
<sec sec-type="supplementary-material" id="sec23">
<title>Supplementary material</title>
<p>The Supplementary material for this article can be found online at: <ext-link xlink:href="https://www.frontiersin.org/articles/10.3389/fmicb.2024.1320845/full#supplementary-material" ext-link-type="uri">https://www.frontiersin.org/articles/10.3389/fmicb.2024.1320845/full#supplementary-material</ext-link></p>
<supplementary-material xlink:href="Table_1.xlsx" id="SM1" mimetype="application/vnd.openxmlformats-officedocument.spreadsheetml.sheet" xmlns:xlink="http://www.w3.org/1999/xlink"/>
<supplementary-material xlink:href="Table_2.xlsx" id="SM2" mimetype="application/vnd.openxmlformats-officedocument.spreadsheetml.sheet" xmlns:xlink="http://www.w3.org/1999/xlink"/>
<supplementary-material xlink:href="Table_3.xlsx" id="SM3" mimetype="application/vnd.openxmlformats-officedocument.spreadsheetml.sheet" xmlns:xlink="http://www.w3.org/1999/xlink"/>
<supplementary-material xlink:href="Table_4.xlsx" id="SM4" mimetype="application/vnd.openxmlformats-officedocument.spreadsheetml.sheet" xmlns:xlink="http://www.w3.org/1999/xlink"/>
<supplementary-material xlink:href="Image_1.TIF" id="SM5" mimetype="image/tiff" xmlns:xlink="http://www.w3.org/1999/xlink"/>
<supplementary-material xlink:href="Image_2.TIF" id="SM6" mimetype="image/tiff" xmlns:xlink="http://www.w3.org/1999/xlink"/>
</sec>
<ref-list>
<title>References</title>
<ref id="ref1"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Abu Shaqra</surname> <given-names>Q. M.</given-names></name></person-group> (<year>2000</year>). <article-title>Epidemiological aspects of brucellosis in Jordan</article-title>. <source>Eur. J. Epidemiol.</source> <volume>16</volume>, <fpage>581</fpage>&#x2013;<lpage>584</lpage>. doi: <pub-id pub-id-type="doi">10.1023/a:1007688925027</pub-id>, PMID: <pub-id pub-id-type="pmid">11049102</pub-id></citation></ref>
<ref id="ref2"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Akar</surname> <given-names>K.</given-names></name> <name><surname>Erganis</surname> <given-names>O.</given-names></name></person-group> (<year>2022</year>). <article-title>Evaluation of the genetic profiles of Brucella melitensis strain from Turkey using multilocus variable number tandem repeat analysis (MLVA) and multilocus sequence typing (MLST) techniques</article-title>. <source>Vet. Microbiol.</source> <volume>269</volume>:<fpage>109423</fpage>. doi: <pub-id pub-id-type="doi">10.1016/j.vetmic.2022.109423</pub-id>, PMID: <pub-id pub-id-type="pmid">35462118</pub-id></citation></ref>
<ref id="ref3"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>An</surname> <given-names>C. H.</given-names></name> <name><surname>Liu</surname> <given-names>Z. G.</given-names></name> <name><surname>Nie</surname> <given-names>S. M.</given-names></name> <name><surname>Sun</surname> <given-names>Y. X.</given-names></name> <name><surname>Fan</surname> <given-names>S. P.</given-names></name> <name><surname>Luo</surname> <given-names>B. Y.</given-names></name> <etal/></person-group>. (<year>2021</year>). <article-title>Changes in the epidemiological characteristics of human brucellosis in Shaanxi Province from 2008 to 2020</article-title>. <source>Sci. Rep.</source> <volume>11</volume>:<fpage>17367</fpage>. doi: <pub-id pub-id-type="doi">10.1038/s41598-021-96774-x</pub-id>, PMID: <pub-id pub-id-type="pmid">34462483</pub-id></citation></ref>
<ref id="ref4"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>An</surname> <given-names>C. H.</given-names></name> <name><surname>Nie</surname> <given-names>S. M.</given-names></name> <name><surname>Sun</surname> <given-names>Y. X.</given-names></name> <name><surname>Fan</surname> <given-names>S. P.</given-names></name> <name><surname>Luo</surname> <given-names>B. Y.</given-names></name> <name><surname>Li</surname> <given-names>Z.</given-names></name> <etal/></person-group>. (<year>2022</year>). <article-title>Seroprevalence trend of human brucellosis and MLVA genotyping characteristics of Brucella melitensis in Shaanxi Province, China, during 2008&#x2013;2020</article-title>. <source>Transbound. Emerg. Dis.</source> <volume>69</volume>, <fpage>e423</fpage>&#x2013;<lpage>e434</lpage>. doi: <pub-id pub-id-type="doi">10.1111/tbed.14320</pub-id></citation></ref>
<ref id="ref5"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Araj</surname> <given-names>G. F.</given-names></name></person-group> (<year>2010</year>). <article-title>Update on laboratory diagnosis of human brucellosis</article-title>. <source>Int. J. Antimicrob. Agents</source> <volume>36</volume>, <fpage>S12</fpage>&#x2013;<lpage>S17</lpage>. doi: <pub-id pub-id-type="doi">10.1016/j.ijantimicag.2010.06.014</pub-id>, PMID: <pub-id pub-id-type="pmid">20692128</pub-id></citation></ref>
<ref id="ref6"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bagheri Nejad</surname> <given-names>R.</given-names></name> <name><surname>Krecek</surname> <given-names>R. C.</given-names></name> <name><surname>Khalaf</surname> <given-names>O. H.</given-names></name> <name><surname>Hailat</surname> <given-names>N.</given-names></name> <name><surname>Arenas-Gamboa</surname> <given-names>A. M.</given-names></name></person-group> (<year>2020</year>). <article-title>Brucellosis in the Middle East: current situation and a pathway forward</article-title>. <source>PLoS Negl. Trop. Dis.</source> <volume>14</volume>:<fpage>e0008071</fpage>. doi: <pub-id pub-id-type="doi">10.1371/journal.pntd.0008071</pub-id>, PMID: <pub-id pub-id-type="pmid">32437346</pub-id></citation></ref>
<ref id="ref9001"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bricker</surname> <given-names>B. J.</given-names></name> <name><surname>Halling</surname> <given-names>S. M.</given-names></name></person-group> (<year>1994</year>). <article-title>Differentiation of Brucella abortus bv. 1, 2, and 4, Brucella melitensis, Brucella ovis, and Brucella suis bv. 1 by PCR</article-title>. <source>J Clin Microbiol</source> <volume>32</volume>, <fpage>2660</fpage>&#x2013;<lpage>2666</lpage>. doi: <pub-id pub-id-type="doi">10.1128/jcm.32.11.2660-2666.1994</pub-id></citation></ref>
<ref id="ref7"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Chen</surname> <given-names>Q.</given-names></name> <name><surname>Lai</surname> <given-names>S.</given-names></name> <name><surname>Yin</surname> <given-names>W.</given-names></name> <name><surname>Zhou</surname> <given-names>H.</given-names></name> <name><surname>Li</surname> <given-names>Y.</given-names></name> <name><surname>Mu</surname> <given-names>D.</given-names></name> <etal/></person-group>. (<year>2016</year>). <article-title>Epidemic characteristics, high-risk townships and space-time clusters of human brucellosis in Shanxi Province of China, 2005-2014</article-title>. <source>BMC Infect. Dis.</source> <volume>16</volume>:<fpage>760</fpage>. doi: <pub-id pub-id-type="doi">10.1186/s12879-016-2086-x</pub-id></citation></ref>
<ref id="ref8"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Dadar</surname> <given-names>M.</given-names></name> <name><surname>Fakhri</surname> <given-names>Y.</given-names></name> <name><surname>Shahali</surname> <given-names>Y.</given-names></name> <name><surname>Mousavi Khaneghah</surname> <given-names>A.</given-names></name></person-group> (<year>2020</year>). <article-title>Contamination of milk and dairy products by Brucella species: a global systematic review and meta-analysis</article-title>. <source>Food Res. Int.</source> <volume>128</volume>:<fpage>108775</fpage>. doi: <pub-id pub-id-type="doi">10.1016/j.foodres.2019.108775</pub-id></citation></ref>
<ref id="ref9"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>De Massis</surname> <given-names>F.</given-names></name> <name><surname>Ancora</surname> <given-names>M.</given-names></name> <name><surname>Atzeni</surname> <given-names>M.</given-names></name> <name><surname>Rolesu</surname> <given-names>S.</given-names></name> <name><surname>Bandino</surname> <given-names>E.</given-names></name> <name><surname>Danzetta</surname> <given-names>M. L.</given-names></name> <etal/></person-group>. (<year>2015</year>). <article-title>MLVA as an epidemiological tool to trace Back Brucella melitensis Biovar 1 re-emergence in Italy</article-title>. <source>Transbound. Emerg. Dis.</source> <volume>62</volume>, <fpage>463</fpage>&#x2013;<lpage>469</lpage>. doi: <pub-id pub-id-type="doi">10.1111/tbed.12397</pub-id>, PMID: <pub-id pub-id-type="pmid">26194658</pub-id></citation></ref>
<ref id="ref10"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Deqiu</surname> <given-names>S.</given-names></name> <name><surname>Donglou</surname> <given-names>X.</given-names></name> <name><surname>Jiming</surname> <given-names>Y.</given-names></name></person-group> (<year>2002</year>). <article-title>Epidemiology and control of brucellosis in China</article-title>. <source>Vet. Microbiol.</source> <volume>90</volume>, <fpage>165</fpage>&#x2013;<lpage>182</lpage>. doi: <pub-id pub-id-type="doi">10.1016/s0378-1135(02)00252-3</pub-id>, PMID: <pub-id pub-id-type="pmid">12414142</pub-id></citation></ref>
<ref id="ref11"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Di Bonaventura</surname> <given-names>G.</given-names></name> <name><surname>Angeletti</surname> <given-names>S.</given-names></name> <name><surname>Ianni</surname> <given-names>A.</given-names></name> <name><surname>Petitti</surname> <given-names>T.</given-names></name> <name><surname>Gherardi</surname> <given-names>G.</given-names></name></person-group> (<year>2021</year>). <article-title>Microbiological laboratory diagnosis of human brucellosis: an overview</article-title>. <source>Pathogens</source> <volume>10</volume>:<fpage>121623</fpage>. doi: <pub-id pub-id-type="doi">10.3390/pathogens10121623</pub-id>, PMID: <pub-id pub-id-type="pmid">34959578</pub-id></citation></ref>
<ref id="ref15"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Duan</surname> <given-names>H.</given-names></name> <name><surname>Wu</surname> <given-names>X.</given-names></name> <name><surname>Ma</surname> <given-names>T.</given-names></name> <name><surname>Sun</surname> <given-names>J.</given-names></name> <name><surname>Qi</surname> <given-names>R.</given-names></name> <name><surname>Fang</surname> <given-names>Y.</given-names></name></person-group> (<year>2023</year>). <article-title>Epidemiological analysis of human brucellosis in Ningxia Hui autonomous region from 2004 to 2019</article-title>. <source>Chin. J. Endemiol.</source> <volume>41</volume>, <fpage>580</fpage>&#x2013;<lpage>584</lpage>. doi: <pub-id pub-id-type="doi">10.3760/cma.j.cn231583-20210429-00145</pub-id></citation></ref>
<ref id="ref12"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Franco</surname> <given-names>M. P.</given-names></name> <name><surname>Mulder</surname> <given-names>M.</given-names></name> <name><surname>Gilman</surname> <given-names>R. H.</given-names></name> <name><surname>Smits</surname> <given-names>H. L.</given-names></name></person-group> (<year>2007</year>). <article-title>Human brucellosis</article-title>. <source>Lancet Infect. Dis.</source> <volume>7</volume>, <fpage>775</fpage>&#x2013;<lpage>786</lpage>. doi: <pub-id pub-id-type="doi">10.1016/s1473-3099(07)70286-4</pub-id>, PMID: <pub-id pub-id-type="pmid">18045560</pub-id></citation></ref>
<ref id="ref13"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Gong</surname> <given-names>W.</given-names></name> <name><surname>Sun</surname> <given-names>P.</given-names></name> <name><surname>Zhai</surname> <given-names>C.</given-names></name> <name><surname>Yuan</surname> <given-names>J.</given-names></name> <name><surname>Chen</surname> <given-names>Y.</given-names></name> <name><surname>Chen</surname> <given-names>Q.</given-names></name> <etal/></person-group>. (<year>2023</year>). <article-title>Accessibility of the three-year comprehensive prevention and control of brucellosis in Ningxia: a mathematical modeling study</article-title>. <source>BMC Infect. Dis.</source> <volume>23</volume>:<fpage>292</fpage>. doi: <pub-id pub-id-type="doi">10.1186/s12879-023-08270-4</pub-id>, PMID: <pub-id pub-id-type="pmid">37147629</pub-id></citation></ref>
<ref id="ref14"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Gwida</surname> <given-names>M.</given-names></name> <name><surname>Al Dahouk</surname> <given-names>S.</given-names></name> <name><surname>Melzer</surname> <given-names>F.</given-names></name> <name><surname>R&#x00F6;sler</surname> <given-names>U.</given-names></name> <name><surname>Neubauer</surname> <given-names>H.</given-names></name> <name><surname>Tomaso</surname> <given-names>H.</given-names></name></person-group> (<year>2010</year>). <article-title>Brucellosis - regionally emerging zoonotic disease?</article-title> <source>Croat. Med. J.</source> <volume>51</volume>, <fpage>289</fpage>&#x2013;<lpage>295</lpage>. doi: <pub-id pub-id-type="doi">10.3325/cmj.2010.51.289</pub-id>, PMID: <pub-id pub-id-type="pmid">20718081</pub-id></citation></ref>
<ref id="ref9002"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Hunter</surname> <given-names>P. R.</given-names></name></person-group> (<year>1990</year>). <article-title>Reproducibility and indices of discriminatory power of microbial typing methods</article-title>. <source>J Clin Microbiol</source> <volume>28</volume>, <fpage>1903</fpage>&#x2013;<lpage>1905</lpage>. doi: <pub-id pub-id-type="doi">10.1128/jcm.28.9.1903-1905.1990</pub-id></citation></ref>
<ref id="ref16"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Jamil</surname> <given-names>T.</given-names></name> <name><surname>Khan</surname> <given-names>A. U.</given-names></name> <name><surname>Saqib</surname> <given-names>M.</given-names></name> <name><surname>Hussain</surname> <given-names>M. H.</given-names></name> <name><surname>Melzer</surname> <given-names>F.</given-names></name> <name><surname>Rehman</surname> <given-names>A.</given-names></name> <etal/></person-group>. (<year>2021</year>). <article-title>Animal and human brucellosis in Pakistan</article-title>. <source>Front. Public Health</source> <volume>9</volume>:<fpage>660508</fpage>. doi: <pub-id pub-id-type="doi">10.3389/fpubh.2021.660508</pub-id>, PMID: <pub-id pub-id-type="pmid">34395357</pub-id></citation></ref>
<ref id="ref17"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kili&#x00E7;</surname> <given-names>S.</given-names></name> <name><surname>Ivanov</surname> <given-names>I. N.</given-names></name> <name><surname>Durmaz</surname> <given-names>R.</given-names></name> <name><surname>Bayraktar</surname> <given-names>M. R.</given-names></name> <name><surname>Ayaslioglu</surname> <given-names>E.</given-names></name> <name><surname>Uyanik</surname> <given-names>M. H.</given-names></name> <etal/></person-group>. (<year>2011</year>). <article-title>Multiple-locus variable-number tandem-repeat analysis genotyping of human Brucella isolates from Turkey</article-title>. <source>J. Clin. Microbiol.</source> <volume>49</volume>, <fpage>3276</fpage>&#x2013;<lpage>3283</lpage>. doi: <pub-id pub-id-type="doi">10.1128/jcm.02538-10</pub-id></citation></ref>
<ref id="ref18"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kydyshov</surname> <given-names>K.</given-names></name> <name><surname>Usenbaev</surname> <given-names>N.</given-names></name> <name><surname>Berdiev</surname> <given-names>S.</given-names></name> <name><surname>Dzhaparova</surname> <given-names>A.</given-names></name> <name><surname>Abidova</surname> <given-names>A.</given-names></name> <name><surname>Kebekbaeva</surname> <given-names>N.</given-names></name> <etal/></person-group>. (<year>2022</year>). <article-title>First record of the human infection of Brucella melitensis in Kyrgyzstan: evidence from whole-genome sequencing-based analysis</article-title>. <source>Infect. Dis. Poverty</source> <volume>11</volume>:<fpage>120</fpage>. doi: <pub-id pub-id-type="doi">10.1186/s40249-022-01044-1</pub-id>, PMID: <pub-id pub-id-type="pmid">36482466</pub-id></citation></ref>
<ref id="ref19"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Lai</surname> <given-names>S.</given-names></name> <name><surname>Zhou</surname> <given-names>H.</given-names></name> <name><surname>Xiong</surname> <given-names>W.</given-names></name> <name><surname>Gilbert</surname> <given-names>M.</given-names></name> <name><surname>Huang</surname> <given-names>Z.</given-names></name> <name><surname>Yu</surname> <given-names>J.</given-names></name> <etal/></person-group>. (<year>2017</year>). <article-title>Changing epidemiology of human brucellosis, China, 1955-2014</article-title>. <source>Emerg. Infect. Dis.</source> <volume>23</volume>, <fpage>184</fpage>&#x2013;<lpage>194</lpage>. doi: <pub-id pub-id-type="doi">10.3201/eid2302.151710</pub-id>, PMID: <pub-id pub-id-type="pmid">28098531</pub-id></citation></ref>
<ref id="ref20"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Laine</surname> <given-names>C. G.</given-names></name> <name><surname>Scott</surname> <given-names>H. M.</given-names></name> <name><surname>Arenas-Gamboa</surname> <given-names>A. M.</given-names></name></person-group> (<year>2022</year>). <article-title>Human brucellosis: widespread information deficiency hinders an understanding of global disease frequency</article-title>. <source>PLoS Negl. Trop. Dis.</source> <volume>16</volume>:<fpage>e0010404</fpage>. doi: <pub-id pub-id-type="doi">10.1371/journal.pntd.0010404</pub-id>, PMID: <pub-id pub-id-type="pmid">35580076</pub-id></citation></ref>
<ref id="ref9003"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Le Fl&#x00E8;che</surname> <given-names>P.</given-names></name> <name><surname>Jacques</surname> <given-names>I.</given-names></name> <name><surname>Grayon</surname> <given-names>M.</given-names></name> <name><surname>Al Dahouk</surname> <given-names>S.</given-names></name> <name><surname>Bouchon</surname> <given-names>P.</given-names></name> <name><surname>Denoeud</surname> <given-names>F.</given-names></name> <etal/></person-group>. (<year>2006</year>). <article-title>Evaluation and selection of tandem repeat loci for a Brucella MLVA typing assay</article-title>. <source>BMC Microbiol</source> <volume>6</volume>:<fpage>9</fpage>. doi: <pub-id pub-id-type="doi">10.1186/1471-2180-6-9</pub-id></citation></ref>
<ref id="ref9004"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Li</surname> <given-names>Z.</given-names></name> <name><surname>Wang</surname> <given-names>X. M.</given-names></name> <name><surname>Zhu</surname> <given-names>X.</given-names></name> <name><surname>Wang</surname> <given-names>M.</given-names></name> <name><surname>Cheng</surname> <given-names>H.</given-names></name> <name><surname>Li</surname> <given-names>D.</given-names></name> <etal/></person-group> (<year>2020</year>). <article-title>Molecular Characteristics of Brucella Isolates Collected From Humans in Hainan Province</article-title>, <source>China. Front Microbiol</source> <volume>11</volume>:<fpage>452</fpage>. doi: <pub-id pub-id-type="doi">10.3389/fmicb.2020.00452</pub-id></citation></ref>
<ref id="ref21"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Liu</surname> <given-names>Z. G.</given-names></name> <name><surname>Wang</surname> <given-names>M.</given-names></name> <name><surname>Zhao</surname> <given-names>H. Y.</given-names></name> <name><surname>Piao</surname> <given-names>D. R.</given-names></name> <name><surname>Jiang</surname> <given-names>H.</given-names></name> <name><surname>Li</surname> <given-names>Z. J.</given-names></name></person-group> (<year>2019</year>). <article-title>Investigation of the molecular characteristics of Brucella isolates from Guangxi Province</article-title>. <source>China. BMC Microbiol</source> <volume>19</volume>:<fpage>292</fpage>. doi: <pub-id pub-id-type="doi">10.1186/s12866-019-1665-6</pub-id>, PMID: <pub-id pub-id-type="pmid">31842756</pub-id></citation></ref>
<ref id="ref9005"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Luc</surname> <given-names>A.</given-names></name></person-group> (<year>1995</year>). <article-title>Local indicators of spatial association-LISA</article-title>. <source>Geographical Analysis</source> <volume>27</volume>, <fpage>93</fpage>&#x2013;<lpage>115</lpage>. doi: <pub-id pub-id-type="doi">10.1111/j.1538-4632.1995.tb00338.x</pub-id></citation></ref>
<ref id="ref22"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Matle</surname> <given-names>I.</given-names></name> <name><surname>Ledwaba</surname> <given-names>B.</given-names></name> <name><surname>Madiba</surname> <given-names>K.</given-names></name> <name><surname>Makhado</surname> <given-names>L.</given-names></name> <name><surname>Jambwa</surname> <given-names>K.</given-names></name> <name><surname>Ntushelo</surname> <given-names>N.</given-names></name></person-group> (<year>2021</year>). <article-title>Characterisation of Brucella species and biovars in South Africa between 2008 and 2018 using laboratory diagnostic data</article-title>. <source>Vet Med Sci</source> <volume>7</volume>, <fpage>1245</fpage>&#x2013;<lpage>1253</lpage>. doi: <pub-id pub-id-type="doi">10.1002/vms3.483</pub-id>, PMID: <pub-id pub-id-type="pmid">33974356</pub-id></citation></ref>
<ref id="ref23"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Nascimento</surname> <given-names>M.</given-names></name> <name><surname>Sousa</surname> <given-names>A.</given-names></name> <name><surname>Ramirez</surname> <given-names>M.</given-names></name> <name><surname>Francisco</surname> <given-names>A. P.</given-names></name> <name><surname>Carri&#x00E7;o</surname> <given-names>J. A.</given-names></name> <name><surname>Vaz</surname> <given-names>C.</given-names></name></person-group> (<year>2016</year>). <article-title>PHYLOViZ 2.0: providing scalable data integration and visualization for multiple phylogenetic inference methods</article-title>. <source>Bioinformatics</source> <volume>33</volume>, <fpage>128</fpage>&#x2013;<lpage>129</lpage>. doi: <pub-id pub-id-type="doi">10.1093/bioinformatics/btw582</pub-id>, PMID: <pub-id pub-id-type="pmid">27605102</pub-id></citation></ref>
<ref id="ref24"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Pappas</surname> <given-names>G.</given-names></name> <name><surname>Papadimitriou</surname> <given-names>P.</given-names></name> <name><surname>Akritidis</surname> <given-names>N.</given-names></name> <name><surname>Christou</surname> <given-names>L.</given-names></name> <name><surname>Tsianos</surname> <given-names>E. V.</given-names></name></person-group> (<year>2006</year>). <article-title>The new global map of human brucellosis</article-title>. <source>Lancet Infect. Dis.</source> <volume>6</volume>, <fpage>91</fpage>&#x2013;<lpage>99</lpage>. doi: <pub-id pub-id-type="doi">10.1016/s1473-3099(06)70382-6</pub-id>, PMID: <pub-id pub-id-type="pmid">16439329</pub-id></citation></ref>
<ref id="ref25"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Sayour</surname> <given-names>A. E.</given-names></name> <name><surname>Elbauomy</surname> <given-names>E.</given-names></name> <name><surname>Abdel-Hamid</surname> <given-names>N. H.</given-names></name> <name><surname>Mahrous</surname> <given-names>A.</given-names></name> <name><surname>Carychao</surname> <given-names>D.</given-names></name> <name><surname>Cooley</surname> <given-names>M. B.</given-names></name> <etal/></person-group>. (<year>2020</year>). <article-title>MLVA fingerprinting of Brucella melitensis circulating among livestock and cases of sporadic human illness in Egypt</article-title>. <source>Transbound. Emerg. Dis.</source> <volume>67</volume>, <fpage>2435</fpage>&#x2013;<lpage>2445</lpage>. doi: <pub-id pub-id-type="doi">10.1111/tbed.13581</pub-id>, PMID: <pub-id pub-id-type="pmid">32304280</pub-id></citation></ref>
<ref id="ref26"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Wei</surname> <given-names>K. F.</given-names></name> <name><surname>Zhang</surname> <given-names>H.</given-names></name> <name><surname>He</surname> <given-names>J.</given-names></name> <name><surname>Gou</surname> <given-names>F. X.</given-names></name> <name><surname>Cheng</surname> <given-names>Y.</given-names></name> <name><surname>Liu</surname> <given-names>X. F.</given-names></name></person-group> (<year>2019</year>). <article-title>Spatial distribution of brucellosis in Gansu province, 2013&#x2013;2018</article-title>. <source>Zhonghua Liu Xing Bing Xue Za Zhi</source> <volume>40</volume>, <fpage>1099</fpage>&#x2013;<lpage>1105</lpage>. doi: <pub-id pub-id-type="doi">10.3760/cma.j.issn.0254-6450.2019.09.015</pub-id>, PMID: <pub-id pub-id-type="pmid">31594153</pub-id></citation></ref>
<ref id="ref27"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Yagupsky</surname> <given-names>P.</given-names></name> <name><surname>Morata</surname> <given-names>P.</given-names></name> <name><surname>Colmenero</surname> <given-names>J. D.</given-names></name></person-group> (<year>2019</year>). <article-title>Laboratory diagnosis of human brucellosis</article-title>. <source>Clin. Microbiol. Rev.</source> <volume>33</volume>:<fpage>e19</fpage>. doi: <pub-id pub-id-type="doi">10.1128/cmr.00073-19</pub-id>, PMID: <pub-id pub-id-type="pmid">31722888</pub-id></citation></ref>
<ref id="ref28"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Yang</surname> <given-names>H.</given-names></name> <name><surname>Zhang</surname> <given-names>S.</given-names></name> <name><surname>Wang</surname> <given-names>T.</given-names></name> <name><surname>Zhao</surname> <given-names>C.</given-names></name> <name><surname>Zhang</surname> <given-names>X.</given-names></name> <name><surname>Hu</surname> <given-names>J.</given-names></name> <etal/></person-group>. (<year>2020</year>). <article-title>Epidemiological characteristics and spatiotemporal trend analysis of human brucellosis in China, 1950&#x2013;2018</article-title>. <source>Int. J. Environ. Res. Public Health</source> <volume>17</volume>:<fpage>72382</fpage>. doi: <pub-id pub-id-type="doi">10.3390/ijerph17072382</pub-id></citation></ref>
<ref id="ref29"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Zhang</surname> <given-names>J.</given-names></name> <name><surname>Yin</surname> <given-names>F.</given-names></name> <name><surname>Zhang</surname> <given-names>T.</given-names></name> <name><surname>Yang</surname> <given-names>C.</given-names></name> <name><surname>Zhang</surname> <given-names>X.</given-names></name> <name><surname>Feng</surname> <given-names>Z.</given-names></name> <etal/></person-group>. (<year>2014</year>). <article-title>Spatial analysis on human brucellosis incidence in mainland China: 2004&#x2013;2010</article-title>. <source>BMJ Open</source> <volume>4</volume>:<fpage>e004470</fpage>. doi: <pub-id pub-id-type="doi">10.1136/bmjopen-2013-004470</pub-id>, PMID: <pub-id pub-id-type="pmid">24713215</pub-id></citation></ref>
<ref id="ref30"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Zhao</surname> <given-names>Y.</given-names></name> <name><surname>Guo</surname> <given-names>Z. Q.</given-names></name> <name><surname>Zhao</surname> <given-names>J. H.</given-names></name> <name><surname>Liang</surname> <given-names>P. F.</given-names></name></person-group> (<year>2020</year>). <article-title>Spatial-temporal distribution of human brucellosis and its correlation with the number of livestock in Ningxia Hui autonomous region, 2012&#x2013;2018</article-title>. <source>Zhonghua Liu Xing Bing Xue Za Zhi</source> <volume>41</volume>, <fpage>872</fpage>&#x2013;<lpage>876</lpage>. doi: <pub-id pub-id-type="doi">10.3760/cma.j.cn112338-20190925-00701</pub-id>, PMID: <pub-id pub-id-type="pmid">32564552</pub-id></citation></ref>
<ref id="ref31"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Zhu</surname> <given-names>X.</given-names></name> <name><surname>Zhao</surname> <given-names>Z.</given-names></name> <name><surname>Ma</surname> <given-names>S.</given-names></name> <name><surname>Guo</surname> <given-names>Z.</given-names></name> <name><surname>Wang</surname> <given-names>M.</given-names></name> <name><surname>Li</surname> <given-names>Z.</given-names></name> <etal/></person-group>. (<year>2020</year>). <article-title>Brucella melitensis, a latent travel bacterium, continual spread and expansion from northern to southern China and its relationship to worldwide lineages</article-title>. <source>Emerg Microbes Infect</source> <volume>9</volume>, <fpage>1618</fpage>&#x2013;<lpage>1627</lpage>. doi: <pub-id pub-id-type="doi">10.1080/22221751.2020.1788995</pub-id>, PMID: <pub-id pub-id-type="pmid">32594852</pub-id></citation></ref>
</ref-list>
</back>
</article>
