Rapid, broad-scale gene expression evolution in experimentally harvested fish populations

Directional selection, not the direction of selection, affects telomere length and copy number at ribosomal RNA loci

Many fisheries exert directional selection on traits such as body size and growth rate. Whether directional selection impacts regions of the genome associated with traits related to growth is unknown. To address this issue, we characterised copy number variation in three regions of the genome associated with cell division, (1) telomeric DNA, (2) loci transcribed as ribosomal RNA (rDNA), and (3) mitochondrial DNA (mtDNA), in three selection lines of zebrafish reared at three temperatures (22 °C, 28 °C, and 34 °C). Selection lines differed in (1) the direction of selection (two lines experienced directional selection for large or small body size) and (2) whether they experienced any directional selection itself. Lines that had experienced directional selection were smaller, had lower growth rate, shorter telomeres, and lower rDNA copy number than the line that experiencing no directional selection. Neither telomere length nor rDNA copy number were affected by temperature. In contrast, mtDNA content increased at elevated temperature but did not differ among selection lines. Though directional selection impacts rDNA and telomere length, direction of such selection did not matter, whereas mtDNA acts as a stress marker for temperature. Future work should examine the consequences of these genomic changes in natural fish stocks.

Does size-selective harvesting erode adaptive potential to thermal stress?

Overharvesting is a serious threat to many fish populations. High mortality and directional selection on body size can cause evolutionary change in exploited populations via selection for a specific phenotype and a potential reduction in phenotypic diversity. Whether the loss of phenotypic diversity that accompanies directional selection impairs response to environmental stress is not known. To address this question, we exposed three zebrafish selection lines to thermal stress. Two lines had experienced directional selection for (1) large and (2) small body size, and one was (3) subject to random removal of individuals with respect to body size (i.e. line with no directional selection). Selection lines were exposed to three temperatures (elevated, 34°C; ambient, 28°C; low, 22°C) to determine the response to an environmental stressor (thermal stress). We assessed differences among selection lines in their life history (growth and reproduction), physiological traits (metabolic rate and critical thermal max) and behaviour (activity and feeding behaviour) when reared at different temperatures. Lines experiencing directional selection (i.e. size selected) showed reduced growth rate and a shift in average phenotype in response to lower or elevated thermal stress compared with fish from the random-selected line. Our data indicate that populations exposed to directional selection can have a more limited capacity to respond to thermal stress compared with fish that experience a comparable reduction in population size (but without directional selection). Future studies should aim to understand the impacts of environmental stressors on natural fish stocks.

Environment and genotype predict the genomic nature of domestication of salmonids as revealed by gene expression

Billions of salmonids are produced annually by artificial reproduction for harvest and conservation. Morphologically, behaviourally and physiologically these fish differ from wild-born fish, including in ways consistent with domestication. Unlike most studied domesticates, which diverged from wild ancestors millennia ago, salmonids offer a tractable model for early-stage domestication. Here, we review a fundamental mechanism for domestication-driven differences in early-stage domestication, differentially expressed genes (DEGs), in salmonids. We found 34 publications examining DEGs under domestication driven by environment and genotype, covering six species, over a range of life-history stages and tissues. Three trends emerged. First, domesticated genotypes have increased expression of growth hormone and related metabolic genes, with differences magnified under artificial environments with increased food. Regulatory consequences of these DEGs potentially drive overall DEG patterns. Second, immune genes are often DEGs under domestication and not simply owing to release from growth-immune trade-offs under increased food. Third, domesticated genotypes exhibit reduced gene expression plasticity, with plasticity further reduced in low-complexity environments typical of production systems. Recommendations for experimental design improvements, coupled with tissue-specific expression and emerging analytical approaches for DEGs present tractable avenues to understand the evolution of domestication in salmonids and other species.

Evolutionary implications of size-selective mortality on the ontogenetic development of shoal cohesion: a neurochemical approach using a zebrafish, Danio rerio, harvest selection experiment

Abstract

