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.     

Gene expression remodelling and immune response during adaptive divergence in an African cichlid fish

Variation in gene expression contributes to ecological speciation by facilitating population persistence in novel environments. Likewise, immune responses can be of relevance in speciation driven by adaptation to different environments. Previous studies examining gene expression differences between recently diverged ecotypes have often relied on only one pair of populations, targeted the expression of only a subset of genes or used wild‐caught individuals. Here, we investigated the contribution of habitat‐specific parasites and symbionts and the underlying immunological abilities of ecotype hosts to adaptive divergence in lake–river population pairs of the cichlid fish Astatotilapia burtoni. To shed light on the role of phenotypic plasticity in adaptive divergence, we compared parasite and microbiota communities, immune response, and gene expression patterns of fish from natural habitats and a lake‐like pond set‐up. In all investigated population pairs, lake fish were more heavily parasitized than river fish, in terms of both parasite taxon composition and infection abundance. The innate immune response in the wild was higher in lake than in river populations and was elevated in a river population exposed to lake parasites in the pond set‐up. Environmental differences between lake and river habitat and their distinct parasite communities have shaped differential gene expression, involving genes functioning in osmoregulation and immune response. Most changes in gene expression between lake and river samples in the wild and in the pond set‐up were based on a plastic response. Finally, gene expression and bacterial communities of wild‐caught individuals and individuals acclimatized to lake‐like pond conditions showed shifts underlying adaptive phenotypic plasticity.     

Selection from parasites favours immunogenetic diversity but not divergence among locally adapted host populations

The unprecedented polymorphism in the major histocompatibility complex (MHC) genes is thought to be maintained by balancing selection from parasites. However, do parasites also drive divergence at MHC loci between host populations, or do the effects of balancing selection maintain similarities among populations? We examined MHC variation in populations of the livebearing fish Poecilia mexicana and characterized their parasite communities. Poecilia mexicana populations in the Cueva del Azufre system are locally adapted to darkness and the presence of toxic hydrogen sulphide, representing highly divergent ecotypes or incipient species. Parasite communities differed significantly across populations, and populations with higher parasite loads had higher levels of diversity at class II MHC genes. However, despite different parasite communities, marked divergence in adaptive traits and in neutral genetic markers, we found MHC alleles to be remarkably similar among host populations. Our findings indicate that balancing selection from parasites maintains immunogenetic diversity of hosts, but this process does not promote MHC divergence in this system. On the contrary, we suggest that balancing selection on immunogenetic loci may outweigh divergent selection causing divergence, thereby hindering host divergence and speciation. Our findings support the hypothesis that balancing selection maintains MHC similarities among lineages during and after speciation (trans‐species evolution).     

Differentiation in parasitism among ecotypes of whitefish segregating along depth gradients

Parasitism is a potential mechanism initiating or facilitating ecotypic differentiation and speciation in freshwater fish. While recent studies have begun to explore this question, there are no empirical studies of parasitism in evolutionary replicates of ecotype‐pairs at variable stages of speciation. Such comparative studies of parasitism along continuums of host differentiation are needed as a first step towards testing the role of parasites in ecological speciation. We explored parasitism of whitefish Coregonus lavaretus in four pre‐alpine lakes in Switzerland that hold replicate species radiations of whitefish. We sampled shallow and deep‐spawning ecotypes on their breeding grounds. We found significant and consistent differences in infection between the ecotypes so that the shallow‐spawning fish had more trematode infections, whereas the deepspawning fish had more cestodes. The magnitude of these differences correlated positively with the degree of the genetic differentiation among the ecotypes and negatively with the extent of eutrophication of the lakes. Although the overall diversity of infections was low, some parasite species with potential effects on fish showed marked differences in infection between the ecotypes, suggesting that parasitism may have a role in maintaining ecotype differentiation in this system. Our results also indicate previously unknown habitat segregation of the better differentiated ecotypes, i.e. species, along the depth gradient outside the breeding season. Moreover, oligotrophic lakes tended to have higher parasite species richness and higher abundances of infection, than mesotrophic and eutrophic lakes, suggesting that the history of eutrophication affects parasite diversity.     

Divergent Macroparasite Infections in Parapatric Swiss Lake-Stream Pairs of Threespine Stickleback (Gasterosteus aculeatus)

