Genomic signatures of relaxed disruptive selection associated with speciation reversal in whitefish

Background: Speciation reversal: the erosion of species differentiation via an increase in introgressive hybridization due to the weakening of previously divergent selection regimes, is thought to be an important, yet poorly understood, driver of biodiversity loss. Our study system, the Alpine whitefish (Coregonus spp.) species complex is a classic example of a recent postglacial adaptive radiation: forming an array of endemic lake flocks, with the independent origination of similar ecotypes among flocks. However, many of the lakes of the Alpine radiation have been seriously impacted by anthropogenic nutrient enrichment, resulting in a collapse in neutral genetic and phenotypic differentiation within the most polluted lakes. Here we investigate the effects of eutrophication on the selective forces that have shaped this radiation, using population genomics. We studied eight sympatric species assemblages belonging to five independent parallel adaptive radiations, and one species pair in secondary contact. We used AFLP markers, and applied FST outlier (BAYESCAN, DFDIST) and logistic regression analyses (MATSAM), to identify candidate regions for disruptive selection in the genome and their associations with adaptive traits within each lake flock. The number of outlier and adaptive trait associated loci identified per lake were then regressed against two variables (historical phosphorus concentration and contemporary oxygen concentration) representing the strength of eutrophication. Results: Whilst we identify disruptive selection candidate regions in all lake flocks, we find similar trends, across analysis methods, towards fewer disruptive selection candidate regions and fewer adaptive trait/candidate loci associations in the more polluted lakes. Conclusions: Weakened disruptive selection and a concomitant breakdown in reproductive isolating mechanisms in more polluted lakes has lead to increased gene flow between coexisting Alpine whitefish species. We hypothesize that the resulting higher rates of interspecific recombination reduce either the number or extent of genomic islands of divergence surrounding loci evolving under disruptive natural selection. This produces the negative trend seen in the number of selection candidate loci recovered during genome scans of whitefish species flocks, with increasing levels of anthropogenic eutrophication: as the likelihood decreases that AFLP restriction sites will fall within regions of heightened genomic divergence and therefore be classified as FST outlier loci. This study explores for the first time the potential effects of human-mediated relaxation of disruptive selection on heterogeneous genomic divergence between coexisting species.

Widespread phenotypic and genetic divergence along altitudinal gradients in animals

Altitudinal gradients offer valuable study systems to investigate how adap- tive genetic diversity is distributed within and between natural populations and which factors promote or prevent adaptive differentiation. The environ- mental clines along altitudinal gradients tend to be steep relative to the dispersal distance of many organisms, providing an opportunity to study the joint effects of divergent natural selection and gene flow. Temperature is one variable showing consistent altitudinal changes, and altitudinal gradi- ents can therefore provide spatial surrogates for some of the changes antici- pated under climate change. Here, we investigate the extent and patterns of adaptive divergence in animal populations along altitudinal gradients by sur- veying the literature for (i) studies on phenotypic variation assessed under common garden or reciprocal transplant designs and (ii) studies looking for signatures of divergent selection at the molecular level. Phenotypic data show that significant between-population differences are common and taxo- nomically widespread, involving traits such as mass, wing size, tolerance to thermal extremes and melanization. Several lines of evidence suggest that some of the observed differences are adaptively relevant, but rigorous tests of local adaptation or the link between specific phenotypes and fitness are sorely lacking. Evidence for a role of altitudinal adaptation also exists for a number of candidate genes, most prominently haemoglobin, and for anony- mous molecular markers. Novel genomic approaches may provide valuable tools for studying adaptive diversity, also in species that are not amenable to experimentation.

Repeated and predictable patterns of ecotypic differentiation during a biological invasion: lake-stream divergence in parapatric Swiss stickleback

