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

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.

Three recent sex chromosome-to-autosome fusions in a Drosophila virilis strain with high satellite DNA content

The karyotype, or number and arrangement of chromosomes, has varying levels of stability across both evolution and disease. Karyotype changes often originate from DNA breaks near the centromeres of chromosomes, which generally contain long arrays of tandem repeats or satellite DNA. Drosophila virilis possesses among the highest relative satellite abundances of studied species, with almost half its genome composed of three related 7 bp satellites. We discovered a strain of D. virilis that we infer recently underwent three independent chromosome fusion events involving the X and Y chromosomes, in addition to one subsequent fission event. Here, we isolate and characterize the four different karyotypes we discovered in this strain which we believe demonstrates remarkable genome instability. We discovered that one of the substrains with an X-autosome fusion has an X-to-Y chromosome nondisjunction rate 20 × higher than the D. virilis reference strain (21% vs 1%). Finally, we found an overall higher rate of DNA breakage in the substrain with higher satellite DNA compared to a genetically similar substrain with less satellite DNA. This suggests that satellite DNA abundance may play a role in the risk of genome instability. Overall, we introduce a novel system consisting of a single strain with four different karyotypes, which we believe will be useful for future studies of genome instability, centromere function, and sex chromosome evolution.

Genomic evidence of recent European introgression into North American farmed and wild Atlantic salmon

Gene flow between wild and domestic populations has been repeatedly demonstrated across a diverse range of taxa. Ultimately, the genetic impacts of gene flow from domestic into wild populations depend both on the degree of domestication and the original source of the domesticated population. Atlantic salmon, Salmo salar, used in North American aquaculture are ostensibly of North American origin. However, evidence of European introgression into North American aquaculture salmon has accumulated in recent decades, even though the use of diploid European salmon has never been approved in Canada. The full extent of such introgression as well as the potential impacts on wild salmon in the Northwest Atlantic remains uncertain. Here, we extend previous work comparing North American and European wild salmon (n = 5799) using a 220 K SNP array to quantify levels of recent European introgression into samples of domestic salmon, aquaculture escapees, and wild salmon collected throughout Atlantic Canada. Analysis of North American farmed salmon (n = 403) and escapees (n = 289) displayed significantly elevated levels of European ancestry by comparison with wild individuals (p < 0.001). Of North American farmed salmon sampled between 2011 and 2018, ~17% had more than 10% European ancestry and several individuals exceeded 40% European ancestry. Samples of escaped farmed salmon similarly displayed elevated levels of European ancestry, with two individuals classified as 100% European. Analysis of juvenile salmon collected in rivers proximate to aquaculture locations also revealed evidence of elevated European ancestry and larger admixture tract in comparison to individuals collected at distance from aquaculture. Overall, our results demonstrate that even though diploid European salmon have never been approved for use in Canada, individuals of full and partial European ancestry have been in use over the last decade, and that some of these individuals have escaped and hybridized in the wild.

Environmentally associated chromosomal structural variation influences fine-scale population structure of Atlantic Salmon (Salmo salar)

Chromosomal rearrangements (e.g., inversions, fusions, and translocations) have long been associated with environmental variation in wild populations. New genomic tools provide the opportunity to examine the role of these structural variants in shaping adaptive differences within and among wild populations of non-model organisms. In Atlantic Salmon (Salmo salar), variations in chromosomal rearrangements exist across the species natural range, yet the role and importance of these structural variants in maintaining adaptive differences among wild populations remains poorly understood. We genotyped Atlantic Salmon (n = 1429) from 26 populations within a highly genetically structured region of southern Newfoundland, Canada with a 220K SNP array. Multivariate analysis, across two independent years, consistently identified variation in a structural variant (translocation between chromosomes Ssa01 and Ssa23), previously associated with evidence of trans-Atlantic secondary contact, as the dominant factor influencing population structure in the region. Redundancy analysis suggested that variation in the Ssa01/Ssa23 chromosomal translocation is strongly correlated with temperature. Our analyses suggest environmentally mediated selection acting on standing genetic variation in genomic architecture introduced through secondary contact may underpin fine-scale local adaptation in Placentia Bay, Newfoundland, Canada, a large and deep embayment, highlighting the importance of chromosomal structural variation as a driver of contemporary adaptive divergence.

