ALL Metrics
-
Views
-
Downloads
Get PDF
Get XML
Cite
Export
Track
Research Article

Modeling enzootic raccoon rabies from land use patterns - Georgia (USA) 2006-2010

[version 1; peer review: 2 approved]
PUBLISHED 27 Dec 2013
Author details Author details
OPEN PEER REVIEW
REVIEWER STATUS

This article is included in the Neglected Tropical Diseases collection.

Abstract

We analyzed how land-use patterns and changes in urbanization influence reported rabid raccoons in Georgia from 2006 - 2010.  Using Geographical Information Systems and rabies surveillance data, multivariate analysis was conducted on 15 land-use variables that included natural topography, agricultural development, and urbanization to model positive raccoon rabies cases while controlling for potential raccoon submission bias associated with higher human population densities.  Low intensity residential development was positively associated with reported rabid raccoons while a negative association was found with evergreen forest.  Evergreen forests may offer a barrier effect where resources are low and raccoon populations are not supported.  Areas with pure stands of upland evergreen forest might be utilized in baiting strategies for oral rabies vaccination programs where fewer or no baits may be needed.  Their use as a barrier should be considered carefully in a cost-effective strategy for oral rabies vaccination (ORV) programs to contain the western spread of this important zoonotic disease.

Keywords

enzootic, environmental heterogeneity, GIS, land-use, raccoon, urbanization

Introduction

Here we use Geographical Information Systems (GIS) analysis to exhibit how both natural topography and urbanization potentially affect the dynamics of a zoonotic disease in its reservoir species. Today, the primary reservoir for rabies in Georgia, and, in fact, the entire eastern part of the USA, is the raccoon (Procyon lotor). Although dogs previously served as the main reservoir for rabies throughout the USA, dog rabies virus variants in the USA have been eliminated through persistent vaccination programs and policies1. Rabies virus has evolved to adapt to a wildlife meso-carnivore in a population that can be defined by geographical regions where transmission is sustained (Figure 1)2. Over the second half of the last century, the eastern raccoon rabies variant has spread northward and, in Georgia, the raccoon rabies enzootic has been established for over 50 years3.

ac25883b-23de-4af9-a2ca-df559a09f817_figure1.gif

Figure 1. Reservoir species for rabies across the USA.

Infection from the reservoir species to another mammalian host can, and still does occur, but the non-reservoir animal most likely becomes a dead-end host4. Cases of rabies in Georgia in non-reservoir wildlife, such as foxes, coyotes and skunks, and in domesticated animals, such as dogs, cats, and farm animals, can almost always be traced either to the raccoon variant, or to one of the circulating variants among several bat species, through DNA molecular typing5. Using Georgia Department of Public Health laboratory data on raccoon specimens submitted for rabies testing, and the US Geological Survey Land Cover Database, we report a model to predict rabid raccoon cases from land use patterns.

Raccoon rabies epidemiology

It is well documented that higher densities of raccoon populations exist in urbanized areas than in more natural habitats, and, that the highest densities are found at the urban-rural interface69. An opportunistic omnivore, eating anything from crayfish in creeks to acorns, bird eggs, garden vegetables and even garbage, raccoons have been essentially subsidized by urbanization by the creation of an edge ecosystem. An edge ecosystem is created when habitat is broken up. This can occur naturally, such as disturbance by fire, or can be part of the geomorphology, such as the demarcation between uplands and wetlands. Or it can occur when development creates a new transition between forest and open space. Although some of the open space is replaced with impervious surfaces, i.e. roads and rooftops, higher raccoon populations can be sustained because they have increased access to food resources, and fewer predators, while still using the adjacent natural habitat for protective cover10.

This raccoon ecological theory, and the way it relates to land use patterns, has been supported in numerous studies where raccoon rabies epidemics occur during the initial introduction of rabies into naïve raccoon populations. These studies have described “sensational” epizootics so fulminating that natural barriers such as rivers can only slow them down11,12. Once the infection becomes enzootic in a raccoon population, large outbreaks rarely occur and the disease becomes incessant (Figure 2)13.

ac25883b-23de-4af9-a2ca-df559a09f817_figure2.gif

