EXPLOSIVE RADIATION OR CRYPTIC MASS EXTINCTION? INTERPRETING SIGNATURES IN MOLECULAR PHYLOGENIES

The Major Features of Macroevolution

Evolutionary dynamics operating across deep time leave footprints in the shapes of phylogenetic trees. For the last several decades, researchers have used increasingly large and robust phylogenies to study the evolutionary history of individual clades and to investigate the causes of the glaring disparities in diversity among groups. Whereas typically not the focal point of individual clade-level studies, many researchers have remarked on recurrent patterns that have been observed across many different groups and at many different time scales. Whereas previous studies have documented various such regularities in topology and branch length distributions, they have typically focused on a single pattern and used a disparate collection (oftentimes, of quite variable reliability) of trees to assess it. Here we take advantage of modern megaphylogenies and unify previous disparate observations about the shapes embedded in the Tree of Life to create a catalog of the "major features of macroevolution." By characterizing such a large swath of subtrees in a consistent way, we hope to provide a set of phenomena that process-based macroevolutionary models of diversification ought to seek to explain.

Interplay of Ecological Opportunities and Functional Traits Drives the Evolution and Diversification of Millettiod Legumes (Fabaceae)

Understanding the striking diversity of the angiosperms is a paramount issue in biology and of interest to biologists. The Millettiod legumes is one of the most hyper-diverse groups of the legume family, containing many economically important medicine, furniture and craft species. In the present study, we explore how the interplay of past climate change, ecological opportunities and functional traits' evolution may have triggered diversification of the Millettiod legumes. Using a comprehensive species-level phylogeny from three plastid markers, we estimate divergence times, infer habit shifts, test the phylogenetic and temporal diversification heterogeneity, and reconstruct ancestral biogeographical ranges. We found that three dramatic accumulations of the Millettiod legumes occurred during the Miocene. The rapid diversification of the Millettiod legumes in the Miocene was driven by ecological opportunities created by the emergence of new niches and range expansion. Additionally, habit shifts and the switch between biomes might have facilitated the rapid diversification of the Millettiod legumes. The Millettiod legumes provide an excellent case for supporting the idea that the interplay of functional traits, biomes, and climatic and geographic factors drives evolutionary success.

Phylogenomics and continental biogeographic disjunctions: insight from the Australian starflowers (Calytrix)

PREMISE

Continental-scale disjunctions and associated drivers are core research interests in biogeographic studies. Here, we selected a species-rich Australian plant genus (Calytrix; Myrtaceae) as a case study to investigate these patterns. Species of this endemic Australian starflower genus have a disjunct distribution across the mesic fringes of the continent and are largely absent from the arid center.

METHODS

We used high-throughput sequencing to generate unprecedented resolution and near complete species-level nuclear and plastid phylogenies for Calytrix. BioGeoBEARS and biogeographic stochastic mapping were used to infer ancestral areas, the relative contributions of vicariance and dispersal events, and directionality of dispersal.

RESULTS

Present-day disjunctions in Calytrix are explained by a combination of scenarios: (1) retreat of multiple lineages from the continental center to the more mesic fringes as Australia became progressively more arid, with subsequent extinction in the center as well as (2) origination of ancestral lineages in southwestern Australia (SWA) for species-rich clades. The SWA biodiversity hotspot is a major diversification center and the most common source area of dispersals, with multiple lineages originating in SWA and subsequently spreading to the adjacent arid Eremaean region.

CONCLUSIONS

Our results suggest that major extinction, as a result of cooling and drying of the Australian continent in the Eocene-Miocene, shaped the present-day biogeography of Calytrix. We hypothesize that this peripheral vicariance pattern, which is similar to the African Rand flora, may explain the disjunctions of many other Australian plant groups. Further studies with densely sampled phylogenies are required to test this hypothesis.

Food Plant Shifts Drive the Diversification of Sack-Bearer Moths

AbstractLepidoptera are a highly diverse group of herbivorous insects; however, some superfamilies have relatively few species. Two alternative hypotheses for drivers of Lepidoptera diversity are shifts in food plant use or shifts from concealed to external feeding as larvae. Many studies address the former hypothesis but with bias toward externally feeding taxa. One of the most striking examples of species disparity between sister lineages in Lepidoptera is between the concealed-feeding sack-bearer moths (Mimallonoidea), which contain about 300 species, and externally feeding Macroheterocera, which have over 74,000 species. We provide the first dated tree of Mimallonidae to understand the diversification dynamics of these moths in order to fill a knowledge gap pertaining to drivers of diversity within an important concealed-feeding clade. We find that Mimallonidae is an ancient Lepidoptera lineage that originated in the Cretaceous ∼105 million years ago and has had a close association with the plant order Myrtales for the past 40 million years. Diversification dynamics are tightly linked with food plant usage in this group. Reliance on Myrtales may have influenced diversification of Mimallonidae because clades that shifted away from the ancestral condition of feeding on Myrtales have the highest speciation rates in the family.

Pulled Diversification Rates, Lineages-Through-Time Plots and Modern Macroevolutionary Modelling

Estimating time-dependent rates of speciation and extinction from dated phylogenetic trees of extant species (timetrees), and determining how and why they vary, is key to understanding how ecological and evolutionary processes shape biodiversity. Due to an increasing availability of phylogenetic trees, a growing number of process-based methods relying on the birth–death model have been developed in the last decade to address a variety of questions in macroevolution. However, this methodological progress has regularly been criticized such that one may wonder how reliable the estimations of speciation and extinction rates are. In particular, using lineages-through-time (LTT) plots, a recent study has shown that there are an infinite number of equally likely diversification scenarios that can generate any timetree. This has led to questioning whether or not diversification rates should be estimated at all. Here, we summarize, clarify, and highlight technical considerations on recent findings regarding the capacity of models to disentangle diversification histories. Using simulations, we illustrate the characteristics of newly proposed “pulled rates” and their utility. We recognize that the recent findings are a step forward in understanding the behavior of macroevolutionary modeling, but they in no way suggest we should abandon diversification modeling altogether. On the contrary, the study of macroevolution using phylogenetic trees has never been more exciting and promising than today. We still face important limitations in regard to data availability and methods, but by acknowledging them we can better target our joint efforts as a scientific community. [Birth–death models; extinction; phylogenetics; speciation.]

Conserved ancestral tropical niche but different continental histories explain the latitudinal diversity gradient in brush-footed butterflies

