Skip to main content

Epidemiological characteristics and spatio-temporal analysis of brucellosis in Shandong province, 2015–2021

Abstract

Background

Brucellosis is one of the major public health problems in China, it not only causes huge economic losses to the society, but also threatens the human’s physical and mental health. The reported cases of brucellosis in Shandong province were at a high level, therefore, it is necessary for us to understand the epidemic characteristics and distribution trend of Brucellosis in Shandong province. This study aims to describe the epidemiological characteristics and spatial clustering characteristics of brucellosis in Shandong Province, provide a reference for the scientific prevention and control.

Methods

Human brucellosis data in Shandong province from 2015 to 2021 were obtained from the China Information System for Disease Control and Prevention, the data were analyzed by descriptive epidemiological methods, spatial autocorrelation analysis and spatial-temporal cluster analysis methods use ArcGIS and SaTScan software, the results were presented in ArcMap.

Results

A total of 22,251 human cases of brucellosis were reported, the annual incidence ranged between 2.41/100,000 and 4.07/100,000 from 2015 to 2021 in Shandong province, incidence has been decreasing year by year, while there was a significant increase in 2021. The distribution of brucellosis was of a seasonal trend, mainly concentrating during March to August. The age of the cases was mainly concentrated in the 30–74 age ranges, the average annual incidence rate was significantly higher in males than in females. The spatial analysis showed that the epidemics were mainly concentrated in the north and southwest. For the spatial autocorrelation analysis, a high global autocorrelation was observed at the county level, and the high–high clusters mainly distributed in the north and southwest region. For the spatio-temporal scanning, the most likely cluster areas mainly distributed in the north area, and then gradually moved southward, and the radius of clustered narrowed.

Conclusions

Human brucellosis remains a common challenge, particularly in northern region in spring and summer. More disease prevention and control measures should be taken in high-risk populations, and such higher-risk susceptible areas to reduce the incidence of brucellosis and ensure the health of the people.

Peer Review reports

Background

Brucellosis is the most important zoonotic disease caused by Brucella species comprising Gram negative, facultative, intracellular pathogens, [1, 2]. Domestic animals such as cattle, goats, sheep, pigs, camel, buffalo and dogs serve as a reservoir hosts [3]. The disease is transmitted to humans by direct/indirect contact with infected animals or through the consumption of raw meat and dairy products. The main transmission routes are digestive tract, skin, and mucosal and respiratory tract contact with blood body fluids and aerosols [4]. The impact of brucellosis on human health is a major issue, it cause fever, nausea, muscular pain, abdominal pain, sweating, weakness, decreased appetite, weight loss, and liver inflammation [1]. Brucellosis not only causes huge economic losses to the society, but also threatens the human’s physical and mental health [5]. Human brucellosis is one of the major public health problems in China.

Brucellosis is prevalent worldwide, Cases of brucellosis have been reported in 170 countries and regions. There are about 50,000 to 6 million people with more than 500,000 human cases reported annually worldwide, the Middle East, Asia, the Mediterranean Basin, Africa, Central America and the Caribbean being the main endemic regions [6, 7]. The evolution of the “global village” through international tourism. new endemic foci have emerged [8], So far, China has 30 provinces (municipalities, autonomous regions) have different degrees of brucellosis epidemic, there are now 300,000 to 500,000 patients [9]. 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 [10]. The reported cases of brucellosis in Shandong province are at a high level, with the number of reported cases averaging more than 2700 cases per year in the past three years. The situation of brucellosis prevention and control is grim. This study aims to describe the epidemiological characteristics and spatial clustering characteristics of brucellosis in Shandong Province, provide a reference for the scientific prevention and control of brucellosis.

Methods

Research data

The data of human brucellosis cases in Shandong province between 2015 and 2021were obtained from the China Information System for Disease Control and Prevention. The main contents of the database of human brucellosis include the number of cases, incidence rate, number of deaths and mortality rate in multiple dimensions by the region, age group, gender and occupation. The diagnostic criteria were based on the national health industry standard. The demographic information were obtained from the Shandong Statistical Bureau. Geographic space information was acquired from the National Geographic Information Public Service Platform.

