Genome-wide association identifies genomic regions influencing fillet color in Northwest Atlantic salmon (Salmo salar Linnaeus 1758)

Atlantic salmon (Salmo salar) is an important source of food globally; however, fillet color can significantly affect consumer purchasing, leading to potential food waste. Fish diets can be supplemented with astaxanthin to increase the organic pigment, carotenoid, responsible for flesh coloration; however, there is variation in the amount of overall fillet coloration in response to feeding astaxanthin. The uptake of this pigment is influenced by the environment and genetics and has been shown to be heritable. Therefore, we set out to determine the genomic associations of two separate year classes of farmed North American Atlantic salmon with measured Minolta Chroma Meter (lightness, redness, and yellowness) and SalmoFan phenotypic traits. Using ASReml-R genome-wide association, two genetic markers on chromosome 26 were significantly associated with almost all color traits, and these two markers explained between 6.0% and 12.5% of the variances. The genomic region on chromosome 26 was importantly found to be associated with the beta-carotene oxygenase 1 (bco1) gene, which is essential in the conversion of beta-carotenoids to vitamin A, implying that this gene may also play an important role in flesh coloration in North American Atlantic salmon. Additionally, there were several genomic regions significantly associated with color traits, in which the accompanying genes had functions in line with thermogenesis, immune function, and pathogenic responses. Understanding how environmental and genetic factors work together to affect fillet quality traits will help inform genetic improvement.

Genomic selection for survival under naturally occurring Saprolegnia oomycete infection in farmed European whitefish Coregonus lavaretus

Saprolegnia oomycete infection causes serious economic losses and reduces fish health in aquaculture. Genomic selection based on thousands of DNA markers is a powerful tool to improve fish traits in selective breeding programs. Our goal was to develop a single nucleotide polymorphism (SNP) marker panel and to test its use in genomic selection for improved survival against Saprolegnia infection in European whitefish Coregonus lavaretus, the second most important farmed fish species in Finland. We used a double digest restriction site associated DNA (ddRAD) genotyping by sequencing method to produce a SNP panel, and we tested it analyzing data from a cohort of 1,335 fish, which were measured at different times for mortality to Saprolegnia oomycete infection and weight traits. We calculated the genetic relationship matrix (GRM) from the genome-wide genetic data, integrating it in multivariate mixed models used for the estimation of variance components and genomic breeding values (GEBVs), and to carry out Genome-Wide Association Studies for the presence of quantitative trait loci (QTL) affecting the phenotypes in analysis. We identified one major QTL on chromosome 6 affecting mortality to Saprolegnia infection, explaining 7.7% to 51.3% of genetic variance, and a QTL for weight on chromosome 4, explaining 1.8% to 5.4% of genetic variance. Heritability for mortality was 0.20 to 0.43 on the liability scale, and heritability for weight was 0.44 to 0.53. The QTL for mortality showed an additive allelic effect. We tested whether integrating the QTL for mortality as a fixed factor, together with a new GRM calculated excluding the QTL from the genetic data, would improve the accuracy estimation of GEBVs. This test was done through a cross-validation approach, which indicated that the inclusion of the QTL increased the mean accuracy of the GEBVs by 0.28 points, from 0.33 to 0.61, relative to the use of full GRM only. The area under the curve of the receiver-operator curve for mortality increased from 0.58 to 0.67 when the QTL was included in the model. The inclusion of the QTL as a fixed effect in the model increased the correlation between the GEBVs of early mortality with the late mortality, compared to a model that did not include the QTL. These results validate the usability of the produced SNP panel for genomic selection in European whitefish and highlight the opportunity for modeling QTLs in genomic evaluation of mortality due to Saprolegnia infection.

Quantitative Genetics of Growth Rate and Filet Quality Traits in Atlantic Salmon Inferred From a Longitudinal Bayesian Model for the Left-Censored Gaussian Trait Growth Rate

