Population transcriptomics of life-history variation in the genus Salmo

Assessing Outlier Probabilities in Transcriptomics Data When Evaluating a Classifier

Outliers in the training or test set used to fit and evaluate a classifier on transcriptomics data can considerably change the estimated performance of the model. Hence, an either too weak or a too optimistic accuracy is then reported and the estimated model performance cannot be reproduced on independent data. It is then also doubtful whether a classifier qualifies for clinical usage. We estimate classifier performances in simulated gene expression data with artificial outliers and in two real-world datasets. As a new approach, we use two outlier detection methods within a bootstrap procedure to estimate the outlier probability for each sample and evaluate classifiers before and after outlier removal by means of cross-validation. We found that the removal of outliers changed the classification performance notably. For the most part, removing outliers improved the classification results. Taking into account the fact that there are various, sometimes unclear reasons for a sample to be an outlier, we strongly advocate to always report the performance of a transcriptomics classifier with and without outliers in training and test data. This provides a more diverse picture of a classifier's performance and prevents reporting models that later turn out to be not applicable for clinical diagnoses.

Alternative migratory tactics in brown trout (Salmo trutta) are underpinned by divergent regulation of metabolic but not neurological genes

The occurrence of alternative morphs within populations is common, but the underlying molecular mechanisms remain poorly understood. Many animals, for example, exhibit facultative migration, where two or more alternative migratory tactics (AMTs) coexist within populations. In certain salmonid species, some individuals remain in natal rivers all their lives, while others (in particular, females) migrate to sea for a period of marine growth. Here, we performed transcriptional profiling ("RNA-seq") of the brain and liver of male and female brown trout to understand the genes and processes that differentiate between migratory and residency morphs (AMT-associated genes) and how they may differ in expression between the sexes. We found tissue-specific differences with a greater number of genes expressed differentially in the liver (n = 867 genes) compared with the brain (n = 10) between the morphs. Genes with increased expression in resident livers were enriched for Gene Ontology terms associated with metabolic processes, highlighting key molecular-genetic pathways underlying the energetic requirements associated with divergent migratory tactics. In contrast, smolt-biased genes were enriched for biological processes such as response to cytokines, suggestive of possible immune function differences between smolts and residents. Finally, we identified evidence of sex-biased gene expression for AMT-associated genes in the liver (n = 12) but not the brain. Collectively, our results provide insights into tissue-specific gene expression underlying the production of alternative life histories within and between the sexes, and point toward a key role for metabolic processes in the liver in mediating divergent physiological trajectories of migrants versus residents.

Investigating Diadromy in Fishes and Its Loss in an -Omics Era

Diadromy, the predictable movements of individuals between marine and freshwater environments, is biogeographically and phylogenetically widespread across fishes. Thus, despite the high energetic and potential fitness costs involved in moving between distinct environments, diadromy appears to be an effective life history strategy. Yet, the origin and molecular mechanisms that underpin this migratory behavior are not fully understood. In this review, we aim first to summarize what is known about diadromy in fishes; this includes the phylogenetic relationship among diadromous species, a description of the main hypotheses regarding its origin, and a discussion of the presence of non-migratory populations within diadromous species. Second, we discuss how recent research based on -omics approaches (chiefly genomics, transcriptomics, and epigenomics) is beginning to provide answers to questions on the genetic bases and origin(s) of diadromy. Finally, we suggest future directions for -omics research that can help tackle questions on the evolution of diadromy.

Anadromy, potamodromy and residency in brown trout Salmo trutta: the role of genes and the environment

