Phylogenetic Paleoecology: Tree-Thinking and Ecology in Deep Time.

The new and emerging field of phylogenetic paleoecology leverages the evolutionary relationships among species to explain temporal and spatial changes in species diversity, abundance, and distribution in deep time. This field is poised for rapid progress as knowledge of the evolutionary relationships among fossil species continues to expand. In particular, this approach will lend new insights to many of the longstanding questions in evolutionary biology, such as: the relationships among character change, ecology, and evolutionary rates; the processes that determine the evolutionary relationships among species within communities and along environmental gradients; and the phylogenetic signal underlying ecological selectivity in background and mass extinctions and in major evolutionary radiations.

Phylogenetic Clustering of Origination and Extinction across the Late Ordovician Mass Extinction.

Mass extinctions can have dramatic effects on the trajectory of life, but in some cases the effects can be relatively small even when extinction rates are high. For example, the Late Ordovician mass extinction is the second most severe in terms of the proportion of genera eliminated, yet is noted for the lack of ecological consequences and shifts in clade dominance. By comparison, the end-Cretaceous mass extinction was less severe but eliminated several major clades while some rare surviving clades diversified in the Paleogene. This disconnect may be better understood by incorporating the phylogenetic relatedness of taxa into studies of mass extinctions, as the factors driving extinction and recovery are thought to be phylogenetically conserved and should therefore promote both origination and extinction of closely related taxa. Here, we test whether there was phylogenetic selectivity in extinction and origination using brachiopod genera from the Middle Ordovician through the Devonian. Using an index of taxonomic clustering (RCL) as a proxy for phylogenetic clustering, we find that A) both extinctions and originations shift from taxonomically random or weakly clustered within families in the Ordovician to strongly clustered in the Silurian and Devonian, beginning with the recovery following the Late Ordovician mass extinction, and B) the Late Ordovician mass extinction was itself only weakly clustered. Both results stand in stark contrast to Cretaceous-Cenozoic bivalves, which showed significant levels of taxonomic clustering of extinctions in the Cretaceous, including strong clustering in the mass extinction, but taxonomically random extinctions in the Cenozoic. The contrasting patterns between the Late Ordovician and end-Cretaceous events suggest a complex relationship between the phylogenetic selectivity of mass extinctions and the long-term phylogenetic signal in origination and extinction patterns.

Origination and immigration drive latitudinal gradients in marine functional diversity.

Global patterns in the functional attributes of organisms are critical to understanding biodiversity trends and predicting biotic responses to environmental change. In the first global marine analysis, we find a strong decrease in functional richness, but a strong increase in functional evenness, with increasing latitude using intertidal-to-outer-shelf bivalves as a model system (N = 5571 species). These patterns appear to be driven by the interplay between variation in origination rates among functional groups, and latitudinal patterns in origination and range expansion, as documented by the rich fossil record of the group. The data suggest that (i) accumulation of taxa in spatial bins and functional categories has not impeded continued diversification in the tropics, and (ii) extinctions will influence ecosystem function differentially across latitudes.

Out of the tropics, but how? Fossils, bridge species, and thermal ranges in the dynamics of the marine latitudinal diversity gradient.

