Genomic organization and evolution of the Atlantic salmon hemoglobin repertoire

BACKGROUND

The genomes of salmonids are considered pseudo-tetraploid undergoing reversion to a stable diploid state. Given the genome duplication and extensive biological data available for salmonids, they are excellent model organisms for studying comparative genomics, evolutionary processes, fates of duplicated genes and the genetic and physiological processes associated with complex behavioral phenotypes. The evolution of the tetrapod hemoglobin genes is well studied; however, little is known about the genomic organization and evolution of teleost hemoglobin genes, particularly those of salmonids. The Atlantic salmon serves as a representative salmonid species for genomics studies. Given the well documented role of hemoglobin in adaptation to varied environmental conditions as well as its use as a model protein for evolutionary analyses, an understanding of the genomic structure and organization of the Atlantic salmon α and β hemoglobin genes is of great interest.

RESULTS

We identified four bacterial artificial chromosomes (BACs) comprising two hemoglobin gene clusters spanning the entire α and β hemoglobin gene repertoire of the Atlantic salmon genome. Their chromosomal locations were established using fluorescence in situ hybridization (FISH) analysis and linkage mapping, demonstrating that the two clusters are located on separate chromosomes. The BACs were sequenced and assembled into scaffolds, which were annotated for putatively functional and pseudogenized hemoglobin-like genes. This revealed that the tail-to-tail organization and alternating pattern of the α and β hemoglobin genes are well conserved in both clusters, as well as that the Atlantic salmon genome houses substantially more hemoglobin genes, including non-Bohr β globin genes, than the genomes of other teleosts that have been sequenced.

CONCLUSIONS

We suggest that the most parsimonious evolutionary path leading to the present organization of the Atlantic salmon hemoglobin genes involves the loss of a single hemoglobin gene cluster after the whole genome duplication (WGD) at the base of the teleost radiation but prior to the salmonid-specific WGD, which then produced the duplicated copies seen today. We also propose that the relatively high number of hemoglobin genes as well as the presence of non-Bohr β hemoglobin genes may be due to the dynamic life history of salmon and the diverse environmental conditions that the species encounters.Data deposition: BACs S0155C07 and S0079J05 (fps135): GenBank GQ898924; BACs S0055H05 and S0014B03 (fps1046): GenBank GQ898925.

Risk-based analysis of polychlorinated biphenyl toxicity in harbor seals

Persistent organic pollutants (POPs) have been associated with adverse health effects in marine mammals. However, the complex mixtures to which free-ranging populations are exposed constrain the elucidation of cause-and-effect relationships between specific POPs and the observed health risks. In this study, we 1) assembled data from studies showing polychlorinated biphenyl (PCB)-associated effects on the health of free-ranging harbor seals in the northeastern Pacific Ocean, 2) carried out additional POP analyses on seal samples to broaden the available data on contaminant residues, and 3) derived estimates of individual POPs and their toxic risks. Taken together, these components were used to generate a new toxicity reference value (TRV) for the protection of marine mammal health. In this case study of seals in British Columbia, Canada, and Washington State, USA, PCBs were the single most abundant POP and were correlated with several adverse health effects. PCB exposures consistently exceeded regulatory toxicity thresholds for fish-eating wildlife. Nursing seal pups were at particular risk, reflecting their greatly increased dietary intake of PCBs and their sensitivity to developmental toxicity. Based on the collective evidence obtained, we propose TRVs (consisting of 5% tissue residue concentration and dose) of 1.3 mg/kg lipid weight tissue residue in blubber and 0.05 mg/kg lipid weight tolerable daily intake in prey. Insofar as the TRVs are lower than previously established TRVs and regulatory guidelines, our study highlights the current underestimation of risks associated with PCBs in high-trophic-level wildlife.

Sequencing the genome of the Atlantic salmon (Salmo salar)

The International Collaboration to Sequence the Atlantic Salmon Genome (ICSASG) will produce a genome sequence that identifies and physically maps all genes in the Atlantic salmon genome and acts as a reference sequence for other salmonids.

Evolution of duplicated IgH loci in Atlantic salmon, Salmo salar

BACKGROUND

