DNA methylation reprogramming, TE derepression, and postzygotic isolation of nascent animal species.

The genomic shock hypothesis stipulates that the stress associated with divergent genome admixture can cause transposable element (TE) derepression, which could act as a postzygotic isolation mechanism. TEs affect gene structure, expression patterns, and chromosome organization and may have deleterious consequences when released. For these reasons, they are silenced by heterochromatin formation, which includes DNA methylation. Here, we show that a significant proportion of TEs are differentially methylated between the "dwarf" (limnetic) and the "normal" (benthic) whitefish, two nascent species that diverged some 15,000 generations ago within the Coregonus clupeaformis species complex. Moreover, TEs are overrepresented among loci that were demethylated in hybrids, indicative of their transcriptional derepression. These results are consistent with earlier studies in this system that revealed TE transcriptional derepression causes abnormal embryonic development and death of hybrids. Hence, this supports a role of DNA methylation reprogramming and TE derepression in postzygotic isolation of nascent animal species.

The role of ecotype-environment interactions in intraspecific trophic niche partitioning subsequent to stocking.

Worldwide, stocking of fish represents a valuable tool for conservation and maintenance of species exploited by recreational fishing. Releases of hatchery-reared fish are more and more recognized to have numerous demographic, ecological, and genetic impacts on wild populations. However, consequences on intraspecific trophic relationships have rarely been investigated. In this study, we assessed the impacts of supplementation stocking and resulting introgressive hybridization on the trophic niches occupied by stocked, local, and hybrid lake trout (Salvelinus namaycush) within populations of piscivorous and planktivorous ecotypes stocked from a wild piscivorous source population. We compared trophic niches using stable isotope analysis (δ13 C and δ15 N) and trophic position among the three genetic origins. Putative genetic effects were tested with phenotype-genotype association of "life history" ecological traits (body size, growth rate, condition index, and trophic niche) and genotypes (RADseq SNP markers) using redundant discriminant analysis (RDA). Results showed that sympatry resulting from the stocking of contrasting ecotypes is a risk factor for niche partitioning. Planktivorous populations are more susceptible to niche partitioning, by competitive exclusion of the local fish from a littoral niche to an alternative pelagic/profundal niche. Observed niche partitioning is probably a manifestation of competitive interactions between ecotypes. Our results emphasize that ecotypic variation should be considered for more efficient management and conservation practices and in order to mitigate negative impact of supplementation stocking.

Divergence in physiological factors affecting swimming performance between anadromous and resident populations of brook charr Salvelinus fontinalis.

In this study, an anadromous strain (L) and a freshwater-resident (R) strain of brook charr Salvelinus fontinalis as well as their reciprocal hybrids, were reared in a common environment and submitted to swimming tests combined with salinity challenges. The critical swimming speeds (Ucrit ) of the different crosses were measured in both fresh (FW) and salt water (SW) and the variations in several physiological traits (osmotic, energetic and metabolic capacities) that are predicted to influence swimming performance were documented. Anadromous and resident fish reached the same Ucrit in both FW and SW, with Ucrit being 14% lower in SW compared with FW. The strains, however, seemed to use different underlying strategies: the anadromous strain relied on its streamlined body shape and higher osmoregulatory capacity, while the resident strain had greater citrate synthase (FW) and lactate dehydrogenase (FW, SW) capacity and either greater initial stores or more efficient use of liver (FW, SW) and muscle (FW) glycogen during exercise. Compared with R♀ L♂ hybrids, L♀ R♂ hybrids had a 20% lower swimming speed, which was associated with a 24% smaller cardio-somatic index and higher physiological costs. Thus swimming performance depends on cross direction (i.e. which parental line was used as dam or sire). The study thus suggests that divergent physiological factors between anadromous and resident S. fontinalis may result in similar swimming capacities that are adapted to their respective lifestyles.

On the maintenance of genetic variation and adaptation to environmental change: considerations from population genomics in fishes.

