Search

COVID-19 brought back to the attention of the scientific community that males are more susceptible to infectious diseases. What is clear for other infections-that sex and gender differences influence both risk of infection and mortality-is not yet fully elucidated for plague, particularly bubonic plague, although this knowledge can help find specific defences against a disease for which a vaccine is not yet available. To address this question, we analysed data on plague from hospitals in different parts of the world since the early eighteenth century, which provide demographic information on individual patients, diagnosis and course of the disease in the pre-antibiotic era. Assuming that the two sexes were equally represented, we observe a worldwide prevalence of male cases hospitalized at any age, a result which seems better explained by gender-biased (thus cultural) behaviours than biological sex-related factors. Conversely, case fatality rates differ among countries and geographic macro-areas, while globally, lethality appears slightly prevalent in young females and older adults (regardless of sex). Logistic regression models confirm that the main risk factor for bubonic plague death was the geographical location of the cases and being older than 50 years, whereas sex only showcased a slight trend.

Biological processes and physical oceanography are often integrated in numerical modelling of marine fish larvae, but rarely in statistical analyses of spatio-temporal observation data. Here, we examine the relative contribution of inter-annual variability in spawner distribution, advection by ocean currents, hydrography and climate in modifying observed distribution patterns of cod larvae in the Lofoten-Barents Sea. By integrating predictions from a particle-tracking model into a spatially explicit statistical analysis, the effects of advection and the timing and locations of spawning are accounted for. The analysis also includes other environmental factors: temperature, salinity, a convergence index and a climate threshold determined by the North Atlantic Oscillation (NAO). We found that the spatial pattern of larvae changed over the two climate periods, being more upstream in low NAO years. We also demonstrate that spawning distribution and ocean circulation are the main factors shaping this distribution, while temperature effects are different between climate periods, probably due to a different spatial overlap of the fish larvae and their prey, and the consequent effect on the spatial pattern of larval survival. Our new methodological approach combines numerical and statistical modelling to draw robust inferences from observed distributions and will be of general interest for studies of many marine fish species.

Social parasites may exploit their hosts by mimicking other organisms that the hosts normally benefit from investing in or responding to in some other way. Some parasites exaggerate key characters of the organisms they mimic, possibly in order to increase the response from the hosts. The huge gape and extreme begging intensity of the parasitic common cuckoo chick (Cuculus canorus) may be an example. In this paper, the evolutionary stability of manipulating hosts through exaggerated signals is analysed using game theory. Our model indicates that a parasite's signal intensity must be below a certain threshold in order to ensure acceptance and that this threshold depends directly on the rate of parasitism. The only evolutionarily stable strategy (ESS) combination is when hosts accept all signallers and parasites signal at their optimal signal intensity, which must be below the threshold. Supernormal manipulation by parasites is only evolutionarily stable under sufficiently low rates of parasitism. If the conditions for the ESS combination are not satisfied, rejector hosts can invade using signal intensity as a cue for identifying parasites. These qualitative predictions are discussed with respect to empirical evidence from parasitic mimicry systems that have been suggested to involve supernormal signalling, including evicting avian brood parasites and insect-mimicking Ophrys orchids.

Here we present, to the authors' knowledge for the very first time for a small marsupial, a thorough analysis of the demography and population dynamics of the mouse opossum (Thylamys elegans) in western South America. We test the relative importance of feedback structure and climatic factors (rainfall and the Southern Oscillation Index) in explaining the temporal variation in the demography of the mouse opossum. The demographic information was incorporated into a stage-structured population dynamics model and the model's predictions were compared with observed patterns. The mouse opossum's capture rates showed seasonal (within-year) and between-year variability, with individuals having higher capture rates during late summer and autumn and lower capture rates during winter and spring. There was also a strong between-year effect on capture probabilities. The reproductive (the fraction of reproductively active individuals) and recruitment rates showed a clear seasonal and a between-year pattern of variation with the peak of reproductive activity occuring during winter and early spring. In addition, the fraction of reproductive individuals was positively related to annual rainfall, while population density and annual rainfall positively influenced the recruitment rate. The survival rates were negatively related to annual rainfall. The average finite population growth rate during the study period was estimated to be 1.011 +/- 0.0019 from capture-recapture estimates. While the annual growth rate estimated from the seasonal linear matrix models was 1.026, the subadult and adult survival and maturation rates represent between 54% (winter) and 81% (summer) of the impact on the annual growth rate.

