Rapid vertebrate speciation via isolation, bottlenecks, and drift.

Speciation is often driven by selective processes like those associated with viability, mate choice, or local adaptation, and "speciation genes" have been identified in many eukaryotic lineages. In contrast, neutral processes are rarely considered as the primary drivers of speciation, especially over short evolutionary timeframes. Here, we describe a rapid vertebrate speciation event driven primarily by genetic drift. The White Sands pupfish (Cyprinodon tularosa) is endemic to New Mexico's Tularosa Basin where the species is currently managed as two Evolutionarily significant units (ESUs) and is of international conservation concern (Endangered). Whole-genome resequencing data from each ESU showed remarkably high and uniform levels of differentiation across the entire genome (global FST ≈ 0.40). Despite inhabiting ecologically dissimilar springs and streams, our whole-genome analysis revealed no discrete islands of divergence indicative of strong selection, even when we focused on an array of candidate genes. Demographic modeling of the joint allele frequency spectrum indicates the two ESUs split only ~4 to 5 kya and that both ESUs have undergone major bottlenecks within the last 2.5 millennia. Our results indicate the genome-wide disparities between the two ESUs are not driven by divergent selection but by neutral drift due to small population sizes, geographic isolation, and repeated bottlenecks. While rapid speciation is often driven by natural or sexual selection, here we show that isolation and drift have led to speciation within a few thousand generations. We discuss these evolutionary insights in light of the conservation management challenges they pose.

Central African dwarf crocodiles found in syntopy are comparably divergent to South American dwarf caimans.

Recent molecular taxonomic advancements have expanded our understanding of crocodylian diversity, revealing the existence of previously overlooked species, including the Congo dwarf crocodile (Osteolaemus osborni) in the central Congo Basin rainforests. This study explores the genomic divergence between O. osborni and its better-known relative, the true dwarf crocodile (Osteolaemus tetraspis), shedding light on their evolutionary history. Field research conducted in the northwestern Republic of the Congo uncovered a locality where both species coexist in sympatry/syntopy. Genomic analysis of sympatric individuals reveals a level of divergence comparable to that between ecologically similar South American dwarf caimans (Paleosuchus palpebrosus and Paleosuchus trigonatus), suggesting parallel speciation in the Afrotropics and Neotropics during the Middle to Late Miocene, 10-12 Ma. Comparison of the sympatric and allopatric dwarf crocodiles indicates no gene flow between the analysed sympatric individuals of O. osborni and O. tetraspis. However, a larger sample will be required to answer the question of whether or to what extent these species hybridize. This study emphasizes the need for further research on the biology and conservation status of the Congo dwarf crocodile, highlighting its significance in the unique biodiversity of the Congolian rainforests and thus its potential as a flagship species.

Extreme elevational migration spurred cryptic speciation in giant hummingbirds.

The ecoevolutionary drivers of species niche expansion or contraction are critical for biodiversity but challenging to infer. Niche expansion may be promoted by local adaptation or constrained by physiological performance trade-offs. For birds, evolutionary shifts in migratory behavior permit the broadening of the climatic niche by expansion into varied, seasonal environments. Broader niches can be short-lived if diversifying selection and geography promote speciation and niche subdivision across climatic gradients. To illuminate niche breadth dynamics, we can ask how "outlier" species defy constraints. Of the 363 hummingbird species, the giant hummingbird (Patagona gigas) has the broadest climatic niche by a large margin. To test the roles of migratory behavior, performance trade-offs, and genetic structure in maintaining its exceptional niche breadth, we studied its movements, respiratory traits, and population genomics. Satellite and light-level geolocator tracks revealed an >8,300-km loop migration over the Central Andean Plateau. This migration included a 3-wk, ~4,100-m ascent punctuated by upward bursts and pauses, resembling the acclimatization routines of human mountain climbers, and accompanied by surging blood-hemoglobin concentrations. Extreme migration was accompanied by deep genomic divergence from high-elevation resident populations, with decisive postzygotic barriers to gene flow. The two forms occur side-by-side but differ almost imperceptibly in size, plumage, and respiratory traits. The high-elevation resident taxon is the world's largest hummingbird, a previously undiscovered species that we describe and name here. The giant hummingbirds demonstrate evolutionary limits on niche breadth: when the ancestral niche expanded due to evolution (or loss) of an extreme migratory behavior, speciation followed.

