Skip to main content

The spatiotemporal distribution of human brucellosis in mainland China from 2007-2016



Despite the considerable efforts made to address the issue of brucellosis worldwide, its prevalence in dairy products continues to be difficult to estimate and represents a key public health issue around the world today. The aim of the present study was to better understand the epidemiology of this disease in mainland China. We set out to investigate the yearly spatial distribution and possible hotspots of the disease.


Human brucellosis data from mainland China between 2007 and 2016 were collected from the China Information System for Disease Control and Prevention. A geographic information system ArcGIS10.3 (ESRI, Redlands) was used to identify potential changes in the spatial and temporal distribution of human brucellosis in mainland China during the study period. These distributions were evaluated using three-dimensional trend analysis and spatial autocorrelation analyse. A gravity-center was used to analyse the migration track of human brucellosis.


A total of 399,578 cases of human brucellosis were reported during the 10-year study period. The monthly incidence of brucellosis in China demonstrates clear seasonality. Spring and summer are the peak seasons, while May is the peak month for brucellosis. Three-dimensional trend analysis suggests that brucellosis is on the rise from south to north, and that the epidemic situation in northern China is more severe. Between 2007 and 2016, the overall migration distance of the brucellosis incidence gravity-center was 906.43 km, and the direction was southwest. However, the overall gravity center of brucellosis was still in the northern part of China. In the global autocorrelation analysis, brucellosis in China demonstrated a non-random distribution between 2013 and 2014, with spatial autocorrelation (Z > 1.96, P < 0.05) and a clustering trend, while no clustering trend was found from 2007 to 2012 or from 2015 to 2016. In the local autocorrelation analysis, a Low-Low cluster phenomenon was found in the south of China in 2013 and 2014.


Human brucellosis remains a widespread challenge, particularly in northern China. The hotspots highlight potential high-risk areas which may require special plans and resources for monitoring and controlling the disease.

Peer Review reports


Brucellosis, which is also known as undulant fever, is one of the most common zoonoses in the world [1,2,3,4]. It is an infectious disease which can be acute or chronic. It is caused by brucella, which is a natural epidemic disease [5]. Brucellosis is globally distributed. More than 170 countries and regions around the world have recorded brucellosis epidemics, accounting for approximately one-sixth to one-fifth of the world’s population [6]. There are between five and six million patients who are living with brucellosis worldwide, the number of cases reported to the WHO each year is 500,000 [7,8,9,10]. The Mediterranean region as well as Asia, Central and South America are the areas with the highest incidences [11,12,13,14,15].

In China, the first incidence of brucellosis occurred for two foreign patients in Shanghai in 1905 and was recorded as Malta fever [16]. However, several patients in China with similar clinical symptoms had been observed in the 10 years before that report. Following the 1905 report, three more cases were reported in Chongqing in 1906 [17]. By the 1950s and 1960s, brucellosis had become more prevalent for Chinese people and animals. Subsequently, from the 1980s to the early 1990s, activities for the prevention and control of brucellosis have gradually been introduced in mainland China, meaning that the brucellosis epidemic had been basically controlled [18]. During this time, the human infection rate was just 3%, the incidence rate was only 0.02/100,000 [16]. However, brucellosis has rebounded since 1993, with sporadic cases increasing at a rate of between 30 and 50%. In 2015, the incidence of human cases reached 4.18/100,000. In the ranking of statutory infectious diseases, it rose from 16th place in 2000 to 6th place in 2014. A strong upward trend such as this is extremely rare among all such reported infectious diseases.

Previous studies have demonstrated that the global epidemiology of brucellosis has changed dramatically over the past decades, particularly in industrialized countries [19]. The spatial distribution characteristics of a brucellosis outbreak also have changed distinctly [20], although whether this change has temporal and spatial covariance requires further investigation.

The aim of our study was to use a spatial distribution model to analyze the spatial distribution characteristics and correlations of human brucellosis in mainland China between 2007 and 2016, and so to provide theoretical support for the prevention and control of human brucellosis in the context of local conditions.


Data source

Human brucellosis cases from 31 provinces (municipalities and autonomous regions) in mainland China between 2007 and 2016 were collected through National Population and Health Science Data Sharing Platform, the Data-center of China Public Health Science. The platform was established in 2004 and is more integrated, effective and reliable than the previous case-reporting system [21]. The copyright of data sharing platform belongs to Chinese center for disease control and prevention (CDC). The data on this platform is publicly available, anyone can download the data from platform after the account application being approved. In mainland China, suspected and confirmed cases of human brucellosis must be reported to local and provincial CDC and then to the Chinese CDC (CCDC). In order to meet the case definition, a diagnosis of human brucellosis must not only be accompanied by clinical signs, but must also be confirmed by a serological test or in isolation, in accordance with the case definition from the WHO [22].

Statistical analysis

Descriptive analysis

