MHC polymorphism and disease resistance in Atlantic salmon (Salmo salar); facing pathogens with single expressed major histocompatibility class I and class II loci

The dominantly expressed class I molecule of the chicken MHC is explained by coevolution with the polymorphic peptide transporter (TAP) genes

In most mammals, the MHC class I molecules are polymorphic and determine the specificity of peptide presentation, whereas the transporter associated with antigen presentation (TAP) heterodimers are functionally monomorphic. In chickens, there are two classical class I genes but only one is expressed at a high level, which can result in strong MHC associations with resistance to particular infectious pathogens. However, the basis for having a single dominantly expressed class I molecule has been unclear. Here we report TAP1 and TAP2 sequences from 16 chicken lines, and show that both genes have high allelic polymorphism and moderate sequence diversity, with variation in positions expected for peptide binding. We analyze peptide translocation in two MHC haplotypes, showing that chicken TAPs specify translocation at three peptide positions, matching the peptide motif of the single dominantly expressed class I molecule. These results show that coevolution between class I and TAP genes can explain the presence of a single dominantly expressed class I molecule in common chicken MHC haplotypes. Moreover, such coevolution in the primordial MHC may have been responsible for the appearance of the antigen presentation pathways at the birth of the adaptive immune system.

Detecting population structure in a high gene-flow species, Atlantic herring (Clupea harengus): direct, simultaneous evaluation of neutral vs putatively selected loci

In many marine fish species, genetic population structure is typically weak because populations are large, evolutionarily young and have a high potential for gene flow. We tested whether genetic markers influenced by natural selection are more efficient than the presumed neutral genetic markers to detect population structure in Atlantic herring (Clupea harengus), a migratory pelagic species with large effective population sizes. We compared the spatial and temporal patterns of divergence and statistical power of three traditional genetic marker types, microsatellites, allozymes and mitochondrial DNA, with one microsatellite locus, Cpa112, previously shown to be influenced by divergent selection associated with salinity, and one locus located in the major histocompatibility complex class IIA (MHC-IIA) gene, using the same individuals across analyses. Samples were collected in 2002 and 2003 at two locations in the North Sea, one location in the Skagerrak and one location in the low-saline Baltic Sea. Levels of divergence for putatively neutral markers were generally low, with the exception of single outlier locus/sample combinations; microsatellites were the most statistically powerful markers under neutral expectations. We found no evidence of selection acting on the MHC locus. Cpa112, however, was highly divergent in the Baltic samples. Simulations addressing the statistical power for detecting population divergence showed that when using Cpa112 alone, compared with using eight presumed neutral microsatellite loci, sample sizes could be reduced by up to a tenth while still retaining high statistical power. Our results show that the loci influenced by selection can serve as powerful markers for detecting population structure in high gene-flow marine fish species.

Contrasting responses to selection in class I and class IIα major histocompatibility-linked markers in salmon

Comparison of levels and patterns of genetic variation in natural populations either across loci or against neutral expectation can yield insight into locus-specific differences in the strength and direction of evolutionary forces. We used both approaches to test the hypotheses on patterns of selection on major histocompatibility (MH)-linked markers. We performed temporal analyses of class I and class IIα MH-linked markers and eight microsatellite loci in two Atlantic salmon populations in Ireland on two temporal scales: over six decades and 9 years in the rivers Burrishoole and Delphi, respectively. We also compared contemporary Burrishoole and Delphi samples with nearby populations for the same loci. On comparing patterns of temporal and spatial differentiation among classes of loci, the class IIα MH-linked marker was consistently identified as an outlier compared with patterns at the other microsatellite loci or neutral expectation. We found higher levels of temporal and spatial heterogeneity in heterozygosity (but not in allelic richness) for the class IIα MH-linked marker compared with microsatellites. Tests on both within- and among-population differentiation are consistent with directional selection acting on the class IIα-linked marker in both temporal and spatial comparisons, but only in temporal comparisons for the class I-linked marker. Our results indicate a complex pattern of selection on MH-linked markers in natural populations of Atlantic salmon. These findings highlight the importance of considering selection on MH-linked markers when using these markers for management and conservation purposes.

Allelic polymorphism, gene duplication and balancing selection of the MHC class II DAB gene of Cynoglossus semilaevis (Cynoglossidae)

Major histocompatibility complex (MHC) genes play an important role in the immune response of vertebrates. Allelic polymorphism and evolutionary mechanism of MHC genes have been investigated in many mammals, but much less is known in teleosts. We examined the polymorphism, gene duplication and balancing selection of the MHC class II DAB gene of the half-smooth tongue sole (Cynoglossus semilaevis); 23 alleles were found in this species. Gene duplication manifested as three to six distinct sequences at each domain in the same individuals. Non-synonymous substitutions occurred at a significantly higher frequency than synonymous substitutions in the PBR domain, suggesting balancing selection for maintaining polymorphisms at the MHC II DAB locus. Many positive selection sites were found to act very intensely on antigen-binding sites of MHC class II DAB gene.

Geographic variation of the major histocompatibility complex in Eastern Atlantic grey seals (Halichoerus grypus)

Pathogen-driven balancing selection maintains high genetic diversity in many vertebrates, particularly in the major histocompatibility complex (MHC) immune system gene family, which is often associated with disease susceptibility. In large natural populations where subpopulations face different pathogen pressures, the MHC should show greater genetic differentiation within a species than neutral markers. We examined genetic diversity at the MHC-DQB locus and nine putatively neutral microsatellite markers in grey seals (Halichoerus grypus) from eight United Kingdom (UK) colonies, the Faeroe Islands and Sable Island, Canada. Five DQB alleles were identified in grey seals, which varied in prevalence across the grey seal range. Among the seal colonies, significant differences in DQB allele and haplotype frequencies and in average DQB heterozygosity were observed. Additionally, the DQB gene exhibited greater differentiation among colonies compared with neutral markers, yet a weaker pattern of isolation by distance (IBD). After correcting for the underlying IBD pattern, subpopulations breeding in similar habitats were more similar to one another in DQB allele frequencies than populations breeding in different habitats, but the same did not hold true for microsatellites, suggesting that habitat-specific pathogen pressure influences MHC evolution. Overall, the data are consistent with selection at MHC-DQB loci in grey seals with both varying selective pressures and geographic population structure appearing to influence the DQB genetic composition of breeding colonies.

