Genome assembly, structural variants, and genetic differentiation between lake whitefish young species pairs (Coregonus sp.) with long and short reads.

Nascent pairs of ecologically differentiated species offer an opportunity to get a better glimpse at the genetic architecture of speciation. Of particular interest is our recent ability to consider a wider range of genomic variants, not only single-nucleotide polymorphisms (SNPs), thanks to long-read sequencing technology. We can now identify structural variants (SVs) such as insertions, deletions and other rearrangements, allowing further insights into the genetic architecture of speciation and how different types of variants are involved in species differentiation. Here, we investigated genomic patterns of differentiation between sympatric species pairs (Dwarf and Normal) belonging to the lake whitefish (Coregonus clupeaformis) species complex. We assembled the first reference genomes for both C. clupeaformis sp. Normal and C. clupeaformis sp. Dwarf, annotated the transposable elements and analysed the genomes in the light of related coregonid species. Next, we used a combination of long- and short-read sequencing to characterize SVs and genotype them at the population scale using genome-graph approaches, showing that SVs cover five times more of the genome than SNPs. We then integrated both SNPs and SVs to investigate the genetic architecture of species differentiation in two different lakes and highlighted an excess of shared outliers of differentiation. In particular, a large fraction of SVs differentiating the two species correspond to insertions or deletions of transposable elements (TEs), suggesting that TE accumulation may represent a key component of genetic divergence between the Dwarf and Normal species. Together, our results suggest that SVs may play an important role in speciation and that, by combining second- and third-generation sequencing, we now have the ability to integrate SVs into speciation genomics.

Environmental selection, rather than neutral processes, best explain regional patterns of diversity in a tropical rainforest fish.

To conserve the high functional and genetic variation in hotspots such as tropical rainforests, it is essential to understand the forces driving and maintaining biodiversity. We asked to what extent environmental gradients and terrain structure affect morphological and genomic variation across the wet tropical distribution of an Australian rainbowfish, Melanotaenia splendida splendida. We used an integrative riverscape genomics and morphometrics framework to assess the influence of these factors on both putative adaptive and non-adaptive spatial divergence. We found that neutral genetic population structure was largely explainable by restricted gene flow among drainages. However, environmental associations revealed that ecological variables had a similar power to explain overall genetic variation, and greater power to explain body shape variation, than the included neutral covariables. Hydrological and thermal variables were the strongest environmental predictors and were correlated with traits previously linked to heritable habitat-associated dimorphism in rainbowfishes. In addition, climate-associated genetic variation was significantly associated with morphology, supporting heritability of shape variation. These results support the inference of evolved functional differences among localities, and the importance of hydroclimate in early stages of diversification. We expect that substantial evolutionary responses will be required in tropical rainforest endemics to mitigate local fitness losses due to changing climates.

Locally Adaptive Inversions Modulate Genetic Variation at Different Geographic Scales in a Seaweed Fly.

Across a species range, multiple sources of environmental heterogeneity, at both small and large scales, create complex landscapes of selection, which may challenge adaptation, particularly when gene flow is high. One key to multidimensional adaptation may reside in the heterogeneity of recombination along the genome. Structural variants, like chromosomal inversions, reduce recombination, increasing linkage disequilibrium among loci at a potentially massive scale. In this study, we examined how chromosomal inversions shape genetic variation across a species range and ask how their contribution to adaptation in the face of gene flow varies across geographic scales. We sampled the seaweed fly Coelopa frigida along a bioclimatic gradient stretching across 10° of latitude, a salinity gradient, and a range of heterogeneous, patchy habitats. We generated a chromosome-level genome assembly to analyze 1,446 low-coverage whole genomes collected along those gradients. We found several large nonrecombining genomic regions, including putative inversions. In contrast to the collinear regions, inversions and low-recombining regions differentiated populations more strongly, either along an ecogeographic cline or at a fine-grained scale. These genomic regions were associated with environmental factors and adaptive phenotypes, albeit with contrasting patterns. Altogether, our results highlight the importance of recombination in shaping adaptation to environmental heterogeneity at local and large scales.

Landscape genomics of the American lobster (Homarus americanus).