Human brucellosis cases which occurred between 2007 and 2016 were collected by month and year, and Excel software was used to draw a time series diagram. Data were sorted and described according to time and region.

Three-dimensional trend analysis

Both the spatial distribution and incidence trend of human brucellosis were evaluated and visualized via the three-dimensional trend analysis, utilizing ArcGIS10.3 (ESRI, Redlands) software [23]. The three-dimensional trend analysis uses the least square method to establish the binary polynomial regression model, which excludes both accidental and local variation, and displays the general distribution rule of disease in the study area [24]. The X-axis and Y-axis represent the geometric center of specific study region, while the Z-axis represents the human brucellosis. The higher the Z-axis, the higher the incidence [25]. This means that a point in three-dimensional space (X, Y, Z) represents human brucellosis incidences in a study region. The points in the three-dimensional space are projected onto both the XZ plane and the YZ plane respectively. Following this, an optimal fitting line is created according to the scatter plot, fitting polynomial onto the projection plane, while the curve is used to simulate the spatial trend of the disease [26].

Gravity-center migration

A common method for studying the migration of matter and energy in geo-distribution is through the use of the ‘gravity-center’. A ‘geographic gravity-center’ can represent the spatial and temporal distribution characteristics of geographical elements. The moving direction of a gravity-center points to the ‘high density’ element of spatial phenomena, while the distance reflects the spatial difference of the variation range of the geographical elements [27].


Using spatial statistical tools (such as the mean center) to calculate the incidence gravity-center of brucellosis each year, we assumed that a region is composed of several sub-regions. Each sub-region has a particle Pi (Xi, Yi) and its attribute value is Mi, meaning the center coordinate formula of an attribute is [28]:

$$ X=\frac{\varSigma_{i=1}^n{M}_i{X}_i}{\varSigma_{i=1}^n{M}_i}\kern2.12em Y=\frac{\varSigma_{i=1}^n{M}_i{Y}_i}{\varSigma_{i=1}^n{M}_i} $$

X and Y respectively represent the latitude and longitude of the gravity-center of a region’s attribute, while (Xi, Yi) is the geographic coordinate of the central city of the i-th secondary region and Mi is the magnitude of an attribute of the i-th secondary region.

Migration track

When a certain attribute value in the study area accounts for a large proportion or a rapid growth rate in the population, due to the unbalanced distribution of the spatial attributes, the regional gravity-center will move to this area. This means that a ‘gravity-center offset’ will occur. The gravity-center deviation demonstrates that the attribute develops rapidly in this direction, while the strength of the overall social activity attribute is enhanced. Through the regional gravity-center migration, we can see the migration track and spatial distribution difference of this region’s attribute, which is calculated as follows:

$$ {D}_{ab}=C\times \sqrt{{\left({X}_a-{X}_b\right)}^2+{\left({Y}_a-{Y}_b\right)}^2} $$

Here, a and b represent different years, D represents the distance of the gravity center offset between different years, (Xa,Xb) and (Ya, Yb) represent center coordinates of an attribute of the a and b years, respectively. C is a constant, representing the coefficient from the coordinate unit of the earth’s surface to the plane distance; C = 111.111.

Spatial autocorrelation analysis

Tobler’s first law of geography states that geographic objects or attributes are interrelated in spatial distribution. This means that the degree of interaction between things with a close geographical distance is higher than that of things with a long distance [29]. Spatial autocorrelation in GIS helps understand the degree to which one object is similar to other nearby objects. Moran’s I (Index) measures spatial autocorrelation. Spatial autocorrelation measures how much close objects are in comparison with other close objects. Moran’s I can be classified as positive, negative and no spatial auto-correlation. Positive spatial autocorrelation is when similar values cluster together in a map. Negative spatial autocorrelation is when dissimilar values cluster together in a map. This concept was extended to spatial epidemiology, which can be used to demonstrate the spatial distribution of diseases. Therefore, this analysis has its own advantages for recognizing the spatial distribution characteristics and pattern causes of diseases. Spatial autocorrelation can be divided into global and local spatial autocorrelation.

Global spatial autocorrelation can be used to describe the overall spatial distribution pattern of brucellosis in China [30]. Here, global Moran’s I value was used as a measure, which ranges from [− 1,1]. I > 0 indicates a positive spatial correlation; I < 0 indicates a negative spatial correlation; if I value is close to 0, no spatial correlation exists. The larger the absolute value of I, the stronger the correlation. When |Z| > 1.96, P < 0.05 was considered statistically significant and there was spatial autocorrelation, this is written in mathematical notation as:

$$ I=\frac{n\sum \limits_{i=1}^n\sum \limits_{j=1}^n{w}_{ij}\left({x}_i-\overline{x}\right)\left({x}_j-\overline{x}\right)}{\sum \limits_{i=1}^n\sum \limits_{j=1}^n{w}_{ij}\sum \limits_{i=1}^n{\left({x}_i-\overline{x}\right)}^2} $$

