Juvenile divergence in adaptive traits among seven sympatric fish ecomorphs arises before moving to different lacustrine habitats

Identifying the mechanisms initiating sympatric diversification in vertebrates has remained a conceptual challenge. Here, we analyse an assemblage of sympatric charr (Salvelinus malma) morphs from landlocked Lake Kronotskoe basin as a model to uncover the divergence pathways in freshwater fishes during the early life history stages. All morphs have distinct developmental biology, but a similar developmental rate retardation compared to the ancestor. Our study reveals that adult morphological differences, which acquire functionality at maturation, originate in the early juvenile stages due to heterochrony in skeletogenesis and allometric changes triggered by variation in metabolic activity. The craniofacial differences among the morphs result from asynchronous development of several skeletal modules. The accelerated ossification of teeth‐armed bones occurs in predatory feeding morphs, whereas cranial cover ossification is promoted in benthivorous morphs. These contrasting growth patterns have led to seven phenotypes that span a range far beyond the ancestral variability. The most distinct morphs are a riverine spawning, epilimnetic predator and a lacustrine spawning, profundal benthic feeder. Taken together, we argue that the adaptive morphological differentiation in these sympatric freshwater fishes is driven by diverging patterns in ossification rate and metabolic activity against a background of uneven somatic growth. This divergence is primarily associated with basic environmental differences on the nursery grounds that might be unrelated to resource use. This nonheritable phenotype divergence is then exposed to natural selection that could result in further adaptive genetic changes.     

Flickers of speciation: Sympatric colour morphs of the arc‐eye hawkfish,Paracirrhites arcatus,reveal key elements of divergence with gene flow

One of the primary challenges of evolutionary research is to identify ecological factors that favour reproductive isolation. Therefore, studying partially isolated taxa has the potential to provide novel insight into the mechanisms of evolutionary divergence. Our study utilizes an adaptive colour polymorphism in the arc‐eye hawkfish (Paracirrhites arcatus) to explore the evolution of reproductive barriers in the absence of geographic isolation. Dark and light morphs are ecologically partitioned into basaltic and coral microhabitats a few metres apart. To test whether ecological barriers have reduced gene flow among dark and light phenotypes, we evaluated genetic variation at 30 microsatellite loci and a nuclear exon (Mc1r) associated with melanistic coloration. We report low, but significant microsatellite differentiation among colour morphs and stronger divergence in the coding region of Mc1r indicating signatures of selection. Critically, we observed greater genetic divergence between colour morphs on the same reefs than that between the same morphs in different geographic locations. We hypothesize that adaptation to the contrasting microhabitats is overriding gene flow and is responsible for the partial reproductive isolation observed between sympatric colour morphs. Combined with complementary studies of hawkfish ecology and behaviour, these genetic results indicate an ecological barrier to gene flow initiated by habitat selection and enhanced by assortative mating. Hence, the arc‐eye hawkfish fulfil theoretical expectations for the earliest phase of speciation with gene flow.     

Parallels,nonparallels, and plasticity in population differentiation of threespine stickleback within a lake

The frequent occurrence of parallel phenotypic divergence in similar habitats is often evoked when emphasizing the role of ecology in adaptive radiation and speciation. However, because phenotypic plasticity can contribute to the observed pattern of divergence, confirmation of divergence at loci underlying phenotypic traits is important for confirming adaptive divergence. In the present study, we examine parallel morphological, neutral, and potentially adaptive genetic divergence of threespine stickleback inhabiting different habitats within a lake. Three genetic clusters best explained the neutral genetic structure within the lake; however, morphological differences were only weakly connected to genetic clusters and there was considerable phenotypic variation within clusters. Among the factors that could contribute to the observed pattern of morphological and genetic divergence are phenotypic plasticity, selective mortality of hybrids, and habitat choice based on morphology. Several loci are identified as outliers indicating divergent selection between the morphs and some parallels in morphological and adaptive genetic divergence are found in stickleback spawning at two lava sites. However, neutral genetic structure indicates considerable genetic connectivity among the two lava sites, and the parallels in morphology may therefore represent selective distribution of phenotypes rather than parallel divergence. © 2009 The Linnean Society of London, Biological Journal of the Linnean Society, 2009, 98 , 803–813.     