Statistical analysis

The SPSS software (version 23.0, IBM company, NewYork, USA) was used for the descriptive epidemiological analysis. χ2 test was used for comparison of count data, and the test level was α = 0.05. Spatial mapping and spatial autocorrelation analysis were based on Arcgis software(Version, 10.6, Environmental Systems Research Institute, Redlands, USA). SaTScan software(version 9.5, Boston, MA, USA) was used for spatial-temporal statistical.

Spatial autocorrelation analysis

As a spatial statistical method, global spatial autocorrelation and local spatial autocorrelation are used to describe the relationship between study areas and measure the degree of clustered or dispersion [11,12,13]. The Moran’s I index is used to measure overall spatial autocorrelation and spatial distribution of the study areas while the local one can be further used to reflects the local spatial autocorrelation and the specific clustering areas [14]. The value of Moran’s I ranges between [− 1, 1], I > 0 indicates a positive spatial correlation, while I < 0 indicates a negative spatial correlation, if the I value is close to 0, indicates that the cases are randomly distributed in space, no spatial correlation exists [15, 16]. The larger the absolute value of I, the stronger the correlation. When Z > 1.96, P < 0.05 was considered to be statistically significant and a spatial autocorrelation existed [17, 18]. LISA is used to reflect a geographical phenomenon on a regional unit or the degree of correlation between an attribute value and the same geographical phenomenon or attribute value on the adjacent unit [19]. The spatial correlation patterns obtained from the local Moran’s I index can be shown by LISA, which are classified into four types ,low-high cluster (L-H, which indicated that the low cluster areas were surrounded by high cluster areas), high–low cluster (H-L, which indicated that the high cluster areas were surrounded by other low cluster areas), low-low cluster (L-L, which indicated the cold spot), and high-high cluster (H-H, which indicated the hot spot) [20, 21]. The spatial autocorrelation analysis was conducted by ArcGIS software.

Spatio-temporal cluster analysis

We used the spatio-temporal scan statistics to detect the center and radius of the aggregation area [22, 23], and verify whether the time and geographic clustering of human brucellosis was caused by random variation or not [24]. The basic principles of spatio-temporal scan is based on a discrete Poisson model [25]. The theoretical incidence number of each scanning window is calculated and compared with the actual incidence number to construct the log likelihood ratio (LLR) [26], By calculating the LLR of the spatial unit attributes within and outside the dynamic window area under different centers and radiuses, SaTScan makes the statistical inference and explores the maximum possible clustering area using Monte Carlo for the statistical significance evaluation [27].For each possible spatial-temporal clustering area, when P < 0.05, the larger the log-likelihood ratio (LLR) value, the more the likelihood that the area covered by the scanning dynamic window represents the clustering area [28, 29]. Finally, the window with the largest LLR value is selected as the maximum possible clustering area, which represents this high-risk region [30], while other statistically significant windows were the secondary possible clustering area [31]. In spatio-temporal scan analysis, the spatial unit was set as county, the temporal unit was set as year. Circular moving windows were set to scan the study area. The default size settings for the window and time are usually set to 50% [32], The number of Monte Carlo randomization was set to 999.SatScan software was used for spatial-temporal statistical, and Arc GIS software was used for visual presentation of the scanning results.

Results

Demographic characteristics

From 2015 to 2021, the average annual incidence rate was 4.42/100,000 in males and 1.84/100,000 in females, which were significantly higher in males than in females (χ2 = 3142.31, P<0.01). The male-to-female ratio ranging from2.43:1 to 2.65:1. The incidence is generally higher in the age range of 30–74 years ranging from1.57/100,000 to 9.88/100,000, we found that the occupation distribution is mainly farmers, from 2015 to 2021, a total of 18,975 cases were reported, accounting for 85.28% of the total cases. (Table 1).

