Age truncation due to disease shrinks metapopulation viability for amphibians

Metapopulations often exist in a fragile balance between local extinctions and (re)colonisations, in which case emerging threats that alter species vital rates may drastically increase metapopulation extinction risk. We combined empirical data with metapopulation simulations to examine how demographic shifts associated with amphibian chytrid fungus (Batrachochytrium dendrobatidis, Bd) have altered metapopulation viability for threatened amphibians in Australia. Comparing the ages of museum specimens collected before Bd emerged in Australia with individuals from geographically matched remnant populations revealed significant truncation of age structures post-Bd, with a halving of annual adult survival probabilities. Spatially realistic metapopulation modelling demonstrated that reduced adult survival led to major reductions in the parameter space over which persistence was possible for the focal species, with contractions to landscapes with higher landscape connectivity, lower environmental stochasticity and considerably higher recruitment rates. Metapopulation persistence post-Bd required greater landscape connectivity than pre-Bd. This arises from a landscape-level analogue of compensatory recruitment at the population level, in which higher (re)colonisation rates can offset more frequent local extinctions, enabling persistence of amphibians susceptible to Bd. Interactions between recruitment rate, environmental stochasticity and landscape connectivity were also more important for metapopulation persistence post-Bd. Higher recruitment was required to mitigate the impacts of environmental stochasticity, and higher landscape connectivity was required to mitigate the impacts of environmental stochasticity and poor recruitment. Increased reliance on these interdependencies shrunk the parameter space over which metapopulations could persist post-Bd. Our study demonstrates that emerging threats that alter species vital rates can drastically reduce the capacity of certain environments to support metapopulations. For our focal species, reductions in adult survival rates due to Bd produced major reductions in the conditions under which persistence was possible, providing a mechanistic insight into the processes underpinning observed range and niche contractions of amphibians impacted by this pathogen. More broadly, our study illustrates how environmentally mediated host resilience can enable persistence following the emergence of novel pathogens. This pathway to persistence is worthy of greater attention on both conceptual and applied grounds.

Front propagation in a spatial system of weakly interacting networks

We consider the spread of epidemic in a spatial metapopulation system consisting of weakly interacting patches. Each local patch is represented by a network with a certain node degree distribution and individuals can migrate between neighboring patches. Stochastic particle simulations of the SIR model show that after a short transient, the spatial spread of epidemic has a form of a propagating front. A theoretical analysis shows that the speed of front propagation depends on the effective diffusion coefficient and on the local proliferation rate similarly to fronts described by the Fisher-Kolmogorov equation. To determine the speed of front propagation, first, the early-time dynamics in a local patch is computed analytically by employing degree based approximation for the case of a constant disease duration. The resulting delay differential equation is solved for early times to obtain the local growth exponent. Next, the reaction diffusion equation is derived from the effective master equation and the effective diffusion coefficient and the overall proliferation rate are determined. Finally, the fourth order derivative in the reaction diffusion equation is taken into account to obtain the discrete correction to the front propagation speed. The analytical results are in a good agreement with the results of stochastic particle simulations.

Intracellular common gardens reveal niche differentiation in transposable element community during bacterial adaptive evolution

The distribution and abundance of transposable elements across the tree of life have significantly shaped the evolution of cellular organisms, but the underlying mechanisms shaping these ecological patterns remain elusive. Here we establish a "common garden" approach to study causal ecological interactions between a xenogeneic conditional lethal sacB gene and the community of transposable insertion sequences (ISs) in a multipartite prokaryote genome. Xenogeneic sacB of low, medium, or high GC content was individually inserted into three replicons of a model bacterium Sinorhizobium fredii, and exhibited replicon- and GC-dependent variation in genetic stability. This variation was largely attributable to multidimensional niche differentiation for IS community members. The transposition efficiency of major active ISs depended on the nucleoid-associated xenogeneic silencer MucR. Experimentally eliminating insertion activity of specific ISs by deleting MucR strongly demonstrated a dominant role of niche differentiation among ISs. This intracellular common garden approach in the experimental evolution context allows not only for evaluating genetic stability of natural and synthetic xenogeneic genes of different sequence signatures in host cells but also for tracking and testing causal relationships in unifying ecological principles in genome ecology.

Exploring links between climatic predictability and the evolution of within- and transgenerational plasticity

In variable environments, phenotypic plasticity can increase fitness by providing tight environment-phenotype matching. However, adaptive plasticity is expected to evolve only when the future selective environment can be predicted based on the prevailing conditions. That is, the juvenile environment should be predictive of the adult environment (within-generation plasticity) or the parental environment should be predictive of the offspring environment (transgenerational plasticity). Moreover, the environmental predictability can also shape transient responses such as stress response in an adaptive direction. Here, we test links between environmental predictability and the evolution of adaptive plasticity by combining time series analyses and a common garden experiment using temperature as a stressor in a temperate butterfly (Melitaea cinxia). Time series analyses revealed that across season fluctuations in temperature over 48 years are overall predictable. However, within the growing season, temperature fluctuations showed high heterogeneity across years with low autocorrelations and the timing of temperature peaks were asynchronous. Most life-history traits showed strong within-generation plasticity for temperature and traits such as body size and growth rate broke the temperature-size rule. Evidence for transgenerational plasticity, however, was weak and detected for only two traits each in an adaptive and non-adaptive direction. We suggest that the low predictability of temperature fluctuations within the growing season likely disfavors the evolution of adaptive transgenerational plasticity but instead favors strong within-generation plasticity.

