Persistence of within-species lineages: a neglected control of speciation rates

Genomic, morphological and physiological data support fast ecotypic differentiation and incipient speciation in an alpine diving beetle

An intricate interplay between evolutionary and demographic processes has frequently resulted in complex patterns of genetic and phenotypic diversity in alpine lineages, posing serious challenges to species delimitation and biodiversity conservation planning. Here we integrate genomic data, geometric morphometric analyses and thermal tolerance experiments to explore the role of Pleistocene climatic changes and adaptation to alpine environments on patterns of genomic and phenotypic variation in diving beetles from the taxonomically complex Agabus bipustulatus species group. Genetic structure and phylogenomic analyses revealed the presence of three geographically cohesive lineages, two representing trans-Palearctic and Iberian populations of the elevation-generalist A. bipustulatus and another corresponding to the strictly-alpine A. nevadensis, a narrow-range endemic taxon from the Sierra Nevada mountain range in southeastern Iberia. The best-supported model of lineage divergence, along with the existence of pervasive genetic introgression and admixture in secondary contact zones, is consistent with a scenario of population isolation and connectivity linked to Quaternary climatic oscillations. Our results suggest that A. nevadensis is an alpine ecotype of A. bipustulatus, whose genotypic, morphological and physiological differentiation likely resulted from an interplay between population isolation and local altitudinal adaptation. Remarkably, within the Iberian Peninsula, such ecotypic differentiation is unique to Sierra Nevada populations and has not been replicated in other alpine populations of A. bipustulatus. Collectively, our study supports fast ecotypic differentiation and incipient speciation processes within the study complex and points to Pleistocene glaciations and local adaptation along elevational gradients as key drivers of biodiversity generation in alpine environments.

Drivers of genomic differentiation landscapes in populations of disparate ecological and geographical settings within mainland Apis cerana

Elucidating the evolutionary processes that drive population divergence can enhance our understanding of the early stages of speciation and inform conservation management decisions. The honeybee Apis cerana displays extensive population divergence, providing an informative natural system for exploring these processes. The mainland lineage A. cerana includes several peripheral subspecies with disparate ecological and geographical settings radiated from a central ancestor. Under this evolutionary framework, we can explore the patterns of genome differentiation and the evolutionary models that explain them. We can also elucidate the contribution of non-genomic spatiotemporal mechanisms (extrinsic features) and genomic mechanisms (intrinsic features) that influence these genomic differentiation landscapes. Based on 293 whole genomes, a small part of the genome is highly differentiated between central-peripheral subspecies pairs, while low and partial parallelism partly reflects idiosyncratic responses to environmental differences. Combined elements of recurrent selection and speciation-with-gene-flow models generate the heterogeneous genome landscapes. These elements weight differently between central-island and other central-peripheral subspecies pairs, influenced by glacial cycles superimposed on different geomorphologies. Although local recombination rates exert a significant influence on patterns of genomic differentiation, it is unlikely that low-recombination rates regions were generated by structural variation. In conclusion, complex factors including geographical isolation, divergent ecological selection and non-uniform genome features have acted concertedly in the evolution of reproductive barriers that could reduce gene flow in part of the genome and facilitate the persistence of distinct populations within mainland lineage of A. cerana.

Quantitative evaluation of the drivers of species richness in a Mediterranean ecosystem (Cape, South Africa)

BACKGROUND AND AIMS

Mediterranean ecosystems have a high vascular plant species richness (SR) relative to their surface area. This SR, representing the balance between speciation and extinction, has been attributed to multiple mechanisms that result in both high rates of speciation and/or low rates of extinction. An abiding question is, however, what is special about Mediterranean ecosystems that enables this high SR? Apart from the long-term climatic stability of the region, SR has also been related to resource availability, the many individuals hypothesis, resource spatial heterogeneity, temporal heterogeneity and biotic feedbacks.

METHODS

Spatial patterns of species richness were related to climatic, edaphic and biotic variables and to spatial variability within the Greater Cape Floristic Region (GCFR) of South Africa. Boosted regression tree models were used to explore the strength of relationships between SR and environmental predictors related to each hypothesized mechanism.

KEY RESULTS

Water availability (i.e. precipitation) was a stronger predictor of SR than potential evapotranspiration or temperature. Scarcity of nutrients was also related to SR. There was no indication that SR was related to the density of individuals and only temporal heterogeneity induced by fire was related to SR. Spatial heterogeneities of climatic, edaphic and biotic variables were strongly associated with SR. Biotic interactions remain difficult to assess, although we have some evidence for a putative role in regulating SR.

CONCLUSIONS

While the lack of ecosystem-resetting disturbances (e.g. glaciation) is undoubtedly a key requirement for high species accumulation, predictably, no one explanation holds the key to understanding SR. In the GCFR high SR is the product of a combination of adequate water, nutrient scarcity, spatial and temporal heterogeneity, and possibly biotic feedbacks.

An Elevational Phylogeographic Diversity Gradient in Neotropical Birds Is Decoupled from Speciation Rates

AbstractA key question about macroevolutionary speciation rates is whether they are controlled by microevolutionary processes operating at the population level. For example, does spatial variation in population genetic differentiation underlie geographical gradients in speciation rates? Previous work suggests that speciation rates increase with elevation in Neotropical birds, but underlying population-level gradients remain unexplored. Here, we characterize elevational phylogeographic diversity between montane and lowland birds in the megadiverse Andes-Amazonian system and assess its relationship to speciation rates to evaluate the link between population-level differentiation and species-level diversification. We aggregated and georeferenced nearly 7,000 mitochondrial DNA sequences across 103 species or species complexes in the Andes and Amazonia and used these sequences to describe phylogeographic differentiation across both regions. Our results show increased levels of both discrete and continuous metrics of population structure in the Andean mountains compared with the Amazonian lowlands. However, higher levels of population differentiation do not predict higher rates of speciation in our dataset. Multiple potential factors may lead to our observed decoupling of initial population divergence and speciation rates, including the ephemerality of incipient species and the multifaceted nature of the speciation process, as well as methodological challenges associated with estimating rates of population differentiation and speciation.

