Genome-Wide Association Analysis for Resistance to Infectious Pancreatic Necrosis Virus Identifies Candidate Genes Involved in Viral Replication and Immune Response in Rainbow Trout (Oncorhynchus mykiss)

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.

G2P Provides an Integrative Environment for Multi-model genomic selection analysis to improve genotype-to-phenotype prediction

Genotype-to-phenotype (G2P) prediction has become a mainstream paradigm to facilitate genomic selection (GS)-assisted breeding in the seed industry. Many methods have been introduced for building GS models, but their prediction precision may vary depending on species and specific traits. Therefore, evaluation of multiple models and selection of the appropriate one is crucial to effective GS analysis. Here, we present the G2P container developed for the Singularity platform, which not only contains a library of 16 state-of-the-art GS models and 13 evaluation metrics. G2P works as an integrative environment offering comprehensive, unbiased evaluation analyses of the 16 GS models, which may be run in parallel on high-performance computing clusters. Based on the evaluation outcome, G2P performs auto-ensemble algorithms that not only can automatically select the most precise models but also can integrate prediction results from multiple models. This functionality should further improve the precision of G2P prediction. Another noteworthy function is the refinement design of the training set, in which G2P optimizes the training set based on the genetic diversity analysis of a studied population. Although the training samples in the optimized set are fewer than in the original set, the prediction precision is almost equivalent to that obtained when using the whole set. This functionality is quite useful in practice, as it reduces the cost of phenotyping when constructing training population. The G2P container and source codes are freely accessible at https://g2p-env.github.io/.

Assessment of listing and categorisation of animal diseases within the framework of the Animal Health Law (Regulation (EU) No 2016/429): infectious pancreatic necrosis (IPN)

Infectious pancreatic necrosis (IPN) was assessed according to the criteria of the Animal Health Law (AHL), in particular, the criteria of Article 7 on disease profile and impacts, Article 5 on its eligibility to be listed, Annex IV for its categorisation according to disease prevention and control rules as in Article 9, and Article 8 for listing animal species related to IPN. The assessment was performed following a methodology previously published. The outcome reported is the median of the probability ranges provided by the experts, which indicates whether each criterion is fulfilled (lower bound ≥ 66%) or not (upper bound ≤ 33%), or whether there is uncertainty about fulfilment. Reasoning points are reported for criteria with an uncertain outcome. According to the assessment here performed, it is uncertain whether IPN can be considered eligible to be listed for Union intervention according to Article 5 of the AHL (50-90% probability). According to the criteria in Annex IV, for the purpose of categorisation related to the level of prevention and control as in Article 9 of the AHL, the AHAW Panel concluded that IPN does not meet the criteria in Section 1 (Category A; 0-1% probability of meeting the criteria) and it is uncertain whether it meets the criteria in Sections 2, 3, 4 and 5 (Categories B, C, D and E; 33-66%, 33-66%, 50-90% and 50-99% probability of meeting the criteria, respectively). The animal species to be listed for IPN according to Article 8 criteria are provided.

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.

GWAS on Imputed Whole-Genome Sequence Variants Reveal Genes Associated with Resistance to Piscirickettsia salmonis in Rainbow Trout (Oncorhynchus mykiss)

Genome-wide association studies (GWAS) allow the identification of associations between genetic variants and important phenotypes in domestic animals, including disease-resistance traits. Whole Genome Sequencing (WGS) data can help increase the resolution and statistical power of association mapping. Here, we conduced GWAS to asses he facultative intracellular bacterium Piscirickettsia salmonis, which affects farmed rainbow trout, Oncorhynchus mykiss, in Chile using imputed genotypes at the sequence level and searched for candidate genes located in genomic regions associated with the trait. A total of 2130 rainbow trout were intraperitoneally challenged with P. salmonis under controlled conditions and genotyped using a 57K single nucleotide polymorphism (SNP) panel. Genotype imputation was performed in all the genotyped animals using WGS data from 102 individuals. A total of 488,979 imputed WGS variants were available in the 2130 individuals after quality control. GWAS revealed genome-wide significant quantitative trait loci (QTL) in Omy02, Omy03, Omy25, Omy26 and Omy27 for time to death and in Omy26 for binary survival. Twenty-four (24) candidate genes associated with P. salmonis resistance were identified, which were mainly related to phagocytosis, innate immune response, inflammation, oxidative response, lipid metabolism and apoptotic process. Our results provide further knowledge on the genetic variants and genes associated with resistance to intracellular bacterial infection in rainbow trout.

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.

Genome-wide association study of VHSV-resistance trait in Paralichthys olivaceus

