Molecular Mechanisms of Sexual Dimorphism in Threespine Stickleback

Sexual dimorphism is commonly understood as differences in external features, such as morphological features or coloration. However, it can more broadly encompass behavior and physiology and at the core of these differences is the genetic mechanism – mRNA and protein expression. How, and which, molecular mechanisms influence sexually dimorphic features is not well understood thus far. DNA, RNA and proteins are the template required to create the phenotype of an individual, and they are connected to each other via processes of transcription and translation. As the genome of males and females are almost identical with the exception of the few genes on the sex chromosome or the sex-determining alleles (in the case of organisms without sex chromosomes), it is likely that many of the downstream processes resulting in sexual dimorphism are produced by changes in gene regulation and result from a regulatory cascade and not from a vastly different gene composition. Thus, in this thesis a systems biology approach is used to understand sexual dimorphism at all molecular levels and how different genomic features, e.g. sex chromosome evolution, can affect the interplay of these molecules. The threespine stickleback, Gasterosteus aculeatus, is used as the model to investigate molecular mechanisms of sexual dimorphism. It has well-characterized ecology and behavior, especially in the breeding season when sexual dimorphism is high. Moreover, threespine stickleback has a recently evolved Y chromosome in the early stages of sex chromosome evolution, characterized by a lack of recombination leading to degeneration (i.e. gene loss). The aim of my thesis is to investigate how the genotype links to the molecular phenotype and relates to differences in molecular expression between males and females. Based on previous research on sex differences in mRNA expression, I investigated sex-biased protein expression in adult fish outside the breeding season to see if differences persisted after translation. As sex-biased expression also prevailed in the proteome and previous transcription expression seemed to be related to the sex chromosomes, I investigated the genome level with a particular focus on the sex-chromosomes. I characterized the status of Y chromosome degeneration in the threespine stickleback and its effects on gene function. Furthermore, since the degeneration process leaves genes in a single copy in males, I examined whether the resulting dosage difference of messenger RNA for hemizygous genes is compensated as it is in other organisms. In addition, threespine sticklebacks have wellcharacterized behavioral differences related to the male’s social status during the breeding season. To understand the connection between the genotype and behavior, I examined gene expression patterns related to breeding behavior using dominant and subordinate males as well as female

Partitioning the relative fitness effects of diet and trophic morphology in the threespine stickleback

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

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.

Distinctive insular forms of threespine stickleback (Gasterosteus aculeatus) from western Mediterranean islands

Neutral and adaptive variation among populations within a species is a major component of biological diversity and may be pronounced among insular populations due to geographical isolation and island specific evolutionary forces at work. Detecting and preserving potential evolutionary significant units below the species rank has become a crucial task for conservation biology. Combining genetic, phenotypic and ecological data, we investigated evolutionary patterns among the enigmatic threespine stickleback populations from western Mediterranean islands, all of which are threatened by habitat deterioration and climate change. We find indications that these populations derive from different genetic lineages, being genetically highly distinct from the stickleback of mainland Europe and the northern Atlantic as well as from each other. Mediterranean island stickleback populations are also phenotypically distinct from mainland populations but interestingly stickleback from Iceland have converged on a similar phenotype. This distinctive island stickleback phenotype seems to be driven by distinct selective regimes on islands versus continents. Overall, our results reveal the status of western Mediterranean island stickleback as evolutionarily distinct units, important for conservation of biodiversity.

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.

Cryptic invasion drives phenotypic changes in central European threespine stickleback

Cryptic invasions are commonly associated with genetic changes of the native species or genetic lineage that the invaders replace. Phenotypic shifts resulting from cryptic invasions are less commonly reported given the relative paucity of historical specimens that document such phenotypic changes. Here, I study such a case in two populations of threespine stickleback from central Europe, comparing contemporary patterns of gene flow with phenotypic changes between historical and contemporary population samples. I find gene flow from an invasive lineage to be associated with significant phenotypic changes, where the degree of phenotypic change corresponds with the level of gene flow that a population receives. These findings underline the utility of combining genetic approaches with phenotypic data to estimate the impact of gene flow in systems where anthropogenic alterations have removed former geographic barriers promoting cryptic invasions.

