A single major chromosomal region controls natural killer cell-like activity in rainbow trout

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.

Production and verification of a 2nd generation clonal group of Japanese flounder, Paralichthys olivaceus

Clonal fishes are useful tools in biology and aquaculture studies due to their isogenicity. In Japanese flounder (Paralichthys olivaceus), a group of homozygous clones was created by inducing meiogynogenesis in eggs from a mitogynogenetic homozygous diploid. As the clones reached sexual maturity, meiogynogenesis was again induced in order to produce a 2^nd generation clonal group of Japanese flounder. After 3 months, there were 611 healthy, surviving individuals. Twenty-four microsatellite markers, that covered all the linkage groups of Japanese flounder, were used to identify the homozygosity of the 2^nd generation clones; no heterozygous locus was detected. This indicates that the production of a 2^nd generation clonal group of Japanese flounder was successful. Restriction-site DNA associated sequencing at the genomic level also confirmed the homozygosity and clonality of the 2^nd generation clonal group. Furthermore, these 2^nd generation clones had a small coefficient of variation for body shape indices at 210 days of age and showed a high degree of similarity in body characteristics among individuals. The successful production of 2^nd generation clones has laid the foundation for the large-scale production of clonal Japanese flounder.

Putative apolipoprotein A-I, natural killer cell enhancement factor and lysozyme g are involved in the early immune response of brown-marbled grouper, Epinephelus fuscoguttatus, Forskal, to Vibrio alginolyticus

The gram-negative bacterium, Vibrio alginolyticus, has frequently been identified as the pathogen responsible for the infectious disease called vibriosis. This disease is one of the major challenges facing brown-marbled grouper aquaculture, causing fish farmers globally to suffer substantial economic losses. The objective of this study was to investigate the proteins involved in the immune response of brown-marbled grouper fingerlings during their initial encounter with pathogenic organisms. To achieve this objective, a challenge experiment was performed, in which healthy brown-marbled grouper fingerlings were divided into two groups. Fish in the treated group were subjected to intraperitoneal injection with an infectious dose of V. alginolyticus suspended in phosphate-buffered saline (PBS), and those in the control group were injected with an equal volume of PBS. Blood samples were collected from a replicate number of fish from both groups at 4 h post-challenge and analysed for immune response-related serum proteins via two-dimensional gel electrophoresis. The results showed that 14 protein spots were altered between the treated and control groups; these protein spots were further analysed to determine the identity of each protein via MALDI-TOF/TOF. Among the altered proteins, three were clearly overexpressed in the treated group compared with the control; these were identified as putative apolipoprotein A-I, natural killer cell enhancement factor and lysozyme g. Based on these results, these three highly expressed proteins participate in immune response-related reactions during the initial exposure (4 h) of brown-marbled grouper fingerling to V. alginolyticus infection.

Transfer of serum and cells from Yersinia ruckeri vaccinated doubled-haploid hot creek rainbow trout into outcross F1 progeny elucidates mechanisms of vaccine-induced protection

Yersinia ruckeri is a well-established bacterial pathogen for many salmonid species, against which a formalin-killed bacterin vaccine has been effective in reducing disease outbreaks. Previous studies have reported conflicting results about the protective value of the systemic humoral response to Y. ruckeri vaccination. Here we directly demonstrate that plasma contains the long-term protective component elicited by both immersion and intraperitoneal injection vaccination of rainbow trout. A total of 0.5 μL of plasma from vaccinated fish provided almost complete protection against experimental challenge. Conversely, the cells obtained from peripheral blood conferred little or no protection in naïve recipients. The protective component of immune sera was IgM based on size exclusion chromatography and recognition by monoclonal antibody Warr 1-14. Immune plasma generated against a Y. ruckeri biotype 1 strain protected equally against challenges with Y. ruckeri biotype 1 and 2 strains. These results illustrate the importance of the humoral IgM response against Y. ruckeri and the use of doubled haploid rainbow trout (Oncorhynchus mykiss) and transfer of plasma/serum and cells into F1 outcross progeny as a model system for dissection of the mechanism(s) of vaccine-induced protection.