Size-selective mortality may evolutionarily alter life-history as well as individual behavioral and physiological traits. Moreover, size-selective mortality can affect group behavioral traits, such as shoaling and collective properties (e.g., shoal cohesion), which are relevant for finding food and reducing risk of predation. Here, we present experimental evidence using selection lines of zebrafish (Danio rerio) that were exposed to positive (large-harvested), negative (small-harvested), and random (control) size-selective mortality for five generations, followed by eight generations during which harvesting was halted to remove maternal effects and to study evolutionarily fixed outcomes. We investigated changes in shoal cohesion and turnover in monoamines in zebrafish through ontogeny. To that end, we repeatedly measured inter-individual distance in groups of eight fish and the turnovers of dopamine and serotonin in brains of fish from juvenile to the adult stage at 40-day intervals. We, firstly, found that shoal cohesion was overall consistent through ontogeny at group levels suggesting the presence of collective personality. Secondly, we found a decrease in shoal cohesion through ontogeny in the small-harvested and control lines, while the large-harvested line did not show any ontogenetic change. Thirdly, the selection lines did not differ among each other in shoal cohesion at any ontogenetic stage. Fourthly, dopamine turnover increased through ontogeny in a similar way for all lines while the serotonin turnover decreased in the large-harvested and control lines, but not in the small-harvested line. The large-harvested line also had higher serotonin turnover than controls at specific time periods. In conclusion, intensive size-selective mortality left an evolutionary legacy of asymmetric selection responses in the ontogeny of shoal cohesion and the underlying physiological mechanisms in experimentally harvested zebrafish in the laboratory.

Significant statement

The evolution of animal behavior can be affected by human activities both at behavioral and physiological levels, but causal evidence is scarce and mostly focusing on single life-stages. We studied whether and to what extent size-selective harvesting, a common selection pattern in fisheries, can be an evolutionary driver of the development of shoal cohesion during ontogeny. We used a multi-generation experiment with zebrafish to study cause-and-effects of opposing size-selection patterns. We quantified shoal cohesion, and serotonin and dopamine turnover in the brain. We found that shoal cohesion emerged as a collective personality trait and that behavioral and physiological responses were asymmetrical with respect to the opposing selection patterns.

Size-selective mortality fosters ontogenetic changes in collective risk-taking behaviour in zebrafish, Danio rerio

Size-selective mortality is common in fish populations and can operate either in a positive size-selective fashion by harvesting larger-than-average fish or be negatively size-selective by harvesting smaller-than-average fish. Through various mechanisms (like genetic correlations among behaviour and life-history traits or direct selection on behaviour co-varying with growth rate or size-at-maturation), size-selection can result in evolutionary changes in behavioural traits. Theory suggests that both positive and negative size-selection without additional selection on behaviour favours boldness, while evolution of shyness is possible if the largest fish are harvested. Here we examined the impact of size-selective mortality on collective boldness across ontogeny using three experimental lines of zebrafish (Daniorerio) generated through positive (large-harvested), negative (small-harvested) and random (control line) size-selective mortality for five generations and then relaxed selection for 10 generations to examine evolutionarily fixed outcomes. We measured collective risk-taking during feeding (boldness) under simulated aerial predation threat, and across four contexts in presence/absence of a cichlid. Boldness decreased across ontogeny under aerial predation threat, and the small-harvested line was consistently bolder than controls. The large and small-harvested lines showed higher behavioural plasticity as larvae and developed personality earlier compared to the controls. The large-harvested line showed increased variability and plasticity in boldness throughout ontogeny. In the presence of a live predator, fish did not differ in boldness in three contexts compared to the controls, but the large-harvested line showed reduced behavioural plasticity across contexts than controls. Our results confirmed theory by demonstrating that size-selective harvesting evolutionarily alters collective boldness and its variability and plasticity.

Captivity induces large and population-dependent brain transcriptomic changes in wild-caught cane toads (Rhinella marina)

Gene expression levels are key molecular phenotypes at the interplay between genotype and environment. Mounting evidence suggests that short-term changes in environmental conditions, such as those encountered in captivity, can substantially affect gene expression levels. Yet, the exact magnitude of this effect, how general it is, and whether it results in parallel changes across populations are not well understood. Here, we take advantage of the well-studied cane toad, Rhinella marina, to examine the effect of short-term captivity on brain gene expression levels, and determine whether effects of captivity differ between long-colonized and vanguard populations of the cane toad's Australian invasion range. We compared the transcriptomes of wild-caught toads immediately assayed with those from toads captured from the same populations but maintained in captivity for seven months. We found large differences in gene expression levels between captive and wild-caught toads from the same population, with an over-representation of processes related to behaviour and the response to stress. Captivity had a much larger effect on both gene expression levels and gene expression variability in toads from vanguard populations compared to toads from long-colonized areas, potentially indicating an increased plasticity in toads at the leading edge of the invasion. Overall, our findings indicate that short-term captivity can induce large and population-specific transcriptomic changes, which has significant implications for studies comparing phenotypic traits of wild-caught organisms from different populations that have been held in captivity.