Centromere identification and inactivation on a neo-Y chromosome fusion in threespine stickleback fish (Gasterosteus aculeatus)

Thesis (Ph.D.)--University of Washington, 2016-05

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.

Habitat Choice and Female Preference in a Polymorphic Stickleback Population

Background: A small pond, c. 90 years old, near Bern, Switzerland contains a population of threespine stickleback (Gasterosteus aculeatus) with two distinct male phenotypes. Males of one type are large, and red, and nest in the shallow littoral zone. The males of the other are small and orange, and nest offshore at slightly greater depth. The females in this population are phenotypically highly variable but cannot easily be assigned to either male type. Question: Is the existence of two sympatric male morphs maintained by substrate-associated male nest site choice and facilitated by female mate preferences? Organisms: Male stickleback caught individually at their breeding sites. Females caught with minnow traps. Methods: In experimental tanks, we simulated the slope and substrate of the two nesting habitats. We then placed individual males in a tank and observed in which habitat the male would build his nest. In a simultaneous two-stimulus choice design, we gave females the choice between a large, red male and a small, orange one. We measured female morphology and used linear mixed effect models to determine whether female preference correlated with female morphology. Results: Both red and orange males preferred nesting in the habitat that simulated the slightly deeper offshore condition. This is the habitat occupied by the small, orange males in the pond itself. The proportion of females that chose a small orange male was similar to that which chose a large red male. Several aspects of female phenotype correlated with the male type that a female preferred.

Ecosystem size matters: The dimensionality of intralacustrine diversification in Icelandic stickleback is predicted by lake size

Cases of evolutionary diversification can be characterized along a continuum from weak to strong genetic and phenotypic differentiation. Several factors may facilitate or constrain the differentiation process. Comparative analyses of replicates of the same taxon at different stages of differentiation can be useful to identify these factors. We estimated the number of distinct phenotypic groups in threespine stickleback populations from nine lakes in Iceland and in one marine population. Using the inferred number of phenotypic groups in each lake, genetic divergence from the marine population, and physical lake and landscape variables, we tested if ecosystem size, approximated by lake size and depth, or isolation from the ancestral marine gene pool predict the occurrence and the extent of phenotypic and genetic diversification within lakes. We find intralacustrine phenotypic diversification to be the rule rather than the exception, occurring in all but the youngest lake population and being manifest in ecologically important phenotypic traits. Neutral genetic data further indicates non-random mating in four out of nine studied lakes, and restricted gene flow between sympatric phenotypic groups in two. Although neither the phenotypic variation nor the number of intralacustrine phenotypic groups were associated with any of our environmental variables, the number of phenotypic traits that were differentiated was significantly positively related to lake size, and evidence for restricted gene flow between sympatric phenotypic groups was only found in the largest lakes where trait specific phenotypic differentiation was highest.

Evidence of Adaptive Evolutionary Divergence during Biological Invasion

Rapid phenotypic diversification during biological invasions can either arise by adaptation to alternative environments or by adaptive phenotypic plasticity. Where experimental evidence for adaptive plasticity is common, support for evolutionary diversification is rare. Here, we performed a controlled laboratory experiment using full-sib crosses between ecologically divergent threespine stickleback populations to test for a genetic basis of adaptation. Our populations are from two very different habitats, lake and stream, of a recently invaded range in Switzerland and differ in ecologically relevant morphological traits. We found that in a lake-like food treatment lake fish grow faster than stream fish, resembling the difference among wild type individuals. In contrast, in a stream-like food treatment individuals from both populations grow similarly. Our experimental data suggest that genetically determined diversification has occurred within less than 140 years after the arrival of stickleback in our studied region.

