Genetic architecture and genomic selection of fatty acid composition predicted by Raman spectroscopy in rainbow trout.

BACKGROUND: In response to major challenges regarding the supply and sustainability of marine ingredients in aquafeeds, the aquaculture industry has made a large-scale shift toward plant-based substitutions for fish oil and fish meal. But, this also led to lower levels of healthful n-3 long-chain polyunsaturated fatty acids (PUFAs)-especially eicosapentaenoic (EPA) and docosahexaenoic (DHA) acids-in flesh. One potential solution is to select fish with better abilities to retain or synthesise PUFAs, to increase the efficiency of aquaculture and promote the production of healthier fish products. To this end, we aimed i) to estimate the genetic variability in fatty acid (FA) composition in visceral fat quantified by Raman spectroscopy, with respect to both individual FAs and groups under a feeding regime with limited n-3 PUFAs; ii) to study the genetic and phenotypic correlations between FAs and processing yields- and fat-related traits; iii) to detect QTLs associated with FA composition and identify candidate genes; and iv) to assess the efficiency of genomic selection compared to pedigree-based BLUP selection. RESULTS: Proportions of the various FAs in fish were indirectly estimated using Raman scattering spectroscopy. Fish were genotyped using the 57 K SNP Axiom™ Trout Genotyping Array. Following quality control, the final analysis contained 29,652 SNPs from 1382 fish. Heritability estimates for traits ranged from 0.03 ± 0.03 (n-3 PUFAs) to 0.24 ± 0.05 (n-6 PUFAs), confirming the potential for genomic selection. n-3 PUFAs are positively correlated to a decrease in fat deposition in the fillet and in the viscera but negatively correlated to body weight. This highlights the potential interest to combine selection on FA and against fat deposition to improve nutritional merit of aquaculture products. Several QTLs were identified for FA composition, containing multiple candidate genes with indirect links to FA metabolism. In particular, one region on Omy1 was associated with n-6 PUFAs, monounsaturated FAs, linoleic acid, and EPA, while a region on Omy7 had effects on n-6 PUFAs, EPA, and linoleic acid. When we compared the effectiveness of breeding programmes based on genomic selection (using a reference population of 1000 individuals related to selection candidates) or on pedigree-based selection, we found that the former yielded increases in selection accuracy of 12 to 120% depending on the FA trait. CONCLUSION: This study reveals the polygenic genetic architecture for FA composition in rainbow trout and confirms that genomic selection has potential to improve EPA and DHA proportions in aquaculture species.

Correction to: Quantitative trait loci for resistance to Flavobacterium psychrophilum in rainbow trout: effect of the mode of infection and evidence of epistatic interactions.

After publication of this work [1], we noted that there was an error in Table 3 Line 4.

Genome-wide estimates of genetic diversity, inbreeding and effective size of experimental and commercial rainbow trout lines undergoing selective breeding.

BACKGROUND: Selective breeding is a relatively recent practice in aquaculture species compared to terrestrial livestock. Nevertheless, the genetic variability of farmed salmonid lines, which have been selected for several generations, should be assessed. Indeed, a significant decrease in genetic variability due to high selection intensity could have occurred, potentially jeopardizing the long-term genetic progress as well as the adaptive capacities of populations facing change(s) in the environment. Thus, it is important to evaluate the impact of selection practices on genetic diversity to limit future inbreeding. The current study presents an analysis of genetic diversity within and between six French rainbow trout (Oncorhynchus mykiss) experimental or commercial lines based on a medium-density single nucleotide polymorphism (SNP) chip and various molecular genetic indicators: fixation index (FST), linkage disequilibrium (LD), effective population size (Ne) and inbreeding coefficient derived from runs of homozygosity (ROH). RESULTS: Our results showed a moderate level of genetic differentiation between selected lines (FST ranging from 0.08 to 0.15). LD declined rapidly over the first 100 kb, but then remained quite high at long distances, leading to low estimates of Ne in the last generation ranging from 24 to 68 depending on the line and methodology considered. These results were consistent with inbreeding estimates that varied from 10.0% in an unselected experimental line to 19.5% in a commercial line, and which are clearly higher than corresponding estimates in ruminants or pigs. In addition, strong variations in LD and inbreeding were observed along the genome that may be due to differences in local rates of recombination or due to key genes that tended to have fixed favorable alleles for domestication or production. CONCLUSIONS: This is the first report on ROH for any aquaculture species. Inbreeding appeared to be moderate to high in the six French rainbow trout lines, due to founder effects at the start of the breeding programs, but also likely to sweepstakes reproductive success in addition to selection for the selected lines. Efficient management of inbreeding is a major goal in breeding programs to ensure that populations can adapt to future breeding objectives and SNP information can be used to manage the rate at which inbreeding builds up in the fish genome.

