PLoS OnePLoS ONEplosplosonePLoS ONE1932-6203Public Library of ScienceSan Francisco, USA243491823861328PONE-D-13-0680010.1371/journal.pone.0082048Research ArticleGeotemporal Analysis of Neisseria meningitidis Clones in the United States: 2000–2005Geotemporal Analysis of N. meningitidisWiringaAnn E.1ShuttKathleen A.2MarshJane W.2CohnAmanda C.3MessonnierNancy E.3ZanskyShelley M.4PetitSusan5FarleyMonica M.6GershmanKen7LynfieldRuth8ReingoldArthur9SchaffnerWilliam10ThompsonJamie11BrownShawn T.12LeeBruce Y.13HarrisonLee H.214*Department of Epidemiology, University of Pittsburgh Graduate School of Public Health, Pittsburgh, Pennsylvania, United States of AmericaInfectious Diseases Epidemiology Research Unit, University of Pittsburgh Graduate School of Public Health and School of Medicine, Pittsburgh, Pennsylvania, United States of AmericaMeningitis and Vaccine Preventable Diseases Branch, Division of Bacterial Diseases, Centers for Disease Control and Prevention, Atlanta, Georgia, United States of AmericaNew York State Department of Health, Albany, New York, United States of AmericaConnecticut Department of Public Health, Hartford, Connecticut, United States of AmericaEmory University and VA Medical Center, Atlanta, Georgia, United States of AmericaColorado Department of Public Health and Environment, Denver, Colorado, United States of AmericaMinnesota Department of Health, St. Paul, Minnesota, United States of AmericaUniversity of California, Berkeley, Berkeley, California, United States of AmericaVanderbilt University School of Medicine, Nashville, Tennessee, United States of AmericaOregon Public Health Division, Portland, Oregon, United States of AmericaPittsburgh Supercomputing Center, Carnegie Mellon University, Pittsburgh, Pennsylvania, United States of AmericaPublic Health Computational and Operations Research (PHICOR), University of Pittsburgh School of Medicine and Graduate School of Public Health, Pittsburgh, Pennsylvania, United States of AmericaDepartment of International Health, Johns Hopkins Bloomberg School of Public Health, Baltimore, Maryland, United States of AmericaBereswillStefanEditorCharité-University Medicine Berlin, Germany* E-mail: lharriso@edc.pitt.edu

Competing Interests: This study was funded in part by the Centers for Disease Control and Prevention. Dr. Harrison has received research support and lecture fees from Sanofi Pasteur; lecture fees from Novartis Vaccines; and has served as a consultant to GlaxoSmithKline, Merck, Novartis Vaccines, Sanofi Pasteur, and Pfizer. Dr. Harrison’s financial ties with industry were terminated before he became a voting member of the Advisory Committee on Immunization Practices in July 2012. Dr. Schaffner is a member of data safety monitoring boards for Merck and has served as a consultant to Dynavax, Pfizer, GSK and Sanofi Pasteur. Dr. Marsh has received research support from Sanofi Pasteur, Merck and ViroPharma Inc. Martin Kulldorff provided technical assistance with the SaTScan software. There are no patents, products in development or marketed products to declare. This does not alter our adherence to all the PLOS ONE policies on sharing data and materials.

Conceived and designed the experiments: AW KS JM AC NM BL LH. Performed the experiments: AW KS. Analyzed the data: AW KS. Contributed reagents/materials/analysis tools: SZ SP MF KG RL AR WS JT SB. Wrote the paper: AW KS AC NM LH.

201312122013312014812e82048222013291020132013This is an open-access article, free of all copyright, and may be freely reproduced, distributed, transmitted, modified, built upon, or otherwise used by anyone for any lawful purpose. The work is made available under the Creative Commons CC0 public domain dedication.Background

The detection of meningococcal outbreaks relies on serogrouping and epidemiologic definitions. Advances in molecular epidemiology have improved the ability to distinguish unique Neisseria meningitidis strains, enabling the classification of isolates into clones. Around 98% of meningococcal cases in the United States are believed to be sporadic.

Methods

Meningococcal isolates from 9 Active Bacterial Core surveillance sites throughout the United States from 2000 through 2005 were classified according to serogroup, multilocus sequence typing, and outer membrane protein (porA, porB, and fetA) genotyping. Clones were defined as isolates that were indistinguishable according to this characterization. Case data were aggregated to the census tract level and all non-singleton clones were assessed for non-random spatial and temporal clustering using retrospective space-time analyses with a discrete Poisson probability model.

Results

Among 1,062 geocoded cases with available isolates, 438 unique clones were identified, 78 of which had ≥2 isolates. 702 cases were attributable to non-singleton clones, accounting for 66.0% of all geocoded cases. 32 statistically significant clusters comprised of 107 cases (10.1% of all geocoded cases) were identified. Clusters had the following attributes: included 2 to 11 cases; 1 day to 33 months duration; radius of 0 to 61.7 km; and attack rate of 0.7 to 57.8 cases per 100,000 population. Serogroups represented among the clusters were: B (n = 12 clusters, 45 cases), C (n = 11 clusters, 27 cases), and Y (n = 9 clusters, 35 cases); 20 clusters (62.5%) were caused by serogroups represented in meningococcal vaccines that are commercially available in the United States.

