Punctuational evolution is pervasive in distal site metastatic colonization.

The evolution of metastasis represents a lethal stage of cancer progression. Yet, the evolutionary kinetics of metastatic disease remain unresolved. Here, using single cell CRISPR-Cas9 lineage tracing data, we show that in metastatic disease, gradual molecular evolution is punctuated by episodes of rapid evolutionary change associated with lineage divergence. By measuring punctuational effects across the metastatic cascade, we show that punctuational effects contribute more to the molecular diversity at distal site metastases compared to the paired primary tumor, suggesting qualitatively different modes of evolution may drive primary and metastatic tumor progression. This is the first empirical evidence for distinct patterns of molecular evolution at early and late stages of metastasis and demonstrates the complex interplay of cell intrinsic and extrinsic factors that shape lethal cancer.

Deep-sea origin and depth colonization associated with phenotypic innovations in scleractinian corals.

The deep sea (>200 m) is home to a surprisingly rich biota, which in some cases compares to that found in shallow areas. Scleractinian corals are an example of this - they are key species in both shallow and deep ecosystems. However, what evolutionary processes resulted in current depth distribution of the marine fauna is a long-standing question. Various conflicting hypotheses have been proposed, but few formal tests have been conducted. Here, we use global spatial distribution data to test the bathymetric origin and colonization trends across the depth gradient in scleractinian corals. Using a phylogenetic approach, we infer the origin and historical trends in directionality and speed of colonization during the diversification in depth. We also examine how the emergence of photo-symbiosis and coloniality, scleractinian corals' most conspicuous phenotypic innovations, have influenced this process. Our results strongly support an offshore-onshore pattern of evolution and varying dispersion capacities along depth associated with trait-defined lineages. These results highlight the relevance of the evolutionary processes occurring at different depths to explain the origin of extant marine biodiversity and the consequences of altering these processes by human impact, highlighting the need to include this overlooked evolutionary history in conservation plans.

Terrestrial reproduction and parental care drive rapid evolution in the trade-off between offspring size and number across amphibians.

The trade-off between offspring size and number is central to life history strategies. Both the evolutionary gain of parental care or more favorable habitats for offspring development are predicted to result in fewer, larger offspring. However, despite much research, it remains unclear whether and how different forms of care and habitats drive the evolution of the trade-off. Using data for over 800 amphibian species, we demonstrate that, after controlling for allometry, amphibians with direct development and those that lay eggs in terrestrial environments have larger eggs and smaller clutches, while different care behaviors and adaptations vary in their effects on the trade-off. Specifically, among the 11 care forms we considered at the egg, tadpole and juvenile stage, egg brooding, male egg attendance, and female egg attendance increase egg size; female tadpole attendance and tadpole feeding decrease egg size, while egg brooding, tadpole feeding, male tadpole attendance, and male tadpole transport decrease clutch size. Unlike egg size that shows exceptionally high rates of phenotypic change in just 19 branches of the amphibian phylogeny, clutch size has evolved at exceptionally high rates in 135 branches, indicating episodes of strong selection; egg and tadpole environment, direct development, egg brooding, tadpole feeding, male tadpole attendance, and tadpole transport explain 80% of these events. By explicitly considering diversity in parental care and offspring habitat by stage of offspring development, this study demonstrates that more favorable conditions for offspring development promote the evolution of larger offspring in smaller broods and reveals that the diversity of parental care forms influences the trade-off in more nuanced ways than previously appreciated.

Strong support for a heterogeneous speciation decline model in Dinosauria: a response to claims made by Bonsor et al. (2020).

Through phylogenetic modelling, we previously presented strong support for diversification decline in the three major subclades of dinosaurs (Sakamoto et al. 2016 Proc. Natl Acad. Sci. USA 113, 5036-5040. (doi:10.1073/pnas.1521478113)). Recently, our support for this model has been criticized (Bonsor et al. 2020 R. Soc. Open Sci. 7, 201195. (doi:10.1098/rsos.201195)). Here, we highlight that these criticisms seem to largely stem from a misunderstanding of our study: contrary to Bonsor et al.'s claims, our model accounts for heterogeneity in diversification dynamics, was selected based on deviance information criterion (DIC) scores (not parameter significance), and intercepts were estimated to account for uncertainties in the root age of the phylogenetic tree. We also demonstrate that their new analyses are not comparable to our models: they fit simple, Dinosauria-wide models as a direct comparison to our group-wise models, and their additional trees are subclades that are limited in taxonomic coverage and temporal span, i.e. severely affected by incomplete sampling. We further present results of new analyses on larger, better-sampled trees (N = 961) of dinosaurs, showing support for the time-quadratic model. Disagreements in how we interpret modelled diversification dynamics are to be expected, but criticisms should be based on sound logic and understanding of the model under discussion.

