Unifying latitudinal gradients in range size and richness across marine and terrestrial systems

Geographic patterns of taxonomic and phylogenetic β-diversity of aquatic angiosperms in China

China covers a vast territory harbouring a large number of aquatic plants. Although there are many studies on the β-diversity of total, herbaceous or woody plants in China and elsewhere, few studies have focused on aquatic plants. Here, we analyse a comprehensive data set of 889 aquatic angiosperm species in China, and explore the geographic patterns and climatic correlates of total taxonomic and phylogenetic β-diversity as well as their turnover and nestedness components. Our results show that geographic patterns of taxonomic and phylogenetic β-diversity are highly congruent for aquatic angiosperms, and taxonomic β-diversity is consistently higher than phylogenetic β-diversity. The ratio between the nestedness component and total β-diversity is high in northwestern China and low in southeastern China. The geographic patterns of taxonomic and phylogenetic β-diversity of aquatic angiosperms in China are obviously affected by geographic and climatic distances, respectively. In conclusion, the geographic patterns of taxonomic and phylogenetic β-diversity of aquatic angiosperms are consistent across China. Climatic and geographic distances jointly affect the geographic patterns of β-diversity of aquatic angiosperms. Overall, our work provides insight into understanding the large-scale patterns of aquatic angiosperm β-diversity, and is a critical addition to previous studies on the macroecological patterns of terrestrial organisms.

Evolutionary time and species diversity in aquatic ecosystems worldwide

The latitudinal diversity gradient (LDG) is frequently described as the most dramatic biodiversity pattern on Earth, yet ecologists and biogeographers have failed to reach consensus on its primary cause. A key problem in explaining the LDG involves collinearity between multiple factors that are predicted to affect species richness in the same direction. In terrestrial systems, energy input, geographic area, and evolutionary time for species accumulation tend to covary positively with species richness at the largest spatial scales, such that their individual contributions to the LDG are confounded in global analyses. I review three diversity patterns from marine and freshwater systems that break this collinearity and which may thus provide stronger tests of the influence of time on global richness gradients. Specifically, I contrast biodiversity patterns along oceanic depth gradients, in geologically young versus ancient lakes, and in the north versus south polar marine biomes. I focus primarily on fishes due to greater data availability but synthesize patterns for invertebrates where possible. I find that regional-to-global species richness generally declines with depth in the oceans, despite the great age and stability of the deep-sea biome. Geologically ancient lakes generally do not contain more species than young lakes, and the Antarctic marine biome is not appreciably more species rich than the much younger Arctic marine biome. However, endemism is consistently higher in older systems. Patterns for invertebrate groups are less clear than for fishes and reflect a critical need for primary biodiversity data. In summary, the available data suggest that species richness is either decoupled from or only weakly related to the amount of time for diversification. These results suggest that energy, productivity, or geographic area are the primary drivers of large-scale diversity gradients. To the extent that marine and terrestrial diversity gradients result from similar processes, these examples provide evidence against a primary role for evolutionary time as the cause of the LDG.

Evolutionary drivers of the hump-shaped latitudinal gradient of benthic polychaete species richness along the Southeastern Pacific coast