Spatial heterogeneity in diversity and intensity of parasitism is a typical feature of most host-parasite interactions, but understanding of the evolutionary implications of such variation is limited. One possible outcome of infection heterogeneities is parasite-mediated divergent selection between host populations, ecotypes or species which may facilitate the process of ecological speciation. However, very few studies have described infections in population-pairs along the speciation continuum from low to moderate or high degree of genetic differentiation that would address the possibility of parasite-mediated divergent selection in the early stages of the speciation process. Here we provide an example of divergent parasitism in freshwater fish ecotypes by examining macroparasite infections in threespine stickleback (Gasterosteus aculeatus) of four Swiss lake systems each harbouring parapatric lake-stream ecotype pairs. We demonstrate significant differences in infections within and between the pairs that are driven particularly by the parasite taxa transmitted to fish from benthic invertebrates. The magnitude of the differences tended to correlate positively with the extent of neutral genetic differentiation between the parapatric lake and stream populations of stickleback, whereas no such correlation was found among allopatric populations from similar or contrasting habitats. This suggests that genetic differentiation is unrelated to the magnitude of parasite infection contrasts when gene flow is constrained by geographical barriers while in the absence of physical barriers, genetic differentiation and the magnitude of differences in infections tend to be positively correlated.     

Lack of parallel genetic patterns underlying the repeated ecological divergence of beach and stream‐spawning kokanee salmon

Recent progress in methods for detecting adaptive population divergence in situ shows promise for elucidating the conditions under which selection acts to generate intraspecific diversity. Rapid ecological diversification is common in fishes; however, the role of phenotypic plasticity and adaptation to local environments is poorly understood. It is now possible to investigate genetic patterns to make inferences regarding phenotypic traits under selection and possible mechanisms underlying ecotype divergence, particularly where similar novel phenotypes have arisen in multiple independent populations. Here, we employed a bottom‐up approach to test for signatures of directional selection associated with divergence of beach‐ and stream‐spawning kokanee, the obligate freshwater form of sockeye salmon (Oncorhynchus nerka). Beach‐ and stream‐spawners co‐exist in many post‐glacial lakes and exhibit distinct reproductive behaviours, life‐history traits and spawning habitat preferences. Replicate ecotype pairs across five lakes in British Columbia, Canada were genotyped at 57 expressed sequence tag‐linked and anonymous microsatellite loci identified in a previous genome scan. Fifteen loci exhibited signatures of directional selection (high F_ST outliers), four of which were identified in multiple lakes. However, the lack of parallel genetic patterns across all lakes may be a result of: 1) an inability to detect loci truly under selection; 2) alternative genetic pathways underlying ecotype divergence in this system; and/or 3) phenotypic plasticity playing a formative role in driving kokanee spawning habitat differences. Gene annotations for detected outliers suggest pathogen resistance and energy metabolism as potential mechanisms contributing to the divergence of beach‐ and stream‐spawning kokanee, but further study is required.     

Long‐term coevolution between avian brood parasites and their hosts

Coevolutionary theory predicts that the most common long‐term outcome of the relationships between brood parasites and their hosts should be coevolutionary cycles based on a dynamic change selecting the currently least‐defended host species, given that when well‐defended hosts are abandoned, hosts will be selected to decrease their defences as these are usually assumed to be costly. This is assumed to be the case also in brood parasite‐host systems. Here I examine the frequency of the three potential long‐term outcomes of brood parasite–host coevolution (coevolutionary cycles, lack of rejection, and successful resistance) in 182 host species. The results of simple exploratory comparisons show that coevolutionary cycles are very scarce while the lack of rejection and successful resistance, which are considered evolutionary enigmas, are much more frequent. I discuss these results considering (i) the importance of different host defences at all stages of the breeding cycle, (ii) the role of phenotypic plasticity in long‐term coevolution, and (iii) the evolutionary history of host selection. I suggest that in purely antagonistic coevolutionary interactions, such as those involving brood parasites and their hosts, that although cycles will exist during an intermediate phase of the interactions, the arms race will end with the extinction of the host or with the host acquiring successful resistance. As evolutionary time passes, this resistance will force brood parasites to use previously less suitable host species. Furthermore, I present a model that represents the long‐term trajectories and outcomes of coevolutionary interactions between brood parasites and their hosts with respect to the evolution of egg‐rejection defence. This model suggests that as an increasing number of species acquire successful resistance, other unparasitized host species become more profitable and their parasitism rate and the costs imposed by brood parasitism at the population level will increase, selecting for the evolution of host defences. This means that although acceptance is adaptive when the parasitism rate and the costs of parasitism are very low, this cannot be considered to represent an evolutionary equilibrium, as conventional theory has done to date, because it is not stable.     

Parasitism, immune response and reproductive success in the house martin Delichonurbica