I is Moran’s I value, n refers to the number of spatial elements, xi and xj are the observed values at cell i and j, \( \overline{x} \) is the mean value of attribute values of all spatial units, while wij refers to the spatial weight between elements i and j. The expression of wij is as follows:

$$ {w}_{ij}=\left\{\begin{array}{l}1\ \mathrm{If}\ \mathrm{the}\ i-\mathrm{th}\ \mathrm{and}\ j-\mathrm{th}\ \mathrm{spatial}\ \mathrm{units}\ \mathrm{are}\ \mathrm{adjacent}\\ {}0\ \mathrm{If}\ \mathrm{the}\ i-\mathrm{th}\ \mathrm{and}\ j-\mathrm{th}\ \mathrm{spatial}\ \mathrm{units}\ \mathrm{are}\ \mathrm{not}\ \mathrm{adjacent}\\ {}0\ \mathrm{If}\ \mathrm{the}\ i-\mathrm{th}\ \mathrm{and}\ j-\mathrm{th}\ \mathrm{are}\ \mathrm{equal}\end{array}\right\} $$

In mathematical notation, this is:

$$ {z}_1=\frac{\left(I-E(I)\right)}{\sqrt{Var(I)}} $$

The corresponding variance is:

$$ Var(I)=\frac{1}{S_0^2\left({n}^2-1\right)}\left({n}^2{S}_1-n{S}_2+3{S}_0^2\right)-\frac{1}{{\left(n-1\right)}^2} $$

Among them:

$$ {S}_0={\sum}_{ij}{w}_{ij},{S}_1=2{\sum}_{ij}{w^2}_{ij},{S}_2=4{\sum}_i{w^2}_i,{w}_i={\sum}_j{w}_{ij} $$
$$ E(I)=-1/n-1 $$

Where E(I) refers to the mean value of spatial Morn’s I under the completely random assumption theoretical.

Local spatial autocorrelation: Local cluster analysis was performed using the local Moran’s I statistic [31]. Anselin’s local Moran’s I (local indicators of spatial association or LISA) statistic indicates the location of local clusters and spatial outliers [32, 33]. There are four cluster types in a local spatial autocorrelation analysis cluster map: High-High, Low-Low, High-Low and Low-High values, which are relative to neighboring districts [34, 35]. High-High and Low-Low clusters indicate that the observed values of an area are similar to those of the surrounding areas. Clusters of areas which have high crude brucellosis incidence rates (High-High or HH) are considered ‘hotspots’, whereas clusters of areas with low crude brucellosis incidence rates (Low-Low or LL) are ‘cold spots’. High-Low and Low-High clusters indicate that the observed values of an area are larger than (or less than) those of the surrounding areas. Clustering and outlier analysis were used to study local autocorrelation, as they can detect correlations between the disease occurrence of the provinces (municipalities and autonomous regions) as well as neighboring provinces (municipalities and autonomous regions) in mainland China and so identify areas where the disease is concentrated. The formula is as follows:

$$ I=\frac{n^2\left({x}_i-\overline{x}\right){\varSigma}_{j=1}^n{w}_{ij}\left({x}_j-x\right)}{\varSigma_{i=1}^n\left({x}_i-\overline{x}\right)2} $$

Here, the letters represent the same meanings as in the global Moran’s I coefficients.


Epidemic situation of human brucellosis

Time distribution

During the study period of 2007 to 2016, there were 399,578 reported cases of human brucellosis in mainland China. Annual case totals ranged from 19,721 in 2007 to 57,222 in 2014, with a mean number of approximately 34,498 (95% CI: 25610, 43,386) cases reported per year. The annual mean incidence of brucellosis was 2.97/100,000, with a low of 1.50/100,000 in 2007 and a high of 4.22/100,000 in 2014. As a whole, annual incidences follow a growth trend, with a linear relationship of y = 0.0211x + 1.7003 from 2007 to 2016 (Fig. 1). Furthermore, annual and monthly incidences fluctuated during this 10-year period. Between 2007 and 2009, the annual average incidence of brucellosis in China demonstrated an upward trend. Following that, a short-term decline occurred in 2010, after which incidences of the disease began to rise again, continuing to rise until 2014, and then decline again in 2015 and 2016. A purely temporal analysis between 2007 and 2016 demonstrated that high-incidence seasons of human brucellosis were the time between late spring and early summer, while the incidence was highest in May.

Fig. 1

Time distribution of human brucellosis incidences in China from 2007 to 2016. (The histogram shows monthly incidences and the dots show annual incidences of brucellosis. Annual incidences have a y = 0.0211x + 1.7003 relationship. Results of a purely temporal analysis between 2007 and 2016 demonstrated that high-incidence seasons of human brucellosis occurred between late spring and early summer, with incidences highest in May)

Spatial distribution

