Resistance to amoebic gill disease (AGD) is characterised by the transcriptional dysregulation of immune and cell cycle pathways

Resistance to Saprolegnia parasitica infection: A heritable trait in Atlantic salmon

A controlled Saprolegnia parasitica infection model was used to challenge 1158 fish representing 105 pedigreed Atlantic salmon families to evaluate the possibility of selecting for Saprolegnia resistance in a commercial breeding programme. Fish were infected in five study tanks and observed for 40 days post-infection for lesion score and survival. Survival analysis of the top 10 resistant and bottom 10 susceptible families indicated that the hazard of dying following Saprolegnia infection was 1509% higher in susceptible families. In all fish, a 10 g increase in weight correlated with a 7.8% increase in the hazard of dying while sex did not affect mortality. Resistance to Saprolegnia was estimated to have a heritability of 0.25, indicating that selection is possible. Genetic and phenotypic correlations indicated that the 11-point scoring system, developed in this study to quantify Saprolegnia infection severity, had a high negative correlation with survival as days to mortality at ≥-0.922(±0.005), suggesting that the scoring method could help assess lesion development in studies where mortality is not the primary biological endpoint.

Diseases Caused by Amoebae in Fish: An Overview

Parasitic and amphizoic amoebae are ubiquitous and can affect a huge variety of hosts, from invertebrates to humans, and fish are not an exception. Most of the relationships between amoebae and fish are based on four different types: ectocommensals, ectoparasites, endocommensals and endoparasites, although the lines between them are not always clear. As ectocommensals, they are located specially on the gills and particularly the amphizoic Neoparamoeba perurans is the most relevant species, being a real pathogenic parasite in farmed salmon. It causes amoebic gill disease, which causes a progressive hyperplasia of epithelial cells in the gill filaments and lamellae. Nodular gill disease is its analogue in freshwater fish but the causative agent is still not clear, although several amoebae have been identified associated to the lesions. Other species have been described in different fish species, affecting not only gills but also other organs, even internal ones. In some cases, species of the genera Naegleria or Acanthamoeba, which also contain pathogenic species affecting humans, are usually described affecting freshwater fish species. As endocommensals, Entamoebae species have been described in the digestive tract of freshwater and marine fish species, but Endolimax nana can reach other organs and cause systemic infections in farmed Solea senegalensis. Other systemic infections caused by amoebae are usually described in wild fish, although in most cases these are isolated cases without clinical signs or significance.

Bacteriomic Profiling of Branchial Lesions Induced by Neoparamoeba perurans Challenge Reveals Commensal Dysbiosis and an Association with Tenacibaculum dicentrarchi in AGD-Affected Atlantic Salmon (Salmo salar L.)

Amoebic gill disease is a parasitic condition that commonly affects marine farmed Atlantic salmon. The causative agent, Neoparamoeba perurans, induces a marked proliferation of the gill mucosa and focal superficial necrosis upon branchial lesions. The effect that amoebic branchialitis has upon gill associated commensal bacteria is unknown. A 16S rRNA sequencing approach was employed to profile changes in bacterial community composition, within amoebic gill disease (AGD)-affected and non-affected gill tissue. The bacterial diversity of biopsies with and without diseased tissue was significantly lower in the AGD-affected fish compared to uninfected fish. Furthermore, within the AGD-affected tissue, lesions appeared to contain a significantly higher abundance of the Flavobacterium, Tenacibaculum dicentrarchi compared to adjunct unaffected tissues. Quantitative PCR specific to both N. perurans and T. dicentrarchi was used to further examine the co-abundance of these known fish pathogens. A moderate positive correlation between these pathogens was observed. Taken together, the present study sheds new light on the complex interaction between the host, parasite and bacterial communities during AGD progression. The role that T. dicentrarchi may play in this complex relationship requires further investigation.

Integration of Transcriptome, Gross Morphology and Histopathology in the Gill of Sea Farmed Atlantic Salmon (Salmo salar): Lessons From Multi-Site Sampling

