Genomics reveal population structure, evolutionary history, and signatures of selection in the northern bottlenose whale, Hyperoodon ampullatus.

Information on wildlife population structure, demographic history, and adaptations are fundamental to understanding species evolution and informing conservation strategies. To study this ecological context for a cetacean of conservation concern, we conducted the first genomic assessment of the northern bottlenose whale, Hyperoodon ampullatus, using whole-genome resequencing data (n = 37) from five regions across the North Atlantic Ocean. We found a range-wide pattern of isolation-by-distance with a genetic subdivision distinguishing three subgroups: the Scotian Shelf, western North Atlantic, and Jan Mayen regions. Signals of elevated levels of inbreeding in the Endangered Scotian Shelf population indicate this population may be more vulnerable than the other two subgroups. In addition to signatures of inbreeding, evidence of local adaptation in the Scotian Shelf was detected across the genome. We found a long-term decline in effective population size for the species, which poses risks to their genetic diversity and may be exacerbated by the isolating effects of population subdivision. Protecting important habitat and migratory corridors should be prioritized to rebuild population sizes that were diminished by commercial whaling, strengthen gene flow, and ensure animals can move across regions in response to environmental changes.

Resistance and resilience of genetic and phenotypic diversity to "black swan" flood events: A retrospective analysis with historical samples of guppies.

Rare extreme "black swan" disturbances can impact ecosystems in many ways, such as destroying habitats, depleting resources, and causing high mortality. In rivers, for instance, exceptional floods that occur infrequently (e.g., so-called "50-year floods") can strongly impact the abundance of fishes and other aquatic organisms. Beyond such ecological effects, these floods could also impact intraspecific diversity by elevating genetic drift or dispersal and by imposing strong selection, which could then influence the population's ability to recover from disturbance. And yet, natural systems might be resistant (show little change) or resilient (show rapid recovery) even to rare extreme events - perhaps as a result of selection due to past events. We considered these possibilities in two rivers where native guppies experienced two extreme floods - one in 2005 and another in 2016. For each river, we selected four sites and used archived "historical" samples to compare levels of genetic and phenotypic diversity before vs. after floods. Genetic diversity was represented by 33 neutral microsatellite markers, and phenotypic diversity was represented by body length and male melanic (black) colour. We found that genetic diversity and population structure was mostly "resistant" to even these extreme floods; whereas the larger impacts on phenotypic diversity were short-lived, suggesting additional "resilience". We discuss the determinants of these two outcomes for guppies facing floods, and then consider the general implications for the resistance and resilience of intraspecific variation to black swan disturbances.

Multiple decades of stocking has resulted in limited hatchery introgression in wild brook trout (Salvelinus fontinalis) populations of Nova Scotia.

Many populations of freshwater fishes are threatened with losses, and increasingly, the release of hatchery individuals is one strategy being implemented to support wild populations. However, stocking of hatchery individuals may pose long-term threats to wild populations, particularly if genetic interactions occur between wild and hatchery individuals. One highly prized sport fish that has been heavily stocked throughout its range is the brook trout (Salvelinus fontinalis). In Nova Scotia, Canada, hatchery brook trout have been stocked since the early 1900s, and despite continued stocking efforts, populations have suffered declines in recent decades. Before this study, the genetic structure of brook trout populations in the province was unknown; however, given the potential negative consequences associated with hatchery stocking, it is possible that hatchery programs have adversely affected the genetic integrity of wild populations. To assess the influence of hatchery supplementation on wild populations, we genotyped wild brook trout from 12 river systems and hatchery brook trout from two major hatcheries using 100 microsatellite loci. Genetic analyses of wild trout revealed extensive population genetic structure among and within river systems and significant isolation-by-distance. Hatchery stocks were genetically distinct from wild populations, and most populations showed limited to no evidence of hatchery introgression (<5% hatchery ancestry). Only a single location had a substantial number of hatchery-derived trout and was located in the only river where a local strain is used for supplementation. The amount of hatchery stocking within a watershed did not influence the level of hatchery introgression. Neutral genetic structure of wild populations was influenced by geography with some influence of climate and stocking indices. Overall, our study suggests that long-term stocking has not significantly affected the genetic integrity of wild trout populations, highlighting the variable outcomes of stocking and the need to evaluate the consequences on a case-by-case basis.

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.