The first goal of this paper was to overview modern approaches to local adaptation, with a focus on the use of population genomics data to detect signals of natural selection in fishes. Several mechanisms are discussed that may enhance the maintenance of genetic variation and evolutionary potential, which have been overlooked and should be considered in future theoretical development and predictive models: the prevalence of soft sweeps, polygenic basis of adaptation, balancing selection and transient polymorphisms, parallel evolution, as well as epigenetic variation. Research on fish population genomics has provided ample evidence for local adaptation at the genome level. Pervasive adaptive evolution, however, seems to almost never involve the fixation of beneficial alleles. Instead, adaptation apparently proceeds most commonly by soft sweeps entailing shifts in frequencies of alleles being shared between differentially adapted populations. One obvious factor contributing to the maintenance of standing genetic variation in the face of selective pressures is that adaptive phenotypic traits are most often highly polygenic, and consequently the response to selection should derive mostly from allelic co-variances among causative loci rather than pronounced allele frequency changes. Balancing selection in its various forms may also play an important role in maintaining adaptive genetic variation and the evolutionary potential of species to cope with environmental change. A large body of literature on fishes also shows that repeated evolution of adaptive phenotypes is a ubiquitous evolutionary phenomenon that seems to occur most often via different genetic solutions, further adding to the potential options of species to cope with a changing environment. Moreover, a paradox is emerging from recent fish studies whereby populations of highly reduced effective population sizes and impoverished genetic diversity can apparently retain their adaptive potential in some circumstances. Although more empirical support is needed, several recent studies suggest that epigenetic variation could account for this apparent paradox. Therefore, epigenetic variation should be fully integrated with considerations pertaining to role of soft sweeps, polygenic and balancing selection, as well as repeated adaptation involving different genetic basis towards improving models predicting the evolutionary potential of species to cope with a changing world.

Making sense of the relationships between Ne, Nb and Nc towards defining conservation thresholds in Atlantic salmon (Salmo salar).

Effective population size over a generation (Ne) or over a reproductive cycle (Nb) and the adult census size (Nc) are important parameters in both conservation and evolutionary biology. Ne provides information regarding the rate of loss of genetic diversity and can be tracked back in time to infer demographic history of populations, whereas Nb may often be more easily quantified than Nc for short-term abundance monitoring. In this study, we propose (1) an empirical context to Waples et al. (2014) who introduced a correction to bias due to overlapping generations, and (2) a mathematical relationship between Ne and Nb for direct application in Atlantic salmon populations in Québec, Canada. To achieve this, we investigate the relationships between Ne, Nb and Nc in 10 Atlantic salmon populations, Canada, for which we genotyped 100 randomly sampled young-of-the year individuals for 5 consecutive years. The results show a positive correlation between Ne, Nb and Nc, suggesting that Nb is an indicative parameter for tracking effective population size and abundance of Atlantic salmon. However, our model allows predicting Nc from Nb values at 27% that can be partly explained by high variance in Nb/Nc both among populations (37%) and among years (19%). This result illustrates the need for thorough calibration of Nb/Nc before using Nb in monitoring programs, as well as a full understanding of the limits of such an approach. Finally, we discuss the importance of these results for the management of wild populations.

Investigating genomic and phenotypic parallelism between piscivorous and planktivorous lake trout (Salvelinus namaycush) ecotypes by means of RADseq and morphometrics analyses.

Repeated adaptive ecological diversification has commonly been reported in fish and has often been associated with trophic niche diversity. The main goal of this study was to investigate the extent of parallelism in the genomic and phenotypic divergence between piscivorous and planktivorous lake trout ecotypes from Laurentian Shield lakes, Canada. This was achieved by documenting the extent of morphological differentiation using geometric morphometrics and linear measurements as well as the pattern of genomic divergence by means of RADseq genotyping (3925 filtered SNPs) in 12 lakes. Our results indicate that the two ecotypes evolved distinct body shape and several linear measurements in parallel. Neutral genetic differentiation was pronounced between all isolated populations (Mean FST  = 0.433), indicating no or very limited migration and pronounced genetic drift. Significant genetic differentiation also suggested partial reproductive isolation between ecotypes in the two lakes where they are found in sympatry. Combining different outlier detection methods, we identified 48 SNPs putatively under divergent selection between ecotypes, among which 10 could be annotated and related to functions such as developmental processes and ionic regulation. Finally, our results indicate that parallel morphological divergence is accompanied by both parallel and nonparallel genomic divergence, which is associated with the use of different trophic niches between ecotypes. The results are also discussed in the context of management and conservation of this highly exploited species throughout northern North America.

