Decelerating spread of West Nile virus by percolation in a heterogeneous urban landscape

Vector-borne diseases are emerging and re-emerging in urban environments throughout the world, presenting an increasing challenge to human health and a major obstacle to development. Currently, more than half of the global population is concentrated in urban environments, which are highly heterogeneous in the extent, degree, and distribution of environmental modifications. Because the prevalence of vector-borne pathogens is so closely coupled to the ecologies of vector and host species, this heterogeneity has the potential to significantly alter the dynamical systems through which pathogens propagate, and also thereby affect the epidemiological patterns of disease at multiple spatial scales. One such pattern is the speed of spread. Whereas standard models hold that pathogens spread as waves with constant or increasing speed, we hypothesized that heterogeneity in urban environments would cause decelerating travelling waves in incipient epidemics. To test this hypothesis, we analysed data on the spread of West Nile virus (WNV) in New York City (NYC), the 1999 epicentre of the North American pandemic, during annual epizootics from 2000–2008. These data show evidence of deceleration in all years studied, consistent with our hypothesis. To further explain these patterns, we developed a spatial model for vector-borne disease transmission in a heterogeneous environment. An emergent property of this model is that deceleration occurs only in the vicinity of a critical point. Geostatistical analysis suggests that NYC may be on the edge of this criticality. Together, these analyses provide the first evidence for the endogenous generation of decelerating travelling waves in an emerging infectious disease. Since the reported deceleration results from the heterogeneity of the environment through which the pathogen percolates, our findings suggest that targeting control at key sites could efficiently prevent pathogen spread to remote susceptible areas or even halt epidemics.

Current theory of the spatial spread of pathogens predicts travelling waves at constant or increasing speed in homogeneous environments. However, in urban environments, increasing and often unregulated development produces a highly heterogeneous landscape. Such heterogeneity affects pathogens spread by insect vectors particularly, which typically have short dispersal distances. We hypothesized that high levels of heterogeneity can slow the spread of such pathogens, resulting in decelerating epidemic waves. We analysed the annual spread of West Nile virus (WNV) in New York City (NYC), using a dataset containing >1,000,000 records since the origin of the North American pandemic in 1999. Our analysis provides the first evidence of endogenous decelerating travelling waves in an emerging infectious disease. We found that WNV spread with decreasing speed in each season and rejected four alternative hypotheses to explain this deceleration. A mathematical model shows that high levels of heterogeneity can lead to such decelerating travelling waves. Interestingly, the level of heterogeneity in land-cover types associated with WNV-positive dead birds in NYC is of the order of magnitude required to produce decelerating travelling waves in the model. Consequently, we propose that control strategies targeting key sites may be effective at slowing WNV spread in NYC.

Utility of frameless stereotaxy in the resection of skull base and Basal cerebral lesions in children

Since 1991, we have performed nearly 300 stereotactic procedures using the ISG viewing wand on a variety of cranial lesions in patients under 22 years of age. Of these, 38 procedures were performed on 34 patients for basal cerebral and skull base lesions. Our patients ranged in age from 3.5 months to 22 years with a mean age of 9.45 years. There were 18 females and 16 males. Twenty-one patients had basal cerebral lesions located in the thalamus (10), basal ganglia (2), third ventricle (2), and hypothalamus (7). Thirteen patients had skull base lesions located within the anterior optic apparatus (3), sella turcica (4), middle and posterior cranial fossae (4), and craniocervical region (2). Preoperative CT and/or MRI scan images were taken as a volume acquisition and transferred to the computer workstation utilizing the ISG Wand software. This workstation was transferred to the operating room where it was calibrated to a faro Surgicom arm which interfaces with the patient and the three-dimensional radiological image. The ISG Wand was utilized to plan the scalp and bone flaps and to select the optional trajectory to lesion. The surgical approaches which were specifically used in this series with the ISG Wand included transcallosal (15), pterional (5), frontal (3), subtemporal (4), transsphenoidal (3), temporal (3), tumor cyst shunt insertion (1), burr hole drainage (1), transoral (2), bifrontal (1), bifrontal mid facial (1), and transnasal (1). Although brain shift occurred following craniotomy and with brain retraction, the relative immobility of these lesions at the skull or cerebral base permitted an accurate targeting of all lesions with an error range of 1.0-2.5 mm throughout the entire procedure. This relatively precise intraoperative feedback led to more accurate recognition of tumor landmarks. It is the authors' impression that a more aggressive resection of these lesions was achieved than could be without the device. We conclude that a frameless stereotactic device such as the ISG Wand is particularly valuable in the approach to skull base and basal cerebral tumors in children.