Brown trout Salmo trutta is endemic to Europe, western Asia and north-western Africa; it is a prominent member of freshwater and coastal marine fish faunas. The species shows two resident (river-resident, lake-resident) and three main facultative migratory life histories (downstream-upstream within a river system, fluvial-adfluvial potamodromous; to and from a lake, lacustrine-adfluvial (inlet) or allacustrine (outlet) potamodromous; to and from the sea, anadromous). River-residency v. migration is a balance between enhanced feeding and thus growth advantages of migration to a particular habitat v. the costs of potentially greater mortality and energy expenditure. Fluvial-adfluvial migration usually has less feeding improvement, but less mortality risk, than lacustrine-adfluvial or allacustrine and anadromous, but the latter vary among catchments as to which is favoured. Indirect evidence suggests that around 50% of the variability in S. trutta migration v. residency, among individuals within a population, is due to genetic variance. This dichotomous decision can best be explained by the threshold-trait model of quantitative genetics. Thus, an individual's physiological condition (e.g., energy status) as regulated by environmental factors, genes and non-genetic parental effects, acts as the cue. The magnitude of this cue relative to a genetically predetermined individual threshold, governs whether it will migrate or sexually mature as a river-resident. This decision threshold occurs early in life and, if the choice is to migrate, a second threshold probably follows determining the age and timing of migration. Migration destination (mainstem river, lake, or sea) also appears to be genetically programmed. Decisions to migrate and ultimate destination result in a number of subsequent consequential changes such as parr-smolt transformation, sexual maturity and return migration. Strong associations with one or a few genes have been found for most aspects of the migratory syndrome and indirect evidence supports genetic involvement in all parts. Thus, migratory and resident life histories potentially evolve as a result of natural and anthropogenic environmental changes, which alter relative survival and reproduction. Knowledge of genetic determinants of the various components of migration in S. trutta lags substantially behind that of Oncorhynchus mykiss and other salmonines. Identification of genetic markers linked to migration components and especially to the migration-residency decision, is a prerequisite for facilitating detailed empirical studies. In order to predict effectively, through modelling, the effects of environmental changes, quantification of the relative fitness of different migratory traits and of their heritabilities, across a range of environmental conditions, is also urgently required in the face of the increasing pace of such changes.

Role of genomics and transcriptomics in selection of reintroduction source populations

The use and importance of reintroduction as a conservation tool to return a species to its historical range from which it has been extirpated will increase as climate change and human development accelerate habitat loss and population extinctions. Although the number of reintroduction attempts has increased rapidly over the past 2 decades, the success rate is generally low. As a result of population differences in fitness-related traits and divergent responses to environmental stresses, population performance upon reintroduction is highly variable, and it is generally agreed that selecting an appropriate source population is a critical component of a successful reintroduction. Conservation genomics is an emerging field that addresses long-standing challenges in conservation, and the potential for using novel molecular genetic approaches to inform and improve conservation efforts is high. Because the successful establishment and persistence of reintroduced populations is highly dependent on the functional genetic variation and environmental stress tolerance of the source population, we propose the application of conservation genomics and transcriptomics to guide reintroduction practices. Specifically, we propose using genome-wide functional loci to estimate genetic variation of source populations. This estimate can then be used to predict the potential for adaptation. We also propose using transcriptional profiling to measure the expression response of fitness-related genes to environmental stresses as a proxy for acclimation (tolerance) capacity. Appropriate application of conservation genomics and transcriptomics has the potential to dramatically enhance reintroduction success in a time of rapidly declining biodiversity and accelerating environmental change.

Exploring transcriptomic diversity in muscle revealed that cellular signaling pathways mainly differentiate five Western porcine breeds

Background

Among transcriptomic studies, those comparing species or populations can increase our understanding of the impact of the evolutionary forces on the differentiation of populations. A particular situation is the one of short evolution time with breeds of a domesticated species that underwent strong selective pressures. In this study, the gene expression diversity across five pig breeds has been explored in muscle. Samples came from: 24 Duroc, 33 Landrace, 41 Large White dam line, 10 Large White sire line and 39 Piétrain. From these animals, 147 muscle samples obtained at slaughter were analyzed using the porcine Agilent 44 K v1 microarray.

Results

