Generation of a reference transcriptome for evaluating rainbow trout responses to various stressors

Effect of Dietary Soy Protein Source on Effluent Water Quality and Growth Performance of Rainbow Trout Reared in a Serial Reuse Water System

Juvenile rainbow trout (125 ± 0.8 g) were fed a fishmeal control diet (C), a high soy protein concentrate diet (SP), a high soybean meal diet (HiS), or a diet with high levels of fermented soy protein concentrate (fSP) for 12 weeks in a tank system capable of receiving 1st and 3rd use water from a serial-reuse production hatchery. Water quality was generally lower in 3rd use compared to 1st use water and after passing through tanks (inflow vs. outflow). Total dissolved solids were significantly higher (p = 0.003) for 3rd use compared to 1st use water, and values were also higher (p < 0.001) for the fSP diet. Turbidity and ammonia were highest in tanks for trout fed the HiS and fSP and SP and fSP diets, respectively, but were characterized by high variation, which likely prevented the detection of significant differences. Weight gain (p < 0.001) and survival (p = 0.008) were significantly lower for trout in 3rd use compared to 1st use water. Trout fed the HiS diet were generally in poorer physiological condition with lower body fat stores (p = 0.05) and lower growth rate (p < 0.001) and survival (p = 0.05) compared to the other diets, which were similar. The expression of several stress-associated genes (FK506, DIO2, REGPS, Cyp1a, G6PH, GADD45a, and IRF-1) in the liver and gill showed that diet and water source affected their regulation. Replacement of FM by SP providing 50% of dietary protein promotes acceptable growth performance compared to an FM diet and was superior to HiS. The impacts of soy protein concentrate on water quality under commercial production conditions, however, require further study.

Comparative transcriptome sequencing analysis of female and male Decapterus macrosoma

Sexual growth dimorphism is a common phenomenon in teleost fish and has led to many reproductive strategies. Growth- and sex-related gene research in teleost fish would broaden our understanding of the process. In this study, transcriptome sequencing of shortfin scad Decapterus macrosoma was performed for the first time, and a high-quality reference transcriptome was constructed. After identification and assembly, a total of 58,475 nonredundant unigenes were obtained with an N50 length of 2,266 bp, and 28,174 unigenes were successfully annotated with multiple public databases. BUSCO analysis determined a level of 92.9% completeness for the assembled transcriptome. Gene expression analysis revealed 2,345 differentially expressed genes (DEGs) in the female and male D. macrosoma, 1,150 of which were female-biased DEGs, and 1,195 unigenes were male-biased DEGs. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses showed that the DEGs were mainly involved in biological processes including protein synthesis, growth, rhythmic processes, immune defense, and vitellogenesis. Then, we identified many growth- and sex-related genes, including Igf, Fabps, EF-hand family genes, Zp3, Zp4 and Vg. In addition, a total of 19,573 simple sequence repeats (SSRs) were screened and identified from the transcriptome sequences. The results of this study can provide valuable information on growth- and sex-related genes and facilitate further exploration of the molecular mechanism of sexual growth dimorphism.

Development of a High-Density 665 K SNP Array for Rainbow Trout Genome-Wide Genotyping

Single nucleotide polymorphism (SNP) arrays, also named « SNP chips », enable very large numbers of individuals to be genotyped at a targeted set of thousands of genome-wide identified markers. We used preexisting variant datasets from USDA, a French commercial line and 30X-coverage whole genome sequencing of INRAE isogenic lines to develop an Affymetrix 665 K SNP array (HD chip) for rainbow trout. In total, we identified 32,372,492 SNPs that were polymorphic in the USDA or INRAE databases. A subset of identified SNPs were selected for inclusion on the chip, prioritizing SNPs whose flanking sequence uniquely aligned to the Swanson reference genome, with homogenous repartition over the genome and the highest Minimum Allele Frequency in both USDA and French databases. Of the 664,531 SNPs which passed the Affymetrix quality filters and were manufactured on the HD chip, 65.3% and 60.9% passed filtering metrics and were polymorphic in two other distinct French commercial populations in which, respectively, 288 and 175 sampled fish were genotyped. Only 576,118 SNPs mapped uniquely on both Swanson and Arlee reference genomes, and 12,071 SNPs did not map at all on the Arlee reference genome. Among those 576,118 SNPs, 38,948 SNPs were kept from the commercially available medium-density 57 K SNP chip. We demonstrate the utility of the HD chip by describing the high rates of linkage disequilibrium at 2–10 kb in the rainbow trout genome in comparison to the linkage disequilibrium observed at 50–100 kb which are usual distances between markers of the medium-density chip.

Aquaculture omics: An update on the current status of research and data analysis

Aquaculture is the fast-growing agricultural sector and has the ability to meet the growing demand for protein nutritional security for future population. In future aquaculture is going to be the major source of fish proteins as capture fisheries reached at its maximum. However, several challenges need to overcome such as lack of genetically improved strains/varieties, lack of species-specific feed/functional feed, round the year availability of quality fish seed, pollution of ecosystems and increased frequencies of disease occurrence etc. In recent years, the continuous development of high throughput sequencing technology has revolutionized the biological sciences and provided necessary tools. Application of 'omics' in aquaculture research have been successfully used to resolve several productive and reproductive issues and thus ensure its sustainability and profitability. To date, high quality draft genomes of over fifty fish species have been generated and successfully used to develop large number of single nucleotide polymorphism markers (SNPs), marker panels and other genomic resources etc in several aquaculture species. Similarly, transcriptome profiling and miRNAs analysis have been used in aquaculture research to identify key transcripts and expression analysis of candidate genes/miRNAs involved in reproduction, immunity, growth, development, stress toxicology and disease. Metagenome analysis emerged as a promising scientific tool to analyze the complex genomes contained within microbial communities. Metagenomics has been successfully used in the aquaculture sector to identify novel and potential pathogens, antibiotic resistance genes, microbial roles in microcosms, microbial communities forming biofloc, probiotics etc. In the current review, we discussed application of high-throughput technologies (NGS) in the aquaculture sector.

