A major effect quantitative trait locus for whirling disease resistance identified in rainbow trout (Oncorhynchus mykiss)

Mapping of a major QTL for increased robustness and detection of genome assembly errors in Asian seabass (Lates calcarifer)

BACKGROUND

For Asian seabass (Lates calcarifer, Bloch 1790) cultured at sea cages various aquatic pathogens, complex environmental and stress factors are considered as leading causes of disease, causing tens of millions of dollars of annual economic losses. Over the years, we conducted farm-based challenges by exposing Asian seabass juveniles to complex natural environmental conditions. In one of these challenges, we collected a total of 1,250 fish classified as either 'sensitive' or 'robust' individuals during the 28-day observation period.

RESULTS

We constructed a high-resolution linkage map with 3,089 SNPs for Asian seabass using the double digest Restriction-site Associated DNA (ddRAD) technology and a performed a search for Quantitative Trait Loci (QTL) associated with robustness. The search detected a major genome-wide significant QTL for increased robustness in pathogen-infected marine environment on linkage group 11 (ASB_LG11; 88.9 cM to 93.6 cM) with phenotypic variation explained of 81.0%. The QTL was positioned within a > 800 kb genomic region located at the tip of chromosome ASB_LG11 with two Single Nucleotide Polymorphism markers, R1-38468 and R1-61252, located near to the two ends of the QTL. When the R1-61252 marker was validated experimentally in a different mass cross population, it showed a statistically significant association with increased robustness. The majority of thirty-six potential candidate genes located within the QTL have known functions related to innate immunity, stress response or disease. By utilizing this ddRAD-based map, we detected five mis-assemblies corresponding to four chromosomes, namely ASB_LG8, ASB_LG9, ASB_LG15 and ASB_LG20, in the current Asian seabass reference genome assembly.

CONCLUSION

According to our knowledge, the QTL associated with increased robustness is the first such finding from a tropical fish species. Depending on further validation in other stocks and populations, it might be potentially useful for selecting robust Asian seabass lines in selection programs.

Control of parasitic diseases in aquaculture

Finfish aquaculture in freshwater and marine environments is continuously expanding globally, and the potential for a substantial further increase is well documented. The industry is supplying fish products for human consumption to the same extent as capture fisheries, and new fish species for domestication are still being selected by the industry. The challenge faced by all aquacultured species, classical and novel, is the range of pathogens associated with each new fish type. A fish host in its natural environment carries a series of more or less specific parasites (specialists and generalists). Some of these show a marked ability to propagate in aquaculture settings. They may then elicit disease when infection intensities in the confined aquaculture environment reach high levels. In addition, the risk of transmission of parasites from aquaculture enterprises to wild fish stocks adds to the parasitic challenge. Control programmes of various kinds are needed and these may include chemotherapeutants and medicines as the farmer's first and convenient choice, but mechanical, biological, immunological and genetic control methods are available solutions. New methods are still to be developed by scrutinizing the life cycle of each particular parasite species and pin-pointing the vulnerable stage to be targeted. As parasites exhibit a huge potential for adaptation to environmental changes, one must realize that only one approach rarely is sufficient. The present work therefore elaborates on and advocates for implementation of integrated control strategies for diseases caused by protozoan and metazoan parasites.

Genetic and epigenetic regulation of growth, reproduction, disease resistance and stress responses in aquaculture

Major progress has been made with genomic and genetic studies in aquaculture in the last decade. However, research on epigenetic regulation of aquaculture traits is still at an early stage. It is apparent that most, if not all, aquaculture traits are regulated at both genetic and epigenetic levels. This paper reviews recent progress in understanding of genetic and epigenetic regulation of important aquaculture traits such as growth, reproduction, disease resistance, and stress responses. Although it is challenging to make generalized statements, DNA methylation is mostly correlated with down-regulation of gene expression, especially when at promoters and enhancers. As such, methylation of growth factors and their receptors is negatively correlated with growth; hypomethylation of genes important for stress tolerance is correlated with increased stress tolerance; hypomethylation of genes important for male or female sex differentiation leads to sex differentiation into males or females, respectively. It is apparent that environmental regulation of aquaculture traits is mediated at the level of epigenetic regulation, and such environment-induced epigenetic changes appeared to be intergenerationally inherited, but evidences for transgenerational inheritance are still limited.

Responses towards eyefluke (Diplostomum pseudospathaceum) in different genetic lineages of rainbow trout

Marker-assisted selective breeding of fish with higher levels of resistance towards specific pathogens may improve fish health, but the impact of host genotype on susceptibility to multiple pathogen infections is still poorly investigated. This study examined the resistance in rainbow trout Oncorhynchus mykiss towards infection with the eye fluke Diplostomum pseudospathaceum. We used genetically selected rainbow trout, carrying SNPs associated with resistance towards the parasitic ciliate Ichthyophthirius multifiliis, and exposed the fish to eye fluke cercariae. We showed that fish partly resistant to I. multifiliis were more susceptible to eye fluke invasion. The expression of immune relevant genes (encoding innate and adaptive factors) was also affected as these genotypes responded less strongly to a secondary fluke infection. The complexity of genome architecture in disease resistance towards multiple pathogens is discussed.

Modulation of local and systemic immune responses in brown trout (Salmo trutta) following exposure to Myxobolus cerebralis

Myxobolus cerebralis, the etiological agent of Whirling Disease (WD), is a freshwater myxozoan parasite with considerable economic and ecological relevance for salmonids. There are differences in disease susceptibility between species and strains of salmonids. Recently, we have reported that the suppressor of cytokine signaling SOCS1 and SOCS3 are key in modulating rainbow trout (Oncorhynchus mykiss) immune responses and that resistant fish apparently exhibit effective Th17 cell response after exposure to M. cerebralis. It is unclear whether such molecules and pathways are also involved in the immune response of M. cerebralis infected brown trout (Salmo trutta). Hence, this study aimed to explore their role during immune modulation in infected brown trout, which is considered resistant to this parasite. Fish were exposed to the triactinomyxon (TAM) stages of M. cerebralis and quantitative real-time PCR (RT-qPCR) was carried out to examine local (caudal fin) and systemic (head kidney, spleen) immune transcriptional changes associated with WD over time in infected and control fish. All of the immune genes in the three tissues studied were differentially expressed in infected fish at multiple time points. Brown trout reduced the parasite load and demonstrated effective immune responses, likely by keeping pro-inflammatory and anti-inflammatory cytokines in balance whilst stimulating efficient Th17-mediated immunity. This study increases knowledge on the brown trout immune response to M. cerebralis and helps us to understand the underlying mechanisms of WD resistance.

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.

STAT3/SOCS3 axis contributes to the outcome of salmonid whirling disease

