Search

In genomics-scale datasets, loci are closely packed within chromosomes and hence provide correlated information. Averaging across loci as if they were independent creates pseudoreplication, which reduces the effective degrees of freedom (df’) compared to the nominal degrees of freedom, df. This issue has been known for some time, but consequences have not been systematically quantified across the entire genome. Here we measured pseudoreplication (quantified by the ratio df’/df) for a common metric of genetic differentiation (FST) and a common measure of linkage disequilibrium between pairs of loci (r2). Based on data simulated using models (SLiM and msprime) that allow efficient forward-in-time and coalescent simulations while precisely controlling population pedigrees, we estimated df’ and df’/df by measuring the rate of decline in the variance of mean FST and mean r2 as more loci were used. For both indices, df’ increases with Ne and genome size, as expected. However, even for large Ne and large genomes, df’ for mean r2 plateaus after a few thousand loci, and a variance components analysis indicates that the limiting factor is uncertainty associated with sampling individuals rather than genes. Pseudoreplication is less extreme for FST, but df’/df ≤0.01 can occur in datasets using tens of thousands of loci. Commonly-used block-jackknife methods consistently overestimated var(FST), producing very conservative confidence intervals. Predicting df’ based on our modeling results as a function of Ne, L, S, and genome size provides a robust way to quantify precision associated with genomics-scale datasets.

Knowledge of how individuals are related is important in many areas of research, and numerous methods for inferring pairwise relatedness from genetic data have been developed. However, the majority of these methods were not developed for situations where data are limited. Specifically, most methods rely on the availability of population allele frequencies, the relative genomic position of variants and accurate genotype data. But in studies of non‐model organisms or ancient samples, such data are not always available. Motivated by this, we present a new method for pairwise relatedness inference, which requires neither allele frequency information nor information on genomic position. Furthermore, it can be applied not only to accurate genotype data but also to low‐depth sequencing data from which genotypes cannot be accurately called. We evaluate it using data from a range of human populations and show that it can be used to infer close familial relationships with a similar accuracy as a widely used method that relies on population allele frequencies. Additionally, we show that our method is robust to SNP ascertainment and applicable to low‐depth sequencing data generated using different strategies, including resequencing and RADseq, which is important for application to a diverse range of populations and species.

Whole genome duplications have occurred in the recent ancestors of many plants, fish, and amphibians, resulting in a pervasiveness of paralogous loci and the potential for both disomic and tetrasomic inheritance in the same genome. Paralogs can be difficult to reliably genotype and are often excluded from genotyping-by-sequencing (GBS) analyses; however, removal requires paralogs to be identified which is difficult without a reference genome. We present a method for identifying paralogs in natural populations by combining two properties of duplicated loci: 1) the expected frequency of heterozygotes exceeds that for singleton loci, and 2) within heterozygotes, observed read ratios for each allele in GBS data will deviate from the 1:1 expected for singleton (diploid) loci. These deviations are often not apparent within individuals, particularly when sequence coverage is low; but, we postulated that summing allele reads for each locus over all heterozygous individuals in a population would provide sufficient power to detect deviations at those loci. We identified paralogous loci in three species: Chinook salmon (Oncorhynchus tshawytscha) which retains regions with ongoing residual tetrasomy on eight chromosome arms following a recent whole genome duplication, mountain barberry (Berberis alpina) which has a large proportion of paralogs that arose through an unknown mechanism, and dusky parrotfish (Scarus niger) which has largely re-diploidized following an ancient whole genome duplication. Importantly, this approach only requires the genotype and allele-specific read counts for each individual, information which is readily obtained from most GBS analysis pipelines. This article is protected by copyright. All rights reserved.

Effective population size (Ne) is among the most important metrics in evolutionary biology. In natural populations, it is often difficult to collect adequate demographic data to calculate Ne directly. Consequently, genetic methods to estimate Ne have been developed. Two Ne estimators based on sibship reconstruction using multi-locus genotype data have been developed in recent years; sibship assignment and parentage analysis without parents. In this study, we evaluated the accuracy of sibship reconstruction using a large empirical dataset from five hatchery steelhead populations with known pedigrees and using 95 single nucleotide polymorphism (SNP) markers. We challenged the software COLONY with 2,599,961 known relationships and demonstrated that reconstruction of full-sib and unrelated pairs were greater than 95% and 99% accurate, respectively. However, reconstruction of half-sib pairs was poor (< 5% accurate). Despite poor half-sib reconstruction, both estimators provided accurate estimates of the effective number of breeders (Nb) when sample sizes were near or greater than the true Nb and when assuming a monogamous mating system. We further demonstrated that both methods provide roughly equivalent estimates of Nb. Our results indicate that sibship reconstruction and current SNP panels provide promise for estimating Nb in steelhead populations in the region. This article is protected by copyright. All rights reserved.