Latitudinal diversity gradients (LDG) and their explanatory factors are among the most challenging topics in macroecology and biogeography. Despite of its apparent generality, a growing body of evidence shows that 'anomalous' LDG (i.e., inverse or hump-shaped trends) are common among marine organisms along the Southeastern Pacific (SEP) coast. Here, we evaluate the shape of the LDG of marine benthic polychaetes and its underlying causes using a dataset of 643 species inhabiting the continental shelf (<200 m depth), using latitudinal bands with a spatial resolution of 0.5°, along the SEP (3-56° S). The explanatory value of six oceanographic (Sea Surface Temperature (SST), SST range, salinity, salinity range, primary productivity and shelf area), and one macroecological proxy (median latitudinal range of species) were assessed using a random forest model. The taxonomic structure was used to estimate the degree of niche conservatism of predictor variables and to estimate latitudinal trends in phylogenetic diversity, based on three indices (phylogenetic richness (PD_SES), mean pairwise distance (MPD_SES), and variation of pairwise distances (VPD)). The LDG exhibits a hump-shaped trend, with a maximum peak of species richness at ca. 42° S, declining towards northern and southern areas of SEP. The latitudinal pattern was also evident in local samples controlled by sampling effort. The random forest model had a high accuracy (pseudo-r² = 0.95) and showed that the LDG could be explained by four variables (median latitudinal range, SST, salinity, and SST range), yet the functional relationship between species richness and these predictors was variable. A significant degree of phylogenetic conservatism was detected for the median latitudinal range and SST. PD_SES increased toward the southern region, whereas VPD showed the opposite trend, both statistically significant. MPD_SES has the same trend as PD_SES, but it is not significant. Our results reinforce the idea that the south Chile fjord area, particularly the Chiloé region, was likely the evolutionary source of new species of marine polychaetes along SEP, creating a hotspot of diversity. Therefore, in the same way as the canonical LDG shows a decline in diversity while moving away from the tropics; on this case the decline occurs while moving away from Chiloé Island. These results, coupled with a strong phylogenetic signal of the main predictor variables suggest that processes operating mainly at evolutionary timescales govern the LDG.

The species range-size patterns for vascular plants of Seorak Mountain (Korea): Relationship between group of life forms and phytogeography affinity along the elevational gradient

Research on species richness patterns and the advanced elevational Rapoport rule (ERR) has been widespread in recent years; however, there is a lack of such research for the temperate mountainous regions in northeast Asia. Here, we collected plant species from the Seorak Mountain in northeast Asia through field surveys. The species were divided into 11 groups according to the life-form types and phytogeography affinities of each species. The ERR was evaluated using Steven's method and by examining the species richness patterns of each group. The species richness patterns revealed a positive multimodal pattern along the elevation gradient, but phytogeography affinities (increasing trend) and life-form analysis (unimodal) exhibited different patterns. The elevation gradients (1,350 m for the mean elevation-range relationships), which are affected by the boundary effect and different life forms, did not consistently support the ERR. However, herbs as well as rare, endemic, and red list species showed consistent support for the ERR, which could be attributed to the influence by phytogeography affinities. Therefore, the results from Seorak Mountain showed that the ERR was not consistent for different plant life forms in the same area; however, phytogeography affinities could support and explain ERR.

Scale-dependent effects of niche specialisation: The disconnect between individual and species ranges

One of the most general expectations of species range dynamics is that widespread species tend to have broader niches. However, it remains unclear how this relationship is expressed at different levels of biological organisation, which involve potentially distinctive processes operating at different spatial and temporal scales. Here, we show that range sizes of terrestrial non-volant mammals at the individual and species level show contrasting relationships with two ecological niche dimensions: diet and habitat breadth. While average individual home range size appears to be mainly shaped by the interplay of diet niche breadth and body mass, species geographical range size is primarily related to habitat niche breadth but not to diet niche breadth. Our findings suggest that individual home range size is shaped by the trade-off between energetic requirements, movement capacity and trophic specialisation, whereas species geographical range size is related to the ability to persist under various environmental conditions.

Global patterns and a latitudinal gradient of flower disparity: perspectives from the angiosperm order Ericales

Morphological diversity (disparity) is an essential but often neglected aspect of biodiversity. Hence, it seems timely and promising to re-emphasize morphology in modern evolutionary studies. Disparity is a good proxy for the diversity of functions and interactions with the environment of a group of taxa. In addition, geographical and ecological patterns of disparity are crucial to understand organismal evolution and to guide biodiversity conservation efforts. Here, we analyse floral disparity across latitudinal intervals, growth forms, climate types, types of habitats, and regions for a large and representative sample of the angiosperm order Ericales. We find a latitudinal gradient of floral disparity and a decoupling of disparity from species richness. Other factors investigated are intercorrelated, and we find the highest disparity for tropical trees growing in African and South American forests. Explanations for the latitudinal gradient of floral disparity may involve the release of abiotic constraints and the increase of biotic interactions towards tropical latitudes, allowing tropical lineages to explore a broader area of the floral morphospace. Our study confirms the relevance of biodiversity parameters other than species richness and is consistent with the importance of species interactions in the tropics, in particular with respect to angiosperm flowers and their pollinators.

