Cenozoic history of the tropical marine biodiversity hotspot.

The region with the highest marine biodiversity on our planet is known as the Coral Triangle or Indo-Australian Archipelago (IAA)1,2. Its enormous biodiversity has long attracted the interest of biologists; however, the detailed evolutionary history of the IAA biodiversity hotspot remains poorly understood3. Here we present a high-resolution reconstruction of the Cenozoic diversity history of the IAA by inferring speciation-extinction dynamics using a comprehensive fossil dataset. We found that the IAA has exhibited a unidirectional diversification trend since about 25 million years ago, following a roughly logistic increase until a diversity plateau beginning about 2.6 million years ago. The growth of diversity was primarily controlled by diversity dependency and habitat size, and also facilitated by the alleviation of thermal stress after 13.9 million years ago. Distinct net diversification peaks were recorded at about 25, 20, 16, 12 and 5 million years ago, which were probably related to major tectonic events in addition to climate transitions. Key biogeographic processes had far-reaching effects on the IAA diversity as shown by the long-term waning of the Tethyan descendants versus the waxing of cosmopolitan and IAA taxa. Finally, it seems that the absence of major extinctions and the Cenozoic cooling have been essential in making the IAA the richest marine biodiversity hotspot on Earth.

The Angiosperm Terrestrial Revolution buffered ants against extinction.

With ~14,000 extant species, ants are ubiquitous and of tremendous ecological importance. They have undergone remarkable diversification throughout their evolutionary history. However, the drivers of their diversity dynamics are not well quantified or understood. Previous phylogenetic analyses have suggested patterns of diversity dynamics associated with the Angiosperm Terrestrial Revolution (ATR), but these studies have overlooked valuable information from the fossil record. To address this gap, we conducted a comprehensive analysis using a large dataset that includes both the ant fossil record (~24,000 individual occurrences) and neontological data (~14,000 occurrences), and tested four hypotheses proposed for ant diversification: co-diversification, competitive extinction, hyper-specialization, and buffered extinction. Taking into account biases in the fossil record, we found three distinct diversification periods (the latest Cretaceous, Eocene, and Oligo-Miocene) and one extinction period (Late Cretaceous). The competitive extinction hypothesis between stem and crown ants is not supported. Instead, we found support for the co-diversification, buffered extinction, and hyper-specialization hypotheses. The environmental changes of the ATR, mediated by the angiosperm radiation, likely played a critical role in buffering ants against extinction and favoring their diversification by providing new ecological niches, such as forest litter and arboreal nesting sites, and additional resources. We also hypothesize that the decline and extinction of stem ants during the Late Cretaceous was due to their hyper-specialized morphology, which limited their ability to expand their dietary niche in changing environments. This study highlights the importance of a holistic approach when studying the interplay between past environments and the evolutionary trajectories of organisms.

Convergence in sympatric swallowtail butterflies reveals ecological interactions as a key driver of worldwide trait diversification.

Ecological interactions can promote phenotypic diversification in sympatric species. While competition can enhance trait divergence, other ecological interactions may promote convergence in sympatric species. Within butterflies, evolutionary convergences in wing color patterns have been reported between distantly related species, especially in females of palatable species, where mimetic color patterns are promoted by predator communities shared with defended species living in sympatry. Wing color patterns are also often involved in species recognition in butterflies, and divergence in this trait has been reported in closely related species living in sympatry as a result of reproductive character displacement. Here, we investigate the effect of sympatry between species on the convergence vs. divergence of their wing color patterns in relation to phylogenetic distance, focusing on the iconic swallowtail butterflies (family Papilionidae). We developed an unsupervised machine learning-based method to estimate phenotypic distances among wing color patterns of 337 species, enabling us to finely quantify morphological diversity at the global scale among species and allowing us to compute pairwise phenotypic distances between sympatric and allopatric species pairs. We found phenotypic convergence in sympatry, stronger among distantly related species, while divergence was weaker and restricted to closely related males. The convergence was stronger among females than males, suggesting that differential selective pressures acting on the two sexes drove sexual dimorphism. Our results highlight the significant effect of ecological interactions driven by predation pressures on trait diversification in Papilionidae and provide evidence for the interaction between phylogenetic proximity and ecological interactions in sympatry, acting on macroevolutionary patterns of phenotypic diversification.

