A bit of sex stabilizes host-parasite dynamics

To date, only a few studies have focused on the effects of sex on population dynamics. Previous models have typically found that sexual reproduction dampens population fluctuations. Although asexual and sexual reproduction are just the two endpoints along a continuum of varying rates of sex, previous work has ignored the effects of intermediate degrees of sex on population dynamics. Here we study the effects of partial sexual reproduction (i.e. sex occurs only every few generations or with small probability in each generation) on the coupled population dynamics of a Nicholson-Bailey host-parasite model. We show that complex dynamics are simplified for high host population growth rates if the frequency of sex is sufficiently high in both host and parasite: sex decreases fluctuations in population density, and leads to non-chaotic dynamics for population growth rates that would result in chaotic dynamics in the absence of sexual reproduction. However, the simplification does not increase gradually with an increasing frequency of sex but appears abruptly at low-to-intermediate frequencies of sex. For some parameter settings, intermediate frequencies of sexual reproduction can simplify the dynamics more than lower or higher frequencies. Thus, in agreement with earlier results, sexual reproduction typically stabilizes complex population dynamics in our models. Additionally, our results suggest that low-to-intermediate frequencies of sex may often be as (or even more) stabilizing as high frequencies.     

Why spacing behavior does not stabilize density in cyclic populations of microtine rodents

Summary There are several published hypotheses that consider spacing behavior to be a significant factor causing the multiannual density fluctuations characteristic of some microtine rodent populations. Recent modeling efforts have concluded, however, that spacing behavior should have a stabilizing rather than a destabilizing effect on population dynamics. Why doesn't spacing behavior stabilize these cyclic populations? We argue that while spacing behavior does have a stabilizing influence on population dynamics by limiting the number of breeding individuals, reproduction continues and population size is not limited in an asymptotic manner. Rather, microtine social organization produces demographic changes within a population that allow density cycles to occur under certain conditions. Using a simulation model, we demonstrate that in a strongly seasonal environment populations with low density dependence in reproduction will cycle whereas populations with high density dependence in reproduction will have relatively stable densities. Given such complicating factors as the annual species nature of microtine rodents, occasionally intense predation, and the tendency for territoriality to break down during the non-breeding season, individuals with low density dependence in reproduction will always be able to invade and eventually dominate populations with high density dependence in reproduction, regardless of the resulting destabilization of population dynamics.     

Density,parasitism, and sexual reproduction are strongly correlated in lake Daphnia populations

Many organisms can reproduce both asexually and sexually. For cyclical parthenogens, periods of asexual reproduction are punctuated by bouts of sexual reproduction, and the shift from asexual to sexual reproduction has large impacts on fitness and population dynamics. We studied populations of Daphnia dentifera to determine the amount of investment in sexual reproduction as well as the factors associated with variation in investment in sex. To do so, we tracked host density, infections by nine different parasites, and sexual reproduction in 15 lake populations of D. dentifera for 3 years. Sexual reproduction was seasonal, with male and ephippial female production beginning as early as late September and generally increasing through November. However, there was substantial variation in the prevalence of sexual individuals across populations, with some populations remaining entirely asexual throughout the study period and others shifting almost entirely to sexual females and males. We found strong relationships between density, prevalence of infection, parasite species richness, and sexual reproduction in these populations. However, strong collinearity between density, parasitism, and sexual reproduction means that further work will be required to disentangle the causal mechanisms underlying these relationships.     

Demography can favour female-advantageous alleles

When female fecundity is relatively independent of male abundance, while male reproduction is proportional to female abundance, females have a larger effect on population dynamics than males (i.e. female demographic dominance). This population dynamic phenomenon might not appear to influence evolution, because male and female genomes still contribute equally much to the next generation. However, here we examine two evolutionary scenarios to provide a proof of principle that spatial structure can make female demographic dominance matter. Our two simulation models combine dispersal evolution with local adaptation subjected to intralocus sexual conflict and environmentally driven sex ratio biases, respectively. Both models have equilibria where one environment (without being intrinsically poorer) has so few reproductive females that trait evolution becomes disproportionately determined by those environments where females survive better (intralocus sexual conflict model), or where daughters are overproduced (environmental sex determination model). Surprisingly, however, the two facts that selection favours alleles that benefit females, and population growth is improved when female fitness is high, together do not imply that all measures of population performance are improved. The sex-specificity of the source–sink dynamics predicts that populations can evolve to fail to persist in habitats where alleles do poorly when expressed in females.     

