Global assessment of effective population sizes: Consistent taxonomic differences in meeting the 50/500 rule

Effective population size (Ne) is a particularly useful metric for conservation as it affects genetic drift, inbreeding and adaptive potential within populations. Current guidelines recommend a minimum Ne of 50 and 500 to avoid short-term inbreeding and to preserve long-term adaptive potential respectively. However, the extent to which wild populations reach these thresholds globally has not been investigated, nor has the relationship between Ne and human activities. Through a quantitative review, we generated a dataset with 4610 georeferenced Ne estimates from 3829 populations, extracted from 723 articles. These data show that certain taxonomic groups are less likely to meet 50/500 thresholds and are disproportionately impacted by human activities; plant, mammal and amphibian populations had a <54% probability of reachingN ̂ e $$ {\hat{N}}_e $$ = 50 and a <9% probability of reachingN ̂ e $$ {\hat{N}}_e $$ = 500. Populations listed as being of conservation concern according to the IUCN Red List had a smaller medianN ̂ e $$ {\hat{N}}_e $$ than unlisted populations, and this was consistent across all taxonomic groups.N ̂ e $$ {\hat{N}}_e $$ was reduced in areas with a greater Global Human Footprint, especially for amphibians, birds and mammals, however relationships varied between taxa. We also highlight several considerations for future works, including the role that gene flow and subpopulation structure plays in the estimation ofN ̂ e $$ {\hat{N}}_e $$ in wild populations, and the need for finer-scale taxonomic analyses. Our findings provide guidance for more specific thresholds based on Ne and help prioritise assessment of populations from taxa most at risk of failing to meet conservation thresholds.

Adaptation to seasonal reproduction and environment-associated factors drive temporal and spatial differentiation in northwest Atlantic herring despite gene flow

Understanding how marine organisms adapt to local environments is crucial for predicting how populations will respond to global climate change. The genomic basis, environmental factors and evolutionary processes involved in local adaptation are however not well understood. Here we use Atlantic herring, an abundant, migratory and widely distributed marine fish with substantial genomic resources, as a model organism to evaluate local adaptation. We examined genomic variation and its correlation with environmental variables across a broad environmental gradient, for 15 spawning aggregations in Atlantic Canada and the United States. We then compared our results with available genomic data of northeast Atlantic populations. We confirmed that population structure lies in a fraction of the genome including likely adaptive genetic variants of functional importance. We discovered 10 highly differentiated genomic regions distributed across four chromosomes. Nine regions show strong association with seasonal reproduction. One region, corresponding to a known inversion on chromosome 12, underlies a latitudinal pattern discriminating populations north and south of a biogeographic transition zone on the Scotian Shelf. Genome-environment associations indicate that winter seawater temperature best correlates with the latitudinal pattern of this inversion. The variation at two so-called 'islands of divergence' related to seasonal reproduction appear to be private to the northwest Atlantic. Populations in the northwest and northeast Atlantic share variation at four of these divergent regions, simultaneously displaying significant diversity in haplotype composition at another four regions, which includes an undescribed structural variant approximately 7.7 Mb long on chromosome 8. Our results suggest that the timing and geographic location of spawning and early development may be under diverse selective pressures related to allelic fitness across environments. Our study highlights the role of genomic architecture, ancestral haplotypes and selection in maintaining adaptive divergence in species with large population sizes and presumably high gene flow.

Epigenetic patterns in Atlantic herring (Clupea harengus): Temperature and photoperiod as environmental stressors during larval development

Understanding the molecular mechanisms underlying individual responses to environmental changes is crucial for species conservation and management. Pelagic fishes including Atlantic herring (Clupea harengus) are of particular interest because of their key ecological and economic roles and their susceptibility to a changing ocean from global warming. Temperature and photoperiod have been linked with spawning time and location in adult herring, but no study has thus far investigated the role of environmental factors on gene regulation during the vulnerable early developmental stages. Here, we examine DNA methylation patterns of larval herring bred under two temperatures (11°C and 13°C) and photoperiod (6 and 12 h) regimes in a 2 × 2 factorial design. We found consistently high levels of global methylation across all individuals and a decline in global methylation with increased developmental stage that was more pronounced at 13°C (p ≤ 0.007) than at 11°C (p ≥ 0.21). Most of the differentially methylated sites were in exon and promoter regions for genes linked to metabolism and development, some of which were hypermethylated at higher temperature. These results demonstrate the important role of DNA methylation during larval development and suggest that this molecular mechanism might be key in regulating early-stage responses to environmental stressors in Atlantic herring.

SNP Panel and Genomic Sex Identification in Atlantic Halibut (Hippoglossus hippoglossus)

The ability to identify sex is necessary in population biology for a proper understanding of the dynamics of a population. In Atlantic halibut, phenotypic sex identification is not possible due to the lack of significant external morphological differences. We developed an Illumina SNP panel for Atlantic halibut with 4000 SNPs spread evenly throughout the genome with a minor allele frequency MAF ≥ 0.4, except for N = 249 SNPs located in a sex-determining region on chromosome 12, N = 176 of these SNPs were selected to genetically identify male and female individuals using a DAPC analysis. The genomic identification of sex allows for non-lethal sex determination and validation of sex identification in the field. The SNP panel is a new genomic resource for Atlantic halibut that will make it possible to generate the genotypic data for the large number of individuals needed to estimate population abundance using genomics and the Close Kin Mark Recapture (CKMR) approach, an emerging component of fisheries management and stock monitoring.

Geography, environment, and colonization history interact with morph type to shape genomic variation in an Arctic fish

