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

The last of their kind: Is the genus Scutiger (Anura: Megophryidae) a relict element of the paleo-transhimalaya biota?

The orographic evolution of the Himalaya-Tibet Mountain system continues to be a subject of controversy, leading to considerable uncertainty regarding the environment and surface elevation of the Tibetan Plateau during the Cenozoic era. As many geoscientific (but not paleontological) studies suggest, elevations close to modern heights exist in vast areas of Tibet since at least the late Paleogene, implicating the presence of large-scale alpine environments for more than 30 million years. To explore a recently proposed alternative model that assumes a warm temperate environment across paleo-Tibet, we carried out a phylogeographic survey using genomic analyses of samples covering the range of endemic lazy toads (Scutiger) across the Himalaya-Tibet orogen. We identified two main clades, with several, geographically distinct subclades. The long temporal gap between the stem and crown age of Scutiger may suggest high extinction rates. Diversification within the crown group, depending on the calibration, occurred either from the Mid-Miocene or Late-Miocene and continued until the Holocene. The present-day Himalayan Scutiger fauna could have evolved from lineages that existed on the southern edges of the paleo-Tibetan area (the Transhimalaya = Gangdese Shan), while extant species living on the eastern edge of the Plateau originated probably from the eastern edges of northern parts of the ancestral Tibetan area (Hoh Xil, Tanggula Shan). Based on the Mid-Miocene divergence time estimation and ancestral area reconstruction, we propose that uplift-associated aridification of a warm temperate Miocene-Tibet, coupled with high extirpation rates of ancestral populations, and species range shifts along drainage systems and epigenetic transverse valleys of the rising mountains, is a plausible scenario explaining the phylogenetic structure of Scutiger. This hypothesis aligns with the fossil record but conflicts with geoscientific concepts of high elevated Tibetan Plateau since the late Paleogene. Considering a Late-Miocene/Pliocene divergence time, an alternative scenario of dispersal from SE Asia into the East, Central, and West Himalaya cannot be excluded, although essential evolutionary and biogeographic aspects remain unresolved within this model.

Reliable biogeography requires fossils: insights from a new species-level phylogeny of extinct and living carnivores

A central objective of historical biogeography is to understand where clades originated and how they moved across space and over time. However, given the dynamic history of ecosystem changes in response to climate change and geological events, the manifold long-distance dispersals over evolutionary timescales, and regional and global extinctions, it remains uncertain how reliable inferences based solely on extant taxa can be achieved. Using a novel species-level phylogeny of all known extant and extinct species of the mammalian order Carnivora and related extinct groups, we show that far more precise and accurate ancestral areas can be estimated by fully integrating extinct species into the analyses, rather than solely relying on extant species or identifying ancestral areas only based on the geography of the oldest fossils. Through a series of simulations, we further show that this conclusion is robust under realistic scenarios in which the unknown extinct taxa represent a biased subset of all extinct species. Our results highlight the importance of integrating fossil taxa into a phylogenetic framework to further improve our understanding of historical biogeography and reveal the dynamic dispersal and diversification history of carnivores.

A story of vicariance? How the geology of oceanic archipelagos influenced the evolutionary history of endemic damselflies

South Pacific islands provide an ideal study system to explore patterns of speciation, specifically examining the role of dispersal versus vicariance. Dispersal is often the suggested mechanism of diversification in the South Pacific, specifically among remote island chains. Here, we provide a phylogeny of several related genera of Coenagrionidae (Odonata: Zygoptera) from the South Pacific, based on five molecular loci, in order to examine patterns of speciation in the region. We used the endemic damselfly genera Nesobasis, Nikoulabasis, and Vanuatubasis found across both Fiji and Vanuatu. Knowledge of the geologic history of the region was used to inform our understanding of the evolution of these genera. Both archipelagos used to be part of the Vitiaz arc which spanned from the Solomon Islands to Tonga and began to break apart 10-12 Ma. Results of our divergence-time estimations and biogeographic reconstructions support that the breakup of this arc acted as a significant vicariance event in the evolution of these taxa. Specifically, it led to the extant generic diversity seen in these damselflies. We find that within the archipelago of Vanuatu, that Espiritu Santo served as an important source for dispersal to other islands with Malekula acting as a stepping stone to Efate.

Origin and evolutionary history of Populus (Salicaceae): Further insights based on time divergence and biogeographic analysis

INTRODUCTION

Populus (Salicaceae) species harbour rich biodiversity and are widely distributed throughout the Northern Hemisphere. However, the origin and biogeography of Populus remain poorly understood.

METHODS

We infer the divergence times and the historical biogeography of the genus Populus through phylogenetic analysis of 34 chloroplast fragments based on a large sample.

RESULTS AND DISCUSSION