Rise, fall and hope for the Sicilian endemic plant Muscari gussonei: A story of survival in the face of narrow germination optimum, climate changes, desertification and habitat fragmentation

Muscari gussonei is an endangered endemic plant growing on fragmented Mediterranean coastal dunes. This study focused on the germination performance of M. gussonei at two fixed temperatures, 10 and 15 °C, and at an alternating one, 10/20 °C, and on the multi-temporal trends of temperature and rainfall during 1931-2020, as well as on the patterns of desertification and land-cover changes over the last 60 years. High and similar germinability was found for different populations of M. gussonei, in particular, the final germination percentage (FGP) was ≥95 % for the three treatments. The general pattern was the lower the temperature the higher and faster the germination. However, germination speed varied significantly among populations. This intraspecific variability of germination behavior may suggest a certain level of ecophysiological plasticity in M. gussonei, thus raising hopes on the capacity of M. gussonei to respond better to the ongoing severe environmental changes. In the period 1931-2020, indeed, the average temperature rose by 1.5 °C, from 16.8 to 18.3 °C, which is equivalent to the enormous increase of 0.17 °C per decade. Similarly, the average rainfall declined by 100 mm, from 600 to 500 mm. Another serious stressor was desertification, which affects >90 % of the distributional area of M. gussonei. A further factor of ecological degradation is a considerably altered landscape, where the agricultural component accounts for c. 85 %, whereas natural and seminatural areas were only c. 10 %. Increasing temperature and dryness will inevitably reduce the germinability of M. gussonei, characterized by a narrow germination optimum of 10-15 °C. The future of M. gussonei looks even more dramatic if we consider its small and scattered populations distributed in an agricultural matrix affected by high levels of desertification. Only multivariate information at different space-time scales can provide an exhaustive picture for implementing effective conservation strategies.

Discerning structure versus speciation in phylogeographic analysis of Seepage Salamanders (Desmognathus aeneus) using demography, environment, geography, and phenotype

Numerous mechanisms can drive speciation, including isolation by adaptation, distance, and environment. These forces can promote genetic and phenotypic differentiation of local populations, the formation of phylogeographic lineages, and ultimately, completed speciation. However, conceptually similar mechanisms may also result in stabilizing rather than diversifying selection, leading to lineage integration and the long-term persistence of population structure within genetically cohesive species. Processes that drive the formation and maintenance of geographic genetic diversity while facilitating high rates of migration and limiting phenotypic differentiation may thereby result in population genetic structure that is not accompanied by reproductive isolation. We suggest that this framework can be applied more broadly to address the classic dilemma of "structure" versus "species" when evaluating phylogeographic diversity, unifying population genetics, species delimitation, and the underlying study of speciation. We demonstrate one such instance in the Seepage Salamander (Desmognathus aeneus) from the southeastern United States. Recent studies estimated up to 6.3% mitochondrial divergence and four phylogenomic lineages with broad admixture across geographic hybrid zones, which could potentially represent distinct species supported by our species-delimitation analyses. However, while limited dispersal promotes substantial isolation by distance, microhabitat specificity appears to yield stabilizing selection on a single, uniform, ecologically mediated phenotype. As a result, climatic cycles promote recurrent contact between lineages and repeated instances of high migration through time. Subsequent hybridization is apparently not counteracted by adaptive differentiation limiting introgression, leaving a single unified species with deeply divergent phylogeographic lineages that nonetheless do not appear to represent incipient species.

Genetic signatures of lineage fusion closely resemble population decline

Accurate interpretation of the genetic signatures of past demographic events is crucial for reconstructing evolutionary history. Lineage fusion (complete merging, resulting in a single panmictic population) is a special case of secondary contact that is seldom considered. Here, the circumstances under which lineage fusion can be distinguished from population size constancy, growth, bottleneck, and decline were investigated. Multi-locus haplotype data were simulated under models of lineage fusion with different divergence versus sampling lag times (D:L ratios). These pseudo-observed datasets also differed in their allocation of a fixed amount of sequencing resources (number of sampled alleles, haplotype length, number of loci). Distinguishability of lineage fusion versus each of 10 untrue non-fusion scenarios was quantified based on six summary statistics (neutrality tests). Some datasets were also analyzed using extended Bayesian skyline plots. Results showed that signatures of lineage fusion very closely resemble those of decline-high distinguishability was generally limited to the most favorable scenario (D:L = 9), using the most sensitive summary statistics (F S and Z nS), coupled with the optimal sequencing resource allocation (maximizing number of loci). Also, extended Bayesian skyline plots often erroneously inferred population decline. Awareness of the potential for lineage fusion to carry the hallmarks of population decline is critical.

Conceptual and empirical bridges between micro- and macroevolution

Explaining broad molecular, phenotypic and species biodiversity patterns necessitates a unifying framework spanning multiple evolutionary scales. Here we argue that although substantial effort has been made to reconcile microevolution and macroevolution, much work remains to identify the links between biological processes at play. We highlight four major questions of evolutionary biology whose solutions require conceptual bridges between micro and macroevolution. We review potential avenues for future research to establish how mechanisms at one scale (drift, mutation, migration, selection) translate to processes at the other scale (speciation, extinction, biogeographic dispersal) and vice versa. We propose ways in which current comparative methods to infer molecular evolution, phenotypic evolution and species diversification could be improved to specifically address these questions. We conclude that researchers are in a better position than ever before to build a synthesis to understand how microevolutionary dynamics unfold over millions of years.

Integrative characterization of genetic and phenotypic differentiation in an ant species complex with strong hierarchical population structure and low dispersal abilities

