Skip to main content

The characteristics of spatial-temporal distribution and cluster of tuberculosis in Yunnan Province, China, 2005–2018

Abstract

Background

Tuberculosis (TB) makes a big challenge to public health, especially in high TB burden counties of China and Greater Mekong Subregion (GMS). The aim of this study was to identify the spatial-temporal dynamic process and high-risk region of notified pulmonary tuberculosis (PTB), sputum smear-positive tuberculosis (SSP-TB) and sputum smear-negative tuberculosis (SSN-TB) cases in Yunnan, the south-western of China between years of 2005 to 2018. Meanwhile, to evaluate the similarity of prevalence pattern for TB among GMS.

Methods

Data for notified PTB were extracted from the China Information System for Disease Control and Prevention (CISDCP) correspond to population information in 129 counties of Yunnan between 2005 to 2018. Seasonally adjusted time series defined the trend cycle and seasonality of PTB prevalence. Kulldorff’s space-time scan statistics was applied to identify temporal, spatial and spatial-temporal PTB prevalence clusters at county-level of Yunnan. Pearson correlation coefficient and hierarchical clustering were applied to define the similarity of TB prevalence among borders with GMS.

Result

There were a total of 381,855 notified PTB cases in Yunnan, and the average prevalence was 59.1 per 100,000 population between 2005 to 2018. A declined long-term trend with seasonality of a peak in spring and a trough in winter for PTB was observed. Spatial-temporal scan statistics detected the significant clusters of PTB prevalence, the most likely cluster concentrated in the northeastern angle of Yunnan between 2011 to 2015 (RR = 2.6, P < 0.01), though the most recent cluster for PTB and spatial cluster for SSP-TB was in borders with GMS. There were six potential TB prevalence patterns among GMS.

Conclusion

This study detected aggregated time interval and regions for PTB, SSP-TB, and SSN-TB at county-level of Yunnan province. Similarity prevalence pattern was found in borders and GMS. The localized prevention strategy should focus on cross-boundary transmission and SSN-TB control.

Peer Review reports

Background

Tuberculosis (TB) is a communicable disease caused by the agent of Mycobacterium tuberculosis, which mainly invade lung tissues and lead to airborne infectious disease of pulmonary tuberculosis (PTB).

TB made a big challenge to public health especially in high disease burden counties [1, 2]. Globally, high TB burden counties were primarily concentrated in Asia and Africa. World Health Organization (WHO) estimated there were 10.0 million TB cases worldwide, TB incidence was 889,000 and a rate of 63 per 100,000 population for China in 2017 [3]. Especially, Five out of six Greater Mekong Subregion (GMS) countries in southeast Asia were defined as high epidemiological TB burden, thus, Myanmar, Lao People’s Democratic Republic (Laos) and Vietnam shared not only national boundaries but also TB burden with China.

Although nation-wide longitudinal TB prevalence surveys showed that prevalence of PTB and the most infectious sputum smear-positive tuberculosis (SSP-TB) substantially decreased through two decades efforts and intensive directly observed treatment, short-course strategy (DOTS) programme [4], China still has the second-largest burden of TB cases in the world. Previous statistical modeling analysis by applying time series method reported that the trend for notified TB cases decreased in China between 2005 to 2012 [5], but long-term trend various in different provinces in the nation [6, 7]. Yunnan Province followed the National Tuberculosis Control Program in China (China NTP), yet implemented non-routine strategy of active case finding to detect TB cases in hard-to-reach population among few communities [8]. This public health action may have an impact on local TB epidemic status and made it a challenge to describe and understand the trend of TB prevalence.

As an airborne disease, PTB epidemics influenced the transmission in geographical neighborhoods, disease hotspots were defined as high-risk clusters. Studies reported hotspots regions or high-risk clusters of TB in China by using spatial-temporal scan analysis [9,10,11,12,13,14], spatial-temporal distribution characteristics were illustrated at the national, provincial, prefectural, county-level or individual level in diverse time frames, all of which have shown the geographical and temporal heterogeneity of TB epidemic.

The geospatial difference of TB prevalence distribution was observed for Yunnan in 2018, the highest notification rate of Lanping county was 28 times higher than the lowest rate of Tonghai county [15]. Yunnan province was in a unique geographical location, as one part of GMS and southwestern gateway of China, moreover, surrounded by GMS high TB burden counties, making it indispensable to clarify the spatial-temporal heterogeneous distribution of TB in the province.

Our aim of this study was, first, to detect TB spatial-temporal clusters at county-level of Yunnan; second, to understand the temporal trend of notified TB in Yunnan; in addition, to explore correlated pattern of TB prevalence among GMS and border counties of Yunnan. Knowledge of long-term trend and spatial-temporal distribution of TB prevalence was crucial to understand the dynamic transmission of TB and to provide local evidence of classified TB prevention and control strategies.

Methods

General setting

Yunnan locate in the far southwest corner of China, between longitude 97°31′ and 106°11′ East, latitude of 21°80′ and 29°15′ North, with a total area of 397,100 km2 while mountains occupy 84% of the territory, while the altitude various in tremendous range from 76 to 6740 m. There are 16 prefectures and 129 counties in Yunnan, the population of which was 48 million in 2018 [16]. Yunnan province compose GMS and adjoin Myanmar, Laos, and Vietnam by land (Fig. 1a). There are 25 border counties in Yunnan; 18, 2 and 7 counties of them adjoin with Myanmar, Laos, and Vietnam respectively. Mengla county adjoin by Myanmar and Laos, Jiangcheng county bordering on Laos and Vietnam simultaneously (Fig. 1b).

