High prevalence and concomitant infection of Ranavirus and Eustrongylides sp. in the invasive American Bullfrog in Brazil

American Bullfrogs, Aquarana catesbeiana, are invasive anuran species distributed worldwide. One of the adverse impacts that this species causes in native communities is as a reservoir host for pathogens and parasites. Here, we report the coinfection of two pathogenic organisms in A. catesbeiana: Ranavirus and the nematode Eustrongylides. Bullfrogs were collected in the wild in a pond close to the urban area of São Paulo, Brazil. The prevalence of both pathogens was high: 77% were infected with ranavirus with a mean viral load of 1010.3 viral copies, and 100% of the bullfrogs were infected by Eustrongylides sp. with a mean intensity of infection of 13.4 specimens per host. Four host specimens (31%) presented pathological signs that seemed to be related to the Eustrongylides sp. infection, such as internal organs adhered to each other due to high intensity and large size of the nematodes, ulcers, and raw flesh wounds caused by the nematode. The pathogenic and concomitant infections have potential zoonotic implications and raise concerns about human infection risks for Eustrongylides infection. Moreover, such infections may represent an additional level of threat to native communities through the potential shifts in patterns of parasite and pathogen transmission. Future research involving the native anuran community is essential to ascertain whether invasive bullfrogs are attenuating or exacerbating diseases such as ranavirosis and eustrongylidiosis.

Ranavirus and helminth parasite co-infection in invasive American bullfrogs in the Atlantic forest, Brazil

Emerging infectious diseases threaten amphibian species across the globe. In Brazil, the American bullfrog (Aquarana catesbeiana) is a highly invasive species that can potentially transmit parasites and pathogens to native amphibians. This is the first assessment of co-infection of Ranavirus and helminth macroparasites in invasive populations of bullfrogs in South America. We collected, measured, and euthanized 65 specimens of A. catesbeiana sampled from 9 sites across three states of Brazil in the Atlantic Forest biome. We collected and identified helminth macroparasites and sampled host liver tissue to test for the presence and load of Ranavirus with quantitative PCR. We documented patterns of prevalence, parasite load, and co-infection with generalized linear mixed models, generalized logistic regressions, and randomization tests. Most individual bullfrogs did not exhibit clinical signs of infection, but the overall Ranavirus prevalence was 27% (95% confidence interval, [CI 17-38]). Bullfrogs were infected with helminth macroparasites from 5 taxa. Co-infection of helminth macroparasites and Ranavirus was also common (21% CI [12-31]). Bullfrog size was positively correlated with total macroparasite abundance and richness, and the best-fitting model included a significant interaction between bullfrog size and Ranavirus infection status. We observed a negative correlation between Ranavirus viral load and nematode abundance (slope = -0.22, P = 0.03). Invasive bullfrogs (A. catesbeiana) in Brazil were frequently infected with both Ranavirus and helminth macroparasites, so adult bullfrogs could serve as reservoir hosts for both pathogens and parasites. However, many macroparasites collected were encysted and not developing. Coinfection patterns suggest a potential interaction between Ranavirus and macroparasites because helminth abundance increased with bullfrog size but was lower in Ranavirus infected individuals. Future studies of bullfrogs in the Atlantic Forest should investigate their potential role in pathogen and parasite transmission to native anurans.

Synthesis reveals approximately balanced biotic differentiation and homogenization

It is commonly thought that the biodiversity crisis includes widespread declines in the spatial variation of species composition, called biotic homogenization. Using a typology relating homogenization and differentiation to local and regional diversity changes, we synthesize patterns across 461 metacommunities surveyed for 10 to 91 years, and 64 species checklists (13 to 500+ years). Across all datasets, we found that no change was the most common outcome, but with many instances of homogenization and differentiation. A weak homogenizing trend of a 0.3% increase in species shared among communities/year on average was driven by increased numbers of widespread (high occupancy) species and strongly associated with checklist data that have longer durations and large spatial scales. At smaller spatial and temporal scales, we show that homogenization and differentiation can be driven by changes in the number and spatial distributions of both rare and common species. The multiscale perspective introduced here can help identify scale-dependent drivers underpinning biotic differentiation and homogenization.

Looking back on biodiversity change: lessons for the road ahead