Figure 2. Schematic of the temporal stages of a typical county in New York during the raccoon rabies epizootic from 1992 to 2000.

Stage a: Pre-raccoon variant – rare cases in raccoons; cases of raccoon rabies might spill over from other wildlife. b1: initial epizootic of raccoon variant rabies as it moves through the county. c1: sequential epizootic. d1: sequential inter-epizootics. (Adapted from Gordon et al., 200413).

Mitigating the passive public health database

The number of raccoons submitted for rabies testing listed by county and month was acquired from the Georgia Department of Public Health for a five-year period from 2006 to 2010. This database only represents results from raccoons involved in human or domestic animal exposure incidents that are submitted to the Georgia Public Health Laboratories (GPHL). Testing is done solely to inform medical professionals how to proceed with rabies prophylaxis for potential human exposure cases. The presence of rabies virus in the raccoon specimens is tested by direct fluorescent-antibody staining of brain tissue. [http://www.cdc.gov/rabies/diagnosis/direct_fluorescent_antibody.html].

The mosaic of urban and suburban development poses challenges to this type of surveillance of the raccoon rabies reservoir because higher human population densities bias reporting to passive public health surveillance systems (Figure 3); detection of positive cases in developed areas is likely increased due to the probability of more human-raccoon and domestic animal-raccoon interactions14. Due to potential confidentiality issues concerning the person exposed, this data is aggregated and reported at the county level. The number of raccoons submitted for testing closely parallels the human population densities.

ac25883b-23de-4af9-a2ca-df559a09f817_figure3.gif

Figure 3. The population densities of 159 Georgia counties by quartile from the 2010 census as compared to the number of raccoons submitted per county for rabies testing by quartile from 2006–2010.

Prior to 2006, only positive cases were reported by GPHL. Because the 2006–2010 GPHL dataset reported both positive and negative submissions, there was opportunity to possibly mitigate the surveillance bias to some extent by incorporating a standardizing variable. The Getis-Ord Gi statistical tool in GIS (ESRI ArcMap 10) was used to exhibit spatial clustering of counties with higher or lower than expected number of positive cases, positive cases per person per square kilometer, and the number of submissions per person per square kilometer (Figure 4)15.

ac25883b-23de-4af9-a2ca-df559a09f817_figure4.gif

Figure 4. Significant clustering of positive cases occurs mostly in the metro Atlanta area where there are also high submission rates.

The clustering of positive cases is mostly away from major population centers when analyzed as a rate of population density. Clustering of submissions/density closely parallels positive cases/density.

When the number of positive cases is put at a rate of population density, the clustering disperses away from the metropolitan areas. However, because our model counted the number of positive cases as the outcome variable, we wanted to use some form of the submissions data as the “weighted” standardizing variable on the right side of the regression model to normalize counties.

Table 1 summarizes the human demographic and GPHL testing data. There were forty-six counties that were inclusive in the upper quartile for positive cases. Forty counties were inclusive in the upper quartile for the number of submissions for testing, as well as all other independent land use variables. The high odds ratio for both submissions, and population density, confirms that as there is more testing, there are more positive cases; it is a matter of more intense surveillance, and does not necessarily reflect the pattern of the rabies reservoir in the raccoon. However, when density is used as the denominator to analyze the effect of rate of submissions on disease predictability, the odds ration falls (3.1). The drop in the odds ratio, though still significant, could reflect two possibilities: that there truly are higher raccoon densities in more populated areas and/or that the database bias has been mitigated to some extent.

Table 1. Human demographic data, number of submissions for testing, and positive raccoon cases from 2006–2010 for 159 Georgia counties.

Sub/density is the weighted variable used in the regression model. High enzootic is the number of counties in the upper quartile of positive cases (≥ 8).

EnzooticEnzootic
VariableMedianUpper quartileRangehighlowOR95% CIp value
Densitya24.9583.2–986.12317 5.6 2.6–12.3 <0.001
Submissions8150–148355 68.7 22–211 <0.001
Positive480–50 --- --- --- --- ---
Sub/densityb0.2380.4320–1.661921 3.1 1.5–6.6 0.003

apeople per km2

bnumber of submitted raccoons for testing per person per km2

OR: odds ration

CI: confidence interval

Land use/land cover and human demographic analysis

Using GIS, we extracted fifteen land use/land cover categories from the National Land Cover Database 2006 [http://www.mrlc.gov/nlcd06_leg.php]. Variable construction was based on the method of Jones et al. 200316, to define low intensity residential, high intensity residential and commercial land use. These three development variables are determined from measuring the impervious surface, and constitute less than 50% (combining “open space” and “low intensity”), 50 – 79% (“medium intensity”), and 80% or greater impervious surface area (“high intensity”), respectively. The agricultural variable was also constructed on the method of Jones et al. 2003, and is a combination of pasture/hay and cultivated crops. Both, demographic data, and land area of each county, were obtained from the 2010 US Census [http://quickfacts.census.gov/qfd/states/13000.html].

There were 1011 rabid raccoons from a total of 2064 raccoon specimens submitted from 150 counties over the five year period. A univariate analysis was conducted for all variables, including land-use data, human population density, and the number of rabid raccoons (outcome variable), by dichotomizing at the upper quartile to determine crude odds ratios (OR).

Urbanization variables (“Low intensity res”, “High intensity res”, and “Com/Ind/Trans”) had the highest crude odds ratios among all the land use variables examined (Table 2). Shrub and the woody wetlands had similar protective odds ratios (0.2–0.3; CI 0.1–0.9) while evergreen forest had the lowest crude odds ratio at 0.04 (CI 0.006–0.3).

Table 2. Crude odds ratios for land-use patterns associated with relative risk of high enzootic raccoon rabies (≥8 positive cases) within 159 Georgia counties.

EnzooticEnzootic
VariableMedianaUpper quartileaRangeahighlowOR95% CIp value
Open water0.791.660.4–29.614261.50.7–3.10.329
Low intensity res 58.1982.8–44.52614 9.2 4.1–20.6 <0.001
High intensity res 0.280.700.1–8.02515 7.8 3.5–17.2 <0.001
Com/Ind/Trans 0.72.3.03–4.82317 5.6 2.6–12.3 <0.001
Barren 0.220.57.08–1.11723 2.3 1.1–4.9 0.029
Shrub 2.844.720.8–9.1535 0.3 0.1–0.7 0.008
Deciduous forest21.3530.174.1–46.113271.30.6–2.70.566
Evergreen forest 21.5431.0511.5–12.6139 0.4 .006–0.3 <0.001
Mixed forest3.025.081.1–9.610300.80.3–1.70.528
Agricultural16.7225.329.0–29.716242.00.9–4.20.075
Grasslands6.158.593.9–4.111290.90.4–2.00.818
Woody wetlands 6.5316.032.9–29.3436 0.2 0.07–0.6 0.002
Herbaceous wetlands0.271.21.01–23.76340.30.1–0.90.025

apercentage of land use

OR: odds ratio

CI: confidence interval

The independent variables with a p-value less than 0.1 were selected to build a regression model for predicting the number of rabid raccoons. Variables were eliminated through a manual backwards stepwise regression. Although the selection of independent variables (percentage land-use) to begin the model was determined with dichotomized data (odds ratios), the regression model predicted the number of rabid cases from continuous data. Potential autocorrelation between variables in the model were considered in the backwards stepwise regression by first removing co-variables with a correlation coefficient greater than 0.6; the co-variable with the highest p-value was eliminated. Because of high variance among the counties, including zeros, in the number of rabid raccoons, both Poisson and negative binomial distributions were evaluated for the regression analysis17.

A negative binomial distribution was chosen for the final regression model after comparing the goodness of fit statistics (AIC=795, deviance/degrees of freedom=1.1, Pearson chi-square/degrees of freedom=0.96) to a model using the Poisson distribution (961, 3.2, and 3.2, respectively). The final negative binomial regression model included:

+ 1.9(Sub/density)[p<0.001] – 1.8(Evergreen forest)[p=0.02]– 52.2(Com/Ind/Trans)[p<0.001] + 11.5(Low intensity res)[p<0.001] + 43.5(Barren)[p=0.019] as predictors of total positive raccoon rabies cases.

Conclusion

We have presented one example of how anthropogenic changes to the landscape impacts the natural ecology of the earth. This study supports the hypothesis that the majority of rabid raccoons are associated with low intensity urbanized areas. If raccoon populations were not “subsidized” and the raccoon densities were allowed to remain at normal carrying capacities, rabies as a reservoir in raccoons might possibly eradicate itself. Our finding that upland evergreen forests are inversely associated with the presence of rabid raccoons supports this possibility too, and, coincides with epidemiological theory that a population threshold is needed to maintain transmission of any infectious disease18. For example, pine forests that are managed for timber production generally have lower raccoon utilization as compared to other forested habitats19. Therefore, compelling evidence exists that the pure stands of managed pine forests in western and southern Alabama have been attributed as being the major barrier to the spread of the raccoon enzootic further west into Mississippi20.

The finding that commercial development was a positive predictor of raccoon rabies in the univariate analysis but had a negative coefficient in the final model suggests that, alone, commercial development is a co-associate of the low intensity urbanization phenomenon. Since the commercial development variable remained in the regression model with a negative coefficient however, it seems that as mixed land patterns give way to highly impacted landscapes, with high levels of impervious surface, there reaches a point at which protective cover falls below a minimal threshold.

The finding that barren land-cover is a predictor of positive raccoon rabies cases has no precedence in the literature, and contributed very little to the model, even with a high coefficient, due to the low proportion of this land type within our entire study area (0.34%). Barren land can mean rock that is exposed either naturally, or from strip mines and gravel pits. There seems to be no ecological reason for these land types to be associated with raccoons. However, “barren” also consists of land with less than 15% vegetative cover that has been temporarily cleared in preparation for development (future urbanization) and can include landfills (food resources) [http://www.hq.nasa.gov/iwgsdi/Barren_Land.html].

The submissions/population density variable was an attempt to normalize all counties in the multivariate model. We added this variable into the regression equation because the crude odds ratio for submissions/population (alone) was higher (4.1). The lower odds ratio using density as the denominator indicates that testing bias is better controlled. The odds of a county being in the upper quartile of percent evergreen forest coverage and having submitted fewer raccoon specimens for testing per person per square kilometer (the standardized variable) as compared to any other county is 1.1 (CI: 0.5–2.6 CI). This insignificant odds ratio means that although there are fewer submissions for testing in those counties with high coverage of evergreen forest as compared to other counties there is no difference in submission rate among the counties, and surveillance bias has been mitigated.

Finally, our study model has implications for the enzootic areas along an oral rabies vaccination (ORV) zone. This vaccination zone lies along the entire western edge of the eastern raccoon rabies reservoir (see Figure 1). Oral vaccination baits are dropped along roadways and from planes along this front to prevent the westward spread of the raccoon rabies variant. Evidence here suggests support for possibly using few or no baits in a pure stand of upland evergreen forest, at least in the southern US. In a similar fashion to dropping baits along a river in a riparian area, baits could be dropped up to the edge of an upland evergreen forest, as those raccoons in a pure upland evergreen forest may primarily use adjacent habitats for their resource needs21. The baits saved through this procedure could be distributed elsewhere with no loss in control effort. Conversely, bait density could be increased in areas that have been fragmented due to low intensity development.

Comments on this article Comments (0)

Version 2
VERSION 2 PUBLISHED 27 Dec 2013
Comment
Author details Author details
Competing interests
Grant information
Copyright
Download
 
Export To
metrics
Views Downloads
F1000Research - -
PubMed Central
Data from PMC are received and updated monthly.
- -
Citations
CITE
how to cite this article
Duke JE, Blanton JD, Ivey M and Rupprecht C. Modeling enzootic raccoon rabies from land use patterns - Georgia (USA) 2006-2010 [version 1; peer review: 2 approved]. F1000Research 2013, 2:285 (https://doi.org/10.12688/f1000research.2-285.v1)
NOTE: If applicable, it is important to ensure the information in square brackets after the title is included in all citations of this article.
track
receive updates on this article
Track an article to receive email alerts on any updates to this article.

Open Peer Review

Current Reviewer Status: ?
Key to Reviewer Statuses VIEW
ApprovedThe paper is scientifically sound in its current form and only minor, if any, improvements are suggested
Approved with reservations A number of small changes, sometimes more significant revisions are required to address specific details and improve the papers academic merit.
Not approvedFundamental flaws in the paper seriously undermine the findings and conclusions
Version 1
VERSION 1
PUBLISHED 27 Dec 2013
Views
46
Cite
Reviewer Report 06 Mar 2014
Nicole Gottdenker, Department of Pathology, University of Georgia, Athens, GA, USA 
Approved
VIEWS 46
General comments:

This article provides an interesting analysis of the available data on raccoon rabies from a public health database in order to evaluate how land use impacts raccoon rabies transmission. In particular, results show that evergreen forests (likely managed pine ... Continue reading
CITE
CITE
HOW TO CITE THIS REPORT
Gottdenker N. Reviewer Report For: Modeling enzootic raccoon rabies from land use patterns - Georgia (USA) 2006-2010 [version 1; peer review: 2 approved]. F1000Research 2013, 2:285 (https://doi.org/10.5256/f1000research.3243.r3926)
NOTE: it is important to ensure the information in square brackets after the title is included in all citations of this article.
  • Author Response 10 Mar 2014
    John Duke, Institute of Public Health, Georgia State University, Atlanta, Georgia 30303, USA
    10 Mar 2014
    Author Response
    Thank you so much for your comments.

    While the different statistical methodologies were discussed among the authors and others, the final decision was that, in our opinion, backwards stepwise regression was ... Continue reading
COMMENTS ON THIS REPORT
  • Author Response 10 Mar 2014
    John Duke, Institute of Public Health, Georgia State University, Atlanta, Georgia 30303, USA
    10 Mar 2014
    Author Response
    Thank you so much for your comments.

    While the different statistical methodologies were discussed among the authors and others, the final decision was that, in our opinion, backwards stepwise regression was ... Continue reading
Views
34
Cite
Reviewer Report 03 Mar 2014
Colleen Jonsson, Department of Microbiology and Immunology, University of Louisville, Louisville, KY, USA 
Approved
VIEWS 34
I find the content and data of this article acceptable. I would add a few more sentences ... Continue reading
CITE
CITE
HOW TO CITE THIS REPORT
Jonsson C. Reviewer Report For: Modeling enzootic raccoon rabies from land use patterns - Georgia (USA) 2006-2010 [version 1; peer review: 2 approved]. F1000Research 2013, 2:285 (https://doi.org/10.5256/f1000research.3243.r3369)
NOTE: it is important to ensure the information in square brackets after the title is included in all citations of this article.
  • Author Response 10 Mar 2014
    John Duke, Institute of Public Health, Georgia State University, Atlanta, Georgia 30303, USA
    10 Mar 2014
    Author Response
    Thank you so much for your comments.

    We felt that a reader could quickly look up the term "crude odds ratio" to have enough understanding to proceed.  I hope we made ... Continue reading
COMMENTS ON THIS REPORT
  • Author Response 10 Mar 2014
    John Duke, Institute of Public Health, Georgia State University, Atlanta, Georgia 30303, USA
    10 Mar 2014
    Author Response
    Thank you so much for your comments.

    We felt that a reader could quickly look up the term "crude odds ratio" to have enough understanding to proceed.  I hope we made ... Continue reading

Comments on this article Comments (0)

Version 2
VERSION 2 PUBLISHED 27 Dec 2013
Comment
Alongside their report, reviewers assign a status to the article:
Approved - the paper is scientifically sound in its current form and only minor, if any, improvements are suggested
Approved with reservations - A number of small changes, sometimes more significant revisions are required to address specific details and improve the papers academic merit.
Not approved - fundamental flaws in the paper seriously undermine the findings and conclusions
Sign In
If you've forgotten your password, please enter your email address below and we'll send you instructions on how to reset your password.

The email address should be the one you originally registered with F1000.

Email address not valid, please try again

You registered with F1000 via Google, so we cannot reset your password.

To sign in, please click here.

If you still need help with your Google account password, please click here.

You registered with F1000 via Facebook, so we cannot reset your password.

To sign in, please click here.

If you still need help with your Facebook account password, please click here.

Code not correct, please try again
Email us for further assistance.
Server error, please try again.