Emerg Infect DisEmerging Infect. DisEIDEmerging Infectious Diseases1080-60401080-6059Centers for Disease Control and Prevention23621956364767212-145310.3201/eid1903.121453DispatchDispatchPrioritizing Tuberculosis Clusters by Genotype for Public Health Action, Washington, USATuberculosis Clusters, WashingtonLindquistScottAllenSheanneFieldKimGhoshSmitaHaddadMaryam B.NaritaMasahiroOrenEyalAuthor affiliations: Kitsap County Health District, Bremerton, Washington, USA (S. Lindquist); Washington State Department of Health, Olympia, Washington, USA (S. Lindquist, S. Allen, K. Field); Centers for Disease Control and Prevention, Atlanta, Georgia, USA (S. Ghosh, M.B. Haddad); Public Health–Seattle & King County Tuberculosis Control Program, Seattle, Washington, USA (M. Narita, E. Oren).Address for correspondence: Scott Lindquist, Kitsap County Health District, Kitsap Public Health District, 345 6th St, Suite 300, Bremerton, WA 98337, USA; email: scott.lindquist@kitsappublichealth.org32013193493496

Groups of tuberculosis cases with indistinguishable Mycobacterium tuberculosis genotypes (clusters) might represent recent transmission. We compared geospatial concentration of genotype clusters with independent priority rankings determined by local public health officials; findings were highly correlated. Routine use of geospatial statistics could help health departments identify recent disease transmission.

Keywords: tuberculosisgenotypedisease outbreaksclusterstuberculosis and other mycobacteriaUnited StatesWashingtonTBgenotypingMycobacterium tuberculosis

Mycobacterium tuberculosis genotyping has been applied to tuberculosis (TB) control activities for >2 decades, and epidemiologic or genotyping data can confirm or disprove outbreaks (14). Investigation of genotype clusters can identify unrecognized transmission and lead to interventions that interrupt further transmission (5,6). However, cluster investigations are complex, requiring patient interviews and field observations. Focusing resources on clusters that most likely represent recent TB transmission could reduce the number of unnecessary investigations.

Geospatial statistics can identify higher-than-expected concentrations of TB cases with indistinguishable genotypes (7). We describe a comparison of a quantitative geospatial statistic analysis with qualitative expert opinion for prioritizing TB cluster investigations in Washington, USA, a state with moderate TB incidence (3.5 cases/100,000 persons) (8). The comparison was performed for initial and follow-up 3-year periods, 2005–2007 (period 1) and 2008–2010 (period 2).

The Study

TB genotype clusters were defined as groups of >3 TB case-patients whose isolates had matching spoligotyping and 12-locus mycobacterial interspersed repetitive unit–variable number tandem repeat (MIRU-VNTR) (9) genotyping results that were reported in the same county within Washington. A log-likelihood ratio (LLR) was calculated for each genotype cluster identified during each of the two 3-year periods (Figure). The larger the LLR, the greater the possibility the cluster represented geographically concentrated TB cases, a proxy for recent TB transmission. The cutoff point for the LLR was set to 5.0, based on the value used by the national TB Genotyping Information Management System (10).

Formula used to calculate geospatial statistic (a modified log-likelihood ratio [LLR]) on the basis of geographic distribution of Mycobacterium tuberculosis genotype clusters, Washington, USA. Variables are classified as follows: a = number of tuberculosis (TB) cases with the genotype of interest in the selected county; b = number of cases with the genotype of interest in the United States; c = number of cases without the genotype of interest in the selected county; d = number of cases without the genotype of interest in the United States; N = total number of TB cases.

Qualitative analysis came from a 5-member expert panel of TB public health officials in Washington. In 2008, the panel participated in a discussion of all county-level TB clusters, ranking each as high or low priority for additional investigation. Priority was determined on the basis of a review of patient characteristics, epidemiologic links from field investigations, and maps of genotype distributions. The panel also had information from enhanced contact investigations from local public health investigation teams that included the ability to order IS6110 restriction fragment-length polymorphism (IS6110 RFLP) and 24-locus MIRU-VNTR testing for clusters of concern, but results from these tests were not universally available. The ranking exercise with the same 5-member panel was repeated after period 1 for clusters from period 1. The expert panel was blinded to the LLR.

LLRs were compared with the expert opinion ranking to assess concordance. With expert opinion as the standard, negative and positive predictive values (NPV and PPV, respectively) were calculated for period 1 using a cutoff point of LLR >5.0. Alternative cutoff points were evaluated to maximize NPV and PPV. Sensitivity and specificity of the >5.0 LLR cutoff point and exact binomial 95% CIs were calculated for period 1 clusters. An alternative cutoff point to maximize sensitivity and specificity was also determined.

