PacBio Iso-Seq Improves the Rainbow Trout Genome Annotation and Identifies Alternative Splicing Associated With Economically Important Phenotypes

Rainbow trout is an important model organism that has received concerted international efforts to study the transcriptome. For this purpose, short-read sequencing has been primarily used over the past decade. However, these sequences are too short of resolving the transcriptome complexity. This study reported a first full-length transcriptome assembly of the rainbow trout using single-molecule long-read isoform sequencing (Iso-Seq). Extensive computational approaches were used to refine and validate the reconstructed transcriptome. The study identified 10,640 high-confidence transcripts not previously annotated, in addition to 1,479 isoforms not mapped to the current Swanson reference genome. Most of the identified lncRNAs were non-coding variants of coding transcripts. The majority of genes had multiple transcript isoforms (average ∼3 isoforms/locus). Intron retention (IR) and exon skipping (ES) accounted for 56% of alternative splicing (AS) events. Iso-Seq improved the reference genome annotation, which allowed identification of characteristic AS associated with fish growth, muscle accretion, disease resistance, stress response, and fish migration. For instance, an ES in GVIN1 gene existed in fish susceptible to bacterial cold-water disease (BCWD). Besides, under five stress conditions, there was a commonly regulated exon in prolyl 4-hydroxylase subunit alpha-2 (P4HA2) gene. The reconstructed gene models and their posttranscriptional processing in rainbow trout provide invaluable resources that could be further used for future genetics and genomics studies. Additionally, the study identified characteristic transcription events associated with economically important phenotypes, which could be applied in selective breeding.

A long reads-based de-novo assembly of the genome of the Arlee homozygous line reveals chromosomal rearrangements in rainbow trout

Currently, there is still a need to improve the contiguity of the rainbow trout reference genome and to use multiple genetic backgrounds that will represent the genetic diversity of this species. The Arlee doubled haploid line was originated from a domesticated hatchery strain that was originally collected from the northern California coast. The Canu pipeline was used to generate the Arlee line genome de-novo assembly from high coverage PacBio long-reads sequence data. The assembly was further improved with Bionano optical maps and Hi-C proximity ligation sequence data to generate 32 major scaffolds corresponding to the karyotype of the Arlee line (2 N = 64). It is composed of 938 scaffolds with N50 of 39.16 Mb and a total length of 2.33 Gb, of which ∼95% was in 32 chromosome sequences with only 438 gaps between contigs and scaffolds. In rainbow trout the haploid chromosome number can vary from 29 to 32. In the Arlee karyotype the haploid chromosome number is 32 because chromosomes Omy04, 14 and 25 are divided into six acrocentric chromosomes. Additional structural variations that were identified in the Arlee genome included the major inversions on chromosomes Omy05 and Omy20 and additional 15 smaller inversions that will require further validation. This is also the first rainbow trout genome assembly that includes a scaffold with the sex-determination gene (sdY) in the chromosome Y sequence. The utility of this genome assembly is shown through the improved annotation of the duplicated genome loci that harbor the IGH genes on chromosomes Omy12 and Omy13.

Identification of High-Confidence Structural Variants in Domesticated Rainbow Trout Using Whole-Genome Sequencing

Genomic structural variants (SVs) are a major source of genetic and phenotypic variation but have not been investigated systematically in rainbow trout (Oncorhynchus mykiss), an important aquaculture species of cold freshwater. The objectives of this study were 1) to identify and validate high-confidence SVs in rainbow trout using whole-genome re-sequencing; and 2) to examine the contribution of transposable elements (TEs) to SVs in rainbow trout. A total of 96 rainbow trout, including 11 homozygous lines and 85 outbred fish from three breeding populations, were whole-genome sequenced with an average genome coverage of 17.2×. Putative SVs were identified using the program Smoove which integrates LUMPY and other associated tools into one package. After rigorous filtering, 13,863 high-confidence SVs were identified. Pacific Biosciences long-reads of Arlee, one of the homozygous lines used for SV detection, validated 98% (3,948 of 4,030) of the high-confidence SVs identified in the Arlee homozygous line. Based on principal component analysis, the 85 outbred fish clustered into three groups consistent with their populations of origin, further indicating that the high-confidence SVs identified in this study are robust. The repetitive DNA content of the high-confidence SV sequences was 86.5%, which is much higher than the 57.1% repetitive DNA content of the reference genome, and is also higher than the repetitive DNA content of Atlantic salmon SVs reported previously. TEs thus contribute substantially to SVs in rainbow trout as TEs make up the majority of repetitive sequences. Hundreds of the high-confidence SVs were annotated as exon-loss or gene-fusion variants, and may have phenotypic effects. The high-confidence SVs reported in this study provide a foundation for further rainbow trout SV studies.

Publisher Correction: Sex-dependent dominance maintains migration supergene in rainbow trout

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Sex-dependent dominance maintains migration supergene in rainbow trout

Males and females often differ in their fitness optima for shared traits that have a shared genetic basis, leading to sexual conflict. Morphologically differentiated sex chromosomes can resolve this conflict and protect sexually antagonistic variation, but they accumulate deleterious mutations. However, how sexual conflict is resolved in species that lack differentiated sex chromosomes is largely unknown. Here we present a chromosome-anchored genome assembly for rainbow trout (Oncorhynchus mykiss) and characterize a 55-Mb double-inversion supergene that mediates sex-specific migratory tendency through sex-dependent dominance reversal, an alternative mechanism for resolving sexual conflict. The double inversion contains key photosensory, circadian rhythm, adiposity and sex-related genes and displays a latitudinal frequency cline, indicating environmentally dependent selection. Our results show sex-dependent dominance reversal across a large autosomal supergene, a mechanism for sexual conflict resolution capable of protecting sexually antagonistic variation while avoiding the homozygous lethality and deleterious mutations associated with typical heteromorphic sex chromosomes.

A New Single Nucleotide Polymorphism Database for Rainbow Trout Generated Through Whole Genome Resequencing