Eurasia is the likely location of the early divergences of Salicaceae after the Cretaceous-Paleogene (K-Pg) mass extinction, followed by recurrent spread to the remainder of the Old World and the New World beginning in the Eocene; the extant Populus species began to diversity during the early Oligocene (approximately 27.24 Ma), climate changes during the Oligocene may have facilitated the diversification of modern poplar species; three separate lineages of Populus from Eurasia colonized North America in the Cenozoic via the Bering Land Bridges (BLB); We hypothesize that the present day disjunction in Populus can be explained by two scenarios: (i) Populus likely originated in Eurasia and subsequently colonized other regions, including North America; and (ii) the fact that the ancestor of the genus Populus that was once widely distributed in the Northern Hemisphere and eventually wiped out due to the higher extinction rates in North America, similar to the African Rand flora. We hypothesize that disparities in extinction across the evolutionary history of Populus in different regions shape the modern biogeography of Populus. Further studies with dense sampling and more evidence are required to test these hypotheses. Our research underscores the significance of combining phylogenetic analyses with biogeographic interpretations to enhance our knowledge of the origin, divergence, and distribution of biodiversity in temperate plant floras.

Mapping species diversification metrics in macroecology: Prospects and challenges

The intersection of macroecology and macroevolution is one of today’s most active research in biology. In the last decade, we have witnessed a steady increment of macroecological studies that use metrics attempting to capture macroevolutionary processes to explain present-day biodiversity patterns. Evolutionary explanations of current species richness gradients are fundamental for understanding how diversity accumulates in a region. Although multiple hypotheses have been proposed to explain the patterns we observe in nature, it is well-known that the present-day diversity patterns result from speciation, extinction, colonization from nearby areas, or a combination of these macroevolutionary processes. Whether these metrics capture macroevolutionary processes across space is unknown. Some tip-rate metrics calculated directly from a phylogenetic tree (e.g., mean root distance -MRD-; mean diversification rate -mDR-) seem to return very similar geographical patterns regardless of how they are estimated (e.g., using branch lengths explicitly or not). Model-based tip-rate metrics —those estimated using macroevolutionary mixtures, e.g., the BAMM approach— seem to provide better net diversification estimates than only speciation rates. We argue that the lack of appropriate estimates of extinction and dispersal rates in phylogenetic trees may strongly limit our inferences about how species richness gradients have emerged at spatial and temporal scales. Here, we present a literature review about this topic and empirical comparisons between select taxa with several of these metrics. We implemented a simple null model approach to evaluate whether mapping of these metrics deviates from a random sampling process. We show that phylogenetic metrics by themselves are relatively poor at capturing speciation, extinction, and dispersal processes across geographical gradients. Furthermore, we provide evidence of how parametric biogeographic methods can improve our inference of past events and, therefore, our conclusions about the evolutionary processes driving biodiversity patterns. We recommend that further studies include several approaches simultaneously (e.g., spatial diversification modeling, parametric biogeographic methods, simulations) to disentangle the relative role of speciation, extinction, and dispersal in the generation and maintenance of species richness gradients at regional and global scales.

Estimating ancient biogeographic patterns with statistical model discrimination

The geographic ranges in which species live is a function of many factors underlying ecological and evolutionary contingencies. Observing the geographic range of an individual species provides valuable information about these historical contingencies for a lineage, determining the distribution of many distantly related species in tandem provides information about large-scale constraints on evolutionary and ecological processes generally. We present a linear regression method that allows for the discrimination of various hypothetical biogeographical models for determining which landscape distributional pattern best matches data from the fossil record. The linear regression models used in the discrimination rely on geodesic distances between sampling sites (typically geologic formations) as the independent variable and three possible dependent variables: Dice/Sorensen similarity; Euclidean distance; and phylogenetic community dissimilarity. Both the similarity and distance measures are useful for full-community analyses without evolutionary information, whereas the phylogenetic community dissimilarity requires phylogenetic data. Importantly, the discrimination method uses linear regression residual error to provide relative measures of support for each biogeographical model tested, not absolute answers or p-values. When applied to a recently published dataset of Campanian pollen, we find evidence that supports two plant communities separated by a transitional zone of unknown size. A similar case study of ceratopsid dinosaurs using phylogenetic community dissimilarity provided no evidence of a biogeographical pattern, but this case study suffers from a lack of data to accurately discriminate and/or too much temporal mixing. Future research aiming to reconstruct the distribution of organisms across a landscape has a statistical-based method for determining what biogeographic distributional model best matches the available data.

Wallacean and Melanesian Islands Promote Higher Rates of Diversification within the Global Passerine radiation Corvides

The complex island archipelagoes of Wallacea and Melanesia have provided empirical data behind integral theories in evolutionary biology, including allopatric speciation and island biogeography. Yet, questions regarding the relative impact of the layered biogeographic barriers, such as deep-water trenches and isolated island systems, on faunal diversification remain underexplored. One such barrier is Wallace's Line, a significant biogeographic boundary that largely separates Australian and Asian biodiversity. To assess the relative roles of biogeographic barriers-specifically isolated island systems and Wallace's Line-we investigated the tempo and mode of diversification in a diverse avian radiation, Corvides (Crows and Jays, Birds-of-paradise, Vangas, and allies). We combined a genus-level dataset of thousands of ultraconserved elements (UCEs) and a species-level, 12-gene Sanger sequence matrix to produce a well-resolved supermatrix tree that we leveraged to explore the group's historical biogeography and effects of biogeographic barriers on their macroevolutionary dynamics. The tree is well-resolved and differs substantially from what has been used extensively for past comparative analyses within this group. We confirmed that Corvides, and its major constituent clades, arose in Australia and that a burst of dispersals west across Wallace's Line occurred after the uplift of Wallacea during the mid-Miocene. We found that dispersal across this biogeographic barrier were generally rare, though westward dispersals were two times more frequent than eastward dispersals. Wallacea's central position between Sundaland and Sahul no doubt acted as a bridge for island-hopping dispersal out of Australia, across Wallace's Line, to colonize the rest of Earth. In addition, we found that the complex island archipelagoes east of Wallace's Line harbor the highest rates of net diversification and are a substantial source of colonists to continental systems on both sides of this biogeographic barrier. Our results support emerging evidence that island systems, particularly the geologically complex archipelagoes of the Indo-pacific, are drivers of species diversification.

