A genetic linkage map for the salmon louse (Lepeophtheirus salmonis): evidence for high male:female and inter-familial recombination rate differences

A salmon louse (Lepeophtheirus salmonis salmonis) genetic linkage map was constructed to serve as a genomic resource for future investigations into the biology of this important marine parasitic copepod species, and to provide insights into the inheritance patterns of genetic markers in this species. SNP genotyping of 8 families confirmed the presence of 15 linkage groups based upon the assignment of 93,773 markers. Progeny sample size weight adjusted map sizes in males (with the exception of SL12 and SL15) ranged in size from 96.50 cM (SL11) to 134.61 cM (SL06), and total combined map steps or bins ranged from 143 (SL09) to 203 (SL13). The SL12 male map was the smallest linkage group with a weight-averaged size of 3.05 cM with 6 recombination bins. Male:female specific recombination rate differences are 10.49:1 and represent one of the largest reported sex-specific differences for any animal species. Recombination ratio differences (M:F) ranged from 1.0 (SL12) to 29:1 (SL15). The number of markers exhibiting normal Mendelian segregation within the sex linkage group SL15 was extremely low (N = 80) in comparison to other linkage groups genotyped [range: 1459 (SL12)-10206 markers (SL05)]. Re-evaluation of Mendelian inheritance patterns of markers unassigned to any mapping parent according to hemizygous segregation patterns (models presented) identified matches for many of these markers to hemizygous patterns. The greatest proportion of these markers assigned to SL15 (N increased to 574). Inclusion of the hemizygous markers revised SL15 sex-specific recombination rate differences to 28:1. Recombination hot- and coldspots were identified across all linkage groups with all linkage groups possessing multiple peaks. Nine of 13 linkage groups evaluated possessed adjacent domains with hot-coldspot transitional zones. The most common pattern was for one end of the linkage to show elevated recombination in addition to internal regions. For SL01 and SL06, however, a terminal region with high recombination was not evident while a central domain possessing extremely high-recombination levels was present. High levels of recombination were weakly coupled to higher levels of SNP variation within domains, but this association was very strong for the central domains of SL01 and SL06. From the pooled paternal half-sib lots (several virgin females placed with 1 male), only 1 or two surviving family lots were obtained. Surviving families possessed parents where both the male and female possessed either inherently low or high recombination rates. This study provides insight into the organization of the sea louse genome, and describes large differences in recombination rate that exist among individuals of the same sex, and between the sexes. These differences in recombination rate may be coupled to the capabilities of this species to adapt to environmental and pharmaceutical treatments, given that family survivorship appears to be enhanced when parents have similar recombination levels.

The Arctic charr (Salvelinus alpinus) genome and transcriptome assembly

Arctic charr have a circumpolar distribution, persevere under extreme environmental conditions, and reach ages unknown to most other salmonids. The Salvelinus genus is primarily composed of species with genomes that are structured more like the ancestral salmonid genome than most Oncorhynchus and Salmo species of sister genera. It is thought that this aspect of the genome may be important for local adaptation (due to increased recombination) and anadromy (the migration of fish from saltwater to freshwater). In this study, we describe the generation of a new genetic map, the sequencing and assembly of the Arctic charr genome (GenBank accession: GCF_002910315.2) using the newly created genetic map and a previous genetic map, and present several analyses of the Arctic charr genes and genome assembly. The newly generated genetic map consists of 8,574 unique genetic markers and is similar to previous genetic maps with the exception of three major structural differences. The N50, identified BUSCOs, repetitive DNA content, and total size of the Arctic charr assembled genome are all comparable to other assembled salmonid genomes. An analysis to identify orthologous genes revealed that a large number of orthologs could be identified between salmonids and many appear to have highly conserved gene expression profiles between species. Comparing orthologous gene expression profiles may give us a better insight into which genes are more likely to influence species specific phenotypes.

A 200K SNP chip reveals a novel Pacific salmon louse genotype linked to differential efficacy of emamectin benzoate

Antiparasitic drugs such as emamectin benzoate (EMB) are relied upon to reduce the parasite load, particularly of the sea louse Lepeophtheirus salmonis, on farmed salmon. The decline in EMB treatment efficacy for this purpose is an important issue for salmon producers around the world, and particularly for those in the Atlantic Ocean where widespread EMB tolerance in sea lice is recognized as a significant problem. Salmon farms in the Northeast Pacific Ocean have not historically experienced the same issues with treatment efficacy, possibly due to the relatively large population of endemic salmonid hosts that serve to both redistribute surviving lice and dilute populations potentially under selection by introducing naïve lice to farms. Frequent migration of lice among farmed and wild hosts should limit the effect of farm-specific selection pressures on changes to the overall allele frequencies of sea lice in the Pacific Ocean. A previous study using microsatellites examined L. salmonis oncorhynchi from 10 Pacific locations from wild and farmed hosts and found no population structure. Recently however, a farm population of sea lice was detected where EMB bioassay exposure tolerance was abnormally elevated. In response, we have developed a Pacific louse draft genome that complements the previously-released Atlantic louse sequence. These genomes were combined with whole-genome re-sequencing data to design a highly sensitive 201,279 marker SNP array applicable for both subspecies (90,827 validated Pacific loci; 153,569 validated Atlantic loci). Notably, kmer spectrum analysis of the re-sequenced samples indicated that Pacific lice exhibit a large within-individual heterozygosity rate (average of 1 in every 72 bases) that is markedly higher than that of Atlantic individuals (1 in every 173 bases). The SNP chip was used to produce a high-density map for Atlantic sea louse linkage group 5 that was previously shown to be associated with EMB tolerance in Atlantic lice. Additionally, 478 Pacific louse samples from farmed and wild hosts obtained between 2005 and 2014 were also genotyped on the array. Clustering analysis allowed us to detect the apparent emergence of an otherwise rare genotype at a high frequency among the lice collected from two farms in 2013 that had reported elevated EMB tolerance. This genotype was not observed in louse samples collected from the same farm in 2010, nor in any lice sampled from other locations prior to 2013. However, this genotype was detected at low frequencies in louse samples from farms in two locations reporting elevated EMB tolerance in 2014. These results suggest that a rare genotype present in Pacific lice may be locally expanded in farms after EMB treatment. Supporting this hypothesis, 437 SNPs associated with this genotype were found to be in a region of linkage group 5 that overlaps the region associated with EMB resistance in Atlantic lice. Finally, five of the top diagnostic SNPs within this region were used to screen lice that had been subjected to an EMB survival assay, revealing a significant association between these SNPs and EMB treatment outcome. To our knowledge this work is the first report to identify a genetic link to altered EMB efficacy in L. salmonis in the Pacific Ocean.

