Identifying the factors which favour the establishment and spread of non-native species in novel environments is one of the keys to predict - and hence prevent or control - biological invasions. This includes biological factors (i.e. factors associated with the invasive species themselves), and one of the prevailing hypotheses is that some species traits may explain their impressive success to establish and spread in novel environments [1]. In animals, most research studies have focused on traits associated with fecundity, age at maturity, level of affiliation to humans or dispersal ability for instance. The “composite picture” of the perfect (i.e. successful) invader that has gradually emerged is a small-bodied animal strongly affiliated to human activities with high fecundity, high dispersal ability and a super high level of plasticity. Of course, the story is not that simple, and actually a perfect invader sometimes – if not often- takes another form… Carrying on to identify what makes a species a successful invader or not is hence still an important research axis with major implications.
In this manuscript, Charbonnel and collaborators [2] provide an interesting opportunity to gain novel insights into our understanding of (the) traits underlying invasion success. They nicely combine the power of Next-Generation Sequencing (NGS) with a clever comparative approach of two closely-related invasive rodents (the house mouse Mus musculus and the black rat Rattus rattus) in a common environment. They use this experimental design to test the appealing hypothesis that pathogens may be actors of the story, and may indirectly explain why some non-native species are so successful in invading novel habitats.
It is generally assumed that the community of pathogens encountered by non-native species in novel environments is different from that of their native area. On the one hand (the enemy-release hypothesis), it can be hypothesized that non-native species, when they arrive into a novel environment, will be relaxed from the pressure imposed by their native pathogens because local pathogens are not adapted (and hence do not infect) to this novel host. Because immune defence against pathogens is highly costly, non-native species establishing into a novel environment could hence reallocate these costs to other functions such as fecundity or dispersal apparatus. This scenario has been termed the “evolution of increased competitive ability” (EICA) hypothesis [3]. On the other hand (the EICA-refined hypothesis [4]), one can assume that invaders will encounter new pathogens in newly established areas, and will allocate energy toward cost-effective immune pathways to permit allocating a non-negligible amount of energy toward other functions. Finally, a last hypothesis (the “immune protection” hypothesis) assumes major changes in pathogen composition between native and invaded areas, which should lead to an overall increase in immune investment by the native species to successfully invade novel environments [4]. This last hypothesis suggests that only non-native species being able to take up the associated costs of immunity will be successful invaders.
The role of immunity in invasion success has yet been poorly investigated, mainly because of the difficulty to simultaneously analyse multiple immune pathways [4]. Charbonnel and collaborators [2] overpass this difficulty by screening all genes expressed (using a whole RNA sequencing approach) in an immune tissue: the spleen. They do so along the invasion routes of two sympatric invasive rodents in Africa and compare anciently and newly invaded areas (respectively). For one of the two species (the house mouse), they found a high number of immune-related genes to be up-regulated in newly invaded areas compared to anciently invaded areas. All categories of immune pathways (costly and cost-effective) were up-regulated, suggesting an overall increase in immune investment in the mouse, which corroborates the “immune protection” hypothesis. For the black rat, patterns of gene expression were somewhat different, with much less pronounced differentiation in gene expression between newly and anciently invaded areas. Among the few differentiated genes, a few were associated to immune responses and some of theses genes were even down-regulated in the newly invaded areas. This pattern may actually corroborate the EICA hypothesis, although it could alternatively suggest that stochastic processes (drift) associated to recent decrease in population size (which is expected during a colonisation event) are more important than selection imposed by pathogens in shaping patterns of immune gene expression.
Overall, this study [2] suggests (i) that immune-related traits are important in predicting invasion success and (ii) that two successful species with a similar invasion history and living in similar environments can use different life-history strategies to reach the same success. This later finding is particularly relevant and intriguing as it suggests that the traits and strategies deployed by species to colonise new habitats might actually be idiosyncratic, and that, if general trends actually emerge in regards of traits predicting the success of invaders, the devil might actually be into the details. Comparative studies are extremely important to identify the general rules and the specificities sustaining actual patterns, but these approaches are yet poorly used in biological invasions (at least empirically). The work presented by Charbonnel and colleagues [2] calls for future comparative studies performed at multiple spatial scales (native vs. non-native areas, anciently vs. recently invaded areas), multiple taxonomic resolutions and across multiple traits (to search for trade-offs), so that the success of invasive species can be properly understood and predicted.

