Search

Brook charr is a cold-water species which is highly sensitive to increased water temperatures, such as those associated with climate change. Environmental variation can potentially induce phenotypic changes that are inherited across generations, for instance, via epigenetic mechanisms. Here, we tested whether parental thermal regimes (intergenerational plasticity) and offspring-rearing temperatures (within-generational plasticity) modify the brain transcriptome of Brook charr progeny (fry stage). Parents were exposed to either cold or warm temperatures during final gonad maturation and their progeny were reared at 5 or 8 °C during the first stages of development. Illumina Novaseq6000 was used to sequence the brain transcriptome at the yolk sac resorption stage. The number of differentially expressed genes was very low when comparing fry reared at different temperatures (79 differentially expressed genes). In contrast, 9,050 differentially expressed genes were significantly differentially expressed between fry issued from parents exposed to either cold or warm temperatures. There was a significant downregulation of processes related to neural and synaptic activity in fry originating from the warm parental group vs fry from the cold parental one. We also observed significant upregulation of DNA methylation genes and of the most salient processes associated with compensation to warming, such as metabolism, cellular response to stress, and adaptive immunity. [ABSTRACT FROM AUTHOR]

Associations between host genotype and host‐associated microbiomes have been shown in a variety of animal clades, but studies on teleosts mostly show weak associations. Our study aimed to explore these relationships in four sympatric Serrasalmidae (i.e., piranha) teleosts from an Amazonian lake, using data sets from the hosts genomes (single nucleotide polymorphisms from genotyping by sequencing), skin and gut microbiomes (16S rRNA gene metataxonomics) and diets (COI metabarcoding) from the same fish individuals. First, we investigated whether there were significant covariations of microbiome and fish genotypes at the inter‐ and intraspecific levels. We also assessed the extent of covariation between Serrasalmidae diet and microbiome, to isolate genotypic from dietary effects on community structure. We observed a significant covariation of skin microbiomes and host genotypes at interspecific (R2 = 24.4%) and intraspecific (R2 = 6.2%) levels, whereas gut microbiomes correlated poorly with host genotypes. Serrasalmidae diet composition was significantly correlated to fish genotype only at the interspecific level (R2 = 5.4%), but did not covary with gut microbiome composition (Mantel R = −.04). Second, we investigated whether the study of interspecific differentiation could benefit from considering host‐associated microbial communities in addition to host genotypes. By using a nonmetric multidimensional scaling (NMDS) ordination‐based approach, we observed that ordinations from skin‐ and gut species‐specific bacterial biomarkers identified through a random forest algorithm could significantly increase the average interspecific differentiation detected through host genotype data alone. Although future studies encompassing additional species and environments are needed, our results suggest Serrasalmidae microbiomes could constitute an insightful trait to be considered when studying the interspecific differences between members of this clade. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

While Atlantic salmon (Salmo salar) of the northernmost American populations is alimentary, economically, and culturally important for Ungava Inuit communities (Nunavik, Canada) and might play a key role in the persistence of the species in a global warming context, many mysteries remain about those remote and atypical populations. Thus, our first aim was to document the genomic structure of the Nunavik populations. The second objective was to determine whether salmon only migrating to the estuary without reaching the sea, apparently unique to those populations, represent distinct populations from the typical anadromous salmons and subsequently explore the genetic basis of migratory life‐history tactics in the species. Finally, the third goal was to quantify the contribution of each genetically distinct population and life‐history tactic in the mixed‐stock subsistence fishery of the Koksoak R. estuary. We used Genotyping‐by‐Sequencing to genotype 14,061 single nucleotide polymorphisms in the genome of 248 individuals from 8 source populations and 280 individuals from the Koksoak estuary mixed‐stock fishery. Life‐history tactics were identified by a visual assessment of scales. Results show a hierarchical structure mainly influenced by isolation‐by‐distance with 7 populations out of the 8 studied rivers. While no obvious structure was detected between marine and estuarine salmon within the population, we have identified genomic regions putatively associated with those migration tactics. Finally, all salmon captured in the Koksoak estuary originated from the Koksoak drainage and mostly from 2 tributaries, but no inter‐annual variation in the contribution of these tributaries was found. Our results indicate, however, that both marine and estuarine salmon contribute substantially to estuarine fisheries and that there is inter‐annual variation in this contribution. These findings provide crucial information for the conservation of salmon populations in a rapidly changing ecosystem, as well as for fishery management to improve the food security of Inuit communities. [ABSTRACT FROM AUTHOR]

Differences between sexes in trait fitness optima can generate intralocus sexual conflicts that have the potential to maintain genetic diversity through balancing selection. However, these differences are unlikely to be associated with strong selective coefficients and are challenging to detect. Additionally, recent studies have highlighted that duplications on sexual chromosomes can create artifactual signals of intralocus sexual conflicts. Thus, testing the relationship between intralocus sexual conflicts and balancing selection requires stringent filtering of duplicated regions, and dedicated methods to detect loci with low levels of intersex differentiation. In this study, we investigated intralocus sexual conflicts in the three-spined stickleback using whole-genome sequencing (mean coverage = 12×) of 50 females and 49 males from an anadromous population in the St. Lawrence River, Québec, Canada. After stringent filtering of duplications from the sex chromosomes, we compared three methods to detect intralocus sexual conflicts. We found only two genomic regions under potential intralocus sexual conflict that also showed signals of balancing selection. Overall, our results suggest that most intralocus sexual conflicts do not drive long-term balancing selection and are most likely transient. [ABSTRACT FROM AUTHOR]

