Identification of Multiple QTL Hotspots in Sockeye Salmon (Oncorhynchus nerka) Using Genotyping-by-Sequencing and a Dense Linkage Map

The genomic basis of reproductive and migratory behaviour in a polymorphic salmonid

Recent ecotypic differentiation provides unique opportunities to investigate the genomic basis and architecture of local adaptation, while offering insights into how species form and persist. Sockeye salmon (Oncorhynchus nerka) exhibit migratory and resident ("kokanee") ecotypes, which are further distinguished into shore-spawning and stream-spawning reproductive ecotypes. Here, we analysed 36 sockeye (stream-spawning) and kokanee (stream- and shore-spawning) genomes from a system where they co-occur and have recent common ancestry (Okanagan Lake/River in British Columbia, Canada) to investigate the genomic basis of reproductive and migratory behaviour. Examination of the genomic landscape of differentiation, differences in allele frequencies and genotype-phenotype associations revealed three main blocks of sequence differentiation on chromosomes 7, 12 and 20, associated with migratory behaviour, spawning location and spawning timing. Structural variants identified in these same areas suggest they could contribute to ecotypic differentiation directly as causal variants or via maintenance of their genomic architecture through recombination suppression mechanisms. Genes in these regions were related to spatial memory and swimming endurance (SYNGAP, TPM3), as well as eye and brain development (including SIX6), potentially associated with differences in migratory behaviour and visual habitats across spawning locations, respectively. Additional genes (GREB1L, ROCK1) identified here have been associated with timing of migration in other salmonids and could explain variation in timing of O. nerka spawning. Together, these results based on the joint analysis of sequence and structural variation represent a significant advance in our understanding of the genomic landscape of ecotypic differentiation at different stages in the speciation continuum.

Evolution of thermal physiology alters the projected range of threespine stickleback under climate change

Species distribution models (SDMs) are widely used to predict range shifts but could be unreliable under climate change scenarios because they do not account for evolution. The thermal physiology of a species is a key determinant of its range and thus incorporating thermal trait evolution into SDMs might be expected to alter projected ranges. We identified a genetic basis for physiological and behavioural traits that evolve in response to temperature change in natural populations of threespine stickleback (Gasterosteus aculeatus). Using these data, we created geographical range projections using a mechanistic niche area approach under two climate change scenarios. Under both scenarios, trait data were either static ("no evolution" models), allowed to evolve at observed evolutionary rates ("evolution" models) or allowed to evolve at a rate of evolution scaled by the trait variance that is explained by quantitative trait loci (QTL; "scaled evolution" models). We show that incorporating these traits and their evolution substantially altered the projected ranges for a widespread panmictic marine population, with over 7-fold increases in area under climate change projections when traits are allowed to evolve. Evolution-informed SDMs should improve the precision of forecasting range dynamics under climate change, and aid in their application to management and the protection of biodiversity.

Constraint and divergence in the evolution of male and female recombination rates in fishes

Recombination is a fundamental feature of sexual reproduction across eukaryotes, yet recombination rates are highly variable both within and between species. In particular, sex differences in recombination rate between males and females (heterochiasmy) are more often the rule than the exception, but despite the prevalence of heterochiasmy the ultimate causes of global patterns of heterochiasmy remain unclear. Here, we assemble a comprehensive dataset of sex-specific recombination rate estimates for 61 fish species, and combine this with information on sex determination, fertilization mode, and sexual dimorphism to test competing theories for the causes and evolution of heterochiasmy. We find that sex differences in recombination rate are evolutionary labile, with frequent shifts in the direction and magnitude of heterochiasmy. This rapid turnover does not appear to be driven by simple neutral processes and is inconsistent with nonadaptive explanations for heterochiasmy, including biological sex differences in meiosis. Although patterns of heterochiasmy across the phylogeny indicate a potential role for adaptive processes, we are unable to directly link variation in heterochiasmy with proxies for sexual selection or sexual conflict across species, indicating that these effects-if present-are either subtle or complex. Finally, we show evidence for correlated rates of recombination rate evolution between males and females, indicating the potential for genetic constraints and sexual conflict over the recombination landscape.

Dissection of quantitative trait loci in the Lachancea waltii yeast species highlights major hotspots

