Warm and arid regions of the world are hotspots of superorganism complexity.

Biologists have long been fascinated by the processes that give rise to phenotypic complexity of organisms, yet whether there exist geographical hotspots of phenotypic complexity remains poorly explored. Phenotypic complexity can be readily observed in ant colonies, which are superorganisms with morphologically differentiated queen and worker castes analogous to the germline and soma of multicellular organisms. Several ant species have evolved 'worker polymorphism', where workers in a single colony show quantifiable differences in size and head-to-body scaling. Here, we use 256 754 occurrence points from 8990 ant species to investigate the geography of worker polymorphism. We show that arid regions of the world are the hotspots of superorganism complexity. Tropical savannahs and deserts, which are typically species-poor relative to tropical or even temperate forests, harbour the highest densities of polymorphic ants. We discuss the possible adaptive advantages that worker polymorphism provides in arid environments. Our work may provide a window into the environmental conditions that promote the emergence of highly complex phenotypes.

Sodium co-limits and catalyzes macronutrients in a prairie food web.

Nitrogen and phosphorus frequently limit terrestrial plant production, but have a mixed record in regulating the abundance of terrestrial invertebrates. We contrasted four ways that Na could interact with an NP fertilizer to shape the plants and invertebrates of an inland prairie. We applied NP and Na to m2 plots in a factorial design. Aboveground invertebrate abundance was independently co-limited by NaCl and NP, but with +NP plots supporting more individuals. We suggest the disparity arises because NP enhanced plant height by 35% (1 SD) over controls, providing both food and habitat, whereas NaCl provides only food. Belowground invertebrates showed evidence of serial co-limitation, where NaCl additions alone were ineffectual, but catalyzed access to NP. This suggests the increased belowground food availability in NP plots increased Na demand. Na and NP supply rates vary with climate, land use, and with inputs like urine. The co-limitation and catalysis of N and P by Na thus has the potential for predicting patterns of abundance and diversity across spatial scales.

Biogeochemistry drives diversity in the prokaryotes, fungi, and invertebrates of a Panama forest.

Humans are both fertilizing the world and depleting its soils, decreasing the diversity of aquatic ecosystems and terrestrial plants in the process. We know less about how nutrients shape the abundance and diversity of the prokaryotes, fungi, and invertebrates of Earth's soils. Here we explore this question in the soils of a Panama forest subject to a 13-yr fertilization with factorial combinations of nitrogen (N), phosphorus (P), and potassium (K) and a separate micronutrient cocktail. We contrast three hypotheses linking biogeochemistry to abundance and diversity. Consistent with the Stress Hypothesis, adding N suppressed the abundance of invertebrates and the richness of all three groups of organisms by ca. 1 SD or more below controls. Nitrogen addition plots were 0.8 pH units more acidic with 18% more exchangeable aluminum, which is toxic to both prokaryotes and eukaryotes. These stress effects were frequently reversed, however, when N was added with P (for prokaryotes and invertebrates) and with added K (for fungi). Consistent with the Abundance Hypothesis, adding P generally increased prokaryote and invertebrate diversity, and adding K enhanced invertebrate diversity. Also consistent with the Abundance Hypothesis, increases in invertebrate abundance generated increases in richness. We found little evidence for the Competition Hypothesis: that single nutrients suppressed diversity by favoring a subset of high nutrient specialists, and that nutrient combinations suppressed diversity even more. Instead, combinations of nutrients, and especially the cation/micronutrient treatment, yielded the largest increases in richness in the two eukaryote groups. In sum, changes in soil biogeochemistry revealed a diversity of responses among the three dominant soil groups, positive synergies among nutrients, and-in contrast with terrestrial plants-the frequent enhancement of soil biodiversity.

The latitudinal species richness gradient in New World woody angiosperms is consistent with the tropical conservatism hypothesis.

