Revealing the evolutionary history and contemporary population structure of Pacific salmon in the Fraser River through genome resequencing

The Fraser River once supported massive salmon returns. However, over the last century, the largest returns have consistently been less than half of the recorded historical maximum. There is substantial interest from surrounding communities and governments to increase salmon returns for both human use and functional ecosystems. To generate resources for this endeavor, we resequenced genomes of Chinook (Oncorhynchus tshawytscha), coho (O. kisutch), and sockeye salmon (O. nerka) from the Fraser River at moderate coverage (∼16x). A total of 954 resequenced genomes were analyzed, with 681 collected specifically for this study from tissues sampled between 1997 and 2021. An additional 273 were collected from previous studies. At the species level, Chinook salmon appeared to have 1.6-2.1x more SNPs than coho or sockeye salmon, respectively. This difference may be attributable to large historical declines of coho and sockeye salmon. At the population level, three Fraser River genetic groups were identified for each species using principal component and admixture analyses, which is consistent with previous research and supports the continued use of these groups in conservation and management efforts. Environmental factors and a migration barrier were identified as major factors influencing the boundaries of these genetic groups. Additionally, 20 potentially adaptive loci were identified among the genetic groups. This information may be valuable in new management and conservation efforts. Furthermore, the resequenced genomes are an important resource for contemporary genomics research on Fraser River salmon and have been made publicly available.

Insights from Hi-C data regarding the Pacific salmon louse (Lepeophtheirus salmonis) sex chromosomes

Salmon lice, Lepeophtheirus salmonis (family Caligidae), are ectoparasites that have negatively impacted the salmon aquaculture industry and vulnerable wild salmon populations. Researchers have studied salmon lice to better understand their biology to develop effective control strategies. In this study, we updated the chromosome-level reference genome assembly of the Pacific subspecies of L. salmonis using Hi-C data. The previous version placed contigs/scaffolds using an Atlantic salmon louse genetic map. By utilizing Hi-C data from Pacific salmon lice, we were able to properly assign locations to contigs/scaffolds previously unplaced or misplaced. This resulted in a more accurate genome assembly and a more comprehensive characterization of the sex chromosome unique to females (W). We found evidence that the same ZW-ZZ mechanism is common in both Atlantic and Pacific subspecies of salmon lice using PCR assays. The W chromosome was approximately 800 kb in size, which is ∼30 times smaller than the Z chromosome (24 Mb). The W chromosome contained 61 annotated genes, including 32 protein-coding genes, 27 long noncoding RNA (lncRNA) genes, and 2 pseudogenes. Among these 61 genes, 39 genes shared homology to genes found on other chromosomes, while 20 were unique to the W chromosome. Two genes of interest on the W chromosome, prohibitin-2 and kinase suppressor of ras-2, were previously identified as potential sex-linked markers in the salmon louse. However, we prioritized the 20 unique genes on the W chromosome as sex-determining candidates. This information furthers our understanding of the biology of this ectoparasite and may help in the development of more effective management strategies.

Insights from a chum salmon (Oncorhynchus keta) genome assembly regarding whole-genome duplication and nucleotide variation influencing gene function

Chum salmon are ecologically important to Pacific Ocean ecosystems and commercially important to fisheries. To improve the genetic resources available for this species, we sequenced and assembled the genome of a male chum salmon using Oxford Nanopore read technology and the Flye genome assembly software (contig N50: ∼2 Mbp, complete BUSCOs: ∼98.1%). We also resequenced the genomes of 59 chum salmon from hatchery sources to better characterize the genome assembly and the diversity of nucleotide variants impacting phenotype variation. With genomic sequences from a doubled haploid individual, we were able to identify regions of the genome assembly that have been collapsed due to high sequence similarity between homeologous (duplicated) chromosomes. The homeologous chromosomes are relics of an ancient salmonid-specific genome duplication. These regions were enriched with genes whose functions are related to the immune system and responses to toxins. From analyzing nucleotide variant annotations of the resequenced genomes, we were also able to identify genes that have increased levels of variants thought to moderately impact gene function. Genes related to the immune system and the detection of chemical stimuli (olfaction) had increased levels of these variants based on a gene ontology enrichment analysis. The tandem organization of many of the enriched genes raises the question of why they have this organization.

