Spatial patterns in the contribution of biotic and abiotic factors to the population dynamics of three freshwater fish species

BACKGROUND

Population dynamics are driven by a number of biotic (e.g., density-dependence) and abiotic (e.g., climate) factors whose contribution can greatly vary across study systems (i.e., populations). Yet, the extent to which the contribution of these factors varies across populations and between species and whether spatial patterns can be identified has received little attention.

METHODS

Here, we used a long-term (1982-2011), broad scale (182 sites distributed across metropolitan France) dataset to study spatial patterns in the population's dynamics of three freshwater fish species presenting contrasted life-histories and patterns of elevation range shifts in recent decades. We used a hierarchical Bayesian approach together with an elasticity analysis to estimate the relative contribution of a set of biotic (e.g., strength of density dependence, recruitment rate) and abiotic (mean and variability of water temperature) factors affecting the site-specific dynamic of two different size classes (0⁺ and >0⁺ individuals) for the three species. We then tested whether the local contribution of each factor presented evidence for biogeographical patterns by confronting two non-mutually exclusive hypotheses: the "range-shift" hypothesis that predicts a gradient along elevation or latitude and the "abundant-center" hypothesis that predicts a gradient from the center to the edge of the species' distributional range.

RESULTS

Despite contrasted life-histories, the three species displayed similar large-scale patterns in population dynamics with a much stronger contribution of biotic factors over abiotic ones. Yet, the contribution of the different factors strongly varied within distributional ranges and followed distinct spatial patterns. Indeed, while abiotic factors mostly varied along elevation, biotic factors-which disproportionately contributed to population dynamics-varied along both elevation and latitude.

CONCLUSIONS

Overall while our results provide stronger support for the range-shift hypothesis, they also highlight the dual effect of distinct factors on spatial patterns in population dynamics and can explain the overall difficulty to find general evidence for geographic gradients in natural populations. We propose that considering the separate contribution of the factors affecting population dynamics could help better understand the drivers of abundance-distribution patterns.

Inter-specific variability in demographic processes affects abundance-occupancy relationships

Species with large local abundances tend to occupy more sites. One of the mechanisms proposed to explain this widely reported inter-specific relationship is a cross-scale hypothesis based on dynamics at the population level. Called the vital rates mechanism; it uses within-population demographic processes of population growth and density dependence to predict when inter-specific abundance-occupancy relationships can arise and when these relationships can weaken and even turn negative. Even though the vital rates mechanism is mathematically simple, its predictions has never been tested directly because of the difficulty estimating the demographic parameters involved. Here, using a recently introduced mark-recapture analysis method, we show that there is no relationship between abundance and occupancy among 17 bird species. Our results are consistent with the predictions of the vital rate mechanism regarding the demographic processes that are expected to weaken this relationship. Specifically, we find that intrinsic growth rate and local abundance are not correlated, and density dependence strength shows considerable variation across species. Variability in density dependence strength is related to variability in species-level local average abundance and intrinsic growth rate; species with lower growth rate have higher abundance and are strongly regulated by density dependent processes, especially acting on survival rates. More generally, our findings support a cross-scale mechanism of macroecological abundance-occupancy relationship emerging from density-dependent dynamics at the population level.

Gut microbiota reflect the crowding stress of space shortage, physical and non-physical contact in Brandt's voles (Lasiopodomys brandtii)

Density-dependence plays a critical role in behavior and population regulation of small mammals, which is likely mediated by hormones and gut microbiota. High density-induced crowding effects often cause a combination of various social stresses including space shortage, physical contact and non-physical contact, but their distinct effects on gut microbiota in animals have not been investigated. In this study, we examined the crowding effects of space shortage and physical or non-physical contact stress on serum corticosterone and gut microbiota of Brandt's voles in both laboratory and field conditions. Our results demonstrated that the space shortage stress showed a more predominant impact on serum corticosterone and gut microbiota of voles than physical or non-physical contact stress; the crowding effects of non-physical contact stress became stronger in high density conditions, while physical contact stress was stronger in a larger group without density effects. High density or group size treatments under both laboratory and semi-natural enclosure conditions significantly increased the relative abundance of key differential taxa, including Bacteroidetes, TM7, S24_7, Streptococcus, and Lactobacillus; while high density or group size treatments decreased the relative abundance of Firmicutes, Staphylococcaceae, Bacteroides, Faecalibacterium, and Adlercreutzia. Our study suggests that high density-induced space shortage and physical contact or non-physical contact stress may play a significant role in behavior and population regulation through altering gut microbiota in small mammals. Our results may also have significant implications in rodent control or health management for livestock.

Long-term patterns in winter habitat selection, breeding and predation in a density-fluctuating, high Arctic lemming population

Habitat selection is expected to balance benefits and costs that maximizes fitness. Using a rare data set on collared lemming (Dicrostonyx groenlandicus) winter nest location spanning more than two decades, we show that lemmings actively select for Salix snow beds, likely due to its favorable micro-climate, and that lemming habitat selection was density-dependent. Lemmings nevertheless exhibited some flexibility in their habitat selection, which appeared to be influenced by the year-to-year variation in snow conditions. The likelihood of both lemming breeding and nest predation by stoats (Mustela erminea) was not directly linked to habitat despite a delicate interplay between habitat, nest size, breeding, and predation. Hence, the larger lemming nests were found in Salix snow beds, and these were more often used for breeding, but both larger nests and nests used for breeding were also predated more often than other nests. Our study provides a clear example of how density-dependent habitat selection acts to balance fitness in the various habitats utilized by collared lemmings.

