Roger L. Vallejo

Single nucleotide polymorphism discovery in rainbow trout by deep sequencing of a reduced representation library.

BackgroundTo enhance capabilities for genomic analyses in rainbow trout, such as genomic selection, a large suite of polymorphic markers that are amenable to high-throughput genotyping protocols must be identified. Expressed Sequence Tags (ESTs) have been used for single nucleotide polymorphism (SNP) discovery in salmonids. In those strategies, the salmonid semi-tetraploid genomes often led to assemblies of paralogous sequences and therefore resulted in a high rate of false positive SNP identification. Sequencing genomic DNA using primers identified from ESTs proved to be an effective but time consuming methodology of SNP identification in rainbow trout, therefore not suitable for high throughput SNP discovery. In this study, we employed a high-throughput strategy that used pyrosequencing technology to generate data from a reduced representation library constructed with genomic DNA pooled from 96 unrelated rainbow trout that represent the National Center for Cool and Cold Water Aquaculture (NCCCWA) broodstock population.ResultsThe reduced representation library consisted of 440 bp fragments resulting from complete digestion with the restriction enzyme Hae III; sequencing produced 2,000,000 reads providing an average 6 fold coverage of the estimated 150,000 unique genomic restriction fragments (300,000 fragment ends). Three independent data analyses identified 22,022 to 47,128 putative SNPs on 13,140 to 24,627 independent contigs. A set of 384 putative SNPs, randomly selected from the sets produced by the three analyses were genotyped on individual fish to determine the validation rate of putative SNPs among analyses, distinguish apparent SNPs that actually represent paralogous loci in the tetraploid genome, examine Mendelian segregation, and place the validated SNPs on the rainbow trout linkage map. Approximately 48% (183) of the putative SNPs were validated; 167 markers were successfully incorporated into the rainbow trout linkage map. In addition, 2% of the sequences from the validated markers were associated with rainbow trout transcripts.ConclusionThe use of reduced representation libraries and pyrosequencing technology proved to be an effective strategy for the discovery of a high number of putative SNPs in rainbow trout; however, modifications to the technique to decrease the false discovery rate resulting from the evolutionary recent genome duplication would be desirable.

A second generation genetic map for rainbow trout (Oncorhynchus mykiss)

BackgroundGenetic maps characterizing the inheritance patterns of traits and markers have been developed for a wide range of species and used to study questions in biomedicine, agriculture, ecology and evolutionary biology. The status of rainbow trout genetic maps has progressed significantly over the last decade due to interest in this species in aquaculture and sport fisheries, and as a model research organism for studies related to carcinogenesis, toxicology, comparative immunology, disease ecology, physiology and nutrition. We constructed a second generation genetic map for rainbow trout using microsatellite markers to facilitate the identification of quantitative trait loci for traits affecting aquaculture production efficiency and the extraction of comparative information from the genome sequences of model fish species.ResultsA genetic map ordering 1124 microsatellite loci spanning a sex-averaged distance of 2927.10 cM (Kosambi) and having 2.6 cM resolution was constructed by genotyping 10 parents and 150 offspring from the National Center for Cool and Cold Water Aquaculture (NCCCWA) reference family mapping panel. Microsatellite markers, representing pairs of loci resulting from an evolutionarily recent whole genome duplication event, identified 180 duplicated regions within the rainbow trout genome. Microsatellites associated with genes through expressed sequence tags or bacterial artificial chromosomes produced comparative assignments with tetraodon, zebrafish, fugu, and medaka resulting in assignments of homology for 199 loci.ConclusionThe second generation NCCCWA genetic map provides an increased microsatellite marker density and quantifies differences in recombination rate between the sexes in outbred populations. It has the potential to integrate with cytogenetic and other physical maps, identifying paralogous regions of the rainbow trout genome arising from the evolutionarily recent genome duplication event, and anchoring a comparative map with the zebrafish, medaka, tetraodon, and fugu genomes. This resource will facilitate the identification of genes affecting traits of interest through fine mapping and positional cloning of candidate genes.