The Atlantic salmon (Salmo salar) immunoglobulin heavy chain (IgH) locus possesses two parallel IgH isoloci (IGH-A and IGH-B), that are related to the genomic duplication event in the family Salmonidae. These duplicated IgH loci in Atlantic salmon provide a unique opportunity to examine the mechanisms of genome diversity and genome evolution of the IgH loci in vertebrates. In this study, we defined the structure of these loci in Atlantic salmon, and sequenced 24 bacterial artificial chromosome (BAC) clones that were assembled into the IGH-A (1.1 Mb) and IGH-B (0.9 Mb) loci. In addition, over 7,000 cDNA clones from the IgH variable (VH) region have been sequenced and analyzed.

RESULTS

The present study shows that the genomic organization of the duplicated IgH loci in Atlantic salmon differs from that in other teleosts and other vertebrates. The loci possess multiple Cτ genes upstream of the Cμ region, with three of the Cτ genes being functional. Moreover, the duplicated loci possess over 300 VH segments which could be classified into 18 families. This is the largest number of VH families currently defined in any vertebrate. There were significant structural differences between the two loci, indicating that both IGH-A and -B loci have evolved independently in the short time after the recent genome duplication approximately 60 mya.

CONCLUSIONS

Our results indicate that the duplication of the IgH loci in Atlantic salmon significantly contributes to the increased diversity of the antibody repertoire, as compared with the single IgH locus in other vertebrates.

Zonadhesin is essential for species specificity of sperm adhesion to the egg zona pellucida

Interaction of rapidly evolving molecules imparts species specificity to sperm-egg recognition in marine invertebrates, but it is unclear whether comparable interactions occur during fertilization in any vertebrate species. In mammals, the sperm acrosomal protein zonadhesin is a rapidly evolving molecule with species-specific binding activity for the egg zona pellucida (ZP). Here we show using null mice produced by targeted disruption of Zan that zonadhesin confers species specificity to sperm-ZP adhesion. Sperm capacitation selectively exposed a partial von Willebrand D domain of mouse zonadhesin on the surface of living, motile cells. Antibodies to the exposed domain inhibited adhesion of wild-type spermatozoa to the mouse ZP but did not inhibit adhesion of spermatozoa lacking zonadhesin. Zan(-/-) males were fertile, and their spermatozoa readily fertilized mouse eggs in vitro. Remarkably, however, loss of zonadhesin increased adhesion of mouse spermatozoa to pig, cow, and rabbit ZP but not mouse ZP. We conclude that zonadhesin mediates species-specific ZP adhesion, and Zan(-/-) males are fertile because their spermatozoa retain adhesion capability that is not species-specific. Mammalian sperm-ZP adhesion is therefore molecularly robust, and species-specific egg recognition by a protein in the sperm acrosome is conserved between invertebrates and vertebrates, even though the adhesion molecules themselves are unrelated.

High gene expression of inflammatory markers and IL-17A correlates with severity of injection site reactions of Atlantic salmon vaccinated with oil-adjuvanted vaccines

Background

Two decades after the introduction of oil-based vaccines in the control of bacterial and viral diseases in farmed salmonids, the mechanisms of induced side effects manifested as intra-abdominal granulomas remain unresolved. Side effects have been associated with generation of auto-antibodies and autoimmunity but the underlying profile of inflammatory and immune response has not been characterized. This study was undertaken with the aim to elucidate the inflammatory and immune mechanisms of granuloma formation at gene expression level associated with high and low side effect (granuloma) indices.

Groups of Atlantic salmon parr were injected intraperitoneally with oil-adjuvanted vaccines containing either high or low concentrations of Aeromonas salmonicida or Moritella viscosa antigens in order to induce polarized (severe and mild) granulomatous reactions. The established granulomatous reactions were confirmed by gross and histological methods at 3 months post vaccination when responses were known to have matured. The corresponding gene expression patterns in the head kidneys were profiled using salmonid cDNA microarrays followed by validation by real-time quantitative PCR (qPCR). qPCR was also used to examine the expression of additional genes known to be important in the adaptive immune response.

Results

Granulomatous lesions were observed in all vaccinated fish. The presence of severe granulomas was associated with a profile of up-regulation of innate immunity-related genes such as complement factors C1q and C6, mannose binding protein, lysozyme C, C-type lectin receptor, CD209, Cathepsin D, CD63, LECT-2, CC chemokine and metallothionein. In addition, TGF-β (p = 0.001), IL-17A (p = 0.007) and its receptor (IL-17AR) (p = 0.009) representing T_H17 were significantly up-regulated in the group with severe granulomas as were arginase and IgM. None of the genes directly reflective of T_H1 T cell lineage (IFN-γ, CD4) or T_H2 (GATA-3) responses were differentially expressed.