Drivers of the artiodactyl turnover in insular western Europe at the Eocene-Oligocene Transition.

Simultaneously investigating the effects of abiotic and biotic factors on diversity dynamics is essential to understand the evolutionary history of clades. The Grande Coupure corresponds to a major faunal turnover at the Eocene-Oligocene transition (EOT) (~34.1 to 33.55 Mya) and is defined in western Europe as an extinction of insular European mammals coupled with the arrival of crown clades from Asia. Here, we focused on the species-rich group of endemic European artiodactyls to determine the drivers of the Grande Coupure during the major environmental disruptions at the EOT. Using Bayesian birth-death models, we analyzed an original high-resolution fossil dataset (90 species, >2,100 occurrences) from southwestern France (Quercy area) and estimated the regional diversification and diversity dynamics of endemic and immigrant artiodactyls. We show that the endemic artiodactyl radiation was mainly related to the Eocene tropical conditions, combined with biotic controls on speciation and clade-related diversity dependence. We further highlight that the major environmental changes at the transition (77% of species became extinct) and the concurrent increase in seasonality in Europe during the Oligocene were likely the main drivers of their decline. Surprisingly, our results do not support the widely-held hypothesis of active competition between endemic and immigrant artiodactyls but rather suggest a passive or opportunistic replacement by immigrants, which is further supported by morphological clustering of specific ecological traits across the Eocene-Oligocene transition. Our analyses provide insights into the evolutionary and ecological processes driving the diversification and decline of mammalian clades during a major biological and climatic crisis.

An eosimiid primate of South Asian affinities in the Paleogene of Western Amazonia and the origin of New World monkeys.

Recent fossil discoveries in Western Amazonia revealed that two distinct anthropoid primate clades of African origin colonized South America near the Eocene/Oligocene transition (ca. 34 Ma). Here, we describe a diminutive fossil primate from Brazilian Amazonia and suggest that, surprisingly, a third clade of anthropoids was involved in the Paleogene colonization of South America by primates. This new taxon, Ashaninkacebus simpsoni gen. et sp. nov., has strong dental affinities with Asian African stem anthropoids: the Eosimiiformes. Morphology-based phylogenetic analyses of early Old World anthropoids and extinct and extant New World monkeys (platyrrhines) support relationships of both Ashaninkacebus and Amamria (late middle Eocene, North Africa) to the South Asian Eosimiidae. Afro-Arabia, then a mega island, played the role of a biogeographic stopover between South Asia and South America for anthropoid primates and hystricognathous rodents. The earliest primates from South America bear little adaptive resemblance to later Oligocene-early Miocene platyrrhine monkeys, and the scarcity of available paleontological data precludes elucidating firmly their affinities with or within Platyrrhini. Nonetheless, these data shed light on some of their life history traits, revealing a particularly small body size and a diet consisting primarily of insects and possibly fruit, which would have increased their chances of survival on a natural floating island during this extraordinary over-water trip to South America from Africa. Divergence-time estimates between Old and New World taxa indicate that the transatlantic dispersal(s) could source in the intense flooding events associated with the late middle Eocene climatic optimum (ca. 40.5 Ma) in Western Africa.

The origin and speciation of orchids.

Orchids constitute one of the most spectacular radiations of flowering plants. However, their origin, spread across the globe, and hotspots of speciation remain uncertain due to the lack of an up-to-date phylogeographic analysis. We present a new Orchidaceae phylogeny based on combined high-throughput and Sanger sequencing data, covering all five subfamilies, 17/22 tribes, 40/49 subtribes, 285/736 genera, and c. 7% (1921) of the 29 524 accepted species, and use it to infer geographic range evolution, diversity, and speciation patterns by adding curated geographical distributions from the World Checklist of Vascular Plants. The orchids' most recent common ancestor is inferred to have lived in Late Cretaceous Laurasia. The modern range of Apostasioideae, which comprises two genera with 16 species from India to northern Australia, is interpreted as relictual, similar to that of numerous other groups that went extinct at higher latitudes following the global climate cooling during the Oligocene. Despite their ancient origin, modern orchid species diversity mainly originated over the last 5 Ma, with the highest speciation rates in Panama and Costa Rica. These results alter our understanding of the geographic origin of orchids, previously proposed as Australian, and pinpoint Central America as a region of recent, explosive speciation.

