Evidence for ecological speciation and its alternative

The Ecology of Hybrid Incompatibilities

Ecologically mediated selection against hybrids, caused by hybrid phenotypes fitting poorly into available niches, is typically viewed as distinct from selection caused by epistatic Dobzhansky-Muller hybrid incompatibilities. Here, we show how selection against transgressive phenotypes in hybrids manifests as incompatibility. After outlining our logic, we summarize current approaches for studying ecology-based selection on hybrids. We then quantitatively review QTL-mapping studies and find traits differing between parent taxa are typically polygenic. Next, we describe how verbal models of selection on hybrids translate to phenotypic and genetic fitness landscapes, highlighting emerging approaches for detecting polygenic incompatibilities. Finally, in a synthesis of published data, we report that trait transgression-and thus possibly extrinsic hybrid incompatibility in hybrids-escalates with the phenotypic divergence between parents. We discuss conceptual implications and conclude that studying the ecological basis of hybrid incompatibility will facilitate new discoveries about mechanisms of speciation.

Climatic niche shifts and ecological sky-island dynamics in Mesoamerican montane birds

An ongoing challenge in evolutionary and ecological research focuses on testing biogeographic hypotheses for the understanding of both species' distributional patterns and of the factors influencing range limits. In this study, we described the climatic niches of Neotropical humid montane forest birds through the analysis of factors driving their evolution at inter- and intraspecific levels; and tested for differences among allopatric lineages within Aulacorhynchus, Chlorospingus, Cardellina, and Eupherusa. We employed ecological niche models (ENMs) along with an ordination approach with kernel smoothing to perform niche overlap analyses and test hypotheses of niche equivalence/similarity among lineages. In addition, we described the potential distributions of each lineage during the Late Pleistocene climate fluctuations, identifying historical range expansions, connectivity, and stability. Overall, we observed differences in environmental variables influencing climatic requirements and distributional patterns for our selected species. We detected the highest values of niche overlap mainly between Eupherusa and some Chlorospingus lineages. At both interspecific and intraspecific levels, sister lineages showed non-identical environmental niches. Our results offer weak support to a moist forest model, in which populations followed the expansion and contraction cycles of montane forests, leading to a lack of niche conservatism among lineages (they tend to occupy not identical climatic environments) throughout Mesoamerica. Therefore, historical climatic conditions may act as ecological barriers determining the distributional ranges of these species.

The interplay between climatic niche evolution, polyploidy and reproductive traits explains plant speciation in the Mediterranean Basin: a case study in Centaurium (Gentianaceae)

Speciation and diversification patterns in angiosperms are frequently shaped by niche evolution. Centaurium Hill is a Mediterranean genus with ca. 25 species, of which 60% are polyploids (tetra- and hexaploids), distributed mainly in the Mediterranean Basin and in areas with temperate and arid climates of Asia, Europe, North-Central Africa and North America. The evolutionary history of this genus has been studied using morphological, biogeographical and molecular approaches, but its climatic niche characterization and its relation with genome evolution (chromosome number and ploidy level) has not been addressed yet. Thus, this study aims to identify the role of the evolution of climatic niche, ploidy level, life cycle and floral traits in the diversification of Centaurium. Climatic niche characterization involved estimating present climate preferences using quantitative data and reconstructing ancestral niches to evaluate climatic niche shifts. The evolution of climatic niche towards selective optima determined by ploidy level (three ploidy levels) and different binary traits (polyploidy, floral size, floral display, herkogamy and life cycle) was addressed under the Ornstein-Uhlenbeck model. Chromosome number evolution was inferred using the ChromoSSE model, testing if changes are clado- or anagenetic. Chromosome number evolution and its link with cladogenesis, life cycle and floral traits was modeled on the phylogeny. The reconstruction of the ancestral niches shows that Centaurium originated in a mild climate and diversified to both humid and cold as well as to dry and warmer climates. Niche conservatism was estimated in the climatic niche of the ancestors, while the climatic niche of the current taxa experienced transitions from their ancestors' niche. Besides, the evolution of climatic niche towards multiple selective optima determined by the studied traits was supported, life cycle optima receiving the highest support. The reconstruction of chromosome number transitions shows that the rate of speciation process resulting from chromosomal changes (chromosomal cladogenesis) is similar to that of non-chromosomal cladogenesis. Additionally, dependent evolution of floral size, floral display and herkogamy with chromosome number variation was supported. In conclusion, polyploidization is a crucial process in the Mediterranean region that assisted speciation and diversification into new areas with different climates, entailing niche shifts and evolution of reproductive strategies.

Coupling of Barriers to Gene Exchange: Causes and Consequences

Coupling has emerged as a concept to describe the transition from differentiated populations to newly evolved species through the strengthening of reproductive isolation. However, the term has been used in multiple ways, and relevant processes have sometimes not been clearly distinguished. Here, we synthesize existing uses of the concept of coupling and find three main perspectives: (1) coupling as the build-up of linkage disequilibrium among loci underlying barriers to gene exchange, (2) coupling as the build-up of genome-wide linkage disequilibrium, and (3) coupling as the process generating a coincidence of distinct barrier effects. We compare and contrast these views, show the diverse processes involved and the complexity of the relationships among recombination, linkage disequilibrium, and reproductive isolation, and, finally, we emphasize how each perspective can guide new directions in speciation research. Although the importance of coupling for evolutionary divergence and speciation is well established, many theoretical and empirical questions remain unanswered.

The role of deep hybridization in fern speciation: Examples from the Thelypteridaceae

PREMISE

Hybridization is recognized as an important mechanism in fern speciation, with many allopolyploids known among congeners, as well as evidence of ancient genome duplications. Several contemporary instances of deep (intergeneric) hybridization have been noted, invariably resulting in sterile progeny. We chose the christelloid lineage of the family Thelypteridaceae, recognized for its high frequency of both intra- and intergeneric hybrids, to investigate recent hybrid speciation between deeply diverged lineages. We also seek to understand the ecological and evolutionary outcomes of resulting lineages across the landscape.

METHODS

By phasing captured reads within a phylogenomic data set of GoFlag 408 nuclear loci using HybPhaser, we investigated candidate hybrids to identify parental lineages. We estimated divergence ages by inferring a dated phylogeny using fossil calibrations with treePL. We investigated ecological niche conservatism between one confirmed intergeneric allotetraploid and its diploid progenitors using the centroid, overlap, unfilling, and expansion (COUE) framework.

RESULTS

We provide evidence for at least six instances of intergeneric hybrid speciation within the christelloid clade and estimate up to 45 million years of divergence between progenitors. The niche quantification analysis showed moderate niche overlap between an allopolyploid species and its progenitors, with significant divergence from the niche of one progenitor and conservatism to the other.

CONCLUSIONS

The examples provided here highlight the overlooked role that allopolyploidization following intergeneric hybridization may play in fern diversification and range and niche expansions. Applying this approach to other fern taxa may reveal a similar pattern of deep hybridization resulting in highly successful novel lineages.

Sex-specific life-history trait expression in hybrids of a cave- and surface-dwelling fish (Poecilia mexicana, Poeciliidae)