RNA-Seq Identifies SNP Markers for Growth Traits in Rainbow Trout

Fast growth is an important and highly desired trait, which affects the profitability of food animal production, with feed costs accounting for the largest proportion of production costs. Traditional phenotype-based selection is typically used to select for growth traits; however, genetic improvement is slow over generations. Single nucleotide polymorphisms (SNPs) explain 90% of the genetic differences between individuals; therefore, they are most suitable for genetic evaluation and strategies that employ molecular genetics for selective breeding. SNPs found within or near a coding sequence are of particular interest because they are more likely to alter the biological function of a protein. We aimed to use SNPs to identify markers and genes associated with genetic variation in growth. RNA-Seq whole-transcriptome analysis of pooled cDNA samples from a population of rainbow trout selected for improved growth versus unselected genetic cohorts (10 fish from 1 full-sib family each) identified SNP markers associated with growth-rate. The allelic imbalances (the ratio between the allele frequencies of the fast growing sample and that of the slow growing sample) were considered at scores >5.0 as an amplification and <0.2 as loss of heterozygosity. A subset of SNPs (n = 54) were validated and evaluated for association with growth traits in 778 individuals of a three-generation parent/offspring panel representing 40 families. Twenty-two SNP markers and one mitochondrial haplotype were significantly associated with growth traits. Polymorphism of 48 of the markers was confirmed in other commercially important aquaculture stocks. Many markers were clustered into genes of metabolic energy production pathways and are suitable candidates for genetic selection. The study demonstrates that RNA-Seq at low sequence coverage of divergent populations is a fast and effective means of identifying SNPs, with allelic imbalances between phenotypes. This technique is suitable for marker development in non-model species lacking complete and well-annotated genome reference sequences.

Rainbow trout resistance to bacterial cold-water disease is moderately heritable and is not adversely correlated with growth.

The objectives of this study were to estimate the heritabilities for and genetic correlations among resistance to bacterial cold-water disease and growth traits in a population of rainbow trout (Oncorhynchus mykiss). Bacterial cold-water disease, a chronic disease of rainbow trout, is caused by Flavobacterium psychrophilum. This bacterium also causes acute losses in young fish, known as rainbow trout fry syndrome. Selective breeding for increased disease resistance is a promising strategy that has not been widely used in aquaculture. At the same time, improving growth performance is critical for efficient production. At the National Center for Cool and Cold Water Aquaculture, reducing the negative impact of diseases on rainbow trout culture and improving growth performance are primary objectives. In 2005, when fish averaged 2.4 g, 71 full-sib families were challenged with F. psychrophilum and evaluated for 21 d. Overall survival was 29.3% and family rates of survival varied from 1.5 to 72.5%. Heritability of postchallenge survival, an indicator of disease resistance, was estimated to be 0.35 +/- 0.09. Body weights at 9 and 12 mo posthatch and growth rate from 9 to 12 mo were evaluated on siblings of the fish in the disease challenge study. Growth traits were moderately heritable, from 0.32 for growth rate to 0.61 for 12-mo BW. Genetic and phenotypic correlations between growth traits and resistance to bacterial cold-water disease were not different from zero. These results suggest that genetic improvement can be made simultaneously for growth and bacterial cold-water disease resistance in rainbow trout by using selective breeding.

Response to selection for bacterial cold water disease resistance in rainbow trout

