Conservation genomics of Atlantic salmon: SNPs associated with QTLs for adaptive traits in parr from four trans-Atlantic backcrosses

Risk assessment of antimony-arsenic contaminated soil remediated using zero-valent iron at different pH values combined with freeze-thaw cycles

Soil in mining wastelands is seriously polluted with heavy metals. Zero-valent iron (ZVI) is widely used for remediation of heavy metal-polluted soil because of its excellent adsorption properties; however, the remediation process is affected by complex environmental conditions, such as acid rain and freeze-thaw cycles. In this study, the effects of different pH values and freeze-thaw cycles on remediation of antimony (Sb)- and arsenic (As)-contaminated soil by ZVI were investigated in laboratory simulation experiments. The stability and potential human health risks associated with the remediated soil were evaluated. The results showed that ZVI has a significant stabilizing effect on Sb and As in both acidic and alkaline soils contaminated with dual levels of Sb and As, and the freeze-thaw process in different pH value solution systems further enhances the ability of ZVI to stabilize Sb and As, especially in acidic soils. However, it should be noted that apart from the pH=1.0 solution environment, ZVI's ability to stabilize As is attenuated under other circumstances, potentially leading to leaching of its unstable form and thereby increasing contamination risks. This indicates that the F1 (2% ZVI+pH=1 solution+freeze-thaw cycle) processing exhibits superior effectiveness. After F1 treatment, the bioavailability of Sb and As in both soils also significantly decreased during the gastric and intestinal stages (about 60.00%), the non-carcinogenic and carcinogenic risks of Sb and As in alkaline soils are eliminated for children and adults, with a decrease ranging from 60.00% to 70.00%, while in acidic soil, the non-carcinogenic and carcinogenic risks of As to adults and children is acceptable, but Sb still poses non-carcinogenic risks to children, despite reductions of about 65.00%. These findings demonstrate that soil pH is a crucial factor influencing the efficacy of ZVI in stabilizing Sb and As contaminants during freeze-thaw cycles. This provides a solid theoretical foundation for utilizing ZVI in the remediation of Sb- and As-contaminated soils, emphasizing the significance of considering both pH levels and freeze-thaw conditions to ensure effective and safe treatment.

Genomic underpinnings of head and body shape in Arctic charr ecomorph pairs

Across its Holarctic range, Arctic charr (Salvelinus alpinus) populations have diverged into distinct trophic specialists across independent replicate lakes. The major aspect of divergence between ecomorphs is in head shape and body shape, which are ecomorphological traits reflecting niche use. However, whether the genomic underpinnings of these parallel divergences are consistent across replicates was unknown but key for resolving the substrate of parallel evolution. We investigated the genomic basis of head shape and body shape morphology across four benthivore-planktivore ecomorph pairs of Arctic charr in Scotland. Through genome-wide association analyses, we found genomic regions associated with head shape (89 SNPs) or body shape (180 SNPs) separately and 50 of these SNPs were strongly associated with both body and head shape morphology. For each trait separately, only a small number of SNPs were shared across all ecomorph pairs (3 SNPs for head shape and 10 SNPs for body shape). Signs of selection on the associated genomic regions varied across pairs, consistent with evolutionary demography differing considerably across lakes. Using a comprehensive database of salmonid QTLs newly augmented and mapped to a charr genome, we found several of the head- and body-shape-associated SNPs were within or near morphology QTLs from other salmonid species, reflecting a shared genetic basis for these phenotypes across species. Overall, our results demonstrate how parallel ecotype divergences can have both population-specific and deeply shared genomic underpinnings across replicates, influenced by differences in their environments and demographic histories.

Is shape in the eye of the beholder? Assessing landmarking error in geometric morphometric analyses on live fish

Geometric morphometrics is widely used to quantify morphological variation between biological specimens, but the fundamental influence of operator bias on data reproducibility is rarely considered, particularly in studies using photographs of live animals taken under field conditions. We examined this using four independent operators that applied an identical landmarking scheme to replicate photographs of 291 live Atlantic salmon (Salmo salar L.) from two rivers. Using repeated measures tests, we found significant inter-operator differences in mean body shape, suggesting that the operators introduced a systematic error despite following the same landmarking scheme. No significant differences were detected when the landmarking process was repeated by the same operator on a random subset of photographs. Importantly, in spite of significant operator bias, small but statistically significant morphological differences between fish from the two rivers were found consistently by all operators. Pairwise tests of angles of vectors of shape change showed that these between-river differences in body shape were analogous across operator datasets, suggesting a general reproducibility of findings obtained by geometric morphometric studies. In contrast, merging landmark data when fish from each river are digitised by different operators had a significant impact on downstream analyses, highlighting an intrinsic risk of bias. Overall, we show that, even when significant inter-operator error is introduced during digitisation, following an identical landmarking scheme can identify morphological differences between populations. This study indicates that operators digitising at least a sub-set of all data groups of interest may be an effective way of mitigating inter-operator error and potentially enabling data sharing.

Redundancy analysis, genome-wide association studies and the pigmentation of brown trout (Salmo trutta L.)

