The genetic basis of divergent melanic pigmentation in benthic and limnetic threespine stickleback.

Pigmentation is an excellent trait to examine patterns of evolutionary change because it is often under natural selection. Benthic and limnetic threespine stickleback (Gasterosteus aculeatus) exhibit distinct pigmentation phenotypes, likely an adaptation to occupation of divergent niches. The genetic architecture of pigmentation in vertebrates appears to be complex. Prior QTL mapping of threespine stickleback pigmentation phenotypes has identified several candidate loci. However-relative to other morphological phenotypes (e.g., spines or lateral plates)-the genetic architecture of threespine stickleback pigmentation remains understudied. Here, we performed QTL mapping for two melanic pigmentation traits (melanophore density and lateral barring) using benthic-limnetic F2 crosses. The two traits mapped to different chromosomes, suggesting a distinct genetic basis. The resulting QTLs were additive, but explained a relatively small fraction of the total variance (~6%). QTLs maps differed by F1 family, suggesting variation in genetic architecture or ability to detect loci of small effect. Functional analysis identified enriched pathways for candidate loci. Several of the resulting candidate loci for pigmentation, including three loci in enriched pathways (bco1, sulf1, and tyms) have been previously indicated to affect pigmentation in other vertebrates. These findings add to a growing body of evidence suggesting pigmentation is often polygenic.

Effects of urban-induced mutations on ecology, evolution and health.

Increasing evidence suggests that urbanization is associated with higher mutation rates, which can affect the health and evolution of organisms that inhabit cities. Elevated pollution levels in urban areas can induce DNA damage, leading to de novo mutations. Studies on mutations induced by urban pollution are most prevalent in humans and microorganisms, whereas studies of non-human eukaryotes are rare, even though increased mutation rates have the potential to affect organisms and their populations in contemporary time. Our Perspective explores how higher mutation rates in urban environments could impact the fitness, ecology and evolution of populations. Most mutations will be neutral or deleterious, and higher mutation rates associated with elevated pollution in urban populations can increase the risk of cancer in humans and potentially other species. We highlight the potential for urban-driven increased deleterious mutational loads in some organisms, which could lead to a decline in population growth of a wide diversity of organisms. Although beneficial mutations are expected to be rare, we argue that higher mutation rates in urban areas could influence adaptive evolution, especially in organisms with short generation times. Finally, we explore avenues for future research to better understand the effects of urban-induced mutations on the fitness, ecology and evolution of city-dwelling organisms.

Gut Microbiota Uniqueness Is Associated with Lake Size, a Proxy for Diet Diversity, in Stickleback Fish.

AbstractOrganismal divergence can be driven by differential resource use and adaptation to different trophic niches. Variation in diet is a major factor shaping the gut microbiota, which is crucial for many aspects of their hosts' biology. However, it remains largely unknown how host diet diversity affects the gut microbiota, and it could be hypothesized that trophic niche width is positively associated with gut microbiota diversity. To test this idea, we sequenced the 16S ribosomal RNA gene from intestinal tissue of 14 threespine stickleback populations from lakes of varying size on Vancouver Island, Canada, that have been shown to differ in trophic niche width. Using lake size as a proxy for trophic ecology, we found evidence for higher gut microbiota uniqueness among individuals from populations with broader trophic niches. While these results suggest that diet diversity might promote gut microbiota diversity, additional work investigating diet and gut microbiota variation of the same host organisms will be necessary. Yet our results motivate the question of how host population diversity (e.g., ecological, morphological, genetic) might interact with the gut microbiota during the adaptation to ecological niches.

The effects of host ecology and phylogeny on gut microbiota (non)parallelism across birds and mammals.