The global increase in species richness toward the tropics across continents and taxonomic groups, referred to as the latitudinal diversity gradient, stimulated the formulation of many hypotheses to explain the underlying mechanisms of this pattern. We evaluate several of these hypotheses to explain spatial diversity patterns in a butterfly family, the Nymphalidae, by assessing the contributions of speciation, extinction, and dispersal, and also the extent to which these processes differ among regions at the same latitude. We generate a time-calibrated phylogeny containing 2,866 nymphalid species (~45% of extant diversity). Neither speciation nor extinction rate variations consistently explain the latitudinal diversity gradient among regions because temporal diversification dynamics differ greatly across longitude. The Neotropical diversity results from low extinction rates, not high speciation rates, and biotic interchanges with other regions are rare. Southeast Asia is also characterized by a low speciation rate but, unlike the Neotropics, is the main source of dispersal events through time. Our results suggest that global climate change throughout the Cenozoic, combined with tropical niche conservatism, played a major role in generating the modern latitudinal diversity gradient of nymphalid butterflies.

The explosive radiation, intense host-shifts and long-term failure to speciate in the evolutionary history of the feather mite genus Analges (Acariformes: Analgidae) from European passerines

Mites of the genus Analges (Acariformes: Analgidae) inhabit the down feathers of passeriform birds. The evolutionary history of Analges and the co-phylogentic relationships between these mites and their hosts are unknown. Our phylogenetic analysis supported the monophyly of the genus, but it did not support previous taxonomic hypotheses subdividing the genus into the subgenera Analges and Analgopsis or arranging some species into the A. chelopus and A. passerinus species groups. Molecular data reveal seven new species inhabiting Eurasian passerines and support the existence of several multi-host species. According to molecular dating, the origin of the Analges (c. 41 Mya) coincided with the Eocene diversification of Passerida into Sylvioidea and Muscicapoidea–Passeroidea. The initial diversification of Analges took place on the Muscicapoidea clade, while remaining passerine superfamilies appear to have been colonized because of host-switching. Co-speciation appears to be relatively common among Analges species and their hosts, but the most striking pattern in the co-phylogenetic scenario involves numerous complete host-switches, spreads and several failures to speciate. The mechanism of long-term gene-flow among different populations of multi-host Analges species is enigmatic and difficult to resolve. Probably, in some cases mites could be transferred between birds via feathers used as nest material.

Understanding Diversity and Systematics in Australian Fabaceae Tribe Mirbelieae

Australia has a very diverse pea-flowered legume flora with 1715 native and naturalised species currently recognised. Tribe Mirbelieae s.l. includes 44% of Australia’s peas in 24 genera with 756 recognised species. However, several genera within the Pultenaea alliance in tribe Mirbelieae are considered to be non-monophyletic and two main options have been proposed: option one is to merge ca. 18 genera containing ca. 540 species (the largest genus, Pultenaea has nomenclatural priority); and option two is to re-circumscribe some genera and describe new genera as required to form monophyletic groups. At the species level, option one would require 76% of names to be changed; whereas based on available data, option two is likely to require, at most, 8.3% of names to change. Option two therefore provides the least nomenclatural disruption but cannot be implemented without a robust phylogenetic framework to define new generic limits. Here we present novel analyses of available plastid DNA data (trnL-F) which suggest that option two would be feasible once sufficient data are generated to resolve relationships. However, the reticulate evolutionary histories or past rapid speciation suggested for this group may prevent the resolution of all nodes. We propose targeted use of Next-Generation Sequencing technology as the best way to resolve relationships between the key clades in the tribe and present a framework for such a study. An overview of current taxonomy in the tribe is presented, along with the state of taxonomic knowledge and availability of published descriptions for electronic flora treatments. Several new combinations and typifications are published in an appendix.

Adaptive radiation and speciation in Rhipicephalus ticks: A medley of novel hosts, nested predator-prey food webs, off-host periods and dispersal along temperature variation gradients

Rhipicephalus are a species-diverse genus of ticks, mainly distributed in the Afrotropics with some species in the Palearctic and Oriental regions. Current taxonomic consensus comprise nine informal species groups/lineages based on immature morphology. This work integrates biogeographic, ecological and molecular lines of evidence to better understand Rhipicephalus evolution. Phylogenetic analysis based on four genes (12S, 16S, 28S-D2 and COI) recovered five distinct clades with nine descendant clades that are generally congruent with current taxonomy, with some exceptions. Historical biogeography is inferred from molecular divergence times, ancestral distribution areas, host-use and climate niches of four phylogenetically significant bioclimatic variables (isothermality, annual, seasonal and diurnal temperature range). Novel hosts enabled host-linked dispersal events into new environments, and ticks exploited new hosts through nested predator-prey connections in food webs. Diversification was further induced by climate niche partitioning along gradients in temperature range during off-host periods. Ancestral climate niche estimates corroborated dispersal events by indicating hypothetical ancestors moved into environments with different annual and seasonal temperature ranges along latitudinal gradients. Host size for immature and adult life stages was important for dispersal and subsequent diversification rates. Clades that utilise large, mobile hosts (ungulates and carnivores) early in development have wider geographic ranges but slower diversification rates, and those utilising small, less mobile hosts (rodents, lagomorphs and afroinsectivores) early in development have smaller ranges but higher diversification rates. These findings suggest diversification is driven by a complex set of factors linked to both host-associations (host size, ranges and mobility) and climate niche partitioning along annual and seasonal temperature range gradients that vary with latitude. Moreover, competitive interactions can reinforce these processes and drive speciation. Off-host periods facilitate adaptive radiation by enabling host switches along nested predator-prey connections in food webs, but at the cost of environmental exposure that partitions niches among dispersing progenitors, disrupting geneflow and driving diversification. As such, the evolution and ecological niches of Rhipicephalus are characterised by trade-offs between on- and off-host periods, and these trade-offs interact with nested predator-prey connections in food webs, host-use at different life stages, as well as gradients in annual and seasonal temperature ranges to drive adaptive radiation and speciation.

Tectonics, climate and the diversification of the tropical African terrestrial flora and fauna