Polymorphic species are useful models for investigating the evolutionary processes driving diversification. Such processes include colonization history as well as contemporary selection, gene flow, and genetic drift, which can vary between intraspecific morphs as a function of their distinct life histories. The interactive and relative influence of such evolutionary processes on morph differentiation critically informs morph-specific management decisions and our understanding of incipient speciation. We therefore investigated how geographic distance, environmental conditions, and colonization history interacted with morph migratory capacity in the highly polymorphic fish species, Arctic Charr (Salvelinus alpinus). Using an 87k SNP chip we genetically characterized recently evolved anadromous, resident, and landlocked charr collected from 45 locations across a secondary contact zone of three charr glacial lineages in eastern Canada. A strong pattern of isolation by distance across all populations suggested geographic distance principally shaped genetic structure. Landlocked populations had lower genetic diversities and higher genetic differentiation than anadromous populations. However, effective population size was generally temporally stable in landlocked populations in comparison to anadromous populations. Genetic diversity positively correlated with latitude, potentially indicating southern anadromous populations' vulnerability to climate change and greater introgression between the Arctic and Atlantic glacial lineages in northern Labrador. Local adaptation was suggested by the observation of several environmental variables strongly associating with functionally relevant outlier genes including a region on chromosome AC21 potentially associated with anadromy. Our results demonstrate that gene flow, colonization history, and local adaptation uniquely interact to influence the genetic variation and evolutionary trajectory of populations.

Does the loss of diadromy imply the loss of salinity tolerance? A gene expression study with replicate nondiadromous populations of Galaxias maculatus

The recurrent colonization of freshwater habitats and subsequent loss of diadromy is a major ecological transition that has been reported in many ancestrally diadromous fishes. Such residency is often accompanied by a loss of tolerance to seawater. The amphidromous Galaxias maculatus has repeatedly colonized freshwater streams with evidence that freshwater-resident populations exhibit stark differences in their tolerance to higher salinities. Here, we used transcriptomics to gain insight into the mechanisms contributing to reduced tolerance to higher salinities in freshwater resident populations. We conducted an acute salinity challenge (0 ppt to 23-25 ppt) and measured osmoregulatory ability (muscle water content) over 48 h in three populations: diadromous, saltwater intolerant resident (Toltén), and saltwater tolerant resident (Valdivia). RNA sequencing of the gills identified genes that were differentially expressed in association with the salinity change and associated with the loss of saltwater tolerance in the Toltén population. Key genes associated with saltwater acclimation were characterized in diadromous G. maculatus individuals, some of which were also expressed in the saltwater tolerant resident population (Valdivia). We found that some of these "saltwater acclimation" genes, including the cystic fibrosis transmembrane conductance regulator gene (CFTR), were not significantly upregulated in the saltwater intolerant resident population (Toltén), suggesting a potential mechanism for the loss of tolerance to higher salinities. As the suite of differentially expressed genes in the diadromous-resident comparison differed between freshwater populations, we hypothesize that diadromy loss results in unique evolutionary trajectories due to drift, so the loss of diadromy does not necessarily lead to a loss in upper salinity tolerance.

Commerson's dolphin population structure: evidence for female phylopatry and male dispersal

A key in species conservation is understanding the amount and distribution of genetic diversity and how environmental changes that occurred in the recent past may have influenced current patterns of population structure. Commerson's dolphin, Cephalorhynchus commersonii, has two subspecies, one of which is endemic to South America (C. commersonii commersonii) and little is known about its population genetics. Our objective was to investigate the population genetics of this subspecies throughout its distribution. Using 70 skin samples and information available in GenBank, 308 mitochondrial DNA sequences and 28 species-specific microsatellites were analyzed. The species presented low genetic diversity when compared to other dolphin species, but was consistent with other species within the genus. Strong population structure based on mitochondrial DNA was exhibited throughout its entire distribution, a pattern consistent with female philopatry. However, this pattern was not detected when using microsatellites, suggesting male-mediated gene flow. Demographic tests suggested a population expansion beginning approximately 15,000 years ago, after the Last Glacial Maximum. In a climate change scenario, we recommended considering each sampling location as an independent population management unit in order to evaluate the impact of possible environmental changes on the distribution of genetic information within the species.

Low STR variability in the threatened marsh deer, Blastocerus dichotomus, detected through amplicon sequencing in non-invasive samples

Blastocerus dichotomus is the largest deer in South America. We have used 25 microsatellite markers detected and genotyped by Next Generation Sequencing to estimate the genetic variability of B. dichotomus in Argentina, where most of its populations are threatened. Primer design was based on the sequence of a shallow partial genome (15,967,456 reads; 16.66% genome coverage, mean depth 1.64) of a single individual. From the thousands of microsatellite loci found, even under high stringency selection, we chose and tested a set of 80 markers on 30 DNA samples extracted from tissue and feces from three Argentinean populations. Heterozygosity levels were low across all loci in all populations (H=0.31 to 0.40). Amplicon sequencing is a fast, easy, and affordable technique that can be very useful for the characterization of microsatellite marker sets for the conservation genetics of non-model organisms. This work is also one of the first ones to use amplicon sequencing in non-invasive samples and represents an important development for the study of threatened species.

Genetic Causes and Consequences of Sympatric Morph Divergence in Salmonidae: A Search for Mechanisms

Repeatedly and recently evolved sympatric morphs exhibiting consistent phenotypic differences provide natural experimental replicates of speciation. Because such morphs are observed frequently in Salmonidae, this clade provides a rare opportunity to uncover the genomic mechanisms underpinning speciation. Such insight is also critical for conserving salmonid diversity, the loss of which could have significant ecological and economic consequences. Our review suggests that genetic differentiation among sympatric morphs is largely nonparallel apart from a few key genes that may be critical for consistently driving morph differentiation. We discuss alternative levels of parallelism likely underlying consistent morph differentiation and identify several factors that may temper this incipient speciation between sympatric morphs, including glacial history and contemporary selective pressures. Our synthesis demonstrates that salmonids are useful for studying speciation and poses additional research questions to be answered by future study of this family. Expected final online publication date for the Annual Review of Animal Biosciences, Volume 10 is February 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Genomic basis of deep-water adaptation in Arctic Charr (Salvelinus alpinus) morphs