High Preservation Potential of Paleogeographic Range Size Distributions in Deep Time

AbstractReconstructing geographic range sizes from fossil data is a crucial tool in paleoecology, elucidating macroecological and macroevolutionary processes. Studies examining links between range size and extinction risk may also offer a predictive tool for identifying species most vulnerable in the "sixth mass extinction." However, the extent to which paleogeographic ranges can be recorded reliably in the fossil record is unknown. We perform simulation-based extinction experiments to examine (1) the fidelity of paleogeographic range size preservation in deep time, (2) the relative performance of different methods for reconstructing range size, and (3) the reliability of detecting patterns of extinction "selectivity" on range size. Our results suggest both that relative paleogeographic range size can be consistently reconstructed and that selectivity patterns on range size can be preserved under many extinction intensities, even when sedimentary rocks are scarce. By identifying patterns of selectivity across Earth's history, paleontologists can thus augment neontological work that aims to predict and prevent extinctions of living species. Last, we find that introducing "false extinctions" in the fossil record can produce spurious range-selectivity signals. Errors in the temporal ranges of species may pose a larger barrier to reconstructing range size-extinction risk signals than the spatial distribution of fossiliferous sediments.

Ecological drivers of global gradients in avian dispersal inferred from wing morphology

An organism’s ability to disperse influences many fundamental processes, from speciation and geographical range expansion to community assembly. However, the patterns and underlying drivers of variation in dispersal across species remain unclear, partly because standardised estimates of dispersal ability are rarely available. Here we present a global dataset of avian hand-wing index (HWI), an estimate of wing shape widely adopted as a proxy for dispersal ability in birds. We show that HWI is correlated with geography and ecology across 10,338 (>99%) species, increasing at higher latitudes and with migration, and decreasing with territoriality. After controlling for these effects, the strongest predictor of HWI is temperature variability (seasonality), with secondary effects of diet and habitat type. Finally, we also show that HWI is a strong predictor of geographical range size. Our analyses reveal a prominent latitudinal gradient in HWI shaped by a combination of environmental and behavioural factors, and also provide a global index of avian dispersal ability for use in community ecology, macroecology, and macroevolution.

In birds, the hand-wing index is a morphological trait that can be used as a proxy for flight efficiency. Here the authors examine variation of hand-wing index in over 10,000 bird species, finding that it is higher in migratory and non-territorial species, and lower in the tropics.

Common latitudinal gradients in functional richness and functional evenness across marine and terrestrial systems

Functional diversity is an important aspect of biodiversity, but its relationship to species diversity in time and space is poorly understood. Here we compare spatial patterns of functional and taxonomic diversity across marine and terrestrial systems to identify commonalities in their respective ecological and evolutionary drivers. We placed species-level ecological traits into comparable multi-dimensional frameworks for two model systems, marine bivalves and terrestrial birds, and used global species-occurrence data to examine the distribution of functional diversity with latitude and longitude. In both systems, tropical faunas show high total functional richness (FR) but low functional evenness (FE) (i.e. the tropics contain a highly skewed distribution of species among functional groups). Functional groups that persist toward the poles become more uniform in species richness, such that FR declines and FE rises with latitude in both systems. Temperate assemblages are more functionally even than tropical assemblages subsampled to temperate levels of species richness, suggesting that high species richness in the tropics reflects a high degree of ecological specialization within a few functional groups and/or factors that favour high recent speciation or reduced extinction rates in those groups.

Loss of Biodiversity Dimensions through Shifting Climates and Ancient Mass Extinctions