Genetic and transcriptomic analyses provide new insights on the early antiviral response to VHSV in resistant and susceptible rainbow trout.

BACKGROUND: The viral hemorrhagic septicemia virus (VHSV) is a major threat for salmonid farming and for wild fish populations worldwide. Previous studies have highlighted the importance of innate factors regulated by a major quantitative trait locus (QTL) for the natural resistance to waterborne VHSV infection in rainbow trout. The aim of this study was to analyze the early transcriptomic response to VHSV inoculation in cell lines derived from previously described resistant and susceptible homozygous isogenic lines of rainbow trout to obtain insights into the molecular mechanisms responsible for the resistance to the viral infection. RESULTS: We first confirmed the presence of the major QTL in a backcross involving a highly resistant fish isogenic line (B57) and a highly susceptible one (A22), and were able to define the confidence interval of the QTL and to identify its precise position. We extended the definition of the QTL since it controls not only resistance to waterborne infection but also the kinetics of mortality after intra-peritoneal injection. Deep sequencing of the transcriptome of B57 and A22 derived cell lines exposed to inactivated VHSV showed a stronger response to virus inoculation in the resistant background. In line with our previous observations, an early and strong induction of interferon and interferon-stimulated genes was correlated with the resistance to VHSV, highlighting the major role of innate immune factors in natural trout resistance to the virus. Interestingly, major factors of the antiviral innate immunity were much more expressed in naive B57 cells compared to naive A22 cells, which likely contributes to the ability of B57 to mount a fast antiviral response after viral infection. These observations were further extended by the identification of several innate immune-related genes localized close to the QTL area on the rainbow trout genome. CONCLUSIONS: Taken together, our results improve our knowledge in virus-host interactions in vertebrates and provide novel insights in the molecular mechanisms explaining the resistance to VHSV in rainbow trout. Our data also provide a collection of potential markers for resistance and susceptibility of rainbow trout to VHSV infection.

Quantitative trait loci for resistance to Flavobacterium psychrophilum in rainbow trout: effect of the mode of infection and evidence of epistatic interactions.

BACKGROUND: Bacterial cold-water disease, which is caused by Flavobacterium psychrophilum, is one of the major diseases that affect rainbow trout (Oncorhynchus mykiss) and a primary concern for trout farming. Better knowledge of the genetic basis of resistance to F. psychrophilum would help to implement this trait in selection schemes and to investigate the immune mechanisms associated with resistance. Various studies have revealed that skin and mucus may contribute to response to infection. However, previous quantitative trait loci (QTL) studies were conducted by using injection as the route of infection. Immersion challenge, which is assumed to mimic natural infection by F. psychrophilum more closely, may reveal different defence mechanisms. RESULTS: Two isogenic lines of rainbow trout with contrasting susceptibilities to F. psychrophilum were crossed to produce doubled haploid F2 progeny. Fish were infected with F. psychrophilum either by intramuscular injection (115 individuals) or by immersion (195 individuals), and genotyped for 9654 markers using RAD-sequencing. Fifteen QTL associated with resistance traits were detected and only three QTL were common between the injection and immersion. Using a model that accounted for epistatic interactions between QTL, two main types of interactions were revealed. A "compensation-like" effect was detected between several pairs of QTL for the two modes of infection. An "enhancing-like" interaction effect was detected between four pairs of QTL. Integration of the QTL results with results of a previous transcriptomic analysis of response to F. psychrophilum infection resulted in a list of potential candidate immune genes that belong to four relevant functional categories (bacterial sensors, effectors of antibacterial immunity, inflammatory factors and interferon-stimulated genes). CONCLUSIONS: These results provide new insights into the genetic determinism of rainbow trout resistance to F. psychrophilum and confirm that some QTL with large effects are involved in this trait. For the first time, the role of epistatic interactions between resistance-associated QTL was evidenced. We found that the infection protocol used had an effect on the modulation of defence mechanisms and also identified relevant immune functional candidate genes.