The gill of teleost fish is a multifunctional organ involved in many physiological processes such as gas exchange, osmotic and ionic regulation, acid-base balance and excretion of nitrogenous waste. Due to its extensive interface with the environment, the gill plays a key role as a primary mucosal defense tissue against pathogens, as manifested by the presence of the gill-associated lymphoid tissue (GIALT). In recent years, the prevalence of multifactorial gill pathologies has increased significantly, causing substantial losses in Atlantic salmon aquaculture. The transition from healthy to unhealthy gill phenotypes and the progression of multifactorial gill pathologies, such as proliferative gill disease (PGD), proliferative gill inflammation (PGI) and complex gill disorder (CGD), are commonly characterized by epithelial hyperplasia, lamellar fusion and inflammation. Routine monitoring for PGD relies on visual inspection and non-invasive scoring of the gill tissue (gross morphology), coupled with histopathological examination of gill sections. To explore the underlying molecular events that are associated with the progression of PGD, we sampled Atlantic salmon from three different marine production sites in Scotland and examined the gill tissue at three different levels of organization: gross morphology with the use of PGD scores (macroscopic examination), whole transcriptome (gene expression by RNA-seq) and histopathology (microscopic examination). Our results strongly suggested that the changes in PGD scores of the gill tissue were not associated with the changes in gene expression or histopathology. In contrast, integration of the gill RNA-seq data with the gill histopathology enabled us to identify common gene expression patterns associated with multifactorial gill disease, independently from the origin of samples. We demonstrated that the gene expression patterns associated with multifactorial gill disease were dominated by two processes: a range of immune responses driven by pro-inflammatory cytokines and the events associated with tissue damage and repair, driven by caspases and angiogenin.

Genome-wide association mapping and accuracy of predictions for amoebic gill disease in Atlantic salmon (Salmo salar)

Amoebic gill disease (AGD) is a parasitic disease caused by the amoeba Paramoeba perurans, which colonizes the gill tissues and causes distress for the host. AGD can cause high morbidity and mortalities in salmonid and non-salmonid fish species. To understand the genetic basis of AGD and improve health status of farmed A. salmon, a population of ~ 6,100 individuals belonging to 150 full-sib families was monitored for development of AGD in the sea of Ireland. The population was followed for two rounds of AGD infections, and fish were gill scored to identify severity of disease in first (N = 3,663) and the second (N = 3,511) infection with freshwater treatment after the first gill-scoring. A subset of this gill-scored population (N = 1,141) from 119 full-sib families were genotyped with 57,184 SNPs using custom-made Affymetrix SNP-chip. GWAS analyses were performed which resulted in five significantly associated SNP variants distributed over chromosome 1, 2 and 5. Three candidate genes; c4, tnxb and slc44a4 were found within QTL region of chromosome 2. The tnxb and c4 genes are known to be a part of innate immune system, and may play a role in resistance to AGD. The gain in prediction accuracy obtained by involving genomic information was 9–17% higher than using traditional pedigree information.

Characterising the mechanisms underlying genetic resistance to amoebic gill disease in Atlantic salmon using RNA sequencing

Background

Gill health is one of the main concerns for Atlantic salmon aquaculture, and Amoebic Gill Disease (AGD), attributable to infection by the amoeba Neoparamoeba perurans, is a frequent cause of morbidity. In the absence of preventive measures, increasing genetic resistance of salmon to AGD via selective breeding can reduce the incidence of the disease and mitigate gill damage. Understanding the mechanisms leading to AGD resistance and the underlying causative genomic features can aid in this effort, while also providing critical information for the development of other control strategies. AGD resistance is considered to be moderately heritable, and several putative QTL have been identified. The aim of the current study was to improve understanding of the mechanisms underlying AGD resistance, and to identify putative causative genomic factors underlying the QTL. To achieve this, RNA was extracted from the gill and head kidney of AGD resistant and susceptible animals following a challenge with N. perurans, and sequenced.

Results

Comparison between resistant and susceptible animals primarily highlighted differences mainly in the local immune response in the gill, involving red blood cell genes and genes related to immune function and cell adhesion. Differentially expressed immune genes pointed to a contrast in Th2 and Th17 responses, which is consistent with the increased heritability observed after successive challenges with the amoeba. Five QTL-region candidate genes showed differential expression, including a gene connected to interferon responses (GVINP1), a gene involved in systemic inflammation (MAP4K4), and a positive regulator of apoptosis (TRIM39). Analyses of allele-specific expression highlighted a gene in the QTL region on chromosome 17, cellular repressor of E1A-stimulated genes 1 (CREG1), showing allelic differential expression suggestive of a cis-acting regulatory variant.

Conclusions

In summary, this study provides new insights into the mechanisms of resistance to AGD in Atlantic salmon, and highlights candidate genes for further functional studies that can further elucidate the genomic mechanisms leading to resistance and contribute to enhancing salmon health via improved genomic selection.

In vitro gill cell monolayer successfully reproduces in vivo Atlantic salmon host responses to Neoparamoeba perurans infection

