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.

Behavioural responses of threespine stickleback with lateral line asymmetries to experimental mechanosensory stimuli

Behavioural asymmetry, typically referred to as laterality, is widespread among bilaterians and is often associated with asymmetry in brain structure. However, the influence of sensory receptor asymmetry on laterality has undergone limited investigation. Here we use threespine stickleback (Gasterosteus aculeatus) to investigate the influence of lateral line asymmetry on laterality during lab simulations of three mechanosensation-dependent behaviours: predator evasion, prey localization and rheotaxis. We recorded the response of stickleback to impacts at the water surface and water flow in photic conditions and low-frequency oscillations in the dark, across four repeat trials. We then compared individuals’ laterality to asymmetry in the number of neuromasts on either side of their body. Stickleback hovered with their right side against the arena wall 57% of the time (P<0.001) in illuminated surface impact trials and 56% of the time in (P=0.085) dark low-frequency stimulation trials. Light regime modulated the effect of neuromast count on laterality, as fish with more neuromasts were more likely to hover with the wall on their right during illumination (P=0.007) but were less likely to do so in darkness (P=0.025). Population level laterality diminished in later trials across multiple behaviours and individuals did not show a consistent side bias in any behaviours. Our results demonstrate a complex relationship between sensory structure asymmetry and laterality, suggesting that laterality is modulated multiple sensory modalities and temporally dynamic.

Single-cell RNA sequencing reveals micro-evolution of the stickleback immune system

Pathogenic infection is an important driver of many ecological processes. Furthermore, variability in immune function is an important driver of differential infection outcomes. New evidence would suggest that immune variation extends to broad cellular structure of immune systems. However, variability at such broad levels is traditionally difficult to detect in non-model systems. Here we leverage single cell transcriptomic approaches to document signatures of microevolution of immune system structure in a natural system, the three-spined stickleback (Gasterosteus aculeatus). We sampled nine adult fish from three populations with variability in resistance to a cestode parasite, Schistocephalus solidus, to create the first comprehensive immune cell atlas for G. aculeatus. Eight major immune cell types, corresponding to major vertebrate immune cells, were identified. We were also able to document significant variation in both abundance and expression profiles of the individual immune cell types, among the three populations of fish. This variability may contribute to observed variability in parasite susceptibility. Finally, we demonstrate that identified cell type markers can be used to reinterpret traditional transcriptomic data. Combined our study demonstrates the power of single cell sequencing to not only document evolutionary phenomena (i.e. microevolution of immune cells), but also increase the power of traditional transcriptomic datasets.

Insights Into Amino Acid Metabolism and Incorporation by Compound- Specific Carbon Isotope Analysis of three-Spined Sticklebacks

Interpretation of isotope data is of upmost importance in ecology to build sound models for the study of animal diets, migration patterns and physiology. However, our understanding of isotope fractionation and incorporation into different consumer tissues is still limited. We therefore measured the δ13C values of individual amino acids from muscle and liver tissue of three-spined sticklebacks (Gasterosteus aculeatus) over time in a controlled feeding experiment. The isotope signatures of amino acids in the liver quickly responded to even small shifts in dietary isotope compositions, whereas muscle tissue remained isotopically constant over time. No significant fractionation between diet and fish tissues was observed for half of the amino acids in both liver and muscle tissue, supporting the idea of direct nutrient routing on a protein rich diet. Small isotopic differences were observed for alanine, glycine and serine between diet and liver, indicating that metabolic processes such as glycolysis or gluconeogenesis can be tracked by the isotope signatures of their main substitutes. We demonstrate that compound-specific isotope analysis has great potential to investigate the main metabolic pathways of organisms and suggest further investigations using isotopically enriched materials to facilitate the correct interpretation of ecological field data.

Evolution of stickleback spines through independent cis-regulatory changes at HOXDB

Understanding the genetic mechanisms leading to new traits is a fundamental goal of evolutionary biology. We show that HOXDB regulatory changes have been used repeatedly in different stickleback fish species to alter the length and number of bony dorsal spines. In Gasterosteus aculeatus, a variant HOXDB allele is genetically linked to shortening an existing spine and adding a spine. In Apeltes quadracus, a variant allele is associated with lengthening an existing spine and adding a spine. The alleles alter the same conserved non-coding HOXDB enhancer by diverse molecular mechanisms, including SNPs, deletions, and transposable element insertions. The independent cis-acting regulatory changes are linked to anterior expansion or contraction of HOXDB expression. Our findings support the long-standing hypothesis that natural Hox gene variation underlies key morphological patterning changes in wild populations and illustrate how different mutational mechanisms affecting the same region may produce opposite gene expression changes with similar phenotypic outcomes.

Isolating and growing fibroblast cells from threespine stickleback (Gasterosteus aculeatus) v2