Conclusions

Around 10% of meningococcal disease cases in the U.S. could be assigned to a geotemporal cluster. Molecular characterization of isolates, combined with geotemporal analysis, is a useful tool for understanding the spread of virulent meningococcal clones and patterns of transmission in populations.

This study was supported by the US Centers for Disease Control and Prevention grant DHM 99–006042. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. CDC coauthors participated in various aspects of the study and manuscript reparation.
Introduction

Neisseria meningitidis is an important cause of meningitis and other serious bacterial infections globally. [1] In the U.S., over 98% of meningococcal disease cases are considered to be sporadic, that is, unrelated to other cases, but outbreaks also occur. [2], [3] Identification of meningococcal outbreaks typically relies on serogrouping and epidemiologic definitions. [3], [4] The incidence of meningococcal disease in the U.S. is at historically low levels [2].

Advances in objective, DNA sequence-based molecular epidemiologic tools for N. meningitidis have enhanced the ability to characterize this organism. Multilocus sequence typing (MLST) is a standard molecular subtyping approach for determining genetic lineage. [5] DNA sequencing of genes that encode outer membrane proteins (OMPs) provides additional discriminatory power among strains belonging to the same sequence type (STs). [6], [7] Together, MLST and OMP genotyping allow for the classification of meningococcal isolates into specific clones, which can be used to detect outbreaks and study patterns of spread within populations. [8][10] In one study a spatial scan statistic was used to identify 26 clusters of invasive meningococcal disease in Germany using a clone definition based on serogroup, porA and fetA typing [8].

We recently reported the population structure of invasive meningococcal isolates throughout the United States. [11] The purpose of the present study was to assess geotemporal clustering patterns of specific meningococcal clones among the isolates reported in that study to determine whether this approach could identify both known and previously undetected clusters of meningococcal disease.

MethodsEthics Statement

The study was approved by the Johns Hopkins Bloomberg School of Public Health, University of Pittsburgh, Emory University, Georgia Department of Public Health, Vanderbilt University School of Medicine and Tennessee Department of Health institutional review boards. Institutional review board approval was not required from the California, Colorado, Connecticut, Minnesota, New York or Oregon Active Bacterial Core surveillance sites because meningococcal disease is a reportable condition and the activities fall under routine disease surveillance authority.

Study Isolates and Determination of Serogroup

Invasive meningococcal study isolates were obtained through active surveillance during the period of January 1, 2000-December 31, 2005 from 9 Active Bacterial Core surveillance (ABCs) sites. ABCs, an active laboratory- and population-based surveillance program for invasive bacterial pathogens, is a core component of the Centers for Disease Control and Prevention (CDC) Emerging Infections Programs Network. [12] ABCs defines a case as isolation of N. meningitidis from a normally sterile site, such as blood or cerebrospinal fluid, in a resident of an ABCs surveillance area. [13], [14] The CDC case definition of a serogroup C N. meningitidis outbreak is ≥3 confirmed or probable cases in ≤3 months, resulting in a primary attack rate of ≥10 cases per 100,000 population among persons with a common organizational affiliation or who live in the same community. [4] Although the Advisory Committee on Immunization Practices (ACIP) definitions and guidelines for the control and prevention of meningococcal disease were initially developed for serogroup C meningitis, the same principles are relevant for the control of cases attributable to other vaccine-preventable N. meningitidis serogroups including A, Y and W-135 [4].

Participating ABCs sites included the following areas: California (three counties in the San Francisco Bay area), Colorado (5 counties in the Denver area), Connecticut, Georgia, Maryland, Minnesota, New York (7 counties in the Rochester area and 8 in the vicinity of Albany), Oregon, and Tennessee (11 urban counties in the Nashville, Memphis, Knoxville and Chattanooga areas). The population under surveillance in 2005 was approximately 38 million persons [15].

Laboratory work for this study was performed at the CDC and the University of Pittsburgh. Serogrouping, MLST, and OMP genotyping of porA VR1 and VR2, porB, and fetA VR were performed as previously described. [7], [11] Meningococcal clones were classified according to serogroup, ST, and OMP (porA, porB, and fetA) genotyping. To be considered as belonging to the same clone, isolates had to be indistinguishable based on the results of these assays.

Geotemporal Cluster Analysis

The methodology of cluster detection using SaTScan version 9.1.1 (Information Management Services, Inc., Silver Spring, MD, and Martin Kulldorff, Boston, MA) is available on the Technical Documentation page of the developer’s web site (www.satscan.org/techdoc.html). [16][21] All analyses were conducted using a dedicated personal computer with 4 GB of RAM. Case data were aggregated to the census tract level and all non-singleton clones were assessed for geotemporal clustering using retrospective space-time analyses with a discrete Poisson probability model. Analyses tested the hypothesis that the expected number of cases in each census tract was proportional to its population size. The scan identified clusters with high rates of N. meningitidis indicative of non-random spatial and/or temporal distribution of cases.

