Taking the beat of the Arctic: are lemming population cycles changing due to winter climate?

Reports of fading vole and lemming population cycles and persisting low populations in some parts of the Arctic have raised concerns about the spread of these fundamental changes to tundra food web dynamics. By compiling 24 unique time series of lemming population fluctuations across the circumpolar region, we show that virtually all populations displayed alternating periods of cyclic/non-cyclic fluctuations over the past four decades. Cyclic patterns were detected 55% of the time (n = 649 years pooled across sites) with a median periodicity of 3.7 years, and non-cyclic periods were not more frequent in recent years. Overall, there was an indication for a negative effect of warm spells occurring during the snow onset period of the preceding year on lemming abundance. However, winter duration or early winter climatic conditions did not differ on average between cyclic and non-cyclic periods. Analysis of the time series shows that there is presently no Arctic-wide collapse of lemming cycles, even though cycles have been sporadic at most sites during the last decades. Although non-stationary dynamics appears a common feature of lemming populations also in the past, continued warming in early winter may decrease the frequency of periodic irruptions with negative consequences for tundra ecosystems.

Exploring the influence of host community composition on the outbreak potential of Anaplasma phagocytophilum and Borrelia burgdorferi s.l

In large parts of the northern hemisphere, multiple deer species coexist, and management actions can strongly influence wild deer communities. Such changes may also indirectly influence other species in the community, such as small mammals and birds, because deer can have strong effects on their habitats and resources. Deer, small mammals and birds play an important role in the dynamics of tick-borne zoonotic diseases. It is, however, relatively underexplored how the abundance and composition of vertebrate communities may affect the outbreak potential, maintenance and circulation of tick-borne pathogens. In this study we focus on the outbreak potential by exploring how the basic reproduction number R0 for different tick-borne pathogens depends on host community composition. We used published data on co-varying roe deer (Capreolus capreolus) and fallow deer (Dama dama) densities following a hunting ban, and different small mammal and bird densities, to investigate how the change in host community influences the R0 of four tick-borne pathogens: one non-zoonotic, namely Anaplasma phagocytophilum ecotype 2, and three zoonotic, namely A. phagocytophilum ecotype 1, Borrelia afzelii and Borrelia garinii. We calculated R0 using a next generation matrix approach, and used elasticities to quantify the contributions to R0 of the different groups of host species. The value of R0 for A. phagocytophilum ecotype 1 was higher with high fallow deer density and low roe deer density, while it was the other way round for A. phagocytophilum ecotype 2. For B. afzelii, R0 was mostly related to the density of small mammals and for B. garinii it was mostly determined by bird density. Our results show that the effect of species composition is substantial in the outbreak potential of tick-borne pathogens. This implies that also management actions that change this composition, can (indirectly and unintentionally) affect the outbreak potential of tick-borne diseases.

Urban forest soils harbour distinct and more diverse communities of bacteria and fungi compared to less disturbed forest soils

Anthropogenic changes to land use drive concomitant changes in biodiversity, including that of the soil microbiota. However, it is not clear how increasing intensity of human disturbance is reflected in the soil microbial communities. To address this issue, we used amplicon sequencing to quantify the microbiota (bacteria and fungi) in the soil of forests (n = 312) experiencing four different land uses, national parks (set aside for nature conservation), managed (for forestry purposes), suburban (on the border of an urban area) and urban (fully within a town or city), which broadly represent a gradient of anthropogenic disturbance. Alpha diversity of bacteria and fungi increased with increasing levels of anthropogenic disturbance, and was thus highest in urban forest soils and lowest in the national parks. The forest soil microbial communities were structured according to the level of anthropogenic disturbance, with a clear urban signature evident in both bacteria and fungi. Despite notable differences in community composition, there was little change in the predicted functional traits of urban bacteria. By contrast, urban soils exhibited a marked loss of ectomycorrhizal fungi. Soil pH was positively correlated with the level of disturbance, and thus was the strongest predictor of variation in alpha and beta diversity of forest soil communities, indicating a role of soil alkalinity in structuring urban soil microbial communities. Hence, our study shows how the properties of urban forest soils promote an increase in microbial diversity and a change in forest soil microbiota composition.

Population fluctuations and synanthropy explain transmission risk in rodent-borne zoonoses

Population fluctuations are widespread across the animal kingdom, especially in the order Rodentia, which includes many globally important reservoir species for zoonotic pathogens. The implications of these fluctuations for zoonotic spillover remain poorly understood. Here, we report a global empirical analysis of data describing the linkages between habitat use, population fluctuations and zoonotic reservoir status in rodents. Our quantitative synthesis is based on data collated from papers and databases. We show that the magnitude of population fluctuations combined with species' synanthropy and degree of human exploitation together distinguish most rodent reservoirs at a global scale, a result that was consistent across all pathogen types and pathogen transmission modes. Our spatial analyses identified hotspots of high transmission risk, including regions where reservoir species dominate the rodent community. Beyond rodents, these generalities inform our understanding of how natural and anthropogenic factors interact to increase the risk of zoonotic spillover in a rapidly changing world.

Puumala Orthohantavirus Infection Does Not Affect the Trapping Success of Its Reservoir Host

Pathogens might affect behavior of infected reservoir hosts and hence their trappability, which could bias population estimates of pathogen prevalence. In this study, we used snap-trapping data on Puumala orthohantavirus (PUUV)-infected (n = 1619) and noninfected (n = 6940) bank voles (Myodes glareolus) from five vole cycles, normally representing increase, peak, and decline phase, to evaluate if infection status affected trapping success. If PUUV infection, as previously suggested, increases activity and/or mobility, we would expect a higher proportion of infected than noninfected specimens in the first trapping night. However, the proportion of PUUV-infected voles did not differ across the three trapping nights. We conclude that PUUV infection did not affect trapping success, confirming snap trapping as an appropriate trapping method for studies on PUUV prevalence and likely other orthohantaviruses.

Alphacoronavirus in a Daubenton’s Myotis Bat (Myotis daubentonii) in Sweden

