A simple kit to use computational notebooks for more openness, reproducibility, and productivity in research.

The ubiquitous use of computational work for data generation, processing, and modeling increased the importance of digital documentation in improving research quality and impact. Computational notebooks are files that contain descriptive text, as well as code and its outputs, in a single, dynamic, and visually appealing file that is easier to understand by nonspecialists. Traditionally used by data scientists when producing reports and informing decision-making, the use of this tool in research publication is not common, despite its potential to increase research impact and quality. For a single study, the content of such documentation partially overlaps with that of classical lab notebooks and that of the scientific manuscript reporting the study. Therefore, to minimize the amount of work required to manage all the files related to these contents and optimize their production, we present a starter kit to facilitate the implementation of computational notebooks in the research process, including publication. The kit contains the template of a computational notebook integrated into a research project that employs R, Python, or Julia. Using examples of ecological studies, we show how computational notebooks also foster the implementation of principles of Open Science, such as reproducibility and traceability. The kit is designed for beginners, but at the end we present practices that can be gradually implemented to develop a fully digital research workflow. Our hope is that such minimalist yet effective starter kit will encourage researchers to adopt this practice in their workflow, regardless of their computational background.

'Keeping the kids at home' can limit the persistence of contagious pathogens in social animals.

Social networks are considered to be 'highly modular' when individuals within one module are more connected to each other than they are to individuals in other modules. It is currently unclear how highly modular social networks influence the persistence of contagious pathogens that generate lifelong immunity in their hosts when between-group interactions are age dependent. This trait occurs in social species with communal nurseries, where juveniles are reared together for a substantial period in burrows or similar forms of containment and are thus in isolation from contact with individuals in other social groups. Our main objective was to determine whether, and to what extent, such age-dependent patterns of between-group interactions consistently increased the fade-out probability of contagious pathogens that generate lifelong immunity in their hosts. We hypothesised that in populations of species where juveniles are raised in communal nurseries, a high proportion of recovered adults in a group would form a 'protective barrier' around susceptible juveniles against pathogen transmission, thereby increasing the probability of epidemic fade-out in the population. To test this idea, we implemented a spatially implicit individual-based susceptible-infected-recovered (SIR) model for a large range of generic host and pathogen traits. Our results indicated that (a) the probability of epidemic fade-out was consistently higher in populations with communal nurseries, especially for highly contagious pathogens (high basic reproduction number, R0 ) and (b) communal nurseries can counteract the cost of group living in terms of infection risk to a greater extent than variation in other traits. We discuss our findings in relation to herd immunity and outline the importance of considering the network structure of a given host population before implementing management measures such as vaccinations, since interventions focused on individuals with high between-group contact should be particularly effective for controlling pathogen spread in hosts with communal nurseries.

The dimensionality of stability depends on disturbance type.

Ecosystems respond in various ways to disturbances. Quantifying ecological stability therefore requires inspecting multiple stability properties, such as resistance, recovery, persistence and invariability. Correlations among these properties can reduce the dimensionality of stability, simplifying the study of environmental effects on ecosystems. A key question is how the kind of disturbance affects these correlations. We here investigated the effect of three disturbance types (random, species-specific, local) applied at four intensity levels, on the dimensionality of stability at the population and community level. We used previously parameterized models that represent five natural communities, varying in species richness and the number of trophic levels. We found that disturbance type but not intensity affected the dimensionality of stability and only at the population level. The dimensionality of stability also varied greatly among species and communities. Therefore, studying stability cannot be simplified to using a single metric and multi-dimensional assessments are still to be recommended.

Seasonal host life-history processes fuel disease dynamics at different spatial scales.

Understanding the drivers underlying disease dynamics is still a major challenge in disease ecology, especially in the case of long-term disease persistence. Even though there is a strong consensus that density-dependent factors play an important role for the spread of diseases, the main drivers are still discussed and, more importantly, might differ between invasion and persistence periods. Here, we analysed long-term outbreak data of classical swine fever, an important disease in both wild boar and livestock, prevalent in the wild boar population from 1993 to 2000 in Mecklenburg-Vorpommern, Germany. We report outbreak characteristics and results from generalized linear mixed models to reveal what factors affected infection risk on both the landscape and the individual level. Spatiotemporal outbreak dynamics showed an initial wave-like spread with high incidence during the invasion period followed by a drop of incidence and an increase in seroprevalence during the persistence period. Velocity of spread increased with time during the first year of outbreak and decreased linearly afterwards, being on average 7.6 km per quarter. Landscape- and individual-level analyses of infection risk indicate contrasting seasonal patterns. During the persistence period, infection risk on the landscape level was highest during autumn and winter seasons, probably related to spatial behaviour such as increased long-distance movements and contacts induced by rutting and escaping movements. In contrast, individual-level infection risk peaked in spring, probably related to the concurrent birth season leading to higher densities, and was significantly higher in piglets than in reproductive animals. Our findings highlight that it is important to investigate both individual- and landscape-level patterns of infection risk to understand long-term persistence of wildlife diseases and to guide respective management actions. Furthermore, we highlight that exploring different temporal aggregation of the data helps to reveal important seasonal patterns, which might be masked otherwise.

Resource asynchrony and landscape homogenization as drivers of virulence evolution: The case of a directly transmitted disease in a social host.