Dissecting the genetic basis of complex trait remains a real challenge. The budding yeast Saccharomyces cerevisiae has become a model organism for studying quantitative traits, successfully increasing our knowledge in many aspects. However, the exploration of the genotype-phenotype relationship in non-model yeast species could provide a deeper insight into the genetic basis of complex traits. Here, we have studied this relationship in the Lachancea waltii species which diverged from the S. cerevisiae lineage prior to the whole-genome duplication. By performing linkage mapping analyses in this species, we identified 86 quantitative trait loci (QTL) impacting the growth in a large number of conditions. The distribution of these loci across the genome has revealed two major QTL hotspots. A first hotspot corresponds to a general growth QTL, impacting a wide range of conditions. By contrast, the second hotspot highlighted a trade-off with a disadvantageous allele for drug-free conditions which proved to be advantageous in the presence of several drugs. Finally, a comparison of the detected QTL in L. waltii with those which had been previously identified for the same trait in a closely related species, Lachancea kluyveri was performed. This analysis clearly showed the absence of shared QTL across these species. Altogether, our results represent a first step toward the exploration of the genetic architecture of quantitative trait across different yeast species.

The sockeye salmon genome, transcriptome, and analyses identifying population defining regions of the genome

Sockeye salmon (Oncorhynchus nerka) is a commercially and culturally important species to the people that live along the northern Pacific Ocean coast. There are two main sockeye salmon ecotypes—the ocean-going (anadromous) ecotype and the fresh-water ecotype known as kokanee. The goal of this study was to better understand the population structure of sockeye salmon and identify possible genomic differences among populations and between the two ecotypes. In pursuit of this goal, we generated the first reference sockeye salmon genome assembly and an RNA-seq transcriptome data set to better annotate features of the assembly. Resequenced whole-genomes of 140 sockeye salmon and kokanee were analyzed to understand population structure and identify genomic differences between ecotypes. Three distinct geographic and genetic groups were identified from analyses of the resequencing data. Nucleotide variants in an immunoglobulin heavy chain variable gene cluster on chromosome 26 were found to differentiate the northwestern group from the southern and upper Columbia River groups. Several candidate genes were found to be associated with the kokanee ecotype. Many of these genes were related to ammonia tolerance or vision. Finally, the sex chromosomes of this species were better characterized, and an alternative sex-determination mechanism was identified in a subset of upper Columbia River kokanee.

A high-density SNP-based genetic map and several economic traits-related loci in Pelteobagrus vachelli

BACKGROUND

A high-density genetic linkage map is essential for QTL fine mapping, comparative genome analysis, identification of candidate genes and marker-assisted selection in aquaculture species. Pelteobagrus vachelli is a very popular commercial species in Asia. However, some specific characters hindered achievement of the traditional selective breeding based on phenotypes, such as lack of large-scale genomic resource and short of markers tightly associated with growth, sex determination and hypoxia tolerance related traits.

RESULTS

By making use of 5059 ddRAD markers in P. vachelli, a high-resolution genetic linkage map was successfully constructed. The map' length was 4047.01 cM by using an interval of 0.11 cm, which is an average marker standard. Comparative genome mapping revealed that a high proportion (83.2%) of markers with a one-to-one correspondence were observed between P. vachelli and P. fulvidraco. Based on the genetic map, 8 significant genome-wide QTLs for 4 weight, 1 body proportion, 2 sex determination, and 1 hypoxia tolerance related traits were detected on 4 LGs. Some SNPs from these significant genome-wide QTLs were observably associated with these phenotypic traits in other individuals by Kompetitive Allele Specific PCR. In addition, two candidate genes for weight, Sipa1 and HSD11B2, were differentially expressed between fast-, medium- and slow-growing P. vachelli. Sema7a, associated with hypoxia tolerance, was induced after hypoxia exposure and reoxygenation.

CONCLUSIONS

We mapped a set of suggestive and significant QTLs as well as candidate genes for 12 growth, 1 sex determination and 1 hypoxia tolerance related traits based on a high-density genetic linkage map by making use of SNP markers for P. fulvidraco. Our results have offered a valuable method about the much more efficient production of all-male, fast growth and hypoxia tolerance P. vachelli for the aquaculture industry.

