Genomic architecture of ecologically divergent body shape in a pair of sympatric crater lake cichlid fishes

Determining the genetic bases of adaptations and their roles in speciation is a prominent issue in evolutionary biology. Cichlid fish species flocks are a prime example of recent rapid radiations, often associated with adaptive phenotypic divergence from a common ancestor within a short period of time. In several radiations of freshwater fishes, divergence in ecomorphological traits — including body shape, colour, lips and jaws — is thought to underlie their ecological differentiation, specialization and, ultimately, speciation. The Midas cichlid species complex (Amphilophus spp.) of Nicaragua provides one of the few known examples of sympatric speciation where species have rapidly evolved different but parallel morphologies in young crater lakes. This study identified significant QTL for body shape using SNPs generated via ddRAD sequencing and geometric morphometric analyses of a cross between two ecologically and morphologically divergent, sympatric cichlid species endemic to crater Lake Apoyo: an elongated limnetic species (Amphilophus zaliosus) and a high‐bodied benthic species (Amphilophus astorquii). A total of 453 genome‐wide informative SNPs were identified in 240 F₂ hybrids. These markers were used to construct a genetic map in which 25 linkage groups were resolved. Seventy‐two segregating SNPs were linked to 11 QTL. By annotating the two most highly supported QTL‐linked genomic regions, genes that might contribute to divergence in body shape along the benthic–limnetic axis in Midas cichlid sympatric adaptive radiations were identified. These results suggest that few genomic regions of large effect contribute to early stage divergence in Midas cichlids.     

Evidence for ecological speciation via a host shift in the holly leaf miner,Phytomyza glabricola (Diptera: Agromyzidae)

Evolutionary radiations have been well documented in plants and insects, and natural selection may often underly these radiations. If radiations are adaptive, the diversity of species could be due to ecological speciation in these lineages. Agromyzid flies exhibit patterns of repeated host‐associated radiations. We investigated whether host‐associated population divergence and evidence of divergent selection exist in the leaf miner Phytomyza glabricola on its sympatric host plants, the holly species, Ilex coriacea and I. glabra. Using AFLPs and nuclear sequence data, we found substantial genetic divergence between host‐associated populations of these flies throughout their geographic range. Genome scans using the AFLP data identified 13 loci under divergent selection, consistent with processes of ecological speciation. EF‐1α data suggest that I. glabra is the original host of P. glabricola and that I. coriacea is the novel host, but the AFLP data are ambiguous with regard to directionality of the host shift.     

Weak disruptive selection and incomplete phenotypic divergence in two classic examples of sympatric speciation: cameroon crater lake cichlids

Abstract Recent documentation of a few compelling examples of sympatric speciation led to a proliferation of theoretical models. Unfortunately, plausible examples from nature have rarely been used to test model predictions, such as the initial presence of strong disruptive selection. Here I estimated the form and strength of selection in two classic examples of sympatric speciation: radiations of Cameroon cichlids restricted to Lakes Barombi Mbo and Ejagham. I measured five functional traits and relative growth rates in over 500 individuals within incipient species complexes from each lake. Disruptive selection was prevalent in both groups on single and multivariate trait axes but weak relative to stabilizing selection on other traits and most published estimates of disruptive selection. Furthermore, despite genetic structure, assortative mating, and bimodal species-diagnostic coloration, trait distributions were unimodal in both species complexes, indicating the earliest stages of speciation. Long waiting times or incomplete sympatric speciation may result when disruptive selection is initially weak. Alternatively, I present evidence of additional constraints in both species complexes, including weak linkage between coloration and morphology, reduced morphological variance aligned with nonlinear selection surfaces, and minimal ecological divergence. While other species within these radiations show complete phenotypic separation, morphological and ecological divergence in these species complexes may be slow or incomplete outside optimal parameter ranges, in contrast to rapid divergence of their sexual coloration.     

Divergence in coloration and ecological speciation in the Anolis marmoratus species complex