Evolutionary impacts differ between two exploited populations of northern bottlenose whale (Hyperoodon ampullatus).

Interpretation of conservation status should be informed by an appreciation of genetic diversity, past demography, and overall trends in population size, which contribute to a species' evolutionary potential and resilience to genetic risks. Low genetic diversity can be symptomatic of rapid demographic declines and impose genetic risks to populations, but can also be maintained by natural processes. The northern bottlenose whale Hyperoodon ampullatus has the lowest known mitochondrial diversity of any cetacean and was intensely whaled in the Northwest Atlantic over the last century, but whether exploitation imposed genetic risks that could limit recovery is unknown. We sequenced full mitogenomes and genotyped 37 novel microsatellites for 128 individuals from known areas of abundance in the Scotian Shelf, Northern and Southern Labrador, Davis Strait, and Iceland, and a newly discovered group off Newfoundland. Despite low diversity and shared haplotypes across all regions, both markers supported the Endangered Scotian Shelf population as distinct from the combined northern regions. The genetic affinity of Newfoundland was uncertain, suggesting an area of mixing with no clear population distinction for the region. Demographic reconstruction using mitogenomes suggests that the northern region underwent population expansion following the last glacial maximum, but for the peripheral Scotian Shelf population, a stable demographic trend was followed by a drastic decline over a temporal scale consistent with increasing human activity in the Northwest Atlantic. Low connectivity between the Scotian Shelf and the rest of the Atlantic likely compounded the impact of intensive whaling for this species, potentially imposing genetic risks affecting recovery of this population. We highlight how the combination of historical environmental conditions and modern exploitation of this species has had very different evolutionary impacts on structured populations of northern bottlenose whales across the western North Atlantic.

megasat: automated inference of microsatellite genotypes from sequence data.

megasat is software that enables genotyping of microsatellite loci using next-generation sequencing data. Microsatellites are amplified in large multiplexes, and then sequenced in pooled amplicons. megasat reads sequence files and automatically scores microsatellite genotypes. It uses fuzzy matches to allow for sequencing errors and applies decision rules to account for amplification artefacts, including nontarget amplification products, replication slippage during PCR (amplification stutter) and differential amplification of alleles. An important feature of megasat is the generation of histograms of the length-frequency distributions of amplification products for each locus and each individual. These histograms, analogous to electropherograms traditionally used to score microsatellite genotypes, enable rapid evaluation and editing of automatically scored genotypes. megasat is written in Perl, runs on Windows, Mac OS X and Linux systems, and includes a simple graphical user interface. We demonstrate megasat using data from guppy, Poecilia reticulata. We genotype 1024 guppies at 43 microsatellites per run on an Illumina MiSeq sequencer. We evaluated the accuracy of automatically called genotypes using two methods, based on pedigree and repeat genotyping data, and obtained estimates of mean genotyping error rates of 0.021 and 0.012. In both estimates, three loci accounted for a disproportionate fraction of genotyping errors; conversely, 26 loci were scored with 0-1 detected error (error rate ≤0.007). Our results show that with appropriate selection of loci, automated genotyping of microsatellite loci can be achieved with very high throughput, low genotyping error and very low genotyping costs.

Complete mitochondrial genomes for Icelus spatula, Aspidophoroides olrikii and Leptoclinus maculatus: pan-Arctic marine fishes from Canadian waters.

Three Arctic marine fishes Icelus spatula, Aspidophoroides olrikii and Leptoclinus maculatus have been identified as target species for investigating the effects of ocean warming on population patterns in high-latitude marine habitats around Canada. In preparation for this research, we have resolved whole mitochondrial genome sequences of 16 384, 17 200 and 16 384 bp for each species, respectively. GC content for each species was 47.5%, 44.2% and 45.3%, respectively. Mitogenome gene composition included 13 protein-encoding genes, 2 rRNA and 22 tRNA genes, for I. spatula and L. maculatus, consistent with other teleosts. Only 20 tRNA genes were annotated for A. olrikii, because tRNA-Pro and tRNA-Thr are poorly characterized and aberrantly located in this species.