Low dispersal, occurrence of asexual reproduction and geographic discontinuity increase genetic differentiation between populations, which ultimately can lead to speciation. In this work, we used a multidisciplinary framework to characterize the genetic and phenotypic differentiation between and within two cryptic ant species with restricted dispersal, Cataglyphis cursor and C. piliscapa and used behavioral experiments to test for reproductive isolation. Their distribution is segregated by the Rhône River and they have been traditionally distinguished only by hair numbers, although a statistical assessment is still lacking. We found strong genetic (microsatellites, nuclear and mitochondrial sequences), morphological (number of hairs, tibia length, male genitalia) and chemical (cuticular hydrocarbons) differentiation not only between species but also among localities within species. However, inter-specific differentiation was slightly higher than intra-specific differentiation for most markers. Overall, this pattern could either reflect reproductive isolation or could result from a longer period of geographic isolation between species than among localities within species without necessarily involving reproductive isolation. Interestingly, our behavioral experiments showed an absence of mating between species associated to a higher aggressiveness of workers towards heterospecific males. This suggests that sexual selection may, at least partially, fuel reproductive isolation. We also showed that cuticular hydrocarbons, mtDNA sequences and number of hairs provide reliable criteria allowing species discrimination. Overall, this species complex offers a case study to further investigate varying stages of a speciation continuum by estimating reproductive isolation between pairs of localities varying by their level of genetic differentiation.

Evolutionary time and species diversity in aquatic ecosystems worldwide

The latitudinal diversity gradient (LDG) is frequently described as the most dramatic biodiversity pattern on Earth, yet ecologists and biogeographers have failed to reach consensus on its primary cause. A key problem in explaining the LDG involves collinearity between multiple factors that are predicted to affect species richness in the same direction. In terrestrial systems, energy input, geographic area, and evolutionary time for species accumulation tend to covary positively with species richness at the largest spatial scales, such that their individual contributions to the LDG are confounded in global analyses. I review three diversity patterns from marine and freshwater systems that break this collinearity and which may thus provide stronger tests of the influence of time on global richness gradients. Specifically, I contrast biodiversity patterns along oceanic depth gradients, in geologically young versus ancient lakes, and in the north versus south polar marine biomes. I focus primarily on fishes due to greater data availability but synthesize patterns for invertebrates where possible. I find that regional-to-global species richness generally declines with depth in the oceans, despite the great age and stability of the deep-sea biome. Geologically ancient lakes generally do not contain more species than young lakes, and the Antarctic marine biome is not appreciably more species rich than the much younger Arctic marine biome. However, endemism is consistently higher in older systems. Patterns for invertebrate groups are less clear than for fishes and reflect a critical need for primary biodiversity data. In summary, the available data suggest that species richness is either decoupled from or only weakly related to the amount of time for diversification. These results suggest that energy, productivity, or geographic area are the primary drivers of large-scale diversity gradients. To the extent that marine and terrestrial diversity gradients result from similar processes, these examples provide evidence against a primary role for evolutionary time as the cause of the LDG.

The phoenix hypothesis of speciation

Genetic divergence among allopatric populations builds reproductive isolation over time. This process is accelerated when populations face a changing environment that allows large-effect mutational differences to accumulate, but abrupt change also places populations at risk of extinction. Here we use simulations of Fisher's geometric model with explicit population dynamics to explore the genetic changes that occur in the face of environmental changes. Because evolutionary rescue leads to the fixation of mutations whose phenotypic effects are larger on average compared with populations not at risk of extinction, these mutations are thus more likely to lead to reproductive isolation. We refer to the formation of new species from the ashes of populations in decline as the phoenix hypothesis of speciation. The phoenix hypothesis predicts more substantial hybrid fitness breakdown among populations surviving a higher extinction risk. The hypothesis was supported when many loci underlie adaptation. With only a small number of potential rescue mutations, however, mutations that fixed in different populations were more likely to be identical, with such parallel changes reducing isolation. Consequently, reproductive isolation builds fastest in populations subject to an intermediate extinction risk, given a limited number of mutations available for adaptation.

Genomic evidence of speciation by fusion in a recent radiation of grasshoppers

Postdivergence gene flow can trigger a number of creative evolutionary outcomes, ranging from the transfer of beneficial alleles across species boundaries (i.e., adaptive introgression) to the formation of new species (i.e., hybrid speciation). Although neutral and adaptive introgression has been broadly documented in nature, hybrid speciation is assumed to be rare and the evolutionary and ecological context facilitating this phenomenon still remains controversial. Through combining genomic and phenotypic data, we evaluate the hypothesis that the dual feeding regime (based on both scrub legumes and gramineous herbs) of the taxonomically controversial grasshopper Chorthippus saulcyi algoaldensis resulted from hybridization between the sister taxa C. binotatus (that exclusively feeds on scrub legumes) and C. saulcyi (that only feeds on gramineous herbs). Genetic clustering analyses and inferences from coalescent-based demographic simulations confirm that C. s. algoaldensis represents an independently evolving lineage and support the ancient hybrid origin of this taxon (about 1.4 Ma), which sheds light on its uncertain phylogenetic position and might explain its broader trophic niche. We propose a Pleistocene hybrid speciation model where range shifts resulting from climatic oscillations can promote the formation of hybrid swarms and facilitate their long-term persistence through geographic isolation from parental forms in topographically complex landscapes.

Soil seed bank responses to edge effects in temperate European forests

Aim

The amount of forest edges is increasing globally due to forest fragmentation and land-use changes. However, edge effects on the soil seed bank of temperate forests are still poorly understood. Here, we assessed edge effects at contrasting spatial scales across Europe and quantified the extent to which edges can preserve the seeds of forest specialist plants.

Location

Temperate European deciduous forests along a 2,300-km latitudinal gradient.

Time period

2018-2021.

Major taxa studied

Vascular plants.

Methods

Through a greenhouse germination experiment, we studied how edge effects alter the density, diversity, composition and functionality of forest soil seed banks in 90 plots along different latitudes, elevations and forest management types. We also assessed which environmental conditions drive the seed bank responses at the forest edge versus interior and looked at the relationship between the seed bank and the herb layer species richness.

Results

Overall, 10,108 seedlings of 250 species emerged from the soil seed bank. Seed density and species richness of generalists (species not only associated with forests) were higher at edges compared to interiors, with a negative influence of C : N ratio and litter quality. Conversely, forest specialist species richness did not decline from the interior to the edge. Also, edges were compositionally, but not functionally, different from interiors. The correlation between the seed bank and the herb layer species richness was positive and affected by microclimate.