An in vitro model to study the host response to Neoparamoeba perurans, the causative agent of amoebic gill disease (AGD), was evaluated. The rainbow trout gill derived cell line, RTgill-W1, was seeded onto permeable cell culture supports and maintained asymmetrically with apical seawater. Cells were inoculated with either a passage attenuated or a recent wild clone of N. perurans. Amoebae, loaded with phagocytosed fluorescent beads, were observed associated with host cells within 20 min post inoculation (pi). By 6 h small foci of cytopathic effect appeared and at 72 h cytolysis was observed, with total disruption of the cell monolayer at 96 h pi. Due to cell monolayer disruption, the platform could not support proliferation of amoebae, which showed a 3-log reduction in parasite 18S rRNA mRNA after 72 h (10⁶ copies at 1 h to 10³ at 72 h pi). SEM observations showed amoebae-like cells with either short pseudopodia and a malleiform shape, or, long pseudopodia embedded within the gill cells and erosion of the cell monolayer. To study the host immune response, inoculated gill cells were harvested from triplicate inserts at 0, 1, 3, 6, 24 and 48 h pi, and expression of 12 genes involved in the Atlantic salmon response to AGD was compared between infected and uninfected cells and between amoebic clones. Both clones induced similar host inmate immune responses, with the up-regulation of proinflammatory cytokine IL1β, complement C3 and cell receptor MHC-1. The Th2 pathway was up-regulated, with increased gene expression of the transcription factor GATA3, and Th2 cytokines IL10, IL6 and IL4/13A. PCNA and AG-2 were also up-regulated. The wild clone induced significantly higher up-regulation of IL1β, MHC-1, PCNA, lysozyme and IL10 than the attenuated clone for at least some exposure times, but AG-2 gene expression was higher in cells inoculated with the attenuated one. A principal component analysis showed that AG-2 and IL10 were key genes in the in vitro host response to N. perurans. This in vitro model has proved to be a promising tool to study host responses to amoebae and may therefore reduce the requirement for in vivo studies when evaluating alternative therapeutants to AGD control.

Genomic and Transcriptomic Analysis of Amoebic Gill Disease Resistance in Atlantic Salmon (Salmo salar L.)

Amoebic gill disease (AGD) is one of the most important parasitic diseases of farmed Atlantic salmon. It is a source of major economic loss to the industry and poses significant threats to animal welfare. Previous studies have shown that resistance against this disease has a moderate, heritable genetic component, although the genes and the genetic pathways that contribute to this process have yet to be elucidated. In this study, to identify the genetic mechanisms of AGD resistance, we first investigated the molecular signatures of AGD infection in Atlantic salmon through a challenge model, where we compared the transcriptome profiles of the naïve and infected animals. We then conducted a genome-wide association analysis with 1,333 challenged tested fish to map the AGD resistance genomic regions, supported by the results from the transcriptomic data. Further, we investigated the potential of incorporating gene expression analysis results in genomic prediction to improve prediction accuracy. Our data suggest thousands of genes have modified their expression following infection, with a significant increase in the transcription of genes with functional properties in cell adhesion and a sharp decline in the abundance of various components of the immune system genes. From the genome-wide association analysis, QTL regions on chromosomes ssa04, ssa09, and ssa13 were detected to be linked with AGD resistance. In particular, we found that QTL region on ssa04 harbors members of the cadherin gene family. These genes play a critical role in target recognition and cell adhesion. The QTL region on ssa09 also is associated with another member of the cadherin gene family, protocadherin Fat 4. The associated genetic markers on ssa13 span a large genomic region that includes interleukin-18-binding protein, a gene with function essential in inhibiting the proinflammatory effect of cytokine IL18. Incorporating gene expression information through a weighted genomic relationship matrix approach decreased genomic prediction accuracy and increased bias of prediction. Together, these findings help to improve our breeding programs and animal welfare against AGD and advance our knowledge of the genetic basis of host-pathogen interactions.

Omics in fish mucosal immunity

The mucosal immune system of fish is a complex network of immune cells and molecules that are constantly surveilling the environment and protecting the host from infection. A number of "omics" tools are now available and utilized to understand the complexity of mucosal immune systems in non-traditional animal models. This review summarizes recent advances in the implementation of "omics" tools pertaining to the four mucosa-associated lymphoid tissues in teleosts. Genomics, transcriptomics, proteomics, and "omics" in microbiome research require interdisciplinary collaboration and careful experimental design. The data-rich datasets generated are proving really useful at discovering new innate immune players in fish mucosal secretions, identifying novel markers of specific mucosal immune responses, unraveling the diversity of the B and T cell repertoires and characterizing the diversity of the microbial communities present in teleost mucosal surfaces. Bioinformatics, data analysis and storage platforms should be developed to facilitate rapid processing of large datasets, especially when mammalian tools such as bioinformatics analysis software are not available in fishes.

