National Center for Cool and Cold Water Aquaculture

The whole-genome duplication 80 million years ago of the common ancestor of salmonids (salmonid-specific fourth vertebrate whole-genome duplication, Ss4R) provides unique opportunities to learn about the evolutionary fate of a duplicated vertebrate genome in 70 extant lineages. Here we present a high-quality genome assembly for Atlantic salmon (Salmo salar), and show that large genomic reorganizations, coinciding with bursts of transposon-mediated repeat expansions, were crucial for the post-Ss4R rediploidization process. Comparisons of duplicate gene expression patterns across a wide range of tissues with orthologous genes from a pre-Ss4R outgroup unexpectedly demonstrate far more instances of neofunctionalization than subfunctionalization. Surprisingly, we find that genes that were retained as duplicates after the teleost-specific whole-genome duplication 320 million years ago were not more likely to be retained after the Ss4R, and that the duplicate retention was not influenced to a great extent by the nature of the predicted protein interactions of the gene products. Finally, we demonstrate that the Atlantic salmon assembly can serve as a reference sequence for the study of other salmonids for a range of purposes.

Infectious diseases of fishes occur when susceptible fishes are exposed to virulent pathogens under certain environmental stress conditions. Very little research has been carried out to show the effect of pollution on outbreaks of infectious diseases of fishes. Therefore, examples taken from the literature were selected and reviewed to show the coincidence of infectious diseases with stress caused by temperature, eutrophication, sewage, metabolic products of fishes, industrial pollution, and pesticides.

Brown marmorated stink bug, <it>Halyomorpha halys</it> Stål, is an invasive, herbivorous insect species that was accidentally introduced to the United States from Asia. First discovered in Allentown, PA, in 1996, <it>H. halys</it> has now been reported from at least 40 states in the United States. Additional invasions have been detected in Canada, Switzerland, France, Germany, Italy, and Lichtenstein, suggesting this invasive species could emerge as a cosmopolitan pest species. In its native range, <it>H. halys</it> is classified as an outbreak pest; however, in North America, <it>H. halys</it> has become a major agricultural pest across a wide range of commodities. <it>H. halys</it> is a generalist herbivore, capable of consuming >100 different species of host plants, often resulting in substantial economic damage; its feeding damage resulted in US$37 million of losses in apple in 2010, but this stink bug species also attacks other fruit, vegetable, field crop, and ornamental plant species. <it>H. halys</it> has disrupted integrated pest management programs for multiple cropping systems. Pesticide applications, including broad-spectrum insecticides, have increased in response to <it>H. halys</it> infestations, potentially negatively influencing populations of beneficial arthropods and increasing secondary pest outbreaks. <it>H. halys</it> is also challenging because it affects homeowners as a nuisance pest; the bug tends to overwinter in homes and outbuildings. Although more research is required to better understand the ecology and biology of <it>H. halys</it>, we present its life history, host plant damage, and the management options available for this invasive pest species.

We describe the organization of a nascent international effort, the Functional Annotation of Animal Genomes (FAANG) project, whose aim is to produce comprehensive maps of functional elements in the genomes of domesticated animal species.

Antimicrobial resistance in Yersinia pestis is rare, yet constitutes a significant international public health and biodefense threat. In 1995, the first multidrug resistant (MDR) isolate of Y. pestis (strain IP275) was identified, and was shown to contain a self-transmissible plasmid (pIP1202) that conferred resistance to many of the antimicrobials recommended for plague treatment and prophylaxis. Comparative analysis of the DNA sequence of Y. pestis plasmid pIP1202 revealed a near identical IncA/C plasmid backbone that is shared by MDR plasmids isolated from Salmonella enterica serotype Newport SL254 and the fish pathogen Yersinia ruckeri YR71. The high degree of sequence identity and gene synteny between the plasmid backbones suggests recent acquisition of these plasmids from a common ancestor. In addition, the Y. pestis pIP1202-like plasmid backbone was detected in numerous MDR enterobacterial pathogens isolated from retail meat samples collected between 2002 and 2005 in the United States. Plasmid-positive strains were isolated from beef, chicken, turkey and pork, and were found in samples from the following states: California, Colorado, Connecticut, Georgia, Maryland, Minnesota, New Mexico, New York and Oregon. Our studies reveal that this common plasmid backbone is broadly disseminated among MDR zoonotic pathogens associated with agriculture. This reservoir of mobile resistance determinants has the potential to disseminate to Y. pestis and other human and zoonotic bacterial pathogens and therefore represents a significant public health concern.