Incidences of human brucellosis in each province vary significantly, with an interquartile range of 4.25/100,000 (0.045~4.297/100,000) and a range of 70.10/100,000 (0~70.10/100,000). We defined regions with an incidence greater than 10/100,000 as high incidence areas, while incidences lower than 0.999/100,000 were defined as low incidence areas. Over the total study period, high incidence areas increased from two to five, with the number of high incidence provinces and cities increasing year by year. In 2007, there were six provinces (6/31, 19.35%) with no cases of human brucellosis. These were mainly located in central and southern regions such as Guizhou, Hunan, Jiangxi and Hainan. As time passed, the epidemic area gradually expanded into all provinces and cities in China. The spatial and temporal distribution map (Fig. 2) shows that 31 provinces (municipalities and autonomous regions) had incidences between 2007 and 2016. The epidemic mainly occurred in the north of China, such as the Inner Mongolia and Shanxi provinces. After 2009, the high incidence area grandly expanded in other northern areas, such as the Heilongjiang, Xinjiang and Ningxia provinces. In Inner Mongolia, disease incidence was the highest in history, exceeding the 50/100,000 population between 2009 and 2011. Between 2014 and 2016, incidences increased considerably in areas which had been historically low areas such as Gansu, Qinghai and Shandong. This was most marked in Ningxia, where incidences increased sharply from 0.4636/100,000 in 2007 to 43.6557/100,000 in 2015.

Fig. 2