The post-glacial colonization of Gander Lake in Newfoundland, Canada, by Arctic Charr (Salvelinus alpinus) provides the opportunity to study the genomic basis of adaptation to extreme deep-water environments. Colonization of deep-water (>50 m) habitats often requires extensive adaptation to cope with novel environmental challenges from high hydrostatic pressure, low temperature, and low light, but the genomic mechanisms underlying evolution in these environments are rarely known. Here, we compare genomic divergence between a deep-water morph adapted to depths of up to 288 m and a larger, piscivorous pelagic morph occupying shallower depths. Using both a SNP array and resequencing of whole nuclear and mitochondrial genomes, we find clear genetic divergence (F_ST  = 0.11-0.15) between deep and shallow water morphs, despite an absence of morph divergence across the mitochondrial genome. Outlier analyses identified many diverged genomic regions containing genes enriched for processes such as gene expression and DNA repair, cardiac function, and membrane transport. Detection of putative copy number variants (CNVs) uncovered 385 genes with CNVs distinct to piscivorous morphs, and 275 genes with CNVs distinct to deep-water morphs, enriched for processes associated with synapse assembly. Demographic analyses identified evidence for recent and local morph divergence, and ongoing reductions in diversity consistent with postglacial colonization. Together, these results show that Arctic Charr morph divergence has occurred through genome-wide differentiation and elevated divergence of genes underlying multiple cellular and physiological processes, providing insight into the genomic basis of adaptation in a deep-water habitat following postglacial recolonization.

Population abundance in arctic grayling using genetics and close-kin mark-recapture

Arctic Grayling (Thymallus arcticus) are among the most widely distributed and abundant freshwater fish in the Yukon Territory of Canada, yet little information exists regarding their broad and fine-scale population structures or the number and size of these populations. The estimation of population abundance is fundamental for robust management and conservation, yet estimating abundance in the wild is often difficult. Here, we estimated abundance of an Arctic Grayling population using multiple genetic markers and the close-kin mark-recapture (CKMR) method. A total of N = 1,104 Arctic Grayling collected from two systems in Yukon were genotyped at 38 sequenced microsatellites. We first identified structure and assessed genetic diversity (effective population size,N ^ e ). Collections from one of the systems (Lubbock River) comprised adults and young-of-the-year sampled independently allowing the identification of parent-offspring pairs (POPs), and thus, the estimation of abundance using CKMR. We used COLONY and CKMRsim to identify POPs and both provided similar results leading to indistinguishable estimates (95% CI) of census size, that is,N ^ c ( C O L O N Y ) = 1858 (1259-2457) andN ^ c ( C K M R s i m ) = 1812 (1229-2389). The accuracy of the population abundance estimates can in the future be improved with temporal sampling and more precise age or size-specific fecundity estimates for Arctic Grayling. Our study demonstrates that the method can be used to inform management and conservation policy for Arctic Grayling and likely also for other fish species for which the assumption of random and independent sampling of adults and offspring can be assured.

Chromosome level reference of Atlantic halibut Hippoglossus hippoglossus provides insight into the evolution of sexual determination systems

Changes in the genetic mechanisms that control sexual determination have occurred independently across the tree of life, and with exceptional frequency in teleost fishes. To investigate the genomic changes underlying the evolution of sexual determination, we sequenced a chromosome-level genome, multitissue transcriptomes, and reduced representation population data for the Atlantic halibut (Hippoglossus hippoglossus), which has an XY/XX sex determination mechanism and has recently diverged (0.9-3.8 Ma) from the Pacific halibut (Hippoglossus stenolepis), which has a ZZ/ZW system. We used frequency and coverage-based population approaches to identify a putative sex-determining factor, GSDF. We characterized regions with elevated heterozygosity and linkage disequilibrium indicating suppression of recombination across a nascent sex chromosome. We detected testis-specific expression of GSDF, the sequence of which is highly conserved across flatfishes. Based on evidence from genome-wide association, coverage, linkage disequilibrium, testis and brain transcriptomes, and sequence conservation with other flatfishes, we propose a mechanism for the recent evolution of an XY sex-determination mechanism in Atlantic halibut. Changes to the ancestral sex-determining gene DMRT1 in regulating the downstream gene GSDF probably coincided with GSDF, or a proximal regulatory element of it, becoming the primary sex-determining factor. Our results suggest changes to a small number of elements can have drastic repercussions for the genomic substrate available to sex-specific evolutionary forces, providing insight into how certain elements repeatedly evolve to control sex across taxa. Our chromosome-level assembly, multitissue transcriptomes, and population genomic data provide a valuable resource and understanding of the evolution of sexual systems in fishes.

Functional genetic diversity in an exploited marine species and its relevance to fisheries management

The timing of reproduction influences key evolutionary and ecological processes in wild populations. Variation in reproductive timing may be an especially important evolutionary driver in the marine environment, where the high mobility of many species and few physical barriers to migration provide limited opportunities for spatial divergence to arise. Using genomic data collected from spawning aggregations of Pacific herring (Clupea pallasii) across 1600 km of coastline, we show that reproductive timing drives population structure in these pelagic fish. Within a specific spawning season, we observed isolation by distance, indicating that gene flow is also geographically limited over our study area. These results emphasize the importance of considering both seasonal and spatial variation in spawning when delineating management units for herring. On several chromosomes, we detected linkage disequilibrium extending over multiple Mb, suggesting the presence of chromosomal rearrangements. Spawning phenology was highly correlated with polymorphisms in several genes, in particular SYNE2, which influences the development of retinal photoreceptors in vertebrates. SYNE2 is probably within a chromosomal rearrangement in Pacific herring and is also associated with spawn timing in Atlantic herring (Clupea harengus). The observed genetic diversity probably underlies resource waves provided by spawning herring. Given the ecological, economic and cultural significance of herring, our results support that conserving intraspecific genetic diversity is important for maintaining current and future ecosystem processes.

Investigating Diadromy in Fishes and Its Loss in an -Omics Era

Diadromy, the predictable movements of individuals between marine and freshwater environments, is biogeographically and phylogenetically widespread across fishes. Thus, despite the high energetic and potential fitness costs involved in moving between distinct environments, diadromy appears to be an effective life history strategy. Yet, the origin and molecular mechanisms that underpin this migratory behavior are not fully understood. In this review, we aim first to summarize what is known about diadromy in fishes; this includes the phylogenetic relationship among diadromous species, a description of the main hypotheses regarding its origin, and a discussion of the presence of non-migratory populations within diadromous species. Second, we discuss how recent research based on -omics approaches (chiefly genomics, transcriptomics, and epigenomics) is beginning to provide answers to questions on the genetic bases and origin(s) of diadromy. Finally, we suggest future directions for -omics research that can help tackle questions on the evolution of diadromy.

