Detection and confirmation of a major QTL affecting resistance to infectious pancreatic necrosis (IPN) in Atlantic salmon (Salmo salar)

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.

Estimating individuals' genetic and non-genetic effects underlying infectious disease transmission from temporal epidemic data

Individuals differ widely in their contribution to the spread of infection within and across populations. Three key epidemiological host traits affect infectious disease spread: susceptibility (propensity to acquire infection), infectivity (propensity to transmit infection to others) and recoverability (propensity to recover quickly). Interventions aiming to reduce disease spread may target improvement in any one of these traits, but the necessary statistical methods for obtaining risk estimates are lacking. In this paper we introduce a novel software tool called SIRE (standing for "Susceptibility, Infectivity and Recoverability Estimation"), which allows for the first time simultaneous estimation of the genetic effect of a single nucleotide polymorphism (SNP), as well as non-genetic influences on these three unobservable host traits. SIRE implements a flexible Bayesian algorithm which accommodates a wide range of disease surveillance data comprising any combination of recorded individual infection and/or recovery times, or disease diagnostic test results. Different genetic and non-genetic regulations and data scenarios (representing realistic recording schemes) were simulated to validate SIRE and to assess their impact on the precision, accuracy and bias of parameter estimates. This analysis revealed that with few exceptions, SIRE provides unbiased, accurate parameter estimates associated with all three host traits. For most scenarios, SNP effects associated with recoverability can be estimated with highest precision, followed by susceptibility. For infectivity, many epidemics with few individuals give substantially more statistical power to identify SNP effects than the reverse. Importantly, precise estimates of SNP and other effects could be obtained even in the case of incomplete, censored and relatively infrequent measurements of individuals' infection or survival status, albeit requiring more individuals to yield equivalent precision. SIRE represents a new tool for analysing a wide range of experimental and field disease data with the aim of discovering and validating SNPs and other factors controlling infectious disease transmission.

Genome-Wide Association Mapping and Gene Expression Analyses Reveal Genetic Mechanisms of Disease Resistance Variations in Cynoglossus semilaevis

The sustainable development of aquaculture has been impeded by infectious diseases worldwide. However, the genomic architecture and the genetic basis underlying the disease resistance remain poorly understood, which severely hampers both the understanding of the evolution of fish disease resistance traits and the prevention of these diseases in the aquaculture community. Cynoglossus semilaevis is a representative and commercially-important flatfish species. Here we combined genome-wide association study and Fst and nucleotide diversity filtration to identify loci important for the disease resistance. Based on 1,016,774 single-nucleotide polymorphisms (SNPs) identified from 650 Gb genome resequencing data of 505 individuals, we detected 33 SNPs significantly associated with disease resistance and 79 candidate regions after filtration steps. Both the allele frequencies and genotype frequencies of the associated loci were significantly different between the resistant and susceptible fish, suggesting a role in the genetic basis of disease resistance. The SNP with strongest association with disease resistance was located in Chr 17, at 145 bp upstream of fblx19 gene, and overlapped with the major quantitative trait locus previously identified. Several genes, such as plekha7, nucb2, and fgfr2, were also identified to potentially play roles in the disease resistance. Furthermore, the expression of some associating genes were likely under epigenetic regulations between the bacterial resistant and susceptible families. These results provide insights into the mechanism that enable variation of disease resistance to bacterial pathogen infection. The identified polymorphisms and genes are valuable targets and molecular resources for disease resistance and other traits, and for advanced breeding practice for superior germplasm in fish aquaculture.

Genome-Wide Association Analysis With a 50K Transcribed Gene SNP-Chip Identifies QTL Affecting Muscle Yield in Rainbow Trout

