Suggestive association of major histocompatibility IB genetic markers with resistance to bacterial cold water disease in rainbow trout (Oncorhynchus mykiss)

Heat stress induces oxidative stress and weakens the immune system in catfish Clarias magur: Evidence from physiological, histological, and transcriptomic analyses

Climate change is unequivocal, causing a rise in the Earth's temperature, which ultimately impacts all ecosystems. However, aquatic ecosystems are most severely affected by rising temperatures resulting in huge losses to aquaculture industry. The present study investigated the oxidative stress, histopathological changes, and transcriptomic responses in a freshwater catfish Clarias magur subjected to acute heat stress. Fish were exposed to four different temperatures, i.e., 28, 30, 32, and 34 °C, for 96 h to assess their heat tolerance and adaptation behavior. Fish kept at 26 °C were considered the control group. Elevated levels of key antioxidative enzymes such as catalase, glutathione reductase, and superoxide dismutase, were recorded in vital organs (gills, kidney, liver, and rosette). High rates of lipid peroxidation were also observed in the gills, kidney, liver, and rosette. An analysis of the top 25 differentially expressed genes of the gill transcriptome revealed that 72 percent of the transcripts were represented by innate and adaptive immune response genes. Downregulation of BOLA class I and MHC class I molecules indicated impaired immunity whereas, upregulation of MHC class II beta chain and GTPase IMAP8 suggested a compensatory immune response. These findings were also supported by the observed histoarchitectural alterations, such as disintegration of the skin barrier, hepatic and nephrotic apoptosis, tissue hyperplasia, macrophage infiltration, and development of splenic granulomas. This study provides important insights into physiological and molecular mechanisms underlying acute heat stress responses. Understanding these mechanisms is important for developing mitigation strategies to improve the sustainability and resilience of commercially important catfish under continuously changing climatic conditions.

Survival, serum biochemical parameters, hepatic antioxidant status, and gene expression of three Nile tilapia strains under pathogenic Streptococcus agalactiae challenge

Streptococcosis is the leading bacterial disease impacting Nile tilapia (Oreochromis niloticus) and causes substantial economic losses in China. This study assessed the resistance and tolerance of three Nile tilapia (88, 99, and NG) to Streptococcus agalactiae infection. Survival rates were monitored, and samples were collected from blood, liver, and spleen at 0, 1, 2, 5, 7, and 14 days post-infection. Serum biochemical parameters, hepatic antioxidant enzymes, and the expression of related antioxidant and immune genes were measured. Results showed that strain 88 had superior resistance, with the highest survival rate, reduced liver damage, and lower alanine aminotransferase and alkaline phosphatase levels at 1 and 5 days post-infection. This strain also had higher serum cholesterol and triglyceride levels between days 5 and 14. Antioxidant activities increased in all strains post-infection, with strain 88 showing significantly higher superoxide dismutase (SOD) activity and elevated catalase (CAT) and SOD gene expressions; this indicated enhanced control of reactive oxygen species and reduced tissue damage. Additionally, strain 88 exhibited lower levels of inflammatory markers (TNF-α, IL-1β, IL-10) with higher IgM levels and superior MHC-II expression during early infection, which demonstrated a strong immune response. These results indicate that strain 88 is a promising candidate for selective breeding to improve streptococcosis resistance in tilapia.

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.

What Can Genetics Do for the Control of Infectious Diseases in Aquaculture?

Simple Summary

Infectious diseases place an economic burden on aquaculture and a limitation to its growth. This state-of-the-art review describes the application of genetics and genomics as novel tools to control infectious disease in aquaculture.

Abstract

Infectious diseases place an economic burden on aquaculture and a limitation to its growth. An innovative approach to mitigate their impact on production is breeding for disease resistance: selection for domestication, family-based selection, marker-assisted selection, and more recently, genomic selection. Advances in genetics and genomics approaches to the control of infectious diseases are key to increasing aquaculture efficiency, profitability, and sustainability and to reducing its environmental footprint. Interaction and co-evolution between a host and pathogen can, however, turn breeding to boost infectious disease resistance into a potential driver of pathogenic change. Parallel molecular characterization of the pathogen and its virulence and antimicrobial resistance genes is therefore essential to understand pathogen evolution over time in response to host immunity, and to apply appropriate mitigation strategies.

Toward Genome-Based Selection in Asian Seabass: What Can We Learn From Other Food Fishes and Farm Animals?