In flounder aquaculture, selective breeding plays a vital role in the development of disease-resistant traits and animals with high growth rates. Moreover, superior animals are required to achieve high profits. Unlike growth-related traits, disease-resistant experiments need to be conducted in a controlled environment, as the improper measurement of traits often leads to low genetic correlation and incorrect estimation of breeding values. In this study, viral hemorrhagic septicemia virus (VHSV) resistance was studied using a genome-wide association study (GWAS), and the genetic parameters were estimated. Genotyping was performed using a high-quality 70 K single nucleotide polymorphism (SNP) Affymetrix® Axiom® myDesign™ Genotyping Array of olive flounder. A heritability of ∼0.18 for resistance to VHSV was estimated using genomic information of the fish. According to the GWAS, significant SNPs were detected in chromosomes 21, 24, and contig AGQT02032065.1. Three SNPs showed significance at the genome-wide level (p < 1 × 10^-6), while others showed significance above the suggestive cutoff (p < 1 × 10^-4). The 3% phenotypic variation was explained by the highest significant SNP, named AX-419319631. Of the important genes for disease resistance, SNPs were associated with plcg1, epha4, clstn2, pik3cb, hes6, meis3, prx6, cep164, siae, and kirrel3b. Most of the genes associated with these SNPs have been previously reported with respect to viral entry, propagation, and immune mechanisms. Therefore, our study provides helpful information regarding VHSV resistance in olive flounder, which can be used for breeding applications.

Differential Transcriptomic Response of Rainbow Trout to Infection with Two Strains of IPNV

The IPN virus (IPNV) causes a highly contagious disease that affects farmed salmonids. IPNV isolates have been phylogenetically classified into seven genogroups, of which two are present in Chile, genogroups 1 and 5. This study aimed to compare the transcriptomic response of rainbow trout fry challenged with two Chilean isolates of IPNV, RTTX (genogroup 1), and ALKA (genogroup 5). Tissue samples from challenged individuals and controls were taken at 1, 7, and 20 days post-challenge and analyzed by RNA-Seq. The results revealed that infection with RTTX elicited a greater modulation of the trout transcriptome compared to ALKA infection, generating a greater number of highly differentially expressed genes in relation to the control fish. Gene Ontology enrichment indicated that functions related to the inflammatory and immune responses were modulated in fish challenged with both isolates throughout the trial, but with different regulation patterns. On day 1 post challenge, these functions were activated in those challenged with ALKA, but suppressed in RTTX-challenged fish. These results suggest that rainbow trout exhibit a differential transcriptomic response to infection with the two genetically distinct IPNV isolates, especially at early times post-infection.

The nedd-8 activating enzyme gene underlies genetic resistance to infectious pancreatic necrosis virus in Atlantic salmon

Genetic resistance to infectious pancreatic necrosis virus (IPNV) in Atlantic salmon is a rare example of a trait where a single locus (QTL) explains almost all of the genetic variation. Genetic marker tests based on this QTL on salmon chromosome 26 have been widely applied in selective breeding to markedly reduce the incidence of the disease. In the current study, whole genome sequencing and functional annotation approaches were applied to characterise genes and variants in the QTL region. This was complemented by an analysis of differential expression between salmon fry of homozygous resistant and homozygous susceptible genotypes challenged with IPNV. These analyses pointed to the NEDD-8 activating enzyme 1 (nae1) gene as a putative functional candidate underlying the QTL effect. The role of nae1 in IPN resistance was further assessed via CRISPR-Cas9 knockout of the nae1 gene and chemical inhibition of the nae1 protein activity in Atlantic salmon cell lines, both of which resulted in highly significant reduction in productive IPNV replication. In contrast, CRISPR-Cas9 knockout of a candidate gene previously purported to be a cellular receptor for the virus (cdh1) did not have a major impact on productive IPNV replication. These results suggest that nae1 is the causative gene underlying the major QTL affecting resistance to IPNV in salmon, provide further evidence for the critical role of neddylation in host-pathogen interactions, and highlight the value in combining high-throughput genomics approaches with targeted genome editing to understand the genetic basis of disease resistance.

Whole-genome association study searching for QTL for Aeromonas salmonicida resistance in rainbow trout

