Search

Captive rearing in salmon hatcheries can have considerable impacts on both fish phenotype and fitness within a single generation, even in the absence of genetic change. Evidence for hatchery-induced changes in DNA methylation is becoming abundant, though questions remain on the sex-specificity of these effects, their persistence until spawning, and potential for transmission to future generations. Here we performed whole genome methylation sequencing of fin tissue for 16 hatchery and 16 wild Atlantic salmon (Salmo salar) returning to spawn in the Rimouski River, Québec. We identified two cohorts of hatchery-reared salmon through methylation analysis, one of which was epigenetically similar to wild fish, suggesting that supplementation efforts may be able to minimize the epigenetic effects of hatchery rearing. We found considerable sex-specific effects of hatchery rearing, with few genomic regions being affected in both males and females. We also analysed the methylome of 32 F1 offspring from four groups (pure wild, pure hatchery origin, and reciprocal hybrids). We found that few epigenetic changes due to parental hatchery rearing persisted in the F1 offspring though the patterns of inheritance appear to be complex, involving nonadditive effects. Our results suggest that the epigenetic effects of hatchery rearing can be minimal in F0. There may also be minimal epigenetic inheritance and rapid loss of epigenetic changes associated with hatchery rearing. However, due to sex-specificity and nonadditive patterns of inheritance, methylation changes due to captive rearing are rather complex and the field would benefit from further research on minimizing the epigenetic effects of captive rearing in conservation efforts.

Human activities induce direct or indirect selection pressure on natural population and may ultimately affect population's integrity. While numerous conservation programs aimed to minimize human‐induced genomic variation, human‐induced environmental variation may generate epigenomic variation potentially affecting fitness through phenotypic modifications. Major questions remain pertaining to how much epigenomic variation arises from environmental heterogeneity, whether this variation can persist throughout life, and whether it can be transmitted across generations. We performed whole genome bisulfite sequencing (WGBS) on the sperm of genetically indistinguishable hatchery and wild‐born migrating adults of Coho salmon (Oncorhynchus kisutch) from two geographically distant rivers at different epigenome scales. Our results showed that coupling WGBS with fine‐scale analyses (local and chromosomal) allows the detection of parallel early‐life hatchery‐induced epimarks that differentiate wild from hatchery‐reared salmon. Four chromosomes and 183 differentially methylated regions (DMRs) displayed a significant signal of methylation differentiation between hatchery and wild‐born Coho salmon. Moreover, those early‐life epimarks persisted in germ line cells despite about 1.5 year spent in the ocean following release from hatchery, opening the possibility for transgenerational inheritance. Our results strengthen the hypothesis that epigenomic modifications environmentally induced during early‐life development persist in germ cells of adults until reproduction, which could potentially impact their fitness.

Here, we present an annotated, chromosome-anchored, genome assembly for Lake Trout (Salvelinus namaycush) - a highly diverse salmonid species of notable conservation concern and an excellent model for research on adaptation and speciation. We leveraged Pacific Biosciences long-read sequencing, paired-end Illumina sequencing, proximity ligation (Hi-C) sequencing, and a previously published linkage map to produce a highly contiguous assembly composed of 7378 contigs (contig N50 = 1.8 Mb) assigned to 4120 scaffolds (scaffold N50 = 44.975 Mb). Long read sequencing data were generated using DNA from a female double haploid individual. 84.7% of the genome was assigned to 42 chromosome-sized scaffolds and 93.2% of Benchmarking Universal Single Copy Orthologues were recovered, putting this assembly on par with the best currently available salmonid genomes. Estimates of genome size based on k-mer frequency analysis were highly similar to the total size of the finished genome, suggesting that the entirety of the genome was recovered. A mitochondrial genome assembly was also produced. Self-versus-self synteny analysis allowed us to identify homeologs resulting from the salmonid specific autotetraploid event (Ss4R) as well as regions exhibiting delayed rediploidization. Alignment with three other salmonid genomes and the Northern Pike (Esox lucius) genome also allowed us to identify homologous chromosomes in related taxa. We also generated multiple resources useful for future genomic research on Lake Trout, including a repeat library and a sex-averaged recombination map. A novel RNA sequencing data set for liver tissue was also generated in order to produce a publicly available set of annotations for 49,668 genes and pseudogenes. Potential applications of these resources to population genetics and the conservation of native populations are discussed.