The ongoing COVID-19 pandemic has stimulated a search for reservoirs and species potentially involved in back and forth transmission. Studies have postulated bats as one of the key reservoirs of coronaviruses (CoVs), and different CoVs have been detected in bats. So far, CoVs have not been found in bats in Sweden and we therefore tested whether they carry CoVs. In summer 2020, we sampled a total of 77 adult bats comprising 74 Myotis daubentonii, 2 Pipistrellus pygmaeus, and 1 M. mystacinus bats in southern Sweden. Blood, saliva and feces were sampled, processed and subjected to a virus next-generation sequencing target enrichment protocol. An Alphacoronavirus was detected and sequenced from feces of a M. daubentonii adult female bat. Phylogenetic analysis of the almost complete virus genome revealed a close relationship with Finnish and Danish strains. This was the first finding of a CoV in bats in Sweden, and bats may play a role in the transmission cycle of CoVs in Sweden. Focused and targeted surveillance of CoVs in bats is warranted, with consideration of potential conflicts between public health and nature conservation required as many bat species in Europe are threatened and protected.

Climate change accelerates winter transmission of a zoonotic pathogen

Many zoonotic diseases are weather sensitive, raising concern how their distribution and outbreaks will be affected by climate change. At northern high latitudes, the effect of global warming on especially winter conditions is strong. By using long term monitoring data (1980-1986 and 2003-2013) from Northern Europe on temperature, precipitation, an endemic zoonotic pathogen (Puumala orthohantavirus, PUUV) and its reservoir host (the bank vole, Myodes glareolus), we show that early winters have become increasingly wet, with a knock-on effect on pathogen transmission in its reservoir host population. Further, our study is the first to show a climate change effect on an endemic northern zoonosis, that is not induced by increased host abundance or distribution, demonstrating that climate change can also alter transmission intensity within host populations. Our results suggest that rainy early winters accelerate PUUV transmission in bank voles in winter, likely increasing the human zoonotic risk in the North.

Low-level environmental metal pollution is associated with altered gut microbiota of a wild rodent, the bank vole (Myodes glareolus)

Mining and related industries are a major source of metal pollution. In contrast to the well-studied effects of exposure to metals on animal physiology and health, the impacts of environmental metal pollution on the gut microbiota of wild animals are virtually unknown. As the gut microbiota is a key component of host health, it is important to understand whether metal pollution can alter wild animal gut microbiota composition. Using a combination of 16S rRNA amplicon sequencing and quantification of metal levels in kidneys, we assessed whether multi-metal exposure (the sum of normalized levels of fifteen metals) was associated with changes in gut microbiota of wild bank voles (Myodes glareolus) from two locations in Finland. Exposure to increased metal load was associated with higher gut microbiota species diversity (α-diversity) and altered community composition (β-diversity), but not dispersion. Multi-metal exposure and increased levels of several metals (Cd, Hg, Pb and Se) were associated with differences in the abundance of microbial taxa, especially those within the families Clostridiales vadinBB60 group, Desulfovibrionaceae, Lachnospiraceae, Muribaculaceae and Ruminococcaceae. Our data indicate that even low-level metal pollution can affect the diversity of microbiota and be associated with deterministic differences in composition of host gut microbiota in wild animal populations. These findings highlight the need to study a broader range of metals and their cocktails that are more representative of the types of environmental exposure experienced by wild animals.

Wild ungulate species differ in their contribution to the transmission of Ixodes ricinus-borne pathogens

BACKGROUND

Several ungulate species are feeding and propagation hosts for the tick Ixodes ricinus as well as hosts to a wide range of zoonotic pathogens. Here, we focus on Anaplasma phagocytophilum and Borrelia burgdorferi (s.l.), two important pathogens for which ungulates are amplifying and dilution hosts, respectively. Ungulate management is one of the main tools to mitigate human health risks associated with these tick-borne pathogens. Across Europe, different species of ungulates are expanding their ranges and increasing in numbers. It is currently unclear if and how the relative contribution to the life-cycle of I. ricinus and the transmission cycles of tick-borne pathogens differ among these species. In this study, we aimed to identify these relative contributions for five European ungulate species.

METHODS

We quantified the tick load and collected ticks and spleen samples from hunted fallow deer (Dama dama, n = 131), moose (Alces alces, n = 15), red deer (Cervus elaphus, n = 61), roe deer (Capreolus capreolus, n = 30) and wild boar (Sus scrofa, n = 87) in south-central Sweden. We investigated the presence of tick-borne pathogens in ticks and spleen samples using real-time PCR. We determined if ungulate species differed in tick load (prevalence and intensity) and in infection prevalence in their tissue as well as in the ticks feeding on them.

RESULTS

Wild boar hosted fewer adult female ticks than any of the deer species, indicating that deer are more important as propagation hosts. Among the deer species, moose had the lowest number of female ticks, while there was no difference among the other deer species. Given the low number of infected nymphs, the relative contribution of all ungulate species to the transmission of B. burgdorferi (s.l.) was low. Fallow deer, red deer and roe deer contributed more to the transmission of A. phagocytophilum than wild boar.

CONCLUSIONS

The ungulate species clearly differed in their role as a propagation host and in the transmission of B. burgdorferi and A. phagocytophilum. This study provides crucial information for ungulate management as a tool to mitigate zoonotic disease risk and argues for adapting management approaches to the local ungulate species composition and the pathogen(s) of concern.

Geographical Distribution and Genetic Diversity of Bank Vole Hepaciviruses in Europe

The development of new diagnostic methods resulted in the discovery of novel hepaciviruses in wild populations of the bank vole (Myodes glareolus, syn. Clethrionomys glareolus). The naturally infected voles demonstrate signs of hepatitis similar to those induced by hepatitis C virus (HCV) in humans. The aim of the present research was to investigate the geographical distribution of bank vole-associated hepaciviruses (BvHVs) and their genetic diversity in Europe. Real-time reverse transcription polymerase chain reaction (RT-qPCR) screening revealed BvHV RNA in 442 out of 1838 (24.0%) bank voles from nine European countries and in one of seven northern red-backed voles (Myodes rutilus, syn. Clethrionomys rutilus). BvHV RNA was not found in any other small mammal species (n = 23) tested here. Phylogenetic and isolation-by-distance analyses confirmed the occurrence of both BvHV species (Hepacivirus F and Hepacivirus J) and their sympatric occurrence at several trapping sites in two countries. The broad geographical distribution of BvHVs across Europe was associated with their presence in bank voles of different evolutionary lineages. The extensive geographical distribution and high levels of genetic diversity of BvHVs, as well as the high population fluctuations of bank voles and occasional commensalism in some parts of Europe warrant future studies on the zoonotic potential of BvHVs.