Adaptive divergence in coloration is expected to produce reproductive isolation in species that use colourful signals in mate choice and species recognition. Indeed, many adaptive radiations are characterized by differentiation in colourful signals, suggesting that divergent selection acting on coloration may be an important component of speciation. Populations in the Anolis marmoratus species complex from the Caribbean island of Guadeloupe display striking divergence in the colour and pattern of adult males that occurs over small geographic distances, suggesting strong divergent selection. Here we test the hypothesis that divergence in coloration results in reduced gene flow among populations. We quantify variation in adult male coloration across a habitat gradient between mesic and xeric habitats, use a multilocus coalescent approach to infer historical demographic parameters of divergence, and examine gene flow and population structure using microsatellite variation. We find that colour variation evolved without geographic isolation and in the face of gene flow, consistent with strong divergent selection and that both ecological and sexual selection are implicated. However, we find no significant differentiation at microsatellite loci across populations, suggesting little reproductive isolation and high levels of contemporary gene exchange. Strong divergent selection on loci affecting coloration probably maintains clinal phenotypic variation despite high gene flow at neutral loci, supporting the notion of a porous genome in which adaptive portions of the genome remain fixed whereas neutral portions are homogenized by gene flow and recombination. We discuss the impact of these findings for studies of colour evolution and ecological speciation.     

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.     

Genomic signatures of divergent selection and speciation patterns in a ‘natural experiment’, the young parallel radiations of Nicaraguan crater lake cichlid fishes

Divergent selection is the main driving force in sympatric ecological speciation and may also play a strong role in divergence between allopatric populations. Characterizing the genome‐wide impact of divergent selection often constitutes a first step in unravelling the genetic bases underlying adaptation and ecological speciation. The Midas cichlid fish (Amphilophus citrinellus) species complex in Nicaragua is a powerful system for studying evolutionary processes. Independent colonizations of isolated young crater lakes by Midas cichlid populations from the older and great lakes of Nicaragua resulted in the repeated evolution of adaptive radiations by intralacustrine sympatric speciation. In this study we performed genome scans on two repeated radiations of crater lake species and their great lake source populations (1030 polymorphic AFLPs, n ～ 30 individuals per species). We detected regions under divergent selection (0.3% in the crater lake Xiloá flock and 1.7% in the older crater lake Apoyo radiation) that might be responsible for the sympatric diversifications. We find no evidence that the same genomic regions have been involved in the repeated evolution of parallel adaptations across crater lake flocks. However, there is some genetic parallelism apparent (seven out of 51 crater lake to great lake outlier loci are shared; 13.7%) that is associated with the allopatric divergence of both crater lake flocks. Interestingly, our results suggest that the number of outlier loci involved in sympatric and allopatric divergence increases over time. A phylogeny based on the AFLP data clearly supports the monophyly of both crater lake species flocks and indicates a parallel branching order with a primary split along the limnetic‐benthic axis in both radiations.     

Sexual selection,reproductive behavior,and speciation in the mbuna species flock of Lake Malawi (Pisces: Cichlidae)

Synopsis The cichlid species flocks of the African Great Lakes represent the most extreme case of adaptive radiation among vertebrates. Recently, attention has focused on the potential for sexual selection to drive or accelerate speciation in these fishes. Cichlids as a whole are social in nature and display complex behavior, particularly during courtship and spawning; however, the extent to which changes in species recognition cues may account for species diversity among haplochromine lineages has remained speculative. Our investigations have indicated that oral incubating haplochromines show a reduction in diversity and extent of courtship relative to substrate brooding cichlids, and apparently retain aspects of a primitive specific mate recognition system. Laboratory observations of courtship in the Malawian endemicPseudotropheus zebra suggest that organization of the spawning bout is loose, and lacking in any well defined stimulus-response chain. Interspecific comparisons of in situ courtship behavior among male mbuna, lithophilous haplochromines of Lake Malawi, revealed only one potential example of species-specific behavior, and indicated that mate choice occurred prior to the onset of intense courtship. Courtship display by male mbuna does not appear critical to species recognition and may represent an evolutionary relict. Alternatively, male courtship display may have an indirect effect on successful reproduction or may be maintained through pleiotropic effects. The mbuna display no evidence of behavioral innovation and show limited interspecific differentiation in behavioral expression. More generally, there is no evidence at present to suggest that epigamic sexual selection, acting on courtship behavior, has been a major mechanism in the diversification of the haplochromine species flocks.     

