Joint species distribution modelling with the r-package Hmsc

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

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

Land-use homogenization reduces the occurrence and diversity of frugivorous birds in a tropical biodiversity hotspot

Understanding how human-modified landscapes maintain biodiversity and provide ecosystem services is crucial for establishing conservation practices. Given that responses to land-use are species-specific, it is crucial to understand how land-use changes may shape patterns of species diversity and persistence in human-modified landscapes. Here, we used a comprehensive data set on bird distribution from the Brazilian Atlantic Forest to understand how species richness and individual occurrences of frugivorous bird species responded to land-use spatial predictors and, subsequently, assess how ecological traits and phylogeny modulated these responses. Using Bayesian hierarchical modeling, we reveal that the richness of frugivorous birds was positively associated with the amount of native forest and negatively with both agriculture and pasture amount at the landscape scale. Conversely, the effect of these predictors on species occurrence and ecological traits was highly variable and presented a weak phylogenetic signal. Furthermore, land-use homogenization (i.e., the conversion of forest to pasture or agriculture) led to pervasive consequences for forest-dependent bird species, whereas several generalist species thrived in deforested areas, replacing those sensitive to habitat disturbances.

Global Spore Sampling Project: A global, standardized dataset of airborne fungal DNA

Novel methods for sampling and characterizing biodiversity hold great promise for re-evaluating patterns of life across the planet. The sampling of airborne spores with a cyclone sampler, and the sequencing of their DNA, have been suggested as an efficient and well-calibrated tool for surveying fungal diversity across various environments. Here we present data originating from the Global Spore Sampling Project, comprising 2,768 samples collected during two years at 47 outdoor locations across the world. Each sample represents fungal DNA extracted from 24 m3 of air. We applied a conservative bioinformatics pipeline that filtered out sequences that did not show strong evidence of representing a fungal species. The pipeline yielded 27,954 species-level operational taxonomic units (OTUs). Each OTU is accompanied by a probabilistic taxonomic classification, validated through comparison with expert evaluations. To examine the potential of the data for ecological analyses, we partitioned the variation in species distributions into spatial and seasonal components, showing a strong effect of the annual mean temperature on community composition.

Predicting the responses of European grassland communities to climate and land cover change

European grasslands are among the most species-rich ecosystems on small spatial scales. However, human-induced activities like land use and climate change pose significant threats to this diversity. To explore how climate and land cover change will affect biodiversity and community composition in grassland ecosystems, we conducted joint species distribution models (SDMs) on the extensive vegetation-plot database sPlotOpen to project distributions of 1178 grassland species across Europe under current conditions and three future scenarios. We further compared model accuracy and computational efficiency between joint SDMs (JSDMs) and stacked SDMs, especially for rare species. Our results show that: (i) grassland communities in the mountain ranges are expected to suffer high rates of species loss, while those in western, northern and eastern Europe will experience substantial turnover; (ii) scaling anomalies were observed in the predicted species richness, reflecting regional differences in the dominant drivers of assembly processes; (iii) JSDMs did not outperform stacked SDMs in predictive power but demonstrated superior efficiency in model fitting and predicting; and (iv) incorporating co-occurrence datasets improved the model performance in predicting the distribution of rare species. This article is part of the theme issue 'Ecological novelty and planetary stewardship: biodiversity dynamics in a transforming biosphere'.

A successional shift enhances stability in ant symbiont communities

Throughout succession, communities undergo structural shifts, which can alter the relative abundances of species and how they interact. It is frequently asserted that these alterations beget stability, i.e. that succession selects for communities better able to resist perturbations. Yet, whether and how alterations of network structure affect stability during succession in complex communities is rarely studied in natural ecosystems. Here, we explore how network attributes influence stability of different successional stages of a natural network: symbiotic arthropod communities forming food webs inside red wood ant nests. We determined the abundance of 16 functional groups within the symbiont community across 51 host nests in the beginning and end stages of succession. Nest age was the main driver of the compositional shifts: symbiont communities in old nests contained more even species abundance distributions and a greater proportion of specialists. Based on the abundance data, we reconstructed interaction matrices and food webs of the symbiont community for each nest. We showed that the enhanced community evenness in old nests leads to an augmented food web stability in all but the largest symbiont communities. Overall, this study demonstrates that succession begets stability in a natural ecological network by making the community more even.

Transposable elements regulate thymus development and function

Transposable elements (TEs) are repetitive sequences representing ~45% of the human and mouse genomes and are highly expressed by medullary thymic epithelial cells (mTECs). In this study, we investigated the role of TEs on T-cell development in the thymus. We performed multiomic analyses of TEs in human and mouse thymic cells to elucidate their role in T-cell development. We report that TE expression in the human thymus is high and shows extensive age- and cell lineage-related variations. TE expression correlates with multiple transcription factors in all cell types of the human thymus. Two cell types express particularly broad TE repertoires: mTECs and plasmacytoid dendritic cells (pDCs). In mTECs, transcriptomic data suggest that TEs interact with transcription factors essential for mTEC development and function (e.g., PAX1 and REL), and immunopeptidomic data showed that TEs generate MHC-I-associated peptides implicated in thymocyte education. Notably, AIRE, FEZF2, and CHD4 regulate small yet non-redundant sets of TEs in murine mTECs. Human thymic pDCs homogenously express large numbers of TEs that likely form dsRNA, which can activate innate immune receptors, potentially explaining why thymic pDCs constitutively secrete IFN ɑ/β. This study highlights the diversity of interactions between TEs and the adaptive immune system. TEs are genetic parasites, and the two thymic cell types most affected by TEs (mTEcs and pDCs) are essential to establishing central T-cell tolerance. Therefore, we propose that orchestrating TE expression in thymic cells is critical to prevent autoimmunity in vertebrates.

Unravelling the impact of climate change on honey bees: An ensemble modelling approach to predict shifts in habitat suitability in Queensland, Australia

Honey bees play a vital role in providing essential ecosystem services and contributing to global agriculture. However, the potential effect of climate change on honey bee distribution is still not well understood. This study aims to identify the most influential bioclimatic and environmental variables, assess their impact on honey bee distribution, and predict future distribution. An ensemble modelling approach using the biomod2 package in R was employed to develop three models: a climate-only model, an environment-only model, and a combined climate and environment model. By utilising bioclimatic data (radiation of the wettest and driest quarters and temperature seasonality) from 1990 to 2009, combined with observed honey bee presence and pseudo absence data, this model predicted suitable locations for honey bee apiaries for two future time spans: 2020-2039 and 2060-2079. The climate-only model exhibited a true skill statistic (TSS) value of 0.85, underscoring the pivotal role of radiation and temperature seasonality in shaping honey bee distribution. The environment-only model, incorporating proximity to floral resources, foliage projective cover, and elevation, demonstrated strong predictive performance, with a TSS of 0.88, emphasising the significance of environmental variables in determining habitat suitability for honey bees. The combined model had a higher TSS of 0.96, indicating that the combination of climate and environmental variables enhances the model's performance. By the 2020-2039 period, approximately 88% of highly suitable habitats for honey bees are projected to transition from their current state to become moderate (14.84%) to marginally suitable (13.46%) areas. Predictions for the 2060-2079 period reveal a concerning trend: 100% of highly suitable land transitions into moderately (0.54%), marginally (17.56%), or not suitable areas (81.9%) for honey bees. These results emphasise the critical need for targeted conservation efforts and the implementation of policies aimed at safeguarding honey bees and the vital apiary industry.

Functional redundancy of weed seed predation is reduced by intensified agriculture

