Defining the speciation continuum

Chromosomal inversions from an initial ecotypic divergence drive a gradual repeated radiation of Galápagos beetles

Island faunas exhibit some of the most iconic examples where similar forms repeatedly evolve within different islands. Yet, whether these deterministic evolutionary trajectories within islands are driven by an initial, singular divergence and the subsequent exchange of individuals and adaptive genetic variation between islands remains unclear. Here, we study a gradual, repeated evolution of low-dispersive highland ecotypes from a dispersive lowland ecotype of Calosoma beetles along the island progression of the Galápagos. We show that repeated highland adaptation involved selection on multiple shared alleles within extensive chromosomal inversions that originated from an initial adaptation event on the oldest island. These highland inversions first spread through dispersal of highland individuals. Subsequent admixture with the lowland ecotype resulted in polymorphic dispersive populations from which the highland populations evolved on the youngest islands. Our findings emphasize the significance of an ancient divergence in driving repeated evolution and highlight how a mixed contribution of inter-island colonization and within-island evolution can shape parallel species communities.

A multidimensional selective landscape drives adaptive divergence between and within closely related Phlox species

Selection causes local adaptation across populations within species and simultaneously divergence between species. However, it is unclear if either the force of or the response to selection is similar across these scales. We show that natural selection drives divergence between closely related species in a pattern that is distinct from local adaptation within species. We use reciprocal transplant experiments across three species of Phlox wildflowers to characterize widespread adaptive divergence. Using provenance trials, we also find strong local adaptation between populations within a species. Comparing divergence and selection between these two scales of diversity we discover that one suite of traits predicts fitness differences between species and that an independent suite of traits predicts fitness variation within species. Selection drives divergence between species, contributing to speciation, while simultaneously favoring extensive diversity that is maintained across populations within a species. Our work demonstrates how the selection landscape is complex and multidimensional.

Embedded Complexity of Evolutionary Sequences

Multiple pathways and outcomes are common in evolutionary sequences for biological and other environmental systems due to nonlinear complexity, historical contingency, and disturbances. From any starting point, multiple evolutionary pathways are possible. From an endpoint or observed state, multiple possibilities exist for the sequence of events that created it. However, for any observed historical sequence-e.g., ecological or soil chronosequences, stratigraphic records, or lineages-only one historical sequence actually occurred. Here, a measure of the embedded complexity of historical sequences based on algebraic graph theory is introduced. Sequences are represented as system states S(t), such that S(t - 1) ≠ S(t) ≠ S(t + 1). Each sequence of N states contains nested subgraph sequences of length 2, 3, …, N - 1. The embedded complexity index (which can also be interpreted in terms of embedded information) compares the complexity (based on the spectral radius λ1) of the entire sequence to the cumulative complexity of the constituent subsequences. The spectral radius is closely linked to graph entropy, so the index also reflects information in the sequence. The analysis is also applied to ecological state-and-transition models (STM), which represent observed transitions, along with information on their causes or triggers. As historical sequences are lengthened (by the passage of time and additional transitions or by improved resolutions or new observations of historical changes), the overall complexity asymptotically approaches λ1 = 2, while the embedded complexity increases as N2.6. Four case studies are presented, representing coastal benthic community shifts determined from biostratigraphy, ecological succession on glacial forelands, vegetation community changes in longleaf pine woodlands, and habitat changes in a delta.

Heterochronic shift in gene expression leads to ontogenetic morphological divergence between two closely related polyploid species

Heterochrony-alteration to the rate or timing of development-is an important mechanism of trait differentiation associated with speciation. Heterochrony may explain the morphological divergence between two polyploid species, June sucker (Chasmistes liorus) and Utah sucker (Catostomus ardens). The larvae of both species have terminal mouths; however, as adults, June sucker and Utah sucker develop subterminal and ventral mouths, respectively. We document a difference in the timing of shape development and a corresponding change in the timing of gene expression, suggesting the distinctive mouth morphology in June suckers may result from paedomorphosis. Specifically, adult June suckers exhibit an intermediate mouth morphology between the larval (terminal) and ancestral (ventral) states. Endemic and sympatric Chasmistes/Catostomus pairs in two other lakes also are morphologically divergent, but genetically similar. These species pairs could have resulted from the differential expression of genes and corresponding divergence in trait development. Paedomorphosis may lead to adaptive diversification in Catostomids.

Diverse pathways to speciation revealed by marine snails

Speciation is a key evolutionary process that is not yet fully understood. Combining population genomic and ecological data from multiple diverging pairs of marine snails (Littorina) supports the search for speciation mechanisms. Placing pairs on a one-dimensional speciation continuum, from undifferentiated populations to species, obscured the complexity of speciation. Adding multiple axes helped to describe either speciation routes or reproductive isolation in the snails. Divergent ecological selection repeatedly generated barriers between ecotypes, but appeared less important in completing speciation while genetic incompatibilities played a key role. Chromosomal inversions contributed to genomic barriers, but with variable impact. A multidimensional (hypercube) approach supported framing of questions and identification of knowledge gaps and can be useful to understand speciation in many other systems.

Species delimitation with limited sampling: An example from rare trapdoor spider genus Cyclocosmia (Mygalomorphae, Halonoproctidae)

The outcome of species delimitation depends on many factors, including conceptual framework, study design, data availability, methodology employed and subjective decision making. Obtaining sufficient taxon sampling in endangered or rare taxa might be difficult, particularly when non-lethal tissue collection cannot be utilized. The need to avoid overexploitation of the natural populations may thus limit methodological framework available for downstream data analyses and bias the results. We test species boundaries in rare North American trapdoor spider genus Cyclocosmia Ausserer (1871) inhabiting the Southern Coastal Plain biodiversity hotspot with the use of genomic data and two multispecies coalescent model methods. We evaluate the performance of each methodology within a limited sampling framework. To mitigate the risk of species over splitting, common in taxa with highly structured populations, we subsequently implement a species validation step via genealogical diversification index (gdi), which accounts for both genetic isolation and gene flow. We delimited eight geographically restricted lineages within sampled North American Cyclocosmia, suggesting that major river drainages in the region are likely barriers to dispersal. Our results suggest that utilizing BPP in the species discovery step might be a good option for datasets comprising hundreds of loci, but fewer individuals, which may be a common scenario for rare taxa. However, we also show that such results should be validated via gdi, in order to avoid over splitting.