Stability of patch-turnover relationships under equilibrium and nonequilibrium metapopulation dynamics driven by biogeography

Two controversial tenets of metapopulation biology are whether patch quality and the surrounding matrix are more important to turnover (colonisation and extinction) than biogeography (patch area and isolation) and whether factors governing turnover during equilibrium also dominate nonequilibrium dynamics. We tested both tenets using 18 years of surveys for two secretive wetland birds, black and Virginia rails, during (1) a period of equilibrium with stable occupancy and (2) after drought and arrival of West Nile Virus (WNV), which resulted in WNV infections in rails, increased extinction and decreased colonisation probabilities modified by WNV, nonequilibrium dynamics for both species and occupancy decline for black rails. Area (primarily) and isolation (secondarily) drove turnover during both stable and unstable metapopulation dynamics, greatly exceeding the effects of patch quality and matrix conditions. Moreover, slopes between turnover and patch characteristics changed little between equilibrium and nonequilibrium, confirming the overriding influences of biogeographic factors on turnover.

Effects of environment and genotype on dispersal differ across departure, transfer and settlement in a butterfly metapopulation

Active dispersal is driven by extrinsic and intrinsic factors at the three stages of departure, transfer and settlement. Most empirical studies capture only one stage of this complex process, and knowledge of how much can be generalized from one stage to another remains unknown. Here we use genetic assignment tests to reconstruct dispersal across 5 years and 232 habitat patches of a Glanville fritillary butterfly (Melitaea cinxia) metapopulation. We link individual dispersal events to weather, landscape structure, size and quality of habitat patches, and individual genotype to identify the factors that influence the three stages of dispersal and post-settlement survival. We found that nearly all tested factors strongly affected departure probabilities, but that the same factors explained very little variation in realized dispersal distances. Surprisingly, we found no effect of dispersal distance on post-settlement survival. Rather, survival was influenced by weather conditions, quality of the natal habitat patch, and a strong interaction between genotype and occupancy status of the settled habitat patch, with more mobile genotypes having higher survival as colonists rather than as immigrants. Our work highlights the multi-causality of dispersal and that some dispersal costs can only be understood by considering extrinsic and intrinsic factors and their interaction across the entire dispersal process.

Improved chromosome-level genome assembly of the Glanville fritillary butterfly (Melitaea cinxia) integrating Pacific Biosciences long reads and a high-density linkage map

BACKGROUND

The Glanville fritillary (Melitaea cinxia) butterfly is a model system for metapopulation dynamics research in fragmented landscapes. Here, we provide a chromosome-level assembly of the butterfly's genome produced from Pacific Biosciences sequencing of a pool of males, combined with a linkage map from population crosses.

RESULTS

The final assembly size of 484 Mb is an increase of 94 Mb on the previously published genome. Estimation of the completeness of the genome with BUSCO indicates that the genome contains 92-94% of the BUSCO genes in complete and single copies. We predicted 14,810 genes using the MAKER pipeline and manually curated 1,232 of these gene models.

CONCLUSIONS

The genome and its annotated gene models are a valuable resource for future comparative genomics, molecular biology, transcriptome, and genetics studies on this species.

The effect of summer drought on the predictability of local extinctions in a butterfly metapopulation

The ecological impacts of extreme climatic events on population dynamics and community composition are profound and predominantly negative. Using extensive data of an ecological model system, we tested whether predictions from ecological models remain robust when environmental conditions are outside the bounds of observation. We observed a 10-fold demographic decline of the Glanville fritillary butterfly (Melitaea cinxia) metapopulation on the Åland islands, Finland in the summer of 2018 and used climatic and satellite data to demonstrate that this year was an anomaly with low climatic water balance values and low vegetation productivity indices across Åland. Population growth rates were strongly associated with spatiotemporal variation in climatic water balance. Covariates shown previously to affect the extinction probability of local populations in this metapopulation were less informative when populations were exposed to severe drought during the summer months. Our results highlight the unpredictable responses of natural populations to extreme climatic events.

Earth and field observations underpin metapopulation dynamics in complex landscapes: Near-term study on carabids

Understanding risks to biodiversity requires predictions of the spatial distribution of species adapting to changing ecosystems and, to that end, Earth observations integrating field surveys prove essential as they provide key numbers for assessing landscape-wide biodiversity scenarios. Here, we develop, and apply to a relevant case study, a method suited to merge Earth/field observations with spatially explicit stochastic metapopulation models to study the near-term ecological dynamics of target species in complex terrains. Our framework incorporates the use of species distribution models for a reasoned estimation of the initial presence of the target species and accounts for imperfect and incomplete detection of the species presence in the study area. It also uses a metapopulation fitness function derived from Earth observation data subsuming the ecological niche of the target species. This framework is applied to contrast occupancy of two species of carabids (Pterostichus flavofemoratus, Carabus depressus) observed in the context of a large ecological monitoring program carried out within the Gran Paradiso National Park (GPNP, Italy). Results suggest that the proposed framework may indeed exploit the hallmarks of spatially explicit ecological approaches and of remote Earth observations. The model reproduces well the observed in situ data. Moreover, it projects in the near term the two species' presence both in space and in time, highlighting the features of the metapopulation dynamics of colonization and extinction, and their expected trends within verifiable timeframes.