Skip to main content

Community transmission of multidrug-resistant tuberculosis is associated with activity space overlap in Lima, Peru



Transmission of multidrug-resistant tuberculosis (MDRTB) requires spatial proximity between infectious cases and susceptible persons. We assess activity space overlap among MDRTB cases and community controls to identify potential areas of transmission.


We enrolled 35 MDRTB cases and 64 TB-free community controls in Lima, Peru. Cases were whole genome sequenced and strain clustering was used as a proxy for transmission. GPS data were gathered from participants over seven days. Kernel density estimation methods were used to construct activity spaces from GPS locations and the utilization distribution overlap index (UDOI) was used to quantify activity space overlap.


Activity spaces of controls (median = 35.6 km2, IQR = 25.1–54) were larger than cases (median = 21.3 km2, IQR = 17.9–48.6) (P = 0.02). Activity space overlap was greatest among genetically clustered cases (mean UDOI = 0.63, sd = 0.67) and lowest between cases and controls (mean UDOI = 0.13, sd = 0.28). UDOI was positively associated with genetic similarity of MDRTB strains between case pairs (P < 0.001). The odds of two cases being genetically clustered increased by 22% per 0.10 increase in UDOI (OR = 1.22, CI = 1.09–1.36, P < 0.001).


Activity space overlap is associated with MDRTB clustering. MDRTB transmission may be occurring in small, overlapping activity spaces in community settings. GPS studies may be useful in identifying new areas of MDRTB transmission.

Peer Review reports


Tuberculosis (TB) is commonly transmitted outside of the home in community-based settings during social contact between infectious cases and susceptible community members [1, 2]. Prior studies have used questionnaires to identify epidemiologic links between unrelated TB cases and have found spatial links between cases in community settings [3,4,5,6]. Identifying spatial areas of transmission is important for contact tracing and infection control, as epidemiologic links between TB cases are often unclear [7, 8]. Moreover, targeting infection control in geographic hotspots of TB transmission may reduce overall levels of community transmission [9].

Activity spaces are used in epidemiologic studies to represent geographic spaces wherein people spend their time during regular daily activities [10,11,12]. Past studies of human activity spaces have relied predominantly on place-tracing questionnaires to delineate activity spaces [3,4,5, 13]. More recent approaches now leverage global positioning systems (GPS) [10, 14, 15]. The widespread availability of GPS technology has facilitated fine scale study of human movement patterns and are less prone to recall and measurement errors common in retrospective place-tracing interviews [16,17,18,19].

The goal of this study was to compare the activity spaces of multidrug resistant tuberculosis (MDRTB) cases and healthy (TB-free) community controls, identify areas of activity space overlap among clustered cases to identify areas of potential transmission, and to quantify the association between activity space overlap and genetic clustering of MDR Mycobacterium tuberculosis (Mtb) strains.


Study design and setting

Study participants were recruited from the areas of Callao and Lima Sur located to the north and south of Lima, Peru, respectively. These two regions report the greatest proportion of incident MDRTB cases in Peru [20]. Callao has an area of 147 km2 with nearly all one million residents living in urban areas. Lima Sur encompasses 11 districts with a total area of 852 km2 and 1.5 million residents.

Between February 2016 and May 2017, patients were recruited from a completed parent study that enrolled a household-based cohort of MDRTB patients between 2010 and 2013 [21]. Sputum samples from cases were taken at diagnosis and processed on liquid microscopic observation drug susceptibility assays (MODS) and solid Ogawa media. Aliquots of positive sputum samples were reserved for DNA extraction and genotyping via whole genome sequencing [22]. The single nucleotide polymorphism (SNP) calling analysis was performed on an Illumina HiSeq2000 with paired-end reads of length of 100 bp [23]. A pairwise matrix of MDRTB cases and the number of single nucleotide polymorphism (SNP) differences in their Mtb strains was assembled; cases with Mtb strains within ≤5 SNP differences were considered genetically clustered [24]. This threshold was our working definition for MDRTB transmission. Exclusion criteria included genetically clustered pairs from the same household. Where a case was genetically clustered to multiple cases from the same household, only one member of that household was enrolled.