Single-nucleotide polymorphisms (SNPs) are highly abundant markers, which are broadly distributed in animal genomes. For rainbow trout (Oncorhynchus mykiss), SNP discovery has been previously done through sequencing of restriction-site associated DNA (RAD) libraries, reduced representation libraries (RRL) and RNA sequencing. Recently we have performed high coverage whole genome resequencing with 61 unrelated samples, representing a wide range of rainbow trout and steelhead populations, with 49 new samples added to 12 aquaculture samples from AquaGen (Norway) that we previously used for SNP discovery. Of the 49 new samples, 11 were double-haploid lines from Washington State University (WSU) and 38 represented wild and hatchery populations from a wide range of geographic distribution and with divergent migratory phenotypes. We then mapped the sequences to the new rainbow trout reference genome assembly (GCA_002163495.1) which is based on the Swanson YY doubled haploid line. Variant calling was conducted with FreeBayes and SAMtools mpileup, followed by filtering of SNPs based on quality score, sequence complexity, read depth on the locus, and number of genotyped samples. Results from the two variant calling programs were compared and genotypes of the double haploid samples were used for detecting and filtering putative paralogous sequence variants (PSVs) and multi-sequence variants (MSVs). Overall, 30,302,087 SNPs were identified on the rainbow trout genome 29 chromosomes and 1,139,018 on unplaced scaffolds, with 4,042,723 SNPs having high minor allele frequency (MAF > 0.25). The average SNP density on the chromosomes was one SNP per 64 bp, or 15.6 SNPs per 1 kb. Results from the phylogenetic analysis that we conducted indicate that the SNP markers contain enough population-specific polymorphisms for recovering population relationships despite the small sample size used. Intra-Population polymorphism assessment revealed high level of polymorphism and heterozygosity within each population. We also provide functional annotation based on the genome position of each SNP and evaluate the use of clonal lines for filtering of PSVs and MSVs. These SNPs form a new database, which provides an important resource for a new high density SNP array design and for other SNP genotyping platforms used for genetic and genomics studies of this iconic salmonid fish species.

Increased mitochondrial DNA diversity in ancient Columbia River basin Chinook salmon Oncorhynchus tshawytscha

The Columbia River and its tributaries provide essential spawning and rearing habitat for many salmonid species, including Chinook salmon (Oncorhynchus tshawytscha). Chinook salmon were historically abundant throughout the basin and Native Americans in the region relied heavily on these fish for thousands of years. Following the arrival of Europeans in the 1800s, salmon in the basin experienced broad declines linked to overfishing, water diversion projects, habitat destruction, connectivity reduction, introgression with hatchery-origin fish, and hydropower development. Despite historical abundance, many native salmonids are now at risk of extinction. Research and management related to Chinook salmon is usually explored under what are termed "the four H's": habitat, harvest, hatcheries, and hydropower; here we explore a fifth H, history. Patterns of prehistoric and contemporary mitochondrial DNA variation from Chinook salmon were analyzed to characterize and compare population genetic diversity prior to recent alterations and, thus, elucidate a deeper history for this species. A total of 346 ancient and 366 contemporary samples were processed during this study. Species was determined for 130 of the ancient samples and control region haplotypes of 84 of these were sequenced. Diversity estimates from these 84 ancient Chinook salmon were compared to 379 contemporary samples. Our analysis provides the first direct measure of reduced genetic diversity for Chinook salmon from the ancient to the contemporary period, as measured both in direct loss of mitochondrial haplotypes and reductions in haplotype and nucleotide diversity. However, these losses do not appear equal across the basin, with higher losses of diversity in the mid-Columbia than in the Snake subbasin. The results are unexpected, as the two groups were predicted to share a common history as parts of the larger Columbia River Basin, and instead indicate that Chinook salmon in these subbasins may have divergent demographic histories.

Migration-related phenotypic divergence is associated with epigenetic modifications in rainbow trout

Migration is essential for the reproduction and survival of many animals, yet little is understood about its underlying molecular mechanisms. We used the salmonid Oncorhynchus mykiss to gain mechanistic insight into smoltification, which is a morphological, physiological and behavioural transition undertaken by juveniles in preparation for seaward migration. O. mykiss is experimentally tractable and displays intra- and interpopulation variation in migration propensity. Migratory individuals can produce nonmigratory progeny and vice versa, indicating a high degree of phenotypic plasticity. One potential way that phenotypic plasticity might be linked to variation in migration-related life history tactics is through epigenetic regulation of gene expression. To explore this, we quantitatively measured genome-scale DNA methylation in fin tissue using reduced representation bisulphite sequencing of F2 siblings produced from a cross between steelhead (migratory) and rainbow trout (nonmigratory) lines. We identified 57 differentially methylated regions (DMRs) between smolt and resident O. mykiss juveniles. DMRs were high in magnitude, with up to 62% differential methylation between life history types, and over half of the gene-associated DMRs were in transcriptional regulatory regions. Many of the DMRs encode proteins with activity relevant to migration-related transitions (e.g. circadian rhythm pathway, nervous system development, protein kinase activity). This study provides the first evidence of a relationship between epigenetic variation and life history divergence associated with migration-related traits in any species.

Differential Expression of Genes that Control Respiration Contribute to Thermal Adaptation in Redband Trout (Oncorhynchus mykiss gairdneri)