Genomewide single nucleotide polymorphism discovery in Atlantic salmon (Salmo salar): validation in wild and farmed American and European populations.

A considerable number of single nucleotide polymorphisms (SNPs) are required to elucidate genotype-phenotype associations and determine the molecular basis of important traits. In this work, we carried out de novo SNP discovery accounting for both genome duplication and genetic variation from American and European salmon populations. A total of 9 736 473 nonredundant SNPs were identified across a set of 20 fish by whole-genome sequencing. After applying six bioinformatic filtering steps, 200 K SNPs were selected to develop an Affymetrix Axiom(®) myDesign Custom Array. This array was used to genotype 480 fish representing wild and farmed salmon from Europe, North America and Chile. A total of 159 099 (79.6%) SNPs were validated as high quality based on clustering properties. A total of 151 509 validated SNPs showed a unique position in the genome. When comparing these SNPs against 238 572 markers currently available in two other Atlantic salmon arrays, only 4.6% of the SNP overlapped with the panel developed in this study. This novel high-density SNP panel will be very useful for the dissection of economically and ecologically relevant traits, enhancing breeding programmes through genomic selection as well as supporting genetic studies in both wild and farmed populations of Atlantic salmon using high-resolution genomewide information.

RAD sequencing reveals within-generation polygenic selection in response to anthropogenic organic and metal contamination in North Atlantic Eels.

Measuring the effects of selection on the genome imposed by human-altered environment is currently a major goal in ecological genomics. Given the polygenic basis of most phenotypic traits, quantitative genetic theory predicts that selection is expected to cause subtle allelic changes among covarying loci rather than pronounced changes at few loci of large effects. The goal of this study was to test for the occurrence of polygenic selection in both North Atlantic eels (European Eel, Anguilla anguilla and American Eel, A. rostrata), using a method that searches for covariation among loci that would discriminate eels from 'control' vs. 'polluted' environments and be associated with specific contaminants acting as putative selective agents. RAD-seq libraries resulted in 23 659 and 14 755 filtered loci for the European and American Eels, respectively. A total of 142 and 141 covarying markers discriminating European and American Eels from 'control' vs. 'polluted' sampling localities were obtained using the Random Forest algorithm. Distance-based redundancy analyses (db-RDAs) were used to assess the relationships between these covarying markers and concentration of 34 contaminants measured for each individual eel. PCB153, 4'4'DDE and selenium were associated with covarying markers for both species, thus pointing to these contaminants as major selective agents in contaminated sites. Gene enrichment analyses suggested that sterol regulation plays an important role in the differential survival of eels in 'polluted' environment. This study illustrates the power of combining methods for detecting signals of polygenic selection and for associating variation of markers with putative selective agents in studies aiming at documenting the dynamics of selection at the genomic level and particularly so in human-altered environments.

Signatures of natural selection between life cycle stages separated by metamorphosis in European eel.