IMPORTANCE: What are the roles of determinism and contingency in evolution? The paleontologist and evolutionary biologist Stephen J. Gould raised this question in his famous thought experiment of "replaying life's tape." Settings where independent lineages have repeatedly adapted to similar ecological niches (i.e., parallel evolution) are well suited to address this question. Here, we quantified whether repeated ecological shifts across 53 mammalian and 50 avian host species are associated with parallel gut microbiota changes. Our results indicate that parallel shifts in host diet are associated with greater gut microbiota parallelism (i.e., more deterministic). While further research will be necessary to obtain a comprehensive picture of the circumstances under which deterministic gut microbiota changes might be expected, our study can be instrumental in motivating the use of more quantitative methods in microbiota research. This, in turn, can help us better understand microbiota dynamics during adaptive evolution of their hosts.

Gut microbiota parallelism and divergence associated with colonisation of novel habitats.

An organism's gut microbiota can change in response to novel environmental conditions, in particular when colonisation of new habitats is accompanied by shifts in the host species' ecology. Here, we investigated the gut microbiota of three lizard species (A. inornata, H. maculata and S. cowlesi) from their ancestral-like habitat in the Chihuahuan desert and two colonised habitats with contrasting geological and ecological compositions: the White Sands and Carrizozo lava flow. The host species and the lizards' environment both shape gut microbiota composition, but host effects were overall stronger. Further, we found evidence that colonisation of the same environment by independent host species led to parallel changes of the gut microbiota, whereas the colonisation of two distinct environments by the same host species led to gut microbiota divergence. Some of the gut microbiota changes that accompanied the colonisation of the White Sands were associated with shifts in diet (based on diet information from previous studies), which is congruent with the general observation that trophic ecology has a strong effect on gut microbiota composition. Our study provides insights into how shifts in host ecology accompanying colonisation of novel environments can affect gut microbiota composition and diversity.

Assessing the validity of fecal sampling for characterizing variation in threespine stickleback's gut microbiota.

The gut microbiota is crucial for many aspects of their hosts' biology, and it has been characterized for many species across the animal kingdom. Yet, we still don't have a good understanding of whether non-lethal sampling can accurately capture the diversity of gut-associated bacterial communities, as estimated from lethal sampling of intestinal tissue. We further lack knowledge on whether non-lethal sampling methods are suitable for detecting gut microbiota shifts associated with changes in environmental factors (e.g., diet). We addressed these questions in threespine stickleback fish, a model system for evolutionary ecology, by comparing bacterial communities from intestinal tissue and feces. Despite some differences in community composition between the two sample types and considerable temporal variation among fecal samples, bacterial communities appear to largely overlap. Further, we detected consistent and significant changes of fecal bacterial communities associated with an experimental diet manipulation. This suggests that fecal sampling can represent an adequate non-lethal method to characterize the gut microbiota of threespine stickleback, but additional studies will be necessary before drawing general conclusions regarding the validity of fecal sampling for gut microbiota studies. To this end, we give recommendations to improve the characterization of the gut microbiota via fecal sampling. Fecal sampling allows studying temporal gut microbiota shifts associated with environmental change at the individual level, which increases opportunities for future experimental gut microbiota research.

Repeated genetic divergence plays a minor role in repeated phenotypic divergence of lake-stream stickleback.

Recent studies have shown that the repeated evolution of similar phenotypes in response to similar ecological conditions (here "parallel evolution") often occurs through mutations in the same genes. However, many previous studies have focused on known candidate genes in a limited number of systems. Thus, the question of how often parallel phenotypic evolution is due to parallel genetic changes remains open. Here, we used quantitative trait locus (QTL) mapping in F2 intercrosses between lake and stream threespine stickleback (Gasterosteus aculeatus) from four independent watersheds on Vancouver Island, Canada to determine whether the same QTL underlie divergence in the same phenotypes across, between, and within watersheds. We find few parallel QTL, even in independent crosses from the same watershed or for phenotypes that have diverged in parallel. These findings suggest that different mutations can lead to similar phenotypes. The low genetic repeatability observed in these lake-stream systems contrasts with the higher genetic repeatability observed in other stickleback systems. We speculate that differences in evolutionary history, gene flow, and/or the strength and direction of selection might explain these differences in genetic parallelism and emphasize that more work is needed to move beyond documenting genetic parallelism to identifying the underlying causes.