Evaluating the fitness of hybrids can provide important insights into genetic differences between species or diverging populations. We focused on surface- and cave-ecotypes of the widespread Atlantic molly Poecilia mexicana and raised F1 hybrids of reciprocal crosses to sexual maturity in a common-garden experiment. Hybrids were reared in a fully factorial 2 × 2 design consisting of lighting (light vs. darkness) and resource availability (high vs. low food). We quantified survival, ability to realize their full reproductive potential (i.e., completed maturation for males and 3 consecutive births for females) and essential life-history traits. Compared to the performance of pure cave and surface fish from a previous experiment, F1s had the highest death rate and the lowest proportion of fish that reached their full reproductive potential. We also uncovered an intriguing pattern of sex-specific phenotype expression, because male hybrids expressed cave molly life histories, while female hybrids expressed surface molly life histories. Our results provide evidence for strong selection against hybrids in the cave molly system, but also suggest a complex pattern of sex-specific (opposing) dominance, with certain surface molly genes being dominant in female hybrids and certain cave molly genes being dominant in male hybrids.

Local adaptation, recombination, and the fate of neopolyploids

Polyploidy is widely recognized as an important speciation mechanism because it isolates tetraploids from their diploid progenitors. Polyploidy also provides new genetic material that may facilitate adaptive evolution. However, new mutations are more likely to arise after a neopolyploid already has successfully invaded a population. Thus, the role of adaptive forces in establishing a polyploid remains unclear. One solution to this apparent paradox may lie in the capacity of polyploids to suppress recombination among preexisting locally adapted alleles. The local adaptation mechanism requires that spatially heterogeneous selection acts on multiple loci and that gene flow introduces maladapted alleles to the population where the polyploid forms. The mechanism requires neither strong genetic drift nor any intrinsic benefit of genome doubling and can accommodate any mode of gene action. A unique prediction of the mechanism is that adaptive alleles should predate polyploidization, a pattern consistent with observations from a few well-studied polyploids. The mechanism is also consistent with the coexistence of both diploid and tetraploid cytotypes, fitness heterogeneity among independently derived polyploids, and the prevalence of outcrossing among older polyploids. The local adaptation mechanism also makes novel predictions about circumstances favoring polyploid invasions that can be tested using molecular genetic or comparative approaches.

Roosting ecology drives the evolution of diverse bat landing maneuvers

Bats and birds are the only living vertebrates capable of powered flight. However, bats differ from birds in that their flight required the evolution of ascending landing maneuvers that achieve their iconic head-under-heels roosting posture. We examined the evolution of landing flight in bats and tested its association with the physical properties of roosts. Bats performed four maneuvers, each correlated with patterns of peak impact force, impulse, and roosting ecology, a critical aspect of bat biology. Our findings indicate that the common ancestor of bats performed simple, four-limbed landings, similar to extant gliding mammals, and that rotationally complex landings enhancing control over impact forces coevolved multiple times with shifts to stiff, horizontal roosts. These results suggest landing biomechanics is central to bat biology: it was critical to flight adaptation in the past, mediates roost use in the present, and may affect bats' ability to respond to deforestation in the future.

Museum genomics approach to study the taxonomy and evolution of Woolly-necked storks using historic specimens

The accessibility of genomic tools in evolutionary biology has allowed for a thorough exploration of various evolutionary processes associated with adaptation and speciation. However, genomic studies in natural systems present numerous challenges, reflecting the inherent complexities of studying organisms in their native habitats. The utilization of museum specimens for genomics research has received increased attention in recent times, facilitated by advancements in ancient DNA techniques. In this study, we have utilized a museum genomics approach to analyze historic specimens of Woolly-necked storks (Ciconia spp.) and examine their genetic composition and taxonomic status and explore the evolutionary and adaptive trajectories of populations over the years. The Woolly-necked storks are distributed in Asia and Africa with a taxonomic classification that has been a matter of ambiguity. Asian and African Woollynecks were recently recognized as different species based on their morphological differences; however, their genomic validation was lacking. In this study, we have used ∼70-year-old museum samples for whole-genome population-scale sequencing. Our study has revealed that Asian and African Woollynecks are genetically distinct, consistent with the current taxonomic classification based on morphological features. However, we also found a high genetic divergence between the Asian subspecies Ciconia episcopus neglecta and Ciconia episcopus episcopus, suggesting this classification requires a detailed examination to explore processes of ongoing speciation. Because taxonomic classification directly impacts conservation efforts, and there is evidence of declining populations of Asian Woollynecks in Southeast Asia, our results highlight that population-scale studies are urgent to determine the genetic, ecological, and phylogenetic diversity of these birds.

Comparative Genomic Analysis of the Pattern of Evolution of Male and Female Reproductive Proteins in Seed Beetles

Male seminal fluid proteins often show signs of positive selection and divergent evolution, believed to reflect male-female coevolution. Yet, our understanding of the predicted concerted evolution of seminal fluid proteins and female reproductive proteins is limited. We sequenced, assembled, and annotated the genome of two species of seed beetles allowing a comparative analysis of four closely related species of these herbivorous insects. We compare the general pattern of evolution in genes encoding seminal fluid proteins and female reproductive proteins with those in digestive protein genes and well-conserved reference genes. We found that female reproductive proteins showed an overall ratio of nonsynonymous to synonymous substitutions (ω) similar to that of conserved genes, while seminal fluid proteins and digestive proteins exhibited higher overall ω values. Further, seminal fluid proteins and digestive proteins showed a higher proportion of sites putatively under positive selection, and explicit tests showed no difference in relaxed selection between protein types. Evolutionary rate covariation analyses showed that evolutionary rates among seminal fluid proteins were on average more closely correlated with those in female reproductive proteins than with either digestive or conserved genes. Gene expression showed the expected negative covariation with ω values, except for male-biased genes where this negative relationship was reversed. In conclusion, seminal fluid proteins showed relatively rapid evolution and signs of positive selection. In contrast, female reproductive proteins evolved at a lower rate under selective constraints, on par with genes known to be well conserved. Although our findings provide support for concerted evolution of seminal fluid proteins and female reproductive proteins, they also suggest that these two classes of proteins evolve under partly distinct selective regimes.

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

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

Mate choice in the brain: species differ in how male traits 'turn on' gene expression in female brains

Mate choice plays a fundamental role in speciation, yet we know little about the molecular mechanisms that underpin this crucial decision-making process. Stickleback fish differentially adapted to limnetic and benthic habitats are reproductively isolated and females of each species use different male traits to evaluate prospective partners and reject heterospecific males. Here, we integrate behavioural data from a mate choice experiment with gene expression profiles from the brains of females actively deciding whether to mate. We find substantial gene expression variation between limnetic and benthic females, regardless of behavioural context, suggesting general divergence in constitutive gene expression patterns, corresponding to their genetic differentiation. Intriguingly, female gene co-expression modules covary with male display traits but in opposing directions for sympatric populations of the two species, suggesting male displays elicit a dynamic neurogenomic response that reflects known differences in female preferences. Furthermore, we confirm the role of numerous candidate genes previously implicated in female mate choice in other species, suggesting evolutionary tinkering with these conserved molecular processes to generate divergent mate preferences. Taken together, our study adds important new insights to our understanding of the molecular processes underlying female decision-making critical for generating sexual isolation and speciation.

Short-term fluctuating and long-term divergent selection on sympatric Monkeyflowers: insights from decade-spanning reciprocal transplants