Detection of coding/functional SNPs that change the biological function of a gene may lead to identification of putative causative alleles within QTL regions and discovery of genetic markers with large effects on phenotypes. This study has two-fold objectives, first to develop, and validate a 50K transcribed gene SNP-chip using RNA-Seq data. To achieve this objective, two bioinformatics pipelines, GATK and SAMtools, were used to identify ~21K transcribed SNPs with allelic imbalances associated with important aquaculture production traits including body weight, muscle yield, muscle fat content, shear force, and whiteness in addition to resistance/susceptibility to bacterial cold-water disease (BCWD). SNPs ere identified from pooled RNA-Seq data collected from ~620 fish, representing 98 families from growth- and 54 families from BCWD-selected lines with divergent phenotypes. In addition, ~29K transcribed SNPs without allelic-imbalances were strategically added to build a 50K Affymetrix SNP-chip. SNPs selected included two SNPs per gene from 14K genes and ~5K non-synonymous SNPs. The SNP-chip was used to genotype 1728 fish. The average SNP calling-rate for samples passing quality control (QC; 1,641 fish) was ≥ 98.5%. The second objective of this study was to test the feasibility of using the new SNP-chip in GWA (Genome-wide association) analysis to identify QTL explaining muscle yield variance. GWA study on 878 fish (representing 197 families from 2 consecutive generations) with muscle yield phenotypes and genotyped for 35K polymorphic markers (passing QC) identified several QTL regions explaining together up to 28.40% of the additive genetic variance for muscle yield in this rainbow trout population. The most significant QTLs were on chromosomes 14 and 16 with 12.71 and 10.49% of the genetic variance, respectively. Many of the annotated genes in the QTL regions were previously reported as important regulators of muscle development and cell signaling. No major QTLs were identified in a previous GWA study using a 57K genomic SNP chip on the same fish population. These results indicate improved detection power of the transcribed gene SNP-chip in the target trait and population, allowing identification of large-effect QTLs for important traits in rainbow trout.

Teleosts Genomics: Progress and Prospects in Disease Prevention and Control

Genome wide studies based on conventional molecular tools and upcoming omics technologies are beginning to gain functional applications in the control and prevention of diseases in teleosts fish. Herein, we provide insights into current progress and prospects in the use genomics studies for the control and prevention of fish diseases. Metagenomics has emerged to be an important tool used to identify emerging infectious diseases for the timely design of rational disease control strategies, determining microbial compositions in different aquatic environments used for fish farming and the use of host microbiota to monitor the health status of fish. Expounding the use of antimicrobial peptides (AMPs) as therapeutic agents against different pathogens as well as elucidating their role in tissue regeneration is another vital aspect of genomics studies that had taken precedent in recent years. In vaccine development, prospects made include the identification of highly immunogenic proteins for use in recombinant vaccine designs as well as identifying gene signatures that correlate with protective immunity for use as benchmarks in optimizing vaccine efficacy. Progress in quantitative trait loci (QTL) mapping is beginning to yield considerable success in identifying resistant traits against some of the highly infectious diseases that have previously ravaged the aquaculture industry. Altogether, the synopsis put forth shows that genomics studies are beginning to yield positive contribution in the prevention and control of fish diseases in aquaculture.

Genomic Prediction Accuracy for Resistance Against Piscirickettsia salmonis in Farmed Rainbow Trout

Salmonid rickettsial syndrome (SRS), caused by the intracellular bacterium Piscirickettsia salmonis, is one of the main diseases affecting rainbow trout (Oncorhynchus mykiss) farming. To accelerate genetic progress, genomic selection methods can be used as an effective approach to control the disease. The aims of this study were: (i) to compare the accuracy of estimated breeding values using pedigree-based best linear unbiased prediction (PBLUP) with genomic BLUP (GBLUP), single-step GBLUP (ssGBLUP), Bayes C, and Bayesian Lasso (LASSO); and (ii) to test the accuracy of genomic prediction and PBLUP using different marker densities (0.5, 3, 10, 20, and 27 K) for resistance against P. salmonis in rainbow trout. Phenotypes were recorded as number of days to death (DD) and binary survival (BS) from 2416 fish challenged with P. salmonis. A total of 1934 fish were genotyped using a 57 K single-nucleotide polymorphism (SNP) array. All genomic prediction methods achieved higher accuracies than PBLUP. The relative increase in accuracy for different genomic models ranged from 28 to 41% for both DD and BS at 27 K SNP. Between different genomic models, the highest relative increase in accuracy was obtained with Bayes C (∼40%), where 3 K SNP was enough to achieve a similar accuracy to that of the 27 K SNP for both traits. For resistance against P. salmonis in rainbow trout, we showed that genomic predictions using GBLUP, ssGBLUP, Bayes C, and LASSO can increase accuracy compared with PBLUP. Moreover, it is possible to use relatively low-density SNP panels for genomic prediction without compromising accuracy predictions for resistance against P. salmonis in rainbow trout.

Microsatellites Associated with Growth Performance and Analysis of Resistance to Aeromonas hydrophila in Tambaqui Colossoma macropomum