Divergent and linked selection shape patterns of genomic differentiation between European and North American Atlantic salmon (Salmo salar)

As populations diverge many processes can shape genomic patterns of differentiation. Regions of high differentiation can arise due to divergent selection acting on selected loci, genetic hitchhiking of nearby loci, or through repeated selection against deleterious alleles (linked background selection); this divergence may then be further elevated in regions of reduced recombination. Atlantic salmon (Salmo salar) from Europe and North America diverged >600,000 years ago and despite some evidence of secondary contact, the majority of genetic data indicate substantial divergence between lineages. This deep divergence with potential gene flow provides an opportunity to investigate the role of different mechanisms that shape the genomic landscape during early speciation. Here, using 184,295 single nucleotide polymorphisms (SNPs) and 80 populations, we investigate the genomic landscape of differentiation across the Atlantic Ocean with a focus on highly differentiated regions and the processes shaping them. We found evidence of high (mean F_ST  = 0.26) and heterogeneous genomic differentiation between continents. Genomic regions associated with high trans-Atlantic differentiation ranged in size from single loci (SNPs) within important genes to large regions (1-3 Mbp) on four chromosomes (Ssa06, Ssa13, Ssa16 and Ssa19). These regions showed signatures consistent with selection, including high linkage disequilibrium, despite no significant reduction in recombination. Genes and functional enrichment of processes associated with differentiated regions may highlight continental differences in ocean navigation and parasite resistance. Our results provide insight into potential mechanisms underlying differences between continents, and evidence of near-fixed and potentially adaptive trans-Atlantic differences concurrent with a background of high genome-wide differentiation supports subspecies designation in Atlantic salmon.

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.

Genetic Differentiation in Hatchery and Stocked Populations of Sea Trout in the Southern Baltic: Selection Evidence at SNP Loci

The impacts and interactions between hatchery-bred fish and wild fish populations has been a topic of active investigation in recent decades. In some instances, the benefits of stocking can be overshadowed by negative effects such as genetic introgression with natural populations, loss of genetic diversity, and dilution of local adaptations. Methods that facilitate the identification of stocked fish enable us to estimate not only the effectiveness of stocking but also the level of natural reproduction and the degree of hybridization. The longest Baltic river, the Vistula, also has the second highest discharge. Historically, it hosted numerous populations of the anadromous form of brown trout (sea trout); however, dam construction has since interfered with and reduced spawning migration to a rate that is much lower than before. Reduced spawning has resulted in a population collapse and a negative flow-on effect on commercial catches. In response, Poland (along with many other Baltic countries) initiated an intensive stocking program which continues today and which sees the average annual release of 700,000 smolts. As a consequence, today’s main-river and inshore catches come from stock-enhanced populations. High-throughput single-nucleotide polymorphism (SNP) genotyping was performed on samples of sea trout from southern Baltic populations; results suggest that a significant portion of the sea trout catches in the Vistula mouth region have direct hatchery origin and indicate the presence of Pomeranian specimens. SNP loci identified as outliers indicate a potential selection pressure that may be related with effects of hatchery breeding and mixing with natural populations. The brown trout SNP array applied in this study showed high effectiveness not only for population differentiation, but more importantly, it emerged as a sensitive tool to provide evidence of detection selection.

Identification of genomic regions regulating sex determination in Atlantic salmon using high density SNP data

Background

A complete understanding of the genetic basis for sexual determination and differentiation is necessary in order to implement efficient breeding schemes at early stages of development. Atlantic salmon belongs to the family Salmonidae of fishes and represents a species of great commercial value. Although the species is assumed to be male heterogametic with XY sex determination, the precise genetic basis of sexual development remains unclear. The complexity is likely associated to the relatively recent salmonid specific whole genome duplication that may be responsible for certain genome instability. This instability together with the capacity of the sex-determining gene to move across the genome as reported by previous studies, may explain that sexual development genes are not circumscribed to the same chromosomes in all members of the species. In this study, we have used a 220 K SNP panel developed for Atlantic salmon to identify the chromosomes explaining the highest proportion of the genetic variance for sex as well as candidate regions and genes associated to sexual development in this species.