Sexual reproduction and diversity: Connection between sexual selection and biological communities via population dynamics

Sexual reproduction is a mysterious phenomenon. Most animals and plants invest in sexual reproduction, even though it is more costly than asexual reproduction. Theoretical studies suggest that occasional or conditional use of sexual reproduction, involving facultative switching between sexual and asexual reproduction, is the optimal reproductive strategy. However, obligate sexual reproduction is common in nature. Recent studies suggest that the evolution of facultative sexual reproduction is prevented by males that coerce females into sexual fertilization; thus, sexual reproduction has the potential to enforce costs on a given species. Here, the effect of sex on biodiversity is explored by evaluating the reproductive costs arising from sex. Sex provides atypical selection pressure that favors traits that increase fertilization success, even at the expense of population growth rates, that is, sexual selection. The strength of sexual selection depends on the density of a given species. Sexual selection often causes strong negative effects on the population growth rates of species that occur at high density. Conversely, a species that reduces its density is released from this negative effect, and so increases its growth rate. Thus, this negative density-dependent effect on population growth that arises from sexual selection could be used to rescue endangered species from extinction, prevent the overgrowth of common species and promote the coexistence of competitive species. Recent publications on sexual reproduction provide several predictions related to the evolution of reproductive strategies, which is an important step toward integrating evolutionary dynamics, demographic dynamics and community dynamics.     

The genetics and evolution of obligate reproductive parasitism in Trichogramma pretiosum infected with parthenogenesis-inducing Wolbachia

Parthenogenesis-inducing (PI) Wolbachia belong to a class of intracellular symbionts that distort the offspring sex ratio of their hosts toward a female bias. In many PI Wolbachia-infected species sex ratio distortion has reached its ultimate expression-fixation of infection and all-female populations. This is only possible with thelytokous PI symbionts as they provide an alternative form of reproduction and remove the requirement for males and sexual reproduction. Many populations fixed for PI Wolbachia infection have lost the ability to reproduce sexually, even when cured of the infection. We examine one such population in the species Trichogramma pretiosum. Through a series of backcrossing experiments with an uninfected Trichogramma pretiosum population we were able to show that the genetic basis for the loss of female sexual function could be explained by a dominant nuclear effect. Male sexual function had not been completely lost, though some deterioration of male sexual function was also evident when males from the infected population (created through antibiotic curing of infected females) were mated to uninfected females. We discuss the dynamics of sex ratio selection in PI Wolbachia-infected populations and the evolution of non-fertilizing mutations.     

Sperm storage reduces the strength of the mate‐finding Allee effect

Mate searching is a key component of sexual reproduction that can have important implications for population viability, especially for the mate‐finding Allee effect. Interannual sperm storage by females may be an adaptation that potentially attenuates mate limitation, but the demographic consequences of this functional trait have not been studied. Our goal is to assess the effect of female sperm storage durability on the strength of the mate‐finding Allee effect and the viability of populations subject to low population density and habitat alteration. We used an individual‐based simulation model that incorporates realistic representations of the demographic and spatial processes of our model species, the spur‐thighed tortoise (Testudo graeca). This allowed for a detailed assessment of reproductive rates, population growth rates, and extinction probabilities. We also studied the relationship between the number of reproductive males and the reproductive rates for scenarios combining different levels of sperm storage durability, initial population density, and landscape alteration. Our results showed that simulated populations parameterized with the field‐observed demographic rates collapsed for short sperm storage durability, but were viable for a durability of one year or longer. In contrast, the simulated populations with a low initial density were only viable in human‐altered landscapes for sperm storage durability of 4 years. We find that sperm storage is an effective mechanism that can reduce the strength of the mate‐finding Allee effect and contribute to the persistence of low‐density populations. Our study highlights the key role of sperm storage in the dynamics of species with limited movement ability to facilitate reproduction in patchy landscapes or during population expansion. This study represents the first quantification of the effect of sperm storage durability on population dynamics in different landscapes and population scenarios.     

Density dependence in group dynamics of a highly social mongoose, Suricata suricatta

