Species delimitation with abc and other coalescent-based methods: a test of accuracy with simulations and an empirical example with lizards of the liolaemus darwinii complex (squamata: liolaemidae)

Species delimitation is a major research focus in evolutionary biology because accurate species boundaries are a prerequisite for the study of speciation. New species delimitation methods (SDMs) can accommodate nonmonophyletic species and gene tree discordance as a result of incomplete lineage sorting via the coalescent model, but do not explicitly accommodate gene flow after divergence. Approximate Bayesian computation (ABC) can incorporate gene flow and estimate other relevant parameters of the speciation process while testing alternative species delimitation hypotheses. We evaluated the accuracy of BPP, SpeDeSTEM, and ABC for delimiting species using simulated data and applied these methods to empirical data from lizards of the Liolaemus darwinii complex. Overall, BPP was the most accurate, ABC showed an intermediate accuracy, and SpeDeSTEM was the least accurate under most simulated conditions. All three SDMs showed lower accuracy when speciation occurred despite gene flow, as found in previous studies, but ABC was the method with the smallest decrease in accuracy. All three SDMs consistently supported the distinctness of southern and northern lineages within L. darwinii. These SDMs based on genetic data should be complemented with novel SDMs based on morphological and ecological data to achieve truly integrative and statistically robust approaches to species discovery.     

Predictive modeling for allopatric Strix (Strigiformes: Strigidae) owls in South America: determinants of their distributions and ecological niche‐based processes

Strix (Strigidae) is a worldwide genus of 17 owl species typical of forested habitats, including Rusty‐barred Owls (S. hylophila), Chaco Owls (S. chacoensis), and Rufous‐legged Owls (S. rufipes) in South America. These species are distributed allopatrically, but the ecological traits that determine their distributions remain largely unknown and their phylogenetic relationships are unclear. We used species distribution models (SDMs) to identify variables explaining their distribution patterns and test hypotheses about ecological divergence and conservatism based on niche overlap analysis. For Rusty‐barred Owls and Chaco Owls, climatic factors related to temperature played a major role, whereas a rainfall variable was more important for Rufous‐legged Owls. When niche overlaps were compared, accounting for regional similarities in the habitat available to each species, an ecological niche divergence process was supported for Chaco Owl‐Rusty‐barred Owl and Chaco Owl‐Rufous‐legged Owl, whereas a niche conservatism process was supported for Rusty‐barred Owl‐Rufous‐legged Owl. Different ecological requirements support current species delimitation, but they are in disagreement with the two main hypotheses currently envisaged about their phylogenetic relationships (Chaco Owls as the sister taxa of either Rufous‐legged Owls or Rusty‐barred Owls) and support a new phylogenetic hypothesis (Rufous‐legged Owls as sister taxa of Rusty‐barred Owls). Our findings suggest that speciation of Rusty‐barred Owls and Rufous‐legged Owls was a vicariant event resulting from Atlantic marine transgressions in southern South America in the Miocene, but their niche was conserved because habitat changed little in their respective ranges. In contrast, Chaco Owls diverged ecologically from the other two species as a result of their adaptations to the habitat they currently occupy. Ecological and historical approaches in biogeography can be embedded to explain distribution patterns, and results provided by SDMs can be used to infer historical and ecological processes in an integrative way.     

Ecological niche modelling as an exploratory tool for identifying species limits: an example based on Mexican muroid rodents

Niche conservatism theory suggests that recently diverged sister species share the same ecological niche. However, if the ecological niche evolves as part of the speciation process, the ecological pattern could be useful for recognizing cryptic species. In a broad sense systematists agree that the niche characters could be used for species differentiation. However, to date such characters have been ignored. We used the genetic algorithm for rule‐set production for modelling the ecological niche as a means of inferring ecological divergence in allopatric populations of muroid rodents for which taxonomic identity is uncertain. Our results show that niche differentiation is significant in most of the identified phylogroups. The differentiation is likely associated with natural evolutionary units, which can be identified by applying species concepts based on phylogenetic and ecological patterns (e.g. phylogenetic, cohesive, evolutionary). Even so, the role of the niche partition within phylogenetic reconstruction may be a limited one.     