Genetic structure is correlated with phenotypic divergence rather than geographic isolation in the highly polymorphic strawberry poison‐dart frog

Phenotypic and genetic divergence can be influenced by a variety of factors, including sexual and natural selection, genetic drift and geographic isolation. Investigating the roles of these factors in natural systems can provide insight into the relative influences of allopatric and ecological modes of biological diversification in nature. The strawberry poison frog, Dendrobates pumilio, presents an excellent opportunity for this kind of research, displaying a diverse array of colour morphs and inhabiting a heterogeneous landscape that includes oceanic islands, fragmented rainforest patches and wide expanses of suitable habitat. In this study, we use 15 highly polymorphic microsatellite loci to estimate population structure and gene flow among populations from across the range of D. pumilio and a causal modelling framework to statistically test 12 hypotheses regarding the geographic and phenotypic variables that explain genetic differentiation within this system. Our results demonstrate that the genetic distance between populations is most strongly associated with differences in dorsal coloration. Previous experimental studies have shown that phenotypic differences can result in sexual and natural selection against non‐native phenotypes, and our results now show that these forces lead to genetic isolation between different colour morphs in the wild, presenting a potential case of incipient speciation through selection.     

Recent speciation between sympatric Tanganyikan cichlid colour morphs

Lake Tanganyika, Africa's oldest lake, harbours an impressive diversity of cichlid fishes. Although diversification in its radiating groups is thought to have been initially rapid, cichlids from Lake Tanganyika show little evidence for ongoing speciation. In contrast, examples of recent divergence among sympatric colour morphs are well known in haplochromine cichlids from Lakes Malawi and Victoria. Here, we report genetic evidence for recent divergence between two sympatric Tanganyikan cichlid colour morphs. These Petrochromis morphs share mitochondrial haplotypes, yet microsatellite loci reveal that their sympatric populations form distinct genetic groups. Nuclear divergence between the two morphs is equivalent to that which arises geographically within one of the morphs over short distances and is substantially smaller than that among other sympatric species in this genus. These patterns suggest that these morphs diverged only recently, yet that barriers to gene flow exist which prevent extensive admixture despite their sympatric distribution. The morphs studied here provide an unusual example of active diversification in Lake Tanganyika's generally ancient cichlid fauna and enable comparisons of speciation processes between Lake Tanganyika and other African lakes.     

The interaction of resource use and gene flow on the phenotypic divergence of benthic and pelagic morphs of Icelandic Arctic charr (Salvelinus alpinus)

Conceptual models of adaptive divergence and ecological speciation in sympatry predict differential resource use, phenotype–environment correlations, and reduced gene flow among diverging phenotypes. While these predictions have been assessed in past studies, connections among them have rarely been assessed collectively. We examined relationships among phenotypic, ecological, and genetic variation in Arctic charr (Salvelinus alpinus) from six Icelandic localities that have undergone varying degrees of divergence into sympatric benthic and pelagic morphs. We characterized morphological variation with geometric morphometrics, tested for differential resource use between morphs using stable isotopes, and inferred the amount of gene flow from single nucleotide polymorphisms. Analysis of stable isotopic signatures indicated that sympatric morphs showed similar difference in resource use across populations, likely arising from the common utilization of niche space within each population. Carbon isotopic signature was also a significant predictor of individual variation in body shape and size, suggesting that variation in benthic and pelagic resource use is associated with phenotypic variation. The estimated percentage of hybrids between sympatric morphs varied across populations (from 0% to 15.6%) but the majority of fish had genotypes (ancestry coefficients) characteristic of pure morphs. Despite evidence of reduced gene flow between sympatric morphs, we did not detect the expected negative relationship between divergence in resource use and gene flow. Three lakes showed the expected pattern, but morphs in the fourth showed no detectable hybridization and had relatively low differences in resource use between them. This coupled with the finding that resource use and genetic differentiation had differential effects on body shape variation across populations suggests that reproductive isolation maintains phenotypic divergence between benthic and pelagic morphs when the effects of resource use are relatively low. Our ability to assess relationships between phenotype, ecology, and genetics deepens our understanding of the processes underlying adaptive divergence in sympatry.     

