Evolution of Cranial Shape in Caecilians (Amphibia: Gymnophiona)

Insights into morphological diversification can be obtained from the ways the species of a clade occupy morphospace. Projecting a phylogeny into morphospace provides estimates of evolutionary trajectories as lineages diversified information that can be used to infer the dynamics of evolutionary processes that produced patterns of morphospace occupation. We present here a large-scale investigation into evolution of morphological variation in the skull of caecilian amphibians, a major clade of vertebrates. Because caecilians are limbless, predominantly fossorial animals, diversification of their skull has occurred within a framework imposed by the functional demands of head-first burrowing. We examined cranial shape in 141 species, over half of known species, using X-ray computed tomography and geometric morphometrics. Mapping an existing phylogeny into the cranial morphospace to estimate the history of morphological change (phylomorphospace), we find a striking pattern: most species occupy distinct clusters in cranial morphospace that closely correspond to the main caecilian clades, and each cluster is separated by unoccupied morphospace. The empty spaces in shape space are unlikely to be caused entirely by extinction or incomplete sampling. The main caecilian clades have different amounts of morphological disparity, but neither clade age nor number of species account for this variation. Cranial shape variation is clearly linked to phyletic divergence, but there is also homoplasy, which is attributed to extrinsic factors associated with head-first digging: features of caecilian crania that have been previously argued to correlate with differential microhabitat use and burrowing ability, such as subterminal and terminal mouths, degree of temporal fenestration (stegokrotaphy/zygokrotaphy), and eyes covered by bone, have evolved and many combinations occur in modern species. We find evidence of morphological convergence in cranial shape, among species that have eyes covered by bone, resulting in a narrow bullet-shaped head. These results reveal a complex history, including early expansion of morphospace and both divergent and convergent evolution resulting in the diversity we observe today.     

MORPHOLOGICAL RATES OF ANGIOSPERM SEED SIZE EVOLUTION

The evolution of seed size among angiosperms reflects their ecological diversification in a complex fitness landscape of life‐history strategies. The lineages that have evolved seeds beyond the upper and lower boundaries that defined nonflowering seed plants since the Paleozoic are more dispersed across the angiosperm phylogeny than would be expected under a neutral model of phenotypic evolution. Morphological rates of seed size evolution estimated for 40 clades based on 17,375 species ranged from 0.001 (Garryales) to 0.207 (Malvales). Comparative phylogenetic analysis indicated that morphological rates are not associated with the clade's seed size but are negatively correlated with the clade's position in the overall distribution of angiosperm seed sizes; clades with seed sizes closer to the angiosperm mean had significantly higher morphological rates than clades with extremely small or extremely large seeds. Likewise, per‐clade taxonomic diversification rates are not associated with the seed size of the clade but with where the clade falls within the angiosperm seed size distribution. These results suggest that evolutionary rates (morphological and taxonomic) are elevated in densely occupied regions of the seed morphospace relative to lineages whose ecophenotypic innovations have moved them toward the edges.     

Species and shape diversification are inversely correlated among gobies and cardinalfishes (Teleostei: Gobiiformes)

Gobies and their relatives are significant components of nearshore marine, estuarine, and freshwater fish faunas in both tropical and temperate habitats worldwide. They are remarkable for their ability to adapt to and diversify in a wide range of environments. Among gobiiform clades, species diversities vary widely, ranging from two species in Kurtidae to more than 1,000 species in Gobiidae. There is also great variation in head and body shape and in environmental preferences (fresh, brackish, or marine habitats). In this study, I used a time-calibrated molecular phylogeny, coupled with morphometric and comparative analyses, to examine evolutionary rates of both speciation and morphological diversification among gobiiform lineages. Projection of the phylogeny onto a shape-derived morphospace shows that Gobioidei is morphometrically distinct from its sister taxon Apogonoidei, but that families within Gobioidei overlap in morphospace. Analysis of species diversification rates indicates that three rate shifts have occurred over the evolutionary history of Gobiiformes. Relative to the other lineages, Kurtidae has exhibited a slowdown in speciation, whereas both Apogonidae and Gobiidae?+?Gobionellidae have experienced an increase in diversification. Comparative analyses show that in Apogonidae and Gobiidae?+?Gobionellidae, increased speciation is correlated with diminished rates of morphological diversification, differently manifested in either clade and among the various sublineages. The elevation in speciation rates and diminishment in rates of morphological change in both Apogonidae and the clade Gobiidae?+?Gobionellidae are correlated with shifts to oceanic habitats from freshwater. This pattern is the complement to that seen across the global radiation of acanthomorph fishes in which a decrease in species diversification is associated with an increase in morphological disparity.     