Tropical Africa is home to an astonishing biodiversity occurring in a variety of ecosystems. Past climatic change and geological events have impacted the evolution and diversification of this biodiversity. During the last two decades, around 90 dated molecular phylogenies of different clades across animals and plants have been published leading to an increased understanding of the diversification and speciation processes generating tropical African biodiversity. In parallel, extended geological and palaeoclimatic records together with detailed numerical simulations have refined our understanding of past geological and climatic changes in Africa. To date, these important advances have not been reviewed within a common framework. Here, we critically review and synthesize African climate, tectonics and terrestrial biodiversity evolution throughout the Cenozoic to the mid-Pleistocene, drawing on recent advances in Earth and life sciences. We first review six major geo-climatic periods defining tropical African biodiversity diversification by synthesizing 89 dated molecular phylogeny studies. Two major geo-climatic factors impacting the diversification of the sub-Saharan biota are highlighted. First, Africa underwent numerous climatic fluctuations at ancient and more recent timescales, with tectonic, greenhouse gas, and orbital forcing stimulating diversification. Second, increased aridification since the Late Eocene led to important extinction events, but also provided unique diversification opportunities shaping the current tropical African biodiversity landscape. We then review diversification studies of tropical terrestrial animal and plant clades and discuss three major models of speciation: (i) geographic speciation via vicariance (allopatry); (ii) ecological speciation impacted by climate and geological changes, and (iii) genomic speciation via genome duplication. Geographic speciation has been the most widely documented to date and is a common speciation model across tropical Africa. We conclude with four important challenges faced by tropical African biodiversity research: (i) to increase knowledge by gathering basic and fundamental biodiversity information; (ii) to improve modelling of African geophysical evolution throughout the Cenozoic via better constraints and downscaling approaches; (iii) to increase the precision of phylogenetic reconstruction and molecular dating of tropical African clades by using next generation sequencing approaches together with better fossil calibrations; (iv) finally, as done here, to integrate data better from Earth and life sciences by focusing on the interdisciplinary study of the evolution of tropical African biodiversity in a wider geodiversity context.

Phylogenomic Resolution of Sea Spider Diversification through Integration of Multiple Data Classes

Despite significant advances in invertebrate phylogenomics over the past decade, the higher-level phylogeny of Pycnogonida (sea spiders) remains elusive. Due to the inaccessibility of some small-bodied lineages, few phylogenetic studies have sampled all sea spider families. Previous efforts based on a handful of genes have yielded unstable tree topologies. Here, we inferred the relationships of 89 sea spider species using targeted capture of the mitochondrial genome, 56 conserved exons, 101 ultraconserved elements, and 3 nuclear ribosomal genes. We inferred molecular divergence times by integrating morphological data for fossil species to calibrate 15 nodes in the arthropod tree of life. This integration of data classes resolved the basal topology of sea spiders with high support. The enigmatic family Austrodecidae was resolved as the sister group to the remaining Pycnogonida and the small-bodied family Rhynchothoracidae as the sister group of the robust-bodied family Pycnogonidae. Molecular divergence time estimation recovered a basal divergence of crown group sea spiders in the Ordovician. Comparison of diversification dynamics with other marine invertebrate taxa that originated in the Paleozoic suggests that sea spiders and some crustacean groups exhibit resilience to mass extinction episodes, relative to mollusk and echinoderm lineages.

Genome-wide macroevolutionary signatures of key innovations in butterflies colonizing new host plants

The mega-diversity of herbivorous insects is attributed to their co-evolutionary associations with plants. Despite abundant studies on insect-plant interactions, we do not know whether host-plant shifts have impacted both genomic adaptation and species diversification over geological times. We show that the antagonistic insect-plant interaction between swallowtail butterflies and the highly toxic birthworts began 55 million years ago in Beringia, followed by several major ancient host-plant shifts. This evolutionary framework provides a valuable opportunity for repeated tests of genomic signatures of macroevolutionary changes and estimation of diversification rates across their phylogeny. We find that host-plant shifts in butterflies are associated with both genome-wide adaptive molecular evolution (more genes under positive selection) and repeated bursts of speciation rates, contributing to an increase in global diversification through time. Our study links ecological changes, genome-wide adaptations and macroevolutionary consequences, lending support to the importance of ecological interactions as evolutionary drivers over long time periods.

Species Selection Regime and Phylogenetic Tree Shape

Species selection, the effect of heritable traits in generating between-lineage diversification rate differences, provides a valuable conceptual framework for understanding the relationship between traits, diversification, and phylogenetic tree shape. An important challenge, however, is that the nature of real diversification landscapes—curves or surfaces which describe the propensity of species-level lineages to diversify as a function of one or more traits—remains poorly understood. Here, we present a novel, time-stratified extension of the QuaSSE model in which speciation/extinction rate is specified as a static or temporally shifting Gaussian or skewed-Gaussian function of the diversification trait. We then use simulations to show that the generally imbalanced nature of real phylogenetic trees, as well as their generally greater than expected frequency of deep branching events, are typical outcomes when diversification is treated as a dynamic, trait-dependent process. Focusing on four basic models (Gaussian-speciation with and without background extinction; skewed-speciation; Gaussian-extinction), we also show that particular features of the species selection regime produce distinct tree shape signatures and that, consequently, a combination of tree shape metrics has the potential to reveal the species selection regime under which a particular lineage diversified. We evaluate this idea empirically by comparing the phylogenetic trees of plant lineages diversifying within climatically and geologically stable environments of the Greater Cape Floristic Region, with those of lineages diversifying in environments that have experienced major change through the Late Miocene-Pliocene. Consistent with our expectations, the trees of lineages diversifying in a dynamic context are less balanced, show a greater concentration of branching events close to the present, and display stronger diversification rate-trait correlations. We suggest that species selection plays an important role in shaping phylogenetic trees but recognize the need for an explicit probabilistic framework within which to assess the likelihoods of alternative diversification scenarios as explanations of a particular tree shape. [Cape flora; diversification landscape; environmental change; gamma statistic; species selection; time-stratified QuaSSE model; trait-dependent diversification; tree imbalance.]

A fossil-calibrated phylogeny reveals the biogeographic history of the Cladrastis clade, an amphi-Pacific early-branching group in papilionoid legumes

