RNA-seq analysis reveals extensive transcriptional plasticity to temperature stress in a freshwater fish species
© Smith et al.; licensee BioMed Central Ltd. 2013
Received: 29 January 2013
Accepted: 27 May 2013
Published: 5 June 2013
Identifying genes of adaptive significance in a changing environment is a major focus of ecological genomics. Such efforts were restricted, until recently, to researchers studying a small group of model organisms or closely related taxa. With the advent of next generation sequencing (NGS), genomes and transcriptomes of virtually any species are now available for studies of adaptive evolution. We experimentally manipulated temperature conditions for two groups of crimson spotted rainbowfish (Melanotaenia duboulayi) and measured differences in RNA transcription between them. This non-migratory species is found across a latitudinal thermal gradient in eastern Australia and is predicted to be negatively impacted by ongoing environmental and climatic change.
Using next generation RNA-seq technologies on an Illumina HiSeq2000 platform, we assembled a de novo transcriptome and tested for differential expression across the treatment groups. Quality of the assembly was high with a N50 length of 1856 bases. Of the 107,749 assembled contigs, we identified 4251 that were differentially expressed according to a consensus of four different mapping and significance testing approaches. Once duplicate isoforms were removed, we were able to annotate 614 up-regulated transfrags and 349 that showed reduced expression in the higher temperature group.
Annotated blast matches reveal that differentially expressed genes correspond to critical metabolic pathways previously shown to be important for temperature tolerance in other fish species. Our results indicate that rainbowfish exhibit predictable plastic regulatory responses to temperature stress and the genes we identified provide excellent candidates for further investigations of population adaptation to increasing temperatures.
KeywordsRainbowfish Melanotaenia duboulayi Transcriptomes Climate change Thermal adaptation
The ability of species and populations to adapt to environmental change is the cornerstone of the emerging field of ecological genomics [1, 2]. Until recently, genome-wide studies of genetic adaptation in non-model organisms were not possible. With the advent of massively parallel next generation sequencing technologies (NGS), these types of studies have become a reality  and while many of the challenges and preferred strategies are still being addressed [4–6], empirical studies are now starting to be reported [7–14]. Studies of transcriptome level responses to environmental change offer an opportunity to understand the underlying genetic basis for adaptation. Such studies represent a powerful approach to assessing the genes involved in adaptation to a changing climate, particularly increasing temperatures. By profiling transcriptional changes induced by temperature stress, it is possible to identify the gene regions or pathways that are likely to be the targets of selection. This information is crucial to enable researchers to assess levels of variation across these gene regions, at a landscape scale, to predict the capacity of organisms to adapt to a warming climate.
Genes involved in physiological adaptation to temperature stress have been uncovered in many species. Heat shock proteins , alcohol dehydrogenase  and lactate dehydrogenase genes  have all been shown to be related to heat tolerance. In fish, the list of candidates also includes many from other gene regions related to respiration and protein binding [18–20]. Apart from differences in coding regions, transcriptional regulation is also a source of variation that can potentially contribute to adaptive evolutionary change, particularly in the early stages of divergence. Studies in natural populations of gobies (Gillichthys mirabilis) have shown that short term exposure (<8 hours) to a temperature of 32°C induces a strong upregulation of heat shock proteins (Hsps) in both gill and muscle tissues . Many other transcripts related to a wide variety of biological processes including protein homeostasis, cell cycle control, cytoskeletal reorganisations, metabolic regulation, and signal transduction were differentially expressed in treatment and control groups. The majority of these genes displayed tissue-specific responses presumably related to the differing molecular functions associated with each tissue type. Logan and Somero  found that, with long-term acclimation to increased temperature (up to 28°C), there was no upregulation of stress-related proteins and only slight, although detectable, differences in expression of genes involved in protein biosynthesis, transport and various metabolic categories. This they suggest indicates evidence of long-term acclimation showing a steady state condition involving relative energy costs for different processes. They later showed however, that stress related genes (HSP70, UBIQ, and CDKN1B) were induced in long-term acclimatised fish subsequently exposed to acute heating conditions (4°C/hour) and that the onset temperature for significant expression change varied according to acclimation temperature . Quinn et al. also found increased expression of HSPs and Ubiquitin in Arctic charr (Salvelinus alpinus) exposed to temperature stress and reported a down regulation of haemoglobin genes in fish that showed tolerance to increased temperatures. Another cold climate fish, Trematomus bernacchii, has been shown to be unable to mount a heat shock response despite retaining the heat shock gene Hsp70 and the regulation factor HSF1 . Further work showed that many other genes associated with the cellular stress response were induced by heat stress. The inability to mount a heat shock response however, highlights the susceptibility of this species to global warming and raises the question as to how this and other species will be able to adapt to increasing temperatures.
Buckley and Hofmann  examined the extensive plasticity in Hsp induction in gobies acclimatised to different thermal backgrounds (13°C, 21°C, and 28°C). They found that the activation temperature of the transcriptional regulator HSF1 was positively associated with the acclimatisation temperature indicating that plasticity in heat shock response is linked to plasticity in the regulatory framework governing Hsps. While adaptive plasticity is often seen as a mechanism that can slow or dampen divergent selection, it has been argued that it can also lead to rapid speciation if there are strong correlations between phenotype and environment combined with significant population structure . By examining the transcriptomic response to temperature stress we can develop a better understanding of the genes and biochemical pathways that are fundamental to physiological acclimatisation to a warming environment and gain insights into the regulatory changes that accompany adaptation over evolutionary timescales .
Australian rainbowfish are an ideal species group to test hypotheses about the genetic responses to increasing temperatures. In particular, the crimson-spotted rainbowfish (Melanotaenia duboulayi) is a subtropical freshwater species found along a north–south temperature gradient in eastern Australia. Their distribution ranges over several ecoregions which, coupled with a strong population structure and local abundance [29–31], makes them a well suited model for studying local adaptation. The ease of maintaining captive populations of rainbowfish also make them amenable to a range of laboratory-based experimental studies [32–34]. In this study, we maintained groups of M. duboulayi at ambient and elevated temperature levels and then used an RNA-seq approach to assess transcriptome level changes related to temperature stress. Our aim is to provide an initial investigation of the transcriptomic response to thermal stress in rainbowfish. As such, this will allow for the screening of many more individuals via genotyping of candidate SNPs. In addition we present the first annotated transcriptome and gene catalogue for the order Atheriniformes. Our goal is to identify key candidate genes and make a first step towards understanding the important biochemical pathways on which selection is likely to act in a warming climate.
Source of fish and design of temperature trial
Crimson spotted rainbowfish were collected using a hand-net from a location in the upper reaches of the Brisbane River, near the township of Fernvale (27°26'37.39"S, 152°40'12.76"E). Water monitoring data from the Queensland Department of Environment and Resource Monitoring (DERM) show the average daily mean temperatures for this location ranged between 12.2°C in winter and 28.3°C in summer from January 1st 2004 to January 1st 2011 (http://watermonitoring.derm.qld.gov.au). Fish were transported live to Flinders University animal rearing facility and acclimatised at a temperature of 21°C for a period of 30 days prior to the start of the temperature trials. For the trials we used only adult male rainbowfish of about the same length (a proxy for age), since gender and age can affect expression responses . These individuals were randomly assigned to a treatment or a control group (n = 6 per group). Temperature in the treatment group was increased by 2°C per day over a period of six days towards a target of 33°C. This target represents the projected average summer temperature for this region in 2070 based on a high emission scenario of the International Panel on Climate Change: http://www.climatechangeinaustralia.gov.au/qldtemp15.php. This temperature condition was then maintained for 14 days. The control group was kept at 21°C for the duration of the experiment. All animal handling was performed in accordance with the Australian Code of Practice for the Care and Use of Animals for Scientific Purposes, 2004 and approved by the Flinders University Animal Welfare Committee (AWC E342).
RNA extraction, Illumina library preparation and sequencing
Upon completion of the temperature trial, fish were sacrificed using AQUI-S® solution  and dissected immediately to remove their livers. Although increased temperature has been shown to differentially induce expression changes in different tissue types [21, 37], we were restricted to examining just one tissue type due to logistical constraints. We selected liver due to previous research linking this tissue type to heat stress responses [38–40]. Total RNAs were individually extracted using the Ambion Magmax™-96 total RNA isolation kit (Life Sciences) according to the manufacturer’s instructions. Briefly, 5 mg of tissue was placed in the lysis solution and homogenised in Qiagen Tissuelyzer™ for a period of 30 sec. Nucleic acids were captured onto magnetic beads, washed and treated with DNase. Total RNA was then eluted in 50 μl elution buffer. RNA quality and concentration was measured using an RNA Pico chip on an Agilent Bioanalyzer. Normalised starting quantities of total RNA were then used to prepare 12 separate Illumina sequencing libraries with the TruSeq™ RNA sample preparation kit (Illumina). Library preparation was performed as per the manufacturer’s instructions. In the final step before sequencing, all 12 individual libraries were normalised and pooled together using the adapter indices supplied by the manufacturer (Illumina MID tags 2, 4–7, 12–16, 18, 19). Pooled sequencing was then performed as 101 bp, paired-end reads in a single lane of an Illumina HiSeq2000 instrument housed at the Ramaciotti Centre for Gene Function Analysis, University of New South Wales.
Quality control and de novo assembly
Sequence data were sorted by individual and adapters were trimmed by the service provider prior to analysis. Quality filtering was performed using the FastX-toolkit suite of pre-processing tools (http://hannonlab.cshl.edu/fastx_toolkit/index.html) in a Galaxy setting . Based on the FastX quality statistics, the first two and last 5 bases were trimmed from each read as they had consistently low phred scores (<Q15). Paired reads were then joined and a quality filter applied such that any combined reads having <90% of bases with a phred score of Q20 or higher were discarded. Paired reads were then split and interleaved to suit the input style of the de novo assembly program. Transcriptome assembly was performed de novo with the program Velvet/Oases . This program reconstructs independent assemblies based on different k-mer values used to build a de Bruijn graph. The program then uses dynamic error removal adapted to RNA-seq data and implements a robust scaffolding method to predict full length transfrags. Multiple single k-mer assemblies are then merged to cover genes at different expression levels without redundancy. Two individuals from each of the treatment and control groups were pooled as input for the assembly. Assemblies were compiled for a k-mer range of 19 to 49 with an expected insert size between paired ends of 300 bp and a coverage cut-off value set to 4.2. We tested different merged assembly ranges based on the summary statistics for each individual k-mer assembly . The outcome of each merge was assessed with respect to the optimal assembly parameters . The optimal assembly should achieve the best balance between large median, mean and N50 contig lengths while minimising the total number of contigs but maintaining a large summed contig length. As Oases is vulnerable to mis-assembly at low k-mer values, we adopted a conservative approach of merging k-mer values > k = 19. Optimal assembly was achieved with a k-mer range of 19 to 41.
Mapping of sequence reads and differential expression analysis
To test for differential expression (DE), individual sequence reads for each sample were mapped back to the assembled transcriptome with the alignment program Bowtie . Bowtie was implemented in the –v alignment mode with the maximum number of mismatches set to 3. Paired end reads were aligned to the transcriptome with both read pairs needing a valid alignment within a given locus to be counted as a match. If more than one alignment was possible the best match was reported according to the least number of mismatches for each read and overall for the pair. The reproducibility of the alignment approach was tested by performing the mapping step with BWA, an alternative alignment program . The number of reads aligning to each transfrag for each sample was calculated with the IdxStats command of Samtools . Count data was then used as input for the program DESeq  which estimates variance-mean dependence in the data and tests for differential expression based on the negative binomial distribution. The six samples from each treatment were used to generate mean expression levels with associated variances. Differential expression was tested at a significance level of α= 0.05 adjusted to match a 5% false discovery rate using the Benjamini-Hochberg procedure. The threshold for fold-change differences is determined by the significance testing as the power to detect significant differential expression depends on the expression strength. For weakly expressed genes, stronger changes are required for the gene to be called significantly expressed. We also compared DE methodology by running the EdgeR program to assess significant differences in the count data. A consensus list of DE genes was then generated from the four analysis approaches adopted (i.e. Bowtie-DESeq, Bowtie-EdgeR, BWA-DESeq, BWA-EdgeR). Significantly up and down regulated transfrags were selected and blasted against the NCBI database using blastx in the program Blast2GO . Blastx was performed against the NCBI nucleotide database with the minimum E-value score set to 1.0E-06. To assign gene ontology terms to each annotated sequence, successful blast hits were mapped and annotated using Blast2GO for the entire assembled transcriptome with the annotation cut-off threshold set to 55 and the GO level weighting set to 5.
Results and discussion
Raw sequencing data and quality statistics
Assembly statistics for k-mer lengths 19–49 and different k-mer merge ranges from the Oases de novo assembly program
Min sequence length
Max sequence length
Average sequence length
Median sequence length
(A + T)s
(G + C)s
Differential expression analyses
The Blast2GO program was able to find sequence similarities for 2,740 of the DE transfrags but could not find mapping or annotation information for a further 634 of them, leaving 2,106 DE transfrags which were successfully annotated. The top 15 matching species from the BLAST query were all fish species with the most BLAST hits being for the Nile tilapia Oreochromis niloticus with 583 matches. Duplicate gene isoforms were detected by matching identical annotated gene names from the Blast2GO output. These isoforms were then combined and reported as single “genes”. Once isoforms were combined, there were 614 genes that were up-regulated in the high temperature treatment with 349 genes being down-regulated (see Additinal file 1: Table S2a and b). For significantly down-regulated transfrags, the mean fold-change between ambient and high-temperature conditions was 4.0-fold, with a range from 55.6-fold for g2/m phase specific e3 ubiquitin-protein ligase to 2.2-fold for the Phytanoyl-peroxisomal-like protein. The mean fold-change for significantly up-regulated transfrags was 11.13, ranging from 1.98 (for the cyclin-dependent kinase 2 interacting protein) to 259-fold (for the Heat shock protein Hsp-90-like).
Ontology of differentially expressed genes
While the Hsp genes are commonly identified as overexpressed in short-term temperature manipulation experiments [24, 37], they are less likely to be targets for selection during gradual temperature shifts associated with climate change [22, 53]. Hsp genes represent a physiological response to sudden stressors and therefore plasticity in these traits is unlikely to be adaptive over longer timescales . The more likely candidates for an adaptive genetic response are those genes involved in what has been termed the “cellular homeostatis response” to long-term temperature stress . Unlike stress response genes that provide an immediate early response to macromolecular damage and sudden changes in cellular redox potential, the cellular homeostatasis response involves effector proteins mediating parameter specific adaptation to environmental change.
Responses associated with prolonged exposure to heat stress
Prolonged exposure to increased temperatures has previously been associated with gene ontologies related to protein folding, oxidative stress and immune function [18, 19]. Similarly, we detected significant upregulation of genes with these ontologies in the high temperature treatment such as Calnexin (2.8 x), NADH dehydrogenase (2.5 x), and glutathione S-transferase (5.1 x) suggesting long-term reallocation of energy resources. Plasticity in the expression of these genes is more likely to be adaptive and allow localised populations to survive in a changing environment, eventually leading to divergent selection. Kassahn et al. grouped stress-response transcripts into four different clusters according to the pattern of regulation detected under short versus long-term exposure to heat stress. They suggested that long-term exposure to heat stress in a coral reef fish (31°C for five days) induces expression of genes involved in development and immune function whereas genes related to metabolic function are suppressed. Our data, from long-term exposure to heat stress in rainbowfish (33°C for 14 days), support those findings. Developmental processes and metabolic processes accounted for 48% of dysregulated transfrags (Figure 3a). Immune function seems less important in our dataset and is covered by the “response to stimuli” category representing 9% of DE transfrags including the natural killer cell enhancement factor (upregulated 2.8 x). It is possible that the longer exposure to heat stress in our study allowed recovery from the immediate activation of the immune function genes.
Candidate genes for broad scale studies of temperature response in the crimson spotted rainbowfish, Melanotaenia duboulayi
min. e value
Primary gene Ontology
C:endoplasmic reticulum membrane
P:isoprenoid biosynthetic process
smooth muscle cell-specific protein sm22 alpha
P:muscle organ development
ap-2 complex subunit sigma
acyl carrier mitochondrial precursor
nadh-cytochrome b5 reductase 2
F:cytochrome-b5 reductase activity
sterol-4-alpha-carboxylate 3- decarboxylating-like
P:steroid biosynthetic process
retinol-binding protein cellular
ribosomal rna processing protein 36 homolog
protein cdv3 homolog
nadh dehydrogenase 1 alpha subcomplex subunit 6
C:mitochondrial inner membrane
y chain e2~ubiquitin-hect
P:catecholamine metabolic process
small nuclear ribonucleoprotein sm d2
P:ncRNA metabolic process
rho-class glutathione s-transferase
thioredoxin domain-containing protein 14 precursor
P:cell redox homeostasis
ubiquitin-conjugating enzyme e2 variant 2
F:acid-amino acid ligase activity
coiled-coil domain-containing protein 47 precursor
ubiquinol-cytochrome c reductase core protein ii
3-hydroxy-3-methylglutaryl-coenzyme a synthase 1
P:response to tellurium ion
nuclear factor erythroid derived 2-like 1
P:heme biosynthetic process
F:glutathione transferase activity
cyclin-dependent kinase inhibitor 1
P:cellular response to stimulus
catechol-o-methyltransferase domain-containing protein 1
period homolog 3
P:respiratory gaseous exchange
scinderin like a
5-aminolevulinate erythroid- mitochondrial-like
P:response to hypoxia
55 kda erythrocyte membrane protein
C:intracellular non-membrane-bounded organelle
cbp p300-interacting transactivator 3b
actin-binding lim protein 1 long isoform isoform cra_a
udp-glucuronosyltransferase 2a2-like isoform 2
F:transferase activity, transferring hexosyl groups
glucose-fructose oxidoreductase domain-containing protein 1-like
dual specificity tyrosine-phosphorylation-regulated kinase 1b
P:protein amino acid autophosphorylation
synaptobrevin homolog ykt6
serine--pyruvate mitochondrial precursor
transmembrane protein 192
ras-related protein rab-13-like
c-jun-amino-terminal kinase-interacting protein 4 isoform partial
histone-lysine n-methyltransferase setd3-like
peroxisome proliferator-activated receptor alpha
P:steroid hormone mediated signaling pathway
cytochrome p450 1a
C:endoplasmic reticulum membrane
thyrotrophic embryonic factor
P:cellular response to light stimulus
nuclear receptor subfamily 1 group d member 2
P:steroid hormone mediated signaling pathway
F:lipid transporter activity
RNA-seq recommendations for non-model taxa
The results of this study highlight the appropriateness of an RNA-seq approach for studies of adaptation (including adaptive plasticity) in non-model organisms. With the paucity of genomic resources available for most wildlife species, NGS technologies offer the best hope for unravelling the processes of evolutionary adaptation in a natural setting. Rainbowfish are evolutionarily very different from their nearest genome-enabled species, Oryzias latipes, yet in this study we were able to generate a substantial list of candidate genes involved in a response to increasing temperatures. Over the past few years, the proliferation of software resources and validated pipelines for RNA-seq means that virtually any organism can now be the focus of ecological genomic research and this is reflected in the rapid increase in publications reporting RNA-seq analyses in non-human taxa. The limiting factors that remain now are bioinformatic expertise and incomplete reference data. Over half of the dysregulated transfrags identified in our study were unable to be identified or were of unknown function. This continues to be a major challenge for studies of ecological and evolutionary genomics . Interpretation of genomic data lags well behind the current ability to generate that data. The limitation stems from the fact that annotation of genes of ecological interest still relies upon inferring homologies with genomic features established and developed in a few model species for non-ecological purposes. Better data integration is needed to facilitate the association of gene transcripts with specific natural conditions or phenotypic responses. Further work to characterise the function of these unknown genes via experimental studies of non-model organisms will enhance our understanding of the important biological pathways involved in responses to temperature stress and other environmental changes. We have shown that differing mapping and DE analysis approaches lead to very different outcomes in terms of the DE genes identified. While a combination of all available approaches is preferable to identify overlap in the candidate genes detected, we found that combining output from just Bowtie mapping and DESeq significance testing with BWA mapping and DESeq significance testing delivered just 21 more DE genes than combining all four approaches tested in our study (see Figure 1). This conservative approach is an efficient way to avoid large numbers of false positives being detected in RNA-seq studies.
Temperature increases predicted over the coming decades suggests species with limited dispersal abilities will need substantial adaptive potential to avoid extinction. That adaptive potential will likely come from a number of sources including adaptive phenotypic plasticity, standing genetic variation, and newly-derived mutations. Regardless of the source, adaptation will be most important in those processes related to heat tolerance. We have presented a first insight into which processes are likely to be important in the rainbowfish, M. duboulayi. This provides a foundation for future research into temperature-driven adaptive responses in freshwater species but also invites more detailed study of the phenome-genome interaction under conditions of temperature stress.
We identified a predictable suite of heat shock genes that responded sharply to increased temperatures in the treatment group. However, we also identified transfrags related to regulation of metabolic functions and developmental processes that showed mid-range levels of dysregulation and may be stronger candidates as genes for long-term adaptation to a warming environment. We present these candidate genes as targets for ongoing research into populations representing different thermal environments throughout the species range. We also expect that these candidates will be useful targets for studies of other freshwater species experiencing long-term thermal challenges.
The expression level changes we have presented may be an example of a plastic response. To check for an adaptive component it is necessary to repeat the temperature trial on other geographically distant populations and/or sister taxa. Parallel expression level changes in these populations would indicate plasticity whereas altered responses would be suggestive of adaptation at the genome level. Such “common garden” experiments allow the disentangling of pure plastic vs. genetic responses and are ideal approaches for future research. Other avenues to explore evolutionary adaptation to increased temperatures include investigating if DNA polymorphisms are present within and between populations at the gene regions we have identified in this study. Extensions of this research to include adaptive traits from other important environmental impacts will enable a much broader understanding of how freshwater species are likely to cope with human-induced habitat and climatic change.
Availability of supporting data
Raw sequencing data is available through the NCBI Sequence Read Archive under Project ID PRJNA205235 (http://trace.ncbi.nlm.nih.gov/Traces/sra/). All samples were sequenced as 101 bp paired-end reads on an Illumina HiSeq2000 sequencer.
This study was funded by the Discovery Program of the Australian Research Council (ARC grant DP110101207 to L. Beheregaray and L. Bernatchez). We thank Leo O’Reilly and Simon Westergaard for assistance with sampling and fish husbandry, respectively. Animal ethical approval was received from Flinders University (AWC E342).
- Shi Y, Li J, Jin Z: Advances in ecological genomics. Shengtai Xuebao/ Acta Ecologica Sinica. 2012, 32 (18): 5846-5858. 10.5846/stxb201108041143.View ArticleGoogle Scholar
- Ungerer MC, Johnson LC, Herman MA: Ecological genomics: understanding gene and genome function in the natural environment. Heredity. 2007, 100 (2): 178-183.View ArticlePubMedGoogle Scholar
- Ekblom R, Galindo J: Applications of next generation sequencing in molecular ecology of non-model organisms. Heredity. 2011, 107 (1): 1-15. 10.1038/hdy.2010.152.PubMed CentralView ArticlePubMedGoogle Scholar
- Hornett EA, Wheat CW: Quantitative RNA-Seq analysis in non-model species: assessing transcriptome assemblies as a scaffold and the utility of evolutionary divergent genomic reference species. BMC Genomics. 2012, 13 (1): 361-10.1186/1471-2164-13-361.PubMed CentralView ArticlePubMedGoogle Scholar
- Vijay N, Poelstra JW, Künstner A, Wolf JBW: Challenges and strategies in transcriptome assembly and differential gene expression quantification. 2012, Mol Ecol: A comprehensive in silico assessment of RNA-seq experimentsGoogle Scholar
- Pavey SA, Bernatchez L, Aubin-Horth N, Landry CR: What is needed for next-generation ecological and evolutionary genomics?. Ecology and Evolution. 2012, 27 (12): 673-678. 10.1016/j.tree.2012.07.014.View ArticlePubMedGoogle Scholar
- de Ong W, Voo LYC, Kumar VS: De Novo Assembly, Characterization and Functional Annotation of Pineapple Fruit Transcriptome through Massively Parallel Sequencing. PLoS One. 2012, 7 (10): e46937-10.1371/journal.pone.0046937.View ArticleGoogle Scholar
- Gagnaire PA, Normandeau E, Cote C, Hansen MM, Bernatchez L: The Genetic Consequences of Spatially Varying Selection in the Panmictic American Eel (Anguilla rostrata). Genetics. 2012, 190 (2): 725-U703. 10.1534/genetics.111.134825.PubMed CentralView ArticlePubMedGoogle Scholar
- Gao X, Han J, Lu Z, Li Y, He C: Characterization of the spotted seal Phoca largha transcriptome using Illumina paired-end sequencing and development of SSR markers. Comp Biochem Physiol Part D Genomics Proteomics. 2012, 7 (3): 277-284. 10.1016/j.cbd.2012.05.001.View ArticlePubMedGoogle Scholar
- Garcia TI, Shen Y, Crawford D, Oleksiak MF, Whitehead A, Walter RB: RNA-Seq reveals complex genetic response to deepwater horizon oil release in Fundulus grandis. BMC Genomics. 2012, 13 (1): 474-10.1186/1471-2164-13-474.PubMed CentralView ArticlePubMedGoogle Scholar
- Liu M, Qiao G, Jiang J, Yang H, Xie L, Xie J, Zhuo R: Transcriptome Sequencing and De Novo Analysis for Ma Bamboo (Dendrocalamus latiflorus Munro) Using the Illumina Platform. PLoS One. 2012, 7 (10): e46766-10.1371/journal.pone.0046766.PubMed CentralView ArticlePubMedGoogle Scholar
- Miller HC, Biggs PJ, Voelckel C, Nelson NJ: De novo sequence assembly and characterisation of a partial transcriptome for an evolutionarily distinct reptile, the tuatara (Sphenodon punctatus). BMC Genomics. 2012, 13 (1): 439-10.1186/1471-2164-13-439.PubMed CentralView ArticlePubMedGoogle Scholar
- St-Cyr J, Derome N, Bernatchez L: The transcriptomics of life-history trade-offs in whitefish species pairs (Coregonus sp.). Mol Ecol. 2008, 17 (7): 1850-1870. 10.1111/j.1365-294X.2008.03696.x.View ArticlePubMedGoogle Scholar
- Cooke G, Chao N, Beheregaray L: Natural selection in the water: freshwater invasion and adaptation by water colour in the Amazonian pufferfish. J Evol Biol. 2012, 25 (7): 1305-1320. 10.1111/j.1420-9101.2012.02514.x.View ArticlePubMedGoogle Scholar
- Lindquist S, Craig E: The heat-shock proteins. Annu Rev Genet. 1988, 22 (1): 631-677. 10.1146/annurev.ge.22.120188.003215.View ArticlePubMedGoogle Scholar
- Pipkin S, Rhodes C, Williams N: Influence of temperature on Drosophila alcohol dehydrogenase polymorphism. J Hered. 1973, 64 (4): 181-185.PubMedGoogle Scholar
- Crawford DL, Powers DA: Molecular basis of evolutionary adaptation at the lactate dehydrogenase-B locus in the fish Fundulus heteroclitus. Proc Natl Acad Sci. 1989, 86 (23): 9365-9369. 10.1073/pnas.86.23.9365.PubMed CentralView ArticlePubMedGoogle Scholar
- Kassahn KS, Caley MJ, Ward AC, Connolly AR, Stone G, Crozier RH: Heterologous microarray experiments used to identify the early gene response to heat stress in a coral reef fish. Mol Ecol. 2007, 16 (8): 1749-1763. 10.1111/j.1365-294X.2006.03178.x.View ArticlePubMedGoogle Scholar
- Lewis JM, Hori TS, Rise ML, Walsh PJ, Currie S: Transcriptome responses to heat stress in the nucleated red blood cells of the rainbow trout (Oncorhynchus mykiss). Physiol Genomics. 2010, 42 (3): 361-373. 10.1152/physiolgenomics.00067.2010.View ArticlePubMedGoogle Scholar
- Long Y, Li L, Li Q, He X, Cui Z: Transcriptomic Characterization of Temperature Stress Responses in Larval Zebrafish. PLoS One. 2012, 7 (5): e37209-10.1371/journal.pone.0037209.PubMed CentralView ArticlePubMedGoogle Scholar
- Buckley BA, Gracey AY, Somero GN: The cellular response to heat stress in the goby Gillichthys mirabilis: A cDNA microarray and protein-level analysis. J Exp Biol. 2006, 209 (14): 2660-2677. 10.1242/jeb.02292.View ArticlePubMedGoogle Scholar
- Logan CA, Somero GN: Transcriptional responses to thermal acclimation in the eurythermal fish Gillichthys mirabilis (Cooper 1864). Am J Physiol Regul Integr Comp Physiol. 2010, 299 (3): R843-R852. 10.1152/ajpregu.00306.2010.View ArticlePubMedGoogle Scholar
- Logan CA, Somero GN: Effects of thermal acclimation on transcriptional responses to acute heat stress in the eurythermal fish Gillichthys mirabilis (Cooper). Am J Physiol Regul Integr Comp Physiol. 2011, 300 (6): 1373-1383. 10.1152/ajpregu.00689.2010.View ArticleGoogle Scholar
- Quinn NL, McGowan CR, Cooper GA, Koop BF, Davidson WS: Identification of genes associated with heat tolerance in arctic charr exposed to acute thermal stress. Physiol Genomics. 2011, 43 (11): 685-696. 10.1152/physiolgenomics.00008.2011.View ArticlePubMedGoogle Scholar
- Buckley BA, Place SP, Hofmann GE: Regulation of heat shock genes in isolated hepatocytes from an Antarctic fish. Trematomus bernacchii. J Exp Biol. 2004, 207 (21): 3649-3656. 10.1242/jeb.01219.View ArticlePubMedGoogle Scholar
- Buckley BA, Hofmann GE: Thermal acclimation changes DNA-binding activity of heat shock factor 1 (HSF1) in the goby Gillichthys mirabilis: Implications for plasticity in the heat-shock response in natural populations. J Exp Biol. 2002, 205 (20): 3231-3240.PubMedGoogle Scholar
- West-Eberhard MJ: Phenotypic Plasticity and the Origins of Diversity. Annu Rev Ecol Evol Syst. 1989, 20: 249-278. 10.1146/annurev.es.20.110189.001341. ArticleType: research-article / Full publication date: 1989 / Copyright © 1989 Annual ReviewsView ArticleGoogle Scholar
- Cheviron ZA, Whitehead A, Brumfield RT: Transcriptomic variation and plasticity in rufous-collared sparrows (Zonotrichia capensis) along an altitudinal gradient. Mol Ecol. 2008, 17 (20): 4556-4569. 10.1111/j.1365-294X.2008.03942.x.View ArticlePubMedGoogle Scholar
- Allen GR, Midgley SH, Allen M: Field guide to the freshwater fishes of. 2002, Western Australian Museum: AustraliaGoogle Scholar
- Hattori A, Warburton K: Microhabitat use by the rainbowfish, Melanotaenia duboulayi, in a subtropical Australian stream. J Ethol. 2003, 21 (1): 15-22.Google Scholar
- McGuigan K, Zhu D, Allen GR, Moritz C: Phylogenetic relationships and historical biogeography of melanotaeniid fishes in Australia and New Guinea. Mar Freshwat Res. 2000, 51 (7): 713-723. 10.1071/MF99159.View ArticleGoogle Scholar
- Brown C: Do female rainbowfish (Melanotaenia spp.) prefer to shoal with familiar individuals under predation pressure?. J Ethol. 2002, 20 (2): 89-94. 10.1007/s10164-002-0059-6.View ArticleGoogle Scholar
- Holdway DA, Hefferman J, Smith A: Multigeneration assessment of nonylphenol and endosulfan using a model Australian freshwater fish. Melanotaenia fluviatilis. Environ Toxicol. 2008, 23 (2): 253-262. 10.1002/tox.20329.View ArticlePubMedGoogle Scholar
- Ferris R, Wilson RS: The physiological arms race: Exploring thermal acclimation among interacting species. J Therm Biol. 2012, 37 (3): 236-242. 10.1016/j.jtherbio.2012.01.006.View ArticleGoogle Scholar
- McCairns R, Bernatchez L: Landscape genetic analyses reveal cryptic population structure and putative selection gradients in a large‒scale estuarine environment. Mol Ecol. 2008, 17 (17): 3901-3916. 10.1111/j.1365-294X.2008.03884.x.View ArticlePubMedGoogle Scholar
- Young MJ: The efficacy of the aquatic anaesthetic AQUI-S® for anaesthesia of a small freshwater fish. Melanotaenia australis. J Fish Biol. 2009, 75 (7): 1888-1894. 10.1111/j.1095-8649.2009.02431.x.View ArticlePubMedGoogle Scholar
- Flanagan SW, Ryan AJ, Gisolfi CV, Moseley PL: Tissue-specific HSP70 response in animals undergoing heat stress. Am J Physiol Regul Integr Comp Physiol. 1995, 268 (1): R28-R32.Google Scholar
- Hall DM, Xu L, Drake VJ, Oberley LW, Oberley TD, Moseley PL, Kregel KC: Aging reduces adaptive capacity and stress protein expression in the liver after heat stress. J Appl Physiol. 2000, 89 (2): 749-759.PubMedGoogle Scholar
- Kew M, Bersohn I, Seftel H, Kent G: Liver damage in heatstroke. Am J Med. 1970, 49 (2): 192-202. 10.1016/S0002-9343(70)80075-4.View ArticlePubMedGoogle Scholar
- Rabergh CM, Airaksinen S, Soitamo A, Bjorklund HV, Johansson T, Nikinmaa M, Sistonen L: Tissue-specific expression of zebrafish (Danio rerio) heat shock factor 1 mRNAs in response to heat stress. J Exp Biol. 2000, 203 (12): 1817-1824.PubMedGoogle Scholar
- Giardine B, Riemer C, Hardison RC, Burhans R, Elnitski L, Shah P, Zhang Y, Blankenberg D, Albert I, Taylor J, et al: Galaxy: A platform for interactive large-scale genome analysis. Genome Res. 2005, 15 (10): 1451-1455. 10.1101/gr.4086505.PubMed CentralView ArticlePubMedGoogle Scholar
- Schulz MH, Zerbino DR, Vingron M, Birney E: 2012, Oases: Robust de novo RNA-seq assembly across the dynamic range of expression levels. BioinformaticsGoogle Scholar
- Surget-Groba Y, Montoya-Burgos JI: Optimization of de novo transcriptome assembly from next-generation sequencing data. Genome Res. 2010, 20 (10): 1432-1440. 10.1101/gr.103846.109.PubMed CentralView ArticlePubMedGoogle Scholar
- Langmead B, Trapnell C, Pop M, Salzberg S: Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009, 10 (3): R25-10.1186/gb-2009-10-3-r25.PubMed CentralView ArticlePubMedGoogle Scholar
- Li H, Durbin R: Fast and accurate short read alignment with Burrows–Wheeler transform. Bioinformatics. 2009, 25 (14): 1754-1760. 10.1093/bioinformatics/btp324.PubMed CentralView ArticlePubMedGoogle Scholar
- Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R, Subgroup GPDP: The Sequence Alignment/Map format and SAMtools. Bioinformatics. 2009, 25 (16): 2078-2079. 10.1093/bioinformatics/btp352.PubMed CentralView ArticlePubMedGoogle Scholar
- Anders S, Huber W: Differential expression analysis for sequence count data. Genome Biol. 2010, 11 (10): R106-10.1186/gb-2010-11-10-r106.PubMed CentralView ArticlePubMedGoogle Scholar
- Conesa A, Götz S, Garcia-Gomez JM, Terol J, Talon M, Robles M: Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005, 21: 3674-3676. 10.1093/bioinformatics/bti610.View ArticlePubMedGoogle Scholar
- Li H, Homer N: A survey of sequence alignment algorithms for next-generation sequencing. Brief Bioinform. 2010, 11 (5): 473-483. 10.1093/bib/bbq015.PubMed CentralView ArticlePubMedGoogle Scholar
- Robles JA, Qureshi SE, Stephen SJ, Wilson SR, Burden CJ, Taylor JM: Efficient experimental design and analysis strategies for the detection of differential expression using RNA-Sequencing. BMC Genomics. 2012, 13 (1): 484-10.1186/1471-2164-13-484.PubMed CentralView ArticlePubMedGoogle Scholar
- Kvam VM, Liu P, Si Y: A comparison of statistical methods for detecting differentially expressed genes from RNA-seq data. Am J Bot. 2012, 99 (2): 248-256. 10.3732/ajb.1100340.View ArticlePubMedGoogle Scholar
- Airaksinen S, Jokilehto T, Råbergh CM, Nikinmaa M: Heat- and cold-inducible regulation of HSP70 expression in zebrafish ZF4 cells. Comp Biochem Physiol B Biochem Mol Biol. 2003, 136 (2): 275-282. 10.1016/S1096-4959(03)00205-7.View ArticlePubMedGoogle Scholar
- Kassahn K, Crozier R, Ward A, Stone G, Caley MJ: From transcriptome to biological function: environmental stress in an ectothermic vertebrate, the coral reef fish Pomacentrus moluccensis. BMC Genomics. 2007, 8 (1): 358-10.1186/1471-2164-8-358.PubMed CentralView ArticlePubMedGoogle Scholar
- Keller JM, Escara-Wilke JF, Keller ET: Heat stress-induced heat shock protein 70 expression is dependent on ERK activation in zebrafish (Danio rerio) cells. Comparative Biochemistry and Physiology Part A: Molecular & Integrative Physiology. 2008, 150 (3): 307-314. 10.1016/j.cbpa.2008.03.021.View ArticleGoogle Scholar
- Krone PH, Lele Z, Sass JB: Heat shock genes and the heat shock response in zebrafish embryos. Biochem Cell Biol. 1997, 75 (5): 487-497. 10.1139/o97-083.View ArticlePubMedGoogle Scholar
- Murtha JM, Keller ET: Characterization of the heat shock response in mature zebrafish (Danio rerio). Exp Gerontol. 2003, 38 (6): 683-691. 10.1016/S0531-5565(03)00067-6.View ArticlePubMedGoogle Scholar
- Quinn NL, McGowan CR, Cooper GA, Koop BF, Davidson WS: Ribosomal genes and heat shock proteins as putative markers for chronic, sublethal heat stress in arctic charr: Applications for aquaculture and wild fish. Physiol Genomics. 2011, 43 (18): 1056-1064. 10.1152/physiolgenomics.00090.2011.View ArticlePubMedGoogle Scholar
- Fitzpatrick BM: Underappreciated consequences of phenotypic plasticity for ecological speciation. International Journal of Ecology. 2012, 2012: Art. ID 25601-View ArticleGoogle Scholar
- Kültz D: Molecular and evolutionary basis of the cellular stress response. Annu Rev Physiol. 2005, 67: 225-257. 10.1146/annurev.physiol.67.040403.103635.View ArticlePubMedGoogle Scholar
- Thibert-Plante X, Hendry A: The consequences of phenotypic plasticity for ecological speciation. J Evol Biol. 2011, 24 (2): 326-342. 10.1111/j.1420-9101.2010.02169.x.View ArticlePubMedGoogle Scholar
This article is published under license to 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.