Aquaculture Europe 2023

September 18 - 21, 2023


Add To Calendar 21/09/2023 09:30:0021/09/2023 09:45:00Europe/ViennaAquaculture Europe 2023GENOME WIDE ASSOCIATION STUDIES OF PRODUCTION TRAITS IN RAINBOW TROUT USING WHOLE GENOME SEQUENCINGSchubert 5The European Aquaculture Societywebmaster@aquaeas.orgfalseDD/MM/YYYYaaVZHLXMfzTRLzDrHmAi181982


V. Thomas1 , F. Enez2 , P. Haffray2, B. Imarazene3, M. Dupont-Nivet1, P. Phocas1


1 Université Paris-Saclay, INRAE, AgroParisTech, GABI, 78 350 Jouy-en-Josas, France

2 SYSAAF, Station LPGP-INRAE, 35 042 Rennes, France

3 Les Sources de l’Avance, Pissos 40 410, France

E-mail :



 Rainbow trout is  one of the most commonly farmed salmonid species  across the world and  the most farmed fish species in France. Trout are reared in France for the production of large fillets which are consumed fresh or smoked . Thus, the main objectives of genetic selection are to improve growth and fillet’s yield. To assess growth and fillet’s quality various traits are measured, such as body weight, carcass  and fillet yields or fat content. Though some QTLs have been found in previous experiments

 , we still lack precise information about genes and biologicals mechanisms linked to those traits. Our aim was to estimate genetic parameters of production traits in a commercial line of rainbow trout and to detect and accurately localize  associated QTLs.


Materials and methods

 The stock was established from 3 generations of  a commercial selected line of  Les Sources de l’Avance  breeding company (Aqualande G roup, France) whose breeding program started 12 generations ago. The stock was reared under commercial conditions until harvesting. From the 9th, the 10th and the 12th generations, 2198 , 1410 and 824  fish  were sampled  between 503 and  539 days post fecundation. Fish were measured for fork length (FL, cm), body weight (BW, g), carcass yield (CY), head gutted carcass weight (HGCW, g), head gutted carcass yield (HGCY), viscera weight (VW, g), gut yield (GY), fat content in the muscle (FAT, %) recorded using a Fish Torry Fat-meter®.

Fin samples of  3390  phenotyped fish were  genotyped for 57,501 SNPs using the Axiom® Trout Genotyping array by the INRAE genotyping platform Gentyane (Clermont-Ferrand,  France). All the parents of phenotyped fish were also genotyped with the 57K array (184, 183 and 80 parents of the 9th generation,10th and 12th generations).

 After SNPs quality control including filtering out SNPs with  minor allele frequency (MAF) < 1%, and SNP call rate < 97%; 31,968 SNPs were retained for the analysis.

 The  99 sires  of the phenotyped fish in 10th  generation  were  also sequenced with the NovaSeq6000® paired-end technology (Illumina 2×150bp) at INRAE sequencing platform Get-PlaGe (Toulouse, France). We used the nfcore/sarek 2.7.1 pipeline to call variants . We retained 10,520,443  bi-allelic SNPs  with  MAF > 1%.   The imputation of the 32K genotypes of phenotyped fish into 10,520K genotypes was performed using  FIMPUTE3  software using pedigree information.  After imputation we kept 1,231, 034 SNPs with a MAF above 10%, a mendelian error rate below 3% and linkage disequilibrium r2 between SNPs < 0.9 in sliding 100kb-windows.

 Genetic parameters were estimated using BLUPF90 software

 . GWAS was performed for the 8  traits of interest with a B ayesian S parse L inear M ixed M odel (BSLMM) using GEMMA - 0.98.5 software


Results and discussion

Heritability varied from moderate to high values (TABLE 1). The highest heritability was estimated for FAT (0,63 ± 0,03) and the lowest value for FL (0,16 ± 0,02). P roduction traits were polygenic, the proportion of genetic variance accounted for by the 150 to 500 SNPs with the largest effects ranged from 75 to 95% depending on the trait.

 We detected 17 QTLs (TABLE 2) with very strong evidence based on posterior inclusion probability (PIP) of SNPs in BSLMM and associated Bayes Factor (2lnBF > 14). We did not re-detect seven of the ten QTLs found in the previous study based on generations 9th and 10th of the same rainbow trout line (Blay et al., 2021).

 For FL, 3 QTLs were identified, among them the QTL at 7.1 Mb on  chr 22 was also detected for BW and entirely located in the slc39a10  gene. F or BW, 4 other QTLs were detected on chr 4, 5 and 8. The QTL on chr 4 (with peak SNP at 20.2 Mb ) was located within rcan2.  The 2  QTLs on chr 8 and 22 were previously identified in the same trout line

 . The QTL on chr 8 spanned the region  of map3k7 and bach2b genes.  For HGCW and HGCY a common QTL was identified on chr 6 i n  a region where top3b and ppm1f genes are annotated .  All the other QTLs were in intergenic regions. 


Thanks to imputation on over 1 million SNPs, we significantly refined  the location of  2 QTLs previously identified, as well as we detected  15 new QTLs. These results will be validated in an undergoing study by genotyping the 13th  generation of the Aqualande’s line to confirm the associations between genotypes and production traits.


The European Maritime and Fisheries Fund supported this study (Omega-Truite PFEA470017FA1000008, Phenomix PFEA470019FA1000013, RedOut PFEA470019FA1000007).


Blay et al (2021). Frontiers in Genetics, 12.

Misztal et al (2014). Manual for BLUPF90 family of programs.

Zhou et al (2013). PLoS Genetics, 9(2), e1003264.