Sablefish (Anoplopoma fimbria) chromosome-level genome assembly

Sablefish (Anoplopoma fimbria) are in the suborder Cottioidei, which also includes stickleback and lumpfish. This species inhabits coastal regions of the northeastern and northwestern Pacific Ocean from California to Japan. A commercial fishery for sablefish began to flourish in the 1960s, though a downward trend in stock biomass and landings has been observed since 2010. Aquaculture protocols have been developed for sablefish; eggs and sperm from wild-caught and hatchery-reared captive broodstock are used to generate offspring that reach market size in about two years. Parentage analyses show that survival in aquaculture varies among families. Growth rate and disease resistance also vary among individuals and cohorts, but the extent to which genetics and the environment contribute to this variation is unclear. The sablefish genome assembly reported here will form the foundation for SNP-based surveys designed to detect genetic markers associated with survival, growth rate, and pathogen resistance. Beyond its contribution to sablefish domestication, the sablefish genome can be a resource for the management of the wild sablefish fishery. The assembly generated in this study had a length of 653 Mbp, a scaffold N50 of 26.74 Mbp, a contig N50 of 2.57 Mbp, and contained more than 98% of the 3640 Actinopterygii core genes. We placed 620.9 Mbp (95% of the total) onto 24 chromosomes using a genetic map derived from six full-sib families and Hi-C contact data.

Population-size history inferences from the coho salmon (Oncorhynchus kisutch) genome

Coho salmon (Oncorhynchus kisutch) are a culturally and economically important species that return from multiyear ocean migrations to spawn in rivers that flow to the Northern Pacific Ocean. Southern stocks of coho salmon in Canada and the United States have significantly declined over the past quarter century, and unfortunately, conservation efforts have not reversed this trend. To assist in stock management and conservation efforts, we generated a chromosome-level genome assembly. We also resequenced the genomes of 83 coho salmon across their North American range to identify nucleotide variants, and understand the demographic histories of these salmon by modeling effective population size from genome-wide data. From demographic history modeling, we observed reductions in effective population sizes between 3750-8000 years ago for several northern sampling sites, which may correspond to bottleneck events during recolonization after glacial retreat.

Long-distance migration is a major factor driving local adaptation at continental scale in Coho salmon

Inferring the genomic basis of local adaptation is a long-standing goal of evolutionary biology. Beyond its fundamental evolutionary implications, such knowledge can guide conservation decisions for populations of conservation and management concern. Here, we investigated the genomic basis of local adaptation in the Coho salmon (Oncorhynchus kisutch) across its entire North American range. We hypothesized that extensive spatial variation in environmental conditions and the species' homing behaviour may promote the establishment of local adaptation. We genotyped 7829 individuals representing 217 sampling locations at more than 100,000 high-quality RADseq loci to investigate how recombination might affect the detection of loci putatively under selection and took advantage of the precise description of the demographic history of the species from our previous work to draw accurate population genomic inferences about local adaptation. The results indicated that genetic differentiation scans and genetic-environment association analyses were both significantly affected by variation in recombination rate as low recombination regions displayed an increased number of outliers. By taking these confounding factors into consideration, we revealed that migration distance was the primary selective factor driving local adaptation and partial parallel divergence among distant populations. Moreover, we identified several candidate single nucleotide polymorphisms associated with long-distance migration and altitude including a gene known to be involved in adaptation to altitude in other species. The evolutionary implications of our findings are discussed along with conservation applications.

Re-evaluating Coho salmon (Oncorhynchus kisutch) conservation units in Canada using genomic data

