Comparative genomic analysis of Atlantic salmon, Salmo salar, from Europe and North America

The generation of the first chromosome-level de novo genome assembly and the development and validation of a 50K SNP array for the St. John River aquaculture strain of North American Atlantic salmon

Atlantic salmon (Salmo salar) in Northeastern US and Eastern Canada has high economic value for the sport fishing and aquaculture industries. Large differences exist between the genomes of Atlantic salmon of European origin and North American (N.A.) origin. Given the genetic and genomic differences between the 2 lineages, it is crucial to develop unique genomic resources for N.A. Atlantic salmon. Here, we describe the resources that we recently developed for genomic and genetic research in N.A. Atlantic salmon aquaculture. Firstly, a new single nucleotide polymorphism (SNP) database for N.A. Atlantic salmon consisting of 3.1 million putative SNPs was generated using data from whole-genome resequencing of 80 N.A. Atlantic salmon individuals. Secondly, a high-density 50K SNP array enriched for the genic regions of the genome and containing 3 sex determination and 61 putative continent of origin markers was developed and validated. Thirdly, a genetic map composed of 27 linkage groups with 36K SNP markers was generated from 2,512 individuals in 141 full-sib families. Finally, a chromosome-level de novo genome assembly from a male N.A. Atlantic salmon from the St. John River aquaculture strain was generated using PacBio long reads. Information from Hi-C proximity ligation sequences and Bionano optical mapping was used to concatenate the contigs into scaffolds. The assembly contains 1,755 scaffolds and only 1,253 gaps, with a total length of 2.83 Gb and N50 of 17.2 Mb. A BUSCO analysis detected 96.2% of the conserved Actinopterygii genes in the assembly, and the genetic linkage information was used to guide the formation of 27 chromosome sequences. Comparative analysis with the reference genome assembly of the European Atlantic salmon confirmed that the karyotype differences between the 2 lineages are caused by a fission in chromosome Ssa01 and 3 chromosome fusions including the p arm of chromosome Ssa01 with Ssa23, Ssa08 with Ssa29, and Ssa26 with Ssa28. The genomic resources we have generated for Atlantic salmon provide a crucial boost for genetic research and for management of farmed and wild populations in this highly valued species.

Population genetic analysis reveals a geographically limited transition zone between two genetically distinct Atlantic salmon lineages in Norway

Atlantic salmon is characterized by a high degree of population genetic structure throughout its native range. However, while populations inhabiting rivers in Norway and Russia make up a significant proportion of salmon in the Atlantic, thus far, genetic studies in this region have only encompassed low to modest numbers of populations. Here, we provide the first "in-depth" investigation of population genetic structuring in the species in this region. Analysis of 18 microsatellites on >9,000 fish from 115 rivers revealed highly significant population genetic structure throughout, following a hierarchical pattern. The highest and clearest level of division separated populations north and south of the Lofoten region in northern Norway. In this region, only a few populations displayed intermediate genetic profiles, strongly indicating a geographically limited transition zone. This was further supported by a dedicated cline analysis. Population genetic structure was also characterized by a pattern of isolation by distance. A decline in overall genetic diversity was observed from the south to the north, and two of the microsatellites showed a clear decrease in number of alleles across the observed transition zone. Together, these analyses support results from previous studies, that salmon in Norway originate from two main genetic lineages, one from the Barents-White Sea refugium that recolonized northern Norwegian and adjacent Russian rivers, and one from the eastern Atlantic that recolonized the rest of Norway. Furthermore, our results indicate that local conditions in the limited geographic transition zone between the two observed lineages, characterized by open coastline with no obvious barriers to gene flow, are strong enough to maintain the genetic differentiation between them.

Genome-wide nucleotide diversity of hatchery-reared Atlantic and Mediterranean strains of brown trout Salmo trutta compared to wild Mediterranean populations

