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

Genetic diversity patterns in farmed rainbow trout (Oncorhynchus mykiss) populations using genome-wide SNP and haplotype data

Rainbow trout is one of the most popular aquaculture species worldwide, with a long history of domestication. However, limited information exists about the genetic diversity of farmed rainbow trout populations globally, with most available reports relying on low-throughput genotyping technologies. Notably, no information exists about the genetic diversity status of farmed rainbow trout in Sweden. Double-digest restriction-site-associated DNA sequencing was performed on more than 500 broodfish from two leading producers in Sweden and from the country's national breeding program. Following the detection of single nucleotide polymorphisms (SNPs), genetic diversity was studied by using either individual SNPs (n = 8680; one SNP retained per 300 bp sequence reads) or through SNP haplotypes (n = 20 558; all SNPs retained in 300 bp sequence reads). Similar amounts of genetic diversity were found amongst the three populations when individual SNPs were used. Furthermore, principal component analysis and discriminant analysis of principal components suggested two genetic clusters with the two industry populations grouped together. Genetic differentiation based on the FST fixation index was ~0.01 between the industry populations and ~0.05 when those were compared with the breeding program. Preliminary estimates of effective population size (Ne ) and inbreeding (based on runs of homozygosity; FROH ) were similar amongst the three populations (Ne  ≈ 50-80; median FROH  ≈ 0.11). Finally, the haplotype-based analysis suggested that animals from the breeding program had higher shared coancestry levels than those from the other two populations. Overall, our study provides novel insights into the genetic diversity and structure of Sweden's three main farmed rainbow trout populations, which could guide their future management.

Plasma proteomic profiling of bacterial cold water disease-resistant and -susceptible rainbow trout lines and biomarker discovery

Genetic variation for disease resistance is present in salmonid fish; however, the molecular basis is poorly understood, and biomarkers of disease susceptibility/resistance are unavailable. Previously, we selected a line of rainbow trout for high survival following standardized challenge with Flavobacterium psychrophilum (Fp), the causative agent of bacterial cold water disease. The resistant line (ARS-Fp-R) exhibits over 60 percentage points higher survival compared to a reference susceptible line (ARS-Fp-S). To gain insight into the differential host response between genetic lines, we compared the plasma proteomes from day 6 after intramuscular challenge. Pooled plasma from unhandled, PBS-injected, and Fp-injected groups were simultaneously analyzed using a TMT 6-plex label, and the relative abundance of 513 proteins was determined. Data are available via ProteomeXchange, with identifier PXD041308, and the relative protein abundance values were compared to mRNA measured from a prior, whole-body RNA-seq dataset. Our results identified a subset of differentially abundant intracellular proteins was identified, including troponin and myosin, which were not transcriptionally regulated, suggesting that these proteins were released into plasma following pathogen-induced tissue damage. A separate subset of high-abundance, secreted proteins were transcriptionally regulated in infected fish. The highest differentially expressed protein was a C1q family member (designated complement C1q-like protein 3; C1q-LP3) that was upregulated over 20-fold in the infected susceptible line while only modestly upregulated, 1.8-fold, in the infected resistant line. Validation of biomarkers was performed using immunoassays and C1q-LP3, skeletal muscle troponin C, cathelcidin 2, haptoglobin, leptin, and growth and differentiation factor 15 exhibited elevated concentration in susceptible line plasma. Complement factor H-like 1 exhibited higher abundance in the resistant line compared to the susceptible line in both control and challenged fish and thus was a baseline differentiator between lines. C1q-LP3 and STNC were elevated in Atlantic salmon plasma following experimental challenge with Fp. In summary, these findings further the understanding of the differential host response to Fp and identifies salmonid biomarkers that may have use for genetic line evaluation and on-farm health monitoring.

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.

Age-Dependent Growth-Related QTL Variations in Pacific Abalone, Haliotis discus hannai