Transcriptome Profiling Reveals a Divergent Adaptive Response to Hyper- and Hypo-Salinity in the Yellow Drum, Nibea albiflora

Simple Summary

Global warming and certain climate disasters (typhoon, tsunami, etc.) can lead to fluctuation in seawater salinity that causes salinity stress in fish. The aim of this study was to investigate the functional genes and relevant pathways in response to salinity stress in the yellow drum. Genes and pathways related to signal transduction, osmoregulation, and metabolism may be involved in the adaptive regulation to salinity in the yellow drum. Additionally, the genes under salinity stress were mainly divided into three expression trends. Our results provided novel insights into further study of the salinity adaptability of euryhaline fishes.

Abstract

The yellow drum (Nibea albiflora) is an important marine economic fish that is widely distributed in the coastal waters of the Northwest Pacific. In order to understand the molecular regulatory mechanism of the yellow drum under salinity stress, in the present study, transcriptome analysis was performed under gradients with six salinities (10, 15, 20, 25, 30, and 35 psu). Compared to 25 psu, 907, 1109, 1309, 18, and 243 differentially expressed genes (DEGs) were obtained under 10, 15, 20, 30, and 35 psu salinities, respectively. The differential gene expression was further validated by quantitative real-time PCR (qPCR). The results of the tendency analysis showed that all DEGs of the yellow drum under salinity fluctuation were mainly divided into three expression trends. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis showed that the PI3K-Akt signaling pathway, Jak-STAT signaling pathway as well as the glutathione metabolism and steroid biosynthesis pathways may be the key pathways for the salinity adaptive regulation mechanism of the yellow drum. G protein-coupled receptors (GPCRs), the solute carrier family (SLC), the transient receptor potential cation channel subfamily V member 6 (TRPV6), isocitrate dehydrogenase (IDH1), and fructose-bisphosphate aldolase C-B (ALDOCB) may be the key genes in the response of the yellow drum to salinity stress. This study explored the transcriptional patterns of the yellow drum under salinity stress and provided fundamental information for the study of salinity adaptability in this species.

Comprehensive analysis of transcriptomics and metabolomics to understand the flesh quality regulation of crucian carp (Carassius auratus) treated with short term micro-flowing water system

The short term micro-flowing purification system (STMFPS) has been shown to improve the flesh quality of freshwater fish. However, few studies have focused on the involved underlying mechanisms. This study explored the effect of STMFPS on the flesh quality of market-size freshwater fish based on the combination of metabolomics and transcriptomics methods. The UPLC-QTOF/MS based metabolomics method was utilized to screen metabolites and predict the possible major metabolic pathways during different STMFPS treatment periods (0 d, 1 d, 5 d and 9 d). Furthermore, the transcriptomic data demonstrated that the differentially expressed genes detected in crucian carp muscle were 2915, 7852 and 7183 after 1 d, 5 d and 9 d STMFPS treatment. Results showed that the TCA cycle, ornithine cycle, purine metabolism and amino acid catabolism play important roles in improving the flesh quality of crucian carp. This study may help to understand the mechanism of improving the flesh quality of aquatic products using STMFPS.

Structure and regulation of the NK-lysin (1-4) and NK-lysin like (a and b) antimicrobial genes in rainbow trout (Oncorhynchus mykiss)

Nk-lysin (Nkl), an antimicrobial peptide (AMP) product of natural killer cells and cytotoxic T cells in mammals, has recently been characterized in a number of finfish species. In this study, we identified six genes with sequence homology to Nkl and characterized their patterns of mRNA expression and abundances in rainbow trout (Oncorhynchus mykiss). The cDNA sequences for the six Nkls encoded precursor peptides of 128-133 aa in length, and mature peptides of 109-111 aa in length. Genomic DNA of the nkl1-4 genes consisted of five exons and four introns, whereas the nkl-like a & b genes consisted of four exons and three introns. Chromosomal locations of these peptides show that nkl1 was located on chromosome arm 25q, whereas the other five nkl genes were clustered on chromosome arm 19q. Phylogenetic analysis revealed a conserved structure of Nkls among the teleosts and further protein sequence analyses suggests that all six nkl genes fall within the Nkl sub-family of the Saposin family of proteins. Patterns of tissue-specific mRNA expression were asymmetric among the six trout Nkl homologues, with nkl1, nkl3, and nkl-like a & b occurring in immune competent organs such as spleen, gill, intestine and kidney, as well as pineal gland, brain and oocytes. However, nkl2 and nkl4, showed primary abundances in brain, pineal gland and oocyte tissues. Using mRNA sequencing, in whole-body pools of juvenile trout fry (1 g bw) exposed to Flavobacterium psychrophilum infection, we observed modest up-regulation (2-3 fold) of five (nkl 2-4 and nkl-like a & b) of the six nkl mRNAs over the five-day post-challenge time-course. However, no upregulation could be recorded in spleen tissue measured by qPCR in juvenile trout (270 g bw). Using mRNA sequencing again, mRNA abundances were determined in gill of juvenile trout (~57.7 g bw) exposed to various aquaculture stressors. The results indicated that all six nkls (nkl1-4 and nkl-like a and nkl-like b) were downregulated when exposed to high temperature, and that nkl1 was significantly downregulated following salinity challenge. Overall, these newly characterized AMPs may contribute to host innate immunity as they are modulated following pathogen challenge and by physiological stressors.

Isolation and characterization of SNP markers of rainbow trout (Onchorhynchus mykiss Walbaum, 1792) from transcriptomic sequences

