Skip to main content

The positive externalities of migrant-based TB control strategy in a Chinese urban population with internal migration: a transmission-dynamic modeling study

Abstract

Background

Large-scale rural-to-urban migration has changed the epidemiology of tuberculosis (TB) in large Chinese cities. We estimated the contribution of TB importation, reactivation of latent infection, and local transmission to new TB cases in Shanghai, and compared the potential impact of intervention options.

Methods

We developed a transmission dynamic model of TB for Songjiang District, Shanghai, which has experienced high migration over the past 25 years. We calibrated the model to local demographic data, TB notifications, and molecular epidemiologic studies. We estimated epidemiological drivers as well as future outcomes of current TB policies and compared this base-case scenario with scenarios describing additional targeted interventions focusing on migrants or vulnerable residents.

Results

The model captured key demographic and epidemiological features of TB among migrant and resident populations in Songjiang District, Shanghai. Between 2020 and 2035, we estimate that over 60% of TB cases will occur among migrants and that approximately 43% of these cases will result from recent infection. While TB incidence will decline under current policies, we estimate that additional interventions—including active screening and preventive treatment for migrants—could reduce TB incidence by an additional 20% by 2035.

Conclusions

Migrant-focused TB interventions could produce meaningful health benefits for migrants, as well as for young residents who receive indirect protection as a result of reduced TB transmission in Shanghai. Further studies to measure cost-effectiveness are needed to evaluate the feasibility of these interventions in Shanghai and similar urban centers experiencing high migration volumes.

Peer Review reports

Background

China has an estimated 4.5 million prevalent tuberculosis (TB) cases, of which 80% arise among the rural population [1]. TB incidence in China has declined in recent decades, with national scale-up of directly observed treatment strategies and ongoing improvements in living standards. Currently, rural-to-urban migration is a major impediment for TB control in China. Most internal migrants are men who leave rural areas to join the urban wage economy, and there will be an estimated 240–260 million urban migrants by 2030. Due to China’s household registration system, most migrants typically face difficulties in becoming formal urban residents [2]. Consequently, they may not be entitled to subsidized housing or education, have poor access to social security and medical benefits, and often live and work in circumstances that promote Mycobacterium tuberculosis (Mtb) transmission and impede prompt diagnosis [3,4,5]. Studies have shown that urban migrants who develop TB have reduced access to TB care, and worse treatment outcomes than urban residents [6,7,8].

In China, identification of TB cases is primarily through passive case-finding, with symptomatic individuals voluntarily seeking diagnosis and treatment. High volume internal migration has been linked to elevated infectious disease incidence in Chinese cities, and challenges for disease control [6,7,8]. In recent years, rising TB case numbers associated with urban migration have been reported in many major cities—including Beijing, Shanghai, Shenzhen, Guangzhou, Wuhan and Hangzhou [4, 9,10,11,12,13]—yet there was no specific guideline for TB screening among urban migrants, and the underlying epidemiological mechanisms have not been completely described. The optimal response from urban TB control agencies will likely depend on the sources of new TB cases—whether elevated TB case notifications among urban migrants are due to reactivation of latent TB infection (LTBI) acquired before migration, transmission following migration, or prevalent TB among new arrivals.

In this study, we used a transmission-dynamic model to investigate the epidemiology of migration-linked TB in Shanghai. We conducted our study in Songjiang, a district of Shanghai in which the majority of the 2016 population were internal migrants. For Songjiang, previous genotypic analyses of clinical isolates have described transmission linkages between TB cases among residents and migrants [14]. Using these data, alongside detailed demographic and clinical data for 2000–2018, we estimated the contribution of different epidemiological mechanisms to TB incidence among migrants and residents and assessed the potential impact of new interventions focused on these population groups.

Methods

Study setting

Shanghai has experienced major rural-to-urban migration since the 1990s (Additional file 1: Figure S1) [15]. By 2016 the Shanghai population totaled 24 million, 41% of whom were internal migrants [16]. Songjiang has a population of 1.77 million, of which 62% are internal migrants, defined as people without a Shanghai household registration status through Chinese hukou system. Most originate from western and central China. In a previous study during 2006 to 2008, migrants had higher tuberculosis rates (38.9/100,000 population) compared to local residents (27.8/100,000 population) in Songjiang (Additional file 1: Figure S2) [6]. Since the 1990s, Songjiang has participated in a comprehensive TB surveillance and reporting system involving Songjiang District Centers for Diseases Control and Prevention, a TB-designated hospital, and community-based clinics. This network provides TB services including passive case finding, diagnosis, treatment, and follow-up. Data on TB cases, including demographic, epidemiological, and clinical information, is recorded in an electronic TB Information Management System.

Currently, Songjiang applies a referral system based on TB-like symptoms to prompt TB diagnosis and notification. Although there was no specific national guideline for TB control among internal migrants, since late 2004, an enhanced program was implemented in Shanghai to extend free TB services (reimbursement for diagnosis and anti-TB drugs) to the migrant population. Despite this expansion of standard TB care policy, there was no specific routine TB screening among migrants in Songjiang. Many local manufactory industries request a general pre-employment physical examination for both migrant and resident applicants, including chest X-ray. Residents aged 65 years and above can receive a voluntary annual physical examination, including a chest X-ray examination (not specific for TB screening).