Phylogeny of Merlin's grass (Isoetaceae): revealing an "Amborella syndrome" and the importance of geographic distribution for understanding current and historical diversity

Background

Merlin’s grass (Isoetes, Isoetaceae, Lycopsida), is the extant remnant of the isoetalean wood-producing lycopsids that originated during the Paleozoic, possibly in aquatic or boggy habitats. Modern day species are aquatic, semi-aquatic or terrestrial and occur almost worldwide. They display little morphological variation; the lobed corm has helically arranged leaves with internal air channels and basal sporangia. Genetic variation has also proven limited, which has hampered phylogenetic inference. We investigate evolutionary relationships in Isoetes, using molecular data and an extended sample of species compared to previous work, adding species that have never before been included in a phylogenetic study.

Results

Our results reveal an unexpected discovery of an “Amborella syndrome” in Isoetaceae: a single poorly known species is sister to the remaining family. The species, Isoetes wormaldii, is a rare endemic to the Eastern Cape of South Africa. Its leaves are flattened with a rounded point, which sharply contrasts with the awl-shaped leaves of most other species of Isoetes. The remaining species of Isoetes are resolved in five major clades, also indicated in previous work. While the phylogeny shows geographic structure, the patterns are complex. For example, tropical-southern African species occur in at least five clades, and Indian, Australian and Mediterranean species in at least three clades each.

Conclusion

The evolutionary and biogeographical history of Isoetes is not easily explained, and may conceivably include ample extinction and a mixture of ancient and more recent processes. Previously shown difficulties with node age estimation increase the problem. The here demonstrated sister-relationship between the phylogenetically, morphologically and genetically distinct Isoetes wormaldii and the remaining family appears to bridge the morphological gap between Isoetes and its extinct relatives, although further studies are needed. Moreover, it shortens the branch length to its living sister genus Selaginella, and may enhance node age estimation in future studies. Isoetes wormaldii is critically endangered, known only from one (to a few) minor populations. Immediate actions need to be taken if we want to prevent this unique species from going extinct.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-022-01988-w.

Cyphoderia ampulla (Cyphoderiidae: Rhizaria), a tale of freshwater sailors. The causes and consequences of ecological transitions through the salinity barrier in a family of benthic protists

The salinity barrier that separates marine and freshwater biomes is probably the most important division in biodiversity on Earth. Those organisms that successfully performed this transition had access to new ecosystems while undergoing changes in selective pressure, which often led to major shifts in diversification rates. While these transitions have been extensively investigated in animals, the tempo, mode, and outcome of crossing the salinity barrier have been scarcely studied in other eukaryotes. Here, we reconstructed the evolutionary history of the species complex Cyphoderia ampulla (Euglyphida: Cercozoa: Rhizaria) based on DNA sequences from the nuclear SSU rRNA gene and the mitochondrial cytochrome oxidase subunit I gene, obtained from publicly available environmental DNA data (GeneBank, EukBank) and isolated organisms. A tree calibrated with euglyphid fossils showed that four independent transitions towards freshwater systems occurred from the Mid Miocene onwards, coincident with important fluctuations in sea level. Ancestral trait reconstructions indicated that the whole family Cyphoderiidae had a marine origin and suggest that ancestors of the freshwater forms were euryhaline and lived in environments with fluctuating salinity. Diversification rates did not show any obvious increase concomitant with ecological transitions, but morphometric analyses indicated that species increased in size and homogenized their morphology after colonizing the new environments. This suggests adaptation to changes in selective pressure exerted by life in freshwater sediments.

Mitogenomics and hidden-trait models reveal the role of phoresy and host shifts in the diversification of parasitoid blister beetles (Coleoptera: Meloidae)

Changes in life history traits are often considered speciation triggers and can have dramatic effects on the evolutionary history of a lineage. Here, we examine the consequences of changes in two life history traits, host‐type and phoresy, in the hypermetamorphic blister beetles, Meloidae. Subfamilies Nemognathinae and Meloinae exhibit a complex life cycle involving multiple metamorphoses and parasitoidism. Most genera and tribes are bee‐parasitoids, and include phoretic or nonphoretic species, while two tribes feed on grasshopper eggs. These different life strategies are coupled with striking differences in species richness among clades. We generated a mitogenomic phylogeny for Nemognathinae and Meloinae, confirming the monophyly of these two clades, and used the dated phylogeny to explore the association between diversification rates and changes in host specificity and phoresy, using state‐dependent speciation and extinction (SSE) models that include the effect of hidden traits. To account for the low taxon sampling, we implemented a phylogenetic‐taxonomic approach based on birth‐death simulations, and used a Bayesian framework to integrate parameter and phylogenetic uncertainty. Results show that the ancestral hypermetamorphic Meloidae was a nonphoretic bee‐parasitoid, and that transitions towards a phoretic bee‐parasitoid and grasshopper parasitoidism occurred multiple times. Nonphoretic bee‐parasitoid lineages exhibit significantly higher relative extinction and lower diversification rates than phoretic bee‐and grasshopper‐parasitoids, but no significant differences were found between the latter two strategies. This suggests that Orthopteran host shifts and phoresy contributed jointly to the evolutionary success of the parasitoid meloidae. We also demonstrate that SSE models can be used to identify hidden traits coevolving with the focal trait in driving a lineage's diversification dynamics.