Genomic basis of the loss of diadromy in Galaxias maculatus: Insights from reciprocal transplant experiments

Diadromy is known for having major effects on the distribution and richness of aquatic species, and so does its loss. The loss of diadromy has led to the diversification of many species, yet research focusing on understanding its molecular basis and consequences are limited. This is particularly true for amphidromous species despite being the most abundant group of diadromous species. Galaxias maculatus, an amphidromous species and one of the most widely distributed fishes in the Southern Hemisphere, exhibits many instances of nonmigratory or resident populations. The existence of naturally replicated resident populations in Patagonia can serve as an ideal system for the study of the mechanisms that lead to the loss of the diadromy and its ecological and evolutionary consequences. Here, we studied two adjacent river systems in which resident populations are genetically differentiated yet derived from the same diadromous population. By combining a reciprocal transplant experiment with genomic data, we showed that the two resident populations followed different evolutionary pathways by exhibiting a differential response in their capacity to survive in salt water. While one resident population was able to survive salt water, the other was not. Genomic analyses provided insights into the genes that distinguished (a) migratory from nonmigratory populations; (b) populations that can vs those that cannot survive a saltwater environment; and (c) between these resident populations. This study demonstrates that the loss of diadromy can be achieved by different pathways and that environmental (selection) and random (genetic drift) forces shape this dynamic evolutionary process.

Multiple drainage reversal episodes and glacial refugia in a Patagonian fish revealed by sequenced microsatellites

The rise of the southern Andes and the Quaternary glacial cycles influenced the landscape of Patagonia, affecting the phylogeographic and biogeographic patterns of its flora and fauna. Here, we examine the phylogeography of the freshwater fish, Percichthys trucha, using 53 sequenced microsatellite DNA markers. Fish (n= 835) were collected from 16 river systems (46 locations) spanning the species range on both sides of the Andes. Eleven watersheds drain to the Pacific, five of which are trans-Andean (headwaters east of Andes). The remaining five drainages empty into the Atlantic. Three analytical approaches (neighbour-joining tree, hierarchical AMOVAs, Structure) revealed evidence of historic drainage reversals: fish from four of the five trans-Andean systems (Puelo, Futalaufquen/Yelcho, Baker, Pascua) exhibited greater genetic similarity with Atlantic draining systems than with Pacific systems with headwaters west of Andes. Present-day drainage (Pacific versus Atlantic) explained only 5% of total genetic variance, while ancestral drainage explained nearly 27% of total variance. Thus, the phylogeographic structure of P. trucha is consistent with episodes of drainage reversal in multiple systems and suggests a major role for deglaciation in the genetic and indeed the geographical distribution of P. trucha in Patagonia. The study emphasizes the significant role of historical processes in the current pattern of genetic diversity and differentiation in a fish from a southern temperate region.

Resolving fine-scale population structure and fishery exploitation using sequenced microsatellites in a northern fish

The resiliency of populations and species to environmental change is dependent on the maintenance of genetic diversity, and as such, quantifying diversity is central to combating ongoing widespread reductions in biodiversity. With the advent of next-generation sequencing, several methods now exist for resolving fine-scale population structure, but the comparative performance of these methods for genetic assignment has rarely been tested. Here, we evaluate the performance of sequenced microsatellites and a single nucleotide polymorphism (SNP) array to resolve fine-scale population structure in a critically important salmonid in north eastern Canada, Arctic Charr (Salvelinus alpinus). We also assess the utility of sequenced microsatellites for fisheries applications by quantifying the spatial scales of movement and exploitation through genetic assignment of fishery samples to rivers of origin and comparing these results with a 29-year tagging dataset. Self-assignment and simulation-based analyses of 111 genome-wide microsatellite loci and 500 informative SNPs from 28 populations of Arctic Charr in north-eastern Canada identified largely river-specific genetic structure. Despite large differences (~4X) in the number of loci surveyed between panels, mean self-assignment accuracy was similar with the microsatellite loci and the SNP panel (>90%). Subsequent analysis of 996 fishery-collected samples using the microsatellite panel revealed that larger rivers contribute greater numbers of individuals to the fishery and that coastal fisheries largely exploit individuals originating from nearby rivers, corroborating results from traditional tagging experiments. Our results demonstrate the efficacy of sequence-based microsatellite genotyping to advance understanding of fine-scale population structure and harvest composition in northern and understudied species.

Human-induced habitat fragmentation effects on connectivity, diversity, and population persistence of an endemic fish, Percilia irwini, in the Biobío River basin (Chile)

An understanding of how genetic variability is distributed in space is fundamental for the conservation and maintenance of diversity in spatially fragmented and vulnerable populations. While fragmentation can occur from natural barriers, it can also be exacerbated by anthropogenic activities such as hydroelectric power plant development. Whatever the source, fragmentation can have significant ecological effects, including disruptions of migratory processes and gene flow among populations. In Chile, the Biobío River basin exhibits a high degree of habitat fragmentation due to the numerous hydroelectric power plants in operation, the number of which is expected to increase following new renewable energy use strategies. Here, we assessed the effects of different kinds of barriers on the genetic structure of the endemic freshwater fish Percilia irwini, knowledge that is critically needed to inform conservation strategies in light of current and anticipated further fragmentation initiatives in the system. We identified eight genetic units throughout the entire Biobío system with high effective sizes. A reduced effective size estimate was, however, observed in a single population located between two impassable barriers. Both natural waterfalls and human-made dams were important drivers of population differentiation in this system; however, dams affect genetic diversity differentially depending on their mode of operation. Evidence of population extirpation was found in two river stretches limited by upstream and downstream dams. Significant gene flow in both directions was found among populations not separated by natural or anthropogenic barriers. Our results suggest a significant vulnerability of P. irwini populations to future dam development and demonstrate the importance of studying basin-wide data sets with genetic metrics to understand the strength and direction of anthropogenic impacts on fish populations.

Invasive species and postglacial colonization: their effects on the genetic diversity of a Patagonian fish