Study design and model structure

We developed an age-structured transmission-dynamic model of TB epidemiology in Songjiang District, adapted from published models (Fig. 1) [17, 18]. The model stratifies the Songjiang population into four dimensions, with a core dimension representing TB transmission, natural history, diagnosis, and treatment (Fig. 1a and Fig. S4). Additional dimensions represented differences in population demography, TB burden and susceptibility, mixing patterns, and TB care access by age, sex, and migrant status (Fig. 1b–d). Migrant populations were stratified according to the duration of stay in Shanghai (Fig. 1d), with immigration and emigration modeled explicitly. Resident individuals entering Songjiang from other districts of Shanghai were not considered “migrants” for the analysis, as their epidemiological profile is expected to be similar to those of current residents and they held the same household registration status as the Songjiang local residents. We did not simulate transmission dynamics in rural areas and calibrated TB burden in incoming migrants to match observed disease patterns in Songjiang. The full model structure results from the crossing of the four dimensions and is represented with many distinct compartments and processes and the model code was detailed in Additional file 2.

Fig. 1
figure1

Schematic of TB simulation model. The schematic shows compartments within each dimension of the model: a Main model: core TB dimensions; b Age group dimension; c Gender dimension; d Risk group dimension by migration status. Arrows identify entries to the model and possible transitions between compartments. Solid arrows indicate state transitions, dotted arrows indicate model entry for resident individuals, and dashed arrows indicate model entry for internal migrant individuals. Migrants were divided into different status based on their duration of staying in Shanghai. Migrants exits due to death are not shown

Data inputs

Demographic data

We extracted population data from the most recent National Census (2010), Shanghai 1% Population Sampling Survey (2015), and annual statistical yearbooks of Songjiang (2000–2018) [15, 16, 19]. These data included annual population by age and migrant status, distribution of the migrant population by reported years since entry and occupational category, birth rates, and mortality rates (Additional file 1:Table S1) [15]. Future migration estimates were based on Shanghai government population policy, which targets limited future population growth [20].

Observed TB cases and outcomes

From the TB surveillance system, we extracted TB case counts (available from 2000) stratified by age, gender, and migrant status. These data were used to build historical time-series of TB notification rates. We compiled additional evidence on TB care (including TB treatment discontinuation and treatment outcomes) from administrate data sources [6].

Transmission patterns

We used genotyping data that utilized an optimized MIRU-VNTR loci set (including the hypervariable loci) to estimate the extent of shared respiratory contacts among residents and migrants in Songjiang [14]. These data represent culture-confirmed TB patients reported in Songjiang from 2009 to 2015 and were used to estimate the probability that a clustered TB case in one of the two migrant strata (residents vs. migrants) would be due to transmission from an individual of the other group. Age-based mixing patterns were based on published contact matrices [21], and we allowed for increased respiratory contacts among 20–30-year-olds to account for elevated occupational exposure.

Natural history of latent infection and active disease

We specified parameters describing the probability of infection given Mtb exposure, rates of progression to TB for infected individuals, partial immunity of previously-exposed individuals, progression from early to late-stage disease, smear positivity (an indicator of infectiousness and mortality risk), infectiousness, TB-specific mortality and self-cure rates for individuals with active disease based on the published literature, and stratified by age and sex. Passive case detection was assumed to identify individuals with advanced TB for treatment, who could subsequently die, be cured, or return to untreated TB following treatment discontinuation or failure. We assumed that all diagnosed TB cases are reported, and that there was no treatment of LTBI under the base-case scenario. Apart from differential infection risks and TB treatment access, we assumed that natural history parameters would be the same between migrants and residents, conditional on age and sex. Tables S2-S3 (Additional file 1) provide parameter definitions, values, and data sources [6, 15, 22,23,24,25,26,27,28,29,30,31,32].

Scenarios for future TB control options

The base-case scenario assumed continuation of current TB control activities, represented by current rates of treatment-seeking and current TB care outcomes. We compared this to hypothetical scenarios describing potential expansions in TB control. Table 1 describes model scenarios and operationalization [33,34,35,36,37,38,39,40]. In brief, scenarios one and two focus on improving TB prevention and care for elderly residents; scenarios three, four, and five target improvements in TB prevention and care for migrants. Scenario six represents the combination of migrant-focus interventions, and scenario seven represents the combination of all interventions included in other scenarios. We assumed any changes to current services would be introduced over 3 years, 2020–2022, then remain in place for the rest of the analysis period.

Table 1 Detailed model scenarios and operationalization

Modeled outcomes

Outcomes included the number of effective Mtb transmissions by individuals with TB, the number of new infections resulting from transmission (these two outcomes respectively describing who transmits TB, and who is infected), incident TB cases, TB deaths, and life-years lost due to TB. We also decomposed new TB cases into those due to reactivation of established LTBI (infected > 2 years ago), recent transmission (infected < 2 years), and entry of individuals into Songjiang with prevalent disease. We reported outcomes per 100,000 of population, stratified by age and migrant status. We estimated outcomes until 2035 and reported results for 2020, 2025, and 2035.