The channel catfish genome sequence provides insights into the evolution of scale formation in teleosts

Catfish represent 12% of teleost or 6.3% of all vertebrate species, and are of enormous economic value. Here we report a high-quality reference genome sequence of channel catfish (Ictalurus punctatus), the major aquaculture species in the US. The reference genome sequence was validated by genetic mapping of 54,000 SNPs, and annotated with 26,661 predicted protein-coding genes. Through comparative analysis of genomes and transcriptomes of scaled and scaleless fish and scale regeneration experiments, we address the genomic basis for the most striking physical characteristic of catfish, the evolutionary loss of scales and provide evidence that lack of secretory calcium-binding phosphoproteins accounts for the evolutionary loss of scales in catfish. The channel catfish reference genome sequence, along with two additional genome sequences and transcriptomes of scaled catfishes, provide crucial resources for evolutionary and biological studies. This work also demonstrates the power of comparative subtraction of candidate genes for traits of structural significance.

Transcriptome profiling in fast versus slow-growing rainbow trout across seasonal gradients

Background

Circannual rhythms in vertebrates can influence a wide variety of physiological processes. Some notable examples include annual reproductive cycles and for poikilotherms, seasonal changes modulating growth. Increasing water temperature elevates growth rates in fishes, but increases in photoperiod regime can have similar influences even at constant temperature. Therefore, in order to understand the dynamics of growth in fish it is important to consider the background influence of photoperiod regime on gene expression differences. This study examined the influence of a declining photoperiod regime (winter solstice) compared to an increasing photoperiod regime (spring equinox) on white muscle transcriptome profiles in fast and slow-growing rainbow trout from a commercial aquaculture strain.

Results

Slow-growing fish could be characterized as possessing transcriptome profiles that conform in many respects to an endurance training regime in humans. They have elevated mitochondrial and cytosolic creatine kinase expression levels and appear to suppress mTOR-signaling as evidenced by elevated TSC2 expression, and they also have elevated p53 levels. Large fish display a physiological repertoire that may be consistent with strength/resistance physiology having elevated cytoskeletal gene component expression and glycogen metabolism cycling along with higher PI3K levels. In many respects small vs. large fish match eccentric vs. concentric muscle expression patterns, respectively. Lipid metabolic genes are also more elevated in larger fish, the most notable being the G0S2 switch gene. M and Z-line sarcomere remodelling appears to be more prevalent in large fish. Twenty-three out of 26 gene families with previously reported significant SNP-based growth differences were detected as having significant expression differences.

Conclusions

Larger fish display a broader array of genes showing higher expression, and their profiles are more similar to those observed in December lot fish (i.e., an accelerated growth period). Conversely, small fish display gene profiles more similar to seasonal growth decline phases (i.e., September lot fish). Overall, seasonal timing was coupled to greater differences in gene expression compared to differences associated with fish size.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2363-5) contains supplementary material, which is available to authorized users.

Co-localization of growth QTL with differentially expressed candidate genes in rainbow trout

We tested whether genes differentially expressed between large and small rainbow trout co-localized with familial QTL regions for body size. Eleven chromosomes, known from previous work to house QTL for weight and length in rainbow trout, were examined for QTL in half-sibling families produced in September (1 XY male and 1 XX neomale) and December (1 XY male). In previous studies, we identified 108 candidate genes for growth expressed in the liver and white muscle in a subset of the fish used in this study. These gene sequences were BLASTN aligned against the rainbow trout and stickleback genomes to determine their location (rainbow trout) and inferred location based on synteny with the stickleback genome. Across the progeny of all three males used in the study, 63.9% of the genes with differential expression appear to co-localize with the QTL regions on 6 of the 11 chromosomes tested in these males. Genes that co-localized with QTL in the mixed-sex offspring of the two XY males primarily showed up-regulation in the muscle of large fish and were related to muscle growth, metabolism, and the stress response.

Comparison of growth-related traits and gene expression profiles between the offspring of neomale (XX) and normal male (XY) rainbow trout

All-female lines of fish are created by crossing sex reversed (XX genotype) males with normal females. All-female lines avoid the deleterious phenotypic effects that are typical of precocious maturation in males. To determine whether all-female and mixed sex populations of rainbow trout (Oncorhynchus mykiss) differ in performance, we compared the growth and gene expression profiles in progeny groups produced by crossing a XX male and a XY male to the same five females. Body weight and length were measured in the resulting all-female (XX) and mixed sex (XX/XY) offspring groups. Microarray experiments with liver and white muscle were used to determine if the gene expression profiles of large and small XX offspring differ from those in large and small XX/XY offspring. We detected no significant differences in body length and weight between offspring groups but XX offspring were significantly less variable in the value of these traits. A large number of upregulated genes were shared between the large XX and large XX/XY offspring; the small XX and small XX/XY offspring also shared similar expression profiles. No GO category differences were seen in the liver or between the large XX and large XX/XY offspring in the muscle. The greatest differences between the small XX and small XX/XY offspring were in the genes assigned to the "small molecule metabolic process" and "cellular metabolic process" GO level 3 categories. Similarly, genes within these categories as well as the category "macromolecule metabolic process" were more highly expressed in small compared to large XX fish.

Construction of a high-density, high-resolution genetic map and its integration with BAC-based physical map in channel catfish

Construction of genetic linkage map is essential for genetic and genomic studies. Recent advances in sequencing and genotyping technologies made it possible to generate high-density and high-resolution genetic linkage maps, especially for the organisms lacking extensive genomic resources. In the present work, we constructed a high-density and high-resolution genetic map for channel catfish with three large resource families genotyped using the catfish 250K single-nucleotide polymorphism (SNP) array. A total of 54,342 SNPs were placed on the linkage map, which to our knowledge had the highest marker density among aquaculture species. The estimated genetic size was 3,505.4 cM with a resolution of 0.22 cM for sex-averaged genetic map. The sex-specific linkage maps spanned a total of 4,495.1 cM in females and 2,593.7 cM in males, presenting a ratio of 1.7 : 1 between female and male in recombination fraction. After integration with the previously established physical map, over 87% of physical map contigs were anchored to the linkage groups that covered a physical length of 867 Mb, accounting for ∼90% of the catfish genome. The integrated map provides a valuable tool for validating and improving the catfish whole-genome assembly and facilitates fine-scale QTL mapping and positional cloning of genes responsible for economically important traits.