Conclusions

Granulomatous reactions following vaccination with oil-based vaccines in Atlantic salmon have the profile of strong expression of genes related to innate immune responses. The expression of TGF-β, IL-17A and its receptor suggests an involvement of T_H17 T cell lineage and is in conformity with strong infiltration of neutrophils and macrophages into inflamed areas. Arginase upregulation shows that macrophages in these reactions are alternatively activated, indicating also a T_H2-profile. To what extent the expression of IL-17A and its receptor reflects an autoimmune vaccine-based reaction remains elusive but would be in conformity with previous observations of autoimmune reactions in salmon when vaccinated with oil-based vaccines.

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.

Regulation, expression and characterization of aromatase (cyp19b1) transcripts in ovary and testis of rainbow trout (Oncorhynchus mykiss)

Cytochrome P450 aromatase is the key enzyme in the pathway that converts androgens to estrogens. The enzyme functions in the smooth endoplasmic reticulum in a complex with NADPH-cytochrome P450 reductase. In teleost fish, at least two separate loci, cyp19a and cyp19b, encode distinct aromatase isoforms. The activity of cyp19a and cyp19b are predominantly associated with the ovary and brain, respectively, although their expression is not confined solely to these tissues. We found that at least five cyp19b1 transcripts with different 5'-UTRs are generated in the ovary and testis of rainbow trout. Regulation for selection of these variants may be through signals present in exon 2 that recruit alternative splicing factors. Also, binding elements for FOXL2 and SF-1 located within the cyp19b1 intron 1 may influence formation of transcripts that contain the 3'-end of the intron. Another transcript devoid of the exon 2 methionine initiator codon may utilize other downstream in-frame start codons. Less developed stages of ovarian and testicular tissues express only the intron-containing transcripts whereas precocious and more mature gonads express all five cyp19b1 messages. The function of these different 5'-UTRs may be for regulation of cyp19b1 at particular developmental stages or to specify control in distinct gonadal cell-types.

Comprehensive analysis of MHC class I genes from the U-, S-, and Z-lineages in Atlantic salmon

BACKGROUND

We have previously sequenced more than 500 kb of the duplicated MHC class I regions in Atlantic salmon. In the IA region we identified the loci for the MHC class I gene Sasa-UBA in addition to a soluble MHC class I molecule, Sasa-ULA. A pseudolocus for Sasa-UCA was identified in the nonclassical IB region. Both regions contained genes for antigen presentation, as wells as orthologues to other genes residing in the human MHC region.

RESULTS

The genomic localisation of two MHC class I lineages (Z and S) has been resolved. 7 BACs were sequenced using a combination of standard Sanger and 454 sequencing. The new sequence data extended the IA region with 150 kb identifying the location of one Z-lineage locus, ZAA. The IB region was extended with 350 kb including three new Z-lineage loci, ZBA, ZCA and ZDA in addition to a UGA locus. An allelic version of the IB region contained a functional UDA locus in addition to the UCA pseudolocus. Additionally a BAC harbouring two MHC class I genes (UHA) was placed on linkage group 14, while a BAC containing the S-lineage locus SAA (previously known as UAA) was placed on LG10. Gene expression studies showed limited expression range for all class I genes with exception of UBA being dominantly expressed in gut, spleen and gills, and ZAA with high expression in blood.

CONCLUSION

Here we describe the genomic organization of MHC class I loci from the U-, Z-, and S-lineages in Atlantic salmon. Nine of the described class I genes are located in the extension of the duplicated IA and IB regions, while three class I genes are found on two separate linkage groups. The gene organization of the two regions indicates that the IB region is evolving at a different pace than the IA region. Expression profiling, polymorphic content, peptide binding properties and phylogenetic relationship show that Atlantic salmon has only one MHC class Ia gene (UBA), in addition to a multitude of nonclassical MHC class I genes from the U-, S- and Z-lineages.

Grayling (Thymallinae) phylogeny within salmonids: complete mitochondrial DNA sequences of Thymallus arcticus and Thymallus thymallus