Estimating biodiversity change across the planet in the context of widespread human modification is a critical challenge. Here, we review how biodiversity has changed in recent decades across scales and taxonomic groups, focusing on four diversity metrics: species richness, temporal turnover, spatial beta-diversity and abundance. At local scales, change across all metrics includes many examples of both increases and declines and tends to be centred around zero, but with higher prevalence of declining trends in beta-diversity (increasing similarity in composition across space or biotic homogenization) and abundance. The exception to this pattern is temporal turnover, with changes in species composition through time observed in most local assemblages. Less is known about change at regional scales, although several studies suggest that increases in richness are more prevalent than declines. Change at the global scale is the hardest to estimate accurately, but most studies suggest extinction rates are probably outpacing speciation rates, although both are elevated. Recognizing this variability is essential to accurately portray how biodiversity change is unfolding, and highlights how much remains unknown about the magnitude and direction of multiple biodiversity metrics at different scales. Reducing these blind spots is essential to allow appropriate management actions to be deployed. This article is part of the theme issue 'Detecting and attributing the causes of biodiversity change: needs, gaps and solutions'.

Time-travelling pathogens and their risk to ecological communities

Permafrost thawing and the potential 'lab leak' of ancient microorganisms generate risks of biological invasions for today's ecological communities, including threats to human health via exposure to emergent pathogens. Whether and how such 'time-travelling' invaders could establish in modern communities is unclear, and existing data are too scarce to test hypotheses. To quantify the risks of time-travelling invasions, we isolated digital virus-like pathogens from the past records of coevolved artificial life communities and studied their simulated invasion into future states of the community. We then investigated how invasions affected diversity of the free-living bacteria-like organisms (i.e., hosts) in recipient communities compared to controls where no invasion occurred (and control invasions of contemporary pathogens). Invading pathogens could often survive and continue evolving, and in a few cases (3.1%) became exceptionally dominant in the invaded community. Even so, invaders often had negligible effects on the invaded community composition; however, in a few, highly unpredictable cases (1.1%), invaders precipitated either substantial losses (up to -32%) or gains (up to +12%) in the total richness of free-living species compared to controls. Given the sheer abundance of ancient microorganisms regularly released into modern communities, such a low probability of outbreak events still presents substantial risks. Our findings therefore suggest that unpredictable threats so far confined to science fiction and conjecture could in fact be powerful drivers of ecological change.

Six decades of museum collections reveal disruption of native ant assemblages by introduced species

There is a looming environmental crisis characterized by widespread declines in global biodiversity,^1,2,3,4,5,6 coupled with the establishment of introduced species at accelerated rates.^{7,8,9,10,11,12,13,14} We quantified how multi-species invasions affect litter ant communities in natural ecosystems by leveraging museum records and contemporary collections to assemble a large (18,990 occurrences, 6,483 sampled local communities, and 177 species) 54-year (1965-2019) dataset for the entire state of Florida, USA. Nine of ten species that decreased most strongly in relative abundance ("losers") were native, while nine of the top ten "winners" were introduced species. These changes led to shifts in the composition of rare and common species: in 1965, only two of the ten most common ants were introduced, whereas by 2019, six of ten were introduced species. Native losers included seed dispersers and specialist predators, suggesting a potential loss of ecosystem function through time, despite no obvious loss of phylogenetic diversity. We also examined the role of species-level traits as predictors of invasion success. Introduced species were more likely to be polygynous than native species. The tendency to form supercolonies, where workers from separate nests integrate, also differed between native and introduced species and was correlated with the degree to which species increased in their rank abundances over 50 years. In Florida, introduced ants now account for 30% of occurrence records, and up to 70% in southern Florida. If current trends continue, introduced species will account for over half of occurrence records in all Florida's litter ant communities within the next 50 years.

Regional occupancy increases for widespread species but decreases for narrowly distributed species in metacommunity time series

While human activities are known to elicit rapid turnover in species composition through time, the properties of the species that increase or decrease their spatial occupancy underlying this turnover are less clear. Here, we used an extensive dataset of 238 metacommunity time series of multiple taxa spread across the globe to evaluate whether species that are more widespread (large-ranged species) differed in how they changed their site occupancy over the 10-90 years the metacommunities were monitored relative to species that are more narrowly distributed (small-ranged species). We found that on average, large-ranged species tended to increase in occupancy through time, whereas small-ranged species tended to decrease. These relationships were stronger in marine than in terrestrial and freshwater realms. However, in terrestrial regions, the directional changes in occupancy were less extreme in protected areas. Our findings provide evidence for systematic decreases in occupancy of small-ranged species, and that habitat protection could mitigate these losses in the face of environmental change.