Plant diversity, like that of most other taxonomic groups, peaks in the tropics, where climatic conditions are warm and wet, and it declines toward the temperate and polar zones as conditions become colder and drier, with more seasonally variable temperatures. Climate and evolutionary history are often considered competing explanations for the latitudinal gradient, but they are linked by the evolutionarily conserved environmental adaptations of species and the history of Earth's climate system. The tropical conservatism hypothesis (TCH) invokes niche conservatism, climatic limitations on establishment and survival, and paleoclimatic history to explain the latitudinal diversity gradient. Here, we use latitudinal distributions for over 12,500 woody angiosperm species, a fossil-calibrated supertree, and null modeling to test predictions of the TCH. Regional assemblages in the northern and southern temperate zones are less phylogenetically diverse than expected based on their species richness, because temperate taxa are clustered into relatively few clades. Moreover, lineages with temperate affinities are generally younger and nested within older, more tropical lineages. As predicted by the TCH, the vast majority of temperate lineages have arisen since global cooling began at the Eocene-Oligocene boundary (34 Mya). By linking physiological tolerances of species to evolutionary and biogeographic processes, phylogenetic niche conservatism may provide a theoretical framework for a generalized explanation for Earth's predominant pattern of biodiversity.

Constancy in Functional Space across a Species Richness Anomaly.

The relationship between large-scale gradients in species richness and functional diversity provides important information regarding the mechanisms driving patterns of biodiversity. A classic hypothesis in ecology is that strong interspecific interactions should result in an increase in the functional volume of assemblages as the species richness increases, whereas climatic constraints may result in no change in functional volume. Most research of this kind examines latitudinal gradients in species richness, but the results are likely confounded by underlying gradients in climate and phylogenetic composition. We take an alternative approach that examines functional richness across a tree species richness anomaly where species richness doubles from Europe to eastern North America. The results demonstrate that the functional richness on both continents saturates at a similar point as species richness increases and that the packing of functional space becomes tighter. Further, the species richness anomaly is driven primarily by genera unique to North America, but those genera contribute less than expected functional richness to the region, indicating a high level of redundancy with genera shared between the continents. Taken together, the results indicate that the species richness anomaly is associated with diversification within a climatically constrained trait space. More generally, the work demonstrates the power of utilizing species richness anomalies in biodiversity research, particularly when they are coupled with information regarding organismal function.

Biogeographic patterns of soil diazotrophic communities across six forests in North America.

Soil diazotrophs play important roles in ecosystem functioning by converting atmospheric N2 into biologically available ammonium. However, the diversity and distribution of soil diazotrophic communities in different forests and whether they follow biogeographic patterns similar to macroorganisms still remain unclear. By sequencing nifH gene amplicons, we surveyed the diversity, structure and biogeographic patterns of soil diazotrophic communities across six North American forests (126 nested samples). Our results showed that each forest harboured markedly different soil diazotrophic communities and that these communities followed traditional biogeographic patterns similar to plant and animal communities, including the taxa-area relationship (TAR) and latitudinal diversity gradient. Significantly higher community diversity and lower microbial spatial turnover rates (i.e. z-values) were found for rainforests (~0.06) than temperate forests (~0.1). The gradient pattern of TARs and community diversity was strongly correlated (r(2)  > 0.5) with latitude, annual mean temperature, plant species richness and precipitation, and weakly correlated (r(2)  < 0.25) with pH and soil moisture. This study suggests that even microbial subcommunities (e.g. soil diazotrophs) follow general biogeographic patterns (e.g. TAR, latitudinal diversity gradient), and indicates that the metabolic theory of ecology and habitat heterogeneity may be the major underlying ecological mechanisms shaping the biogeographic patterns of soil diazotrophic communities.

Global phylogenetic structure of the hyperdiverse ant genus Pheidole reveals the repeated evolution of macroecological patterns.

