RNA-seq reveals differential gene expression in the brains of juvenile resident and migratory smolt rainbow trout (Oncorhynchus mykiss)

Comparative genomics of rainbow trout (Oncorhynchus mykiss): Is the genetic architecture of migratory behavior conserved among populations?

Rainbow trout (Oncorhynchus mykiss) are a partially migratory species wherein some individuals undergo long-distance anadromous migrations, and others stay as residents in their native freshwater streams. The decision to migrate is known to be highly heritable, and yet, the underlying genes and alleles associated with migration are not fully characterized. Here we used a pooled approach of whole-genome sequence data from migratory and resident trout of two native populations-Sashin Creek, Alaska and Little Sheep Creek, Oregon-to obtain a genome-wide perspective of the genetic architecture of resident and migratory life history. We calculated estimates of genetic differentiation, genetic diversity, and selection between the two phenotypes to locate regions of interest and then compared these associations between populations. We identified numerous genes and alleles associated with life history development in the Sashin Creek population with a notable area on chromosome 8 that may play a critical role in the development of the migratory phenotype. However, very few alleles appeared to be associated with life history development in the Little Sheep Creek system, suggesting population-specific genetic effects are likely important in the development of anadromy. Our results indicate that a migratory life history is not controlled by a singular gene or region but supports the idea that there are many independent ways for a migratory phenotype to emerge in a population. Therefore, conserving and promoting genetic diversity in migratory individuals is paramount to conserving these populations. Ultimately, our data add to a growing body of literature that suggests that population-specific genetic effects, likely mediated through environmental variation, contribute to life history development in rainbow trout.

Whole-Genome Resequencing to Evaluate Life History Variation in Anadromous Migration of Oncorhynchus mykiss

Anadromous fish experience physiological modifications necessary to migrate between vastly different freshwater and marine environments, but some species such as Oncorhynchus mykiss demonstrate variation in life history strategies with some individuals remaining exclusively resident in freshwater, whereas others undergo anadromous migration. Because there is limited understanding of genes involved in this life history variation across populations of this species, we evaluated the genomic difference between known anadromous (n = 39) and resident (n = 78) Oncorhynchus mykiss collected from the Klickitat River, WA, USA, with whole-genome resequencing methods. Sequencing of these collections yielded 5.64 million single-nucleotide polymorphisms that were tested for significant differences between resident and anadromous groups along with previously identified candidate gene regions. Although a few regions of the genome were marginally significant, there was one region on chromosome Omy12 that provided the most consistent signal of association with anadromy near two annotated genes in the reference assembly: COP9 signalosome complex subunit 6 (CSN6) and NACHT, LRR, and PYD domain–containing protein 3 (NLRP3). Previously identified candidate genes for anadromy within the inversion region of chromosome Omy05 in coastal steelhead and rainbow trout were not informative for this population as shown in previous studies. Results indicate that the significant region on chromosome Omy12 may represent a minor effect gene for male anadromy and suggests that this life history variation in Oncorhynchus mykiss is more strongly driven by other mechanisms related to environmental rearing such as epigenetic modification, gene expression, and phenotypic plasticity. Further studies into regulatory mechanisms of this trait are needed to understand drivers of anadromy in populations of this protected species.

Back to the basics? Transcriptomics offers integrative insights into the role of space, time and the environment for gene expression and behaviour

Fuelled by the ongoing genomic revolution, broadscale RNA expression surveys are fast replacing studies targeting one or a few genes to understand the molecular basis of behaviour. Yet, the timescale of RNA-sequencing experiments and the dynamics of neural gene activation are insufficient to drive real-time switches between behavioural states. Moreover, the spatial, functional and transcriptional complexity of the brain (the most commonly targeted tissue in studies of behaviour) further complicates inference. We argue that a Central Dogma-like 'back-to-basics' assumption that gene expression changes cause behaviour leaves some of the most important aspects of gene-behaviour relationships unexplored, including the roles of environmental influences, timing and feedback from behaviour-and the environmental shifts it causes-to neural gene expression. No perfect experimental solutions exist but we advocate that explicit consideration, exploration and discussion of these factors will pave the way toward a richer understanding of the complicated relationships between genes, environments, brain gene expression and behaviour over developmental and evolutionary timescales.

Differential gene expression associated with behavioral variation in ecotypes of Lake Superior brook trout (Salvelinus fontinalis)

