Analysis of Transmission and Control of Tuberculosis in Mainland China, 2005-2016, Based on the Age-Structure Mathematical Model.

Tuberculosis (TB), an air-borne infectious disease, is a major public-health problem in China. The reported number of the active tuberculosis cases is about one million each year. The morbidity data for 2005-2012 reflect that the difference in morbidity based on age group is significant, thus the role of age-structure on the transmission of TB needs to be further developed. In this work, based on the reported data and the observed morbidity characteristics, we propose a susceptible-exposed-infectious-recovered (SEIR) epidemic model with age groupings, involving three categories: children, the middle-aged, and senior to investigate the role of age on the transmission of tuberculosis in Mainland China from 2005 to 2016. Then, we evaluated the parameters by the Least Square method and simulated the model and it had good alignment with the reported infected TB data in Mainland China. Furthermore, we estimated the basic reproduction number R₀ of 1.7858, with an obtained 95% confidence interval for R₀ of (1.7752, 1.7963) by Latin hypercube sampling, and we completed a sensitivity analysis of R₀ in terms of some parameters. Our study demonstrates that diverse age groups have different effects on TB. Two effective measures were found that would help reach the goals of theWorld Health Organization (WHO) End TB Strategy: an increase in the recovery rate and the reduction in the infectious rate of the senior age group.

1. Introduction

Tuberculosis (TB) is an air-borne infective disease caused by the slowly-replicating bacterium Mycobacterium tuberculosis (Mtb). Person-to-person transmission of Mtb occurs via the respiratory system, which can happen through both close contact between people and through infectious bacilli being carried throughout buildings by air currents [1]. According to the World Health Organization’s (WHO) Global Tuberculosis Report 2013 [2], an estimated 8.6 million new cases of TB and 1.3 million deaths, including 320,000 deaths among HIV-positive people, were recorded in 2012. Approximately 80% of all new TB cases in the world occur in 22 high burden countries that have incidence rates from 59 to 1003 per 100,000 people. India and China have the largest number of cases at 26% and 12% of the global total, respectively. Despite widespread implementation of control measures, including the Bacillus Calmette Guerin (BCG) vaccination, antiretroviral therapy, antimicrobial chemotherapy, and the WHO’s short course chemotherapy strategy, TB remains a major cause of illness and death in most high-incidence countries.

As in many high-TB burden countries, TB is a major public health problem in China. Attributed to the emergence of drug-resistant strains of Mtb and the spread of HIV/AIDS, the TB burden in China has increased over the past two decades [3]. Meanwhile, according to the United Nations, as a result of the one-child policy and low mortality, China is ageing more rapidly than almost any other country [4]. From the morbidity statistics from 2005 to 2012, shown in Figure 1, we know that the difference in morbidity between different age groups is remarkable: the morbidity of the 0–14 age group is the lowest and that of the 60 and older age group is the highest, while the morbidity of the 15–59 age group is in between. However, the role of this age grouping pattern in the transmission of TB in China is not yet clear. Therefore, we wanted to identify the potential effect of age on the morbidity for the transmission of TB in China, which could help provide effective control strategies for diverse sub-populations.

Figure 1

The tuberculosis (TB) morbidity based on age structure from 2005 to 2012 in China (source: CDC, China).

Mathematical modeling has become a powerful tool for analyzing epidemiological characteristics [5,6,7,8,9]. Different models have been developed for defining target sub-populations for treating latent TB infections and incorporating certain factors, such as drug-resistant strains, co-infection with HIV, relapse, re-infection, and vaccination, to study the transmission dynamics of TB [7,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24]. In particular, Blower et al. [10] proposed a simple TB transmission model and presented a theoretical framework for assessing the intrinsic TB transmission dynamics. Bhunuet et al. [14] considered a TB model incorporated the treatment of infectives and chemoprophylaxis. Liu et al. [19] studied a TB model incorporating seasonality. Huynh et al. [24] developed an individual-based computational model to explore the trajectory of the TB burden if the DOTS strategy is maintained or if new interventions are introduced. A more detailed discussion on different TB models was completed by White and Garnett [25]. However, few works have used mathematical models with age groupings to study the transmission of TB in China. In this paper, based on the reported data and the observed morbidity characteristics, we created a susceptible-exposed-infectious-recovered (SEIR) model with age groups of childhood, middle-aged, and senior, to investigate the role of age on the transmission process and evaluate feasible control strategies to reach the goals outlined in the WHO End TB Strategy. We estimated the basic reproduction number , analyzed the globally dynamic behavior of the model, and used the model to simulate the annual data of infected TB cases reported by the Center for Disease Control (CDC) from 2005 to 2016. Finally, we completed uncertainty and sensitivity analysis of , and explored some effective and targeted control measures for the transmission of TB in China. The rest of this paper is organized as follows. In Section 2, we present the data collection; formulate the TB model; obtain the theoretical results, such as existence and uniqueness of the solution; and define the basic reproduction ratio and global stability of disease-free equilibrium. In Section 3, data fitting and sensitivity analysis of are shown, and the feasibility of the WHO End TB Strategy is assessed. A brief discussion ensues in Section 4.