AN ADAPTIVE RADIATION OF FROGS IN A SOUTHEAST ASIAN ISLAND ARCHIPELAGO

Living amphibians exhibit a diversity of ecologies, life histories, and species‐rich lineages that offers opportunities for studies of adaptive radiation. We characterize a diverse clade of frogs (Kaloula, Microhylidae) in the Philippine island archipelago as an example of an adaptive radiation into three primary habitat specialists or ecotypes. We use a novel phylogenetic estimate for this clade to evaluate the tempo of lineage accumulation and morphological diversification. Because species‐level phylogenetic estimates for Philippine Kaloula are lacking, we employ dense population sampling to determine the appropriate evolutionary lineages for diversification analyses. We explicitly take phylogenetic uncertainty into account when calculating diversification and disparification statistics and fitting models of diversification. Following dispersal to the Philippines from Southeast Asia, Kaloula radiated rapidly into several well‐supported clades. Morphological variation within Kaloula is partly explained by ecotype and accumulated at high levels during this radiation, including within ecotypes. We pinpoint an axis of morphospace related directly to climbing and digging behaviors and find patterns of phenotypic evolution suggestive of ecological opportunity with partitioning into distinct habitat specialists. We conclude by discussing the components of phenotypic diversity that are likely important in amphibian adaptive radiations.     

Disparity and Evolutionary Rate Do Not Explain Diversity Patterns in Muroid Rodents (Rodentia: Muroidea)

A positive correlation between diversity and disparity/evolutionary rate is predicted by multiple evolutionary theories. However, recent empirical studies in various taxa do not always find such an association. Similarly, we find no correlation between these two levels of variation, based on cranial morphometric data and molecular phylogenetic data from 317 muroid rodent species and dipodoid outgroups, analyzed using three-dimensional geometric morphometrics. This disassociation was found using both phylogenetic and non-phylogenetic approaches, indicating that an increase in clade richness is not necessarily followed by an increase in morphological divergence and vice versa. Furthermore, the distribution of muroid families in morphospace is highly overlapping suggesting greater variation within than between clades. Taken together with the observation that families with the most distinctive cranial morphologies (nesomyids, dipodids, and spalacids) are the least diverse, indicates that evolution of new cranial morphologies may not play an important role in the diversification of muroid rodents.     

Ecological variation in South American geophagine cichlids arose during an early burst of adaptive morphological and functional evolution

Diversity and disparity are unequally distributed both phylogenetically and geographically. This uneven distribution may be owing to differences in diversification rates between clades resulting from processes such as adaptive radiation. We examined the rate and distribution of evolution in feeding biomechanics in the extremely diverse and continentally distributed South American geophagine cichlids. Evolutionary patterns in multivariate functional morphospace were examined using a phylomorphospace approach, disparity-through-time analyses and by comparing Brownian motion (BM) and adaptive peak evolutionary models using maximum likelihood. The most species-rich and functionally disparate clade (CAS) expanded more efficiently in morphospace and evolved more rapidly compared with both BM expectations and its sister clade (GGD). Members of the CAS clade also exhibited an early burst in functional evolution that corresponds to the development of modern ecological roles and may have been related to the colonization of a novel adaptive peak characterized by fast oral jaw mechanics. Furthermore, reduced ecological opportunity following this early burst may have restricted functional evolution in the GGD clade, which is less species-rich and more ecologically specialized. Patterns of evolution in ecologically important functional traits are consistent with a pattern of adaptive radiation within the most diverse clade of Geophagini.     

Testing for unequal rates of morphological diversification in the absence of a detailed phylogeny: a case study from Characiform fishes