Pacific abalone is a high-value, commercially important marine invertebrate. It shows low growth as well as individual and yearly growth variation in aquaculture. Marker-assisted selection breeding could potentially resolve the problem of low and variable growth and increase genetic gain. Expression of quantitative trait loci (QTLs) for growth-related traits, viz., body weight, shell length, and shell width were analyzed at the first, second, and third year of age using an F1 cross population. A total of 37 chromosome-wide QTLs were identified in linkage groups 01, 02, 03, 04, 06, 07, 08, 10, 11, 12, and 13 at different ages. None of the QTLs detected at any one age were expressed in all three age groups. This result suggests that growth-related traits at different ages are influenced by different QTLs in each year. However, multiple-trait QTLs (where one QTL affects all three traits) were detected each year that are also age-specific. Eleven multiple-trait QTLs were detected at different ages: two QTLs in the first year; two QTLs in the second year; and seven QTLs in the third year. As abalone hatcheries use three-year-old abalone for breeding, QTL-linked markers that were detected at the third year of age could potentially be used in marker-assisted selection breeding programs.

Potential of low-density genotype imputation for cost-efficient genomic selection for resistance to Flavobacterium columnare in rainbow trout (Oncorhynchus mykiss)

BACKGROUND

Flavobacterium columnare is the pathogen agent of columnaris disease, a major emerging disease that affects rainbow trout aquaculture. Selective breeding using genomic selection has potential to achieve cumulative improvement of the host resistance. However, genomic selection is expensive partly because of the cost of genotyping large numbers of animals using high-density single nucleotide polymorphism (SNP) arrays. The objective of this study was to assess the efficiency of genomic selection for resistance to F. columnare using in silico low-density (LD) panels combined with imputation. After a natural outbreak of columnaris disease, 2874 challenged fish and 469 fish from the parental generation (n = 81 parents) were genotyped with 27,907 SNPs. The efficiency of genomic prediction using LD panels was assessed for 10 panels of different densities, which were created in silico using two sampling methods, random and equally spaced. All LD panels were also imputed to the full 28K HD panel using the parental generation as the reference population, and genomic predictions were re-evaluated. The potential of prioritizing SNPs that are associated with resistance to F. columnare was also tested for the six lower-density panels.

RESULTS

The accuracies of both imputation and genomic predictions were similar with random and equally-spaced sampling of SNPs. Using LD panels of at least 3000 SNPs or lower-density panels (as low as 300 SNPs) combined with imputation resulted in accuracies that were comparable to those of the 28K HD panel and were 11% higher than the pedigree-based predictions.

CONCLUSIONS

Compared to using the commercial HD panel, LD panels combined with imputation may provide a more affordable approach to genomic prediction of breeding values, which supports a more widespread adoption of genomic selection in aquaculture breeding programmes.

Validation of a QTL associated with resistance to Vibrio anguillarum in rainbow trout (Oncorhynchus mykiss)

Vibriosis is a bacterial disease in fish caused by the Gram negative bacterium Vibrio anguillarum with severe impact on rainbow trout (Oncorhynchus mykiss) farming. Sustainable control methods should be developed and we here show that marker assisted selective breeding of fish naturally resistant to the disease is feasible. We have validated the use of a single nucleotide polymorphism (SNP) marker SNP AX-89,945,921 (QTL on chromosome 21). The QTL was previously found associated with resistance to vibriosis and described following a genome wide association analysis (GWAS) of trout exposed to the bacterium. For this validation spawners were genotyped by use of the 57 K Axiom®Trout Microarray (Affymetrix) and homozygous male fish carrying the allele with the SNP AX-89,945,921 were then selected and used to fertilize eggs from outbred female trout resulting in fish all carrying the SNP (QTL-fish). Control fish (non-QTL fish) were produced by fertilizing the same batch of eggs by use of male parents negative for the SNP. The fish were exposed in freshwater to V. anguillarum (water bath infection) at 19 C°. A total of 900 fish were challenged in a common garden set-up in triplicate. A bacterial solution of V. anguillarum (serotype O1) was added to each of three freshwater fish tanks, each with 150 QTL and 150 non-QTL fish. Fish were tagged by tail fin cut (upper/lower) to discern the two groups, whereafter fish were monitored around the clock to detect disease signs and remove moribund fish. Clinical vibriosis developed within two days in non-QTL-fish (overall morbidity of 70%). QTL fish developed clinical signs later and the morbidity was significantly lower and did not reach 50%. Rainbow trout farming may benefit from using the QTL associated with higher resistance towards vibriosis. The effect may be optimized in the future by use of both male and female parents homozygous for the marker allele.

