This article has Open Peer Review reports available.
Assessment of tuberculosis spatial hotspot areas in Antananarivo, Madagascar, by combining spatial analysis and genotyping
© The Author(s). 2017
Received: 24 March 2017
Accepted: 30 July 2017
Published: 14 August 2017
Tuberculosis (TB) remains a public health problem in Madagascar. A crucial element of TB control is the development of an easy and rapid method for the orientation of TB control strategies in the country. Our main objective was to develop a TB spatial hotspot identification method by combining spatial analysis and TB genotyping method in Antananarivo.
Sputa of new pulmonary TB cases from 20 TB diagnosis and treatment centers (DTCs) in Antananarivo were collected from August 2013 to May 2014 for culture. Mycobacterium tuberculosis complex (MTBC) clinical isolates were typed by spoligotyping on a Luminex® 200 platform. All TB patients were respectively localized according to their neighborhood residence and the spatial distribution of all pulmonary TB patients and patients with genotypic clustered isolates were scanned respectively by the Kulldorff spatial scanning method for identification of significant spatial clustering. Areas exhibiting spatial clustering of patients with genotypic clustered isolates were considered as hotspot TB areas for transmission.
Overall, 467 new cases were included in the study, and 394 spoligotypes were obtained (84.4%). New TB cases were distributed in 133 of the 192 Fokontany (administrative neighborhoods) of Antananarivo (1 to 15 clinical patients per Fokontany) and patients with genotypic clustered isolates were distributed in 127 of the 192 Fokontany (1 to 13 per Fokontany). A single spatial focal point of epidemics was detected when ignoring genotypic data (p = 0.039). One Fokontany of this focal point and three additional ones were detected to be spatially clustered when taking genotypes into account (p < 0.05). These four areas were declared potential TB transmission hotspots in Antananarivo and will be considered as priority targets for surveillance in the future.
This method, combining spatial analysis and TB genotyping will now be used for further focused clinical and epidemiological studies in Madagascar and will allow better TB control strategies by public health authorities.
KeywordsMycobacterium tuberculosis Geographic Information System Genotyping Spatial cluster Antananarivo
Tuberculosis (TB), caused by Mycobacterium tuberculosis complex (M. tuberculosis), remains one of the deadliest infectious diseases worldwide. In 2014, 9.6 million people contracted TB and 1.5 million died from the disease . The number of deaths due to TB slowly declined between 2000 and 2013 due to effective diagnosis and treatment, but remains unacceptably high .
Spoligotyping is a method investigating the diversity at a highly variable CRISPR locus evolving by deletion in M. tuberculosis complex . Spoligotyping has been widely used to classify M. tuberculosis clinical isolates by family and subfamily [4, 5] that were later found to be largely concordant with lineages as defined by an whole genome sequencing (WGS) approach, allowing single nucleotide polymorphisms (SNPs) to be identified . In addition and despite a relatively low discrimination level of spoligotyping [7, 8] it can be used for first genetic identification of patient’s clinical isolates and suggest or exclude recent transmission cases which should be investigated further using more discriminatory methods [7, 9, 10]. Previous studies demonstrated that spatial clustering of TB data when associated to genetic clustering of TB cases more easily allows to focus on adequate settings to distinguish most vulnerable populations and reactivation versus recent transmission cases [10–13].
Hence, geospatial tools may be helpful to study the TB dynamics of urban areas with high prevalence of TB. Geospatial tools incorporating Geographic Information Systems (GIS) enable the identification and mapping of spatiotemporal clustering of disease or patients . The GIS method has been used to study the spatial distribution of human TB cases and has identified the heterogeneity of epidemic areas [10, 15]. To confirm recent TB transmission, isolates must be found to be clonal. Clonality must then be investigated by more discriminatory genotyping methods. In Antananarivo, an exhaustive thorough genomic characterization of clinical isolates either by 24 MIRU-VNTR or by WGS remains out of reach for economic reasons for the time-being. For this reason, we chose a classical spoligotyping approach, which remains a first-line method to characterize clinical isolates in resource-limited countries. Such a combination of approaches (spatial and genetic clustering) is interesting to locate spatial clusters of TB and attempt to assess where recent TB transmission cases may occur.
In Madagascar, the incidence of TB in 2013 was estimated to be approximately 233/100,000 inhabitants . TB prevalence distribution in Madagascar is likely heterogeneous with particularly high rates in specific areas driven by uncontrolled transmission, as in most resource-limited countries with either remote or isolated settings [15, 17, 18]. A crucial element of TB control efforts is the identification of “TB Hotspot areas” for the orientation of TB control strategies given the lack of resources. Therefore it would be beneficial for the TB program to have a tool for targeting areas of high transmission risk where interventions should be concentrated. Previous studies using genotyping techniques on clinical isolates have shown a large diversity of circulating M. tuberculosis genotypes in Madagascar [19–21]. A preliminary study based on TB notification rates identified spatial aggregation of TB cases in Antananarivo [15, 17]. This aggregation could be due either to actual transmission and/or to reactivation cases. Our aim here was to accumulate evidence concerning the possibility that previously identified TB hotspot areas in Antananarivo could be linked to transmission events.
We used a combination of spatial tools and genotyping to identify potential high-prevalence and likely higher recent transmission risk, TB areas in Antananarivo, Madagascar, as done recently in Brazil and Japan [11, 12]. These methods will be used in Madagascar for further understanding of the TB epidemiology in Antananarivo and the identification of priority targets for TB control strategies.
Fokontany were identified by conventional alphanumeric codes and were geolocalized by their centroid. Data of the population census of the 192 Fokontany were provided by the Development Office of Antananarivo (DOA) and poorest neighborhoods were described previously .
New confirmed TB cases for individuals were included without age limitation between August 2013 and May 2014 (n = 523). Only patients living in Antananarivo for at least one year have been recruited since they are more likely to have got TB in this study area. All 17 DTCs in the study area and three reference DTCs in the suburban zone of Antananarivo (Fenoarivo, Ambohibao, and Itaosy), totaling 20 DTCs were included in the study. The DTC of the central jail was excluded because of the potential transmission bias associated with this last area.
For each recruited patient, the objectives and nature of the study were explained by a DTC agent. Information written in Malagasy language was given to any patient or parents for children under 18 years. He/she was then asked if he/she agrees to participate in the study. Study subjects were allowed to refuse to participate in the study or withdraw from the study at any time without prior justification. When the subject agreed to participate in the study, he/she signed the informed consent form. An interview was conducted and the information was collected on a form. The residence address and Fokontany were recorded for every patient. The residence address provided by the patients in the questionnaire was confirmed by the DTC registers.
Fresh sputa were collected, stored, and cultured on Lowenstein Jensen (LJ) solid medium. Raw DNA was obtained from the clinical isolates by heating and killing a suspension of a culture colony in a hot dry bath at 80 °C for 30 min, After centrifugation at 13,000 * G for 5 min, supernatant was transferred in a new tube and conserved at −20 °C. Raw DNA was shipped at −20 °C to Orsay, France where the spoligotyping was performed directly using a Luminex 200 platform as described by Zhang . The spoligotype profiles were identified by their Shared-type (SIT – Spoligotyping International Types) and lineage designation by using the SpolDB4/SITVITWEB rules and classification . Moreover, to fit with current genome-based lineage designations, L1-L7 lineage labels were added to spoligotype families [24, 25]. Given the modest benefit of spoligotyping value when used alone in molecular epidemiological studies, we did not infer any recent transmission rate based on spoligotyping-based clustering, but used the genotyping information to correlate spatial and genetic clustering as a clue to define “Hot Spot TB area”.
TB patients were localized according to their Fokontany of residence. All data on the population denominators per Fokontany and metadata used for mapping were provided by the DOA. All TB patients and patients with genotypic clustered isolates (patients associated with isolates with repeated spoligotype or PRS) were scanned separately using the Kulldorff spatial scanning method. The populations per Fokontany with SaTScan™ were considered for the spatial analyses (http://www.satscan.org) . For Satscan® parameters, a maximum radius of 1 km was defined for the spatial scanning. A p-value <0.05 determined by conducting Monte Carlo replications was considered to be statistically significant. We assumed that the number of TB patients in each Fokontany fits with the Poisson distribution. The identified Fokontany associated with spatial clusters representing TB focal points of epidemics were mapped using QuantumGis 2.8® (QGIS Development Team, 2013). Areas representing focal points of epidemics with highly genetically related isolates (clustered by spoligotyping) were considered to be potential TB hotspot transmission areas.
Genotypic diversity of the clinical isolates
We obtained 88 individual spoligotype profiles from 394 clinical isolates (84.4% of the 467 recruited patients). We could not type 73 (16%) clinical isolates as 40 cultures were negative and 33 spoligotypes were non-interpretable. Of the individual 88 spoligotype profiles, 47 were unique and 41 were clustered, totaling 347 clustered isolates, found in clusters of 2 to 39 clinical isolates (Additional file 1: Table S1).
Residences of the included patients (n = 467) were distributed in 133 of the 192 Fokontany (Fig. 1c). Spatial observations showed that the highest number of TB cases was observed in Fokontany Andohatapenaka II in the 1st district (n = 15) whereas the district with the highest number of TB cases was the 5th district (112/467; 24.0%). The relatively most populous Fokontany at the periphery of Antananarivo had the higher TB cases number compared to the center of the city (3rd district) (Fig. 1b, c).
This study aimed to determine the spatial signature of TB by using a combination of genotyping and geospatial tools across the urban city of Antananarivo, Madagascar. The complementarity of both the genotyping and spatial analyses approaches has been used for the determination of TB transmission areas and risk studies in Brazil and Japan [11, 12]. These two methods are not new but they had never been combined for the detection of potential high risk areas of TB transmission in Antananarivo.
The combination of GIS and spoligotype data identified four hotspots with potential TB transmission in the city of Antananarivo. Analysis of all TB cases (without distinction of genotypic clustering of strains) identified a TB disease focal point overlapping with one of the potential transmission area (constituted by the Fokontany of Andohatapenaka II). TB hotspots such as that detected by the analysis taking only spatial data into account is likely constituted simultaneously by patients linked by indirect and relatively ancient TB transmission and by patients linked by local and recent TB transmission. Spatial clustering of patients with genotypically clustered isolates may therefore concern more patients associated with local and recent transmission. Additionally, our combination of approaches permitted the detection of spatial clusters of TB patients which were not detected with only the spatial scan of TB cases. TB transmissions might have occurred in these areas. For definitive proof of TB transmission, more discriminatory genotyping tools should however be used.
While previous studies determined risk factors associated with spatial clustering of TB cases in Antananarivo [15, 18], risk factors associated with transmission were poorly investigated. The first potential transmission hotspot is constituted by the Fokontany of Antohomadinika Afovoany that is localized in one of the poorest neighborhood in the 1st urban district of Antananarivo. Most of the houses in this area are made of precarious wooden hovels and the majority of the local population does not have standard health care access. While life under fragile conditions and environmental factors [14, 19] are known factors that contribute to TB reactivation , they could also foster local transmission as suggested by this study. The relatively high diversity of SITs seen in this poor area supports the previous findings that environmental factors also favor reactivation of latent TB. This area, containing the highest rate of TB, constitutes the first high risk area of TB transmission in Antananarivo.
The three other areas of potential TB transmission (#2, #3 and #4) host or lie close to three public markets (Anosibe, Isotry and Ambanidia market) where there is an important flow of people. Similar studies have shown that a flow of persons promotes TB transmission [11, 12]. This study further adds evidence for potential transmission in markets.
The discriminatory level of spoligotyping is relatively low, and the genotypic clustering cannot be used to estimate TB recent transmission rates. Thus isolates with the same spoligotype may be coming from different chains of transmission. However, when compared with isolates with single genotype, that might represent reactivation or relatively far recent transmission, these clustered isolates are more likely to have originated from recent spatially localized transmission. The aim of the study was to identify high risk areas, this limitation does therefore not invalidate this study and should be taken as a first step towards identification of real TB transmission hotspots.
From 1994 until 2000, studies of the M. tuberculosis genotype profiles present in Antananarivo reported six lineages and sub-lineages of M. tuberculosis clinical isolates: T, LAM, and H (Lineage 4), EAI (Lineage 1), CAS (Lineage 3), and Beijing (Lineage 2) . The 6 lineages were shown to be present in the capital of Madagascar in 2004/2005, with the appearance of other minor sub-lineages such as S and X along with some unknowns (U) . The distribution of the major circulating clinical isolates observed in these studies may have evolved slightly, but we noted no significant change since then, suggesting that all the lineages were being continuously transmitted in this area.
In this study, the population density in every studied Fokontany was taken into account, increasing the accuracy of spatial clustering identification. The use of spoligotyping as a first screen is a relatively simple and inexpensive method. It allowed us to identify clusters that were overlooked using only spatial information. We thus plan to keep this strategy as a first-line detection of potential transmission areas. A possible improvement in our approach would be the use of methods providing both spoligotyping and resistance data such as TB-SPRINT . This high-throughput assay tests for rifampicin and isoniazid resistance simultaneously with spoligotyping on Luminex. It could be a useful strategy in the near future to survey and prevent the spread of MDR-TB cases.
Another limit of this study was the short recruitment time (9 months). However, this recruitment period was sufficient to achieve a large sample size, and given the stability of the population, such duration should not have created too much bias. Some of the enrolled patients did not agree to give consent to participate in the study, those residing in Antananarivo Renivohitra consulting with DTCs outside Antananarivo, and those not consulting with a DTC had to be removed from the inclusion. Finally, we chose the patient residence to perform the spatial analysis although the patient residence is clearly not the unique site of possible TB transmission. Activity areas have been linked to transmission in other studies [11, 12]. We still chose to locate TB cases according to the residential Fokontany since people with most poor living condition and without work circulate most of their time in the same Fokontany in Antananarivo. A further investigation using the combined spatial methods on both working and residential areas could still be useful for TB epidemiological surveillance.
Areas of high risk of TB transmission were suggested by the combination of M. tuberculosis spoligotyping and spatial analysis in the Urban District of Antananarivo. This method may be helpful for TB epidemiological surveillance in Madagascar and developing countries for guiding TB control strategies by the identification of priority target areas. Massive active diagnosis for children and susceptible persons can be suggested in these areas, as well as intense of TB awareness.
The “Projet interne” of the Institut Pasteur de Madagascar funded the work. Our thanks go to the “Programme National de Lutte Contre la Tuberculose” represented by the Dr. Andrianantenaina RAKOTOSON for the authorization to conduct this study, the officials and staff of the various “Centre de Diagnostic et de Traitement de la Tuberculose”, the Mycobacteria Unit of the IPM, the SIG cell of the Epidemiology Unit of IPM, and the IGEPE team of the Institute of Integrative Cell Biology, UMR9198, CEA-CNRS-University Paris Saclay, France, for their collaboration. We are also grateful to the TB patients who accepted to participate and provided biological samples to perform this study and the Development Office of Antananarivo (DOA) for the different datas used for the study.
Availability of data and materials
All data generated or analyzed during this study are included in this published article and in the Additional file 1: Table S1.
NHR: Initiated the research, wrote the study design, conducted the research, did data entry and analysis and wrote the manuscript. NR: Was involved in the data analysis and in the write up of the manuscript. SLR: Conducted a part of the research and participated in data analysis. MSR: Participated in the study design and in the conducting of the research. GR: Initiated the research. CS: Served as co-supervisor, initiated the research and write up the manuscript. FR: Served as co-supervisor, initiated the research involved in data analysis and write up the manuscript. VR: Served as main supervisor, initiated the research, involved in data analysis and write up the manuscript. All authors read and approved the final manuscript.
Ethics approval and consent to participate
The study was approved by the ethics committee of the Malagasy Ministry of Public Health (Authorization No: 057MSANP/EC – 04th July 2013).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.
- WHO. Global Tuberculosis report 2015; 2015.Google Scholar
- WHO. WHO Report 2014. Global tuberculosis control 2014; 2014.Google Scholar
- Kamerbeek J, Schouls L, Kolk A, van Agterveld M, van Soolingen D, Kuijper S, et al. Simultaneous detection and strain differentiation of Mycobacterium tuberculosis for diagnosis and epidemiology. J Clin Microbiol. 1997;35(4):907–14.PubMedPubMed CentralGoogle Scholar
- Demay C, Liens B, Burguiere T, Hill V, Couvin D, Millet J, et al. SITVITWEB--a publicly available international multimarker database for studying Mycobacterium tuberculosis genetic diversity and molecular epidemiology. Infect Genet Evol. 2012;12(4):755–66.View ArticlePubMedGoogle Scholar
- Brudey K, Driscoll JR, Rigouts L, Prodinger WM, Gori A, Al-Hajoj SA, et al. Mycobacterium tuberculosis complex genetic diversity: mining the fourth international spoligotyping database (SpolDB4) for classification, population genetics and epidemiology. BMC Microbiol. 2006;6:23.View ArticlePubMedPubMed CentralGoogle Scholar
- Comas I, Homolka S, Niemann S, Gagneux S. Genotyping of genetically monomorphic bacteria: DNA sequencing in Mycobacterium tuberculosis highlights the limitations of current methodologies. PLoS One. 2009;4(11):e7815.Google Scholar
- Stucki D, Ballif M, Egger M, Furrer H, Altpeter E, Battegay M, et al. Standard genotyping overestimates transmission of Mycobacterium tuberculosis among immigrants in a low incidence country. J Clin Microbiol. 2016;54(7):1862-70.Google Scholar
- Supply P, Allix C, Lesjean S, Cardoso-Oelemann M, Rusch-Gerdes S, Willery E, et al. Proposal for standardization of optimized mycobacterial interspersed repetitive unit-variable-number tandem repeat typing of Mycobacterium tuberculosis. J Clin Microbiol. 2006;44(12):4498–510.View ArticlePubMedPubMed CentralGoogle Scholar
- Small PM, Hopewell PC, Singh SP, Paz A, Parsonnet J, Ruston DC, et al. The epidemiology of tuberculosis in San Francisco-a population-based study using conventional and molecular methods. N Engl J Med. 1994;330(24):1703–9.View ArticlePubMedGoogle Scholar
- Moonan PK, Ghosh S, Oeltmann JE, Kammerer JS, Cowan LS, Navin TR. Using genotyping and geospatial scanning to estimate recent Mycobacterium tuberculosis transmission, United States. Emerg Infect Dis. 2012;18(3):458–65.View ArticlePubMedPubMed CentralGoogle Scholar
- 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(9):e0138831.View ArticlePubMedPubMed CentralGoogle Scholar
- Ribeiro FK, Pan W, Bertolde A, Vinhas SA, Peres RL, Riley L, et al. Genotypic and spatial analysis of Mycobacterium tuberculosis transmission in a high-incidence urban setting. Clin Infect Dis. 2015;61(5):758-66.Google Scholar
- Gurjav U, Burneebaatar B, Narmandakh E, Tumenbayar O, Ochirbat B, Hill-Cawthorne GA, et al. Spatiotemporal evidence for cross-border spread of MDR-TB along the Trans-Siberian Railway line. Int J Tuberc Lung Dis. 2015;19(11):1376–82.View ArticlePubMedGoogle Scholar
- Kulldorff M, Nagarwalla N. Spatial disease clusters: detection and inference. Stat Med. 1995;14(8):799–810.View ArticlePubMedGoogle Scholar
- Randremanana RV, Sabatier P, Rakotomanana F, Randriamanantena A, Richard V. Spatial clustering of pulmonary tuberculosis and impact of the care factors in Antananarivo City. Tropical Med Int Health. 2009;14(4):429–37.View ArticleGoogle Scholar
- WHO. Tuberculose: profil de pays/Madagascar. 2014.Google Scholar
- Rakotosamimanana S, Mandrosovololona V, Rakotonirina J, Ramamonjisoa J, Ranjalahy JR, Randremanana RV, et al. Spatial analysis of pulmonary tuberculosis in Antananarivo Madagascar: tuberculosis-related knowledge, attitude and practice. PLoS ONE. 2014;9(11):e110471.View ArticlePubMedPubMed CentralGoogle Scholar
- Randremanana RV, Richard V, Rakotomanana F, Sabatier P, Bicout DJ. Bayesian mapping of pulmonary tuberculosis in Antananarivo, Madagascar. BMC Infect Dis. 2010;10(1):21.View ArticlePubMedPubMed CentralGoogle Scholar
- Ferdinand S, Sola C, Chanteau S, Ramarokoto H, Rasolonavalona T, Rasolofo-Razanamparany V, et al. A study of spoligotyping-defined Mycobacterium tuberculosis clades in relation to the origin of peopling and the demographic history in Madagascar. Infect Genet Evol. 2005;5(4):340–8.View ArticlePubMedGoogle Scholar
- Razanamparany VR, Ménard D, Aurégan G, Gicquel B, Chanteau S. Extrapulmonary and pulmonary tuberculosis in Antananarivo (Madagascar): high clustering rate in female patients. J Clin Microbiol. 2002;40(11):3964–9.View ArticleGoogle Scholar
- Rasolofo-Razanamparany V, Auregan G, Ratsirahonana O, Raharimanana R, Ramarokoto H, Gicquel B, et al. Genetic polymorphism of M. tuberculosis strains in Antananaviro. Arch Inst Pasteur Madagascar. 1994;62(1):59–64.Google Scholar
- Wachsberger J-M. Les quartiers pauvres à Antananarivo. Autrepart. 2009(3):117–37. http://www.cairn.info/revue-autrepart-2009-3-page-117.htm. Accessed 3 June.
- Zhang J, Abadia E, Refregier G, Tafaj S, Boschiroli ML, Guillard B, et al. Mycobacterium tuberculosis complex CRISPR genotyping: improving efficiency, throughput and discriminative power of ‘spoligotyping’ with new spacers and a microbead-based hybridization assay. J Med Microbiol. 2010;59(Pt 3):285–94.View ArticlePubMedGoogle Scholar
- Coll F, McNerney R, Guerra-Assunçào JA, Glynn JR, Perdigào J, Viveiros M, et al. A robust SNP barcode for typing Mycobacterium tuberculosis complex strains. Nat Commun. 2014;5:4812.View ArticlePubMedPubMed CentralGoogle Scholar
- Coll F, Preston M, Guerra-Assunçào JA, Hill-Cawthorn G, Harris D, Perdigào J, et al. PolyTB: a genomic variation map for Mycobacterium tuberculosis. Tuberculosis. 2014;94(3):346–54.View ArticlePubMedPubMed CentralGoogle Scholar
- Kulldorff M. A spatial scan statistic. Communications in statistics. 1997;26:1481–96.View ArticleGoogle Scholar
- Rakotosamimanana N, Raharimanga V, Andriamandimby SF, Soares JL, Doherty TM, Ratsitorahina M, et al. Variation in gamma interferon responses to different infecting strains of Mycobacterium tuberculosis in acid-fast bacillus smear-positive patients and household contacts in Antananarivo, Madagascar. Clin Vaccine Immunol. 2010;17(7):1094–103.View ArticlePubMedPubMed CentralGoogle Scholar
- Gomgnimbou MK, Hernandez-Neuta I, Panaiotov S, Bachiyska E, Palomino JC, Martin A, et al. Tuberculosis-spoligo-rifampin-isoniazid typing: an all-in-one assay technique for surveillance and control of multidrug-resistant tuberculosis on Luminex devices. J Clin Microbiol. 2013;51(11):3527–34.View ArticlePubMedPubMed CentralGoogle Scholar