Rainbow trout (Onchorhynchus mykiss) is one of the most important freshwater aquaculture fish in Iran. It is necessary to develop available molecular marker such as SNPs, which represent a useful tool in detecting adaptive signals in populations and also parentage assignment for O. mykiss. Genetic architecture of broodstock populations is important for breeding programs, as it enables decisions on broodstock screening and genomic selection. In this study, 52 novel single nucleotide polymorphism (SNP) markers for O. mykiss were discovered and validated based on transcriptome sequencing, by means of paired-end sequencing in an Illumina HiSeq 2500 platform. The SNPs were identified through liver transcriptome sequencing from fifteen samples. The observed and expected heterozygosities ranged from 0.177 to 1.000 and 0.239 to 0.638, respectively. The minimum allele frequency (MAF) ranged from 0.166 to 0.489. Among these SNP loci, twenty-two loci showed significant departures from the Hardy-Weinberg equilibrium after Bonferroni correction (p < 0.05) and significant linkage disequilibrium was found. The SNP markers identified in this research could be useful for novel studies, such as those related to associations between high-resolution molecular markers and quantitative traits studies. Moreover, these SNP markers would be used in genetic studies helping economic performance improvement and management of this species.

The Transcription Factor ArcA Modulates Salmonella's Metabolism in Response to Neutrophil Hypochlorous Acid-Mediated Stress

Salmonella Typhimurium, a bacterial pathogen with high metabolic plasticity, can adapt to different environmental conditions; these traits enhance its virulence by enabling bacterial survival. Neutrophils play important roles in the innate immune response, including the production of microbicidal reactive oxygen species (ROS). In addition, the myeloperoxidase in neutrophils catalyzes the formation of hypochlorous acid (HOCl), a highly toxic molecule that reacts with essential biomolecules, causing oxidative damage including lipid peroxidation and protein carbonylation. The bacterial response regulator ArcA regulates adaptive responses to oxygen levels and influences the survival of Salmonella inside phagocytic cells. Here, we demonstrate by whole transcriptomic analyses that ArcA regulates genes related to various metabolic pathways, enabling bacterial survival during HOCl-stress in vitro. Also, inside neutrophils, ArcA controls the transcription of several metabolic pathways by downregulating the expression of genes related to fatty acid degradation, lysine degradation, and arginine, proline, pyruvate, and propanoate metabolism. ArcA also upregulates genes encoding components of the oxidative pathway. These results underscore the importance of ArcA in ATP generation inside the neutrophil phagosome and its participation in bacterial metabolic adaptations during HOCl stress.

Transcriptome profiling of the low-salinity stress responses in the gills of the juvenile Pseudopleuronectes yokohamae

Salinity is an important environmental factor that affects the life cycle of fish, including their growth, development and reproduction. The marbled flounder, Pseudopleuronectes yokohamae, is an important economic resource and serves as a good model to investigate osmoregulation, as it can adapt to a wide range of salinity levels. However, the lack of genomic resources for this species has hampered the understanding of the mechanisms underlying its salinity tolerance. In this study, RNA-Seq analysis was conducted to identify genes related to salinity adaptation and osmotic regulation in the gill tissue of marbled flounder exposed to different concentrations of environmental salinity (6 and 30 ppt). After de novo assembly, 19,265 genes were annotated by the Nr database. A comparison of expression between the two salinity groups revealed 673 differentially expressed genes, of which 180 were upregulated and 493 were downregulated in the low salinity group relative to the high salinity group. The related molecular biological processes were explored from several important perspectives, and potential functions were determined by enrichment analyses, including those of metabolites in ion transportation, energy metabolism and protein synthesis, and immune responses. This study is the first transcriptomic study conducted on marbled flounder, and it revealed many novel sequences for further biological analyses. In addition, the candidate genes identified in the gene expression analysis provided insights into responses to salinity change and molecular mechanisms underlying osmoregulation in the gills of marbled flounder.

Integrated transcriptome analysis of immunological responses in the pearl sac of the triangle sail mussel (Hyriopsis cumingii) after mantle implantation

For pearl culture of bivalve Hyriopsis cumingii, implantation of the sabio may cause nucleus discharge and increased host death rates. We performed a transcriptome analysis of the pearl sac of H. cumingii for 30 days after mantle implantation; 293863 unigenes were obtained, and 27176 unigenes were identified using nr, nt, KO, Swiss-Prot, Pfam, GO, and KOG databases. We detected 4878 differentially expressed genes (DEGs) through pairwise comparisons. We speculated that the physical condition of the recipient mussels returned to normal in about one month; the period was divided into six vital phases (0, 2 h-6 h, 12 h-24 h, 48 h to 7 days, 14 days and 30 days) on the basis of the overall similarities in DEGs. We compared the DEGs between time points and identified key immune-related genes. Our findings provide information on the immunological reactions induced by implantation in pearl mussels.

Label-free quantification of protein expression in the rainbow trout (Oncorhynchus mykiss) in response to short-term exposure to heat stress

Rainbow trout (Oncorhynchus mykiss) are a cold-water salmonid species that is highly susceptible to heat stress. Summer temperature stress is a common issue in trout aquaculture. To better understand the molecular mechanisms of the heat-stress response in the trout, we used label-free quantitative proteome techniques to identify differentially expressed proteins in the livers of rainbow trout exposed to heat stress. We identified 3362 proteins and 152 differentially expressed proteins (p < 0.05; fold-change >2). Of these, 37 were uniquely expressed in the heat-stress group and 35 were uniquely expressed in the control group. In addition, 42 proteins were significantly upregulated (fold-change >2) and 38 proteins were significantly downregulated (fold-change >2). GO (Gene Ontology) analysis indicated that these differentially expressed proteins were primarily expressed in the nucleus, extracellular matrix, and cytoplasm, and were associated with a variety of functions, including protein binding/bridging and enzyme facilitation. KEGG (Kyoto Encyclopedia of Genes and Genomes) analysis of the differentially expressed proteins showed that, during high temperature stress, many biological processes were extensively altered, particularly the estrogen signaling pathway, the complement and coagulation cascades, and the platelet activation pathway. Our study focused on the identification of a systematic approach for the characterization of regulatory networks. Our results provide a framework for further studies of the heat-stress response in fish.

Transcriptomic analysis reveal an efficient osmoregulatory system in Siberian sturgeon Acipenser baeri in response to salinity stress