Sympatric species often adapt to distinct microhabitats, leading to reproductive isolation and phenotypic diversity. However, temporal variation in selection may cause local maladaptation and species boundary breakdown, especially in years with large climatic events leading to episodic selection. Repeated reciprocal transplants can reveal short and long-term patterns of natural selection. To examine evolutionary trajectories of sympatric Monkeyflowers adapted to different niches, Mimulus guttatus and M. laciniatus, we performed three replicated transplants and combined them with previous experiments to leverage a dataset of five transplants spanning 10 years. We performed phenotypic selection analyses on parents and hybrids in parental habitats in Yosemite NP, CA during years of differing snowpack. If there is ecological isolation, then we predicted local adaptation and divergent phenotypic selection between habitats in line with species' differences. We found interannual fluctuations in selection, often not in predicted directions. Episodic selection due to extreme high snowpack caused a reversal of local adaptation and contributed to overall maladaptation of M. guttatus . However, a combined-year analysis detected longer-term divergent selection on flowering time, a key temporally isolating and adaptative trait, in agreement with species' differences. In conclusion, even with annual fluctuations, longer-term divergent selection may still promote species boundaries.

Ecological transcriptomics reveals stress response pathways of a ground-herb species in a waterlogging gradient of Amazonian riparian forests

Environmental stress is a fundamental facet of life and a significant driver of natural selection in the wild. Gene expression diversity may facilitate adaptation to environmental changes, without necessary genetic change, but its role in adaptive divergence remains largely understudied in Neotropical systems. In Amazonian riparian forests, species distribution is predominantly influenced by species' waterlogging tolerance. The flooding gradient delineates distinct wetland forest types, shaping habitats and species characteristics. Here we investigated the molecular basis of environmental stress response in a tropical ground-herb species (Ischnosiphon puberulus) to environmental variation in Amazonian riparian forests. We compared environmental variables and gene expression profiles from individuals collected in two forest types: Igapó and Terra firme in the Amazonian riparian forests. Predictable seasonal flooding poses a significant challenge in Igapó compared to Terra firme environments, with the former presenting higher water column height and longer flooding duration. Our findings suggest that contrasting environmental conditions related to flooding regimes are important drivers of population genetic differentiation and differential gene expression in I. puberulus. Enriched gene ontology terms highlight associations with environmental stresses, such as defence response, water transport, phosphorylation, root development, response to auxin, salicylic acid and oxidative stress. By uncovering key environmental stress response pathways conserved across populations, I. puberulus offers novel genetic insights into the molecular basis of plant reactions to environmental constraints found in flooded areas of this highly biodiverse neotropical ecosystem.

A historical stepping-stone path for an island-colonizing cactus across a submerged "bridge" archipelago

Here we use population genomic data (ddRAD-Seq) and ecological niche modeling to test biogeographic hypotheses for the divergence of the island-endemic cactus species Cereus insularis Hemsl. (Cereeae; Cactaceae) from its sister species C. fernambucensis Lem. The Cereus insularis grows in the Fernando de Noronha Islands (FNI), a Neotropical archipelago located 350 km off the Brazilian Atlantic Forest (BAF) coast. Phylogeographic reconstructions support a northward expansion by the common ancestor of C. insularis and C. fernambucensis along the mainland BAF coast, with C. insularis diverging from the widespread mainland taxon C. fernambucensis after colonizing FNI in the late Pleistocene. The morphologically distinct C. insularis is monophyletic and nested within C. fernambucensis, as expected from a progenitor-derivative speciation model. We tested alternative biogeographic and demographic hypotheses for the colonization of the FNI using Approximate Bayesian Computation. We found the greatest support for a stepping-stone path that emerged during periods of decreased sea level (the "bridge" hypothesis), in congruence with historical ecological niche modeling that shows highly suitable habitats on stepping-stone islands during glacial periods. The outlier analyses reveal signatures of selection in C. insularis, suggesting a putative role of adaptation driving rapid anagenic differentiation of this species in FNI.

Recurrent gene flow events occurred during the diversification of clownfishes of the skunk complex

Clownfish (subfamily Amphiprioninae) are an iconic group of coral reef fish that evolved a mutualistic interaction with sea anemones, which triggered the adaptive radiation of the clade. Within clownfishes, the "skunk complex" is particularly interesting. Besides ecological speciation, interspecific gene flow and hybrid speciation are thought to have shaped the evolution of the group. We investigated the mechanisms characterizing the diversification of this complex. By taking advantage of their disjunct geographical distribution, we obtained whole-genome data of sympatric and allopatric populations of the three main species of the complex (Amphiprion akallopisos, A. perideraion and A. sandaracinos). We examined population structure, genomic divergence and introgression signals and performed demographic modelling to identify the most realistic diversification scenario. We excluded scenarios of strict isolation or hybrid origin of A. sandaracinos. We discovered moderate gene flow from A. perideraion to the ancestor of A. akallopisos + A. sandaracinos and weak gene flow between the species in the Indo-Australian Archipelago throughout the diversification of the group. We identified introgressed regions in A. sandaracinos and detected in A. perideraion two large regions of high divergence from the two other species. While we found that gene flow has occurred throughout the species' diversification, we also observed that recent admixture was less pervasive than initially thought, suggesting a role of host repartition or behavioural barriers in maintaining the genetic identity of the species in sympatry.

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.

Reproductive isolation arises during laboratory adaptation to a novel hot environment

BACKGROUND

Reproductive isolation can result from adaptive processes (e.g., ecological speciation and mutation-order speciation) or stochastic processes such as "system drift" model. Ecological speciation predicts barriers to gene flow between populations from different environments, but not among replicate populations from the same environment. In contrast, reproductive isolation among populations independently adapted to the same/similar environment can arise from both mutation-order speciation or system drift.

RESULTS

In experimentally evolved populations adapting to a hot environment for over 100 generations, we find evidence for pre- and postmating reproductive isolation. On one hand, an altered lipid metabolism and cuticular hydrocarbon composition pointed to possible premating barriers between the ancestral and replicate evolved populations. On the other hand, the pronounced gene expression differences in male reproductive genes may underlie the postmating isolation among replicate evolved populations adapting to the same environment with the same standing genetic variation.

CONCLUSION

Our study confirms that replicated evolution experiments provide valuable insights into the mechanisms of speciation. The rapid emergence of the premating reproductive isolation during temperature adaptation showcases incipient ecological speciation. The potential evidence of postmating reproductive isolation among replicates gave rise to two hypotheses: (1) mutation-order speciation through a common selection on early fecundity leading to an inherent inter-locus sexual conflict; (2) system drift with genetic drift along the neutral ridges.

Interspecific selection in a diverse mycorrhizal symbiosis

Coevolution describes evolutionary change in which two or more interacting species reciprocally drive each other's evolution, potentially resulting in trait diversification and ecological speciation. Much progress has been made in analysis of its dynamics and consequences, but relatively little is understood about how coevolution works in multispecies interactions, i.e., those with diverse suites of species on one or both sides of an interaction. Interactions among plant hosts and their mutualistic ectomycorrhizal fungi (ECM) may provide an ecologically unique arena to examine the nature of selection in multispecies interactions. Using native genotypes of Monterey pine (Pinus radiata), we performed a common garden experiment at a field site that contains native stands to investigate selection from ECM fungi on pine traits. We planted seedlings from all five native populations, as well as inter-population crosses to represent intermediate phenotypes/genotypes, and measured seedling traits and ECM fungal traits to evaluate the potential for evolution in the symbiosis. We then combined field estimates of selection gradients with estimates of heritability and genetic variance-covariance matrices for multiple traits of the mutualism to determine which fungal traits drive plant fitness variation. We found evidence that certain fungal operational taxonomic units, families and species-level morphological traits by which ECM fungi acquire and transport nutrients exert selection on plant traits related to growth and allocation patterns. This work represents the first field-based, community-level study measuring multispecific coevolutionary selection in nutritional symbioses.