Bayesian total-evidence dating revisits sloth phylogeny and biogeography: a cautionary tale on morphological clock analyses.

Combining morphological and molecular characters through Bayesian total-evidence dating allows inferring the phylogenetic and timescale framework of both extant and fossil taxa, while accounting for the stochasticity and incompleteness of the fossil record. Such an integrative approach is particularly needed when dealing with clades such as sloths (Mammalia: Folivora), for which developmental and biomechanical studies have shown high levels of morphological convergence whereas molecular data can only account for a limited percentage of their total species richness. Here, we propose an alternative hypothesis of sloth evolution that emphasizes the pervasiveness of morphological convergence and the importance of considering the fossil record and an adequate taxon sampling in both phylogenetic and biogeographic inferences. Regardless of different clock models and morphological datasets, the extant sloth Bradypus is consistently recovered as a megatherioid, and Choloepus as a mylodontoid, in agreement with molecular-only analyses. The recently extinct Caribbean sloths (Megalocnoidea) are found to be a monophyletic sister-clade of Megatherioidea, in contrast to previous phylogenetic hypotheses. Our results contradict previous morphological analyses and further support the polyphyly of "Megalonychidae", whose members were found in five different clades. Regardless of taxon sampling and clock models, the Caribbean colonization of sloths is compatible with the exhumation of islands along Aves Ridge and its geological time frame. Overall, our total-evidence analysis illustrates the difficulty of positioning highly incomplete fossils, although a robust phylogenetic framework was recovered by an a posteriori removal of taxa with high percentages of missing characters. Elimination of these taxa improved topological resolution by reducing polytomies and increasing node support. However, it introduced a systematic and geographic bias because most of these incomplete specimens are from northern South America. This is evident in biogeographic reconstructions, which suggest Patagonia as the area of origin of many clades when taxa are underrepresented, but Amazonia and/or Central and Southern Andes when all taxa are included. More generally, our analyses demonstrate the instability of topology and divergence time estimates when using different morphological datasets and clock models, and thus caution against making macroevolutionary inferences when node support is weak or when uncertainties in the fossil record are not considered.

Reconciling fossils with phylogenies reveals the origin and macroevolutionary processes explaining the global cycad biodiversity.

The determinants of biodiversity patterns can be understood using macroevolutionary analyses. The integration of fossils into phylogenies offers a deeper understanding of processes underlying biodiversity patterns in deep time. Cycadales are considered a relict of a once more diverse and globally distributed group but are restricted to low latitudes today. We still know little about their origin and geographic range evolution. Combining molecular data for extant species and leaf morphological data for extant and fossil species, we study the origin of cycad global biodiversity patterns through Bayesian total-evidence dating analyses. We assess the ancestral geographic origin and trace the historical biogeography of cycads with a time-stratified process-based model. Cycads originated in the Carboniferous on the Laurasian landmass and expanded in Gondwana in the Jurassic. Through now-vanished continental connections, Antarctica and Greenland were crucial biogeographic crossroads for cycad biogeography. Vicariance is an essential speciation mode in the deep and recent past. Their latitudinal span increased in the Jurassic and restrained toward subtropical latitudes in the Neogene in line with biogeographic inferences of high-latitude extirpations. We show the benefits of integrating fossils into phylogenies to estimate ancestral areas of origin and to study evolutionary processes explaining the global distribution of present-day relict groups.

Phylogeny and diversification of planthoppers (Hemiptera: Fulgoromorpha) based on a comprehensive molecular dataset and large taxon sampling.

