Using the multivariate Hawkes process to study interactions between multiple species from camera trap data.

Interspecific interactions can influence species' activity and movement patterns. In particular, species may avoid or attract each other through reactive responses in space and/or time. However, data and methods to study such reactive interactions have remained scarce and were generally limited to two interacting species. At this time, the deployment of camera traps opens new opportunities but adapted statistical techniques are still required to analyze interaction patterns with such data. We present the multivariate Hawkes process (MHP) and show how it can be used to analyze interactions between several species using camera trap data. Hawkes processes use flexible pairwise interaction functions, allowing us to consider asymmetries and variations over time when depicting reactive temporal interactions. After describing the theoretical foundations of the MHP, we outline how its framework can be used to study interspecific interactions with camera trap data. We design a simulation study to evaluate the performance of the MHP and of another existing method to infer interactions from camera trap-like data. We also use the MHP to infer reactive interactions from real camera trap data for five species from South African savannas (impala Aepyceros melampus, greater kudu Tragelaphus strepsiceros, lion Panthera leo, blue wildebeest Connochaetes taurinus and Burchell's zebra Equus quagga burchelli). The simulation study shows that the MHP can be used as a tool to benchmark other methods of interspecific interaction inference and that this model can reliably infer interactions when enough data are considered. The analysis of real data highlights evidence of predator avoidance by prey and herbivore-herbivore attraction. Lastly, we present the advantages and limits of the MHP and discuss how it can be improved to infer attraction/avoidance patterns more reliably. As camera traps are increasingly used, the multivariate Hawkes process provides a promising framework to decipher the complexity of interactions structuring ecological communities.

PhylteR: Efficient Identification of Outlier Sequences in Phylogenomic Datasets.

In phylogenomics, incongruences between gene trees, resulting from both artifactual and biological reasons, can decrease the signal-to-noise ratio and complicate species tree inference. The amount of data handled today in classical phylogenomic analyses precludes manual error detection and removal. However, a simple and efficient way to automate the identification of outliers from a collection of gene trees is still missing. Here, we present PhylteR, a method that allows rapid and accurate detection of outlier sequences in phylogenomic datasets, i.e. species from individual gene trees that do not follow the general trend. PhylteR relies on DISTATIS, an extension of multidimensional scaling to 3 dimensions to compare multiple distance matrices at once. In PhylteR, these distance matrices extracted from individual gene phylogenies represent evolutionary distances between species according to each gene. On simulated datasets, we show that PhylteR identifies outliers with more sensitivity and precision than a comparable existing method. We also show that PhylteR is not sensitive to ILS-induced incongruences, which is a desirable feature. On a biological dataset of 14,463 genes for 53 species previously assembled for Carnivora phylogenomics, we show (i) that PhylteR identifies as outliers sequences that can be considered as such by other means, and (ii) that the removal of these sequences improves the concordance between the gene trees and the species tree. Thanks to the generation of numerous graphical outputs, PhylteR also allows for the rapid and easy visual characterization of the dataset at hand, thus aiding in the precise identification of errors. PhylteR is distributed as an R package on CRAN and as containerized versions (docker and singularity).

The global spectrum of plant form and function: enhanced species-level trait dataset.

Here we provide the 'Global Spectrum of Plant Form and Function Dataset', containing species mean values for six vascular plant traits. Together, these traits -plant height, stem specific density, leaf area, leaf mass per area, leaf nitrogen content per dry mass, and diaspore (seed or spore) mass - define the primary axes of variation in plant form and function. The dataset is based on ca. 1 million trait records received via the TRY database (representing ca. 2,500 original publications) and additional unpublished data. It provides 92,159 species mean values for the six traits, covering 46,047 species. The data are complemented by higher-level taxonomic classification and six categorical traits (woodiness, growth form, succulence, adaptation to terrestrial or aquatic habitats, nutrition type and leaf type). Data quality management is based on a probabilistic approach combined with comprehensive validation against expert knowledge and external information. Intense data acquisition and thorough quality control produced the largest and, to our knowledge, most accurate compilation of empirically observed vascular plant species mean traits to date.