Parasites often exert strong selection pressures on their hosts that have evolved anti-parasite defences to counter the negative effects of parasites. We studied the relationship between intensity of parasitism, one aspect of host immune response, and host reproductive success, using the house martin bug Oeciacushirundinis and its house martin Delichonurbica host as a model system. Experimental manipulation of parasite load of nests during laying of the first clutch altered the intensity of parasitism. Parasites reduced the reproductive success of their hosts measured in terms of body condition and survival of nestlings. Host immune response, measured as the concentration of gammaglobulins and total plasma proteins, was positively associated with parasite reproduction, estimated as the number of juvenile parasites, but was only weakly related to the intensity of adult parasites. The concentration of gammaglobulins was negatively related to nestling body mass, implying a trade-off between immune function and body condition. Parasite reproduction thus exerts a cost on hosts by increasing the immune response. Received: 25 August 1997 / Accepted: 3 November 1997     

The evolution of acceptance and tolerance in hosts of avian brood parasites

Avian brood parasites lay their eggs in the nests of their hosts, which rear the parasite's progeny. The costs of parasitism have selected for the evolution of defence strategies in many host species. Most research has focused on resistance strategies, where hosts minimize the number of successful parasitism events using defences such as mobbing of adult brood parasites or rejection of parasite eggs. However, many hosts do not exhibit resistance. Here we explore why some hosts accept parasite eggs in their nests and how this is related to the virulence of the parasite. We also explore the extent to which acceptance of parasites can be explained by the evolution of tolerance; a strategy in which the host accepts the parasite but adjusts its life history or other traits to minimize the costs of parasitism. We review examples of tolerance in hosts of brood parasites (such as modifications to clutch size and multi‐broodedness), and utilize the literature on host–pathogen interactions and plant herbivory to analyse the prevalence of each type of defence (tolerance or resistance) and their evolution. We conclude that (i) the interactions between brood parasites and their hosts provide a highly tractable system for studying the evolution of tolerance, (ii) studies of host defences against brood parasites should investigate both resistance and tolerance, and (iii) tolerance and resistance can lead to contrasting evolutionary scenarios.     

Genomics of Rapid Incipient Speciation in Sympatric Threespine Stickleback

Ecological speciation is the process by which reproductively isolated populations emerge as a consequence of divergent natural or ecologically-mediated sexual selection. Most genomic studies of ecological speciation have investigated allopatric populations, making it difficult to infer reproductive isolation. The few studies on sympatric ecotypes have focused on advanced stages of the speciation process after thousands of generations of divergence. As a consequence, we still do not know what genomic signatures of the early onset of ecological speciation look like. Here, we examined genomic differentiation among migratory lake and resident stream ecotypes of threespine stickleback reproducing in sympatry in one stream, and in parapatry in another stream. Importantly, these ecotypes started diverging less than 150 years ago. We obtained 34,756 SNPs with restriction-site associated DNA sequencing and identified genomic islands of differentiation using a Hidden Markov Model approach. Consistent with incipient ecological speciation, we found significant genomic differentiation between ecotypes both in sympatry and parapatry. Of 19 islands of differentiation resisting gene flow in sympatry, all were also differentiated in parapatry and were thus likely driven by divergent selection among habitats. These islands clustered in quantitative trait loci controlling divergent traits among the ecotypes, many of them concentrated in one region with low to intermediate recombination. Our findings suggest that adaptive genomic differentiation at many genetic loci can arise and persist in sympatry at the very early stage of ecotype divergence, and that the genomic architecture of adaptation may facilitate this.     

Interaction of a host plant and its holoparasite: effects of previous selection by the parasite

If parasites decrease the fitness of their hosts one could expect selection for host traits (e.g. resistance and tolerance) that decrease the negative effects of parasitic infection. To study selection caused by parasitism, we used a novel study system: we grew host plants (Urtica dioica) that originated from previously parasitized and unparasitized natural populations (four of each) with or without a holoparasitic plant (Cuscuta europaea). Infectivity of the parasite (i.e. qualitative resistance of the host) did not differ between the two host types. Parasites grown with hosts from parasitized populations had lower performance than parasites grown with hosts from unparasitized populations, indicating host resistance in terms of parasite’s performance (i.e. quantitative resistance). However, our results suggest that the tolerance of parasitic infection was lower in hosts from parasitized populations compared with hosts from unparasitized populations as indicated by the lower above‐ground vegetative biomass of the infected host plants from previously parasitized populations.     

Flexible mate choice may contribute to ecotype assortative mating in pumpkinseed sunfish (Lepomis gibbosus)