Results

Results from regional heritability analysis showed that the chromosomes explaining the highest proportion of variance in these populations were Ssa02 (heritability = 0.42, SE = 0.12) and Ssa21 (heritability = 0.26, SE = 0.11). After pruning by linkage disequilibrium, genome-wide association analyses revealed 114 SNPs that were significantly associated with sex, being Ssa02 the chromosome containing a greatest number of regions. Close examination of the candidate regions evidenced important genes related to sex in other species of Class Actinopterygii, including SDY, genes from family SOX, RSPO1, ESR1, U2AF2A, LMO7, GNRH-R, DND and FIGLA.

Conclusions

The combined results from regional heritability analysis and genome-wide association have provided new advances in the knowledge of the genetic regulation of sex determination in Atlantic salmon, supporting that Ssa02 is the candidate chromosome for sex in this species and suggesting an alternative population lineage in Spanish wild populations according to the results from Ssa21.

Exploring a Pool-seq-only approach for gaining population genomic insights in nonmodel species

Developing genomic insights is challenging in nonmodel species for which resources are often scarce and prohibitively costly. Here, we explore the potential of a recently established approach using Pool-seq data to generate a de novo genome assembly for mining exons, upon which Pool-seq data are used to estimate population divergence and diversity. We do this for two pairs of sympatric populations of brown trout (Salmo trutta): one naturally sympatric set of populations and another pair of populations introduced to a common environment. We validate our approach by comparing the results to those from markers previously used to describe the populations (allozymes and individual-based single nucleotide polymorphisms [SNPs]) and from mapping the Pool-seq data to a reference genome of the closely related Atlantic salmon (Salmo salar). We find that genomic differentiation (F ST) between the two introduced populations exceeds that of the naturally sympatric populations (F ST = 0.13 and 0.03 between the introduced and the naturally sympatric populations, respectively), in concordance with estimates from the previously used SNPs. The same level of population divergence is found for the two genome assemblies, but estimates of average nucleotide diversity differ ( π ¯  ≈ 0.002 and π ¯  ≈ 0.001 when mapping to S. trutta and S. salar, respectively), although the relationships between population values are largely consistent. This discrepancy might be attributed to biases when mapping to a haploid condensed assembly made of highly fragmented read data compared to using a high-quality reference assembly from a divergent species. We conclude that the Pool-seq-only approach can be suitable for detecting and quantifying genome-wide population differentiation, and for comparing genomic diversity in populations of nonmodel species where reference genomes are lacking.

Genomic signatures and correlates of widespread population declines in salmon

Global losses of biodiversity are occurring at an unprecedented rate, but causes are often unidentified. Genomic data provide an opportunity to isolate drivers of change and even predict future vulnerabilities. Atlantic salmon (Salmo salar) populations have declined range-wide, but factors responsible are poorly understood. Here, we reconstruct changes in effective population size (Ne) in recent decades for 172 range-wide populations using a linkage-based method. Across the North Atlantic, Ne has significantly declined in >60% of populations and declines are consistently temperature-associated. We identify significant polygenic associations with decline, involving genomic regions related to metabolic, developmental, and physiological processes. These regions exhibit changes in presumably adaptive diversity in declining populations consistent with contemporary shifts in body size and phenology. Genomic signatures of widespread population decline and associated risk scores allow direct and potentially predictive links between population fitness and genotype, highlighting the power of genomic resources to assess population vulnerability.

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.

Comparing genomic signatures of domestication in two Atlantic salmon (Salmo salar L.) populations with different geographical origins

