Resistance to viral nervous necrosis in European sea bass (Dicentrarchus labrax L.): heritability and relationships with body weight, cortisol concentration, and antibody titer

Aleutian disease: Risk factors and ImmunAD strategy for genetic improvement of tolerance in American mink (Neogale vison)

Aleutian disease (AD) is a devastating infectious disease in American mink (Neogale vison) industry caused by Aleutian mink disease virus (AMDV). Two crucial steps toward controlling infectious diseases in farm animals are: (i) assessment of the infection risk factors to minimize the likelihood of infection and (ii) selection of animals with superior immune responses against pathogens to build tolerant farms. This study aimed to investigate AD risk factors and evaluate a novel "ImmunAD" approach for genetic improvement of AD tolerance. Phenotypic records and pedigree information of 1,366 and 24,633 animals were included in this study. The risk of animal's age, sex, color type, and year of sampling on AMDV infection was assessed using a logistic regression model and counter immune-electrophoresis (CIEP) test results. ImmunAD phenotype was calculated based on AMDVG enzyme-linked immunosorbent assay (ELISA) and CIEP test results, and breeding values for ImmunAD were estimated using an animal model. Animals were classified into high-coordinated (HCIR), average-coordinated (ACIR), and low-coordinated immune responders (LCIR) using ImmunAD's breeding values, and the impact of selection of HCIR on live grade of pelt quality (PQ), harvest weight (HW), and harvest length (HL) breeding values were evaluated. Age of > 1 year, male sex, and year of sampling were identified as significant risk factors of AD (p < 0.05). A moderate-to-high heritability (0.55±0.07) was estimated for ImmunAD, while a higher heritability was observed among the CIEP-positive animals (0.76±0.06). Significantly higher breeding values were observed for PQ and HL among HCIR than those for LCIR and ACIR (p < 0.05). Our findings indicate the critical role of male breeders in AD distribution within mink farms. Regular screening of AD in male breeders before pairing them with females during breeding seasons can help disease control. ImmunAD strategy can be applied to genetic improvement of AD tolerance, with favorable impacts on some growth and production traits. Higher genetic gains can be achieved in populations with higher AD seroprevalences.

Identification of candidate SNPs and genes associated with resistance to nervous necrosis virus in leopard coral grouper (Plectropomus leopardus) using GWAS

The leopard coral grouper (Plectropomus leopardus), which has become increasingly popular in consumption due to its bright body color and great nutritional, holds a high economic and breeding potential. However, in recent years, the P.leopardus aquaculture industry has been impeded by the nervous necrosis virus (NNV) outbreak, leading to widespread mortality among fry and juvenile grouper. However, the genetic basis of resistance to NNV in P. leopardus remains to be investigated. In the present study, we conducted a genome-wide association analysis (GWAS) on 100 resistant and 100 susceptible samples to discover variants and potential genes linked with NNV resistance. For this study, 157,926 high-quality single nucleotide polymorphisms (SNPs) based on whole genome resequencing were discovered, and eighteen SNPs loci linked to disease resistance were discovered. We annotated six relevant candidate genes, including sik2, herc2, pip5k1c, npr1, mybpc3, and arhgap9, which showed important roles in lipid metabolism, oxidative stress, and neuronal survival. In the brain tissues of resistant and susceptible groups, candidate genes against NNV infection showed significant differential expression. The results indicate that regulating neuronal survival or pathways involved in lipid metabolism may result in increased resistance to NNV. Understanding the molecular mechanisms that lead to NNV resistance will be beneficial for the growth of the P. leopardus breeding sector. Additionally, the identified SNPs could be employed as biomarkers of disease resistance in P. leopardus, which will facilitate the selective breeding of grouper.

Whole-genome sequencing identifies interferon-induced protein IFI6/IFI27-like as a strong candidate gene for VNN resistance in European sea bass

BACKGROUND

Viral nervous necrosis (VNN) is a major disease that affects European sea bass, and understanding the biological mechanisms that underlie VNN resistance is important for the welfare of farmed fish and sustainability of production systems. The aim of this study was to identify genomic regions and genes that are associated with VNN resistance in sea bass.