A total of 12,358 genes were identified as expressed in muscle after normalization and 1,703 genes were declared differential for at least one breed (FDR < 0.001). The functional analysis highlighted that gene expression diversity is mainly linked to cellular signaling pathways such as the PI3K (phosphoinositide 3-kinase) pathway. The PI3K pathway is known to be involved in the control of development of the skeletal muscle mass by affecting extracellular matrix - receptor interactions, regulation of actin cytoskeleton pathways and some metabolic functions. This study also highlighted 228 spots (171 unique genes) that differentiate the breeds from each other. A common subgroup of 15 genes selected by three statistical methods was able to differentiate Duroc, Large White and Piétrain breeds.

Conclusions

This study on transcriptomic differentiation across Western pig breeds highlighted a global picture: mainly signaling pathways were affected. This result is consistent with the selection objective of increasing muscle mass. These transcriptional changes may indicate selection pressure or simply breed differences which may be driven by human selection. Further work aiming at comparing genetic and transcriptomic diversities would further increase our understanding of the consequences of human impact on livestock species.

Electronic supplementary material

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

Divergence in expression of candidate genes for the smoltification process between juvenile resident rainbow and anadromous steelhead trout

Rainbow and steelhead trout (Oncorhynchus mykiss), among other salmonid fishes, exhibit tremendous life history diversity, foremost of which is variation in migratory propensity. While some individuals possess the ability to undertake an anadromous marine migration, others remain resident in freshwater throughout their life cycle. Those that will migrate undergo tremendous physiological, morphological, and behavioral transformations in a process called smoltification which transitions freshwater-adapted parr to marine-adapted smolts. While the behavior, ecology, and physiology of smoltification are well described, our understanding of the proximate genetic mechanisms that trigger the process are not well known. Quantitative genetic analyses have identified several genomic regions associated with smoltification and migration-related traits within this species. Here we investigate the divergence in gene expression of 18 functional and positional candidate genes for the smoltification process in the brain, gill, and liver tissues of migratory smolts, resident parr, and precocious mature male trout at the developmental stage of out-migration. Our analysis reveals several genes differentially expressed between life history classes and validates the candidate nature of several genes in the parr-smolt transformation including Clock1α, FSHβ, GR, GH2, GHR1, GHR2, NDK7, p53, SC6a7, Taldo1, THRα, THRβ, and Vdac2.

There and back again: migration in freshwater fishes

Animal migration is an amazing phenomenon that has fascinated humans for long. Many freshwater fishes also show remarkable migrations, whereof the spectacular mass migrations of salmonids from the spawning streams are the most well known and well studied. However, recent studies have shown that migration occurs in a range of freshwater fish taxa from many different habitats. In this review we focus on the causes and consequences of migration in freshwater fishes. We start with an introduction of concepts and categories of migration, and then address the evolutionary causes that drive individuals to make these migratory journeys. The basis for the decision of an individual fish to migrate or stay resident is an evaluation of the costs and benefits of different strategies to maximize its lifetime reproductive effort. We provide examples by discussing our own work on the causes behind seasonal migration in a cyprinid fish, roach (Rutilus rutilus (L., 1758)), within this framework. We then highlight different adaptations that allow fish to migrate over sometimes vast journeys across space, including capacity for orientation, osmoregulation, and efficient energy expenditure. Following this we consider the consequences of migration in freshwater fish from ecological, evolutionary, and conservation perspectives, and finally, we detail some of the recent developments in the methodologies used to collect data on fish migration and how these could be used in future research.

The physiological basis of the migration continuum in brown trout (Salmo trutta)

Partial migration is common in many animal taxa; however, the physiological variation underpinning migration strategies remains poorly understood. Among salmonid fishes, brown trout (Salmo trutta) is one of the species that exhibits the most complex variation in sympatric migration strategies, expressed as a migration continuum, ranging from residency to anadromy. In looking at brown trout, our objective with this study was to test the hypothesis that variation in migration strategies is underpinned by physiological variation. Prior to migration, physiological samples were taken from fish in the stream and then released at the capture site. Using telemetry, we subsequently classified fish as resident, short-distance migrants (potamodromous), or long-distance migrants (potentially anadromous). Our results revealed that fish belonging to the resident strategy differed from those exhibiting any of the two migratory strategies. Gill Na,K-ATPase activity, condition factor, and indicators of nutritional status suggested that trout from the two migratory strategies were smoltified and energetically depleted before leaving the stream, compared to those in the resident strategy. The trout belonging to the two migratory strategies were generally similar; however, lower triacylglycerides levels in the short-distance migrants indicated that they were more lipid depleted prior to migration compared with the long-distance migrants. In the context of migration cost, we suggest that additional lipid depletion makes migrants more inclined to terminate migration at the first given feeding opportunity, whereas individuals that are less lipid depleted will migrate farther. Collectively, our data suggest that the energetic state of individual fish provides a possible mechanism underpinning the migration continuum in brown trout.

