Q(ST)-F(ST) comparisons: evolutionary and ecological insights from genomic heterogeneity.

Comparative studies of the divergence of quantitative traits and neutral molecular markers, known as Q(ST)-F(ST) comparisons, provide a means for researchers to distinguish between natural selection and genetic drift as causes of population differentiation in complex polygenic traits. The use of Q(ST)-F(ST) comparisons has increased rapidly in the last few years, highlighting the utility of this approach for addressing a wide range of questions that are relevant to evolutionary and ecological genetics. These studies have also provided lessons for the design of future Q(ST)-F(ST) comparisons. Methods based on the Q(ST)-F(ST) approach could also be used to analyse various types of 'omics' data in new and revealing ways.

Sticklebacks adapted to divergent osmotic environments show differences in plasticity for kidney morphology and candidate gene expression.

Novel physiological challenges in different environments can promote the evolution of divergent phenotypes, either through plastic or genetic changes. Environmental salinity serves as a key barrier to the distribution of nearly all aquatic organisms, and species diversification is likely to be enabled by adaptation to alternative osmotic environments. The threespine stickleback (Gasterosteus aculeatus) is a euryhaline species with populations found both in marine and freshwater environments. It has evolved both highly plastic and locally adapted phenotypes due to salinity-derived selection, but the physiological and genetic basis of adaptation to salinity is not fully understood. We integrated comparative cellular morphology of the kidney, a key organ for osmoregulation, and candidate gene expression to explore the underpinnings of evolved variation in osmotic plasticity within two populations of sticklebacks from distinct salinity zones in the Baltic Sea: the high salinity Kattegat, representative of the ancestral marine habitat; and the low salinity Bay of Bothnia. A common-garden experiment revealed that kidney morphology in the ancestral high-salinity population had a highly plastic response to salinity conditions whereas this plastic response was reduced in the low-salinity population. Candidate gene expression in kidney tissue revealed a similar pattern of population-specific differences, with a higher degree of plasticity in the native high-salinity population. Together these results suggest that renal cellular morphology has become canalized to low salinity, and that these structural differences may have functional implications for osmoregulation.

The evolution and adaptive potential of transcriptional variation in sticklebacks--signatures of selection and widespread heritability.

Evidence implicating differential gene expression as a significant driver of evolutionary novelty continues to accumulate, but our understanding of the underlying sources of variation in expression, both environmental and genetic, is wanting. Heritability in particular may be underestimated when inferred from genetic mapping studies, the predominant "genetical genomics" approach to the study of expression variation. Such uncertainty represents a fundamental limitation to testing for adaptive evolution at the transcriptomic level. By studying the inheritance of expression levels in 10,495 genes (10,527 splice variants) in a threespine stickleback pedigree consisting of 563 individuals, half of which were subjected to a thermal treatment, we show that 74-98% of transcripts exhibit significant additive genetic variance. Dominance variance is also prevalent (41-99% of transcripts), and genetic sources of variation seem to play a more significant role in expression variance in the liver than a key environmental variable, temperature. Among-population comparisons suggest that the majority of differential expression in the liver is likely due to neutral divergence; however, we also show that signatures of directional selection may be more prevalent than those of stabilizing selection. This predominantly aligns with the neutral model of evolution for gene expression but also suggests that natural selection may still act on transcriptional variation in the wild. As genetic variation both within- and among-populations ultimately defines adaptive potential, these results indicate that broad adaptive potential may be found within the transcriptome.

Variation of Tolerance to Isothiazolinones Among Daphnia pulex Clones.

Isothiazolinones are a family of broad-spectrum biocides widely used in industry and consumer products. Chloro- and methyl-isothiazolinones (CMIT and MIT) are documented as strong irritants, yet they are still used in a wide variety of applications, including cosmetics, cleansers, hygienic products, and various industrial applications. The subsequent substantial release of these molecules from urban sources into freshwater environments, and their potential impacts on aquatic species, have nevertheless received little attention so far, with few studies reporting on the toxicity of either CMIT or MIT to nontarget organisms. The present study addresses this current knowledge gap by evaluating the acute toxicity to Daphnia pulex (Cladocera) of CMIT/MIT (3:1) and MIT, the two formulations most commonly used by manufacturers. In addition, genetic diversity is known to be a major component of variability in phenotypic responses, although it is largely overlooked in typical toxicity tests. Thus the potential range of responses inherent to genetic diversity is rarely considered. Therefore, to account for intraspecific variations in sensitivity, our design involved eight clonal lines of D. pulex stemming from distinct natural populations or commercial strains. Clones exhibited strong variation in their responses, with median lethal concentration (LC50) values ranging from 0.10 to 1.84 mg/L for the mixture CMIT/MIT, and from 0.68 to 2.84 mg/L for MIT alone. These intraspecific ranges of LC50 values challenge the use of single clones of daphnids in standard ecotoxicological tests and the predictions based on their results. The present study brings new evidence that assessing ecological risk of chemicals while ignoring genotype diversity is neither ecologically relevant, nor a representative evaluation of the diversity of potential adverse outcomes. Environ Toxicol Chem 2023;42:805-814. © 2023 SETAC.

