Molecular markers for variation in spawning date in a hatchery population of rainbow trout (Oncorhynchus mykiss).

We examined the distribution of alleles at 63 microsatellite loci distributed across 29 linkage groups in broodstock females from a commercial population of rainbow trout spawning on different dates throughout the season (August to January). A total of 368 females, 184 and 117 females from each of the tail-ends of the spawning distribution and a subsample of 67 females spawning in the middle, were used to detect marker-trait associations. Twenty-one loci in a subset of genomic regions (RT-5, 7, 8, 10, 12, 14, 15, 22, 23, 24, 25, 29, 30, and 31) were significantly associated with variation in spawning date. Many of these markers localize to regions with known spawning date quantitative trait loci based on previous studies. An individual assignment analysis was used to test how well the molecular data could be used to assign individuals to their correct spawning group, and markers were given a ranking reflecting their contribution to the accuracy of assignment. The top 15 ranked markers were successful at assigning the majority of females to the correct spawning group based on genotype with an average accuracy of 76 %. The most likely genes that could contribute to these differences in spawning date are discussed. Together, these data indicate that the loci could be incorporated into a selection index with phenotype data to increase the accuracy of selection for spawning date.

Comparative analysis of genetic parameters and quantitative trait loci for growth traits in Fraser strain Arctic charr (Salvelinus alpinus) reared in freshwater and brackish water environments.

To determine the potential for genetic improvement in Fraser strain Arctic charr (AC, Salvelinus alpinus), we calculated genetic parameters for BW and condition factor (K) and tested if previously identified QTL for these traits were detectable across a commercial broodstock reared in both freshwater (FRW) and brackish water (BRW). Individuals from 30 full-sib families were reared up to 29 mo of age in FRW and BRW tanks at a commercial facility. Heritability for BW and K was moderate in FRW (0.29 to 0.38) but lower in BRW (0.14 to 0.17). Genetic correlations for BW across environments were positive and moderate (0.33 to 0.67); however, equivalent K correlations were very weak (0.24 to 0.37). We identified a single BW QTL with experimentwide effects on linkage group AC-8, 4 BW QTL (AC-4, -13, -14, and -19), and 3 K QTL (AC-4, -5, and -20) with chromosomewide effects across families. Notably, the QTL on AC-8 had significant effects with BW at 3 out of 4 sampling dates in FRW and had significant allelic phase disequilibrium with BW across families, suggesting a tight coupling of the marker region to the QTL in this population. Body weight QTL were identified on AC-4 in both FRW and BRW environments and AC-4 was the only linkage group with a detectable QTL for both K and BW. Modest consistency of some QTL effects as well as moderate heritability in both environments suggests that there is some potential for genetic improvement of growth in this species even though gene × environment interactions are high.

Whole genome duplication: challenges and considerations associated with sequence orthology assignment in Salmoninae.

To illustrate some of the challenges and considerations in assigning correct orthology necessary for any comparative genomic investigation among salmonids, sequence data from the non-coding regions of different chromosomes in three members of the subfamily Salmoninae, rainbow trout Oncorhynchus mykiss, Atlantic salmon Salmo salar and Arctic charr Salvelinus alpinus, were compared. By analysing c. 55 distinct loci, corresponding to c. 142 kbp sequence information per species, 18 duplicated patterns representative of the two sequential rounds of teleost-specific whole genome duplications (i.e. 3R and 4R WGD) were identified. Sequence similarities between the 4R paralogues were c. 90%, which was slightly lower than those of the 4R orthologues and c. 60% for the 3R products. Through careful examination of the sequence data, however, only 14 loci could reliably be assigned as true orthologues. Locus-specific trees were constructed through maximum parsimony, maximum likelihood and neighbour-joining methods and were rooted using the information from a close relative, lake whitefish Coregonus clupeaformis. All approaches generated congruent trees supporting the {Coregonus [Salmo (Oncorhynchus, Salvelinus)]} topology. The general phenotypic characteristics of sequences, however, were highly suggestive of the basal position of Oncorhynchus, raising the hypothesis of an accelerated rate of nucleotide evolution in this species.

