Gene expression in the brain and kidney of rainbow trout in response to handling stress
BMC Genomics volume 6, Article number: 3 (2005)
Microarray technologies are rapidly becoming available for new species including teleost fishes. We constructed a rainbow trout cDNA microarray targeted at the identification of genes which are differentially expressed in response to environmental stressors. This platform included clones from normalized and subtracted libraries and genes selected through functional annotation. Present study focused on time-course comparisons of stress responses in the brain and kidney and the identification of a set of genes which are diagnostic for stress response.
Fish were stressed with handling and samples were collected 1, 3 and 5 days after the first exposure. Gene expression profiles were analysed in terms of Gene Ontology categories. Stress affected different functional groups of genes in the tissues studied. Mitochondria, extracellular matrix and endopeptidases (especially collagenases) were the major targets in kidney. Stress response in brain was characterized with dramatic temporal alterations. Metal ion binding proteins, glycolytic enzymes and motor proteins were induced transiently, whereas expression of genes involved in stress and immune response, cell proliferation and growth, signal transduction and apoptosis, protein biosynthesis and folding changed in a reciprocal fashion. Despite dramatic difference between tissues and time-points, we were able to identify a group of 48 genes that showed strong correlation of expression profiles (Pearson r > |0.65|) in 35 microarray experiments being regulated by stress. We evaluated performance of the clone sets used for preparation of microarray. Overall, the number of differentially expressed genes was markedly higher in EST than in genes selected through Gene Ontology annotations, however 63% of stress-responsive genes were from this group.
1. Stress responses in fish brain and kidney are different in function and time-course. 2. Identification of stress-regulated genes provides the possibility for measuring stress responses in various conditions and further search for the functionally related genes.
Until recently multiple gene expression profiling was applied almost exclusively to human and a few model organisms. At present cDNA microarrays are being constructed for new species including teleost fishes [1–6]. Since EST sequencing projects are carried out with a large number of species, continuous development of new platforms can be expected in the future. We designed a salmonid fish cDNA microarray primarily to characterize responses to stress, toxicity and pathogens. This paper focuses on time-course comparisons of stress responses in rainbow trout and the usage of functional annotation to conduct analyses of gene expression data.
Functional annotation of genes, especially Gene Ontology  is increasingly being used for analyses and interpretation of microarray results [8–13]. We applied Gene Ontology in several modes to facilitate implementation of our research tasks. Furthermore, experimental results generated guidelines for the development of specialized microarrays. Well designed platforms are expected to ensure identification of differentially expressed genes while containing representative coverage from important functional groups. Custom made microarrays include clones from cDNA libraries and/or selected genes, which have advantages and drawbacks. Indiscriminant spotting of EST may result in under representation of many functional classes. On the other hand selection of genes fully relies on annotations and hypotheses, which can be misleading and limit possibilities for nontrivial findings. We used clones from normalized and subtracted cDNA libraries as well as genes selected by the functional categories of Gene Ontology for inclusion onto a microarray targeted at characterizing transcriptome responses to environmental stressors. Designing a new platform requires balancing a large number of genes versus multiple replications of spots, which enhances statistical analyses of data. The rainbow trout microarray was prepared by spotting of relatively small number of genes (1300) in 6 replicates. We show that multiple replications combined with the dye-swap design of hybridization [14, 15] allows for accurate detection of relatively small alterations in expression levels, which is important for the functional interpretation of results.
Stress is closely associated with many diverse issues in fish biology and environmental research (reviewed in ). Stress is generally defined as the reaction to external forces and abnormal conditions that tend to disturb an organism's homeostasis. To illustrate the major trends in the studies of stress in fish, we performed a computer-assisted analysis of Medline abstracts covering this area (Table 1). Salmonids have been studied more extensively than any other fish species. Research has focused on various biotic and abiotic factors including toxicity, environmental parameters (oxygen, temperature, salinity, acidosis), diseases, social interactions (crowding, aggressiveness) and farming manipulations. Analysis of Medline abstracts indicated physiological processes, cellular structure and selected proteins that have been the major foci of previous fish stress studies. This provided an outline for interpretation of our results. We analyzed the effects of stress on the transcriptome in the brain and kidney, which are considered important target tissues along with muscle, blood cells, liver and epithelia. We report a profound difference of stress response in these tissues and the identification of a diagnostic set of genes.
1 Design of cDNA microarray
The rainbow trout cDNA microarray was composed of EST and selected genes. The cDNA libraries were prepared from tissues of stressed fish using suppression subtractive hybridization, SSH  and a modification of the cap-finder method  supplemented with enzymatic normalization . We sequenced 2000 clones and redundancy of the subtracted libraries was markedly greater than that of the normalized (306% and 134% respectively). In addition to EST we selected rainbow trout transcripts from the normalized multi-tissue cDNA library  based on their assignment to functional categories of Gene Ontology (stress and defense response, regulation of cell cycle, signal transduction, chaperone activity and apoptosis). The selected genes substantially improved the coverage of many functional classes (Table 2), though the number of differentially expressed genes in this group was markedly inferior to EST (Figure 1). Subtraction cloning enriched genes that showed strong alteration of expression at response to stress (p < 0.01 or lower, Figure 1A), however the SSH clone set did not provide any advantage when microarray was used for the related research tasks (Figure 1B).
2 Stress response in the brain and kidney of rainbow trout
2.1 Differentially expressed genes
Fish were stressed with netting and samples were collected 1, 3 and 5 days after the first exposure. We used plasma cortisol as a stress marker . The hormone levels increased 7.6-fold after 1 day and did not change significantly to the end of experiment (Figure 2).
Microarray results were submitted to GEO (GSM22355). Two genes were up-regulated in both tissues at all time-points (Figure 3). One is a putative homolog to the mammalian N-myc regulated genes, which are induced with steroid hormones in the brain  and kidney . Mitochondrial ADP, ATP carrier can be implicated to both normal functions and cell death . Metallothionein-IL, a classical stress marker was induced to the end of experiment and a similar profile was seen in midkine precursor (growth factor), histone H1.0 and B-cell translocation protein 1. In kidney we observed consistent up-regulation of genes related to energy metabolism, such as mitochondrial proteins (cytochromes b and c, cytochrome oxidases), enzymes (glyceraldehyde 3-phosphate dehydrogenase, fructose-bisphosphate aldolase, serine-pyruvate aminotransferase) and similar profiles were seen in two heat shock proteins and two signal transducers (cytohesin binding protein and GRB2-adaptor). The repressed genes were related to actin binding (coronin and profilin) and immune response (meprin, immunoglobulin epsilon receptor, thymosin and lysozyme).
Rapid alteration of gene expression was a remarkable feature of stress response in the brain. Only one gene, aquaporin, was up-regulated for the duration of the experiment. Water channel aquaporin plays a key role in water homeostasis being implicated in various physiological processes and pathological conditions . A panel of genes which showed markedly increased expression after 1 day was also suppressed after 5 days. Surprisingly, this group included mainly genes that are predominantly expressed in skeletal or cardiac muscle (myosin light chain 1 and 2, skeletal and cardiac isoforms, myosin heavy chain, troponin I, T and C) or are involved in regulation of muscle contraction (parvalbumin alpha and sarcoplasmic reticulum calcium ATPase). An opposite tendency was shown by a large group of genes however the magnitude of expression changes was smaller. We analysed 5 differentially expressed genes with qPCR and the results were in close concordance with the microarray data (not shown).
2.2 Functional classes
The search for enriched Gene Ontology functional categories in the lists of differentially expressed genes found almost no overlap between the tissues (Table 3). In the brain stress affected binding and transport of metal ions, especially calcium and manganese, chaperones and heat shock proteins, cytoskeleton and microtubules and a number of signaling pathways; whereas, mitochondrion, extracellular structures and peptidases appeared the primary targets in the kidney.
Comparison of the differentially expressed genes by the Gene Ontology categories suggested coordinated regulation of various cellular functions in the brain. Early stress response was marked with transient induction of the cytoskeleton proteins and similar profiles were observed in the metal binding proteins and enzymes of carbohydrate metabolism (Figure 4). An opposite expression pattern was shown by a large group of genes involved in stress and immune response, regulation of growth and cell cycle, apoptosis, signal transduction and cell to cell signaling. This was in parallel with enhancement of transcription and translation, ubiquitin-dependent protein catabolism and protein folding. In the kidney the temporal alterations were much weaker. Expression of metal binding proteins increased slowly in parallel with peptidases. Strong induction of collagenases coincided with decrease of collagen expression. At the same time a number of metabolic functions were suppressed (oxidative phosphorylation and oxidoreductase activity, amine metabolism and RNA binding).
3 Stress-responsive genes
Microarray design included genes from functional categories which were expected to be affected by stress (Table 2). Overall observations of differences in gene expression from this group in response to handling stress were minimal; however, this could be accounted for by its heterogeneity. Therefore we searched for the subgroups of genes with correlated expression profiles within the functional classes using results of 35 microarray experiments conducted by our laboratory. Both factorial and cluster analyses revealed 9 defense response genes that showed tightly coordinated expression being induced with stress. We continued search using the consensus profile of this subgroup and found 47 positively and 1 negatively correlated genes (Pearson r > |0.65|). Of these 29 were identified by the protein products (Figure 5A), 19 being from the set of selected clones. Expression of the stress-responsive genes changed significantly in several experiments including this study (Figure 5B). They were up-regulated in kidney with stress and injection of cortisol, combination of these treatments showed an additive effect (Figure 5C). These genes also responded to the model water contaminants, being induced with low and medium and repressed with high doses (Figure 5D).
I Stress response in rainbow trout
Our study aimed at comparison of time-course of stress response in rainbow trout brain and kidney and finding of a diagnostic set of genes. These tasks were implemented with an aid of Gene Ontology annotation, which was used in several modes. The most straightforward and commonly used approach is counting of Gene Ontology classes in the lists of differentially expressed genes. Statistical inference of enrichment and depletion is made with Z-score of hypergeometric distribution, exact Fisher's test or its modifications. Such analyses helped us to interpret differences of stress responses in the brain and kidney (Table 3). In the brain handling stress mainly affected expression of transcripts for structural proteins (especially cytoskeleton), signal transduction, and binding of metal ions, whereas mitochondria, extracellular structures and peptidases appeared the key targets in the kidney. Computer-assisted analysis of Medline abstracts suggested that most of these themes have not been addressed in the studies of fish stress (Table 1).
Searches of the enriched Gene Ontology categories associated with differentially expressed transcripts is useful for rapid screens of microarray data; however, it presumes coordinated expression of functionally related genes. This assumption is not valid for many classes, especially large and heterogenous groups, such as stress and defense response. Because the gene composition of microarray is used as a reference, uneven presentation of functional categories can distort the results. Finally, this analysis does not take into account direction and magnitude of differential expression. To overcome these problems, enrichment of Gene Ontology classes is analysed in groups of genes with similar expression profiles revealed with cluster or factorial analyses. In this study we preferred straight comparison of Gene Ontology classes by the mean log expression ratios which helped in interpretation of the time-course of stress response.
In the kidney temporal alterations were relatively weak though significant. Expression of peptidases (especially collagenases) increased steadily, which implied possible degradation of tissue with prolonged stress. We could expect abrupt fluctuations in the rainbow trout brain, since transient induction and up-regulation of gene expression was observed in response to cold in the brain of channel catfish . In our study most differentially expressed genes fell into two groups with distinct temporal profiles which showed remarkable coherence of the functional classes. Early phase was associated with dramatic up-regulation of structural and metal binding proteins, which were repressed in later phases. Expression of genes involved in stress and defense response, apoptosis and signal transduction, cell cycle and growth changed in a reciprocal fashion. Activation of metal binding proteins could be accounted for the role of ions (particularly calcium) in multiple pathways of gene expression regulation in the brain . Motor proteins of cytoskeleton play key roles in the transport of vesicles and the establishment and rearrangement of neuronal networks [27–29] which also could be implicated to the stress response in fish. However, in mammals these functions are associated with non-muscle isoforms and therefore differential expression of the sarcomeric proteins was unexpected. Additional experiments confirmed induction of these proteins at early phase of stress response. Previously we observed high activity of skeletal α-actin and myosin light chain 2 promoters in the neural tissues of rainbow trout embryos . Sequencing of salmonid fish cDNA libraries provided evidence for the brain expression of sarcomeric proteins, but their role remains fully unknown. At present there is sparse evidence for differential expression of structural muscle proteins in the mammalian brain. For example regulation of troponin I with dextromethorphan (antagonist of excitatory amino acid receptors) was reported in the rat hyppocampus and cortex .
Grouping of individual differentially expressed genes by the functional classes reduced noise and enhanced cluster and factorial analyses. This helped to identify stress-responsive genes that showed correlated expression in 35 microarray experiments (22 experiments are shown in Figure 5A). Association with stress is well established for most of these proteins and some are used as stress markers. The list of enriched Gene Ontology categories (stress, defense and humoral immune response, signal transduction and response to oxidative stress, p < 0.05) suggested biological relevance of this group. Computer analysis of Medline abstracts (Table 1) showed that immunity and metabolism of reactive oxidative species are prioritized in studies of fish stress and these functional categories were enriched in the list of stress-responsive genes. Thus Gene Ontology provided a useful starting point for search of functionally related genes and results of these analyses can be used further for the revision of annotations.
2 Construction of microarrays
Results of our experiments helped to evaluate the strategy used in construction of the rainbow trout microarray. Researchers developing microarrays for new species are commonly choosing between specially selected genes and clones from normalized and subtracted cDNA libraries. We used SSH, which is at present probably the most popular method of subtraction. Though proven efficient in many studies, this method has a number of drawbacks. Subtraction requires re-association of tester and formation of double-stranded DNA, hence many rare transcripts are not cloned and variations in concentrations of cDNA and hybridization conditions may have strong impacts on library composition. High redundancy is a common feature of the SSH libraries. Apart from these problems, rapid alterations of gene expression observed in this study and many other microarray experiments make the advantages of subtractive cloning ambiguous. Subtraction achieves enrichment of the transcripts, which are over or under represented in the test sample. In many cases one sample will not provide coverage of differentially expressed genes for the whole series, whereas pooling of samples may reduce fluctuations. Furthermore, we observed relatively high ratio of differentially expressed genes among the clones from the unsubtracted cDNA libraries, which are easier for construction and much less redundant. The advantage of subtractive cloning becomes negligible when microarrays are used for different, though related research tasks (Fig. 1B).
At present selection of genes for microarrays is facilitated with advances of functional annotation. This helped us to improve presentation of many functional categories (Table 1) and enhanced interpretation of results. Most of the selected genes did not show differential expression in our studies, however 63% of stress-responsive genes were from this group. In our view, this finding is a strong argument for utilizing Gene Ontology in the development of specialized platforms.
Given the limited number of spots on slides, microarray design requires a careful balance between the number of genes and replication of spots. Apparent advantages of genome-wide platforms are compromised with the problems associated with identifying significantly differentially expressed genes. We preferred combination of multiple spotting and dye-swap normalization, which ensured robust normalization and accurate detection of differential expression at low ratios. Coordinated expression of functionally related genes suggested biological relevance of relatively small alterations in the transcription levels. Selection of differentially expressed genes by the cutoff values would result in loss of valuable information in our experiments. For instance, most of the stress-responsive genes showed small or moderate expression changes, the identification of this group would not be likely without multiple replications.
1. Combination of EST and selected genes appears a reasonable way for construction of cDNA microarrays. Multiple replications of spots and dye swap design of hybridization ensure robust normalization and high power of statistical analyses. Finding of differential expression at small ratios is essential for the functional interpretation of microarray data.
2. Stress response in fish brain and kidney is different both by the target functions and time-course. In brain slow progression of adaptive response was preceded with dramatic transient induction of motor and metal ion binding proteins. Prolonged stress was likely to result in slow degradation of extracellular matrix in kidney.
3. Finding of stress-responsive genes provides possibility for measurement of stress in various conditions and search for the functionally related genes.
1. Computer-assisted analysis of Medline abstracts
Search of Medline was made with queries: "fishes AND stress" (1060 abstracts) and "fishes NOT stress" (10069 abstracts). Abstracts were split into separate words and a list of non-redundant terms was composed. The numbers of abstracts including each term were estimated. The terms were ranked by the Z-scores of hypergeometric distribution and enrichment was analysed with exact Fisher's test (p < 0.05).
2. Experiments with fish, exposure and sampling
One year old rainbow trout were stressed with netting for 2 min, this treatment was repeated once a day for a duration of 5 days. Fish were killed with over-dose of anaesthetic (MS-222) and blood was taken from the caudal vein. The kidneys and brains were snap-frozen in liquid nitrogen. Plasma cortisol was determined with RIA using Orion Spectra Cortisol kit.
3. Preparation of microarrays
RNA was extracted with Trizol reagent (Invitrogen) and mRNA was purified with Dynabeads kit (Dynal). SSH cloning was performed as described . For preparation of normalized libraries, synthesis of cDNA with PowerScript reverse transcriptase (Clontech) was primed with oligonucleotides including EcoRI and NotI sites: 5'-ACGAGGCGAATTC ACAGAGAGGAG(T)VN-3', 5'-GAGAGAGAGTGGTGCGGCCGC GGTGTATGGGG-3'). Double-stranded cDNA was generated using Advantage DNA polymerase mix (Clontech) and PCR primers: 5'-ACGAGGCGAATTC ACAGAGAGGAG-3' and 5'-GAGAGTGGTGCGGCCGC GGTGTA-3'. The PCR products were purified with QIAquick kit (Qiagen), precipitated with ethanol and dissolved to 1 μg/μl in hybridization buffer (1 M NaCl, 50 mM HEPES (pH 8.3), 1 mM EDTA). DNA was denaturated for 5 min at 94°C. Following re-association at 72°C for 16 hours, DNA was ethanol precipitated and digested with 150 U of exonuclease III (MBI Fermentas) for 15 min at 37°C. This treatment eliminates re-associated double-stranded DNA . Single-stranded DNA was PCR amplified, size separated with agarose gel electrophoresis and cloned into pGEM®-11Zf (+) (Promega). Normalized and subtracted cDNA libraries were prepared from the stressed fish (whole fry, brain, kidney and spleen of 1-year old fish). A number of clones were from the rainbow trout and Baltic salmon cDNA libraries constructed in University of Turku. The sequences were analysed with stand-alone blastn and blastx  Microarray incldued 315 genes selected by the Gene Ontology functional categories. Of these, 282 clones were from the normalized multi-tissue library  and the rest were produced with RT PCR. The cDNA inserts were amplified with PCR using universal primers and purified with Millipore Montage PCR96 Cleanup Kit. DNA was spotted onto poly-(L) lysine-coated slides and each clone was printed in 6 replicates.
4. Microarray analyses
Total RNA was extracted with Trizol reagent (Invitrogen) and 4 individuals were pooled in each sample. Stressed fish was compared with time-matched control. Labeling with Cy3- and Cy5-dCTP (Amersham Pharmacia) was made using SuperScript III (Invitrogen) and oligo(dT) primer; cDNA was purified with Microcon YM30 (Millipore). We used a dye swap experimental design [14, 15] and each sample was hybridized to two microarrays. For the first slide, test and control cDNA were labeled with Cy5 and Cy3 respectively, and for the second array dye assignments were reversed. The slides were pretreated with 1% BSA, fraction V, 5 x SSC, 0.1% SDS (30 min at 50°C) and washed with 2 x SSC (3 min) and 0.2 x SSC (3 min) and hybridized overnight in cocktail containing 1.3 x Denhardt's, 3 x SSC 0.3% SDS, 0.67 μg/μl polyadenylate and 1.4 μg/μl yeast tRNA. All chemicals were from Sigma-Aldrich. Scanning was performed with ScanArray 5000 and images were processed with QuantArray (GSI Luminomics). The measurements in spots were filtered by criteria I/B ≥ 3 and (I-B)/(S I +S B ) ≥ 0.6, where I and B are the mean signal and background intensities and S I , S B are the standard deviations. After subtraction of mean background, lowess normalization  was performed. Differential expression was analysed with Student's t-test (p < 0.01) and the genes were ranked by the log(p-level).
5. Quantitative RT PCR
Primers (Table 4) were designed to amplify 194–305 b fragments. RNA was processed with Rnase-free Dnase (Promega). Synthesis of cDNA with Superscript III reverse transcriptase (Invitrogen) was primed with oligo(dT). Analyses were carried out using Dynamo SYBR Green kit (Finnzymes) and ABI Prism 7700 (Amersham-Pharmacia).
AK designed microarray, carried out experiments with fish and data analyses and drafted the manuscript. HK conducted the microarray analyses. PP developed software for annotation of genes and analyses of Medline abstracts. CR constructed the multi-tissue cDNA library and provided the selected genes. SA developed software for management of microarray data and performed the statistical analyses. HM coordinated research. All authors read and approved the final manuscript.
Gracey AY, Troll JV, Somero GN: Hypoxia-induced gene expression profiling in the euryoxic fish Gillichthys mirabilis. Proc Natl Acad Sci U S A. 2001, 98: 1993-1998. 10.1073/pnas.98.4.1993.
Rise ML, von Schalburg KR, Brown GD, Mawer MA, Devlin RH, Kuipers N, Busby M, Beetz-Sargent M, Alberto R, Gibbs AR, Hunt P, Shukin R, Zeznik JA, Nelson C, Jones SR, Smailus DE, Jones SJ, Schein JE, Marra MA, Butterfield YS, Stott JM, Ng SH, Davidson WS, Koop BF: Development and application of a salmonid EST database and cDNA microarray: data mining and interspecific hybridization characteristics. Genome Res. 2004, 14: 478-490. 10.1101/gr.1687304.
Williams TD, Gensberg K, Minchin SD, Chipman JK: A DNA expression array to detect toxic stress response in European flounder (Platichthys flesus). Aquat Toxicol. 2003, 65: 141-157. 10.1016/S0166-445X(03)00119-X.
Ju Z, Dunham RA, Liu Z: Differential gene expression in the brain of channel catfish (Ictalurus punctatus) in response to cold acclimation. Mol Genet Genomics. 2002, 268: 87-95. 10.1007/s00438-002-0727-9.
Oleksiak MF, Churchill GA, Crawford DL: Variation in gene expression within and among natural populations. Nat Genet. 2002, 32: 261-266. 10.1038/ng983.
Ton C, Stamatiou D, Dzau VJ, Liew CC: Construction of a zebrafish cDNA microarray: gene expression profiling of the zebrafish during development. Biochem Biophys Res Commun. 2002, 296: 1134-1142. 10.1016/S0006-291X(02)02010-7.
Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, Harris MA, Hill DP, Issel-Tarver L, Kasarskis A, Lewis S, Matese JC, Richardson JE, Ringwald M, Rubin GM, Sherlock G: Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet. 2000, 25: 25-29. 10.1038/75556.
Doniger SW, Salomonis N, Dahlquist KD, Vranizan K, Lawlor SC, Conklin BR: MAPPFinder: using Gene Ontology and GenMAPP to create a global gene-expression profile from microarray data. Genome Biol. 2003, 4: R7-10.1186/gb-2003-4-1-r7.
Draghici S, Khatri P, Martins RP, Ostermeier GC, Krawetz SA: Global functional profiling of gene expression. Genomics. 2003, 81: 98-104. 10.1016/S0888-7543(02)00021-6.
Herrero J, Al-Shahrour F, Diaz-Uriarte R, Mateos A, Vaquerizas JM, Santoyo J, Dopazo J: GEPAS: A web-based resource for microarray gene expression data analysis. Nucleic Acids Res. 2003, 31: 3461-3467. 10.1093/nar/gkg591.
Zeeberg BR, Feng W, Wang G, Wang MD, Fojo AT, Sunshine M, Narasimhan S, Kane DW, Reinhold WC, Lababidi S, Bussey KJ, Riss J, Barrett JC, Weinstein JN: GoMiner: a resource for biological interpretation of genomic and proteomic data. Genome Biol. 2003, 4: R28-10.1186/gb-2003-4-4-r28.
Robinson PN, Wollstein A, Bohme U, Beattie B: Ontologizing gene-expression microarray data: characterizing clusters with Gene Ontology. Bioinformatics. 2004, 20: 979-981. 10.1093/bioinformatics/bth040.
McCarroll SA, Murphy CT, Zou S, Pletcher SD, Chin CS, Jan YN, Kenyon C, Bargmann CI, Li H: Comparing genomic expression patterns across species identifies shared transcriptional profile in aging. Nat Genet. 2004, 36: 197-204. 10.1038/ng1291.
Kerr MK, Churchill GA: Statistical design and the analysis of gene expression microarray data. Genet Res. 2001, 77: 123-128. 10.1017/S0016672301005055.
Yang YH, Dudoit S, Luu P, Speed T: Normalization for cDNA Microarray Data. In Microarray Optical Technologies and Informatics. Edited by: Bittner ML, Chen Y, Dorsel AN, Dougherty ER. 2001, San Jose: SPIE, Society for Optical Engineering
Wendelaar Bonga SE: The stress response in fish. Physiol Rev. 1997, 77: 591-625.
Diatchenko L, Chenchik A, Siebert P: Suppression subtractive hybridization: A method for generating subtracted cDNA libraries starting from poly (A+) or total RNA. In RT-PCR Methods for Gene Cloning and Analysis. Edited by: Siebert P, Larrick J. 1998, BioTechniques Books, 213-239.
Chenchik A, Zhu Y, Diatchenko L, Li R, Hill J, Siebert P: Generation and use of high quality cDNA from small amounts of total RNA by SMART PCR. In RT-PCR Methods for Gene Cloning and Analysis. Edited by: Siebert P, Larrick J. 1998, BioTechniques Books, 305-319.
Zeng J, Gorski RA, Hamer D: Differential cDNA cloning by enzymatic degrading subtraction (EDS). Nucleic Acids Res. 1994, 22: 4381-4385.
Rexroad CE, Keele JW, Gahr SA, Palti Y, Lee Y, Karamycheva S, Quackenbush J: Sequence analysis of a rainbow trout cDNA library and creation of a Gene Index. Cytogenetics and Genome Research. 2003, 102: 347-354. 10.1159/000075773.
Barton BA, Iwama GK: Physiological changes in fish from stress in aquaculture with emphasis on the response and effects of corticosteroids. Annual Review of Fish Diseases. 1991, 1: 3-26. 10.1016/0959-8030(91)90019-G.
Nichols NR: Ndrg2, a novel gene regulated by adrenal steroids and antidepressants, is highly expressed in astrocytes. Ann N Y Acad Sci. 2003, 1007: 349-356. 10.1196/annals.1286.034.
Boulkroun S, Fay M, Zennaro MC, Escoubet B, Jaisser F, Blot-Chabaud M, Farman N, Courtois-Coutry N: Characterization of rat NDRG2 (N-Myc downstream regulated gene 2), a novel early mineralocorticoid-specific induced gene. J Biol Chem. 2002, 277: 31506-31515. 10.1074/jbc.M200272200.
Halestrap AP, Brennerb C: The adenine nucleotide translocase: a central component of the mitochondrial permeability transition pore and key player in cell death. Curr Med Chem. 2003, 10: 1507-1525.
Amiry-Moghaddam M, Ottersen OP: The molecular basis of water transport in the brain. Nat Rev Neurosci. 2003, 4: 991-1001.
West AE, Chen WG, Dalva MB, Dolmetsch RE, Kornhauser JM, Shaywitz AJ, Takasu MA, Tao X, Greenberg ME: Calcium regulation of neuronal gene expression. Proc Natl Acad Sci U S A. 2001, 98: 11024-11031. 10.1073/pnas.191352298.
Brown J, Bridgman PC: Role of myosin II in axon outgrowth. J Histochem Cytochem. 2003, 51: 421-428.
Bridgman PC: Myosin-dependent transport in neurons. J Neurobiol. 2004, 58: 164-174. 10.1002/neu.10320.
Brown ME, Bridgman PC: Myosin function in nervous and sensory systems. J Neurobiol. 2004, 58: 118-130. 10.1002/neu.10285.
Krasnov A, Teerijoki H, Gorodilov Y, Mölsa H: Cloning of rainbow trout (Oncorhynchus mykiss) alpha-actin, myosin regulatory light chain genes and the 5'-flanking region of alpha-tropomyosin. Functional assessment of promoters. J Exp Biol. 2003, 206: 601-608. 10.1242/jeb.00116.
Lee KH, Ahn JI, Yu DH, Koh HC, Kim SH, Yang BH, Lee YS, Lee YS: Dextromethorphan alters gene expression in rat brain hippocampus and cortex. Int J Mol Med. 2003, 11: 559-568.
Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25: 3389-3402. 10.1093/nar/25.17.3389.
Cleveland WS, Grosse E, Shyu WM: Local regression models. In Statistical Models in S. Edited by: Chambers JM, Hastie TJ. 1992, Pacific Grove: Wadsworth & Brooks/Cole, 309-376.
Koskinen H, Pehkonen P, Vehniäinen E, Krasnov A, Rexroad C, Afanasyev S, Mölsa H, Oikari A: Response of rainbow trout transcriptome to model chemical contaminants. Biochem Biophys Res Communs. 2004, 320: 745-753. 10.1016/j.bbrc.2004.06.024.
This study was supported by the National Agency of Technology, Finland. We wish to thank Rolf Sara (CBT, University of Turku) for preparation of cDNA microarrays and Seppo Kukkonen for analyses of plasma cortisol.
Heikki Koskinen, Petri Pehkonen contributed equally to this work.
About this article
Cite this article
Krasnov, A., Koskinen, H., Pehkonen, P. et al. Gene expression in the brain and kidney of rainbow trout in response to handling stress. BMC Genomics 6, 3 (2005). https://doi.org/10.1186/1471-2164-6-3