Statistics and operationalization

The model was programmed in R and C++ [41]. We used Bayesian evidence synthesis [42] to combine data and calibrate parameters to reproduce observed time-series of population and TB cases by age, sex, and migrant status. Calibration was implemented using Incremental Mixture Importance Sampling [43, 44], which produces a large number of epidemiological trajectories consistent with prior parameter distributions and calibration data (Additional file 1: Tables S4–5) [14]. We used the fitted model to simulate future outcomes under each analytic scenario. We calculated point estimates as the mean value across all simulated epidemiological trajectories and calculated 95% uncertainty intervals of the distribution of simulation results.

Sensitivity analyses

In sensitivity analyses, we tested the impact of alternative assumptions of a 5% annual increase and 5% annual decrease in future migrant population size, which was held flat in the main analysis. Secondly, we tested alternative assumptions of 0% and 6% annual decline in LTBI and TB prevalence among future entering migrants, as compared to the 3% decline assumed in the main analysis. Thirdly, we assessed the impact of an extended 30-year time horizon (to 2050), as compared to 15 years used for the main analysis. Table S6 in the additional file 1 describes the additional parameters used for model sensitivity analysis [33,34,35,36,37,38,39,40].

Results

Comparison of model estimates to data on recent population and TB case trends

Recent population data show rapid population growth (from 0.1 to 1.1 million between 2000 and 2014, Fig. 2a), and a much younger age distribution (Fig. 2b) for migrants in Songjiang, relative to the resident population. Similar trends were observed in the distribution of TB notifications over time (Fig. 2c) and age (Fig. 2e, f), with migrant cases clustered in the 20–35 age group while resident cases demonstrated peaks at 20–35 and 60–80 years-old. By 2016 migrants represented approximately 60% of the overall population and 74% of all reported TB cases in Songjiang, had a higher notification rate (43.5 per 100,000) compared to that among resident (24.3 per 100,000), which were consistent with previous reports [6, 14]. Additional detail on these demographic and epidemiological trends is shown in Figs. S1–3 [6, 15, 16]. The model was able to reproduce these and other recent epidemiological data, including the TB case fatality estimates by age (Fig. 2d).

Fig. 2
figure2

Comparison of model estimates and data describing population and TB outcomes by age and years: migrant and resident population (a, b), the annual number of diagnosed TB cases (c, e, and f), and case fatality rate (d). Model results are plotted for 1000 calibrated parameter sets to represent uncertainty in modeled results. The population size and age distribution data were abstracted from the national census and yearbook statistic records. TB case and facility data were drawn from the routine TB case reporting system

Model estimates of TB epidemiology

Figure 3 presents estimates for five key outcomes in 2020—effective Mtb transmissions by individuals with TB, new infections, incident TB cases, TB deaths, and years of life lost (YLL) from TB—stratified by age group and migrant status. Migrants in the 20–35-year-old age group represent a major fraction of each of these outcomes, reflecting their share of the population as well as elevated TB burden. Migrants were estimated to contribute 87.0% [95% uncertainty interval: 83.7, 89.9] of all instanced of TB transmission, and represent a smaller fraction (78.6% [74.4, 82.3]) of those newly infected, reflecting greater migrant-to-resident than resident-to-migrant transmission. Transmission from migrants was estimated to produce 91.7% [89.6, 93.9] of new infections among migrants and 69.9% [65.0, 73.6] of new infections among residents. The estimated fraction of all new infections caused by an individual of the other group (migrant to resident or resident to migrant) fitted well with the observed value (model 20.5% [18.8, 22.1], data 21.2%, supplementary materials). Residents represented an increased share of incident TB cases and 47.4% [42.3, 52.1] of TB deaths, reflecting elevated rates of reactivation TB and higher TB case fatality in the elderly resident population. Due to the younger average age, there was an estimated 35.2 [33.8, 36.7] YLL per TB death among migrants, compared to 16.5 [15.6, 17.4] YLL among residents. Migrants represented 70.2% [66.0, 74.4] of all life-years lost due to TB death.

Fig. 3
figure3

Model estimated key epidemiological outcomes in 2020 by age groups and migration status. Line represents the mean value and shaded areas represent 95% credible intervals. a The number of instances of Mtb transmission. b The number of new Mtb infections. c The number of new TB cases. d Deaths due to TB. e Life-years lost to TB

Figure 4 reports the fraction of new TB cases in 2020 due to LTBI reactivation, recent Mtb infection, and prevalent TB among individuals migrating into Songjiang, stratified by age and migrant status. For both migrants and residents, LTBI reactivation was estimated to represent the large majority of TB cases for individuals > 50 years old and was an important source of new TB cases, representing 72.2% [67.9, 76.8] of resident TB cases and 48.0% [42.9, 54.3] of migrant TB cases. Recent infection contributed 26.9% [23.2, 32.1] of resident TB cases and 42.9% [38.6, 45.8] of migrant TB cases, and was the primary cause of TB among individuals < 30 years old. Prevalent TB at entry to Songjiang contributed only 0.8% [0.5, 1.4] and 9.1% [7.0, 11.3] of TB cases among residents and migrants respectively. Overall, recent infection was estimated to cause 38.0% [33.3, 42.3] of incident TB cases overall.

