Continuous and arrested morphological diversification in sister clades of characiform fishes: a phylomorphospace approach

Understanding how and why certain clades diversify greatly in morphology whereas others do not remains a major theme in evolutionary biology. Projecting families of phylogenies into multivariate morphospaces can distinguish two scenarios potentially leading to unequal morphological diversification: unequal magnitude of change per phylogenetic branch, and unequal efficiency in morphological innovation. This approach is demonstrated using a case study of skulls in sister clades within the South American fish superfamily Anostomoidea. Unequal morphological diversification in this system resulted not from the morphologically diverse clade changing more on each phylogenetic branch, but from that clade distributing an equal amount of change more widely through morphospace and innovating continually. Although substantial morphological evolution occurred throughout the less diverse clade's history, most of that clade's expansion in morphospace occurred in the most basal branches, and more derived portions of that radiation oscillated within previously explored limits. Because simulations revealed that there is a maximum 2.7% probability of generating two clades that differ so greatly in the density of lineages within morphospace under a null Brownian model, the observed difference in pattern likely reflects a difference in the underlying evolutionary process. Clade-specific factors that may have promoted or arrested morphological diversification are discussed.     

The molecular phylogenetic signature of clades in decline

Molecular phylogenies have been used to study the diversification of many clades. However, current methods for inferring diversification dynamics from molecular phylogenies ignore the possibility that clades may be decreasing in diversity, despite the fact that the fossil record shows this to be the case for many groups. Here we investigate the molecular phylogenetic signature of decreasing diversity using the most widely used statistic for inferring diversity dynamics from molecular phylogenies, the γ statistic. We show that if a clade is in decline its molecular phylogeny may show evidence of the decrease in the diversification rate that occurred between its diversification and decline phases. The ability to detect the change in diversification rate depends largely on the ratio of the speciation rates of the diversification and decline phases, the higher the ratio the stronger the signal of the change in diversification rate. Consequently, molecular phylogenies of clades in relative rapid decline do not carry a signature of their decreasing diversification. Further, the signal of the change in diversification rate, if present, declines as the diversity drop. Unfortunately, the molecular signature of clades in decline is the same as the signature produced by diversity dependent diversification. Given this similarity, and the inability of current methods to detect declining diversity, it is likely that some of the extant clades that show a decrease in diversification rate, currently interpreted as evidence for diversity dependent diversification, are in fact in decline. Unless methods can be developed that can discriminate between the different modes of diversification, specifically diversity dependent diversification and declining diversity, we will need the fossil record, or data from some other source, to distinguish between these very different diversity trajectories.     

Adaptive radiation and the evolution of nectarivory in a large songbird clade

The accumulation of exceptional ecological diversity within a lineage is a key feature of adaptive radiation resulting from diversification associated with the subdivision of previously underutilized resources. The invasion of unoccupied niche space is predicted to be a key determinant of adaptive diversification, and this process may be particularly important if the diversity of competing lineages within the area, in which the radiation unfolds, is already high. Here, we test whether the evolution of nectarivory resulted in significantly higher rates of morphological evolution, more extensive morphological disparity, and a heightened build‐up of sympatric species diversity in a large adaptive radiation of passerine birds (the honeyeaters, about 190 species) that have diversified extensively throughout continental and insular settings. We find that a large increase in rates of body size evolution and general expansion in morphological space followed an ancestral shift to nectarivory, enabling the build‐up of large numbers of co‐occurring species that vary greatly in size, compared to related and co‐distributed nonnectarivorous clades. These results strongly support the idea that evolutionary shifts into novel areas of niche space play a key role in promoting adaptive radiation in the presence of likely competing lineages.     

Relationships of diversity,disparity, and their evolutionary rates in squirrels (Sciuridae)

Several theories predict that rapidly diversifying clades will also rapidly diverge phenotypically; yet, there are also reasons for suspecting that diversification and divergence might not be correlated. In the widely distributed squirrel clade (Sciuridae), we test for correlations between per lineage speciation rates, species richness, disparity, and a time‐invariant measure of disparity that allows for comparing rates when evolutionary modes differ, as they do in squirrels. We find that species richness and speciation rates are not correlated with clade age or with each other. Disparity appears to be positively correlated with clade age because young, rapidly diversifying Nearctic grassland clades are strongly pulled to a single stable optimum but older, slowly diversifying Paleotropical forest clades contain lineages that diverge along multiple ecological and morphological lines. That contrast is likely due to both the environments they inhabit and their phylogenetic community structure. Our results argue against a shared explanation for diversity and disparity in favor of geographically mediated modes of speciation and ecologically mediated modes of phenotypic evolution.     