The early-branching Cladrastis clade of papilionoid legumes (Leguminosae, Papilionoideae) has an intriguing amphi-Pacific disjunct distribution in eastern Asia and temperate-tropical Americas. Here we used nuclear and three plastid regions to reconstruct the phylogenetic relationships and divergence times in the Cladrastis clade, as well as the evolution of morphological characters that might have been key in its biogeographic history. The ancestral character state estimation revealed that the most recent common ancestor of the Cladrastis clade was deciduous trees possessing compressed, winged fruits. The Cladrastis clade was inferred to have originated in the mid-latitude thermophilic forests of North America in the early Eocene, followed by the split between ancestors of wing-fruited Platyosprion and the non-wing-fruited group, and later the divergence of Cladrastis s.s. from the non-wing-fruited group in middle Eocene. Platyosprion and Cladrastis s.s. display an "out-of-North-America" biogeographic pattern and might have migrated to Asia via the Bering land bridge (BLB) or the North Atlantic land bridges (NALB) during middle to late Eocene. Our results, coupled with the relatively well documented fossil record for the clade, suggest that Platyosprion experienced an extinction event in North America caused by climatic cooling around the Eocene-Oligocene transition, which drove a major vegetation shift in western North America, in turn serving as a barrier for the vicariance of Pickeringia and Styphnolobium. The evolution of shrubby habit and sclerophyllous leaves in the former might be adaption to the chaparral vegetation in southwestern North America; the latter gained the trait of moniliform, succulent fruit. Styphnolobium further dispersed southward to tropical North America in the Oligocene, and eastward to Asia through BLB during middle Miocene. Subsequent sundering of BLB facilitated the vicariance of St. affine and St. japonicum.

The early wasp plucks the flower: disparate extant diversity of sawfly superfamilies (Hymenoptera: ‘Symphyta’) may reflect asynchronous switching to angiosperm hosts

The insect order Hymenoptera originated during the Permian nearly 300 Mya. Ancestrally herbivorous hymenopteran lineages today make up the paraphyletic suborder ‘Symphyta’, which encompasses c. 8200 species with very diverse host-plant associations. We use phylogeny-based statistical analyses to explore the drivers of diversity dynamics within the ‘Symphyta’, with a particular focus on the hypothesis that diversification of herbivorous insects has been driven by the explosive radiation of angiosperms during and after the Cretaceous. Our ancestral-state estimates reveal that the first symphytans fed on gymnosperms, and that shifts onto angiosperms and pteridophytes – and back – have occurred at different time intervals in different groups. Trait-dependent analyses indicate that average net diversification rates do not differ between symphytan lineages feeding on angiosperms, gymnosperms or pteridophytes, but trait-independent models show that the highest diversification rates are found in a few angiosperm-feeding lineages that may have been favoured by the radiations of their host taxa during the Cenozoic. Intriguingly, lineages-through-time plots show signs of an early Cretaceous mass extinction, with a recovery starting first in angiosperm-associated clades. Hence, the oft-invoked assumption of herbivore diversification driven by the rise of flowering plants may overlook a Cretaceous global turnover in insect herbivore communities during the rapid displacement of gymnosperm- and pteridophyte-dominated floras by angiosperms.

Exploring the effects of character construction and choice, outgroups and analytical method on phylogenetic inference from discrete characters in extant crocodilians

Phylogenies for fossil taxa must be inferred from morphology, but accuracy of inference is questionable. Here, morphological characters for extant crocodilians are investigated to assess how to improve inference accuracy. The homoplasy of characters is assessed against a DNA-based phylogenetic tree. Cranial characters are significantly less homoplastic, but this result is perhaps confounded by research effort. Meristic characters are significantly more homoplastic and should be used with caution. Characters were reassessed first hand and documented. Those characters passing tests of robust construction are significantly less homoplastic. Suggestions are made for means to improve coding of discrete characters. Phylogenies inferred using only robust characters and a reassessed matrix, including corrected scorings, were not overall closer to the DNA tree, but did often place the gharial (Gavialis) in a position agreeing with or closer to it. The effects of the choice of analytical method were modest, but Bayesian analysis of the reassessed matrix placed Gavialis and Mecistops (slender-snouted crocodile) in DNA-concordant positions. Use of extant rather than extinct outgroups, even with the original matrix, placed Gavialis in a more DNA-concordant position, as did factoring out 3D skull shape. The morphological case for placement of Gavialis outside other extant crocodilians is arguably overstated, with many characters linked to skull shape.

Turnover of southern cypresses in the post-Gondwanan world: extinction, transoceanic dispersal, adaptation and rediversification

Cupressaceae subfamily Callitroideae has been an important exemplar for vicariance biogeography, but its history is more than just disjunctions resulting from continental drift. We combine fossil and molecular data to better assess its extinction and, sometimes, rediversification after past global change. Key fossils were reassessed and their phylogenetic placement for calibration was determined using trait mapping and Bayes Factors. Five vicariance hypotheses were tested by comparing molecular divergence times with the timing of tectonic rifting. The role of adaptation to fire (serotiny) in its spread across a drying Australia was tested for Callitris. Our findings suggest that three transoceanic disjunctions within the Callitroideae probably arose from long-distance dispersal. A signature of extinction, centred on the end-Eocene global climatic chilling and drying, is evident in lineages-through-time plots and in the fossil record. Callitris, the most diverse extant callitroid genus, suffered extinctions but surviving lineages adapted and re-radiated into dry, fire-prone biomes that expanded in the Neogene. Serotiny, a key adaptation to fire, likely evolved in Callitris coincident with the biome shift. Both extinction and adaptive shifts have probably played major roles in this chronicle of turnover and renewal, but better understanding of biogeographical history requires improved taxonomy of fossils.

Analysis of temporal diversification of African Cyprinidae (Teleostei, Cypriniformes)

Recent evidence that freshwater fishes diversify faster than marine fishes signifies that the evolutionary history of biodiversity in freshwater system is of particular interest. Here, the evolutionary diversification events of African Cyprinidae, a freshwater fish family with wide geographic distribution, were reconstructed and analysed. The overall diversification rate of African Cyprinidae is 0.08 species per million year (when extinction rate is very high, i.e., ε = 0.9) and 0.11 species per million year (when ε = 0). This overall rate is lower than the rate reported for African Cichlids, suggesting that African freshwaters might be less conducive for a rapid diversification of Cyprinidae. However, the observed diversification events of African Cyprinidae occurred in the last 10 million years. The temporal pattern of these events follows a non-constant episodic birth-death model (Bayes Factor > 28) and the rate-constant model never outperformed any of the non-constant models tested. The fact that most diversification events occurred in the last 10 million years supports the pattern reported for Cyprinidae in other continent, e.g., Asia, perhaps pointing to concomitant diversification globally. However, the diversification events coincided with major geologic and paleo-climatic events in Africa, suggesting that geological and climatic events may have mediated the diversification patterns of Cyprinidae on the continent.