The phylogenetic relationships among the three subfamilies (Salmoninae, Coregoninae and Thymallinae) in the Salmonidae have not been addressed extensively at the molecular level. In this study, the whole mitochondrial genomes of two Thymallinae species, Thymallus arcticus and Thymallus thymallus were sequenced, and the published mitochondrial genome sequences of other salmonids were used for Bayesian and maximum-likelihood phylogenetic analyses. These results support an ancestral Coregoninae, branching within the Salmonidae, with Thymallinae as the sister group to Salmoninae.

Assignment of Atlantic salmon (Salmo salar) linkage groups to specific chromosomes: conservation of large syntenic blocks corresponding to whole chromosome arms in rainbow trout (Oncorhynchus mykiss)

Background

Most teleost species, especially freshwater groups such as the Esocidae which are the closest relatives of salmonids, have a karyotype comprising 25 pairs of acrocentric chromosomes and 48–52 chromosome arms. After the common ancestor of salmonids underwent a whole genome duplication, its karyotype would have 100 chromosome arms, and this is reflected in the modal range of 96–104 seen in extant salmonids (e.g., rainbow trout). The Atlantic salmon is an exception among the salmonids as it has 72–74 chromosome arms and its karyotype includes 12 pairs of large acrocentric chromosomes, which appear to be the result of tandem fusions. The purpose of this study was to integrate the Atlantic salmon's linkage map and karyotype and to compare the chromosome map with that of rainbow trout.

Results

The Atlantic salmon genetic linkage groups were assigned to specific chromosomes in the European subspecies using fluorescence in situ hybridization with BAC probes containing genetic markers mapped to each linkage group. The genetic linkage groups were larger for metacentric chromosomes compared to acrocentric chromosomes of similar size. Comparison of the Atlantic salmon chromosome map with that of rainbow trout provides strong evidence for conservation of large syntenic blocks in these species, corresponding to entire chromosome arms in the rainbow trout.

Conclusion

It had been suggested that some of the large acrocentric chromosomes in Atlantic salmon are the result of tandem fusions, and that the small blocks of repetitive DNA in the middle of the arms represent the sites of chromosome fusions. The finding that the chromosomal regions on either side of the blocks of repetitive DNA within the larger acrocentric chromosomes correspond to different rainbow trout chromosome arms provides support for this hypothesis.

The sex determining loci and sex chromosomes in the family salmonidae

Salmonids are descended from a common ancestor that underwent an autotetraploidization event. After a whole genome duplication species could deal with sex determination by deleting one copy of SEX, the sex determining locus, or by recruiting a duplicated transcription factor to become a novel sex determining gene. It is not known which if any of these strategies salmonids adopted, but it appears that they all have primarily a genetic mechanism of sex determination with male heterogamety. The sharing of sex-linked markers on the X and Y chromosomes and the difficulty in identifying Y-specific markers indicate that X and Y chromosomes in salmonids have a large pseudoautosomal region and a small sex determining region. Linkage analyses suggest that either SEX differs in different lineages or else has remained the same and moved by transposition to different chromosomes. The identification of the sex chromosomes in salmonid species has not resolved this issue. It is clear that salmonids are at an early stage in sex chromosome differentiation and therefore provide a wonderful opportunity to study the evolution of sex determination. The availability of a reference salmonid genome sequence would provide an important resource for research in this area.

Identification of a molecular marker for type A spermatogonia by microarray analysis using gonadal cells from pvasa-GFP transgenic rainbow trout (Oncorhynchus mykiss)

The spermatogonia of fish can be classified as being either undifferentiated type A spermatogonia or differentiated type B spermatogonia. Although type A spermatogonia, which contain spermatogonial stem cells, have been demonstrated to be a suitable material for germ cell transplantation, no molecular markers for distinguishing between type A and type B spermatogonia in fish have been developed to date. We therefore sought to develop a molecular marker for type A spermatogonia in rainbow trout. Using GFP-dependent flow cytometry (FCM), enriched fractions of type A and type B spermatogonia, testicular somatic cells, and primordial germ cells were prepared from rainbow trout possessing the green fluorescent protein (GFP) gene driven by trout vasa regulatory regions (pvasa-GFP rainbow trout). The gene-expression profiles of each cell fraction were then compared with a microarray containing cDNAs representing 16,006 genes from several salmonid species. Genes exhibiting high expression for type A spermatogonia relative to above-mentioned other types of gonadal cells were identified and subjected to RT-PCR and quatitative PCR analysis. Since only the rainbow trout notch1 homologue showed significantly high expression in the type A spermatogonia-enriched fraction, we propose that notch1 may be a useful molecular marker for type A spermatogonia. The combination of GFP-dependent FCM and microarray analysis of pvasa-GFP rainbow trout can therefore be applied to the identification of potentially useful molecular markers of germ cells in fish.