Conservation units (CUs) are important tools for supporting the implementation of standardized management practices for exploited species. Following the adoption of the Wild Salmon Policy in Canada, CUs were defined for Pacific salmon based on characteristics related to ecotype, life history and genetic variation using microsatellite markers as indirect measures of local adaptation. Genomic data sets have the potential to improve the definition of CUs by reducing variance around estimates of population genetic parameters, thereby increasing the power to detect more subtle patterns of population genetic structure and by providing an opportunity to incorporate adaptive information more directly with the identification of variants putatively under selection. We used one of the largest genomic data sets recently published for a nonmodel species, comprising 5662 individual Coho salmon (Oncorhynchus kisutch) from 149 sampling locations and a total of 24,542 high-quality SNPs obtained using genotyping-by-sequencing and mapped to the Coho salmon reference genome to (1) evaluate the current delineation of CUs for Coho in Canada and (2) compare patterns of population structure observed using neutral and outlier loci from genotype-environment association analyses to determine whether separate CUs that capture adaptive diversity are needed. Our results reflected CU boundaries on the whole, with the majority of sampling locations managed in the same CU clustering together within genetic groups. However, additional groups that are not currently represented by CUs were also uncovered. We observed considerable overlap in the genetic clusters identified using neutral or candidate loci, indicating a general congruence in patterns of genetic variation driven by local adaptation and gene flow in this species. Consequently, we suggest that the current CU boundaries for Coho salmon are largely well-suited for meeting the Canadian Wild Salmon Policy's objective of defining biologically distinct groups, but we highlight specific areas where CU boundaries may be refined.

An update of the salmon louse (Lepeophtheirus salmonis) reference genome assembly

Salmon lice have plagued the salmon farming industry and have negatively impacted salmon populations in the wild. In response, researchers have generated high density genetic maps, genome assemblies, transcriptomes, and whole-genome resequencing data to better understand this parasite. In this study, we used long-read sequencing technology to update the previous genome assemblies of Atlantic Ocean salmon lice with a more contiguous assembly and a more comprehensive gene catalog of Pacific Ocean salmon lice. We were also able to further characterize genomic features previously identified from other studies by using published resequenced genomes of 25 Atlantic and 15 Pacific salmon lice. One example was further characterizing the ZW sex chromosomes. For both the Atlantic and Pacific Ocean salmon lice subspecies, we found that the female W-chromosome is only a small fraction of the Z-chromosome and that the vast majority of the W and Z-chromosome do not contain conserved regions (i.e. pseudoautosomal regions). However, conserved orthologous protein sequences can still be identified between the W- and Z-chromosomes.

The pink salmon genome: Uncovering the genomic consequences of a two-year life cycle

Pink salmon (Oncorhynchus gorbuscha) adults are the smallest of the five Pacific salmon native to the western Pacific Ocean. Pink salmon are also the most abundant of these species and account for a large proportion of the commercial value of the salmon fishery worldwide. A two-year life history of pink salmon generates temporally isolated populations that spawn either in even-years or odd-years. To uncover the influence of this genetic isolation, reference genome assemblies were generated for each year-class and whole genome re-sequencing data was collected from salmon of both year-classes. The salmon were sampled from six Canadian rivers and one Japanese river. At multiple centromeres we identified peaks of Fst between year-classes that were millions of base-pairs long. The largest Fst peak was also associated with a million base-pair chromosomal polymorphism found in the odd-year genome near a centromere. These Fst peaks may be the result of a centromere drive or a combination of reduced recombination and genetic drift, and they could influence speciation. Other regions of the genome influenced by odd-year and even-year temporal isolation and tentatively under selection were mostly associated with genes related to immune function, organ development/maintenance, and behaviour.

The salmon louse genome: Copepod features and parasitic adaptations

Copepods encompass numerous ecological roles including parasites, detrivores and phytoplankton grazers. Nonetheless, copepod genome assemblies remain scarce. Lepeophtheirus salmonis is an economically and ecologically important ectoparasitic copepod found on salmonid fish. We present the 695.4 Mbp L. salmonis genome assembly containing ≈60% repetitive regions and 13,081 annotated protein-coding genes. The genome comprises 14 autosomes and a ZZ-ZW sex chromosome system. Assembly assessment identified 92.4% of the expected arthropod genes. Transcriptomics supported annotation and indicated a marked shift in gene expression after host attachment, including apparent downregulation of genes related to circadian rhythm coinciding with abandoning diurnal migration. The genome shows evolutionary signatures including loss of genes needed for peroxisome biogenesis, presence of numerous FNII domains, and an incomplete heme homeostasis pathway suggesting heme proteins to be obtained from the host. Despite repeated development of resistance against chemical treatments L. salmonis exhibits low numbers of many genes involved in detoxification.

Genomic basis of deep-water adaptation in Arctic Charr (Salvelinus alpinus) morphs