Lack of correlation between the resistances to two rhabdovirus infections in rainbow trout.

The Viral Hemorrhagic Septicemia Virus (VHSV) and the Infectious Hematopoietic Necrosis Virus (IHNV) are two rhabdoviruses responsible for serious outbreaks in salmonid farms. To date, little is known about the variability of host response to these viruses. Using gynogenetic clonal lines of rainbow trout exhibiting a wide range of resistance to viral infections, we showed that there was no correlation between the resistance to VHSV and IHNV. We also confirmed the importance of fish weight for its susceptibility to IHNV infection. Finally, using a chimeric recombinant IHNV expressing the VHSV glycoprotein, we showed that the glycoprotein plays a key role in the virulence and in the level of resistance observed in different genetic backgrounds. Taken together, our results provide new prospects for a better understanding of host responses to rhabdovirus infections in salmonids.

Interplay between a bacterial pathogen and its host in rainbow trout isogenic lines with contrasted susceptibility to cold water disease.

Infectious diseases are a major constraint on aquaculture. Genetic lines with different susceptibilities to diseases are useful models to identify resistance mechanisms to pathogens and to improve prophylaxis. Bacterial cold-water disease (BCWD) caused by Flavobacterium psychrophilum represents a major threat for freshwater salmonid farming worldwide. A collection of rainbow trout (Oncorhynchus mykiss) isogenic lines was previously produced from a French domestic population. Here, we compared BCWD resistance phenotypes using a subset of isogenic lines chosen for their contrasted susceptibilities to F. psychrophilum. We applied individual monitoring to document the infection process, including time-course quantification of bacteremia and innate immune response. Strikingly, BCWD resistance was correlated with a lower bacterial growth rate in blood. Several immune genes were expressed at higher levels in resistant fish regardless of infection: the Type II arginase (arg2), a marker for M2 macrophages involved in anti-inflammatory responses and tissue repair, and two Toll-like receptors (tlr2/tlr7), responsible for pathogen detection and inflammatory responses. This study highlights the importance of innate and intrinsic defense mechanisms in determining the outcome of F. psychrophilum infections, and illustrates that non-lethal time-course blood sampling for individual monitoring of bacteremia is a powerful tool to resolve within-host pathogen behavior in bacterial fish diseases.

Detection of QTL with effects on osmoregulation capacities in the rainbow trout (Oncorhynchus mykiss).

BACKGROUND: There is increasing evidence that the ability to adapt to seawater in teleost fish is modulated by genetic factors. Most studies have involved the comparison of species or strains and little is known about the genetic architecture of the trait. To address this question, we searched for QTL affecting osmoregulation capacities after transfer to saline water in a nonmigratory captive-bred population of rainbow trout. RESULTS: A QTL design (5 full-sib families, about 200 F2 progeny each) was produced from a cross between F0 grand-parents previously selected during two generations for a high or a low cortisol response after a standardized confinement stress. When fish were about 18 months old (near 204 g body weight), individual progeny were submitted to two successive hyper-osmotic challenges (30 ppt salinity) 14 days apart. Plasma chloride and sodium concentrations were recorded 24 h after each transfer. After the second challenge, fish were sacrificed and a gill index (weight of total gill arches corrected for body weight) was recorded. The genome scan was performed with 196 microsatellites and 85 SNP markers. Unitrait and multiple-trait QTL analyses were carried out on the whole dataset (5 families) through interval mapping methods with the QTLMap software. For post-challenge plasma ion concentrations, significant QTL (P < 0.05) were found on six different linkage groups and highly suggestive ones (P < 0.10) on two additional linkage groups. Most QTL affected concentrations of both chloride and sodium during both challenges, but some were specific to either chloride (2 QTL) or sodium (1 QTL) concentrations. Six QTL (4 significant, 2 suggestive) affecting gill index were discovered. Two were specific to the trait, while the others were also identified as QTL for post-challenge ion concentrations. Altogether, allelic effects were consistent for QTL affecting chloride and sodium concentrations but inconsistent for QTL affecting ion concentrations and gill morphology. There was no systematic lineage effect (grand-parental origin of QTL alleles) on the recorded traits. CONCLUSIONS: For the first time, genomic loci associated with effects on major physiological components of osmotic adaptation to seawater in a nonmigratory fish were revealed. The results pave the way for further deciphering of the complex regulatory mechanisms underlying seawater adaptation and genes involved in osmoregulatory physiology in rainbow trout and other euryhaline fishes.