A long evolutionary reach for fishing nets

Earlier maturation of Atlantic salmon is linked to indirect effects of fisheries on its prey.

The Effects of Selective Harvest on Japanese Spanish Mackerel (Scomberomorus niphonius) Phenotypic Evolution

Japanese Spanish mackerel (Scomberomorus niphonius) is an important fish species in the China Seas with wide distribution, extensive migration, and high economic value. This species has been yielding high fisheries production despite experiencing continuously high fishing pressure and the conversion from gillnet to trawl harvesting. Meanwhile, changes in life-history traits have been observed, including earlier maturation and smaller size at age. Here, we build an individual-based eco-genetic model parameterized for Japanese Spanish mackerel to investigate the population’s response to different fishing scenarios (fishing by trawl or by gillnet). The model allows evolution of life-history processes including maturation, reproduction and growth. It also incorporates environmental variability, phenotypic plasticity, and density-dependent feedbacks. Our results show that different gear types can result in different responses of life-history traits and altered population dynamics. The population harvested by gillnet shows weaker response to fishing than that by trawl. When fishing ceases, gillnet-harvested population can recover to the pre-harvest level more easily than that harvested by trawl. The different responses of population growth rate and evolution to different fishing gears demonstrated in this study shed light on the sustainable management and utilization of Japanese Spanish mackerel in the over-exploited China Seas.

Capture vulnerability of sea turtles on recreational fishing piers

Capture vulnerability of commercial and recreational fishes has been associated with behavioral, morphological, and life-history traits; however, relationships with non-target species, such as sea turtles, have not been adequately studied. We examined species composition, timing of captures, morphological variables including body size and head width, and body condition of sea turtles captured from a recreational fishing pier in the northern Gulf of Mexico and of sea turtles captured in the waters adjacent to the pier. From 2014 to 2019, 148 net captures and 112 pier captures of three sea turtle species were documented. Green turtles were captured most frequently in the net and on the pier. Turtles captured from the pier were larger than those captured in the net. There was no difference in head width between net-caught and pier-caught turtles; however, small sample sizes limited those comparisons. The body condition index was lower for pier-caught than net-caught Kemp';s ridleys but did not differ with green turtles or loggerheads. Differences were also observed in the timing of capture on the pier as compared to in the net. Finally, the relationship between size, body condition, and pier-capture vulnerability suggests these are complex interactions. Mortality of sea turtles captured from fishing piers could be selecting against bolder individuals, which may result in changes in sea turtle population demographics over a long time period.

A cis-Acting Mutation in the PxABCG1 Promoter Is Associated with Cry1Ac Resistance in Plutella xylostella (L.)

The molecular mechanisms of insect resistance to Cry toxins generated from the bacterium Bacillus thuringiensis (Bt) urgently need to be elucidated to enable the improvement and sustainability of Bt-based products. Although downregulation of the expression of midgut receptor genes is a pivotal mechanism of insect resistance to Bt Cry toxins, the underlying transcriptional regulation of these genes remains elusive. Herein, we unraveled the regulatory mechanism of the downregulation of the ABC transporter gene PxABCG1 (also called Pxwhite), a functional midgut receptor of the Bt Cry1Ac toxin in Plutella xylostella. The PxABCG1 promoters of Cry1Ac-susceptible and Cry1Ac-resistant strains were cloned and analyzed, and they showed clear differences in activity. Subsequently, a dual-luciferase reporter assay, a yeast one-hybrid (Y1H) assay, and RNA interference (RNAi) experiments demonstrated that a cis-mutation in a binding site of the Hox transcription factor Antennapedia (Antp) decreased the promoter activity of the resistant strain and eliminated the binding and regulation of Antp, thereby enhancing the resistance of P. xylostella to the Cry1Ac toxin. These results advance our knowledge of the roles of cis- and trans-regulatory variations in the regulation of midgut Cry receptor genes and the evolution of Bt resistance, contributing to a more complete understanding of the Bt resistance mechanism.

Genomic stability through time despite decades of exploitation in cod on both sides of the Atlantic