References

[1] Jeschke, J. M., & Strayer, D. L. (2006). Determinants of vertebrate invasion success in Europe and North America. Global Change Biology, 12(9), 1608-1619. doi: 10.1111/j.1365-2486.2006.01213.x
[2] Blossey, B., & Notzold, R. (1995). Evolution of increased competitive ability in invasive nonindigenous plants: a hypothesis. Journal of Ecology, 83(5), 887-889. doi: 10.2307/2261425
[3] Charbonnel, N., Galan, M., Tatard, C., Loiseau, A., Diagne, C. A., Dalecky, A., Parrinello, H., Rialle, S., Severac, D., & Brouat, C. (2019). Differential immune gene expression associated with contemporary range expansion of two invasive rodents in Senegal. bioRxiv, 442160, ver. 5 peer-reviewed and recommended by PCI Ecology. doi: 10.1101/442160
[4] Lee, K. A., & Klasing, K. C. (2004). A role for immunology in invasion biology. Trends in Ecology & Evolution, 19(10), 523-529. doi: 10.1016/j.tree.2004.07.012

How do you estimate the biodiversity of a whole community, or the distribution of abundances and ranges of its species, from presence/absence data in scattered samples?
It all starts with the collector's dilemma: if you double the number of samples, you will not get double the number of species, since you will find many of the same common species, and only a few new rare ones.
This non-additivity has prompted many ecologists to study the Species-Area Relationship. A common theoretical approach has been to connect this spatial pattern to the overall distribution of how common or rare a species can be. At least since Fisher's celebrated log-series [1], ecologists have been trying to, first, infer the shape of the Species Abundance Distribution, and then, use it to predict how many species should be found in a given area or a given number of samples. This has found many applications, from microbial communities to tropical forests, from estimating the number of yet-unknown species to predicting how much biodiversity may be lost if a fraction of the habitat is removed.
In this elegant work, Tovo et al. [2] propose a method that starts only from presence/absence data over a number of samples, and provides the community's diversity, as well as its abundance and range size distributions. This method is simple, analytically explicit, and accurate: the authors test it on the classic Pasoh and Barro Colorado Island tropical forest datasets, and on simulated data. They make a very laudable effort in both explaining its theoretical underpinnings, and proposing a straightforward step-by-step guide to applying it to data.
The core of Tovo et al's method is a simple property: the scale invariance of the Negative Binomial (NB) distribution. Subsampling from a NB gives another NB, where a single parameter has changed. Therefore, if the Species Abundance Distribution is close enough to some NB (which is flexible enough to accommodate all the data here), we can estimate how this parameter changes when going from (1) a single sample to (2) all the available samples, and from there, extrapolate to (3) the entire community.
This principle was first applied by the authors in a previous study [3] that required abundance data in the samples, rather than just presence/absence. Given that binary occurrence data is far more available in a variety of empirical settings, this extension is worthwhile (including its new predictions on range size distributions), and it deserves to be widely known and tested.

ADDITIONAL COMMENTS

1) To explain the novelty of the authors' contribution, it is useful to look at competing techniques.
Some ""parametric"" approaches try to infer the whole-community Species Abundance Distribution (SAD) by guessing its functional form (Gaussian, power-law, log-series...) and fitting its parameters from sampled data. The issue is that this distribution shape may not remain in the same family as we increase the sampling effort or area, so the regression problem may not be well-defined. This is where the Negative Binomial's scale invariance is useful.
Other ""non-parametric"" approaches have renounced guessing the whole SAD: they simply try to approximate of its tail of rare species, by looking at how many species are found in only one (or a few) samples. From this, they derive an estimate of biodiversity that is agnostic to the rest of the SAD. Tovo et al. [2] show the issue with these approaches: they extrapolate from the properties of individual samples to the whole community, but do not properly account for the bias introduced by the amount of sampling (the intermediate scale (2) in the summary above).