Due to the steadily increasing need for seafood and the plateauing output of fisheries, more fish need to be produced by aquaculture production. In parallel with the improvement of farming methods, elite food fish lines with superior traits for production must be generated by selection programs that utilize cutting-edge tools of genomics. The purpose of this review is to provide a historical overview and status report of a selection program performed on a catadromous predator, the Asian seabass (Lates calcarifer, Bloch 1790) that can change its sex during its lifetime. We describe the practices of wet lab, farm and lab in detail by focusing onto the foundations and achievements of the program. In addition to the approaches used for selection, our review also provides an inventory of genetic/genomic platforms and technologies developed to (i) provide current and future support for the selection process; and (ii) improve our understanding of the biology of the species. Approaches used for the improvement of terrestrial farm animals are used as examples and references, as those processes are far ahead of the ones used in aquaculture and thus they might help those working on fish to select the best possible options and avoid potential pitfalls.

Salmonid Antibacterial Immunity: An Aquaculture Perspective

The aquaculture industry is continuously threatened by infectious diseases, including those of bacterial origin. Regardless of the disease burden, aquaculture is already the main method for producing fish protein, having displaced capture fisheries. One attractive sector within this industry is the culture of salmonids, which are (a) uniquely under pressure due to overfishing and (b) the most valuable finfish per unit of weight. There are still knowledge gaps in the understanding of fish immunity, leading to vaccines that are not as effective as in terrestrial species, thus a common method to combat bacterial disease outbreaks is the use of antibiotics. Though effective, this method increases both the prevalence and risk of generating antibiotic-resistant bacteria. To facilitate vaccine design and/or alternative treatment efforts, a deeper understanding of the teleost immune system is essential. This review highlights the current state of teleost antibacterial immunity in the context of salmonid aquaculture. Additionally, the success of current techniques/methods used to combat bacterial diseases in salmonid aquaculture will be addressed. Filling the immunology knowledge gaps highlighted here will assist in reducing aquaculture losses in the future.

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.

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.

A Genome-Wide Association Study Reveals That Genes with Functions for Bone Development Are Associated with Body Conformation in Catfish

Body conformation is of great scientific and commercial interest for aquaculture fish species because it affects biological adaptation of the organism to environments, and is of economic importance to the aquaculture industry considering its direct effect on fillet yield. Catfish is the primary aquaculture species in the USA. Two major species used in the aquaculture industry, channel catfish and blue catfish, differ in body shape and therefore the backcross progenies serve as a good model for quantitative trait locus (QTL) analysis. Here, a genome-wide association study (GWAS) with hybrid catfish was conducted to identify the QTL for body conformation, including deheaded body length (DBL), body length (BL), body depth (BD), and body breadth (BB), which were all standardized by cubic root of body weight. Overall, the results indicate that the traits are polygenic. For DBL, linkage group (LG) 2 and LG 24 contain significant QTL, and LG 13 and LG 26 contain suggestively associated QTL (-log₁₀(P value) > 4.5). Compared with DBL, additional SNPs were identified to be associated with body length on LG 2, LG 7, and LG 18. Although no significant QTL for body depth was found, three suggestively associated QTLs were identified on LG 5, LG 13, and LG 14. No SNP for body breadth reached the threshold for suggestive association. Genes close to the associated SNPs were determined, many of which are known to be involved in bone development. This work therefore provides the basis for future identification of causal genes for the control of body conformation.

Similar Genetic Architecture with Shared and Unique Quantitative Trait Loci for Bacterial Cold Water Disease Resistance in Two Rainbow Trout Breeding Populations

Bacterial cold water disease (BCWD) causes significant mortality and economic losses in salmonid aquaculture. In previous studies, we identified moderate-large effect quantitative trait loci (QTL) for BCWD resistance in rainbow trout (Oncorhynchus mykiss). However, the recent availability of a 57 K SNP array and a reference genome assembly have enabled us to conduct genome-wide association studies (GWAS) that overcome several experimental limitations from our previous work. In the current study, we conducted GWAS for BCWD resistance in two rainbow trout breeding populations using two genotyping platforms, the 57 K Affymetrix SNP array and restriction-associated DNA (RAD) sequencing. Overall, we identified 14 moderate-large effect QTL that explained up to 60.8% of the genetic variance in one of the two populations and 27.7% in the other. Four of these QTL were found in both populations explaining a substantial proportion of the variance, although major differences were also detected between the two populations. Our results confirm that BCWD resistance is controlled by the oligogenic inheritance of few moderate-large effect loci and a large-unknown number of loci each having a small effect on BCWD resistance. We detected differences in QTL number and genome location between two GWAS models (weighted single-step GBLUP and Bayes B), which highlights the utility of using different models to uncover QTL. The RAD-SNPs detected a greater number of QTL than the 57 K SNP array in one population, suggesting that the RAD-SNPs may uncover polymorphisms that are more unique and informative for the specific population in which they were discovered.

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.