A family-based selection program was initiated at the National Center for Cool and Cold Water Aquaculture in 2005 to improve resistance to bacterial cold water disease (BCWD) in rainbow trout. The objective of this study was to estimate response to 2 generations of selection. A total of 14,841 juvenile fish (BW = 3.1 g; SD = 1.1 g) from 230 full-sib families and 3 randomly mated control lines were challenged intraperitoneally with Flavobacterium psychrophilum, the bacterium that causes BCWD, and mortalities were observed for 21 d. Selection was applied to family EBV derived from a proportional-hazards frailty (animal) model while constraining rate of inbreeding to <or=1% per generation. After adjusting for nongenetic effects, survival rate of select-line families increased by 24.6 +/- 6.8 and 44.7 +/- 6.7 (cumulative) percentage points after 1 and 2 generations of selection, respectively (P < 0.01). Genetic trend, estimated from a linear animal model that fit genetic group effects, was 19.0 +/- 4.1 percentage points per generation and approached significance (P = 0.07). Heritability estimates from the proportional-hazards frailty model and linear animal model were similar (0.22 and 0.23, respectively), and family EBV from both models were highly correlated (-0.92). Accuracy of selection, estimated as the correlation between mid-parent EBV and progeny survival rate, was 0.20 (P < 0.01) for the proportional-hazards frailty model and 0.18 (P = 0.01) for the linear animal model. Accuracy estimates were not different (P = 0.81) between the models. This study demonstrates that selective breeding can be effective for improving resistance to experimental BCWD challenge in rainbow trout.

A first generation integrated map of the rainbow trout genome.

BackgroundRainbow 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. An integrated physical and genetic map is needed to facilitate fine mapping of QTL and the selection of positional candidate genes for incorporation in marker-assisted selection (MAS) programs for improving rainbow trout aquaculture production.ResultsThe first generation integrated map of the rainbow trout genome is composed of 238 BAC contigs anchored to chromosomes of the genetic map. It covers more than 10% of the genome across segments from all 29 chromosomes. Anchoring of 203 contigs to chromosomes of the National Center for Cool and Cold Water Aquaculture (NCCCWA) genetic map was achieved through mapping of 288 genetic markers derived from BAC end sequences (BES), screening of the BAC library with previously mapped markers and matching of SNPs with BES reads. In addition, 35 contigs were anchored to linkage groups of the INRA (French National Institute of Agricultural Research) genetic map through markers that were not informative for linkage analysis in the NCCCWA mapping panel. The ratio of physical to genetic linkage distances varied substantially among chromosomes and BAC contigs with an average of 3,033 Kb/cM.ConclusionsThe integrated map described here provides a framework for a robust composite genome map for rainbow trout. This resource is needed for genomic analyses in this research model and economically important species and will facilitate comparative genome mapping with other salmonids and with model fish species. This resource will also facilitate efforts to assemble a whole-genome reference sequence for rainbow trout.

Suggestive association of major histocompatibility IB genetic markers with resistance to bacterial cold water disease in rainbow trout (Oncorhynchus mykiss).

Genes within the major histocompatibility complex (MHC) are important for both innate and adaptive immune responses in mammals; however, much less is known regarding their contribution in teleost fishes. We examined the involvement of four major histocompatibility (MH) genomic regions in rainbow trout in resistance to the causative agent of bacterial coldwater disease (BCWD), Flavobacterium psychrophilum. Fish from the 2005 NCCCWA brood-year (71 full-sib families) were challenged with F. psychrophilum strain CSF 259–93. The overall mortality rate was 70%, with large variation in mortality between families. Disease resistance was quantified as post-challenge days to death. Phenotypic variation and additive genetic variation were estimated using mixed models of survival analysis. To examine association, eight microsatellite markers were isolated from MH gene-containing BAC clones and mapped onto the rainbow trout genetic linkage map. The parents and grandparents of the 2005 brood-year families were genotyped with these eight markers and another two markers tightly linked to the MH-IB region to assess the extent of linkage disequilibrium (LD) of MH genomic regions MH-IA, MH-IB, TAP1, and MH-II with survival post-challenge. MH-IB and MH-II markers were linked to BCWD survivability when data were analyzed by family. Tests for disease association at the population level substantiated the involvement of MH-IB, but not MH-II, with disease resistance. The impact of selective breeding for disease resistance on MH sequence variation is discussed in the context of aquaculture production.

Detection and Validation of QTL Affecting Bacterial Cold Water Disease Resistance in Rainbow Trout Using Restriction-Site Associated DNA Sequencing