Our understanding of the evolution of Fulgoromorpha (Insects, Hemiptera) has relied on molecular studies that have only considered either a limited number of taxa where all the families were not represented simultaneously, or a reduced number of genes.The absence of a global analysis comparing all the available data has thus led to significant biases in the analyzes, as evidenced by the incongruence of the results reported for planthopper phylogeny. Here we provide a phylogenetic and dating analysis of the Fulgoromorpha with a large sampling of 531 ingroup taxa, representing about 80% of the currently described suprageneric taxonomic diversity in this group. This study is based on most of the molecular sequences available to date and duly verified, for a set of nuclear and mitochondrial genes from a taxonomic sampling as complete as possible. The most significant results of our study are: (1) the unexpected paraphyly of Delphacidae whose Protodelphacida seem more related to Cixiidae than to other Delphacidae;(2) the group Meenoplidae-Kinnaridae recovered sister to the remaining Fulgoroidea families; (3) the early branching node of Tettigometridae sister of all the other families;(4) the Achilidae-Derbidae clade with Achilidae Plectoderini including Achilixiidae recovered as monophyletic as well as theFulgoridae-Dictyopharidae clade; and (5) the Tropiduchidae placed sister to the other so called 'higher' families (sec. Shcherbakov, 2006).Our divergence times analysis, calibrated with a set of duly verified fossils, suggests that the first diversification of planthoppers occurred in the Early Triassic around 240 Mya and those of the superfamilies Delphacoidea and Fulgoroidea in the Middle-Late Triassic around 210 Mya and 230 Mya, respectively. By the end of the Jurassic, all major planthopper lineages were originated, and all families, around 125 Mya, might havebeen driven in their distribution and evolution (in their first subfamilial divisions) by the geographical constraints of the Gondwanan break-up.Rapid evolutionary radiations occurred particularly in Fulgoridae around 125-130 Mya. Our results stress the importance of the good quality of the sequences used in the molecular analyzes and the primordial importance of a large sampling when analyzing the phylogeny of the group.

A comprehensive phylogeny and revised taxonomy illuminate the origin and diversification of the global radiation of Papilio (Lepidoptera: Papilionidae).

The swallowtail genus Papilio (Lepidoptera: Papilionidae) is species rich, distributed worldwide, and has broad morphological habits and ecological niches. Because of its elevated species richness, it has been historically difficult to reconstruct a densely sampled phylogeny for this clade. Here we provide a taxonomic working list for the genus, resulting in 235 Papilio species, and assemble a molecular dataset of seven gene fragments representing ca. 80% of the currently described diversity. Phylogenetic analyses reconstructed a robust tree with highly supported relationships within subgenera, although a few nodes in the early history of the Old World Papilio remain unresolved. Contrasting with previous results, we found that Papilio alexanor is sister to all Old World Papilio and that the subgenus Eleppone is no longer monotypic. The latter includes the recently described Fijian Papilio natewa with the Australian Papilio anactus and is sister to subgenus Araminta (formerly included in subgenus Menelaides) occurring in Southeast Asia. Our phylogeny also includes rarely studied (P. antimachus, P. benguetana) or endangered species (P. buddha, P. chikae). Taxonomic changes resulting from this study are elucidated. Molecular dating and biogeographic analyses indicate that Papilio originated ca. 30 million years ago (Oligocene), in a northern region centered on Beringia. A rapid early Miocene radiation in the Paleotropics is revealed within Old World Papilio, potentially explaining their low early branch support. Most subgenera originated in the early to middle Miocene followed by synchronous southward biogeographic dispersals and repeated local extirpations in northern latitudes. This study provides a comprehensive phylogenetic framework for Papilio with clarification of subgeneric systematics and species taxonomic changes enumerated, which will facilitate further studies to address questions on their ecology and evolutionary biology using this model clade.

Phylogenomic and Macroevolutionary Evidence for an Explosive Radiation of a Plant Genus in the Miocene.

Mountain systems harbor a substantial fraction of global biodiversity and, thus, provide excellent opportunities to study rapid diversification and to understand the historical processes underlying the assembly of biodiversity hotspots. The rich biodiversity in mountains is widely regarded as having arisen under the influence of geological and climatic processes as well as the complex interactions among them. However, the relative contribution of geology and climate in driving species radiation is seldom explored. Here, we studied the evolutionary radiation of Oreocharis (Gesneriaceae), which has diversified extensively throughout East Asia, especially within the Hengduan Mountains (HDM), using transcriptomic data and a time calibrated phylogeny for 88% (111/126) of all species of the genus. In particular, we applied phylogenetic reconstructions to evaluate the extent of incomplete lineage sorting accompanying the early and rapid radiation in the genus. We then fit macroevolutionary models to explore its spatial and diversification dynamics in Oreocharis and applied explicit birth-death models to investigate the effects of past environmental changes on its diversification. Evidence from 574 orthologous loci suggest that Oreocharis underwent an impressive early burst of speciation starting ca. 12 Ma in the Miocene, followed by a drastic decline in speciation toward the present. Although we found no evidence for a shift in diversification rate across the phylogeny of Oreocharis, we showed a difference in diversification dynamics between the HDM and non-HDM lineages, with higher diversification rates in the HDM. The diversification dynamic of Oreocharis is most likely positively associated with temperature-dependent speciation and dependency on the Asian monsoons. We suggest that the warm and humid climate of the mid-Miocene was probably the primary driver of the rapid diversification in Oreocharis, while mountain building of the HDM might have indirectly affected species diversification of the HDM lineage. This study highlights the importance of past climatic changes, combined with mountain building, in creating strong environmental heterogeneity and driving diversification of mountain plants, and suggests that the biodiversity in the HDM cannot directly be attributed to mountain uplift, contrary to many recent speculations.[East Asian monsoons; environmental heterogeneity; Hengduan Mountains; incomplete lineage sorting; Oreocharis; past climate change; rapid diversification; transcriptome.].

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

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.].