The curse of the Pharaoh revisited: evolutionary bi-stability in environmentally transmitted pathogens

It is increasingly evident that for a number of high-profile pathogens, transmission involves both direct and environmental pathways. Much of the distinguished evolutionary theory has, however, focused on each of transmission component separately. Herein, we use the framework of adaptive dynamics to study the evolutionary consequences of mixed transmission. We find that environmental transmission can select for increased virulence when direct transmission is low. Increasing the efficiency of direct transmission gives rise to an evolutionary bi-stability, with coexistence of different levels of virulence. We conclude that the overlooked contribution of environmental transmission may explain the curious appearance of high virulence in pathogens that are typically only moderately pathogenic, as observed for avian influenza viruses and cholera.

Why are some plant genera more invasive than others?

Determining how biological traits are related to the ability of groups of organisms to become economically damaging when established outside of their native ranges is a major goal of population biology, and important in the management of invasive species. Little is known about why some taxonomic groups are more likely to become pests than others among plants. We investigated traits that discriminate vascular plant genera, a level of taxonomic generality at which risk assessment and screening could be more effectively performed, according to the proportion of naturalized species which are pests. We focused on the United States and Canada, and, because our purpose is ultimately regulatory, considered species classified as weeds or noxious. Using contingency tables, we identified 11 genera of vascular plants that are disproportionately represented by invasive species. Results from boosted regression tree analyses show that these categories reflect biological differences. In summary, approximately 25% of variation in genus proportions of weeds or noxious species was explained by biological covariates. Key explanatory traits included genus means for wetland habitat affinity, chromosome number, and seed mass.

Experimental demonstration of a two-phase population extinction hazard

Population extinction is a fundamental biological process with applications to ecology, epidemiology, immunology, conservation biology and genetics. Although a monotonic relationship between initial population size and mean extinction time is predicted by virtually all theoretical models, attempts at empirical demonstration have been equivocal. We suggest that this anomaly is best explained with reference to the transient properties of ensembles of populations. Specifically, we submit that under experimental conditions, many populations escape their initially vulnerable state to reach quasi-stationarity, where effects of initial conditions are erased. Thus, extinction of populations initialized far from quasi-stationarity may be exposed to a two-phase extinction hazard. An empirical prediction of this theory is that the fit Cox proportional hazards regression model for the observed survival time distribution of a group of populations will be shown to violate the proportional hazards assumption early in the experiment, but not at later times. We report results of two experiments with the cladoceran zooplankton Daphnia magna designed to exhibit this phenomenon. In one experiment, habitat size was also varied. Statistical analysis showed that in one of these experiments a transformation occurred so that very early in the experiment there existed a transient phase during which the extinction hazard was primarily owing to the initial population size, and that this was gradually replaced by a more stable quasi-stationary phase. In the second experiment, only habitat size unambiguously displayed an effect. Analysis of data pooled from both experiments suggests that the overall extinction time distribution in this system results from the mixture of extinctions during the initial rapid phase, during which the effects of initial population size can be considerable, and a longer quasi-stationary phase, during which only habitat size has an effect. These are the first results, to our knowledge, of a two-phase population extinction process.

Regional differences in the association between land cover and West Nile virus disease incidence in humans in the United States