Evolution repeats itself in replicate long-term studies in the wild

The extent to which evolution is repeatable remains debated. Here, we study changes over time in the frequency of cryptic color-pattern morphs in 10 replicate long-term field studies of a stick insect, each spanning at least a decade (across 30 years of total data). We find predictable "up-and-down" fluctuations in stripe frequency in all populations, representing repeatable evolutionary dynamics based on standing genetic variation. A field experiment demonstrates that these fluctuations involve negative frequency-dependent natural selection (NFDS). These fluctuations rely on demographic and selective variability that pushes populations away from equilibrium, such that they can reliably move back toward it via NFDS. Last, we show that the origin of new cryptic forms is associated with multiple structural genomic variants such that which mutations arise affects evolution at larger temporal scales. Thus, evolution from existing variation is predictable and repeatable, but mutation adds complexity even for traits evolving deterministically under natural selection.

Convergent genetic adaptation of Escherichia coli in minimal media leads to pleiotropic divergence

Adaptation in an environment can either be beneficial, neutral or disadvantageous in another. To test the genetic basis of pleiotropic behaviour, we evolved six lines of E. coli independently in environments where glucose and galactose were the sole carbon sources, for 300 generations. All six lines in each environment exhibit convergent adaptation in the environment in which they were evolved. However, pleiotropic behaviour was observed in several environmental contexts, including other carbon environments. Genome sequencing reveals that mutations in global regulators rpoB and rpoC cause this pleiotropy. We report three new alleles of the rpoB gene, and one new allele of the rpoC gene. The novel rpoB alleles confer resistance to Rifampicin, and alter motility. Our results show how single nucleotide changes in the process of adaptation in minimal media can lead to wide-scale pleiotropy, resulting in changes in traits that are not under direct selection.

Trade-offs, trade-ups, and high mutational parallelism underlie microbial adaptation during extreme cycles of feast and famine

Microbes are evolutionarily robust organisms capable of rapid adaptation to complex stress, which enables them to colonize harsh environments. In nature, microbes are regularly challenged by starvation, which is a particularly complex stress because resource limitation often co-occurs with changes in pH, osmolarity, and toxin accumulation created by metabolic waste. Often overlooked are the additional complications introduced by eventual resource replenishment, as successful microbes must withstand rapid environmental shifts before swiftly capitalizing on replenished resources to avoid invasion by competing species. To understand how microbes navigate trade-offs between growth and survival, ultimately adapting to thrive in environments with extreme fluctuations, we experimentally evolved 16 Escherichia coli populations for 900 days in repeated feast/famine conditions with cycles of 100-day starvation before resource replenishment. Using longitudinal population-genomic analysis, we found that evolution in response to extreme feast/famine is characterized by narrow adaptive trajectories with high mutational parallelism and notable mutational order. Genetic reconstructions reveal that early mutations result in trade-offs for biofilm and motility but trade-ups for growth and survival, as these mutations conferred positively correlated advantages during both short-term and long-term culture. Our results demonstrate how microbes can navigate the adaptive landscapes of regularly fluctuating conditions and ultimately follow mutational trajectories that confer benefits across diverse environments.

Decoding populations in the ocean microbiome

Understanding the characteristics and structure of populations is fundamental to comprehending ecosystem processes and evolutionary adaptations. While the study of animal and plant populations has spanned a few centuries, microbial populations have been under scientific scrutiny for a considerably shorter period. In the ocean, analyzing the genetic composition of microbial populations and their adaptations to multiple niches can yield important insights into ecosystem function and the microbiome's response to global change. However, microbial populations have remained elusive to the scientific community due to the challenges associated with isolating microorganisms in the laboratory. Today, advancements in large-scale metagenomics and metatranscriptomics facilitate the investigation of populations from many uncultured microbial species directly from their habitats. The knowledge acquired thus far reveals substantial genetic diversity among various microbial species, showcasing distinct patterns of population differentiation and adaptations, and highlighting the significant role of selection in structuring populations. In the coming years, population genomics is expected to significantly increase our understanding of the architecture and functioning of the ocean microbiome, providing insights into its vulnerability or resilience in the face of ongoing global change. Video Abstract.

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.

Environmental correlates of adaptive diversification in postglacial freshwater fishes

Determining how environmental conditions contribute to divergence among populations and drive speciation is fundamental to resolving mechanisms and understanding outcomes in evolutionary biology. Postglacial freshwater fish species in the Northern Hemisphere are ideal biological systems to explore the effects of environment on diversification in morphology, ecology, and genetics (ecomorph divergences) within lakes. To date, various environmental factors have been implicated in the presence of multiple ecomorphs within particular lakes or regions. However, concerted evidence for generalizable patterns in environmental variables associated with speciation across geographical regions and across species and genera has been lacking. Here, we aimed to identify key biotic and abiotic factors associated with ecological divergence of postglacial freshwater fish species into multiple sympatric ecomorphs, focusing on species in the well-studied, widespread, and co-distributed genera Gasterosteus, Salvelinus, and Coregonus (stickleback, charr, and whitefish, respectively). We found that the presence of multiple sympatric ecomorphs tended to be associated with increasing lake surface area, maximum depth, and nutrient availability. In addition, predation, competition, and prey availability were suggested to play a role in divergence into multiple ecomorphs, but the effects of biotic factors require further study. Although we identified several environmental factors correlated with the presence of multiple ecomorphs, there were substantial data gaps across species and regions. An improved understanding of these systems may provide insight into both generalizable environmental factors involved in speciation in other systems, and potential ecological and evolutionary responses of species complexes when these variables are altered by environmental change.

Historical biogeography and evolutionary diversification of Lilium (Liliaceae): New insights from plastome phylogenomics

Here, we infer the historical biogeography and evolutionary diversification of the genus Lilium. For this purpose, we used the complete plastomes of 64 currently accepted species in the genus Lilium (14 plastomes were newly sequenced) to recover the phylogenetic backbone of the genus and a time-calibrated phylogenetic framework to estimate biogeographical history scenarios and evolutionary diversification rates of Lilium. Our results suggest that ancient climatic changes and geological tectonic activities jointly shaped the distribution range and drove evolutionary radiation of Lilium, including the Middle Miocene Climate Optimum (MMCO), the late Miocene global cooling, as well as the successive uplift of the Qinghai-Tibet Plateau (QTP) and the strengthening of the monsoon climate in East Asia during the late Miocene and the Pliocene. This case study suggests that the unique geological and climatic events in the Neogene of East Asia, in particular the uplift of QTP and the enhancement of monsoonal climate, may have played an essential role in formation of uneven distribution of plant diversity in the Northern Hemisphere.

Shades of white: The Petunia long corolla tube clade evolutionary history

Delimiting species is challenging in recently diverged species, and adaptive radiation is fundamental to understanding the evolutionary processes because it requires multiple ecological opportunities associated with adaptation to biotic and abiotic environments. The young Petunia genus (Solanaceae) is an excellent opportunity to study speciation because of its association with pollinators and unique microenvironments. This study evaluated the phylogenetic relationships among a Petunia clade species with different floral syndromes that inhabit several environments. We based our work on multiple individuals per lineage and employed nuclear and plastid phylogenetic markers and nuclear microsatellites. The phylogenetic tree revealed two main groups regarding the elevation of the distribution range, whereas microsatellites showed high polymorphism-sharing splitting lineages into three clusters. Isolation by distance, migration followed by new environment colonization, and shifts in floral syndrome were the motors for lineage differentiation, including infraspecific structuring, which suggests the need for taxonomic revision in the genus.