Effects of different mating scenarios on embryo viability in brown trout

Mating with attractive or dominant males is often predicted to offer indirect genetic benefits to females, but it is still largely unclear how important such non-random mating can be with regard to embryo viability. We sampled a natural population of adult migratory brown trout (Salmo trutta), bred them in vitro in a half-sib breeding design to separate genetic from maternal environmental effects, raised 2098 embryos singly until hatching, and exposed them experimentally to different levels of pathogen stress at a late embryonic stage. We found that the embryos' tolerance to the induced pathogen stress was linked to the major histocompatibility complex (MHC) of their parents, i.e. certain MHC genotypes appeared to provide better protection against infection than others. We also found significant additive genetic variance for stress tolerance. Melanin-based dark skin patterns revealed males with 'good genes', i.e. embryos fathered by dark coloured males had a high tolerance to infection. Mating with large and dominant males would, however, not improve embryo viability when compared to random mating. We used simulations to provide estimates of how mate choice based on MHC or melanin-based skin patterns would influence embryos' tolerance to the experimentally induced pathogen stress.

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.

MHC class IIalpha gene polymorphism and its association with resistance/susceptibility to Vibrio anguillarum in Japanese flounder (Paralichthys olivaceus)

Association between polymorphism of the major histocompatibility complex (MHC) genes and disease resistance has been documented for few teleosts. In this study, we first investigate the genetic variation at the MHC IIalpha gene in Japanese flounder (Paralichthys olivaceus) for survival after challenge with bacterial infection. To explore the specific allele associated with disease resistance, about 6000 individuals from 60 families challenged with Vibrio anguillarum, which causes significantly different mortality in flounder families. 15-20 individuals from each of six high-resistance (HR) and six low-resistance (LR) families were screened for their MHC class IIalpha genotypes using sequence analysis. High polymorphism of MHC IIalpha gene and at least two loci were discovered in Japanese flounder and the rate of d(N) occurred at a significantly higher frequency than that of d(S) in PBR and non-PBR, the balancing selection hypothesis could explain the high polymorphism of MHC IIalpha gene in flounder. From the 212 individuals, a total of 55 MHC class IIalpha alleles were identified, and 9 alleles were used to study association between alleles and resistance to disease. Three alleles, Paol-DAA*1301, Paol-DAA*1401 and Paol-DAA*2201 were significantly associated with resistance against V. anguillarum, and Paol-DAA*1001 and Paol-DAA*1501 alleles were significantly associated with increased susceptibility to V. anguillarum. This study confirmed the association between alleles of MHC class IIalpha gene and disease resistance or susceptibility to bacterial infection in flounder, and the disease resistance-related MHC markers could be used for molecular marker-assisted selective breeding in the flounder.

Gene duplication and evidence for balancing selection acting on MHC class II DAA gene of the half-smooth tongue sole (Cynoglossus semilaevis)

Allelic polymorphism and evolution mechanism of major histocompatibility complex (MHC) genes has been investigated in many mammals, however, much less is known in teleost. In order to investigate the mechanisms creating and maintaining variability at the MHC class II DAA locus, we examined the polymorphism, gene duplication and balancing selection of MHC class II DAA gene of the half-smooth tongue sole (Cynoglossus semilaevis). We described 33 alleles in the C. semilaevis, recombination and gene duplication seems to play more important roles in the origin of new alleles. The rate of non-synonymous substitutions (d(N)) occurred at a significantly higher frequency than that of synonymous substitutions (d(S)) in peptide-binding region (PBR) and non-PBR, suggesting balancing selection for maintaining polymorphisms at the MHC II DAA locus. Many positive selection sites were found to act very intensively on antigen-binding sites. Our founding suggests a snapshot in an evolutionary process of MHC-DAA gene evolution of the C. semilaevis.

MHC evolution in three salmonid species: a comparison between class II alpha and beta genes

The genes of the major histocompatibility complex (MHC) are amongst the most variable in vertebrates and represent some of the best candidates to study processes of adaptive evolution. However, despite the number of studies available, most of the information on the structure and function of these genes come from studies in mammals and birds in which the MHC class I and II genes are tightly linked and class II alpha exhibits low variability in many cases. Teleost fishes are among the most primitive vertebrates with MHC and represent good organisms for the study of MHC evolution because their class I and class II loci are not physically linked, allowing for independent evolution of both classes of genes. We have compared the diversity and molecular mechanisms of evolution of classical MH class II alpha and class II beta loci in farm populations of three salmonid species: Oncorhynchus kisutch, Oncorhynchus mykiss and Salmo salar. We found single classical class II loci and high polymorphism at both class II alpha and beta genes in the three species. Mechanisms of evolution were common for both class II genes, with recombination and point mutation involved in generating diversity and positive selection acting on the peptide-binding residues. These results suggest that the maintenance of variability at the class IIalpha gene could be a mechanism to increase diversity in the MHC class II in salmonids in order to compensate for the expression of one single classical locus and to respond to a wider array of parasites.

Polymorphism of the superoxide dismutase gene family in the bay scallop (Argopecten irradians) and its association with resistance/susceptibility to Vibrio anguillarum