BACKGROUND: Species showing complex life cycles provide excellent opportunities to study the genetic associations between life cycle stages, as selective pressures may differ before and after metamorphosis. The European eel presents a complex life cycle with two metamorphoses, a first metamorphosis from larvae into glass eels (juvenile stage) and a second metamorphosis into silver eels (adult stage). We tested the hypothesis that different genes and gene pathways will be under selection at different life stages when comparing the genetic associations between glass eels and silver eels. RESULTS: We used two sets of markers to test for selection: first, we genotyped individuals using a panel of 80 coding-gene single nucleotide polymorphisms (SNPs) developed in American eel; second, we investigated selection at the genome level using a total of 153,423 RAD-sequencing generated SNPs widely distributed across the genome. Using the RAD approach, outlier tests identified a total of 2413 (1.57%) potentially selected SNPs. Functional annotation analysis identified signal transduction pathways as the most over-represented group of genes, including MAPK/Erk signalling, calcium signalling and GnRH (gonadotropin-releasing hormone) signalling. Many of the over-represented pathways were related to growth, while others could result from the different conditions that eels inhabit during their life cycle. CONCLUSIONS: The observation of different genes and gene pathways under selection when comparing glass eels vs. silver eels supports the adaptive decoupling hypothesis for the benefits of metamorphosis. Partitioning the life cycle into discrete morphological phases may be overall beneficial since it allows the different life stages to respond independently to their unique selection pressures. This might translate into a more effective use of food and niche resources and/or performance of phase-specific tasks (e.g. feeding in the case of glass eels, migrating and reproducing in the case of silver eels).

How does salinity influence habitat selection and growth in juvenile American eels Anguilla rostrata?

The influence of salinity on habitat selection and growth in juvenile American eels Anguilla rostrata captured in four rivers across eastern Canada was assessed in controlled experiments in 2011 and 2012. Glass eels were first categorized according to their salinity preferences towards fresh (FW), salt (SW) or brackish water (BW) and the growth rate of each group of elvers was subsequently monitored in controlled FW and BW environments for 7 months. Most glass eels (78-89%) did not make a choice, i.e. they remained in BW. Salinity preferences were not influenced by body condition, although a possible role of pigmentation could not be ruled out. Glass eels that did make a choice displayed a similar preference for FW (60-75%) regardless of their geographic origin but glass eels from the St Lawrence Estuary displayed a significantly higher locomotor activity than those from other regions. Neither the salinity preferences showed by glass eels in the first experiment nor the rearing salinities appeared to have much influence on growth during the experiments. Elvers from Nova Scotia, however, reached a significantly higher mass than those from the St Lawrence Estuary thus supporting the hypothesis of genetically (or epigenetically) based differences for growth between A. rostrata from different origins. These results provide important ecological knowledge for the sustained exploitation and conservation of this threatened species.

Speciation and demographic history of Atlantic eels (Anguilla anguilla and A. rostrata) revealed by mitogenome sequencing.

Processes leading to speciation in oceanic environments without obvious physical barriers remain poorly known. European and American eel (Anguilla anguilla and A. rostrata) spawn in partial sympatry in the Sargasso Sea. Larvae are advected by the Gulf Stream and other currents towards the European/North African and North American coasts, respectively. We analyzed 104 mitogenomes from the two species along with mitogenomes of other Anguilla and outgroup species. We estimated divergence time between the two species to identify major events involved in speciation. We also considered two previously stated hypotheses: one where the ancestral species was present in only one continent but was advected across the Atlantic by ocean current changes and another where population declines during Pleistocene glaciations led to increasing vicariance, facilitating speciation. Divergence time was estimated to ∼3.38 Mya, coinciding with the closure of the Panama Gateway that led to reinforcement of the Gulf Stream. This could have advected larvae towards European/North African coasts, in which case American eel would be expected to be the ancestral species. This scenario could, however, not be unequivocally confirmed by analyses of dN/dS, nucleotide diversity and effective population size estimates. Extended bayesian skyline plots showed fluctuations of effective population sizes and declines during glaciations, and thus also lending support to the importance of vicariance during speciation. There was evidence for positive selection at the ATP6 and possibly ND5 genes, indicating a role in speciation. The findings suggest an important role of ocean current changes in speciation of marine organisms.

Microbiome investigation in the ecological speciation context of lake whitefish (Coregonus clupeaformis) using next-generation sequencing.