Disruptive natural selection predicts divergence between the sexes during adaptive radiation

Evolution of sexual dimorphism in ecologically relevant traits, for example, via resource competition between the sexes, is traditionally envisioned to stall the progress of adaptive radiation. An alternative view is that evolution of ecological sexual dimorphism could in fact play an important positive role by facilitating sex‐specific adaptation. How competition‐driven disruptive selection, ecological sexual dimorphism, and speciation interact during real adaptive radiations is thus a critical and open empirical question. Here, we examine the relationships between these three processes in a clade of salamanders that has recently radiated into divergent niches associated with an aquatic life cycle. We find that morphological divergence between the sexes has occurred in a combination of head shape traits that are under disruptive natural selection within breeding ponds, while divergence among species means has occurred independently of this disruptive selection. Further, we find that adaptation to aquatic life is associated with increased sexual dimorphism across taxa, consistent with the hypothesis of clade‐wide character displacement between the sexes. Our results suggest the evolution of ecological sexual dimorphism may play a key role in niche divergence among nascent species and demonstrate that ecological sexual dimorphism and ecological speciation can and do evolve concurrently in the early stages of adaptive radiation.     

Adaptive divergence and the evolution of reproductive isolation in the wild: an empirical demonstration using introduced sockeye salmon

Populations exposed to different ecological environments should diverge for phenotypic traits that influence survival and reproduction. This adaptive divergence should reduce gene flow between populations because immigrants become less fit than residents and because hybrids perform poorly in either environment (i.e., ecologically-dependent reproductive isolation). Here I demonstrate adaptive divergence and the evolution of reproductive isolation in populations of sockeye salmon (Oncorhynchus nerka) introduced from a common ancestral source into a new lake system (Lake Washington, Washington). The introduced fish founded several new populations, two of which experience very different environments during breeding and early development (Cedar River v.s. Pleasure Point beach). Over 13 generations, the two populations diverged for adult traits (female body size, male body depth; measured in the wild) and embryo traits (survival to hatching, development rate, size at emergence; measured in a common environment). The rates of divergence for these characters were similar to those observed in other examples of rapid evolution, and can best be attributed to natural selection. Partial reproductive isolation has evolved in concert with adaptive divergence: the rate of exchange of adults between the populations (determined using natural tags) is higher than the rate of gene flow (determined using DNA microsatellites). The demonstration that adaptive divergence can initiate reproductive isolation in less than 13 generations suggests that the first signs of ecological speciation may appear soon after new environments are first colonized.     

Sexual isolation promotes divergence between parapatric lake and stream stickleback

Speciation can be initiated by adaptive divergence between populations in ecologically different habitats, but how sexually based reproductive barriers contribute to this process is less well understood. We here test for sexual isolation between ecotypes of threespine stickleback fish residing in adjacent lake and stream habitats in the Lake Constance basin, Central Europe. Mating trials exposing females to pairings of territorial lake and stream males in outdoor mesocosms allowing for natural reproductive behaviour reveal that mating occurs preferentially between partners of the same ecotype. Compared to random mating, this sexual barrier reduces gene flow between the ecotypes by some 36%. This relatively modest strength of sexual isolation is surprising because comparing the males between the two ecotypes shows striking differentiation in traits generally considered relevant to reproductive behaviour (body size, breeding coloration, nest size). Analysing size differences among the individuals in the mating trials further indicates that assortative mating is not related to ecotype differences in body size. Overall, we demonstrate that sexually based reproductive isolation promotes divergence in lake–stream stickleback along with other known reproductive barriers, but we also caution against inferring strong sexual isolation from the observation of strong population divergence in sexually relevant traits.     