Coral reefs within the Coral Triangle (CT) are home to the greatest marine diversity on the globe and are an important supplier of marine resources to densely populated coastal regions. Many coral reefs within the CT and around the world are under threat from over-exploitation. Marine protected areas (MPAs) have been proven to be effective tools in restoring fish stocks. However, the role that MPAs play in promoting connectivity at greater distances through larval dispersal is still unknown. RADseq was used to discover 4253 SNPs from 81 individuals of the dusky parrotfish (Scarus niger) collected from three sites within the Philippines. A lack of population structure suggested a high rate of gene flow (F ST = 0.007). Estimates of Ne from linkage disequilibrium are relatively large, ranging from 1200 to 2000. A sibling analysis revealed one pair of well-supported full siblings (r = 0.773) and one pair of putative half siblings (r = 0.191) between sites separated by more than 500 km. The low F ST values indicate a high degree of gene flow between the reefs within the sampling area while the sibling analysis suggests demographic connectivity between the Sibuyan Sea and the Sulu Sea. The Mindoro–Panay throughflow is a likely vector by which larvae are transported between these sites, suggesting that reefs in Romblon are sources for reefs near Basay, 400 km to the south. Given the reliance of a vast majority of coral reef fishes on larval dispersal, this study reveals that MPAs established within the central Philippines can supply varying levels of larvae to overfished reefs.

Whole-genome duplications have occurred in the recent ancestors of many plants, fish and amphibians. Signals of these whole-genome duplications still exist in the form of paralogous loci. Recent advances have allowed reliable identification of paralogs in genotyping-by-sequencing (GBS) data such as that generated from restriction-site-associated DNA sequencing (RADSeq); however, excluding paralogs from analyses is still routine due to difficulties in genotyping. This exclusion of paralogs may filter a large fraction of loci, including loci that may be adaptively important or informative for population genetic analyses. We present a maximum-likelihood method for inferring allele dosage in paralogs and assess its accuracy using simulated GBS, empirical RADSeq and amplicon sequencing data from Chinook salmon. We accurately infer allele dosage for some paralogs from a RADSeq data set and show how accuracy is dependent upon both read depth and allele frequency. The amplicon sequencing data set, using RADSeq-derived markers, achieved sufficient depth to infer allele dosage for all paralogs. This study demonstrates that RADSeq locus discovery combined with amplicon sequencing of targeted loci is an effective method for incorporating paralogs into population genetic analyses.

Whole-genome duplications are major evolutionary events with a lasting impact on genome structure. Duplication events complicate genetic analyses as paralogous sequences are difficult to distinguish; consequently, paralogs are often excluded from studies. The effects of an ancient whole-genome duplication (approximately 88 MYA) are still evident in salmonids through the persistence of numerous paralogous gene sequences and partial tetrasomic inheritance. We use restriction site-associated DNA sequencing on 10 collections of chum salmon from the Salish Sea in the USA and Canada to investigate genetic diversity and population structure in both tetrasomic and rediploidized regions of the genome. We use a pedigree and high-density linkage map to identify paralogous loci and to investigate genetic variation across the genome. By applying multivariate statistical methods, we show that it is possible to characterize paralogous loci and that they display similar patterns of population structure as the diploidized portion of the genome. We find genetic associations with the adaptively important trait of run-timing in both sets of loci. By including paralogous loci in genome scans, we can observe evolutionary signals in genomic regions that have routinely been excluded from population genetic studies in other polyploid-derived species.

Salmonids are an important cultural and ecological resource exhibiting near worldwide distribution between their native and introduced range. Previous research has generated linkage maps and genomic resources for several species as well as genome assemblies for two species. We first leveraged improvements in mapping and genotyping methods to create a dense linkage map for Chinook salmon Oncorhynchus tshawytscha by assembling family data from different sources. We successfully mapped 14 620 SNP loci including 2336 paralogs in subtelomeric regions. This improved map was then used as a foundation to integrate genomic resources for gene annotation and population genomic analyses. We anchored a total of 286 scaffolds from the Atlantic salmon genome to the linkage map to provide a framework for the placement 11 728 Chinook salmon ESTs. Previously identified thermotolerance QTL were found to colocalize with several candidate genes including HSP70, a gene known to be involved in thermal response, as well as its inhibitor. Multiple regions of the genome with elevated divergence between populations were also identified, and annotation of ESTs in these regions identified candidate genes for fitness related traits such as stress response, growth and behaviour. Collectively, these results demonstrate the utility of combining genomic resources with linkage maps to enhance evolutionary inferences.