A nonclassical MHC class I U lineage locus in zebrafish with a null haplotypic variant

Three sequence lineages of MHC class I genes have been described in zebrafish (Danio rerio): U, Z, and L. The U lineage genes encoded on zebrafish chromosome 19 are predicted to provide the classical function of antigen presentation. This MHC class I locus displays significant haplotypic variation and is the only MHC class I locus in zebrafish that shares conserved synteny with the core mammalian MHC. Here, we describe two MHC class I U lineage genes, mhc1ula and mhc1uma, that map to chromosome 22. Unlike the U lineage proteins encoded on chromosome 19, Ula and Uma likely play a nonclassical role as they lack conservation of key peptide binding residues, display limited polymorphic variation, and exhibit tissue-specific expression. We also describe a null haplotype at this chromosome 22 locus in which the mhc1ula and mhc1uma genes are absent due to a ~30 kb deletion with no other MHC class I sequences present. Functional and non-functional transcripts of mhc1ula and mhc1uma were identified; however, mhc1uma transcripts were often not amplified or amplified at low levels from individuals possessing an apparently bona fide gene. These distinct U lineage genes may be restricted to the superorder Ostariophysi as similar sequences only could be identified from the blind cavefish (Astyanax mexicanus), fathead minnow (Pimephales promelas), goldfish (Carassius auratus), and grass carp (Ctenopharyngodon idella).

Influence of immune-relevant genotype on the reproductive success of a salmonid alternative mating strategy

Major histocompatibility complex (MHC) and immune-relevant gene markers were used to evaluate differences in reproductive success (RS) among naturally spawning coho salmon Oncorhynchus kisutch mate pairs involving an alternative male reproductive phenotype, known as jacks. These mate pairs included both hatchery-reared and wild origin fish such that three classes were evaluated in two consecutive years (2005 and 2006) using a previously constructed multigenerational genetic pedigree: wild × wild (W × W), hatchery × hatchery (H × H) and wild × hatchery (W × H). Oncorhynchus kisutch jack mate pairs mated randomly based on immune-relevant genotype in both years; a result consistent with the opportunistic mating strategy of jacks. An association between greater number of alleles shared at three immune-relevant gene markers and increased RS was found for: W × H mate pairs in 2005 (BHMS429), W × H pairs in 2006 (SsalR016TKU) and W × W pairs in 2006 (OMM3085). No correlation between immune gene diversity and RS was found for H × H pairs in either year. The results suggest that the influence of immune-relevant genotype on mating success may be different for jacks when compared with previous studies of large adult male O. kisutch.

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.

A genome scan for quantitative trait loci associated with Vibrio anguillarum infection resistance in Japanese flounder (Paralichthys olivaceus) by bulked segregant analysis

A recent genetic linkage map was employed to detect quantitative trait loci (QTLs) associated with Vibrio anguillarum resistance in Japanese flounder. An F1 family established and challenged with V. anguillarum in 2009 was used for QTL mapping. Of the 221 simple sequence repeat (SSR) markers used to detect polymorphisms in the parents of F1, 170 were confirmed to be polymorphic. The average distance between the markers was 10.6 cM. Equal amounts of genomic DNA from 15 fry that died early and from 15 survivors were pooled separately to constitute susceptible bulk and resistance bulk DNA. Bulked segregant analysis and QTL mapping were combined to detect candidate SSR markers and regions associated with the disease. A genome scan identified four polymorphic SSR markers, two of which were significantly different between susceptible and resistance bulk (P=0.008). These two markers were located in linkage group (LG) 7; therefore, all the SSR markers in LG7 were genotyped in all the challenged fry by single marker analysis. Using two different models, 11-17 SSR markers were detected with different levels of significance. To confirm the associations of these markers with the disease, composite interval mapping was employed to genotype all the challenged individuals. One and three QTLs, which explained more than 60 % of the phenotypic variance, were detected by the two models. Two of the QTLs were located at 48.6 cM. The common QTL may therefore be a major candidate region for disease resistance against V. anguillarum infection.

Teleost T and NK cell immunity