Pleiotropy facilitates parallel adaptation in sticklebacks.

Highly pleiotropic genes are predicted to be used less often during adaptation, as mutations in these loci are more likely to have negative fitness consequences. Following this logic, we tested whether pleiotropy impacts the probability that a locus will be used repeatedly in adaptation. We used two proxies to estimate pleiotropy: number of phenotypic traits affected by a given genomic region and gene connectivity. We first surveyed 16 independent stream-lake and three independent benthic-limnetic ecotype pairs of threespine stickleback to estimate genome-wide patterns in parallel genomic differentiation. Our analysis revealed parallel divergence across the genome; 30%-37% of outlier regions were shared between at least two independent pairs in either the stream-lake or benthic-limnetic comparisons. We then tested whether parallel genomic regions are less pleiotropic than nonparallel regions. Counter to our a priori prediction, parallel genomic regions contained genes with significantly more pleiotropy; that is, influencing a greater number of traits and more highly connected. The increased pleiotropy of parallel regions could not be explained by other genomic factors, as there was no significant difference in mean gene count, mutation or recombination rates between parallel and nonparallel regions. Interestingly, although nonparallel regions contained genes that were less connected and influenced fewer mapped traits on average than parallel regions, they also tended to contain the genes that were predicted to be the most pleiotropic. Taken together, our findings are consistent with the idea that pleiotropy only becomes constraining at high levels and that low or intermediate levels of pleiotropy may be beneficial for adaptation.

Analysis of ancestry heterozygosity suggests that hybrid incompatibilities in threespine stickleback are environment dependent.

Hybrid incompatibilities occur when interactions between opposite ancestry alleles at different loci reduce the fitness of hybrids. Most work on incompatibilities has focused on those that are "intrinsic," meaning they affect viability and sterility in the laboratory. Theory predicts that ecological selection can also underlie hybrid incompatibilities, but tests of this hypothesis using sequence data are scarce. In this article, we compiled genetic data for F2 hybrid crosses between divergent populations of threespine stickleback fish (Gasterosteus aculeatus L.) that were born and raised in either the field (seminatural experimental ponds) or the laboratory (aquaria). Because selection against incompatibilities results in elevated ancestry heterozygosity, we tested the prediction that ancestry heterozygosity will be higher in pond-raised fish compared to those raised in aquaria. We found that ancestry heterozygosity was elevated by approximately 3% in crosses raised in ponds compared to those raised in aquaria. Additional analyses support a phenotypic basis for incompatibility and suggest that environment-specific single-locus heterozygote advantage is not the cause of selection on ancestry heterozygosity. Our study provides evidence that, in stickleback, a coarse-albeit indirect-signal of environment-dependent hybrid incompatibility is reliably detectable and suggests that extrinsic incompatibilities can evolve before intrinsic incompatibilities.

Parallel shifts of visual sensitivity and body coloration in replicate populations of extremophile fish.

Visual sensitivity and body pigmentation are often shaped by both natural selection from the environment and sexual selection from mate choice. One way of quantifying the impact of the environment is by measuring how traits have changed after colonization of a novel habitat. To do this, we studied Poecilia mexicana populations that have repeatedly adapted to extreme sulphidic (H2 S-containing) environments. We measured visual sensitivity using opsin gene expression, as well as body pigmentation, for populations in four independent drainages. Both visual sensitivity and body pigmentation showed significant parallel shifts towards greater medium-wavelength sensitivity and reflectance in sulphidic populations. Altogether we found that sulphidic habitats select for differences in visual sensitivity and pigmentation. Shifts between habitats may be due to both differences in the water's spectral properties and correlated ecological changes.

Quantifying (non)parallelism of gut microbial community change using multivariate vector analysis.

