As far back as Clements [1] and Gleason [2], the historical schism between deterministic and stochastic perspectives has divided ecologists. Deterministic theories tend to emphasize niche-based processes such as environmental filtering and species interactions as the main drivers of species distribution in nature, while stochastic theories mainly focus on chance colonization, random extinctions and ecological drift [3]. Although the old days when ecologists were fighting fiercely over null models and their adequacy to capture niche-based processes is over [4], the ghost of that debate between deterministic and stochastic perspectives came back to haunt ecologists in the form of the ‘environment versus space’ debate with the development of metacommunity theory [5]. While interest in that question led to meaningful syntheses of metacommunity dynamics in natural systems [6], it also illustrated how context-dependant the answer was [7]. One of the next frontiers in metacommunity ecology is to identify the underlying drivers of this observed context-dependency in the relative importance of ecological processus [7, 8].
Reflecting on seminal work by Robert MacArthur emphasizing different processes at different spatial scales [9, 10] (the so-called ‘MacArthur paradox’), Chase and Myers proposed in 2011 that a key in solving the deterministic versus stochastic debate was probably to turn our attention to how the relative importance of local processes changes across spatial scales [3]. Scale-dependance is a well-acknowledged challenge in ecology, hampering empirical syntheses and comparisons between studies [11-14]. Embracing the scale-dependance of ecological processes would not only lead to stronger syntheses and consolidation of current knowledge, it could also help resolve many current debates or apparent contradictions [11, 15, 16].
The timely study by Siqueira et al. [17] fits well within this historical context by exploring the relative importance of ecological drift and selection across a gradient of community size (number of individuals in a given community). More specifically, they tested the hypothesis that small communities are more dissimilar among each other because of ecological drift compared to large communities, which are mainly structured by niche selection [17]. That smaller populations or communities should be more affected by drift is a mathematical given [18], but the main questions are i) for a given community size how important is ecological drift relative to other processes, and ii) how small does a community have to be before random assembly dominates? The authors answer these questions using an extensive stream dataset with a community size gradient sampled from 200 streams in two climatic regions (Brazil and Finland). Combining linear models with recent null model approaches to measure deviations from random expectations [19], they show that, as expected based on theory and recent experimental work, smaller communities tend to have higher β-diversity, and that those β-diversity patterns could not be distinguished from random assembly processes [17]. Spatial turnover among larger communities is mainly driven by niche-based processes related to species sorting or dispersal dynamics [17]. Given the current environmental context, with many anthropogenic perturbations leading to reduced community size, it is legitimate to wonder, as the authors do, whether we are moving toward a more stochastic and thus less predictable world with obvious implications for the conservation of biodiversity [17].
The real strength of the study by Siqueira et al. [17], in my opinion, is in the inclusion of stream data from boreal and tropical regions. Interestingly and most importantly, the largest communities in the tropical streams are as large as the smallest communities in the boreal streams. This is where the study should really have us reflect on the notions of context-dependency in observed patterns because the negative relationship between community size and β-diversity was only observed in the tropical streams, but not in the boreal streams [17]. This interesting nonlinearity in the response means that a study that would have investigated the drift versus niche-based question only in Finland would have found very different results from the same study in Brazil. Only by integrating such a large scale gradient of community sizes together could the authors show the actual shape of the relationship, which is the first step toward building a comprehensive synthesis on a debate that has challenged ecologists for almost a century.

References