Mycobacterium tuberculosis genetic features associated with pulmonary tuberculosis severity.

OBJECTIVES: Mycobacterium tuberculosis (Mtb) infections result in a wide spectrum of clinical presentations but without proven Mtb genetic determinants. Herein, we hypothesized that the genetic features of Mtb clinical isolates, such as specific polymorphisms or microdiversity, may be linked to tuberculosis (TB) severity. METHODS: A total of 234 patients with pulmonary TB (including 193 drug-susceptible and 14 monoresistant cases diagnosed between 2017 and 2020 and 27 multidrug-resistant cases diagnosed between 2010 and 2020) were stratified according to TB disease severity, and Mtb genetic features were explored using whole genome sequencing, including heterologous single-nucleotide polymorphism (SNP), calling to explore microdiversity. Finally, we performed a structural equation modeling analysis to relate TB severity to Mtb genetic features. RESULTS: The clinical isolates from patients with mild TB carried mutations in genes associated with host-pathogen interaction, whereas those from patients with moderate/severe TB carried mutations associated with regulatory mechanisms. Genome-wide association study identified an SNP in the promoter of the gene coding for the virulence regulator espR, statistically associated with moderate/severe disease. Structural equation modeling and model comparisons indicated that TB severity was associated with the detection of Mtb microdiversity within clinical isolates and to the espR SNP. CONCLUSION: Taken together, these results provide a new insight to better understand TB pathophysiology and could provide a new prognosis tool for pulmonary TB severity.

Genetic and species-level biodiversity patterns are linked by demography and ecological opportunity.

The processes that give rise to species richness gradients are not well understood, but may be linked to resource-based limits on the number of species a region can support. Ecological limits placed on regional species richness should also affect population demography, suggesting that these processes could also generate genetic diversity gradients. If true, we might better understand how broad-scale biodiversity patterns are formed by identifying the common causes of genetic diversity and species richness. We develop a hypothetical framework based on the consequences of regional variation in ecological limits set by resource availability and heterogeneity to simultaneously explain spatial patterns of species richness and neutral genetic diversity. Repurposing raw genotypic data spanning 38 mammal species sampled across 801 sites in North America, we show that estimates of genome-wide genetic diversity and species richness share spatial structure. Notably, species richness hotspots tend to harbor lower levels of within-species genetic variation. A structural equation model encompassing eco-evolutionary processes related to resource availability, habitat heterogeneity, and contemporary human disturbance supports the spatial patterns we detect. These results suggest broad-scale patterns of species richness and genetic diversity could both partly be caused by intraspecific demographic and evolutionary processes acting simultaneously across species.

Are human natal sex ratio differences across the world adaptive? A test of Fisher's principle.

Fisher's principle states that natural selection favours an equal number of male and female births at the population level, unless there are sex differences in rearing costs or sex differences in mortality before the end of the period of parental investment. Sex differences in rearing costs should be more pronounced in low- than in high-resource settings. We, therefore, examined whether human development index and sex differences in child mortality contribute to the natural variation in human sex ratio at birth across the globe. As predicted by Fisher's principle, the proportion of male births increased with both increasing male-biased childhood mortality and level of development of each country. However, these relationships were absent after accounting for spatial autocorrelation in the residuals, which our inference is conditioned on. This work shows how the failure to account for residual spatial autocorrelation can lead to incorrect conclusions regarding support for predictions from sex allocation theory.

Can an herbivore affect where a top predator kills its prey by modifying woody vegetation structure?