Oxidative stress among SOD-1 genotypes in rainbow trout (Oncorhynchus mykiss)

Natural variation in the antioxidant-enzyme SOD-1 (superoxide dismutase) is known to alter the impacts of oxidative damage at both the cellular and organismal levels. Using three homozygous clonal lines of rainbow trout [Hot Creek (n=30), Arlee (n=21), and Swanson (n=10)], which differ for single nucleotide polymorphisms (SNPs) and amino acid substitutions at the SOD-1 locus, we investigated the functional effects of this variation on SOD-1 activity during ozone stress and subsequent levels of oxidative damage to DNA and cell membranes. Fish from each line were subjected to either control conditions or 24h of ozone stress, after which tissues were analyzed for antioxidant status and oxidative damage. Liver SOD-1 activity was lower in ozonated vs. control fish in the Hot Creek line, and among ozonated fish, Hot Creek was lower than Arlee. Total erythrocyte SOD activity was not significantly impacted by ozonation; however significant differences in total erythrocyte SOD activity were measured among clonal lines, driven primarily by lower activity in the Hot Creek line. Ozone had a significant treatment effect in all oxidative damage parameters assessed: it increased DNA lesions in erythrocytes and levels of lipid peroxidation in gill tissue and plasma. Among lines, Swanson showed higher lipid peroxidation levels in gill tissue after ozonation than Arlee or Hot Creek. Conversely, Swanson control and treatment fish had significantly lower plasma lipid peroxidation levels than did fish from the other lines. Overall, the among-line differences in SOD and SOD-1 activity and oxidative damage provide evidence that SOD-1 genotypes differ functionally under both oxidative stress and control conditions; however, other genetic differences among lines should be investigated in order to further explain the phenotypic differences in SOD enzyme activity and oxidative damage described here.

Characterization of the OmyY1 Region on the Rainbow Trout Y Chromosome

We characterized the male-specific region on the Y chromosome of rainbow trout, which contains both sdY (the sex-determining gene) and the male-specific genetic marker, OmyY1. Several clones containing the OmyY1 marker were screened from a BAC library from a YY clonal line and found to be part of an 800 kb BAC contig. Using fluorescence in situ hybridization (FISH), these clones were localized to the end of the short arm of the Y chromosome in rainbow trout, with an additional signal on the end of the X chromosome in many cells. We sequenced a minimum tiling path of these clones using Illumina and 454 pyrosequencing. The region is rich in transposons and rDNA, but also appears to contain several single-copy protein-coding genes. Most of these genes are also found on the X chromosome; and in several cases sex-specific SNPs in these genes were identified between the male (YY) and female (XX) homozygous clonal lines. Additional genes were identified by hybridization of the BACs to the cGRASP salmonid 4x44K oligo microarray. By BLASTn evaluations using hypothetical transcripts of OmyY1-linked candidate genes as query against several EST databases, we conclude at least 12 of these candidate genes are likely functional, and expressed.

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.

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

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

Analysis of BAC-end sequences in rainbow trout: content characterization and assessment of synteny between trout and other fish genomes

BACKGROUND

Rainbow trout (Oncorhynchus mykiss) are cultivated worldwide for aquaculture production and are widely used as a model species to gain knowledge of many aspects of fish biology. The common ancestor of the salmonids experienced a whole genome duplication event, making extant salmonids such as the rainbow trout an excellent model for studying the evolution of tetraploidization and re-diploidization in vertebrates. However, the lack of a reference genome sequence hampers research progress for both academic and applied purposes. In order to enrich the genomic tools already available in this species and provide further insight on the complexity of its genome, we sequenced a large number of rainbow trout BAC-end sequences (BES) and characterized their contents.

RESULTS