The mode and extent of rapid evolution and genomic change in response to human harvesting are key conservation issues. Although experiments and models have shown a high potential for both genetic and phenotypic change in response to fishing, empirical examples of genetic responses in wild populations are rare. Here, we compare whole-genome sequence data of Atlantic cod (Gadus morhua) that were collected before (early 20th century) and after (early 21st century) periods of intensive exploitation and rapid decline in the age of maturation from two geographically distinct populations in Newfoundland, Canada, and the northeast Arctic, Norway. Our temporal, genome-wide analyses of 346,290 loci show no substantial loss of genetic diversity and high effective population sizes. Moreover, we do not find distinct signals of strong selective sweeps anywhere in the genome, although we cannot rule out the possibility of highly polygenic evolution. Our observations suggest that phenotypic change in these populations is not constrained by irreversible loss of genomic variation and thus imply that former traits could be reestablished with demographic recovery.

Size Selective Harvesting Does Not Result in Reproductive Isolation among Experimental Lines of Zebrafish, Danio rerio: Implications for Managing Harvest-Induced Evolution

Simple Summary

Mortality in fish populations is commonly size-selective. In fisheries, larger fish are preferentially caught while natural predators preferentially consume smaller fish. Removal of certain sized fish from populations and elevated fishing mortality constitute a selection pressure which may change life-history, behaviour and reduce adult body-size. Because behaviour and body-size are related and influence mating preferences and reproductive output, size-selective mortality may favour subpopulations that less readily mate with each other. Our aim is to test this possibility using three experimental lines of zebrafish (Danio rerio) generated in laboratory by removing large-sized, small-sized and random-sized fish for five generations. We tested mating preferences among males and females and tested if they spawned together. We found males and females of all subpopulations to reproduce among themselves. Females generally preferred large-sized males. Females of all lines spawned with males, and males of all lines fertilised eggs of females independent of the subpopulation origin. Our study shows that size-selective mortality typical of fisheries or in populations facing heavy predation does not result in evolution of reproductive barriers. Thus, when populations adapted to fishing pressure come in contact with populations unexposed to such pressures, interbreeding may happen thereby helping exploited populations recover from harvest-induced evolution.

Abstract

Size-selective mortality is common in fish stocks. Positive size-selection happens in fisheries where larger size classes are preferentially targeted while gape-limited natural predation may cause negative size-selection for smaller size classes. As body size and correlated behavioural traits are sexually selected, harvest-induced trait changes may promote prezygotic reproductive barriers among selection lines experiencing differential size-selective mortality. To investigate this, we used three experimental lines of zebrafish (Danio rerio) exposed to positive (large-harvested), negative (small-harvested) and random (control line) size-selective mortality for five generations. We tested prezygotic preferences through choice tests and spawning trials. In the preference tests without controlling for body size, we found that females of all lines preferred males of the generally larger small-harvested line. When the body size of stimulus fish was statistically controlled, this preference disappeared and a weak evidence of line-assortative preference emerged, but only among large-harvested line fish. In subsequent spawning trials, we did not find evidence for line-assortative reproductive allocation in any of the lines. Our study suggests that size-selection due to fisheries or natural predation does not result in reproductive isolation. Gene flow between wild-populations and populations adapted to size-selected mortality may happen during secondary contact which can speed up trait recovery.

Evolutionary responses of marine organisms to urbanized seascapes

Many of the world's major cities are located in coastal zones, resulting in urban and industrial impacts on adjacent marine ecosystems. These pressures, which include pollutants, sewage, runoff and debris, temperature increases, hardened shorelines/structures, and light and acoustic pollution, have resulted in new evolutionary landscapes for coastal marine organisms. Marine environmental changes influenced by urbanization may create new selective regimes or may influence neutral evolution via impacts on gene flow or partitioning of genetic diversity across seascapes. While some urban selective pressures, such as hardened surfaces, are similar to those experienced by terrestrial species, others, such as oxidative stress, are specific to aquatic environments. Moreover, spatial and temporal scales of evolutionary responses may differ in the ocean due to the spatial extent of selective pressures and greater capacity for dispersal/gene flow. Here, we present a conceptual framework and synthesis of current research on evolutionary responses of marine organisms to urban pressures. We review urban impacts on genetic diversity and gene flow and examine evidence that marine species are adapting, or are predicted to adapt, to urbanization over rapid evolutionary time frames. Our findings indicate that in the majority of studies, urban stressors are correlated with reduced genetic diversity. Genetic structure is often increased in urbanized settings, but artificial structures can also act as stepping stones for some hard-surface specialists, promoting range expansion. Most evidence for rapid adaptation to urban stressors comes from studies of heritable tolerance to pollutants in a relatively small number of species; however, the majority of marine ecotoxicology studies do not test directly for heritability. Finally, we highlight current gaps in our understanding of evolutionary processes in marine urban environments and present a framework for future research to address these gaps.