Sturgeons are euryhaline fish species that have developed specific mechanisms of osmotic and ion regulation to adapt to waters of varying salinity. For the aim to elucidate the osmoregulation strategy behind its high salinity tolerance of sturgeons, the transcriptomes of gills in Siberian sturgeon Acipenser baeri under salinity stress (30 ppt) were sequenced using deep-sequencing platform Illumina/HiSeq-2500 and differential expression genes (DEGs) were identified. A total of 167, 501, 278 clean reads were obtained and 280, 238 unigenes were composed of those clean reads with the mean length of 520nt, and the N50 of 630 bp. Unigenes Sequence alignment was implemented via KEGG, KOG, NT, NR, PFAM, Swiss-Prot, and GO databases. 62, 242 unigenes (22.21%) were annoated in at least one database. 11380 significantly differentially expressed unigenes were found, 6969 of which were up-regulated and 4411 were down-regulated by salinity stress. Amongst the top 20 KEGG pathways with the most amount of annotation sequences, some pathways such as glycerophospholipid metabolism, fatty-acid biosynthesis, glycolysis/gluconeogenesis, oxidative phosphorylation have been comprehensively proved to be relevant to osmoregulation. Despite of these, three possible osmoregulation-related signaling pathways as lipid metabolism related pathways, tight junction pathway and thyroid hormone signaling pathway have been widely analyzed in the current study. In all DEGs, some of the typical genes involved in osmoregulation, including calcium-transporting ATPase 4 (ATP2B4), Na⁺/K⁺-ATPase alpha subunit (α-NKA), potassium-transporting ATPase alpha chain 1 (ATP4A) and Ras GTPase-activating protein (RasGAP) etc were also identified. RNA-seq results were validated with quantitative real-time PCR (qPCR), the 12 selected genes showed a consistent direction in both DGE library and qPCR analysis, proving that the RNA-seq results are reliable. The present results would be helpful to elucidate the osmoregulation mechanism of aquatic animals adapting to salinity challenge.

Population genetic analysis of aquaculture salmonid populations in China using a 57K rainbow trout SNP array

Various salmonid species are cultivated in cold water aquaculture. However, due to limited genomic data resources, specific high-throughput genotyping tools are not available to many of the salmonid species. In this study, a 57K single nucleotide polymorphism (SNP) array for rainbow trout (Oncorhynchus mykiss) was utilized to detect polymorphisms in seven salmonid species, including Hucho taimen, Oncorhynchus masou, Salvelinus fontinalis, Brachymystax lenok, Salvelinus leucomaenis, O. kisutch, and O. mykiss. The number of polymorphic markers per population ranged from 3,844 (O. kisutch) to 53,734 (O. mykiss), indicating that the rainbow trout SNP array was applicable as a universal genotyping tool for other salmonid species. Among the six other salmonid populations from four genera, 28,882 SNPs were shared, whereas 525 SNPs were polymorphic in all four genera. The genetic diversity and population relationships of the seven salmonid species were studied by principal component analysis (PCA). The phylogenetic relationships among populations were analyzed using the maximum likelihood method, which indicated that the shared SNP markers provide reliable genomic information for population genetic analyses in common aquaculture salmonid fishes. Furthermore, this obtained genomic information may be applicable for population genetic evaluation, marker-assisted breeding, and propagative parent selection in fry production.

Pathologic features of mycobacteriosis in naturally infected Syngnathidae and novel transcriptome assembly in association with disease

Syngnathidae (seahorses, seadragons and pipefish) suffer significant losses from non-tuberculous mycobacteria. However, they produce markedly different lesions in response to the disease compared to other teleost species, notably infrequent granuloma formation. This study evaluated 270 syngnathid fish, from which 92 were diagnosed with mycobacteriosis by histopathology, culture or both. Microscopic lesions variably consisted of random foci of coagulative necrosis in multiple organs, containing high numbers of free bacteria and large aggregates or sheets of macrophages with cytoplasm laden with acid-fast bacilli. Mycobacterial associated granulomas were identified in only six seahorses. Five fish had positive cultures with no observed microscopic changes. RNA-seq of the head kidney was performed to investigate the transcriptome of two infected and six non-infected lined seahorses Hippocampus erectus. Assembled and annotated putative transcripts serve to enrich the database for this species, as well as provide baseline data for understanding the pathogenesis of mycobacteriosis in seahorses. Putative components of the innate immune system (IL-1β, IL-6, TNF, NOS, Toll-like receptor 1, MHC Class I, NF-κβ, transforming growth factor beta, MyD88) were identified in the RNA-seq data set. However, a homolog for a key component in the TH1 adaptive immune response, interferon-gamma, was not identified and may underlie the unique pathologic presentation.

RNA-Seq analysis of salinity stress-responsive transcriptome in the liver of spotted sea bass (Lateolabrax maculatus)

Salinity is one of the most prominent abiotic factors, which greatly influence reproduction, development, growth, physiological and metabolic activities of fishes. Spotted sea bass (Lateolabrax maculatus), as a euryhaline marine teleost, has extraordinary ability to deal with a wide range of salinity changes. However, this species is devoid of genomic resources, and no study has been conducted at the transcriptomic level to determine genes responsible for salinity regulation, which impedes the understanding of the fundamental mechanism conferring tolerance to salinity fluctuations. Liver, as the major metabolic organ, is the key source supplying energy for iono- and osmoregulation in fish, however, little attention has been paid to its salinity-related functions but which should not be ignored. In this study, we perform RNA-Seq analysis to identify genes involved in salinity adaptation and osmoregulation in liver of spotted sea bass, generating from the fishes exposed to low and high salinity water (5 vs 30ppt). After de novo assembly, annotation and differential gene expression analysis, a total of 455 genes were differentially expressed, including 184 up-regulated and 271 down-regulated transcripts in low salinity-acclimated fish group compared with that in high salinity-acclimated group. A number of genes with a potential role in salinity adaptation for spotted sea bass were classified into five functional categories based on the gene ontology (GO) and enrichment analysis, which include genes involved in metabolites and ion transporters, energy metabolism, signal transduction, immune response and structure reorganization. The candidate genes identified in L. maculates liver provide valuable information to explore new pathways related to fish salinity and osmotic regulation. Besides, the transcriptomic sequencing data supplies significant resources for identification of novel genes and further studying biological questions in spotted sea bass.