Phenotypic and QTL allelic associations among embryonic developmental rate, body size, and precocious maturation in male rainbow trout

We examined associations among embryonic developmental rate (EDR) as measured by hatching time, juvenile body weight (BW) and propensity for precocial sexual maturation (PM) at two years in two sets of diallel crosses of rainbow trout produced in two spawning seasons (September and December) at both the phenotypic and genotypic levels. Dams and sires had highly significant effects on the body weight of their male juvenile progeny on three measurement dates where parental effects remained consistent through time. Dams spawning earlier in the season produced a greater number of mature male progeny (56.7%) than did later spawning females (25.6%). The families from the December lot showed the expected associations among traits in that earlier hatching fish were significantly heavier on all three measurement dates than later hatching fish and were more likely to mature earlier when families were combined. Moreover, earlier maturing fish were significantly heavier on the third measurement date than those that did not mature. In the September lot, mature fish were significantly heavier as juveniles on all three measurement dates than immature fish as predicted but no significant associations were detected between EDR and BW or between PM and EDR. Significant QTL were detected for all three traits but the linkage group location varied depending on the trait and half-sib group analyzed (across dams and sires in each lot). A strong QTL for EDR with genome-wide effects was detected on linkage group RT-8 in all four half-sib analyses. None of the four linkage groups analyzed had QTL for all three traits. However, the phenotypic association between EDR and BW observed in the December lot was supported by the co-localization of QTL to linkage group RT-8 and a positive coupling of allelic effects. RT-8 marker alleles significantly associated with faster EDR were also associated with larger BW and this was observed in numerous families on all three measurement dates. Linkage group RT-24 had weaker QTL for all three traits in the September lot but these were not detected in the same half-sib group simultaneously. At the allelic level, marker alleles for faster EDR were also associated with BW but only at the third measurement date and the progeny of one male. Similarly, RT-30 had weaker QTL for EDR and PM in the December paternal half-sib analysis but no associations were evident at the allelic level. The detection of associations between life history traits and growth at both the phenotypic and genotypic levels has significant implications to aquaculture breeding programs where selection for a desirable trait may lead to unwanted alterations of other traits. Furthermore, the differences between spawning season lots emphasize the complex interaction between environment and genotype on economically important traits and the resulting challenges for aquaculture.

Hidden genetic variation evolves with ecological specialization: the genetic basis of phenotypic plasticity in Arctic charr ecomorphs

The genetic variance that determines phenotypic variation can change across environments through developmental plasticity and in turn play a strong role in evolution. Induced changes in genotype-phenotype relationships should strongly influence adaptation by exposing different sets of heritable variation to selection under some conditions, while also hiding variation. Therefore, the heritable variation exposed or hidden from selection is likely to differ among habitats. We used ecomorphs from two divergent populations of Arctic charr (Salvelinus alpinus) to test the prediction that genotype-phenotype relationships would change in relation to environment. If present over several generations this should lead to divergence in genotype-phenotype relationships under common conditions, and to changes in the amount and type of hidden genetic variance that can evolve. We performed a common garden experiment whereby two ecomorphs from each of two Icelandic lakes were reared under conditions that mimicked benthic and limnetic prey to induce responses in craniofacial traits. Using microsatellite based genetic maps, we subsequently detected QTL related to these craniofacial traits. We found substantial changes in the number and type of QTL between diet treatments and evidence that novel diet treatments can in some cases provide a higher number of QTL. These findings suggest that selection on phenotypic variation, which is both genetically and environmentally determined, has shaped the genetic architecture of adaptive divergence in Arctic charr. However, while adaptive changes are occurring in the genome there also appears to be an accumulation of hidden genetic variation for loci not expressed in the contemporary environment.

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.

Differential gene expression in small and large rainbow trout derived from two seasonal spawning groups

Background

Growth in fishes is regulated via many environmental and physiological factors and is shaped by the genetic background of each individual. Previous microarray studies of salmonid growth have examined fish experiencing either muscle wastage or accelerated growth patterns following refeeding, or the influence of growth hormone and transgenesis. This study determines the gene expression profiles of genetically unmanipulated large and small fish from a domesticated salmonid strain reared on a typical feeding regime. Gene expression profiles of white muscle and liver from rainbow trout (Oncorhynchus mykiss) from two seasonal spawning groups (September and December lots) within a single strain were examined when the fish were 15 months of age to assess the influence of season (late fall vs. onset of spring) and body size (large vs. small).

Results

Although IGFBP1 gene expression was up-regulated in the livers of small fish in both seasonal lots, few expression differences were detected in the liver overall. Faster growing Dec. fish showed a greater number of differences in white muscle expression compared to Sept. fish. Significant differences in the GO Generic Level 3 categories ‘response to external stimulus’, ‘establishment of localization’, and ‘response to stress’ were detected in white muscle tissue between large and small fish. Larger fish showed up-regulation of cytoskeletal component genes while many genes related to myofibril components of muscle tissue were up-regulated in small fish. Most of the genes up-regulated in large fish within the ‘response to stress’ category are involved in immunity while in small fish most of these gene functions are related to apoptosis.

Conclusions

A higher proportion of genes in white muscle compared to liver showed similar patterns of up- or down-regulation within the same size class across seasons supporting their utility as biomarkers for growth in rainbow trout. Differences between large and small Sept. fish in the ‘response to stress’ and ‘response to external stimulus’ categories for white muscle tissue, suggests that smaller fish have a greater inability to handle stress compared to the large fish. Sampling season had a significant impact on the expression of genes related to the growth process in rainbow trout.

An integrated transcriptomic and comparative genomic analysis of differential gene expression in Arctic charr (Salvelinus alpinus) following seawater exposure

High-throughput RNA sequencing was employed to compare expression profiles in two Arctic charr (Salvelinus alpinus) families post seawater exposure to identify genes and biological processes involved in hypo-osmoregulation and regulation of salinity tolerance. To further understand the genetic architecture of hypo-osmoregulation, the genomic organization of differentially expressed (DE) genes was also analysed. Using a de novo gill transcriptome assembly we found over 2300 contigs to be DE. Major transporters from the seawater mitochondrion-rich cell (MRC) complex were up-regulated in seawater. Expression ratios for 257 differentially expressed contigs were highly correlated between families, suggesting they are strictly regulated. Based on expression profiles and known molecular pathways we inferred that seawater exposure induced changes in methylation states and elevated peroxynitrite formation in gill. We hypothesized that concomitance between DE immune genes and the transition to a hypo-osmoregulatory state could be related to Cl- sequestration by antimicrobial defence mechanisms. Gene Ontology analysis revealed that cell division genes were up-regulated, which could reflect the proliferation of ATP1α1b-type seawater MRCs. Comparative genomics analyses suggest that hypo-osmoregulation is influenced by the relative proximities among a contingent of genes on Arctic charr linkage groups AC-4 and AC-12 that exhibit homologous affinities with a region on stickleback chromosome Ga-I. This supports the hypothesis that relative gene location along a chromosome is a property of the genetic architecture of hypo-osmoregulation. Evidence of non-random structure between hypo-osmoregulation candidate genes was found on AC-1/11 and AC-28, suggesting that interchromosomal rearrangements played a role in the evolution of hypo-osmoregulation in Arctic charr.