Selective breeding and genetic improvement have left detectable signatures on the genomes of domestic species. The elucidation of such signatures is fundamental for detecting genomic regions of biological relevance to domestication and improving management practices. In aquaculture, domestication was carried out independently in different locations worldwide, which provides opportunities to study the parallel effects of domestication on the genome of individuals that have been selected for similar traits. In this study, we aimed to detect potential genomic signatures of domestication in two independent pairs of wild/domesticated Atlantic salmon populations of Canadian and Scottish origins, respectively. Putative genomic regions under divergent selection were investigated using a 200K SNP array by combining three different statistical methods based either on allele frequencies (LFMM, Bayescan) or haplotype differentiation (Rsb). We identified 337 and 270 SNPs potentially under divergent selection in wild and hatchery populations of Canadian and Scottish origins, respectively. We observed little overlap between results obtained from different statistical methods, highlighting the need to test complementary approaches for detecting a broad range of genomic footprints of selection. The vast majority of the outliers detected were population-specific but we found four candidate genes that were shared between the populations. We propose that these candidate genes may play a role in the parallel process of domestication. Overall, our results suggest that genetic drift may have override the effect of artificial selection and/or point toward a different genetic basis underlying the expression of similar traits in different domesticated strains. Finally, it is likely that domestication may predominantly target polygenic traits (e.g., growth) such that its genomic impact might be more difficult to detect with methods assuming selective sweeps.

Chromosomal fusion and life history-associated genomic variation contribute to within-river local adaptation of Atlantic salmon

Chromosomal inversions have been implicated in facilitating adaptation in the face of high levels of gene flow, but whether chromosomal fusions also have similar potential remains poorly understood. Atlantic salmon are usually characterized by population structure at multiple spatial scales; however, this is not the case for tributaries of the Miramichi River in North America. To resolve genetic relationships between populations in this system and the potential for known chromosomal fusions to contribute to adaptation, we genotyped 728 juvenile salmon using a 50 K SNP array. Consistent with previous work, we report extremely weak overall population structuring (Global F_ST  = 0.004) and failed to support hierarchical structure between the river's two main branches. We provide the first genomic characterization of a previously described polymorphic fusion between chromosomes 8 and 29. Fusion genomic characteristics included high LD, reduced heterozygosity in the fused homokaryotes, and strong divergence between the fused and the unfused rearrangement. Population structure based on fusion karyotype was five times stronger than neutral variation (F_ST  = 0.019), and the frequency of the fusion was associated with summer precipitation supporting a hypothesis that this rearrangement may contribute local adaptation despite weak neutral differentiation. Additionally, both outlier variation among populations and a polygenic framework for characterizing adaptive variation in relation to climate identified a 250-Kb region of chromosome 9, including the gene six6 that has previously been linked to age-at-maturity and run-timing for this species. Overall, our results indicate that adaptive processes, independent of major river branching, are more important than neutral processes for structuring these populations.

Genetic evaluation of Atlantic salmon for growth traits incorporating SNP markers

Genetic evaluations of individual fish were calculated for growth traits in North American Atlantic salmon with and without inclusion of genetic markers. The number of SNP markers was reduced to 6,000 and further to 270 in order to reduce the problem of overparameterization. SNP genotypes were predicted for all ungenotyped animals in the pedigree. Analysis of traits used a model with polygenic effects and SNP markers together. Polygenic effects refer to the additive genetic effects that remain after accounting for SNP genotypes. SNP marker genotypes were included as covariates to evaluate fish for growth traits (weight and length) in different environments (freshwater and seawater) with genders separated. Including regressions on SNP marker genotypes reduced the sum of squares of residuals by 2.7%-12.5% and increased the variability of Mendelian sampling effects (i.e., within-family variation) compared to traditional animal model evaluations. Genetic evaluations may be carried out with a few hundred markers which may be more affordable for genotyping large numbers of fish.

Extensive hybridization following a large escape of domesticated Atlantic salmon in the Northwest Atlantic

Domestication is rife with episodes of interbreeding between cultured and wild populations, potentially challenging adaptive variation in the wild. In Atlantic salmon, Salmo salar, the number of domesticated individuals far exceeds wild individuals, and escape events occur regularly, yet evidence of the magnitude and geographic scale of interbreeding resulting from individual escape events is lacking. We screened juvenile Atlantic salmon using 95 single nucleotide polymorphisms following a single, large aquaculture escape in the Northwest Atlantic and report the landscape-scale detection of hybrid and feral salmon (27.1%, 17/18 rivers). Hybrids were reproductively viable, and observed at higher frequency in smaller wild populations. Repeated annual sampling of this cohort revealed decreases in the presence of hybrid and feral offspring over time. These results link previous observations of escaped salmon in rivers with reports of population genetic change, and demonstrate the potential negative consequences of escapes from net-pen aquaculture on wild populations.