The present distribution of Patagonian species is the result of a complex history involving Quaternary refugial populations, Holocene range expansions and demographic changes occurring during the Anthropocene. Invasive salmonids were introduced in Patagonia during the last century, occupying most rivers and lakes, preying on and competing with native species, including the fish Galaxias platei. Here, we used G. platei as a case study to understand how long-term (i.e. population differentiation during the Holocene) and short-term historical processes (salmonid introductions) affect genetic diversity. Using a suite of microsatellite markers, we found that the number of alleles is negatively correlated with the presence of salmonids (short-term processes), with G. platei populations from lakes with salmonids exhibiting significantly lower genetic diversity than populations from lakes without salmonids. Simulations (100 years backwards) showed that this difference in genetic diversity can be explained by a 99% reduction in population size. Allelic richness and observed heterozygosities were also negatively correlated with the presence of salmonids, but also positively correlated with long-term processes linked to Quaternary glaciations. Our results show how different genetic parameters can help identify processes taking place at different scales and their importance in terms of conservation.

Extensive secondary contact among three glacial lineages of Arctic Char (Salvelinus alpinus) in Labrador and Newfoundland

AIM

The Pleistocene glaciation event prompted the allopatric divergence of multiple glacial lineages of Arctic char (Salvelinus alpinus), some of which have come into secondary contact upon their recolonization of the Holarctic. While three glacial lineages (Arctic, Atlantic, and Acadian) are known to have recolonized the western Atlantic, the degree of overlap of these three lineages is largely unknown. We sought to determine the distribution of these three glacial lineages in Labrador and Newfoundland at a fine spatial scale to assess their potential for introgression and their relative contribution to local fisheries.

LOCATION

Labrador and Newfoundland, Canada.

METHODS

We sequenced a portion of the D-loop region in over 1,000 Arctic char (S. alpinus) samples from 67 locations across Labrador and Newfoundland.

RESULTS

Within Labrador, the Arctic and Atlantic lineages were widespread. Two locations (one landlocked and one with access to the sea) also contained individuals of the Acadian lineage, constituting the first record of this lineage in Labrador. Atlantic and Acadian lineage individuals were found in both eastern and western Newfoundland. Multiple sampling locations in Labrador and Newfoundland contained fish of two or more different glacial lineages, implying their introgression. Glacial lineage did not appear to dictate contemporary genetic divergence between the pale and dark morph of char present in Gander Lake, Newfoundland. Both were predominately of the Atlantic lineage, suggesting the potential for their divergence in sympatry.

MAIN CONCLUSIONS

Our study reveals Labrador and Newfoundland to be a unique junction of three glacial lineages which have likely hybridized extensively in this region.

Temporal stability and assignment power of adaptively divergent genomic regions between herring (Clupea harengus) seasonal spawning aggregations

Atlantic herring (Clupea harengus), a vital ecosystem component and target of the largest Northwest Atlantic pelagic fishery, undergo seasonal spawning migrations that result in elusive sympatric population structure. Herring spawn mostly in fall or spring, and genomic differentiation was recently detected between these groups. Here we used a subset of this differentiation, 66 single nucleotide polymorphisms (SNPs) to analyze the temporal dynamics of this local adaptation and the applicability of SNP subsets in stock assessment. We showed remarkable temporal stability of genomic differentiation corresponding to spawning season, between samples taken a decade apart (2005 N = 90 vs. 2014 N = 71) in the Gulf of St. Lawrence, and new evidence of limited interbreeding between spawning components. We also examined an understudied and overexploited herring population in Bras d'Or lake (N = 97); using highly reduced SNP panels (N SNPs > 6), we verified little-known sympatric spawning populations within this unique inland sea. These results describe consistent local adaptation, arising from asynchronous reproduction in a migratory and dynamic marine species. Our research demonstrates the efficiency and precision of SNP-based assessments of sympatric subpopulations; and indeed, this temporally stable local adaptation underlines the importance of such fine-scale management practices.

Past, present, and future of a freshwater fish metapopulation in a threatened landscape

It is well documented that hydropower plants can affect the dynamics of fish populations through landscape alterations and the creation of new barriers. Less emphasis has been placed on the examination of the genetic consequences for fish populations of the construction of dams. The relatively few studies that focus on genetics often do not consider colonization history and even fewer tend to use this information for conservation purposes. As a case study, we used a 3-pronged approach to study the influence of historical processes, contemporary landscape features, and potential future anthropogenic changes in landscape on the genetic diversity of a fish metapopulation. Our goal was to identify the metapopulation's main attributes, detect priority areas for conservation, and assess the consequences of the construction of hydropower plants for the persistence of the metapopulation. We used microsatellite markers and coalescent approaches to examine historical colonization processes, traditional population genetics, and simulations of future populations under alternate scenarios of population size reduction and gene flow. Historical gene flow appeared to have declined relatively recently and contemporary populations appeared highly susceptible to changes in landscape. Gene flow is critical for population persistence. We found that hydropower plants could lead to a rapid reduction in number of alleles and to population extirpation 50-80 years after their construction. More generally, our 3-pronged approach for the analyses of empirical genetic data can provide policy makers with information on the potential impacts of landscape changes and thus lead to more robust conservation efforts.

Whole-genome sequencing approaches for conservation biology: Advantages, limitations and practical recommendations

Whole-genome resequencing (WGR) is a powerful method for addressing fundamental evolutionary biology questions that have not been fully resolved using traditional methods. WGR includes four approaches: the sequencing of individuals to a high depth of coverage with either unresolved or resolved haplotypes, the sequencing of population genomes to a high depth by mixing equimolar amounts of unlabelled-individual DNA (Pool-seq) and the sequencing of multiple individuals from a population to a low depth (lcWGR). These techniques require the availability of a reference genome. This, along with the still high cost of shotgun sequencing and the large demand for computing resources and storage, has limited their implementation in nonmodel species with scarce genomic resources and in fields such as conservation biology. Our goal here is to describe the various WGR methods, their pros and cons and potential applications in conservation biology. WGR offers an unprecedented marker density and surveys a wide diversity of genetic variations not limited to single nucleotide polymorphisms (e.g., structural variants and mutations in regulatory elements), increasing their power for the detection of signatures of selection and local adaptation as well as for the identification of the genetic basis of phenotypic traits and diseases. Currently, though, no single WGR approach fulfils all requirements of conservation genetics, and each method has its own limitations and sources of potential bias. We discuss proposed ways to minimize such biases. We envision a not distant future where the analysis of whole genomes becomes a routine task in many nonmodel species and fields including conservation biology.