The association of molecular variants with phenotypic variation is a main issue in biology, often tackled with genome-wide association studies (GWAS). GWAS are challenging, with increasing, but still limited, use in evolutionary biology. We used redundancy analysis (RDA) as a complimentary ordination approach to single- and multitrait GWAS to explore the molecular basis of pigmentation variation in brown trout (Salmo trutta) belonging to wild populations impacted by hatchery fish. Based on 75,684 single nucleotide polymorphic (SNP) markers, RDA, single- and multitrait GWAS allowed the extraction of 337 independent colour patterning loci (CPLs) associated with trout pigmentation traits, such as the number of red and black spots on flanks. Collectively, these CPLs (i) mapped onto 35 out of 40 brown trout linkage groups indicating a polygenic genomic architecture of pigmentation, (ii) were found to be associated with 218 candidate genes, including 197 genes formerly mentioned in the literature associated to skin pigmentation, skin patterning, differentiation or structure notably in a close relative, the rainbow trout (Onchorhynchus mykiss), and (iii) related to functions relevant to pigmentation variation (e.g., calcium- and ion-binding, cell adhesion). Annotated CPLs include genes with well-known pigmentation effects (e.g., PMEL, SLC45A2, SOX10), but also markers associated with genes formerly found expressed in rainbow or brown trout skins. RDA was also shown to be useful to investigate management issues, especially the dynamics of trout pigmentation submitted to several generations of hatchery introgression.

Quantitative Genetics of Smoltification Status at the Time of Seawater Transfer in Atlantic Salmon (Salmo Salar)

High mortality during grow out in the sea is a challenge for farmed Atlantic salmon production in Norway and globally, which is partly attributed to suboptimal smolt quality. In this study, two groups of pre-smolts were put on a standard light smoltification regime with alternating 12L:12D per day for 6 weeks (Phase I), followed by 24L:0D per day for 6 weeks (Phase II); one group was 0 + smolt (EXP1) and the other 1 + smolt (EXP2). To monitor the smoltification status of the fish, 100 (EXP1) and 60 (EXP2) fish were randomly sampled per week during Phase II. The following phenotypes for smoltification status were studied: RT-qPCR relative mRNA expression of values of two alpha catalytic subunits of the variants of the Na⁺K⁺ATPase (NKA) expressed in the sampled gill tissues of each fish. The first variant, alpha1a with increased expression in freshwater (FW) and the second variant alpha1b with increased expression in seawater variant (SW), as well as their ratio SW/FW. At the optimal time for seawater transfer based on the SW/FW trait, 1,000 (at sixth sampling of EXP1) and 1,500 (at fifth sampling of EXP2) fish were sampled for genetic parameter estimation. The individual variation in FW, SW, and SW/FW was very large at each of the seven samplings indicating a large variation among individuals in the optimum time of transfer to seawater. SW/FW showed significant genetic variation in both 0+ and 1+ smolts, which indicates the possibility for selection for improved synchronization of smoltification status of Atlantic salmon at the time where the largest proportion of the fish is considered to be smolt. However, the genetic correlation between SW/FW of 0+ and 1+ was not significantly different from zero indicating very little shared genetic variation in SW/FW in 0+ and 1+ fish. Smoltification phenotypes showed temporal progression over the smoltification period, and this progression varied between 0+ and 1+ smolt highlighting the importance of correctly timing the major sampling point, and when cohorts are transferred to seawater. This also highlighted the need for further research into noninvasive methods of objectively measuring individual smoltification through time and subsequent smolt survival and growth rate at sea.

Epistatic regulation of growth in Atlantic salmon revealed: a QTL study performed on the domesticated-wild interface

Background

Quantitative traits are typically considered to be under additive genetic control. Although there are indications that non-additive factors have the potential to contribute to trait variation, experimental demonstration remains scarce. Here, we investigated the genetic basis of growth in Atlantic salmon by exploiting the high level of genetic diversity and trait expression among domesticated, hybrid and wild populations.

Results

After rearing fish in common-garden experiments under aquaculture conditions, we performed a variance component analysis in four mapping populations totaling ~ 7000 individuals from six wild, two domesticated and three F1 wild/domesticated hybrid strains. Across the four independent datasets, genome-wide significant quantitative trait loci (QTLs) associated with weight and length were detected on a total of 18 chromosomes, reflecting the polygenic nature of growth. Significant QTLs correlated with both length and weight were detected on chromosomes 2, 6 and 9 in multiple datasets. Significantly, epistatic QTLs were detected in all datasets.

Discussion

The observed interactions demonstrated that the phenotypic effect of inheriting an allele deviated between half-sib families. Gene-by-gene interactions were also suggested, where the combined effect of two loci resulted in a genetic effect upon phenotypic variance, while no genetic effect was detected when the two loci were considered separately. To our knowledge, this is the first documentation of epistasis in a quantitative trait in Atlantic salmon. These novel results are of relevance for breeding programs, and for predicting the evolutionary consequences of domestication-introgression in wild populations.

Genome-wide association analyses reveal polygenic genomic architecture underlying divergent shell morphology in Spanish Littorina saxatilis ecotypes

Gene flow between diverging populations experiencing dissimilar ecological conditions can theoretically constrain adaptive evolution. To minimize the effect of gene flow, alleles underlying traits essential for local adaptation are predicted to be located in linked genome regions with reduced recombination. Local reduction in gene flow caused by selection is expected to produce elevated divergence in these regions. The highly divergent crab-adapted and wave-adapted ecotypes of the marine snail Littorina saxatilis present a model system to test these predictions. We used genome-wide association (GWA) analysis of geometric morphometric shell traits associated with microgeographic divergence between the two L. saxatilis ecotypes within three separate sampling sites. A total of 477 snails that had individual geometric morphometric data and individual genotypes at 4,066 single nucleotide polymorphisms (SNPs) were analyzed using GWA methods that corrected for population structure among the three sites. This approach allowed dissection of the genomic architecture of shell shape divergence between ecotypes across a wide geographic range, spanning two glacial lineages. GWA revealed 216 quantitative trait loci (QTL) with shell size or shape differences between ecotypes, with most loci explaining a small proportion of phenotypic variation. We found that QTL were evenly distributed across 17 linkage groups, and exhibited elevated interchromosomal linkage, suggesting a genome-wide response to divergent selection on shell shape between the two ecotypes. Shell shape trait-associated loci showed partial overlap with previously identified outlier loci under divergent selection between the two ecotypes, supporting the hypothesis of diversifying selection on these genomic regions. These results suggest that divergence in shell shape between the crab-adapted and wave-adapted ecotypes is produced predominantly by a polygenic genomic architecture with positive linkage disequilibrium among loci of small effect.

