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.     

Diversification and asymmetrical gene flow across time and space: lineage sorting and hybridization in polytypic barking frogs

Young species complexes that are widespread across ecologically disparate regions offer important insights into the process of speciation because of their relevance to how local adaptation and gene flow influence diversification. We used mitochondrial DNA and up to 28 152 genomewide single nucleotide polymorphisms from polytypic barking frogs (Craugastor augusti complex) to infer phylogenetic relationships and test for the signature of introgressive hybridization among diverging lineages. Our phylogenetic reconstructions suggest (i) a rapid Pliocene–Pleistocene radiation that produced at least nine distinct lineages and (ii) that geographic features of the arid Central Mexican Plateau contributed to two independent northward expansions. Despite clear lineage differentiation (many private alleles and high between‐lineage F_ST scores), D‐statistic tests, which differentiate introgression from ancestral polymorphism, allowed us to identify two putative instances of reticulate gene flow. Partitioned D‐statistics provided evidence that these events occurred in the same direction between clades but at different points in time. After correcting for geographic distance, we found that lineages involved in hybrid gene flow interactions had higher levels of genetic variation than independently evolving lineages. These findings suggest that the nature of hybrid compatibility can be conserved overlong periods of evolutionary time and that hybridization between diverging lineages may contribute to standing levels of genetic variation.     

Tracing the recombination and colonization history of hybrid species in space and time

Hybrid speciation has long fascinated evolutionary biologists and laymen alike, presumably because it challenges our classical view of evolution as a ‘one‐way street’ leading to strictly tree‐like patterns of ancestry and descent. Homoploid hybrid speciation (HHS) has been a particularly interesting puzzle, as it appears to occur extremely rapidly, perhaps within less than 50 generations ( McCarthy et al. 1995 ; Buerkle et al. 2000 ). Nevertheless, HHS may sometimes involve extended or repeated periods of recombination and gene exchange between populations subject to strong divergent natural selection ( Buerkle & Rieseberg 2008 ). Thus, HHS provides a highly interesting setting for understanding the drivers and tempo of adaptive divergence and speciation in the face of gene flow ( Arnold 2006 ; Rieseberg & Willis 2007 ; Nolte & Tautz 2009). In the present issue of Molecular Ecology, Wang et al. (2011) explore a particularly challenging issue connected to HHS: they attempt to trace the colonization and recombination history of an ancient (several MYA) hybrid species, from admixture and recombination in the ancestral hybrid zone to subsequent range shifts triggered by tectonic events (uplift of the Tibetan plateau) and climatic shifts (Pleistocene ice ages). This work is important because it addresses key issues related to the origin of the standing genetic variation available for adaptive responses (e.g. to climate change) and speciation in temperate species, which are topics of great current interest ( Rieseberg et al. 2003 ; Barrett & Schluter 2008 ; de Carvalho et al. 2010 ).     

An invasive lineage of sculpins, Cottus sp. (Pisces, Teleostei) in the Rhine with new habitat adaptations has originated from hybridization between old phylogeographic groups

Fish abundance surveys in the Rhine system have shown in the past two decades that there is a rapid upriver invasion of a freshwater sculpin of the genus Cottus. These fish are found in habitats that are atypical for the known species Cottus gobio, which is confined to small cold streams within the Rhine drainage. Phylogeographic analysis based on mitochondrial haplotypes and diagnostic single nucleotide polymorphisms indicates that the invasive sculpins are hybrids between two old lineages from the River Scheldt drainage and the River Rhine drainage, although it is morphologically more similar to the Scheldt sculpins. Most importantly, however, the invasive population possesses a unique ecological potential that does not occur in either of the source populations from the Rhine or the Scheldt, which allows the colonization of new habitats that have previously been free of sculpins. Microsatellite analysis shows that the new lineage is genetically intermediate between the old lineages and that it forms a distinct genetic group across its whole expansion range. We conclude that hybridization between long separated groups has lead to the fast emergence of a new, adaptationally distinct sculpin lineage.     

Rapid formation of distinct hybrid lineages after secondary contact of two fish species (Cottus sp.)