Plastome sequences and exploration of tree-space help to resolve the phylogeny of riceflowers (Thymelaeaceae: Pimelea)

Data sets comprising small numbers of genetic markers are not always able to resolve phylogenetic relationships. This has frequently been the case in molecular systematic studies of plants, with many analyses being based on sequence data from only two or three chloroplast genes. An example of this comes from the riceflowers Pimelea Banks & Sol. ex Gaertn. (Thymelaeaceae), a large genus of flowering plants predominantly distributed in Australia. Despite the considerable morphological variation in the genus, low sequence divergence in chloroplast markers has led to the phylogeny of Pimelea remaining largely uncertain. In this study, we resolve the backbone of the phylogeny of Pimelea in comprehensive Bayesian and maximum-likelihood analyses of plastome sequences from 41 taxa. However, some relationships received only moderate to poor support, and the Pimelea clade contained extremely short internal branches. By using topology-clustering analyses, we demonstrate that conflicting phylogenetic signals can be found across the trees estimated from individual chloroplast protein-coding genes. A relaxed-clock dating analysis reveals that Pimelea arose in the mid-Miocene, with most divergences within the genus occurring during a subsequent rapid diversification. Our new phylogenetic estimate offers better resolution and is more strongly supported than previous estimates, providing a platform for future taxonomic revisions of both Pimelea and the broader subfamily. Our study has demonstrated the substantial improvements in phylogenetic resolution that can be achieved using plastome-scale data sets in plant molecular systematics.

Evolutionary persistence in Gunnera and the contribution of southern plant groups to the tropical Andes biodiversity hotspot

Several studies have demonstrated the contribution of northern immigrants to the flora of the tropical Andes—the world’s richest and most diverse biodiversity hotspot. However, much less is known about the biogeographic history and diversification of Andean groups with southern origins, although it has been suggested that northern and southern groups have contributed roughly equally to the high Andean (i.e., páramo) flora. Here we infer the evolutionary history of the southern hemisphere plant genus Gunnera, a lineage with a rich fossil history and an important ecological role as an early colonising species characteristic of wet, montane environments. Our results show striking contrasts in species diversification, where some species may have persisted for some 90 million years, and whereas others date to less than 2 Ma since origination. The outstanding longevity of the group is likely linked to a high degree of niche conservatism across its highly disjunct range, whereby Gunnera tracks damp and boggy soils in cool habitats. Colonisation of the northern Andes is related to Quaternary climate change, with subsequent rapid diversification appearing to be driven by their ability to take advantage of environmental opportunities. This study demonstrates the composite origin of a mega-diverse biota.

Diversification in tropics and subtropics following the mid-Miocene climate change: A case study of the spider genus Nesticella

Caves may offer suitable refugia for troglophilic invertebrates during periods of unfavourable climatic conditions because of their stable microclimates. As a consequence, allopatric divergence from their epigean counterparts may occur, leading to formation of truly hypogean communities (the Climatic Relict Hypothesis). Unlike the well-studied effects of Pleistocene glaciations, we know little about how ancient climate changes drove the development of cave-dwelling organisms living at both middle and lower latitudes. We investigate the evolutionary history of the troglophilic spider genus Nesticella (Araneae, Nesticidae) in relation to Asian Neogene (23-2.6 Ma) climatic changes. Our analyses discern clear differences in the evolution of the two main clades of Nesticella, which occur in temperate/subtropical and tropical latitudes. Eastern Asian Nesticella gradually evolved greater sedentariness and a strict subterranean lifestyle starting from the middle Miocene Epoch (~15-14 Ma) in conjunction with the progressive deterioration of the climate and vegetational shifts. Caves appear to have acted as refugia because of their internally uniform temperature and humidity, which allowed these spiders to survive increasing external seasonality and habitat loss. In contrast, a uniform accumulation of lineages, long-lasting times for dispersals and the lack of a comparable habitat shifting characterized the tropical lineage. This difference in pattern likely owes to the mild effects of climate change at low latitudes and the consequent lack of strong climatic drivers in tropical environments. Thus, the mid-Miocene climatic shift appears to be the major evolutionary force shaping the ecological differences between Asian troglophilic invertebrates and the driver of the permanent hypogean communities in middle latitudes.

Paleoclimate determines diversification patterns in the fossorial snake family Uropeltidae Cuvier, 1829

Understanding how and why diversification rates vary across evolutionary time is central to understanding how biodiversity is generated and maintained. Recent mathematical models that allow estimation of diversification rates across time from reconstructed phylogenies have enabled us to make inferences on how biodiversity copes with environmental change. Here, we explore patterns of temporal diversification in Uropeltidae, a diverse fossorial snake family. We generate a time-calibrated phylogenetic hypothesis for Uropeltidae and show a significant correlation between diversification rate and paleotemperature during the Cenozoic. We show that the temporal diversification pattern of this group is punctuated by one rate shift event with a decrease in diversification and turnover rate between ca. 11Ma to present, but there is no strong support for mass extinction events. The analysis indicates higher turnover during periods of drastic climatic fluctuations and reduced diversification rates associated with contraction and fragmentation of forest habitats during the late Miocene. Our study highlights the influence of environmental fluctuations on diversification rates in fossorial taxa such as uropeltids, and raises conservation concerns related to present rate of climate change.

The tree balance signature of mass extinction is erased by continued evolution in clades of constrained size with trait-dependent speciation

The kind and duration of phylogenetic topological “signatures” left in the wake of macroevolutionary events remain poorly understood. To this end, we examined a broad range of simulated phylogenies generated using trait-biased, heritable speciation probabilities and mass extinction that could be either random or selective on trait value, but also using background extinction and diversity-dependence to constrain clade sizes. In keeping with prior results, random mass extinction increased imbalance of clades that recovered to pre-extinction size, but was a relatively weak effect. Mass extinction that was selective on trait values tended to produce clades of similar or greater balance compared to random extinction or controls. Allowing evolution to continue past the point of clade-size recovery resulted in erosion and eventual erasure of this signal, with all treatments converging on similar values of imbalance, except for very intense extinction regimes targeted at taxa with high speciation rates. Return to a more balanced state with extended post-extinction evolution was also associated with loss of the previous phylogenetic root in most treatments. These results further demonstrate that while a mass extinction event can produce a recognizable phylogenetic signal, its effects become increasingly obscured the further an evolving clade gets from that event, with any sharp imbalance due to unrelated evolutionary factors.