Main conclusions

Our results underpin how edge effects shape species diversity and composition of soil seed banks in ancient forests, especially increasing the proportion of generalist species and thus potentially favouring a shift in community composition. However, the presence of many forest specialists suggests that soil seed banks still play a key role in understorey species persistence and could support the resilience of our fragmented forests.

Ehrlich and Raven escape and radiate coevolution hypothesis at different levels of organization: Past and future perspectives

The classic paper by Ehrlich and Raven on coevolution will soon be 60 years old. Although they were not the first to develop the idea of coevolution, their thought-provoking paper certainly popularized this idea and inspired several generations of scientists interested in coevolution. Here, we describe some of their main contributions, quantitatively measure the impact of their seminal paper on different fields of research, and discuss how ideas related to their original paper might push the study of coevolution forward. To guide our discussion, we explore their original hypothesis into three research fields that are associated with distinct scales/levels of organization: (1) the genetic mechanisms underlying coevolutionary interactions; (2) the potential association between coevolutionary diversification and the organization of ecological networks; and (3) the micro- and macroevolutionary mechanisms and expected patterns under their hypothesis. By doing so, we discuss potentially overlooked aspects and future directions for the study of coevolutionary dynamics and diversification.

Niche breadth and biodiversity change derived from marine Amphipoda species off Iceland

Understanding the ecological requirements and thresholds of individual species is crucial to better predict potential outcomes of climate change on species distribution. In particular, species optima and lower and upper limits along resource gradients require attention. Based on Huisman-Olff-Fresco (HOF) models, we determined species-specific responses along gradients of nine environmental parameters including depth in order to estimate niche attributes of 30 deep-sea benthic amphipods occurring around Iceland. We, furthermore, examined the relationships between niche breadth, occupancy, and geographic range assuming that species with a wider niche are spatially more widely dispersed and vice versa. Overall, our results reveal that species react very differently to environmental gradients, which is independent of the family affiliation of the respective species. We could infer a strong relationship between occupancy and geographic range and also relate this to differences in niche breadth; that is specialist species with a narrow niche had a more limited distribution and may thus be more threatened by changing environmental conditions than generalist species, which are more widespread. Given the preponderance of rare species in the deep sea, this implies that many species could be at risk. However, this must be carefully weighed against geographical data gaps in this area, given that many deep-sea areas are severely undersampled and the true distribution of most species is unknown. After all, our results underline that an accurate taxonomic classification is of crucial importance, without which ecological niche properties cannot be determined and which is hence fundamental for the assessment and understanding of changes in biodiversity in the face of increasing human perturbations.

Species persistence with hybridization in toad-headed lizards driven by divergent selection and low recombination

Speciation plays a central role in evolutionary studies, and particularly how reproductive isolation (RI) evolves. The origins and persistence of RI are distinct processes that require separate evaluations. Treating them separately clarifies the drivers of speciation and then it is possible to link the processes to understand large-scale patterns of diversity. Recent genomic studies have focused predominantly on how species or RI originate. However, we know little about how species persist in face of gene flow. Here, we evaluate a contact zone of two closely related toad-headed lizards (Phrynocephalus) using a chromosome-level genome assembly and population genomics. To some extent, recent asymmetric introgression from Phrynocephalus putjatai to P. vlangalii reduces their genomic differences. However, their highly divergent regions (HDRs) have heterogeneous distributions across the genomes. Functional gene annotation indicates that many genes within HDRs are involved in reproduction and RI. Compared with allopatric populations, contact areas exhibit recent divergent selection on the HDRs and a lower population recombination rate. Taken together, this implies that divergent selection and low genetic recombination help maintain RI. This study provides insights into the genomic mechanisms that drive RI and two species persistence in the face of gene flow during the late stage of speciation.

No link between population isolation and speciation rate in squamate reptiles

Rates of species formation vary widely across the tree of life and contribute to massive disparities in species richness among clades. This variation can emerge from differences in metapopulation-level processes that affect the rates at which lineages diverge, persist, and evolve reproductive barriers and ecological differentiation. For example, populations that evolve reproductive barriers quickly should form new species at faster rates than populations that acquire reproductive barriers more slowly. This expectation implicitly links microevolutionary processes (the evolution of populations) and macroevolutionary patterns (the profound disparity in speciation rate across taxa). Here, leveraging extensive field sampling from the Neotropical Cerrado biome in a biogeographically controlled natural experiment, we test the role of an important microevolutionary process-the propensity for population isolation-as a control on speciation rate in lizards and snakes. By quantifying population genomic structure across a set of codistributed taxa with extensive and phylogenetically independent variation in speciation rate, we show that broad-scale patterns of species formation are decoupled from demographic and genetic processes that promote the formation of population isolates. Population isolation is likely a critical stage of speciation for many taxa, but our results suggest that interspecific variability in the propensity for isolation has little influence on speciation rates. These results suggest that other stages of speciation-including the rate at which reproductive barriers evolve and the extent to which newly formed populations persist-are likely to play a larger role than population isolation in controlling speciation rate variation in squamates.

Faster evolution of a premating reproductive barrier is not associated with faster speciation rates in New World passerine birds

Why are speciation rates so variable across the tree of life? One hypothesis is that this variation is explained by how rapidly reproductive barriers evolve. We tested this hypothesis by conducting a comparative study of the evolution of bird song, a premating barrier to reproduction. Speciation in birds is typically initiated when geographically isolated (allopatric) populations evolve reproductive barriers. We measured the strength of song as a premating barrier between closely related allopatric populations by conducting 2339 field experiments to measure song discrimination for 175 taxon pairs of allopatric or parapatric New World passerine birds, and estimated recent speciation rates from molecular phylogenies. We found evidence that song discrimination is indeed an important reproductive barrier: taxon pairs with high song discrimination in allopatry did not regularly interbreed in parapatry. However, evolutionary rates of song discrimination were not associated with recent speciation rates. Evolutionary rates of song discrimination were also unrelated to latitude or elevation, but species with innate song (suboscines) evolved song discrimination much faster than species with learned song (oscines). We conclude that song is a key premating reproductive barrier in birds, but faster evolution of this reproductive barrier between populations does not consistently result in faster diversification between species.