There are differences in disease susceptibility to whirling disease (WD) among strains of rainbow trout. The North American strain Trout Lodge (TL) is highly susceptible, whereas the German Hofer (HO) strain is more resistant. The suppressor of cytokine signaling (SOCS) proteins are key in inhibiting cytokine signaling. Their role in modulating the immune response against whirling disease is not completely clear. This study aimed at investigating the transcriptional response of SOCS1 and SOCS3 genes to Myxobolus cerebralis along with that of several upstream regulators and immune response genes. M. cerebralis induced the expression of SOCS1, the IL-6-dependent SOCS3, the anti-inflammatory cytokine IL-10 and the Treg associated transcription factor FOXP3 in TL fish at multiple time points, which likely caused a restricted STAT1 and STAT3 activity affecting the Th17/Treg17 balance. The expression of SOCS1 and the IL-6-dependent SOCS3 was induced constraining the activation of STAT1 and STAT3 in TL fish, thereby causing Th17/Treg17 imbalance and leaving the fish unable to establish a protective immune response against M. cerebralis or control inflammatory reactions increasing susceptibility to WD. Conversely, in HO fish, the expression of SOCS1 and SOCS3 was restrained, whereas the expression of STAT1 and IL-23-mediated STAT3 was induced potentially enabling more controlled immune responses, accelerating parasite clearance and elevating resistance. The induced expression of STAT1 and IL-23-mediated STAT3 likely maintained a successful Th17/Treg17 balance and enabled fish to promote effective immune responses favouring resistance against WD. The results provide insights into the role of SOCS1 and SOCS3 in regulating the activation and magnitude of host immunity in rainbow trout, which may help us understand the mechanisms that underlie the variation in resistance to WD.

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.

GWAS Analysis Indicated Importance of NF-κB Signaling Pathway in Host Resistance Against Motile Aeromonas Septicemia Disease in Catfish

Motile Aeromonas septicemia (MAS) disease caused by a bacterial pathogen, Aeromonas hydrophila, is an emerging but severe disease of catfish. Genetic enhancement of disease resistance is considered to be effective to control the disease. To provide an insight into the genomic basis of MAS disease resistance, in this study, we conducted a genome-wide association study (GWAS) to identify quantitative trait loci (QTL). A total of 1820 interspecific backcross catfish of 7 families were challenged with A. hydrophila, and 382 phenotypic extremes were selected for genotyping with the catfish 690 K SNP arrays. Three QTL on linkage group (LG) 2, 26 and 29 were identified to be significantly associated with MAS resistance. Within these regions, a total of 24 genes had known functions in immunity, 10 of which were involved in NF-κB signaling pathway, suggesting the importance of NF-κB signaling pathway in MAS resistance. In addition, three suggestively significant QTL were identified on LG 11, 17, and 20. The limited numbers of QTL involved in MAS resistance suggests that marker-assisted selection may be a viable approach for catfish breeding.

Kinetics of local and systemic immune cell responses in whirling disease infection and resistance in rainbow trout

Background

Whirling disease (WD), caused by the myxozoan parasite Myxobolus cerebralis, is responsible for high mortalities in rainbow trout hatcheries and natural populations. To elucidate how resistant and susceptible rainbow trout strains respond to early invasion, a well-established model of WD was used to demonstrate the kinetics of local and systemic immune responses in two rainbow trout strains, the susceptible American Trout Lodge (TL) and the more resistant German Hofer strain (HO).

Methods

Parasite load and cellular immune responses were compared across several time points after M. cerebralis exposure to elucidate the kinetics of immune cells in resistant and susceptible rainbow trout in response to early invasion. In the course of the 20 days following exposure, leukocyte kinetics was monitored by flow cytometry in the caudal fin (CF), head kidney (HK) and spleen (SP). For the analysis of the leukocyte composition, cells were stained using a set of monoclonal antibodies with known specificity for distinct subpopulations of rainbow trout leukocytes.

Results

Experiments indicated general increases of CF, HK and SP myeloid cells, while decreases of B cells and T cells in the SP and HK were observed at several time points in the TL strain. On the other hand, in the HO strain, increases of T cells were dominant in CF, HK and SP at multiple time points. The differences between HO and TL were most distinct at 2, 4, 12 and 48 hours post-exposure (hpe) as well as at 4 days post-exposure (dpe), with the vast majority of innate immune response cells having higher values in the susceptible TL strain. Alteration of the leukocyte populations with augmented local cellular responses and excessive immune reactions likely lead to subsequent host tissue damage and supports parasite invasion and development in TL.

Conclusions

The findings of this study highlight the significance of effective local and systemic immune reaction and indicate proper activation of T lymphocytes critical for host resistance during M. cerebralis infection. The present study provides insights into the cellular basis of protective immune responses against M. cerebralis and can help us to elucidate the mechanisms underlying the variation in resistance to WD.

Major Histocompatibility Complex (MHC) Genes and Disease Resistance in Fish

Fascinating about classical major histocompatibility complex (MHC) molecules is their polymorphism. The present study is a review and discussion of the fish MHC situation. The basic pattern of MHC variation in fish is similar to mammals, with MHC class I versus class II, and polymorphic classical versus nonpolymorphic nonclassical. However, in many or all teleost fishes, important differences with mammalian or human MHC were observed: (1) The allelic/haplotype diversification levels of classical MHC class I tend to be much higher than in mammals and involve structural positions within but also outside the peptide binding groove; (2) Teleost fish classical MHC class I and class II loci are not linked. The present article summarizes previous studies that performed quantitative trait loci (QTL) analysis for mapping differences in teleost fish disease resistance, and discusses them from MHC point of view. Overall, those QTL studies suggest the possible importance of genomic regions including classical MHC class II and nonclassical MHC class I genes, whereas similar observations were not made for the genomic regions with the highly diversified classical MHC class I alleles. It must be concluded that despite decades of knowing MHC polymorphism in jawed vertebrate species including fish, firm conclusions (as opposed to appealing hypotheses) on the reasons for MHC polymorphism cannot be made, and that the types of polymorphism observed in fish may not be explained by disease-resistance models alone.

Genome-wide association analysis of intra-specific QTL associated with the resistance for enteric septicemia of catfish

Disease resistance is one of the most important traits for aquaculture industry. For catfish industry, enteric septicemia of catfish (ESC), caused by the bacterial pathogen Edwardsiella ictaluri, is the most severe disease, causing enormous economic losses every year. In this study, we used three channel catfish families with 900 individuals (300 fish per family) and the 690K catfish SNP array, and conducted a genome-wide association study to detect the quantitative trait loci (QTL) associated with ESC resistance. Three significant QTL, with two of located on LG1 and one on LG26, and three suggestive QTL located on LG1, LG3, and LG21, respectively, were identified to be associated with ESC resistance. With a well-assembled- and -annotated reference genome sequence, genes around the involved QTL regions were identified. Among these genes, 37 genes had known functions in immunity, which may be involved in ESC resistance. Notably, nlrc3 and nlrp12 identified here were also found in QTL regions of ESC resistance in the channel catfish × blue catfish interspecific hybrid system, suggesting this QTL was operating within both intra-specific channel catfish populations and interspecific hybrid backcross populations. Many of the genes of the Class I MHC pathway, for mediated antigen processing and presentation, were found in the QTL regions. The positional correlation found in this study and the expressional correlation found in previous studies indicated that Class I MHC pathway was significantly associated with ESC resistance. This study validated one QTL previously identified using the second and fourth generation of the interspecific hybrid backcross progenies, and identified five additional QTL among channel catfish families. Taken together, it appears that there are only a few major QTL for ESC disease resistance, making marker-assisted selection an effective approach for genetic improvements of ESC resistance.