Power of a dual-use SNP panel for pedigree reconstruction and population assignment

The use of high-throughput, low-density sequencing approaches has dramatically increased in recent years in studies of eco-evolutionary processes in wild populations and domestication in commercial aquaculture. Most of these studies focus on identifying panels of SNP loci for a single downstream application, whereas there have been few studies examining the trade-offs for selecting panels of markers for use in multiple applications. Here, we detail the use of a bioinformatic workflow for the development of a dual-purpose SNP panel for parentage and population assignment, which included identifying putative SNP loci, filtering for the most informative loci for the two tasks, designing effective multiplex PCR primers, optimizing the SNP panel for performance, and performing quality control steps for downstream applications. We applied this workflow to two adjacent Alaskan Sockeye Salmon populations and identified a GTseq panel of 142 SNP loci for parentage and 35 SNP loci for population assignment. Only 50-75 panel loci were necessary for >95% accurate parentage, whereas population assignment success, with all 172 panel loci, ranged from 93.9% to 96.2%. Finally, we discuss the trade-offs and complexities of the decision-making process that drives SNP panel development, optimization, and testing.

Comparative Genomic Analyses and a Novel Linkage Map for Cisco (Coregonus artedi) Provide Insights into Chromosomal Evolution and Rediploidization Across Salmonids

Whole-genome duplication (WGD) is hypothesized to be an important evolutionary mechanism that can facilitate adaptation and speciation. Genomes that exist in states of both diploidy and residual tetraploidy are of particular interest, as mechanisms that maintain the ploidy mosaic after WGD may provide important insights into evolutionary processes. The Salmonidae family exhibits residual tetraploidy, and this, combined with the evolutionary diversity formed after an ancestral autotetraploidization event, makes this group a useful study system. In this study, we generate a novel linkage map for cisco (Coregonus artedi), an economically and culturally important fish in North America and a member of the subfamily Coregoninae, which previously lacked a high-density haploid linkage map. We also conduct comparative genomic analyses to refine our understanding of chromosomal fusion/fission history across salmonids. To facilitate this comparative approach, we use the naming strategy of protokaryotype identifiers (PKs) to associate duplicated chromosomes to their putative ancestral state. The female linkage map for cisco contains 20,292 loci, 3,225 of which are likely within residually tetraploid regions. Comparative genomic analyses revealed that patterns of residual tetrasomy are generally conserved across species, although interspecific variation persists. To determine the broad-scale retention of residual tetrasomy across the salmonids, we analyze sequence similarity of currently available genomes and find evidence of residual tetrasomy in seven of the eight chromosomes that have been previously hypothesized to show this pattern. This interspecific variation in extent of rediploidization may have important implications for understanding salmonid evolutionary histories and informing future conservation efforts.

Mapping of Adaptive Traits Enabled by a High-Density Linkage Map for Lake Trout

Understanding the genomic basis of adaptative intraspecific phenotypic variation is a central goal in conservation genetics and evolutionary biology. Lake trout (Salvelinus namaycush) are an excellent species for addressing the genetic basis for adaptive variation because they express a striking degree of ecophenotypic variation across their range; however, necessary genomic resources are lacking. Here we utilize recently-developed analytical methods and sequencing technologies to (1) construct a high-density linkage and centromere map for lake trout, (2) identify loci underlying variation in traits that differentiate lake trout ecophenotypes and populations, (3) determine the location of the lake trout sex determination locus, and (4) identify chromosomal homologies between lake trout and other salmonids of varying divergence. The resulting linkage map contains 15,740 single nucleotide polymorphisms (SNPs) mapped to 42 linkage groups, likely representing the 42 lake trout chromosomes. Female and male linkage group lengths ranged from 43.07 to 134.64 centimorgans, and 1.97 to 92.87 centimorgans, respectively. We improved the map by determining coordinates for 41 of 42 centromeres, resulting in a map with 8 metacentric chromosomes and 34 acrocentric or telocentric chromosomes. We use the map to localize the sex determination locus and multiple quantitative trait loci (QTL) associated with intraspecific phenotypic divergence including traits related to growth and body condition, patterns of skin pigmentation, and two composite geomorphometric variables quantifying body shape. Two QTL for the presence of vermiculations and spots mapped with high certainty to an arm of linkage group Sna3, growth related traits mapped to two QTL on linkage groups Sna1 and Sna12, and putative body shape QTL were detected on six separate linkage groups. The sex determination locus was mapped to Sna4 with high confidence. Synteny analysis revealed that lake trout and congener Arctic char (Salvelinus alpinus) are likely differentiated by three or four chromosomal fissions, possibly one chromosomal fusion, and 6 or more large inversions. Combining centromere mapping information with putative inversion coordinates revealed that the majority of detected inversions differentiating lake trout from other salmonids are pericentric and located on acrocentric and telocentric linkage groups. Our results suggest that speciation and adaptive divergence within the genus Salvelinus may have been associated with multiple pericentric inversions occurring primarily on acrocentric and telocentric chromosomes. The linkage map presented here will be a critical resource for advancing conservation oriented genomic research on lake trout and exploring chromosomal evolution within and between salmonid species.

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.