Prevalence of Orthohantavirus-Reactive Antibodies in Humans and Peri-Domestic Rodents in Northern Ethiopia

In 2012, Tigray orthohantavirus was discovered in Ethiopia, but its seasonal infection in small mammals, and whether it poses a risk to humans was unknown. The occurrence of small mammals, rodents and shrews, in human inhabitations in northern Ethiopia is affected by season and presence of stone bunds. We sampled small mammals in two seasons from low- and high-density stone bund fields adjacent to houses and community-protected semi-natural habitats in Atsbi and Hagere Selam, where Tigray orthohantavirus was first discovered. We collected blood samples from both small mammals and residents using filter paper. The presence of orthohantavirus-reactive antibodies in blood was then analyzed using immunofluorescence assay (human samples) and enzyme linked immunosorbent assays (small mammal samples) with Puumala orthohantavirus as antigen. Viral RNA was detected by RT-PCR using small mammal blood samples. Total orthohantavirus prevalence (antibodies or virus RNA) in the small mammals was 3.37%. The positive animals were three Stenocephalemys albipes rats (prevalence in this species = 13.04%). The low prevalence made it impossible to determine whether season and stone bunds were associated with orthohantavirus prevalence in the small mammals. In humans, we report the first detection of orthohantavirus-reactive IgG antibodies in Ethiopia (seroprevalence = 5.26%). S. albipes lives in close proximity to humans, likely increasing the risk of zoonotic transmission.

Global patterns and determinants of lake macrophyte taxonomic, functional and phylogenetic beta diversity

Documenting the patterns of biological diversity on Earth has always been a central challenge in macroecology and biogeography. However, for the diverse group of freshwater plants, such research program is still in its infancy. Here, we examined global variation in taxonomic, functional and phylogenetic beta diversity patterns of lake macrophytes using regional data from six continents. A data set of ca. 480 lake macrophyte community observations, together with climatic, geographical and environmental variables, was compiled across 16 regions worldwide. We (a) built the very first phylogeny comprising most freshwater plant lineages; (b) exploited a wide array of functional traits that are important to macrophyte autoecology or that relate to lake ecosystem functioning; (c) assessed if different large-scale beta diversity patterns show a clear latitudinal gradient from the equator to the poles using null models; and (d) employed evolutionary and regression models to first identify the degree to which the studied functional traits show a phylogenetic signal, and then to estimate community-environment relationships at multiple spatial scales. Our results supported the notion that ecological niches evolved independently of phylogeny in macrophyte lineages worldwide. We also showed that taxonomic and phylogenetic beta diversity followed the typical global trend with higher diversity in the tropics. In addition, we were able to confirm that species, multi-trait and lineage compositions were first and foremost structured by climatic conditions at relatively broad spatial scales. Perhaps more importantly, we showed that large-scale processes along latitudinal and elevational gradients have left a strong footprint in the current diversity patterns and community-environment relationships in lake macrophytes. Overall, our results stress the need for an integrative approach to macroecology, biogeography and conservation biology, combining multiple diversity facets at different spatial scales.

Geographical Distribution of Ljungan Virus in Small Mammals in Europe

Ljungan virus (LV), which belongs to the Parechovirus genus in the Picornaviridae family, was first isolated from bank voles (Myodes glareolus) in Sweden in 1998 and proposed as a zoonotic agent. To improve knowledge of the host association and geographical distribution of LV, tissues from 1685 animals belonging to multiple rodent and insectivore species from 12 European countries were screened for LV-RNA using reverse transcriptase (RT)-PCR. In addition, we investigated how the prevalence of LV-RNA in bank voles is associated with various intrinsic and extrinsic factors. We show that LV is widespread geographically, having been detected in at least one host species in nine European countries. Twelve out of 21 species screened were LV-RNA PCR positive, including, for the first time, the red vole (Myodes rutilus) and the root or tundra vole (Alexandromys formerly Microtus oeconomus), as well as in insectivores, including the bicolored white-toothed shrew (Crocidura leucodon) and the Valais shrew (Sorex antinorii). Results indicated that bank voles are the main rodent host for this virus (overall RT-PCR prevalence: 15.2%). Linear modeling of intrinsic and extrinsic factors that could impact LV prevalence showed a concave-down relationship between body mass and LV occurrence, so that subadults had the highest LV positivity, but LV in older animals was less prevalent. Also, LV prevalence was higher in autumn and lower in spring, and the amount of precipitation recorded during the 6 months preceding the trapping date was negatively correlated with the presence of the virus. Phylogenetic analysis on the 185 base pair species-specific sequence of the 5' untranslated region identified high genetic diversity (46.5%) between 80 haplotypes, although no geographical or host-specific patterns of diversity were detected.

Selective Predation by Owls on Infected Bank Voles (Myodes glareolus) as a Possible Sentinel of Tularemia Outbreaks

Tularemia is a widely spread zoonotic disease in the northern hemisphere, caused by the bacterium Francisella tularensis. In humans, tularemia is an acute febrile illness with incidence peaks in late summer to early autumn of outbreak years, but there is no early warning system in place that can reduce the impact of disease by providing timely risk information. In this study, we revisit previously unpublished data on F. tularensis in water, sediment, soil, and small mammals from 1984 in northern Sweden. In addition, we used human case data from the national surveillance system for tularemia in the same year. In the environmental and small mammal material, bank vole (Myodes glareolus) samples from urine and bladder were the only samples that tested positive for F. tularensis. The prevalence of F. tularensis among trapped bank voles was 13.5%, although all six bank voles that were retrieved from owl nest boxes in early May tested positive. Forty-two human tularemia cases were reported from August to December in 1984. Based on these results, we encourage investigating the potential role of tularemia-infected bank voles retrieved from owl nest boxes in spring as an early warning for outbreaks of tularemia among humans in summer and autumn of the same year.

Spatio-temporal variation of metals and organic contaminants in bank voles (Myodes glareolus)