Associations between behaviors and the development of different life history tactics have been documented in several species of salmon, trout, and charr. While it is well known that such behaviors are heritable the genes and molecular pathways connected to these behaviors remain unknown. We used an RNA-seq approach to identify genes and molecular pathways differentially regulated in brain tissue between "shy" and "bold" brook trout (Salvelinus fontinalis). A small number of genes were differentially expressed between the behavioral types at several months after hatching and two years of age. Pathway analysis revealed that EIF2 signaling differed consistently between shy and bold individuals suggesting large-scale differences in protein synthesis between behavioral types in the brain. Additionally, the RNA-seq data were used to find polymorphisms within the brook trout genome and a GWAS approach was used to test for statistical associations between genetic variants and behavior type. One allele located in a transcription factor (TSHZ3) contained a protein-coding non-synonymous SNP suggesting that functional variation within TSHZ3 is connected to the development of different behaviors. These results suggest that the molecular basis of behavioral development is complex and due to the differential expression of many genes involved in a wide-range of different molecular pathways.

PacBio Iso-Seq Improves the Rainbow Trout Genome Annotation and Identifies Alternative Splicing Associated With Economically Important Phenotypes

Rainbow trout is an important model organism that has received concerted international efforts to study the transcriptome. For this purpose, short-read sequencing has been primarily used over the past decade. However, these sequences are too short of resolving the transcriptome complexity. This study reported a first full-length transcriptome assembly of the rainbow trout using single-molecule long-read isoform sequencing (Iso-Seq). Extensive computational approaches were used to refine and validate the reconstructed transcriptome. The study identified 10,640 high-confidence transcripts not previously annotated, in addition to 1,479 isoforms not mapped to the current Swanson reference genome. Most of the identified lncRNAs were non-coding variants of coding transcripts. The majority of genes had multiple transcript isoforms (average ∼3 isoforms/locus). Intron retention (IR) and exon skipping (ES) accounted for 56% of alternative splicing (AS) events. Iso-Seq improved the reference genome annotation, which allowed identification of characteristic AS associated with fish growth, muscle accretion, disease resistance, stress response, and fish migration. For instance, an ES in GVIN1 gene existed in fish susceptible to bacterial cold-water disease (BCWD). Besides, under five stress conditions, there was a commonly regulated exon in prolyl 4-hydroxylase subunit alpha-2 (P4HA2) gene. The reconstructed gene models and their posttranscriptional processing in rainbow trout provide invaluable resources that could be further used for future genetics and genomics studies. Additionally, the study identified characteristic transcription events associated with economically important phenotypes, which could be applied in selective breeding.

Alternative migratory tactics in brown trout (Salmo trutta) are underpinned by divergent regulation of metabolic but not neurological genes

The occurrence of alternative morphs within populations is common, but the underlying molecular mechanisms remain poorly understood. Many animals, for example, exhibit facultative migration, where two or more alternative migratory tactics (AMTs) coexist within populations. In certain salmonid species, some individuals remain in natal rivers all their lives, while others (in particular, females) migrate to sea for a period of marine growth. Here, we performed transcriptional profiling ("RNA-seq") of the brain and liver of male and female brown trout to understand the genes and processes that differentiate between migratory and residency morphs (AMT-associated genes) and how they may differ in expression between the sexes. We found tissue-specific differences with a greater number of genes expressed differentially in the liver (n = 867 genes) compared with the brain (n = 10) between the morphs. Genes with increased expression in resident livers were enriched for Gene Ontology terms associated with metabolic processes, highlighting key molecular-genetic pathways underlying the energetic requirements associated with divergent migratory tactics. In contrast, smolt-biased genes were enriched for biological processes such as response to cytokines, suggestive of possible immune function differences between smolts and residents. Finally, we identified evidence of sex-biased gene expression for AMT-associated genes in the liver (n = 12) but not the brain. Collectively, our results provide insights into tissue-specific gene expression underlying the production of alternative life histories within and between the sexes, and point toward a key role for metabolic processes in the liver in mediating divergent physiological trajectories of migrants versus residents.

Continent-wide genomic signatures of adaptation to urbanisation in a songbird across Europe