A total of 806 TB cases were reported in Washington during period 1. Of 659 culture-positive cases, 642 (97.4%) had genotyped isolates; of these, 318 cases formed 21 clusters. Five of these clusters had a high LLR; the expert panel ranked all 5 of these clusters high priority and identified them as clusters of concern. Of the 16 clusters with LLR <5.0, the expert panel ranked 12 (75.0%) as low priority (Table).

Comparison of geospatial analysis results and expert panel priority status rankings for county-level genotype clusters of TB cases, Washington, USA, 2005–2010*
CountySpoliogtype12-locus MIRU-VNTRPeriod 1†
Period 2†
Key epidemiologic features‡
No. casesLLRExpert priorityNo. casesLLRExpert priority
A0000000000037712233211536433231.8High1918.3HighH, SA
000000000003771223325173533170.1High80.4LowFB
677777477413771254326223432140.2Low140.1LowFB
70337740000177122742511343441.2Low62.9LowFB, H
00000000000377122332515353330.5Low1<0.01LowFB
00000000000377122332516353330.2Low71.7LowFB, SA
67777747741377125432622342261.3Low40.4LowFB
77700037776077122512511332230.6Low31.0LowFB
00000000000377122232517354392.4High30.1HighH, SA
77777777776077122512511332230.7Low20.4LowFB
00000000000377122232517353372.6High82.3High
00000000000377122332517343331.5Low10.2LowFB
00000000000377122332513353331.4Low0<0.01LowFB
00000000000000022332512353449.2High0<0.01LowH, SA

777776777620601
224325153323
3
4.4
High

0
0.03
Low
H, SA
B
000000000003771
223425173563
13
19.5
High

8
13
High
H, SA
C67777747741377125432622343261.3Low3<0.01LowFB, SA

777776777760771
125325153225
7
9.8
High

2
1.8
High
H, SA
D67777747741377125432622343250.8Low40.4LowFB

777776757760771
223325143324
6
8.7
High

1
1.2
Low

E67777747741377125432622343231.0Low20.9LowFB

*LLRs >5.0 (boldface) indicate greater possibility that the cluster represents geographically concentrated TB cases, a proxy for recent TB transmission). TB, tuberculosis; MIRU-VNTR, mycobacterial interspersed repetitive unit–variable number tandem repeat; LLR, log-likelihood ratio, a measure of geospatial concentration; H, homelessness; FB, foreign born; SA, substance abuse. 
†Period 1, January 1, 2005–December 31, 2007; period 2, January 1, 2008–December 31, 2010.
‡Blank cells indicate none present.

A total of 723 TB cases were reported in Washington during period 2. Of 592 culture-positive cases, 576 (97.3%) had genotyped isolates. The expert panel reexamined the 21 clusters identified during period 1 and focused on new activity within those areas during period 2. Two clusters with a high LLR during period 1 continued to have a high LLR during period 2; the expert panel continued to rank these high priority. The 3 other clusters that had a high LLR during period 1 had a low LLR for period 2; one of those was still considered a high priority by the panel. Of the remaining 16 clusters, which continued to have a LLR <5, 14 (87.5%) were still ranked low priority by the panel.

Two clusters in the same county, PCR00309 and PCR00803, had low LLRs but were considered high priority by the expert panel. For cluster PCR00309, the panel cited high levels of homelessness among case-patients as reason to rank it high priority. For cluster PCR00803, the panel cited a highly mobile population from a TB-endemic country that regularly traveled into and out of the United States as reason to rank it high priority. However, the travel history among case-patients in this cluster made it difficult for investigators to determine whether transmission was occurring within Washington or abroad.

The NPV and PPV for a LLR cutoff point of 5.0 were 75% and 100%, respectively. Lowering the cutoff point to a LLR >2.0 increased the NPV to 92.3%, but the PPV remained at 100%. For period 1 clusters only, a cutoff point of LLR >5.0 generated a sensitivity of 55.6% (95% CI 21.2%–86.3%) and specificity of 100% (95% CI 73.5%–100.0%) for identifying clusters for further investigation. Decreasing the cutoff point to >2.0 increased the sensitivity to 88.9% (95% CI 51.8%–99.7%) but did not change the specificity (100%; 95% CI 73.5%–100.0%).

Conclusions

The geospatial statistic in this study was highly correlated with experts’ perceived need for public health action. This finding indicates that automated alerts generated on the basis of geospatial concentration of TB cases might help the state TB program identify clusters that would benefit from additional investigation. Automated alerts can be generated by using routinely collected surveillance data and are currently part of the national TB Genotyping Information Management System (10).