Transcriptional response to heat shock in liver of snow trout (Schizothorax richardsonii)--a vulnerable Himalayan Cyprinid fish

The snow trout (Schizothorax richardsonii) belonging to family Cyprinidae, is an endemic fish of the Himalayan region. The species is tagged as vulnerable species in the IUCN red list of threatened species. The fish thrives well in snowmelt water of several streams and rivers in the region but are occasionally exposed to more than 20 °C during the summer season. Therefore, we have used deep RNA sequencing to decipher the transcriptome of snow trout and characterize the genes and molecular pathways involved in heat shock response. In this study 72,601,298 and 65,428,283 raw reads for heat-shocked and control, respectively, were obtained by Illumina paired-end sequencing technology. The de novo assembled transcriptome was tested for differential gene expression across the treatment groups. The quality of assembly was high with N75 and N50 lengths of 461 and 1274 bases, respectively. A total of 65 unique transcripts were differentially expressed in liver under heat shock and control. Annotated blast matches reveal that differentially expressed transcripts correspond to critical chaperones and molecular pathways, previously shown to be important for thermal stress in other fish species. Eight randomly selected heat-stressed responsive transcripts were also observed to be upregulated during qRT-PCR analysis. This study is the preliminary step to understanding the responses during sudden environmental changes like heat shock. The reference transcriptome database would also aid further studies on biological and physiological aspects of the snow trout under abiotic stresses.

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

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

Transcriptomic response to heat stress among ecologically divergent populations of redband trout

Background

As ectothermic organisms have evolved to differing aquatic climates, the molecular basis of thermal adaptation is a key area of research. In this study, we tested for differential transcriptional response of ecologically divergent populations of redband trout (Oncorhynchus mykiss gairdneri) that have evolved in desert and montane climates. Each pure strain and their F1 cross were reared in a common garden environment and exposed over four weeks to diel water temperatures that were similar to those experienced in desert climates within the species’ range. Gill tissues were collected from the three strains of fish (desert, montane, F1 crosses) at the peak of heat stress and tested for mRNA expression differences across the transcriptome with RNA-seq.

Results

Strong differences in transcriptomic response to heat stress were observed across strains confirming that fish from desert environments have evolved diverse mechanisms to cope with stressful environments. As expected, a large number of total transcripts (12,814) were differentially expressed in the study (FDR ≤ 0.05) with 2310 transcripts in common for all three strains, but the desert strain had a larger number of unique differentially expressed transcripts (2875) than the montane (1982) or the F1 (2355) strain. Strongly differentiated genes (>4 fold change and FDR ≤ 0.05) were particularly abundant in the desert strain (824 unique contigs) relative to the other two strains (montane = 58; F1 = 192).

Conclusions

This study demonstrated patterns of acclimation (i.e., phenotypic plasticity) within strains and evolutionary adaptation among strains in numerous genes throughout the transcriptome. Key stress response genes such as molecular chaperones (i.e., heat shock proteins) had adaptive patterns of gene expression among strains, but also a much higher number of metabolic and cellular process genes were differentially expressed in the desert strain demonstrating these biological pathways are critical for thermal adaptation to warm aquatic climates. The results of this study further elucidate the molecular basis for thermal adaptation in aquatic ecosystems and extend the potential for identifying genes that may be critical for adaptation to changing climates.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1246-5) contains supplementary material, which is available to authorized users.

Transcriptome analysis of northern elephant seal (Mirounga angustirostris) muscle tissue provides a novel molecular resource and physiological insights

Background

The northern elephant seal, Mirounga angustirostris, is a valuable animal model of fasting adaptation and hypoxic stress tolerance. However, no reference sequence is currently available for this and many other marine mammal study systems, hindering molecular understanding of marine adaptations and unique physiology.

Results

We sequenced a transcriptome of M. angustirostris derived from muscle sampled during an acute stress challenge experiment to identify species-specific markers of stress axis activation and recovery. De novo assembly generated 164,966 contigs and a total of 522,699 transcripts, of which 68.70% were annotated using mouse, human, and domestic dog reference protein sequences. To reduce transcript redundancy, we removed highly similar isoforms in large gene families and produced a filtered assembly containing 336,657 transcripts. We found that a large number of annotated genes are associated with metabolic signaling, immune and stress responses, and muscle function. Preliminary differential expression analysis suggests a limited transcriptional response to acute stress involving alterations in metabolic and immune pathways and muscle tissue maintenance, potentially driven by early response transcription factors such as Cebpd.

Conclusions

We present the first reference sequence for Mirounga angustirostris produced by RNA sequencing of muscle tissue and cloud-based de novo transcriptome assembly. We annotated 395,102 transcripts, some of which may be novel isoforms, and have identified thousands of genes involved in key physiological processes. This resource provides elephant seal-specific gene sequences, complementing existing metabolite and protein expression studies and enabling future work on molecular pathways regulating adaptations such as fasting, hypoxia, and environmental stress responses in marine mammals.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1253-6) contains supplementary material, which is available to authorized users.

Digital gene expression analysis in hemocytes of the white shrimp Litopenaeus vannamei in response to low salinity stress

The white shrimp Litopenaeus vannamei has been greatly impacted by low salinity stress. To gain knowledge on the immune response in L. vannamei under such stress, we investigated digital gene expression (DEG) in L. vannamei hemocytes using the deep-sequencing platform Illumina HiSeq 2000. In total, 38,155 high quality unigenes with average length 770 bp were generated; 145 and 79 genes were identified up- or down-regulated, respectively. Functional categorization and pathways of the differentially expressed genes revealed that immune signaling pathways, cellular immunity, humoral immunity, apoptosis, cellular protein synthesis, lipid transport and energy metabolism were the differentially regulated processes occurring during low salinity stress. These results will provide a resource for subsequent gene expression studies regarding environmental stress and a valuable gene information for a better understanding of immune mechanisms of L. vannamei under low salinity stress.