In marine species experiencing intense fishing pressures, knowledge of genetic structure and local adaptation represent a critical information to assist sustainable management. In this study, we performed a landscape genomics analysis in the American lobster to investigate the issues pertaining to the consequences of making use of putative adaptive loci to reliably infer population structure and thus more rigorously delineating biological management units in marine exploited species. Toward this end, we genotyped 14,893 single nucleotide polymorphism (SNPs) in 4190 lobsters sampled across 96 sampling sites distributed along 1000 km in the northwest Atlantic in both Canada and the USA. As typical for most marine species, we observed a weak, albeit highly significant genetic structure. We also found that adaptive genetic variation allows detecting fine-scale population structure not resolved by neutral genetic variation alone. Using the recent genome assembly of the American lobster, we were able to map and annotate several SNPs located in functional genes potentially implicated in adaptive processes such as thermal stress response, salinity tolerance and growth metabolism pathways. Taken together, our study indicates that weak population structure in high gene flow systems can be resolved at various spatial scales, and that putatively adaptive genetic variation can substantially enhance the delineation of biological management units of marine exploited species.

Thermal adaptation rather than demographic history drives genetic structure inferred by copy number variants in a marine fish.

Increasing evidence shows that structural variants represent an overlooked aspect of genetic variation with consequential evolutionary roles. Among those, copy number variants (CNVs), including duplicated genomic regions and transposable elements (TEs), may contribute to local adaptation and/or reproductive isolation among divergent populations. Those mechanisms suppose that CNVs could be used to infer neutral and/or adaptive population genetic structure, whose study has been restricted to microsatellites, mitochondrial DNA and Amplified fragment length polymorphism markers in the past and more recently the use of single nucleotide polymorphisms (SNPs). Taking advantage of recent developments allowing CNV analysis from RAD-seq data, we investigated how variation in fitness-related traits, local environmental conditions and demographic history are associated with CNVs, and how subsequent copy number variation drives population genetic structure in a marine fish, the capelin (Mallotus villosus). We collected 1538 DNA samples from 35 sampling sites in the north Atlantic Ocean and identified 6620 putative CNVs. We found associations between CNVs and the gonadosomatic index, suggesting that six duplicated regions could affect female fitness by modulating oocyte production. We also detected 105 CNV candidates associated with water temperature, among which 20% corresponded to genomic regions located within the sequence of protein-coding genes, suggesting local adaptation to cold water by means of gene sequence amplification. We also identified 175 CNVs associated with the divergence of three previously defined parapatric glacial lineages, of which 24% were located within protein-coding genes, making those loci potential candidates for reproductive isolation. Lastly, our analyses unveiled a hierarchical, complex CNV population structure determined by temperature and local geography, which was in stark contrast to that inferred based on SNPs in a previous study. Our findings underline the complementarity of those two types of genomic variation in population genomics studies.

Demography, genetics, and decline of a spatially structured population of lekking bird.

Understanding the mechanisms underlying population decline is a critical challenge for conservation biologists. Both deterministic (e.g. habitat loss, fragmentation, and Allee effect) and stochastic (i.e. demographic and environmental stochasticity) demographic processes are involved in population decline. Simultaneously, a decrease of population size has far-reaching consequences for genetics of populations by increasing the risk of inbreeding and the strength of genetic drift, which together inevitably results in a loss of genetic diversity and a reduced effective population size ([Formula: see text]). These genetic factors may retroactively affect vital rates (a phenomenon coined 'inbreeding depression'), reduce population growth, and accelerate demographic decline. To date, most studies that have examined the demographic and genetic processes driving the decline of wild populations have neglected their spatial structure. In this study, we examined demographic and genetic factors involved in the decline of a spatially structured population of a lekking bird, the western capercaillie (Tetrao urogallus). To address this issue, we collected capture-recapture and genetic data over a 6-years period in the Vosges Mountains (France). Our study showed that the population of T. urogallus experienced a severe decline between 2010 and 2015. We did not detect any Allee effect on survival and recruitment. By contrast, individuals of both sexes dispersed to avoid small subpopulations, thus suggesting a potential behavioral response to a mate finding Allee effect. In parallel to this demographic decline, the population showed low levels of genetic diversity, high inbreeding and low effective population sizes at both subpopulation and population levels. Despite this, we did not detect evidence of inbreeding depression: neither adult survival nor recruitment were affected by individual inbreeding level. Our study underlines the benefit from combining demographic and genetic approaches to investigate processes that are involved in population decline.

Copy number variants outperform SNPs to reveal genotype-temperature association in a marine species.

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 account 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 nonmodel species and highlights the importance of considering structural variants alongside SNPs to enhance our understanding of ecological and evolutionary processes shaping adaptive population structure.