Fig. 1
figure1

Location of Yunnan Province and Greater Mekong Subregion with border counties, and the county-level average tuberculosis prevalence of Yunnan, 2005–2018. Location and prefectures of Yunnan Province and China, Myanmar, Laos and Vietnam (a), the border and inland counties of China with GMS counties (b), The average prevalence of PTB in Yunnan, 2005–2018 (c), the average prevalence of SSP-TB in Yunnan, 2005–2018 (d), the average prevalence of SSN-TB in Yunnan, 2005–2018 (e). GMS Greater Mekong Subregion, PTB Pulmonary tuberculosis, SSP-TB Sputum smear-positive tuberculosis, SSN-TB Sputum smear-negative tuberculosis. World geographic database were extracted from the GADM database (https://www.gadm.org/), China and Yunnan geographic database provided by National Geomatics Center of China (http://www.ngcc.cn/ngcc/) at a 1:1,000,000 scale as the layer’s attribute

Data source

TB diagnosis was based on patients’ suspicious symptoms plus results of Chest X-ray and sputum smear microscopy in TB designated hospital or county/district level Center for Disease Control and Prevention (CDC) TB clinic in Yunnan. Clinic diagnosed or laboratory confirmed PTB cases were notified in the China Information System for Disease Control and Prevention (CISDCP). Monthly TB cases notification data of 129 counties/districts in Yunnan from January 2005 to December 2018 were extracted from CISDCP. Annually demographic data of 129 counties/districts were collected from Yunnan statistical yearbooks from 2005 to 2018 [16]. The average TB notification rate was calculated by the numerator of summing county-level notified cases and the denominator of the summing population of counties in 14 years (Fig. 1 c, d, e). The raw counts and prevalence for PTB, SSP-TB, and SSN-TB (sputum smear-negative tuberculosis) cases were applied for analysis.

National-level estimated TB incidence and reported SSP-TB cases from 2005 to 2017 in GMS countries of Myanmar, Laos, and Vietnam were extracted from WHO’s global tuberculosis database (https://www.who.int/tb/data/en/).

Statistical methods

Time series analysis

Monthly reported TB cases counts were aggregated into provincial level then analyzed by X-12-ARIMA (autoregressive integrated moving average) seasonal adjustment [5, 17], X-12-ARIMA was developed by US Census Bureau to define seasonal adjustment for time series [18]. By applying the X-12-ARIMA process, series of notified TB cases were decomposed into three components: seasonal variation, long-term trend cycle, and random irregular noise. Seasonal factors between 2005 to 2018 were calculated to explore seasonality of TB.

Temporal, spatial and space-time scan statistic

Kulldorff’s scan statistics was applied to detect clusters of TB cases in either a purely temporal, purely spatial or space-time setting. Spatial-temporal scan method created the infinite number of scanning windows, scan statistics gradually scanning a window across time and/or space, recorded the number of observed and expected observations inside the window at each location [19]. The scanning window was a time interval in the purely temporal scan for one-dimensional line, a circle in the purely spatial scan for bidimensional surface, or a cylinder in spatial-temporal scan with a circular base and height of time period for three-dimensional space [20, 21].

The surveillance data of notified TB case number against the population at risk in the same county/district was presumptively followed the Poisson probability distribution. Our study was a retrospective space-time scan analysis based on the discrete Poisson model. Under the Poisson distribution assumption, for each location and size of scanning windows, the alternative hypothesis was that there was an elevated risk within the window as compared to outside. The likelihood function calculated to define clusters, likelihood function was maximized overall scanned windows, maximum log-likelihood ratio (LLR) correspond the most likely cluster, which means the least likely to have occurred by chance; meanwhile, other ordinal statistically significant LLRs were matched to secondary clusters. P-value of maximum likelihood test was obtained through Monte Carlo hypothesis testing by randomly replicated simulations tests for comparing the rank of the maximum LLR from real data with from random data. The relative risk (RR) defined with the risk within scanned window compared to risk outside the scan window, RR representing how much more common disease was in this location and time period compared to the baseline [21].

Spatial-temporal scan parameters were selected by following the principle of reducing geographical overlap in clusters [10, 13]. The length of scanning time window covered 30% of the entire study period, and the scanning space window was set to cover 13% of the population at risk. Spatial scan for PTB prevalence in each year of study was applied to detect dynamic of geo-clusters, then scan spatial-temporal clusters for PTB, SSP-TB and SSN-TB prevalence by aggregating data of 14 years with the same parameters. The limit number of Monte Carlo replications was set to 999 times to detect the most likely and secondary clusters.

TB epidemic intensity in spatial-temporal clusters was measured by monthly average TB prevalence, which was calculated by the numerator of observed cases in cluster and denominator of population in the region multiply time interval (month) of the cluster.

Correlation and hierarchical clustering analysis

Chi-square tests were conducted to compare the proportion of PTB and SSP-TB purely spatial or spatial-temporal clusters between border and inland counties, then Pearson correlation coefficients were calculated by applying the PTB and SSP-TB prevalence of border counties or border counties within purely spatial clusters, Yunnan province, Myanmar, Laos, and Vietnam between 2005 to 2017. For the hierarchical clustering analysis, the first step was scaling and centering prevalence data; then define the similarity of prevalence by Euclidean distance; between-class distance and class agglomeration were defined by maximum between-class distance and complete linkage method; finally, the dendrogram presented the hierarchical clustering results.

The SaTScan™ software 9.6 (https://www.satscan.org) was applied for temporal, spatial, and spatial-temporal analysis. ArcGIS 10.2 (ESRI Inc) demonstrated geographical visualization for average TB prevalence and TB high-risk clusters. R software 3.5.2 (http://www.R-project.org) was applied for other statistical analysis. P < 0.05 indicates a statistical significance.

Results

Temporal distribution and clusters

From 2005 to 2018, a total of 381,855 PTB cases were notified in Yunnan. The average PTB prevalence was 59.1 per 100,000 population, meanwhile, the average prevalence of SSP-TB and SSN-TB was 20.3 and 27.1 per 100,000 population respectively. Figure 2 showed the X-12-ARIMA seasonal decomposition of raw notified PTB counts, original series of PTB cases (Fig. 2a) with trend cycle (Fig. 2b), segregated seasonal factors (Fig. 2c), and irregular components (Fig. 2d). PTB trend cycle could be stated with 5 stages, beginning with a peak in 2005, a sharp downward trend between 2006 to 2011, followed by a rapid climb in 2011 to 2012, a steady fluctuation between 2013 to 2016, and a recent upward trend in 2017 to 2018. There was an obvious seasonality of notified PTB cases in isolated seasonal distribution, in each year the trough appeared in December, then sharply reached the peak in January. Though the seasonal factors declined in recent years, it periodically appeared across the entire study time frame of 14 years (Fig. 2e).

Fig. 2
figure2

Seasonal time series decomposition of monthly tuberculosis cases in Yunnan, China, 2005–2018. Original series (a), Trend cycle (b), Seasonal factors (c), and irregular components (d). Irregular components: the residual that removed trend cycle and seasonal factor from original series

The annually temporal clusters included January except 2017 (Table 1), range of TB prevalence clusters concentrated in spring and/or summer. In the entire study time frame, temporal aggregation interval of PTB prevalence was in January 2005 to August 2006, for SSP-TB and SSN-TB the interval was from January 2008 to October 2011, and from January 2013 to February 2017 respectively (P < 0.01).

Table 1 Temporal clustering of pulmonary tuberculosis cases in Yunnan, 2005–2018

Spatial distribution and clusters

The average prevalence of PTB, SSP-TB, and SSN-TB showed explicitly geographical heterogeneity in 129 counties of Yunnan from 2005 to 2018 (Fig. 1c,d,e). Annually spatial clustering analysis showed the dynamic process of PTB aggregations in space, the time-dependent most likely clusters were mainly concentrated in three regions and were shown in Fig. 3. In the year of 2005, 2006 and 2008, the cluster was in Linxiang county of Lincang prefecture. Another hotspot located in the northeast angle of Yunnan, which Zhenxiong and Weixin counties of Zhaotong prefecture were the center of clusters in 2007, and 2009 to 2015. Recently, the most likely cluster was in southwestern borders of Pu’er and Xishuangbanna prefectures between 2016 to 2018.

Fig. 3
figure3

Spatial clustering of notified pulmonary tuberculosis cases in Yunnan, China, 2005–2018

The spatial clustering in the entire study period was similar to the annual aggregation, the most likely cluster of PTB and SSN-TB prevalence was Zhenxiong county between 2005 and 2018 (Fig. 3 and Fig. 4b), but the remarkable aggregation of SSP-TB was in the southwest frontier counties (Fig. 4a). Other secondary clusters scattered in the central, southeast and northeast of Yunnan.

Fig. 4
figure4

Spatial clustering of notified sputum smear-positive and smear-negative tuberculosis cases in Yunnan, China, 2005–2018. SSP-TB clusters (a) and SSN-TB clusters (b). SSP-TB Sputum smear-positive tuberculosis, SSN-TB Sputum smear-negative tuberculosis

Spatial-temporal distribution and clusters

The spatial-temporal distribution showed a consistent pattern with purely spatial scan for TB prevalence. Figure 5 and Table 2 showed the 17 spatial-temporal clusters for PTB, 9 clusters for SSP-TB, and 15 clusters for SSN-TB in 129 counties from 2005 to 2018. The most likely cluster for PTB was in northeast angle of Yunnan with coordinates of 27.90 N and 105.00 E, contained Zhenxiong and Weixin county in Zhaotong prefecture with a cluster circle radius of 43.1 Kilometers, which the high-risk period from February 2011 to March 2015 (LLR = 3657.6, P < 0.001). The risk was 2.6 times higher of developing active TB among residents within cluster compared with outside. A total of 10,963 PTB cases reported within the cluster and the average monthly prevalence was 13.4/100,000. Besides, there were the maximum reported cases (20080) and coverage counties (18 counties) in secondary cluster 1; the highest average monthly TB prevalence (14.2 per 100,000 population) in secondary cluster 2; specifically, both these two secondary clusters were in border region and were the latest clusters (Fig. 5a,d).

Fig. 5
figure5

Spatial-temporal clustering of notified tuberculosis cases and monthly average tuberculosis prevalence of clusters in Yunnan, China, 2005–2018. Spatial-temporal clustering of notified PTB (a), Spatial-temporal clustering of notified SSP-TB (b), Spatial-temporal clustering of notified SSN-TB (c); Monthly average TB prevalence of PTB clusters (d), Monthly average TB prevalence of SSP-TB clusters (e), Monthly average TB prevalence of SSN-TB clusters (f). TB tuberculosis, PTB Pulmonary tuberculosis, SSP-TB Sputum smear-positive tuberculosis, SSN-TB Sputum smear-negative tuberculosis

Table 2 Spatial-temporal clusters of pulmonary tuberculosis cases in Yunnan, 2005–2018

Moreover, the spatial-temporal most likely cluster of SSP-TB and SSN-TB centered in northeastern Zhenxiong and Weixin county, particularly, the time frame of SSP-TB clusters aggregated before the year of 2012, though most of the SSN-TB clusters were detected after 2012 (Fig. 5b,c,e,f). More detail of these clusters could be found in additional files (Additional files 1 and 2).

Correlation and hierarchical clustering analysis

The association between regions and scanned clusters were shown in Table 3. For purely spatial scanning, a higher proportion of border counties were defined within PTB and SSP-TB clusters than inland counties(P < 0.01), meanwhile, for spatial-temporal scanning, there was no significant difference for the proportion of counties within and outside clusters between border and inland counties. Pairwise correlation coefficients between 25 border counties or 17 border counties within purely spatial clusters and GMS were shown in Fig. 6.

Table 3 Association between geographical regions and pulmonary tuberculosis clusters in Yunnan
Fig. 6
figure6

Correlation of tuberculosis prevalence among Greater Mekong Subregion and borders of Yunnan. All border counties of Yunnan (a), and border counties within spatial clusters (b). GMS Greater Mekong Subregion, TB Tuberculosis, SSP-TB Sputum smear-positive tuberculosis, MTB Myanmar TB incidence, MTBSSP Myanmar SSP-TB incidence, LTB Laos TB incidence, LTBSSP Laos SSP-TB incidence, VTB Vietnam TB incidence, VTBSSP Vietnam SSP-TB incidence, BMCTB Myanmar-China border TB prevalence, BMCTBSSP Myanmar-China border SSP-TB prevalence, BLCTB Laos-China border TB prevalence, BLCTBSSP Laos-China border SSP-TB prevalence, BVCTB Vietnam-China border TB prevalence, BVCTBSSP Vietnam-China border SSP-TB prevalence, YNTB Yunnan Province TB prevalence, YNTBSSP Yunnan Province SSP-TB prevalence

Similarity patterns of TB prevalence defined by hierarchical clustering for 25 border counties or 17 border counties within purely spatial clusters with GMS were showed in Fig. 7. There were 6 possible patterns of TB prevalence among borders and GMS regardless of whether the borders were within purely spatial clusters. PTB incidence of Myanmar, Vietnam, Laos were independent categories; border counties PTB prevalence clustered with Myanmar SSP-TB prevalence; meanwhile, border counties SSP-TB prevalence were similar to Yunnan SSP-TB pattern; moreover, the pattern indicated that PTB prevalence of Yunnan highly correlated with Vietnam and Laos SSP-TB prevalence.

Fig. 7
figure7

Hierarchical clustering dendrogram of tuberculosis prevalence pattern among Greater Mekong Subregion and borders of Yunnan. All border counties of Yunnan (a), and border counties within spatial clusters (b). GMS Greater Mekong Subregion, TB Tuberculosis, SSP-TB Sputum smear-positive tuberculosis, MTB Myanmar TB incidence, MTBSSP Myanmar SSP-TB incidence, LTB Laos TB incidence, LTBSSP Laos SSP-TB incidence, VTB Vietnam TB incidence, VTBSSP Vietnam SSP-TB incidence, BMCTB Myanmar-China border TB prevalence, BMCTBSSP Myanmar-China border SSP-TB prevalence, BLCTB Laos-China border TB prevalence, BLCTBSSP Laos-China border SSP-TB prevalence, BVCTB Vietnam-China border TB prevalence, BVCTBSSP Vietnam-China border SSP-TB prevalence, YNTB Yunnan Province TB prevalence, YNTBSSP Yunnan Province SSP-TB prevalence

Discussion

This study reveals temporal, spatial and spatial-temporal TB prevalence distribution in southwestern China at county-level in Yunnan province from 2005 to 2018. In brief, we detected TB high-risk time interval and high-epidemic areas for scanning spatial-temporal characteristics, yet identified similar patterns of TB prevalence among GMS. This study also presented the dynamic perspective of spatial-temporal PTB, SSP-TB and SSN-TB epidemic in Yunnan province at county-level.

The annual notification of PTB was 59.6 per 100,000 population of Yunnan in 2018. Although significant efforts made an annual decline rate of 1.5% from 2005, the high TB burden with absolute number of 28,618 cases reported in 2018, made a big challenge to achieve the goal of End TB in 2035 without the breakthrough of vaccine or new drug [22].

Time series decomposed secular trend showed that after PTB notification peak at 2005, the prevalence decreased first and then increased in recent years. Between 2003 and 2005, the detection of SSP-TB cases by the public health system more than doubled, from 30% of new cases to 80% [23]. The reason for the notification peak in 2005 was that China launched the direct internet-based reporting system for infectious diseases in 2004, thus greatly increased the TB reporting cases in 2005. After 2003, plenty efforts such as intensive DOTS implementation coverage, increased government commitment and improved public-health funding, all these measures focused on TB control lead to an acceleration of prevalence decline in the following decade [4, 24]. Recent years, active case finding strategy implementation increased the number of TB suspects with symptoms, which were 184,618 in 2018, almost doubled from 104,960 in 2015 [15, 25, 26], result in the current PTB prevalence upward in Yunnan between 2016 and 2018.

Seasonality was observed in different counties for TB notification. Interestingly, the peak month and trough month of TB notification were constant in the North Hemisphere regardless of the locations’ longitude. The peak months were roughly the same in American (spring, March) [27], South Korea (spring and summer) [28], Indian (spring, March to May) [29], Singapore (spring and summer, March and July) [30], China (spring, April) [5], Wuhan city (spring, March) [6] and Xingjiang autonomous prefecture of China (spring, March) [7]. Our study was consistent with these researches in the North Hemisphere, seasonal factors were observed and the peak in January and the secondary peak in May. The hypothesis of TB seasonality was related to the lack of sunshine and the lower temperature in winter. Vitamin D deficiency due to shorter daylight hours in winter [31], the temperature was inversely and lagged associated with TB incidence [32], all of which caused seasonality disease for the peak in spring and summer. In China, the Spring Festival effect should also be considered. Which means TB notification significant reduced during Spring Festival holidays, consequently, seasonal factors sharply declined in habitual Spring Festival month of February. Meanwhile, the purely temporal scan revealed the temporal clusters were concentrated in spring and summer in Yunnan each year. In the whole study time frame, the cluster interval for SSN-TB was from 2008 to 2011, though for SSN-TB was more recently from 2013 to 2017, which suggested the ongoing TB control policy should focus on SSN-TB in Yunnan.

Kulldorff’s scan statistics method was developed to evaluate temporal and geospatial distribution, it was applied to detect communicable disease, vector-borne diseases and cancer geospatial aggregation [33,34,35,36], meanwhile, the sensitivity of spatial-temporal statistics prompted early detection of disease outbreak and emergency disease from surveillance system [37, 38]. This powerful method showed the strength of statistical robustness and interpretability of analyzed results. Scan statistics were widely applied in study topic related to TB [39,40,41,42,43], whereas, data aggregated into large scales of administrative regions may ignore the disease variation in small size of population, information lose lead to inaccurate and insensitive conclusion [44], these national-level researches could not preciously detect localized cluster on the resolution of province or prefecture [9, 10, 45]. Meanwhile, due to the stochastic scan statistics sensitive to parameters, the analytical results on high-resolution scan of county-level may not stable. Small changes on the algorithm parameters lead to different results, especially in small size of population [13]. The fitness of setting the parameters is crucial to the analysis as a whole.

The purely spatial scan showed that the PTB in Yunnan were not randomly distributed, and the dynamic prevalence of PTB revealed three mainly aggregated regions, the hotspot of Yunnan north-eastern angle in Zhaotong prefecture was high frequently in clusters and hold 8 of 14 years in study interval. Previous studies examined the PTB clusters for Zhaotong prefecture were in towns of Zhengxiong and Weixin county [14, 46].

Spatial-temporal cluster pattern was in line with pure spatial scanning. Unexpectedly, spatial-temporal scan detected two clusters were implemented active cases finding, one was secondary cluster 2 (Lanping county) and another was cluster 11 (Dongchuan county) in PTB clusters. The time frame of clustering matched with activities of active cases finding [15, 47]. This suggested that by considering the cluster time interval, higher sensitivity and closer to reality outcome for the spatial-temporal scan. Besides, the time frame for SSP-TB clusters concentrated before 2012, though most of the SSN-TB clusters were defined after 2012, which indicated the decline of SSP-TB and the progress and achievement for tuberculosis control in Yunnan.

Our study found that the most recent cluster of PTB and the SSP-TB spatial cluster for the whole interval was in southwestern borders neighbored with Myanmar, Laos, and Vietnam. Furthermore, the correlation of TB prevalence among borders and GMS were relatively high. Strikingly, hierarchical clustering indicated that there were 6 subclasses for TB epidemic pattern among GMS, thus the borders’ TB prevalence was similar to Myanmar TB epidemic pattern. Based on the consistency of traditional and molecular epidemiology evidence which confirmed the relatively lower prevalence of Beijing genotype in the border region of Pu’er, Xishuangbanna, as well as Vietnam and Myanmar [48,49,50], we speculated that the residents living in the border region moved across the national boundary for livelihood while the air-borne disease of TB was carried beyond frontier. Recent high TB prevalence and high-risk temporal-spatial clusters in the GMS region suggested that cross-boundary intervention and international control policy should be implemented in these clusters.

Our study has some limitations. Firstly, the surveillance data did not contain covariates of patients’ demographic information for sex, age, etc., yet we did not introduce ecology factors like geographical, meteorological and economic situation, all of which could be possible indicators of TB incidence and prevalence. Secondly, we do not take account for unreported cases when using notifications data, since there was a risk of underestimated prevalence regardless of missing or unreported cases. Thirdly, it is difficult to collect the genetic and lower-level detailed geospatial information for TB among GMS, although it will advance the understanding of TB transmission among GMS. Further study should address these points.

Conclusions

This study defined the long-term trend and seasonality of PTB prevalence for Yunnan province. Besides, we applied scan statistics to detect the temporal, spatial and spatial-temporal clusters of PTB prevalence between 2005 to 2018. The most likely cluster for PTB concentrated in the northeastern angle of Yunnan in distant past, the most recent spatial-temporal cluster of PTB and the spatial cluster of SSP-TB was in southwestern borders with GMS, meanwhile, the SSN-TB clusters were aggregated between 2012 to 2018. The similarity prevalence pattern of PTB among GMS suggested that priority of localized preventing implement should focus on cross-boundary intervention, yet the controlling strategy and resource allocation for the whole province should take account of the SSN-TB prevention.

Availability of data and materials

Data of the study was not publicly available, the datasets used and analysed during the current study are available from the corresponding author on reasonable request.

Abbreviations

ARIMA:

Autoregressive integrated moving average

CDC:

Center for Disease Control and Prevention

China NTP:

National Tuberculosis Control Program in China

CISDCP:

China Information System for Disease Control and Prevention

DOTS:

Directly observed treatment, short-course strategy

GMS:

Greater Mekong Subregion

LLR:

Log-likelihood ratio

RR:

Relative risk

PTB:

Pulmonary tuberculosis

SSN-TB:

Sputum smear-negative tuberculosis

SSP-TB:

Sputum smear-positive tuberculosis

TB:

Tuberculosis

WHO:

World Health Organization

References

  1. 1.

    Stanaway JD, Afshin A, Gakidou E, Lim SS, Abate D, Abate KH, et al. Global, regional, and national comparative risk assessment of 84 behavioural, environmental and occupational, and metabolic risks or clusters of risks for 195 countries and territories, 1990–2017: a systematic analysis for the global burden of disease study 2017. Lancet. 2018;392:1923–94.

    Article  Google Scholar 

  2. 2.

    Kyu HH, Maddison ER, Henry NJ, Mumford JE, Barber R, Shields C, et al. The global burden of tuberculosis: results from the global burden of disease study 2015. Lancet Infect Dis. 2018;18:261–84.

    Article  Google Scholar 

  3. 3.

    World Health Organization. Global tuberculosis report 2018. Geneva: World Health Organization; 2018.

    Google Scholar 

  4. 4.

    Wang L, Zhang H, Ruan Y, Chin DP, Xia Y, Cheng S, et al. Tuberculosis prevalence in China, 1990–2010; a longitudinal analysis of national survey data. Lancet. 2014;383:2057–64.

    PubMed  Article  PubMed Central  Google Scholar 

  5. 5.

    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:2005–12.

    Google Scholar 

  6. 6.

    Yang X, Duan Q, Wang J, Zhang Z, Jiang G. Seasonal variation of newly notified pulmonary tuberculosis cases from 2004 to 2013 in Wuhan. China. PLoS One. 2014;9.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  7. 7.

    Wubuli A, Li Y, Xue F, Yao X, Upur H, Wushouer Q. Seasonality of active tuberculosis notification from 2005 to 2014 in Xinjiang. China PLoS One. 2017;12:e0180226.

    PubMed  Article  CAS  PubMed Central  Google Scholar 

  8. 8.

    Chen J-O, Qiu Y-B, Rueda ZV, Hou J-L, Lu K-Y, Chen L-P, et al. Role of community-based active case finding in screening tuberculosis in Yunnan province of China. Infect Dis Poverty. 2019;8:92.

    PubMed  PubMed Central  Article  Google Scholar 

  9. 9.

    Zhao F, Cheng S, He G, Huang F, Zhang H, Xu B, et al. Space-time clustering characteristics of tuberculosis in China, 2005-2011. PLoS One. 2013;8:2005–11.

    Google Scholar 

  10. 10.

    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:1–10.

    Article  Google Scholar 

  11. 11.

    Ge E, Zhang X, Wang X, Wei X. Spatial and temporal analysis of tuberculosis in Zhejiang Province, China, 2009-2012. Infect Dis Poverty. 2016;5:1–10.

    Article  Google Scholar 

  12. 12.

    Li L, Xi Y, Ren F. Spatio-temporal distribution characteristics and trajectory similarity analysis of tuberculosis in Beijing, China. Int J Environ Res Public Health. 2016;13.

  13. 13.

    Rao H, Shi X, Zhang X. Using the kulldorff’s scan statistical analysis to detect spatio-temporal clusters of tuberculosis in Qinghai Province, China, 2009-2016. BMC Infect Dis. 2017;17:1–11.

    Article  Google Scholar 

  14. 14.

    Huang L, Li XX, Abe EM, Xu L, Ruan Y, Cao CL, et al. Spatial-temporal analysis of pulmonary tuberculosis in the northeast of the Yunnan province, People’s republic of China. Infect Dis Poverty. 2017;6:1–11.

    Article  Google Scholar 

  15. 15.

    Yunnan Centre For Disease Control and Prevention. Yunnan tuberculosis control report 2018. Kunming: Yunnan Centre For Disease Control and Prevention; 2018.

    Google Scholar 

  16. 16.

    Statistics Bureau of Yunnan. Yunnan statistical yearbooks. 2018. http://www.stats.yn.gov.cn/. Accessed 26 Jun 2019.

  17. 17.

    Jaganath D, Wobudeya E, Sekadde MP, Nsangi B, Haq H, Cattamanchi A. Seasonality of childhood tuberculosis cases in Kampala, Uganda, 2010-2015. PLoS One. 2019;14:e0214555.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  18. 18.

    Census Bureau US. X-12-ARIMA reference manual. Washington DC: U.S. Census Bureau; 2011.

    Google Scholar 

  19. 19.

    Kulldorff M. A spatial scan statistic. Commun Stat - Theory Methods. 1997;26:1481–96.

    Article  Google Scholar 

  20. 20.

    Kulldorff M, Athas WF, Feurer EJ, Miller BA, Key CR. Evaluating cluster alarms: a space-time scan statistic and brain cancer in Los Alamos. New Mexico Am J Public Health. 1998;88:1377–80.

    PubMed  Article  CAS  PubMed Central  Google Scholar 

  21. 21.

    Kulldorff M. SaTScan™ user guide for version 9.6. 2018. https://www.satscan.org/techdoc.html. .

  22. 22.

    World Health Organization. WHO End TB Strategy. Geneva: World Health Organization. 2015. http://www.who.int/tb/post2015_strategy/en/. Accessed 11 Jul 2018.

  23. 23.

    Wang L, Liu J, Chin DP. Progress in tuberculosis control and the evolving public-health system in China. Lancet. 2007;369:691–6.

    PubMed  Article  PubMed Central  Google Scholar 

  24. 24.

    Wang Y. Data of the 5th National Prevalence Survey of tuberculosis. Beijing: The Chinese Military Medical Publishing House; 2011.

    Google Scholar 

  25. 25.

    Yunnan Centre For Disease Control and Prevention. Yunnan tuberculosis control report 2016. Kunming: Yunnan Centre For Disease Control and Prevention; 2016.

    Google Scholar 

  26. 26.

    Yunnan Centre For Disease Control and Prevention. Yunnan tuberculosis control report 2017. Kunming: Yunnan Centre For Disease Control and Prevention; 2017.

    Google Scholar 

  27. 27.

    Willis MD, Winston CA, Heilig CM, Cain KP, Walter ND, Mac Kenzie WR. Seasonality of tuberculosis in the United States, 1993-2008. Clin Infect Dis. 2012;54:1553–60.

    PubMed  PubMed Central  Article  Google Scholar 

  28. 28.

    Kim EH, Bae J-M. Seasonality of tuberculosis in the Republic of Korea, 2006-2016. Epidemiol Health. 2018;40:e2018051.

    PubMed  PubMed Central  Article  Google Scholar 

  29. 29.

    Kumar V, Singh A, Adhikary M, Daral S, Khokhar A, Singh S. Seasonality of tuberculosis in Delhi. India: A Time Ser Anal Tuberc Res Treat. 2014;2014:1–5.

    Google Scholar 

  30. 30.

    Wah W, Das S, Earnest A, Lim LKY, Chee CBE, Cook AR, et al. Time series analysis of demographic and temporal trends of tuberculosis in Singapore. BMC Public Health. 2014;14:1121.

    PubMed  PubMed Central  Article  Google Scholar 

  31. 31.

    Koh GCKW, Hawthorne G, Turner AM, Kunst H, Dedicoat M. Tuberculosis incidence correlates with sunshine: an ecological 28-year time series study. PLoS One. 2013;8:1–5.

    Google Scholar 

  32. 32.

    Xiao Y, He L, Chen Y, Wang Q, Meng Q, Chang W, et al. The influence of meteorological factors on tuberculosis incidence in Southwest China from 2006 to 2015. Sci Rep. 2018;8:1–8.

    Article  CAS  Google Scholar 

  33. 33.

    Fischer E, Pahan D, Chowdhury S, Oskam L, Richardus J. The spatial distribution of leprosy in four villages in Bangladesh: an observational study. BMC Infect Dis. 2008;8:125.

    PubMed  PubMed Central  Article  Google Scholar 

  34. 34.

    Pinchoff J, Chipeta J, Banda GC, Miti S, Shields T, Curriero F, et al. Spatial clustering of measles cases during endemic (1998–2002) and epidemic (2010) periods in Lusaka. Zambia BMC Infect Dis. 2015;15:121.

    PubMed  Article  PubMed Central  Google Scholar 

  35. 35.

    González R, Augusto OJ, Munguambe K, Pierrat C, Pedro EN, Sacoor C, et al. HIV incidence and spatial clustering in a rural area of southern Mozambique. PLoS One. 2015;10:e0132053.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  36. 36.

    Lal A, Hales S. Heterogeneity in hotspots: spatio-temporal patterns in neglected parasitic diseases. Epidemiol Infect. 2015;143:631–9.

    PubMed  Article  CAS  PubMed Central  Google Scholar 

  37. 37.

    Yih WK, Cocoros NM, Crockett M, Klompas M, Kruskal BA, Kulldorff M, et al. Automated influenza-like illness reporting—an efficient adjunct to traditional sentinel surveillance. Public Health Rep. 2014;129:55–63.

    PubMed  PubMed Central  Article  Google Scholar 

  38. 38.

    Greene SK, Peterson ER, Kapell D, Fine AD, Kulldorff M. Daily reportable disease spatiotemporal cluster detection, New York City, New York, USA, 2014–2015. Emerg Infect Dis. 2016;22:1808–12.

    PubMed  PubMed Central  Article  Google Scholar 

  39. 39.

    Tadesse S, Enqueselassie F, Hagos S. Spatial and space-time clustering of tuberculosis in Gurage zone. Southern Ethiopia PLoS One. 2018;13:1–14.

    Google Scholar 

  40. 40.

    Roth D, Otterstatter M, Wong J, Cook V, Johnston J, Mak S. Identification of spatial and cohort clustering of tuberculosis using surveillance data from British Columbia, Canada, 1990–2013. Soc Sci Med. 2016;168:214–22.

    PubMed  Article  PubMed Central  Google Scholar 

  41. 41.

    Dangisso MH, Datiko DG, Lindtjørn B. Spatio-temporal analysis of smear-positive tuberculosis in the Sidama zone. Southern Ethiopia. PLoS One. 2015;10:1–19.

    Google Scholar 

  42. 42.

    Gurjav U, Burneebaatar B, Narmandakh E, Tumenbayar O, Ochirbat B, Hill-Cawthorne GA, et al. Spatiotemporal evidence for cross-border spread of MDR-TB along the trans-Siberian railway line. Int J Tuberc Lung Dis. 2015;19:1376–82.

    PubMed  Article  CAS  PubMed Central  Google Scholar 

  43. 43.

    Higgs BW, Mohtashem M, Grinsdale J, Kawamura LM. Early detection of tuberculosis outbreaks among the San Francisco homeless: trade-offs between spatial resolution and temporal scale. PLoS One. 2007;2.

    PubMed  PubMed Central  Article  Google Scholar 

  44. 44.

    Zhou X, Wang X, Yang K, Yang G, Lv S. Spatial epidemiology. Beijing: Science Press; 2009.

    Google Scholar 

  45. 45.

    Mao Q, Zeng C, Zheng D, Yang Y. Analysis on spatial-temporal distribution characteristics of smear positive pulmonary tuberculosis in China, 2004–2015. Int J Infect Dis. 2019;80:S36–44.

    Article  Google Scholar 

  46. 46.

    Huang L, Abe EM, Li XX, Bergquist R, Xu L, Xue JB, et al. Space-time clustering and associated risk factors of pulmonary tuberculosis in Southwest China. Infect Dis Poverty. 2018;7:1–11.

    Article  Google Scholar 

  47. 47.

    Zhang C, Ruan Y, Cheng J, Zhao F, Xia Y, Zhang H, et al. Comparing yield and relative costs of WHO TB screening algorithms in selected risk groups among people aged 65 years and over in China, 2013. PLoS One. 2017;12:e0176581.

    PubMed  PubMed Central  Article  CAS  Google Scholar 

  48. 48.

    Chen L, Pang Y, Ma L, Yang H, Ru H, Yang X, et al. First insight into the molecular epidemiology of mycobacterium tuberculosis isolates from the minority enclaves of southwestern China. Biomed Res Int. 2017;2017:1–9.

    Google Scholar 

  49. 49.

    Phyu S, Stavrum R, Lwin T, Svendsen OS, Ti T, Grewal HMS. Predominance of mycobacterium tuberculosis EAI and Beijing lineages in Yangon. Myanmar J Clin Microbiol. 2009;47:335–44.

    PubMed  Article  CAS  PubMed Central  Google Scholar 

  50. 50.

    Buu TN, Huyen MN, Lan NTN, Quy HT, Hen NV, Zignol M, et al. The Beijing genotype is associated with young age and multidrug-resistant tuberculosis in rural Vietnam. Int J Tuberc Lung Dis. 2009;13:900–6.

    PubMed  CAS  PubMed Central  Google Scholar 

Download references

Acknowledgments

This study was supported by the Yunnan Center for Disease Control and Prevention. We thank TB department colleagues for their assistance in data collection.

Funding

This work was funded by the Yunnan Provincial Health Bureau. Funding from the Yunnan Provincial Health Bureau was provided through the Yunnan Provincial Center for Tuberculosis Research (NO. 2014NS354). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Author information

Affiliations

Authors

Contributions

JC and LX conceived and designed the study. JC, YQ, RY, LL, JH, and KL provided the source data of the study. JC prepared software and performed the statistical analysis. JC and LX prepared the manuscript and interpreted the data. YQ, RY, LL, JH and KL assisted with the editing of the paper and provided critical comments. LL, JH, KL and LX revised it critically for important intellectual content. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Lin Xu.

Ethics declarations

Ethical approval and consent to participate

This retrospective study was consulted to the ethics committees of Yunnan Center for Disease Control and Prevention. Ethics approval was not available in this study because we did not include any data of patients’ personal information, including name, identity information, address, telephone number, etc. This study only reviewed of secondary aggregated data on county-level, therefore, waiver of written informed consent was taken and waived off ethical approval.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

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

Supplementary information

Additional file 1.

Spatial-temporal clusters of sputum smear-positive tuberculosis cases in Yunnan, 2005–2018.

Additional file 2.

Spatial-temporal clusters of sputum smear-negative tuberculosis cases in Yunnan, 2005–2018.

Rights and permissions

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

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Chen, J., Qiu, Y., Yang, R. et al. The characteristics of spatial-temporal distribution and cluster of tuberculosis in Yunnan Province, China, 2005–2018. BMC Public Health 19, 1715 (2019). https://doi.org/10.1186/s12889-019-7993-5

Download citation

Keywords

  • Tuberculosis
  • Spatial-temporal cluster
  • Scan statistics
  • Yunnan
  • Greater Mekong subregion