Eco-evolutionary dynamics driven by fishing: From single species models to dynamic evolution within complex food webs

Evidence of contemporary evolution across ecological time scales stimulated research on the eco-evolutionary dynamics of natural populations. Aquatic systems provide a good setting to study eco-evolutionary dynamics owing to a wealth of long-term monitoring data and the detected trends in fish life-history traits across intensively harvested marine and freshwater systems. In the present study, we focus on modelling approaches to simulate eco-evolutionary dynamics of fishes and their ecosystems. Firstly, we review the development of modelling from single species to multispecies approaches. Secondly, we advance the current state-of-the-art methodology by implementing evolution of life-history traits of a top predator into the context of complex food web dynamics as described by the allometric trophic network (ATN) framework. The functioning of our newly developed eco-evolutionary ATNE framework is illustrated using a well-studied lake food web. Our simulations show how both natural selection arising from feeding interactions and size-selective fishing cause evolutionary changes in the top predator and how those feed back to its prey species and further cascade down to lower trophic levels. Finally, we discuss future directions, particularly the need to integrate genomic discoveries into eco-evolutionary projections.

Size-selective mortality induces evolutionary changes in group risk-taking behaviour and the circadian system in a fish

Intensive and trait-selective mortality of fish and wildlife can cause evolutionary changes in a range of life-history and behavioural traits. These changes might in turn alter the circadian system due to co-evolutionary mechanisms or correlated selection responses both at behavioural and molecular levels, with knock-on effects on daily physiological processes and behavioural outputs. We examined the evolutionary impact of size-selective mortality on group risk-taking behaviour and the circadian system in a model fish species. We exposed zebrafish Danio rerio to either large or small size-selective harvesting relative to a control over five generations, followed by eight generations during which harvesting was halted to remove maternal effects. Size-selective mortality affected fine-scale timing of behaviours. In particular, small size-selective mortality, typical of specialized fisheries and gape-limited predators targeting smaller size classes, increased group risk-taking behaviuor during feeding and after simulated predator attacks. Moreover, small size-selective mortality increased early peaks of daily activity as well as extended self-feeding daily activity to the photophase compared to controls. By contrast large size-selective mortality, typical of most wild capture fisheries, only showed an almost significant effect of decreasing group risk-taking behaviour during the habituation phase and no clear changes in fine-scale timing of daily behavioural rhythms compared to controls. We also found changes in the molecular circadian core clockwork in response to both size-selective mortality treatments. These changes disappeared in the clock output pathway because both size-selected lines showed similar transcription profiles. This switch downstream to the molecular circadian core clockwork also resulted in similar overall behavioural rhythms (diurnal swimming and self-feeding in the last hours of darkness) independent of the underlying molecular clock. To conclude, our experimental harvest left an asymmetrical evolutionary legacy in group risk-taking behaviour and in fine-scale daily behavioural rhythms. Yet, the overall timing of activity showed evolutionary resistance probably maintained by a molecular switch. Our experimental findings suggest that size-selective mortality can have consequences for behaviour and physiological processes.

Echolocation call frequency variation in horseshoe bats: molecular basis revealed by comparative transcriptomics

Recently diverged taxa with contrasting phenotypes offer opportunities for unravelling the genetic basis of phenotypic variation in nature. Horseshoe bats are a speciose group that exhibit a derived form of high-duty cycle echolocation in which the inner ear is finely tuned to echoes of the narrowband call frequency. Here, by focusing on three recently diverged subspecies of the intermediate horseshoe bat (Rhinolophus affinis) that display divergent echolocation call frequencies, we aim to identify candidate loci putatively involved in hearing frequency variation. We used de novo transcriptome sequencing of two mainland taxa (himalayanus and macrurus) and one island taxon (hainanus) to compare expression profiles of thousands of genes. By comparing taxa with divergent call frequencies (around 15 kHz difference), we identified 252 differentially expressed genes, of which six have been shown to be involved in hearing or deafness in human/mouse. To obtain further validation of these results, we applied quantitative reverse transcription-PCR to the candidate gene FBXL15 and found a broad association between the level of expression and call frequency across taxa. The genes identified here represent strong candidate loci associated with hearing frequency variation in bats.