Parallel evolution of phenotypic traits is regarded as strong evidence for natural selection and has been studied extensively in a variety of taxa. However, we have limited knowledge of whether parallel evolution of host organisms is accompanied by parallel changes of their associated microbial communities (i.e., microbiotas), which are crucial for their hosts' ecology and evolution. Determining the extent of microbiota parallelism in nature can improve our ability to identify the factors that are associated with (putatively adaptive) shifts in microbial communities. While it has been emphasized that (non)parallel evolution is better considered as a quantitative continuum rather than a binary phenomenon, quantitative approaches have rarely been used to study microbiota parallelism. We advocate using multivariate vector analysis (i.e., phenotypic change vector analysis) to quantify direction and magnitude of microbiota changes and discuss the applicability of this approach for studying parallelism, and we compiled an R package for multivariate vector analysis of microbial communities ('multivarvector'). We exemplify its use by reanalyzing gut microbiota data from multiple fish species that exhibit parallel shifts in trophic ecology. We found that multivariate vector analysis results were largely consistent with other statistical methods, parallelism estimates were not affected by the taxonomic level at which the microbiota is studied, and parallelism might be stronger for gut microbiota function compared to taxonomic composition. This approach provides an analytical framework for quantitative comparisons across host lineages, thereby providing the potential to advance our capacity to predict microbiota changes. Hence, we emphasize that the development and application of quantitative measures, such as multivariate vector analysis, should be further explored in microbiota research in order to better understand the role of microbiota dynamics during their hosts' adaptive evolution, particularly in settings of parallel evolution.

The genomic signature of ecological divergence along the benthic-limnetic axis in allopatric and sympatric threespine stickleback.

The repeated occurrence of similar phenotypes in independent lineages (i.e., parallel evolution) in response to similar ecological conditions can provide compelling insights into the process of adaptive evolution. An intriguing question is to what extent repeated phenotypic changes are underlain by repeated changes at the genomic level and whether patterns of genomic divergence differ with the geographic context in which populations evolve. Here, we combined genomic, morphological and ecological data sets to investigate the genomic signatures of divergence across populations of threespine stickleback (Gasterosteus aculeatus) that adapted to contrasting ecological niches (benthic or limnetic) in either sympatry or allopatry. We found that genome-wide differentiation (FST ) was an order of magnitude higher and substantially more repeatable for sympatric benthic and limnetic specialists compared to allopatric populations with similar levels of ecological divergence. We identified genomic regions consistently differentiated between sympatric ecotypes that were also differentiated between or associated with benthic vs. limnetic niche in allopatric populations. These candidate regions were enriched on three chromosomes known to be involved in the benthic-limnetic divergence of threespine stickleback. Some candidate regions overlapped with QTL for body shape and trophic traits such as gill raker number, traits that strongly differ between benthic and limnetic ecotypes. In summary, our study shows that magnitude and repeatability of genomic signatures of ecological divergence in threespine stickleback highly depend on the geographic context. The identified candidate regions provide starting points to identify functionally important genes for the adaptation to benthic and limnetic niches.

Shared Patterns of Genome-Wide Differentiation Are More Strongly Predicted by Geography Than by Ecology.

Closely related populations often display similar patterns of genomic differentiation, yet it remains an open question which ecological and evolutionary forces generate these patterns. The leading hypothesis is that this similarity in divergence is driven by parallel natural selection. However, several recent studies have suggested that these patterns may instead be a product of the depletion of genetic variation that occurs as result of background selection (i.e., linked negative selection). To date, there have been few direct tests of these competing hypotheses. To determine the relative contributions of background selection and parallel selection to patterns of repeated differentiation, we examined 24 independently derived populations of freshwater stickleback occupying a variety of niches and estimated genomic patterns of differentiation in each relative to their common marine ancestor. Patterns of genetic differentiation were strongly correlated across pairs of freshwater populations adapting to the same ecological niche, supporting a role for parallel natural selection. In contrast to other recent work, our study comparing populations adapting to the same niche produced no evidence signifying that similar patterns of genomic differentiation are generated by background selection. We also found that overall patterns of genetic differentiation were considerably more similar for populations found in closer geographic proximity. In fact, the effect of geography on the repeatability of differentiation was greater than that of parallel selection. Our results suggest that shared selective landscapes and ancestral variation are the key drivers of repeated patterns of differentiation in systems that have recently colonized novel environments.