How does variation in total and relative abundance contribute to gradients of species diversity?

Patterns of biodiversity provide insights into the processes that shape biological communities around the world. Variation in species diversity along biogeographical or ecological gradients, such as latitude or precipitation, can be attributed to variation in different components of biodiversity: changes in the total abundance (i.e., more-individual effects) and changes in the regional species abundance distribution (SAD). Rarefaction curves can provide a tool to partition these sources of variation on diversity, but first must be converted to a common unit of measurement. Here, we partition species diversity gradients into components of the SAD and abundance using the effective number of species (ENS) transformation of the individual-based rarefaction curve. Because the ENS curve is unconstrained by sample size, it can act as a standardized unit of measurement when comparing effect sizes among different components of biodiversity change. We illustrate the utility of the approach using two data sets spanning latitudinal diversity gradients in trees and marine reef fish and find contrasting results. Whereas the diversity gradient of fish was mostly associated with variation in abundance (86%), the tree diversity gradient was mostly associated with variation in the SAD (59%). These results suggest that local fish diversity may be limited by energy through the more-individuals effect, while species pool effects are the larger determinant of tree diversity. We suggest that the framework of the ENS-curve has the potential to quantify the underlying factors influencing most aspects of diversity change.

Source-sink behavioural dynamics limit institutional evolution in a group-structured society

Social change in any society entails changes in both behaviours and institutions. We model a group-structured society in which the transmission of individual behaviour occurs in parallel with the selection of group-level institutions. We consider a cooperative behaviour that generates collective benefits for groups but does not spread between individuals on its own. Groups exhibit institutions that increase the diffusion of the behaviour within the group, but also incur a group cost. Groups adopt institutions in proportion to their fitness. Finally, the behaviour may also spread globally. We find that behaviour and institutions can be mutually reinforcing. But the model also generates behavioural source-sink dynamics when behaviour generated in institutionalized groups spreads to non-institutionalized groups and boosts their fitness. Consequently, the global diffusion of group-beneficial behaviour creates a pattern of institutional free-riding that limits the evolution of group-beneficial institutions. Our model suggests that, in a group-structured society, large-scale beneficial social change can be best achieved when the relevant behaviour and institutions remain correlated.

Random placement models explain species richness and dissimilarity of frog assemblages within Atlantic Forest fragments

Understanding the effects of random versus niche-based processes on biodiversity patterns is a central theme in ecology, and an important tool for predicting effects of habitat loss and fragmentation on biodiversity. We investigated the predictive power of random processes to explain species richness and species dissimilarity of amphibian assemblages in a fragmented tropical landscape of the Atlantic Forest of South America. We analyzed a large database of amphibian abundance and occupancy, sampled in 21 forest fragments ranging in size from 1.9 to 619 ha. We compared observed species richness and species dissimilarity with the outcomes of two null (random placement) models: 1- the traditional Coleman's area-based model and 2- an abundance-based model (based on the number of individuals observed in each fragment). We applied these models for all species combined, and separately for forest-dependent and habitat-generalist species. The abundance-based model fitted the observed species richness data better than the area-based model for all species, forest-dependent species, and generalist species. The area-based and the abundance-based models were also able to significantly explain species dissimilarity for all species and for generalists, but not for forest dependent species. The traditional area-based model assigned too many individuals to large fragments, thus failing to accurately explain species richness within patches across the landscape. Although niche-based processes may be important to structuring the regional pool of species in fragmented landscapes, our results suggest that part of the variation in species richness and species dissimilarity can be successfully explained by random placement models, especially for generalist species. Evaluating which factors cause variation in the number of individuals among patches should be a focus in future studies aiming to understand biodiversity patterns in fragmented landscapes.

Long-term changes in temperate marine fish assemblages are driven by a small subset of species