2) The main condition for all such approaches to work is well-mixedness: each sample should be sufficiently like a lot drawn from the same skewed lottery. As long as that condition applies, finding the best approach is a theoretical matter of probabilities and combinatorics that may, in time, be given a definite answer.
The authors also show that ""well-mixed"" is not as restrictive as it sounds: the method works both on real data (which is never perfectly mixed) and on simulations where species are even more spatially clustered than the empirical data. In addition, the Negative Binomial's scale invariance entails that, if it works well enough at some spatial scale, it will also work at all higher scales (until one reaches the edges of the sufficiently-well-mixed community)

3) One may ask: why the Negative Binomial as a Species Abundance Distribution?
If one wishes for some dynamical explanation, the Negative Binomial can be derived from neutral birth and death process with immigration, as shown by the authors in [3]. But to be applied to data, it should only be able to approximate the empirical distribution well enough (at all relevant scales). Depending on one's taste, this type of probabilistic approaches can be interpreted as:
- purely phenomenological, describing only the observational process of sampling from an existing state of affairs, not the ecological processes that gave rise to that state.
- a null model, from which everything in practice is expected to deviate to some extent.
- or a way to capture the statistical forces that tend to induce stable relationships between different patterns (as long as no ecological process opposes them strongly enough).

References

[1] Fisher, R. A., Corbet, A. S., & Williams, C. B. (1943). The relation between the number of species and the number of individuals in a random sample of an animal population. The Journal of Animal Ecology, 42-58. doi: 10.2307/1411
[2] Tovo, A., Formentin, M., Suweis, S., Stivanello, S., Azaele, S., & Maritan, A. (2019). Inferring macro-ecological patterns from local species' occurrences. bioRxiv, 387456, ver. 2 peer-reviewed and recommended by PCI Ecol. doi: 10.1101/387456
[3] Tovo, A., Suweis, S., Formentin, M., Favretti, M., Volkov, I., Banavar, J. R., Azaele, S., & Maritan, A. (2017). Upscaling species richness and abundances in tropical forests. Science Advances, 3(10), e1701438. doi: 10.1126/sciadv.1701438