The superoxide dismutases (SODs) are a family of enzymes that function as the first line of antioxidant defense against highly reactive superoxide radicals. The bay scallop Argopecten irradians contains three unique superoxide dismutases: Ai-icCuZnSOD, Ai-MnSOD and Ai-ecCuZnSOD, which were characterized in our previous studies. qRT-PCR was also performed to characterize the temporal expression of SODs in the hemocytes of bay scallops injected with the bacterium Vibrio anguillarum. To characterize the SOD family in A. irradians completely, we evaluated the polymorphism in the SOD genes and investigated the association of this polymorphism with resistance/susceptibility to V. anguillarum. Fifty-nine SNPs were identified in the promoter, exon and partial intron sequences of the three SOD genes. AiECSOD contained the most SNPs, as compared to AiCuZnSOD and AiMnSOD, and the majority of these were located in the promoter. Among them, the genotypes of -1739 T-C SNP in the AiCuZnSOD promoter and alleles of the -498 A-T and -267 G-A SNPs in the AiECSOD promoter showed a significant association with resistance/susceptibility to V. angullarum (P<0.05). The only non-synonymous SNP that was identified, E1-38 C-A in Ai-ecCuZnSOD, was a dimorphism caused by a C to A transition that resulted in a Thr to Lys substitution at position 13 in the signal peptide. The Thr allele was associated with increased susceptibility to V. anguillarum (P<0.05). To confirm the presumption, another independent challenge experiment was performed, in which the cumulative mortality of Ai-icCuZnSOD Q-1739 genotype TT was significantly lower than TC (P<0.05). Ai-ecCuZnSOD Q-498 genotype AA and AT were significantly lower than TT (P<0.05), Ai-ecCuZnSOD E1-3 genotype AA was significantly higher than CA and CC (P<0.05). The results suggested that these three polymorphic loci could be potential gene markers for the future molecular selection of strains that are resistant to diseases caused by V. anguillarum.

Salmo salar and Esox lucius full-length cDNA sequences reveal changes in evolutionary pressures on a post-tetraploidization genome

Background

Salmonids are one of the most intensely studied fish, in part due to their economic and environmental importance, and in part due to a recent whole genome duplication in the common ancestor of salmonids. This duplication greatly impacts species diversification, functional specialization, and adaptation. Extensive new genomic resources have recently become available for Atlantic salmon (Salmo salar), but documentation of allelic versus duplicate reference genes remains a major uncertainty in the complete characterization of its genome and its evolution.

Results

From existing expressed sequence tag (EST) resources and three new full-length cDNA libraries, 9,057 reference quality full-length gene insert clones were identified for Atlantic salmon. A further 1,365 reference full-length clones were annotated from 29,221 northern pike (Esox lucius) ESTs. Pairwise d_N/d_S comparisons within each of 408 sets of duplicated salmon genes using northern pike as a diploid out-group show asymmetric relaxation of selection on salmon duplicates.

Conclusions

9,057 full-length reference genes were characterized in S. salar and can be used to identify alleles and gene family members. Comparisons of duplicated genes show that while purifying selection is the predominant force acting on both duplicates, consistent with retention of functionality in both copies, some relaxation of pressure on gene duplicates can be identified. In addition, there is evidence that evolution has acted asymmetrically on paralogs, allowing one of the pair to diverge at a faster rate.

Temperature and length-dependent modulation of the MH class II beta gene expression in brook charr (Salvelinus fontinalis) by a cis-acting minisatellite

It is widely recognized that the variation in gene regulation is an important factor from which evolutionary changes in diverse aspects of phenotype can be observed in all organisms. Distinctive elements with functional roles on gene regulation have been identified within the non-coding part of the genome, including repeated elements. Major histocompatibility complex (MHC) genes have been the subject of an abundant literature which made them unique candidates for studies of adaptation in natural populations. Yet, the vast majority of studies on MHC genes have dealt with patterns of polymorphism in sequence variation while very few paid attention to the possible implication of differential expression in adaptive responses. In this paper, we report the identification of a polymorphic minisatellite formed of a 32 nucleotides motif (38% G+C) involved in regulation of the major histocompatibility class II beta gene (MHII beta) of brook charr (Salvelinus fontinalis). Our main objectives were: to analyze the variability of this minisatellite found in the second intron of the MHII beta gene and to document its effect to the variation of expression level of this gene under different environmental conditions. Distinctive number of the minisatellite repeats were associated with each different MHII beta alleles identified from exon 2 sequences. Relative expression levels of specific alleles in heterozygous individuals were determined from fish lymphocytes in different genotypes. We found that alleles carrying the longest minisatellite showed a significant 1.67-2.56-fold reduction in the transcript expression relatively to the shortest one. Results obtained in three different genotypes also indicated that the repressive activity associated to the longest minisatellite was more effective at 18 degrees C compared to 6 degrees C. In contrast, no significant difference was observed in transcript levels between alleles with comparable minisatellite length at both temperatures. We also depicted a significant up-regulation of the total MHII beta transcript at 6 degrees C relative to 18 degrees C. These results reveal for the first time that a temperature-sensitive minisatellite could potentially play an important role in the gene regulation of the adaptive immune response in fishes.

Genetic architecture of rainbow trout survival from egg to adult

Survival from birth to a reproductive adult is a challenge that only robust individuals resistant to a variety of mortality factors will overcome. To assess whether survival traits share genetic architecture throughout the life cycle, we estimated genetic correlations for survival within fingerling stage, and across egg, fingerling and grow-out stages in farmed rainbow trout. Genetic parameters of survival at three life cycle stages were estimated for 249 166 individuals originating from ten year classes of a pedigreed population. Despite being an important fitness component, survival traits harboured significant but modest amount of genetic variation (h2=0.07-0.27). Weak associations between survival during egg-fry and fingerling periods, between early and late fingerling periods (rG=0.30) and generally low genetic correlations between fingerling and grow-out survival (mean rG=0.06) suggested that life-stage specific survival traits are best regarded as separate traits. However, in the sub-set of data with detailed time of death records, positive genetic correlations between early and late fingerling survival (rG=0.89) showed that during certain years the best genotypes in the early period were also among the best in the late period. That survival across fingerling period can be genetically the same, trait was indicated also by only slightly higher heritability (h2=0.15) estimated with the survival analysis of time to death during fingerling period compared to the analysis treating fingerling survival as a binary character (h2=0.11). The results imply that (1) inherited resistance against unknown mortality factors exists, but (2) ranking of genotypes changes across life stages.

Identification of a novel allele of peroxisome proliferator-activated receptor gamma (PPARG) and its association with resistance to Aeromonas salmonicida in Atlantic salmon (Salmo salar)