Throughout the last decades, the emergence of zoonotic diseases and the frequency of disease outbreaks have increased substantially, fuelled by habitat encroachment and vectors overlapping with more hosts due to global change. The virulence of pathogens is one key trait for successful invasion. In order to understand how global change drivers such as habitat homogenization and climate change drive pathogen virulence evolution, we adapted an established individual-based model of host-pathogen dynamics. Our model simulates a population of social hosts affected by a directly transmitted evolving pathogen in a dynamic landscape. Pathogen virulence evolution results in multiple strains in the model that differ in their transmission capability and lethality. We represent the effects of global change by simulating environmental changes both in time (resource asynchrony) and space (homogenization). We found an increase in pathogenic virulence and a shift in strain dominance with increasing landscape homogenization. Our model further indicated that lower virulence is dominant in fragmented landscapes, although pulses of highly virulent strains emerged under resource asynchrony. While all landscape scenarios favoured co-occurrence of low- and high-virulent strains, the high-virulence strains capitalized on the possibility for transmission when host density increased and were likely to become dominant. With asynchrony likely to occur more often due to global change, our model showed that a subsequent evolution towards lower virulence could lead to some diseases becoming endemic in their host populations.

Host weight, seasonality and anthropogenic factors contribute to parasite community differences between urban and rural foxes.

Pathogens often occur at different prevalence along environmental gradients. This is of particular importance for gradients of anthropogenic impact such as rural-urban transitions presenting a changing interface between humans and wildlife. The assembly of parasite communities is affected by both the external environmental conditions and individual host characteristics. Hosts with low body weight (smaller individuals or animals with poor body condition) might be more susceptible to infection. Furthermore, parasites' mode of transmission might affect their occurrence: rural environments with better availability of intermediate hosts might favour trophic transmission, while urban environments, typically with dense definitive host populations, might favour direct transmission. We here study helminth communities (141 intestinal samples) within the red fox (Vulpes vulpes), a synanthropic host, using DNA metabarcoding of multiple marker genes. We analysed the effect of urbanisation, seasonality and host-intrinsic (weight, sex) variables on helminth communities. Helminth species richness increased in foxes with lower body weight and in winter and spring. Season and urbanisation, however, had strong effects on the community composition, i.e., on the identity of the detected species. Surprisingly, transmission in two-host life cycles (trophic transmission) was more pronounced in urban Berlin than in rural Brandenburg. This disagrees with the prevailing hypothesis that trophically transmitted helminths are less prevalent in urban areas than in rural areas. Generally, co-infestations with multiple helminths and high infection intensity are associated with lighter (younger, smaller or low body condition) animals. Both host-intrinsic traits and environmental drivers together shape parasite community composition and turnover along urban-rural gradients.

Movement can mediate temporal mismatches between resource availability and biological events in host-pathogen interactions.

Global change is shifting the timing of biological events, leading to temporal mismatches between biological events and resource availability. These temporal mismatches can threaten species' populations. Importantly, temporal mismatches not only exert strong pressures on the population dynamics of the focal species, but can also lead to substantial changes in pairwise species interactions such as host-pathogen systems. We adapted an established individual-based model of host-pathogen dynamics. The model describes a viral agent in a social host, while accounting for the host's explicit movement decisions. We aimed to investigate how temporal mismatches between seasonal resource availability and host life-history events affect host-pathogen coexistence, that is, disease persistence. Seasonal resource fluctuations only increased coexistence probability when in synchrony with the hosts' biological events. However, a temporal mismatch reduced host-pathogen coexistence, but only marginally. In tandem with an increasing temporal mismatch, our model showed a shift in the spatial distribution of infected hosts. It shifted from an even distribution under synchronous conditions toward the formation of disease hotspots, when host life history and resource availability mismatched completely. The spatial restriction of infected hosts to small hotspots in the landscape initially suggested a lower coexistence probability due to the critical loss of susceptible host individuals within those hotspots. However, the surrounding landscape facilitated demographic rescue through habitat-dependent movement. Our work demonstrates that the negative effects of temporal mismatches between host resource availability and host life history on host-pathogen coexistence can be reduced through the formation of temporary disease hotspots and host movement decisions, with implications for disease management under disturbances and global change.

Movement-mediated community assembly and coexistence.

Organismal movement is ubiquitous and facilitates important ecological mechanisms that drive community and metacommunity composition and hence biodiversity. In most existing ecological theories and models in biodiversity research, movement is represented simplistically, ignoring the behavioural basis of movement and consequently the variation in behaviour at species and individual levels. However, as human endeavours modify climate and land use, the behavioural processes of organisms in response to these changes, including movement, become critical to understanding the resulting biodiversity loss. Here, we draw together research from different subdisciplines in ecology to understand the impact of individual-level movement processes on community-level patterns in species composition and coexistence. We join the movement ecology framework with the key concepts from metacommunity theory, community assembly and modern coexistence theory using the idea of micro-macro links, where various aspects of emergent movement behaviour scale up to local and regional patterns in species mobility and mobile-link-generated patterns in abiotic and biotic environmental conditions. These in turn influence both individual movement and, at ecological timescales, mechanisms such as dispersal limitation, environmental filtering, and niche partitioning. We conclude by highlighting challenges to and promising future avenues for data generation, data analysis and complementary modelling approaches and provide a brief outlook on how a new behaviour-based view on movement becomes important in understanding the responses of communities under ongoing environmental change.