Abstract
Background:
In the United States, spillover of West Nile virus (WNV) into wild mammal populations has been reported since the introduction of the virus into the New World in 1999. Eastern gray squirrels (Sciurus carolinensis) exhibit a high seroprevalence for WNV in urban settings where high virus circulation and human spillover have been reported. In Atlanta, Georgia, human cases of WNV are uncommon despite high infection rates in birds and mosquitoes. In this study, we evaluated WNV exposure of eastern gray squirrels in a WNV hot spot in Atlanta.
Materials and Methods:
Gray squirrels were live-trapped in Grant Park, Atlanta, during July–October, 2012, and a census was conducted to estimate squirrel density in the study site. Sera from trapped animals were tested for circulating virus-by-virus isolation in cell culture and for WNV-specific antibodies by enzyme-linked immunosorbent assay and plaque reduction neutralization test. Mosquitoes were collected at the same location and tested for virus isolation.
Results:
Among the 69 collected squirrels, 25 (36.2%) tested positive for WNV antibodies, although none were viremic. Seroprevalence was lower in juveniles (18.8%) than in adults (37.5%), but this difference was not statistically significant. Gender and squirrel density had no effect on seroprevalence. Seasonality of squirrel seroprevalence and of mosquito infection were significantly associated, both peaking in August. No difference in squirrel exposure was detected across the collection sites.
Conclusions:
We report a high degree of WNV exposure in squirrels in Grant Park that was correlated with seasonality of mosquito infection. The detection of antibodies in juveniles suggests that circulation of WNV in the surveyed population is ongoing. Eastern gray squirrels may be suitable indicators for virus amplification and for risk of human spillover on a local scale in urban settings.
Introduction
W
WNV exposure has been detected in many other mammals, such as bats, raccoons, opossum, rabbit, mice, skunks, and tree squirrels (Dietrich et al. 2005, Root et al. 2005, Gómez et al. 2008a, Blitvich et al. 2009). Most mammalian species are considered dead-end hosts because the viremia levels following exposure are not high enough to infect the vector (Bunning et al. 2002, Ratterree et al. 2004, Read et al. 2005). Over the last decade, several studies performed in both laboratory conditions and natural settings re-examined whether specific mammals are noncompetent hosts. Studies of viremia in cottontail rabbits, eastern chipmunks, golden hamsters, and tree squirrels demonstrated that these species presented viremia above the threshold required to infect Culex mosquitoes (Tesh et al. 2005, Tiawsirisup et al. 2005, Root et al. 2006, Platt et al. 2007, Gómez et al. 2008b).
Field studies have documented that tree squirrels and chipmunks are commonly exposed to WNV with mean population seroprevalence comparable to levels reported in bird populations (Dietrich et al. 2005, Root et al. 2005, Gómez et al. 2008a, Blitvich et al. 2009). In squirrels, laboratory experiments detected virus in blood up to 5 days postinoculation, with a viremia peak of 104.98 plaque-forming units (pfu)/mL. WNV has also been isolated from fecal and urine samples obtained from infected squirrels up to 20 days after initial infection (Root et al. 2006, Gómez et al. 2008b).
Urban green spaces are favorable habitats for sciurids (e.g., tree squirrels and chipmunks), which typically form big colonies in these settings, sharing this urban environment with mosquito vectors, birds, and humans. The home range of squirrels is limited compared to that of birds, and, thus, they may indicate the potential for virus spillover to mammals on a local scale (Don 1983, Koprowski 1994, Whitaker and Hamilton 1998). Eastern gray squirrels may be a useful indicator of local virus transmission because they breed twice annually, and seroprevalence in juveniles can be used to detect virus transmission during the WNV transmission season (Whitaker and Hamilton 1998). Furthermore, the levels of viremia and seroprevalence in sciurids suggest that they also may indicate risk of WNV spillover to humans and may even play a role in WNV transmission.
In Georgia, although rates of seroprevalence in bird reservoirs and of mosquito infection are similar to those observed in areas with high levels of human spillover (Gibbs et al. 2006, Vazquez-Prokopec et al. 2010), the number of reported human cases is relatively low (Centers for Disease Control and Prevention 2014). Recently, WNV hot spots have been identified in the Atlanta metropolitan area (Vazquez-Prokopec et al. 2010) in and around wooded areas (the city of Atlanta is known as the “city in the forest” because of its high tree coverage). This type of landscape favors high abundance of tree squirrels (Sciurus carolinensis), which, given their arboreal behavior (squirrels usually build their nests at the fork of tree branches or inside tree holes), share their habitat with birds. There is also a positive correlation between tree coverage and a high abundance of Culex (Cx.) quinquefasciatus in Atlanta (Vazquez-Prokopec et al. 2010), which may increase the amount of contact between mosquitoes and tree squirrels and therefore, the risk of WNV spillover.
Here, we report the results of a study targeting squirrels to help explain the relatively low number of reported human cases in Georgia (Centers for Disease Control and Prevention 2014, Ruiz et al. 2004, 2007). We focused on the possible role of squirrels in WNV transmission in Grant Park, an urban park in Atlanta. In this area, there has been reported evidence of spillover of WNV to humans (Vazquez-Prokopec et al. 2010), and WNV has been detected in birds and mosquitoes in and around Grant Park annually since 2001 (D.G. Mead, personal communication). We also examined whether WNV seroprevalence in tree squirrels could be used as an indicator for local-scale WNV transmission risk to humans in an urban setting.
Materials and Methods
Study area
Between June and October of 2012, we captured squirrels and collected mosquitoes in Grant Park, Atlanta, GA. This park has high tree density and is heavily frequented by humans. The study area has described earlier (Levine et al. 2013), and WNV infection in birds has been documented (Levine et al. 2013).
Squirrel trapping and processing
We collected specimens in seven sites with a mean distance of 150 meters between sites. Trapping was performed biweekly using five Tomahawk traps (48.33×15.23×15.2 cm; Tomahawk Live Trap, Tomahawk, WI) in an area of ∼50 meters2 per site for a total of 280 trap/days. Traps were baited with peanut butter and a mixture of seeds. At each site, the five traps were set approximately 5 meters apart, and opened during mornings (6:00 to 12:00) or afternoons (16:00 to 20:00), periods that correspond to the daily peaks of squirrel activity. Because the park is heavily used, we checked the traps every 45 min to ensure that they were not disturbed.
All captured squirrels were anesthetized according to published guidelines (Parker et al. 2008) by inhalation of isoflurane (Iso Flo1, Abbott Laboratories, North Chicago, IL). We aged each animal as adult or juvenile based on size, fur color, and breeding attributes (Koprowski 1994). Blood samples were collected from the femoral vein, transported on ice to the laboratory, and centrifuged at 10,000 rpm for 30 min. After centrifugation, the serum was removed and stored at −80°C for further processing. All procedures were approved by Emory University IACUC committee (DAR-2001178-051515B, approved 5/15/2012).
Squirrel density estimation method
To evaluate squirrel density, two squirrel censuses were conducted at the beginning (early July, 2012) and the end (late September, 2012) of the trapping season. We applied area counts to estimate squirrel density in our study (Hein 1997, Steele and Koprowski 2003). We divided the study area into 30 equal-sized quadrants of 100 meters2 and established a vantage point within each quadrant, from which each operator started the surveillance facing north and covered the quadrant turning clockwise. Each quadrant was observed for 15 min, and the operator recorded the distance and direction of each sighted squirrel from the vantage point.
The sampling order of each area was predetermined using random number generation (Flyger 1959, Parker and Nilon 2008). Counts were conducted in each area from sunrise to 2 h after sunrise and 2 h prior to sunset until sunset, which corresponds to peak squirrel activity in urban areas (Manski et al. 1981, Gustafson and VanDruff 1990). To reduce the probability of count bias due to disturbance factors, we visited each quadrant both in mornings and evenings on different days. Density was calculated as described by Flyger (1959).
Mosquito collection and processing
Adult mosquitoes were collected using CDC gravid traps (Reiter et al. 1986) baited with a mixture of dog food and hay infusion. One trap per site was set night (from 18:00 to 8:00) before the squirrel-trapping session. Collected female mosquitoes were identified to genus (Slaff and Apperson 1989), sorted into pools of up to 25 (according to species, trap site, and collection date), and stored at −80°C. Mosquito pools were homogenized using a Qiagen Mixer Mill 300 (Valencia, CA) and clarified by centrifugation (10 min at 9000 rpm) for virus isolation performed at the Southeastern Cooperative Wildlife Disease Study at the University of Georgia (Athens, GA).
Enzyme-linked immunosorbent assay
Sera were screened at a dilution of 1:10 for antibodies to flaviviruses by epitope-blocking enzyme-linked immunosorbent assay (ELISA), as described previously (Blitvich et al. 2003, 2009). Briefly, ELISAs were performed using the WNV-specific monoclonal antibody (mAb) 3.1112G (Chemicon International, Temecula, CA) or the flavivirus group-reactive mAb 6B6C-1 (InBios International, Seattle, WA). Antigen was prepared from WNV-infected Aedes albopictus (C6/36) cell cultures. The ability of the test sera to block the binding of the mAbs to WNV antigen was compared with the blocking ability of control serum without antibody to flaviviruses. Data were expressed as relative percentages, and inhibition values ≥30% were considered as indicating the presence of viral antibodies.
Plaque reduction neutralization test
All sera positive for flavivirus antibodies by blocking ELISA were further tested by plaque reduction neutralization test (PRNT) to identify the infecting virus. PRNTs were conducted according to standard methods in the Biosafety Level 3 facilities at Iowa State University (Beaty et al. 1995) using WNV (strain NY99–35261–11) and St. Louis encephalitis virus (SLEV; strain TBH-28). SLEV was included because it has been identified in Georgia and is known to react with antibodies to WNV (Calisher et al. 1989, Gubler et al. 2007). Sera were initially tested at a dilution of 1:10, and those that tested positive were further diluted and tested to determine their end-point titers. Titers were expressed as the reciprocal of highest serum dilutions yielding ≥90% reduction in the number of plaques (PRNT90). For etiologic diagnosis, the PRNT90 antibody titer to the respective virus needed to be at least four-fold greater than that shown by the other tested flavivirus.
Virus isolation in cell culture
All squirrel sera and mosquito homogenates were tested for cytopathic virus by performing virus isolation in Vero Middle America Research Unit (MARU, Vero M) cells. Cell cultures were examined regularly for evidence of cytopathic effect (CPE). CPE-positive cultures were tested for WNV using VecTest® strips (Medical Analysis Systems, Inc., Camarillo, CA). Reverse transcription polymerase chain reaction (RT-PCR) was performed to confirm the presence of WNV RNA in VecTest-positive samples using WNV-specific primers as described by Levine et al. (2013). The minimum infection rate (MIR) was calculated using maximum likelihood (ML) (Bustamante and Lord 2010) and is expressed as the number of positive mosquito pool per 1000 tested mosquitoes.
Statistical analysis
The Fisher least significant difference (LSD) test (Conover 1998) was used to evaluate the differences between the proportion of seropositive squirrels by study site, age group (adult/juvenile), and gender. The Fisher LSD test was also applied to test for significant differences between mosquito MIR by site. Differences in abundance of collected mosquitoes between months were compared using a Mann–Whitney–Wilcoxon test.
The effect of squirrel gender, age, study site, and monthly MIR on animal seropositivity was tested using a logistic regression based on a generalized linear model (GLM). Starting from a full model formula:
We applied a multimodel selection approach to find the best model(s) using Akaike information criteria (AIC) (Burnham and Anderson 2002). We calculated the ΔAIC as the difference between each model and the one with the lowest AIC value. All models with a ΔAIC <2 were included in the set of best models (Burnham and Anderson 2002). A cumulative link mixed models (CLMM) procedure was performed to investigate the relationship of antibody titer (ordered categorical data) in seropositive squirrels with specimen characteristics and capture month. In the CLMM model, sampling sites were set as a random effect to adjust for environmental components and pseudo-replication.
Results
Squirrel census and collection
The density of squirrels for the study area was determined to be 18.9/ha, with a range of 14.8–22/ ha (Table 1). During our study, we collected 69 individuals with three (4.2%) recaptures. All three recaptured squirrels were retrapped at their original trapping site, suggesting a limited movement range. Of the 69 individuals, 44 were males (61.4%) and 16 were juveniles (22.8%). Thirty-one squirrels (44.8%) showed signs of breeding attributes (e.g., testicular dimension, vagina shape, and pregnancy). The lowest number of squirrels was captured in July (Fig. 1), when the collections were affected by consecutive rain events within a short time period (data not shown), and we were not able to set traps until the 18th of the month.