1.?For social species, the link between individual behaviour and population dynamics is mediated by group-level demography. 2.?Populations of obligate cooperative breeders are structured into social groups, which may be subject to inverse density dependence (Allee effects) that result from a dependence on conspecific helpers, but evidence for population-wide Allee effects is rare. 3.?We use field data from a long-term study of cooperative meerkats (Suricata suricatta; Schreber, 1776) - a species for which local Allee effects are not reflected in population-level dynamics - to empirically model interannual group dynamics. 4.?Using phenomenological population models, modified to incorporate environmental conditions and potential Allee effects, we first investigate overall patterns of group dynamics and find support only for conventional density dependence that increases after years of low rainfall. 5.?To explain the observed patterns, we examine specific demographic rates and assess their contributions to overall group dynamics. Although per-capita meerkat mortality is subject to a component Allee effect, it contributes relatively little to observed variation in group dynamics, and other (conventionally density dependent) demographic rates - especially emigration - govern group dynamics. 6.?Our findings highlight the need to consider demographic processes and density dependence in subpopulations before drawing conclusions about how behaviour affects population processes in socially complex systems.     

Demographic heterogeneity impacts density-dependent population dynamics

Among-individual variation in vital parameters such as birth and death rates that is unrelated to age, stage, sex, or environmental fluctuations is referred to as demographic heterogeneity. This kind of heterogeneity is prevalent in ecological populations, but is almost always left out of models. Demographic heterogeneity has been shown to affect demographic stochasticity in small populations and to increase growth rates for density-independent populations. The latter is due to ??cohort selection,?? where the most frail individuals die out first, lowering the cohort??s average mortality as it ages. The importance of cohort selection to population dynamics has only recently been recognized. We use a continuous-time model with density dependence, based on the logistic equation, to study the effects of demographic heterogeneity in mortality and reproduction. Reproductive heterogeneity is introduced in three ways: parent fertility, offspring viability, and parent?Coffspring correlation. We find that both the low-density growth rate and the equilibrium population size increase as the magnitude of mortality heterogeneity increases or as parent?Coffspring phenotypic correlation increases. Population dynamics are affected by complex interactions among the different types of heterogeneity, and trade-off scenarios are examined which can sometimes reverse the effect of increased heterogeneity. We show that there are a number of different homogeneous approximations to heterogeneous models, but all fail to capture important parts of the dynamics of the full model.     

Sex in space: population dynamic consequences

Sex, so important in the reproduction of bigametic species, is nonetheless often ignored in explorations of the dynamics of populations. Using a growth model of dispersal-coupled populations we can keep track of fluctuations in numbers of females and males. The sexes may differ from each other in their ability to disperse and their sensitivity to population density. As a further complication, the breeding system is either monogamous or polygamous. We use the harmonic mean birth function to account for sex-ratio-dependent population growth in a Moran–Ricker population renewal process. Incorporating the spatial dimension stabilizes the dynamics of populations with monogamy as the breeding system, but does not stabilize the population dynamics of polygamous species. Most notably, in populations coupled with dispersal, where the sexes differ in their dispersal ability there are rarely stable and equal sex ratios. Rather, a two-point cycle, four-point cycle and eventually complex behaviour of sex-ratio dynamics will emerge with increasing birth rates. Monogamy often leads to less noisy sex-ratio dynamics than polygamy. In our model, the sex-ratio dynamics of coupled populations differ from those of an isolated population system, where a stable 50:50 sex ratio is achievable with equal density-dependence costs for females and males. When sexes match in their dispersal ability, population dynamics and sex-ratio dynamics of coupled populations collapse to those of isolated populations.     

Current and historical factors drive variation of reproductive traits in unisexual mosses in Europe: A case study

Unisexual bryophytes provide excellent models to study the mechanisms that regulate the frequency of sexual versusasexual reproduction in plants, and their ecological and evolutionary implications. Here, we determined sex expression, phenotypic sex ratio, and individual shoot traits in 242 populations of the cosmopolitan moss Pseudoscleropodium purum spanning its whole distributional range. We tested whether niche differentiation, sex-specific differences in shoot size, and biogeographical history explained the spatial variation of reproductive traits. We observed high levels of sex expression and predominantly female-biased populations, although both traits showed high intraspecific variation among populations. Sex expression and sex ratio were partly explained by current macroscale environmental variation, with male shoots being less frequent at the higher end of the environmental gradients defined by the current distribution of the species. Female bias in population sex ratio was significantly lower in areas recolonized after the last glacial maximum (recent populations) than in glacial refugia (long-term persistent populations). We demonstrated that reproductive trait variation in perennial unisexual mosses is partially driven by macroscale and historical environmental variation. Based on our results, we hypothesize that sexual dimorphism in environmental tolerance and vegetative growth contribute to sex ratio bias over time, constraining the chances of sexual reproduction, especially in long-term persistent populations. Further studies combining genetic analyses and population monitoring should improve our understanding of the implications of the intraspecific variation in the frequency of sexual versusasexual reproduction in bryophyte population fitness and eco-evolutionary dynamics.     