The effect of divergent and parallel selection on the genomic landscape of divergence

While the role of selection in divergence along the speciation continuum is theoretically well understood, defining specific signatures of selection in the genomic landscape of divergence is empirically challenging. Modelling approaches can provide insight into the potential role of selection on the emergence of a heterogenous genomic landscape of divergence. Here, we extend and apply an individual-based approach that simulates the phenotypic and genotypic distributions of two populations under a variety of selection regimes, genotype-phenotype maps, modes of migration, and genotype-environment interactions. We show that genomic islands of high differentiation and genomic valleys of similarity may respectively form under divergent and parallel selection between populations. For both types of between-population selection, negative and positive frequency-dependent selection within populations generated genomic islands of higher magnitude and genomic valleys of similarity, respectively. Divergence rates decreased under strong dominance with divergent selection, as well as in models including genotype-environment interactions under parallel selection. For both divergent and parallel selection models, divergence rate was higher under an intermittent migration regime between populations, in contrast to a constant level of migration across generations, despite an equal number of total migrants. We highlight that interpreting a particular evolutionary history from an observed genomic pattern must be done cautiously, as similar patterns may be obtained from different combinations of evolutionary processes. Modelling approaches such as ours provide an opportunity to narrow the potential routes that generate the genomic patterns of specific evolutionary histories.

Ecological determinants of Cope's rule and its inverse

Cope's rule posits that evolution gradually increases the body size in lineages. Over the last decades, two schools of thought have fueled a debate on the applicability of Cope's rule by reporting empirical evidence, respectively, for and against Cope's rule. The apparent contradictions thus documented highlight the need for a comprehensive process-based synthesis through which both positions of this debate can be understood and reconciled. Here, we use a process-based community-evolution model to investigate the eco-evolutionary emergence of Cope's rule. We report three characteristic macroevolutionary patterns, of which only two are consistent with Cope's rule. First, we find that Cope's rule applies when species interactions solely depend on relative differences in body size and the risk of lineage extinction is low. Second, in environments with higher risk of lineage extinction, the recurrent evolutionary elimination of top predators induces cyclic evolution toward larger body sizes, according to a macroevolutionary pattern we call the recurrent Cope's rule. Third, when interactions between species are determined not only by their body sizes but also by their ecological niches, the recurrent Cope's rule may get inverted, leading to cyclic evolution toward smaller body sizes. This recurrent inverse Cope's rule is characterized by highly dynamic community evolution, involving the diversification of species with large body sizes and the extinction of species with small body sizes. To our knowledge, these results provide the first theoretical foundation for reconciling the contrasting empirical evidence reported on body-size evolution.

A new species of Amolops (Amphibia, Anura, Ranidae) from Guizhou Province, China

The Torrent frogs of the genus Amolops are widely distributed in Nepal and northern India eastwards to southern China and southwards to Malaysia. The genus currently contains 84 species. Previous studies indicated underestimated species diversity in the genus. In the context, a new species occurring from the mountains in the northwestern Guizhou Province, China is found and described based on morphological comparisons and molecular phylogenetic analyses, Amolopsdafangensissp. nov. Phylogenetic analyses based on DNA sequences of the mitochondrial 16S rRNA and COI genes supported the new species as an independent lineage. The uncorrected genetic distances between the 16S rRNA and COI genes in the new species and its closest congener were 0.7% and 2.6%, respectively, which are higher than or at the same level as those among many pairs of congeners. Morphologically, the new species can be distinguished from its congeners by a combination of the following characters: body size moderate (SVL 43.2-46.8 mm in males); head length larger than head width slightly; tympanum distinct, oval; vocal sacs absent; vomerine teeth present; dorsolateral folds weak formed by series of glands; nuptial pads present on the base of finger I; heels overlapping when thighs are positioned at right angles to the body; tibiotarsal articulation reaching the level far beyond the tip of the snout when leg stretched forward.

Genetic drift promotes and recombination hinders speciation on holey fitness landscapes

Dobzhansky and Muller proposed a general mechanism through which microevolution, the substitution of alleles within populations, can cause the evolution of reproductive isolation between populations and, therefore, macroevolution. As allopatric populations diverge, many combinations of alleles differing between them have not been tested by natural selection and may thus be incompatible. Such genetic incompatibilities often cause low fitness in hybrids between species. Furthermore, the number of incompatibilities grows with the genetic distance between diverging populations. However, what determines the rate and pattern of accumulation of incompatibilities remains unclear. We investigate this question by simulating evolution on holey fitness landscapes on which genetic incompatibilities can be identified unambiguously. We find that genetic incompatibilities accumulate more slowly among genetically robust populations and identify two determinants of the accumulation rate: recombination rate and population size. In large populations with abundant genetic variation, recombination selects for increased genetic robustness and, consequently, incompatibilities accumulate more slowly. In small populations, genetic drift interferes with this process and promotes the accumulation of genetic incompatibilities. Our results suggest a novel mechanism by which genetic drift promotes and recombination hinders speciation.

Exploring coral speciation: Multiple sympatric Stylophora pistillata taxa along a divergence continuum on the Great Barrier Reef

Understanding how biodiversity originates and is maintained are fundamental challenge in evolutionary biology. Speciation is a continuous process and progression along this continuum depends on the interplay between evolutionary forces driving divergence and forces promoting genetic homogenisation. Coral reefs are broadly connected yet highly heterogeneous ecosystems, and divergence with gene flow at small spatial scales might therefore be common. Genomic studies are increasingly revealing the existence of closely related and sympatric taxa within taxonomic coral species, but the extent to which these taxa might still be exchanging genes and sharing environmental niches is unclear. In this study, we sampled extensively across diverse habitats at multiple reefs of the Great Barrier Reef (GBR) and comprehensively examined genome-wide diversity and divergence histories within and among taxa of the Stylophora pistillata species complex. S. pistillata is one of the most abundant and well-studied coral species, yet we discovered five distinct taxa, with wide geographic ranges and extensive sympatry. Demographic modelling showed that speciation events have occurred with gene flow and that taxa are at different stages along a divergence continuum. We found significant correlations between genetic divergence and specific environmental variables, suggesting that niche partitioning may have played a role in speciation and that S. pistillata taxa might be differentially adapted to different environments. Conservation actions rely on estimates of species richness, population sizes and species ranges, which are biased if divergent taxa are lumped together. As coral reefs are rapidly degrading due to climate change, our study highlights the importance of recognising evolutionarily distinct and differentially adapted coral taxa to improve conservation and restoration efforts aiming at protecting coral genetic diversity.

Anthropogenic Change and the Process of Speciation