Organisms can adapt to local environmental conditions as a plastic response or become adapted through natural selection on genetic variation. The ability to adapt to increased water temperatures will be of paramount importance for many fish species as the climate continues to warm and water resources become limited. Because increased water temperatures will reduce the dissolved oxygen available for fish, we hypothesized that adaptation to low oxygen environments would involve improved respiration through oxidative phosphorylation (OXPHOS). To test this hypothesis, we subjected individuals from two ecologically divergent populations of inland (redband) rainbow trout (Oncorhynchus mykiss gairdneri) with historically different temperature regimes (desert and montane) and their F1 progeny to diel cycles of temperature stress and then examined gene expression data for 80 nuclear- and mitochondrial-encoded OXPHOS subunits that participate in respiration. Of the 80 transcripts, 7 showed ≥ 2-fold difference in expression levels in gill tissue from desert fish under heat stress whereas the montane fish had none and the F1 only had one differentially expressed gene. A structural analysis of the proteins encoded by those genes suggests that the response could coordinate the formation of supercomplexes and oligomers. Supercomplexes may increase the efficiency of respiration because complexes I, III, and IV are brought into close proximity and oligomerization of complex V alters the macrostructure of mitochondria to improve respiration. Significant differences in gene expression patterns in response to heat stress in a common environment indicate that the response was not due to plasticity but had a genetic basis.

Transcriptome assembly, gene annotation and tissue gene expression atlas of the rainbow trout

Efforts to obtain a comprehensive genome sequence for rainbow trout are ongoing and will be complemented by transcriptome information that will enhance genome assembly and annotation. Previously, transcriptome reference sequences were reported using data from different sources. Although the previous work added a great wealth of sequences, a complete and well-annotated transcriptome is still needed. In addition, gene expression in different tissues was not completely addressed in the previous studies. In this study, non-normalized cDNA libraries were sequenced from 13 different tissues of a single doubled haploid rainbow trout from the same source used for the rainbow trout genome sequence. A total of ~1.167 billion paired-end reads were de novo assembled using the Trinity RNA-Seq assembler yielding 474,524 contigs > 500 base-pairs. Of them, 287,593 had homologies to the NCBI non-redundant protein database. The longest contig of each cluster was selected as a reference, yielding 44,990 representative contigs. A total of 4,146 contigs (9.2%), including 710 full-length sequences, did not match any mRNA sequences in the current rainbow trout genome reference. Mapping reads to the reference genome identified an additional 11,843 transcripts not annotated in the genome. A digital gene expression atlas revealed 7,678 housekeeping and 4,021 tissue-specific genes. Expression of about 16,000–32,000 genes (35–71% of the identified genes) accounted for basic and specialized functions of each tissue. White muscle and stomach had the least complex transcriptomes, with high percentages of their total mRNA contributed by a small number of genes. Brain, testis and intestine, in contrast, had complex transcriptomes, with a large numbers of genes involved in their expression patterns. This study provides comprehensive de novo transcriptome information that is suitable for functional and comparative genomics studies in rainbow trout, including annotation of the genome.

Aggressive behavior, brain size and domestication in clonal rainbow trout lines

Domestication causes behavior and brain size changes in many species. We addressed three questions using clonal rainbow trout lines: What are the mirror-elicited aggressive tendencies in lines with varying degrees of domestication? How does brain size relate to genotype and domestication level? Finally, is there a relationship between aggressive behavior and brain size? Clonal lines, although sampling a limited subset of the species variation, provide us with a reproducible experimental system with which we can develop hypotheses for further research. We performed principal component analyses on 12 continuous behavior and brain/body size variables and one discrete behavioral variable ("yawn") and detected several aggression syndromes. Two behaviors, "freeze" and "escape", associated with high domestication; "display" and "yawn" behavior associated with wild lines and "swim against the mirror" behavior associated with semi-wild and domestic lines. Two brain size traits, total brain and olfactory volume, were significantly related to domestication level when taking total body size into account, with domesticated lines having larger total brain volume and olfactory regions. The aggression syndromes identified indicate that future QTL mapping studies on domestication-related traits would likely be fruitful.

Characterization of the rainbow trout spleen transcriptome and identification of immune-related genes

Resistance against diseases affects profitability of rainbow trout. Limited information is available about functions and mechanisms of teleost immune pathways. Immunogenomics provides powerful tools to determine disease resistance genes/gene pathways and develop genetic markers for genomic selection. RNA-Seq sequencing of the rainbow trout spleen yielded 93,532,200 reads (100 bp). High quality reads were assembled into 43,047 contigs. 26,333 (61.17%) of the contigs had hits to the NR protein database and 7024 (16.32%) had hits to the KEGG database. Gene ontology showed significant percentages of transcripts assigned to binding (51%), signaling (7%), response to stimuli (9%) and receptor activity (4%) suggesting existence of many immune-related genes. KEGG annotation revealed 2825 sequences belonging to "organismal systems" with the highest number of sequences, 842 (29.81%), assigned to immune system. A number of sequences were identified for the first time in rainbow trout belonging to Toll-like receptor signaling (35), B cell receptor signaling pathway (44), T cell receptor signaling pathway (56), chemokine signaling pathway (73), Fc gamma R-mediated phagocytosis (52), leukocyte transendothelial migration (60) and NK cell mediated cytotoxicity (42). In addition, 51 transcripts were identified as spleen-specific genes. The list includes 277 full-length cDNAs. The presence of a large number of immune-related genes and pathways similar to other vertebrates suggests that innate and adaptive immunity in fish are conserved. This study provides deep-sequence data of rainbow trout spleen transcriptome and identifies many new immune-related genes and full-length cDNAs. This data will help identify allelic variations suitable for genomic selection and genetic manipulation in aquaculture.

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

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

Mapping and Expression of Candidate Genes for Development Rate in Rainbow Trout (Oncorhynchus mykiss)

Development rate has important implications for individual fitness and physiology. In salmonid fishes, development rate correlates with many traits later in life, including life-history diversity, growth, and age and size at sexual maturation. In rainbow trout (Oncorhynchus mykiss), a quantitative trait locus for embryonic development rate has been detected on chromosome 5 across populations. However, few candidate genes have been identified within this region. In this study, we use gene mapping, gene expression, and quantitative genetic methods to further identify the genetic basis of embryonic developmental rate in O. mykiss Among the genes located in the region of the major development rate quantitative trait locus (GHR1, Clock1a, Myd118-1, and their paralogs), all were expressed early in embryonic development (fertilization through hatch), but none were differentially expressed between individuals with the fast- or slow-developing alleles for a major embryonic development rate quantitative trait locus. In a follow-up study of migratory and resident rainbow trout from natural populations in Alaska, we found significant additive variation in development rate and, moreover, found associations between development rate and allelic variation in all 3 candidate genes within the quantitative trait locus for embryonic development. The mapping of these genes to this region and associations in multiple populations provide positional candidates for further study of their roles in growth, development, and life-history diversity in this model salmonid.