Approximation of a physiologically structured population model with seasonal reproduction by a stage-structured biomass model

Seasonal reproduction causes, due to the periodic inflow of young small individuals in the population, seasonal fluctuations in population size distributions. Seasonal reproduction furthermore implies that the energetic body condition of reproducing individuals varies over time. Through these mechanisms, seasonal reproduction likely affects population and community dynamics. While seasonal reproduction is often incorporated in population models using discrete time equations, these are not suitable for size-structured populations in which individuals grow continuously between reproductive events. Size-structured population models that consider seasonal reproduction, an explicit growing season and individual-level energetic processes exist in the form of physiologically structured population models. However, modeling large species ensembles with these models is virtually impossible. In this study, we therefore develop a simpler model framework by approximating a cohort-based size-structured population model with seasonal reproduction to a stage-structured biomass model of four ODEs. The model translates individual-level assumptions about food ingestion, bioenergetics, growth, investment in reproduction, storage of reproductive energy, and seasonal reproduction in stage-based processes at the population level. Numerical analysis of the two models shows similar values for the average biomass of juveniles, adults, and resource unless large-amplitude cycles with a single cohort dominating the population occur. The model framework can be extended by adding species or multiple juvenile and/or adult stages. This opens up possibilities to investigate population dynamics of interacting species while incorporating ontogenetic development and complex life histories in combination with seasonal reproduction.     

Eco‐evolutionary feedback promotes Red Queen dynamics and selects for sex in predator populations

Although numerous hypotheses exist to explain the overwhelming presence of sexual reproduction across the tree of life, we still cannot explain its prevalence when considering all inherent costs involved. The Red Queen hypothesis states that sex is maintained because it can create novel genotypes with a selective advantage. This occurs when the interactions between species induce frequent environmental change. Here, we investigate whether coevolution and eco‐evolutionary feedback dynamics in a predator‐prey system allows for indirect selection and maintenance of sexual reproduction in the predator. Combining models and chemostat experiments of a rotifer‐algae system we show a continuous feedback between population and trait change along with recurrent shifts from selection by predation and competition for a limited resource. We found that a high propensity for sex was indirectly selected and was maintained in rotifer populations within environments containing these eco‐evolutionary dynamics; whereas within environments under constant conditions, predators evolved rapidly to lower levels of sex. Thus, our results indicate that the influence of eco‐evolutionary feedback dynamics on the overall evolutionary change has been underestimated.     

Sex and space destabilize intransitive competition within and between species

Organisms ranging from bacteria and corals to plants and vertebrates can form intransitive competitive networks, in which coexistence can be maintained because no one species or genotype is superior to all others. However, in the simplest case with three competing types, the long-term outcome may not be so clear if two of the three represent the ends of a continuous heritable trait distribution within one species, as has been recently demonstrated empirically in a short-term experiment with plants. Using simulation models of this scenario, results with asexual reproduction confirm previous studies which showed that local interactions promote coexistence. However, with sexual reproduction, genetic variance is reduced because selection fluctuates between favouring the two extremes during population cycles, while sex continually produces intermediates. Sex thus slows the response to selection when it is the strongest and therefore slows the recovery from extreme abundances, creating larger abundance fluctuations. Local interactions do not stabilize dynamics with sex because the resultant spatial patches of one species are genetically heterogeneous, such that particular phenotypes do not benefit from spatial refuges. In sharp contrast to previous models suggesting that sex or local interactions stabilize population dynamics, here sex and local interactions destabilize dynamics and increase extinction risk.     

Population and evolutionary dynamics in spatially structured seasonally varying environments