Testing for differences in rates of speciation, extinction, and morphological evolution in four tribes of cichlids endemic to Lake Tanganyika, East Africa

Patterns of morphological disparity yield important insight into the causes of diversification and adaptive radiation in East African cichlids. However, comparisons of cichlid disparity have often failed to consider the effects that differing clade ages or stochasticity may have on disparity before making interpretations. Here, a model of branching morphological evolution allows assessment of the relative contributions of differing turnover and morphological change rates, clade ages, and stochastic variation to the observed patterns of disparity in four endemic tribes of Lake Tanganyika cichlids. Simulations compare the likelihood of generating the observed disparity of the four tribes using 200-parameter combinations and four model conditioning variations, which allows inference of evolutionary rate differences among clades. The model is generally robust to model conditioning, the approach to data analysis, and model assumptions. Disparity differences among the first three cichlid tribes, Ectodini, Lamprologini, and Tropheini, can be explained entirely by stochasticity and age, whereas the fourth tribe, Cyprichromini, has likely experienced lower rates of turnover and morphological change. This rate difference is likely related to the low dietary diversity of the Cyprichromini. These results highlight the importance of considering both clade age and stochastic variation when interpreting morphological diversity and evolutionary processes.     

Cryptic diversity in a widespread North American songbird: phylogeography of the Brown Creeper (Certhia americana)

The identification of species via morphological characteristics has traditionally left cryptic species undescribed in taxa under selection for morphological conservation (or a lack of selection for morphological change). Treecreepers (Genus: Certhia) have a conserved morphological appearance, making it difficult to ascertain relationships in the genus based on morphology alone. Recent genetic and song structure studies of Eurasian Treecreepers identified cryptic species within Old World Certhia that were previously undescribed using morphological characteristics. Here, we use mtDNA to investigate cryptic diversity and patterns of diversification in the Brown Creeper (Certhia americana), the single described Certhia species in the Americas. Phylogenetic analyses identified six well-supported geographically-structured clades; the basal divergence separates a northern and a southern lineage in the Brown Creeper, likely cryptic species previously characterized as many subspecies. Sympatry is prevalent between clades in western North America, where possible contact zones warrant further investigation. Allopatry appears to be the primary driver of deep phylogeographic structure within the Brown Creeper; however, within clade diversity is highly correlated with the life history traits of the populations that comprise the geographically structured phylogroups.     

Detecting shifts in diversity limits from molecular phylogenies: what can we know?

Large complete species-level molecular phylogenies can provide the most direct information about the macroevolutionary history of clades having poor fossil records. However, extinction will ultimately erode evidence of pulses of rapid speciation in the deep past. Assessment of how well, and for how long, phylogenies retain the signature of such pulses has hitherto been based on a--probably untenable--model of ongoing diversity-independent diversification. Here, we develop two new tests for changes in diversification 'rules' and evaluate their power to detect sudden increases in equilibrium diversity in clades simulated with diversity-dependent speciation and extinction rates. Pulses of diversification are only detected easily if they occurred recently and if the rate of species turnover at equilibrium is low; rates reported for fossil mammals suggest that the power to detect a doubling of species diversity falls to 50 per cent after less than 50 Myr even with a perfect phylogeny of extant species. Extinction does eventually draw a veil over past dynamics, suggesting that some questions are beyond the limits of inference, but sudden clade-wide pulses of speciation can be detected after many millions of years, even when overall diversity is constrained. Applying our methods to existing phylogenies of mammals and angiosperms identifies intervals of elevated diversification in each.     

ABSOLUTE DIVERSIFICATION RATES IN ANGIOSPERM CLADES