Tambaqui, Colossoma macropomum, is the main native fish species produced in Brazil, and is an important species for genetic improvement in aquaculture. In addition, breeding studies on this species can be optimized with the use of molecular markers associated with productive phenotypes. The objective of the present study was to test the performance of growth traits and resistance to the bacteria, Aeromonas hydrophila, in association with microsatellite markers in C. macropomum. In this study, three full-sib families were subjected to bacterial challenge and morphometric growth assessments. Tambaqui families subjected to the bacterial challenge differed significantly in death time and mortality rate. There was, however, no association between resistance to bacteria and microsatellite markers. In relation to growth traits, we observed a marker/phenotype association in two microsatellites. The marker in the 6b isoform x5 gene (TNCRC6b) was associated with length, whereas an anonymous marker was associated with height. The present study highlighted the evaluation of molecular markers associated with growth traits, and can serve as the basis for future marker-assisted selection (MAS) of tambaqui.

Infectious Pancreatic Necrosis Virus Causing Clinical and Subclinical Infections in Atlantic Salmon Have Different Genetic Fingerprints

Infectious pancreatic necrosis virus (IPNV) is the causative agent of IPN, an important disease of salmonids. IPNV infections result in either sub-clinical or overt disease and the basis of this difference is not well-understood. The objective of the present study was to determine the VP2 gene of the virus associated with the different forms of clinical manifestation. Groups of Atlantic salmon (Salmo salar L.) reared in farms located in different IPN disease pressures were monitored from brood stock until grow-out over a 3 year period. Hatcheries A1 and B1 as well as cooperating seawater farms were located in a low disease risk area while hatcheries A2 and B2 as well as their cooperating seawater farms were in high IPN risk areas. Samples including eggs, milt, whole fry, kidney depending on the stage of production were collected during outbreaks or in apparently healthy populations where no outbreaks occurred. The virus was re-isolated in CHSE cells and the VP2 gene amplified by RT-PCR followed by sequencing. During the freshwater stage, there were no disease outbreaks at hatcheries A1, A2, and B1 (except in one fish group that originated from hatchery B2), although IPNV was isolated from some of the fish groups at all 3 hatcheries. By contrast, all fish groups at hatchery B2 suffered IPN outbreaks. In seawater, only groups of fish originating from hatchery A1 had no IPN outbreaks albeit virus being isolated from the fish. On the other hand, fish originating from hatcheries A2, B1, and B2 experienced outbreaks in seawater. The VP2 amino acid fingerprint of the virus associated with subclinical infections from A1 and co-operating seawater sites was V64A137P217T221A247N252S281D282E319. By contrast, all virus isolates associated with clinical infections had the motif I64T137T217A221T247V252T281N282A319, where underlined amino acids represent the avirulent and highly virulent motif, respectively. Phylogenetic analysis of amino acid sequences showed 2 clades, one of isolates associated with subclinical infections (from A1 and cooperating seawater farms) and the other of isolates from fish with overt disease (all other sites). Furthermore, the clustering pattern of isolates suggests more circulation of virus within fish groups rather than between them.

Fine mapping QTL for resistance to VNN disease using a high-density linkage map in Asian seabass

Asian seabass has suffered from viral nervous necrosis (VNN) disease. Our previous study has mapped quantitative trait loci (QTL) for resistance to VNN disease. To fine map these QTL and identify causative genes, we identified 6425 single nucleotide polymorphisms (SNPs) from 85 dead and 94 surviving individuals. Combined with 155 microsatellites, we constructed a genetic map consisting of 24 linkage groups (LGs) containing 3000 markers, with an average interval of 1.27 cM. We mapped one significant and three suggestive QTL with phenotypic variation explained (PVE) of 8.3 to 11.0%, two significant and two suggestive QTL with PVE of 7.8 to 10.9%, for resistance in three LGs and survival time in four LGs, respectively. Further analysis one QTL with the largest effect identified protocadherin alpha-C 2-like (Pcdhac2) as the possible candidate gene. Association study in 43 families with 1127 individuals revealed a 6 bp insertion-deletion was significantly associated with disease resistance. qRT-PCR showed the expression of Pcdhac2 was significantly induced in the brain, muscle and skin after nervous necrosis virus (NNV) infection. Our results could facilitate marker-assisted selection (MAS) for resistance to NNV in Asian seabass and set up the basis for functional analysis of the potential causative gene for resistance.

Gene expression comparison of resistant and susceptible Atlantic salmon fry challenged with Infectious Pancreatic Necrosis virus reveals a marked contrast in immune response

BACKGROUND