Phylogenetic and paleobotanical evidence for late Miocene diversification of the Tertiary subtropical lineage of ivies (Hedera L., Araliaceae)

BACKGROUND

Hedera (ivies) is one of the few temperate genera of the primarily tropical Asian Palmate group of the Araliaceae, which extends its range out of Asia to Europe and the Mediterranean basin. Phylogenetic and phylogeographic results suggested Asia as the center of origin and the western Mediterranean region as one of the secondary centers of diversification. The bird-dispersed fleshy fruits of ivies suggest frequent dispersal over long distances (e.g. Macaronesian archipelagos), although reducing the impact of geographic barriers to gene flow in mainland species. Genetic isolation associated with geographic barriers and independent polyploidization events have been postulated as the main driving forces of diversification. In this study we aim to evaluate past and present diversification patterns in Hedera within a geographic and temporal framework to clarify the biogeographic history of the genus.

RESULTS

Phylogenetic (biogeographic, time divergence and diversification) and phylogeographic (coalescence) analyses using four DNA regions (nrITS, trnH-psbA, trnT-trnL, rpl32) revealed a complex spatial pattern of lineage divergence. Scarce geographic limitation to gene flow and limited diversification are observed during the early-mid Miocene, followed by a diversification rate increase related to geographic divergence from the Tortonian/Messinian. Genetic and palaeobotanical evidence points the origin of the Hedera clade in Asia, followed by a gradual E-W Asian extinction and the progressive E-W Mediterranean colonization. The temporal framework for the E Asia - W Mediterranean westward colonization herein reported is congruent with the fossil record. Subsequent range expansion in Europe and back colonization to Asia is also inferred. Uneven diversification among geographic areas occurred from the Tortonian/Messinian onwards with limited diversification in the newly colonized European and Asian regions. Eastern and western Mediterranean regions acted as refugia for Miocene and post-Miocene lineages, with a similar role as consecutive centers of centrifugal dispersal (including islands) and speciation.

CONCLUSIONS

The Miocene Asian extinction and European survival of Hedera question the general pattern of Tertiary regional extinction of temperate angiosperms in Europe while they survived in Asia. The Tortonian/Messinian diversification increase of ivies in the Mediterranean challenges the idea that this aridity period was responsible for the extinction of the Mediterranean subtropical Tertiary flora. Differential responses of Hedera to geographic barriers throughout its evolutionary history, linked to spatial isolation related to historical geologic and climatic constraints may have shaped diversification of ivies in concert with recurrent polyploidy.

Phylogenetic perspectives on reef fish functional traits

Functional traits have been fundamental to the evolution and diversification of entire fish lineages on coral reefs. Yet their relationship with the processes promoting speciation, extinction and the filtering of local species pools remains unclear. We review the current literature exploring the evolution of diet, body size, water column use and geographic range size in reef-associated fishes. Using published and new data, we mapped functional traits on to published phylogenetic trees to uncover evolutionary patterns that have led to the current functional diversity of fishes on coral reefs. When examining reconstructed patterns for diet and feeding mode, we found examples of independent transitions to planktivory across different reef fish families. Such transitions and associated morphological alterations may represent cases in which ecological opportunity for the exploitation of different resources drives speciation and adaptation. In terms of body size, reconstructions showed that both large and small sizes appear multiple times within clades of mid-sized fishes and that extreme body sizes have arisen mostly in the last 10 million years (Myr). The reconstruction of range size revealed many cases of disparate range sizes among sister species. Such range size disparity highlights potential vicariant processes through isolation in peripheral locations. When accounting for peripheral speciation processes in sister pairs, we found a significant relationship between labrid range size and lineage age. The diversity and evolution of traits within lineages is influenced by trait-environment interactions as well as by species and trait-trait interactions, where the presence of a given trait may trigger the development of related traits or behaviours. Our effort to assess the evolution of functional diversity across reef fish clades adds to the burgeoning research focusing on the evolutionary and ecological roles of functional traits. We argue that the combination of a phylogenetic and a functional approach will improve the understanding of the mechanisms of species assembly in extraordinarily rich coral reef communities.

Phylogeny and biogeography of the amphi-Pacific genus Aphananthe

Aphananthe is a small genus of five species showing an intriguing amphi-Pacific distribution in eastern, southern and southeastern Asia, Australia, and Mexico, also with one species in Madagascar. The phylogenetic relationships of Aphananthe were reconstructed with two nuclear (ITS & ETS) and two plastid (psbA-trnH & trnL-trnF) regions. Clade divergence times were estimated with a Bayesian approach, and the ancestral areas were inferred using the dispersal-extinction-cladogenesis and Bayesian Binary MCMC analyses. Aphananthe was supported to be monophyletic, with the eastern Asian A. aspera resolved as sister to a clade of the remaining four species. Aphananthe was inferred to have originated in the Late Cretaceous (71.5 mya, with 95% HPD: 66.6-81.3 mya), and the crown age of the genus was dated to be in the early Miocene (19.1 mya, with 95% HPD: 12.4-28.9 mya). The fossil record indicates that Aphananthe was present in the high latitude thermophilic forests in the early Tertiary, and experienced extinctions from the middle Tertiary onwards. Aphananthe originated in Europe based on the inference that included fossil and extant species, but eastern Asia was estimated to be the ancestral area of the clade of the extant species of Aphananthe. Both the West Gondwanan vicariance hypothesis and the boreotropics hypothesis could be excluded as explanation for its amphi-Pacific distribution. Long-distance dispersals out of eastern Asia into North America, southern and southeastern Asia and Australia, and Madagascar during the Miocene account for its wide intercontinental disjunct distribution.

The Evolutionary Legacy of Diversification Predicts Ecosystem Function

