- Research Article
- Open Access
Spatial-temporal analysis of pulmonary tuberculosis in the northeast of the Yunnan province, People’s Republic of China
Infectious Diseases of Povertyvolume 6, Article number: 53 (2017)
The number of pulmonary tuberculosis (PTB) cases in China ranks third in the world. A continuous increase in cases has recently been recorded in Zhaotong prefecture-level city, which is located in the northeastern part of Yunnan province. This study explored the space-time dynamics of PTB cases in Zhaotong to provide useful information that will help guide policymakers to formulate effective regional prevention and control strategies.
The data on PTB cases were extracted from the nationwide tuberculosis online registration system. Time series and spatial cluster analyses were applied to detect PTB temporal trends and spatial patterns at the town level between 2011 and 2015 in Zhaotong. Three indicators of PTB treatment registration history were used: initial treatment registration rate, re-treatment registration rate, and total PTB registration rate.
Seasonal trends were detected with an apparent symptom onset peak during the winter season and a registration peak during the spring season. A most likely cluster and six secondary clusters were identified for the total PTB registration rate, one most likely cluster and five secondary clusters for the initial treatment registration rate, and one most likely cluster for the re-treatment registration rate. The most likely cluster of the three indicators had a similar spatial distribution and size in Zhenxiong County, which is characterised by a poor socio-economic level and the largest population in Yunnan.
This study identified temporal and spatial distribution of PTB in a high PTB burden area using existing health data. The results of the study provide useful information on the prevailing epidemiological situation of PTB in Zhaotong and could be used to develop strategies for more effective PTB control at the town level. The cluster that overlapped the three PTB indicators falls within the geographic areas where PTB control efforts should be prioritised.
Tuberculosis (TB) is an infectious disease caused by the bacteria Mycobacterium tuberculosis. Pulmonary tuberculosis (PTB) is a result of an attack on the lungs by Mycobacterium tuberculosis. Tuberculosis continues to be a major public health problem in China with an estimated one million incident cases reported which alone contributed 11% to the global TB incident in 2010 . Therefore, the Chinese national authorities established the 12th five-year National TB Control Programme in 2011 to effectively control TB. The implementation of the integrated control strategies has helped control the burden of PTB in China during the last 5 years. However, a continuous increase of the PTB burden has recently been observed in Yunnan province, especially in Zhaotong prefecture-level city.
The number of active PTB cases in Zhaotong is ranked highest in Yunnan and contributes almost one-quarter of the provincial registered cases yearly . The epidemic disparity of PTB between Zhaotong and other cities in Yunnan appears to follow an increasing trend yearly. The notification rate of PTB in Zhaotong is twice and four times higher than the rate reported for the province and for Yuxi city, which has the lowest notification rate in Yunnan, respectively .
The use of spatial-temporal cluster analysis has been applied for the detection of infectious disease hotspots in recent years [3,4,5]. It has also been successfully used for the identification of TB clusters, with meaningful results yielded [6,7,8]. Previous studies have shown that PTB is an airborne infectious disease with spatial and temporal heterogeneous distribution, and it is believed that a better understanding of the spatial epidemiology of PTB will help guide policymakers and different stakeholders to formulate effective regional prevention and control strategies [9,10,11].
Until now, few studies have been conducted to explore the spatial epidemiology of TB in Yunnan, China. Therefore, we conducted a spatial-temporal cluster analysis using town-level PTB registered data from 2011 to 2015 in Zhaotong to determine the clustering areas of the PTB epidemic and provide evidence to the local TB control programme for strategy development and intensified intervention.
Yunnan Province is located in southwest China and has a relatively poor socio-economic status. In 2015, the average gross domestic product per capita was 29 015 Yuan, and it was only 13 060 Yuan in Zhaotong City. There are 47 million inhabitants in Yunnan . The province covers an area of 394 100 km2, of which 94% is mountainous. There are 16 prefecture-level cities in Yunnan with a total of 129 counties. Zhaotong, located in the northeast of Yunnan, is strategically located in the province as it borders both the Guizhou and Sichuan provinces. There are 5.43 million inhabitants in Zhaotong, covering an area of 23 021 km2, containing 11 counties and 143 towns, with 23 ethnic minorities living in the mountainous and hilly areas.
Data collection and management
The data on TB cases were obtained from the database of the Yunnan provincial TB programme at the Yunnan Provincial Center for Disease Control and Prevention (Yunnan CDC). All the cases were reported to the nationwide TB online registration system. The epidemiological data of each case were collected, including: age, gender, current address, date of TB symptom onset, date of diagnosis, result of smear microscopy, therapeutic category and result of treatment. Town-level PTB registration data in Zhaotong from 2011 to 2015 were extracted from the database system. Three indicators of PTB treatment registration history were used: initial treatment registration rate, re-treatment registration rate and total PTB registration rate. All the cases were geocoded and matched to the town-level polygon maps from the geographic information system (Geographic database from China CDC) at a 1:1 000 000 scale as the layer’s attribute table by the same identified number.
Furthermore, the town-level point layer containing latitudes and longitudes of central points for each town was created using the ArcGIS v10.2 software (ESRI, Redlands, CA, USA). The geocoding process was applied, as described in previous studies [13,14,15]. Demographic information and the administrative code of Zhaotong were obtained from the Annual Statistical Report and integrated into the study database.
Time series analysis and descriptive analysis
All the PTB cases were aggregated by month in relation to their registration date and their symptom onset date, to identify the temporal patterns of the disease. The time series included 60 months in total from January 2011 to December 2015, and was examined using the EXCELL 2007 (see Fig. 1). The gender, age and the year of infection of the PTB cases were also aggregated to identify the demographic characteristics of the disease yearly. Time series and descriptive analyses were conducted at the town level as all the PTB cases registered in Zhaotong during the period between 2011 and 2015 were included in our analyses.
Spatial and space-time scan statistic
The spatial and space-time scan statistic was performed using SaTScanTM software v9.4.2 . The purely spatial scan statistic was applied to identify geographic areas with high PTB registration rates that statistically significantly exceeded nearby areas. The standard purely spatial scan statistic imposes a circular window on the map. The space-time scan statistic is defined by a cylindrical window with a circular geographic base and with height corresponding to time. The base is defined as the purely spatial scan statistic, while the height reflects the time period of potential clusters. The cylindrical window is then moved in space and time, so that for each possible geographical location and size, it also visits each possible time period. For this analysis, a Poisson based model was used, where the number of events in an area is Poisson distributed according to a known underlying population at risk. PTB cases registered in each town were used and recorded with the population in the same town which was assumed as the population in the Poisson probability model. Under the null hypothesis, and when there are no covariates, the expected number of cases in each area is proportional to its population size, or to the person-years in that area. For this analysis, the geographic size of the window was limited to half the expected number of cases and the time period was also limited to half the total time period. The test of significance of the identified clusters was based on comparing the likelihood ratio test statistics against a null distribution obtained from a Monte Carlo simulation. The number of permutation was set to 999 and the significance level was set as 0.05. The space-time scan window with the maximum likelihood value was the most likely cluster and other significant windows were called secondary clusters .
Descriptive analysis of PTB cases
A total of 23 605 confirmed active PTB cases were included in this study, of which 21 lived in other provinces. Table 1 shows that the number of male cases was twice that of female cases in any given year. In addition, significant shares of PTB infections were found in the age groups of 15–30 years (over 30%) and 30–45 years (around 25%). Among the registered cases, the percentage of sputum smear-positive (SS+) PTB cases declined yearly (see Table 1).
Temporal patterns of PTB cases
Figure 1 shows the monthly time series of the numbers of PTB cases registered from 2011 to 2015. There was an evident trend variation; with the maximum number of PTB cases registered in the spring each year and a volatile declining trend the following summer and autumn, and then a minimum number of cases were recorded in winter. Peaks of the PTB symptom onset date were observed every January from 2011 to 2015, and then declined gradually after March. A non-linear correlation was observed both in the registration date (χ 2 = 336.831, P < 0.001) and symptom onset date (χ 2 = 6.703, P = 0.010), when comparing the different seasons amalgamated across all years (see Fig. 1).
Spatial patterns of PTB cases
The spatial variations of total PTB registration rates between 2011 and 2015 in Zhaotong shows that the annualised average rate at the town level ranged from 8.09 to 175.50 per 100 000 people. The highest registration rates were found in the towns of Zhenxiong County and border towns in Weixin County, such as Shuitian, Linfeng and Shuanghe (Fig. 2). In addition, Yongfeng and Xiaolongdong towns in Zhaoyang district; Yanjing, Dousha and Shizi towns in Yanjin County; and Cuihua town in Daguan County had the highest registration rates. The rates were relatively low in Yongshan County, Qiaojia County and Yiliang County. The cut-off values were the same, which is a visual display to show the registration rate of each town changes over time. The registration rates of PTB cases were significantly different during the years under study, and the registration rate in Wufeng town in Zhenxiong County appeared to increase from 148.94 in 2011 to 191.03 per 100 000 people in 2015. The disease demonstrated high heterogeneity in space and time.
The spatial variations of initial treatment registration rates between 2011 and 2015 in Zhaotong shows that the annual cases per 100 000 people were 76.8. The highest registration rates were found in Wufeng, Yudong, Wanchang and Chishuiyuan towns in Zhenxiong County and Yongfeng town in Zhaoyang County (Fig. 3). The registration rates of initial treatment cases were significantly different between 2011 and 2015, while the rates in Wufeng town and Yudong town in Zhenxiong County increased continuously between 2011 and 2015.
Figure 4 shows the spatial variations of PTB re-treatment registration rates between 2011 and 2015 in Zhaotong. Annual cases per 100 000 people were 4.3. The highest registration rates were found in Wanchang and Yanyuan towns of Zhenxiong County, Yanjing town in Yanjin County and Shanggaoqiao town in Daguan County. The registration rates of re-treatment cases were significantly different in the studied years. There was an obvious decrease in the registration of re-treatment rates in Wanchang and Yanyuan. However, the rate in Shanggaoqiao continuously increased from 2011 to 2015.
And the Fig. 5 shows the space-time scan analysis results regarding the total PTB cases using the Poisson model, including one most likely cluster and six secondary clusters. Two secondary clusters were statistically insignificant with P > 0.05. The clusters are distinguished in different colours on Fig. 5. The main cluster occupies large areas of the Zhenxiong County, covering 29 towns and 5 994 PTB cases. The cluster is significant, with long-term persistence observed from early spring 2013 to the end of summer 2015. The important finding is that people within this cluster were facing more than twice risk of acquiring PTB infection, compared to those outside the hotspot (relative risk, RR = 2.01). In addition, about three-quarters of the secondary clusters occurred in the spring and persisted for about 6 months. The other quarter of secondary clusters included Zhongming and LongAn towns in Yiliang County (secondary cluster 4) located in the south of Zhaotong, which emerged in December 2015 but lasted for a short period of only 1 month. The relative RR of these secondary clusters varied from 3.91 to 1.70.
The space-time cluster for initial treatment and re-treatment cases was detected using a Poisson-based model. One most likely cluster and five secondary clusters for a high registration rate of initial treatment, and one most likely cluster for a high registration rate of re-treatment were detected, with the results depicted in Fig. 6. The most likely cluster of the PTB initial treatment registration rate had the same clustering centre and size as the total PTB registration rate in the Wufeng-central regions of Zhenxiong County (radius 45.51 km, covering 27 towns in Zhenxiong and two towns in Weixin County), between 2014 and 2015 (RR = 1.86, P < 0.001), with 4 383 observed cases and 2 586 expected cases. In addition, the significant secondary clusters of the initial treatment registration rate had the same clustering centres as the total PTB registration rate in Fenghuang, Dousha and Cuihua towns. The remaining two secondary clusters were statistically insignificant with P > 0.05. The most likely cluster for re-treatment cases was found to exist at Zhongtun town in Zhenxiong County (radius 50.58 km, covering 27 towns in Zhenxiong and Miaogou town in Weixin), between 2011 and 2012 (RR = 5.78, P < 0.001), with 134 observed cases and four expected cases.
This study analysed the spatial-temporal status of PTB in Zhaotong city, northeast Yunnan, using the smallest administrative unit in the city-wide scale. Hence, it is the first study to analyse PTB risk at the town level in a highly epidemic area of Yunnan. The identification of TB risk areas using surveillance data based on the geographic unit is advantageous for policymakers . Systemic utilisation of cluster detection techniques for regular surveillance of TB is useful to effectively carry out TB control programmes . Therefore, it is important to equip policymakers with information on the spatial-temporal status of PTB to help guide the provision, planning and optimisation of PTB control strategies in Zhaotong.
In this study, a total of 23 605 PTB cases were registered between 2011 and 2015 in Zhaotong. Findings show that there may be seasonal trends, with an apparent symptom onset peak in the winter and a registration peak during the spring of the following year. This reflects the delay between the time that PTB symptoms become evident in patients and the period that patients seek medical attention. It is also important to stress that most TB patients experience delay in referral due to their lack of awareness that there are designated PTB clinics and seek treatment in general hospitals. This might also lead to late PTB registration in our TB online registration system . The possible reason of the symptom onset peak during the winter is a result of the reduction of exposure to ultraviolet rays from the sun and more indoor activities, which increases the possible susceptibility of people to infections with TB bacteria in the winter . Similar trends in seasonal patterns had been reported elsewhere . Resources should be distributed and planned for rationally in order to cope with the arrival of the incidence peak period and increase awareness-raising and running of focus groups during this period. It is of great significance to control the PTB epidemic situation and protect the elderly, patients with diabetes and chronic diseases and other susceptible groups.
The results obtained in this study show that the trend of the smear positive PTB cases percentage is very dramatic. The total number of smear positive cases decreased during the years under study, but there was a substantial increase in the number of smear-negative cases. The decrease in the number of smear-positive patients was due to the efforts for early case finding during the 12th five-year national TB integrated control plan that was initiated in 2011. Increased case finding efforts can result in a lower proportion of cases with bacteriological confirmation (that is smear positive), as some cases are found before the bacterial load increases to a detectable level . Estimates from the World Health Organization Western Pacific Region lead to a similar conclusion: The proportion of cases bacteriologically confirmed among all TB cases decreased from 52% in 2005 to 37% in 2013 in this region .
This study identified total PTB space-time clusters, which differ in their socio-economic status, demography structures and natural environment. The most likely cluster included 27 towns in Zhenxiong County and two neighbouring towns in Weixin County, which have the highest risk of PTB infection (RR = 2.01). This region in Zhaotong is vast and densely populated, in which mostly poor people with very low incomes or without serious means of livelihood live. Therefore, a high prevalence of PTB infections in these areas is associated with a low socio-economic status, as poverty is reported to be one of the factors responsible for a high PTB prevalence. The prevalence of PTB infection in Zhaoyang district, which is ranked second to Zhenxiong County in terms of PTB burden in Zhaotong, as shown in the cluster, might be associated with the migration of people into the district. Zhaoyang is the capital of Zhaotong and also serves as the commercial centre of the region. The movement of PTB patients from the different counties could facilitate PTB transmission in the district.
The PTB spatial-temporal characteristics were identified using different indicators. As an indicator of PTB burden, the total PTB registration rate lumps together the initial treatment registration rate and the re-treatment registration rate. The initial treatment registration rate is one indicator of severe epidemic situation and the retreatment registration rate is one indicator of poor supervision and follow- up. The statistically significant clusters of initial treatment cases were in line with the total PTB cases. It is obvious that a high TB epidemic situation might trigger the development of initial treatment cases That is, people within the situation are more likely infected with TB bacteria than those of outside the situation. Thus, to achieve the target of reducing PTB cases, adequate attention should be given to the infected population in these clusters, and it is also important to optimise the use of available resources for effective PTB control. It is advised that the control of PTB should be incorporated into the public health services with standardised service contents and methods. These include: encouraging the role of health workers to disseminate information to patients in order to shorten the time interval between incidence and registration; and carrying out patients’ active screening in this cluster for early detection (case finding) – this should be incorporated into the control strategy plan for the next five years. In addition, more detailed survey-based investigations are needed within and outside the identified clusters to evaluate the factors influencing PTB distribution and purposeful intervention measures to control PTB prevalence in this region should be put forward.
The most likely cluster of the three indicators chosen for this study has similar spatial distribution and size in Zhenxiong County. This shows that the high registration rate of total PTB cases is not only limited to the high registration rate of initial treatment, but also the high registration rate of re-treatment. This area carries a disproportional burden of PTB cases. The overlap of clusters of the three indicators shows that there is a high active PTB incidence in this region, but patients do not receive adequate attention, desired care or satisfactory management. Therefore, there is a large number of re-treatment cases. One of the major challenges that might arise in this area is that if the management of patients cannot keep up with the numbers of PTB patients found, this might strongly influence the clustering of PTB re-treatment cases within a few years. Therefore, it is important to enhance treatment and management quality of active PTB cases under effective control strategy programmatic conditions in Zhenxiong County.
Though our study demonstrated the usefulness of spatial and temporal clustering analysis, it also had a few limitations. First, the estimated risk of PTB infection might be underestimated in some areas because cases may not be recorded in the official registration system. Second, the power of cylindrical scan statistics could be limited in some irregular geography. Third, we cannot judge the recent transmission or reactivation of PTB due to a lack of laboratory evidence, but we hope to improve this in future studies by accessing more funding to carry out molecular-based studies on TB transmission. Fourth, the potential risk factors that could be associated with clustering were not assessed in this study. The high PTB prevalence is frequently associated with individual and socio-economic factors [25,26,27,28,29]. Fifth, we did not get each patient’s geographic coordinates in this study. However, we can consider using individual coordinates in future studies when we analyse on village scale.
Tuberculosis remains a public health burden in Zhaotong and our findings show that there are significant spatial-temporal characteristics of PTB at town scale in the region. Therefore, the findings of this study have provided useful information on the prevailing epidemiological situation of PTB in Zhaotong using existing health data and could be used to develop strategies for more effective PTB control at the town level.
Center of disease control and prevention
Wang L. Status analysis of tuberculosis prevention and control, China. Chin J Public Health. 2012;28:413–4.
Yunnan CDC. Yunnan tuberculosis control report 2015. Kunming: Yunnan Centre For Disease Control and Prevention; 2015.
Metras R, Porphyre T, Pfeiffer DU, Kemp A, Thompson PN, Collins LM, White RG. Exploratory space-time analyses of rift valley fever in South Africa in 2008–2011. PLoSNegl Trop Dis. 2012;6(8):e1808.
Chong S, Nelson M, Byun R, Harris L, Eastwood J, Jalaludin B. Geospatial analyses to identify clusters of adverse antenatal factors for targeted interventions. Int J Health Geogr. 2013;12:46.
Chimoyi LA, Musenge E. Spatial analysis of factors associated with HIV infection among young people in Uganda, 2011. Bmc Public Health. 2014;14:555.
Baker BJ, Moonan PK. Characterizing tuberculosis genotype clusters along the United States–Mexico border [Short communication]. Int J Tuberc Lung Dis. 2014;18(3):289–91.
Hassarangsee S, Tripathi N, Souris M. Spatial pattern detection of tuberculosis: a case study of Si Sa Ket Province, Thailand. Int j Env res Pub he. 2015;12(12):16005–18.
Dangisso MH, Datiko DG, Lindtjørn B. Spatio-temporal analysis of smear-positive tuberculosis in the sidama zone,Southern Ethiopia. Plos One. 2015;10(6):e126369.
Liu Y, Li X, Wang W, Li Z, Hou M, He Y, Wu W, Wang H, Liang H, Guo X. 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.
Maciel EL, Pan W, Dietze R, Peres RL, Vinhas SA, Ribeiro FK, Palaci M, Rodrigues RR, Zandonade E, Golub JE. Spatial patterns of pulmonary tuberculosis incidence and their relationship to socio-economic status in Vitoria. Brazil Int J Tuberc Lung Dis. 2010;14(11):1395–402.
Musenge E, Vounatsou P, Collinson M, Tollman S, Kahn K. The contribution of spatial analysis to understanding HIV/TB mortality in children: a structural equation modelling approach. Glob Health Action. 2013;6:38–47. http://dx.doi.org/10.3402/gha.v6i0.19266.
Statistics Bureau of Yunnan. Yunnan national economic and social development statistical bullet. 2015. http://www.sei.gov.cn/ShowArticle.asp?ArticleID = 263398. Accessed Mar 2016.
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:11.
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.
Zhao F, Cheng S, He G, Huang F, Zhang H, Xu B, Murimwa TC, Cheng J, Hu D, Wang L. Space-time clustering characteristics of tuberculosis in China, 2005–2011. PLoS One. 2013;8(12):e83605.
Kulldorff M. SaTScanTM User Guide for version 9.4.4. http://www.satscan.org/. Accessed 29 Jan 2016.
Deng T, Huang Y, Yu S, Gu J, Huang C, Xiao G, Hao Y. Spatial-temporal clusters and risk factors of hand, foot, and mouth disease at the district level in Guangdong Province, China. PLoS One. 2013;8(2):e56943.
Dietz NA, Sherman R, Mackinnon J, Fleming L, Arheart KL, Wohler B, Lee DJ. Toward the identification of communities with increased tobacco-associated cancer burden: application of spatial modeling techniques. J Carcinog. 2011;10:22.
Touray K, Adetifa IM, Jallow A, Rigby J, Jeffries D, Cheung YB, Donkor S, Adegbola RA, Hill PC. Spatial analysis of tuberculosis in an urban west African setting: is there evidence of clustering? Trop Med Int Health. 2010;15(6):664–72.
Wei X, Zou G, Yin J, Walley J, Sun Q. Comparing patient care seeking pathways in three models of hospital and TB programme collaboration in China. BMC Infect Dis. 2013;13:93.
Thorpe LE, Frieden TR, Laserson KF, Wells C, Khatri GR. Seasonality of tuberculosis in India: is it real and what does it tell us? Lancet. 2004;364(9445):1613–4.
Li XX, Wang LX, Zhang H, Du X, Jiang SW, Shen T, Zhang YP, Zeng G. Seasonal variations in notification of active tuberculosis cases in China, 2005–2012. PLoS One. 2013;8(7):e68102.
Systematic screening for active tuberculosis: principles and recommendations. Geneva, World Health Organization, 2013 (http://www.who.int/tb/tbscreening/en/. Accessed 29 March 2016.
Tom Hiatt, Nobuyuki Nishikiori: Epidemiology and control of tuberculosis in the Western Pacific Region: update with 2013 case notification data. World Health Organization Regional Office for the Western Pacific.
Harling G, Castro MC. A spatial analysis of social and economic determinants of tuberculosis in Brazil. Health Place. Manila. 2014;25:56–67. doi:10.5365/wpsar.2015.6.010.
Pereira AGL, Medronho RDA, Escosteguy CC, Valencia LIO, Magalhães MDAF. Spatial distribution and socioeconomic context of tuberculosis in Rio de Janeiro, Brazil. Revista de SaúdePública. 2015;49:1–8.
Randremanana RV, Sabatier P, Rakotomanana F, Randriamanantena A, Richard V. Spatial clustering of pulmonary tuberculosis and impact of the care factors in Antananarivo City. Trop Med Int Health. 2009;14(4):429–37.
Li XX, Wang LX, Zhang H, Jiang SW, Fang Q, Chen JX, Zhou XN. Spatial variations of pulmonary tuberculosis prevalence co-impacted by socio-economic and geographic factors in People’s Republic of China, 2010. Bmc Public Health. 2014;14:257.
Rakotosamimanana S, Mandrosovololona V, Rakotonirina J, Ramamonjisoa J, Ranjalahy JR, Randremanana RV, Rakotomanana F. Spatial analysis of pulmonary tuberculosis in Antananarivo Madagascar: tuberculosis-related knowledge, attitude and practice. PLoS One. 2014;9(11):e110471.
The authors would like to thank colleagues of the TB programme in the Yunnan CDC for their assistance in the data collection.
This study was supported by the National Special Science and Technology Project for Major Infectious Diseases of China (Grant No. 2012ZX10004-220, 2016ZX10004222-006) and the China-UK Global Health Support Programme (Grant No.: GHSP-CS-OP1-01). The Forth Round of Three-Year Public Health Action Plan of Shanghai, China(No. 15GWZK0101, GWIV-29). The funders had no role in the study design, data collection and analysis, decision to publish, or preparation of the paper.
Availability of data and materials
In order to protect the privacy of TB patients, we will not share the original copies of information from the database. We would like to share statistical results of this study. If anyone needs these data, please contact the corresponding author for a soft copy.
LH and S-ZL designed the research and tools, and supervised the study. LH analysed the data and wrote the paper. X-XL and S-ZL revised the paper and provided critical comments. EMA assisted in the editing of the paper. LX, YR and C-LC revised and approved the final paper and accept accountability for all aspects of the work. All authors approve the paper for publication.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
In this study, TB data were collected by routine TB surveillance and control activities. Therefore ethical consent was not essential.