Causes and Consequences of Apparent Timescaling Across All Estimated Evolutionary Rates

Evolutionary rates play a central role in connecting micro- and macroevolution. All evolutionary rate estimates, including rates of molecular evolution, trait evolution, and lineage diversification, share a similar scaling pattern with time: The highest rates are those measured over the shortest time interval. This creates a disconnect between micro- and macroevolution, although the pattern is the opposite of what some might expect: Patterns of change over short timescales predict that evolution has tremendous potential to create variation and that potential is barely tapped by macroevolution. In this review, we discuss this shared scaling pattern across evolutionary rates. We break down possible explanations for scaling into two categories, estimation error and model misspecification, and discuss how both apply to each type of rate. We also discuss the consequences of this ubiquitous pattern, which can lead to unexpected results when comparing ratesover different timescales. Finally, after addressing purely statistical concerns, we explore a few possibilities for a shared unifying explanation across the three types of rates that results from a failure to fully understand and account for how biological processes scale over time.

Incipient speciation and its impact on taxonomic decision: a case study using a sky island sister-species pair of stag beetles (Lucanidae: Lucanus)

Species delimitation can be difficult when the divergence between focal taxa is in the incipient stage of speciation, because conflicting results are expected among different data sets, and the species limits can differ depending on the species concept applied. We studied speciation history and investigated the impact on taxonomic decision-making when using different types of data in a Taiwanese endemic sister-species pair of stag beetles, Lucanus miwai and Lucanus yulaoensis, from sky island habitats. We showed that the two geographical taxa can be diagnosed by male mandibular shape. We found two mitochondrial co1 lineages with pairwise sequence divergence > 3%; however, L. miwai might not be monophyletic. The result of our multispecies coalescent-based species delimitation using five nuclear loci supported the evolutionary independence of the two sister species, but the calculated values of the genealogical divergence index (gdi) corresponded to the ambiguous zone of species delimitation. We also showed that post-divergence gene flow is unlikely. Our study demonstrates challenges in the delineation of incipient species, but shows the importance of understanding the speciation history and adopting integrative approaches to reconcile seemingly conflicting results before making evolutionarily relevant taxonomic decisions.

Geographical isolation versus dispersal: Relictual alpine grasshoppers support a model of interglacial diversification with limited hybridization

Alpine biotas are paradigmatic of the countervailing roles of geographical isolation and dispersal during diversification. In temperate regions, repeated distributional shifts driven by Pleistocene climatic oscillations produced both recurrent pulses of population fragmentation and opportunities for gene flow during range expansions. Here, we test whether a model of divergence in isolation vs. with gene flow is more likely in the diversification of flightless alpine grasshoppers of the genus Podisma from the Iberian Peninsula. The answer to this question can also provide key insights about the pace of evolution. Specifically, if the data fit a divergence in isolation model, this suggests rapid evolution of reproductive isolation. Genomic data confirm a Pleistocene origin of the species complex, and multiple analytical approaches revealed limited asymmetric historical hybridization between two taxa. Genomic-based demographic reconstructions, spatial patterns of genetic structure and range shifts inferred from palaeodistribution modelling suggest severe range contraction accompanied by declines in effective population sizes during interglacials (i.e., contemporary populations confined to sky islands are relicts) and expansions during the coldest stages of the Pleistocene in each taxon. Although limited hybridization during secondary contact leads to phylogenetic uncertainty if gene flow is not accommodated when estimating evolutionary relationships, all species exhibit strong genetic cohesiveness. Our study lends support to the notion that the accumulation of incipient differences during periods of isolation were sufficient to lead to lineage persistence, but also that the demographic changes, dispersal constraints and spatial distribution of the sky islands themselves mediated species diversification in temperate alpine biotas.

Plotting for change: an analytical framework to aid decisions on which lineages are candidate species in phylogenomic species discovery

A recent study argued that coalescent-based models of species delimitation mostly delineate population structure, not species, and called for the validation of candidate species using biological information additional to the genetic information, such as phenotypic or ecological data. Here, we introduce a framework to interrogate genomic datasets and coalescent-based species trees for the presence of candidate species in situations where additional biological data are unavailable, unobtainable or uninformative. For de novo genomic studies of species boundaries, we propose six steps: (1) visualize genetic affinities among individuals to identify both discrete and admixed genetic groups from first principles and to hold aside individuals involved in contemporary admixture for independent consideration; (2) apply phylogenetic techniques to identify lineages; (3) assess diagnosability of those lineages as potential candidate species; (4) interpret the diagnosable lineages in a geographical context (sympatry, parapatry, allopatry); (5) assess significance of difference or trends in the context of sampling intensity; and (6) adopt a holistic approach to available evidence to inform decisions on species status in the difficult cases of allopatry. We adopt this approach to distinguish candidate species from within-species lineages for a widespread species complex of Australian freshwater fishes (Retropinna spp.). Our framework addresses two cornerstone issues in systematics that are often not discussed explicitly in genomic species discovery: diagnosability and how to determine it, and what criteria should be used to decide whether diagnosable lineages are conspecific or represent different species.

Phytoplankton biodiversity and the inverted paradox

Earth's aquatic food webs are overwhelmingly supported by planktonic microalgae that live in the sunlit water column where only a minimum number of physical niches are readily identifiable. Despite this paucity of environmental differentiation, these "phytoplankton" populations exhibit a rich biodiversity, an observation not easily reconciled with broadly accepted rules of resource-based competitive exclusion. This conundrum is referred to as the "Paradox of the Plankton". Consideration of physical distancing between nutrient depletion zones around individual phytoplankton, however, suggests a competition-neutral resource landscape. Application of neutral theory to the sheer number of phytoplankton in physically-mixed water masses yields a prediction of astronomical biodiversity, suggesting the inverted paradox: Why are there so few phytoplankton species? Here, we introduce a trophic constraint on phytoplankton that, when combined with stochastic principals of ecological drift, predicts only modest levels of diversity in an otherwise competition-neutral landscape. Our "trophic exclusion" principle predicts diversity to be independent of population size and yields a species richness across cell-size classes that is consistent with broad oceanographic survey observations.