The future spread and impact of an introduced species will depend on how it adapts to the abiotic and biotic conditions encountered in its new range, so the potential for rapid evolution subsequent to species introduction is a critical, evolutionary dimension of invasion biology. Using a resurrection approach, we provide a direct test for change over time within populations in a species' introduced range, in the Asian shade annual Polygonum cespitosum. We document, over an 11-year period, the evolution of increased reproductive output as well as greater physiological and root-allocational plasticity in response to the more open, sunny conditions found in the North American range in which the species has become invasive. These findings show that extremely rapid adaptive modifications to ecologically-important traits and plastic expression patterns can evolve subsequent to a species' introduction, within populations established in its introduced range. This study is one of the first to directly document evolutionary change in adaptive plasticity. Such rapid evolutionary changes can facilitate the spread of introduced species into novel habitats and hence contribute to their invasive success in a new range. The data also reveal how evolutionary trajectories can differ among populations in ways that can influence invasion dynamics.

An important use of genetic parentage analysis is the ability to directly calculate the number of offspring produced by each parent (ki) and hence effective population size, Ne. But what if parental genotypes are not available? In theory, given enough markers, it should be possible to reconstruct parental genotypes based entirely on a sample of progeny, and if so the vector of parental ki values. However, this would provide information only about parents that actually contributed offspring to the sample. How would ignoring the ‘null’ parents (those that produced no offspring) affect an estimate of Ne? The surprising answer is that null parents have no effect at all. We show that: (i) The standard formula for inbreeding Ne can be rewritten so that it is a function only of sample size and P k 2 �� ; it is not necessary to know the total number of parents (N). This same relationship does not hold for variance Ne. (ii) This novel formula provides an unbiased estimate of Ne even if only a subset of progeny is available, provided the parental contributions are accurately determined, in which case precision is also high compared to other single-sample estimators of Ne. (iii) It is not necessary to actually reconstruct parental genotypes; from a matrix of pairwise relationships (as can be estimated by some current software programs), it is possible to construct the vector of ki values and estimate Ne. The new method based on parentage analysis without parents (PwoP) can potentially be useful as a single-sample estimator of contemporary Ne, provided that either (i) relationships can be accurately determined, or (ii) P k 2 �� can be estimated directly.

Studying the effect of similar environments on diverse genetic backgrounds has long been a goal of evolutionary biologists with studies typically relying on experimental approaches. Pink salmon, a highly abundant and widely ranging salmonid, provide a naturally occurring opportunity to study the effects of similar environments on divergent genetic backgrounds due to a strict two-year semelparous life history. The species is composed of two reproductively isolated lineages with overlapping ranges that share the same spawning and rearing environments in alternate years. We used restriction-site-associated DNA (RAD) sequencing to discover and genotype approximately 8000 SNP loci in three population pairs of even- and odd-year pink salmon along a latitudinal gradient in North America. We found greater differentiation within the odd-year than within the even-year lineage and greater differentiation in the southern pair from Puget Sound than in the northern Alaskan population pairs. We identified 15 SNPs reflecting signatures of parallel selection using both a differentiation-based method (BAYESCAN) and an environmental correlation method (BAYENV). These SNPs represent genomic regions that may be particularly informative in understanding adaptive evolution in pink salmon and exploring how differing genetic backgrounds within a species respond to selection from the same natural environment.

Recent advances in population genomics have made it possible to detect previously unidentified structure, obtain more accurate estimates of demographic parameters, and explore adaptive divergence, potentially revolutionizing the way genetic data are used to manage wild populations. Here, we identified 10 944 single-nucleotide polymorphisms using restriction-site-associated DNA (RAD) sequencing to explore population structure, demography, and adaptive divergence in five populations of Chinook salmon (Oncorhynchus tshawytscha) from western Alaska. Patterns of population structure were similar to those of past studies, but our ability to assign individuals back to their region of origin was greatly improved (>90% accuracy for all populations). We also calculated effective size with and without removing physically linked loci identified from a linkage map, a novel method for nonmodel organisms. Estimates of effective size were generally above 1000 and were biased downward when physically linked loci were not removed. Outlier tests based on genetic differentiation identified 733 loci and three genomic regions under putative selection. These markers and genomic regions are excellent candidates for future research and can be used to create high-resolution panels for genetic monitoring and population assignment. This work demonstrates the utility of genomic data to inform conservation in highly exploited species with shallow population structure.