Conserved QTL and chromosomal inversion affect resistance to columnaris disease in 2 rainbow trout (Oncorhyncus mykiss) populations

We present a comparative genetic analysis of the quantitative trait loci underlying resistance to warm water columnaris disease in 2 farmed rainbow trout (Oncorhynchus mykiss) populations. We provide evidence for the conservation of a major quantitative trait loci on Omy03, and the putative role played by a chromosomal rearrangement on Omy05. A total of 3,962 individuals from the 2 populations experienced a natural Flavobacterium columnare outbreak. Data for 25,823 genome-wide SNPs were generated for both cases (fatalities) and controls (survivors). F_ST and pairwise additive genetic relationships suggest that, despite being currently kept as separate broodstocks, the 2 populations are closely related. Association analyses identified a major quantitative trait loci on chromosome Omy03 and a second smaller quantitative trait loci on Omy05. Quantitative trait loci on Omy03 consistently explained 3–11% of genetic variation in both populations, whereas quantitative trait loci on Omy05 showed different degree of association across populations and sexes. The quantitative trait loci on Omy05 was found within a naturally occurring, 54.84 cM long inversion which is easy to tag due to a strong linkage disequilibrium between the 375 tagging SNPs. The ancestral haplotype on Omy05 was associated with decreased mortality. Genetic correlation between mortality in the 2 populations was estimated at 0.64, implying that the genetic basis of resistance is partly similar in the 2 populations. Our quantitative trait loci validation identifies markers that can be potentially used to complement breeding value evaluations to increase resistance against columnaris disease, and help to mitigate effects of climate change on aquaculture.

Identification of Haplotypes Associated With Resistance to Bacterial Cold Water Disease in Rainbow Trout Using Whole-Genome Resequencing

Bacterial cold water disease (BCWD) is an important disease in rainbow trout aquaculture. Previously, we have identified and validated two major QTL (quantitative trait loci) for BCWD resistance, located on chromosomes Omy08 and Omy25, in the odd-year Troutlodge May spawning population. We also demonstrated that marker-assisted selection (MAS) for BCWD resistance using the favorable haplotypes associated with the two major QTL is feasible. However, each favorable haplotype spans a large genomic region of 1.3–1.6 Mb. Recombination events within the haplotype regions will result in new haplotypes associated with BCWD resistance, which will reduce the accuracy of MAS for BCWD resistance over time. The objectives of this study were 1) to identify additional SNPs (single nucleotide polymorphisms) associated with BCWD resistance using whole-genome sequencing (WGS); 2) to validate the SNPs associated with BCWD resistance using family-based association mapping; 3) to refine the haplotypes associated with BCWD resistance; and 4) to evaluate MAS for BCWD resistance using the refined QTL haplotypes. Four consecutive generations of the Troutlodge May spawning population were evaluated for BCWD resistance. Parents and offspring were sequenced as individuals and in pools based on their BCWD phenotypes. Over 12 million SNPs were identified by mapping the sequences from the individuals and pools to the reference genome. SNPs with significantly different allele frequencies between the two BCWD phenotype groups were selected to develop SNP assays for family-based association mapping in three consecutive generations of the Troutlodge May spawning population. Among the 78 SNPs derived from WGS, 77 SNPs were associated with BCWD resistance in at least one of the three consecutive generations. The additional SNPs associated with BCWD resistance allowed us to reduce the physical sizes of haplotypes associated with BCWD resistance to less than 0.5 Mb. We also demonstrated that the refined QTL haplotypes can be used for MAS in the Troutlodge May spawning population. Therefore, the SNPs and haplotypes reported in this study provide additional resources for improvement of BCWD resistance in rainbow trout.