We enrolled community controls, who verbally confirmed that they had never received TB treatment or diagnosis, as a comparison sample. Controls included community health workers and nurses that worked at community health posts serving case neighborhoods. Additional controls were referred by community health workers and sourced from churches, restaurants, communal kitchens and education centers located in case neighborhoods. Controls were frequency matched to cases on age (± 5 years), sex (male/female) and study region (Callao/Lima Sur) to ensure comparability across these variables.

Informed consent was obtained from participants prior to data collection. The study protocol, consent forms and data collection instruments were reviewed and approved by the Institutional Committee of Ethics for Humans at La Universidad Peruana Cayetano Heredia.

Data sources and measurements

Questionnaires were used to collect demographic information from participants during face-to-face interviews. We used Qstarz BT-Q1000XT (Qstarz International, Taipei, Taiwan) GPS loggers to gather data on participant’s movements over seven days of observation. The units were configured to log participant’s locations (i.e., geocoordinates) every minute. Consenting participants were given a GPS logger and instructed to keep the logger powered on, carried with them at all times, and recharged nightly. Study nurses called participants every other day during the 7-day data collection period to remind them to carry GPS loggers.

Constructing activity spaces

Spatial ecologists have developed a suite of methods to study movement patterns of wildlife using GPS technology to delineate areas of regular space use called ‘home ranges’ and are analogous to activity spaces in human research [25]. These methods account for non-uniform space use by representing home ranges (i.e., activity spaces) as spatial probability density functions of space use called utilization distributions (UD) [26]. Instead of assuming space-use is uniform across an activity space, UDs highlight areas of concentrated activity with probability contours and are better representations of space use. In this study, we use home ranges to represent participant activity spaces.

Kernel density estimation (KDE) was used to construct participant activity spaces [25, 27]. We used iterative visualizations of GPS kernel densities generated with a Gaussian kernel and bandwidths between 100 m and 1200 m to identify UDs that provided adequate smoothing of raw GPS locations while highlighting distinct “peaks” of areas where locations were concentrated. The final chosen bandwidth was 950 m.

Home range (i.e., activity space) sizes were estimated at 50, 95 and 99% contours of each participant’s UD. These percentages correspond to the smallest home range area encompassing 50, 95 and 99% of a participant’s GPS locations. The 95% contour is the standard used in home range studies, while the 50% contour is considered the “core area” of activity [28, 29]. The 99% contour is the most inclusive contour, containing areas of sparse activity.

Measuring spatial overlap

After estimating the UDs to represent each participant’s activity space, we calculated the utilization distribution overlap index (UDOI) to quantify the amount of activity space spatial overlap between participants [26]. The UDOI is estimated as the cumulative sum of the cell-by-cell product of two participant UDs multiplied by the intersecting area (i.e., product of two UDs) [26, 28]. The UDOI of two participants is high when their GPS locations are concentrated within the same space [26]. The UDOI ranges from 0 (no spatial overlap) to 1 (complete spatial overlap) and can take on values > 1 if the UDs are non-uniformly distributed and have a high degree of overlap. We estimated the UDOI’s at each home range contour level (50, 95, and 99%) to examine the magnitude of association between MDRTB transmission and activity space overlap. We created UDs for each participant and estimated the UDOIs for all pairs of participants using the ‘adehabitatHR’ package [30] in R (version 3.6.1, The R Foundation).

Statistical analyses

T-tests, Wilcoxon rank sum tests (when appropriate), and chi-squared tests were used to compare cases and controls by demographics, home range size and mean UDOI.

We compared the mean pairwise UDOI of cases and controls to determine the degree of spatial overlap within and between groups. The mean UDOI of case dyads (i.e., case-case pairs), control dyads (i.e., control-control pairs), and case and control dyads (i.e., case-control pairs) were evaluated. Bonferroni adjusted P-values were reported to account for multiple comparisons.

We used logistic regression to estimate the odds of being genetically clustered as a function of spatial overlap (among cases). Linear regression was used to assess the relationship between the UDOI of case pairs and degree of genetic strain similarity (i.e., SNP differences). SNP difference values were log (base 10) transformed and logit transformed UDOI values were used to meet normality assumptions for modelling.


