Open Access

Spatial transmission and meteorological determinants of tuberculosis incidence in Qinghai Province, China: a spatial clustering panel analysis

  • Hua-Xiang Rao1,
  • Xi Zhang2,
  • Lei Zhao1,
  • Juan Yu3,
  • Wen Ren1,
  • Xue-Lei Zhang1,
  • Yong-Cheng Ma4,
  • Yan Shi4,
  • Bin-Zhong Ma4,
  • Xiang Wang1,
  • Zhen Wei1,
  • Hua-Fang Wang1 and
  • Li-Xia Qiu1Email author
Infectious Diseases of Poverty20165:45

DOI: 10.1186/s40249-016-0139-4

Received: 25 September 2015

Accepted: 26 April 2016

Published: 2 June 2016

Abstract

Background

Tuberculosis (TB) is the notifiable infectious disease with the second highest incidence in the Qinghai province, a province with poor primary health care infrastructure. Understanding the spatial distribution of TB and related environmental factors is necessary for developing effective strategies to control and further eliminate TB.

Methods

Our TB incidence data and meteorological data were extracted from the China Information System of Disease Control and Prevention and statistical yearbooks, respectively. We calculated the global and local Moran’s I by using spatial autocorrelation analysis to detect the spatial clustering of TB incidence each year. A spatial panel data model was applied to examine the associations of meteorological factors with TB incidence after adjustment of spatial individual effects and spatial autocorrelation.

Results

The Local Moran’s I method detected 11 counties with a significantly high-high spatial clustering (average annual incidence: 294/100 000) and 17 counties with a significantly low-low spatial clustering (average annual incidence: 68/100 000) of TB annual incidence within the examined five-year period; the global Moran’s I values ranged from 0.40 to 0.58 (all P-values < 0.05). The TB incidence was positively associated with the temperature, precipitation, and wind speed (all P-values < 0.05), which were confirmed by the spatial panel data model. Each 10 °C, 2 cm, and 1 m/s increase in temperature, precipitation, and wind speed associated with 9 % and 3 % decrements and a 7 % increment in the TB incidence, respectively.

Conclusions

High TB incidence areas were mainly concentrated in south-western Qinghai, while low TB incidence areas clustered in eastern and north-western Qinghai. Areas with low temperature and precipitation and with strong wind speeds tended to have higher TB incidences.

Keywords

Tuberculosis incidence Meteorological factors Spatial clustering Spatial panel data model

Multilingual abstracts

Please see Additional file 1 for translations of the abstract into the six official working languages of the United Nations.

Background

Tuberculosis (TB) remains a major public health burden in many developing countries [13]. According to the World Health Organization TB annual report, in 2013, the number of new reported TB cases in the world was estimated at 11.4 million. Among the 22 high TB burden countries and regions, mainland China ranks second, with 1.3 million cases, giving an incidence of 98/100 000 [4]. Based on the data from the fifth TB epidemiological sampling survey in China, the TB number was higher in the western regions compared with the central and eastern regions [5]. The occurrence of TB in China has obvious periodic and seasonal features, more frequently occurring in the winter and spring, which suggests that TB might be associated with meteorological factors [6]. Current evidence has suggested that, besides the traditional factors such as genetic susceptibility [7], sex [8], education, ethnicity, drinking [9], smoking [10], and related diseases, several ecological factors, including geographic, climatic, and socioeconomic factors, also have critical impacts on the prevalence of TB [1113].

Understanding the spatial variations in TB prevalence and its determinants is crucial for improved targeting of interventions and resources. Many geospatial analytical methods, such as spatial autocorrelation analysis (Moran’s I and Getis-Ord G) [1416] and space-time scan statistic (SaTScan) methods [1719] have been used for understanding TB and other public health problems [20, 21]. In China, county-level studies have used various spatial epidemiological methods to identify clustering of health conditions, including notifiable pandemic influenza A in Hong Kong [22], and TB in Linyi [2], Beijing [5], and Xinjiang [23]. However, there are currently no related studies on geospatial distribution of TB in Qinghai province.

Moreover, while several studies have shown that the TB incidence might be related with the temperature, precipitation, and wind speed [24], few studies have considered the modifier effects from time and spatial factors in the relationship between meteorological factors and TB incidence. With further research of spatial econometrics and the rapid development of computer technology, spatial autocorrelation and spatial panel data models are becoming useful tools in the analysis of spatiotemporal data, and are gradually being applied to the research of infectious diseases [2426]. Therefore, in this study, we aimed to understand the distribution of TB and to explore the associations between meteorological factors and TB incidence using spatial autocorrelation analysis and a spatial panel data model based on surveillance data from the Qinghai Center for Disease Control and Prevention.

Methods

Qinghai province, which comprises 8 cities, including a total of 46 counties, is located between longitude 89°35′ and 103°04′ East, and latitude 31°40′ and 39°19′ North (Fig. 1). As an underdeveloped region in north-western China, it has a higher annual incidence of TB than other regions of the country. The average altitude is 3 000-5 000 m, with a typical plateau continental climate, which includes little rain, low temperatures, and long sunshine hours.
Fig. 1

Location of the study areas, Qinghai Province, China. The map was created using the ArcGIS software (version 10.0, ESRI Inc., Redlands, CA, USA)

Tuberculosis incidence data