The species composition of plant and animal assemblages across the globe has changed substantially over the past century. How do the dynamics of individual species cause this change? We classified species into seven unique categories of temporal dynamics based on the ordered sequence of presences and absences that each species contributes to an assemblage time series. We applied this framework to 14,434 species trajectories comprising 280 assemblages of temperate marine fishes surveyed annually for 20 or more years. Although 90% of the assemblages diverged in species composition from the baseline year, this compositional change was largely driven by only 8% of the species' trajectories. Quantifying the reorganization of assemblages based on species shared temporal dynamics should facilitate the task of monitoring and restoring biodiversity. We suggest ways in which our framework could provide informative measures of compositional change, as well as leverage future research on pattern and process in ecological systems.

Environment-host-microbial interactions shape the Sarracenia purpurea microbiome at the continental scale

The importance of climate, habitat structure, and higher trophic levels on microbial diversity is only beginning to be understood. Here, we examined the influence of climate variables, plant morphology, and the abundance of aquatic invertebrates on the microbial biodiversity of the northern pitcher plant Sarracenia purpurea. The plant's cup-shaped leaves fill with rainwater and support a miniature, yet full-fledged, ecosystem with a diverse microbiome that decomposes captured prey and a small network of shredding and filter-feeding aquatic invertebrates that feed on microbes. We characterized pitcher microbiomes of 108 plants sampled at 36 sites from Florida to Quebec. Structural equation models revealed that annual precipitation and temperature, plant size, and midge abundance had direct effects on microbiome taxonomic and phylogenetic diversity. Climate variables also exerted indirect effects through plant size and midge abundance. Further, spatial structure and climate influenced taxonomic composition, but not phylogenetic composition. Our results suggest that direct effects of midge abundance and climate and indirect effects of climate through its effect on plant-associated factors lead to greater richness of microbial phylotypes in warmer, wetter sites.

A multiscale framework for disentangling the roles of evenness, density, and aggregation on diversity gradients

Disentangling the drivers of diversity gradients can be challenging. The Measurement of Biodiversity (MoB) framework decomposes scale‐dependent changes in species diversity into three components of community structure: species abundance distribution (SAD), total community abundance, and within‐species spatial aggregation. Here we extend MoB from categorical treatment comparisons to quantify variation along continuous geographic or environmental gradients. Our approach requires sites along a gradient, each consisting of georeferenced plots of abundance‐based species composition data. We demonstrate our method using a case study of ants sampled along an elevational gradient of 28 sites in a mixed deciduous forest of the Great Smoky Mountains National Park, USA. MoB analysis revealed that decreases in ant species richness along the elevational gradient were associated with decreasing evenness and total number of species, which counteracted the modest increase in richness associated with decreasing spatial aggregation along the gradient. Total community abundance had a negligible effect on richness at all but the finest spatial grains, SAD effects increased in importance with sampling effort, and the aggregation effect had the strongest effect at coarser spatial grains. These results do not support the more‐individuals hypothesis, but they are consistent with a hypothesis of stronger environmental filtering at coarser spatial grains. Our extension of MoB has the potential to elucidate how components of community structure contribute to changes in diversity along environmental gradients and should be useful for a variety of assemblage‐level data collected along gradients.

Using Climatic Credits to Pay the Climatic Debt

Many organisms are accumulating climatic debt as they respond more slowly than expected to rising global temperatures, leading to disequilibrium of species diversity with contemporary climate. The resulting transient dynamics are complex and may cause overoptimistic biodiversity assessments. We propose a simple budget framework to integrate climatic debt with two classes of intervention: (i) climatic credits that pay some of the debt, reducing the overall biological change required to reach a new equilibrium; and (ii) options to adjust the debt repayment rate, either making a system more responsive by increasing the rate or temporarily reducing the rate to buy more time for local adaptation and credit implementation. We illustrate how this budget can be created and highlight limitations and challenges.

Clockwise and counterclockwise hysteresis characterize state changes in the same aquatic ecosystem