The relative importance of ecological selection and geographical isolation in promoting and constraining genetic and phenotypic differentiation among populations is not always obvious. Interacting with divergent selection, restricted opportunity for gene flow may in some cases be as much a cause as a consequence of adaptation, with the latter being a hallmark of ecologi- cal speciation. Ecological speciation is well studied in parts of the native range of the three-spined stickleback. Here, we study this process in a recently invaded part of its range. Switzerland was colonized within the past 140 years from at least three different colonization events involving differ- ent stickleback lineages. They now occupy diverse habitats, ranging from small streams to the pelagic zone of large lakes. We use replicated systems of parapatric lake and stream populations, some of which trace their origins to different invasive lineages, to ask (i) whether phenotypic divergence occurred among populations inhabiting distinct habitats, (ii) whether trajec- tories of phenotypic divergence follow predictable parallel patterns and (iii) whether gene flow constrains divergent adaptation or vice versa. We find consistent phenotypic divergence between populations occupying distinct habitats. This involves parallel evolution in several traits with known eco- logical relevance in independent evolutionary lineages. Adaptive divergence supersedes homogenizing gene flow even at a small spatial scale. We find evidence that adaptive phenotypic divergence places constraints on gene flow over and above that imposed by geographical distance, signalling the early onset of ecological speciation.

Adaptive plasticity and genetic divergence in feeding efficiency during parallel adaptive radiation of whitefish (Coregonus spp.)

Parallel phenotypic divergence in replicated adaptive radiations could either result from parallel genetic divergence in response to similar divergent selec- tion regimes or from equivalent phenotypically plastic response to the repeated occurrence of contrasting environments. In post-glacial fish, repli- cated divergence in phenotypes along the benthic-limnetic habitat axis is commonly observed. Here, we use two benthic-limnetic species pairs of whitefish from two Swiss lakes, raised in a common garden design, with reciprocal food treatments in one species pair, to experimentally measure whether feeding efficiency on benthic prey has a genetic basis or whether it underlies phenotypic plasticity (or both). To do so, we offered experimental fish mosquito larvae, partially burried in sand, and measured multiple feed- ing efficiency variables. Our results reveal both, genetic divergence as well as phenotypically plastic divergence in feeding efficiency, with the pheno- typically benthic species raised on benthic food being the most efficient forager on benthic prey. This indicates that both, divergent natural selection on genetically heritable traits and adaptive phenotypic plasticity, are likely important mechanisms driving phenotypic divergence in adaptive radiation.

Experimental evidence for trait utility of gill raker number in adaptive radiation of a north temperate fish

North temperate fish in post-glacial lakes are textbook examples for rapid parallel adaptive radiation into multiple trophic specialists within individual lakes. Speciation repeatedly proceeded along the benthic – limnetic habitat axis, and benthic – limnetic sister species diverge in the number of gill rakers. Yet, the utility of different numbers of gill rakers for consuming benthic vs. limnetic food has only very rarely been experimentally demonstrated. We bred and raised families of a benthic – limnetic species pair of whitefish under common garden conditions to test whether these species (i) show heritable differentiation in feeding efficiency on zooplankton, and (ii) whether varia- tion in feeding efficiency is predicted by variation in gill raker numbers. We used zooplankton of three different size classes to investigate prey size dependency of divergence in feeding efficiency and to investigate the effect strength of variation in the number of gill rakers. Our results show strong interspecific differences in feeding efficiency. These differences are largest when fish were tested with the smallest zooplankton. Importantly, feeding efficiency is significantly positively correlated with the number of gill rakers when using small zooplankton, also when species identity is statistically controlled for. Our results support the hypothesis that a larger number of gill rakers are of adaptive significance for feeding on zooplankton and pro- vide one of the first experimental demonstrations of trait utility of gill raker number when fish feed on zooplankton. These results are consistent with the suggested importance of divergent selection driven feeding adaptation during adaptive radiation of fish in post-glacial lakes.

Sexual dimorphism dominates divergent host plant use in stick insect trophic morphology

Background: Clear examples of ecological speciation exist, often involving divergence in trophic morphology. However, substantial variation also exists in how far the ecological speciation process proceeds, potentially linked to the number of ecological axes, traits, or genes subject to divergent selection. In addition, recent studies highlight how differentiation might occur between the sexes, rather than between populations. We examine variation in trophic morphology in two host-plant ecotypes of walking-stick insects (Timema cristinae), known to have diverged in morphological traits related to crypsis and predator avoidance, and to have reached an intermediate point in the ecological speciation process. Here we test how host plant use, sex, and rearing environment affect variation in trophic morphology in this species using traditional multivariate, novel kernel density based and Bayesian morphometric analyses. Results: Contrary to expectations, we find limited host-associated divergence in mandible shape. Instead, the main predictor of shape variation is sex, with secondary roles of population of origin and rearing environment. Conclusion: Our results show that trophic morphology does not strongly contribute to host-adapted ecotype divergence in T. cristinae and that traits can respond to complex selection regimes by diverging along different intraspecific lines, thereby impeding progress toward speciation.

