Thermal adaptation and ecological speciation

Ecological speciation is defined as the emergence of reproductive isolation as a direct or indirect consequence of divergent ecological adaptation. Several empirical examples of ecological speciation have been reported in the literature which very often involve adaptation to biotic resources. In this review, we investigate whether adaptation to different thermal habitats could also promote speciation and try to assess the importance of such processes in nature. Our survey of the literature identified 16 animal and plant systems where divergent thermal adaptation may underlie (partial) reproductive isolation between populations or may allow the stable coexistence of sibling taxa. In many of the systems, the differentially adapted populations have a parapatric distribution along an environmental gradient. Isolation often involves extrinsic selection against locally maladapted parental or hybrid genotypes, and additional pre- or postzygotic barriers may be important. Together, the identified examples strongly suggest that divergent selection between thermal environments is often strong enough to maintain a bimodal genotype distribution upon secondary contact. What is less clear from the available data is whether it can also be strong enough to allow ecological speciation in the face of gene flow through reinforcement-like processes. It is possible that intrinsic features of thermal gradients or the genetic basis of thermal adaptation make such reinforcement-like processes unlikely but it is equally possible that pertinent systems are understudied. Overall, our literature survey highlights (once again) the dearth of studies that investigate similar incipient species along the continuum from initial divergence to full reproductive isolation and studies that investigate all possible reproductive barriers in a given system.

The Role of Parasitism in Adaptive Radiations - When Might Parasites Promote and When Might They Constrain Ecological Speciation?

Research on speciation and adaptive radiation has flourished during the past decades, yet factors underlying initiation of reproductive isolation often remain unknown. Parasites represent important selective agents and have received renewed attention in speciation research. We review the literature on parasite-mediated divergent selection in context of ecological speciation and present empirical evidence for three nonexclusive mechanisms by which parasites might facilitate speciation: reduced viability or fecundity of immigrants and hybrids, assortative mating as a pleiotropic by-product of host adaptation, and ecologically-based sexual selection. We emphasise the lack of research on speciation continuums, which is why no study has yet made a convincing case for parasite driven divergent evolution to initiate the emergence of reproductive isolation. We also point interest towards selection imposed by single vs. multiple parasite species, conceptually linking this to strength and multifariousness of selection. Moreover, we discuss how parasites, by manipulating behaviour or impairing sensory abilities of hosts, may change the form of selection that underlies speciation. We conclude that future studies should consider host populations at variable stages of the speciation process, and explore recurrent patterns of parasitism and resistance that could pinpoint the role of parasites in imposing the divergent selection that initiates ecological speciation.

Adaptive divergence of ancient gene duplicates in the avian MHC class II beta.

Gene duplication and neofunctionalization are known to be important processes in the evolution of phenotypic complexity. They account for important evolutionary novelties that confer ecological adaptation, such as the major histocompatibility complex (MHC), a multigene family crucial to the vertebrate immune system. In birds, two MHC class II β (MHCIIβ) exon 3 lineages have been recently characterized, and two hypotheses for the evolutionary history of MHCIIβ lineages were proposed. These lineages could have arisen either by 1) an ancient duplication and subsequent divergence of one paralog or by 2) recent parallel duplications followed by functional convergence. Here, we compiled a data set consisting of 63 MHCIIβ exon 3 sequences from six avian orders to distinguish between these hypotheses and to understand the role of selection in the divergent evolution of the two avian MHCIIβ lineages. Based on phylogenetic reconstructions and simulations, we show that a unique duplication event preceding the major avian radiations gave rise to two ancestral MHCIIβ lineages that were each likely lost once later during avian evolution. Maximum likelihood estimation shows that following the ancestral duplication, positive selection drove a radical shift from basic to acidic amino acid composition of a protein domain facing the α-chain in the MHCII α β-heterodimer. Structural analyses of the MHCII α β-heterodimer highlight that three of these residues are potentially involved in direct interactions with the α-chain, suggesting that the shift following duplication may have been accompanied by coevolution of the interacting α- and β-chains. These results provide new insights into the long-term evolutionary relationships among avian MHC genes and open interesting perspectives for comparative and population genomic studies of avian MHC evolution.

Case studies and mathematical models of ecological speciation. 1. Cichlids in a crater lake.