A sequential diversification with Miocene extinction and Pliocene speciation linked to mountain uplift explains the diversity of the African rain forest clade Monodoreae (Annonaceae).

BACKGROUND AND AIMS: Throughout the Cenozoic, Africa underwent several climatic and geological changes impacting the evolution of tropical rain forests (TRF). African TRF are thought to have extended from East to West in a 'pan-African' TRF, followed by several events of fragmentation during drier climate periods. During the Miocene, climate cooling and mountain uplift led to the aridification of tropical Africa and open habitats expanded at the expense of TRF, which likely experienced local extinctions. However, in plants, these drivers were previously inferred using limited taxonomic and molecular data. Here, we tested the impact of climate and geological changes on diversification within the diverse clade Monodoreae (Annonaceae) composed of 90 tree species restricted to African TRF. METHODS: We reconstructed a near complete phylogenetic tree, based on 32 nuclear genes, and dated using relaxed clocks and fossil calibrations in a Bayesian framework. We inferred the biogeographic history and the diversification dynamics of the clade using multiple birth-death models. KEY RESULTS: Monodoreae originated in East African TRF ca. 25 million years ago (Ma) and expanded toward Central Africa during the Miocene. We inferred range contractions during the middle Miocene and document important connections between East and West African TRF after 15-13 Ma. Our results indicated a sudden extinction event during the late Miocene, followed by an increase in speciation rates. Birth-death models suggested that African elevation change (orogeny) is positively linked to speciation in this clade. CONCLUSION: East Africa is inferred as an important source of Monodoreae species, and possibly for African plant diversity in general. Our results support a "sequential scenario of diversification" where increased aridification triggered extinction of TRF species in Monodoreae. This was quickly followed by rain forests fragmentation, subsequently enhancing lagged speciation resulting from vicariance and improved climate conditions. In contrast to previous ideas, the uplift of East Africa is shown to have played a positive role in Monodoreae diversification.

The rise of angiosperms pushed conifers to decline during global cooling.

Competition among species and entire clades can impact species diversification and extinction, which can shape macroevolutionary patterns. The fossil record shows successive biotic turnovers such that a dominant group is replaced by another. One striking example involves the decline of gymnosperms and the rapid diversification and ecological dominance of angiosperms in the Cretaceous. It is generally believed that angiosperms outcompeted gymnosperms, but the macroevolutionary processes and alternative drivers explaining this pattern remain elusive. Using extant time trees and vetted fossil occurrences for conifers, we tested the hypotheses that clade competition or climate change led to the decline of conifers at the expense of angiosperms. Here, we find that both fossil and molecular data show high congruence in revealing 1) low diversification rates, punctuated by speciation pulses, during warming events throughout the Phanerozoic and 2) that conifer extinction increased significantly in the Mid-Cretaceous (100 to 110 Ma) and remained high ever since. Their extinction rates are best explained by the rise of angiosperms, rejecting alternative models based on either climate change or time alone. Our results support the hypothesis of an active clade replacement, implying that direct competition with angiosperms increased the extinction of conifers by pushing their remaining species diversity and dominance out of the warm tropics. This study illustrates how entire branches on the Tree of Life may actively compete for ecological dominance under changing climates.

A phylogenetic study to assess the link between biome specialization and diversification in swallowtail butterflies.