Recent local adaptation to pollution has been evidenced in several organisms inhabiting environments heavily contaminated by metals. Nevertheless, the molecular mechanisms underlying adaptation to high metal concentrations are poorly understood, especially in fishes. Yellow perch ( Perca flavescens ) populations from lakes in the mining area of Rouyn-Noranda (QC, Canada) have been faced with metal contamination for about 90 years. Here, we examine gene transcription patterns of fish reciprocally transplanted between a reference and a metal-contaminated lake and also fish caged in their native lake. After four weeks, 111 genes were differentially transcribed in metal-naïve fish transferred to the metal-contaminated lake, revealing a plastic response to metal exposure. Genes involved in the citric cycle and beta-oxidation pathways were under-transcribed, suggesting a potential strategy to mitigate the effects of metal stress by reducing energy turnover. However, metal-contaminated fish transplanted to the reference lake did not show any transcriptomic response, indicating a reduced plastic response capability to sudden reduction in metal concentrations. Moreover, the transcription of other genes, especially ones involved in energy metabolism, was affected by caging. Overall, our results highlight environmental stress response mechanisms in yellow perch at the transcriptomic level and support a rapid adaptive response to metal exposure through genetic assimilation. [ABSTRACT FROM AUTHOR]

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 r² ~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. [ABSTRACT FROM AUTHOR]

Next-generation sequencing (NGS) is revolutionising marker development and the rapidly increasing amount of transcriptomes published across a wide variety of taxa is providing valuable sequence databases for the identification of genetic markers without the need to generate new sequences. Microsatellites are still the most important source of polymorphic markers in ecology and evolution. Motivated by our long-term interest in the adaptive radiation of a non-model species complex of whitefishes (Coregonus spp.), in this study, we focus on microsatellite characterisation and multiplex optimisation using transcriptome sequences generated by Illumina® and Roche-454, as well as online databases of Expressed Sequence Tags (EST) for the study of whitefish evolution and demographic history. We identified and optimised 40 polymorphic loci in multiplex PCR reactions and validated the robustness of our analyses by testing several population genetics and phylogeographic predictions using 494 fish from five lakes and 2 distinct ecotypes. [ABSTRACT FROM AUTHOR]

Highlights: [•] The transcriptional responses of Perca flavescens to both metal and non metal stressors were measured with a 1000 candidate-gene microarray. [•] 475, 287 and 176 genes were differentially transcribed depending on temperature, Ni and Cd concentrations, respectively. [•] Genes involved in iron metabolism, transcriptional and translational processes, vitamin metabolism, blood coagulation, and calcium transport were impacted by metals. [•] The developed microarray contributes to a better characterization of the impact of different stressors on the transcriptome. [ABSTRACT FROM AUTHOR]

The evolution of reproductive isolation in an ecological context may involve multiple facets of species divergence on which divergent selection may operate. These include variation in quantitative phenotypic traits, regulation of gene expression, and differential transmission of particular allelic combinations. Thus, an integrative approach to the speciation process involves identifying the genetic basis of these traits, in order to understand how they are affected by divergent selection in nature and how they ultimately contribute to reproductive isolation. In the Lake Whitefish ( Coregonus clupeaformis), dwarf and normal species pairs sympatrically occur in several North American postglacial lakes. The limnetic dwarf whitefish distinguishes from its normal benthic relative by numerous life history, behavioural, morphological and gene expression traits, in relation with the exploitation of distinct ecological niches. Here, we have applied the RAD-Sequencing method to a hybrid backcross family to reconstruct a high-density genetic linkage map and perform QTL mapping in the Lake Whitefish. The 3061 cM map encompassed 3438 segregating RAD markers distributed over 40 linkage groups, for an average resolution of 0.89 cM. We mapped phenotypic and expression QTL underlying ecologically important traits as well as transmission ratio distortion QTL, and identified genomic regions harbouring clusters of such QTL. A narrow genomic region strongly associated with sex determination was also evidenced. Positional and functional information revealed in this study will be useful in ongoing population genomic studies to illuminate our understanding of the genomic architecture of reproductive isolation between whitefish species pairs. [ABSTRACT FROM AUTHOR]

Escaped domesticated individuals can introduce disadvantageous traits into wild populations due to both adaptive differences between population ancestors and human-induced changes during domestication. In contrast to their domesticated counterparts, some endangered wild Atlantic salmon populations encounter during their marine stage large amounts of suspended sediments, which may act as a selective agent. We used microarrays to elucidate quantitative transcriptional differences between a domesticated salmon strain, a wild population and their first-generation hybrids during their marine life stage, to describe transcriptional responses to natural suspended sediments, and to test for adaptive genetic variation in plasticity relating to a history of natural exposure or nonexposure to suspended sediments. We identified 67 genes differing in transcription level among salmon groups. Among these genes, processes related to energy metabolism and ion homoeostasis were over-represented, while genes contributing to immunity and actin-/myosin-related processes were also involved in strain differentiation. Domestic-wild hybrids exhibited intermediate transcription patterns relative to their parents for two-thirds of all genes that differed between their parents; however, genes deviating from additivity tended to have similar levels to those expressed by the wild parent. Sediments induced increases in transcription levels of eight genes, some of which are known to contribute to external or intracellular damage mitigation. Although genetic variation in plasticity did not differ significantly between groups after correcting for multiple comparisons, two genes (metallothionein and glutathione reductase) tended to be more plastic in response to suspended sediments in wild and hybrid salmon, and merit further examination as candidate genes under natural selection. [ABSTRACT FROM AUTHOR]

