BackgroundThere are large differences between the sexes at the genetic level; these differences include heterogametic sex chromosomes and/or differences in expression of genes between the sexes. In rainbow trout (Oncorhynchus mykiss) qRT-PCR studies have found significant differences in expression of several candidate sex determining genes. However, these genes represent a very small fraction of the genome and research in other species suggests there are large portions of the transcriptome that are differentially expressed between the sexes. These differences are especially noticeable once gonad differentiation and maturation has occurred, but less is known at earlier stages of development. Here we use data from a microarray and qRT-PCR to identify genes differentially expressed between the sexes at three time points in pre-hatch embryos, prior to the known timing of sexual differentiation in this species.ResultsThe microarray study revealed 883 differentially expressed features between the sexes with roughly equal numbers of male and female upregulated features across time points. Most of the differentially expressed genes on the microarray were not related to sex function, suggesting large scale differences in gene expression between the sexes are present early in development. Candidate gene analysis revealed sox9, DMRT1, Nr5a1 and wt1 were upregulated in males at some time points and foxl2, ovol1, fst and cyp19a1a were upregulated in females at some time points.ConclusionThis is the first study to identify sexual dimorphism in expression of the genome during embryogenesis in any fish and demonstrates that transcriptional differences are present before the completion of gonadogenesis.

In rainbow trout (Oncorhynchus mykiss) and other fishes, embryonic development rate is an ecologically and evolutionarily important trait that is closely associated with survival and physiological performance later in life. To identify genes differentially regulated in fast and slow-developing embryos of rainbow trout, we examined gene expression across developmental time points in rainbow trout embryos possessing alleles linked to a major quantitative trait loci (QTL) for fast versus slow embryonic development rate. Whole genome expression microarray analyses were conducted using embryos from a fourth generation backcross family, whereby each backcross generation involved the introgression of the fast-developing alleles for a major development rate QTL into a slow-developing clonal line of rainbow trout. Embryos were collected at 15, 19, and 28 days post-fertilization; sex and QTL genotype were determined using molecular markers, and cDNA from 48 embryos were used for microarray analysis. A total of 183 features were identified with significant differences between embryonic development rate genotypes. Genes associated with cell cycle growth, muscle contraction and protein synthesis were expressed significantly higher in embryos with the fast-developing allele (Clearwater) than those with the slow-developing allele (Oregon State University), which may associate with fast growth and early body mass construction in embryo development. Across time points, individuals with the fast-developing QTL allele appeared to have earlier onset of these developmental processes when compared to individuals with the slow development alleles, even as early as 15 days post-fertilization. Differentially expressed candidate genes chosen for linkage mapping were localized primarily to regions outside of the major embryonic development rate QTL, with the exception of a single gene (very low-density lipoprotein receptor precursor).

Development rate has important implications for individual fitness and physiology. In salmonid fishes, development rate correlates with many traits later in life, including life-history diversity, growth, and age and size at sexual maturation. In rainbow trout (Oncorhynchus mykiss), a quantitative trait locus for embryonic development rate has been detected on chromosome 5 across populations. However, few candidate genes have been identified within this region. In this study, we use gene mapping, gene expression, and quantitative genetic methods to further identify the genetic basis of embryonic developmental rate in O. mykiss Among the genes located in the region of the major development rate quantitative trait locus (GHR1, Clock1a, Myd118-1, and their paralogs), all were expressed early in embryonic development (fertilization through hatch), but none were differentially expressed between individuals with the fast- or slow-developing alleles for a major embryonic development rate quantitative trait locus. In a follow-up study of migratory and resident rainbow trout from natural populations in Alaska, we found significant additive variation in development rate and, moreover, found associations between development rate and allelic variation in all 3 candidate genes within the quantitative trait locus for embryonic development. The mapping of these genes to this region and associations in multiple populations provide positional candidates for further study of their roles in growth, development, and life-history diversity in this model salmonid.

Fitness related traits often show spatial variation across populations of widely distributed species. Comparisons of genetic variation among populations in putatively neutral DNA markers and in phenotypic traits susceptible to selection (QST FST analysis) can be used to determine to what degree differentiation among populations can be attributed to selection or genetic drift. Traditionally, QST FST analyses require a large number of populations to achieve sufficient statistical power; however, new methods have been developed that allow QST FST comparisons to be conducted on as few as two populations if their pedigrees are informative. This study compared genetic and morphological divergence in three strains of brook trout Salvelinus fontinalis that were historically or currently used for stocking in the Lake Superior Basin. Herein we examined if morphological divergence among populations showed temporal variation, and if divergence could be attributed to selection or was indistinguishable from genetic drift. Multivariate QST FST analysis showed evidence for divergent selection between populations. Univariate analyses suggests that the pattern observed in the multivariate analyses was largely driven by divergent selection for length and weight, and moreover by divergence between the Assinica strain and each of the Iron River and Siskiwit strains rather than divergent selection between each population pair. While it could not be determined if divergence was due to natural selection or inadvertent artificial selection in hatcheries, selected differences were consistent with patterns of domestication commonly found in salmonids.