The main function of the immune system is to maintain the organism's homeostasis when invaded by foreign material or organisms. Prior to successful elimination of the invader it is crucial to distinguish self from non-self. Most pathogens and altered cells can be recognized by immune cells through expressed pathogen- or danger-associated molecular patterns (PAMPS or DAMPS, respectively), through non-self (e.g. allogenic or xenogenic cells) or missing major histocompatibility (MHC) class I molecules (some virus-infected target cells), and by presenting foreign non-self peptides of intracellular (through MHC class I-e.g. virus-infected target cells) or extracellular (through MHC class II-e.g. from bacteria) origin. In order to eliminate invaders directly or by destroying their ability to replicate (e.g. virus-infected cells) specialized immune cells of the innate and adaptive responses appeared during evolution. The first line of defence is represented by the evolutionarily ancient macrophages and natural killer (NK) cells. These innate mechanisms are well developed in bony fish. Two types of NK cell homologues have been described in fish: non-specific cytotoxic cells and NK-like cells. Adaptive cell-mediated cytotoxicity (CMC) requires key molecules expressed on cytotoxic T lymphocytes (CTLs) and target cells. CTLs kill host cells harbouring intracellular pathogens by binding of their T cell receptor (TCR) and its co-receptor CD8 to a complex of MHC class I and bound peptide on the infected host cell. Alternatively, extracellular antigens are taken up by professional antigen presenting cells such as macrophages, dendritic cells and B cells to process those antigens and present the resulting peptides in association with MHC class II to CD4(+) T helper cells. During recent years, genes encoding MHC class I and II, TCR and its co-receptors CD8 and CD4 have been cloned in several fish species and antibodies have been developed to study protein expression in morphological and functional contexts. Functional assays for innate and adaptive lymphocyte responses have been developed in only a few fish species. This review summarizes and discusses recent results and developments in the field of T and NK cell responses with focus on economically important and experimental model fish species in the context of vaccination.

Three vibrio-resistance related EST-SSR markers revealed by selective genotyping in the clam Meretrix meretrix

The clam Meretrix meretrix is an important commercial bivalve distributed in the coastal areas of South and Southeast Asia. In this study, marker-trait association analyses were performed based on the stock materials of M. meretrix with different vibrio-resistance profile obtained by selective breeding. Forty-eight EST-SSR markers were screened and 27 polymorphic SSRs of them were genotyped in the clam stocks with different resistance to Vibrio parahaemolyticus (11-R and 11-S) and to Vibrio harveyi (09-R and 09-C). Allele frequency distributions of the SSRs among different stocks were compared using Pearson's Chi-square test, and three functional EST-SSR markers (MM959, MM4765 and MM8364) were found to be associated with vibrio-resistance trait. The 140-bp allele of MM959 and 128-bp allele of MM4765 had significantly higher frequencies in resistant groups (11-R and 09-R) than in susceptive/control groups (11-S and 09-C) (P < 0.01), which suggested that the clams carrying these two alleles have stronger resistance against vibrio. Clam individuals of 11-S were divided into three subgroups based on the survival time post-challenge, and the multi-dimensional scaling (MDS) analysis showed that clusters generated by genetic similarity revealed by the three SSR markers were consistent with the three subgroups distinctions. The putative functions of contig959, contig4765 and contig8364 also suggested that the three SSR-involved genes might play important roles in immunity of M. meretrix. All these results supported that EST-SSR markers MM959, MM4765 and MM8364 were associated with vibrio-resistance and would be useful for marker-assisted selection (MAS) in M. meretrix genetic breeding.

Mapping rainbow trout immune genes involved in inflammation reveals conserved blocks of immune genes in teleosts

We report the genetic map location of 14 genes involved in the inflammatory response to salmonid bacterial and viral pathogens, which brings the total number of immune genes mapped in rainbow trout (RT, Oncorhynchus mykiss) to 61. These genes were mapped as candidate genes that may be involved in resistance to bacterial kidney disease, as well as candidates for known QTL for resistance to infectious hematopoietic necrosis virus, infectious pancreatic necrosis virus and Ceratomyxa shasta. These QTL map to one or more of the linkage groups containing immune genes. The combined analysis of these linkage results and those of previously mapped immune genes in RT shows that many immune genes are found in syntenic blocks of genes that have been retained in teleosts despite species divergence and genome duplication events.

Linkage mapping of toll-like receptors (TLRs) in Japanese flounder, Paralichthys olivaceus

Toll-like receptors (TLRs) are responsible for the recognition of specific pathogen-associated molecular patterns and consequently activate signal pathways leading to inflammatory and interferon responses. The region surrounding several TLRs was previously found to be associated with resistance to specific disease. Hence, we determined the location of 11 TLRs in Japanese flounder (Paralichthys olivaceus) using polymorphic microsatellite markers. TLR1 and TLR3 were located on linkage group (LG) 21 and 7, respectively. Membrane TLR5 and soluble TLR5 were mapped to LG22. TLR7 and TLR8 were mapped to LG3. TLR9 was found on LG1 and TLR14 and TLR21 were located on the same linkage group, LG10. TLR22 was found on LG8. Interestingly, TLR2 was mapped with the previously reported Poli9-8TUF microsatellite marker which is tightly associated with lymphocystis virus disease resistance. Therefore, TLR2 is a candidate gene for resistance to lymphocystis disease. These results imply that the location of a TLR associated with a particular disease may be valuable for the research on the relationship between host immune response and disease resistance.

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.