The Effects of Quantitative Trait Architecture on Detection Power in Short-Term Artificial Selection Experiments

Evolve and resequence (E&R) experiments, in which artificial selection is imposed on organisms in a controlled environment, are becoming an increasingly accessible tool for studying the genetic basis of adaptation. Previous work has assessed how different experimental design parameters affect the power to detect the quantitative trait loci (QTL) that underlie adaptive responses in such experiments, but so far there has been little exploration of how this power varies with the genetic architecture of the evolving traits. In this study, we use forward simulation to build a more realistic model of an E&R experiment in which a quantitative polygenic trait experiences a short, but strong, episode of truncation selection. We study the expected power for QTL detection in such an experiment and how this power is influenced by different aspects of trait architecture, including the number of QTL affecting the trait, their starting frequencies, effect sizes, clustering along a chromosome, dominance, and epistasis patterns. We show that all of these parameters can affect allele frequency dynamics at the QTL and linked loci in complex and often unintuitive ways, and thus influence our power to detect them. One consequence of this is that existing detection methods based on models of independent selective sweeps at individual QTL often have lower detection power than a simple measurement of allele frequency differences before and after selection. Our findings highlight the importance of taking trait architecture into account when designing and interpreting studies of molecular adaptation with temporal data. We provide a customizable modeling framework that will enable researchers to easily simulate E&R experiments with different trait architectures and parameters tuned to their specific study system, allowing for assessment of expected detection power and optimization of experimental design.

Genetic analysis redraws the management boundaries for the European sprat

Sustainable fisheries management requires detailed knowledge of population genetic structure. The European sprat is an important commercial fish distributed from Morocco to the Arctic circle, Baltic, Mediterranean, and Black seas. Prior to 2018, annual catch advice on sprat from the International Council for the Exploration of the Sea (ICES) was based on five putative stocks: (a) North Sea, (b) Kattegat-Skagerrak and Norwegian fjords, (c) Baltic Sea, (d) West of Scotland-southern Celtic Seas, and (e) English Channel. However, there were concerns that the sprat advice on stock size estimates management plan inadequately reflected the underlying biological units. Here, we used ddRAD sequencing to develop 91 SNPs that were thereafter used to genotype approximately 2,500 fish from 40 locations. Three highly distinct and relatively homogenous genetic groups were identified: (a) Norwegian fjords; (b) Northeast Atlantic including the North Sea, Kattegat-Skagerrak, Celtic Sea, and Bay of Biscay; and (c) Baltic Sea. Evidence of genetic admixture and possibly physical mixing was detected in samples collected from the transition zone between the North and Baltic seas, but not between any of the other groups. These results have already been implemented by ICES with the decision to merge the North Sea and the Kattegat-Skagerrak sprat to be assessed as a single unit, thus demonstrating that genetic data can be rapidly absorbed to align harvest regimes and biological units.

Size reductions and genomic changes within two generations in wild walleye populations: associated with harvest?

The extent and rate of harvest‐induced genetic changes in natural populations may impact population productivity, recovery, and persistence. While there is substantial evidence for phenotypic changes in harvested fishes, knowledge of genetic change in the wild remains limited, as phenotypic and genetic data are seldom considered in tandem, and the number of generations needed for genetic changes to occur is not well understood. We quantified changes in size‐at‐age, sex‐specific changes in body size, and genomic metrics in three harvested walleye (Sander vitreus) populations and a fourth reference population with low harvest levels over a 15‐year period in Mistassini Lake, Quebec. We also collected Indigenous knowledge (IK) surrounding concerns about these populations over time. Using ~9,000 SNPs, genomic metrics included changes in population structure, neutral genomic diversity, effective population size, and signatures of selection. Indigenous knowledge revealed overall reductions in body size and number of fish caught. Smaller body size, a small reduction in size‐at‐age, nascent changes to population structure (population differentiation within one river and homogenization between two others), and signatures of selection between historical and contemporary samples reflected coupled phenotypic and genomic change in the three harvested populations in both sexes, while no change occurred in the reference population. Sex‐specific analyses revealed differences in both body size and genomic metrics but were inconclusive about whether one sex was disproportionately affected. Although alternative explanations cannot be ruled out, our collective results are consistent with the hypothesis that genetic changes associated with harvesting may arise within 1–2.5 generations in long‐lived wild fishes. This study thus demonstrates the need to investigate concerns about harvest‐induced evolution quickly once they have been raised.