Our understanding of the genetic basis of local adaptation has recently benefited from the increased power to identify functional variants associated with environmental variables at the genome scale. However, it often remains challenging to determine whether locally adaptive alleles are actively maintained at intermediate frequencies by spatially varying selection. Here, we evaluate the extent to which this particular type of balancing selection explains the retention of adaptive genetic variation in the extreme situation of perfect panmixia, using the American eel (Anguilla rostrata) as a model. We first conducted a genome scan between two samples from opposite ends of a latitudinal environmental gradient using 454 sequencing of individually tagged cDNA libraries. Candidate SNPs were then genotyped in 992 individuals from 16 sampling sites at different life stages of the same cohort (including larvae from the Sargasso Sea, glass eels, and 1-year-old individuals) as well as in glass eels of the following cohort. Evidence for spatially varying selection was found at 13 loci showing correlations between allele frequencies and environmental variables across the entire species range. Simulations under a multiple-niche Levene's model using estimated relative fitness values among genotypes rarely predicted a stable polymorphic equilibrium at these loci. Our results suggest that some genetic-by-environment interactions detected in our study arise during the progress toward fixation of a globally advantageous allele with spatially variable effects on fitness. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Epigenetic inheritance can result in plastic responses to changing environments being faithfully transmitted to offspring. However, it remains unclear how epigenetic mechanisms such as DNA methylation can contribute to multigenerational acclimation and adaptation to environmental stressors. Brook charr (Salvelinus fontinalis), an economically important salmonid, is highly sensitive to thermal stress and is of conservation concern in the context of climate change. We studied the effects of temperature during parental sexual maturation and offspring rearing on whole-genome DNA methylation in brook charr juveniles (fry). Parents were split between warm and cold temperatures during sexual maturation, mated in controlled breeding designs, then offspring from each family were split between warm (8°C) and cold (5°C) rearing environments. Using whole-genome bisulfite sequencing, we found 188 differentially methylated regions (DMRs) due to parental maturation temperature after controlling for family structure. By contrast, offspring rearing temperature had a negligible effect on offspring methylation. Stable intergenerational inheritance of DNA methylation and minimal plasticity in progeny could result in the transmission of acclimatory epigenetic states to offspring, priming them for a warming environment. Our findings have implications pertaining to the role of intergenerational epigenetic inheritance in response to ongoing climate change. [ABSTRACT FROM AUTHOR]

Despite the commercial importance of Greenland Halibut (Reinhardtius hippoglossoides), important gaps still persist in our knowledge of this species, including its reproductive biology and sex determination mechanism. Here, we combined single-molecule sequencing of long reads (Pacific Sciences) with chromatin conformation capture sequencing (Hi-C) data to assemble the first chromosome-level reference genome for this species. The high-quality assembly encompassed more than 598 Megabases (Mb) assigned to 1594 scaffolds (scaffold N50?25Mb) with 96% of its total length distributed among 24 chromosomes. Investigation of the syntenic relationship with other economically important flatfish species revealed a high conservation of synteny blocks among members of this phylogenetic clade. Sex determination analysis revealed that similar to other teleost fishes, flatfishes also exhibit a high level of plasticity and turnover in sex determination mechanisms. A low-coverage whole-genome sequence analysis of 198 individuals revealed that Greenland Halibut possesses a male heterogametic XY system and several putative candidate genes implied in the sex determination of this species. Our study also suggests for the first time in flatfishes that a putative Y-autosomal fusion could be associated with a reduction of recombination typical of the early steps of sex chromosome evolution. [ABSTRACT FROM AUTHOR]

Environmental DNA (eDNA) metabarcoding has gained growing attention as a strategy for monitoring biodiversity in ecology. However, taxa identifications produced through metabarcoding require sophisticated processing of high‐throughput sequencing data from taxonomically informative DNA barcodes. Various sets of universal and taxon‐specific primers have been developed, extending the usability of metabarcoding across archaea, bacteria and eukaryotes. Accordingly, a multitude of metabarcoding data analysis tools and pipelines have also been developed. Often, several developed workflows are designed to process the same amplicon sequencing data, making it somewhat puzzling to choose one among the plethora of existing pipelines. However, each pipeline has its own specific philosophy, strengths and limitations, which should be considered depending on the aims of any specific study, as well as the bioinformatics expertise of the user. In this review, we outline the input data requirements, supported operating systems and particular attributes of thirty‐two amplicon processing pipelines with the goal of helping users to select a pipeline for their metabarcoding projects. [ABSTRACT FROM AUTHOR]