A genome-wide assessment of diversity is provided for wild Mediterranean brown trout Salmo trutta populations from headwater tributaries of the Orb River and from Atlantic and Mediterranean hatchery-reared strains that have been used for stocking. Double-digest restriction-site-associated DNA sequencing (dd-RADseq) was performed and the efficiency of de novo and reference-mapping approaches to obtain individual genotypes was compared. Large numbers of single nucleotide polymorphism (SNP) markers with similar genome-wide distributions were discovered using both approaches (196 639 v. 121 016 SNPs, respectively), with c. 80% of the loci detected de novo being also found with reference mapping, using the Atlantic salmon Salmo salar genome as a reference. Lower mapping density but larger nucleotide diversity (π) was generally observed near extremities of linkage groups, consistent with regions of residual tetrasomic inheritance observed in salmonids. Genome-wide diversity estimates revealed reduced polymorphism in hatchery strains (π = 0·0040 and π = 0·0029 in Atlantic and Mediterranean strains, respectively) compared to wild populations (π = 0·0049), a pattern that was congruent with allelic richness estimated from microsatellite markers. Finally, pronounced heterozygote deficiency was found in hatchery strains (Atlantic F_IS = 0·18; Mediterranean F_IS = 0·42), indicating that stocking practices may affect the genetic diversity in wild populations. These new genomic resources will provide important tools to define better conservation strategies in S. trutta.

Genomic Instability of the Sex-Determining Locus in Atlantic Salmon (Salmo salar)

Atlantic salmon and rainbow trout, like other members of the subfamily Salmoninae, are gonochoristic with male heterogamety. The finding that sex-linked genetic markers varied between species suggested that the sex-determining gene differs among salmonid species, or that there is one sex-determining gene that has the capacity to move around the genome. The discovery of sdY, the sex-determining gene in rainbow trout, and its presence in many male salmonids gave support to the latter. Additional evidence for a salmonid-specific, sex-determining jumping gene came from the mapping of the sex-determining locus to three different chromosomes in Tasmanian male Atlantic salmon lineages. To characterize the sex-determining region, we isolated three sdY containing BACs from an Atlantic salmon male library. Sequencing of these BACs yielded two contigs, one of which contained the sdY gene. Sequence analysis of the borders of male-specific and female/male common regions revealed highly repetitive sequences associated with mobile elements, which may allow an sdY cassette to jump around the genome. FISH analysis using a BAC or a plasmid containing the sdY gene showed that the sdY gene did indeed localize to the chromosomes where SEX had been mapped in different Tasmanian Atlantic salmon families. Moreover, the plasmid sdY gene probe hybridized primarily to one of the sex chromosomes as would be expected of a male-specific gene. Our results suggest that a common salmonid sex-determining gene (sdY) can move between three specific loci on chromosomes 2, 3, and 6, giving the impression that there are multiple SEX loci both within and between salmonid species.

Population genomic analyses of early-phase Atlantic Salmon (Salmo salar) domestication/captive breeding

Domestication can have adverse genetic consequences, which may reduce the fitness of individuals once released back into the wild. Many wild Atlantic salmon (Salmo salarL.) populations are threatened by anthropogenic influences, and they are supplemented with captively bred fish. The Atlantic salmon is also widely used in selective breeding programs to increase the mean trait values for desired phenotypic traits. We analyzed a genomewide set of SNPs in three domesticated Atlantic salmon strains and their wild conspecifics to identify loci underlying domestication. The genetic differentiation between domesticated strains and wild populations was low (F_ST < 0.03), and domesticated strains harbored similar levels of genetic diversity compared to their wild conspecifics. Only a few loci showed footprints of selection, and these loci were located in different linkage groups among the different wild population/hatchery strain comparisons. Simulated scenarios indicated that differentiation in quantitative trait loci exceeded that in neutral markers during the early phases of divergence only when the difference in the phenotypic optimum between populations was large. This study indicates that detecting selection using standard approaches in the early phases of domestication might be challenging unless selection is strong and the traits under selection show simple inheritance patterns.

Comparative mapping between Coho Salmon (Oncorhynchus kisutch) and three other salmonids suggests a role for chromosomal rearrangements in the retention of duplicated regions following a whole genome duplication event