A recent study of a pair of sympatric species of cichlids in Lake Apoyo in Nicaragua is viewed as providing probably one of the most convincing examples of sympatric speciation to date. Here, we describe and study a stochastic, individual-based, explicit genetic model tailored for this cichlid system. Our results show that relatively rapid (<20,000 generations) colonization of a new ecological niche and (sympatric or parapatric) speciation via local adaptation and divergence in habitat and mating preferences are theoretically plausible if: (i) the number of loci underlying the traits controlling local adaptation, and habitat and mating preferences is small; (ii) the strength of selection for local adaptation is intermediate; (iii) the carrying capacity of the population is intermediate; and (iv) the effects of the loci influencing nonrandom mating are strong. We discuss patterns and timescales of ecological speciation identified by our model, and we highlight important parameters and features that need to be studied empirically to provide information that can be used to improve the biological realism and power of mathematical models of ecological speciation.

Selection on a genetic polymorphism counteracts ecological speciation in a stick insect.

The interplay between selection and aspects of the genetic architecture of traits (such as linkage, dominance, and epistasis) can either drive or constrain speciation [1-3]. Despite accumulating evidence that speciation can progress to "intermediate" stages-with populations evolving only partial reproductive isolation-studies describing selective mechanisms that impose constraints on speciation are more rare than those describing drivers. The stick insect Timema cristinae provides an example of a system in which partial reproductive isolation has evolved between populations adapted to different host plant environments, in part due to divergent selection acting on a pattern polymorphism [4, 5]. Here, we demonstrate how selection on a green/melanistic color polymorphism counteracts speciation in this system. Specifically, divergent selection between hosts does not occur on color phenotypes because melanistic T. cristinae are cryptic on the stems of both host species, are resistant to a fungal pathogen, and have a mating advantage. Using genetic crosses and genome-wide association mapping, we quantify the genetic architecture of both the pattern and color polymorphism, illustrating their simple genetic control. We use these empirical results to develop an individual-based model that shows how the melanistic phenotype acts as a "genetic bridge" that increases gene flow between populations living on different hosts. Our results demonstrate how variation in the nature of selection acting on traits, and aspects of trait genetic architecture, can impose constraints on both local adaptation and speciation.

Establishing Winter Origins of Migrating Lesser Snow Geese Using Stable Isotopes

Increases in Snow Goose (Chen caerulescens) populations and large-scale habitat changes in North America have contributed to the concentration of migratory waterfowl on fewer wetlands, reducing resource availability, and enhancing risks of disease transmission. Predicting wintering locations of migratory individuals is critical to guide wildlife population management and habitat restoration. We used stable carbon (δ13C), nitrogen (δ15N), and hydrogen (δ2H) isotope ratios in muscle tissue of wintering Snow Geese to discriminate four major wintering areas, the Playa Lake Region, Texas Gulf Coast, Louisiana Gulf Coast, and Arkansas, and infer the wintering locations of individuals collected later during the 2007 and 2008 spring migrations in the Rainwater Basin (RWB) of Nebraska. We predicted the wintering ground derivation of migrating Snow Geese using a likelihood-based approach. Our three-isotope analysis provided an efficient discrimination of the four wintering areas. The assignment model predicted that 53% [95% CI: 37-69] of our sample of Snow Geese from the RWB in 2007 had most likely originated in Louisiana, 38% [23-54] had wintered on Texas Gulf Coast, and 9% [0-20] in Arkansas; the assessment suggested that 89% [73-100] of our 2008 sample had most likely come from Texas Gulf Coast, 9% [0-27] from Louisiana Gulf Coast, and 2% [0-9] from Arkansas. Further segregation of wintering grounds and additional sampling of spring migrating Snow Geese would refine overall assignment and help explain interannual variations in migratory connectivity. The ability to distinguish origins of northbound geese can support the development of spatially-adaptive management strategies for the midcontinent Snow Goose population. Establishing migratory connectivity using isotope assignment techniques can be extended to other waterfowl species to determine critical habitat, evaluate population energy requirements, and inform waterfowl conservation and management strategies.

Adaptive divergence and speciation among sexual and pseudoviviparous populations of Festuca