Ultra-deep Illumina sequencing accurately identifies MHC class IIb alleles and provides evidence for copy number variation in the guppy (Poecilia reticulata).

We address the bioinformatic issue of accurately separating amplified genes of the major histocompatibility complex (MHC) from artefacts generated during high-throughput sequencing workflows. We fit observed ultra-deep sequencing depths (hundreds to thousands of sequences per amplicon) of allelic variants to expectations from genetic models of copy number variation (CNV). We provide a simple, accurate and repeatable method for genotyping multigene families, evaluating our method via analyses of 209 b of MHC class IIb exon 2 in guppies (Poecilia reticulata). Genotype repeatability for resequenced individuals (N = 49) was high (100%) within the same sequencing run. However, repeatability dropped to 83.7% between independent runs, either because of lower mean amplicon sequencing depth in the initial run or random PCR effects. This highlights the importance of fully independent replicates. Significant improvements in genotyping accuracy were made by greatly reducing type I genotyping error (i.e. accepting an artefact as a true allele), which may occur when using low-depth allele validation thresholds used by previous methods. Only a small amount (4.9%) of type II error (i.e. rejecting a genuine allele as an artefact) was detected through fully independent sequencing runs. We observed 1-6 alleles per individual, and evidence of sharing of alleles across loci. Variation in the total number of MHC class II loci among individuals, both among and within populations was also observed, and some genotypes appeared to be partially hemizygous; total allelic dosage added up to an odd number of allelic copies. Collectively, observations provide evidence of MHC CNV and its complex basis in natural populations.

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.

Beyond lifetime reproductive success: the posthumous reproductive dynamics of male Trinidadian guppies.

In semelparous populations, dormant germ banks (e.g. seeds) have been proposed as important in maintaining genotypes that are adaptive at different times in fluctuating environments. Such hidden storage of genetic diversity need not be exclusive to dormant banks. Genotype diversity may be preserved in many iteroparous animals through sperm-storage mechanisms in females. This allows males to reproduce posthumously and increase the effective sizes of seemingly female-biased populations. Although long-term sperm storage has been demonstrated in many organisms, the understanding of its importance in the wild is very poor. We here show the prevalence of male posthumous reproduction in wild Trinidadian guppies, through the combination of mark-recapture and pedigree analyses of a multigenerational individual-based dataset. A significant proportion of the reproductive population consisted of dead males, who could conceive up to 10 months after death (the maximum allowed by the length of the dataset), which is more than twice the estimated generation time. Demographic analysis shows that the fecundity of dead males can play an important role in population growth and selection.

Evaluating SNP ascertainment bias and its impact on population assignment in Atlantic cod, Gadus morhua.

The increasing use of single nucleotide polymorphisms (SNPs) in studies of nonmodel organisms accentuates the need to evaluate the influence of ascertainment bias on accurate ecological or evolutionary inference. Using a panel of 1641 expressed sequence tag-derived SNPs developed for northwest Atlantic cod (Gadus morhua), we examined the influence of ascertainment bias and its potential impact on assignment of individuals to populations ranging widely in origin. We hypothesized that reductions in assignment success would be associated with lower diversity in geographical regions outside the location of ascertainment. Individuals were genotyped from 13 locations spanning much of the contemporary range of Atlantic cod. Diversity, measured as average sample heterozygosity and number of polymorphic loci, declined (c. 30%) from the western (H(e) = 0.36) to eastern (H(e) = 0.25) Atlantic, consistent with a signal of ascertainment bias. Assignment success was examined separately for pools of loci representing differing degrees of reductions in diversity. SNPs displaying the largest declines in diversity produced the most accurate assignment in the ascertainment region (c. 83%) and the lowest levels of correct assignment outside the ascertainment region (c. 31%). Interestingly, several isolated locations showed no effect of assignment bias and consistently displayed 100% correct assignment. Contrary to expectations, estimates of accurate assignment range-wide using all loci displayed remarkable similarity despite reductions in diversity. Our results support the use of large SNP panels in assignment studies of high geneflow marine species. However, our evidence of significant reductions in assignment success using some pools of loci suggests that ascertainment bias may influence assignment results and should be evaluated in large-scale assignment studies.