Paedomorphosis and retention of juvenile diet lead speciation in a group of Neotropical snakes (Colubroides-Philodryadini).

Dipsadidae is one of the largest clades of extant reptiles, showing an impressive morphological and ecological diversity. Despite this fact, the developmental processes behind its diversity are still largely unknown. In this study, we used 3D reconstructions based on micro-CT data and geometric morphometrics to evaluate the skull morphology of Philodryas agassizii, a small, surface-dwelling dipsadid that consume spiders. Adult individuals of P. agassizii exhibit a cranial morphology frequently observed in juveniles of other surface-dwelling colubroideans, represented in our analysis by its close relative Philodryas patagoniensis. Large orbits, gibbous neurocranial roof and a relatively short jaw complex are features present in juveniles of the latter species. Furthermore, we performed an extensive survey about diet of P. patagoniensis in which we detected an ontogenetic dietary shift, indicating that arthropods are more frequently consumed by juveniles of this dietary generalist. Thus, we infer that P. agassizzii retained not only the ancestral juvenile skull morphology but also dietary preferences. This study reveals that morphological changes driven by heterochronic changes, specifically paedomorphosis, influenced the retention of ancestral life history traits in P. agassizii, and therefore promoted cladogenesis. In this way, we obtained first evidence that heterochronic processes lead speciation in the snake megadiverse clade Dipsadidae.

Speciation patterns of Aedes mosquitoes in the Scutellaris Group: a mitochondrial perspective.

The Scutellaris Group of Aedes comprises 47 mosquito species, including Aedes albopictus. While Ae. albopictus is widely distributed, the other species are mostly found in the Asia-Pacific region. Evolutionary history researches of Aedes species within the Scutellaris Group have mainly focused on Ae. albopictus, a species that raises significant public health concerns, neglecting the other species. In this study, we aimed to assess genetic diversity and estimate speciation times of several species within the Scutellaris Group. Mosquitoes were therefore collected from various Asia-Pacific countries. Their mitochondrial cytochrome c oxidase subunit 1 (cox1) and subunit 3 (cox3) sequences were analyzed alongside those of other Scutellaris Group species available in the GenBank database. To estimate the divergence time, we analyzed 1849 cox1 gene sequences from 21 species, using three species (Aedes aegypti, Aedes notoscriptus and Aedes vigilax) as outgroups. We found that most of the speciation dates occurred during the Paleogene and the Neogene periods. A separation between the Scutellaris Subgroup and the Albopictus Subgroup occurred approximately 64-61 million years ago (MYA). We also identified a split between species found in Asia/Micronesia and those collected in Melanesia/Polynesia approximately 36-35 MYA. Our findings suggest that the speciation of Aedes species within the Scutellaris Group may be driven by diversity in mammalian hosts, climate and environmental changes, and geological dynamics rather than human migration.

A new species of Populus and the extensive hybrid speciation arising from it on the Qinghai-Tibet Plateau.

While the diversity of species formation is broadly acknowledged, significant debate exists regarding the universal nature of hybrid species formation. Through an 18-year comprehensive study of all Populus species on the Qinghai-Tibet Plateau, 23 previously recorded species and 8 new species were identified. Based on morphological characteristics, these can be classified into three groups: species in section Leucoides, species with large leaves, and species with small leaves in section Tacamahaca. By conducting whole-genome re-sequencing of 150 genotypes from these 31 species, 2.28 million single nucleotide polymorphisms (SNPs) were identified. Phylogenetic analysis utilizing these SNPs not only revealed a highly intricate evolutionary network within the large-leaf species of section Tacamahaca but also confirmed that a new species, P. curviserrata, naturally hybridized with P. cathayana, P. szechuanica, and P. ciliata, resulting in 11 hybrid species. These findings indicate the widespread occurrence of hybrid species formation within this genus, with hybridization serving as a key evolutionary mechanism for Populus on the plateau. A novel hypothesis, "Hybrid Species Exterminating Their Ancestral Species (HSEAS)," is introduced to explain the mechanisms of hybrid species formation at three different scales: the entire plateau, the southeastern mountain region, and individual river valleys.