Most population genomic tools rely on accurate single nucleotide polymorphism (SNP) calling and filtering to meet their underlying assumptions. However, genomic complexity, resulting from structural variants, paralogous sequences, and repetitive elements, presents significant challenges in assembling contiguous reference genomes. Consequently, short-read resequencing studies can encounter mismapping issues, leading to SNPs that deviate from Mendelian expected patterns of heterozygosity and allelic ratio. In this study, we employed the ngsParalog software to identify such deviant SNPs in whole-genome sequencing (WGS) data with low (1.5×) to intermediate (4.8×) coverage for four species: Arctic Char (Salvelinus alpinus), Lake Whitefish (Coregonus clupeaformis), Atlantic Salmon (Salmo salar), and the American Eel (Anguilla rostrata). The analyses revealed that deviant SNPs accounted for 22% to 62% of all SNPs in salmonid datasets and approximately 11% in the American Eel dataset. These deviant SNPs were particularly concentrated within repetitive elements and genomic regions that had recently undergone rediploidization in salmonids. Additionally, narrow peaks of elevated coverage were ubiquitous along all four reference genomes, encompassed most deviant SNPs, and could be partially associated with transposons and tandem repeats. Including these deviant SNPs in genomic analyses led to highly distorted site frequency spectra, underestimated pairwise FST values, and overestimated nucleotide diversity. Considering the widespread occurrence of deviant SNPs arising from a variety of sources, their important impact in estimating population parameters, and the availability of effective tools to identify them, we propose that excluding deviant SNPs from WGS datasets is required to improve genomic inferences for a wide range of taxa and sequencing depths. [ABSTRACT FROM AUTHOR]