Homoploid hybridization after secondary contact between related species can lead to mixtures of genotypes which have the potential for rapid adaptation to new environmental conditions. Here, we focus on a case where anthropogenic changes within the past 200 years have allowed the hybridization between two fish species (Cottus rhenanus and Cottus perifretum) in the Netherlands. Specifically, we address the question of the dynamics of the emergence of these hybrids and invasion of the river systems. Using a set of 81 mostly ancestry-informative SNP markers, as well as broad sample coverage in and around the area of the initial contact, we find a structured hybrid swarm with at least three distinct hybrid lineages that have emerged out of this secondary contact situation. We show that genetically coherent groups can occur at geographically distant locations, while geographically adjacent groups can be genetically different, indicating that some form of reproductive isolation between the lineages is already effective. Using a newly developed modelling approach, we test the relative influence of founding admixture, drift and migration on the allele compositions of the sampling sites. We find that the allele frequency distributions can best be explained if continued gene flow between the parental species and the hybrid lineages is invoked. Genome mapping of the invasive lineage in the Rhine shows that major chromosomal rearrangements were not involved in creating this distinct lineage. Our results show that hybridization after secondary contact can quickly lead to multiple independent new lineages that have the capacity to form hybrid species.     

Evolutionary history of freshwater sculpins, genus Cottus (Teleostei; Cottidae) and related taxa, as inferred from mitochondrial DNA phylogeny

The freshwater sculpins, genus Cottus (Teleostei; Cottidae), comprise bottom-dwelling fishes that exhibit various life-history styles, having radiated throughout Northern Hemisphere freshwater habitats. The phylogenetic relationships among Cottus and related taxa were estimated from mitochondrial DNA 12S rRNA and control region (CR) sequences, the freshwater sculpins examined falling into five lineages (A-E). Lineage A consisted of Trachidermus fasciatus and C. kazika, both having a catadromous life-history. The remaining species (lineages B-E) spawn in freshwater habitats regardless of life-history (amphidromous, lacustrine or fluvial), suggesting that the various life-history types post-dated a common ancestor of lineages B-E. Molecular clock estimates suggested a Pliocene-Pleistocene radiation (or Miocene-Pliocene from the alternative clock) of lineages B-E. In eastern Eurasia, speciation with life-history changes to amphidromous or fluvial styles has apparently occurred independently in some lineages, as a general pattern. Mitochondrial DNA CR phylogeny showed the monophyletic Baikalian cottoids (Cottoidei) to be nested within Cottus and Trachidermus, suggesting that the former ecologically and morphologically divergent cottoids may have originated from a single lineage which invaded the ancient lake.     

Fuzzy species limits in Mediterranean gorgonians (Cnidaria,Octocorallia): inferences on speciation processes

The study of the interplay between speciation and hybridization is of primary importance in evolutionary biology. Octocorals are ecologically important species whose shallow phylogenetic relationships often remain to be studied. In the Mediterranean Sea, three congeneric octocorals can be observed in sympatry: Eunicella verrucosa, Eunicella cavolini and Eunicella singularis. They display morphological differences and E. singularis hosts photosynthetic Symbiodinium, contrary to the two other species. Two nuclear sequence markers were used to study speciation and gene flow between these species, through network analysis and Approximate Bayesian Computation (ABC). Shared sequences indicated the possibility of hybridization or incomplete lineage sorting. According to ABC, a scenario of gene flow through secondary contact was the best model to explain these results. At the intraspecific level, neither geographical nor ecological isolation corresponded to distinct genetic lineages in E. cavolini. These results are discussed in the light of the potential role of ecology and genetic incompatibilities in the persistence of species limits.     

Repeated Reticulate Evolution in North American Papilio machaon Group Swallowtail Butterflies