Fig. 4
figure4

Model estimated fraction of new TB cases in 2020. The color represents the new TB cases due to LTBI reactivation (green), recent Mtb infection (orange), and prevalent TB among individuals migrating into Songjiang (blue), by age and migrant status (panel a, resident TB; panel b, migrant TB). The blue region in panel a represents residents who enter Songjiang from other districts in Shanghai. The polygon shows the distribution of new TB cases by age and the histogram presents the average of new TB cases from the three contributors. Faded color and border line represent the 95% uncertainty intervals of fraction

Epidemiological projections

Figure 5 shows base-case projections of population and TB incidence. While the total population is projected to stabilize following rapid growth during 2000–2015 (Fig. 5a), TB cases are projected to decline. TB cases among the migrant population are projected to decline to 18.4 [17.1, 23.2] per 100,000 by 2035, a 49% [36–53] reduction from 2016. By 2035 resident TB cases are projected to decline to 11.6 [10.6, 14.4] per 100,000 (Fig. 5c), a 51% [39–55] reduction from 2016. These declines are consistent with recent TB case trends (Fig. 5, points). Migrants are projected to represent a relatively stable share of total TB cases, equivalent to 63% of total TB cases while representing 53% of the 2035 population.

Fig. 5
figure5

Projected trends in population size and TB incidence under base-case scenario in Songjiang, 2010–2035. Lines represent best estimate projection for a given population group. Dots represent the reported data. Shaded areas represent 95% credible intervals. a Migrant and resident population per thousand. b TB new cases. c Number of new TB cases per million

Table 2 provides base-case projections for other TB outcomes. These outcomes are projected to follow similar trends to new TB cases, with TB transmission, new infections, TB deaths, and YLL due to TB all declining log-linearly over the projection period.

Table 2 Model estimated TB epidemiologic outcomes for the total population, migrants, and residents under base-case scenarios in 2015, 2020, 2025, and 2035

Alternative scenarios for interventions among residents and migrants

Figure 6 shows incremental outcomes for the targeted interventions described by scenarios 1–5 compared with the base-case, summed over 2020–2035, and stratified by age and migrant status. Among these scenarios, interventions directed at the migrant population (particularly scenarios 4 and 5) are projected to have the greatest impact on most outcomes, with benefits concentrated in younger age groups consistent with the migrant age distribution. An exception to this pattern is for TB deaths, where scenario 2 is projected to have a major impact on TB deaths in the elderly resident population. In all cases, interventions targeted at migrants are projected to produce beneficial secondary effect for residents. In contrast, interventions targeted at residents are projected to have minimal impact on migrant TB outcomes.

Fig. 6
figure6

Model estimated incremental changes of key epidemiological outcomes for interventions scenarios 1–5 compared with the base-case from 2020 to 2035, stratified by age and migrant status. Resident outcomes (a, c, e, g, i) and migrant outcomes (b, d, f, h, j). Lines represent best estimate incremental changes for a given population group, scenario, and age group

Table 3 summarizes incremental reductions in TB outcomes for each scenario by 2035. In general, outcomes related to TB transmission (Mtb transmission, new Mtb infections) are more sensitive to intervention than other outcomes, with the combination of all interventions (scenario 7) projected to reduce these outcomes by 20–23% by 2035 compared to the base-case, while estimated reductions in other outcomes ranged from 10 to 12%.

Table 3 Model estimated cumulated incremental scenarios of TB epidemiologic outcomes by migration status as a percentage (95% CI) of reduction compared to their base-case value

Sensitivity analysis

Sensitivity analyses representing higher future migrant populations resulted in a greater effect of migrant-focused interventions relative to resident-focused interventions, and similar results were seen for sensitivity analyses representing slower declines in TB burden among future entering migrants. Even for sensitivity analyses representing faster future declines in migration volume and migrant TB burden, interventions focused on migrants still produced the greatest reductions in TB incidence and beneficial spillovers for resident populations. Results observed over a longer time horizon were consistent with the main analysis, with ongoing log-linear declines in TB burden when outcomes were projected to 2050. Over the long-term, intervention scenarios were shown to produce a step change in TB transmission, incidence and death rates, with this percentage difference with the base case increasing slowly over the extended simulation period. Detailed results for sensitivity analyses are shown in Fig. S5 and Table S7.

Discussion

Over the past decades, there has been substantial growth of cities in China, driven by massive rural-to-urban internal migration [2]. This large-scale migration of young adults seeking economic opportunities has created opportunities for infectious disease transmission and has been associated with increases in TB case notifications. This poses new challenges for TB control in many urban areas, particularly in populous mega-cities like Shanghai, Beijing, Guangzhou, Wuhan, and Hangzhou [9,10,11,12,13]. In this study, we developed a mathematical model of TB that accounts for the demographic changes produced by high volume migration, and the differential risks of TB transmission, progression, and mortality across population groups. Using detailed local data on demography, case notifications, and molecular epidemiology, we calibrated the model to observed TB epidemiology in Shanghai, which is similar to other Chinese cities in terms of migration-linked population growth and concomitant increases in TB notifications [9,10,11,12,13]. We used this model to investigate the mechanisms underlying TB epidemiology in cities like Shanghai, and estimate the health impacts of alternative intervention strategies.

