Skip to main content
Download PDF
Download ePub
  • Research article
  • Open access
  • Published:

Spatial-temporal distribution of human brucellosis in mainland China from 2004 to 2017 and an analysis of social and environmental factors

Environmental Health and Preventive Medicine volume 25, Article number: 1 (2020) Cite this article

Abstract

Background

This study aimed to describe the changing distribution of human brucellosis between 2004 and 2017 in mainland China and seek scientific evidence of the relationship between socio-economic, environmental, and ecological factors and human brucellosis incidence.

Methods

The annual numbers of brucellosis cases and incidence rates from 31 provinces in mainland China between 2004 and 2017 were obtained from the Data-Center for China Public Health Science. The number of monthly brucellosis cases in 2018 was obtained from the Chinese Center for Disease Control and Prevention. The electronic map of the People’s Republic of China was downloaded from the National Earth System Science Data Sharing Platform. Human population density, gross domestic product (GDP), and an inventory of cattle and sheep at the end of each year from 2004 to 2017 were obtained from the National Bureau of Statistics of China. Annual rainfall data from 31 provinces in the People’s Republic of China from 2004 to 2017 were collected from the China Meteorological Data Service Center. The risk distribution and changing trends of human brucellosis were mapped with ArcGIS. A cluster analysis was employed to identify geographical areas and periods with statistically significant incidence rates. Multivariate linear regression was used to determine possible factors that were significantly correlated with the presence of human brucellosis cases.

Results

Human brucellosis cases have spread throughout the whole country. Human brucellosis cases occurred mostly from March to August and were concentrated from April to July. The inventory of sheep, GDP, and climate were significantly correlated with the presence of brucellosis cases in mainland China.

Conclusions

The geographical expansion of human brucellosis in mainland China was observed, so did the high-incidence clusters between 2004 and 2017. Most of the cases were reported during the early spring to early summer (February–August). Results from the multivariate linear regression suggested that the inventory of sheep, GDP, and climate were significantly associated with the incidence of human brucellosis in mainland China.

Background

Brucellosis, caused by Brucella species, has been recognized as one of the most common worldwide zoonotic diseases [1]. It not only causes substantial direct economic losses in animal production but also presents a great threat to human health [2]. In many countries, there has been a marked decrease in the incidence of human brucellosis, which can be attributed to successful domestic animal brucellosis control or eradication programs [3].

However, brucellosis is endemic in China [4, 5]. Human brucellosis remains a serious public health concern in China, with a 7.8% annual increase in the number of reported cases from 2007 to 2017 and more widespread natural foci [6, 7]. From the perspective of geographical distribution, the affected regions in China gradually expanded from the northern traditional pasturing regions to the agricultural areas and finally to the southern coastal and southwestern areas [7]. From the perspective of occupations at risk, human brucellosis cases have been identified in all classifications of occupations listed in China’s National Noticeable Infectious Disease Reporting (NIDR) system; notably, brucellosis cases have been reported in students, children, and retirees who were not directly exposed to livestock [8]. Low and insufficient awareness and knowledge about brucellosis have also been reported among the Chinese population [9, 10]. The abovementioned changing epidemiology of human brucellosis requires additional focus on the spatial-temporal distribution over time and a better understanding of potential drivers or socio-economic predictors to adjust brucellosis control strategies and allocate health resources.

Thus, the present study was conducted to visualize the annual and monthly human brucellosis incidence rates of each province in the People’s Republic of China from 2004 to 2017 and to explore the impacts of the gross domestic product (GDP), the number of cattle and sheep, and precipitation on the human brucellosis incidence rates in each province from 2004 to 2017, with the aim of providing information for the prevention and control of human brucellosis.

Methods

Data source