Genetic linkage mapping and analysis of muscle fiber-related QTLs in common carp (Cyprinus carpio L.)

A genetic linkage map of common carp (Cyprinus carpio L.) was constructed using Type I and Type II microsatellite markers and a pseudo-testcross mapping strategy. The microsatellite markers were isolated from microsatellite-enriched genomic libraries and tested for their segregation in a full-sib mapping panel containing 92 individuals. A total of 161 microsatellite loci were mapped into 54 linkage groups. The total lengths of the female, male and consensus maps were 2,000, 946, and 1,852 cM, with an average marker spacing of approximately 13, 7, and 11 cM, respectively. Muscle fiber-related traits, including muscle fiber cross-section area and muscle fiber density, were mapped to the genetic map. Three QTLs for muscle fiber cross-section area and two QTLs for muscle fiber density were identified when considering both significant and suggestive QTL effects. The QTLs with largest effects for muscle fiber cross-section area and muscle fiber density were 21.9% and 18.9%, and they were located in LG3, respectively.

Winning the invasion roulette: escapes from fish farms increase admixture and facilitate establishment of non-native rainbow trout

Aquaculture is a major source of invasive aquatic species, despite the fact that cultured organisms often have low genetic diversity and tend to be maladapted to survive in the wild. Yet, to what extent aquaculture escapees become established by means of high propagule pressure and multiple origins is not clear. We analysed the genetic diversity of 15 established populations and four farmed stocks of non-native rainbow trout in Chile, a species first introduced for recreational fishing around 1900, but which has in recent decades escaped in large numbers from fish farms and become widespread. Aquaculture propagule pressure was a good predictor of the incidence of farm escapees, which represented 16% of all free-ranging rainbow trout and were present in 80% of the study rivers. Hybrids between farm escapes and established trout were present in all rivers at frequencies ranging between 7 and 69%, and population admixture was positively correlated with genetic diversity. We suggest that non-native salmonids introduced into the Southern Hemisphere could benefit from admixture because local adaptations may not have yet developed, and there may be initially little fitness loss resulting from outbreeding depression.

A first generation integrated map of the rainbow trout genome

Background

Rainbow trout (Oncorhynchus mykiss) are the most-widely cultivated cold freshwater fish in the world and an important model species for many research areas. Coupling great interest in this species as a research model with the need for genetic improvement of aquaculture production efficiency traits justifies the continued development of genomics research resources. Many quantitative trait loci (QTL) have been identified for production and life-history traits in rainbow trout. An integrated physical and genetic map is needed to facilitate fine mapping of QTL and the selection of positional candidate genes for incorporation in marker-assisted selection (MAS) programs for improving rainbow trout aquaculture production.

Results

The first generation integrated map of the rainbow trout genome is composed of 238 BAC contigs anchored to chromosomes of the genetic map. It covers more than 10% of the genome across segments from all 29 chromosomes. Anchoring of 203 contigs to chromosomes of the National Center for Cool and Cold Water Aquaculture (NCCCWA) genetic map was achieved through mapping of 288 genetic markers derived from BAC end sequences (BES), screening of the BAC library with previously mapped markers and matching of SNPs with BES reads. In addition, 35 contigs were anchored to linkage groups of the INRA (French National Institute of Agricultural Research) genetic map through markers that were not informative for linkage analysis in the NCCCWA mapping panel. The ratio of physical to genetic linkage distances varied substantially among chromosomes and BAC contigs with an average of 3,033 Kb/cM.

Conclusions

The integrated map described here provides a framework for a robust composite genome map for rainbow trout. This resource is needed for genomic analyses in this research model and economically important species and will facilitate comparative genome mapping with other salmonids and with model fish species. This resource will also facilitate efforts to assemble a whole-genome reference sequence for rainbow trout.

Quantitative trait loci for resistance to fish pasteurellosis in gilthead sea bream (Sparus aurata)