Mapping of Quantitative Trait Loci for Growth and Carcass-Related Traits in Chickens Using a Restriction-Site Associated DNA Sequencing Method

In the present study, quantitative trait loci (QTLs) analysis was performed to identify the chromosomal positions of growth and carcass-related trait QTLs using 319 F₂ chickens obtained from intercrosses of an Oh-Shamo male and four White Plymouth Rock females. Body weight was measured weekly until the birds were 7 weeks old. Carcass-related traits were also measured at this timepoint. A genetic linkage map was constructed using 545 single nucleotide polymorphism (SNP) markers that were developed using a restriction-site associated DNA sequencing method. The linkage map included the 23 autosomes and the Z chromosome. Using simple interval QTL mapping, we were able to identify 10 significant and suggestive main-effect QTLs for growth and carcass-related traits present on chromosomes 1, 2, 3, 5, 8, 19, 24, and Z. These loci explained 5.60–16.52% of the phenotypic variances. The chromosomal positions of the 10 QTLs overlapped with those of previously reported QTLs, whereas the targeted traits varied. Our QTLs will aid future breeding programs in improving growth and meat yield of chickens (e.g., via marker-assisted selection), particularly in the Japanese brand chicken industry.

Applications of genotyping by sequencing in aquaculture breeding and genetics

Selective breeding is increasingly recognized as a key component of sustainable production of aquaculture species. The uptake of genomic technology in aquaculture breeding has traditionally lagged behind terrestrial farmed animals. However, the rapid development and application of sequencing technologies has allowed aquaculture to narrow the gap, leading to substantial genomic resources for all major aquaculture species. While high-density single-nucleotide polymorphism (SNP) arrays for some species have been developed recently, direct genotyping by sequencing (GBS) techniques have underpinned many of the advances in aquaculture genetics and breeding to date. In particular, restriction-site associated DNA sequencing (RAD-Seq) and subsequent variations have been extensively applied to generate population-level SNP genotype data. These GBS techniques are not dependent on prior genomic information such as a reference genome assembly for the species of interest. As such, they have been widely utilized by researchers and companies focussing on nonmodel aquaculture species with relatively small research communities. Applications of RAD-Seq techniques have included generation of genetic linkage maps, performing genome-wide association studies, improvements of reference genome assemblies and, more recently, genomic selection for traits of interest to aquaculture like growth, sex determination or disease resistance. In this review, we briefly discuss the history of GBS, the nuances of the various GBS techniques, bioinformatics approaches and application of these techniques to various aquaculture species.

Using Linkage Maps as a Tool To Determine Patterns of Chromosome Synteny in the Genus Salvelinus