Local adaptation at the transcriptome level in brown trout: evidence from early life history temperature genomic reaction norms

Local adaptation and its underlying molecular basis has long been a key focus in evolutionary biology. There has recently been increased interest in the evolutionary role of plasticity and the molecular mechanisms underlying local adaptation. Using transcriptome analysis, we assessed differences in gene expression profiles for three brown trout (Salmo trutta) populations, one resident and two anadromous, experiencing different temperature regimes in the wild. The study was based on an F2 generation raised in a common garden setting. A previous study of the F1 generation revealed different reaction norms and significantly higher Q_ST than F_ST among populations for two early life-history traits. In the present study we investigated if genomic reaction norm patterns were also present at the transcriptome level. Eggs from the three populations were incubated at two temperatures (5 and 8 degrees C) representing conditions encountered in the local environments. Global gene expression for fry at the stage of first feeding was analysed using a 32k cDNA microarray. The results revealed differences in gene expression between populations and temperatures and population × temperature interactions, the latter indicating locally adapted reaction norms. Moreover, the reaction norms paralleled those observed previously at early life-history traits. We identified 90 cDNA clones among the genes with an interaction effect that were differently expressed between the ecologically divergent populations. These included genes involved in immune- and stress response. We observed less plasticity in the resident as compared to the anadromous populations, possibly reflecting that the degree of environmental heterogeneity encountered by individuals throughout their life cycle will select for variable level of phenotypic plasticity at the transcriptome level. Our study demonstrates the usefulness of transcriptome approaches to identify genes with different temperature reaction norms. The responses observed suggest that populations may vary in their susceptibility to climate change.

Genetic architecture of migration-related traits in rainbow and steelhead trout, Oncorhynchus mykiss

Although migration plays a critical role in the evolution and diversification of species, relatively little is known of the genetic architecture underlying this life history in any species. Rainbow and steelhead trout (Oncorhynchus mykiss) naturally segregate for both resident and migratory life-history types, respectively, as do other members of the salmonid family of fishes. Using an experimental cross derived from wild resident rainbow and wild migratory steelhead trout from Southeast Alaska and high throughput restriction-site associated DNA (RAD) tag sequencing, we perform a quantitative trait locus (QTL) analysis to identify the number, position, and relative contribution of genetic effects on a suite of 27 physiological and morphological traits associated with the migratory life history in this species. In total, 37 QTL are localized to 19 unique QTL positions, explaining 4-13.63% of the variation for 19 of the 27 migration-related traits measured. Two chromosomal positions, one on chromosome Omy12 and the other on Omy14 each harbor 7 QTL for migration-related traits, suggesting that these regions could harbor master genetic controls for the migratory life-history tactic in this species. Another QTL region on Omy5 has been implicated in several studies of adaptive life histories within this species and could represent another important locus underlying the migratory life history. We also evaluate whether loci identified in this out-crossed QTL study colocalize to genomic positions previously identified for associations with migration-related traits in a doubled haploid mapping family.

Brain transcriptome variation among behaviorally distinct strains of zebrafish (Danio rerio)

BACKGROUND

Domesticated animal populations often show profound reductions in predator avoidance and fear-related behavior compared to wild populations. These reductions are remarkably consistent and have been observed in a diverse array of taxa including fish, birds, and mammals. Experiments conducted in common environments indicate that these behavioral differences have a genetic basis. In this study, we quantified differences in fear-related behavior between wild and domesticated zebrafish strains and used microarray analysis to identify genes that may be associated with this variation.

