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

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.

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.

Applying genetic technologies to combat infectious diseases in aquaculture

Disease and parasitism cause major welfare, environmental and economic concerns for global aquaculture. In this review, we examine the status and potential of technologies that exploit genetic variation in host resistance to tackle this problem. We argue that there is an urgent need to improve understanding of the genetic mechanisms involved, leading to the development of tools that can be applied to boost host resistance and reduce the disease burden. We draw on two pressing global disease problems as case studies-sea lice infestations in salmonids and white spot syndrome in shrimp. We review how the latest genetic technologies can be capitalised upon to determine the mechanisms underlying inter- and intra-species variation in pathogen/parasite resistance, and how the derived knowledge could be applied to boost disease resistance using selective breeding, gene editing and/or with targeted feed treatments and vaccines. Gene editing brings novel opportunities, but also implementation and dissemination challenges, and necessitates new protocols to integrate the technology into aquaculture breeding programmes. There is also an ongoing need to minimise risks of disease agents evolving to overcome genetic improvements to host resistance, and insights from epidemiological and evolutionary models of pathogen infestation in wild and cultured host populations are explored. Ethical issues around the different approaches for achieving genetic resistance are discussed. Application of genetic technologies and approaches has potential to improve fundamental knowledge of mechanisms affecting genetic resistance and provide effective pathways for implementation that could lead to more resistant aquaculture stocks, transforming global aquaculture.

Genetic and epigenetic regulation of growth, reproduction, disease resistance and stress responses in aquaculture

Major progress has been made with genomic and genetic studies in aquaculture in the last decade. However, research on epigenetic regulation of aquaculture traits is still at an early stage. It is apparent that most, if not all, aquaculture traits are regulated at both genetic and epigenetic levels. This paper reviews recent progress in understanding of genetic and epigenetic regulation of important aquaculture traits such as growth, reproduction, disease resistance, and stress responses. Although it is challenging to make generalized statements, DNA methylation is mostly correlated with down-regulation of gene expression, especially when at promoters and enhancers. As such, methylation of growth factors and their receptors is negatively correlated with growth; hypomethylation of genes important for stress tolerance is correlated with increased stress tolerance; hypomethylation of genes important for male or female sex differentiation leads to sex differentiation into males or females, respectively. It is apparent that environmental regulation of aquaculture traits is mediated at the level of epigenetic regulation, and such environment-induced epigenetic changes appeared to be intergenerationally inherited, but evidences for transgenerational inheritance are still limited.

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.

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

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

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.

Optimization of Genomic Selection to Improve Disease Resistance in Two Marine Fishes, the European Sea Bass (Dicentrarchus labrax) and the Gilthead Sea Bream (Sparus aurata)

Disease outbreaks are a major threat to the aquaculture industry, and can be controlled by selective breeding. With the development of high-throughput genotyping technologies, genomic selection may become accessible even in minor species. Training population size and marker density are among the main drivers of the prediction accuracy, which both have a high impact on the cost of genomic selection. In this study, we assessed the impact of training population size as well as marker density on the prediction accuracy of disease resistance traits in European sea bass (Dicentrarchus labrax) and gilthead sea bream (Sparus aurata). We performed a challenge to nervous necrosis virus (NNV) in two sea bass cohorts, a challenge to Vibrio harveyi in one sea bass cohort and a challenge to Photobacterium damselae subsp. piscicida in one sea bream cohort. Challenged individuals were genotyped on 57K-60K SNP chips. Markers were sampled to design virtual SNP chips of 1K, 3K, 6K, and 10K markers. Similarly, challenged individuals were randomly sampled to vary training population size from 50 to 800 individuals. The accuracy of genomic-based (GBLUP model) and pedigree-based estimated breeding values (EBV) (PBLUP model) was computed for each training population size using Monte-Carlo cross-validation. Genomic-based breeding values were also computed using the virtual chips to study the effect of marker density. For resistance to Viral Nervous Necrosis (VNN), as one major QTL was detected, the opportunity of marker-assisted selection was investigated by adding a QTL effect in both genomic and pedigree prediction models. As training population size increased, accuracy increased to reach values in range of 0.51-0.65 for full density chips. The accuracy could still increase with more individuals in the training population as the accuracy plateau was not reached. When using only the 6K density chip, accuracy reached at least 90% of that obtained with the full density chip. Adding the QTL effect increased the accuracy of the PBLUP model to values higher than the GBLUP model without the QTL effect. This work sets a framework for the practical implementation of genomic selection to improve the resistance to major diseases in European sea bass and gilthead sea bream.

Genome-Wide Association Analysis Reveals the Genetic Architecture of Parasite (Cryptocaryon irritans) Resistance in Large Yellow Croaker (Larimichthys crocea)

Large yellow croaker is an important marine culture species in China. Recently, the large yellow croaker industry is threatened by various disease problems, especially for the white spot disease, which is caused by parasite Cryptocaryon irritans. In the current study, we conducted a genome-wide association study (GWAS) for C. irritans resistance in two large yellow croaker populations (n = 264 and n = 480, respectively). We identified 15 QTL with explained genetic variance ranging from 1 to 8% in the two populations. One QTL on chromosome 23 was shared by the two populations, and three QTL had been reported in the previous study. We identified a lot of biological pathways associated with C. irritans resistance, such as hormone transport, response to bacterium, apoptotic process, acute inflammatory response to antigenic stimulus, and NF-kappa B signaling pathway. The genes casp8 and traf6 involved in regulatory network for apoptosis and inflammation were identified to be candidate genes for C. irritans resistance. Our results showed the complex polygenic architecture of resistance of large yellow croaker against C. irritans. These results would be helpful for the researches of the molecular mechanism of C. irritans resistance and genome-assisted breeding of large yellow croaker.

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.