Next generation sequencing techniques have revolutionized the collection of genome and transcriptome data from nonmodel organisms. This manuscript details the application of restriction site-associated DNA sequencing (RADseq) to generate a marker-dense genetic map for Brook Trout (Salvelinus fontinalis). The consensus map was constructed from three full-sib families totaling 176 F₁ individuals. The map consisted of 42 linkage groups with a total female map size of 2502.5 cM, and a total male map size of 1863.8 cM. Synteny was confirmed with Atlantic Salmon for 38 linkage groups, with Rainbow Trout for 37 linkage groups, Arctic Char for 36 linkage groups, and with a previously published Brook Trout linkage map for 39 linkage groups. Comparative mapping confirmed the presence of 8 metacentric and 34 acrocentric chromosomes in Brook Trout. Six metacentric chromosomes seem to be conserved with Arctic Char suggesting there have been at least two species-specific fusion and fission events within the genus Salvelinus. In addition, the sex marker (sdY; sexually dimorphic on the Y chromosome) was mapped to Brook Trout BC35, which is homologous with Atlantic Salmon Ssa09qa, Rainbow Trout Omy25, and Arctic Char AC04q. Ultimately, this linkage map will be a useful resource for studies on the genome organization of Salvelinus, and facilitates comparisons of the Salvelinus genome with Salmo and Oncorhynchus.

Genomic Changes Associated with Reproductive and Migratory Ecotypes in Sockeye Salmon (Oncorhynchus nerka)

Mechanisms underlying adaptive evolution can best be explored using paired populations displaying similar phenotypic divergence, illuminating the genomic changes associated with specific life history traits. Here, we used paired migratory [anadromous vs. resident (kokanee)] and reproductive [shore- vs. stream-spawning] ecotypes of sockeye salmon (Oncorhynchus nerka) sampled from seven lakes and two rivers spanning three catchments (Columbia, Fraser, and Skeena) in British Columbia, Canada to investigate the patterns and processes underlying their divergence. Restriction-site associated DNA sequencing was used to genotype this sampling at 7,347 single nucleotide polymorphisms, 334 of which were identified as outlier loci and candidates for divergent selection within at least one ecotype comparison. Sixty-eight of these outliers were present in two or more comparisons, with 33 detected across multiple catchments. Of particular note, one locus was detected as the most significant outlier between shore and stream-spawning ecotypes in multiple comparisons and across catchments (Columbia, Fraser, and Snake). We also detected several genomic islands of divergence, some shared among comparisons, potentially showing linked signals of differential selection. The single nucleotide polymorphisms and genomic regions identified in our study offer a range of mechanistic hypotheses associated with the genetic basis of O. nerka life history variation and provide novel tools for informing fisheries management.

Significant Synteny and Colocalization of Ecologically Relevant Quantitative Trait Loci Within and Across Species of Salmonid Fishes

The organization of functional regions within genomes has important implications for evolutionary potential. Considerable research effort has gone toward identifying the genomic basis of phenotypic traits of interest through quantitative trait loci (QTL) analyses. Less research has assessed the arrangement of QTL in the genome within and across species. To investigate the distribution, extent of colocalization, and the synteny of QTL for ecologically relevant traits, we used a comparative genomic mapping approach within and across a range of salmonid species. We compiled 943 QTL from all available species [lake whitefish (Coregonus clupeaformis), coho salmon (Oncorhynchus kisutch), rainbow trout (O. mykiss), Chinook salmon (O. tshawytscha), Atlantic salmon (Salmo salar), and Arctic charr (Salvelinus alpinus)]. We developed a novel analytical framework for mapping and testing the distribution of these QTL. We found no correlation between QTL density and gene density at the chromosome level but did at the fine-scale. Two chromosomes were significantly enriched for QTL. We found multiple synteny blocks for morphological, life history, and physiological traits across species, but only morphology and physiology had significantly more than expected. Two or three pairs of traits were significantly colocalized in three species (lake whitefish, coho salmon, and rainbow trout). Colocalization and fine-scale synteny suggest genetic linkage between traits within species and a conserved genetic basis across species. However, this pattern was weak overall, with colocalization and synteny being relatively rare. These findings advance our understanding of the role of genomic organization in the renowned ecological and phenotypic variability of salmonid fishes.

Sex Chromosome Evolution, Heterochiasmy, and Physiological QTL in the Salmonid Brook Charr Salvelinus fontinalis