Identification of single nucleotide polymorphisms from the transcriptome of an organism with a whole genome duplication

BACKGROUND

The common ancestor of salmonid fishes, including rainbow trout (Oncorhynchus mykiss), experienced a whole genome duplication between 20 and 100 million years ago, and many of the duplicated genes have been retained in the trout genome. This retention complicates efforts to detect allelic variation in salmonid fishes. Specifically, single nucleotide polymorphism (SNP) detection is problematic because nucleotide variation can be found between the duplicate copies (paralogs) of a gene as well as between alleles.

RESULTS

We present a method of differentiating between allelic and paralogous (gene copy) sequence variants, allowing identification of SNPs in organisms with multiple copies of a gene or set of genes. The basic strategy is to: 1) identify windows of unique cDNA sequences with homology to each other, 2) compare these unique cDNAs if they are not shared between individuals (i.e. the cDNA is homozygous in one individual and homozygous for another cDNA in the other individual), and 3) give a "SNP score" value between zero and one to each candidate sequence variant based on six criteria. Using this strategy we were able to detect about seven thousand potential SNPs from the transcriptomes of several clonal lines of rainbow trout. When directly compared to a pre-validated set of SNPs in polyploid wheat, we were also able to estimate the false-positive rate of this strategy as 0 to 28% depending on parameters used.

CONCLUSIONS

This strategy has an advantage over traditional techniques of SNP identification because another dimension of sequencing information is utilized. This method is especially well suited for identifying SNPs in polyploids, both outbred and inbred, but would tend to be conservative for diploid organisms.

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

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

Domestication shapes morphology in rainbow trout Oncorhynchus mykiss

In this study, clonal lines from North American resident and migratory populations of rainbow trout Oncorhynchus mykiss adapted to different geographical conditions and with different domestication histories were characterized morphologically. Lines reared in a common-garden experiment were characterized for external shape and meristic values, searching for a general pattern of morphological variation due to exposure to captive conditions. A sharp distinction was identified between wild and captive lines. The body profile was deeper in captive lines, with longer dorsal and anal fins and shorter and deeper caudal peduncles. Highly significant differences were also identified in meristic values among the lines but no consistent relation between meristic values and domestication status was detected. This morphological characterization will facilitate the selection of lines with divergent phenotypes for subsequent quantitative trait loci analysis, aimed at identifying genome regions linked with morphological adaptive response to captive conditions.

Y chromosome phylogeny for cutthroat trout (Oncorhynchus clarkii) subspecies is generally concordant with those of other markers

Sequence divergence was evaluated in the non-recombining, male-specific OmyY1 region of the Y chromosome among the subspecies of cutthroat trout (Oncorhynchus clarkii) in the western United States. This evaluation identified subspecies-discriminating OmyY1-haplotypes within a ∼1200bp region of the OmyY1 locus and localized the region to the end of the Y chromosome by FISH analysis. OmyY1 sequences were aligned and used to reconstruct a phylogeny of the cutthroat trout subspecies and related species via maximum-parsimony and Bayesian analyses. In the Y-haplotype phylogeny, clade distributions generally corresponded to the geographic distributions of the recognized subspecies. This phylogeny generally corresponded to a mitochondrial tree obtained for these subspecies in a previous study. Both support a clade of trout vs. Pacific salmon, of rainbow trout, and of a Yellowstone cutthroat group within the cutthroat trout. In our OmyY1 tree, however, the cutthroat "clade", although present topologically, was not statistically significant. Some key differences were found between trees obtained from the paternally-inherited OmyY1 vs. maternally-inherited mitochondrial haplotypes in cutthroat trout compared to rainbow trout. Other findings are: The trout OmyY1 region evolves between 3 and 13 times slower than the trout mitochondrial regions that have been studied. The Lahontan cutthroat trout had a fixed OmyY1 sequence throughout ten separate populations, suggesting this subspecies underwent a severe population bottleneck prior to its current dispersal throughout the Great Basin during the pluvial phase of the last ice age. The Yellowstone group is the most derived among the cutthroat trout and consists of the Yellowstone, Bonneville, Colorado, Rio Grande and greenback subspecies. Identification of subspecies and sex with this Y-chromosome marker may prove useful in conservation efforts.

Analysis of host genetic diversity and viral entry as sources of between-host variation in viral load

Little is known about the factors that drive the high levels of between-host variation in pathogen burden that are frequently observed in viral infections. Here, two factors thought to impact viral load variability, host genetic diversity and stochastic processes linked with viral entry into the host, were examined. This work was conducted with the aquatic vertebrate virus, Infectious hematopoietic necrosis virus (IHNV), in its natural host, rainbow trout. It was found that in controlled in vivo infections of IHNV, a suggestive trend of reduced between-fish viral load variation was observed in a clonal population of isogenic trout compared to a genetically diverse population of out-bred trout. However, this trend was not statistically significant for any of the four viral genotypes examined, and high levels of fish-to-fish variation persisted even in the isogenic trout population. A decrease in fish-to-fish viral load variation was also observed in virus injection challenges that bypassed the host entry step, compared to fish exposed to the virus through the natural water-borne immersion route of infection. This trend was significant for three of the four virus genotypes examined and suggests host entry may play a role in viral load variability. However, high levels of viral load variation also remained in the injection challenges. Together, these results indicate that although host genetic diversity and viral entry may play some role in between-fish viral load variation, they are not major factors. Other biological and non-biological parameters that may influence viral load variation are discussed.

