Genetic architecture of acute hyperthermia resistance in juvenile rainbow trout (Oncorhynchus mykiss) and genetic correlations with production traits.

BACKGROUND: Selective breeding is a promising solution to reduce the vulnerability of fish farms to heat waves, which are predicted to increase in intensity and frequency. However, limited information about the genetic architecture of acute hyperthermia resistance in fish is available. Two batches of sibs from a rainbow trout commercial line were produced: the first (N = 1382) was phenotyped for acute hyperthermia resistance at nine months of age and the second (N = 1506) was phenotyped for main production traits (growth, body length, muscle fat content and carcass yield) at 20 months of age. Fish were genotyped on a 57 K single nucleotide polymorphism (SNP) array and their genotypes were imputed to high-density based on the parent's genotypes from a 665 K SNP array. RESULTS: The heritability estimate of resistance to acute hyperthermia was 0.29 ± 0.05, confirming the potential of selective breeding for this trait. Since genetic correlations of acute hyperthermia resistance with the main production traits near harvest age were all close to zero, selecting for acute hyperthermia resistance should not impact the main production traits, and vice-versa. A genome-wide association study revealed that resistance to acute hyperthermia is a highly polygenic trait, with six quantitative trait loci (QTL) detected, but explaining less than 5% of the genetic variance. Two of these QTL, including the most significant one, may explain differences in acute hyperthermia resistance across INRAE isogenic lines of rainbow trout. Differences in mean acute hyperthermia resistance phenotypes between homozygotes at the most significant SNP was 69% of the phenotypic standard deviation, showing promising potential for marker-assisted selection. We identified 89 candidate genes within the QTL regions, among which the most convincing functional candidates are dnajc7, hsp70b, nkiras2, cdk12, phb, fkbp10, ddx5, cygb1, enpp7, pdhx and acly. CONCLUSIONS: This study provides valuable insight into the genetic architecture of acute hyperthermia resistance in juvenile rainbow trout. We show that the selection potential for this trait is substantial and selection for this trait should not be too detrimental to improvement of other traits of interest. Identified functional candidate genes provide new knowledge on the physiological mechanisms involved in acute hyperthermia resistance, such as protein chaperoning, oxidative stress response, homeostasis maintenance and cell survival.

Development of a High-Density 665 K SNP Array for Rainbow Trout Genome-Wide Genotyping.

Single nucleotide polymorphism (SNP) arrays, also named « SNP chips ¼, enable very large numbers of individuals to be genotyped at a targeted set of thousands of genome-wide identified markers. We used preexisting variant datasets from USDA, a French commercial line and 30X-coverage whole genome sequencing of INRAE isogenic lines to develop an Affymetrix 665 K SNP array (HD chip) for rainbow trout. In total, we identified 32,372,492 SNPs that were polymorphic in the USDA or INRAE databases. A subset of identified SNPs were selected for inclusion on the chip, prioritizing SNPs whose flanking sequence uniquely aligned to the Swanson reference genome, with homogenous repartition over the genome and the highest Minimum Allele Frequency in both USDA and French databases. Of the 664,531 SNPs which passed the Affymetrix quality filters and were manufactured on the HD chip, 65.3% and 60.9% passed filtering metrics and were polymorphic in two other distinct French commercial populations in which, respectively, 288 and 175 sampled fish were genotyped. Only 576,118 SNPs mapped uniquely on both Swanson and Arlee reference genomes, and 12,071 SNPs did not map at all on the Arlee reference genome. Among those 576,118 SNPs, 38,948 SNPs were kept from the commercially available medium-density 57 K SNP chip. We demonstrate the utility of the HD chip by describing the high rates of linkage disequilibrium at 2-10 kb in the rainbow trout genome in comparison to the linkage disequilibrium observed at 50-100 kb which are usual distances between markers of the medium-density chip.

Estimates of recent and historical effective population size in turbot, seabream, seabass and carp selective breeding programmes.

BACKGROUND: The high fecundity of fish species allows intense selection to be practised and therefore leads to fast genetic gains. Based on this, numerous selective breeding programmes have been started in Europe in the last decades, but in general, little is known about how the base populations of breeders have been built. Such knowledge is important because base populations can be created from very few individuals, which can lead to small effective population sizes and associated reductions in genetic variability. In this study, we used genomic information that was recently made available for turbot (Scophthalmus maximus), gilthead seabream (Sparus aurata), European seabass (Dicentrarchus labrax) and common carp (Cyprinus carpio) to obtain accurate estimates of the effective size for commercial populations. METHODS: Restriction-site associated DNA sequencing data were used to estimate current and historical effective population sizes. We used a novel method that considers the linkage disequilibrium spectrum for the whole range of genetic distances between all pairs of single nucleotide polymorphisms (SNPs), and thus accounts for potential fluctuations in population size over time. RESULTS: Our results show that the current effective population size for these populations is small (equal to or less than 50 fish), potentially putting the sustainability of the breeding programmes at risk. We have also detected important drops in effective population size about five to nine generations ago, most likely as a result of domestication and the start of selective breeding programmes for these species in Europe. CONCLUSIONS: Our findings highlight the need to broaden the genetic composition of the base populations from which selection programmes start, and suggest that measures designed to increase effective population size within all farmed populations analysed here should be implemented in order to manage genetic variability and ensure the sustainability of the breeding programmes.

Optimizing Genomic Prediction of Host Resistance to Koi Herpesvirus Disease in Carp.