GWAS analysis using interspecific backcross progenies reveals superior blue catfish alleles responsible for strong resistance against enteric septicemia of catfish

Infectious diseases pose significant threats to the catfish industry. Enteric septicemia of catfish (ESC) caused by Edwardsiella ictaluri is the most devastating disease for catfish aquaculture, causing huge economic losses annually. Channel catfish and blue catfish exhibit great contrast in resistance against ESC, with channel catfish being highly susceptible and blue catfish being highly resistant. As such, the interspecific backcross progenies provide an ideal system for the identification of quantitative trait locus (QTL). We previously reported one significant QTL on linkage group (LG) 1 using the third-generation backcrosses, but the number of founders used to make the second- and third-generation backcross progenies was very small. Although the third-generation backcross progenies provided a greater power for fine mapping than the first-generation backcrosses, some major QTL for disease resistance may have been missing due to the small numbers of founders used to produce the higher generation backcrosses. In this study, we performed a genome-wide association study using first-generation backcrosses with the catfish 690 K SNP arrays to identify additional ESC disease resistance QTL, especially those at the species level. Two genomic regions on LG1 and LG23 were determined to be significantly associated with ESC resistance as revealed by a mixed linear model and family-based association test. Examination of the resistance alleles indicated their origin from blue catfish, indicating that at least two major disease resistance loci exist among blue catfish populations. Upon further validation, markers linked with major ESC disease resistance QTL should be useful for marker-assisted introgression, allowing development of highly ESC resistant breeds of catfish.

Genetically influenced resistance to stress and disease in salmonids in relation to present-day breeding practice - a short review

While intensive fish production has many advantages, it also has a number of drawbacks as regards disease and stress. To date, there has been no conclusive review of disease resistance at Czech fish farms. The aim of the study was to describe briefly the existing salmonid breeding practice in the Czech Republic and to point out the trends and new possibilities gaining ground around Europe. However, the present situation in the Czech stocks is not rare at all and therefore it is used here as a model example representing numerous breeding practices in Europe. Stress and disease resistance in fish is polygenic and quantitative, making selection for such traits difficult. In recent years, however, fish breeding methods have developed rapidly, with the use of genetic analysis tools, for example, now allowing much greater selection accuracy. Gradual progress in understanding the importance of individual genetic markers offers many new options that can be utilised in breeding practice. New selection methods, such as quantitative trait loci (QTLs) and genomic selection, are increasingly employed in European aquaculture. Next generation sequencing techniques now help in the finding of new and promising QTLs that can be used in assisted selection. This review maps the current progress in improving salmonid resistance to stress and disease in aquaculture and at the same time provides the breeders with a short overview of the latest tools of genetically controlled breeding and of the newest products available at the European market.

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.

Genome-Wide Association and Genomic Selection for Resistance to Amoebic Gill Disease in Atlantic Salmon

Amoebic gill disease (AGD) is one of the largest threats to salmon aquaculture, causing serious economic and animal welfare burden. Treatments can be expensive and environmentally damaging, hence the need for alternative strategies. Breeding for disease resistance can contribute to prevention and control of AGD, providing long-term cumulative benefits in selected stocks. The use of genomic selection can expedite selection for disease resistance due to improved accuracy compared to pedigree-based approaches. The aim of this work was to quantify and characterize genetic variation in AGD resistance in salmon, the genetic architecture of the trait, and the potential of genomic selection to contribute to disease control. An AGD challenge was performed in ∼1,500 Atlantic salmon, using gill damage and amoebic load as indicator traits for host resistance. Both traits are heritable (h2 ∼0.25-0.30) and show high positive correlation, indicating they may be good measurements of host resistance to AGD. While the genetic architecture of resistance appeared to be largely polygenic in nature, two regions on chromosome 18 showed suggestive association with both AGD resistance traits. Using a cross-validation approach, genomic prediction accuracy was up to 18% higher than that obtained using pedigree, and a reduction in marker density to ∼2,000 SNPs was sufficient to obtain accuracies similar to those obtained using the whole dataset. This study indicates that resistance to AGD is a suitable trait for genomic selection, and the addition of this trait to Atlantic salmon breeding programs can lead to more resistant stocks.

A novel C-type lectin gene is a strong candidate gene for Benedenia disease resistance in Japanese yellowtail, Seriola quinqueradiata

Little is known about mechanisms of resistance to parasitic diseases in marine finfish. Benedenia disease is caused by infection by the monogenean parasite Benedenia seriolae. Previous quantitative trait locus (QTL) analyses have identified a major QTL associated with resistance to Benedenia disease in linkage group Squ2 of the Japanese yellowtail/amberjack Seriola quinqueradiata. To uncover the bioregulatory mechanism of Benedenia disease resistance, complete Illumina sequencing of BAC clones carrying genomic DNA for the QTL region in linkage group Squ2 was performed to reveal a novel C-type lectin in this region. Expression of the mRNA of this C-type lectin was detected in skin tissue parasitized by B. seriolae. Scanning for single nucleotide polymorphisms (SNPs) uncovered a SNP in the C-type lectin/C-type lectin-like domain that was significantly associated with B. seriolae infection levels. These results strongly suggest that the novel C-type lectin gene controls resistance to Benedenia disease in Japanese yellowtails.

QTL Mapping for Resistance to Iridovirus in Asian Seabass Using Genotyping-by-Sequencing

Identifying quantitative trait loci (QTL) for viral disease resistance is of particular importance in selective breeding programs of fish species. Genetic markers linked to QTL can be useful in marker-assisted selection (MAS) for elites resistant to specific pathogens. Here, we conducted a genome scan for QTL associated with Singapore grouper iridovirus (SGIV) resistance in an Asian seabass (Lates calcarifer) family, using a high-density linkage map generated with genotyping-by-sequencing. One genome-wide significant and three suggestive QTL were detected at LG21, LG6, LG13, and LG15, respectively. The phenotypic variation explained (PVE) by the four QTL ranged from 7.5 to 15.6%. The position of the most significant QTL at LG21 was located between 31.88 and 36.81 cM. The SNP marker (SNP130416) nearest to the peak of this QTL was significantly associated with SGIV resistance in an unrelated multifamily population. One candidate gene, MECOM, close to the peak of this QTL region, was predicted. Evidence of alternative splicing was observed for MECOM and one specific category of splicing variants was differentially expressed at 5 days post-SGIV infection. The QTL detected in this study are valuable resources and can be used in the selective breeding programs of Asian seabass with regard to resistance to SGIV.

Quantitative Genetic Variation in, and Environmental Effects on, Pathogen Resistance and Temperature-Dependent Disease Severity in a Wild Trout