On the scent of speciation: the chemosensory system and its role in premating isolation

Chemosensory speciation is characterized by the evolution of barriers to genetic exchange that involve chemosensory systems and chemical signals. Here, we review some representative studies documenting chemosensory speciation in an attempt to evaluate the importance and the different aspects of the process in nature and to gain insights into the genetic basis and the evolutionary mechanisms of chemosensory trait divergence. Although most studies of chemosensory speciation concern sexual isolation mediated by pheromone divergence, especially in Drosophila and moth species, other chemically based behaviours (habitat choice, pollinator attraction) can also play an important role in speciation and are likely to do so in a wide range of invertebrate and vertebrate species. Adaptive divergence of chemosensory traits in response to factors such as pollinators, hosts and conspecifics commonly drives the evolution of chemical prezygotic barriers. Although the genetic basis of chemosensory speciation remains largely unknown, genomic approaches to chemosensory gene families and to enzymes involved in biosynthetic pathways of signal compounds now provide new opportunities to dissect the genetic basis of these complex traits and of their divergence among taxa.     

The importance (or lack thereof) of niche divergence to the maintenance of a northern species complex: the case of the long‐toed salamander (Ambystoma macrodactylum Baird)

The relative importance of ecological vs. nonecological factors for the origin and maintenance of species is an open question in evolutionary biology. Young lineages – such as the distinct genetic groups that make up the ranges of many northern species – represent an opportunity to study the importance of ecological divergence during the early stages of diversification. Yet, few studies have examined the extent of niche divergence between lineages in previously glaciated regions and the role of ecology in maintaining the contact zones between them. In this study, we used tests of niche overlap in combination with ecological niche models to explore the extent of niche divergence between lineages of the long‐toed salamander (Ambystoma macrodactylum Baird) species complex and to determine whether contact zones correspond to (divergent) niche limits. We found limited evidence for niche divergence between the different long‐toed salamander lineages, substantial overlap in the predicted distribution of suitable climatic space for all lineages and range limits that are independent of niche limits. These results raise questions as to the importance of ecological divergence to the development of this widespread species complex and highlight the potential for non‐ecological factors to play a more important role in the maintenance of northern taxa.     

Correlated patterns of genetic diversity and differentiation across an avian family

Comparative studies of closely related taxa can provide insights into the evolutionary forces that shape genome evolution and the prevalence of convergent molecular evolution. We investigated patterns of genetic diversity and differentiation in stonechats (genus Saxicola), a widely distributed avian species complex with phenotypic variation in plumage, morphology and migratory behaviour, to ask whether similar genomic regions have become differentiated in independent, but closely related, taxa. We used whole‐genome pooled sequencing of 262 individuals from five taxa and found that levels of genetic diversity and divergence are strongly correlated among different stonechat taxa. We then asked whether these patterns remain correlated at deeper evolutionary scales and found that homologous genomic regions have become differentiated in stonechats and the closely related Ficedula flycatchers. Such correlation across a range of evolutionary divergence and among phylogenetically independent comparisons suggests that similar processes may be driving the differentiation of these independently evolving lineages, which in turn may be the result of intrinsic properties of particular genomic regions (e.g. areas of low recombination). Consequently, studies employing genome scans to search for areas important for reproductive isolation or adaptation should account for corresponding regions of differentiation, as these regions may not necessarily represent speciation islands or evidence of local adaptation.     

Cryptic niche conservatism among evolutionary lineages of an invasive lizard