The genetic architecture of embryonic developmental rate and genetic covariation with age at maturation in rainbow trout Oncorhynchus mykiss.

The genetic architecture underlying variation in embryonic developmental rate (DR) and genetic covariation with age of maturation (MAT) was investigated in rainbow trout Oncorhynchus mykiss. Highly significant additive parental effects and more limited evidence of epistatic effects on progeny hatching time were detected in three diallel sets of families. Genome scans with an average of 142 microsatellite loci from all 29 linkage groups in two families detected significant quantitative trait loci (QTL) for developmental rate on RT-8 and RT-30 with genome-wide and chromosome-wide effects, respectively. The QTL on linkage group RT-8 explained 23·7% of the phenotypic variation and supports results from previous studies. The co-localization of QTL for both DR and MAT to several linkage groups and the observation that alleles associated with faster developmental rate were found significantly more often in early maturing rather than typical and later maturing male ancestors supports the hypothesis of genetic covariation between DR and MAT. The maturation background and schedule of additional sires, however, did not have a consistent association with their progeny hatching times, suggesting that other genetic, environmental and physiological effects contribute to variation in these life-history traits.

Genetic mapping of the ornithine decarboxylase (odc) gene complex in rainbow trout (Oncorhynchus mykiss).

Ornithine decarboxylase (odc) and its associated gene complex, including antizymes of ODC and inhibitors of the antizymes, play a key role in regulating the overall polyamine levels (i.e., putrescine) in cells. Polyamine production levels are intricately coupled to mitotic and cellular turnover rates. Hence, these genes may be important candidates of growth regulation in vertebrates, if their chromosomal locations coincide with known quantitative trait locus (QTL) regions influencing growth traits. Here we report the genetic mapping of multiple duplicated forms of genes within this complex to previously known life-history and growth QTL regions in rainbow trout (Oncorhynchus mykiss). Specifically, duplicated copies of the 2 antizyme genes map to Om-8/9/24, while antizyme inhibitor and odc duplicates map to Om-27/31 and Om-14/25, respectively. All sets of paralogous mapping locations correspond to ancestrally duplicated synteny regions within the genome of rainbow trout (i.e., ancestral A, B, GH/I and M linkage groups).

Genomic organization of the IGF1, IGF2, MYF5, MYF6 and GRF/PACAP genes across Salmoninae genera.

Whole-genome duplication in the ancient ray-finned fish and subsequent tetraploidization in the ancestor to the salmonids have complicated genomic and candidate gene studies in these organisms as many genes with multiple copies are present throughout their genomes. In an attempt to identify genes with a potential influence on growth and development, we investigated the genomic positions of insulin-like growth factors 1 and 2 (IGF1, IGF2), myogenic factors 5 and 6 (MYF5, MYF6) and growth hormone-releasing factor/pituitary adenylate cyclase-activating polypeptide (GRF/PACAP) in three salmonid species: rainbow trout (Oncorhynchus mykiss), Atlantic salmon (Salmo salar) and Arctic charr (Salvelinus alpinus). Our results suggest a tight association between the IGF1, MYF5 and MYF6 genes in all three species. We further localized the duplicated copies of IGF1 to the homeologous linkage groups RT-7/15 in rainbow trout and AC-3/24 in Arctic charr, and the two copies of MYF6 to homeologous linkage groups AS-22/24 in Atlantic salmon. Localization of GRF/PACAP to RT-7, AS-31 and AC-27 and IGF2 to RT-27, AS-2 and AC-4 in rainbow trout, Atlantic salmon and Arctic charr respectively is consistent with previously reported homologies among these chromosomal segments identified using other genetic markers. However, localization of the second copy of GRF/PACAP to RT-19 and AC-14 and the duplicated copy of IGF2 to AC-19 suggest a possible new homology/homeology between these chromosomes. These results might also be an indication of a more ancient polyploidization event that occurred deep in the ray-finned fish lineage.