Increasingly imperative objectives in ecology are to understand and forecast population dynamic and evolutionary responses to seasonal environmental variation and change. Such population and evolutionary dynamics result from immediate and lagged responses of all key life‐history traits, and resulting demographic rates that affect population growth rate, to seasonal environmental conditions and population density. However, existing population dynamic and eco‐evolutionary theory and models have not yet fully encompassed within‐individual and among‐individual variation, covariation, structure and heterogeneity, and ongoing evolution, in a critical life‐history trait that allows individuals to respond to seasonal environmental conditions: seasonal migration. Meanwhile, empirical studies aided by new animal‐tracking technologies are increasingly demonstrating substantial within‐population variation in the occurrence and form of migration versus year‐round residence, generating diverse forms of ‘partial migration’ spanning diverse species, habitats and spatial scales. Such partially migratory systems form a continuum between the extreme scenarios of full migration and full year‐round residence, and are commonplace in nature. Here, we first review basic scenarios of partial migration and associated models designed to identify conditions that facilitate the maintenance of migratory polymorphism. We highlight that such models have been fundamental to the development of partial migration theory, but are spatially and demographically simplistic compared to the rich bodies of population dynamic theory and models that consider spatially structured populations with dispersal but no migration, or consider populations experiencing strong seasonality and full obligate migration. Second, to provide an overarching conceptual framework for spatio‐temporal population dynamics, we define a ‘partially migratory meta‐population’ system as a spatially structured set of locations that can be occupied by different sets of resident and migrant individuals in different seasons, and where locations that can support reproduction can also be linked by dispersal. We outline key forms of within‐individual and among‐individual variation and structure in migration that could arise within such systems and interact with variation in individual survival, reproduction and dispersal to create complex population dynamics and evolutionary responses across locations, seasons, years and generations. Third, we review approaches by which population dynamic and eco‐evolutionary models could be developed to test hypotheses regarding the dynamics and persistence of partially migratory meta‐populations given diverse forms of seasonal environmental variation and change, and to forecast system‐specific dynamics. To demonstrate one such approach, we use an evolutionary individual‐based model to illustrate that multiple forms of partial migration can readily co‐exist in a simple spatially structured landscape. Finally, we summarise recent empirical studies that demonstrate key components of demographic structure in partial migration, and demonstrate diverse associations with reproduction and survival. We thereby identify key theoretical and empirical knowledge gaps that remain, and consider multiple complementary approaches by which these gaps can be filled in order to elucidate population dynamic and eco‐evolutionary responses to spatio‐temporal seasonal environmental variation and change.     

Stage-structured ontogeny in resource populations generates non-additive stabilizing and de-stabilizing forces in populations and communities

Demographic heterogeneity influences how populations respond to density dependent intraspecific competition and trophic interactions. Distinct stages across an organism's development, or ontogeny, are an important example of demographic heterogeneity. In consumer populations, ontogenetic stage structure has been shown to produce categorical differences in population dynamics, community dynamics and even species coexistence compared to models lacking explicit ontogeny. The study of consumer–resource interactions must also consider the ontogenetic stage structure of the resource itself, particularly plants, given their fundamental role at the basis of terrestrial food webs. We incorporate distinct ontogenetic stages of plants into an adaptable multi-stage consumer–resource modeling framework that facilitates studying how stage specific consumers shape trophic dynamics at low trophic levels. We describe the role of density dependent demographic rates in mediating the dynamics of stage-structured plant populations. We then investigate how these demographic rates interact with consumer pressure to influence stability and coexistence in multiple stage-specific consumer–resource interactions. Results detail how density dependent effects across distinct ontogenetic stages in plant development produce non-additivity in the drivers of dynamic stability both in single populations and in consumer–resource settings, challenging the ubiquity of certain traditional ecological dynamic paradigms. We also find categorical differences in the population variability induced by herbivores consuming separate plant stages. Consumer–resource models, such as plant–herbivore interactions, often average out demographic heterogeneity in populations. Here, we show that explicitly including plant demographic heterogeneity through ontogeny yields distinct dynamic expectations for both plants and herbivores compared to traditional consumer–resource formulations. Our results indicate that efforts to understand the demographic effect of herbivores on plant populations may need to also consider the effects of plant demographics on herbivores and the reciprocal relationship between them.     

Soil-related variations in the population dynamics of six dipterocarp tree species with strong habitat preferences