Comparing cases and controls

A total of 99 participants were enrolled, including 35 MDRTB cases (35%) and 64 healthy community controls (65%). Only 35 participants of the original study could be contacted and consented, as the majority were no longer contactable, had moved house or had died. Sixteen (46%) cases were genetically clustered with matching Mycobacterium tuberculosis (Mtb) strains (within ≤5 SNPs). Frequency matching between cases and controls was achieved by region, age and sex (Table 1). The median age, gender and regional distribution of participants in both case and control groups were comparable, though not found to be statistically significant. Participants in both groups had similar levels of employment, but participants in the case group generally had a lower level of educational attainment compared to controls (31.4% vs 6.3% had primary education or less). The case group had a lower monthly income compared to the control group (77.1% vs 20.3% earned less than 1000 PEN/month). The differences in education and income were statistically supported. The sublineages of Mtb strains included in our study cohort were compared with those from the parent study population (Table S1). There was a similar distribution of sublineages represented between the two studies.

Table 1 Summary of study participants and comparison of multidrug resistant tuberculosis (MDRTB) cases and TB-free community controls enrolled in this study

A median of 8253 GPS locations (IQR 5634-9828) were collected per participant, representing about 138 h of data (IQR 94–164) per participant. A median of 7756 GPS locations (IQR 4813-9493), equivalent to 129 h (IQR 80–158), were collected for cases and a median of 8364 GPS locations (IQR 6766-9887), equivalent to 139 h (IQR 113–165), were collected for controls. There were, however, no statistically significant differences in the amount of GPS locations collected between cases and controls.

Home ranges of cases and controls

The aggregated 50 and 95% home ranges of all cases in Callao and Lima Sur are shown in Fig. 1.

Fig. 1

District maps of Lima, Peru illustrating home ranges in Callao and Lima Sur. a home ranges of all MDRTB cases in Callao and b) home ranges of all MDRTB cases in Lima Sur. The 50% home range contours are drawn with a thick solid border and 95% home range contours are drawn with a dashed line

In Callao, the 50% home range of cases were located on either side of Jorge Chávez International Airport in high traffic corridors and densely populated neighborhoods; another 50% core was located in a densely commercial area of Ventanilla. In Lima Sur, there was one 50% home range core located in the district of Villa Maria Del Triunfo. Case home ranges were significantly smaller than control home ranges at all contour levels (Table 1).

Spatial overlap between cases and controls

Table 2 summarizes the mean UDOI by dyad type (i.e., pairs of participant types) at the 95% home range contour. Mean UDOI was greatest among clustered case dyads (mean = 0.29, sd = 0.55) in all regions whereas mean UDOI was lowest among non-clustered cases (mean = 0.05, sd = 0.20) (P < 0.01). Mean UDOI was higher among controls (mean = 0.14, sd = 0.36) than among cases (mean = 0.06, sd = 0.22) (P < 0.001). Figure 2 provides three illustrative examples of low, medium and high UDOI values for three dyads of cases. There was a high proportion of non-overlapping dyads (between 35 to 47% of possible dyads that did not have overlapping activity spaces). As a sensitivity analysis, we analyzed the mean UDOI of dyad types with only overlapping dyads and found that mean UDOI increased across all dyad types, but there was no change in inferences (Table 2).

Table 2 Comparison of mean 95% utilization distribution overlap index (UDOI) values by dyad type and study regions (Callao and Lima Sur)
Fig. 2

Three examples of home ranges with differing 95% home range overlaps. Low (Panel a), moderate (Panel b) and high (Panel c) 95% home range overlaps measured by their utilization distribution overlap index (UDOI). The small circles are GPS relocation trails and the larger outlines are the 95% home ranges of each case. The blue and red coloring are used to differentiate the two cases. Map outline is of districts in Lima. Panel a shows two (non-clustered) cases with very low overlap (UDOI< 0.01); Panel b shows two (non-clustered) cases with moderate overlap (UDOI = 0.63) and Panel c shows two (clustered) cases with very high overlap (UDOI> 1.0)

Association between spatial overlap and genetic clustering