Date of isolate culture was used to define the onset date of illness. The units of time precision and time aggregation for all analyses were days and months, respectively. The maximum temporal window size for all analyses was 50% of the study period (36 months), and the maximum spatial window 50% of the population in each ABCs site. We chose these wide windows to gain insight into the persistence of invasive meningococcal clones over space and time. Our study aimed to assess the geotemporal distribution of molecularly related clusters and, accordingly, the selection of this long temporal window allowed identification of clusters demonstrating a persistence of specific clones with durations exceeding those of classically defined meningococcal outbreaks. When the maximum window is set to 50% of the population, both small and large clusters can be identified. The statistical significance of each cluster was determined by 999 replications of Monte Carlo hypothesis testing and interpreted as significant when p≤0.05. Nine ABCs sites were analyzed independently. To detect possible clustering across contiguous ABCs sites, the following combined locations were also analyzed: (1) the census tracts comprising metropolitan Chattanooga, TN, and the state of Georgia; (2) the census tracts in the vicinity of Albany, NY, and the state of Connecticut. ABCs sites were asked to determine whether the clusters reported herein were identified at the time they occurred, and whether the current study failed to detect any previously known clusters.

Population and Geographic Data

Census 2000 geographic boundary files for states, counties and census tracts participating in ABCs were obtained from the Census Cartographic Boundary Files Collection of the U.S. Census Bureau’s Geography Division. [22] Population and land area data for ABCs surveillance sites were retrieved from the U.S. Census Bureau Census 2000 Summary File (SF 1) 100-Percent Data Set. [23] Using ArcGIS version 9.3.1 (ESRI, Inc., Redlands, CA), the location of each case was geocoded using a perturbation algorithm. The perturbation distance was designed to be an inverse function of the population density in the respective census tract. The purpose of the perturbation is to de-identify the location of individual cases on a map. As a result, geocoded case locations may exhibit displacement into an adjacent census tract in figures. Analyses were conducted at the level of census tract with cases aggregated to their true tract of residence. Census tract population density (residents per km2) and attack rate (per 100,000 population) were calculated based on tract population and case count.

Results

A total of 1,159 N. meningitidis isolates were obtained from cases meeting the ABCs case definition between January 1, 2000 and December 31, 2005. 1,062 (91.7%) of these cases were successfully geocoded and included. The case counts by ABCs site were California, 113 (10.6% of the total); Colorado, 49 (4.6%); Connecticut, 82 (7.7%); Georgia 163 (15.4%); Maryland, 117 (11.0%); Minnesota, 126 (11.9%); New York, 75 (7.1%); Oregon, 267 (25.1%); and Tennessee, 70 (6.6%). Isolates were characterized by molecular subtyping into 438 unique clones, 78 (17.8%) of which were non-singletons. Overall there were 702 cases attributable to non-singleton clones, accounting for 66.0% of all geocoded isolates. Molecular characteristics of the non-singletons, ranked in decreasing order of frequency, are shown in Table 1.