Density-dependent condition of juvenile penaeid shrimps in seagrass-dominated aquatic vegetation beds located at different distance from a tidal inlet

Seagrasses are critical habitats for the recruitment and growth of juvenile penaeid shrimps within estuaries and coastal lagoons. The location of a seagrass bed within the lagoon can determine the value of a particular bed for shrimp populations. Consequently, differences in the abundance of shrimp can be found in seagrasses depending on their location. As shrimp density increases, density-dependent effects on biological parameters are more likely to occur. However, knowledge about density-dependent processes on shrimp populations in nursery habitats remains limited. The present investigation was undertaken to examine the effects of population density on shrimp condition in two selected seagrass beds, located at different distance from a tidal inlet, one 25 km away (distant) and the other 1 km away (nearby), in a subtropical coastal lagoon. The study was based on monthly samplings during one year in Laguna Madre (Mexico), performing a total of 36 shrimp trawls (100 m2 each one) within each seagrass bed (n = 3 trawls per bed per month for 12 months). Shrimp density was related to the proximity to the tidal inlet (higher density was consistently observed in the nearby seagrass bed), which in turn adversely affected the condition of both species studied (Penaeus aztecus and P. duorarum). In this regard, the present study provides the first evidence of density-dependent effects on shrimp condition inhabiting a nursery habitat. Both shrimp species exhibited a negative relationship between condition and shrimp density. However, this pattern differed depending on the proximity to the tidal inlet, suggesting that shrimp populations inhabiting the nearby seagrass bed are exposed to density-dependent effects on condition; whereas, such effects were not detected in the distant seagrass bed. Shrimp density within the distant seagrass bed was probably below carrying capacity, which is suggested by the better shrimp condition observed in that area of the lagoon. Intra and interspecific competition for food items is surmised to occur, predominantly within seagrass beds near the tidal inlet. However, this hypothesis needs to be tested in future studies.

Intraspecific interactions affect outcomes of pulse toxicity at different Daphnia magna population phases

Traditional toxicity tests assess stressor effects on individuals, while protection goals are focused on the population-level and above. Additionally, these tests ignore common ecological factors such as resource levels and population growth phase. The objective of this research was to explore effects of - and interactions between - resource availability and stress response at the individual and population levels using Daphnia magna as a model. We hypothesized that density-dependent changes in resources at various phases of population growth would cause different population responses to the same toxicant stress. Laboratory populations of Daphnia magna were exposed to a 48-h pulse of 20 or 30 μg/l pyraclostrobin in one of four distinct phases of laboratory population cycles: growth, peak, decline, and stable. Population size and recovery were observed throughout the 51-day study. Populations exposed to pyraclostrobin during the growth phase had the least mortality and fastest recovery, while populations in the peak phase had the greatest mortality and slowest recovery. These data suggested that high density and low food at the peak phase resulted in more sensitive daphnids. To further test this hypothesis, a resource-amended acute toxicity study was conducted to quantify the effects of food resource on pyraclostrobin toxicity to Daphnia magna. Three age classes of Daphnia magna (neonate, subadult, adult) were fed low or high food levels and exposed to pyraclostrobin for 48 h. Toxicity was greater, as shown by lower 48 h LC50s, for smaller Daphnia magna age classes and lower food levels comporting results in the population study. Importantly, the acute toxicity studies generally yielded lower effect levels than the population studies suggesting that while the standard acute studies are ecologically unrealistic, they may be protective of toxicity under some circumstances. Collectively, these data point to the importance of population phase and the resource environment in modulating toxicity.

A computational study of density-dependent individual movement and the formation of population clusters in two-dimensional spatial domains

The patterns of collective behaviour in a population emerging from individual animal movement have long been of interest to ecologists, as has the emergence of heterogeneous patterns among a population. In this paper we will consider these phenomena by using an individual-based modelling approach to simulate a population whose individuals undergo density-dependent movement in 2D spatial domains. We first show that the introduction of density-dependent movement in the form of two parameters, a perception radius and a probability of directed movement, leads to the formation of clusters. We then show that the properties of the clusters and their stability over time are different between populations of Brownian and non-Brownian walkers and are also dependent on the choice of parameters. Finally, we consider the effect of the probability of directed movement on the temporal stability of clusters and show that while clusters formed by Brownian and non-Brownian walkers may have similar properties with certain parameter sets, the spatio-temporal dynamics remain different.

High housing density increases stress hormone- or disease-associated fecal microbiota in male Brandt's voles (Lasiopodomys brandtii)