Hybridization, polyploidization, and morphological convergence make dozens of taxa into one chaotic genetic pool: a phylogenomic case of the Ficus erecta species complex (Moraceae)

The Ficus erecta complex, characterized by its morphological diversity and frequent interspecific overlap, shares pollinating fig wasps among several species. This attribute, coupled with its intricate phylogenetic relationships, establishes it as an exemplary model for studying speciation and evolutionary patterns. Extensive researches involving RADseq (Restriction-site associated DNA sequencing), complete chloroplast genome data, and flow cytometry methods were conducted, focusing on phylogenomic analysis, genetic structure, and ploidy detection within the complex. Significantly, the findings exposed a pronounced nuclear-cytoplasmic conflict. This evidence, together with genetic structure analysis, confirmed that hybridization within the complex is a frequent occurrence. The ploidy detection revealed widespread polyploidy, with certain species exhibiting multiple ploidy levels, including 2×, 3×, and 4×. Of particular note, only five species (F. abelii, F. erecta, F. formosana, F. tannoensis and F. vaccinioides) in the complex were proved to be monophyletic. Species such as F. gasparriniana, F. pandurata, and F. stenophylla were found to encompass multiple phylogenetically distinct lineages. This discovery, along with morphological comparisons, suggests a significant underestimation of species diversity within the complex. This study also identified F. tannoensis as an allopolyploid species originating from F. vaccinioide and F. erecta. Considering the integration of morphological, molecular systematics, and cytological evidences, it is proposed that the scope of the F. erecta complex should be expanded to the entire subsect. Frutescentiae. This would redefine the complex as a continuously evolving group comprising at least 33 taxa, characterized by blurred species boundaries, frequent hybridization and polyploidization, and ambiguous genetic differentiation.

Recent genetic, phenetic and ecological divergence across the Mesoamerican highlands: a study case with Diglossa baritula (Aves: Thraupidae)

The topographical, geological, climatic and biodiversity complexity of Mesoamerica has made it a primary research focus. The Mesoamerican highlands is a region with particularly high species richness and within-species variation. The Cinnamon-bellied Flowerpiercer, Diglossa baritula (Wagler, 1832), is a species endemic to the Mesoamerican highlands, with three allopatric subspecies currently recognized. To characterize divergence within this species, we integrated genomics, morphology, coloration and ecological niche modeling approaches, obtained from sampling individuals across the entire geographic distribution of the species. Our results revealed a clear genomic divergence between the populations to the east versus the west of the Isthmus of Tehuantepec. In contrast to the genomic results, morphology and coloration analyses showed intermediate levels of differentiation, indicating that population groups within D. baritula have probably been under similar selective pressures. Our morphology results indicated that the only sexually dimorphic morphological variable is the wing chord, with males having a longer wing chord than females. Finally, ecological data indicated that there are differences in ecological niche within D. baritula. Our data suggest that D. baritula could contain two or more incipient species at the intermediate phase of the speciation continuum. These results highlight the importance of the geographical barrier of the Isthmus of Tehuantepec and Pleistocene climatic events in driving isolation and population divergence in D. baritula. The present investigation illustrates the speciation potential of the D. baritula complex and the capacity of Mesoamerican highlands to create cryptic biodiversity and endemism.

The global speciation continuum of the cyanobacterium Microcoleus

Speciation is a continuous process driven by genetic, geographic, and ecological barriers to gene flow. It is widely investigated in multicellular eukaryotes, yet we are only beginning to comprehend the relative importance of mechanisms driving the emergence of barriers to gene flow in microbial populations. Here, we explored the diversification of the nearly ubiquitous soil cyanobacterium Microcoleus. Our dataset consisted of 291 genomes, of which 202 strains and eight herbarium specimens were sequenced for this study. We found that Microcoleus represents a global speciation continuum of at least 12 lineages, which radiated during Eocene/Oligocene aridification and exhibit varying degrees of divergence and gene flow. The lineage divergence has been driven by selection, geographical distance, and the environment. Evidence of genetic divergence and selection was widespread across the genome, but we identified regions of exceptional differentiation containing candidate genes associated with stress response and biosynthesis of secondary metabolites.

Experimental evolution of hybrid populations to identify Dobzhansky-Muller incompatibility loci

Epistatic interactions between loci that reduce fitness in interspecies hybrids, Dobzhansky-Muller incompatibilities (DMIs), contribute genetically to the inviability and infertility within hybrid populations. It remains a challenge, however, to identify the loci that contribute to DMIs as causes of reproductive isolation between species. Here, we assess through forward simulation the power of evolve-and-resequence (E&R) experimental evolution of hybrid populations to map DMI loci. We document conditions under which such a mapping strategy may be most feasible and demonstrate how mapping power is sensitive to biologically relevant parameters such as one-way versus two-way incompatibility type, selection strength, recombination rate, and dominance interactions. We also assess the influence of parameters under direct control of an experimenter, including duration of experimental evolution and number of replicate populations. We conclude that an E&R strategy for mapping DMI loci, and other cases of epistasis, can be a viable option under some circumstances for study systems with short generation times like Caenorhabditis nematodes.

Genomic differentiation and niche divergence in the Hetaerina americana (Odonata) cryptic species complex

The evolution of reproductive barriers, that is, the speciation process, implies the limitation of gene flow between populations. Different patterns of genomic differentiation throughout the speciation continuum may provide insights into the causal evolutionary forces of species divergence. In this study, we analysed a cryptic species complex of the genus Hetaerina (Odonata). This complex includes H. americana and H. calverti; however, in H. americana two highly differentiated genetic groups have been previously detected, which, we hypothesize, may correspond to different species with low morphological variation. We obtained single nucleotide polymorphism (SNP) data for 90 individuals belonging to the different taxa in the complex and carried out differentiation tests to identify genetic isolation. The results from STRUCTURE and discriminant analysis of principal components (DAPC), based on almost 5000 SNPs, confirmed the presence of three highly differentiated taxa. Also, we found FST values above 0.5 in pairwise comparisons, which indicates a considerable degree of genetic isolation among the suggested species. We also found low climatic niche overlap among all taxa, suggesting that each group occurs at specific conditions of temperature, precipitation and elevation. We propose that H. americana comprises two cryptic species, which may be reproductively isolated by ecological barriers related to niche divergence, since the morphological variation is minimal and, therefore, mechanical barriers are probably less effective compared to other related species such as H. calverti.

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.

Hybridization in agricultural weeds: A review from ecological, evolutionary, and management perspectives