Genomic organization and evolution of the vomeronasal type 2 receptor-like (OlfC) gene clusters in Atlantic salmon, Salmo salar

There are three major multigene superfamilies of olfactory receptors (OR, V1R, and V2R) in mammals. The ORs are expressed in the main olfactory organ, whereas the V1Rs and V2Rs are located in the vomeronasal organ. Fish only possess one olfactory organ in each nasal cavity, the olfactory rosette; therefore, it has been proposed that their V2R-like genes be classified as olfactory C family G protein-coupled receptors (OlfC). There are large variations in the sizes of OR gene repertoires. Previous studies have shown that fish have between 12 and 46 functional V2R-like genes, whereas humans have lost all functional V2Rs, and frog sp. have more than 240. Pseudogenization of V2R genes is a prevalent event across species. In the mouse and frog genomes, there are approximately double the number of pseudogenes compared with functional genes. An oligonucleotide probe was designed from a conserved sequence from four Atlantic salmon OlfC genes and used to screen the Atlantic salmon bacterial artificial chromosome (BAC) library. Hybridization-positive BACs were matched to fingerprint contigs, and representative BACs were shotgun cloned and sequenced. We identified 55 OlfC genes. Twenty-nine of the OlfC genes are classified as putatively functional genes and 26 as pseudogenes. The OlfC genes are found in two genomic clusters on chromosomes 9 and 20. Phylogenetic analysis revealed that the OlfC genes could be divided into 10 subfamilies, with nine of these subfamilies corresponding to subfamilies found in other teleosts and one being salmon specific. There is also a large expansion in the number of OlfC genes in one subfamily in Atlantic salmon. Subfamily gene expansions have been identified in other teleosts, and these differences in gene number reflect species-specific evolutionary requirements for olfaction. Total RNA was isolated from the olfactory epithelium and other tissues from a presmolt to examine the expression of the odorant genes. Several of the putative OlfC genes that we identified are expressed only in the olfactory epithelium, consistent with these genes encoding odorant receptors.

A salmonid EST genomic study: genes, duplications, phylogeny and microarrays

Background

Salmonids are of interest because of their relatively recent genome duplication, and their extensive use in wild fisheries and aquaculture. A comprehensive gene list and a comparison of genes in some of the different species provide valuable genomic information for one of the most widely studied groups of fish.

Results

298,304 expressed sequence tags (ESTs) from Atlantic salmon (69% of the total), 11,664 chinook, 10,813 sockeye, 10,051 brook trout, 10,975 grayling, 8,630 lake whitefish, and 3,624 northern pike ESTs were obtained in this study and have been deposited into the public databases. Contigs were built and putative full-length Atlantic salmon clones have been identified. A database containing ESTs, assemblies, consensus sequences, open reading frames, gene predictions and putative annotation is available. The overall similarity between Atlantic salmon ESTs and those of rainbow trout, chinook, sockeye, brook trout, grayling, lake whitefish, northern pike and rainbow smelt is 93.4, 94.2, 94.6, 94.4, 92.5, 91.7, 89.6, and 86.2% respectively. An analysis of 78 transcript sets show Salmo as a sister group to Oncorhynchus and Salvelinus within Salmoninae, and Thymallinae as a sister group to Salmoninae and Coregoninae within Salmonidae. Extensive gene duplication is consistent with a genome duplication in the common ancestor of salmonids. Using all of the available EST data, a new expanded salmonid cDNA microarray of 32,000 features was created. Cross-species hybridizations to this cDNA microarray indicate that this resource will be useful for studies of all 68 salmonid species.

Conclusion

An extensive collection and analysis of salmonid RNA putative transcripts indicate that Pacific salmon, Atlantic salmon and charr are 94–96% similar while the more distant whitefish, grayling, pike and smelt are 93, 92, 89 and 86% similar to salmon. The salmonid transcriptome reveals a complex history of gene duplication that is consistent with an ancestral salmonid genome duplication hypothesis. Genome resources, including a new 32 K microarray, provide valuable new tools to study salmonids.

