Data from: Design of a 9K SNP chip for polar bears (Ursus maritimus) from RAD and transcriptome sequencing

Single-nucleotide polymorphisms (SNPs) offer numerous advantages over anonymous markers such as microsatellites, including improved estimation of population parameters, finer-scale resolution of population structure and more precise genomic dissection of quantitative traits. However, many SNPs are needed to equal the resolution of a single microsatellite, and reliable large-scale genotyping of SNPs remains a challenge in nonmodel species. Here, we document the creation of a 9K Illumina Infinium BeadChip for polar bears (Ursus maritimus), which will be used to investigate: (i) the fine-scale population structure among Canadian polar bears and (ii) the genomic architecture of phenotypic traits in the Western Hudson Bay subpopulation. To this end, we used restriction-site associated DNA (RAD) sequencing from 38 bears across their circumpolar range, as well as blood/fat transcriptome sequencing of 10 individuals from Western Hudson Bay. Six-thousand RAD SNPs and 3000 transcriptomic SNPs were selected for the chip, based primarily on genomic spacing and gene function respectively. Of the 9000 SNPs ordered from Illumina, 8042 were successfully printed, and – after genotyping 1450 polar bears – 5441 of these SNPs were found to be well clustered and polymorphic. Using this array, we show rapid linkage disequilibrium decay among polar bears, we demonstrate that in a subsample of 78 individuals, our SNPs detect known genetic structure more clearly than 24 microsatellites genotyped for the same individuals and that these results are not driven by the SNP ascertainment scheme. Here, we present one of the first large-scale genotyping resources designed for a threatened species

Short reads, circular genome: skimming SOLiD sequence to construct the bighorn sheep mitochondrial genome

As sequencing technology improves, an increasing number of projects aim to generate full genome sequence, even for nonmodel taxa. These projects may be feasibly conducted at lower read depths if the alignment can be aided by previously developed genomic resources from a closely related species. We investigated the feasibility of constructing a complete mitochondrial (mt) genome without preamplification or other targeting of the sequence. Here we present a full mt genome sequence (16,463 nucleotides) for the bighorn sheep (Ovis canadensis) generated though alignment of SOLiD short-read sequences to a reference genome. Average read depth was 1240, and each base was covered by at least 36 reads. We then conducted a phylogenomic analysis with 27 other bovid mitogenomes, which placed bighorn sheep firmly in the Ovis clade. These results show that it is possible to generate a complete mitogenome by skimming a low-coverage genomic sequencing library. This technique will become increasingly applicable as the number of taxa with some level of genome sequence rises

SNP and microsatellite data

RAD, transcriptomic, and microsatellite data used in the manuscript, as well as inputs/outputs of analyses. This archive contains the complete RAD .vcf file. For the complete transcriptomic .vcf file, see doi:10.5061/dryad.606j6, which is described in doi:10.1111/1755-0998.12190

Circumpolar Genetic Structure and Recent Gene Flow of Polar Bears: A Reanalysis.

Recently, an extensive study of 2,748 polar bears (Ursus maritimus) from across their circumpolar range was published in PLOS ONE, which used microsatellites and mitochondrial haplotypes to apparently show altered population structure and a dramatic change in directional gene flow towards the Canadian Archipelago-an area believed to be a future refugium for polar bears as their southernmost habitats decline under climate change. Although this study represents a major international collaborative effort and promised to be a baseline for future genetics work, methodological shortcomings and errors of interpretation undermine some of the study's main conclusions. Here, we present a reanalysis of this data in which we address some of these issues, including: (1) highly unbalanced sample sizes and large amounts of systematically missing data; (2) incorrect calculation of FST and of significance levels; (3) misleading estimates of recent gene flow resulting from non-convergence of the program BayesAss. In contrast to the original findings, in our reanalysis we find six genetic clusters of polar bears worldwide: the Hudson Bay Complex, the Western and Eastern Canadian Arctic Archipelago, the Western and Eastern Polar Basin, and-importantly-we reconfirm the presence of a unique and possibly endangered cluster of bears in Norwegian Bay near Canada's expected last sea-ice refugium. Although polar bears' abundance, distribution, and population structure will certainly be negatively affected by ongoing-and increasingly rapid-loss of Arctic sea ice, these genetic data provide no evidence of strong directional gene flow in response to recent climate change

Data from: Temporal dynamics of linkage disequilibrium in two populations of bighorn sheep