Anthropogenic impacts on the environment alter speciation processes by affecting both geographical contexts and selection patterns on a worldwide scale. Here we review evidence of these effects. We find that human activities often generate spatial isolation between populations and thereby promote genetic divergence but also frequently cause sudden secondary contact and hybridization between diverging lineages. Human-caused environmental changes produce new ecological niches, altering selection in diverse ways that can drive diversification; but changes also often remove niches and cause extirpations. Human impacts that alter selection regimes are widespread and strong in magnitude, ranging from local changes in biotic and abiotic conditions to direct harvesting to global climate change. Altered selection, and evolutionary responses to it, impacts early-stage divergence of lineages, but does not necessarily lead toward speciation and persistence of separate species. Altogether, humans both promote and hinder speciation, although new species would form very slowly relative to anthropogenic hybridization, which can be nearly instantaneous. Speculating about the future of speciation, we highlight two key conclusions: (1) Humans will have a large influence on extinction and "despeciation" dynamics in the short term and on early-stage lineage divergence, and thus potentially speciation in the longer term, and (2) long-term monitoring combined with easily dated anthropogenic changes will improve our understanding of the processes of speciation. We can use this knowledge to preserve and restore ecosystems in ways that promote (re-)diversification, increasing future opportunities of speciation and enhancing biodiversity.

Divergent selection and climate adaptation fuel genomic differentiation between sister species of Sphagnum (peat moss)

BACKGROUND AND AIMS

New plant species can evolve through the reinforcement of reproductive isolation via local adaptation along habitat gradients. Peat mosses (Sphagnaceae) are an emerging model system for the study of evolutionary genomics and have well-documented niche differentiation among species. Recent molecular studies have demonstrated that the globally distributed species Sphagnum magellanicum is a complex of morphologically cryptic lineages that are phylogenetically and ecologically distinct. Here, we describe the architecture of genomic differentiation between two sister species in this complex known from eastern North America: the northern S. diabolicum and the largely southern S. magniae.

METHODS

We sampled plant populations from across a latitudinal gradient in eastern North America and performed whole genome and restriction-site associated DNA sequencing. These sequencing data were then analyzed computationally.

KEY RESULTS

Using sliding-window population genetic analyses we find that differentiation is concentrated within 'islands' of the genome spanning up to 400 kb that are characterized by elevated genetic divergence, suppressed recombination, reduced nucleotide diversity and increased rates of non-synonymous substitution. Sequence variants that are significantly associated with genetic structure and bioclimatic variables occur within genes that have functional enrichment for biological processes including abiotic stress response, photoperiodism and hormone-mediated signalling. Demographic modelling demonstrates that these two species diverged no more than 225 000 generations ago with secondary contact occurring where their ranges overlap.

CONCLUSIONS

We suggest that this heterogeneity of genomic differentiation is a result of linked selection and reflects the role of local adaptation to contrasting climatic zones in driving speciation. This research provides insight into the process of speciation in a group of ecologically important plants and strengthens our predictive understanding of how plant populations will respond as Earth's climate rapidly changes.

Candidate genes under selection in song sparrows co-vary with climate and body mass in support of Bergmann's Rule

Ecogeographic rules denote spatial patterns in phenotype and environment that may reflect local adaptation as well as a species' capacity to adapt to change. To identify genes underlying Bergmann's Rule, which posits that spatial correlations of body mass and temperature reflect natural selection and local adaptation in endotherms, we compare 79 genomes from nine song sparrow (Melospiza melodia) subspecies that vary ~300% in body mass (17 - 50 g). Comparing large- and smaller-bodied subspecies revealed 9 candidate genes in three genomic regions associated with body mass. Further comparisons to the five smallest subspecies endemic to California revealed eight SNPs within four of the candidate genes (GARNL3, RALGPS1, ANGPTL2, and COL15A1) associated with body mass and varying as predicted by Bergmann's Rule. Our results support the hypothesis that co-variation in environment, body mass and genotype reflect the influence of natural selection on local adaptation and a capacity for contemporary evolution in this diverse species.

Parallel evolution in Crassostrea oysters along the latitudinal gradient is associated with variation in multiple genes involved in adipogenesis

Parallel diversification provides a proper framework for studying the role of natural selection in evolution. Yet, empirical studies from ecological 'non-model' species of invertebrates are limited at the whole genome level. Here, we presented a chromosome-scale genome assembly for Crassostrea angulata and investigated the parallel genomic evolution in oysters. Specifically, we used population genomics approaches to compare two southern-northern oyster species pairs (C. angulata-C. gigas and southern-northern C. ariakensis) along the coast of China. The estimated divergence time of C. angulata and C. gigas is earlier than that of southern and northern C. ariakensis, which aligns with the overall elevated genome-wide divergence. However, the southern-northern C. ariakensis FST profile represented more extremely divergent "islands". Combined with recent reciprocal hybridization studies, we proposed that they are currently at an early stage of speciation. These two southern-northern oyster species pairs exhibited significant repeatability in patterns of genome-wide differentiation, especially in genomic regions with extremely high and low divergence. This suggested that divergent and purifying selection has contributed to the genomic parallelism between southern and northern latitudes. Top differentiated genomic regions shared in these two oyster species pairs contained candidate genes enriched for functions in energy metabolism, especially adipogenesis, which are closely related to reproductive behaviours. These genes might be good candidates for further investigation in vivo. In conclusion, our results suggest that similar divergent selection and shared genomic features could predictably transform standing genetic variation within one species pair into differences in another.

Accelerated body size evolution in upland environments is correlated with recent speciation in South American freshwater fishes

Speciation rates vary greatly among taxa and regions and are shaped by both biotic and abiotic factors. However, the relative importance and interactions of these factors are not well understood. Here we investigate the potential drivers of speciation rates in South American freshwater fishes, the most diverse continental vertebrate fauna, by examining the roles of multiple biotic and abiotic factors. We integrate a dataset on species geographic distribution, phylogenetic, morphological, climatic, and habitat data. We find that Late Neogene-Quaternary speciation events are strongly associated with body-size evolution, particularly in lineages with small body sizes that inhabit higher elevations near the continental periphery. Conversely, the effects of temperature, area, and diversity-dependence, often thought to facilitate speciation, are negligible. By evaluating multiple factors simultaneously, we demonstrate that habitat characteristics associated with elevation, as well as body size evolution, correlate with rapid speciation in South American freshwater fishes. Our study emphasizes the importance of integrative approaches that consider the interplay of biotic and abiotic factors in generating macroecological patterns of species diversity.

Genomics of plant speciation

Studies of plants have been instrumental for revealing how new species originate. For several decades, botanical research has complemented and, in some cases, challenged concepts on speciation developed via the study of other organisms while also revealing additional ways in which species can form. Now, the ability to sequence genomes at an unprecedented pace and scale has allowed biologists to settle decades-long debates and tackle other emerging challenges in speciation research. Here, we review these recent genome-enabled developments in plant speciation. We discuss complications related to identification of reproductive isolation (RI) loci using analyses of the landscape of genomic divergence and highlight the important role that structural variants have in speciation, as increasingly revealed by new sequencing technologies. Further, we review how genomics has advanced what we know of some routes to new species formation, like hybridization or whole-genome duplication, while casting doubt on others, like population bottlenecks and genetic drift. While genomics can fast-track identification of genes and mutations that confer RI, we emphasize that follow-up molecular and field experiments remain critical. Nonetheless, genomics has clarified the outsized role of ancient variants rather than new mutations, particularly early during speciation. We conclude by highlighting promising avenues of future study. These include expanding what we know so far about the role of epigenetic and structural changes during speciation, broadening the scope and taxonomic breadth of plant speciation genomics studies, and synthesizing information from extensive genomic data that have already been generated by the plant speciation community.

What Darwin could not see: island formation and historical sea levels shape genetic divergence and island biogeography in a coastal marine species