This study develops the random phylogenies rate test (RAPRATE), a likelihood method that simulates morphological evolution along randomly generated phylogenies, and uses it to determine whether a considerable difference in morphological diversity between two sister clades of South American fishes should be taken as evidence of differing rates of morphological change or lineage turnover. Despite identical ages of origin, similar species richness, and sympatric geographic distributions, the morphological and ecological diversity of the superfamily Anostomoidea exceeds that of the Curimatoidea. The test shows with 90% confidence (using variance among species as the measure of morphological diversity) or 99% confidence (using volume of occupied morphospace) that the rate of morphological change per unit time in the Anostomoidea likely exceeded that of the Curimatoidea. Variation in the rate of lineage turnover (speciation and extinction rates) is not found to affect greatly the morphological diversity of simulated clades and is not a likely explanation of the observed difference in morphological diversity in this case study. Though a 17% or greater delay in the onset of diversification in the Curimatoidea remains a possible alternative explanation of unequal morphological diversification, further simulations suggest that two clades drawn from the possible treespace of the Anostomoidea and Curimatoidea will rarely differ so greatly in the onset of diversification. Several uniquely derived morphological and ecological features of the Anostomoidea and Curimatoidea may have accelerated or decelerated their rate of morphological change, including a marked lengthening of the quadrate that may have relaxed structural constraints on the evolution of the anostomoid jaw.     

A major shift in diversification rate helps explain macroevolutionary patterns in primate species diversity

Primates represent one of the most species rich, wide ranging, and ecologically diverse clades of mammals. What major macroevolutionary factors have driven their diversification and contributed to the modern distribution of primate species remains widely debated. We employed phylogenetic comparative methods to examine the role of clade age and evolutionary rate heterogeneity in the modern distribution of species diversity of Primates. Primate diversification has accelerated since its origin, with decreased extinction leading to a shift to even higher evolutionary rates in the most species rich family (Cercopithecidae). Older primate clades tended to be more diverse, however a shift in evolutionary rate was necessary to adequately explain the imbalance in species diversity. Species richness was also poorly explained by geographic distribution, especially once clade age and evolutionary rate shifts were accounted for, and may relate instead to other ecological factors. The global distribution of primate species diversity appears to have been strongly impacted by heterogeneity in evolutionary rates.     

Discerning evolutionary processes in patterns of tamarin (genus Saguinus) craniofacial variation

Quantitative genetic theory specifies evolutionary expectations for morphological diversification by genetic drift in a monophyletic clade. If genetic drift is responsible for the evolutionary morphological diversification of a clade, patterns of within- and between-taxon morphological variance/covariance should be proportional. We tested for proportionality of within- and between-species craniofacial morphological variation in 12 species of tamarins (genus Saguinus). We found that within- and between-taxon morphological variations across the entire genus were not proportional, and hence not likely to be due to genetic drift alone. The primary deviation from proportionality is that size and size-related shape in the cranium is more variable relative to other aspects of cranial morphology than expected under genetic drift, suggesting differential size selection between the two major clades, the small-bodied and large-bodied tamarins. Within each of these major clades, most of the interspecific variation is consistent with the pattern expected under genetic drift, although specific contrasts may indicate the involvement of differential selection. Morphological distances among taxa do not correspond very closely to the phylogeny derived from mtDNA. In particular, S. oedipus and S. geoffroyi are very distinct morphologically from the rest of the tamarins, although they are phylogenetically the sister clade to a clade containing S. midas and S. bicolor. Morphological similarity is not a good guide to phylogenetic affinity in the tamarins, especially with regard to deeper nodes in the phylogenetic tree.     

Do different disparity proxies converge on a common signal? Insights from the cranial morphometrics and evolutionary history of Pterosauria (Diapsida: Archosauria)