Incremental increases in a driver variable, such as nutrients or detritus, can trigger abrupt shifts in aquatic ecosystems that may exhibit hysteretic dynamics and a slow return to the initial state. A model system for understanding these dynamics is the microbial assemblage that inhabits the cup-shaped leaves of the pitcher plant Sarracenia purpurea. With enrichment of organic matter, this system flips within three days from an oxygen-rich state to an oxygen-poor state. In a replicated greenhouse experiment, we enriched pitcher-plant leaves at different rates with bovine serum albumin (BSA), a molecular substitute for detritus. Changes in dissolved oxygen (DO) and undigested BSA concentration were monitored during enrichment and recovery phases. With increasing enrichment rates, the dynamics ranged from clockwise hysteresis (low), to environmental tracking (medium), to novel counter-clockwise hysteresis (high). These experiments demonstrate that detrital enrichment rate can modulate a diversity of hysteretic responses within a single aquatic ecosystem, and suggest different management strategies may be needed to mitigate the effects of high vs. low rates of detrital enrichment.

Investigating Biotic Interactions in Deep Time

Recent renewed interest in using fossil data to understand how biotic interactions have shaped the evolution of life is challenging the widely held assumption that long-term climate changes are the primary drivers of biodiversity change. New approaches go beyond traditional richness and co-occurrence studies to explicitly model biotic interactions using data on fossil and modern biodiversity. Important developments in three primary areas of research include analysis of (i) macroevolutionary rates, (ii) the impacts of and recovery from extinction events, and (iii) how humans (Homo sapiens) affected interactions among non-human species. We present multiple lines of evidence for an important and measurable role of biotic interactions in shaping the evolution of communities and lineages on long timescales.

Reorganization of surviving mammal communities after the end-Pleistocene megafaunal extinction

Large mammals are at high risk of extinction globally. To understand the consequences of their demise for community assembly, we tracked community structure through the end-Pleistocene megafaunal extinction in North America. We decomposed the effects of biotic and abiotic factors by analyzing co-occurrence within the mutual ranges of species pairs. Although shifting climate drove an increase in niche overlap, co-occurrence decreased, signaling shifts in biotic interactions. Furthermore, the effect of abiotic factors on co-occurrence remained constant over time while the effect of biotic factors decreased. Biotic factors apparently played a key role in continental-scale community assembly before the extinctions. Specifically, large mammals likely promoted co-occurrence in the Pleistocene, and their loss contributed to the modern assembly pattern in which co-occurrence frequently falls below random expectations.

Trade-Offs in Cold Resistance at the Northern Range Edge of the Common Woodland Ant Aphaenogaster picea (Formicidae)

Geographic variation in low temperatures at poleward range margins of terrestrial species often mirrors population variation in cold resistance, suggesting that range boundaries may be set by evolutionary constraints on cold physiology. The northeastern woodland ant Aphaenogaster picea occurs up to approximately 45°N in central Maine. We combined presence/absence surveys with classification tree analysis to characterize its northern range limit and assayed two measures of cold resistance operating on different timescales to determine whether and how marginal populations adapt to environmental extremes. The range boundary of A. picea was predicted primarily by temperature, but low winter temperatures did not emerge as the primary correlate of species occurrence. Low summer temperatures and high seasonal variability predicted absence above the boundary, whereas high mean annual temperature (MAT) predicted presence in southern Maine. In contrast, assays of cold resistance across multiple sites were consistent with the hypothesis of local cold adaptation at the range edge: among populations, there was a 4-min reduction in chill coma recovery time across a 2° reduction in MAT. Baseline resistance and capacity for additional plastic cold hardening shifted in opposite directions, with hardening capacity approaching zero at the coldest sites. This trade-off between baseline resistance and cold-hardening capacity suggests that populations at range edges may adapt to colder temperatures through genetic assimilation of plastic responses, potentially constraining further adaptation and range expansion.

Proportional mixture of two rarefaction/extrapolation curves to forecast biodiversity changes under landscape transformation

Progressive habitat transformation causes global changes in landscape biodiversity patterns, but can be hard to quantify. Rarefaction/extrapolation approaches can quantify within-habitat biodiversity, but may not be useful for cases in which one habitat type is progressively transformed into another habitat type. To quantify biodiversity patterns in such transformed landscapes, we use Hill numbers to analyse individual-based species abundance data or replicated, sample-based incidence data. Given biodiversity data from two distinct habitat types, when a specified proportion of original habitat is transformed, our approach utilises a proportional mixture of two within-habitat rarefaction/extrapolation curves to analytically predict biodiversity changes, with bootstrap confidence intervals to assess sampling uncertainty. We also derive analytic formulas for assessing species composition (i.e. the numbers of shared and unique species) for any mixture of the two habitat types. Our analytical and numerical analyses revealed that species unique to each habitat type are the most important determinants of landscape biodiversity patterns.