In selective breeding programs for Atlantic salmon, test fish are slaughtered at an average body weight where growth rate and carcass traits as filet fat (F F), filet pigment (F P) and visceral fat index (F F) are recorded. The objective of this study was to obtain estimates of genetic correlations between growth rate (GR), and the three carcass quality traits when fish from the same 206 families (offspring of 120 sires and 206 dams from 2 year-classes) were recorded both at the same age (SA) and about the same body weight (SW). In the SW group, the largest fish were slaughtered at five different slaughter events and the remaining fish at the sixth slaughter event over 6 months. Estimates of genetic parameters for the traits were obtained from a Bayesian multivariate model for (potentially) truncated Gaussian traits through a Gibbs sampler procedure in which phantom GR values were obtained for the unslaughtered, and thus censored SW group fish at each slaughter event. The heritability estimates for the same trait in each group was similar; about 0.2 for FF, 0.15 for FP and 0.35 for VF and GR. The genetic correlation between the same traits in the two groups was high for growth rate (0.91 ± 0.05) visceral index (0.86 ± 0.05), medium for filet fat (0.45 ± 0.17) and low for filet pigment (0.13 ± 0.27). Within the two groups, the genetic correlation between growth rate and filet fat changed from positive (0.59 ± 0.14) for the SA group to negative (-0.45 ± 0.17) for the SW group, while the genetic correlation between growth rate and filet pigment changed from negative (-0.33 ± 0.22) for the SA group to positive (0.62 ± 0.16) for the SW group. The genetic correlation of growth rate with FF and FP is sensitive to whether the latter traits are measured at the same age or the same body weight. The results indicate that selection for increased growth rate is not expected to have a detrimental effect on the quality traits if increased growth potential is realized through a reduced production time.

Do whitefish (Coregonus lavaretus) larvae show adaptive variation in the avoidance of microplastic ingestion?

The presence of microplastics in aquatic ecosystems has recently received increased attention. Small plastic particles may resemble natural food items of larval fish and other aquatic organisms, and create strong selective pressures on the feeding traits in exposed populations. Here, we examined if larval ingestion of 90 μm polystyrene microspheres, in the presence of zooplankton (Artemia nauplii, mean length = 433 μm), shows adaptive variation in the European whitefish (Coregonus lavaretus). A full-factorial experimental breeding design allowed us to estimate the relative contributions of male (sire) and female (dam) parents and full-sib family variance in early feeding traits, and also genetic (co)variation between these traits. We also monitored the magnitude of intake and elimination of microplastics from the alimentary tracts of the larvae. In general, larval whitefish ingested small numbers of microplastics (mean = 1.8, range = 0-26 particles per larva), but ingestion was marginally affected by the dam, and more strongly by the full-sib family variation. Microsphere ingestion showed no statistically significant additive genetic variation, and thus, no heritability. Moreover, microsphere ingestion rate covaried positively with the ingestion of Artemia, further suggesting that larvae cannot adaptively avoid microsphere ingestion. Together with the detected strong genetic correlation between food intake and microplastic intake, the results suggest that larval fish do not readily possess additive genetic variation that would help them to adapt to the increasing pollution by microplastics. The conflict between feeding on natural food and avoiding microplastics deserves further attention.

Depuration of geosmin- and 2-methylisoborneol-induced off-flavors in recirculating aquaculture system (RAS) farmed European whitefish Coregonus lavaretus

European whitefish Coregonus lavaretus has increasingly become an important species for aquatic food production, especially in the Nordic countries. Whitefish is produced in traditional cage and pond operations, and in recirculating aquaculture system (RAS) in which, unfortunately, off-flavors and odors, mostly caused by geosmin (GSM) and 2-methylisoborneol (MIB), can accumulate in fish flesh from the circulating water. GSM and MIB have very low human sensory detection limits and, therefore, often disliked by consumers even at low concentrations. In this study, concentrations of GSM and MIB in RAS farmed European whitefish were studied by an analytical method based on headspace solid phase microextraction and gas chromatography–mass spectrometry. Concentrations were determined in different parts of fish: fillet, neck, belly, and tail during a depuration period and in depuration water. The highest initial concentrations were on average 32 ng g⁻¹ (GSM) and 24 ng g⁻¹ (MIB) in European whitefish fillet and 128 ng L⁻¹ (GSM) and 94 ng L⁻¹ (MIB) in water, respectively. After a depuration period of 16 days, concentrations decreased to below the detection limits, indicating the importance of the depuration period.