Abstract The extraordinary contemporary species richness and ecological predominance of flowering plants (angiosperms) are even more remarkable when considering the relatively recent onset of their evolutionary diversification. We examine the evolutionary diversification of angiosperms and the observed differential distribution of species in angiosperm clades by estimating the rate of diversification for angiosperms as a whole and for a large set of angiosperm clades. We also identify angiosperm clades with a standing diversity that is either much higher or lower than expected, given the estimated background diversification rate. Recognition of angiosperm clades, the phylogenetic relationships among them, and their taxonomic composition are based on an empirical compilation of primary phylogenetic studies. By making an integrative and critical use of the paleobotanical record, we obtain reasonably secure approximations for the age of a large set of angiosperm clades. Diversification was modeled as a stochastic, time‐homogeneous birth‐and‐death process that depends on the diversification rate (r) and the relative extinction rate (∈). A statistical analysis of the birth and death process was then used to obtain 95% confidence intervals for the expected number of species through time in a clade that diversifies at a rate equal to that of angiosperms as a whole. Confidence intervals were obtained for stem group and for crown group ages in the absence of extinction (∈= 0.0) and under a high relative extinction rate (∈= 0.9). The standing diversity of angiosperm clades was then compared to expected species diversity according to the background rate of diversification, and, depending on their placement with respect to the calculated confidence intervals, exceedingly species‐rich or exceedingly species‐poor clades were identified. The rate of diversification for angiosperms as a whole ranges from 0.077 (∈= 0.9) to 0.089 (∈= 0.0) net speciation events per million years. Ten clades fall above the confidence intervals of expected species diversity, and 13 clades were found to be unexpectedly species poor. The phylogenetic distribution of clades with an exceedingly high number of species suggests that traits that confer high rates of diversification evolved independently in different instances and do not characterize the angiosperms as a whole.     

Large scale mitochondrial sequencing in Mexican Americans suggests a reappraisal of Native American origins

Background

The family Pteropodidae comprises bats commonly known as megabats or Old World fruit bats. Molecular phylogenetic studies of pteropodids have provided considerable insight into intrafamilial relationships, but these studies have included only a fraction of the extant diversity (a maximum of 26 out of the 46 currently recognized genera) and have failed to resolve deep relationships among internal clades. Here we readdress the systematics of pteropodids by applying a strategy to try to resolve ancient relationships within Pteropodidae, while providing further insight into subgroup membership, by 1) increasing the taxonomic sample to 42 genera; 2) increasing the number of characters (to >8,000 bp) and nuclear genomic representation; 3) minimizing missing data; 4) controlling for sequence bias; and 5) using appropriate data partitioning and models of sequence evolution.

Results

Our analyses recovered six principal clades and one additional independent lineage (consisting of a single genus) within Pteropodidae. Reciprocal monophyly of these groups was highly supported and generally congruent among the different methods and datasets used. Likewise, most relationships within these principal clades were well resolved and statistically supported. Relationships among the 7 principal groups, however, were poorly supported in all analyses. This result could not be explained by any detectable systematic bias in the data or incongruence among loci. The SOWH test confirmed that basal branches' lengths were not different from zero, which points to closely-spaced cladogenesis as the most likely explanation for the poor resolution of the deep pteropodid relationships. Simulations suggest that an increase in the amount of sequence data is likely to solve this problem.

Conclusions

The phylogenetic hypothesis generated here provides a robust framework for a revised cladistic classification of Pteropodidae into subfamilies and tribes and will greatly contribute to the understanding of character evolution and biogeography of pteropodids. The inability of our data to resolve the deepest relationships of the major pteropodid lineages suggests an explosive diversification soon after origin of the crown pteropodids. Several characteristics of pteropodids are consistent with this conclusion, including high species diversity, great morphological diversity, and presence of key innovations in relation to their sister group.     

The ecological dynamics of clade diversification and community assembly

Clades diversify in an ecological context, but most macroevolutionary models do not directly encapsulate ecological mechanisms that influence speciation and extinction. A data set of 245 chordate, arthropod, mollusk, and magnoliophyte phylogenies had a majority of clades that showed rapid lineage accumulation early with a slowing more recently, whereas a small but significant minority showed accelerated lineage accumulation in their recent histories. Previous analyses have demonstrated that macroevolutionary birth-death models can replicate the pattern of slowing lineage accumulation only by a strong decrease in speciation rate with increasing species richness and extinction rate held extremely low or absent. In contrast, the metacommunity model presented here could generate the full range of patterns seen in the real phylogenies by simply manipulating the degree of ecological differentiation of new species at the time of speciation. Specifically, the metacommunity model predicts that clades showing decelerating lineage accumulation rates are those that have diversified by ecological modes of speciation, whereas clades showing accelerating lineage accumulation rates are those that have diversified primarily by modes of speciation that generate little or no ecological diversification. A number of testable predictions that integrate data from molecular systematics, community ecology, and biogeography are also discussed.     