Bacterial and viral diseases are major problems in Atlantic salmon aquaculture, but may be challenged through selection of brood stock with enhanced survival to diseases. Today's selection strategy is based on controlled challenge tests using siblings of the breeding candidates, and is thus indirect. Direct trait records on breeding candidates can potentially be provided through identification of genetic variation linked to the susceptibility to the disease. Peroxisome proliferator-activated receptor gamma (PPARG) is a lipid-sensing transcription factor primarily known for inducing fat-accumulation in adipocytes, but also in lipid-accumulating macrophages, in mammalian species. Here we report a novel allele of PPARG, pparg-2, in Atlantic salmon. pparg-2 has an insertion of sixty nucleotides that encodes two additional copies of the almost perfect decapeptide motif, (F/C/Y)NHSPDR(S/N)HS, compared to the previously described pparg-1. pparg-1 contains six copies of this repeat unit whereas eight copies are present in the novel pparg-2 allele. pparg-2 mRNA was detectable in kidney and spleen of random Atlantic salmon samples. Here, we studied the effect of pparg-1 and pparg-2 on survival upon challenge to a highly virulent bacterium, Aeromonas salmonicida, causing furunculosis, and the virus causing infectious salmon anaemia (ISA), respectively, in a Norwegian aquaculture population of Atlantic salmon. ppar alleles were found to be significantly associated with survival upon challenge to A. salmonicida, but not to ISA. pparg-2 was the better allele in terms of survival in the challenge test for furunculosis, survival rates being 0.32, 0.40 and 0.42 for animals with the pparg-1,-1, pparg-1, -2 and pparg-2, -2 genotypes, respectively. We conclude that pparg-2 is in linkage disequilibrium (LD) with, or identical to, a locus contributing to different susceptibility to furunculosis in Atlantic salmon. PPARG was mapped to linkage group eight (LG8) but could only be positioned on the male linkage map since all the informative parents in the mapping families were males. This is the first report showing an association between pparg alleles and an enhanced immune response in fish.

MHC class IIB additive and non-additive effects on fitness measures in the guppy Poecilia reticulata

The genetic architecture of fitness at the class IIB gene of the major histocompatibility complex (MHC) in the guppy Poecilia reticulata was analysed. Diversity at the MHC is thought to be maintained by some form of balancing selection; heterozygote advantage, frequency-dependent selection or spatially and temporally fluctuating selection. Here these hypotheses are evaluated by using an algorithm that partitions the effect of specific MHC allele and genotypes on fitness measures. The effect of MHC genotype on surrogate measures of fitness was tested, including growth rate (at high and low bulk food diets), parasite load following a parasite challenge and survival. The number of copies of the Pore_a132 MHC allele was inversely related to infection by Gyrodactylus flukes and it appeared to be positively related to faster growth. Also, genotypes combining the Pore_a132 or other relatively common alleles paired with rare MHC alleles produced both advantageous and detrimental non-additive effects. Thus, the genetic architecture underlying fitness at the MHC is complex in the P. reticulata.

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.

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.

MHC standing genetic variation and pathogen resistance in wild Atlantic salmon

Pathogens are increasingly emerging in human-altered environments as a serious threat to biodiversity. In this context of rapid environmental changes, improving our knowledge on the interaction between ecology and evolution is critical. The objective of this study was to evaluate the influence of an immunocompetence gene, the major histocompatibility complex (MHC) class IIbeta, on the pathogen infection levels in wild Atlantic salmon populations, Salmo salar, and identify selective agents involved in contemporary coevolution. MHC variability and bacterial infection rate were determined throughout the summer in juvenile salmon from six rivers belonging to different genetic and ecological regions in Québec, Canada. A total of 13 different pathogens were identified in kidney by DNA sequence analysis, including a predominant myxozoa, most probably recently introduced in North America. Infection rates were the highest in southern rivers at the beginning of the summer (average 47.6+/-6.3% infected fish). One MHC allele conferred a 2.9 times greater chance of being resistant to myxozoa, while another allele increased susceptibility by 3.4 times. The decrease in frequency of the susceptibility allele but not other MHC or microsatellite alleles during summer was suggestive of a mortality event from myxozoa infection. These results supported the hypothesis of pathogen-driven selection in the wild by means of frequency-dependent selection or change in selection through time and space rather than heterozygous advantage, and underline the importance of MHC standing genetic variation for facing pathogens in a changing environment.

Confirmation and fine-mapping of a major QTL for resistance to infectious pancreatic necrosis in Atlantic salmon (Salmo salar): population-level associations between markers and trait

BACKGROUND

Infectious pancreatic necrosis (IPN) is one of the most prevalent and economically devastating diseases in Atlantic salmon (Salmo salar) farming worldwide. The disease causes large mortalities at both the fry- and post-smolt stages. Family selection for increased IPN resistance is performed through the use of controlled challenge tests, where survival rates of sib-groups are recorded. However, since challenge-tested animals cannot be used as breeding candidates, within-family selection is not performed and only half of the genetic variation for IPN resistance is being exploited. DNA markers linked to quantitative trait loci (QTL) affecting IPN resistance would therefore be a powerful selection tool. The aim of this study was to identify and fine-map QTL for IPN-resistance in Atlantic salmon, for use in marker-assisted selection to increase the rate of genetic improvement for this trait.

RESULTS

A genome scan was carried out using 10 large full-sib families of challenge-tested Atlantic salmon post-smolts and microsatellite markers distributed across the genome. One major QTL for IPN-resistance was detected, explaining 29% and 83% of the phenotypic and genetic variances, respectively. This QTL mapped to the same location as a QTL recently detected in a Scottish Atlantic salmon population. The QTL was found to be segregating in 10 out of 20 mapping parents, and subsequent fine-mapping with additional markers narrowed the QTL peak to a 4 cM region on linkage group 21. Challenge-tested fry were used to show that the QTL had the same effect on fry as on post-smolt, with the confidence interval for QTL position in fry overlapping the confidence interval found in post-smolts. A total of 178 parents were tested for segregation of the QTL, identifying 72 QTL-heterozygous parents. Genotypes at QTL-heterozygous parents were used to determine linkage phases between alleles at the underlying DNA polymorphism and alleles at single markers or multi-marker haplotypes. One four-marker haplotype was found to be the best predictor of QTL alleles, and was successfully used to deduce genotypes of the underlying polymorphism in 72% of the parents of the next generation within a breeding nucleus. A highly significant population-level correlation was found between deduced alleles at the underlying polymorphism and survival of offspring groups in the fry challenge test, parents with the three deduced genotypes (QQ, Qq, qq) having mean offspring mortality rates of 0.13, 0.32, and 0.49, respectively. The frequency of the high-resistance allele (Q) in the population was estimated to be 0.30. Apart from this major QTL, one other experiment-wise significant QTL for IPN-resistance was detected, located on linkage group 4.