Phylogeographic and phenotypic outcomes of brown anole colonization across the Caribbean provide insight into the beginning stages of an adaptive radiation

Some of the most important insights into the ecological and evolutionary processes of diversification and speciation have come from studies of island adaptive radiations, yet relatively little research has examined how these radiations initiate. We suggest that Anolis sagrei is a candidate for understanding the origins of the Caribbean Anolis adaptive radiation and how a colonizing anole species begins to undergo allopatric diversification, phenotypic divergence and, potentially, speciation. We undertook a genomic and morphological analysis of representative populations across the entire native range of A. sagrei, finding that the species originated in the early Pliocene, with the deepest divergence occurring between western and eastern Cuba. Lineages from these two regions subsequently colonized the northern Caribbean. We find that at the broadest scale, populations colonizing areas with fewer closely related competitors tend to evolve larger body size and more lamellae on their toepads. This trend follows expectations for post‐colonization divergence from progenitors and convergence in allopatry, whereby populations freed from competition with close relatives evolve towards common morphological and ecological optima. Taken together, our results show a complex history of ancient and recent Cuban diaspora with populations on competitor‐poor islands evolving away from their ancestral Cuban populations regardless of their phylogenetic relationships, thus providing insight into the original diversification of colonist anoles at the beginning of the radiation. Our research also supplies an evolutionary framework for the many studies of this increasingly important species in ecological and evolutionary research.     

Divergence in host use ability of a marine herbivore from two habitat types

We studied ecological divergence of host use ability in a generalist marine herbivore living in two distinct host plant assemblages. We collected Idotea balthica isopods from three populations dominated by the brown alga Fucus vesiculosus and three dominated by the seagrass Zostera marina. In two reciprocal common garden feeding experiments for adult and laboratory‐born juvenile isopods, we found that isopods from both assemblages performed better with their sympatric dominant host species than did isopods allopatric to this host. This indicates parallel divergence of populations according to the sympatric host plant assemblage. Furthermore, initial body size and body size‐dependent mortality differed between populations from the two assemblages. In nature, this may result in lower fitness of immigrants compared with that of residents and consequently reinforce divergence of the populations. Finally, we discuss how phenotypic plasticity and maternal and random effects may associate with the results.     

Does sympatry predict life history and morphological diversification in the Mexican livebearing fish Poeciliopsis baenschi?

Understanding why some species coexist and others do not remains one of the fundamental challenges of ecology. Although there is evidence to suggest that closely‐related species are unlikely to occupy the same habitat because of competitive exclusion, there are many cases where closely‐related species do co‐occur. Research comparing sympatric and allopatric populations of co‐occurring species provides a framework for understanding the role of phenotypic diversification in species coexistence. In the present study, we compare phenotypic divergence between sympatric and allopatric populations of the livebearing fish, Poeciliopsis baenschi. We focus on life‐history traits and body shape, comprising two sets of integrated traits likely to diverge in response to varying selective pressures. Given that males and females can express different phenotypic traits, we also test for patterns of divergence among sexes by comparing size at maturity and sexual dimorphism in body shape between males and females in each population type. We take advantage of a natural experiment in western Mexico where, in some locations, P. baenschi co‐occur with a closely‐related species, Poeciliopsis turneri (sympatric populations) and, in other locations, they occur in isolation (allopatric populations). The results obtained in the present study show that sympatric populations of P. baenschi differed significantly in life‐history traits and in body shape compared to their allopatric counterparts. Additionally, males and females showed different responses for size at maturity in sympatric conditions versus allopatric conditions. However, the amount of sexual dimorphism did not differ between sympatric and allopatric populations of P. baenschi. Hence, we conclude that not all traits show similar levels of phenotypic divergence in response to sympatric conditions. © 2010 The Linnean Society of London, Biological Journal of the Linnean Society, 2010, 100 , 608–618.     