RESULTS

Compared to wild zebrafish, domesticated zebrafish spent more time near the water surface and were more likely to occupy the front of the aquarium nearest a human observer. Microarray analysis of the brain transcriptome identified high levels of population variation in gene expression, with 1,749 genes significantly differentially expressed among populations. Genes that varied among populations belonged to functional categories that included DNA repair, DNA photolyase activity, response to light stimulus, neuron development and axon guidance, cell death, iron-binding, chromatin reorganization, and homeobox genes. Comparatively fewer genes (112) differed between domesticated and wild strains with notable genes including gpr177 (wntless), selenoprotein P1a, synaptophysin and synaptoporin, and acyl-CoA binding domain containing proteins (acbd3 and acbd4).

CONCLUSIONS

Microarray analysis identified a large number of genes that differed among zebrafish populations and may underlie behavioral domestication. Comparisons with similar microarray studies of domestication in rainbow trout and canids identified sixteen evolutionarily or functionally related genes that may represent components of shared molecular mechanisms underlying convergent behavioral evolution during vertebrate domestication. However, this conclusion must be tempered by limitations associated with comparisons among microarray studies and the low level of population-level replication inherent to these studies.

Ecological transcriptomics of lake-type and riverine sockeye salmon (Oncorhynchus nerka)

BACKGROUND

There are a growing number of genomes sequenced with tentative functions assigned to a large proportion of the individual genes. Model organisms in laboratory settings form the basis for the assignment of gene function, and the ecological context of gene function is lacking. This work addresses this shortcoming by investigating expressed genes of sockeye salmon (Oncorhynchus nerka) muscle tissue. We compared morphology and gene expression in natural juvenile sockeye populations related to river and lake habitats. Based on previously documented divergent morphology, feeding strategy, and predation in association with these distinct environments, we expect that burst swimming is favored in riverine population and continuous swimming is favored in lake-type population. In turn we predict that morphology and expressed genes promote burst swimming in riverine sockeye and continuous swimming in lake-type sockeye.

RESULTS

We found the riverine sockeye population had deep, robust bodies and lake-type had shallow, streamlined bodies. Gene expression patterns were measured using a 16 k microarray, discovering 141 genes with significant differential expression. Overall, the identity and function of these genes was consistent with our hypothesis. In addition, Gene Ontology (GO) enrichment analyses with a larger set of differentially expressed genes found the "biosynthesis" category enriched for the riverine population and the "metabolism" category enriched for the lake-type population.

CONCLUSIONS

This study provides a framework for understanding sockeye life history from a transcriptomic perspective and a starting point for more extensive, targeted studies determining the ecological context of genes.

Genomic investigation of year-long and multigenerational exposures of fathead minnow to the munitions compound RDX

We assessed the impacts of exposure to an environmentally representative concentration (0.83 mg/L) of the explosive cyclotrimethylenetrinitramine (RDX) on fathead minnows (Pimephales promelas) in one-year and multigenerational bioassays. In the one-year bioassay, impacts were assessed by statistical comparisons of females from breeding groups reared in control or RDX-exposure conditions. The RDX had no significant effect on gonadosomatic index or condition factor assayed at 1 d and at one, three, six, nine, and 12 months. The liver-somatic index was significantly increased versus controls only at the 12-month timepoint. RDX had no significant effect on live-prey capture rates, egg production, or fertilization. RDX caused minimal differential-transcript expression with no consistent discernable effect on gene-functional categories for either brain or liver tissues in the one-year exposure. In the multigenerational assay, the effects of acute (96 h) exposure to RDX were compared in fish reared to the F(2) generation in either control or RDX-exposure conditions. Enrichment of gene functions including neuroexcitatory glutamate metabolism, sensory signaling, and neurological development were observed comparing control-reared and RDX-reared fish. Our results indicated that exposure to RDX at a concentration representing the highest levels observed in the environment (0.83 mg/L) had limited impacts on genomic, individual, and population-level endpoints in fathead minnows in a one-year exposure. However, multigenerational exposures altered transcript expression related to neural development and function. Environ.