Aim There is increasing evidence that the quality and breadth of ecological niches vary among individuals, populations, evolutionary lineages and therefore also across the range of a species. Sufficient knowledge about niche divergence among clades might thus be crucial for predicting the invasion potential of species. We tested for the first time whether evolutionary lineages of an invasive species vary in their climate niches and invasive potential. Furthermore, we tested whether lineage‐specific models show a better performance than combined models. Location Europe. Methods We used species distribution models (SDMs) based on climatic information at native and invasive ranges to test for intra‐specific niche divergence among mitochondrial DNA (mtDNA) clades of the invasive wall lizard Podarcis muralis. Using DNA barcoding, we assigned 77 invasive populations in Central Europe to eight geographically distinct evolutionary lineages. Niche similarity among lineages was assessed and the predictive power of a combination of clade‐specific SDMs was compared with a combined SDM using the pooled records of all lineages. Results We recorded eight different invasive mtDNA clades in Central Europe. The analysed clades had rather similar realized niches in their native and invasive ranges, whereas inter‐clade niche differentiation was comparatively strong. However, we found only a weak correlation between geographic origin (i.e. mtDNA clade) and invasive occurrences. Clades with narrow realized niches still became successful invaders far outside their native range, most probably due to broader fundamental niches. The combined model using data for all invasive lineages achieved a much better prediction of the invasive potential. Conclusions Our results indicate that the observed niche differentiation among evolutionary lineages is mainly driven by niche realization and not by differences in the fundamental niches. Such cryptic niche conservatism might hamper the success of clade‐specific niche modelling. Cryptic niche conservatism may in general explain the invasion success of species in areas with apparently unsuitable climate.     

Speciation with gene flow: Evidence from a complex of alpine butterflies (Coenonympha,Satyridae)

Until complete reproductive isolation is achieved, the extent of differentiation between two diverging lineages is the result of a dynamic equilibrium between genetic isolation and mixing. This is especially true for hybrid taxa, for which the degree of isolation in regard to their parental species is decisive in their capacity to rise as a new and stable entity. In this work, we explored the past and current patterns of hybridization and divergence within a complex of closely related butterflies in the genus Coenonympha in which two alpine species, C. darwiniana and C. macromma, have been shown to result from hybridization between the also alpine C. gardetta and the lowland C. arcania. By testing alternative scenarios of divergence among species, we show that gene flow has been uninterrupted throughout the speciation process, although leading to different degrees of current genetic isolation between species in contact zones depending on the pair considered. Nonetheless, at broader geographic scale, analyses reveal a clear genetic differentiation between hybrid lineages and their parental species, pointing out to an advanced stage of the hybrid speciation process. Finally, the positive correlation observed between ecological divergence and genetic isolation among these butterflies suggests a potential role for ecological drivers during their speciation processes.     

Impact of Geography and Climate on the Genetic Differentiation of the Subtropical Pine Pinus yunnanensis

Southwest China is a biodiversity hotspot characterized by complex topography, heterogeneous regional climates and rich flora. The processes and driving factors underlying this hotspot remain to be explicitly tested across taxa to gain a general understanding of the evolution of biodiversity and speciation in the region. In this study, we examined the role played by historically neutral processes, geography and environment in producing the current genetic diversity of the subtropical pine Pinus yunnanensis. We used genetic and ecological methods to investigate the patterns of genetic differentiation and ecological niche divergence across the distribution range of this species. We found both continuous genetic differentiation over the majority of its range, and discrete isolated local clusters. The discrete differentiation between two genetic groups in the west and east peripheries is consistent with niche divergence and geographical isolation of these groups. In the central area of the species’ range, population structure was shaped mainly by neutral processes and geography rather than by ecological selection. These results show that geographical and environmental factors together created stronger and more discrete genetic differentiation than isolation by distance alone, and illustrate the importance of ecological factors in forming or maintaining genetic divergence across a complex landscape. Our findings differ from other phylogenetic studies that identified the historical drainage system in the region as the primary factor shaping population structure, and highlight the heterogeneous contributions that geography and environment have made to genetic diversity among taxa in southwest China.     

Pliocene–Pleistocene ecological niche evolution shapes the phylogeography of a Mediterranean plant group