[1] Clements, F. E. (1936). Nature and structure of the climax. Journal of ecology, 24(1), 252-284. doi: 10.2307/2256278
[2] Gleason, H. A. (1917). The structure and development of the plant association. Bulletin of the Torrey Botanical Club, 44(10), 463-481. doi: 10.2307/2479596
[3] Chase, J. M., and Myers, J. A. (2011). Disentangling the importance of ecological niches from stochastic processes across scales. Philosophical transactions of the Royal Society B: Biological sciences, 366(1576), 2351-2363. doi: 10.1098/rstb.2011.0063
[4] Diamond, J. M., and Gilpin, M. E. (1982). Examination of the “null” model of Connor and Simberloff for species co-occurrences on islands. Oecologia, 52(1), 64-74. doi: 10.1007/BF00349013
[5] Leibold M. A., et al. (2004). The metacommunity concept: a framework for multi‐scale community ecology. Ecology letters, 7(7), 601-613. doi: 10.1111/j.1461-0248.2004.00608.x
[6] Cottenie, K. (2005). Integrating environmental and spatial processes in ecological community dynamics. Ecology letters, 8(11), 1175-1182. doi: 10.1111/j.1461-0248.2005.00820.x
[7] Leibold, M. A. and Chase, J. M. (2018). Metacommunity Ecology. Monographs in Population Biology, vol. 59. Princeton University Press. [8] Vellend, M. (2010). Conceptual synthesis in community ecology. The Quarterly review of biology, 85(2), 183-206. doi: 10.1086/652373
[9] MacArthur, R. H., and Wilson, E. O. (1963). An equilibrium theory of insular zoogeography. Evolution, 17(4), 373-387. doi: 10.1111/j.1558-5646.1963.tb03295.x
[10] MacArthur, R. H., and Levins, R. (1967). The limiting similarity, convergence, and divergence of coexisting species. The American Naturalist, 101(921), 377-385. doi: 10.1086/282505
[11] Viana, D. S., and Chase, J. M. (2019). Spatial scale modulates the inference of metacommunity assembly processes. Ecology, 100(2), e02576. doi: 10.1002/ecy.2576
[12] Chave, J. (2013). The problem of pattern and scale in ecology: what have we learned in 20 years?. Ecology letters, 16, 4-16. doi: 10.1111/ele.12048
[13] Patrick, C. J., and Yuan, L. L. (2019). The challenges that spatial context present for synthesizing community ecology across scales. Oikos, 128(3), 297-308. doi: 10.1111/oik.05802
[14] Chase, J. M., and Knight, T. M. (2013). Scale‐dependent effect sizes of ecological drivers on biodiversity: why standardised sampling is not enough. Ecology letters, 16, 17-26. doi: 10.1111/ele.12112
[15] Horváth, Z., Ptacnik, R., Vad, C. F., and Chase, J. M. (2019). Habitat loss over six decades accelerates regional and local biodiversity loss via changing landscape connectance. Ecology letters, 22(6), 1019-1027. doi: 10.1111/ele.13260
[16] Chase, J. M, Gooriah, L., May, F., Ryberg, W. A, Schuler, M. S, Craven, D., and Knight, T. M. (2019). A framework for disentangling ecological mechanisms underlying the island species–area relationship. Frontiers of Biogeography, 11(1). doi: 10.21425/F5FBG40844.
[17] Siqueira T., Saito V. S., Bini L. M., Melo A. S., Petsch D. K. , Landeiro V. L., Tolonen K. T., Jyrkänkallio-Mikkola J., Soininen J. and Heino J. (2019). Community size affects the signals of ecological drift and niche selection on biodiversity. bioRxiv 515098, ver. 4 peer-reviewed and recommended by PCI Ecology. doi: 10.1101/515098
[18] Hastings A., Gross L. J. eds. (2012). Encyclopedia of theoretical ecology (University of California Press, Berkeley).
[19] Chase, J. M., Kraft, N. J., Smith, K. G., Vellend, M., and Inouye, B. D. (2011). Using null models to disentangle variation in community dissimilarity from variation in α‐diversity. Ecosphere, 2(2), 1-11. doi: 10.1890/ES10-00117.1

Worldwide, city expansion is happening at a fast rate and at the same time, urbanists are more and more required to make place for biodiversity. Choices have to be made regarding the area and spatial arrangement of suitable spaces for non-human living organisms, that will favor the long-term survival of their populations. To guide those choices, it is necessary to understand the mechanisms driving the effects of land management on biodiversity.

Research results on the effects of urbanization on genetic diversity have been very diverse, with studies showing higher genetic diversity in rural than in urban populations (e.g. Delaney et al. 2010), the contrary (e.g. Miles et al. 2018) or no difference (e.g. Schoville et al. 2013). The same is true for studies investigating genetic differentiation. The reasons for these differences probably lie in the relative intensities of gene flow and genetic drift in each case study, which are hard to disentangle and quantify in empirical datasets.

In their paper, Savary et al. (2024) used an elegant and powerful simulation approach to better understand the diversity of observed patterns and investigate the effects of dispersal limitation on genetic patterns (diversity and differentiation). Their simulations involved the landscapes of 325 real European cities, each under three different scenarios mimicking 3 virtual urban tolerant species with different abilities to move within cities while genetic drift intensity was held constant across scenarios. The cities were chosen so that the proportion of artificial areas was held constant (20%) but their location and shape varied. This design allowed the authors to investigate the effect of connectivity and spatial configuration of habitat on the genetic responses to spatial variations in dispersal in cities. 