Genome organization of glutamine synthetase genes in rainbow trout (Oncorhynchus mykiss).

Unlike mammals, bony fish appear to possess multiple genes encoding glutamine synthetase (GS), the nitrogen metabolism enzyme responsible for the conversion of glutamate and ammonia into glutamine at the expense of ATP. This study reports on the development of genetic markers for each of the four isoforms identified thus far in rainbow trout (Oncorhynchus mykiss) and their genome localization by linkage mapping. We found that genes coding for GS01, GS02, GS03, and GS04 map to four different linkage groups in the trout genome, namely RT-24, RT-23, RT-08, and RT-13, respectively. Linkage groups RT-23 and RT-13 appear to represent distinct chromosomes sharing duplicated marker regions, which lends further support to the previous suggestion that GS02 and GS04 may be duplicate gene copies that evolved from a whole-genome duplication in the trout ancestor. In contrast, there is at present no further evidence that RT-24 and RT-08 share ancestrally homologous segments and additional genomic studies will be needed to clarify the evolutionary origin of genes coding for GS01 and GS03.

The candidate gene, Clock, localizes to a strong spawning time quantitative trait locus region in rainbow trout.

We applied a candidate gene mapping approach to an existing quantitative trait loci (QTL) data set for spawning date in rainbow trout (Oncorynchus mykiss) to ascertain whether these genes could potentially account for any observed QTL effects. Several genes were chosen for their known or suspected roles in reproduction, circadian, or circannual timing, including salmon-type gonadotropin-releasing hormone 3A and 3B (GnRH3A and GnRH3B), Clock, Period1, and arylalkylamine N-acetlytransferase-1 and -2 (AANAT-1 and AANAT-2). Genes were sequenced, and polymorphisms were identified in parents of two rainbow trout mapping families, one of which was used previously to detect spawn timing QTL. Interval mapping was used to identify associations between genetic markers and spawning date effects. Using a genetic map that was updated with 574 genetic markers (775 total), we found evidence for 11 significant or suggestive QTL regions. Most QTL were only localized within one of the parents; however, a strong QTL region was identified in both female and male parents on linkage group RT-8 that explained 20% and 50% of trait variance, respectively. The Clock gene mapped to this region. Period1 mapped to a region in the female parent associated with a marginal effect (P = .056) on spawn timing. Other candidate genes were not associated with significant QTL effects.

Identification of the sex-determining locus of Atlantic salmon (Salmo salar) on chromosome 2.

We have integrated data from linkage mapping, physical mapping and karyotyping to gain a better understanding of the sex-determining locus, SEX, in Atlantic salmon (Salmo salar). SEX has been mapped to Atlantic salmon linkage group 1 (ASL1) and is associated with several microsatellite markers. We have used probes designed from the flanking regions of these sex-linked microsatellite markers to screen a bacterial artificial chromosome (BAC) library, representing an 11.7x coverage of the Atlantic salmon genome, which has been HindIII fingerprinted and assembled into contigs. BACs containing sex-linked microsatellites and their related contigs have been identified and representative BACs have been placed on the Atlantic salmon chromosomes by fluorescent in situ hybridization (FISH). This identified chromosome 2, a large metacentric, as the sex chromosome. By positioning several BACs on this chromosome by FISH, it was possible to orient ASL1 with respect to chromosome 2. The region containing SEX appears to lie on the long arm between marker Ssa202DU and a region of heterochromatin identified by DAPI staining. BAC end-sequencing of clones within sex-linked contigs revealed five hitherto unmapped genes along the sex chromosome. We are using an in silico approach coupled with physical probing of the BAC library to extend the BAC contigs to provide a physical map of ASL1, with a view to sequencing chromosome 2 and, in the process, identifying the sex-determining gene.

Characterization and comparison of microsatellites derived from repeat-enriched libraries and expressed sequence tags.