Estimating species ability to adapt to environmental changes is crucial to understand their past and future response to climate change. The Mediterranean Basin has experienced remarkable climatic changes since the Miocene, which have greatly influenced the evolution of the Mediterranean flora. Here, we examine the evolutionary history and biogeographic patterns of two sedge sister species (Carex, Cyperaceae) restricted to the western Mediterranean Basin, but with Pliocene fossil record in central Europe. In particular, we estimated the evolution of climatic niches through time and its influence in lineage differentiation. We carried out a dated phylogenetic–phylogeographic study based on seven DNA regions (nDNA and ptDNA) and fingerprinting data (AFLPs), and modelled ecological niches and species distributions for the Pliocene, Pleistocene and present. Phylogenetic and divergence time analyses revealed that both species form a monophyletic lineage originated in the late Pliocene–early Pleistocene. We detected clear genetic differentiation between both species with distinct genetic clusters in disjunct areas, indicating the predominant role of geographic barriers limiting gene flow. We found a remarkable shift in the climatic requirements between Pliocene and extant populations, although the niche seems to have been relatively conserved since the Pleistocene split of both species. This study highlights how an integrative approach combining different data sources and analyses, including fossils, allows solid and robust inferences about the evolutionary history of a plant group since the Pliocene.     

Habitat availability and gene flow influence diverging local population trajectories under scenarios of climate change: a place‐based approach

Ecological niche theory holds that species distributions are shaped by a large and complex suite of interacting factors. Species distribution models (SDMs) are increasingly used to describe species’ niches and predict the effects of future environmental change, including climate change. Currently, SDMs often fail to capture the complexity of species’ niches, resulting in predictions that are generally limited to climate‐occupancy interactions. Here, we explore the potential impact of climate change on the American pika using a replicated place‐based approach that incorporates climate, gene flow, habitat configuration, and microhabitat complexity into SDMs. Using contemporary presence–absence data from occupancy surveys, genetic data to infer connectivity between habitat patches, and 21 environmental niche variables, we built separate SDMs for pika populations inhabiting eight US National Park Service units representing the habitat and climatic breadth of the species across the western United States. We then predicted occurrence probability under current (1981–2010) and three future time periods (out to 2100). Occurrence probabilities and the relative importance of predictor variables varied widely among study areas, revealing important local‐scale differences in the realized niche of the American pika. This variation resulted in diverse and – in some cases – highly divergent future potential occupancy patterns for pikas, ranging from complete extirpation in some study areas to stable occupancy patterns in others. Habitat composition and connectivity, which are rarely incorporated in SDM projections, were influential in predicting pika occupancy in all study areas and frequently outranked climate variables. Our findings illustrate the importance of a place‐based approach to species distribution modeling that includes fine‐scale factors when assessing current and future climate impacts on species’ distributions, especially when predictions are intended to manage and conserve species of concern within individual protected areas.     

DOES NICHE DIVERGENCE ACCOMPANY ALLOPATRIC DIVERGENCE IN APHELOCOMA JAYS AS PREDICTED UNDER ECOLOGICAL SPECIATION?: INSIGHTS FROM TESTS WITH NICHE MODELS

The role of ecology in the origin of species has been the subject of long‐standing interest to evolutionary biologists. New sources of spatially explicit ecological data allow for large‐scale tests of whether speciation is associated with niche divergence or whether closely related species tend to be similar ecologically (niche conservatism). Because of the confounding effects of spatial autocorrelation of environmental variables, we generate null expectations for niche divergence for both an ecological‐niche modeling and a multivariate approach to address the question: do allopatrically distributed taxa occupy similar niches? In a classic system for the study of niche evolution—the Aphelocoma jays—we show that there is little evidence for niche divergence among Mexican Jay (A. ultramarina) lineages in the process of speciation, contrary to previous results. In contrast, Aphelocoma species that exist in partial sympatry in some regions show evidence for niche divergence. Our approach is widely applicable to the many cases of allopatric lineages in the beginning stages of speciation. These results do not support an ecological speciation model for Mexican Jay lineages because, in most cases, the allopatric environments they occupy are not significantly more divergent than expected under a null model.     