2. Materials and Methods

2.1. Data Collection

The reported annual and cumulative Tuberculosis cases in Mainland China from 2005 to 2016 were obtained from the National Notifiable Disease Surveillance System (NNDSS) (Table 1). Between 2005 and 2012, a total of 14,221,317 active TB cases were reported. The number of new infected TB cases declined 4.479% per year, from 1,259,308 in 2005 to 836,236 in 2016.

Table 1

Reported infected TB cases from 2005 to 2016 in China.

Year	2005	2006	2007	2008	2009	2010
Cases	1,259,308	1,127,571	1,163,959	1,169,540	953,275	951,508
Year	2011	2012	2013	2014	2015	2016
Cases	1,076,938	991,350	889,381	864,015	864,015	836,236

2.2. Model Formulation

In this section, we introduce a deterministic TB model incorporating age grouping with control measures. The entire population is classified into four classes: susceptible (S), latency (E), infectious (I) and recovered (R). Based on the observation that the morbidity among diverse age groups is significantly different (Figure 1), to explore the role of age on the infection pattern between susceptible and infectious classes, the susceptible class was further divided into three age groups: childhood (), middle-aged (), and senior (). We also assumed that the latent, infectious, and recovered classes are the same for different age groups. Since the latent TB cases, which are individuals who have been infected by TB bacteria but are asymptomatic, and cured TB cases may not directly cause death [26], we assumed that the death rate of the latent and recovered classes were related to the natural death rate d. Additionally, for infectious class, we added the term , based on natural death rate d, to describe the deaths caused by TB infection. Our assumptions for the dynamic transmission of TB in China with age groupings are demonstrated in Figure 2.

Figure 2

Flow chart of TB transmission with age-structure.

The model we created has the compartmental structure of the classical SEIR epidemic model, and is described by the following differential equations: where all the parameters are positive. A is the annual birth rate of the population; and are the conversion rates from the susceptible children to the susceptible middle-aged group, and from the susceptible middle-aged group to the susceptible senior group, respectively; and are the morbidities of children, middle-aged, and senior susceptible age groups, respectively; p is the fraction of fast-developing infectious cases; v is the re-activation rate of the latent TB patients; and are, respectively, the mortalities of the adolescent, the middle-aged and the elderly susceptible age groups; d is the natural death rate; is the disease-induced death rate; is the recovery rate; and is the recurrence rate of successfully treated TB cases.

Due to the severity of the transmission situation, China developed and implemented two five-year national plans in the 1980s and one 10-year national plan in the 1990s to control TB. After implementing these national TB control programs, the modern TB control strategy was implemented. Subsequently, China increased high-quality directly observed treatment, short course chemotherapy (DOTS) [24], and a compulsory Bacillus Calmette Guerin (BCG) immunization program for newborns [27]. These actions helped to effectively control the increase of TB in China. Given this, and based on Model (1), we considered two kinds of control strategies for TB in China: the incremental recovery rate per year due to DOTS, (), and the immunity rate of the BCG vaccine, (). By assuming the newborns that received the BCG vaccine remain in the susceptible compartment, Model (1) becomes the following:

2.3. Theoretical Results of Model (2)

In epidemiology, the basic reproduction number (denoted ) of an infection can be viewed as the number of cases one case generates on average over the course of its infectious period [28]. This is one of the most important indexes in evaluating the risk of an infectious disease. The asymptotical dynamic behavior of infectious diseases can be reflected by the steady state, which implies the disease will die out or persist in the future. Therefore, we first provided some mathematical analysis results of Model (2), whose proofs are shown in Appendix A.

Model (2) has the following positively invariant set:

Making use of the next generation matrix (see [29]), we obtained the basic reproduction number of Model (2) as follows:

This model has a disease-free equilibrium , where

and the endemic equilibrium , which is determined by the following equations

If , the disease-free equilibrium is globally asymptotically stable.