Intensified agriculture, a driver of biodiversity loss, can diminish ecosystem functions and their stability. Biodiversity can increase functional redundancy and is expected to stabilize ecosystem functions. Few studies, however, have explored how agricultural intensity affects functional redundancy and its link with ecosystem function stability. Here, within a continental-wide study, we assess how functional redundancy of seed predation is affected by agricultural intensity and landscape simplification. By combining carabid abundances with molecular gut content data, functional redundancy of seed predation was quantified for 65 weed genera across 60 fields in four European countries. Across weed genera, functional redundancy was reduced with high field management intensity and simplified crop rotations. Moreover, functional redundancy increased the spatial stability of weed seed predation at the field scale. We found that ecosystem functions are vulnerable to disturbances in intensively managed agroecosystems, providing empirical evidence of the importance of biodiversity for stable ecosystem functions across space.

Examining the co-occurrences of human threats within terrestrial protected areas

Human threats to biodiversity are prevalent within protected areas (PAs), undermining their effectiveness in halting biodiversity loss. Certain threats tend to co-occur, resulting in amplified cumulative impact through synergistic effects. However, it remains unclear which threats are related the most. We analyzed a dataset of 71 human threats in 18 013 terrestrial PAs of the European Union's Natura 2000 network, using a Joint Species Distribution Modelling approach, to assess the threats' co-occurrence patterns and potential drivers. Overall, threats were more frequently correlated positively than negatively. Threats related to agriculture and urbanization were correlated strongly with most other threats. Approximately 70% of the variance in our model was explained by country-specific factors, indicating the importance of local drivers. Minimizing the negative impact of key threats can likely reduce the impact of related threats. However, more research is needed to understand better the relationships among threats and, importantly, their combined impact on biodiversity.

Rewilding soil and litter invertebrates and fungi increases decomposition rates and alters detritivore communities

Habitat degradation and associated reductions in ecosystem functions can be reversed by reintroducing or 'rewilding' keystone species. Rewilding projects have historically targeted restoration of processes such as grazing regimes or top-down predation effects. Few projects focus on restoring decomposition efficiency, despite the pivotal role decomposition plays in global carbon sequestration and nutrient cycling. Here, we tested whether rewilding entire communities of detritivorous invertebrates and fungi can improve litter decomposition efficiency and restore detritivore communities during ecological restoration. Rewilding was conducted by transplanting leaf litter and soil, including associated invertebrate and fungal communities from species-rich remnant sites into species-poor, and geographically isolated, revegetated farmland sites in a temperate woodland region of southeastern Australia. We compared communities in sites under the following treatments: remnant (conservation area and source of litter transplant), rewilded revegetation (revegetated farmland site with litter transplant) and control revegetation (revegetated site, no transplant). In one 'before' and three 'after' sampling periods, we measured litter decomposition and the abundance and diversity of detritivorous invertebrates and fungi. We quantified the effect of detritivores on the rate of litter decomposition using piecewise Structural Equation Modelling. Decomposition was significantly faster in rewilding sites than in both control and remnant areas and was largely driven by a greater abundance of invertebrate detritivores. Similarly, the abundance of invertebrate detritivores in rewilding revegetation sites exceeded the level of remnant communities, whereas there was little difference between control and remnant sites. In contrast, rewilding did not increase saprotrophic fungi relative abundance/diversity and there was no strong relationship between decomposition and fungal diversity. Our findings suggest the relatively simple act of transplanting leaf litter and soil can increase functional efficiency during restoration and alter community composition. Our methods may prove important across a range of contexts where other restoration methods have failed to restore ecosystem processes to pre-degradation levels.

How non-target chironomid communities respond to mosquito control: Integrating DNA metabarcoding and joint species distribution modelling

The conservation and management of riparian ecosystems rely on understanding the ecological consequences of anthropogenic stressors that impact natural communities. In this context, studies investigating the effects of anthropogenic stressors require reliable methods capable of mapping the relationships between taxa occurrence or abundance and environmental predictors within a spatio-temporal framework. Here, we present an integrative approach using DNA metabarcoding and Hierarchical Modelling of Species Communities (HMSC) to unravel the intricate dynamics and resilience of chironomid communities exposed to Bacillus thuringiensis var. israelensis (Bti). Chironomid emergence was sampled from a total of 12 floodplain pond mesocosms, half of which received Bti treatment, during a 16-week period spanning spring and summer of 2020. Subsequently, we determined the community compositions of chironomids and examined their genus-specific responses to the Bti treatment, considering their phylogenetic affiliations and ecological traits of the larvae. Additionally, we investigated the impact of the Bti treatment on the body size distribution of emerging chironomids. Our study revealed consistent responses to Bti among different chironomid genera, indicating that neither phylogenetic affiliations nor larval feeding strategies significantly contributed to the observed patterns. Both taxonomic and genetic diversity were positively correlated with the number of emerged individuals. Furthermore, our findings demonstrated Bti-related effects on chironomid body size distribution, which could have relevant implications for size-selective terrestrial predators. Hence, our study highlights the value of employing a combination of DNA metabarcoding and HMSC to unravel the complex dynamics of Bti-related non-target effects on chironomid communities. The insights gained from this integrated framework contribute to our understanding of the ecological consequences of anthropogenic stressors and provide a foundation for informed decision-making regarding the conservation and management of riparian ecosystems.

Spatially explicit estimation of freshwater fish stock biomass with limited data: A case study of an endangered endemic fish on the Tibetan Plateau, China

Accurate evaluation of fish stock biomass is essential for effective conservation management and targeted species enhancement efforts. However, this remains challenging owing to limited data availability. Therefore, we present an integrated modeling framework combining catch per unit effort with ensemble species distribution modeling called CPUESDM, which explicitly assesses and validates the spatial distribution of stock biomass for freshwater fish species with limited data, applied to Herzensteinia microcephalus. The core algorithm incorporates the Leslie regression model, ensemble species distribution modeling, and exploratory spatial interpolation techniques. We found that H. microcephalus biomass in the Yangtze River source area yielded an initial estimate of 113.52 tons. Our validation results demonstrate high accuracy with a Cohen's kappa coefficient of 0.78 and root mean square error of 0.05. Furthermore, our spatially-explicit, global, absolute biomass density map effectively identified areas with high and low concentrations of biomass distribution centers. Additionally, this study offers access to the source code, example raw data, and a step-by-step instruction manual for other researchers using field data to explore the application of this model. Our findings can help inform for future conservation efforts around fish stock biomass estimation, especially for endangered species.

Whether including exotic species alters conservation prioritization: a case study in the Min River in southeastern China

Conservation practices from the perspective of functional diversity (FD) and conservation prioritization need to account for the impacts of exotic species in freshwater ecosystems. This work first simulated the influence of exotic species on the values of FD in a schemed mechanistic model, and then a practical case study of conservation prioritization was performed in the Min River, the largest river in southeastern China, to discuss whether including exotic species alters prioritization. The mechanistic model revealed that exotic species significantly altered the expected FD if the number of exotic species occupied 2% of the community. Joint species distribution modelling indicated that the highest FD occurred in the west, northwest and north upstreams of the Min River. Values of FD in 64.69% of the basin decreased after the exotic species were removed from calculation. Conservation prioritization with the Zonation software proved that if first the habitats of exotic species were removed during prioritization, 62.75% of the highest prioritized areas were shifted, average species representation of the endemic species was improved and mean conservation efficiency was increased by 7.53%. Existence of exotic species will significantly alter the metrics of biodiversity and the solution for conservation prioritization, and negatively weighting exotic species in the scope of conservation prioritization is suggested to better protect endemic species. This work advocates a thorough estimate of the impacts of exotic species on FD and conservation prioritization, providing complementary evidence for conservation biology and valuable implications for local freshwater fish conservation.

Functional alterations and predictive capacity of gut microbiome in type 2 diabetes