Hybridization between distinct populations or species is increasingly recognized as an important process for generating biodiversity. However, the interaction between hybridization and speciation is complex, and the diverse evolutionary outcomes of hybridization are difficult to differentiate. Here we characterize potential hybridization in a species group of swallowtail butterflies using microsatellites, DNA sequences, and morphology, and assess whether adaptive introgression or homoploid hybrid speciation was the primary process leading to each putative hybrid lineage. Four geographically separated hybrid populations were identified in the Papilio machaon species group. One distinct mitochondrial DNA clade from P. machaon was fixed in three hybrid taxa (P. brevicauda, P. joanae, and P. m. kahli), while one hybrid swarm (P. zelicaon x machaon) exhibited this hybrid mtDNA clade as well as widespread parental mtDNA haplotypes from both parental species. Microsatellite markers and morphology showed variable admixture and intermediacy, ranging from signatures of prolonged differential introgression from the paternal species (P. polyxenes/P. zelicaon) to current gene flow with both parental species. Divergences of the hybrid lineages dated to early- to mid-Pleistocene, suggesting that repeated glaciations and subsequent range shifts of parental species, particularly P. machaon hudsonianus, facilitated initial hybridization. Although each lineage is distinct, P. joanae is the only taxon with sufficient evidence (ecological separation from parental species) to define it as a homoploid hybrid species. The repetition of hybridization in this group provides a valuable foundation for future research on hybridization, and these results emphasize the potential for hybridization to drive speciation in diverse ways.     

A phylogeographic investigation of the hybrid origin of a species of swordtail fish from Mexico

Hybrid speciation may contribute significantly to generating biodiversity, but only a few well‐documented examples for it exist so far that do not involve polyploidization as a mechanism. The swordtail fish, Xiphophorus clemenciae, shows common hallmarks of a hybrid origin and still overlaps in its current geographic distribution with its putative ancestral species (Xiphophorus hellerii and Xiphophorus maculatus). Xiphophorus clemenciae provides an ideal system for investigating the possible continued genetic interactions between a hybrid and its parental species. Here, we use microsatellite and mitochondrial markers to investigate the population structure of these species of swordtails and search for signs of recent hybridization. Individuals were sampled from 21 localities across the known range of X. clemenciae– the Isthmus of Tehuantepec (IT) Mexico, and several environmental parameters that might represent barriers to dispersal were recorded. The hybridization event that gave rise to X. clemenciae appears to be rather ancient, and a single origin is likely. We find negligible evidence for ongoing hybridization and introgression between the putative ancestral species, because they now occupy distinct ecological niches, and a common haplotype is shared by most populations of X. clemenciae. The population structure within these species shows an isolation‐by‐distance (IBD) pattern and genetic differentiation between most populations is significant and high. We infer that tectonic evolution in the Isthmus has greatly restricted gene flow between the southern and central IT populations of X. clemenciae and X. helleriii and provide preliminary information to aid in conservation management of this geographically restricted hybrid species, X. clemenciae.     

Hybridization and asymmetric introgression across a narrow zone of contact between Neotoma fuscipes and N. macrotis (Rodentia: Cricetidae)

When lineages diverge in allopatry and come into secondary contact, we have a unique opportunity to examine the degree to which they have become reproductively isolated from one another and the mechanisms that contribute to rates of interspecific gene flow. If hybridization and introgression have occurred in the past or are ongoing, examining patterns of variation in differentially inherited genetic markers can provide insight into underlying mechanisms determining interspecific reproductive interactions. We investigated genetic variation in a recently discovered contact zone between two species of woodrats (Neotoma fuscipes and N. macrotis) in central coastal California. Previous studies have found evidence of historic hybridization between these species, but an active site of sympatry had yet to be discovered. Here, we describe the first known area of sympatry between the two species and present evidence of ongoing hybridization at this site. We intensively sampled throughout the narrow area of interspecific contact and continuously into the adjacent areas where each taxon existed largely in isolation of the other. We genotyped 851 individuals sampled in a single year at 15 nuclear microsatellite loci and sequenced a portion of the cytochrome b gene in a subset of these to examine patterns of introgression. We find that the area of sympatry is less than 1 km in length, but hybridization appeared to be fairly common, with 15% of the woodrats being of hybrid origin. At least some hybrids are fertile because introgression was evident. However, introgression appeared to be asymmetric with backcrosses toward N. macrotis being more common, a pattern we suspect is due in part to directional mate choice related to differential body size. This study adds to our growing understanding of the nature of species boundaries, especially between lineages that are far along the continuum towards speciation. © 2015 The Linnean Society of London, Biological Journal of the Linnean Society, 2015, 115 , 162–172.     