Comparative transcriptome analysis reveals potential evolutionary differences in adaptation of temperature and body shape among four Percidae species

Considering the divergent temperature habitats and morphological traits of four Percidae species: yellow perch (Perca flavescens), Eurasian perch (Perca fluviatilis), pike perch (Sander lucioperca), and ruffe (Gymnocephalus cernua), we stepped into the transcriptome level to discover genes and mechanisms that drive adaptation to different temperature environments and evolution in body shape. Based on 93,566 to 181,246 annotated unigenes of the four species, we identified 1,117 one-to-one orthologous genes and subsequently constructed the phylogenetic trees that are consistent with previous studies. Together with the tree, the ratios of nonsynonymous to synonymous substitutions presented decreased evolutionary rates from the D. rerio branch to the sub-branch clustered by P. flavescens and P. fluviatilis. The specific 93 fast-evolving genes and 57 positively selected genes in P. flavescens, compared with 22 shared fast-evolving genes among P. fluviatilis, G. cernua, and S. lucioperca, showed an intrinsic foundation that ensure its adaptation to the warmer Great Lakes and farther south, especially in functional terms like “Cul4-RING E3 ubiquitin ligase complex.” Meanwhile, the specific 78 fast-evolving genes and 41 positively selected genes in S. lucioperca drew a clear picture of how it evolved to a large and elongated body with camera-type eyes and muscle strength so that it could occupy the highest position in the food web. Overall, our results uncover genetic basis that support evolutionary adaptation of temperature and body shape in four Percid species, and could furthermore assist studies on environmental adaptation in fishes.

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.

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.

Genetic linkage of distinct adaptive traits in sympatrically speciating crater lake cichlid fish

Our understanding of how biological diversity arises is limited, especially in the case of speciation in the face of gene flow. Here we investigate the genomic basis of adaptive traits, focusing on a sympatrically diverging species pair of crater lake cichlid fishes. We identify the main quantitative trait loci (QTL) for two eco-morphological traits: body shape and pharyngeal jaw morphology. These traits diverge in parallel between benthic and limnetic species in the repeated adaptive radiations of this and other fish lineages. Remarkably, a single chromosomal region contains the highest effect size QTL for both traits. Transcriptomic data show that the QTL regions contain genes putatively under selection. Independent population genomic data corroborate QTL regions as areas of high differentiation between the sympatric sister species. Our results provide empirical support for current theoretical models that emphasize the importance of genetic linkage and pleiotropy in facilitating rapid divergence in sympatry.

The genetic architecture of growth and fillet traits in farmed Atlantic salmon (Salmo salar)

BACKGROUND

Performance and quality traits such as harvest weight, fillet weight and flesh color are of economic importance to the Atlantic salmon aquaculture industry. The genetic factors underlying these traits are of scientific and commercial interest. However, such traits are typically polygenic in nature, with the number and size of QTL likely to vary between studies and populations. The aim of this study was to investigate the genetic basis of several growth and fillet traits measured at harvest in a large farmed salmon population by using SNP markers. Due to the marked heterochiasmy in salmonids, an efficient two-stage mapping approach was applied whereby QTL were detected using a sire-based linkage analysis, a sparse SNP marker map and exploiting low rates of recombination, while a subsequent dam-based analysis focused on the significant chromosomes with a denser map to confirm QTL and estimate their position.

RESULTS

The harvest traits all showed significant heritability, ranging from 0.05 for fillet yield up to 0.53 for the weight traits. In the sire-based analysis, 1695 offspring with trait records and their 20 sires were successfully genotyped for the SNPs on the sparse map. Chromosomes 13, 18, 19 and 20 were shown to harbor genome-wide significant QTL affecting several growth-related traits. The QTL on chr. 13, 18 and 20 were detected in the dam-based analysis using 512 offspring from 10 dams and explained approximately 6-7 % of the within-family variation in these traits.

CONCLUSIONS

We have detected several QTL affecting economically important complex traits in a commercial salmon population. Overall, the results suggest that the traits are relatively polygenic and that QTL tend to be pleiotropic (affecting the weight of several components of the harvested fish). Comparison of QTL regions across studies suggests that harvest trait QTL tend to be relatively population-specific. Therefore, the application of marker or genomic selection for improvement in these traits is likely to be most effective when the discovery population is closely related to the selection candidates (e.g. within-family genomic selection).

Genome-wide association study (GWAS) for growth rate and age at sexual maturation in Atlantic salmon (Salmo salar)