Isolation, characterization and comparison of Atlantic and Chinook salmon growth hormone 1 and 2

Background

Growth hormone (GH) is an important regulator of skeletal growth, as well as other adapted processes in salmonids. The GH gene (gh) in salmonids is represented by duplicated, non-allelic isoforms designated as gh1 and gh2. We have isolated and characterized gh-containing bacterial artificial chromosomes (BACs) of both Atlantic and Chinook salmon (Salmo salar and Oncorhynchus tshawytscha) in order to further elucidate our understanding of the conservation and regulation of these loci.

Results

BACs containing gh1 and gh2 from both Atlantic and Chinook salmon were assembled, annotated, and compared to each other in their coding, intronic, regulatory, and flanking regions. These BACs also contain the genes for skeletal muscle sodium channel oriented in the same direction. The sequences of the genes for interferon alpha-1, myosin alkali light chain and microtubule associated protein Tau were also identified, and found in opposite orientations relative to gh1 and gh2. Viability of each of these genes was examined by PCR. We show that transposon insertions have occurred differently in the promoters of gh, within and between each species. Other differences within the promoters and intronic and 3'-flanking regions of the four gh genes provide evidence that they have distinct regulatory modes and possibly act to function differently and/or during different times of salmonid development.

Conclusion

A core proximal promoter for transcription of both gh1 and gh2 is conserved between the two species of salmon. Nevertheless, transposon integration and regulatory element differences do exist between the promoters of gh1 and gh2. Additionally, organization of transposon families into the BACs containing gh1 and for the BACs containing gh2, are very similar within orthologous regions, but much less clear conservation is apparent in comparisons between the gh1- and gh2-containing paralogous BACs for the two fish species. This is consistent with the hypothesis that a burst of transposition activity occurred during the speciation events which led to Atlantic and Pacific salmon. The Chinook and other Oncorhynchus GH1s are strikingly different in comparison to the other GHs and this change is not apparent in the surrounding non-coding sequences.

EST and mitochondrial DNA sequences support a distinct Pacific form of salmon louse, Lepeophtheirus salmonis

Nuclear deoxyribonucleic acid sequences from approximately 15,000 salmon louse expressed sequence tags (ESTs), the complete mitochondrial genome (16,148bp) of salmon louse, and 16S ribosomal ribonucleic acid (rRNA) and cytochrome oxidase subunit I (COI) genes from 68 salmon lice collected from Japan, Alaska, and western Canada support a Pacific lineage of Lepeophtheirus salmonis that is distinct from that occurring in the Atlantic Ocean. On average, nuclear genes are 3.2% different, the complete mitochondrial genome is 7.1% different, and 16S rRNA and COI genes are 4.2% and 6.1% different, respectively. Reduced genetic diversity within the Pacific form of L. salmonis is consistent with an introduction into the Pacific from the Atlantic Ocean. The level of divergence is consistent with the hypothesis that the Pacific form of L. salmonis coevolved with Pacific salmon (Onchorhynchus spp.) and the Atlantic form coevolved with Atlantic salmonids (Salmo spp.) independently for the last 2.5-11 million years. The level of genetic divergence coincides with the opportunity for migration of fish between the Atlantic and Pacific Ocean basins via the Arctic Ocean with the opening of the Bering Strait, approximately 5 million years ago. The genetic differences may help explain apparent differences in pathogenicity and environmental sensitivity documented for the Atlantic and Pacific forms of L. salmonis.

Assessing the feasibility of GS FLX Pyrosequencing for sequencing the Atlantic salmon genome

Background

With a whole genome duplication event and wealth of biological data, salmonids are excellent model organisms for studying evolutionary processes, fates of duplicated genes and genetic and physiological processes associated with complex behavioral phenotypes. It is surprising therefore, that no salmonid genome has been sequenced. Atlantic salmon (Salmo salar) is a good representative salmonid for sequencing given its importance in aquaculture and the genomic resources available. However, the size and complexity of the genome combined with the lack of a sequenced reference genome from a closely related fish makes assembly challenging. Given the cost and time limitations of Sanger sequencing as well as recent improvements to next generation sequencing technologies, we examined the feasibility of using the Genome Sequencer (GS) FLX pyrosequencing system to obtain the sequence of a salmonid genome. Eight pooled BACs belonging to a minimum tiling path covering ~1 Mb of the Atlantic salmon genome were sequenced by GS FLX shotgun and Long Paired End sequencing and compared with a ninth BAC sequenced by Sanger sequencing of a shotgun library.