Genomic selection (GS) is increasingly applied in breeding programs of major aquaculture species, enabling improved prediction accuracy and genetic gain compared to pedigree-based approaches. Koi Herpesvirus disease (KHVD) is notifiable by the World Organization for Animal Health and the European Union, causing major economic losses to carp production. GS has potential to breed carp with improved resistance to KHVD, thereby contributing to disease control. In the current study, Restriction-site Associated DNA sequencing (RAD-seq) was applied on a population of 1,425 common carp juveniles which had been challenged with Koi herpes virus, followed by sampling of survivors and mortalities. GS was tested on a wide range of scenarios by varying both SNP densities and the genetic relationships between training and validation sets. The accuracy of correctly identifying KHVD resistant animals using GS was between 8 and 18% higher than pedigree best linear unbiased predictor (pBLUP) depending on the tested scenario. Furthermore, minor decreases in prediction accuracy were observed with decreased SNP density. However, the genetic relationship between the training and validation sets was a key factor in the efficacy of genomic prediction of KHVD resistance in carp, with substantially lower prediction accuracy when the relationships between the training and validation sets did not contain close relatives.

Mapping and Sequencing of a Significant Quantitative Trait Locus Affecting Resistance to Koi Herpesvirus in Common Carp.

Cyprinids are the most highly produced group of fishes globally, with common carp being one of the most valuable species of the group. Koi herpesvirus (KHV) infections can result in high levels of mortality, causing major economic losses, and is listed as a notifiable disease by the World Organization for Animal Health. Selective breeding for host resistance has the potential to reduce morbidity and losses due to KHV. Therefore, improving knowledge about host resistance and methods of incorporating genomic data into breeding for resistance may contribute to a decrease in economic losses in carp farming. In the current study, a population of 1,425 carp juveniles, originating from a factorial cross between 40 sires and 20 dams was challenged with KHV. Mortalities and survivors were recorded and sampled for genotyping by sequencing using Restriction Site-Associated DNA sequencing (RADseq). Genome-wide association analyses were performed to investigate the genetic architecture of resistance to KHV. A genome-wide significant QTL affecting resistance to KHV was identified on linkage group 44, explaining approximately 7% of the additive genetic variance. Pooled whole genome resequencing of a subset of resistant (n = 60) and susceptible animals (n = 60) was performed to characterize QTL regions, including identification of putative candidate genes and functional annotation of associated polymorphisms. The TRIM25 gene was identified as a promising positional and functional candidate within the QTL region of LG 44, and a putative premature stop mutation in this gene was discovered.

Potential for Genetic Improvement of the Main Slaughter Yields in Common Carp With in vivo Morphological Predictors.

Common carp is a major aquaculture species worldwide, commonly sold alive but also as processed headless carcass or filets. However, recording of processing yields is impossible on live breeding candidates, and alternatives for genetic improvement are either sib selection based on slaughtered fish, or indirect selection on correlated traits recorded in vivo. Morphological predictors that can be measured on live fish and that correlate with real slaughter yields hence remain a possible alternative. To quantify the power of morphological predictors for genetic improvement of yields, we estimated genetic parameters of slaughter yields and various predictors in 3-year-old common carp reared communally under semi-intensive pond conditions. The experimental stock was established by a partial factorial design of 20 dams and 40 sires, and 1553 progenies were assigned to their parents using 12 microsatellites. Slaughter yields were highly heritable (h2 = 0.46 for headless carcass yield, 0.50 for filet yield) and strongly genetically correlated with each other (rg = 0.96). To create morphological predictors, external (phenotypes, 2D digitization) and internal measurements (ultrasound imagery) were recorded and combined by multiple linear regression to predict slaughter yields. The accuracy of the phenotypic prediction was high for headless carcass yield (R2 = 0.63) and intermediate for filet yield (R2 = 0.49). Interestingly, heritability of predicted slaughter yields (0.48-0.63) was higher than that of the real yields to predict, and had high genetic correlations with the real yields (rg = 0.84-0.88). In addition, both predicted yields were highly phenotypically and genetically correlated with each other (0.95 for both), suggesting that using predicted headless carcass yield in a breeding program would be a good way to also improve filet yield. Besides, two individual predictors (P1 and P2) included in the prediction models and two simple internal measurements (E4 and E23) exhibited intermediate to high heritability estimates (h2 = 0.34 - 0.72) and significant genetic correlations to the slaughter yields (rg = |0.39 - 0.83|). The results show that there is a solid potential for genetic improvement of slaughter yields by selecting for predictor traits recorded on live breeding candidates of common carp.

Accuracy of Genomic Evaluations of Juvenile Growth Rate in Common Carp (Cyprinus carpio) Using Genotyping by Sequencing.

Cyprinids are the most important group of farmed fish globally in terms of production volume, with common carp (Cyprinus carpio) being one of the most valuable species of the group. The use of modern selective breeding methods in carp is at a formative stage, implying a large scope for genetic improvement of key production traits. In the current study, a population of 1,425 carp juveniles, originating from a partial factorial cross between 40 sires and 20 dams, was used for investigating the potential of genomic selection (GS) for juvenile growth, an exemplar polygenic production trait. RAD sequencing was used to identify and genotype SNP markers for subsequent parentage assignment, construction of a medium density genetic map (12,311 SNPs), genome-wide association study (GWAS), and testing of GS. A moderate heritability was estimated for body length of carp at 120 days (as a proxy of juvenile growth) of 0.33 (s.e. 0.05). No genome-wide significant QTL was identified using a single marker GWAS approach. Genomic prediction of breeding values outperformed pedigree-based prediction, resulting in 18% improvement in prediction accuracy. The impact of reduced SNP densities on prediction accuracy was tested by varying minor allele frequency (MAF) thresholds, with no drop in prediction accuracy until the MAF threshold is set <0.3 (2,744 SNPs). These results point to the potential for GS to improve economically important traits in common carp breeding programs.

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 (h2 = 0.36-1.00), except survival in both periods and fat change (both patterns) during overwintering (h2 = 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 (rg = -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.