Hybrid speciation driven by multilocus introgression of ecological traits.

Hybridization allows adaptations to be shared among lineages and may trigger the evolution of new species1,2. However, convincing examples of homoploid hybrid speciation remain rare because it is challenging to demonstrate that hybridization was crucial in generating reproductive isolation3. Here we combine population genomic analysis with quantitative trait locus mapping of species-specific traits to examine a case of hybrid speciation in Heliconius butterflies. We show that Heliconius elevatus is a hybrid species that is sympatric with both parents and has persisted as an independently evolving lineage for at least 180,000 years. This is despite pervasive and ongoing gene flow with one parent, Heliconius pardalinus, which homogenizes 99% of their genomes. The remaining 1% introgressed from the other parent, Heliconius melpomene, and is scattered widely across the H. elevatus genome in islands of divergence from H. pardalinus. These islands contain multiple traits that are under disruptive selection, including colour pattern, wing shape, host plant preference, sex pheromones and mate choice. Collectively, these traits place H. elevatus on its own adaptive peak and permit coexistence with both parents. Our results show that speciation was driven by introgression of ecological traits, and that speciation with gene flow is possible with a multilocus genetic architecture.

A new nuclear phylogeny of the tea family (Theaceae) unravels rapid radiations in genus Camellia.

Molecular analyses of rapidly radiating groups often reveal incongruence between gene trees. This mainly results from incomplete lineage sorting, introgression, and gene tree estimation error, which complicate the estimation of phylogenetic relationships. In this study, we reconstruct the phylogeny of Theaceae using 348 nuclear loci from 68 individuals and two outgroup taxa. Sequence data were obtained by target enrichment using the recently released Angiosperm 353 universal probe set applied to herbarium specimens. The robustness of the topologies to variation in data quality was established under a range of different filtering schemes, using both coalescent and concatenation approaches. Our results confirmed most of the previously hypothesized relationships among tribes and genera, while clarifying additional interspecific relationships within the rapidly radiating genus Camellia. We recovered a remarkably high degree of gene tree heterogeneity indicative of rapid radiation in the group and observed cytonuclear conflicts, especially within Camellia. This was especially pronounced around short branches, which we primarily associate with gene tree estimation error. Our analysis also indicates that incomplete lineage sorting (ILS) contributed to gene-tree conflicts and accounted for approximately 14 % of the explained variation, whereas inferred introgression levels were low. Our study advances the understanding of the evolution of this important plant family and provides guidance on the application of target capture methods and the evaluation of key processes that influence phylogenetic discordances.

Marine introgressions and Andean uplift have driven diversification in neotropical Monkey tree frogs (Anura, Phyllomedusinae).

The species richness in the Neotropics has been linked to environmental heterogeneity and a complex geological history. We evaluated which biogeographic processes were associated with the diversification of Monkey tree frogs, an endemic clade from the Neotropics. We tested two competing hypotheses: the diversification of Phyllomedusinae occurred either in a "south-north" or a "north-south" direction in the Neotropics. We also hypothesized that marine introgressions and Andean uplift had a crucial role in promoting their diversification. We used 13 molecular markers in a Bayesian analysis to infer phylogenetic relationships among 57 species of Phyllomedusinae and to estimate their divergence times. We estimated ancestral ranges based on 12 biogeographic units considering the landscape modifications of the Neotropical region. We found that the Phyllomedusinae hypothetical ancestor range was probably widespread throughout South America, from Western Amazon to Southern Atlantic Forest, at 29.5 Mya. The Phyllomedusines' ancestor must have initially diverged through vicariance, generally followed by jump-dispersals and sympatric speciation. Dispersal among areas occurred mostly from Western Amazonia towards Northern Andes and the South American diagonal of dry landscapes, a divergent pattern from both "south-north" and "north-south" diversification hypotheses. Our results revealed a complex diversification process of Monkey tree frogs, occurring simultaneously with the orogeny of Northern Andes and the South American marine introgressions in the last 30 million years.