Genome-Wide Association Study Identifies Genomic Loci Affecting Filet Firmness and Protein Content in Rainbow Trout

Filet quality traits determine consumer satisfaction and affect profitability of the aquaculture industry. Soft flesh is a criterion for fish filet downgrades, resulting in loss of value. Filet firmness is influenced by many factors, including rate of protein turnover. A 50K transcribed gene SNP chip was used to genotype 789 rainbow trout, from two consecutive generations, produced in the USDA/NCCCWA selective breeding program. Weighted single-step GBLUP (WssGBLUP) was used to perform genome-wide association (GWA) analyses to identify quantitative trait loci affecting filet firmness and protein content. Applying genomic sliding windows of 50 adjacent SNPs, 212 and 225 SNPs were associated with genetic variation in filet shear force and protein content, respectively. Four common SNPs in the ryanodine receptor 3 gene (RYR3) affected the aforementioned filet traits; this association suggests common mechanisms underlying filet shear force and protein content. Genes harboring SNPs were mostly involved in calcium homeostasis, proteolytic activities, transcriptional regulation, chromatin remodeling, and apoptotic processes. RYR3 harbored the highest number of SNPs (n = 32) affecting genetic variation in shear force (2.29%) and protein content (4.97%). Additionally, based on single-marker analysis, a SNP in RYR3 ranked at the top of all SNPs associated with variation in shear force. Our data suggest a role for RYR3 in muscle firmness that may be considered for genomic- and marker-assisted selection in breeding programs of rainbow trout.

Co-inheritance of sea age at maturity and iteroparity in the Atlantic salmon vgll3 genomic region

Co-inheritance in life-history traits may result in unpredictable evolutionary trajectories if not accounted for in life-history models. Iteroparity (the reproductive strategy of reproducing more than once) in Atlantic salmon (Salmo salar) is a fitness trait with substantial variation within and among populations. In the Teno River in northern Europe, iteroparous individuals constitute an important component of many populations and have experienced a sharp increase in abundance in the last 20 years, partly overlapping with a general decrease in age structure. The physiological basis of iteroparity bears similarities to that of age at first maturity, another life-history trait with substantial fitness effects in salmon. Sea age at maturity in Atlantic salmon is controlled by a major locus around the vgll3 gene, and we used this opportunity demonstrate that these two traits are co-inherited around this genome region. The odds ratio of survival until second reproduction was up to 2.4 (1.8-3.5 90% CI) times higher for fish with the early-maturing vgll3 genotype (EE) compared to fish with the late-maturing genotype (LL). The L allele was dominant in individuals remaining only one year at sea before maturation, but the dominance was reversed, with the E allele being dominant in individuals maturing after two or more years at sea. Post hoc analysis indicated that iteroparous fish with the EE genotype had accelerated growth prior to first reproduction compared to first-time spawners, across all age groups, whereas this effect was not detected in fish with the LL genotype. These results broaden the functional link around the vgll3 genome region and help us understand constraints in the evolution of life-history variation in salmon. Our results further highlight the need to account for genetic correlations between fitness traits when predicting demographic changes in changing environments.

Sperm pre-fertilization thermal environment shapes offspring phenotype and performance