Transcriptomics of salinity tolerance capacity in Arctic charr (Salvelinus alpinus): a comparison of gene expression profiles between divergent QTL genotypes

Osmoregulatory capabilities have played an important role in the evolution, dispersal, and diversification of vertebrates. To better understand the genetic architecture of hypo-osmoregulation in fishes and to determine which genes and biological processes affect intraspecific variation in salinity tolerance, we used mRNA sequence libraries from Arctic charr gill tissue to compare gene expression profiles in fish exhibiting divergent salinity tolerance quantitative trait locus (QTL) genotypes. We compared differentially expressed genes with QTL positions to gain insight about the nature of the underlying polymorphisms and examined gene expression within the context of genome organization to gain insight about the evolution of hypo-osmoregulation in fishes. mRNA sequencing of 18 gill tissue libraries produced 417 million reads, and the final reduced de novo transcriptome assembly consisted of 92,543 contigs. Families contained a similar number of differentially expressed contigs between high and low salinity tolerance capacity groups, and log2 expression ratios ranged from 10.4 to -8.6. We found that intraspecific variation in salinity tolerance capacity correlated with differential expression of immune response genes. Some differentially expressed genes formed clusters along linkage groups. Most clusters comprised gene pairs, though clusters of three, four, and eight genes were also observed. We postulated that conserved synteny of gene clusters on multiple ancestral and teleost chromosomes may have been preserved via purifying selection. Colocalization of QTL with differentially expressed genes suggests that polymorphisms in cis-regulatory elements are part of a majority of QTL.

Influence of developmental stage and genotype on liver mRNA levels among wild, domesticated, and hybrid rainbow trout (Oncorhynchus mykiss)

BACKGROUND

Release of domesticated strains of fish into nature may pose a threat to wild populations with respect to their evolved genetic structure and fitness. Understanding alterations that have occurred in both physiology and genetics as a consequence of domestication can assist in evaluating the risks posed by introgression of domesticated genomes into wild genetic backgrounds, however the molecular causes of these consequences are currently poorly defined. The present study has examined levels of mRNA in fast-growing pure domesticated (D), slow-growing age-matched pure wild (Wa), slow-growing size-matched pure wild (Ws), and first generation hybrid cross (W/D) rainbow trout (Oncorhynchus mykiss) to investigate the influence of genotype (domesticated vs. wild, and their interactions in hybrids) and developmental stage (age- or size-matched animals) on genetic responses (i.e. dominant vs. recessive) and specific physiological pathways.

RESULTS

Highly significant differences in mRNA levels were found between domesticated and wild-type rainbow trout genotypes (321 mRNAs), with many mRNAs in the wild-domesticated hybrid progeny showing intermediate levels. Differences were also found between age-matched and size-matched wild-type trout groups (64 mRNAs), with unique mRNA differences for each of the wild-type groups when compared to domesticated trout (Wa: 114 mRNAs, Ws: 88 mRNAs), illustrating an influence of fish developmental stage affecting findings when used as comparator groups to other genotypes. Analysis of differentially expressed mRNAs (found for both wild-type trout to domesticated comparisons) among the genotypes indicates that 34.8% are regulated consistent with an additive genetic model, whereas 39.1% and 26.1% show a recessive or dominant mode of regulation, respectively. These molecular data are largely consistent with phenotypic data (growth and behavioural assessments) assessed in domesticated and wild trout strains.

CONCLUSIONS

The present molecular data are concordant with domestication having clearly altered rainbow trout genomes and consequent phenotype from that of native wild populations. Although mainly additive responses were noted in hybrid progeny, the prevalence of dominant and non-additive responses reveals that introgression of domesticated and wild genotypes alters the type of genetic control of mRNA levels from that of wild-type, which may lead to disruption of gene regulation systems important for developing phenotypes for optimal fitness in nature. A clear influence of both fish age and size (developmental stage) on mRNA levels was also noted in this study, which highlights the importance of examining multiple control samples to provide a comprehensive understanding of changes observed between strains possessing differences in growth rate.

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.

Genomic arrangement of salinity tolerance QTLs in salmonids: a comparative analysis of Atlantic salmon (Salmo salar) with Arctic charr (Salvelinus alpinus) and rainbow trout (Oncorhynchus mykiss)

Background

Quantitative trait locus (QTL) studies show that variation in salinity tolerance in Arctic charr and rainbow trout has a genetic basis, even though both these species have low to moderate salinity tolerance capacities. QTL were observed to localize to homologous linkage group segments within putative chromosomal regions possessing multiple candidate genes. We compared salinity tolerance QTL in rainbow trout and Arctic charr to those detected in a higher salinity tolerant species, Atlantic salmon. The highly derived karyotype of Atlantic salmon allows for the assessment of whether disparity in salinity tolerance in salmonids is associated with differences in genetic architecture. To facilitate these comparisons, we examined the genomic synteny patterns of key candidate genes in the other model teleost fishes that have experienced three whole-genome duplication (3R) events which preceded a fourth (4R) whole genome duplication event common to all salmonid species.

Results

Nine linkage groups contained chromosome-wide significant QTL (AS-2, -4p, -4q, -5, -9, -12p, -12q, -14q -17q, -22, and −23), while a single genome-wide significant QTL was located on AS-4q. Salmonid genomes shared the greatest marker homology with the genome of three-spined stickleback. All linkage group arms in Atlantic salmon were syntenic with at least one stickleback chromosome, while 18 arms had multiple affinities. Arm fusions in Atlantic salmon were often between multiple regions bearing salinity tolerance QTL. Nine linkage groups in Arctic charr and six linkage group arms in rainbow trout currently have no synteny alignments with stickleback chromosomes, while eight rainbow trout linkage group arms were syntenic with multiple stickleback chromosomes. Rearrangements in the stickleback lineage involving fusions of ancestral arm segments could account for the 21 chromosome pairs observed in the stickleback karyotype.