West Nile virus (WNV) is generally considered to be an urban pathogen in the United States, but studies associating land cover and disease incidence, seroprevalence, or infection rate in humans, birds, domesticated and wild mammals, and mosquitoes report varying and sometimes contradictory results at an array of spatial extents. Human infection can provide insight about basic transmission activity; therefore, we analyzed data on the incidence of WNV disease in humans to obtain a comprehensive picture of how human disease and land cover type are associated across the United States. Human WNV disease incidence in Northeastern regions was positively associated with urban land covers, whereas incidence in the Western United States was positively associated with agricultural land covers. We suggest that these regional associations are explained by the geographic distributions of prominent WNV vectors: Culex pipiens complex (including Cx. pipiens and Cx. quinquefasciatus) in the Northeast and Cx. tarsalis in the Western United States.

Time since introduction, seed mass, and genome size predict successful invaders among the cultivated vascular plants of Hawaii

Extensive economic and environmental damage has been caused by invasive exotic plant species in many ecosystems worldwide. Many comparative studies have therefore attempted to predict, from biological traits, which species among the pool of naturalized non-natives become invasive. However, few studies have investigated which species establish and/or become pests from the larger pool of introduced species and controlled for time since introduction. Here we present results from a study aimed at quantifying predicting three classes of invasive species cultivated in Hawaii. Of 7,866 ornamental species cultivated in Hawaii between 1840 and 1999, 420 (5.3%) species naturalized, 141 (1.8%) have been classified as weeds, and 39 (0.5%) were listed by the state of Hawaii as noxious. Of the 815 species introduced >80 years ago, 253 (31%) have naturalized, 90 (11%) are classed as weeds, and 22 (3%) as noxious by the state of Hawaii. Using boosted regression trees we classified each group with nearly 90% accuracy, despite incompleteness of data and the low proportion of naturalized or pest species. Key biological predictors were seed mass and highest chromosome number standardized by genus which, when data on residence time was removed, were able to predict all three groups with 76-82% accuracy. We conclude that, when focused on a single region, screening for potential weeds or noxious plants based on a small set of biological traits can be achieved with sufficient accuracy for policy and management purposes.

Critical patch size generated by Allee effect in gypsy moth, Lymantria dispar (L.)

Allee effects are important dynamical mechanisms in small-density populations in which per capita population growth rate increases with density. When positive density dependence is sufficiently severe (a ‘strong’ Allee effect), a critical density arises below which populations do not persist. For spatially distributed populations subject to dispersal, theory predicts that the occupied area also exhibits a critical threshold for population persistence, but this result has not been confirmed in nature. We tested this prediction in patterns of population persistence across the invasion front of the European gypsy moth (Lymantria dispar) in the United States in data collected between 1996 and 2008. Our analysis consistently provided evidence for effects of both population area and density on persistence, as predicted by the general theory, and confirmed here using a mechanistic model developed for the gypsy moth system. We believe this study to be the first empirical documentation of critical patch size induced by an Allee effect.

Development of a sarcoidosis murine lung granuloma model using Mycobacterium superoxide dismutase A peptide

Sarcoidosis is characterized by noncaseating granulomas containing CD4(+) T cells with a Th1 immunophenotype. Although the causative antigens remain unknown, independent studies noted molecular and immunologic evidence of mycobacterial virulence factors in sarcoidosis specimens. A major limiting factor in discovering new insights into the pathogenesis of sarcoidosis is the lack of an animal model. Using a distinct superoxide dismutase A peptide (sodA) associated with sarcoidosis granulomas, we developed a pulmonary model of sarcoidosis granulomatous inflammation. Mice were sensitized by a subcutaneous injection of sodA, incorporated in incomplete Freund's adjuvant (IFA). Control subjects consisted of mice with no sensitization (ConNS), sensitized with IFA only (ConIFA), or with Schistosoma mansoni eggs. Fourteen days later, sensitized mice were challenged by tail-vein injection of naked beads, covalently coupled to sodA peptides or to schistosome egg antigens (SEA). Histologic analysis revealed hilar lymphadenopathy and noncaseating granulomas in the lungs of sodA-treated or SEA-treated mice. Flow cytometry of bronchoalveolar lavage (BAL) demonstrated CD4(+) T-cell responses against sodA peptide in the sodA-sensitized mice only. Cytometric bead analysis revealed significant differences in IL-2 and IFN-γ secretion in the BAL fluid of sodA-treated mice, compared with mice that received SEA or naked beads (P = 0.008, Wilcoxon rank sum test). ConNS and ConIFA mice demonstrated no significant formation of granuloma, and no Th1 immunophenotype. The use of microbial peptides distinct for sarcoidosis reveals a histologic and immunologic profile in the murine model that correlates well with those profiles noted in human sarcoidosis, providing the framework to investigate the molecular basis for the progression or resolution of sarcoidosis.