Few studies have applied NGS methods to investigate the microbiome of vertebrates in their natural environment and in freshwater fishes in particularly. Here, we used pyrosequencing of the 16S gene rRNA to (i) test for differences in kidney bacterial communities (i.e. microbiota) of dwarf and normal whitefish found as sympatric pairs, (ii) test the hypothesis of higher bacterial diversity in normal compared with dwarf whitefish and (iii) test for the occurrence of parallelism with the presence and composition of bacterial communities across species pairs inhabiting different lakes. The kidney microbiota of 253 dwarf and normal whitefish from five lakes was analysed combining a double-nested PCR approach with 454 pyrosequencing. Bacteria were detected in 52.6% of the analysed whitefish. There was no overall significant difference among lakes and forms, although the lake × form interaction was found significant. We identified 579 bacterial genera, which is substantially more than previous descriptions using less sensitive techniques of fish bacterial diversity in kidney, pathogenic or not. Ten of these genera contained eighteen pathogenic species. Differences in bacteria composition between whitefish forms were not parallel among lakes. In accordance with the higher diversity of prey types, normal whitefish kidney tissue consistently had a more diverse bacterial community and this pattern was parallel among lakes. These results add to building evidence from previous studies on this system that the adaptive divergence of dwarf, and normal whitefish has been driven by both parallel and nonparallel ecological conditions across lakes.

Genome-wide single-generation signatures of local selection in the panmictic European eel.

Next-generation sequencing and the collection of genome-wide data allow identifying adaptive variation and footprints of directional selection. Using a large SNP data set from 259 RAD-sequenced European eel individuals (glass eels) from eight locations between 34 and 64(o) N, we examined the patterns of genome-wide genetic diversity across locations. We tested for local selection by searching for increased population differentiation using F(ST) -based outlier tests and by testing for significant associations between allele frequencies and environmental variables. The overall low genetic differentiation found (F(ST) = 0.0007) indicates that most of the genome is homogenized by gene flow, providing further evidence for genomic panmixia in the European eel. The lack of genetic substructuring was consistent at both nuclear and mitochondrial SNPs. Using an extensive number of diagnostic SNPs, results showed a low occurrence of hybrids between European and American eel, mainly limited to Iceland (5.9%), although individuals with signatures of introgression several generations back in time were found in mainland Europe. Despite panmixia, a small set of SNPs showed high genetic differentiation consistent with single-generation signatures of spatially varying selection acting on glass eels. After screening 50 354 SNPs, a total of 754 potentially locally selected SNPs were identified. Candidate genes for local selection constituted a wide array of functions, including calcium signalling, neuroactive ligand-receptor interaction and circadian rhythm. Remarkably, one of the candidate genes identified is PERIOD, possibly related to differences in local photoperiod associated with the >30° difference in latitude between locations. Genes under selection were spread across the genome, and there were no large regions of increased differentiation as expected when selection occurs within just a single generation due to panmixia. This supports the conclusion that most of the genome is homogenized by gene flow that removes any effects of diversifying selection from each new generation.

Assessing patterns of hybridization between North Atlantic eels using diagnostic single-nucleotide polymorphisms.

The two North Atlantic eel species, the European eel (Anguilla anguilla) and the American eel (Anguilla rostrata), spawn in partial sympatry in the Sargasso Sea, providing ample opportunity to interbreed. In this study, we used a RAD (Restriction site Associated DNA) sequencing approach to identify species-specific diagnostic single-nucleotide polymorphisms (SNPs) and design a low-density array that combined with screening of a diagnostic mitochondrial DNA marker. Eels from Iceland (N=159) and from the neighboring Faroe Islands (N=29) were genotyped, along with 94 larvae (49 European and 45 American eel) collected in the Sargasso Sea. Our SNP survey showed that the majority of Icelandic eels are pure European eels but there is also an important contribution of individuals of admixed ancestry (10.7%). Although most of the hybrids were identified as F1 hybrids from European eel female × American eel male crosses, backcrosses were also detected, including a first-generation backcross (F1 hybrid × pure European eel) and three individuals identified as second-generation backcrosses originating from American eel × F1 hybrid backcrosses interbreeding with pure European eels. In comparison, no hybrids were observed in the Faroe Islands, the closest bodies of land to Iceland. It is possible that hybrids show an intermediate migratory behaviour between the two parental species that ultimately brings hybrid larvae to the shores of Iceland, situated roughly halfway between the Sargasso Sea and Europe. Only two hybrids were observed among Sargasso Sea larvae, both backcrosses, but no F1 hybrids, that points to temporal variation in the occurrence of hybridization.