Environmental contamination with metals and organic compounds is of increasing concern for ecosystem and human health. Still, our knowledge about spatial distribution, temporal changes and ecotoxicological fate of metals and organic contaminants in wildlife is limited. We studied concentrations of 69 elements and 50 organic compounds in 300 bank voles (Myodes glareolus), Europe's most common mammal, sampled in spring and autumn 2017-2018 in five monitoring areas, representing three biogeographic regions. In addition, we compared measured concentrations with previous results from bank voles sampled within the same areas in 1995-1997 and 2001. In general, our results show regional differences, but no consistent patterns among contaminants and study areas. The exception was for the lowest concentrations of organic contaminants (e.g. perfluorooctane sulfonate, PFOS), which were generally found in the northern Swedish mountain area. Concentrations of metals and organic contaminants in adults varied seasonally with most organic contaminants being higher in spring; likely induced by diet shifts but potentially also related to age differences. In addition, metal concentrations varied between organs (liver vs. kidney), age classes (juveniles vs. adults; generally higher in adults) as well as between males and females. Concentrations of chromium and nickel in kidney and liver in the northernmost mountain area were lower in 2017-2018 than in 1995-1997 and in three of four areas, lead concentrations were lower in 2017-2018 than in 2001. Current metal concentrations (except mercury) are not expected to negatively affect the voles. Concentrations of hexachlorobenzene displayed highest concentrations in 2001 in the mountains, while it was close to detection limit in 2017-2018. Likewise, PFOS concentrations decreased in the mountains and in south-central lowland forests between 2001 and 2017-2018. Our results suggest that season, age class and sex need to be considered when designing and interpreting results from monitoring programs targeting inorganic and organic contaminants in wildlife.

Documenting lemming population change in the Arctic: Can we detect trends?

Lemmings are a key component of tundra food webs and changes in their dynamics can affect the whole ecosystem. We present a comprehensive overview of lemming monitoring and research activities, and assess recent trends in lemming abundance across the circumpolar Arctic. Since 2000, lemmings have been monitored at 49 sites of which 38 are still active. The sites were not evenly distributed with notably Russia and high Arctic Canada underrepresented. Abundance was monitored at all sites, but methods and levels of precision varied greatly. Other important attributes such as health, genetic diversity and potential drivers of population change, were often not monitored. There was no evidence that lemming populations were decreasing in general, although a negative trend was detected for low arctic populations sympatric with voles. To keep the pace of arctic change, we recommend maintaining long-term programmes while harmonizing methods, improving spatial coverage and integrating an ecosystem perspective.

Correction to: Documenting lemming population change in the Arctic: Can we detect trends?

In the original published article, some of the symbols in figure 1A were modified incorrectly during the typesetting and publication process. The correct version of the figure is provided in this correction.

Effect of spatial scale and latitude on diversity-disease relationships

Natural ecosystems provide humans with different types of ecosystem services, often linked to biodiversity. The dilution effect (DE) predicts a negative relationship between biodiversity and risk of infectious diseases of humans, other animals, and plants. We hypothesized that a stronger DE would be observed in studies conducted at smaller spatial scales, where biotic drivers may predominate, compared to studies at larger spatial scales where abiotic drivers may more strongly affect disease patterns. In addition, we hypothesized a stronger DE in studies from temperate regions at mid latitudes than in those from subtropical and tropical regions, due to more diffuse species interactions at low latitudes. To explore these hypotheses, we conducted a meta‐analysis of observational studies of diversity–disease relationships for animals across spatial scales and geographic regions. Negative diversity–disease relationships were significant at small (combined site and local), intermediate (combined landscape and regional), and large (combined continental and global) scales and the effect did not differ depending on size of the study areas. For the geographic region analysis, a strongly negative diversity–disease relationship was found in the temperate region while no effect was found in the subtropical and tropical regions. However, no overall effect of absolute latitude on the strength of the dilution effect was detected. Our results suggest that a negative diversity–disease relationship occurs across scales and latitudes and is especially strong in the temperate region. These findings may help guide future management efforts in lowering disease risk.

Catchment properties and the photosynthetic trait composition of freshwater plant communities

Unlike in land plants, photosynthesis in many aquatic plants relies on bicarbonate in addition to carbon dioxide (CO2) to compensate for the low diffusivity and potential depletion of CO2 in water. Concentrations of bicarbonate and CO2 vary greatly with catchment geology. In this study, we investigate whether there is a link between these concentrations and the frequency of freshwater plants possessing the bicarbonate use trait. We show, globally, that the frequency of plant species with this trait increases with bicarbonate concentration. Regionally, however, the frequency of bicarbonate use is reduced at sites where the CO2 concentration is substantially above the air equilibrium, consistent with this trait being an adaptation to carbon limitation. Future anthropogenic changes of bicarbonate and CO2 concentrations may alter the species compositions of freshwater plant communities.

Wildfire-induced short-term changes in a small mammal community increase prevalence of a zoonotic pathogen?

Natural disturbances like droughts and fires are important determinants of wildlife community structure and are suggested to have important implications for prevalence of wildlife-borne pathogens. After a major wildfire affecting >1,600 ha of boreal forest in Sweden in 2006, we took the rare opportunity to study the short-term response (2007-2010 and 2015) of small mammal community structure, population dynamics, and prevalence of the Puumala orthohantavirus (PUUV) hosted by bank voles (Myodes glareolus). We performed snap-trapping in permanent trapping plots in clear-cuts (n = 3), unburnt reference forests (n = 7), and the fire area (n = 7) and surveyed vegetation and habitat structure. Small mammal species richness was low in all habitats (at maximum three species per trapping session), and the bank vole was the only small mammal species encountered in the fire area after the first postfire year. In autumns of years of peak rodent densities, the trapping index of bank voles was lowest in the fire area, and in two of three peak-density years, it was highest in clear-cuts. Age structure of bank voles varied among forest types with dominance of overwintered breeders in the fire area in the first postfire spring. PUUV infection probability in bank voles was positively related to vole age. Infection probability was highest in the fire area due to low habitat complexity in burnt forests, which possibly increased encounter rate among bank voles. Our results suggest that forest fires induce cascading effects, including fast recovery/recolonization of fire areas by generalists like bank voles, impoverished species richness of small mammals, and altered prevalence of a rodent-borne zoonotic pathogen. Our pilot study suggests high human infection risk upon encountering a bank vole in the fire area, however, with even higher overall risk in unburnt forests due to their higher vole numbers.

OPEN RESEARCH BADGES

This article has earned an Open Data Badge for making publicly available the digitally-shareable data necessary to reproduce the reported results. The data is available at https://osf.io/6fsy3/.

Population Dynamics of Bank Voles Predicts Human Puumala Hantavirus Risk