There was a statistically significant positive association between the UDOI of cases and genetic clustering (i.e. being part of a recent transmission chain). At the 95% contour, the odds of two cases being clustered increased by 22% per 0.10 increase in their shared UDOI (OR = 1.22, CI = 1.09–1.36, P < 0.001). The effect was larger at the 50% contour where the odds of two cases being clustered was increased by four-fold for every 0.10 increase in UDOI (OR = 4.25, CI = 1.85–9.73, P = 0.001) (Table 3).

Table 3 Association between genetic similarity of Mycobacterium strains of two case pairs (measured in number of different single nucleotide polymorphisms (SNP) and their utilization distribution overlap index (UDOI)

There was a statistically significant association between spatial overlap and Mtb genetic similarity (i.e., SNP differences) among cases (Fig. 3a and b). At the 95% home range contour, for every log increase in the number of SNP difference, there was a 2.6-unit reduction in the logit transformed UDOI (coef = − 2.6, CI = -3.8–1.4, P < 0.001) or approximately 0.0025 reduction in UDOI (Table 3). In other words, as genetic similarity of two Mtb strains increased, their level of home range overlap increased.

Fig. 3

Two-way scatter plots and linear regression predictions of 95% UDOI and SNP difference. a Two-way scatter plot of 95% UDOI and SNP difference between case pairs, and b) linear regression prediction of 95% UDOI and SNP difference with confidence interval. There is a statistically significant negative association between UDOI and SNP difference (Coef = − 2.59, P < 0.001). Note UDOI is logit transformed to reduce variable skew for regression modelling. Panel c and d exclude case pairs with no overlap

Since there was a large proportion of non-overlapping dyads that might have influenced the negative association between SNP difference and UDOI, we conducted a sensitivity analysis excluding non-overlapping dyads and found the association between Mtb genetic dissimilarity and UDOI of two cases attenuated but was still statistically significant (Fig. 3c and d). During post-hoc review, we found one clustered dyad of related family members, resulting in a very high UDOI that may be considered an outlier. After excluding this dyad from analysis, the mean UDOI of clustered cases was reduced, but remained greater than all other dyad pairs (Table S2). The association between the genetic similarity of two Mtb strains and UDOI was slightly attenuated, but remained statistically significant at all home range levels (Table S3).


In this study, we compared the activity spaces of MDRTB cases to healthy community controls. We found that cases had significantly smaller activity spaces than controls, that the activity space overlap was greatest among genetically clustered cases, and that there was a statistically significant association between activity space overlap and Mtb strain genetic similarity. These findings suggest that MDRTB contact, exposure and transmission may be occurring among cases in relatively small, overlapping activity spaces in community settings.

The demonstration of high overlap amongst genetically clustered cases and lower overlap between cases and controls in the dyad analysis suggests that spatial segregation between the groups may be occurring, as found in a previous study [31]. Spatial segregation of MDRTB cases may support the use of spatially targeted screening interventions to improve local control and indirectly reduce MDRTB prevalence [32]. The high spatial overlap among cases may explain why spatial clustering of MDRTB genotypes has been previously found in localized hotspots in Lima [2]. The association between Mtb genetic similarity and spatial overlap was observed at all home range contours but was particularly pronounced at the smaller 50% home range contour (median of 5km2), suggesting that transmission of MDRTB may be occurring at very local levels near case residences. Yang et al. found that genetic similarity of Mtb strains of MDRTB case pairs in China increased as their residential proximity increased [33]. A phenomenon of MDRTB spillover from a prison in Lima to the surrounding population was demonstrated by Warren et al. and proposes a mechanism whereby local transmission may be occurring in the community [34]. Moreover, as clustered cases tended to have substantially smaller activity spaces than non-clustered cases, this suggests that movement and spatial overlap in small geographic hotspots may be driving MDRTB transmission in this population. Focusing infection control in those areas of high overlap may reduce MDRTB transmission and result in community-wide benefits as suggested by Dowdy et al. [9]. It is important to recognize, however, that drivers of transmission hotspots are likely to vary between rural and urban environments. Nelson et al. demonstrated that transmission of extensively drug resistant TB (XDRTB) was likely occurring far away from cases’ residences in rural South Africa. The authors suggest that cyclical rural-to-urban migration for work was an important determinant of transmission (Nelson et al. 2018). Shift of XDRTB transmission towards workplaces in urban centers in Durban was also highlighted by Peterson et al. [13].