Distinguishing neutral and adaptive genetic variation is one of the main challenges in investigating processes shaping population structure in the wild, and landscape genomics can help identify signatures of adaptation to contrasting environments. Arctic Char (Salvelinus alpinus) is an anadromous salmonid and the most harvested fish species by Inuit people, including in Nunavik (Québec, Canada), one of the most recently deglaciated regions in the world. Unlike many other anadromous salmonids, Arctic Char occupy coastal habitats near their natal rivers during their short marine phase restricted to the summer ice‐free period. Our main objective was to document putatively neutral and adaptive genomic variation in anadromous Arctic Char populations from Nunavik and bordering regions to inform local fisheries management. We used genotyping by sequencing (GBS) to genotype 18,112 filtered single nucleotide polymorphisms (SNP) in 650 individuals from 23 sampling locations along >2000 km of coastline. Our results reveal a hierarchical genetic structure, whereby neighboring hydrographic systems harbor distinct populations grouped by major oceanographic basins: Hudson Bay, Hudson Strait, Ungava Bay, and Labrador Sea. We found genetic diversity and differentiation to be consistent both with the expected postglacial recolonization history and with patterns of isolation‐by‐distance reflecting contemporary gene flow. Results from three gene–environment association methods supported the hypothesis of local adaptation to both freshwater and marine environments (strongest associations with sea surface and air temperatures during summer and salinity). Our results support a fisheries management strategy at a regional scale, and other implications for hatchery projects and adaptation to climate change are discussed.

Protecting freshwater biodiversity is considered an ultimate challenge but depends on reliable surveys of species distribution and abundance which eDNA metabarcoding (environmental DNA metabarcoding) may offer. To do so, a better understanding of the sources of temporal variation among species eDNA abundance and of data transformation in eDNA metabarcoding studies is needed. Here, we show that transformation based on relative abundance is critical to suitable analyses of eDNA metabarcoding data and that Hellinger transformation performed slightly better than other methods. Furthermore, we show that site localities significantly explain eDNA metabarcoding variation, while no variation is explained by time of sampling. This indicates that species communities vary more spatially than temporally within a dendritic system composed of small rivers. We then further documented the community structure in the St. Charles River (Québec City, Canada) and six of its tributaries. This revealed the existence of eight species communities explaining 82.1% of eDNA read variation within this river network. Moreover, variation in environmental variables among sites explained 53.0% of eDNA reads, while sampling events and temporal environmental variation explained no eDNA metabarcoding variation. Altogether, this supports the claim that eDNA metabarcoding is a powerful tool to document and monitor fish communities in watersheds composed of small river dendritic systems.

Large rivers and their estuaries are structurally complex and comprise a diversity of habitats supporting a rich biodiversity. As a result, identifying and monitoring fish communities using traditional methods in such systems may often be logistically challenging. Using the mitochondrial DNA 12S MiFish primers, we performed an eDNA metabarcoding analysis to assess the effect of spatial and environmental factors on the variation of the fish community structure along most of the St. Lawrence River/Estuary/Gulf (Québec Canada), a transect spanning 1300 km across a diversity of habitats from a fluviatile non‐tidal section to a marine environment. A total of 129 species were identified including freshwater and marine species. For the freshwater sectors, eDNA identified 80 species compared with the 85 species previously reported based on conventional sampling. eDNA also revealed similar species diversity and communities in the fluviatile section of the St. Lawrence River. Furthermore, our study improved current knowledge about the brackish and marine sections by describing community transition between freshwater and marine fish communities in association with a drastic shift in environmental conditions observed between the end of the fluvial estuary and the beginning of the middle (brackish) estuary. Altogether, this study exemplifies how eDNA metabarcoding is a powerful tool to document fish community shifts in large temperate lotic ecosystems.