RESULTS

We generated a dataset of 838,451 single nucleotide polymorphisms (SNPs) identified from whole-genome sequencing (WGS) in the parental generation of two commercial populations (A: 2371 individuals and B: 3428 individuals) of European sea bass with phenotypic records for binary survival in a VNN challenge. For each population, three cohorts were submitted to a red-spotted grouper nervous necrosis virus (RGNNV) challenge by immersion and genotyped on a 57K SNP chip. After imputation of WGS SNPs from their parents, quantitative trait loci (QTL) were mapped using a Bayesian sparse linear mixed model (BSLMM). We found several QTL regions that were specific to one of the populations on different linkage groups (LG), and one 127-kb QTL region on LG12 that was shared by both populations and included the genes ZDHHC14, which encodes a palmitoyltransferase, and IFI6/IFI27-like, which encodes an interferon-alpha induced protein. The most significant SNP in this QTL region was only 1.9 kb downstream of the coding sequence of the IFI6/IFI27-like gene. An unrelated population of four large families was used to validate the effect of the QTL. Survival rates of susceptible genotypes were 40.6% and 45.4% in populations A and B, respectively, while that of the resistant genotype was 66.2% in population B and 78% in population A.

CONCLUSIONS

We have identified a genomic region that carries a major QTL for resistance to VNN and includes the ZDHHC14 and IFI6/IFI27-like genes. The potential involvement of the interferon pathway, a well-known anti-viral defense mechanism in several organisms (chicken, human, or fish), in survival to VNN infection is of particular interest. Our results can lead to major improvements for sea bass breeding programs through marker-assisted genomic selection to obtain more resistant fish.

A Systematic Review and Meta-Analysis of Basal and Post-Stress Circulating Cortisol Concentration in an Important Marine Aquaculture Fish Species, European Sea Bass, Dicentrarchus labrax

BACKGROUND

European sea bass is a species characterized by high and dispersed cortisol levels. The aim of the present study was to analyze all published data on basal and post-acute stress cortisol levels in this species.

METHODS

For this systematic review and meta-analysis the Web of Science and Scopus databases were searched for papers reporting plasma or serum cortisol levels in E. sea bass, without language or date restrictions. Data were extracted directly for the reported results and were analyzed separately for basal and post-acute stress levels, as well their standardized mean differences (SMD) using random-effects meta-analyses.

RESULTS

Of 407 unique records identified, 69 were eligible. Basal cortisol levels had a pooled effect of 88.7 ng mL^-1 (n = 57), while post-acute stress levels were 385.9 ng mL^-1 (n = 34). The average SMD between basal and post-stress was calculated to be 3.02 (n = 22). All analyses had a high between-study heterogeneity. Results for basal and post-stress levels were affected by the assay type and anesthesia prior to blood sampling.

CONCLUSIONS

Cortisol levels in E. sea bass are higher than most studied fish species and display large heterogeneity. Application of stress led to elevated cortisol levels in all studies examined. In all cases, sources of between-studies heterogeneity were identified.

Viral nervous necrosis resistance in gilthead sea bream (Sparus aurata) at the larval stage: heritability and accuracy of genomic prediction with different training and testing settings

BACKGROUND

The gilthead sea bream (Sparus aurata) has long been considered resistant to viral nervous necrosis (VNN), until recently, when significant mortalities caused by a reassortant nervous necrosis virus (NNV) strain were reported. Selective breeding to enhance resistance against NNV might be a preventive action. In this study, 972 sea bream larvae were subjected to a NNV challenge test and the symptomatology was recorded. All the experimental fish and their parents were genotyped using a genome-wide single nucleotide polymorphism (SNP) array consisting of over 26,000 markers.

RESULTS