The gut microbiome plays a significant role in the development of Type 2 Diabetes Mellitus (T2DM), but the functional mechanisms behind this association merit deeper investigation. Here, we used the nanopore sequencing technology for metagenomic analyses to compare the gut microbiome of individuals with T2DM from the United Arab Emirates (n = 40) with that of control (n = 44). DMM enterotyping of the cohort resulted concordantly with previous results, in three dominant groups Bacteroides (K1), Firmicutes (K2), and Prevotella (K3) lineages. The diversity analysis revealed a high level of diversity in the Firmicutes group (K2) both in terms of species richness and evenness (Wilcoxon rank-sum test, p value < 0.05 vs. K1 and K3 groups), consistent with the Ruminococcus enterotype described in Western populations. Additionally, functional enrichment analyses of KEGG modules showed significant differences in abundance between individuals with T2DM and controls (FDR < 0.05). These differences include modules associated with the degradation of amino acids, such as arginine, the degradation of urea as well as those associated with homoacetogenesis. Prediction analysis with the Predomics approach suggested potential biomarkers for T2DM, including a balance between a depletion of Enterococcus faecium and Blautia lineages with an enrichment of Absiella spp or Eubacterium limosum in T2DM individuals, highlighting the potential of metagenomic analysis in predicting predisposition to diabetic cardiomyopathy in T2DM patients.

Seasonal distribution patterns of Scomberomorus commerson in the Taiwan Strait in relation to oceanographic conditions: An ensemble modeling approach

The decline in the stock of the narrow-barred Spanish mackerel in the Taiwan Strait has sparked interest in conservation efforts. To optimize conservation and restoration efforts, it is crucial to understand their habitat preference in response to changing environments. In this study, ensemble modeling was used to investigate the seasonal distribution patterns of Spanish mackerel. Winter was identified as the most productive season, followed by fall; productivity was the lowest in summer. Five single-algorithm models were developed, and on the basis of their performance, four were selected for inclusion in an ensemble species distribution model. The spatial distribution of Spanish mackerel was primarily along the latitudinal range 23°-25°N in spring and summer. However, in fall and winter, the geographical range increased toward the southern region. The findings of this study will contribute to the understanding of this specific species and the approach used in this study may be applicable to other fisheries stocks also.

A protocol for evaluating the entomological impact of larval source reduction on mosquito vectors at hotel compounds in Zanzibar

There is an increasing awareness of the association between tourism activity and risks of emerging mosquito-borne diseases (MBDs) worldwide. In previous studies we showed that hotels in Zanzibar may play an important role in maintaining residual foci of mosquito vectors populations of public health concern. These findings indicated larval sources removal (LSR) interventions may have a significant negative impact on vector communities. However, a thorough analysis of the response vector species to potential LSM strategies must be evaluated prior to implementation of a large-scale area-wide control campaign. Here we propose a protocol for evaluation of the impact of LSR against mosquito vectors at hotel settings in Zanzibar. This protocol is set to determine the efficacy of LSR in a randomized control partial cross-over experimental design with four hotel compounds representing the unit of randomization for allocation of interventions. However, the protocol can be applied to evaluate the impact of LRS in more than four sites. Proposed interventions are active removal of disposed containers, and installation of water dispenser to replace single use discarded plastic water bottles, which were identified as the most important source of mosquitoes studied hotels. The ideal time for allocating intervention to the intervention arms the dry season, when the mosquito abundance is predictably lower. The possible impact of interventions on mosquito occurrence and abundance risks is then evaluated throughout subsequent rainy and dry seasons. If an appreciable reduction in mosquito abundance and occurrence risks is observed during the trial period, intervention could be extended to the control arm to determine whether any potential reduction of mosquito density is reproducible. A rigorous evaluation of the proposed LRS interventions will inspire large scale trials and provide support for evidence-based mosquito management at hotel facilities in Zanzibar and similar settings.

Interspecific Host Variation and Biotic Interactions Drive Pathogen Community Assembly in Chinese Bumblebees

Bumblebees have been considered one of the most important pollinators on the planet. However, recent reports of bumblebee decline have raised concern about a significant threat to ecosystem stability. Infectious diseases caused by multiple pathogen infections have been increasingly recognized as an important mechanism behind this decline worldwide. Understanding the determining factors that influence the assembly and composition of pathogen communities among bumblebees can provide important implications for predicting infectious disease dynamics and making effective conservation policies. Here, we study the relative importance of biotic interactions versus interspecific host resistance in shaping the pathogen community composition of bumblebees in China. We first conducted a comprehensive survey of 13 pathogens from 22 bumblebee species across China. We then applied joint species distribution modeling to assess the determinants of pathogen community composition and examine the presence and strength of pathogen-pathogen associations. We found that host species explained most of the variations in pathogen occurrences and composition, suggesting that host specificity was the most important variable in predicting pathogen occurrences and community composition in bumblebees. Moreover, we detected both positive and negative associations among pathogens, indicating the role of competition and facilitation among pathogens in determining pathogen community assembly. Our research demonstrates the power of a pluralistic framework integrating field survey of bumblebee pathogens with community ecology frameworks to understand the underlying mechanisms of pathogen community assembly.

Investigating Ecological Interactions in the Tumor Microenvironment using Joint Species Distribution Models for Point Patterns

The tumor microenvironment (TME) is a complex and dynamic ecosystem that involves interactions between different cell types, such as cancer cells, immune cells, and stromal cells. These interactions can promote or inhibit tumor growth and affect response to therapy. Multitype Gibbs point process (MGPP) models are statistical models used to study the spatial distribution and interaction of different types of objects, such as the distribution of cell types in a tissue sample. Such models are potentially useful for investigating the spatial relationships between different cell types in the tumor microenvironment, but so far studies of the TME using cell-resolution imaging have been largely limited to spatial descriptive statistics. However, MGPP models have many advantages over descriptive statistics, such as uncertainty quantification, incorporation of multiple covariates and the ability to make predictions. In this paper, we describe and apply a previously developed MGPP method, the saturated pairwise interaction Gibbs point process model , to a publicly available multiplexed imaging dataset obtained from colorectal cancer patients. Importantly, we show how these methods can be used as joint species distribution models (JSDMs) to precisely frame and answer many relevant questions related to the ecology of the tumor microenvironment.

Quantitative Amplicon Sequencing Is Necessary to Identify Differential Taxa and Correlated Taxa Where Population Sizes Differ

High-throughput, multiplexed-amplicon sequencing has become a core tool for understanding environmental microbiomes. As researchers have widely adopted sequencing, many open-source analysis pipelines have been developed to compare microbiomes using compositional analysis frameworks. However, there is increasing evidence that compositional analyses do not provide the information necessary to accurately interpret many community assembly processes. This is especially true when there are large gradients that drive distinct community assembly processes. Recently, sequencing has been combined with Q-PCR (among other sources of total quantitation) to generate "Quantitative Sequencing" (QSeq) data. QSeq more accurately estimates the true abundance of taxa, is a more reliable basis for inferring correlation, and, ultimately, can be more reliably related to environmental data to infer community assembly processes. In this paper, we use a combination of published data sets, synthesis, and empirical modeling to offer guidance for which contexts QSeq is advantageous. As little as 5% variation in total abundance among experimental groups resulted in more accurate inference by QSeq than compositional methods. Compositional methods for differential abundance and correlation unreliably detected patterns in abundance and covariance when there was greater than 20% variation in total abundance among experimental groups. Whether QSeq performs better for beta diversity analysis depends on the question being asked, and the analytic strategy (e.g., what distance metric is being used); for many questions and methods, QSeq and compositional analysis are equivalent for beta diversity analysis. QSeq is especially useful for taxon-specific analysis; QSeq transformation and analysis should be the default for answering taxon-specific questions of amplicon sequence data. Publicly available bioinformatics pipelines should incorporate support for QSeq transformation and analysis.

Mammals show distinct functional gut microbiome dynamics to identical series of environmental stressors

IMPORTANCE