Predicting risk of zoonotic diseases, i.e., diseases shared by humans and animals, is often complicated by the population ecology of wildlife host(s). We here demonstrate how ecological knowledge of a disease system can be used for early prediction of human risk using Puumala hantavirus (PUUV) in bank voles (Myodes glareolus), which causes Nephropathia epidemica (NE) in humans, as a model system. Bank vole populations at northern latitudes exhibit multiannual fluctuations in density and spatial distribution, a phenomenon that has been studied extensively. Nevertheless, existing studies predict NE incidence only a few months before an outbreak. We used a time series on cyclic bank vole population density (1972-2013), their PUUV infection rates (1979-1986; 2003-2013), and NE incidence in Sweden (1990-2013). Depending on the relationship between vole density and infection prevalence (proportion of infected animals), either overall density of bank voles or the density of infected bank voles may be used to predict seasonal NE incidence. The density and spatial distribution of voles at density minima of a population cycle contribute to the early warning of NE risk later at its cyclic peak. When bank voles remain relatively widespread in the landscape during cyclic minima, PUUV can spread from a high baseline during a cycle, culminating in high prevalence in bank voles and potentially high NE risk during peak densities.

Diet shift in bank voles induced by competition from grey-sided voles?

Grey-sided voles (Myodes rufocanus) and bank voles (Myodes glareolus) co-exist in boreal forests in northern Scandinavia. Previous studies suggest that the 2 species interact interspecifically, the grey-sided vole being the dominant species. We tested the hypothesis that bank voles shift their diet due to competition with the dominant grey-sided vole by studying stable isotope ratios in both species. Muscle samples were taken from voles in patches of old forest occupied by only bank voles and patches of old forest occupied by both grey-sided voles and bank voles. We found that: (i) stable isotope ratios of bank voles differed in areas with and without grey-sided voles; and that (ii) the stable isotope ratios of bank voles were more similar to those of grey-sided voles in areas where grey-sided voles were absent. Our data suggests that grey-sided voles forced bank voles to change their diet due to interspecific competition.

Global patterns in the metacommunity structuring of lake macrophytes: regional variations and driving factors

We studied community-environment relationships of lake macrophytes at two metacommunity scales using data from 16 regions across the world. More specifically, we examined (a) whether the lake macrophyte communities respond similar to key local environmental factors, major climate variables and lake spatial locations in each of the regions (i.e., within-region approach) and (b) how well can explained variability in the community-environment relationships across multiple lake macrophyte metacommunities be accounted for by elevation range, spatial extent, latitude, longitude, and age of the oldest lake within each metacommunity (i.e., across-region approach). In the within-region approach, we employed partial redundancy analyses together with variation partitioning to investigate the relative importance of local variables, climate variables, and spatial location on lake macrophytes among the study regions. In the across-region approach, we used adjusted R2 values of the variation partitioning to model the community-environment relationships across multiple metacommunities using linear regression and commonality analysis. We found that niche filtering related to local lake-level environmental conditions was the dominant force structuring macrophytes within metacommunities. However, our results also revealed that elevation range associated with climate (increasing temperature amplitude affecting macrophytes) and spatial location (likely due to dispersal limitation) was important for macrophytes based on the findings of the across-metacommunities analysis. These findings suggest that different determinants influence macrophyte metacommunities within different regions, thus showing context dependency. Moreover, our study emphasized that the use of a single metacommunity scale gives incomplete information on the environmental features explaining variation in macrophyte communities.

Rewilding wetlands: beaver as agents of within-habitat heterogeneity and the responses of contrasting biota

Ecosystem engineers can increase biodiversity by creating novel habitat supporting species that would otherwise be absent. Their more routine activities further influence the biota occupying engineered habitats. Beavers are well-known for transforming ecosystems through dam building and are therefore increasingly being used for habitat restoration, adaptation to climate extremes and in long-term rewilding. Abandoned beaver ponds (BP) develop into meadows or forested wetlands that differ fundamentally from other terrestrial habitats and thus increase landscape diversity. Active BP, by contrast, are superficially similar to other non-engineered shallow wetlands, but ongoing use and maintenance might affect how BP contribute to aquatic biodiversity. We explored the 'within-habitat' effect of an ecosystem engineer by comparing active BP in southern Sweden with coexisting other wetlands (OW), using sedentary (plants) and mobile (water beetles) organisms as indicators. BP differed predictably from OW in environmental characteristics and were more heterogeneous. BP supported more plant species at plot (+15%) and site (+33%) scales, and plant beta diversity, based on turnover between plots, was 17% higher than in OW, contributing to a significantly larger species pool in BP (+17%). Beetles were not differentiated between BP and OW based on diversity measures but were 26% more abundant in BP. Independent of habitat creation beaver are thus significant agents of within-habitat heterogeneity that differentiates BP from other standing water habitat; as an integral component of the rewilding of wetlands re-establishing beaver should benefit aquatic biodiversity across multiple scales.This article is part of the theme issue 'Trophic rewilding: consequences for ecosystems under global change'.

Mixtures of macrophyte growth forms promote nitrogen cycling in wetlands

The importance of aquatic plant diversity in regulating nutrient cycling in wetlands remains poorly understood. We investigated how variation in macrophyte growth form (emerging, submerged and bryophyte) combinations and species mixtures affect nitrogen (N) removal from the water and N accumulation in plant biomass. We conducted a wetland mesocosm experiment for 100 days during July-September 2015. Twelve species were grown in mono- and in two-species mixed cultures for a total of 32 single and two-growth form combinations. Nitrogen removal from the water was quantified on three occasions during the experiment, while N accumulation in plant biomass was determined following termination of the experiment. The number of species and growth forms present increased N removal and accumulation. The growth form combinations of emerging and bryophyte species showed the highest N accumulation and N removal from water, followed by combinations of emerging species. By contrast, submerged species growing in the presence of emerging or other submerged species showed the lowest levels of N accumulation and N removal. Temporal variation in N removal also differed among growth form combinations: N removal was highest for emerging-bryophyte combinations in July, but peaked for the emerging-submerged and emerging-bryophyte combinations in August. Indeed, the occurrence of complementarity among macrophyte species, particularly in combinations of bryophyte and emerging species, enhanced N removal and uptake during the entire growing season. Our study highlights the importance of bryophytes, which have been neglected in research on nutrient cycling in wetlands, for aquatic N cycling, especially given their worldwide distribution across biomes. Overall, our findings point towards the potential important role of the diversity of macrophyte growth forms in regulating key ecosystem processes related to N cycling in wetlands.

Seasonal shift of diet in bank voles explains trophic fate of anthropogenic osmium?