In this study, we focused on the cases of pulmonary TB. Our TB data were based on the China Information System for Disease Control and Prevention (CISDCP, http://1.202.129.170/UVSSERVER2.0), which was established in 2005. TB cases were diagnosed using X-ray, pathogen detection, and pathologic diagnosis according to the diagnosis criteria recommended by the National Health and Family Planning Commission of the People’s Republic of China (Former Ministry of Health) in 2008 [6]. The relevant information, including age, sex, occupation, diagnostic category, and diagnostic date, was collected to analyse the epidemic characteristics of TB in Qinghai province.

TB is a notifiable infectious disease in China; all cases must be reported online within 24 h after diagnosis in the hospital. We collected the county- and city-level data from January 2009 to December 2013 in Qinghai and randomly selected 261 of all 771 medical institutions in Qinghai province and checked all medical records of these selected institutions to confirm whether there were missing TB cases or not. During the period from 2009 to 2013, in all 46 counties, a total of 27 665 TB cases and 51 TB-related deaths were identified; no missing cases or outbreaks were declared.

Environmental data

Our meteorological data were based on the statistical data from the Qinghai statistics office, which is reported yearly by the meteorological bureau of each city. In this study, we focused on four main meteorological factors, including the monthly average temperature (MAT, °C), precipitation (MP, mm), total sunshine hours (MSH, hours), and wind speed (MAWS, m/s). Considering a time lag between infectious disease development and meteorological factors, the meteorological data were collected from July 2008 to December 2013.

Statistical methods

Spatial autocorrelation analysis method

Spatial autocorrelation analysis was conducted by using Open GeoDa software (GeoDa Center for Geospatial Analysis and Computation, Arizona State University, AZ, USA) and used to identify the spatial clustering of the annual TB incidence of all 46 counties [2, 5, 27, 28]. The row standardized first-order contiguity Rook neighbours were used as the criterion for identifying neighbours in this paper. In this rule, if regions i and j are neighbours and share a boundary, w ij  = 1; otherwise, w ij  = 0. Global Moran’s I was calculated to test the spatial autocorrelation of all counties in Qinghai province, ranging from -1 to +1 [5]. Positive/negative spatial autocorrelation occurs when Moran’s I is close to +1/-1, which indicates that areas with similar (high-high or low-low)/dissimilar (high-low or low-high) incidence of TB are clustered together [29]. Monte Carlo randomization (9999 permutations) was employed to assess the significance of Moran’s I, with the null hypothesis being that the distribution of TB incidence in Qinghai province is completely spatially random [27]; in other words, that the counties with high and low TB incidence are randomly distributed across the study area [30]. If the test is significant (P ≤ 0.05), this suggests a clustering/dispersing of the TB incidence [2, 16]. Subsequently, we used local indicators of spatial association (LISA; Local Moran’ I) analysis and a Moran scatter plot to examine the spatial autocorrelation of each county in Qinghai province and to determine the locations of the clusters [31]. Moran’s plot shows high-high and low-low clustering in the upper right and lower left quadrants, respectively. Statistically significant high-high, low-low, and outlier local clusters (high-low and low-high) were visualized using a cluster map with county boundaries [15].

Spatial panel analysis method

In our study, we used data of TB incidence and meteorological factors, which were collected from different cities monthly, to explore the associations of meteorological factors with TB incidence. The data were collected at different times (monthly) and in different areas (cities), and were hence considered repeated observations data, referred to as panel data (also known as pooled time series and cross-section data). The spatial heterogeneity and spatial dependence in different cities need to be considered in the analysis [5, 32]. A spatial panel data model is able to address data with spatial dependence and also enables researchers to consider spatial heterogeneity, which typically refers to data containing continuous observations of a number of spatial units [2426]. Compared to traditional methods based on cross-sectional or time series models, spatial panel data models are more informative and contain more variation and less collinearity among the variables [26, 33]. The use of panel data results in a greater availability of degrees of freedom, and hence increases the efficiency of the estimation [26, 34]. Therefore, we used a spatial panel data model to examine the association between the TB incidence and meteorological factors after adjustment for the spatial confounders. As the distribution of the TB incidence rate is highly skewed, log transformation of the TB incidence was used in the analyses, according to the following formula: log incidence = lg (TB incidence). The unit root test and co-integration test were conducted to confirm the stationary of the data [35]. The spatial panel data model analyses were conducted by using Matlab R2009a (Mathworks Inc., Natick, MA, USA), and the significant level was 0.05.

Results

Annually, the incidences of TB in Qinghai province were 93, 87, 93, 112, and 106 per 100 000 people from 2009 to 2013, accounting for 12.24, 13.96, 17.13, 19.00, and 16.84 % of all reported infectious diseases, respectively (Table 1). The high incidence areas were mainly concentrated in the cities of Yushu and Guoluo, while the top three counties of TB annual incidence were Maduo (656/100 000), Jiuzhi (461/100 000), and Zaduo (393/100 000) in the southwest of Qinghai. The low incidence areas were mainly concentrated in the cities of Haixi and Xining, with the bottom three counties being Mangya (15/100 000), Delingha (39/100 000), and Dachaidan (42/100 000) in the northwest of Qinghai (Fig. 2). The TB incidence rate showed significant periodicity and seasonality, reaching a seasonal peak around April and then decreasing to a trough in December (Fig. 3). The significant meteorological characteristics in Qinghai province included the strong sunlight and relative low temperature throughout the year (Table 2).
Table 1

Characteristics of tuberculosis cases in Qinghai Province, China, 2009-2013

Variables

2009

2010

2011

2012

2013

Report cases (incidence rate, 1/100 000)

Death cases (mortality rate, 1/100 000)

Report cases (incidence rate, 1/100 000)

Death cases (mortality rate, 1/100 000)

Report cases (incidence rate, 1/100 000)

Death cases (mortality rate, 1/100 000)

Report cases (incidence rate, 1/100 000)

Death cases (mortality rate, 1/100 000)

Report cases (incidence rate, 1/100 000)

Death cases (mortality rate, 1/100 000)

Infectious disease

42 006 (757.82)

50 (0.90)

34 865 (625.60)

29 (0.52)

30 541 (542.79)

22 (0.39)

33 514 (589.86)

42 (0.74)

35 954 (627.28)

38 (0.66)

Pulmonary tuberculosis

5 141 (92.75)

21 (0.38)

4 868 (87.35)

6 (0.11)

5 232 (92.99)

10 (0.18)

6 369 (112.10)

7 (0.12)

6 055 (105.64)

7 (0.12)

Age (years)

0-

12 (15.57)

0 (0.00)

11 (14.15)

1 (1.29)

16 (21.07)

0 (0.00)

14 (19.57)

0 (0.00)

12 (16.99)

0 (0.00)

1-

1 (1.32)

0 (0.00)

5 (6.58)

0 (0.00)

5 (6.73)

0 (0.00)

14 (19.07)

0 (0.00)

5 (6.94)

0 (0.00)

2-

2 (2.62)

0 (0.00)

0 (0.00)

0 (0.00)

2 (2.70)

0 (0.00)

9 (12.37)

0 (0.00)

7 (9.87)

0 (0.00)

3-

1 (1.30)

0 (0.00)

1 (1.32)

0 (0.00)

5 (6.72)

0 (0.00)

10 (13.76)

0 (0.00)

7 (9.93)

0 (0.00)

4-

3 (3.83)

0 (0.00)

3 (3.90)

0 (0.00)

3 (3.99)

0 (0.00)

7 (9.53)

0 (0.00)

7 (9.85)

0 (0.00)

5-

3 (3.76)

0 (0.00)

7 (8.95)

0 (0.00)

7 (9.14)

0 (0.00)

7 (9.48)

0 (0.00)

3 (4.20)

0 (0.00)

6-

8 (9.89)

0 (0.00)

4 (5.04)

0 (0.00)

10 (12.88)

0 (0.00)

5 (6.68)

0 (0.00)

6 (8.28)

0 (0.00)

7-

3 (3.66)

0 (0.00)

8 (9.95)

0 (0.00)

9 (11.44)

0 (0.00)

6 (7.69)

0 (0.00)

8 (10.22)

0 (0.00)

8-

4 (5.10)

0 (0.00)

5 (6.11)

0 (0.00)

11 (13.75)

0 (0.00)

15 (16.25)

0 (0.00)

8 (7.34)

0 (0.00)

9-

6 (9.49)

0 (0.00)

7 (8.90)

0 (0.00)

6 (7.69)

0 (0.00)

16 (19.60)

0 (0.00)

11 (12.73)

0 (0.00)

10-

72 (16.16)

0 (0.00)

82 (19.68)

0 (0.00)

107 (25.96)

0 (0.00)

147 (36.70)

0 (0.00)

116 (30.73)

0 (0.00)

15-

320 (62.73)

2 (0.39)

346 (69.38)

0 (0.00)

396 (79.19)

0 (0.00)

570 (118.71)

0 (0.00)

497 (106.16)

0 (0.00)

20-

444 (97.99)

0 (0.00)

492 (104.25)

0 (0.00)

506 (107.63)

2 (0.43)

662 (131.56)

0 (0.00)

577 (111.56)

0 (0.00)

25-

432 (97.20)

1 (0.23)

438 (101.70)

0 (0.00)

490 (112.60)

0 (0.00)

511 (124.20)

1 (0.24)

467 (118.75)

0 (0.00)

30-

495 (99.50)

1 (0.20)

452 (94.60)

0 (0.00)

464 (95.63)

0 (0.00)

493 (105.84)

0 (0.00)

486 (110.15)

0 (0.00)

35-

514 (91.63)

2 (0.36)

475 (84.59)

0 (0.00)

472 (81.85)

0 (0.00)

536 (95.05)

0 (0.00)

526 (94.52)

0 (0.00)

40-

465 (89.57)

0 (0.00)

451 (86.40)

0 (0.00)

471 (87.35)

0 (0.00)

533 (94.34)

0 (0.00)

501 (87.92)

1 (0.18)

45-

365 (120.16)

3 (0.99)

338 (92.96)

0 (0.00)

391 (102.58)

0 (0.00)

475 (94.71)

1 (0.20)

497 (83.64)

2 (0.34)

50-

308 (111.30)

2 (0.72)

236 (90.25)

1 (0.38)

205 (75.55)

1 (0.37)

323 (136.33)

0 (0.00)

329 (141.60)

0 (0.00)

55-

317 (134.66)

2 (0.85)

314 (129.21)

0 (0.00)

335 (134.05)

1 (0.40)

396 (164.15)

0 (0.00)

378 (151.86)

1 (0.40)

60-

356 (202.39)

0 (0.00)

345 (189.88)

1 (0.55)

376 (201.76)

1 (0.54)

438 (242.60)

1 (0.55)

473 (259.34)

2 (1.10)

65-

379 (257.16)

1 (0.68)

316 (207.94)

1 (0.66)

336 (226.41)

0 (0.00)

448 (306.43)

0 (0.00)

432 (293.55)

0 (0.00)

70-

341 (306.56)

5 (4.50)

292 (259.21)

1 (0.89)

320 (290.80)

4 (3.64)

399 (355.38)

1 (0.89)

350 (304.40)

0 (0.00)

75-

201 (326.43)

2 (3.25)

170 (254.85)

1 (1.50)

195 (299.84)

1 (1.54)

218 (314.61)

2 (2.89)

239 (319.44)

1 (1.34)

80-

75 (292.16)

0 (0.00)

53 (204.76)

0 (0.00)

79 (312.09)

0 (0.00)

86 (314.02)

0 (0.00)

85 (298.13)

0 (0.00)

85-

14 (251.57)

0 (0.00)

17 (291.75)

0 (0.00)

15 (262.97)

0 (0.00)

31 (254.86)

1 (8.50)

28 (232.00)

0 (0.00)

Sex

Men

3 179 (111.66)

7 (0.25)

2 952 (103.18)

2 (0.07)

3 123 (107.18)

7 (0.24)

3 787 (131.70)

7 (0.24)

3 547 (120.26)

3 (0.10)

Women

1 962 (72.77)

14 (0.52)

1 916 (70.65)

4 (0.15)

2 109 (77.74)

3 (0.11)

2 582 (92.01)

0 (0.00)

2 508 (90.14)

4 (0.14)

Occupation

Farmers and herdsmen

3 441 (-)

13 (-)

3 329 (-)

5 (-)

3 558 (-)

7 (-)

4 328 (-)

6 (-)

4 249 (-)

4 (-)

Student

354 (-)

1 (-)

383 (-)

0 (-)

475 (-)

0 (-)

685 (-)

0 (-)

592 (-)

0 (-)

Worker

333 (-)

1 (-)

317 (-)

0 (-)

302 (-)

0 (-)

264 (-)

0 (-)

228 (-)

0 (-)

Attendant

79 (-)

0 (-)

60 (-)

0 (-)

66 (-)

0 (-)

59 (-)

0 (-)

36 (-)

0 (-)

Teacher

60 (-)

1 (-)

34 (-)

0 (-)

36 (-)

0 (-)

39 (-)

0 (-)

30 (-)

0 (-)

Medical personnel

19 (-)

0 (-)

16 (-)

0 (-)

29 (-)

0 (-)

29 (-)

0 (-)

23 (-)

0 (-)

Unemployed and retirees

855 (-)

5 (-)

729 (-)

1 (-)

766 (-)

3 (-)

965 (-)

1 (-)

897 (-)

3 (-)

Diagnostic category

Sputum smear positive

2 710 (-)

12 (-)

2 633 (-)

5 (-)

2 670 (-)

5 (-)

2 665 (-)

6 (-)

2 326 (-)

5 (-)

Bacterium negative

1 431 (-)

6 (-)

1 328 (-)

1 (-)

1 582 (-)

4 (-)

2 164 (-)

0 (-)

2 374 (-)

2 (-)

No detection in sputum

974 (-)

3 (-)

900 (-)

0 (-)

958 (-)

1 (-)

1 528 (-)

1 (-)

1 333 (-)

0 (-)

Only germiculture positive

26 (-)

0 (-)

7 (-)

0 (-)

22 (-)

0 (-)

12 (-)

0 (-)

22 (-)

0 (-)

Fig. 2

Annual incidence of tuberculosis in Qinghai Province, China, 2009-2013. The areas of high annual incidence of TB were mainly concentrated in south-western Qinghai, with the top three counties being Maduo, Jiuzhi, and Zaduo

Fig. 3

Monthly incidence rates of TB in Qinghai Province, China, from January 2009 to December 2013. The TB incidence rate showed significant periodicity and seasonality, reaching a seasonal peak around April and decreasing to a trough in December. TB, tuberculosis

Table 2

Descriptive statistics for meteorological variables in Qinghai Province, China, from July 2008 to December 2013

District, city

Meteorological variables

Mean

Standard deviation

Minimum

Maximum

Total

MAT (°C)

5.0

8.9

-13.4

20.4

MP (mm)

35.3

39.4

0.0

195.1

MSH (hours)

218.9

31.8

99.6

307.9

MAWS (m/s)

1.5

0.5

0.1

3.5

Xining

MAT (°C)

6.1

9.1

-11.2

18.8

MP (mm)

35.4

36.8

0.0

147.7

MSH (hours)

217.7

29.5

151.9

289.6

MAWS (m/s)

1.0

0.2

0.7

1.8

Haidong

MAT (°C)

7.8

9.1

-9.9

20.4

MP (mm)

28.6

30.4

0.0

128.9

MSH (hours)

222.2

27.9

156.6

282.2

MAWS (m/s)

2.0

0.4

1.0

2.6

Hainan

MAT (°C)

5.6

8.8

-9.8

18.8

MP (mm)

28.0

31.3

0.0

101.4

MSH (hours)

238.3

28.1

176.6

303.8

MAWS (m/s)

1.6

0.5

0.7

3.5

Haibei

MAT (°C)

2.1

9.1

-13.4

14.9

MP (mm)

41.0

42.7

0.0

195.1

MSH (hours)

206.2

24.4

133.1

262.3

MAWS (m/s)

1.5

0.4

0.6

2.5

Haixi

MAT (°C)

5.0

9.7

-11.8

19.5

MP (mm)

21.0

30.0

0.0

124.2

MSH (hours)

243.3

28.1

192.1

307.9

MAWS (m/s)

1.6

0.5

0.9

2.6

Huangnan

MAT (°C)

7.1

8.3

-8.7

19.2

MP (mm)

36.1

37.7

0.0

144.2

MSH (hours)

209.7

30.0

121.3

278.6

MAWS (m/s)

1.1

0.3

0.1

1.6

Guoluo

MAT (°C)

0.9

7.8

-12.8

12.4

MP (mm)

47.3

50.7

0.0

169.5

MSH (hours)

210.8

34.6

99.6

285.5

MAWS (m/s)

1.8

0.4

1.2

2.7

Yushu

MAT (°C)

5.1

7.8

-7.7

19.2

MP (mm)

44.9

44.8

0.0

159.8

MSH (hours)

202.9

26.3

132.9

274.1

MAWS (m/s)

1.3

0.4

0.8

2.5

MAT monthly average temperature, MP monthly precipitation, MSH monthly total sunshine hours, MAWS monthly average wind speed

Spatial autocorrelation analysis of TB incidence

The global Moran’s I values of each year at the county level were high, ranging from 0.40 to 0.58 (all P-values < 0.05), which indicated that the counties with high TB incidence tended to be adjacent to the districts with high TB incidence, and that the counties with low TB incidence tended to be adjacent to the districts with low TB incidence (Fig. 4). LISA analysis revealed 11 counties with a significantly high-high spatial clustering and 17 counties with a significantly low-low spatial clustering of TB annual incidence in the five-year period. The high-high clustering areas were mainly concentrated around the cities of Yushu and Guoluo, such as Chengduo, Maduo, Qumalai, and Dari counties, with an average annual incidence of 294/100 000. The low-low clustering areas were concentrated in Xining city and surrounding areas, as well as in several counties of Haixi city, such as Mangya, Lenghu, and Dachaidan, with an average incidence of 68/100 000 (Fig. 5).
Fig. 4

Moran scatter plot for the annual incidence of tuberculosis in Qinghai Province, China, 2009-2013. The horizontal axis shows the standardized incidence of the counties, and the vertical axis indicates the spatial lag factors; the linear slope is the Moran’s I

Fig. 5

LISA significance map and cluster map for annual tuberculosis incidence in Qinghai Province, China, 2009-2013. The high risk areas were mainly concentrated in the cities of Yushu and Guoluo, while the low incidence districts were mainly distributed in the cities of Xining and Haixi. LISA, local indicators of spatial association

Spatial panel analysis of meteorological factors

TB is a chronic infectious disease with a certain amount of time lag between the influencing factors and the disease. We used the meteorological factors with a 0- to 6-month lag from the incidence of TB to fit the simple linear model and classical panel data models. The associations between TB incidence and meteorological factors with a 3-month lag were found to have the best goodness of fit. Subsequently, we examined the individual effects of spatial cities by using the Hausman-test and F-test; as a result, a significant fixed effect of each city was found (P < 0.001). The Durbin-Watson statistic and Moran’s I (P < 0.001) indicated a spatial autocorrelation in error term. The Lagrange multiplier (LM) test showed that the spatial lag effect was more significant than the spatial error effect (Table 3). Therefore, we finally used the spatial lag fixed effects panel data model to examine the associations between TB incidence and meteorological factors.
Table 3

Results of the classical panel data models for the log of TB incidence with meteorological factors of a 3-month lag

Factors

Simple linear regression

Fixed effects model

Random effects model

Coefficient

t-value

P-value

Coefficient

t-value

P-value

Coefficient

t-value

P-value

Constant

1.0097

10.42

<0.0001

   

1.0367

12.27

<0.0001

MAT (°C)

-0.0195

-10.03

<0.0001

-0.0059

-4.96

<0.0001

-0.0060

-5.12

<0.0001

MP (mm)

0.0024

5.07

<0.0001

-0.0008

-2.77

0.0060

-0.0007

-2.63

0.0090

MSH (hours)

-0.0011

-2.56

0.0110

-0.0003

-1.10

0.2720

-0.0003

-1.14

0.2570

MAWS (m/s)

0.1441

5.94

<0.0001

0.0320

1.93

0.0550

0.0339

2.05

0.0410

Log likelihood

1.92

  

328.57

     

AIC

0.01

  

-1.32

     

SC

0.06

  

-1.21

     

F-statistic

193.90

 

<0.0001

      

H-statistic

      

10.41

 

0.0340

LM lag

14.14

 

<0.0001

39.40

 

<0.0001

   

Robust LM lag

161.22

 

<0.0001

21.00

 

<0.0001

   

LM error

34.13

 

<0.0001

35.93

 

<0.0001

   

Robust LM error

181.21

 

<0.0001

17.54

 

<0.0001

   

Moran’s I

0.20

6.07

<0.0001

0.20

6.20

<0.0001

   

TB tuberculosis, MAT monthly average temperature, MP monthly precipitation, MSH monthly total sunshine hours, MAWS monthly average wind speed, AIC Akaike information criterion, SC Schwarz Criterion, H-statistic Hausman-statistic, LM Lagrange multiplier

The result showed that the background incidence of each city was different; the highest and lowest background incidences were found in Guoluo city (12.88/100 000) and Haixi city (2.75/100 000), respectively (Table 4). The spatial autocorrelation coefficient was 0.30, meaning that a spatial spillover phenomenon existed; the TB incidence of the study city increased 1 time when the TB incidence of adjacent cities increased 9 times and other influencing factors were kept constant. The regression coefficients of MAT, MP, and MAWS were all statistically significant (all P-values < 0.05), with each 10 °C, 2 cm, and 1 m/s increase in temperature, precipitation, and wind speed being associated with 9 % and 3 % decrements and a 7 % increment in the TB incidence, respectively (Table 5).
Table 4

Results for spatial individual effects of each city by using the spatial panel data model

Cities

Intercept term (μ i )

Background incidence (1/100 000)

Cities

Intercept term (μ i )

Background incidence (1/100 000)

Haixi

0.44

2.75

Hainan

0.70

5.01

Xining

0.51

3.24

Haibei

0.79

6.17

Haidong

0.57

3.72

Yushu

1.09

12.30

Huangnan

0.64

4.37

Guoluo

1.11

12.88

Table 5

Results of the spatial panel data model for the log of TB incidence with meteorological factors of a 3-month lag

Factors

Coefficient

95 % CIs of coefficients

P-value

MAT (°C)

-0.0040

-0.0063, -0.0017

<0.0001

MP (mm)

-0.0006

-0.0011, -0.0001

0.0250

MSH (hours)

-0.0002

-0.0007, 0.0002

0.2680

MAWS (m/s)

0.0309

0.0001, 0.0618

0.0490

ρ

0.2997

0.2003, 0.3992

<0.0001

TB tuberculosis, CI confidence interval, MAT monthly average temperature, MP monthly precipitation, MSH monthly total sunshine hours, MAWS monthly average wind speed, ρ spatial autocorrelation coefficient

Discussion

In our study, we used the data of TB from the CISDCP to analyse the characteristics of the TB incidence in Qinghai province. The incidence of TB in Qinghai province showed clear periodic and seasonal features in the colder winter and spring months, similar to the patterns of epidemic regularity reported in Beijing [5] and Shandong [2], China, and in India [36, 37] and Mongolia [38]. In our study, it was shown that the TB incidence was higher in elders than in young people, with farmers and herdsmen at particularly high risks (Table 1), which supported the notion of TB pathogenesis regularity [39]. Moreover, the TB incidence in regions with relatively poor economic conditions and high altitudes has been shown to be high, potentially owning to differences in the ethnicity, medical and health conditions, and economic and educational levels of the residents. Additionally, the local environment and climatic conditions may also influence the incidence of TB [37, 4043]. Thus, in Qinghai, the farmers and herdsmen living in the regions with cold climate and high altitudes should be made aware of the high risk of TB occurrence.

The incidence and spatial autocorrelation differ, but are associated in some regards. A higher incidence reflects higher epidemic strength. Spatial autocorrelation describes the relationship of incidence between one region and the surrounding areas and is based on incidence and overall consideration of regional geographical, human, and environmental factors. In our study, we chose Moran’s I to analyse the spatial autocorrelation. It showed a positive correlation within regions, and the correlation displayed an increasing tendency year by year, indicating that the distribution of TB in Qinghai is not random, with obvious spatial clusters. These results are consistent with previously published studies [2, 11, 37, 40].

Spatial clustering analysis has suggested that classical multiple linear regression analysis is not suitable for exploration of TB risk factors at the ecological level. Thus, this method should be combined with a spatial statistical model to explore the risk factors [25]. A panel data model can be directly used to assess the differences between regions, control for individual heterogeneity, and present more reasonable results. However, this model does not consider the spatial autocorrelation in the study area. Accordingly, a spatial panel data model combining spatial metrology with the panel data model was used in our study. This model takes the individual effects and spatial autocorrelation into consideration, which could result in better use of the spatiotemporal information from infectious diseases surveillance data [24, 26]. In a spatial panel data model, the introduction of a spatial individual effect can correct deviations caused by unobserved variables, and the use of a spatial weight matrix can illustrate the spatial correlation and reflect the interaction between regions [34]. Moreover, the modelling estimation is more effective and could help us research unknown variables in depth.

Some studies have reported that the average incubation period of TB ranges from four to eight weeks, with a two-month interval from the symptom appearance to medical diagnosis [6]. Accordingly, we created a fitting model with a 3-month lag [4446]. The results showed that the combined data using the traditional method would overstate the influence of the meteorological factors, because the effect is different in different regions. Conversely, the spatial panel data model reduced error occurrence and increased the goodness of fit. The spatial autocorrelation coefficient was 0.30, which indicated that the introduction of the spatial lag dependent variable could reasonably explain the spatial autocorrelation. Further, the regression coefficients of MAT and MP were negative, while that of MAWS was positive, suggesting that with increasing MAT and MP, the incidence of TB decreases exponentially. Conversely, with increasing MAWS, the incidence increases exponentially. These result are in accordance with the findings of previous studies, such as those by Li et al. [6] and Naranbat et al. [38].

In cold condition, especially in the winter, most people stay indoors for a long time, hence, once someone releases bacteria into the environment, elderly and immunocompromised populations are at particularly high risk of infection because of the poor ventilation. Moreover, some studies have shown a relationship between vitamin D levels and TB incidence [38]. Fewer outdoor activities and less exposure time to sunlight result in decreased synthesis of vitamin D in the human body. In turn, this may increase the risk of TB infection [6, 47]. In the traditional view, TB infection is mainly caused by bacteria in the sputum. When the sputum gets dry, bacteria are expelled into the air and enter other people’s respiratory systems. With increasing rainfalls, the air humidity increases and the transmission of aerosols decreases. Hence, the risk of bacteria entering the respiratory system is decreased, as is the incidence of TB.

In summary, we found that the spread of TB in Qinghai province is not random, but rather present obvious spatial clustering. The use of spatial statistic methods may offer necessary feedback in terms of the prevalence nature and epidemic characteristics of TB in various regions, and may consequently result in public health officials providing TB control and developing novel prevention strategies. In this study, we quantified the relationship between TB and meteorological factors using a spatial panel data model on the basis of panel data for the first time. A spatial panel data model is appropriate if longitudinal data of multiple units are available and if spatial autocorrelation exists. This model has a better model fitting and provides more precise effect size estimation. Additionally, as meteorological factors obviously affect the TB incidence in Qinghai province, future strategies of TB control and prevention should consider climate variations.

However, there are some limitations in the present study that are worth mentioning. First, as it is difficult to collect meteorological data of each county, our analysis was initiated at the city level. Second, the temperature, precipitation, sunshine hours, and wind speed are not the only meteorological factors affecting TB distribution. Therefore, additional county-level factors, such as atmospheric pressure, average vapour pressure, and average relative humidity, will be taken into consideration in our future study.

Conclusions

Our study found that high-high clustering areas of TB incidence were mainly concentrated in the southwest, while low-low clustering areas were found mainly in eastern and north-western Qinghai. The TB incidence was positively associated with the temperature, precipitation, and wind speed after adjusting for spatial heterogeneity and spatial correlation in Qinghai province.

Ethics and consent

The Ethics Review Board of the Qinghai Center for Disease Control and Prevention was consulted, since the project used only surveillance data, confirmed that ethics approval was not required because the data did not contain any personal or health information that could be connected back to the original identifiers.

Abbreviations

LISA: 

local indicators of spatial association

MAT: 

monthly average temperature

MAWS: 

monthly average wind speed

MP: 

monthly precipitation

MSH: 

monthly total sunshine hours

TB: 

tuberculosis

Declarations

Acknowledgements

The authors have no external support of funding to report. This study was supported by the Qinghai Center for Disease Control and Prevention (CDC). The authors gratefully acknowledge the staff involved in TB surveillance at all participating levels in Qinghai province, China.

Open AccessThis 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.

Authors’ Affiliations

(1)
Department of Health Statistics, School of Public Health, Shanxi Medical University
(2)
Department of Epidemiology, Richard M. Fairbanks School of Public Health, Indiana University
(3)
Department of Biochemistry and Molecular Biology, School of Basic Medical Sciences, Shanxi Medical University
(4)
Institute for Communicable Disease Control and Prevention, Qinghai Center for Disease Control and Prevention

References

  1. Addis Z, Adem E, Alemu A, Birhan W, Mathewos B, Tachebele B, et al. Prevalence of smear positive pulmonary tuberculosis in Gondar prisoners, North West Ethiopia. Asian Pac J Trop Med. 2015;8(2):127–31. doi:10.1016/S1995-7645(14)60302-3.View ArticlePubMedGoogle Scholar
  2. Wang T, Xue F, Chen Y, Ma Y, Liu Y. The spatial epidemiology of tuberculosis in Linyi City, China, 2005-2010. BMC Public Health. 2012;12:885. doi:10.1186/1471-2458-12-885.View ArticlePubMedPubMed CentralGoogle Scholar
  3. Hu H, Chen J, Sato KD, Zhou Y, Jiang H, Wu P, et al. Factors that associated with TB patient admission rate and TB inpatient service cost: a cross-sectional study in China. Infect Dis Poverty. 2016;5(1):4. doi:10.1186/s40249-016-0097-x.View ArticlePubMedPubMed CentralGoogle Scholar
  4. Tang S. Tackling challenges of TB/MDRTB in China: concerted actions are imperative. Infect Dis Poverty. 2015;4:19. doi:10.1186/s40249-015-0050-4.View ArticlePubMedPubMed CentralGoogle Scholar
  5. Liu Y, Li X, Wang W, Li Z, Hou M, He Y, et al. Investigation of space-time clusters and geospatial hot spots for the occurrence of tuberculosis in Beijing. Int J Tuberc Lung Dis. 2012;16(4):486–91. doi:10.5588/ijtld.11.0255.View ArticlePubMedGoogle Scholar
  6. Li XX, Wang LX, Zhang H, Du X, Jiang SW, Shen T, et al. Seasonal variations in notification of active tuberculosis cases in China, 2005-2012. PLoS One. 2013;8(7):e68102. doi:10.1371/journal.pone.0068102.View ArticlePubMedPubMed CentralGoogle Scholar
  7. Ben-Selma W, Harizi H, Letaief M, Boukadida J. Age- and gender-specific effects on NRAMP1 gene polymorphisms and risk of the development of active tuberculosis in Tunisian populations. Int J Infect Dis. 2012;16(7):e543–50. doi:10.1016/j.ijid.2011.11.016.View ArticlePubMedGoogle Scholar
  8. Rhines AS. The role of sex differences in the prevalence and transmission of tuberculosis. Tuberculosis (Edinb). 2013;93(1):104–7. doi:10.1016/j.tube.2012.10.012.View ArticleGoogle Scholar
  9. Cois A, Ehrlich R. Problem drinking as a risk factor for tuberculosis: a propensity score matched analysis of a national survey. BMC Public Health. 2013;13:871. doi:10.1186/1471-2458-13-871.View ArticlePubMedPubMed CentralGoogle Scholar
  10. Singh PN, Yel D, Kheam T, Hurd G, Job JS. Cigarette smoking and tuberculosis in Cambodia: findings from a national sample. Tob Induc Dis. 2013;11(1):8. doi:10.1186/1617-9625-11-8.View ArticlePubMedPubMed CentralGoogle Scholar
  11. Sun W, Gong J, Zhou J, Zhao Y, Tan J, Ibrahim AN, et al. A spatial, social and environmental study of tuberculosis in China using statistical and GIS technology. Int J Environ Res Public Health. 2015;12(2):1425–48. doi:10.3390/ijerph120201425.View ArticlePubMedPubMed CentralGoogle Scholar
  12. Li XX, Wang LX, Zhang J, Liu YX, Zhang H, Jiang SW, et al. Exploration of ecological factors related to the spatial heterogeneity of tuberculosis prevalence in P. R. China. Glob Health Action. 2014;7:23620. doi:10.3402/gha.v7.23620.PubMedGoogle Scholar
  13. Ukwaja KN, Alobu I, Abimbola S, Hopewell PC. Household catastrophic payments for tuberculosis care in Nigeria: incidence, determinants, and policy implications for universal health coverage. Infect Dis Poverty. 2013;2(1):21. doi:10.1186/2049-9957-2-21.View ArticlePubMedPubMed CentralGoogle Scholar
  14. Zou B, Peng F, Wan N, Mamady K, Wilson GJ. Spatial cluster detection of air pollution exposure inequities across the United States. PLoS One. 2014;9(3):e91917. doi:10.1371/journal.pone.0091917.View ArticlePubMedPubMed CentralGoogle Scholar
  15. Varga C, Pearl DL, McEwen SA, Sargeant JM, Pollari F, Guerin MT. Area-level global and local clustering of human Salmonella Enteritidis infection rates in the city of Toronto, Canada, 2007-2009. BMC Infect Dis. 2015;15:359. doi:10.1186/s12879-015-1106-6.View ArticlePubMedPubMed CentralGoogle Scholar
  16. Abbas T, Younus M, Muhammad SA. Spatial cluster analysis of human cases of Crimean Congo hemorrhagic fever reported in Pakistan. Infect Dis Poverty. 2015;4:9. doi:10.1186/2049-9957-4-9.View ArticlePubMedPubMed CentralGoogle Scholar
  17. Selemani M, Mrema S, Shamte A, Shabani J, Mahande MJ, Yeates K, et al. Spatial and space- time clustering of mortality due to malaria in rural Tanzania: evidence from Ifakara and Rufiji Health and Demographic Surveillance System sites. Malar J. 2015;14:369. doi:10.1186/s12936-015-0905-y.View ArticlePubMedPubMed CentralGoogle Scholar
  18. Saman DM, Cole HP, Odoi A, Myers ML, Carey DI, Westneat SC. A spatial cluster analysis of tractor overturns in Kentucky from 1960 to 2002. PLoS One. 2012;7(1):e30532. doi:10.1371/journal.pone.0030532.View ArticlePubMedPubMed CentralGoogle Scholar
  19. Ngamini Ngui A, Apparicio P, Moltchanova E, Vasiliadis HM. Spatial analysis of suicide mortality in Quebec: spatial clustering and area factor correlates. Psychiatry Res. 2014;220(1-2):20–30. doi:10.1016/j.psychres.2014.07.033.View ArticlePubMedGoogle Scholar
  20. 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. doi:10.1017/S0031182013001972.View ArticlePubMedGoogle Scholar
  21. Zulu LC, Kalipeni E, Johannes E. Analyzing spatial clustering and the spatiotemporal nature and trends of HIV/AIDS prevalence using GIS: the case of Malawi, 1994-2010. BMC Infect Dis. 2014;14:285. doi:10.1186/1471-2334-14-285.View ArticlePubMedPubMed CentralGoogle Scholar
  22. Lee SS, Wong NS. The clustering and transmission dynamics of pandemic influenza A (H1N1) 2009 cases in Hong Kong. J Infect. 2011;63(4):274–80. doi:10.1016/j.jinf.2011.03.011.View ArticlePubMedGoogle Scholar
  23. Wubuli A, Xue F, Jiang D, Yao X, Upur H, Wushouer Q. Socio-Demographic Predictors and Distribution of Pulmonary Tuberculosis (TB) in Xinjiang, China: A Spatial Analysis. PLoS One. 2015;10(12):e0144010. doi:10.1371/journal.pone.0144010.View ArticlePubMedPubMed CentralGoogle Scholar
  24. Wang H, Du Z, Wang X, Liu Y, Yuan Z, Liu Y, et al. Detecting the association between meteorological factors and hand, foot, and mouth disease using spatial panel data models. Int J Infect Dis. 2015;34:66–70. doi:10.1016/j.ijid.2015.03.007.View ArticlePubMedGoogle Scholar
  25. Ehlert A, Oberschachtsiek D. Does managed care reduce health care expenditure? Evidence from spatial panel data. Int J Health Care Finance Econ. 2014;14(3):207–27. doi:10.1007/s10754-014-9145-x.View ArticlePubMedGoogle Scholar
  26. Barufi AM, Haddad E, Paez A. Infant mortality in Brazil, 1980-2000: a spatial panel data analysis. BMC Public Health. 2012;12:181. doi:10.1186/1471-2458-12-181.View ArticlePubMedPubMed CentralGoogle Scholar
  27. Wu W, Guo J, Guan P, Sun Y, Zhou B. Clusters of spatial, temporal, and space-time distribution of hemorrhagic fever with renal syndrome in Liaoning Province, Northeastern China. BMC Infect Dis. 2011;11:229. doi:10.1186/1471-2334-11-229.View ArticlePubMedPubMed CentralGoogle Scholar
  28. Chen Y. New approaches for calculating Moran’s index of spatial autocorrelation. PLoS One. 2013;8(7):e68336. doi:10.1371/journal.pone.0068336.View ArticlePubMedPubMed CentralGoogle Scholar
  29. Xia J, Cai S, Zhang H, Lin W, Fan Y, Qiu J, et al. Spatial, temporal, and spatiotemporal analysis of malaria in Hubei Province, China from 2004-2011. Malar J. 2015;14:145. doi:10.1186/s12936-015-0650-2.View ArticlePubMedPubMed CentralGoogle Scholar
  30. Yazdani-Charati J, Siamian H, Kazemnejad A, Mohammad V. Spatial clustering of tuberculosis incidence in the North of Iran. Glob J Health Sci. 2014;6(6):288–94. doi:10.5539/gjhs.v6n6p288.PubMedPubMed CentralGoogle Scholar
  31. Li Y, Zhou YZ, Zhang CB, Dou L, Du HY, Lin XM, et al. Application of local Moran’s I and GIS to identify hotspots of Ni, Cr of vegetable soils in high-incidence area of liver cancer from the Pearl River Delta, South China. Huan Jing Ke Xue (in Chinese). 2010;31(6):1617–23.Google Scholar
  32. Wei Y, Huang Y, Luo L, Xiao X, Liu L, Yang Z. Rapid increase of scrub typhus: an epidemiology and spatial-temporal cluster analysis in Guangzhou City, Southern China, 2006–2012. PLoS One. 2014;9(7):e101976. doi:10.1371/journal.pone.0101976.View ArticlePubMedPubMed CentralGoogle Scholar
  33. Kelejian HH, Prucha IR. Specification and Estimation of Spatial Autoregressive Models with Autoregressive and Heteroskedastic Disturbances. J Econom. 2010;157(1):53–67. doi:10.1016/j.jeconom.2009.10.025.View ArticlePubMedPubMed CentralGoogle Scholar
  34. Ponicki WR, Gruenewald PJ, Remer LG. Spatial panel analyses of alcohol outlets and motor vehicle crashes in California: 1999-2008. Accid Anal Prev. 2013;55:135–43. doi:10.1016/j.aap.2013.03.001.View ArticlePubMedPubMed CentralGoogle Scholar
  35. Leonard WR, Reyes-Garcia V, Tanner S, Rosinger A, Schultz A, Vadez V, et al. The Tsimane’ Amazonian Panel Study (TAPS): Nine years (2002-2010) of annual data available to the public. Econ Hum Biol. 2015;19:51–61. doi:10.1016/j.ehb.2015.07.004.View ArticlePubMedGoogle Scholar
  36. Kumar V, Singh A, Adhikary M, Daral S, Khokhar A, Singh S. Seasonality of tuberculosis in delhi, India: a time series analysis. Tuberc Res Treat. 2014;2014:514093. doi:10.1155/2014/514093.PubMedPubMed CentralGoogle Scholar
  37. Narula P, Sihota P, Azad S, Lio P. Analyzing seasonality of tuberculosis across Indian states and union territories. J Epidemiol Glob Health. 2015;5(4):337–46. doi:10.1016/j.jegh.2015.02.004.View ArticlePubMedGoogle Scholar
  38. Naranbat N, Nymadawa P, Schopfer K, Rieder HL. Seasonality of tuberculosis in an Eastern- Asian country with an extreme continental climate. Eur Respir J. 2009;34(4):921–5. doi:10.1183/09031936.00035309.View ArticlePubMedGoogle Scholar
  39. Huynh GH, Klein DJ, Chin DP, Wagner BG, Eckhoff PA, Liu R, et al. Tuberculosis control strategies to reach the 2035 global targets in China: the role of changing demographics and reactivation disease. BMC Med. 2015;13:88. doi:10.1186/s12916-015-0341-4.View ArticlePubMedPubMed CentralGoogle Scholar
  40. Harling G, Castro MC. A spatial analysis of social and economic determinants of tuberculosis in Brazil. Health Place. 2014;25:56–67. doi:10.1016/j.healthplace.2013.10.008.View ArticlePubMedGoogle Scholar
  41. Prussing C, Castillo-Salgado C, Baruch N, Cronin WA. Geo-epidemiologic and molecular characterization to identify social, cultural, and economic factors where targeted tuberculosis control activities can reduce incidence in Maryland, 2004-2010. Public Health Rep. 2013;128 Suppl 3:104–14.PubMedPubMed CentralGoogle Scholar
  42. Scott C, Kirking HL, Jeffries C, Price SF, Pratt R. Tuberculosis trends--United States, 2014. MMWR Morb Mortal Wkly Rep. 2015;64(10):265–9.PubMedGoogle Scholar
  43. Mamaev IA, Khachirov Dzh G, Guseinov GK. Tuberculosis morbidity in a male population by the climatic, geographic and administrative rural areas of Dagestan. Probl Tuberk (in Russian). 2000;2:13–7.Google Scholar
  44. Tattevin P, Che D, Fraisse P, Gatey C, Guichard C, Antoine D, et al. Factors associated with patient and health care system delay in the diagnosis of tuberculosis in France. Int J Tuberc Lung Dis. 2012;16(4):510–5. doi:10.5588/ijtld.11.0420.View ArticlePubMedGoogle Scholar
  45. Cheng S, Chen W, Yang Y, Chu P, Liu X, Zhao M, et al. Effect of Diagnostic and Treatment Delay on the Risk of Tuberculosis Transmission in Shenzhen, China: An Observational Cohort Study, 1993-2010. PLoS One. 2013;8(6):e67516. doi:10.1371/journal.pone.0067516.View ArticlePubMedPubMed CentralGoogle Scholar
  46. Belay M, Bjune G, Ameni G, Abebe F. Diagnostic and treatment delay among Tuberculosis patients in Afar Region, Ethiopia: a cross-sectional study. BMC Public Health. 2012;12:369. doi:10.1186/1471-2458-12-369.View ArticlePubMedPubMed CentralGoogle Scholar
  47. Fares A. Seasonality of tuberculosis. J Glob Infect Dis. 2011;3(1):46–55. doi:10.4103/0974-777X.77296.View ArticlePubMedPubMed CentralGoogle Scholar

Copyright

© Rao et al. 2016

Advertisement