Health after pathogen contact varies among individuals because of differences in pathogen load (which is limited by resistance) and disease severity in response to pathogen load (which is limited by tolerance). To understand pathogen-induced host evolution, it is critical to know not only the relative contributions of nongenetic and genetic variation to resistance and tolerance but also how they change environmentally. We quantified nongenetic and genetic variation in parasite load and the associated temperature-dependent disease among trout siblings from two rivers. We detected a genetic variance for parasite load 6.6 times as large in the colder river. By contrast, genetic variance for disease traits tended to be larger in the warmer river, where the disease was manifested more severely. The relationships between disease severity and pathogen load (tolerance) exhibited plateaus at low pathogen load and stronger steepening slopes at high pathogen load in the warmer river. Our study demonstrates the environmental influence on disease severity, nongenetic and genetic variance for health-damage-limiting host abilities, and the shape of tolerance curves. Environmental variability is predicted to govern the presence and intensity of selection, change the relative contributions of nongenetic and genetic variance, and therefore hamper evolution toward more resistant and tolerant hosts.

A RNAi-based therapeutic proof of concept targets salmonid whirling disease in vivo

Myxobolus cerebralis is a cnidarian-myxozoan parasite that causes salmonid whirling disease. M. cerebralis alternates between two hosts: (1) a vertebrate salmonid and (2) an invertebrate oligochaete, Tubifex tubifex. There is no successful treatment for salmonid whirling disease. MyxSP-1 is a M. cerebralis serine protease implicated in whirling disease pathogenesis. We hypothesized that short-interfering RNA (siRNA)-induced RNA interference (RNAi) can silence MyxSP-1 in the invertebrate host and abrogate the M. cerebralis life cycle. This would preclude whirling disease infection in the salmonid host. To test this hypothesis, we first developed a siRNA delivery protocol in T. tubifex. Second, we determined the effective dose for siRNA treatment of M. cerebralis-infected T. tubifex. M. cerebralis-infected T. tubifex were treated with different concentrations of MyxSP-1 or negative control siRNAs (1μM, 2μM, 5μM or 7μM) at 15°C for 24h, 48h, 72h and 96h, respectively. We monitored MyxSP-1 knockdown using real-time quantitative PCR (qPCR). siRNA treatment with MyxSP-1 siRNA at 2μM concentration for 24h at 15°C showed maximum significant MyxSP-1 knockdown in T. tubifex. Third, we determined the time points in the M. cerebralis life cycle in T. tubifex at which siRNA treatment was most effective. M. cerebralis-infected T. tubifex were treated with MyxSP-1 or negative control siRNAs (2μM concentration for 24h at 15°C) at 24 hours post-infection (24hpi), 48hpi, 72hpi, 96hpi, 1 month post-infection (1mpi), 2mpi and 3mpi, respectively. We observed that siRNA treatment of T. tubifex was most effective at 1mpi, 2mpi and 3mpi. Fourth, we immersed specific-pathogen-free rainbow trout fry in water inhabited by MyxSP-1 siRNA-treated T. tubifex (at 1mpi, 2mpi and 3mpi). The salmonids did not develop whirling disease and showed significant MyxSP-1 knockdown. We also observed long-term RNAi in T. tubifex. Together these results demonstrate a novel RNAi-based therapeutic proof of concept in vivo against salmonid whirling disease.

Studying the Genetics of Resistance to CyHV-3 Disease Using Introgression from Feral to Cultured Common Carp Strains

Sustainability and further development of aquaculture production are constantly challenged by outbreaks of fish diseases, which are difficult to prevent or control. Developing fish strains that are genetically resistant to a disease is a cost-effective and a sustainable solution to address this challenge. To do so, heritable genetic variation in disease resistance should be identified and combined together with other desirable production traits. Aquaculture of common carp has suffered substantial losses from the infectious disease caused by the cyprinid herpes virus type 3 (CyHV-3) virus and the global spread of outbreaks indicates that many cultured strains are susceptible. In this research, CyHV-3 resistance from the feral strain “Amur Sassan” was successfully introgressed into two susceptible cultured strains up to the first backcross (BC₁) generation. Variation in resistance of families from F₁ and BC₁ generations was significantly greater compared to that among families of any of the susceptible parental lines, a good starting point for a family selection program. Considerable additive genetic variation was found for CyHV-3 resistance. This phenotype was transferable between generations with contributions to resistance from both the resistant feral and the susceptible cultured strains. Reduced scale coverage (mirror phenotype) is desirable and common in cultured strains, but so far, cultured mirror carp strains were found to be susceptible. Here, using BC₁ families ranging from susceptible to resistant, no differences in resistance levels between fully scaled and mirror full-sib groups were found, indicating that CyHV-3 resistance was successfully combined with the desirable mirror phenotype. In addition, the CyHV-3 viral load in tissues throughout the infection of susceptible and resistant fish was followed. Although resistant fish get infected, viral loads in tissues of these fish are significantly lesser than in those of susceptible fish, allowing them to survive the disease. Taken together, in this study we have laid the foundation for breeding CyHV-3-resistant strains and started to address the mechanisms underlying the phenotypic differences in resistance to this disease.

Aquaculture genomics, genetics and breeding in the United States: current status, challenges, and priorities for future research

Advancing the production efficiency and profitability of aquaculture is dependent upon the ability to utilize a diverse array of genetic resources. The ultimate goals of aquaculture genomics, genetics and breeding research are to enhance aquaculture production efficiency, sustainability, product quality, and profitability in support of the commercial sector and for the benefit of consumers. In order to achieve these goals, it is important to understand the genomic structure and organization of aquaculture species, and their genomic and phenomic variations, as well as the genetic basis of traits and their interrelationships. In addition, it is also important to understand the mechanisms of regulation and evolutionary conservation at the levels of genome, transcriptome, proteome, epigenome, and systems biology. With genomic information and information between the genomes and phenomes, technologies for marker/causal mutation-assisted selection, genome selection, and genome editing can be developed for applications in aquaculture. A set of genomic tools and resources must be made available including reference genome sequences and their annotations (including coding and non-coding regulatory elements), genome-wide polymorphic markers, efficient genotyping platforms, high-density and high-resolution linkage maps, and transcriptome resources including non-coding transcripts. Genomic and genetic control of important performance and production traits, such as disease resistance, feed conversion efficiency, growth rate, processing yield, behaviour, reproductive characteristics, and tolerance to environmental stressors like low dissolved oxygen, high or low water temperature and salinity, must be understood. QTL need to be identified, validated across strains, lines and populations, and their mechanisms of control understood. Causal gene(s) need to be identified. Genetic and epigenetic regulation of important aquaculture traits need to be determined, and technologies for marker-assisted selection, causal gene/mutation-assisted selection, genome selection, and genome editing using CRISPR and other technologies must be developed, demonstrated with applicability, and application to aquaculture industries.Major progress has been made in aquaculture genomics for dozens of fish and shellfish species including the development of genetic linkage maps, physical maps, microarrays, single nucleotide polymorphism (SNP) arrays, transcriptome databases and various stages of genome reference sequences. This paper provides a general review of the current status, challenges and future research needs of aquaculture genomics, genetics, and breeding, with a focus on major aquaculture species in the United States: catfish, rainbow trout, Atlantic salmon, tilapia, striped bass, oysters, and shrimp. While the overall research priorities and the practical goals are similar across various aquaculture species, the current status in each species should dictate the next priority areas within the species. This paper is an output of the USDA Workshop for Aquaculture Genomics, Genetics, and Breeding held in late March 2016 in Auburn, Alabama, with participants from all parts of the United States.