RAPID BEHAVIORAL CHANGES IN MEDAKA (ORYZIAS LATIPES) CAUSED BY SELECTION FOR COMPETITIVE AND NONCOMPETITIVE GROWTH

Behavior is a major component of growth rate variation in fish, but the genetic relationship between behavior and growth is largely unknown. In particular, it is not known how behavior responds to natural or artificial selection on growth. It has been argued that artificial selection for fast growth in competitive environments might lead to higher levels of overall aggression, and therefore there would be no net gain in assimilation efficiency or growth. The contrary has also been argued. It is not immediately obvious, therefore, whether fish that avoid interacting with conspecifics should grow faster, or whether more aggressive, dominant individuals should grow faster. In the present study, we artificially selected fish on the basis of growth in two environments that differed in the intensity of social interactions. In the high interaction environment (HI), food was provided to excess inside a floating cork ring, which limited access to the food and allowed fish to attempt to monopolize the food supply. In the low interaction environment (LI), an equivalent amount of food was spread over the container's surface. In each social environment, fish were selected within family for fast and for slow growth rate during two generations (G₁ and G₂ ). The response was measured in the third generation (G₃ ) on the directly selected trait (growth) and on a correlated trait (agonistic behavior). The magnitude of the direct response to selection on growth in medaka was directly related to population density, i.e., selection was most effective at high densities. Agonistic behavior was inversely related to growth when individuals were raised and selected in an environment where enforced social interaction took place, so long as food was not limited in quantity.

Parallel adaptive evolution of geographically distant herring populations on both sides of the North Atlantic Ocean

Atlantic herring is an excellent species for studying the genetic basis of adaptation in geographically distant populations because of its characteristically large population sizes and low genetic drift. In this study we compared whole-genome resequencing data of Atlantic herring populations from both sides of the Atlantic Ocean. An important finding was the very low degree of genetic differentiation among geographically distant populations (fixation index = 0.026), suggesting lack of reproductive isolation across the ocean. This feature of the Atlantic herring facilitates the detection of genetic factors affecting adaptation because of the sharp contrast between loci showing genetic differentiation resulting from natural selection and the low background noise resulting from genetic drift. We show that genetic factors associated with timing of reproduction are shared between genetically distinct and geographically distant populations. The genes for thyroid-stimulating hormone receptor (TSHR), the SOX11 transcription factor (SOX11), calmodulin (CALM), and estrogen receptor 2 (ESR2A), all with a significant role in reproductive biology, were among the loci that showed the most consistent association with spawning time throughout the species range. In fact, the same two SNPs located at the 5' end of TSHR showed the most significant association with spawning time in both the east and west Atlantic. We also identified unexpected haplotype sharing between spring-spawning oceanic herring and autumn-spawning populations across the Atlantic Ocean and the Baltic Sea. The genomic regions showing this pattern are unlikely to control spawning time but may be involved in adaptation to ecological factor(s) shared among these populations.

A portrait of a sucker using landscape genetics: how colonization and life history undermine the idealized dendritic metapopulation

Dendritic metapopulations have been attributed unique properties by in silico studies, including an elevated genetic diversity relative to a panmictic population of equal total size. These predictions have not been rigorously tested in nature, nor has there been full consideration of the interacting effects among contemporary landscape features, colonization history and life history traits of the target species. We tested for the effects of dendritic structure as well as the relative importance of life history, environmental barriers and historical colonization on the neutral genetic structure of a longnose sucker (Catostomus catostomus) metapopulation in the Kogaluk watershed of northern Labrador, Canada. Samples were collected from eight lakes, genotyped with 17 microsatellites, and aged using opercula. Lakes varied in differentiation, historical and contemporary connectivity, and life history traits. Isolation by distance was detected only by removing two highly genetically differentiated lakes, suggesting a lack of migration-drift equilibrium and the lingering influence of historical factors on genetic structure. Bayesian analyses supported colonization via the Kogaluk's headwaters. The historical concentration of genetic diversity in headwaters inferred by this result was supported by high historical and contemporary effective sizes of the headwater lake, T-Bone. Alternatively, reduced allelic richness in headwaters confirmed the dendritic structure's influence on gene flow, but this did not translate to an elevated metapopulation effective size. A lack of equilibrium and upstream migration may have dampened the effects of dendritic structure. We suggest that interacting historical and contemporary factors prevent the achievement of the idealized traits of a dendritic metapopulation in nature.

Scorched mussels (Brachidontes spp., Bivalvia: Mytilidae) from the tropical and warm-temperate southwestern Atlantic: the role of the Amazon River in their speciation

Antitropicality is a distribution pattern where closely related taxa are separated by an intertropical latitudinal gap. Two potential examples include Brachidontes darwinianus (south eastern Brazil to Uruguay), considered by some authors as a synonym of B. exustus (Gulf of Mexico and the Caribbean), and B. solisianus, distributed along the Brazilian coast with dubious records north of the intertropical zone. Using two nuclear (18S and 28S rDNA) and one mitochondrial gene (mtDNA COI), we aimed to elucidate the phylogeographic and phylogenetic relationships among the scorched mussels present in the warm-temperate region of the southwest Atlantic. We evaluated a divergence process mediated by the tropical zone over alternative phylogeographic hypotheses. Brachidontes solisianus was closely related to B. exustus I, a species with which it exhibits an antitropical distribution. Their divergence time was approximately 2.6 Ma, consistent with the intensification of Amazon River flow. Brachidontes darwinianus, an estuarine species is shown here not to be related to this B. exustus complex. We suspect ancestral forms may have dispersed from the Caribbean to the Atlantic coast via the Trans-Amazonian seaway (Miocene). The third species, B rodriguezii is presumed to have a long history in the region with related fossil forms going back to the Miocene. Although scorched mussels are very similar in appearance, their evolutionary histories are very different, involving major historical contingencies as the formation of the Amazon River, the Panama Isthmus, and the last marine transgression.