Water quality improvements offset the climatic debt for stream macroinvertebrates over twenty years

Many species are accumulating climatic debt as they fail to keep pace with increasing global temperatures. In theory, concomitant decreases in other stressors (e.g. pollution, fragmentation) could offset some warming effects, paying climatic debt with accrued environmental credit. This process may be occurring in many western European rivers. We fit a Markov chain model to ~20,000 macroinvertebrate samples from England and Wales, and demonstrate that despite large temperature increases 1991-2011, macroinvertebrate communities remained close to their predicted equilibrium with environmental conditions. Using a novel analysis of multiple stressors, an accumulated climatic debt of 0.64 (±0.13 standard error) °C of warming was paid by a water-quality credit equivalent to 0.89 (±0.04)°C of cooling. Although there is finite scope for mitigating additional climate warming in this way, water quality improvements appear to have offset recent temperature increases, and the concept of environmental credit may be a useful tool for communicating climate offsetting.

Diversity-disease relationships and shared species analyses for human microbiome-associated diseases

Diversity indices have been routinely computed in the study of human microbiome-associated diseases (MADs). However, it is still unclear whether there is a consistent diversity-disease relationship (DDR) for the human MADs, and whether there are consistent differences in the taxonomic composition of microbiomes sampled from healthy versus diseased individuals. Here we reanalyzed raw data and used a meta-analysis to compare the microbiome diversity and composition of healthy versus diseased individuals in 41 comparisons extracted from 27 previously published studies of human MADs. In the DDR analysis, the average effect size across studies did not differ from zero for a comparison of healthy versus diseased individuals. In 30 of 41 comparisons (73%) there was no significant difference in microbiome diversity of healthy versus diseased individuals, or of different disease classes. For the species composition analysis (shared species analysis), the effect sizes were significantly different from zero. In 33 of 41 comparisons (80%), there were fewer OTUs (operational taxonomic units) shared between healthy and diseased individuals than expected by chance, but with 49% (20 of 41 comparisons) statistically significant. These results imply that the taxonomic composition of disease-associated microbiomes is often distinct from that of healthy individuals. Because species composition changes with disease state, some microbiome OTUs may serve as potential diagnostic indicators of disease. However, the overall species diversity of human microbiomes is not a reliable indicator of disease.

Ecological drift and competitive interactions predict unique patterns in temporal fluctuations of population size

Recent studies have highlighted the importance of higher-order competitive interactions in stabilizing population dynamics in multi-species communities. But how does the structure of competitive hierarchies affect population dynamics and extinction processes? We tackled this important question by using spatially explicit simulations of ecological drift (10 species in a homogeneous landscape of 64 patches) in which birth rates were influenced by interspecific competition. Specifically, we examined how transitive (linear pecking orders) and intransitive (pecking orders with loops) competitive hierarchies affected extinction rates and population dynamics in simulated communities through time. In comparison to a pure neutral model, an ecological drift model including transitive competition increased extinction rates, caused synchronous density-dependent population fluctuations, and generated a white-noise distribution of population sizes. In contrast, the drift model with intransitive competitive interactions decreased extinctions rates, caused asynchronous (compensatory) density-dependent population fluctuations, and generated a brown noise distribution of population sizes. We also explored the effect on community stability of more complex patterns of competitive interactions in which pairwise competitive relationships were assigned probabilistically. These probabilistic competition models also generated density-dependent trajectories and a brown noise distribution of population sizes. However, extinction rates and the degree of population synchrony were comparable to those observed in purely neutral communities. Collectively, our results confirm that intransitive competition has a strong and stabilizing effect on local populations in species-poor communities. This effect wanes with increasing species richness. Empirical assemblages characterized by brown spectral noise, density-dependent regulation, and asynchronous (compensatory) population fluctuations may indicate a signature of intransitive competitive interactions.

Draft Aphaenogaster genomes expand our view of ant genome size variation across climate gradients

