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 reaching N ̂ e $$ {\hat{N}}_e $$ = 50 and a <9% probability of reaching N ̂ e $$ {\hat{N}}_e $$ = 500. Populations listed as being of conservation concern according to the IUCN Red List had a smaller median N ̂ 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 of N ̂ 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.

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 87 k 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.

La colonización recurrente de hábitats de agua dulce y la subsecuente pérdida de diadromía es una transición ecológica importante que ha sido reportada en varias especies de peces con ancestros diádromos. Esta residencia está acompañada frecuentemente por la pérdida de tolerancia a ambientes de agua salada. Galaxias maculatus, especie anfídroma, ha colonizado ríos repetidamente y existe evidencia que las poblaciones residentes presentan diferencias respecto a la tolerancia al agua salada. En este estudio, usamos transcriptómica para dilucidar los mecanismos que contribuyen a la reducida tolerancia a altas salinidades en las poblaciones residentes de agua dulce. Realizamos un desafío agudo de salinidad (0 ppt a 23-2 ppt) y medimos la habilidad osmoreguladora (contenido de agua en músculo) por 48 horas en individuos de tres poblaciones: una diádroma, una intolerante a agua salada (Toltén) y una tolerante a agua salada (Valdivia). Con el secuenciamiento de ARN de las branquias identificamos los genes expresados diferencialmente al cambio de salinidad y cuales están asociados a la pérdida de tolerancia a agua salada en la población de Toltén. Genes claves asociados a la aclimatación al agua salada fueron caracterizados en individuos diádromos, algunos de ellos también se expresaron en la población residente tolerante al agua salada (Valdivia). Sin embargo, algunos genes involucrados en la aclimatación al agua salada, incluyendo el gen regulador de la conductancia transmembrana de la fibrosis quística (CFTR), no se diferenciaron significativamente en la población residente intolerante al agua salada (Toltén), sugiriendo un mecanismo potencial de la pérdida de tolerancia a ambientes con salinidad elevada. Como el conjunto de genes expresados difiere entre las dos poblaciones residentes al compararse con la población diádroma, hipotetizamos que la pérdida de diadromía resulta en trayectorias evolutivas únicas debido a deriva génica, por lo que la pérdida de la diadromía no necesariamente conlleva a la pérdida de la tolerancia a aguas saladas.

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.

Demographic resilience of brook trout populations subjected to experimental size-selective harvesting.

Sustainable management of exploited populations benefits from integrating demographic and genetic considerations into assessments, as both play a role in determining harvest yields and population persistence. This is especially important in populations subject to size-selective harvest, because size selective harvesting has the potential to result in significant demographic, life-history, and genetic changes. We investigated harvest-induced changes in the effective number of breeders ( N ^ b ) for introduced brook trout populations (Salvelinus fontinalis) in alpine lakes from western Canada. Three populations were subject to 3 years of size-selective harvesting, while three control populations experienced no harvest. The N ^ c decreased consistently across all harvested populations (on average 60.8%) but fluctuated in control populations. There were no consistent changes in N ^ b between control or harvest populations, but one harvest population experienced a decrease in N ^ b of 63.2%. The N ^ b / N ^ c ratio increased consistently across harvest lakes; however we found no evidence of genetic compensation (where variance in reproductive success decreases at lower abundance) based on changes in family evenness ( FE ^ ) and the number of full-sibling families ( N ^ fam ). We found no relationship between FE ^ and N ^ c or between N ^ fam / N ^ c and FE ^ . We posit that change in N ^ b was buffered by constraints on breeding habitat prior to harvest, such that the same number of breeding sites were occupied before and after harvest. These results suggest that effective size in harvested populations may be resilient to considerable changes in Nc in the short-term, but it is still important to monitor exploited populations to assess the risk of inbreeding and ensure their long-term survival.

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.

The Genomic Consistency of the Loss of Anadromy in an Arctic Fish (Salvelinus alpinus).