In our manuscript, we report the first interspecific comparative study about the plasticity of the gut microbiota. We conducted a captivity experiment that exposed wild-captured mammals to a series of environmental challenges over 45 days. We characterized their gut microbial communities using genome-resolved metagenomics and modeled how the taxonomic, phylogenetic, and functional microbial dynamics varied across a series of disturbances in both species. Our results indicate that the intrinsic properties (e.g., diversity and functional redundancy) of microbial communities coupled with physiological attributes (e.g., thermal plasticity) of hosts shape the taxonomic, phylogenetic, and functional response of gut microbiomes to environmental stressors, which might influence their contribution to the acclimation and adaptation capacity of animal hosts.

Accounting for trait variability and coordination in predictions of drought-induced range shifts in woody plants

Functional traits offer a promising avenue to improve predictions of species range shifts under climate change, which will entail warmer and often drier conditions. Although the conceptual foundation linking traits with plant performance and range shifts appears solid, the predictive ability of individual traits remains generally low. In this review, we address this apparent paradox, emphasizing examples of woody plants and traits associated with drought responses at the species' rear edge. Low predictive ability reflects the fact not only that range dynamics tend to be complex and multifactorial, as well as uncertainty in the identification of relevant traits and limited data availability, but also that trait effects are scale- and context-dependent. The latter results from the complex interactions among traits (e.g. compensatory effects) and between them and the environment (e.g. exposure), which ultimately determine persistence and colonization capacity. To confront this complexity, a more balanced coverage of the main functional dimensions involved (stress tolerance, resource use, regeneration and dispersal) is needed, and modelling approaches must be developed that explicitly account for: trait coordination in a hierarchical context; trait variability in space and time and its relationship with exposure; and the effect of biotic interactions in an ecological community context.

Alternative measures of trait-niche relationships: A test on dispersal traits in saproxylic beetles

Functional trait approaches are common in ecology, but a lack of clear hypotheses on how traits relate to environmental gradients (i.e., trait-niche relationships) often makes uncovering mechanisms difficult. Furthermore, measures of community functional structure differ in their implications, yet inferences are seldom compared among metrics. Community-weighted mean trait values (CWMs), a common measure, are largely driven by the most common species and thus do not reflect community-wide trait-niche relationships per se. Alternatively, trait-niche relationships can be estimated across a larger group of species using hierarchical joint species distribution models (JSDMs), quantified by a parameter Γ. We investigated how inferences about trait-niche relationships are affected by the choice of metric. Using deadwood-dependent (saproxylic) beetles in fragmented Finnish forests, we followed a protocol for investigating trait-niche relationships by (1) identifying environmental filters (climate, forest age, and deadwood volume), (2) relating these to an ecological function (dispersal ability), and (3) identifying traits related to this function (wing morphology). We tested 18 hypothesized dispersal relationships using both CWM and Γ estimates across these environmental gradients. CWMs were more likely than Γ to show support for trait-niche relationships. Up to 13% of species' realized niches were explained by dispersal traits, but the directions of effects were consistent with fewer than 11%-39% of our 18 trait-niche hypotheses (depending on the metric used). This highlights the difficulty in connecting morphological traits and ecological functions in insects, despite the clear conceptual link between landscape connectivity and flight-related traits. Caution is thus warranted in hypothesis development, particularly where apparent trait-function links are less clear. Inferences differ when CWMs versus Γ estimates are used, necessitating the choice of a metric that reflects study questions. CWMs help explain the effects of environmental gradients on community trait composition, whereas the effects of traits on species' niches are better estimated using hierarchical JSDMs.

The first arriving virus shapes within-host viral diversity during natural epidemics

Viral diversity has been discovered across scales from host individuals to populations. However, the drivers of viral community assembly are still largely unknown. Within-host viral communities are formed through co-infections, where the interval between the arrival times of viruses may vary. Priority effects describe the timing and order in which species arrive in an environment, and how early colonizers impact subsequent community assembly. To study the effect of the first-arriving virus on subsequent infection patterns of five focal viruses, we set up a field experiment using naïve Plantago lanceolata plants as sentinels during a seasonal virus epidemic. Using joint species distribution modelling, we find both positive and negative effects of early season viral infection on late season viral colonization patterns. The direction of the effect depends on both the host genotype and which virus colonized the host early in the season. It is well established that co-occurring viruses may change the virulence and transmission of viral infections. However, our results show that priority effects may also play an important, previously unquantified role in viral community assembly. The assessment of these temporal dynamics within a community ecological framework will improve our ability to understand and predict viral diversity in natural systems.

Mixed effects of a national protected area network on terrestrial and freshwater biodiversity

Protected areas are considered fundamental to counter biodiversity loss. However, evidence for their effectiveness in averting local extinctions remains scarce and taxonomically biased. We employ a robust counterfactual multi-taxon approach to compare occupancy patterns of 638 species, including birds (150), mammals (23), plants (39) and phytoplankton (426) between protected and unprotected sites across four decades in Finland. We find mixed impacts of protected areas, with only a small proportion of species explicitly benefiting from protection-mainly through slower rates of decline inside protected areas. The benefits of protection are enhanced for larger protected areas and are traceable to when the sites were protected, but are mostly unrelated to species conservation status or traits (size, climatic niche and threat status). Our results suggest that the current protected area network can partly contribute to slow down declines in occupancy rates, but alone will not suffice to halt the biodiversity crisis. Efforts aimed at improving coverage, connectivity and management will be key to enhance the effectiveness of protected areas towards bending the curve of biodiversity loss.

Statistical methods to identify mechanisms in studies of eco-evolutionary dynamics

While the reciprocal effects of ecological and evolutionary dynamics are increasingly recognized as an important driver for biodiversity, detection of such eco-evolutionary feedbacks, their underlying mechanisms, and their consequences remains challenging. Eco-evolutionary dynamics occur at different spatial and temporal scales and can leave signatures at different levels of organization (e.g., gene, protein, trait, community) that are often difficult to detect. Recent advances in statistical methods combined with alternative hypothesis testing provides a promising approach to identify potential eco-evolutionary drivers for observed data even in non-model systems that are not amenable to experimental manipulation. We discuss recent advances in eco-evolutionary modeling and statistical methods and discuss challenges for fitting mechanistic models to eco-evolutionary data.

Demography-environment relationships improve mechanistic understanding of range dynamics under climate change

Species respond to climate change with range and abundance dynamics. To better explain and predict them, we need a mechanistic understanding of how the underlying demographic processes are shaped by climatic conditions. Here, we aim to infer demography-climate relationships from distribution and abundance data. For this, we developed spatially explicit, process-based models for eight Swiss breeding bird populations. These jointly consider dispersal, population dynamics and the climate-dependence of three demographic processes-juvenile survival, adult survival and fecundity. The models were calibrated to 267 nationwide abundance time series in a Bayesian framework. The fitted models showed moderate to excellent goodness-of-fit and discriminatory power. The most influential climatic predictors for population performance were the mean breeding-season temperature and the total winter precipitation. Contemporary climate change benefitted the population trends of typical mountain birds leading to lower population losses or even slight increases, whereas lowland birds were adversely affected. Our results emphasize that generic process-based models embedded in a robust statistical framework can improve our predictions of range dynamics and may allow disentangling of the underlying processes. For future research, we advocate a stronger integration of experimental and empirical studies in order to gain more precise insights into the mechanisms by which climate affects populations. This article is part of the theme issue 'Detecting and attributing the causes of biodiversity change: needs, gaps and solutions'.

Hurdles and opportunities in implementing marine biosecurity systems in data-poor regions

Managing marine nonindigenous species (mNIS) is challenging, because marine environments are highly connected, allowing the dispersal of species across large spatial scales, including geopolitical borders. Cross-border inconsistencies in biosecurity management can promote the spread of mNIS across geopolitical borders, and incursions often go unnoticed or unreported. Collaborative surveillance programs can enhance the early detection of mNIS, when response may still be possible, and can foster capacity building around a common threat. Regional or international databases curated for mNIS can inform local monitoring programs and can foster real-time information exchange on mNIS of concern. When combined, local species reference libraries, publicly available mNIS databases, and predictive modeling can facilitate the development of biosecurity programs in regions lacking baseline data. Biosecurity programs should be practical, feasible, cost-effective, mainly focused on prevention and early detection, and be built on the collaboration and coordination of government, nongovernment organizations, stakeholders, and local citizens for a rapid response.