Relaxed trait covariance in interspecific cichlid hybrids predicts morphological diversity in adaptive radiations

The process of adaptive radiation involves multiple events of speciation in short succession, associated with ecological diversification. Understanding this process requires identifying the origins of heritable phenotypic variation that allows adaptive radiation to progress. Hybridization is one source of genetic and morphological variation that may spur adaptive radiation. We experimentally explored the potential role of hybridization in facilitating the onset of adaptive radiation. We generated first- and second-generation hybrids of four species of African cichlid fish, extant relatives of the putative ancestors of the adaptive radiations of Lakes Victoria and Malawi. We com- pared patterns in hybrid morphological variation with the variation in the lake radiations. We show that significant fractions of the interspecific mor- phological variation and the major trajectories in morphospace that charac- terize whole radiations can be generated in second-generation hybrids. Furthermore, we show that covariation between traits is relaxed in second- generation hybrids, which may facilitate adaptive diversification. These results support the idea that hybridization can provide the heritable pheno- typic diversity necessary to initiate adaptive radiation.

Habitat-dependent geographical variation in ontogenetic allometry of the shiner perch Cymatogaster aggregata Gibbons (Teleostei: Embiotocidae)

Studies of intraspecific morphological variation in fishes have traditionally focused on freshwater rather than marine species. In addition, such studies typically focus on adults, although causes and intensities of selective pressures most likely vary through an individual’s lifetime. In this study, body and head shape of a marine species, shiner perch Cymatogaster aggregata Gibbons were compared among localities along the Pacific Northwest coast of North America. Evidence was found for intraspecific variation in ontogenetic allometry, and for a closer correlation of body shape with environment rather than geographical proximity. This correlation with environment was more evident in younger fish, thereby demonstrating the importance of analysing multiple life stages. A common garden experiment suggests both environmental and genetic bases for the observed differences. Recognizing intraspecific ecomorphological complexity and its specificity to habitat and/or life stage can have important consequences for understanding the role of local adaptation and population dynamics in macroecology.

Scanning electron microscope study of the developing microvasculature in the postnatal rat lung.

During postnatal growth the parenchymal septa of rat lung undergo an impressive restructuring. While immature septa are thick and contain two capillary layers, mature septa are slender and contain a single microvascular network. Using the Mercox casting technique and scanning electron microscopy, we investigated the mode and the timing of the transformation of the pulmonary capillary bed. During the third postnatal week the parenchymal septa rapidly mature to match adult morphology. Even in adult lungs, however, remnants of the immature status are present: A capillary bilayer is regularly found at the base and the tip of the septa. Our observations support the concept that reduction of intervening tissue, partial fusion of the two capillary networks, and preferential growth lead to the mature vascular arrangement. The fact that true mature interalveolar septa show a denser capillary network than alveolar walls abutting onto pleura, bronchi, or larger vessels is consonant with the fusion theory. Towards the nonparenchyma, the capillary network surrounding every airspace had no counterpart to fuse with. From quantitative data it can be calculated that owing to lung growth, mesh size should increase more than four times between birth and adult age. The adult lung network, however, is denser than the one in young animals. This means that new meshes must be added during growth. We propose that small holes observed in sheet-like regions of the microvasculature enlarge to form new capillary meshes. With this mechanism of in-itself or intussusceptional growth, sprouting of individual capillary segments to increase network size is no longer needed.

Positive effect of the yellow morph on female reproductive success in the flower colour polymorphic Iris lutescens (Iridaceae), a deceptive species

The deceptive Iris lutescens (Iridaceae) shows a heritable and striking flower colour polymorphism, with both yellow- and purple-flowered individuals growing sympatrically. Deceptive species with flower colour polymorphism are mainly described in the family Orchidaceae and rarely found in other families. To explain the maintenance of flower colour polymorphism in I.lutescens, we investigated female reproductive success in natural populations of southern France, at both population and local scales (within populations). Female reproductive success was positively correlated with yellow morph frequency, at both the population scale and the local scale. Therefore, we failed to observe negative frequency-dependent selection (NFDS), a mechanism commonly invoked to explain flower colour polymorphism in deceptive plant species. Flower size and local flower density could also affect female reproductive success in natural populations. Pollinator behaviour could explain the positive effect of the yellow morph, and our results suggest that flower colour polymorphism might not persist in I.lutescens, but alternative explanations not linked to pollinator behaviour are discussed. In particular, NFDS, although an appealingly simple explanation previously demonstrated in orchids, may not always contribute to maintaining flower colour polymorphism, even in deceptive species.