Density-dependence is an important mechanism in the population regulation of small mammals. Stressors induced by high-density (e.g., crowding and aggression) can cause physiological and neurological disorders, and are hypothesized to be associated with alterations in gut microbiota, which may in turn reduce the fitness of animals by increasing stress- or disease-associated microbes. In this study, we examined the effects of housing density on the hormone levels, immunity, and composition of gut microbiota in male Brandt's voles (Lasiopodomys brandtii) by conducting two specific housing density experiments with or without physical contact between voles. Voles in high density groups exhibited higher serum corticosterone (CORT), serotonin (5-HT), and immunoglobulin G (IgG) levels, as well as higher testosterone (T) levels only in the experiment with physical contact. Meanwhile, high-density treatments induced significant changes in the composition of gut microbiota by increasing disease-associated microbes. The levels of hormones and immunity (i.e., CORT, 5-HT, and IgG) elevated by the high density treatment were significantly correlated with some specific microbes. These results imply that high-density-induced stress may shape the fitness of animals under natural conditions by altering their gut microbiota. Our study provides novel insights into the potential roles of gut microbiota in the density-dependent population regulation of small rodents as well as the potential mechanisms underlying psychological disorders in humans and animals under crowded conditions.

How parasite exposure and time interact to determine Australapatemon burti (Trematoda: Digenea) infections in second intermediate hosts (Erpobdella microstoma) (Hirudinea: Erpodellidae)

Australapatemon spp. are cosmopolitan trematodes that infect freshwater snails, aquatic leeches, and birds. Despite their broad geographic distribution, relatively little is known about interactions between Australapatemon spp. and their leech hosts, particularly under experimental conditions and in natural settings. We used experimental exposures to determine how Australapatemon burti cercariae dosage (number administered to leech hosts, Erpobdella microstoma) affected infection success (fraction to encyst as metacercariae), infection abundance, host survival, and host size over the 100 days following exposure. Interestingly, infection success was strongly density-dependent, such that there were no differences in metacercariae load even among hosts exposed to a 30-fold difference in cercariae. This relationship suggests that local processes (e.g., resource availability, interference competition, or host defenses) may play a strong role in parasite transmission. Our results also indicated that metacercariae did not become evident until ~4 weeks post exposure, with average load climbing until approximately 13 weeks. There was no evidence of metacercariae death or clearance over the census period. Parasite exposure had no detectable effects on leech size or survival, even with nearly 1,000 cercariae. Complementary surveys of leeches in California revealed that 11 of 14 ponds supported infection by A. burti (based on morphology and molecular sequencing), with an average prevalence of 32% and similar metacercariae intensity as in our experimental exposures. The extended development time and extreme density dependence of A. burti has implications for studying naturally occurring host populations, for which detected infections may represent only a fraction of cercariae to which animals have been exposed. Future investigation of these underlying mechanisms would be benefical in understanding host-parasite relationships.

Post-dispersal factors influence recruitment patterns but do not override the importance of seed limitation in populations of a native thistle

Whether plant populations are limited by seed or microsite availability is a long-standing debate. However, since both can be important, increasing emphasis is placed on disentangling their relative importance and how they vary through space and time. Although uncommon, seed addition studies that include multiple levels of seed augmentation, and follow plants through to the adult stage, are critical to achieving this goal. Such data are also vital to understanding when biotic pressures, such as herbivory, influence plant abundance. In this study, we experimentally added seeds of a native thistle, Cirsium canescens, at four augmentation densities to plots at two long-term study sites and quantified densities of seedlings and reproductive adults over 9 years. Recruitment to both seedling and adult stages was strongly seed-limited at both sites; however, the relative strength of seed limitation decreased with plant age. Fitting alternative recruitment functions to our data indicated that post-dispersal mortality factors were important as well. Strong density-dependent mortality limited recruitment at one site, while density-independent limitation predominated at the other. Overall, our experimental seed addition demonstrates that the environment at these sites remains suitable for C. canescens survival to reproduction and that seed availability limits adult densities. The results thus provide support for the hypothesis that seed losses due to the invasive weevil, Rhinocyllus conicus, rather than shifting microsite conditions, are driving C. canescens population declines. Shifts in the importance of density-dependent recruitment limitation between sites highlights that alternate strategies may be necessary to recover plant populations at different locations.

Estimating tsetse fertility: daily averaging versus periodic larviposition

When computing mean daily fertility in adult female tsetse, the common practice of taking the reciprocal of the interlarval period (called averaged fertility) was compared with the method of taking the sum of the products of daily fertility and adult survivorship divided by the sum of daily survivorships (called periodic fertility). The latter method yielded a consistently higher measure of fertility (approximately 10% for tsetse) than the former method. A conversion factor was calculated to convert averaged fertility to periodic fertility. A feasibility criterion was determined for the viability of a tsetse population. Fertility and survivorship data from tsetse populations on Antelope Is. and Redcliff Is., both in Zimbabwe, were used to illustrate the feasibility criterion, as well as the limitations imposed by survivorship and fertility on the viability of tsetse populations. The 10% difference in fertility between the two methods of calculation makes the computation of population feasibility with some parameter combinations sometimes result in a wrong answer. It also underestimates both sterile male release rates required to eradicate a pest population, as well as the speed of resurgence if an eradication attempt fails.

Population dynamics of small endotherms under global change: Greater white-toothed shrews Crocidura russula in Mediterranean habitats