Gene flow is expected to limit adaptive divergence, but the ecological and behavioural factors that govern gene flow are still poorly understood, particularly at the earliest stages of population divergence. Reduced gene flow through mate choice (sexual isolation) can evolve even under conditions of subtle population divergence if intermediate phenotypes have reduced fitness. We indirectly tested the hypothesis that mate choice has evolved between coexisting littoral and pelagic ecotypes of polyphenic pumpkinseed sunfish (Lepomis gibbosus) that have diverged in morphology and resource use and where intermediate phenotypes have reduced performance. We assessed the ecotype of nesting males and females using stable isotope estimates of diet and a divergent male morphological trait, oral jaw width. We found positive assortative mating between ecotypes in a common spawning habitat along exposed lake shorelines, but contrary to expectations, assortative mating was variably expressed between two sampling years. Although the factors that influence variable assortative mating remain unclear, our results are consistent with mate choice being expressed by ecotypes. Despite being variably expressed, mate choice will reduce gene flow between ecotypes and could contribute to further adaptive divergence depending on its frequency and strength in the population. Our findings add to a growing body of evidence indicating mate choice behaviour can be a plastic trait, an idea that should be more explicitly considered in empirical studies of mate choice as well as conceptual frameworks of mate choice evolution and adaptive divergence.     

Ecology of a diplozoon parasite on the gills of the African cyprinid Barbus neumayeri

In oxygen‐deficient waters, the difficulties of oxygen uptake in gill parasites and their fish hosts may influence host and parasite densities, site selection by the parasite, and effects of the parasite on host condition. This study quantified the prevalence and intensity of the gill monogenean Neodiplozoon polycotyleus in the African cyprinid fish Barbus neumayeri from an intermittent forest stream in western Uganda. Oxygen levels were low in the stream over the 12‐month study, averaging only 2.5 mg litre^?1 (monthly range = 1.2–4.3 mg litre^?1). However, parasite prevalence was high (47.2%), suggesting high tolerance to low oxygen in N. polycotyleus. The prevalence of parasites varied with host body size, with the highest frequency of occurrence in the middle size classes. Prevalence also varied over the year; seasonal peaks of rainfall coincided with a lower frequency of N. polycotyleus. The significantly nonrandom frequency distribution of parasites among hosts suggests regulation of parasite numbers. Of the hosts infected, 37.1% harboured one N. polycotyleus parasite, and 62.9% harboured two parasites. No fish were infected with more than two diplozoons. There was evidence for strong site specificity by N. polycotyleus within hosts; 77.7% of the parasites were located on the filaments of the second gill arch, which may relate to increased oxygen availability. In addition, only one of the 178 infected fish had more than one parasite on one side of the branchial basket. Although N. polycotyleus is undoubtedly parasitic, we found no evidence of a negative parasitic effect on the condition or reproductive status of B. neumayeri.     

Higher parasite resistance in Daphnia populations with recent epidemics

Natural populations often show genetic variation in parasite resistance, forming the basis for evolutionary response to selection imposed by parasitism. We investigated whether previous epidemics selected for higher resistance to novel parasite isolates in a Daphnia galeata–microparasite system by comparing susceptibility of host clones from populations with varying epidemic history. We manipulated resource availability to evaluate whether diet influences Daphnia susceptibility as epidemics are common in nutrient‐rich lakes. Exposing clones from 10 lakes under two food treatments to an allopatric protozoan parasite, we found that Daphnia originating from lakes (mainly nutrient rich) with previous epidemics better resist infection. Despite this result, there was a tendency of higher susceptibility in the low food treatment, suggesting that higher resistance of clones from populations with epidemic background is not directly caused by lake nutrient level. Rather, our results imply that host populations respond to parasite‐mediated selection by evolving higher parasite resistance.     

Genomic architecture of habitat‐related divergence and signature of directional selection in the body shapes of Gnathopogon fishes