A general multi-strain model with environmental transmission: invasion conditions for the disease-free and endemic states

Although many infectious diseases of humans and wildlife are transmitted via an environmental reservoir, the theory of environmental transmission remains poorly elaborated. Here we introduce an SIR-type multi-strain disease transmission model with perfect cross immunity where environmental transmission is broadly defined by three axioms. We establish the conditions under which a multi-strain endemic state is invaded by another strain which is both directly and environmentally transmitted. We discuss explicit forms for environmental transmission terms and apply our newly derived invasion conditions to a two-strain system. Then, we consider the case of two strains with matching basic reproduction numbers (i.e., R(0)), one directly transmitted only and the other both directly and environmentally transmitted, invading each other's endemic state. We find that the strain which is only directly transmitted can invade the endemic state of the strain with mixed transmission. However, the endemic state of the first strain is neutrally stable to invasion by the second strain. Thus, our results suggest that environmental transmission makes the endemic state less resistant to invasion.

Experimental demonstration of population extinction due to a predator-driven Allee effect

1. Allee effects may result in negative growth rates at low population density, with important implications for conservation and management of exploited populations. Theory predicts prey populations will exhibit Allee effects when their predator exhibits a Type II functional response, but empirical evidence linking this positively density-dependent variation in predator-induced individual mortality to population growth rate and probability of extinction is lacking. 2. Here, we report a demonstration of extinction due to predator-driven Allee effects in an experimental Daphnia-Chaoborus system. A component Allee effect caused by higher predation rates at low Daphnia density led to positive density dependence in per capita growth rate and accelerated extinction rate at low density. 3. A stochastic model of the process revealed how the critical density below which population growth is negative depends on the mechanistic details of the predator-prey interaction. 4. The ubiquity of predator-prey interactions and saturating functional responses suggests predator-driven Allee effects are potentially important in determining extinction risk of a large number of species.

Evolutionary relationships among human-isolated and wildlife-isolated West Nile viruses

The evolutionary relationships among pathogen lineages in multi-host systems are often the only observable signature of unobserved ecological and epidemiological processes. The evolution of viruses infecting humans, particularly, is of interest because of the public health importance of understanding the relationship of virus exposure to disease risk. Here I report results of two analyses of the evolutionary relationships among West Nile viruses in North America. These analyses suggest that (1) assortative mixing occurs between virus groups and human vs. non-human hosts and (2) human-derived isolates are related to each other. The ecological processes generating these viruses and the epidemiological consequences of West Nile virus host preference are unknown.

Environment, but not migration rate, influences extinction risk in experimental metapopulations

Ecological theory suggests that several demographic factors influence metapopulation extinction risk, including synchrony in population size between subpopulations, metapopulation size and the magnitude of fluctuations in population size. Theoretically, each of these is influenced by the rate of migration between subpopulations. Here we report on an experiment where we manipulated migration rate within metapopulations of the freshwater zooplankton Daphnia magna to examine how migration influenced each of these demographic variables, and subsequent effects on metapopulation extinction. In addition, our experimental procedures introduced unplanned but controlled differences between metapopulations in light intensity, enabling us to examine the relative influences of environmental and demographic factors. We found that increasing migration rate increased subpopulation synchrony. We failed to detect effects of migration on population size and fluctuations in population size at the metapopulation or subpopulation level, however. In contrast, light intensity did not influence synchrony, but was positively correlated with population size and negatively correlated with population fluctuation. Finally, synchrony did not influence time to extinction, while population size and the magnitude of fluctuations did. We conclude that environmental factors had a greater influence on extinction risk than demographic factors, and that metapopulation size and fluctuation were more important to extinction risk than metapopulation synchrony.