Mapping seasonal migration in a songbird hybrid zone -- heritability, genetic correlations, and genomic patterns linked to speciation.

Seasonal migration is a widespread behavior relevant for adaptation and speciation, yet knowledge of its genetic basis is limited. We leveraged advances in tracking and sequencing technologies to bridge this gap in a well-characterized hybrid zone between songbirds that differ in migratory behavior. Migration requires the coordinated action of many traits, including orientation, timing, and wing morphology. We used genetic mapping to show these traits are highly heritable and genetically correlated, explaining how migration has evolved so rapidly in the past and suggesting future responses to climate change may be possible. Many of these traits mapped to the same genomic regions and small structural variants indicating the same, or tightly linked, genes underlie them. Analyses integrating transcriptomic data indicate cholinergic receptors could control multiple traits. Furthermore, analyses integrating genomic differentiation further suggested genes underlying migratory traits help maintain reproductive isolation in this hybrid zone.

Species Delimitation and Exploration of Species Partitions with ASAP and LIMES.

DNA barcoding plays an important role in exploring undescribed biodiversity and is increasingly used to delimit lineages at the species level (see Chap. 4 by Miralles et al.). Although several approaches and programs have been developed to perform species delimitation from datasets of single-locus DNA sequences, such as DNA barcodes, most of these were not initially provided as user-friendly GUI-driven executables. In spite of their differences, most of these tools share the same goal, i.e., inferring de novo a partition of subsets, potentially each representing a distinct species. More recently, a proposed common exchange format for the resulting species partitions (SPART) has been implemented by several of these tools, paving the way toward developing an interoperable digital environment entirely dedicated to integrative and comparative species delimitation. In this chapter, we provide detailed protocols for the use of two bioinformatic tools, one for single locus molecular species delimitation (ASAP) and one for statistical comparison of species partitions resulting from any kind of species delimitation analyses (LIMES).

Rapid clade divergence and phyletic gradualism in an interacting particle model of sympatric speciation.

The coexistence of cladogenesis, i.e., the branching of lineages along an evolutionary tree as observed in the fossil record, and anagenesis, which is the progressive evolution within populations, lacks a clear explanation. In this study, we examine a simple model that simulates the evolutionary changes occurring within populations inhabiting the same environment in sympatry, and driven by ecological competition. Our model characterizes populations through a set of evolving morphological traits represented by mathematical points within a two-dimensional morphospace. Such points may reproduce or die due to overcrowding, implying competition in morphospace as suggested by the ecological phenomenon of character displacement. By focusing on the morphospace rather than physical space, the model effectively captures the simultaneous evolution of coexisting populations. Central to the model is the delicate balance between the range of competition and the range of reproduction within the morphospace. Interesting patterns emerge when the ratio between the competition to reproducetion ranges, referred to as CR ratio, changes from values slightly smaller to significantly larger than unity. When competition acts over short distances relative to the reproduction range (low CR), the phylogenetic tree takes on a nearly uniform appearance, gradually transforming into a more bush-like structure for slightly higher CR values. With further increases in CR, evolutionary lineages become more discernible, and the morphogenetic pattern shifts from a bush-like shape to a more tree-like arrangement and few branches for very large CRs. At specific time sections, the synthetic phylogenetic tree appears as an assembly of clusters of individuals within the morphospace. These clusters, interpretable as simulated models of species, exhibit distinct separation within the morphospace and are subject to dynamic inter-cluster repulsion. Notably, clusters tend to be resistant to change. They maintain relatively constant abundances while gradually shifting their positions within the morphospace-a phase that aligns with the concept of phyletic gradualism. However, this predictable pattern is occasionally upset by the abrupt divisions into multiple groups, interpreted as cladogenesis events. The intricacies of the splitting process are explored, revealing that in scenarios with large CR values, the splitting can emerge much more rapidly than phyletic changes. This accelerated process of splitting is initiated by one or few individuals at the fringes of a cluster, where competition is minimal. The newly generated cluster then undergoes deformation, swiftly followed by divergence and splitting (seen as branching in the synthetic phylogenetic tree), as if an inherent "repulsion" triggered the division between species. The simple rules implied in the interacting-particle model may provide insight into the coexistence of gradualism and cladogenesis along lineages, illustrating the capacity for rapid shifts during cladogenesis and the more gradual process of anagenesis.