Decabromodiphenyl ether (BDE 209), the major congener in the high volume industrial flame retardant mixture "DecaBDE", has recently been shown to be metabolized by carp. To further explore this phenomenon, juvenile rainbow trout were exposed to BDE 209 via the diet for a five month period. Analysis of the whole body homogenate, liver, serum, and intestinal tissues revealed that BDE 209 accumulated in rainbow trout tissues and was most concentrated in the liver. In addition to BDE 209, several hepta-, octa-, and nonaBDE congeners also accumulated in rainbow trout tissues over the same period as a result of BDE 209 debromination. Based on the total body burden of the hepta- through decaBDE congeners, uptake of BDE 209 was estimated at 3.2%. Congener profiles were different among whole body homogenate, liver, and serum, with the whole body homogenates having a greater contribution of the debrominated biotransformation products. Extracts of the rainbow trout whole body homogenates were compared with extracts from a previous experiment with common carp. This comparison revealed that BDE 202 (2,2',3,3',5,5',6,6'-octabromodiphenyl ether) was a dominant debromination product in both studies. To determine whether the observed debromination was metabolically driven, liver microsomal fractions were prepared from both common carp and rainbow trout. Analysis of the microsomal fractions following incubation with BDE 209 revealed that rainbow trout biotransformed as much as 22% of the BDE 209 mass, primarily to octa- and nonaBDE congeners. In contrast, carp liver microsomes biotransformed up to 65% of the BDE 209 mass, primarily down to hexaBDE congeners. These microsomal incubations confirm a metabolic pathway for BDE 209 debromination.

Males and females often differ in their fitness optima for shared traits that have a shared genetic basis, leading to sexual conflict. Morphologically differentiated sex chromosomes can resolve this conflict and protect sexually antagonistic variation, but they accumulate deleterious mutations. However, how sexual conflict is resolved in species that lack differentiated sex chromosomes is largely unknown. Here we present a chromosome-anchored genome assembly for rainbow trout (Oncorhynchus mykiss) and characterize a 55-Mb double-inversion supergene that mediates sex-specific migratory tendency through sex-dependent dominance reversal, an alternative mechanism for resolving sexual conflict. The double inversion contains key photosensory, circadian rhythm, adiposity and sex-related genes and displays a latitudinal frequency cline, indicating environmentally dependent selection. Our results show sex-dependent dominance reversal across a large autosomal supergene, a mechanism for sexual conflict resolution capable of protecting sexually antagonistic variation while avoiding the homozygous lethality and deleterious mutations associated with typical heteromorphic sex chromosomes.

The development and refinement of amphibian medicine comprise an ongoing science that reflects the unique life history of these animals and our growing knowledge of amphibian diseases. Amphibians are notoriously fastidious in terms of captive care requirements, and the majority of diseases of amphibians maintained in captivity will relate directly or indirectly to husbandry and management. Investigators have described many infectious and noninfectious diseases that occur among various species of captive and wild amphibians, and there is considerable overlap in the diseases of captive versus free-ranging populations. In this article, some of the more commonly reported infectious and noninfectious diseases as well as their etiological agents and causative factors are reviewed. Some of the more common amphibian diseases with bacterial etiologies include bacterial dermatosepticemia or "red leg syndrome," flavobacteriosis, mycobacteriosis, and chlamydiosis. The most common viral diseases of amphibians are caused by the ranaviruses, which have an impact on many species of anurans and caudates. Mycotic and mycotic-like organisms cause a number of diseases among amphibians, including chytridiomycosis, zygomycoses, chromomycoses, saprolegniasis, and ichthyophoniasis. Protozoan parasites of amphibians include a variety of amoeba, ciliates, flagellates, and sporozoans Common metazoan parasites include various myxozoans, helminths (particularly trematodes and nematodes), and arthropods. Commonly encountered noninfectious disease etiologies for amphibians include neoplasia, absolute or specific nutritional deficiencies or overloads, chemical toxicities, and inadequate husbandry or environmental management.