Adaptive radiations are of particular interest owing to what they reveal about the ecological and evolutionary regulation of biodiversity. This applies to localized island radiations such as Darwin's finches, and also to rapid radiations occurring on a global scale. Here we analyse the macroevolution and macroecology of Pheidole, a famously hyperdiverse and ecologically dominant ant genus. We generate and analyse four novel datasets: (i) a robust global phylogeny including 285 Pheidole species, (ii) a global database on regional Pheidole richness in 365 political areas summarizing over 97 000 individual records from more than 6500 studies, (iii) a global database of Pheidole richness from 3796 local communities and (iv) a database of Pheidole body sizes across species. Analysis of the potential climate drivers of richness revealed that the patterns are statistically very similar across different biogeographic regions, with both regional and local richness associated with the same coefficients of temperature and precipitation. This similarity occurs even though phylogenetic analysis shows that Pheidole reached dominance in communities through serial localized radiations into different biomes within different continents and islands. Pheidole body size distributions have likewise converged across geographical regions. We propose these cases of convergence indicate that the global radiation of Pheidole is structured by deterministic factors regulating diversification and diversity.

Global models of ant diversity suggest regions where new discoveries are most likely are under disproportionate deforestation threat.

Most of the described and probably undescribed species on Earth are insects. Global models of species diversity rarely focus on insects and none attempt to address unknown, undescribed diversity. We assembled a database representing about 13,000 records for ant generic distribution from over 350 regions that cover much of the globe. Based on two models of diversity and endemicity, we identified regions where our knowledge of ant diversity is most limited, regions we have called "hotspots of discovery." A priori, such regions might be expected to be remote and untouched. Instead, we found that the hotspots of discovery are also the regions in which biodiversity is the most threatened by habitat destruction. Our results not only highlight the immediate need for conservation of the remaining natural habitats in these regions, but also the extent to which, by focusing on well-known groups such as vertebrates, we may fail to conserve the far greater diversity of the smaller species yet to be found.

Capacity and establishment rules govern the number of nonnative species in communities of ground-dwelling invertebrates.

Nonnative species are a key agent of global change. However, nonnative invertebrates remain understudied at the community scales where they are most likely to drive local extirpations. We use the North American NEON pitfall trapping network to document the number of nonnative species from 51 invertebrate communities, testing four classes of drivers. We sequenced samples using the eDNA from the sample's storage ethanol. We used AICc informed regression to evaluate how native species richness, productivity, habitat, temperature, and human population density and vehicular traffic account for continent-wide variation in the number of nonnative species in a local community. The percentage of nonnatives varied 3-fold among habitat types and over 10-fold (0%-14%) overall. We found evidence for two types of constraints on nonnative diversity. Consistent with Capacity rules (i.e., how the number of niches and individuals reflect the number of species an ecosystem can support) nonnatives increased with existing native species richness and ecosystem productivity. Consistent with Establishment Rules (i.e., how the dispersal rate of nonnative propagules and the number of open sites limits nonnative species richness) nonnatives increased with automobile traffic-a measure of human-generated propagule pressure-and were twice as common in pastures than native grasslands. After accounting for drivers associated with a community's ability to support native species (native species richness and productivity), nonnatives are more common in communities that are regularly seasonally disturbed (pastures and, potentially deciduous forests) and those experiencing more vehicular traffic. These baseline values across the US North America will allow NEON's monitoring mission to document how anthropogenic change-from disturbance to propagule transport, from temperature to trends in local extinction-further shape biotic homogenization.

Clinal variation in colony breeding structure and level of inbreeding in the subterranean termites Reticulitermes flavipes and R. grassei.

