Figures
Abstract
Tuberculosis (TB) disease can be caused by either recent transmission from infectious patients or reactivation of remote latent infection. Spatial dependency (correlation between nearby geographic areas) in tuberculosis incidence is a signature for chains of recent transmission with geographic diffusion. To understand the contribution of recent transmission in the TB endemic in Taiwan, where reactivation has been assumed to be the predominant mode of pathogenesis, we used spatial regression analysis to examine whether there was spatial dependency between the TB incidence in each township and in its neighbors. A total of 90,661 TB cases from 349 townships in 2003–2008 were included in this analysis. After adjusting for the effects of confounding socioeconomic variables, including the percentages of aboriginals and average household income, the results show that the spatial lag parameter remains positively significant (0.43, p<0.001), which indicates that the TB incidences of neighboring townships had an effect on the TB incidence in each township. Townships with substantial spatial spillover effects were mainly located in the northern, western and eastern parts of Taiwan. Spatial dependency implies that recent transmission plays a significant role in the pathogenesis of TB in Taiwan. Therefore, in addition to the current focus on improving the cure rate under directly observed therapy programs, more resource need to be allocated to active case finding in order to break the chain of transmission.
Citation: Ng I-C, Wen T-H, Wang J-Y, Fang C-T (2012) Spatial Dependency of Tuberculosis Incidence in Taiwan. PLoS ONE 7(11): e50740. https://doi.org/10.1371/journal.pone.0050740
Editor: Madhukar Pai, McGill University, Canada
Received: September 9, 2012; Accepted: October 23, 2012; Published: November 30, 2012
Copyright: © 2012 Ng et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: The authors acknowledge the financial support provided by the Infectious Diseases Research and Education Center, Department of Health and the National Taiwan University. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Human tuberculosis (TB) is an airborne infectious disease caused by Mycobacterium tuberculosis. The risk of progressing to active disease is highest in the first 2 years after infection, during which half of the symptomatic TB cases occur [1]. Active TB disease can be the result of either recent transmission from infectious patients or reactivation of remote latent infection [1], [2]. Genotyping and geospatial scanning investigations by the Centers for Disease Control and Prevention (Atlanta, Georgia) have shown that approximately 1 in 4 active TB cases reported in the United States may be attributed to recent transmission [2]. Limited data are available for the role of recent transmission in nationwide TB incidence in other countries.
Taiwan is a middle-burden country with an annual TB incidence remaining around 70 per 100,000 people from 1997 through 2005 [3], despite BCG vaccination and anti-TB drug therapy. National Directly Observed Treatment (DOT) programs were started in 2006, and the annual TB incidence gradually decreased to 57 per 100,000 people in 2010 [3]. Reactivation has been assumed to be the predominant mode of pathogenesis because of age transition in tuberculosis patients – from a disease of young adults during 1957–1961 to a disease of elderly (≥65 years) people during 1997–2001 [4]. Although several outbreaks of active TB cases, which occurred in a family or within a hospital, were identified using genotyping techniques [5]–[9], there remains a lack of nationwide genotyping or geospatial investigations on the role of recent transmission in the TB endemic in Taiwan.
Because the cumulative effect of local TB transmission among communities will cause geographic diffusion, we hypothesize that, if recent transmission plays a significant role in the TB endemic in Taiwan, we should be able to observe the presence of spatial dependency (the correlation between nearby geographical areas) in TB incidence between neighboring townships after adjusting for the spatial autocorrelation of the underlying sociodemographic and ethnic factors that influence the incidence of TB and TB reactivation (i.e. age, economic status, human immunodeficiency virus (HIV) infection, and aborigines) [10], [11].
To understand the role of recent transmission in TB endemic in Taiwan, we applied spatial regression analyses to examine whether spatial dependency exists for the TB incidence at the township-level, after adjusting for the effects of socioeconomic geography.
Methods
Data Sources
Pulmonary TB is a notifiable disease that must be reported in Taiwan. Anonymized data on TB cases were obtained from the Notifiable Infectious Disease Statistics System [3] of Taiwan Centers for Disease Control (Taipei, Taiwan). Cases occurring from 2003–2008 were included in this study. The townships where TB cases occurred were mapped according to the patients’ residential addresses. Cases from the outlying islands (including Penghu County, Kinmen County, Lienchiang County, Green Island, Orchid Island, and Liu Chiu Island) were excluded. A total of 90,661 TB cases from 349 townships were included. The TB incidence of each township was estimated using the number of TB cases during a period divided by total population of the township. Demographic and socioeconomic data were obtained from the 2000 Taiwan Census. Anonymized data on HIV infection cases were also obtained from the Notifiable Infectious Disease Statistics System [12].
Ethics Statement
Taiwan Centers for Disease Control (Taipei, Taiwan) approved the use of data for the present study. The study procedure was reviewed and approved by the Institutional Review Board (IRB) of National Taiwan University Hospital (Taipei, Taiwan). The IRB approved the exemption of informed consent because the data on TB and HIV cases had been anonymized by the Notifiable Infectious Disease Statistics System.
Socioeconomic Variables
Taiwan Census data included the population density, average household income, average number of persons per household, average years of education, and percentages of the population that were elderly (>60 years), aboriginal, Southeast Asian brides, and Southeast Asian laborers for each township. The average household income and average years of education were analyzed by quartiles using dummy variables (see Table 1 for details).
Spatial Autocorrelation and the Spatial Weight Matrix
Spatial autocorrelation identifies the patterns of spatial dependency by calculating the correlation of a variable with itself within a geographic space, meaning that the value of a variable is associated with those of the same variable in nearby areas. If spatial autocorrelation exists, general statistical methods that assume values of observations are independent may be invalid for further analysis. Spatial autocorrelation can occur in two directions: positive and negative. Positive spatial autocorrelation implies that the values of neighboring areas are similar to one another, while negative autocorrelation implies they are opposed to each other. The statistic used in this study to measure spatial autocorrelation is Moran’s I. This measure is used for variables at interval or ratio scales. The value of Moran’s I is calculated based on the deviation from the mean of two neighboring values [13]. The mathematical formula is as follows:(1)
where N is the sample size, is the mean of the variable, Xi is the value of the variable at a particular location i, Xj is the variable value at location j, and Wij is a spatial weight indexing the location of i relative to j. The value of this statistic is scored between −1 and 1. A score close to 1 represents positive autocorrelation and townships that may be hot spots. A score near −1 shows negative autocorrelation, indicating that the values of neighboring areas are opposite that of the township being examined. The significance of Moran’s I is evaluated by using a Z score and p-value generated by random permutation. The null hypothesis states that there is no spatial autocorrelation for the variable within the geographic area.
Spatial neighbors can be defined by a spatial weight matrix that is created in accordance with the neighbor definition chosen. We first calculated the mean distances between population centers of townships. Townships with shorter distances between their population centers were defined as neighbors. The geospatial relationships between pairs of the 349 townships were stored in a 349×349 matrix. The weight of each cell was the inverse of the distance between the two neighbors.
Spatial Lag Model
We inspected the residuals from the ordinary least squares (OLS) regression model to identify the spatial dependency of the residuals. If spatial dependency exists, it violates the assumption that the error terms of individual observations are independent of each other in the OLS regression; therefore a model that considers spatial autocorrelation is necessary [14].
We used a spatially lagged y model [15] that incorporated a spatially lagged dependent variable (y) on the right side of the regression equation. This regression model in matrix notation is represented as follows:(2)where ρ is the spatial lag parameter, W is the spatial weight matrix, X is a matrix of explanatory variables with an associated vector of regression coefficients β, and ε is a vector of normally distributed, random error terms. When the parameter associated with spatial lag (ρ) was positive, it indicated that, for townships where TB incidence was high, their neighbors also had a high TB incidence.
Because y is recruited in both sides of the regression equation, spatial dynamics creates a feedback effect between townships, in which a township’s level of TB incidence has an effect on its neighbors’, and the neighbors’ neighbors are also affected, throughout all connected townships [14]. This phenomenon leads to a chain reaction that finally returns to influence the initial TB incidence via the spatially lagged y term. In equilibrium, the expected value for y is calculated as follows:(3)
The spatial multiplier, (I-ρW)−1, shows how much the change in independent variable x in one township “spills over” onto other surrounding townships. This “spillover” then affects y through the effect of its spatial lag [16].
We used maps and a histogram to illustrate the variability in the spatial spillover (diffusion) of each township. These figures present the spillover (diffusion) at equilibrium of TB incidence into the surrounding townships with one-unit changes in the explanatory variables.
We were also interested in determining whether a neighbors’ previous TB incidence could be associated with that township’s future TB incidence. The space-time model appears as yt = ρWyt−1+Xβ+ε, where we set yt as the TB incidence from 2006–2008 and yt−1 as the TB incidence from 2003–2005.
Statistical Analysis
The associations between socioeconomic variables and TB incidences were analyzed using a linear regression model. Natural logarithmic transformations were used for TB incidence to accommodate the assumption of normal distribution. Stepwise regression modeling was conducted using SAS version 9.2 (SAS Institute, Cary, North Carolina). Moran’s I statistic calculation, the permutation process, and the spatial regression analysis were performed using Geoda® version 0.9.5-I [17]. The spatial multiplier was calculated using R version 2.9.0.
Results
Spatial Distribution of TB Incidence
Figure 1 shows the TB incidences of the 349 townships from 2003–2008. The Moran’s I statistic for the TB incidences of the 349 townships was positive (0.37) and statistically significant, indicating the presence of spatial clustering of TB incidences.
Univariate Analysis and OLS Model
We performed linear regression to identify the socioeconomic variables associated with higher TB incidence. Univariate analyses were performed for each independent variable, and they showed that most socioeconomic variables were significant, except for lower middle education (EDU1), the percentage of elderly (ELDER_P), and lower middle income (INCOME1) (Table 1). There was also a correlation between the socioeconomic variables (Table 2). We subsequently conducted stepwise multiple regression analysis, which showed that only the percentages of aborigines (ABOR_P), middle income (INCOME2), and high income (INCOME3) were independent factors (Table 3). The variance inflation factor of these variables remained below 2, excluding multicollinearity.
We again used Moran’s I statistic to test if there was still spatial autocorrelation for the residuals of OLS regression. The Moran’s I statistic was 0.18, indicating that the independent variables in the OLS model did not account for all spatial dependence in the outcome variable. These results confirmed the need to conduct spatial regression.
Spatial Lag Model
Spatial lag regression was conducted using the distance between population centers of polygons as the spatial weight. These results are shown in Table 3. Both the percentage of aborigines and high household income remain significant in the spatial lag model. A high percentage of aborigines was associated with higher TB incidence, while a high average household income was associated with lower TB incidence. These associations became smaller in the spatial lag model. Middle income (INCOME2) was significant in the OLS model but not in the spatial lag model. The ρ coefficient for the spatial parameter was significant and positive (0.43, p<0.001), which implies a positive correlation between the TB incidences of neighboring townships.
The log likelihood and Akaike’s information criterion (AIC) showed that the spatial lag model had a better fit than the OLS model. The Moran’s I statistic for the residuals of the spatial lag model was 0.05, which was very close to 0. This demonstrated that the spatial parameter could eliminate the effect of spatial autocorrelation in the regression model.
Spatial Multiplier
The spatial multiplier for the spatial lag model was calculated for each township and presented in Figure 2. This multiplier represented the interdependence of TB incidence for adjacent townships and had a minimum value of 1.06, which indicated that the independent variables for every township had a certain degree of spillover. The average value was 1.74 and the standard deviation was 0.15. Townships with high spatial multiplier values were mainly located in the northern, western and eastern parts of Taiwan.
Spatial-Time Lag Model
We further consider a spatial-time lag model. The model appeared as yt = Xβ+ρWyt−1+ε, where we set yt as the log-transformed TB incidence from 2006–2008 (under national DOT programs) and yt−1 as the log-transformed TB incidence from 2003–2005 (before national DOT programs). Using yt as the dependent variable, univariate analyses were performed; these analyses are presented in Table 1. The result was similar to using the log transformed TB incidence from 2003–2008 as the dependent variable. The term Wyt−1 was then calculated along with the percentages of aborigines, middle income, and high income and put into the model as independent variables. The spatial-time lag parameter, percentage of aborigines, and high average household income remained significant (Table 3). Therefore, the neighboring townships’ TB incidences from 2003–2005 were associated with a township’s TB incidence from 2006–2008. Moran’s I for the residuals of the spatial-time lag model (0.13) were much smaller than that of the log transformed TB incidence from 2006–2008 (0.29). Thus, the spatial-time lag parameter could partially eliminate the effect of spatial autocorrelation.
Discussion
Our geospatial analysis of the countrywide TB data for Taiwan indicated that the TB incidence in a township was significantly affected by the TB incidence in neighboring townships, which implies that recent transmission plays a significant role in TB endemic in Taiwan. Therefore, in addition to the current focus on improving the cure rate under DOT programs, more resource need to be allocated to active case finding in order to break the chain of transmission.
Using spatial regression modeling, we demonstrated that there exists a spatial dependency of township-level TB incidences in Taiwan, after adjusting for the effects of confounding socioeconomic variables, including the percentages of aboriginals and average household income. Furthermore, when we considered the temporality of the infectious processes, the spatial-time lag model indicated that a town’s TB incidences from 2006–2008 were affected by their neighbors’ TB incidences from 2003–2005, as would be expected from the cumulative effects of local TB transmission with contagious diffusion among the community.
The geospatial findings in the present study are consistent with molecular epidemiologic findings [5]–[9]. A large nosocomial TB outbreak in 2003 involving 66 health care workers at a district hospital in Taipei was traced to an index case who was hospitalized in February 2002 by matched DNA fingerprints [5]. Mycobacterial interspersed repetitive-unit-variable-number tandem-repeat (MIRU-VNTR) typing and spoligotyping was performed on TB isolates from 365 patients treated at a hospital in Taipei from 2002–2004; these results showed that 236 (65%) were clustered by genotype [7]. Another study in Hualien County in eastern Taiwan showed that 45 (62%) of 73 multidrug-resistant TB isolates were clustered [9]. These clustering rates were significantly higher than those reported in San Francisco (40%) [18], New York (37.5%) [19], the Netherlands (46%) [20], and Denmark (57%) [21], but lower than those reported from Malawi (72%) [22] and South Africa (67% [23], 72% [24]). Molecular genotyping further revealed that at least 51% of recurrent TB cases in Taiwan were caused by re-infection by a different strain, rather than by relapse or re-activation [25]. In keeping with previous molecular epidemiologic findings, our geospatial analyses provide necessary, complementary evidence on the significant role of recent TB transmission in Taiwan.
Our analysis showed that the percentage of aborigines is an independent risk factor for higher TB incidence after adjusting for the effects of spatial dependency and household income. This finding was in agreement with previous studies of TB incidence in Taiwan [26]–[29]. It has been shown that aboriginal areas have a TB incidence that is 3–5 times higher than non-aboriginal areas [26] and that the socioeconomic and health statuses of people living in aboriginal areas were generally lower than the national average [27], [28]. Poor compliance with anti-TB treatment might lengthen infectious period, thus increasing transmission [29].
Consistent with previous observations that TB is a disease of the deprived and the poor [30]–[32], our analysis found that a high average household income (the highest quartile) is an independent factor for lower TB incidence. TB is related to poverty in a number of ways, including higher contact rates due to crowded and poorly ventilated environments, reduced immunity status, and decreased odds of receiving proper treatment [30]. From Table 2, we can see that high average household income (INCOME3) was significantly correlated with all other variables, which may be the reason that most other socioeconomic variables were insignificant in the stepwise multiple regression analysis.
HIV infection weakens the immunity of patients and increases the risk of rapid progression to active TB disease after infection [33]. High HIV prevalence in the population may increase TB incidence [33]–[35]. In univariate analysis, we found that the cumulative HIV incidence from 1984–2002 was a significant risk factor for higher TB incidences in 2003–2005, as well as 2003–2008 (Table 1). Nevertheless, HIV infection status did not remain an independent factor in the multiple regression model. One probable reason for this change is the low HIV prevalence in Taiwan: there were only 4,145 adult HIV cases at the end of 2002 out of a population of 23 million. In addition, the positive correlation between HIV and high average income (Table 2) could mask the potential effect of HIV infection on TB incidence.
The resolution of the geospatial analysis in the present study was limited to the township level because further details on the residential addresses of TB patients were kept confidential by the Notifiable Infectious Disease Statistics System. Therefore, we were unable to use spatial point analysis methods to identify localized spatial clustering of TB cases. Another limitation of this study is the lack of data on the molecular genotype of clinical isolates and the host factors of individual persons, as well as the social network data, which restricts our inferences to the ecological level. The last limitation is that, if a spatially autocorrelated determinant of reactivated latent TB cases has been overlooked, our conclusions could be incorrect. We do take into consideration a range of important socioeconomic factors, but it is still possible that an important variable is missing. Our findings justify further large-scale genotyping-geospatial correlation studies to provide more insight on TB epidemiology in Taiwan.
In conclusion, our results add to the evidence that recent transmission plays a significant role in TB incidence in Taiwan, as well as highlighting the importance of taking a geospatial perspective in TB epidemiology.
Author Contributions
Conceived and designed the experiments: THW CTF. Analyzed the data: ICN. Contributed reagents/materials/analysis tools: THW CTF. Wrote the paper: ICN THW CTF. Interpreted data and revised the manuscript: JYW.
References
- 1. Horsburgh CR (2004) Priorities for the treatment of latent tuberculosis infection in the United States. N Engl J Med 350: 2060–2067.
- 2. Moonan PK, Ghosh S, Oeltmann JE, Kammerer JS, Cowan LS, et al. (2012) Using genotyping and geospatial scanning to estimate recent Mycobacterium tuberculosis transmission, United States. Emerg Infect Dis 18: 458–465.
- 3.
Centers for Disease Control T Notifiable Infectious Diseases Statistical System - Tuberculosis. Available: http://nidss.cdc.gov.tw/SingleDisease.aspx?dt=3&dc=1&disease=010&d=1&i=0&s=determined_cnt&RK=Y. Accessed 2009 January 3.
- 4. Yu MC, Bai KJ, Chang JH, Lee CN (2006) Age transition of tuberculosis patients in Taiwan, 1957–2001. J Formos Med Assoc 105: 25–30.
- 5. Huang WL, Jou R, Yeh PF, Huang A, Outbreak Investigation T (2007) Laboratory investigation of a nosocomial transmission of tuberculosis at a district general hospital. J Formos Med Assoc 106: 520–527.
- 6. Yang CY, Jou R, Chuang PC, Chang JT, Lee JJ, et al. (2007) Transmission of Mycobacterium tuberculosis in a family proved by genotyping. J Formos Med Assoc 106: 808–814.
- 7. Dou HY, Tseng FC, Lin CW, Chang JR, Sun JR, et al. (2008) Molecular epidemiology and evolutionary genetics of Mycobacterium tuberculosis in Taipei. BMC Infect Dis 8: 170.
- 8. Chen TC, Lu PL, Yang CJ, Lin WR, Lin CY, et al. (2010) Management of a nosocomial outbreak of Mycobacterium tuberculosis Beijing/W genotype in Taiwan: an emphasis on case tracing with high-resolution computed tomography. Jpn J Infect Dis 63: 199–203.
- 9. Hsu AH, Lin CB, Lee YS, Chiang CY, Chen LK, et al. (2010) Molecular epidemiology of multidrug-resistant Mycobacterium tuberculosis in Eastern Taiwan. Int J Tuberc Lung Dis 14: 924–926.
- 10. Horsburgh CR, O’Donnell M, Chamblee S, Moreland JL, Johnson J, et al. (2010) Revisiting rates of reactivation tuberculosis: a population-based approach. Am J Respir Crit Care Med 182: 420–425.
- 11. Johnson IL, Thomson M, Manfreda J, Hershfield ES (1985) Risk factors for reactivation of tuberculosis in Manitoba. CMAJ 133: 1221–1224.
- 12.
Centers for Disease Control T Notifiable Infectious Diseases Statistical System - HIV. Available: http://nidss.cdc.gov.tw/SingleDisease.aspx?Pt=s&Dc=1&Dt=3&disease=044. Accessed 2009 January 3.
- 13. Moran PA (1950) Notes on continuous stochastic phenomena. Biometrika 37: 17–23.
- 14.
Ward MD, Gleditsch KS (2008) Spatial Regression Models. Thousand Oaks, CA.: Sage.
- 15.
Anselin L (1988) Spatial econometrics : methods and models. Dordrecht; Boston: Kluwer Academic Publishers.
- 16. Anselin L (2003) Spatial externalities, spatial multipliers, and spatial econometrics. Int Regional Science Rev 26: 153–166.
- 17.
Anselin L. GeoDa 0.9.5-i. GeoDa Center for Geospatial Analysis and Computation. https://geodacenter.asu.edu/.
- 18. Small PM, Hopewell PC, Singh SP, Paz A, Parsonnet J, et al. (1994) The epidemiology of tuberculosis in San Francisco. A population-based study using conventional and molecular methods. N Engl J Med 330: 1703–1709.
- 19. Alland D, Kalkut GE, Moss AR, McAdam RA, Hahn JA, et al. (1994) Transmission of tuberculosis in New York City. An analysis by DNA fingerprinting and conventional epidemiologic methods. N Engl J Med 330: 1710–1716.
- 20. van Soolingen D, Borgdorff MW, de Haas PE, Sebek MM, Veen J, et al. (1999) Molecular epidemiology of tuberculosis in the Netherlands: a nationwide study from 1993 through 1997. J Infect Dis 180: 726–736.
- 21. Bauer J, Yang Z, Poulsen S, Andersen AB (1998) Results from 5 years of nationwide DNA fingerprinting of Mycobacterium tuberculosis complex isolates in a country with a low incidence of M. tuberculosis infection. J Clin Microbiol 36: 305–308.
- 22. Glynn JR, Crampin AC, Yates MD, Traore H, Mwaungulu FD, et al. (2005) The importance of recent infection with Mycobacterium tuberculosis in an area with high HIV prevalence: a long-term molecular epidemiological study in Northern Malawi. J Infect Dis 192: 480–487.
- 23. Godfrey-Faussett P, Sonnenberg P, Shearer SC, Bruce MC, Mee C, et al. (2000) Tuberculosis control and molecular epidemiology in a South African gold-mining community. Lancet 356: 1066–1071.
- 24. Verver S, Warren RM, Munch Z, Vynnycky E, van Helden PD, et al. (2004) Transmission of tuberculosis in a high incidence urban community in South Africa. Int J Epidemiol 33: 351–357.
- 25. Wang JY, Lee LN, Lai HC, Hsu HL, Liaw YS, et al. (2007) Prediction of the tuberculosis reinfection proportion from the local incidence. J Infect Dis 196: 281–288.
- 26. Yu MC, Bai KJ, Chang JH, Lee CN (2004) Tuberculosis incidence and mortality in aboriginal areas of Taiwan, 1997–2001. J Formos Med Assoc 103: 817–823.
- 27. Ko YC, Hsieh SF (1994) Leading causes of death in the aborigines in Taiwan. Kaohsiung J Med Sci 10: 352–366.
- 28. Ko YC, Liu BH, Hsieh SF (1994) Issues on aboriginal health in Taiwan. Kaohsiung J Med Sci 10: 337–351.
- 29. Chang YM, Tsai BY, Wu YC, Yang SY, Chen CH (2006) Risk of Mycobacterium tuberculosis transmission in an aboriginal village, Taiwan. Southeast Asian J Trop Med Public Health 37 Suppl 3161–164.
- 30. Waaler HT (2002) Tuberculosis and poverty. Int J Tuberc Lung Dis 6: 745–746.
- 31. Figueroa-Munoz JI, Ramon-Pardo P (2008) Tuberculosis control in vulnerable groups. Bull World Health Organ 86: 733–735.
- 32. Bishai WR, Graham NM, Harrington S, Pope DS, Hooper N, et al. (1998) Molecular and geographic patterns of tuberculosis transmission after 15 years of directly observed therapy. JAMA 280: 1679–1684.
- 33. Corbett EL, Watt CJ, Walker N, Maher D, Williams BG, et al. (2003) The growing burden of tuberculosis: global trends and interactions with the HIV epidemic. Arch Intern Med 163: 1009–1021.
- 34. Houben RM, Crampin AC, Mallard K, Mwaungulu JN, Yates MD, et al. (2009) HIV and the risk of tuberculosis due to recent transmission over 12 years in Karonga District, Malawi. Trans R Soc Trop Med Hyg 103: 1187–1189.
- 35. Nava-Aguilera E, Andersson N, Harris E, Mitchell S, Hamel C, et al. (2009) Risk factors associated with recent transmission of tuberculosis: systematic review and meta-analysis. Int J Tuberc Lung Dis 13: 17–26.