Estimates of pedigree-based and genomic heritabilities of VNN symptomatology were consistent with each other (0.21, highest posterior density interval at 95% (HPD95%): 0.1-0.4; 0.19, HPD95%: 0.1-0.3, respectively). The genome-wide association study suggested one genomic region, i.e., in linkage group (LG) 23 that might be involved in sea bream VNN resistance, although it was far from the genome-wide significance threshold. The accuracies (r) of the predicted estimated breeding values (EBV) provided by three Bayesian genomic regression models (Bayes B, Bayes C, and Ridge Regression) were consistent and on average were equal to 0.90 when assessed in a set of cross-validation (CV) procedures. When genomic relationships between training and testing sets were minimized, accuracy decreased greatly (r = 0.53 for a validation based on genomic clustering, r = 0.12 for a validation based on a leave-one-family-out approach focused on the parents of the challenged fish). Classification of the phenotype using the genomic predictions of the phenotype or using the genomic predictions of the pedigree-based, all data included, EBV as classifiers was moderately accurate (area under the ROC curve 0.60 and 0.66, respectively).

CONCLUSIONS

The estimate of the heritability for VNN symptomatology indicates that it is feasible to implement selective breeding programs for increased resistance to VNN of sea bream larvae/juveniles. Exploiting genomic information offers the opportunity of developing prediction tools for VNN resistance, and genomic models can be trained on EBV using all data or phenotypes, with minimal differences in classification performance of the trait phenotype. In a long-term view, the weakening of the genomic ties between animals in the training and test sets leads to decreased genomic prediction accuracies, thus periodical update of the reference population with new data is mandatory.

Genetic Basis for Resistance Against Viral Nervous Necrosis: GWAS and Potential of Genomic Prediction Explored in Farmed European Sea Bass (Dicentrarchus labrax)

Viral nervous necrosis (VNN) is an infectious disease caused by the red-spotted grouper nervous necrosis virus (RGNNV) in European sea bass and is considered a serious concern for the aquaculture industry with fry and juveniles being highly susceptible. To understand the genetic basis for resistance against VNN, a survival phenotype through the challenge test against the RGNNV was recorded in populations from multiple year classes (YC2016 and YC2017). A total of 4,851 individuals from 181 families were tested, and a subset (n∼1,535) belonging to 122 families was genotyped using a ∼57K Affymetrix Axiom array. The survival against the RGNNV showed low to moderate heritability with observed scale estimates of 0.18 and 0.25 obtained using pedigree vs. genomic information, respectively. The genome-wide association analysis showed a strong signal of quantitative trait loci (QTL) at LG12 which explained ∼33% of the genetic variance. The QTL region contained multiple genes (ITPK1, PLK4, HSPA4L, REEP1, CHMP2, MRPL35, and SCUBE) with HSPA4L and/or REEP1 genes being highly relevant with a likely effect on host response in managing disease-associated symptoms. The results on the accuracy of predicting breeding values presented 20–43% advantage in accuracy using genomic over pedigree-based information which varied across model types and applied validation schemes.

Transmission Pathways of the VNN Introduced in Croatian Marine Aquaculture

Due to the insufficient capacity of Croatian hatcheries, marine aquaculture depends on the importation of fry from different countries in the Mediterranean basin. Importation enables a risk of spreading pathogenic agents. Viral nervous necrosis (VNN), caused by betanodavirus is devastating for the farming of European sea bass. We described a VNN outbreak that occurred in Croatia in 2014. After the diagnosis of VNN in sea bass fry introduced from the same hatchery to five unconnected marine farms at the Adriatic Coast, we performed surveillance within one of the affected farms. It resulted in proven horizontal spreading of the virus within the farm and to feral fish around farm cages. Real-time RT-PCR tested samples showed the dependence of the virus’ proliferation to the water temperature and the fish age. The highest mortality rates were noted during higher sea temperatures. Phylogenetic analysis of partial sequences of RNA1 and RNA2 supported the hypothesis that the virus was introduced to all studied farms from the same hatchery. Moreover, phylogenetic analysis of the whole genome sequences of infected farmed sea bass and thicklip mullet showed high similarity and it is unlikely that infection in Croatian sea bass farms has originated from wild reservoirs, as the first positive record in wild mullet was recorded after the disease outbreak.