Oceanic islands play a central role in the study of evolution and island biogeography. The Galapagos Islands are one of the most studied oceanic archipelagos but research has almost exclusively focused on terrestrial organisms compared to marine species. Here we used the Galapagos bullhead shark (Heterodontus quoyi) and single nucleotide polymorphisms (SNPs) to examine evolutionary processes and their consequences for genetic divergence and island biogeography in a shallow-water marine species without larval dispersal. The sequential separation of individual islands from a central island cluster gradually established different ocean depths between islands that pose barriers to dispersal in H. quoyi. Isolation by resistance analysis suggested that ocean bathymetry and historical sea level fluctuations modified genetic connectivity. These processes resulted in at least three genetic clusters that exhibit low genetic diversity and effective population sizes that scale with island size and the level of geographic isolation. Our results exemplify that island formation and climatic cycles shape genetic divergence and biogeography of coastal marine organisms with limited dispersal comparable to terrestrial taxa. Because similar scenarios exist in oceanic islands around the globe our research provides a new perspective on marine evolution and biogeography with implications for the conservation of island biodiversity.

Allometry of Defense: Predator Shift Alters Ontogenetic Growth Patterns in an Antipredator Trait

Predation is a major factor driving prey trait diversification and promoting ecological speciation. Consequently, antipredator traits are widely studied among prey species. However, comparative studies that examine how different predators shape the ontogenetic growth of antipredator traits are scarce. In larval dragonflies, abdominal spines are effective traits against predatory fish in fish lakes, which prefer larger prey. However, defensive spines increase mortality in habitats dominated by invertebrate predators (invertebrate lakes), which prefer smaller prey. Thus, species from fish lakes may accelerate spine growth at a later body size compared to species from invertebrate lakes when growing into the preferred prey size range of predatory fish. In this study, we constructed the allometric relationship between spine length and body size and compared the inflexion point of those growth curves in five species of Leucorrhinia dragonfly larvae. We found that fish-lake Leucorrhinia species accelerated spine growth at a larger body size than congenerics from invertebrate lakes. Further, rather than extending spine length constantly through development, fish-lake species rapidly accelerated spine growth at a larger body size. This is likely to be adaptive for avoiding invertebrate predation at an early life stage, which are also present in fish lakes, though in smaller numbers. Our results highlight that comparative studies of ontogenetic patterns in antipredator traits might be essential to develop an integrated understanding of predator-prey interactions.

Beringia as a high-latitude engine of avian speciation

Beringia is a biogeographically dynamic region that extends from northeastern Asia into northwestern North America. This region has affected avian divergence and speciation in three important ways: (i) by serving as a route for intercontinental colonisation between Asia and the Americas; (ii) by cyclically splitting (and often reuniting) populations, subspecies, and species between these continents; and (iii) by providing isolated refugia through glacial cycles. The effects of these processes can be seen in taxonomic splits of shallow to increasing depths and in the presence of regional endemics. We review the taxa involved in the latter two processes (splitting-reuniting and isolation), with a focus on three research topics: avian diversity, time estimates of the generation of that diversity, and the regions within Beringia that might have been especially important. We find that these processes have generated substantial amounts of avian diversity, including 49 pairs of avian subspecies or species whose breeding distributions largely replace one another across the divide between the Old World and the New World in Beringia, and 103 avian species and subspecies endemic to this region. Among endemics, about one in three is recognised as a full biological species. Endemic taxa in the orders Charadriiformes (shorebirds, alcids, gulls, and terns) and Passeriformes (perching birds) are particularly well represented, although they show very different levels of diversity through evolutionary time. Endemic Beringian Charadriiformes have a 1.31:1 ratio of species to subspecies. In Passeriformes, endemic taxa have a 0.09:1 species-to-subspecies ratio, suggesting that passerine (and thus terrestrial) endemism might be more prone to long-term extinction in this region, although such 'losses' could occur through their being reconnected with wider continental populations during favourable climatic cycles (e.g. subspecies reintegration with other populations). Genetic evidence suggests that most Beringian avian taxa originated over the past 3 million years, confirming the importance of Quaternary processes. There seems to be no obvious clustering in their formation through time, although there might be temporal gaps with lower rates of diversity generation. For at least 62 species, taxonomically undifferentiated populations occupy this region, providing ample potential for future evolutionary diversification.

How subcultures emerge

Sympatric speciation is typically presented as a rare phenomenon, but urban subcultures frequently emerge even in the absence of geographic isolation. Is there perhaps something that culture has but biological inheritance does not that would account for this difference? We present a novel model that combines assortative interaction and multidimensional inheritance. Our computer simulations show that assortment alone can lead to the formation of cohesive clusters of individuals with low within-group and large between-group variability even in the absence of a spatial separation or disruptive natural selection. All it takes is a proportionality between the variance of inputs (cultural 'parents') and outputs (cultural 'offspring'). We argue that variability-dependent inheritance cannot be easily accomplished by genes alone, but it may be the norm, not the exception, in the transmission of culture between humans. This model explains the frequent emergence of subcultures and behavioural clustering in our species and possibly also other cultural animals.

The genetic basis of conspicuous coloration in the Guadeloupean anole: Evolution by sexual and ecological selection

Understanding how natural selection acts on the genome and contributes to the process of speciation is a primary aim of the study of evolution. Here we used natural variation in two subspecies of the Guadeloupean anole (Anolis marmoratus ssp.), from the island of Guadeloupe in the Lesser Antilles, to explore the genomic basis of adaptation and speciation in Anolis lizards. These subspecies inhabit distinct ecological environments and display marked differences in adult male color and pattern. We sequenced the complete genomes of 20 anoles, 10 from each subspecies, at 1.4× coverage. We used genome-wide scans of population differentiation, allele frequency spectrum, and linkage disequilibrium to characterize the genomic architecture within and between the subspecies. While most of the genome was undifferentiated, we observed five large divergent regions. Within these regions we identified blocks, 5 kb pairs in length, enriched for fixed single nucleotide polymorphisms. These blocks encompass 97 genes, two of which are candidate pigmentation genes. One is melanophilin (mlph), which helps transport melanosomes within melanocytes. The other is a cluster of differentiation 36 (cd36), which regulates carotenoid pigment sequestration. We used high-pressure liquid chromatography to confirm that carotenoid pigments are significantly more abundant in the conspicuous orange-pigmented skin of male A. m. marmoratus suggesting that cd36 may be regulating pigment deposition in this tissue. We identified for the first time a carotenoid gene that is a potential target of divergent sexual selection and may be contributing to the early stages of speciation in Anolis lizards.

Population genomics indicate three different modes of divergence and speciation with gene flow in the green-winged teal duck complex