Quantification and comparison of gene expression associated with iron regulation and metabolism in a virulent and attenuated strain of Flavobacterium psychrophilum

Iron is essential for growth and virulence in most pathogenic bacterial strains. In some cases, the hosts for these pathogenic bacteria develop specialized strategies to sequester iron and limit infectivity. This in turn may result in the invading pathogens utilizing high-affinity iron transport mechanisms, such as the use of iron-chelating siderophores, to extend beyond the host defences. Flavobacterium psychrophilum, the causative agent of bacterial coldwater disease (BCWD) in salmonids, relies on iron metabolism for infectivity, and the genome of the model CSF-259-93 strain has recently been made available. Further, this strain serves as a parent strain for a live-attenuated vaccine strain, B.17, which has been shown to provide rainbow trout with protection against BCWD. To elucidate specific gene expression responses to iron metabolism and compare strain differences, both F. psychrophilum strains were grown under iron-limiting conditions and 26 genes related to iron metabolism were mapped for 96 hr in culture via qPCR analyses. Results indicate increased production of the ferrous iron transport protein B (FITB; p =.008), and ferric receptor CfrA (FR 1; p =.012) in the wild-type CSF-259-93 strain at 72 hr and 96 hr post-exposure to iron-limiting media. In the B.17 vaccine strain, siderophore synthase (SS) expression was found to be downregulated at 72 hr, in comparison with 0h (p =.018). When strains were compared, FITB (p =.021), FR1 (p =.009) and SS (p =.016) were also elevated in B.17 at 0 hr and TonB outer protein membrane receptor 1 (TBomr1; p =.005) had a lower expression at 96 hr. Overall, this study demonstrated strain-related gene expression changes in only a fraction of the iron metabolism genes tested; however, results provide insight on potential virulence mechanisms and clarification on iron-related gene expression for F. psychrophilum.

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.

A Major QTL for Resistance to Vibrio anguillarum in Rainbow Trout

Genetic selection of disease resistant fish is a major strategy to improve health, welfare and sustainability in aquaculture. Mapping of single nucleotide polymorphisms (SNP) in the fish genome may be a fruitful tool to define relevant quantitative trait loci (QTL) and we here show its use for characterization of Vibrio anguillarum resistant rainbow trout (Oncorhynchus mykiss). Fingerlings were exposed to the pathogen V. anguillarum serotype O1 in a solution of 1.5 × 10⁷ cfu/ml and observed for 14 days. Disease signs appeared 3 days post exposure (dpe) whereafter mortality progressed exponentially until 6 dpe reaching a total mortality of 55% within 11 days. DNA was sampled from all fish – including survivors – and analyzed on a 57 k Affymetrix SNP platform whereby it was shown that disease resistance was associated with a major QTL on chromosome 21 (Omy 21). Gene expression analyses showed that diseased fish activated genes associated with innate and adaptive immune responses. The possible genes associated with resistance are discussed.

The Type IX Secretion System Is Required for Virulence of the Fish Pathogen Flavobacterium psychrophilum