Advancing the production efficiency and profitability of aquaculture is dependent upon the ability to utilize a diverse array of genetic resources. The ultimate goals of aquaculture genomics, genetics and breeding research are to enhance aquaculture production efficiency, sustainability, product quality, and profitability in support of the commercial sector and for the benefit of consumers. In order to achieve these goals, it is important to understand the genomic structure and organization of aquaculture species, and their genomic and phenomic variations, as well as the genetic basis of traits and their interrelationships. In addition, it is also important to understand the mechanisms of regulation and evolutionary conservation at the levels of genome, transcriptome, proteome, epigenome, and systems biology. With genomic information and information between the genomes and phenomes, technologies for marker/causal mutation-assisted selection, genome selection, and genome editing can be developed for applications in aquaculture. A set of genomic tools and resources must be made available including reference genome sequences and their annotations (including coding and non-coding regulatory elements), genome-wide polymorphic markers, efficient genotyping platforms, high-density and high-resolution linkage maps, and transcriptome resources including non-coding transcripts. Genomic and genetic control of important performance and production traits, such as disease resistance, feed conversion efficiency, growth rate, processing yield, behaviour, reproductive characteristics, and tolerance to environmental stressors like low dissolved oxygen, high or low water temperature and salinity, must be understood. QTL need to be identified, validated across strains, lines and populations, and their mechanisms of control understood. Causal gene(s) need to be identified. Genetic and epigenetic regulation of important aquaculture traits need to be determined, and technologies for marker-assisted selection, causal gene/mutation-assisted selection, genome selection, and genome editing using CRISPR and other technologies must be developed, demonstrated with applicability, and application to aquaculture industries.Major progress has been made in aquaculture genomics for dozens of fish and shellfish species including the development of genetic linkage maps, physical maps, microarrays, single nucleotide polymorphism (SNP) arrays, transcriptome databases and various stages of genome reference sequences. This paper provides a general review of the current status, challenges and future research needs of aquaculture genomics, genetics, and breeding, with a focus on major aquaculture species in the United States: catfish, rainbow trout, Atlantic salmon, tilapia, striped bass, oysters, and shrimp. While the overall research priorities and the practical goals are similar across various aquaculture species, the current status in each species should dictate the next priority areas within the species. This paper is an output of the USDA Workshop for Aquaculture Genomics, Genetics, and Breeding held in late March 2016 in Auburn, Alabama, with participants from all parts of the United States.

As global aquaculture fish production continues to expand, an improved understanding of how environmental factors interact in fish health and production is needed. Significant advances have been made toward economical alternatives to costly fishmeal-based diets, such as grain-based formulations, and toward defining the effect of rearing density on fish health and production. Little research, however, has examined the effects of fishmeal- and grain-based diets in combination with alterations in rearing density. Moreover, it is unknown whether interactions between rearing density and diet impact the composition of the fish intestinal microbiota, which might in turn impact fish health and production. We fed aquacultured adult rainbow trout (Oncorhynchus mykiss) fishmeal- or grain-based diets, reared them under high- or low-density conditions for 10 months in a single aquaculture facility, and evaluated individual fish growth, production, fin indices, and intestinal microbiota composition using 16S rRNA gene sequencing. We found that the intestinal microbiotas were dominated by a shared core microbiota consisting of 52 bacterial lineages observed across all individuals, diets, and rearing densities. Variations in diet and rearing density resulted in only minor changes in intestinal microbiota composition despite significant effects of these variables on fish growth, performance, fillet quality, and welfare. Significant interactions between diet and rearing density were observed only in evaluations of fin indices and the relative abundance of the bacterial genus Staphylococcus. These results demonstrate that aquacultured rainbow trout can achieve remarkable consistency in intestinal microbiota composition and suggest the possibility of developing novel aquaculture strategies without overtly altering intestinal microbiota composition.