The succession of ecological divergence and reproductive isolation in adaptive radiations.

Adaptive radiation is a major source of biodiversity but the way in which known components of ecological opportunity, ecological differentiation, and reproductive isolation underpin such biodiversity patterns remains elusive. Much is known about the evolution of ecological differentiation and reproductive isolation during single speciation events, but exactly how those processes scale up to complete adaptive radiations is less understood. Do we expect complete reproductive barriers between newly formed species before the ecological differentiation continues, or does proper species formation occur much later, long after the ecological diversification? Our goal is to improve our mechanistic understanding of adaptive radiations by analyzing an individual-based model that includes a suite of mechanisms that are known to contribute to biodiversity. The model includes variable biogeographic settings, ecological opportunities, and types of mate choice, which makes several different scenarios of an adaptive radiation possible. We find that evolving clades tend to exploit ecological opportunities early whereas reproductive barriers evolve later, demonstrating a decoupling of ecological differentiation and species formation. In many cases, we also find a long-term trend where assortative mating associated with ecological traits is replaced by sexual selection of neutral display traits as the primary mechanism for reproductive isolation. Our results propose that reticulate phylogenies are likely common and stem from initially low reproductive barriers, rather than the previously suggested idea of repeated hybridization events between well-separated species.

Chromosomal Inversions and the Demography of Speciation in Drosophila montana and Drosophila flavomontana.

Chromosomal inversions may play a central role in speciation given their ability to locally reduce recombination and therefore genetic exchange between diverging populations. We analyzed long- and short-read whole-genome data from sympatric and allopatric populations of 2 Drosophila virilis group species, Drosophila montana and Drosophila flavomontana, to understand if inversions have contributed to their divergence. We identified 3 large alternatively fixed inversions on the X chromosome and one on each of the autosomes 4 and 5. A comparison of demographic models estimated for inverted and noninverted (colinear) chromosomal regions suggests that these inversions arose before the time of the species split. We detected a low rate of interspecific gene flow (introgression) from D. montana to D. flavomontana, which was further reduced inside inversions and was lower in allopatric than in sympatric populations. Together, these results suggest that the inversions were already present in the common ancestral population and that gene exchange between the sister taxa was reduced within inversions both before and after the onset of species divergence. Such ancestrally polymorphic inversions may foster speciation by allowing the accumulation of genetic divergence in loci involved in adaptation and reproductive isolation inside inversions early in the speciation process, while gene exchange at colinear regions continues until the evolving reproductive barriers complete speciation. The overlapping X inversions are particularly good candidates for driving the speciation process of D. montana and D. flavomontana, since they harbor strong genetic incompatibilities that were detected in a recent study of experimental introgression.

Integrative taxonomic analyses reveal that rapid genetic divergence drives Acropora speciation.

Reef-building corals provide the structural basis for one of Earth's most spectacular and diverse but increasingly threatened ecosystems. The reef-building coral genus Acropora may have undergone substantial speciation during the Pleistocene climate and sea-level changes. Here, we aimed to evaluate the speciation history of four morphologically similar tabular Acropora species (Acropora aff. hyacinthus, A. cf. bifurcata, A. cf. cytherea, and A. cf. subulata) using an integrative approach with morphology, genetic, and reproduction methodology. Extensive morphological analyses showed that these four species are distinct and exhibited high gamete incompatibility, preventing hybridization. Furthermore, population structure and principal component analyses with SNPs (>60,000) indicated that these species were genetically distinct, and the ABBA-BABA test did not support introgression among these species. Many of their coding and noncoding RNA sequences showed high genetic variance at loci with high Fst values along the genome. Comparison of these orthologs with those of other Acropora species suggested that many of these genes are under positive selection, which could be associated with spawning time, gamete, and morphological divergence. Our findings show that the speciation of tabular Acropora occurred without hybridization, and the divergence accompanying the rapid evolution of genes in species-rich Acropora could be associated with speciation.

