Previous Article |
Table of Contents
| Next Article
Ecology
Satellite imagery characterizes local animal reservoir populations of Sin Nombre virus in the southwestern United States








*W. Harry Feinstone Department of Molecular Microbiology and
Immunology and ||Department of Environmental Health
Sciences, The Johns Hopkins Bloomberg School of Public Health,
Baltimore, MD 21205;
Department of Biology and Museum of
Southwestern Biology, University of New
Mexico, Albuquerque, NM 87131;
Special Pathogens Branch, Division of Viral and
Rickettsial Diseases, Centers for Disease Control and Prevention,
Atlanta, GA 30333; and ¶IBM T.J. Watson Research Center,
Yorktown Heights, NY 10598
Communicated by M. Gordon Wolman, The Johns Hopkins University, Baltimore, MD and approved October 11, 2002 (received for review February 7, 2002)
| Abstract |
|---|
|
|
|---|
Abbreviations: ENSO, El Niño Southern Oscillation; TM, Thematic Mapper; HPS, hantaviral pulmonary syndrome; SNV, Sin Nombre virus
Recently, Glass et al. (5) used Landsat Thematic Mapper (TM) satellite data to identify environments associated with human risk of hantaviral disease caused by Sin Nombre virus (SNV). This retrospective study of human disease in the southwestern United States introduced an epidemiologic risk algorithm that combined several TM bands, taken nearly a year before the outbreak, and elevation, to estimate human risk in 1993, when disease was reported from nearly 30 locations. When the study was repeated using TM imagery from 1995 to predict risk in 1996 the extent of high-risk areas was substantially reduced and only a single human case was observed (5).
The hantaviral pulmonary syndrome (HPS) outbreak led to a series of studies of hantavirusrodent reservoir systems examining SNV and its main reservoir, Peromyscus maniculatus, the deer mouse (3, 68). Three key observations emerged from this research. First, as with other hantaviruses, SNV was maintained within a single reservoir species. Although spillover to other mammalian species occurred, viral persistence in the ecosystem relied on the deer mouse (6). Second, SNV appeared to be transmitted horizontally, so that the prevalence of infection increased with age class (or its surrogate, body mass) (3, 68). Thus, the prevalence of infection in deer mouse populations tended to peak with an aging population rather than with the absolute size of the population (4, 8, 9). Consequently, HPS risk could be higher in areas where populations showed a preponderance of adult animals (4, 9). Third, there was a sex bias in infection; more adult male than adult female mice were infected (3, 6). Mills et al. (3) surveyed >50 locations in the southwestern United States and found a crude prevalence of infection in female deer mice of 6%, compared with 14% in males. Among adult deer mice the difference was more striking with <10% of adult females with antibody to SNV, whereas nearly 30% of adult males were similarly positive.
The underlying dynamics that generate these patterns in local populations remain a matter of investigation (4, 811). However, Parmenter and colleagues (12, 13) hypothesized that, for several zoonotic agents in the region, local environmental conditions were major drivers that affected rodent populations and led to the increased risk of human disease. These events were triggered by meteorological conditions that affected vegetation and invertebrate populations that, in turn, affected the rodent reservoir populations. The changes in the rodent populations altered SNV transmission within the rodent populations (12) and, subsequently, human risk at a time remote from the triggering events themselves. Parmenter and colleagues proposed the trophic cascade hypothesis and they hypothesized that El Niño events could alter precipitation patterns in the southwestern United States sufficiently to trigger the trophic cascade.
In the spring of 1997, sea surface temperature patterns in the equatorial Pacific indicated a new El Niño Southern Oscillation (ENSO) event. This ENSO continued through mid-1998 (NOAA Climate Prediction Center, www.cpc.noaa.gov/products/analysis_monitoring/ensostuff/ensoyears.html) and provided an opportunity to examine the trophic cascade hypothesis and to evaluate the relationship between the epidemiologic risk algorithm by using satellite data and the ecology of the hantavirus system. We developed a series of predictions relating human risk estimated from satellite imagery with key aspects of the SNVdeer mouse system. High-risk sites were expected to have a faunal structure that was predominated by members of the genus Peromyscus, especially P. maniculatus, and that these local populations would differ in their age and sex structure from lower-risk sites. A consequence would be that levels of SNV infection within local populations in low- and high-risk sites should differ so that human populations in high-risk areas would more likely acquire infection. We further predicted that if the ENSO was an environmental trigger driving the trophic cascade leading to HPS outbreaks, then the dynamics of the SNVdeer mouse system during the ENSO should differ between high- and low-risk areas.
We tested the hypotheses by classifying human risk by using the risk algorithm (5) with satellite imagery from mid-1997 and mid-1998. A series of locations was then selected and sampled for small mammals by using a standardized protocol in 1998 and 1999. We show that high-risk sites were ecologically distinct from low-risk sites and that, by the end of the study, SNV infection was especially prevalent among deer mice in high-risk areas.
| Materials and Methods |
|---|
|
|
|---|
The images were processed and annual predicted HPS risk maps were generated by using described methods (5). Logistic regression was used to estimate risk for each year by using the digital numbers in each of three TM bands (bands 1, 5, and 7) and elevations at 28 HPS cases and 170 randomly selected control sites from the original 1993 outbreak. Risk at each pixel within the study region was assigned by weighting the digital numbers for each band, and elevation by the respective coefficients from the logistic models.
Pixels were categorized into three risk categories. Areas of low,
moderate, or high human risk were aggregated by using a receiver
operator characteristic (ROC) function (5). Low-risk areas were those
with risk values accounting for
10% of cases, moderate-risk areas
had values with >10% but <25% of cases, and high-risk areas had
signatures associated with
25% of cases. These categories were used
to sample rodent-trapping sites with similar levels of human risk.
Rodent Faunal Sampling. In 1998, 40 sites encompassing a wide range of predicted HPS risk (0.050.86 from the 1997 images) within the study area were selected. Sites were chosen by an individual unfamiliar with the region to reduce selection bias. The criteria used in site selection were that the sites were within relatively homogeneous patches of risk categories and that they were within 300 m of a road to allow access. For the field surveys, sites were dichotomized into those above and below the high-risk threshold, and two sites above and two sites below predicted high risk were surveyed at a time. This methodology was devised to control for environmental conditions (e.g., meteorological and lunar conditions) that could affect trap success. The field crews were blinded to the predicted human risk level to further reduce sampling bias. These sites were resampled in the spring of 1999 by using the same protocols. The 1998 TM images used to generate the risk map for 1999 did not include two of the original survey sites in moderate-risk areas, and these were excluded from further analyses.
A standardized trapping protocol was used to sample the local small
mammal fauna (3, 14). One hundred single-capture small mammal live
traps were set for two nights as four parallel lines of 25 traps each.
Traps were spaced at
10-m intervals. Longitude and latitude were
recorded at the contralateral corners of the two outer trap lines by
using global positioning system (GPS) receivers. The data were used to
locate the sampling sites on the satellite images. The major vegetation
in the overstory, understory, and herbaceous layers and the extent of
cover, as well as the amount of cover by litter and bare soil or rock,
were recorded by using a standardized protocol (14). These data were
used to qualitatively compare sites for similarities in plant community
structure.
Standard measures of faunal composition and abundance were generated from the trapping data. Trap success for selected species was measured as the number of animals captured divided by the number of trap nights (maximally 200) at each study site and used as an indication of abundance. Corrections were made for traps that were lost or sprung but empty (14). The age, sex, body mass, and species identification were recorded for each animal captured. These data were used to characterize the age and sex structure of populations for selected species (3). Collected animals were processed according to protocols for SNV detection (15). All animals were accessioned to the Museum of Southwestern Biology (MSB), University of New Mexico, and retained as voucher specimens [NK7796978300 (1998) and NK8491985500 (1999)]. Frozen tissues and blood samples from voucher specimens were deposited in the Division of Genomic Resources, MSB.
Serological Testing. Sera from captured rodents were tested by ELISA for antibodies to hantaviruses by using described methods (16). Briefly, recombinant SNV nucleocapsid antigen cloned into a vector was expressed and coated onto ELISA plates. As control antigen, protein from the vector alone was expressed and coated onto plates. Sera were screened at a dilution of 1:100 in PBSTween-20 (0.1%) and 5% skim milk. Sera were incubated in duplicate wells and detected by using a combination of goat anti-Peromyscus IgG and goat anti-Rattus IgG conjugated with horseradish peroxidase followed with 2,2'-azinodi(3-ethyl)benzthiazoline-6-sulfonic acid (ABTS). Optical densities were read at 410 nm and the difference between the pair of SNV antigen wells and the pair of control wells for each sample was calculated.
Statistical Analysis. The human risk for HPS at each trapping site was determined by averaging the risk for all of the pixels incorporated within the area enclosed by the trapping lines. Prevalence of SNV antibody at the small mammal sampling locations was estimated by dividing the number of individuals that were antibody positive by the numbers of individuals tested. Previously identified risk factors for hantavirus infection in small mammals (species, sex, and age composition) were used to compare risk sites (3).
We tested the hypothesis that the predicted HPS risk derived from the satellite imagery identified biologically relevant features associated with hantavirus transmission by predicting characteristics of local rodent populations. First, we examined the relationship between P. maniculatus abundance at the 38 trapping sites and the HPS risk by using multiple linear regression. The number of deer mice captured (logarithmically transformed captures + 1) was the dependent variable. The HPS risk obtained from the previous June's satellite image and the change in the risk from the prior year, for 1999, were used as the independent variables. Satellite images from 1996 were not available to allow us to calculate changes in risk from 1997 to 1998, so the analysis for 1998 was performed as a simple linear regression. We predicted that there should be a significant relationship between the abundance of deer mice captured at a site and the predicted risk of HPS because deer mice are the principal reservoir of SNV.
As a second test, the age and sex structures of deer mouse populations at high-risk sites were compared with low-risk sites. Adult male deer mice are disproportionately infected with SNV compared with females and juveniles of either sex. Thus, we hypothesized that high-risk sites would have a population sampled that was biased toward adult male deer mice.
Finally, we examined the prevalence of SNV infection, as measured by the presence of specific antibodies in individual animals, in high- and low-risk sites. We predicted that the prevalence of SNV infection should be higher in high-risk sites, especially in the year after the ENSO, as this was the time lag associated with the previous HPS outbreak in 1993.
| Results |
|---|
|
|
|---|
|
Trapping Results. A total of 15,042 trap nights were completed during sampling at the 38 trapping sites (7,527 trap nights and 7,515 trap nights in 1998 and 1999, respectively). Thirteen species of small mammals were captured, with the preponderance (77.8%; n = 826) being members of the genus Peromyscus. Small mammal captures were more common in 1999 (494 animals = 6.4% trap success), the year after the ENSO, than during 1998 (332 animals = 4.41% trap success; P < 0.0001). P. maniculatus was the most common species captured and was found in both low- to moderate- and high-risk areas (Fig. 2). Peromyscus boylii and Peromyscus truei were more commonly observed in high-risk areas than in low- to moderate-risk areas.
|
![]() |
![]() |
Overall, the abundance of deer mice varied between years and was
related to the predicted human risk (Fig. 2). At least one deer mouse
was captured at 10/11 lowmoderate-risk sites in 1998
(
= 8.0 mice per site; SD = 16.06; range,
056 mice). However, the average abundance was skewed by results from
a single location. When that one moderate-risk site was excluded,
capture statistics were substantially lower (
=
3.2; SD = 2.20; range, 06 mice). Deer mice were captured at only
4/7 sites that changed from high risk in 1998 to lowmoderate risk
in the following year (
= 2.0 mice per site;
SD = 2.08; range, 05 mice). Deer mice were trapped at 17/20
sites in persistently high-risk areas in 1998 (
= 5.4 mice per site; SD = 5.19; range, 017 mice). The
differences in average local deer mouse abundance by risk category in
1998 were not statistically significant (KruskalWallis;
H = 1.37; df = 2,35; P = 0.27).
In 1999, the average abundance of deer mice in lowmoderate risk sites
was not significantly different from that observed in 1998
(
= 2.1; SD = 4.50; range, 015). However,
deer mice were captured at fewer lowmoderate risk sites (4/11 sites
vs. 10/11 sites). Sites that had changed from high risk to low risk
in 1999 had marginally lower numbers of deer mice
(
= 0.9; SD = 0.90; range, 02) and deer
mice were captured at only four of seven sites. In contrast, all
high-risk sites had at least two deer mice and the numbers of deer mice
increased in high-risk areas (
= 9.4; SD =
0.5.79; range, 224). The average local abundance of deer mice at
high-risk sites was significantly greater than at other sites
(KruskalWallis; H = 24.42; df = 2,35;
P < 0.001).
In addition to the increase in the abundance of reservoir populations
in high-risk areas in 1999, there were differences in the deer mouse
population structure depending on the level of risk. The sex ratios
varied depending on the year and the risk categories of sites (Fig.
3). In both years the lowmoderate-risk
sites and sites where risk varied between years had sex ratios that did
not differ from 1:1. Populations in high-risk sites, however, were
biased toward males (2.97:1 in 1998 and 1.24:1 in 1999). In 1998, this
difference was highly significant (
2 = 15.11;
df = 1; P = 0.0001), and, while the trend
persisted in 1999, it was not statistically significant
(
2 = 1.07; df = 1; P =
0.30). Populations in high-risk areas also had a biased age structure
(Fig. 4). Populations from high-risk
areas during both 1998 and 1999 had a disproportionate fraction of
adult animals compared with other risk areas
(
2 = 12.74; df = 5; P =
0.026). By 1999, subadults represented the largest class of individuals
in populations from lowmoderate-risk sites, whereas adults continued
to represent the largest age class in high-risk areas (Fig. 4).
|
|
2 = 0.32; df = 1; P =
0.570). The crude prevalence of infection (6.6%; n =
181) did not differ significantly between high- and low-risk areas.
Infection was predominantly in male (91.7%; n = 12),
adult (75%; n = 12) mice. Three subadult P.
maniculatus also were SNV-antibody-positive. Only one animal of
another species (1/26; P. truei) had antibodies to SNV in
1998.
With the end of the El Niño in 1999 more high-risk sites had
infected mice and the prevalence of infection in deer mice from the
high-risk areas was markedly increased (Fig.
5). Overall, SNV-antibody-positive deer
mice were captured at more sites than in 1998
(
2 = 11.94; df = 1; P =
0.0006). In 1999, 19/38 (50.0%) of the sites had deer mice that were
antibody positive. However, infected deer mice were found nearly
exclusively in high-risk areas. Only 1/11 lowmoderate-risk sites
had antibody-positive animals, whereas infected deer mice were found at
18/20 high-risk sites (
2 = 20.74; df =
1; P < 0.0001). The crude prevalence of infection
among mice differed between lowmoderate- and high-risk sites (Fig.
5). The prevalence (8.3%; n = 24) in
lowmoderate-risk sites did not differ from the prevalence in those
sites in 1998, whereas the prevalence of SNV infection among deer mice
in high-risk sites (30.8%; n = 185) increased from
1998 (
2 = 32.97; df = 1;
P < 0.0001). During 1999, infected deer mice were
predominantly male (76.3%; n = 59), adult
(67.8%; n = 59) mice, although 1 juvenile female and
18 subadults (3 females and 15 males) also were antibody positive. The
prevalence of SNV in adult deer mice within high-risk sites was 36.4%
(40/110). In addition to antibody-positive P. maniculatus,
5/36 (13.9%) P. boylii and 10/68 (14.7%) and P.
truei also were positive to SNV during 1999. All positive
P. boylii and P. truei were captured in high-risk
areas.
|
| Discussion |
|---|
|
|
|---|
Predicted high-risk areas for humans were characterized by a small mammal fauna dominated by the known reservoir species, P. maniculatus, and local rodent populations in high-risk areas were dominated by those subclasses of the population that were most likely to carry SNV. Deer mice occurred in persistently low-risk areas as well as sites that showed significant year-to-year fluctuations in risk. However, populations in these areas were subject to local extinction after the ENSO, as evidenced by the increasing number of these sites that failed to have any deer mice in 1999 (10/18 in 1999 vs. 4/18 in 1998) and decreases in the numbers of animals caught in the sites that did have mice (Fig. 2). The characteristics of local reservoir populations in high-risk areas differed from those of low-risk areas so much that, after the 1998 ENSO the chance of exposure to SNV was more likely in high-risk areas (Figs. 2 and 5). Populations of deer mice in high-risk areas continued to increase in abundance during 1999, and mice were captured in all high-risk sites. Most significantly, by 1999, the prevalence of infection with SNV exceeded 36% in adult deer mice from high-risk areas, while infection occurred only sporadically (1/18 sites) and at low prevalence (8.3%) in low-risk areas.
The precise ecological factors that link the increased levels of SNV in some local populations of deer mice, the rodent population dynamics, and the signature associated with the satellite imagery remain to be determined. We postulate that many local populations become established when animals colonize habitats as conditions ameliorate during an ENSO. The colonists are derived from animals that disperse from locations with persistently suitable environmental conditions, but colonizing populations are subject to periodic extinctions as local conditions become more stringent after an ENSO. The likelihood of extinction in the colonized areas depends on annual changes in environmental conditions that may be related to soil moisture, as indicated by rapid changes in predicted risk at some sites between 1998 and 1999 (Fig. 1). Low-risk sites represent environmental conditions where populations are able to survive only for short periods of time (<2 years), whereas environmental conditions in high-risk sites favor the establishment and persistence of deer mouse populations for longer periods (9, 17). Because SNV is horizontally transmitted, populations must persist for a sufficient time to allow for the transmission among local population members to reach a level where human risk is elevated. Empirically, this appears to require at least 2 years (Fig. 5). Rodent populations in sites where extinctions are common do not survive long enough (or reach sufficient population densities) for SNV to become widespread. The behavioral ecology of deer mice, with male-dominant dispersal and territoriality (18), may account for the observed demographic characteristics of high-risk areas (Figs. 35). Sites that favor population persistence and SNV transmission are predominantly colonized by male mice. These mice defend their territories from other males through aggressive interactions (18). Thus, when suitable sites for SNV persistence are initially colonized, they have sex- and age-biased populations, but later they evolve toward a more equitable sex ratio because of in situ reproduction (as seen in 1999). However, behavioral interactions would cause most juveniles and subadults to disperse from these more suitable habitats, resulting in persistent age-biased populations. In addition, aggressive interactions associated with SNV transmission would cause dispersing subadult and other adult male population members to spread SNV, leading to increased human HPS risk, and the dispersal of SNV to other areas. In contrast, low-risk areas, being less suitable over the long run, may have similar population characteristics to higher-risk sites during colonization but serve as population sinks, as environmental conditions became more severe in subsequent years. These sites have lower population densities and a higher likelihood of extinction, and would be primarily occupied by juvenile and subadult animals (Fig. 4) with associated lower levels of SNV infection.
This hypothesis is consistent with several surveys of deer mice that indicate either the presence of infected mice, or the occurrence of human disease, was associated with specific vegetational zones, and, presumably, environmental conditions (3, 7, 19, 20). Studies of the interplay between local habitat characteristics and the dynamics of reservoir populations show they affect SNV persistence and transmission in reservoir populations (9, 17) much as our analysis proposes. A primary difference is that we propose that vegetation community is a surrogate (albeit a highly correlated one) for the underlying environmental risk factors affecting mouse populations. Our qualitative characterization of habitat indicates that predicted high-risk areas are generally those associated with higher elevation and more mesic vegetation, whereas low-risk areas tended to be vegetation communities associated with lower elevation and xeric plant communities, as the earlier studies indicated. However, the subset of high-risk sites in 1998 that were reclassified as low-risk sites in 1999 had vegetation structure consistent with lower-elevation plant communities (Fig. 1). This finding indicates that the algorithm using the TM imagery recognized something in the environment, rather than vegetation as associated with human risk. Glass et al. (5) suggested that soil moisture may be a factor recognized by the classification algorithm which would account for the loss of some predicted high-risk areas after the end of El Niño in 1998.
Much of the high-risk area identified from the satellite imagery lay at elevations that are only rarely occupied by people. Areas that have persistently remained at high risk since HPS was recognized in 1993 correspond to a relatively small fraction of the environment (data not shown) associated with arroyos and valleys leading from higher-elevation areas. These areas may serve as refugia from which deer mice disperse to lower elevations as environmental conditions become more suitable.
The high prevalence of SNV infection in rodent populations in high-risk areas after an ENSO is consistent with observations from sites where human cases have occurred (6). Prevalences of 30% in those survey samples were not uncommon. Practically, this result suggests the reported prevalences in rodents during case investigations do not represent selection bias. This study also shows that a substantial elevation in prevalence occurs over a year's time in high-risk areas. The lag time between the end of the 19971998 ENSO and the increase in SNV prevalence in high-risk areas mirrored the 1993 outbreak of HPS, which occurred within a year of the 19911992 ENSO. This time lag also is consistent with the trophic cascade hypothesis (12, 13). These results also demonstrate that increased SNV prevalence occurs locally, while prevalence remains low or absent in nearby locations, emphasizing that local conditions significantly impact virus transmission. Equally important, areas of high SNV prevalence can be detected by satellite.
Future studies, using satellite imagery to identify areas that remain high risk, could characterize the ecological dynamics of local environmental conditions and monitor SNV transmission in deer mouse populations and compare these sites with lower-risk locations. We believe this approach will help identify the environmental determinants of SNV persistence in the environment. These sites may provide important insights into identifying the environmental conditions that lead to increased levels of SNV in reservoir populations and the subsequent increased risk of human disease.
| Acknowledgements |
|---|
| Footnotes |
|---|
To whom correspondence should be addressed at: The
Johns Hopkins Bloomberg School of Public Health, 615 North Wolfe
Street, Baltimore, MD 21205. E-mail: ggurrigl{at}jhsph.edu. | References |
|---|
|
|
|---|
This article has been cited by other articles in HighWire Press-hosted journals:
![]() |
J. R. Levine, C. L. Fritz, and M. G. Novak Occupational Risk of Exposure to Rodent-borne Hantavirus at US Forest Service Facilities in California Am J Trop Med Hyg, February 1, 2008; 78(2): 352 - 357. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. S. CARROLL, J. N. MILLS, J. M. MONTGOMERY, D. G. BAUSCH, P. J. BLAIR, J. P. BURANS, V. FELICES, A. GIANELLA, N. IIHOSHI, S. T. NICHOL, et al. HANTAVIRUS PULMONARY SYNDROME IN CENTRAL BOLIVIA: RELATIONSHIPS BETWEEN RESERVOIR HOSTS, HABITATS, AND VIRAL GENOTYPES Am J Trop Med Hyg, January 1, 2005; 72(1): 42 - 46. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||