Predominantly Eastward Long-Distance Dispersal in Pantropical Ochnaceae Inferred From Ancestral Range Estimation and Phylogenomics

Ochnaceae is a pantropical family with multiple transoceanic disjunctions at deep and shallow levels. Earlier attempts to unravel the processes that led to such biogeographic patterns suffered from insufficient phylogenetic resolution and unclear delimitation of some of the genera. In the present study, we estimated divergence time and ancestral ranges based on a phylogenomic framework with a well-resolved phylogenetic backbone to tackle issues of the timing and direction of dispersal that may explain the modern global distribution of Ochnaceae. The nuclear data provided the more robust framework for divergence time estimation compared to the plastome-scale data, although differences in the inferred clade ages were mostly small. While Ochnaceae most likely originated in West Gondwana during the Late Cretaceous, all crown-group disjunctions are inferred as dispersal-based, most of them as transoceanic long-distance dispersal (LDD) during the Cenozoic. All LDDs occurred in an eastward direction except for the SE Asian clade of Sauvagesieae, which was founded by trans-Pacific dispersal from South America. The most species-rich clade by far, Ochninae, originated from either a widespread neotropical-African ancestor or a solely neotropical ancestor which then dispersed to Africa. The ancestors of this clade then diversified in Africa, followed by subsequent dispersal to the Malagasy region and tropical Asia on multiple instances in three genera during the Miocene-Pliocene. In particular, Ochna might have used the South Arabian land corridor to reach South Asia. Thus, the pantropical distribution of Ochnaceae is the result of LDD either transoceanic or via land bridges/corridors, whereas vicariance might have played a role only along the stem of the family.

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.

Biogeography Meets Niche Modeling: Inferring the Role of Deep Time Climate Change When Data Is Limited

Geographic range shifts are one major organism response to climate change, especially if the rate of climate change is higher than that of species adaptation. Ecological niche models (ENM) and biogeographic inferences are often used in estimating the effects of climatic oscillations on species range dynamics. ENMs can be used to track climatic suitable areas over time, but have often been limited to shallow timescales; biogeographic inference can reach greater evolutionary depth, but often lacks spatial resolution. Here, we present a simple approach that treats them as independent and complementary sources of evidence, which, when used in partnership, can be employed to reconstruct geographic range shifts over deep evolutionary timescales. For testing this, we chose two extreme African disjunctions: Camptoloma (Scrophulariaceae) and Canarina (Campanulaceae), each comprising of three species disjunctly distributed in Macaronesia and eastern/southern Africa. Using inferred ancestral ranges in tandem with preindustrial and paleoclimate ENM hindcastings, we show that the disjunct pattern was the result of fragmentation and extinction events linked to Neogene aridification cycles. Our results highlight the importance of considering temporal resolution when building ENMs for rare endemics with small population sizes and restricted climatic tolerances such as Camptoloma, for which models built on averaged monthly variables were more informative than those based on annual bioclimatic variables. Additionally, we show that biogeographic information can be used as truncation threshold criteria for building ENMs in the distant past. Our approach is suitable when there is sparse sampling on species occurrences and associated patterns of genetic variation, such as in the case of ancient endemics with widely disjunct distributions as a result of climate change.

Rare and widespread: integrating Bayesian MCMC approaches, Sanger sequencing and Hyb-Seq phylogenomics to reconstruct the origin of the enigmatic Rand Flora genus Camptoloma

PREMISE

Genera that are widespread, with geographically discontinuous distributions and represented by few species, are intriguing. Is their achieved disjunct distribution recent or ancient in origin? Why are they species-poor? The Rand Flora is a continental-scale pattern in which closely related species appear codistributed in isolated regions over the continental margins of Africa. Genus Camptoloma (Scrophulariaceae) is the most notable example, comprising three species isolated from each other on the northwest, eastern, and southwest Africa.

METHODS

We employed Sanger sequencing of nuclear and plastid markers, together with genomic target sequencing of 2190 low-copy nuclear genes, to infer interspecies relationships and the position of Camptoloma within Scrophulariaceae by using supermatrix and multispecies-coalescent approaches. Lineage divergence times and ancestral ranges were inferred with Bayesian Markov chain Monte Carlo (MCMC) approaches. The population history was estimated with phylogeographic coalescent methods.

RESULTS

Camptoloma rotundifolium, restricted to Southern Africa, was shown to be a sister species to the disjunct clade formed by C. canariense, endemic to the Canary Islands, and C. lyperiiflorum, distributed in the Horn of Africa-Southern Arabia. Camptoloma was inferred to be sister to the mostly South African tribes Teedieae and Buddlejeae. Stem divergence was dated in the Late Miocene, while the origin of the extant disjunction was inferred as Early Pliocene.