Agricultural weeds frequently hybridize with each other or with related crop species. Some hybrid weeds exhibit heterosis (hybrid vigor), which may be stabilized through mechanisms like genome duplication or vegetative reproduction. Even when heterosis is not stabilized, hybridization events diversify weed gene pools and often enable adaptive introgression. Consequently, hybridization may promote weed evolution and exacerbate weed-crop competition. However, hybridization does not always increase weediness. Even when viable and fertile, hybrid weeds sometimes prove unsuccessful in crop fields. This review provides an overview of weed hybridization and its management implications. We describe intrinsic and extrinsic factors that influence hybrid fitness in agroecosystems. We also survey the rapidly growing literature on crop-weed hybridization and the link between hybridization and invasiveness. These topics are increasingly relevant in this era of genetic tools for crop improvement, intensive and simplified cropping systems, and globalized trade. The review concludes with suggested research priorities, including hybridization in the context of climate change, plant-insect interactions, and redesigned weed management programs. From a weed management perspective, hybridization is one of many reasons that researchers and land managers must diversify their weed control toolkits.

Genomic divergence and introgression between cryptic species of a widespread North American songbird

Analysis of genomic variation among related populations can sometimes reveal distinct species that were previously undescribed due to similar morphological appearances, and close examination of such cases can provide much insight regarding speciation. Genomic data can also reveal the role of reticulate evolution in differentiation and speciation. White-breasted nuthatches (Sitta carolinensis) are widely distributed North American songbirds that are currently classified as a single species but have been suspected to represent a case of cryptic speciation. Previous genetic analyses suggested four divergent groups, but it was unclear whether these represented multiple reproductively isolated species. Using extensive genomic sampling of over 350 white-breasted nuthatches from across North America and a new chromosome-level reference genome, we asked if white-breasted nuthatches are comprised of multiple species and whether introgression has occurred between divergent populations. Genomic variation of over 300,000 loci revealed four highly differentiated populations (Pacific, n = 45; Eastern, n = 23; Rocky Mountains North, n = 138; and Rocky Mountains South, n = 150) with geographic ranges that are adjacent. We observed a moderate degree of admixture between Rocky Mountain populations but only a small number of hybrids between the Rockies and the Eastern population. The rarity of hybrids together with high levels of differentiation between populations is supportive of populations having some level of reproductive isolation. Between populations, we show evidence for introgression from a divergent ghost lineage of white-breasted nuthatches into the Rocky Mountains South population, which is otherwise closely related to Rocky Mountains North. We conclude that white-breasted nuthatches are best considered at least three species and that ghost lineage introgression has contributed to differentiation between the two Rocky Mountain populations. White-breasted nuthatches provide a dramatic case of morphological similarity despite high genomic differentiation, and the varying levels of reproductive isolation among the four groups provide an example of the speciation continuum.

Genetic and karyotype divergence between parents affect clonality and sterility in hybrids

Asexual reproduction can be triggered by interspecific hybridization, but its emergence is supposedly rare, relying on exceptional combinations of suitable genomes. To examine how genomic and karyotype divergence between parental lineages affect the incidence of asexual gametogenesis, we experimentally hybridized fishes (Cobitidae) across a broad phylogenetic spectrum, assessed by whole exome data. Gametogenic pathways generally followed a continuum from sexual reproduction in hybrids between closely related evolutionary lineages to sterile or inviable crosses between distant lineages. However, most crosses resulted in a combination of sterile males and asexually reproducing females. Their gametes usually experienced problems in chromosome pairing, but females also produced a certain proportion of oocytes with premeiotically duplicated genomes, enabling their development into clonal eggs. Interspecific hybridization may thus commonly affect cell cycles in a specific way, allowing the formation of unreduced oocytes. The emergence of asexual gametogenesis appears tightly linked to hybrid sterility and constitutes an inherent part of the extended speciation continuum.

Host-Shift Speciation Proceeded with Gene Flow in Algae Covering Shells

AbstractHost shifts represent the advancement of a novel niche and often lead to speciation in symbionts. However, its mechanisms are not well understood. Here, we focused on the alga Pseudocladophora conchopheria growing on the shells of intertidal snails. Previous surveys have shown that the alga has host specificity-only attaching to the shell of Lunella correensis-but we discovered that the alga attaches to the shells of multiple sympatric snails. A genome-wide single-nucleotide polymorphism analysis (MIG-seq) was performed to determine whether host-associated speciation occurred in the algae. As a result, there was no gene flow or limited gene flow among the algae from different hosts, and some algae were genetically differentiated among hosts. In addition, the demographic estimate revealed that speciation with gene flow occurred between the algae from different hosts. Therefore, these results support the idea that host-shift speciation gradually proceeded with gene flow in the algae, providing insight into the early evolution of host shifts.

Acknowledging more biodiversity without more species

Delimiting and naming biodiversity is a vital step toward wildlife conservation and research. However, species delimitation must be consistent across biota so that the limited resources available for nature protection can be spent effectively and objectively. To date, newly discovered lineages typically are either left undescribed and thus remain unprotected or are being erroneously proposed as new species despite mixed evidence for completed speciation, in turn contributing to the emerging problem of taxonomic inflation. Inspired by recent conceptual and methodological progress, we propose a standardized workflow for species delimitation that combines phylogenetic and hybrid zone analyses of genomic datasets ("genomic taxonomy"), in which phylogeographic lineages that do not freely admix are ranked as species, while those that have remained fully genetically compatible are ranked as subspecies. In both cases, we encourage their formal taxonomic naming, diagnosis, and description to promote social awareness toward biodiversity. The use of loci throughout the genome overcomes the unreliability of widely used barcoding genes when phylogeographic patterns are complex, while the evaluation of divergence and reproductive isolation unifies the long-opposed concepts of lineage species and biological species. We suggest that a shift in conservation assessments from a single level (species) toward a two-level hierarchy (species and subspecies) will lead to a more balanced perception of biodiversity in which both intraspecific and interspecific diversity are valued and more adequately protected.

RAD-sequencing improves the genetic characterization of a threatened tree peony (Paeonia ludlowii) endemic to China: Implications for conservation