The resource-use hypothesis, proposed by E.S. Vrba, states that habitat fragmentation caused by climatic oscillations would affect particularly biome specialists (species inhabiting only one biome), which might show higher speciation and extinction rates than biome generalists. If true, lineages would accumulate biome-specialist species. This effect would be particularly exacerbated for biomes located at the periphery of the global climatic conditions, namely, biomes that have high/low precipitation and high/low temperature such as rainforest (warm-humid), desert (warm-dry), steppe (cold-dry) and tundra (cold-humid). Here, we test these hypotheses in swallowtail butterflies, a clade with more than 570 species, covering all the continents but Antarctica, and all climatic conditions. Swallowtail butterflies are among the most studied insects, and they are a model group for evolutionary biology and ecology studies. Continental macroecological rules are normally tested using vertebrates, this means that there are fewer examples exploring terrestrial invertebrate patterns at global scale. Here, we compiled a large Geographic Information System database on swallowtail butterflies' distribution maps and used the most complete time-calibrated phylogeny to quantify diversification rates (DRs). In this paper, we aim to answer the following questions: (1) Are there more biome-specialist swallowtail butterflies than biome generalists? (2) Is DR related to biome specialization? (3) If so, do swallowtail butterflies inhabiting extreme biomes show higher DRs? (4) What is the effect of species distribution area? Our results showed that swallowtail family presents a great number of biome specialists which showed substantially higher DRs compared to generalists. We also found that biome specialists are unevenly distributed across biomes. Overall, our results are consistent with the resource-use hypothesis, species climatic niche and biome fragmentation as key factors promoting isolation.

Impact of Pleistocene Eustatic Fluctuations on Evolutionary Dynamics in Southeast Asian Biodiversity Hotspots.

Pleistocene climatic fluctuations (PCF) are frequently highlighted as important evolutionary engines that triggered cycles of biome expansion and contraction. Although there is ample evidence of the impact of PCF on biodiversity of continental biomes, the consequences in insular systems depend on the geology of the islands and the ecology of the taxa inhabiting them. The idiosyncratic aspects of insular systems are exemplified by the islands of the Sunda Shelf in Southeast Asia (Sundaland), where PCF-induced eustatic fluctuations had complex interactions with the geology of the region, resulting in high species diversity and endemism. Emergent land in Southeast Asia varied drastically with sea-level fluctuations during the Pleistocene. Climate-induced fluctuations in sea level caused temporary connections between insular and continental biodiversity hotspots in Southeast Asia. These exposed lands likely had freshwater drainage systems that extended between modern islands: the Paleoriver Hypothesis. Built upon the assumption that aquatic organisms are among the most suitable models to trace ancient river boundaries and fluctuations of landmass coverage, the present study aims to examine the evolutionary consequences of PCF on the dispersal of freshwater biodiversity in Southeast Asia. Time-calibrated phylogenies of DNA-delimited species were inferred for six species-rich freshwater fish genera in Southeast Asia (Clarias, Channa, Glyptothorax, Hemirhamphodon, Dermogenys, Nomorhamphus). The results highlight rampant cryptic diversity and the temporal localization of most speciation events during the Pleistocene, with 88% of speciation events occurring during this period. Diversification analyses indicate that sea-level-dependent diversification models poorly account for species proliferation patterns for all clades excepting Channa. Ancestral area estimations point to Borneo as the most likely origin for most lineages, with two waves of dispersal to Sumatra and Java during the last 5 myr. Speciation events are more frequently associated with boundaries of the paleoriver watersheds, with 60%, than islands boundaries, with 40%. In total, one-third of speciation events are inferred to have occurred within paleorivers on a single island, suggesting that habitat heterogeneity and factors other than allopatry between islands substantially affected diversification of Sundaland fishes. Our results suggest that species proliferation in Sundaland is not wholly reliant on Pleistocene sea-level fluctuations isolating populations on different islands. [Dispersal; diversification; eustatic fluctuations; freshwater fishes; insular systems; Milankovitch cycles; paleoenvironments; vicariance.].

Climate cooling and clade competition likely drove the decline of lamniform sharks.