Environmental Risk Assessment (ERA) is a strategic conceptual framework to characterize the nature and magnitude of risks, to humans and biodiversity, of the release of chemical contaminants in the environment. Several measures have been suggested to enhance the science and application of ERA, including the identification and acknowledgment of uncertainties that potentially influence the outcome of risk assessments, and the appropriate consideration of temporal scale and its linkage to assessment endpoints [1].
Baudrot & Charles [2] proposed to approach these questions by coupling toxicokinetics-toxicodynamics models, which describe the time-course of processes leading to the adverse effects of a toxicant, with Bayesian inference. TKTD models separate processes influencing an organismal internal exposure (´toxicokinetics´, i.e., the uptake, bioaccumulation, distribution, biotransformation and elimination of a toxicant) from processes leading to adverse effects and ultimately its death (´toxicodynamics´) [3]. Although species and substance specific, the mechanistic nature of TKTD models facilitates the comparison of different toxicants, species, life stages, environmental conditions and endpoints [4].
Baudrot & Charles [2] investigated the use of a Bayesian framework to assess the uncertainties surrounding the calibration of General Unified Threshold Models of Survival (a category of TKTD) with data from standard toxicity tests, and their propagation to predictions of regulatory toxicity endpoints such as LC(x,t) [the lethal concentration affecting any x% of the population at any given exposure duration of time t] and MF(x,t) [an exposure multiplication factor leading to any x% effect reduction due to the contaminant at any time t].
Once calibrated with empirical data, GUTS models were used to explore individual survival over time, and under untested exposure conditions. Lethal concentrations displayed a strong curvilinear decline with time of exposure. For a given total amount of contaminant, pulses separated by short time intervals yielded higher mortality than pulses separated by long time intervals, as did few pulses of high amplitude when compared to multiple pulses of low amplitude. The response to a pulsed contaminant exposure was strongly influenced by contaminant depuration times. These findings highlight one important contribution of TKTD modelling in ecotoxicology: they represent just a few of the hundreds of exposure scenarios that could be mathematically explored, and that would be unfeasible or even unethical to conduct experimentally.
GUTS models were also used for interpolations or extrapolations of assessment endpoints, and their marginal distributions. A case in point is the incipient lethal concentration. The responses of model organisms to contaminants in standard toxicity tests are typically assessed at fixed times of exposure (e.g. 24h or 48h in the Daphnia magna acute toxicity test). However, because lethal concentrations are strongly time-dependent, it has been suggested that a more meaningful endpoint would be the incipient (i.e. asymptotic) lethal concentration when time of exposure increases to infinity. The authors present a mathematical solution for calculating the marginal distribution of such incipient lethal concentration, thereby providing both more relevant information and a way of comparing experiments, compounds or species tested for different periods of time.
Uncertainties were found to change drastically with time of exposure, being maximal at extreme values of x for both LC(x,t) and MF(x,t). In practice this means that assessment endpoints estimated when the effects of the contaminant are weak (such as LC10, the contaminant concentration resulting in the mortality of 10% of the experimental population), a commonly used assessment value in ERA, are prone to be highly variable.
The authors end with recommendations for improved experimental design, including (i) using assessment endpoints at intermediate values of x (e.g., LC50 instead of LC10) (ii) prolonging exposure and recording mortality over the course of the experiment (iii) experimenting one or few peaks of high amplitude close to each other when assessing pulsed exposure. Whereas these recommendations are not that different from current practices, they are based on a more coherent mechanistic grounding.
Overall, this and other contributions from Charles, Baudrot and their research group contribute to turn TKTD models into a real tool for Environmental Risk Assessment. Further enhancement of ERA´s science and application could be achieved by extending the use of TKTD models to sublethal rather than lethal effects, and to chronic rather than acute exposure, as these are more controversial issues in decision-making regarding contaminated sites.

References

[1] Dale, V. H., Biddinger, G. R., Newman, M. C., Oris, J. T., Suter, G. W., Thompson, T., ... & Chapman, P. M. (2008). Enhancing the ecological risk assessment process. Integrated environmental assessment and management, 4(3), 306-313. doi: 10.1897/IEAM_2007-066.1
[2] Baudrot, V., & Charles, S. (2018). Recommendations to address uncertainties in environmental risk assessment using toxicokinetics-toxicodynamics models. bioRxiv, 356469, ver. 3 peer-reviewed and recommended by PCI Ecol. doi: 10.1101/356469
[3] EFSA Panel on Plant Protection Products and their Residues (PPR), Ockleford, C., Adriaanse, P., Berny, P., Brock, T., Duquesne, S., ... & Kuhl, T. (2018). Scientific Opinion on the state of the art of Toxicokinetic/Toxicodynamic (TKTD) effect models for regulatory risk assessment of pesticides for aquatic organisms. EFSA Journal, 16(8), e05377. doi: 10.2903/j.efsa.2018.5377
[4] Jager, T., Albert, C., Preuss, T. G., & Ashauer, R. (2011). General unified threshold model of survival-a toxicokinetic-toxicodynamic framework for ecotoxicology. Environmental science & technology, 45(7), 2529-2540. doi: 10.1021/es103092a