Environmental transmission of low pathogenicity avian influenza viruses and its implications for pathogen invasion

Understanding the transmission dynamics and persistence of avian influenza viruses (AIVs) in the wild is an important scientific and public health challenge because this system represents both a reservoir for recombination and a source of novel, potentially human-pathogenic strains. The current paradigm locates all important transmission events on the nearly direct fecal/oral bird-to-bird pathway. In this article, on the basis of overlooked evidence, we propose that an environmental virus reservoir gives rise to indirect transmission. This transmission mode could play an important epidemiological role. Using a stochastic model, we demonstrate how neglecting environmentally generated transmission chains could underestimate the explosiveness and duration of AIV epidemics. We show the important pathogen invasion implications of this phenomenon: the nonnegligible probability of outbreak even when direct transmission is absent, the long-term infectivity of locations of prior outbreaks, and the role of environmental heterogeneity in risk.

Scaling rules for the final decline to extinction

Space-time scaling rules are ubiquitous in ecological phenomena. Current theory postulates three scaling rules that describe the duration of a population's final decline to extinction, although these predictions have not previously been empirically confirmed. We examine these scaling rules across a broader set of conditions, including a wide range of density-dependent patterns in the underlying population dynamics. We then report on tests of these predictions from experiments using the cladoceran Daphnia magna as a model. Our results support two predictions that: (i) the duration of population persistence is much greater than the duration of the final decline to extinction and (ii) the duration of the final decline to extinction increases with the logarithm of the population's estimated carrying capacity. However, our results do not support a third prediction that the duration of the final decline scales inversely with population growth rate. These findings not only support the current standard theory of population extinction but also introduce new empirical anomalies awaiting a theoretical explanation.

The role of environmental transmission in recurrent avian influenza epidemics

Avian influenza virus (AIV) persists in North American wild waterfowl, exhibiting major outbreaks every 2-4 years. Attempts to explain the patterns of periodicity and persistence using simple direct transmission models are unsuccessful. Motivated by empirical evidence, we examine the contribution of an overlooked AIV transmission mode: environmental transmission. It is known that infectious birds shed large concentrations of virions in the environment, where virions may persist for a long time. We thus propose that, in addition to direct fecal/oral transmission, birds may become infected by ingesting virions that have long persisted in the environment. We design a new host-pathogen model that combines within-season transmission dynamics, between-season migration and reproduction, and environmental variation. Analysis of the model yields three major results. First, environmental transmission provides a persistence mechanism within small communities where epidemics cannot be sustained by direct transmission only (i.e., communities smaller than the critical community size). Second, environmental transmission offers a parsimonious explanation of the 2-4 year periodicity of avian influenza epidemics. Third, very low levels of environmental transmission (i.e., few cases per year) are sufficient for avian influenza to persist in populations where it would otherwise vanish.

Complex population dynamics and control of the invasive biennial Alliaria petiolata (garlic mustard)