Urbanisation is increasing worldwide, and there is now ample evidence of phenotypic changes in wild organisms in response to this novel environment. Yet, the genetic changes and genomic architecture underlying these adaptations are poorly understood. Here, we genotype 192 great tits (Parus major) from nine European cities, each paired with an adjacent rural site, to address this major knowledge gap in our understanding of wildlife urban adaptation. We find that a combination of polygenic allele frequency shifts and recurrent selective sweeps are associated with the adaptation of great tits to urban environments. While haplotypes under selection are rarely shared across urban populations, selective sweeps occur within the same genes, mostly linked to neural function and development. Collectively, we show that urban adaptation in a widespread songbird occurs through unique and shared selective sweeps in a core-set of behaviour-linked genes.

Investigating Diadromy in Fishes and Its Loss in an -Omics Era

Diadromy, the predictable movements of individuals between marine and freshwater environments, is biogeographically and phylogenetically widespread across fishes. Thus, despite the high energetic and potential fitness costs involved in moving between distinct environments, diadromy appears to be an effective life history strategy. Yet, the origin and molecular mechanisms that underpin this migratory behavior are not fully understood. In this review, we aim first to summarize what is known about diadromy in fishes; this includes the phylogenetic relationship among diadromous species, a description of the main hypotheses regarding its origin, and a discussion of the presence of non-migratory populations within diadromous species. Second, we discuss how recent research based on -omics approaches (chiefly genomics, transcriptomics, and epigenomics) is beginning to provide answers to questions on the genetic bases and origin(s) of diadromy. Finally, we suggest future directions for -omics research that can help tackle questions on the evolution of diadromy.

Genomic basis of the loss of diadromy in Galaxias maculatus: Insights from reciprocal transplant experiments

Diadromy is known for having major effects on the distribution and richness of aquatic species, and so does its loss. The loss of diadromy has led to the diversification of many species, yet research focusing on understanding its molecular basis and consequences are limited. This is particularly true for amphidromous species despite being the most abundant group of diadromous species. Galaxias maculatus, an amphidromous species and one of the most widely distributed fishes in the Southern Hemisphere, exhibits many instances of nonmigratory or resident populations. The existence of naturally replicated resident populations in Patagonia can serve as an ideal system for the study of the mechanisms that lead to the loss of the diadromy and its ecological and evolutionary consequences. Here, we studied two adjacent river systems in which resident populations are genetically differentiated yet derived from the same diadromous population. By combining a reciprocal transplant experiment with genomic data, we showed that the two resident populations followed different evolutionary pathways by exhibiting a differential response in their capacity to survive in salt water. While one resident population was able to survive salt water, the other was not. Genomic analyses provided insights into the genes that distinguished (a) migratory from nonmigratory populations; (b) populations that can vs those that cannot survive a saltwater environment; and (c) between these resident populations. This study demonstrates that the loss of diadromy can be achieved by different pathways and that environmental (selection) and random (genetic drift) forces shape this dynamic evolutionary process.

Transcriptome analysis of the brain provides insights into the regulatory mechanism for Coilia nasus migration

Background

Coilia nasus (C. nasus) is an important anadromous fish species that resides in the Yangtze River of China, and has high ecological and economical value. However, wild resources have suffered from a serious reduction in population, attributed to the over-construction of water conservancy projects, overfishing, and environmental pollution. The Ministry of Agriculture and Rural Affairs of the People’s Republic of China has issued a notice banning the commercial fishing of wild C. nasus in the Yangtze River. Wild C. nasus populations urgently need to recover. A better understanding of C. nasus migration patterns is necessary to maximize the efficiency of conservation efforts. Juvenile C. nasus experience a simultaneous effect of increasing salinity and cold stress during seaward migration, and the brain plays a comprehensive regulatory role during this process. Therefore, to explore the early seaward migration regulation mechanism of juvenile C. nasus, we performed a comparative transcriptome analysis on the brain of juvenile C. nasus under salinity and cold stress simultaneously.

Results

Relevant neurotransmitters, receptors, and regulatory proteins from three categories of regulatory pathway play synergistic regulatory roles during the migration process: neuronal signaling, the sensory system, and environmental adaptation. The significant differential expression of growth-related hormones, thyroid receptors, haptoglobin, and prolactin receptors was similar to the results of relevant research on salmonids and steelhead trout.

Conclusions

This study revealed a regulatory network that the brain of juvenile C. nasus constructs during migration, thereby providing basic knowledge on further studies could build on. This study also revealed key regulatory genes similar to salmonids and steelhead trout, thus, this study will lay a theoretical foundation for further study on migration regulation mechanism of anadromous fish species.