Early sexual maturation is considered a serious drawback for Atlantic salmon aquaculture as it retards growth, increases production times and affects flesh quality. Although both growth and sexual maturation are thought to be complex processes controlled by several genetic and environmental factors, selection for these traits has been continuously accomplished since the beginning of Atlantic salmon selective breeding programs. In this genome-wide association study (GWAS) we used a 6.5K single-nucleotide polymorphism (SNP) array to genotype ∼480 individuals from the Cermaq Canada broodstock program and search for SNPs associated with growth and age at sexual maturation. Using a mixed model approach we identified markers showing a significant association with growth, grilsing (early sexual maturation) and late sexual maturation. The most significant associations were found for grilsing, with markers located in Ssa10, Ssa02, Ssa13, Ssa25 and Ssa12, and for late maturation with markers located in Ssa28, Ssa01 and Ssa21. A lower level of association was detected with growth on Ssa13. Candidate genes, which were linked to these genetic markers, were identified and some of them show a direct relationship with developmental processes, especially for those in association with sexual maturation. However, the relatively low power to detect genetic markers associated with growth (days to 5 kg) in this GWAS indicates the need to use a higher density SNP array in order to overcome the low levels of linkage disequilibrium observed in Atlantic salmon before the information can be incorporated into a selective breeding program.

Mapping QTLs for swimming ability related traits in Cyprinus carpio L

Body height (BH), head length (HL), snout length (SL), and tail length (TL) are important traits related with swimming ability of fish. Therefore, improving these traits will increase the production which is the basic goal of aquaculture breeding. To understand the genetic basis of swimming ability related traits in Cyprinus carpio L., a high-density linkage map spanning 3,301 cM in 50 linkage groups was utilized for quantitative trait locus (QTL) mapping. Mapping family comprised 190 offspring and 627 molecular markers were genotyped with average distance of 5.6 cM. A total of 15 QTLs including four (qBH13, qBH30, qBH33, qBH48) for BH, four (qHL10, qHL18, qHL29, qHL48) for HL, three (qSL24, qSL27, qSL45) for SL, and four (qTL15, qTL17, qTL18, qTL44) for TL were detected on 13 linkage groups LG10, LG13, LG15, LG17, LG18, LG24, LG27, LG29, LG30, LG33, LG44, LG45, and LG48. Each LG consisted on single QTL except LG18 and LG48. LG18 was found with two QTLs associated with HL and TL. While LG48 was comprised, the QTLs related with BH and HL. The phenotype variance was recorded from 12.6 to 40.6 %. Five QTLs, qHL48, qSL45, qTL15, qTL18, and qTL44, explained phenotype variance of >20 % with a significant levels of 0.047, 0.049, 0.037, 0.025, and 0.023, respectively. The neighbored loci of these QTLs were considered as main region of chromosomes controlling the traits. These identified genetic regions will be the main source of discovering gene(s) associated with swimming ability related traits in C. carpio L.

Single nucleotide polymorphisms in the insulin-like growth factor 1 (IGF1) gene are associated with growth-related traits in farmed Atlantic salmon

Understanding the genetic basis of variation in traits related to growth and fillet quality in Atlantic salmon is of importance to the aquaculture industry. Several growth-related QTL have been identified via the application of genetic markers. The IGF1 gene is considered a highly conserved and crucial growth-regulating gene in salmonid species. However, the association between polymorphisms in the IGF1 gene and growth-related traits in Atlantic salmon is unknown. Therefore, in this study, regions of the Atlantic salmon IGF1 gene were sequenced, aligned and compared across individuals. Three SNPs were identified in the putative promoter (SNP1, g.5763G>T; GenBank no. AGKD01012745), intron 1 (SNP2, g.7292C>T; GenBank no. AGKD01012745) and intron 3 (SNP3, g.4671A>C; GenBank no. AGKD01133398) regions respectively. These SNPs were genotyped in a population of 4800 commercial Atlantic salmon with data on several weight and fillet traits measured at harvest (at approximately 3 years of age). In a mixed model, association analysis of individual SNPs, SNP1 and SNP3 were both significantly associated with several weight traits (P < 0.05). The estimated additive effect on overall harvest weight was approximately 35 and 110 g for SNPs 1 and 3 respectively. A haplotype analysis confirmed the association between genetic variation in the IGF1 gene with overall body weight (P < 0.05) and fillet component traits (P < 0.05). Our findings suggest the identified nucleotide polymorphisms of the IGF1 gene may either affect farmed Atlantic salmon growth directly or be in population-wide linkage disequilibrium with causal variation, highlighting their possible utility as candidates for marker-assisted selection in the aquaculture industry.

A candidate gene association study for growth performance in an improved giant freshwater prawn (Macrobrachium rosenbergii ) culture line

A candidate gene approach using type I single nucleotide polymorphism (SNP) markers can provide an effective method for detecting genes and gene regions that underlie phenotypic variation in adaptively significant traits. In the absence of available genomic data resources, transcriptomes were recently generated in Macrobrachium rosenbergii to identify candidate genes and markers potentially associated with growth. The characterisation of 47 candidate loci by ABI re-sequencing of four cultured and eight wild samples revealed 342 putative SNPs. Among these, 28 SNPs were selected in 23 growth-related candidate genes to genotype in 200 animals selected for improved growth performance in an experimental GFP culture line in Vietnam. The associations between SNP markers and individual growth performance were then examined. For additive and dominant effects, a total of three exonic SNPs in glycogen phosphorylase (additive), heat shock protein 90 (additive and dominant) and peroxidasin (additive), and a total of six intronic SNPs in ankyrin repeats-like protein (additive and dominant), rolling pebbles (dominant), transforming growth factor-β induced precursor (dominant), and UTP-glucose-1-phosphate uridylyltransferase 2 (dominant) genes showed significant associations with the estimated breeding values in the experimental animals (P =0.001-0.031). Individually, they explained 2.6-4.8 % of the genetic variance (R²=0.026-0.048). This is the first large set of SNP markers reported for M. rosenbergii and will be useful for confirmation of associations in other samples or culture lines as well as having applications in marker-assisted selection in future breeding programs.

Genomic architecture of ecologically divergent body shape in a pair of sympatric crater lake cichlid fishes