The main results of this simulation study demonstrate that variations in dispersal spatial patterns, for a given level of genetic drift, trigger variations in genetic patterns. Genetic diversity was lower and genetic differentiation was larger when species had more difficulties to move through the more hostile components of the urban environment. The increase of the relative importance of drift over gene flow when dispersal was spatially more constrained was visible through the associated disappearance of the pattern of isolation by resistance. Forest patches (usually located at the periphery of the cities) usually exhibited larger genetic diversity and were less differentiated than urban green spaces. But interestingly, the presence of habitat patches at the interface between forest and urban green spaces lowered those differences through the promotion of gene flow. 

One other noticeable result, from a landscape genetic method point of view, is the fact that there might be a limit to the detection of barriers to genetic clusters through clustering analyses because of the increased relative effect of genetic drift. This result needs to be confirmed, though, as genetic structure has only been investigated with a recent approach based on spatial graphs. It would be interesting to also analyze those results with the usual Bayesian genetic clustering approaches. 

Overall, this study addresses an important scientific question about the mechanisms explaining the diversity of observed genetic patterns in cities. But it also provides timely cues for connectivity conservation and restoration applied to cities.  
 
References

Delaney, K. S., Riley, S. P., and Fisher, R. N. (2010). A rapid, strong, and convergent genetic response to urban habitat fragmentation in four divergent and widespread vertebrates. PLoS ONE, 5(9):e12767.
https://doi.org/10.1371/journal.pone.0012767
 
Miles, L. S., Dyer, R. J., and Verrelli, B. C. (2018). Urban hubs of connectivity: Contrasting patterns of gene flow within and among cities in the western black widow spider. Proceedings of the Royal Society B, 285(1884):20181224. https://doi.org/10.1098/rspb.2018.1224
 
Savary P., Tannier C., Foltête J.-C., Bourgeois M., Vuidel G., Khimoun A., Moal H., and Garnier S. (2024). How does dispersal shape the genetic patterns of animal populations in European cities? A simulation approach. EcoEvoRxiv, ver. 3 peer-reviewed and recommended by Peer Community in Ecology. https://doi.org/10.32942/X2JS41.
 
Schoville, S. D., Widmer, I., Deschamps-Cottin, M., and Manel, S. (2013). Morphological clines and weak drift along an urbanization gradient in the butterfly, Pieris rapae. PLoS ONE, 8(12):e83095.
https://doi.org/10.1371/journal.pone.0083095

