JCM Accepts, published online ahead of print on 29 February 2012 J. Clin. Microbiol. doi:10.1128/jcm.00132-12 Copyright 2012, American Society for Microbiology. All Rights Reserved. 1 2 Real-time molecular epidemiology of tuberculosis by direct genotyping of smear-positive clinical specimens 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 María Alonso 1,2,3, Marta Herranz 1,2,3, Miguel Martínez Lirola 4,, Milagros González-Rivera 2,5, Emilio Bouza 1,2,3, Darío García de Viedma* 1,2,3 1 Servicio de Microbiología y Enfermedades Infecciosas, Hospital General Universitario Gregorio Marañón, Madrid, Spain 2 Instituto de Investigación Sanitaria Gregorio Marañón, Madrid, Spain 3 CIBER Enfermedades Respiratorias-CIBERES, Spain 4 Complejo Hospitalario Torrecárdenas, Almería, Spain. On behalf of INDAL-TB group 5 Línea Instrumental Secuenciación, Hospital Gregorio Marañón, Madrid, Spain Running title: Real-time molecular epidemiology in tuberculosis *Corresponding author: Darío García de Viedma Servicio de Microbiología y Enfermedades Infecciosas Hospital General Universitario Gregorio Marañón C/ Dr Esquerdo, 46, 28007 Madrid, Spain Fax: 91 5044906 Email: dgviedma2@gmail.com 1
24 25 26 27 28 29 30 31 Abstract We applied MIRU-VNTR to directly analyze the bacilli present in 61 stainpositive specimens from tuberculosis patients. A complete MIRU-type (24 loci) was obtained from all but one (no amplification in one locus) of the specimens (98.4%), and fully correlated with those obtained from the corresponding cultures. Our study is the first to demonstrate that real-time genotyping of Mycobacterium tuberculosis can be achieved, fully transforming the way in which molecular epidemiology can be integrated into control programs. Downloaded from http://jcm.asm.org/ on December 18, 2018 by guest 2
32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 The application of fingerprinting tools has proven crucial in the identification of cases infected by the same Mycobacterium tuberculosis strain (clustered cases), which are considered to result from recent transmission events and constitute a key indicator when evaluating the efficiency of tuberculosis control programs (1). The time to identify clusters from cultured samples has been greatly reduced by a new, faster PCRbased technique, MIRU-VNTR (7). However, one challenge remains unresolved, namely, how to obtain the MTB genotype sufficiently quickly for it to be integrated in the survey of tuberculosis transmission, while contact tracing is still ongoing. To get this target the only option would be to switch from a culture-based genotyping approach to one based on analyzing directly the bacilli present in respiratory specimens and this is the aim of this study. We selected all cases (January 2005-December 2010) of culture-confirmed TB with a stain-positive sputum identified in a population-based molecular epidemiology survey on transmission of tuberculosis performed in Almeria (3). We established the following criteria of inclusion: a) to have frozen stain-positive sputum available; b) to have the complete MIRU-VNTR genotype from their corresponding cultured MTB isolate available. Among them, we selected representatives of clustered cases (infected by the same strain) and orphan cases, attempting to maintain a proportion between clustered and orphan cases similar to that found during the survey in Almeria (around 35:65). The estimated sample size was based on the following assumptions: expected proportion of complete agreement (absolute concordance) =.90; width of the confidence interval: +/- 10% (.10); alpha two-sided error = 0.05. With these assumptions, the estimated sample size (estimation of one proportion) was 35. In order to increase the precision of our estimation we decided to extend the recruitment to all the cases (N: 61) who fulfilled the inclusion criteria mentioned above. 3
57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 Specimens were retrieved from the Complejo Hospitalario Torrecardenas in Almeria and sent blind to Hospital Gregorio Marañon in Madrid for analysis. The final genotyping data obtained from the analysis on specimens were sent back to Almeria to perform a concordance analysis with the data previously obtained from the cultured MTB isolates. DNA was extracted from 1 ml of decontaminated sputum using a column-based purification method (QIAamp DNA Mini Kit protocol, Qiagen, Courtaboeuf, France) and eluted in 70 μl of buffer AE. The MIRU-VNTR protocol followed was a modified version of the original one (eight triplex PCRs) developed to analyze cultured isolates (6). The final reaction mixture (50 μl) included 25 μl of PCR Master Mix (QIAGEN Multiplex PCR kit, Qiagen, Courtaboeuf, France), 5 μl of Q solution, and 0.25 μm of each unlabeled and labeled oligonucleotide for mixes 1-3, 5, 6 and 8. Mixes 4 and 7 were analyzed with the PuReTaq Ready-To-Go PCR Beads System (GE Healthcare, Buckinghamshire, UK), and the final reaction mixtures (25 μl) included 0.25 μm of each unlabeled and labeled oligonucleotide (0.5 μm for loci QUB-4156 and MIRU20) and 1.5 μl of dimethyl sulfoxide. The primers used for PCR amplification were described by Supply et al (6) and amplification profiles were as described elsewhere (5), except for the number of cycles in multiplex PCR assays (35 cycles). PCR products were purified with the High Pure PCR product Purification Kit (Roche Diagnostics, Mannheim, Germany), and their concentration was measured and diluted, to rule out signal saturation, at 2.5 ng/μl for all the mixes except 4 and 7, in which were fixed at 10 ng/μl. PCR products were analyzed by capillary electrophoresis (3130xl GA- POP-7 polymer - 1200 LIZ ISS - Applied Biosystems, Foster City, California, USA). The PCR 4
82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 fragments were sized using the GeneMapper v 4.0 software package. Throughout the study, technical performance was evaluated by monitoring a blank and two fixed PCR samples as quality controls. The retention time data from the ISS peaks and the amplified product peaks of the two quality controls were used to rule out deviations in the correct allelic assignment of the unknown samples. We performed MIRU-VNTR analysis on 61 stain-positive specimens from 61 patients. Based on the number of acid-fast bacilli observed (4), clinical specimens were classified as having >90 bacilli/field (4+; N=30), 10-90 b/f (3+; N=15), 1-9 b/f (2+; N=11) and 1-9 b/10f (1+, N=5). In 7 specimens (two 4+, three 3+ and two 2+), amplification was not detected in one/two of the loci of one of the mixes (three specimens in locus Mtub39, two in locus QUB-4156, one in loci Mtub39/ QUB-4156 and one in loci ETRA/QUB-11b; Table). In these cases the determination was repeated applying a single PCR for each specific locus. The simplex-pcr reaction was run at a final reaction mixture (50 μl) composed of 1 U of HotStart Taq DNA Polymerase (Qiagen), 5 μl of PCR buffer plus 3 μl of MgCl 2 with a final Mg 2+ concentration of 2.25 mm, 10 μl of Q-solution, 0.4 μl of dntp Mix, and 0.1 μm of each unlabeled and labeled oligonucleotide. After applying simplex-pcr, the locus QUB11b from patient 48 still failed to be amplified but all but one of the remaining pendant loci rendered a result, which meant the obtention of a complete genotype for all but one of the 61 specimens (98,36%; CI95%: 91,2-99,9) (Table), even from the specimens with the lowest bacterial load (1-9 b/10 f). When we compared our data with the MIRUtype obtained by analyzing the corresponding cultured isolates, the results were identical for all of them. The analysis of all patterns enabled us to identify 23 patients grouped in nine clusters (Table, C1-C9, 5
106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 including 2-4 cases) and 38 orphan cases (cases 24-61; Table), which coincided with the clustered/orphan distribution obtained from the analysis of the cultured isolates. Our study is the first evaluation, with molecular-epidemiology purposes, of the feasibility of a tuberculosis genotyping scheme applied directly to clinical specimens. If MIRU-VNTR had been applied directly on specimens, we would have succeeded in classifying isolates as clustered or orphan, around 30 days in advance, without waiting for the MTB culture results. However we must admit that could there be a gap between our study and the routine performance of applying MIRU-VNTR on clinical samples. The size of this gap should be measured after final evaluation of the true potential of our proposal in forthcoming studies which evaluate the feasibility to introduce this scheme into a prospective routine practice Genotyping MTB directly from uncultured clinical specimens would also make it possible to extend the applicability of molecular epidemiology approaches to developing areas where MTB culture is not systematically performed and standard genotyping cannot be performed, and, consequently, precise information on tuberculosis transmission dynamics is lacking. Traditional retrospective descriptive epidemiology will soon be replaced by interventionist molecular epidemiology (2), a new format which attempts to rapidly integrate available molecular cluster data to optimize tuberculosis control programs by offering first-line information on transmission routes and dynamics. Our data show that a revolution in the way we understand the molecular epidemiology of tuberculosis is now possible. 6
129 Acknowledgments 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 María Alonso is supported by a contract (REF CA09/00054) from the Instituto de Salud Carlos III (Fondo de Investigaciones Sanitarias) and provides technical support in the Unidad Central de Análisis Molecular of Hospital General Universitario Gregorio Marañón. LightCycler 2.0 was acquired with a grant (REF IF08/36173) from Instituto de Salud Carlos III (Fondo de Investigaciones Sanitarias). The 3130xl Genetic Analyzer was partially financed by grants from Fondo de Investigaciones Sanitarias (IF01-3624, IF08-36173). The study was partially financed by Fondo de Investigaciones Sanitarias (S09/02205), Junta de Andalucía (PI-0444/2008 and PI-0306-2009), and SEPAR (763-09). We are grateful to Ainhoa Simón Zárate, the sequencer technician, who holds a contract from the Fondo de Investigaciones Sanitarias (CA08/00160). We would like to thank Thomas O Boyle for proofreading the manuscript. 7
146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 References 1. Cattamanchi, A., P. C. Hopewell, L. C. Gonzalez, D. H. Osmond, L. Masae Kawamura, C. L. Daley, and R. M. Jasmer. 2006. A 13-year molecular epidemiological analysis of tuberculosis in San Francisco. Int J Tuberc Lung Dis 10:297-304. 2. Garcia de Viedma, D., I. Mokrousov, and N. Rastogi. 2011. Innovations in the molecular epidemiology of tuberculosis. Enferm Infecc Microbiol Clin 29 Suppl 1:8-13. 3. Martinez-Lirola, M., N. Alonso-Rodriguez, M. L. Sanchez, M. Herranz, S. Andres, T. Penafiel, M. C. Rogado, T. Cabezas, J. Martinez, M. A. Lucerna, M. Rodriguez, C. Bonillo Mdel, E. Bouza, and D. Garcia de Viedma. 2008. Advanced survey of tuberculosis transmission in a complex socioepidemiologic scenario with a high proportion of cases in immigrants. Clin Infect Dis 47:8-14. 4. Nolte, F. 1995. Mycobacterium. Manual of Clinical Microbiology 6ed. Washington, DC. 5. Oelemann, M. C., R. Diel, V. Vatin, W. Haas, S. Rusch-Gerdes, C. Locht, S. Niemann, and P. Supply. 2007. Assessment of an optimized mycobacterial interspersed repetitive- unit-variable-number tandem-repeat typing system combined with spoligotyping for population-based molecular epidemiology studies of tuberculosis. J Clin Microbiol 45:691-7. 6. Supply, P., C. Allix, S. Lesjean, M. Cardoso-Oelemann, S. Rusch-Gerdes, E. Willery, E. Savine, P. de Haas, H. van Deutekom, S. Roring, P. Bifani, N. Kurepina, B. Kreiswirth, C. Sola, N. Rastogi, V. Vatin, M. C. Gutierrez, M. Fauville, S. Niemann, R. Skuce, K. Kremer, C. Locht, and D. van Soolingen. 2006. Proposal for standardization of optimized mycobacterial interspersed 8
171 172 173 174 repetitive unit-variable-number tandem repeat typing of Mycobacterium tuberculosis. J Clin Microbiol 44:4498-510. 7. Supply, P., S. Lesjean, E. Savine, K. Kremer, D. van Soolingen, and C. Locht. 2001. Automated high-throughput genotyping for study of global 175 176 177 178 epidemiology of Mycobacterium tuberculosis based on interspersed repetitive units. J Clin Microbiol 39:3563-71. mycobacterial Downloaded from http://jcm.asm.org/ on December 18, 2018 by guest 9
Table. MIRU-VNTR data from analysis of respiratory specimens MIRU-VNTR loci aliases a Clustered strains Cluster status Patient Bacterial load MIRU2 MIRU20 MIRU23 MIRU24 MIRU27 MIRU39 MIRU4 MIRU26 MIRU40 MIRU10 MIRU16 MIRU31 Mtub04 ETR C ETR A Mtub30 Mtub39 QUB-4156 QUB-11b Mtub21 QUB-26 Mtub29 ETR B Mtub34 Full agreement with cultured isolates C1 1 +++ 2 2 6 1 3 2 2 3 4 4 3 3 3 4 2 2 1 2 4 2 5 4 2 1 + C1 2 +++ 2 2 6 1 3 2 2 3 4 4 3 3 3 4 2 2 1 2 4 2 5 4 2 1 + C2 3 ++ 2 2 6 1 3 2 2 4 1 4 2 3 2 4 2 1 2 2 2 3 4 4 1 5 + C2 4 +++ 2 2 6 1 3 2 2 4 1 4 2 3 2 4 2 1 2 2 2 3 4 4 1 5 + C2 5 +++ 2 2 6 1 3 2 2 4 1 4 2 3 2 4 2 1 2 2 2 3 4 4 1 5 + C3 6 ++++ 1 2 6 1 3 2 2 5 2 5 3 2 3 2 2 1 2 2 2 3 4 4 2 3 + C3 7 ++++ 1 2 6 1 3 2 2 5 2 5 3 2 3 2 2 1 2 2 2 3 4 4 2 3 + C4 8 ++++ 2 1 3 1 3 2 2 5 1 5 3 3 2 3 3 4 3 3 6 3 4 2 2 3 + C4 9 +++ 2 1 3 1 3 2 2 5 1 5 3 3 2 3 3 4 3 3 6 3 4 2 2 3 + C4 10 ++++ 2 1 3 1 3 2 2 5 1 5 3 3 2 3 3 4 3 3 6 3 4 2 2 3 + C4 11 + 2 1 3 1 3 2 2 5 1 5 3 3 2 3 3 4 3 3 6 3 4 2 2 3 + C5 12 +++ 2 2 3 1 4 2 3 6 5 4 3 3 3 4 2 2 3 2 4 3 5 5 2 3 + C5 13 ++++ 2 2 3 1 4 2 3 6 5 4 3 3 3 4 2 2 3 2 4 3 5 5 2 3 + C6 14 ++++ 2 2 5 1 3 2 2 4 1 3 3 3 2 4 2 2 5 2 4 2 4 4 2 3 + C6 15 ++ 2 2 5 1 3 2 2 4 1 3 3 3 2 4 2 2 5 2 4 2 4 4 2 3 + C7 16 ++++ 2 2 5 1 3 2 2 5 2 5 3 3 2 3 3 4 3 1 5 3 5 4 2 3 + C7 17 ++++ 2 2 5 1 3 2 2 5 2 5 3 3 2 3 3 4 3 1 5 3 5 4 2 3 + C8 18 ++++ 2 2 5 1 3 2 2 5 3 6 3 3 2 3 3 4 2 3 6 2 6 4 2 3 + C8 19 ++++ 2 2 5 1 3 2 2 5 3 6 3 3 2 3 3 4 2 3 6 2 6 4 2 3 + C9 20 + 2 2 5 1 3 2 2 5 4 3 2 3 2 4 3 2 3 2 4 2 5 4 2 3 + C9 21 ++++ 2 2 5 1 3 2 2 5 4 3 2 3 2 4 3 2 3 2 4 2 5 4 2 3 + C9 22 ++++ 2 2 5 1 3 2 2 5 4 3 2 3 2 4 3 2 3 2 4 2 5 4 2 3 + C9 23 +++ 2 2 5 1 3 2 2 5 4 3 2 3 2 4 3 2 3 2 4 2 5 4 2 3 + O1 24 ++++ 3 2 4 2 3 2 2 4 1 5 2 5 2 5 4 4 4 3 3 4 6 3 4 3 + O2 25 ++++ 2 2 5 1 3 2 3 5 4 3 2 3 2 4 3 2 3 2 5 2 5 4 2 3 + O3 26 ++++ 2 2 5 1 3 2 2 5 3 3 3 3 2 3 3 4 4 2 4 3 7 4 2 3 + O4 27 +++ 2 2 3 1 3 2 2 5 3 5 3 3 2 3 3 4 3* 3* 5 3 5 2 2 3 + O5 28 +++ 2 2 5 1 4 2 2 5 3 5 3 3 3 3 3 4 3 3 6 2 4 4 2 2 + O6 29 ++++ 1 2 6 1 3 2 2 5 1 3 3 3 2 4 3 2 6* 2 5 2 5 4 2 3 + O7 30 ++++ 2 2 4 2 3 2 2 4 1 6 2 5 2 5 7 4 4 3 3 4 6 3 5 3 + O8 31 ++++ 1 2 3 1 3 2 2 5 5 4 3 2 3 2 2 1 2 2 2 3 6 4 2 3 + O9 32 ++++ 2 2 5 1 3 2 2 5 3 5 3 3 2 3 3 4 3 3 5 3 5 4 2 3 + O10 33 ++++ 2 2 6 1 3 2 2 4 1 3 2 3 4 4 2 1 2 2 4 3 8 4 1 5 + O11 34 +++ 2 1 5 1 3 3 1 5 3 3 1 4 2 4 3 2 3 2 3 2 5 4 3 3 + O12 35 ++ 2 2 5 1 3 2 2 4 0 3 4 3 2 4 2 2 3 2* 5 2 5 4 2 3 + O13 36 ++ 2 2 5 1 3 3 2 5 3 3 3 5 4 4 3 4 3 2 6 5 8 4 2 3 + O14 37 ++ 2 2 5 1 3 2 2 5 3 4 3 3 2 3 3 4 3 3 6 3 7 4 2 3 + O15 38 ++++ 2 2 3 1 3 2 2 6 3 5 3 3 2 3 3 3 3 3 4 2 7 2 2 3 + O16 39 ++ 1 2 6 1 3 2 1 5 5 2 3 3 5 4 2 1 2 3* 2 3 8 4 2 3 + O17 40 +++ 2 2 5 1 1 2 2 5 2 3 3 3 2 4 2 2 3 2 4 2 5 4 2 3 + O18 41 ++++ 2 2 3 1 3 2 2 5 3 5 3 3 2 3 3 4 3 3 4 3 5 2 2 3 + O19 42 ++++ 2 2 6 1 1 1 2 5 1 3 3 3 2 4 3 2 6 2 4 2 5 4 2 3 + O20 43 + 2 2 5 1 3 2 2 5 3 4 3 3 2 3 3 4 4 3 5 3 7 4 2 3 + O21 44 ++++ 2 2 3 1 3 2 2 5 4 5 3 3 2 3 3 4 3 3 6 3 7 2 2 3 + O22 45 ++++ 1 2 5 1 3 2 2 5 5 4 3 2 3 2 2 1 2 2 2 3 6 4 2 3 + O23 46 ++++ 2 2 5 1 3 2 2 6 4 4 3 3 2 3 3 4 3 3 3 2 7 4 2 3 + O24 47 ++ 2 1 3 1 3 2 2 5 3 5 3 3 2 3 3 4 3 3 6 3 7 2 2 3 + O25 48 ++++ 2 2 5 1 3 3 2 7 3 5 4 5 5 2 TA 2 3 4 4 8 4 2 3 NA O26 49 +++ 2 2 5 1 3 2 2 5 4 5 2 3 2 3 3 4 3 3 4 3 5 4 2 3 + O27 50 + 2 2 5 1 3 2 2 6 3 4 1 3 2 3 3 4 3 1 3 2 6 4 2 2 + O28 51 ++++ 2 2 3 1 3 2 2 5 4 5 2 3 2 3 3 3 4 3 3 3 7 2 2 3 + O29 52 ++ 2 2 5 1 3 3 2 7 3 3 3 5 4 4 4 4 3 2 5 5 8 4 2 3 + O30 53 + 2 1 5 1 3 2 2 5 4 5 3 3 2 3 3 4 3 3 6 3 7 4 2 3 + O31 54 ++++ 2 2 5 1 3 2 2 5 3 4 3 3 2 3 4 4 3 3 5 3 7 4 2 2 + O32 55 ++ 1 2 7 1 1 2 2 5 6 4 3 3 3 4 2 1 2 3 2 3 7 4 2 3 + O33 56 ++ 1 2 6 1 3 2 2 5 4 4 2 4 3 4 2 1 2 3 2 3 6 4 2 3 + O34 57 +++ 2 1 5 1 3 2 2 5 3 3 3 3 2 4 4 2 6* 2 6 3 5 4 2 3 + O35 58 +++ 2 1 3 1 3 2 2 5 3 5 3 3 2 3 3 4 5* 3 5 3 7 2 2 3 + O36 59 +++ 1 2 6 1 3 2 2 6 4 4 3 3 3 4 2 1 2 3 2 3 5 4 2 3 + O37 60 ++ 2 2 6 1 3 2 2 3 4 5 2 3 3 4 TA 2 1 2 2 2 7 4 2 1 + O38 61 ++++ 2 2 6 1 3 2 2 4 1 4 2 3 2 4 2 1 2 2 3 3 3 4 1 5 + TA: Target Absence. NA : No amplification of the alelle 2 in locus QUB-11b. *: Locus amplification required simplex-pcr format. a : Loci order according to MIRU-VNTRplus database (http://www.miru-vntrplus.org). CN: Cluster number, ON: Orphan strain number Orphan strains