BACKGROUND: Previously, we have shown that bacterial cold water disease (BCWD) resistance in rainbow trout can be improved using traditional family-based selection, but progress has been limited to exploiting only between-family genetic variation. Genomic selection (GS) is a new alternative that enables exploitation of within-family genetic variation. METHODS: We compared three GS models [single-step genomic best linear unbiased prediction (ssGBLUP), weighted ssGBLUP (wssGBLUP), and BayesB] to predict genomic-enabled breeding values (GEBV) for BCWD resistance in a commercial rainbow trout population, and compared the accuracy of GEBV to traditional estimates of breeding values (EBV) from a pedigree-based BLUP (P-BLUP) model. We also assessed the impact of sampling design on the accuracy of GEBV predictions. For these comparisons, we used BCWD survival phenotypes recorded on 7893 fish from 102 families, of which 1473 fish from 50 families had genotypes [57 K single nucleotide polymorphism (SNP) array]. Naïve siblings of the training fish (n = 930 testing fish) were genotyped to predict their GEBV and mated to produce 138 progeny testing families. In the following generation, 9968 progeny were phenotyped to empirically assess the accuracy of GEBV predictions made on their non-phenotyped parents. RESULTS: The accuracy of GEBV from all tested GS models were substantially higher than the P-BLUP model EBV. The highest increase in accuracy relative to the P-BLUP model was achieved with BayesB (97.2 to 108.8%), followed by wssGBLUP at iteration 2 (94.4 to 97.1%) and 3 (88.9 to 91.2%) and ssGBLUP (83.3 to 85.3%). Reducing the training sample size to n = ~1000 had no negative impact on the accuracy (0.67 to 0.72), but with n = ~500 the accuracy dropped to 0.53 to 0.61 if the training and testing fish were full-sibs, and even substantially lower, to 0.22 to 0.25, when they were not full-sibs. CONCLUSIONS: Using progeny performance data, we showed that the accuracy of genomic predictions is substantially higher than estimates obtained from the traditional pedigree-based BLUP model for BCWD resistance. Overall, we found that using a much smaller training sample size compared to similar studies in livestock, GS can substantially improve the selection accuracy and genetic gains for this trait in a commercial rainbow trout breeding population.

Androgenetic doubled haploid progeny produced from a cross between the Oregon State University and Arlee clonal rainbow trout (Oncorhynchus mykiss) lines, used for a previous published rainbow trout map, were used to update the map with the addition of more amplified fragment length polymorphic (AFLP) markers, microsatellites, type I and allozyme markers. We have added more than 900 markers, bringing the total number to 1359 genetic markers and the sex phenotype including 799 EcoRI AFLPs, 174 PstI AFLPs, 226 microsatellites, 72 VNTR, 38 SINE markers, 29 known genes, 12 minisatellites, five RAPDs, and four allozymes. Thirty major linkage groups were identified. Synteny of linkage groups in our map with the outcrossed microsatellite map has been established for all except one linkage group in this doubled haploid cross. Putative homeologous relationships among linkage groups, resulting from the autotetraploid nature of the salmonid genome, have been revealed based on the placement of duplicated microsatellites and type I loci.

Salmonid fishes exhibit extensive local adaptations owing to abundant environmental variation and precise natal homing. This extensive local adaptation makes conservation and restoration of salmonids a challenge. For example, defining unambiguous units of conservation is difficult, and restoration attempts often fail owing to inadequate adaptive matching of translocated populations. A better understanding of the genetic architecture of local adaptation in salmonids could provide valuable information to assist in conserving and restoring natural populations of these important species. Here, we use a combination of laboratory crosses and next-generation sequencing to investigate the genetic architecture of the parallel adaptation of rapid development rate in two geographically and genetically distant populations of rainbow trout (Oncorhynchus mykiss). Strikingly, we find that not only is a parallel genetic mechanism used but that a conserved haplotype is responsible for this intriguing adaptation. The repeated use of adaptive genetic variation across distant geographical areas could be a general theme in salmonids and have important implications for conservation and restoration.

Multidrug resistance (MDR) plasmids belonging to the IncA/C plasmid family are widely distributed among Salmonella and other enterobacterial isolates from agricultural sources and have, at least once, also been identified in a drug-resistant Yersinia pestis isolate (IP275) from Madagascar. Here, we present the complete plasmid sequences of the IncA/C reference plasmid pRA1 (143,963 bp), isolated in 1971 from the fish pathogen Aeromonas hydrophila, and of the cryptic IncA/C plasmid pRAx (49,763 bp), isolated from Escherichia coli transconjugant D7-3, which was obtained through pRA1 transfer in 1980. Using comparative sequence analysis of pRA1 and pRAx with recent members of the IncA/C plasmid family, we show that both plasmids provide novel insights into the evolution of the IncA/C MDR plasmid family and the minimal machinery necessary for stable IncA/C plasmid maintenance. Our results indicate that recent members of the IncA/C plasmid family evolved from a common ancestor, similar in composition to pRA1, through stepwise integration of horizontally acquired resistance gene arrays into a conserved plasmid backbone. Phylogenetic comparisons predict type IV secretion-like conjugative transfer operons encoded on the shared plasmid backbones to be closely related to a group of integrating conjugative elements, which use conjugative transfer for horizontal propagation but stably integrate into the host chromosome during vegetative growth. A hipAB toxin-antitoxin gene cluster found on pRA1, which in Escherichia coli is involved in the formation of persister cell subpopulations, suggests persistence as an early broad-spectrum antimicrobial resistance mechanism in the evolution of IncA/C resistance plasmids.