Spatial distribution of human brucellosis incidences in mainland China from 2007 to 2016 (1/100,000). (The map data comes from National Geographic Center of China. After the website application, the data can be downloaded for free. The website link is

Three-dimensional trend analysis

The overall three-dimensional trend for 2007–2016 indicated an uneven distribution of human brucellosis in mainland China (Fig. 3). The curve on the XZ plane represents the trend of East-West direction, while the curve on the YZ plane represents the change of North-South direction. Each point represents the projection of incidence. It can be seen that incidences of human brucellosis in the northern region were higher than in the southern region. However, in the east-west direction, incidences of brucellosis were evenly distributed in mainland China.

Fig. 3

Three-dimensional trend analysis of mean incidences of brucellosis in mainland China between 2007 and 2016. (The curve on the XZ plane represents the East-West direction trend, while the curve on the YZ plane represents change in the North-South direction. Overall, we can see that incidences of brucellosis in northern China are relatively high, although in recent years, incidences in southern China have gradually increased)

Gravity center migration

Most cases of human brucellosis were reported in northern China during 2007 to 2016, while most provinces in northern China experienced a serious epidemic (Fig. 4). Between 2007 and 2012, the gravity-center of brucellosis was concentrated in the central area of Inner Mongolia. After 2012, there was a trend of gradual migration to the southwest. However, the distance and direction of annual migration were not obvious. Results showed that brucellosis was still prevalent in the northern part of China. During the study period, the gravity-center of brucellosis migrated 906.43 km southwest. The distance and direction of annual migration can be seen in Table 1.

Fig. 4

The center of gravity migration track of human brucellosis in China from 2007 to 2016. (The small black dots in the picture show the gravity center of brucellosis between 2007 and 2016. The map data comes from National Geographic Center of China. The website link is, Plan approval number is GS (2019)1815)

Table 1 Distance and direction of center of gravity migration of brucellosis between 2007 and 2016

Spatial autocorrelation analysis

Global spatial autocorrelation

Using provincial units to carry out the global autocorrelation analysis, we obtained Moran’s I value, variance, Z score and P value from 2007 to 2016, respectively (see Table 2). The values of Moran’s I were 0.1179 and 0.1181 respectively for 2013 and 2014, while Z values were greater than 1.96 (all P < 0.05), indicating that incidences of brucellosis in China between 2013 to 2014 had a non-random distribution, while there was no spatial autocorrelation in either 2007–2012 or 2015–2016 (Z < 1.96, P > 0.05, meaning there was no clustering trend.

Table 2 Global autocorrelation of Moran’s I values of brucellosis in mainland China from 2007 to 2016

Local spatial autocorrelation

While global clustering indices evaluate the distribution of events, they do not specify the location of hotspots. To identify areas characterized by a higher rate of human brucellosis incidences, a local spatial autocorrelation analysis was carried out for the years which had a statistical significance of Moran’s I value, and a local spatial autocorrelation aggregation distribution map was generated (Fig. 5). There were no High-High, High-Low or Low-High cluster areas in China in 2013, while Low-Low cluster areas tended to be concentrated in the southern provinces of China, including the provinces of Yunnan, Guizhou, Guangxi, Hunan, Guangdong, Hubei, Jiangxi, Fujian, Zhejiang and Anhui, with an average incidence of 0.08/100,000. The new High-Low cluster area in 2014 was Xinjiang, while the Low-Low cluster area decreased compared to 2013, now including Yunnan, Guizhou, Guangxi, Hunan, Guangdong, Jiangxi, Fujian, Zhejiang and Anhui provinces. The mean incidence of Low-Low clusters was 0.14/100,000 in 2014, meaning the incidence of Low-Low clusters in 2014 was higher than in 2013. Incidences in the Hubei Province increased rapidly from 0.08/100,000 in 2013 to 0.34/100,000 in 2014, thus reducing Hubei Province in Low-Low cluster areas in 2014. Cluster maps for 2013 and 2014 demonstrate that the cluster area in China was relatively stable, and that compared to the spatial distribution map of incidence, the Low-Low cluster area in each year tends to be the same as the low incidence area, all of which are located in the south of China.

Fig. 5

Local spatial autocorrelation analysis of incidences of human brucellosis in mainland China between 2013 and 2014. (Bright red indicates a High-Low association, bright blue indicates Low-High, light red indicates High-High and light blue indicates Low-Low. The map data comes from National Geographic Center of China. The website link is


Exploring the spatial distribution of a disease is beneficial for guiding the development and implementation of disease surveillance and a control plan. The application of ArcGIS10.3 (ESRI, Redlands) in disease surveillance has yielded good results. Based on surveillance data about human brucellosis between 2007 to 2016, the temporal and spatial distribution of incidences in mainland China suggested that brucellosis in mainland China has unique spatial and temporal clustering characteristics.

In recent years, incidences of brucellosis in mainland China have been relatively high and growing rapidly. Incidences in 2014 are about 2.8 times of that 2007. In addition, brucellosis has obvious seasonality. Incidences were higher in spring and summer and reached the highest in May every year. This is largely consistent with findings from related studies [11, 36,37,38]. The high incidence may be related to the vigorous development of animal husbandry in mainland China in the past 10 years, with a sharp increase in the number of livestock raised, frequent transactions and frequent transfer of livestock. However, grass-roots animal epidemic prevention departments are weak, and existing resources are not enough to meet the current increasing demand for epidemic prevention, making it difficult to deal with prevention and control measures such as immunization, quarantine and treatment of dead animals. High incidence rates of the disease in spring and summer are mainly due to the fact that it is the calving season, meaning that the chances of exposure to placenta increase. At this time, the contact probability of the relevant occupational personnel with abortion products and secretions of the sick livestock is greatly increased. If there are no effective protective measures for these occupational personnel, incidences will continue to increase.

It was found that brucellosis in China tends to show an upward trend from south to north, all of which began to spread from highly affected areas before gradually spreading outward. The gravity-center of brucellosis in China in recent years gradually moved southwest. However, overall migration distance and direction were small, and the gravity-center was always located in Inner Mongolia. The direction of migration was related to the location of the four major pastoral areas in China, which is similar to previously published research results [39]. High incidences of brucellosis were in the north of mainland China, where some areas showed much higher incidences than in the south of China. This may be because the north is the main gathering area of pastoral areas in China, meaning there are greater numbers of livestock being raised there as well as frequent livestock exchange.

A global spatial autocorrelation study demonstrated that brucellosis in China had a positive spatial correlation and clustering between 2013 and 2014, while there was no spatial autocorrelation between 2007 and 2012 or 2015 and 2016. Further local autocorrelation analysis found that the Low-Low concentration areas of brucellosis in China in 2013 were mainly in the southern part of the country, including the provinces of Yunnan, Guizhou, Guangxi, Hunan, Guangdong, Hubei, Jiangxi, Fujian, Zhejiang and Anhui, which is basically consistent with the low incidence areas. In 2014, the Low-Low cluster areas remained in the south of China. Understanding the spatial distribution of brucellosis in China is of great significance for controlling its occurrence and reducing the burden it causes, as well as for optimizing the allocation of resources. A prevention policy of key monitoring should be adopted for High-High cluster areas, while favorable factors for disease control should be found in time for Low-Low cluster areas, so as to provide guidance for infectious disease surveillance. In recent years, incidences of brucellosis in Low-Low areas in China have been clustered in fixed places, meaning the influencing factors of this clustering need further investigation.

At the same time, since brucellosis is a natural epidemic disease, it is not infectious between people, meaning the natural and social environment have a great impact on it. Some natural factors in neighbouring areas may make brucellosis spread more easily, leading to a clustering distribution of the affected areas. It was found that the areas with a high incidence of brucellosis tended to be concentrated in the north of China (Fig. 4). Previous studies [40,41,42] have found that besides raising livestock, brucellosis may also be connected to indirect environmental conditions such as altitude, rainfall, temperature, wind speed and other climatic factors. Li found that areas with a high incidence also tended to be located at middle and slightly higher altitudes (1600-1800 m), while incidences of brucellosis at other altitudes were relatively low [43]. Low incidence areas tend to be concentrated in the southern part of China. As well as fewer livestock, there may be other influencing factors. These factors require further study.

Our study is not without limitations. Firstly, incidences of human brucellosis were underestimated to some extent, as our data was passively collected by a monitoring system, while surveillance data quality was influenced by comprehensive factors, such as the capacity of local health workers, availability of laboratory diagnostics, levels of awareness about the need to visit doctors and so on, all of which may have affected the study’s accuracy. Secondly, our analysis was based on provincial-level data. Smaller spatial units may provide more location-specific information about the design and implementation stages of public health programs. Thirdly, this study only explored the clustering areas of brucellosis and did not analyze factors affecting the clustering distribution, therefore neglecting the impact of human and natural environment factors. In the future, we will therefore further explore the temporal and spatial distribution of brucellosis as well as the factors affecting its distribution.

In summary, given developments in animal husbandry in China, the prevention and detection of brucellosis urgently need to be strengthened. This spatial clustering study of brucellosis is helpful for identifying high-risk areas for brucellosis, and to a certain extent provides a basis for the decision-making of relevant departments. In terms of the distribution characteristics of brucellosis epidemic in this study, we suggest that the detection of brucellosis in northern China should be further strengthened, and that effective methods for controlling brucellosis should be found in southern China, where incidence of brucellosis is low. Governments at all levels should attach importance to the establishment of joint prevention and control mechanism among high-incidence areas. Relevant departments in different regions should strengthen both prevention and control throughout the year. More resources for prevention and control should be appropriately increased in these areas to curb the spread of brucellosis.


It may be concluded that human brucellosis continues to be a widespread challenge in mainland China, especially in the northwest provinces, which are high risk areas. This study utilized ArcGIS10.3 (ESRI, Redlands) to analyze the temporal and spatial distribution of brucellosis in China, thus contributing to the study of high incidence seasons and areas for brucellosis. However, further research should focus on an analysis of environmental, humanistic and socio-economic factors in order to determine risk factors affecting the occurrence and transmission of brucellosis. Such information has the potential to provide critical guidelines for policy makers to initiate prevention measures and control strategies, aimed at susceptible areas that might be high risk, and to prevent or lessen the incidence of human brucellosis in these areas.

Availability of data and materials

The datasets analyzed during the current study are available from the corresponding author on reasonable request. The URL link of the dataset is:



Centers for Disease Control and Prevention


World Health Organization


Geographic Information System


Local indicators of spatial association


  1. 1.

    Mcdermott J, Grace D, Zinsstag J. Economics of brucellosis impact and control in low-income countries. Rev Sci Tech. 2013;32(1):249–61.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  2. 2.

    Farina F, Fuser R, Rossi M, Scotton PG. Brucellosis outbreak in Treviso province caused by infected cheese from an endemic area. Infez Med. 2008;16(3):154–7.

    CAS  PubMed  PubMed Central  Google Scholar 

  3. 3.

    Park SH, Lee YH, Chu H, Hwang SD, Hwang KJ, Choi HY, et al. Application of the microagglutination test for serologic diagnosis of human brucellosis. Osong Public Health Res Perspect. 2012;3(1):19–23.

    PubMed  PubMed Central  Article  Google Scholar 

  4. 4.

    Luo J, Zeng Z, Song Y, He H. Brucellosis in Takins, China. Emerg Infect Dis. 2012;18(9):1527–9.

    PubMed  PubMed Central  Article  Google Scholar 

  5. 5.

    Ahmadkhani M, Alesheikh A. Space-time analysis of human brucellosis considering environmental factors in Iran. Asian Pac J Trop Dis. 2017;7(5):257–65.

    Article  Google Scholar 

  6. 6.

    Haerry TE, Gehring WJ. Intron of the mouse Hoxa-7 gene contains conserved Homeodomain binding sites that can function as an enhancer element in drosophila. Proc Natl Acad Sci U S A. 1996;93(24):13884–9.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  7. 7.

    Al-Ballaa SR, Al-Balla SR, Al-Aska A, Kambal A, Al-Hedaithy MA. Seasonal variation of culture positive brucellosis at a major teaching hospital. Ann Saudi Med. 1994;14(1):12–5.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  8. 8.

    Yumuk Z, O’Callaghan D. Brucellosis in Turkey-an overview. Int J Infect Dis. 2012;16(4):e228–35.

    PubMed  Article  PubMed Central  Google Scholar 

  9. 9.

    Georgios P. The changing Brucella ecology: novel reservoirs, new threats. Int J Antimicrob Agents. 2010;36(Suppl 1):S8–11.

    Google Scholar 

  10. 10.

    Chen JD, Ke CW, Deng X, Jiang S, Liang W, Ke BX, et al. Brucellosis in Guangdong Province, People’s Republic of China, 2005–2010. Emerg Infect Dis. 2013;19(5):817.

    PubMed  PubMed Central  Article  Google Scholar 

  11. 11.

    Park MY, Lee CS, Choi YS, Park SJ, Lee JS, Lee HB. A sporadic outbreak of human brucellosis in Korea. J Korean Med Sci. 2005;20(6):941–6..

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  12. 12.

    Gwida M, Al Dahouk S, Melzer F, Rösler U, Neubauer H, Tomaso H. Brucellosis-regionally emerging zoonotic disease? Croat Med J. 2010;51(4):289–95.

    PubMed  PubMed Central  Article  Google Scholar 

  13. 13.

    Hotez PJ, Savioli L, Fenwick A. Neglected tropical diseases of the Middle East and North Africa: review of their prevalence, distribution, and opportunities for control. PLoS Negl Trop Dis. 2012;6(2):e1475.

    PubMed  PubMed Central  Article  Google Scholar 

  14. 14.

    Dean AS, Crump L, Greter H, Schelling E, Zinsstag J. Global burden of human brucellosis: a systematic review of disease frequency. PLoS Negl Trop Dis. 2012;6(10):e1865.

    PubMed  PubMed Central  Article  Google Scholar 

  15. 15.

    Sprague LD, Al-Dahouk S, Neubauer H. A review on camel brucellosis: a zoonosis sustained by ignorance and indifference. Pathog Glob Health. 2012;106(3):144–9.

    PubMed  PubMed Central  Article  Google Scholar 

  16. 16.

    Shengjie L, Hang Z, Weiyi X, Marius G, Zhuojie H, Jianxing Y, et al. Changing epidemiology of human brucellosis, China, 1955-2014. Emerg Infect Dis. 2017;23(2):184–94..

    Article  Google Scholar 

  17. 17.

    Zhong Z, Yu S, Wang X, Dong S, Xu J, Wang Y, et al. Human brucellosis in the people’s republic of China during 2005-2010. Int J Infect Dis. 2013;17(5):e289–92.

    PubMed  Article  Google Scholar 

  18. 18.

    Peng G, Wei W, Desheng H. Trends of reported human brucellosis cases in mainland China from 2007 to 2017: an exponential smoothing time series analysis. Environ Health Prev Med. 2018;23(1):23.

    Article  Google Scholar 

  19. 19.

    Seleem MN, Boyle SM, Sriranganathan N. Brucellosis: a re-emerging zoonosis. Vet Microbiol. 2010;140(3–4):392–8.

    PubMed  Article  PubMed Central  Google Scholar 

  20. 20.

    Pappas G, Papadimitriou P, Akritidis N, Christou L, Tsianos EV. The new global map of human brucellosis. Lancet Infect Dis. 2006;6(2):91–9.

    PubMed  Article  PubMed Central  Google Scholar 

  21. 21.

    Cao M, Feng Z, Zhang J, Ma J, Li X. Contextual risk factors for regional distribution of Japanese encephalitis in the People’s Republic of China. Trop Med Int Health. 2010;15(8):918–23.

    PubMed  Article  PubMed Central  Google Scholar 

  22. 22.

    Zhang WY, Guo WD, Sun SH, et al. Human brucellosis, Inner Mongolia, China. Emerg Infect Dis. 2010;16(12):2001–3.

    PubMed  PubMed Central  Article  Google Scholar 

  23. 23.

    Scott LM, Janikas MV. Spatial statistics in ArcGIS// handbook of applied spatial analysis; 2010.

    Google Scholar 

  24. 24.

    Du PL, Wu KS, Fang JY, Zeng Y, Xu ZX, Tang WR, et al. Cervical Cancer Mortality Trends in China, 1991-2013, and Predictions for the Future. Asian Pac J Cancer Prew. 2015;16(15):6391–6.

    Article  Google Scholar 

  25. 25.

    Rao HX, Zhang X, Zhao L, Yu J, Ren W, Zhang XL, et al. Spatial transmission and meteorological determinants of tuberculosis incidence in Qinghai Province, China: a spatial clustering panel analysis. Infect Dis Poverty. 2016;5(1):45.

    PubMed  PubMed Central  Article  Google Scholar 

  26. 26.

    Fang JY, Dong HL, Wu KS, Du PL, Xu ZX. Characteristics and prediction of lung Cancer mortality in China from 1991 to 2013. Asian Pac J Cancer Prev. 2015;16(14):5829–34.

    PubMed  Article  PubMed Central  Google Scholar 

  27. 27.

    Lei S, Yingui C, Wei Z, et al. Study on the variation of arable land use and management countermeasures under rapid urbanization: the application of a gravity model in a regional perspective. Environ Monit Assess. 2019;191(2):120.

    Article  Google Scholar 

  28. 28.

    Xiaopeng Q, Liang W, Barker L, Lekiachvili A, Xingyou Z. Comparison of ArcGIS and SAS Geostatistical analyst to estimate population-weighted monthly temperature for US counties. J Resour Ecol. 2012;3(3):220–9.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  29. 29.

    Tobler W. A computer movie simulating urban growth in the Detroit region. Econ Geogr. 1970;46:234–40.

    Article  Google Scholar 

  30. 30.

    Chen Y. New Approaches for Calculating Moran’s Index of Spatial Autocorrelation. PLoS One. 2013;8(7):e68336.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  31. 31.

    Anselin L. Local indicators of spatial association-LISA. Geogr Anal. 1995;27(2):93–115.

    Article  Google Scholar 

  32. 32.

    Khalid AA, Ali AZ. Spatial autocorrelation of Cancer incidence in Saudi Arabia. Int J Environ Res Public Health. 2013;10(12):7207–28.

    Article  Google Scholar 

  33. 33.

    Zhao X, Huang X, Liu Y. Spatial autocorrelation analysis of Chinese inter-provincial industrial chemical oxygen demand discharge. Int J Environ Res Public Health. 2012;9(6):2031–44.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  34. 34.

    Xia J, Cai S, Zhang H, Lin W, Fan Y, Qiu J, et al. Spatial, temporal, and spatiotemporal analysis of malaria in Hubei Province, China from 2004-2011. Malar J. 2015;14(1):145.

    PubMed  PubMed Central  Article  Google Scholar 

  35. 35.

    Zhang Y, Shen Z, Ma C, Jiang C, Feng C, Shankar N, et al. Cluster of human infections with avian influenza a (H7N9) cases: a temporal and spatial analysis. Int J Environ Res Public Health. 2015;12(1):816–28.

    PubMed  PubMed Central  Article  Google Scholar 

  36. 36.

    Karagiannis I, Mellou K, Gkolfinopoulou K, Dougas G, Theocharopoulos G, Vourvidis D, et al. Outbreak investigation of brucellosis in Thassos, Greece, 2008. Euro Surveill. 2012;17(11):20116.

    PubMed  PubMed Central  Google Scholar 

  37. 37.

    De Massis F, Di Girolamo A, Petrini A, Pizzigallo E, Giovannini A. Correlation between animal and human brucellosis in Italy during the period 1997-2002. Clin Microbiol Infect. 2005;11(8):632–6.

    PubMed  Article  PubMed Central  Google Scholar 

  38. 38.

    Golshani M, Buozari S. A review of brucellosis in Iran: epidemiology, risk factors, diagnosis, control, and prevention. Iran Biomed J. 2017;21(6):349–59.

    PubMed  PubMed Central  Google Scholar 

  39. 39.

    Zhang J, Yin F, Zhang T, Yang C, Zhang X, Feng Z, et al. Spatial analysis on human brucellosis incidence in mainland China: 2004-2010. BMJ Open. 2014;4(4):e004470.

    PubMed  PubMed Central  Article  Google Scholar 

  40. 40.

    Mai HM, Irons PC, Kabir J, Thompson PN. Herd-level risk factors for campylobacter fetus infection, Brucella seropositivity and within-herd seroprevalence of brucellosis in cattle in northern Nigeria. Prev Vet Med. 2013;111(3–4):256–67.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  41. 41.

    Ma JY, Wang H, Zhang XF, Xu LQ, Hu GY, Jiang H, et al. MLVA and MLST typing of Brucella from Qinghai, China. Inf Dis Poverty. 2016;5(1):26.

    Article  Google Scholar 

  42. 42.

    Silva I, Dangolla A, Kulachelvy K. Seroepidemiology of Brucella abortus infection in bovids in Sri Lanka. Prev Vet Med. 2000;46(1):51–9.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  43. 43.

    Li YJ, Li XL, Liang S, Fang LQ, Cao WC. Epidemiological features and risk factors associated with the spatial and temporal distribution of human brucellosis in China. BMC Infect Dis. 2013;13(1):547.

    PubMed  PubMed Central  Article  Google Scholar 

Download references


We are very grateful to referees and editors for their careful reading and valuable comments. We would also like to thank Profs. Yu Zhao, Jiangping Li for their valuable discussion on this work. The authors would like to express their gratitude to EditSprings ( for the expert linguistic services provided.


This work was supported by the Grant from Ningxia Provincial Natural Science Foundation (NO. NZ17187), by the Basic research project of Northwest University for Nationalities (NO. 31920180088).

Author information




PFL participated in the design of the study, YZ drafted the manuscript, JHZ and DFP revised manuscript critically for important intellectual content. ZQG managed the project implementation. All authors read and approved the final version of the manuscript.

Corresponding author

Correspondence to Zhongqin Guo.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Liang, P., Zhao, Y., Zhao, J. et al. The spatiotemporal distribution of human brucellosis in mainland China from 2007-2016. BMC Infect Dis 20, 249 (2020).

Download citation


  • Human brucellosis
  • Mainland China
  • Geographic information system
  • Spatial autocorrelation analysis