To date, studies investigating transmission sites for TB have relied on questionnaires about frequented locations [6], which generally underestimate spatial mobility and are subject to information bias [16, 18]. This study uses GPS tracking to obtain objective spatial information on participants and does not rely on participants’ recall. Activity spaces of individuals, calculated using GPS logging data, were shown to be larger than those derived from geotagging venues reported in questionnaires [35]. Moreover, place-tracing questionnaires often focus only on community venues and do not take into account routes travelled between them [35]. This suggests that GPS methods are superior in acquiring a greater amount of spatial information. This methodology is easily reproducible and demonstrates the utility of GPS tracking in combination with whole genome sequencing to identify potential transmission sites.

There are several limitations to our study that should be noted. Firstly, it is likely that selection bias was introduced through our non-random selection of controls. As this was an exploratory study, our sample of controls were convenience-based and were often health workers or their family and friends. As a result, they generally had a higher level of education and income and were not necessarily representative of the general population. Moreover, given the higher socioeconomic status (SES) of controls, it is likely that controls did not live in the same neighborhoods nor frequent the same shops and venues that cases did, which may have resulted in the observed low UDOI between case and control dyads (but the relatively higher UDOI among control-dyads). We attempted to address this source of bias through matching controls by age, sex and study region, but this was not sufficient to control for confounding caused by SES. Our healthy controls were also recruited on the basis that they had never been treated or diagnosed with previous TB. However, this was not confirmed by medical records and may represent another potential source of bias. Our analysis of genetic clusters relied on a SNP difference threshold to determine which cases were genetically clustered, representative of recent transmission [24]. While a small SNP difference between two M. tuberculosis strains is generally regarded as evidence of transmission, the appropriate clustering threshold depends on the environment and setting [36]. While GPS monitoring is considered to be more precise than structured interviews at identifying activity locations, GPS locations alone do not provide context for locations (i.e., types of places or reasons for visiting areas) that interviews could elicit; combining the two forms of data collection is preferable [15, 16]. On average, study participants provided slightly less than seven days of GPS location data, so these movement patterns might not be representative of typical activity. However, prior studies have found that human movement patterns tend to be regular and stable, particularly in urban settings where routines are structured and people tend to spend significant amounts of time in few, regularly visited locations [18, 25, 37].

The small sample size was another limitation to the study. This was mainly due to the large numbers of potential participants who had moved or had unfortunately died. The findings of this study can therefore only be regarded as exploratory in nature.

A larger prospective study with a bigger sample size would be useful to confirm our findings and determine whether the differences between cases and controls, in terms of activity space size and overlap likelihood at all activity space contour sizes, are statistically meaningful. Potential sources of bias could be reduced through the recruitment of confirmed TB-free controls from cases’ neighbourhoods. A detailed investigation into community venues within overlapping activity spaces of genetically-clustered participants is essential to isolate specific areas where transmission is occurring. Additionally, follow-up questionnaires conducted alongside GPS telemetry would be useful to characterise the context of visits to community venues. Data collection could be made simpler and more effective by using alternative sources of GPS data, such as Google maps location history on participants’ smartphones. Data on the movement patterns of TB patients during the period of transmission itself may provide greater insight into specific locations where transmission may have occurred.


In Lima, Peru, activity space overlap is associated with genetic clustering of MDRTB cases and case activity spaces are relatively small. Case finding activities should focus on areas within 5 km of case residences which is the core area of movement for cases and where community transmission may be most likely.

Availability of data and materials

The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.



Global positioning systems


Kernel density estimation


Multidrug resistant


Microscopic observation drug susceptibility assays


Mycobacterium tuberculosis


Single nucleotide polymorphism




Utilization distribution


Utilization distribution overlap index