Understanding the processes shaping population structure and reproductive isolation of marine organisms can improve their management and conservation. Using genomic markers combined with estimation of individual ancestries, assignment tests, spatial ecology, and demographic modeling, we (i) characterized the contemporary population structure, (ii) assessed the influence of space, fishing depth, and sampling years on contemporary distribution, and (iii) reconstructed the speciation history of two cryptic redfish species, Sebastes mentella and S. fasciatus. We genotyped 860 individuals in the Northwest Atlantic Ocean using 24,603 filtered single nucleotide polymorphisms (SNPs). Our results confirmed the clear genetic distinctiveness of the two species and identified three ecotypes within S. mentella and five populations in S. fasciatus. Multivariate analyses highlighted the influence of spatial distribution and depth on the overall genomic variation, while demographic modeling revealed that secondary contact models best explained inter‐ and intragenomic divergence. These species, ecotypes, and populations can be considered as a rare and wide continuum of genomic divergence in the marine environment. This acquired knowledge pertaining to the evolutionary processes driving population divergence and reproductive isolation will help optimizing the assessment of demographic units and possibly to refine fishery management units.

Biodiversity assessment is an important part of conservation management that ideally can be accomplished with noninvasive methods without influencing the structure and functioning of ecosystems. Environmental DNA (eDNA) metabarcoding has provided a promising tool to enable fast and comprehensive monitoring of entire ecosystems, but widespread adoption of this technique requires performance evaluations that compare it with conventional surveys. We compared eDNA metabarcoding and trawling data to evaluate their efficiency to characterize demersal fish communities in the Estuary and Gulf of Saint‐Lawrence, Canada. Seawater and bottom trawling samples were collected in parallel at 84 stations. For a subset of 30 of these stations, water was also collected at three different depths (15, 50, and 250 m) across the water column. An eDNA metabarcoding assay based on the 12S mitochondrial gene using the MiFish‐U primers was applied to detect fish eDNA. We detected a total of 88 fish species with both methods combined, with 72 species being detected by eDNA, 64 species detected by trawl, and 47 species (53%) overlapped between both methods. eDNA was more efficient for quantifying species richness, mainly because it detected species known to be less vulnerable to trawling gear. Our results indicated that the relative abundance estimated by eDNA and trawl is significantly correlated for species detected by both methods, while the relationship was also influenced by environmental variables (temperature, depth, salinity, and oxygen). Integrating eDNA metabarcoding to bottom trawling surveys could provide additional information on vertical fish distribution in the water column. Environmental DNA metabarcoding thus appears to be a reliable and complementary approach to trawling surveys for documenting fish biodiversity, including for obtaining relative quantitative estimates in the marine environment.