A new geographical gradient in the dynamics of small rodents is demonstrated by analysing 29 time series of density indices of the common vole (Microtus arvalis) from Poland, the Czech Republic and the Slovak Republic. This gradient extends from more stable northerly populations in coastal Poland to more variable and cyclic populations in the southernmost parts of the Slovak Republic, and is hence a reversal of the Fennoscandian gradient. All studied variables (such as mean density, cycle amplitude, density variability and the coefficients in a second-order autoregressive model) exhibit consistent latitudinal variation. Possible underlying factors are discussed. In particular, we suggest that seasonality may be a key element in explaining the observed new gradient.

There are only a few recent studies that have demonstrated senescence in ungulates and nothing is known regarding how patterns of senescence may vary as a function of density Senescence in general is linked to the cost of reproduction, which probably increases with density in ungulates and may differ between the sexes. Further, senescence in ungulates is also regarded to be a function of tooth wear rates. As density dependence and sexual differences in food choice have been well documented, this may lead to different tooth wear rates and, thus, possibly density-dependent and sex-specific patterns of senescence. We therefore investigated the effects of age, sex, density and their possible interactions on the variability of body weight in 29,047 red deer harvested during 1965-1998 from Norway, out of which 380 males and 1452 females were eight years or older. There was clear evidence that spatio-temporal variation in density correlated negatively with body weights. In addition, there was evidence of senescence in both male and female red deer. Age at onset of senescence in females was after 20 years of age and independent of population density. In males, the onset and rate of senescence increased with increasing population density. The onset of senescence for males was at ca. 12 years of age at low density, but decreased to approximately ten years of age at high density. The pattern of density-dependent senescence in males, but not that in females, can be explained if (i) the cost of reproduction increases with density more strongly in male than in female red deer, and/or (ii) tooth wear rates are density dependent in males, but not in females. We discuss the ability of these two different, not mutually exclusive hypotheses in explaining the observed pattern of senescence.

In order to evaluate the biogeographical hypothesis that the Norwegian lemming (Lemmus lemmus) survived the last glacial period in some Scandinavian refugia, we examined variation in the nucleotide sequence of the mitochondrial control region (402 base pairs (bp)) and the cytochrome b (cyt b) region (633 bp) in Norwegian and Siberian (Lemmus sibiricus) lemmings. The phylogenetic distinction and cyt b divergence estimate of 1.8% between the Norwegian and Siberian lemmings suggest that their separation pre-dated the last glaciation and imply that the Norwegian lemming is probably a relic of the Pleistocene populations from Western Europe. The star-like control region phylogeny and low mitochondrial DNA diversity in the Norwegian lemming indicate a reduction in its historical effective size followed by population expansion. The average estimate of post-bottleneck time (19-21 kyr) is close to the last glacial maximum (18-22 kyr BP). Taking these findings and the fossil records into consideration, it seems likely that, after colonization of Scandinavia in the Late Pleistocene, the Norwegian lemming suffered a reduction in its population effective size and survived the last glacial maximum in some local Scandinavian refugia, as suggested by early biogeographical work.