Identification of candidate genes and physiological pathways involved in gonad deformation in whitefish (Coregonus spp.) from Lake Thun, Switzerland

In 2000, fishermen reported the appearance of deformed reproductive organs in whitefish (Coregonus spp.) from Lake Thun, Switzerland. Despite intensive investigations, the causes of these abnormalities remain unknown. Using gene expression profiling, we sought to identify candidate genes and physiological processes possibly associated with the observed gonadal deformations, in order to gain insights into potential causes. Using in situ-synthesized oligonucleotide arrays, we compared the expression levels at 21,492 unique transcript probes in liver and head kidney tissue of male whitefish with deformed and normally developed gonads, respectively. The fish had been collected on spawning sites of two genetically distinct whitefish forms of Lake Thun. We contrasted the gene expression profiles of 56 individuals, i.e., 14 individuals of each phenotype and of each population. Gene-by-gene analysis revealed weak expression differences between normal and deformed fish, and only one gene, ictacalcin, was found to be up-regulated in head kidney tissue of deformed fish from both whitefish forms, However, this difference could not be confirmed with quantitative real-time qPCR. Enrichment analysis on the level of physiological processes revealed (i) the involvement of immune response genes in both tissues, particularly those linked to complement activation in the liver, (ii) proteolysis in the liver and (iii) GTPase activity and Ras protein signal transduction in the head kidney. In comparison with current literature, this gene expression pattern signals a chronic autoimmune disease in the testes. Based on the recent observations that gonad deformations are induced through feeding of zooplankton from Lake Thun we hypothesize that a xenobiotic accumulated in whitefish via the plankton triggering autoimmunity as the likely cause of gonad deformations. We propose several experimental strategies to verify or reject this hypothesis.

Extent and scale of local adaptation in salmonid fishes: review and meta-analysis

What is the extent and scale of local adaptation (LA)? How quickly does LA arise? And what is its underlying molecular basis? Our review and meta-analysis on salmonid fishes estimates the frequency of LA to be ∼55-70%, with local populations having a 1.2 times average fitness advantage relative to foreign populations or to their performance in new environments. Salmonid LA is evident at a variety of spatial scales (for example, few km to>1000 km) and can manifest itself quickly (6-30 generations). As the geographic scale between populations increases, LA is generally more frequent and stronger. Yet the extent of LA in salmonids does not appear to differ from that in other assessed taxa. Moreover, the frequency with which foreign salmonid populations outperform local populations (∼23-35%) suggests that drift, gene flow and plasticity often limit or mediate LA. The relatively few studies based on candidate gene and genomewide analyses have identified footprints of selection at both small and large geographical scales, likely reflecting the specific functional properties of loci and the associated selection regimes (for example, local niche partitioning, pathogens, parasites, photoperiodicity and seasonal timing). The molecular basis of LA in salmonids is still largely unknown, but differential expression at the same few genes is implicated in the convergent evolution of certain phenotypes. Collectively, future research will benefit from an integration of classical and molecular approaches to understand: (i) species differences and how they originate, (ii) variation in adaptation across scales, life stages, population sizes and environmental gradients, and (iii) evolutionary responses to human activities.

Gene expression analysis for the identification of selection and local adaptation in fishes

In recent years, variation in gene expression has been recognized as an important component of environmental adaptation in multiple model species, including a few fish species. There is, however, still little known about the genetic basis of adaptation in gene expression resulting from variation in the aquatic environment (e.g. temperature, salinity and oxygen) and the physiological effect and costs of such differences in gene expression. This review presents and discusses progress and pitfalls of applying gene expression analyses to fishes and suggests simple frameworks to get started with gene expression analysis. It is emphasized that well-planned gene expression studies can serve as an important tool for the identification of selection in local populations of fishes, even for non-traditional model species where limited genomic information is available. Recent studies focusing on gene expression variation among natural fish populations are reviewed, highlighting the latest applications that combine genetic evidence from neutral markers and gene expression data.