In large mammal communities, little is known about modification of interspecific interactions through habitat structure changes. We assessed the effects of African elephants (Loxodonta africana) on features of woody habitat structure that can affect predator-prey interactions. We then explored how this can influence where African lions (Panthera leo) kill their prey. Indeed, lions are stalk-and-ambush predators and habitat structure and concealment opportunities are assumed to influence their hunting success. During 2 years, in Hwange National Park, Zimbabwe, kill sites (n = 167) of GPS-collared lions were characterized (visibility distance for large mammals, distance to a potential ambush site and presence of elephant impacts). We compared characteristics of lion kill sites with characteristics of random sites (1) at a large scale (i.e. in areas intensively used by lions, n = 418) and (2) at the microhabitat scale (i.e. in the direct surrounding available habitat, < 150 m, n = 167). Elephant-impacted sites had a slightly higher visibility and a longer distance to a potential ambush site than non-impacted sites, but these relationships were characterized by a high variability. At large scale, kill sites were characterized by higher levels of elephant impacts compared to random sites. At microhabitat scale, compared to the direct nearby available habitat, kill sites were characterized by a reduced distance to a potential ambush site. We suggest a conceptual framework whereby the relative importance of habitat features and prey abundance could change upon the scale considered.

Investigating microbial associations from sequencing survey data with co-correspondence analysis.

Microbial communities, which drive major ecosystem functions, consist of a wide range of interacting species. Understanding how microbial communities are structured and the processes underlying this is crucial to interpreting ecosystem responses to global change but is challenging as microbial interactions cannot usually be directly observed. Multiple efforts are currently focused to combine next-generation sequencing (NGS) techniques with refined statistical analysis (e.g., network analysis, multivariate analysis) to characterize the structures of microbial communities. However, most of these approaches consider a single table of sequencing data measured for several samples. Technological advances now make it possible to collect NGS data on different taxonomic groups simultaneously for the same samples, allowing us to analyse a pair of tables. Here, an analytical framework based on co-correspondence analysis (CoCA) is proposed to study the distributions, assemblages and interactions between two microbial communities. We show the ability of this approach to highlight the relationships between two microbial communities, using two data sets exhibiting various types of interactions. CoCA identified strong association patterns between autotrophic and heterotrophic microbial eukaryote assemblages, on the one hand, and between microalgae and viruses, on the other. We demonstrate also how CoCA can be used, complementary to network analysis, to reorder co-occurrence networks and thus investigate the presence of patterns in ecological networks.

Flower phenology as a disruptor of the fruiting dynamics in temperate oak species.

Many perennial plants display masting, that is, fruiting with strong interannual variations, irregular and synchronized between trees within the population. Here, we tested the hypothesis that the early flower phenology in temperate oak species promotes stochasticity into their fruiting dynamics, which could play a major role in tree reproductive success. From a large field monitoring network, we compared the pollen phenology between temperate and Mediterranean oak species. Then, focusing on temperate oak species, we explored the influence of the weather around the time of budburst and flowering on seed production, and simulated with a mechanistic model the consequences that an evolutionary shifting of flower phenology would have on fruiting dynamics. Temperate oak species release pollen earlier in the season than do Mediterranean oak species. Such early flowering in temperate oak species results in pollen often being released during unfavorable weather conditions and frequently results in reproductive failure. If pollen release were delayed as a result of natural selection, fruiting dynamics would exhibit much reduced stochastic variation. We propose that early flower phenology might be adaptive by making mast-seeding years rare and unpredictable, which would greatly help in controlling the dynamics of seed consumers.

Phenotypic plasticity in the invasive pest Drosophila suzukii: activity rhythms and gene expression in response to temperature.

Phenotypic plasticity may contribute to the invasive success of an alien species in a new environment. A highly plastic species may survive and reproduce in more diverse environments, thereby supporting establishment and colonization. We focused on plasticity in the circadian rhythm of activity, which can favour species coexistence in invasion, for the invasive species Drosophila suzukii, which is expected to be a weaker direct competitor than other Drosophila species of the resident community. We compared the circadian rhythms of the locomotor activity in adults and the expression of clock genes in response to temperature in the invasive D. suzukii and the resident Drosophila melanogaster. We showed that D. suzukii is active in a narrower range of temperatures than D. melanogaster and that the activities of the two species overlap during the day, regardless of the temperature. Both species are diurnal and exhibit rhythmic activity at dawn and dusk, with a much lower activity at dawn for D. suzukii females. Our results show that the timeless and clock genes are good candidates to explain the plastic response that is observed in relation to temperature. Overall, our results suggest that thermal phenotypic plasticity in D. suzukii activity is not sufficient to explain the invasive success of D. suzukii and call for testing other hypotheses, such as the release of competitors and/or predators.