CONCLUSION

The QTL confirmed in this study represents a case of a major gene explaining the bulk of genetic variation for a presumed complex trait. QTL genotypes were deduced within most parents of the 2005 generation of a major breeding company, providing a solid framework for linkage-based MAS within the whole population in subsequent generations. Since haplotype-trait associations valid at the population level were found, there is also a potential for MAS based on linkage disequilibrium (LD). However, in order to use MAS across many generations without reassessment of linkage phases between markers and the underlying polymorphism, the QTL needs to be positioned with even greater accuracy. This will require higher marker densities than are currently available.

The polymorphism of lysozyme gene in Zhikong scallop (Chlamys farreri) and its association with susceptibility/resistance to Listonella anguillarum

Lysozyme functions as a crucial biodefence effector against the infection of bacterial pathogens in innate immunity. The nucleotide sequence polymorphisms in promoter region of a nuclear goose type lysozyme gene from Zhikong scallop Chlamys farreri (designated as CFLysG) were investigated to explore their association with susceptibility/resistance to Listonella anguillarum infection. Eight sites of single nucleotide polymorphisms (SNPs) and two sites of insert-deletion (ins-del) polymorphisms were identified in the promoter region of CFLysG. Two of them, -753 TATCTCGATCAGG ins-del polymorphism and -391 A-G SNP were selected to analyze their distribution in the susceptible and resistant stocks, which were identified according to the survival time after L. anguillarum challenge. Using polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP), two genotypes were found at each site, which were ins/del and ins/ins at locus -753, and A/A and A/G at locus -391, respectively. The -753 ins/del genotype was more prevalent in the resistant stock than that in the susceptible stock, 30% vs 16.67% in frequency, but there was no significant difference in the frequency distribution between these two stocks (P=0.15). In contrast, the frequency of -391A/G genotype in the resistant stock was significantly higher (30%) than that in the susceptible stock (7.14%) (P=0.007), indicating a significant association with the resistance of Zhikong scallop to L. anguillarum. To confirm the presumption, another independent challenge experiment was performed, in which the cumulative mortality of scallops with -391 A/A genotype (96.8%) was significantly higher than those with -391 A/G genotype (64.5%) (P=0.001), which further validate the association between -391 A/G genotype and the resistance of Zhikong scallop to L. anguillarum. These results suggested that the -391 A/G could be a potential marker applied in future selection of Zhikong scallop with enhanced resistance to L. anguillarum.

Molecular cloning, genomic structure, polymorphism and expression analysis of major histocompatibility complex class IIA and IIB genes of half-smooth tongue sole (Cynoglossus semilaevis)

Major histocompatibility complex (MHC) genes play an important role in the immune response of vertebrates. Its function is to present foreign peptide to the T-cell. In order to study the function and molecular polymorphism of class II genes in teleost, the full lengths of MHC class IIA and IIB cDNA were cloned from half-smooth tongue sole by homology cloning and rapid amplification of cDNA ends polymerase chain reaction (RACE-PCR). Genomic organizations, molecular polymorphism, and expression profiles of class IIA and IIB were examined to study the function in fish. As in other teleost, four exons and three introns were identified in half-smooth tongue sole class IIA gene, five exons and four introns were identified in class IIB gene. The deduced amino acid sequence of class IIA had 27.3-69.8% identity with those of mammal and teleost. Nine class IIA alleles were identified from four individuals. Four different alleles observed in a single individual may infer the existence of two loci at least. The deduced amino acid sequence of class IIB had 7.9-71.9% identity with those of other species. Fifteen class IIB alleles were identified. Six different alleles observed in a single individual may suggest that there are at least three loci in class IIB genes. Real-time quantitative RT-PCR demonstrated that the MHC class IIA and IIB were ubiquitously expressed in twelve normal tissues. Challenge of half-smooth tongue sole with the pathogenic bacteria, Vibrio anguillarum, resulted in significant changes in the expression of MHC IIA and IIB mRNA in three tissues.

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.

Lifetime reproductive success is maximized with optimal major histocompatibility complex diversity

Individual diversity at the major histocompatibility complex (MHC) is predicted to be optimal at intermediate rather than at maximal levels. We showed previously in sticklebacks that an intermediate MHC diversity is predominant in natural populations and provides maximal resistance in experimental multiple parasite infections in the laboratory. However, what counts ultimately is the lifetime reproductive success (LRS). Here, we measured LRS of six laboratory-bred sib-groups-to minimize the influence of non-MHC genes-three-spined sticklebacks (Gasterosteus aculeatus) during their entire breeding period, each in a seminatural enclosure in the lake of their parents, where they were exposed to the natural spectrum of parasites. We collected developing clutches at regular intervals and determined parenthood for a representative number of eggs (2279 in total) per clutch with 18 microsatellites. Both males and females with an intermediate MHC class IIB variant number had the highest LRS. The mechanistic link of MHC diversity and LRS differed between the sexes: in females, we found evidence for a trade-off between number of eggs and immunocompentence, whereas in males this correlation was concealed by different timing strategies of reproduction.

A review of the need and possible uses for genetically standardized Atlantic salmon (Salmo salar) in research

Large numbers of Atlantic salmon ( Salmo salar) are used as research animals in basic research and to solve challenges related to the fish-farming industry. Most of this research is performed on farmed animals provided by local breeders or national breeding companies. The genetic constitution of these animals is usually unknown and highly variable. As a result, large numbers of fish are often needed to produce significant results, and results from one study are often impossible to reproduce in another facility. The production of standardized salmon could in many cases reduce the number of animals used in research and at the same time provide more reproducible results. This paper provides an overview of the methods available for the production of standardized Atlantic salmon, and discusses the pros and cons of each technique. The use of zebrafish and other well-defined laboratory fish species as a model for salmon is also discussed. Access to genetically defined fish would greatly benefit the scientific community, in the same way as genetically defined lines of rodents have revolutionized mammalian research.