RNA profiling identifies novel, photoperiod-history dependent markers associated with enhanced saltwater performance in juvenile Atlantic salmon

Atlantic salmon migrate to sea following completion of a developmental process known as smolting, which establishes a seawater (SW) tolerant phenotype. Smolting is stimulated by exposure to long photoperiod or continuous light (LL) following a period of exposure to short photoperiod (SP), and this leads to major changes in gill ion exchange and osmoregulatory function. Here, we performed an RNAseq experiment to discover novel genes involved in photoperiod-dependent remodeling of the gill. This revealed a novel cohort of genes whose expression rises dramatically in fish transferred to LL following SP exposure, but not in control fish maintained continuously on LL or on SP. A follow-up experiment revealed that the SP-history dependence of LL induction of gene expression varies considerably between genes. Some genes were inducible by LL exposure after only 2 weeks exposure to SP, while others required 8 weeks prior SP exposure for maximum responsiveness to LL. Since subsequent SW growth performance is also markedly improved following 8 weeks SP exposure, these photoperiodic history-dependent genes may be useful predictive markers for full smolt development.

Plasma proteome profiling of freshwater and seawater life stages of rainbow trout (Oncorhynchus mykiss)

The sea-run phenotype of rainbow trout (Oncorhynchus mykiss), like other anadromous salmonids, present a juvenile stage fully adapted to life in freshwater known as parr. Development in freshwater is followed by the smolt stage, where preadaptations needed for seawater life are developed making fish ready to migrate to the ocean, after which event they become post-smolts. While these three life stages have been studied using a variety of approaches, proteomics has never been used for such purpose. The present study characterised the blood plasma proteome of parr, smolt and post-smolt rainbow trout using a gel electrophoresis liquid chromatography tandem mass spectrometry approach alone or in combination with low-abundant protein enrichment technology (combinatorial peptide ligand library). In total, 1,822 proteins were quantified, 17.95% of them being detected only in plasma post enrichment. Across all life stages, the most abundant proteins were ankyrin-2, DNA primase large subunit, actin, serum albumin, apolipoproteins, hemoglobin subunits, hemopexin-like proteins and complement C3. When comparing the different life stages, 17 proteins involved in mechanisms to cope with hyperosmotic stress and retinal changes, as well as the downregulation of nonessential processes in smolts, were significantly different between parr and smolt samples. On the other hand, 11 proteins related to increased growth in post-smolts, and also related to coping with hyperosmotic stress and to retinal changes, were significantly different between smolt and post-smolt samples. Overall, this study presents a series of proteins with the potential to complement current seawater-readiness assessment tests in rainbow trout, which can be measured non-lethally in an easily accessible biofluid. Furthermore, this study represents a first in-depth characterisation of the rainbow trout blood plasma proteome, having considered three life stages of the fish and used both fractionation alone or in combination with enrichment methods to increase protein detection.

A large-scale chromosomal inversion is not associated with life history development in rainbow trout from Southeast Alaska

In studying the causative mechanisms behind migration and life history, the salmonids-salmon, trout, and charr-are an exemplary taxonomic group, as life history development is known to have a strong genetic component. A double inversion located on chromosome 5 in rainbow trout (Oncorhynchus mykiss) is associated with life history development in multiple populations, but the importance of this inversion has not been thoroughly tested in conjunction with other polymorphisms in the genome. To that end, we used a high-density SNP chip to genotype 192 F1 migratory and resident rainbow trout and focused our analyses to determine whether this inversion is important in life history development in a well-studied population of rainbow trout from Southeast Alaska. We identified 4,994 and 436 SNPs-predominantly outside of the inversion region-associated with life history development in the migrant and resident familial lines, respectively. Although F1 samples showed genomic patterns consistent with the double inversion on chromosome 5 (reduced observed and expected heterozygosity and an increase in linkage disequilibrium), we found no statistical association between the inversion and life history development. Progeny produced by crossing resident trout and progeny produced by crossing migrant trout both consisted of a mix of migrant and resident individuals, irrespective of the individuals' inversion haplotype on chromosome 5. This suggests that although the inversion is present at a low frequency, it is not strongly associated with migration as it is in populations of Oncorhynchus mykiss from lower latitudes.

Habitat Partitioning and its Possible Genetic Background Between Two Sympatrically Distributed Eel Species in Taiwan