Diversity indices for ecological networks: a unifying framework using Hill numbers.

Describing how ecological interactions change over space and time and how they are shaped by environmental conditions is crucial to understand and predict ecosystem trajectories. However, it requires having an appropriate framework to measure network diversity locally, regionally and between samples (α-, γ- and ß-diversity). Here, we propose a unifying framework that builds on Hill numbers and accounts both for the probabilistic nature of biotic interactions and the abundances of species or groups. We emphasise the importance of analysing network diversity across different species aggregation levels (e.g. from species to trophic groups) to get a better understanding of network structure. We illustrate our framework with a simulation experiment and an empirical analysis using a global food-web database. We discuss further usages of the framework and show how it responds to recent calls on comparing ecological networks and analysing their variation across environmental gradients and time.

Integrating spatial and phylogenetic information in the fourth-corner analysis to test trait-environment relationships.

The fourth-corner analysis aims to quantify and test for relationships between species traits and site-specific environmental variables, mediated by site-specific species abundances. Since there is no common unit of observation, the significance of the relationships is tested using a double permutation procedure (site based and species based). This method implies that all species and sites are independent of each other. However, this fundamental hypothesis might be flawed because of phylogenetic relatedness between species and spatial autocorrelation in the environmental data. Here, using a simulation-based experiment, we demonstrate how the presence of spatial and phylogenetic autocorrelations can, in some circumstances, lead to inflated type I error rates, suggesting that significant associations can be misidentified. As an alternative, we propose a new randomization approach designed to avoid this issue, based on Moran's spectral randomization. In this approach, standard permutations are replaced by constrained randomizations so that the distribution of the statistic under the null hypothesis is built with additional constraints to preserve the phylogenetic and spatial structures of the observed data. The inclusion of this new randomization approach provides total control over type I error rates and should be used in real studies where spatial and phylogenetic autocorrelations often occur.

Optimizing the choice of a spatial weighting matrix in eigenvector-based methods.

Eigenvector-mapping methods such as Moran's eigenvector maps (MEM) are derived from a spatial weighting matrix (SWM) that describes the relations among a set of sampled sites. The specification of the SWM is a crucial step, but the SWM is generally chosen arbitrarily, regardless of the sampling design characteristics. Here, we compare the statistical performances of different types of SWMs (distance-based or graph-based) in contrasted realistic simulation scenarios. Then, we present an optimization method and evaluate its performances compared to the arbitrary choice of the most-widely used distance-based SWM. Results showed that the distance-based SWMs generally had lower power and accuracy than other specifications, and strongly underestimated spatial signals. The optimization method, using a correction procedure for multiple tests, had a correct type I error rate, and had higher power and accuracy than an arbitrary choice of the SWM. Nevertheless, the power decreased when too many SWMs were compared, resulting in a trade-off between the gain of accuracy and the loss of power. We advocate that future studies should optimize the choice of the SWM using a small set of appropriate candidates. R functions to implement the optimization are available in the adespatial package and are detailed in a tutorial.

Beyond neutrality: disentangling the effects of species sorting and spurious correlations in community analysis.

The methods of direct gradient analysis and variation partitioning are the most widely used frameworks to evaluate the contributions of species sorting to metacommunity structure. In many cases, however, species are also driven by spatial processes that are independent of environmental heterogeneity (e.g., neutral dynamics). As such, spatial autocorrelation can occur independently in both species (due to limited dispersal) and the environmental data, leading to spurious correlations between species distributions and the spatialized (i.e., spatially autocorrelated) environment. In these cases, the method of variation partitioning may present high Type I error rates (i.e., reject the null hypothesis more often than the pre-established critical level) and inflated estimates regarding the environmental component that is used to estimate the importance of species sorting. In this paper, we (1) demonstrate that metacommunities driven by neutral dynamics (via limited dispersal) alone or in combination with species sorting leads to inflated estimates and Type I error rates when testing for the importance of species sorting; and (2) propose a general and flexible new variation partitioning procedure to adjust for spurious contributions due to spatial autocorrelation from the environmental fraction. We used simulated metacommunity data driven by pure neutral, pure species sorting, and mixed (i.e., neutral + species sorting dynamics) processes to evaluate the performances of our new methodological framework. We also demonstrate the utility of the proposed framework with an empirical plant dataset in which we show that half of the variation initially due to the environment by the standard variation partitioning framework was due to spurious correlations.