Local eukaryotic and bacterial stream community assembly is shaped by regional land use effects

With anticipated expansion of agricultural areas for food production and increasing intensity of pressures stemming from land-use, it is critical to better understand how species respond to land-use change. This is particularly true for microbial communities which provide key ecosystem functions and display fastest responses to environmental change. However, regional land-use effects on local environmental conditions are often neglected, and, hence, underestimated when investigating community responses. Here we show that the effects stemming from agricultural and forested land use are strongest reflected in water conductivity, pH and phosphorus concentration, shaping microbial communities and their assembly processes. Using a joint species distribution modelling framework with community data based on metabarcoding, we quantify the contribution of land-use types in determining local environmental variables and uncover the impact of both, land-use, and local environment, on microbial stream communities. We found that community assembly is closely linked to land-use type but that the local environment strongly mediates the effects of land-use, resulting in systematic variation of taxon responses to environmental conditions, depending on their domain (bacteria vs. eukaryote) and trophic mode (autotrophy vs. heterotrophy). Given that regional land-use type strongly shapes local environments, it is paramount to consider its key role in shaping local stream communities.

Influence of ecological traits on spatio-temporal dynamics of an elasmobranch community in a heavily exploited basin

Elasmobranchs, which include sharks and batoids, play critical roles in maintaining the integrity and stability of marine food webs. However, these cartilaginous fish are among the most threatened vertebrate lineages due to their widespread depletion. Consequently, understanding dynamics and predicting changes of elasmobranch communities are major research topics in conservation ecology. Here, we leverage long-term catch data from a standardized bottom trawl survey conducted from 1996 to 2019, to evaluate the spatio-temporal dynamics of the elasmobranch community in the heavily exploited Adriatic Sea, where these fish have historically been depleted. We use joint species distribution modeling to quantify the responses of the species to environmental variation while also including important traits such as species age at first maturity, reproductive mode, trophic level, and phylogenetic information. We present spatio-temporal changes in the species community and associated modification of the trait composition, highlighting strong spatial and depth-mediated patterning. We observed an overall increase in the abundance of the dominant elasmobranch species, except for spurdog, which has shown a continued decline. However, our results showed that the present community displays lower age at first maturity and a smaller fraction of viviparous species compared to the earlier observed community due to changes in species' relative abundance. The selected traits contributed considerably to explaining community patterns, suggesting that the integration of trait-based approaches in elasmobranch community analyses can aid efforts to conserve this important lineage of fish.

Tick-borne pathogens and body condition of cattle in smallholder rural livestock production systems in East and West Africa

BACKGROUND

The majority of the African population lives in rural areas where they heavily depend on crop and livestock production for their livelihoods. Given their socio-economic importance, we initiated a standardized multi-country (Benin, Burkina Faso, Ghana, Nigeria, Ethiopia Tanzania and Uganda) surveillance study to assess the current status of important tick-borne haemoparasites (TBHPs) of cattle.

METHODS

We assessed pathogen prevalences (Anaplasma marginale, Anaplasma centrale, Babesia bigemina, Babesia bovis, Ehrlichia ruminantium, and Theileria parva) in the blood of 6447 animals spread over fourteen districts (two districts per country). In addition, we screened for intrinsic (sex, weight, body condition) and extrinsic (husbandry, tick exposure) risk factors as predictors of infections with TBHPs.

RESULTS

There was a large macro-geographic variation observed in A. marginale, B. bigemina, B. bovis and E. ruminantium prevalences. Most correlated with the co-occurrence of their specific sets of vector-competent ticks. Highest numbers of infected cattle were found in Ghana and Benin, and lowest in Burkina Faso. While T. parva was seldomly found (Uganda only: 3.0%), A. marginale was found in each country with a prevalence of at least 40%. Babesia bovis infected individuals had lower body condition scores. Age (as estimated via body weight) was higher in A. marginale infected cattle, but was negatively correlated with B. bigemina and E. ruminantium prevalences. Ehrlichia ruminantium infection was more often found in males, and A. marginale more often in transhumance farming. High levels of co-infection, especially the combination A. marginale × B. bigemina, were observed in all countries, except for Uganda and Burkina Faso. Babesia bigemina was more or less often observed than expected by chance, when cattle were also co-infected with E. ruminantium or A. marginale, respectively.

CONCLUSIONS

Tick-borne pathogens of cattle are ubiquitous in African's smallholder cattle production systems. Our standardized study will help a wide range of stakeholders to provide recommendations for TBHP surveillance and prevention in cattle, especially for B. bovis which heavily impacts production and continues its spread over the African continent via the invasive Rhipicephalus microplus tick.

Biological Crust Diversity Related to Elevation and Soil Properties at Local Scale in a Montane Scrub of Ecuador

The montane shrublands of southern Ecuador represent one of the least studied ecosystems, which in the last decade have been seriously threatened by increasing wildfires, deforestation, overgrazing, and conversion to forest plantations. Our main objective was to determine, at the local scale, the diversity of species composing the biological soil crust (BSC) at three elevations (2100, 2300, and 2500 m.a.s.l.) and their possible relationships with soil physical and chemical properties in montane shrublands. For this purpose, three monitoring plots of 100 m² were established at each elevation, and within each plot, 20 subplots were established (180 subplots sampled in total). In addition, composite soil samples were collected at a depth of 0 to 10 cm, and some physical and biochemical parameters (e.g., bulk density, texture, pH, organic matter, soil organic carbon, total nitrogen, available phosphorus, and potassium) of the soil were analyzed. The results show 35 species (23 lichens, 10 bryophytes and 2 cyanobacteria) at three elevations with a bell-shaped or hump-shaped distribution pattern. This allowed us to point out that the species richness was higher at the intermediate elevations and that the composition showed significant differences in the three elevations related to soil factors. Elevation and soil drivers may help to better chose the more suitable biological soil crust (lichen-dominated and bryophyte-dominated BSC) for the management and conservation of the montane scrub of Ecuador, which is strongly threatened by human activities.

Habitat partitioning, co-occurrence patterns, and mixed-species group formation in sympatric delphinids

Numerous species have been reported to form mixed-species groups, however, little is known about the interplay between niche partitioning and mixed-species group formation. Furthermore, it is often unclear whether species come together by chance due to overlapping habitat preferences, by shared attraction to resources, or by attraction between them. We assessed habitat partitioning, co-occurrence patterns, and mixed-species group formation of sympatric Australian humpback (Sousa sahulensis) and Indo-Pacific bottlenose dolphins (Tursiops aduncus) around the North West Cape, Western Australia, with a joint species distribution model and temporal analyses of sighting data. Australian humpback dolphins preferred shallower and more nearshore waters than Indo-Pacific bottlenose dolphins, yet these species co-occurred more often than expected by chance given shared responses to environmental variables. Indo-Pacific bottlenose dolphins were sighted more often than Australian humpback dolphins during the afternoon, however, we did not find any temporal patterns in the occurrence of mixed-species groups. We propose that the positive association in the species' occurrence indicates the active formation of mixed-species groups. By evaluating habitat partitioning and co-occurrence patterns, this study provides direction for future work which should proceed to investigate the benefits that these species may gain from grouping with each other.

Environmental filtering drives the assembly of mammal communities in a heterogeneous Mediterranean region