The geographical distributions of the Japanese eel (Anguilla japonica) and Giant-mottled eel (A. marmorata) overlap in many regions in East Asia and therefore suffer from interspecific competition in the same rivers. After a long period of adaptation, the Japanese eel and Giant-mottled eel may exhibit habitat partitioning in the rivers to diminish the interspecific competition between them. In this study, we conducted a field investigation in the Fengshan River in Taiwan to survey the habitat distributions of the Japanese eel and Giant-mottled eel throughout a river. Moreover, we investigated whether their habitat distributions are related to their swimming and upstream migration. Thus, the mRNA expression levels of several candidate genes that may be associated with the swimming and upstream migration of eel were examined in the glass eels of the Japanese eel and Giant-mottled eel. Field investigation indicated that the Japanese eel mainly inhabited the lower and middle reaches of the Fengshan River, but the Giant- mottled eel was distributed over the middle to upper reaches. The mRNA expression levels of fMYH, dio2, gria3, and neurod1 were higher in the Giant-mottled eel than in the Japanese eel, implying that Giant- mottled eels might have better swimming bursts and more active upstream migration than Japanese eels. These results suggest that there is a habitat partition at which these two eel species coexist in a river, and their habitat distributions may be linked to their swimming bursts and upstream migration. Determining the habitat distributions of freshwater eels is important for developing applicable plans for eel conservation and resource management.

Long-Term Conservation of Ohnologs Through Partial Tetrasomy Following Whole-Genome Duplication in Salmonidae

Whole-genome duplications (WGDs) have occurred repeatedly and broadly throughout the evolutionary history of eukaryotes. However, the effects of WGD on genome function and evolution remain unclear. The salmonid WGD that occurred approximately 88 million years ago presents an excellent opportunity for studying the effects of WGD as ∼10-15% of each salmonid genome still exhibits tetrasomic inheritance. Herein, we utilized the rainbow trout (Oncorhynchus mykiss) genome assembly and brain transcriptome data to examine the fate of gene pairs (ohnologs) following the salmonid whole-genome duplication. We find higher sequence identity between ohnologs located within known tetrasomic regions than between ohnologs found in disomic regions, and that tetrasomically inherited ohnologs showed greater similarity in patterns of gene expression and per ohnolog were lower expressed, than disomically inherited ohnologs. Enrichment testing for Gene Ontology terms identified 49 over-represented terms in tetrasomically inherited ohnologs compared to disomic ohnologs. However, why these ohnologs are retained as tetrasomic is difficult to answer. It could be that we have identified salmonid specific "dangerous duplicates", that is, genes that cannot take on new roles following WGD. Alternatively, there may be adaptive advantages for retaining genes as functional duplicates in tetrasomic regions, as presumably, movement of these genes into disomic regions would affect both their sequence identity and their gene expression patterns.

Temporal Dynamics of DNA Methylation Patterns in Response to Rearing Juvenile Steelhead (Oncorhynchus mykiss) in a Hatchery versus Simulated Stream Environment

Genetic selection is often implicated as the underlying cause of heritable phenotypic differences between hatchery and wild populations of steelhead trout (Oncorhynchus mykiss) that also differ in lifetime fitness. Developmental plasticity, which can also affect fitness, may be mediated by epigenetic mechanisms such as DNA methylation. Our previous study identified significant differences in DNA methylation between adult hatchery- and natural-origin steelhead from the same population that could not be distinguished by DNA sequence variation. In the current study, we tested whether hatchery-rearing conditions can influence patterns of DNA methylation in steelhead with known genetic backgrounds, and assessed the stability of these changes over time. Eyed-embryos from 22 families of Methow River steelhead were split across traditional hatchery tanks or a simulated stream-rearing environment for 8 months, followed by a second year in a common hatchery tank environment. Family assignments were made using a genetic parentage analysis to account for relatedness among individuals. DNA methylation patterns were examined in the liver, a relatively homogeneous organ that regulates metabolic processes and somatic growth, of juveniles at two time points: after eight months of rearing in either a tank or stream environment and after a subsequent year of rearing in a common tank environment. Further, we analyzed DNA methylation in the sperm of mature 2-year-old males from the earlier described treatments to assess the potential of environmentally-induced changes to be passed to offspring. Hepatic DNA methylation changes in response to hatchery versus stream-rearing in yearling fish were substantial, but few persisted after a second year in the tank environment. However, the early rearing environment appeared to affect how fish responded to developmental and environmental signals during the second year since novel DNA methylation differences were identified in the livers of hatchery versus stream-reared fish after a year of common tank rearing. Furthermore, we found profound differences in DNA methylation due to age, irrespective of rearing treatment. This could be due to smoltification associated changes in liver physiology after the second year of rearing. Although few rearing-treatment effects were observed in the sperm methylome, strong family effects were observed. These data suggest limited potential for intergenerational changes, but highlight the importance of understanding the effects of kinship among studied individuals in order to properly analyze and interpret DNA methylation data in natural populations. Our work is the first to study family effects and temporal dynamics of DNA methylation patterns in response to hatchery-rearing.