Isolation and differentiation of Rivulus hartii across Trinidad and neighboring islands.

Diversification of freshwater fishes on islands is considered unlikely because the traits that enable successful colonization-specifically, broad salinity tolerances and the potential for oceanic dispersal-may also constrain post-colonization genetic differentiation. Some secondary freshwater fish, however, exhibit pronounced genetic differentiation and geographic structure on islands, whereas others do not. It is unclear what conditions give rise to contrasting patterns of differentiation because few comparative reconstructions of population history have been carried out for insular freshwater fishes. In this study, we examined the phylogeography of Hart's killifish (Rivulus hartii) across Trinidad, with reference to neighboring islands and northern South America, to test hypotheses of colonization and differentiation derived from comparable work on co-occurring guppies (Poecilia reticulata). Geographic patterns of mitochondrial DNA haplotype variation and microsatellite genotype variation provide evidence of genetic differentiation of R. hartii among islands and across Trinidad. Our findings are largely consistent with patterns of geographically structured ancestry and admixture found in Trinidadian guppies, which suggests that both species share a history of colonization and differentiation and that post-colonization diversification may be more common among members of insular freshwater fish assemblages than has been previously thought.

Parallel adaptive evolution of Atlantic cod on both sides of the Atlantic Ocean in response to temperature.

Despite the enormous economic and ecological importance of marine organisms, the spatial scales of adaptation and biocomplexity remain largely unknown. Yet, the preservation of local stocks that possess adaptive diversity is critical to the long-term maintenance of productive stable fisheries and ecosystems. Here, we document genomic evidence of range-wide adaptive differentiation in a broadcast spawning marine fish, Atlantic cod (Gadus morhua), using a genome survey of single nucleotide polymorphisms. Of 1641 gene-associated polymorphisms examined, 70 (4.2%) tested positive for signatures of selection using a Bayesian approach. We identify a subset of these loci (n=40) for which allele frequencies show parallel temperature-associated clines (p<0.001, r2=0.89) in the eastern and western north Atlantic. Temperature associations were robust to the statistical removal of geographic distance or latitude effects, and contrasted 'neutral' loci, which displayed no temperature association. Allele frequencies at temperature-associated loci were significantly correlated, spanned three linkage groups and several were successfully annotated supporting the involvement of multiple independent genes. Our results are consistent with the evolution and/or selective sweep of multiple genes in response to ocean temperature, and support the possibility of a new conservation paradigm for non-model marine organisms based on genomic approaches to resolving functional and adaptive diversity.

Ten novel microsatellite loci characterized for a remarkably widespread fish: Galaxias maculatus (Galaxiidae).

Ten polymorphic microsatellite markers (five tetra-, one compound tetra-, one octa- and three dinucleotides) were isolated and characterized for Galaxias maculatus, a fish species widely distributed in the Southern Hemisphere. Markers were tested in 89 individual samples from a single location and the number of alleles ranged between 2 and 28. Observed and expected heterozygosities ranged from 0.103 to 0.910 and 0.098 to 0.935 respectively. No evidence was detected for either linkage disequilibrium (P-values > 0.05 for each locus pair) or deviations from HWE (P-values > 0.05 for every loci).

Isolation and characterization of 13 microsatellite loci for Percichthys trucha (Percichthyidae).

Thirteen polymorphic microsatellite loci are described for the South American freshwater fish Percichthys trucha. Number of alleles per locus ranged from two to 21 and observed heterozygosities ranged from 0.304 to 0.915 in a sample of 47 individuals from four different sampling locations.