Accurate data characterizing species distribution and abundance are critical for conservation and management of aquatic resources. Inventory methods, such as gillnet surveys, are widely used to estimate distribution and abundance of fish. However, gillnet surveys can be costly in terms of material and human resources, may cause unwanted mortality in the fish communities being studied, and is subject to size and species selection bias. Detecting allochthonous DNA released by species in their environment (i.e., environmental DNA, hereafter eDNA) could be used as a noninvasive and less costly alternative. In this study, we directly compare eDNA metabarcoding and gillnets for monitoring freshwater fish communities in terms of species richness and relative species abundance. Metabarcoding was performed with the 12S Mifish primers. We also used species‐specific quantitative PCR (qPCR) for the most abundant species, the walleye (Sander vitreus), to compare estimated relative abundance with metabarcoding and gillnet captures. Water sample collection, prior to gillnet assessment, was performed on 17 sites in the hydroelectric impoundment of the Rupert River (James Bay, Canada), comparing two water filtration methods. After controlling for amplification biases and repeatability, we show that fish communities’ complexity is better represented using eDNA metabarcoding than previously recorded gillnet data and that metabarcoding read count correlates with qPCR (r = 0.78, p

Large‐scale biomonitoring of Arctic coastal marine communities is essential to track temporal changes in ecosystems. Despite the potential of environmental DNA (eDNA) as an innovative coastal biomonitoring tool, important questions remain pertaining to its temporal and spatial variation and how this may affect the evaluation of ecosystem changes over time in hydrodynamic ecosystems. In this study, we used eDNA metabarcoding of coastal water samples in two Canadian Arctic ports to evaluate the potential of eDNA to detect temporal transition in marine coastal communities. We sequenced eDNA from approximately 20 surface water samples collected each month (N ≈ 150 samples) covering the transition period between summer and late fall using four different universal primer pairs (two pairs of COI mitochondrial genes and two pairs of 18S rRNA genes). Our results from both primer pairs highlighted a significant transition from the summer to the fall marine community. We also observed a putative link between eDNA peaks of read abundance and timing for different life stages (e.g., spawning and larvae) of several species with the most abundant sequence reads. As such, our results show that temporal variation must be considered in ensuring comprehensive coastal biomonitoring with eDNA. Although much remains to be investigated about the ecology of eDNA, our results contribute to fundamental knowledge on the origin of eDNA and highlight the importance of considering temporal variation in developing guidance for coastal biomonitoring with this approach.

Conservation of freshwater biodiversity requires being able to track the presence and abundance of entire fish communities. However, studying fish community composition within rivers remains a technical challenge because of high spatial and temporal physico‐chemical variability, anthropic activities and connections with other river catchments, which may all contribute to important variations in local ecology and communities. Here, we used environmental DNA metabarcoding to document spatial variation in fish communities at a small geographic scale in a large river system. The study was conducted in the Contrecoeur sector (5.5 km long and approximately 1–1.5 km wide) of the St. Lawrence River (Québec, Canada), where two water masses with different physico‐chemical properties, known as "brown waters" and "green waters," flow in parallel with limited admixing. Water samples were collected during two consecutive days at 53 stations located in both water masses. Using universal PCR MiFish 12S primers, Illumina MiSeq sequencing, and the Barque (www.github.com/enormandeau/barque) eDNA analysis software developed by our group, a total of 67 fish species were detected. PERMANOVA and redundancy analyses (RDA) performed on relative read abundance revealed that each water mass comprised distinct communities that depended on turbidity, depth, and to a lesser extent on the upstream versus downstream position along the study area. eDNA metabarcoding results were compared with those of traditional surveys conducted previously in the sector and up to 40 km upstream of it. As previously reported, higher species diversity was detected by eDNA and with substantially lower sampling effort. Our results represent one of the few studies documenting the potential of eDNA metabarcoding to investigate small‐scale variation in 2D spatial patterns of distribution of whole fish communities associated with habitat characteristics within a lotic system.

Copy number variants (CNVs) are a major component of genotypic and phenotypic variation in genomes. To date, our knowledge of genotypic variation and evolution has largely been acquired by means of single nucleotide polymorphism (SNPs) analyses. Until recently, the adaptive role of structural variants (SVs) and particularly that of CNVs has been overlooked in wild populations, partly due to their challenging identification. Here, we document the usefulness of Rapture, a derived reduced-representation shotgun sequencing approach, to detect and investigate copy number variants (CNVs) alongside SNPs in American lobster (Homarus americanus) populations. We conducted a comparative study to examine the potential role of SNPs and CNVs in local adaptation by sequencing 1,141 lobsters from 21 sampling sites within the southern Gulf of St. Lawrence, which experiences the highest yearly thermal variance of the Canadian marine coastal waters. Our results demonstrated that CNVs accounts for higher genetic differentiation than SNP markers. Contrary to SNPs, for which no significant genetic-environment association was found, 48 CNV candidates were significantly associated with the annual variance of sea surface temperature, leading to the genetic clustering of sampling locations despite their geographic separation. Altogether, we provide a strong empirical case that CNVs putatively contribute to local adaptation in marine species and unveil stronger spatial signal of population structure than SNPs. Our study provides the means to study CNVs in non-model species and highlights the importance of considering structural variants alongside SNPs to enhance our understanding of ecological and evolutionary processes shaping adaptive population structure.

Investigating the relative importance of neutral versus selective processes governing the accumulation of genetic variants is a key goal in both evolutionary and conservation biology. This is particularly true in the context of small populations, where genetic drift can counteract the effect of selection. Using Brook Charr (Salvelinus fontinalis) from Québec, Canada, as a case study, we investigated the importance of demographic versus selective processes governing the accumulation of both adaptive and maladaptive mutations in closed versus open and connected populations to assess gene flow effect. This was achieved by using 14,779 high-quality filtered SNPs genotyped among 1,416 fish representing 50 populations from three life history types: lacustrine (closed populations), riverine and anadromous (connected populations). Using the PROVEAN algorithm, we observed a considerable accumulation of putative deleterious mutations across populations. The absence of correlation between the occurrence of putatively beneficial or deleterious mutations and local recombination rate supports the hypothesis that genetic drift might be the main driver of the accumulation of such variants. However, despite a lower genetic diversity observed in lacustrine than in riverine or anadromous populations, lacustrine populations do not exhibit more deleterious mutations than the two other history types, suggesting that the negative effect of genetic drift in lacustrine populations may be mitigated by that of relaxed purifying selection. Moreover, we also identified genomic regions associated with anadromy, as well as an overrepresentation of transposable elements associated with variation in environmental variables, thus supporting the importance of transposable elements in adaptation.

Unraveling genetic population structure is challenging in species potentially characterized by large population size and high dispersal rates, often resulting in weak genetic differentiation. Genotyping a large number of samples can improve the detection of subtle genetic structure, but this may substantially increase sequencing cost and downstream bioinformatics computational time. To overcome this challenge, alternative, cost‐effective sequencing approaches, namely Pool‐seq and Rapture, have been developed. We empirically measured the power of resolution and congruence of these two methods in documenting weak population structure in nonmodel species with high gene flow comparatively to a conventional genotyping‐by‐sequencing (GBS) approach. For this, we used the American lobster (Homarus americanus) as a case study. First, we found that GBS, Rapture, and Pool‐seq approaches gave similar allele frequency estimates (i.e., correlation coefficient over 0.90) and all three revealed the same weak pattern of population structure. Yet, Pool‐seq data showed FST estimates three to five times higher than GBS and Rapture, while the latter two methods returned similar FST estimates, indicating that individual‐based approaches provided more congruent results than Pool‐seq. We conclude that despite higher costs, GBS and Rapture are more convenient approaches to use in the case of species exhibiting very weak differentiation. While both GBS and Rapture approaches provided similar results with regard to estimates of population genetic parameters, GBS remains more cost‐effective in project involving a relatively small numbers of genotyped individuals (e.g.

Because significant global changes are currently underway in the Arctic, creating a large-scale standardized database for Arctic marine biodiversity is particularly pressing. This study evaluates the potential of aquatic environmental DNA (eDNA) metabarcoding to detect Arctic coastal biodiversity changes and characterizes the local spatio-temporal distribution of eDNA in two locations. We extracted and amplified eDNA using two COI primer pairs from ~80 water samples that were collected across two Canadian Arctic ports, Churchill and Iqaluit, based on optimized sampling and preservation methods for remote regions surveys. Results demonstrate that aquatic eDNA surveys have the potential to document large-scale Arctic biodiversity change by providing a rapid overview of coastal metazoan biodiversity, detecting nonindigenous species, and allowing sampling in both open water and under the ice cover by local northern-based communities. We show that DNA sequences of ~50% of known Canadian Arctic species and potential invaders are currently present in public databases. A similar proportion of operational taxonomic units was identified at the species level with eDNA metabarcoding, for a total of 181 species identified at both sites. Despite the cold and well-mixed coastal environment, species composition was vertically heterogeneous, in part due to river inflow in the estuarine ecosystem, and differed between the water column and tide pools. Thus, COI-based eDNA metabarcoding may quickly improve large-scale Arctic biomonitoring using eDNA, but we caution that aquatic eDNA sampling needs to be standardized over space and time to accurately evaluate community structure changes., Ecology and Evolution, 8 (16), ISSN:2045-7758

Marine populations are typically characterized by weak genetic differentiation due to the potential for long-distance dispersal favouring high levels of gene flow. However, strong directional advection of water masses or retentive hydrodynamic forces can influence the degree of genetic exchange among marine populations. To determine the oceanographic drivers of genetic structure in a highly dispersive marine invertebrate, the giant California sea cucumber (Parastichopus californicus), we first tested for the presence of genetic discontinuities along the coast of North America in the northeastern Pacific Ocean. Then, we tested two hypotheses regarding spatial processes influencing population structure: (i) isolation by distance (IBD: genetic structure is explained by geographic distance) and (ii) isolation by resistance (IBR: genetic structure is driven by ocean circulation). Using RADseq, we genotyped 717 individuals from 24 sampling locations across 2,719 neutral SNPs to assess the degree of population differentiation and integrated estimates of genetic variation with inferred connectivity probabilities from a biophysical model of larval dispersal mediated by ocean currents. We identified two clusters separating north and south regions, as well as significant, albeit weak, substructure within regions (FST = 0.002, p = .001). After modelling the asymmetric nature of ocean currents, we demonstrated that local oceanography (IBR) was a better predictor of genetic variation (R2 = .49) than geographic distance (IBD) (R2 = .18), and directional processes played an important role in shaping fine-scale structure. Our study contributes to the growing body of literature identifying significant population structure in marine systems and has important implications for the spatial management of P. californicus and other exploited marine species.

Age at maturity is a key life‐history trait of most organisms. In anadromous salmonid fishes such as Atlantic Salmon (Salmo salar), age at sexual maturity is associated with sea age, the number of years spent at sea before the spawning migration. For the first time, we investigated the presence of two nonsynonymous vgll3 polymorphisms in North American Atlantic Salmon populations that relate to sea age in European salmon and quantified the natural variation at these and two additional candidate SNPs from two other genes. A targeted resequencing assay was developed and 1,505 returning adult individuals of size‐inferred sea age and sex from four populations were genotyped. Across three of four populations sampled in Québec, Canada, the late‐maturing component (MSW) of the population of a given sex exhibited higher proportions of SNP genotypes 54Thrvgll3 and 323Lysvgll3 compared to early‐maturing fish (1SW), for example, 85% versus 53% of females from Trinité River carried 323Lysvgll3 (n MSW = 205 vs. n 1SW = 30; p 66%) to be female. In summary, two nonsynonymous vgll3 polymorphisms were confirmed in North American populations of Atlantic Salmon and our results suggest that variation at those loci correlates with sea age and sex. Our results also suggest that this correlation varies among populations. Future work would benefit from a more balanced sampling and from adding data on juvenile riverine life stages to contrast our data.

Genomic GC content can vary locally, and GC-rich regions are usually associated with increased DNA thermostability in thermophilic prokaryotes and warm-blooded eukaryotes. Among vertebrates, fish and amphibians appeared to possess a distinctly less heterogeneous AT/GC organization in their genomes, whereas cytogenetically detectable GC heterogeneity has so far only been documented in mammals and birds. The subject of our study is the gar, an ancient "living fossil" of a basal ray-finned fish lineage, known from the Cretaceous period. We carried out cytogenomic analysis in two gar genera (Atractosteus and Lepisosteus) uncovering a GC chromosomal pattern uncharacteristic for fish. Bioinformatic analysis of the spotted gar (Lepisosteus oculatus) confirmed a GC compartmentalization on GC profiles of linkage groups. This indicates a rather mammalian mode of compositional organization on gar chromosomes. Gars are thus the only analyzed extant ray-finned fishes with a GC compartmentalized genome. Since gars are cold-blooded anamniotes, our results contradict the generally accepted hypothesis that the phylogenomic onset of GC compartmentalization occurred near the origin of amniotes. Ecophysiological findings of other authors indicate a metabolic similarity of gars with mammals. We hypothesize that gars might have undergone convergent evolution with the tetrapod lineages leading to mammals on both metabolic and genomic levels. Their metabolic adaptations might have left footprints in their compositional genome evolution, as proposed by the metabolic rate hypothesis. The genome organization described here in gars sheds new light on the compositional genome evolution in vertebrates generally and contributes to better understanding of the complexities of the mechanisms involved in this process.

During speciation-with-gene-flow, effective migration varies across the genome as a function of several factors, including proximity of selected loci, recombination rate, strength of selection, and number of selected loci. Genome scans may provide better empirical understanding of the genome-wide patterns of genetic differentiation, especially if the variance due to the previously mentioned factors is partitioned. In North American lake whitefish (Coregonus clupeaformis), glacial lineages that diverged in allopatry about 60,000 years ago and came into contact 12,000 years ago have independently evolved in several lakes into two sympatric species pairs (a normal benthic and a dwarf limnetic). Variable degrees of reproductive isolation between species pairs across lakes offer a continuum of genetic and phenotypic divergence associated with adaptation to distinct ecological niches. To disentangle the complex array of genetically based barriers that locally reduce the effective migration rate between whitefish species pairs, we compared genome-wide patterns of divergence across five lakes distributed along this divergence continuum. Using restriction site associated DNA (RAD) sequencing, we combined genetic mapping and population genetics approaches to identify genomic regions resistant to introgression and derive empirical measures of the barrier strength as a function of recombination distance. We found that the size of the genomic islands of differentiation was influenced by the joint effects of linkage disequilibrium maintained by selection on many loci, the strength of ecological niche divergence, as well as demographic characteristics unique to each lake. Partial parallelism in divergent genomic regions likely reflected the combined effects of polygenic adaptation from standing variation and independent changes in the genetic architecture of postzygotic isolation. This study illustrates how integrating genetic mapping and population genomics of multiple sympatric species pairs provide a window on the speciation-with-gene-flow mechanism.

The unique life-history characteristics of North Atlantic catadromous eels have long intrigued evolutionary biologists, especially with respect to mechanisms that could explain their persistence as two ecologically very similar but reproductively and geographically distinct species. Differential developmental schedules during young larval stages have commonly been hypothesized to represent such a key mechanism. We performed a comparative analysis of gene expression by means of microarray experiments with American and European eel leptocephali collected in the Sargasso Sea in order to test the alternative hypotheses of (1) differential timing of gene expression regulation during early development versus (2) species-specific differences in expression of particular genes. Our results provide much stronger support for the former hypothesis since no gene showed consistent significant differences in expression levels between the two species. In contrast, 146 genes showed differential timings of expression between species, although the observed expression level differences between the species were generally small. Consequently, species-specific gene expression regulation seems to play a minor role in species differentiation. Overall, these results show that the basis of the early developmental divergence between the American and European eel is probably influenced by differences in the timing of gene expression regulation for genes involved in a large array of biological functions.

Because of intrinsic differences in their genetic architectures, wild populations invaded by domesticated individuals could experience population-specific consequences following introgression by genetic material of domesticated origin. Expression levels of 16 000 transcripts were quantified by microarrays in liver tissue from farm, wild, and farm-wild backcross (i.e. F1 farm-wild hybrid × wild; total n = 50) Atlantic salmon (Salmo salar) raised under common environmental conditions. The wild populations and farm strain originated from three North American rivers in eastern Canada (Stewiacke, Tusket, and Saint John rivers, respectively). Analysis of variance revealed 177 transcripts with different expression levels among the five strains compared. Five times more of these transcripts were differentiated between farmed parents and Tusket backcrosses (n = 53) than between Stewiacke backcrosses and their farmed parents (n = 11). Altered biological processes in backcrosses also differed between populations both in number and in the type of processes impacted (metabolism vs immunity). Over-dominant gene expression regulation in backcrosses varied considerably between populations (23% in Stewiacke vs 44% in Tusket). Hence, the consequences of introgression of farm genetic material on gene expression depended on population-specific genetic architectures. These results support the need to evaluate impacts of farm-wild genetic interactions at the population scale.