Whole genome duplication has been implicated in evolutionary innovation and rapid diversification. In salmonid fishes, however, whole genome duplication significantly pre-dates major transitions across the family, and re-diploidization has been a gradual process between genomes that have remained essentially collinear. Nevertheless, pairs of duplicated chromosome arms have diverged at different rates from each other, suggesting that the retention of duplicated regions through occasional pairing between homeologous chromosomes may have played an evolutionary role across species pairs. Extensive chromosomal arm rearrangements have been a key mechanism involved in re-dipliodization of the salmonid genome; therefore, we investigated their influence on degree of differentiation between homeologs across salmon species. We derived a linkage map for coho salmon and performed comparative mapping across syntenic arms within the genus Oncorhynchus, and with the genus Salmo, to determine the phylogenetic relationship between chromosome arrangements and the retention of undifferentiated duplicated regions. A 6596.7 cM female coho salmon map, comprising 30 linkage groups with 7415 and 1266 nonduplicated and duplicated loci, respectively, revealed uneven distribution of duplicated loci along and between chromosome arms. These duplicated regions were conserved across syntenic arms across Oncorhynchus species and were identified in metacentric chromosomes likely formed ancestrally to the divergence of Oncorhynchus from Salmo. These findings support previous studies in which observed pairings involved at least one metacentric chromosome. Re-diploidization in salmon may have been prevented or retarded by the formation of metacentric chromosomes after the whole genome duplication event and may explain lineage-specific innovations in salmon species if functional genes are found in these regions.

Spatial and temporal genetic structure of a river-resident Atlantic salmon (Salmo salar) after millennia of isolation

The river-resident Salmo salar ("småblank") has been isolated from other Atlantic salmon populations for 9,500 years in upper River Namsen, Norway. This is the only European Atlantic salmon population accomplishing its entire life cycle in a river. Hydropower development during the last six decades has introduced movement barriers and changed more than 50% of the river habitat to lentic conditions. Based on microsatellites and SNPs, genetic variation within småblank was only about 50% of that in the anadromous Atlantic salmon within the same river. The genetic differentiation (F ST) between småblank and the anadromous population was 0.24. This is similar to the differentiation between anadromous Atlantic salmon in Europe and North America. Microsatellite analyses identified three genetic subpopulations within småblank, each with an effective population size Ne of a few hundred individuals. There was no evidence of reduced heterozygosity and allelic richness in contemporary samples (2005-2008) compared with historical samples (1955-56 and 1978-79). However, there was a reduction in genetic differentiation between sampling localities over time. SNP data supported the differentiation of småblank into subpopulations and revealed downstream asymmetric gene flow between subpopulations. In spite of this, genetic variation was not higher in the lower than in the upper areas. The meta-population structure of småblank probably maintains genetic variation better than one panmictic population would do, as long as gene flow among subpopulations is maintained. Småblank is a unique endemic island population of Atlantic salmon. It is in a precarious situation due to a variety of anthropogenic impacts on its restricted habitat area. Thus, maintaining population size and avoiding further habitat fragmentation are important.

Insight into genomic changes accompanying divergence: genetic linkage maps and synteny of Lucania goodei and L. parva reveal a Robertsonian fusion

Linkage maps are important tools in evolutionary genetics and in studies of speciation. We performed a karyotyping study and constructed high-density linkage maps for two closely related killifish species, Lucania parva and L. goodei, that differ in salinity tolerance and still hybridize in their contact zone in Florida. Using SNPs from orthologous EST contigs, we compared synteny between the two species to determine how genomic architecture has shifted with divergence. Karyotyping revealed that L. goodei possesses 24 acrocentric chromosomes (1N) whereas L. parva possesses 23 chromosomes (1N), one of which is a large metacentric chromosome. Likewise, high-density single-nucleotide polymorphism-based linkage maps indicated 24 linkage groups for L. goodei and 23 linkage groups for L. parva. Synteny mapping revealed two linkage groups in L. goodei that were highly syntenic with the largest linkage group in L. parva. Together, this evidence points to the largest linkage group in L. parva being the result of a chromosomal fusion. We further compared synteny between Lucania with the genome of a more distant teleost relative medaka (Oryzias latipes) and found good conservation of synteny at the chromosomal level. Each Lucania LG had a single best match with each medaka chromosome. These results provide the groundwork for future studies on the genetic architecture of reproductive isolation and salinity tolerance in Lucania and other Fundulidae.