High-throughput sequencing-based techniques such as DNA metabarcoding are increasingly advocated as providing numerous benefits over morphology‐based identifications for biodiversity inventories and ecosystem biomonitoring [1]. These benefits are particularly apparent for highly-diversified and/or hardly accessible aquatic and marine environments, where simple water or sediment samples could already produce acceptably accurate biodiversity estimates based on the environmental DNA present in the samples [2,3]. However, sequence-based characterization of biodiversity comes with its own challenges. A major one resides in the capacity to disentangle true biological diversity (be it taxonomic or genetic) from artefactual diversity generated by sequence-errors accumulation during PCR and sequencing processes, or from the amplification of non-target genes (i.e. pseudo-genes). On one hand, the stringent elimination of sequence variants might lead to biodiversity underestimation through the removal of true species, or the clustering of closely-related ones. On the other hand, a more permissive sequence filtering bears the risks of biodiversity inflation. Recent studies have outlined an excellent methodological framework for addressing this issue by proposing bioinformatic tools that allow the amplicon-specific error-correction as alternative or as complement to the more arbitrary approach of clustering into Molecular Taxonomic Units (MOTUs) based on sequence dissimilarity [4,5]. But to date, the relevance of amplicon-specific error-correction tools has been demonstrated only for a limited set of taxonomic groups and gene markers.
The study of Brandt et al. [6] successfully builds upon existing methodological frameworks for filling this gap in current literature. By proposing a bioinformatic pipeline combining Amplicon Sequence Variants (ASV) curation with MOTU clustering and additional post-clustering curation, the authors show that contrary to previous recommendations, ASV-based curation alone does not represent an adequate approach for DNA metabarcoding-based inventories of metazoans. Metazoans indeed, do exhibit inherently higher intra-specific and intra-individual genetic variability, necessarily leading to biased biodiversity estimates unbalanced in favor of species with higher intraspecific diversity in the absence of MOTU clustering. Interestingly, the positive effect of additional clustering showed to be dependent on the target gene region. Additional clustering had proportionally higher effect on the more polymorphic mitochondrial COI region (as compared to the 18S ribosomal gene). Thus, the major advantage of the study lies in the provision of optimal curation parameters that reflect the best possible balance between minimizing the impact of PCR/sequencing errors and the loss of true biodiversity across markers with contrasting levels of intragenomic variation. This is important as combining multiple markers is increasingly considered for improving the taxonomic coverage and resolution of data in DNA metabarcoding studies.
Another critical aspect of the study is the taxonomic assignation of curated OTUs (which is also the case for the majority of DNA metabarcoding-based biodiversity assessments). Facing the double challenge of focusing on taxonomic groups that are both highly diverse and poorly represented in public sequence reference databases, the authors failed to obtain high-resolution taxonomic assignments for several of the most closely-related species. As a result, taxa with low divergence levels were clustered as single taxonomic units, subsequently leading to underestimation of true biodiversity present. This finding adds to the argument that in order to be successful, sequence-based techniques still require the availability of comprehensive, high-quality reference databases.
Perhaps the only regret we might have with the study is the absence of mock community validation for the prokaryotes compartment. Even though the analyses of natural samples seem to suggest a positive effect of the curation pipeline, the concept of intra- versus inter-species variation in naturally occurring prokaryote communities remains at best ambiguous. Of course, constituting a representative sample of taxonomically-resolved prokaryote taxa from deep-sea habitats does not come without difficulties but has the benefit of opening opportunities for further studies on the matter.

References

[1] Porter, T. M., and Hajibabaei, M. (2018). Scaling up: A guide to high-throughput genomic approaches for biodiversity analysis. Molecular Ecology, 27(2), 313–338. doi: 10.1111/mec.14478
[2] Valentini, A., Taberlet, P., Miaud, C., Civade, R., Herder, J., Thomsen, P. F., … Dejean, T. (2016). Next-generation monitoring of aquatic biodiversity using environmental DNA metabarcoding. Molecular Ecology, 25(4), 929–942. doi: 10.1111/mec.13428
[3] Leray, M., and Knowlton, N. (2015). DNA barcoding and metabarcoding of standardized samples reveal patterns of marine benthic diversity. Proceedings of the National Academy of Sciences, 112(7), 2076–2081. doi: 10.1073/pnas.1424997112
[4] Callahan, B. J., McMurdie, P. J., and Holmes, S. P. (2017). Exact sequence variants should replace operational taxonomic units in marker-gene data analysis. The ISME Journal, 11(12), 2639–2643. doi: 10.1038/ismej.2017.119
[5] Edgar, R. C. (2016). UNOISE2: improved error-correction for Illumina 16S and ITS amplicon sequencing. BioRxiv, 081257. doi: 10.1101/081257
[6] Brandt, M. I., Trouche, B., Quintric, L., Wincker, P., Poulain, J., and Arnaud-Haond, S. (2020). A flexible pipeline combining clustering and correction tools for prokaryotic and eukaryotic metabarcoding. BioRxiv, 717355, ver. 3 peer-reviewed and recommended by PCI Ecology. doi: 10.1101/717355