Fish pasteurellosis is a bacterial disease causing important losses in farmed fish, including gilthead sea bream, a teleost fish of great relevance in marine aquaculture. We report in this study a QTL analysis for resistance to fish pasteurellosis in this species. An experimental population of 500 offspring originating from eight sires and six dams in a single mass-spawning event was subjected to a disease challenge with Photobacterium damselae subsp. piscicida (Phdp), the causative agent of fish pasteurellosis. A total of 151 microsatellite loci were genotyped in the experimental population, and half-sib regression QTL analysis was carried out on two continuous traits, body length at time of death and survival, and for two binary traits, survival at day 7 and survival at day 15, when the highest peaks of mortality were observed. Two significant QTLs were detected for disease resistance. The first one was located on linkage group LG3 affecting late survival (survival at day 15). The second one, for overall survival, was located on LG21, which allowed us to highlight a potential marker (Id13) linked to disease resistance. A significant QTL was also found for body length at death on LG6 explaining 5-8% of the phenotypic variation.

Relationships between growth and disease resistance in rainbow trout, Oncorhynchus mykiss (Walbaum)

Rainbow trout from 23 families were evaluated for growth and resistance to the bacterial coldwater disease (BCWD) caused by Flavobacterium psychrophilum and infectious haematopoietic necrosis (IHN) caused by IHN virus. Average family weights were between 161 and 263 g with an average of 225 g at 213 days post-fertilization with specific growth rates ranging from 2.37 to 2.88. Per cent survival of fish challenged with F. psychrophilum was between 18% and 100%, while for those challenged with IHNV, the range was between 12% and 93%. Significant positive correlations were found for end body weight and resistance to IHN (P < 0.05) and for early body weight and resistance to BCWD (P < 0.1). However, no significant correlations were detected between resistance to both pathogens or disease resistance and overall genetic diversity or diversity within the major histocompatibility locus.

Construction of a genetic linkage map and mapping of a female-specific DNA marker in half-smooth tongue sole (Cynoglossus semilaevis)

The half-smooth tongue sole (Cynoglossus semilaevis, hereafter, "tongue sole") is a marine flatfish with great commercial importance for fisheries and aquaculture in China. It has also been a promising model for the study of sex determination mechanisms in fish. Here, we report the construction of a genetic linkage map for the tongue sole, based on 137 markers including 103 AFLP markers, 33 microsatellite markers, and one female-specific DNA marker. Twenty-six linkage groups (LGs) were found. The total map length was 934.6 cM (Kosambi), with an average spacing of 8.4 cM, covering 64.4% of the estimated genome size. Furthermore, a female-specific SCAR marker, CseF-382, was mapped on LG5. This study represents the first genetic linkage map in the tongue sole. This map has great potential in the identification of quantitative traits loci and sex-related genes and marker-assisted selection in the tongue sole. Meanwhile, the new set of polymorphic microsatellite markers developed in this study is not only useful for genetic mapping but also of critical importance for studies on genetic diversity and broodstock management in tongue sole.

Major histocompatibility complex heterozygote advantage and widespread bacterial infections in populations of Chinook salmon (Oncorhynchus tshawytscha)

Despite growing evidence for parasite-mediated selection on the vertebrate major histocompatibility complex (MHC), little is known about variation in the bacterial parasite community within and among host populations or its influence on MHC evolution. In this study, we characterize variation in the parasitic bacterial community associated with Chinook salmon (Oncorhynchus tshawytscha) fry in five populations in British Columbia (BC), Canada across 2 years, and examine whether bacterial infections are a potential source of selection on the MHC. We found an unprecedented diversity of bacteria infecting fry with a total of 55 unique bacteria identified. Bacterial infection rates varied from 9% to 29% among populations and there was a significant isolation by distance relationship in bacterial community phylogenetic similarity across the populations. Spatial variation in the frequency of infections and in the phylogenetic similarity of bacterial communities may result in differential parasite-mediated selection at the MHC across populations. Across all populations, we found evidence of a heterozygote advantage at the MHC class II, which may be a source of balancing selection on this locus. Interestingly, a co-inertia analysis indicated only susceptibility associations between a few of the MHC class I and II alleles and specific bacterial parasites; there was no evidence that any of the alleles provided resistance to the bacteria. Our results reveal a complex bacterial community infecting populations of a fish and underscore the importance of considering the role of multiple pathogens in the evolution of host adaptations.

A first generation BAC-based physical map of the rainbow trout genome

Background

Rainbow trout (Oncorhynchus mykiss) are the most-widely cultivated cold freshwater fish in the world and an important model species for many research areas. Coupling great interest in this species as a research model with the need for genetic improvement of aquaculture production efficiency traits justifies the continued development of genomics research resources. Many quantitative trait loci (QTL) have been identified for production and life-history traits in rainbow trout. A bacterial artificial chromosome (BAC) physical map is needed to facilitate fine mapping of QTL and the selection of positional candidate genes for incorporation in marker-assisted selection (MAS) for improving rainbow trout aquaculture production. This resource will also facilitate efforts to obtain and assemble a whole-genome reference sequence for this species.