The genetics and epigenetics of animal migration and orientation: birds, butterflies and beyond

Migration is a complex behavioural adaptation for survival that has evolved across the animal kingdom from invertebrates to mammals. In some taxa, closely related migratory species, or even populations of the same species, exhibit different migratory phenotypes, including timing and orientation of migration. In these species, a significant proportion of the phenotypic variance in migratory traits is genetic. In others, the migratory phenotype and direction is triggered by seasonal changes in the environment, suggesting an epigenetic control of their migration. The genes and epigenetic changes underpinning migratory behaviour remain largely unknown. The revolution in (epi)genomics and functional genomic tools holds great promise to rapidly move the field of migration genetics forward. Here, we review our current understanding of the genetic and epigenetic architecture of migratory traits, focusing on two emerging models: the European blackcap and the North American monarch butterfly. We also outline a vision of how technical advances and integrative approaches could be employed to identify and functionally validate candidate genes and cis-regulatory elements on these and other migratory species across both small and broad phylogenetic scales to significantly advance the field of genetics of animal migration.

Transcriptional shifts during juvenile Coho salmon (Oncorhynchus kisutch) life stage changes in freshwater and early marine environments

There is a paucity of information on the physiological changes that occur over the course of salmon early marine migration. Here we aim to provide insight on juvenile Coho salmon (Oncorhynchus kisutch) physiology using the changes in gene expression (cGRASP 44K microarray) of four tissues (brain, gill, muscle, and liver) across the parr to smolt transition in freshwater and through the first eight months of ocean residence. We also examined transcriptome changes with body size as a covariate. The strongest shift in the transcriptome for brain, gill, and muscle occurred between summer and fall in the ocean, representing physiological changes that we speculate may be associated with migration preparation to feeding areas. Metabolic processes in the liver were positively associated with body length, generally consistent with enhanced feeding opportunities. However, a notable exception to this metabolic pattern was for spring post-smolts sampled soon after entry into the ocean, which showed a pattern of gene expression more likely associated with depressed feeding or recent fasting. Overall, this study has revealed life stages that may be the most critical developmentally (fall post-smolt) and for survival (spring post-smolt) in the early marine environment. These life stages may warrant further investigation.

Evidence of sex-bias in gene expression in the brain transcriptome of two populations of rainbow trout (Oncorhynchus mykiss) with divergent life histories

Sex-bias in gene expression is a mechanism that can generate phenotypic variance between the sexes, however, relatively little is known about how patterns of sex-bias vary during development, and how variable sex-bias is between different populations. To that end, we measured sex-bias in gene expression in the brain transcriptome of rainbow trout (Oncorhynchus mykiss) during the first two years of development. Our sampling included from the fry stage through to when O. mykiss either migrate to the ocean or remain resident and undergo sexual maturation. Samples came from two F1 lines: One from migratory steelhead trout and one from resident rainbow trout. All samples were reared in a common garden environment and RNA sequencing (RNA-seq) was used to estimate patterns of gene expression. A total of 1,716 (4.6% of total) genes showed evidence of sex-bias in gene expression in at least one time point. The majority (96.7%) of sex-biased genes were differentially expressed during the second year of development, indicating that patterns of sex-bias in expression are tied to key developmental events, such as migration and sexual maturation. Mapping of differentially expressed genes to the O. mykiss genome revealed that the X chromosome is enriched for female upregulated genes, and this may indicate a lack of dosage compensation in rainbow trout. There were many more sex-biased genes in the migratory line than the resident line suggesting differences in patterns of gene expression in the brain between populations subjected to different forces of selection. Overall, our results suggest that there is considerable variation in the extent and identity of genes exhibiting sex-bias during the first two years of life. These differentially expressed genes may be connected to developmental differences between the sexes, and/or between adopting a resident or migratory life history.