3. Numerical Simulations and Sensitivity Analysis

Despite the central government completing two 10-year control plans, many difficulties still exist elsewhere in the country’s TB control programs. The spread of severe acute respiratory syndrome (SARS) in 2003 revealed substantial weaknesses in the country’s public health system. After the SARS epidemic was controlled, the government made better efforts to tackle public health problems, and increased public health funding, revised laws that concerned the control of infectious diseases, implemented the world’s largest Internet-based disease reporting system, and started a program to rebuild local public health facilities. These measures contributed to an acceleration in the efforts to control tuberculosis [30,31]. Because the data quality for TB is higher after 2004, we decided to fit the data for the infected TB cases for 2005–2016 in China using Model (2). The data from, 2005–2015 were used to fit and those of 2016 were used to check the predictive power by residual and statistic.

3.1. Parameter Estimation

To perform the numerical simulations, we first needed to estimate the model parameters. According to the existing literature and related results of the Chinese population statistic yearbook, we estimated the parameters. The values of the parameters are listed in Table 2, and the detailed estimation process of the parameter values are as follows.

Table 2

The parameter values of Model (2).

Parameter	Value	Confidence Interval	Source
A	1.623×10^7	(1.606×10^7,1.643×10^7)	[32]
d	0.007	(0.00688, 0.00702)	[32]
d_1	0.0017	(0.0013, 0.0021)	[32]
d_2	0.0023	(0.0023, 0.0024)	[32]
d_3	0.0367	(0.0352,0.0382)	[32]
m_1	0.079	−	Fitting
m_2	0.0067	−	Fitting
μ	0.0025	−	[2]
p	0.05	(0,0.3)	[10]
γ	0.496	−	[10]
ξ	0.51	−	[3]
v	6	−	[33]
η	0.00341	−	[30]
φ	0.9	−	[3]
λ_1	1.268×10^‒10	−	Fitting
λ_2	5.11×10^‒10	−	Fitting
λ_3	2.553×10^‒9	−	Fitting

(a) From the results in the China population statistic yearbook from 2005 to 2015 [32], we obtained the natural death rate of the entire population and those of the three age groups. Hence, the mean and the 95% confidence interval are , CI , , CI , , CI , , CI and , CI .

(b) Based on the proportions of age groups in 2005, we calculated the initial value , . The initial number of people infected with TB is , which is the number of people infected with TB in 2005. Moreover, the percent ages of people with a TB bacteria infection but asymptomatic and those successfully treated for TB are 12.1% and 80%, respectively [3]. In 2004, the number of people infected with TB was , so we assumed the initial value of and .

(c) Using the following system and the census data of total population in China from 2005 to 2015, we estimated the parameters and by nonlinear Least-Square method (see Figure 3).

Figure 3

The fitting curve of parameters and by L-S method using the total population from 2005 to 2015. The blue cross point and red dashed line represent the total number of population and the fitted curve, respectively.

(d) The latent period of TB is about two months [33], thus we calculated the re-activation rate of latent TB patients annually. From the 2013 WHO global tuberculosis report [2], we obtained the disease-induced death rate , and from Blower et al. [10], we knew the fraction of fast-developing infectious cases p is and the recovery rate is . According to the Fifth national TB epidemiological survey [3], we knew that the incremental recovery rate of TB is and is .

(e) By use of Model (2), we simulated the cumulative number of people infected with TB from 2005 to 2016. The infection rate values were obtained by the nonlinear Least-Square method. First, we let denote the cumulative number of people infected with TB at time t. According to the flow chart of TB transmission by age grouping (Figure 2), we knew that three parts contributed to the number of infectious compartments: the number of infected people from the three susceptible age group , the latency, and the TB recurrence from recovery: where represents the cumulative number of people infected with TB at time t, and denotes the number in compartment I at time t, which includes the newly-infected TB cases and recovery TB cases at time t. Thus, to estimate the newly-infected TB cases, we had represent the newly-infected TB cases. In the following, we used to simulate the reported TB infected cases per year.

3.2. Numerical Simulations from 2005 to 2016

The decrease in infected TB cases may be due to the current control strategies not being fully effective [31], which aligns with the dynamic behaviors of Model (2). China developed and implemented two five-year national plans in the 1980s and one 10-year national plan in the 1990s to control TB. After implementing these national TB control programs, the full modern TB control strategy was implemented. The increase of TB in China has since been effectively controlled.