Determining the genetic bases of adaptations and their roles in speciation is a prominent issue in evolutionary biology. Cichlid fish species flocks are a prime example of recent rapid radiations, often associated with adaptive phenotypic divergence from a common ancestor within a short period of time. In several radiations of freshwater fishes, divergence in ecomorphological traits - including body shape, colour, lips and jaws - is thought to underlie their ecological differentiation, specialization and, ultimately, speciation. The Midas cichlid species complex (Amphilophus spp.) of Nicaragua provides one of the few known examples of sympatric speciation where species have rapidly evolved different but parallel morphologies in young crater lakes. This study identified significant QTL for body shape using SNPs generated via ddRAD sequencing and geometric morphometric analyses of a cross between two ecologically and morphologically divergent, sympatric cichlid species endemic to crater Lake Apoyo: an elongated limnetic species (Amphilophus zaliosus) and a high-bodied benthic species (Amphilophus astorquii). A total of 453 genome-wide informative SNPs were identified in 240 F2 hybrids. These markers were used to construct a genetic map in which 25 linkage groups were resolved. Seventy-two segregating SNPs were linked to 11 QTL. By annotating the two most highly supported QTL-linked genomic regions, genes that might contribute to divergence in body shape along the benthic-limnetic axis in Midas cichlid sympatric adaptive radiations were identified. These results suggest that few genomic regions of large effect contribute to early stage divergence in Midas cichlids.

A SNP in the 5' flanking region of the myostatin-1b gene is associated with harvest traits in Atlantic salmon (Salmo salar)

Background

Myostatin (MSTN) belongs to the transforming growth factor-β superfamily and is a potent negative regulator of skeletal muscle development and growth in mammals. Most teleost fish possess two MSTN paralogues. However, as a consequence of a recent whole genome-duplication event, salmonids have four: MSTN-1 (−1a and -1b) and MSTN-2 (−2a and -2b). Evidence suggests that teleost MSTN plays a role in the regulation of muscle growth. In the current study, the MSTN-1b gene was re-sequenced and screened for SNP markers in a commercial population of Atlantic salmon. After genotyping 4,800 progeny for the discovered SNPs, we investigated their association with eight harvest traits - four body-weight traits, two ratios of weight traits, flesh colour and fat percentage - using a mixed model association analysis.

Results

Three novel SNPs were discovered in the MSTN-1b gene of Atlantic salmon. One of the SNPs, located within the 5′ flanking region (g.1086C > T), had a significant association with harvest traits (p < 0.05), specifically for: Harvest Weight (kg), Gutted Weight (kg), Deheaded Weight (kg) and Fillet Weight (kg). The haplotype-based association analysis was consistent with this result because the two haplotypes that showed a significant association with body-weight traits, hap4 and hap5 (p < 0.05 and p < 0.01, respectively), differ by a single substitution at the g.1086C > T locus. The alleles at g.1086C > T act in an additive manner and explain a small percentage of the genetic variation of these phenotypes.

Conclusions

The association analysis revealed that g.1086C > T had a significant association with all body-weight traits under study. Although the SNP explains a small percentage of the variance, our results indicate that a variation in the 5′ flanking region of the myostatin gene is associated with the genetic regulation of growth in Atlantic salmon.

Experimental tests for heritable morphological color plasticity in non-native brown trout (Salmo trutta) populations

The success of invasive species is frequently attributed to phenotypic plasticity, which facilitates persistence in novel environments. Here we report on experimental tests to determine whether the intensity of cryptic coloration patterns in a global invader (brown trout, Salmo trutta) was primarily the result of plasticity or heritable variation. Juvenile F₁ offspring were created through experimental crosses of wild-caught parents and reared for 30 days in the laboratory in a split-brood design on either light or dark-colored gravel substrate. Skin and fin coloration quantified with digital photography and image analysis indicated strong plastic effects in response to substrate color; individuals reared on dark substrate had both darker melanin-based skin color and carotenoid-based fin colors than other members of their population reared on light substrate. Slopes of skin and fin color reaction norms were parallel between environments, which is not consistent with heritable population-level plasticity to substrate color. Similarly, we observed weak differences in population-level color within an environment, again suggesting little genetic control on the intensity of skin and fin colors. Taken as whole, our results are consistent with the hypothesis that phenotypic plasticity may have facilitated the success of brown trout invasions and suggests that plasticity is the most likely explanation for the variation in color intensity observed among these populations in nature.

Application of single nucleotide polymorphisms to non-model species: a technical review

Single nucleotide polymorphisms (SNPs) have gained wide use in humans and model species and are becoming the marker of choice for applications in other species. Technology that was developed for work in model species may provide useful tools for SNP discovery and genotyping in non-model organisms. However, SNP discovery can be expensive, labour intensive, and introduce ascertainment bias. In addition, the most efficient approaches to SNP discovery will depend on the research questions that the markers are to resolve as well as the focal species. We discuss advantages and disadvantages of several past and recent technologies for SNP discovery and genotyping and summarize a variety of SNP discovery and genotyping studies in ecology and evolution.

Cost-effective genome-wide estimation of allele frequencies from pooled DNA in Atlantic salmon (Salmo salar L.)

Background

New sequencing technologies have tremendously increased the number of known molecular markers (single nucleotide polymorphisms; SNPs) in a variety of species. Concurrently, improvements to genotyping technology have now made it possible to efficiently genotype large numbers of genome-wide distributed SNPs enabling genome wide association studies (GWAS). However, genotyping significant numbers of individuals with large number of SNPs remains prohibitively expensive for many research groups. A possible solution to this problem is to determine allele frequencies from pooled DNA samples, such ‘allelotyping’ has been presented as a cost-effective alternative to individual genotyping and has become popular in human GWAS. In this article we have tested the effectiveness of DNA pooling to obtain accurate allele frequency estimates for Atlantic salmon (Salmo salar L.) populations using an Illumina SNP-chip.