10.1371/journal.pone.0082048.t001Case counts of non-singleton meningococcal clones, ranked by decreasing frequency.
Clone IDSerogroupClone (CC: ST: porB: porA VR1, porA VR2: fetA)Count (% of total)ABCs site(s) with clone
1Y23∶23∶3-36: P1.5-2,10-1: F4-1142 (13.4)CA,CO,CT,GA,MD,MN,NY,OR,TN
2B32∶32∶3-24: P1.7,16: F3-3113 (10.6)CA,CO,CT,GA,MD,MN,OR,TN
3Y23∶23∶2-55: P1.5-1,2-2: F5-859 (5.6)CA,CO,CT,GA,MD,MN,NY,OR,TN
4C11∶11∶2-2: P1.5,2: F3-640 (3.8)CA,CO,CT,GA,MD,MN,NY,OR,TN
5B32∶32∶3-24: P1.7,16-20: F3-327 (2.5)CA,CO,GA,MD,MN,OR
6C11∶11∶2-2: P1.5,2: F1-3022 (2.1)CO,CT,GA,MD,MN,OR,TN
7B162∶162∶3-73: P1.22,14: F5-916 (1.5)CA,GA,MD,MN,NY,TN
8B32∶32∶3-36: P1.7,16: F3-314 (1.3)OR
9C11∶2962∶2-75: P1.5-1,10-4: F3-614 (1.3)CA,GA,MD,NY,OR,TN
10C11∶11∶2-2: P1.22-1,14: F3-613 (1.2)GA,MD,MN,NY
11B32∶3584∶3-1: P1.7,16: F3-310 (0.9)OR
12C103∶2006∶2-110: P1.5-1,10-4: F3-910 (0.9)MN
13C11∶11∶2-2: P1.5-1,10-8: F3-610 (0.9)GA,MD,MN,NY
14B32∶32∶3-1: P1.7,16: F3-39 (0.8)CT,GA,MD,NY,OR
15C11∶2962∶2-88: P1.5-1,10-4: F3-69 (0.8)CO,GA,MD,MN,OR,TN
16B32∶32∶3-84: P1.7,16: F3-38 (0.8)GA,TN
17Y23∶1625∶2-55: P1.5-1,2-2: F5-87 (0.7)CT,GA,MD,MN,OR,TN
18Y23∶23∶3-36: P1.5,2: F4-17 (0.7)CA,CT,MN,OR,TN
19B41/44∶136∶3-107: P1.17,16-3: F5-56 (0.6)CA,CO,GA,NY,TN
20C11∶2961∶2-48: P1.5,2: F1-306 (0.6)MD,NY
21B35∶35∶3-39: P1.22-1,14: F4-15 (0.5)GA,MD,OR
22B41/44∶154∶3-1: P1.7-2,4: F1-55 (0.5)CA,CT,OR
23C11∶11∶2-73: P1.5,2: F3-65 (0.5)GA,MD,MN,NY
24C32∶32∶3-24: P1.7,16: F3-35 (0.5)MN,OR
25W-13522∶22∶2-23: P1.18-1,3: F4-15 (0.5)CO,MN,OR,TN
26Y167∶1624∶2-55: P1.5-1,10-4: F3-45 (0.5)CT,GA,NY,OR
27B32∶32∶3-133: P1.7,16-20: F3-34 (0.4)OR
28B41/44∶4682∶3-71: P1.22-1,14: F5-24 (0.4)CO,OR
29C11∶11∶2-2: P1.5,2: F1-54 (0.4)CA
30C11∶11∶2-2: P1.5-1,2-2: F3-64 (0.4)GA,NY,OR
31Y23∶1625∶2-141: P1.5-1,2-2: F5-84 (0.4)OR
32Y23∶3582∶2-55: P1.5-1,2-2: F5-84 (0.4)OR
33B32∶32∶3-1: P1.7,16-20: F3-33 (0.3)CT,MD
34B41/44∶136∶3-107: P1.17,16-23: F5-53 (0.3)CA,GA,TN
35B41/44∶170∶3-138: P1.21,16: F1-53 (0.3)OR
36B41/44∶42∶3-1: P1.7-2,4: F1-53 (0.3)CA
37B41/44∶44∶3-38: P1.7-1,1: F1-73 (0.3)NY,OR
38B41/44∶44∶3-45: P1.7-4,1: F1-73 (0.3)OR
39B41/44∶5097∶3-1: P1.7-2,4: F1-53 (0.3)TN
40B41/44∶5111∶3-45: P1.21,16: F1-73 (0.3)CT
41CNo CC:2048∶3-16: P1.5,2: F3-63 (0.3)GA
42C11∶11∶2-2: P1.17,16-3: F3-63 (0.3)MD
43C11∶11∶2-2: P1.5,2: F3-33 (0.3)CA,OR
44C11∶11∶2-2: P1.5,2: F5-53 (0.3)OR
45C11∶11∶2-2: P1.5-1,2: F5-363 (0.3)MN
46Y23∶1621∶3-36: P1.5-2,10-1: F4-13 (0.3)MD
47BNo CC:2048∶3-16: P1.12-1,16-8: F3-62 (0.2)CO,OR
48BNo CC:2875∶2-136: porA Deletion: F4-12 (0.2)GA
49B254∶254∶3-223: P1.19,15: F1-52 (0.2)GA
50B269∶2974∶3-113: P1.7-2,13-1: F5-72 (0.2)MD
51B32∶1364∶3-24: P1.7,16: F3-32 (0.2)CA,OR
52B32∶32∶3-107: P1.22-1,14: F3-32 (0.2)GA
53B32∶32∶3-1: P1.7,16-33: F3-32 (0.2)NY,OR
54B32∶32∶3-24: P1.7,16: F1-72 (0.2)OR
55B41/44∶318∶3-1: P1.7-2,4: F1-72 (0.2)MN
56B41/44∶409∶3-82: P1.18-1,34-2: F1-52 (0.2)CA,GA
57B41/44∶41∶3-172: P1.7-2,4: F1-52 (0.2)GA,TN
58B41/44∶41∶3-1: P1.7-2,4: F1-52 (0.2)CT
59B41/44∶437∶3-114: P1.22-1,14: F5-22 (0.2)MD
60B41/44∶437∶3-71: P1.22-1,14: F5-22 (0.2)GA
61B41/44∶43∶3-16: P1.19,15-1: F1-52 (0.2)GA,NY
62B41/44∶4489∶3-1: P1.7-2,4: F1-52 (0.2)GA
63B60∶60∶3-8: P1.22-1,14: F3-92 (0.2)MD
64C103∶5837∶2-22: P1.17,16-3: F1-182 (0.2)TN
65C11∶11∶2-2: P1.5,2: F3-12 (0.2)MD
66C11∶11∶2-2: P1.5,2: F4-122 (0.2)CA,OR
67C11∶11∶2-2: P1.5-2,10-2: F5-362 (0.2)MN
68C11∶11∶2-85: P1.22-1,14: F3-62 (0.2)MN
69C41/44∶41∶3-1: P1.7-2,4: F1-52 (0.2)CA
70C8∶8∶2-3: P1.5,2: F5-82 (0.2)MN
71W-13522∶1476∶2-23: P1.18-1,3: F4-12 (0.2)CO,OR
72W-13522∶22∶2-109: P1.18-1,3: F4-12 (0.2)CO
73Y22∶1265∶2-23: P1.18-1,3: F1-72 (0.2)GA
74Y23∶183∶3-53: P1.5-2,10-2: F4-12 (0.2)NY
75Y23∶23∶3-36: P1.5-2,10-12: F4-12 (0.2)MD,TN
76Y23∶23∶3-36: P1.5-2,10-29: F4-12 (0.2)GA,NY
77Y23∶23∶3-53: P1.5-2,10-2: F4-12 (0.2)GA,NY
78Y23∶3587∶3-36: P1.5-2,10-2: F4-12 (0.2)GA,OR