Mapping of quantitative trait loci associated with size, shape, and parr mark traits using first- and second-generation backcrosses between European and North American Atlantic salmon (Salmo salar)

Little is known about the genetic architecture of traits important for salmonid restoration ecology. We mapped quantitative trait loci (QTL) using single nucleotide polymorphisms (SNPs) for juvenile body length, weight, shape, and vertical skin pigmentation patterns (parr marks) within three hybrid backcross families between European and North American subspecies of Atlantic salmon. Amounts of variation in skin colour and pattern quantified in the two second-generation transAtlantic families exceeded the ranges seen in purebred populations. GridQTL analyses using low-density female-specific linkage maps detected QTL showing experiment-wide significance on Ssa02, Ssa03, Ssa09, Ssa11, Ssa19, and Ssa26/28 for both length and weight; on Ssa04 and Ssa23 for parr mark number; on Ssa09 and Ssa13 for parr mark contrast; and on Ssa05, Ssa07, Ssa10, Ssa11, Ssa18, Ssa23, and Ssa26/28 for geometric morphometric shape coordinates. Pleiotrophic QTL on Ssa11 affected length, weight, and shape. No QTL was found that explained more than 10% of the phenotypic variance in pigmentation or shape traits. Each QTL was approximately positioned on the physical map of the Atlantic salmon genome. Some QTL locations confirmed previous studies but many were new. Studies like ours may increase the success of salmon restoration projects by enabling better phenotypic and genetic matching between introduced and extirpated strains.

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.

Range-wide parallel climate-associated genomic clines in Atlantic salmon

Clinal variation across replicated environmental gradients can reveal evidence of local adaptation, providing insight into the demographic and evolutionary processes that shape intraspecific diversity. Using 1773 genome-wide single nucleotide polymorphisms we evaluated latitudinal variation in allele frequency for 134 populations of North American and European Atlantic salmon (Salmo salar). We detected 84 (4.74%) and 195 (11%) loci showing clinal patterns in North America and Europe, respectively, with 12 clinal loci in common between continents. Clinal single nucleotide polymorphisms were evenly distributed across the salmon genome and logistic regression revealed significant associations with latitude and seasonal temperatures, particularly average spring temperature in both continents. Loci displaying parallel clines were associated with several metabolic and immune functions, suggesting a potential basis for climate-associated adaptive differentiation. These climate-based clines collectively suggest evidence of large-scale environmental associated differences on either side of the North Atlantic. Our results support patterns of parallel evolution on both sides of the North Atlantic, with evidence of both similar and divergent underlying genetic architecture. The identification of climate-associated genomic clines illuminates the role of selection and demographic processes on intraspecific diversity in this species and provides a context in which to evaluate the impacts of climate change.

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 genome scan for selection signatures comparing farmed Atlantic salmon with two wild populations: Testing colocalization among outlier markers, candidate genes, and quantitative trait loci for production traits

Comparative genome scans can be used to identify chromosome regions, but not traits, that are putatively under selection. Identification of targeted traits may be more likely in recently domesticated populations under strong artificial selection for increased production. We used a North American Atlantic salmon 6K SNP dataset to locate genome regions of an aquaculture strain (Saint John River) that were highly diverged from that of its putative wild founder population (Tobique River). First, admixed individuals with partial European ancestry were detected using STRUCTURE and removed from the dataset. Outlier loci were then identified as those showing extreme differentiation between the aquaculture population and the founder population. All Arlequin methods identified an overlapping subset of 17 outlier loci, three of which were also identified by BayeScan. Many outlier loci were near candidate genes and some were near published quantitative trait loci (QTLs) for growth, appetite, maturity, or disease resistance. Parallel comparisons using a wild, nonfounder population (Stewiacke River) yielded only one overlapping outlier locus as well as a known maturity QTL. We conclude that genome scans comparing a recently domesticated strain with its wild founder population can facilitate identification of candidate genes for traits known to have been under strong artificial selection.

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.

