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.

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.

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

Background

Susceptibility of European sea bass (Dicentrarchus labrax L.) to viral nervous necrosis (VNN) is well-known. Interest towards selective breeding as a tool to enhance genetic resistance in this species has increased sharply due to the major threat represented by VNN for farmed sea bass and limitations concerning specific therapeutical measures. A sea bass experimental population (N = 650) was challenged with nervous necrosis virus (NNV) to investigate genetic variation in VNN mortality. In addition, relationships of this trait with serum cortisol concentration after stress exposure, antibody titer against NNV antigens, and body weight at a fixed age were studied to identify potential indicator traits of VNN resistance.

Results

The estimate of heritability for VNN mortality was moderate and ranged from 0.15 (HPD95%, 95% highest posterior density interval: 0.02, 0.31) to 0.23 (HPD95%: 0.06, 0.47). Heritability estimates for cortisol concentration, antibody titer, and body weight were 0.19 (HPD95%: 0.07, 0.34), 0.36 (HPD95%: 0.16, 0.59) and 0.57 (HPD95%: 0.33, 0.84), respectively. Phenotypic relationships between traits were trivial and not statistically significant, except for the estimated correlation between antibody titer and body weight (0.24). Genetic correlations of mortality with body weight or antibody titer (− 0.39) exhibited a 0.89 probability of being negative. A negligible genetic correlation between mortality and cortisol concentration was detected. Antibody titer was estimated to be positively correlated with body weight (0.49).

Conclusions

Antibody titer against NNV offers the opportunity to use indirect selection to enhance resistance, while the use of cortisol concentration as an indicator trait in breeding programs for VNN resistance is questionable. The estimate of heritability for VNN mortality indicates the feasibility of selective breeding to enhance resistance to NNV and raises attention to the development of genomic prediction tools to simplify testing procedures for selection candidates.

Strong selection pressures maintain divergence on genomic islands in Atlantic cod (Gadus morhua L.) populations

Background

Two distinct populations have been extensively studied in Atlantic cod (Gadus morhua L.): the Northeast Arctic cod (NEAC) population and the coastal cod (CC) population. The objectives of the current study were to identify genomic islands of divergence and to propose an approach to quantify the strength of selection pressures using whole-genome single nucleotide polymorphism (SNP) data. After applying filtering criteria, information on 93 animals (9 CC individuals, 50 NEAC animals and 34 CC × NEAC crossbred individuals) and 3,123,434 autosomal SNPs were used.

Results

Four genomic islands of divergence were identified on chromosomes 1, 2, 7 and 12, which were mapped accurately based on SNP data and which extended in size from 11 to 18 Mb. These regions differed considerably between the two populations although the differences in the rest of the genome were small due to considerable gene flow between the populations. The estimates of selection pressures showed that natural selection was substantially more important than genetic drift in shaping these genomic islands. Our data confirmed results from earlier publications that suggested that genomic islands are due to chromosomal rearrangements that are under strong selection and reduce recombination between rearranged and non-rearranged segments.

Conclusions

Our findings further support the hypothesis that selection and reduced recombination in genomic islands may promote speciation between these two populations although their habitats overlap considerably and migrations occur between them.

Characterization of a male specific region containing a candidate sex determining gene in Atlantic cod

The genetic mechanisms determining sex in teleost fishes are highly variable and the master sex determining gene has only been identified in few species. Here we characterize a male-specific region of 9 kb on linkage group 11 in Atlantic cod (Gadus morhua) harboring a single gene named zkY for zinc knuckle on the Y chromosome. Diagnostic PCR test of phenotypically sexed males and females confirm the sex-specific nature of the Y-sequence. We identified twelve highly similar autosomal gene copies of zkY, of which eight code for proteins containing the zinc knuckle motif. 3D modeling suggests that the amino acid changes observed in six copies might influence the putative RNA-binding specificity. Cod zkY and the autosomal proteins zk1 and zk2 possess an identical zinc knuckle structure, but only the Y-specific gene zkY was expressed at high levels in the developing larvae before the onset of sex differentiation. Collectively these data suggest zkY as a candidate master masculinization gene in Atlantic cod. PCR amplification of Y-sequences in Arctic cod (Arctogadus glacialis) and Greenland cod (Gadus macrocephalus ogac) suggests that the male-specific region emerged in codfishes more than 7.5 million years ago.