Infectious Pancreatic Necrosis (IPN) is a highly contagious birnavirus disease of farmed salmonid fish, which often causes high levels of morbidity and mortality. A large host genetic component to resistance has been previously described for Atlantic salmon (Salmo salar L.), which mediates high mortality rates in some families and zero mortality in others. However, the molecular and immunological basis for this resistance is not yet fully known. This manuscript describes a global comparison of the gene expression profiles of resistant and susceptible Atlantic salmon fry following challenge with the IPN virus.

RESULTS

Salmon fry from two IPNV-resistant and two IPNV-susceptible full sibling families were challenged with the virus and sampled at 1 day, 7 days and 20 days post-challenge. Significant viral titre was observed in both resistant and susceptible fish at all timepoints, although generally at higher levels in susceptible fish. Gene expression profiles combined with gene ontology and pathway analyses demonstrated that while a clear immune response was observed in both resistant and susceptible fish, there were striking differences between the two phenotypes. The susceptible fish showed marked up-regulation of genes related to cytokine activity and inflammatory response that evidently failed to protect against the virus. In contrast, the resistant fish demonstrated a less pronounced immune response including up-regulation of genes relating to the M2 macrophage system.

CONCLUSIONS

While only the susceptible phenotype shows appreciable mortality levels, both resistant and susceptible fish can become infected with IPNV. Susceptible fish are characterized by a much larger, yet ineffective, immune response, largely related to cytokine and inflammatory systems. Resistant fish demonstrate a more moderate, putative macrophage-mediated inflammatory response, which may contribute to their survival.

A genome-wide association study in catfish reveals the presence of functional hubs of related genes within QTLs for columnaris disease resistance

Background

Columnaris causes severe mortalities among many different wild and cultured freshwater fish species, but understanding of host resistance is lacking. Catfish, the primary aquaculture species in the United States, serves as a great model for the analysis of host resistance against columnaris disease. Channel catfish in general is highly resistant to the disease while blue catfish is highly susceptible. F2 generation of hybrids can be produced where phenotypes and genotypes are segregating, providing a useful system for QTL analysis. To identify genes associated with columnaris resistance, we performed a genome-wide association study (GWAS) using the catfish 250 K SNP array with 340 backcross progenies derived from crossing female channel catfish (Ictalurus punctatus) with male F1 hybrid catfish (female channel catfish I. punctatus × male blue catfish I. furcatus).

Results

A genomic region on linkage group 7 was found to be significantly associated with columnaris resistance. Within this region, five have known functions in immunity, including pik3r3b, cyld-like, adcyap1r1, adcyap1r1-like, and mast2. In addition, 3 additional suggestively associated QTL regions were identified on linkage groups 7, 12, and 14. The resistant genotypes on the QTLs of linkage groups 7 and 12 were found to be homozygous with both alleles being derived from channel catfish. The paralogs of the candidate genes in the suggestively associated QTL of linkage group 12 were found on the QTLs of linkage group 7. Many candidate genes on the four associated regions are involved in PI3K pathway that is known to be required by many bacteria for efficient entry into the host.

Conclusion

The GWAS revealed four QTLs associated with columnaris resistance in catfish. Strikingly, the candidate genes may be arranged as functional hubs; the candidate genes within the associated QTLs on linkage groups 7 and 12 are not only co-localized, but also functionally related, with many of them being involved in the PI3K signal transduction pathway, suggesting its importance for columnaris resistance.

Genetics and genomics of disease resistance in salmonid species

Infectious and parasitic diseases generate large economic losses in salmon farming. A feasible and sustainable alternative to prevent disease outbreaks may be represented by genetic improvement for disease resistance. To include disease resistance into the breeding goal, prior knowledge of the levels of genetic variation for these traits is required. Furthermore, the information from the genetic architecture and molecular factors involved in resistance against diseases may be used to accelerate the genetic progress for these traits. In this regard, marker assisted selection and genomic selection are approaches which incorporate molecular information to increase the accuracy when predicting the genetic merit of selection candidates. In this article we review and discuss key aspects related to disease resistance in salmonid species, from both a genetic and genomic perspective, with emphasis in the applicability of disease resistance traits into breeding programs in salmonids.

Assessment of genetic correlation between bacterial cold water disease resistance and spleen index in a domesticated population of rainbow trout: identification of QTL on chromosome Omy19