Transcriptome and expression profiling analysis of the hemocytes reveals a large number of immune-related genes in mud crab Scylla paramamosain during Vibrio parahaemolyticus infection

BACKGROUND

Mud crab Scylla paramamosain is an economically important marine species in China. However, frequent outbreaks of infectious diseases caused by marine bacteria, such as Vibrio parahaemolyticus, result in great economic losses.

METHODOLOGY/PRINCIPAL FINDINGS

Comparative de novo transcriptome analysis of S. paramamosain infected with V. parahaemolyticus was carried out to investigate the molecular mechanisms underlying the immune response to pathogenic bacteria by using the Illumina paired-end sequencing platform. A total of 52,934,042 clean reads from the hemocytes of V. parahaemolyticus-infected mud crabs and controls were obtained and assembled into 186,193 contigs. 59,120 unigenes were identified from 81,709 consensus sequences of mud crabs and 48,934 unigenes were matched proteins in the Nr or Swissprot databases. Among these, 10,566 unigenes belong to 3 categories of Gene Ontology, 25,349 to 30 categories of KEGG, and 15,191 to 25 categories of COG database, covering almost all functional categories. By using the Solexa/Illumina's DGE platform, 1213 differentially expressed genes (P<0.05), including 538 significantly up-regulated and 675 down-regulated, were detected in V. parahaemolyticus-infected crabs as compared to that in the controls. Transcript levels of randomly-chosen genes were further measured by quantitative real-time PCR to confirm the expression profiles. Many differentially expressed genes are involved in various immune processes, including stimulation of the Toll pathway, Immune Deficiency (IMD) pathway, Ras-regulated endocytosis, and proPO-activating system.

CONCLUSIONS/SIGNIFICANCE

Analysis of the expression profile of crabs under infection provides invaluable new data for biological research in S. paramamosain, such as the identification of novel genes in the hemocytes during V. parahaemolyticus infection. These results will facilitate our comprehensive understanding of the mechanisms involved in the immune response to bacterial infection and will be helpful for diseases prevention in crab aquaculture.

The development and characterization of a 57K single nucleotide polymorphism array for rainbow trout

In this study, we describe the development and characterization of the first high-density single nucleotide polymorphism (SNP) genotyping array for rainbow trout. The SNP array is publically available from a commercial vendor (Affymetrix). The SNP genotyping quality was high, and validation rate was close to 90%. This is comparable to other farm animals and is much higher than previous smaller scale SNP validation studies in rainbow trout. High quality and integrity of the genotypes are evident from sample reproducibility and from nearly 100% agreement in genotyping results from other methods. The array is very useful for rainbow trout aquaculture populations with more than 40 900 polymorphic markers per population. For wild populations that were confounded by a smaller sample size, the number of polymorphic markers was between 10 577 and 24 330. Comparison between genotypes from individual populations suggests good potential for identifying candidate markers for populations' traceability. Linkage analysis and mapping of the SNPs to the reference genome assembly provide strong evidence for a wide distribution throughout the genome with good representation in all 29 chromosomes. A total of 68% of the genome scaffolds and contigs were anchored through linkage analysis using the SNP array genotypes, including ~20% of the genome assembly that has not been previously anchored to chromosomes.

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

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

Qualitative de novo analysis of full length cDNA and quantitative analysis of gene expression for common marmoset (Callithrix jacchus) transcriptomes using parallel long-read technology and short-read sequencing

The common marmoset (Callithrix jacchus) is a non-human primate that could prove useful as human pharmacokinetic and biomedical research models. The cytochromes P450 (P450s) are a superfamily of enzymes that have critical roles in drug metabolism and disposition via monooxygenation of a broad range of xenobiotics; however, information on some marmoset P450s is currently limited. Therefore, identification and quantitative analysis of tissue-specific mRNA transcripts, including those of P450s and flavin-containing monooxygenases (FMO, another monooxygenase family), need to be carried out in detail before the marmoset can be used as an animal model in drug development. De novo assembly and expression analysis of marmoset transcripts were conducted with pooled liver, intestine, kidney, and brain samples from three male and three female marmosets. After unique sequences were automatically aligned by assembling software, the mean contig length was 718 bp (with a standard deviation of 457 bp) among a total of 47,883 transcripts. Approximately 30% of the total transcripts were matched to known marmoset sequences. Gene expression in 18 marmoset P450- and 4 FMO-like genes displayed some tissue-specific patterns. Of these, the three most highly expressed in marmoset liver were P450 2D-, 2E-, and 3A-like genes. In extrahepatic tissues, including brain, gene expressions of these monooxygenases were lower than those in liver, although P450 3A4 (previously P450 3A21) in intestine and P450 4A11- and FMO1-like genes in kidney were relatively highly expressed. By means of massive parallel long-read sequencing and short-read technology applied to marmoset liver, intestine, kidney, and brain, the combined next-generation sequencing analyses reported here were able to identify novel marmoset drug-metabolizing P450 transcripts that have until now been little reported. These results provide a foundation for mechanistic studies and pave the way for the use of marmosets as model animals for drug development in the future.

Transcriptome analysis of the Portunus trituberculatus: de novo assembly, growth-related gene identification and marker discovery

Background

The swimming crab, Portunus trituberculatus, is an important farmed species in China, has been attracting extensive studies, which require more and more genome background knowledge. To date, the sequencing of its whole genome is unavailable and transcriptomic information is also scarce for this species. In the present study, we performed de novo transcriptome sequencing to produce a comprehensive transcript dataset for major tissues of Portunus trituberculatus by the Illumina paired-end sequencing technology.

Results