The construction of high-density linkage maps for use in identifying loci underlying important traits requires the development of large numbers of polymorphic genetic markers spanning the entire genome at regularly spaced intervals. As part of our efforts to develop markers for rainbow trout (Oncorhynchus mykiss), we performed a comparison of allelic variation between microsatellite markers developed from expressed sequence tag (EST) data and anonymous markers identified from repeat-enriched libraries constructed from genomic DNA. A subset of 70 markers (37 from EST databases and 33 from repeat enriched libraries) was characterized with respect to polymorphism information content (PIC), number of alleles, repeat number, locus duplication within the genome and ability to amplify in other salmonid species. Higher PIC was detected in dinucleotide microsatellites derived from ESTs than anonymous markers (72.7% vs. 54.0%). In contrast, dinucleotide repeat numbers were higher for anonymous microsatellites than for EST derived microsatellites (27.4 vs.18.1). A higher rate of cross-species amplification was observed for EST microsatellites. Approximately half of each marker type was duplicated within the genome. Unlike single-copy markers, amplification of duplicated microsatellites in other salmonids was not correlated to phylogenetic distance. Genomic microsatellites proved more useful than EST derived microsatellites in discriminating among the salmonids. In total, 428 microsatellite markers were developed in this study for mapping and population genetic studies in rainbow trout.

Sex-linked quantitative trait loci for thermotolerance and length in the rainbow trout.

We hypothesized that correlation between growth traits and upper thermal tolerance (UTT) in rainbow trout (Oncorhynchus mykiss) might be explained by quantitative trait loci (QTL) localized to the same linkage groups. Microsatellites on three autosomal linkage groups carrying UTT QTL in rainbow trout were tested for associations with fork length (FL) and condition factor (K) in half-sib families of outbred rainbow trout and in backcrosses of trout lines selected on UTT. Additionally, we used a sex-linked microsatellite (OmyFGT19TUF) to test for marker-trait associations at the sex chromosomes. The sex-linked marker OmyFGT19TUF was significantly associated with FL and UTT, accounting for up to 9.6% and 9.7% of variance in these traits, respectively. Male advantages in FL (and, to a lesser extent, UTT) relative to their female sibs were dependent on the origin of the Y chromosome and thus varied among grandsire lines. However, males had higher K in a manner unrelated to Y chromosomal origin, suggesting a partially sex-limited expression of this trait. Omy325UoG was significantly associated with K in one of the outbred half-sib families, but no other significant autosomal marker-trait associations were detected. Our findings illustrate minor evidence that correlation between UTT and FL is partially determined by one or more sex-chromosomal QTL.

QTL for body weight and condition factor in Atlantic salmon (Salmo salar): comparative analysis with rainbow trout (Oncorhynchus mykiss) and Arctic charr (Salvelinus alpinus).

Genotypes at 91 microsatellite loci in three full-sib families were used to search for QTL affecting body weight (BW) and condition factor in North American Atlantic salmon (Salmo salar). More than one informative marker was identified on 16-18 linkage groups in each family, allowing at least one chromosomal interval to be analyzed per linkage group. Two significant QTL for BW on linkage groups AS-8 and AS-11, and four significant QTL for condition factor on linkage groups AS-2, AS-5, AS-11, and AS-14 were identified. QTL for both BW and condition factor were located on linkage groups AS-1, 6, 8, 11, and 14 when considering both significant and suggestive QTL effects. The largest QTL effects for BW (AS-8) and for condition factor (AS-14) accounted for 20.1 and 24.9% of the trait variation, respectively. Three of the QTL for BW occur on linkage groups where similar effects have been detected on the homologous regions in either rainbow trout (Oncorhynchus mykiss) or Arctic charr (Salvelinus alpinus).

Linkage arrangement of Na,K-ATPase genes in the tetraploid-derived genome of the rainbow trout (Oncorhynchus mykiss).