Parallel changes in gut microbiome composition and function during colonization, local adaptation and ecological speciation.

The processes of local adaptation and ecological speciation are often strongly shaped by biotic interactions such as competition and predation. One of the strongest lines of evidence that biotic interactions drive evolution comes from the repeated divergence of lineages in association with repeated changes in the community of interacting species. Yet relatively little is known about the repeatability of changes in gut microbial communities and their role in adaptation and divergence of host populations in nature. Here we use three cases of rapid, parallel adaptation and speciation in freshwater threespine stickleback to test for parallel changes in associated gut microbiomes. We find that features of the gut microbial communities have shifted repeatedly in the same direction in association with parallel divergence and speciation of stickleback hosts. These results suggest that changes to gut microbiomes can occur rapidly and predictably in conjunction with host evolution, and that host-microbe interactions might play an important role in host adaptation and diversification.

Genetics of adaptation: Experimental test of a biotic mechanism driving divergence in traits and genes.

The genes underlying adaptations are becoming known, yet the causes of selection on genes-a key step in the study of the genetics of adaptation-remains uncertain. We address this issue experimentally in a threespine stickleback species pair showing exaggerated divergence in bony defensive armor in association with competition-driven character displacement. We used semi-natural ponds to test the role of a native predator in causing divergent evolution of armor and two known underlying genes. Predator presence/absence altered selection on dorsal spines and allele frequencies at the Msx2a gene across a generation. Evolutionary trajectories of alleles at a second gene, Pitx1, and the pelvic spine trait it controls, were more variable. Our experiment demonstrates how manipulation of putative selective agents helps to identify causes of evolutionary divergence at key genes, rule out phenotypic plasticity as a sole determinant of phenotypic differences, and eliminate reliance on fitness surrogates. Divergence of predation regimes in sympatric stickleback is associated with coevolution in response to resource competition, implying a cascade of biotic interactions driving species divergence. We suggest that as divergence proceeds, an increasing number of biotic interactions generate divergent selection, causing more evolution in turn. In this way, biotic adaptation perpetuates species divergence through time during adaptive radiation in an expanding number of traits and genes.

Ecological factors and morphological traits are associated with repeated genomic differentiation between lake and stream stickleback.

The repeated evolution of similar phenotypes in independent populations (i.e. parallel or convergent evolution) provides an opportunity to identify genetic and ecological factors that influence the process of adaptation. Threespine stickleback fish ( Gasterosteus aculeatus) are an excellent model for such studies, as they have repeatedly adapted to divergent habitats across the Northern hemisphere. Here, we use genomic, ecological and morphological data from 16 independent pairs of stickleback populations adapted to divergent lake and stream habitats. We combine a population genomic approach to identify regions of the genome that are likely under selection in these divergent habitats with an association mapping approach to identify regions of the genome that underlie variation in ecological factors and morphological traits. Over 37% of genomic windows are repeatedly differentiated across lake-stream pairs. Similarly, many genomic windows are associated with variation in abiotic factors, diet items and morphological phenotypes. Both the highly differentiated windows and candidate trait windows are non-randomly distributed across the genome and show some overlap. However, the overlap is not significant on a genome-wide scale. Together, our data suggest that adaptation to divergent food resources and predation regimes are drivers of differentiation in lake-stream stickleback, but that additional ecological factors are also important. This article is part of the theme issue 'Convergent evolution in the genomics era: new insights and directions'.

Early-life manipulation of cortisol and its receptor alters stress axis programming and social competence.