The American eel (Anguilla rostrata) has long been regarded as a panmictic fish and has been confirmed as such in the northern part of its range. In this paper, we tested for the first time whether panmixia extends to the tropical range of the species. To do so, we first assembled a reference genome (975 Mbp, 19 chromosomes) combining long (PacBio and Nanopore and short (Illumina paired‐end) reads technologies to support both this study and future research. To test for population structure, we estimated genotype likelihoods from low‐coverage whole‐genome sequencing of 460 American eels, collected at 21 sampling sites (in seven geographic regions) ranging from Canada to Trinidad and Tobago. We estimated genetic distance between regions, performed ADMIXTURE‐like clustering analysis and multivariate analysis, and found no evidence of population structure, thus confirming that panmixia extends to the tropical range of the species. In addition, two genomic regions with putative inversions were observed, both geographically widespread and present at similar frequencies in all regions. We discuss the implications of lack of genetic population structure for the species. Our results are key for the future genomic research in the American eel and the implementation of conservation measures throughout its geographic range. Additionally, our results can be applied to fisheries management and aquaculture of the species. [ABSTRACT FROM AUTHOR]

Identification of stocks and quantification of their relative contribution to recruitment are major objectives toward improving the management and conservation of marine exploited species. Next‐generation sequencing allows for thousands of genomic markers to be analyzed, which provides the resolution needed to address these questions in marine species with weakly differentiated populations. Greenland Halibut (Reinhardtius hippoglossoides) is one of the most important exploited demersal species throughout the North Atlantic, in particular in the Gulf of St. Lawrence, Canada. There, two nurseries are known, the St. Lawrence Estuary and the northern Anticosti Island, but their contribution to the renewal of stocks remains unknown. The goals of this study were (a) to document the genetic structure and (b) to estimate the contribution of the different identified breeding stocks to nurseries. We sampled 100 juveniles per nursery and 50 adults from seven sites ranging from Saguenay Fjord to offshore Newfoundland, with some sites sampled over two consecutive years in order to evaluate the temporal stability of the contribution. Our results show that after removing sex‐linked markers, the Estuary/Gulf of St. Lawrence represents a single stock which is genetically distinct from the Atlantic around Newfoundland (FST = 0.00146, p‐value =.001). Population assignment showed that recruitment in both nurseries is predominantly associated with the St. Lawrence stock. However, we found that the relative contribution of both stocks to the nurseries is temporally variable with 1% contribution of the Newfoundland stock one year but up to 33% for the second year, which may be caused by year‐to‐year variation in larval transport into the Gulf of St. Lawrence. This study serves as a model for the identification of stocks for fisheries resources in a context where few barriers to dispersal occur, in addition to demonstrating the importance of considering sex‐linked markers and temporal replicates in studies of population genomics. [ABSTRACT FROM AUTHOR]

Abstract: 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. [ABSTRACT FROM AUTHOR]

Coho salmon (Oncorhynchus kisutch) are a culturally and economically important species that return from multiyear ocean migrations to spawn in rivers that flow to the Northern Pacific Ocean. Southern stocks of coho salmon in Canada and the United States have significantly declined over the past quarter century, and unfortunately, conservation efforts have not reversed this trend. To assist in stock management and conservation efforts, we generated a chromosome-level genome assembly. We also resequenced the genomes of 83 coho salmon across the North American range to identify nucleotide variants and understand the demographic histories of these salmon by modeling effective population size from genome-wide data. From demographic history modeling, we observed reductions in effective population sizes between 3,750 and 8,000 years ago for several northern sampling sites, which may correspond to bottleneck events during recolonization after glacial retreat. [ABSTRACT FROM AUTHOR]

Large chromosomal rearrangements are thought to facilitate adaptation to heterogeneous environments by limiting genomic recombination. Indeed, inversions have been implicated in adaptation along environmental clines and in ecotype specialization. Here, we combine classical ecological studies and population genetics to investigate an inversion polymorphism previously documented in Europe among natural populations of the seaweed fly Coelopa frigida along a latitudinal cline in North America. We test if the inversion is present in North America and polymorphic, assess which environmental conditions modulate the inversion karyotype frequencies, and document the relationship between inversion karyotype and adult size. We sampled nearly 2000 flies from 20 populations along several environmental gradients to quantify associations of inversion frequencies to heterogeneous environmental variables. Genotyping and phenotyping showed a widespread and conserved inversion polymorphism between Europe and America. Variation in inversion frequency was significantly associated with environmental factors, with parallel patterns between continents, indicating that the inversion may play a role in local adaptation. The three karyotypes of the inversion are differently favoured across microhabitats and represent life-history strategies likely to be maintained by the collective action of several mechanisms of balancing selection. Our study adds to the mounting evidence that inversions are facilitators of adaptation and enhance within-species diversity. [ABSTRACT FROM AUTHOR]

Freshwater eels ( Anguilla sp.) have large economic, cultural, ecological and aesthetic importance worldwide, but they suffered more than 90% decline in global stocks over the past few decades. Proper genetic resources, such as sequenced, assembled and annotated genomes, are essential to help plan sustainable recoveries by identifying physiological, biochemical and genetic mechanisms that caused the declines or that may lead to recoveries. Here, we present the first sequenced genome of the American eel. This genome contained 305 043 contigs (N50 = 7397) and 79 209 scaffolds (N50 = 86 641) for a total size of 1.41 Gb, which is in the middle of the range of previous estimations for this species. In addition, protein-coding regions, including introns and flanking regions, are very well represented in the genome, as 95.2% of the 458 core eukaryotic genes and 98.8% of the 248 ultra-conserved subset were represented in the assembly and a total of 26 564 genes were annotated for future functional genomics studies. We performed a candidate gene analysis to compare three genes among all three freshwater eel species and, congruent with the phylogenetic relationships, Japanese eel ( A. japanica) exhibited the most divergence. Overall, the sequenced genome presented in this study is a crucial addition to the presently available genetic tools to help guide future conservation efforts of freshwater eels. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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, Canada. 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. [ABSTRACT FROM AUTHOR]

Inferring the genomic basis of local adaptation is a long‐standing goal of evolutionary biology. Beyond its fundamental evolutionary implications, such knowledge can guide conservation decisions for populations of conservation and management concern. Here, we investigated the genomic basis of local adaptation in the Coho salmon (Oncorhynchus kisutch) across its entire North American range. We hypothesized that extensive spatial variation in environmental conditions and the species' homing behaviour may promote the establishment of local adaptation. We genotyped 7829 individuals representing 217 sampling locations at more than 100,000 high‐quality RADseq loci to investigate how recombination might affect the detection of loci putatively under selection and took advantage of the precise description of the demographic history of the species from our previous work to draw accurate population genomic inferences about local adaptation. The results indicated that genetic differentiation scans and genetic–environment association analyses were both significantly affected by variation in recombination rate as low recombination regions displayed an increased number of outliers. By taking these confounding factors into consideration, we revealed that migration distance was the primary selective factor driving local adaptation and partial parallel divergence among distant populations. Moreover, we identified several candidate single nucleotide polymorphisms associated with long‐distance migration and altitude including a gene known to be involved in adaptation to altitude in other species. The evolutionary implications of our findings are discussed along with conservation applications. see also the Perspective by Rebekah L. Horn and Shawn R. Narum [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Whole genome duplication (WGD) can provide material for evolutionary innovation. Family Salmonidae is ideal for studying the effects of WGD as the ancestral salmonid underwent WGD relatively recently, ~65 Ma, then rediploidized and diversified. Extensive synteny between homologous chromosome arms occurs in extant salmonids, but each species has both conserved and unique chromosome arm fusions and fissions. Assembly of large, outbred eukaryotic genomes can be difficult, but structural rearrangements within such taxa can be investigated using linkage maps. RAD sequencing provides unprecedented ability to generate high-density linkage maps for nonmodel species, but can result in low numbers of homologous markers between species due to phylogenetic distance or differences in library preparation. Here, we generate a high-density linkage map (3,826 markers) for the Salvelinus genera (Brook Charr S. fontinalis), and then identify corresponding chromosome arms among the other available salmonid high-density linkage maps, including six species of Oncorhynchus, and one species for each of Salmo, Coregonus, and the nonduplicated sister group for the salmonids, Northern Pike Esoxlucius for identifying post-duplicated homeologs. To facilitate this process, we developed MAPCOMP to identify identical and proximate (i.e. nearby) markers between linkage maps using a reference genome of a related species as an intermediate, increasing the number of comparable markers between linkage maps by 5-fold. This enabled a characterization of the most likely history of retained chromosomal rearrangements post-WGD, and several conserved chromosomal inversions. Analyses of RADseq-based linkage maps from other taxa will also benefit from MAPCOMP, available at: https://github.com/enormandeau/mapcomp [ABSTRACT FROM AUTHOR]

Among vertebrates, herpetofauna has the highest proportion of declining species. Detection of environmental DNA (eDNA) is a promising method towards significantly increasing large-scale herpetological conservation efforts. However, the integration of eDNA results within a management framework requires an evaluation of the efficiency of the method in large natural environments and the calibration of eDNA surveys with the quantitative monitoring tools currently used by conservation biologists. Towards this end, we first developed species-specific primers to detect the wood turtle ( Glyptemys insculpta) a species at risk in Canada, by quantitative PCR (qPCR). The rate of eDNA detection obtained by qPCR was also compared to the relative abundance of this species in nine rivers obtained by standardized visual surveys in the Province of Québec (Canada). Second, we developed multi-species primers to detect North American amphibian and reptile species using eDNA metabarcoding analysis. An occurrence index based on the distribution range and habitat type was compared with the eDNA metabarcoding dataset from samples collected in seven lakes and five rivers. Our results empirically support the effectiveness of eDNA metabarcoding to characterize herpetological species distributions. Moreover, detection rates provided similar results to standardized visual surveys currently used to develop conservation strategies for the wood turtle. We conclude that eDNA detection rates may provide an effective semiquantitative survey tool, provided that assay calibration and standardization is performed. [ABSTRACT FROM AUTHOR]

The eastern-most members of the Salmo trutta species complex in the Aralo-Caspian Sea region were studied to infer their population genetic structure and biogeographic origin. A total of 68 individuals collected from Iranian endorheic inland basins (Namak and Urmia lakes), tributaries of the Caspian (Haraz, Kura, Samur, Volga, and Ural river drainages) and Aral (Amu River) seas, and the Baltic Sea basin were genotyped using 26,202 SNPs via Genotyping-by-Sequencing. The data were analyzed using admixture, discriminant analysis of principal components (DAPC), analysis of molecular variance (AMOVA), species tree, genetic differentiation (FST), allele frequency difference (AFD), and neighbor network approaches. Trout in the southern Caspian Sea basin differ from those of the western and northern Caspian Sea. Based on our results, the Lake Namak trout is divergent from the southern and western Caspian trout populations. Aral Sea and Lake Namak trouts likely originated from the northern and southern Caspian Sea populations, respectively. Although only few populations were considered in this study, six conservation/management units of trouts are proposed. [ABSTRACT FROM AUTHOR]

Conservation units (CUs) are important tools for supporting the implementation of standardized management practices for exploited species. Following the adoption of the Wild Salmon Policy in Canada, CUs were defined for Pacific salmon based on characteristics related to ecotype, life history and genetic variation using microsatellite markers as indirect measures of local adaptation. Genomic data sets have the potential to improve the definition of CUs by reducing variance around estimates of population genetic parameters, thereby increasing the power to detect more subtle patterns of population genetic structure and by providing an opportunity to incorporate adaptive information more directly with the identification of variants putatively under selection. We used one of the largest genomic data sets recently published for a nonmodel species, comprising 5662 individual Coho salmon (Oncorhynchus kisutch) from 149 sampling locations and a total of 24,542 high‐quality SNPs obtained using genotyping‐by‐sequencing and mapped to the Coho salmon reference genome to (1) evaluate the current delineation of CUs for Coho in Canada and (2) compare patterns of population structure observed using neutral and outlier loci from genotype–environment association analyses to determine whether separate CUs that capture adaptive diversity are needed. Our results reflected CU boundaries on the whole, with the majority of sampling locations managed in the same CU clustering together within genetic groups. However, additional groups that are not currently represented by CUs were also uncovered. We observed considerable overlap in the genetic clusters identified using neutral or candidate loci, indicating a general congruence in patterns of genetic variation driven by local adaptation and gene flow in this species. Consequently, we suggest that the current CU boundaries for Coho salmon are largely well‐suited for meeting the Canadian Wild Salmon Policy's objective of defining biologically distinct groups, but we highlight specific areas where CU boundaries may be refined. [ABSTRACT FROM AUTHOR]

Summary The two primary ways that species respond to heterogeneous environments is through local adaptation and phenotypic plasticity. The American eel ( Anguilla rostrata ) presents a paradox; despite inhabiting drastically different environments [ 1 ], the species is panmictic [ 2, 3 ]. Spawning takes place only in the southern Sargasso Sea in the Atlantic Ocean [ 1 ]. Then, the planktonic larvae (leptocephali) disperse to rearing locations from Cuba to Greenland, and juveniles colonize either freshwater or brackish/saltwater habitats, where they spend 3–25 years before returning to the Sargasso Sea to spawn as a panmictic species. Depending on rearing habitat, individuals exhibit drastically different ecotypes [ 4–6 ]. In particular, individuals rearing in freshwater tend to grow slowly and mature older and are more likely to be female in comparison to individuals that rear in brackish/saltwater [ 4, 6 ]. The hypothesis that phenotypic plasticity alone can account for all of the differences was not supported by three independent controlled experiments [ 7–10 ]. Here, we present a genome-wide association study that demonstrates a polygenic basis that discriminates these habitat-specific ecotypes belonging to the same panmictic population. We found that 331 co-varying loci out of 42,424 initially considered were associated with the divergent ecotypes, allowing a reclassification of 89.6%. These 331 SNPs are associated with 101 genes that represent vascular and morphological development, calcium ion regulation, growth and transcription factors, and olfactory receptors. Our results are consistent with divergent natural selection of phenotypes and/or genotype-dependent habitat choice by individuals that results in these genetic differences between habitats, occurring every generation anew in this panmictic species. [ABSTRACT FROM AUTHOR]

The capacity of an individual to battle infection is an important fitness determinant in wild vertebrate populations. The major histocompatibility complex ( MHC) genes are crucial for a host's adaptive immune system to detect pathogens. However, anthropogenic activities may disrupt natural cycles of co-evolution between hosts and pathogens. In this study, we investigated the dynamic sequence and expression variation of host parasite interactions in brook charr ( Salvelinus fontinalis) in a context of past human disturbance via population supplementation from domestic individuals. To do so, we developed a new method to examine selection shaping MHC diversity within and between populations and found a complex interplay between neutral and selective processes that varied between lakes that were investigated. We provided evidence for a lower introgression rate of domestic alleles and found that parasite infection increased with domestic genomic background of individuals. We also documented an association between individual MHC alleles and parasite taxa. Finally, longer cis-regulatory minisatellites were positively correlated with MHC II down-regulation and domestic admixture, suggesting that inadvertent selection during domestication resulted in a lower immune response capacity, through a trade-off between growth and immunity, which explained the negative selection of domestic alleles at least under certain circumstances. [ABSTRACT FROM AUTHOR]

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

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Bioinformatic analysis of eDNA metabarcoding data is a crucial step toward rigorously assessing biodiversity. Many programs are now available for each step of the required analyses, but their relative abilities at providing fast and accurate species lists have seldom been evaluated. We used simulated mock communities and real fish eDNA metabarcoding data to evaluate the performance of 13 bioinformatic programs and pipelines to retrieve fish occurrence and read abundance using the 12S mt rRNA gene marker. We used four indices to compare the outputs of each program with the simulated samples: sensitivity, F‐measure, root‐mean‐square error (RMSE) on read relative abundances, and execution time. We found marked differences among programs only for the taxonomic assignment step, both in terms of sensitivity, F‐measure and RMSE. Running time was highly different between programs for each step. The fastest programs with best indices for each step were assembled into a pipeline. We compared this pipeline to pipelines constructed from existing toolboxes (OBITools, Barque, and QIIME 2). Our pipeline and Barque obtained the best performance for all indices and appear to be better alternatives to highly used pipelines for analysing fish eDNA metabarcoding data when a complete reference database is available. Analysis on real eDNA metabarcoding data also indicated differences for taxonomic assignment and execution time only. This study reveals major differences between programs during the taxonomic assignment step. The choice of algorithm for the taxonomic assignment can have a significant impact on diversity estimates and should be made according to the objectives of the study. [ABSTRACT FROM AUTHOR]

Domestication has been practiced for centuries yet directed toward relatively few terrestrial crops and animals. While phenotypic and quantitative genetic changes associated with domestication have been amply documented, little is known about the molecular changes underlying the phenotypic evolution during the process. Here, we have investigated the brook charr (Salvelinus fontinalis) responses to artificial selection by means of transcriptional analysis of ∼32,000 cDNA features performed in both selected and control populations reared under identical environmental conditions during four generations. Our results indicate that selective breeding led to significant changes in the transcription of genes at the juvenile stage, where we observed 4.16% (156/3750) of differentially expressed genes between the two lines. No significant genes were revealed at the earlier life stage. Moreover, when comparing our results to those of previous studies on Atlantic salmon that compared lines that were selected for five to seven generations for similar traits (e.g., growth), genes with similar biological functions were found to be under selection in both studies. These observations indicate that (1) four generations of selection caused substantial changes in regulation of gene transcription between selected and control populations and (2) selective breeding for improving the same phenotypic traits (e.g., rapid growth) in brook charr and Atlantic salmon tended to select for the same changes in transcription profiles as the expression of a small and similar set of genes was affected by selection. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

The genus Rutilus is widespread in the western and central Palearctic region. In the Caspian Sea, the taxonomic status of different populations of Rutilus lacustris has been debated due to different sub-specific names attributed to each population. We genotyped 7,984 single nucleotide polymorphisms and sequenced the mitochondrial cytochrome C oxidase subunit I gene of 37 R. lacustris and Rutilus frisii from the southeast and southwest Caspian Sea and the Aras River in the Kura River drainage. We analysed data using clustering, Bayes factor delimitation, introgression, assessment of migration rate, and phylogenetic analyses. The results showed that the southeast and southwest Caspian Sea populations of R. lacustris were closely related, but highly differentiated from R. lacustris in the Aras River. The Bayes factor delimitation test supported the existence of three populations of R. lacustris in the studied area. Three hybrid individuals with mtDNA from Abramis brama or R. frisii and nuclear DNA from R. lacustris were detected. To protect R. lacustris in the southern Caspian Sea, we propose that the Aras River and sea-run R. lacustris be treated as two separate conservation units and the southern Caspian Sea R. lacustris populations should be viewed as two potentially separate management units. [ABSTRACT FROM AUTHOR]

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. see also the Perspective by Anna Tigano. [ABSTRACT FROM AUTHOR]

Many populations of Caspian Sea trout (Salmo caspius)—a nationally endangered species in Iran—have been extirpated or depleted due to anthropogenic impacts. The Lar National Park hosts large populations of Caspian Sea trout, which have not been subject to fisheries management programs before, but the population/s also face different human-related threats that may endanger their sustainability. A total of 357 Caspian Sea trout collected from different streams in Lar National Park were genotyped at 7978 filtered SNP using Genotyping-By-Sequencing to document population genetic structure and the contribution of each population/habitat to lake-run trout fisheries. Our results revealed a fine-scale population genetic structure, which is probably a product of factors including natural and artificial barriers to gene flow, geographic distance, and behavioral differences between resident and lake-run trout. Mixed-Stock Analyses revealed a high contribution from four panmictic populations of the national park to lake-run fish and almost no contribution from streams located in upper reaches or from streams with hydro-chemical or physical barriers. Our results highlighted the necessity for a more serious conservation plan for both the populations contributing greatly to lake fisheries and the highly diverged upstream populations due to their uniqueness. [ABSTRACT FROM AUTHOR]

A thorough reconstruction of historical processes is essential for a comprehensive understanding of the mechanisms shaping patterns of genetic diversity. Indeed, past and current conditions influencing effective population size have important evolutionary implications for the efficacy of selection, increased accumulation of deleterious mutations, and loss of adaptive potential. Here, we gather extensive genome-wide data that represent the extant diversity of the Coho salmon (Oncorhynchus kisutch) to address two objectives. We demonstrate that a single glacial refugium is the source of most of the present-day genetic diversity, with detectable inputs from a putative secondary micro-refugium. We found statistical support for a scenario whereby ancestral populations located south of the ice sheets expanded recently, swamping out most of the diversity from other putative micro-refugia. Demographic inferences revealed that genetic diversity was also affected by linked selection in large parts of the genome. Moreover, we demonstrate that the recent demographic history of this species generated regional differences in the load of deleterious mutations among populations, a finding that mirrors recent results from human populations and provides increased support for models of expansion load. We propose that insights from these historical inferences should be better integrated in conservation planning of wild organisms, which currently focuses largely on neutral genetic diversity and local adaptation, with the role of potentially maladaptive variation being generally ignored. Author summary: Reconstruction of a species' past demographic history from genetic data can highlight historical factors that have shaped the distribution of genetic diversity along its genome and its geographic range. Here, we combine genotyping-by-sequencing with demographic modelling to address these issues in the Coho salmon, a Pacific salmon of conservation concern in some parts of its range, notably in the south. Our demographic reconstructions reveal a linear decrease in genetic diversity toward the north of the species range, supporting the hypothesis of a northern route of postglacial recolonization from a single major southern refugium. As predicted by theory, we also observed a higher proportion of deleterious mutations in the most distant populations from this refugium. Beyond this general pattern, among-site variation in the proportion of deleterious mutations is consistent with different local trends in effective population sizes. Our results highlight the potential importance of understanding historical factors that have shaped geographic patterns of the distribution of deleterious mutations in order to implement effective management programs for the conservation of wild populations. Such fundamental knowledge of human historical demography is now having major impacts on health sciences, and we argue it is time to integrate such approaches in conservation science as well. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Fig 2 is incorrectly cropped, removing four points in the dataset. The grey shaded area along the regression line corresponds to the 95% confidence level interval obtained from the linear models applied to each dataset. [Extracted from the article]

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., <1,000). Overall, this study illustrates the complexity of estimating genetic differentiation and other summary statistics in complex biological systems characterized by large population size and migration rates. [ABSTRACT FROM AUTHOR]