Genetic Parameters and Genome-Wide Association Studies of Quality Traits Characterised Using Imaging Technologies in Rainbow Trout, Oncorhynchus mykiss.

One of the top priorities of the aquaculture industry is the genetic improvement of economically important traits in fish, such as those related to processing and quality. However, the accuracy of genetic evaluations has been hindered by a lack of data on such traits from a sufficiently large population of animals. The objectives of this study were thus threefold: (i) to estimate genetic parameters of growth-, yield-, and quality-related traits in rainbow trout (Oncorhynchus mykiss) using three different phenotyping technologies [invasive and non-invasive: microwave-based, digital image analysis, and magnetic resonance imaging (MRI)], (ii) to detect quantitative trait loci (QTLs) associated with these traits, and (iii) to identify candidate genes present within these QTL regions. Our study collected data from 1,379 fish on growth, yield-related traits (body weight, condition coefficient, head yield, carcass yield, headless gutted carcass yield), and quality-related traits (total fat, percentage of fat in subcutaneous adipose tissue, percentage of fat in flesh, flesh colour); genotypic data were then obtained for all fish using the 57K SNP Axiom® Trout Genotyping array. Heritability estimates for most of the 14 traits examined were moderate to strong, varying from 0.12 to 0.67. Most traits were clearly polygenic, but our genome-wide association studies (GWASs) identified two genomic regions on chromosome 8 that explained up to 10% of the genetic variance (cumulative effects of two QTLs) for several traits (weight, condition coefficient, subcutaneous and total fat content, carcass and headless gutted carcass yields). For flesh colour traits, six QTLs explained 1-4% of the genetic variance. Within these regions, we identified several genes (htr1, gnpat, ephx1, bcmo1, and cyp2x) that have been implicated in adipogenesis or carotenoid metabolism, and thus represent good candidates for further functional validation. Finally, of the three techniques used for phenotyping, MRI demonstrated particular promise for measurements of fat content and distribution, while the digital image analysis-based approach was very useful in quantifying colour-related traits. This work provides new insights that may aid the development of commercial breeding programmes in rainbow trout, specifically with regard to the genetic improvement of yield and flesh-quality traits as well as the use of invasive and/or non-invasive technologies to predict such traits.

Sources of variation of DNA methylation in rainbow trout: combined effects of temperature and genetic background.

Phenotypic plasticity is a key component of the ability of organisms to respond to changing environmental conditions. In this study, we aimed to study the establishment of DNA methylation marks in response to an environmental stress in rainbow trout and to assess whether these marks depend on the genetic background. The environmental stress chosen here was temperature, a known induction factor of epigenetic marks in fish. To disentangle the role of epigenetic mechanisms such as DNA methylation in generating phenotypic variations, nine rainbow trout isogenic lines with no genetic variability within a line were used. For each line, half of the eggs were incubated at standard temperature (11°C) and the other half at high temperature (16°C), from eyed-stage to hatching. In order to gain a first insight into the establishment of DNA methylation marks in response to an early temperature regime (control 11°C vs. heated 16°C), we have studied the expression of 8 dnmt3 (DNA methyltransferase) genes, potentially involved in de novo methylation, and analysed global DNA methylation in the different rainbow trout isogenic lines using LUMA (LUminometric Methylation Assay). Finally, finer investigation of genome-wide methylation patterns was performed using EpiRADseq, a reduced-representation library approach based on the ddRADseq (Double Digest Restriction Associated DNA) protocol, for six rainbow trout isogenic lines. We have demonstrated that thermal history during embryonic development alters patterns of DNA methylation, but to a greater or lesser extent depending on the genetic background.

Combining Multiple Approaches and Models to Dissect the Genetic Architecture of Resistance to Infections in Fish.