Controlling species invasions is a leading problem for applied ecology. While controlling populations expanding linearly or exponentially is straightforward, intervention in systems with complex dynamics can have complicated, and sometimes counterintuitive, consequences. Most invasive plant populations are stage-structured and density-dependent--a recipe for complex dynamics--and yet few population models have been created to explore the effects of control efforts on such species. We examined the demography of the invasive biennial plant Alliaria petiolata (garlic mustard) on the front of its spread into a natural area and found evidence of strong density dependence in vital rates of first-year rosette and second-year adult stage classes. We parameterized a density-dependent, stage-structured projection model using field-collected data. This model produces two-point cycles with alternating years in which adults vs. rosettes are more prevalent. Such population dynamics match observations in natural populations, suggesting that these complicated population dynamics may result from deterministic rules. We used this model to evaluate simulated management strategies, including herbicide treatment of rosettes and clipping or pulling of adult plants. Management of A. petiolata by inducing mortality of either rosettes or adults will not be effective at reducing population density unless the induced mortality is very high (>95% for rosettes and >85% for adults) and repeated every year. Indeed, induced mortality of rosettes can be counterproductive, causing increases in the stationary distribution of A. petiolata density. This species is typical of many invasive plants (stage-structured, short-lived, high fertility) and exhibits common forms of density dependence. Thus, the management implications of our study should apply broadly to other species with similar life histories. We suggest that management should focus on managing adults rather than rosettes, and on creating efficient control in targeted areas of the population, rather than spreading less efficient efforts widely.

Environmental variation, stochastic extinction, and competitive coexistence

Understanding how environmental fluctuations affect population persistence is essential for predicting the ecological impacts of expected future increases in climate variability. However, two bodies of theory make opposite predictions about the effect of environmental variation on persistence. Single-species theory, common in conservation biology and population viability analyses, suggests that environmental variation increases the risk of stochastic extinction. By contrast, coexistence theory has shown that environmental variation can buffer inferior competitors against competitive exclusion through a storage effect. We reconcile these two perspectives by showing that in the presence of demographic stochasticity, environmental variation can increase the chance of extinction while simultaneously stabilizing coexistence. Our stochastic simulations of a two-species storage effect model reveal a unimodal relationship between environmental variation and coexistence time, implying maximum coexistence at intermediate levels of environmental variation. The unimodal pattern reflects the fact that the stabilizing influence of the storage effect accumulates rapidly at low levels of environmental variation, whereas the risk of extinction due to the combined effects of environmental variation and demographic stochasticity increases most rapidly at higher levels of variation. Future increases in environmental variation could either increase or decrease an inferior competitor's expected persistence time, depending on the distance between the present level of environmental variation and the optimal level anticipated by this theory.

Pandora's Box Contained Bait: The Global Problem of Introduced Earthworms

Introduced exotic earthworms now occur in every biogeographic region in all but the driest or coldest habitat types on Earth. The global distribution of a few species (e.g., Pontoscolex corethrurus) was noted by early naturalists, but now approximately 120 such peregrine species are recognized to be widespread from regional to global scales, mainly via human activities. Species adapted to human transport and to colonization of disturbed habitats are most widespread and are the principal invasive species. We identify a number of endogenous and exogenous factors that may contribute to the successful establishment and spread of peregrine species. Quantification of these factors may help to determine why certain species become invasive while others do not. Recent advances in theory and modeling of biological invasions and in molecular techniques should prove fruitful in improving our understanding of invasive earthworms, as well as in predicting their impacts on ecosystems.

Effects of habitat quality and size on extinction in experimental populations

Stochastic population theory makes clear predictions about the effects of reproductive potential and carrying capacity on characteristic time-scales of extinction. At the same time, the effects of habitat size and quality on reproduction and regulation have been hotly debated. To trace the causal relationships among these factors, we looked at the effects of habitat size and quality on extinction time in experimental populations of Daphnia magna. Replicate model systems representative of a broad-spectrum consumer foraging on a continuously supplied resource were established under crossed treatments of habitat size (two levels) and habitat quality (three levels) and monitored until eventual extinction of all populations. Using statistically derived estimates of key parameters, we related experimental treatments to persistence time through their effect on carrying capacity and the population growth rate. We found that carrying capacity and the intrinsic rate of increase were each influenced similarly by habitat size and quality, and that carrying capacity and the intrinsic rate of increase were in turn both correlated with time to population extinction. We expected habitat quality to have a greater influence on extinction. However, owing to an unexpected effect of habitat size on reproductive potential, habitat size and quality were similarly important for population persistence. These results support the idea that improving the population growth rate or carrying capacity will reduce extinction risk and demonstrate that both are possible by improving habitat quality or increasing habitat size.