CONCLUSIONS

The current disjunct distribution of Camptoloma across Africa was likely the result of fragmentation and extinction and/or population bottlenecking events associated with historical aridification cycles during the Neogene; the pattern of species divergence, from south to north, is consistent with the "climatic refugia" Rand Flora hypothesis.

Reconstruction of the evolutionary biogeography reveal the origins and diversification of oysters (Bivalvia: Ostreidae)

Oysters (Bivalvia: Ostreidae Rafinesque, 1815) live in the intertidal and shallow subtidal areas worldwide. Despite their long evolutionary histories, abundant fossil records, global distribution, and ecological significance, a systematic time-dependent biogeographical analysis of this family is still lacking. Using combined mitochondrial (COI and 16S rRNA) and nuclear (18S rRNA, 28S rRNA, H3 and ITS2) gene makers for 80% (70/88) of the recognized extant Ostreidae, we reconstructed the global phylogenetic and biogeographical relationships throughout the evolutionary history of oysters. The result provided a holistic view of the origin, migration and dispersal patterns of Ostreidae. The phylogenetic results and fossil evidence indicated that Ostreidae originated from the circum-Arctic region in the Early Jurassic. The widening of the Atlantic Ocean and changes in the Tethys Ocean further facilitated their subsequent diversification during the Cretaceous and the Palaeogene periods. In particular, Crassostrea and Saccostrea exhibited relatively low dispersal abilities and their major diversifications were consistent with the tectonic events. Environmental adaptations and reproductive patterns, therefore, should play key roles in the formation of oyster distribution patterners, rather than the dispersal ability of their planktonic larvae. The diversity dynamics inferred by standard phylogenetic are consistent with the fossil record, however, further systematic classification, especially for fossil genus Ostrea, would enhance our understanding on extant and fossil oysters. The present study of the historical biogeography of oysters provides new insights into the evolution and speciation of oysters. Our findings also provide a foundation for the assessment of evolutionary patterns and ecological processes in intertidal and inshore life.

Phytogeographic history of the Tea family inferred through high-resolution phylogeny and fossils

The tea family (Theaceae) has a highly unusual amphi-Pacific disjunct distribution: most extant species in the family are restricted to subtropical evergreen broadleaf forests in East Asia, while a handful of species occur exclusively in the subtropical and tropical Americas. Here we used an approach that integrates the rich fossil evidence of this group with phylogenies in biogeographic analysis to study the processes behind this distribution pattern. We first combined genome-skimming sequencing with existing molecular data to build a robust species-level phylogeny for c.140 Theaceae species, resolving most important unclarified relationships. We then developed an empirical Bayesian method to incorporate distribution evidence from fossil specimens into historical biogeographic analyses and used this method to account for the spatiotemporal history of Theaceae fossils. We compared our method with an alternative Bayesian approach and show that it provides consistent results while significantly reduces computational demands which allows analyses of much larger datasets. Our analyses revealed a circumboreal distribution of the family from the early Cenozoic to the Miocene and inferred repeated expansions and retractions of the modelled distribution in the Northern Hemisphere, suggesting that the current Theaceae distribution could be the remnant of a larger continuous distribution associated with the boreotropical forest that has been hypothesized to occupy most of the northern latitudes in the early Cenozoic. These results contradict with studies that only considered current species distributions and showcase the necessity of integrating fossil and molecular data in phylogeny-based parametric biogeographic models to improve the reliability of inferred biogeographical events.

Incorporating palaeogeography into ancestral area estimation can explain the disjunct distribution of land snails in Macaronesia and the Balearic Islands (Helicidae: Allognathini)

The systematics and biogeographical history of the Eastern Mediterranean and Macaronesian land snail tribe Allognathini (Helicidae: Helicinae) is investigated based on mitochondrial and nuclear DNA sequence data. Our molecular phylogenetic analyses indicate that the genus-group systematics of the tribe needs to be revised. We show for the first time that the narrow-range endemics Lampadia and Idiomela from the Madeira Archipelago belong to Allognathini and represent together the sister group of the diverse Canary Island Hemicycla radiation. We therefore suggest synonymising Lampadiini with Allognathini. Sister to these Macaronesian genera was the Balearic Island Allognathus radiation. Pseudotachea was not recovered as a monophyletic group and the two currently recognised species clustered in Iberus. Similarly, Adiverticula was not recovered as a monophyletic group and clustered in Hemicycla. We therefore suggest synonymising Pseudotachea with Iberus and Adiverticula with Hemicycla. The six genera in Allognathini, which we distinguish here (Cepaea, Iberus, Allognathus, Hemicycla, Idiomela and Lampadia), originated in Western to South-western Europe according to our ancestral area estimation and the fossil record. The disjunct distribution of the Balearic Islands and Macaronesian sister clades and the mainly Iberian Iberus clade that separated earlier can be explained by the separation of the Betic-Rif System from the Iberian Peninsula during the late Oligocene to early Miocene, along with independent Miocene dispersals to the Balearic Islands and Macaronesia from the Iberian Peninsula, where the ancestral lineage became extinct.

Mitogenomes Reveal Multiple Colonization of Mountains by Rattus in Sundaland