Genes underlying the evolution of tetrapod testes size.

BACKGROUND: Testes vary widely in mass relative to body mass across species, but we know very little about which genes underlie and contribute to such variation. This is partly because evidence for which genes are implicated in testis size variation tends to come from investigations involving just one or a few species. Contemporary comparative phylogenetic methods provide an opportunity to test candidate genes for their role in phenotypic change at a macro-evolutionary scale-across species and over millions of years. Previous attempts to detect genotype-phenotype associations across species have been limited in that they can only detect where genes have driven directional selection (e.g. brain size increase). RESULTS: Here, we introduce an approach that uses rates of evolutionary change to overcome this limitation to test whether any of twelve candidate genes have driven testis size evolution across tetrapod vertebrates-regardless of directionality. We do this by seeking a relationship between the rates of genetic and phenotypic evolution. Our results reveal five genes (Alkbh5, Dmrtb1, Pld6, Nlrp3, Sp4) that each have played unique and complex roles in tetrapod testis size diversity. In all five genes, we find strong significant associations between the rate of protein-coding substitutions and the rate of testis size evolution. Such an association has never, to our knowledge, been tested before for any gene or phenotype. CONCLUSIONS: We describe a new approach to tackle one of the most fundamental questions in biology: how do individual genes give rise to biological diversity? The ability to detect genotype-phenotype associations that have acted across species has the potential to build a picture of how natural selection has sculpted phenotypic change over millions of years.

Temperature thresholds of ecosystem respiration at a global scale.

Ecosystem respiration is a major component of the global terrestrial carbon cycle and is strongly influenced by temperature. The global extent of the temperature-ecosystem respiration relationship, however, has not been fully explored. Here, we test linear and threshold models of ecosystem respiration across 210 globally distributed eddy covariance sites over an extensive temperature range. We find thresholds to the global temperature-ecosystem respiration relationship at high and low air temperatures and mid soil temperatures, which represent transitions in the temperature dependence and sensitivity of ecosystem respiration. Annual ecosystem respiration rates show a markedly reduced temperature dependence and sensitivity compared to half-hourly rates, and a single mid-temperature threshold for both air and soil temperature. Our study indicates a distinction in the influence of environmental factors, including temperature, on ecosystem respiration between latitudinal and climate gradients at short (half-hourly) and long (annual) timescales. Such climatological differences in the temperature sensitivity of ecosystem respiration have important consequences for the terrestrial net carbon sink under ongoing climate change.

Nothing Wrong with the Analysis of Clades in Comparative Evolutionary Studies: A Reply to Poe et al.

In a recent paper, Poe et al. assert that scientists should abandon clade-based approaches, particularly those using named taxonomic ranks. Poe et al. attempt to demonstrate that clade selection can have effects on the results of evolutionary analyses but unfortunately fall short of making any robust conclusions. Here, we demonstrate that the assertions made by Poe et al. have two important flaws: (i) an erroneous view of modern phylogenetic comparative methods; and (ii) a lack of statistical rigor in their analyses. We repeat Poe et al.'s analysis but using appropriate phylogenetic comparative approaches. We demonstrate that results remain consistent regardless of the clade definition. We go on to discuss the value of taxonomic groupings and how they can provide meaningful units of comparison in evolutionary study. Unlike the disheartening suggestion to abandon the use of clades, scientists can instead continue to use phylogenetic " corrections" that are already the standard for most comparative evolutionary analyses. [Comparative methods; evolution; phylogeny; taxonomy.].

150 million years of sustained increase in pterosaur flight efficiency.