Integrated Management Strategies for Viral Nervous Necrosis (VNN) Disease Control in Marine Fish Farming in the Mediterranean

Viral nervous necrosis (VNN) is the most important viral disease affecting farmed fish in the Mediterranean. VNN can affect multiple fish species in all production phases (broodstock, hatchery, nursery and ongrowing) and sizes, but it is especially severe in larvae and juvenile stages, where can it cause up to 100% mortalities. European sea bass has been and is still the most affected species, and VNN in gilthead sea bream has become an emerging problem in recent years affecting larvae and juveniles and associated to the presence of new nervous necrosis virus (NNV) reassortants. The relevance of this disease as one of the main biological hazards for Mediterranean finfish farming has been particularly addressed in two recent H2020 projects: PerformFISH and MedAID. The presence of the virus in the environment and in the farming systems poses a serious menace for the development of the Mediterranean finfish aquaculture. Several risks associated to the VNN development in farms have been identified in the different phases of the farming system. The main risks concerning VNN affecting gilthead seabream and European seabass have been identified as restocking from wild fish in broodstock facilities, the origin of eggs and juveniles, quality water supply and live food in hatcheries and nurseries, and infected juveniles and location of farms in endemic areas for on-growing sites. Due to the potential severe impact, a holistic integrated management approach is the best strategy to control VNN in marine fish farms. This approach should include continuous surveillance and early and accurate diagnosis, essential for an early intervention when an outbreak occurs, the implementation of biosecurity and disinfection procedures in the production sites and systematic vaccination with effective vaccines. Outbreak management practices, clinical aspects, diagnostic techniques, and disinfections methods are reviewed in detail in this paper. Additionally, new strategies are becoming more relevant, such as the use of genetic resistant lines and boosting the fish immune system though nutrition.

Genomic Predictions of Phenotypes and Pseudo-Phenotypes for Viral Nervous Necrosis Resistance, Cortisol Concentration, Antibody Titer and Body Weight in European Sea Bass

Simple Summary

Selective breeding programs based on genomic data are still not a common practice in aquaculture, although genomic selection has been widely demonstrated to be advantageous when trait phenotyping is a difficult task. In this study, we investigated the accuracy of predicting the phenotype and the estimated breeding value (EBV) of three Bayesian models and a Random Forest algorithm exploiting the information of a genome-wide SNP panel for European sea bass. The genomic predictions were developed for mortality caused by viral nervous necrosis, post-stress cortisol concentration, antibody titer against nervous necrosis virus and body weight. Selective breeding based on genomic data is a possible option for improving these traits while overcoming difficulties related to individual phenotyping of the investigated traits. Our results evidenced that the EBV used as a pseudo-phenotype enhances the predictive performances of genomic models, and that EBV can be predicted with satisfactory accuracy. The genomic prediction of the EBV for mortality might also be used to classify the phenotype for the same trait.

Abstract

In European sea bass (Dicentrarchus labrax L.), the viral nervous necrosis mortality (MORT), post-stress cortisol concentration (HC), antibody titer (AT) against nervous necrosis virus and body weight (BW) show significant heritability, which makes selective breeding a possible option for their improvement. An experimental population (N = 650) generated by a commercial broodstock was phenotyped for the aforementioned traits and genotyped with a genome-wide SNP panel (16,075 markers). We compared the predictive accuracies of three Bayesian models (Bayes B, Bayes C and Bayesian Ridge Regression) and a machine-learning method (Random Forest). The prediction accuracy of the EBV for MORT was approximately 0.90, whereas the prediction accuracies of the EBV and the phenotype were 0.86 and 0.21 for HC, 0.79 and 0.26 for AT and 0.71 and 0.38 for BW. The genomic prediction of the EBV for MORT used to classify the phenotype for the same trait showed moderate classification performance. Genome-wide association studies confirmed the polygenic nature of MORT and demonstrated a complex genetic structure for HC and AT. Genomic predictions of the EBV for MORT could potentially be used to classify the phenotype of the same trait, though further investigations on a larger experimental population are needed.