A conserved haplotype controls parallel adaptation in geographically distant salmonid populations

Salmonid fishes exhibit extensive local adaptations owing to abundant environmental variation and precise natal homing. This extensive local adaptation makes conservation and restoration of salmonids a challenge. For example, defining unambiguous units of conservation is difficult, and restoration attempts often fail owing to inadequate adaptive matching of translocated populations. A better understanding of the genetic architecture of local adaptation in salmonids could provide valuable information to assist in conserving and restoring natural populations of these important species. Here, we use a combination of laboratory crosses and next-generation sequencing to investigate the genetic architecture of the parallel adaptation of rapid development rate in two geographically and genetically distant populations of rainbow trout (Oncorhynchus mykiss). Strikingly, we find that not only is a parallel genetic mechanism used but that a conserved haplotype is responsible for this intriguing adaptation. The repeated use of adaptive genetic variation across distant geographical areas could be a general theme in salmonids and have important implications for conservation and restoration.

Differential gene expression in male and female rainbow trout embryos prior to the onset of gross morphological differentiation of the gonads

BACKGROUND

There are large differences between the sexes at the genetic level; these differences include heterogametic sex chromosomes and/or differences in expression of genes between the sexes. In rainbow trout (Oncorhynchus mykiss) qRT-PCR studies have found significant differences in expression of several candidate sex determining genes. However, these genes represent a very small fraction of the genome and research in other species suggests there are large portions of the transcriptome that are differentially expressed between the sexes. These differences are especially noticeable once gonad differentiation and maturation has occurred, but less is known at earlier stages of development. Here we use data from a microarray and qRT-PCR to identify genes differentially expressed between the sexes at three time points in pre-hatch embryos, prior to the known timing of sexual differentiation in this species.

RESULTS

The microarray study revealed 883 differentially expressed features between the sexes with roughly equal numbers of male and female upregulated features across time points. Most of the differentially expressed genes on the microarray were not related to sex function, suggesting large scale differences in gene expression between the sexes are present early in development. Candidate gene analysis revealed sox9, DMRT1, Nr5a1 and wt1 were upregulated in males at some time points and foxl2, ovol1, fst and cyp19a1a were upregulated in females at some time points.

CONCLUSION

This is the first study to identify sexual dimorphism in expression of the genome during embryogenesis in any fish and demonstrates that transcriptional differences are present before the completion of gonadogenesis.

Transcriptome profiling of embryonic development rate in rainbow trout advanced backcross introgression lines

In rainbow trout (Oncorhynchus mykiss) and other fishes, embryonic development rate is an ecologically and evolutionarily important trait that is closely associated with survival and physiological performance later in life. To identify genes differentially regulated in fast and slow-developing embryos of rainbow trout, we examined gene expression across developmental time points in rainbow trout embryos possessing alleles linked to a major quantitative trait loci (QTL) for fast versus slow embryonic development rate. Whole genome expression microarray analyses were conducted using embryos from a fourth generation backcross family, whereby each backcross generation involved the introgression of the fast-developing alleles for a major development rate QTL into a slow-developing clonal line of rainbow trout. Embryos were collected at 15, 19, and 28 days post-fertilization; sex and QTL genotype were determined using molecular markers, and cDNA from 48 embryos were used for microarray analysis. A total of 183 features were identified with significant differences between embryonic development rate genotypes. Genes associated with cell cycle growth, muscle contraction and protein synthesis were expressed significantly higher in embryos with the fast-developing allele (Clearwater) than those with the slow-developing allele (Oregon State University), which may associate with fast growth and early body mass construction in embryo development. Across time points, individuals with the fast-developing QTL allele appeared to have earlier onset of these developmental processes when compared to individuals with the slow development alleles, even as early as 15 days post-fertilization. Differentially expressed candidate genes chosen for linkage mapping were localized primarily to regions outside of the major embryonic development rate QTL, with the exception of a single gene (very low-density lipoprotein receptor precursor).

Characterization of the rainbow trout transcriptome using Sanger and 454-pyrosequencing approaches

Background

Rainbow trout are important fish for aquaculture and recreational fisheries and serves as a model species for research investigations associated with carcinogenesis, comparative immunology, toxicology and evolutionary biology. However, to date there is no genome reference sequence to facilitate the development of molecular technologies that utilize high-throughput characterizations of gene expression and genetic variation. Alternatively, transcriptome sequencing is a rapid and efficient means for gene discovery and genetic marker development. Although a large number (258,973) of EST sequences are publicly available, the nature of rainbow trout duplicated genome hinders assembly and complicates annotation.

Results

High-throughput deep sequencing of the Swanson rainbow trout doubled-haploid transcriptome using 454-pyrosequencing technology yielded ~1.3 million reads with an average length of 344 bp, a total of 447 million bases. De novo assembly of the sequences yielded 151,847 Tentative Consensus (TC) sequences (average length of 662 bp) and 224,391 singletons. A combination assembly of both the 454-pyrosequencing ESTs and the pre-existing sequences resulted in 161,818 TCs (average length of 758 bp) and 261,071 singletons. Gene Ontology analysis of the combination assembly showed high similarities to transcriptomes of other fish species with known genome sequences.

Conclusion

The 454 library significantly increased the suite of ESTs available for rainbow trout, allowing improved assembly and annotation of the transcriptome. Furthermore, the 454 sequencing enables functional genome research in rainbow trout, providing a wealth of sequence data to serve as a reference transcriptome for future studies including identification of paralogous sequences and/or allelic variation, digital gene expression and proteomic research.

Deep divergence and apparent sex-biased dispersal revealed by a Y-linked marker in rainbow trout