A total of 176,485 high quality BES, were generated, representing approximately 4% of the trout genome. BES analyses identified 6,848 simple sequence repeats (SSRs), of which 3,854 had high quality flanking sequences for PCR primers design. The first rainbow trout repeat elements database (INRA RT rep1.0) containing 735 putative repeat elements was developed, and identified almost 59.5% of the BES database in base-pairs as repetitive sequence. Approximately 55% of the BES reads (97,846) had more than 100 base pairs of contiguous non-repetitive sequences. The fractions of the 97,846 non-repetitive trout BES reads that had significant BLASTN hits against the zebrafish, medaka and stickleback genome databases were 15%, 16.2% and 17.9%, respectively, while the fractions of the non-repetitive BES reads that had significant BLASTX hits against the zebrafish, medaka, and stickleback protein databases were 10.7%, 9.5% and 9.5%, respectively. Comparative genomics using paired BAC-ends revealed several regions of conserved synteny across all the fish species analyzed in this study.

CONCLUSIONS

The characterization of BES provided insights on the rainbow trout genome. The discovery of specific repeat elements will facilitate analyses of sequence content (e.g. for SNPs discovery and for transcriptome characterization) and future genome sequence assemblies. The numerous microsatellites will facilitate integration of the linkage and physical maps and serve as valuable resource for fine mapping QTL and positional cloning of genes affecting aquaculture production traits. Furthermore, comparative genomics through BES can be used for identifying positional candidate genes from QTL mapping studies, aid in future assembly of a reference genome sequence and elucidating sequence content and complexity in the rainbow trout genome.

Detection of QTL with effects on osmoregulation capacities in the rainbow trout (Oncorhynchus mykiss)

Background

There is increasing evidence that the ability to adapt to seawater in teleost fish is modulated by genetic factors. Most studies have involved the comparison of species or strains and little is known about the genetic architecture of the trait. To address this question, we searched for QTL affecting osmoregulation capacities after transfer to saline water in a nonmigratory captive-bred population of rainbow trout.

Results

A QTL design (5 full-sib families, about 200 F2 progeny each) was produced from a cross between F0 grand-parents previously selected during two generations for a high or a low cortisol response after a standardized confinement stress. When fish were about 18 months old (near 204 g body weight), individual progeny were submitted to two successive hyper-osmotic challenges (30 ppt salinity) 14 days apart. Plasma chloride and sodium concentrations were recorded 24 h after each transfer. After the second challenge, fish were sacrificed and a gill index (weight of total gill arches corrected for body weight) was recorded. The genome scan was performed with 196 microsatellites and 85 SNP markers. Unitrait and multiple-trait QTL analyses were carried out on the whole dataset (5 families) through interval mapping methods with the QTLMap software. For post-challenge plasma ion concentrations, significant QTL (P < 0.05) were found on six different linkage groups and highly suggestive ones (P < 0.10) on two additional linkage groups. Most QTL affected concentrations of both chloride and sodium during both challenges, but some were specific to either chloride (2 QTL) or sodium (1 QTL) concentrations. Six QTL (4 significant, 2 suggestive) affecting gill index were discovered. Two were specific to the trait, while the others were also identified as QTL for post-challenge ion concentrations. Altogether, allelic effects were consistent for QTL affecting chloride and sodium concentrations but inconsistent for QTL affecting ion concentrations and gill morphology. There was no systematic lineage effect (grand-parental origin of QTL alleles) on the recorded traits.

Conclusions

For the first time, genomic loci associated with effects on major physiological components of osmotic adaptation to seawater in a nonmigratory fish were revealed. The results pave the way for further deciphering of the complex regulatory mechanisms underlying seawater adaptation and genes involved in osmoregulatory physiology in rainbow trout and other euryhaline fishes.

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.

Production of viable homozygous, doubled haploid channel catfish (Ictalurus punctatus)

Production of doubled haploids via mitotic gynogenesis is a useful tool for the creation of completely inbred fish. In order to produce viable doubled haploid channel catfish, we utilized hydrostatic pressure or thermal treatments on eggs fertilized with sperm that had been exposed to ultraviolet light. At 1.5 h post-fertilization, the embryos were exposed to either 590 kg/cm(2) hydrostatic pressure for 3 min, 37 degrees C for 5 min, or 41 degrees C for 3 min. In the pressure-treated group, only 21 offspring hatched from five spawns with family sizes of one, two, two, four, and 12 offspring each. Eight embryos from the 37 degrees C treatment and 32 embryos from the 41 degrees C treatment survived to hatch. Genotype analysis using microsatellite loci demonstrated all 21 offspring resulting from pressure treatment were homozygous at the 64 loci tested, and none contained alleles unique to the donor male. Eleven of 32 offspring from the 41 degrees C treatment were homozygous at the 18 loci tested, while 21 offspring were heterozygous at six to 12 of these loci. Again, no offspring contained alleles unique to the donor male. However, all eight offspring from the 37 degrees C treatment were heterozygous at multiple loci, and one contained unambiguous paternal alleles. These experiments demonstrated our ability to produce viable homozygous, doubled haploid channel catfish. Doubled haploid catfish can be used to create completely inbred populations for genetic analyses, and homozygous genomic templates will be useful in gene identification and genome characterization.