Tropical mountains are cradles of biodiversity and endemism. Sundaland, tropical Southeast Asia, hosts 3 species of Rattus endemic to elevations above 2000 m with an apparent convergence in external morphology: Rattus korinchi and R. hoogerwerfi from Sumatra, and R. baluensis from Borneo. A fourth one, R. tiomanicus, is restricted to lowland elevations across the whole region. The origins of these endemics are little known due to the absence of a robust phylogenetic framework. We use complete mitochondrial genomes from the 3 high altitude Rattus, and several related species to determine their relationships, date divergences, reconstruct their history of colonization, and test for selection on the mitochondrial DNA. We show that mountain colonization happened independently in Borneo (<390 Kya) and Sumatra (~1.38 Mya), likely from lowland lineages. The origin of the Bornean endemic R. baluensis is very recent and its genetic diversity is nested within the diversity of R. tiomanicus. We found weak evidence of positive selection in the high-elevation lineages and attributed the greater nonsynonymous mutations on these branches (specially R. baluensis) to lesser purifying selection having acted on the terminal branches in the phylogeny.

Role of diversification rates and evolutionary history as a driver of plant naturalization success

Human introductions of species beyond their natural ranges and their subsequent establishment are defining features of global environmental change. However, naturalized plants are not uniformly distributed across phylogenetic lineages, with some families contributing disproportionately more to the global alien species pool than others. Additionally, lineages differ in diversification rates, and high diversification rates have been associated with characteristics that increase species naturalization success. Here, we investigate the role of diversification rates in explaining the naturalization success of angiosperm plant families. We use five global data sets that include native and alien plant species distribution, horticultural use of plants, and a time-calibrated angiosperm phylogeny. Using phylogenetic generalized linear mixed models, we analysed the effect of diversification rate, different geographical range measures, and horticultural use on the naturalization success of plant families. We show that a family's naturalization success is positively associated with its evolutionary history, native range size, and economic use. Investigating interactive effects of these predictors shows that native range size and geographic distribution additionally affect naturalization success. High diversification rates and large ranges increase naturalization success, especially of temperate families. We suggest this may result from lower ecological specialization in temperate families with large ranges, compared with tropical families with smaller ranges.

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.

An escape-to-radiate model for explaining the high plant diversity and endemism in campos rupestres†

With extraordinary levels of plant diversity and endemism, the Brazilian campos rupestres across the Espinhaço Range have a species/area ratio 40 times higher than the lowland Amazon. Although diversification drivers in campos rupestres remain a matter of debate, the Pleistocene refugium hypothesis (PRH) is often adopted as the most plausible explanation for their high diversity. The PRH has two main postulates: highland interglacial refugia and a species pump mechanism catalysed by climatic changes. We critically assessed studies on campos rupestres diversification at different evolutionary levels and conclude that most of them are affected by sampling biases, unrealistic assumptions or inaccurate results that do not support the PRH. By modelling the palaeo-range of campos rupestres based on the distribution of 1123 species of vascular plants endemic to the Espinhaço Range and using climate and edaphic variables, we projected a virtually constant suitable area for campos rupestres across the last glacial cycle. We challenge the great importance placed on Pleistocene climatic oscillations in campos rupestres plant diversification and offer an alternative explanation named escape-to-radiate model, which emphasizes niche shifts. Under this biogeographic model of diversification, the long-term fragmentation of campos rupestres combined with recurrent extinctions after genetic drift and sporadic events of adaptive radiation may provide an explanation for the current diversity and endemism in the Espinhaço Range. We conclude that long-term diversification dynamics in campos rupestres are mainly driven by selection, while most endemic diversity is ephemeral, extremely fragile and mainly driven by drift.

Ancient tropical extinctions at high latitudes contributed to the latitudinal diversity gradient

Global biodiversity currently peaks at the equator and decreases toward the poles. Growing fossil evidence suggest this hump-shaped latitudinal diversity gradient (LDG) has not been persistent through time, with similar diversity across latitudes flattening out the LDG during past greenhouse periods. However, when and how diversity declined at high latitudes to generate the modern LDG remains an open question. Although diversity-loss scenarios have been proposed, they remain mostly undemonstrated. We outline the "asymmetric gradient of extinction and dispersal" framework that contextualizes previous ideas behind the LDG under a time-variable scenario. Using phylogenies and fossils of Testudines, Crocodilia, and Lepidosauria, we find that the hump-shaped LDG could be explained by (1) disproportionate extinctions of high-latitude tropical-adapted clades when climate transitioned from greenhouse to icehouse, and (2) equator-ward biotic dispersals tracking their climatic preferences when tropical biomes became restricted to the equator. Conversely, equivalent diversification rates across latitudes can account for the formation of an ancient flat LDG. The inclusion of fossils in macroevolutionary studies allows revealing time-dependent extinction rates hardly detectable from phylogenies only. This study underscores that the prevailing evolutionary processes generating the LDG during greenhouses differed from those operating during icehouses.

Fossil and phylogenetic analyses reveal recurrent periods of diversification and extinction in dictyopteran insects

