Email updates

Keep up to date with the latest news and content from BMC Genomics and BioMed Central.

Open Access Research article

Development and characterisation of an expressed sequence tags (EST)-derived single nucleotide polymorphisms (SNPs) resource in rainbow trout

Mekki Boussaha1*, René Guyomard1, Cédric Cabau2, Diane Esquerré3 and Edwige Quillet1

Author affiliations

1 INRA, UMR 1313 Génétique Animale et Biologie Intégrative, 78350, Jouy-en-Josas, France

2 INRA, SIGENAE UR83 Recherches Avicoles, 37380, Nouzilly, France

3 INRA, UMR 444 Laboratoire de Génétique Cellulaire Plateforme GET, Castanet, Tolosan, France

For all author emails, please log on.

Citation and License

BMC Genomics 2012, 13:238  doi:10.1186/1471-2164-13-238


The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2164/13/238


Received:20 February 2012
Accepted:13 June 2012
Published:13 June 2012

© 2012 Boussaha et al.; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License ( http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

There is considerable interest in developing high-throughput genotyping with single nucleotide polymorphisms (SNPs) for the identification of genes affecting important ecological or economical traits. SNPs are evenly distributed throughout the genome and are likely to be functionally relevant. In rainbow trout, in silico screening of EST databases represents an attractive approach for de novo SNP identification. Nevertheless, EST sequencing errors and assembly of EST paralogous sequences can lead to the identification of false positive SNPs which renders the reliability of EST-derived SNPs relatively low. Further validation of EST-derived SNPs is therefore required. The objective of this work was to assess the quality of and to validate a large number of rainbow trout EST-derived SNPs.

Results

A panel of 1,152 EST-derived SNPs was selected from the INRA Sigenae SNP database and was genotyped in standard and double haploid individuals from several populations using the Illumina GoldenGate BeadXpress assay. High-quality genotyping data were obtained for 958 SNPs representing a genotyping success rate of 83.2 %, out of which, 350 SNPs (36.5 %) were polymorphic in at least one population and were designated as true SNPs. They also proved to be a potential tool to investigate genetic diversity of the species, as the set of SNP successfully sorted individuals into three main groups using STRUCTURE software. Functional annotations revealed 28 non-synonymous SNPs, out of which four substitutions were predicted to affect protein functions. A subset of 223 true SNPs were polymorphic in the two INRA mapping reference families and were integrated into the INRA microsatellite-based linkage map.

Conclusions

Our results represent the first study of EST-derived SNPs validation in rainbow trout, a species whose genome sequences is not yet available. We designed several specific filters in order to improve the genotyping yield. Nevertheless, our selection criteria should be further improved in order to reduce the observed high rate of false positive SNPs which results from the occurrence of whole genome duplications.

Background

International genome initiatives have resulted in draft sequences of the genome of several farm animals (cattle, pig, chicken, and horse) and of model fish species (zebrafish (Danio rerio), medaka (Oryzias latipes), stickleback (Gasterosteus aculeatus), takifugu (Takifugu rubripes), and tetraodon (Tetraodon nigroviridis)). Whole genome sequencing are currently underway for a number of aquaculture species: rainbow trout (Oncorhynchus mykiss), Atlantic salmon (Salmo salar), Nile tilapia (Oreochromis niloticus), Asian seabass (Lates calcarifer), European seabass (Dicentrarchus labrax), channel catfish (Ictalurus punctatus) and common carp (Cyprinus carpio). At the same time, high-throughput genomic tools have been developed, improving the description of genomic structure and function.

Projects associated with genome sequencing activities using different breeds from the same species have provided the opportunity to discover hundreds of thousands of potential single-base changes, also known as single nucleotide polymorphisms (SNPs) or short insertion/deletion mutations (indels). The bi-allelic nature of SNPs makes them less informative than microsatellites. Nevertheless, SNPs are considered as a highly reliable and valuable molecular marker system for genotyping and selective breeding because of their omnipresence throughout the entire genome, both within gene coding and non-coding regions.

SNPs in gene coding sequences can be either synonymous (silent polymorphism) or non-synonymous (replacement polymorphism). They are of particular interest to study the genetics of expressed genes and to map functional traits. Synonymous SNPs may alter RNA secondary structures and can affect protein conformation and function [1]. Non-synonymous SNPs can potentially have deleterious functional effects because they lead to changes in amino acid sequences and possibly affect protein structure and function [2,3].

SNPs in non-coding regions can occur in introns, promoters, intergenic sequences, and in 5'- or 3'-untranslated regions. They may alter gene expression by affecting gene splicing, transcription factor binding, mRNA degradation, or non-coding RNA sequences.

Over the last decades, large-scale SNP production initiatives have been associated with the development of high-throughput genotyping technologies that facilitate the simultaneous analysis of hundreds of thousands of SNPs. These low-cost but highly reliable assays have permitted fine-scale gene mapping and candidate gene association studies for complex traits in several species such as humans [4], mouse [5] chicken [6], cattle [7] and sheep [8].

In species whose complete genome sequences are not yet accessible, the increasing availability of expressed sequence tags (ESTs) represents an alternative in silico strategy for de novo SNP identification. This approach does not require any additional bench work, offers a low cost source of SNPs, and has been recently used in a few aquaculture species such as blue and channel catfish species [9], and salmonids [10-13]. Moreover, EST-derived SNPs are considered as gene-derived SNPs since they are located within gene coding and 3′-UTR regions and they can lead to the identification of quantitative trait nucleotides (QTN) [14].

However, the usefulness of EST-derived SNPs remains putative until their true informativity (sequence polymorphism) and duplication status have been checked with genomic DNA in the populations of interest. Although it is possible to use base quality values to discern true allelic variations from sequencing errors, validation is a key step for detection of true SNPs [15]. This is generally carried out by genotyping several population samples with a subset of the EST-derived SNPs [10].

Rainbow trout is the most widely cultivated cold freshwater fish in the world. It has great potential for aquaculture and recreational sport fisheries. In addition to its commercial interest, rainbow trout is also a model species for a wide range of genome-related research activities [16].

The rainbow trout haploid genome size was estimated to be between 2.4 and 3.0 × 109 bp [17,18]. A common ancestor of rainbow trout and other salmonids has undergone a fourth whole-genome duplication (4R WGD) event about 25 to 100 million years ago, which was followed by a period of re-diploidization resulting in a semi-tetraploid state [19]. It has been estimated that up to half of the loci are still duplicated [20]. Although the tetraploidization event increases the genome complexity, it also makes the salmonids an attractive model to study the mechanisms behind the whole-genome duplication event and the subsequent reduction of one of the two copies of the duplicated gene(s).

Both the interest brought into rainbow trout as a research model and the need for its genetic improvement for aquaculture production efficiency and product quality led to the development of several genomic resources for this species. Meanwhile, great efforts have been and are still devoted to the development of SNP genetic markers [21-23].

Previous efforts using reduced representation libraries [22] and reference transcriptome datasets [24] resulted in the production of up to 47,000 and 58,000 putative SNPs, respectively. A subset of 384 randomly selected SNPs were genotyped on individual fish and 184 (48 %) were validated [22]. The observed low validation rate could be partly explained by the presence of paralogous sequences with allelic variation which resulted in the production of false positive SNPs.

Finally, these putative SNPs were not yet publicly available. Therefore, EST-derived SNPs could represent an alternative and complementary in silico approach to assess the quality and to validate larger numbers of SNPs. These resources will add to the already available 184 SNPs validated from the reduced representation libraries study.

Miller and co-authors [23] have also used the RAD (Restriction site Associated DNA) sequencing technology for low density SNP genotyping and reported the construction of a high-resolution linkage map containing 4,563 markers. However, the flanking sequences for these SNPs were only 68 nucleotides long and thus may not be suitable for the design of high-throughput genotyping assays, such as the Illumina assays. Retrieving longer flanking sequences suitable for high-throughput genotyping studies using these RAD-associated markers will need additional information on the whole genome sequence. Efforts are in progress in France and USA [25,26] to provide a rainbow trout reference genome sequence in the near future. Nevertheless, in both cases, aiming at facilitating the assembly step, the sequencing was performed using a doubled haploid homozygous DNA sample which hinders the identification of new SNPs.

Mining EST datasets remains an attractive alternative approach for in silico SNP identification in rainbow trout. Up to 31,121 in silico EST-derived SNPs are currently available at the INRA Sigenae database ( http://www.sigenae.org/ webcite). However, they do not provide any information neither on their true informativity nor on their duplication status. Therefore, it is necessary to validate the status of these markers. Validation of rainbow trout EST-derived SNPs in a large number of populations will not only allow to identify fully informative true SNPs but also will highlight the proportion of informative SNPs shared across different populations, a crucial information to efficiently design future rainbow trout specific SNP chips. These new tools will contribute to studies on population genetics and will facilitate quantitative trait loci (QTL) identification, and marker assisted selection.

In the present study, a panel of 1,152 EST-derived SNPs was selected from the Sigenae SNP database and were subsequently assayed for allelic variation in several rainbow trout population samples using the Illumina GoldenGate assays. Successfully validated EST-derived SNPs were used to analyse the genetic diversity in three bisexually reproducing experimental stocks and a collection of doubled haploid (DH) clones and to update the INRA linkage map by integrating 223 new markers.

Methods

Selection of SNPs for validation

The INRA Sigenae rainbow trout EST-derived SNP database ( http://www.sigenae.org/ webcite; restricted access) was used to select a validation SNP panel. A public version of this release will be available in the near future. Almost 31,121 SNPs were produced by assembling EST sequences. Briefly, several stringent filters were used to improve the quality of predicted SNPs: (1) the value of the local depth at the polymorphic position must be at least equal to 7; (2) the 4 bases flanking regions around the SNP position need to be exactly conserved within the aligned sequences; (3) the minimal number of sequences with the lowest represented base must be at least equal to 3; (4) gaps on consensus sequences were ignored; and (5) N or gaps on sequences were ignored.

Several selection filters (Figure 1) were applied in order to select a panel of 1,152 EST-derived SNPs for validation: (1) in order to meet the requirements for probe design constraints for the Illumina genotyping platform, all SNPs with less than 60 nucleotides between two neighbouring SNPs and with flanking sequences less than 100 nucleotides long were removed; (2) in order to overcome problems due to exon-intron junctions, the SNP flanking sequences were aligned against rainbow trout BAC-end sequences [27] using megablast tools and against zebrafish, medaka, and stickleback genomic sequences using blastn tools. All SNP sequences with an alignment length equal to the flanking sequence length were selected for further analysis. The filtered SNP sequences were then submitted to Illumina to assess their design quality. Only those showing a minimum quality score of 0.6 were further filtered against sequence similarities between each other and against the presence of repetitive sequences. After applying the above filters, a panel of 1,152 EST-derived SNPs was constructed and was used to genotype a large number of rainbow trout individuals.

thumbnailFigure 1. Selection of the validation panel. Filters used to select EST-derived SNPs for validation from the INRA Sigenae SNP database release som10 were summarized.

DNA sources

Two hundred and fifty seven DNA individual rainbow trout were genotyped for each of the 1,152 SNPs using the Illumina GoldenGate assay. These include 37 INRA doubled haploid (DH) individuals, DNA from 10 DH individuals from various origins provided by Dr Gary Thoorgard (ARS), 20 individuals from the INRA synthetic reference strain (INRA-SY), 20 individuals from the INRA spring spawning strain (INRA-SP), and DNA from 44 individuals from five NCCCWA mapping families [28] provided by Dr Yniv Palti (ARS). The two INRA reference mapping families (two parents with four grandparents and 120 DH progeny) were also included. DNA was isolated from fin clips stored in 95 % ethanol, according to the protocols previously described [29].

INRA rainbow trout fin clips were collected from euthanized and/or anesthetized fish elevated at the INRA fish farm facilities. Under French regulation, the INRA facilities are authorized for experimental activities and both the staff of the facilities and scientists have personal authorization to conduct animal experimentations. All animal manipulations were done according to the good animal practice as defined by the French Direction of Veterinary Services.

Genotyping

High-throughput genotyping reactions were performed at the INRA genomics GET PlaGe core facility, using the Illumina GoldenGate BeadXpress systems, according to the manufacturer's protocol [30]. SNPs with an Illumina design score above 0.6 were retained for further analysis. Oligonucleotides were designed, synthesised, and assembled into three custom oligo pooled assays (OPA) by Illumina Inc.

Genotype clustering was performed using the GenomeStudio software (Illumina Inc.). GenCall and GenTrain quality scores for each genotype were generated. A GenCall score cutoff of 0.25 was used to determine valid genotypes at each SNP and the retained SNPs had to have a minimum GenTrain score of 0.25 (a stringent criterion that is used in human genetic studies) [31]. Clusters were visually inspected to ensure high quality data. Genotype calls were exported as spreadsheets from the GenomeStudio software for further analysis.

Population structure

The STRUCTURE software [32] was used to assess the population structure. This program implements a model-based clustering method to infer population structure using genotype data of unlinked markers. We used the admixture model and correlated allele frequency version of the STRUCTURE program [33]. To choose the most likely number of clusters modelling the data, several analyses were performed, for a number of fixed subgroups K (number of populations) from 1 to 5. Each analysis involved five independent runs with a burn-in period of 50,000 and 200,000 iterations for the likelihood estimation. The best K value which corresponds to the K with the highest Delta K score was determined using a non parametric test as previously described [34]. This test uses an ad hoc quantity (delta K) calculated based on the second order rate of change of the likelihood (delta K).

Functional annotations of polymorphic SNPs

Both contig and SNP allele sequences were analysed for gene content by blastx using the ENSEMBL non redundant protein databases for zebrafish (Danio_rerio.Zv9.64.pep.all.fa).

Blastx searches were carried out using an e-value cut off of 1e-5. The blastx search results were filtered to remove non specific homologies using the following filtration steps: (1) the Ensembl protein ID in the blastx results were renamed by their corresponding Ensembl gene ID (since each gene may encode several peptides due to alternative splicing), (2) for each sequence read (query ID) with a gene hit (subject ID), results were filtered to keep only the hits with the minimal e-value score; and (3) sequence reads with several hits having the same minimal e-value were further filtered to keep the hits with the highest HSP (high-scoring segment pairs; calculated as the product of % identity multiplied by alignment length). Only SNP sequences and their corresponding contig sequences having a gene hit were used for further analysis.

For each contig read, query start and query end positions were used to retrieve corresponding contig sequences between these two values. DNA sequences were then translated and the translation product was used to construct a specific RT peptide database. Both SNP allele sequences were then analyzed for synonymous/non-synonymous SNPs by blastx using the produced RT peptide database. For synonymous SNPs, both allele sequences should result in a perfect match with a given peptide sequence (100 % identity). For non-synonymous SNPs, one allele sequence should result in a perfect match and the other should present only one amino acid mismatch.

Finally, we assessed the deleterious effect of non-synonymous SNPs using SIFT (Sorting Intolerant From Tolerant) program (sift.jcvi.org/). Prediction was carried out using the SIFT sequence tool through PSI-Blast searches against UniProt - SwissProt databases (release 57.15, April 2011). Median conservation of sequences was fixed to 3.0 and hits showing more than 90 % identity to the query sequence were removed.

Linkage map construction

Linkage groups were constructed with CARTHAGENE [35] and optimized with the annealing option (argument values: 15, 300, 0.1, 0.5) (see Carthagene help for argument meaning). Since interference is close to 1 in salmonids, we used the percentage of recombination as mapping function. Graphical representations were obtained with MAPCHART [36].

Results and Discussions

Sigenae SNP database characterization and selection of a subset for validation

The Sigenae rainbow trout EST-derived SNP database ( http://www.sigenae.org webcite) contains 31,121 putative SNPs identified in 13,374 EST contigs (Table 1). The total length of contig sequences was estimated to be 2.23 Mb with an average contig length of 1,889 bp ranging from 134 to 9,913 bp. This corresponds to one SNP every 716 bp which is slightly higher than previously reported frequencies from a panel of SNPs obtained using the RAD sequencing approach [23]. The average sequence coverage was estimated to be 12.7 sequence/contig and ranged from 7 to 466 sequences/contig.

Table 1. Distribution of SNPs in EST contigs

Total number of EST contigs containing one or more SNPs were indicated for both the initial database and the validation panel.

Almost 83 % of the EST-derived SNPs were identified from contigs containing one to five SNPs (Table 1). The mean minor sequence frequency among all SNPs was 0.37 ± 0.1 (SD), while the mean observed heterozygosity based on sequence coverage at the polymorphic site was 0.45 ± 0.07, and the mean PIC (polymorphism information content) was 0.34 ± 0.04 (Additional file 1: Sheet 1).

Additional file 1. Minor sequence frequencies, heterozygosity and PIC values were indicated for all SNPs in the INRA Sigenae database (worksheet 1) and for the validation panel (worksheet 2). Allele frequencies were calculated by direct count of sequence depth for each allele, expected heterozygosities per locus (He) were estimated by He = 2pq and PIC values were estimated by PIC = H-2p2q2, where p and q are the sequence frequencies of the two alternate alleles.

Format: XLS Size: 8.9MB Download file

This file can be viewed with: Microsoft Excel ViewerOpen Data

After application of several selective filters (Figure 1) designed to improve the expected yield of genotyping, a subset of 1,152 SNPs was selected for our study (Additional file 2). Three OPA (Oligo Pool Assays) each comprising 384 SNPs were designed and were called the validation panel. SNPs from the validation panel were identified in 863 contigs, of which 66 % contain one to five SNPs (Table 1). The mean minor sequence frequency among the validation panel was estimated to 0.36 ± 0.1 (SD), while the mean observed heterozygosity based on sequence coverage at the polymorphic site was 0.44 ± 0.08 and the mean PIC was 0.34 ± 0.05 (Additional file 1: Sheet 2) which were very close to those calculated from the Sigenae SNPs database.

Additional file 2. Illumina quality design scores for the validation panel.

Format: XLS Size: 466KB Download file

This file can be viewed with: Microsoft Excel ViewerOpen Data

SNP validation

The efficiency of the selection approach and the relevance of the resulting SNPs were assessed by genotyping the validation panel in a number of rainbow trout doubled haploid (DH) and standard individuals from three different domestic populations.

Assays were developed for 1,152 putative EST-derived SNPs, out of which 958 (83 %) were successfully genotyped (Table 2 and Additional file 3: Sheets 1–5) while genotyping failed for 194 SNPs (17 %). These did not either cluster well according to genotype or failed to amplify most probably because of the sequence complexity or the presence of polymorphisms within flanking sequences or failed manufacture with Illumina. These were considered "failed assays". Out of the 958 successfully genotyped SNPs, 55 % were selected from contigs containing no more than 5 SNPs and the overall proportion of successfully genotyped SNPs over those from the validated panel did not depend on the SNP content in EST contigs (Table 3).

Additional file 3. SNP allele frequencies. Allele frequencies for both alleles and minor allele frequencies (MAF) were shown for DH clonal lines (worksheet 1), INRA-SY (worksheet 2), INRA-SP (worksheet 3), NCCCWA (worksheet 4), all 3 populations (worksheet 5) and for true SNPs (worksheet 6). Allele frequencies were calculated by direct count in the raw genotype data file, expected heterozygosities per locus (He) were estimated by He = 2pq and PIC values were estimated by PIC = H-2p2q2, where p and q are the allele frequencies of the two alternate alleles.

Format: XLS Size: 548KB Download file

This file can be viewed with: Microsoft Excel ViewerOpen Data

Table 2. Minor allele frequency of validated SNPs

Table 3. Distribution of the genotyped SNPs in contigs

Almost 36 % of the successfully genotyped SNPs were homozygous in all samples (i.e. only one SNP variant was identified in all individuals). These were incorrectly identified as SNPs since EST sequencing presents a high rate of sequencing errors resulting in the identification of pseudo-SNPs. Some of these SNPs may also correspond to rare polymorphisms that were not present in the population samples genotyped in this study. ESTs are issued from a wide variety of tissues usually collected from a limited number of individuals and genetically constrained populations and may represent a bias in true allelic variations. The effectiveness of EST resources to detect in silico SNPs highly depends on the collection of tissues used and the diversity of the target samples as well as on how well this diversity is represented within the EST databases used for SNP identification [37,38].

Two hundred and sixty two SNPs (27 %) revealed paralogous sequences as all samples, including DH individuals, were heterozygous. Out of these, 63 % were identified in contigs containing at least six SNPs/contig (Table 3). Since up to 50 % of the rainbow trout genome could have retained duplicated regions, the high proportion is most probably due to the assembly of duplicated gene sequences which could result in the production of paralogous site variants (PSVs). PSVs are sequence differences between two paralogous loci but the substitution does not segregate within either locus and were considered false positive SNPs. Similar observations were obtained with Atlantic salmon [10,12,13].

Finally, 37 % (350) of the successfully genotyped SNPs were polymorphic and reliably scored, and thus were considered as true SNPs (Table 2). They were identified in 321 contigs. The yield of true SNPs decreased with the number of SNPs/contig and almost 75 % of those true SNPs were identified from contigs containing no more than five SNPs/contig (Table 3). The mean of observed minor allele frequency (MAF) among true SNPs was 0.27 ± 0.14 (SD), while the mean observed heterozygosity (Figure 2) across loci was 0.35 ± 0.14, and the mean PIC (Figure 3) was 0.28 ± 0.1 (see also Additional file 3: Sheet 6). Since observed heterozygosity and PIC rates are near the maximum theoretical values for a bi-allelic marker, we can conclude that the validation panel is highly informative for this type of markers. These SNPs are of particular interest for linkage analysis since we can easily follow up their segregation from one generation to another.

thumbnailFigure 2. Distribution of observed heterozygosity for true SNPs in three populations. SNPs were clustered into categories based on their observed heterozygosity values.

thumbnailFigure 3. Frequency distribution of PIC values across the three genotyped populations. SNPs were clustered into categories based on their observed PIC values.

Population assignment

Three domestic populations of different origins (INRA SP and SY strains and NCCCWA population) were used in this study. This offers the opportunity to determine whether a Bayesian clustering software such as STRUCTURE could detect the underlying genetic populations among all analysed samples using the observed SNP genotypes only.

We first used the non parametric approach [34] to infer the optimal number of populations (true K value). Inference of the best K using the delta K method revealed a clear peak at K = 3 (Figure 4 and Additional file 4) which corresponds to the true number of populations used in the present study. For this K value, STRUCTURE software successfully sorted individuals into three main groups which corresponded entirely to the discrete three main populations sampled in the study (Figure 5). These results should assist the design of SNP- rather than microsatellite-based studies to detect population structure in a larger collection of rainbow trout.

Additional file 4. Best k analysis. Results were shown for the five independent runs using a bur-in of 50,000 and 200,000 iterations. The best K value was determined using the non parametric test as previously described [29].

Format: XLS Size: 35KB Download file

This file can be viewed with: Microsoft Excel ViewerOpen Data

thumbnailFigure 4. Prediction of the best value of K. Delta K analysis was performed as previously described (Evanno, 2005) in order to predict the best value of K.

thumbnailFigure 5. Genetic population structure predicted by STRUCTURE software. Genetic population structure was inferred with the structure software using K = 3, a burn-in period of 50,000 200,000 iterations for the likelihood estimation. Individuals 1 to 20 correspond to the INRA-SP population, individuals 21 to 40 correspond to the INRA-SY population and individuals 41 to 84 correspond to the NCCCWA population.

Even though microsatellites have higher allele diversity, the frequent occurrence of mutations following a stepwise mutation model within these markers may lead to homoplastic alleles which may represent a significant problem in population genetics [39]. SNPs are considered biallelic and individual SNP loci have lower information content than microsatellite. However, they are highly frequent within genomes, have low mutation rates and these features allow reconstituting highly informative and non homoplasic haplotypes. With the advent of high throughput genotyping strategies, SNPs open new avenues in population genetics such as association studies in natural populations. From a practical point of view, they offer more rapid, highly automated and more reliable genotyping which are also useful properties for population inference analysis.

Functional annotation of true SNPs

To assign putative functions to the 350 true SNPs, we performed blastx searches of both SNPs and corresponding EST contig sequences against the ENSEMBL zebrafish non redundant peptide database.

Blastx search results made it possible to assign putative functions to 321 contig sequences (Additional file 5: Sheet 1). Of these, 279 contig sequences showed unique gene hits, 35 contig sequences showed unique gene hits but multiple alignment positions and seven contig sequences had multiple paralogous gene hits.

Additional file 5. SNP functional annotations. BlastX search results using EST contig (workseet 1) and SNP (worksheet 2) sequences as well as synonymous/non-synonymous SNP prediction results (workseet 3) were summerized.

Format: XLS Size: 278KB Download file

This file can be viewed with: Microsoft Excel ViewerOpen Data

Blastx searches using true SNP sequences revealed that 339 markers resulted in the same gene hits as those found with their corresponding contig sequences (Additional file 5: Sheet 2). Among these, 11 SNPs did not show any homology search results and 12 SNPs matched with target regions in gene hits different from those found with the corresponding contig sequences. These 23 markers were excluded from the SNP panel used for synonymous/non-synonymous prediction analyses. Out of the remaining 327 SNPs, 28 were identified as non-synonymous and four of these substitutions were predicted to affect protein functions (Additional file 5: Sheet 3). This is particularly important since they could be considered as valuable sources of candidate gene polymorphisms underlying important traits leading to the identification of causative genes. However, these predictions were conducted using computational tools and functional data analyses are therefore needed to validate the most likely effect of these substitutions on the protein functions.

Non-synonymous SNPs are of particular interest because they are more likely to alter the biological function of a protein. They are suitable markers for comparative genome mapping and for marker-assisted selection of economically important traits [40,41].

Even though synonymous SNPs have long been considered as silent substitutions, they are also of particular interest since they can alter RNA secondary structures and affect regulation of gene expression [1].

Transitions/transversions ratio

About 74.1 % of SNPs from the validation panel were A ↔ G and C ↔ T transitions representing 28.2 % and 45.9 of the total SNPs, respectively (Table 4). For both true and synonymous SNPs, observed transition over transversion (Ts/Tv) ratios were (3.48) and (4.07), respectively. In addition, the Ts/Tv ratio was found to be higher in synonymous (4.07) than in non-synonymous SNPs (1.15).

Table 4. Observed occurrence of the different SNP types (transitions and transversions)

On average, an excess of transitions was observed in this study, which is believed to be attributable to the abundant hypermutable methylated dinucleotide 5′-CpG-3 [42]. One probable explanation would be the high spontaneous rate of deamination of 5'-methylated cytosines (5mC) at CpG di-nucleotides to thymidine (C↔T) SNPs and (G↔A) on the complementary strand [43].

Synonymous SNPs are more often transitions than transversions. It is generally agreed that degeneracy in the genetic code and the stronger selection pressure against non-synonymous substitutions account for the observed increase in the relative frequencies of transitions over transversions [44].

It is difficult to draw general conclusion from the smallnumber of non-synonymous substitutions recorded here. However, transition and transversion classes tended to occur at similar levels (Table 5). Both transition and transversion classes can result in amino acid substitutions, but the biochemical differences in the corresponding protein products tend to be greater for transversions [45].

Table 5. Distribution of SNPs within INRA rainbow trout microsatellite-based linkage map

Integration of EST-derived SNPs in the INRA RT linkage map

Two unrelated reference families (two parents with four grandparents and 120 progeny) were successfully genotyped with the 350 true SNPs. Out of these, 120 SNPs were not polymorphic in the two reference families and were therefore excluded from the final panel for linkage mapping. Segregation pattern analyses of the remaining 230 SNPs revealed seven SNPs showing non Mendelian inheritance. These were excluded from the final panel for linkage mapping. The remaining 223 SNPs were successfully positioned on the microsatellite-based linkage map (Additional file 6). The number of SNPs that were assigned to each linkage group (RT) varied from 0 (RT17) to 17 (RT27). Linkage groups RT04 and RT25 were unlinked, but were artificially merged to form a metacentric linkage group, as previously reported [28,46].

Additional file 6. INRA rainbow trout linkage maps. Newly integrated SNPs were highlighted in red.

Format: PDF Size: 730KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

One hundred markers (29 %) were identified within duplicated regions or adjacent to duplicated markers delimiting known duplicated regions of rainbow trout (Additional file 7). Four SNPs were found duplicated, out of which three pairs (snp_BX855659_100/1 and /2, snp_CX719996_684/1 and /2, snp_F3GS49K01ARLCI_353/1 and /2) were successfully mapped to the expected homeologous groups (Additional file 6). The fourth pair (snp_BX078786_199/1 and /2) were assigned to linkage groups RT01 and RT19, which did not show any homology in previous studies [28,29,46]. Additionally, three pairs of SNPs (snp_BX869590_239 and snp_CA387137_459, snp_GBPNQUD01ALF2N_185 and snp_GBPNQUD01BTTCT_137, snp_F3GS49K01CUTTZ_151 and snp_F3GS49K02I3USJ_458) were originally identified within the same EST contigs. For each pair, one SNP was assigned to a linkage group and the second one to the homeologous group (RT06 and RT30, RT12 and RT16 and RT07 and RT15, respectively). These seven markers were true but duplicated SNPs and were considered multisite variants (MSVs). Since up to 50 % of the rainbow trout pseudotetraploid genome still retains duplicated regions, this can lead to the production of MSVs. These are SNPs within paralogous loci with one locus fixed (monomorphic) and the other one polymorphic. MSVs can be used for both linkage and association studies and they can provide information on two paralogous loci resulting from the latest WGD event.

Additional file 7. List of duplicated markers and markers assigned in duplicated region.

Format: XLS Size: 113KB Download file

This file can be viewed with: Microsoft Excel ViewerOpen Data

Conclusions

One of the main objectives of rainbow trout genomic research is the development of rapid, accurate and automated genotyping systems for sequence variations influencing economically important traits. Developing larger sets of SNP markers for genome analyses in rainbow trout will facilitate fine QTL mapping and will improve the identification and exploitation of genes affecting important traits and enable selective breeding through genomic selection.

In the lack of a rainbow trout reference genome, EST resources represent an attractive approach for in silico SNP identification. In this study, we have assessed the quality of 1,152 rainbow trout EST-derived SNPs by genotyping several population samples using the Illumina GoldenGate BeadXpress assays. High-quality genotype data were obtained from 958 SNPs representing a genotyping success rate of 83.2 %. Polymorphism information was validated for 350 (36.5 %) rainbow trout EST-derived SNPs.

Almost two thirds of the successfully genotyped SNPs were considered as either false SNP (homozygous in all individuals) or duplicated heterozygous in double haploid individuals, expected to be all homozygous. This higher failure rate is most probably due to EST sequencing errors or the assembly of duplicated gene sequences, which can lead to the production of MSVs (true but duplicated SNPs) and PSVs (false positives). Although MSVs and PSVs can provide information on paralogous loci resulting from the latest WGD event in rainbow trout, their presence can hinder the use of SNP arrays in genome-wide association and should be avoided in population genetics studies. MSVs are of particular interest since they can be used for both mapping, linkage, and for association studies. However, since these markers are not easily interpreted, they were either discarded or their genotypes were incorrectly assigned. Since up to half of the rainbow trout genome still retains duplicated regions, calling and assigning MSVs to the correct paralogous loci become difficult and hinder genotype calling and genetic mapping. Moreover, currently available standard genotyping software programs have not been designed to automatically deal with this type of polymorphisms. It is noteworthy that nearly all the SNPs which were mapped in the duplicated regions were assigned to the same linkage group suggesting that, in most cases, only one of the two duplicated sites could be genotyped. This is particularly important for most of populations genetics studies where the same set of markers must be analysed.

The Illumina GenomeStudio software was recently updated and provides the GT module v2010/v1.8 which allows assigning genotypes in tetraploid genomes. However, users need to perform manual clustering of genotype calls for polyploid loci. This is unrealistic to accomplish when thousands of genotyped markers are analysed simultaneously. Consequently, an accurate and automated genotyping system for calling polyploid genotypes is highly desirable. A recent genotyping software, beadarrayMSV, has been designed and was used for genotype calling and mapping of MSVs in Atlantic salmon [47]. Nevertherless, this tool was rather useful for manipulating data genotyped using Illumina Infinium BeadArrays (Illumina Inc.) rather than GoldenGate data and therefore could not be used in the present study.

Combining the identification and validation of MSV markers with regular SNP polymorphisms in rainbow trout will allow the identification of fully informative genetic markers in this species and will provide valuable tools to efficiently design a future rainbow trout specific SNP chip.

Results reported here represent the first study of EST-derived SNP validation in rainbow trout and demonstrate the utility of EST databases as an alternative approach for de novo SNP identification in species whose genome sequences are not yet available.

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

MB conceived and supervised the project, and drafted the manuscript. RG carried out the linkage maps and genetic diversity analysis, and revised the manuscript. CC performed EST-derived SNP identification and Sigenae SNP database management. DE carried out genotyping work. EQ contributed in designing the study, helped with interpretation of data analyses and revised the manuscript. All authors read and approved the final manuscript.

Acknowledgements

This study has been financially supported by the INRA Department of Animal Genetics. Genotyping work was carried out at the INRA genomics GET PlaGe core facility at Toulouse. DNA samples from 10 DH individuals and NCCCWA population were generously provided by Dr Thorgaard and Dr Yniv Palti (ARS), respectively. Sample DNA extraction for genotyping was performed by Aurélie Charlet and Céline Ciobotaru. We are grateful to Franck Samson from the INRA MIGALE bioinformatics platform ( http://migale.jouy.inra.fr) for providing help and support.

References

  1. Sauna ZE, Kimchi-Sarfaty C, Ambudkar SV, Gottesman MM: The sounds of silence: synonymous mutations affect function.

    Pharmacogenomics 2007, 8(6):527-532. PubMed Abstract | Publisher Full Text OpenURL

  2. Kim H, Schmidt CJ, Decker KS, Emara MG: A double-screening method to identify reliable candidate non-synonymous SNPs from chicken EST data.

    Anim Genet 2003, 34(4):249-254. PubMed Abstract | Publisher Full Text OpenURL

  3. Garg K, Green P, Nickerson DA: Identification of candidate coding region single nucleotide polymorphisms in 165 human genes using assembled expressed sequence tags.

    Genome Res 1999, 9(11):1087-1092. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  4. Consortium WTCC: Genome-wide association study of 14,000 cases of seven common diseases and 3,000 shared controls.

    Nature 2007, 447(7145):661-678. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  5. Valdar W, Solberg LC, Gauguier D, Burnett S, Klenerman P, Cookson WO, Taylor MS, Rawlins JN, Mott R, Flint J: Genome-wide genetic association of complex traits in heterogeneous stock mice.

    Nat Genet 2006, 38(8):879-887. PubMed Abstract | Publisher Full Text OpenURL

  6. Muir WM, Wong GK, Zhang Y, Wang J, Groenen MA, Crooijmans RP, Megens HJ, Zhang H, Okimoto R, Vereijken A, et al.: Genome-wide assessment of worldwide chicken SNP genetic diversity indicates significant absence of rare alleles in commercial breeds.

    Proc Natl Acad Sci U S A 2008, 105(45):17312-17317. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  7. Kim ES, Berger PJ, Kirkpatrick BW: Genome-wide scan for bovine twinning rate QTL using linkage disequilibrium.

    Anim Genet 2009, 40(3):300-307. PubMed Abstract | Publisher Full Text OpenURL

  8. Kijas JW, Townley D, Dalrymple BP, Heaton MP, Maddox JF, McGrath A, Wilson P, Ingersoll RG, McCulloch R, McWilliam S, et al.: A genome wide survey of SNP variation reveals the genetic structure of sheep breeds.

    PLoS One 2009, 4(3):e4668. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  9. Wang S, Peatman E, Abernathy J, Waldbieser G, Lindquist E, Richardson P, Lucas S, Wang M, Li P, Thimmapuram J, et al.: Assembly of 500,000 inter-specific catfish expressed sequence tags and large scale gene-associated marker development for whole genome association studies.

    Genome Biol 2010, 11(1):R8. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  10. Hayes BJ, Nilsen K, Berg PR, Grindflek E, Lien S: SNP detection exploiting multiple sources of redundancy in large EST collections improves validation rates.

    Bioinformatics 2007, 23(13):1692-1693. PubMed Abstract | Publisher Full Text OpenURL

  11. Moen T, Hayes B, Baranski M, Berg PR, Kjoglum S, Koop BF, Davidson WS, Omholt SW, Lien S: A linkage map of the Atlantic salmon (Salmo salar) based on EST-derived SNP markers.

    BMC Genomics 2008, 9:223. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  12. Ryynanen HJ, Primmer CR: Single nucleotide polymorphism (SNP) discovery in duplicated genomes: intron-primed exon-crossing (IPEC) as a strategy for avoiding amplification of duplicated loci in Atlantic salmon (Salmo salar) and other salmonid fishes.

    BMC Genomics 2006, 7:192. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  13. Smith CT, Elfstrom CM, Seeb LW, Seeb JE: Use of sequence data from rainbow trout and atlantic salmon for SNP detection in pacific salmon.

    Mol Ecol 2005, 14(13):4193-4203. PubMed Abstract | Publisher Full Text OpenURL

  14. Jalving R, Oost BA, Slot R: Chicken single nucleotide polymorphism identification and selection for genetic mapping.

    Poult Sci 2004, 83(12):1925-1931. PubMed Abstract | Publisher Full Text OpenURL

  15. Marth GT, Korf I, Yandell MD, Yeh RT, Gu Z, Zakeri H, Stitziel NO, Hillier L, Kwok PY, Gish WR: A general approach to single-nucleotide polymorphism discovery.

    Nat Genet 1999, 23(4):452-456. PubMed Abstract | Publisher Full Text OpenURL

  16. Thorgaard GH, Bailey GS, Williams D, Buhler DR, Kaattari SL, Ristow SS, Hansen JD, Winton JR, Bartholomew JL, Nagler JJ, et al.: Status and opportunities for genomics research with rainbow trout.

    Comp Biochem Physiol B Biochem Mol Biol 2002, 133(4):609-646. PubMed Abstract | Publisher Full Text OpenURL

  17. Ng SH, Artieri CG, Bosdet IE, Chiu R, Danzmann RG, Davidson WS, Ferguson MM, Fjell CD, Hoyheim B, Jones SJ, et al.: A physical map of the genome of Atlantic salmon, Salmo salar.

    Genomics 2005, 86(4):396-404. PubMed Abstract | Publisher Full Text OpenURL

  18. Young WP, Wheeler PA, Coryell VH, Keim P, Thorgaard GH: A detailed linkage map of rainbow trout produced using doubled haploids.

    Genetics 1998, 148(2):839-850. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  19. Allendorf FW: TG. Turner BJ: Evolutionary Genetics of Fishes; 1984. OpenURL

  20. Palti Y, Gahr SA, Hansen JD, Rexroad CE: Characterization of a new BAC library for rainbow trout: evidence for multi-locus duplication.

    Anim Genet 2004, 35(2):130-133. PubMed Abstract | Publisher Full Text OpenURL

  21. Le Bras Y, Dechamp N, Krieg F, Filangi O, Guyomard R, Boussaha M, Bovenhuis H, Pottinger TG, Prunet P, Le Roy P, et al.: Detection of QTL with effects on osmoregulation capacities in the rainbow trout (Oncorhynchus mykiss).

    BMC Genet 2011, 12:46. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  22. Sanchez CC, Smith TP, Wiedmann RT, Vallejo RL, Salem M, Yao J, Rexroad CE: Single nucleotide polymorphism discovery in rainbow trout by deep sequencing of a reduced representation library.

    BMC Genomics 2009, 10:559. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  23. Miller MR, Brunelli JP, Wheeler PA, Liu S, Rexroad CE, Palti Y, Doe CQ, Thorgaard GH: A conserved haplotype controls parallel adaptation in geographically distant salmonid populations.

    Mol Ecol 2011, 21(2):237. PubMed Abstract | Publisher Full Text OpenURL

  24. Sanchez CC, Weber GM, Gao G, Cleveland BM, Yao J, Rexroad CE: Generation of a reference transcriptome for evaluating rainbow trout responses to various stressors.

    BMC Genomics 2011, 12(1):626. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  25. Genotrout Séquençage du génome de la truite arc-en-ciel par l’utilisation des nouvelles techniques de séquençage à ultra haut déb.

    http://www4.rennes.inra.fr/scribe/Les-recherches/Programmes-nationaux-et-internationaux-en-cours/Genotrout webcite

    PubMed Abstract OpenURL

  26. Production of a Draft Reference Genome Sequence for Rainbow Trout.

    http://www.ars.usda.gov/research/projects/projects.htm?ACCN_NO = 419700 webcite

    OpenURL

  27. Genet C, Dehais P, Palti Y, Gao G, Gavory F, Wincker P, Quillet E, Boussaha M: Analysis of BAC-end sequences in rainbow trout: content characterization and assessment of synteny between trout and other fish genomes.

    BMC Genomics 2011, 12:314. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  28. Rexroad CE, Palti Y, Gahr SA, Vallejo RL: A second generation genetic map for rainbow trout (Oncorhynchus mykiss).

    BMC Genet 2008, 9:74. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  29. Guyomard R, Mauger S, Tabet-Canale K, Martineau S, Genet C, Krieg F, Quillet E: A type I and type II microsatellite linkage map of rainbow trout (Oncorhynchus mykiss) with presumptive coverage of all chromosome arms.

    BMC Genomics 2006, 7:302. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  30. Lin CH, Yeakley JM, McDaniel TK, Shen R: Medium- to high-throughput SNP genotyping using VeraCode microbeads.

    Methods Mol Biol 2009, 496:129-142. PubMed Abstract | Publisher Full Text OpenURL

  31. Fan JB, Oliphant A, Shen R, Kermani BG, Garcia F, Gunderson KL, Hansen M, Steemers F, Butler SL, Deloukas P, et al.: Highly parallel SNP genotyping.

    Cold Spring Harb Symp Quant Biol 2003, 68:69-78. PubMed Abstract | Publisher Full Text OpenURL

  32. Pritchard JK, Stephens M, Donnelly P: Inference of population structure using multilocus genotype data.

    Genetics 2000, 155(2):945-959. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  33. Falush D, Stephens M, Pritchard JK: Inference of population structure using multilocus genotype data: linked loci and correlated allele frequencies.

    Genetics 2003, 164(4):1567-1587. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  34. Evanno G, Regnaut S, Goudet J: Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study.

    Mol Ecol 2005, 14(8):2611-2620. PubMed Abstract | Publisher Full Text OpenURL

  35. de Givry S, Bouchez M, Chabrier P, Milan D, Schiex T: Carhta Gene: multipopulation integrated genetic and radiation hybrid mapping.

    Bioinformatics 2005, 21(8):1703-1704. PubMed Abstract | Publisher Full Text OpenURL

  36. Voorrips RE: MapChart: software for the graphical presentation of linkage maps and QTLs.

    J Hered 2002, 93(1):77-78. PubMed Abstract | Publisher Full Text OpenURL

  37. Rafalski A: Applications of single nucleotide polymorphisms in crop genetics.

    Curr Opin Plant Biol 2002, 5(2):94-100. PubMed Abstract | Publisher Full Text OpenURL

  38. Picoult-Newberg L, Ideker TE, Pohl MG, Taylor SL, Donaldson MA, Nickerson DA, Boyce-Jacino M: Mining SNPs from EST databases.

    Genome Res 1999, 9(2):167-174. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  39. Estoup A, Jarne P, Cornuet JM: Homoplasy and mutation model at microsatellite loci and their consequences for population genetics analysis.

    Mol Ecol 2002, 11(9):1591-1604. PubMed Abstract | Publisher Full Text OpenURL

  40. Sarropoulou E, Nousdili D, Magoulas A, Kotoulas G: Linking the genomes of nonmodel teleosts through comparative genomics.

    Mar Biotechnol (NY) 2008, 10(3):227-233. Publisher Full Text OpenURL

  41. Wang S, Sha Z, Sonstegard TS, Liu H, Xu P, Somridhivej B, Peatman E, Kucuktas H, Liu Z: Quality assessment parameters for EST-derived SNPs from catfish.

    BMC Genomics 2008, 9:450. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  42. Holliday R, Grigg GW: DNA methylation and mutation.

    Mutat Res 1993, 285(1):61-67. PubMed Abstract | Publisher Full Text OpenURL

  43. Duncan BK, Miller JH: Mutagenic deamination of cytosine residues in DNA.

    Nature 1980, 287(5782):560-561. PubMed Abstract | Publisher Full Text OpenURL

  44. Moriyama EN, Powell JR: Intraspecific nuclear DNA variation in drosophila.

    Mol Biol Evol 1996, 13(1):261-277. PubMed Abstract | Publisher Full Text OpenURL

  45. Zhang J: Rates of conservative and radical nonsynonymous nucleotide substitutions in mammalian nuclear genes.

    J Mol Evol 2000, 50(1):56-68. PubMed Abstract | Publisher Full Text OpenURL

  46. Danzmann RG, Davidson EA, Ferguson MM, Gharbi K, Koop BF, Hoyheim B, Lien S, Lubieniecki KP, Moghadam HK, Park J, et al.: Distribution of ancestral proto-actinopterygian chromosome arms within the genomes of 4R-derivative salmonid fishes (Rainbow trout and Atlantic salmon).

    BMC Genomics 2008, 9:557. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  47. Gidskehaug L, Kent M, Hayes BJ, Lien S: Genotype calling and mapping of multisite variants using an atlantic salmon iSelect SNP array.

    Bioinformatics 2011, 27(3):303-310. PubMed Abstract | Publisher Full Text OpenURL