Factors shaping community assemblages and species co-occurrence of different trophic levels.

Species assemblages are the results of various processes, including dispersion and habitat filtering. Disentangling the effects of these different processes is challenging for statistical analysis, especially when biotic interactions should be considered. In this study, we used plants (producers) and leafhoppers (phytophagous) as model organisms, and we investigated the relative importance of abiotic versus biotic factors that shape community assemblages, and we infer on their biotic interactions by applying three-step statistical analysis. We applied a novel statistical analysis, that is, multiblock Redundancy Analysis (mbRA, step 1) and showed that 51.8% and 54.1% of the overall variation in plant and leafhopper assemblages are, respectively, explained by the two multiblock models. The most important blocks of variables to explain the variations in plant and leafhopper assemblages were local topography and biotic factors. Variation partitioning analysis (step 2) showed that pure abiotic filtering and pure biotic processes were relatively less important than their combinations, suggesting that biotic relationships are strongly structured by abiotic conditions. Pairwise co-occurrence analysis (step 3) on generalist leafhoppers and the most common plants identified 40 segregated species pairs (mainly between plant species) and 16 aggregated pairs (mainly between leafhopper species). Pairwise analysis on specialist leafhoppers and potential host plants clearly revealed aggregated patterns. Plant segregation suggests heterogeneous resource availability and competitive interactions, while leafhopper aggregation suggests host feeding differentiation at the local level, different feeding microhabitats on host plants, and similar environmental requirements of the species. Using the novel mbRA, we disentangle for the first time the relative importance of more than five distinct groups of variables shaping local species communities. We highlighted the important role of abiotic processes mediated by bottom-up effects of plants on leafhopper communities. Our results revealed that in-field structure diversification and trophic interactions are the main factors causing the co-occurrence patterns observed.

A critical issue in model-based inference for studying trait-based community assembly and a solution.

Statistical testing of trait-environment association from data is a challenge as there is no common unit of observation: the trait is observed on species, the environment on sites and the mediating abundance on species-site combinations. A number of correlation-based methods, such as the community weighted trait means method (CWM), the fourth-corner correlation method and the multivariate method RLQ, have been proposed to estimate such trait-environment associations. In these methods, valid statistical testing proceeds by performing two separate resampling tests, one site-based and the other species-based and by assessing significance by the largest of the two p-values (the pmax test). Recently, regression-based methods using generalized linear models (GLM) have been proposed as a promising alternative with statistical inference via site-based resampling. We investigated the performance of this new approach along with approaches that mimicked the pmax test using GLM instead of fourth-corner. By simulation using models with additional random variation in the species response to the environment, the site-based resampling tests using GLM are shown to have severely inflated type I error, of up to 90%, when the nominal level is set as 5%. In addition, predictive modelling of such data using site-based cross-validation very often identified trait-environment interactions that had no predictive value. The problem that we identify is not an "omitted variable bias" problem as it occurs even when the additional random variation is independent of the observed trait and environment data. Instead, it is a problem of ignoring a random effect. In the same simulations, the GLM-based pmax test controlled the type I error in all models proposed so far in this context, but still gave slightly inflated error in more complex models that included both missing (but important) traits and missing (but important) environmental variables. For screening the importance of single trait-environment combinations, the fourth-corner test is shown to give almost the same results as the GLM-based tests in far less computing time.

Interspecific interference competition at the resource patch scale: do large herbivores spatially avoid elephants while accessing water?