In many vertebrate species, early social experience generates long-term effects on later life social behaviour. These effects are accompanied by persistent modifications in the expression of genes implicated in the stress axis. It is unknown, however, whether stress axis programming can affect the development of social competence, and if so, by which mechanism(s). Here, we used pharmacological manipulations to persistently reprogramme the hypothalamic-pituitary-interrenal axis of juvenile cooperatively breeding cichlids, Neolamprologus pulcher. During the first two months of life, juveniles were repeatedly treated with cortisol, mifepristone or control treatments. Three months after the last manipulation, we tested for treatment effects on (i) social competence, (ii) the expression of genes coding for corticotropin-releasing factor ( crf), glucocorticoid receptor ( gr1) and mineralocorticoid receptor ( mr) in the telencephalon and hypothalamus and (iii) cortisol levels. Social competence in a social challenge was reduced in cortisol-treated juveniles, which is in accordance with previous work applying early-life manipulations using different social experiences. During early life, both cortisol and mifepristone treatments induced a persistent downregulation of crf and upregulation of mr in the telencephalon. We suggest that these persistent changes in stress gene expression may represent an effective physiological mechanism for coping with stress. This article is part of the theme issue 'Developing differences: early-life effects and evolutionary medicine'.

Differential predation alters pigmentation in threespine stickleback (Gasterosteus aculeatus).

Animal pigmentation plays a key role in many biological interactions, including courtship and predator avoidance. Sympatric benthic and limnetic ecotypes of threespine stickleback (Gasterosteus aculeatus) exhibit divergent pigment patterns. To test whether differential predation by cutthroat trout contributes to the differences in pigmentation seen between the ecotypes, we used a within-generation selection experiment on F2 benthic-limnetic hybrids. After 10 months of differential selection, we compared the pigmentation of fish under trout predation to control fish not exposed to trout predation. We found that stickleback exhibited more lateral barring in ponds with trout predation. Ponds with trout were also less turbid, and a greater degree of barring was negatively correlated with the magnitude of turbidity across pond replicates. A more benthic diet, a proxy for habitat use, was also correlated with greater lateral barring and green dorsal pigmentation. These patterns suggest that differential exposure to cutthroat trout predation may explain the divergence in body pigmentation between benthic and limnetic ecotypes.

Gene flow and selection interact to promote adaptive divergence in regions of low recombination.

Adaptation to new environments often occurs in the face of gene flow. Under these conditions, gene flow and recombination can impede adaptation by breaking down linkage disequilibrium between locally adapted alleles. Theory predicts that this decay can be halted or slowed if adaptive alleles are tightly linked in regions of low recombination, potentially favouring divergence and adaptive evolution in these regions over others. Here, we compiled a global genomic data set of over 1,300 individual threespine stickleback from 52 populations and compared the tendency for adaptive alleles to occur in regions of low recombination between populations that diverged with or without gene flow. In support of theory, we found that putatively adaptive alleles (FST and dXY outliers) tend to occur more often in regions of low recombination in populations where divergent selection and gene flow have jointly occurred. This result remained significant when we employed different genomic window sizes, controlled for the effects of mutation rate and gene density, controlled for overall genetic differentiation, varied the genetic map used to estimate recombination and used a continuous (rather than discrete) measure of geographic distance as proxy for gene flow/shared ancestry. We argue that our study provides the first statistical evidence that the interaction of gene flow and selection biases divergence toward regions of low recombination.

Evolutionary ecology of opsin gene sequence, expression and repertoire.

Linking molecular evolution to biological function is a long-standing challenge in evolutionary biology. Some of the best examples of this involve opsins, the genes that encode the molecular basis of light reception. In this issue of Molecular Ecology, three studies examine opsin gene sequence, expression and repertoire to determine how natural selection has shaped the visual system. First, Escobar-Camacho et al. () use opsin repertoire and expression in three Amazonian cichlid species to show that a shift in sensitivity towards longer wavelengths is coincident with the long-wavelength-dominated Amazon basin. Second, Stieb et al. () explore opsin sequence and expression in reef-dwelling damselfish and find that UV- and long-wavelength vision are both important, but likely for different ecological functions. Lastly, Suvorov et al. () study an expansive opsin repertoire in the insect order Odonata and find evidence that copy number expansion is consistent with the permanent heterozygote model of gene duplication. Together these studies emphasize the utility of opsin genes for studying both the local adaptation of sensory systems and, more generally, gene family evolution.