Results

An initial assembly using only GS FLX shotgun sequences (average read length 248.5 bp) with ~30× coverage allowed gene identification, but was incomplete even when 126 Sanger-generated BAC-end sequences (~0.09× coverage) were incorporated. The addition of paired end sequencing reads (additional ~26× coverage) produced a final assembly comprising 175 contigs assembled into four scaffolds with 171 gaps. Sanger sequencing of the ninth BAC (~10.5× coverage) produced nine contigs and two scaffolds. The number of scaffolds produced by the GS FLX assembly was comparable to Sanger-generated sequencing; however, the number of gaps was much higher in the GS FLX assembly.

Conclusion

These results represent the first use of GS FLX paired end reads for de novo sequence assembly. Our data demonstrated that this improved the GS FLX assemblies; however, with respect to de novo sequencing of complex genomes, the GS FLX technology is limited to gene mining and establishing a set of ordered sequence contigs. Currently, for a salmonid reference sequence, it appears that a substantial portion of sequencing should be done using Sanger technology.

A linkage map of the Atlantic salmon (Salmo salar) based on EST-derived SNP markers

Background

The Atlantic salmon is a species of commercial and ecological significance. Like other salmonids, the species displays residual tetrasomy and a large difference in recombination rate between sexes. Linkage maps with full genome coverage, containing both type I and type II markers, are needed for progress in genomics. Furthermore, it is important to estimate levels of linkage disequilibrium (LD) in the species. In this study, we developed several hundred single nucleotide polymorphism (SNP) markers for the Atlantic salmon, and constructed male and female linkage maps containing SNP and microsatellite markers. We also investigated further the distribution of male and female recombination events across the genome, and estimated levels of LD between pairs of markers.

Results

The male map had 29 linkage groups and was 390 cM long. The female map had 30 linkage groups as was 1983 cM long. In total, the maps contained 138 microsatellite markers and 304 SNPs located within genes, most of which were successfully annotated. The ratio of male to female recombination events was either close to zero or very large, indicating that there is little overlap between regions in which male and female crossovers occur. The female map is likely to have close to full genome coverage, while the majority of male linkage groups probably lack markers in telomeric regions where male recombination events occur. Levels of r² increased with decreasing inter-marker distance in a bimodal fashion; increasing slowly from ~60 cM, and more rapidly more from ~12 cM. Long-ranging LD may be consequence of recent admixture in the population, the population being a 'synthetic' breeding population with contributions from several distinct rivers. Levels of r² dropped to half its maximum value (above baseline) within 15 cM, and were higher than 0.2 above baseline for unlinked markers ('useful LD') at inter-marker distances less than 5 cM.

Conclusion

The linkage map presented here is an important resource for genetic, comparative, and physical mapping of the Atlantic salmon. The female map is likely to have a map coverage that is not far from complete, whereas the male map length is likely to be significantly shorter than the true map, due to suboptimal marker coverage in the apparently small physical regions where male crossovers occur. 'Useful LD' was found at inter-marker distances less than 5 cM.

Genomic organization and characterization of two vomeronasal 1 receptor-like genes (ora1 and ora2) in Atlantic salmon Salmo salar

Olfactory receptors are encoded by three large multigene superfamilies (OR, V1R and V2R) in mammals. Fish do not possess a vomeronasal system; therefore, it has been proposed that their V1R-like genes be classified as olfactory receptors related to class A G protein-coupled receptors (ora). Unlike mammalian genomes, which contain more than a hundred V1R genes, the five species of teleost fish that have been investigated to date appear to have six ora genes (ora1-6) except for pufferfish that have lost ora1. The common ancestor of salmonid fishes is purported to have undergone a whole genome duplication. As salmonids have a life history that requires the use of olfactory cues to navigate back to their natal habitats to spawn, we set out to determine if ora1 or ora2 is duplicated in a representative species, Atlantic salmon (Salmo salar). We used an oligonucleotide probe designed from a conserved sequence of several teleost ora2 genes to screen an Atlantic salmon BAC library (CHORI-214). Hybridization-positive BACs belonged to a single fingerprint contig of the Atlantic salmon physical map. All were also positive for ora2 by PCR. One of these BACs was chosen for further study, and shotgun sequencing of this BAC identified two V1R-like genes, ora1 and ora2, that are in a head-to-head conformation as is seen in some other teleosts. The gene products, ora1 and ora2, are highly conserved among teleosts. We only found evidence for a single ora1-2 locus in the Atlantic salmon genome, which was mapped to linkage group 6. Fluorescent in situ hybridization (FISH) analysis placed ora1-2 on chromosome 12. Conserved synteny was found surrounding the ora1 and ora2 genes in Atlantic salmon, medaka and three-spined stickleback, but not zebrafish.