Many aspects of climate affect the deployment of biodiversity in time and space, and so changes in climate might be expected to drive regional and global extinction of both taxa and their ecological functions. Here we examine the association of past climate changes with extinction in marine bivalves, which are increasingly used as a model system for macroecological and macroevolutionary analysis. Focusing on the Cenozoic Era (66 Myr ago to the present), we analyze extinction patterns in shallow-water marine bivalve genera relative to temperature dynamics as estimated from isotopic data in microfossils. When the entire Cenozoic timeseries is considered, extinction intensity is not significantly associated with the mean temperature or the detrended variance in temperature within a given time interval (stratigraphic stage). However, extinction increases significantly with both the rate of temperature change within the stage of extinction and the absolute change in mean temperature from the preceding stage to the stage of extinction. Thus, several extinction events, particularly the extinction pulse near the Pliocene-Pleistocene boundary, do appear to have climatic drivers. Further, the latitudinal diversity gradient today and the Cenozoic history of polar faunas suggest that long-term, regional extinctions associated with cooling removed not just taxa but a variety of ecological functions from high-latitude seas. These dynamics of biodiversity loss contrast with the two mass extinctions bracketing the Mesozoic Era, which had negligible effects on the diversity of ecological functions despite removing nearly as many taxa as the latitudinal gradient does today. Thus, the fossil record raises a key issue: whether the biotic consequences of present-day stresses will more closely resemble the long-term effects of past climate changes or those that cascaded from the mass extinctions.

Extinction risk in extant marine species integrating palaeontological and biodistributional data

Extinction risk assessments of marine invertebrate species remain scarce, which hinders effective management of marine biodiversity in the face of anthropogenic impacts. To help close this information gap, in this paper we provide a metric of relative extinction risk that combines palaeontological data, in the form of extinction rates calculated from the fossil record, with two known correlates of risk in the modern day: geographical range size and realized thermal niche. We test the performance of this metric-Palaeontological Extinction Risk In Lineages (PERIL)-using survivorship analyses of Pliocene bivalve faunas from California and New Zealand, and then use it to identify present-day hotspots of extinction vulnerability for extant shallow-marine Bivalvia. Areas of the ocean where concentrations of bivalve species with higher PERIL scores overlap with high levels of climatic or anthropogenic stressors should be considered of most immediate concern for both conservation and management.

Reliable estimates of beta diversity with incomplete sampling

Beta diversity, the compositional variation among communities or assemblages, is crucial to understanding the principles of diversity assembly. The mean pairwise proportional dissimilarity expresses overall heterogeneity of samples in a data set and is among the most widely used and most robust measures of beta diversity. Obtaining a complete list of taxa and their abundances requires substantial taxonomic expertise and is time consuming. In addition, the information is generally incomplete due to sampling biases. Based on the concept of the ecological significance of dominant taxa, we explore whether determining proportional dissimilarity can be simplified based on dominant species. Using simulations and six case studies, we assess the correlation between complete community compositional data and reduced subsets of a varying number of dominant species. We find that gross beta diversity is usually depicted accurately when only the 80th percentile or five of the most abundant species of each site is considered. In data sets with very high evenness, at least the 10 most abundant species should be included. Focusing on dominant species also maintains the rank-order of beta diversity among sites. Our new approach will allow ecologists and paleobiologists to produce a far greater amount of data on diversity patterns with less time and effort, supporting conservation studies and basic science.

Geographic range velocity and its association with phylogeny and life history traits in North American woody plants