Linking transcriptomic and genomic variation to growth in brook charr hybrids (Salvelinus fontinalis, Mitchill).

Hybridization can lead to phenotypic differences arising from changes in gene expression patterns or new allele combinations. Variation in gene expression is thought to be controlled by differences in transcription regulation of parental alleles, either through cis- or trans-regulatory elements. A previous study among brook charr hybrids from different populations (Rupert, Laval, and domestic) showing distinct length at age during early life stages also revealed different patterns in transcription regulation inheritance of transcript abundance. In the present study, transcript abundance using RNA-sequencing and quantitative real-time PCR, single-nucleotide polymorphism (SNP) genotypes and allelic imbalance were assessed in order to understand the molecular mechanisms underlying the observed transcriptomic and differences in length at age among domestic × Rupert hybrids and Laval × domestic hybrids. We found 198 differentially expressed genes between the two hybrid crosses, and allelic imbalance could be analyzed for 69 of them. Among these 69 genes, 36 genes exhibited cis-acting regulatory effects in both of the two crosses, thus confirming the prevalent role of cis-acting regulatory elements in the regulation of differentially expressed genes among intraspecific hybrids. In addition, we detected a significant association between SNP genotypes of three genes and length at age. Our study is thus one of the few that have highlighted some of the molecular mechanisms potentially involved in the differential phenotypic expression in intraspecific hybrids for nonmodel species.

The influence of parental effects on transcriptomic landscape during early development in brook charr (Salvelinus fontinalis, Mitchill).

Parental effects represent an important source of variation in offspring phenotypes. Depending on the specific mechanisms involved, parental effects may be caused to different degrees by either the maternal or the paternal parent, and these effects may in turn act at different stages of development. To detect parental effects acting on gene transcription regulation and length phenotype during ontogeny, the transcriptomic profiles of two reciprocal hybrids from Laval × Rupert and Laval × Domestic populations of brook charr were compared at hatching, yolk sac resorption and 15 weeks after exogenous feeding. Using a salmonid cDNA microarray, our results show that parental effects modulated gene expression among reciprocal hybrids only at the yolk sac resorption stage. In addition, Laval × Domestic and Laval × Rupert reciprocal hybrids differed in the magnitude of theses parental effects, with 199 and 630 differentially expressed transcripts, respectively. This corresponds to a maximum of 18.5% of the analyzed transcripts. These transcripts are functionally related to cell cycle, nucleic acid metabolism and intracellular protein traffic, which is consistent with observed differences associated with embryonic development and growth differences in other fish species. Our results thus illustrate how parental effects on patterns of gene transcription seem dependent on the genetic architecture of the parents. In addition, in absence of transcriptional differences, non-transcript deposits in the yolk sac could contribute to the observed length differences among the reciprocal hybrids before yolk sac resorption.

Phylogenetic status of brown trout Salmo trutta populations in five rivers from the southern Caspian Sea and two inland lake basins, Iran: a morphogenetic approach.