Mapping QTL for Resistance Against Viral Nervous Necrosis Disease in Asian Seabass

Viral nervous necrosis disease (VNN), caused by nervous necrosis virus (NNV), leads to mass mortality in mariculture. However, phenotypic selection for resistance against VNN is very difficult. To facilitate marker-assisted selection (MAS) for resistance against VNN and understanding of the genetic architecture underlying the resistance against this disease, we mapped quantitative trait loci (QTL) for resistance against VNN in Asian seabass. We challenged fingerlings at 37 days post-hatching (dph), from a single back-cross family, with NNV at a concentration of 9 × 10(6) TCID50/ml for 2 h. Daily mortalities were recorded and collected. A panel of 330 mortalities and 190 surviving fingerlings was genotyped using 149 microsatellites with 145 successfully mapped markers covering 24 linkage groups (LGs). Analysis of QTL for both resistance against VNN and survival time was conducted using interval mapping. Five significant QTL located in four LGs and eight suggestive QTL in seven LGs were identified for resistance. Another five significant QTL in three LGs and five suggestive QTL in three LGs were detected for survival time. One significant QTL, spanning 3 cM in LG20, was identified for both resistance and survival time. These QTL explained 2.2-4.1% of the phenotypic variance for resistance and 2.2-3.3% of the phenotypic variance for survival time, respectively. Our results suggest that VNN resistance in Asian seabass is controlled by many loci with small effects. Our data provide information for fine mapping of QTL and identification of candidate genes for a better understanding of the mechanism of disease resistance.

Whirling disease revisited: pathogenesis, parasite biology and disease intervention

Whirling disease (WD) is an ecologically and economically debilitating disease of rainbow trout Oncorhynchus mykiss caused by the actinosporean spores of the parasite Myxobolus cerebralis. M. cerebralis has a complex, 2-host life cycle alternating between salmonid fish and the oligochaete host Tubifex tubifex. The parasite alternates between 2 spore forms as transmission stages: an actinosporean triactinomyxon spore that is produced in the oligochaete host and a myxosporean spore that develops in the salmonid host. Waterborne triactinomyxon spores released from infected T. tubifex oligochaetes attach to the salmonid host by polar filament extrusion elicited by chemical (nucleoside) and mechanical (thigmotropy) stimuli-a process which is rapidly followed by active penetration of the sporoplasms into the fish epidermis. Upon penetration, sporoplasms multiply and migrate via peripheral nerves and the central nervous system to reach the cartilage where they form trophozoites which undergo further multiplication and subsequent sporogenesis. M. cerebralis myxospores are released into the aquatic environment when infected fish die and autolyse, or when they are consumed and excreted by predators. Myxospores released into the water are ingested by susceptible T. tubifex where they develop intercellularly in the intestine over a period of 3 mo through 4 developmental stages to give rise to mature actinospores. In this article, we review our current understanding of WD-the parasite and its alternate hosts, life cycle and development of the parasite in either host, disease distribution, susceptibility and resistance mechanisms in salmonid host and strategies involved in diagnosis, prevention and control of WD.

Parental genetic diversity of brown trout (Salmo trutta m. fario) brood stock affects offspring susceptibility to whirling disease

BACKGROUND

Whirling disease, caused by the myxozoan parasite Myxobolus cerebralis, has high economical and ecological importance worldwide. Susceptibility to the disease varies considerably among salmonid species. In brown trout (Salmo trutta) the infection is usually subclinical with low mortality, which increases the risk of parasite dissemination, especially when farm fish are used for stocking natural habitats. The influence of intraspecific genetic differences (especially the level of homozygosity) on susceptibility is unknown. Therefore, we examined the possible correlations between parental genetic diversity and offspring susceptibility of brown trout stocks to whirling disease.

METHODS

Two brown trout brood stocks from a German and a Hungarian fish farm were genetically characterized using microsatellite and lineage-specific genetic markers. The individual inbreeding coefficient f and pairwise relatedness factor r were estimated based on eight microsatellite markers. Brood stock populations were divided into groups according to low and high f and r value estimates and subjected to selective fertilization. The offspring from these separate groups were exposed to M. cerebralis actinospores, and the infection prevalence and intensity was measured and statistically analysed.

RESULTS

The analysis of phylogeographic lineage heritage revealed high heterogeneity in the Hungarian brood stock since > 50% of individuals were Atlantic-Danubian hybrids, while only pure Atlantic-descending specimens were detected in the German population. Based on f msat and r msat estimations, classified non-inbred (NIB), inbred (IB) and a group of closely related fish (REL) were created. The susceptibility of their offspring varied considerably. Although there was no significant difference in the prevalence of M. cerebralis infection, the mean intensity of infection differed significantly between NIB and IB groups. In REL and IB groups, a high variability was observed in infection intensity. No external clinical signs were observed in the exposed brown trout groups.

CONCLUSIONS

Our findings indicate that the allelic diversity of brown trout brood stock may constitute a significant factor in disease susceptibility, i.e. the intensity of parasite infection in the subsequent generation.

Genetic and genomic analyses for economically important traits and their applications in molecular breeding of cultured fish

The traits of cultured fish must continually be genetically improved to supply high-quality animal protein for human consumption. Economically important fish traits are controlled by multiple gene quantitative trait loci (QTL), most of which have minor effects, but a few genes may have major effects useful for molecular breeding. In this review, we chose relevant studies on some of the most intensively cultured fish and concisely summarize progress on identifying and verifying QTLs for such traits as growth, disease and stress resistance and sex in recent decades. The potential applications of these major-effect genes and their associated markers in marker-assisted selection and molecular breeding, as well as future research directions are also discussed. These genetic and genomic analyses will be valuable for elucidating the mechanisms modulating economically important traits and to establish more effective molecular breeding techniques in fish.

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.

The development and characterization of a 57K single nucleotide polymorphism array for rainbow trout

In this study, we describe the development and characterization of the first high-density single nucleotide polymorphism (SNP) genotyping array for rainbow trout. The SNP array is publically available from a commercial vendor (Affymetrix). The SNP genotyping quality was high, and validation rate was close to 90%. This is comparable to other farm animals and is much higher than previous smaller scale SNP validation studies in rainbow trout. High quality and integrity of the genotypes are evident from sample reproducibility and from nearly 100% agreement in genotyping results from other methods. The array is very useful for rainbow trout aquaculture populations with more than 40 900 polymorphic markers per population. For wild populations that were confounded by a smaller sample size, the number of polymorphic markers was between 10 577 and 24 330. Comparison between genotypes from individual populations suggests good potential for identifying candidate markers for populations' traceability. Linkage analysis and mapping of the SNPs to the reference genome assembly provide strong evidence for a wide distribution throughout the genome with good representation in all 29 chromosomes. A total of 68% of the genome scaffolds and contigs were anchored through linkage analysis using the SNP array genotypes, including ~20% of the genome assembly that has not been previously anchored to chromosomes.

Biochemical reference intervals and pathophysiological changes in Flavobacterium psychrophilum-resistant and -susceptible rainbow trout lines