Putative causes and consequences of MHC variation within and between locally adapted stickleback demes.

Genes of the major histocompatibility complex (MHC) have been a source of considerable research interest, owing in large part to the growing body of evidence that they may be subject to both natural and sexual selection. However, much remains to be learned about the dynamics of MHC genes in subdivided populations, particularly those characterized by divergent ecological pressures. In this study, we attempt to disentangle the relative roles of both parasite-mediated selection and MHC-mediated mate choice in an open estuarine system inhabited by two parapatric, adaptively divergent threespine stickleback (Gasterosteus aculeatus) demes. We sequenced the putative peptide-binding region (PBR) of an estimated four Class IIß loci from 127 individuals, identifying 329 sequence variants (276 translated amino acid sequences). Demes differed significantly both in the frequency of MHC alleles and in the communities of helminth parasites infecting resident sticklebacks. Strong signatures of natural selection were inferred from analyses of codon substitutions, particularly in the derived (freshwater) rather than the ancestral (marine) deme. Relationships between parasite load and MHC diversity were indicative of balancing selection, but only within the freshwater deme. Signals of MHC-mediated mate choice were weak and differed significantly between demes. Moreover, MHC-mediated mate choice was significantly influenced by environmental salinity and appeared of secondary importance to tendencies towards assortative mating. We discuss the implications of these findings in respect to ecological adaptation and the potential demographic consequences of possible outcomes of MHC-mediated mate choice.

Intraspecific divergence in the lateral line system in the nine-spined stickleback (Pungitius pungitius).

The mechanosensory lateral line system of fishes is an important organ system conveying information crucial to individual fitness. Yet, our knowledge of lateral line diversity is almost exclusively based on interspecific studies, whereas intraspecific variability and possible population divergence have remained largely unexplored. We investigated lateral line system variability in four marine and five pond populations of nine-spined stickleback (Pungitius pungitius). We found significant differences in neuromast number between pond and marine fish. In particular, three of seventeen lateral line regions (viz. caudal peduncle superficial neuromasts; canal neuromasts from the anterior trunk and caudal peduncle) showed strong divergence between habitats. Similar results were obtained with laboratory-reared individuals from a subset of populations, suggesting that the patterns found in nature likely have a genetic basis. Interestingly, we also found habitat-dependent population divergence in neuromast variability, with pond populations showing greater heterogeneity than marine populations, although only in wild-caught fish. A comparison of neutral genetic (F(ST)) and phenotypic (P(ST)) differentiation suggested that natural selection is likely associated with habitat-dependent divergence in neuromast counts. Hence, the results align with the conclusion that the mechanosensory lateral line system divergence among marine and pond nine-spined sticklebacks is adaptive.

FishResp: R package and GUI application for analysis of aquatic respirometry data.

Intermittent-flow respirometry is widely used to measure oxygen uptake rates and subsequently estimate aerobic metabolic rates of aquatic animals. However, the lack of a standard quality-control software to detect technical problems represents a potential impediment to comparisons across studies in the field of evolutionary and conservation physiology. Here, we introduce 'FishResp', a versatile R package and its graphical implementation for quality-control and filtering of raw respirometry data. Our goal is to provide a straightforward, cross-platform and free software to help improve the quality and comparability of metabolic rate estimates for reducing methodological fragmentation in the field of aquatic respirometry. FishResp accepts data from various respirometry systems, allows users to detect potential mechanical problems which can occur during oxygen uptake measurements (e.g. chamber leaking, poor water circulation), and offers six options to correct raw data for microbial oxygen consumption. The software performs filtering of raw data based on user criteria, and produces accurate and unbiased estimates of absolute and mass-specific metabolic rates. Using data from three-spined sticklebacks (Gasterosteus aculeatus) and Trinidadian guppies (Poecilia reticulata), we demonstrate the virtues of FishResp, highlighting the importance of detecting mechanical problems and correcting measurements for background respiration.

Landscape genetic analyses reveal cryptic population structure and putative selection gradients in a large-scale estuarine environment.