ST, sequence type; CC, clonal complex; VR, variable region; CA, California; CO, Colorado; CT, Connecticut; GA, Georgia; MD, Maryland; MN, Minnesota; NY, New York; OR, Oregon; TN, Tennessee.

Given the large number of distinct non-singleton clones identified (116 distinct site/clone parings, with 7 separate analysis files for each–one for each of the 6 study years, and a composite 2000–2005 file), parameterizing and executing the scans to generate the results reported herein was time and computationally intensive. Individual scans took anywhere from seconds to upwards of 24 hours.

Thirty-two statistically significant clusters involving 107 cases (10.1% of all geocoded isolates) attributable to 23 distinct clones were identified (Table 2 and Figure 1) by independent analysis of each ABCs site. Clusters were identified in all sites except Colorado. Clusters ranged in duration from 1 day to 33 months. 19 clusters (59.4%) were composed of two cases and 13 (40.6%) included 3 or more cases. Incidence ranged from 0.7 to 57.8 per 100,000, with 10 clusters (31.3%) having an attack rate of ≥10 cases per 100,000 population over the cluster time period. Annualized incidence ranged from 1.1 to 693.6 cases per 100,000 population. The range of cluster radii was 0 to 61.7 km. Five clusters (15.6%) had radius of 0 km, indicating that the involved cases occurred within a single census tract. Three of the five zero-radius clusters (Clusters C, G and M) were previously identified by the respective ABCs sites, and all were noted to be case pairs among household members. Two zero-radius clusters (Clusters S and CC), also comprised of case pairs, were previously unidentified. No purely temporal clusters (those encompassing an entire ABCs surveillance area) or purely spatial clusters (those with a temporal window spanning the entire 72 month study period) were identified.

10.1371/journal.pone.0082048.g001Meningococcal clusters, by Active Bacterial Core surveillance site, 2000–2005.

(A) California with San Francisco Bay area inset; (B) Connecticut; (C) Georgia; (D) Maryland; (E) Minnesota; (F) New York; (G) Oregon; (H) Tennessee; (I) Georgia and Tennessee.

10.1371/journal.pone.0082048.t002Characteristics of statistically significant meningococcal disease clusters.
Cluster IDClone IDASerogroupSiteFirst CaseLast CaseRadius (km)PopulationNo. CasesAttack Rate (per 100,000)Expected CasesObserved/ExpectedAnnualized Incidence (per 100,000)p-value
A1YCA1/9/20003/30/20002.576,37745.20.021190.520.90.003
B3YCT11/26/20014/17/20024.5105,97343.80.04981.69.10.027
C3YCT2/19/20053/15/20050.0B3561256.20.000533773.6674.00.022
D1YGA2/27/20009/15/200138.8399,87892.30.5317.01.40.003
E16BGA1/22/20057/7/200519.395,67744.20.0045888.98.40.001
F48BGA11/24/20009/16/200114.320,51129.80.000762631.611.70.014
G62BGA6/11/20026/11/20020.0B3956250.60.000013153,846.2606.70.001
H1YMD12/21/20002/28/20018.5361,66951.40.06280.65.50.004
I10CMD3/26/20016/19/200131.7262,09620.80.0055363.63.10.049
J46YMD3/6/20003/7/20003.9169,84521.20.00141428.614.10.023
K59BMD2/8/20014/2/20011.89588220.90.0001513,333.3125.20.003
L4CMN8/10/20008/19/20005.5111,01021.80.00131538.521.60.038
M13CMN10/30/200211/1/20020.0B3460257.80.0001216,666.7693.60.004
N67CMN8/21/20023/15/20034.416,616212.00.000752666.720.60.011
O70CMN3/31/20024/17/20024.423,59728.50.000277407.4101.70.005
P9CNY12/3/200212/11/200220.4138,88421.40.00191052.617.30.010
Q10CNY2/4/20018/1/200127.9382,31941.00.07057.12.10.007
R20CNY11/3/20015/14/20028.37818225.60.00181111.143.90.027
S26YNY4/8/200012/12/20020.0B8549223.40.0037540.58.50.017
T37BNY9/25/20012/9/20045.210,547219.00.0042476.27.80.022
U1YOR8/10/200011/2/200046.9717,34450.70.1533.32.80.033
V2BOR12/3/20004/18/200125.8300,817103.30.5916.98.00.001
W4COR9/21/200210/11/200261.733,85138.90.000833614.5106.30.001
X8BOR1/24/200012/31/200130.6138,831117.90.1957.94.00.001
Y14BOR2/13/20012/23/200143.492,92622.20.00141428.625.80.025
Z24COR12/30/20041/2/20052.320,76729.60.000345882.4115.60.006
AA35BOR1/25/20007/31/200056.290,87133.30.0077389.65.70.003
BB54BOR3/10/20045/3/20044.649,12124.10.00121666.724.40.012
CC1YTN6/28/20019/16/20010.0B3863251.80.00111818.2207.10.032
DD9CTN9/25/20018/22/200211.0397,45641.00.09442.61.10.012
EE16BTN10/9/200510/20/20056.12467028.10.000504000.097.30.011
FF39BTN4/22/20033/2/20056.625,579311.70.0091329.76.10.004
GG16BGA/TN1/22/200510/20/200521.9341,09272.10.045155.62.70.001