We found that large influxes of young migrants (adults predominantly 20–40 years old) has resulted in a high number of TB cases, new infections, and transmissions attributable to young adults in Shanghai. While the number of TB-related deaths is approximately equal between residents and migrants, because migrant cases are typically much younger than resident cases, the YLL due to TB is substantially larger for migrants. Our analysis suggests distinct mechanisms driving TB incidence in residents and migrants: for residents, approximately 70% of disease is due to reactivation of established latent infection, whereas for migrants, almost half of all TB cases are attributable to recent infection occurring after arrival in Shanghai. This was consistent with the findings from epidemiological observation studies in Shanghai [6, 14, 45]. We estimate that importation of prevalent disease among migrants represents a relatively small contribution to the TB burden in Shanghai. Importantly, our model also suggests very similar age-specific contributions of reactivation and recent infection. The difference in the relative importance of these mechanisms for the two population groups arises as a result of the fact that migrants are a much younger subgroup than residents, but within individual age strata, the relative contribution of each mechanism is similar.

Our analyses suggest that with the current TB control strategy, TB incidence in Shanghai will continue to decline slowly, and this slow downward trend is robust to realistic changes in migration volume. Our intervention scenarios were chosen to reflect current policy alternatives being discussed in China [46]. These included interventions focused on elderly residents, who are at high risk of developing TB (scenarios 1 and 2), and interventions focused on migrants (scenarios 3–6). Interventions focused on the elderly were estimated to have little incremental impact on the overall epidemic, but could produce sizable reductions in TB cases, deaths and YLL among the elderly. In contrast, pre-employment screening for disease (scenario 3) and treatment of latent Mtb infection (scenario 4) for recent migrants were projected to reduce TB cases and deaths in young migrants and produce secondary benefits for young residents with whom these migrants are likely to mix. Similar findings were noted in a population cohort study in the UK, that pre-screening of migrants could reduce later transmission risks [47]. If routine screening of factory workers can be brought to scale (scenario 5), we estimate that similar improvements in TB outcomes could be obtained as compared with pre-employment screening. Routine screening in factories may have advantages, as these factories are geographically concentrated in industrial zones and this might facilitate implementation of these interventions [13].

Several limitations in our study should be noted. Our projections are accompanied by uncertainty arising from incomplete knowledge of the natural history of TB, the current epidemiology of TB in Shanghai, and possible future changes in migration policy and disease burden. While we used detailed population, TB notification, and molecular epidemiological data from Shanghai to constrain model behavior to reproduce observed age- and migrant-status related TB epidemiology, uncertainty remains around serval model inputs, and as long-term model projections are difficult to validate, model design choices can bias future projections while still allowing good fit to current data. An influential model input was the MIRU-VNTR genotyping data, used to calibrate mixing assumptions. While the MIRU-VNTR loci set was optimized for local Mtb transmission in Songjiang, we cannot fully exclude risks of bias caused by low discrimination power of this method, or missed migrant cases due to emigration before disease detection. In addition, our analysis focused on the potential health benefits of interventions, yet other factors must be considered before adoption. Meanwhile, the current pandemic of COVID-19 could impact on the continuity of TB care and the future projection of TB transmission [48]. Further work is needed to investigate implementation strategies and costs for these interventions, to generate estimates of cost-effectiveness and examine feasibility.

Conclusions

In summary, this analysis reveals that a substantial fraction of new TB cases in Shanghai are the result of local transmission in the city, driven by transmission involving young adults migrating for work. While this analysis was set in a relatively confined urban area—based on the rich data available in Songjiang for implementing the analysis—this epidemiological context is relevant for many large cities in China, with young adults migrating from rural areas with higher disease burden, and facing urban living conditions that allow the spread of respiratory pathogens. In this epidemiological context, we found that the introduction of interventions designed to interrupt local transmission related to migrants are likely to have greater population benefits than interventions directed toward elderly residents, the other main risk group for TB disease in this setting. This is related to both the greater transmission estimated for migrant groups and their large share of the population in this setting. While these new interventions could reduce TB incidence and death, we found that the magnitude of these benefits would be modest, with the most ambitious migrant-focused combination interventions achieving incremental reductions in TB incidence of approximately 10–15% by 2035. This suggests that current TB control strategies may be largely effective at achieving prompt treatment and limiting transmission, and greater innovation will be needed to meet high-level targets for TB control and elimination.

Availability of data and materials

All data generated or analyzed during this study are included in this published article and its supplementary information files.

Abbreviations

TB:

Tuberculosis

Mtb:

Mycobacterium tuberculosis

LTBI:

Latent tuberculosis infection

MIRU-VNTR:

Mycobacterial interspersed repetitive unit-variable number tandem repeat

YLL:

Years of life lost