Compared with traditional genetic markers, genomic approaches have proved valuable to the conservation of endangered species. Paeonia ludlowii having rarely and pure yellow flowers, is one of the world's most famous tree peonies. However, only several wild populations remain in the Yarlung Zangbo Valley (Nyingchi and Shannan regions, Xizang) in China due to increasing anthropogenic impact on the natural habitats. We used genome-wide single nucleotide polymorphisms to elucidate the spatial pattern of genetic variation, population structure and demographic history of P. ludlowii from the fragmented region comprising the entire range of this species, aiming to provide a basis for conserving the genetic resources of this species. Unlike genetic uniformity among populations revealed in previous studies, we found low but varied levels of intra-population genetic diversity, in which lower genetic diversity was detected in the population in Shannan region compared to those in Nyingzhi region. These spatial patterns may be likely associated with different population sizes caused by micro-environment differences in these two regions. Additionally, low genetic differentiation among populations (Fst = 0.0037) were detected at the species level. This line of evidence, combined with the result of significant genetic differentiation between the two closest populations and lack of isolation by distance, suggested that shared ancestry among now remnant populations rather than contemporary genetic connectivity resulted in subtle population structure. Demographic inference suggested that P. ludlowii probably experienced a temporal history of sharp population decline during the period of Last Glacial Maximum, and a subsequent bottleneck event resulting from prehistoric human activities on the Qinghai-Tibet Plateau. All these events, together with current habitat fragment and excavation might contribute to the endangered status of P. ludlowii. Our study improved the genetic characterization of the endangered tree peony (P. ludlowii) in China, and these genetic inferences should be considered when making different in situ and ex situ conservation actions for P. ludlowii in this evolutionary hotspot region.

Speciation and repeated origins of hypertrophied lips in parallel adaptive radiations of cyprinid fish from East Africa

The evolution of convergent phenotypes is one of the most interesting phenomena of repeated adaptive radiations. Here, we examined the repeated patterns of thick-lipped or "rubberlip" phenotype of cyprinid fish of the genus Labeobarbus discovered in riverine environments of the Ethiopian Highlands, East Africa. To test the adaptive value of thickened lips, identify the ecological niche of the thick-lipped ecomorphs, and test whether these ecomorphs are the products of adaptive divergence, we studied six sympatric pairs of ecomorphs with hypertrophied lips and the normal lip structure from different riverine basins. Trophic morphology, diet, stable isotope (δ15N and δ13C) signatures, as well as mtDNA markers and genome-wide SNP variation, were analyzed. Our results show that thick-lipped ecomorphs partition trophic resources with generalized ecomorphs in only one-half of the examined sympatric pairs despite the pronounced divergence in lip structure. In these thick-lipped ecomorphs that were trophically diverged, the data on their diet along with the elevated 15N values suggest an invertivorous specialization different from the basal omnivorous-detritivouros feeding mode of the generalized ecomorphs. Genetic data confirmed an independent and parallel origin of all six lipped ecomorphs. Yet, only one of those six thick-lipped ecomorphs had a notable genetic divergence with sympatric non-lipped ecomorphs based on nuclear SNPs data (F ST = 0.21). Sympatric pairs can be sorted by combinations of phenotypic, ecological, and genetic divergence from an ecologically non-functional mouth polymorphism via ecologically functional polymorphism to a matured speciation stage via divergent evolution.

Complex Patterns of Diversification in the Gray Zone of Speciation: Model-Based Approaches Applied to Patagonian Liolaemid Lizards (Squamata: Liolaemus kingii clade)

In this study we detangled the evolutionary history of the Patagonian lizard clade Liolaemus kingii, coupling dense geographic sampling and novel computational analytical approaches. We analyzed nuclear and mitochondrial data (restriction site-associated DNA sequencing and cytochrome b) to hypothesize and evaluate species limits, phylogenetic relationships, and demographic histories. We complemented these analyses with posterior predictive simulations to assess the fit of the genomic data to the multispecies coalescent model. We also employed a novel approach to time-calibrate a phylogenetic network. Our results show several instances of mito-nuclear discordance and consistent support for a reticulated history, supporting the view that the complex evolutionary history of the kingii clade is characterized by extensive gene flow and rapid diversification events. We discuss our findings in the contexts of the "gray zone" of speciation, phylogeographic patterns in the Patagonian region, and taxonomic outcomes. [Model adequacy; multispecies coalescent; multispecies network coalescent; phylogenomics; species delimitation.].

Introgression between highly divergent fungal sister species

To understand how species evolve and adapt to changing environments, it is important to study gene flow and introgression due to their influence on speciation and radiation events. Here, we apply a novel experimental system for investigating these mechanisms using natural populations. The system is based on two fungal sister species with morphological and ecological similarities occurring in overlapping habitats. We examined introgression between these species by conducting whole genome sequencing of individuals from populations in North America and Europe. We assessed genome-wide nucleotide divergence and performed crossing experiments to study reproductive barriers. We further used ABBA-BABA statistics together with a network analysis to investigate introgression, and conducted demographic modelling to gain insight into divergence times and introgression events. The results revealed that the species are highly divergent and incompatible in vitro. Despite this, small regions of introgression were scattered throughout the genomes and one introgression event likely involves a ghost population (extant or extinct). This study demonstrates that introgression can be found among divergent species and that population histories can be studied without collections of all the populations involved. Moreover, the experimental system is shown to be a useful tool for research on reproductive isolation in natural populations.

Do sex-linked male meiotic drivers contribute to intrinsic hybrid incompatibilities? Recent empirical studies from flies and rodents

Intrinsic hybrid incompatibility is one of the important isolating barriers between species. In organisms with sex chromosomes, intrinsic hybrid incompatibility often follows two rules: Haldane's rule and large-X effects. One explanation for these two rules is that sex chromosomes are hotspots for meiotic drivers that can cause intrinsic hybrid incompatibility between geographically isolated populations. Although this hypothesis seems plausible and several empirical data are consistent with it, we are still unsure whether such mechanisms occur in nature, particularly with respect to speciation with gene flow. Here, we review empirical studies that have investigated the roles of meiotic drive in sex-chromosome evolution and speciation and propose future studies necessary for testing this hypothesis.

Speciation and development

Understanding general principles about the origin of species remains one of the foundational challenges in evolutionary biology. The genomic divergence between groups of individuals can spawn hybrid inviability and hybrid sterility, which presents a tantalizing developmental problem. Divergent developmental programs may yield either conserved or divergent phenotypes relative to ancestral traits, both of which can be responsible for reproductive isolation during the speciation process. The genetic mechanisms of developmental evolution involve cis- and trans-acting gene regulatory change, protein-protein interactions, genetic network structures, dosage, and epigenetic regulation, all of which also have roots in population genetic and molecular evolutionary processes. Toward the goal of demystifying Darwin's "mystery of mysteries," this review integrates microevolutionary concepts of genetic change with principles of organismal development, establishing explicit links between population genetic process and developmental mechanisms in the production of macroevolutionary pattern. This integration aims to establish a more unified view of speciation that binds process and mechanism.