Disparity, or morphological diversity, is often quantified by evolutionary biologists investigating the macroevolutionary history of clades over geological timescales. Disparity is typically quantified using proxies for morphology, such as measurements, discrete anatomical characters, or geometric morphometrics. If different proxies produce differing results, then the accurate quantification of disparity in deep time may be problematic. However, despite this, few studies have attempted to examine disparity of a single clade using multiple morphological proxies. Here, as a case study for this question, we examine the disparity of the volant Mesozoic fossil reptile clade Pterosauria, an intensively studied group that achieved substantial morphological, ecological and taxonomic diversity during their 145+ million-year evolutionary history. We characterize broadscale patterns of cranial morphological disparity for pterosaurs for the first time using landmark-based geometric morphometrics and make comparisons to calculations of pterosaur disparity based on alternative metrics. Landmark-based disparity calculations suggest that monofenestratan pterosaurs were more diverse cranially than basal non-monofenestratan pterosaurs (at least when the aberrant anurognathids are excluded), and that peak cranial disparity may have occurred in the Early Cretaceous, relatively late in pterosaur evolution. Significantly, our cranial disparity results are broadly congruent with those based on whole skeleton discrete character and limb proportion data sets, indicating that these divergent approaches document a consistent pattern of pterosaur morphological evolution. Therefore, pterosaurs provide an exemplar case demonstrating that different proxies for morphological form can converge on the same disparity signal, which is encouraging because often only one such proxy is available for extinct clades represented by fossils. Furthermore, mapping phylogeny into cranial morphospace demonstrates that pterosaur cranial morphology is significantly correlated with, and potentially constrained by, phylogenetic relationships.     

Ecological limits on diversification of the himalayan core corvoidea

Within regions, differences in the number of species among clades must be explained by clade age, net diversification rate, or immigration. We examine these alternatives by assessing historical causes of the low diversity of a bird parvorder in the Himalayas (the core Corvoidea, 57 species present), relative to its more species rich sister clade (the Passerida, ～400 species present), which together comprise the oscine passerines within this region. The core Corvoidea contain ecologically diverse species spanning a large range of body sizes and elevations. Despite this diversity, on the basis of ecological, morphological, and phylogenetic information, we infer that the best explanation for the low number of species within the Himalayan core Corvoidea is one in which ecology limits diversification and/or immigration. Within the core Corvoidea, body size is correlated with elevation: large species are found at high elevations, and small species at lower elevations. This contrasts with the presence of many small-bodied species spanning all elevations in the Passerida and many large bodied species at low elevations in the other orders of birds (the nonpasserines). Cladogenetic events leading to ecological differences between species in body size and shape mostly occurred millions of years ago, and the rate of evolutionary change has declined toward the present. Elevational distributions have been evolutionarily more labile, but are also associated with ancient cladogenetic events. We suggest the core Corvoidea occupy a restricted volume of ecological space in competition with other bird species, and this has limited in situ diversification and/or immigration.     

The role of immigration and in-situ radiation in explaining blood parasite assemblages in an island bird clade

Parasite communities on islands are assembled through multiple immigrations and/or in-situ diversification. In this study, we used a phylogenetic approach to investigate the role of such processes in shaping current patterns of diversity in Leucocytozoon, a group of haemosporidian blood parasites infecting whites eyes (Zosterops) endemic to the Mascarene archipelago (south-western Indian Ocean). We found that this parasite community arose through a combination of multiple immigrations and in-situ diversification, highlighting the importance of both processes in explaining island diversity. Specifically, two highly diverse parasite clades appear to have been present in the Mascarenes for most of their evolutionary history and have diversified within the archipelago, while another lineage apparently immigrated more recently, probably with human-introduced birds. Interestingly, the evolutionary histories of one clade of parasites and Indian Ocean Zosterops seem tightly associated with a significant signal for phylogenetic congruence, suggesting that host-parasite co-divergence may have occurred in this system.     

A morphospace‐based test for competitive exclusion among flying vertebrates: did birds,bats and pterosaurs get in each other's space?

Three vertebrate groups – birds, bats and pterosaurs – have evolved flapping flight over the past 200 million years. This innovation allowed each clade access to new ecological opportunities, but did the diversification of one of these groups inhibit the evolutionary radiation of any of the others? A related question is whether having the wing attached to the hindlimbs in bats and pterosaurs constrained their morphological diversity relative to birds. Fore‐ and hindlimb measurements from 894 specimens were used to construct a morphospace to assess morphological overlap and range, a possible indicator of competition, among the three clades. Neither birds nor bats entered pterosaur morphospace across the Cretaceous–Paleogene (Tertiary) extinction. Bats plot in a separate area from birds, and have a significantly smaller morphological range than either birds or pterosaurs. On the basis of these results, competitive exclusion among the three groups is not supported.     

Diversity Dynamics in Nymphalidae Butterflies: Effect of Phylogenetic Uncertainty on Diversification Rate Shift Estimates