Estimates of linkage disequilibrium and effective population size in rainbow trout

Background

The use of molecular genetic technologies for broodstock management and selective breeding of aquaculture species is becoming increasingly more common with the continued development of genome tools and reagents. Several laboratories have produced genetic maps for rainbow trout to aid in the identification of loci affecting phenotypes of interest. These maps have resulted in the identification of many quantitative/qualitative trait loci affecting phenotypic variation in traits associated with albinism, disease resistance, temperature tolerance, sex determination, embryonic development rate, spawning date, condition factor and growth. Unfortunately, the elucidation of the precise allelic variation and/or genes underlying phenotypic diversity has yet to be achieved in this species having low marker densities and lacking a whole genome reference sequence. Experimental designs which integrate segregation analyses with linkage disequilibrium (LD) approaches facilitate the discovery of genes affecting important traits. To date the extent of LD has been characterized for humans and several agriculturally important livestock species but not for rainbow trout.

Results

We observed that the level of LD between syntenic loci decayed rapidly at distances greater than 2 cM which is similar to observations of LD in other agriculturally important species including cattle, sheep, pigs and chickens. However, in some cases significant LD was also observed up to 50 cM. Our estimate of effective population size based on genome wide estimates of LD for the NCCCWA broodstock population was 145, indicating that this population will respond well to high selection intensity. However, the range of effective population size based on individual chromosomes was 75.51 - 203.35, possibly indicating that suites of genes on each chromosome are disproportionately under selection pressures.

Conclusions

Our results indicate that large numbers of markers, more than are currently available for this species, will be required to enable the use of genome-wide integrated mapping approaches aimed at identifying genes of interest in rainbow trout.

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.

Recombination is suppressed over a large region of the rainbow trout Y chromosome

The previous genetic mapping data have suggested that most of the rainbow trout sex chromosome pair is pseudoautosomal, with very small X-specific and Y-specific regions. We have prepared an updated genetic and cytogenetic map of the male rainbow trout sex linkage group. Selected sex-linked markers spanning the X chromosome of the female genetic map have been mapped cytogenetically in normal males and genetically in crosses between the OSU female clonal line and four different male clonal lines as well as in outcrosses involving outbred OSU and hybrids between the OSU line and the male clonal lines. The cytogenetic maps of the X and Y chromosomes were very similar to the female genetic map for the X chromosome. Five markers on the male maps are genetically very close to the sex determination locus (SEX), but more widely spaced on the female genetic map and on the cytogenetic map, indicating a large region of suppressed recombination on the Y chromosome surrounding the SEX locus. The male map is greatly extended at the telomere. A BAC clone containing the SCAR (sequence characterized amplified region) Omy-163 marker, which maps close to SEX, was subjected to shotgun sequencing. Two carbonyl reductase genes and a gene homologous to the vertebrate skeletal ryanodine receptor were identified. Carbonyl reductase is a key enzyme involved in production of trout ovarian maturation hormone. This brings the number of type I genes mapped to the sex chromosome to six and has allowed us to identify a region on zebrafish chromosome 10 and medaka chromosome 13 which may be homologous to the distal portion of the long arm of the rainbow trout Y chromosome.

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.

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.

Comparative genome mapping reveals evidence of gene conversion between Sox9 paralogs of rainbow trout (Oncorhynchus mykiss)