Change in sexual signalling traits outruns morphological divergence across an ecological gradient in the post‐glacial radiation of the songbird genus Junco

The relative roles of natural and sexual selection in promoting evolutionary lineage divergence remains controversial and difficult to assess in natural systems. Local adaptation through natural selection is known to play a central role in promoting evolutionary divergence, yet secondary sexual traits can vary widely among species in recent radiations, suggesting that sexual selection may also be important in the early stages of speciation. Here, we compare rates of divergence in ecologically relevant traits (morphology) and sexually selected signalling traits (coloration) relative to neutral structure in genome‐wide molecular markers and examine patterns of variation in sexual dichromatism to explore the roles of natural and sexual selection in the diversification of the songbird genus Junco (Aves: Passerellidae). Juncos include divergent lineages in Central America and several dark‐eyed junco (J. hyemalis) lineages that diversified recently as the group recolonized North America following the last glacial maximum (ca. 18,000 years ago). We found an accelerated rate of divergence in sexually selected characters relative to ecologically relevant traits. Moreover, sexual dichromatism measurements suggested a positive relationship between the degree of colour divergence and the strength of sexual selection when controlling for neutral genetic distance. We also found a positive correlation between dichromatism and latitude, which coincides with the geographic axis of decreasing lineage age in juncos but also with a steep ecological gradient. Finally, we found significant associations between genome‐wide variants linked to functional genes and proxies of both sexual and natural selection. These results suggest that the joint effects of sexual and ecological selection have played a prominent role in the junco radiation.     

STRONG ASSORTATIVE MATING BY DIET,COLOR, SIZE,AND MORPHOLOGY BUT LIMITED PROGRESS TOWARD SYMPATRIC SPECIATION IN A CLASSIC EXAMPLE: CAMEROON CRATER LAKE CICHLIDS

Models predict that sympatric speciation depends on restrictive parameter ranges, such as sufficiently strong disruptive selection and assortative mating, but compelling examples in nature have rarely been used to test these predictions. I measured the strength of assortative mating within a species complex of Tilapia in Lake Ejagham, Cameroon, a celebrated example of incipient sympatric adaptive radiation. This species complex is in the earliest stages of speciation: morphological and ecological divergence are incomplete, species differ primarily in breeding coloration, and introgression is common. I captured 27 mated pairs in situ and measured the diet, color, size, and morphology of each individual. I found strong assortative mating by color, size, head depth, and dietary source of benthic or pelagic prey along two independent dimensions of assortment. Thus, Ejagham Tilapia showed strong assortative mating most conducive to sympatric speciation. Nonetheless, in contrast to a morphologically bimodal Sarotherodon cichlid species pair in the lake, Ejagham Tilapia show more limited progress toward speciation, likely due to insufficient strength of disruptive selection on morphology estimated in a previous study (γ = 0.16). This supports the predicted dependence of sympatric speciation on strong assortment and strong disruptive selection by examining a potentially stalled example in nature.     

Causes of male sexual trait divergence in introduced populations of guppies

Males from different populations of the same species often differ in their sexually selected traits. Variation in sexually selected traits can be attributed to sexual selection if phenotypic divergence matches the direction of sexual selection gradients among populations. However, phenotypic divergence of sexually selected traits may also be influenced by other factors, such as natural selection and genetic constraints. Here, we document differences in male sexual traits among six introduced Australian populations of guppies and untangle the forces driving divergence in these sexually selected traits. Using an experimental approach, we found that male size, area of orange coloration, number of sperm per ejaculate and linear sexual selection gradients for male traits differed among populations. Within populations, a large mismatch between the direction of selection and male traits suggests that constraints may be important in preventing male traits from evolving in the direction of selection. Among populations, however, variation in sexual selection explained more than half of the differences in trait variation, suggesting that, despite within‐population constraints, sexual selection has contributed to population divergence of male traits. Differences in sexual traits were also associated with predation risk and neutral genetic distance. Our study highlights the importance of sexual selection in trait divergence in introduced populations, despite the presence of constraining factors such as predation risk and evolutionary history.     