Rapid evolution of osmoregulatory function by modification of gene transcription in steelhead trout

Populations experiencing sudden environmental change must be capable of rapidly evolving to survive. Here we explore changes in gene transcription as a mechanism for rapid adaptation at four osmoregulatory genes (CFTR I, NaK ATPase1αa, NaK ATPase1αb and GHRII) in anadromous steelhead trout versus a derived land-locked population after 14 generations. Transcription was measured before and after a 24-h saltwater challenge in pure and reciprocal hybrid offspring of fish from both populations reared in a common environment for two generations. Significant differences between the landlocked and migratory populations were observed, particularly in fresh water at the NaK ATPase1αa and GHRII genes, indicating rapid evolutionary change, possibly associated with reduced energy expenditure in the landlocked lake system. Phenotypic divergence analysis (Q (ST)) shows that the observed transcriptional differences deviate from neutral expectations. Some reciprocal crosses exhibited anomalous transcription consistent with sex-linked epistatic or genetic imprinting effects. Our results highlight unpredictable phenotypic outcomes of hybridization among locally adapted populations and the need to exercise caution when interbreeding populations for conservation purposes.

Promise and progress in environmental genomics: a status report on the applications of gene expression-based microarray studies in ecologically relevant fish species

The advent of any new technology is typically met with great excitement. So it was a few years ago, when the combination of advances in sequencing technology and the development of microarray technology made measurements of global gene expression in ecologically relevant species possible. Many of the review papers published around that time promised that these new technologies would revolutionize environmental biology as they had revolutionized medicine and related fields. A few years have passed since these technological advancements have been made, and the use of microarray studies in non-model fish species has been adopted in many laboratories internationally. Has the relatively widespread adoption of this technology really revolutionized the fields of environmental biology, including ecotoxicology, aquaculture and ecology, as promised? Or have these studies merely become a novelty and a potential distraction for scientists addressing environmentally relevant questions? In this review, the promises made in early review papers, in particular about the advances that the use of microarrays would enable, are summarized; these claims are compared to the results of recent studies to determine whether the forecasted changes have materialized. Some applications, as discussed in the paper, have been realized and have led to advances in their field, others are still under development.

Testes and brain gene expression in precocious male and adult maturing Atlantic salmon (Salmo salar)

BACKGROUND

The male Atlantic salmon generally matures in fresh water upon returning after one or several years at sea. Some fast-growing male parr develop an alternative life strategy where they sexually mature before migrating to the oceans. These so called 'precocious' parr or 'sneakers' can successfully fertilise adult female eggs and so perpetuate their line. We have used a custom-built cDNA microarray to investigate gene expression changes occurring in the salmon gonad and brain associated with precocious maturation. The microarray has been populated with genes selected specifically for involvement in sexual maturation (precocious and adult) and in the parr-smolt transformation.

RESULTS