Effects of diesel on survival, growth, and gene expression in rainbow trout (Oncorhynchus mykiss) fry

Diesel spills are all too frequent disturbances of freshwater ecosystems, largely as a result of the quantities transported and consumed. Assessing the risk that such events may pose to aquatic life remains a difficult process, because of the complexity of this hydrocarbon mixture and our limited knowledge of its toxicity. A diesel spike experiment with rainbow trout (Oncorhynchus mykiss) fry was carried out to fill this knowledge gap. Survival, growth, and gene expression changes were assessed and toxicity thresholds were determined. Whereas the biological end points were consistent in the determination of (sub)lethal doses, microarrays supplied additional information on the mechanism of toxicity (oxygen deprivation) and potential long-term effects (feminization, immune system alterations) of diesel exposure on salmonids. Hemoglobins, prostaglandins, cytochromes, and gluthathion-S-transferases were among the molecular biomarkers proposed for use in future risk assessments based on microarray results. By bridging traditional toxicity testing with recent microarray technologies, this study shows the potential of genomics tools in ecotoxicity studies as well as industrial applications, including risk assessment, in the near future.

Functional adaptive diversity of the Atlantic salmon T-cell receptor gamma locus

The Atlantic salmon (Salmo salar) T-cell receptor gamma (TCRgamma) gene has extensive diversity in its capacity for antigen recognition due to the V-J-C gene segments recombinational possibilities, and N-region diversity. This is the first report completely characterizing and annotating a TCRgamma gene locus in teleosts. We identified two different TCRgamma loci in Atlantic salmon. The first locus, TCRgamma 1, spans 260 kbp and contains four tandemly repeated clusters each of which consists of 1-4 variable (V) segments, 1-2 sets of a joining (J) segment and a constant (C) region. In total, 10 V segments, 5 J segments and 5 C regions were found in locus 1. In the second locus, TCRgamma 2, a single non-expressed V-J-C cluster was found. Surprisingly, the Atlantic salmon TCRgamma loci have a larger number of C regions relative to other teleosts, mouse and human. In addition, each Atlantic salmon TCRgamma C region has a different connecting peptide region that may result in distinct cellular responses. Expression data confirm the diverse repertoire found at the genomic level. At least 6 out of 7 functional V segments, all 5 J segments and all 5 C regions found in TCRgamma 1 were identified in TCRgamma transcription analysis. The identification and characterization of the functional TCRgamma 1 along with the associated TCR alpha/delta locus suggest that Atlantic salmon have a functional gammadelta T-cell immune component.

Striking antigen recognition diversity in the Atlantic salmon T-cell receptor alpha/delta locus

The complete TCR alpha/delta locus of Atlantic salmon (Salmo salar) has been characterized and annotated. In the 900 kb TCR alpha/delta locus, 292 Valpha/delta segments and 123 Jalpha/delta segments were identified. Of these, 128 Valpha/delta, 113 Jalpha, and a Jdelta segment appeared to be functional as they lacked frame shifts or stop codons. This represents the largest repertoire of Valpha/delta and Jalpha segments of any organism to date. The 128 functional Valpha/delta segments could be grouped into 29 subgroups based upon 70% nucleotide similarity. Expression data confirmed the usage of the diverse repertoire found at the genomic level. At least 99 Valpha, 13 Vdelta 86 Jalpha, 1 Jdelta, and 2 Ddelta segments were used in TCR alpha or delta transcription, and 652 unique genes were identified from a sample of 759 TCRalpha cDNA clones. Cumulatively, the genomic and expression data suggest that the Atlantic salmon T-cell receptor has enormous capacity to recognize a wide diversity of antigens.