Horn Growth and Reproduction in a Long-Lived Male Mammal: No Compensation for Poor Early-Life Horn Growth

Secondary sexual traits in males of polygynous species are important determinants of reproductive success. It is, however, unknown if and how the development of continuously growing traits at different life-stages is related to reproduction in long-lived male mammals. In this study, we evaluated the relationship of early and late horn growth on social status and reproduction in long-lived male Alpine ibex (Capra ibex). For this, we analysed individual horn growth and assessed its effect on dominance and reproduction. No evidence was detected for compensatory horn growth, as late-life horn growth positively depended on early-life horn growth in males. Still, individuals with longer horn segments grown during early adulthood experienced a stronger age-dependent length decline in annual horn growth during the late development. Accordingly, a divergence between individual growth potential and realized horn growth late in life has to be assumed. Residual age-specific horn length and length of early grown horn segments both positively affected dominance and reproductive success, whereas, contrary to our expectation, no significant effect of the length of horn segments grown during the late development was detected. Suspected higher somatic costs incurred by high-quality males during their late development might at least partly be responsible for this finding. Overall, our study suggests that the total length of horns and their early development in long-lived male Alpine ibex is a reliable indicator of reproductive success and that individuals may be unable to compensate for poor early-life growth performance at a later point in life.

Contemporary ecotypic divergence during a recent range expansion was facilitated by adaptive introgression

Although rapid phenotypic evolution during range expansion associated with colonization of contrasting habitats has been documented in several taxa, the evolutionary mechanisms that underlie such phenotypic divergence have less often been investigated. A strong candidate for rapid ecotype formation within an invaded range is the three-spine stickleback in the Lake Geneva region of central Europe. Since its introduction only about 140 years ago, it has undergone a significant expansion of its range and its niche, now forming phenotypically differentiated parapatric ecotypes that occupy either the pelagic zone of the large lake or small inlet streams, respectively. By comparing museum collections from different times with contemporary population samples, we here reconstruct the evolution of parapatric phenotypic divergence through time. Using genetic data from modern samples, we infer the underlying invasion history. We find that parapatric habitat-dependent phenotypic divergence between the lake and stream was already present in the first half of the twentieth century, but the magnitude of differentiation increased through time, particularly in antipredator defence traits. This suggests that divergent selection between the habitats occurred and was stable through much of the time since colonization. Recently, increased phenotypic differentiation in antipredator defence traits likely results from habitat-dependent selection on alleles that arrived through introgression from a distantly related lineage from outside the Lake Geneva region. This illustrates how hybridization can quickly promote phenotypic divergence in a system where adaptation from standing genetic variation was constrained.

Quick divergence but slow convergence during ecotype formation in lake and stream stickleback pairs of variable age

When genetic constraints restrict phenotypic evolution, diversification can be predicted to evolve along so-called lines of least resistance. To address the importance of such constraints and their resolution, studies of parallel phenotypic divergence that differ in their age are valuable. Here, we investigate the parapatric evolution of six lake and stream threespine stickleback systems from Iceland and Switzerland, ranging in age from a few decades to several millennia. Using phenotypic data, we test for parallelism in ecotypic divergence between parapatric lake and stream populations and compare the observed patterns to an ancestral-like marine population. We find strong and consistent phenotypic divergence, both among lake and stream populations and between our freshwater populations and the marine population. Interestingly, ecotypic divergence in low-dimensional phenotype space (i.e. single traits) is rapid and seems to be often completed within 100 years. Yet, the dimensionality of ecotypic divergence was highest in our oldest systems and only there parallel evolution of unrelated ecotypes was strong enough to overwrite phylogenetic contingency. Moreover, the dimensionality of divergence in different systems varies between trait complexes, suggesting different constraints and evolutionary pathways to their resolution among freshwater systems.