Immature and mature parr collected from a hatchery-reared stock in January were significantly different in weight, length and condition factor. Changes in brain expression were small - never more than 2-fold on the microarray, and down-regulation of genes was much more pronounced than up-regulation. Significantly changing genes included isotocin, vasotocin, cathepsin D, anamorsin and apolipoprotein E. Much greater changes in expression were seen in the testes. Among those genes in the testis with the most significant changes in expression were anti-Mullerian hormone, collagen 1A, and zinc finger protein (Zic1), which were down-regulated in precocity and apolipoproteins E and C-1, lipoprotein lipase and anti-leukoproteinase precursor which were up-regulated in precocity. Expression changes of several genes were confirmed in individual fish by quantitative PCR and several genes (anti-Mullerian hormone, collagen 1A, beta-globin and guanine nucleotide binding protein (G protein) beta polypeptide 2-like 1 (GNB2L1) were also examined in adult maturing testes. Down-regulation of anti-Mullerian hormone was judged to be greater than 160-fold for precocious males and greater than 230-fold for November adult testes in comparison to July testes by this method. For anti-Mullerian hormone and guanine nucleotide binding protein beta polypeptide 2-like 1 expression changes in precocious males mirrored mature adults (November) but for collagen 1A and beta-globin the pattern was more complex.

CONCLUSIONS

Expression changes in the fish brain during the process of precocious sexual maturation were small compared to those in the testes. Microarray analysis suggested down-regulation of housekeeping functions and up-regulation of a small number of specific processes. Transcriptional changes in the testes were much more pronounced with anti-Mullerian hormone playing a major role. Expression profiles for mature parr and maturing adult testes indicate subtle differences in gene expression between these two related groups.

Key considerations for measuring allelic expression on a genomic scale using high-throughput sequencing

Differences in gene expression are thought to be an important source of phenotypic diversity, so dissecting the genetic components of natural variation in gene expression is important for understanding the evolutionary mechanisms that lead to adaptation. Gene expression is a complex trait that, in diploid organisms, results from transcription of both maternal and paternal alleles. Directly measuring allelic expression rather than total gene expression offers greater insight into regulatory variation. The recent emergence of high-throughput sequencing offers an unprecedented opportunity to study allelic transcription at a genomic scale for virtually any species. By sequencing transcript pools derived from heterozygous individuals, estimates of allelic expression can be directly obtained. The statistical power of this approach is influenced by the number of transcripts sequenced and the ability to unambiguously assign individual sequence fragments to specific alleles on the basis of transcribed nucleotide polymorphisms. Here, using mathematical modelling and computer simulations, we determine the minimum sequencing depth required to accurately measure relative allelic expression and detect allelic imbalance via high-throughput sequencing under a variety of conditions. We conclude that, within a species, a minimum of 500-1000 sequencing reads per gene are needed to test for allelic imbalance, and consequently, at least five to 10 millions reads are required for studying a genome expressing 10 000 genes. Finally, using 454 sequencing, we illustrate an application of allelic expression by testing for cis-regulatory divergence between closely related Drosophila species.

Gene-expression signatures of Atlantic salmon's plastic life cycle

How genomic expression differs as a function of life history variation is largely unknown. Atlantic salmon exhibits extreme alternative life histories. We defined the gene-expression signatures of wild-caught salmon at two different life stages by comparing the brain expression profiles of mature sneaker males and immature males, and early migrants and late migrants. In addition to life-stage-specific signatures, we discovered a surprisingly large gene set that was differentially regulated-at similar magnitudes, yet in opposite direction-in both life history transitions. We suggest that this co-variation is not a consequence of many independent cellular and molecular switches in the same direction but rather represents the molecular equivalent of a physiological shift orchestrated by one or very few master regulators.

Major disruption of gene expression in hybrids between young sympatric anadromous and resident populations of brook charr (Salvelinus fontinalis Mitchill)

Genome-wide analyses of the transcriptome have suggested that male-biased genes are the first targets of genomic incompatibilities (g.i.) in inter-specific hybrids. However, those studies have almost invariably focused on Drosophila species that diverged at least 0.9 Ma, and with sterile male hybrids. Here, we use microarrays to analyse patterns of gene expression in very closely related (divergence <12,000 years), sympatric, but ecologically divergent anadromous and resident populations of brook charr (Salvelinus fontinalis) and their F(1) hybrids. Our results show a dramatic breakdown of gene expression patterns in hybrids compared with their parental relatives. Several disrupted genes are related to energetic metabolism, immune response, osmoregulation and protection against oxidative stress, and none has sex-biased functions. Besides, pure individuals show no expression differences at most of the genes disrupted in hybrids, which may suggest the operation of some form of stabilizing selection. Taken together, these results both confirm the idea that perturbations of regulatory networks represent a significant source of g.i. and support the suggestion that developmental pathways can diverge through time without any manifest change in the phenotypic outcome. While the role of other evolutionary forces (e.g. genetic drift) cannot be ruled out, this study suggests that ecological selective processes may provide the initial driving force behind disruption of gene expression in inter-specific hybrids.