Coral reefs and their associated fauna are largely impacted by ongoing climate change. Unravelling species responses to past climatic variations might provide clues on the consequence of ongoing changes. Here, we tested the relationship between changes in sea surface temperature and sea levels during the Quaternary and present-day distributions of coral reef fish species. We investigated whether species-specific responses are associated with life-history traits. We collected a database of coral reef fish distribution together with life-history traits for the Indo-Pacific Ocean. We ran species distribution models (SDMs) on 3,725 tropical reef fish species using contemporary environmental factors together with a variable describing isolation from stable coral reef areas during the Quaternary. We quantified the variance explained independently by isolation from stable areas in the SDMs and related it to a set of species traits including body size and mobility. The variance purely explained by isolation from stable coral reef areas on the distribution of extant coral reef fish species largely varied across species. We observed a triangular relationship between the contribution of isolation from stable areas in the SDMs and body size. Species, whose distribution is more associated with historical changes, occurred predominantly in the Indo-Australian archipelago, where the mean size of fish assemblages is the lowest. Our results suggest that the legacy of habitat changes of the Quaternary is still detectable in the extant distribution of many fish species, especially those with small body size and the most sedentary. Because they were the least able to colonize distant habitats in the past, fish species with smaller body size might have the most pronounced lags in tracking ongoing climate change.

Liu et al. (2018) used a virtual species approach to test the effects of outliers on species distribution models. In their simulations, they applied a threshold value over the simulated suitabilities to generate the species distributions, suggesting that using a probabilistic simulation approach would have been more complex and yield the same results. Here, we argue that using a probabilistic approach is not necessarily more complex and may significantly change results. Although the threshold approach may be justified under limited circumstances, the probabilistic approach has multiple advantages. First, it is in line with ecological theory, which largely assumes non-threshold responses. Second, it is more general, as it includes the threshold as a limiting case. Third, it allows a better separation of the relevant intervening factors that influence model performance. Therefore, we argue that the probabilistic simulation approach should be used as a general standard in virtual species studies.

Habitat modelling is increasingly relevant in biodiversity and conservation studies. A typical application is to predict potential zones of specific conservation interest. With many environmental covariates, a large number of models can he investigated but multi-model inference may become impractical. Shrinkage regression overcomes this issue by dealing with the identification and accurate estimation of effect size for prediction. In a Bayesian framework we investigated the use of a shrinkage prior, the Horseshoe, for variable selection in spatial generalized linear models (GLM). As study cases, we considered 5 datasets on small pelagic fish abundance in the Gulf of Lion (Mediterranean Sea, France) and 9 environmental inputs. We compared the predictive performances of a simple kriging model, a full spatial GLM model with independent normal priors for regression coefficients, a full spatial GLM model with a Horseshoe prior for regression coefficients and 2 zero-inflated models (spatial and non-spatial) with a Horseshoe prior. Predictive performances were evaluated by cross validation on a hold-out subset of the data: models with a Horseshoe prior performed best, and the full model with independent normal priors worst. With an increasing number of inputs, extrapolation quickly became pervasive as we tried to predict from novel combinations of covariate values. By shrinking regression coefficients with a Horseshoe prior, only one model needed to be fitted to the data in order to obtain reasonable and accurate predictions, including extrapolations.

Predictive models are central to many scientific disciplines and vital for informing management in a rapidly changing world. However, limited understanding of the accuracy and precision of models transferred to novel conditions (their 'transferability') undermines confidence in their predictions. Here, 50 experts identified priority knowledge gaps which, if filled, will most improve model transfers. These are summarized into six technical and six fundamental challenges, which underlie the combined need to intensify research on the determinants of ecological predictability, including species traits and data quality, and develop best practices for transferring models. Of high importance is the identification of a widely applicable set of transferability metrics, with appropriate tools to quantify the sources and impacts of prediction uncertainty under novel conditions.

Factors shaping the distribution of mesophotic octocorals (30-200 m depth) remain poorly understood, potentially leaving overlooked coral areas, particularly near their bathymetric and geographic distributional limits. Yet, detailed knowledge about habitat requirements is crucial for conservation of sensitive gorgonians. Here we use Ecological Niche Modelling (ENM) relating thirteen environmental predictors and a highly comprehensive presence dataset, enhanced by SCUBA diving surveys, to investigate the suitable habitat of an important structuring species, Paramuricea clavata, throughout its distribution (Mediterranean and adjacent Atlantic). Models showed that temperature (11.5-25.5 degrees C) and slope are the most important predictors carving the niche of P. clavata. Prediction throughout the full distribution (TSS 0.9) included known locations of P. clavata alongside with previously unknown or unreported sites along the coast of Portugal and Africa, including seamounts. These predictions increase the understanding of the potential distribution for the northern Mediterranean and indicate suitable hard bottom areas down to > 150 m depth. Poorly sampled habitats with predicted presence along Algeria, Alboran Sea and adjacent Atlantic coasts encourage further investigation. We propose that surveys of target areas from the predicted distribution map, together with local expert knowledge, may lead to discoveries of new P. clavata sites and identify priority conservation areas.