Collections of female mosquitoes and squirrels, June–October, 2012. (
MIR, mnimum infection rate; CI, confidence interval.
We collected a total of 72 sera from the 69 squirrels. Twenty-five (36.2%; 95% confidence interval [CI] 24.9–48.6%) sera tested positive for WNV antibodies (Table 1). Among the seropositive individuals, three (12%) were juveniles. Seroprevalence of the adult group (37.5%; 95% CI 24.9–51.5%) was higher than in juveniles (18.8%; 95% CI 4.1–45.6%), but this difference was not significant (odds ratio [OR]=2.7; 95% CI 0.61–15.68; p>0.05). Seroprevalence in males (32.6%; 95% CI 19.1–48.5%) and females (34.5%; 95% CI 17.9–54.3%) was not significantly different (OR=1.09; 95% CI 0.35–3.29, p>0.05). No significant differences (OR=2.33; 95% CI 0.77–7.25; p>0.05) were detected between the seroprevalence of squirrels showing breeding signs (44.8%; 95% CI 26.4–64.3%) and nonbreeding individuals (25.6%; 95% CI 13.5–41.2%). A significant difference (Fisher test, p<0.05) was found between the number of seropositive squirrels in August, when prevalence peaked (70.6%; 95% CI 44–89.6%), and July when none of the individuals captured were seropositive (0%; 95% CI 0–33.6%).
WNV PRNT90 titers of positive squirrels ranged from 40 to 1280 (Fig. 2). Only squirrels collected in the summer had antibody titers higher than 320 (Fig. 2), and the highest value of antibody titer (1280) was recorded in an adult squirrel in June. Juvenile squirrels had antibody titers ranging from 160 to 640, and the highest value for this age group was recorded in August (Fig. 2). Viremia was not detected in any squirrels.

Seroprevalence and antibody titers of collected squirrels and minimum infection rate (MIR) of mosquitoes. (
Mosquito collection
We collected a total of 4276 female mosquitoes, of which the majority (n=4131; 96.5%) were Culex (Cx.) spp., which were grouped into 287 pools. The remaining 146 specimens belonged to other genera and were grouped into 45 pools (unpublished data). In June, we collected the highest number of female mosquitoes per gravid trap (Fig. 1), and the number of collected mosquitoes dropped significantly (Mann–Whitney–Wilcoxon test, p<0.05) in July, and increased slightly in August and September (Fig. 1). The lowest number of collected mosquitoes per gravid trap was recorded in October (Fig. 1).
WNV was only detected in pools containing Culex spp. The WNV MIR for the entire study area was 5.8 (95% CI 3.7; 8.6) and ranged from 2.7 (95% CI 0.8; 6.2) to 9.8 (95% CI 0.2; 53.4) across sites. Mosquito infection peaked in August, and WNV was not detected in pools collected during September or October (Fig. 2). The lower WNV MIR recorded during July was not significantly different from June and August (Fisher LSD test, p=0.23). No significant difference in mosquito MIRs was observed between the sites (Fisher LSD test, p=0.34).
Model results
Two models were selected as the best candidates (Table 2). The variables included in these two best models were the monthly mosquito MIR (Σω=0.70) and squirrel age (Σω=0.38). The monthly MIR was significant and positively associated with squirrel seropositivity. Young squirrels had a lower probability of showing seropositivity than adults, but this difference was not significant. Antibody titers in September and October were significantly lower than in August according to the CLMM (Table 3). Gender and age were not correlated with antibody titer.
The Σωi for each variable is included.
p<0.05.
August was chosen as the reference month because it had the highest seroprevalence.
p<0.05.
p<0.01.
CI, confidence interval.
Discussion
Several studies have shown that spillover of WNV from the mosquito–bird cycle can be detected in many mammal species, but its extent and impact have been studied in depth only for humans. Previous studies reporting WNV seroprevalence in wild mammal communities did not address seasonality of virus spillover. In our study, we collected data on a fine temporal and spatial resolution, which allowed us to quantify monthly spatial and temporal spillover dynamics for a squirrel population living in a WNV hot spot in Atlanta.
Our findings point to WNV spillover into squirrels living in Grant Park (overall seroprevalence of 37.5%). This is similar to the seroprevalence levels recorded in other populations of S. carolinensis in the eastern and midwestern United States, including New York (Kramer and Bernard 2001, Marfin et al. 2001, Root et al. 2005), Illinois (Heinz-Taheny et al. 2004), Pennsylvania (Root et al. 2005), Louisiana (Dietrich et al. 2005), Maryland, and Washington DC (Gómez et al. 2008a).
The highest proportion of seropositive squirrels was reported in Grant Park in August, when the WNV MIR in Culex spp. mosquitoes and seroprevalence in birds captured in the study area also peaked in both 2010 and in 2011 (unpublished data). Thus, the high seroprevalence of WNV in squirrels in Grant Park corresponded to the period of maximum WNV transmission between vectors and avian reservoirs.
Spillover to squirrel populations may be facilitated by them sharing their habitat with avian reservoirs. Moreover, Cx. quinquefasciatus is one of the most abundant host-seeking mosquito species in canopy habitats (Savage et al. 2008), and a high prevalence of WNV infection was reported in mosquitoes collected in tree canopies (Anderson et al. 2004). This co-habitation contributes not only to the mosquito–bird WNV enzootic cycle, but, because Culex spp. females are opportunistic and also feed on several mammalian species (Apperson et al. 2002, Molaei et al. 2006, Hamer et al. 2008), it could explain the high seroprevalence in squirrels recorded here and in other studies.
The antibody titers of specimens collected in fall were significantly lower than in summer. This was consistent with absence of positive mosquito pools from this period. We also reported a low squirrel seroprevalence (no positive specimens) and mosquito MIR in July. In our study year, during July, 2012, we recorded the highest level of rainfall (8.2 cm) with 15 days of rain, which may have reduced the contact rate between vectors and hosts, limiting transmission. High levels of precipitation within a short period negatively impact survival of mosquito larvae (e.g., by flushing away egg rafts and larvae, and diluting nutrients) (Koenraadt and Harrington 2008) and may also reduce female host-seeking activity (Roiz et al. 2010), resulting in a negative impact on WNV circulation (Jones et al. 2012, Crowder et al. 2013). Our sample size in July was not sufficient to test the recorded reduction in seroprevalence for significance, but our data provide further indication that infection levels in squirrels may be associated with corresponding infection levels in vectors.
The correspondence between seasonality and level of seroprevalence of squirrels and birds collected in Grant Park could be linked to feeding behavior of mosquitoes circulating in the area. As shown by Kilpatrick et al. (2006), WNV spillover is driven by mosquito feeding preference, which determines contact rates between host and vectors. Our findings highlight the feeding of Atlanta mosquito vectors on mammals, a likely route for human spillover.
Seroprevalence in squirrels collected in our study may also be linked partially to nonvector transmission routes. It has been shown that WNV can be transmitted by oral inoculation or contact in birds (Komar et al. 2003, Bowen and Nemeth 2007). Although the diet of eastern gray squirrels consists primarily of nuts, this species is opportunistic and sometimes preys on nestling and adult passerine birds (Callahan 1993), a behavior that becomes more frequent during summer. Several studies have reported presence of virus in samples of saliva, feces, urine, and skin from infected tree squirrels, which can persist for weeks (Root et al. 2005, 2006, Padgett et al. 2007, Gómez et al. 2008b, Tiawsirisup et al. 2010), suggesting that some squirrels may be infected by contact with infected individuals that are shedding virus (e.g., through grooming, territorial fighting, communal nests).
Model selection showed that age was an important factor predicting seropositivity. Juvenile squirrels are less likely to be seropositive compared with adults. Although this result was not significant, it highlights that seroprevalence in young individuals is affected by the reduced exposure to the virus compared with adults. Low contact with virus of juveniles has been shown in eastern gray squirrels (Gómez et al. 2008a).
We did not detect significant differences in seroprevalence among squirrels collected in different sites, although the distance between trapping locations of our study was higher than the small home range (usually <5 ha) (Don 1983) of squirrels in urban areas, where food availability is high. Furthermore, the distance between trap locations was smaller than the home range of primary and secondary vectors for WNV (Brust 1980, Hamer et al. 2014), indicating that mosquitoes can transmit the virus to the squirrel population of Grant Park without creating spatial heterogeneity. Indeed, evidence of spatial heterogeneity in seroprevalence of squirrels has been detected only when collection sites were distributed at distances greater than mosquito flight range (Gómez et al. 2008a).
In summary, our findings show that the seasonal pattern of squirrel exposure to WNV is similar to that one experienced by birds living in Grant Park. Given the evidence of WNV circulation in their population and their small home range compared to birds, and possibly even to mosquitoes, tree squirrels may serve to identify WNV spillover in mammals at local scale. Furthermore, seroprevalence in tree squirrels may also be used to detect urban areas with high risk for human spillover due to vector feeding behavior.
Footnotes
Acknowledgments
This research was funded by research support from Emory University to U.K. and to E.N.V.S. students. Plaque-reduction neutralization tests were performed using intramural funds provided by Iowa State University.
We thank the students of the Vector Ecology Lab at Emory University who participated in the squirrel trapping and census. We also thank Jamie Allen, Stewart Haddock, Nat Slaughter, Laura Straub, and the rest of the Inman Park Squirrel Census team for their technical support during planning of collections and the squirrel census.
Author Disclosure Statement
No competing financial interests exist.