Diversification of the African legless skinks in the subfamily Acontinae (Family Scincidae)

Cladogenic diversification is often explained by referring to climatic oscillations and geomorphic shifts that cause allopatric speciation. In this regard, southern Africa retains a high level of landscape heterogeneity in vegetation, geology, and rainfall patterns. The legless skink subfamily Acontinae occurs broadly across the southern African subcontinent and therefore provides an ideal model group for investigating biogeographic patterns associated with the region. A robust phylogenetic study of the Acontinae with comprehensive coverage and adequate sampling of each taxon has been lacking up until now, resulting in unresolved questions regarding the subfamily's biogeography and evolution. In this study, we used multi-locus genetic markers (three mitochondrial and two nuclear) with comprehensive taxon coverage (all currently recognized Acontinae species) and adequate sampling (multiple specimens for most taxa) of each taxon to infer a phylogeny for the subfamily. The phylogeny retrieved four well-supported clades in Acontias and supported the monophyly of Typhlosaurus. Following the General Lineage Concept (GLC), many long-standing phylogenetic enigmas within Acontias occidentalis and the A. kgalagadi, A. lineatus and A. meleagris species complexes, and within Typhlosaurus were resolved. Our species delimitation analyses suggest the existence of hidden taxa in the A. occidentalis, A. cregoi and A. meleagris species groups, but also suggest that some currently recognized species in the A. lineatus and A. meleagris species groups, and within Typhlosaurus, should be synonymised. We also possibly encountered "ghost introgression" in A. occidentalis. Our inferred species tree revealed a signal of gene flow, which implies possible cross-over in some groups. Fossil evidence calibration dating results showed that the divergence between Typhlosaurus and Acontias was likely influenced by cooling and increasing aridity along the southwest coast in the mid-Oligocene caused by the opening of the Drake Passage. Further cladogenesis observed in Typhlosaurus and Acontias was likely influenced by Miocene cooling, expansion of open habitat, uplifting of the eastern Great Escarpment (GE), and variation in rainfall patterns, together with the effect of the warm Agulhas Current since the early Miocene, the development of the cold Benguela Current since the late Miocene, and their co-effects. The biogeographic pattern of the Acontinae bears close resemblance to that of other herpetofauna (e.g., rain frogs and African vipers) in southern Africa.

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.

Speciation without gene-flow in hybridizing deer

Under the ecological speciation model, divergent selection acts on ecological differences between populations, gradually creating barriers to gene flow and ultimately leading to reproductive isolation. Hybridisation is part of this continuum and can both promote and inhibit the speciation process. Here, we used white-tailed (Odocoileus virginianus) and mule deer (O. hemionus) to investigate patterns of speciation in hybridizing sister species. We quantified genome-wide historical introgression and performed genome scans to look for signatures of four different selection scenarios. Despite ample modern evidence of hybridisation, we found negligible patterns of ancestral introgression and no signatures of divergence with gene flow, rather localized patterns of allopatric and balancing selection were detected across the genome. Genes under balancing selection were related to immunity, MHC and sensory perception of smell, the latter of which is consistent with deer biology. The deficiency of historical gene-flow suggests that white-tailed and mule deer were spatially separated during the glaciation cycles of the Pleistocene and genome wide differentiation accrued via genetic drift. Dobzhansky-Muller incompatibilities and selection against hybrids are hypothesised to be acting, and diversity correlations to recombination rates suggests these sister species are far along the speciation continuum.

Males and females contribute differently to the evolution of habitat segregation driven by hybridization

Costly heterospecific mating interactions, such as hybridization, select for prezygotic reproductive isolation. One of the potential traits responding to the selection arising from maladaptive hybridization is habitat preference, whose divergence results in interspecific habitat segregation. Theoretical studies have so far assumed that habitat preference is a sexually shared trait. However, male and female habitat preferences can experience different selection pressures. Here, by combining analytical and simulation approaches, we theoretically examine the evolution of sex-specific habitat preferences. Habitat segregation can have demographic consequences, potentially generating eco-evolutionary dynamics. We thus explicitly consider demography in the simulation model. We also vary the degrees of species discrimination to examine how mate choice influences the evolution of habitat preferences. Results show that both sexes can reduce hybridisation risk by settling in the habitats where abundant conspecific mates reside. However, when females can discriminate species, excess conspecific male aggregation intensifies male-male competition for mating opportunities, posing an obstacle to conspecific aggregation. Meanwhile, conspecific female aggregation attracts conspecific males, by offering the mating opportunity. Therefore, under effective species discrimination, females play a leading role in initiating habitat use divergence. Simulations typically result in either the coexistence with established habitat segregation or the extinction of one of the species. The former result is especially likely when the species differ to some extent in habitat preferences upon secondary contact. Our results disentangle the selection pressures acting on male and female habitat preferences, deepening our understanding of the evolutionary process of habitat segregation due to hybridization.

Picks in the Fabric of a Polyploidy Complex: Integrative Species Delimitation in the Tetraploid Leucanthemum Mill. (Compositae, Anthemideae) Representatives

Based on the results of a preceding species-delimitation analysis for the diploid representatives of the genus Leucanthemum (Compositae, Anthemideae), the present study aims at the elaboration of a specific and subspecific taxonomic treatment of the tetraploid members of the genus. Following an integrative taxonomic approach, species-level decisions on eight predefined morphotaxon hypotheses were based on genetic/genealogical, morphological, ecological, and geographical differentiation patterns. ddRADseq fingerprinting and SNP-based clustering revealed genetic integrity for six of the eight morphotaxa, with no clear differentiation patterns observed between the widespread L. ircutianum subsp. ircutianum and the N Spanish (Cordillera Cantábrica) L. cantabricum and the S French L. delarbrei subsp. delabrei (northern Massif Central) and L. meridionale (western Massif Central). The inclusion of differentiation patterns in morphological (leaf dissection and shape), ecological (climatological and edaphic niches), and geographical respects (pair-wise tests of sympatry vs. allopatry) together with the application of a procedural protocol for species-rank decisions (the 'Wettstein tesseract') led to the proposal of an acknowledgement of the eight predefined morphotaxon hypotheses as six species (two of them with two subspecies). Nomenclatural consequences following from these results are drawn and lead to the following new combinations: Leucanthemum delarbrei subsp. meridionale (Legrand) Oberpr., T.Ott & Vogt, comb. nov. and Leucanthemum ruscinonense (Jeanb. & Timb.-Lagr.) Oberpr., T.Ott & Vogt, comb. et stat. nov.