Effective number of breeders, effective population size and their relationship with census size in an iteroparous species, Salvelinus fontinalis

The relationship between the effective number of breeders (Nb) and the generational effective size (Ne) has rarely been examined empirically in species with overlapping generations and iteroparity. Based on a suite of 11 microsatellite markers, we examine the relationship between Nb, Ne and census population size (Nc) in 14 brook trout (Salvelinus fontinalis) populations inhabiting 12 small streams in Nova Scotia and sampled at least twice between 2009 and 2015. Unbiased estimates of Nb obtained with individuals of a single cohort, adjusted on the basis of age at first maturation (α) and adult lifespan (AL), were from 1.66 to 0.24 times the average estimates of Ne obtained with random samples of individuals of mixed ages (i.e. [Formula: see text]). In turn, these differences led to adjusted Ne estimates that were from nearly five to 0.7 times the estimates derived from mixed-aged individuals. These differences translate into the same range of variation in the ratio of effective to census population size [Formula: see text] within populations. Adopting [Formula: see text] as the more precise and unbiased estimates, we found that these brook trout populations differ markedly in their effective to census population sizes (range approx. 0.3 to approx. 0.01). Using AgeNe, we then showed that the variance in reproductive success or reproductive skew varied among populations by a factor of 40, from Vk/k ≈ 5 to 200. These results suggest wide differences in population dynamics, probably resulting from differences in productivity affecting the intensity of competition for access to mates or redds, and thus reproductive skew. Understanding the relationship between Ne, Nb and Nc, and how these relate to population dynamics and fluctuations in population size, are important for the design of robust conservation strategies in small populations with overlapping generations and iteroparity.

The complete mitochondrial genome of the freshwater fish Galaxias platei and a comparison with other species of the genus Galaxias (faraway, so close?)

We present the whole mitochondrial genome for Galaxias platei, a freshwater fish widely distributed throughout the Patagonian Andes, and compare it with the mitochondrial genome of three congeneric species. The position of G. platei in the phylogenetic reconstruction differs from that shown by earlier studies using multiple markers. We discuss the results in terms of the phylogenetic position of G. platei and the use of whole mitochondrial genomes versus specific regions of multiple organelles.

Echoes of a distant time: effects of historical processes on contemporary genetic patterns in Galaxias platei in Patagonia

Interpreting the genetic structure of a metapopulation as the outcome of gene flow over a variety of timescales is essential for the proper understanding of how changes in landscape affect biological connectivity. Here we contrast historical and contemporary connectivity in two metapopulations of the freshwater fish Galaxias platei in northern and southernmost Patagonia where paleolakes existed during the Holocene and Pleistocene, respectively. Contemporary gene flow was mostly high and asymmetrical in the northern system while extremely reduced in the southernmost system. Historical migration patterns were high and symmetric in the northern system and high and largely asymmetric in the southern system. Both systems showed a moderate structure with a clear pattern of isolation by distance (IBD). Effective population sizes were smaller in populations with low contemporary gene flow. An approximate Bayesian computation (ABC) approach suggests a late Holocene colonization of the lakes in the northern system and recent divergence of the populations from refugial populations from east and west of the Andes. For the southern system, the ABC approach reveals that some of the extant G. platei populations most likely derive from an ancestral population inhabiting a large Pleistocene paleolake while the rest derive from a higher-altitude lake. Our results suggest that neither historical nor contemporary processes individually fully explain the observed structure and geneflow patterns and both are necessary for a proper understanding of the factors that affect diversity and its distribution. Our study highlights the importance of a temporal perspective on connectivity to analyse the diversity of spatially complex metapopulations.

Contemporary effective population and metapopulation size (N e and meta-N e): comparison among three salmonids inhabiting a fragmented system and differing in gene flow and its asymmetries

We estimated local and metapopulation effective sizes ([Formula: see text] and meta-[Formula: see text]) for three coexisting salmonid species (Salmo salar, Salvelinus fontinalis, Salvelinus alpinus) inhabiting a freshwater system comprising seven interconnected lakes. First, we hypothesized that [Formula: see text] might be inversely related to within-species population divergence as reported in an earlier study (i.e., FST: S. salar> S. fontinalis> S. alpinus). Using the approximate Bayesian computation method implemented in ONeSAMP, we found significant differences in [Formula: see text] ([Formula: see text]) between species, consistent with a hierarchy of adult population sizes ([Formula: see text]). Using another method based on a measure of linkage disequilibrium (LDNE: [Formula: see text]), we found more finite [Formula: see text] values for S. salar than for the other two salmonids, in line with the results above that indicate that S. salar exhibits the lowest [Formula: see text] among the three species. Considering subpopulations as open to migration (i.e., removing putative immigrants) led to only marginal and non-significant changes in [Formula: see text], suggesting that migration may be at equilibrium between genetically similar sources. Second, we hypothesized that meta-[Formula: see text] might be significantly smaller than the sum of local [Formula: see text]s (null model) if gene flow is asymmetric, varies among subpopulations, and is driven by common landscape features such as waterfalls. One 'bottom-up' or numerical approach that explicitly incorporates variable and asymmetric migration rates showed this very pattern, while a number of analytical models provided meta-[Formula: see text] estimates that were not significantly different from the null model or from each other. Our study of three species inhabiting a shared environment highlights the importance and utility of differentiating species-specific and landscape effects, not only on dispersal but also in the demography of wild populations as assessed through local [Formula: see text]s and meta-[Formula: see text]s and their relevance in ecology, evolution and conservation.

Trophic interference by Salmo trutta on Aplochiton zebra and Aplochiton taeniatus in southern Patagonian lakes

