Front. Genet. Frontiers in Genetics Front. Genet. 1664-8021 Frontiers Media S.A. 1394656 10.3389/fgene.2024.1394656 Genetics Original Research Genome-wide association analysis of the resistance to infectious hematopoietic necrosis virus in two rainbow trout aquaculture lines confirms oligogenic architecture with several moderate effect quantitative trait loci Palti et al. 10.3389/fgene.2024.1394656 Palti Yniv 1 * Vallejo Roger L. 1 Purcell Maureen K. 2 † Gao Guangtu 1 Shewbridge Kristy L. 1 Long Roseanna L. 1 Setzke Christopher 3 Fragomeni Breno O. 4 Cheng Hao 5 Martin Kyle E. 6 Naish Kerry A. 3 1 National Center for Cool and Cold Water Aquaculture, USDA-ARS, Kearneysville, WV, United States 2 US Geological Survey, Western Fisheries Research Center, Seattle, WA, United States 3 School of Aquatic and Fishery Sciences, University of Washington, Seattle, WA, United States 4 Department of Animal Science, University of Connecticut, Storrs, CT, United States 5 Department of Animal Science, University of California, Davis, Davis, CA, United States 6 Hendrix Genetics, Sumner, WA, United States

Edited by: Gustavo Núñez-Acuña, University of Concepcion, Chile

Reviewed by: Clémence Fraslin, University of Edinburgh, United Kingdom

Marie Lillehammer, Norwegian Institute of Food, Fisheries and Aquaculture Research (Nofima), Norway

 *Correspondence: Yniv Palti, yniv.palti@usda.gov

Present address: Maureen K. Purcell, US Geological Survey, Forest and Rangeland Ecosystem Science Center, Corvallis, OR, United States

 24 05 2024 2024 15 1394656 01 03 2024 30 04 2024 Copyright © 2024 Palti, Vallejo, Purcell, Gao, Shewbridge, Long, Setzke, Fragomeni, Cheng, Martin and Naish. 2024 Palti, Vallejo, Purcell, Gao, Shewbridge, Long, Setzke, Fragomeni, Cheng, Martin and Naish

This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

Infectious hematopoietic necrosis (IHN) is a disease of salmonid fish that is caused by the IHN virus (IHNV), which can cause substantial mortality and economic losses in rainbow trout aquaculture and fisheries enhancement hatchery programs. In a previous study on a commercial rainbow trout breeding line that has undergone selection, we found that genetic resistance to IHNV is controlled by the oligogenic inheritance of several moderate and many small effect quantitative trait loci (QTL). Here we used genome wide association analyses in two different commercial aquaculture lines that were naïve to previous exposure to IHNV to determine whether QTL were shared across lines, and to investigate whether there were major effect loci that were still segregating in the naïve lines. A total of 1,859 and 1,768 offspring from two commercial aquaculture strains were phenotyped for resistance to IHNV and genotyped with the rainbow trout Axiom 57K SNP array. Moderate heritability values (0.15–0.25) were estimated. Two statistical methods were used for genome wide association analyses in the two populations. No major QTL were detected despite the naïve status of the two lines. Further, our analyses confirmed an oligogenic architecture for genetic resistance to IHNV in rainbow trout. Overall, 17 QTL with notable effect (≥1.9% of the additive genetic variance) were detected in at least one of the two rainbow trout lines with at least one of the two statistical methods. Five of those QTL were mapped to overlapping or adjacent chromosomal regions in both lines, suggesting that some loci may be shared across commercial lines. Although some of the loci detected in this GWAS merit further investigation to better understand the biological basis of IHNV disease resistance across populations, the overall genetic architecture of IHNV resistance in the two rainbow trout lines suggests that genomic selection may be a more effective strategy for genetic improvement in this trait.

 GWAS-genome-wide association study IHNV infection rainbow trout (oncorhynchus mykiss) QTL (loci of quantitative traits) aquaculture heritability section-at-acceptance Livestock Genomics Introduction

Infectious hematopoietic necrosis virus (IHNV; Novirhabdovirus salmonid) is endemic to western North America, and has been spread to parts of Asia, Europe, and the Middle East (Kurath et al., 2003; Enzmann et al., 2005; Adel et al., 2016; Wang et al., 2016). In the western U.S., IHNV causes morbidity and mortality in the commercial rainbow trout (Oncorhynchus mykiss) and Atlantic salmon (Salmo salar) aquaculture industries, as well as federal, state and tribal salmon and trout hatcheries (Kurath et al., 2003; Troyer and Kurath, 2003; Saksida, 2006; Bootland and Leong, 2009; Breyta et al., 2013). The disease caused by IHNV, infectious hematopoietic necrosis (IHN), is associated with acute mortality with losses as high as 90% (LaPatra, 1998). Rainbow trout fry are the most susceptible life stage but larger fish can still suffer disease (LaPatra et al., 1990). An effective IHNV DNA vaccine exists, but the need to inject the vaccine makes it impractical for immunizing large numbers of rainbow trout fry (Anderson et al., 1996; Lorenzen and LaPatra, 2005). Novel vaccine and delivery methods remain an active area of research, given the importance of IHNV to the trout industry (Plant and LaPatra, 2011; LaPatra et al., 2015; Larragoite et al., 2016). However, at present no cost-effective vaccine is commercially available for rainbow trout fry, and additional approaches are needed to control IHN disease in aquaculture. Genetic improvement for enhanced disease resistance represents such approach. The use of selectively bred resistant animals integrates well into an overall viral disease control strategy that incorporates good animal husbandry, biosecurity, use of specific pathogen free eggs and vaccination if practical, cost-effective and available.

Additive genetic variance is the basis for selective breeding in agricultural animals. Additive genetic variance for resistance to IHNV infection in rainbow trout has been estimated to be moderate based on heritability estimates for IHNV survival status (h2 = 0.23–0.55) and survival days (h2 = 0.02–0.20) in a steelhead trout (O. mykiss) population (Brieuc et al., 2015) and for survival status (h2 = 0.25) in the commercial breeding population of Clear Springs Food (CSF) (Vallejo et al., 2019). These results suggest that resistant rainbow trout lines can be developed using family-based selective breeding methods. Furthermore, selective breeding of a rainbow trout line for resistance to IHNV has been implemented at the CSF breeding program since the year 2000 (Campbell et al., 2014), and recently we have shown that genomic-enabled selection can accelerate the improvement of disease resistance in that commercial rainbow trout strain (Vallejo et al., 2020). Here, we examine the opportunity for similar practises in two commercial aquaculture lines that were naïve to previous exposure to IHNV.

Genomic technologies have improved breeding predictions accuracy in agriculture by identifying DNA markers linked to complex phenotypic traits (Barabaschi et al., 2016). Genomic selection (GS) is a selective breeding strategy that examines together the association between all genetic markers genotypes in a population with the trait or traits of interest to predict the breeding value of an individual animal from the population (Meuwissen et al., 2001; Goddard et al., 2011). For example, in rainbow trout aquaculture it was shown that GS can double the accuracy of breeding value predictions for resistance to bacterial cold-water disease (Vallejo et al., 2017a), and in recent years the technology has been widely adopted by the salmonids aquaculture industry and in other aquaculture species (Song et al., 2023; Yáñez et al., 2023). Marker assisted selection (MAS) is a simplified strategy of genomic selection that can be implemented for traits in which most of the genetic variance is controlled a single locus. For example, MAS for infectious pancreatic necrosis (IPN) resistance in Atlantic salmon has resulted in a 75% reduction in the occurrence of that viral disease in the European Atlantic salmon industry (Moen et al., 2007; Houston et al., 2008; Houston et al., 2010; Moen et al., 2015). A key first step in conducting such studies is to identify molecular markers linked to the trait by studying their co-inheritance in segregating populations, using analytical approaches such as genome-wide association studies (GWAS) (Goldstein et al., 2010; Visscher et al., 2012; Fernando et al., 2014; Misztal et al., 2014; Vallejo et al., 2017b; Silva et al., 2019; Vallejo et al., 2019). Such studies are more powerful if they are conducted across populations, because the overall number and overlapping genomic regions identified across populations indicate the potential for genomic selection program and marker assisted selection across different broodstocks. These approaches are particularly effective when attempting to initiate a breeding program for a trait in a population that has not undergone previous selection for that trait.

Troutlodge Inc. maintains a year-round production of eggs though the use of four distinct broodstock populations, with peak spawning in February, May, August and November. The two-year spawning cycle has separated these lines further into even and odd year groups. Phylogenetic analyses based on genotypes from 96 Fluidigm SNP assays (Liu et al., 2016) indicated that all four lines are genetically distinct (Liu et al., 2017). Here, we studied the May 2019 and November 2018 populations because they represent distinct genotypic differences within the breeding program. Identification of conserved regions between these two lines may indicate a higher likelihood of identifying these loci in other commercial trout lines. While these lines have undergone selection for multiple production traits, there has been no intentional selection on IHNV resistance. Therefore, these lines are considered naïve for IHNV and may be polymorphic for genes with large effects on IHNV resistance. The CSF breeding population that has already been subject to IHNV survival selection may be fixed for large effect loci, making it more challenging to identify large effect QTL (Campbell et al., 2014; Vallejo et al., 2019).