Disentangling the relative contributions of selective and neutral processes underlying phenotypic and genetic variation under natural, environmental conditions remains a central challenge in evolutionary ecology. However, much of the variation that could be informative in this area of research is likely to be cryptic in nature; thus, the identification of wild populations suitable for study may be problematic. We use a landscape genetics approach to identify such populations of three-spined stickleback inhabiting the Saint Lawrence River estuary. We sampled 1865 adult fish over multiple years. Individuals were genotyped for nine microsatellite loci, and georeferenced multilocus data were used to infer population groupings, as well as locations of genetic discontinuities, under a Bayesian model framework (geneland). We modelled environmental data using nonparametric multiple regression to explain genetic differentiation as a function of spatio-ecological effects. Additionally, we used genotype data to estimate dispersal and gene flow to parameterize a simple model predicting adaptive vs. plastic divergence between demes. We demonstrate a bipartite division of the genetic landscape into freshwater and maritime zones, independent of geographical distance. Moreover, we show that the greatest proportion of genetic variation (31.5%) is explained by environmental differences. However, the potential for either adaptive or plastic divergence between demes is highly dependent upon the strength of migration and selection. Consequently, we highlight the utility of landscape genetics as a tool for hypothesis generation and experimental design, to identify focal populations and putative selection gradients, in order to distinguish between phenotypic plasticity and local adaptation.

Selection on the morphology-physiology-performance nexus: Lessons from freshwater stickleback morphs.

Conspecifics inhabiting divergent environments frequently differ in morphology, physiology, and performance, but the interrelationships amongst traits and with Darwinian fitness remains poorly understood. We investigated population differentiation in morphology, metabolic rate, and swimming performance in three-spined sticklebacks (Gasterosteus aculeatus L.), contrasting a marine/ancestral population with two distinct freshwater morphotypes derived from it: the "typical" low-plated morph, and a unique "small-plated" morph. We test the hypothesis that similar to plate loss in other freshwater populations, reduction in lateral plate size also evolved in response to selection. Additionally, we test how morphology, physiology, and performance have evolved in concert as a response to differences in selection between marine and freshwater environments. We raised pure-bred second-generation fish originating from three populations and quantified their lateral plate coverage, burst- and critical swimming speeds, as well as standard and active metabolic rates. Using a multivariate QST-FST framework, we detected signals of directional selection on metabolic physiology and lateral plate coverage, notably demonstrating that selection is responsible for the reduction in lateral plate coverage in a small-plated stickleback population. We also uncovered signals of multivariate selection amongst all bivariate trait combinations except the two metrics of swimming performance. Divergence between the freshwater and marine populations exceeded neutral expectation in morphology and in most physiological and performance traits, indicating that adaptation to freshwater habitats has occurred, but through different combinations of traits in different populations. These results highlight both the complex interplay between morphology, physiology and performance in local adaptation, and a framework for their investigation.

Challenges and advances for transcriptome assembly in non-model species.

Analyses of high-throughput transcriptome sequences of non-model organisms are based on two main approaches: de novo assembly and genome-guided assembly using mapping to assign reads prior to assembly. Given the limits of mapping reads to a reference when it is highly divergent, as is frequently the case for non-model species, we evaluate whether using blastn would outperform mapping methods for read assignment in such situations (>15% divergence). We demonstrate its high performance by using simulated reads of lengths corresponding to those generated by the most common sequencing platforms, and over a realistic range of genetic divergence (0% to 30% divergence). Here we focus on gene identification and not on resolving the whole set of transcripts (i.e. the complete transcriptome). For simulated datasets, the transcriptome-guided assembly based on blastn recovers 94.8% of genes irrespective of read length at 0% divergence; however, assignment rate of reads is negatively correlated with both increasing divergence level and reducing read lengths. Nevertheless, we still observe 92.6% of recovered genes at 30% divergence irrespective of read length. This analysis also produces a categorization of genes relative to their assignment, and suggests guidelines for data processing prior to analyses of comparative transcriptomics and gene expression to minimize potential inferential bias associated with incorrect transcript assignment. We also compare the performances of de novo assembly alone vs in combination with a transcriptome-guided assembly based on blastn both via simulation and empirically, using data from a cyprinid fish species and from an oak species. For any simulated scenario, the transcriptome-guided assembly using blastn outperforms the de novo approach alone, including when the divergence level is beyond the reach of traditional mapping methods. Combining de novo assembly and a related reference transcriptome for read assignment also addresses the bias/error in contigs caused by the dependence on a related reference alone. Empirical data corroborate these findings when assembling transcriptomes from the two non-model organisms: Parachondrostoma toxostoma (fish) and Quercus pubescens (plant). For the fish species, out of the 31,944 genes known from D. rerio, the guided and de novo assemblies recover respectively 20,605 and 20,032 genes but the performance of the guided assembly approach is much higher for both the contiguity and completeness metrics. For the oak, out of the 29,971 genes known from Vitis vinifera, the transcriptome-guided and de novo assemblies display similar performance, but the new guided approach detects 16,326 genes where the de novo assembly only detects 9,385 genes.