SNP Arrays for Species Identification in Salmonids

The use of SNP genotyping microarrays, developed in one species to analyze a closely related species for which genomic sequence information is scarce, enables the rapid development of a genomic resource (SNP information) without the need to develop new species-specific markers. Using large numbers of microarray SNPs offers the best chance to detect informative markers in nontarget species, markers that can very often be assayed using a lower throughput platform as is described in this paper.

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.

Effects of crossovers between homeologs on inheritance and population genomics in polyploid-derived salmonid fishes

A whole genome duplication occurred in the ancestor of all salmonid fishes some 50-100 million years ago. Early inheritance studies with allozymes indicated that loci in the salmonid genome are inherited disomically in females. However, some pairs of duplicated loci showed patterns of inheritance in males indicating pairing and recombination between homeologous chromosomes. Nearly 20% of loci in the salmonid genome are duplicated and share the same alleles (isoloci), apparently due to homeologous recombination. Half-tetrad analysis revealed that isoloci tend to be telomeric. These results suggested that residual tetrasomic inheritance of isoloci results from homeologous recombination near chromosome ends and that continued disomic inheritance resulted from homologous pairing of centromeric regions. Many current genetic maps of salmonids are based on single nucleotide polymorphisms and microsatellites that are no longer duplicated. Therefore, long sections of chromosomes on these maps are poorly represented, especially telomeric regions. In addition, preferential multivalent pairing of homeologs from the same species in F1 hybrids results in an excess of nonparental gametes (so-called pseudolinkage). We consider how not including duplicated loci has affected our understanding of population and evolutionary genetics of salmonids, and we discuss how incorporating these loci will benefit our understanding of population genomics.

Comparative mapping for bighead carp (Aristichthys nobilis) against model and non-model fishes provides insights into the genomic evolution of cyprinids

Comparative mapping provides an efficient method to connect genomes of non-model and model fishes. In this study, we used flanking sequences of the 659 microsatellites on a genetic map of bighead carp (Aristichthys nobilis) to comprehensively study syntenic relationships between bighead carp and nine model and non-model fishes. Of the five model and two food fishes with whole genome data, Cyprinus carpio showed the highest rate of positive BLAST hits (95.3 %) with bighead carp map, followed by Danio rerio (70.9 %), Oreochromis niloticus (21.7 %), Tetraodon nigroviridis (6.4 %), Gasterosteus aculeatus (5.2 %), Oryzias latipes (4.7 %) and Fugu rubripes (3.5 %). Chromosomal syntenic analyses showed that inversion was the basic chromosomal rearrangement during genomic evolution of cyprinids, and the extent of inversions and translocations was found to be positively correlated with evolutionary relationships among fishes studied. Among the five investigated cyprinids, linkage groups (LGs) of bighead carp, Hypophthalmichthys molitrix and Ctenopharyngodon idella exhibited a one-to-one relationship. Besides, LG 9 of bighead carp and homologous LGs of silver carp and grass carp all corresponded to the chromosomes 10 and 22 of zebrafish, suggesting that chromosomal fission may have occurred in the ancestor of zebrafish. On the other hand, LGs of bighead carp and common carp showed an approximate one-to-two relationship with extensive translocations, confirming the occurrence of a 4th whole genome duplication in common carp. This study provides insights into the understanding of genome evolution among cyprinids and would aid in transferring positional and functional information of genes from model fish like zebrafish to non-model fish like bighead carp.

Footprints of directional selection in wild Atlantic salmon populations: evidence for parasite-driven evolution?