Through genome-wide association analyses in the two commercial rainbow trout lines, the objectives of this study were to (1) determine whether the opportunity for genomic selection is consistent across different broodstock lines by comparing the genetic architecture of the trait in the distinct May and November lines; (2) identify whether there were shared genomic regions across lines that are significantly associated with resistance to IHNV, to assess the potential for improving unrelated lines through marker-assisted selection; and (3) identify candidate genes within the QTL regions that will enable future functional studies aimed at identifying the host genes involved in IHNV resistance, which in turn can accelerate genetic improvement through selective breeding and development of more effective treatments or vaccination models.

 Materials and methods Ethics statement

All fish work was conducted in accordance with national and international guidelines. The protocol for this study was approved by the Institutional Animal Care and Use Committee of the University of Washington, Seattle, WA (Protocol number 4456-01).

 Fish rearing and IHNV challenge

In February 2019, rainbow trout fry (age of ∼80 days post hatching) from 100 nucleus families from the Troutlodge, Inc. November 2018 spawning population (TLUN 2018) were transported from the Sumner, Washington hatchery to the fish rearing facility of the United States Geological Survey (USGS) Western Fisheries Research Center (WFRC) in Seattle, Washington. The 100 full-sib (FS) families were generated from 100 dams to 62 sires, and they included 37 paternal half-sib (HS) families. The fish were reared at the USGS facility for 90 days until the disease challenge with IHNV. Similarly, in August 2019 fry from 103 nucleus families from the Troutlodge May 2019 spawning population (TLUM 2019) were transported to the WFRC facility and reared there for 30 days prior to the disease challenge. The 103 FS TLUM2019 families were generated using 75 sires and 103 dams, and they included 28 paternal half-sib families.

For the TLUN2018 population, the fish were approximately 5.0 g in size and ∼170 days post hatching at the start of the disease challenge, and ∼21 fish per family were placed in three replicated tanks (final count N = 2,025 total) and immersed in IHNV (strain 220-90) using standard methods (Brieuc et al., 2015). A mock control tank with additional 102 fish was immersed in virus-free cell culture medium and tank water volumes were adjusted to provide a similar biological density as the 3 challenge tanks. IHNV was cultured in the fathead minnow EPC cell line (Winton et al., 2010) at 15°C and viral amount was determined by plaque assay (Batts and Winton, 1989). A pilot trial was conducted to determine optimal dosage that was expected to cause ∼50% mortality over 21 days. Confirmed IHNV free trout fry were immersed in a dose of 2,000 plaque forming units (PFU) per ml in static water with aeration for 1 h. Total water volume in the challenge tanks was kept at ∼153 L with fish density of ∼23 g/L, and for the static immersion challenge the volume was dropped to ∼80 L (∼44 g/L). After 1 h, the water flow was resumed, and the challenge tanks were monitored for moribund fish or mortalities daily for a 21-day period and all survivors were euthanized at 21 days. All experiments were conducted at water temperature of ∼15°C to mimic water temperatures common in the rainbow trout aquaculture industry in Southern Idaho.

For the TLUM2019 population, the fish were approximately 2.0 g in size and ∼110 days post hatching at the start of the disease challenge, and 30 fish per family were placed in three replicated tanks (final count N = 3,049 total) and immersed in IHNV (strain 220-90) water. Challenge conditions were as described above, but the rainbow trout fry were immersed in a dose of 1,000 PFU per ml based on finding of slightly higher susceptibility in the pre-challenge pilot study. Total water volumes were similar to the TLUN2018 population challenge, with a fish density of ∼13 g/L for acclimation and survival trial and ∼25 g/L for the 1-h static immersion challenge.

One main difference in this study compared to Brieuc et al., 2015 was the pooling of individuals from all families and the nearly even distributing the pooled families across three ‘common garden’ replicated tanks. In contrast to individual family tanks where the fish from susceptible families are exposed to higher viral load shedding from their siblings, the ‘common garden’ experimental design has the advantage that all individuals are exposed to the same infection pressure regardless of family. Fin tissues were excised from the parents of the nucleus families at spawning and from all virus-challenged individuals and preserved in 95% ethanol for DNA analysis. After the challenge, identification of an individual to family was performed by genotype analyses. A schematic illustration of the experimental design of the disease challenges is presented in Figure 1.

Schematic representation of the experimental design for disease challenges. The Troutlodge November line is represented in this study by year class TLUN2018 and the May line is represented by TLUM 2019.

 IHNV resistance phenotypes

Two phenotypes were recorded; survival days (days to Death; DAYS) and status at the end of the challenge (alive or dead at 21 days; STATUS). However, only the survival days phenotype records were used in the data analyses because survival rate was lower than 10% in the two challenges. Due to the low phenotypic variation in the survival status, and resulting low heritability, this phenotype was excluded from further data analyses.

 Fish populations used for GWAS

TLUN 2018: Eight generations including 7,142 individuals with pedigree records were available for this population. The fish genotyped represent 100 FS families, including 1,859 offspring that were sampled from the three challenge tanks. Due to logistic limitations we had to exclude 150 fish (∼8%) from genotyping. The 50 fish that were excluded from genotyping and pedigree assignments per tank were from the peak mortality days that also represented the median distribution of the mortality per day curve. This exclusion of ∼8% of the fish from the middle of the phenotypic distribution does represent small selective genotyping bias, but we believe that the effect on over estimation of variance components must have been very small.

TLUM 2019: Only two generations of pedigree records (parents and offspring) were available for this population which was used to combine parents from the May spawning odd and even years. The fish genotyped represent 103 full-sib (FS) families, including all the offspring (N = 1,768) that were sampled from the two tanks with the lower mortality rate.

 SNP genotyping and pedigree assignments of offspring

The DNA samples from the two populations were genotyped with the Rainbow Trout Axiom 57K SNP array (Thermo Fisher) following previously described procedures (Vallejo et al., 2019). The samples were genotyped by a commercial service provider (Center for Aquaculture Technologies, San Diego, CA). Following genotype calls, all the SNP markers and samples with call rate (CR) greater than 95% were retained for pedigree assignments. Parentage assignment of individual offspring to all known parental pairs used to produce the full-sib families was conducted using a customized Perl script (Supplementary Material S1). Assignment required at least a 99% Mendelian segregation match between each parent pair and offspring trio.

Before conducting GWAS, the raw marker genotype dataset was filtered using the software BLUPF90 (Misztal et al., 2015). The QC retained samples and SNPs with a genotype calling rate ≥0.97, minor allele frequency ≥0.05, and departures from Hardy-Weinberg equilibrium ≤0.15. Parent-progeny pairs were tested again with the BLUPF90 script for Mendelian error rate (MER). Samples and SNPs were discarded from further analysis if they had MER ≥1% and ≥3%, respectively. Next, we determined the physical map position (Genome Accession GCA_013265735.3) (Gao et al., 2021) of each of the QC filtered markers. A small fraction of the filtered markers (<3%) that were not mapped to chromosome sequences did not have a physical map location, and those markers were excluded from further data analysis. The maps generated for the two populations can be found in Supplementary Material S2, S3.

After genotype data QC and filtering of animals and markers in the TLUN2018 population, a total of 1,796 offspring and 162 parents were used in the association analysis with genotype data from 33,715 informative SNPs. In the TLUM2019 population, 1,742 offspring and 178 parents were successfully genotyped with 35,135 informative SNPs after QC and filtering of animals and markers.

 Estimation of genetic variance parameters

Three methods were used to estimate genetic variance parameters, two of which were based on genetic markers and one used only pedigree and phenotype records.

The DAYS records from TLUN 2018 n = 1796 and TLUM 2019 n = 1742 were fit to an animal linear model to estimate genetic variance parameters for survival DAYS in each population dataset, separately. We fitted a linear mixed model for DAYS records using this animal model: y = 1 μ + X d d + Za + Wc + e , where y is the vector of phenotypic records, 1 is a vector of 1s, μ is the overall mean of phenotypic records, d is a vector of fixed effects, X d is an incidence matrix relating records to fixed effects in d , a is a vector of random individual animal effects, c is a vector of random common environmental effects (i.e., nested effect of families within challenge tank), e is a vector of residual effects, and Z and W are incidence matrices relating records to random animal and random common environmental effects in a and c , respectively.

For the single-step genomic best linear unbiased prediction method (ssGBLUP), the variances of a , c and e are: v a r a c e = H σ a 2 0 0 0 I σ c 2 0 0 0 I σ e 2 where σ a 2 , σ c 2 and σ e 2 are additive genetic, common environment and residual variances, respectively, and H is a matrix that combines pedigree ( A ) and genomic ( G ) relationship matrices as in Aguilar et al. (Aguilar et al., 2010), and the variances of a , c and e are estimated by replacing H with A matrix. We used two and three tanks in the TLUM2019 and TLUN2018 populations, respectively, and the effect of tank on DAYS was significant (p < 0.05) only in TLUM2019 but not in TLUN 2018. Thus, the effect of tank was included only in the vector of fixed effects d for TLUM 2019.

The genetic parameters were estimated using the pedigree-based BLUP (PBLUP) and PBLUP with genomic data using single-step method (ssGBLUP) under a Bayesian framework, using the software gibbsf90+ from the computer application BLUPF90 (Misztal et al., 2015). Our Gibbs sampling scheme included 250,000 iterations, of which the first 50,000 iterations were discarded; from the remaining 200,000 iterations one sample was saved from every 40 iterations, such that results from 5000 samples were used in the genetic analysis. The proper mixing and convergence of the Markov chain Monte Carlo (MCMC) Gibbs sampling approach was evaluated with the R package CODA (Plummer et al., 2006).