The species rich butterfly family Nymphalidae has been used to study evolutionary interactions between plants and insects. Theories of insect-hostplant dynamics predict accelerated diversification due to key innovations. In evolutionary biology, analysis of maximum credibility trees in the software MEDUSA (modelling evolutionary diversity using stepwise AIC) is a popular method for estimation of shifts in diversification rates. We investigated whether phylogenetic uncertainty can produce different results by extending the method across a random sample of trees from the posterior distribution of a Bayesian run. Using the MultiMEDUSA approach, we found that phylogenetic uncertainty greatly affects diversification rate estimates. Different trees produced diversification rates ranging from high values to almost zero for the same clade, and both significant rate increase and decrease in some clades. Only four out of 18 significant shifts found on the maximum clade credibility tree were consistent across most of the sampled trees. Among these, we found accelerated diversification for Ithomiini butterflies. We used the binary speciation and extinction model (BiSSE) and found that a hostplant shift to Solanaceae is correlated with increased net diversification rates in Ithomiini, congruent with the diffuse cospeciation hypothesis. Our results show that taking phylogenetic uncertainty into account when estimating net diversification rate shifts is of great importance, as very different results can be obtained when using the maximum clade credibility tree and other trees from the posterior distribution.     

Large-scale phylogenetic analysis provides insights into the diversification and evolution of sessilid peritrich ciliates (Protista: Ciliophora)

The Peritrichia is a speciose and morphologically distinctive assemblage of ciliated protists that was first observed by Antonie van Leeuwenhoek over 340 years ago. In the last two decades, the phylogenetic relationships of this group have been increasingly debated as morphological and molecular analyses have generated contrasting conclusions, mainly owing to limited sampling. In the present study, we performed expanded phylogenetic analyses of 152 sessilid peritrichs collected from 14 different provinces of China and 141 SSU rDNA peritrich sequences from GenBank. The results of the analyses revealed new divergent relationships between and within major clades that challenge the morphological classification of this group including, (1) the recovery of four major phylogenetically divergent clades in the monophyletic order Sessilida, (2) aboral structures such as the stalk and spasmoneme were evolutionary labile, (3) the stalk or/and spasmoneme was lost in each divergent clade indicating that parallel evolution occurred in sessilid peritrichs and (4) the life cycle and habit drive the diversity of aboral structures as well as diversification and evolution in peritrichs.     

Simultaneous Quaternary radiations of three damselfly clades across the Holarctic

If climate change during the Quaternary shaped the macroevolutionary dynamics of a taxon, we expect to see three features in its history: elevated speciation or extinction rates should date to this time, more northerly distributed clades should show greater discontinuities in these rates, and similar signatures of those effects should be evident in the phylogenetic and phylodemographic histories of multiple clades. In accordance with the role of glacial cycles, speciation rates increased in the Holarctic Enallagma damselflies during the Quaternary, with a 4.25x greater increase in a more northerly distributed clade as compared with a more southern clade. Finer-scale phylogenetic analyses of three radiating clades within the northern clade show similar, complex recent histories over the past 250,000 years to produce 17 Nearctic and four Palearctic extant species. All three are marked by nearly synchronous deep splits that date to approximately 250,000 years ago, resulting in speciation in two. This was soon followed by significant demographic expansions in at least two of the three clades. In two, these expansions seem to have preceded the radiations that have given rise to most of the current biodiversity. Each also produced species at the periphery of the clade's range. In spite of clear genetic support for reproductive isolation among almost all species, mtDNA signals of past asymmetric hybridization between species in different clades also suggest a role for the evolution of mate choice in generating reproductive isolation as species recolonized the landscape following deglaciation. These analyses suggest that recent climate fluctuations resulted in radiations driven by similar combinations of speciation processes acting in different lineages.     

A preliminary classification of evolutionary radiations