Linkage disequilibrium (LD) is the nonrandom association of alleles at two markers. Patterns of LD have biological implications as well as practical ones when designing association studies or conservation programs aimed at identifying the genetic basis of fitness differences within and among populations. However, the temporal dynamics of LD in wild populations has received little empirical attention. In this study, we examined the overall extent of LD, the effect of sample size on the accuracy and precision of LD estimates, and the temporal dynamics of LD in two populations of bighorn sheep (Ovis canadensis) with different demographic histories. Using over 200 microsatellite loci, we assessed two metrics of multi-allelic LD, D′, and χ′2. We found that both populations exhibited high levels of LD, although the extent was much shorter in a native population than one that was founded via translocation, experienced a prolonged bottleneck post founding, followed by recent admixture. In addition, we observed significant variation in LD in relation to the sample size used, with small sample sizes leading to depressed estimates of the extent of LD but inflated estimates of background levels of LD. In contrast, there was not much variation in LD among yearly cross-sections within either population once sample size was accounted for. Lack of pronounced interannual variability suggests that researchers may not have to worry about interannual variation when estimating LD in a population and can instead focus on obtaining the largest sample size possible

RM_Data

Microsatellite genotypes for the Ram Mountain population

Data from: "Polar bear (Ursus maritimus) transcriptome assembly and SNP discovery" in Genomic Resources Notes accepted 1 August 2013-30 September 2013

Polar bears (Ursus maritimus) in the Western Hudson Bay subpopulation have been declining in size and body condition for decades, as climate change causes earlier sea ice breakup, reduced hunting time on the ice, and an increasingly long fasting season. As Western Hudson Bay females have decreased in size, rates of litter production and average litter size have also decreased, while cub mortality and average time to independence have increased. Although these changes have potential evolutionary consequences, little is yet known about the adaptive genetic variation in body size or fat accumulation that would have to underlie any such change. In this study, we used high-throughput Illumina sequencing to develop SNPs from pooled blood and fat transcriptomes, using samples from five adult female polar bears and five (unrelated) dependent cubs. In total, we generated 371,258 transcripts of which 36,755 were deemed to be “full length” (i.e., covered more than 90% of their best BLAST hit), and we identified 63,020 SNPs. Since this study was conducted, we have used a subset of these SNPs to develop an Illumina BeadArray for quantitative genetics research in Western Hudson Bay

Circumpolar Genetic Structure and Recent Gene Flow of Polar Bears: A Reanalysis

<div><p>Recently, an extensive study of 2,748 polar bears (<i>Ursus maritimus</i>) from across their circumpolar range was published in PLOS ONE, which used microsatellites and mitochondrial haplotypes to apparently show altered population structure and a dramatic change in directional gene flow towards the Canadian Archipelago—an area believed to be a future refugium for polar bears as their southernmost habitats decline under climate change. Although this study represents a major international collaborative effort and promised to be a baseline for future genetics work, methodological shortcomings and errors of interpretation undermine some of the study’s main conclusions. Here, we present a reanalysis of this data in which we address some of these issues, including: (1) highly unbalanced sample sizes and large amounts of systematically missing data; (2) incorrect calculation of <i>F</i>_<i>ST</i> and of significance levels; (3) misleading estimates of recent gene flow resulting from non-convergence of the program BayesAss. In contrast to the original findings, in our reanalysis we find six genetic clusters of polar bears worldwide: the Hudson Bay Complex, the Western and Eastern Canadian Arctic Archipelago, the Western and Eastern Polar Basin, and—importantly—we reconfirm the presence of a unique and possibly endangered cluster of bears in Norwegian Bay near Canada’s expected last sea-ice refugium. Although polar bears’ abundance, distribution, and population structure will certainly be negatively affected by ongoing—and increasingly rapid—loss of Arctic sea ice, these genetic data provide no evidence of strong directional gene flow in response to recent climate change.</p></div

Hierarchical analysis of molecular variance (AMOVA) for mitochondrial DNA among management units within the six genetic clusters identified in this paper and shown in Table 2.

<p>Note that many management units (incl. the entire Norwegian Bay cluster) were excluded entirely from this AMOVA because of inadequate sampling. Because we lacked sample location information for downloaded haplotypes, we were unable to split Davis Strait or the Laptev Sea into northern/southern or eastern/western samples; therefore, these MUs were removed for this calculation in addition to the MUs that were removed for low sample sizes in <a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0148967#pone.0148967.t002" target="_blank">Table 2</a>.</p