The heritability for DAYS was estimated as: h 2 = σ a 2 / σ a 2 + σ c 2 + σ e 2 , where h 2 is the estimated narrow-sense heritability, σ a 2 is additive genetic variance, σ c 2 is variance due to nested effect of families within challenge tank (i.e., common environment effect) and σ e 2 is residual error variance.

We also estimated genetic variance parameters using Bayesian multiple regression with BayesB model (BMR-BayesB) with software JWAS (Cheng et al., 2018) run with option ‘single-step = false’. We fitted a linear mixed model for DAYS records using this animal model: y = 1 μ + X d d + Xb + Za + Wc + e , where X is an n × k matrix of observed genotype covariates for k total number of SNPs across the genome for genotyped n individuals, b is a vector of k additive SNP effects, a is a vector of random polygenic effects. The remaining linear model terms were already defined in the analyses based on PBLUP and ssGBLUP, as described earlier. Scaled inverse chi-squared distributions were used for genetic variance and residual variance as described by Fernando et al., 2014; Fernando et al., 2014); in these priors, the degree of freedom was four and the scaled parameters were estimated by assuming that the proportion of variance of the phenotypic data explained by the regression is 0.5.

The BMR-BayesB method fits a mixture model to estimate marker effects, which assumes that there are two types of SNPs: a fraction of SNPs with non-zero effects 1 − π that are drawn from distributions with a marker-specific variance σ α 2 , and another known fraction of SNPs π that a-priori have zero effect on the quantitative trait (Meuwissen et al., 2001). In our study, the mixture parameter π was assumed to be known and defined to meet the condition k ≤ n ; where n is the number of fish with genotype records, p is the effective number of SNPs, and k = 1 − π p is the number of markers sampled as having a non-zero effect that are fitted simultaneously in the Bayesian multiple regression model (Garrick and Fernando, 2013). In this study, we used π = 0.999 which enabled sampling about 34 and 35 non-zero effect SNPs and fitted in the multiple regression model at each MCMC iteration in the BMR-BayesB analysis performed with TLUN2018 and TLUM2019 datasets, respectively. The MCMC Gibbs sampling scheme and the assessment of correct mixing and convergence of the MCMC iterations were like those described in the section of estimation of genetic parameters with PBLUP and ssGBLUP.

 GWAS with wssGBLUP

The November and May spawning lines are two separate and distinct genetic lines with no shared pedigree. Therefore, the GWAS was performed separately for the two populations using DAYS phenotype and genotype data records from TLUM2019 and TLUN 2018. We conducted GWAS with the wssGBLUP method using 1-Mb sliding SNP windows (Wang et al., 2012; Misztal et al., 2015). Briefly, the effects were calculated for individual SNPs in the first step, as shown below. Afterward, the effects of all SNPs within a 1-Mb distance were added and recorded for each sliding window. Briefly, the 1-Mb window slides by one SNP at a time from the first SNP until the last SNP on each chromosome and the results for SNPs that are included in the window are jointly summarized; thus, the estimates for SNP effects is a moving average of n adjacent SNPs included in the 1-Mb window (Misztal et al., 2015). The choice of a 1-Mb window size was based on our recent estimate of strong LD (r 2 ≥ 0.25) spanning on average over 1 Mb in the rainbow trout genome (Vallejo et al., 2018).

In GWAS with wssGBLUP, the weights for each SNP are one for the first iteration, which indicates that all SNPs have the same weight (i.e., single-step GBLUP). For the subsequent iterations (2nd, 3rd, etc.), the weights are SNP-specific variances that are calculated using the estimate of the SNP allele-substitution effect from the previous iteration and the corresponding SNP allele frequencies (Wang et al., 2012). The estimates of SNP effects were calculated using a weighted relationship matrix, using the following equation: u ^ = DM ′ MDM ′ − 1 a ^ g , where u ^ is the vector of the estimated SNP effects; D is a diagonal matrix of weights for variances of SNP effects; M is a matrix linking genotypes of each SNP to each individual; and a ^ g is the estimate of the additive genetic effect for genotyped animals. The individual variance of SNP effects, which corresponds to the diagonal elements of D , was estimated as suggested by Zhang et al. (Zhang et al., 2021): σ ^ u , i 2 = u ^ i 2 2 p i 1 − p i , where: u ^ i 2 is the square of the effect at SNP i , and p i is the observed allele frequency for the second allele of SNP i . In this GWAS, we used results from the second iteration of wssGBLUP, because generally they provide the highest accuracy of genomic predictions (Vallejo et al., 2016) and marker effects (Wang et al., 2012; Irano et al., 2016; Melo et al., 2016; Vallejo et al., 2016). The linear model used for GWAS followed that of the ssGBLUP model previously described for the estimation of genetic variance parameters. It was conducted using Gibbs sampling methods implemented in the software gibbsf90+ from the computer application BLUPF90 (Misztal et al., 2015).

 GWAS with bayesian multiple regression

We conducted GWAS for DAYS with BMR-BayesB model using 1-Mb non-overlapping SNP windows (Fernando et al., 2014; Fernando et al., 2016). The BMR-BayesB model uses the same pedigree information and all animals that had phenotype and genotype records, following the wssGBLUP method. The 1-Mb window’s posterior probability of association (WPPA) with the phenotype was used to estimate the window’s proportion of false positive as P F P = 1 − W P P A (Fernando et al., 2017).

The GWAS for DAYS was performed with BMR-BayesB method (Cheng et al., 2018) implemented in the software JWAS (Cheng et al., 2018). The linear model used for GWAS conformed with the BMR-BayesB model previously described in the section on estimation of genetic variance parameters. The BMR-BayesB uses MCMC Gibbs sampling in the GWAS analysis (Garrick and Fernando, 2013). The MCMC Gibbs sampling scheme and the assessment of correct mixing and convergence of the MCMC iterations have also been previously described (section on estimation of genetic parameters). We did not perform GWAS using sliding SNP windows with BMR-BayesB because this method has not been implemented in the software JWAS.

 Detection of QTL

Quantitative trait loci (QTL) associated with resistance to IHNV were defined as 1-MB SNP windows that explained additive genetic variance (AGV) higher than 1.9% (Vallejo et al., 2022). All QTL windows that were mapped to the same chromosome and were located within less than 20 Mb from each other were defined as belonging to the same QTL region. We have previously used these conservative criteria for defining QTL in GWA analyses in rainbow trout to reduce the type I error rate as much as possible.

For comparison of the QTL genome positions with those previously identified in the Clear Springs Foods (CSF) line (Vallejo et al., 2019), the positions of the flanking SNPs from each QTL in that study were identified on the current version of the rainbow trout reference genome (Gao et al., 2021).

 Identification of protein coding genes in the QTL regions