A dense linkage map for Chinook salmon (Oncorhynchus tshawytscha) reveals variable chromosomal divergence after an ancestral whole genome duplication event

Comparisons between the genomes of salmon species reveal that they underwent extensive chromosomal rearrangements following whole genome duplication that occurred in their lineage 58−63 million years ago. Extant salmonids are diploid, but occasional pairing between homeologous chromosomes exists in males. The consequences of re-diploidization can be characterized by mapping the position of duplicated loci in such species. Linkage maps are also a valuable tool for genome-wide applications such as genome-wide association studies, quantitative trait loci mapping or genome scans. Here, we investigated chromosomal evolution in Chinook salmon (Oncorhynchus tshawytscha) after genome duplication by mapping 7146 restriction-site associated DNA loci in gynogenetic haploid, gynogenetic diploid, and diploid crosses. In the process, we developed a reference database of restriction-site associated DNA loci for Chinook salmon comprising 48528 non-duplicated loci and 6409 known duplicated loci, which will facilitate locus identification and data sharing. We created a very dense linkage map anchored to all 34 chromosomes for the species, and all arms were identified through centromere mapping. The map positions of 799 duplicated loci revealed that homeologous pairs have diverged at different rates following whole genome duplication, and that degree of differentiation along arms was variable. Many of the homeologous pairs with high numbers of duplicated markers appear conserved with other salmon species, suggesting that retention of conserved homeologous pairing in some arms preceded species divergence. As chromosome arms are highly conserved across species, the major resources developed for Chinook salmon in this study are also relevant for other related species.

Linkage maps of the Atlantic salmon (Salmo salar) genome derived from RAD sequencing

BACKGROUND

Genetic linkage maps are useful tools for mapping quantitative trait loci (QTL) influencing variation in traits of interest in a population. Genotyping-by-sequencing approaches such as Restriction-site Associated DNA sequencing (RAD-Seq) now enable the rapid discovery and genotyping of genome-wide SNP markers suitable for the development of dense SNP linkage maps, including in non-model organisms such as Atlantic salmon (Salmo salar). This paper describes the development and characterisation of a high density SNP linkage map based on SbfI RAD-Seq SNP markers from two Atlantic salmon reference families.

RESULTS

Approximately 6,000 SNPs were assigned to 29 linkage groups, utilising markers from known genomic locations as anchors. Linkage maps were then constructed for the four mapping parents separately. Overall map lengths were comparable between male and female parents, but the distribution of the SNPs showed sex-specific patterns with a greater degree of clustering of sire-segregating SNPs to single chromosome regions. The maps were integrated with the Atlantic salmon draft reference genome contigs, allowing the unique assignment of ~4,000 contigs to a linkage group. 112 genome contigs mapped to two or more linkage groups, highlighting regions of putative homeology within the salmon genome. A comparative genomics analysis with the stickleback reference genome identified putative genes closely linked to approximately half of the ordered SNPs and demonstrated blocks of orthology between the Atlantic salmon and stickleback genomes. A subset of 47 RAD-Seq SNPs were successfully validated using a high-throughput genotyping assay, with a correspondence of 97% between the two assays.

CONCLUSIONS

This Atlantic salmon RAD-Seq linkage map is a resource for salmonid genomics research as genotyping-by-sequencing becomes increasingly common. This is aided by the integration of the SbfI RAD-Seq SNPs with existing reference maps and the draft reference genome, as well as the identification of putative genes proximal to the SNPs. Differences in the distribution of recombination events between the sexes is evident, and regions of homeology have been identified which are reflective of the recent salmonid whole genome duplication.

Meiotic maps of sockeye salmon derived from massively parallel DNA sequencing

Background

Meiotic maps are a key tool for comparative genomics and association mapping studies. Next-generation sequencing and genotyping by sequencing are speeding the processes of SNP discovery and the development of new genetic tools, including meiotic maps for numerous species. Currently there are limited genetic resources for sockeye salmon, Oncorhynchus nerka. We develop the first dense meiotic map for sockeye salmon using a combination of novel SNPs found in restriction site associated DNA (RAD tags) and SNPs available from existing expressed sequence tag (EST) based assays.