Mechanisms of host-parasite co-adaptation have long been of interest in evolutionary biology; however, determining the genetic basis of parasite resistance has been challenging. Current advances in genome technologies provide new opportunities for obtaining a genome-scale view of the action of parasite-driven natural selection in wild populations and thus facilitate the search for specific genomic regions underlying inter-population differences in pathogen response. European populations of Atlantic salmon (Salmo salar L.) exhibit natural variance in susceptibility levels to the ectoparasite Gyrodactylus salaris Malmberg 1957, ranging from resistance to extreme susceptibility, and are therefore a good model for studying the evolution of virulence and resistance. However, distinguishing the molecular signatures of genetic drift and environment-associated selection in small populations such as land-locked Atlantic salmon populations presents a challenge, specifically in the search for pathogen-driven selection. We used a novel genome-scan analysis approach that enabled us to i) identify signals of selection in salmon populations affected by varying levels of genetic drift and ii) separate potentially selected loci into the categories of pathogen (G. salaris)-driven selection and selection acting upon other environmental characteristics. A total of 4631 single nucleotide polymorphisms (SNPs) were screened in Atlantic salmon from 12 different northern European populations. We identified three genomic regions potentially affected by parasite-driven selection, as well as three regions presumably affected by salinity-driven directional selection. Functional annotation of candidate SNPs is consistent with the role of the detected genomic regions in immune defence and, implicitly, in osmoregulation. These results provide new insights into the genetic basis of pathogen susceptibility in Atlantic salmon and will enable future searches for the specific genes involved.

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.

Consolidation of the genetic and cytogenetic maps of turbot (Scophthalmus maximus) using FISH with BAC clones

Bacterial artificial chromosomes (BAC) have been widely used for fluorescence in situ hybridization (FISH) mapping of chromosome landmarks in different organisms, including a few in teleosts. In this study, we used BAC-FISH to consolidate the previous genetic and cytogenetic maps of the turbot (Scophthalmus maximus), a commercially important pleuronectiform. The maps consisted of 24 linkage groups (LGs) but only 22 chromosomes. All turbot LGs were assigned to specific chromosomes using BAC probes obtained from a turbot 5× genomic BAC library. It consisted of 46,080 clones with inserts of at least 100 kb and <5 % empty vectors. These BAC probes contained gene-derived or anonymous markers, most of them linked to quantitative trait loci (QTL) related to productive traits. BAC clones were mapped by FISH to unique marker-specific chromosomal positions, which showed a notable concordance with previous genetic mapping data. The two metacentric pairs were cytogenetically assigned to LG2 and LG16, and the nucleolar organizer region (NOR)-bearing pair was assigned to LG15. Double-color FISH assays enabled the consolidation of the turbot genetic map into 22 linkage groups by merging LG8 with LG18 and LG21 with LG24. In this work, a first-generation probe panel of BAC clones anchored to the turbot linkage and cytogenetical map was developed. It is a useful tool for chromosome traceability in turbot, but also relevant in the context of pleuronectiform karyotypes, which often show small hardly identifiable chromosomes. This panel will also be valuable for further integrative genomics of turbot within Pleuronectiformes and teleosts, especially for fine QTL mapping for aquaculture traits, comparative genomics, and whole-genome assembly.

CCR7 is mainly expressed in teleost gills, where it defines an IgD+IgM- B lymphocyte subset

Chemokine receptor CCR7, the receptor for both CCL19 and CCL21 chemokines, regulates the recruitment and clustering of circulating leukocytes to secondary lymphoid tissues, such as lymph nodes and Peyer's patches. Even though teleost fish do not have either of these secondary lymphoid structures, we have recently reported a homolog to CCR7 in rainbow trout (Oncorhynchus mykiss). In the present work, we have studied the distribution of leukocytes bearing extracellular CCR7 in naive adult tissues by flow cytometry, observing that among the different leukocyte populations, the highest numbers of cells with membrane (mem)CCR7 were recorded in the gill (7.5 ± 2% CCR7(+) cells). In comparison, head kidney, spleen, thymus, intestine, and peripheral blood possessed <5% CCR7(+) cells. When CCR7 was studied at early developmental stages, we detected a progressive increase in gene expression and protein CCR7 levels in the gills throughout development. Surprisingly, the majority of the CCR7(+) cells in the gills were not myeloid cells and did not express membrane CD8, IgM, nor IgT, but expressed IgD on the cell surface. In fact, most IgD(+) cells in the gills expressed CCR7. Intriguingly, the IgD(+)CCR7(+) population did not coexpress memIgM. Finally, when trout were bath challenged with viral hemorrhagic septicemia virus, the number of CCR7(+) cells significantly decreased in the gills while significantly increased in head kidney. These results provide evidence of the presence of a novel memIgD(+)memIgM(-) B lymphocyte subset in trout that expresses memCCR7 and responds to viral infections. Similarities with IgD(+)IgM(-) subsets in mammals are discussed.