Standard template queries were used on AquaMine v1.2 (https://aquamine.rnet.missouri.edu/aquamine/) to identify protein coding genes in the most significant QTL SNP windows and to search for the genes description, the GO terms and predicted pathways.

 Results IHNV survival phenotypes

The mortality rate in both challenges was higher than predicted based on the pilot dosing experiments. In the TLUN2018 population the average mortality per tank was 86% with SD of 2% compared to 2% mortality in the mock trial control tank. In the TLUM2019 population the average mortality per tank was 92% with SD of 11% compared to no mortality in the mock trial. Daily mortality peaked on day 5 in the TLUN2018 challenge and on day 7 in the TLUM2019 challenge (Figure 2).

Distribution of daily mortality post IHNV immersion challenge in the commercial rainbow trout TLUN2018 population (A) and in the TLUM2019 population (B). The bars represent standard deviation (SD) between tanks.

 Heritability of resistance to IHNV

The estimated heritability for survival days had a range of 0.08–0.25 and 0.15–0.23 across three methods of estimation in TLUM2019 and TLUN 2018, respectively (Table 1). These estimates of low to moderate heritability for the phenotype of survival days underline the potential for genetic improvement of IHNV resistance through selective breeding in these two rainbow trout breeding populations. Due to the very low survival rate in both experiments, heritability for the binary trait of survival on day 21 post challenge was not different from zero (data not shown), and this trait was excluded from genome-wide association analysis.

Genetic parameter estimates for IHNV resistance (survival days, DAYS) in two rainbow trout aquaculture strains.

Strain a Method b Genetic parameter c
σ g 2 σ c 2 σ e 2 h 2
TLUM2019 PBLUP 1.97 ± 1.02 1.79 ± 0.53 19.45 ± 0.88 0.08 (±0.04)
ssGBLUP 3.89 ± 1.04 1.45 ± 0.46 18.34 ± 0.83 0.16 (±0.04)
BMR-BayesB 5.36 ± 0.76 1.60 ± 0.36 16.47 ± 0.80 0.25 (±0.03)
TLUN2018 PBLUP 8.75 ± 2.17 0.63 ± 0.46 28.21 ± 1.59 0.23 (±0.05)
ssGBLUP 5.55 ± 1.29 1.01 ± 0.56 30.43 ± 1.21 0.15 (±0.03)
BMR-BayesB 8.22 ± 1.22 1.86 ± 0.45 27.57 ± 1.36 0.23 (±0.03)

Variance components analysis was conducted using fish from the year-class (YC) 2019 of the Troutlodge May-spawning strain (TLUM, 2019) and from YC, 2018 of the November-spawning strain (TLUN, 2018).

Genetic parameters were estimated using pedigree-based BLUP (PBLUP), single-step GBLUP (ssGBLUP), and Bayesian multi regression with Bayes B (BMR-BayesB) methods. The survival days records were analyzed using animal linear models.

Genetic parameter estimate (± standard error): σ g 2 is the additive genetic variance; σ c 2 is the variance due to the nested effects of families within challenge tanks; σ e 2 is the residual error variance; and h 2 is the estimate of narrow-sense heritability for survival days.

 QTL associated with resistance to IHNV

In the TLUM2019 population GWAS, we detected 14 SNP windows of 1-Mb that were associated with resistance to IHNV (AGV≥1.9%) with either the wssGBLUP or the BMR-BayesB GWAS methods (Table 2). The 14 QTL windows were determined to be located within 11 chromosomal regions. Eight of those QTL windows were detected with wssGBLUP and six were detected with BMR-BayesB. Jointly, the 11 QTL regions explained up to 49.7% of the additive genetic variance when accounting only for the largest effect SNPs window detected in each QTL region. As shown in Figure 3, the QTL windows on chromosomes 1 (AGV = 11.4%), 7 (AGV = 7.0%) and 23 (AGV = 7.6%) had the highest proportions of explained additive genetic variance. However, only one QTL window, the window on chromosome 1, had a false positive probability (PFP) lower than 0.05.

Characteristics of quantitative trait loci associated with resistance to IHNV in the TLUM2019 aquaculture strain of rainbow trout a .

Chr b AGV (%) c PFP d GWAS method e Physical map (bp) f Window flanking SNPs SNPs per window
Start End Start End
1 11.4 0.03 BMR-BayesB 48,049,099 48,967,892 Affx-88955618 Affx-88930703 17
1 3.5 N/A wssGBLUP 48,632,135 49,623,888 Affx-88952382 Affx-88938488 24
3 3.6 N/A wssGBLUP 51,430,431 52,429,704 Affx-88933807 Affx-88905010 25
5 4.0 0.38 BMR-BayesB 19,009,057 19,953,516 Affx-88906315 Affx-88940046 23
6 1.9 N/A wssGBLUP 25,083,663 26,066,808 Affx-88919141 Affx-88953291 29
7 7.0 0.26 BMR-BayesB 21,017,550 21,868,247 Affx-88941998 Affx-88958253 23
11 3.9 0.44 BMR-BayesB 57,030,022 57,990,833 Affx-88911746 Affx-88928288 20
12 2.3 N/A wssGBLUP 66,510,299 67,487,681 Affx-88913164 Affx-88923715 27
19 3.1 N/A wssGBLUP 15,012,572 15,954,248 Affx-88940815 Affx-88911258 15
23 7.6 0.09 BMR-BayesB 13,029,157 13,930,763 Affx-88914566 Affx-88906601 23
23 2.8 N/A wssGBLUP 28,704,800 29,704,625 Affx-88958162 Affx-88907862 25
30 3.0 0.54 BMR-BayesB 14,136,821 14,976,738 Affx-88904384 Affx-88943429 8
30 2.8 N/A wssGBLUP 20,284,236 21,267,747 Affx-88907684 Affx-88904888 27
32 1.9 N/A wssGBLUP 33,372,622 34,329,512 Affx-88926444 Affx-88945805 31

GWAS, performed using fish from year-class 2019 families from the Troutlodge May-spawning nucleus breeding population.

Chromosome numbers are based on the USDA_OmykA_1.1 genome assembly (Gao et al., 2021; GenBank Assembly Accession GCA_013265735.3).

Explained additive genetic variance (AGV) by tested 1-Mb window (%). The 1-Mb windows with AGV≥1.9% were defined as associated with IHNV, resistance.

PFP is the probability of false positives defined as P F P = 1 − W P P A , where WPPA is the window posterior probability of association. WPPA was estimated only with the BMR-BayesB model.

GWAS, conducted using weighted single-step GBLUP (wssGBLUP) and Bayesian multiple regression with BayesB (BMR-BayesB) methods.

SNP, positions in base pairs (bp) based on rainbow trout reference genome sequence (Gao et al., 2021; GenBank Assembly Accession GCA_013265735.3).

Manhattan plots showing the association of 1-Mb SNP windows with resistance to IHNV in year-class 2019 families from the May-spawning nucleus breeding population (TLUM) using two GWAS models: (A) Weighted singe-step GBLUP (wssGBLUP); and (B) Bayesian multiple regression with BayesB (BMR-BayesB).

In the TLUN18 population, we detected 13 SNP windows of 1-Mb in 12 chromosomal regions that were associated with resistance to IHNV (Table 3). From these windows, six were detected with wssGBLUP and seven were detected with BMR-BayesB. The 12 QTL regions jointly explained up to 69.23% of the additive genetic variance when accounting only for the largest effect SNP window detected in each QTL region. As shown in Figure 4, the QTL windows on chromosomes 13 (AGV = 13.2%), 21 (AGV = 7.9%) and 29 (AGV = 12.2%) had the highest proportions of explained additive genetic variance. The windows detected with BMR-BayesB on chromosomes 13 and 29 had a false positive probability (PFP) lower than 0.05.

Characteristics of quantitative trait loci associated with resistance to IHNV in the TLUN2018 aquaculture strain of rainbow trout a .

Chr b AGV (%) c PFP d GWAS method e Physical map (bp) f Window flanking SNP SNPs per window
Start End Start End
1 5.9 0.34 BMR-BayesB 37,024,181 37,985,996 Affx-88940706 Affx-88919761 19
3 2.6 N/A wssGBLUP 51,104,028 52,087,029 Affx-88955207 Affx-88951474 26
5 7.0 0.23 BMR-BayesB 43,002,688 43,987,163 Affx-88954192 Affx-88917739 19
6 2.6 0.57 BMR-BayesB 25,050,389 25,985,630 Affx-88944756 Affx-88920694 31
8 2.2 0.68 BMR-BayesB 53,052,503 53,969,345 Affx-88937384 Affx-88925776 28
10 5.5 N/A wssGBLUP 35,810,116 36,734,512 Affx-88928023 Affx-88930040 40
13 13.2 0.02 BMR-BayesB 41,096,300 41,961,648 Affx-88930116 Affx-88942374 16
17 5.4 N/A wssGBLUP 26,575,732 27,437,130 Affx-88955005 Affx-88913439 21
19 2.8 0.63 BMR-BayesB 59,071,937 59,743,776 Affx-88919240 Affx-88958570 17
21 7.9 N/A wssGBLUP 25,853,803 26,848,204 Affx-88919092 Affx-88925348 34
23 1.9 N/A wssGBLUP 15,680,409 16,648,321 Affx-88911729 Affx-88934547 25
29 (Y) 2.6 N/A wssGBLUP 4,884,989 5,764,531 Affx-88912960 Affx-88948694 13
29 (Y) 12.2 0.03 BMR-BayesB 6,021,854 6,963,833 Affx-88953281 Affx-88950556 11

GWAS, performed using fish from year-class 2018 families from the Troutlodge November-spawning nucleus breeding population.

Chromosome numbers are based on the USDA_OmykA_1.1 genome assembly (Gao et al., 2021; GenBank Assembly Accession GCA_013265735.3).

Explained additive genetic variance (AGV) by tested 1-Mb window (%). The 1-Mb windows with AGV≥1.9% were defined as associated with IHNV, resistance.

PFP is the probability of false positives defined as P F P = 1 − W P P A , where WPPA is the window posterior probability of association. WPPA was estimated only with BMR-BayesB model.

GWAS, conducted using weighted single-step GBLUP (wssGBLUP) and Bayesian multiple regression with BayesB (BMR-BayesB) methods.

SNP, positions in base pairs (bp) based on rainbow trout reference genome sequence (Gao et al., 2021; GenBank Assembly Accession GCA_013265735.3).

Manhattan plots showing the association of 1-Mb SNP windows with resistance to IHNV in year-class 2018 families from the November-spawning nucleus breeding population (TLUN) using two GWAS models: (A) Weighted singe-step GBLUP (wssGBLUP); and (B) Bayesian multiple regression with BayesB (BMR-BayesB).

Two 1-Mb QTL windows were detected within chromosomes 1, 23 and 30 in the TLUM2019 population (Table 2) and chromosome 29 in the TLUN2018 population (Table 3). In each of those four cases we defined the two windows to be part of the same QTL region since those pair of windows were detected within a region smaller than 20 Mb. It is also important to note that within each of those four pairs of 1-MB QTL windows each window was detected by a different GWAS detection algorithm, which further support the findings of true QTL in those four chromosomal regions.

 Predicted protein coding genes identified in QTL regions

Gene content from the most significant QTL was identified by exploring protein coding genes that were annotated by the NCBI Refseq and by the EBI Ensembl. We confined our analyses to the QTL windows that were detected in TLUM2019 on chromosomes 1 (PFP<0.05) and 23 (detected by two GWAS methods and PFP<0.1), and in TLUN2018 on chromosomes 13 (PFP<0.05), 21 (strongest QTL detected with the wssGBLUP method) and chromosome Y or 29 (PFP<0.05). The list of annotated genes from those five QTL regions is provided in Supplementary Material S4, including the start and end positions on the chromosome, and where available, the gene description and the metabolic pathways that were linked to those genes in Refseq. The list of all the Gene Ontology (GO) terms that were assigned to the genes from each of the five QTL regions is presented in Supplementary Material S5.

 Discussion

Overall, we identified oligogenic structure for resistance to IHNV in the two commercial rainbow trout breeding populations that were evaluated in this study, with several moderate-effect QTLs accounting for large portion of the total additive genetic variance (AGV). In the May line (TLUM 2019) we detected 11 QTL regions that jointly explained nearly 50% of the AGV, and in the November line (TLUN 2018) we detected 12 QTL regions that jointly explained nearly 70% of the AGV. Only the QTL 1-Mb windows detected in TLUN2018 on chromosome 21 (25.8-26.8 Mb) and TLUM2019 on chromosome 30 (14.1-14.9 Mb) in this study are co-localized near QTL regions that were detected in previously in a third aquaculture breeding population, in which selection for resistance has occurred (Supplementary Material S6) (Vallejo et al., 2019). Due to the lack of major QTL for this trait in the commercial rainbow trout populations that were analysed, it is likely that genomic enabled selection, rather than marker assisted selection, is the more effective approach for improving the accuracy of breeding predictions for IHNV disease resistance. A similar conclusion was reached in our previous study on the third aquaculture breeding population (Vallejo et al., 2020). Therefore, our hypothesis that a major QTL for resistance to IHNV may be found segregating in the two aquaculture breeding populations - since they do not have history of exposure or selective breeding for IHNV - has not been supported by the results of this study.

We identified overlapping QTL regions located on chromosomes 3 and 6 between the two populations and neighboring QTL regions that may be in linkage disequilibrium (LD) with the same causative gene or genes on chromosomes 1, 5 and 23 (Figure 5). Rainbow trout aquaculture strains have been shown to have long-range LD, particularly on chromosome 5 (Vallejo et al., 2018; Pearse et al., 2019; Vallejo et al., 2019), and thus we cannot rule out the possibility that QTL windows that are even within 20 Mb from each other are in LD with the same causative gene or genes. QTL that are detected in multiple populations are more likely to be true QTL and are also of greater interest for further investigation. The QTL from chromosomes 1 and 23 were also the most significant QTL that were detected in TLUM 2019, making them strong candidates for further evaluation and analyses. Seven QTL located on chromosomes 8, 10, 13, 17, 19, 21 and 29 were unique to the TLUN2018 population and six QTL located on chromosomes 7, 11, 12, 19, 30 and 32 were unique to TLUM 2019. Interestingly, the most significant QTL detected in TLUN2018 from chromosomes 13, 21 and 29 were not shared with TLUM 2019. In comparison with the Clear Springs Food (CSF) population QTL detected in our previous study (Vallejo et al., 2019) we determined that the TLUN2018 QTL on chromosomes 21 and the TLUM2019 QTL on chromosome 30 are co-localized with the CSF QTL that were detected on the same chromosomes (Supplementary Material S6), making QTL-21 that was also one of the strongest QTL in the TLUN2018 population another top candidate for further evaluation.

Schematic illustration of chromosome positions of co-localized and unique QTL windows associated with resistance to IHNV (AGV ≥1.9%) in the rainbow trout November 2018 and May 2019 spawning nucleus breeding populations.

The age and size of the fish in the two disease challenge experiments was substantially different. The fish from the TLUN2018 population were challenged at the age of ∼170 days post hatching (dph) with an average weight of 5 g and from the TLUM2019 population at the age of ∼110 dph with an average weight of 2 g. Rainbow trout were previously found to remain susceptible to IHNV over a range of sizes but show decreasing susceptibility with increasing size (LaPatra et al., 1990; LaPatra, 1998). For example, rainbow trout ranging from 1.7, 3.4 and 7.4 g exposed to IHNV (1 × 105 pfu/mL) had cumulative percent mortality of 55%, 44% and 18%, respectively (LaPatra, 1998). Previous studies implicate defence barriers or early innate responses as being important for genetic differences in IHNV susceptibility (Quillet et al., 2007; Purcell et al., 2010; Bledsoe et al., 2022). Adaptive immune response such as neutralizing antibodies is important for long-term immunity to IHNV but is unlikely to contribute to survival of naïve fish in an acute challenge (Lorenzen and Lapatra, 1999). Given that the mechanism underlying the impact of age and size on acute IHNV susceptibility has not been characterized, we cannot rule out differential immune responses. However, IHNV susceptibility at two different sizes was found to be a stable trait in both clonal rainbow trout lines (Quillet et al., 2007) and full-sibling rainbow trout families (Purcell et al., 2010). Therefore, while overall susceptibility may change with size, the line or family susceptibility was not found in previous studies to have a transitory effect related to size or age. Taking all this information together, we cannot rule out that the meaningful difference in age and size had an impact on the trait of resistance to IHNV between the two disease challenges in this study and might have contributed in part to the detection of different QTL between the two populations.

The heritability estimates for IHNV resistance in our study were of low to moderate magnitude (0.08-0.25) and somewhat lower than those reported previously in two different rainbow trout population (Brieuc et al., 2015; Vallejo et al., 2019). The different genetic background of those populations may be a factor in the differences in heritability estimates. In addition, one notable difference from those two previous studies was that here we used a common garden design for the disease challenge compared to individual family tanks that were used in the previous studies. Hence, in the current study design the effect of the common shared environment on the family-based response to IHNV infection and survival was reduced, which might have contributed to less inflated heritability estimates compared to the previous studies. Still, the estimated low-moderate heritability in the current study underlines the potential for genetic improvement through selective breeding for IHNV resistance in the two commercial rainbow trout populations.

As described in the methods section, 150 fish (∼8%) were excluded from genotyping in the TLUN2018 population. This exclusion of ∼8% of the fish from the middle of the phenotypic distribution does represent small selective genotyping bias, but we believe that the effect of that small bias on over estimation of variance components must have been very small.

The heritability estimates based on genotype data were comparable across populations when using the same algorithm (0.15–0.16 with ssGBLUP and 0.23–0.25 with BMR-BayesB). Clearly the estimates with BMR-BayesB were higher than ssGBLUP. The difference is caused by the way that the two algorithms estimate the additive genetic variance. In ssGBLUP, the additive genetic variance is calculated using the estimated random individual animal effects. In comparison, with BMR-BayesB the estimate of additive genetic variance is the sum of the additive genetic variance due to random animal polygenic effects (i.e., breeding values of animals) and random marker effects (Fernando et al., 2016).

A larger difference in heritability values was observed using the pedigree and phenotype records (PBLUP; 0.08 and 0.23 for TLUM2019 and TLUN 2018, respectively). This difference may be caused by the reference population that is used by the different models. The reference for the genomic model is the same-generation population from which the genotypes and phenotypes were obtained (also called the training population for genomic selection models), whereas for PBLUP it is the initial generation in the pedigree records (Vallejo et al., 2021; Vallejo et al., 2022). The pedigree records for TLUN2018 trace back through eight generations to a single base population, providing a large dataset of pedigree records to estimate the genetic relatedness of the offspring with phenotypes from this population. In contrast, the TLUM2019 population was formed from a merger of two year-classes that were established from two different base-populations. Thus, the pedigree relies on only one generation of parents and offspring to calculate relatedness among the offspring. Therefore, we believe that the PBLUP heritability estimate for TLUN2018 is more reliable than that for TLUM 2019. The true heritability values for the two populations are likely within the genomic-based estimates of 0.15–0.25.

Two multiple regression GWAS methods were used to identify QTL associated with resistance to IHNV in this study. Such methods are powerful at the population level since they account for linkage disequilibrium between neighbouring loci, and utilize information from all available pedigree, genotype, and phenotype data (Garrick and Fernando, 2013; Misztal et al., 2014). We have previously shown that using two high quality and widely used algorithms for GWAS reveals more information about the genetic architecture of complex disease resistance traits in rainbow trout (Vallejo et al., 2017b; Vallejo et al., 2019; Vallejo et al., 2022). The detection of QTL with both algorithms increases the confidence in the detection of true QTL. In the current study, the QTL on chromosomes 1, 23 and 30 in TLUM2019 and on chromosome 29 in TLUN2018 were detected by both algorithms using the threshold of AGV >1.9. GWAS methods that are based on windows are sensitive to the SNP density and allele frequency (Li et al., 2021). In this study we used 1 Mb windows based on our experience with rainbow trout aquaculture populations where 1 Mb is approximately the average size of the LD extent of r 2 > 0.25 (Vallejo et al., 2018; Vallejo et al., 2020). We agree that very high or very low number of SNPs per window might affect the signal intensity and that LD is not evenly distributed across the genome.

Several of the genes located within the detected QTL regions may be implicated to be involved in immune response to infection based on their GO term and metabolic pathways predictions. However, it is important to caution that at this point we do not know if the true causative variant for each QTL is found in a protein coding gene. Therefore, our discussion here is limited to the gene functions that are predicted by the gene models that are based on orthology and sequence similarity to genes from model species. Furthermore, for practical reasons we restricted the list of annotated genes to those that were found within the 1-Mb QTL window boundaries, although it is possible that due to linkage disequilibrium with the causative genome variant, the causative genes or sequence variants may be located near rather than inside the SNP windows that were found to have the strongest association with the survival days phenotype. Here, we explore the genes from the five most supported QTL regions that may be involved in the innate or adaptive antiviral immunity of the rainbow trout based on the genome annotation and gene prediction models used in Refseq and Ensembl.

The shiftless antiviral inhibitor of ribosomal frameshifting (shfl) gene is located in the QTL region on chromosome 13. It is an interferon-stimulated gene that in human has been shown to function as a broad-spectrum antiviral factor with suppressive activity against various types of RNA and DNA viruses (Wang et al., 2019; Suzuki and Murakawa, 2022). Also located within the chromosome 13 QTL region is the notch receptor 3 gene (notch3), which is involved in the notch signalling pathway. Among its multiple functions, the notch signalling pathway was also implicated to be involved in innate lymphoid cell fate decision and immune response (Golub, 2021). IL15 receptor alpha chain (il15ra) is located within the QTL region on chromosome 21. The interleukin-15-mediated signalling pathway is involved in STAT activation of natural killing (NK) cells and has been implicated in the suppression of HIV latency (Macedo et al., 2022). Also in this QTL region, suppressor of cytokine signalling 2 (socs2) is a negative regulator of receptor signalling pathway via JAK-STAT and of type I interferon-mediated signalling pathway (Zhang et al., 2021; Kausar et al., 2022). Hence, by their known role in feed-back regulation of signalling pathways associated with viral suppression, both il15ra and socs2 are possible positional candidates for being the causative gene or genes responsible for the association of this chromosomal region with IHNV resistance in rainbow trout. The IKAROS family zinc finger one (ikzf1) is found in the QTL on chromosome 23. The structure and expression pattern of this gene confirm that it is a master switch of hematopoiesis and a key regulator of gene expression in lymphoid cells in rainbow trout (Hansen et al., 1997). In mammals, ikzf1 is required for B-cell development, and has been shown to regulate the development and function of dendritic cells and monocytes in human (Cytlak et al., 2018). Also located in this QTL region is the B cell linker protein, which is a molecular scaffold essential for the B cell receptor signalling pathway and required to promote B cell development (Pappu et al., 1999). The wingless-type MMTV integration site family member 2 (wnt2) is located in the QTL region on chromosome one and in cattle it was found through linkage disequilibrium analysis to be associated with resistance to a pathogenic bacterium and antibody response (Pauciullo et al., 2015). The QTL on chromosome Y (29) does not contain genes that can be implicated directly in antiviral immunity based on their functional annotation. However, it contains multiple genes encoding for the spliceosome sub-units 2 and 5. Alternative splicing has been recently implicated as a component of host-pathogen interactions in human and other mammals (Boudreault et al., 2019; Chauhan et al., 2019; Rotival et al., 2019; Yamaguchi et al., 2022; Lyu et al., 2023), including direct interactions between viral proteins and sub-unit 5 (U5) (Boudreault et al., 2022a; Boudreault et al., 2022b). Thus, it is plausible that variants of the spliceosome subunits may have an impact on the host-pathogen interactions after IHNV infection, and potentially may be the causative variant or variants for resistance to IHNV infection that was detected on chromosome Y QTL.

 Conclusion

The genetic architecture of resistance to IHNV in two commercial rainbow trout aquaculture strains was found to be oligogenic. None of the QTL detected had a large enough effect on the genetic variance for the trait that would merit marker assisted selection. Four of the five most supported QTL contain genes that may encode proteins that can be involved in antiviral function or that regulate components of the innate or adaptive immune response of the host. The fifth contain multiple sub-units of the spliceosome, which has also been recently implicated to be involved in host-pathogen interactions in human and other mammals. Further evaluation to refine those QTL regions with more genetic markers will bring us closer to the identification of the causative genome variants that are associated with IHNV resistance in the two commercially important rainbow trout aquaculture strains.

 Data availability statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.

 Ethics statement

The animal study was approved by the Institutional Animal Care and Use Committee of the University of Washington, Seattle, WA (Protocol number 4456-01). The study was conducted in accordance with the local legislation and institutional requirements.

 Author contributions

YP: Conceptualization, Data curation, Formal Analysis, Investigation, Methodology, Project administration, Resources, Supervision, Visualization, Writing–original draft, Writing–review and editing. RV: Data curation, Formal Analysis, Investigation, Methodology, Software, Validation, Visualization, Writing–original draft, Writing–review and editing. MP: Conceptualization, Investigation, Methodology, Resources, Supervision, Visualization, Writing–review and editing. GG: Data curation, Methodology, Software, Writing–review and editing. KS: Data curation, Methodology, Writing–review and editing. RL: Data curation, Methodology, Writing–review and editing. CS: Investigation, Validation, Writing–review and editing. BF: Software, Writing–review and editing. HC: Software, Writing–review and editing. KM: Conceptualization, Data curation, Funding acquisition, Methodology, Resources, Supervision, Writing–review and editing. KN: Conceptualization, Funding acquisition, Project administration, Resources, Supervision, Writing–review and editing.

 Funding

The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This project was supported by the Western Regional Aquaculture Center [grant numbers 2016- 38500-25755 and 2018-38500-28889] from the USDA National Institute of Food and Agriculture and the USDA Agricultural Research Service in house project number 8082-31000- 013. Mention of trade names or commercial products in this publication is solely for the purpose of providing specific information and does not imply recommendation or endorsement by the U.S. government. USDA is an equal opportunity provider and employer.

 Conflict of interest

KM who is one of the co-authors of the manuscript is an employee of Troutlodge, Inc., a daughter company of Hendrix Genetics. Rainbow trout fish from the germplasm developed by Troutlodge were used in this study and their resistance to IHNV disease is being characterized in the study including the genetic architecture of this commercial germplasm. However, there is no endorsement in the manuscript of a commercial product from Troutlodge or Hendrix Genetics. Mention of trade names or commercial products in this publication is solely for the purpose of providing specific information and does not imply recommendation or endorsement by the U.S. government.

The author(s) declared that they were an editorial board member of Frontiers, at the time of submission. This had no impact on the peer review process and the final decision.

 Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

 Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2024.1394656/full#supplementary-material

Perl script to assign best trio matching of parental pair genotypes to offspring genotypes.

Map of markers’ genome positions for TLUM 2019.

Map of markers’ genome positions for TLUN 2018.

Genes located within the most supportive 1-Mb QTL windows.

GO terms assigned to the genes in the QTL windows.

Overlapping QTL regions with the Clear Springs Food study (Vallejo et al., 2019).

 References Adel M. Amiri A. B. Dadar M. Breyta R. Kurath G. Laktarashi B. (2016). Phylogenetic relationships of Iranian infectious hematopoietic necrosis virus of rainbow trout (Oncorhynchus mykiss) based on the glycoprotein gene. Arch. Virol. 161, 657–663. 10.1007/s00705-015-2684-8 Aguilar I. Misztal I. Johnson D. L. Legarra A. Tsuruta S. Lawlor T. J. (2010). Hot topic: a unified approach to utilize phenotypic, full pedigree, and genomic information for genetic evaluation of Holstein final score. J. Dairy Sci. 93, 743–752. 10.3168/jds.2009-2730 Anderson E. D. Mourich D. V. Fahrenkrug S. C. La Patra S. E. Shepard J. Leong J. C. (1996). Genetic immunization of rainbow trout (Oncorhynchus mykiss) against infectious hematopoietic necrosis virus. Mol. Mar. Biol. Biotechnol. 5, 114–122. Barabaschi D. Tondelli A. Desiderio F. Volante A. Vaccino P. Vale G. (2016). Next generation breeding. Plant Sci. 242, 3–13. 10.1016/j.plantsci.2015.07.010 Batts W. N. Winton J. R. (1989). Enhanced detection of infectious hematopoietic necrosis virus and other fish viruses by pretreatment of cell monolayers with polyethylene glycol. J. Aquatic Animal Health 1, 284–290. 10.1577/1548-8667(1989)001<0284:edoihn>2.3.co;2 Bledsoe J. W. Ma J. Cain K. Bruce T. J. Rawles A. Abernathy J. (2022). Multi-tissue RNAseq reveals genetic and temporal differences in acute response to viral (IHNV) infection among three selected lines of rainbow trout with varying resistance. Fish Shellfish Immunol. 124, 343–361. 10.1016/j.fsi.2022.03.034 Bootland L. M. Leong J. C. (2009). “Infectious haematopoietic necrosis virus,” in Fish diseases and disorders. Editors Woo P. T. K. Bruno D. W. (Cambridge: CAB International), 66–109. Boudreault S. Durand M. Martineau C.-A. Perreault J.-P. Lemay G. Bisaillon M. (2022b). Reovirus μ2 protein modulates host cell alternative splicing by reducing protein levels of U5 snRNP core components. Nucleic Acids Res. 50, 5263–5281. 10.1093/nar/gkac272 Boudreault S. Lemay G. Bisaillon M. (2022a). U5 snRNP core proteins are key components of the defense response against viral infection through their roles in programmed cell Death and interferon induction. Viruses 14, 2710. 10.3390/v14122710 Boudreault S. Roy P. Lemay G. Bisaillon M. (2019). Viral modulation of cellular RNA alternative splicing: a new key player in virus–host interactions? WIREs RNA 10, e1543. 10.1002/wrna.1543 Breyta R. Jones A. Stewart B. Brunson R. Thomas J. Kerwin J. (2013). Emergence of MD type infectious hematopoietic necrosis virus in Washington State coastal steelhead trout. Dis. Aquat. Organ 104, 179–195. 10.3354/dao02596 Brieuc M. S. Purcell M. K. Palmer A. D. Naish K. A. (2015). Genetic variation underlying resistance to infectious hematopoietic necrosis virus in a steelhead trout (Oncorhynchus mykiss) population. Dis. Aquat. Organ. 117, 77–83. 10.3354/dao02933 Campbell N. R. LaPatra S. E. Overturf K. Towner R. Narum S. R. (2014). Association mapping of disease resistance traits in rainbow trout using restriction site associated DNA sequencing. G3 Genes.|Genomes|Genetics 4, 2473–2481. 10.1534/g3.114.014621 Chauhan K. Kalam H. Dutt R. Kumar D. (2019). RNA splicing: a new paradigm in host–pathogen interactions. J. Mol. Biol. 431, 1565–1575. 10.1016/j.jmb.2019.03.001 Cheng H. Fernando R. Garrick D. (2018) JWAS: julia implementation of whole-genome analysis software. Auckland, New Zealand: World Congress on Genetics Applied to Livestock Production. Cytlak U. Resteu A. Bogaert D. Kuehn H. S. Altmann T. Gennery A. (2018). Ikaros family zinc finger 1 regulates dendritic cell development and function in humans. Nat. Commun. 9, 1239. 10.1038/s41467-018-02977-8 Enzmann P. J. Kurath G. Fichtner D. Bergmann S. M. (2005). Infectious hematopoietic necrosis virus: monophyletic origin of European isolates from North American genogroup M. Dis. Aquatic Org. 66, 187–195. 10.3354/dao066187 Fernando R. Toosi A. Wolc A. Garrick D. Dekkers J. (2017). Application of whole-genome prediction methods for genome-wide association studies: a bayesian approach. J. Agric. Biol. Environ. Statistics 22, 172–193. 10.1007/s13253-017-0277-6 Fernando R. L. Cheng H. Golden B. L. Garrick D. J. (2016). Computational strategies for alternative single-step Bayesian regression models with large numbers of genotyped and non-genotyped animals. Genet. Sel. Evol. 48, 96. 10.1186/s12711-016-0273-2 Fernando R. L. Dekkers J. C. M. Garrick D. J. (2014). A class of Bayesian methods to combine large numbers of genotyped and non-genotyped animals for whole-genome analyses. Genet. Sel. Evol. 46, 50. 10.1186/1297-9686-46-50 Gao G. Magadan S. Waldbieser G. C. Youngblood R. C. Wheeler P. A. Scheffler B. E. (2021). A long reads-based de-novo assembly of the genome of the Arlee homozygous line reveals chromosomal rearrangements in rainbow trout. G3 Genes.|Genomes|Genetics 11, jkab052. 10.1093/g3journal/jkab052 Garrick D. J. Fernando R. L. (2013). “Implementing a QTL detection study (GWAS) using genomic prediction methodology,” in Genome-wide association studies and genomic prediction. Editors Gondro C. Werf J. Hayes B. (New York: Humana Press). Goddard M. E. Hayes B. J. Meuwissen T. H. E. (2011). Using the genomic relationship matrix to predict the accuracy of genomic selection. J. Animal Breed. Genet. 128, 409–421. 10.1111/j.1439-0388.2011.00964.x Goldstein B. A. Hubbard A. E. Cutler A. Barcellos L. F. (2010). An application of Random Forests to a genome-wide association dataset: methodological considerations & new findings. BMC Genet. 11, 49. 10.1186/1471-2156-11-49 Golub R. (2021). The Notch signaling pathway involvement in innate lymphoid cell biology. Biomed. J. 44, 133–143. 10.1016/j.bj.2020.12.004 Hansen J. D. Strassburger P. Du Pasquier L. (1997). Conservation of a master hematopoietic switch gene during vertebrate evolution: isolation and characterization of Ikaros from teleost and amphibian species. Eur. J. Immunol. 27, 3049–3058. 10.1002/eji.1830271143 Houston R. D. Haley C. S. Hamilton A. Guy D. R. Mota-Velasco J. C. Gheyas A. A. (2010). The susceptibility of Atlantic salmon fry to freshwater infectious pancreatic necrosis is largely explained by a major QTL. Heredity 105, 318–327. 10.1038/hdy.2009.171 Houston R. D. Haley C. S. Hamilton A. Guy D. R. Tinch A. E. Taggard J. B. (2008). Major quantitative trait loci affect resistance to infectious pancreatic necrosis in Atlantic salmon (Salmo salar). Genetics 178, 1109–1115. 10.1534/genetics.107.082974 Irano N. de Camargo G. M. Costa R. B. Terakado A. P. Magalhaes A. F. Silva R. M. (2016). Genome-wide association study for indicator traits of sexual precocity in nellore cattle. PLoS One 11, e0159502. 10.1371/journal.pone.0159502 Kausar S. Gul I. Liu R. Ke X.-X. Dong Z. Abbas M. N. (2022). Antheraea pernyi suppressor of cytokine signaling 2 negatively modulates the JAK/STAT pathway to attenuate microbial infection. Int. J. Mol. Sci. 23, 10389. 10.3390/ijms231810389 Kurath G. Garver K. A. Troyer R. M. Emmenegger E. J. Einer-Jensen K. Anderson E. D. (2003). Phylogeography of infectious haematopoietic necrosis virus in North America. J. Gen. Virol. 84, 803–814. 10.1099/vir.0.18771-0 LaPatra S. Kao S. Erhardt E. B. Salinas I. (2015). Evaluation of dual nasal delivery of infectious hematopoietic necrosis virus and enteric red mouth vaccines in rainbow trout (Oncorhynchus mykiss). Vaccine 33, 771–776. 10.1016/j.vaccine.2014.12.055 LaPatra S. E. (1998). Factors affecting pathogenicity of infectious hematopoietic necrosis virus (IHNV) for salmonid fish. J. Aquatic Animal Health 10, 121–131. 10.1577/1548-8667(1998)010<0121:fapoih>2.0.co;2 LaPatra S. E. Groberg W. J. Rohovec J. S. Fryer J. L. (1990). Size-related susceptibility of salmonids to two strains of infectious hematopoietic necrosis virus. Trans. Am. Fish. Soc. 119, 25–30. 10.1577/1548-8659(1990)119<0025:ssostt>2.3.co;2 Larragoite E. T. Tacchi L. LaPatra S. E. Salinas I. (2016). An attenuated virus vaccine appears safe to the central nervous system of rainbow trout (Oncorhynchus mykiss) after intranasal delivery. Fish. Shellfish Immunol. 49, 351–354. 10.1016/j.fsi.2016.01.006 Li J. Wang Z. Fernando R. Cheng H. (2021). Tests of association based on genomic windows can lead to spurious associations when using genotype panels with heterogeneous SNP densities. Genet. Sel. Evol. 53, 45. 10.1186/s12711-021-00638-x Liu S. Palti Y. Gao G. Rexroad Iii C. E. (2016). Development and validation of a SNP panel for parentage assignment in rainbow trout. Aquaculture 452, 178–182. 10.1016/j.aquaculture.2015.11.001 Liu S. Palti Y. Martin K. E. Parsons J. E. Rexroad Iii C. E. (2017). Assessment of genetic differentiation and genetic assignment of commercial rainbow trout strains using a SNP panel. Aquaculture 468 (Part 1), 120–125. 10.1016/j.aquaculture.2016.10.004 Lorenzen N. Lapatra S. E. (1999). Immunity to rhabdoviruses in rainbow trout: the antibody response. Fish Shellfish Immunol. 9, 345–360. 10.1006/fsim.1999.0194 Lorenzen N. LaPatra S. E. (2005). Vacunas de ADN para peces de vivero: -EN- DNA vaccines for aquacultured fish -FR- Vaccins à ADN destinés aux poissons d’élevage -ES-. Revue Sci. Technique-Office Int. Des Epizooties 24, 201–213. 10.20506/rst.24.1.1565 Lyu M. Lai H. Wang Y. Zhou Y. Chen Y. Wu D. (2023). Roles of alternative splicing in infectious diseases: from hosts, pathogens to their interactions. Chin. Med. J. 136, 767–779. 10.1097/CM9.0000000000002621 Macedo A. B. Levinger C. Nguyen B. N. Richard J. Gupta M. Cruz C. R. Y. (2022). The HIV latency reversal agent HODHBt enhances NK cell effector and memory-like functions by increasing interleukin-15-mediated STAT activation. J. Virology 96, 00372222–e100322. 10.1128/jvi.00372-22 Melo T. P. Takada L. Baldi F. Oliveira H. N. Dias M. M. Neves H. H. (2016). Assessing the value of phenotypic information from non-genotyped animals for QTL mapping of complex traits in real and simulated populations. BMC Genet. 17, 89. 10.1186/s12863-016-0394-1 Meuwissen T. H. E. Hayes B. J. Goddard M. E. (2001). Prediction of total genetic value using genome-wide dense marker maps. Genetics 157, 1819–1829. 10.1093/genetics/157.4.1819 Misztal I. Tsuruta S. Lourenco D. Aguilar I. Legarra A. Vitezica Z. (2015) Manual for BLUPF90 family of programs. Athens, GA, USA: University of Georgia. Misztal I. Wang H. Aguilar I. Legarra A. Tsuruta S. Lourenco D. A. L. (2014) GWAS using ssGBLUP. Proceedings, 10th world congress of genetics applied to Livestock production. Moen T. Sonesson A. K. Hayes B. Lien S. Munck H. Meuwissen T. H. E. (2007). Mapping of a quantitative trait locus for resistance against infectious salmon anaemia in Atlantic salmon (Salmo Salar): comparing survival analysis with analysis on affected/resistant data. BMC Genet. 8, 53. 10.1186/1471-2156-8-53 Moen T. Torgersen J. Santi N. Davidson W. S. Baranski M. Odegard J. (2015). Epithelial cadherin determines resistance to infectious pancreatic necrosis virus in atlantic salmon. Genetics 200, 1313–1326. 10.1534/genetics.115.175406 Pappu R. Cheng A. M. Li B. Gong Q. Chiu C. Griffin N. (1999). Requirement for B cell linker protein (BLNK) in B cell development. Science. 286, 1949–1954. 10.1126/science.286.5446.1949 Pauciullo A. Küpper J. Brandt H. Donat K. Iannuzzi L. Erhardt G. (2015). Wingless-type MMTV integration site family member 2 (WNT2) gene is associated with resistance to MAP in faecal culture and antibody response in Holstein cattle. Anim. Genet. 46, 122–132. 10.1111/age.12261 Pearse D. E. Barson N. J. Nome T. Gao G. Campbell M. A. Abadía-Cardoso A. (2019). Sex-dependent dominance maintains migration supergene in rainbow trout. Nat. Ecol. Evol. 3, 1731–1742. 10.1038/s41559-019-1044-6 Plant K. P. LaPatra S. E. (2011). Advances in fish vaccine delivery. Dev. Comp. Immunol. 35, 1256–1262. 10.1016/j.dci.2011.03.007 Plummer M. Best N. Cowles K. Vines K. (2006). CODA: convergence diagnosis and output analysis for MCMC. R. News 6. Purcell M. K. LaPatra S. E. Woodson J. C. Kurath G. Winton J. R. (2010). Early viral replication and induced or constitutive immunity in rainbow trout families with differential resistance to Infectious hematopoietic necrosis virus (IHNV). Fish Shellfish Immunol. 28, 98–105. 10.1016/j.fsi.2009.10.005 Quillet E. Dorson M. Le Guillou S. Benmansour A. Boudinot P. (2007). Wide range of susceptibility to rhabdoviruses in homozygous clones of rainbow trout. Fish Shellfish Immunol. 22, 510–519. 10.1016/j.fsi.2006.07.002 Rotival M. Quach H. Quintana-Murci L. (2019). Defining the genetic and evolutionary architecture of alternative splicing in response to infection. Nat. Commun. 10, 1671. 10.1038/s41467-019-09689-7 Saksida S. M. (2006). Infectious haematopoietic necrosis epidemic (2001 to 2003) in farmed Atlantic salmon Salmo salar in British Columbia. Dis. Aquat. Organ 72, 213–223. 10.3354/dao072213 Silva R. M. O. Evenhuis J. P. Vallejo R. L. Gao G. Martin K. E. Leeds T. D. (2019). Whole-genome mapping of quantitative trait loci and accuracy of genomic predictions for resistance to columnaris disease in two rainbow trout breeding populations. Genet. Sel. Evol. 51, 42. 10.1186/s12711-019-0484-4 Song H. Dong T. Yan X. Wang W. Tian Z. Sun A. (2023). Genomic selection and its research progress in aquaculture breeding. Rev. Aquac. 15, 274–291. 10.1111/raq.12716 Suzuki Y. Murakawa T. (2022). Restriction of flaviviruses by an interferon-stimulated gene SHFL/C19orf66. Int. J. Mol. Sci. 23, 12619. 10.3390/ijms232012619 Troyer R. M. Kurath G. (2003). Molecular epidemiology of infectious hematopoietic necrosis virus reveals complex virus traffic and evolution within southern Idaho aquaculture. Dis. Aquatic Org. 55, 175–185. 10.3354/dao055175 Vallejo R. L. Cheng H. Fragomeni B. O. Gao G. Silva R. M. O. Martin K. E. (2021). The accuracy of genomic predictions for bacterial cold water disease resistance remains higher than the pedigree-based model one generation after model training in a commercial rainbow trout breeding population. Aquaculture 545, 737164. 10.1016/j.aquaculture.2021.737164 Vallejo R. L. Cheng H. Fragomeni B. O. Shewbridge K. L. Gao G. MacMillan J. R. (2019). Genome-wide association analysis and accuracy of genome-enabled breeding value predictions for resistance to infectious hematopoietic necrosis virus in a commercial rainbow trout breeding population. Genet. Sel. Evol. 51, 47. 10.1186/s12711-019-0489-z Vallejo R. L. Evenhuis J. P. Cheng H. Fragomeni B. O. Gao G. Liu S. (2022). Genome-wide mapping of quantitative trait loci that can be used in marker-assisted selection for resistance to bacterial cold water disease in two commercial rainbow trout breeding populations. Aquaculture 738574. 10.1016/j.aquaculture.2022.738574 Vallejo R. L. Fragomeni B. O. Cheng H. Gao G. Long R. L. Shewbridge K. L. (2020). Assessing accuracy of genomic predictions for resistance to infectious hematopoietic necrosis virus with progeny testing of selection candidates in a commercial rainbow trout breeding population. Front. Veterinary Sci. 7, 590048. 10.3389/fvets.2020.590048 Vallejo R. L. Leeds T. D. Fragomeni B. O. Gao G. Hernandez A. G. Misztal I. (2016). Evaluation of genome-enabled selection for bacterial cold water disease resistance using progeny performance data in rainbow trout: insights on genotyping methods and genomic prediction models. Front. Genet. 7, 96. 10.3389/fgene.2016.00096 Vallejo R. L. Leeds T. D. Gao G. Parsons J. E. Martin K. E. Evenhuis J. P. (2017a). Genomic selection models double the accuracy of predicted breeding values for bacterial cold water disease resistance compared to a traditional pedigree-based model in rainbow trout aquaculture. Genet. Sel. Evol. 49, 17. 10.1186/s12711-017-0293-6 Vallejo R. L. Liu S. Gao G. Fragomeni B. O. Hernandez A. G. Leeds T. D. (2017b). Similar genetic architecture with shared and unique quantitative trait loci for bacterial cold water disease resistance in two rainbow trout breeding populations. Front. Genet. 8, 156. 10.3389/fgene.2017.00156 Vallejo R. L. Silva R. M. O. Evenhuis J. P. Gao G. Liu S. Parsons J. E. (2018). Accurate genomic predictions for BCWD resistance in rainbow trout are achieved using low-density SNP panels: evidence that long-range LD is a major contributing factor. J. Animal Breed. Genet. 135, 263–274. 10.1111/jbg.12335 Visscher P. M. Brown M. A. McCarthy M. I. Yang J. (2012). Five years of GWAS discovery. Am. J. Hum. Genet. 90, 7–24. 10.1016/j.ajhg.2011.11.029 Wang C. Zhao L. L. Li Y. J. Tang L. J. Qiao X. Y. Jiang Y. P. (2016). Analysis of the genome sequence of infectious hematopoietic necrosis virus HLJ-09 in China. Virus Genes. 52, 29–37. 10.1007/s11262-015-1263-0 Wang H. Misztal I. Aguilar I. Legarra A. Muir W. M. (2012). Genome-wide association mapping including phenotypes from relatives without genotypes. Genet. Res. 94, 73–83. 10.1017/S0016672312000274 Wang X. Xuan Y. Han Y. Ding X. Ye K. Yang F. (2019). Regulation of HIV-1 gag-pol expression by shiftless, an inhibitor of programmed -1 ribosomal frameshifting. Cell. 176, 625–635. 10.1016/j.cell.2018.12.030 Winton J. Batts W. DeKinkelin P. LeBerre M. Bremont M. Fijan N. (2010). Current lineages of the epithelioma papulosum cyprini (EPC) cell line are contaminated with fathead minnow, Pimephales promelas, cells. J. Fish Dis. 33, 701–704. 10.1111/j.1365-2761.2010.01165.x Yamaguchi K. Ishigaki K. Suzuki A. Tsuchida Y. Tsuchiya H. Sumitomo S. (2022). Splicing QTL analysis focusing on coding sequences reveals mechanisms for disease susceptibility loci. Nat. Commun. 13, 4659. 10.1038/s41467-022-32358-1 Yáñez J. M. Barría A. López M. E. Moen T. Garcia B. F. Yoshida G. M. (2023). Genome-wide association and genomic selection in aquaculture. Rev. Aquac. 15, 645–675. 10.1111/raq.12750 Zhang X. Zhuang Y. Qin T. Chang M. Ji X. Wang N. (2021). Suppressor of cytokine signalling-2 controls hepatic gluconeogenesis and hyperglycemia by modulating JAK2/STAT5 signalling pathway. Metabolism 122, 154823. 10.1016/j.metabol.2021.154823