Considerable evidence suggests that one genome duplication event preceded the divergence of teleost fishes and a second genome duplication event occurred before the radiation of teleosts of the family Salmonidae. Two Sox9 genes have been isolated from a number of teleosts and are called Sox9a and Sox9b. Two Sox9 gene copies have also been isolated from rainbow trout, a salmonid fish and are called Sox9 and Sox9?2. Previous evaluations of the evolutionary history of rainbow trout Sox9 gene copies using phylogenetic reconstructions of their coding regions indicated that they both belong to the Sox9b clade. In this study, we determine the true evolutionary history of Sox9 gene copies in rainbow trout. We show that the locus referred to as Sox9 in rainbow trout is itself duplicated. Mapping of the duplicated Sox9 gene copies indicates that they are co-orthologs of Sox9b while mapping of Sox9?2 indicates that it is an ortholog of Sox9a. This relationship is supported by phylogenetic reconstruction of Sox9 gene copies in teleosts using their 3? untranslated regions. The conflicting phylogenetic topology of Sox9 genes in rainbow trout indicates the occurrence of gene conversion events between Sox9 and Sox9?2 which is supported by a number of recombination analyses.

Mapping of five candidate sex-determining loci in rainbow trout (Oncorhynchus mykiss)

BACKGROUND

Rainbow trout have an XX/XY genetic mechanism of sex determination where males are the heterogametic sex. The homology of the sex-determining gene (SDG) in medaka to Dmrt1 suggested that SDGs evolve from downstream genes by gene duplication. Orthologous sequences of the major genes of the mammalian sex determination pathway have been reported in the rainbow trout but the map position for the majority of these genes has not been assigned.

RESULTS

Five loci of four candidate genes (Amh, Dax1, Dmrt1 and Sox6) were tested for linkage to the Y chromosome of rainbow trout. We exclude the role of all these loci as candidates for the primary SDG in this species. Sox6i and Sox6ii, duplicated copies of Sox6, mapped to homeologous linkage groups 10 and 18 respectively. Genotyping fishes of the OSU x Arlee mapping family for Sox6i and Sox6ii alleles indicated that Sox6i locus might be deleted in the Arlee lineage.

CONCLUSION

Additional candidate genes should be tested for their linkage to the Y chromosome. Mapping data of duplicated Sox6 loci supports previously suggested homeology between linkage groups 10 and 18. Enrichment of the rainbow trout genomic map with known gene markers allows map comparisons with other salmonids. Mapping of candidate sex-determining loci is important for analyses of potential autosomal modifiers of sex-determination in rainbow trout.

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.

Avian NK activities, cells and receptors

Natural killer (NK) activity has been examined in birds for over 30 years, but evidence that avian NK activity plays crucial roles in disease is only suggestive. In chickens, NK activity is mediated by TCR0 cells in the intestinal epithelium, but elsewhere subsets of alphabeta and gammadelta T cells (NKT cells) may be more important. There are few lectin-like NK receptor genes, located in the genomic region syntenic with the natural killer complex (NKC) as well as the major histocompatibility complex (MHC). In contrast, a huge number of Ig-like receptor genes are located in a region syntenic with the leukocyte receptor complex (LRC).

NK gene complex dynamics and selection for NK cell receptors

Natural killer (NK) cells play important roles in innate defense against infectious agents particularly viruses and also tumors. They mediate their effects through direct cytolysis, release of cytokines and regulation of subsequent adaptive immune responses. NK cells are equipped with sophisticated arrays of inhibitory and activation receptors that regulate their function. In this review we illustrate some of the major evolutionary relationships between NK cell receptors among different animal species and what some of the major mechanisms are that give rise to this diversity in receptor families, including the potential roles of pathogens such as viruses in driving receptor evolution.

A microsatellite genetic map of the turbot (Scophthalmus maximus)