Host genetic resistance against disease-causing pathogens can be enhanced through family-based selective breeding. At present, there is an incomplete understanding of how artificial selection of fish alters host physiology and response following pathogen exposure. We previously reported the generation of selectively-bred rainbow trout Oncorhynchus mykiss lines with either increased resistance (ARS-Fp-R) or susceptibility (ARS-Fp-S) to bacterial cold water disease (BCWD). This study (1) determined baseline reference-range intervals for packed cell volume (PCV) and 18 plasma biochemistry analytes, and (2) examined pathophysiological changes following infection between the genetic lines. PCV and biochemistry reference-range intervals did not significantly differ between genetic lines; thus data were pooled into a single reference-range population (n = 85). ARS-Fp-R and ARS-Fp-S line fish were intraperitoneally challenged with Flavobacterium psychrophilum, and plasma was collected on Days 1, 3, 6, and 9 post-challenge. Splenic bacterial load was measured using an F. psychrophilum-specific qPCR assay. In both genetic lines, changes were observed in mean PCV, total protein, albumin, glucose, cholesterol, chloride, and calcium, falling outside the established reference intervals and significantly differing from phosphate-buffered saline challenged fish, on at least 1d post-challenge. Mean PCV, total protein, and calcium significantly differed between ARS-Fp-R and ARS-Fp-S line fish on Day 9 post-infection, with values in the ARS-Fp-S line deviating most from the reference interval. PCV, total protein, cholesterol, and calcium negatively correlated with bacterial load. These findings identify divergent pathophysiological responses between ARS-Fp-R and ARS-Fp-S line fish following laboratory challenge that are likely associated with differential survival.

Histopathologic Changes in Disease-Resistant-Line and Disease-Susceptible-Line Juvenile Rainbow Trout Experimentally Infected with Flavobacterium psychrophilum

Abstract A challenge to selectively breeding fish populations for improved disease resistance in aquaculture is an incomplete understanding of how artificial selection affects innate immunity at the host-pathogen level. The objective of this study was to determine whether Rainbow Trout Oncorhynchus mykiss bred for differential susceptibility to bacterial cold-water disease exhibited altered tissue damage and cellular inflammatory response following experimental challenge with Flavobacterium psychrophilum. Fish from disease-resistant (ARS-Fp-R) and disease-susceptible (ARS-Fp-S) lines were experimentally challenged as juveniles, and mortalities, as well as survivors, were sampled for histopathology during the acute phase of the disease. Microscopic lesions were quantified or semiquantified and statistically compared for changes over time and between genetic lines. Significant progression in the degree of perisplenitis, splenic necrosis, splenic inflammatory infiltrates, average splenic ellipsoid area, total splenic ellipsoid area, and peritonitis was present over time in both genetic lines on at least one postinfection time point. No differences were found between renal inflammatory infiltrates and renal hematopoietic cell depletion over time. Perisplenitis was significantly lower in fish from the ARS-Fp-R line on day 9 postinfection than in fish from the ARS-Fp-S line. The ARS-Fp-R line demonstrated a trend towards reduced splenic necrosis compared with the ARS-Fp-S line that approached significance, and fish from the ARS-Fp-S line were 3.6 times more likely than fish from the ARS-Fp-R line to have a higher splenic necrosis lesion score after day 3 postinfection. These findings support the hypothesis that differential survival is a result of divergence in disease magnitude and not altered disease course between genetic lines. Characterization of histopathologic changes between genetic lines and over time helps elucidate mechanisms of disease resistance and contributes to our understanding of disease pathogenesis in fish infected with F. psychrophilum. Received January 7, 2014; accepted March 10, 2014.

Photobacteriosis: prevention and diagnosis

Photobacteriosis or fish pasteurellosis is a bacterial disease affecting wild and farm fish. Its etiological agent, the gram negative bacterium Photobacterium damselae subsp. piscicida, is responsible for important economic losses in cultured fish worldwide, in particular in Mediterranean countries and Japan. Efforts have been focused on gaining a better understanding of the biology of the pathogenic microorganism and its natural hosts with the aim of developing effective vaccination strategies and diagnostic tools to control the disease. Conventional vaccinology has thus far yielded unsatisfactory results, and recombinant technology has been applied to identify new antigen candidates for the development of subunit vaccines. Furthermore, molecular methods represent an improvement over classical microbiological techniques for the identification of P. damselae subsp. piscicida and the diagnosis of the disease. The complete sequencing, annotation, and analysis of the pathogen genome will provide insights into the pathogen laying the groundwork for the development of vaccines and diagnostic methods.

Survival and reproduction of Myxobolus cerebralis-resistant rainbow trout introduced to the Colorado river and increased resistance of age-0 progeny

Myxobolus cerebralis caused severe declines in rainbow trout populations across Colorado following its introduction in the 1980s. One promising approach for the recovery of Colorado's rainbow trout populations has been the production of rainbow trout that are genetically resistant to the parasite. We introduced one of these resistant crosses, known as the GR×CRR (cross between the German Rainbow [GR] and Colorado River Rainbow [CRR] trout strains), to the upper Colorado River. The abundance, survival, and growth of the stocked GR×CRR population was examined to determine if GR×CRRs had contributed offspring to the age-0 population, and determine whether these offspring displayed increased resistance and survival characteristics compared to their wild CRR counterparts. Apparent survival of the introduced GR×CRR over the entire study period was estimated to be 0.007 (±0.001). Despite low survival of the GR×CRRs, age-0 progeny of the GR×CRR were encountered in years 2008 through 2011. Genetic assignments revealed a shift in the genetic composition of the rainbow trout fry population over time, with CRR fish comprising the entirety of the fry population in 2007, and GR-cross fish comprising nearly 80% of the fry population in 2011. A decrease in average infection severity (myxospores fish-1) was observed concurrent with the shift in the genetic composition of the rainbow trout fry population, decreasing from an average of 47,708 (±8,950) myxospores fish-1 in 2009 to 2,672 (±4,379) myxospores fish-1 in 2011. Results from this experiment suggest that the GR×CRR can survive and reproduce in rivers with a high prevalence of M. cerebralis. In addition, reduced myxospore burdens in age-0 fish indicated that stocking this cross may ultimately lead to an overall reduction in infection prevalence and severity in the salmonid populations of the upper Colorado River.

Biology and mucosal immunity to myxozoans

Myxozoans are among the most abundant parasites in nature. Their life cycles involve two hosts: an invertebrate, usually an annelid, and a vertebrate, usually a fish. They affect fish species in their natural habitats but also constitute a menace for fish aquaculture. Using different strategies they are able to parasitize and cause damage in multiple organs, including mucosal tissues, which they use also as portals of entry. In fish, the main mucosal sites include the intestine, skin and gills. Recently the finding of a specific mucosal immunoglobulin in teleost (IgT), analogous to mammalian IgA, and the capacity of fish to develop a specific mucosal immune response against different pathogens, has highlighted the importance of studying immune responses at mucosal sites. In this review, we describe the major biological characteristics of myxozoan parasites and present the data available regarding immune responses for species that infect mucosal sites. As models for mucosal immunity we review the responses to Enteromyxum spp. and Ceratomyxa shasta, both of which parasitize the intestine. The immune response at the skin and gills is also described, as these mucosal tissues are used by myxozoans as attaching surfaces and portal of entry, and some species also parasitize these sites. Finally, the development of immunoprophylactic strategies is discussed.

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.