Selective breeding of animals for increased disease resistance is an effective strategy to reduce mortality in aquaculture. However, implementation of selective breeding programs is limited by an incomplete understanding of host resistance traits. We previously reported results of a rainbow trout selection program that demonstrated increased survival following challenge with Flavobacterium psychrophilum, the causative agent of bacterial cold water disease (BCWD). Mechanistic study of disease resistance identified a positive phenotypic correlation between post-challenge survival and spleen somatic-index (SI). Herein, we investigated the hypothesis of a genetic correlation between the two traits influenced by colocalizing QTL. We evaluated the inheritance and calculated the genetic correlation in five year-classes of odd- and even-year breeding lines. A total of 322 pedigreed families (n = 25,369 fish) were measured for disease resistance, and 251 families (n = 5,645 fish) were evaluated for SI. Spleen index was moderately heritable in both even-year (h(2)  = 0.56±0.18) and odd-year (h(2)  = 0.60±0.15) lines. A significant genetic correlation between SI and BCWD resistance was observed in the even-year line (rg  = 0.45±0.20, P = 0.03) but not in the odd-year line (rg  = 0.16±0.12, P = 0.19). Complex segregation analyses of the even-year line provided evidence of genes with major effect on SI, and a genome scan of a single family, 2008132, detected three significant QTL on chromosomes Omy19, 16 and 5, in addition to ten suggestive QTL. A separate chromosome scan for disease resistance in family 2008132 identified a significant BCWD QTL on Omy19 that was associated with time to death and percent survival. In family 2008132, Omy19 microsatellite alleles that associated with higher disease resistance also associated with increased spleen size raising the hypothesis that closely linked QTL contribute to the correlation between these traits. To our knowledge, this is the first estimation of spleen size heritability and evidence for genetic linkage with specific disease resistance in a teleost fish.

Characterisation of QTL-linked and genome-wide restriction site-associated DNA (RAD) markers in farmed Atlantic salmon

BACKGROUND

Restriction site-associated DNA sequencing (RAD-Seq) is a genome complexity reduction technique that facilitates large-scale marker discovery and genotyping by sequencing. Recent applications of RAD-Seq have included linkage and QTL mapping with a particular focus on non-model species. In the current study, we have applied RAD-Seq to two Atlantic salmon families from a commercial breeding program. The offspring from these families were classified into resistant or susceptible based on survival/mortality in an Infectious Pancreatic Necrosis (IPN) challenge experiment, and putative homozygous resistant or susceptible genotype at a major IPN-resistance QTL. From each family, the genomic DNA of the two heterozygous parents and seven offspring of each IPN phenotype and genotype was digested with the SbfI enzyme and sequenced in multiplexed pools.

RESULTS

Sequence was obtained from approximately 70,000 RAD loci in both families and a filtered set of 6,712 segregating SNPs were identified. Analyses of genome-wide RAD marker segregation patterns in the two families suggested SNP discovery on all 29 Atlantic salmon chromosome pairs, and highlighted the dearth of male recombination. The use of pedigreed samples allowed us to distinguish segregating SNPs from putative paralogous sequence variants resulting from the relatively recent genome duplication of salmonid species. Of the segregating SNPs, 50 were linked to the QTL. A subset of these QTL-linked SNPs were converted to a high-throughput assay and genotyped across large commercial populations of IPNV-challenged salmon fry. Several SNPs showed highly significant linkage and association with resistance to IPN, and population linkage-disequilibrium-based SNP tests for resistance were identified.

CONCLUSIONS

We used RAD-Seq to successfully identify and characterise high-density genetic markers in pedigreed aquaculture Atlantic salmon. These results underline the effectiveness of RAD-Seq as a tool for rapid and efficient generation of QTL-targeted and genome-wide marker data in a large complex genome, and its possible utility in farmed animal selection programs.

Transcriptomic analysis of the host response to early stage salmonid alphavirus (SAV-1) infection in Atlantic salmon Salmo salar L