Major histocompatibility complex (MHC) heterozygote superiority to natural multi-parasite infections in the water vole (Arvicola terrestris)

The fundamental role of the major histocompatibility complex (MHC) in immune recognition has led to a general consensus that the characteristically high levels of functional polymorphism at MHC genes is maintained by balancing selection operating through host-parasite coevolution. However, the actual mechanism by which selection operates is unclear. Two hypotheses have been proposed: overdominance (or heterozygote superiority) and negative frequency-dependent selection. Evidence for these hypotheses was evaluated by examining MHC-parasite relationships in an island population of water voles (Arvicola terrestris). Generalized linear mixed models were used to examine whether individual variation at an MHC class II DRB locus explained variation in the individual burdens of five different parasites. MHC genotype explained a significant amount of variation in the burden of gamasid mites, fleas (Megabothris walkeri) and nymphs of sheep ticks (Ixodes ricinus). Additionally, MHC heterozygotes were simultaneously co-infected by fewer parasite types than homozygotes. In each case where an MHC-dependent effect on parasite burden was resolved, the heterozygote genotype was associated with fewer parasites, and the heterozygote outperformed each homozygote in two of three cases, suggesting an overall superiority against parasitism for MHC heterozygote genotypes. This is the first demonstration of MHC heterozygote superiority against multiple parasites in a natural population, a mechanism that could help maintain high levels of functional MHC genetic diversity in natural populations.

Parental relatedness and major histocompatibility effects on early embryo survivorship in Atlantic salmon

Salmon have provided key insights into the relative influence of natural and sexual selection on major histocompatibility complex (MHC) variation. Natural selection on salmon MHC genes has been demonstrated in pathogen studies, and there is evidence of MHC-based mate choice (sexual selection). We tested whether parental MHC genes affect survivorship of juvenile Atlantic salmon (Salmo salar) by quantifying the influence of parental genome-wide relatedness and MHC genotype on survivorship to the swim-up stage. Thirteen microsatellite loci were used to estimate the influence of genome-wide relatedness between parents on offspring survivorship and MHC genotypes were determined by sequencing part of the class IIbeta gene. Our results revealed no significant relationship between early offspring survivorship and genome-wide relatedness, predicted MHC heterozygosity, or MHC allelic similarity. Overall, our data are consistent with the contention that excess MHC heterozygosity in Atlantic salmon juveniles is due to sexual selection as well as differential survival of offspring due to MHC genotype.

Does intra-individual major histocompatibility complex diversity keep a golden mean?

An adaptive immune response is usually initiated only if a major histocompatibility complex (MHC) molecule presents pathogen-derived peptides to T-cells. Every MHC molecule can present only peptides that match its peptide-binding groove. Thus, it seems advantageous for an individual to express many different MHC molecules to be able to resist many different pathogens. However, although MHC genes are the most polymorphic genes of vertebrates, each individual has only a very small subset of the diversity at the population level. This is an evolutionary paradox. We provide an overview of the current data on infection studies and mate-choice experiments and conclude that overall evidence suggests that intermediate intra-individual MHC diversity is optimal. Selective forces that may set an upper limit to intra-individual MHC diversity are discussed. An updated mathematical model based on recent findings on T-cell selection can predict the natural range of intra-individual MHC diversity. Thus, the aim of our review is to evaluate whether the number of MHC alleles usually present in individuals may be optimal to balance the advantages of presenting an increased range of peptides versus the disadvantages of an increased loss of T-cells.

Cell markers and determinants in fish immunology

Despite the impressive increase in the cloning and expression of genes encoding fish immunoregulatory molecules, the knowledge on "in vivo" and "in vitro" functional immunology of the corresponding peptide products is still at an initial stage. This is partly due to the lacking of specific markers for immunoregulatory peptides, that represent an indispensible tool to dissect immune reactions and to trace the fate of cellular events downstream of the activation. In this review we summarise the available information on functional immune activities of some teleost species and discuss the obtained data in an evolutionary and applied context.

Peptide-binding motif prediction by using phage display library for SasaUBA*0301, a resistance haplotype of MHC class I molecule from Atlantic Salmon (Salmo salar)

The structure of the peptide-binding specificity of major histocompatibility complex (MHC) class I has been analyzed extensively in human and mouse. For fish, there are no crystallographic models of MHC molecules, neither are there data on the peptide-binding specificity. In this study, we describe for the first time the identification of a fish class I peptide-MHC ligand binding motif. Phage display technology using both 7 mer and 12 mer libraries enabled us to identify peptide ligands with unique specificity that interacts with the recombinant Salmon MHC class I molecule. The recombinant proteins, beta 2m/SasaUBA*0301, were produced in Escherichia coli, in which the carboxyl terminus of beta 2-microglobulin is joined together with a flexible (GGGGS)3 linker to the amino terminus of the heavy chain. One hundred and seven individual phages bound to beta 2m/SasaUBA*0301 were isolated after four rounds of panning from the 7 mer random-peptide library. The peptide encoding sequences were determined and peptide alignment led to the prediction of position-specific anchor residue. A prominent proline at position 2 was observed and we predict that it might be one of the anchors at the N-terminus. Meanwhile, phage display peptide library encoding random 12 mer peptides was also screened against beta 2m/SasaUBA*0301. Eighty-five percentages of the corresponding peptides have an enrichment of leucine, methionine, valine, or isoleucine at the C-terminus. We predict that this particular allele of Salmon class I molecule might have a very similar binding motif at the C-terminus compared with a known mouse class I molecule H2-Kb which has L, or I, V, M at p8. Previous work showed that Atlantic Salmon carrying the allele SasaUBA*0301 are resistant to infectious Salmon aneamia virus and there is a significant association between MHC polymorphism and the disease resistance. Therefore, our study might contribute to designing a peptide vaccine against this viral disease.

MHC-mediated mate choice increases parasite resistance in salmon

