In some locations, LD incidence recently has increased dramatically; in Virginia, the number of confirmed cases nearly tripled from 2006 to 2007 (http://www.vdh.virginia.gov/epidemiology/surveillance/surveillancedata/index.htm) to ≈12.4 cases per 100,000 residents, well above the 1998–2006 average of 2.2 per 100,000 (1). A 1990 report of LD cases in Virginia noted that the disease was rare in the early 1980s but apparently increased in incidence and geographic distribution through the late 1980s, leading the authors to conclude that the disease was expanding southward (11). Before 2006, most studies of I. scapularis ticks in Virginia focused on the eastern and southeastern parts of the state and found that densities of I. scapularisticks declined, as did their rate of infection with B. burgdorferi, with distance from the coast (12,13). Several early surveys for I. scapularis ticks in Virginia’s neighboring states of North Carolina and Maryland also found them to be most abundant on the Coastal Plain but absent or less common in the Piedmont and Appalachian Mountains. During 1983–1987, Apperson et al. surveyed 1,629 hunter-killed deer from the Coastal Plain, Piedmont, and Appalachian Mountain regions of North Carolina and found I. scapularis ticks only on deer from the Coastal Plain (14). Amerasinghe et al. surveyed 1,281, and 922 hunter-killed deer in 1989 and 1991, respectively, at sites from the Coastal Plain to the Appalachian Mountains of Maryland and found I. scapularis ticks on 59%–70% of deer on the Coastal Plain, fewer on deer in the Piedmont Region, and on only 1%–5% of deer in the Appalachian Mountains (15,16).

Although I. scapularis ticks exist in the southeastern United States (6–9), they are most easily detected by drag sampling, a method used as a proxy for risk to tick exposure (5), in areas associated with highest LD incidence, i.e., the Northeast (New Jersey through Massachusetts) and upper Midwest (Wisconsin and Minnesota) (3,5,17). The difference in apparent abundance of I. scapularis ticks and risk for LD between the northern and southeastern United States has been the subject of much discussion and debate (18) and might be related, either through behavioral or physiologic mechanisms, to genetic differences between I. scapularis populations in these regions (7,19–22). Population genetic structure of I. scapularis ticks has shown that dynamic range shifts are likely to have occurred in recent evolutionary history (19–22) and that 2 distinct lineages within this species can be identified; a relatively genetically uniform “American clade” exists in the northern United States (although this lineage has also been detected in the South), and a genetically diverse “southern clade,” members of which have been found only in the South (20). Although other nomenclatures have been proposed for these 2 lineages (e.g., clades A and B for northern and southern lineages, respectively [19]), we follow the terminology established by Norris et al.: “American” describes the widely distributed yet less diverse clade and “southern” describes the geographically restricted yet more diverse mtDNA clade of I. scapularis ticks (20).

Range expansion of I. scapularis ticks over relatively short periods has been observed (23,24). Moreover, recent environmental modeling, based on extensive field collections of host-seeking I. scapularis ticks, suggests that this species suggests that the range of this species is expanding widely and its occurrence in a given area depends on the lack of abiotic drivers, vapor pressure deficit and elevation (5,25). In Virginia, studies found that I. scapularisticks were concentrated in in northern sites; very few ticks were reported in other parts of the state (5,17,25). In contrast, human LD cases at inland, higher-elevation locations have increased in recent years in Virginia (http://www.vdh.virginia.gov/epidemiology/surveillance/surveillancedata/index.htm). The incongruity between human case and vector abundance datasets might be explained by recent (i.e., since 2007) spatial and/or numerical expansion of I. scapularis populations. We hypothesized that density ofB. burgdorferi–infected ticks would be highest in counties associated with high incidence of human disease if epidemiologic data represent cases in tick-endemic areas. In contrast, low numbers of infected ticks in areas of high human disease might indicate either misdiagnosis or allochthonous exposure.

Methods

Data from Cases in Humans

We compiled LD cases reported to the Virginia Department of Health (http://www.vdh.virginia.gov/epidemiology/surveillance/surveillancedata/index.htm) directly by physicians or identified through follow-up of positive laboratory results by county public health department personnel during 2000–2011. All LD cases counted in Virginia State Reportable Disease Reports met clinical and laboratory criteria specified in the National Surveillance Case Definition (SCD) for LD. We assessed cases counted in Virginia during 2000–2007 using criteria in the 1996 SCD and cases counted during 2008–2011 using criteria in the 2008 SCD. The 1996 SCD enabled states to liberally interpret what constituted laboratory evidence of infection; an IgM-positive Western immunoblot (WB IgM) test result could be counted as laboratory evidence of infection even though a more specific 2-tier test that used an enzyme-linked immunoassay (EIA) and the WB IgM was recommended. The Virginia Department of Health used the less restrictive interpretation of laboratory evidence in its LD surveillance from 1996 through the end of 2007. However, given that single-tier positive results from either the EIA or the WB IgM are less specific than a positive 2-tier result from both the EIA and the WB IgM (26,27), laboratory evidence of infection in the 2008 SCD required, at a minimum, a positive 2-tier test result on blood collected during the acute phase of illness (i.e., within 30 days after illness onset). The more stringent laboratory criteria adopted in the 2008 SCD were designed to minimize the number of false cases counted by state surveillance programs.

We analyzed all data at the county level, which required us to reclassify cases reported in cities to the counties in which they are situated because cities and counties in Virginia are often separate administrative entities. We estimated LD incidence per county for each year during 2000–2011 by dividing the annual number of counted cases by the estimated population size in 2007 (28). To characterize annual change in incidence per county, we calculated the difference in cases between successive years and then averaged these values across years. We analyzed the spatial distribution of human LD cases at the state level by identifying the centroid, or geometric center, of county-level LD incidence for each year, starting in 2000 using ArcMAP 10.0 (ESRI, Redlands, CA, USA). We then used weighted linear regression to determine the effect of year on latitude and longitude of that year’s centroid position weighted by annual number of cases.

Study Sites and Field Collections

In May and June 2011, we sampled ticks at 4 closed-canopy deciduous forest sites along an east-west elevational gradient: Crawfords State Forest (CR) (30 m), a University of Richmond–owned tract in Goochland County (GR) (80 m), Appomattox-Buckingham State Forest (AB) (170–200 m), and Lesesne State Forest (LE) (380–450 m) (Figure 1). We collected ticks at all sites by drag sampling (29) whereby a 1-m² piece of corduroy was dragged along both sides of 5 haphazardly selected 100-m transects (1,000 m² total), stopping every 20 m to remove ticks (17,25). We visited each site 4 times during May–July 2011 with at least 10 days separating visits. All ticks were speciated by light microscopy using dichotomous keys (30), and density of I. scapularisticks was calculated as the average number of ticks collected per transect. Difference in density of I. scapualrisnymphs among sites and visits was determined by analysis of variance of square root–transformed count data. We compared infection prevalence in ticks among sites by Gtest and by creating log-likelihood estimates of 95% CIs with a binomial probability function (31).

Molecular and Phylogenetic Methods

To extract total DNA, individual ticks were dried and flash-frozen by using liquid nitrogen, crushed by using a sterilized pestle, and processed with Qiagen DNeasy Blood and Animal Tissue Kit (QIAGEN, Valencia, CA, USA) by using manufacturer’s protocols. We tested for B. burgdorferi DNA by PCR amplification of the outer surface protein C (ospC) gene and the intergenic spacer region of 16S–23S rRNA genes (32). Presence of amplified DNA was determined by gel electrophoresis, and samples that produced amplicons were purified with a QIAquick PCR Purification Kit (QIAGEN) and submitted for sequencing at the Nucleic Acids Research Facility at Virginia Commonwealth University (Richmond, VA, USA). We also performed PCR to amplify and subsequently sequence an ≈460-bp portion of the I. scapularis 16S rRNA gene using primers 16S +1 and 16S –1 (20). Bidirectional chromatograms from all sequence data were assembled and initially analyzed with Sequencher 4.10.1 (Gene Codes, Ann Arbor, MI, USA). B. burgdorferi sequences were blasted by using GenBank (http://blast.ncbi.nlm.nih.gov/Blast.cgi) to confirm species identification. Sequence data from I. scapularis 16S samples were aligned with reference sequences (33) by using ClustalW (http://www.clustal.org) implemented in MEGA 5.0 (http://www.megasoftware.net/), which was also used to select among models of evolution and to reconstruct phylogeny.

Results

During 1995–1998, the Virginia Department of Health counted 55–73 LD cases per year. The number increased to 122 cases in 1999, and cases continued to increase through the early 2000s. Although Virginia’s LD activity during 2000–2005 was focused primarily on northern Virginia and the Eastern Shore of Virginia (a peninsula extending south from Maryland on the eastern side of the Chesapeake Bay), small numbers of LD cases were recorded in counties across Virginia, including counties in the most southern and southwestern parts (Figure 2). During 2006–2007 the incidence of LD increased substantially in counties throughout the Appalachian Mountains (Figure 2). After the change in the SCD in 2008, many of the most southern and southwestern counties that had recorded LD cases before 2008 ceased to report cases, and the geographic progression of LD appeared as a compact front that progressed from county to county from northeast to southwest. LD cases were not observed again in any of the far southwestern counties until 2011, by which time LD was considered endemic to many of the counties immediately to their northwest (Figure 2).

We collected 2,549 ticks from the field: 2,192 Amblyomma americanum (1 larva, 1,917 nymphs, 274 adults), 306I. scapularis (304 nymphs, 2 adults), 50 Dermacentor variabilis (all adults), and 1 I. dentatus (nymph). Sampling site was a major determinant of I. scapularis density (F = 71.07, p<0.0001, degree of freedom [df] = 3), as was sampling date (F = 6.85, p=0.024, df = 1). Post hoc comparisons indicated that tick density at the highest elevation site (9.55 nymphs/200 m²) was significantly greater than at any other site and that tick density at GR (1.66 nymphs/200 m²) was significantly higher than at site AB (0.25 nymphs/200 m²) (Table). We detected B. burgdorferi DNA in 48 I. scapularis nymphs, 45 of which produced unambiguous sequence reads for at least 1 locus (ospC or intergenic spacer region. Infection prevalence varied significantly among sites (likelihood ratio test, G = 16.3, p<0.0001, df = 3); the prevalence of infection was significantly higher at site LE (0.2) than at sites CR (0.00) and GR (0.04). Because of low sample size, site AB did not yield a reliable estimate of infection prevalence (Figure 3).

Analysis of I. scapularis 16S sequences yielded 17 haplotypes (GenBank accession nos. KF146631–47) from 85 individual nymphs (14 haplotypes from 44 ticks at LE, 1 from 2 ticks at AB, 6 from 21 ticks at GR, and 4 from 18 ticks at CR). Maximum-likelihood phylogenetic reconstruction using Tamura 3-parameter model (34) indicated that all haplotypes detected fall within the American clade; none of the ticks we sampled were phylogenetically identified as southern clade I. scapularis (Figure 4). In addition to an overall increase in human LD cases (from 136 in 2000 to an average of >1,000 in 2010 and 2011), we observed a significant spatial shift of the geometric center of LD incidence in Virginia. The longitude value associated with the centroid of each year’s LD incidence depended significantly on year from 2000 to 2011 (F = 12.48, p = 0.005, r² = 0.56) (Figure 5). Latitude values did not change significantly over time (F = 0.14, p = 0.71, r² = 0.01). We also calculated the average LD incidence per county for 2000–2006 (before the dramatic spike in cases in Virginia) and for 2007–2011 to identify counties in which the largest increases in cases occurred (Table).

Discussion

Our results indicate that 1) human LD incidence in Virginia has increased since 2000 and that the spatial distribution of cases has changed significantly, 2) abundance of I. scapularis nymphs and prevalence of B. burgdorferi infection are consistent with recent changes in human disease data, and 3) I. scapularis populations detected in central and western Virginia are dominated by American-clade haplotypes. Taken together, these results suggest recent spatial and/or demographic expansion of I. scapularis ticks in Virginia, resulting in increased human exposure to B. burgdorferi; the most notable increases in ticks and disease risk are at higher elevations in the western part of Virginia. More generally, our results indicate a dynamic pattern of LD risk. The spatial trends we identified through acarologic sampling are consistent with observed changes in disease incidence and are of paramount public health importance; the observed changes LD epidemiology in Virginia most likely reflect a spatial increase in disease endemicity (Table). We propose that the increase in LD in Virginia is caused by either increasing abundance of I. scapularis ticks, increasing prevalence of B. burgdorferiinfection in the vector, or both. Our data suggest that this vector species may be more abundant than it was before 2007; during widespread collections during 2004–2007, I. scapularis ticks existed throughout most of Virginia, and no infected I. scapularis ticks were detected in central or western Virginia (17,25). Similar range expansion of I. scapularis ticks has been described in Wisconsin and Michigan (23,24).

The extent to which the spatial distribution of LD cases in Virginia will continue to change is unclear. Environmental variables previously identified as important drivers of I. scapularis abundance may not have uniform effects throughout the range of this species. For example, on the basis of extensive sampling in the eastern United States over several years, Diuk-Wasser et al. estimated an elevational threshold of 510 m for this species (25), and Rosen et al. detected more I. scapularis on deer at low elevation than high elevation sites in Tennessee (35). However, our sampling showed the highest density of host-seeking I. scapularis nymphs at elevations approaching this threshold, and we have subsequently collected host-seeking nymphs at >1,000 m in Nelson County in west-central Virginia (R.J. Brinkerhoff, unpub. data). In 2007, a growing focus of LD incidence was observed in southwestern Virginia in Pulaski, Floyd, and Montgomery Counties. These counties have continued to have that region’s highest incidence of LD through 2011 (http://www.vdh.virginia.gov/epidemiology/surveillance/surveillancedata/index.htm) and mostly occupy high mountain valleys with average elevations of 584–762 m. An elevational threshold that limits tick populations at northern latitudes, where high elevation sites experience extreme cold during winter months, would not be expected where equivalent elevations are associated with more moderate climatic conditions.

Our analysis of LD data from humans indicates that the largest increases in LD incidence since 2007 has occurred in higher-elevation counties in western Virginia; the correspondence between these data and acarologic sampling suggests that the cases reported in these locations most likely are locally acquired and indicate recent spatial and/or numerical expansion of human disease. Our results are notably inconsistent with the findings of surveys of I. scapularis ticks on hunter-killed deer in North Carolina and Maryland during 1987–1992, which indicated that I. scapularis ticks were most abundant on the Coastal Plain and absent or uncommon in the Appalachian Mountains (14–16). When human LD surveillance began in Virginia in 1989, the highest incidence was on the state’s Eastern Shore (http://www.vdh.virginia.gov/epidemiology/surveillance/surveillancedata/index.htm). This finding was consistent with early surveys of ticks indicating that I. scapularis was the most common species in the Coastal Plain and much less common at higher elevations to the west (14). A logical conclusion at that time was that LD would continue to spread southward along the state’s Coastal Plain. However, during 2000–2011, LD became more prevalent in Virginia’s upper Piedmont and Appalachian Mountain zones than in the lower Piedmont and Coastal Plain. The results of older surveys of ticks and recent environmental models are not consistent with the current geographic incidence of LD or our field data. This discrepancy suggests a southwestward spatial expansion of northern tick populations into the upper Piedmont and mountain regions of Virginia or demographic expansion of persons into areas of previously low tick density in western localities. We do not have acarologic data from each county in which LD incidence has increased, nor do we have long-term systematic sampling data, and thus we cannot directly attribute local changes in LD to changes in tick densities.

Analysis of single-nucleotide polymorphisms in the I. scapularis genome reinforces the hypothesis that these ticks recolonized northern North America after the most recent glaciation event and that northern populations are genetically less diverse than southern populations (21). Moreover, analyses of single-nucleotide polymorphism data are consistent with south-to-north postglaciation gene flow, whereby northern American-clade populations are a subset of the genetic variation found in southern-clade populations (21) resulting from founder effects when ticks recolonized northern latitudes (22). Tick populations within both LD-endemic foci show signs of genetic isolation from one another and from southern populations (22), and evidence exists for similar lack of gene flow among populations within regions (19). Identification of American-clade I. scapularis ticks in the southeastern United States (19,33) might reflect remnant American-clade lineages in the South or might indicate southward dispersal of American-clade ticks. Qiu et al. noted that coastal sites in southern states were associated with strictly American-clade populations, whereas a mix of American- and southern-clade ticks was detected at inland sites (19). With respect to our study, we point to the recent lack of detection of I. scapularis ticks at high-elevation sites in western or central Virginia (17,25) and the presence of exclusively American-clade I. scapularis ticks in the current study as possible evidence consistent with the population expansion of American-clade ticks from northern population foci. However, we cannot exclude the possibility that the distribution of endemic American-clade ticks simply has expanded in Virginia.

Although American- and southern-clade I. scapularis ticks are now considered 1 species, apparent differences exist in host-seeking behavior, biting behavior, and duration of attachment to different host types (9,36,37). Genetic differences between the major I. scapularis lineages have been well documented (7,19–22), and if American-clade ticks are more likely to feed on humans, the emergence of LD in Virginia would be consistent with increased relative abundance of this variant. In the South, immature I. scapularis ticks feed predominantly on low-competence or noncompetent lizard species and are relatively uncommon on rodents (8,36–38). Southern-clade nymphs may have questing behavior that makes them unlikely to be collected on cloth drags or to bite humans (9); thus, nymphal ticks are difficult to collect, even in places where adult ticks are common. LD risk should be very low in areas where I. scapularis nymphs are unlikely to bite humans, and immature ticks are more likely to feed on reptiles than on competent vertebrate reservoirs. However, data from a single mitochondrial gene, albeit one that has been widely characterized for this species, do not necessarily reflect patterns of differentiation found in nuclear markers (21) and probably are not useful for delineating among behavioral phenotypes. Moreover, we sampled in daytime hours during the presumed peak period of nymphal activity (late spring, early summer) and thus would not have detected ticks exhibiting different host-seeking behaviors. It is possible that multilocus genomic analysis or year-round sampling would yield different insights from those reached in this study.

The latitudinal gradient in LD risk in the eastern United States is not easily explained and probably is driven by demographic and environmental factors (5,26,39). However, our data suggest that the boundary between regions to which I. scapularis ticks are and are not endemic is moving and that B. burgdorferi–infected ticks might be expanding in or into areas from which they historically have been absent. As a result, clinicians and epidemiologists need to be vigilant in the face of changing spatial distributions of risk, especially in transition zones where patterns of disease are rapidly changing (40).

Dr Brinkerhoff is an assistant professor at the University of Richmond and holds an honorary senior lectureship in the School of Life Sciences, University of KwaZulu-Natal, Pietermaritzburg, South Africa. His research focuses on the ecology, evolution, and epidemiology of bacterial pathogens transmitted by arthropod vectors.

Acknowledgments

We thank M. Massaro, R. Kelly, and L. Spicer for assistance with sample collection and processing of laboratory samples. We also thank 2 anonymous reviewers for their thoughtful comments on an earlier draft of this manuscript.

Funding for this project was provided by a Thomas F. and Kate Miller Jeffress Memorial Trust Award (J-1036) to R.J.B. and a summer research fellowship to W.F.G. provided through the University of Richmond School of Arts and Sciences.

References

1. Bacon RM, Kugeler KJ, Mead PS. Surveillance for LD—United States, 1992–2006. MMWR Surveill Summ. 2008;57:1–9 .PubMed
2. Barbour AG, Fish D. The biological and social phenomenon of LD. Science. 1993;260:1610–6. DOIPubMed
3. Pepin KM, Eisen RJ, Mead PS, Piesman J, Fish D, Hoen AG, Geographic variation in the relationship between human LD incidence and density of infected host–seeking Ixodes scapularis nymphs in the eastern United States. Am J Trop Med Hyg. 2012;86:1062–71. DOIPubMed
4. Centers for Disease Control and Prevention. LD incidence rates by state, 2002–2011 [cited 2013 May 7].http://www.cdc.gov/lyme/stats/chartstables/incidencebystate.html
5. Diuk-Wasser MA, Hoen AG, Cislo P, Brinkerhoff R, Hamer SA, Rowland M, Human risk of infection with Borrelia burgdorferi, the LD agent, in eastern United States. Am J Trop Med Hyg. 2012;86:320–7. DOIPubMed
6. Goddard J. Ecological studies of adult Ixodes scapularis in central Mississippi: questing activity in relation to time of year, vegetation type, and meteorologic conditions. J Med Entomol. 1992;29:501–6 .PubMed
7. Rich SM, Caporale DA, Telford SR, Kocher TD, Hartl DL, Spielman A. Disctribution of the Ixodes ricinus–like ticks of eastern North America. Proc Natl Acad Sci U S A. 1995;92:6284–8. DOIPubMed
8. Durden LA, Oliver JH Jr, Banks CW, Vogel GN. Parasitism of lizards by immature stages of the blacklegged tick, Ixodes scapularis (Acari, Ixodidae).Exp Appl Acarol. 2002;26:257–66. DOIPubMed
9. Goddard J, Piesman J. New records of immature Ixodes scapularis from Mississippi. J Vector Ecol. 2006;31:421–2. DOIPubMed
10. Centers for Disease Control and Prevention. Incidence by state, 2005–2010 [cited 2012 Jul 1].http://www.cdc.gov/lyme/stats/chartstables/incidencebystate.html
11. Heimberger T, Jenkins S, Russell H, Duma R. Epidemiology of LD in Virginia. Am J Med Sci. 1990;300:283–7. DOIPubMed
12. Casteel M, Sonenshine D. Abundance of adult Ixodes scapularis and infection with Borrelia burgdorferi in eastern Virginia. Va J Sci.1996;47:293–300.
13. Sonenshine DE, Ratzlaff RE, Troyer J, Demmerle S, Demmerle ER, Austin WE, Borrelia burgdorferi in eastern Virginia: comparison between a coastal and inland locality. Am J Trop Med Hyg. 1995;53:123–33 .PubMed
14. Apperson CS, Levine JF, Nicholson WL. Geographic occurrence of Ixodes scapularis and Amblyomma americanum (Acari: Ixodidae) infesting white-tailed deer in North Carolina. J Wildl Dis. 1990;26:550–3. DOIPubMed
15. Amerasinghe FP, Breisch NL, Azad AF, Gimpel WF, Greko M, Neidhardt MK, Distribution, density and LD spirochete infection in Ixodes dammini(Acari: Ixodidae) on white-tailed deer in Maryland. J Med Entomol. 1992;29:54–61 .PubMed
16. Amerasinghe FP, Breisch NL, Neidhardt K, Pagac B, Scott TW. Increasing density and Borrelia burgdorferi infection of deer-infesting Ixodes dammini (Acari:Ixodidae) in Maryland. J Med Entomol. 1993;30:858–64 .PubMed
17. Diuk-Wasser MA, Hoen AG, Cortinas R, Yaremych-Hamer S, Tsao J, Kitron U, Spatiotemporal patterns of host-seeking Ixodes scapularis nymphs (Acari: Ixodidae) in the United States. J Med Entomol. 2006;43:166–76. DOIPubMed
18. Stromdahl EY, Hickling GJ. Beyond Lyme: aetiology of tick-borne human diseases with emphasis on the south-eastern United States. Zoonoses Public Health. 2012;59(Suppl. 2):48–64. DOIPubMed
19. Qiu WG, Dykhuizen DE, Acosta MS, Luf BJ. Geographic uniformity of the LD spirochete (Borrelia burgdorferi) and its shared history with tick vector (Ixodes scapularis) in the northeastern United States. Genetics. 2002;160:833–49 .PubMed
20. Norris DE, Klompen JSH, Keirans JE, Black WIC. Population genetics of Ixodes scapularis (Acari: Ixodidae) based on mitochondrial 16S and 12S genes. J Med Entomol. 1996;33:78–89 .PubMed
21. Van Zee J, Black WC, Levin M, Goddard J, Smith J, Piesman J. High SNP density in the black-legged tick, Ixodes scapularis, the principal vector of LD spirochetes. Ticks Tick Borne Dis. 2013;4:63–71.
22. Humphrey PT, Caporale DA, Brisson D. Uncoordinated phylogeography of Borrelia burgdorferi and its tick vector Ixodes scapularis. Evolution.2010;64:2653–63. DOIPubMed
23. Hamer SA, Tsao JI, Walker ED, Hickling GJ. Invasion of the LD vector Ixodes scapularis: implications for Borrelia burgdorferi endemicity.EcoHealth. 2010;7:47–63. DOIPubMed
24. Lee X, Hardy K, Johnson DH, Paskewitz SM. Hunter-killed deer surveillance to assess changes in the prevalence and distribution of Ixodes scapularis (Acari:Ixodiae) in Wisconsin. J Med Entomol. 2013;50:632–9. DOIPubMed
25. Diuk-Wasser MA, Vourc’h G, Cislo P, Hoen AG, Melton F, Hamer SA, Field and climate-based model for predicting the density of host-seeking nymphal Ixodes scapularis, an important vector of tick-borne diseases agents in the eastern United States. Glob Ecol Biogeogr. 2010;19:504–14.
26. Centers for Disease Control and Prevention. Recommendations for test performance and interpretation from the Second National Conference on Serologic Diagnosis of LD. MMWR Morb Mortal Wkly Rep. 1995;44:590–1 .PubMed
27. Engstrom SM, Shoop E, Johnson RC. Immunoblot interpretation criteria for serodiagnosis of early LD. J Clin Microbiol. 1995;33:419–27.PubMed
28. ESRI. ESRI demographic updates: 2012/2017, and ESRI white paper. ESRI. 2012; http://www.esri.com/library/whitepapers/pdfs/demographic-update-methodology-2012.pdf
29. Falco RC, Fish D. A comparison of methods for sampling the deer tick, Ixodes dammini, in a LD-endemic area. Exp Appl Acarol. 1992;14:165–73.DOIPubMed
30. Sonenshine DE. Ticks of Virginia (Acari, Metastigmata). Insects of Virginia series, no. 13. Blacksburg (VA): Virginia Polytechnic Institute and State University Press. 1979. 35 pp.
31. Brinkerhoff RJ, Collinge SK, Bai Y, Ray C. Are carnivores universally good sentinels of plague? Vector Borne Zoonotic Dis. 2009;9:491–7.DOIPubMed
32. Bunikis J, Garpmo U, Tsao J, Berglund J, Fish D, Barbour AG. Sequence typing reveals extensive strain diversity of the Lyme borreliosis agentsBorrelia burgdorferi in North America and Borrelia afzelii in Europe. Microbiology. 2004;150:1741–55. DOIPubMed
33. Trout RT, Steelman CD, Szalanski AL. Population genetics and phylogeography of Ixodes scapularis from canines and deer in Arkansas. Southwest Entomologist. 2009;34:273–87. DOI
34. Tamura K. Estimation of the number of nucleotide substitutions when there are strong transition-transversion and G + C content biases. Mol Biol Evol. 1992;9:678–87 .PubMed
35. Rosen ME, Hamer SA, Gerhardt RR, Jones CJ, Muller LI, Scott MC, Borrelia burgdorferi not detected in widespread Ixodes scapularis (Acari: Ixodidae) collected from white-tailed deer in Tennessee. J Med Entomol. 2012;49:1473–80. DOIPubMed
36. Oliver JH Jr, Cummins GA, Joiner MS. Ixodes scapularis (Acari: Ixodidae) parasitizing lizards from the southeastern USA. J Parasitol.1993;79:684–9. DOIPubMed
37. Oliver JH Jr, Owsley MR, Hutcheson HJ, James AM, Chen C, Irby WS, Conspecificity of the ticks Ixodes scapularis and I. dammini (Acari: Ixodidae).J Med Entomol. 1993;30:54–63 .PubMed
38. Spielman A, Levine JF, Wilson ML. Vectorial capacity of North American Ixodes ticks. Yale J Biol Med. 1984;57:507–13 .PubMed
39. Estrada-Peña A. Increasing habitat suitability in the United States for the tick that transmits LD: a remote sensing approach. Environ Health Perspect. 2002;110:635–40. DOIPubMed
40. Lantos PM, Brinkerhoff RJ, Wormser GP, Clemen R. Empiric antibiotic treatment of erythema migrans–like skin lesions as a function of geography: a clinical and cost effectiveness modeling study. Vector Borne Zoonotic Dis. 2013;13:877–83. DOIPubMed

Figures

• Figure 1. Locations of 4 field sites at which ticks were sampled, Virginia, May-July 2011. Circles indicate sampling areas. LE, Lesesne State Forest; AB, Appomattox-Buckingham State Forest; GR, University of Richmond–owned field site; CR, Crawfords...
• Figure 2. Progressive geographic spread of human Lyme disease across Virginia, 2001–2011. Data were reported by the Virginia Department of Health http://www.vdh.virginia.gov/epidemiology/surveillance/surveillancedata/index.htm. Cases per 100,000 population were calculated by county or...
• Figure 3. Variation in estimated prevalence of Borrelia burgdorferi infection in Ixodes scapularis nymphs at 4 field sites in Virginia. Sites are arranged west to east from left to right. LE,...
• Figure 4. Maximum-likelihood phylogenetic reconstruction of Ixodes scapularis lineages based on 16S rRNA gene sequences using Tamura 3-parameter model (35). All samples beginning with IS were collected during this study; reference......
• Figure 5. Centroids of annual incidence, by county, Virginia, 2000–2011. The size of each circle represents the annual number of cases reported by the Virginia Department of Health and is proportional...

Table

• Table. Average density of host-seeking Ixodes scapularis nymphs, average prevalence of Borrelia burgdorferi–infected nymphs, and average annual LD incidence in 5 counties, Virginia

Suggested citation for this article: Brinkerhoff RJ, Gilliam WF, Gaines D. Lyme disease, Virginia, USA, 2000–2011. Emerg Infect Dis [Internet]. 2014 Oct [date cited]. http://dx.doi.org/10.3201/eid2010.130782

DOI: 10.3201/eid2010.130782

The reference has no authors. Please proof carefully. (in reference 6 "Goddard J. Ecological studies of, 1992").

Medline indexes "Va J Sci" but cannot find a listing for reference 12 "Casteel, Sonenshine, 1996". Please check the reference for accuracy.

Medline indexes "Glob Ecol Biogeogr" but cannot find a listing for reference 25 "Diuk-Wasser, Vourc’h, Cislo, Hoen, Melton, Hamer, et al., 2010". Please check the reference for accuracy.

Medline cannot find the journal "Southwest Entomologist" (in reference 33 "Trout, Steelman, Szalanski, 2009"). Please check the journal name.

Table of Contents – Volume 20, Number 10—October 2014

Link Here