Results

In total, 56 Atlantic salmon DNA pools from 14 populations were analyzed on an Atlantic salmon SNP-chip containing probes for 5568 SNP markers, 3928 of which were bi-allelic. We developed an efficient quality control filter which enables exclusion of loci showing high error rate and minor allele frequency (MAF) close to zero. After applying multiple quality control filters we obtained allele frequency estimates for 3631 bi-allelic loci. We observed high concordance (r > 0.99) between allele frequency estimates derived from individual genotyping and DNA pools. Our results also indicate that even relatively small DNA pools (35 individuals) can provide accurate allele frequency estimates for a given sample.

Conclusions

Despite of higher level of variation associated with array replicates compared to pool construction, we suggest that both sources of variation should be taken into account. This study demonstrates that DNA pooling allows fast and high-throughput determination of allele frequencies in Atlantic salmon enabling cost-efficient identification of informative markers for discrimination of populations at various geographical scales, as well as identification of loci controlling ecologically and economically important traits.

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.

The genetic architecture of growth rate in juvenile Takifugu species

Closely related species have often evolved dramatic differences in body size. Takifugu rubripes (fugu) is a large marine pufferfish whose genome has been sequenced, whereas T. niphobles is the smallest species among Takifugu. We show that, unsurprisingly, the juvenile growth rate of T. rubripes is higher than that of T. niphobles in a laboratory setting. We produced F(2) progenies of their F(1) hybrids and found one quantitative trait locus (QTL) significantly associated with variation in juvenile body size. This QTL region (3.5 Mb) contains no known genes directly related to growth phenotype (such as IGFs) except Fgf21, which inhibits growth hormone signaling in mouse. The QTL in Takifugu spp. is distinct from the region previously known to control body size variations in stickleback or tilapia. Our results suggest that in the fish tested herein, genomic regions underlying body size evolution might have different genetic origins. They also suggest that many diverse traits in Takifugu spp. are amenable to genetic mapping.

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.

Discovery and characterization of single nucleotide polymorphisms in Chinook salmon, Oncorhynchus tshawytscha

Molecular population genetics of non-model organisms has been dominated by the use of microsatellite loci over the last two decades. The availability of extensive genomic resources for many species is contributing to a transition to the use of single nucleotide polymorphisms (SNPs) for the study of many natural populations. Here we describe the discovery of a large number of SNPs in Chinook salmon, one of the world's most important fishery species, through large-scale Sanger sequencing of expressed sequence tag (EST) regions. More than 3 Mb of sequence was collected in a survey of variation in almost 132 kb of unique genic regions, from 225 separate ESTs, in a diverse ascertainment panel of 24 salmon. This survey yielded 117 TaqMan (5' nuclease) assays, almost all from separate ESTs, which were validated in population samples from five major stocks of salmon from the three largest basins on the Pacific coast of the contiguous United States: the Sacramento, Klamath and Columbia Rivers. The proportion of these loci that was variable in each of these stocks ranged from 86.3% to 90.6% and the mean minor allele frequency ranged from 0.194 to 0.236. There was substantial differentiation between populations with these markers, with a mean F(ST) estimate of 0.107, and values for individual loci ranging from 0 to 0.592. This substantial polymorphism and population-specific differentiation indicates that these markers will be broadly useful, including for both pedigree reconstruction and genetic stock identification applications.

Outlier SNPs show more genetic structure between two Bay of Fundy metapopulations of Atlantic salmon than do neutral SNPs

Atlantic salmon of Eastern Canada were once of considerable importance to aboriginal, recreational, and commercial fisheries, yet many populations are now in decline, particularly those of the inner Bay of Fundy (iBoF), which were recently listed as endangered. We investigated whether nonneutral SNPs could be used to assign individual Atlantic salmon accurately to either the iBoF or the outer Bay of Fundy (oBoF) metapopulations because this has been difficult with existing neutral markers. We first searched for markers under diversifying selection by genotyping eight captively bred Bay of Fundy (BoF) populations for 320 SNP loci with the Sequenom MassARRAY™ system and then analysed the data set with four different F(ST) outlier detection programs. Three outlier loci were identified by both BayesFST and BayeScan whereas seven outlier loci, including the three previously mentioned, were identified by both Fdist and Arlequin. A subset of 14 nonneutral SNPs was more accurate (85% accuracy) than a subset of 67 neutral SNPs (75% accuracy) at assigning individual salmon back to their metapopulation. We then chose a subset of nine outlier SNP markers and used them to inexpensively genotype archived DNA samples from seven wild BoF populations using Invader™ chemistry. Hierarchical AMOVA of these independent wild samples corroborated our previous findings of significant genetic differentiation between iBoF and oBoF salmon metapopulations. Our research shows that identifying and using outlier loci is an important step towards achieving the goal of consistently and accurately distinguishing iBoF from oBoF Atlantic salmon, which will aid in their conservation.

The genetic basis of divergent pigment patterns in juvenile threespine sticklebacks