Temporal expression patterns of rainbow trout immune-related genes in response to Myxobolus cerebralis exposure

Infection of salmonids by the myxozoan parasite Myxobolus cerebralis can cause whirling disease, which is responsible for high mortalities in rainbow trout hatcheries and natural populations in the United States. Although considerable research has provided insight into disease pathology, host invasion, and inheritance patterns of resistance, the causal genetic variants and molecular mechanisms underlying host resistance or susceptibility remain elusive. A previous study found that expression changes of specific metallothionein genes following M. cerebralis infection are implicated in whirling disease resistance. The present study examines the dynamic transcriptional response to infection of several upstream regulators of the metallothionein gene family (IL-1β, KLF2, STAT3, STAT5), along with innate immune response genes (IFN-γ, IRF1 and iNOS). Pathogen loads and gene expression were compared across multiple time points after M. cerebralis exposure to elucidate how resistant and susceptible rainbow trout strains transcriptionally respond to early invasion. IL-1β, IFN-γ, IRF1, and iNOS all showed increased expression following M. cerebralis exposure for one or both strains across multiple time points. The interferon-related genes IFN-γ and IRF1 had consistently increased expression in the susceptible strain in comparison to the resistant strain, likely due to a less effective initial immune response. STAT3 was the only gene with consistently increased expression in the resistant strain following infection while remaining unchanged in the susceptible strain. Given its pleiotropic effects on immune response, STAT3 is an excellent candidate for future research of whirling disease resistance mechanisms.

Quantitative trait loci (QTL) associated with resistance to a monogenean parasite (Benedenia seriolae) in yellowtail (Seriola quinqueradiata) through genome wide analysis

Benedenia infections caused by the monogenean fluke ectoparasite Benedenia seriolae seriously impact marine finfish aquaculture. Genetic variation has been inferred to play a significant role in determining the susceptibility to this parasitic disease. To evaluate the genetic basis of Benedenia disease resistance in yellowtail (Seriola quinqueradiata), a genome-wide and chromosome-wide linkage analyses were initiated using F₁ yellowtail families (n = 90 per family) based on a high-density linkage map with 860 microsatellite and 142 single nucleotide polymorphism (SNP) markers. Two major quantitative trait loci (QTL) regions on linkage groups Squ2 (BDR-1) and Squ20 (BDR-2) were identified. These QTL regions explained 32.9–35.5% of the phenotypic variance. On the other hand, we investigated the relationship between QTL for susceptibility to B. seriolae and QTL for fish body size. The QTL related to growth was found on another linkage group (Squ7). As a result, this is the first genetic evidence that contributes to detailing phenotypic resistance to Benedenia disease, and the results will help resolve the mechanism of resistance to this important parasitic infection of yellowtail.

Quantitative trait loci for morphometric body measurements of the hybrids of silver carp (Hypophthalmichthys molitrix) and bighead carp (H. nobilis)

Quantitative trait loci (QTL) for 11 morphometric body measurements of the hybrids of silver (Hypophthalmichthys molitrix) and bighead carp (H. nobilis) including body weight (BW), standard length (SL), body depth (BD), body thickness (BT), head length (HL), head depth (HD), length of ventral keel (LVK), length of pectoral fin (Lpec), length of pelvic fin (Lpel), length of caudal fin (Lcau) and space between pectoral and pelvic fins (SPP) were located on the sex average microsatellite linkage map constructed using the hybrids of a female bighead and a male silver carp, on which 15 microsatellites were newly mapped. One locus was found to be responsible for BW, LV K and SPP, respectively. As many as 6 loci were found to be responsible for HD. The variances of remaining traits were partitioned by different numbers of loci varying between 2 and 5. The variance explained each locus ranged from 9.1% to 23.8% of the total. The variance explained by all loci responsible for each measurement ranged from 17.7% to 75.1%. It was noted that multiple measurements were mapped on the same locus. For example, a region bounded by Hym435 and Hym145 was found to be responsible for all the measurements analyzed.

Resistance to a rhabdovirus (VHSV) in rainbow trout: identification of a major QTL related to innate mechanisms

Health control is a major issue in animal breeding and a better knowledge of the genetic bases of resistance to diseases is needed in farm animals including fish. The detection of quantitative trait loci (QTL) will help uncovering the genetic architecture of important traits and understanding the mechanisms involved in resistance to pathogens. We report here the detection of QTL for resistance to Viral Haemorrhagic Septicaemia Virus (VHSV), a major threat for European aquaculture industry. Two induced mitogynogenetic doubled haploid F2 rainbow trout (Oncorhynchus mykiss) families were used. These families combined the genome of susceptible and resistant F0 breeders and contained only fully homozygous individuals. For phenotyping, fish survival after an immersion challenge with the virus was recorded, as well as in vitro virus replication on fin explants. A bidirectional selective genotyping strategy identified seven QTL associated to survival. One of those QTL was significant at the genome-wide level and largely explained both survival and viral replication in fin explants in the different families of the design (up to 65% and 49% of phenotypic variance explained respectively). These results evidence the key role of innate defence in resistance to the virus and pave the way for the identification of the gene(s) responsible for resistance. The identification of a major QTL also opens appealing perspectives for selective breeding of fish with improved resistance.

Genetic basis of differences in myxospore count between whirling disease-resistant and susceptible strains of rainbow trout

We used a quantitative genetics approach and estimated broad sense heritability (h2b) of myxospore count and the number of genes involved in myxospore formation to gain a better understanding of how resistance to Myxobolus cerebralis, the parasite responsible for whirling disease, is inherited in rainbow trout Oncorhynchus mykiss. An M. cerebralis-resistant strain of rainbow trout, the German Rainbow (GR), and a wild, susceptible strain of rainbow trout, the Colorado River Rainbow (CRR), were spawned to create 3 intermediate crossed populations (an F1 cross, F2 intercross, and a B2 backcross between the F1 and the CRR). Within each strain or cross, h2b was estimated from the between-family variance of myxospore counts using full-sibling families. Estimates of h2b and average myxospore counts were lowest in the GR strain, F1 cross, and F2 intercross (h2b = 0.34, 0.42, and 0.34; myxospores fish-1 = 275, 9566, and 45780, respectively), and highest in the B2 backcross and CRR strain (h2b = 0.93 and 0.89; myxospores fish-1 = 97865 and 187595, respectively). Comparison of means and a joint-scaling test suggest that resistance alleles arising from the GR strain are dominant to susceptible alleles from the CRR strain. Resistance was retained in the intermediate crosses but decreased as filial generation number increased (F2) or backcrossing occurred (B2). The estimated number of segregating loci responsible for differences in myxospore count in the parental strains was 9 ± 5. Our results indicate that resistance to M. cerebralis is a heritable trait within these populations and would respond to either artificial selection in hatcheries or natural selection in the wild.

Genomics in eels--towards aquaculture and biology