A three‐way contact zone between forms of Patella rustica (Mollusca: Patellidae) in the central Mediterranean Sea

Previous studies have reported the occurrence of three differentiated mtDNA lineages within Patella rustica in the Mediterranean Sea. Two hypotheses have been proposed to explain these observations: (1) the maintenance of ancestral polymorphism within a single species; (2) the occurrence of cryptic species not identified previously. To distinguish between these hypotheses, we screened the genetic variability at nine allozyme loci, an intron from the α‐amylase gene and a mitochondrial gene for 187 individuals of P. rustica sampled from seven Mediterranean localities. Eight additional localities were screened for the last two markers to place the differentiated lineages in a clear geographic context. Our results demonstrate that the three mtDNA lineages correspond to three distinct nuclear genotype clusters and provide further details on their distribution: the cluster corresponding to the mtDNA lineage from the Atlantic and western Mediterranean extends as far as the south coast of Italy, whereas the remaining two clusters occur in sympatry in the eastern Mediterranean. One of the eastern Mediterranean clusters is highly differentiated and seems to be reproductively isolated from the codistributed form; we therefore suggest that it corresponds to a new species. The remaining two clusters are less differentiated and form a contact zone across south Italian shores. This three‐way contact zone constitutes an interesting model for the study of speciation in the marine realm. © 2010 The Linnean Society of London, Biological Journal of the Linnean Society, 2010, 100 , 154–169.     

Niche segregation and genetic structure of Campylobacter jejuni populations from wild and agricultural host species

Bacterial populations can display high levels of genetic structuring but the forces that influence this are incompletely understood. Here, by combining modelling approaches with multilocus sequence data for the zoonotic pathogen Campylobacter, we investigated how ecological factors such as niche (host) separation relate to population structure. We analysed seven housekeeping genes from published C. jejuni and C. coli isolate collections from a range of food and wild animal sources as well as abiotic environments. By reconstructing genetic structure and the patterns of ancestry, we quantified C. jejuni host association, inferred ancestral populations, investigated genetic admixture in different hosts and determined the host origin of recombinant C. jejuni alleles found in hybrid C. coli lineages. Phylogenetically distinct C. jejuni lineages were associated with phylogenetically distinct wild birds. However, in the farm environment, phylogenetically distant host animals shared several C. jejuni lineages that could not be segregated according to host origin using these analyses. Furthermore, of the introgressed C. jejuni alleles found in C. coli lineages, 73% were attributed to genotypes associated with food animals. Our results are consistent with an evolutionary scenario where distinct Campylobacter lineages are associated with different host species but the ecological factors that maintain this are different in domestic animals such that phylogenetically distant animals can harbour closely related strains.     

Diversification along a benthic to pelagic gradient contributes to fish diversity in the world's largest lake (Lake Baikal,Russia)

Insights into the generation of diversity in both plants and animals have relied heavily on studying speciation in adaptive radiations. Russia's Lake Baikal has facilitated a putative adaptive radiation of cottid fishes (sculpins), some of which are highly specialized to inhabit novel niches created by the lake's unique geology and ecology. Here, we test evolutionary relationships and novel morphological adaptation in a piece of this radiation: the Baikal cottid genus, Cottocomephorus, a morphologically derived benthopelagic genus of three described species. We used a combination of mitochondrial DNA and restriction site associated DNA sequencing from all Cottocomephorus species. Analysis of mitochondrial cytochrome b haplotypes was only able to two resolve two lineages: C. grewingkii and C. comephoroides/inermis. Phylogenetic inference, principal component analysis, and faststructure of genome-wide SNPs uncovered three lineages within Cottocomephorus: C. comephoroides, C. inermis and C. grewingkii. We found recent divergence and admixture between C. comephoroides and C. inermis and deep divergence between these two species and C. grewingkii. Contrasting other fish radiations, we found no evidence of ancient hybridization among Cottocomephorus species. Digital morphology revealed highly derived pelagic phenotypes that reflect divergence by specialization to the benthopelagic niche in Cottocomephorus. Among Cottocomephorus species, we found evidence of ongoing adaptation to the pelagic zone. This pattern highlights the importance of speciation along a benthic-pelagic gradient seen in Cottocomephorus and across other adaptive fish radiations.     