The processes leading to divergence and speciation can differ broadly among taxa with different life histories. We examine these processes in a small clade of ducks with historically uncertain relationships and species limits. The green-winged teal (Anas crecca) complex is a Holarctic species of dabbling duck currently categorized as three subspecies (Anas crecca crecca, A. c. nimia, and A. c. carolinensis) with a close relative, the yellow-billed teal (Anas flavirostris) from South America. A. c. crecca and A. c. carolinensis are seasonal migrants, while the other taxa are sedentary. We examined divergence and speciation patterns in this group, determining their phylogenetic relationships and the presence and levels of gene flow among lineages using both mitochondrial and genome-wide nuclear DNA obtained from 1,393 ultraconserved element (UCE) loci. Phylogenetic relationships using nuclear DNA among these taxa showed A. c. crecca, A. c. nimia, and A. c. carolinensis clustering together to form one polytomous clade, with A. flavirostris sister to this clade. This relationship can be summarized as (crecca, nimia, carolinensis)(flavirostris). However, whole mitogenomes revealed a different phylogeny: (crecca, nimia)(carolinensis, flavirostris). The best demographic model for key pairwise comparisons supported divergence with gene flow as the probable speciation mechanism in all three contrasts (crecca-nimia, crecca-carolinensis, and carolinensis-flavirostris). Given prior work, gene flow was expected among the Holarctic taxa, but gene flow between North American carolinensis and South American flavirostris (M ∼0.1-0.4 individuals/generation), albeit low, was not expected. Three geographically oriented modes of divergence are likely involved in the diversification of this complex: heteropatric (crecca-nimia), parapatric (crecca-carolinensis), and (mostly) allopatric (carolinensis-flavirostris). Our study shows that ultraconserved elements are a powerful tool for simultaneously studying systematics and population genomics in systems with historically uncertain relationships and species limits.

Geographic-genomic and geographic-phenotypic differentiation of the Aquilegia viridiflora complex

How species diverge into different lineages is a central issue in evolutionary biology. Despite the increasing evidence indicating that such divergences do not need geographic isolation, the correlation between lineage divergence and the adaptive ecological divergence of phenotype corresponding to distribution is still unknown. In addition, gene flow has been widely detected during and through such diverging processes. We used one widely distributed Aquilegia viridiflora complex as a model system to examine genomic differentiation and corresponding phenotypic variations along geographic gradients. Our phenotypic analyses of 20 populations from northwest to northeast China identified two phenotypic groups along the geographic cline. All examined traits are distinct from each other, although a few intermediate individuals occur in their contacting regions. We further sequenced the genomes of representative individuals of each population. However, four distinct genetic lineages were detected based on nuclear genomes. In particular, we recovered numerous genetic hybrids in the contact regions of four lineages. Gene flow is widespread and continuous between four lineages but much higher between contacting lineages than geographically isolated lineages. Gene flow and natural selection might result in inconsistency between heredity and phenotype. Moreover, many genes with fast lineage-specific mutations were identified to be involved in local adaptation. Our results suggest that both geographic isolation and local selection exerted by the environment and pollinators may together create geographic distributions of phenotypic variations as well as the underlying genomic divergences in numerous lineages.

The Contribution of Epigenetics to Evolutionary Adaptation in Zingiber kawagoii Hayata (Zingiberaceae) Endemic to Taiwan

We epigenotyped 211 individuals from 17 Zingiber kawagoii populations using methylation-sensitive amplification polymorphism (MSAP) and investigated the associations of methylated (mMSAP) and unmethylated (uMSAP) loci with 16 environmental variables. Data regarding genetic variation based on amplified fragment length polymorphism (AFLP) were obtained from an earlier study. We found a significant positive correlation between genetic and epigenetic variation. Significantly higher mean mMSAP and uMSAP uH_E (unbiased expected heterozygosity: 0.223 and 0.131, respectively, p < 0.001) per locus than that estimated based on AFLP (uH_E = 0.104) were found. Genome scans detected 10 mMSAP and 9 uMSAP F_ST outliers associated with various environmental variables. A significant linear fit for 11 and 12 environmental variables with outlier mMSAP and uMSAP ordination, respectively, generated using full model redundancy analysis (RDA) was found. When conditioned on geography, partial RDA revealed that five and six environmental variables, respectively, were the most important variables influencing outlier mMSAP and uMSAP variation. We found higher genetic (average F_ST = 0.298) than epigenetic (mMSAP and uMSAP average F_ST = 0.044 and 0.106, respectively) differentiation and higher genetic isolation-by-distance (IBD) than epigenetic IBD. Strong epigenetic isolation-by-environment (IBE) was found, particularly based on the outlier data, controlling either for geography (mMSAP and uMSAP β_E = 0.128 and 0.132, respectively, p = 0.001) or for genetic structure (mMSAP and uMSAP β_E = 0.105 and 0.136, respectively, p = 0.001). Our results suggest that epigenetic variants can be substrates for natural selection linked to environmental variables and complement genetic changes in the adaptive evolution of Z. kawagoii populations.

A clue to the evolutionary history of modern East Asian flora: insights from phylogeography and diterpenoid alkaloid distribution pattern of the Spiraea japonica complex

Each subkingdom of East Asian flora (EAF) has a unique evolutionary history, but which has rarely been described based on phylogeographic studies of EAF species. The Spiraea japonica L. complex, which is widespread in East Asia (EA), has received considerable attention because of the presence of diterpenoid alkaloids (DAs). It provides a proxy for understanding the genetic diversity and DA distribution patterns of species under various environmental conditions associated with the geological background in EA. In the present study, the plastome and chloroplast/nuclear DNA of 71 populations belonging to the S. japonica complex and its congeners were sequenced, combined with DA identification, environmental analyses, and ecological niche modelling, to investigate their phylogenetic relationships, genetic and DAs distribution patterns, biogeography, and demographic dynamics. An "ampliative" S. japonica complex was put forward, comprising all species of Sect. Calospira Ser. Japonicae, of which three evolutionary units carrying their respective unique types of DAs were identified and associated with the regionalization of EAF (referring to the Hengduan Mountains, central China, and east China). Moreover, a transition belt in central China with its biogeographic significance was revealed by genetic and DA distribution patterns from the perspective of ecological adaptation. The origin and onset differentiation of the "ampliative" S. japonica complex was estimated in the early Miocene (22.01/19.44 Ma). The formation of Japanese populations (6.75 Ma) was facilitated by the land bridge, which subsequently had a fairly stable demographic history. The populations in east China have undergone a founder effect after the Last Glacial Maximum, which may have been promoted by the expansion potential of polyploidization. Overall, the in-situ origin and diversification of the "ampliative" S. japonica complex since the early Miocene is a vertical section of the formation and development of modern EAF and was shaped by the geological history of each subkingdom.

The roles of recombination and selection in shaping genomic divergence in an incipient ecological species complex

Speciation genomic studies have revealed that genomes of diverging lineages are shaped jointly by the actions of gene flow and selection. These evolutionary forces acting in concert with processes such as recombination and genome features such as gene density shape a mosaic landscape of divergence. We investigated the roles of recombination and gene density in shaping the patterns of differentiation and divergence between the cyclically parthenogenetic ecological sister-taxa, Daphnia pulicaria and Daphnia pulex. First, we assembled a phased chromosome-scale genome assembly using trio-binning for D. pulicaria and constructed a genetic map using an F2-intercross panel to understand sex-specific recombination rate heterogeneity. Finally, we used a ddRADseq data set with broad geographic sampling of D. pulicaria, D. pulex, and their hybrids to understand the patterns of genome-scale divergence and demographic parameters. Our study provides the first sex-specific estimates of recombination rates for a cyclical parthenogen, and unlike other eukaryotic species, we observed male-biased heterochiasmy in D. pulicaria, which may be related to this somewhat unique breeding mode. Additionally, regions of high gene density and recombination are generally more divergent than regions of suppressed recombination. Outlier analysis indicated that divergent genomic regions are probably driven by selection on D. pulicaria, the derived lineage colonizing a novel lake habitat. Together, our study supports a scenario of selection acting on genes related to local adaptation shaping genome-wide patterns of differentiation despite high local recombination rates in this species complex. Finally, we discuss the limitations of our data in light of demographic uncertainty.