Table 1 Demographic characteristics of Brucellosis in Shandong province, 2015–2021

Temporal characteristics

A total of 22,251 cases of brucellosis were reported in Shandong province from 2015 to 2021, and the annual reported incidence ranged between 2.41/100,000 and 4.07/100,000. Since 2015, the incidence of brucellosis in Shandong province has firstly increased and then decreased, and reached the peak in 2016, with reported incidence 4.07/100,000 From 2017 to 2020, brucellosis has been decreasing year by year, while there was a significant increase in 2021(Fig. 1). Brucellosis cases were reported every month throughout the year, with clear seasonal distribution characteristics. The peak of brucellosis incidence was concentrated from March to August, accounting for 66.08% of the total cases (Fig. 2).

Fig. 1
figure 1

Reported cases and incidence rate of Brucellosis in Shandong province, 2015–2021

Fig. 2
figure 2

Monthly distribution of Brucellosis cases in Shandong province, 2015–2021

Spatial characteristics

Brucellosis cases were reported in 16 cities of Shandong province from 2015 to 2021, and the number of counties (districts) affected by the epidemic was 123, 122, 122, 140, 137, 140 and 140 respectively. The top five counties (districts) with average annual incidence were Lijin County (29.32/100 000), Kenli County (12.39/100 000), Sishui County (10.92/100 000), Wudi County (9.85/100 000) and Zhanhua County (9.63/100 000), mainly located in the north and southwest of Shandong Province (Fig. 3). In space, there is a gradual narrowing of the central epidemic area. In general, the incidence rates in most areas have been decreasing year by year, especially in the central area (Fig. 4).

Fig. 3
figure 3

Regional distribution of annual incidence of Brucellosis in Shandong Province, 2015–2021

Fig. 4
figure 4

The yearly incidence of Brucellosis at the county level in Shandong province, 2015–2021

Spatial autocorrelation analysis

The global spatial autocorrelation analysis found that the annual global Moran’s I values ranged from 0.205 to 1.547, while the Z values were greater than 1.96 (all p < 0.001), which indicates that the incidence of brucellosis between 2015 and 2021 had a positive spatial autocorrelation and obvious spatial clustered distribution at the county level, the highest clustering degree was in 2020 (Moran’s I = 1.547) and the lowest in 2016 (Moran’s I = 0. 205) (Table 2). Local autocorrelation analysis results were shown in Fig. 5. The LISA clustered distribution map showed that there were15, 12,10, 8, 6,8 and9 counties (districts) in the high-high clustered area from 2015 to 2021, and the “high-high” clustered area was mainly concentrated in the northern part of the province. Located in Distributed in Hekou District, Lijin County, Kenli District, Wudi County, Yangxin County, Huimin County, Zhanhua County, Bincheng District, Qingyun County, Shanghe County, Pingyi County, Shouguang city, Sishui County, Qufu City, Xintai City, Ningyang County and other regions. Hekou District, Lijin County and Zhanhua County are the “high-high” gathering areas over the years. Over time, the number of counties (districts) covered by high-high clustered areas showed a downward trend. In recent years, Pingyi County of Linyi City, Ningyang County of Taian city and Qufu City of Jining became the new “high-high” clustered areas. The “low-low” clustered area area became smaller year by year and its distribution changed to concentrated in the eastern region.

Table 2 Results of the spatial autocorrelation test on Brucellosis cases in Shandong province, 2015–2021
Fig. 5
figure 5

Results of local spatial autocorrelation analysis of Brucellosis at the county level in Shandong province, 2015–2021

Spatio-temporal cluster analysis