Latitudinal diversity gradients are underlain by complex combinations of origination, extinction, and shifts in geographic distribution and therefore are best analyzed by integrating paleontological and neontological data. The fossil record of marine bivalves shows, in three successive late Cenozoic time slices, that most clades (operationally here, genera) tend to originate in the tropics and then expand out of the tropics (OTT) to higher latitudes while retaining their tropical presence. This OTT pattern is robust both to assumptions on the preservation potential of taxa and to taxonomic revisions of extant and fossil species. Range expansion of clades may occur via "bridge species," which violate climate-niche conservatism to bridge the tropical-temperate boundary in most OTT genera. Substantial time lags (∼5 Myr) between the origins of tropical clades and their entry into the temperate zone suggest that OTT events are rare on a per-clade basis. Clades with higher diversification rates within the tropics are the most likely to expand OTT and the most likely to produce multiple bridge species, suggesting that high speciation rates promote the OTT dynamic. Although expansion of thermal tolerances is key to the OTT dynamic, most latitudinally widespread species instead achieve their broad ranges by tracking widespread, spatially-uniform temperatures within the tropics (yielding, via the nonlinear relation between temperature and latitude, a pattern opposite to Rapoport's rule). This decoupling of range size and temperature tolerance may also explain the differing roles of species and clade ranges in buffering species from background and mass extinctions.

Global environmental predictors of benthic marine biogeographic structure.

Analyses of how environmental factors influence the biogeographic structure of biotas are essential for understanding the processes underlying global diversity patterns and for predicting large-scale biotic responses to global change. Here we show that the large-scale geographic structure of shallow-marine benthic faunas, defined by existing biogeographic schemes, can be predicted with 89-100% accuracy by a few readily available oceanographic variables; temperature alone can predict 53-99% of the present-day structure along coastlines. The same set of variables is also strongly correlated with spatial changes in species compositions of bivalves, a major component of the benthic marine biota, at the 1° grid-cell resolution. These analyses demonstrate the central role of coastal oceanography in structuring benthic marine biogeography and suggest that a few environmental variables may be sufficient to model the response of marine biogeographic structure to past and future changes in climate.

Genus age, provincial area and the taxonomic structure of marine faunas.

Species are unevenly distributed among genera within clades and regions, with most genera species-poor and few species-rich. At regional scales, this structure to taxonomic diversity is generated via speciation, extinction and geographical range dynamics. Here, we use a global database of extant marine bivalves to characterize the taxonomic structure of climate zones and provinces. Our analyses reveal a general, Zipf-Mandelbrot form to the distribution of species among genera, with faunas from similar climate zones exhibiting similar taxonomic structure. Provinces that contain older taxa and/or encompass larger areas are expected to be more species-rich. Although both median genus age and provincial area correlate with measures of taxonomic structure, these relationships are interdependent, nonlinear and driven primarily by contrasts between tropical and extra-tropical faunas. Provincial area and taxonomic structure are largely decoupled within climate zones. Counter to the expectation that genus age and species richness should positively covary, diverse and highly structured provincial faunas are dominated by young genera. The marked differences between tropical and temperate faunas suggest strong spatial variation in evolutionary rates and invasion frequencies. Such variation contradicts biogeographic models that scale taxonomic diversity to geographical area.

Differential extinction and the contrasting structure of polar marine faunas.

BACKGROUND: The low taxonomic diversity of polar marine faunas today reflects both the failure of clades to colonize or diversify in high latitudes and regional extinctions of once-present clades. However, simple models of polar evolution are made difficult by the strikingly different faunal compositions and community structures of the two poles. METHODOLOGY/PRINCIPAL FINDINGS: A comparison of early Cenozoic Arctic and Antarctic bivalve faunas with modern ones, within the framework of a molecular phylogeny, shows that while Arctic losses were randomly distributed across the tree, Antarctic losses were significantly concentrated in more derived families, resulting in communities dominated by basal lineages. Potential mechanisms for the phylogenetic structure to Antarctic extinctions include continental isolation, changes in primary productivity leading to turnover of both predators and prey, and the effect of glaciation on shelf habitats. CONCLUSIONS/SIGNIFICANCE: These results show that phylogenetic consequences of past extinctions can vary substantially among regions and thus shape regional faunal structures, even when due to similar drivers, here global cooling, and provide the first phylogenetic support for the "retrograde" hypothesis of Antarctic faunal evolution.

A macroevolutionary perspective on species range limits.

Understanding the factors that determine the geographic range limits of species is important for many questions in ecology, evolution and conservation biology. These limits arise from complex interactions among ecology and dispersal ability of species and the physical environment, but many of the underlying traits can be conserved among related species and clades. Thus, the range limits of species are likely to be influenced by their macroevolutionary history. Using palaeontological and biogeographic data for marine bivalves, we find that the range limits of genera are significantly related to their constituent species richness, but the effects of age are weak and inconsistent. In addition, we find a significant phylogenetic signal in the range limits at both genus and family levels, although the strength of this effect shows interoceanic variation. This phylogenetic conservatism of range limits gives rise to an evolutionary pattern where wide-ranging lineages have clusters of species within the biogeographic provinces, with a few extending across major boundaries.

Signature of the end-Cretaceous mass extinction in the modern biota.

The long-term effects of mass extinctions on spatial and evolutionary dynamics have been poorly studied. Here we show that the evolutionary consequences of the end-Cretaceous [Cretaceous/Paleogene (K/Pg)] mass extinction persist in present-day biogeography. The geologic ages of genera of living marine bivalves show a significant break from a smooth exponential distribution, corresponding to the K/Pg boundary. The break reflects a permanent increase in origination rates, intermediate between the Mesozoic rate and the post-extinction recovery pulse. This global rate shift is most clearly seen today in tropical bioprovinces and weakens toward the poles. Coupled with the modern geographic distributions of taxa originating before and after the K/Pg boundary, this spatial pattern indicates that tropical origination rates after the K/Pg event have left a permanent mark on the taxonomic and biogeographic structure of the modern biota, despite the complex Cenozoic history of marine environments.

Generation of Earth's first-order biodiversity pattern.

The first-order biodiversity pattern on Earth today and at least as far back as the Paleozoic is the latitudinal diversity gradient (LDG), a decrease in richness of species and higher taxa from the equator to the poles. LDGs are produced by geographic trends in origination, extinction, and dispersal over evolutionary timescales, so that analyses of static patterns will be insufficient to reveal underlying processes. The fossil record of marine bivalve genera, a model system for the analysis of biodiversity dynamics over large temporal and spatial scales, shows that an origination and range-expansion gradient plays a major role in generating the LDG. Peak origination rates and peak diversities fall within the tropics, with range expansion out of the tropics the predominant spatial dynamic thereafter. The origination-diversity link occurs even in a "contrarian" group whose diversity peaks at midlatitudes, an exception proving the rule that spatial variations in origination are key to latitudinal diversity patterns. Extinction rates are lower in polar latitudes (> or =60 degrees ) than in temperate zones and thus cannot create the observed gradient alone. They may, however, help to explain why origination and immigration are evidently damped in higher latitudes. We suggest that species require more resources in higher latitudes, for the seasonality of primary productivity increases by more than an order of magnitude from equatorial to polar regions. Higher-latitude species are generalists that, unlike potential immigrants, are adapted to garner the large share of resources required for incumbency in those regions. When resources are opened up by extinctions, lineages spread chiefly poleward and chiefly through speciation.

Phanerozoic trends in the global diversity of marine invertebrates.

It has previously been thought that there was a steep Cretaceous and Cenozoic radiation of marine invertebrates. This pattern can be replicated with a new data set of fossil occurrences representing 3.5 million specimens, but only when older analytical protocols are used. Moreover, analyses that employ sampling standardization and more robust counting methods show a modest rise in diversity with no clear trend after the mid-Cretaceous. Globally, locally, and at both high and low latitudes, diversity was less than twice as high in the Neogene as in the mid-Paleozoic. The ratio of global to local richness has changed little, and a latitudinal diversity gradient was present in the early Paleozoic.

Species-genus ratios reflect a global history of diversification and range expansion in marine bivalves.

The distribution of marine bivalve species among genera and higher taxa takes the form of the classic hollow curve, wherein few lineages are species rich and many are species poor. The distribution of species among genera (S/G ratio) varies with latitude, with temperate S/G's falling within the null expectation, and tropical and polar S/G's exceeding it. Here, we test several hypotheses for this polar overdominance in the species richness of small numbers of genera. We find a significant positive correlation between the latitudinal range of a genus and its species richness, both globally and within regions. Genus age and species richness are also positively related, but this relationship breaks down when the analysis is limited to genera endemic to climate zones or with narrow latitudinal ranges. The data suggest a link between speciation and range-expansion, with genera expanding out of the tropical latitudinal bins tending to speciate more prolifically, both globally and regionally. These genera contain more species within climate zones than taxa endemic to that zone. Range expansion thus appears to be fundamentally coupled with speciation, producing the skewed distribution of species among genera, both globally and regionally, whereas clade longevity is achieved through extinction -- resistance conferred by broad geographical ranges.

Contrarian clade confirms the ubiquity of spatial origination patterns in the production of latitudinal diversity gradients.

The latitudinal diversity gradient (LDG), wherein the number of species and higher taxa peaks in the tropics and decreases toward the poles, is the best-documented large-scale diversity pattern on Earth, but hypotheses explaining the standard LDG must also account for rare "contrarian" taxa that show diversity maxima outside of the tropics. For marine bivalves, one of the few groups that provide spatially explicit temporal data on a global scale, we show that a major contrarian group, the Anomalodesmata, unexpectedly exhibits the same large-scale dynamics as related clades having normal LDGs in two key respects. First, maxima in standing genus diversity and genus origination rates coincide spatially. Second, the strength of a clade's present-day LDG is significantly related to the proportion of its living genera that originated in the tropics during the late Cenozoic, with the contrarian gradient strength at both species and genus level predicted quantitatively by the values for the other clades. Geologic age distributions indicate that the anomalous LDG results from origination that is damped in the tropics rather than heightened in the temperate zones. The pervasive role of spatial origination patterns in shaping LDGs, regardless of the position of their diversity maxima, corroborates hypotheses based on clades showing standard gradients and underscores the insights that contrarian groups can provide into general principles of diversity dynamics.