Diet shifts are common in mammals and birds, but little is known about how such shifts along the food web affect contaminant exposure. Voles are staple food for many mammalian and avian predators. There is therefore a risk of transfer of contaminants accumulated in voles within the food chain. Osmium is one of the rarest earth elements with osmium tetroxide (OsO₄) as the most toxic vapor-phase airborne contaminant. Anthropogenic OsO₄ accumulates in fruticose lichens that are important winter food of bank voles (Myodes glareolus). Here, we test if a) anthropogenic osmium accumulates in bank voles in winter, and b) accumulation rates and concentrations are lower in autumn when the species is mainly herbivorous. Our study, performed in a boreal forest impacted by anthropogenic osmium, supported the hypotheses for all studied tissues (kidney, liver, lung, muscle and spleen) in 50 studied bank voles. In autumn, osmium concentrations in bank voles were even partly similar to those in the graminivorous field vole (Microtus agrestis; n=14). In autumn but not in late winter/early spring, osmium concentrations were generally negatively correlated with body weight and root length of the first mandible molar, i.e. proxies of bank vole age. Identified negative correlations between organ-to-body weight ratios and osmium concentrations in late winter/early spring indicate intoxication. Our results suggest unequal accumulation risk for predators feeding on different cohorts of bank voles.

Spatial prediction and validation of zoonotic hazard through micro-habitat properties: where does Puumala hantavirus hole - up?

Background

To predict the risk of infectious diseases originating in wildlife, it is important to identify habitats that allow the co-occurrence of pathogens and their hosts. Puumala hantavirus (PUUV) is a directly-transmitted RNA virus that causes hemorrhagic fever in humans, and is carried and transmitted by the bank vole (Myodes glareolus). In northern Sweden, bank voles undergo 3–4 year population cycles, during which their spatial distribution varies greatly.

Methods

We used boosted regression trees; a technique inspired by machine learning, on a 10 – year time-series (fall 2003–2013) to develop a spatial predictive model assessing seasonal PUUV hazard using micro-habitat variables in a landscape heavily modified by forestry. We validated the models in an independent study area approx. 200 km away by predicting seasonal presence of infected bank voles in a five-year-period (2007–2010 and 2015).

Results

The distribution of PUUV-infected voles varied seasonally and inter-annually. In spring, micro-habitat variables related to cover and food availability in forests predicted both bank vole and infected bank vole presence. In fall, the presence of PUUV-infected voles was generally restricted to spruce forests where cover was abundant, despite the broad landscape distribution of bank voles in general. We hypothesize that the discrepancy in distribution between infected and uninfected hosts in fall, was related to higher survival of PUUV and/or PUUV-infected voles in the environment, especially where cover is plentiful.

Conclusions

Moist and mesic old spruce forests, with abundant cover such as large holes and bilberry shrubs, also providing food, were most likely to harbor infected bank voles. The models developed using long-term and spatially extensive data can be extrapolated to other areas in northern Fennoscandia. To predict the hazard of directly transmitted zoonoses in areas with unknown risk status, models based on micro-habitat variables and developed through machine learning techniques in well-studied systems, could be used.

Electronic supplementary material

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

Dampening of population cycles in voles affects small mammal community structure, decreases diversity, and increases prevalence of a zoonotic disease

Long-term decline and depression of density in cyclic small rodents is a recent widespread phenomenon. These observed changes at the population level might have cascading effects at the ecosystem level. Here, we assessed relationships between changing boreal landscapes and biodiversity changes of small mammal communities. We also inferred potential effects of observed community changes for increased transmission risk of Puumala virus (PUUV) spread, causing the zoonotic disease nephropatica epidemica in humans. Analyses were based on long-term (1971-2013) monitoring data of shrews and voles representing 58 time series in northern Sweden. We calculated richness, diversity, and evenness at alpha, beta, and gamma level, partitioned beta diversity into turnover (species replacement) and nestedness (species addition/removal), used similarity percentages (SIMPER) analysis to assess community structure, and calculated the cumulated number of PUUV-infected bank voles and average PUUV prevalence (percentage of infected bank voles) per vole cycle. Alpha, beta, and gamma richness and diversity of voles, but not shrews, showed long-term trends that varied spatially. The observed patterns were associated with an increase in community contribution of bank vole (Myodes glareolus), a decrease of gray-sided vole (M. rufocanus) and field vole (Microtus agrestis) and a hump-shaped variation in contribution of common shrew (Sorex araneus). Long-term biodiversity changes were largely related to changes in forest landscape structure. Number of PUUV-infected bank voles in spring was negatively related to beta and gamma diversity, and positively related to turnover of shrews (replaced by voles) and to community contribution of bank voles. The latter was also positively related to average PUUV prevalence in spring. We showed that long-term changes in the boreal landscape contributed to explain the decrease in biodiversity and the change in structure of small mammal communities. In addition, our results suggest decrease in small mammal diversity to have knock-on effects on dynamics of infectious diseases among small mammals with potential implications for disease transmission to humans.

Sublethal Lead Exposure Alters Movement Behavior in Free-Ranging Golden Eagles

Lead poisoning of animals due to ingestion of fragments from lead-based ammunition in carcasses and offal of shot wildlife is acknowledged globally and raises great concerns about potential behavioral effects leading to increased mortality risks. Lead levels in blood were correlated with progress of the moose hunting season. Based on analyses of tracking data, we found that even sublethal lead concentrations in blood (25 ppb, wet weight), can likely negatively affect movement behavior (flight height and movement rate) of free-ranging scavenging Golden Eagles (Aquila chrysaetos). Lead levels in liver of recovered post-mortem analyzed eagles suggested that sublethal exposure increases the risk of mortality in eagles. Such adverse effects on animals are probably common worldwide and across species, where game hunting with lead-based ammunition is widespread. Our study highlights lead exposure as a considerably more serious threat to wildlife conservation than previously realized and suggests implementation of bans of lead ammunition for hunting.

Declining ecosystem health and the dilution effect