Integrating top-down and bottom-up approaches to understand the genetic architecture of speciation across a monkeyflower hybrid zone

Understanding the phenotypic and genetic architecture of reproductive isolation is a long-standing goal of speciation research. In several systems, large-effect loci contributing to barrier phenotypes have been characterized, but such causal connections are rarely known for more complex genetic architectures. In this study, we combine "top-down" and "bottom-up" approaches with demographic modelling toward an integrated understanding of speciation across a monkeyflower hybrid zone. Previous work suggests that pollinator visitation acts as a primary barrier to gene flow between two divergent red- and yellow-flowered ecotypes of Mimulus aurantiacus. Several candidate isolating traits and anonymous single nucleotide polymorphism loci under divergent selection have been identified, but their genomic positions remain unknown. Here, we report findings from demographic analyses that indicate this hybrid zone formed by secondary contact, but that subsequent gene flow was restricted by widespread barrier loci across the genome. Using a novel, geographic cline-based genome scan, we demonstrate that candidate barrier loci are broadly distributed across the genome, rather than mapping to one or a few "islands of speciation." Quantitative trait locus (QTL) mapping reveals that most floral traits are highly polygenic, with little evidence that QTL colocalize, indicating that most traits are genetically independent. Finally, we find little evidence that QTL and candidate barrier loci overlap, suggesting that some loci contribute to other forms of reproductive isolation. Our findings highlight the challenges of understanding the genetic architecture of reproductive isolation and reveal that barriers to gene flow other than pollinator isolation may play an important role in this system.

Genomic evidence for three distinct species in the Erebia manto complex in Central Europe (Lepidoptera, Nymphalidae)

A problem to implement conservation strategies is that in many cases recognized taxa are in fact complexes of several cryptic species. Failure to properly delineate species may lead to misplaced priorities or to inadequate conservation measures. One such species complex is the yellow-spotted ringlet Erebia manto, which comprises several phenotypically distinct lineages, whose degree of genomic isolation has so far not been assessed. Some of these lineages are geographically restricted and thus possibly represent distinct units with conservation priorities. Using several thousand nuclear genomic markers, we evaluated to which degree the bubastis lineage from the Alps and the vogesiaca lineage from the Vosges, are genetically isolated from the widespread manto lineage. Our results suggest that both lineages are genetically as strongly differentiated from manto as other taxonomically well separated sibling species in this genus from each other, supporting a delineation of bubastis and vogesiaca as independent species. Given the restricted and isolated range of vogesiaca as well as the disjunct distribution of bubastis, our findings have significant implication for future conservation efforts on these formerly cryptic species and highlight the need to investigate the genomic identity within species complexes.

Supplementary Information

The online version contains supplementary material available at 10.1007/s10592-023-01501-w.

Intraspecific variation in reproductive barriers between two closely related Arabidopsis sister species

Reproductive isolation (RI) is a critical component of speciation and varies strongly in timing and strength among different sister taxa, depending on, for example the geography of speciation and divergence time. However, these factors may also produce variation in timing and strength among populations within species. Here we tested for variation in the expression of RI among replicate population pairs between the sister taxa Arabidopsis lyrata subsp. lyrata and A. arenicola. While the former is predominantly outcrossing, the latter is predominantly selfing. We focused on intrinsic prezygotic and postzygotic RI as both species occur largely in allopatry. We assessed RI by performing within-population crosses and interspecific between-population crosses, and by raising offspring. RI was generally high between all interspecific population pairs, but it varied in timing and strength depending on population history. Prezygotic isolation was strongest between the closest-related population pair, while early postzygotic isolation was high for all other population pairs. Furthermore, the timing and strength of RI depended strongly on cross direction. Our study provides empirical support that reproductive barriers between species are highly variable among population pairs and asymmetric within population pairs, and this variation seems to follow patterns typically described across species pairs.

Single gene mutation in a plant MYB transcription factor causes a major shift in pollinator preference

Understanding the molecular basis of reproductive isolation and speciation is a key goal of evolutionary genetics. In the South American genus Petunia, the R2R3-MYB transcription factor MYB-FL regulates the biosynthesis of UV-absorbing flavonol pigments, a major determinant of pollinator preference. MYB-FL is highly expressed in the hawkmoth-pollinated P. axillaris, but independent losses of its activity in sister taxa P. secreta and P. exserta led to UV-reflective flowers and associated pollinator shifts in each lineage (bees and hummingbirds, respectively). We created a myb-fl CRISPR mutant in P. axillaris and studied the effect of this single gene on innate pollinator preference. The mutation strongly reduced the expression of the two key flavonol-related biosynthetic genes but only affected the expression of few other genes. The mutant flowers were UV reflective as expected but additionally contained low levels of visible anthocyanin pigments. Hawkmoths strongly preferred the wild-type P. axillaris over the myb-fl mutant, whereas both social and solitary bee preference depended on the level of visible color of the mutants. MYB-FL, with its specific expression pattern, small number of target genes, and key position at the nexus of flavonol and anthocyanin biosynthetic pathways, provides a striking example of evolution by single mutations of large phenotypic effect.

Ecology and the origin of non-ephemeral species

AbstractResearch over the past three decades has shown that ecology-based extrinsic reproductive barriers can rapidly arise to generate incipient species-but such barriers can also rapidly dissolve when environments change, resulting in incipient species collapse. Understanding the evolution of unconditional, "intrinsic" reproductive barriers is therefore important for understanding the longer-term buildup of biodiversity. In this article, we consider ecology's role in the evolution of intrinsic reproductive isolation. We suggest that this topic has fallen into a gap between disciplines: while evolutionary ecologists have traditionally focused on the rapid evolution of extrinsic isolation between co-occurring ecotypes, speciation geneticists studying intrinsic isolation in other taxa have devoted little attention to the ecological context in which it evolves. We argue that for evolutionary ecology to close this gap, the field will have to expand its focus beyond rapid adaptation and its traditional model systems. Synthesizing data from several subfields, we present circumstantial evidence for and against different forms of ecological adaptation as promoters of intrinsic isolation and discuss alternative forces that may be significant. We conclude by outlining complementary approaches that can better address the role of ecology in the evolution of nonephemeral reproductive barriers and, by extension, less ephemeral species.