Total RNA was isolated from eyestalk, gill, heart, hepatopancreas and muscle. Equal quantities of RNA from each tissue were pooled to construct a cDNA library. Using the Illumina paired-end sequencing technology, we generated a total of 120,137 transcripts with an average length of 1037 bp. Further assembly analysis showed that all contigs contributed to 87,100 unigenes, of these, 16,029 unigenes (18.40% of the total) can be matched in the GenBank non-redundant database. Potential genes and their functions were predicted by GO, KEGG pathway mapping and COG analysis. Based on our sequence analysis and published literature, many putative genes with fundamental roles in growth and muscle development, including actin, myosin, tropomyosin, troponin and other potentially important candidate genes were identified for the first time in this specie. Furthermore, 22,673 SSRs and 66,191 high-confidence SNPs were identified in this EST dataset.

Conclusion

The transcriptome provides an invaluable new data for a functional genomics resource and future biological research in Portunus trituberculatus. The data will also instruct future functional studies to manipulate or select for genes influencing growth that should find practical applications in aquaculture breeding programs. The molecular markers identified in this study will provide a material basis for future genetic linkage and quantitative trait loci analyses, and will be essential for accelerating aquaculture breeding programs with this species.

RNA-seq analysis of early hepatic response to handling and confinement stress in rainbow trout

Fish under intensive rearing conditions experience various stressors which have negative impacts on survival, growth, reproduction and fillet quality. Identifying and characterizing the molecular mechanisms underlying stress responses will facilitate the development of strategies that aim to improve animal welfare and aquaculture production efficiency. In this study, we used RNA-seq to identify transcripts which are differentially expressed in the rainbow trout liver in response to handling and confinement stress. These stressors were selected due to their relevance in aquaculture production. Total RNA was extracted from the livers of individual fish in five tanks having eight fish each, including three tanks of fish subjected to a 3 hour handling and confinement stress and two control tanks. Equal amount of total RNA of six individual fish was pooled by tank to create five RNA-seq libraries which were sequenced in one lane of Illumina HiSeq 2000. Three sequencing runs were conducted to obtain a total of 491,570,566 reads which were mapped onto the previously generated stress reference transcriptome to identify 316 differentially expressed transcripts (DETs). Twenty one DETs were selected for qPCR to validate the RNA-seq approach. The fold changes in gene expression identified by RNA-seq and qPCR were highly correlated (R² = 0.88). Several gene ontology terms including transcription factor activity and biological process such as glucose metabolic process were enriched among these DETs. Pathways involved in response to handling and confinement stress were implicated by mapping the DETs to reference pathways in the KEGG database.

Accession Numbers

Raw RNA-seq reads have been submitted to the NCBI Short Read Archive under accession number SRP022881.

Customized Perl Scripts

All customized scripts described in this paper are available from Dr. Guangtu Gao or the corresponding author.

Transcriptome analysis of Portunus trituberculatus in response to salinity stress provides insights into the molecular basis of osmoregulation

Background

The swimming crab, Portunus trituberculatus, which is naturally distributed in the coastal waters of Asia-Pacific countries, is an important farmed species in China. Salinity is one of the most important abiotic factors that influence not only the distribution and abundance of crustaceans, it is also an important factor for artificial propagation of the crab. To better understand the interaction between salinity stress and osmoregulation, we performed a transcriptome analysis in the gills of Portunus trituberculatus challenged with salinity stress, using the Illumina Deep Sequencing technology.

Results

We obtained 27,696,835, 28,268,353 and 33,901,271 qualified Illumina read pairs from low salinity challenged (LC), non-challenged (NC), and high salinity challenged (HC) Portunus trituberculatus cDNA libraries, respectively. The overall de novo assembly of cDNA sequence data generated 94,511 unigenes, with an average length of 644 bp. Comparative genomic analysis revealed that 1,705 genes differentially expressed in salinity stress compared to the controls, including 615 and 1,516 unigenes in NC vs LC and NC vs HC respectively. GO functional enrichment analysis results showed some differentially expressed genes were involved in crucial processes related to osmoregulation, such as ion transport processes, amino acid metabolism and synthesis processes, proteolysis process and chitin metabolic process.

Conclusion

This work represents the first report of the utilization of the next generation sequencing techniques for transcriptome analysis in Portunus trituberculatus and provides valuable information on salinity adaptation mechanism. Results reveal a substantial number of genes modified by salinity stress and a few important salinity acclimation pathways, which will serve as an invaluable resource for revealing the molecular basis of osmoregulation in Portunus trituberculatus. In addition, the most comprehensive sequences of transcripts reported in this study provide a rich source for identification of novel genes in the crab.

Transcriptome of Atlantic cod (Gadus morhua L.) early embryos from farmed and wild broodstocks

Significant efforts have been made to elucidate factors affecting egg quality in fish. Recently, we have shown that eggs originating from wild broodstock (WB) of Atlantic cod (Gadus morhua L.) are of superior quality to those derived from farmed broodstock (FB), and this is associated with differences in the chemical composition of egg yolk. However, maternal transcripts, accumulated during oogenesis, have not been studied extensively in fish. The aim of the present study was to characterize putative maternal mRNA transcriptome in fertilized eggs of Atlantic cod and to compare transcript pools between WB and FB in order to investigate the relation between egg developmental potential and putative maternal mRNA deposits. We performed high-throughput 454 pyrosequencing. For each WB and FB group, five cDNA libraries were individually tagged and sequenced, resulting in 98,687 (WB) and 119,333 (FB) average reads per library. Sequencing reads were de novo assembled, annotated, and mapped. Out of 13,726 identified isotigs, 238 were differentially expressed between WB and FB, with 155 isotigs significantly upregulated in WB. The sequence reads were mapped to 11,340 different Atlantic cod transcripts and 158 sequences were differentially expressed between the 2 groups. Important transcripts involved in fructose metabolism, fatty acid metabolism, glycerophospholipid metabolism, and oxidative phosphorylation were differentially represented between the two broodstock groups, showing potential as biomarkers of egg quality in teleosts. Our findings contribute to the hypothesis that maternal mRNAs affect egg quality and, consequently, the early development of fish.

Basal polarization of the mucosal compartment in Flavobacterium columnare susceptible and resistant channel catfish (Ictalurus punctatus)