Animal pigment patterns are important for a range of functions, including camouflage and communication. Repeating pigment patterns, such as stripes, bars and spots have been of particular interest to developmental and theoretical biologists, but the genetic basis of natural variation in such patterns is largely unexplored. In this study, we identify a difference in a periodic pigment pattern among juvenile threespine sticklebacks (Gasterosteus aculeatus) from different environments. Freshwater sticklebacks exhibit prominent vertical bars that visually break up the body shape, but sticklebacks from marine populations do not. We hypothesize that these distinct pigment patterns are tuned to provide crypsis in different habitats. This phenotypic difference is widespread and appears in most of the freshwater populations that we sampled. We used quantitative trait locus (QTL) mapping in freshwater-marine F2 hybrids to elucidate the genetic architecture underlying divergence in this pigmentation pattern. We identified two QTL that were significantly associated with variation in barring. Interestingly, these QTL were associated with two distinct aspects of the pigment pattern: melanophore number and overall pigment level. We compared the QTL locations with positions of known pigment candidate genes in the stickleback genome. We also identified two major QTL for juvenile body size, providing new insights into the genetic basis of juvenile growth rates in natural populations. In summary, although there is a growing literature describing simple genetic bases for adaptive coloration differences, this study emphasizes that pigment patterns can also possess a more complex genetic architecture.

Mapping QTL for an adaptive trait: the length of caudal fin in Lates calcarifer

The caudal fin represents a fundamental design feature of fishes and plays an important role in locomotor dynamics in fishes. The shape of caudal is an important parameter in traditional systematics. However, little is known about genes involved in the development of different forms of caudal fins. This study was conducted to identify and map quantitative trait loci (QTL) affecting the length of caudal fin and the ratio between tail length and standard body length in Asian seabass (Lates calcarifer). One F1 family containing 380 offspring was generated by crossing two unrelated individuals. One hundred and seventeen microsatellites almost evenly distributed along the whole genome were genotyped. Length of caudal fin at 90 days post-hatch was measured. QTL analysis detected six significant (genome-wide significant) and two suggestive (linkage-group-wide significant) QTL on seven linkage groups. The six significant QTL explained 5.5-16.6% of the phenotypic variance, suggesting these traits were controlled by multiple genes. Comparative genomics analysis identified several potential candidate genes for the length of caudal fin. The QTL for the length of caudal fin detected for the first time in marine fish may provide a starting point for the future identification of genes involved in the development of different forms of caudal fins in fishes.

Temporal change in genetic integrity suggests loss of local adaptation in a wild Atlantic salmon (Salmo salar) population following introgression by farmed escapees

In some wild Atlantic salmon populations, rapid declines in numbers of wild returning adults has been associated with an increase in the prevalence of farmed salmon. Studies of phenotypic variation have shown that interbreeding between farmed and wild salmon may lead to loss of local adaptation. Yet, few studies have attempted to assess the impact of interbreeding at the genome level, especially among North American populations. Here, we document temporal changes in the genetic makeup of the severely threatened Magaguadavic River salmon population (Bay of Fundy, Canada), a population that might have been impacted by interbreeding with farmed salmon for nearly 20 years. Wild and farmed individuals caught entering the river from 1980 to 2005 were genotyped at 112 single-nucleotide polymorphisms (SNPs), and/or eight microsatellite loci, to scan for potential shifts in adaptive genetic variation. No significant temporal change in microsatellite-based estimates of allele richness or gene diversity was detected in the wild population, despite its precipitous decline in numbers over the last two decades. This might reflect the effect of introgression from farmed salmon, which was corroborated by temporal change in linkage-disequilibrium. Moreover, SNP genome scans identified a temporal decrease in candidate loci potentially under directional selection. Of particular interest was a SNP previously shown to be strongly associated with an important quantitative trait locus for parr mark number, which retained its genetic distinctiveness between farmed and wild fish longer than other outliers. Overall, these results indicate that farmed escapees have introgressed with wild Magaguadavic salmon resulting in significant alteration of the genetic integrity of the native population, including possible loss of adaptation to wild conditions.

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

Background

Several lines of evidence including allozyme analysis, restriction digest patterns and sequencing of mtDNA as well as mini- and micro-satellite allele frequencies indicate that Atlantic salmon (Salmo salar) from North America and Europe are genetically distinct. These observations are supported by karyotype analysis, which revealed that North American Atlantic salmon have 27 pairs of chromosomes whereas European salmon have 29 pairs. We set out to construct a linkage map for a North American Atlantic salmon family and to compare this map with the well developed map for European Atlantic salmon.

Results

We used microsatellite markers, which had previously been mapped in the two Atlantic salmon SALMAP mapping families from the River Tay, Scotland, to carry out linkage analysis in an Atlantic salmon family (NB1) whose parents were derived from the Saint John River stock in New Brunswick, Canada. As large differences in recombination rates between female and male Atlantic salmon have been noted, separate genetic maps were constructed for each sex. The female linkage map comprises 218 markers in 37 linkage groups while the male map has 226 markers in 28 linkage groups. We combined 280 markers from the female and male maps into 27 composite linkage groups, which correspond to the haploid number of chromosomes in Atlantic salmon from the Western Atlantic.

Conclusions

A comparison of the composite NB1 and SALMAP linkage maps revealed the reason for the difference in the chromosome numbers between European and North American Atlantic salmon: Linkage groups AS-4 and AS-32 in the Scottish salmon, which correspond to chromosomes Ssa-6 and Ssa-22, are combined into a single NB1 linkage group as are linkage groups AS-21 and AS-33 (corresponding to chromosomes Ssa-26 and Ssa-28). The comparison of the linkage maps also suggested some additional chromosomal rearrangements, but it will require finer mapping, potentially using SNPs, to test these predictions. Our results provide the first comparison of the genomic architecture of Atlantic salmon from North America and Europe with respect to chromosome organization.

Growth-related quantitative trait loci in domestic and wild rainbow trout (Oncorhynchus mykiss)

Background