Global climate fluctuated considerably throughout the Pliocene-Pleistocene period, influencing the evolutionary history of a wide array of species. Using the phylogeographic patterns within the hartebeest (Alcelaphus buselaphus (Pallas, 1766)) complex, we evaluated the evolutionary consequences of such environmental change for a typical large mammal ranging on the African savannah. Our results, as generated from two mitochondrial DNA markers (the D-loop and cytochrome b), suggest an origin of the hartebeest in eastern Africa from where the species has colonized other parts of the continent. Phylogenetic analyses revealed an early diversification into southern and northern hartebeest lineages, an event that may be related to the formation of the Rift Valley lakes. The northern lineage has further diverged into eastern and western lineages, most probably as a result of the expanding central African rainforest belt and subsequent contraction of savannah habitats during a period of global warming. The diversification events appear to have coincided with major climatic changes and are highly correlated in time. These observations strongly suggest that large-scale climatic fluctuations have been a major determinant for the species' evolutionary history and that hartebeest evolution has mainly taken place in isolated yet environmentally favourable refugia during periods of global warming. Indications of sudden population expansion for two putative ancestral hartebeest populations provide further support for a refugia-based explanation of the diversification events. Reciprocal monophyly between southern and northern lineages may suggest that reproductive barriers exist and that the hartebeest complex comprises two different species.

Although empirical and theoretical studies suggest that climate influences the timing of life-history events in animals and plants, correlations between climate and the timing of events such as egg-laying, migration or flowering do not reveal the mechanisms by which natural selection operates on life-history events. We present a general autoregressive model of the timing of life-history events in relation to variation in global climate that, like autoregressive models of population dynamics, allows for a more mechanistic understanding of the roles of climate, resources and competition. We applied the model to data on 50 years of annual dates of first flowering by three species of plants in 26 populations covering 4 degrees of latitude in Norway. In agreement with earlier studies, plants in most populations and all three species bloomed earlier following warmer winters. Moreover, our model revealed that earlier blooming reflected increasing influences of resources and density-dependent population limitation under climatic warming. The insights available from the application of this model to phenological data in other taxa will contribute to our understanding of the roles of endogenous versus exogenous processes in the evolution of the timing of life-history events in a changing climate.

Typically, in many studies in ecology, epidemiology, biomedicine and others, we are confronted with panels of short time-series of which we are interested in obtaining a biologically meaningful grouping. Here, we propose a bootstrap approach to test whether the regression functions or the variances of the error terms in a family of stochastic regression models are the same. Our general setting includes panels of time-series models as a special case. We rigorously justify the use of the test by investigating its asymptotic properties, both theoretically and through simulations. The latter confirm that for finite sample size, bootstrap provides a better approximation than classical asymptotic theory. We then apply the proposed tests to the mink-muskrat data across 81 trapping regions in Canada. Ecologically interpretable groupings are obtained, which serve as a necessary first step before a fuller biological and statistical analysis of the food chain interaction.

We analysed two 26-year long (1970-1995) time-series on annual population growth rates of Norwegian lemmings (Lemmus lemmus) from Finse, south Norway, using a threshold autoregressive (TAR) approach. We demonstrate that the population dynamics is both phase- and density-dependent. The phase-dependence accounts for the observed nonlinearity. We used the deduced stochastic model structure as a basis for evaluating the dynamic properties of this system. The dynamics is characterized either by limit cycles or chaos (the latter with a strong semi-periodic component). Stochasticity is seen to play an important role in the determination of the periodicity. The ecological implications of these statistical and mathematical results are discussed.

We present a model on plant-deer climate interactions developed for improving our understanding of the temporal dynamics of deer abundance and, in particular, how intrinsic (density-dependent) and extrinsic (plants, climate) factors influence these dynamics. The model was tested statistically by analysing the dynamics of five Norwegian red deer populations between 1964 and 1993. Direct and delayed density-dependence significantly influenced the development of the populations: delayed density-dependence primarily operated through female density, whereas direct density-dependence acted through both female and male densities. Furthermore, population dynamics of Norwegian red deer were significantly affected by climate (as measured by the global weather phenomenon, the North Atlantic Oscillation: NAO). Warm, snowy winters (high NAO) were associated with decreased deer abundance, whereas the delayed (two-year) effect of warm, snowy winters had a positive effect on deer abundance. Our analyses are argued to have profound implications for the general understanding of climate change and terrestrial ecosystem functioning.