The sperm pre-fertilization environment has recently been suggested to mediate remarkable transgenerational consequences for offspring phenotype (transgenerational plasticity, TGB), but the adaptive significance of the process has remained unclear. Here, we studied the transgenerational effects of sperm pre-fertilization thermal environment in a cold-adapted salmonid, the European whitefish (Coregonus lavaretus). We used a full-factorial breeding design where the eggs of five females were fertilized with the milt of 10 males that had been pre-incubated at two different temperatures (3.5°C and 6.5°C) for 15 h prior to fertilization. Thermal manipulation did not affect sperm motility, cell size, fertilization success or embryo mortality. However, offspring that were fertilized with 6.5°C-exposed milt were smaller and had poorer swimming performance than their full-siblings that had been fertilized with the 3.5°C-exposed milt. Furthermore, the effect of milt treatment on embryo mortality varied among different females (treatment×female interaction) and male-female combinations (treatment×female×male interaction). Together, these results indicate that sperm pre-fertilization thermal environment shapes offspring phenotype and post-hatching performance and modifies both the magnitude of female (dam) effects and the compatibility of the gametes. Generally, our results suggest that short-term changes in sperm thermal conditions may have negative impact for offspring fitness. Thus, sperm thermal environment may have an important role in determining the adaptation potential of organisms to climate change. Detailed mechanism(s) behind our findings require further attention.

Gradual and Acute Temperature Rise Induces Crossing Endocrine, Metabolic, and Immunological Pathways in Maraena Whitefish (Coregonus maraena)

The complex and still poorly understood nature of thermoregulation in various fish species complicates the determination of the physiological status on the basis of diagnostic marker genes and indicative molecular pathways. The present study aimed to compare the physiological impacts of both gradual and acute temperature rise from 18 to 24°C on maraena whitefish in aquaculture. Microarray-based transcriptome profiles in the liver, spleen and kidney of heat-stressed maraena whitefish revealed the modulation of a significantly higher number of genes in those groups exposed to gradually rising temperatures compared with the acutely stressed groups, which might reflect early adaptation mechanisms. Moreover, we suggest a common set of 11 differentially expressed genes that indicate thermal stress induced by gradual or acute temperature rise in the three selected tissues. Besides the two pathways regulated in both data sets unfolded protein response and aldosterone signaling in epithelial cells, we identified unique tissue- and stress type-specific pathways reflecting the crossroads between signal transduction, metabolic and immunologic pathways to cope with thermal stress. In addition, comparing lists of differentially regulated genes with meta-analyzed published data sets revealed that “acute temperature rise”-responding genes that encode members of the HSP70, HSP90, and HSP40 families; their functional homologs; co-chaperones and stress-signal transducers are well-conserved across different species, tissues and/or cell types and experimental approaches.

The genetics of overwintering performance in two-year old common carp and its relation to performance until market size

Using farmed common carp, we investigated the genetic background of the second year overwintering performance and its relation to the performance during the third growing season and at market size. The experimental stock was established by partial factorial design with a series of 4 factorial matings of 5 dams and 10 sires each. The families were reared communally and pedigree was re-constructed with 93.6% success using 12 microsatellites on 2008 offspring. Three successive recordings (second autumn, third spring, and third autumn—market size) covering two periods (second overwintering, third growing season) were included. Body weight, Fulton’s condition factor and percent muscle fat content were recorded at all times and headless carcass yield and fillet yield were recorded at market size. Specific growth rate, absolute and relative fat change and overall survival were calculated for each period. Heritability estimates were significantly different from zero and almost all traits were moderately to highly heritable (h² = 0.36–1.00), except survival in both periods and fat change (both patterns) during overwintering (h² = 0.12–0.15). Genetic and phenotypic correlations imply that selection against weight loss and fat loss during overwintering is expected to lead to a better winter survival, together with a positive effect on growth in the third growing season. Interestingly, higher muscle fat content was genetically correlated to lower survival in the following period (r_g = -0.59; -0.53, respectively for winter and the third summer). On the other hand, higher muscle fat was also genetically linked to better slaughter yields. Moreover, selection for higher condition factor would lead to better performance during winter, growing season and at market size.

Genetic improvement of feed conversion ratio via indirect selection against lipid deposition in farmed rainbow trout (Oncorhynchus mykiss Walbaum)