Social insects exhibit remarkable variation in their colony breeding structures, both within and among species. Ecological factors are believed to be important in shaping reproductive traits of social insect colonies, yet there is little information linking specific environmental variables with differences in breeding structure. Subterranean termites (Rhinotermitidae) show exceptional variation in colony breeding structure, differing in the number of reproductives and degree of inbreeding; colonies can be simple families headed by a single pair of monogamous reproductives (king and queen) or they can be extended families headed by multiple inbreeding neotenic reproductives (wingless individuals). Using microsatellite markers, we characterized colony breeding structure and levels of inbreeding in populations over large parts of the range of the subterranean termites Reticulitermes flavipes in the USA and R. grassei in Europe. Combining these new data with previous results on populations of both species, we found that latitude had a strong effect on the proportion of extended-family colonies in R. flavipes and on levels of inbreeding in both species. We examined the effect of several environmental variables that vary latitudinally; while the degree of inbreeding was greatest in cool, moist habitats in both species, seasonality affected the species differently. Inbreeding in R. flavipes was most strongly associated with climatic variables (mean annual temperature and seasonality), whereas nonclimatic variables, including the availability of wood substrate and soil composition, were important predictors of inbreeding in R. grassei. These results are the first showing that termite breeding structure is shaped by local environmental factors and that species can vary in their responses to these factors.

Temperature-habitat interactions constrain seasonal activity in a continental array of pitfall traps.

Activity density (AD), the rate at which animals collectively move through their environment, emerges as the product of a taxon's local abundance and its velocity. We analyze drivers of seasonal AD using 47 localities from the National Ecological Observatory Network (NEON) both to better understand variation in ecosystem rates like pollination and seed dispersal as well as the constraints of using AD to monitor invertebrate populations. AD was measured as volume from biweekly pitfall trap arrays (ml trap-1 14 days-1 ). Pooled samples from 2017 to 2018 revealed AD extrema at most temperatures but with a strongly positive overall slope. However, habitat types varied widely in AD's seasonal temperature sensitivity, from negative in wetlands to positive in mixed forest, grassland, and shrub habitats. The temperature of maximum AD varied threefold across the 47 localities; it tracked the threefold geographic variation in maximum growing season temperature with a consistent gap of ca. 3°C across habitats, a novel macroecological result. AD holds potential as an effective proxy for investigating ecosystem rates driven by activity. However, our results suggest that its use for monitoring insect abundance is complicated by the many ways that both abundance and velocity are constrained by a locality's temperature and plant physiognomy.

Strong influence of regional species pools on continent-wide structuring of local communities.

There is a long tradition in ecology of evaluating the relative contribution of the regional species pool and local interactions on the structure of local communities. Similarly, a growing number of studies assess the phylogenetic structure of communities, relative to that in the regional species pool, to examine the interplay between broad-scale evolutionary and fine-scale ecological processes. Finally, a renewed interest in the influence of species source pools on communities has shown that the definition of the source pool influences interpretations of patterns of community structure. We use a continent-wide dataset of local ant communities and implement ecologically explicit source pool definitions to examine the relative importance of regional species pools and local interactions for shaping community structure. Then we assess which factors underlie systematic variation in the structure of communities along climatic gradients. We find that the average phylogenetic relatedness of species in ant communities decreases from tropical to temperate regions, but the strength of this relationship depends on the level of ecological realism in the definition of source pools. We conclude that the evolution of climatic niches influences the phylogenetic structure of regional source pools and that the influence of regional source pools on local community structure is strong.

Testing the role of body size and litter depth on invertebrate diversity across six forests in North America.

Ecologists search for rules by which traits dictate the abundance and distribution of species. Here we search for rules that apply across three common taxa of litter invertebrates in six North American forests from Panama to Oregon. We use image analysis to quantify the abundance and body size distributions of mites, springtails, and spiders in 21 1-m2 plots per forest. We contrast three hypotheses: two of which focus on trait-abundance relationships and a third linking abundance to species richness. Despite three orders of magnitude variation in size, the predicted negative relationship between mean body size and abundance per area occurred in only 18% of cases, never for large bodied taxa like spiders. We likewise found only 18% of tests supported our prediction that increasing litter depth allows for high abundance; two-thirds of which occurred at a single deciduous forest in Massachusetts. In contrast, invertebrate abundance constrained species richness 76% of the time. Our results suggest that body size and habitat volume in brown food webs are rarely good predictors of variation in abundance, but that variation in diversity is generally well predicted by abundance.