Results

The physical map was constructed from DNA fingerprinting of 192,096 BAC clones using the 4-color high-information content fingerprinting (HICF) method. The clones were assembled into physical map contigs using the finger-printing contig (FPC) program. The map is composed of 4,173 contigs and 9,379 singletons. The total number of unique fingerprinting fragments (consensus bands) in contigs is 1,185,157, which corresponds to an estimated physical length of 2.0 Gb. The map assembly was validated by 1) comparison with probe hybridization results and agarose gel fingerprinting contigs; and 2) anchoring large contigs to the microsatellite-based genetic linkage map.

Conclusion

The production and validation of the first BAC physical map of the rainbow trout genome is described in this paper. We are currently integrating this map with the NCCCWA genetic map using more than 200 microsatellites isolated from BAC end sequences and by identifying BACs that harbor more than 300 previously mapped markers. The availability of an integrated physical and genetic map will enable detailed comparative genome analyses, fine mapping of QTL, positional cloning, selection of positional candidate genes for economically important traits and the incorporation of MAS into rainbow trout breeding programs.

MHC genetic structure and divergence across populations of Chinook salmon (Oncorhynchus tshawytscha)

The major histocompatibility complex (MHC) is thought to be under strong selection pressure because of its integral role in pathogen recognition. Consequently, patterns of MHC genetic variation should reflect selection pressures across the landscape. We examined genetic variation and population genetic structure at the MHC class I-A1 and class II-B1 exons in five Chinook salmon (Oncorhynchus tshawytscha) populations from two geographic regions in British Columbia, Canada. We then compared estimates of population structure at the MHC genes with neutral estimates based on microsatellites to examine the potential for local adaptation at the MHC. Chinook salmon are in decline throughout much of their native range and understanding the degree of local adaptation exhibited by the MHC may be important in conservation planning. Comparisons among populations yielded higher G'(ST) estimates for the MHC class I than expected under neutrality based on the microsatellites. In contrast, the MHC class II tended to exhibit lower G'(ST) values than did the microsatellites. These results suggest that across populations unique selection pressures are driving allele frequency differences at the MHC class I but that the MHC class II may be the subject of homogenizing selection. Rates of nonsynonymous versus synonymous substitutions found in codons associated within the MHC class I and II peptide-binding regions provided strong evidence of positive selection. Together, these results support the hypothesis that selection is influencing genetic variation at the MHC, but suggest that selection pressures may vary at the two classes of loci both at the sequence and population levels.

Resistance of common carp (Cyprinus carpio L.) to Cyprinid herpesvirus-3 is influenced by major histocompatibility (MH) class II B gene polymorphism

The role of MH class II B (Cyca-DAB1-like) genes in resistance of common carp (Cyprinus carpio L.) to Cyprinid herpesvirus-3 (CyHV-3), also known as koi herpesvirus (KHV) was analysed. The material consisted of 934 fish from six carp crosses. Fish were challenged with CyHV-3 at an age of 7 and 10 months. During challenge experiments the peak of mortality caused by CyHV-3 was observed at days 8-12 p.i. and the overall cumulative mortality reached 79.9%. Among six Cyca-DAB1-like genotypes, revealed by PCR-RF-SSCP analysis, one genotype (E) was found associated with higher resistance to CyHV-3. Three other genotypes (B, H and J) could be linked to higher susceptibility to CyHV-3. Analysis of the alleles that compose the Cyca-DAB1-like genotypes linked one particular allele (Cyca-DAB1*05) to significantly increased, and two alleles (Cyca-DAB1*02 and Cyca-DAB1*06) to significantly decreased resistance to CyHV-3. Our data indicate that MH class II B genes could be used as potential genetic markers in breeding of common carp for resistance to this virus.

Rainbow trout resistance to bacterial cold-water disease is moderately heritable and is not adversely correlated with growth