Repeatability,ephemerality and inconvenient truths in the speciation process

Everyone appreciates the happy fiction that species conform to the simple theoretical convenience of a single panmictic population. In speciation genetics, a further standard simplification is that it is only those genetic differences that are fixed between diverging populations that need concern us in order to understand the accumulation of intrinsic barriers to reproduction. To a first approximation, of course, both of these assumptions are appropriate and theory based on them provides compelling insights into diverse evolutionary phenomena (Orr & Turelli 2001 ). But what else can we learn about the begetting of biodiversity, speciation, by considering explicitly some less convenient realities of natural populations? Specifically, how does genetic variation at incompatibility loci within a species influence interspecies hybridization upon secondary contact? And, in nature, how repeatable among distinct bouts of secondary contact are the genomic outcomes of hybridization? Mandeville et al. ( 2015 ) tackle exactly this question in their new study in Molecular Ecology on five species of suckers, fish of the genus Catostomus, that overlap sympatrically in different portions of their subdivided ranges that occupy different rivers. They document substantial genomic heterogeneity in realized hybridization in nature, both among species pairs and among the source populations for hybrids of a given species pair. This imperfect repeatability of episodes of hybridization implies greater permeability of species barriers in some parts of their range, with intriguing consequences for how the integrity of species as independently evolving units could be susceptible to collapse.     

Empirical demonstration of hybrid chromosomal races in house mice

Western house mice (Mus musculus domesticus) and common shrews (Sorex araneus) are important models for study of chromosomal speciation. Both had ancestral karyotypes consisting of telocentric chromosomes, and each is subdivided into numerous chromosomal races many of which have resulted from fixation of new mutations (Robertsonian fusions and whole‐arm reciprocal translocations). However, some chromosomal races in both species may alternatively have originated through hybridization, with particular homozygous recombinant products reaching fixation. Here, we demonstrate the process of generation of hybrid chromosomal races for the first time in either species using molecular markers. Analysis of centromeric microsatellite markers show that the Mid Valtellina (IMVA) and Upper Valtellina (IUVA) chromosomal races of the house mouse are recombinant products of hybridization of the Lower Valtellina (ILVA) and Poschiavo (CHPO) chromosomal races, supporting earlier theoretical analysis. IMVA and IUVA occupy a small area of the Italian Alps where ILVA makes contact with CHPO. IUVA and CHPO have previously been shown to be reproductively isolated in one village, emphasizing that hybrid chromosomal races in small mammals, as in plants, have the potential to be part of the speciation process.     

Parallel speciation or long‐distance dispersal? Lessons from seaweeds (Fucus) in the Baltic Sea

Parallel evolution has been invoked as a forceful mechanism of ecotype and species formation in many animal taxa. However, parallelism may be difficult to separate from recently monophyletically diverged species that are likely to show complex genetic relationships as a result of considerable shared ancestral variation and secondary hybridization in local areas. Thus, species' degrees of reproductive isolation, barriers to dispersal and, in particular, limited capacities for long‐distance dispersal will affect demographical structures underlying mechanisms of divergent evolution. Here, we used nine microsatellite DNA markers to study intra‐ and interspecific genetic diversity of two recently diverged species of brown macroalgae, Fucus radicans (L. Bergström & L. Kautsky) and F. vesiculosus (Linnaeus), in the Baltic Sea. We further performed biophysical modelling to identify likely connectivity patterns influencing the species' genetic structures. For each species, we found intraspecific contrasting patterns of clonality incidence and population structure. In addition, strong genetic differentiation between the two species within each locality supported the existence of two distinct evolutionary lineages (F_ST = 0.15–0.41). However, overall genetic clustering analyses across both species' populations revealed that all populations from one region (Estonia) were more genetically similar to each other than to their own taxon from the other two regions (Sweden and Finland). Our data support a hypothesis of parallel speciation. Alternatively, Estonia may be the ancestral source of both species, but is presently isolated by oceanographic barriers to dispersal. Thus, a limited gene flow in combination with genetic drift could have shaped the seemingly parallel structure.     