Evidence for multiple sex-determining loci in Tasmanian Atlantic salmon (Salmo salar)

Phenotypic sex in salmonids is determined primarily by a genetic male heterogametic system; yet, sex reversal can be accomplished via hormonal treatment. In Tasmanian Atlantic salmon aquaculture, to overcome problems associated with early sexual maturation in males, sex-reversed females are crossed with normal females to produce all female stock. However, phenotypic distinction of sex-reversed females (neo-males) from true males is problematic. We set out to identify genetic markers that could make this distinction. Microsatellite markers from chromosome 2 (Ssa02), to which the sex-determining locus (SEX) has been mapped in two Scottish Atlantic salmon families, did not predict sex in a pilot study of seven families. A TaqMan 64 SNP genome-wide scan suggested SEX was on Ssa06 in these families, and this was confirmed by microsatellite markers. A survey of 58 families in total representing 38 male lineages in the SALTAS breeding program found that 34 of the families had SEX on Ssa02, in 22 of the families SEX was on Ssa06, and two of the families had a third SEX locus, on Ssa03. A PCR test using primers designed from the recently published sdY gene is consistent with Tasmanian Atlantic salmon having a single sex-determining gene that may be located on at least three linkage groups.

Perspective on rhabdomyolysis-induced acute kidney injury and new treatment options

AIM

The purpose of this study was to assess the role of caspase-dependent apoptosis, caspase 1, calpain 1, inducible nitric oxide synthase (iNOS) and endothelial nitric oxide synthase (eNOS) and the protective effect of grape seed proanthocyanidin extract (GSPE) in the development of rhabdomyolysis-induced acute kidney injury (AKI).

MATERIALS AND METHODS

Twenty-one rats were divided into 3 groups - control, rhabdomyolysis and rhabdomyolysis + GSPE. Rhabdomyolysis was induced in the rhabdomyolysis and rhabdomyolysis + GSPE groups with the injection into both hind limbs of 10 ml/kg hypertonic (50%) glycerol following 24-hour dehydration on the 6th day. The rhabdomyolysis + GSPE group was given GSPE at 100 mg/kg by gavage for 7 days. The experiment was concluded 48 h after glycerol injection. Blood specimens were collected, and kidney tissues were extracted for histopathological examination.

RESULTS

We identified an increase in blood urea nitrogen, creatinine, histopathological score, iNOS, caspase 3, caspase 1 and calpain 1 expression in the rhabdomyolysis group compared to the controls and a decrease in eNOS expression. In the rhabdomyolysis + GSPE group, however, there was a decrease in these mediators, together with an increase in eNOS expression.

CONCLUSION

This study shows for the first time in the literature that calpain 1 is involved in the pathogenesis of rhabdomyolysis-induced AKI, and that GSPE may have a renoprotective effect.

Screen for Footprints of Selection during Domestication/Captive Breeding of Atlantic Salmon

Domesticated animals provide a unique opportunity to identify genomic targets of artificial selection to the captive environment. Here, we screened three independent domesticated/captive Atlantic salmon (Salmo salar) strains and their wild progenitor populations in an effort to detect potential signals of domestication selection by typing of 261 SNPs and 70 microsatellite loci. By combining information from four different neutrality tests, in total ten genomic regions showed signs of directional selection based on multiple sources of evidence. Most of the identified candidate regions were rather small ranging from zero to a few centimorgans (cM) in the female Atlantic salmon linkage map. We also evaluated how adaptation from standing variation affects adjacent SNP and microsatellite variation along the chromosomes and, by using forward simulations with strong selection, we were able to generate genetic differentiation patterns comparable to the observed data. This study highlights the significance of standing genetic variation during the early stages of adaptation and represents a useful step towards identifying functional variants involved in domestication of Atlantic salmon.