Given the abundance, broad distribution, and diversity of roles that ants play in many ecosystems, they are an ideal group to serve as ecosystem indicators of climatic change. At present, only a few whole-genome sequences of ants are available (19 of >16,000 species), mostly from tropical and sub-tropical species. To address this limited sampling, we sequenced genomes of temperate-latitude species from the genus Aphaenogaster, a genus with important seed dispersers. In total, we sampled seven colonies of six species: Aphaenogaster ashmeadi, Aphaenogaster floridana, Aphaenogaster fulva, Aphaenogaster miamiana, Aphaenogaster picea, and Aphaenogaster rudis. The geographic ranges of these species collectively span eastern North America from southern Florida to southern Canada, which encompasses a latitudinal gradient in which many climatic variables are changing rapidly. For the six genomes, we assembled an average of 271,039 contigs into 47,337 scaffolds. The Aphaenogaster genomes displayed high levels of completeness with 96.1% to 97.6% of Hymenoptera BUSCOs completely represented, relative to currently sequenced ant genomes which ranged from 88.2% to 98.5%. Additionally, the mean genome size was 370.5 Mb, ranging from 310.3 to 429.7, which is comparable to that of other sequenced ant genomes (212.8-396.0 Mb) and flow cytometry estimates (210.7-690.4 Mb). In an analysis of currently sequenced ant genomes and the new Aphaenogaster sequences, we found that after controlling for both spatial autocorrelation and phylogenetics ant genome size was marginally correlated with sample site climate similarity. Of all examined climate variables, minimum temperature, and annual precipitation had the strongest correlations with genome size, with ants from locations with colder minimum temperatures and higher levels of precipitation having larger genomes. These results suggest that climate extremes could be a selective force acting on ant genomes and point to the need for more extensive sequencing of ant genomes.

Deciphering the enigma of undetected species, phylogenetic, and functional diversity based on Good-Turing theory

Estimating the species, phylogenetic, and functional diversity of a community is challenging because rare species are often undetected, even with intensive sampling. The Good-Turing frequency formula, originally developed for cryptography, estimates in an ecological context the true frequencies of rare species in a single assemblage based on an incomplete sample of individuals. Until now, this formula has never been used to estimate undetected species, phylogenetic, and functional diversity. Here, we first generalize the Good-Turing formula to incomplete sampling of two assemblages. The original formula and its two-assemblage generalization provide a novel and unified approach to notation, terminology, and estimation of undetected biological diversity. For species richness, the Good-Turing framework offers an intuitive way to derive the non-parametric estimators of the undetected species richness in a single assemblage, and of the undetected species shared between two assemblages. For phylogenetic diversity, the unified approach leads to an estimator of the undetected Faith's phylogenetic diversity (PD, the total length of undetected branches of a phylogenetic tree connecting all species), as well as a new estimator of undetected PD shared between two phylogenetic trees. For functional diversity based on species traits, the unified approach yields a new estimator of undetected Walker et al.'s functional attribute diversity (FAD, the total species-pairwise functional distance) in a single assemblage, as well as a new estimator of undetected FAD shared between two assemblages. Although some of the resulting estimators have been previously published (but derived with traditional mathematical inequalities), all taxonomic, phylogenetic, and functional diversity estimators are now derived under the same framework. All the derived estimators are theoretically lower bounds of the corresponding undetected diversities; our approach reveals the sufficient conditions under which the estimators are nearly unbiased, thus offering new insights. Simulation results are reported to numerically verify the performance of the derived estimators. We illustrate all estimators and assess their sampling uncertainty with an empirical dataset for Brazilian rain forest trees. These estimators should be widely applicable to many current problems in ecology, such as the effects of climate change on spatial and temporal beta diversity and the contribution of trait diversity to ecosystem multi-functionality.

Environmental proteomics reveals taxonomic and functional changes in an enriched aquatic ecosystem