Pseudovivipary is an environmentally induced flowering abnormality in which vegetative shoots replace seminiferous (sexual) inflorescences. Pseudovivipary is usually retained in transplantation experiments, indicating that the trait is not solely induced by the growing environment. Pseudovivipary is the defining characteristic of Festuca vivipara, and arguably the only feature separating this species from its closest seminiferous relative, Festuca ovina. We performed phylogenetic and population genetic analysis on sympatric F. ovina and F. vivipara samples to establish whether pseudovivipary is an adaptive trait that accurately defines the separation of genetically distinct Festuca species. Chloroplast and nuclear marker-based analyses revealed that variation at a geographical level can exceed that between F. vivipara and F. ovina. We deduced that F. vivipara is a recent species that frequently arises independently within F. ovina populations and has not accumulated significant genetic differentiation from its progenitor. We inferred local gene flow between the species. We identified one amplified fragment length polymorphism marker that may be linked to a pseudovivipary-related region of the genome, and several other markers provide evidence of regional local adaptation in Festuca populations. We conclude that F. vivipara can only be appropriately recognized as a morphologically and ecologically distinct species; it lacks genetic differentiation from its relatives. This is the first report of a ‘failure in normal flowering development’ that repeatedly appears to be adaptive, such that the trait responsible for species recognition constantly reappears on a local basis.

Ecological speciation and phenotypic plasticity affect ecosystems

Phenotypic differences among closely related populations and species can cause contrasting effects on ecosystems; however, it is unknown whether such effects result from genetic divergence, phenotypic plasticity, or both. To test this, we reared sympatric limnetic and benthic species of whitefish from a young adaptive radiation in a common garden, where the benthic species was raised on two distinct food types. We then used these fish in a mesocosm experiment to test for contrasting ecosystem effects of closely related species and of plastically induced differences within a species. We found that strong contrasting ecosystem effects resulted more frequently from genetic divergence, although they were not stronger overall than those resulting from phenotypic plasticity. Overall, our results provide evidence that genetically based differences among closely related species that evolved during a young adaptive radiation can affect ecosystems, and that phenotypic plasticity can modify the ecosystem effects of such species.

Integrative analyses of speciation and divergence in Psammodromus hispanicus (Squamata: Lacertidae).

ABSTRACT: BACKGROUND: Genetic, phenotypic and ecological divergence within a lineage is the result of past and ongoing evolutionary processes, which lead ultimately to diversification and speciation. Integrative analyses allow linking diversification to geological, climatic, and ecological events, and thus disentangling the relative importance of different evolutionary drivers in generating and maintaining current species richness. RESULTS: Here, we use phylogenetic, phenotypic, geographic, and environmental data to investigate diversification in the Spanish sand racer (Psammodromus hispanicus). Phylogenetic, molecular clock dating, and phenotypic analyses show that P. hispanicus consists of three lineages. One lineage from Western Spain diverged 8.3 (2.9-14.7) Mya from the ancestor of Psammodromus hispanicus edwardsianus and P. hispanicus hispanicus Central lineage. The latter diverged 4.8 (1.5-8.7) Mya. Molecular clock dating, together with population genetic analyses, indicate that the three lineages experienced northward range expansions from southern Iberian refugia during Pleistocene glacial periods. Ecological niche modelling shows that suitable habitat of the Western lineage and P. h. edwardsianus overlap over vast areas, but that a barrier may hinder dispersal and genetic mixing of populations of both lineages. P. h. hispanicus Central lineage inhabits an ecological niche that overlaps marginally with the other two lineages. CONCLUSIONS: Our results provide evidence for divergence in allopatry and niche conservatism between the Western lineage and the ancestor of P. h. edwardsianus and P. h. hispanicus Central lineage, whereas they suggest that niche divergence is involved in the origin of the latter two lineages. Both processes were temporally separated and may be responsible for the here documented genetic and phenotypic diversity of P. hispanicus. The temporal pattern is in line with those proposed for other animal lineages. It suggests that geographic isolation and vicariance played an important role in the early diversification of the group, and that lineage diversification was further amplified through ecological divergence.

Trophic dynamics in the northern Humboldt Current system: insights from stable isotopes and stomach content analyses