Small endotherms would be especially exposed to main global change drivers (habitat and climate changes) but would also be able to withstand them by adjusting population dynamics locally to changing climate- and habitat-driven food and predation conditions. We analyse the relative importance of changes in climate (mean and variability, including relevant time-lags) and habitat conditions on the abundance, age structure and growth rate of Mediterranean populations of a small endotherm, the greater white-toothed shrew Crocidura russula, along a 10-year period (2008-2017). Habitat type and season were the key factors shaping shrew population dynamics, which showed consistent peak numbers in open habitats in autumn, after the spring-summer reproductive period. Significant increases in aridity (increasing temperature and decreasing rainfall) along the study period did not explain variation in shrew numbers, although short-term variations in abundance were negatively related to relative air humidity and temperature over three last months prior to the surveys. Overall, ongoing climate change have not yet affected shrew population dynamics in its core areas of the Mediterranean region, in spite of expectations based on climate change rate in this region and small endotherm sensitivity to these changes. Reliance on open habitats with lower predation pressure would explain the resilience of shrew populations to climate change. However, current trends of land use change (land abandonment and afforestation) threaten Mediterranean open habitats, so that resilience would not last for long if these trends are not counteracted.

Smallest Anopheles farauti occur during the peak transmission season in the Solomon Islands

BACKGROUND

Malaria transmission varies in intensity amongst Solomon Island villages where Anopheles farauti is the only vector. This variation in transmission intensity might be explained by density-dependent processes during An. farauti larval development, as density dependence can impact adult size with associated fitness costs and daily survivorship.

METHODS

Adult anophelines were sampled from six villages in Western and Central Provinces, Solomon Islands between March 2014 and February 2017. The size of females was estimated by measuring wing lengths, and then analysed for associations with biting densities and rainfall.

RESULTS

In the Solomon Islands, three anopheline species, An. farauti, Anopheles hinesorum and Anopheles lungae, differed in size. The primary malaria vector, An. farauti, varied significantly in size among villages. Greater rainfall was directly associated with higher densities of An. farauti biting rates, but inversely associated with body size with the smallest mean sized mosquitoes present during the peak transmission period. A measurable association between body size and survivorship was not found.

CONCLUSIONS

Density dependent effects are likely impacting the size of adult An. farauti emerging from a range of larval habitats. The data suggest that rainfall increases An. farauti numbers and that these more abundant mosquitoes are significantly smaller in size, but without any reduced survivorship being associated with smaller size. The higher malaria transmission rate in a high malaria focus village appears to be determined more by vector numbers than size or survivorship of the vectors.

Effect of density-dependent individual movement on emerging spatial population distribution: Brownian motion vs Levy flights

Individual animal movement has been a focus of intense research and considerable controversy over the last two decades, however the understanding of wider ecological implications of various movement behaviours is lacking. In this paper, we consider this issue in the context of pattern formation. Using an individual-based modelling approach and computer simulations, we first show that density dependence ("auto-taxis") of the individual movement in a population of random walkers typically results in the formation of a strongly heterogeneous population distribution consisting of clearly defined animal clusters or patches. We then show that, when the movement takes place in a large spatial domain, the properties of the clusters are significantly different in the populations of Brownian and non-Brownian walkers. Whilst clusters tend to be stable in the case of Brownian motion, in the population of Levy walkers clusters are dynamical so that the number of clusters fluctuates in the course of time. We also show that the population dynamics of non-Brownian walkers exhibits two different time scales: a short time scale of the relaxation of the initial condition and a long time scale when one type of dynamics is replaced by another. Finally, we show that the distribution of sample values in the populations of Brownian and non-Brownian walkers is significantly different.

Modeling effects of ecological factors on evolution of polygenic pesticide resistance

Widespread use of pesticides has resulted in the evolution of resistance in many insect pests worldwide limiting their use in pest control. Effective pest and resistance management practices require understanding of the genetics of resistance and of the life history of the pest. Most models for pesticide resistance assume that resistance is monogenic, conferred by a single gene. However, resistance could evolve as a polygenic quantitative trait resulting from the action of several genes, especially when pesticide dose is low. Further, fitness of the pest could be density-dependent and might depend upon abiotic factors such as temperature. It is not known how these factors affect the evolution of polygenic resistance or pest population dynamics when resistance evolves. We use the western corn rootworm, Diabrotica virgifera virgifera, as a case study and use data on density-dependent survival, heritability and survival rates on the transgenic Cry3Bb1 toxin and corresponding LC₅₀ values, to model polygenic resistance to Cry3Bb1. We found that LC₅₀ increased rapidly even at doses that produced a mortality of less than 99.9%. However, survival reached 100% only when mortality was as high as 99.9%. Fast response to high selection pressure produced cyclical larval densities while lower selection pressures produced equilibrium densities. Interestingly we found that a relatively low density observed in a population may not be evidence for a low survival to the pesticide. As a consequence we found that larger refuges might not necessarily help in reducing pest densities especially when pesticide mortality is low. This effect, arising from the tradeoff between response to selection and density dependence, calls for careful assessment of the evolution of resistance based on change in survival as well as on pest densities. When selection pressure is low a lower initial density resulted in a larger response to selection. Finally, we showed that populations with shorter developmental times developed resistance faster initially irrespective of selection pressure. However, when selection pressure is low survival eventually became higher in populations with longer developmental times. Since developmental time depends on degree days spatio-temporal variation in temperature could be an important factor in polygenic resistance evolution.

Population response of an apex Antarctic consumer to its prey and climate fluctuations