The length and mass ratio, diet and isotopic composition of Aplochiton zebra and Aplochiton taeniatus inhabiting a Salmo trutta-invaded and a S. trutta-free lake in southern Patagonia were compared. Results indicate that S. trutta exercises important trophic interference over A. zebra and A. taeniatus, causing changes in their dietary composition by reducing the consumption of winged Diptera through changes in feeding behaviours that involve jumping out of the water. This effect is significantly higher in A. zebra than in A. taeniatus a species that has a highly specialized diet. The dietary changes of A. zebra and A. taeniatus in sympatry with S. trutta lead to an impoverishment of their isotopic nitrogen signals (δ(15)N), suggesting a reduction of their trophic position. In the case of A. zebra, this translates into a significant decrease in its body condition factor. Such interference could lead to a population decline of this species and would explain the current distribution range decline and allopatry with S. trutta in fluvial systems.

Genomic islands of divergence and their consequences for the resolution of spatial structure in an exploited marine fish

As populations diverge, genomic regions associated with adaptation display elevated differentiation. These genomic islands of adaptive divergence can inform conservation efforts in exploited species, by refining the delineation of management units, and providing genomic tools for more precise and effective population monitoring and the successful assignment of individuals and products. We explored heterogeneity in genomic divergence and its impact on the resolution of spatial population structure in exploited populations of Atlantic cod, Gadus morhua, using genome wide expressed sequence derived single nucleotide polymorphisms in 466 individuals sampled across the range. Outlier tests identified elevated divergence at 5.2% of SNPs, consistent with directional selection in one-third of linkage groups. Genomic regions of elevated divergence ranged in size from a single position to several cM. Structuring at neutral loci was associated with geographic features, whereas outlier SNPs revealed genetic discontinuities in both the eastern and western Atlantic. This fine-scale geographic differentiation enhanced assignment to region of origin, and through the identification of adaptive diversity, fundamentally changes how these populations should be conserved. This work demonstrates the utility of genome scans for adaptive divergence in the delineation of stock structure, the traceability of individuals and products, and ultimately a role for population genomics in fisheries conservation.

Mature male parr contribution to the effective size of an anadromous Atlantic salmon (Salmo salar) population over 30 years

We describe temporal changes in the genetic composition of a small anadromous Atlantic salmon (Salmo salar) population from South Newfoundland, an area where salmon populations are considered threatened (COSEWIC 2010). We examined the genetic variability (13 microsatellite loci) in 869 out-migrating smolt and post-spawning kelt samples, collected from 1985 to 2011 for a total of 22 annual collections and a 30 year span of assigned cohorts. We estimated the annual effective number of breeders (Nb) and the generational effective population size (Ne) through genetic methods and demographically using the adult sex ratio. Comparisons between genetic and demographic estimates show that the adult spawners inadequately explain the observed Ne estimates, suggesting that mature male parr are significantly increasing Nb and Ne over the study period. Spawning as parr appears to be a viable and important strategy in the near absence of adult males.

Estuarine survival and migratory behaviour of Atlantic salmon Salmo salar smolts

To estimate mortality rates, assess the spatio-temporal dynamics of natural mortality and examine migratory behaviour during the fresh to saltwater transition, 185 wild Atlantic salmon Salmo salar smolts were implanted with coded acoustic transmitters. Seaward migration of tagged S. salar from four river systems in an area of Nova Scotia, Canada known as the Southern Upland was monitored using fixed receivers and active telemetry over 3 years. Cumulative survival through the river, inner estuary, outer estuary and bay habitats averaged 59·6% (range = 39·4-73·5%). When standardized to distance travelled, survival rates followed two patterns: (1) constant rates of survival independent of habitat or (2) low survival most frequently associated with inner estuary habitats. In rivers where survival was independent of habitat, residency periods were also independent of habitat, post-smolts exhibited few upstream movements, took a more direct route to the ocean and reached the ocean rapidly. Alternatively, in rivers where survival was habitat specific, residency was also habitat specific with overall increased residency, more frequent upstream movements and delayed arrival to the open ocean. The sudden disappearance of most (75-100%) smolts and post-smolts assumed dead during the course of this study warrants further examination into the role of avian predators as a mortality vector.

Permanent genetic resources added to Molecular Ecology Resources Database 1 December 2011-31 January 2012

This article documents the addition of 473 microsatellite marker loci and 71 pairs of single-nucleotide polymorphism (SNP) sequencing primers to the Molecular Ecology Resources Database. Loci were developed for the following species: Barteria fistulosa, Bombus morio, Galaxias platei, Hematodinium perezi, Macrocentrus cingulum Brischke (a.k.a. M. abdominalis Fab., M. grandii Goidanich or M. gifuensis Ashmead), Micropogonias furnieri, Nerita melanotragus, Nilaparvata lugens Stål, Sciaenops ocellatus, Scomber scombrus, Spodoptera frugiperda and Turdus lherminieri. These loci were cross-tested on the following species: Barteria dewevrei, Barteria nigritana, Barteria solida, Cynoscion acoupa, Cynoscion jamaicensis, Cynoscion leiarchus, Cynoscion nebulosus, Cynoscion striatus, Cynoscion virescens, Macrodon ancylodon, Menticirrhus americanus, Nilaparvata muiri and Umbrina canosai. This article also documents the addition of 116 sequencing primer pairs for Dicentrarchus labrax.

Climate-induced changes to the ancestral population size of two Patagonian galaxiids: the influence of glacial cycling

Patagonia is one of the few areas in the Southern Hemisphere to have been directly influenced by Quaternary glaciers. In this study, we evaluate the influence that Quaternary glacial ice had on the genetic diversity of two congeneric fish species, the diadromous Galaxias maculatus and the nondiadromous Galaxias platei, using multilocus estimates of effective population size through time. Mid-Quaternary glaciations had far-reaching consequences for both species. Galaxias maculatus and G. platei each experienced severe genetic bottlenecks during the period when Patagonia ice sheet advance reached its maximum positions c. 1.1-0.6 Ma. Concordant drops in effective size during this time suggest that range sizes were under similar constraints. It is therefore unlikely that coastal (brackish/marine) environments served as a significant refuge for G. maculatus during glacial periods. An earlier onset of population declines for G. platei suggests that this species was vulnerable to modest glacial advances. Declines in effective sizes were continuous for both species and lasted into the late-Pleistocene. However, G. maculatus exhibited a strong population recovery during the late-Quaternary (c. 400,000 bp). Unusually long and warm interglacials associated with the late-Quaternary may have helped to facilitate a strong population rebound in this primarily coastal species.