The long-term accumulation of biodiversity has been punctuated by remarkable evolutionary transitions that allowed organisms to exploit new ecological opportunities. Mesozoic flying reptiles (the pterosaurs), which dominated the skies for more than 150 million years, were the product of one such transition. The ancestors of pterosaurs were small and probably bipedal early archosaurs1, which were certainly well-adapted to terrestrial locomotion. Pterosaurs diverged from dinosaur ancestors in the Early Triassic epoch (around 245 million years ago); however, the first fossils of pterosaurs are dated to 25 million years later, in the Late Triassic epoch. Therefore, in the absence of proto-pterosaur fossils, it is difficult to study how flight first evolved in this group. Here we describe the evolutionary dynamics of the adaptation of pterosaurs to a new method of locomotion. The earliest known pterosaurs took flight and subsequently appear to have become capable and efficient flyers. However, it seems clear that transitioning between forms of locomotion2,3-from terrestrial to volant-challenged early pterosaurs by imposing a high energetic burden, thus requiring flight to provide some offsetting fitness benefits. Using phylogenetic statistical methods and biophysical models combined with information from the fossil record, we detect an evolutionary signal of natural selection that acted to increase flight efficiency over millions of years. Our results show that there was still considerable room for improvement in terms of efficiency after the appearance of flight. However, in the Azhdarchoidea4, a clade that exhibits gigantism, we test the hypothesis that there was a decreased reliance on flight5-7 and find evidence for reduced selection on flight efficiency in this clade. Our approach offers a blueprint to objectively study functional and energetic changes through geological time at a more nuanced level than has previously been possible.

Rapid decreases in relative testes mass among monogamous birds but not in other vertebrates.

Larger testes produce more sperm and therefore improve reproductive success in the face of sperm competition. Adaptation to social mating systems with relatively high and low sperm competition are therefore likely to have driven changes in relative testes size in opposing directions. Here, we combine the largest vertebrate testes mass dataset ever collected with phylogenetic approaches for measuring rates of morphological evolution to provide the first quantitative evidence for how relative testes mass has changed over time. We detect explosive radiations of testes mass diversity distributed throughout the vertebrate tree of life: bursts of rapid change have been frequent during vertebrate evolutionary history. In socially monogamous birds, there have been repeated rapid reductions in relative testes mass. We see no such pattern in other monogamous vertebrates; the prevalence of monogamy in birds may have increased opportunities for investment in alternative behaviours and physiologies allowing reduced investment in expensive testes.

The decoupled nature of basal metabolic rate and body temperature in endotherm evolution.

The origins of endothermy in birds and mammals are important events in vertebrate evolution. Endotherms can maintain their body temperature (Tb) over a wide range of ambient temperatures primarily using the heat that is generated continuously by their high basal metabolic rate (BMR)1. There is also an important positive feedback loop as Tb influences BMR1-3. Owing to this interplay between BMRs and Tb, many ecologists and evolutionary physiologists posit that the evolution of BMR and Tb must have been coupled during the radiation of endotherms3-5, changing with similar trends6-8. However, colder historical environments might have imposed strong selective pressures on BMR to compensate for increased rates of heat loss and to keep Tb constant9-12. Thus, adaptation to cold ambient temperatures through increases in BMR could have decoupled BMR from Tb and caused different evolutionary routes to the modern diversity in these traits. Here we show that BMR and Tb were decoupled in approximately 90% of mammalian phylogenetic branches and 36% of avian phylogenetic branches. Mammalian BMRs evolved with rapid bursts but without a long-term directional trend, whereas Tb evolved mostly at a constant rate and towards colder bodies from a warmer-bodied common ancestor. Avian BMRs evolved predominantly at a constant rate and without a long-term directional trend, whereas Tb evolved with much greater rate heterogeneity and with adaptive evolution towards colder bodies. Furthermore, rapid shifts that lead to both increases and decreases in BMRs were linked to abrupt changes towards colder ambient temperatures-although only in mammals. Our results suggest that natural selection effectively exploited the diversity in mammalian BMRs under diverse, often-adverse historical thermal environments.

Extreme and rapid bursts of functional adaptations shape bite force in amniotes.

Adaptation is the fundamental driver of functional and biomechanical evolution. Accordingly, the states of biomechanical traits (absolute or relative trait values) have long been used as proxies for adaptations in response to direct selection. However, ignoring evolutionary history, in particular ancestry, passage of time and the rate of evolution, can be misleading. Here, we apply a recently developed phylogenetic statistical approach using significant rate shifts to detect instances of exceptional rates of adaptive changes in bite force in a large group of terrestrial vertebrates, the amniotes. Our results show that bite force in amniotes evolved through multiple bursts of exceptional rates of adaptive changes, whereby whole groups-including Darwin's finches, maniraptoran dinosaurs (group of non-avian dinosaurs including birds), anthropoids and hominins (fossil and modern humans)-experienced significant rate increases compared to the background rate. However, in most parts of the amniote tree of life, we find no exceptional rate increases, indicating that coevolution with body size was primarily responsible for the patterns observed in bite force. Our approach represents a template for future studies in functional morphology and biomechanics, where exceptional rates of adaptive changes can be quantified and potentially linked to specific ecological factors underpinning major evolutionary radiations.