Regional genetic diversity for NNV grouper viruses across the Indo-Asian region - implications for selecting virus resistance in farmed groupers

Grouper aquaculture around Asia is impacted by the nervous necrosis virus (NNV) and, in response, host resistance to this infection is being considered as a trait for selection. However efficient selection may be confounded if there are different genetic strains of NNV within and between regions and over years. This study uses statistical approaches and assessment of “characteristic attributes” (i.e. nucleotide positions that discriminate among strains) to assess whether published and new NNV RNA2 cds sequences show genetic differentiation over geography, host species and years. Rather clear evidence was found for regional strains of NNV. Interestingly, most of the geographic defining “characteristic attributes” were in codon position three, and not translated into differences for the protein capsid (i.e. they were synonymous variations), suggesting that while NNV strains were geographically isolated and had diverged in different regions for RNA sequences, selection had largely conserved the protein sequences among regions. The apparent selection constraint on the capsid protein may mitigate the risk that despite geographic subdivision, NNV strain variability will confound genetic selection for host resistance. The existence of regional Asian NNV strains may suggest that hatcheries are at risk from NNV not only from imported material but also from endemic reservoirs.

Viral encephalopathy and retinopathy in aquaculture: a review

Viral encephalopathy and retinopathy (VER), otherwise known as viral nervous necrosis (VNN), is a major devastating threat for aquatic animals. Betanodaviruses have been isolated in at least 70 aquatic animal species in marine and in freshwater environments throughout the world, with the notable exception of South America. In this review, the main features of betanodavirus, including its diversity, its distribution and its transmission modes in fish, are firstly presented. Then, the existing diagnosis and detection methods, as well as the different control procedures of this disease, are reviewed. Finally, the potential of selective breeding, including both conventional and genomic selection, as an opportunity to obtain resistant commercial populations, is examined.

Genomic predictions can accelerate selection for resistance against Piscirickettsia salmonis in Atlantic salmon (Salmo salar)

BACKGROUND

Salmon Rickettsial Syndrome (SRS) caused by Piscirickettsia salmonis is a major disease affecting the Chilean salmon industry. Genomic selection (GS) is a method wherein genome-wide markers and phenotype information of full-sibs are used to predict genomic EBV (GEBV) of selection candidates and is expected to have increased accuracy and response to selection over traditional pedigree based Best Linear Unbiased Prediction (PBLUP). Widely used GS methods such as genomic BLUP (GBLUP), SNPBLUP, Bayes C and Bayesian Lasso may perform differently with respect to accuracy of GEBV prediction. Our aim was to compare the accuracy, in terms of reliability of genome-enabled prediction, from different GS methods with PBLUP for resistance to SRS in an Atlantic salmon breeding program. Number of days to death (DAYS), binary survival status (STATUS) phenotypes, and 50 K SNP array genotypes were obtained from 2601 smolts challenged with P. salmonis. The reliability of different GS methods at different SNP densities with and without pedigree were compared to PBLUP using a five-fold cross validation scheme.

RESULTS

Heritability estimated from GS methods was significantly higher than PBLUP. Pearson's correlation between predicted GEBV from PBLUP and GS models ranged from 0.79 to 0.91 and 0.79-0.95 for DAYS and STATUS, respectively. The relative increase in reliability from different GS methods for DAYS and STATUS with 50 K SNP ranged from 8 to 25% and 27-30%, respectively. All GS methods outperformed PBLUP at all marker densities. DAYS and STATUS showed superior reliability over PBLUP even at the lowest marker density of 3 K and 500 SNP, respectively. 20 K SNP showed close to maximal reliability for both traits with little improvement using higher densities.

CONCLUSIONS

These results indicate that genomic predictions can accelerate genetic progress for SRS resistance in Atlantic salmon and implementation of this approach will contribute to the control of SRS in Chile. We recommend GBLUP for routine GS evaluation because this method is computationally faster and the results are very similar with other GS methods. The use of lower density SNP or the combination of low density SNP and an imputation strategy may help to reduce genotyping costs without compromising gain in reliability.

The polymorphisms of a MIF gene and their association with Vibrio resistance in the clam Meretrix meretrix