Diversification of unicellular eukaryotes: cryptomonad colonizations of marine and fresh waters inferred from revised 18S rRNA phylogeny

The cryptomonads is a well-defined lineage of unicellular eukaryotes, composed of several marine and freshwater groups. However, the evolutionary relationships among these groups are unclear due to conflicting inferences between morphological and molecular phylogenies. Here, we have inferred the evolutionary relationships among marine and freshwater species in order to better understand the importance of the marine-freshwater boundary on the historical diversification patterns of cryptomonads. We have constructed improved molecular phylogenies by taking into account rate variation both across sites and across sequences (covarion substitutions), and by analysing the vast majority of publicly available cryptomonad 18S rRNA sequences and related environmental phylotypes. The resulting phylogenies included 55 sequences, and revealed two novel freshwater cryptomonad clades (CRY1 and CRY2) and a large hidden diversity of cryptomonads. CRY1 was placed deeply within the cryptomonad phylogeny together with all the major freshwater lineages (i.e. Goniomonas and Cryptomonas), while CRY2 was placed within a lineage of marine species identified as Plagioselmis-like with the aid of a new sequence generated from a cultured species. The inferred phylogenies suggest only few successful marine-freshwater transitions over the history of cryptomonads. Most of the transitions seem to have occurred from marine to fresh waters, but re-colonizations of marine habitats have also taken place. This implies that the differences in the biogeophysical conditions between marine and fresh waters constitute a substantial barrier for the cross-colonization of these environments by cryptomonads.     

Vertebral evolution and the diversification of squamate reptiles

Taxonomic, morphological, and functional diversity are often discordant and independent components of diversity. A fundamental and largely unanswered question in evolutionary biology is why some clades diversify primarily in some of these components and not others. Dramatic variation in trunk vertebral numbers (14 to >300) among squamate reptiles coincides with different body shapes, and snake-like body shapes have evolved numerous times. However, whether increased evolutionary rates or numbers of vertebrae underlie body shape and taxonomic diversification is unknown. Using a supertree of squamates including 1375 species, and corresponding vertebral and body shape data, we show that increased rates of evolution in vertebral numbers have coincided with increased rates and disparity in body shape evolution, but not changes in rates of taxonomic diversification. We also show that the evolution of many vertebrae has not spurred or inhibited body shape or taxonomic diversification, suggesting that increased vertebral number is not a key innovation. Our findings demonstrate that lineage attributes such as the relaxation of constraints on vertebral number can facilitate the evolution of novel body shapes, but that different factors are responsible for body shape and taxonomic diversification.     

Lineage diversification and morphological evolution in a large-scale continental radiation: the neotropical ovenbirds and woodcreepers (aves: Furnariidae)

Patterns of diversification in species-rich clades provide insight into the processes that generate biological diversity. We tested different models of lineage and phenotypic diversification in an exceptional continental radiation, the ovenbird family Furnariidae, using the most complete species-level phylogenetic hypothesis produced to date for a major avian clade (97% of 293 species). We found that the Furnariidae exhibit nearly constant rates of lineage accumulation but show evidence of constrained morphological evolution. This pattern of sustained high rates of speciation despite limitations on phenotypic evolution contrasts with the results of most previous studies of evolutionary radiations, which have found a pattern of decelerating diversity-dependent lineage accumulation coupled with decelerating or constrained phenotypic evolution. Our results suggest that lineage accumulation in tropical continental radiations may not be as limited by ecological opportunities as in temperate or island radiations. More studies examining patterns of both lineage and phenotypic diversification are needed to understand the often complex tempo and mode of evolutionary radiations on continents.     

MORPHOLOGICAL DIVERSIFICATION AND ADAPTIVE RADIATION: A COMPARISON OF TWO DIVERSE LIZARD CLADES