Whole-genome duplication (WGD) can have large impacts on genome evolution, and much remains unknown about these impacts. This includes the mechanisms of coping with a duplicated sex determination system and whether this has an impact on increasing the diversity of sex determination mechanisms. Other impacts include sexual conflict, where alleles having different optimums in each sex can result in sequestration of genes into nonrecombining sex chromosomes. Sex chromosome development itself may involve sex-specific recombination rate (i.e., heterochiasmy), which is also poorly understood. The family Salmonidae is a model system for these phenomena, having undergone autotetraploidization and subsequent rediploidization in most of the genome at the base of the lineage. The salmonid master sex determining gene is known, and many species have nonhomologous sex chromosomes, putatively due to transposition of this gene. In this study, we identify the sex chromosome of Brook Charr Salvelinus fontinalis and compare sex chromosome identities across the lineage (eight species and four genera). Although nonhomology is frequent, homologous sex chromosomes and other consistencies are present in distantly related species, indicating probable convergence on specific sex and neo-sex chromosomes. We also characterize strong heterochiasmy with 2.7-fold more crossovers in maternal than paternal haplotypes with paternal crossovers biased to chromosome ends. When considering only rediploidized chromosomes, the overall heterochiasmy trend remains, although with only 1.9-fold more recombination in the female than the male. Y chromosome crossovers are restricted to a single end of the chromosome, and this chromosome contains a large interspecific inversion, although its status between males and females remains unknown. Finally, we identify quantitative trait loci (QTL) for 21 unique growth, reproductive, and stress-related phenotypes to improve knowledge of the genetic architecture of these traits important to aquaculture and evolution.

Congruent population structure across paralogous and nonparalogous loci in Salish Sea chum salmon (Oncorhynchus keta)

Whole-genome duplications are major evolutionary events with a lasting impact on genome structure. Duplication events complicate genetic analyses as paralogous sequences are difficult to distinguish; consequently, paralogs are often excluded from studies. The effects of an ancient whole-genome duplication (approximately 88 MYA) are still evident in salmonids through the persistence of numerous paralogous gene sequences and partial tetrasomic inheritance. We use restriction site-associated DNA sequencing on 10 collections of chum salmon from the Salish Sea in the USA and Canada to investigate genetic diversity and population structure in both tetrasomic and rediploidized regions of the genome. We use a pedigree and high-density linkage map to identify paralogous loci and to investigate genetic variation across the genome. By applying multivariate statistical methods, we show that it is possible to characterize paralogous loci and that they display similar patterns of population structure as the diploidized portion of the genome. We find genetic associations with the adaptively important trait of run-timing in both sets of loci. By including paralogous loci in genome scans, we can observe evolutionary signals in genomic regions that have routinely been excluded from population genetic studies in other polyploid-derived species.

A Dense Brown Trout (Salmo trutta) Linkage Map Reveals Recent Chromosomal Rearrangements in the Salmo Genus and the Impact of Selection on Linked Neutral Diversity

High-density linkage maps are valuable tools for conservation and eco-evolutionary issues. In salmonids, a complex rediploidization process consecutive to an ancient whole genome duplication event makes linkage maps of prime importance for investigating the evolutionary history of chromosome rearrangements. Here, we developed a high-density consensus linkage map for the brown trout (Salmo trutta), a socioeconomically important species heavily impacted by human activities. A total of 3977 ddRAD markers were mapped and ordered in 40 linkage groups using sex- and lineage-averaged recombination distances obtained from two family crosses. Performing map comparison between S. trutta and its sister species, S. salar, revealed extensive chromosomal rearrangements. Strikingly, all of the fusion and fission events that occurred after the S. salar/S. trutta speciation happened in the Atlantic salmon branch, whereas the brown trout remained closer to the ancestral chromosome structure. Using the strongly conserved synteny within chromosome arms, we aligned the brown trout linkage map to the Atlantic salmon genome sequence to estimate the local recombination rate in S. trutta at 3721 loci. A significant positive correlation between recombination rate and within-population nucleotide diversity (π) was found, indicating that selection constrains variation at linked neutral sites in brown trout. This new high-density linkage map provides a useful genomic resource for future aquaculture, conservation, and eco-evolutionary studies in brown trout.

A SNP Based Linkage Map of the Arctic Charr (Salvelinus alpinus) Genome Provides Insights into the Diploidization Process After Whole Genome Duplication