Increased knowledge of what factors are determining species interactions are of major importance for our understanding of dynamics and functionality of ecological communities [1]. Currently, when ongoing temperature modifications lead to changes in species temporal and spatial limits the subject gets increasingly topical. A species phenology determines whether it thrive or survive in its environment. However, as the phenologies of different species are not necessarily equally affected by environmental changes, temporal or spatial mismatches can occur and affect the species-species interactions in the network [2] and as such the full network structure.
In this preprint by Manincor et al. [3] the authors explore the effect of phenology overlap on a large network of species interactions in calcareous grasslands in France. They analyze if and how this effect varies along a latitudinal gradient using empirical data on six plant-hoverfly networks. When comparing ecological network along gradients a well-known problem is that the network metrics is dependent on network size [4]. Therefore, instead of focusing on complete network structure the authors here focus on the factors that determine the probability of interactions and interaction frequency (number of visits). The authors use Bayesian Structural Equation Models (SEM) to link the interaction probability and number of visits to phenology overlap and species abundance. SEM is a multivariate technique that can be used to test several hypotheses and evaluate multiple causal relationships using both observed and latent variables to explain some other observed variables. The authors provide a nice description of the approach for this type of study system. In addition, the study also tests whether phenology affects network compartmentalization, by analyzing species subgroups using a latent block model (LBM) which is a clustering method particularly well-suited for weighted networks.
The authors identify phenology overlap as an important determinant of plant-pollinator interactions, but also conclude this factor alone is not sufficient to explain the species interactions. Species abundances was important for number of visits. Plant phenology drives the duration of the phenology overlap between plant and hoverflies in the studied system. This in turn influences either the probability of interaction or the expected number of visits, as well as network compartmentalization. Longer phenologies correspond to lower modularity inferring less constrained interactions, and shorter phenologies correspond to higher modularity inferring more constrained interactions.
What make this study particularly interesting is the presentation of SEMs as an innovative approach to compare networks of different sizes along environmental gradients. The authors show that these methods can be a useful tool when the aim is to understand the structure of plant-pollinator networks and data is varying in complexities. During the review process the authors carefully addressed to the comments from the two reviewers and the manuscript improved during the process. Both reviewers have expertise highly relevant for the research performed and the development of the manuscript. In my opinion this is a highly interesting and valuable piece of work both when it comes to the scientific question and the methodology. I look forward to further follow this research.

References

[1] Pascual, M., and Dunne, J. A. (Eds.). (2006). Ecological networks: linking structure to dynamics in food webs. Oxford University Press.
[2] Parmesan, C. (2007). Influences of species, latitudes and methodologies on estimates of phenological response to global warming. Global Change Biology, 13(9), 1860-1872. doi: 10.1111/j.1365-2486.2007.01404.x
[3] de Manincor, N., Hautekeete, N., Piquot, Y., Schatz, B., Vanappelghem, C. and Massol, F. (2019). Does phenology explain plant-pollinator interactions at different latitudes? An assessment of its explanatory power in plant-hoverfly networks in French calcareous grasslands. Zenodo, 2543768, ver. 4 peer-reviewed and recommended by PCI Ecology. doi: 10.5281/zenodo.2543768
[4] Staniczenko, P. P., Kopp, J. C., and Allesina, S. (2013). The ghost of nestedness in ecological networks. Nature communications, 4, 1391. doi: 10.1038/ncomms2422