A fundamental endeavor in population ecology is to identify the drivers of population dynamics. A few empirical studies included the effect of prey abundance when investigating simultaneously the effects of density-dependence and climate factors on marine top-predator population dynamics. Our aim was to unravel the mechanisms forcing population dynamics of an apex consumer seabird, the south polar skua, using long-term climatic and population time series of the consumer and its prey in Terre Adélie, Antarctica. Influences of density-dependence, climatic factors, and prey abundance with lag effects were tested on the breeding population dynamics with a Bayesian multi-model inference approach. We evidenced a negative trend in breeding population growth rate when density increased. Lagged effects of sea-ice concentration and air temperature in spring and a contemporary effect of prey resources were supported. Remarkably, results outline a reverse response of the south polar skua and one of its main preys to the same environmental factor (sea-ice concentration), suggesting a strong link between skua and penguin dynamics. The causal mechanisms may involve competition for food and space through territorial behavior as well as local climate and prey availability, which probably operate on breeding parameters (breeding propensity, breeding success, or recruitment) rather than on adult survival. Our results provide new insights on the relative importance of factors forcing the population dynamics of an apex consumer including density-dependence, local climate conditions, and direct and indirect effects of prey abundance.

Extrinsic Mortality Can Shape Life-History Traits, Including Senescence

The Williams' hypothesis is one of the most widely known ideas in life history evolution. It states that higher adult mortality should lead to faster and/or earlier senescence. Theoretically derived gradients, however, do not support this prediction. Increased awareness of this fact has caused a crisis of misinformation among theorists and empirical ecologists. We resolve this crisis by outlining key issues in the measurement of fitness, assumptions of density dependence, and their effect on extrinsic mortality. The classic gradients apply only to a narrow range of ecological contexts where density-dependence is either absent or present but with unrealistic stipulations. Re-deriving the classic gradients, using a more appropriate measure of fitness and incorporating density, shows that broad ecological contexts exist where Williams' hypothesis is supported.

Nonlinear population dynamics in a bounded habitat

A key issue in ecology is whether a population will survive long term or go extinct. This is the question we address in this paper for a population in a bounded habitat. We will restrict our study to the case of a single species in a one-dimensional habitat of length L. The evolution of the population density distribution ρ(x, t), where x is the position and t the time, is governed by elementary processes such as growth and dispersal, which, in standard models, are typically described by a constant per capita growth rate and normal diffusion, respectively. However, feedbacks in the regulatory mechanisms and external factors can produce density-dependent rates. Therefore, we consider a generalization of the standard evolution equation, which, after dimensional scaling and assuming large carrying capacity, becomes ∂_tρ=∂_x(ρ^ν-1∂_xρ)+ρ^μ, where μ,ν∈R. This equation is complemented by absorbing boundaries, mimicking adverse conditions outside the habitat. For this nonlinear problem, we obtain, analytically, exact expressions of the critical habitat size L_c for population survival, as a function of the exponents and initial conditions. We find that depending on the values of the exponents (ν, μ), population survival can occur for either L ≥ L_c, L ≤ L_c or for any L. This generalizes the usual statement that L_c represents the minimum habitat size. In addition, nonlinearities introduce dependence on the initial conditions, affecting L_c.

Phenotypic and genotypic monitoring of Schistosoma mansoni in Tanzanian schoolchildren five years into a preventative chemotherapy national control programme

BACKGROUND

Schistosoma mansoni is a parasite of profound medical importance. Current control focusses on mass praziquantel (PZQ) treatment of populations in endemic areas, termed Preventative Chemotherapy (PC). Large-scale PC programmes exert prolonged selection pressures on parasites with the potential for, direct and/or indirect, emergence of drug resistance. Molecular methods can help monitor genetic changes of schistosome populations over time and in response to drug treatment, as well as estimate adult worm burdens through parentage analysis. Furthermore, methods such as in vitro drug sensitivity assays help phenotype in vivo parasite genotypic drug efficacy.

METHODS

We conducted combined in vitro PZQ efficacy testing with population genetic analyses of S. mansoni collected from children from two schools in 2010, five years after the introduction of a National Control Programme. Children at one school had received four annual PZQ treatments and the other school had received two mass treatments in total. We compared genetic differentiation, indices of genetic diversity, and estimated adult worm burden from parasites collected in 2010 with samples collected in 2005 (before the control programme began) and in 2006 (six months after the first PZQ treatment). Using 2010 larval samples, we also compared the genetic similarity of those with high and low in vitro sensitivity to PZQ.

RESULTS

We demonstrated that there were individual parasites with reduced PZQ susceptibility in the 2010 collections, as evidenced by our in vitro larval behavioural phenotypic assay. There was no evidence, however, that miracidia showing phenotypically reduced susceptibility clustered together genetically. Molecular analysis also demonstrated a significant reduction of adult worm load over time, despite little evidence of reduction in parasite infection intensity, as measured by egg output. Genetic diversity of infections did not reduce over time, despite changes in the genetic composition of the parasite populations.

CONCLUSIONS