Agricultural expansion and intensification are major drivers of global change. Quantifying the importance of different processes governing the assembly of local communities in agroecosystems is essential to guide the conservation effort allocated to enhancing habitat connectivity, improving habitat quality or managing species interactions. We used multiple detection methods to record the occurrence of medium-sized and large-sized mammals in three managed landscapes of a heterogeneous Mediterranean region. Then we used a joint species distribution model to evaluate the relative influence of dispersal limitation, environmental filtering, and interspecific interactions on the local assembly of mammal communities in 4-km² plots. The partitioning of the explained variation in species occurrence was attributed on average 99% to environmental filters and 1% to dispersal filters. No role was attributed to biotic filters, in agreement with the scarce support for strong competition or other negative interactions found after a literature review. Four principal environmental factors explained on average 63% of variance in species occurrence and operated mainly at the landscape scale. The amount of shrub cover in the neighboring landscape was the most influential factor favoring mammal occurrence and accounted for nearly one-third of the total variance. The proportion of intensively managed croplands and proxies of human activity within landscape samples limited mammal presence. At the microhabitat scale (~80 m² plots) the mean percentage area deprived of woody vegetation also had a negative effect. Functional traits such as body mass or social behavior accounted for a substantial fraction of the variation attributed to environmental factors. We concluded that multiscale environmental filtering governed local community assembly, whereas the role of dispersal limitation and interspecific interactions was negligible. Our results suggest that further removal of shrubland, the expansion of intensive agriculture, and the increase of human activity are expected to result in species losses. The fact that community integrity responds to a single type of ecological process simplifies practical recommendations. Management strategies should focus on the conservation and restoration of shrubland, adopting alternatives to intensive schemes of agricultural production, and minimizing recreational and other human activities in remnant natural habitats within agroecosystems or mosaic landscapes.

Butterfly foraging is remarkably synchronous in an experimental tropical macrocosm

Diel patterns in foraging activity are dictated by a combination of abiotic, biotic and endogenous limits. Understanding these limits is important for insects because ectotherm taxa will respond more pronouncedly to ongoing climatic change, potentially affecting crucial ecosystem services. We leverage an experimental macrocosm, the Montreal Insectarium Grand Vivarium, to test the importance of endogenous mechanisms in determining temporal patterns in foraging activity of butterflies. Specifically, we assessed the degree of temporal niche partitioning among 24 butterfly species originating from the Earth's tropics within controlled environmental conditions. We found strong niche overlap, with the frequency of foraging events peaking around solar noon for 96% of the species assessed. Our models suggest that this result was not due to the extent of cloud cover, which affects radiational heating and thus limits body temperature in butterflies. Together, these findings suggest that an endogenous mechanism evolved to regulate the timing of butterfly foraging activity within suitable environmental conditions. Understanding similar mechanisms will be crucial to forecast the effects of climate change on insects, and thus on the many ecosystem services they provide.

Whole-of-community invertebrate rewilding: Leaf litter transplants rapidly increase beetle diversity during restoration

Restoration of degraded areas is now a central tool in humanity's response to continued species-loss. However, restoration projects often report exceedingly slow or failed recolonization of fauna, especially dispersal-constrained groups such as invertebrates. Active interventions via reintroducing or "rewilding" invertebrates may assist recolonization and speed up restoration of communities toward a desired target. However, invertebrate rewilding is rarely implemented during ecological restoration. Here, we studied the efficacy of invertebrate rewilding as a means of reintroducing dispersal-constrained species and improving diversity and compositional similarities to remnant communities during restoration. Rewilding was conducted by transplanting leaf litter and soil, including associated communities of invertebrates from species rich remnant sites into species poor, and geographically isolated, revegetated farmland sites. We sampled pre- and post-rewilding invertebrate communities in remnant, rewilded revegetation, and control revegetation sites. We analyzed morphospecies richness, abundance, community composition, and modeled morphospecies traits (dispersal method/trophic guild) using a Hierarchical Modelling of Species Communities approach to determine which biological properties facilitated establishment. Beetle (Coleoptera) morphospecies richness increased rapidly in rewilded sites and was indistinguishable from remnant communities as early as 7 months post-rewilding. Beetle community similarity in the rewilding sites significantly deviated from the control sites 27 months post-rewilding, however remnant communities remained distinct over the study timeframe. Establishment success varied as other taxa did not respond as consistently as beetles within the study timeframe. Furthermore, there were no discernible shifts in dispersal traits in rewilded sites. However, predatory morphospecies were more likely to establish post-rewilding than other trophic groups. Our results demonstrate that the relatively simple act of transplanting leaf litter can result in comparatively large increases in morphospecies richness during restoration in a short timeframe. We advocate methodologies such as ours should be adopted more frequently to address failed community restoration as they are cost-effective and can be easily applied by practitioners in various restoration settings. However, further efficacy tests (e.g., varying the number of rewilding events) and longer study timeframes are needed to ensure effectiveness for a broader range of invertebrate taxa and ecosystems.

Macrosystem community change in lake phytoplankton and its implications for diversity and function

Aim

We use lake phytoplankton community data to quantify the spatio-temporal and scale-dependent impacts of eutrophication, land-use and climate change on species niches and community assembly processes while accounting for species traits and phylogenetic constraints.

Location

Finland.

Time period

1977-2017.

Major taxa

Phytoplankton.

Methods

We use hierarchical modelling of species communities (HMSC) to model metacommunity trajectories at 853 lakes over four decades of environmental change, including a hierarchical spatial structure to account for scale-dependent processes. Using a "region of common profile" approach, we evaluate compositional changes of species communities and trait profiles and investigate their temporal development.

Results

We demonstrate the emergence of novel and widespread community composition clusters in previously more compositionally homogeneous communities, with cluster-specific community trait profiles, indicating functional differences. A strong phylogenetic signal of species responses to the environment implies similar responses among closely related taxa. Community cluster-specific species prevalence indicates lower taxonomic dispersion within the current dominant clusters compared with the historically dominant cluster and an overall higher prevalence of smaller species sizes within communities. Our findings denote profound spatio-temporal structuring of species co-occurrence patterns and highlight functional differences of lake phytoplankton communities.

Main conclusions

Diverging community trajectories have led to a nationwide reshuffling of lake phytoplankton communities. At regional and national scales, lakes are not single entities but metacommunity hubs in an interconnected waterscape. The assembly mechanisms of phytoplankton communities are strongly structured by spatio-temporal dynamics, which have led to novel community types, but only a minor part of this reshuffling could be linked to temporal environmental change.

Factors at multiple scales drive parasite community structure

Understanding how ecological communities are assembled remains a key goal of ecosystem ecology. Because communities are hierarchical, factors acting at multiple scales can contribute to patterns of community structure. Parasites provide a natural system to explore this idea, as they exist as discrete communities within host individuals, which are themselves part of a community and metacommunity. We aimed to understand the relative contribution of multi-scale drivers in parasite community assembly and assess how patterns at one level may mask those occurring at another. Specifically, we wanted to disentangle patterns caused by passive sampling from those determined by ecological drivers, and how these vary with scale. We applied a Markov Random Fields model and assessed measures of β-diversity and nestedness for 420 replicate parasite infracommunities (parasite assemblages in host individuals) across two freshwater mussel host species, three sites and two time periods, comparing our results to simulations from four different ecologically relevant null models. We showed that β-diversity between sites (explaining 25% of variation in parasite distribution) and host species (41%) is greater than expected, and β-diversity between individual hosts is smaller than expected, even after accounting for parasite prevalence and characteristics of host individuals. Furthermore, parasite communities were significantly less nested than expected once parasite prevalence and host characteristics were both accounted for, but more nested than expected otherwise, suggesting a degree of modularity at the within-host level that is masked if underlying host and parasite characteristics are not taken into account. The Markov Random Fields model provided evidence for possible competitive within-host parasite interactions, providing a mechanism for the observed infracommunity modularity. An integrative approach that examines factors at multiple scales is necessary to understand the composition of ecological communities. Furthermore, patterns at one level can alter the interpretation of ecologically important drivers at another if variation at higher scales is not accounted for.