Hybridization between two parapatric ranid frog species in the northern Sierra Nevada,California, USA

Contact zones between species provide a unique opportunity to test whether taxa can hybridize or not. Cross‐breeding or hybridization between closely related taxa can promote gene flow (introgression) between species, adaptation, or even speciation. Though hybridization events may be short‐lived and difficult to detect in the field, genetic data can provide information about the level of introgression between closely related taxa. Hybridization can promote introgression between species, which may be an important evolutionary mechanism for either homogenization (reversing initial divergence between species) or reproductive isolation (potentially leading to speciation). Here, we used thousands of genetic markers from nuclear DNA to detect hybridization between two parapatric frog species (Rana boylii and Rana sierrae) in the Sierra Nevada of California. Based on principal components analysis, admixture, and analysis of heterozygosity at species diagnostic SNPs, we detected two F1 hybrid individuals in the Feather River basin, as well as a weak signal of introgression and gene flow between the frog species compared with frog populations from two other adjacent watersheds. This study provides the first documentation of hybridization and introgression between these two species, which are of conservation concern.     

Snails in the desert: Species diversification of Theba (Gastropoda: Helicidae) along the Atlantic coast of NW Africa

The spatial subdivision of species often plays a pivotal role in speciation. Across their entire range, species are rarely panmictic and crucial consequences of spatial subdivision are (1) random genetic drift including historical factors, (2) uniform selection, and (3) divergent selection. Each of these consequences may result in geographic variation and eventually reproductive isolation, but their relative importance in speciation is still unclear. In this study, we used a combination of genetic, morphological, and climatic data to obtain a comprehensive picture of differentiation among three closely related, parapatrically distributed taxa of the land snail genus Theba occurring along the Atlantic coasts of South Morocco and Western Sahara. We conducted Mantel and partial Mantel tests to relate phenotypic and genotypic variation of these species to geography and/or climate. As null hypothesis for an evolutionary scenario, we assumed nonadaptive speciation and expected a pattern of isolation by distance among taxa. Rejection of the null hypothesis would indicate isolation by environment due to adaptation. Generally, genetic drift plays an important role but is rarely considered as sole driver of speciation. It is the combination of drift and selection that predominantly drives speciation. This study, however, provides a potential example, in which nonadaptive speciation, that is, genetic drift, is apparently the main driver of shaping the diversity of Theba in NW Africa. Restriction of gene flow between populations caused by geographic isolation probably has played an important role. Climate oscillations during the Plio‐ and Pleistocene may have led to repeated ecological changes in NW Africa and disruptions of habitats promoting differentiation by geographic isolation. The inferred evolutionary scenario, however, did not fully explain the incongruence between the AFLP‐ and mtDNA‐tree topologies. This incongruence might indicate past hybridization among the studied Theba forms.     

Influence of environmental heterogeneity on genetic diversity and structure in an endemic southern Californian oak

Understanding how specific environmental factors shape gene flow while disentangling their importance relative to the effects of geographical isolation is a major question in evolutionary biology and a specific goal of landscape genetics. Here, we combine information from nuclear microsatellite markers and ecological niche modelling to study the association between climate and spatial genetic structure and variability in Engelmann oak (Quercus engelmannii), a wind-pollinated species with high potential for gene flow. We first test whether genetic diversity is associated with climatic niche suitability and stability since the Last Glacial Maximum (LGM). Second, we use causal modelling to analyse the potential influence of climatic factors (current and LGM niche suitability) and altitude in the observed patterns of genetic structure. We found that genetic diversity is negatively associated with local climatic stability since the LGM, which may be due to higher immigration rates in unstable patches during favourable climatic periods and/or temporally varying selection. Analyses of spatial genetic structure revealed the presence of three main genetic clusters, a pattern that is mainly driven by two highly differentiated populations located in the northern edge of the species distribution range. After controlling for geographic distance, causal modelling analyses showed that genetic relatedness decreases with the environmental divergence among sampling sites estimated as altitude and current and LGM niche suitability. Natural selection against nonlocal genotypes and/or asynchrony in reproductive phenology may explain this pattern. Overall, this study suggests that local environmental conditions can shape patterns of genetic structure and variability even in species with high potential for gene flow and relatively small distribution ranges.     