Rapid increases in taxonomic diversity are generally described as adaptive or evolutionary radiations. Such radiations differ widely in the rate and extent of morphologic innovation, taxonomic diversification and phylogenetic breadth, suggesting that several patterns, and likely processes, are involved. At least four distinct patterns of evolutionary radiation can be identified: novelty events, which generate new morphological complexity (altering the body plan of the group under consideration) but not necessarily with the associated production of many lower taxa; broad diversification events involving many independent lineages that undergo diversification, generate many new species and are driven by new ecological opportunities; economic radiations of a limited group of ecologically (but not necessarily phylogenetically) related clades exploiting a limited new ecologic opportunity; and adaptive radiations that may occur at any taxonomic level, but involve a rapid increase in diversity within a single clade, including “true”; adaptive radiations. Many events produce simple diversity increases with no corresponding increase in genetic/developmental/morphological/behavioral sophistication, but the most evolutionarily interesting events add new levels of complexity.     

Bats, clocks, and rocks: diversification patterns in Chiroptera

Identifying nonrandom clade diversification is a critical first step toward understanding the evolutionary processes underlying any radiation and how best to preserve future phylogenetic diversity. However, differences in diversification rates have not been quantitatively assessed for the majority of groups because of the lack of necessary analytical tools (e.g., complete species-level phylogenies, estimates of divergence times, and robust statistics which incorporate phylogenetic uncertainty and test appropriate null models of clade growth). Here, for the first time, we investigate diversification rate heterogeneity in one of the largest groups studied thus far, the bats (Mammalia: Chiroptera). We use a recent, robust statistical approach (whole-tree likelihood-based relative rate tests) on complete dated species-level supertree phylogenies. As has been demonstrated previously for most other groups, among-lineage diversification rate within bats has not been constant. However, we show that bat diversification is more heterogeneous than in other mammalian clades thus far studied. The whole-tree likelihood-based relative rates tests suggest that clades within the families Phyllostomidae and Molossidae underwent a number of significant changes in relative diversification rate. There is also some evidence for rate shifts within Pteropodidae, Emballonuridae, Rhinolophidae, Hipposideridae, and Vespertilionidae, but the significance of these shifts depends on polytomy resolution within each family. Diversification rate in bats has also not been constant, with the largest diversification rate shifts occurring 30-50 million years ago, a time overlapping with the greatest number of shifts in flowering plant diversification rates.     

Divergent patterns of diversification in courtship and genitalic characters of Timema walking‐sticks

Understanding the patterns of diversification in sexual traits and the selection underlying such diversification represents a major unresolved question in evolutionary biology. We examined the phylogenetic diversification for courtship and external genitalic characters across ten species of Timema walking‐sticks, to infer the tempos and modes of character change in these sexual traits and to draw inferences regarding the selective pressures underlying speciation and diversification in this clade. Rates of inferred change in male courtship behaviours were proportional to speciation events, but male external genitalic structures showed a pattern of continuous change across evolutionary time, with divergence proportional to branch lengths. These findings suggest that diversification of courtship behaviour is mediated by processes that occur in association with speciation, whereas diversification of genitalia occurs more or less continuously, most likely driven by forces of sexual selection.     

Trait evolution rates shape continental patterns of species richness in North America's most diverse angiosperm genus (Carex,Cyperaceae)

Ecological opportunity has been associated with increases in diversification rates across the tree of life. Under an ecological diversification model, the emergence of novel environments is hypothesized to promote morpho- and ecospace evolution. Whether this model holds at the clade level within the most species-rich angiosperm genus found in North America (Carex, Cyperaceae) is yet to be tested. Recent works demonstrate a temporal coupling of climate cooling and widespread colonization of Carex in North America, implicating ecological diversification. In addition, research has consistently found asymmetric patterns of lineage-level diversification in the genus. Why does variation in clade sizes exist in the genus? Is ecological diversification involved? In this study, we tested whether rates of morphological and ecological trait evolution are correlated with clade-level species richness in Carex of North America north of Mexico. We constructed a phylogeny of 477 species—an almost complete regional sample. We estimated rates of evolution of morphological traits, habitat, and climatic niche and assessed whether differences in rates of evolution correlate with species richness differences in replicate non-nested sister clades. Our work demonstrates significant positive correlations between climatic niche rates, habitat and reproductive morphological evolution, and species richness. This coupling of trait and niche evolution and species richness in a diverse, continental clade sample strongly suggests that the ability of clades to explore niche and functional space has shaped disparities in richness and functional diversity across the North American flora region. Our findings highlight the importance of the evolutionary history of trait and niche evolution in shaping continental and regional floras.