Y-chromosome and mitochondrial DNA markers can reveal phylogenetic patterns by allowing tracking of male and female lineages, respectively. We used sequence data from a recently discovered Y-linked marker and a mitochondrial marker to examine phylogeographic structure in the widespread and economically important rainbow trout (Oncorhynchus mykiss). Two distinct geographic groupings that generally correspond to coastal and inland subspecies were evident within the Y-marker network while the mtDNA haplotype network showed little geographic structure. Our results suggest that male-specific behavior has prevented widespread admixture of Y haplotypes and that gene flow between the coastal and inland subspecies has largely occurred through females. This new Y marker may also aid conservation efforts by genetically identifying inland populations that have not hybridized with widely stocked coastal-derived hatchery fish.

A first generation BAC-based physical map of the rainbow trout genome

Background

Rainbow trout (Oncorhynchus mykiss) are the most-widely cultivated cold freshwater fish in the world and an important model species for many research areas. Coupling great interest in this species as a research model with the need for genetic improvement of aquaculture production efficiency traits justifies the continued development of genomics research resources. Many quantitative trait loci (QTL) have been identified for production and life-history traits in rainbow trout. A bacterial artificial chromosome (BAC) physical map is needed to facilitate fine mapping of QTL and the selection of positional candidate genes for incorporation in marker-assisted selection (MAS) for improving rainbow trout aquaculture production. This resource will also facilitate efforts to obtain and assemble a whole-genome reference sequence for this species.

Results

The physical map was constructed from DNA fingerprinting of 192,096 BAC clones using the 4-color high-information content fingerprinting (HICF) method. The clones were assembled into physical map contigs using the finger-printing contig (FPC) program. The map is composed of 4,173 contigs and 9,379 singletons. The total number of unique fingerprinting fragments (consensus bands) in contigs is 1,185,157, which corresponds to an estimated physical length of 2.0 Gb. The map assembly was validated by 1) comparison with probe hybridization results and agarose gel fingerprinting contigs; and 2) anchoring large contigs to the microsatellite-based genetic linkage map.

Conclusion

The production and validation of the first BAC physical map of the rainbow trout genome is described in this paper. We are currently integrating this map with the NCCCWA genetic map using more than 200 microsatellites isolated from BAC end sequences and by identifying BACs that harbor more than 300 previously mapped markers. The availability of an integrated physical and genetic map will enable detailed comparative genome analyses, fine mapping of QTL, positional cloning, selection of positional candidate genes for economically important traits and the incorporation of MAS into rainbow trout breeding programs.

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

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

Universal mtDNA primers for species identification of degraded bony fish samples

We developed primers for amplifying and sequencing highly degraded mtDNA from diverse fish species. The primers flank a variable 148-bp fragment within the 12S region of mtDNA. We screened and sequenced 82 samples of bony fishes representing 17 families to confirm cross-species amplification and identification. Salmonid species were analysed and demonstrate 13 species-specific SNPs within this region. Based on alignments of additional deposited sequences, these primers are conserved in many other species, making them useful for species identification using degraded DNA samples such as archaeological specimens.

Characterization of the interferon genes in homozygous rainbow trout reveals two novel genes, alternate splicing and differential regulation of duplicated genes

The genes encoding the type I and type II interferons (IFNs) have previously been identified in rainbow trout and their proteins partially characterized. These previous studies reported a single type II IFN (rtIFN-gamma) and three rainbow trout type I IFN genes that are classified into either group I (rtIFN1, rtIFN2) or group II (rtIFN3). In this present study, we report the identification of a novel IFN-gamma gene (rtIFN-gamma2) and a novel type I group II IFN (rtIFN4) in homozygous rainbow trout and predict that additional IFN genes or pseudogenes exist in the rainbow trout genome. Additionally, we provide evidence that short and long forms of rtIFN1 are actively and differentially transcribed in homozygous trout, and likely arose due to alternate splicing of the first exon. Quantitative reverse transcriptase PCR (qRT-PCR) assays were developed to systematically profile all of the rainbow trout IFN transcripts, with high specificity at an individual gene level, in naïve fish and after stimulation with virus or viral-related molecules. Cloned PCR products were used to ensure the specificity of the qRT-PCR assays and as absolute standards to assess transcript abundance of each gene. All IFN genes were modulated in response to Infectious hematopoietic necrosis virus (IHNV), a DNA vaccine based on the IHNV glycoprotein, and poly I:C. The most inducible of the type I IFN genes, by all stimuli tested, were rtIFN3 and the short transcript form of rtIFN1. Gene expression of rtIFN-gamma1 and rtIFN-gamma2 was highly up-regulated by IHNV infection and DNA vaccination but rtIFN-gamma2 was induced to a greater magnitude. The specificity of the qRT-PCR assays reported here will be useful for future studies aimed at identifying which cells produce IFNs at early time points after infection.

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

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

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

BACKGROUND

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

RESULTS

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

CONCLUSION

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

Genome assembly anchored QTL map of bovine chromosome 14

Bovine chromosome 14 (BTA14) has been widely explored for quantitative trait loci (QTL) and genes related to economically important traits in both dairy and beef cattle. We reviewed more than 40 investigations and anchored 126 QTL to the current genome assembly (Btau 4_0). Using this anchored QTL map, we observed that, in dairy cattle, the region spanning 0 – 10 Mb on BTA14 has the highest density QTL map with a total of 56 QTL, mainly for milk production traits. It is very likely that both somatic cell score (SCS) and clinical mastitis share some common QTL in two regions: 61.48 Mb - 73.84 Mb and 7.86 Mb – 39.55 Mb, respectively. As well, both ovulation rate and twinning rate might share a common QTL region from 34.16 Mb to 65.38 Mb. However, there are no common QTL locations in three pregnancy related phenotypes: non-return rate, pregnancy rate and daughter pregnancy rate. In beef cattle, the majority of QTL are located in a broad region of 15 Mb – 45 Mb on the chromosome. Functional genes, such as CRH, CYP11B1, DGAT1, FABP4 and TG, as potential candidates for some of these QTL, were also reviewed. Therefore, our review provides a standardized QTL map anchored within the current genome assembly, which would enhance the process of selecting positional and physiological candidate genes for many important traits in cattle.