Natural (parasite-driven) and sexual selection are thought to maintain high polymorphism in the genes of the major histocompatibility complex (MHC), but support for a link between mate choice, MHC variation and increased parasite resistance is circumstantial. We compared MHC diversity and Anisakis loads among anadromous Atlantic salmon (Salmo salar L.) returning to four rivers to spawn, which had originated from natural spawning (parents allowed to mate freely) or artificial crosses (parents deprived from the potential benefits of mate choice). We found that the offspring of artificially bred salmon had higher parasite loads and were almost four times more likely to be infected than free-mating salmon, despite having similar levels of MHC diversity. Moreover, the offspring of wild salmon were more MHC dissimilar than the offspring of artificially crossed salmon, and uninfected fish were more dissimilar for MHC than infected fish. Thus, our results suggest a link between disassortative mating and offspring benefits and indicate that MHC-mediated mate choice and natural (parasite-driven) selection act in combination to maintain MHC diversity, and hence fitness. Therefore, artificial breeding programmes that negate the potential genetic benefits of mate choice may result in inherently inferior offspring, regardless of population size, rearing conditions or genetic diversity.

Major histocompatibility complex based resistance to a common bacterial pathogen of amphibians

Given their well-developed systems of innate and adaptive immunity, global population declines of amphibians are particularly perplexing. To investigate the role of the major histocompatibility complex (MHC) in conferring pathogen resistance, we challenged Xenopus laevis tadpoles bearing different combinations of four MHC haplotypes (f, g, j, and r) with the bacterial pathogen Aeromonas hydrophila in two experiments. In the first, we exposed ff, fg, gg, gj, and jj tadpoles, obtained from breeding MHC homozygous parents, to one of three doses of A. hydrophila or heat-killed bacteria as a control. In the second, we exposed ff, fg, fr, gg, rg, and rr tadpoles, obtained from breeding MHC heterozygous parents and subsequently genotyped by PCR, to A. hydrophila, heat-killed bacteria or media alone as controls. We thereby determined whether the same patterns of MHC resistance emerged within as among families, independent of non-MHC heritable differences. Tadpoles with r or g MHC haplotypes were more likely to die than were those with f or j haplotypes. Growth rates varied among MHC types, independent of exposure dose. Heterozygous individuals with both susceptible and resistant haplotypes were intermediate to either homozygous genotype in both size and survival. The effect of the MHC on growth and survival was consistent between experiments and across families. MHC alleles differentially confer resistance to, or tolerance of, the bacterial pathogen, which affects tadpoles' growth and survival.

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

Genes within the major histocompatibility complex (MHC) are important for both innate and adaptive immune responses in mammals; however, much less is known regarding their contribution in teleost fishes. We examined the involvement of four major histocompatibility (MH) genomic regions in rainbow trout in resistance to the causative agent of bacterial coldwater disease (BCWD), Flavobacterium psychrophilum. Fish from the 2005 NCCCWA brood-year (71 full-sib families) were challenged with F. psychrophilum strain CSF 259-93. The overall mortality rate was 70%, with large variation in mortality between families. Disease resistance was quantified as post-challenge days to death. Phenotypic variation and additive genetic variation were estimated using mixed models of survival analysis. To examine association, eight microsatellite markers were isolated from MH gene-containing BAC clones and mapped onto the rainbow trout genetic linkage map. The parents and grandparents of the 2005 brood-year families were genotyped with these eight markers and another two markers tightly linked to the MH-IB region to assess the extent of linkage disequilibrium (LD) of MH genomic regions MH-IA, MH-IB, TAP1, and MH-II with survival post-challenge. MH-IB and MH-II markers were linked to BCWD survivability when data were analyzed by family. Tests for disease association at the population level substantiated the involvement of MH-IB, but not MH-II, with disease resistance. The impact of selective breeding for disease resistance on MH sequence variation is discussed in the context of aquaculture production.

Identification of immune-related genes from kidney and spleen of turbot, Psetta maxima (L.), by suppression subtractive hybridization following challenge with Vibrio harveyi

Suppression subtractive hybridization was used to investigate the response of turbot, Psetta maxima (L.), to Vibrio harveyi, by using a cDNA library constructed from artificially infected turbot kidney and spleen mRNA. Forty-nine expressed sequence tags were obtained. Several immune system genes were identified, including a major histocompatibility complex (MHC) class Ia gene and a heat shock protein 70 gene. Some signalling molecules were also present in the cDNA libraries, including src-family tyrosine kinase SCK, sgk-1 serine-threonine protein kinase and amyloid precursor-like protein 2. The full length of MHC class Ia cDNA was cloned from turbot cDNA by rapid amplification of cDNA ends polymerase chain reaction. The nucleotide sequence of turbot MHC class Ia has been submitted to GenBank with accession number EF032639. The turbot MHC class Ia cDNA has an open reading frame encoding 354 amino acids, and the deduced amino acid sequence of turbot MHC class Ia has 68%, 54%, 51%, 52%, 57%, 33%, 29% and 29% identities to those of olive flounder, medaka, rainbow trout, Atlantic cod, tiger puffer, chicken, mouse and human, respectively. Quantitative reverse transcriptase-PCR was performed for the MHC class Ia gene, and it was revealed that the expression level of the MHC class Ia gene in V. harveyi-challenged turbot increased to fourfold that of the controls.

The MHC and non-random mating in a captive population of Chinook salmon

Detailed analysis of variation in reproductive success can provide an understanding of the selective pressures that drive the evolution of adaptations. Here, we use experimental spawning channels to assess phenotypic and genotypic correlates of reproductive success in Chinook salmon (Oncorhynchus tshawytscha). Groups of 36 fish in three different sex ratios (1:2, 1:1 and 2:1) were allowed to spawn and the offspring were collected after emergence from the gravel. Microsatellite genetic markers were used to assign parentage of each offspring, and the parents were also typed at the major histocompatibility class IIB locus (MHC). We found that large males, and males with brighter coloration and a more green/blue hue on their lateral integument sired more offspring, albeit only body size and brightness had independent effects. There was no similar relationship between these variables and female reproductive success. Furthermore, there was no effect of sex ratio on the strength or significance of any of the correlations. Females mated non-randomly at the MHC, appearing to select mates that produced offspring with greater genetic diversity as measured by amino-acid divergence. Females mated randomly with respect to male genetic relatedness and males mated randomly with respect to both MHC and genetic relatedness. These results indicate that sexual selection favours increased body size and perhaps integument coloration in males as well as increases genetic diversity at the MHC by female mate choice.