The geographic ranges of taxa change in response to environmental conditions. Yet whether rates of range movement (biotic velocities) are phylogenetically conserved is not well known. Phylogenetic conservatism of biotic velocities could reflect similarities among related lineages in climatic tolerances and dispersal-associated traits. We assess whether late Quaternary biotic velocities were phylogenetically conserved and whether they correlate with climatic tolerances and dispersal-associated traits. We used phylogenetic regression and nonparametric correlation to evaluate associations between biotic velocities, dispersal-associated traits, and climatic tolerances for 28 woody plant genera and subgenera in North America. The velocities with which woody plant taxa shifted their core geographic range limits were positively correlated from time step to time step between 16 and 7 ka. The strength of this correlation weakened after 7 ka as the pace of climate change slowed. Dispersal-associated traits and climatic tolerances were not associated with biotic velocities. Although the biotic velocities of some genera were consistently fast and others consistently slow, biotic velocities were not phylogenetically conserved. The rapid late Quaternary range shifts of plants lacking traits that facilitate frequent long-distance dispersal has long been noted (i.e., Reid's Paradox). Our results are consistent with this paradox and show that it remains robust when phylogenetic information is taken into account. The lack of association between biotic velocities, dispersal-associated traits, and climatic tolerances may reflect several, nonmutually exclusive processes, including rare long-distance dispersal, biotic interactions, and cryptic refugia. Because late Quaternary biotic velocities were decoupled from dispersal-associated traits, trait data for genera and subgenera cannot be used to predict longer-term (millennial-scale) floristic responses to climate change.

Contrasting responses of functional diversity to major losses in taxonomic diversity

Taxonomic diversity of benthic marine invertebrate shelf species declines at present by nearly an order of magnitude from the tropics to the poles in each hemisphere along the latitudinal diversity gradient (LDG), most steeply along the western Pacific where shallow-sea diversity is at its tropical maximum. In the Bivalvia, a model system for macroevolution and macroecology, this taxonomic trend is accompanied by a decline in the number of functional groups and an increase in the evenness of taxa distributed among those groups, with maximum functional evenness (FE) in polar waters of both hemispheres. In contrast, analyses of this model system across the two era-defining events of the Phanerozoic, the Permian-Triassic and Cretaceous-Paleogene mass extinctions, show only minor declines in functional richness despite high extinction intensities, resulting in a rise in FE owing to the persistence of functional groups. We hypothesize that the spatial decline of taxonomic diversity and increase in FE along the present-day LDG primarily reflect diversity-dependent factors, whereas retention of almost all functional groups through the two mass extinctions suggests the operation of diversity-independent factors. Comparative analyses of different aspects of biodiversity thus reveal strongly contrasting biological consequences of similarly severe declines in taxonomic diversity and can help predict the consequences for functional diversity among different drivers of past, present, and future biodiversity loss.

Approaches to Macroevolution: 2. Sorting of Variation, Some Overarching Issues, and General Conclusions

Approaches to macroevolution require integration of its two fundamental components, within a hierarchical framework. Following a companion paper on the origin of variation, I here discuss sorting within an evolutionary hierarchy. Species sorting-sometimes termed species selection in the broad sense, meaning differential origination and extinction owing to intrinsic biological properties-can be split into strict-sense species selection, in which rate differentials are governed by emergent, species-level traits such as geographic range size, and effect macroevolution, in which rates are governed by organism-level traits such as body size; both processes can create hitchhiking effects, indirectly causing the proliferation or decline of other traits. Several methods can operationalize the concept of emergence, so that rigorous separation of these processes is increasingly feasible. A macroevolutionary tradeoff, underlain by the intrinsic traits that influence evolutionary dynamics, causes speciation and extinction rates to covary in many clades, resulting in evolutionary volatility of some clades and more subdued behavior of others; the few clades that break the tradeoff can achieve especially prolific diversification. In addition to intrinsic biological traits at multiple levels, extrinsic events can drive the waxing and waning of clades, and the interaction of traits and events are difficult but important to disentangle. Evolutionary trends can arise in many ways, and at any hierarchical level; descriptive models can be fitted to clade trajectories in phenotypic or functional spaces, but they may not be diagnostic regarding processes, and close attention must be paid to both leading and trailing edges of apparent trends. Biotic interactions can have negative or positive effects on taxonomic diversity within a clade, but cannot be readily extrapolated from the nature of such interactions at the organismic level. The relationships among macroevolutionary currencies through time (taxonomic richness, morphologic disparity, functional variety) are crucial for understanding the nature of evolutionary diversification. A novel approach to diversity-disparity analysis shows that taxonomic diversifications can lag behind, occur in concert with, or precede, increases in disparity. Some overarching issues relating to both the origin and sorting of clades and phenotypes include the macroevolutionary role of mass extinctions, the potential differences between plant and animal macroevolution, whether macroevolutionary processes have changed through geologic time, and the growing human impact on present-day macroevolution. Many challenges remain, but progress is being made on two of the key ones: (a) the integration of variation-generating mechanisms and the multilevel sorting processes that act on that variation, and (b) the integration of paleontological and neontological approaches to historical biology.

Probabilistic models of species discovery and biodiversity comparisons

Inferring large-scale processes that drive biodiversity hinges on understanding the phylogenetic and spatial pattern of species richness. However, clades and geographic regions are accumulating newly described species at an uneven rate, potentially affecting the stability of currently observed diversity patterns. Here, we present a probabilistic model of species discovery to assess the uncertainty in diversity levels among clades and regions. We use a Bayesian time series regression to estimate the long-term trend in the rate of species description for marine bivalves and find a distinct spatial bias in the accumulation of new species. Despite these biases, probabilistic estimates of future species richness show considerable stability in the currently observed rank order of regional diversity. However, absolute differences in richness are still likely to change, potentially modifying the correlation between species numbers and geographic, environmental, and biological factors thought to promote biodiversity. Applied to scallops and related clades, we find that accumulating knowledge of deep-sea species will likely shift the relative richness of these three families, emphasizing the need to consider the incomplete nature of bivalve taxonomy in quantitative studies of its diversity. Along with estimating expected changes to observed patterns of diversity, the model described in this paper pinpoints geographic areas and clades most urgently requiring additional systematic study-an important practice for building more complete and accurate models of biodiversity dynamics that can inform ecological and evolutionary theory and improve conservation practice.

Shaping the Latitudinal Diversity Gradient: New Perspectives from a Synthesis of Paleobiology and Biogeography

An impediment to understanding the origin and dynamics of the latitudinal diversity gradient (LDG)-the most pervasive large-scale biotic pattern on Earth-has been the tendency to focus narrowly on a single causal factor when a more synthetic, integrative approach is needed. Using marine bivalves as a model system and drawing on other systems where possible, we review paleobiologic and biogeographic support for two supposedly opposing views, that the LDG is shaped primarily by (a) local environmental factors that determine the number of species and higher taxa at a given latitude (in situ hypotheses) or (b) the entry of lineages arising elsewhere into a focal region (spatial dynamics hypotheses). Support for in situ hypotheses includes the fit of present-day diversity trends in many clades to such environmental factors as temperature and the correlation of extinction intensities in Pliocene bivalve faunas with net regional temperature changes. Support for spatial dynamics hypotheses includes the age-frequency distribution of bivalve genera across latitudes, which is consistent with an out-of-the-tropics dynamic, as are the higher species diversities in temperate southeastern Australia and southeastern Japan than in the tropical Caribbean. Thus, both in situ and spatial dynamics processes must shape the bivalve LDG and are likely to operate in other groups as well. The relative strengths of the two processes may differ among groups showing similar LDGs, but dissecting their effects will require improved methods of integrating fossil data with molecular phylogenies. We highlight several potential research directions and argue that many of the most dramatic biotic patterns, past and present, are likely to have been generated by diverse, mutually reinforcing drivers.

Support for the elevational Rapoport's rule among seed plants in Nepal depends on biogeographical affinities and boundary effects

As one of the most important hypotheses on biogeographical distribution, Rapoport's rule has attracted attention around the world. However, it is unclear whether the applicability of the elevational Rapoport's Rule differs between organisms from different biogeographical regions. We used Stevens' method, which uses species diversity and the averaged range sizes of all species within each (100 m) elevational band to explore diversity-elevation, range-elevation, and diversity-range relationships. We compared support for the elevational Rapoport's rule between tropical and temperate species of seed plants in Nepal. Neither tropical nor temperate species supported the predictions of the elevational Rapoport's rule along the elevation gradient of 100-6,000 m a.s.l. for any of the studied relationships. However, along the smaller 1,000-5,000 m a.s.l. gradient (4,300 m a.s.l. for range-elevation relationships) which is thought to be less influenced by boundary effects, we observed consistent support for the rule by tropical species, although temperate species did not show consistent support. The degree of support for the elevational Rapoport's rule may not only be influenced by hard boundary effects, but also by the biogeographical affinities of the focal taxa. With ongoing global warming and increasing variability of temperature in high-elevation regions, tropical taxa may shift upward into higher elevations and expand their elevational ranges, causing the loss of temperate taxa diversity. Relevant studies on the elevational Rapoport's rule with regard to biogeographical affinities may be a promising avenue to further our understanding of this rule.