Extensively drug resistant


  1. 1.

    Yates TA, Khan PY, Knight GM, Taylor JG, McHugh TD, Lipman M, et al. The transmission of Mycobacterium tuberculosis in high burden settings. Lancet Infect Dis. 2016;16:227–38.

    Article  PubMed  Google Scholar 

  2. 2.

    Zelner JL, Murray MB, Becerra MC, Galea J, Lecca L, Calderon R, et al. Identifying hotspots of multidrug-resistant tuberculosis transmission using spatial and molecular genetic data. J Infect Dis. 2015;213:287–94.

    Article  PubMed  PubMed Central  Google Scholar 

  3. 3.

    Chamie G, Wandera B, Marquez C, Kato-Maeda M, Kamya MR, Havlir DV, et al. Identifying locations of recent TB transmission in rural Uganda: a multidisciplinary approach. Trop Med Int Heal. 2015;20:537–45.

    Article  Google Scholar 

  4. 4.

    Chamie G, Kato-Maeda M, Emperador DM, Wandera B, Mugagga O, Crandall J, et al. Spatial overlap links seemingly unconnected genotype-matched TB cases in rural Uganda. PLoS One. 2018;13:e0192666.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  5. 5.

    Izumi K, Ohkado A, Uchimura K, Murase Y, Tatsumi Y, Kayebeta A, et al. Detection of tuberculosis infection hotspots using activity spaces based spatial approach in an urban Tokyo, from 2003 to 2011. PLoS One. 2015;10:e0138831.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  6. 6.

    Shaweno D, Karmakar M, Alene KA, Ragonnet R, Clements AC, Trauer JM, et al. Methods used in the spatial analysis of tuberculosis epidemiology: a systematic review. BMC Med. 2018;16:193.

    Article  PubMed  PubMed Central  Google Scholar 

  7. 7.

    Gardy JL, Johnston JC, Sui SJH, Cook VJ, Shah L, Brodkin E, et al. Whole-genome sequencing and social-network analysis of a tuberculosis outbreak. N Engl J Med. 2011;364:730–9.

    CAS  Article  PubMed  Google Scholar 

  8. 8.

    Klovdahl AS, Graviss EA, Yaganehdoost A, Ross MW, Wanger A, Adams GJ, et al. Networks and tuberculosis: an undetected community outbreak involving public places. Soc Sci Med. 2001;52:681–94

    CAS  Article  Google Scholar 

  9. 9.

    Dowdy DW, Golub JE, Chaisson RE, Saraceni V. Heterogeneity in tuberculosis transmission and the role of geographic hotspots in propagating epidemics. Proc Natl Acad Sci U S A. 2012;109:9557–62.

    Article  PubMed  PubMed Central  Google Scholar 

  10. 10.

    Hirsch JA, Winters M, Clarke P, McKay H. Generating GPS activity spaces that shed light upon the mobility habits of older adults: a descriptive analysis. Int J Health Geogr. 2014;13:51.

    Article  PubMed  PubMed Central  Google Scholar 

  11. 11.

    Sherman JE, Spencer J, Preisser JS, Gesler WM, Arcury TA. A suite of methods for representing activity space in a healthcare accessibility study. Int J Health Geogr. 2005;4:24.

    Article  PubMed  PubMed Central  Google Scholar 

  12. 12.

    Worrell MC, Kramer M, Yamin A, Ray SM, Goswami ND. Use of activity space in a tuberculosis outbreak: bringing homeless persons into spatial analyses. Open Forum Infect Dis. 2017;4.

  13. 13.

    Peterson ML, Gandhi NR, Clennon J, Nelson KN, Morris N, Ismail N, et al. Extensively drug-resistant tuberculosis “hotspots” and sociodemographic associations in Durban, South Africa. Int J Tuberc Lung Dis. 2019;23:720–7.

    CAS  Article  Google Scholar 

  14. 14.

    Lee NC, Voss C, Frazer AD, Hirsch JA, McKay HA, Winters M. Does activity space size influence physical activity levels of adolescents?—a GPS study of an urban environment. Prev Med Reports. 2016;3:75–8.

    Article  Google Scholar 

  15. 15.

    Shareck M, Kestens Y, Gauvin L. Examining the spatial congruence between data obtained with a novel activity location questionnaire, continuous GPS tracking, and prompted recall surveys. Int J Health Geogr. 2013;12:40.

    Article  PubMed  PubMed Central  Google Scholar 

  16. 16.

    Paz-Soldan VA, Reiner RC, Morrison AC, Stoddard ST, Kitron U, Scott TW, et al. Strengths and weaknesses of global positioning system (GPS) data-loggers and semi-structured interviews for capturing fine-scale human mobility: findings from Iquitos. Peru PLoS Negl Trop Dis. 2014;8:e2888.

    Article  PubMed  Google Scholar 

  17. 17.

    Thierry B, Chaix B, Kestens Y. Detecting activity locations from raw GPS data: a novel kernel-based algorithm. Int J Health Geogr. 2013;12:14.

    Article  PubMed  PubMed Central  Google Scholar 

  18. 18.

    Vazquez-Prokopec GM, Bisanzio D, Stoddard ST, Paz-Soldan V, Morrison AC, Elder JP, et al. Using GPS technology to quantify human mobility, Dynamic Contacts and Infectious Disease Dynamics in a Resource-Poor Urban Environment. PLoS One. 2013;8:e58802.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  19. 19.

    Vazquez-Prokopec GM, Stoddard ST, Paz-Soldan V, Morrison AC, Elder JP, Kochel TJ, et al. Usefulness of commercially available GPS data-loggers for tracking human movement and exposure to dengue virus. Int J Health Geogr. 2009;8:68.

    Article  PubMed  PubMed Central  Google Scholar 

  20. 20.

    Alarcon V, Alarcon E, Figueroa C, Mendoza-Ticona A. Tubersulosis en el Perú: situación epidemiológica, avances y desafíos para su control. Rev Peru Med Exp Salud Publica. 2017;34:299–310.

    Article  Google Scholar 

  21. 21.

    Grandjean L, Gilman RH, Martin L, Soto E, Castro B, Lopez S, et al. Transmission of multidrug-resistant and drug-susceptible tuberculosis within households: a prospective cohort study. PLoS Med. 2015;12:e1001843.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  22. 22.

    Grandjean L, Gilman RH, Iwamoto T, Köser CU, Coronel J, Zimic M, et al. Convergent evolution and topologically disruptive polymorphisms among multidrug-resistant tuberculosis in Peru. PLoS One. 2017;12:e0189838.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  23. 23.

    Grandjean L, Monteserin J, Gilman R, Pauschardt J, Bonilla C, Ritacco V, et al. Association between bacterial homoplastic variants and radiological pathology in tuberculosis Tuberculosis. Thorax. 2020.

  24. 24.

    Meehan CJ, Moris P, Kohl TA, Pečerska J, Akter S, Merker M, et al. The relationship between transmission time and clustering methods in Mycobacterium tuberculosis epidemiology. EBioMedicine. 2018;37:410–6.

    Article  PubMed  PubMed Central  Google Scholar 

  25. 25.

    Powell RA, Mitchell MS. What is a home range? J Mammal. 2012;93:948–58.

    Article  Google Scholar 

  26. 26.

    Fieberg J, Kochanny CO. Quantifying home-range overlap: the importance of the utilization distribution. J Wildl Manag. 2005;69:1346–59.

    Article  Google Scholar 

  27. 27.

    Worton BJ. Kernel methods for estimating the utilization distribution in home-range studies. Ecology. 1989;70:164–8.

    Article  Google Scholar 

  28. 28.

    Clapp JG, Beck JL. Evaluating distributional shifts in home range estimates. Ecol Evol. 2015;5:3869–78.

    Article  PubMed  PubMed Central  Google Scholar 

  29. 29.

    Walter WD, Onorato DP, Fischer JW. Is there a single best estimator?Selection of home range estimators using area-under-the-curve. Mov Ecol. 2015;3:10.

    Article  PubMed  PubMed Central  Google Scholar 

  30. 30.

    Calenge C. Home range estimation in R: the adehabitatHR package. 2015.

    Google Scholar 

  31. 31.

    Bui DP, Oren E, Roe DJ, Brown HE, Harris RB, Knight GM, et al. A Case Control Study to Identify Community Venues Associated with Genetically Clustered Multidrug-Resistant Tuberculosis Disease in Lima, Peru. Clin Infect Dis. 2018:ciy746.

  32. 32.

    Cudahy PGT, Andrews JR, Bilinski A, Dowdy DW, Mathema B, Menzies NA, et al. Spatially targeted screening to reduce tuberculosis transmission in high-incidence settings. Lancet Infect Dis. 2019;19:e89–95.

    Article  PubMed  Google Scholar 

  33. 33.

    Yang C, Lu L, Warren JL, Wu J, Jiang Q, Zuo T, et al. Internal migration and transmission dynamics of tuberculosis in Shanghai, China: an epidemiological, spatial, genomic analysis. Lancet Infect Dis. 2018;18:788–95.

    Article  PubMed  PubMed Central  Google Scholar 

  34. 34.

    Warren JL, Grandjean L, Moore DAJ, Lithgow A, Coronel J, Sheen P, et al. Investigating spillover of multidrug-resistant tuberculosis from a prison: a spatial and molecular epidemiological analysis. BMC Med. 2018;16:122.

    Article  PubMed  PubMed Central  Google Scholar 

  35. 35.

    Kestens Y, Thierry B, Shareck M, Steinmetz-Wood M, Chaix B. Integrating activity spaces in health research: comparing the VERITAS activity space questionnaire with 7-day GPS tracking and prompted recall. Spat Spatiotemporal Epidemiol. 2018;25:1–9.

    Article  PubMed  Google Scholar 

  36. 36.

    Hatherell H-A, Colijn C, Stagg HR, Jackson C, Winter JR, Abubakar I. Interpreting whole genome sequencing for investigating tuberculosis transmission: a systematic review. BMC Med. 2016;14:21.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  37. 37.

    Siła-Nowicka K, Vandrol J, Oshan T, Long JA, Demšar U, Fotheringham AS. Analysis of human mobility patterns from GPS trajectories and contextual information. Int J Geogr Inf Sci. 2016;30:881–906.

    Article  Google Scholar 