Convergent evolutionary processes driven by foraging opportunity in two sympatric morph pairs of Arctic charr with contrasting post‐glacial origins

The expression of two or more discrete phenotypes amongst individuals within a species (morphs) provides multiple modes upon which selection can act semi‐independently, and thus may be an important stage in speciation. In the present study, we compared two sympatric morph systems aiming to address hypotheses related to their evolutionary origin. Arctic charr in sympatry in Loch Tay, Scotland, exhibit one of two discrete, alternative body size phenotypes at maturity (large or small body size). Arctic charr in Loch Awe segregate into two temporally segregated spawning groups (breeding in either spring or autumn). Mitochondrial DNA restriction fragment length polymorphism analysis showed that the morph pairs in both lakes comprise separate gene pools, although segregation of the Loch Awe morphs is more subtle than that of Loch Tay. We conclude that the Loch Awe morphs diverged in situ (within the lake), whereas Loch Tay morphs most likely arose through multiple invasions by different ancestral groups that segregated before post‐glacial invasion (i.e. in allopatry). Both morph pairs showed clear trophic segregation between planktonic and benthic resources (measured by stable isotope analysis) but this was significantly less distinct in Loch Tay than in Loch Awe. By contrast, both inter‐morph morphological and life‐history differences were more subtle in Loch Awe than in Loch Tay. The strong ecological but relatively weak morphological and life‐history divergence of the in situ derived morphs compared to morphs with allopatric origins indicates a strong link between early ecological and subsequent genetic divergence of sympatric origin emerging species pairs. The emergence of parallel specialisms despite distinct genetic origins of these morph pairs suggests that the effect of available foraging opportunities may be at least as important as genetic origin in structuring sympatric divergence in post‐glacial fishes with high levels of phenotypic plasticity. © 2012 The Linnean Society of London, Biological Journal of the Linnean Society, 2012, ?? , ??–??.     

Effects of population,family, and diet on craniofacial morphology of Icelandic Arctic charr (Salvelinus alpinus)

We evaluated hypotheses of intralacustrine diversification and plastic responses to two diet environments in Icelandic Arctic charr (Salvelinus alpinus). Full‐sib families of progeny of wild polymorphic charr from two lakes where morphs vary in their degree of phenotypic and ecological divergence were split, with half of the offspring reared on a benthic and half on a limnetic type of diet to estimate family norms of reaction. We focused on variation in craniofacial traits because they are probably functionally related to diet and complement a previous study of body shape in these charr. A hierarchical analysis of phenotypic variation between lakes, pairs of morphs within each lake, and two families within each morph found that phenotypic variation partitioned between families relative to morphs was reduced in the more ecologically diversified population, which is consistent with adaptive diversification. The effect size of plastic responses between lake populations was similar, suggesting little difference in the degree of canalization in contrast to a previous analysis of body form plasticity. Thus, the role that plastic morphological responses play in the adaptive diversification of morphs and different lake populations of Arctic charr may depend on the trait. © 2013 The Linnean Society of London     

DETERMINISM IN THE DIVERSIFICATION OF HISPANIOLAN TRUNK‐GROUND ANOLES (ANOLIS CYBOTES SPECIES COMPLEX)