Somatic growth is a complex process that involves the action and interaction of genes and environment. A number of quantitative trait loci (QTL) previously identified for body weight and condition factor in rainbow trout (Oncorhynchus mykiss), and two other salmonid species, were used to further investigate the genetic architecture of growth-influencing genes in this species. Relationships among previously mapped candidate genes for growth and their co-localization to identified QTL regions are reported. Furthermore, using a comparative genomic analysis of syntenic rainbow trout linkage group clusters to their homologous regions within model teleost species such as zebrafish, stickleback and medaka, inferences were made regarding additional possible candidate genes underlying identified QTL regions.

Results

Body weight (BW) QTL were detected on the majority of rainbow trout linkage groups across 10 parents from 3 strains. However, only 10 linkage groups (i.e., RT-3, -6, -8, -9, -10, -12, -13, -22, -24, -27) possessed QTL regions with chromosome-wide or genome-wide effects across multiple parents. Fewer QTL for condition factor (K) were identified and only six instances of co-localization across families were detected (i.e. RT-9, -15, -16, -23, -27, -31 and RT-2/9 homeologs). Of note, both BW and K QTL co-localize on RT-9 and RT-27. The incidence of epistatic interaction across genomic regions within different female backgrounds was also examined, and although evidence for interaction effects within certain QTL regions were evident, these interactions were few in number and statistically weak. Of interest, however, was the fact that these predominantly occurred within K QTL regions. Currently mapped growth candidate genes are largely congruent with the identified QTL regions. More QTL were detected in male, compared to female parents, with the greatest number evident in an F₁ male parent derived from an intercross between domesticated and wild strain of rainbow trout which differed strongly in growth rate.

Conclusions

Strain background influences the degree to which QTL effects are evident for growth-related genes. The process of domestication (which primarily selects faster growing fish) may largely reduce the genetic influences on growth-specific phenotypic variation. Although heritabilities have been reported to be relatively high for both BW and K growth traits, the genetic architecture of K phenotypic variation appears less defined (i.e., fewer major contributing QTL regions were identified compared with BW QTL regions).

Mapping of quantitative trait loci for flesh colour and growth traits in Atlantic salmon (Salmo salar)

BACKGROUND

Flesh colour and growth related traits in salmonids are both commercially important and of great interest from a physiological and evolutionary perspective. The aim of this study was to identify quantitative trait loci (QTL) affecting flesh colour and growth related traits in an F2 population derived from an isolated, landlocked wild population in Norway (Byglands Bleke) and a commercial production population.

METHODS

One hundred and twenty-eight informative microsatellite loci distributed across all 29 linkage groups in Atlantic salmon were genotyped in individuals from four F2 families that were selected from the ends of the flesh colour distribution. Genotyping of 23 additional loci and two additional families was performed on a number of linkage groups harbouring putative QTL. QTL analysis was performed using a line-cross model assuming fixation of alternate QTL alleles and a half-sib model with no assumptions about the number and frequency of QTL alleles in the founder populations.

RESULTS

A moderate to strong phenotypic correlation was found between colour, length and weight traits. In total, 13 genome-wide significant QTL were detected for all traits using the line-cross model, including three genome-wide significant QTL for flesh colour (Chr 6, Chr 26 and Chr 4). In addition, 32 suggestive QTL were detected (chromosome-wide P < 0.05). Using the half-sib model, six genome-wide significant QTL were detected for all traits, including two for flesh colour (Chr 26 and Chr 4) and 41 suggestive QTL were detected (chromosome-wide P < 0.05). Based on the half-sib analysis, these two genome-wide significant QTL for flesh colour explained 24% of the phenotypic variance for this trait.

CONCLUSIONS

A large number of significant and suggestive QTL for flesh colour and growth traits were found in an F2 population of Atlantic salmon. Chr 26 and Chr 4 presented the strongest evidence for significant QTL affecting flesh colour, while Chr 10, Chr 5, and Chr 4 presented the strongest evidence for significant QTL affecting growth traits (length and weight). These QTL could be strong candidates for use in marker-assisted selection and provide a starting point for further characterisation of the genetic components underlying flesh colour and growth.

Perspectives on the genetic architecture of divergence in body shape in sticklebacks

The body shape of fishes encompasses a number of morphological traits that are intrinsically linked to functional systems and affect various measures of performance, including swimming, feeding, and avoiding predators. Changes in shape can allow a species to exploit a new ecological niche and can lead to ecological speciation. Body shape results from the integration of morphological, behavioral and physiological traits. It has been well established that functional interdependency among traits plays a large role in constraining the evolution of shape, affecting both the speed and the repeated evolution of particular body shapes. However, it is less clear what role genetic or developmental constraints might play in biasing the rate or direction of the evolution of body shape. Here, we suggest that the threespine stickleback (Gasterosteus aculeatus) is a powerful model system in which to address the extent to which genetic or developmental constraints play a role in the evolution of body shape in fishes. We review the existing data that begins to address these issues in sticklebacks and provide suggestions for future areas of research that will be particularly fruitful for illuminating the mechanisms that contribute to the evolution of body shape in fishes.

Detection of QTL affecting harvest traits in a commercial Atlantic salmon population

Genetic variation in performance and quality traits measured at harvest has previously been demonstrated in Atlantic salmon aquaculture populations. To map major loci underlying this variation, we utilized data from 10 families from a commercial breeding programme. Significant QTL were detected affecting harvest weight and length traits on linkage group 1, and affecting waste weight on linkage group 5. In total, 11 of the 29 linkage groups examined showed at least suggestive evidence for a QTL. These data suggest that major loci affecting economically important harvest characteristics are segregating in commercial salmon populations.