In farmed fish, selective breeding for feed conversion ratio (FCR) may be possible via indirectly selecting for easily-measured indicator traits correlated with FCR. We tested the hypothesis that rainbow trout with low lipid% have genetically better FCR, and that lipid% may be genetically related to retention efficiency of macronutrients, making lipid% a useful indicator trait. A quantitative genetic analysis was used to quantify the benefit of replacing feed intake in a selection index with one of three lipid traits: body lipid%, muscle lipid% or viscera% weight of total body weight (reflecting visceral lipid). The index theory calculations showed that simultaneous selection for weight gain and against feed intake (direct selection to improve FCR) increased the expected genetic response in FCR by 1·50-fold compared with the sole selection for growth. Replacing feed intake in the selection index with body lipid%, muscle lipid% or viscera% increased genetic response in FCR by 1·29-, 1·49- and 1·02-fold, respectively, compared with the sole selection for growth. Consequently, indirect selection for weight gain and against muscle lipid% was almost as effective as direct selection for FCR. Fish with genetically low body and muscle lipid% were more efficient in turning ingested protein into protein weight gain. Both physiological and genetic mechanisms promote the hypothesis that low-lipid% fish are more efficient. These results highlight that in breeding programmes of rainbow trout, control of lipid deposition improves not only FCR but also protein-retention efficiency. This improves resource efficiency of aquaculture and reduces nutrient load to the environment.

Toll-like receptors in maraena whitefish: Evolutionary relationship among salmonid fishes and patterns of response to Aeromonas salmonicida

Toll-like receptors (TLRs) interact directly with particular pathogenic structures and are thus highly important to innate immunity. The present manuscript characterises a suite of 14 TLRs in maraena whitefish (Coregonus maraena), a salmonid species with increasing importance for aquaculture. Whitefish TLRs were structurally and evolutionary analysed. The results revealed a close relationship with TLRs from salmonid fish species rainbow trout and Atlantic salmon. Profiling the baseline expression of TLR genes in whitefish indicated that mainly members of the TLR11 family were highly expressed across all investigated tissues. A stimulation model with inactivated Aeromonas salmonicida was used to induce inflammation in the peritoneal cavity of whitefish. This bacterial challenge induced the expression of pro-inflammatory cytokine genes and evoked a strong influx of granulated cells of myeloid origin into the peritoneal cavity. As a likely consequence, the abundance of TLR-encoding transcripts increased moderately in peritoneal cells, with the highest levels of transcripts encoding non-mammalian TLR22a and a soluble TLR5 variant. In the course of inflammation, the proportion of granulated cells increased in peripheral blood accompanied by elevated TLR copy numbers in spleen and simultaneously reduced TLR copy numbers in head kidney at day 3 post-stimulation. Altogether, the present study provides in-vivo evidence for relatively modest TLR response patterns, but marked trafficking of myeloid cells as an immunophysiological consequence of A. salmonicida inflammation in whitefish. The present results contribute to improved understanding of the host-pathogen interaction in salmonid fish.

Identification and de novo sequencing of housekeeping genes appropriate for gene expression analyses in farmed maraena whitefish (Coregonus maraena) during crowding stress

Maraena whitefish (Coregonus maraena; synonym Coregonus lavaretus f. balticus) is a high-quality food fish in the Southern Baltic Sea belonging to the group of salmonid fishes. Coregonus sp. is successfully kept in aquaculture throughout northern Europe (e.g. in Finland, Germany, Russia) and North America. In this regard, the molecular and immunological characterisation of stress response in maraena whitefish contributes to the development of robust and fast-growing maraena whitefish breeding strains for aquaculture. Thus, in the present study, the potential housekeeping genes beta actin (ACTB), elongation factor 1 alpha (EEF1A1), glyceraldehydes-3-phosphate dehydrogenase (GAPDH), ribosomal protein 9 (RPL9), ribosomal protein 32 (RPL32) and ribosomal protein S20 (RPS20) were de novo sequenced and tested concerning their applicability as reference genes in quantitative real-time PCR (qPCR) in maraena whitefish under different stocking densities. For this purpose, tissue samples of liver, kidney, gills, head kidney, skin, adipose tissue, heart and dorsal fin were investigated. qPCR data were analysed with Normfinder tool to determine gene expression stability. DNA sequencing exposed transcribed paralogous EEF1A1A and EEF1A1B genes differing in their putative protein structure. Normfinder analysis revealed RPL9 and RPL32 as most stable, GAPDH and ACTB as least stable genes for qPCR analyses, respectively. This is the first study that provides a subset of seven de novo sequenced housekeeping genes usable as reference genes in studies of stress response in maraena whitefish.