The evolutionary processes that produce adaptive radiations are enigmatic. They can only be studied after the fact, once a radiation has occurred and been recognized, rather than while the processes are ongoing. One way to connect pattern to process is to study the processes driving divergence today among populations of species that belong to an adaptive radiation, and compare the results to patterns observed at a deeper, macroevolutionary level. We tested whether evolution is a deterministic process with similar outcomes during different stages of the adaptive radiation of Anolis lizards. Using a clade of terrestrial–scansorial lizards in the genus Anolis, we inferred the adaptive basis of spatial variation among contemporary populations and tested whether axes of phenotypic differentiation among them mirror known axes of diversification at deeper levels of the anole radiation. Nonparallel change associated with genetic divergence explains the vast majority of geographic variation. However, we found phenotypic variation to be adaptive as confirmed by convergence in populations occurring in similar habitats in different mountain ranges. Morphological diversification among populations recurs deterministically along two axes of diversification previously identified in the anole radiation, but the characters involved differ from those involved in adaptation at higher levels of anole phylogeny.     

Genomic regions repeatedly involved in divergence among plant‐specialized pea aphid biotypes

Understanding the genetic bases of biological diversification is a long‐standing goal in evolutionary biology. Here, we investigate whether replicated cases of adaptive divergence involve the same genomic regions in the pea aphid, Acyrthosiphon pisum, a large complex of genetically differentiated biotypes, each specialized on different species of legumes. A previous study identified genomic regions putatively involved in host‐plant adaptation and/or reproductive isolation by performing a hierarchical genome scan in three biotypes. This led to the identification of 11 F_ST outliers among 390 polymorphic microsatellite markers. In this study, the outlier status of these 11 loci was assessed in eight biotypes specialized on other host plants. Four of the 11 previously identified outliers showed greater genetic differentiation among these additional biotypes than expected under the null hypothesis of neutral evolution (α < 0.01). Whether these hotspots of genomic divergence result from adaptive events, intrinsic barriers or reduced recombination is discussed.     

Rampant host- and defensive phenotype-associated diversification in a goldenrod gall midge

Natural selection can play an important role in the genetic divergence of populations and their subsequent speciation. Such adaptive diversification, or ecological speciation, might underlie the enormous diversity of plant-feeding insects that frequently experience strong selection pressures associated with host plant use as well as from natural enemies. This view is supported by increasing documentation of host-associated (genetic) differentiation in populations of plant-feeding insects using alternate hosts. Here, we examine evolutionary diversification in a single nominal taxon, the gall midge Asteromyia carbonifera (O.S.), with respect to host plant use and gall phenotype. Because galls can be viewed as extended defensive phenotypes of the midges, gall morphology is likely to be a reflection of selective pressures by enemies. Using phylogenetic and comparative analyses of mtDNA and nuclear sequence data, we find evidence that A. carbonifera populations are rapidly diversifying along host plant and gall morphological lines. At a broad scale, geography explains surprisingly little genetic variation, and there is little evidence of strict co-cladogenesis with their Solidago hosts. Gall morphology is relatively labile, distinct gall morphs have evolved repeatedly and colonized multiple hosts, and multiple genetically and morphologically distinct morphs frequently coexist on a single host plant species. These results suggest that Asteromyia carbonifera is in the midst of an adaptive radiation driven by multitrophic selective pressures. Similar complex community pressures are likely to play a role in the diversification of other herbivorous insect groups.     

Evidence for divergence and adaptive isolation in post‐glacially derived bimodal allopatric and sympatric rainbow smelt populations

The roles of ecological speciation and reinforcement in the formation of contemporary diversity remain contentious. In the present study, we contrast phenotypic and molecular divergence within morphologically diverged bimodal sympatric and allopatric pairs of rainbow smelt, Osmerus mordax. We hypothesize that, in sympatry, evidence of selection associated with resource partitioning will be visible through strong divergence, reinforcement, and greater character displacement. Parallel morphological divergence was observed between the two trophic forms (macrophagous and microphagous), with several examples of greater trait divergence in sympatry than allopatry. Mitochondrial DNA sequence analysis indicated no association between historical clades and morphology; however, Bayesian clustering using microsatellites supported the isolation of these morphs under both allopatry and sympatry. Estimates of genetic isolation were one order of magnitude lower than measures of morphological divergence, consistent with a hypothesis of strong contemporary selection. Using experimental crosses, we obtained similar rates of fertilization success among the allopatric hybrid and pure crosses; whereas, in the sympatric hybrid crosses, fertilization rates dropped by 30–50%, suggesting a clear role for reinforcement through prezygotic incompatibilities. The present study supports the hypothesis that processes of post‐glacial radiation and diversification differ between sympatry and allopatry, and indicates a role for reinforcement and ecological processes in recent sympatric diversification. © 2010 The Linnean Society of London, Biological Journal of the Linnean Society, 2010, 101 , 583–594.     