Freshwater eels (genus Anguilla), especially the species inhabiting the temperate areas such as the European, American and Japanese eels, are important aquaculture species. Although artificial reproduction has been attempted since the 1930s and large numbers of studies have been conducted, it has not yet fully succeeded. Problems in eel artificial breeding are highly diverse, for instance, lack of basic information about reproduction in nature, no appropriate food for larvae, high mortality, and high individual variation in adults in response to maturation induction. Over the last decade, genomic data have been obtained for a variety of aquatic organisms. Recent technological advances in sequencing and computation now enable the accumulation of genomic information even for non-model species. The draft genome of the European eel Anguilla anguilla has been recently determined using Illumina technology and transcriptomic data based on next generation sequencing have been emerging. Extensive genomic information will facilitate many aspects of the artificial reproduction of eels. Here, we review the progress in genome-wide studies of eels, including additional analysis of the European eel genome data, and discuss future directions and implications of genomic data for aquaculture.

Biostatic activity of piscine serum and mucus on myxozoan fish infective stages

Since the basis of host specificity in Myxozoa, i.e. the differential disposition and extinction of erroneously penetrated myxozoan infective stages in non-susceptible fish hosts, remains puzzling, we aimed to explore the role of the innate immune system in this issue. In a comparative incubation challenge of actinospore sporoplasms of the freshwater parasite species Myxobolus cerebralis, Henneguya nuesslini and Myxobolus pseudodispar to isolates of host and non-host muci and blood sera, we measured cellular disintegration proportions and times by means of a double staining viability assay utilizing fluorescent dyes. After their activation, emerging primary and secondary sporoplasm cells were evaluated microscopically for physical integrity and onset of cell death due to exposure. Impairment by any mucus used was not detected up to 100 min of exposure. All parasites showed significantly increased cellular breakdown in non-susceptible host serum compared to the respective substrates from susceptible host fish. Except for M. cerebralis, the serum of the susceptible host was considerably less effective over time. In this species, both the primary and the secondary cells were affected in much shorter times than in the other two representatives. Inhibition of protease activity did not affect carp serum effect on M. cerebralis stages. We suggest the active components to be complement or complement induced factors since heat inactivation and withdrawal of bivalent metal ions lowered serum activity significantly. The study marks the first in vitro viability challenge of activated myxozoan transmission stages with teleost derived immune factors.

Gene expression profiles of spleen, liver, and head kidney in turbot (Scophthalmus maximus) along the infection process with Philasterides dicentrarchi using an immune-enriched oligo-microarray

We evaluated the expression profiles of turbot in spleen, liver, and head kidney across five temporal points of the Philasterides dicentrarchi infection process using an 8x15K Agilent oligo-microarray. The microarray included 2,176 different fivefold replicated gene probes designed from a turbot 3' sequenced EST database. We were able to identify 221 differentially expressed (DE) genes (8.1% of the whole microarray), 113 in spleen, 83 in liver, and 90 in head kidney, in at least 1 of the 5 temporal points sampled for each organ. Most of these genes could be annotated (83.0%) and functionally categorized using GO terms (69.1%) after the additional sequencing of DE genes from the 5' end. Many DE genes were related to innate and acquired immune functions. A high proportion of DE genes were organ-specific (70.6%), although their associated GO functions showed notable similarities in the three organs. The most striking difference in functional distribution was observed between the up- and downregulated gene groups. Upregulated genes were mostly associated to immune functions, while downregulated ones mainly involved metabolism-related genes. Genetic response appeared clustered in a few groups of genes with similar expression profiles along the temporal series. The information obtained will aid to understand the turbot immune response and will specifically be valuable to develop strategies of defense to P. dicentrarchi to achieve more resistant broodstocks for turbot industry.

A second generation integrated map of the rainbow trout (Oncorhynchus mykiss) genome: analysis of conserved synteny with model fish genomes

DNA fingerprints and end sequences from bacterial artificial chromosomes (BACs) from two new libraries were generated to improve the first generation integrated physical and genetic map of the rainbow trout (Oncorhynchus mykiss) genome. The current version of the physical map is composed of 167,989 clones of which 158,670 are assembled into contigs and 9,319 are singletons. The number of contigs was reduced from 4,173 to 3,220. End sequencing of clones from the new libraries generated a total of 11,958 high quality sequence reads. The end sequences were used to develop 238 new microsatellites of which 42 were added to the genetic map. Conserved synteny between the rainbow trout genome and model fish genomes was analyzed using 188,443 BAC end sequence (BES) reads. The fractions of BES reads with significant BLASTN hits against the zebrafish, medaka, and stickleback genomes were 8.8%, 9.7%, and 10.5%, respectively, while the fractions of significant BLASTX hits against the zebrafish, medaka, and stickleback protein databases were 6.2%, 5.8%, and 5.5%, respectively. The overall number of unique regions of conserved synteny identified through grouping of the rainbow trout BES into fingerprinting contigs was 2,259, 2,229, and 2,203 for stickleback, medaka, and zebrafish, respectively. These numbers are approximately three to five times greater than those we have previously identified using BAC paired ends. Clustering of the conserved synteny analysis results by linkage groups as derived from the integrated physical and genetic map revealed that despite the low sequence homology, large blocks of macrosynteny are conserved between chromosome arms of rainbow trout and the model fish species.

QTL detection for Aeromonas salmonicida resistance related traits in turbot (Scophthalmus maximus)

Background

Interactions between fish and pathogens, that may be harmless under natural conditions, often result in serious diseases in aquaculture systems. This is especially important due to the fact that the strains used in aquaculture are derived from wild strains that may not have had enough time to adapt to new disease pressures. The turbot is one of the most promising European aquaculture species. Furunculosis, caused by the bacterium Aeromonas salmonicida, produces important losses to turbot industry. An appealing solution is to achieve more robust broodstock, which can prevent or diminish the devastating effects of epizooties. Genomics strategies have been developed in turbot to look for candidate genes for resistance to furunculosis and a genetic map with appropriate density to screen for genomic associations has been also constructed. In the present study, a genome scan for QTL affecting resistance and survival to A. salmonicida in four turbot families was carried out. The objectives were to identify consistent QTL using different statistical approaches (linear regression and maximum likelihood) and to locate the tightest associated markers for their application in genetic breeding strategies.

Results

Significant QTL for resistance were identified by the linear regression method in three linkage groups (LGs 4, 6 and 9) and for survival in two LGs (6 and 9). The maximum likelihood methodology identified QTL in three LGs (5, 6 and 9) for both traits. Significant association between disease traits and genotypes was detected for several markers, some of them explaining up to 17% of the phenotypic variance. We also identified candidate genes located in the detected QTL using data from previously mapped markers.

Conclusions

Several regions controlling resistance to A. salmonicida in turbot have been detected. The observed concordance between different statistical methods at particular linkage groups gives consistency to our results. The detected associated markers could be useful for genetic breeding strategies. A finer mapping will be necessary at the detected QTL intervals to narrow associations and around the closely associated markers to look for candidate genes through comparative genomics or positional cloning strategies. The identification of associated variants at specific genes will be essential, together with the QTL associations detected in this study, for future marker assisted selection programs.