The results of the spatio-temporal cluster analysis for the human brucellosis cases in Shandong province from 2015 to 2021 included one most likely cluster, one secondary cluster and one 2nd secondary cluster. The most likely cluster was located in the northern area, including 37 counties and districts, of which the cluster center was (37.88 N,118.52E) and the cluster radius was 167.91 km. and the time frame was from January 2015 to December 2017, with a relative risk (RR) value of 2.52 (P < 0.0001). The secondary clusters were relatively smaller, mainly distributed in the central area, including Qufu City, Sishui County, of which the cluster center was (35.66 N, 117.25E) and the cluster radius was25.56 km. and the cluster time frames mainly ranged from January 2019 to December 2021 ,with a relative risk (RR) value of 4.17 (P < 0.0001). The 2nd secondary clusters were mainly distributed in the southwest area of Shandong province, including Dingtao District, Mudan district, Cao County, and the cluster time frames mainly ranged from January 2016 to December 2018 with a relative risk (RR) value of 2.59 (P < 0.0001) (Table 3; Fig. 6).

Table 3 Spatio-temporal clustering of human Brucellosis in Shandong province, 2015–2021
Fig. 6
figure 6

Spatio-temporal clustering of brucellosis cases in Shandong province from 2015to 2021. (The red circles represent the clustering areas, such that the largest represents the most likely cluster area.)

Discussion

Brucellosis is one of the most widespread zoonoses worldwide [7, 33], In China brucellosis has been increasingly causing huge economic loss, and it has been a public health problem in recent years [4]. Public health researchers are usually interested in using data visualization methods to describe the distribution of disease. It can guide public health officers and medical providers to fast response and make decision about containment, management, and prevention of diseases, to improve the surveillance and control strategies [34, 35]. In this study, the incidence of brucellosis at the county level in Shandong province from 2015 to 2021 was used to discuss the epidemiological characteristics and investigate its spatial and temporal distribution rules and possible clustered areas. Better display of incidence and spatial clustered through map visualization, this study provides the theoretical support for the prevention and control of human brucellosis.

From 2005 to 2014, the incidence of brucellosis in China was on the rise, with the reported incidence rising from 1.40/100,000 to 4.22/100,000, and then the incidence was on the decline from 2015 to 2017 [36]. The brucellosis epidemic in Shandong province is slightly different from that in the whole country, showing an upward trend and then a downward trend since 2015 to 2020, and an upward trend in 2021. It peaked in 2016 with incidence rate of 4.07 per 100,000. Shandong province issued the corresponding strengthening of prevention and control policies, to suppress brucellosis epidemic rising trend played a positive role. The increase in 2021 May be associated with abundant rain in 2021, which is more conducive to the propagation and spread of pathogens in the high-incidence season, the spread of infectious diseases with natural foci is affected by rainfall [37].

Our study found that the higher incidence was concentrated in the 30–74 age ranges, most of the people working with livestock are middle-aged, within this age range. In rural Chinese families, most middle-aged and elderly people aged between 40 and 65 are raising livestock [38], so the incidence of brucellosis among this age group is higher. Drinking unpasteurized goat milk or cow’s milk can also lead to brucellosis infection [39, 40], infection in other age groups may be associated with drinking raw milk or other means of exposure to infected livestock or their products. More males than females are affected, because of there are more men in jobs that have more contact with livestock, such as breeding and slaughtering.

The results of our temporal cluster analysis show that there is a distinct seasonality in human brucellosis. The main cluster time of the reported cases is from March to August, the peak incidence occurred in early spring to early summer, and the lowest incidence usually occurred in winter. This result is consistent with the results from other studies in China [41,42,43] and some other Asian countries [44, 45]. The reason might be that these months represent the time of lambing in the agricultural and pastoral areas. Since Brucella exists in tissues and body fluids in the placenta, fetal membranes,

amniotic fluid, breasts, and lymph nodes. during the breeding seasons, the risk of exposure to Brucella increases correspondingly [41]. It may also be due to warm temperatures suitable for the transmission of zoonosis, according to some reported research findings, the factors of temperature, sunshine, wind, speed, altitude and rainfall will affect the introduction of brucellosis [46, 47]. Identifying the peak of seasonal incidence can be used to strengthen prevention and control strategies before the incidence increases.