BACKGROUND: Fast, efficiently growing animals have increased protein synthesis and/or reduced protein degradation relative to slow, inefficiently growing animals. Consequently, minimizing the energetic cost of protein turnover is a strategic goal for enhancing animal growth. Characterization of gene expression profiles associated with protein turnover would allow us to identify genes that could potentially be used as molecular biomarkers to select for germplasm with improved protein accretion. RESULTS: We evaluated changes in hepatic global gene expression in response to 3-week starvation in rainbow trout (Oncorhynchus mykiss). Microarray analysis revealed a coordinated, down-regulated expression of protein biosynthesis genes in starved fish. In addition, the expression of genes involved in lipid metabolism/transport, aerobic respiration, blood functions and immune response were decreased in response to starvation. However, the microarray approach did not show a significant increase of gene expression in protein catabolic pathways. Further studies, using real-time PCR and enzyme activity assays, were performed to investigate the expression of genes involved in the major proteolytic pathways including calpains, the multi-catalytic proteasome and cathepsins. Starvation reduced mRNA expression of the calpain inhibitor, calpastatin long isoform (CAST-L), with a subsequent increase in the calpain catalytic activity. In addition, starvation caused a slight but significant increase in 20S proteasome activity without affecting mRNA levels of the proteasome genes. Neither the mRNA levels nor the activities of cathepsin D and L were affected by starvation. CONCLUSION: These results suggest a significant role of calpain and 20S proteasome pathways in protein mobilization as a source of energy during fasting and a potential association of the CAST-L gene with fish protein accretion.

BACKGROUND: Salmonids are of interest because of their relatively recent genome duplication, and their extensive use in wild fisheries and aquaculture. A comprehensive gene list and a comparison of genes in some of the different species provide valuable genomic information for one of the most widely studied groups of fish. RESULTS: 298,304 expressed sequence tags (ESTs) from Atlantic salmon (69% of the total), 11,664 chinook, 10,813 sockeye, 10,051 brook trout, 10,975 grayling, 8,630 lake whitefish, and 3,624 northern pike ESTs were obtained in this study and have been deposited into the public databases. Contigs were built and putative full-length Atlantic salmon clones have been identified. A database containing ESTs, assemblies, consensus sequences, open reading frames, gene predictions and putative annotation is available. The overall similarity between Atlantic salmon ESTs and those of rainbow trout, chinook, sockeye, brook trout, grayling, lake whitefish, northern pike and rainbow smelt is 93.4, 94.2, 94.6, 94.4, 92.5, 91.7, 89.6, and 86.2% respectively. An analysis of 78 transcript sets show Salmo as a sister group to Oncorhynchus and Salvelinus within Salmoninae, and Thymallinae as a sister group to Salmoninae and Coregoninae within Salmonidae. Extensive gene duplication is consistent with a genome duplication in the common ancestor of salmonids. Using all of the available EST data, a new expanded salmonid cDNA microarray of 32,000 features was created. Cross-species hybridizations to this cDNA microarray indicate that this resource will be useful for studies of all 68 salmonid species. CONCLUSION: An extensive collection and analysis of salmonid RNA putative transcripts indicate that Pacific salmon, Atlantic salmon and charr are 94-96% similar while the more distant whitefish, grayling, pike and smelt are 93, 92, 89 and 86% similar to salmon. The salmonid transcriptome reveals a complex history of gene duplication that is consistent with an ancestral salmonid genome duplication hypothesis. Genome resources, including a new 32 K microarray, provide valuable new tools to study salmonids.