Using massively parallel sequencing data from two species with different life history traits, American lobster ( Homarus americanus) and Arctic Char ( Salvelinus alpinus), we highlight how an unbalanced sex ratio in the samples and a few sex-linked markers may lead to false interpretations of population structure and thus to potentially erroneous management recommendations. Here, multivariate analyses revealed two genetic clusters separating samples by sex instead of by expected spatial variation: inshore and offshore locations in lobster, or east and west locations in Arctic Char. To further investigate this, we created several subsamples artificially varying the sex ratio in the inshore/offshore and east/west groups and then demonstrated that significant genetic differentiation could be observed despite panmixia in lobster, and that FST values were overestimated in Arctic Char. This pattern was due to 12 and 94 sex-linked markers driving differentiation for lobster and Arctic Char, respectively. Removing sex-linked markers led to nonsignificant genetic structure in lobster and a more accurate estimation of FST in Arctic Char. The locations of these markers and putative identities of genes containing or nearby the markers were determined using available transcriptomic and genomic data, and this provided new information related to sex determination in both species. Given that only 9.6% of all marine/diadromous population genomic studies to date have reported sex information, we urge researchers to collect and consider individual sex information. Sex information is therefore relevant for avoiding unexpected biases due to sex-linked markers as well as for improving our knowledge of sex determination systems in nonmodel species. [ABSTRACT FROM AUTHOR]

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 54Thr vgll3 and 323Lys vgll3 compared to early-maturing fish (1SW), for example, 85% versus 53% of females from Trinité River carried 323Lys vgll3 ( nMSW = 205 vs. n1SW = 30; p < .001). However, the association between vgll3 polymorphism and sea age was more pronounced in females than in males in the rivers we studied. Logistic regression analysis of vgll3 SNP genotypes revealed increased probabilities of exhibiting higher sea age for 54Thr vgll3 and 323Lys vgll3 genotypes compared to alternative genotypes, depending on population and sex. Moreover, individuals carrying the heterozygous vgll3 SNP genotypes were more likely (>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. [ABSTRACT FROM AUTHOR]