Aeromonas salmonicida subsp. salmonicida, the causative agent of furunculosis, has extensive negative effects on wild and farmed salmonids worldwide. Vaccination induces some protection under certain conditions but disease outbreaks occur even in vaccinated fish. Therefore, alternative disease control approaches are required to ensure the sustainable expansion of rainbow trout aquaculture. Selective breeding can be applied to enhance host resistance to pathogens. The present work used genome-wide association study (GWAS) to identify quantitative trait loci (QTL) associated with A. salmonicida resistance in rainbow trout. A total 798 rainbow trout exposed to A. salmonicida by bath challenge revealed 614 susceptible and 138 resistant fish. Genotyping was conducted using the 57 K single nucleotide polymorphism (SNP) array and the GWAS was performed for survival and time to death phenotypes. We identified a QTL on chromosome 16 and located positional candidate genes in the proximity of the most significant SNPs. In addition, samples from exposed fish were examined for expression of 24 immune-relevant genes indicating a systematic immune response to the infection. The present work demonstrated that resistance to A. salmonicida is moderately heritable with oligogenic architecture. These result will be useful for the future breeding programs for improving the natural resistance of rainbow trout against furunculosis.

Systematic review of the status of veterinary epidemiological research in two species regarding the FAIR guiding principles

Background

The FAIR (Findable, Accessible, Interoperable, Reusable) principles were proposed in 2016 to set a path towards reusability of research datasets. In this systematic review, we assessed the FAIRness of datasets associated with peer-reviewed articles in veterinary epidemiology research published since 2017, specifically looking at salmonids and dairy cattle. We considered the differences in practices between molecular epidemiology, the branch of epidemiology using genetic sequences of pathogens and hosts to describe disease patterns, and non-molecular epidemiology.

Results

A total of 152 articles were included in the assessment. Consistent with previous assessments conducted in other disciplines, our results showed that most datasets used in non-molecular epidemiological studies were not available (i.e., neither findable nor accessible). Data availability was much higher for molecular epidemiology papers, in line with a strong repository base available to scientists in this discipline. The available data objects generally scored favourably for Findable, Accessible and Reusable indicators, but Interoperability was more problematic.

Conclusions

None of the datasets assessed in this study met all the requirements set by the FAIR principles. Interoperability, in particular, requires specific skills in data management which may not yet be broadly available in the epidemiology community. In the discussion, we present recommendations on how veterinary research could move towards greater reusability according to FAIR principles. Overall, although many initiatives to improve data access have been started in the research community, their impact on the availability of datasets underlying published articles remains unclear to date.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12917-021-02971-1.

Increased accuracy of genomic predictions for growth under chronic thermal stress in rainbow trout by prioritizing variants from GWAS using imputed sequence data

Through imputation of genotypes, genome‐wide association study (GWAS) and genomic prediction (GP) using whole‐genome sequencing (WGS) data are cost‐efficient and feasible in aquaculture breeding schemes. The objective was to dissect the genetic architecture of growth traits under chronic heat stress in rainbow trout (Oncorhynchus mykiss) and to assess the accuracy of GP based on imputed WGS and different preselected single nucleotide polymorphism (SNP) arrays. A total of 192 and 764 fish challenged to a heat stress experiment for 62 days were genotyped using a customized 1 K and 26 K SNP panels, respectively, and then, genotype imputation was performed from a low‐density chip to WGS using 102 parents (36 males and 66 females) as the reference population. Imputed WGS data were used to perform GWAS and test GP accuracy under different preselected SNP scenarios. Heritability was estimated for body weight (BW), body length (BL) and average daily gain (ADG). Estimates using imputed WGS data ranged from 0.33 ± 0.05 to 0.55 ± 0.05 for growth traits under chronic heat stress. GWAS revealed that the top five cumulatively SNPs explained a maximum of 0.94%, 0.86% and 0.51% of genetic variance for BW, BL and ADG, respectively. Some important functional candidate genes associated with growth‐related traits were found among the most important SNPs, including signal transducer and activator of transcription 5B and 3 (STAT5B and STAT3, respectively) and cytokine‐inducible SH2‐containing protein (CISH). WGS data resulted in a slight increase in prediction accuracy compared with pedigree‐based method, whereas preselected SNPs based on the top GWAS hits improved prediction accuracies, with values ranging from 1.2 to 13.3%. Our results support the evidence of the polygenic nature of growth traits when measured under heat stress. The accuracies of GP can be improved using preselected variants from GWAS, and the use of WGS marginally increases prediction accuracy.

Transcriptomic Analysis of Inbred Chicken Lines Reveals Infectious Bursal Disease Severity Is Associated with Greater Bursal Inflammation In Vivo and More Rapid Induction of Pro-Inflammatory Responses in Primary Bursal Cells Stimulated Ex Vivo