Conclusions

Salinity tolerance in salmonids from three genera is to some extent controlled by the same loci. Synteny between QTL in salmonids and candidate genes in stickleback suggests genetic variation at candidate gene loci could affect salinity tolerance in all three salmonids investigated. Candidate genes often occur in pairs on chromosomes, and synteny patterns indicate these pairs are generally conserved in 2R, 3R, and 4R genomes. Synteny maps also suggest that the Atlantic salmon genome contains three larger syntenic combinations of candidate genes that are not evident in any of the other 2R, 3R, or 4R genomes examined. These larger synteny tracts appear to have resulted from ancestral arm fusions that occurred in the Atlantic salmon ancestor. We hypothesize that the superior hypo-osmoregulatory efficiency that is characteristic of Atlantic salmon may be related to these clusters.

Sex chromosome polymorphisms in Arctic charr and their evolutionary origins

Current data on the Y-specific sex-determining region of salmonid fishes from genera Salvelinus, Salmo, and Oncorhynchus indicate variable polymorphisms in the homologous chromosomal locations of the sex-specific determining region. In the majority of the Atlantic lineage Arctic charr, including populations from the Fraser River, in Labrador Canada, as well as Swedish and Norwegian strains, the sex-determining locus maps to linkage group AC-4. Previously, sex-linked polymorphisms (i.e., variation in the associated sex-linked markers on AC-4) have been described in Arctic charr. Here, we report further evidence for intraspecific sex linkage group polymorphisms in Arctic charr (i.e., the detection of the SEX locus on either the AC-1 or AC-21 linkage group) and a possible conservation of a sex linkage arrangement in Icelandic Arctic charr and Atlantic salmon, involving sex-linked markers on the AC-1/21 homeologs and the European AS-1/6 homeologous linkage groups in Atlantic salmon. The evolutionary origins for the multiple sex-determining regions within the salmonid family are discussed. We also relate the variable sex-determining regions in salmonids to their ancestral proto-teleost karyotypic origins and compare these findings with what has been observed in other teleost species in general.

The genetic basis of salinity tolerance traits in Arctic charr (Salvelinus alpinus)

BACKGROUND

The capacity to maintain internal ion homeostasis amidst changing conditions is particularly important for teleost fishes whose reproductive cycle is dependent upon movement from freshwater to seawater. Although the physiology of seawater osmoregulation in mitochondria-rich cells of fish gill epithelium is well understood, less is known about the underlying causes of inter- and intraspecific variation in salinity tolerance. We used a genome-scan approach in Arctic charr (Salvelinus alpinus) to map quantitative trait loci (QTL) correlated with variation in four salinity tolerance performance traits and six body size traits. Comparative genomics approaches allowed us to infer whether allelic variation at candidate gene loci (e.g., ATP1α1b, NKCC1, CFTR, and cldn10e) could have underlain observed variation.

RESULTS

Combined parental analyses yielded genome-wide significant QTL on linkage groups 8, 14 and 20 for salinity tolerance performance traits, and on 1, 19, 20 and 28 for body size traits. Several QTL exhibited chromosome-wide significance. Among the salinity tolerance performance QTL, trait co-localizations occurred on chromosomes 1, 4, 7, 18 and 20, while the greatest experimental variation was explained by QTL on chromosomes 20 (19.9%), 19 (14.2%), 4 (14.1%) and 12 (13.1%). Several QTL localized to linkage groups exhibiting homeologous affinities, and multiple QTL mapped to regions homologous with the positions of candidate gene loci in other teleosts. There was no gene × environment interaction among body size QTL and ambient salinity.

CONCLUSIONS

Variation in salinity tolerance capacity can be mapped to a subset of Arctic charr genomic regions that significantly influence performance in a seawater environment. The detection of QTL on linkage group 12 was consistent with the hypothesis that variation in salinity tolerance may be affected by allelic variation at the ATP1α1b locus. IGF2 may also affect salinity tolerance capacity as suggested by a genome-wide QTL on linkage group 19. The detection of salinity tolerance QTL in homeologous regions suggests that candidate loci duplicated from the salmonid-specific whole-genome duplication may have retained their function on both sets of homeologous chromosomes. Homologous affinities suggest that loci affecting salinity tolerance in Arctic charr may coincide with QTL for smoltification and salinity tolerance traits in rainbow trout. The effects of body size QTL appear to be independent of changes in ambient salinity.

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.

Genome evolution in the fish family salmonidae: generation of a brook charr genetic map and comparisons among charrs (Arctic charr and brook charr) with rainbow trout

BACKGROUND

Salmonids are regarded as 4R derivative species, having experienced 4 whole genome duplication events in their ancestry. Many duplicated chromosome regions still share extensive homology with one another which is maintained primarily through male-based homeologous chromosome pairings during meiosis. The formation of quadrivalents during meiosis leads to pseudolinkage. This phenomenon is more prevalent within 5 of the 12 ancestral teleost linkage groups in salmonids.

RESULTS

We constructed a genetic linkage map for brook charr and used this in combination with the genetic map from Arctic charr, to make comparisons with the genetic map of rainbow trout. Although not all chromosome arms are currently mapped, some homologous chromosome rearrangements were evident between Arctic charr and brook charr. Notably, 10 chromosome arms in brook charr representing 5 metacentric chromosomes in Arctic charr have undergone rearrangements. Three metacentrics have one arm translocated and fused with another chromosome arm in brook charr to a make a new metacentrics while two metacentrics are represented by 4 acrocentric pairs in brook charr. In two cases (i.e., BC-4 and BC-16), an apparent polymorphism was observed with the identification of both a putative metacentric structure (similar to metacentric AC-4 = BC-4 and a joining of acrocentric AC-16 + one arm of AC-28 = BC-16), as well as two separate acrocentric linkage groups evident in the mapping parents. Forty-six of the expected 50 karyotypic arms could be inter-generically assigned. SEX in brook charr (BC-4) was localized to the same homologous linkage group region as in Arctic charr (AC-4). The homeologous affinities detected in the two charr species facilitated the identification of 20 (expected number = 25) shared syntenic regions with rainbow trout, although it is likely that some of these regions were partial or overlapping arm regions.

CONCLUSIONS