Genotypic and phenotypic monitoring did not indicate a selective sweep, as may be expected if PZQ treatment was selecting a small number of related "resistant" parasites, but there was evidence of genetic changes at the population level over time. Genetic data were used to estimate adult worm burdens, which unlike parasite infection intensity, showed reductions over time, suggesting the relaxation of negative density-dependent constraints on parasite fecundity with PZQ treatment. We thereby demonstrated that density-dependence in schistosome populations may complicate evaluation and monitoring of control programmes.

Influence of soil pathogens on early regeneration success of tropical trees varies between forest edge and interior

Soil fungi are key mediators of negative density-dependent mortality in seeds and seedlings, and the ability to withstand pathogens in the shaded understory of closed-canopy forests could reinforce light gradient partitioning by tree species. For four species of tropical rainforest trees-two shade-tolerant and two shade-intolerant-we conducted a field experiment to examine the interactive effects of fungal pathogens, light, and seed density on germination and early seedling establishment. In a fully factorial design, seeds were sown into 1 m² plots containing soil collected from underneath conspecific adult trees, with plots assigned to forest edge (high light) or shaded understory, high or low density, and fungicide or no fungicide application. We monitored total seed germination and final seedling survival over 15 weeks. Shade-intolerant species were strongly constrained by light; their seedlings survived only at the edge. Fungicide application significantly improved seedling emergence and/or survival for three of the four focal species. There were no significant interactions between fungicide and seed density, suggesting that pathogen spread with increased aggregation of seeds and seedlings did not contribute to pathogen-mediated mortality. Two species experienced significant edge-fungicide interactions, but fungicide effects in edge vs. interior forest varied with species and recruitment stage. Our results suggest that changes to plant-pathogen interactions could affect plant recruitment in human-impacted forests subject to fragmentation and edge-effects.

Density-dependent sex ratio and sex-specific preference for host traits in parasitic bat flies

Background

Deviation of sex ratios from unity in wild animal populations has recently been demonstrated to be far more prevalent than previously thought. Ectoparasites are prominent examples of this bias, given that their sex ratios vary from strongly female- to strongly male-biased both among hosts and at the metapopulation level. To date our knowledge is very limited on how and why these biased sex ratios develop. It was suggested that sex ratio and sex-specific aggregation of ectoparasites might be shaped by the ecology, behaviour and physiology of both hosts and their parasites. Here we investigate a highly specialised, hematophagous bat fly species with strong potential to move between hosts, arguably limited inbreeding effects, off-host developmental stages and extended parental care.

Results

We collected a total of 796 Nycteribia kolenatii bat flies from 147 individual bats using fumigation and subsequently determined their sex. We report a balanced sex ratio at the metapopulation level and a highly variable sex ratio among infrapopulations ranging from 100% male to 100% female. We show that infrapopulation sex ratio is not random and is highly correlated with infrapopulation size. Sex ratio is highly male biased in small and highly female biased in large infrapopulations. We show that this pattern is most probably the result of sex-specific preference in bat flies for host traits, most likely combined with a higher mobility of males. We demonstrate that female bat flies exert a strong preference for high host body condition and female hosts, while the distribution of males is more even.

Conclusions

Our results suggest that locally biased sex ratios can develop due to sex-specific habitat preference of parasites. Moreover, it is apparent that the sex of both hosts and parasites need to be accounted for when a better understanding of host-parasite systems is targeted.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-017-2340-0) contains supplementary material, which is available to authorized users.

Stochastic population growth in spatially heterogeneous environments: the density-dependent case

This work is devoted to studying the dynamics of a structured population that is subject to the combined effects of environmental stochasticity, competition for resources, spatio-temporal heterogeneity and dispersal. The population is spread throughout n patches whose population abundances are modeled as the solutions of a system of nonlinear stochastic differential equations living on \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$[0,\infty )^n$$\end{document}[0,∞)n. We prove that r, the stochastic growth rate of the total population in the absence of competition, determines the long-term behaviour of the population. The parameter r can be expressed as the Lyapunov exponent of an associated linearized system of stochastic differential equations. Detailed analysis shows that if \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ r>0$$\end{document}r>0, the population abundances converge polynomially fast to a unique invariant probability measure on \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$(0,\infty )^n$$\end{document}(0,∞)n, while when \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ r<0$$\end{document}r<0, the population abundances of the patches converge almost surely to 0 exponentially fast. This generalizes and extends the results of Evans et al. (J Math Biol 66(3):423–476, 2013) and proves one of their conjectures. Compared to recent developments, our model incorporates very general density-dependent growth rates and competition terms. Furthermore, we prove that persistence is robust to small, possibly density dependent, perturbations of the growth rates, dispersal matrix and covariance matrix of the environmental noise. We also show that the stochastic growth rate depends continuously on the coefficients. Our work allows the environmental noise driving our system to be degenerate. This is relevant from a biological point of view since, for example, the environments of the different patches can be perfectly correlated. We show how one can adapt the nondegenerate results to the degenerate setting. As an example we fully analyze the two-patch case, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$n=2$$\end{document}n=2, and show that the stochastic growth rate is a decreasing function of the dispersion rate. In particular, coupling two sink patches can never yield persistence, in contrast to the results from the non-degenerate setting treated by Evans et al. which show that sometimes coupling by dispersal can make the system persistent.

A density-dependent model of competition for one resource in the chemostat