As part of our efforts to characterize Na,K-ATPase isoforms in salmonid fish, we investigated the linkage arrangement of genes coding for the alpha and beta-subunits of the enzyme complex in the tetraploid-derived genome of the rainbow trout (Oncorhynchus mykiss). Genetic markers were developed from four of five previously characterized alpha-subunit isoforms (alpha1b, alpha1c, alpha2 and alpha3) and four expressed sequence tags derived from yet undescribed beta-subunit isoforms (beta1a, beta1b, beta3a and beta3b). Sex-specific linkage analysis of polymorphic loci in a reference meiotic panel revealed that Na,K-ATPase genes are generally dispersed throughout the rainbow trout genome. A notable exception was the colocalization of two alpha-subunit genes and one beta-subunit gene on linkage group RT-12, which may thus share a conserved orthologous segment with linkage group 1 in zebrafish (Danio rerio). Consistent with previously reported homeologous relationships among the chromosomes of the rainbow trout, primers designed from the alpha3-isoform detected a pair of duplicated genes on linkage groups RT-27 and RT-31. Similarly, the evolutionary conservation of homeologous regions on linkage groups RT-12 and RT-16 was further supported by the map localization of gene duplicates for the beta1b isoform. The detection of homeologs within each gene family also raises the possibility that novel isoforms may be discovered as functional duplicates.

Population genetic structure of Arctic charr, Salvelinus alpinus from northwest Europe on large and small spatial scales.

To examine the population genetic structure of lake-resident Arctic charr, Salvelinus alpinus from northwest Europe on multiple spatial scales, 2367 individuals from 43 lakes located in three geographical regions (Iceland, the British Isles and Scandinavia) were genotyped at six microsatellite loci. On a large scale, data provided little evidence to support clustering of populations according to geographical region. Hierarchical analysis of molecular variance indicated that, although statistically significant, only 2.17% of the variance in allelic frequencies was partitioned at the among-region level. Within regions, high levels of genetic differentiation were typically found between lakes regardless of the geographical distance separating them. These results are consistent with the hypothesis of rapid postglacial recolonization of all of northwest Europe from a single charr lineage, with subsequent restriction of gene flow. On a smaller scale, there was evidence for close genetic relationships among lakes from within common drainage basins in Scotland. Thus, interlake genetic structure reflects localized patterns of recent (or contemporary) gene flow superimposed onto a larger scale structure that is largely a result of historical processes. There was also evidence for widespread genetic structuring at the within-lake level, with sympatric populations detected in 10 lakes, and multilocus heterozygote deficits found in 23 lakes. This evidence of the Wahlund effect was found in all lakes known to contain discrete phenotypic morphs, as well as many others, suggesting that morphs may often represent separate breeding populations, and also that the phenomenon of polymorphism in this species may be more widespread than is currently realized.

A genetic linkage map for Arctic char (Salvelinus alpinus): evidence for higher recombination rates and segregation distortion in hybrid versus pure strain mapping parents.

We constructed a genetic linkage map for Arctic char (Salvelinus alpinus) using two backcrosses between genetically divergent strains. Forty-six linkage groups (expected = 39-41) and 19 homeologous affinities (expected = 25) were identified using 184 microsatellites, 129 amplified fragment length polymorphisms (AFLPs), 13 type I gene markers, and one phenotypic marker, SEX. Twenty-six markers remain unlinked. Female map distance (9.92 Morgans) was substantially higher than male map distance (3.90 Morgans) based on the most complete parental information (i.e., the F1 hybrids). Female recombination rates were often significantly higher than those of males across all pairwise comparisons within homologous chromosomal segments (average female to male ratios within families was 1.69:1). The female hybrid parent had significantly higher recombination rates than the pure strain female parent. Segregation distortion was detected in four linkage groups (4, 8, 13, 20) for both families. In family 3, only the largest fish were sampled for genotyping, suggesting that segregation distortion may represent regions possessing influences on growth. In family 2, almost all cases showing segregation distortion involved markers in the female hybrid parent.