With help from the MATLAB (The Mathworks, Inc., Natick, MA, USA) tool fminsearch, which is part of the optimization toolbox, we estimated the optimal parameters for Model (2). Then, using the fourth-order and five-order Runge-Kutta algorithm (ode45 function), which is a powerful tool for solving ordinary differential equations, according to the corresponding parameters of Model (2) listed in Table 2, we simulated the data of the cumulative number and reported cases of TB infection from 2005 to 2016. Meanwhile, by random sampling of the 95% confidence interval (CI) of the parameters, we further plotted the 95% CI of the trajectories of the TB infection data, both cumulative and newly-infected TB cases, based on 2000 independent repeated simulations of Model (2) (see Figure 4).

Figure 4

The fitting curve and observed data from 2005 to 2016 in Mainland China with Model (2). (a) The simulation result for the cumulative number of people infected with TB; (b) The simulation result for newly-infected TB cases, where the red dashed line is the simulated curve of Model (2), the light grey area is the 95% confidence interval (CI) for all 2000 simulations, and the blue circles are the reported TB infection cases.

Figure 4 shows both the time evolution of infection cases and a comparison with the empirical records of TB infection cases, and also shows the 95% percent interval for all 3000 passing simulation trajectories. Moreover, we calculated the residual of 2016 as 235 and R-square (R) statistic to show goodness of fit [34], where the R-square value is 0.9812. We also observed that the actual reported TB infection data almost fell into the 95% CI of our simulation trajectories. Thus, our simulation results are in good accordance with the reported TB infection data, both cumulative and newly-infected TB cases, from the CDC in China from 2005 to 2016. Model (2) had a better predictive performance.

In addition, to evaluate the TB burden of China based on our model, according to the definition of incidence that the number of new and relapse cases of TB arising in a given time period, usually one year, we can further translate the reported infected TB cases into the incidence rate of TB. For comparison, we also plot the global TB incidence and the estimated TB incidence of WHO from 2005 to 2015. Figure 5 shows that after 2008, the TB incidence is lower than that of estimated value by WHO, which may implies that China substantial accelerate the control effects of TB. Moreover, we can observed that the TB incidence of China is far below the global level.

Figure 5

The comparison of TB incidence from 2005 to 2015.

3.3. Uncertainty and Sensitivity Analysis of

Due to the uncertainty in the initial parameter estimates, we performed a Latin hypercube sampling (LHS) on the estimated parameters (see, e.g., [8,35,36]). Since the LHS requires assigning a probability density function (PDF) to each of the parameters, we stratified the PDFs into 3000 equiprobability areas and then independently randomly sampled 3000 times without replacement, forming 3000 input parameter vectors [21].These input parameter vectors were then used to calculate the numerical distribution of the basic reproduction number . With the simulated parameter values, we obtained the numerical distribution of the basic reproduction number (see Figure 6), and estimated the basic reproduction number from 2005 to 2016 is and the confidence interval of is .

Figure 6

The numerical distribution of the basic reproduction number .

For the sensitivity analysis of , we can calculate partial rank correlation coefficient (PRCC), which reflects the correlation between parameters and . The PRCC of the estimated parameters with respect to are listed in Table 3. It follows from Table 3 that there exist a positive correlation between and , and a negative correlation between and . Furthermore, we can obtain that , namely, play the most important role to determine .

Table 3

The PRCC of the estimated parameters with respect to .

Parameters	PRCC	p-Value
A	0.9961	0.0005
γ	−0.7482	0.0009
λ_3	0.6141	0.0017
m_2	−0.5265	0.0001
λ_1	0.5195	0.004
λ_2	0.4820	0.0022
η	0.1808	0.0036
m_1	0.0175	0.0012

3.4. Feasibility Assessment of Reaching WHO End TB Strategy

Significant progress in controlling TB has been made during the last two decades, however, the WHO proposed a post-2015 global End TB Strategy in 2014 [37]. This strategy aims to end the global TB epidemic, with targets to cut new cases by 90% by 2035 and a milestones of 50% reduction in TB incidence rate in 2025.

In the above analysis, and are the most important risk factors for TB control. To examine the TB controlling effects with respect to and , we examined if reaching the WHO End TB Strategy would be feasible based on the current different control strategies. We used the parameter values listed in Table 2 as a baseline to compare the control effects. First, we only considered the single intervention scenario including , and, as shown in Figure 7a, we would not be able to reach the goal of WHO End TB Strategy under the current plan, even with decreasing by 50%. Then, we considered the single intervention scenario of , and as shown in Figure 7b, 15% increasing of the baseline would allow us to reach the WHO target. Finally, we considered an integrated control strategy including both and simultaneously. Figure 7c shows that if we can reduce the morbidity in the senior group by 15%, and increase the recovery rate by 10%, then we will meet the TB End Target. Therefore, we concluded that, by using the current TB control interventions, China may not reach the WHO End TB Strategy in 2025. To achieve the WHO End TB Strategy goal, China will need to pay more attention to enhance their combination TB interventions and further explore the feasibility of additional control strategies.