Aquatic ecosystem enrichment can lead to distinct and irreversible changes to undesirable states. Understanding changes in active microbial community function and composition following organic-matter loading in enriched ecosystems can help identify biomarkers of such state changes. In a field experiment, we enriched replicate aquatic ecosystems in the pitchers of the northern pitcher plant, Sarracenia purpurea. Shotgun metaproteomics using a custom metagenomic database identified proteins, molecular pathways, and contributing microbial taxa that differentiated control ecosystems from those that were enriched. The number of microbial taxa contributing to protein expression was comparable between treatments; however, taxonomic evenness was higher in controls. Functionally active bacterial composition differed significantly among treatments and was more divergent in control pitchers than enriched pitchers. Aerobic and facultative anaerobic bacteria contributed most to identified proteins in control and enriched ecosystems, respectively. The molecular pathways and contributing taxa in enriched pitcher ecosystems were similar to those found in larger enriched aquatic ecosystems and are consistent with microbial processes occurring at the base of detrital food webs. Detectable differences between protein profiles of enriched and control ecosystems suggest that a time series of environmental proteomics data may identify protein biomarkers of impending state changes to enriched states.

Community-level regulation of temporal trends in biodiversity

Many theoretical models of community dynamics predict that species richness (S) and total abundance (N) are regulated in their temporal fluctuations. We present novel evidence for widespread regulation of biodiversity. For 59 plant and animal assemblages from around the globe monitored annually for a decade or more, the majority exhibited regulated fluctuations compared to the null hypothesis of an unconstrained random walk. However, there was little evidence for statistical artifacts, regulation driven by correlations with average annual temperature, or local-scale compensatory fluctuations in S or N. In the absence of major environmental perturbations, such as urbanization or cropland transformation, species richness and abundance may be buffered and exhibit some resilience in their temporal trajectories. These results suggest that regulatory processes are occurring despite unprecedented environmental change, highlighting the need for community-level assessment of biodiversity trends, as well as extensions of existing theory to address open source pools and shifting environmental conditions.

GEOGRAPHIC VARIATION IN LIFE-HISTORY TRAITS OF THE ANT LION, MYRMELEON IMMACULATUS: EVOLUTIONARY IMPLICATIONS OF BERGMANN'S RULE

In eastern North America, body size of the larval ant lion Myrmeleon immaculatus increases from south to north, following Bergmann's rule. We used a common-garden experiment and a reciprocal-transplant experiment to evaluate the effects of food and temperature on ant lion growth, body size, and survivorship. In the laboratory common-garden experiment, first-instar larvae from two southern (Georgia, South Carolina) and two northern (Connecticut, Rhode Island) populations were reared in incubators under high- and low-food and high- and low-temperature regimes. For all populations, high food increased final body mass and growth rate and decreased development time. Growth rates were higher at low temperatures, but temperature did not affect larval or adult body mass. Survivorship was highest in high-food and low-temperature treatments. Across all food and temperature treatments, northern populations exhibited a larger final body mass, shorter development time, faster growth rate, and greater survivorship than did southern populations. Results were similar for a field reciprocal-transplant experiment of third-instar larvae between populations in Connecticut and Oklahoma: Connecticut larvae grew faster than Oklahoma larvae, regardless of transplant site. Conversely, larvae transplanted to Oklahoma grew faster than larvae transplanted to Connecticut, regardless of population source. These results suggest that variation in food availability, not temperature, may account for differences in growth and body size of northern and southern ant lions. Although northern larvae grew faster and reached a larger body size in both experiments, northern environments should suppress growth because of reduced food availability and a limited growing season. This study provides the first example of countergradient selection causing Bergmann's rule in an ectotherm.

Functional trait diversity maximizes ecosystem multifunctionality

Understanding the relationship between biodiversity and ecosystem functioning has been a core ecological research topic over the last decades. Although a key hypothesis is that the diversity of functional traits determines ecosystem functioning, we do not know how much trait diversity is needed to maintain multiple ecosystem functions simultaneously (multifunctionality). Here, we uncovered a scaling relationship between the abundance distribution of two key plant functional traits (specific leaf area, maximum plant height) and multifunctionality in 124 dryland plant communities spread over all continents except Antarctica. For each trait, we found a strong empirical relationship between the skewness and the kurtosis of the trait distributions that cannot be explained by chance. This relationship predicted a strikingly high trait diversity within dryland plant communities, which was associated with a local maximization of multifunctionality. Skewness and kurtosis had a much stronger impact on multifunctionality than other important multifunctionality drivers such as species richness and aridity. The scaling relationship identified here quantifies how much trait diversity is required to maximize multifunctionality locally. Trait distributions can be used to predict the functional consequences of biodiversity loss in terrestrial ecosystems.