Flavobacterium psychrophilum causes bacterial cold-water disease in wild and aquaculture-reared fish and is a major problem for salmonid aquaculture. The mechanisms responsible for cold-water disease are not known. It was recently demonstrated that the related fish pathogen, Flavobacterium columnare, requires a functional type IX protein secretion system (T9SS) to cause disease. T9SSs secrete cell surface adhesins, gliding motility proteins, peptidases, and other enzymes, any of which may be virulence factors. The F. psychrophilum genome has genes predicted to encode components of a T9SS. Here, we used a SacB-mediated gene deletion technique recently adapted for use in the Bacteroidetes to delete a core F. psychrophilum T9SS gene, gldN The ΔgldN mutant cells were deficient for secretion of many proteins in comparison to wild-type cells. Complementation of the mutant with wild-type gldN on a plasmid restored secretion. Compared to wild-type and complemented strains, the ΔgldN mutant was deficient in adhesion, gliding motility, and extracellular proteolytic and hemolytic activities. The ΔgldN mutant exhibited reduced virulence in rainbow trout and complementation restored virulence, suggesting that the T9SS plays an important role in the disease.IMPORTANCE Bacterial cold-water disease, caused by F. psychrophilum, is a major problem for salmonid aquaculture. Little is known regarding the virulence factors involved in this disease, and control measures are inadequate. A targeted gene deletion method was adapted to F. psychrophilum and used to demonstrate the importance of the T9SS in virulence. Proteins secreted by this system are likely virulence factors and targets for the development of control measures.

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.

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.

Major Histocompatibility Complex (MHC) Genes and Disease Resistance in Fish

Fascinating about classical major histocompatibility complex (MHC) molecules is their polymorphism. The present study is a review and discussion of the fish MHC situation. The basic pattern of MHC variation in fish is similar to mammals, with MHC class I versus class II, and polymorphic classical versus nonpolymorphic nonclassical. However, in many or all teleost fishes, important differences with mammalian or human MHC were observed: (1) The allelic/haplotype diversification levels of classical MHC class I tend to be much higher than in mammals and involve structural positions within but also outside the peptide binding groove; (2) Teleost fish classical MHC class I and class II loci are not linked. The present article summarizes previous studies that performed quantitative trait loci (QTL) analysis for mapping differences in teleost fish disease resistance, and discusses them from MHC point of view. Overall, those QTL studies suggest the possible importance of genomic regions including classical MHC class II and nonclassical MHC class I genes, whereas similar observations were not made for the genomic regions with the highly diversified classical MHC class I alleles. It must be concluded that despite decades of knowing MHC polymorphism in jawed vertebrate species including fish, firm conclusions (as opposed to appealing hypotheses) on the reasons for MHC polymorphism cannot be made, and that the types of polymorphism observed in fish may not be explained by disease-resistance models alone.

Rainbow trout resistance to bacterial cold water disease: two new quantitative trait loci identified after a natural disease outbreak on a French farm

In rainbow trout farming, Flavobacterium psychrophilum, the causative agent of bacterial cold water disease, is responsible for important economic losses. Resistance to F. psychrophilum is heritable, and several quantitative trait loci (QTL) with moderate effects have been detected, opening up promising perspectives for the genetic improvement of resistance. In most studies however, resistance to F. psychrophilum was assessed in experimental infectious challenges using injection as the infection route, which is not representative of natural infection. Indeed, injection bypasses external barriers, such as mucus and skin, that likely play a protective role against the infection. In this study, we aimed at describing the genetic architecture of the resistance to F. psychrophilum after a natural disease outbreak. In a 2000-fish cohort, reared on a French farm, 720 fish were sampled and genotyped using the medium-throughput Axiom™ Trout Genotyping Array. Overall mortality at the end of the outbreak was 25%. Genome-wide association studies were performed under two different models for time to death measured on 706 fish with validated genotypes for 30 060 SNPs. This study confirms the polygenic inheritance of resistance to F. psychrophilum with a few QTL with moderate effects and a large polygenic background, the heritability of the trait being estimated at 0.34. Two new chromosome-wide significant QTL and three suggestive QTL were detected, each of them explaining between 1% and 4% of genetic variance.