References

  1. 1.

    Technical Guidance Group of the Fifth National TB Epidemiological Survey; The Office of the Fifth National TB Epidemiological Survey. The fifth national tuberculosis epidemiological survey in 2010. Chin J Antituberculosis. 2012;34(8):485–508.

    Google Scholar 

  2. 2.

    Peng X. China's demographic history and future challenges. Science. 2011;333(6042):581–7.

    CAS  Article  Google Scholar 

  3. 3.

    Wei X, Chen J, Chen P, Newell JN, Li H, Sun C, et al. Barriers to TB care for rural-to-urban migrant TB patients in Shanghai: a qualitative study. Tropical Med Int Health. 2009;14(7):754–60.

    Article  Google Scholar 

  4. 4.

    Wang W, Jiang Q, Abdullah AS, Xu B. Barriers in accessing to tuberculosis care among non-residents in Shanghai: a descriptive study of delays in diagnosis. Eur J Pub Health. 2007;17(5):419–23.

    Article  Google Scholar 

  5. 5.

    Wei X, Zou G, Yin J, et al. Providing financial incentives to rural-to-urban tuberculosis migrants in Shanghai: an intervention study. Infect Dis Poverty. 2012;1(1):9.

    Article  Google Scholar 

  6. 6.

    Shen X, Xia Z, Li XQ, et al. Tuberculosis in an urban area in China: differences between urban migrants and local residents. PLoS One. 2012;7(11):e51133.

    CAS  Article  Google Scholar 

  7. 7.

    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.

    Article  Google Scholar 

  8. 8.

    Tobe RG, Xu L, Song P, Huang Y. The rural-to-urban migrant population in China: gloomy prospects for tuberculosis control. Biosci Trends. 2011;5(6):226–30.

    Article  Google Scholar 

  9. 9.

    Li T, He XX, Chang ZR, et al. Impact of new migrant populations on the spatial distribution of tuberculosis in Beijing. Int J Tuberc Lung Dis. 2011;15(2):163–8 i-iii.

    CAS  PubMed  Google Scholar 

  10. 10.

    Li X, Yang Q, Feng B, et al. Tuberculosis infection in rural labor migrants in Shenzhen, China: emerging challenge to tuberculosis control during urbanization. Sci Rep. 2017;7(1):4457.

    Article  Google Scholar 

  11. 11.

    Li X, Li T, Tan S. Males, ages >/= 45 years, businessperson, floating population, and rural residents may be considered high-risk groups for tuberculosis infection in Guangzhou, China: a review of 136,394 tb confirmed cases. Rev Inst Med Trop Sao Paulo. 2013;55(5):366–8.

    Article  Google Scholar 

  12. 12.

    Zhu M, Wang J, Dib HH, Wang Z. Enhancing the management of cross-regional transfer of floating tuberculosis cases by active follow-up and communication. Eur J Pub Health. 2012;22(4):577–82.

    Article  Google Scholar 

  13. 13.

    Jia ZW, Jia XW, Liu YX, et al. Spatial analysis of tuberculosis cases in migrants and permanent residents, Beijing, 2000-2006. Emerg Infect Dis. 2008;14(9):1413–9.

    Article  Google Scholar 

  14. 14.

    Yang C, Lu L, Warren JL, et al. Internal migration and transmission dynamics of tuberculosis in Shanghai, China: an epidemiological, spatial, genomic analysis. Lancet Infect Dis. 2018;18(7):788–95.

    Article  Google Scholar 

  15. 15.

    HuiMing He HL. Almanac of Songjiang. Shanghai: Shanghai Lexicographical Publishing House; 2000-2018.

  16. 16.

    National Bureau of Statistics of China. 2010 Population Census. (http://www.stats.gov.cn/english/Statisticaldata/CensusData/) 2012. Accessed 23 Apr 2013.

  17. 17.

    Menzies NA, Cohen T, Hill AN, et al. Prospects for tuberculosis elimination in the United States: results of a transmission dynamic model. Am J Epidemiol. 2018;187(9):2011–20.

    Article  Google Scholar 

  18. 18.

    Menzies NA, Cohen T, Lin HH, Murray M, Salomon JA. Population health impact and cost-effectiveness of tuberculosis diagnosis with Xpert MTB/RIF: a dynamic simulation and economic evaluation. PLoS Med. 2012;9(11):e1001347.

    Article  Google Scholar 

  19. 19.

    Shanghai Municipal Statistics Bureau. Shanghai 1% Population Sampling Survey (2015). China Statistics Press 2017. ISBN: 978-7-5037-8063-9. http://cdi.cnki.net/Titles/SingleNJ?NJCode=N2017060217.

  20. 20.

    Shanghai Urban Planning and Land Resource Administration Bureau. Shanghai Master Plan 2017-2035: Striving For The Excellent Global City. Shanghai 2018; (http://www.shanghai.gov.cn/newshanghai/xxgkfj/2035004.pdf). Accessed Jan 2018.

  21. 21.

    Prem K, Cook AR, Jit M. Projecting social contact matrices in 152 countries using contact surveys and demographic data. PLoS Comput Biol. 2017;13(9):e1005697. https://doi.org/10.1371/journal.pcbi.1005697.

  22. 22.

    Dowdy DW, Chaisson RE. The persistence of tuberculosis in the age of DOTS: reassessing the effect of case detection. Bull World Health Organ. 2009;87(4):296–304.

    Article  Google Scholar 

  23. 23.

    Behr MA, Warren SA, Salamon H, et al. Transmission of Mycobacterium tuberculosis from patients smear-negative for acid-fast bacilli. Lancet. 1999;353(9151):444–9.

    CAS  Article  Google Scholar 

  24. 24.

    Tostmann A, Kik SV, Kalisvaart NA, et al. Tuberculosis transmission by patients with smear-negative pulmonary tuberculosis in a large cohort in the Netherlands. Clin Infect Dis. 2008;47(9):1135–42.

    Article  Google Scholar 

  25. 25.

    Tiemersma EW, van der Werf MJ, Borgdorff MW, Williams BG, Nagelkerke NJD. Natural History of Tuberculosis: Duration and Fatality of Untreated Pulmonary Tuberculosis in HIV Negative Patients: A Systematic Review. PLoS ONE. 2011;6(4):e17601. https://doi.org/10.1371/journal.pone.0017601.

  26. 26.

    Sutherland I. The ten-year incidence of clinical tuberculosis following “conversion” in 2550 individuals aged 14 to 19 years. TSRU Progress Report 1968; The Hague: International Union Against Tuberculosis.

  27. 27.

    Ferebee S, Mount F. Tuberculosis morbidity in a controlled trial of the prophylactic use of isoniazid among household contacts. Am Rev Respir Dis. 1962;85:490–510.

    CAS  PubMed  Google Scholar 

  28. 28.

    Andrews JR, Noubary F, Walensky RP, Cerda R, Losina E, Horsburgh CR. Risk of progression to active tuberculosis following reinfection with Mycobacterium tuberculosis. Clin Infect Dis. 2012;54(6):784–91.

    Article  Google Scholar 

  29. 29.

    Hong Kong Chest Service/Tuberculosis Research Centre, Madras/British Medical Research Council. A controlled trial of 2-month, 3-month, and 12-month regimens of chemotherapy for sputum-smear-negative pulmonary tuberculosis. Results at 60 months. Am Rev Respir Dis 1984;130(1):23–8. https://doi.org/10.1371/journal.pone.0017601.

  30. 30.

    Menzies D, Benedetti A, Paydar A, et al. Effect of duration and intermittency of rifampin on tuberculosis treatment outcomes: a systematic review and meta-analysis. PLoS Med. 2009;6(9):e1000146.

    Article  Google Scholar 

  31. 31.

    Espinal MA, Kim SJ, Suarez PG, et al. Standard short-course chemotherapy for drug-resistant tuberculosis: treatment outcomes in 6 countries. JAMA. 2000;283(19):2537–45.

    CAS  Article  Google Scholar 

  32. 32.

    Orenstein EW, Basu S, Shah NS, et al. Treatment outcomes among patients with multidrug-resistant tuberculosis: systematic review and meta-analysis. Lancet Infect Dis. 2009;9(3):153–61.

    Article  Google Scholar 

  33. 33.

    Campbell JR, Krot J, Elwood K, Cook V, Marra F. A systematic review on TST and IGRA tests used for diagnosis of LTBI in immigrants. Mol Diagn Ther. 2015;19(1):9–24.

    Article  Google Scholar 

  34. 34.

    Zhao G, Lu X, Zhou L, et al. Study on INH preventive therapy for patients with HIV/TB co-infection. J Public Health Prev Med. 2016;27(02):42–5.

    Google Scholar 

  35. 35.

    Yang Y, Zhao Q, Ning H, Huang H, Han P. Observation on the effectiveness of preventive treatment for tuberculosis infection in family close contacts. Chin J Modern Drug Appl. 2009;3(14):209–10.

    CAS  Google Scholar 

  36. 36.

    Qi G, Ma X, Zhou C, Li C, Liu F. Analysis of tuberculosis preventive therapy implementation among freshmen in universities of Tongzhou district in Beijing. China J Antituberculosis. 2011;33(01):25–7.

    Google Scholar 

  37. 37.

    Liu Y, Tu D, An Y, ZHANG L. Control in university students in Beijing: preventive therapy for tuberculosis-infected persons. China J Antituberculosis. 2005;03:139–42.

    Google Scholar 

  38. 38.

    Sterling TR, Villarino ME, Borisov AS, et al. Three months of rifapentine and isoniazid for latent tuberculosis infection. N Engl J Med. 2011;365(23):2155–66.

    CAS  Article  Google Scholar 

  39. 39.

    Smieja MJ, Marchetti CA, Cook DJ, Smaill FM. Isoniazid for preventing tuberculosis in non-HIV infected persons. Cochrane Database Syst Rev. 2000;2:CD001363.

    Google Scholar 

  40. 40.

    WHO. Systematic screening for active tuberculosis: principles and recommendations. Geneva: World Health Organazation; 2013.

    Google Scholar 

  41. 41.

    Eddelbuettel D, Francois R. Rcpp: seamless R and C plus plus integration. J Stat Softw. 2011;40(8):1–18.

    Article  Google Scholar 

  42. 42.

    Menzies NA, Soeteman DI, Pandya A, Kim JJ. Bayesian methods for calibrating health policy models: a tutorial. Pharmacoeconomics. 2017;35(6):613–24.

    Article  Google Scholar 

  43. 43.

    Raftery AE, Bao L. Estimating and projecting trends in HIV/AIDS generalized epidemics using incremental mixture importance sampling. Biometrics. 2010;66(4):1162–73.

    Article  Google Scholar 

  44. 44.

    Poole D, Raftery AE. Inference for deterministic simulation models: the Bayesian melding approach. J Am Stat Assoc. 2000;95(452):1244–55.

    Article  Google Scholar 

  45. 45.

    Yang C, Shen X, Peng Y, et al. Transmission of Mycobacterium tuberculosis in China: a population-based molecular epidemiologic study. Clin Infect Dis. 2015;61(2):219–27.

    Article  Google Scholar 

  46. 46.

    Zhang CY, Zhao F, Xia YY, et al. Prevalence and risk factors of active pulmonary tuberculosis among elderly people in China: a population based cross-sectional study. Infectious Diseases of Poverty. 2019;8(1):7.

    Article  Google Scholar 

  47. 47.

    Aldridge RW, Zenner D, White PJ, et al. Tuberculosis in migrants moving from high-incidence to low-incidence countries: a population-based cohort study of 519 955 migrants screened before entry to England, Wales, and Northern Ireland. Lancet. 2016;388(10059):2510–8.

    Article  Google Scholar 

  48. 48.

    Hogan AB, Jewell BL, Sherrard-Smith E, et al. Potential impact of the COVID-19 pandemic on HIV, tuberculosis, and malaria in low-income and middle-income countries: a modelling study. Lancet Glob Health. 2020. https://doi.org/10.1016/S2214-109X(20)30288-6.

Download references

Acknowledgements

We thank the health workers in the tuberculosis control department of Shanghai Municipal Center for Disease Control and Prevention (Shanghai CDC) and Songjiang District CDC for their general support and valuable suggestions for the model scenarios of interventions.

Funding

This work was supported by the National Natural Science Foundation of China (81872679 to XS) and grants provided by the National Institutes of Health (NIH) (R01AI146555 to NAM, R01AI13058 to TC). This work is also supported by the National Major Science and Technology Project of China (2018ZX10715-012-002 to XS, LL, and XG) and the Shanghai Municipal Project for Academic Leaders in Public Health (grant GWV-10.2-XD23 to XS). CY was supported by the Robert E. Leet and Clara Guthrie Patterson Trust Mentored Research Award and CTSA grant (UL1 TR001863, KL2 TR001862) from the National Center for Advancing Translational Science (NCATS), components of the NIH, and NIH roadmap for Medical Research. Its contents are solely the responsibility of the authors and do not necessarily represent the official view of NIH. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Author information

Affiliations

Authors

Contributions

CY, JK, TC, and NAM conceived the study. All authors contributed to the study design. JK, CY, TC, and NAM developed the model. CY, LL, XG, and XS collected and synthesized the data from Songjiang and Shanghai statistics records. JK, CY, and NAM did the data analysis. CY wrote the first draft of the manuscript. All authors contributed to the interpretation of the results, revised the manuscript for important intellectual content, and approved the final version.

Corresponding authors

Correspondence to Xin Shen or Nicolas A. Menzies.

Ethics declarations

Ethics approval and consent to participate

The study and the use of data were reviewed and approved by the Ethical Review Committee at the Shanghai Municipal Center for Disease Control and Prevention (2011630).

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

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

Correspondence: Nicolas A Menzies, nmenzies@hsph.harvard.edu, Harvard T.H. Chan School of Public Health, Boston, MA, USA; Xin Shen, shenxin@scdc.sh.cn, Shanghai Municipal Center for Disease Control and Prevention, Shanghai, China.

Supplementary Information

Additional file 1: Fig. S1.

Population by migrant status in Shanghai (1978–2016). Fig. S2. TB notification by migrant status and age group in Shanghai (2008–2016). Fig. S3. Occupations of migrants in Songjiang District, Shanghai, 2010. Fig. S4. Model structure and mechanisms. Fig. S5. Sensitivity of the estimated change of TB notification and incidence. Table S1. Birth and mortality rates of Songjiang District, Shanghai 2001–2018. Table S2. Model initial parameters – Demographics. Table S3. Model initial parameters – TB epidemiology and natural history. Table S4. Prior and posterior distributions for calibrated model parameters. Table S5. Distribution of cluster size of different cluster types, 2009–2015. Table S6. Additional parameters used for model sensitivity analysis. Table S7. Long-term horizon sensitivity analysis

Additional file 2.

model data and code.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Yang, C., Kang, J., Lu, L. et al. The positive externalities of migrant-based TB control strategy in a Chinese urban population with internal migration: a transmission-dynamic modeling study. BMC Med 19, 95 (2021). https://doi.org/10.1186/s12916-021-01968-9

Download citation

Keywords

  • Tuberculosis
  • Internal migration
  • Transmission-dynamic modeling
  • China
  • Public health policy