This paper deals with a two-microbial species model in competition for a single-resource in the chemostat including general intra- and interspecific density-dependent growth rates with distinct removal rates for each species. In order to understand the effects of intra- and interspecific interference, this general model is first studied by determining the conditions of existence and local stability of steady states. With the same removal rate, the model can be reduced to a planar system and then the global stability results for each steady state are derived. The bifurcations of steady states according to interspecific interference parameters are analyzed in a particular case of density-dependent growth rates which are usually used in the literature. The operating diagrams show how the model behaves by varying the operating parameters and illustrate the effect of the intra- and interspecific interference on the disappearance of coexistence region and the occurrence of bi-stability region. Concerning the small enough interspecific interference terms, we would shed light on the global convergence towards the coexistence steady state for any positive initial condition. When the interspecific interference pressure is large enough this system exhibits bi-stability where the issue of the competition depends on the initial condition.

Decline in territory size and fecundity as a response to carrying capacity in an endangered songbird

Density-dependent processes are fundamental mechanisms for the regulation of populations. Ecological theories differ in their predictions on whether increasing population density leads to individual adjustments of survival and reproductive output or to dominance and monopolization of resources. Here, we use a natural experiment to examine which factors limit population growth in the only remaining population of the endangered pale-headed brush finch (Atlapetes pallidiceps). For three distinct phases (a phase of population suppression, 2001-2002; expansion due to conservation management, 2003-2008; and equilibrium phase, 2009-2014), we estimated demographic parameters with an integrated population model using population size, the proportion of successfully breeding pairs and their productivity, territory size, and mark-recapture data of adult birds. A low proportion of successful breeders due to brood parasitism (0.42, 95% credible interval 0.26-0.59) limited population growth before 2003; subsequent culling of the brood parasite resulted in a two-fold increase of the proportion of successful breeders during the 'expansion phase'. When the population approached the carrying capacity of its habitat, territory size declined by more than 50% and fecundity declined from 1.9 (1.54-2.27) to 1.3 (1.12-1.53) chicks per breeding pair, but the proportion of successful breeders remained constant (expansion phase: 0.85; 0.76-0.93; equilibrium phase: 0.86; 0.79-0.92). This study demonstrates that limiting resources can lead to individual adjustments instead of despotic behavior, and the individual reduction of reproductive output at high population densities is consistent with the slow life-history of many tropical species.

Predicting Wolbachia invasion dynamics in Aedes aegypti populations using models of density-dependent demographic traits

BACKGROUND

Arbovirus transmission by the mosquito Aedes aegypti can be reduced by the introduction and establishment of the endosymbiotic bacteria Wolbachia in wild populations of the vector. Wolbachia spreads by increasing the fitness of its hosts relative to uninfected mosquitoes. However, mosquito fitness is also strongly affected by population size through density-dependent competition for limited food resources. We do not understand how this natural variation in fitness affects symbiont spread, which limits our ability to design successful control strategies.

RESULTS

We develop a mathematical model to predict A. aegypti-Wolbachia dynamics that incorporates larval density-dependent variation in important fitness components of infected and uninfected mosquitoes. Our model explains detailed features of the mosquito-Wolbachia dynamics observed in two independent experimental A. aegypti populations, allowing the combined effects on dynamics of multiple density-dependent fitness components to be characterized. We apply our model to investigate Wolbachia field release dynamics, and show how invasion outcomes can depend strongly on the severity of density-dependent competition at the release site. Specifically, the ratio of released relative to wild mosquitoes required to attain a target infection frequency (at the end of a release program) can vary by nearly an order of magnitude. The time taken for Wolbachia to become established following releases can differ by over 2 years. These effects depend on the relative fitness of field and insectary-reared mosquitoes.

CONCLUSIONS

Models of Wolbachia invasion incorporating density-dependent demographic variation in the host population explain observed dynamics in experimental A. aegypti populations. These models predict strong effects of density-dependence on Wolbachia dynamics in field populations, and can assist in the effective use of Wolbachia to control the transmission of arboviruses such as dengue, chikungunya and zika.

Differences in density-dependence drive dual offspring size strategies in fish

Offspring size reflects the optimal balance between female fecundity and allocation of energy to each offspring. Most fish, in particular teleost species, produce many small eggs, while others, notably elasmobranch species, have low fecundity and large offspring. No general explanation has yet been put forwards to explain these different strategies between species which occupy similar habitats. We approach the problem by (1) examining the differences between life history parameters of teleost fish and elasmobranchs and (2) an evolutionary model. We show that life history parameters characterising growth, mortality and reproductive output are almost similar between teleosts and elasmobranchs. We find that a model which accounts for density-dependence predicts dual offspring size strategies: either invariant with adult size or proportional to adult size. The model predicts that the invariant strategy is associated with an absence of density-dependence in early life whereas proportional offspring are subject to density-dependence throughout life. Parameterising the model using life history data regenerates the observed dual offspring size pattern. We conjecture that the life stage where density-dependent competition occurs is an important factor behind the observed difference in offspring size strategies.

Inter-class competition in stage-structured populations: effects of adult density on life-history traits of adult and juvenile common lizards