The post-glacial colonization of Gander Lake in Newfoundland, Canada, by Arctic Charr (Salvelinus alpinus) provides the opportunity to study the genomic basis of adaptation to extreme deep-water environments. Colonization of deep-water (>50 m) habitats often requires extensive adaptation to cope with novel environmental challenges from high hydrostatic pressure, low temperature, and low light, but the genomic mechanisms underlying evolution in these environments are rarely known. Here, we compare genomic divergence between a deep-water morph adapted to depths of up to 288 m and a larger, piscivorous pelagic morph occupying shallower depths. Using both a SNP array and resequencing of whole nuclear and mitochondrial genomes, we find clear genetic divergence (F_ST  = 0.11-0.15) between deep and shallow water morphs, despite an absence of morph divergence across the mitochondrial genome. Outlier analyses identified many diverged genomic regions containing genes enriched for processes such as gene expression and DNA repair, cardiac function, and membrane transport. Detection of putative copy number variants (CNVs) uncovered 385 genes with CNVs distinct to piscivorous morphs, and 275 genes with CNVs distinct to deep-water morphs, enriched for processes associated with synapse assembly. Demographic analyses identified evidence for recent and local morph divergence, and ongoing reductions in diversity consistent with postglacial colonization. Together, these results show that Arctic Charr morph divergence has occurred through genome-wide differentiation and elevated divergence of genes underlying multiple cellular and physiological processes, providing insight into the genomic basis of adaptation in a deep-water habitat following postglacial recolonization.

Detection of selection signatures in farmed coho salmon (Oncorhynchus kisutch) using dense genome-wide information

Animal domestication and artificial selection give rise to gradual changes at the genomic level in populations. Subsequent footprints of selection, known as selection signatures or selective sweeps, have been traced in the genomes of many animal livestock species by exploiting variation in linkage disequilibrium patterns and/or reduction of genetic diversity. Domestication of most aquatic species is recent in comparison with land animals, and salmonids are one of the most important fish species in aquaculture. Coho salmon (Oncorhynchus kisutch), cultivated primarily in Chile, has been subjected to breeding programs to improve growth, disease resistance traits, and flesh color. This study aimed to identify selection signatures that may be involved in adaptation to culture conditions and traits of productive interest. To do so, individuals of two domestic populations cultured in Chile were genotyped with 200 thousand SNPs, and analyses were conducted using iHS, XP-EHH and CLR. Several signatures of selection on different chromosomal regions were detected across both populations. Some of the identified regions under selection contained genes such anapc2, alad, chp2 and myn, which have been previously associated with body weight in Atlantic salmon, or sec24d and robo1, which have been associated with resistance to Piscirickettsia salmonis in coho salmon. Findings in our study can contribute to an integrated genome-wide map of selection signatures, to help identify the genetic mechanisms of phenotypic diversity in coho salmon.

Comparative regulomics supports pervasive selection on gene dosage following whole genome duplication

BACKGROUND

Whole genome duplication (WGD) events have played a major role in eukaryotic genome evolution, but the consequence of these extreme events in adaptive genome evolution is still not well understood. To address this knowledge gap, we used a comparative phylogenetic model and transcriptomic data from seven species to infer selection on gene expression in duplicated genes (ohnologs) following the salmonid WGD 80-100 million years ago.

RESULTS

We find rare cases of tissue-specific expression evolution but pervasive expression evolution affecting many tissues, reflecting strong selection on maintenance of genome stability following genome doubling. Ohnolog expression levels have evolved mostly asymmetrically, by diverting one ohnolog copy down a path towards lower expression and possible pseudogenization. Loss of expression in one ohnolog is significantly associated with transposable element insertions in promoters and likely driven by selection on gene dosage including selection on stoichiometric balance. We also find symmetric expression shifts, and these are associated with genes under strong evolutionary constraints such as ribosome subunit genes. This possibly reflects selection operating to achieve a gene dose reduction while avoiding accumulation of "toxic mutations". Mechanistically, ohnolog regulatory divergence is dictated by the number of bound transcription factors in promoters, with transposable elements being one likely source of novel binding sites driving tissue-specific gains in expression.

CONCLUSIONS

Our results imply pervasive adaptive expression evolution following WGD to overcome the immediate challenges posed by genome doubling and to exploit the long-term genetic opportunities for novel phenotype evolution.

Assessing the effects of genotype-by-environment interaction on epigenetic, transcriptomic, and phenotypic response in a Pacific salmon