Yosemite toad (Anaxyrus canorus) transcriptome reveals interplay between speciation genes and adaptive introgression.

Genomes are heterogeneous during the early stages of speciation, with small 'islands' of DNA appearing to reflect strong adaptive differences, surrounded by vast seas of relative homogeneity. As species diverge, secondary contact zones between them can act as an interface and selectively filter through advantageous alleles of hybrid origin. Such introgression is another important adaptive process, one that allows beneficial mosaics of recombinant DNA ('rivers') to flow from one species into another. Although genomic islands of divergence appear to be associated with reproductive isolation, and genomic rivers form by adaptive introgression, it is unknown whether islands and rivers tend to be the same or different loci. We examined three replicate secondary contact zones for the Yosemite toad (Anaxyrus canorus) using two genomic data sets and a morphometric data set to answer the questions: (1) How predictably different are islands and rivers, both in terms of genomic location and gene function? (2) Are the adaptive genetic trait loci underlying tadpole growth and development reliably islands, rivers or neither? We found that island and river loci have significant overlap within a contact zone, suggesting that some loci are first islands, and later are predictably converted into rivers. However, gene ontology enrichment analysis showed strong overlap in gene function unique to all island loci, suggesting predictability in overall gene pathways for islands. Genome-wide association study outliers for tadpole development included LPIN3, a lipid metabolism gene potentially involved in climate change adaptation, that is island-like for all three contact zones, but also appears to be introgressing (as a river) across one zone. Taken together, our results suggest that adaptive divergence and introgression may be more complementary forces than currently appreciated.

Digest: Liolaemus lizards as an emerging model system to study hybridization as a driver of rapid radiations.

Recent perspectives on speciation genomics emphasize the pivotal role of hybridization in driving rapid radiations. The Liolaemus lizard genus displays impressive species richness with around 290 species widely distributed across southern South America. Sánchez et al. conducted a comprehensive study on the 5-million-year-old Liolaemus kingii group, which includes 14 species. The research provides new key insights to enhance our understanding of this rapid radiation, including its diversification in space and time and the consequences of hybridization in its morphological evolution and taxonomy.

Evolutionary radiation of the Eurasian Pinus species under pervasive gene flow.

Evolutionary radiation, a pivotal aspect of macroevolution, offers valuable insights into evolutionary processes. The genus Pinus is the largest genus in conifers with c . 90% of the extant species emerged in the Miocene, which signifies a case of rapid diversification. Despite this remarkable history, our understanding of the mechanisms driving radiation within this expansive genus has remained limited. Using exome capture sequencing and a fossil-calibrated phylogeny, we investigated the divergence history, niche diversification, and introgression among 13 closely related Eurasian species spanning climate zones from the tropics to the boreal Arctic. We detected complex introgression among lineages in subsection Pinus at all stages of the phylogeny. Despite this widespread gene exchange, each species maintained its genetic identity and showed clear niche differentiation. Demographic analysis unveiled distinct population histories among these species, which further influenced the nucleotide diversity and efficacy of purifying and positive selection in each species. Our findings suggest that radiation in the Eurasian pines was likely fueled by interspecific recombination and further reinforced by their adaptation to distinct environments. Our study highlights the constraints and opportunities for evolutionary change, and the expectations of future adaptation in response to environmental changes in different lineages.

Divergent dynamics of sexual and habitat isolation at the transition between stick insect populations and species.

Speciation is often viewed as a continuum along which populations diverge until they become reproductively-isolated species. However, such divergence may be heterogeneous, proceeding in fits and bursts, rather than being uniform and gradual. We show in Timema stick insects that one component of reproductive isolation evolves non-uniformly across this continuum, whereas another does not. Specifically, we use thousands of host-preference and mating trials to study habitat and sexual isolation among 42 pairs of taxa spanning a range of genomic differentiation and divergence time. We find that habitat isolation is uncoupled from genomic differentiation within species, but accumulates linearly with it between species. In contrast, sexual isolation accumulates linearly across the speciation continuum, and thus exhibits similar dynamics to morphological traits not implicated in reproductive isolation. The results show different evolutionary dynamics for different components of reproductive isolation and highlight a special relevance for species status in the process of speciation.