This protocol goes through the process of isolating and growing fibroblast cells from threespine stickleback.

Shedding Light on the Circadian Clock of the Threespine Stickleback

The circadian clock is an internal timekeeping system shared by most organisms, and knowledge about its functional importance and evolution in natural environments is still needed. Here, we investigated the circadian clock of wild-caught threespine sticklebacks (Gasterosteus aculeatus) at the behavioural and molecular levels. While their behaviour, ecology, and evolution are well studied, information on their circadian rhythms are scarce. We quantified the daily locomotor activity rhythm under a light-dark cycle (LD) and under constant darkness (DD). Under LD, all fish exhibited significant daily rhythmicity, while under DD, only 18% of individuals remained rhythmic. This interindividual variation suggests that the circadian clock controls activity only in certain individuals. Moreover, under LD, some fish were almost exclusively nocturnal, while others were active around the clock. Furthermore, the most nocturnal fish were also the least active. These results suggest that light masks activity (i.e. suppresses activity without entraining the internal clock) more strongly in some individuals than others. Finally, we quantified the expression of five clock genes in the brain of sticklebacks under DD using qPCR. We did not detect circadian rhythmicity, which could either indicate that the clock molecular oscillator is highly light-dependent, or that there was an oscillation but that we were unable to detect it. Overall, our study suggests that a strong circadian control on behavioural rhythms may not necessarily be advantageous in a natural population of sticklebacks and that the daily phase of activity varies greatly between individuals because of a differential masking effect of light.

Balancing Selection in Structured Populations of Three-Spined Sticklebacks

Balancing selection describes evolutionary processes that maintain genetic diversity. To date, the number of impacted genes and underlying biological functions remain elusive. Using 60 three-spined stickleback genomes (Gasterosteus aculeatus) from five recently diverged lake-river population-pairs, we performed genome-wide scans across two levels of organization: population-pairs and populations. We overlapped Tajima’s D and Watterson’s estimator metrics and verified signals with additional summary statistics, and evaluated alternative explanations: neutral evolution, population structure, associated overdominance, or demographic change. Candidate windows exhibiting signals of balancing selection spanned 2.31% (population-pair) and 3.10% (population) of the autosomes. These candidate windows had extended linkage disequilibrium and were enriched in intergenic and non-synonymous SNPs. We identified 715 (population-pair) and 1,010 (population) candidate genes under balancing selection. Importantly, using conservative thresholds, we found a small proportion of candidate genes overlapped with highly differentiated loci or regions of potential associated overdominance. There was little evidence of confounding effects originating from demographic change. Overall, candidate genes under balancing selection were associated with functions related to interactions with the environment (olfaction and receptor signalling pathways). Our results demonstrate selection that maintains standing genetic variation is common and evolves in response to local environmental pressures, playing an important role in adaptation.

An optimized LC-HRMS untargeted metabolomics workflow for multi-matrices investigations in the three-spined stickleback

Environmental metabolomics has become a growing research field to understand biological and biochemical perturbations of organisms in response to various abiotic or biotic stresses. It focuses on the comprehensive and systematic analysis of a biologic system’s metabolome. This allows the recognition of biochemical pathways impacted by a stressor, and the identification of some metabolites as biomarkers of potential perturbations occurring in a body. In this work, we describe the development and optimization of a complete reliable methodology based on liquid chromatography coupled to high resolution mass spectrometry (LC-HRMS) for untargeted metabolomics studies within a fish model species, the three-spined stickleback (Gasterosteus aculeatus). We evaluated the differences and also the complementarities between four different matrices (brain, gills, liver and whole fish) to obtain metabolome information. To this end, we optimized and compared sample preparation and the analytical method, since the type and number of metabolites detected in any matrix are closely related to these latter. For the sample preparation, a solid-liquid extraction was performed on a low quantity of whole fish, liver, brain, or gills tissues using combinations of methanol/water/heptane. Based on the numbers of features observed in LC-HRMS and on the responses of analytical standards representative of different metabolites groups (amino acids, sugars…), we discuss the influence of the nature, volume, and ratio of extraction solvents, the sample weight, and the reconstitution solvent. Moreover, the analytical conditions (LC columns, pH and additive of mobile phases and ionization modes) were also optimized so as to ensure the maximum metabolome coverages. Thus, two complementary chromatographic procedures were combined in order to cover a broader range of metabolites: a reversed phase separation (RPLC) on a C18 column followed by detection with positive ionization mode (ESI+) and a hydrophilic interaction chromatography (HILIC) on a zwitterionic column followed by detection with negative ionization mode (ESI-). This work provides information on brain, gills, liver, vs the whole body contribution to the stickleback metabolome. These information would help to guide ecotoxicological and biomonitoring studies.