We compared the morphological diversity (i.e., the amount of morphological space occupied) of two similar clades, the lizard genera Anolis and Sceloporus. These species-rich monophyletic clades are similar in body size, age of origin, and many aspects of their natural history. We examined a number of morphological traits whose variation is likely to represent adaptation to different aspects of the environment, including body size, limb proportions, head dimensions, and tail length. Examination of the position of species in multidimensional space, based on a principal components analysis, indicates that the morphological diversity of Anolis, which we refer to as disparity, is significantly greater than that of Sceloporus. One potential explanation for this pattern is that morphological diversification in Anolis was facilitated by the evolution of subdigital toe-pads, which allow anoles to use the environment in ways not available to Sceloporus. The geographic location of diversification (tropical and subtropical for Anolis, arid for Sceloporus) may also have been important.     

TESTING FOR ANCIENT ADAPTIVE RADIATIONS IN NEOTROPICAL CICHLID FISHES

Most contemporary studies of adaptive radiation focus on relatively recent and geographically restricted clades. It is less clear whether diversification of ancient clades spanning entire continents is consistent with adaptive radiation. We used novel fossil calibrations to generate a chronogram of Neotropical cichlid fishes and to test whether patterns of lineage and morphological diversification are congruent with hypothesized adaptive radiations in South and Central America. We found that diversification in the Neotropical cichlid clade and the highly diverse tribe Geophagini was consistent with diversity‐dependent, early bursts of divergence followed by decreased rates of lineage accumulation. South American Geophagini underwent early rapid differentiation in body shape, expanding into novel morphological space characterized by elongate‐bodied predators. Divergence in head shape attributes associated with trophic specialization evolved under strong adaptive constraints in all Neotropical cichlid clades. The South American Cichlasomatini followed patterns consistent with constant rates of morphological divergence. Although morphological diversification in South American Heroini was limited, Eocene invasion of Central American habitats was followed by convergent diversification mirroring variation observed in Geophagini. Diversification in Neotropical cichlids was influenced by the early adaptive radiation of Geophagini, which potentially limited differentiation in other cichlid clades.     

Detecting local diversity‐dependence in diversification

Whether there are ecological limits to species diversification is a hotly debated topic. Molecular phylogenies show slowdowns in lineage accumulation, suggesting that speciation rates decline with increasing diversity. A maximum‐likelihood (ML) method to detect diversity‐dependent (DD) diversification from phylogenetic branching times exists, but it assumes that diversity‐dependence is a global phenomenon and therefore ignores that the underlying species interactions are mostly local, and not all species in the phylogeny co‐occur locally. Here, we explore whether this ML method based on the nonspatial diversity‐dependence model can detect local diversity‐dependence, by applying it to phylogenies, simulated with a spatial stochastic model of local DD speciation, extinction, and dispersal between two local communities. We find that type I errors (falsely detecting diversity‐dependence) are low, and the power to detect diversity‐dependence is high when dispersal rates are not too low. Interestingly, when dispersal is high the power to detect diversity‐dependence is even higher than in the nonspatial model. Moreover, estimates of intrinsic speciation rate, extinction rate, and ecological limit strongly depend on dispersal rate. We conclude that the nonspatial DD approach can be used to detect diversity‐dependence in clades of species that live in not too disconnected areas, but parameter estimates must be interpreted cautiously.     

Correlating early evolution of parasitic platyhelminths to Gondwana breakup

Investigating patterns and processes of parasite diversification over ancient geological periods should involve comparisons of host and parasite phylogenies in a biogeographic context. It has been shown previously that the geographical distribution of host-specific parasites of sarcopterygians was guided, from Palaeozoic to Cainozoic times, mostly by evolution and diversification of their freshwater hosts. Here, we propose phylogenies of neobatrachian frogs and their specific parasites (Platyhelminthes, Monogenea) to investigate coevolutionary processes and historical biogeography of polystomes and further discuss all the possible assumptions that may account for the early evolution of these parasites. Phylogenetic analyses of concatenated rRNA nuclear genes (18S and partial 28S) supplemented by cophylogenetic and biogeographic vicariance analyses reveal four main parasite lineages that can be ascribed to centers of diversity, namely Australia, India, Africa, and South America. In addition, the relationships among these biogeographical monophyletic groups, substantiated by molecular dating, reflect sequential origins during the breakup of Gondwana. The Australian polystome lineage may have been isolated during the first stages of the breakup, whereas the Indian lineage would have arisen after the complete separation of western and eastern Gondwanan components. Next, polystomes would have codiverged with hyloid sensu stricto and ranoid frog lineages before the completion of South American and African plate separation. Ultimately, they would have undergone an extensive diversification in South America when their ancestral host families diversified. Therefore, the presence of polystome parasites in specific anuran host clades and in discrete geographic areas reveals the importance of biogeographic vicariance in diversification processes and supports the occurrence and radiation of amphibians over ancient and recent geological periods.     