In this study, “province” in China was selected as the research unit, and the provinces were divided into two parts, southern and northern provinces, along the Qinling Mountains-Huaihe River line [7]. The map of China showing the province names and the division line is presented in Fig. 1. The annual numbers of brucellosis cases and incidence rates from 31 provinces in China between January 1, 2004, and December 31, 2017, were obtained from the Data-Center for China Public Health Science (http://www.phsciencedata.cn/Share/en/index.jsp). The numbers of monthly brucellosis cases between January 1, 2018, and December 31, 2018, were obtained from the Chinese Center for Disease Control and Prevention (http://www.chinacdc.cn/). An electronic map of China (1:1,000,000) was downloaded from the National Earth System Science Data Sharing Platform (http://www.geodata.cn). The information was unavailable regarding brucellosis incidence in Hong Kong, Macao, and Taiwan; we marked “Data unavailable” when necessary in the figures presented in this study.

Fig. 1
figure 1

Geographical map of mainland China. Provinces were characterized as southern and northern regions according to the Qinling-Huaihe boundary. Areas south of the Qinling-Huaihe boundary are considered southern China, and areas north of the boundary are considered northern China

Full size image

Socio-economic, environmental, and ecological changes have been regarded as triggers for the emergence of emerging infectious diseases (EIDs) [11,12,13]. Pathogens have been indicated to be well adapted to and good at exploiting any social and environmental changes, such as the changes in human demographics or climate, resulting in new chances to spread [14]. In the present study, six factors were selected from the socio-economic (human population density and GDP), environmental (precipitation and climate type), and ecological (the inventory of cattle and the inventory of sheep) perspectives. The human population density in urban areas (persons per km2), GDP, and an inventory of cattle and sheep at the end of each year from 2004 to 2017 were obtained from National Bureau of Statistics of China (http://data.stats.gov.cn/). Annual rainfall data from 173 surface meteorological observation stations in 31 provinces in China (Hong Kong, Macao, and Taiwan were not included) from 2004 to 2017 (average rainfall per year, mm) were collected from China Meteorological Data Service Center (CMDC, http://data.cma.cn/en). The classification of climate type (tropical, subtropical, warm-temperate, mid-temperature, and cold areas) for each province was on the basis of a previous publication [7].

Statistical analysis

The spatial distribution and temporal trend of brucellosis incidence and the relationship between the targeted risk factors and the number of brucellosis cases were explored and visualized using ArcGIS10.6 software (ESRI, Redlands, CA, USA). For the spatial analysis, the data existed in the form of discrete points, and there were no values for the unsampled points. Inversed distance weighted (IDW) interpolation estimates the unknown values by assuming that the weights of the unknown points are greater when they are closer to the known points. Interpolation generates a continuous surface of values. Inversed distance weighted (IDW) interpolation was performed and visualized based on the average incidence between 2004 and 2017 for each province.

Statistically significant high-incidence clusters (hot spots) and low-incidence clusters (cold spots) were identified and visualized with optimized hot spot analysis using ArcGIS 10.6 software. Spatial and temporal statistics were adopted to detect the geographical and temporal clusters between 2004 and 2017 using SaTScan version 9.6 (https://www.satscan.org/), and the results were visualized by ArcGIS10.6. With the assumption that brucellosis incidence followed a discrete Poisson distribution, the maximum spatial cluster size was set to 10% of the population at risk, the minimum temporal cluster size was set to 1 month, and the maximum temporal cluster size was set to 50% of the study period in the Poisson model. The standard Monte Carlo method was employed to determine the P value, and the number of Monte Carlo replications was 999. The most likely cluster, secondary likely clusters, and so on were identified by the parameter log-likelihood ratio (LLR). The LLR compares the difference in the number of cases inside and outside the cluster. The most likely cluster has the greatest value of LLR. Following the most likely cluster, the secondary clusters are outlined by the value of each LLR. A disease cluster was defined when the actual number of cases was greater than the expected number. Assuming that the portion of disease cases inside the scanned cluster was the same as the average portion throughout the whole area, we calculated the value of expected cases as the product of the population of the scanned cluster and the average incidence in the whole country.

Multivariate linear regression was conducted to investigate the association between the inventory of cattle, the inventory of sheep, GDP, human population density, precipitation, climate type, and brucellosis incidence between 2004 and 2017. The model was trained using the annual incidence, number of cattle, number of sheep, GDP, human population density, precipitation, and climate type from 2004 to 2017 for each province. The dataset was shuffled and split into training and testing sets at a ratio of 7 to 3. Before conducting the regression analysis, the values of the GDP, population density, and human brucellosis incidence were log-transformed. An adjusted R2 analysis was selected to evaluate the model performance, and the validation of the model was examined by linearity, a P-P plot, the variance inflation factor (VIF), and a residual versus fitted value scatter plot. A P-P plot was used to check whether the variables in the model were multivariate and normal. The VIF was employed to detect the scale of multicollinearity. A VIF value between 1 and 5 indicated moderate correlation; if the VIF had a value of 10 or more, the correlation was very high, and the model required further manipulation. Multilinear regression was conducted using Python 3.7 software (https://www.python.org/). The correlation coefficients of the chosen variables were obtained through multilinear regression, and we calculated the standardized partial regression coefficients to compare the contribution of each variable to the model.

Results

Brucellosis incidence in China

From 2004 to 2018, there were 530,048 reported cases of human brucellosis in mainland China. Overall, human brucellosis incidence exhibited an increasing trend from 2004 (11,472 cases, 0.88/100,000 persons) and peaked in 2014, with 57,222 cases (4.22/100,000 persons) reported. After that, it decreased slightly, and 40,328 cases were reported in 2018. Most of the brucellosis cases were reported during the period from March to August and concentrated in the period from April to July (Fig. 2). From 2004 to 2017, the years 2004, 2010, and 2016 were selected to compare the quarterly incidence. The second quarter in each year was found to have the highest incidence, and the lowest incidence was found in the fourth quarter (Additional file 1: Figure S1). Other than Tibet, which showed a decreasing trend, human brucellosis incidence in the remaining 30 provinces in China increased over time. Among the 31 provinces, human brucellosis incidence was the highest in the Inner Mongolia Autonomous Region until 2014. The highest incidence was found in Ningxia (43.65/100,000 persons) in 2015 and in Xinjiang (36.56/100,000 persons) in 2016. In 2004, no cases of human brucellosis were reported in Shanghai, Jiangsu, Jiangxi, Hubei, Hunan, Hainan, Chongqing, Sichuan, Guizhou, and Yunnan. After 8 years, human brucellosis cases were reported in all the abovementioned provinces (Fig. 3). Brucellosis hot spots were found in Xinjiang, Inner Mongolia, Heilongjiang, Jilin, Liaoning, Beijing, Tianjin, Hebei, Shanxi, Shaanxi, Ningxia, Gansu, Henan, and Shandong (Additional file 1: Figure S2). The hot spots were all located in northern China, whereas most of the cold spots were in southern China, except Qinghai and Tibet (Additional file 1: Figure S2). The IDW interpolation results also suggested that brucellosis was an endemic disease in northern areas (Fig. 4). With the steadily expanding and increasing trend of epidemics, brucellosis has become a health threat to people living in the southern region.

Fig. 2
figure 2

The number of monthly brucellosis cases from January 2004 to December 2018 in mainland China. The number of brucellosis patients increased between 2004 and 2018. Brucellosis occurred mostly from March to August and was concentrated in April to July in each year

Full size image
Fig. 3
figure 3

Spatial distribution of brucellosis incidence between 2004 and 2017 in mainland China. Since 2013, brucellosis has occurred throughout the country, with an increasing incidence. The incidence in the northern provinces was higher than that in the southern provinces between 2004 and 2017

Full size image
Fig. 4
figure 4

IDW interpolation of average brucellosis incidence between 2004 and 2017 in mainland China. Inverse distance weighting (IDW) estimates the unknown values of the unmeasured points using the weighted average method, which calculates the weight of the unknown points based on their distances to known points. The closer the unknown points are to the known points, the more influence the known will have on the unknown points. The average brucellosis incidence decreased from the north to south between 2004 and 2017. Inner Mongolia was most affected by brucellosis

Full size image

Temporal clusters of human brucellosis

High-risk temporal clusters during 2004 to 2017 were detected from February 2011 to August 2017, which was associated with a 1.91-fold higher risk of contracting human brucellosis than the relative months. High-incidence temporal clusters occurred between February and August of each year; the risk was 1.93- to 3.34-fold higher than that in non-high-incidence temporal clusters (Table 1).

Table 1 Temporal scan of the number of brucellosis cases between 2004 and 2017 in mainland China
Full size table

Seven clusters between 2004 and 2017 were identified using a spatial-temporal scan, including one most likely cluster and six secondary clusters (Fig. 5). Inner Mongolia and Shanxi Province accounted for the most likely cluster, with a radius of 369.72 km; 132,946 cases were observed in the most likely cluster. Most of the cases were reported between April 1, 2008, and March 31, 2015. Between 2004 and 2013, Inner Mongolia and Shanxi Province were both identified as most likely clusters, except in 2010 when Inner Mongolia was the only most likely cluster. An expanding trend of the scale of the most likely clusters was observed, and six provinces (Xinjiang, Tibet, Qinghai, Ningxia, Inner Mongolia, and Gansu) had been identified as most likely clusters since 2014 (Additional file 1: Figure S3, Additional file 1: Table S1).

Fig. 5
figure 5

High-incidence spatial clusters between 2004 and 2017 in mainland China. All of the brucellosis clusters were in the northern provinces, and most of the high-incidence spatial clusters were concentrated around Inner Mongolia and Shanxi Province (most likely cluster), except Tibet (secondary cluster 6)

Full size image

Socio-economic, environmental, and ecological factors that possibly influence brucellosis incidence

The GDP in southeast coastal regions is higher than that in inland regions (Fig. 6). Between 2004 and 2017, Guangdong ($7533), Jiangsu ($6931), and Shandong ($6389) provinces were the top three provinces in terms of GDP (Table 2). Other than Gansu, Guizhou, Hainan, Ningxia, Qinghai, Shanxi, Tibet, and Xinjiang, the GDP in the remaining provinces showed an increasing trend (Additional file 1: Figure S4). The southern provinces with high GDPs had fewer human brucellosis cases between 2004 and 2017 (Fig. 6). The average population densities in Henan (5232.25/km2), Shaanxi (4874.18/km2), and Jiangxi (4367.26/km2) provinces were higher than those in the other provinces in China (Additional file 1: Figure S4, Table 2). The northern provinces are generally characterized as having cold, mid-temperate, or warm-temperate climates. The southern provinces have subtropical and tropical climates (Fig. 7). Tropical, subtropical, and warm-temperate regions had more precipitation than the regions with mid-temperate and cold climates (Additional file 1: Figure S5). The tropical-type, subtropical-type, and high-precipitation areas were associated with lower incidence rates than the other three climate types (Fig. 7). Sheep stocks were larger than cattle stocks (Fig. 8, Table 2). Henan Province had the largest number of cattle stocks, and the Inner Mongolia Autonomous Region had the largest average sheep stock (Additional file 1: Figure S6). Large sheep stocks tended to be associated with more human brucellosis cases than large cattle stocks (Fig. 8).

Fig. 6
figure 6

Density map of the relationship between socio-economic factors (GDP and population density) and average brucellosis cases between 2004 and 2017 in mainland China. In both figures, one dot = 200 human brucellosis cases. The GDP and population density of the southeastern coastal and central provinces were higher than those of the northwest regions. Most of the brucellosis cases were distributed in the northern, northwestern, and northeastern provinces

Full size image
Table 2 Average values of human incidence, the number of sheep, the number of cattle, climate, precipitation, GDP, and population density between 2004 and 2017 in 31 provinces in mainland China
Full size table
Fig. 7
figure 7

Density map of the relationship between environmental factors (precipitation and climate) and average brucellosis cases between 2004 and 2017 in mainland China. In the left figure, one dot = 10. The southern provinces were considered subtropical and tropical regions and had more precipitation than the warm-temperature, mid-temperate, and cold regions located in northern China

Full size image
Fig. 8
figure 8

Density map of the relationship between ecological factors (the number of cattle and the number of sheep) and average brucellosis cases between 2004 and 2017 in mainland China. In both figures, one dot = 10. Cattle stocks were evenly distributed in China; however, sheep stocks (caprine and ovine) were concentrated in the north, which coincided with the reported brucellosis cases

Full size image

By Spearman’s correlation analysis, we found that the correlation between rainfall and cattle and sheep stocks was negative, indicating that where rainfall is high, the number of cattle and sheep might be low. Rainfall positively impacted the GDP. The association between subtropical climate and tropical climate and the number of cattle and sheep was negative. Both mid-temperate and warm-temperate climates had a negative correlation with precipitation, and the mid-temperate climate was also negatively correlated with the GDP, suggesting that the GDP might be relatively low in the mid-temperate areas (Table 3). In subtropical and warm-temperate areas, the regression lines showed that the number of cattle, the number of sheep, precipitation, GDP, and population density barely impacted human brucellosis incidence (Fig. 9). In tropical areas, in addition to population density, which did not play a role in the occurrence of brucellosis, the number of cattle negatively influenced human brucellosis incidence; the number of sheep and GDP had a positive correlation with the incidence between 2004 and 2017 (Fig. 9). Cattle stocks were not associated with the brucellosis incidence in cold areas, but sheep stocks, precipitation, GDP, and population density were associated. Human brucellosis incidence in the mid-temperate areas was higher than that in the areas with the other climate types. Among the selected factors, the number of cattle, the number of sheep, and GDP were positively correlated with human brucellosis incidence, and precipitation and population density were negatively correlated in the mid-temperate areas (Fig. 9).

Table 3 Spearman’s correlation coefficients between the incidence, cattle stocks, sheep stocks, population density, GDP, precipitation, and climate
Full size table
Fig. 9
figure 9

Scatter plots with regression lines representing the association between the number of cattle, the number of sheep, GDP, population density, precipitation, and climate with brucellosis incidence. The blue dots represent those from warm-temperate areas, the yellow squares represent those from mid-temperate areas, the green diamonds represent those from subtropical areas, the red multiplication signs represent those from tropical areas, and the purple upside-down triangles represent those from cold areas. The incidence was higher in the mid-temperate areas than in the other climate areas. In the subtropical areas, the incidence remained at a low level compared with the other climate areas

Full size image

The results from the multivariate linear regression analysis showed that the number of sheep, the GDP, a mid-temperate climate, and a warm-temperate climate were significant independent factors of brucellosis incidence in China between 2004 and 2017 (Table 4). The GDP, the sheep stocks, a mid-temperate climate, and a warm-temperate climate were positively associated with human brucellosis incidence. The value of adjusted R2 was 0.637, meaning that these factors explained 63.7% of the emergence of human brucellosis in China from 2004 to 2017 according to the multivariate linear regression model. The results of the validation of the model (multicollinearity, P-P plot, residual error, linearity check) are compiled in the Additional file (Additional file 1: Table S2, Table S3, and Figure S7).

Table 4 Relationship between socio-economic, environmental, ecological factors and brucellosis incidence in mainland China between 2004 and 2017
Full size table

Discussion

On the basis of the monthly number of nationwide human brucellosis cases during 2004–2018 and the monthly and yearly disease data of each province in mainland China from 2004 to 2017, the present study explored the changing distribution of human brucellosis incidence spatially and temporally and investigated the association between the incidence and certain socio-economic, environmental, and ecological factors. Human brucellosis incidence exhibited an obvious seasonal pattern, gradually increasing in January and peaking in May, after which it gradually decreased. The peak incidence occurred in early spring to early summer, and the lowest incidence usually occurred in winter. This result is similar to the results from other studies in China [15, 16] and other countries, such as Italy [17,18,19]. The identification of seasonal patterns is critical to the success of control and prevention strategies [20]. Sheep breed once or twice a year, mostly in early April and late September; for cattle, the breeding period is from May to July [21]. Since Brucella exists in tissues and body fluids in the placenta, fetal membranes, amniotic fluid, breasts, and lymph nodes [2], during the breeding seasons, the risk of exposure to Brucella increases correspondingly. In the past, brucellosis cases were concentrated in the northern regions of China, and few southern provinces reported any cases [7]. However, since 2014, all provinces have reported brucellosis cases. In the high-incidence cluster scan analysis, we found that the high-incidence clusters expanded from two provinces (Inner Mongolia and Shanxi) in 2004 to six provinces (Xinjiang, Tibet, Qinghai, Ningxia, Inner Mongolia, and Gansu) in 2017 (Additional file 1: Table S1).

With the help of a multivariate linear model, the contributing factors that potentially drove the emergence and spatial expansion of brucellosis in China were explored. Similarly, as previous studies indicated [15, 22], an increased number of sheep were responsible for an increased number of human brucellosis incidence, and our results supported that conclusion. In addition to linear regression, we conducted ridge regression and Poisson regression to compare the values of the coefficients with those calculated by linear regression and verify the relationship between the potential drivers and human brucellosis incidence in China. The results were similar (Additional file 1: Table S4 and Table S5). Thus, we concluded that the number of sheep was positively associated with the increase in human brucellosis incidence. Our results also revealed that cattle’s involvement in human brucellosis incidence was not significant. In China, B. melitensis is the predominant strain. Large-scale serological testing of cattle and sheep found that the positive rate of sheep sera was higher than that of cattle sera [23]. In addition, in high-incidence areas, sheep were the predominant livestock. To cater to the domestic mass demand for meat, the stock of sheep has increased by six times [24]. As sheep are a major influencing factor in the spread of brucellosis, surveillance and corresponding prevention and control measures for sheep populations should be a priority. It has been indicated that changing socio-economic conditions, such as increasing population densities in urban areas and increasing GDPs, may be correlated with the possibility of infection emergence and transmission [13]. Countries with relatively high GDP are normally brucellosis-free [25]. Densely populated areas and areas with poor sanitation and hygiene contribute to the spread of food-borne and soil-borne diseases [26]. With the rapid acceleration of urbanization in developing countries, rural infectious diseases have emerged, posing a significant threat to urban residents [26]. According to the China Statistical Yearbook, the population density in urban areas of China presented an overall increasing trend. In our multivariate linear regression model, no statistically significant correlations between population density and human brucellosis incidence between 2004 and 2017 were found. Urbanization in China does not affect human brucellosis incidence. The Spearman coefficient between the GDP and incidence was − 0.06; however, the GDP was significantly positively correlated with brucellosis incidence in the multivariate analysis. In China, brucellosis cases concentrated in provinces defined as northern areas that have generally lower GDPs and were also categorized in our research as mid-temperate or warm-temperate areas. This positive association was more likely to exist in those provinces that were economically less advantaged in China between 2004 and 2017 than in provinces with high economic advantages. Further analysis should focus on the conditions in which the GDP is positively associated with human brucellosis incidence. We found that among the five climate zones, the mid-temperate climate and warm-temperate climate significantly influenced the emergence of brucellosis in China, with cold climate as the reference climate. There have been some hypotheses regarding rainfall and climate changes as environmental drivers for vector-borne diseases [27, 28]. The present study did not find such a pattern among precipitation and human brucellosis events. This could be because in the model, the variable precipitation was not enough to explain the whole picture, unlike climate [11]. Thus, further study on brucellosis epidemiology regarding ecologic factors should focus on the characteristics and differences in and between each climate type to identify the determinants.

Globally, zoonoses account for 60.3% of EIDs [11], and livestock account for up to 37% of the agricultural GDP [29]. Over the past 14 years, China’s economy has entered a high-speed growth phase, with the per capita GDP increasing from 1828 US dollars in 2004 to 8828 US dollars in 2017 [30]. The increase in the demand for high living standards stimulates the demand for meat, which consequently causes an increase in livestock farming and meat importation [31]; thus, the risk of exposure to animals, either live or slaughtered, is increasing. In China, the amount of food and animal imports is approximately 7–8 times higher than the amount of exports, and importation is increasing at a rate of 15% per year [31]. International trade could facilitate the spread and introduction of new types of Brucella pathogen from other countries, which might increase the difficulty of control and prevention. In light of the expanding scale and increasing incidence of brucellosis in China, we found that the human brucellosis incidence increased more rapidly in certain provinces, such as Heilongjiang, Shandong, Henan, Gansu, Ningxia, Zhejiang, and Guangdong provinces. The driving and influencing factors of the spread should be elucidated. Our results suggested that direct contact with the infectious livestock remains the primary source of human brucellosis. The reasons behind the increasing incidence in the southern provinces in China could be travel and the increasing demand for mutton and beef, which has led to an increase in the transportation of cattle, sheep, and other animal products from the northern pastoral areas, such as Inner Mongolia [32]. It is critical to strengthen the management and supervision of the agricultural market and the international and domestic transportation of livestock to prevent contaminated meat and dairy products from entering the market [3].

Occupational groups, such as slaughterhouse workers, meat-packing employees, veterinarians, and herdsmen, still account for the majority of brucellosis patients. A study in Pakistan showed that inadequate understanding of brucellosis can lead to risky behaviors that could cause infection [10, 33]. With the development of economic globalization, the amount of air travel passengers increased by five times from 2004 to 2017 [34]. There have been several reports about brucellosis infection among returning travelers [35, 36]. Health education for travelers should also be a priority, especially for travelers visiting high-risk regions during high-risk periods. The focus of health education should be avoiding the consumption of undercooked meat and unpasteurized dairy products and maintaining a safe distance from livestock. For medical practitioners, especially those who work in non-brucellosis-endemic areas, if patients complain of fever, sweating, joint pain, and malaise, they should be asked about their travel history and exposure history to avoid a misdiagnosis or late diagnosis of human brucellosis.

“One Health” has been defined as a collaborative, multifaceted, and multidisciplinary approach to address local, regional, national, and worldwide health problems. A strength of the present study was that it was conducted at the livestock-human-ecological interface. We mapped the changing epidemiology of brucellosis between 2004 and 2017 in 31 provinces of China, detected “hot spot” areas and time frames, and quantified the impact of representative socio-economic, environmental, and ecological factors on brucellosis incidence.

There are some limitations to our study. First, data regarding human brucellosis incidence were collected by passive monitoring due to the national policy. Most brucellosis cases occurred in farmers, who often do not have adequate access to a timely diagnosis and effective treatment. Thus, the brucellosis incidence identified in this analysis might underestimate the real extent of the brucellosis situation in China. Second, we macroscopically described the changing trend of human brucellosis incidence throughout the whole country using province as the study unit. Thus, the association between the changing incidence and the potentially driving factors might not accurately describe the local relationship. Future studies should divide the whole country into several categories based on socio-economic or ecological characteristics, such as geographical locations, climate conditions, vegetation coverage, vegetation type, people’s living standards, prices of meat, or meat consumption. In addition, other machine learning models can be used to analyze the changing epidemiology of human brucellosis, explore the influencing factors and their interactions, examine their associations, and further predict the future trends of human brucellosis incidence under these classifications.

Conclusions

In summary, the geographical extent of human brucellosis in China has been expanding with an increasing incidence in each province. The human brucellosis incidence and the high-incidence clusters were concentrated in the northern areas but were also observed to be expanding. Most cases were reported during the early spring to early summer (February–August). The inventory of sheep, GDP, and climate type were significantly associated with the incidence of human brucellosis in China between 2004 and 2017.

Availability of data and materials

The datasets analyzed during the current study are available from the corresponding author on reasonable request.

Abbreviations

CMDC:

China Meteorological Data Service Center

EIDs:

Emerging infectious diseases

GDP:

Gross domestic product

IDW:

Inversed distance weighted

LLR:

Log-likelihood ratio

NIDR:

National Noticeable Infectious Disease Reporting

VIF:

Variance inflation factor

References

  1. Godfroid J. Brucellosis in livestock and wildlife: zoonotic diseases without pandemic potential in need of innovative one health approaches. Arch Public Health. 2017;75:34.

    Article  Google Scholar 

  2. Corbel MJ. Brucellosis: an overview. Emerg Infect Dis. 1997;3(2):213–21.

    Article  CAS  Google Scholar 

  3. Zhang N, Huang D, Wu W, Liu J, Liang F, Zhou B, et al. Animal brucellosis control or eradication programs worldwide: a systematic review of experiences and lessons learned. Prev Vet Med. 2018;160:105–15.

    Article  Google Scholar 

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

    Article  Google Scholar 

  5. Zhang X, Wang Z, Mu G, Wang T. Brucellosis control in northeast China: a long way to go. Public Health. 2015;129(8):1132–4.

    Article  CAS  Google Scholar 

  6. Guan P, Wu W, Huang D. 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 

  7. Lai S, Zhou H, Xiong W, Gilbert M, Huang Z, Yu J, et al. Changing epidemiology of human brucellosis, China, 1955-2014. Emerg Infect Dis. 2017;23(2):184–94.

    Article  Google Scholar 

  8. Cui BY, Jiang H. Surveillance data of brucellosis in China, 2005-2016. Dis Surveill. 2018;33:188–92 (Article in Chinese).

    Google Scholar 

  9. Zeng H, Wang Y, Sun X, Liu P, Xu Q, Huang D, et al. Status and influencing factors of farmers’ private investment in the prevention and control of sheep brucellosis in China: a cross-sectional study. PLoS Negl Trop Dis. 2019;13(3):e0007285.

    Article  Google Scholar 

  10. Zhang N, Zhou H, Huang DS, Guan P. Brucellosis awareness and knowledge in communities worldwide: a systematic review and meta-analysis of 79 observational studies. PLoS Negl Trop Dis. 2019;13(5):e0007366.

    Article  Google Scholar 

  11. Jones KE, Patel NG, Levy MA, Storeygard A, Balk D, Gittleman JL, et al. Global trends in emerging infectious diseases. Nature. 2008;451(7181):990–3.

    Article  CAS  Google Scholar 

  12. Morse SS. Factors in the emergence of infectious diseases. Emerg Infect Dis. 1995;1(1):7–15.

    Article  CAS  Google Scholar 

  13. Weiss RA, McMichael AJ. Social and environmental risk factors in the emergence of infectious diseases. Nat Med. 2004;10(12 Suppl):S70–6.

    Article  CAS  Google Scholar 

  14. Woolhouse MEJ, Gowtage-Sequeria S. Host range and emerging and reemerging pathogens. Emerg Infect Dis. 2005;11(12):1842–7.

    Article  Google Scholar 

  15. Li Y, Li X, Liang S, Fang L, Cao W. Epidemiological features and risk factors associated with the spatial and temporal distribution of human brucellosis in China. BMC Infect Dis. 2013;13:547.

    Article  Google Scholar 

  16. Wang Y, Xu C, Zhang S, Wang Z, Zhu Y, Yuan J. Temporal trends analysis of human brucellosis incidence in mainland China from 2004 to 2018. Sci Rep. 2018;8(1):15901.

    Article  Google Scholar 

  17. Mancini FR, Bella A, Graziani C, Marianelli C. Trends of human brucellosis in Italy, 1998-2010. Epidemiol Infect. 2014;142(6):1188–95.

    Article  CAS  Google Scholar 

  18. Nematollahi S, Ayubi E, Karami M, Khazaei S, Shojaeian M, Zamani R, et al. Epidemiological characteristics of human brucellosis in Hamadan Province during 2009-2015: results from the National Notifiable Diseases Surveillance System. Int J Infect Dis. 2017;61:56–61.

    Article  Google Scholar 

  19. Pakzad R, Pakzad I, Sa S, Reza M. Spatiotemporal analysis of brucellosis incidence in Iran from 2011 to 2014 using GIS. Int J Infect Dis. 2018;67:129–36.

    Article  Google Scholar 

  20. He F, Hu ZJ, Zhang WC, Cai L, Cai GX, Aoyagi K. Construction and evaluation of two computational models for predicting the incidence of influenza in Nagasaki Prefecture, Japan. Sci Rep. 2017;7(1):7192.

    Article  Google Scholar 

  21. Lou P, Wang L, Zhang X, Xu J, Wang K. Modelling seasonal brucellosis epidemics in Bayingolin Mongol Autonomous Prefecture of Xinjiang, China, 2010-2014. Biomed Res Int. 2016;2016:5103718.

    PubMed  PubMed Central  Google Scholar 

  22. Zhou L, Fan M, Hou Q, Jin Z, Sun X. Transmission dynamics and optimal control of brucellosis in Inner Mongolia of China. Math Biosci Eng. 2018;15(2):543–67.

    Article  Google Scholar 

  23. Man T, Wang D, Cui B, Wang Y, Ding F, Li T, et al. Analysis on surveillance data of brucellosis in China, 2009. Dis Surveill. 2010;25:944–6 (Article in Chinese).

    Google Scholar 

  24. National Development and Reform Commission of China. National beef and mutton production development plan (2013-2020). 2013. Available at: http://www.sdpc.gov.cn/zcfb/zcfbghwb/201402/P020140221362074574291.pdf. Accessed 10 Feb 2019.

    Google Scholar 

  25. Cárdenas L, Awada L, Tizzani P, Cáceres P, Casal J. Characterization and evolution of countries affected by bovine brucellosis (1996-2014). Transbound Emerg Dis. 2019;66(3):1280–90.

    Article  Google Scholar 

  26. Neiderud CJ. How urbanization affects the epidemiology of emerging infectious diseases. Infect Ecol Epidemiol. 2015;5(1):27060.

    PubMed  Google Scholar 

  27. Guernier V, Hochberg ME, Guégan JF. Ecology drives the worldwide distribution of human diseases. PLoS Biol. 2004;2(6):e141.

    Article  Google Scholar 

  28. Patz JA, Campbell-Lendrum D, Holloway T, Foley JA. Impact of regional climate change on human health. Nature. 2005;438(7066):310–7.

    Article  CAS  Google Scholar 

  29. Wiethoelter AK, Beltrán-Alcrudo D, Kock R, Mor SM. Global trends in infectious diseases at the wildlife-livestock interface. Proc Natl Acad Sci U S A. 2015;112(31):9662–7.

    Article  CAS  Google Scholar 

  30. National Bureau of Statistics of China. China Statistical Yearbook. 2018. Available at: http://data.stats.gov.cn/english/index.htm. Accessed 10 Feb 2019

    Google Scholar 

  31. Wu T, Perrings C, Kinzig A, Collins JP, Minteer BA, Daszak P. Economic growth, urbanization, globalization, and the risks of emerging infectious diseases in China: a review. Ambio. 2017;46(1):18–29.

    Article  CAS  Google Scholar 

  32. Zhang J, Sun GQ, Sun XD, Hou Q, Li MT, Huang BX, et al. Prediction and control of brucellosis transmission of dairy cattle in Zhejiang Province, China. PLoS One. 2014;9(11):e108592.

    Article  Google Scholar 

  33. Arachchige K, Hemali C, Kothalawala A, Makita K, Kothalawala H. Knowledge, attitudes, and practices (KAP) related to brucellosis and factors affecting knowledge sharing on animal diseases: a cross-sectional survey in the dry zone of Sri Lanka. Trop Anim Health Prod. 2018;50(5):983–9.

    Article  Google Scholar 

  34. World Bank. World development indicators. Washington, D.C: World Bank; 2018. Available at: https://datacatalog.worldbank.org/dataset/world-development-indicators. Accessed 10 Feb 2019

    Google Scholar 

  35. Di Pierdomenico A, Borgia SM, Richardson D, Baqi M. Brucellosis in a returned traveller. CMAJ. 2011;183(10):E690–2.

    Article  Google Scholar 

  36. Rhodes HM, Williams DN, Hansen GT. Invasive human brucellosis infection in travelers to and immigrants from the Horn of Africa related to the consumption of raw camel milk. Travel Med Infect Dis. 2016;14(3):255–60.

    Article  Google Scholar 

Download references

Acknowledgements

This study is based in part on the information from the documents provided by the National Bureau of Statistic the People’s Republic of China, Data-Center of China Public Health, and Chinese Center for Disease Control and Prevention.

Funding

This study was supported by the grant from the National Natural Science Foundation of China (grant nos.71573275 and 71974199) and Science Foundation of Liaoning Provincial Department of Education (LFWK201719).

Author information

Authors and Affiliations

  1. Department of Epidemiology, School of Public Health, China Medical University, Shenyang, 110122, China

    Cheng Peng, Yan-Jun Li, De-Sheng Huang & Peng Guan

  2. Department of Mathematics, School of Fundamental Sciences, China Medical University, Shenyang, 110122, China

    De-Sheng Huang

Authors
  1. Cheng Peng
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Yan-Jun Li
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. De-Sheng Huang
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Peng Guan
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

CP collected and analyzed the statistics, and drafted the manuscript. PG provided the overall leadership to the study and participated in the analysis. DSH and PG secured the funding. YJL and DSH participated in the design and statistical analysis. All authors reviewed and approved the final version of the manuscript.

Corresponding author

Correspondence to Peng Guan.

Ethics declarations

Ethics approval and consent to participate

The study protocol was approved by the Research Institutional Review Board of China Medical University, and it was exempt from institutional review board assessment due to the fact that it was conducted on the basis of the nationwide surveillance of infectious diseases in China and all the data were publicly open, supplied, and analyzed in an anonymous format, and without any access to the personal identifying information.

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.

Supplementary information

Additional file 1: Table

S1. Spatio-temporal scan analysis of high-incidence clusters between 2004 and 2017 in mainland China. Table S2. VIFs of the number of cattle, the number of sheep, GDP, population density, rainfall and climate in the multivariate linear regression model. Table S3. VIF of the number of cattle, the number of sheep, GDP, population density, rainfall, Tropical climate, Mid-temperate climate and Warm-temperate climate in the multivariate linear regression model. Table S4. Coefficients of the number of sheep, GDP, population density, rainfall and climate in the ridge regression model. Table S5. Coefficients of the number of sheep, GDP, population density, rainfall and climate in the Poisson regression model. Figure S1. Spatial distribution of quarterly brucellosis incidence (100,000 persons) in 2004, 2010 and 2016 in mainland China. Figure S2. Hot and cold spots of brucellosis incidence between 2004 and 2017 in mainland China. Figure S3. Spatio-temporal distribution of high-incidence clusters between 2004 and 2017 in mainland China. Figure S4. GDP (left) and population density (right) provincial distribution in mainland China between 2004 and 2017. Figure S5. Climate and Precipitation distribution in provinces of mainland China between 2004 and 2017. Figure S6. Cattle stocks (left) and Sheep stocks (right) provincial distribution in mainland China between 2004 and 2017. Figure S7. Validation of multivariate linear model: (a) P-P plot, (b) Residual error and (c) linearity check.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Peng, C., Li, YJ., Huang, DS. et al. Spatial-temporal distribution of human brucellosis in mainland China from 2004 to 2017 and an analysis of social and environmental factors. Environ Health Prev Med 25, 1 (2020). https://doi.org/10.1186/s12199-019-0839-z

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12199-019-0839-z

Keywords

Environmental Health and Preventive Medicine

ISSN: 1347-4715

Contact us