Animals may anticipate and try to avoid, at some costs, physical encounters with other competitors. This may ultimately impact their foraging distribution and intake rates. Such cryptic interference competition is difficult to measure in the field, and extremely little is known at the interspecific level. We tested the hypothesis that smaller species avoid larger ones because of potential costs of interference competition and hence expected them to segregate from larger competitors at the scale of a resource patch. We assessed fine-scale spatial segregation patterns between three African herbivore species (zebra Equus quagga, kudu Tragelaphus strepsiceros and giraffe Giraffa camelopardalis) and a megaherbivore, the African elephant Loxodonta africana, at the scale of water resource patches in the semi-arid ecosystem of Hwange National Park, Zimbabwe. Nine waterholes were monitored every two weeks during the dry season of a drought year, and observational scans of the spatial distribution of all herbivores were performed every 15 min. We developed a methodological approach to analyse such fine-scale spatial data. Elephants increasingly used waterholes as the dry season progressed, as did the probability of co-occurrence and agonistic interaction with elephants for the three study species. All three species segregated from elephants at the beginning of the dry season, suggesting a spatial avoidance of elephants and the existence of costs of being close to them. However, contrarily to our expectations, herbivores did not segregate from elephants the rest of the dry season but tended to increasingly aggregate with elephants as the dry season progressed. We discuss these surprising results and the existence of a trade-off between avoidance of interspecific interference competition and other potential factors such as access to quality water, which may have relative associated costs that change with the time of the year.

Spatial Distribution of a Large Herbivore Community at Waterholes: An Assessment of Its Stability over Years in Hwange National Park, Zimbabwe.

The spatial structuring of populations or communities is an important driver of their functioning and their influence on ecosystems. Identifying the (in)stability of the spatial structure of populations is a first step towards understanding the underlying causes of these structures. Here we studied the relative importance of spatial vs. interannual variability in explaining the patterns of abundance of a large herbivore community (8 species) at waterholes in Hwange National Park (Zimbabwe). We analyzed census data collected over 13 years using multivariate methods. Our results showed that variability in the census data was mostly explained by the spatial structure of the community, as some waterholes had consistently greater herbivore abundance than others. Some temporal variability probably linked to Park-scale migration dependent on annual rainfall was noticeable, however. Once this was accounted for, little temporal variability remained to be explained, suggesting that other factors affecting herbivore abundance over time had a negligible effect at the scale of the study. The extent of spatial and temporal variability in census data was also measured for each species. This study could help in projecting the consequences of surface water management, and more generally presents a methodological framework to simultaneously address the relative importance of spatial vs. temporal effects in driving the distribution of organisms across landscapes.

Addressing ecological effects of radiation on populations and ecosystems to improve protection of the environment against radiation: Agreed statements from a Consensus Symposium.

This paper reports the output of a consensus symposium organized by the International Union of Radioecology in November 2015. The symposium gathered an academically diverse group of 30 scientists to consider the still debated ecological impact of radiation on populations and ecosystems. Stimulated by the Chernobyl and Fukushima disasters' accidental contamination of the environment, there is increasing interest in developing environmental radiation protection frameworks. Scientific research conducted in a variety of laboratory and field settings has improved our knowledge of the effects of ionizing radiation on the environment. However, the results from such studies sometimes appear contradictory and there is disagreement about the implications for risk assessment. The Symposium discussions therefore focused on issues that might lead to different interpretations of the results, such as laboratory versus field approaches, organism versus population and ecosystemic inference strategies, dose estimation approaches and their significance under chronic exposure conditions. The participating scientists, from across the spectrum of disciplines and research areas, extending also beyond the traditional radioecology community, successfully developed a constructive spirit directed at understanding discrepancies. From the discussions, the group has derived seven consensus statements related to environmental protection against radiation, which are supplemented with some recommendations. Each of these statements is contextualized and discussed in view of contributing to the orientation and integration of future research, the results of which should yield better consensus on the ecological impact of radiation and consolidate suitable approaches for efficient radiological protection of the environment.