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

doi:10.21203/rs.2.16305/v2

Download PDF

Research article

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

https://doi.org/10.21203/rs.2.16305/v2

This work is licensed under a CC BY 4.0 License

Journal Publication

published 21 Dec, 2019

Read the published version in BMC Public Health →

You are reading this older preprint version

Read the latest preprint version →

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.

Infectious Diseases

Health Economics & Outcomes Research

Health Policy

Tuberculosis

Spatial-temporal cluster

Scan statistics

Yunnan

Greater Mekong Subregion

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

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 square kilometers while mountains occupy 84% of the territory, while the altitude various in tremendous range from 76 to 6740 meters. 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 (Figure 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 (Figure 1B).

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 (Figure 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.

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 (Figure 2A) with trend cycle (Figure 2B), segregated seasonal factors (Figure 2C), and irregular components (Figure 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 (Figure 2E).

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

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, 1D and 1E). Annually spatial clustering analysis showed the dynamic process of PTB aggregations in space, the time-dependent most likely clusters were mainly concentrated in 3 regions and were shown in Figure 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.

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 (Figure 3 and Figure 4B), but the remarkable aggregation of SSP-TB was in the southwest frontier counties (Figure 4A). Other secondary clusters scattered in the central, southeast and northeast of Yunnan.

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.10 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 (20 080) 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 (Figure 5A and 5D).

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 (Figure 5B, 5C, 5E, and 5F). More detail of these clusters could be found in additional files (Additional file 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 Figure 6.

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

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

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.

TB: Tuberculosis; PTB: pulmonary tuberculosis; WHO: World Health Organization; GMS: Greater Mekong Subregion; SSP-TB: sputum smear-positive tuberculosis; SSN-TB: sputum smear-negative tuberculosis; DOTS: directly observed treatment, short-course strategy; China NTP: National Tuberculosis Control Program in China; CDC: Center for Disease Control and Prevention; CISDCP: China Information System for Disease Control and Prevention; ARIMA: autoregressive integrated moving average; LLR: log-likelihood ratio; RR: relative risk

Ethical Approval and Consent to participate

Not applicable.

Consent for publication

Not applicable.

Availability of Data and Material

We do not wish to share our data. Data is available upon reasonable request to the corresponding author.

Competing Interests

The authors declare that they have no competing interests.

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.

Authors’ Contributions

JC and LX conceived and designed the study. JC, YQ, RY, LL, JH, KL and YY collected the data. JC performed the statistical analysis. JC and LX prepared the manuscript. All authors read and approved the final manuscript.

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.

Authors’ information

Jinou Chen¹, Yubing Qiu¹, Rui Yang¹, Ling Li¹, Jinglong Hou¹, Kunyun Lu¹, Lin Xu¹

¹Yunnan Center for Disease Control and Prevention, Kunming, China.

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 Stu. Lancet. 2018;392:1923–94.
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.
World Health Organization. Global Tuberculosis Report 2018. Geneva: World Health Organization; 2018.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Yunnan Centre For Disease Control and Prevention. Yunnan tuberculosis control report 2018. Kunming: Yunnan Centre For Disease Control and Prevention; 2018.
Statistics Bureau of Yunnan. Yunnan statistical yearbooks. 2018. http://www.stats.yn.gov.cn/. Accessed 26 Jun 2019.
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.
U.S. Census Bureau. X-12-ARIMA Reference Manual. Washington DC: U.S. Census Bureau; 2011.
Kulldorff M. A spatial scan statistic. Commun Stat - Theory Methods. 1997;26:1481–96.
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.
Kulldorff M. SaTScan^TM user guide for version 9.6. 2018. https://www.satscan.org/techdoc.html. Accessed 9 Jul 2019.
World Health Organization. WHO End TB Strategy. Geneva: World Health Organization. 2015. http://www.who.int/tb/post2015_strategy/en/. Accessed 11 Jul 2018.
Wang L, Liu J, Chin DP. Progress in tuberculosis control and the evolving public-health system in China. Lancet. 2007;369:691–6.
Wang Y. Data of the 5th National Prevalence Survey of Tuberculosis. Beijing: The Chinese Military Medical Publishing House; 2011.
Yunnan Centre For Disease Control and Prevention. Yunnan tuberculosis control report 2016. Kunming: Yunnan Centre For Disease Control and Prevention; 2016.
Yunnan Centre For Disease Control and Prevention. Yunnan tuberculosis control report 2017. Kunming: Yunnan Centre For Disease Control and Prevention; 2017.
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.
Kim EH, Bae J-M. Seasonality of tuberculosis in the Republic of Korea, 2006-2016. Epidemiol Health. 2018;40:e2018051.
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:1–5.
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.
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.
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.
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.
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.
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.
Lal A, Hales S. Heterogeneity in hotspots: spatio-temporal patterns in neglected parasitic diseases. Epidemiol Infect. 2015;143:631–9.
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.
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.
Tadesse S, Enqueselassie F, Hagos S. Spatial and space-time clustering of tuberculosis in Gurage Zone, Southern Ethiopia. PLoS One. 2018;13:1–14.
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.
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.
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.
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.
Zhou X, Wang X, Yang K, Yang G, Lv S. Spatial Epidemiology. Beijing: Science Press; 2009.
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.
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.
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.
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.
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.
Buu TN, Huyen MN, Lan NTN, Quy HT, Hen N V., 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.

Table 1. Temporal clustering of TB cases in Yunnan, 2005-2018

Categories	Year	Cluster time frame (start date to end date)	Observed cases (n)	Expected cases (n)	RR	LLR	P
PTB	2005	2005-1-1 to 2005-8-31	23978	17372	1.41	1181.73	<0.001
PTB	2006	2006-1-1 to 2006-6-30	16386	13015	1.27	418.54	<0.001
PTB	2007	2007-1-1 to 2007-6-30	17331	13114	1.34	639.52	<0.001
PTB	2008	2008-1-1 to 2008-5-31	14656	11087	1.33	538.42	<0.001
PTB	2009	2009-1-1 to 2009-4-30	11055	8808	1.26	271.75	<0.001
PTB	2010	2010-1-1 to 2010-1-31	2998	2288	1.31	100.96	<0.001
PTB	2011	2011-1-1 to 2011-1-31	3022	2301	1.32	103.31	<0.001
PTB	2012	2012-1-1 to 2012-5-31	13594	11372	1.20	210.77	<0.001
PTB	2013	2013-1-1 to 2013-2-28	5487	4439	1.24	116.33	<0.001
PTB	2014	2014-1-1 to 2014-1-31	2997	2346	1.28	83.53	<0.001
PTB	2015	2015-1-1 to 2015-1-31	2938	2359	1.25	66.12	<0.001
PTB	2016	2016-1-1 to 2016-1-31	2761	2374	1.16	30.21	<0.001
PTB	2017	2017-5-1 to 2017-6-30	5020	4710	1.07	10.14	<0.001
PTB	2018	2018-1-1 to 2018-1-31	3004	2403	1.25	70.12	<0.001
PTB	2005-2018	2005-1-1 to 2006-8-31	53814	43601	1.27	1268.18	<0.001
SSP-TB	2005-2018	2008-1-1 to 2011-10-31	45941	35513	1.45	1988.7	<0.001
SSN-TB	2005-2018	2013-1-1 to 2017-2-28	69538	53061	1.52	3494.3	<0.001

TB tuberculosis, PTB Pulmonary tuberculosis, SSP-TB Sputum smear-positive tuberculosis, SSN-TB Sputum smear-negative tuberculosis, RR Relative risk, LLR Log-likelihood ratios

Table 2. Spatial-temporal clusters of PTB cases in Yunnan, 2005-2018

Cluster type	Cluster period	Coordinates/Radius	N	Observed cases (n)	Expected cases (n)	RR	LLR	P
Most likely cluster	2011-02-01 to 2015-03-31	(27.90 N, 105.00 E) / 43.10 km	2	10963	4299	2.60	3657.56	<0.001
Secondary cluster 1	2015-01-01 to 2018-12-31	(21.99 N, 100.31 E) / 240.95 km	18	20080	11739	1.75	2533.15	<0.001
Secondary cluster 2	2015-01-01 to 2018-12-31	(26.03 N, 98.86 E) / 75.94 km	3	4097	1459	2.83	1600.85	<0.001
Secondary cluster 3	2005-01-01 to 2007-06-30	(25.15 N, 100.55 E) / 119.09 km	18	12399	8468	1.48	817.74	<0.001
Secondary cluster 4	2005-01-01 to 2006-07-31	(23.30 N, 103.51 E) / 77.79 km	10	5034	3325	1.52	382.49	<0.001
Secondary cluster 5	2005-03-01 to 2009-03-31	(23.66 N, 105.70 E) / 0 km	1	1953	992	1.97	362.91	<0.001
Secondary cluster 6	2007-01-01 to 2010-06-30	(27.25 N, 103.41 E) / 87.96 km	5	6685	4731	1.42	362.18	<0.001
Secondary cluster 7	2005-02-01 to 2008-08-31	(24.96 N, 104.32 E) / 0 km	1	2021	1160	1.75	261.64	<0.001
Secondary cluster 8	2005-01-01 to 2007-07-31	(25.16 N, 102.66 E) / 14.67 km	2	2072	1247	1.67	228.04	<0.001
Secondary cluster 9	2010-10-01 to 2013-11-30	(25.89 N, 102.58 E) / 0 km	1	1230	763	1.61	120.85	<0.001
Secondary cluster 10	2006-01-01 to 2006-01-31	(25.76 N, 102.20 E) / 34.34 km	2	130	24	5.32	111.79	<0.001
Secondary cluster 11	2009-08-01 to 2013-09-30	(26.15 N, 103.05 E) / 0 km	1	1140	718	1.59	104.95	<0.001
Secondary cluster 12	2008-01-01 to 2008-10-31	(25.03 N, 103.71 E) / 0 km	1	590	307	1.92	102.22	<0.001
Secondary cluster 13	2007-04-01 to 2007-07-31	(26.30 N, 104.13 E) / 0 km	1	467	273	1.71	56.96	<0.001
Secondary cluster 14	2005-01-01 to 2006-09-30	(24.77 N, 103.42 E) / 0 km	1	407	241	1.69	47.25	<0.001
Secondary cluster 15	2012-01-01 to 2012-05-31	(25.79 N, 103.86 E) / 62.11 km	4	709	514	1.38	33.06	<0.001
Secondary cluster 16	2005-04-01 to 2005-06-30	(24.60 N, 102.53 E) / 27.58 km	2	142	85	1.67	15.81	0.012

PTB Pulmonary tuberculosis, N number of counties in the cluster, RR Relative risk,

LLR Log-likelihood ratios

Table 3. Association between geographical regions and TB clusters in Yunnan.

Characteristic of clusters	Regions	Counties within clusters ^†		Counties outside clusters^‡		χ²	P
Characteristic of clusters	Regions	n	(%)	n	(%)	χ²	P
Spatial clustering PTB cases
	Border*	18	72.0	7	28.0	9.16	<0.01
	Inland	40	38.5	64	61.5
Spatial clustering SSP-TB cases
	Border	17	68.0	8	32.0	8.71	<0.01
	Inland	37	35.6	67	64.4
Spatial clustering SSN-TB cases
	Border	15	60.0	10	40.0	3.82	0.051
	Inland	40	38.5	64	61.5
Spatial-temporal clustering PTB cases
	Border	14	56.0	11	44.0	0.01	0.95
	Inland	59	56.7	45	43.3
Spatial-temporal clustering SSP-TB cases
	Border	16	64.0	9	36.0	0.30	0.56
	Inland	60	57.7	44	42.3
Spatial-temporal clustering SSN-TB cases
	Border	13	52.0	12	48.0	0.50	0.48
	Inland	46	44.2	58	55.8

PTB Pulmonary tuberculosis, SSP-TB Sputum smear-positive tuberculosis, SSN-TB Sputum smear-negative tuberculosis,

† Counties within the most likely or secondary clusters

‡ Counties outside the most likely or secondary clusters

* Border: 25 border counties in Yunnan.

Download PDF

Journal Publication

published 21 Dec, 2019

Read the published version in BMC Public Health →

Editorial decision: Minor revision
18 Nov, 2019
Editor assigned by journal
13 Nov, 2019
Submission checks completed at journal
12 Nov, 2019
Editor invited by journal
12 Nov, 2019

You are reading this older preprint version

Read the latest preprint version →

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

Status:

Journal Publication

Version 2

Abstract

Figures

Background

Methods

Results

Discussion

Conclusions

List of Abbreviations

Declarations

References

Tables

Supplementary Files

Status:

Journal Publication

Version 2