This article has Open Peer Review reports available.
Human brucellosis occurrences in inner mongolia, China: a spatio-temporal distribution and ecological niche modeling approach
© Jia and Joyner; licensee BioMed Central. 2015
Received: 19 September 2014
Accepted: 15 January 2015
Published: 3 February 2015
Brucellosis is a common zoonotic disease and remains a major burden in both human and domesticated animal populations worldwide. Few geographic studies of human Brucellosis have been conducted, especially in China. Inner Mongolia of China is considered an appropriate area for the study of human Brucellosis due to its provision of a suitable environment for animals most responsible for human Brucellosis outbreaks.
The aggregated numbers of human Brucellosis cases from 1951 to 2005 at the municipality level, and the yearly numbers and incidence rates of human Brucellosis cases from 2006 to 2010 at the county level were collected. Geographic Information Systems (GIS), remote sensing (RS) and ecological niche modeling (ENM) were integrated to study the distribution of human Brucellosis cases over 1951–2010.
Results indicate that areas of central and eastern Inner Mongolia provide a long-term suitable environment where human Brucellosis outbreaks have occurred and can be expected to persist. Other areas of northeast China and central Mongolia also contain similar environments.
This study is the first to combine advanced spatial statistical analysis with environmental modeling techniques when examining human Brucellosis outbreaks and will help to inform decision-making in the field of public health.
Brucellosis, a common zoonotic disease also referred to as Bang's disease, Crimean fever, Gibraltar fever, Malta fever, Maltese fever, Mediterranean fever, rock fever, and undulant fever, remains a major burden in both human and domesticated animal populations worldwide, especially in the Mediterranean region, Asia, the Middle East, Sub-Saharan Africa, Latin America and the Balkan Peninsula [1-3]. Causing more than 500,000 new human cases worldwide annually , Brucellosis is directly transmitted to humans via contact with animals that carry the pathogenic bacteria Brucella or infectious material during animal husbandry and meat processing, or indirectly through the consumption of unpasteurized dairy products . Infection among different animal species is mostly attributed to different infectious agents in the genus Brucella (B.). The primary causative agents include B. melitensis (sheep and goats), B. ovis (sheep and goat), B. abortus (cattle), B. suis (swine), and B. canis (dog) [4-6], although more have been detected in both domesticated and wildlife species [7,8]. Brucellosis itself yields low mortality rates, but it can cause substantial disabilities and weaknesses to the human immune system . The clinical presentation can be acute, subacute or chronic, varying from joint, muscle and back pain to flu-like symptoms, and even more serious conditions in different organ systems . However, it most commonly targets the reproductive system, resulting in up to a 40% increase in fetal wastage during the early stages of pregnancy and up to 2% of fetal deaths during the later stages of pregnancy for expectant women .
Details of production, multiplication and spread of the Brucella spp. among livestock have been described elsewhere , but a brief summary is given here. Initially, B. abortus replicates itself in regional lymph nodes of cows. Uterine infection occurs during the second trimester and is followed by placental inflammation that develops into placental disruption and endometritis ultimately leading to fetal death, after which the fetus is normally retained 1 to 3 days in utero and delivered with numerous bacteria expelled from the genital tract and spread out by various pathways. Brucellosis is detrimental to human health, and can lead to reductions in livestock production and subsequent economic losses, reduced milk/meat production, time lost by patients from normal daily activities, and increased medical costs [13,14]. Combination therapy with long-acting oxytetracycline and streptomycin has shown to moderate the effectiveness of the disease resulting in a reduction of shedding in most cows and complete elimination of infection in some cows. However, due to low treatment success rates, many countries (e.g., United States, European Union member countries, Australia and New Zealand) adopted eradication programs to slaughter the infected cows and quarantine the exposed herdmates until they could either be slaughtered or recertified as Brucellosis-free [12,15,16].
In China, Brucellosis remains a major public health issue. Previous studies indicated that human Brucellosis was endemic in 25 of 32 provinces (or autonomous regions) of mainland China . Inner Mongolia Autonomous Region has been most severely affected by Brucellosis in both humans and livestock since 1999, reporting the largest number of human Brucellosis cases across China [18,19], and accounting for approximately 40% of the total reported cases in China during 1999–2008 [18,20], 47.2% in 2010 , and almost 50% during 2005–2010 . However, few efforts have been specifically focused on the distribution of human Brucellosis cases in Inner Mongolia. One recent study revealed potential relationships between the incidence of human Brucellosis cases in China and some environmental factors, including temperature, rainfall, hours of sunshine, relative humidity and average wind velocity , but additional research is needed.
The goal of this study is two-fold: 1) examine the spatio-temporal distribution of reported human Brucellosis cases and incidence rates across Inner Mongolia during 1950–2010 and 2) examine the association of environmental factors with human Brucellosis occurrences in Inner Mongolia. Geographic Information Systems (GIS), remote sensing (RS), exploratory spatial data analysis (ESDA) and ecological niche modeling (ENM) were utilized to achieve the goals of this study.
No human subjects work was undertaken in this study. Human Brucellosis case data were extracted from annual reports and published literature. The annual reports provided summarized count data of patients reported by county and year. All data were anonymized.
Inner Mongolia consists of 12 administrative units at the municipal level (nine cities and three leagues), one level finer than province, and 89 subdivisions on the county level. The capital is Hohhot, and Baotou is the largest industrial city. In addition to Hohhot and Baotou, the other ten municipal units include Bayan Nur City, Chifeng City, Hulunbuir City, Ordos City, Tongliao City, Wuhai City, Wulanchabu City, Alxa League, Hinggan League, and Xilin Gol League. The earliest record related to animal Brucellosis in Inner Mongolia dates back to the year 1936 in Ulan Hot City of Hinggan League (Figure 1), called Wangye Temple at that time, where Brucella was found in an aborted bovine fetus.
The aggregated numbers of human Brucellosis cases in each decade from 1951 to 2000 and from 2001 to 2005 at the municipal level were extracted from published literature - the data during 1951–2003 and during 2004–2005 were originally from annual reports by the Inner Mongolia Center for Disease Control and Prevention (IMCDC) and Chinese Center for Disease Control and Prevention (CCDC), respectively . The annual numbers of human Brucellosis cases along with incidence rates during 2006–2010, calculated based on total county population and reported per 100,000, were collected from annual reports by the IMCDC. Each reported case must be confirmed by both clinical signs and serological tests, or isolation in accordance with the case definition of the World Health Organization (WHO) to maintain the same case definitions as cases reported in human Brucellosis studies in China .
The minimum, maximum and mean middle infrared radiation (MIR), land surface temperature (LST), and enhanced vegetation index (EVI) images were produced by Temporal Fourier Analysis [26,27], based on available 1 × 1 km data available between 2001 and 2008 from MODerate-resolution Imaging Spectroradiometer (MODIS) sensors on-board NASA’s Terra and Aqua satellites . EVI is an improved version of the conventional Normalized Difference Vegetation Index (NDVI). While NDVI has been used extensively in successful studies , EVI functions with higher fidelity through reducing saturation of the signal at high vegetation coverage and soil background effects . LST is a measure of surface, or “skin,” temperature . MIR, recorded by MODIS band 7 (2,105-2,155 nm), is often associated with the water content, surface temperature, and structure of vegetation canopies .
Predictor variables (data layers) in the study
Remotely sensed variables
EVI = 2.5 × (NIR-RED)/(NIR + 6.0 × RED–7.5 × BLUE + 1.0)
Maximum LST (°C)
Mean LST (°C)
Minimum LST (°C)
MODIS band 7
Maximum precipitation (mm)
Mean precipitation (mm)
Minimum precipitation (mm)
Annual mean temperature (°C)
Mean diurnal range(°C)
max temp - min temp
(bio2 / bio7) × 100
(std. dev. / mean) × 100
Max temperature of warmest month (°C)
Min temperature of coldest month (°C)
Temperature annual range (°C)
bio5 - bio6
Mean temperature of wettest quarter (°C)
Mean temperature of driest quarter (°C)
Mean temperature of warmest quarter (°C)
Mean temperature of coldest quarter (°C)
Annual precipitation (mm)
Precipitation of wettest month (mm)
Precipitation of driest month (mm)
(std. dev. / mean) × 100
Precipitation of wettest quarter (mm)
Precipitation of driest quarter (mm)
Precipitation of warmest quarter (mm)
Precipitation of coldest quarter (mm)
Travel time to the nearest major settlement of population size 50,000 or more (minutes)
Population density (person/km2) in 2000
Cattle density (heads/km2) in 2005
Small ruminants density (heads/km2) in 2005
Cattle density (heads/km2) in 2010
Goats density (heads /km2) in 2010
Sheep density (heads/km2) in 2010
Since cattle and small ruminants have been responsible for most human Brucellosis occurrences , the densities of cattle and small ruminants in 2005, with a resolution of 5×5 km, were obtained from the Food and Agriculture Organization (FAO) of the United Nations and resampled to 1×1 km. In 2010, the density of cattle was available at a resolution of 1×1 km, as well as the densities of sheep and goats used for representing the density of small ruminants which was not available in 2010 . The Global Rural–urban Mapping Project, Version 1 (GRUMPv1) was used to provide population numbers based on values from 1×1 km grids . A global map of accessibility in 2000, available as a gridded surface at a resolution of 1×1 km with pixel values representing minutes of travel time to the nearest major settlement of population size 50,000 or more, was also utilized for this study . The worldwide Land Use Systems (Version 1.1) was developed by the Land Tenure and Management Unit of the FAO and used to extract livestock-related areas.
The annual numbers of human Brucellosis cases at the municipal level from 1950 to 2010 were temporally aggregated by decade and visually compared. The annual incidence rates of human cases from 2006 to 2010 on the county level were plotted for a visual comparison over time. The χ2 (chi-square) two-sample test was implemented to determine statistically significant differences in incidence rates amongst counties between each subsequent and preceding year. A natural log transformation was performed prior to each t-test to reduce the influence of the skewness of the distribution of incidence rate values.
To reduce the potential impact of inaccurate observed incidence rates caused by an issue of “small case numbers” or incomplete reporting, a method of spatial smoothing was applied and resulted in the replacement of raw rates with adjusted estimates based on incidence rates in neighboring counties. This can be realized by assigning either the radius (distance) of a smoothing window or by defining neighbors and number of neighbors, and is predicated on the first law of geography which states that "everything is related to everything else, but near things are more related than distant things” . Using GeoDa, a spatial statistics and spatial data analysis software package, a first order Rook contiguity weights matrix was created for spatial smoothing, meaning that all polygons sharing a border of some length, more than a point-length border, with the target Thiessen polygon were defined as neighbors .
Local Moran’s I was employed to examine spatial clustering of the counties with high and low numbers of reported cases and incidence rates across Inner Mongolia, where a random permutation procedure (RPP) replicated the statistic 999 times to generate reference distributions . Through comparison of reference distributions, the statistical significance of spatial clusters was tested.
Ecological niche modeling
ENM, similar in many respects to species distribution modeling (SDM), is an innovative approach that associates the occurrence locations of infectious diseases (i.e., presence data) with a set of environmental variables resulting in a geographic range prediction for disease cases and agents. Maxent, used in this study, is a maximum entropy approach to presence-only distribution modeling that has shown a high predictive power even for small sample sizes [43,44]. However, the output from Maxent for evaluating the relative contribution of each independent variable to the model can be affected by high correlation amongst variables . Prior to modeling experiments, variable correlations were examined to eliminate issues of high multicollinearity. The Maxent modeling procedure assesses the importance of the individual and collective contributions that each variable makes to the distribution of human Brucellosis cases through 1) jackknife analysis, 2) average area under the curve (AUC) values for 10 model iterations and 3) average percentage contribution of each variable to the model. The prediction was projected outside of Inner Mongolia to the country of Mongolia to the north and eight adjacent Chinese provinces that border Inner Mongolia, which from east to west include Heilongjiang, Jilin, Liaoning, Hebei, Shanxi, Shannxi, Ningxia and Gansu provinces.
Obtaining point coordinates of Brucellosis presence that best represented conditions where human Brucellosis most likely occurs proved to be challenging since original data were aggregated over administrative units. A priori knowledge was needed to limit and reduce the suitable habitat area. For example, it is well-known that pastoral areas have the highest incidence rates of human Brucellosis [22,46]. Instead of simply using geometric centroids of administrative units as multiple studies have done [47-51], this issue was addressed through a series of dasymetric procedures based on the verity that human Brucellosis is most likely to occur on the grasslands where both people and small ruminants (primarily sheep and goats) inhabit. The density of small ruminants in most of Inner Mongolia ranged from 0 to 58 heads/km2, with abnormally high densities in just a few disjunct locations. Small ruminant (sheep and goat) densities were mapped in Inner Mongolia and the densities were classified into three categories based on a quantile classification method. Goat densities greater than ~22 heads/km2 were classified into the highest quantile, while sheep densities greater than ~16 heads/km2 were classified into the highest quantile. The 16–22 heads/km2 range matched that of a previous study  that indicated that the medium density threshold for small ruminants is ~18 heads/km2. Thus, a threshold of 18 heads/km2 was considered to be a reasonable cutoff value for the medium density of small ruminants in this study. In addition to small ruminant densities, moderate elevation (800–1600 meters) is considered to be more suitable for Brucellosis occurrences . Moreover, potential human Brucellosis occurrences were limited to the counties where total numbers of reported cases from 2006 to 2010 exceeded 50 persons.
Grasslands with low, moderate and high livestock densities were extracted from the land use raster layer and one centroid was created for each raster cell. The values of altitude (800–1600 meters), population density (>1 person/km2), small ruminant density in both 2005 and 2010 (>18 heads/km2) underlying each centroid were then used to further restrain the suspected occurrence locations. The case numbers from 2006 to 2010 were summed for each county and only the suspected locations falling within the counties with more than 50 cases were selected as “presence” points.
Multivariate regression models were built at the county level based on the remaining variables after removing high multicollinearity, to explore the impacts of various environmental and socioeconomic factors on the incidence rates of human Brucellosis over time. The incidence rates over counties in five individual years (2006–2010) and the average incidence rate over five years were separately used as dependent variables in six models, with the same set of independent variables used in the models and a stepwise method applied to each regression procedure.
Spatial-temporal distribution of reported human Brucellosis cases
Descriptive statistics of raw and natural log-transformed incidence rates of Brucellosis occurrences in humans by county during 2006-2010
Log Std. deviation
Log Std. error mean
Results (equal variances assumed) of t -test for the difference of incidence rates of human Brucellosis occurrences between each subsequent and preceding year during 2006-2010
t-test for equality of means
Std. error difference
2007 v.s. 2006
2008 v.s. 2007
2009 v.s. 2008
2010 v.s. 2009
2010 v.s. 2006
Ecological niche modeling
Two ecological niche models (ENM_1 and ENM_2) were produced to predict the potential risk areas for human Brucellosis. The ENM_1 used all suspected occurrence points available based on five years of reported human cases (2006–2010) and 17 environmental and socioeconomic variables after removing some variables that were found to have high multicollinearity with other variables. Of the total suspected points, 25% were randomly excluded from modeling and used for model testing and validation. The ENM_2 used the same 17 variables as ENM_1, and the suspected occurrence points produced based on three years of reported human cases (2006–2008), with the last two years of data (2009–2010) withheld from the model to be used as the testing (validation) dataset.
Analysis of variable contributions and model assessment
Counties (cities, districts or banners) with different extents of probability for human Brucellosis
Probability of presence
Counties (ordered by probability from high to low)
>0.5 – 0.57
West Ujimqin Banner, Plain and Bordered White Banner, Duolun County, Xulun Hoh Banner
>0.4 – 0.5
Bordered Yellow Banner, Xilinhot City, Sonid Right Banner, Abag Banner, Darhan Muminggan United Banner, Tumed Left Banner, Siziwang Banner, Hohhot City, Tumed Right Banner, East Ujimqin Banner, Hexigten Banner, Jungar Banner, Baiyun Mining District
>0.3 – 0.4
Urat Middle Banner, Horinger County, Linxi County, Chahar Right Back Banner, Sonid Left Banner, Urat Front Banner, Guyang County, Shangdu County, Dalad Banner, Baotou City
>0.2 – 0.3
Togtoh County, Qingshuihe County, Huade County, Chahar Right Middle Banner, Harqin Banner, Baarin Right Banner, Xinghe County, Donghe District, Wuchuan County, Hanggin Banner
>0.1 – 0.2
Liangcheng County, Zhuozi County, Taibus Banner, Otog Banner, Chifeng City, Chahar Right Front Banner, Ongniud Banner, Wuyuan County, Baarin Left Banner, Ejin Horo Banner, Dengkou County, Linhe District
>0.01 – 0.1
Hanggin Rear Banner, Uxin Banner, Horqin Right Front Banner, Jining District, Fengzhen City, Jarud Banner, Horqin Right Middle Banner, Ningcheng County, Aohan Banner, Erenhot City, Holingol City, Ar Khorchin Banner, Otog Front Banner, Tuquan County, New Barag Left Banner, Wuhai City, New Barag Right Banner, Urat Rear Banner, Ulan Hot City, Evenk Autonomous Banner
Variables of total numbers (persons or heads) derived from original variables of densities (persons/km 2 or heads/km 2 )
Total number of population in 2005
Total number of cattle in 2005
Total number of small ruminants in 2005
Total number of cattle in 2010
Total number of goats in 2010
Total number of sheep in 2010
Significant independent variables for the incidence rates of human Brucellosis cases in different years (p < 0.05)
Multiple studies have explored the spatial patterns of human Brucellosis in various countries including Azerbaijan , Italy , Germany , United States , Ecuador , and mainland China [14,19,55]. However, the most suitable study area for Brucellosis in China, Inner Mongolia, has not been well-studied. Moreover, limited efforts have primarily focused on the environmental dynamics of human Brucellosis occurrences . To fill the gaps, this study systematically examined the spatial-temporal distribution of reported human Brucellosis cases and incidence rates at the municipal level from 1950 to 2010 and at the county level from 2006 to 2010, as well as the impacts of environmental and socioeconomic factors on human Brucellosis occurrences at point-level and county-level locations.
The study represents the first effort to combine GIS, RS and ENM techniques for the study of human Brucellosis. With the rapid spread of advanced spatial techniques and fast growing volume of high-resolution spatio-temporal data, GIS and RS are playing an increasingly vital role in medical geography and spatial epidemiology over recent years. GIS and RS provide a powerful toolset and diversity of data sources for ENM; ENM in turn has advanced conventional spatio-temporal analysis by adding more mathematical and computational complexity.
The probability of presence for human Brucellosis based on the ENM represents the extent of the similarity of environmental and socioeconomic conditions at each location to those most suitable for human Brucellosis. The locations with high predicted probability of presence do not necessarily have a large number of reported cases in the past. They may just have suitable conditions and, consequently, the potential for future occurrences. Pilot surveillance programs should be launched to determine if any under reporting or reporting biases exists at those locations.
From the perspective of historical variations, although there has been a surprising increase from 1950 to 2010, increasing accessibility to health facilities for diagnosis and reporting may be partly responsible for the perceived increase in human Brucellosis cases, especially in Brucellosis epidemic regions where the incidence rates of human Brucellosis cases increased dramatically between 2001 and 2010 (e.g., Hulunbuir and Xinli Gol). The number of health facilities has also dramatically increased and improved in Inner Mongolia since 2003 when reporting networks were established. Surveillance is typically based on the diagnosis of human cases in the nearest health facilities, which are then communicated to the IMCDC by municipal health departments and subsequently forwarded to the CCDC through the National Notifiable Disease Surveillance System. This procedure may be affected by a lack of accurate diagnosis, imprecise symptoms, uneven case detection efforts, and reporting bias; therefore a spatially smoothed approach was used for reducing spurious outliers. Reporting bias may also, to some extent, result in a mismatch between highly-suspected occurrence points and areas with high risk for Brucellosis predicted by ENM models.
Inner Mongolia is well-known for its extensive prairies and grass-fed livestock including cattle and sheep, which are often the carriers for the bacteria Brucella. The region contained approximately 18.2% of the total sheep population during 2005–2010, larger than any other province in mainland China . Consequently, Inner Mongolia is considered a focal area for the study of human Brucellosis. A descriptive study in Inner Mongolia revealed that occupation (agriculture worker, shepherd, butcher, slaughter-house worker, and cattle dealer), risky practices (handling of ruminant abortions, skinning of stillborn lambs and kids, and crushing the umbilical cord of newborn lambs and kids with teeth), and certain dietary preferences (consuming unpasteurized and unboiled milk and fresh cheese) were correlated with the occurrence of Brucellosis in humans . Larger numbers of reported human cases during 1999–2008 occurred predominantly in a specific gender (male, 70.2%) and age range (30–59 years old, 64.7%) and is likely a result of these groups participating in risky occupations and practices. Occupational features of Brucellosis also help to explain why more cases are reported in the east (humid, grassland-dominant) than west (arid and semiarid, desert-dominant) and in rural areas over urban areas. Additional educational programs should be launched, specifically aimed at enhancing the awareness of work-place Brucellosis prevention for certain groups.
Brucellosis is most often carried by domesticated animals and evidence from a model of animal-human Brucellosis transmission shows that 90% of human Brucellosis cases are small-ruminant derived in Mongolia . In China, higher incidence rates of human Brucellosis were positively associated with the density of sheep and goats  and approximately 84.5% of the total human cases were found to be infected by B. melitensis . Especially in Inner Mongolia, this proportion reached 90.25% or higher during 1996–2010 . Therefore, small ruminant density in 2005 is used as an indicator for the possibility of human occurrences, and sheep and goats are combined to represent small ruminant densities in 2010 when small ruminant density is not available. Cattle was not associated with human Brucellosis occurrences in a previous study in China . In addition to livestock density, nomadicity might play an important role in the transmission of Brucellosis among animals in Inner Mongolia through the mixing of herds, which can affect the transmission in humans. This is a limitation of this study due to the complex movement of livestock.
The vaccination rate for sheep in Inner Mongolia was previously only 31.6% . A study concluded that vaccinating adult sheep alone is insufficient in eradicating Brucellosis based on the prediction that Brucellosis would persist for a long period of time even though all sheep were supposedly vaccinated twice per year . Since 2011, the government in Inner Mongolia began to require vaccinations for all sheep twice per year. Given that infected sheep are still common across Inner Mongolia, especially amongst local herdsman, simultaneous disinfection and vaccination was suggested, as well as regular sheep surveillance which has been realized for some areas in China .
Environmental modeling for human Brucellosis is uncommon because Brucellosis is not as climate-sensitive as some infectious diseases, such as malaria [57,58] and dengue fever [45,59] among other diseases. However, varying environmental conditions can result in a lack of water or grass in the pasturing areas, which increases disease susceptibility in animals with low resistance. Altitude especially can play an important role in the types of vegetation that can grow in certain areas, precipitation levels, temperature ranges, and other environmental conditions important to the survival and transmission of Brucella. Consequently, altitude was the most influential abiotic variable within the ENM. Additionally, Brucellosis is capable of being transmitted by fomites. The Brucella organism can survive for several months in water, aborted fetuses, manure, hay, contaminated equipment and clothes, and especially in conditions with high humidity, low temperature, and no sunlight . Although the theory is on a micro-level, similar environmental impacts have been revealed by a macro-level study, where lower temperature and less sunshine in winter and spring, and a 1–2 month incubation period have been associated with epidemic peaks from March to August, especially in the month of May . Inhabitants of Inner Mongolia also spend more time working outside during epidemic months than during the harsh winter months, increasing their exposure to Brucellosis infection. In addition to environmental and biological determinants of the organism, different cultural and livestock management practices among communities in the highlands might impact differences in risk factors, but these differences warrant additional epidemiological field work before they can be confirmed.
After Inner Mongolia, the five provinces with the most human Brucellosis cases reported during 2005–2010 (Shannxi, Heilongjiang, Hebei, Jilin and Shanxi) all share a border with Inner Mongolia, which was also detected as the primary cluster of incidence rates  and accounted for most major occurrences of human Brucellosis cases in China . Therefore, the ENM prediction was extended outside of Inner Mongolia to all bordering counties. Some human cases were previously reported in non-epidemic parts of China, such as Guangdong and Fujian provinces in southeastern China, which might be attributed to increased transportation of un-quarantined and unvaccinated animals and travelling of people from epidemic areas [11,18,22]. It is recommended to examine travel history as part of a more inclusive epidemiological investigation in non-endemic regions.
A major limitation of the study pertains to the quality of disease data that were available (as is the case with most ENM and SDM research [47-51]) because the true occurrence locations were not available. Brucellosis is not a contagious disease among humans, so predicting the high risk areas of human Brucellosis requires knowledge of the locations where people are most likely infected. The dasymetric procedure provided a much more certain location-to-environment relationship for modeling purposes and was superior to the alternative use of average environmental conditions at the county or municipal levels. According to the study of the occupation of patients during the 1991–2005 period in Inner Mongolia, peasants and herders consisted of 79.36% of all cases . This provided further evidence in selecting appropriate locations for human Brucellosis. Other limitations include the inability to consider the differences in the accessibility to health facilities and diagnostic services between the counties reporting larger and smaller numbers of human cases.
This study provides a new direction for human Brucellosis research and greatly contributes to our knowledge of the roles that various environmental and socioeconomic factors play in the distribution and spread of human Brucellosis occurrences. Moreover, the output of this study will promote comparisons with future research at an international scale and provide a new perspective to inform decision-making in the field of public health in Inner Mongolia. This study also serves as a predecessor to broader and more in-depth studies of human Brucellosis in Inner Mongolia. Expansions in both spatial and temporal scales are warranted and could provide further insight.
The main conclusions of this study are as follows: 1) Human Brucellosis occurrences exhibit a spatial trend gradually changing in the number of cases from east and central Inner Mongolia to the west, with the incidence rates in central Inner Mongolia higher than those of other regions; 2) The major variables contributing to the model include the density of small ruminants (especially sheep), altitude, precipitation seasonality, and temperature annual range; 3) The incidence rates in the highest regions may be skewed by smaller populations; 4) The density of sheep (positively) and the mean of LST (negatively) and MIR (positively) were correlated with the incidence rates of human Brucellosis at the county level; 5) It is expected that parts of Mongolia and northern Hebei, Shanxi and Shannxi, in addition to Inner Mongolia, have some risk for human Brucellosis occurrence, based on their suitable conditions; 6) Counties in Inner Mongolia are classified into a series of categories with differing levels of probability for the presence of human Brucellosis.
This research was partially funded by the Robert C. West Graduate Student Field Research Award from the Department of Geography and Anthropology at Louisiana State University. We thank the East Tennessee State University Office of Research and Sponsored Programs for partial funding support for this project. We would also like to thank the reviewers for their many great comments and suggestions. They have tremendously improved this manuscript.
- Corbel MJ. Brucellosis: an overview. Emerg Infect Dis. 1997;3(2):213–21.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Gwida M, Al Dahouk S, Melzer F, Rosler U, Neubauer H, Tomaso H. Brucellosis - regionally emerging zoonotic disease? Croat Med J. 2010;51(4):289–95.View ArticlePubMedPubMed CentralGoogle Scholar
- Abdullayev R, Kracalik I, Ismayilova R, Ustun N, Talibzade A, Blackburn JK. Analyzing the spatial and temporal distribution of human brucellosis in Azerbaijan (1995–2009) using spatial and spatio-temporal statistics. BMC Infect Dis. 2012;12:185.View ArticlePubMedPubMed CentralGoogle Scholar
- Crawford RP, Huber JD, Adams BS. Epidemiology and surveillance. In: Animal Brucellosis. Boca Raton: CRC Press; 1990. p. 131–51.Google Scholar
- Zinsstag J, Roth F, Orkhon D, Chimed-Ochir G, Nansalmaa M, Kolar J, et al. A model of animal-human brucellosis transmission in Mongolia. Prev Vet Med. 2005;69(1–2):77–95.View ArticlePubMedGoogle Scholar
- Alexander KA, Blackburn JK, Vandewalle ME, Pesapane R, Baipoledi EK, Elzer PH. Buffalo, bush meat, and the zoonotic threat of brucellosis in Botswana. PLoS One. 2012;7(3):e32842.View ArticlePubMedPubMed CentralGoogle Scholar
- Godfroid J, Scholz HC, Barbier T, Nicolas C, Wattiau P, Fretin D, et al. Brucellosis at the animal/ecosystem/human interface at the beginning of the 21st century. Prev Vet Med. 2011;102(2):118–31.View ArticlePubMedGoogle Scholar
- Plumb GE, Olsen SC, Buttke D. Brucellosis: 'One Health' challenges and opportunities. Rev Sci Tech. 2013;32(1):271–8.PubMedGoogle Scholar
- Dean AS, Crump L, Greter H, Hattendorf J, Schelling E, Zinsstag J. Clinical manifestations of human brucellosis: a systematic review and meta-analysis. PLoS Negl Trop Dis. 2012;6(12):e1929.View ArticlePubMedPubMed CentralGoogle Scholar
- Memish ZA, Balkhy HH. Brucellosis and international travel. J Travel Med. 2004;11(1):49–55.View ArticlePubMedGoogle Scholar
- Christensen BW MD, TM H. Female reproductive disorders. In: SB P, editor. Large Animal Internal Medicine. 4th ed. St. Louis, MO: Mosby-Elsevier; 2009. p. 1461.Google Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- 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.View ArticlePubMedPubMed CentralGoogle Scholar
- Godfroid J, Kasbohrer A. Brucellosis in the European Union and Norway at the turn of the twenty-first century. Vet Microbiol. 2002;90(1–4):135–45.View ArticlePubMedGoogle Scholar
- Tweddle NE, Livingstone P. Bovine tuberculosis control and eradication programs in Australia and New Zealand. Vet Microbiol. 1994;40(1–2):23–39.View ArticlePubMedGoogle Scholar
- Shang D, Xiao D, Yin J. Epidemiology and control of brucellosis in China. Vet Microbiol. 2002;90(1–4):165–82.Google Scholar
- Zhang WY, Guo WD, Sun SH, Jiang JF, Sun HL, Li SL, et al. Human brucellosis, Inner Mongolia, China. Emerg Infect Dis. 2010;16(12):2001–3.View ArticlePubMedPubMed CentralGoogle Scholar
- 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 IJID Off Publ Int Soc Infect Dis. 2013;17(5):e289–92.Google Scholar
- Chen Y, Ke Y, Wang Y, Yuan X, Zhou X, Jiang H, et al. Changes of predominant species/biovars and sequence types of Brucella isolates, Inner Mongolia. Chin BMC Infect Dis. 2013;13:514.View ArticleGoogle Scholar
- Hou Q, Sun X, Zhang J, Liu Y, Wang Y, Jin Z. Modeling the transmission dynamics of sheep brucellosis in Inner Mongolia Autonomous Region, China. Math Biosci. 2013;242(1):51–8.View ArticlePubMedGoogle Scholar
- 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:547.View ArticlePubMedPubMed CentralGoogle Scholar
- Angerer J, Han G, Fujisaki I, Havstad K. Climate change and ecosystems of Asia with emphasis on Inner Mongolia and Mongolia. Rangelands. 2008;30(3):46–51.View ArticleGoogle Scholar
- Xue Y. The impact of desertification in the Mongolian and the Inner Mongolian grassland on the regional climate. J Clim. 1996;9(9):2173–89.View ArticleGoogle Scholar
- Chen WBCQ ZHJ, Chi R, Ma Y. Analysis of epidemic characteristics on brucellosis in Inner Mongolia. Chinese J Control of Endemic Dis. 2008;23(1):56–8.Google Scholar
- Rogers DJ. Satellites, space, time and the African trypanosomiases. Adv Parasitol. 2000;47:129–71.View ArticlePubMedGoogle Scholar
- Rogers DJ, Randolph SE, Snow RW, Hay SI. Satellite imagery in the study and forecast of malaria. Nature. 2002;415(6872):710–5.View ArticlePubMedPubMed CentralGoogle Scholar
- Scharlemann JP, Benz D, Hay SI, Purse BV, Tatem AJ, Wint GW, et al. Global data for ecology and epidemiology: a novel algorithm for temporal Fourier processing MODIS data. PLoS One. 2008;3(1):e1408.View ArticlePubMedPubMed CentralGoogle Scholar
- Kitron U, Kazmierczak JJ. Spatial analysis of the distribution of Lyme disease in Wisconsin. Am J Epidemiol. 1997;145(6):558–66.View ArticlePubMedGoogle Scholar
- Huete A, Didan K, Miura T, Rodriguez EP, Gao X, Ferreira LG. Overview of the radiometric and biophysical performance of the MODIS vegetation indices. Remote Sens Environ. 2002;83(1):195–213.View ArticleGoogle Scholar
- Sims DA, Rahman AF, Cordova VD, El-Masri BZ, Baldocchi DD, Bolstad PV, et al. A new model of gross primary productivity for North American ecosystems based solely on the enhanced vegetation index and land surface temperature from MODIS. Remote Sens Environ. 2008;112(4):1633–46.View ArticleGoogle Scholar
- Tatem AJ, Snow RW, Hay SI. Mapping the environmental coverage of the INDEPTH demographic surveillance system network in rural Africa. Tropical Med Int Health TM IH. 2006;11(8):1318–26.View ArticleGoogle Scholar
- Baptista Rosas RC, Hinojosa A, Riquelme M. Ecological niche modeling of Coccidioides spp. in western North American deserts. Ann N Y Acad Sci. 2007;1111:35–46.View ArticlePubMedGoogle Scholar
- Gonzalez C, Wang O, Strutz SE, Gonzalez-Salazar C, Sanchez-Cordero V, Sarkar S. Climate change and risk of leishmaniasis in north america: predictions from ecological niche models of vector and reservoir species. PLoS Negl Trop Dis. 2010;4(1):e585.View ArticlePubMedPubMed CentralGoogle Scholar
- Papeş M, Gaubert P. Modelling ecological niches from low numbers of occurrences: assessment of the conservation status of poorly known viverrids (Mammalia, Carnivora) across two continents. Divers Distrib. 2007;13(6):890–902.View ArticleGoogle Scholar
- Ron L, Benitez W, Speybroeck N, Ron J, Saegerman C, Berkvens D, et al. Spatio-temporal clusters of incident human brucellosis cases in Ecuador. Spatial Spatio-Temporal Epidemiol. 2013;5:1–10.View ArticleGoogle Scholar
- Robinson TP, Wint GW, Conchedda G, Van Boeckel TP, Ercoli V, Palamara E. Mapping the global distribution of livestock. PLoS One. 2014;9(5):e96084.View ArticlePubMedPubMed CentralGoogle Scholar
- Center for International Earth Science Information Network - CIESIN - Columbia University, International Food Policy Research Institute - IFPRI. Global Rural–urban Mapping Project, Version 1 (GRUMPv1): Population Density Grid. In: The World Bank, Centro Internacional de Agricultura Tropical - CIAT. Palisades, NY: NASA Socioeconomic Data and Applications Center (SEDAC); 2011.Google Scholar
- Nelson A. Estimated Travel Time to the Nearest City of 50,000 or More People in Year 2000. Ispra, Italy: In. Edited by Commission GEMU-JRCotE; 2008.Google Scholar
- Tobler WR. A computer movie simulating urban growth in the Detroit region. Econ Geogr. 1970;46:234–40.View ArticleGoogle Scholar
- Anselin L, Syabri I, Kho Y. GeoDa: an introduction to spatial data analysis. Geogr Anal. 2006;38(1):5–22.View ArticleGoogle Scholar
- Anselin L. Local indicators of spatial association—LISA. Geogr Anal. 1995;27(2):93–115.View ArticleGoogle Scholar
- Phillips, S. J., Dudík, M., & Schapire, R. E. A maximum entropy approach to species distribution modeling. In Proceedings of the twenty-first international conference on Machine learning (p. 83). ACM. (2004, July)Google Scholar
- Phillips SJ, Anderson RP, Schapire RE. Maximum entropy modeling of species geographic distributions. Ecol Model. 2006;190(3):231–59.View ArticleGoogle Scholar
- Machado-Machado EA. Empirical mapping of suitability to dengue fever in Mexico using species distribution modeling. Appl Geogr. 2012;33:82–93.View ArticleGoogle Scholar
- McDermott JJ, Arimi SM. Brucellosis in sub-Saharan Africa: epidemiology, control and impact. Vet Microbiol. 2002;90(1–4):111–34.View ArticlePubMedGoogle Scholar
- Mischler P, Kearney M, McCarroll JC, Scholte RG, Vounatsou P, Malone JB. Environmental and socio-economic risk modelling for Chagas disease in Bolivia. Geospat Health. 2012;6(3):S59–66.View ArticlePubMedGoogle Scholar
- Fischer D, Thomas SM, Niemitz F, Reineking B, Beierkuhnlein C. Projection of climatic suitability for Aedes albopictus Skuse (Culicidae) in Europe under climate change conditions. Global Planet Change. 2011;78(1):54–64.View ArticleGoogle Scholar
- Benedict MQ, Levine RS, Hawley WA, Lounibos LP. Spread of the tiger: global risk of invasion by the mosquito Aedes albopictus. Vector Borne Zoonotic Dis. 2007;7(1):76–85.View ArticlePubMedPubMed CentralGoogle Scholar
- Escobar LE, Peterson AT, Favi M, Yung V, Pons DJ, Medina-Vogel G. Ecology and geography of transmission of Two Bat-borne rabies lineages in Chile. PLoS Negl Trop Dis. 2013;7(12):e2577.View ArticlePubMedPubMed CentralGoogle Scholar
- Medley KA. Niche shifts during the global invasion of the Asian tiger mosquito, Aedes albopictus Skuse (Culicidae), revealed by reciprocal distribution models. Glob Ecol Biogeogr. 2010;19(1):122–33.View ArticleGoogle Scholar
- Mancini FR, Bella A, Graziani C, Marianelli C, Mughini-Gras L, Pasquali P, et al. Trends of human brucellosis in Italy, 1998–2010. Epidemiol Infect. 2014;142(6):1188–95.View ArticlePubMedGoogle Scholar
- Al Dahouk S, Neubauer H, Hensel A, Schoneberg I, Nockler K, Alpers K, et al. Changing epidemiology of human brucellosis, Germany, 1962–2005. Emerg Infect Dis. 2007;13(12):1895–900.View ArticlePubMedPubMed CentralGoogle Scholar
- Chomel BB, DeBess EE, Mangiamele DM, Reilly KF, Farver TB, Sun RK, et al. Changing trends in the epidemiology of human brucellosis in California from 1973 to 1992: a shift toward foodborne transmission. J Infect Dis. 1994;170(5):1216–23.View ArticlePubMedGoogle Scholar
- Chen Z, Zhang W, Ke Y, Wang Y, Tian B, Wang D, et al. High-risk regions of human brucellosis in china: implications for prevention and early diagnosis of travel-related infections. Clin infect Dis An Off Publ Infect Dis Soc Am. 2013;57(2):330–2.View ArticleGoogle Scholar
- CAHY. China animal husbandry yearbook. Beijing: China Agriculture Press; 2005.Google Scholar
- Ayala D, Costantini C, Ose K, Kamdem GC, Antonio-Nkondjio C, Agbor JP, et al. Habitat suitability and ecological niche profile of major malaria vectors in Cameroon. Malar J. 2009;8:307.View ArticlePubMedPubMed CentralGoogle Scholar
- Moffett A, Shackelford N, Sarkar S. Malaria in Africa: vector species' niche models and relative risk maps. PLoS One. 2007;2(9):e824.View ArticlePubMedPubMed CentralGoogle Scholar
- Peterson AT, Martinez-Campos C, Nakazawa Y, Martinez-Meyer E. Time-specific ecological niche modeling predicts spatial dynamics of vector insects and human dengue cases. Trans R Soc Trop Med Hyg. 2005;99(9):647–55.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. 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.