Genotype-by-environment (GxE) interactions are non-parallel reaction norms among individuals with different genotypes in response to different environmental conditions. GxE interactions are an extension of phenotypic plasticity and consequently studying such interactions improves our ability to predict effects of different environments on phenotype as well as the fitness of genetically distinct organisms and their capacity to interact with ecosystems. Growth hormone transgenic coho salmon grow much faster than non-transgenics when raised in tank environments, but show little difference in growth when reared in nature-like streams. We used this model system to evaluate potential mechanisms underlying this growth rate GxE interaction, performing RNA-seq to measure gene transcription and whole-genome bisulfite sequencing to measure gene methylation in liver tissue. Gene ontology (GO) term analysis revealed stress as an important biological process potentially influencing growth rate GxE interactions. While few genes with transcription differences also had methylation differences, in promoter or gene regions, many genes were differentially methylated between tank and stream environments. A GO term analysis of differentially methylated genes between tank and stream environments revealed increased methylation in the stream environment of more than 95% of the differentially methylated genes, many with biological processes unrelated to liver function. The lower nutritional condition of the stream environment may cause increased negative regulation of genes less vital for liver tissue function than when fish are reared in tanks with unlimited food availability. These data show a large effect of rearing environment both on gene expression and methylation, but it is less clear that the detected epigenetic marks are responsible for the observed altered growth and physiological responses.

Sexually Dimorphic Growth Stimulation in a Strain of Growth Hormone Transgenic Coho Salmon (Oncorhynchus kisutch)

Growth hormone (GH) transgenic fish often exhibit remarkable transformations in growth rate and other phenotypes relative to wild-type. The 5750A transgenic coho salmon strain exhibits strong sexually dimorphic growth, with females possessing growth stimulation at a level typical of that seen for both sexes in other strains harbouring the same gene construct (e.g. M77), while males display a modest level of growth stimulation. GH mRNA levels were significantly higher in females than in males of the 5750A strain but equivalent in the M77 strain, indicating sex and transgene insertion locus altered transgene expression. We found that acute estradiol treatments did not influence GH expression in either strain (5750A and M77) or the transgene promoter (metallothionein-B), suggesting that estradiol level was not a significant factor influencing transgene activity. The feminization of XX and XY fish of the 5750A and M77 strains generated all-female groups and resulted in equalized growth of the two genetic sexes, suggesting that the presence of the Y chromosome was not directly capable of influencing the GH transgene-mediated growth in a physiological female conditions. These data suggest that the difference in growth rate seen between the sexes in the 5750A strain arises from non-estradiol-mediated sex influences on gene regulation at the transgene locus. This study shows how genetic factors and transgene insertion sites can influence transgene expression with significant consequent effects on phenotype.

Demographic history shaped geographical patterns of deleterious mutation load in a broadly distributed Pacific Salmon

A thorough reconstruction of historical processes is essential for a comprehensive understanding of the mechanisms shaping patterns of genetic diversity. Indeed, past and current conditions influencing effective population size have important evolutionary implications for the efficacy of selection, increased accumulation of deleterious mutations, and loss of adaptive potential. Here, we gather extensive genome-wide data that represent the extant diversity of the Coho salmon (Oncorhynchus kisutch) to address two objectives. We demonstrate that a single glacial refugium is the source of most of the present-day genetic diversity, with detectable inputs from a putative secondary micro-refugium. We found statistical support for a scenario whereby ancestral populations located south of the ice sheets expanded recently, swamping out most of the diversity from other putative micro-refugia. Demographic inferences revealed that genetic diversity was also affected by linked selection in large parts of the genome. Moreover, we demonstrate that the recent demographic history of this species generated regional differences in the load of deleterious mutations among populations, a finding that mirrors recent results from human populations and provides increased support for models of expansion load. We propose that insights from these historical inferences should be better integrated in conservation planning of wild organisms, which currently focuses largely on neutral genetic diversity and local adaptation, with the role of potentially maladaptive variation being generally ignored.