The results of spatial autocorrelation analysis showed that the incidence of brucellosis in Shandong Province from 2015 to 2021 had a certain spatial clustered with positive spatial autocorrelation at the county level. The Moran’s I value of the global spatial autocorrelation index decreased first and then increased, indicating that the degree of global spatial clustered decreased first and then increased, the highest clustering degree was in 2020.Local autocorrelation analysis results were showed that the “high-high” clustered area was mainly concentrated in the northern part of the province. Located in distributed in Hekou District, Lijin County, Kenli District, Wudi County, Yangxin County, Huimin County, Zhanhua County, Bincheng District, Qingyun County, Shanghe County, Pingyi County, Shouguang city, Sishui County, Qufu City, Xintai City, Ningyang County and other regions. Hekou District, Lijin County and Zhanhua County are the “high-high” gathering areas over the years. It may be due to the large number of animal husbandry and animal product processing industries in these areas, there are many workers in these related industries, the backward traditional breeding methods, the poor behavior of family captivity, the mixed living of human and animal, and the weak awareness of protection, leading to more opportunities for infection. Over time, the number of counties (districts) covered by high-high clustered areas showed a downward trend. In recent years, Pingyi County of Linyi City, Ningyang County of Tai ‘an city and Qufu City of Jining became the new “high-high” clustered area. The range of “low-low” clustered had been shrinking year by year, and the new “high-high” clustered appears in the eastern region. Be alert to these new “high-high” clusters and strengthen targeted prevention and control measures.

The Spatio temporal cluster analysis identified one most likely cluster, one secondary cluster and one 2nd secondary cluster of human brucellosis incidence. The most likely cluster was located in the northern areas, including 37 counties and districts, which indicated that the prevention and control measures in these regions still need to be strengthened. The secondary clusters were distributed in the southern area, including Qufu City, Sishui County, mainly ranged from January 2019 to December 2021.This means that brucellosis moves from north to south; thus, prevention and control awareness should also be established in these areas. This Spatio-temporal cluster study of brucellosis helps to identify high-risk areas and time frames for brucellosis, to provide basis for government decision-making. According to the results of our study, we suggest to further strengthen surveillance and prevention and control measures in the southern region.

Direct contact with the infectious livestock remains the primary source of human brucellosis. the increasing demand for mutton and beef with the improvement of living standards, which has led to an increase in the transportation of cattle, sheep, and other animal products from the northern pastoral areas [48]. Inter-provincial movement of livestock products increase exposure to infected animals and contaminated livestock products. It is critical to strengthen the management and supervision of the agricultural market and the interprovincial transportation of livestock to prevent contaminated meat and dairy products from entering the market. It is suggested that the animal husbandry supervision department should strengthen the quarantine of cattle and sheep and the supervision of livestock circulation and transportation to prevent the flow of diseased animals and their livestock products.

Governments at all levels should establish joint prevention and control mechanism involving multiple departments among high-risk areas. Relevant departments in different regions should continue to strengthen both the prevention and control throughout the year, to contain the spread of brucellosis. Brucellosis is a reemerging disease for both humans and animals [49], health and animal husbandry departments should strengthen cooperation, strengthen information exchange and communication, regularly inform each other of the epidemic situation, and ensure that all prevention and control measures are in place. Increase investment in brucellosis prevention and control in key areas, strengthen training of medical personnel, effectively improve the diagnosis and treatment level of brucellosis in medical institutions, earnestly do a good job in the diagnosis and reporting of brucellosis cases, achieve early detection, diagnosis, reporting and treatment, standardize treatment, and reduce the chronic brucellosis.

Strengthen information dissemination and health education of brucellosis among key groups, such as family breeder, intervention packages have been issued to guide livestock breeding and livestock product production and processing occupational groups to do a good job in personal protection. Health education videos and leaflets has been used to increase the intensity and frequency of education and improve the protection awareness of employees. Comprehensive prevention and control measures to reduce brucellosis infection rate, to ensure the health of the people.