Latitudinal variation in climate-associated genes imperils range edge populations.

The ecological impacts of increasing global temperatures are evident in most ecosystems on Earth, but our understanding of how climatic variation influences natural selection and adaptive resilience across latitudes remains largely unknown. Latitudinal gradients allow testing general ecosystem-level theories relevant to climatic adaptation. We assessed differences in adaptive diversity of populations along a latitudinal region spanning highly variable temperate to subtropical climates. We generated and integrated information from environmental mapping, phenotypic variation and genome-wide data from across the geographical range of the rainbowfish Melanotaenia duboulayi, an emerging aquatic system for studies of climate change. We detected, after controlling for spatial population structure, strong interactions between genotypes and environment associated with variation in stream flow and temperature. Some of these hydroclimate-associated genes were found to interact within functional protein networks that contain genes of adaptive significance for projected future climates in rainbowfish. Hydroclimatic selection was also associated with variation in phenotypic traits, including traits known to affect fitness of rainbowfish exposed to different flow environments. Consistent with predictions from the "climatic variability hypothesis," populations exposed to extremes of important environmental variables showed stronger adaptive divergence and less variation in climate-associated genes compared to populations at the centre of the environmental gradient. Our findings suggest that populations that evolved at environmental range margins and at geographical range edges may be more vulnerable to changing climates, a finding with implications for predicting adaptive resilience and managing biodiversity under climate change.

Adaptive and maladaptive genetic diversity in small populations: Insights from the Brook Charr (Salvelinus fontinalis) case study.

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.

Absence of founder effect and evidence for adaptive divergence in a recently introduced insular population of white-tailed deer (Odocoileus virginianus).

Islands are generally colonized by few individuals which could lead to a founder effect causing loss of genetic diversity and rapid divergence by strong genetic drift. Insular conditions can also induce new selective pressures on populations. Here, we investigated the extent of genetic differentiation within a white-tailed deer (Odocoileus virginianus) population introduced on an island and its differentiation with its source mainland population. In response to their novel environmental conditions, introduced deer changed phenotypically from mainland individuals, therefore we investigated the genetic bases of the morphological differentiation. The study was conducted on Anticosti Island (Québec, Canada) where 220 individuals were introduced 120 years ago, resulting in a population size over 160,000 individuals. We used genotyping-by-sequencing (GBS) to generate 8,518 filtered high-quality SNPs and compared patterns of genetic diversity and differentiation between the continental and Anticosti Island populations. Clustering analyses indicated a single panmictic island population and no sign of isolation by distance. Our results revealed a weak, albeit highly significant, genetic differentiation between the Anticosti Island population and its source population (mean FST  = 0.005), which allowed a population assignment success of 93%. Also, the high genetic diversity maintained in the introduced population supports the absence of a strong founder effect due to the large number of founders followed by rapid population growth. We further used a polygenic approach to assess the genetic bases of the divergent phenotypical traits between insular and continental populations. We found loci related to muscular function and lipid metabolism, which suggested that these could be involved in local adaptation on Anticosti Island. We discuss these results in a harvest management context.

Shared ancestral polymorphisms and chromosomal rearrangements as potential drivers of local adaptation in a marine fish.

Gene flow has tremendous importance for local adaptation, by influencing the fate of de novo mutations, maintaining standing genetic variation and driving adaptive introgression. Furthermore, structural variation as chromosomal rearrangements may facilitate adaptation despite high gene flow. However, our understanding of the evolutionary mechanisms impending or favouring local adaptation in the presence of gene flow is still limited to a restricted number of study systems. In this study, we examined how demographic history, shared ancestral polymorphism, and gene flow among glacial lineages contribute to local adaptation to sea conditions in a marine fish, the capelin (Mallotus villosus). We first assembled a 490-Mbp draft genome of M. villosus to map our RAD sequence reads. Then, we used a large data set of genome-wide single nucleotide polymorphisms (25,904 filtered SNPs) genotyped in 1,310 individuals collected from 31 spawning sites in the northwest Atlantic. We reconstructed the history of divergence among three glacial lineages and showed that they probably diverged from 3.8 to 1.8 million years ago and experienced secondary contacts. Within each lineage, our analyses provided evidence for large Ne and high gene flow among spawning sites. Within the Northwest Atlantic lineage, we detected a polymorphic chromosomal rearrangement leading to the occurrence of three haplogroups. Genotype-environment associations revealed molecular signatures of local adaptation to environmental conditions prevailing at spawning sites. Our study also suggests that both shared polymorphisms among lineages, resulting from standing genetic variation or introgression, and chromosomal rearrangements may contribute to local adaptation in the presence of high gene flow.