The freshwater bacterial pathogen, Flavobacterium columnare, infects a variety of ornamental and farmed fish species worldwide through mucosal attachment points on the gill and skin. While previous studies have demonstrated a chemotactic response of F. columnare to fish mucus, little is known about how host gill mucosal molecular and cellular constituents may impact rates of adhesion, tissue invasion, and ultimately, mortality. Here, we describe the use of RNA-seq to profile gill expression differences between channel catfish (Ictalurus punctatus) differing in their susceptibility to F. columnare both basally (before infection) and at three early timepoints post-infection (1 h, 2 h, and 8 h). After sequencing and de novo assembly of over 350 million 100 base-pair transcript reads, between group comparisons revealed 1714 unique genes differentially expressed greater than 1.5-fold at one or more timepoints. In the large dataset, we focused our analysis on basal differential expression between resistant and susceptible catfish as these genes could potentially reveal genetic and/or environmental factors linked with differential rates of infection. A number of critical innate immune components including iNOS2b, lysozyme C, IL-8, and TNF-alpha were constitutively higher in resistant catfish gill, while susceptible fish showed high expression levels of secreted mucin forms, a rhamnose-binding lectin previously linked to susceptibility, and mucosal immune factors such as CD103 and IL-17. Taken together, the immune and mucin profiles obtained by RNA-seq suggest a basal polarization in the gill mucosa, with susceptible fish possessing a putative mucosecretory, toleragenic phenotype which may predispose them to F. columnare infection.

Sequencing of the first ayu (Plecoglossus altivelis) macrophage transcriptome and microarray development for investigation the effect of LECT2 on macrophages

Macrophages play an important role in first-line host defense of innate immune in fishes. However, it is difficult to investigate cellular mechanism of immune response in fish species with little genomic information available. Here we present the first use of RNA-Sequencing to study the macrophage transcriptome of ayu, Plecoglossus altivelis, which is an economically important fish in East Asia. De novo assembly generated 49,808 non-redundant consensus sequences, among which 23,490 transcripts found respective coding sequences. 15,707 transcripts are predicted to be involved in known metabolic or signaling pathways. The sequences were then used to develop a microarray for measurement the effect of recombinant LECT2 on ayu macrophages. LECT2 altered expression of a variety of genes mainly implicated in actin cytoskeleton, pattern recognition receptors and cytokines. Meanwhile, LECT2 enhanced phagocytosis, bacterial killing, and respiratory burst in ayu macrophages, which supported the thought derived from the microarray data that LECT2 activates macrophages. In conclusion, our results contribute to understanding the specific regulation mechanism of LECT2 in macrophage activation, and the combination of transcriptome analysis and microarray assay is a good method for screening a special tissue or cell response to a stimulus or pathogen in non-model fish species.

Stocking impacts the expression of candidate genes and physiological condition in introgressed brook charr (Salvelinus fontinalis) populations

Translocation of plants and animal populations between environments is one of the major forms of anthropogenic perturbation experienced by pristine populations, and consequently, human-mediated hybridization by stocking practices between wild and exogenous conspecifics is of increasing concern. In this study, we compared the expression of seven candidate genes involved in multifactorial traits and regulatory pathways for growth as a function of level of introgressive hybridization between wild and domestic brook charr to test the null hypothesis of no effect of introgression on wild fish. Our analyses revealed that the expression of two of the genes tested, cytochrome c oxidase VIIa and the growth hormone receptor isoform I, was positively correlated with the level of introgression. We also observed a positive relationship between the extent of introgression and physiological status quantified by the Fulton's condition index. The expression of other genes was influenced by other variables, including year of sampling (reflecting different thermal conditions), sampling method and lake of origin. This is the first demonstration in nature that introgression from stocked populations has an impact on the expression of genes playing a role in important biological functions that may be related with fitness in wild introgressed populations.

A moderate increase in ambient temperature modulates the Atlantic cod (Gadus morhua) spleen transcriptome response to intraperitoneal viral mimic injection

BACKGROUND

Atlantic cod (Gadus morhua) reared in sea-cages can experience large variations in temperature, and these have been shown to affect their immune function. We used the new 20K Atlantic cod microarray to investigate how a water temperature change which, simulates that seen in Newfoundland during the spring-summer (i.e. from 10°C to 16°C, 1°C increase every 5 days) impacted the cod spleen transcriptome response to the intraperitoneal injection of a viral mimic (polyriboinosinic polyribocytidylic acid, pIC).

RESULTS

The temperature regime alone did not cause any significant increases in plasma cortisol levels and only minor changes in spleen gene transcription. However, it had a considerable impact on the fish spleen transcriptome response to pIC [290 and 339 significantly differentially expressed genes between 16°C and 10°C at 6 and 24 hours post-injection (HPI), respectively]. Seventeen microarray-identified transcripts were selected for QPCR validation based on immune-relevant functional annotations. Fifteen of these transcripts (i.e. 88%), including DHX58, STAT1, IRF7, ISG15, RSAD2 and IκBα, were shown by QPCR to be significantly induced by pIC.

CONCLUSIONS

The temperature increase appeared to accelerate the spleen immune transcriptome response to pIC. We found 41 and 999 genes differentially expressed between fish injected with PBS vs. pIC at 10°C and sampled at 6HPI and 24HPI, respectively. In contrast, there were 656 and 246 genes differentially expressed between fish injected with PBS vs. pIC at 16°C and sampled at 6HPI and 24HPI, respectively. Our results indicate that the modulation of mRNA expression of genes belonging to the NF-κB and type I interferon signal transduction pathways may play a role in controlling temperature-induced changes in the spleen's transcript expression response to pIC. Moreover, interferon effector genes such as ISG15 and RSAD2 were differentially expressed between fish injected with pIC at 10°C vs. 16°C at 6HPI. These results substantially increase our understanding of the genes and molecular pathways involved in the negative impacts of elevated ambient temperature on fish health, and may also be valuable to our understanding of how accelerated global climate change could impact cold-water marine finfish species.

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

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.