Regulatory Architecture of Gene Expression Variation in the Threespine Stickleback Gasterosteus aculeatus.

Much adaptive evolutionary change is underlain by mutational variation in regions of the genome that regulate gene expression rather than in the coding regions of the genes themselves. An understanding of the role of gene expression variation in facilitating local adaptation will be aided by an understanding of underlying regulatory networks. Here, we characterize the genetic architecture of gene expression variation in the threespine stickleback (Gasterosteus aculeatus), an important model in the study of adaptive evolution. We collected transcriptomic and genomic data from 60 half-sib families using an expression microarray and genotyping-by-sequencing, and located expression quantitative trait loci (eQTL) underlying the variation in gene expression in liver tissue using an interval mapping approach. We identified eQTL for several thousand expression traits. Expression was influenced by polymorphism in both cis- and trans-regulatory regions. Trans-eQTL clustered into hotspots. We did not identify master transcriptional regulators in hotspot locations: rather, the presence of hotspots may be driven by complex interactions between multiple transcription factors. One observed hotspot colocated with a QTL recently found to underlie salinity tolerance in the threespine stickleback. However, most other observed hotspots did not colocate with regions of the genome known to be involved in adaptive divergence between marine and freshwater habitats.

The adaptive potential of subtropical rainbowfish in the face of climate change: heritability and heritable plasticity for the expression of candidate genes.

Whilst adaptation and phenotypic plasticity might buffer species against habitat degradation associated with global climate change, few studies making such claims also possess the necessary and sufficient data to support them. Doing so requires demonstration of heritable variation in traits affecting fitness under new environmental conditions. We address this issue using an emerging aquatic system to study adaptation to climate change, the crimson-spotted rainbowfish (Melanotaenia duboulayi), a freshwater species from a region of eastern Australia projected to be affected by marked temperature increases. Captive born M. duboulayi of known pedigree were used to assess the long-term effects of contemporary and 2070-projected summer temperatures on the expression of genes previously identified in a climate change transcriptomics (RNA-Seq) experiment. Nearly all genes responded to increasing temperature. Significant additive genetic variance explained a moderate proportion of transcriptional variation for all genes. Most genes also showed broad-sense genetic variation in transcriptional plasticity. Additionally, molecular pathways of candidate genes co-occur with genes inferred to be under climate-mediated selection in wild M. duboulayi populations. Together, these results indicate the presence of existing variation in important physiological traits, and the potential for adaptive responses to a changing thermal environment.

Hypoxia and the pharmaceutical diclofenac influence the circadian responses of three-spined stickleback.

Pollution with low concentrations of pharmaceuticals, especially when combined with low-oxygen conditions (hypoxia), is a threat to aquatic ecosystems worldwide. The non-steroidal anti-inflammatory drug diclofenac is commonly detected in wastewater effluents, and has potential to accumulate in the bile of fish. Diclofenac has been shown to activate aryl hydrocarbon receptor (AHR), which induces transcription in the metabolic enzyme cytochrome P450 1a (cyp1a). Previously, crosstalk has been shown to occur between AHR and hypoxia inducible factor 1 (HIF-1). In addition, both of these transcription factors interact with the proteins regulating circadian (24-h) rhythms in vertebrates. Yet little is known about the significance of these interactions during simultaneous exposure to chemicals and hypoxia in fish in vivo. We exposed wild-caught three-spined sticklebacks (Gasterosteus aculeatus) to diclofenac (1 µg/L, 14 days), hypoxia (2.0 mg/L, up to 24h) and the combination of both. We then analyzed markers of chemical biotransformation (EROD activity, cyp1a and ahr mRNA levels), glycolysis (lactate dehydrogenase (LDH) enzyme activity, ldh and enolase 1a mRNA levels), and the transcription of core circadian clock genes clock and period 1 in liver tissue. Samples were taken at three time points during the light period in order to address disturbances in the circadian variation of metabolic processes. The results show that mRNA levels and LDH activity tended to be lowest before the dark period, but this pattern was disturbed by hypoxia and diclofenac. Diclofenac and hypoxia co-exposure induced EROD activity more strongly than diclofenac exposure alone, while cyp1a mRNA level was increased also by hypoxia and diclofenac alone. LDH activity and mRNA expression showed a clear time-dependent response during hypoxia, which is consistent with the previously suggested decreased accumulation of HIF-1 during the dark period. Furthermore, LDH activity and transcription was disturbed by diclofenac, indicating important effects of environmental pollutants in disturbing natural acclimation. This study demonstrates the need for more studies to understand the potential disturbances in endogenous rhythms caused by environmental pollution in natural populations.