Theory suggests that the structure of evolutionary history represented in a species community may affect its functioning, but phylogenetic diversity metrics do not allow for the identification of major differences in this structure. Here we propose a new metric, ELDERness (for Evolutionary Legacy of DivERsity) to estimate evolutionary branching patterns within communities by fitting a polynomial function to lineage-through-time (LTT) plots. We illustrate how real and simulated community branching patterns can be more correctly described by ELDERness and can successfully predict ecosystem functioning. In particular, the evolutionary history of branching patterns can be encapsulated by the parameters of third-order polynomial functions and further measured through only two parameters, the "ELDERness surfaces." These parameters captured variation in productivity of a grassland community better than existing phylogenetic diversity or diversification metrics and independent of species richness or presence of nitrogen fixers. Specifically, communities with small ELDERness surfaces (constant accumulation of lineages through time in LTT plots) were more productive, consistent with increased productivity resulting from complementary lineages combined with niche filling within lineages. Overall, while existing phylogenetic diversity metrics remain useful in many contexts, we suggest that our ELDERness approach better enables testing hypotheses that relate complex patterns of macroevolutionary history represented in local communities to ecosystem functioning.

If Dung Beetles (Scarabaeidae: Scarabaeinae) Arose in Association with Dinosaurs, Did They Also Suffer a Mass Co-Extinction at the K-Pg Boundary?

The evolutionary success of beetles and numerous other terrestrial insects is generally attributed to co-radiation with flowering plants but most studies have focused on herbivorous or pollinating insects. Non-herbivores represent a significant proportion of beetle diversity yet potential factors that influence their diversification have been largely unexamined. In the present study, we examine the factors driving diversification within the Scarabaeidae, a speciose beetle family with a range of both herbivorous and non-herbivorous ecologies. In particular, it has been long debated whether the key event in the evolution of dung beetles (Scarabaeidae: Scarabaeinae) was an adaptation to feeding on dinosaur or mammalian dung. Here we present molecular evidence to show that the origin of dung beetles occurred in the middle of the Cretaceous, likely in association with dinosaur dung, but more surprisingly the timing is consistent with the rise of the angiosperms. We hypothesize that the switch in dinosaur diet to incorporate more nutritious and less fibrous angiosperm foliage provided a palatable dung source that ultimately created a new niche for diversification. Given the well-accepted mass extinction of non-avian dinosaurs at the Cretaceous-Paleogene boundary, we examine a potential co-extinction of dung beetles due to the loss of an important evolutionary resource, i.e., dinosaur dung. The biogeography of dung beetles is also examined to explore the previously proposed “out of Africa” hypothesis. Given the inferred age of Scarabaeinae as originating in the Lower Cretaceous, the major radiation of dung feeders prior to the Cenomanian, and the early divergence of both African and Gondwanan lineages, we hypothesise that that faunal exchange between Africa and Gondwanaland occurred during the earliest evolution of the Scarabaeinae. Therefore we propose that both Gondwanan vicariance and dispersal of African lineages is responsible for present day distribution of scarabaeine dung beetles and provide examples.

Extinction in Phylogenetics and Biogeography: From Timetrees to Patterns of Biotic Assemblage

Global climate change and its impact on biodiversity levels have made extinction a relevant topic in biological research. Yet, until recently, extinction has received less attention in macroevolutionary studies than speciation; the reason is the difficulty to infer an event that actually eliminates rather than creates new taxa. For example, in biogeography, extinction has often been seen as noise, introducing homoplasy in biogeographic relationships, rather than a pattern-generating process. The molecular revolution and the possibility to integrate time into phylogenetic reconstructions have allowed studying extinction under different perspectives. Here, we review phylogenetic (temporal) and biogeographic (spatial) approaches to the inference of extinction and the challenges this process poses for reconstructing evolutionary history. Specifically, we focus on the problem of discriminating between alternative high extinction scenarios using time trees with only extant taxa, and on the confounding effect introduced by asymmetric spatial extinction – different rates of extinction across areas – in biogeographic inference. Finally, we identify the most promising avenues of research in both fields, which include the integration of additional sources of evidence such as the fossil record or environmental information in birth–death models and biogeographic reconstructions, the development of new models that tie extinction rates to phenotypic or environmental variation, or the implementation within a Bayesian framework of parametric non-stationary biogeographic models.

A new time tree reveals Earth history's imprint on the evolution of modern birds

Determining the timing of diversification of modern birds has been difficult. We combined DNA sequences of clock-like genes for most avian families with 130 fossil birds to generate a new time tree for Neornithes and investigated their biogeographic and diversification dynamics. We found that the most recent common ancestor of modern birds inhabited South America around 95 million years ago, but it was not until the Cretaceous-Paleogene transition (66 million years ago) that Neornithes began to diversify rapidly around the world. Birds used two main dispersion routes: reaching the Old World through North America, and reaching Australia and Zealandia through Antarctica. Net diversification rates increased during periods of global cooling, suggesting that fragmentation of tropical biomes stimulated speciation. Thus, we found pervasive evidence that avian evolution has been influenced by plate tectonics and environmental change, two basic features of Earth's dynamics.

Inflation of Molecular Clock Rates and Dates: Molecular Phylogenetics, Biogeography, and Diversification of a Global Cicada Radiation from Australasia (Hemiptera: Cicadidae: Cicadettini)

Dated phylogenetic trees are important for studying mechanisms of diversification, and molecular clocks are important tools for studies of organisms lacking good fossil records. However, studies have begun to identify problems in molecular clock dates caused by uncertainty of the modeled molecular substitution process. Here we explore Bayesian relaxed-clock molecular dating while studying the biogeography of ca. 200 species from the global cicada tribe Cicadettini. Because the available fossils are few and uninformative, we calibrate our trees in part with a cytochrome oxidase I (COI) clock prior encompassing a range of literature estimates for arthropods. We show that tribe-level analyses calibrated solely with the COI clock recover extremely old dates that conflict with published estimates for two well-studied New Zealand subclades within Cicadettini. Additional subclade analyses suggest that COI relaxed-clock rates and maximum-likelihood branch lengths become inflated relative to EF-1[Formula: see text] intron and exon rates and branch lengths as clade age increases. We present corrected estimates derived from: (i) an extrapolated EF-1[Formula: see text] exon clock derived from COI-calibrated analysis within the largest New Zealand subclade; (ii) post hoc scaling of the tribe-level chronogram using results from subclade analyses; and (iii) exploitation of a geological calibration point associated with New Caledonia. We caution that considerable uncertainty is generated due to dependence of substitution estimates on both the taxon sample and the choice of model, including gamma category number and the choice of empirical versus estimated base frequencies. Our results suggest that diversification of the tribe Cicadettini commenced in the early- to mid-Cenozoic and continued with the development of open, arid habitats in Australia and worldwide. We find that Cicadettini is a rare example of a global terrestrial animal group with an Australasian origin, with all non-Australasian genera belonging to two distal clades. Within Australia, we show that Cicadettini is more widely distributed than any other cicada tribe, diverse in temperate, arid and monsoonal habitats, and nearly absent from rainforests. We comment on the taxonomic implications of our findings for thirteen cicada genera.