There are several important topics in the study of infectious diseases that have not been well explored due to technical difficulties. One such topic is pursued by Alizon et al. in “Modelling coinfections to detect within-host interactions from genotype combination prevalences” [1]. Both theory and several important examples have demonstrated that interactions among co-infecting strains can have outsized impacts on disease outcomes, transmission dynamics, and epidemiology. Unfortunately, empirical data on pathogen interactions and their outcomes is often correlational making results difficult to decipher.
The analytical framework developed by Alizon et al. [1] infers the presence and strength of pathogen interactions through their impact on transmission dynamics using a novel application of Approximate Bayesian Computation (ABC)-regression to epidemiological data. Traditional analytic approaches identify pathogen interactions when the observed distribution of pathogens among hosts differ from ‘neutral’ expectations. However, deviations from this expectation are not only a result of inter-strain interactions but can be caused by many ecological interactions, such as heterogeneity in host contact networks. To overcome this difficulty, Alizon et al [1] develop an analytical framework that incorporates explicit epidemiological models to allow inference of interactions among strains of Human Papillomaviruses (HPV) even with other ecological interactions that impact the distribution of strains among hosts. Alizon et al also demonstrate that using more of the available data, including the specific combination of strains present in hosts and knowledge of the connectivity of the hosts (i.e., super-spreaders), leads to more accurate inferences of the strength and direction of within-host interactions among coinfecting strains. This method successfully identified data generated from models with high and moderate inter-strain interaction intensity when the host population was homogeneous and was only slightly less successful when the host population was heterogeneous (super-spreaders present). By comparison, some previously published analytical methods could identify only some inter-strain interactions in datasets generated from models with homogeneous host populations, but host heterogeneity obscured these interactions.
This manuscript makes seamless connections between basic viral biology and its epidemiological consequences by tying them together with realistic models, illustrating the fundamental utility of biological modeling. This analytical framework provides crucial tools for experimentalists, facilitating collaborations with theoreticians to better understand the epidemiological consequences of co-infections. In addition, the method is simple enough to be applied by a broad base of experimentalists to the many pathogens where co-infections are common. Thus, this paper has the potential to impact several research fields and public health practice. Those attempting to apply this method should note the potential limitations noted by the authors. For example, it is not designed to detect the mechanisms of inter-strain interactions (there is no within host component of the models) but to identify the existence of interactions through patterns indicative of these interactions while ruling out other sources that could cause the pattern. This approach is likely to be most accurate when strain identification within hosts is precise and unbiased - which is unlikely in many systems where samples are taken only from symptomatic cases and strain detection is not sufficiently sensitive – and when host contact networks can be reasonably estimated. Importantly, a priori knowledge of the set of possible epidemiological models is needed for accurate parameter estimates, which may be true for several prominent pathogens, but not be so for many other pathogens and symbionts. We look forward to future extensions of this framework where this restriction is relaxed. Alizon et al. [1] have provided a framework that will facilitate theoretical and empirical work on the impact of coinfections on infectious disease and should shape future public health data collection standards.

References

[1] Alizon, S., Murall, C.L., Saulnier, E., & Sofonea, M.T. (2018). Detecting within-host interactions using genotype combination prevalence data. bioRxiv, 256586, ver. 3 peer-reviewed and recommended by PCI Ecology. doi: 10.1101/256586