Parallel epigenetic modifications induced by hatchery rearing in a Pacific salmon.

Wild stocks of Pacific salmonids have experienced sharp declines in abundance over the past century. Consequently, billions of fish are released each year for enhancing abundance and sustaining fisheries. However, the beneficial role of this widely used management practice is highly debated since fitness decrease of hatchery-origin fish in the wild has been documented. Artificial selection in hatcheries has often been invoked as the most likely explanation for reduced fitness, and most studies to date have focused on finding signatures of hatchery-induced selection at the DNA level. We tested an alternative hypothesis, that captive rearing induces epigenetic reprogramming, by comparing genome-wide patterns of methylation and variation at the DNA level in hatchery-reared coho salmon (Oncorhynchus kisutch) with those of their wild counterparts in two geographically distant rivers. We found a highly significant proportion of epigenetic variation explained by the rearing environment that was as high as the one explained by the river of origin. The differentially methylated regions show enrichment for biological functions that may affect the capacity of hatchery-born smolts to migrate successfully in the ocean. Shared epigenetic variation between hatchery-reared salmon provides evidence for parallel epigenetic modifications induced by hatchery rearing in the absence of genetic differentiation between hatchery and natural-origin fish for each river. This study highlights epigenetic modifications induced by captive rearing as a potential explanatory mechanism for reduced fitness in hatchery-reared salmon.

Kin-dependent dispersal influences relatedness and genetic structuring in a lek system.

Kin selection and dispersal play a critical role in the evolution of cooperative breeding systems. Limited dispersal increases relatedness in spatially structured populations (population viscosity), with the result that neighbours tend to be genealogical relatives. Yet the increase in neighbours' fitness-related performance through altruistic interaction may also result in habitat saturation and thus exacerbate local competition between kin. Our goal was to detect the footprint of kin selection and competition by examining the spatial structure of relatedness and by comparing non-effective and effective dispersal in a population of a lekking bird, Tetrao urogallus. For this purpose, we analysed capture-recapture and genetic data collected over a 6-year period on a spatially structured population of T. urogallus in France. Our findings revealed a strong spatial structure of relatedness in males. They also indicated that the population viscosity could allow male cooperation through two non-exclusive mechanisms. First, at their first lek attendance, males aggregate in a lek composed of relatives. Second, the distance corresponding to non-effective dispersal dramatically outweighed effective dispersal distance, which suggests that dispersers incur high post-settlement costs. These two mechanisms result in strong population genetic structuring in males. In females, our findings revealed a lower level of spatial structure of relatedness and genetic structure in respect to males. Additionally, non-effective dispersal and effective dispersal distances in females were highly similar, which suggests limited post-settlement costs. These results indicate that kin-dependent dispersal decisions and costs have a genetic footprint in wild populations and are factors that may be involved in the evolution of cooperative courtship.

Convergence in organ size but not energy metabolism enzyme activities among wild Lake Whitefish (Coregonus clupeaformis) species pairs.

The repeated evolution of similar phenotypes by similar mechanisms can be indicative of local adaptation, constraints or biases in the evolutionary process. Little is known about the incidence of physiological convergence in natural populations, so here we test whether energy metabolism in 'dwarf' and 'normal' Lake Whitefish evolves by similar mechanisms. Prior genomic and transcriptomic studies have found that divergence in energy metabolism is key to local adaptation in whitefish species pairs, but that distinct genetic and transcriptomic changes often underlie phenotypic evolution among lakes. Here, we predicted that traits at higher levels of biological organization, including the activities of energy metabolism enzymes (the product of enzyme concentration and turnover rate) and the relative proportions of metabolically active tissues (heart, liver, skeletal muscle), would show greater convergence than genetic and transcriptomic variation. We compared four whitefish species pairs and found convergence in organ size whereby all dwarf whitefish populations have a higher proportion of red skeletal muscle, three have relatively larger livers and two have relatively larger ventricles than normal fish. On the other hand, hepatic and muscle enzyme activities showed little convergence and were largely dependent on lake of origin. Only the most genetically divergent species pair (Cliff Lake) displayed white muscle enzyme activities matching results from laboratory-reared normal and dwarf whitefish. Overall, these data show convergence in the evolution of organ size, but not in the activities of candidate enzymes of energy metabolism, which may have evolved mainly as a consequence of demographic or ecological differences among lakes.

