This paper is in the following e-collection/theme issue:

Cross-Sectional Studies in Public Health (289) Descriptive Epidemiology and Population Size Estimates (109) Infectious Diseases (non-STD/STI) (1592) Innovative Methods in Public Health and Surveillance (336)

Published on in Vol 9 (2023)

Preprints (earlier versions) of this paper are available at https://preprints.jmir.org/preprint/49859, first published .
Trend of the Tuberculous Pleurisy Notification Rate in Eastern China During 2017-2021: Spatiotemporal Analysis

Trend of the Tuberculous Pleurisy Notification Rate in Eastern China During 2017-2021: Spatiotemporal Analysis

Trend of the Tuberculous Pleurisy Notification Rate in Eastern China During 2017-2021: Spatiotemporal Analysis

Authors of this article:

Ying Zhou1 Author Orcid Image ;   Dan Luo2 Author Orcid Image ;   Kui Liu3, 4 Author Orcid Image ;   Bin Chen3 Author Orcid Image ;   Songhua Chen3 Author Orcid Image ;   Junhang Pan3 Author Orcid Image ;   Zhengwei Liu3 Author Orcid Image ;   Jianmin Jiang3, 5 Author Orcid Image
Article Authors Cited by (4) Tweetations (3) Metrics

    Original Paper

    1School of Public Health, Hangzhou Normal University, Hangzhou, China

    2School of Public Health, Hangzhou Medical College, Hangzhou, China

    3Department of Tuberculosis Control and Prevention, Zhejiang Provincial Center for Disease Control and Prevention, Hangzhou, China

    4National Centre for Tuberculosis Control and Prevention, Chinese Center for Disease Control and Prevention, Beijing, China

    5Key Laboratory of Vaccine, Prevention and Control of Infectious Disease of Zhejiang Province, Zhejiang Provincial Center for Disease Control and Prevention, Hangzhou, China

    *these authors contributed equally

    Corresponding Author:

    Jianmin Jiang, PhD

    Department of Tuberculosis Control and Prevention

    Zhejiang Provincial Center for Disease Control and Prevention

    3399 Binsheng Road

    BinJiang District

    Hangzhou, 310000

    China

    Phone: 86 057187115009

    Fax:86 057187115100

    Email: jmjiang@cdc.zj.cn


    Background: Tuberculous pleurisy (TP) presents a serious allergic reaction in the pleura caused by Mycobacterium tuberculosis; however, few studies have described its spatial epidemiological characteristics in eastern China.

    Objective: This study aimed to determine the epidemiological distribution of TP and predict its further development in Zhejiang Province.

    Methods: Data on all notified cases of TP in Zhejiang Province, China, from 2017 to 2021 were collected from the existing tuberculosis information management system. Analyses, including spatial autocorrelation and spatial-temporal scan analysis, were performed to identify hot spots and clusters, respectively. The prediction of TP prevalence was performed using the seasonal autoregressive integrated moving average (SARIMA), Holt-Winters exponential smoothing, and Prophet models using R (The R Foundation) and Python (Python Software Foundation).

    Results: The average notification rate of TP in Zhejiang Province was 7.06 cases per 100,000 population, peaking in the summer. The male-to-female ratio was 2.18:1. In terms of geographical distribution, clusters of cases were observed in the western part of Zhejiang Province, including parts of Hangzhou, Quzhou, Jinhua, Lishui, Wenzhou, and Taizhou city. Spatial-temporal analysis identified 1 most likely cluster and 4 secondary clusters. The Holt-Winters model outperformed the SARIMA and Prophet models in predicting the trend in TP prevalence.

    Conclusions: The western region of Zhejiang Province had the highest risk of TP. Comprehensive interventions, such as chest x-ray screening and symptom screening, should be reinforced to improve early identification. Additionally, a more systematic assessment of the prevalence trend of TP should include more predictors.

    JMIR Public Health Surveill 2023;9:e49859

    doi:10.2196/49859

    Keywords

    tuberculous pleurisy; spatio-temporal; epidemiology; prediction; time series 


    Humanity has been affected by tuberculosis (TB) for centuries, which is an infectious agent–caused disease considered the leading cause of mortality worldwide since the beginning of recorded history. Among people who contract TB each year, the majority are adults, with the incidence of TB being higher in men than in women. According to the Global Tuberculosis Report of 2022 released by the World Health Organization, the estimated TB incidence rate in 2021 increased by 3.6% from 2020, and China accounted for 7.4% of TB cases worldwide. Furthermore, the number of TB deaths worldwide successively increased from 1.4 million in 2019 to 1.6 million in 2021, only after the onset of the COVID-19 pandemic [1]. Mycobacterium tuberculosis generally affects not only the lungs but also other sites such as the pleura [2].

    Tuberculosis pleurisy (TP) is a common clinical manifestation of TB post primary infection or due to the reactivation of latent TB [3]. Available documentation shows that TP accounts for about 4% of total tuberculosis cases in the United States and Brazil, 20% of cases in South Africa, and 7.3% of cases in Korea [4]. In the United States, the incidence of both TP and TB has been declining, but to a lesser degree for TP [5]. In 2018, it was categorized as pulmonary tuberculosis (PTB) in China [6]. TP can be divided into dry and exudative pleurisy. Dry pleurisy is an early inflammatory reaction of the pleura, usually without obvious imaging abnormalities. Exudative pleurisy is characterized by unilateral or bilateral pleural effusions, which often occur as thickened pleural adhesions with slow absorption [7]. The recommended standard regimen for TP is up to 8 months, which is longer than the usual treatment period for TB [8]. Previous studies have shown that the prevalence of TP is low in women, while the number of patients with TP is higher in eastern China, among farmers and patients with AIDS [9-11]. However, the temporal periodicity and spatial clustering characteristics of TP have not been explored in eastern China.

    Zhejiang Province is located in the eastern region of China and has a developed economy. In recent years, certain parts of the province have experienced a high burden of TB. Understanding the spatial and temporal patterns of TB incidence can provide insights into the clusters of this disease [12]. Spatial-temporal analysis has been widely used to describe the epidemiological distribution and aggregation of infectious diseases in both space and time; however, few studies have explored the spatial epidemiological characteristics of TP. In comparison to traditional epidemiological techniques, tools using geographic information systems have significant advantages in tracking disease, identifying areas of clustering, and analyzing the spread of disease across communities [13]. Spatiotemporal cluster detection has played a crucial role in understanding infectious disease spread, with its dynamic course and dimensions of time and place, and a combination of both [14]. In addition, time-series models were selected to provide scientific clues for predicting the prevalence of diseases such as desiccation syndrome and AIDS [15,16]. However, there are few such studies on TP.

    Therefore, this study aimed to use spatial-temporal analysis to identify the epidemiological characteristics of TP. Furthermore, the spatial differentiation of TP can provide a visual map depicting the diversity of TP occurrence across different regions. Additionally, this study used time-series models to provide scientific clues for the control of and interventions targeting TP infections in the province.


    Overview of the Study Area

    Zhejiang Province is located on the eastern coast of China, with a land area of 101,800 km2, accounting for 1.06% of the country’s land area [17]. The province encompasses plains, mountains, coasts, islands, and lakes. Zhejiang Province comprises 90 counties and cities, including 2 subprovincial cities and 9 prefecture-level cities [17,18].

    Data Source

    All recorded cases of TP in Zhejiang Province from 2017 to 2021 were obtained from the TB Information Management System (TBIMS), which is a first-generation web-based information system designed and built by the National Tuberculosis Control and Prevention Center in 2005 [19]. The system collected data on demographics, diagnosis, laboratory results, and treatment outcomes of all patients with TB during the specified period. In addition, demographic data for Zhejiang Province were obtained from the China Information System for Disease Control and Prevention.

    Definition

    Diagnostic criteria for TP were as follows: (1) confirmed patients: cultures of pleural biopsy specimens were positive, and the strain was identified as M tuberculosis, changes in tuberculosis confirmed through pleural tissue biopsy; and (2) clinically diagnosed patients: imaging examination showing pleural effusion with (i) a moderately positive or strongly positive outcome on the tuberculin pure protein derivative skin test, (ii) a positive outcome on the γ-interferon release test, or (iii) a positive outcome on the M tuberculosis antibody test [20]. Delays in seeking medical treatment were indicated when the interval between the appearance of the first symptom and the visit to the designated hospital was > 2 weeks, and the delay in diagnosis was the time between the first consultation and a confirmed diagnosis [21].

    General Epidemiological Characteristics of Patients With TP

    The features of patients with TP in Zhejiang Province from 2017 to 2021 are described in the context of sex, age group, city, antituberculosis treatment, mode of case finding, treatment history, delay in access to care and final diagnosis, and treatment outcomes.

    Spatial Autocorrelation Analysis

    In this study, the Moran I index and local indicators of spatial association (LISA) were used to analyze global and local spatial autocorrelation. The global Moran I index is a comprehensive measure of spatial autocorrelation in the whole study area, and it can indicate the average degree of spatial differences between each unit and the surrounding area. If partial positive and partial negative spatial correlations coexist in the overall study area, LISA are needed to reveal the possible spatial variability [22].

    Space-Time Scan Statistics

    For the spatiotemporal analysis, we used the Kulldorff spatial-temporal scan statistic to identify clusters based on specific time and location combinations. This statistic uses a Poisson probability distribution mode. The Kulldorff method uses a moving cylindrical window to scan the study area. The height of the cylinder represents the possible clustering time, the bottom represents the clustering area, and the radius serves as the scanning risk population [23]. The significance of the clusters identified at a 95% confidence level was tested using Monte Carlo simulations to identify the clusters.

    Prediction Models

    Three models were used in this study. The Holt-Winters exponential smoothing model eliminated some random fluctuations, which can assign different weights to the data of each period and reasonably predict the future development trend [16]. The seasonal autoregressive integrated moving average (SARIMA) model divided the observed value into 3 parts: residuals, seasonal feature, and true trend [24]. The Holt-Winters exponential smoothing model and SARIMA model were introduced more specifically in a previous study [16]. The Prophet model, developed by Meta, allowed using time as a regressor to fit both linear and nonlinear functions of time [25], and a more specific introduction to this model is provided in Xie et al’s [26] study.

    Statistical Analysis

    Descriptive analysis was conducted using Excel (Microsoft Corp) and SPSS (version 28; IBM Corp). The Moran I index was calculated using GeoDa (version 1.20), and SaTScan software (version 10.1; Harvard Medical School) was used to determine the spatial and temporal aggregation of TP. All visual results are presented using the ArcGIS software (version 10.8; SERI Inc). The prediction of TP prevalence was model-fitted using R (version 4.2.0; The R Foundation) and Python (version 3.8; Python Software Foundation).

    Ethical Considerations

    This study was approved by the Ethics Committee of the Zhejiang Provincial Center for Disease Control and Prevention (audit number 2022-032-01), which waived the requirement for informed consent. All personal information obtained in this study was treated and kept confidential in accordance with the required protocols.


    General Epidemiological Characteristics of Patients With TP

    During the study period, a total of 11,531 patients with TP were identified, with a notification rate of 7.8 cases per 100,000 population in 2017, 8.2 cases per 100,000 population in 2018, 7.7 cases per 100,000 population in 2019, 6.6 cases per 100,000 population in 2020, and 5.0 cases per 100,000 population in 2021, respectively (Figure 1).The male to female ratio was 2.18:1. Among all cases of TP, the top notified cases and notification rate were noted in 2018 (n=2594; 8.2 cases per 100,000 population). The number of cases peaked in the second quarter (n=3354) and plummeted in the fourth quarter (n=2558) during the study period, indicating seasonal distribution characteristics. In terms of age distribution, 74 (0.64%) individuals were aged 0-18 years, 6359 (55.15%) were aged 19-59 years, and 5098 (44.21%) were aged ≥60 years. In terms of geographical distribution, 2319 (20%) cases were notified in Hangzhou, while 1294 (11.22%), 1393 (12.08%), 1452 (12.59%), and 1259 (10.92%) cases were reported in Wenzhou, Ningbo, Jinhua, and Taizhou, respectively. Passive findings was the main source of TP identification in Zhejiang Province. Nearly all patients (n=11,326, 98.22%) were treated for the first time, with nearly 90% of them having completed treatment. Delayed access to medical care was observed among 6317 (54.78%) patients with TP, and diagnosis was delayed among 1462 (12.68%) patients. In addition, 86.26% (n=9949) of patients with TP were treated with a duration of ≥9 months (Table 1).

    ‎
    Figure 1. Change in the trend of tuberculous pleurisy notification rate and case numbers by month.
    Table 1. Epidemiological characteristics of patients with tuberculosis pleurisy in Zhejiang Province from 2017 to 2021 (N=11,531).
    CharacteristicPatients, n (%)
    Sex

    		Male7901 (68.52)

    		Female3630 (31.48)
    Age group (years)

    		0-1874 (0.64)

    		19-596359 (55.15)

    		≥605098 (44.21)
    City

    		Hangzhou2319 (20.11)

    		Jinhua1452 (12.59)

    		Ningbo1393 (12.08)

    		Wenzhou1294 (11.22)

    		Taizhou1259 (10.92)

    		Shaoxing1078 (9.35)

    		Quzhou858 (7.44)

    		Jiaxing674 (5.85)

    		Huzhou582 (5.05)

    		Lishui506 (4.39)

    		Zhoushan116 (1.00)
    Case finding

    		Passive finding11,513 (99.84)

    		Active finding18 (0.16)
    Treatment history

    		Initial treatment11,326 (98.22)

    		Retreatment205 (1.78)
    Anti–tuberculosis pleurisy treatment

    		Yes11,523 (99.93)

    		No8 (0.07)
    Time of delay in seeking medical treatment (days)

    		0-145214 (45.22)

    		15-292968 (25.74)

    		30-441569 (13.61)

    		45-59503 (4.36)

    		≥601261 (10.94)

    		Unknown16 (0.14)
    Time of delay in diagnosis in the health system(days)

    		0-1410,069 (87.32)

    		15-29773 (6.70)

    		30-44302 (2.62)

    		45-59122 (1.06)

    		≥60237 (2.06)

    		Unknown28 (0.24)
    Interval between diagnosis of tuberculosis pleurisy and end of treatment (days)

    		0-180734 (6.37)

    		181-270795 (6.89)

    		≥2719949 (86.28)

    		Unknown53 (0.46)
    Treatment outcome

    		Treatment completed10,422 (90.38)

    		Death203 (1.76)

    		Cure242 (2.10)

    		Changes in diagnosis204 (1.77)

    		Failure14 (0.13)

    		Adverse reactions122 (1.06)

    		Transfer to multidrug-resistant tuberculosis pleurisy treatment10 (0.09)

    		Others314 (2.73)

    Spatial Autocorrelation Analysis

    On spatial autocorrelation analysis, the global Moran I index ranged from 0.187 to 0.362, indicating spatial heterogeneity in the distribution of patients with TP in Zhejiang Province from 2017 to 2021. Furthermore, there was a positive correlation between the notification rate and spatial distribution (P<.01; Table 2). The local autocorrelation LISA cluster plot showed that in Zhejiang Province, hot spot areas were located in parts of Quzhou, Jinhua, and Hangzhou, and cold spot areas were clustered in Jiaxing, Wenzhou, and Lishui. Although the hot and cold spots changed every year, Jiangshan, Qujiang, Changshan in Quzhou City, and Chun'an in Hangzhou City emerged as hot spots throughout the study period. Furthermore, Longwan and Yueqing in Wenzhou City were identified as cold spots in 4 out of the 5 years (Figure 2).

    Table 2. Spatial autocorrelation analysis of the tuberculosis pleurisy notification rate.
    YearMoran I indexz scoreP value
    20170.3625.275.001
    20180.2964.099.001
    20190.1872.632.007
    20200.2093.026.006
    20210.2413.864.003
    ‎
    Figure 2. Local indicators of spatial association of the tuberculosis pleurisy notification rate in Zhejiang Province.

    Spatial-Temporal Scan Statistics

    Spatial-temporal scan statistics of TP in Zhejiang Province from 2010 to 2021 identified 1 most likely cluster and 4 secondary clusters. The most likely cluster was the high-risk cluster, which encompassed 32 counties in western and southwestern Zhejiang Province, including all counties and districts of Quzhou, Jinhua, and Lishui; Wencheng and Taishun in Wenzhou; and Xianju in Lishui. The most likely cluster period ranged from April 2017 to September 2019, with 2383 cases of TP and a relative risk of 1.98 (log likelihood ratio=382.55, P<.01; Table 3 and Figure 3).

    Table 3. Spatial-temporal scanning of the tuberculosis pleurisy notification rate in Zhejiang Province, 2017-2021.
    Cluster typeCluster periodCoordinates/radiusCounties, nNotified (expected) cases, n (n)Log likelihood ratioRelative riskP value
    Most likely clusterJanuary 20, 2017, to September 30, 2019(28.578400° N, 118.595300° E)/199.19 km322383 (1341.84)382.551.98<.01
    Secondary cluster 1January 1, 2017, to June 30, 2018(30.042100° N, 120.663400° E)/51.20 km8785 (488.22)80.061.65<.01
    Secondary cluster 2February 1, 2018, to July 31, 2020(27.885900° N, 120.815300° E)/0 km148 (3.87)76.8412.46<.01
    Secondary cluster 3February 1, 2017, to November 30, 2018(28.602000° N, 121.076900° E)/0 km1124 (48.34)41.402.58<.01
    Secondary cluster 4April 1, 2017, to August 31, 2018(28.652100° N, 121.475800° E)/12.06 km2129 (74.78)16.251.73<.01
    ‎
    Figure 3. Spatial-temporal clustering of the tuberculosis pleurisy notification rate in Zhejiang Province.

    Model Predictions

    On comparing the fitness of predicting TP prevalence trends in 2021, the Holt-Winters exponential smoothing model showed better predictive performance than the SARIMA and Prophet models, with lower mean-square error, root-mean-square error, mean absolute error, mean absolute percentage error, and symmetric mean absolute percentage error values (Table 4). The composition of the original sequence and the predictive results determined using the Holt-Winters exponential smoothing model are presented in Table 4. The exponential smoothing model was used to predict the number of notified cases of TP, and its actual values in 2021 are shown in Table 5.

    Table 4. Critical indices of the Holt-Winters exponential smoothing, SARIMAa, and Prophet models.
    IndicesSARIMAHolt-WintersProphet
    MSEb978.593479.170579.114
    RMSEc31.28221.89024.065
    MAEd25.24717.38017.501
    MAPEe16.27011.130N/Af
    SMAPEg16.88312.080N/A

    aSARIMA: seasonal autoregressive integrated moving average.

    bMSE: mean-square error.

    cRMSE: root-mean-square error.

    dMAE: mean absolute error.

    eMAPE: mean absolute percentage error.

    fN/A: not applicable.

    gSMAPE: symmetric mean absolute percentage error.

    Table 5. Holt-Winters exponential smoothing model used to predict the number of notified cases of tuberculosis pleurisy and its actual values in 2021.
    MonthActual valueEstimate95% CI
    January158141104-178
    February122119105-133
    March191162161-163
    April159176134-217
    May165175140-210
    June155156116-197
    July163156132-181
    August14412999-159
    September145128116-140
    October13312082-157
    November145115108-121
    December160108107-108

    Principal Findings

    TP, as a subgroup of PTB in China, has attracted increasing attention in recent years due to its severity. According to a bibliometric study, there has been an increasing number of TP-related studies in the past decade, and China had the most studies in this field [27]; furthermore, Mexico, the Philippines, Vietnam, China, and India had the highest prevalences of TP worldwide [5]. Zhejiang Province is a developed coastal area of China, and the notification rate of TP accounts for a certain proportion. Therefore, it is necessary to understand cluster epidemics at the spatial and temporal levels.

    In this study, epidemiological characteristics showed that men were more likely to develop TP than women, which is in line with the results of previous studies [28,29]. Given that TP occurs secondary to PTB, we speculated that possible explanations for sex differences were attributed to lifestyle and genetic susceptibility, such as smoking, alcohol use, drug abuse, sex hormones, and genetic variations [30,31]. Therefore, it is necessary to consider men as the key population for prevention and control, such as by strengthening screening and improving the awareness of TP in this specific population. Furthermore, we found that the TP notification rate peaked mainly in the summer. As shown in the literature, when the mean ambient temperature increased from 15.1 °C to 24.5 °C, the cumulative adjusted relative risk of hospitalization for TB was 1.07 (95% CI 1.00-1.14) [32]. Although there is no clear evidence showing a direct correlation between TP and temperature, we postulate that due to the increased temperature, there is a likely increase in the synthesis of mycolic acids, which are vital for the growth of M tuberculosis and contribute to its subverting of and escaping from the immune system [33].

    In our study, 6317 (54.78%) patients with TP experienced delayed access to medical care. A literature analysis of TB cases from 32 provinces in China showed that patient delays were also common among TB cases, with a median of 20 (IQR 6-46) days [21]. We speculated that a comparatively long patient delay may result in further invasion of the pleura. Furthermore, the median delay in PTB diagnosis in China was 1 (IQR 0-8) day, while a delay of TP diagnosis of more than 14 days accounted for 12.68% of cases in our study [21]. To our knowledge, a definite diagnosis of TP was made by detecting M tuberculosis from pleural effusion or pleural tissue [34]. However, culturing M tuberculosis takes 2-8 weeks to obtain results in some regions without access to rapid diagnosis tools, which could influence effective medical interventions [35]. Consequently, further delay in anti-TB treatment for TP might result in pleural thickening or TB empyema that requires surgical resolution [36]. Considering the long duration of the disease, insidious onset, and lack of specificity of diagnostic methods for TP [37,38], it is prone to both abovementioned delays, which may further lead to increased transmission of TB [39]. Therefore, timely detection, such as using x-ray imaging or computed tomography, and adequate health education are needed to identify TP in the early stages and improve awareness of TB information [40]. For patients with normal immune function, the current recommended treatment duration is 6 months for isoniazid and rifampicin with 2 months for ethambutol and pyrazinamide [8]. In China, the National Tuberculosis Program guidelines specifically recommend that the treatment period for TP be ≥9 months, while nearly 14% of patients had received insufficient treatment. Therefore, more attention must be paid to the treatment duration of TP.

    Spatial-temporal analysis is commonly used to explore disease clusters and identify areas at high risk of infectious disease [41]. Spatial analysis indicated spatial heterogeneity among patients with TP, with hot spots located mainly in the western cities of Zhejiang Province. The relatively large agricultural population and relatively low economic status in these areas, along with limited knowledge related to TP, are the common risk factors for the development of TP [42]. Therefore, the early detection strategy of TP in this population should be strengthened by implementing modern approaches, such as large-scale active screening. Furthermore, the spatial-temporal analysis indicated that the most likely clusters were mainly concentrated in the western part of Zhejiang Province, which were similar to the results of spatial analysis. This result unveils the distribution pattern of TP incidence in Zhejiang Province and provides a theoretical basis for the regional rational allocation of health resources related to TP and effective implementation of comprehensive prevention and control measures [41]. Therefore, exploring the causes and mechanisms leading to the regional clustering of TP should be carried out in the future.

    Each predictive model has its own limitations and merits due to the different prerequisites and mechanisms involved in predicting epidemics. In this study, our exponential smoothing model has a good prediction effect, and the actual value of 66.7% falls within the 95% CI of the predicted value. Part of the actual value that is not within the prediction range may be unsatisfactory for clinical applications. The possible reason was the impact of the COVID-19 pandemic, resulting in the deviation of the model prediction.

    Limitations

    However, this study also has some limitations. First, TP is not easy to diagnose, and it was not mandatory to include TBIMS before 2018, which could cause an underestimation of TP prevalence from our available surveillance system. Second, the combination and breakdown of different areas of the panel data makes it impossible to carry out detailed analysis in some areas. Third, the prediction model did not include other meteorological factors such as humidity and wind. Last, some potential influencing factors were not considered, such as HIV coinfection, which would be considered in the future.

    Conclusions

    In this study, the male population had a high notification rate of TP, and a high notification rate of this disease was noted during summer. The western region of Zhejiang Province has the highest relative risk of TP. Comprehensive interventions, such as chest x-ray screening and symptom screening, should be strengthened to improve early identification. Additionally, a more systematic assessment of the prevalence trend of TP should include more predictors.

    Acknowledgments

    We appreciate the local centers for disease control and prevention, community health care centers, and tuberculosis (TB)-designated hospitals in Zhejiang Province to implement TB, including tuberculosis pleurisy, reporting and notification during daily work. This study was supported by the National-Zhejiang Health Commission Major S&T Project (grant WKJ-ZJ-2118), Zhejiang Provincial Medical and Health Project (2021KY618 and 2020KY520).

    Data Availability

    All data and materials were included in this paper. The corresponding author (JJ) can provide data upon reasonable request after all studies and substudies have been completed.

    Authors' Contributions

    JP, BC, and KL conceptualized and designed the study. ZL and JJ managed the study. DL and YZ analyzed the data. SC and BC interpreted the study results. YZ, JJ, and DL drafted the manuscript for important intellectual content. KL, JP, and JJ reviewed and edited the manuscript. YZ, DL, KL, BC, SC, JP, ZL, and JJ approved the final version of the manuscript to be published.

    Conflicts of Interest

    None declared.

    1. Global Tuberculosis Report 2022. World Health Organization. 2022. URL: https://www.who.int/teams/global-tuberculosis-programme/tb-reports/global-tuberculosis-report-2022 [accessed 2023-10-05]
    2. Boonsarngsuk V, Mangkang K, Santanirand P. Prevalence and risk factors of drug-resistant extrapulmonary tuberculosis. Clin Respir J. Jun 2018;12(6):2101-2109. [CrossRef] [Medline]
    3. Ryan H, Yoo J, Darsini P. Corticosteroids for tuberculous pleurisy. Cochrane Database Syst Rev. Mar 14, 2017;3(3):CD001876. [FREE Full text] [CrossRef] [Medline]
    4. Jeon D. Tuberculous pleurisy: an update. Tuberc Respir Dis (Seoul). Apr 2014;76(4):153-159. [FREE Full text] [CrossRef] [Medline]
    5. Baumann MH, Nolan R, Petrini M, Lee YCG, Light RW, Schneider E. Pleural tuberculosis in the United States: incidence and drug resistance. Chest. Apr 2007;131(4):1125-1132. [CrossRef] [Medline]
    6. People’s Republic of China state health and Family Planning Commission. Tuberculosis classification (WS196—2017). Electronic J Emerg Infect Dis. Aug 31, 2018;3:191-192. [CrossRef]
    7. National Health and Family Planning Commission of the People's Republic of China. Diagnostic criteria for tuberculosis. Electronic J Emerg Infect Dis. Feb 28, 2018;3:59-61. [CrossRef]
    8. Shaw JA, Koegelenberg CFN. Pleural Tuberculosis. Clin Chest Med. Dec 2021;42(4):649-666. [CrossRef] [Medline]
    9. Frye MD, Pozsik CJ, Sahn SA. Tuberculous pleurisy is more common in AIDS than in non-AIDS patients with tuberculosis. Chest. Aug 1997;112(2):393-397. [CrossRef] [Medline]
    10. Wen P, Wei M, Han C, He Y, Wang M. Risk factors for tuberculous empyema in pleural tuberculosis patients. Sci Rep. Dec 20, 2019;9(1):19569. [FREE Full text] [CrossRef] [Medline]
    11. Chen S, Wang Y, Zhan Y, Liu C, Wang Q, Feng J, et al. The incidence of tuberculous pleurisy in mainland China from 2005 to 2018. Front Public Health. 2023;11:1180818. [FREE Full text] [CrossRef] [Medline]
    12. Mohamed AS, Levine M, Camp JW, Lund E, Yoder JS, Glickman LT, et al. Temporal patterns of human and canine Giardia infection in the United States: 2003-2009. Prev Vet Med. Feb 01, 2014;113(2):249-256. [CrossRef] [Medline]
    13. Gao S, Mioc D, Anton F, Yi X, Coleman DJ. Online GIS services for mapping and sharing disease information. Int J Health Geogr. Feb 25, 2008;7:8. [FREE Full text] [CrossRef] [Medline]
    14. Rogers DJ, Randolph SE. Studying the global distribution of infectious diseases using GIS and RS. Nat Rev Microbiol. Dec 2003;1(3):231-237. [FREE Full text] [CrossRef] [Medline]
    15. Luo Z, Jia X, Bao J, Song Z, Zhu H, Liu M, et al. A combined model of SARIMA and prophet models in forecasting AIDS incidence in Henan Province, China. Int J Environ Res Public Health. May 12, 2022;19(10) [FREE Full text] [CrossRef] [Medline]
    16. Wang S, Wei F, Li H, Wang Z, Wei P. Comparison of SARIMA model and Holt-Winters model in predicting the incidence of Sjögren's syndrome. Int J Rheum Dis. Nov 2022;25(11):1263-1269. [CrossRef] [Medline]
    17. Liu K, Li T, Vongpradith A, Wang F, Peng Y, Wang W, et al. Identification and prediction of tuberculosis in Eastern China: analyses from 10-year population-based notification data in Zhejiang Province, China. Sci Rep. May 04, 2020;10(1):7425. [FREE Full text] [CrossRef] [Medline]
    18. Liu K, Chen S, Zhang Y, Li T, Xie B, Wang W, et al. Tuberculosis burden caused by migrant population in Eastern China: evidence from notification records in Zhejiang Province during 2013-2017. BMC Infect Dis. Jan 31, 2022;22(1):109. [FREE Full text] [CrossRef] [Medline]
    19. Li T, Du X, Shewade HD, Soe KT, Zhang H. What happens to migrant tuberculosis patients who are transferred out using a web-based system in China? PLoS One. 2018;13(11):e0206580. [FREE Full text] [CrossRef] [Medline]
    20. Huang Y, Zuo L, Ultrasound Professional Committee of Tuberculosis Branch of Chinese Medical Association; Interventional Ultrasound Professional Committee of Interventional Physician Branch of Chinese Medical Doctor Association. Expert consensus on ultrasound diagnosis, classification and interventional therapy of tuberculous pleurisy (2022 Edition). Chinese Journal of Antituberculosis. Sep 10, 2022:880-897. [CrossRef]
    21. Li T, Du X, Kang J, Luo D, Liu X, Zhao Y. Patient, diagnosis, and treatment delays among tuberculosis patients before and during COVID-19 epidemic - China, 2018-2022. China CDC Wkly. Mar 24, 2023;5(12):259-265. [FREE Full text] [CrossRef] [Medline]
    22. Li M, OU J. Spatio-temporal analysis of influenza A(H1N1) in China during 2009-2013 based on GIS. Geogr Res. Nov 18, 2016:2139-2152. [CrossRef]
    23. Yamamura M, de Freitas IM, Santo Neto M, Chiaravalloti Neto F, Popolin MAP, Arroyo LH, et al. Spatial analysis of avoidable hospitalizations due to tuberculosis in Ribeirao Preto, SP, Brazil (2006-2012). Rev Saude Publica. 2016;50:20. [FREE Full text] [CrossRef] [Medline]
    24. Ndeh NT, Tesfaldet YT, Budnard J, Chuaicharoen P. The secondary outcome of public health measures amidst the COVID-19 pandemic in the spread of other respiratory infectious diseases in Thailand. Travel Med Infect Dis. 2022;48:102348. [FREE Full text] [CrossRef] [Medline]
    25. Ogallo W, Wanyana I, Tadesse GA, Wanjiru C, Akinwande V, Kabwama S, et al. Quantifying the impact of COVID-19 on essential health services: a comparison of interrupted time series analysis using Prophet and Poisson regression models. J Am Med Inform Assoc. Mar 16, 2023;30(4):634-642. [CrossRef] [Medline]
    26. Xie C, Wen H, Yang W, Cai J, Zhang P, Wu R, et al. Trend analysis and forecast of daily reported incidence of hand, foot and mouth disease in Hubei, China by Prophet model. Sci Rep. Jan 14, 2021;11(1):1445. [FREE Full text] [CrossRef] [Medline]
    27. Bian Y, Deng M, Zhang Q, Hou G. Global trends of research on tuberculous pleurisy over the past 15 years: a bibliometric analysis. Front Cell Infect Microbiol. 2022;12:937811. [FREE Full text] [CrossRef] [Medline]
    28. Tahseen S, Khanzada FM, Baloch AQ, Abbas Q, Bhutto MM, Alizai AW, et al. Extrapulmonary tuberculosis in Pakistan- A nation-wide multicenter retrospective study. PLoS One. 2020;15(4):e0232134. [FREE Full text] [CrossRef] [Medline]
    29. Luo Y, Yan F, Xue Y, Mao L, Lin Q, Tang G, et al. Diagnostic utility of pleural fluid T-SPOT and interferon-gamma for tuberculous pleurisy: a two-center prospective cohort study in China. Int J Infect Dis. Oct 2020;99:515-521. [FREE Full text] [CrossRef] [Medline]
    30. Khan AH, Sulaiman SAS, Hassali MA, Khan KU, Ming LC, Mateen O, et al. Effect of smoking on treatment outcome among tuberculosis patients in Malaysia; a multicenter study. BMC Public Health. Jun 04, 2020;20(1):854. [FREE Full text] [CrossRef] [Medline]
    31. Kirenga BJ, Ssengooba W, Muwonge C, Nakiyingi L, Kyaligonza S, Kasozi S, et al. Tuberculosis risk factors among tuberculosis patients in Kampala, Uganda: implications for tuberculosis control. BMC Public Health. Jan 21, 2015;15:13. [FREE Full text] [CrossRef] [Medline]
    32. Chong KC, Yeoh EK, Leung CC, Lau SYF, Lam HCY, Goggins WB, et al. Independent effect of weather, air pollutants, and seasonal influenza on risk of tuberculosis hospitalization: an analysis of 22-year hospital admission data. Sci Total Environ. Sep 01, 2022;837:155711. [CrossRef] [Medline]
    33. Takayama K, Armstrong EL, Davidson LA, Kunugi KA, Kilburn JO. Effect of low temperature on growth, viability, and synthesis of mycolic acids of Mycobacterium tuberculosis strain H37Ra. Am Rev Respir Dis. Jul 1978;118(1):113-117. [CrossRef] [Medline]
    34. Jeon D. Tuberculous pleurisy: an update. Tuberc Respir Dis (Seoul). Apr 2014;76(4):153-159. [FREE Full text] [CrossRef] [Medline]
    35. Ritchie SR, Harrison AC, Vaughan RH, Calder L, Morris AJ. New recommendations for duration of respiratory isolation based on time to detect Mycobacterium tuberculosis in liquid culture. Eur Respir J. Sep 2007;30(3):501-507. [FREE Full text] [CrossRef] [Medline]
    36. Barbas CS, Cukier A, de Varvalho CR, Barbas Filho JV, Light RW. The relationship between pleural fluid findings and the development of pleural thickening in patients with pleural tuberculosis. Chest. Nov 1991;100(5):1264-1267. [CrossRef] [Medline]
    37. Liu YJ, Jiang Z, Chen H, Jing H, Cao X, Coia JE, et al. Description of demographic and clinical characteristics of extrapulmonary tuberculosis in Shandong, China. Hippokratia. 2020;24(1):27-32. [FREE Full text] [Medline]
    38. Lee C, Chiu L, Chang C, Chung F, Li S, Chou C, et al. The clinical experience of Mycobacterial culture yield of pleural tissue by pleuroscopic pleural biopsy among tuberculous pleurisy patients. Medicina (Kaunas). Sep 15, 2022;58(9) [FREE Full text] [CrossRef] [Medline]
    39. Zhongwei J, Zuhong L. Effect of delayed diagnosis and treatment on transmission of tuberculosis. J Syst Sci Math Sci. Jan 15, 2016;36(1):28.
    40. Jiang J, Zhang X, Wang F. Influencing factors for the delay in seeking health care of pulmonary tuberculosis in the elderly. Prev Med. 2020;6:612-614. [CrossRef]
    41. Qian Y, Kun C. Spatial distribution patterns of pulmonary tuberculosis incidence in Zhejiang province: spatial autocorrelation analysis. Chin J Public Health. Apr 30, 2013;29(4):485-487.
    42. Gui JJ, Zhang TF, Liu ZF. Epidemiological characteristics and spatial clusters of pulmonary tuberculosis in Zhejiang province,2005-2011. Chin J Public Health. 2016;32(1):11-14. [CrossRef]

    ‎
    LISA: local indicators of spatial association
    PTB: pulmonary tuberculosis
    SARIMA: seasonal autoregressive integrated moving average
    TB: tuberculosis
    TP: tuberculous pleurisy

    Edited by A Mavragani, T Sanchez; submitted 12.06.23; peer-reviewed by T Li, W Tang, Y Liu; comments to author 19.08.23; revised version received 31.08.23; accepted 19.09.23; published 30.10.23.

    Copyright

    ©Ying Zhou, Dan Luo, Kui Liu, Bin Chen, Songhua Chen, Junhang Pan, Zhengwei Liu, Jianmin Jiang. Originally published in JMIR Public Health and Surveillance (https://publichealth.jmir.org), 30.10.2023.

    This is an open-access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work, first published in JMIR Public Health and Surveillance, is properly cited. The complete bibliographic information, a link to the original publication on https://publichealth.jmir.org, as well as this copyright and license information must be included.