Conclusion

Brucellosis continues to be a common challenge in Shandong province, especially for counties in the northern region from March to August. In such higher-risk susceptible areas, more disease prevention and control measures should be taken to reduce the incidence of brucellosis and ensure the health of the people. Further researches should focus more on the influencing factors, such as the environmental and socio-economic factors, to provide more guidelines for policy makers to initiate control strategies.

Data Availability

The datasets generated during and analyzed during the current study are not publicly available but are available from the corresponding author on reasonable request.

Abbreviations

LLR:

Likelihood Rate

RR:

Relative Risk

LISA:

Local indicators of spatial association

References

  1. Sulayman S, Bora RS, Sabir J, et al. Brucellosis: current status of the disease and future perspectives. Postepy Mikrobiologii. 2020;59(4):337–44.

    Google Scholar 

  2. Pappas G, Akritidis N, Bosilkovski M, Tsianos E, Brucellosis. N Engl J Med. 2005;352(22):2325–36.

    CAS  PubMed  Google Scholar 

  3. Haque N, Bari MS, Hossain MA, et al. An overview of brucellosis. Mymensingh Med J. 2011;20(4):742–7.

    CAS  PubMed  Google Scholar 

  4. Zheng R, Xie S, Lu X, et al. A systematic review and Meta-analysis of epidemiology and clinical manifestations of human brucellosis in China. Biomed Res Int. 2018;2018:5712920.

    PubMed  PubMed Central  Google Scholar 

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

    PubMed  Google Scholar 

  6. Godfroid J, Al Dahouk S, Pappas G, Roth F, Matope G, Muma J, Marcotty T, Pfeiffer D, Skjerve E. A “One Health” surveillance and control of brucellosis in developing countries: moving away from improvisation. Comp Immunol Microbiol Infect Dis. 2013;36(3):241–8.

    PubMed  Google Scholar 

  7. 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  Google Scholar 

  8. 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  Google Scholar 

  9. Chen QL. Spatial-temporal analysis of the prevalence of brucellosis in my country and research on the health-related quality of life of patients [D]. Chinese Center for Disease Control and Prevention; 2017.

  10. Lai S, Zhou H, Xiong W, et al. Changing epidemiology of human brucellosis, China, 1955–2014. Emerg Infect Dis. 2017;23(2):184–94.

    PubMed  PubMed Central  Google Scholar 

  11. Flahaut B, Mouchart M, Martin ES, Thomas I. The local spatial autocorrelation and the kernel method for identifying black zones: a comparative approach. Accid Anal Prev. 2003;35(6):991–1004.

    PubMed  Google Scholar 

  12. Mattsson BJ, Zipkin EF, Gardner B, Blank PJ, Sauer JR, Royle JA. Explaining local-scale species distributions: relative contributions of spatial autocorrelation and landscape heterogeneity for an avian assemblage. PLoS ONE. 2013;8:e55097.

    CAS  PubMed  PubMed Central  Google Scholar 

  13. Viladomat J, Mazumder R, Mcinturff A, McCauley DJ, Hastie T. Assessing the significance of global and local correlations under spatial autocorrelation: a nonparametric approach. Biometrics. 2014;70(2):409–18.

    PubMed  PubMed Central  Google Scholar 

  14. Tillé Y, Dickson MM, Espa G, Giuliani D. Measuring the spatial balance of a sample: a new measure based on the Moran’s I index. Spat Stat. 2018;23:182–92.

    Google Scholar 

  15. Mao Y, Zhang N, Zhu B, Liu J, He R. A descriptive analysis of the Spatio-temporal distribution of intestinal infectious diseases in China. BMC Infect Dis. 2019;19(1):766.

    PubMed  PubMed Central  Google Scholar 

  16. Xiao-Ni H, Wei-Wei Z, Dan-Feng S, Hong L, Lian-Di Z, Bao-Guo L. Spatial pattern analysis of heavy metals in Beijing agricultural soils based on spatial autocorrelation statistics. Inter J Env Res Pub Heal. 2011;8:2074–89.

    Google Scholar 

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

    Google Scholar 

  18. Erkihun T, Amsalu T, Yihunie A. Spatial distribution and determinants of acute respiratory infection among under-five children in Ethiopia: ethiopian demographic health survey 2016. PLoS ONE. 2019;14(4):e0215572.

    Google Scholar 

  19. Zhang H, Yang L, Li L, Xu G, Zhang X. The epidemic characteristics and spatial autocorrelation analysis of hand, foot and mouth disease from 2010 to 2015 in Shantou, Guangdong, China. BMC Public Health. 2019;19(1):998.

    PubMed  PubMed Central  Google Scholar 

  20. T Waldhör. The spatial autocorrelation coefficient Moran’s I under heteroscedasticity. Stat Med. 2010;15(7–9):887–92.

    Google Scholar 

  21. Anselin L. Local Indicator of spatial AssociationLISA. GeogrAnal. 2010;27(2):93–115.

    Google Scholar 

  22. Kulldorff M. A spatial scan statistic. Commun Stat TheoryMethods. 1997;26:1481–96.

    Google Scholar 

  23. Kulldorff M, Heffernan R, Hartman J, et al. A space-time permutation scan statistic for Disease Outbreak Detection. PLoS Med. 2005;2(3):e59.

    PubMed  PubMed Central  Google Scholar 

  24. Jones SG, Martin K, Philip B. Influence of spatial resolution on space-time disease cluster detection. PLoS ONE. 2012;7:e48036.

    CAS  PubMed  PubMed Central  Google Scholar 

  25. Kulldorff M, Huang L, Pickle L, Duczmal L. An elliptic spatial scan statistic. Stat Med. 2006;25(22):3929–43.

    PubMed  Google Scholar 

  26. Kleinman KP, Abrams AM, Kulldorff M, Platt R. A model-adjusted space–time scan statistic with an application to syndromic surveillance. Epidemiol Infect. 2005;133(3):409–19.

    CAS  PubMed  PubMed Central  Google Scholar 

  27. Hunachew B, Wakgari D, Abera K, Delia G. Spatial, temporal, and spatio-temporal analysis of under-five diarrhea in southern Ethiopia. Trop Med Health. 2018;46:18.

    Google Scholar 

  28. Sluydts V, Heng S, Coosemans M, Van Roey K, Gryseels C, Canier L, et al. Spatial clustering and risk factors of malaria infections in Ratanakiri Province, Cambodia. Malar J. 2014;13:387.

    PubMed  PubMed Central  Google Scholar 

  29. Zhu B, Liu J, Fu Y, Zhang B, Mao Y. Spatio-temporal epidemiology of viral hepatitis in China (2003–2015): implications for prevention and control policies. Int J Environ Res Public Health. 2018;15(4):661.

    PubMed  PubMed Central  Google Scholar 

  30. Kulldorff M, Feuer EJ, Miller BA, Freedman LS. Breast cancer clusters in the northeast United States: a geographic analysis. Am J Epidemiol. 1997;146(2):161–70.

    CAS  PubMed  Google Scholar 

  31. Liu MY, Li QH, Zhang YJ, Ma Y, Liu Y, Feng W, et al. Spatial and temporal clustering analysis of tuberculosis in the mainland of China at the prefecture level, 2005–2015. Infect Dis Poverty. 2018;7:106.

    PubMed  PubMed Central  Google Scholar 

  32. Wang LY, Zhang WY, Ding F, Hu WB, Magalhaes RJ, Sun HL, et al. Spatio-temporal patterns of japanese encephalitis in China, 2002–2010. PLoS Negl Trop D. 2013;7:e2285.

    Google Scholar 

  33. Rubach MP, Halliday JE, Cleaveland S, Crump JA. Brucellosis in low-income and middle-income countries. Curr Opin Infect Dis. 2013;26(5):404–12.

    PubMed  PubMed Central  Google Scholar 

  34. Li F, He F, Sun J, Zhai Y, Jiang J, Lin J. Spatial and temporal analysis of severe fever with thrombocytopenia syndrome in Zhejiang Province, China, 2011–2015. J Infect Dev Ctries. 2019;13(1):35–43.

    PubMed  Google Scholar 

  35. Caprarelli G, Fletcher S. A brief review of spatial analysis concepts and tools used for mapping, containment and risk modelling of infectious diseases and other illnesses. Parasitology. 2014;141(5):581–601.

    PubMed  Google Scholar 

  36. Cui BY, Jiang H. Epidemiology of brucellosis in China, 2005–2016. Disease Surveillance. 2016;31(12):1018–22.

    Google Scholar 

  37. Shuman EK. Global climate change and infectious diseases. [J]. Int J Occup Environ Med. 2010;2(1):1061–3.

    Google Scholar 

  38. Liang PF, Zhao Y, Zhao JH, Pan DF, Guo ZQ. Human distribution and spatial-temporal clustering analysis of human brucellosis in China from 2012 to 2016. Infect Dis Poverty. 2020;9(1):142.

    PubMed  PubMed Central  Google Scholar 

  39. Maryam D, Youcef S, Adrian M, et al. Human brucellosis caused by raw dairy products: a review on the occurrence, major risk factors and prevention. IntJ Food Microbiol. 2018;292:39.

    Google Scholar 

  40. Glikman D. Human brucellosis in Israel-the Saga continues. Isr Med Assoc J. 2019;21:52–3.

    PubMed  Google Scholar 

  41. Peng C, Li YJ, Huang DS, Guan P. 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. 2020;25(1):1.

    CAS  PubMed  PubMed Central  Google Scholar 

  42. 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.

    PubMed  PubMed Central  Google Scholar 

  43. 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.

    PubMed  PubMed Central  Google Scholar 

  44. 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  Google Scholar 

  45. Pakzad R, Pakzad I, Safiri S, et al. Spatio-temporal analysis of brucellosis incidence in Iran from 2011 to 2014 using GIS. Int J Infect Dis. 2018;67:129–36.

    PubMed  Google Scholar 

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

    Google Scholar 

  47. 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;1:547.

    Google Scholar 

  48. 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.

    PubMed  PubMed Central  Google Scholar 

  49. Zhou K, Wu B, Pan H, et al. ONE Health Approach to address zoonotic brucellosis: a Spatio-temporal Associations Study between Animals and humans. Front Vet Sci. 2020;7:521.

    PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