Detecting patterns of species diversification in the presence of both rate shifts and mass extinctions

BACKGROUND

Recent methodological advances allow better examination of speciation and extinction processes and patterns. A major open question is the origin of large discrepancies in species number between groups of the same age. Existing frameworks to model this diversity either focus on changes between lineages, neglecting global effects such as mass extinctions, or focus on changes over time which would affect all lineages. Yet it seems probable that both lineages differences and mass extinctions affect the same groups.

RESULTS

Here we used simulations to test the performance of two widely used methods under complex scenarios of diversification. We report good performances, although with a tendency to over-predict events with increasing complexity of the scenario.

CONCLUSION

Overall, we find that lineage shifts are better detected than mass extinctions. This work has significance to assess the methods currently used to estimate changes in diversification using phylogenetic trees. Our results also point toward the need to develop new models of diversification to expand our capabilities to analyse realistic and complex evolutionary scenarios.

Three explanations for biodiversity hotspots: small range size, geographical overlap and time for species accumulation. An Australian case study

To understand the generation and maintenance of biodiversity hotspots, we tested three major hypotheses: rates of diversification, ecological limits to diversity, and time for species accumulation. Using dated molecular phylogenies, measures of species' range size and geographical clade overlap, niche modelling, and lineages-through-time plots of Australian Fabaceae, we compared the southwest Australia Floristic Region (SWAFR; a global biodiversity hotspot) with a latitudinally equivalent non-hotspot, southeast Australia (SEA). Ranges of species (real and simulated) were smaller in the SWAFR than in SEA. Geographical overlap of clades was significantly greater for Daviesia in the SWAFR than in SEA, but the inverse for Bossiaea. Lineage diversification rates over the past 10 Myr did not differ between the SWAFR and SEA in either genus. Interaction of multiple factors probably explains the differences in measured diversity between the two regions. Steeper climatic gradients in the SWAFR probably explain the smaller geographical ranges of both genera there. Greater geographical overlap of clades in the SWAFR, combined with a longer time in the region, can explain why Daviesia is far more species-rich there than in SEA. Our results indicate that the time for speciation and ecological limits hypotheses, in concert, can explain the differences in biodiversity.

Confluence, synnovation, and depauperons in plant diversification

We review the empirical phylogenetic literature on plant diversification, highlighting challenges in separating the effects of speciation and extinction, in specifying diversification mechanisms, and in making convincing arguments. In recent discussions of context dependence, key opportunities and landscapes, and indirect effects and lag times, we see a distinct shift away from single-point/single-cause 'key innovation' hypotheses toward more nuanced explanations involving multiple interacting causal agents assembled step-wise through a tree. To help crystalize this emerging perspective we introduce the term 'synnovation' (a hybrid of 'synergy' and 'innovation') for an interacting combination of traits with a particular consequence ('key synnovation' in the case of increased diversification rate), and the term 'confluence' for the sequential coming together of a set of traits (innovations and synnovations), environmental changes, and geographic movements along the branches of a phylogenetic tree. We illustrate these concepts using the radiation of Bromeliaceae. We also highlight the generality of these ideas by considering how rate heterogeneity associated with a confluence relates to the existence of particularly species-poor lineages, or 'depauperons.' Many challenges are posed by this re-purposed research framework, including difficulties associated with partial taxon sampling, uncertainty in divergence time estimation, and extinction.

Origin and diversification of living cycads: a cautionary tale on the impact of the branching process prior in Bayesian molecular dating

Background

Bayesian relaxed-clock dating has significantly influenced our understanding of the timeline of biotic evolution. This approach requires the use of priors on the branching process, yet little is known about their impact on divergence time estimates. We investigated the effect of branching priors using the iconic cycads. We conducted phylogenetic estimations for 237 cycad species using three genes and two calibration strategies incorporating up to six fossil constraints to (i) test the impact of two different branching process priors on age estimates, (ii) assess which branching prior better fits the data, (iii) investigate branching prior impacts on diversification analyses, and (iv) provide insights into the diversification history of cycads.

Results

Using Bayes factors, we compared divergence time estimates and the inferred dynamics of diversification when using Yule versus birth-death priors. Bayes factors were calculated with marginal likelihood estimated with stepping-stone sampling. We found striking differences in age estimates and diversification dynamics depending on prior choice. Dating with the Yule prior suggested that extant cycad genera diversified in the Paleogene and with two diversification rate shifts. In contrast, dating with the birth-death prior yielded Neogene diversifications, and four rate shifts, one for each of the four richest genera. Nonetheless, dating with the two priors provided similar age estimates for the divergence of cycads from Ginkgo (Carboniferous) and their crown age (Permian). Of these, Bayes factors clearly supported the birth-death prior.

Conclusions

These results suggest the choice of the branching process prior can have a drastic influence on our understanding of evolutionary radiations. Therefore, all dating analyses must involve a model selection process using Bayes factors to select between a Yule or birth-death prior, in particular on ancient clades with a potential pattern of high extinction. We also provide new insights into the history of cycad diversification because we found (i) periods of extinction along the long branches of the genera consistent with fossil data, and (ii) high diversification rates within the Miocene genus radiations.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0347-8) contains supplementary material, which is available to authorized users.

Comparative historical biogeography of three groups of Nearctic freshwater fishes across central Mexico

Biogeographic patterns of the three main Nearctic groups of continental fishes inhabiting river drainages in central Mexico (livebearing goodeids, southern Mexican notropins and species of Algansea, the last two representing independent lineages of cyprinids) were obtained and compared by following two approaches: an estimate of divergence times and using a well-defined biogeographic method. Three concordant biogeographic events were identified among the three groups, showing some evidence of a partially congruent evolutionary history. The analysed groups show at least three independent colonization events into central Mexico: two western routes, followed by the Goodeinae and members of Algansea, and an early Plateau route followed by southern notropins. The most recent common ancestor (MRCA) of each of the three freshwater fish groups diversified in central Mexico in the Late Miocene. The lack of a strong congruence in their biogeographic patterns, and the differences in species richness among the three clades might be evidence for distinct patterns of diversification.