Quantitative Trait Loci (QTL) Mapping

Quantitative trait loci (QTL) are genetic regions that influence phenotypic variation of a complex trait, often through genetic interactions with each other and the environment. These are commonly identified through a statistical genetic analysis known as QTL mapping. Here, I present a step-by-step, practical approach to QTL mapping along with a sample data file. I focus on methods commonly used and discoveries that have been made in fishes, and utilize a multiple QTL mapping (MQM) approach in the free software package R/qtl.

Among-family variation in survival and gene expression uncovers adaptive genetic variation in a threatened fish

Variation in among-family transcriptional responses to different environmental conditions can help to identify adaptive genetic variation, even prior to a selective event. Coupling differential gene expression with formal survival analyses allows for the disentanglement of treatment effects, required for understanding how individuals plastically respond to environmental stressors, from the adaptive genetic variation responsible for differential survival. We combined these two approaches to investigate responses to an emerging conservation issue, thiamine (vitamin B₁ ) deficiency, in a threatened population of Atlantic salmon (Salmo salar). Thiamine is an essential vitamin that is increasingly limited in many ecosystems. In Lake Champlain, Atlantic salmon cannot acquire thiamine in sufficient quantities to support natural reproduction; fertilized eggs must be reared in hatcheries and treated with supplemental thiamine. We evaluated transcriptional responses (via RNA sequencing) to thiamine treatment across families and found 3,616 genes differentially expressed between control (no supplemental thiamine) and treatment individuals. Fewer genes changed expression equally across families (i.e., additively) than exhibited genotype × environment interactions in response to thiamine. Differentially expressed genes were related to known physiological effects of thiamine deficiency, including oxidative stress, cardiovascular irregularities and neurological abnormalities. We also identified 1,446 putatively adaptive genes that were strongly associated with among-family survival in the absence of thiamine treatment, many of which related to neurogenesis and visual perception. Our results highlight the utility of coupling RNA sequencing with formal survival analyses to identify candidate genes that underlie the among-family variation in survival required for an adaptive response to natural selection.

Contrasting genomic shifts underlie parallel phenotypic evolution in response to fishing

Humans cause widespread evolutionary change in nature, but we still know little about the genomic basis of rapid adaptation in the Anthropocene. We tracked genomic changes across all protein-coding genes in experimental fish populations that evolved pronounced shifts in growth rates due to size-selective harvest over only four generations. Comparisons of replicate lines show parallel allele frequency shifts that recapitulate responses to size-selection gradients in the wild across hundreds of unlinked variants concentrated in growth-related genes. However, a supercluster of genes also rose rapidly in frequency and dominated the evolutionary dynamic in one replicate line but not in others. Parallel phenotypic changes thus masked highly divergent genomic responses to selection, illustrating how contingent rapid adaptation can be in the face of strong human-induced selection.

Are model organisms representative for climate change research? Testing thermal tolerance in wild and laboratory zebrafish populations

Model organisms can be useful for studying climate change impacts, but it is unclear whether domestication to laboratory conditions has altered their thermal tolerance and therefore how representative of wild populations they are. Zebrafish in the wild live in fluctuating thermal environments that potentially reach harmful temperatures. In the laboratory, zebrafish have gone through four decades of domestication and adaptation to stable optimal temperatures with few thermal extremes. If maintaining thermal tolerance is costly or if genetic traits promoting laboratory fitness at optimal temperature differ from genetic traits for high thermal tolerance, the thermal tolerance of laboratory zebrafish could be hypothesized to be lower than that of wild zebrafish. Furthermore, very little is known about the thermal environment of wild zebrafish and how close to their thermal limits they live. Here, we compared the acute upper thermal tolerance (critical thermal maxima; CT_max) of wild zebrafish measured on-site in West Bengal, India, to zebrafish at three laboratory acclimation/domestication levels: wild-caught, F₁ generation wild-caught and domesticated laboratory AB-WT line. We found that in the wild, CT_max increased with increasing site temperature. Yet at the warmest site, zebrafish lived very close to their thermal limit, suggesting that they may currently encounter lethal temperatures. In the laboratory, acclimation temperature appeared to have a stronger effect on CT_max than it did in the wild. The fish in the wild also had a 0.85-1.01°C lower CT_max compared to all laboratory populations. This difference between laboratory-held and wild populations shows that environmental conditions can affect zebrafish's thermal tolerance. However, there was no difference in CT_max between the laboratory-held populations regardless of the domestication duration. This suggests that thermal tolerance is maintained during domestication and highlights that experiments using domesticated laboratory-reared model species can be appropriate for addressing certain questions on thermal tolerance and global warming impacts.