The “dilution effect” implies that where species vary in susceptibility to infection by a pathogen, higher diversity often leads to lower infection prevalence in hosts. For directly transmitted pathogens, non-host species may “dilute” infection directly (1) and indirectly (2). Competitors and predators may (1) alter host behavior to reduce pathogen transmission or (2) reduce host density. In a well-studied system, we tested the dilution of the zoonotic Puumala hantavirus (PUUV) in bank voles (Myodes glareolus) by two competitors and a predator. Our study was based on long-term PUUV infection data (2003–2013) in northern Sweden. The field vole (Microtus agrestis) and the common shrew (Sorex araneus) are bank vole competitors and Tengmalm’s owl (Aegolius funereus) is a main predator of bank voles. Infection probability in bank voles decreased when common shrew density increased, suggesting that common shrews reduced PUUV transmission. Field voles suppressed bank vole density in meadows and clear-cuts and indirectly diluted PUUV infection. Further, Tengmalm’s owl decline in 1980–2013 may have contributed to higher PUUV infection rates in bank voles in 2003–2013 compared to 1979–1986. Our study provides further evidence for dilution effect and suggests that owls may have an important role in reducing disease risk.

Selective predation on hantavirus-infected voles by owls and confounding effects from landscape properties

It has been suggested that predators may protect human health through reducing disease-host densities or selectively preying on infected individuals from the population. However, this has not been tested empirically. We hypothesized that Tengmalm's owl (Aegolius funereus) selectively preys on hantavirus-infected individuals of its staple prey, the bank vole (Myodes glareolus). Bank voles are hosts of Puumala hantavirus, which causes a form of hemorrhagic fever in humans. Selective predation by owls on infected voles may reduce human disease risk. We compared the prevalence of anti-Puumala hantavirus antibodies (seroprevalence), in bank voles cached by owls in nest boxes to seroprevalence in voles trapped in closed-canopy forest around each nest box. We found no general difference in seroprevalence. Forest landscape structure could partly account for the observed patterns in seroprevalence. Only in more connected forest patches was seroprevalence in bank voles cached in nest boxes higher than seroprevalence in trapped voles. This effect disappeared with increasing forest patch isolation, as seroprevalence in trapped voles increased with forest patch isolation, but did not in cached voles. Our results suggest a complex relationship between zoonotic disease prevalence in hosts, their predators, and landscape structure. Some mechanisms that may have caused the seroprevalence patterns in our results include higher bank vole density in isolated forest patches. This study offers future research potential to shed further light on the contribution of predators and landscape properties to human health.

Relative importance of plant uptake and plant associated denitrification for removal of nitrogen from mine drainage in sub-arctic wetlands

Reactive nitrogen (N) species released from undetonated ammonium-nitrate based explosives used in mining or other blasting operations are an emerging environmental problem. Wetlands are frequently used to treat N-contaminated water in temperate climate, but knowledge on plant-microbial interactions and treatment potential in sub-arctic wetlands is limited. Here, we compare the relative importance of plant uptake and denitrification among five plant species commonly occurring in sub-arctic wetlands for removal of N in nitrate-rich mine drainage in northern Sweden. Nitrogen uptake and plant associated potential denitrification activity and genetic potential for denitrification based on quantitative PCR of the denitrification genes nirS, nirK, nosZI and nosZII were determined in plants growing both in situ and cultivated in a growth chamber. The growth chamber and in situ studies generated similar results, suggesting high relevance and applicability of results from growth chamber experiments. We identified denitrification as the dominating pathway for N-removal and abundances of denitrification genes were strong indicators of plant associated denitrification activity. The magnitude and direction of the effect differed among the plant species, with the aquatic moss Drepanocladus fluitans showing exceptionally high ratios between denitrification and uptake rates, compared to the other species. However, to acquire realistic estimates of N-removal potential of specific wetlands and their associated plant species, the total plant biomass needs to be considered. The species-specific plant N-uptake and abundance of denitrification genes on the root or plant surfaces were affected by the presence of other plant species, which show that both multi- and inter-trophic interactions are occurring. Future studies on N-removal potential of wetland plant species should consider how to best exploit these interactions in sub-arctic wetlands.

Impact of Beaver Pond Colonization History on Methylmercury Concentrations in Surface Water

Elevated concentrations of methylmercury (MeHg) in freshwater ecosystems are of major environmental concern in large parts of the northern hemisphere. Beaver ponds have been identified as a potentially important source of MeHg. The role of beavers might be especially pronounced in large parts of Europe, where beaver populations have expanded rapidly following near-extirpation. This study evaluates the role of the age and colonization history (encompassing patterns of use and reuse) of ponds constructed by the Eurasian beaver Castor fiber in regulating MeHg concentrations in Swedish streams. In 12 beaver systems located in three regions, we quantified MeHg concentrations together with other relevant parameters on five occasions per year in 2012-2013. Five were pioneer systems, inundated for the first time since beaver extirpation, and seven were recolonized, with dams reconstructed by newly recolonizing beavers. MeHg concentrations in pioneer but not in recolonized beaver systems were up to 3.5 fold higher downstream than upstream of the ponds, and varied between seasons and years. Our results show that pioneer inundation by beavers can increase MeHg concentrations in streams, but that this effect is negligible when dams are reconstructed on previously used ponds. We therefore expect that the recovery and expansion of beavers in the boreal system will only have a transitional effect on MeHg in the environment.

Dynamics and drivers of hantavirus prevalence in rodent populations

Human encroachment on wildlife habitats has contributed to the emergence of several zoonoses. Pathogenic hantaviruses are hosted by rodents and cause severe diseases in the Americas and Eurasia. We reviewed several factors that potentially drive prevalence (the proportion of infected rodents) in host populations. These include demography, behavior, host density, small mammal diversity, predation, and habitat and landscape characteristics. This review is the first to include a quantitative summary of the literature investigating hantavirus prevalence in rodents. Demographic structure and density were investigated the most and predation the least. Reported effects of demographic structure and small mammal diversity were consistent, whereby reproductive males were most likely to be infected and prevalence decreased with small mammal diversity. The influences of habitat and landscape properties are often complex and indirect. The relationship between density and prevalence merits more investigation. Most hantavirus hosts are habitat generalists and their control is challenging. Incorporating all potential factors and their interactions is essential to understanding and controlling infection in host populations.

The importance of bank vole density and rainy winters in predicting nephropathia epidemica incidence in Northern Sweden