Diploidization, which follows whole genome duplication events, does not occur evenly across the genome. In salmonid fishes, certain pairs of homeologous chromosomes preserve tetraploid loci in higher frequencies toward the telomeres due to residual tetrasomic inheritance. Research suggests this occurs only in homeologous pairs where one chromosome arm has undergone a fusion event. We present a linkage map for Arctic charr (Salvelinus alpinus), a salmonid species with relatively fewer chromosome fusions. Genotype by sequencing identified 19,418 SNPs, and a linkage map consisting of 4508 markers was constructed from a subset of high quality SNPs and microsatellite markers that were used to anchor the new map to previous versions. Both male- and female-specific linkage maps contained the expected number of 39 linkage groups. The chromosome type associated with each linkage group was determined, and 10 stable metacentric chromosomes were identified, along with a chromosome polymorphism involving the sex chromosome AC04. Two instances of a weak form of pseudolinkage were detected in the telomeric regions of homeologous chromosome arms in both female and male linkage maps. Chromosome arm homologies within the Atlantic salmon (Salmo salar) and rainbow trout (Oncorhynchus mykiss) genomes were determined. Paralogous sequence variants (PSVs) were identified, and their comparative BLASTn hit locations showed that duplicate markers exist in higher numbers on seven pairs of homeologous arms, previously identified as preserving tetrasomy in salmonid species. Homeologous arm pairs where neither arm has been part of a fusion event in Arctic charr had fewer PSVs, suggesting faster diploidization rates in these regions.

Genomic islands of divergence linked to ecotypic variation in sockeye salmon

Regions of the genome displaying elevated differentiation (genomic islands of divergence) are thought to play an important role in local adaptation, especially in populations experiencing high gene flow. However, the characteristics of these islands as well as the functional significance of genes located within them remain largely unknown. Here, we used data from thousands of SNPs aligned to a linkage map to investigate genomic islands of divergence in three ecotypes of sockeye salmon (Oncorhynchus nerka) from a single drainage in southwestern Alaska. We found ten islands displaying high differentiation among ecotypes. Conversely, neutral structure observed throughout the rest of the genome was low and not partitioned by ecotype. One island on linkage group So13 was particularly large and contained six SNPs with F_ST  > 0.14 (average F_ST of neutral SNPs = 0.01). Functional annotation revealed that the peak of this island contained a nonsynonymous mutation in a gene involved in growth in other species (TULP4). The islands that we discovered were relatively small (80-402 Kb), loci found in islands did not show reduced levels of diversity, and loci in islands displayed slightly elevated linkage disequilibrium. These attributes suggest that the islands discovered here were likely generated by divergence hitchhiking; however, we cannot rule out the possibility that other mechanisms may have produced them. Our results suggest that islands of divergence serve an important role in local adaptation with gene flow and represent a significant advance towards understanding the genetic basis of ecotypic differentiation.

Salmonid Chromosome Evolution as Revealed by a Novel Method for Comparing RADseq Linkage Maps

Whole genome duplication (WGD) can provide material for evolutionary innovation. Family Salmonidae is ideal for studying the effects of WGD as the ancestral salmonid underwent WGD relatively recently, ∼65 Ma, then rediploidized and diversified. Extensive synteny between homologous chromosome arms occurs in extant salmonids, but each species has both conserved and unique chromosome arm fusions and fissions. Assembly of large, outbred eukaryotic genomes can be difficult, but structural rearrangements within such taxa can be investigated using linkage maps. RAD sequencing provides unprecedented ability to generate high-density linkage maps for nonmodel species, but can result in low numbers of homologous markers between species due to phylogenetic distance or differences in library preparation. Here, we generate a high-density linkage map (3,826 markers) for the Salvelinus genera (Brook Charr S. fontinalis), and then identify corresponding chromosome arms among the other available salmonid high-density linkage maps, including six species of Oncorhynchus, and one species for each of Salmo, Coregonus, and the nonduplicated sister group for the salmonids, Northern Pike Esox lucius for identifying post-duplicated homeologs. To facilitate this process, we developed MapComp to identify identical and proximate (i.e. nearby) markers between linkage maps using a reference genome of a related species as an intermediate, increasing the number of comparable markers between linkage maps by 5-fold. This enabled a characterization of the most likely history of retained chromosomal rearrangements post-WGD, and several conserved chromosomal inversions. Analyses of RADseq-based linkage maps from other taxa will also benefit from MapComp, available at: https://github.com/enormandeau/mapcomp/