The objectives of this study were to estimate the heritabilities for and genetic correlations among resistance to bacterial cold-water disease and growth traits in a population of rainbow trout (Oncorhynchus mykiss). Bacterial cold-water disease, a chronic disease of rainbow trout, is caused by Flavobacterium psychrophilum. This bacterium also causes acute losses in young fish, known as rainbow trout fry syndrome. Selective breeding for increased disease resistance is a promising strategy that has not been widely used in aquaculture. At the same time, improving growth performance is critical for efficient production. At the National Center for Cool and Cold Water Aquaculture, reducing the negative impact of diseases on rainbow trout culture and improving growth performance are primary objectives. In 2005, when fish averaged 2.4 g, 71 full-sib families were challenged with F. psychrophilum and evaluated for 21 d. Overall survival was 29.3% and family rates of survival varied from 1.5 to 72.5%. Heritability of postchallenge survival, an indicator of disease resistance, was estimated to be 0.35 +/- 0.09. Body weights at 9 and 12 mo posthatch and growth rate from 9 to 12 mo were evaluated on siblings of the fish in the disease challenge study. Growth traits were moderately heritable, from 0.32 for growth rate to 0.61 for 12-mo BW. Genetic and phenotypic correlations between growth traits and resistance to bacterial cold-water disease were not different from zero. These results suggest that genetic improvement can be made simultaneously for growth and bacterial cold-water disease resistance in rainbow trout by using selective breeding.

A second generation genetic map for rainbow trout (Oncorhynchus mykiss)

Background

Genetic maps characterizing the inheritance patterns of traits and markers have been developed for a wide range of species and used to study questions in biomedicine, agriculture, ecology and evolutionary biology. The status of rainbow trout genetic maps has progressed significantly over the last decade due to interest in this species in aquaculture and sport fisheries, and as a model research organism for studies related to carcinogenesis, toxicology, comparative immunology, disease ecology, physiology and nutrition. We constructed a second generation genetic map for rainbow trout using microsatellite markers to facilitate the identification of quantitative trait loci for traits affecting aquaculture production efficiency and the extraction of comparative information from the genome sequences of model fish species.

Results

A genetic map ordering 1124 microsatellite loci spanning a sex-averaged distance of 2927.10 cM (Kosambi) and having 2.6 cM resolution was constructed by genotyping 10 parents and 150 offspring from the National Center for Cool and Cold Water Aquaculture (NCCCWA) reference family mapping panel. Microsatellite markers, representing pairs of loci resulting from an evolutionarily recent whole genome duplication event, identified 180 duplicated regions within the rainbow trout genome. Microsatellites associated with genes through expressed sequence tags or bacterial artificial chromosomes produced comparative assignments with tetraodon, zebrafish, fugu, and medaka resulting in assignments of homology for 199 loci.

Conclusion

The second generation NCCCWA genetic map provides an increased microsatellite marker density and quantifies differences in recombination rate between the sexes in outbred populations. It has the potential to integrate with cytogenetic and other physical maps, identifying paralogous regions of the rainbow trout genome arising from the evolutionarily recent genome duplication event, and anchoring a comparative map with the zebrafish, medaka, tetraodon, and fugu genomes. This resource will facilitate the identification of genes affecting traits of interest through fine mapping and positional cloning of candidate genes.

Phenotypic and genotypic analysis of Flavobacterium psychrophilum isolates from Ontario salmonids with bacterial coldwater disease

Flavobacterium psychrophilum is the causative agent of bacterial coldwater disease (BCWD) and rainbow trout fry syndrome. BCWD has a considerable economic impact on aquaculture operations in Ontario, Canada, and our limited understanding of the population structure and epidemiology of F. psychrophilum isolates is an impediment to the development of improved management strategies. Seventy-five 16S rRNA gene and gyr polymerase chain reaction positive isolates of F. psychrophilum that had been collected over a 16-year period from farmed salmonids with tail rot, necrotic myositis, and osteochondrosis were characterized morphologically, biochemically, and genotypically. Although the isolates were homogeneous by preliminary biochemical and phenotypic characterization, two distinct biovars were found by API ZYM testing. As well, four restriction pattern types were detected by 16S rRNA polymerase chain reaction – restriction fragment length polymorphism analysis and there was a significant (P < 0.001) correlation between biovar I and digestion with MaeIII and between biovar II and digestion with MnlI or no site (P < 0.05). Further heterogenity was detected by sequence analysis of a 194 bp stem loop 3 region of rRNA. Nine sequence types were identified; 40/46 biovar I isolates were sequence type “a”, while 21/32 biovar II isolates belonged to either sequence type “c” or “d”. More than one biovar and genotype was identified among the strains recovered from separate fish sampled from three groups of rainbow trout ( Oncorhynchus mykiss ) experiencing BCWD mortality events. No association was found between genotype or biovar and type of disease. Taken together, these data suggest that F. psychrophilum from Ontario can be grouped into two major lineages based on biovar and 16S rRNA polymorphisms, and although three major strain types were most frequently isolated in this study, it appears that the population of F. psychrophilum with pathogenic potential is quite heterogeneous.