Inter-generic comparisons among 2 species of charr (genus Salvelinus) and a trout (genus Oncorhynchus) have identified that linkage group arm arrangements are largely retained among these species. Previous studies have revealed that up to 7 regions of high duplicate marker retention occur between Salmo species (i.e., Atlantic salmon and brown trout) and rainbow trout, with 5 of these regions exhibiting higher levels of pseudolinkage. Pseudolinkage was detected in the charr species (i.e., BC-1/21, AC-12/27, AC-6/23, = RT-2p/29q, RT-12p/16p, and RT-27p/31p, respectively) consistent with three of the five 'salmonid-specific' pseudolinkage regions. Chromosome arms with the highest number of duplicated markers in rainbow trout are the linkage group arms with the highest retention of duplicated markers in both charr species.

Genetic architecture of body weight, condition factor and age of sexual maturation in Icelandic Arctic charr (Salvelinus alpinus)

The high commercial value from the aquaculture of salmonid fishes has prompted many studies into the genetic architecture of complex traits and the need to identify genomic regions that have repeatable associations with trait variation both within and among species. We searched for quantitative trait loci (QTL) for body weight (BW), condition factor (CF) and age of sexual maturation (MAT) in families of Arctic charr (Salvelinus alpinus) from an Icelandic breeding program. QTL with genome-wide significance were detected for each trait on multiple Arctic charr (AC) linkage groups (BW: AC-4, AC-20; CF: AC-7, AC-20, AC-23, AC-36; MAT: AC-13/34, AC-39). In addition to the genome-wide significant QTL for both BW and CF on AC-20, linkage groups AC-4, AC-7, AC-8, and AC-16 contain QTL for both BW and CF with chromosome-wide significance. These regions had effects (albeit weaker) on MAT with the exception of the region on AC-8. Comparisons with a North American cultured strain of Arctic charr, as well as North American populations of Atlantic salmon (Salmo salar), and rainbow trout (Oncorhynchus mykiss), reveal some conservation in QTL location and structure, particularly with respect to the joint associations of QTL influencing BW and CF. The detection of some differences in genetic architecture between the two aquaculture strains of Arctic charr may be reflective of the differential evolutionary histories experienced by these fishes, and illustrates the importance of including different strains to investigate genetic variation in a species where the intent is to use that variation in selective breeding programs.

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.

Clock genes and their genomic distributions in three species of salmonid fishes: Associations with genes regulating sexual maturation and cell cycling

Background

Clock family genes encode transcription factors that regulate clock-controlled genes and thus regulate many physiological mechanisms/processes in a circadian fashion. Clock1 duplicates and copies of Clock3 and NPAS2-like genes were partially characterized (genomic sequencing) and mapped using family-based indels/SNPs in rainbow trout (RT)(Oncorhynchus mykiss), Arctic charr (AC)(Salvelinus alpinus), and Atlantic salmon (AS)(Salmo salar) mapping panels.

Results

Clock1 duplicates mapped to linkage groups RT-8/-24, AC-16/-13 and AS-2/-18. Clock3/NPAS2-like genes mapped to RT-9/-20, AC-20/-43, and AS-5. Most of these linkage group regions containing the Clock gene duplicates were derived from the most recent 4R whole genome duplication event specific to the salmonids. These linkage groups contain quantitative trait loci (QTL) for life history and growth traits (i.e., reproduction and cell cycling). Comparative synteny analyses with other model teleost species reveal a high degree of conservation for genes in these chromosomal regions suggesting that functionally related or co-regulated genes are clustered in syntenic blocks. For example, anti-müllerian hormone (amh), regulating sexual maturation, and ornithine decarboxylase antizymes (oaz1 and oaz2), regulating cell cycling, are contained within these syntenic blocks.

Conclusions

Synteny analyses indicate that regions homologous to major life-history QTL regions in salmonids contain many candidate genes that are likely to influence reproduction and cell cycling. The order of these genes is highly conserved across the vertebrate species examined, and as such, these genes may make up a functional cluster of genes that are likely co-regulated. CLOCK, as a transcription factor, is found within this block and therefore has the potential to cis-regulate the processes influenced by these genes. Additionally, clock-controlled genes (CCGs) are located in other life-history QTL regions within salmonids suggesting that at least in part, trans-regulation of these QTL regions may also occur via Clock expression.

Growth-related quantitative trait loci in domestic and wild rainbow trout (Oncorhynchus mykiss)

Background

Somatic growth is a complex process that involves the action and interaction of genes and environment. A number of quantitative trait loci (QTL) previously identified for body weight and condition factor in rainbow trout (Oncorhynchus mykiss), and two other salmonid species, were used to further investigate the genetic architecture of growth-influencing genes in this species. Relationships among previously mapped candidate genes for growth and their co-localization to identified QTL regions are reported. Furthermore, using a comparative genomic analysis of syntenic rainbow trout linkage group clusters to their homologous regions within model teleost species such as zebrafish, stickleback and medaka, inferences were made regarding additional possible candidate genes underlying identified QTL regions.

Results

Body weight (BW) QTL were detected on the majority of rainbow trout linkage groups across 10 parents from 3 strains. However, only 10 linkage groups (i.e., RT-3, -6, -8, -9, -10, -12, -13, -22, -24, -27) possessed QTL regions with chromosome-wide or genome-wide effects across multiple parents. Fewer QTL for condition factor (K) were identified and only six instances of co-localization across families were detected (i.e. RT-9, -15, -16, -23, -27, -31 and RT-2/9 homeologs). Of note, both BW and K QTL co-localize on RT-9 and RT-27. The incidence of epistatic interaction across genomic regions within different female backgrounds was also examined, and although evidence for interaction effects within certain QTL regions were evident, these interactions were few in number and statistically weak. Of interest, however, was the fact that these predominantly occurred within K QTL regions. Currently mapped growth candidate genes are largely congruent with the identified QTL regions. More QTL were detected in male, compared to female parents, with the greatest number evident in an F₁ male parent derived from an intercross between domesticated and wild strain of rainbow trout which differed strongly in growth rate.

Conclusions

Strain background influences the degree to which QTL effects are evident for growth-related genes. The process of domestication (which primarily selects faster growing fish) may largely reduce the genetic influences on growth-specific phenotypic variation. Although heritabilities have been reported to be relatively high for both BW and K growth traits, the genetic architecture of K phenotypic variation appears less defined (i.e., fewer major contributing QTL regions were identified compared with BW QTL regions).

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).

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.

Comparative genomics and evolution of conserved noncoding elements (CNE) in rainbow trout

Background