Are biological invasions driven by a few pioneers, running ahead of their conspecifics? Or are these pioneers constantly being caught up by, and folded into, the larger flux of propagules from the established populations behind them?
In ecology and beyond, these two scenarios are known as "pulled" and "pushed" fronts, and they come with different expectations. In a pushed front, invasion speed is not just a matter of how good individuals are at dispersing and settling new locations. It becomes a collective, density-dependent property of population fluxes. And in particular, it can depend on the equilibrium abundance of the established populations inside the range, i.e. the species’ carrying capacity K, factoring in its abiotic environment and biotic interactions.
This realization is especially important because it can flip around our expectations about which species expand fast, and how to manage them. We tend to think of initial colonization and long-term abundance as two independent axes of variation among species or indeed as two ends of a spectrum, in the classic competition-colonization tradeoff [1]. When both play into invasion speed, good dispersers might not outrun good competitors. This is useful knowledge, whether we want to contain an invasion or secure a reintroduction.
In their study "When higher carrying capacities lead to faster propagation", Haond et al [2] combine mathematical analysis, Individual-Based simulations and experiments to show that various mechanisms can cause pushed fronts, whose speed increases with the carrying capacity K of the species. Rather than focus on one particular angle, the authors endeavor to demonstrate that this qualitative effect appears again and again in a variety of settings.
It is perhaps surprising that this notable and general connection between K and invasion speed has managed to garner so little fame in ecology. A large fraction of the literature employs the venerable Fisher-KPP reaction-diffusion model, which combines local logistic growth with linear diffusion in space. This model has prompted both considerable mathematical developments [3] and many applications to modelling real invasions [4]. But it only allows pulled fronts, driven by the small populations at the edge of a species range, with a speed that depends only on their initial growth rate r.
This classic setup is, however, singular in many ways. Haond et al [2] use it as a null model, and introduce three mechanisms or factors that each ensure a role of K in invasion speed, while giving less importance to the pioneers at the border.
Two factors, the Allee effect and demographic stochasticity, make small edge populations slower to grow or less likely to survive. These two factors are studied theoretically, and to make their claims stronger, the authors stack the deck against K. When generalizing equations or simulations beyond the null case, it is easy to obtain functional forms where the parameter K does not only play the role of equilibrium carrying capacity, but also affects dynamical properties such as the maximum or mean growth rate. In that case, it can trivially change the propagation speed, without it meaning anything about the role of established populations behind the front. Haond et al [2] avoid this pitfall by disentangling these effects, at the cost of slightly more peculiar expressions, and show that varying essentially nothing but the carrying capacity can still impact the speed of the invasion front.
The third factor, density-dependent dispersal, makes small populations less prone to disperse. It is well established empirically and theoretically that various biological mechanisms, from collective organization to behavioral switches, can prompt organisms in denser populations to disperse more, e.g. in such a way as to escape competition [5]. The authors demonstrate how this effect induces a link between carrying capacity and invasion speed, both theoretically and in a dispersal experiment on the parasitoid wasp, Trichogramma chilonis.
Overall, this study carries a simple and clear message, supported by valuable contributions from different angles. Although some sections are clearly written for the theoretical ecology crowd, this article has something for everyone, from the stray physicist to the open-minded manager. The collaboration between theoreticians and experimentalists, while not central, is worthy of note. Because the narrative of this study is the variety of mechanisms that can lead to the same qualitative effect, the inclusion of various approaches is not a gimmick, but helps drive home its main message. The work is fairly self-contained, although one could always wish for further developments, especially in the direction of more quantitative testing of these mechanisms.
In conclusion, Haond et al [2] effectively convey the widely relevant message that, for some species, invading is not just about the destination, it is about the many offspring one makes along the way.

References

[1] Levins, R., & Culver, D. (1971). Regional Coexistence of Species and Competition between Rare Species. Proceedings of the National Academy of Sciences, 68(6), 1246–1248. doi: 10.1073/pnas.68.6.1246
[2] Haond, M., Morel-Journel, T., Lombaert, E., Vercken, E., Mailleret, L., & Roques, L. (2018). When higher carrying capacities lead to faster propagation. BioRxiv, 307322. doi: 10.1101/307322
[3] Crooks, E. C. M., Dancer, E. N., Hilhorst, D., Mimura, M., & Ninomiya, H. (2004). Spatial segregation limit of a competition-diffusion system with Dirichlet boundary conditions. Nonlinear Analysis: Real World Applications, 5(4), 645–665. doi: 10.1016/j.nonrwa.2004.01.004
[4] Shigesada, N., & Kawasaki, K. (1997). Biological Invasions: Theory and Practice. Oxford University Press, UK.
[5] Matthysen, E. (2005). Density-dependent dispersal in birds and mammals. Ecography, 28(3), 403–416. doi: 10.1111/j.0906-7590.2005.04073.x