Pathogenic hantaviruses (family Bunyaviridae, genus Hantavirus) are rodent-borne viruses causing hemorrhagic fever with renal syndrome (HFRS) in Eurasia. In Europe, there are more than 10,000 yearly cases of nephropathia epidemica (NE), a mild form of HFRS caused by Puumala virus (PUUV). The common and widely distributed bank vole (Myodes glareolus) is the host of PUUV. In this study, we aim to explain and predict NE incidence in boreal Sweden using bank vole densities. We tested whether the number of rainy days in winter contributed to variation in NE incidence. We forecast NE incidence in July 2013–June 2014 using projected autumn vole density, and then considering two climatic scenarios: 1) rain-free winter and 2) winter with many rainy days. Autumn vole density was a strong explanatory variable of NE incidence in boreal Sweden in 1990–2012 (R² = 79%, p<0.001). Adding the number of rainy winter days improved the model (R² = 84%, p<0.05). We report for the first time that risk of NE is higher in winters with many rainy days. Rain on snow and ground icing may block vole access to subnivean space. Seeking refuge from adverse conditions and shelter from predators, voles may infest buildings, increasing infection risk. In a rainy winter scenario, we predicted 812 NE cases in boreal Sweden, triple the number of cases predicted in a rain-free winter in 2013/2014. Our model enables identification of high risk years when preparedness in the public health sector is crucial, as a rainy winter would accentuate risk.

Predicting grey-sided vole occurrence in northern Sweden at multiple spatial scales

Forestry is continually changing the habitats for many forest-dwelling species around the world. The grey-sided vole (Myodes rufocanus) has declined since the 1970s in forests of northern Sweden. Previous studies suggested that this might partly be caused by reduced focal forest patch size due to clear-cutting. Proximity and access to old pine forest and that microhabitats often contains stones have also been suggested previously but never been evaluated at multiple spatial scales. In a field study in 2010–2011 in northern Sweden, we investigated whether occurrence of grey-sided voles would be higher in (1) large focal patches of >60 years old forest, (2) in patches with high connectivity to surrounding patches, and (3) in patches in proximity to stone fields. We trapped animals in forest patches in two study areas (Västerbotten and Norrbotten). At each trap station, we surveyed structural microhabitat characteristics. Landscape-scale features were investigated using satellite-based forest data combined with geological maps. Unexpectedly, the vole was almost completely absent in Norrbotten. The trap sites in Norrbotten had a considerably lower amount of stone holes compared with sites with voles in Västerbotten. We suggest this might help to explain the absence in Norrbotten. In Västerbotten, the distance from forest patches with voles to stone fields was significantly shorter than from patches without voles. In addition, connectivity to surrounding patches and size of the focal forest patches was indeed related to the occurrence of grey-sided voles, with connectivity being the overall best predictor. Our results support previous findings on the importance of large forest patches, but also highlight the importance of connectivity for occurrence of grey-sided voles. The results further suggest that proximity to stone fields increase habitat quality of the forests for the vole and that the presence of stone fields enhances the voles' ability to move between nearby forest patches through the matrix.

Europe-wide dampening of population cycles in keystone herbivores

Suggestions of collapse in small herbivore cycles since the 1980s have raised concerns about the loss of essential ecosystem functions. Whether such phenomena are general and result from extrinsic environmental changes or from intrinsic process stochasticity is currently unknown. Using a large compilation of time series of vole abundances, we demonstrate consistent cycle amplitude dampening associated with a reduction in winter population growth, although regulatory processes responsible for cyclicity have not been lost. The underlying syndrome of change throughout Europe and grass-eating vole species suggests a common climatic driver. Increasing intervals of low-amplitude small herbivore population fluctuations are expected in the future, and these may have cascading impacts on trophic webs across ecosystems.

Modelling nitrogen transformations in waters receiving mine effluents

This paper presents a biogeochemical model developed for a clarification pond receiving ammonium nitrogen rich discharge water from the Boliden concentration plant located in northern Sweden. Present knowledge about nitrogen (N) transformations in lakes is compiled in a dynamic model that calculates concentrations of the six N species (state variables) ammonium-N (N(am)), nitrate-N (N(ox)), dissolved organic N in water (N(org)), N in phytoplankton (N(pp)), in macrophytes (N(mp)) and in sediment (N(sed)). It also simulates the rate of 16 N transformation processes occurring in the water column and sediment as well as water-sediment and water-atmosphere interactions. The model was programmed in the software Powersim using 2008 data, whilst validation was performed using data from 2006 to 2007. The sensitivity analysis showed that the state variables are most sensitive to changes in the coefficients related to the temperature dependence of the transformation processes. A six-year simulation of N(am) showed stable behaviour over time. The calibrated model rendered coefficients of determination (R(2)) of 0.93, 0.79 and 0.86 for N(am), N(ox) and N(org), respectively. Performance measures quantitatively expressing the deviation between modelled and measured data resulted in values close to zero, indicating a stable model structure. The simulated denitrification rate was on average five times higher than the ammonia volatilisation rate and about three times higher than the permanent burial of N(sed) and, hence, the most important process for the permanent removal of N. The model can be used to simulate possible measures to reduce the nitrogen load and, after some modification and recalibration, it can be applied at other mine sites affected by N rich effluents.

Liquid culture mass production of biocontrol nematodes, Heterorhabditis bacteriophora (Nematoda: Rhabditida): improved timing of dauer juvenile inoculation

Heterorhabditis bacteriophora is used in biological control of soil-borne insect pests in horticulture and turf. Mass production is carried out in monoxenic liquid cultures pre-incubated with the symbiont of the nematodes, the bacterium Photorhabdus luminescens, before nematode dauer juveniles (DJ) are inoculated. As a response to bacterial food signals, the DJ recover from the developmentally arrested dauer stage, grow to adults and produce DJ offspring. Variable DJ recovery after inoculation into cultures of P. luminescens often causes process failure due to low numbers of adult nematodes in the medium. In order to enhance DJ recovery, improve nematode population management and increase yields, the optimal timing for DJ inoculation was sought. The process parameter pH and respiration quotient (RQ) were recorded in order to test whether changes can be used to identify the best moment for DJ inoculation. When DJ were inoculated during the lag and early logarithmic growth phases of P. luminescens cultures, DJ recovery was low and almost no nematode reproduction was obtained. High populations of P. luminescens phase variants were recorded. Recovery and yields increased when DJ were inoculated during the latter log phase during which the RQ dropped to values <0.8 and the pH reached a maximum. The highest DJ recovery and yields were observed in cultures that were inoculated during the late stationary growth phase. This period started with the increase of the pH after its distinct minimum at pH <8.0. Thus optimal timing for DJ inoculation can be defined through monitoring of the pH in the P. luminescens culture.