Evolutionary radiation of Asteromyia carbonifera (Diptera: Cecidomyiidae) gall morphotypes on the goldenrod Solidago altissima (Asteraceae)

Population divergence of phytophagous insects is often coupled to host‐plant shifts and is frequently attributed to the divergent selective environments associated with alternative host‐plants. In some cases, however, divergence is associated with the use of alternative host‐plant organs of a single host species. The basis of within‐host radiations such as these remains poorly understood. In the present stusy, we analysed the radiation of Asteromyia gall midges occurring both within one host plant species and within a single organ on that host. In this system, four morphologically distinct Asteromyia gall forms (morphs) coexist on the leaves of goldenrod Solidago altissima. Our analyses of amplified fragment length polymorphism and DNA sequence data confirm the genetic differentiation among midges from three gall morphs and reveal evidence of a genetically distinct fourth gall morph. The absence of clear gall morph related clades in the mitochondrial DNA derived phylogenies is indicative of incomplete lineage sorting or recent gene flow, suggesting that population divergence among gall forms is recent. We assess the likely history of this radiation and use the results of phylogenetic analyses along with ecological data on phenology and parasitism rates to evaluate potential hypotheses for the mode of differentiation. These preliminary analyses suggest that diversification of the Asteromyia gall morphs is likely shaped by interactions between the midge, a symbiotic fungus, and parasitoid enemies. © 2008 The Linnean Society of London, Biological Journal of the Linnean Society, 2008, 95 , 840–858.     

Gudgeon fish with and without genetically determined countershading coexist in heterogeneous littoral environments of an ancient lake

Countershading, characterized by a darker dorsal surface and lighter ventral surface, is common among many animals. This dorsoventral pigment polarity is often thought to be adaptive coloration for camouflage. By contrast, noncountershaded (melanistic) morphs often occur within a species due to genetic color polymorphism in terrestrial animals. However, the polymorphism with either countershaded or melanistic morphs is poorly known in wild aquatic animals. This study explored the genetic nature of diverged color morphs of a lineage of gudgeon fish (genus Sarcocheilichthys) in the ancient Lake Biwa and propose this system as a novel model for testing hypotheses of functional aspects of countershading and its loss in aquatic environments. This system harbors two color morphs that have been treated taxonomically as separate species; Sarcocheilichthys variegatus microoculus which occurs throughout the littoral zone and Sarcocheilichthys biwaensis which occurs in and around rocky areas. First, we confirmed that the divergence of dorsoventral color patterns between the two morphs is under strict genetic control at the levels of chromatophore distribution and melanin‐related gene expression under common garden rearing. The former morph displayed sharp countershading coloration, whereas the latter morph exhibited a strong tendency toward its loss. The crossing results indicated that this divergence was likely controlled by a single locus in a two‐allele Mendelian inheritance pattern. Furthermore, our population genomic and genome‐wide association study analyses detected no genome‐wide divergence between the two morphs, except for one region near a locus that may be associated with the color divergence. Thus, these morphs are either in a state of intraspecific color polymorphism or two incipient species. Evolutionary forces underlying this polymorphism appear to be associated with heterogeneous littoral environments in this lake. Future ecological genomic research will provide insight into adaptive functions of this widespread coloration, including the eco‐evolutionary drivers of its loss, in the aquatic world.     

Island biology and morphological divergence of the Skyros wall lizard <Emphasis Type="Italic">Podarcis gaigeae</Emphasis>: a combined role for local selection and genetic drift on color morph frequency divergence?

Background  