Disentangling the role of phenotypic plasticity and genetic divergence in contemporary ecotype formation during a biological invasion

The occurrence of contemporary ecotype formation through adaptive divergence of populations within the range of an invasive species typically requires standing genetic variation but can be facilitated by phenotypic plasticity. The relative contributions of both of these to adaptive trait differentiation have rarely been simultaneously quantified in recently diverging vertebrate populations. Here we study a case of intraspecific divergence into distinct lake and stream ecotypes of threespine stickleback that evolved in the past 140 years within the invasive range in Switzerland. Using a controlled laboratory experiment with full-sib crosses and treatments mimicking a key feature of ecotypic niche divergence, we test if the phenotypic divergence that we observe in the wild results from phenotypic plasticity or divergent genetic predisposition. Our experimental groups show qualitatively similar phenotypic divergence as those observed among wild adults. The relative contribution of plasticity and divergent genetic predisposition differs among the traits studied, with traits related to the biomechanics of feeding showing a stronger genetic predisposition, whereas traits related to locomotion are mainly plastic. These results implicate that phenotypic plasticity and standing genetic variation interacted during contemporary ecotype formation in this case.

Hybrid ‘superswarm’ leads to rapid divergence and establishment of populations during a biological invasion

Understanding the genetic background of invading species can be crucial information clarifying why they become invasive. Intraspecific genetic admixture among lineages separated in the native ranges may promote the rate and extent of an invasion by substantially increasing standing genetic variation. Here we examine the genetic relationships among threespine stickleback that recently colonized Switzerland. This invasion results from several distinct genetic lineages that colonized multiple locations and have since undergone range expansions, where they coexist and admix in parts of their range. Using 17 microsatellites genotyped for 634 individuals collected from 17 Swiss and two non-Swiss European sites, we reconstruct the invasion of stickleback and investigate the potential and extent of admixture and hybridization among the colonizing lineages from a population genetic perspective. Specifically we test for an increase in standing genetic variation in populations where multiple lineages coexist. We find strong evidence of massive hybridization early on, followed by what appears to be recent increased genetic isolation and the formation of several new genetically distinguishable populations, consistent with a hybrid ‘superswarm’. This massive hybridization and population formation event(s) occurred over approximately 140 years and likely fuelled the successful invasion of a diverse range of habitats. The implications are that multiple colonizations coupled with hybridization can lead to the formation of new stable genetic populations potentially kick-starting speciation and adaptive radiation over a very short time.

Phenotypic divergence along lines of genetic variance

Natural populations inhabiting the same environment often independently evolve the same phenotype. Is this replicated evolution a result of genetic constraints imposed by patterns of genetic covariation? We looked for associations between directions of morphological divergence and the orientation of the genetic variance-covariance matrix (G) by using an experimental system of morphological evolution in two allopatric nonsister species of rainbow fish. Replicate populations of both Melanotaenia eachamensis and Melanotaenia duboulayi have independently adapted to lake versus stream hydrodynamic environments. The major axis of divergence (z) among all eight study populations was closely associated with the direction of greatest genetic variance (g(max)), suggesting directional genetic constraint on evolution. However, the direction of hydrodynamic adaptation was strongly associated with vectors of G describing relatively small proportions of the total genetic variance, and was only weakly associated with g(max). In contrast, divergence between replicate populations within each habitat was approximately proportional to the level of genetic variance, a result consistent with theoretical predictions for neutral phenotypic divergence. Divergence between the two species was also primarily along major eigenvectors of G. Our results therefore suggest that hydrodynamic adaptation in rainbow fish was not directionally constrained by the dominant eigenvector of G. Without partitioning divergence as a consequence of the adaptation of interest (here, hydrodynamic adaptation) from divergence due to other processes, empirical studies are likely to overestimate the potential for the major eigenvectors of G to directionally constrain adaptive evolution.