Research Article Issue 10 | October 2010 | Vol. 118
The Risk of West Nile Virus Infection Is Associated with Combined Sewer Overflow Streams in Urban Atlanta, Georgia, USA
Gonzalo M. Vazquez-Prokopec1, Jodi L. Vanden Eng2, Rosmarie Kelly3, Daniel G. Mead4, Priti Kolhe5, James Howgate5, Uriel Kitron1,6, Thomas R. Burkot2
1Emory University, Atlanta, Georgia, USA; 2Centers for Disease Control and Prevention, Atlanta, Georgia, USA; 3Georgia Division of Public Health, Atlanta, Georgia, USA; 4University of Georgia, Athens, Georgia, USA; 5Fulton County Department of Health and Wellness, Atlanta, Georgia, USA; 6Fogarty International Center, National Institutes of Health, Bethesda, Maryland, USA
Environ Health Perspect 118:1382-1388 (2010). http://dx.doi.org/10.1289/ehp.1001939 [online 08 June 2010]
Background: At present, the factors favoring transmission and amplification of West Nile Virus (WNV) within urban environments are poorly understood. In urban Atlanta, Georgia, the highly polluted waters of streams affected by combined sewer overflow (CSO) represent significant habitats for the WNV mosquito vector Culex quinquefasciatus. However, their contribution to the risk of WNV infection in humans and birds remains unclear.
Objectives: Our goals were to describe and quantify the spatial distribution of WNV infection in mosquitoes, humans, and corvids, such as blue jays and American crows that are particularly susceptible to WNV infection, and to assess the relationship between WNV infection and proximity to CSO-affected streams in the city of Atlanta, Georgia.
Materials and methods: We applied spatial statistics to human, corvid, and mosquito WNV surveillance data from 2001 through 2007. Multimodel analysis was used to estimate associations of WNV infection in Cx. quinquefasciatus, humans, and dead corvids with selected risk factors including distance to CSO streams and catch basins, land cover, median household income, and housing characteristics.
Results: We found that WNV infection in mosquitoes, corvids, and humans was spatially clustered and statistically associated with CSO-affected streams. WNV infection in Cx. quinquefasciatus was significantly higher in CSO compared with non-CSO streams, and WNV infection rates among humans and corvids were significantly associated with proximity to CSO-affected streams, the extent of tree cover, and median household income.
Conclusions: Our study strongly suggests that CSO-affected streams are significant sources of Cx. quinquefasciatus mosquitoes that may facilitate WNV transmission to humans within urban environments. Our findings may have direct implications for the surveillance and control of WNV in other urban centers that continue to use CSO systems as a waste management practice.
Key words: arbovirus, Culex quinquefasciatus, risk factors, spatial clustering, spatial epidemiology, urban pollution
Address correspondence to G.M. Vazquez-Prokopec, Department of Environmental Studies, Emory University, 400 Dowman Drive, Math and Science Center, 5th floor, Suite E536, Atlanta, GA 30322 USA. Telephone: (404) 712-8728. Fax: (404) 727-4448. E-mail: email@example.com
We gratefully thank A. Brackery for helping with the initial data preparation and geocodification and A. VanNostrand and J. Clennon for helpful comments. We also thank J. Schroeder and three anonymous reviewers for their helpful suggestions for revisions, which greatly improved this article.
This work was supported in part by the Centers for Disease Control and Prevention (CDC) and Emory University (Department of Environmental Studies internal support), Research and Policy for Infectious Disease Dynamics (RAPIDD) program of the Science and Technology Directorate, U.S. Department of Homeland Security, and the Fogarty International Center, National Institutes of Health.
The findings, opinions, and conclusions in this report are those of the authors and do not necessarily represent the views of the CDC or the institutions with which the authors are affiliated.
The authors declare they have no actual or potential competing financial interests.
Received 13 January 2010; accepted 08 June 2010; online 08 June 2010.
West Nile Virus (WNV) is a mosquito-transmitted virus (family Flaviviridae) and a human, equine, and avian neuropathogen (Kramer et al. 2008). The virus is indigenous to the old world (i.e., Africa and Middle East) and is maintained by a bird–mosquito–bird transmission cycle primarily involving Culex sp. mosquitoes; humans, horses, and other mammals are dead-end hosts for the virus (Kramer et al. 2008). Since its emergence in New York City in 1999, WNV has spread over much of North America and the Caribbean and has become a threat to public, equine, and wildlife health (Kramer et al. 2008; Lindsey et al. 2010). Human WNV infections cluster in space and time, with the highest concentrations of cases occurring in urban environments during the summer (Brown et al. 2008; Kramer et al. 2008). Factors that strongly influence WNV transmission dynamics in urban environments throughout the United States include variation in the occurrence and competence of mosquito vectors and vertebrate reservoirs (Kilpatrick et al. 2006a; Reisen et al. 2005, 2008a; Turell et al. 2005), local variation in mosquito feeding patterns (Hamer et al. 2009; Kilpatrick et al. 2006b), bird–herd immunity and infectiousness (Reisen et al. 2008b), WNV genotypes (Bertolotti et al. 2008; Kilpatrick et al. 2008), human population characteristics (Ruiz et al. 2004, 2007), climate, and other environmental variables (Kilpatrick et al. 2008; Soverow et al. 2009).
Information on the distribution and abundance of mosquito vectors in urban environments is paramount for assessing local risks of exposure to WNV. In the United States, members of the Culex pipiens complex (Cx. pipiens and Culex quinquefasciatus), Culex tarsalis, Culex stigmatosoma, Culex salinarius, Culex nigripalpus, and Culex restuans have been implicated as urban vectors of WNV (Hayes et al. 2005; Reisen et al. 2008a; Turell et al. 2005). Most of these species (particularly the Cx. pipiens complex) require organic-rich water for larval development and are commonly found in high numbers in artificial containers, unattended pools, retention ponds, storm drains, catch basins, sewage systems and treatment plants, and pit latrines (Calhoun et al. 2007; Curtis and Feachem 1981; Reisen et al. 2008c).
Combined sewage and street runoff discharges into natural waterways are considered the main sources of urban stream pollution (Paul and Meyer 2001; Tibbetts 2005). According to the U.S. Environmental Protection Agency (EPA), about 850 billion gallons of untreated mixed wastewater and stormwater are discharged annually into U.S. urban waterways, mainly through combined sewer overflow (CSO) systems (Tibbetts 2005; U.S. EPA 2004). In times of low precipitation, CSO systems channel wastewater to a treatment plant before it is discharged into a waterway. However, treatment facilities are by-passed during heavy precipitation or snowmelt, and combined wastewater and stormwater are discharged directly into natural waterways after only minimal chlorine treatment and sieving to remove large physical contaminants. Because of concerns about their potential impacts on both public health and the environment, CSO systems have become the focus of a debate regarding the best techniques to manage growing volumes of sewage and stormwater runoff in many U.S. communities (Tibbetts 2005). About 40 million people in 772 cities from 32 states live in cities with CSO systems; most are located in the mid-Atlantic and midwestern states (U.S. EPA 2004).
An unforeseen effect of CSO discharges into natural streams has been the increased breeding of mosquito species that are considered important vectors of human pathogens. In a study performed in two creeks in the city of Atlanta, Georgia, Calhoun et al. (2007) found that streams that received nutrient pulses from CSO systems had increased numbers of Cx. quinquefasciatus, the main vector of WNV in the area. Other studies have shown that a faster mosquito larval development (Curtis and Feachem 1981), coupled with a reduced larval predation and intraspecific and interspecific competition (Agnew et al. 2000; Beketov and Liess 2007; Mogi and Sota 1996), and a significant tendency of females to lay their eggs in organically enriched waters (Chaves et al. 2009) have contributed to the high Cx. quinquefasciatus productivity in CSO-affected streams. As Cx. quinquefasciatus larval productivity levels in CSO streams seem not to be significantly affected by distance from the source of overflow discharges (Calhoun et al. 2007), the potential impact of such environments as sources of WNV vectors could extend throughout their drainage path, covering a significant proportion of the urban space. However, to our knowledge, an assessment of the potential contribution of CSO-affected streams as sources of WNV infection for humans and birds has never been undertaken.
In this study, our objectives were to describe and quantify the spatial distribution of mosquito abundance and WNV infection in mosquitoes, humans, and birds, specifically corvids (blue jays and American crows), that are considered sentinels of local virus activity (Nielsen and Reisen 2007) in the city of Atlanta, Georgia, from 2001–2007 and to assess the relationship between CSO-affected streams and WNV infection.
Materials and Methods
Study area. From October 2001 through December 2007, the state of Georgia reported 199 human WNV infections and 17 deaths [Centers for Disease Control and Prevention (CDC) 2010]. Of these, 25.1% resided in Fulton County, the most populous county in Georgia (992,472 inhabitants based on 2008 census estimates) and the core of the Atlanta metropolitan area. Fulton County covers 1,385 km² and encompasses populations with a wide range of socioeconomic conditions, from the wealthiest neighborhoods in the state to those with the highest poverty rates in the United States. Seven CSO facilities and the streams associated with them (CSO streams) and about 76 non-CSO streams (excluding tributary creeks) are located within Fulton County.
Data sources and management. Geographic coordinates of each CSO facility and associated streams and data on mosquito abundance and WNV infections in humans, mosquitoes, and bird carcasses from October 2001 through December 2007 were provided by the state of Georgia and the Fulton County Department of Health and Wellness (FCDHW). The protocols for data management, analysis, and dissemination were approved by the CDC Institutional Review Board committee (protocol 990121).
From 21 August 2001 through 31 October 2007, mosquitoes were trapped only in urban environments, and the trapping was concentrated in the city of Atlanta (where most of the population of Fulton County resides). They were collected overnight using paired CO2-baited CDC light traps (Silver 2008) and gravid traps (Reiter 1983) that were placed within 5 m of each other at each sampling location. A total of 455 trap locations were sampled in Fulton County [average number of trap locations per year was 93.7 ± 47.8 (± SD)]. Trap locations were clustered in Fulton County but were randomly distributed within the city of Atlanta [see Supplemental Material, Figure 2 (doi:10.1289/ehp.1001939)]. Every year, mosquito sampling was performed biweekly from 1 May through 31 October by FCDHW staff.
Mosquitoes were classified by species; pooled by date, trap location, species, and trap type; and tested for WNV infection [average number of mosquitoes tested in each group was 13.4 (range, 1–25); see Supplemental Material (doi:10.1289/ehp.1001939) for details on the WNV testing procedures]. To derive WNV infection rates, we used the maximum likelihood (ML) estimation tool (Biggerstaff 2006) to account for variation in the numbers of mosquitoes tested (WNV-positive mosquitoes per 1,000 tested). The intensity of WNV infection was defined as the average ML WNV mosquito infection rate for the period 2001–2007.
We used geocoded data on the 50 human cases of WNV reported in Fulton County from 2001–2007 and human population data from the 2000 U.S. census (U.S. Census Bureau 2000) to estimate annual human WNV infection incidence rates (cases/100,000 persons) for each census tract. Incidence rate estimates for relatively uncommon (or underreported) diseases such as WNV can be imprecise, particularly for small administrative units such as census tracts (Waller and Gotway 2004). Thus, we used spatial Empirical Bayes (EB) smoothing methods that reduce variance by borrowing information from neighboring spatial units. However, because there is a tradeoff between improved precision and the introduction of bias due to the smoothing process (Waller and Gotway 2004), we have reported the results of spatial analyses based on both smoothed (EB) and unsmoothed human WNV incidence rates.
The locations of WNV-positive bird carcasses submitted by citizens to health authorities were geocoded (ArcGis 9.2; ESRI, Redlands, CA). We used data on WNV-positive corvid (i.e., blue jays and American crows) carcasses only, as a consequence of their high WNV mortality and proximity to humans (Nielsen and Reisen 2007). We divided the number of WNV-positive corvid carcasses per census tract by the 2000 population estimate (U.S. Census Bureau 2000) to obtain the number of WNV-positive corvids per 100,000 persons per year in each census tract (i.e., referred to as the WNV-positive corvid death ratio). This metric allowed minimizing bias due to nonrandom variation in the submission of bird carcasses across tracts. We have reported the average WNV-positive corvid death ratios for each tract from August 2001 through December 2004 because bird carcass submissions were highest during this time period.
Environmental factors that might influence mosquito spatial distribution and WNV infection were derived from various data sources. Elevation (in meters above sea level) was derived from the Shuttle Radar Topography Mission (SRTM) data set (90-m pixel size and pixel values describing elevation), generated by the Consultative Group on International Agricultural Research-Consortium for Spatial Information (CGIAR-CSI 2008). Distances (in meters) to the nearest CSO stream and the nearest catch basin (potential Cx. quinquefasciatus urban larval habitats) were estimated by a straight line function in ArcGis 9.2 (expressed as a raster layer with a 30-m pixel size and pixel values indicating the distance of the pixel to the nearest CSO stream or catch basin. The extent (percentage) and variability (range) in tree canopy cover and the percentage of land covered by different land-use classes were obtained from the National Land Cover Database (30-m pixel size and pixel values indicating either percentage of canopy cover or land-cover type) generated from Lansat satellite images collected in 2001 by the Multi-Resolution Land Characteristics Consortium (MRLC 2008).
Predictors of Cx. quinquefasciatus abundance and WNV infection rates (distance to CSO and catch basins, elevation, land use) were estimated at each trap location by averaging values within 1 km of each location (to account for spatial dependence in trap collections due to mosquito flight). Tree canopy range was estimated by calculating the range of canopy cover values within a 1-km radius of each trap. Predictors of human WNV incidence rates and WNV-positive corvid death ratios were estimated at the census tract level. Environmental factors (distance to CSO and catch basins, elevation, land use, and canopy coverage) were averaged per tract, whereas tract estimates of sociodemographic factors (i.e., median household income and percentage of houses built in the 1950s and 1960s) were obtained from the 2000 U.S. census (U.S. Census Bureau 2000). Vegetation, elevation, household income, and age of the houses were identified as significant predictors of WNV human infection in the Greater Chicago area (Ruiz et al. 2004).
Statistical analyses. We used spatial statistics to identify the location and extent of statistically significant clustering (p < 0.05) of high Cx. quinquefasciatus abundance and WNV infection in humans and corvids [see Supplemental Material (doi:10.1289/ehp.1001939) for a detailed description of the statistics].
The spatial distribution of Cx. quinquefasciatus abundance was mapped using a Kernel density function (Silverman 1986) and expressed as mosquitoes per hectare. We used Getis’ Gi*(d) statistic (Ord and Getis 1995) to detect spatial clusters of high Cx. quinquefasciatus abundance per trapping location (the 2001–2007 average number of Cx. quinquefasciatus per trap-night). Significant clustering of Cx. quinquefasciatus WNV infection intensity per trap location (the 2001–2007 average ML WNV infection rate) was also tested by Gi*(d). We used the Bernoulli spatial scan test (Kulldorff 1997) to detect the presence of significant clustering of Cx. quinquefasciatus WNV infection (the presence or absence of WNV infection on a trapping location during 2001–2007). A local Moran’s I local indicator of spatial association (LISA) test (Waller and Gotway 2004) was used to identify significant clustering of human WNV infection rates and WNV-positive corvid death ratios. Statistical inference was based on the difference between observed values and expected values generated from 999 Monte Carlo randomizations.
We assessed relationships of Cx. quinquefasciatus abundance and WNV infection intensity, human EB WNV incidence rates, and WNV-positive corvid death ratios with selected environmental and demographic characteristics (including the average distance to CSO streams and catch basins, land use, percentage of tree cover, median household income, and housing characteristics) using a multimodel selection approach (Burnham and Anderson 2002). Under this analytic framework, we contrasted a set of candidate models (representing different hypotheses) with each other and identified the best model (or the best set of models) based on model fit (Burnham and Anderson 2002), where the best model was the one with the lowest Akaike information criterion (AIC) value. When the lowest AIC value differed from the next best model by < 2 units, we identified a best set of models rather than a single best model (Burnham and Anderson 2002). We estimated the Akaike weight (ωi) for each model as a measure of the probability that a particular model fitted the data better than the alternative set of candidate models (Burnham and Anderson 2002). For each independent variable (j) evaluated, we estimated the sum of Akaike weights (Σωi) as the sum of the ωi from each model in which j was a statistically significant predictor of the outcome (Burnham and Anderson 2002). This metric (bounded between 0 and 1) evaluates the relative importance of each independent variable for predicting the dependent variable (Burnham and Anderson 2002). We contrasted seven models (a full model and six linear combinations of environmental, demographic, and spatial factors) for Cx. quinquefasciatus abundance and WNV infection and five models for human WNV incidence rates and for WNV-positive corvid death ratios.
We performed the statistical analyses using STATA (version 9.0; StataCorp, College Station, TX); for all spatial analyses, we used ClusterSeer2 (TerraSeer, Ann Arbor, MI), GeoDa (http://geodacenter.asu.edu/) and SatScan (http://www.satscan.org/).
Only data from Cx. quinquefasciatus collected from gravid traps (38,120 of the total 71,824 mosquitoes collected) were used in this study [see Supplemental Material (doi:10.1289/ehp.1001939)]. Yearly averages of Cx. quinquefasciatus ML WNV infection rates varied between 2.6 and 16.3 per 1,000 mosquitoes, with no clear temporal pattern (Table 1). Data were available for 321 WNV-positive corvids, which accounted for 98.2% of the WNV-infected dead birds collected in 2001–2006 [see Supplemental Material, Table 1 (doi:10.1289/ehp.1001939)]. The declining temporal trend in WNV-positive corvid death ratios (from 28.0 per 100,000 persons in 2002 to 0.00 per 100,000 in 2006; Table 1) is most likely explained by the decreasing number of dead birds submitted by Atlanta residents to the FCDHW. Of the 50 confirmed cases among Fulton Country residents in 2001–2007, 26 (52%) had symptoms of neuroinvasive WNV. The 2001–2007 average WNV incidence rate among Fulton County residents (0.8 cases per 100,000) was two times higher than the 1999–2008 national average (0.4 cases per 100,000) (Lindsey et al. 2010). With the exception of the years 2001 and 2006, human incidence rates showed little temporal variation, ranging between 0.93 and 1.0 cases per 100,000 (Table 1).
Cx. quinquefasciatus was heterogeneously distributed, with the highest abundances (number per trap-night) and densities (numbers per hectare) concentrated within the boundaries of the city of Atlanta (Figure 1A). Traps with elevated abundance of Cx. quinquefasciatus mosquitoes clustered in three locations within an average radius of 1,305 m (range, 200–2,000) from each trap (Figure 1A) [Gi*(d) > 3.7; p < 0.05]. The three high mosquito abundance clusters encompassed segments of five of the seven CSO streams in Atlanta (Figure 1A). High WNV infection intensity in Cx. quinquefasciatus (WNV-positive mosquitoes/1,000) also clustered in three locations within Atlanta (Figure 1B, Gi*(d) > 3.7; p < 0.05). The largest cluster [average radius = 1,042 m (range, 400–2,000)] surrounded part of the North Avenue CSO creek and encompassed 41 trapping locations (Figure 1B). The remaining two clusters encompassed only two trap locations each, with clustering radii of 320 and 93 m, respectively (Figure 1B). Clusters of WNV infection presence in Cx. quinquefasciatus (presence/absence of WNV infection at a trapping location) also occurred within the city of Atlanta, in proximity [average = 988 m (range, 10–3,118)] to CSO streams (Figure 1B) (Kulldorff Bernoulli test, p < 0.05). K-function analysis (Ripley 1976) confirmed that the observed clustering patterns were not the result of heterogeneous trap distribution, because traps were randomly distributed within the city of Atlanta [Supplemental Material, Figure 2B (doi:10.1289/ehp.1001939)].
(A) Abundance (mosquitoes/trap-night), density distribution (mosquitoes/ha), and local spatial clustering of Cx. quinquefasciatus abundance, 2001–2007. (B) Distribution of WNV infections in Cx. quinquefasciatus and location of clusters of high ML WNV infection intensity [Gi*(d)] and presence (Bernoulli test) in mosquito pools (Cx. quinquefasciatus grouped according to date and location) tested during 2001–2007. Large insets show a detailed view of the city of Atlanta, whereas the upper right inset shows the location of Fulton County.
The proportion of WNV-infected Cx. quinquefasciatus pools (mosquitoes grouped by date and location) was significantly higher within 1 km of a CSO stream (125/1,215) than within 1 km of a non-CSO stream (77/1,044) (p = 0.027). Furthermore, 9 of 12 (75%) trapping locations where WNV-infected Cx. quinquefasciatus mosquitoes were detected during early spring were located within 1 km of a CSO stream. The pattern of ML WNV infection clustering in Cx. quinquefasciatus was consistent across years, with most positive pools located in close proximity to CSO streams [see Supplemental Material, Figure 3 (doi:10.1289/ehp.1001939)].
A regression model that included the average distance between each trapping location to the nearest CSO stream and the range of tree canopy cover within 1 km of each trap location was best supported by the data (ωi, 65.8%) (Table 2). Distance to CSO streams was the best predictor of the abundance of Cx. quinquefasciatus mosquitoes (Σωi = 100%, by being statistically significant in the best models), followed by tree canopy coverage range (Σωi = 66.9%) (Table 2). The remaining environmental factors were not statistically significant. Traps with a high density of Cx. quinquefasciatus (> 20/trap-night) were located at an average distance of 1,156 m (95% CI, 552–1,759 m) from a CSO stream and had an average tree canopy range of 98% (95% CI, 95–99%). Tree canopy range can be interpreted as the degree of canopy variability around a trap location, as high values include highly forested and deforested areas close to each other (within 1 km).
Summary of multiple logistic regression models used to evaluate the average abundance of Cx. quinquefasciatus (mosquitoes/trap-night) and WNV infection presence (presence/absence of WNV-positive pools in a trapping location), urban Atlanta, Georgia, USA, 2001–2007.
The presence of WNV infection in Cx. quinquefasciatus was also predicted by the average distance to CSO streams and by the canopy cover range (Table 2). The two best candidate models included distance to CSO streams as the sole significant predictor of Cx. quinquefasciatus WNV infection presence. Distance to CSO streams (Σωi, 80.2%) followed by canopy cover range (Σωi, 15.2%) were the most important predictors in these models.
Human WNV incidence rates (EB) and WNV-positive corvid death ratios also clustered within the city of Atlanta (Figure 2). One of the two clusters of high human WNV incidence (identified by the Local Moran’s I LISA test) was located near three of the seven CSO streams (Figure 2A). This cluster included six census tracts and six reported cases and had an average annual incidence (± SD) of 48.6 ± 37.4 cases per 100,000 persons. The second (northern) cluster of human WNV incidence was not directly associated with any CSO stream (Figure 2A). This cluster included 10 tracts and 10 cases and had an average annual incidence of 42.6 ± 31.3 cases per 100,000 persons. In contrast, Atlanta census tracts outside the clustering areas had an average annual incidence of 7.5 ± 11.0 human WNV cases per 100,000 persons. Clustering zones were similar but encompassed fewer tracts when we repeated the analysis using unsmoothed (versus EB) human WNV incidence rates [see Supplemental Material, Figure 4 (doi:10.1289/ehp.1001939)]. Spatial clustering of WNV-positive corvid death ratios overlapped with human WNV incidence clusters and Cx. quinquefasciatus infection clusters in the east and north of Atlanta (Figure 2B). Four of the seven CSO streams were partly encompassed by WNV-positive dead corvid ratio clusters (Figure 2B).
Distribution and spatial clustering of (A) EB-smoothed WNV human incidence rate estimates (cases/100,000 persons) and (B) WNV-positive corvid death ratios (number of dead corvids/100,000 persons) in Fulton County. Inset shows a detailed view of the city of Atlanta.
The human WNV incidence rate per census tract (2001–2007, log10 + 1 transformed EB estimates) was best predicted by a model that included the average distance from each tract to the nearest CSO stream, mean tree canopy cover, the proportion of houses built in the 1950s–1960s, median household income, and number of WNV-positive dead corvids identified during 2001–2004 (Table 3) (R2 = 0.21; p < 0.001; ωi = 81.7%). When all variables with p < 0.07 were considered, distance to CSO streams was the best predictor of human WNV incidence (Σωi = 100%), followed by the proportion of houses built in the 1950s–1960s (Σωi = 94.0%) and the median household income (Σωi = 92.0%) (Table 3). The remaining factors were poorly supported by the data. Similarly, log10 + 1 transformed WNV-positive corvid death ratios were best predicted by a model including the average distance to the nearest CSO stream, mean tree canopy cover, proportion of houses built in the 1950s–1960s, and median household income (Table 3) (R2 = 0.31; p < 0.001; ωi = 91.3%). Again, mean distance to CSO streams and mean tree cover (both Σωi = 100%), followed by median household income (Σωi = 91.0%), were the best predictors of high WNV-positive corvid death ratios (Table 3).
The effect of sewage-affected waterways on mosquito dynamics and the pathogens they transmit have been the focus of research in urban environments throughout the world (Curtis and Feachem 1981; El Esnawy and Nagwa 2001; Mogi and Okazawa 1990), but have received limited attention in the United States, even though approximately 746 urban centers are affected by CSOs (U.S. EPA 2004). Surprisingly, the emergence of WNV and other mosquito-transmitted pathogens as a consequence of combined sewer systems is not considered in the 2004 report of the U.S. EPA as a potential human health impact of CSO systems (U.S. EPA 2004). Previous research in the city of Atlanta, Georgia, showed that CSO-affected streams have higher Cx. quinquefasciatus populations and enhance oviposition of female Cx. quinquefasciatus compared with non-CSO streams, where Cx. quinquefasciatus is rarely found (Calhoun et al. 2007; Chaves et al. 2009). Our study builds on that knowledge and shows that WNV infection in Cx. quinquefasciatus is significantly higher in CSO streams than in non-CSO streams and that WNV incidence rates in both humans and corvids are significantly associated with proximity to CSO-affected streams.
According to the U.S. EPA, there is limited information about the human health impacts of CSO systems (Tibbetts 2005; U.S. EPA 2004). Most of the available data focus on the effects of CSO-associated exposures to microbial pathogens (e.g., bacteria, virus, and parasites), toxics (e.g., heavy metals), and biologically active chemicals (e.g., hormones) on human health (Tibbetts 2005; U.S. EPA 2004). Our study shows that WNV infection is significantly higher in close proximity to CSO-affected streams. Similarly, proximity to a sewage plant and its effluent was identified as the only predictor of high WNV human seroprevalence in the city of Hashimiah, Jordan (Batieha et al. 2000). Assuming a Cx. quinquefasciatus urban flight range of approximately 1.0 km (Reisen et al. 1991), the potential influence of CSO streams as sources of WNV infected vectors could extend over significant portions of the urban space. In Atlanta, the area potentially affected by flight-dispersed mosquitoes originating from CSO streams could encompass 23% (78.3/341.1 km2) of the city. These findings may have direct implications for U.S. urban centers that still rely on CSO systems as a primary waste management practice, particularly those located in mid-Atlantic and midwestern states affected by WNV epidemics (Brown et al. 2008).
Spatial clustering of WNV infection in urban environments occurs in areas where favorable larval breeding conditions, competent reservoir hosts, and opportunities for human exposure overlap (Hamer et al. 2009; Nielsen et al. 2008; Ruiz et al. 2004). In Atlanta, such conditions occurred in close proximity to CSO streams that also were characterized by a larger proportion and range of tree canopy cover, low median incomes, and a large proportion of houses built in the 1950s–1960s. Urban tree cover could have facilitated contact between Cx. quinquefasciatus and an abundant bird population (Hamer et al. 2009) and mosquito dispersal away from CSO streams. Because we averaged collections across years, our measures of Cx. quinquefasciatus clustering cannot be considered surrogates of mosquito dispersal. Rather, we consider such distances as measures of the extent of high abundance of Cx. quinquefasciatus populations. On the other hand, people residing in low-income 1950s–1960s housing areas may be at increased risk of exposure to evening-biting Cx. quinquefasciatus compared with residents of high-income areas where air conditioning and protective behaviors such as the use of mosquito repellent or the active avoidance of mosquitoes may be higher. The CSO streams located in the north of Atlanta (one of the wealthiest sectors of the city) overlapped with clusters of WNV infection in mosquitoes and corvids but not human infections, supporting the potential for differential exposure of the human population in such areas. Conversely, in southeast Atlanta (a mid- to low-income area) clustering of WNV infection in humans, corvids, and mosquitoes overlapped. Located within 2.5 km of two CSO streams and connected to them by forest corridors, this area includes a cemetery, recreational areas (e.g., Grant Park, Atlanta Zoo), and commercial and residential spaces, and it has an estimated human WNV incidence rate that is 6.5 times higher than the Atlanta average. Understanding sociodemographic and environmental determinants of virus amplification and transmission dynamics in such transmission hot spots may help explain WNV infection persistence and circulation within urban environments (Bertolotti et al. 2008; Hamer et al. 2009; Ruiz et al. 2004).
Human WNV infection is considered endemic in Atlanta, but the mechanisms explaining WNV persistence in this area are still unclear. The overwintering of the virus in hibernating mosquitoes has been proposed as a mechanism of WNV persistence (Reisen et al. 2006). In the present study, WNV-infected Cx. quinquefasciatus mosquitoes identified in early spring were located in close proximity to CSO streams and facilities, which may have been a consequence of the greater abundance of mosquitoes in CSO versus non-CSO streams. In addition, CSO tunnels and exit pipes may serve as ideal mosquito (and virus) overwintering habitats. Our recent finding of a WNV-infected Cx. quinquefasciatus pool inside CSO tunnels (Vazquez-Prokopec GM, unpublished data) supports such a hypothesis.
In 1999, the city of Atlanta was found to be in violation of both the federal Clean Water Act and the Georgia Water Quality Control Act because of CSO discharges into urban waterways (U.S. EPA 1999). As a part of a major settlement between the city and U.S. EPA, an underground reservoir system to hold excess effluent during heavy precipitation was constructed (at an estimated cost of $3.9 billion) and completed in 2008. The new CSO system was designed to minimize pollution by reducing the number of overflows to < 10. However, this system could provide an ideal environment for Cx. quinquefasciatus mosquitoes and WNV amplification in CSO streams: infrequent flooding with high organic content effluent into streams, a food-rich and predator-free larval habitat, and abundant forest corridors and avian hosts. Further monitoring of mosquito populations in CSO streams in Atlanta is needed to determine if the remediation system will enhance the growth of mosquito populations and further increase the risk of WNV transmission to humans.
We conducted an integrated analysis of information on WNV infections in mosquitoes, dead corvids, and humans and found a strong association between WNV infection and proximity to CSO-affected urban streams. Our findings are consistent with previous research indicating that CSO streams represent major mosquito-breeding habitats and support the need to expand this research to other urban centers that still include CSO systems as part of their wastewater management strategy.
- Agnew P, Haussy C, Michalakis Y. 2000. Effects of density and larval competition on selected life history traits of Culex pipiens quinquefasciatus (Diptera: Culicidae). J Med Entomol 37:732–735.
- Batieha A, Saliba EK, Graham R, Mohareb E, Hijazi Y, Wijeyaratne P. 2000. Seroprevalence of West Nile, Rift Valley, and Sandfly arboviruses in Hashimiah, Jordan. Emerg Infect Dis 6:358–362.
- Beketov MA, Liess M. 2007. Predation risk perception and food scarcity induce alterations of life-cycle traits of the mosquito Culex pipiens. Ecol Entomol 32:405–410.
- Bertolotti L, Kitron UD, Walker ED, Ruiz MO, Brawn JD, Loss SR, et al. 2008. Fine-scale genetic variation and evolution of West Nile Virus in a transmission “hot spot” in suburban Chicago, USA. Virology 374:381–389.
- Biggerstaff BJ. 2006. PooledInfRate, Version 3.0: A Microsoft Excel Add-In to Compute Prevalence Estimates from Pooled Samples. Fort Collins, CO:Centers for Disease Control and Prevention.
- Brown HE, Childs JE, Diuk-Wasser MA, Fish D. 2008. Ecological factors associated with West Nile virus transmission, northeastern United States. Emerg Infect Dis 14:1539–1545.
- Burnham KP, Anderson DR. 2002. Model Selection and Multimodel Inference: A Practical Information-Theoretic Approach. 2nd ed. New York:Springer.
- Calhoun LM, Avery M, Jones L, Gunarto K, King R, Roberts J, et al. 2007. Combined sewage overflows (CSO) are major urban breeding sites for Culex quinquefasciatus in Atlanta, Georgia. Am J Trop Med Hyg 77:478–484.
- CDC (Centers for Disease Control and Prevention). 2010. West Nile Virus Statistics, Surveillance, and Control Archive. Available: http://www.cdc.gov/ncidod/dvbid/westnile/surv&control_archive.htm [accessed 25 August 2010].
- Chaves LF, Keogh CL, Vazquez-Prokopec GM, Kitron U. 2009. Combined sewage overflow enhances oviposition of Culex quinquefasciatus (Diptera:Culicidae) in urban areas. J Med Entomol 46:220–226.
- Consultative Group on International Agricultural Research – Center for Consortium for Spatial Information (CGIAR-CSI). 2008. Level 4 SRTM global elevation data. Available: ftp://srtm.csi.cgiar.org/SRTM_V41/SRTM_Data_GeoTiff/srtm_20_06.zip [accessed 25 August 2010].
- Curtis CF, Feachem RG. 1981. Sanitation and Culex pipiens mosquitoes: a brief review. J Trop Med Hyg (84):17–25.
- El-Esnawy Nagwa A. 2001. Infection by certain arboviruses among workers potentially at risk of infection. J Egypt Public Health Assoc 76:169–182.
- Hamer GL, Kitron UD, Goldberg TL, Brawn JD, Loss SR, Ruiz MO, et al. 2009. Host selection by Culex pipiens mosquitoes and West Nile virus amplification. Am J Trop Med Hyg 80:268–278.
- Hayes EB, Komar N, Nasci RS, Montgomery SP, O’Leary DR, Campbell GL. 2005. Epidemiology and transmission dynamics of West Nile virus disease. Emerg Infect Dis 11:1167–1173.
- Kilpatrick AM, Daszak P, Jones MJ, Marra PP, Kramer LD. 2006. Host heterogeneity dominates West Nile virus transmission. Proc Biol Sci 273:2327–2333.
- Kilpatrick AM, Kramer LD, Jones MJ, Marra PP, Daszak P. 2006. West Nile virus epidemics in North America are driven by shifts in mosquito feeding behavior. PLoS Biol 4:e82; doi:10.1371/journal.pbio.0040082.
- Kilpatrick AM, Meola MA, Moudy RM, Kramer LD. 2008. Temperature, viral genetics, and the transmission of West Nile virus by Culex pipiens mosquitoes. PLoS Pathog 4:e1000092; doi:10.1371/journal.ppat.1000092.
- Kramer LD, Styer LM, Ebel GD. 2008. A global perspective on the epidemiology of West Nile virus. Annu Rev Entomol 53:61–81.
- Kulldorff M. 1997. A spatial scan statistic. Commun Stat Theory Methods 26:1481–1496.
- Lindsey NP, Staples JE, Lehman JA, Fischer M, Centers for Disease Control and Prevention (CDC). 2010. Surveillance for human West Nile virus disease—United States, 1999–2008. MMWR Surveill Summ 59:1–17.
- Mogi M, Okazawa T. 1990. Factors influencing development and survival of Culex pipiens pallens larvae (Diptera, Culicidae) in polluted urban creeks. Res Popul Ecol 32:135–149.
- Mogi M, Sota T. 1996. Physical and biological attributes of water channels utilized by Culex pipiens pallens immatures in Saga City, southwest Japan. J Am Mosq Control Assoc 12:206–214.
- Multi-Resolution Land-cover Consortium (MRLC). 2008. National Land Cover Database. Available: http://gisdata.usgs.net/website/MRLC/viewer.htm [accessed 25 August 2010].
- Nielsen CF, Armijos MV, Wheeler S, Carpenter TE, Boyce WM, Kelley K, et al. 2008. Risk factors associated with human infection during the 2006 West Nile virus outbreak in Davis, a residential community in northern California. Am J Trop Med Hyg 78:53–62.
- Nielsen CF, Reisen WK. 2007. West Nile virus-infected dead corvids increase the risk of infection in Culex mosquitoes (Diptera:Culicidae) in domestic landscapes. J Med Entomol 44:1067–1073.
- Ord JK, Getis A. 1995. Local spatial autocorrelation statistics. Distributional issues and an application. Geogr Anal 27:286–306.
- Paul MJ, Meyer JL. 2001. Streams in the urban landscape. Annu Rev Ecol Syst 32:333–365.
- Reisen WK, Barker CM, Fang Y, Martinez VM. 2008. Does variation in Culex (Diptera:Culicidae) vector competence enable outbreaks of West Nile virus in California? J Med Entomol 45:1126–1138.
- Reisen WK, Fang Y, Lothrop HD, Martinez VM, Wilson J, O’Connor P, et al. 2006. Overwintering of West Nile virus in Southern California. J Med Entomol 43:344–355.
- Reisen WK, Fang Y, Martinez VM. 2005. Avian host and mosquito (Diptera:Culicidae) vector competence determine the efficiency of West Nile and St. Louis encephalitis virus transmission. J Med Entomol 42:367–375.
- Reisen WK, Lothrop HD, Wheeler SS, Kennsington M, Gutierrez A, Fang Y, et al. 2008. Persistent West Nile virus transmission and the apparent displacement St. Louis encephalitis virus in southeastern California, 2003–2006. J Med Entomol 45:494–508.
- Reisen WK, Milby MM, Meyer RP, Pfuntner AR, Spoehel J, Hazelrigg JE, et al. 1991. Mark-release-recapture studies with Culex mosquitoes (Diptera:Culicidae) in Southern California. J Med Entomol 28:: 357–371.
- Reisen WK, Takahashi RM, Carroll BD, Quiring R. 2008. Delinquent mortgages, neglected swimming pools, and West Nile virus, California. Emerg Infect Dis 14:1747–1749.
- Reiter P. 1983. A portable, battery-powered trap for collecting gravid Culex mosquitoes. Mosq News 43:496–498.
- Ripley BD. 1976. Second order analysis of stationary point processes. J Appl Prob 13:255–266.
- Ruiz MO, Tedesco C, McTighe TJ, Austin C, Kitron U. 2004. Environmental and social determinants of human risk during a West Nile virus outbreak in the greater Chicago area, 2002. Int J Health Geogr 3:8.doi:10.1186/1476-072X-3-8
- Ruiz MO, Walker ED, Foster ES, Haramis LD, Kitron UD. 2007. Association of West Nile virus illness and urban landscapes in Chicago and Detroit. Int J Health Geogr 6:10; doi:10.1186/1476-072X-6-10.
- Silver JB. 2008. Sampling adults with light-traps. In: Mosquito Ecology: Field Sampling Methods. New York:Springer, 845–946.
- Silverman BW. 1986. Density Estimation for Statistics and Data Analysis. London; New York:Chapman and Hall.
- Soverow JE, Wellenius GA, Fisman DN, Mittleman MA. 2009. Infectious disease in a warming world: how weather influenced West Nile virus in the United States (2001–2005). Environ Health Perspect 117:1049–1052.
- Tibbetts J. 2005. Combined sewer systems: down, dirty, and out of date. Environ Health Perspect 113:A464–A467.
- Turell MJ, Dohm DJ, Sardelis MR, Oguinn ML, Andreadis TG, Blow JA. 2005. An update on the potential of North American mosquitoes (Diptera:Culicidae) to transmit West Nile Virus. J Med Entomol 42:57–62.
- U.S. Census Bureau. 2000. Census of Population and Housing, Summary File I, American Factfinder. Available: http://factfinder.census.gov [accessed 5 August 2010].
- U.S. EPA (U.S. Environmental Protection Agency). 1999. U.S. Georgia and Atlanta Reach Settlement to Fix City’s Aging Sewer System. Available: http://yosemite.epa.gov/opa/admpress.nsf/016bcfb1deb9fecd85256aca005d74df/1d9bf67474410410852567bd0073d60c!OpenDocument [accessed 5 August 2010].
- U.S. EPA (U.S. Environmental Protection Agency). 2004. Report to Congress. Impacts and Control of CSOs and SSOs. Washington, DC:U.S. Environmental Protection Agency.
- Waller LA, Gotway CA. 2004. Applied Spatial Statistics for Public Health Data. Hoboken, NJ:John Wiley & Sons.
Come visit with EHP’s Science Editor, Jane Schroeder, at the 2016 ISEE Meeting in Rome, Italy, from 1–4 September 2016. This is a great opportunity to meet a member of our team, and to learn more about the journal. A number of our Associate Editors will be in attendance as well. Jane is also co-hosting a half-day writing and publishing workshop on 4 September, following the conference.
EHP is pleased to present the abstracts from the 28th annual meeting of the International Society for Environmental Epidemiology (ISEE), held in Rome, Italy, 1–4 September 2016, and hosted by the Department of Epidemiology Lazio Regional Health Service, ASL Roma 1, and the Italian Epidemiological Association. The focus of this year’s conference is current and future challenges in exposure assessment, study design, and data analyses.
Featured Children’s Health
Birgit Claus Henn, Adrienne S. Ettinger, Marianne R. Hopkins, Rebecca Jim, Chitra Amarasiriwardena, David C. Christiani, Brent A. Coull, David C. Bellinger, and Robert O. Wright
Diane Gilbert-Diamond, Jennifer A. Emond, Emily R. Baker, Susan A. Korrick, and Margaret R. Karagas
Benjamin B. Green, Margaret R. Karagas, Tracy Punshon, Brian P. Jackson, David J. Robbins, E. Andres Houseman, and Carmen J. Marsit