Infectious disease, shifting climates, and opportunistic predators: cumulative factors potentially impacting wild salmon declines

Emerging diseases are impacting animals under high-density culture, yet few studies assess their importance to wild populations. Microparasites selected for enhanced virulence in culture settings should be less successful maintaining infectivity in wild populations, as once the host dies, there are limited opportunities to infect new individuals. Instead, moderately virulent microparasites persisting for long periods across multiple environments are of greatest concern. Evolved resistance to endemic microparasites may reduce susceptibilities, but as barriers to microparasite distributions are weakened, and environments become more stressful, unexposed populations may be impacted and pathogenicity enhanced. We provide an overview of the evolutionary and ecological impacts of infectious diseases in wild salmon and suggest ways in which modern technologies can elucidate the microparasites of greatest potential import. We present four case studies that resolve microparasite impacts on adult salmon migration success, impact of river warming on microparasite replication, and infection status on susceptibility to predation. Future health of wild salmon must be considered in a holistic context that includes the cumulative or synergistic impacts of multiple stressors. These approaches will identify populations at greatest risk, critically needed to manage and potentially ameliorate the shifts in current or future trajectories of wild populations.

Dietary vegetable oils do not alter the intestine transcriptome of gilthead sea bream (Sparus aurata), but modulate the transcriptomic response to infection with Enteromyxum leei

Background

Studies conducted with gilthead sea bream (Sparus aurata L.) have determined the maximum dietary replacement of fish meal and oil without compromising growth or product quality. The present study aimed to analyze the effect of the nutritional background on fish health and fish fed plant protein-based diets with fish oil (FO diet) or a blend of vegetable oils (66VO diet) were exposed for 102 days to the intestinal myxosporean parasite Enteromyxum leei, and the intestine transcriptome was analyzed with a customized oligo-microarray of 7,500 annotated genes.

Results

Infection prevalence was high and similar in the two diet groups, but the outcome of the disease was more pronounced in fish fed the 66VO diet. No differences were found in the transcriptome of both diet control groups, whereas the number of differentially expressed genes in infected groups was considerable. K-means clustering of these differentially expressed genes identified four expression patterns that reflected the progression of the disease with the magnitude of the fold-change being higher in infected 66VO fish. A positive correlation was found between the time of infection and the magnitude of the transcriptional change within the 66VO group, being higher in early infected animals. Within this diet group, a strong up-regulation of many components of the immune specific response was evidenced, whereas other genes related to complement response and xenobiotic metabolism were down-regulated.

Conclusions

The high replacement of fish oil by vegetable oils in practical fish feeds did not modify the intestine transcriptome of gilthead sea bream, but important changes were apparent when fish were exposed to the myxosporean E. leei. The detected changes were mostly a consequence rather than a cause of the different disease progression in the two diet groups. Hence, the developed microarray constitutes an excellent diagnostic tool to address changes associated with the action of intestinal pathogens, but lacks a prognostic value to predict in advance the different susceptibility of growing fish to the current pathogen.

Molecular profiling of the gilthead sea bream (Sparus aurata L.) response to chronic exposure to the myxosporean parasite Enteromyxum leei

The aim of the present work was to investigate the transcriptome response of gilthead sea bream (Sparus aurata) after challenge with the myxosporean Enteromyxum leei, a wide-spread enteric parasite causing heavy economic losses in Mediterranean sparid farms. This parasite causes severe desquamative enteritis which usually leads to death of the fish, and there are no preventative or curative treatments for this enteromyxosis. After 113 days of exposure to parasite-contaminated effluent, fish were classified into three cohorts: control fish not exposed to parasite, those that were exposed and infected, and those that were exposed but not infected. In order to detect target genes that may be candidates for infective status or resistance, a cDNA microarray containing 18,490 cDNA clones enriched in genes differentially expressed after infection was hybridised with head kidney and intestine samples. In infected fish, 371 and 373 genes were differentially regulated at the >1.5-fold level in intestine and head kidney respectively, whereas in non-infected fish 175 and 501 genes were differentially regulated in these tissues, respectively. A global marked gene down-regulation was evident in infected fish, mainly in genes involved in the immune and acute phase response particularly complement and mannose binding lectin. Microarray analysis demonstrated a complex interplay between host and/or parasite derived proteases and protease inhibitors, apoptosis, cell proliferation and antioxidant defence genes in exposed fish. In the head kidney of non-infected fish a marked depression of genes involved in the acute phase response was evident. By contrast, in the intestine of non-infected fish, interferon-stimulated and MHC class II genes involved in antigen processing and presentation were up-regulated, possibly indicating that an active immune response at the local level is important to avoid infection with or proliferation of the parasite.