Macrophage migration inhibitory factor (MIF) is an important proinflammatory cytokine that mediates both innate and adaptive immune responses. In this study, a homolog of MIF was identified in the clam Meretrix meretrix. Ten SNPs in the DNA partial sequence of MmMIF were found to be significantly associated with Vibrio resistance (P < 0.05). Distinct expression patterns of MmMIF among different haplotypes were observed after Vibrio challenge. The results showed that haplotypes did not affect MmMIF expression in the negative control group, while the expression of MmMIF in clams with Hap1 and Hap1/Hap2 was significantly lower than that with Hap2 at 24 h in the PBS-injected group but significantly higher than that with Hap2 in the Vibrio-injected group. The results also indicate that Hap1 and Hap1/Hap2 can specifically respond to mechanical stimulation while Hap2 can specifically respond to Vibrio infection. The effect of a missense mutation was detected by site-directed mutagenesis using fusion expression of the protein, which showed that the SNP g.737 (I > V) has no effect on redox activity and tautomerase activity. These studies identified a potential marker that is enriched in Vibrio-resistant clams that can be used for the genetic breeding of Meretrix meretrix.

Two adjacent inversions maintain genomic differentiation between migratory and stationary ecotypes of Atlantic cod

Atlantic cod is composed of multiple migratory and stationary populations widely distributed in the North Atlantic Ocean. The Northeast Arctic cod (NEAC) population in the Barents Sea undertakes annual spawning migrations to the northern Norwegian coast. Although spawning occurs sympatrically with the stationary Norwegian coastal cod (NCC), phenotypic and genetic differences between NEAC and NCC are maintained. In this study, we resolve the enigma by revealing the mechanisms underlying these differences. Extended linkage disequilibrium (LD) and population divergence were demonstrated in a 17.4-Mb region on linkage group 1 (LG1) based on genotypes of 494 SNPs from 192 parents of farmed families of NEAC, NCC or NEACxNCC crosses. Linkage analyses revealed two adjacent inversions within this region that repress meiotic recombination in NEACxNCC crosses. We identified a NEAC-specific haplotype consisting of 186 SNPs that was fixed in NEAC sampled from the Barents Sea, but segregating under Hardy-Weinberg equilibrium in eight NCC stocks. Comparative genomic analyses determine the NEAC configuration of the inversions to be the derived state and date it to ~1.6-2.0 Mya. The haplotype block harbours 763 genes, including candidates regulating swim bladder pressure, haem synthesis and skeletal muscle organization conferring adaptation to long-distance migrations and vertical movements down to large depths. Our results suggest that the migratory ecotype experiences strong directional selection for the two adjacent inversions on LG1. Despite interbreeding between NEAC and NCC, the inversions are maintaining genetic differentiation, and we hypothesize the co-occurrence of multiple adaptive alleles forming a 'supergene' in the NEAC population.

Transcriptome analysis of medaka following epinecidin-1 and TH1-5 treatment of NNV infection

Nervous necrosis virus (NNV) infects a wide range of larval and juvenile fish species, thereby causing enormous economic losses in the aquaculture industry. Possible solutions to this problem include the use of antimicrobial peptides (AMPs), which directly inhibit bacterial growth, and also modulate host signaling mechanisms. The AMPs epinecidin (Epi)-1 and Tilapia hepcidin (TH) 1-5 have been demonstrated to be effective against Nervous necrosis virus infection in medaka (Oryzias latipes). However, the underlying molecular mechanisms are yet to be explored. Here, microarray analyses were performed to examine how NNV infection and/or epinecidin-1 or TH1-5 treatment affects gene expression in medaka; such analyses enabled the prediction of host signaling pathways affected by virus infection and/or regulated by epinecidin-1 and TH1-5. Transcriptome analysis revealed altered expression of genes involved in B cell activation, T cell activation, adipocytokine signaling, and mast cell activation. We subsequently used real-time PCR to analyze expression of key genes involved in these signaling mechanisms. Medaka infected with NNV exhibited up-regulation of PVALB, CEBPA, IFIM, IFN, IL-6ST, NF-kB2, SOC3, SP1, and TGFB1, and such increases were prevented by pre-treatment with epinecidin-1 or TH1-5. Immunohistochemistry using the anti-NNV antibody to stain brain and eye sections revealed that epinecidin-1 treatment during or after infection clears viral load, while TH1-5 treatment only reduces viral numbers if applied during infection. These observations demonstrate that epinecidin-1 and TH1-5 modulate NNV-induced host signaling mechanisms, thereby preventing viral multiplication in host organisms.