Download references


We’d like to acknowledge the field study team’s support in conducting this research.


This study was supported by a Fulbright-Fogarty Global Health fellowship grant (R25TW009340–5101872) and a Wellcome Trust Early Post-doctoral fellowship grant (201470/Z/16/Z).

Author information




DB and LG designed the study with feedback from EO, HEB, and RBH. DB and LG worked on data collection and statistical analyses with guidance from EO, HEB, and GMK. All authors contributed to interpretation of the data. DB and LG drafted the manuscript; SSC, EO, HEB, RBH, and GMK provided substantive revisions. All authors have read and approved the manuscript.

Corresponding author

Correspondence to Louis Grandjean.

Ethics declarations

Ethics approval and consent to participate

Informed consent was obtained from participants prior to data collection. The study protocol, consent forms and data collection instruments were reviewed and approved by the Institutional Committee of Ethics for Humans at La Universidad Peruana Cayetano Heredia. All methods were performed in accordance with the relevant guidelines and regulations of this institute.

Consent for publication

Not applicable.

Competing interests

None declared.

Additional information

Publisher’s Note

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

Supplementary Information

Additional file 1: Table S1.

Distribution of sublineages in study and parent study. Fisher-exact P-value = 0.83.

Additional file 2: Table S2.

Comparison of mean 95% UDOI by dyad type, excluding clustered cases from the same household.

Additional file 3: Table S3.

Association between genetic similarity of Mycobacterium tuberculosis strains of two case pairs (measured in number of different single nucleotide polymorphisms (SNP)) and their utilization distribution overlap index (UDOI), excluding clustered cases that are household members.

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 The Creative Commons Public Domain Dedication waiver ( 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

Bui, D.P., Chandran, S.S., Oren, E. et al. Community transmission of multidrug-resistant tuberculosis is associated with activity space overlap in Lima, Peru. BMC Infect Dis 21, 275 (2021).

Download citation


  • Tuberculosis
  • Community transmission
  • GPS
  • Home range
  • Activity space