Starvation alters the liver transcriptome of the innate immune response in Atlantic salmon (Salmo salar)

BACKGROUND

The immune response is an energy demanding process, which has effects in many physiological pathways in the body including protein and lipid metabolism. During an inflammatory response the liver is required to produce high levels of acute phase response proteins that attempt to neutralise an invading pathogen. Although this has been extensively studied in both mammals and fish, little is known about how high and low energy reserves modulate the response to an infection in fish which are ectothermic vertebrates. Food withdrawal in fish causes a decrease in metabolic rate so as to preserve protein and lipid energy reserves, which occurs naturally during the life cycle of many salmonids. Here we investigated how the feeding or fasting of Atlantic salmon affected the transcriptional response in the liver to an acute bacterial infection.

RESULTS

Total liver RNA was extracted from four different groups of salmon. Two groups were fed or starved for 28 days. One of each of the fed or starved groups was then exposed to an acute bacterial infection. Twenty four hours later (day 29) the livers were isolated from all fish for RNA extraction. The transcriptional changes were examined by micro array analysis using a 17 K Atlantic salmon cDNA microarray. The expression profiling results showed major changes in gene transcription in each of the groups. Enrichment for particular biological pathways was examined by analysis of gene ontology. Those fish that were starved decreased immune gene transcription and reduced production of plasma protein genes, and upon infection there was a further decrease in genes encoding plasma proteins but a large increase in acute phase response proteins. The latter was greater in magnitude than in the fish that had been fed prior to infection. The expression of several genes that were found altered during microarray analysis was confirmed by real time PCR.

CONCLUSIONS

We demonstrate that both starvation and infection have profound effects on transcription in the liver of salmon. There was a significant effect on the transcriptional response to infection depending on the prior feeding regime of the fish. It is likely that the energy demands on protein synthesis for acute phase response proteins are relatively high in the starved fish which have reduced energy reserves. This has implications for dietary control of fish if an immune response is anticipated.

Differential gene expression during smoltification of Atlantic salmon (Salmo salar L.): a first large-scale microarray study

The life cycle of the Atlantic salmon (Salmo salar) involves a period of 1 to 3 years in freshwater followed by migration to the sea where the salmon undergoes rapid growth. In preparation for the marine environment, while still in freshwater, the salmon undergo a transformation from a freshwater dwelling parr to a saltwater adapted smolt, a process known as smoltification. The Atlantic salmon Transcriptome Analysis of Important Traits of Salmon/Salmon Genome Project (TRAITS/SGP) cDNA microarray was used to investigate how gene expression alters during smoltification. Genes differentially expressed during smoltification were identified by comparing gene expression profiles in smolt brain, gill, and kidney tissue samples with those of parr. Of the three tissues investigated, the number of differentially expressed genes was the greatest in gill. Many of the differentially expressed genes could be assigned to one of four main categories: growth, metabolism, oxygen transport, and osmoregulation. Quantitative polymerase chain reaction successfully confirmed the differential expression of seven of the upregulated genes. The TRAITS/SGP cDNA microarray was used to successfully demonstrate for the first time how gene expression mediates smoltification in the Atlantic salmon. Changes in gene expression observed in this study reflected the physiological and biochemical changes recorded by previous studies describing the parr-smolt transformation. This study significantly increases our knowledge of smoltification and will benefit future studies in this area of research.

Amoebic gill disease resistance is not related to the systemic antibody response of Atlantic salmon, Salmo salar L

Amoebic gill disease (AGD) is a proliferative gill tissue response caused by Neoparamoeba perurans and is the main disease affecting Australian marine farmed Atlantic salmon. We have previously proposed that macroscopic gill health ('gill score') trajectories and challenge survival provide evidence of a change in the nature of resistance to AGD. In order to examine whether the apparent development of resistance was because of an adaptive response, serum was sequentially sampled from the same individuals over the first three rounds of natural AGD infection and from survivors of a subsequent non-intervention AGD survival challenge. The systemic immune reaction to 'wildtype'Neoparamoeba sp. was characterized by Western blot analysis and differentiated to putative carbohydrate or peptide epitopes by periodate oxidation reactions. The proportion of seropositive fish increased from 46% to 77% with each AGD round. Antibody response to carbohydrate epitope(s) was immunodominant, occurring in 43-64% of samples. Antibodies that bound peptide epitope were identified in 16% of the challenge survivors. A 1:50 (single-dilution) enzyme-linked immunosorbent assay confirmed a measurable immune titre in 13% of the survivors. There was no evidence that antibodies recognizing wildtype Neoparamoeba provided significant protection against AGD.