Understanding angiosperm diversification using small and large phylogenetic trees

How will the emerging possibility of inferring ultra-large phylogenies influence our ability to identify shifts in diversification rate? For several large angiosperm clades (Angiospermae, Monocotyledonae, Orchidaceae, Poaceae, Eudicotyledonae, Fabaceae, and Asteraceae), we explore this issue by contrasting two approaches: (1) using small backbone trees with an inferred number of extant species assigned to each terminal clade and (2) using a mega-phylogeny of 55473 seed plant species represented in GenBank. The mega-phylogeny approach assumes that the sample of species in GenBank is at least roughly proportional to the actual species diversity of different lineages, as appears to be the case for many major angiosperm lineages. Using both approaches, we found that diversification rate shifts are not directly associated with the major named clades examined here, with the sole exception of Fabaceae in the GenBank mega-phylogeny. These agreements are encouraging and may support a generality about angiosperm evolution: major shifts in diversification may not be directly associated with major named clades, but rather with clades that are nested not far within these groups. An alternative explanation is that there have been increased extinction rates in early-diverging lineages within these clades. Based on our mega-phylogeny, the shifts in diversification appear to be distributed quite evenly throughout the angiosperms. Mega-phylogenetic studies of diversification hold great promise for revealing new patterns, but we will need to focus more attention on properly specifying null expectation.     

Rates of morphological evolution are correlated with species richness in salamanders

The tempo and mode of species diversification and phenotypic evolution vary widely across the tree of life, yet the relationship between these processes is poorly known. Previous tests of the relationship between rates of phenotypic evolution and rates of species diversification have assumed that species richness increases continuously through time. If this assumption is violated, simple phylogenetic estimates of net diversification rate may bear no relationship to processes that influence the distribution of species richness among clades. Here, we demonstrate that the variation in species richness among plethodontid salamander clades is unlikely to have resulted from simple time-dependent processes, leading to fundamentally different conclusions about the relationship between rates of phenotypic evolution and species diversification. Morphological evolutionary rates of both size and shape evolution are correlated with clade species richness, but are uncorrelated with simple estimators of net diversification that assume constancy of rates through time. This coupling between species diversification and phenotypic evolution is consistent with the hypothesis that clades with high rates of morphological trait evolution may diversify more than clades with low rates. Our results indicate that assumptions about underlying processes of diversity regulation have important consequences for interpreting macroevolutionary patterns.     

A new method for handling missing species in diversification analysis applicable to randomly or nonrandomly sampled phylogenies

Chronograms from molecular dating are increasingly being used to infer rates of diversification and their change over time. A major limitation in such analyses is incomplete species sampling that moreover is usually nonrandom. While the widely used γ statistic with the Monte Carlo constant-rates test or the birth-death likelihood analysis with the δ AICrc test statistic are appropriate for comparing the fit of different diversification models in phylogenies with random species sampling, no objective automated method has been developed for fitting diversification models to nonrandomly sampled phylogenies. Here, we introduce a novel approach, CorSiM, which involves simulating missing splits under a constant rate birth-death model and allows the user to specify whether species sampling in the phylogeny being analyzed is random or nonrandom. The completed trees can be used in subsequent model-fitting analyses. This is fundamentally different from previous diversification rate estimation methods, which were based on null distributions derived from the incomplete trees. CorSiM is automated in an R package and can easily be applied to large data sets. We illustrate the approach in two Araceae clades, one with a random species sampling of 52% and one with a nonrandom sampling of 55%. In the latter clade, the CorSiM approach detects and quantifies an increase in diversification rate, whereas classic approaches prefer a constant rate model; in the former clade, results do not differ among methods (as indeed expected since the classic approaches are valid only for randomly sampled phylogenies). The CorSiM method greatly reduces the type I error in diversification analysis, but type II error remains a methodological problem.