Sexual dimorphism of head morphology in three‐spined stickleback Gasterosteus aculeatus

This study examined sexual dimorphism of head morphology in the ecologically diverse three‐spined stickleback Gasterosteus aculeatus. Male G. aculeatus had longer heads than female G. aculeatus in all 10 anadromous, stream and lake populations examined, and head length growth rates were significantly higher in males in half of the populations sampled, indicating that differences in head size increased with body size in many populations. Despite consistently larger heads in males, there was significant variation in size‐adjusted head length among populations, suggesting that the relationship between head length and body length was flexible. Inter‐population differences in head length were correlated between sexes, thus population‐level factors influenced head length in both sexes despite the sexual dimorphism present. Head shape variation between lake and anadromous populations was greater than that between sexes. The common divergence in head shape between sexes across populations was about twice as important as the sexual dimorphism unique to each population. Finally, much of the sexual dimorphism in head length was due to divergence in the anterior region of the head, where the primary trophic structures were found. It is unclear whether the sexual dimorphism was due to natural selection for niche divergence between sexes or sexual selection. This study improves knowledge of the magnitude, growth rate divergence, inter‐population variation and location of sexual dimorphism in G. aculeatus head morphology.     

Ecomorphological convergence in Eleutherodactylus frogs: a case of replicate radiations in the Caribbean

Replicate radiations, the repeated multiplication of species associated with ecological divergence, have attracted much attention and generated as much debate. Due to the few well‐studied cases, it remains unclear whether replicate radiations are an exceptional result of evolution or a relatively common example of the power of adaptation by natural selection. We examined the case of Eleutherodactylus frogs, which radiated in the Caribbean islands resulting in more than 160 species that occupy very diverse habitats. A time‐calibrated phylogeny revealed that these frogs independently diversified on all larger islands producing species that occupy a broad range of microhabitats in different islands. Using phylogenetic comparative methods, we found an association between morphological traits and particular microhabitats, and for most microhabitats detected significant morphological convergence. Our results indicate Caribbean Eleutherodactylus are a novel example of replicate radiations, and highlight the predictability of evolutionary processes, as similar ecological opportunities can lead to similar outcomes.     

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.     

Variable extent of parallelism in respiratory,circulatory, and neurological traits across lake whitefish species pairs

Parallel adaptive radiation events provide a powerful framework for investigations of ecology's contribution to phenotypic diversification. Ecologically driven divergence has been invoked to explain the repeated evolution of sympatric dwarf and normal lake whitefish (Coregonus clupeaformis) species in multiple lakes in eastern North America. Nevertheless, links between most putatively adaptive traits and ecological variation remain poorly defined within and among whitefish species pairs. Here, we examine four species pairs for variation in gill, heart, and brain size; three traits predicted to show strong phenotypic responses to ecological divergence. In each of the species pairs, normals exhibited larger body size standardized gills compared to dwarfs – a pattern that is suggestive of a common ecological driver of gill size divergence. Within lakes, the seasonal hypoxia experienced in the benthic environment is a likely factor leading to the requirement for larger gills in normals. Interestingly, the morphological pathways used to achieve larger gills varied between species pairs from Québec and Maine, which may imply subtle non‐parallelism in gill size divergence related to differences in genetic background. There was also a non‐significant trend toward larger hearts in dwarfs, the more active species of the two, whereas brain size varied exclusively among the lake populations. Taken together, our results suggest that the diversification of whitefish has been driven by parallel and non‐parallel ecological conditions across lakes. Furthermore, the phenotypic response to ecological variation may depend on genetic background of each population.