Table 1 for detailed information on the molecular characteristics of each clone.A See

B Case data were aggregated to the level of census tract for these analyses. A radius of 0 km indicates that all cases in that cluster occurred within the same census tract.

In analyses of contiguous ABCs sites, one cluster (cluster GG) caused by clone 16 included 7 cases in metropolitan Chattanooga, TN and northwestern Georgia (radius = 21.9 km) that occurred between January 22, 2005 and October 20, 2005 (Figure 1). In single-site analyses, a 4-case cluster of Clone 16 (cluster E, p = 0.001, radius = 19.3 km) was detected in Georgia and a 2-case cluster in Tennessee (cluster EE, p = 0.011, radius = 6.1 km). The cross-site analysis captured an additional case in Tennessee not associated with the independently identified cluster in that state.

Eleven clusters (34.4% of all clusters) had previously been identified by the ABCs sites, while 21 clusters (65.6%) had not been identified. Cluster GG, which spanned the Georgia-Tennessee border, was previously identified by local public health practitioners. A cluster comprised of 2 cases occurring 1 week apart among students at the same university was reported by Connecticut but not detected in our analyses. These isolates were identical by OMP genotyping but differed at a single MLST locus (ST-1374 versus ST-40).

Three serogroups were represented among the clusters: B (n = 12 clusters, 45 cases), C (n = 11 clusters, 27 cases) and Y (n = 9 clusters, 35 cases), indicating that 20 clusters (62.5%) were caused by serogroups represented in meningococcal vaccines that are available in the U.S. Although seven cases attributable to two distinct serogroup W-135 clones were reported, no significant serogroup W-135 N. meningitidis clusters were identified. Eleven distinct serogroup B clones were responsible for 12 significant clusters (37.5% of all clusters) in 5 states: Georgia (3), Maryland (1), New York (1), Oregon (5), Tennessee (2), and an additional cross-site cluster in Georgia/Tennessee. No cluster attributable to a vaccine-preventable serogroup met the CDC definition of an outbreak warranting consideration of vaccination for disease control, and our analyses did not fail to detect any clusters that prompted local health officials to consider vaccination as a control strategy.

Sensitivity analyses were conducted to assess the impact of alternate scanning window parameterizations. Spatial and temporal windows were varied from 25–50% and 25–90%, respectively. The clusters detected differed negligibly for each permutation of settings, so we elected to present the results for an intermediate set of parameterizations with spatial and temporal windows each set at 50%.

Discussion

We identified 32 meningococcal case clusters with non-random spatial and temporal distribution, 21 (65.6%) of which had not been previously identified. To our knowledge, this is the first study of the geotemporal distribution of meningococcal clones causing invasive disease in the United States. Our analyses highlight the large diversity of circulating meningococcal clones in the U.S. and the fact that most invasive meningococcal disease is caused by a limited number of MLST-defined meningococcal lineages. Importantly, 10.1% of cases belonged to a molecularly identified cluster, which is higher than the proportion of cases identified in outbreaks using traditional epidemiologic methods. [3] All five clusters with a radius of 0 km (indicating that both cases in the pair occurred in the same census tract) were associated with case pairs. The attack rate and projected annualized incidence rate in the affected census tracts ranged from 23.4 to 57.8 and 8.5 to 693.6 cases per 100,000 population, respectively. This highly non-random spatial and temporal distribution of cases underscores the importance of assessing the geotemporal distribution of as few as 2 isolates of the same clone. We did not detect any serogroup C clusters that met the CDC definition of an outbreak warranting consideration of vaccination. Assuming the ACIP serogroup C outbreak definition thresholds, no clusters of serogroup Y disease were detected that would prompt consideration of vaccination for disease control. [4] This lack of geographically expansive clusters comprised of large numbers of cases may be, in part, a reflection of the historically low incidence of meningococcal disease in the United States. [2] During the study period the national projection for the annual incidence of meningococcal disease in the U.S. decreased from 0.8 cases per 100,000 population (year 2000) to 0.35 cases per 100,000 population (2005), with a nadir of 0.31 cases per 100,000 in 2004 [24][29].

The finding of a cluster spanning the Georgia-Tennessee border underscores the importance of looking across jurisdictional boundaries and the use of a common definition of clone. Cross-border spread has also been implicated in the increased incidence of cases seen in the Aachen region of Germany and neighboring Netherlands. [30] This edge effect was thought to explain the excess number of districts surrounding Aachen which had incidence rates higher than could be explained by complex space-time conditional intensity modeling of Germany alone. [31] In a study of 3,979 cases of meningococcal disease occurring from 2002–2008 in Niger, 15 clusters were identified. [32] Clusters ranged in size from 9 to 558 cases, and exhibited a geographic predominance in southeastern Niger. Movement across the border between Niger and Nigeria in this region is common and represents another example of this edge effect.