A consensus microsatellite-based linkage map of the turbot (Scophthalmus maximus) was constructed from two unrelated families. The mapping panel was derived from a gynogenetic family of 96 haploid embryos and a biparental diploid family of 85 full-sib progeny with known linkage phase. A total of 242 microsatellites were mapped in 26 linkage groups, six markers remaining unlinked. The consensus map length was 1343.2 cM, with an average distance between markers of 6.5 +/- 0.5 cM. Similar length of female and male maps was evidenced. However, the mean recombination at common intervals throughout the genome revealed significant differences between sexes, approximately 1.6 times higher in the female than in the male. The comparison of turbot microsatellite flanking sequences against the Tetraodon nigroviridis genome revealed 55 significant matches, with a mean length of 102 bp and high sequence similarity (81-100%). The comparative mapping revealed significant syntenic regions among fish species. This study represents the first linkage map in the turbot, one of the most important flatfish in European aquaculture. This map will be suitable for QTL identification of productive traits in this species and for further evolutionary studies in fish and vertebrate species.

Heat shock protein (HSP70) RNA expression differs among rainbow trout (Oncorhynchus mykiss) clonal lines

Heat shock protein 70 (HSP70, 70 kDa) is the most commonly expressed protein in response to thermal stress. The extent of its expression is associated with differences in environmental temperatures. We investigated the heat shock response in red blood cells collected from one-year-old rainbow trout (Oncorhynchus mykiss). Three different clonal lines of rainbow trout (Arlee, OSU and Whale Rock) were utilized, originating from habitats that likely experienced different thermal profile. The relative expression of HSP70 from blood cells treated at 13 degrees C, 16 degrees C, 18 degrees C, 20 degrees C, 22 degrees C, and 24 degrees C was quantified using real-time PCR. The use of red blood cells allows for the control and replication of HSP70 expression patterns. Relative expression of HSP70 differed significantly among the three clonal lines. The Arlee line had the lowest HSP70 response of the three clonal lines at any temperature; indicating a heritable difference. Maximum expression of HSP70 occurred at 22 degrees C in the OSU line and at 24 degrees C in the Whale Rock line. The discovery of variation in HSP70 expression among the clonal lines indicates that future studies to map the genetic control of HSP70 expression differences are possible.

A type I and type II microsatellite linkage map of rainbow trout (Oncorhynchus mykiss) with presumptive coverage of all chromosome arms

Background

The development of large genomic resources has become a prerequisite to elucidate the wide-scale evolution of genomes and the molecular basis of complex traits. Linkage maps represent a first level of integration and utilization of such resources and the primary framework for molecular analyses of quantitative traits. Previously published linkage maps have already outlined the main peculiarities of the rainbow trout meiosis and a correspondance between linkage groups and chromosome arms has been recently established using fluorescent in situ hybridization. The number of chromosome arms which were covered by these maps remained unknown.

Results

We report an updated linkage map based on segregation analysis of more than nine hundred microsatellite markers in two doubled haploid gynogenetic lines. These markers segregated into 31 linkage groups spanning an approximate total map length of 2750 cM. Centromeres were mapped for all the linkage groups using meiogenetic lines. For each of the 31 linkage groups, the meta or acrocentric structure infered from centromere mapping was identical with those recently found with fluorescent in situ hybridization results. The present map is therefore assumed to cover the 52 chromosome arms which constitute the rainbow trout karyotype. Our data confirm the occurrence of a high interference level in this species. Homeologous regions were identified in eleven linkage groups, reflecting the tetraploid nature of the salmonid genome. The data supported the assumption that gene orders are conserved between duplicated groups and that each group is located on a single chromosome arm. Overall, a high congruence with already published rainbow trout linkage maps was found for both gene syntenies and orders.

Conclusion

This new map is likely to cover the whole set of chromosome arms and should provide a useful framework to integrate existing or forthcoming rainbow trout linkage maps and other genomic resources. Since very large numbers of EST containing microsatellite sequences are available in databases, it becomes feasible to construct high-density linkage maps localizing known genes. This will facilitate comparative mapping and, eventually, identification of candidate genes in QTL studies.

Transcriptome profiles of livers and kidneys from three rainbow trout (Oncorhynchus mykiss) clonal lines distinguish stocks from three allopatric populations

