Landscape-Genetic Analysis of Population Structure in the Texas Gray Fox Oral Rabies Vaccination Zone

In west-central Texas, USA, abatement efforts for the gray fox (Urocyon cinereoargenteus) rabies epizootic illustrate the difficulties inherent in large-scale management of wildlife disease. The rabies epizootic has been managed through a cooperative oral rabies vaccination program (ORV) since 1996. Millions of edible baits containing a rabies vaccine have been distributed annually in a 16-km to 24-km zone around the perimeter of the epizootic, which encompasses a geographic area >4 x 105 km2. The ORV program successfully halted expansion of the epizootic into metropolitan areas but has not achieved the ultimate goal of eradication. Rabies activity in gray fox continues to occur periodically outside the ORV zone, preventing ORV zone contraction and dissipation of the epizootic. We employed a landscape-genetic approach to assess gray fox population structure and dispersal in the affected area, with the aim of assisting rabies management efforts. No unique genetic clusters or population boundaries were detected. Instead, foxes were weakly structured over the entire region in an isolation by distance pattern. Local subpopulations appeared to be genetically non-independent over distances >30 km, implying that long-distance movements or dispersal may have been common in the region. We concluded that gray foxes in west-central Texas have a high potential for long-distance rabies virus trafficking. Thus, a 16-km to 24-km ORV zone may be too narrow to contain the fox rabies epizootic. Continued expansion of the ORV zone, although costly, may be critical to the long-term goal of eliminating the Texas fox rabies virus variant from the United States.

Contributions of demography and dispersal parameters to the spatial spread of a stage-structured insect invasion

Stage-structured models that integrate demography and dispersal can be used to identify points in the life cycle with large effects on rates of population spatial spread, information that is vital in the development of containment strategies for invasive species. Current challenges in the application of these tools include: (1) accounting for large uncertainty in model parameters, which may violate assumptions of ‘‘local’’ perturbation metrics such as sensitivities and elasticities, and (2) forecasting not only asymptotic rates of spatial spread, as is usually done, but also transient spatial dynamics in the early stages of invasion. We developed an invasion model for the Diaprepes root weevil (DRW; Diaprepes abbreviatus [Coleoptera: Curculionidae]), a generalist herbivore that has invaded citrus-growing regions of the United States. We synthesized data on DRW demography and dispersal and generated predictions for asymptotic and transient peak invasion speeds, accounting for parameter uncertainty. We quantified the contributions of each parameter toward invasion speed using a ‘‘global’’ perturbation analysis, and we contrasted parameter contributions during the transient and asymptotic phases. We found that the asymptotic invasion speed was 0.02–0.028 km/week, although the transient peak invasion speed (0.03– 0.045 km/week) was significantly greater. Both asymptotic and transient invasions speeds were most responsive to weevil dispersal distances. However, demographic parameters that had large effects on asymptotic speed (e.g., survival of early-instar larvae) had little effect on transient speed. Comparison of the global analysis with lower-level elasticities indicated that local perturbation analysis would have generated unreliable predictions for the responsiveness of invasion speed to underlying parameters. Observed range expansion in southern Florida (1992–2006) was significantly lower than the invasion speed predicted by the model. Possible causes of this mismatch include overestimation of dispersal distances, demographic rates, and spatiotemporal variation in parameter values. This study demonstrates that, when parameter uncertainty is large, as is often the case, global perturbation analyses are needed to identify which points in the life cycle should be targets of management. Our results also suggest that effective strategies for reducing spread during the asymptotic phase may have little effect during the transient phase. Includes Appendix.