Results

We discovered and genotyped putative SNPs in 3,430 RAD tags. We removed paralogous sequence variants leaving 1,672 SNPs; these were combined with 53 EST-based SNP genotypes for linkage mapping. The map contained 29 male and female linkage groups, consistent with the haploid chromosome number expected for sockeye salmon. The female map contains 1,057 loci spanning 4,896 cM, and the male map contains 1,118 loci spanning 4,220 cM. Regions of conservation with rainbow trout and synteny between the RAD based rainbow trout map and the sockeye salmon map were established.

Conclusions

Using RAD sequencing and EST-based SNP assays we successfully generated the first high density linkage map for sockeye salmon.

Chromosomal differences between European and North American Atlantic salmon discovered by linkage mapping and supported by fluorescence in situ hybridization analysis

Background

Geographical isolation has generated a distinct difference between Atlantic salmon of European and North American Atlantic origin. The European Atlantic salmon generally has 29 pairs of chromosomes and 74 chromosome arms whereas it has been reported that the North American Atlantic salmon has 27 chromosome pairs and an NF of 72. In order to predict the major chromosomal rearrangements causing these differences, we constructed a dense linkage map for Atlantic salmon of North American origin and compared it with the well-developed map for European Atlantic salmon.

Results

The presented male and female genetic maps for the North American subspecies of Atlantic salmon, contains 3,662 SNPs located on 27 linkage groups. The total lengths of the female and male linkage maps were 2,153 cM and 968 cM respectively, with males characteristically showing recombination only at the telomeres. We compared these maps with recently published SNP maps from European Atlantic salmon, and predicted three chromosomal reorganization events that we then tested using fluorescence in situ hybridization (FISH) analysis. The proposed rearrangements, which define the differences in the karyotypes of the North American Atlantic salmon relative to the European Atlantic salmon, include the translocation of the p arm of ssa01 to ssa23 and polymorphic fusions: ssa26 with ssa28, and ssa08 with ssa29.

Conclusions

This study identified major chromosomal differences between European and North American Atlantic salmon. However, while gross structural differences were significant, the order of genetic markers at the fine-resolution scale was remarkably conserved. This is a good indication that information from the International Cooperation to Sequence the Atlantic salmon Genome, which is sequencing a European Atlantic salmon, can be transferred to Atlantic salmon from North America.

Three decades of farmed escapees in the wild: a spatio-temporal analysis of Atlantic salmon population genetic structure throughout Norway

Each year, hundreds of thousands of domesticated farmed Atlantic salmon escape into the wild. In Norway, which is the world’s largest commercial producer, many native Atlantic salmon populations have experienced large numbers of escapees on the spawning grounds for the past 15–30 years. In order to study the potential genetic impact, we conducted a spatio-temporal analysis of 3049 fish from 21 populations throughout Norway, sampled in the period 1970–2010. Based upon the analysis of 22 microsatellites, individual admixture, F_ST and increased allelic richness revealed temporal genetic changes in six of the populations. These changes were highly significant in four of them. For example, 76% and 100% of the fish comprising the contemporary samples for the rivers Vosso and Opo were excluded from their respective historical samples at P = 0.001. Based upon several genetic parameters, including simulations, genetic drift was excluded as the primary cause of the observed genetic changes. In the remaining 15 populations, some of which had also been exposed to high numbers of escapees, clear genetic changes were not detected. Significant population genetic structuring was observed among the 21 populations in the historical (global F_ST = 0.038) and contemporary data sets (global F_ST = 0.030), although significantly reduced with time (P = 0.008). This reduction was especially distinct when looking at the six populations displaying temporal changes (global F_ST dropped from 0.058 to 0.039, P = 0.006). We draw two main conclusions: 1. The majority of the historical population genetic structure throughout Norway still appears to be retained, suggesting a low to modest overall success of farmed escapees in the wild; 2. Genetic introgression of farmed escapees in native salmon populations has been strongly population-dependent, and it appears to be linked with the density of the native population.