Adaptation and acclimation of traits associated with swimming capacity in Lake Whitefish (coregonus clupeaformis) ecotypes.

BACKGROUND: Improved performance in a given ecological niche can occur through local adaptation, phenotypic plasticity, or a combination of these mechanisms. Evaluating the relative importance of these two mechanisms is needed to better understand the cause of intra specific polymorphism. In this study, we reared populations of Lake Whitefish (Coregonus clupeaformis) representing the'normal' (benthic form) and the 'dwarf' (derived limnetic form) ecotypes in two different conditions (control and swim-training) to test the relative importance of adaptation and acclimation in the differentiation of traits related to swimming capacity. The dwarf whitefish is a more active swimmer than the normal ecotype, and also has a higher capacity for aerobic energy production in the swimming musculature. We hypothesized that dwarf fish would show changes in morphological and physiological traits consistent with reductions in the energetic costs of swimming and maintenance metabolism. RESULTS: We found differences in traits predicted to decrease the costs of prolonged swimming and standard metabolic rate and allow for a more active lifestyle in dwarf whitefish. Dwarf whitefish evolved a more streamlined body shape, predicted to lead to a decreased drag, and a smaller brain, which may decrease their standard metabolic rate. Contrary to predictions, we also found evidence of acclimation in liver size and metabolic enzyme activities. CONCLUSION: Results support the view that local adaptation has contributed to the genetically-based divergence of traits associated with swimming activity. Presence of post-zygotic barriers limiting gene flow between these ecotype pairs may have favoured repeated local adaptation to the limnetic niches.

Seascape genomics provides evidence for thermal adaptation and current-mediated population structure in American lobster (Homarus americanus).

Investigating how environmental features shape the genetic structure of populations is crucial for understanding how they are potentially adapted to their habitats, as well as for sound management. In this study, we assessed the relative importance of spatial distribution, ocean currents and sea surface temperature (SST) on patterns of putatively neutral and adaptive genetic variation among American lobster from 19 locations using population differentiation (PD) approaches combined with environmental association (EA) analyses. First, PD approaches (using bayescan, arlequin and outflank) found 28 outlier SNPs putatively under divergent selection and 9770 neutral SNPs in common. Redundancy analysis revealed that spatial distribution, ocean current-mediated larval connectivity and SST explained 31.7% of the neutral genetic differentiation, with ocean currents driving the majority of this relationship (21.0%). After removing the influence of spatial distribution, no SST were significant for putatively neutral genetic variation whereas minimum annual SST still had a significant impact and explained 8.1% of the putatively adaptive genetic variation. Second, EA analyses (using Pearson correlation tests, bayescenv and lfmm) jointly identified seven SNPs as candidates for thermal adaptation. Covariation at these SNPs was assessed with a spatial multivariate analysis that highlighted a significant temperature association, after accounting for the influence of spatial distribution. Among the 505 candidate SNPs detected by at least one of the three approaches, we discovered three polymorphisms located in genes previously shown to play a role in thermal adaptation. Our results have implications for the management of the American lobster and provide a foundation on which to predict how this species will cope with climate change.

Epigenetic and Genetic Differentiation Between Coregonus Species Pairs.

Phenotypic diversification is classically associated with genetic differentiation and gene expression variation. However, increasing evidence suggests that DNA methylation is involved in evolutionary processes due to its phenotypic and transcriptional effects. Methylation can increase mutagenesis and could lead to increased genetic divergence between populations experiencing different environmental conditions for many generations, though there has been minimal empirical research on epigenetically induced mutagenesis in diversification and speciation. Whitefish, freshwater members of the salmonid family, are excellent systems to study phenotypic diversification and speciation due to the repeated divergence of benthic-limnetic species pairs serving as natural replicates. Here we investigate whole genome genetic and epigenetic differentiation between sympatric benthic-limnetic species pairs in lake and European whitefish (Coregonus clupeaformis and Coregonus lavaretus) from four lakes (N = 64). We found considerable, albeit variable, genetic and epigenetic differences between species pairs. All SNP types were enriched at CpG sites supporting the mutagenic nature of DNA methylation, though C>T SNPs were most common. We also found an enrichment of overlaps between outlier SNPs with the 5% highest FST between species and differentially methylated loci. This could possibly represent differentially methylated sites that have caused divergent genetic mutations between species, or divergent selection leading to both genetic and epigenetic variation at these sites. Our results support the hypothesis that DNA methylation contributes to phenotypic divergence and mutagenesis during whitefish speciation.