Bacterial cold water disease (BCWD) causes significant economic loss in salmonid aquaculture. Using microsatellite markers in a genome scan, we previously detected significant and suggestive QTL affecting phenotypic variation in survival following challenge with Flavobacterium psychrophilum, the causative agent of BCWD in rainbow trout. In this study, we performed selective genotyping of SNPs from restriction-site associated DNA (RAD) sequence data from two pedigreed families (2009070 and 2009196) to validate the major QTL from the previous work and to detect new QTL. The use of RAD SNPs in the genome scans increased the number of mapped markers from ~300 to ~5,000 per family. The significant QTL detected in the microsatellites scan on chromosome Omy8 in family 2009070 was validated explaining up to 58% of the phenotypic variance in that family, and in addition, a second QTL was also detected on Omy8. Two novel QTL on Omy11 and 14 were also detected, and the previously suggestive QTL on Omy1, 7 and 25 were also validated in family 2009070. In family 2009196, the microsatellite significant QTL on Omy6 and 12 were validated and a new QTL on Omy8 was detected, but none of the previously detected suggestive QTL were validated. The two Omy8 QTL from family 2009070 and the Omy12 QTL from family 2009196 were found to be co-localized with handling and confinement stress response QTL that our group has previously identified in a separate pedigreed family. With the currently available data we cannot determine if the co-localized QTL are the result of genes with pleiotropic effects or a mere physical proximity on the same chromosome segment. The genetic markers linked to BCWD resistance QTL were used to query the scaffolds of the rainbow trout reference genome assembly and the QTL-positive scaffold sequences were found to include 100 positional candidate genes. Several of the candidate genes located on or near the two Omy8 QTL detected in family 2009070 suggest potential linkages between stress response and the regulation of immune response in rainbow trout.

Temporal and pathogen-load dependent changes in rainbow trout (Oncorhynchus mykiss) immune response traits following challenge with biotype 2 Yersinia ruckeri

Rainbow trout infected with Yersinia ruckeri, the causative agent of enteric redmouth disease (ERM), produce a pro-inflammatory and acute-phase response attributed in part to the innate recognition of bacterial-produced flagellin. Recently, variants of Y. ruckeri have been identified that lack flagella and associated motility. These strains are classified as biotype 2 (BT2) Y. ruckeri and these are considered an emerging problem in salmonid aquaculture. Little is known about the salmonid immune response to these variants. Herein, we report temporal and quantitative changes in rainbow trout immune response parameters following a primary challenge with BT2 Y. ruckeri strain YRNC10. Fish were injection-challenged with ten-fold dilutions of viable bacteria and sampled on days 1, 3, 5 and 7 post-challenge. TNFalpha1 and IL1-beta1 transcripts were increased by day 1 post-challenge, and on days 3, 5 and 7 maximal gene transcript up-regulation occurred at a threshold of approximately 64-256CFU per mg spleen tissue. Infection induced robust SAA gene up-regulation that was significantly correlated with increased gene expression of IL-1beta1 (r=0.81, P<0.0001) and TNFalpha1 (r=0.55, P<0.0001). Y. ruckeri infection induced modest changes in INFgamma and Mx-1 gene transcript abundance at intermediate or high challenge doses and the expression patterns of both genes were positively correlated with pro-inflammatory gene and acute-phase gene transcription patterns. TNF superfamily 13b (BAFF) gene expression was significantly down-regulated in response to infection on days 3, 5 and 7 at the highest challenge doses. The spleen somatic index was significantly increased on days 3, 5 and 7 post-infection and positively correlated with spleen colony forming units and abundance of gene transcripts SAA, TNFalpha1, and IL1-beta1. In summary, rainbow trout had a strong innate response following challenge with BT2 Y. ruckeri strain YRNC10 indicating that flagellin expression is not required for production of a robust pro-inflammatory and acute-phase gene transcription response. This study further supports the use of SAA transcript abundance and spleen somatic index as general measures of immunological status and fish health.

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

BackgroundRainbow 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.ResultsThe 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.ConclusionThe 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.