Differences in the density of conspecific tree individuals in response to environmental gradients are well documented for many tree species, but how such density differences are generated and maintained is poorly understood. We examined the segregation of six dipterocarp species among three soil types in the Pasoh tropical forest, Malaysia. We examined how individual performance and population dynamics changed across the soil types using 10-year demographic data to compare tree performance across soil types, and constructed population matrix models to analyze the population dynamics. Species showed only minor changes in mortality and juvenile growth across soil types, although recruitment differed greatly. Clear, interspecific demographic trade-offs between growth and mortality were found in all soil types. The relative trade-offs by a species did not differ substantially among the soil types. Population sizes were projected to remain stable in all soil types for all species with one exception. Our life-table response experiment demonstrated that the population dynamics of a species differed only subtly among soil types. Therefore, species with strong density differences across soil types do not necessarily differ greatly in their population dynamics across the soil types. In contrast, interspecific differences in population dynamics were large. The trade-off between mortality and growth led to a negative correlation between the contributions of mortality and growth to variations in the population growth rate (λ) and thus reduced their net contributions. Recruitment had little impact on the variation in λ. The combination of these factors resulted in little variation in λ among species.     

Population dynamics under parasitic sex ratio distortion.

We analyse the population dynamic effects of sex ratio distortion by vertically transmitted, feminizing parasites. We show that, for diploid hosts, sex ratio distortion may lead to extinction as males become too rare to maintain the host population through reproduction. Feminizers can magnify Allee effects, broadening the range of conditions leading to extinction of small populations. Depending on male mating constraints and strength of density dependence, feminizers may either increase or decrease the equilibrium host density. Under conditions leading to deterministic host extinction, stochastic elimination of the parasite may allow the host population to recover. Hence, infection by parasitic sex ratio distorters may be transient in finite populations. We consider the implications of this process for parasite prevalence, host population regulation, and sex ratio evolution.     

Effect of reproductive modes and environmental heterogeneity in the population dynamics of a geographically widespread clonal desert cactus

The dynamics of plant populations in arid environments are largely affected by the unpredictable environmental conditions and are fine-tuned by biotic factors, such as modes of recruitment. A single species must cope with both spatial and temporal heterogeneity that trigger pulses of sexual and clonal establishment throughout its distributional range. We studied two populations of the clonal, purple prickly pear cactus, Opuntia macrocentra, in order to contrast the factors responsible for the population dynamics of a common, widely distributed species. The study sites were located in protected areas that correspond to extreme latitudinal locations for this species within the Chihuahuan Desert. We studied both populations for four consecutive years and determined the demographic consequences of environmental variability and the mode of reproduction using matrix population models, life table response experiments (LTREs), and loop and perturbation analyses. Although both populations seemed fairly stable (population growth rate, λ∼1), different demographic parameters and different life cycle routes were responsible for this stability in each population. In the southernmost population (MBR) LTRE and loop and elasticity analyses showed that stasis is the demographic process with the highest contributions to λ, followed by sexual reproduction, and clonal propagation contributed the least. The northern population (CR) had both higher elasticities and larger contributions of stasis, followed by clonal propagation and sexual recruitment. Loop analysis also showed that individuals in CR have more paths to complete a life cycle than those in MBR. As a consequence, each population differed in life history traits (e.g., size class structure, size at sexual maturity, and reproductive value). Numerical perturbation analyses showed a small effect of the seed bank on the λ of both populations, while the transition from seeds to seedlings had an important effect mainly in the northern population. Clonal propagation (higher survival and higher contributions to vital rates) seems to be more important for maintaining populations over long time periods than sexual reproduction.     

Reproduction and population density affect humoral immunity in bank voles under field experimental conditions

Density dependence is a common feature in the dynamics of animal populations. Availability of food resources critical to immunity is likely to be one of the mechanisms mediating the effect of population density on individual fitness. The ability to mount an immune response to an antigen is also affected by levels of immunosuppressive hormones associated with reproduction or mediating the response to ecological and social stress. We assessed variation in condition and intensity of humoral immune response to a T-cell-dependent antigen in bank voles (Clethrionomys glareolus) by experimentally altering population density before immunisation. Consistent with our prediction, males had lower humoral immunocompetence in the breeding than in the non-breeding season. Contrary to our expectation, males did not show enhanced immunocompetence and females showed depressed humoral immune response under experimentally lowered population density. Variation of immune response in relation to population density depended on sex, with females but not males showing lower immune response under experimentally reduced density. We conclude that humoral immunity of bank voles was affected by reproduction and social environment rather than by population density. Received: 2 November 1999 / Accepted: 22 March 2000