Genotype-by-environment interaction of growth traits in rainbow trout (Oncorhynchus mykiss): a continental scale study

Rainbow trout is a globally important fish species for aquaculture. However, fish for most farms worldwide are produced by only a few breeding companies. Selection based solely on fish performance recorded at a nucleus may lead to lower-than-expected genetic gains in other production environments when genotype-by-environment (G × E) interaction exists. The aim was to quantify the magnitude of G × E interaction of growth traits (tagging weight; BWT, harvest weight; BWH, and growth rate; TGC) measured across 4 environments, located in 3 different continents, by estimating genetic correlations between environments. A total of 100 families, of at least 25 in size, were produced from the mating 58 sires and 100 dams. In total, 13,806 offspring were reared at the nucleus (selection environment) in Washington State (NUC) and in 3 other environments: a recirculating aquaculture system in Freshwater Institute (FI), West Virginia; a high-altitude farm in Peru (PE), and a cold-water farm in Germany (GER). To account for selection bias due to selective mortality, a multitrait multienvironment animal mixed model was applied to analyze the performance data in different environments as different traits. Genetic correlation (rg) of a trait measured in different environments and rg of different traits measured in different environments were estimated. The results show that heterogeneity of additive genetic variances was mainly found for BWH measured in FI and PE. Additive genetic coefficient of variation for BWH in NUC, FI, PE, and GER were 7.63, 8.36, 8.64, and 9.75, respectively. Genetic correlations between the same trait in different environments were low, indicating strong reranking (BWT: rg = 0.15 to 0.37, BWH: rg = 0.19 to 0.48, TGC: rg = 0.31 to 0.36) across environments. The rg between BWT in NUC and BWH in both FI (0.31) and GER (0.36) were positive, which was also found between BWT in NUC and TGC in both FI (0.10) and GER (0.20). However, rg were negative between BWT in NUC and both BWH (-0.06) and TGC (-0.20) in PE. Correction for selection bias resulted in higher additive genetic variances. In conclusion, strong G × E interaction was found for BWT, BWH, and TGC. Accounting for G × E interaction in the breeding program, either by using sib information from testing stations or environment-specific breeding programs, would increase genetic gains for environments that differ significantly from NUC.

Untangling the positive genetic correlation between rainbow trout growth and survival

Explanations for positive and negative genetic correlations between growth and fitness traits are essential for life-history theory and selective breeding. Here, we test whether growth and survival display genetic trade-off. Furthermore, we assess the potential of third-party traits to explain observed genetic associations. First, we estimated genetic correlations of growth and survival of rainbow trout. We then explored whether these associations are explained by genetic correlations with health, body composition and maturity traits. Analysis included 14 traits across life stages and environments. Data were recorded from 249 166 individuals belonging to 10 year classes of a pedigreed population. The results revealed that rapid growth during grow-out was genetically associated with enhanced survival (mean r_G = 0.17). This resulted because genotypes with less nematode caused cataract grew faster and were more likely to survive. Fingerling survival was not genetically related to weight or to grow-out survival. Instead, rapid fingerling growth made fish prone to deformations (r_G = 0.18). Evolutionary genetics provides a theoretical framework to study variation in genetic correlations. This study demonstrates that genetic correlation patterns of growth and survival can be explained by a set of key explanatory traits recorded at different life stages and that these traits can be simultaneously improved by selective breeding.