Temporal,spatial and ecological dynamics of speciation among amphi‐Beringian small mammals

Aim

Quaternary climate cycles played an important role in promoting diversification across the Northern Hemisphere, although details of the mechanisms driving evolutionary change are still poorly resolved. In a comparative phylogeographical framework, we investigate temporal, spatial and ecological components of evolution within a suite of Holarctic small mammals. We test a hypothesis of simultaneous divergence among multiple taxon pairs, investigating time to coalescence and demographic change for each taxon in response to a combination of climate and geography.

Location

Beringia, the nexus of the northern continents.

Methods

We used approximate Bayesian computation methods to test for simultaneous divergence among eight pairs of taxa, using cytochrome b gene sequences. We calculated coalescence times for eastern and western components of each pair and for the combined pairs, and relate dates to Quaternary climatic periodicity and combinations of environmental events and physical barriers. Population growth and expansion statistics were used to test evolutionary responses among taxa, including range shifts, persistence or periodic extirpation. Species distribution models (SDMs) for each taxon were used to predict their geographical ranges during the present interglacial, Last Glacial Maximum and previous interglacial.

Results

Multiple divergence events across Beringia were primarily coincident with extreme glacial cycles of the late Quaternary. Structure within Beringia is spatially consistent with at least three environmental barriers arising at different times: the Kolyma Uplands, Bering Strait and portions of the Bering Isthmus. Levels of divergence varied substantially, indicating evolutionary processes spanning deep and shallow time‐scales. The different demographics among taxa reflect their distinct ecological responses. SDMs predicted regional distributional changes through time and different spatial responses among taxa.

Main conclusions

Beringia predominantly constituted a dispersal corridor during the early Quaternary and a major centre of endemism in the late Quaternary. Coincident with severe glacial cycles, small mammal species were ‘caught’ in Beringia and diversified over multiple climatic phases. Relative genetic differentiation across Beringia appears to be related to ecological differences reflecting a gradual adaptation to Beringian environments through time. Some methodological constraints associated with resolving recent (late Quaternary) isolation events or drawing inferences from a single locus are discussed.     

Mechanistic niche modelling: combining physiological and spatial data to predict species' ranges

Species distribution models (SDMs) use spatial environmental data to make inferences on species' range limits and habitat suitability. Conceptually, these models aim to determine and map components of a species' ecological niche through space and time, and they have become important tools in pure and applied ecology and evolutionary biology. Most approaches are correlative in that they statistically link spatial data to species distribution records. An alternative strategy is to explicitly incorporate the mechanistic links between the functional traits of organisms and their environments into SDMs. Here, we review how the principles of biophysical ecology can be used to link spatial data to the physiological responses and constraints of organisms. This provides a mechanistic view of the fundamental niche which can then be mapped to the landscape to infer range constraints. We show how physiologically based SDMs can be developed for different organisms in different environmental contexts. Mechanistic SDMs have different strengths and weaknesses to correlative approaches, and there are many exciting and unexplored prospects for integrating the two approaches. As physiological knowledge becomes better integrated into SDMs, we will make more robust predictions of range shifts in novel or non-equilibrium contexts such as invasions, translocations, climate change and evolutionary shifts.     

Along the speciation continuum: Quantifying intrinsic and extrinsic isolating barriers across five million years of evolutionary divergence in California jewelflowers