Y-specific sequences and polymorphisms in rainbow trout and Chinook salmon

Improved methods for genetically sexing salmonids and for characterization of Y-chromosome homologies between species can contribute to understanding the evolution of sex chromosomes and sex-determining mechanisms. In this study we have explored 12.5 kb of Y-chromosome-specific sequence flanking the previously described OtY2 locus in Chinook salmon (Oncorhynchus tshawytscha) and 21 kb of homologous rainbow trout (Oncorhynchus mykiss) Y-chromosome-specific sequence. This is the first confirmed Y-specific sequence for rainbow trout. New Y-specific markers are described for Chinook salmon (OtY3) and rainbow trout (OmyY1), which are readily detected by PCR assays and are advantageous because they also produce autosomal control amplification products. Additionally, AFLP analysis of Chinook salmon yielded another potential Y-chromosome marker. These descriptions will facilitate genotypic sexing and should be useful for population studies of Y-chromosome polymorphisms and for future studies to characterize what appears to be a common sex-determining mechanism between these species.

The cutthroat trout Y chromosome is conserved with that of rainbow trout

Five genetic markers previously shown to be located on the sex chromosomes of rainbow trout (Oncorhynchus mykiss) were tested for linkage with the sex locus of Yellowstone cutthroat trout (Oncorhynchus clarki bouvieri) in a genetic cross created from a rainbow x cutthroat male hybrid. We show that the sex locus of both rainbow and cutthroat trout is on the same homologous linkage group. Fluorescence in situ hybridization (FISH) using a probe for the microsatellite marker Omm1665, which maps close to the sex locus of Yellowstone cutthroat trout, was used to identify the Y chromosome of cutthroat trout in the hybrid. The Y chromosome of cutthroat trout is sub-telocentric and lacks a DAPI band found on the short arm of the Y chromosome of some rainbow trout males.

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

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

A simplified QTL mapping approach for screening and mapping of novel AFLP markers associated with beef marbling

Genome screening of quantitative trait loci (QTL) for a complex trait is usually costly and highly laborious, as it requires a large number of markers spanning the whole genome. Here we present a simplified approach for screening and mapping of QTL-linked markers for beef marbling using a WagyuxLimousin F(2) reference population. This simplified approach involves integration of the amplified fragment length polymorphism (AFLP) with DNA pooling and selective genotyping and comparative bioinformatics tools. AFLP analysis on two high and two low marbling DNA pools yielded ten visually different markers. Among them, four were confirmed based on individual AFLP validation. Sequencing and in silico characterization assigned two of these AFLP markers to bovine chromosomes 1 (BTA1) and 13 (BTA13), which are orthologous to human chromosomes HSA21q22.2 and HSA10p11.23 with both regions harboring QTL for obesity-related phenotypes. Both AFLP markers showed significantly large additive genetic effects (0.28+/-0.11 on BTA1 and 0.54+/-0.21 on BTA13) on beef-marbling score (BMS) (P<0.05). Overall, this approach is less time consuming, inexpensive and in particular, suitable for screening and mapping QTL-linked markers when targeting one or a few complex traits.

Quantitative trait loci x maternal cytoplasmic environment interaction for development rate in Oncorhynchus mykiss

Effects of maternal cytoplasmic environment (MCE) on development rate in rainbow trout were evaluated within a quantitative trait loci (QTL) analysis framework. Previous research had identified QTL for development rate in doubled haploid (DH) progeny produced from a cross between the Oregon State University (OSU) and the Swanson (SW) River rainbow trout clonal lines. In this study, progeny for QTL mapping were produced from a cross between the OSU and Clearwater (CW) River clonal lines. Doubled haploids were produced from the OSU x CW F1 by androgenesis using eggs from different females (or MCEs); with androgenesis, the maternal nuclear genome was destroyed by irradiation and diploidy was restored by blocking the first embryonic cleavage by heat shock. All embryos were incubated at the same temperature and development rate quantified as time to hatch. Using a linkage map constructed primarily with AFLP markers, QTL mapping was performed, including MCE covariates and QTL x MCE effects in models for testing. The major QTL for development rate in the OSU x SW cross overlaps with the major QTL found in this OSU x CW cross; effects at this locus were the same across MCEs. Both MCE and QTL x MCE effects contribute to variability in development rate, but QTL x MCE were minor and detected only at small-effect QTL.

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

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

Assignment of rainbow trout linkage groups to specific chromosomes

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

A novel nuclear-encoded mitochondrial poly(A) polymerase PAPD1 is a potential candidate gene for the extreme obesity related phenotypes in mammals

People with obesity, especially extreme obesity, are at risk for many health problems. However, the responsible genes remain unknown in >95% of severe obesity cases. Our previous genome-wide scan of Wagyu x Limousin F₂ cattle crosses with extreme phenotypes revealed a molecular marker significantly associated with intramuscular fat deposition. Characterization of this marker showed that it is orthologous to the human gene KIAA1462 located on HSA10p11.23, where a major quantitative trait locus for morbid obesity has been reported. The newly identified mitochondrial poly(A) polymerase associated domain containing 1 (PAPD1) gene, which is located near this marker, is particularly interesting because the polymerase is required for the polyadenylation and stabilization of mammalian mitochondrial mRNAs. In the present study, both cDNA and genomic DNA sequences were annotated for the bovine PAPD1 gene and ten genetic markers were detected in the promoter and exon 1 region. Among seven markers assayed on ~ 250 Wagyu x Limousin F₂ animals, two single nucleotide polymorphisms (SNPs) in the promoter region were significantly associated with intramuscular fat (P<0.05). However, there was a significant interaction (P<0.05) between a third SNP, which causes an amino acid change in coding exon 1, and each of these two promoter SNPs on intramuscular fat deposition. In particular, the differences between double heterozygous animals at two polymorphic sites and the slim genotype animals exceeded 2.3 standard deviations for the trait in both cases. Our study provides evidence for a new mechanism – the involvement of compound heterosis in extreme obesity, which warrants further examination.

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

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