The northern Humboldt Current system (NHCS) off Peru is one of the most productive world marine regions. It represents less than 0.1% of the world ocean surface but presently sustains about 10% of the world fish catch, with the Peruvian anchovy or anchoveta Engraulis ringens as emblematic fish resource. Compared with other eastern boundary upwelling systems, the higher fish productivity of the NHCS cannot be explained by a corresponding higher primary productivity. On another hand, the NHCS is the region where El Niño, and climate variability in general, is most notable. Also, surface oxygenated waters overlie an intense and extremely shallow Oxygen Minimum Zone (OMZ). In this context, the main objective of this study is to better understand the trophic flows in the NHCS using both stomach content and stable isotope analyses. The study focuses on a variety of organisms from low trophic levels such as zooplankton to top predators (seabirds and fur seals). The approach combines both long-term and specific studies on emblematic species such as anchoveta, and sardine Sardinops sagax and a more inclusive analysis considering the 'global' food web in the recent years (2008 – 2012) using stable isotope analysis. Revisiting anchovy and sardine we show that whereas phytoplankton largely dominated anchoveta and sardine diets in terms of numerical abundance, the carbon content of prey items indicated that zooplankton was by far the most important dietary component. Indeed for anchovy euphausiids contributed 67.5% of dietary carbon, followed by copepods (26.3%). Selecting the largest prey, the euphausiids, provide an energetic advantage for anchoveta in its ecosystem where oxygen depletion imposes strong metabolic constrain to pelagic fish. Sardine feed on smaller zooplankton than do anchoveta, with sardine diet consisting of smaller copepods and fewer euphausiids than anchoveta diet. Hence, trophic competition between sardine and anchovy in the northern Humboldt Current system is minimized by their partitioning of the zooplankton food resource based on prey size, as has been reported in other systems. These results suggest an ecological role for pelagic fish that challenges previous understanding of their position in the foodweb (zooplanktophagous instead of phytophagous), the functioning and the trophic models of the NHCS. Finally to obtain a more comprehensive vision of the relative trophic position of NHCS main components we used stable isotope analyses. For that purpose we analyzed the δ13C and δ15N stable isotope values of thirteen taxonomic categories collected off Peru from 2008 - 2011, i.e., zooplankton, fish, squids and air-breathing top predators. The δ15N isotope signature was strongly impacted by the species, the body length and the latitude. Along the Peruvian coast, the OMZ get more intense and shallow south of ~7.5ºS impacting the baseline nitrogen stable isotopes. Employing a linear mixed-effects modelling approach taking into account the latitudinal and body length effects, we provide a new vision of the relative trophic position of key ecosystem components. Also we confirm stomach content-based results on anchoveta Engraulis ringens and highlight the potential remarkable importance of an often neglected ecosystem component, the squat lobster Pleuroncodes monodon. Indeed, our results support the hypothesis according to which this species forage to some extent on fish eggs and larvae and can thus predate on the first life stages of exploited species. However, the δ13C values of these two species suggest that anchoveta and squat lobster do not exactly share the same habitat. This would potentially reduce some direct competition and/or predation.

Integrative analyses of speciation and divergence in Psammodromus hispanicus (Squamata: Lacertidae).

BackgroundGenetic, phenotypic and ecological divergence within a lineage is the result of past and ongoing evolutionary processes, which lead ultimately to diversification and speciation. Integrative analyses allow linking diversification to geological, climatic, and ecological events, and thus disentangling the relative importance of different evolutionary drivers in generating and maintaining current species richness.ResultsHere, we use phylogenetic, phenotypic, geographic, and environmental data to investigate diversification in the Spanish sand racer (Psammodromus hispanicus). Phylogenetic, molecular clock dating, and phenotypic analyses show that P. hispanicus consists of three lineages. One lineage from Western Spain diverged 8.3 (2.9-14.7) Mya from the ancestor of Psammodromus hispanicus edwardsianus and P. hispanicus hispanicus Central lineage. The latter diverged 4.8 (1.5-8.7) Mya. Molecular clock dating, together with population genetic analyses, indicate that the three lineages experienced northward range expansions from southern Iberian refugia during Pleistocene glacial periods. Ecological niche modelling shows that suitable habitat of the Western lineage and P. h. edwardsianus overlap over vast areas, but that a barrier may hinder dispersal and genetic mixing of populations of both lineages. P. h. hispanicus Central lineage inhabits an ecological niche that overlaps marginally with the other two lineages.ConclusionsOur results provide evidence for divergence in allopatry and niche conservatism between the Western lineage and the ancestor of P. h. edwardsianus and P. h. hispanicus Central lineage, whereas they suggest that niche divergence is involved in the origin of the latter two lineages. Both processes were temporally separated and may be responsible for the here documented genetic and phenotypic diversity of P. hispanicus. The temporal pattern is in line with those proposed for other animal lineages. It suggests that geographic isolation and vicariance played an important role in the early diversification of the group, and that lineage diversification was further amplified through ecological divergence.