Transcriptome profiling is a powerful means of simultaneously identifying large numbers of genes that respond transcriptionally to stimuli of any sort. Whereas individuality at the level of genomic sequence is readily revealed and can be expected to influence transcriptional responses, knowledge of the global transcriptomic consequences of genomic individuality is in its infancy. Appreciation of the inherent variability of biological systems gives us confidence in predicting that no two individuals in any outbred population will respond identically to a stimulus. More critical for comparative studies, even unstimulated transcriptomes will be distinctive for each individual. To assess the confidence with which inferences may be drawn from transcriptome profiling when genetically identical samples can be assured, we examined the unprovoked transcriptomes of hepatic and pronephric (head kidney) tissues in three clonal lines of Rainbow trout (Oncorhynchus mykiss). Clonal individuals derived from three allopatric populations presented transcriptional profiles for both liver and pronephros that were not statistically significantly different within each clonal line; however each clonal line was distinguished by a subset of genes with constitutively different transcript abundance. Among these, immunologically-relevant genes were over-represented, possibly reflecting evolutionarily recent, pathogen-driven genetic sweeps.

Assignment of rainbow trout linkage groups to specific chromosomes

The rainbow trout genetic linkage groups have been assigned to specific chromosomes in the OSU (2N=60) strain using fluorescence in situ hybridization (FISH) with BAC probes containing genes mapped to each linkage group. There was a rough correlation between chromosome size and size of the genetic linkage map in centimorgans for the genetic maps based on recombination from the female parent. Chromosome size and structure have a major impact on the female:male recombination ratio, which is much higher (up to 10:1 near the centromeres) on the larger metacentric chromosomes compared to smaller acrocentric chromosomes. Eighty percent of the BAC clones containing duplicate genes mapped to a single chromosomal location, suggesting that diploidization resulted in substantial divergence of intergenic regions. The BAC clones that hybridized to both duplicate loci were usually located in the distal portion of the chromosome. Duplicate genes were almost always found at a similar location on the chromosome arm of two different chromosome pairs, suggesting that most of the chromosome rearrangements following tetraploidization were centric fusions and did not involve homeologous chromosomes. The set of BACs compiled for this research will be especially useful in construction of genome maps and identification of QTL for important traits in other salmonid fishes.

Comprehensive gene expression profiling following DNA vaccination of rainbow trout against infectious hematopoietic necrosis virus

The DNA vaccine based on the glycoprotein gene of Infectious hematopoietic necrosis virus induces a non-specific anti-viral immune response and long-term specific immunity against IHNV. This study characterized gene expression responses associated with the early anti-viral response. Homozygous rainbow trout were injected intra-muscularly (I.M.) with vector DNA or the IHNV DNA vaccine. Gene expression in muscle tissue (I.M. site) was evaluated using a 16,008 feature salmon cDNA microarray. Eighty different genes were significantly modulated in the vector DNA group while 910 genes were modulated in the IHNV DNA vaccinate group relative to control group. Quantitative reverse-transcriptase PCR was used to examine expression of selected immune genes at the I.M. site and in other secondary tissues. In the localized response (I.M. site), the magnitudes of gene expression changes were much greater in the vaccinate group relative to the vector DNA group for the majority of genes analyzed. At secondary systemic sites (e.g. gill, kidney and spleen), type I IFN-related genes were up-regulated in only the IHNV DNA vaccinated group. The results presented here suggest that the IHNV DNA vaccine induces up-regulation of the type I IFN system across multiple tissues, which is the functional basis of early anti-viral immunity.

Altered self-erythrocyte recognition and destruction in an inbred line of tilapia (Oreochromis aureus)

Carboxyfluorescein diacetate (cFDA)-stained autologous and syngeneic tilapia (Oreochromis aureus) erythrocytes are recognized by effector peripheral blood leukocytes and lysed after a short culture period of 4 h. The hemolysis level was evaluated by measuring the fluorescence of the released cFDA. The degree of lysis of stained target erythrocytes of 60 individuals revealed a trimodal distribution statistically stratified into three groups of low (LR), intermediate (IR), and high (HR) responders. Depletion of the majority of phagocytes from leukocytes lowered the lysis level of HR to that of LR. A highly significant increase of LR cytotoxicity was obtained after the addition of conditioned medium from HR but only in the presence of phagocytes. Genetic analysis of offspring from four crosses (IR x HR, IR x LR, HR x LR, and LR x LR) revealed a quantitative trait locus (QTL) segregating for the level of response linked to markers UNH207 and UNH231 on linkage group 6 of tilapia. Based on segregation analysis of 58 gynogenetic BIU-1 offspring, the distances from the centromere were estimated as 21.5, 11.5, and 9.0 cM for UNH207, UNH231, and the QTL, respectively. It is suggested that 1) self-target recognition and destruction requires both cFDA-altered self-erythrocyte membrane and membrane structures normally present in autologous, syngeneic, and xenogeneic targets; 2) natural cytotoxic cells and/or macrophages are involved in erythrocyte lysis; and 3) the lysis level is codominantly inherited by a QTL segregating on tilapia linkage group 6.

