Comparing the rates of speciation and extinction between phylogenetic trees

phytools 2.0: an updated R ecosystem for phylogenetic comparative methods (and other things)

Phylogenetic comparative methods comprise the general endeavor of using an estimated phylogenetic tree (or set of trees) to make secondary inferences: about trait evolution, diversification dynamics, biogeography, community ecology, and a wide range of other phenomena or processes. Over the past ten years or so, the phytools R package has grown to become an important research tool for phylogenetic comparative analysis. phytools is a diverse contributed R library now consisting of hundreds of different functions covering a variety of methods and purposes in phylogenetic biology. As of the time of writing, phytools included functionality for fitting models of trait evolution, for reconstructing ancestral states, for studying diversification on trees, and for visualizing phylogenies, comparative data, and fitted models, as well numerous other tasks related to phylogenetic biology. Here, I describe some significant features of and recent updates to phytools, while also illustrating several popular workflows of the phytools computational software.

Macroevolutionary trends and diversification dynamics in Atripliceae (Amaranthaceae s.l., Chenopodioideae): a first approach

BACKGROUND AND AIMS

Atripliceae evolved and diversified by dispersals and radiations across continents in both hemispheres, colonizing similar semi-arid, saline-alkaline environments throughout the world. Meanwhile, its species developed different life forms, photosynthetic pathways, mono- or dioecy, and different morphological features in flowers, fruiting bracteoles and seeds. In this study, we introduce a first approach to the macroevolutionary patterns and diversification dynamics of the Atripliceae to understand how time, traits, speciation, extinction and new habitats influenced the evolution of this lineage.

METHODS

We performed molecular phylogenetic analyses and clade age estimation of Atripliceae to apply time-, trait- and geographic-dependent diversification analyses and ancestral state reconstructions to explore diversification patterns within the tribe.

KEY RESULTS

Opposite diversification dynamics within the two major clades of Atripliceae, the Archiatriplex and Atriplex clades, could explain the unbalanced species richness between them; we found low mean speciation rates in the Archiatriplex clade and one shift to higher speciation rates placed in the branch of the Atriplex core. This acceleration in diversification seems to have started before the transition between C3 and C4 metabolism and before the arrival of Atriplex in the Americas, and matches the Mid-Miocene Climatic Optimum. Besides, the American species of Atriplex exhibit slightly higher net diversification rates than the Australian and Eurasian ones. While time seems not to be associated with diversification, traits such as life form, photosynthetic pathway and plant sex may have played roles as diversification drivers.

CONCLUSIONS

Traits more than time played a key role in Atripliceae diversification, and we could speculate that climate changes could have triggered speciation. The extreme arid or saline environments where Atripliceae species prevail may explain its particular evolutionary trends and trait correlations compared with other angiosperms and highlight the importance of conservation efforts needed to preserve them as genetic resources to deal with climatic changes.

Recent accelerated diversification in rosids occurred outside the tropics

Conflicting relationships have been found between diversification rate and temperature across disparate clades of life. Here, we use a supermatrix comprising nearly 20,000 species of rosids-a clade of ~25% of all angiosperm species-to understand global patterns of diversification and its climatic association. Our approach incorporates historical global temperature, assessment of species' temperature niche, and two broad-scale characterizations of tropical versus non-tropical niche occupancy. We find the diversification rates of most subclades dramatically increased over the last 15 million years (Myr) during cooling associated with global expansion of temperate habitats. Climatic niche is negatively associated with diversification rates, with tropical rosids forming older communities and experiencing speciation rates ~2-fold below rosids in cooler climates. Our results suggest long-term cooling had a disproportionate effect on non-tropical diversification rates, leading to dynamic young communities outside of the tropics, while relative stability in tropical climes led to older, slower-evolving but still species-rich communities.

Estimating rates and patterns of diversification with incomplete sampling: a case study in the rosids

PREMISE

Recent advances in generating large-scale phylogenies enable broad-scale estimation of species diversification. These now common approaches typically are characterized by (1) incomplete species coverage without explicit sampling methodologies and/or (2) sparse backbone representation, and usually rely on presumed phylogenetic placements to account for species without molecular data. We used empirical examples to examine the effects of incomplete sampling on diversification estimation and provide constructive suggestions to ecologists and evolutionary biologists based on those results.

METHODS

We used a supermatrix for rosids and one well-sampled subclade (Cucurbitaceae) as empirical case studies. We compared results using these large phylogenies with those based on a previously inferred, smaller supermatrix and on a synthetic tree resource with complete taxonomic coverage. Finally, we simulated random and representative taxon sampling and explored the impact of sampling on three commonly used methods, both parametric (RPANDA and BAMM) and semiparametric (DR).

RESULTS

We found that the impact of sampling on diversification estimates was idiosyncratic and often strong. Compared to full empirical sampling, representative and random sampling schemes either depressed or inflated speciation rates, depending on methods and sampling schemes. No method was entirely robust to poor sampling, but BAMM was least sensitive to moderate levels of missing taxa.

CONCLUSIONS

We suggest caution against uncritical modeling of missing taxa using taxonomic data for poorly sampled trees and in the use of summary backbone trees and other data sets with high representative bias, and we stress the importance of explicit sampling methodologies in macroevolutionary studies.

Novel phylogenetic methods are needed for understanding gene function in the era of mega-scale genome sequencing

Ongoing large-scale genome sequencing projects are forecasting a data deluge that will almost certainly overwhelm current analytical capabilities of evolutionary genomics. In contrast to population genomics, there are no standardized methods in evolutionary genomics for extracting evolutionary and functional (e.g. gene-trait association) signal from genomic data. Here, we examine how current practices of multi-species comparative genomics perform in this aspect and point out that many genomic datasets are under-utilized due to the lack of powerful methodologies. As a result, many current analyses emphasize gene families for which some functional data is already available, resulting in a growing gap between functionally well-characterized genes/organisms and the universe of unknowns. This leaves unknown genes on the 'dark side' of genomes, a problem that will not be mitigated by sequencing more and more genomes, unless we develop tools to infer functional hypotheses for unknown genes in a systematic manner. We provide an inventory of recently developed methods capable of predicting gene-gene and gene-trait associations based on comparative data, then argue that realizing the full potential of whole genome datasets requires the integration of phylogenetic comparative methods into genomics, a rich but underutilized toolbox for looking into the past.