Use of androgenesis for estimating maternal and mitochondrial genome effects on development and oxygen consumption in rainbow trout, Oncorhynchus mykiss

Chromosome set manipulation was used to produce rainbow trout, Oncorhynchus mykiss, with identical nuclear backgrounds, but different maternal backgrounds to determine mitochondrial effects on development rate and oxygen consumption. Significant differences in development rate and oxygen consumption were observed between groups from different females. Development rates ranged from a mean of 317.97 degree days ( degrees d) to 335.25 degrees d in progeny from different females. Mean oxygen consumption rates ranged from 3.31 micromol O2 g(-1) wet mass h(-1) to 9.66 micromol O2 g(-1) wet mass h(-1). Oxygen consumption and development rate analysis revealed the two slowest developing groups had the highest oxygen consumption rates. Development rate differences between second generation clonal females indicate that mitochondrial genomes play a significant role on early development and are comparable to development rate differences between clonal lines of rainbow trout. These results indicate that selection for mitochondrial genomes could increase growth rates and possibly food conversion ratios in aquaculture species.

Polymorphism and differentiation of rainbow trout Y chromosomes

Most fish species show little morphological differentiation in the sex chromosomes. We have coupled molecular and cytogenetic analyses to characterize the male-determining region of the rainbow trout (Oncorhynchus mykiss) Y chromosome. Four genetically diverse male clonal lines of this species were used for genetic and physical mapping of regions in the vicinity of the sex locus. Five markers were genetically mapped to the Y chromosome in these male lines, indicating that the sex locus was located on the same linkage group in each of the lines. We also confirmed the presence of a Y chromosome morphological polymorphism among these lines, with the Y chromosomes from two of the lines having the more common heteromorphic Y chromosome and two of the lines having Y chromosomes morphologically similar to the X chromosome. The fluorescence in situ hybridization (FISH) pattern of two probes linked to sex suggested that the sex locus is physically located on the long arm of the Y chromosome. Fishes appear to be an excellent group of organisms for studying sex chromosome evolution and differentiation in vertebrates because they show considerable variability in the mechanisms and (or) patterns involved in sex determination.Key words: sex chromosomes, sex markers, cytogenetics, rainbow trout, fish.

Does differential selection on the 5S rDNA explain why the rainbow trout sex chromosome heteromorphism is not linked to the SEX locus?

Many but not all rainbow trout strains have morphologically distinguishable sex chromosomes. In these strains, the short arm of the X has multiple copies of 5S rDNA and a bright DAPI band near the centromere, both of which are missing from the Y chromosome, which has a very small short arm. We examined the presence of these markers using fluorescence in situ hybridization (FISH) in four different YY clonal lines derived from different strains and compared the results with sexed fish of the Donaldson strain with the normal X/Y heteromorphism. The Y chromosome in two of the YY clonal lines (Arlee and Swanson) is indistinguishable from the X chromosome and it is positive for 5S rDNA and the DAPI bright band. On the other hand, both 5S rDNA sequences and the DAPI band were not found on the Y chromosome in Hot Creek and Clearwater which have the normal Y. Thus the presence of these two cytogenetic markers may account for the size difference between the short arm of the X and Y chromosome found in most rainbow trout strains. In fishes the expression of one type of 5S rRNA is restricted to oocytes and previous work suggests that although XX males are fairly common, XY females are rare, implying a selective disadvantage for XY females. A hypothesis is presented to explain why this sex chromosome heteromorphism is not closely linked to the SEX locus, which is found on the long arm of the Y chromosome in rainbow trout.

Behavioral Differences Among Rainbow Trout Clonal Lines

Salmonid fish show considerable geographical variation in morphology, physiology, and behavior. Understanding the genetic mechanisms underlying this variation could be useful for enhancing aquaculture stocks, identifying unique populations for conservation, and determining the genetic factors underlying natural adaptation and domestication. As a first step toward the genetic dissection of salmonid behavioral diversity, we examined variation in behavior patterns among four clonal lines of rainbow trout ( Oncorhynchus mykiss ) derived from geographically diverse source populations with different domestication histories. Clonal lines were crossed with two outbred (i.e., not homozygous) females, and the resulting progenies were reared and tested under identical conditions. Clonal line had significant genetic effects on mean swim level, hiding, foraging, startle response, and aggression level. Multiple comparisons suggest that domestication history of the source populations had a strong influence on these behavior patterns. Progeny of two clonal lines derived from populations reared in captivity for over 100 years exhibited reductions in predator avoidance behavior patterns and increases in aggression compared to progeny of two clonal lines from more recently domesticated populations. These results will facilitate future investigation of the genetic factors underlying population variation in these behavior patterns influenced by domestication.

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

We report the identification of a single major chromosomal region controlling natural killer (NK) cell-like activity in rainbow trout (Oncorhynchus mykiss). A genetic map based on 484 AFLP and 39 microsatellite genotypes from 106 doubled haploid fish was constructed. These fish were produced by androgenesis from a hybrid of two clonal lines divergent in NK-like activity. NK-like activities for 75 of the doubled haploids were quantified by an in vitro chromium release assay utilizing (51)Cr-labeled YAC-1 target cells. Composite interval mapping revealed a single major quantitative trait locus (QTL) associated with NK-like activity in this rainbow trout model. Genetic mapping revealed this QTL to also be unlinked to: fragmented MHC class I and MHC class II regions, the leukocyte receptor cluster, the natural killer cell enhancement factor ( NKEF) gene, the RAG-1 gene, and two QTL associated with resistance to infectious pancreatic necrosis virus in rainbow trout. Collectively, these results extend the utility of rainbow trout as an immunological model and are consistent with the idea that a single chromosomal region homologous to the natural killer cell complex (NKC) located on syntenic portions of mouse chromosome (Chr) 6, human Chr 12, and rat Chr 4 may exist in a lower vertebrate model.