Variations of speciation and extinction rates determine the fate of clades through time. Periods of high diversification and extinction (possibly mass-extinction events) can punctuate the evolutionary history of various clades, but they remain loosely defined for many biological groups, especially nonmarine invertebrates like insects. Here, we examine whether the cockroaches, mantises and termites (altogether included in Dictyoptera) have experienced episodic pulses of speciation or extinction and how these pulses may be associated with environmental fluctuations or mass extinctions. We relied on molecular phylogeny and fossil data to shed light on the times and rates at which dictyopterans diversified. The diversification of Dictyoptera has alternated between (i) periods of high diversification in the late Carboniferous, Early-Middle Triassic, Early Cretaceous and middle Palaeogene, and (ii) periods of high extinction rates particularly at the Permian-Triassic boundary, but not necessarily correlated with the major global biodiversity crises as in the mid-Cretaceous. This study advocates the importance of analyzing, when possible, both molecular phylogeny and fossil data to unveil diversification and extinction periods for a given group. The causes and consequences of extinction must be studied beyond mass-extinction events alone to gain a broader understanding of how clades wax and wane.

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 power of Bayesian birth-death skyline models to detect mass extinction events from phylogenies with only extant taxa

Mass extinction events (MEEs), defined as significant losses of species diversity in significantly short time periods, have attracted the attention of biologists because of their link to major environmental change. MEEs have traditionally been studied through the fossil record, but the development of birth-death models has made it possible to detect their signature based on extant-taxa phylogenies. Most birth-death models consider MEEs as instantaneous events where a high proportion of species are simultaneously removed from the tree ("single pulse" approach), in contrast to the paleontological record, where MEEs have a time duration. Here, we explore the power of a Bayesian Birth-Death Skyline (BDSKY) model to detect the signature of MEEs through changes in extinction rates under a "time-slice" approach. In this approach, MEEs are time intervals where the extinction rate is greater than the speciation rate. Results showed BDSKY can detect and locate MEEs but that precision and accuracy depend on the phylogeny's size and MEE intensity. Comparisons of BDSKY with the single-pulse Bayesian model, CoMET, showed a similar frequency of Type II error and neither model exhibited Type I error. However, while CoMET performed better in detecting and locating MEEs for smaller phylogenies, BDSKY showed higher accuracy in estimating extinction and speciation rates.

Complex interplay of ancient vicariance and recent patterns of geographical speciation in north-western North American temperate rainforests explains the phylogeny of jumping slugs (Hemphillia spp.)

The history of the currently disjunct temperate rainforests of the Pacific Northwest of North America has shaped the evolution and diversity of endemics. This study focuses on how geological and climatic perturbations have driven speciation in the area by isolating lineages. We investigated the phylogenetic relationships and historical biogeography of the endemic jumping slugs (genus Hemphillia) using a multi-locus phylogeny. We evaluated the spatial distribution and divergence times of major lineages, generated ancestral area probabilities and inferred the biogeographical history of the genus. Our study revealed eight genetic lineages that formed three clades: one clade consisting of two Coast/Cascade lineages, and two reciprocally monophyletic clades that each contain a Coast/Cascade and two Rocky Mountains taxa. The results of the biogeographical analysis suggest that the ancestral range of the genus occupied Coast/Cascade habitats and then spread across into Northern Rocky Mountain interior habitats with subsequent fragmentations isolating coastal and inland lineages. Finally, there have been more recent speciation events among three lineage pairs that have shaped shallow structures of all clades. We add to our knowledge of the biogeographical history of the region in that we discovered diversification and speciation events that have occurred in ways more complex than previously thought.

Morphological Innovations and Vast Extensions of Mountain Habitats Triggered Rapid Diversification Within the Species-Rich Irano-Turanian Genus Acantholimon (Plumbaginaceae)

The Irano-Turanian floristic region spans a topographically complex and climatically continental territory, which has served as a source of xerophytic taxa for neighboring regions and is represented by a high percent of endemics. Yet, a comprehensive picture of the abiotic and biotic factors that have driven diversification within this biota remains to be established due to the scarcity of phylogenetic studies. Acantholimon is an important component of the subalpine steppe flora of the Irano-Turanian region, containing c. 200 cushion-forming sub-shrubby pungent-leaved species. Our recent molecular phylogenetic study has led to enlarging the circumscription of this genus to include eight mono- or oligospecific genera lacking the characteristic life-form and leaves. Using the same molecular phylogeny, here we investigate the tempo and mode of diversification as well as the biogeographic patterns in this genus, to test the hypothesis that a combination of key morphological innovations and abiotic factors is behind Acantholimon high species diversity. Molecular dating analysis indicates that Acantholimon s.l. started to diversify between the Late Miocene and the Pliocene and the biogeographic analysis points to an Eastern Iran-Afghanistan origin. Macroevolutionary models support the hypothesis that the high diversity of the genus is explained by accelerated diversification rates in two clades associated with the appearance of morphological key innovations such as a cushion life-form and pungent leaves; this would have favored the colonization of water-stressed, substrate-poor mountainous habitats along the newly uplifted IT mountains during the Mio-Pliocene. Given the apparent similarity of mountain habitats for most species of Acantholimon, we hypothesize that its current high species diversity responds to a scenario of non-adaptive radiation fueled by allopatric speciation rather than evolutionary radiation driven by ecological opportunity. Similar scenarios might underlie the high diversity of other speciose genera in the topographically complex Irano-Turanian landscape, though this remains to be tested with fine-grained distribution and climatic data.