Hyb-Seq provides new insights into the phylogeny and evolution of the Chrysanthemum zawadskii species complex in China

A species complex is an assemblage of closely related species with blurred boundaries, and from which species could arise from different speciation processes and/or a speciation continuum. Such a complex can provide an opportunity to investigate evolutionary mechanisms acting on speciation. The Chrysanthemum zawadskii species complex in China, a monophyletic group of Chrysanthemum, consists of seven species with considerable morphological variation, diverse habitats and different distribution patterns. Here, we used Hyb-Seq data to construct a well-resolved phylogeny of the C. zawadskii complex. Then, we performed comparative analyses of variation patterns in morphology, ecology and distribution to investigate the roles of geography and ecology in this complex's diversification. Lastly, we implemented divergence time estimation, species distribution modelling and ancestral area reconstruction to trace the evolutionary history of this complex. We concluded that the C. zawadskii complex originated in the Qinling-Daba mountains during the early Pliocene and then spread west and northward along the mountain ranges to northern China. During this process, geographical and ecological factors imposing different influences resulted in the current diversification and distribution patterns of this species complex, which is composed of both well-diverged species and diverging lineages on the path of speciation.

Genomic differentiation in Pacific cod using P ool‐ S eq

Patterns of genetic differentiation across the genome can provide insight into selective forces driving adaptation. We used pooled whole genome sequencing, gene annotation, and environmental covariates to evaluate patterns of genomic differentiation and to investigate mechanisms responsible for divergence among proximate Pacific cod (Gadus macrocephalus) populations from the Bering Sea and Aleutian Islands and more distant Washington Coast cod. Samples were taken from eight spawning locations, three of which were replicated to estimate consistency in allele frequency estimation. A kernel smoothing moving weighted average of relative divergence (F_ST) identified 11 genomic islands of differentiation between the Aleutian Islands and Bering Sea samples. In some islands of differentiation, there was also elevated absolute divergence (d_XY) and evidence for selection, despite proximity and potential for gene flow. Similar levels of absolute divergence (d_XY) but roughly double the relative divergence (F_ST) were observed between the distant Bering Sea and Washington Coast samples. Islands of differentiation were much smaller than the four large inversions among Atlantic cod ecotypes. Islands of differentiation between the Bering Sea and Aleutian Island were associated with SNPs from five vision system genes, which can be associated with feeding, predator avoidance, orientation, and socialization. We hypothesize that islands of differentiation between Pacific cod from the Bering Sea and Aleutian Islands provide evidence for adaptive differentiation despite gene flow in this commercially important marine species.

The strength of reproductive isolating barriers in seed plants: Insights from studies quantifying premating and postmating reproductive barriers over the past 15 years

Speciation is driven by the evolution of reproductive isolating barriers that reduce, and ultimately prevent, substantial gene flow between lineages. Despite its central role in evolutionary biology, the process can be difficult to study because it proceeds differently among groups and may occur over long timescales. Due to this complexity, we typically rely on generalizations of empirical data to describe and understand the process. Previous reviews of reproductive isolation (RI) in flowering plants have suggested that prezygotic or extrinsic barriers generally have a stronger effect on reducing gene flow compared to postzygotic or intrinsic barriers. Past conclusions have rested on relatively few empirical estimates of RI; however, RI data have become increasingly abundant over the past 15 years. We analyzed data from recent studies quantifying multiple pre- and postmating barriers in plants and compared the strengths of isolating barriers across 89 taxa pairs using standardized RI metrics. Individual prezygotic barriers were on average stronger than individual postzygotic barriers, and the total strength of prezygotic RI was approximately twice that of postzygotic RI. These findings corroborate that ecological divergence and extrinsic factors, as opposed to solely the accumulation of genetic incompatibilities, are important to speciation and the maintenance of species boundaries in plants. Despite an emphasis in the literature on asymmetric postmating and postzygotic RI, we found that prezygotic barriers acted equally asymmetrically. Overall, substantial variability in the strengths of 12 isolating barriers highlights the great diversity of mechanisms that contribute to plant diversification.

What is reproductive isolation?

Reproductive isolation (RI) is a core concept in evolutionary biology. It has been the central focus of speciation research since the modern synthesis and is the basis by which biological species are defined. Despite this, the term is used in seemingly different ways, and attempts to quantify RI have used very different approaches. After showing that the field lacks a clear definition of the term, we attempt to clarify key issues, including what RI is, how it can be quantified in principle, and how it can be measured in practice. Following other definitions with a genetic focus, we propose that RI is a quantitative measure of the effect that genetic differences between populations have on gene flow. Specifically, RI compares the flow of neutral alleles in the presence of these genetic differences to the flow without any such differences. RI is thus greater than zero when genetic differences between populations reduce the flow of neutral alleles between populations. We show how RI can be quantified in a range of scenarios. A key conclusion is that RI depends strongly on circumstances-including the spatial, temporal and genomic context-making it difficult to compare across systems. After reviewing methods for estimating RI from data, we conclude that it is difficult to measure in practice. We discuss our findings in light of the goals of speciation research and encourage the use of methods for estimating RI that integrate organismal and genetic approaches.

Cell types as species: Exploring a metaphor

The concept of "cell type," though fundamental to cell biology, is controversial. Cells have historically been classified into types based on morphology, physiology, or location. More recently, single cell transcriptomic studies have revealed fine-scale differences among cells with similar gross phenotypes. Transcriptomic snapshots of cells at various stages of differentiation, and of cells under different physiological conditions, have shown that in many cases variation is more continuous than discrete, raising questions about the relationship between cell type and cell state. Some researchers have rejected the notion of fixed types altogether. Throughout the history of discussions on cell type, cell biologists have compared the problem of defining cell type with the interminable and often contentious debate over the definition of arguably the most important concept in systematics and evolutionary biology, "species." In the last decades, systematics, like cell biology, has been transformed by the increasing availability of molecular data, and the fine-grained resolution of genetic relationships have generated new ideas about how that variation should be classified. There are numerous parallels between the two fields that make exploration of the "cell types as species" metaphor timely. These parallels begin with philosophy, with discussion of both cell types and species as being either individuals, groups, or something in between (e.g., homeostatic property clusters). In each field there are various different types of lineages that form trees or networks that can (and in some cases do) provide criteria for grouping. Developing and refining models for evolutionary divergence of species and for cell type differentiation are parallel goals of the two fields. The goal of this essay is to highlight such parallels with the hope of inspiring biologists in both fields to look for new solutions to similar problems outside of their own field.