Rapid Change in Mammalian Eye Shape Is Explained by Activity Pattern.

The rate of morphological evolution along the branches of a phylogeny varies widely [1-6]. Although such rate variation is often assumed to reflect the strength of historical natural selection resulting in adaptation [7-14], this lacks empirical and analytical evidence. One way to demonstrate a relationship between branchwise rates and adaptation would be to show that rapid rates of evolution are linked with ecological shifts or key innovations. Here, we test for this link by determining whether activity pattern, the time of day at which species are active, explains rapid bursts of evolutionary change in eye shape. Using modern approaches to identify shifts in the rate of morphological evolution [7, 13], we find that over 74% of rapid eye-shape change during mammalian evolutionary history is directly explained by distinct selection pressures acting on nocturnal, cathemeral, and diurnal species. Our results reveal how ecological changes occurring along the branches of a phylogeny can manifest in subsequent changes in the rate of morphological evolution. Although selective pressures exerted by different activity patterns have acted uniformly across all mammals, we find differences in the rate of eye-shape evolution among orders. The key to understanding this is in how ecology itself has evolved. We find heterogeneity in how activity pattern has evolved among mammals that ultimately led to differences in the rate of eye-shape evolution among species. Our approach represents an exciting new way to pinpoint factors driving adaptation, enabling a clearer understanding of the factors that drive the evolution of biological diversity.

Transitions between foot postures are associated with elevated rates of body size evolution in mammals.

Terrestrial mammals have evolved various foot postures: flat-footed (plantigrady), tiptoed (digitigrady), and hooved (unguligrady) postures. Although the importance of foot posture on ecology and body size of mammalian species has been widely recognized, its evolutionary trajectory and influence on body size evolution across mammalian phylogeny remain untested. Taking a Bayesian phylogenetic approach combined with a comprehensive dataset of foot postures in 880 extant mammalian species, we investigated the evolutionary history of foot postures and rates of body size evolution, within the same posture and at transitions between postures. Our results show that the common ancestor of mammals was plantigrade, and transitions predominantly occurred only between plantigrady and digitigrady and between digitigrady and unguligrady. At the transitions between plantigrady and digitigrady and between digitigrady and unguligrady, rates of body size evolution are significantly elevated leading to the larger body masses of digitigrade species (∼1 kg) and unguligrade species (∼78 kg) compared with their respective ancestral postures [plantigrady (∼0.75 kg) and digitigrady]. Our results demonstrate the importance of foot postures on mammalian body size evolution and have implications for mammalian body size increase through time. In addition, we highlight a way forward for future studies that seek to integrate morphofunctional and macroevolutionary approaches.

Phylogenetic non-independence in rates of trait evolution.

Statistical non-independence of species' biological traits is recognized in most traits under selection. Yet, whether or not the evolutionary rates of such biological traits are statistically non-independent remains to be tested. Here, we test the hypothesis that phenotypic evolutionary rates are non-independent, i.e. contain phylogenetic signal, using empirical rates of evolution in three separate traits: body mass in mammals, beak shape in birds and bite force in amniotes. Specifically, we test if evolutionary rates are phylogenetically interdependent. We find evidence for phylogenetic signal in evolutionary rates in all three case studies. While phylogenetic signal diminishes deeper in time, this is reflective of statistical power owing to small sample and effect sizes. When effect size is large, e.g. owing to the presence of fossil tips, we detect high phylogenetic signals even in deeper time slices. Thus, we recommend that rates be treated as being non-independent throughout the evolutionary history of the group of organisms under study, and any summaries or analyses of rates through time-including associations of rates with traits-need to account for the undesired effects of shared ancestry.

Dinosaurs reveal the geographical signature of an evolutionary radiation.