Patterns of spatial variation in discrete phenotypic traits can be used to draw inferences about the adaptive significance of traits and evolutionary processes, especially when compared to patterns of neutral genetic variation. Population divergence in adaptive traits such as color morphs can be influenced by both local ecology and stochastic factors such as genetic drift or founder events. Here, we use quantitative color measurements of males and females of Skyros wall lizard, Podarcis gaigeae, to demonstrate that this species is polymorphic with respect to throat color, and the morphs form discrete phenotypic clusters with limited overlap between categories. We use divergence in throat color morph frequencies and compare that to neutral genetic variation to infer the evolutionary processes acting on islet- and mainland populations.     

Machine‐learning‐based detection of adaptive divergence of the stream mayfly Ephemera strigata populations

Adaptive divergence is a key mechanism shaping the genetic variation of natural populations. A central question linking ecology with evolutionary biology is how spatial environmental heterogeneity can lead to adaptive divergence among local populations within a species. In this study, using a genome scan approach to detect candidate loci under selection, we examined adaptive divergence of the stream mayfly Ephemera strigata in the Natori River Basin in northeastern Japan. We applied a new machine‐learning method (i.e., random forest) besides traditional distance‐based redundancy analysis (dbRDA) to examine relationships between environmental factors and adaptive divergence at non‐neutral loci. Spatial autocorrelation analysis based on neutral loci was employed to examine the dispersal ability of this species. We conclude the following: (a) E. strigata show altitudinal adaptive divergence among the populations in the Natori River Basin; (b) random forest showed higher resolution for detecting adaptive divergence than traditional statistical analysis; and (c) separating all markers into neutral and non‐neutral loci could provide full insight into parameters such as genetic diversity, local adaptation, and dispersal ability.     

Parallel and non‐parallel morphological divergence among foraging specialists in European whitefish (Coregonus lavaretus)

Parallel phenotypic evolution occurs when independent populations evolve similar traits in response to similar selective regimes. However, populations inhabiting similar environments also frequently show some phenotypic differences that result from non‐parallel evolution. In this study, we quantified the relative importance of parallel evolution to similar foraging regimes and non‐parallel lake‐specific effects on morphological variation in European whitefish (Coregonus lavaretus). We found evidence for both lake‐specific morphological characteristics and parallel morphological divergence between whitefish specializing in feeding on profundal and littoral resources in three separate lakes. Foraging specialists expressed similar phenotypes in different lakes in both overall body shape and selected measured morphological traits. The morphology of the two whitefish specialists resembled that predicted from other fish species, supporting the conclusion of an adaptive significance of the observed morphological characteristics. Our results indicate that divergent natural selection resulting from foraging specialization is driving and/or maintaining the observed parallel morphological divergence. Whitefish in this study may represent an early stage of divergence towards the evolution of specialized morphs.     

The role of structural genomic variants in population differentiation and ecotype formation in Timema cristinae walking sticks

Theory predicts that structural genomic variants such as inversions can promote adaptive diversification and speciation. Despite increasing empirical evidence that adaptive divergence can be triggered by one or a few large inversions, the degree to which widespread genomic regions under divergent selection are associated with structural variants remains unclear. Here we test for an association between structural variants and genomic regions that underlie parallel host‐plant‐associated ecotype formation in Timema cristinae stick insects. Using mate‐pair resequencing of 20 new whole genomes we find that moderately sized structural variants such as inversions, deletions and duplications are widespread across the genome, being retained as standing variation within and among populations. Using 160 previously published, standard‐orientation whole genome sequences we find little to no evidence that the DNA sequences within inversions exhibit accentuated differentiation between ecotypes. In contrast, a formerly described large region of reduced recombination that harbours genes controlling colour‐pattern exhibits evidence for accentuated differentiation between ecotypes, which is consistent with differences in the frequency of colour‐pattern morphs between host‐associated ecotypes. Our results suggest that some types of structural variants (e.g., large inversions) are more likely to underlie adaptive divergence than others, and that structural variants are not required for subtle yet genome‐wide genetic differentiation with gene flow.