Comparison of GNRH3 genes across salmonid genera.

Genomic sequences of gonadotropin-releasing hormone genes were amplified and examined for sequence divergence among members of three different genera of the subfamily Salmoninae: rainbow trout (Oncorhynchus mykiss), Atlantic salmon (Salmo salar), and Arctic charr (Salvelinus alpinus). Sequences of GNRH3A and GNRH3B (formerly known as sGnRH1 and sGnRH2) were 97-99% similar in coding regions and 94-98% similar in non-coding regions among genera, but comparisons within species between GNRH3A and GNRH3B were only 90-92% similar in coding regions and 83-89% similar in non-coding regions. Polymorphisms in the parents of mapping families for each species allowed for linkage mapping of the GNRH3B gene in all three species and the GNRH3A gene in rainbow trout. GNRH3B maps to linkage group 6 in rainbow trout, linkage group 16 in Atlantic salmon and linkage group 25 in Arctic charr. GNRH3A mapped to linkage group 30 in rainbow trout.

Quantitative trait loci for spawning date and body weight in rainbow trout: testing for conserved effects across ancestrally duplicated chromosomes.

We incorporated 69 microsatellite loci into an existing data set of 132 markers to test for quantitative trait loci (QTLs) affecting spawning date and body weight in a backcross between two outbred strains of rainbow trout (Oncorhynchus mykiss). Twenty-six linkage groups were identified and synteny of duplicated microsatellite markers was used to confirm 13 homeologous chromosome pairs. Gene-centromere data were used to localize the centromeres for 13 linkage groups whose orientations were previously unknown. We applied a combination of interval mapping and single marker analysis to the segregating maternal and paternal alleles at 201 microsatellite loci. Four spawning date QTLs with suggestive evidence for an additional two QTLs were detected in female trout spawning at 3 and 4 years of age. Similarly we detected three QTLs for body weight in females at 2 years of age plus four suggestive QTLs for this trait. We found marginal evidence that three pairs of ancestral homeologues contained detectable QTLs for the same trait. In one of the three pairs of homeologues, the duplicated QTL regions mapped to the same relative chromosomal location, while the exact localization of the QTL position in one of the other pairs was difficult to infer since it was based on data from a male-derived map. The existing data were unable to refute a hypothesis that duplicated functional genes will be maintained within the telomeric regions of salmonids due to preferential male-mediated crossing over in this region. Two of the four spawning date QTLs were detected on linkage groups with unknown homeologous relationships. QTLs with possible pleiotropic effects on both spawning date and body size were localized to two linkage groups.

Effects of genetic sex and genomic background on epistasis in rainbow trout (Oncorhynchus mykiss).

Epistasis among quantitative trait loci (QTL) for survival (upper thermal tolerance, UTT) and morphological (fork length, FL and condition factor, K) traits was detected in purestrain and interstrain rainbow trout (Oncorhynchus mykiss) families. One sex-linked (OmyFGT19TUF) and three autosomal (Omy325UoG, Ssa14DU and Ssa20.19NUIG; linkage groups B, D and S, respectively) microsatellite loci linked to UTT QTL in this species were used. Within half sib families, significant effects of full sib family on epistasis involving Omy325UoG and OmyFGT19TUF were detected at a rate significantly higher than expected for UTT (p < 0.001*) and FL (p < 0.01*), using results significant at comparisonwise significance thresholds derived from permutational analysis. Measured across half sib families, the phenotype of female genotypic classes was more divergent from the family trait mean than that of males where epistasis involved the sex-linked locus OmyFGT19TUF (p = 0.0176*), and also for means over all families (p = 0.0355*). Female genotypic classes were also more divergent (p = 0.0011 **) from the full sib trait mean where three-way interaction between OmyFGT19TUF, one of the autosomal loci and full sib family was significant, and marginally more divergent for trait means of genotypic classes across all full sib families (p = 0.0856+). There was no evidence that these effects were more pronounced in hybrid F1 families than purestrains.