Linking population-level and microevolutionary processes to understand speciation dynamics at the macroevolutionary scale

Although speciation dynamics have been described for several taxonomic groups in distinct geographic regions, most macroevolutionary studies still lack a detailed mechanistic view on how or why speciation rates change. To help partially fill this gap, we suggest that the interaction between the time taken by a species to geographically expand and the time populations take to evolve reproductive isolation should be considered when we are trying to understand macroevolutionary patterns. We introduce a simple conceptual index to guide our discussion on how demographic and microevolutionary processes might produce speciation dynamics at macroevolutionary scales. Our framework is developed under different scenarios: when speciation is mediated by geographical or resource-partitioning opportunities, and when diversity is limited or not. We also discuss how organismal intrinsic properties and different overall geographical settings can influence the tempo and mode of speciation. We argue that specific conditions observed at the microscale might produce a pulse in speciation rates even without a pulse in either climate or physical barriers. We also propose a hypothesis to reconcile the apparent inconsistency between speciation measured at the microscale and macroscale, and emphasize that diversification rates are better seen as an emergent property. We hope to bring the reader's attention to interesting mechanisms to be further studied, to motivate the development of new theoretical models that connect microevolution and macroevolution, and to inspire new empirical and methodological approaches to more adequately investigate speciation dynamics either using neontological or paleontological data.

Incorporating the speciation process into species delimitation

The "multispecies" coalescent (MSC) model that underlies many genomic species-delimitation approaches is problematic because it does not distinguish between genetic structure associated with species versus that of populations within species. Consequently, as both the genomic and spatial resolution of data increases, a proliferation of artifactual species results as within-species population lineages, detected due to restrictions in gene flow, are identified as distinct species. The toll of this extends beyond systematic studies, getting magnified across the many disciplines that rely upon an accurate framework of identified species. Here we present the first of a new class of approaches that addresses this issue by incorporating an extended speciation process for species delimitation. We model the formation of population lineages and their subsequent development into independent species as separate processes and provide for a way to incorporate current understanding of the species boundaries in the system through specification of species identities of a subset of population lineages. As a result, species boundaries and within-species lineages boundaries can be discriminated across the entire system, and species identities can be assigned to the remaining lineages of unknown affinities with quantified probabilities. In addition to the identification of species units in nature, the primary goal of species delimitation, the incorporation of a speciation model also allows us insights into the links between population and species-level processes. By explicitly accounting for restrictions in gene flow not only between, but also within, species, we also address the limits of genetic data for delimiting species. Specifically, while genetic data alone is not sufficient for accurate delimitation, when considered in conjunction with other information we are able to not only learn about species boundaries, but also about the tempo of the speciation process itself.

Reconstructing hotspots of genetic diversity from glacial refugia and subsequent dispersal in Italian common toads (Bufo bufo)

Genetic diversity feeds the evolutionary process and allows populations to adapt to environmental changes. However, we still lack a thorough understanding of why hotspots of genetic diversity are so 'hot'. Here, we analysed the relative contribution of bioclimatic stability and genetic admixture between divergent lineages in shaping spatial patterns of genetic diversity in the common toad Bufo bufo along the Italian peninsula. We combined population genetic, phylogeographic and species distribution modelling (SDM) approaches to map ancestral areas, glacial refugia, and secondary contact zones. We consistently identified three phylogeographic lineages, distributed in northern, central and southern Italy. These lineages expanded from their ancestral areas and established secondary contact zones, before the last interglacial. SDM identified widespread glacial refugia in peninsular Italy, sometimes located under the present-day sea-level. Generalized linear models indicated genetic admixture as the only significant predictor of the levels of population genetic diversity. Our results show that glacial refugia contributed to preserving both levels and patterns of genetic diversity across glacial-interglacial cycles, but not to their formation, and highlight a general principle emerging in Mediterranean species: higher levels of genetic diversity mark populations with substantial contributions from multiple genetic lineages, irrespective of the location of glacial refugia.

On the Origin of Coexisting Species

Speciation is frequently initiated but rarely completed, a phenomenon hypothesized to arise due to the failure of nascent lineages to persist. Although a failure to persist often has ecological causes, key gaps exist between ecological and evolutionary theories that, if filled, would clarify when and why speciation succeeds or fails. Here, we apply ecological coexistence theory to show how the alignment between different forms of niche opportunity and niche use shape the initiation, progression, and completion of speciation. Niche evolution may drive coexistence or competitive exclusion, and an ability to coexist ecologically may help or hinder speciation. Our perspective allows progress towards unifying the origin and maintenance of species diversity across the tree of life.

Is population subdivision different from speciation? From phylogeography to species delimitation

Species-level diversity and the underlying mechanisms that lead to the formation of new species, that is, speciation, have often been confounded with intraspecific diversity and population subdivision. The delineation between intraspecific and interspecific divergence processes has received much less attention than species delimitation. The ramifications of confounding speciation and population subdivision are that the term speciation has been used to describe many different biological divergence processes, rendering the results, or inferences, between studies incomparable. Phylogeographic studies have advanced our understanding of how spatial variation in the pattern of biodiversity can begin, become structured, and persist through time. Studies of species delimitation have further provided statistical and model-based approaches to determine the phylogeographic entities that merit species status. However, without a proper understanding and delineation between the processes that generate and maintain intraspecific and interspecific diversity in a study system, the delimitation of species may still not be biologically and evolutionarily relevant. I argue that variation in the continuity of the divergence process among biological systems could be a key factor leading to the enduring contention in delineating divergence patterns, or species delimitation, meriting future comparative studies to help us better understand the nature of biological species.