Recent advances in the accumulation of genetic mapping and DNA sequence information from several salmonid species support the long standing view of an autopolyploid origin of these fishes (i.e., 4R). However, the paralogy relationships of the chromosomal segments descendent from earlier polyploidization events (i.e., 2R/3R) largely remain unknown, mainly due to an unbalanced pseudogenization of paralogous genes that were once resident on the ancient duplicated segments. Inter-specific conserved noncoding elements (CNE) might hold the key in identifying these regions, if they are associated with arrays of genes that have been highly conserved in syntenic blocks through evolution. To test this hypothesis, we investigated the chromosomal positions of subset of CNE in the rainbow trout genome using a comparative genomic framework.

Results

Through a genome wide analysis, we selected 41 pairs of adjacent CNE located on various chromosomes in zebrafish and obtained their intervening, less conserved, sequence information from rainbow trout. We identified 56 distinct fragments corresponding to about 150 Kbp of sequence data that were localized to 67 different chromosomal regions in the rainbow trout genome. The genomic positions of many duplicated CNE provided additional support for some previously suggested homeologies in this species. Additionally, we now propose 40 new potential paralogous affinities by analyzing the variation in the segregation patterns of some multi-copy CNE along with the synteny association comparison using several model vertebrates. Some of these regions appear to carry signatures of the 1R, 2R or 3R duplications. A subset of these CNE markers also demonstrated high utility in identifying homologous chromosomal segments in the genomes of Atlantic salmon and Arctic charr.

Conclusion

CNE seem to be more efficacious than coding sequences in providing insights into the ancient paralogous affinities within the vertebrate genomes. Such a feature makes these elements extremely attractive for comparative genomics studies, as they can be treated as 'anchor' markers to investigate the association of distally located candidate genes on the homologous genomic segments of closely or distantly related organisms.

Determination of quantitative trait loci (QTL) for early maturation in rainbow trout (Oncorhynchus mykiss)

To identify quantitative trait loci (QTL) influencing early maturation (EM) in rainbow trout (Oncorhynchus mykiss), a genome scan was performed using 100 microsatellite loci across 29 linkage groups. Six inter-strain paternal half-sib families using three inter-strain F(1) brothers (approximately 50 progeny in each family) derived from two strains that differ in the propensity for EM were used in the study. Alleles derived from both parental sources were observed to contribute to the expression of EM in the progeny of the brothers. Four genome-wide significant QTL regions (i.e., RT-8, -17, -24, and -30) were observed. EM QTL detected on RT-8 and -24 demonstrated significant and suggestive QTL effects in both male and female progeny. Furthermore, within both male and female full-sib groupings, QTL on RT-8 and -24 were detected in two or more of the five parents used. Significant genome-wide and several strong chromosome-wide QTL for EM localized to different regions in males and females, suggesting some sex-specific control. Namely, QTL detected on RT-13, -15, -21, and -30 were associated with EM only in females, and those on RT-3, -17, and -19 were associated with EM only in males. Within the QTL regions identified, a comparison of syntenic EST markers from the rainbow trout linkage map with the zebrafish (Danio rerio) genome identified several putative candidate genes that may influence EM.

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.

Quantitative trait loci for body weight, condition factor and age at sexual maturation in Arctic charr (Salvelinus alpinus): comparative analysis with rainbow trout (Oncorhynchus mykiss) and Atlantic salmon (Salmo salar)

In salmonid fishes, life-history changes may often be coupled to early individual growth trajectories. We identified quantitative trait loci (QTL) for body weight (BW), condition factor (K) and age at sexual maturation (MT) in two full-sib families of Arctic charr (Salvelinus alpinus) to ascertain if QTL for MT were confounded with BW QTL intervals. Three significant QTL for BW, three QTL for MT and one significant QTL for K were identified. A BW QTL with major effect was localized to linkage group 8 (AC-8) and explained more than 34% of the phenotypic variation. Markers on AC-8 have previously been identified as being associated with variation in fork length and BW in this species. Similarly, a major QTL (PEV = 23%) with an influence on the female MT was localized to AC-23. Some of these regions are homologous to those in the genomes of rainbow trout (Oncorhynchus mykiss) and Atlantic salmon (Salmo salar), where similar QTL effects have been detected. Our results also suggest the conservation of MT QTL on the homeologous linkage group pair AC-3/24 in Arctic charr. We further identified chromosomal regions that harbor QTL for multiple traits. In particular, markers on AC-4, -20 and -36 had detectable QTL for all traits studied. Significant MT QTL detected on AC-23, -24, and -27 were autonomous of any BW QTL regions, suggesting that the regulation of MT may be more independent of BW control within this species than in other species of salmonids.

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.

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.

A linkage map for brown trout (Salmo trutta): chromosome homeologies and comparative genome organization with other salmonid fish

We report on the construction of a linkage map for brown trout (Salmo trutta) and its comparison with those of other tetraploid-derivative fish in the family Salmonidae, including Atlantic salmon (Salmo salar), rainbow trout (Oncorhynchus mykiss), and Arctic char (Salvelinus alpinus). Overall, we identified 37 linkage groups (2n = 80) from the analysis of 288 microsatellite polymorphisms, 13 allozyme markers, and phenotypic sex in four backcross families. Additionally, we used gene-centromere analysis to approximate the position of the centromere for 20 linkage groups and thus relate linkage arrangements to the physical morphology of chromosomes. Sex-specific maps derived from multiple parents were estimated to cover 346.4 and 912.5 cM of the male and female genomes, respectively. As previously observed in other salmonids, recombination rates showed large sex differences (average female-to-male ratio was 6.4), with male crossovers generally localized toward the distal end of linkage groups. Putative homeologous regions inherited from the salmonid tetraploid ancestor were identified for 10 pairs of linkage groups, including five chromosomes showing evidence of residual tetrasomy (pseudolinkage). Map alignments with orthologous regions in Atlantic salmon, rainbow trout, and Arctic char also revealed extensive conservation of syntenic blocks across species, which was generally consistent with chromosome divergence through Robertsonian translocations.

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.

A physical map of the genome of Atlantic salmon, Salmo salar

A physical map of the Atlantic salmon (Salmo salar) genome was generated based on HindIII fingerprints of a publicly available BAC (bacterial artificial chromosome) library constructed from DNA isolated from a Norwegian male. Approximately 11.5 haploid genome equivalents (185,938 clones) were successfully fingerprinted. Contigs were first assembled via FPC using high-stringency (1e-16), and then end-to-end joins yielded 4354 contigs and 37,285 singletons. The accuracy of the contig assembly was verified by hybridization and PCR analysis using genetic markers. A subset of the BACs in the library contained few or no HindIII recognition sites in their insert DNA. BglI digestion fragment patterns of these BACs allowed us to identify three classes: (1) BACs containing histone genes, (2) BACs containing rDNA-repeating units, and (3) those that do not have BglI recognition sites. End-sequence analysis of selected BACs representing these three classes confirmed the identification of the first two classes and suggested that the third class contained highly repetitive DNA corresponding to tRNAs and related sequences.