Infectious diseases represent a major threat for the sustainable development of fish farming. Efficient vaccines are not available against all diseases, and growing antibiotics resistance limits the use of antimicrobial drugs in aquaculture. It is therefore important to understand the basis of fish natural resistance to infections to help genetic selection and to develop new approaches against infectious diseases. However, the identification of the main mechanisms determining the resistance or susceptibility of a host to a pathogenic microbe is challenging, integrating the complexity of the variation of host genetics, the variability of pathogens, and their capacity of fast evolution and adaptation. Multiple approaches have been used for this purpose: (i) genetic approaches, QTL (quantitative trait loci) mapping or GWAS (genome-wide association study) analysis, to dissect the genetic architecture of disease resistance, and (ii) transcriptomics and functional assays to link the genetic constitution of a fish to the molecular mechanisms involved in its interactions with pathogens. To date, many studies in a wide range of fish species have investigated the genetic determinism of resistance to many diseases using QTL mapping or GWAS analyses. A few of these studies pointed mainly toward adaptive mechanisms of resistance/susceptibility to infections; others pointed toward innate or intrinsic mechanisms. However, in the majority of studies, underlying mechanisms remain unknown. By comparing gene expression profiles between resistant and susceptible genetic backgrounds, transcriptomics studies have contributed to build a framework of gene pathways determining fish responsiveness to a number of pathogens. Adding functional assays to expression and genetic approaches has led to a better understanding of resistance mechanisms in some cases. The development of knock-out approaches will complement these analyses and help to validate putative candidate genes critical for resistance to infections. In this review, we highlight fish isogenic lines as a unique biological material to unravel the complexity of host response to different pathogens. In the future, combining multiple approaches will lead to a better understanding of the dynamics of interaction between the pathogen and the host immune response, and contribute to the identification of potential targets of selection for improved resistance.

Genetic determinism of spontaneous masculinisation in XX female rainbow trout: new insights using medium throughput genotyping and whole-genome sequencing.

Rainbow trout has a male heterogametic (XY) sex determination system controlled by a major sex-determining gene, sdY. Unexpectedly, a few phenotypically masculinised fish are regularly observed in all-female farmed trout stocks. To better understand the genetic determinism underlying spontaneous maleness in XX-rainbow trout, we recorded the phenotypic sex of 20,210 XX-rainbow trout from a French farm population at 10 and 15 months post-hatching. The overall masculinisation rate was 1.45%. We performed two genome-wide association studies (GWAS) on a subsample of 1139 individuals classified as females, intersex or males using either medium-throughput genotyping (31,811 SNPs) or whole-genome sequencing (WGS, 8.7 million SNPs). The genomic heritability of maleness ranged between 0.48 and 0.62 depending on the method and the number of SNPs used for the estimation. At the 31K SNPs level, we detected four QTL on three chromosomes (Omy1, Omy12 and Omy20). Using WGS information, we narrowed down the positions of the two QTL detected on Omy1 to 96 kb and 347 kb respectively, with the second QTL explaining up to 14% of the total genetic variance of maleness. Within this QTL, we detected three putative candidate genes, fgfa8, cyp17a1 and an uncharacterised protein (LOC110527930), which might be involved in spontaneous maleness of XX-female rainbow trout.

Resistance to a rhabdovirus (VHSV) in rainbow trout: identification of a major QTL related to innate mechanisms.

Health control is a major issue in animal breeding and a better knowledge of the genetic bases of resistance to diseases is needed in farm animals including fish. The detection of quantitative trait loci (QTL) will help uncovering the genetic architecture of important traits and understanding the mechanisms involved in resistance to pathogens. We report here the detection of QTL for resistance to Viral Haemorrhagic Septicaemia Virus (VHSV), a major threat for European aquaculture industry. Two induced mitogynogenetic doubled haploid F2 rainbow trout (Oncorhynchus mykiss) families were used. These families combined the genome of susceptible and resistant F0 breeders and contained only fully homozygous individuals. For phenotyping, fish survival after an immersion challenge with the virus was recorded, as well as in vitro virus replication on fin explants. A bidirectional selective genotyping strategy identified seven QTL associated to survival. One of those QTL was significant at the genome-wide level and largely explained both survival and viral replication in fin explants in the different families of the design (up to 65% and 49% of phenotypic variance explained respectively). These results evidence the key role of innate defence in resistance to the virus and pave the way for the identification of the gene(s) responsible for resistance. The identification of a major QTL also opens appealing perspectives for selective breeding of fish with improved resistance.