Terrestrial habitats decouple the relationship between species and subspecies diversification in mammals

Darwin proposed that lineages with higher diversification rates should evidence this capacity at both the species and subspecies level. This should be the case if subspecific boundaries are evolutionary faultlines along which speciation is generally more likely to occur. This pattern has been described for birds, but remains poorly understood in mammals. To investigate the relationship between species richness (SR) and subspecies richness (SSR), we calculated the strength of the correlation between the two across all mammals. Mammalian taxonomic richness correlates positively, but only very weakly, between the species and subspecies level, deviating from the pattern found in birds. However, when mammals are separated by environmental substrate, the relationship between generic SR and average SSR in non-terrestrial taxa is stronger than that reported for birds (Kendall's tau = 0.31, p < 0.001). By contrast, the correlation in terrestrial taxa alone weakens compared to that for all mammals (Kendall's tau = 0.11, p < 0.001). A significant interaction between environmental substrate and SR in phylogenetic regressions confirms a role for terrestrial habitats in disrupting otherwise linked dynamics of diversification across the taxonomic hierarchy. Further, models including species range size as a predictor show that range size affects SSR more in terrestrial taxa. Taken together, these results suggest that the dynamics of diversification of terrestrial mammals are more affected by physical barriers or ecological heterogeneity within ranges than those of non-terrestrial mammals, at two evolutionary levels. We discuss the implication of these results for the equivalence of avian and mammalian subspecies, their potential role in speciation and the broader question of the relationship between microevolution and macroevolution.

Beyond Reproductive Isolation: Demographic Controls on the Speciation Process

Studies of speciation typically investigate the evolution of reproductive isolation between populations, but several other processes can serve as key steps limiting the formation of species. In particular, the probability of successful speciation can be influenced by factors that affect the frequency with which population isolates form as well as their persistence through time. We suggest that population isolation and persistence have an inherently spatial dimension that can be profitably studied using a conceptual framework drawn from metapopulation ecology. We discuss models of speciation that incorporate demographic processes and highlight the need for a broader application of phylogenetic comparative approaches to evaluate the general importance of population isolation, persistence, and reproductive isolation in speciation. We review diverse and nontraditional data sources that can be leveraged to study isolation and persistence in a comparative framework. This incorporation of spatial demographic information facilitates the integration of perspectives on speciation across disciplines and timescales.

Dispersal and local persistence shape the genetic structure of a widespread Neotropical plant species with a patchy distribution

BACKGROUND AND AIMS

Isolated populations constitute an ideal laboratory to study the consequences of intraspecific divergence, because intrinsic incompatibilities are more likely to accumulate under reduced gene flow. Here, we use a widespread bromeliad with a patchy distribution, Pitcairnia lanuginosa, as a model to infer processes driving Neotropical diversification and, thus, to improve our understanding of the origin and evolutionary dynamics of biodiversity in this highly speciose region.

METHODS

We assessed the timing of lineage divergence, genetic structural patterns and historical demography of P. lanuginosa, based on microsatellites, and plastid and nuclear sequence data sets using coalescent analyses and an Approximate Bayesian Computation framework. Additionally, we used species distribution models (SDMs) to independently estimate potential changes in habitat suitability.

KEY RESULTS

Despite morphological uniformity, plastid and nuclear DNA data revealed two distinct P. lanuginosa lineages that probably diverged through dispersal from the Cerrado to the Central Andean Yungas, following the final uplift of the Andes, and passed through long-term isolation with no evidence of migration. Microsatellite data indicate low genetic diversity and high levels of inbreeding within populations, and restricted gene flow among populations, which are likely to be a consequence of bottlenecks (or founder events), and high selfing rates promoting population persistence in isolation. SDMs showed a slight expansion of the suitable range for P. lanuginosa lineages during the Last Glacial Maximum, although molecular data revealed a signature of older divergence. Pleistocene climatic oscillations thus seem to have played only a minor role in the diversification of P. lanuginosa, which probably persisted through adverse conditions in riparian forests.

CONCLUSIONS

Our results imply drift as a major force shaping the evolution of P. lanuginosa, and suggest that dispersal events have a prominent role in connecting Neotropical open and forest biomes.

Exploring the drivers of population structure across desert snakes can help to link micro and macroevolution

To understand the underlying mechanisms generating population genetic divergence and structure is a critical step towards understanding how biodiversity evolves at both micro- and macroevolutionary scales. At the population-level, geographic isolation as well as adaptation to local environmental conditions can generate different patterns of spatial genetic variation among populations. Specific organismal traits as well as the characteristics of the environment might influence the process under which populations become spatially structured. In a From the Cover article in this issue of Molecular Ecology, Myers et al. (2019) present an integrative approach to investigate if the Cochise filter barrier (CFB), lying between the Sonoran and Chihuahuan Deserts, and the surrounding river networks were relevant in driving the population structure of 13 snake species. While local environmental conditions seem to predominantly contribute to lineage divergence, traditionally studied vicariant barriers seem to have played a minor role in shaping population structure across the studied species. This study brings insights into how population-level processes could contribute to the formation of incipient species, which ultimately might affect the speciation rates measured at macroevolutionary scales. Hence, Myers et al. (2019) not only represents an integrative study aiming to understand the drivers of population genetic divergence, but also a potentially important contribution to our ongoing challenge in linking micro- and macroevolution.

Understanding the effect of competition during evolutionary radiations: an integrated model of phenotypic and species diversification

Competition can drive macroevolutionary change, for example during adaptive radiations. However, we still lack a clear understanding of how it shapes diversification processes and patterns. To better understand the macroevolutionary consequences of competition, as well as the signal left on phylogenetic data, we developed a model linking trait evolution and species diversification in an ecological context. We find four main results: first, competition spurs trait diversity but not necessarily species richness; second, competition produces slowdowns in species diversification even in the absence of explicit ecological limits, but not in phenotypic diversification even in the presence of such limits; third, early burst patterns do not provide a reliable way of testing for adaptive radiations; and fourth, looking for phylogenetic signal in trait data and support for phenotypic models incorporating competition is a better alternative. Our results clarify the macroevolutionary consequences of competition and could help design more powerful tests of adaptive radiations in nature.