Reconstruction of a species’ past demographic history from genetic data can highlight historical factors that have shaped the distribution of genetic diversity along its genome and its geographic range. Here, we combine genotyping-by-sequencing with demographic modelling to address these issues in the Coho salmon, a Pacific salmon of conservation concern in some parts of its range, notably in the south. Our demographic reconstructions reveal a linear decrease in genetic diversity toward the north of the species range, supporting the hypothesis of a northern route of postglacial recolonization from a single major southern refugium. As predicted by theory, we also observed a higher proportion of deleterious mutations in the most distant populations from this refugium. Beyond this general pattern, among-site variation in the proportion of deleterious mutations is consistent with different local trends in effective population sizes. Our results highlight the potential importance of understanding historical factors that have shaped geographic patterns of the distribution of deleterious mutations in order to implement effective management programs for the conservation of wild populations. Such fundamental knowledge of human historical demography is now having major impacts on health sciences, and we argue it is time to integrate such approaches in conservation science as well.

Resolving fine-scale population structure and fishery exploitation using sequenced microsatellites in a northern fish

The resiliency of populations and species to environmental change is dependent on the maintenance of genetic diversity, and as such, quantifying diversity is central to combating ongoing widespread reductions in biodiversity. With the advent of next-generation sequencing, several methods now exist for resolving fine-scale population structure, but the comparative performance of these methods for genetic assignment has rarely been tested. Here, we evaluate the performance of sequenced microsatellites and a single nucleotide polymorphism (SNP) array to resolve fine-scale population structure in a critically important salmonid in north eastern Canada, Arctic Charr (Salvelinus alpinus). We also assess the utility of sequenced microsatellites for fisheries applications by quantifying the spatial scales of movement and exploitation through genetic assignment of fishery samples to rivers of origin and comparing these results with a 29-year tagging dataset. Self-assignment and simulation-based analyses of 111 genome-wide microsatellite loci and 500 informative SNPs from 28 populations of Arctic Charr in north-eastern Canada identified largely river-specific genetic structure. Despite large differences (~4X) in the number of loci surveyed between panels, mean self-assignment accuracy was similar with the microsatellite loci and the SNP panel (>90%). Subsequent analysis of 996 fishery-collected samples using the microsatellite panel revealed that larger rivers contribute greater numbers of individuals to the fishery and that coastal fisheries largely exploit individuals originating from nearby rivers, corroborating results from traditional tagging experiments. Our results demonstrate the efficacy of sequence-based microsatellite genotyping to advance understanding of fine-scale population structure and harvest composition in northern and understudied species.

Parallelism in eco-morphology and gene expression despite variable evolutionary and genomic backgrounds in a Holarctic fish

Understanding the extent to which ecological divergence is repeatable is essential for predicting responses of biodiversity to environmental change. Here we test the predictability of evolution, from genotype to phenotype, by studying parallel evolution in a salmonid fish, Arctic charr (Salvelinus alpinus), across eleven replicate sympatric ecotype pairs (benthivorous-planktivorous and planktivorous-piscivorous) and two evolutionary lineages. We found considerable variability in eco-morphological divergence, with several traits related to foraging (eye diameter, pectoral fin length) being highly parallel even across lineages. This suggests repeated and predictable adaptation to environment. Consistent with ancestral genetic variation, hundreds of loci were associated with ecotype divergence within lineages of which eight were shared across lineages. This shared genetic variation was maintained despite variation in evolutionary histories, ranging from postglacial divergence in sympatry (ca. 10-15kya) to pre-glacial divergence (ca. 20-40kya) with postglacial secondary contact. Transcriptome-wide gene expression (44,102 genes) was highly parallel across replicates, involved biological processes characteristic of ecotype morphology and physiology, and revealed parallelism at the level of regulatory networks. This expression divergence was not only plastic but in part genetically controlled by parallel cis-eQTL. Lastly, we found that the magnitude of phenotypic divergence was largely correlated with the genetic differentiation and gene expression divergence. In contrast, the direction of phenotypic change was mostly determined by the interplay of adaptive genetic variation, gene expression, and ecosystem size. Ecosystem size further explained variation in putatively adaptive, ecotype-associated genomic patterns within and across lineages, highlighting the role of environmental variation and stochasticity in parallel evolution. Together, our findings demonstrate the parallel evolution of eco-morphology and gene expression within and across evolutionary lineages, which is controlled by the interplay of environmental stochasticity and evolutionary contingencies, largely overcoming variable evolutionary histories and genomic backgrounds.

Publisher Correction: Sex-dependent dominance maintains migration supergene in rainbow trout

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Sex-dependent dominance maintains migration supergene in rainbow trout