The variability of two fitness-related phenotypic traits (body weight and a mandibular skeletal ratio) was analysed among cohorts and age-classes of red deer in Norway. Phenotypic variation among cohorts was pronounced for calves, yearlings and reproductively mature adults. Fluctuations in cohort-specific mean body weights and skeletal ratios of adults correlated with global climatic variation in winter conditions influenced by the North Atlantic Oscillation while cohorts were in utero. Red deer born following warm winters were smaller than those born after cold winters, and this inter-cohort variability persisted into adulthood. Phenotypic variation among cohorts of red deer influenced by climate change may pose consequences for fitness of cohorts since body size and condition contribute to reproductive success and survival in male and female red deer. In particular, the recent trend of increasingly warm winters in northern Europe and Scandinavia may lead to reduced body size and fecundity of red deer, and perhaps other ungulates, in those areas.

Current ecological information on periodically fluctuating microtine populations are demonstrated to support a hypothesis involving both predation and intrinsic self-regulation as necessary and sufficient factors for explaining the "microtine density cycle'. The structure of the cyclic time series is largely two dimensional with strong delayed density dependence. Together with recent field studies on rodent demography, our modelling suggests that trophic interaction is a likely candidate to generate the dimensionality observed for northern microtine rodent dynamics. It is shown that the trophic interaction must be fairly strong. This suggests that specialist predation is the most likely one among the classes of trophic interactions. We also argue that some - but not too strong - self-regulation must occur to generate the structure of the available time series on northern European microtines.

A total of 31 years of abundance data from 90 populations of the grey-sided vole (Clethrionomys rufocanus) in northern Hokkaido (Japan) were analysed with respect to population dynamic characteristics. Both non-periodic and multi-annually periodic fluctuations occur among the studied populations. The length of the period varies from 2 to 5 years. The dynamics appear approximately linear on a logarithmic scale, and a log-linear stochastic difference model with one time-lag is found to recreate the periodograms of the time series. The deterministic (log-linear) component of these models determines the periodicity while stochasticity sustains the cycles. The dynamics of the Hokkaidian vole populations are found to vary clinally from the western coast eastwards and towards the interior of the study area. This gradient corresponds superficially to the latitudinal gradient seen in microtine populations in Fennoscandia. However, under close scrutiny, these gradients differ greatly: the Hokkaidian gradient is caused by a cline in delayed density dependence. Statistical delayed density dependence is more negative towards the east and the interior.

The patterns of density dependence in Fennoscandian rodents are investigated statistically using a linear autoregressive scheme. Nineteen time series of microtine abundances along a latitudinal gradient in Fennoscandia from 60 degrees N to 69 degrees N are analysed. We provide statistical evidence that there exists a latitudinal gradient in density dependence in Fennoscandian microtines. Southern populations experience significantly stronger direct density dependence than northern populations. Delayed density dependence was significantly negative throughout the region and appeared constant across the latitudinal gradient. The populations consistently exhibit dynamics of second order throughout the region. Together, the clinal direct density dependence and constant delayed density dependence give rise to a cline in cycle period from 3 to 4.5 years. The statistical results are compared to assumptions and predictions made in previous studies on the geographic gradient in the population dynamics of these rodents. The results are in agreement with the predictions of the 'generalist predator hypothesis'.

The dominant Lyapunov exponent, as estimated from time series using the Jacobian-based method, is often used for indicating whether the underlying dynamic system is chaotic or not. The Jacobian-based method together with Response Surface Methodology has been suggested as a method for detecting chaotic dynamics in ecological time series. Besides pointing out that this may not be an appropriate method, we report on estimates of the uncertainty in the estimates of the dominant Lyapunov exponent. For this purpose, we have used time series data on Holarctic microtines. On the basis of our analyses, we are unable to find general evidence for chaotic dynamics in northern microtine populations (north of 60 degrees N) as recently suggested in the ecological literature. It seems, however, that the dynamic properties of the northern and southern populations are different. These patterns are supported by testing for nonlinearity.