Understanding heterogeneity in species richness between closely related clades is a key research question in ecology and evolutionary biology. Multiple hypotheses have been proposed to interpret such diversity contrasts across the tree of life, with most studies focusing on speciation rates to explain clades' evolutionary radiations, while often neglecting extinction rates. Here we study a notorious biological model as exemplified by the sister relationships between mackerel sharks (Lamniformes, 15 extant species) and ground sharks (Carcharhiniformes, ∼290 extant species). Using a comprehensive fossil dataset, we found that the diversity dynamics of lamniforms waxed and waned following repeated cycles of radiation phases and declining phases. Radiation phases peaked up to 3 times the current diversity in the early Late Cretaceous. In the last 20 million years, the group declined to its present-day diversity. Along with a higher extinction risk for young species, we further show that this declining pattern is likely attributed to a combination of abiotic and biotic factors, with a cooling-driven extinction (negative correlation between temperature and extinction) and clade competition with some ground sharks. Competition from multiple clades successively drove the demise and replacement of mackerel sharks due to a failure to originate facing the rise of ground sharks, particularly since the Eocene. These effects came from ecologically similar carcharhiniform species inhibiting diversification of medium- and large-sized lamniforms. These results imply that the interplay between abiotic and biotic drivers had a substantial role in extinction and speciation, respectively, which determines the sequential rise and decline of marine apex predators.

Punctuational ecological changes rather than global factors drive species diversification and the evolution of wing phenotypes in Morpho butterflies.

Assessing the relative importance of geographical and ecological drivers of evolution is paramount to understand the diversification of species and traits at the macroevolutionary scale. Here, we use an integrative approach, combining phylogenetics, biogeography, ecology and quantified phenotypes to investigate the drivers of both species and phenotypic diversification of the iconic Neotropical butterfly genus Morpho. We generated a time-calibrated phylogeny for all known species and inferred historical biogeography. We fitted models of time-dependent (accounting for rate heterogeneity across the phylogeny) and paleoenvironment-dependent diversification (accounting for global effect on the phylogeny). We used geometric morphometrics to assess variation of wing size and shape across the tree and investigated their dynamics of evolution. We found that the diversification of Morpho is best explained when considering variable diversification rates across the tree, possibly associated with lineages occupying different microhabitat conditions. First, a shift from understory to canopy was characterized by an increased speciation rate partially coupled with an increasing rate of wing shape evolution. Second, the occupation of dense bamboo thickets accompanying a major host-plant shift from dicotyledons towards monocotyledons was associated with a simultaneous diversification rate shift and an evolutionary 'jump' of wing size. Our study points to a diversification pattern driven by punctuational ecological changes instead of a global driver or biogeographic history.

Whole Genome Shotgun Phylogenomics Resolves the Pattern and Timing of Swallowtail Butterfly Evolution.

Evolutionary relationships have remained unresolved in many well-studied groups, even though advances in next-generation sequencing and analysis, using approaches such as transcriptomics, anchored hybrid enrichment, or ultraconserved elements, have brought systematics to the brink of whole genome phylogenomics. Recently, it has become possible to sequence the entire genomes of numerous nonbiological models in parallel at reasonable cost, particularly with shotgun sequencing. Here, we identify orthologous coding sequences from whole-genome shotgun sequences, which we then use to investigate the relevance and power of phylogenomic relationship inference and time-calibrated tree estimation. We study an iconic group of butterflies-swallowtails of the family Papilionidae-that has remained phylogenetically unresolved, with continued debate about the timing of their diversification. Low-coverage whole genomes were obtained using Illumina shotgun sequencing for all genera. Genome assembly coupled to BLAST-based orthology searches allowed extraction of 6621 orthologous protein-coding genes for 45 Papilionidae species and 16 outgroup species (with 32% missing data after cleaning phases). Supermatrix phylogenomic analyses were performed with both maximum-likelihood (IQ-TREE) and Bayesian mixture models (PhyloBayes) for amino acid sequences, which produced a fully resolved phylogeny providing new insights into controversial relationships. Species tree reconstruction from gene trees was performed with ASTRAL and SuperTriplets and recovered the same phylogeny. We estimated gene site concordant factors to complement traditional node-support measures, which strengthens the robustness of inferred phylogenies. Bayesian estimates of divergence times based on a reduced data set (760 orthologs and 12% missing data) indicate a mid-Cretaceous origin of Papilionoidea around 99.2 Ma (95% credibility interval: 68.6-142.7 Ma) and Papilionidae around 71.4 Ma (49.8-103.6 Ma), with subsequent diversification of modern lineages well after the Cretaceous-Paleogene event. These results show that shotgun sequencing of whole genomes, even when highly fragmented, represents a powerful approach to phylogenomics and molecular dating in a group that has previously been refractory to resolution.