Multiple evolutionary pathways to decreased lateral plate coverage in freshwater threespine sticklebacks.

Adaptation to novel environments can be based either on standing genetic variation or variation attributable to new mutations. When standing genetic variation for a functional adaptation is lacking, and variation due to new mutations is not yet available, adaptation is possible only through alternative functional solutions. Reduction in the number of bony lateral plates as an adaptation to freshwater colonization by marine threespine sticklebacks (Gasterosteus aculeatus) has occurred in numerous independent cases through allelic substitution in the ectodysplasin-a (Eda) gene. Studying the phenotypic and genetic variation in plate number and size in five marine and six freshwater threespine stickleback populations, we found that when variation in Eda was limiting (i.e., alleles associated with the low-plate morph were missing or in extremely low frequency), plate number reduction did not take place in freshwater populations, but reduced lateral plate coverage was achieved by a reduction in the size of lateral plates. Our results suggest that this phenotypically and genetically discrete "small-plated" threespine stickleback-which is the dominant form in three northern European freshwater populations-may be functionally equivalent to the low-plated morph and hence, serve as an example of convergent evolution toward functional similarity in the face of genetic constraints.

Adaptive divergence between freshwater and marine sticklebacks: insights into the role of phenotypic plasticity from an integrated analysis of candidate gene expression.

Debate surrounding the integration of phenotypic plasticity within the neo-Darwinian paradigm has recently intensified, but is largely dominated by conceptual abstractions. Advances in our capacities to identify candidate genes, and quantify their levels of expression, now facilitate the study of natural variation in inherently plastic traits, and may lead to a more concrete understanding of plasticity's role in adaptive evolution. We present data from parapatric threespine stickleback (Gasterosteus aculeatus) demes inhabiting geologically recent, freshwater and saltwater zones of a large estuary. Reaction norms for survival confirm adaptation to local salinity conditions. Analysis of osmoregulatory candidate gene expression within an ecological quantitative genetics framework suggests putative mechanisms underlying adaptive variation, and provides insights into the role of ancestral trait plasticity in this divergence. A sodium-potassium ATPase (ATP1A1) is identified as a candidate gene for freshwater adaptation. In addition to heritable variation for gene expression, we infer significant correlation between measures of expression and individual fitness. Overall results indicate a loss of plasticity in the freshwater deme. We discuss how this is consistent with adaptation facilitated by ancestral plasticity as a heuristic example that may prove useful for future, explicit tests of the genetic assimilation hypothesis.

Habitat and home range fidelity in a trophically dimorphic pumpkinseed sunfish (Lepomis gibbosus) population.

We investigated habitat selection in a trophically dimorphic population of pumpkinseed sunfish (Lepomis gibbosus) to determine whether littoral and limnetic ecotypes exhibit habitat or site fidelity. A transplant experiment was conducted, in which 998 pumpkinseeds captured from littoral and limnetic sites were marked and released in either the site of capture, the nearest site of the same habitat type, or the nearest site of the opposite habitat type. Daily recapture attempts over the course of the reproductive and growing season provided a 25% recapture rate, 40% of which were recaptured on multiple occasions at the same site. Site fidelity was very high in both ecotypes. Results estimated with a multi-state transition model indicated that the probability of a transplanted pumpkinseed returning to its site of origin ranged from 74% for limnetic pumpkinseeds released into a different limnetic or littoral site, to 93% for littoral pumpkinseeds released into a limnetic site. Furthermore, the probability of a pumpkinseed being recaptured at its site of origin if not transplanted was estimated at 97 and 98% for limnetic individuals and littoral individuals, respectively. Discriminant Function Analysis of helminth parasite loads sampled from littoral and limnetic individuals could classify site of origin with 96-100% accuracy, suggesting that the habitat and site fidelity patterns observed with mark-recapture are indicative of long-term habitat segregation of the two forms. The results of our experiment provide compelling evidence of correlated habitat selection as a function of home range fidelity within both ecotypes of a subtly dimorphic species. Such behaviour could have a significant effect on present or future gene flow.