The authors would like to thank all the reviewers for their review and revision suggestions, Thanks to staff from the city and county CDC for their efforts in brucellosis surveillance. Thanks to the corresponding author for his comments and suggestions, Thanks to all authors for their cooperation.

Funding

This work was supported by the National Science and Technology Major Project (No.2018ZX10714002), by the Shandong Traditional Chinese Medicine Science and Technology Project (No.2021Q001).

Author information

Authors and Affiliations

Authors

Contributions

KZQ participated in the study design and revised the manuscript. YXL carried out data analysis and drafted the manuscript. helped to revise the manuscript. FM, LY, YJM and CLX participated in the revision of the article. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Zengqiang Kou.

Ethics declarations

Competing interests

The authors declare no competing interests.

Ethics approval and consent to participate

We confirm that all methods were carried out in accordance with relevant guidelines and regulations of Helsinki declaration. Informed consent was obtained from all subjects and their legal guardians. All experimental protocols were approved by ethics Committee of Shandong CDC.

Consent for publication

Not applicable.

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 http://creativecommons.org/licenses/by/4.0/. 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 in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Yu, X., Fang, M., Li, Y. et al. Epidemiological characteristics and spatio-temporal analysis of brucellosis in Shandong province, 2015–2021. BMC Infect Dis 23, 669 (2023). https://doi.org/10.1186/s12879-023-08503-6

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12879-023-08503-6

Keywords