Contact zone dynamics during early stages of speciation in a chorus frog (Pseudacris crucifer)

Characterizing the genetic and behavioural consequences of contact between previously geographically isolated lineages provides insights into the mechanisms underlying diversification and ultimately speciation. The spring peeper (Pseudacris crucifer) is a widespread Nearctic chorus frog with six divergent mitochondrial DNA (mtDNA) lineages, many of which came into secondary contact during the Holocene. We examined genetics, morphology, advertisement calls and female preference for two lineages that began diverging in allopatry in the Pliocene and now overlap in southwestern Ontario, Canada. We found non-coincident clines in mtDNA and nuclear DNA, mirroring directionality of premating isolation barriers. We also found divergence in a range of traits between these two lineages, displacement in male call attributes and female preference for calls of their natal lineage in sympatry. Hybrids were morphologically distinct from both parental lineages, but hybrid male calls were acoustically intermediate. Female hybrids showed asymmetrical preference for Eastern male calls. These results considered together provide evidence of either unidirectional hybridization or selection against hybrids, potentially implying reproductive character displacement. Our work demonstrates the utility of integrated, multi-character approaches to understanding the processes of divergence and the nature of speciation.     

Structure and dynamics of hybrid zones at different stages of speciation in the common vole (Microtus arvalis)

The genetic structure and dynamics of hybrid zones provide crucial information for understanding the processes and mechanisms of evolutionary divergence and speciation. In general, higher levels of evolutionary divergence between taxa are more likely to be associated with reproductive isolation and may result in suppressed or strongly restricted hybridization. In this study, we examined two secondary contact zones between three deep evolutionary lineages in the common vole (Microtus arvalis). Differences in divergence times between the lineages can shed light on different stages of reproductive isolation and thus provide information on the ongoing speciation process in M. arvalis. We examined more than 800 individuals for mitochondrial (mtDNA), Y‐chromosome and autosomal markers and used assignment and cline analysis methods to characterize the extent and direction of gene flow in the contact zones. Introgression of both autosomal and mtDNA markers in a relatively broad area of admixture indicates selectively neutral hybridization between the least‐divergent lineages (Central and Eastern) without evidence for partial reproductive isolation. In contrast, a very narrow area of hybridization, shifts in marker clines and the quasi‐absence of Y‐chromosome introgression support a moving hybrid zone and unidirectional selection against male hybrids between the lineages with older divergence (Central and Western). Data from a replicate transect further support non‐neutral processes in this hybrid zone and also suggest a role for landscape history in the movement and shaping of geneflow profiles.     

Molecular investigation of genetic assimilation during the rapid adaptive radiations of East African cichlid fishes

Adaptive radiations are characterized by adaptive diversification intertwined with rapid speciation within a lineage resulting in many ecologically specialized, phenotypically diverse species. It has been proposed that adaptive radiations can originate from ancestral lineages with pronounced phenotypic plasticity in adaptive traits, facilitating ecologically driven phenotypic diversification that is ultimately fixed through genetic assimilation of gene regulatory regions. This study aimed to investigate how phenotypic plasticity is reflected in gene expression patterns in the trophic apparatus of several lineages of East African cichlid fishes, and whether the observed patterns support genetic assimilation. This investigation used a split brood experimental design to compare adaptive plasticity in species from within and outside of adaptive radiations. The plastic response was induced in the crushing pharyngeal jaws through feeding individuals either a hard or soft diet. We find that nonradiating, basal lineages show higher levels of adaptive morphological plasticity than the derived, radiated lineages, suggesting that these differences have become partially genetically fixed during the formation of the adaptive radiations. Two candidate genes that may have undergone genetic assimilation, gif and alas1, were identified, in addition to alterations in the wiring of LPJ patterning networks. Taken together, our results suggest that genetic assimilation may have dampened the inducibility of plasticity related genes during the adaptive radiations of East African cichlids, flattening the reaction norms and canalizing their feeding phenotypes, driving adaptation to progressively more narrow ecological niches.