Ensemble modeling of black pomfret (Parastromateus niger) habitat in the Taiwan Strait based on oceanographic variables

The location, effort, number of captures, and time of fishing were all used in this study to assess the geographic distribution of Parastromateus niger in the Taiwan Strait. Other species distribution models performed worse than generalized linear models (GLMs) based on six oceanographic parameters. The sea surface temperature (SST) was between 26.5 °C and 29.5 °C, the sea surface chlorophyll (SSC) level was between 0.3-0.44 mg/m³, the sea surface salinity (SSS) was between 33.4 °C and 34.4 °C, the mixed layer depth was between 10 °C and 14 °C, the sea surface height was between 0.57 °C and 0.77 °C, and the eddy kinetic energy (EKE) was between 0.603 °C. According to the statistical findings, SST is merely a small effect compared to SSS, SSC level, and EKE in terms of impacting species distribution. By combining four effective single-algorithm models with no obvious bias, an ensemble habitat model was created. The ranges of 117°E-119°E and 22°N-24°N have the highest annual distributions of S.CPUE and nominal CPUE.

The community ecology perspective of omics data

The measurement of uncharacterized pools of biological molecules through techniques such as metabarcoding, metagenomics, metatranscriptomics, metabolomics, and metaproteomics produces large, multivariate datasets. Analyses of these datasets have successfully been borrowed from community ecology to characterize the molecular diversity of samples (ɑ-diversity) and to assess how these profiles change in response to experimental treatments or across gradients (β-diversity). However, sample preparation and data collection methods generate biases and noise which confound molecular diversity estimates and require special attention. Here, we examine how technical biases and noise that are introduced into multivariate molecular data affect the estimation of the components of diversity (i.e., total number of different molecular species, or entities; total number of molecules; and the abundance distribution of molecular entities). We then explore under which conditions these biases affect the measurement of ɑ- and β-diversity and highlight how novel methods commonly used in community ecology can be adopted to improve the interpretation and integration of multivariate molecular data. Video Abstract.

Elevational changes in the bacterial community composition and potential functions in a Tibetan grassland

In the Tibetan grasslands, the distribution patterns of the microbial community structure and function along elevation gradients have attracted considerable attention due to the wide distribution of mountain slopes, but the controlling factors of these patterns are still unclear. Here we investigated the taxonomy and potential functions of soil bacteria along an elevation gradient in a Tibetan mountainous grassland in 2 years, aiming to explore the elevation patterns of the bacterial structure and function and the underlying drivers. High-throughput sequencing and environment attribute measurements were conducted to survey the bacterial and environment characters. Furthermore, PICRUSt2 for prediction of bacterial functions, iCAMP for unraveling the drivers controlling community assembly, and HMSC for variance partitioning of bacterial community composition were used. Elevation did not significantly affect the bacterial diversity but changed their composition, driven by both deterministic and stochastic processes. In addition, elevation did not significantly affect the relative importance of deterministic and stochastic processes. Soil carbon, nitrogen, and temperature were important deterministic factors in driving bacterial community structure. The genus Solirubrobacter in Actinobacteriota was identified as most elevation discriminatory. Based on these observations, the bacterial community in the Tibetan mountainous grasslands was more controlled by edaphic factors than temperature, indicating their relative stability under climate change scenarios.

Among- and within-population variation in morphology, rewards, and scent in a hawkmoth-pollinated plant

PREMISE

Floral scent is a complex trait that mediates many plant-insect interactions, but our understanding of how floral scent variation evolves, either independently or in concert with other traits, remains limited. Assessing variation in floral scent at multiple levels of biological organization and comparing patterns of variation in scent to variation in other floral traits can contribute to our understanding of how scent variation evolves in nature.

METHODS

We used a greenhouse common garden experiment to investigate variation in floral scent at three scales-within plants, among plants, and among populations-and to determine whether scent, alone or in combination with morphology and rewards, contributes to population differentiation in Oenothera cespitosa subsp. marginata. Its range spans most of the biomes in the western United States, such that variation in both the abiotic and biotic environment could contribute to trait variation.

RESULTS

Multiple analytical approaches demonstrated substantial variation among and within populations in compound-specific and total floral scent measures. Overall, populations were differentiated in morphology and reward traits and in scent. Across populations, coupled patterns of variation in linalool, leucine-derived compounds, and hypanthium length are consistent with a long-tongued moth pollination syndrome.

CONCLUSIONS

The considerable variation in floral scent detected within populations suggests that, similar to other floral traits, variation in floral scent may have a heritable genetic component. Differences in patterns of population differentiation in floral scent and in morphology and rewards indicate that these traits may be shaped by different selective pressures.

Measuring, comparing and interpreting phenotypic selection on floral scent

Natural selection on floral scent composition is a key element of the hypothesis that pollinators and other floral visitors drive scent evolution. The measure of such selection is complicated by the high-dimensional nature of floral scent data and uncertainty about the cognitive processes involved in scent-mediated communication. We use dimension reduction through reduced-rank regression to jointly estimate a scent composite trait under selection and the strength of selection acting on this trait. To assess and compare variation in selection on scent across species, time and space, we reanalyse 22 datasets on six species from four previous studies. The results agreed qualitatively with previous analyses in terms of identifying populations and scent compounds subject to stronger selection but also allowed us to evaluate and compare the strength of selection on scent across studies. Doing so revealed that selection on floral scent was highly variable, and overall about as common and as strong as selection on other phenotypic traits involved in pollinator attraction or pollen transfer. These results are consistent with an important role of floral scent in pollinator attraction. Our approach should be useful for further studies of plant-animal communication and for studies of selection on other high-dimensional phenotypes. In particular, our approach will be useful for studies of pollinator-mediated selection on complex scent blends comprising many volatiles, and when no prior information on the physiological responses of pollinators to scent compounds is available.

Generalized Matrix Factorization: efficient algorithms for fitting generalized linear latent variable models to large data arrays

Unmeasured or latent variables are often the cause of correlations between multivariate measurements, which are studied in a variety of fields such as psychology, ecology, and medicine. For Gaussian measurements, there are classical tools such as factor analysis or principal component analysis with a well-established theory and fast algorithms. Generalized Linear Latent Variable models (GLLVMs) generalize such factor models to non-Gaussian responses. However, current algorithms for estimating model parameters in GLLVMs require intensive computation and do not scale to large datasets with thousands of observational units or responses. In this article, we propose a new approach for fitting GLLVMs to high-dimensional datasets, based on approximating the model using penalized quasi-likelihood and then using a Newton method and Fisher scoring to learn the model parameters. Computationally, our method is noticeably faster and more stable, enabling GLLVM fits to much larger matrices than previously possible. We apply our method on a dataset of 48,000 observational units with over 2,000 observed species in each unit and find that most of the variability can be explained with a handful of factors. We publish an easy-to-use implementation of our proposed fitting algorithm.

Hydropeaking impacts on riverine plants downstream from the world's largest hydropower dam, the Three Gorges Dam

Hydropeaking has become a global issue because of extensive hydropower dam construction worldwide. Yet, its ecological impacts on riverine ecosystems are not well studied. We explored the effects of hydropeaking on riverine plants, based on data from a ~300 km reach downstream of the world's largest hydropower dam, the Three Gorges Dam. We tested three hypotheses relating to hydropeaking impacts on species elevational distributions, assemblage structure and species-specific biomass patterns by generalized linear mixed modelling and joint species distribution modelling. We found that, first, hydropeaking greatly shaped species elevational ranges, leading to expansions of herbs to high elevations and shifting species dominance at low elevations. Secondly, we detected contrasting effects of hydropeaking on assemblage-level characteristics of herbs. The inundation induced by hydropeaking had strong effects on assemblage composition and biomass allocation, where more biomass was allocated to belowground part. Hydropeaking blurred the species richness-biomass relationship, although it had little effect on species richness or plot-level biomass. Thirdly, hydropeaking induced inundation was the most important covariate driving species biomass patterns of riverine plants, although complex species-specific effects were identified, and random effects were often large in fitted models. We concluded that hydropeaking likely acted as a major driver of plant community assembly in rivers with a hydropower dam. Conservation and restoration of riverine plants can benefit from the inclusion of water level management in operational schemes of hydropower dams, especially during the early life history stages.