A comparative analysis of the rainbow trout genome with 2 other species of fish (Arctic charr and Atlantic salmon) within the tetraploid derivative Salmonidae family (subfamily: Salmoninae)

We updated the genetic map of rainbow trout (Oncorhynchus mykiss) for 2 outcrossed mapping panels, and used this map to assess the putative chromosome structure and recombination rate differences among linkage groups. We then used the rainbow trout sex-specific maps to make comparisons with 2 other ancestrally polyploid species of salmonid fishes, Arctic charr (Salvelinus alpinus) and Atlantic salmon (Salmo salar) to identify homeologous chromosome affinities within each species and ascertain homologous chromosome relationships among the species. Salmonid fishes exhibit a wide range of sex-specific differences in recombination rate, with some species having the largest differences for any vertebrate species studied to date. Our current estimate of female:male recombination rates in rainbow trout is 4.31:1. Chromosome structure and (or) size is associated with recombination rate differences between the sexes in rainbow trout. Linkage groups derived from presumptive acrocentric type chromosomes were observed to have much lower sex-specific differences in recombination rate than metacentric type linkage groups. Arctic charr is karyotypically the least derived species (i.e., possessing a high number of acrocentric chromosomes) and Atlantic salmon is the most derived (i.e., possessing a number of whole-arm fusions). Atlantic salmon have the largest female:male recombination ratio difference (i.e., 16.81:1) compared with rainbow trout, and Arctic charr (1.69:1). Comparisons of recombination rates between homologous segments of linkage groups among species indicated that when significant experiment-wise differences were detected (7/24 tests), recombination rates were generally higher in the species with a less-derived chromosome structure (6/7 significant comparisons). Greater similarity in linkage group syntenies were observed between Atlantic salmon and rainbow trout, suggesting their closer phylogenetic affinities, and most interspecific linkage group comparisons support a model that suggests whole chromosome arm translocations have occurred in the evolution of this group. However, some possible exceptions were detected and these findings are discussed in relation to their influence on segregation distortion patterns. We also report unusual meiotic segregation patterns in a female parent involving the duplicated (homeologous) linkage group pair 12/16 and discuss several models that may account for these patterns.Key words: linkage analysis, genetic markers, polyploidy, tetrasomic inheritance, segregation distortion, recombination rate.

Evidence for Hox Gene Duplication in Rainbow Trout (Oncorhynchus mykiss): A Tetraploid Model Species

We examined the genomic organization of Hox genes in rainbow trout (Oncorhynchus mykiss), a tetraploid teleost derivative species, in order to test models of presumptive genomic duplications during vertebrate evolution. Thirteen putative clusters were localized in the current rainbow trout genetic map; however, analysis of the sequence data suggests the presence of at least 14 Hox clusters. Many duplicated genes appear to have been retained in the genome and share a high percentage of amino acid similarity with one another. We characterized two Hox genes located within the HoxCb cluster that may have been lost independently in other teleost species studied to date. Finally, we identified conserved syntenic blocks between salmonids and human, and provide data supporting two new linkage group homeologies (i.e., RT-3/16, RT-12/29) and three previously described homeologies (RT-2/9, RT-17/22, and RT-27/31) in rainbow trout.

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.

Evolution of Hox clusters in Salmonidae: a comparative analysis between Atlantic salmon (Salmo salar) and rainbow trout (Oncorhynchus mykiss)

We studied the genomic organization of Hox genes in Atlantic salmon (Salmo salar) and made comparisons to that in rainbow trout (Oncorhynchus mykiss), another member of the family Salmonidae. We used these two species to test the hypothesis that the Hox genes would provide evidence for a fourth round of duplication (4R) of this gene family given the recent polyploid ancestry of the salmonid fish. Thirteen putative Hox clusters were identified and 10 of these complexes were localized to the current Atlantic salmon genetic map. Syntenic regions with the rainbow trout linkage map were detected and further homologies and homeologies are suggested. We propose that the common ancestor of Atlantic salmon and rainbow trout possessed at least 14 clusters of Hox genes, and additional clusters cannot be ruled out. Salmonid Hox cluster complements seem to be more similar to those of zebrafish (Danio rerio) than medaka (Oryzias latipes) or pufferfish (Sphoeroides nephelus and Takifugu rubripes), as both Atlantic salmon and rainbow trout have retained HoxCb ortholog, which has been lost in medaka and pufferfish but not in zebrafish. However, our data suggest that phylogenetically, the homologous genes within each cluster express mosaic relationships among the teleosts tested and, thus, leave unresolved the interfamilial relationships among these taxa.

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).

Characterization of Na, K-ATPase genes in Atlantic salmon (Salmo salar) and comparative genomic organization with rainbow trout (Oncorhynchus mykiss)

A combination of molecular and in silico approaches was employed to assemble a survey of Na, K-ATPase genes contained in the ancestrally tetraploid genome of the Atlantic salmon (Salmo salar). Molecular characterization of genomic clones coding for the alpha subunit revealed two single genes (alpha1a and alpha2) and two pairs of presumably homeologous genes (alpha1b/i-ii and alpha1c/i-ii). Each of the six genes showed high sequence similarity to isoforms previously isolated from rainbow trout and extensive structural differences relative to putative orthologs in the human genome. In silico analysis of expressed sequence tag (EST) collections indicated that at least five alpha (alpha1a, alpha1b, alpha1c, alpha2, and alpha3) and four beta (beta1a, beta1b, beta2, and beta3b) subunit isoforms are expressed in Atlantic salmon. Meiotic linkage analysis further showed that Na, K-ATPase genes are dispersed throughout the salmon genome, with the exception of two multigene clusters on linkage groups AS-22 and AS-28. Duplicate gene copies for the isoform alpha1b were assigned to linkage groups with multiple homeologous anchors (AS-22 and AS-23), while beta2 duplicates suggested a new homeologous affinity between AS-05 and AS-28. In addition, the comparison of linkage arrangements with rainbow trout also showed that the genomic organization of Na, K-ATPase genes is consistent with the evolutionary conservation of syntenic chromosome regions between these species.

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.

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.