The Warps and Wefts of a Polyploidy Complex: Integrative Species Delimitation of the Diploid Leucanthemum (Compositae, Anthemideae) Representatives

Species delimitation—owing to the paramount role of the species rank in evolutionary, ecological, and nature conservation studies—is an essential contribution of taxonomy to biodiversity research. In an ‘integrative taxonomy’ approach to species delimitation on the diploid level, we searched for evolutionary significant units (the warps and wefts) that gave rise to the polyploid complex of European ox-eye daisies (Leucanthemum; Compositae-Anthemideae). Species discovery and validation methods based on genetic, ecological, geographical, and morphometric datasets were applied to test the currently accepted diploid morpho-species, i.e., morphologically delimited species, in Leucanthemum. Novel approaches were taken in the analyses of RADseq data (consensus clustering), morphometrics of reconstructed leaf silhouettes from digitized herbarium specimens, and quantification of species-distribution overlaps. We show that 17 of the 20 Leucanthemum morpho-species are supported by genetic evidence. The taxonomic rank of the remaining three morpho-species was resolved by combining genealogic, ecologic, geographic, and morphologic data in the framework of von Wettstein’s morpho-geographical species concept. We herewith provide a methodological pipeline for the species delimitation in an ‘integrative taxonomy’ fashion using sources of evidence from genealogical, morphological, ecological, and geographical data in the philosophy of De Queiroz’s “Unified Species Concept”.

Ten years of demographic modelling of divergence and speciation in the sea

Understanding population divergence that eventually leads to speciation is essential for evolutionary biology. High species diversity in the sea was regarded as a paradox when strict allopatry was considered necessary for most speciation events because geographical barriers seemed largely absent in the sea, and many marine species have high dispersal capacities. Combining genome-wide data with demographic modelling to infer the demographic history of divergence has introduced new ways to address this classical issue. These models assume an ancestral population that splits into two subpopulations diverging according to different scenarios that allow tests for periods of gene flow. Models can also test for heterogeneities in population sizes and migration rates along the genome to account, respectively, for background selection and selection against introgressed ancestry. To investigate how barriers to gene flow arise in the sea, we compiled studies modelling the demographic history of divergence in marine organisms and extracted preferred demographic scenarios together with estimates of demographic parameters. These studies show that geographical barriers to gene flow do exist in the sea but that divergence can also occur without strict isolation. Heterogeneity of gene flow was detected in most population pairs suggesting the predominance of semipermeable barriers during divergence. We found a weak positive relationship between the fraction of the genome experiencing reduced gene flow and levels of genome-wide differentiation. Furthermore, we found that the upper bound of the 'grey zone of speciation' for our dataset extended beyond that found before, implying that gene flow between diverging taxa is possible at higher levels of divergence than previously thought. Finally, we list recommendations for further strengthening the use of demographic modelling in speciation research. These include a more balanced representation of taxa, more consistent and comprehensive modelling, clear reporting of results and simulation studies to rule out nonbiological explanations for general results.

Epigenetic induction may speed up or slow down speciation with gene flow

Speciation is less likely to occur when there is gene flow between nascent species. Natural selection can oppose gene flow and promote speciation if there is variation in ecological conditions among the nascent species' locations. Previous theory on ecological speciation with gene flow has focused primarily on the role of genetic variation in ecological traits, largely neglecting the role of nongenetic inheritance or transgenerational plasticity. Here, we build and analyze models incorporating both genetic and epigenetic inheritance, the latter representing a form of nongenetic inheritance. We investigate the rate of speciation for a population that inhabits two patches connected by migration, and find that adaptively biased epigenetic induction can speed up or slow down speciation, depending on the form of the map from genotype and epigenotype to phenotype. While adaptively relevant epigenetic variation can speed up speciation by reducing the fitness of migrants and hybrids, it can also slow down speciation. This latter effect occurs when the epialleles are able to achieve adaptation faster than the genetic alleles, thereby weakening selection on the latter.

The Build-Up of Population Genetic Divergence along the Speciation Continuum during a Recent Adaptive Radiation of Rhagoletis Flies

New species form through the evolution of genetic barriers to gene flow between previously interbreeding populations. The understanding of how speciation proceeds is hampered by our inability to follow cases of incipient speciation through time. Comparative approaches examining different diverging taxa may offer limited inferences, unless they fulfill criteria that make the comparisons relevant. Here, we test for those criteria in a recent adaptive radiation of the Rhagoletis pomonella species group (RPSG) hypothesized to have diverged in sympatry via adaptation to different host fruits. We use a large-scale population genetic survey of 1568 flies across 33 populations to: (1) detect on-going hybridization, (2) determine whether the RPSG is derived from the same proximate ancestor, and (3) examine patterns of clustering and differentiation among sympatric populations. We find that divergence of each in-group RPSG taxon is occurring under current gene flow, that the derived members are nested within the large pool of genetic variation present in hawthorn-infesting populations of R. pomonella, and that sympatric population pairs differ markedly in their degree of genotypic clustering and differentiation across loci. We conclude that the RPSG provides a particularly robust opportunity to make direct comparisons to test hypotheses about how ecological speciation proceeds despite on-going gene flow.

Adaptation and ecological speciation in seasonally varying environments at high latitudes: Drosophila virilis group

Living in high latitudes and altitudes sets specific requirements on species’ ability to forecast seasonal changes and to respond to them in an appropriate way. Adaptation into diverse environmental conditions can also lead to ecological speciation through habitat isolation or by inducing changes in traits that influence assortative mating. In this review, we explain how the unique time-measuring systems of Drosophila virilis group species have enabled the species to occupy high latitudes and how the traits involved in species reproduction and survival exhibit strong linkage with latitudinally varying photoperiodic and climatic conditions. We also describe variation in reproductive barriers between the populations of two species with overlapping distributions and show how local adaptation and the reinforcement of prezygotic barriers have created partial reproductive isolation between conspecific populations. Finally, we consider the role of species-specific chromosomal inversions and the X chromosome in the development of reproductive barriers between diverging lineages.