Figure 7

The reported TB cases (a) under control strategies of decreasing from 1 to 0.5; (b) under control strategies of increasing from 1 to 1.15; and (c) under a combination of control strategies of increasing from 1 to 1.1 and decreasing from 1 to 0.85, simultaneously.

4. Discussion

The Millennium Development Goal’s target in China was achieved with the decrease in the reported number of TB cases, however, the aging demographic represents an increasing challenge to TB control as China considers its post-2015 End TB Strategy [24]. Importantly, significant differences exist among different age groups in terms of the morbidity of TB. Taking this into account, and using the reported TB data in China from 2005 to 2016, we proposed a SEIR epidemic model with three age groups, children, middle-aged, and senior, to study the transmission of tuberculosis in China. By means of the Least Square method, we evaluated the parameters and simulated the model, and the model agrees well with the annual reported TB data in China. Furthermore, we calculated the basic reproduction number , and obtain the 95% confidence interval for is about by Latin hypercube sampling. We also assessed the feasibility of reaching the WHO End TB Strategy goal under current China TB control initiatives by using a sensitivity analysis of in terms of the parameters.

(i) Our results demonstrate that taking the age grouping into consideration is reasonable to characterize the transmission and to improve the control strategies of targeting therapy for TB in China. Based on the age-structuring model, more risk factors for different age groups can be identified. Interventions could be targeted toward specific groups, which would be particularly effective as an epidemic control measure [38]. Thus, the age grouping pattern provides a meaningful scheme, based upon the treatment of active cases and the chemoprophylaxis of latently infected individuals, to define targeted sub-populations for treating TB infections. For instance, the BCG vaccine is useful only for younger people but is less effective for the middle- or the senior-aged groups, having an average efficacy of only about 50% for those groups [2,39]. However, with the aging of the Chinese population and high morbidity rate of TB in seniors, perhaps an analogue of the BCG vaccine control strategy should be implemented for the potentially high-risk senior sub-population, which may result in the decreasing the morbidity in that group. In addition, the nationwide DOTS program should be more focused on the senior-aged group, such as providing more financial assistance for this group, who may experience catastrophic costs due to TB [26], and should place more emphasis on the people with latent TB in middle-aged group, who may increase the proportion of the actively infected people in the senior group.

(ii) From the analysis of PRCC of in Table 3, it is shown that , , and are the most effective methods for controlling TB in China. Although the WHO’s target treatment levels may not lead to eradication, these non-eradication treatment levels could significantly reduce morbidity and mortality [11]. Thus, two important indexes must be improved: First, the TB treatment success rate and treatment coverage (increasing ), for example, by providing high-quality TB care to prevent suffering and death from TB. Second, monitoring and detecting the latent TB in the senior population (reducing ) may help prevent the development of active TB in those already infected with Mycobacterium tuberculosis, including further strengthening the public health facilities and providing an isolation policy for those with detected latent TB. For TB infection in children, contact tracing is one of the key components of TB prevention, so educational programming and campaigning can be aimed the youngest age group.

(iii) Our feasibility assessment of reaching WHO End TB Strategy goal for 2015–2025, showed that even with any single intervention or combination of interventions, China may not reach the goal at the country level, as shown by the multi-models result in Houben et al. [5]. Due to the influence of drug-resistant strains, co-infection with other diseases including HIV, diabetes mellitus, etc., and increasing infection opportunities that accompanies world travel, TB will be weakly persistent and should show an overall decreasing trend in the future (see Figure 7). shows that if we can reduce the morbidity of the senior group by 15%, and increase the recovery rate by 10% , then we could potentially achieve the WHO TB End Target. Similar to Wang et al. [31] pointed out, China is not on track, nor does it appear to be currently possible, to reach the required reduction in prevalence. Therefore, there is still a need for sustained improvements in TB control to keep reducing the burden of TB in China.

From a practical viewpoint, clarifying the role of age in the transmission of TB may aid in forecasting the long-term health risks, in proposing a targeted TB control strategy, and in setting objectives and using limited resources more effectively [19].