Salmon pancreas disease, caused by salmonid alphavirus (SAV) of the family Togaviridae, is an economically important disease affecting farmed Atlantic salmon (Salmo salar L.) in Scotland, Norway, and Ireland. The virus causes characteristic lesions in the pancreas, heart, kidney and skeletal muscle of infected fish. The mechanisms responsible for the pathology and the immune responses elicited in infected Atlantic salmon are not fully understood. A microarray-based study was therefore performed to evaluate the host transcriptomic response during the early stages of an experimentally-induced SAV-1 infection. Atlantic salmon parr were injected intra-peritoneally with viral cell culture supernatant or cell culture supernatant without virus. RNA, extracted from head kidney sampled from infected and control fish at 1, 3 and 5 days post-injection (d.p.i.), was interrogated with the 17 k TRAITS/SGP cDNA microarray. The greatest number of significantly differentially expressed genes was recorded at 3 d.p.i., mainly associated with immune and defence mechanisms, including genes involved in interferon I pathways and Major Histocompatibility Complex Class I and II responses. Genes associated with apoptosis and cellular stress were also found to be differentially expressed between infected and uninfected individuals, as were genes involved in inhibiting viral attachment and replication. The microarray results were validated by follow-on analysis of eight genes by real-time PCR. The findings of the study reflect mechanisms used by the host to protect itself during the early stages of SAV-1 infection. In particular, there was evidence of rapid induction of interferon-mediated responses similar to those seen during mammalian alphavirus infections, and also early involvement of an adaptive immune response. This study provides essential knowledge to assist in the development of effective control and management strategies for SAV-1 infection.

Genetic resistance to rhabdovirus infection in teleost fish is paralleled to the derived cell resistance status

Genetic factors of resistance and predisposition to viral diseases explain a significant part of the clinical variability observed within host populations. Predisposition to viral diseases has been associated to MHC haplotypes and T cell immunity, but a growing repertoire of innate/intrinsic factors are implicated in the genetic determinism of the host susceptibility to viruses. In a long-term study of the genetics of host resistance to fish rhabdoviruses, we produced a collection of double-haploid rainbow trout clones showing a wide range of susceptibility to Viral Hemorrhagic Septicemia Virus (VHSV) waterborne infection. The susceptibility of fibroblastic cell lines derived from these clonal fish was fully consistent with the susceptibility of the parental fish clones. The mechanisms determining the host resistance therefore did not associate with specific host immunity, but rather with innate or intrinsic factors. One cell line was resistant to rhabdovirus infection due to the combination of an early interferon IFN induction - that was not observed in the susceptible cells - and of yet unknown factors that hamper the first steps of the viral cycle. The implication of IFN was well consistent with the wide range of resistance of this genetic background to VSHV and IHNV, to the birnavirus IPNV and the orthomyxovirus ISAV. Another cell line was even more refractory to the VHSV infection through different antiviral mechanisms. This collection of clonal fish and isogenic cell lines provides an interesting model to analyze the relative contribution of antiviral pathways to the resistance to different viruses.

IPNV with high and low virulence: host immune responses and viral mutations during infection

Background

Infectious pancreatic necrosis virus (IPNV) is an aquatic member of the Birnaviridae family that causes widespread disease in salmonids. IPNV is represented by multiple strains with markedly different virulence. Comparison of isolates reveals hyper variable regions (HVR), which are presumably associated with pathogenicity. However little is known about the rates and modes of sequence divergence and molecular mechanisms that determine virulence. Also how the host response may influence IPNV virulence is poorly described.

Methods

In this study we compared two field isolates of IPNV (NFH-Ar and NFH-El). The sequence changes, replication and mortality were assessed following experimental challenge of Atlantic salmon. Gene expression analyses with qPCR and microarray were applied to examine the immune responses in head kidney.

Results

Significant differences in mortality were observed between the two isolates, and viral load in the pancreas at 13 days post infection (d p.i.) was more than 4 orders of magnitude greater for NFH-Ar in comparison with NFH-El. Sequence comparison of five viral genes from the IPNV isolates revealed different mutation rates and Ka/Ks ratios. A strong tendency towards non-synonymous mutations was found in the HRV of VP2 and in VP3. All mutations in VP5 produced precocious stop codons. Prior to the challenge, NFH-Ar and NFH-El possessed high and low virulence motifs in VP2, respectively. Nucleotide substitutions were noticed already during passage of viruses in CHSE-214 cells and their accumulation continued in the challenged fish. The sequence changes were notably directed towards low virulence. Co-ordinated activation of anti-viral genes with diverse functions (IFN-a1 and c, sensors - Rig-I, MDA-5, TLR8 and 9, signal transducers - Srk2, MyD88, effectors - Mx, galectin 9, galectin binding protein, antigen presentation - b2-microglobulin) was observed at 13 d p.i. (NFH-Ar) and 29 d p.i. (both isolates).

Conclusions

Mortality and expression levels of the immune genes were directly related to the rate of viral replication, which was in turn associated with sequences of viral genes. Rapid changes in the viral genome that dramatically reduced virus proliferation might indicate a higher susceptibility to protective mechanism employed by the host. Disease outbreak and mortality depend on a delicate balance between host defence, regulation of signalling cascades and virus genomic properties.