Can the pyruvate: ferredoxin oxidoreductase (PFOR) gene be used as an additional marker to discriminate among Blastocystis strains or subtypes?

BACKGROUND

Blastocystis spp. are the most prevalent intestinal eukaryotes identified in humans, with at least 17 genetic subtypes (ST) based on genes coding for the small-subunit ribosomal RNA (18S). It has been argued that the 18S gene should not be the marker of choice to discriminate between STs of these strains because this marker exhibits high intra-genomic polymorphism. By contrast, pyruvate:ferredoxin oxidoreductase (PFOR) is a relevant enzyme involved in the core energy metabolism of many anaerobic microorganisms such as Blastocystis, which, in other protozoa, shows more polymorphisms than the 18S gene and thus may offer finer discrimination when trying to identify Blastocystis ST. Therefore, the objective of the present study was to assess the suitability of the PFOR gene as an additional marker to discriminate among Blastocystis strains or subtypes from symptomatic carrier children.

METHODS

Faecal samples from 192 children with gastrointestinal symptoms from the State of Mexico were submitted for coprological study. Twenty-one of these samples were positive only for Blastocystis spp.; these samples were analysed by PCR sequencing of regions of the 18S and PFOR genes. The amplicons were purified and sequenced; afterwards, both markers were assessed for genetic diversity.

RESULTS

The 18S analysis showed the following frequencies of Blastocystis subtypes: ST3 = 43%; ST1 = 38%; ST2 = 14%; and ST7 = 5%. Additionally, using subtype-specific primer sets, two samples showed mixed Blastocystis ST1 and ST2 infection. For PFOR, Bayesian inference revealed the presence of three clades (I-III); two of them grouped different ST samples, and one grouped six samples of ST3 (III). Nucleotide diversity (π) and haplotype polymorphism (θ) for the 18S analysis were similar for ST1 and ST2 (π = ~0.025 and θ = ~0.036); remarkably, ST3 showed almost 10-fold lower values. For PFOR, a similar trend was found: clade I and II had π = ~0.05 and θ = ~0.05, whereas for clade III, the values were almost 6-fold lower.

CONCLUSIONS

Although the fragment of the PFOR gene analysed in the present study did not allow discrimination between Blastocystis STs, this marker grouped the samples in three clades with strengthened support, suggesting that PFOR may be under different selective pressures and evolutionary histories than the 18S gene. Interestingly, the ST3 sequences showed lower variability with probable purifying selection in both markers, meaning that evolutionary forces drive differential processes among Blastocystis STs.

Conceptual and empirical advances in Neotropical biodiversity research

The unparalleled biodiversity found in the American tropics (the Neotropics) has attracted the attention of naturalists for centuries. Despite major advances in recent years in our understanding of the origin and diversification of many Neotropical taxa and biotic regions, many questions remain to be answered. Additional biological and geological data are still needed, as well as methodological advances that are capable of bridging these research fields. In this review, aimed primarily at advanced students and early-career scientists, we introduce the concept of "trans-disciplinary biogeography," which refers to the integration of data from multiple areas of research in biology (e.g., community ecology, phylogeography, systematics, historical biogeography) and Earth and the physical sciences (e.g., geology, climatology, palaeontology), as a means to reconstruct the giant puzzle of Neotropical biodiversity and evolution in space and time. We caution against extrapolating results derived from the study of one or a few taxa to convey general scenarios of Neotropical evolution and landscape formation. We urge more coordination and integration of data and ideas among disciplines, transcending their traditional boundaries, as a basis for advancing tomorrow's ground-breaking research. Our review highlights the great opportunities for studying the Neotropical biota to understand the evolution of life.

Detecting local diversity-dependence in diversification

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

Why are there so many sedges? Sumatroscirpeae, a missing piece in the evolutionary puzzle of the giant genus Carex (Cyperaceae)

For over a century, the origins and mechanisms underlying the diversification of the enormous temperate genus Carex (>2100 species; Cariceae, Cyperaceae) have remained largely speculative. Characteristics such as its diverse ecology, varied biogeography, and intriguing cytology have made Carex a powerful model for studying plant evolution, but its uncertain sister-group relationships hinder its use in studies that depend on accurate ancestral state estimates and biogeographic inferences. To identify the sister to Carex, we estimated the phylogeny of all genera in the Cariceae-Dulichieae-Scirpeae clade (CDS) using three plastid and two nuclear ribosomal markers. Ancestral state reconstructions of key characters were made, and a time-calibrated tree estimated. Carex is strongly supported as sister to the rare East Asian Sumatroscirpus, sole genus of a new tribe, Sumatroscirpeae trib. nov. Believed to be unique to Carex, the perigynium (prophyllar bract enclosing a flower) is in fact a synapomorphy shared with this small tribe (∼4 species) that appeared 36 Mya. We thus suggest the initial key innovation in the remarkable diversification of Carex is not the perigynium, but could be the release of mechanical constraints on perigynia through spikelet truncation, resulting in novel adaptive morphologies. Monoecy, chromosomal change, and rapid inflorescence development enabling phenological isolation may also be involved. The tiny tribe Sumatroscirpeae will provide unprecedented insights into the inflorescence homology, evolution, diversification, and biogeographic history of its sister-group Carex, one of the world's most diverse plant lineages.