- Data note
- Open Access
De novo assembly of the Brown trout (Salmo trutta m. fario) brain and muscle transcriptome: transcript annotation, tissue differential expression profile and SNP discovery
BMC Research Notes volume 13, Article number: 503 (2020)
The Brown trout is a salmonid species with a high commercial value in Europe. Life history and spawning behaviour include resident (Salmo trutta m. fario) and migratory (Salmo trutta m. trutta) ecotypes. The main objective is to apply RNA-seq technology in order to obtain a reference transcriptome of two key tissues, brain and muscle, of the riverine trout Salmo trutta m. fario. Having a reference transcriptome of the resident form will complement genomic resources of salmonid species.
We generate two cDNA libraries from pooled RNA samples, isolated from muscle and brain tissues of adult individuals of Salmo trutta m. fario, which were sequenced by Illumina technology. Raw reads were subjected to de-novo transcriptome assembly using Trinity, and coding regions were predicted by TransDecoder. A final set of 35,049 non-redundant ORF unigenes were annotated. Tissue differential expression analysis was evaluated by Cuffdiff. A False Discovery Rate (FDR) ≤ 0.01 was considered for significant differential expression, allowing to identify key differentially expressed unigenes. Finally, we have identified SNP variants that will be useful tools for population genomic studies.
Brown trout (Salmo trutta) has been extensively studied by its commercial and biological importance. From the sixty-six species in this family, S. trutta is a species native to Europe with a wide distribution area that includes Atlantic and Mediterranean European basins, as well as northern Africa and western Asia basins [1, 2].The specie has been introduced in North and South America and Australia by its commercial exploitation for sport fishing, as well as farmed for food and game fish, extending their actual geographical distribution as discontinuous populations on all continents except Antarctica .
Life history traits of Brown trout populations include resident forms such as riverine (S. trutta m. fario) and migratory forms such as anadromous (S. trutta m. trutta) ecotype [4, 5]. Anadromous and non-anadromous forms coexist in the same river being apparently genetically indistinguishable [6, 7]. An extended literature on Brown trout research has been produced that includes physiological, ecological and genetic aspects [8,9,10]. As a contribution to this global effort, here we provide a comprehensive transcriptome data set derived from brain and muscle tissues of Salmo trutta m. fario ecotype by using RNA-seq technology. We also evaluated differential transcript expression among these two tissues identifying key differentially expressed unigenes. Finally, we applied an in-silico pipeline that allow us to discover SNP variants useful for population genomic studies. The generated data could provide new valuable genomic resources for population genetic and genomic studies that can help to answer opened questions about the live history traits of riverine S. trutta m. fario as well as differences among S. trutta ecotypes.
Salmo trutta m. fario. brain and muscle tissues were collected from 25 wild type individuals (15 females) captured at the Falmisell river (Lleida, Catalonia). RNA pools from brain (10.2 µg) and muscle (11.4 µg) tissues were obtained with equimolar concentration from each subject. The TruSeq™ RNA sample Prep Kit (Illumina, Madrid, Spain) was used to build cDNA libraries according to manufacturer instructions (Table 1, Data file 1). FASTQ sequence reads were assembled using Trinity  run on the paired end sequences with the fixed default k-mer size of 25 and minimum contig length of 200. Descriptive statistics of assembly and sequencing is found at Table 1 (Data file 2 and Data file 3). Among the 144,984 contigs predicted by Trinity (Table 1, Data file 4 and Data file 8), we identify protein coding regions using TransDecoder package . We retained the longest ORF predicted for each contig sequence with a minimum of 100 amino acids long. Transcript redundancy was further reduced by CD-hit , obtaining a final set of 35,189 non-redundant ORF unigenes as best cluster representatives (Table 1, Data file 5). Size distribution for clustered ORF unigenes is presented in Table 1 (Data file 3). This final set was characterized by homology search to nucleotide and protein databases (Table 1, Data file 10 and Data file 11). Taxonomic representation showed the top hits for a large fraction of unigenes (≈88%) to Neopterigii taxon, with 66% of unigenes assigned to family Salmonidae (Salvelius sp. (1%), Onchorrinchus sp. (14%) and Salmo sp. (51%) (Table 1, Data file 12). A total of 4337 protein motif were assigned to 23,616 ORF unigenes, being the RNA recognition motif (6.4%), Immunoglobulin domain (4.8%), Tetratricopeptide repeat (4.8%) and Protein kinase domain (3.4%) the most prevalent (Table 1, Data file 13).
Similarity search by Blast2GO renders a total of 28,132 (80%) unigenes with GO annotation. GO term were then simplified using a generic GOSlim vocabulary  (Table 1, Data file 14). The ten top GO terms among the Cellular Component (18,071, 64%), Molecular Function (20,691, 74%) and Biological Process (23,954, 85%) ontology at level 2 are shown in Table 1 (Data file 4). Mapping unigenes to the reference canonical pathways in the KEGG database, yields a total of 13,957 (39.8%) ORF unigenes assigned to 3421 KEGG terms (KO) defining a total of 386 pathways (Table 1, Data file 15).
Tissue specific transcriptome expression analysis was performed by normalization of raw reads (FPKM, fragments per kilobase of exon per million fragments) obtained from both tissues (Table 1, Data file 16 and Data file 17). Analysis reveals 1172 ORF unigenes expressed only in muscle, 8595 expressed only in brain and 12,072 expressed in both tissues (Table 1, Data file 5, FigS3). Differentially expressed unigenes at FDR < 0.01 and best homologous sequences are shown at Table 1 (Data file 18 and Data file 19).
Finally, we have identified 73,237 putative SNPs (Table 1, Data file 20) and extracted 150 bp sequence context to each SNP as a source for the design of PCR primers useful for genotyping protocols (Table 1, Data file 21).
The use of pooled RNA samples does not allow us to detect sex or individual specific transcript expression profiles as well as limit our capability to detect transcripts expressed at low level in a specific individual. In addition, pooled samples avoid us to resolve SNP frequency distribution, being this parameter indirectly estimated according to the observed SNP sequence coverage in the pooled sample.
Availability of data and materials
The data described in this Data note can be freely and openly available on Figshare (https://0-doi-org.brum.beds.ac.uk/10.6084/m9.figshare.12902474.v1; https://0-doi-org.brum.beds.ac.uk/10.6084/m9.figshare.12902405.v2; https://0-doi-org.brum.beds.ac.uk/10.6084/m9.figshare.7326464.v1; https://0-doi-org.brum.beds.ac.uk/10.6084/m9.figshare.7712708.v4; https://0-doi-org.brum.beds.ac.uk/10.6084/m9.figshare.12905777.v2; https://0-doi-org.brum.beds.ac.uk/10.6084/m9.figshare.12905747.v1; https://0-doi-org.brum.beds.ac.uk/10.6084/m9.figshare.12905831.v1). Assembly of non-redundant ORF unigene sequences are available from NCBI transcriptome shotgun assembly (TSA) database (https://identifiers.org/ncbi/insdc:GHGR00000000.1). Raw sequence reads are available from the NCBI sequence read archive (SRA) database (https://identifiers.org/insdc.sra:SRP151838). Please see Table 1 and references list [14,15,16,17,18,19,20,21,22] for details and links to the data.
Binary Sequence Alignment/Map
Basic local alignment search tool, bp: base pair
Fragments Per Kilobase of exon model per Million mapped reads
Kyoto Encyclopaedia of Genes and Genomes
Open reading frame
Protein families database
False Discovery Rate
Sequence Alignment/Map tools
Single nucleotide polymorphism
Bagliniere JL. Introduction: the brown trout (Salmo trutta L)—its origin, distribution and economic and scientific significance. Biology and ecology of the Brown and Sea Trout. 3rd ed. London: Springer London; 2000. pp. 1–12. https://0-doi-org.brum.beds.ac.uk/https://0-doi-org.brum.beds.ac.uk/10.1007/978-1-4471-0775-0_1
Klemetsen A, Amundsen PA, Dempson JB, Jonsson B, Jonsson N, O'Connell MF, et al. Atlantic salmon Salmo salar L., brown trout Salmo trutta L. and Arctic charr Salvelinus alpinus (L.): a review of aspects of their life histories. Ecol Freshwater Fish. 2nd ed. 2003;12:1–59. https://0-doi-org.brum.beds.ac.uk/https://0-doi-org.brum.beds.ac.uk/10.1034/j.1600-0633.2003.00010.x
MacCrimmon HR, Marshall TL. World distribution of Brown Trout, Salmo trutta. J Fish Res Board Can. 2011;25:2527–48. https://0-doi-org.brum.beds.ac.uk/10.1139/f68-225.
Elliott JM. Quantitative ecology and the Brown Trout. Oxford University Press, USA; 1994. https://0-doi-org.brum.beds.ac.uk/https://0-doi-org.brum.beds.ac.uk/10.1577/1548-8659-123.6.1006
Poćwierz-Kotus A, Bernaś R, Dębowski P, Kent MP, Lien S, Kesler M, et al. Genetic differentiation of southeast Baltic populations of sea trout inferred from single nucleotide polymorphisms. Anim Genet. 2014;45:96–104. https://0-doi-org.brum.beds.ac.uk/10.1111/age.12095.
Charles K, Guyomard R, Hoyheim B, Ombredane D, Bagliniere JL. Lack of genetic differentiation between anadromous and resident sympatric brown trout (Salmo trutta) in a Normandy population. Aquat Living Resour. 2005;18:65–9. https://0-doi-org.brum.beds.ac.uk/10.1051/alr:2005006.
Charles K, Roussel JM, Lebel JM, Bagliniere JL, Ombredane D. Genetic differentiation between anadromous and freshwater resident brown trout (Salmo trutta L.): insights obtained from stable isotope analysis. Ecol Freshwater Fish. 2006;15:255–63. https://0-doi-org.brum.beds.ac.uk/10.1111/j.1600-0633.2006.00149.x.
Harvey J. Ecology of Atlantic Salmon and Brown Trout: habitat as a template for life histories. Freshw Biol. 2012;57:1531–41. https://0-doi-org.brum.beds.ac.uk/10.1007/978-94-007-1189-1.
Boel M, Aarestrup K, Baktoft H, Larsen T, Søndergaard Madsen S, Malte H, et al. The physiological basis of the migration continuum in brown trout (Salmo trutta). Physiol Biochem Zool. 2014;87:334–45. https://0-doi-org.brum.beds.ac.uk/10.1086/674869.
Oromi N, Jové M, Pascual-Pons M, Royo JL, Rocaspana R, Aparicio E, et al. Differential metabolic profiles associated to movement behaviour of stream-resident brown trout (Salmo trutta). PLoS ONE. 2017;12:e0181697. https://0-doi-org.brum.beds.ac.uk/10.1371/journal.pone.0181697.
Haas BJ, Papanicolaou A, Yassour M, Grabherr M, Blood PD, Bowden J, et al. De novo transcript sequence reconstruction from RNA-seq using the Trinity platform for reference generation and analysis. Nat Protoc. 2013;8:1494–512. https://0-doi-org.brum.beds.ac.uk/10.1038/nprot.2013.084.
Li W, Godzik A. Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics. 2006;22:1658–9. https://0-doi-org.brum.beds.ac.uk/10.1093/bioinformatics/btl158.
McCarthy FM, Bridges SM, Wang N, Magee GB, Williams WP, Luthe DS, et al. AgBase: a unified resource for functional analysis in agriculture. Nucleic Acids Res. 2007;35:D599-603. https://0-doi-org.brum.beds.ac.uk/10.1093/nar/gkl936.
Fibla J, Oromi N, Pascual-Pons M, Royo JL, Palau A and Fibla M. Salmo trutta m. fario Raw sequence reads. NCBI Sequence Read Archive; 2020. https://identifiers.org/insdc.sra:SRP151838
Fibla J, Oromi N, Pascual-Pons M, Royo JL, Palau A and Fibla M. TSA: Salmo trutta fario, transcriptome shotgun assembly. GenBank; 2020. https://identifiers.org/ncbi/insdc:GHGR00000000.1
Fibla J, Oromi N, Pascual-Pons M, Royo JL, Palau A and Fibla M. SNP discovery of Non-redundant ORF unigenes to Gene Ontology, KEGG and Protein Family databases of "De novo assembly of the Brown trout (Salmo trutta m. fario) brain and muscle transcriptome: transcript annotation, tissue differential expression profile and SNP discovery". Figshare; 2020. https://0-doi-org.brum.beds.ac.uk/10.6084/m9.figshare.12905831.v1
Fibla J, Oromi N, Pascual-Pons M, Royo JL, Palau A and Fibla M. Functional annotation of Non-redundant ORF unigenes to Gene Ontology, KEGG and Protein Family databases of "De novo assembly of the Brown trout (Salmo trutta m. fario) brain and muscle transcriptome: transcript annotation, tissue differential expression profile and SNP discovery". Figshare; 2020. https://0-doi-org.brum.beds.ac.uk/10.6084/m9.figshare.12905777.v2
Fibla J, Oromi N, Pascual-Pons M, Royo JL, Palau A and Fibla M. Tissue differential expression profile of "De novo assembly of the Brown trout (Salmo trutta m. fario) brain and muscle transcriptome: transcript annotation, tissue differential expression profile and SNP discovery". Figshare; 2020. https://0-doi-org.brum.beds.ac.uk/10.6084/m9.figshare.12905747.v1
Fibla J, Oromi N, Pascual-Pons M, Royo JL, Palau A and Fibla M. Supplementary Files for "De novo assembly of the Brown trout (Salmo trutta m. fario) brain and muscle transcriptome: transcript annotation, tissue differential expression profile and SNP discovery". Figshare; 2020. https://0-doi-org.brum.beds.ac.uk/10.6084/m9.figshare.12902474.v1
Fibla J, Oromi N, Pascual-Pons M, Royo JL, Palau A and Fibla M. Supplementary Figures for "De novo assembly of the Brown trout (Salmo trutta m. fario) brain and muscle transcriptome: transcript annotation, tissue differential expression profile and SNP discovery". Figshare; 2020. Figure. https://0-doi-org.brum.beds.ac.uk/10.6084/m9.figshare.12902405.v2
Fibla J, Oromi N, Pascual-Pons M, Royo JL, Palau A and Fibla M. Annotation of Non-redundant ORF unigenes to nucleotide and protein databases of "De novo assembly of the Brown trout (Salmo trutta m. fario) brain and muscle transcriptome: transcript annotation, tissue differential expression profile and SNP discovery". Figshare; 2020. https://0-doi-org.brum.beds.ac.uk/10.6084/m9.figshare.7712708.v4
Fibla J, Oromi N, Pascual-Pons M, Royo JL, Palau A, Fibla M. De novo assembled contigs. Figshare. 2020. https://0-doi-org.brum.beds.ac.uk/10.6084/m9.figshare.7326464.v1.
We are grateful to all participants of Gesna Estudis Ambientals S.L. (R. Rocaspana and E. Aparicio) and Eccus Proyectos Técnicos, Medioambientales y Obras SL (MA MarínVitalla), who have participated in the sampling procedure. We would like to dedicate this paper in memory of M.A. Marín Vitalla (Nines), who passed away last year. She is sorely missed.
This study has been supported and financed by the Biodiversity Conservation Plan of ENDESA, S.A. (ENEL Group) to JF. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Ethics approval and consent to participate
Permissions for electrofishing and capture of S. trutta m. fario individuals, was approved by the competent authorities: Departament de Medi Ambient i Habitatge de la Generalitat de Catalunya (current Departament d’Agricultura, Ramaderia, Pesca, Alimentacio i Medi Natural) (SF/602) of the regional authorities of Catalonia.
Consent for publication
The authors did not report any competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Fibla, J., Oromi, N., Pascual-Pons, M. et al. De novo assembly of the Brown trout (Salmo trutta m. fario) brain and muscle transcriptome: transcript annotation, tissue differential expression profile and SNP discovery. BMC Res Notes 13, 503 (2020). https://0-doi-org.brum.beds.ac.uk/10.1186/s13104-020-05351-4
- De novo transcriptome
- Brain & muscle transcriptome
- SNP discovery
- Salmo trutta m. fario