Tick microbial associations at the crossroad of horizontal and vertical transmission pathways

BACKGROUND

Microbial communities can affect disease risk by interfering with the transmission or maintenance of pathogens in blood-feeding arthropods. Here, we investigated whether bacterial communities vary between Ixodes ricinus nymphs which were or were not infected with horizontally transmitted human pathogens.

METHODS

Ticks from eight forest sites were tested for the presence of Borrelia burgdorferi sensu lato, Babesia spp., Anaplasma phagocytophilum, and Neoehrlichia mikurensis by quantitative polymerase chain reaction (qPCR), and their microbiomes were determined by 16S rRNA amplicon sequencing. Tick bacterial communities clustered poorly by pathogen infection status but better by geography. As a second approach, we analysed variation in tick microorganism community structure (in terms of species co-infection) across space using hierarchical modelling of species communities. For that, we analysed almost 14,000 nymphs, which were tested for the presence of horizontally transmitted pathogens B. burgdorferi s.l., A. phagocytophilum, and N. mikurensis, and the vertically transmitted tick symbionts Rickettsia helvetica, Rickettsiella spp., Spiroplasma ixodetis, and Candidatus Midichloria mitochondrii.

RESULTS

With the exception of Rickettsiella spp., all microorganisms had either significant negative (R. helvetica and A. phagocytophilum) or positive (S. ixodetis, N. mikurensis, and B. burgdorferi s.l.) associations with M. mitochondrii. Two tick symbionts, R. helvetica and S. ixodetis, were negatively associated with each other. As expected, both B. burgdorferi s.l. and N. mikurensis had a significant positive association with each other and a negative association with A. phagocytophilum. Although these few specific associations do not appear to have a large effect on the entire microbiome composition, they can still be relevant for tick-borne pathogen dynamics.

CONCLUSIONS

Based on our results, we propose that M. mitochondrii alters the propensity of ticks to acquire or maintain horizontally acquired pathogens. The underlying mechanisms for some of these remarkable interactions are discussed herein and merit further investigation. Positive and negative associations between and within horizontally and vertically transmitted symbionts.

Skin microbiota of oxazolone-induced contact hypersensitivity mouse model

Contact allergy is a common skin allergy, which can be studied utilising contact hypersensitivity (CHS) animal model. However, it is not clear, whether CHS is a suitable model to investigate skin microbiota interactions. We characterised the effect of contact dermatitis on the skin microbiota and studied the biological effects of oxazolone (OXA) -induced inflammation on skin thickness, immune cell numbers and changes of the microbiota in CHS mouse model (n = 72) for 28 days. Through 16S rRNA gene sequencing we defined the composition of bacterial communities and associations of bacteria with inflammation. We observed that the vehicle solution of acetone and olive oil induced bacterial community changes on day 1, and OXA-induced changes were observed mainly on day 7. Many of the notably enriched bacteria present in the OXA-challenged positive group represented the genus Faecalibaculum which were most likely derived from the cage environment. Additionally, skin inflammation correlated negatively with Streptococcus, which is considered a native skin bacterium, and positively with Muribacter muris, which is typical in oral environment. Skin inflammation favoured colonisation of cage-derived faecal bacteria, and additionally mouse grooming transferred oral bacteria on the skin. Due to the observed changes, we conclude that CHS model could be used for certain skin microbiome-related research set-ups. However, since vehicle exposure can alter the skin microbiome as such, future studies should include considerations such as careful control sampling and statistical tests to account for potential confounding factors.

Beyond Basic Diversity Estimates-Analytical Tools for Mechanistic Interpretations of Amplicon Sequencing Data

Understanding microbial ecology through amplifying short read regions, typically 16S rRNA for prokaryotic species or 18S rRNA for eukaryotic species, remains a popular, economical choice. These methods provide relative abundances of key microbial taxa, which, depending on the experimental design, can be used to infer mechanistic ecological underpinnings. In this review, we discuss recent advancements in in situ analytical tools that have the power to elucidate ecological phenomena, unveil the metabolic potential of microbial communities, identify complex multidimensional interactions between species, and compare stability and complexity under different conditions. Additionally, we highlight methods that incorporate various modalities and additional information, which in combination with abundance data, can help us understand how microbial communities respond to change in a typical ecosystem. Whilst the field of microbial informatics continues to progress substantially, our emphasis is on popular methods that are applicable to a broad range of study designs. The application of these methods can increase our mechanistic understanding of the ongoing dynamics of complex microbial communities.

Enriching captivity conditions with natural elements does not prevent the loss of wild‐like gut microbiota but shapes its compositional variation in two small mammals

As continued growth in gut microbiota studies in captive and model animals elucidates the importance of their role in host biology, further pursuit of how to retain a wild‐like microbial community is becoming increasingly important to obtain representative results from captive animals. In this study, we assessed how the gut microbiota of two wild‐caught small mammals, namely Crocidura russula (Eulipotyphla, insectivore) and Apodemus sylvaticus (Rodentia, omnivore), changed when bringing them into captivity. We analyzed fecal samples of 15 A. sylvaticus and 21 C. russula, immediately after bringing them into captivity and 5 weeks later, spread over two housing treatments: a “natural” setup enriched with elements freshly collected from nature and a “laboratory” setup with sterile artificial elements. Through sequencing of the V3–V4 region of the 16S recombinant RNA gene, we found that the initial microbial diversity dropped during captivity in both species, regardless of treatment. Community composition underwent a change of similar magnitude in both species and under both treatments. However, we did observe that the temporal development of the gut microbiome took different trajectories (i.e., changed in different directions) under different treatments, particularly in C. russula, suggesting that C. russula may be more susceptible to environmental change. The results of this experiment do not support the use of microbially enriched environments to retain wild‐like microbial diversities and compositions, yet show that specific housing conditions can significantly affect the drift of microbial communities under captivity.

Early impacts of the largest Amazonian hydropower project on fish communities

Hydropower is a threat to freshwater fishes. Despite a recent boom in dam construction, few studies have assessed their impact on mega-diverse tropical rivers. Using a before-after study design, we investigated the early impacts of the Belo Monte hydroelectric complex, the third-largest hydropower project in the world, on fishes of the Xingu River, a major clear-water tributary of the lower Amazon. We explored impacts across different river sectors (upstream, reservoir, reduced flow sector, and downstream) and spatial scales (individual sectors vs. all sectors combined) using joint species distribution models and different facets of diversity (taxonomic, functional, and phylogenetic). After 5 years of the Belo Monte operation, species richness declined ~12% in lentic and ~16% in lotic environments. Changes in abundance were of less magnitude (<4%). Effects were particularly negative for species of the families Serrasalmidae (mainly pacus), Anostomidae (headstanders), Auchenipteridae, and Pimelodidae (catfishes), whereas no taxonomic group consistently increased in richness or abundance. The reservoir and downstream sectors were the most impacted, with declines of ~24-29% in fish species richness, overall reductions in fish body size and trophic level, and a change in average body shape. Richness and abundance also declined in the reduced river flow, and changes in size, shape, and position of fins were observed. Relatively minor changes were found in the upstream sector. Variation in functional and phylogenetic diversity following river impoundment was subtle; however, across sectors, we found a reduction in functional divergence, indicating a decline in the abundance of species located near the extremities of community functional space. This may be the first sign of an environmental filtering process reducing functional diversity in the region. Greater changes in flow and habitats are expected as hydropower operations ramp up, and continued monitoring is warranted to understand the full scope and magnitude of ecological impacts.