MH class IIalpha polymorphism in local and global adaptation of Arctic charr (Salvelinus alpinus L.)

Arctic charr, a highly plastic salmonid that inhabits the circumpolar region, colonized its current environment after the last glaciation. Recent colonization limits the capacity of many techniques to define and characterize constituent populations. As a novel approach, we used the major histocompatibility (MH) class IIalpha gene polymorphism as a marker that would characterize the genetic divergence of global Arctic charr populations caused by drift and by local adaptation to pathogens. We were able to detect significant isolation of all the lineages previously defined by mitochondrial DNA sequencing and also isolation of some populations within those groups. We found that most of the polymorphism of the class IIalpha gene was distributed globally, which indicates ancestral selection; however, in most cases, distinctive allele frequencies and specific haplotypes distinguished each population suggesting that recent selection has also occurred. Although all studied populations showed similar MH class IIalpha polymorphisms, we also found variation in which particular amino acid positions were polymorphic and which were constant in the different populations studied. This variation provides a greater adaptive capacity for the MH class IIalpha receptors in Arctic charr and is yet another illustration of the extraordinary plasticity of the species.

Mortality selection during the 2003 European heat wave in three-spined sticklebacks: effects of parasites and MHC genotype

Background

Ecological interaction strength may increase under environmental stress including temperature. How such stress enhances and interacts with parasite selection is almost unknown. We studied the importance of resistance genes of the major histocompatibility complex (MHC) class II in 14 families of three-spined sticklebacks Gasterosteus aculeatus exposed to their natural macroparasites in field enclosures in the extreme summer of 2003.

Results

After a mass die-off during the 2003-European heat wave killing 78% of 277 experimental fish, we found strong differences in survival among and within families. In families with higher average parasite load fewer individuals survived. Multivariate analysis revealed that the composition of the infecting parasite fauna was family specific. Within families, individuals with an intermediate number of MHC class IIB sequence variants survived best and had the lowest parasite load among survivors, suggesting a direct functional link between MHC diversity and fitness. The within family MHC effects were, however, small compared to between family effects, suggesting that other genetic components or non-genetic effects were also important.

Conclusion

The correlation between parasite load and mortality that we found at both individual and family level might have appeared only in the extraordinary heatwave of 2003. Due to global warming the frequency of extreme climatic events is predicted to increase, which might intensify costs of parasitism and enhance selection on immune genes.

Gene expression analyses in Atlantic salmon challenged with infectious salmon anemia virus reveal differences between individuals with early, intermediate and late mortality

BACKGROUND

Infectious salmon anemia virus (ISAV) causes a multisystemic disease responsible for severe losses in salmon aquaculture. Better understanding of factors that explain variations in resistance between individuals and families is essential for development of strategies for disease control. To approach this, we compared global gene expression using microarrays in fish dying early and late in the time course following infection from a highly pathogenic ISAV.

RESULTS

Tissues (gill, heart, liver and spleen) from infected Atlantic salmon (cohabitation, ISAV Glesvaer 2/90 isolate) were collected from three stages over the time course of the experiment; early (EM, 0-10% cumulative mortality (CM), 21-25 days post-infection (DPI)), intermediate (IM, 35-55% CM, 28-31 DPI) and late (LM, 75-85% CM, 37-48 DPI) mortality. Viral loads were equal in EM and IM but dropped markedly in LM fish. Gene expression analyses using a 1.8 K salmonid fish cDNA microarray (SFA2.0) and real-time qPCR revealed a large group of genes highly up-regulated across tissues in EM, which were mainly implicated in innate antiviral responses and cellular stress. Despite equal levels of MHC class I in EM and LM, increase of splenic and cardiac expression of immunoglobulin-like genes was found only in LM while a suite of adaptive immunity markers were activated already in IM. The hepatic responses to ISAV were characterized by difference between EM and LM in expression of chaperones and genes involved in eicosanoid metabolism. To develop classification of high and low resistance phenotypes based on a small number of genes, we processed results from qPCR analyses of liver using a linear discriminant analysis. Four genes (5-lipoxygenase activating protein, cytochrome P450 2K4-1, galectin-9 and annexin A1) were sufficient for correct assignment of individuals to EM, LM and uninfected groups, while IM was inseparable from EM. Three of four prognostic markers are involved in metabolism of inflammatory regulators.

CONCLUSION

This study adds to the understanding of molecular determinants for resistance to acute ISAV infection. The most susceptible individuals were characterized by high viral replication and dramatic activation of innate immune responses, which did not provide protection. The ability to endure high levels of infection for sustained periods could be associated with lower inflammatory responses while subsequent protection and viral clearance was most likely conferred by activation of adaptive immunity.

Expression and characterization of recombinant single-chain salmon class I MHC fused with beta2-microglobulin with biological activity

Heterodimeric class I major histocompatibility complex (MHC) molecules consist of a putative 45-kDa heavy chain and a 12-kDa beta2-microglobulin (beta2m) light chain. The knowledge about MHC genes in Atlantic salmon accumulated during the last decade has allowed us to generate soluble and stable MHC class I molecules with biological activity. We report here the use of a bacterial expression system to produce the recombinant single-chain MHC molecules based on a specific allele Sasa-UBA*0301. This particular allele was selected because previous work has shown its association with the resistance to infectious salmon anaemia virus. The single-chain salmon MHC class I molecule has been designed and generated, in which the carboxyl terminus of beta2m is joined together with a flexible 15 or 20 amino acid peptide linker to the amino terminus of the heavy chain (Sasabeta2mUBA*0301). Monoclonal antibodies were successfully produced against both the MHC class I heavy chain and beta(2)m, and showed binding to the recombinant molecule. The recombinant complex Sasabeta2mUBA*0301 was expressed and isolated; the production was scaled up by adjusting to its optimal conditions. Subsequently, the recombinant proteins were purified by affinity chromatography using mAb against beta2m and alpha3. Eluates were analyzed by Western blot and refolded by the removal of denaturant. The correct folding was confirmed by measuring its binding capacity against mAb produced to recognize the native form of MHC molecules by biosensor analysis. This production of sufficient amounts of class I MHC proteins may represent a useful tool to study the peptide-binding specificity of MHC class I molecules, in order to design a peptide vaccine against viral pathogens.