Why are red flowers so rare? Testing the macroevolutionary causes of tippiness

Traits that have arisen multiple times yet still remain rare present a curious paradox. A number of these rare traits show a distinct tippy pattern, where they appear widely dispersed across a phylogeny, are associated with short branches and differ between recently diverged sister species. This phylogenetic pattern has classically been attributed to the trait being an evolutionary dead end, where the trait arises due to some short-term evolutionary advantage, but it ultimately leads species to extinction. While the higher extinction rate associated with a dead end trait could produce such a tippy pattern, a similar pattern could appear if lineages with the trait speciated slower than other lineages, or if the trait was lost more often that it was gained. In this study, we quantify the degree of tippiness of red flowers in the tomato family, Solanaceae, and investigate the macroevolutionary processes that could explain the sparse phylogenetic distribution of this trait. Using a suite of metrics, we confirm that red-flowered lineages are significantly overdispersed across the tree and form smaller clades than expected under a null model. Next, we fit 22 alternative models using HiSSE (Hidden State Speciation and Extinction), which accommodates asymmetries in speciation, extinction and transition rates that depend on observed and unobserved (hidden) character states. Results of the model fitting indicated significant variation in diversification rates across the family, which is best explained by the inclusion of hidden states. Our best fitting model differs between the maximum clade credibility tree and when incorporating phylogenetic uncertainty, suggesting that the extreme tippiness and rarity of red Solanaceae flowers makes it difficult to distinguish among different underlying processes. However, both of the best models strongly support a bias towards the loss of red flowers. The best fitting HiSSE model when incorporating phylogenetic uncertainty lends some support to the hypothesis that lineages with red flowers exhibit reduced diversification rates due to elevated extinction rates. Future studies employing simulations or targeting population-level processes may allow us to determine whether red flowers in Solanaceae or other angiosperms clades are rare and tippy due to a combination of processes, or asymmetrical transitions alone.

Microevolutionary processes impact macroevolutionary patterns

BACKGROUND

Macroevolutionary modeling of species diversification plays important roles in inferring large-scale biodiversity patterns. It allows estimation of speciation and extinction rates and statistically testing their relationships with different ecological factors. However, macroevolutionary patterns are ultimately generated by microevolutionary processes acting at population levels, especially when speciation and extinction are considered protracted instead of point events. Neglecting the connection between micro- and macroevolution may hinder our ability to fully understand the underlying mechanisms that drive the observed patterns.

RESULTS

In this simulation study, we used the protracted speciation framework to demonstrate that distinct microevolutionary scenarios can generate very similar biodiversity patterns (e.g., latitudinal diversity gradient). We also showed that current macroevolutionary models may not be able to distinguish these different scenarios.

CONCLUSIONS

Given the compounded nature of speciation and extinction rates, one needs to be cautious when inferring causal relationships between ecological factors and macroevolutioanry rates. Future studies that incorporate microevolutionary processes into current modeling approaches are in need.

Testing the impact of effective population size on speciation rates - a negative correlation or lack thereof in lichenized fungi

The effect of microevolutionary processes on macroevolutionary patterns, and vice versa, is an important but under-investigated question. Here we present an integrative population genetic and phylogenetic study using molecular sequence data from three lichen-forming fungal lineages to empirically test the potential correlation between effective population size – approximated by the parameter θ – and estimated speciation rates using a phylogenetic tree (λ). A negative association between θ and λ was supported in one lineage of lichen-forming fungi, Melanelixia (Parmeliaceae), while no significant relationships was found for two other genera within the same family, Melanohalea and Xanthoparmelia. We discuss the significance of our results and the importance of considering microevolutionary processes when studying macroevolutionary patterns.

Factors promoting polyploid persistence and diversification and limiting diploid speciation during the K-Pg interlude

The large wave of polyploidization following the Cretaceous-Paleogene (K-Pg) mass extinction has been explained by enhanced polyploid persistence arising from adaptive properties of the polyploids themselves, as well as an increase in unreduced gamete production and diploid hybridization. We propose that the demise of diploids afforded opportunities for polyploid establishment and expansion into novel habitats. Augmented polyploid gene pools from diploid and polyploid relatives, in association with their multiple and independent origins (of both autopolyploids and allopolyploids), facilitated their subsequent diversification. Their ability to recruit genetic variation from their diploid relatives or from products of recurrent origins sharing their genome(s) ostensibly contributed to polyploid persistence. Concomitantly, we propose that the number of congeneric diploid species dramatically contracted disproportionally to polyploids during the K-Pg interval (i.e. a diploid trough), resulting in a reduction in the rate of diploid speciation. Accordingly, the preponderance of neopolyploids was likely autopolyploids.

Commensal associations and benthic habitats shape macroevolution of the bivalve clade Galeommatoidea

The great diversity of marine life has been shaped by the interplay between abiotic and biotic factors. Among different biotic interactions, symbiosis is an important yet less studied phenomenon. Here, we tested how symbiotic associations affected marine diversification, using the bivalve superfamily Galeommatoidea as a study system. This superfamily contains large numbers of obligate commensal as well as free-living species and is therefore amenable to comparative approaches. We constructed a global molecular phylogeny of Galeommatoidea and compared macroevolutionary patterns between free-living and commensal lineages. Our analyses inferred that commensalism/sediment-dwelling is likely to be the ancestral condition of Galeommatoidea and that secondary invasions of hard-bottom habitats linked to the loss of commensalism. One major clade containing most of the free-living species exhibits a 2-4 times higher diversification rate than that of the commensals, likely driven by frequent niche partitioning in highly heterogeneous hard-bottom habitats. However, commensal clades show much higher within-clade morphological disparity, likely promoted by their intimate associations with diverse hosts. Our study highlights the importance of interactions between different ecological factors in shaping marine macroevolution and that biotic factors cannot be ignored if we wish to fully understand processes that generate marine biodiversity.