Activity density at a continental scale: What drives invertebrate biomass moving across the soil surface?

Activity density (AD), the rate that an individual taxon or its biomass moves through the environment, is used both to monitor communities and quantify the potential for ecosystem work. The Abundance Velocity Hypothesis posited that AD increases with aboveground net primary productivity (ANPP) and is a unimodal function of temperature. Here we show that, at continental extents, increasing ANPP may have nonlinear effects on AD: increasing abundance, but decreasing velocity as accumulating vegetation interferes with movement. We use 5 yr of data from the NEON invertebrate pitfall trap arrays including 43 locations and four habitat types for a total of 77 habitat-site combinations to evaluate continental drivers of invertebrate AD. ANPP and temperature accounted for one-third to 92% of variation in AD. As predicted, AD was a unimodal function of temperature in forests and grasslands but increased linearly in open scrublands. ANPP yielded further nonlinear effects, generating unimodal AD curves in wetlands, and bimodal curves in forests. While all four habitats showed no AD trends over 5 yr of sampling, these nonlinearities suggest that trends in AD, often used to infer changes in insect abundance, will vary qualitatively across ecoregions.

Plant geography upon the basis of functional traits: an example from eastern North American trees.

Plant geographers have sought for decades to describe and predict the geographic distribution of vegetation types on the basis of plant function and its relationship with the abiotic environment. Traditionally this has been accomplished using categorical representations such as plant functional types. Increasingly, plant functional ecologists have sought to refine categorical functional types via quantitative functional traits in order to understand the ecological implications of trade-offs in plant form and function. Fewer works have focused upon testing whether commonly measured functional traits enhance our understanding of plant biogeography broadly and the geographic distribution of vegetation types in particular. Here we combine a continental-scale forest inventory data set containing 18 111 plots with a plant functional trait data set to ask: (1) Is there a strong relationship between the abiotic environment and the distribution of functional trait values in forest inventory plots? And (2) can different Holdridge life zones be distinguished upon the basis of their functional trait distributions? The results show geographic patterns of functional trait distributions that are often strongly correlated with climate and also show that the Holdridge life zones in the study area can be differentiated using a combination of functional traits.

More individuals but fewer species: testing the 'more individuals hypothesis' in a diverse tropical fauna.

A positive relationship between species richness and productivity is often observed in nature, but the causes remain contentious. One mechanism, the 'more individuals hypothesis' (MIH), predicts richness increases monotonically with density, as a function of resource flux. To test the MIH, we manipulated resource abundance in a community of tropical rainforest litter ants and measured richness and density responses. A unimodal relationship between richness and density most closely fitted the control and disturbance (resource removal) treatments in contrast to expectations of the MIH. Resource addition resulted in a monotonic increase in richness relative to density, a shift from the pattern in the control. In the disturbance treatment, richness was greater than in the control, opposite to expectations of the MIH. While large-scale correlations between ant diversity and net primary productivity or temperature are reconcilable with the MIH, key elements of the hypothesis are not supported.

Canopy and litter ant assemblages share similar climate-species density relationships.

Tropical forest canopies house most of the globe's diversity, yet little is known about global patterns and drivers of canopy diversity. Here, we present models of ant species density, using climate, abundance and habitat (i.e. canopy versus litter) as predictors. Ant species density is positively associated with temperature and precipitation, and negatively (or non-significantly) associated with two metrics of seasonality, precipitation seasonality and temperature range. Ant species density was significantly higher in canopy samples, but this difference disappeared once abundance was considered. Thus, apparent differences in species density between canopy and litter samples are probably owing to differences in abundance-diversity relationships, and not differences in climate-diversity relationships. Thus, it appears that canopy and litter ant assemblages share a common abundance-diversity relationship influenced by similar but not identical climatic drivers.

Robust and simplified machine learning identification of pitfall trap-collected ground beetles at the continental scale.