Size-selective harvesting fosters adaptations in mating behaviour and reproductive allocation, affecting sexual selection in fish

The role of sexual selection in the context of harvest-induced evolution is poorly understood. However, elevated and trait-selective harvesting of wild populations may change sexually selected traits, which in turn can affect mate choice and reproduction. We experimentally evaluated the potential for fisheries-induced evolution of mating behaviour and reproductive allocation in fish. We used an experimental system of zebrafish (Danio rerio) lines exposed to large, small or random (i.e. control) size-selective mortality. The large-harvested line represented a treatment simulating the typical case in fisheries where the largest individuals are preferentially harvested. We used a full factorial design of spawning trials with size-matched individuals to control for the systematic impact of body size during reproduction, thereby singling out possible changes in mating behaviour and reproductive allocation. Both small size-selective mortality and large size-selective mortality left a legacy on male mating behaviour by elevating intersexual aggression. However, there was no evidence for line-assortative reproductive allocation. Females of all lines preferentially allocated eggs to the generally less aggressive males of the random-harvested control line. Females of the large-harvested line showed enhanced reproductive performance, and males of the large-harvested line had the highest egg fertilization rate among all males. These findings can be explained as an evolutionary adaptation by which individuals of the large-harvested line display an enhanced reproductive performance early in life to offset the increased probability of adult mortality due to harvest. Our results suggest that the large-harvested line evolved behaviourally mediated reproductive adaptations that could increase the rate of recovery when populations adapted to high fishing pressure come into secondary contact with other populations.

Introgressive hybridization between wild and domestic individuals and its relationship with parasitism in brook charr Salvelinus fontinalis

The effects of introgression on parasitism in brook charr Salvelinus fontinalis were investigated in 28 lakes with various levels of stocking in Québec, Canada. No effect of genetic background on parasitism was found at the individual level. Body length seemed to explain most of the variation observed at this level, with largest fish being more infected. However, lakes with the greater average domestic genetic background were found to display significantly lower parasite prevalence and diversity. Since our results indicate no effect of domestic genes at the individual level, the negative association with introgression found at the population level may be mainly attributed to differences in intrinsic environmental quality of lakes (e.g. fishing pressure, availability of food resources, abiotic characteristics).

Toward a mechanistic understanding of vulnerability to hook-and-line fishing: Boldness as the basic target of angling-induced selection

In passively operated fishing gear, boldness-related behaviors should fundamentally affect the vulnerability of individual fish and thus be under fisheries selection. To test this hypothesis, we used juvenile common-garden reared carp (Cyprinus carpio) within a narrow size range to investigate the mechanistic basis of behavioral selection caused by angling. We focused on one key personality trait (i.e., boldness), measured in groups within ponds, two morphological traits (body shape and head shape), and one life-history trait (juvenile growth capacity) and studied mean standardized selection gradients caused by angling. Carp behavior was highly repeatable within ponds. In the short term, over seven days of fishing, total length, not boldness, was the main predictor of angling vulnerability. However, after 20 days of fishing, boldness turned out to be the main trait under selection, followed by juvenile growth rate, while morphological traits were only weakly related to angling vulnerability. In addition, we found juvenile growth rate to be moderately correlated with boldness. Hence, direct selection on boldness will also induce indirect selection on juvenile growth and vice versa, but given that the two traits are not perfectly correlated, independent evolution of both traits is also possible. Our study is among the first to mechanistically reveal that energy-acquisition-related behaviors, and not growth rate per se, are key factors determining the probability of capture, and hence, behavioral traits appear to be the prime targets of angling selection. We predict an evolutionary response toward increased shyness in intensively angling-exploited fish stocks, possibly causing the emergence of a timidity syndrome.