Dinosaurs dominated terrestrial ecosystems across the globe for over 100 million years and provide a classic example of an evolutionary radiation. However, little is known about how these animals radiated geographically to become globally distributed. Here, we use a biogeographical model to reconstruct the dinosaurs' ancestral locations, revealing the spatial mechanisms that underpinned this 170-million-year-long radiation. We find that dinosaurs spread rapidly initially, followed by a significant continuous and gradual reduction in their speed of movement towards the Cretaceous/Tertiary boundary (66 million years ago). This suggests that the predominant mode of dinosaur speciation changed through time with speciation originally largely driven by geographical isolation-when dinosaurs speciated more, they moved further. This was gradually replaced by increasing levels of sympatric speciation (species taking advantage of ecological opportunities within their existing environment) as terrestrial space became a limiting factor. Our results uncover the geographical signature of an evolutionary radiation.

A cautionary note on the use of Ornstein Uhlenbeck models in macroevolutionary studies.

Phylogenetic comparative methods are increasingly used to give new insights into the dynamics of trait evolution in deep time. For continuous traits the core of these methods is a suite of models that attempt to capture evolutionary patterns by extending the Brownian constant variance model. However, the properties of these models are often poorly understood, which can lead to the misinterpretation of results. Here we focus on one of these models - the Ornstein Uhlenbeck (OU) model. We show that the OU model is frequently incorrectly favoured over simpler models when using Likelihood ratio tests, and that many studies fitting this model use datasets that are small and prone to this problem. We also show that very small amounts of error in datasets can have profound effects on the inferences derived from OU models. Our results suggest that simulating fitted models and comparing with empirical results is critical when fitting OU and other extensions of the Brownian model. We conclude by making recommendations for best practice in fitting OU models in phylogenetic comparative analyses, and for interpreting the parameters of the OU model.

Comparison of conservation metrics in a case study of lemurs.

Conservation planning is important to protect species from going extinct now that natural habitats are decreasing owing to human activity and climate change. However, there is considerable controversy in choosing appropriate metrics to weigh the value of species and geographic regions. For example, the added value of phylogenetic conservation-selection criteria remains disputed because high correlations between them and the nonphylogenetic criteria of species richness have been reported. We evaluated the commonly used conservation metrics species richness, endemism, phylogenetic diversity (PD), and phylogenetic endemism (PE) in a case study on lemurs of Madagascar. This enabled us to identify the conservation target of each metric and consider how they may be used in future conservation planning. We also devised a novel metric that uses a phylogeny scaled according to the rate of phenotypic evolution as a proxy for a species' ability to adapt to change. High rates of evolution may indicate generalization or specialization. Both specialization and low rates of evolution may result in an inability to adapt to changing environments. We examined conservation priorities by using the inverse of the rate of body mass evolution to account for species with low rates of evolution. In line with previous work, we found high correlations among species richness and PD (r = 0.96), and endemism and PE (r = 0.82) in Malagasy lemurs. Phylogenetic endemism in combination with rates of evolution and their inverse prioritized grid cells containing highly endemic and specialized lemurs at risk of extinction, such as Avahi occidentalis and Lepilemur edwardsi, 2 endangered lemurs with high rates of phenotypic evolution and low-quality diets, and Hapalemur aureus, a critically endangered species with a low rate of body mass evolution and a diet consisting of very high doses of cyanide.

Dinosaurs in decline tens of millions of years before their final extinction.

Whether dinosaurs were in a long-term decline or whether they were reigning strong right up to their final disappearance at the Cretaceous-Paleogene (K-Pg) mass extinction event 66 Mya has been debated for decades with no clear resolution. The dispute has continued unresolved because of a lack of statistical rigor and appropriate evolutionary framework. Here, for the first time to our knowledge, we apply a Bayesian phylogenetic approach to model the evolutionary dynamics of speciation and extinction through time in Mesozoic dinosaurs, properly taking account of previously ignored statistical violations. We find overwhelming support for a long-term decline across all dinosaurs and within all three dinosaurian subclades (Ornithischia, Sauropodomorpha, and Theropoda), where speciation rate slowed down through time and was ultimately exceeded by extinction rate tens of millions of years before the K-Pg boundary. The only exceptions to this general pattern are the morphologically specialized herbivores, the Hadrosauriformes and Ceratopsidae, which show rapid species proliferations throughout the Late Cretaceous instead. Our results highlight that, despite some heterogeneity in speciation dynamics, dinosaurs showed a marked reduction in their ability to replace extinct species with new ones, making them vulnerable to extinction and unable to respond quickly to and recover from the final catastrophic event.