Males and females often differ in their fitness optima for shared traits that have a shared genetic basis, leading to sexual conflict. Morphologically differentiated sex chromosomes can resolve this conflict and protect sexually antagonistic variation, but they accumulate deleterious mutations. However, how sexual conflict is resolved in species that lack differentiated sex chromosomes is largely unknown. Here we present a chromosome-anchored genome assembly for rainbow trout (Oncorhynchus mykiss) and characterize a 55-Mb double-inversion supergene that mediates sex-specific migratory tendency through sex-dependent dominance reversal, an alternative mechanism for resolving sexual conflict. The double inversion contains key photosensory, circadian rhythm, adiposity and sex-related genes and displays a latitudinal frequency cline, indicating environmentally dependent selection. Our results show sex-dependent dominance reversal across a large autosomal supergene, a mechanism for sexual conflict resolution capable of protecting sexually antagonistic variation while avoiding the homozygous lethality and deleterious mutations associated with typical heteromorphic sex chromosomes.

Standardized IMGT® Nomenclature of Salmonidae IGH Genes, the Paradigm of Atlantic Salmon and Rainbow Trout: From Genomics to Repertoires

In teleost fish as in mammals, humoral adaptive immunity is based on B lymphocytes expressing highly diverse immunoglobulins (IG). During B cell differentiation, IG loci are subjected to genomic rearrangements of V, D, and J genes, producing a unique antigen receptor expressed on the surface of each lymphocyte. During the course of an immune response to infections or immunizations, B cell clones specific of epitopes from the immunogen are expanded and activated, leading to production of specific antibodies. Among teleost fish, salmonids comprise key species for aquaculture. Rainbow trout (Oncorhynchus mykiss) and Atlantic salmon (Salmo salar) are especially important from a commercial point of view and have emerged as critical models for fish immunology. The growing interest to capture accurate and comprehensive antibody responses against common pathogens and vaccines has resulted in recent efforts to sequence the IG repertoire in these species. In this context, a unified and standardized nomenclature of salmonid IG heavy chain (IGH) genes is urgently required, to improve accuracy of annotation of adaptive immune receptor repertoire dataset generated by high-throughput sequencing (AIRRseq) and facilitate comparisons between studies and species. Interestingly, the assembly of salmonids IGH genomic sequences is challenging due to the presence of two large size duplicated IGH loci and high numbers of IG genes and pseudogenes. We used data available for Atlantic salmon to establish an IMGT standardized nomenclature of IGH genes in this species and then applied the IMGT rules to the rainbow trout IGH loci to set up a nomenclature, which takes into account the specificities of Salmonid loci. This unique, consistent nomenclature for Salmonid IGH genes was then used to construct IMGT sequence reference directories allowing accurate annotation of AIRRseq data. The complex issues raised by the genetic diversity of salmon and trout strains are discussed in the context of IG repertoire annotation.

Carotenoid pigmentation in salmon: variation in expression at BCO2-l locus controls a key fitness trait affecting red coloration

Carotenoids are primarily responsible for the characteristic red flesh coloration of salmon. Flesh coloration is an economically and evolutionarily significant trait that varies inter- and intra-specifically, yet the underlying genetic mechanism is unknown. Chinook salmon (Oncorhynchus tshawytscha) represents an ideal system to study carotenoid variation as, unlike other salmonids, they exhibit extreme differences in carotenoid utilization due to genetic polymorphisms. Here, we crossed populations of Chinook salmon with fixed differences in flesh coloration (red versus white) for a genome-wide association study to identify loci associated with pigmentation. Here, the beta-carotene oxygenase 2-like (BCO2-l) gene was significantly associated with flesh colour, with the most significant single nucleotide polymorphism explaining 66% of the variation in colour. BCO2 gene disruption is linked to carotenoid accumulation in other taxa, therefore we hypothesize that an ancestral mutation partially disrupting BCO2-l activity (i.e. hypomorphic mutation) allowed the deposition and accumulation of carotenoids within Salmonidae. Indeed, we found elevated transcript levels of BCO2-l in white Chinook salmon relative to red. The long-standing mystery of why salmon are red, while no other fishes are, is thus probably explained by a hypomorphic mutation in the proto-salmonid at the time of divergence of red-fleshed salmonid genera (approx. 30 Ma).