Genomic organisation and expression of the natural killer cell enhancing factor (NKEF) gene in channel catfish, Ictalurus punctatus (Rafinesque)

The gene encoding Natural Killer Cell Enhancing Factor (NKEF) was identified in clones from a gynogenetic channel catfish BAC library. NKEF gene sequence (5.2 kb) was obtained by direct sequencing of BAC clones. The catfish NKEF gene contains one non-coding exon, five coding exons and five introns. A putative TATA box and other transcription factor binding sites were identified in the promoter region. The NKEF amino acid sequence is highly conserved between fish and mammals. Gene expression, measured by real-time quantitative PCR, was detected in all major tissues with the highest level of expression in stomach and heart and lowest levels in gonad and pituitary gland. Catfish NKEF mRNA levels were slightly upregulated 8 and 24 h after injection of lipopolysaccharide. A TAA-repeat microsatellite was identified in a BAC clone containing the NKEF gene, and this locus contained at least 12 alleles in a random-bred catfish population. Multipoint linkage analysis in two reference families placed the NKEF gene on linkage group U18, 5.1 cM from locus IpCG0177 (r=0.05; LOD=45.8), 6.0 cM from a novel immune type receptor, syntenic (40 cM) with T-cell receptor alpha, and not linked with MHC loci.

Comparative genomics of natural killer cell receptor gene clusters

Many receptors on natural killer (NK) cells recognize major histocompatibility complex class I molecules in order to monitor unhealthy tissues, such as cells infected with viruses, and some tumors. Genes encoding families of NK receptors and related sequences are organized into two main clusters in humans: the natural killer complex on Chromosome 12p13.1, which encodes C-type lectin molecules, and the leukocyte receptor complex on Chromosome 19q13.4, which encodes immunoglobulin superfamily molecules. The composition of these gene clusters differs markedly between closely related species, providing evidence for rapid, lineage-specific expansions or contractions of sets of loci. The choice of NK receptor genes is polarized in the two species most studied, mouse and human. In mouse, the C-type lectin-related Ly49 gene family predominates. Conversely, the single Ly49 sequence is a pseudogene in humans, and the immunoglobulin superfamily KIR gene family is extensive. These different gene sets encode proteins that are comparable in function and genetic diversity, even though they have undergone species-specific expansions. Understanding the biological significance of this curious situation may be aided by studying which NK receptor genes are used in other vertebrates, especially in relation to species-specific differences in genes for major histocompatibility complex class I molecules.

Comparative mapping of expressed sequence tags containing microsatellites in rainbow trout (Oncorhynchus mykiss)

BACKGROUND

Comparative genomics, through the integration of genetic maps from species of interest with whole genome sequences of other species, will facilitate the identification of genes affecting phenotypes of interest. The development of microsatellite markers from expressed sequence tags will serve to increase marker densities on current salmonid genetic maps and initiate in silico comparative maps with species whose genomes have been fully sequenced.

RESULTS

Eighty-nine polymorphic microsatellite markers were generated for rainbow trout of which at least 74 amplify in other salmonids. Fifty-five have been associated with functional annotation and 30 were mapped on existing genetic maps. Homologous sequences were identified for 20 of the EST containing microsatellites to identify comparative assignments within the tetraodon, mouse, and/or human genomes.

CONCLUSION

The addition of microsatellite markers constructed from expressed sequence tag data will facilitate the development of high-density genetic maps for rainbow trout and comparative maps with other salmonids and better studied species.