AbstractThe potentially significant genetic consequences associated with the loss of migratory capacity of diadromous fishes that have become landlocked in freshwater are poorly understood. Consistent selective pressures associated with freshwater residency may drive repeated differentiation both between allopatric landlocked and anadromous populations and within landlocked populations (resulting in sympatric morphs). Alternatively, the strong genetic drift anticipated in isolated landlocked populations could hinder consistent adaptation, limiting genetic parallelism. Understanding the degree of genetic parallelism underlying differentiation has implications for both the predictability of evolution and management practices. We employed an 87k single-nucleotide polymorphism (SNP) array to examine the genetic characteristics of landlocked and anadromous Arctic char (Salvelinus alpinus) populations from five drainages within Labrador, Canada. One gene was detected as an outlier between sympatric, size-differentiated morphs in each of two landlocked lakes. While no single locus differentiated all replicate pairs of landlocked and anadromous populations, several SNPs, genes, and paralogs were consistently detected as outliers in at least 70% of these pairwise comparisons. A significant C-score suggested that the amount of shared outlier SNPs across all paired landlocked and anadromous populations was greater than expected by chance. Our results indicate that despite their isolation, selection due to the loss of diadromy may drive consistent genetic responses in landlocked populations.

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 (FST  = 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) and N ^ 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.

Limited genetic parallelism underlies recent, repeated incipient speciation in geographically proximate populations of an Arctic fish (Salvelinus alpinus).

The genetic underpinnings of incipient speciation, including the genomic mechanisms which contribute to morphological and ecological differentiation and reproductive isolation, remain poorly understood. The repeated evolution of consistently, phenotypically distinct morphs of Arctic Charr (Salvelinus alpinus) within the Quaternary period offer an ideal model to study the repeatability of evolution at the genomic level. Sympatric morphs of Arctic Charr are found across this species' circumpolar distribution. However, the specific genetic mechanisms driving this morph differentiation are largely unknown despite the cultural and economic importance of the anadromous morph. We used a newly designed 87k SNP chip to investigate the character and consistency of the genomic differences among sympatric morphs within three recently deglaciated and geographically proximate lakes in Labrador, Canada. We found genetically distinct small and large morph Arctic Charr in all three lakes consistent with resident and anadromous morphs, respectively. A degree of reproductive isolation among sympatric morphs is likely given genome-wide distributions of outlier SNPs and high genome-wide FST s. Across all lakes, outlier SNPs were largely nonoverlapping suggesting a lack of genetic parallelism driving morph differentiation. Alternatively, several genes and paralogous copies of the same gene consistently differentiated morphs across multiple lakes suggesting their importance to the manifestation of morphs. Our results confirm the utility of Arctic Charr as a model for investigating the predictability of evolution and support the importance of both genetic parallelism and nonparallelism to the incipient speciation of Arctic Charr morphs.

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.

The effects of diadromy and its loss on genomic divergence: The case of amphidromous Galaxias maculatus populations.

Understanding the evolutionary mechanisms that affect the genetic divergence between diadromous and resident populations across heterogeneous environments is a challenging task. While diadromy may promote gene flow leading to a lack of genetic differentiation among populations, resident populations tend to be affected by local adaptation and/or plasticity. Studies on these effects on genomic divergence in nonmodel amphidromous species are scarce. Galaxias maculatus, one of the most widespread fish species in the Southern Hemisphere, exhibits two life histories, an ancestral diadromous, specifically, amphidromous form, and a derived freshwater resident form. We examined the genetic diversity and divergence among 20 estuarine and resident populations across the Chilean distribution of G. maculatus and assessed the extent to which selection is involved in the differentiation among resident populations. We obtained nearly 4,400 SNP markers using a RADcap approach for 224 individuals. As expected, collections from estuarine locations typically consist of diadromous individuals. Diadromous populations are highly differentiated from their resident counterparts by both neutral and putative adaptive markers. While diadromous populations exhibit high gene flow and lack site fidelity, resident populations appear to be the product of different colonization events with relatively low genetic diversity and varying levels of gene flow. In particular, the northernmost resident populations were clearly genetically distinct and reproductively isolated from each other suggesting local adaptation. Our study provides insights into the role of life history differences in the maintenance of genetic diversity and the importance of genetic divergence in species evolution.