Poultry litter applications to land have been based on crop N requirements, resulting in application of P in excess of plant requirements, which may cause degradation of water quality in the Chesapeake Bay watershed. The effect of litter source (the Delmarva Peninsula and Moorefield, West Virginia) and composting of poultry litter on N mineralization and availability of P in two soil types (sandy loam and silt loam) was determined in a controlled environment for 120 d. Nitrogen mineralization (percent total organic N converted to inorganic nitrogen) rates were higher for fresh litter (range of 42 to 64%) than composted litter (range of 1 to 9%). The N mineralization rate of fresh litter from the Delmarva Peninsula was consistently lower than the fresh litter from Moorefield, WV. The N mineralization rate of composted litter from either source was not significantly different for each soil type (7 to 9% in sandy loam and 1 to 5% in silt loam) even though composting conditions were completely different at the two composting facilities. Litter source had a large effect on N mineralization rates of fresh but not composted poultry litter. Composting yielded a more predictable and reliable source of mineralizable N than fresh litter. Water-extractable phosphorus (WEP) was similar in soils amended with composted litter from WV and fresh litter from both sources (approximately 10 to 25 and 2 to 14 mg P kg(-1) for sandy loam and silt loam, respectively). Mehlich 1-extractable phosphorus (MEP) was similar in soils amended with WV fresh litter and composted litter from both sources (approximately 100 to 140 and 60 to 90 mg P kg(-1) for sandy loam and silt loam, respectively). These results suggest that the composting process did not consistently reduce WEP and MEP, and P can be as available in composted poultry litter as in fresh poultry litter.

Fast growth is an important and highly desired trait, which affects the profitability of food animal production, with feed costs accounting for the largest proportion of production costs. Traditional phenotype-based selection is typically used to select for growth traits; however, genetic improvement is slow over generations. Single nucleotide polymorphisms (SNPs) explain 90% of the genetic differences between individuals; therefore, they are most suitable for genetic evaluation and strategies that employ molecular genetics for selective breeding. SNPs found within or near a coding sequence are of particular interest because they are more likely to alter the biological function of a protein. We aimed to use SNPs to identify markers and genes associated with genetic variation in growth. RNA-Seq whole-transcriptome analysis of pooled cDNA samples from a population of rainbow trout selected for improved growth versus unselected genetic cohorts (10 fish from 1 full-sib family each) identified SNP markers associated with growth-rate. The allelic imbalances (the ratio between the allele frequencies of the fast growing sample and that of the slow growing sample) were considered at scores >5.0 as an amplification and <0.2 as loss of heterozygosity. A subset of SNPs (n = 54) were validated and evaluated for association with growth traits in 778 individuals of a three-generation parent/offspring panel representing 40 families. Twenty-two SNP markers and one mitochondrial haplotype were significantly associated with growth traits. Polymorphism of 48 of the markers was confirmed in other commercially important aquaculture stocks. Many markers were clustered into genes of metabolic energy production pathways and are suitable candidates for genetic selection. The study demonstrates that RNA-Seq at low sequence coverage of divergent populations is a fast and effective means of identifying SNPs, with allelic imbalances between phenotypes. This technique is suitable for marker development in non-model species lacking complete and well-annotated genome reference sequences.

Selective breeding of animals for increased innate resistance offers an attractive strategy to control disease in agriculture. However, this approach is limited by an incomplete knowledge of the heritability, duration, and mechanism(s) of resistance, as well as the impact of selection on the immune response to unrelated pathogens. Herein, as part of a rainbow trout broodstock improvement program, we evaluated factors involved in resistance against a bacterial disease agent, Flavobacterium psychrophilum. In 2005, 71 full-sibling crosses, weighing an average of 2.4 g, were screened, and resistant and susceptible crosses were identified. Naive cohorts were evaluated at 10 and 800 g in size, and most maintained their original relative resistant or susceptible phenotypes, indicating that these traits were stable as size increased >300-fold. During the course of these studies, we observed that the normalized spleen weights of the resistant fish crosses were greater than those of the susceptible fish crosses. To test for direct association, we determined the spleen-somatic index of 103 fish crosses; created high, medium, and low spleen-index groups; and determined survival following challenge with F. psychrophilum or Yersinia ruckeri. Consistent with our previous observations, trout with larger spleen indices were significantly more resistant to F. psychrophilum challenge; however, this result was pathogen-specific, as there was no correlation of spleen size with survival following Y. ruckeri challenge. To our knowledge, this is the first report of a positive association between spleen size and disease resistance in a teleost fish. Further evaluation of spleen index as an indirect measure of disease resistance is warranted.

BACKGROUND: 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. RESULTS: 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. CONCLUSIONS: 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.