A previous spatiotemporal analysis of invasive meningococcal disease in Germany used a molecular subtyping scheme similar to ours with the exception that it did not include porB genotyping. [8] This approach, applied to 1,616 cases, resulted in the identification of 383 unique clones. 4.2% of cases were involved in a cluster, and 76.9% of clusters involved only 2 patients. Although the proportion of cases assigned to a cluster was lower in that study than ours (10.1%), a similarly high number of detected clusters included 2 patients (60% in our study).

Clone 2, a serogroup B strain, predominated in Oregon where 97 cases (85.8% of all Clone 2 isolates) were identified during the study period. These isolates represented 36.3% of all invasive meningococcal cases in Oregon from 2000–2005. As reported previously, the Oregon clone also caused disease in seven other ABCs sites, but it was not implicated in any statistically significant clusters in other states, confirming previous reports of its unique epidemiologic behavior in Oregon [11].

The SaTScan methodology has several advantages for the study of clustering among cases of invasive meningococcal disease. The software is highly flexible and can be used for spatial, temporal or space-time scan statistics, in either a retrospective or prospective fashion. Scan parameters, including spatial and temporal windows can be tailored to assess for specific types of endemic or epidemic patterns of activity. We purposely set relatively large windows with the goal of detecting patterns of transmission and disease beyond what would be defined as a community outbreak. [4] Retrospective analyses aid the confirmation of previously detected clusters. Most importantly, this methodology can be used to detect associations between cases without a previously identified epidemiological link. Prospective scans allow for the early detection of related cases and may signal persistence of a clone in space or time.

SaTScan allows the user to define clone characteristics, so parameters can be tailored to complement the level of molecular typing available. In this study, we used a highly conservative definition of a clone, namely, identification by MLST and genotyping of 3 OMP antigens. More clusters would have been detected if porB genotyping results had been restricted to the level of class 2 or 3 genotypes. If we had collapsed porB into class 2 or 3, we would have identified 352 unique clones (a 19.6% decrease), including 80 non-singletons. Using this classification, 790 cases would be attributable to a non-singleton clone, accounting for 74.4% of geocoded isolates. With porB collapsed into class 2 or 3 and all other scan parameters unchanged, we would have detected 50 clusters, as opposed to the 32 detected with the stricter definition used in this study. Consistent with the 2-case cluster reported by Connecticut, in which epidemiologically related isolates were indistinguishable except for a single locus difference by MLST, a less-conservative definition of clone, while leading to the detection of more clusters, would also have likely falsely clustered unrelated cases. Regardless, the specific SaTScan settings and definition of a clone can be tailored to the specific setting in which this methodology is employed.

In summary, the SaTScan methodology is a flexible and practical tool for the surveillance and analysis of meningococcal disease. Geotemporal scan software, in conjunction with molecular subtyping, can be used to study patterns of infectious disease transmission, identify previously undetected clusters, and improve understanding of the epidemiology of circulating molecular clones.