In order to better understand differences in the outcome of infectious bursal disease virus (IBDV) infection, we inoculated a very virulent (vv) strain into White Leghorn chickens of inbred line W that was previously reported to experience over 24% flock mortality, and three inbred lines (15I, C.B4 and 0) that were previously reported to display no mortality. Within each experimental group, some individuals experienced more severe disease than others but line 15I birds experienced milder disease based on average clinical scores, percentage of birds with gross pathology, average bursal lesion scores and average peak bursal virus titre. RNA-Seq analysis revealed that more severe disease in line W was associated with significant up-regulation of pathways involved in inflammation, cytoskeletal regulation by Rho GTPases, nicotinic acetylcholine receptor signaling, and Wnt signaling in the bursa compared to line 15I. Primary bursal cell populations isolated from uninfected line W birds contained a significantly greater percentage of KUL01+ macrophages than cells isolated from line 15I birds (p < 0.01) and, when stimulated ex vivo with LPS, showed more rapid up-regulation of pro-inflammatory gene expression than those from line 15I birds. We hypothesize that a more rapid induction of pro-inflammatory cytokine responses in bursal cells following IBDV infection leads to more severe disease in line W birds than in line 15I.

Genomics to accelerate genetic improvement in tilapia

Selective breeding of tilapia populations started in the early 1990s and over the past three decades tilapia has become one of the most important farmed freshwater species, being produced in more than 125 countries around the globe. Although genome assemblies have been available since 2011, most of the tilapia industry still depends on classical selection techniques using mass spawning or pedigree information to select for growth traits with reported genetic gains of up to 20% per generation. The involvement of international breeding companies and research institutions has resulted in the rapid development and application of genomic resources in the last few years. GWAS and genomic selection are expected to contribute to uncovering the genetic variants involved in economically relevant traits and increasing the genetic gain in selective breeding programs, respectively. Developments over the next few years will probably focus on achieving a deep understanding of genetic architecture of complex traits, as well as accelerating genetic progress in the selection for growth-, quality- and robustness-related traits. Novel phenotyping technologies (i.e. phenomics), lower-cost whole-genome sequencing approaches, functional genomics and gene editing tools will be crucial in future developments for the improvement of tilapia aquaculture.

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.

High-Throughput Single Nucleotide Polymorphism (SNP) Discovery and Validation Through Whole-Genome Resequencing in Nile Tilapia (Oreochromis niloticus)

Nile tilapia (Oreochromis niloticus) is the second most important farmed fish in the world and a sustainable source of protein for human consumption. Several genetic improvement programs have been established for this species in the world. Currently, the estimation of genetic merit of breeders is typically based on genealogical and phenotypic information. Genome-wide information can be exploited to efficiently incorporate traits that are difficult to measure into the breeding goal. Thus, single nucleotide polymorphisms (SNPs) are required to investigate phenotype-genotype associations and determine the genomic basis of economically important traits. We performed de novo SNP discovery in three different populations of farmed Nile tilapia. A total of 29.9 million non-redundant SNPs were identified through Illumina (HiSeq 2500) whole-genome resequencing of 326 individual samples. After applying several filtering steps, including removing SNP based on genotype and site quality, presence of Mendelian errors, and non-unique position in the genome, a total of 50,000 high-quality SNPs were selected for the development of a custom Illumina BeadChip SNP panel. These SNPs were highly informative in the three populations analyzed showing between 43,869 (94%) and 46,139 (99%) SNPs in Hardy-Weinberg Equilibrium; 37,843 (76%) and 45,171(90%) SNPs with a minor allele frequency (MAF) higher than 0.05; and 43,450 (87%) and 46,570 (93%) SNPs with a MAF higher than 0.01. The 50K SNP panel developed in the current work will be useful for the dissection of economically relevant traits, enhancing breeding programs through genomic selection, as well as supporting genetic studies in farmed populations of Nile tilapia using dense genome-wide information.

Comparison of mortality and viral load in rainbow trout (Oncorhynchus mykiss) infected with infectious pancreatic necrosis virus (IPNV) genogroups 1 and 5

Infectious pancreatic necrosis virus (IPNV) is the aetiological agent of a highly contagious disease that affects farmed salmonids. IPNV isolates have been phylogenetically classified into eight genogroups, of which two are present in Chile, genogroups 1 and 5. Here, we compare the mortality rate caused by isolates from both genogroups in rainbow trout (Oncorhynchus mykiss) fry to determine if there is an association between host susceptibility and phylogenetic characterization of IPNV. Fish were challenged by immersion with one of four isolates (two for each genogroup), and mortality curves were assessed after 30 days. Viral load was measured in all mortalities and in live fish sampled at 1, 7 and 20 days post-infection. Although mortality was low throughout the challenge, differences were found between fish infected with different isolates. Both isolates from genogroup 1 caused greater cumulative mortalities than either of the isolates from genogroup 5. When combined, the overall mortality rate of fish challenged with genogroup 1 isolates was significantly higher than those infected with genogroup 5. However, viral load was lower on trout infected with genogroup 1 isolates. These results suggest that rainbow trout are more susceptible to IPNV isolates from genogroup 1 than genogroup 5.