Insect populations are changing rapidly, and monitoring these changes is essential for understanding the causes and consequences of such shifts. However, large-scale insect identification projects are time-consuming and expensive when done solely by human identifiers. Machine learning offers a possible solution to help collect insect data quickly and efficiently.Here, we outline a methodology for training classification models to identify pitfall trap-collected insects from image data and then apply the method to identify ground beetles (Carabidae). All beetles were collected by the National Ecological Observatory Network (NEON), a continental scale ecological monitoring project with sites across the United States. We describe the procedures for image collection, image data extraction, data preparation, and model training, and compare the performance of five machine learning algorithms and two classification methods (hierarchical vs. single-level) identifying ground beetles from the species to subfamily level. All models were trained using pre-extracted feature vectors, not raw image data. Our methodology allows for data to be extracted from multiple individuals within the same image thus enhancing time efficiency, utilizes relatively simple models that allow for direct assessment of model performance, and can be performed on relatively small datasets.The best performing algorithm, linear discriminant analysis (LDA), reached an accuracy of 84.6% at the species level when naively identifying species, which was further increased to >95% when classifications were limited by known local species pools. Model performance was negatively correlated with taxonomic specificity, with the LDA model reaching an accuracy of ~99% at the subfamily level. When classifying carabid species not included in the training dataset at higher taxonomic levels species, the models performed significantly better than if classifications were made randomly. We also observed greater performance when classifications were made using the hierarchical classification method compared to the single-level classification method at higher taxonomic levels.The general methodology outlined here serves as a proof-of-concept for classifying pitfall trap-collected organisms using machine learning algorithms, and the image data extraction methodology may be used for nonmachine learning uses. We propose that integration of machine learning in large-scale identification pipelines will increase efficiency and lead to a greater flow of insect macroecological data, with the potential to be expanded for use with other noninsect taxa.

Climatic drivers of hemispheric asymmetry in global patterns of ant species richness.

Although many taxa show a latitudinal gradient in richness, the relationship between latitude and species richness is often asymmetrical between the northern and southern hemispheres. Here we examine the latitudinal pattern of species richness across 1003 local ant assemblages. We find latitudinal asymmetry, with southern hemisphere sites being more diverse than northern hemisphere sites. Most of this asymmetry could be explained statistically by differences in contemporary climate. Local ant species richness was positively associated with temperature, but negatively (although weakly) associated with temperature range and precipitation. After contemporary climate was accounted for, a modest difference in diversity between hemispheres persisted, suggesting that factors other than contemporary climate contributed to the hemispherical asymmetry. The most parsimonious explanation for this remaining asymmetry is that greater climate change since the Eocene in the northern than in the southern hemisphere has led to more extinctions in the northern hemisphere with consequent effects on local ant species richness.

Species energy and Thermal Performance Theory predict 20-yr changes in ant community abundance and richness.

In an era of rapid climate change, and with it concern over insect declines, we used two theories to predict 20-yr changes in 34 North American ant communities. The ecosystems, from deserts to hardwood forests, were first surveyed in the 1990s. When resurveyed in 2016-2017, they averaged 1°C warmer with 200 g C·m-2 ·yr-1 higher plant productivity. Ant colony abundance changed from -49% to +61%. Consistent with Thermal Performance Theory, colony abundance increased with temperature increases < 1°C, then decreased as a site's mean monthly temperature change increased up to +2.4°C. Consistent with Species Energy Theory, (1) ant abundance tracked changes in a measure of energy availability (net aboveground productivity, g C·m-2 ·yr-1 ) and (2) increases in colony abundance drove increases in local plot- and transect-level species richness but not that of Chao 2, an estimate of the size of the species pool. Even after accounting for these drivers, local species richness was still higher ~20 yr after the original surveys, likely due to the increased activity of ant workers. These results suggest community changes are predictable using theory from geographical ecology, and that warming can first enhance but may ultimately decrease the abundance of this important insect taxon.