Design and validation of a high-density single nucleotide polymorphism array for the Eastern oyster (Crassostrea virginica).

Dense single nucleotide polymorphism (SNP) arrays are essential tools for rapid high-throughput genotyping for many genetic analyses, including genomic selection and high-resolution population genomic assessments. We present a high-density (200 K) SNP array developed for the Eastern oyster (Crassostrea virginica), which is a species of significant aquaculture production and restoration efforts throughout its native range. SNP discovery was performed using low-coverage whole-genome sequencing of 435 F1 oysters from families from 11 founder populations in New Brunswick, Canada. An Affymetrix Axiom Custom array was created with 219,447 SNPs meeting stringent selection criteria and validated by genotyping more than 4,000 oysters across 2 generations. In total, 144,570 SNPs had a call rate >90%, most of which (96%) were polymorphic and were distributed across the Eastern oyster reference genome, with similar levels of genetic diversity observed in both generations. Linkage disequilibrium was low (maximum r2 ∼0.32) and decayed moderately with increasing distance between SNP pairs. Taking advantage of our intergenerational data set, we quantified Mendelian inheritance errors to validate SNP selection. Although most of SNPs exhibited low Mendelian inheritance error rates overall, with 72% of called SNPs having an error rate of <1%, many loci had elevated Mendelian inheritance error rates, potentially indicating the presence of null alleles. This SNP panel provides a necessary tool to enable routine application of genomic approaches, including genomic selection, in C. virginica selective breeding programs. As demand for production increases, this resource will be essential for accelerating production and sustaining the Canadian oyster aquaculture industry.

Genome-wide DNA methylation predicts environmentally driven life history variation in a marine fish.

Epigenetic modifications are thought to be one of the molecular mechanisms involved in plastic adaptive responses to environmental variation. However, studies reporting associations between genome-wide epigenetic changes and habitat-specific variations in life history traits (e.g., lifespan, reproduction) are still scarce, likely due to the recent application of methylome resequencing methods to non-model species. In this study, we examined associations between whole genome DNA methylation and environmentally driven life history variation in 2 lineages of a marine fish, the capelin (Mallotus villosus), from North America and Europe. In both lineages, capelin harbor 2 contrasting life history tactics (demersal vs. beach-spawning). Performing whole genome and methylome sequencing, we showed that life history tactics are associated with epigenetic changes in both lineages, though the effect was stronger in European capelin. Genetic differentiation between the capelin harboring different life history tactics was negligible, but we found genome-wide methylation changes in both lineages. We identified 9,125 European and 199 North American differentially methylated regions (DMRs) due to life history. Gene ontology (GO) enrichment analysis for both lineages revealed an excess of terms related to neural function. Our results suggest that environmental variation causes important epigenetic changes that are associated with contrasting life history tactics in lineages with divergent genetic backgrounds, with variable importance of genetic variation in driving epigenetic variation. Our study emphasizes the potential role of genome-wide epigenetic variation in adaptation to environmental variation.

Epigenetic variations are more substantial than genetic variations in rapid adaptation of oyster to Pacific oyster mortality syndrome.

Disease emergence is accelerating with global changes. Understanding by which mechanisms host populations can rapidly adapt will be crucial for management practices. Pacific oyster mortality syndrome (POMS) imposes a substantial and recurrent selective pressure on oyster populations, and rapid adaptation may arise through genetics and epigenetics. In this study, we used (epi)genome-wide association mapping to show that oysters differentially exposed to POMS displayed genetic and epigenetic signatures of selection. Consistent with higher resistance to POMS, the genes targeted included many genes in several pathways related to immunity. By combining correlation, DNA methylation quantitative trait loci, and variance partitioning, we revealed that a third of phenotypic variation was explained by interactions between the genetic and epigenetic information, ~14% by the genome, and up to 25% by the epigenome alone. Similar to genetically based adaptation, epigenetic mechanisms notably governing immune responses can contribute substantially to the rapid adaptation of hosts to emerging infectious diseases.