Ecological and evolutionary processes in natural populations are largely influenced by the population's stage-structure. Commonly, different classes have different competitive abilities, e.g., due to differences in body size, suggesting that inter-class competition may be important and largely asymmetric. However, experimental evidence states that inter-class competition, which is important, is rare and restricted to marine fish. Here, we manipulated the adult density in six semi-natural populations of the European common lizard, Zootoca vivipara, while holding juvenile density constant. Adult density affected juveniles, but not adults, in line with inter-class competition. High adult density led to lower juvenile survival and growth before hibernation. In contrast, juvenile survival after hibernation was higher in populations with high adult density, pointing to relaxed inter-class competition. As a result, annual survival was not affected by adult density, showing that differences in pre- and post-hibernation survival balanced each other out. The intensity of inter-class competition affected reproduction, performance, and body size in juveniles. Path analyses unravelled direct treatment effects on early growth (pre-hibernation) and no direct treatment effects on the parameters measured after hibernation. This points to allometry of treatment-induced differences in early growth, and it suggests that inter-class competition mainly affects the early growth of the competitively inferior class and thereby their future performance and reproduction. These results are in contrast with previous findings and, together with results in marine fish, suggest that the strength and direction of density dependence may depend on the degree of inter-class competition, and thus on the availability of resources used by the competing classes.

Density-dependent effects of larval dispersal mediated by host plant quality on populations of an invasive insect

The success of invasive species is often thought to be due to release from natural enemies. This hypothesis assumes that species are regulated by top-down forces in their native range and are likely to be regulated by bottom-up forces in the invasive range. Neither of these assumptions has been consistently supported with insects, a group which includes many destructive invasive species. Winter moth (Operophtera brumata) is an invasive defoliator in North America that appears to be regulated by larval mortality. To assess whether regulation was caused by top-down or bottom-up forces, we sought to identify the main causes of larval mortality. We used observational and manipulative field and laboratory studies to demonstrate that larval mortality due to predation, parasitism, and disease were minimal. We measured the response of larval dispersal in the field to multiple aspects of foliar quality, including total phenolics, pH 10 oxidized phenolics, trichome density, total nitrogen, total carbon, and carbon-nitrogen ratio. Tree-level declines in density were driven by density-dependent dispersal of early instars. Late instar larvae dispersed at increased rates from previously damaged as compared to undamaged foliage, and in 2015 field larval dispersal rates were related to proportion of oxidative phenolics. We conclude that larval dispersal is the dominant source of density-dependent larval mortality, may be mediated by induced changes in foliar quality, and likely regulates population densities in New England. These findings suggest that winter moth population densities in New England are regulated by bottom-up forces, aligning with the natural enemy release hypothesis.

Statistical validation of structured population models for Daphnia magna

In this study we use statistical validation techniques to verify density-dependent mechanisms hypothesized for populations of Daphnia magna. We develop structured population models that exemplify specific mechanisms and use multi-scale experimental data in order to test their importance. We show that fecundity and survival rates are affected by both time-varying density-independent factors, such as age, and density-dependent factors, such as competition. We perform uncertainty analysis and show that our parameters are estimated with a high degree of confidence. Furthermore, we perform a sensitivity analysis to understand how changes in fecundity and survival rates affect population size and age-structure.

Post-dispersal seed predation and the establishment of vertebrate dispersed plants in Mediterranean scrublands

The post-dispersal fate of seeds and fruit (diaspores) of three vertebrate-dispersed trees, Crataegus monogyna, Prunus mahaleb and Taxus baccata, was studied in the Andalusian highlands, south-eastern Spain. Exclosures were used to quantify separately the impact of vertebrates and invertebrates on seed removal in relation to diaspore density and microhabitat. The three plant species showed marked differences in the percentage of diaspores removed, ranging from only 5% for C. monogyna to 87% for T. baccata. Although chaffinches (Fringilla coelebs) fed on diaspores, rodents (Apodemus sylvaticus) were the main vertebrate removers of seed and fruit. Two species of ant (Cataglyphis velox and Aphaenogaster iberica) were the only invertebrates observed to remove diaspores. However, the impact of ants was strongly seasonal and they only removed P. mahaleb fruit to any significant extent. While removal of seed by rodents was equivalent to predation, ants were responsible for secondary dispersal. However, their role was limited to infrequent, small-scale redistribution of fruit in the vicinity of parent trees. Rodents and ants differed in their use of different microhabitats. Rodents foraged mostly beneath trees and low shrubs and avoided open areas while the reverse was true of ants. Thus, patterns of post-dispersal seed removal will be contigent on the relative abundance and distribution of ants and rodents. Studies which neglect to quantify separately the impacts of these two guilds of seed removers may fail to elucidate the mechanisms underlying patterns of post-dispersal seed removal. The coincidence of both increased seed deposition by the main avian dispersers (Turdus spp.) and increased seed predation with increasing vegetation height suggested that selection pressures other than post-dispersal seed predation shape the spatial pattern of seed dispersal. Rather than providing a means of escaping post-dispersal seed predators, dispersal appears to direct seeds to microhabitats most suitable for seedling survival. Nevertheless, the reliance of most vertebrate-dispersed trees on regeneration by seed and the absence of persistent soil seed banks imply that post-dispersal seed predators may exert a strong influence on the demography of the plants whose seeds they consume. Even where microsites are limited, the coincidence of the most suitable microhabitats for seedling establishment with those where seed predation is highest provide a means by which selective seed predators can influence community composition.