Evolution of ecomorphologically relevant traits such as body shapes is important to colonize and persist in a novel environment. Habitat‐related adaptive divergence of these traits is therefore common among animals. We studied the genomic architecture of habitat‐related divergence in the body shape of Gnathopogon fishes, a novel example of lake–stream ecomorphological divergence, and tested for the action of directional selection on body shape differentiation. Compared to stream‐dwelling Gnathopogon elongatus, the sister species Gnathopogon caerulescens, exclusively inhabiting a large ancient lake, had an elongated body, increased proportion of the caudal region and small head, which would be advantageous in the limnetic environment. Using an F₂ interspecific cross between the two Gnathopogon species (195 individuals), quantitative trait locus (QTL) analysis with geometric morphometric quantification of body shape and restriction‐site associated DNA sequencing‐derived markers (1622 loci) identified 26 significant QTLs associated with the interspecific differences of body shape‐related traits. These QTLs had small to moderate effects, supporting polygenic inheritance of the body shape‐related traits. Each QTL was mostly located on different genomic regions, while colocalized QTLs were detected for some ecomorphologically relevant traits that are proxy of body and caudal peduncle depths, suggesting different degree of modularity among traits. The directions of the body shape QTLs were mostly consistent with the interspecific difference, and QTL sign test suggested a genetic signature of directional selection in the body shape divergence. Thus, we successfully elucidated the genomic architecture underlying the adaptive changes of the quantitative and complex morphological trait in a novel system.     

Relationship Between Nest Predation Suffered by Hosts and Brown-Headed Cowbird Parasitism: A Comparative Study

There are at least four main hypotheses that may explain how the evolution of host selection by avian brood parasites could be linked to nest predation among their potential hosts. First, selection may have favoured parasite phenotypes discriminating among hosts on the basis of expected nest failure. Second, parasitized nests may be more easily detected by predators and extra costs of parasitism may accelerate the evolution of host defences. Third, selection may have favoured predator phenotypes avoiding parasitized nests because parasitism enhances nest defence. Fourth, female brood parasites may directly or indirectly induce host nesting failures in order to enhance future laying opportunities. We collected data on brood parasitism and nest failure due to predation to test these hypotheses in a comparative approach using North American passerines and their brood parasite, the brown-headed cowbird Molothrus ater. Under the hypotheses 1 or 3 we predicted brood parasitism to be negatively associated with nest predation across species, whereas this relation is expected to be positive if hypotheses 2 or 4 are true. We demonstrate that independent of host suitability, nest location, habitat type, length of the nestling period, body mass and similarity among species due to common ancestry, species experiencing relatively high levels of nest predation suffered lower levels of cowbird parasitism. Our results suggest a previously ignored role for nest predation suffered by hosts on the dynamics of the coevolutionary relationships between hosts and avian brood parasites. Co-ordinating editor: Dr. F. Stuefer     

An experimental test of antagonistic effects of competition and parasitism on host performance in semi‐natural mesocosms

The mechanisms by which parasites can mediate the interactions between species have received increased interest in recent years. Nonetheless, most research has focused on the role of shared parasites as mediators of interspecific competition. Here, we explore the relative effects of Gyrodactylus specialist ectoparasites of Trinidadian guppies Poecilia reticulata on competition between their host and juveniles of the killifish Rivulus hartii. In mesocosms that replicate natural streams, we exposed guppies to only competitors, to only parasites, to both parasites and competitors, or the absence of both. Consistent with previous studies, we found that female guppies grew significantly less where only Gyrodactylus were present, and this was regardless of infection status or parasite load. Surprisingly, this effect of Gyrodactylus on the growth of female guppies was greatly reduced when both parasites and competitors were present in the mesocosms. We conclude that guppies can mediate the effects of Gyrodactylus on competition with Rivulus, by adaptively fine‐tuning their phenotype when simultaneously facing multiple enemies.     

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.     

Genomic evidence for population‐specific responses to co‐evolving parasites in a New Zealand freshwater snail

Reciprocal co‐evolving interactions between hosts and parasites are a primary source of strong selection that can promote rapid and often population‐ or genotype‐specific evolutionary change. These host–parasite interactions are also a major source of disease. Despite their importance, very little is known about the genomic basis of co‐evolving host–parasite interactions in natural populations, especially in animals. Here, we use gene expression and sequence evolution approaches to take critical steps towards characterizing the genomic basis of interactions between the freshwater snail Potamopyrgus antipodarum and its co‐evolving sterilizing trematode parasite, Microphallus sp., a textbook example of natural coevolution. We found that Microphallus‐infected P. antipodarum exhibit systematic downregulation of genes relative to uninfected P. antipodarum. The specific genes involved in parasite response differ markedly across lakes, consistent with a scenario where population‐level co‐evolution is leading to population‐specific host–parasite interactions and evolutionary trajectories. We also used an F_ST‐based approach to identify a set of loci that represent promising candidates for targets of parasite‐mediated selection across lakes as well as within each lake population. These results constitute the first genomic evidence for population‐specific responses to co‐evolving infection in the P. antipodarum‐Microphallus interaction and provide new insights into the genomic basis of co‐evolutionary interactions in nature.