Patient and contact characteristics, transmission venues, and temporality all contribute toward prioritization determination. For example, during period 1, a total of 6 (66.6%) clusters ranked high priority by the expert panel were characterized by homelessness or substance abuse among case-patients, and 8 (88.9%) were characterized by US-born case-patients (Table).

Conversely, 11 (91.7%) clusters ranked low priority were characterized by case-patients who were foreign-born, a known risk factor for latent TB infection (7). None of the period 1 clusters with LLR >5 and only 1 of 9 clusters ranked as high priority by the expert panel were characterized by foreign-born case-patients. These results indicate the need for further study to identify the limitations of the LLR score in detecting localized and recent TB transmission among foreign-born case-patients.

The availability of IS6110 RFLP or 24-locus MIRU-VNTR testing results to the expert panel is the current standard for fieldwork and could have introduced an information bias for the panel in this study. Although this effect is unknown, lack of universal IS6110 RFLP and 24-locus MIRU-VNTR test results is a limitation of this study.

We found that geospatial statistics based on TB genotyping and surveillance data could help identify and prioritize likely recent disease transmission events in Washington. In addition, LLR values should be incorporated into ongoing evaluation by the expert panel; in fact, LLR is now included in routine genotype and cluster reviews. Geospatial statistics are an attractive approach to prioritization, but additional field-based research is needed to assess whether factors such as epidemiologic characteristics could be used to further develop a prioritization algorithm. Integrating these factors and determining ideal cutoff points in different settings will increase predictive value.

Suggested citation for this article: Lindquist S, Allen S, Field K, Ghosh S, Haddad MB, Narita M, et al. Prioritizing tuberculosis clusters by genotype for public health action, Washington, USA. Emerg Infect Dis [Internet]. 2013 Mar [date cited]. http://dx.doi.org/10.3201/eid1903.121453

Acknowledgments

We thank Steve Kammerer, Heidi Soeters, and Juliana Grant for help with conducting the original expert panel and with preparation of this manuscript.

Dr Lindquist is a tuberculosis medical consultant for the Washington State Department of Health, a health officer for the Kitsap County Health District, and a pediatrician at the Port Gamble S’Klallam Tribal Medical Clinic. His research interests include epidemiology of tuberculosis, interferon gamma release assays and public health, and public health delivery systems.

ReferencesCastro KG, Dooley SW, Curran JW. Transmission of HIV-associated tuberculosis to health-care workers. Lancet. 1992;340:10434 10.1016/0140-6736(92)93063-S1357437Centers for Disease Control and Prevention New CDC program for rapid genotyping of Mycobacterium tuberculosis isolates. MMWR Morb Mortal Wkly Rep. 2005;54:47 .16177693Barnes PF, Cave MD. Molecular epidemiology of tuberculosis. N Engl J Med. 2003;349:114956 10.1056/NEJMra02196413679530Buff AM, Moonan PK, Desai MA, McKenna TL, Harris DA, Rogers BJ, South Carolina tuberculosis genotype cluster investigation: a tale of substance abuse and recurrent disease. Int J Tuberc Lung Dis. 2010;14:13479 .20843429Pevzner ES, Robison S, Donovan J, Allis D, Spitters C, Friedman R, Tuberculosis transmission and use of methamphetamines in Snohomish County, WA, 1991–2006. Am J Public Health. 2010;100:24816 10.2105/AJPH.2009.16238820167896Driver CR, Kreiswirth B, Macaraig M, Clark C, Munsiff SS, Driscoll J, Molecular epidemiology of tuberculosis after declining incidence, New York City, 2001–2003. Epidemiol Infect. 2007;135:63443 10.1017/S095026880600727817064454Moonan PK, Ghosh S, Kammerer JS, Oeltmann JE, Cowan LS, Navin TR. Estimating recent transmission of Mycobacterium tuberculosis in the United States. Emerg Infect Dis. 2012;18:45865 10.3201/eid1803.11110722377473Centers for Disease Control and Prevention Reported tuberculosis in the United States, 2010. Atlanta: The Centers; 2011Cowan LS, Diem L, Monson T, Wand P, Temporado D, Oemig TV, Evaluation of a two–step approach for large-scale, prospective genotyping of Mycobacterium tuberculosis isolates in the United States. J Clin Microbiol. 2005;43:68895 10.1128/JCM.43.2.688-695.200515695665Ghosh S, Moonan PK, Cowan L, Grant J, Kammerer S, Navin T. Tuberculosis Genotyping Information Management System: enhancing tuberculosis surveillance in the United States. Infect Genet Evol. 2012;12:7828 10.1016/j.meegid.2011.10.01322044522