Understanding the relative roles of intrinsic and extrinsic reproductive barriers, and their interplay within the geographic context of diverging taxa, remains an outstanding challenge in the study of speciation. We conducted a comparative analysis of reproductive isolation in California Jewelflowers (Streptanthus, s.l., Brassicaceae) by quantifying potential barriers to gene flow at multiple life history stages in 39 species pairs spanning five million years of evolutionary divergence. We quantified nine potential pre‐ and postzygotic barriers and explored patterns of reproductive isolation in relation to genetic distance. Intrinsic postzygotic isolation was initially weak, increased at intermediate genetic distances, and reached a threshold characterized by complete genetic incompatibility. Climatic niche differences were strong at shallow genetic distances, and species pairs with overlapping ranges showed slight but appreciable phenological isolation, highlighting the potential for ecological barriers to contribute to speciation. Geographic analyses suggest that speciation is not regionally allopatric in the California Jewelflowers, as recently diverged taxa occur in relatively close proximity and display substantial range overlap. Young pairs are characterized by incomplete intrinsic postzygotic isolation, suggesting that extrinsic barriers or fine‐scale spatial segregation are more important early in the divergence process than genetic incompatibilities.     

A Comparative Analysis of Genetic Differentiation across Six Shared Willow Host Species in Leaf- and Bud-Galling Sawflies

Genetic divergence and speciation in plant-feeding insects could be driven by contrasting selection pressures imposed by different plant species and taxa. While numerous examples of host-associated differentiation (HAD) have been found, the overall importance of HAD in insect diversification remains unclear, as few studies have investigated its frequency in relation to all speciation events. One promising way to infer the prevalence and repeatability of HAD is to estimate genetic differentiation in multiple insect taxa that use the same set of hosts. To this end, we measured and compared variation in mitochondrial COI and nuclear ITS2 sequences in population samples of leaf-galling Pontania and bud-galling Euura sawflies (Hymenoptera: Tenthredinidae) collected from six Salix species in two replicate locations in northern Fennoscandia. We found evidence of frequent HAD in both species complexes, as individuals from the same willow species tended to cluster together on both mitochondrial and nuclear phylogenetic trees. Although few fixed differences among the putative species were found, hierarchical AMOVAs showed that most of the genetic variation in the samples was explained by host species rather than by sampling location. Nevertheless, the levels of HAD measured across specific pairs of host species were not correlated in the two focal galler groups. Hence, our results support the hypothesis of HAD as a central force in herbivore speciation, but also indicate that evolutionary trajectories are only weakly repeatable even in temporally overlapping radiations of related insect taxa.     

Divergent selection and drift shape the genomes of two avian sister species spanning a saline–freshwater ecotone

The role of species divergence due to ecologically based divergent selection—or ecological speciation—in generating and maintaining biodiversity is a central question in evolutionary biology. Comparison of the genomes of phylogenetically related taxa spanning a selective habitat gradient enables discovery of divergent signatures of selection and thereby provides valuable insight into the role of divergent ecological selection in speciation. Tidal marsh ecosystems provide tractable opportunities for studying organisms' adaptations to selective pressures that underlie ecological divergence. Sharp environmental gradients across the saline–freshwater ecotone within tidal marshes present extreme adaptive challenges to terrestrial vertebrates. Here, we sequence 20 whole genomes of two avian sister species endemic to tidal marshes—the saltmarsh sparrow (Ammospiza caudacutus) and Nelson's sparrow (A. nelsoni)—to evaluate the influence of selective and demographic processes in shaping genome‐wide patterns of divergence. Genome‐wide divergence between these two recently diverged sister species was notably high (genome‐wide F_ST = 0.32). Against a background of high genome‐wide divergence, regions of elevated divergence were widespread throughout the genome, as opposed to focused within islands of differentiation. These patterns may be the result of genetic drift resulting from past tidal march colonization events in conjunction with divergent selection to different environments. We identified several candidate genes that exhibited elevated divergence between saltmarsh and Nelson's sparrows, including genes linked to osmotic regulation, circadian rhythm, and plumage melanism—all putative candidates linked to adaptation to tidal marsh environments. These findings provide new insights into the roles of divergent selection and genetic drift in generating and maintaining biodiversity.