Amongst the many challenges and forms of environmental change that organisms face in our era of global change, climate change is perhaps one of the most straightforward and amenable to investigation. First, measurements of day-to-day temperatures are relatively feasible and accessible, and predictions regarding the expected trends in Earth surface temperature are probably some of the most reliable we have. It appears quite clear, in particular, that beyond the overall increase in average temperature, the heat waves locally experienced by organisms in their natural habitats are bound to become more frequent, more intense, and more long-lasting [1]. Second, it is well appreciated that temperature is a major environmental factor with strong impacts on different facets of organismal development and life-history [2-4]. These impacts have reasonably clear mechanistic underpinnings, with definite connections to biochemistry, physiology, and considerations on energetics. Third, since variation in temperature is a challenge already experienced by natural populations across their current and historical ranges, it is not a completely alien form of environmental change. Therefore, we already learnt quite a lot about it in several species, and so did the species, as they may be expected to have evolved dedicated adaptive mechanisms to respond to elevated temperatures. Last, but not least, temperature is quite amenable to being manipulated as an experimental factor.
For all these reasons, experimental studies of the consequences of increased temperature hit some of a sweetspot and are a source of very nice research, in many different organisms. The work by Leicht and Seppala [5] complements a sequence of earlier studies by this group, using the freshwater snail Lymnaea stagnalis as their model system [6-7].
In the present study, the authors investigate how a heat wave (a period of abnormally elevated temperature, here 25°C versus a normal 15°C) may have indirect effects on the next generation, through maternal effects. They question whether such indirect effects exist, and if they exist, how they compare, in terms of effect size, with the (more straightforward) direct effects observed in individuals that directly experience a heat wave. Transgenerational effects are well-known to occur following periods of physiological stress, and might thus have non negligible contributions to the overall effect of warming.
In this freshwater snail, heat has very strong direct effects: mortality increases at high temperature, but survivors grow much bigger, with a greater propensity to lay eggs and a (spectacular) three-fold increase in the number of eggs laid [6]. Considering that, it is easy to consider that transgenerational effects should be small game. And indeed, the present study also observes the big and obvious direct effects of elevated temperature: higher mortality, but greater propensity to oviposit. However, it was also found that the eggs were smaller if from mothers exposed to high temperature, with a correspondingly smaller size of hatchlings. This suggests that a heat wave causes the snails to lay more eggs, but smaller ones, reminiscent of a size-number trade-off. Unfortunately, clutch size could not be measured in this experiment, so this cannot be investigated any further. For this trait, the indirect effect may indeed be regarded as small game : eggs and hatchlings were about 15 % smaller, an effect size pretty small compared to the mammoth direct positive effect of temperature on shell length (see Figure 4 ; and also [6]). The same is true for developmental time (Figure 3).
However, for some traits the story was different. In particular, it was found that the (smaller) eggs produced from heated mothers were more likely to hatch by almost 10% (Figure 2). Here the indirect effect not only goes against the direct effect (hatching rate is lower at high temperature), but it also has similar effect size. As a consequence, taking into account both the indirect and direct effects, hatching success is essentially the same at 15°C and 25°C (Figure 2). Survival also had comparable effect sizes for direct and indirect effects. Indeed, survival was reduced by about 20% regardless of whom endured the heat stress (the focal individual or her mother; Figure 4). Interestingly, the direct and indirect effects were not quite cumulative: if a mother experienced a heat wave, heating up the offspring did not do much more damage, as though the offspring were ‘adapted’ to the warmer conditions (but keep in mind that, surprisingly, the authors’ stats did not find a significant interaction; Table 2).
At the end of the day, even though at first heat seems a relatively simple and understandable component of environmental change, this study shows how varied its effects can be effects on different components of individual fitness. The overall impact most likely is a mix of direct and indirect effects, of shifts along allocation trade-offs, and of maladaptive and adaptive responses, whose overall ecological significance is not so easy to grasp. That said, this study shows that direct and indirect (maternal) effects can sometimes go against one another and have similar intensities. Indirect effects should therefore not be overlooked in this kind of studies. It also gives a hint of what an interesting challenge it is to understand the adaptive or maladaptive nature of organism responses to elevated temperatures, and to evaluate their ultimate fitness consequences.

References

[1] Meehl, G. A., & Tebaldi, C. (2004). More intense, more frequent, and longer lasting heat waves in the 21st century. Science (New York, N.Y.), 305(5686), 994–997. doi: 10.1126/science.1098704
[2] Adamo, S. A., & Lovett, M. M. E. (2011). Some like it hot: the effects of climate change on reproduction, immune function and disease resistance in the cricket Gryllus texensis. The Journal of Experimental Biology, 214(Pt 12), 1997–2004. doi: 10.1242/jeb.056531
[3] Deutsch, C. A., Tewksbury, J. J., Tigchelaar, M., Battisti, D. S., Merrill, S. C., Huey, R. B., & Naylor, R. L. (2018). Increase in crop losses to insect pests in a warming climate. Science (New York, N.Y.), 361(6405), 916–919. doi: 10.1126/science.aat3466
[4] Sentis, A., Hemptinne, J.-L., & Brodeur, J. (2013). Effects of simulated heat waves on an experimental plant–herbivore–predator food chain. Global Change Biology, 19(3), 833–842. doi: 10.1111/gcb.12094
[5] Leicht, K., & Seppälä, O. (2019). Direct and transgenerational effects of an experimental heat wave on early life stages in a freshwater snail. BioRxiv, 449777, ver. 4 peer-reviewed and recommended by PCI Ecology. doi: 10.1101/449777
[6] Leicht, K., Seppälä, K., & Seppälä, O. (2017). Potential for adaptation to climate change: family-level variation in fitness-related traits and their responses to heat waves in a snail population. BMC Evolutionary Biology, 17(1), 140. doi: 10.1186/s12862-017-0988-x
[7] Leicht, K., Jokela, J., & Seppälä, O. (2013). An experimental heat wave changes immune defense and life history traits in a freshwater snail. Ecology and Evolution, 3(15), 4861–4871. doi: 10.1002/ece3.874