Interrelationships, origin and phylogenetic affinities of brown trout Salmo trutta populations from the southern Caspian Sea basin, Orumieh and Namak Lake basins in Iran were analysed from complete mtDNA control region sequences, 12 microsatellite loci and morphological characters. Among 129 specimens from six populations, seven haplotypes were observed. Based on mtDNA haplotype data, the Orumieh and southern Caspian populations did not differ significantly, but the Namak basin-Karaj population presented a unique haplotype closely related to the haplotypes of the other populations (0·1% Kimura two-parameter, K2P divergence). All Iranian haplotypes clustered as a distinct group within the Danube phylogenetic grouping, with an average K2P distance of 0·41% relative to other Danubian haplotypes. The Karaj haplotype in the Namak basin was related to a haplotype (Da26) formerly identified in the Tigris basin in Turkey, to a Salmo trutta oxianus haplotype from the Aral Sea basin, and to haplotype Da1a with two mutational steps, as well as to other Iranian haplotypes with one to two mutational steps, which may indicate a centre of origin in the Caspian basin. In contrast to results of the mtDNA analysis, more pronounced differentiation was observed among the populations studied in the morphological and microsatellite DNA data, except for the two populations from the Orumieh basin, which were similar, possibly due to anthropogenic causes.

Genetic differentiation between two sympatric morphs of the blind Iran cave barb Iranocypris typhlops.

The phylogenetic relationship between two sympatric morphotypes of the Iran cave barb Iranocypris typhlops, and Garra rufa, was investigated by sequencing the cytochrome c oxidase I (coI) region (788 bp) providing the first molecular evidence of their phylogeny. Consistent with their morphological differences, the mean genetic distance between the two forms of I. typhlops was significantly higher than generally reported for intraspecific divergence in freshwater fishes. They were phylogenetically closer to G. rufa than to any other species.

Oxidative phosphorylation gene transcription in whitefish species pairs reveals patterns of parallel and nonparallel physiological divergence.

Across multiple lakes in North America, lake whitefish (Coregonus clupeaformis) have independently evolved 'dwarf' and 'normal' sympatric species pairs that exhibit pronounced phenotypic and genetic divergence. In particular, traits associated with metabolism have been shown to be highly differentiated between whitefish species. Here, we examine the transcription of genes associated with the five mitochondrial and nuclear genome-encoded oxidative phosphorylation (OXPHOS) complexes, the primary physiological mechanism responsible for the production of ATP, in whitefish species pairs from Cliff Lake and Webster Lake in Maine, USA. We observed OXPHOS gene transcription divergence between dwarf and normal whitefish in each of the two lakes, with the former exhibiting transcription upregulation for genes associated with each of the OXPHOS complexes. We also observed a significant influence of lake on transcription levels for some of the genes, indicating that inter-lake ecological or genetic differences are contributing to variation in OXPHOS gene transcription levels. Together, our results support the hypothesis that metabolic divergence is a critical adaptation involved in whitefish speciation and implicate OXPHOS gene upregulation as a factor involved in meeting the enhanced energetic demands of dwarf whitefish. Further studies are now needed to evaluate the contribution of genetically vs. plasticity driven variation in transcription associated with this critical physiological pathway.

Plasticity and heritability of morphological variation within and between parapatric stickleback demes.

The threespine stickleback (Gasterosteus aculeatus) has emerged as an important model organism in evolutionary ecology, largely due to the repeated, parallel evolution of divergent morphotypes found in populations having colonized freshwater habitats. However, morphological divergence following colonization is not a universal phenomenon. We explore this in a large-scale estuarine ecosystem inhabited by two parapatric stickleback demes, each physiologically adapted to divergent osmoregulatory environments (fresh vs. saline waters). Using geometric morphometric analyses of wild-caught individuals, we detected significant differences between demes, in addition to sexual dimorphism, in body shape. However, rearing full-sib families from each deme under controlled, reciprocal salinity conditions revealed no differences between genotypes and highly significant environmental effects. It is also noteworthy that fish from both demes were fully plated, whether found in the wild or reared under reciprocal salinity conditions. Although we found significant heritability for body shape, we also noted significant direct environmental effects for many latent shape variables. Moreover, we found little evidence for diversifying selection acting on body size and shape (Q(ST) ). Nevertheless, uniform compressive variation did exceed neutral expectations, yet despite evidence of both allometry and genetic correlation with body length, we detected no correlated signatures of selection. Taken together, these results suggest that much of the morphological divergence observed in this system is the result of plastic responses to environmental variation rather than adaptive differentiation.