Search

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.

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.

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.

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

Domesticated populations of rainbow trout (Oncorhynchus mykiss) probably have been selected for high growth rates and large body size, which may have resulted in reduced sprint swimming performance. We hypothesized that more domesticated populations of rainbow trout would have higher growth rates and larger body size, but slower swim speed, relative to semiwild populations. We tested this hypothesis by measuring body mass, body length, and sprint swim speed multiple times over 92 days in progeny from crosses between males from three clonal lines and an outbred female. We found significantly higher body masses and significantly slower swim speeds in the highly domesticated Arlee and Hot Creek progeny groups compared with the semiwild Swanson hybrid progeny group, supporting our hypothesis. Growth rates also differed significantly among groups, but at the ages measured, the Swanson hybrid progeny had an intermediate growth rate. However, given the differences in body mass, either growth rates were higher in the more domesticated progeny groups at young unmeasured ages and (or) they hatched at a significantly larger body mass than the semiwild Swanson hybrid. Les populations domestiquees de truites arc-en-ciel (Oncorhynchus mykiss) ont probablement ete selectionnees pour leur fort taux de croissance et leur importante taille corporelle, ce qui peut avoir conduit a une performance reduite de la nage de pointe. Nous posons l'hypothese selon laquelle les populations plus domestiquees de truites arc-en-ciel devraient avoir un taux de croissance plus eleve et une taille corporelle plus forte, mais une vitesse de nage plus lente, par comparaison aux populations semi-sauvages. Nous avons teste cette hypothese en mesurant la masse corporelle, la longueur du corps et la vitesse de nage de pointe a de nombreuses reprises sur une periode de 92 jours chez des rejetons de croisements entre des males de trois lignees clonales et une femelle exogame. La masse corporelle est significativement plus forte et la vitesse de nage significativement plus lente chez les groupes de rejetons des lignees fortement domestiquees Arlee et Hot Creek, par rapport au groupe de rejetons hybrides de la lignee semi-sauvage Swanson, ce qui confirme notre hypothese. Les taux de croissance different aussi parmi ces groupes, mais aux ages ou nous avons fait les mesures, les rejetons hybrides Swanson ont un taux de croissance intermediaire. Cependant, compte tenu des differences de masse corporelle, il se peut que les taux de croissance aient ete plus eleves dans les groupes de rejetons plus domestiques aux ages anterieurs a nos mesures ou alors (et) que ces rejetons aient eclos avec une masse corporelle significativement plus grande que les hybrides semi-sauvages Swanson. [Traduit par la Redaction], Introduction In natural populations of animals, strong correlations between large body size and increased evolutionary fitness are frequently measured; large body size is related to increased fecundity, reduced predation, and, [...]

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.

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 F2 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 F2 animals, two single nucleotide polymorphisms (SNPs) in the promoter region were significantly associated with intramuscular fat (P

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.