This publication made use of the Neisseria Multi Locus Sequence Typing website (http://pubmlst.org/neisseria/) developed by Keith Jolley and Man-Suen Chan and located at the University of Oxford. [33] The development of this site has been funded by the Wellcome Trust and European Union. We thank Martin Kulldorff for his methodological insight and technical assistance with the SaTScan software. SaTScan™ is a trademark of Martin Kulldorff. The SaTScan™ software was developed under the joint auspices of (i) Martin Kulldorff, (ii) the National Cancer Institute, and (iii) Farzad Mostashari of the New York City Department of Health and Mental Hygiene.

We thank the following Emerging Infections Program (EIP)/Active Bacterial Core surveillance (ABCs) site investigators and staff: Mirasol Apostol, Pam Daily, Joelle Nadle and Gretchen Rothrock (California); Steve Burnite (Colorado); Matt Carter and Michelle Wilson (Connecticut); Kathryn Arnold, Wendy Baughman, Lauren Lorentzson, Paul Malpiedi and Stephanie Thomas (Georgia); Terresa Carter and Rosemary Hollick (Maryland), Brenda Jewell, Christine Lees, Craig Morin and Jean Rainbow (Minnesota); Nancy Bennett, Jillian Karr, Glenda Smith and Nancy Spina (New York); Mark Schmidt and Ann Thomas (Oregon); Brenda Barnes (Tennessee); and Karrie-Ann Toews, Chris Van Beneden, Emily Weston and Carolyn Wright (CDC ABCs program). We also thank participating microbiology laboratory personnel and hospital infection preventionists in ABCs site hospitals and laboratories for identifying the N. meningitidis cases and submitting the bacterial isolates, Melina Lenser for performing molecular subtyping and Elizabeth Mitgang for her assistance with data management.

ReferencesHalperinSA, BettingerJA, GreenwoodB, HarrisonLH, JelfsJ, et al (2012) The changing and dynamic epidemiology of meningococcal disease. Vaccine30 Suppl 2B263622178525CohnAC, MacNeilJR, HarrisonLH, HatcherC, TheodoreJ, et al (2010) Changes in Neisseria meningitidis disease epidemiology in the United States, 1998–2007: implications for prevention of meningococcal disease. Clin Infect Dis50: 18419120001736BrooksR, WoodsCW, BenjaminDKJr, RosensteinNE (2006) Increased case-fatality rate associated with outbreaks of Neisseria meningitidis infection, compared with sporadic meningococcal disease, in the United States, 1994–2002. Clin Infect Dis43: 495416758417Centers for Disease Control and Prevention (1997) Control and prevention of serogroup C meningococcal disease: evaluation and management of suspected outbreaks: recommendations of the Advisory Committee on Immunization Practices (ACIP). MMWR Recomm Rep46: 13219048847BrehonyC, JolleyKA, MaidenMC (2007) Multilocus sequence typing for global surveillance of meningococcal disease. FEMS Microbiol Rev31: 152617168997HarrisonLH (2006) Prospects for vaccine prevention of meningococcal infection. Clin Microbiol Rev19: 14216416418528HarrisonLH, JolleyKA, ShuttKA, MarshJW, O’LearyM, et al (2006) Antigenic shift and increased incidence of meningococcal disease. J Infect Dis193: 1266127416586364EliasJ, HarmsenD, ClausH, HellenbrandW, FroschM, et al (2006) Spatiotemporal analysis of invasive meningococcal disease, Germany. Emerg Infect Dis12: 1689169517283618ReinhardtM, EliasJ, AlbertJ, FroschM, HarmsenD, et al (2008) EpiScanGIS: an online geographic surveillance system for meningococcal disease. Int J Health Geogr7: 3318593474JungI, KulldorffM, RichardOJ (2010) A spatial scan statistic for multinomial data. Stat Med29: 1910191820680984HarrisonLH, ShuttKA, SchminkSE, MarshJW, HarcourtBH, et al (2010) Population structure and capsular switching of invasive Neisseria meningitidis isolates in the pre-meningococcal conjugate vaccine era–United States, 2000–2005. J Infect Dis201: 1208122420199241SchuchatA, HilgerT, ZellE, FarleyMM, ReingoldA, et al (2001) Active bacterial core surveillance of the emerging infections program network. Emerg Infect Dis7: 929911266299HarrisonLH, DwyerDM, MaplesCT, BillmannL (1999) Risk of meningococcal infection in college students. JAMA281: 1906191010349894HarrisonLH, PassMA, MendelsohnAB, EgriM, RosensteinNE, et al (2001) Invasive meningococcal disease in adolescents and young adults. JAMA286: 69469911495619National Center for Health Statistics (2006) Vintage 2005 Bridged-Race Postcensal Population Estimates. Atlanta, GA.Kulldorff M (2010) SaTScan Technical Documentation. Boston, MA.KulldorffM (1997) A spatial scan statistic. Commun Stat-Theor M26: 14811496KulldorffM, NagarwallaN (1995) Spatial disease clusters: detection and inference. Stat Med14: 7998107644860TakahashiK, KulldorffM, TangoT, YihK (2008) A flexibly shaped space-time scan statistic for disease outbreak detection and monitoring. Int J Health Geogr7: 1418402711KulldorffM, HeffernanR, HartmanJ, AssuncaoR, MostashariF (2005) A space-time permutation scan statistic for disease outbreak detection. PLoS medicine2: e5915719066Kulldorff M, Information Management Services Inc. (2011) SaTScan: Software for the spatial and space-time scan statistics. 9.1.1 ed. Boston, MA and Silver Spring, MD,.U.S. Census Bureau Geography Division (2005) U.S. Census Bureau Cartographic Boundary Files Web Site. Washington, DC: U.S. Census Bureau, Geography Division, Cartographic Products Management Branch.U.S. Census Bureau (2000) American FactFinder Census 2000 Summary File 1 (SF 1) 100-Percent Data. Washington, DC: U.S. Census Bureau.Centers for Disease Control and Prevention (2002) Active Bacterial Core Surveillance Report, Emerging Infections Program Network, Neisseria meningitidis, 2000.Centers for Disease Control and Prevention (2003) Active Bacterial Core Surveillance Report, Emerging Infections Program Network, Neisseria meningitidis, 2001.Centers for Disease Control and Prevention (2003) Active Bacterial Core Surveillance Report, Emerging Infections Program Network, Neisseria meningitidis, 2002.Centers for Disease Control and Prevention (2004) Active Bacterial Core Surveillance Report, Emerging Infections Program Network, Neisseria meningitidis, 2003.Centers for Disease Control and Prevention (2005) Active Bacterial Core Surveillance Report, Emerging Infections Program Network, Neisseria meningitidis, 2004.Centers for Disease Control and Prevention (2006) Active Bacterial Core Surveillance Report, Emerging Infections Program Network, Neisseria meningitidis, 2005.EliasJ, SchoulsLM, van de PolI, KeijzersWC, MartinDR, et al (2010) Vaccine preventability of meningococcal clone, Greater Aachen Region, Germany. Emerging infectious diseases16: 46547220202422MeyerS, EliasJ, HohleM (2012) A space-time conditional intensity model for invasive meningococcal disease occurrence. Biometrics68: 60761621981412MainassaraHB, MolinariN, DematteiC, Fabbro-PerayP (2010) The relative risk of spatial cluster occurrence and spatiotemporal evolution of meningococcal disease in Niger, 2002–2008. Geospat Health5: 9310121080324JolleyKA, ChanMS, MaidenMC (2004) mlstdbNet - distributed multi-locus sequence typing (MLST) databases. BMC Bioinformatics5: 8615230973