Find research datasets worth reusing

opencc-zeroDec 2013View details →

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.

opencc-zeroDec 2015View details →

Data from: Discovery of 20,000 RAD–SNPs and development of a 52-SNP array for monitoring river otters

Many North American river otter (Lontra canadensis) populations are threatened or recovering but are difficult to study because they occur at low densities, it is difficult to visually identify individuals, and they inhabit aquatic environments that accelerate degradation of biological samples. Single nucleotide polymorphisms (SNPs) can improve our ability to monitor demographic and genetic parameters of difficult to study species. We used restriction site associated DNA (RAD) sequencing to discover 20,772 SNPs present in Montana, USA, river otter populations, including 14,512 loci that were also variable in at least one other population range-wide. After applying careful filtering criteria meant to minimize ascertainment bias and identify high quality, highly heterozygous (H o = 0.2–0.50) SNPs, we developed and tested 52 independent SNP qPCR genotyping assays, including 41 that performed well with diluted DNA. The 41 loci provided high power for population assignment tests with only 1 misassignment (1.6 %) between closely neighboring populations. Our SNPs showed high power to differentiate individuals and assign them to population of origin, as well as strong concordance of genotypes from high and diluted concentrations of DNA, and between original RAD and the SNP qPCR array.

Data from: RAD sequencing resolves fine-scale population structure in a benthic invertebrate: implications for understanding phenotypic plasticity

The field of molecular ecology is transitioning from the use of small panels of classical genetic markers such as microsatellites to much larger panels of single nucleotide polymorphisms (SNPs) generated by approaches like RAD sequencing. However, few empirical studies have directly compared the ability of these methods to resolve population structure. This could have implications for understanding phenotypic plasticity, as many previous studies of natural populations may have lacked the power to detect genetic differences, especially over micro-geographic scales. We therefore compared the ability of microsatellites and RAD sequencing to resolve fine-scale population structure in a commercially important benthic invertebrate by genotyping great scallops (Pecten maximus) from nine populations around Northern Ireland at 13 microsatellites and 10 539 SNPs. The shells were then subjected to morphometric and colour analysis in order to compare patterns of phenotypic and genetic variation. We found that RAD sequencing was superior at resolving population structure, yielding higher Fst values and support for two distinct genetic clusters, whereas only one cluster could be detected in a Bayesian analysis of the microsatellite dataset. Furthermore, appreciable phenotypic variation was observed in size-independent shell shape and coloration, including among localities that could not be distinguished from one another genetically, providing support for the notion that these traits are phenotypically plastic. Taken together, our results suggest that RAD sequencing is a powerful approach for studying population structure and phenotypic plasticity in natural populations.

opencc-zeroDec 2015View details →

Data from: Bioinformatic processing of RAD-seq data dramatically impacts downstream population genetic inference

Restriction site-associated DNA sequencing (RAD-seq) provides high-resolution population genomic data at low cost, and has become an important component in ecological and evolutionary studies. As with all high-throughput technologies, analytic strategies require critical validation to ensure accurate and unbiased interpretation. To test for the impact of bioinformatic data processing on downstream population genetic inferences, we analysed mammalian RAD-seq data (>100 individuals) with 312 combinations of methodology (de novo vs. mapping to references of increasing divergence) and filtering criteria (missing data, HWE, FIS, coverage, mapping, genotype quality). In an effort to identify commonalities and biases in all pipelines, we computed summary statistics (nr. loci, nr. SNP, π, Hetobs, FIS, FST, Ne, m) and compared the results to independent null expectations (isolation-by-distance correlation, expected transition-to-transversion ratio Ts/Tv, Mendelian mismatch rates of known parent-offspring trios). We observed large differences between reference-based and de novo approaches, the former generally calling more SNPs and reducing FIS and Ts/Tv. Data completion levels showed little impact on most summary statistics, and FST estimates were robust across all pipelines. The site-frequency spectrum (SFS) was highly sensitive to the chosen approach as reflected in large variance of parameter estimates across demographic scenarios (single-population bottlenecks and isolation-with-migration model). Null-expectations were best met by reference-based approaches, though contingent on the specific criteria. We recommend RAD-seq studies employ reference-based approaches to a closely related genome, and due to the high stochasticity associated with the pipeline advocate the use of multiple pipelines to ensure robust population genetic and demographic inferences.

Data from: RAD sequencing reveals genomewide divergence between independent invasions of the European green crab (Carcinus maenas) in the Northwest Atlantic

Genomic studies of invasive species can reveal both invasive pathways and functional differences underpinning patterns of colonization success. The European green crab (Carcinus maenas) was initially introduced to eastern North America nearly 200 years ago where it expanded northwards to eastern Nova Scotia. A subsequent invasion to Nova Scotia from a northern European source allowed further range expansion, providing a unique opportunity to study the invasion genomics of a species with multiple invasions. Here, we use restriction-site-associated DNA sequencing-derived SNPs to explore fine-scale genomewide differentiation between these two invasions. We identified 9137 loci from green crab sampled from 11 locations along eastern North America and compared spatial variation to mitochondrial COI sequence variation used previously to characterize these invasions. Overall spatial divergence among invasions was high (pairwise FST ~0.001 to 0.15) and spread across many loci, with a mean FST ~0.052 and 52% of loci examined characterized by FST values >0.05. The majority of the most divergent loci (i.e., outliers, ~1.2%) displayed latitudinal clines in allele frequency highlighting extensive genomic divergence among the invasions. Discriminant analysis of principal components (both neutral and outlier loci) clearly resolved the two invasions spatially and was highly correlated with mitochondrial divergence. Our results reveal extensive cryptic intraspecific genomic diversity associated with differing patterns of colonization success and demonstrates clear utility for genomic approaches to delineating the distribution and colonization success of aquatic invasive species.

Data from: Sturgeon conservation genomics: SNP discovery and validation using RAD sequencing

Caviar-producing sturgeons belonging to the genus Acipenser are considered to be one of the most endangered species groups in the world. Continued overfishing in spite of increasing legislation, zero catch quotas and extensive aquaculture production have led to the collapse of wild stocks across Europe and Asia. The evolutionary relationships among Adriatic, Russian, Persian and Siberian sturgeons are complex because of past introgression events and remain poorly understood. Conservation management, traceability and enforcement suffer a lack of appropriate DNA markers for the genetic identification of sturgeon at the species, population and individual level. This study employed RAD sequencing to discover and characterize single nucleotide polymorphism (SNP) DNA markers for use in sturgeon conservation in these four tetraploid species over three biological levels, using a single sequencing lane. Four population meta-samples and eight individual samples from one family were barcoded separately before sequencing. Analysis of 14.4 Gb of paired-end RAD data focused on the identification of SNPs in the paired-end contig, with subsequent in silico and empirical validation of candidate markers. Thousands of putatively informative markers were identified including, for the first time, SNPs that show population-wide differentiation between Russian and Persian sturgeons, representing an important advance in our ability to manage these cryptic species. The results highlight the challenges of genotyping-by-sequencing in polyploid taxa, while establishing the potential genetic resources for developing a new range of caviar traceability and enforcement tools.

Data from: Applicability of RAD-tag genotyping for inter-familial comparisons: empirical data from two cetaceans

Restriction site-Associated DNA tag (RAD-tag) sequencing has become a popular approach to generate thousands of SNPs used to address diverse questions in population genomics. Comparatively, the suitability of RAD-tag genotyping to address evolutionary questions across divergent species has been the subject of only a few recent studies. Here, we evaluate the applicability of this approach to conduct genome-wide scans for polymorphisms across two cetacean species belonging to distinct families: the short-beaked common dolphin (Delphinus delphis; n = 5 individuals) and the harbor porpoise (Phocoena phocoena; n = 1 individual). Additionally, we explore the effects of varying two parameters in the Stacks analysis pipeline on the number of loci and level of divergence obtained. We observed a 34% drop in the total number of loci that were present in all individuals when analyzing individuals from the distinct families compared to analyses restricted to intra-specific comparisons (i.e., within D. delphis). Despite relatively stringent quality filters, 3,595 polymorphic loci were retrieved from our inter-familial comparison. Cetaceans have undergone rapid diversification and the estimated divergence time between the two families is relatively recent (14 to 19 My). Thus, our results showed that, for this level of divergence, a large number of orthologous loci can still be genotyped using this approach, which is on par with two recent in silico studies. Our findings constitute one of the first empirical investigations using RAD-tag sequencing at this level of divergence and highlights the great potential of this approach in comparative studies and to address evolutionary questions.

Data from: Mass production of SNP markers in a nonmodel passerine bird through RAD sequencing and contig mapping to the zebra finch genome

Here, we present an adaptation of restriction-site-associated DNA sequencing (RAD-seq) to the Illumina HiSeq2000 technology that we used to produce SNP markers in very large quantities at low cost per unit in the Réunion grey white-eye (Zosterops borbonicus), a nonmodel passerine bird species with no reference genome. We sequenced a set of six pools of 18–25 individuals using a single sequencing lane. This allowed us to build around 600 000 contigs, among which at least 386 000 could be mapped to the zebra finch (Taeniopygia guttata) genome. This yielded more than 80 000 SNPs that could be mapped unambiguously and are evenly distributed across the genome. Thus, our approach provides a good illustration of the high potential of paired-end RAD sequencing of pooled DNA samples combined with comparative assembly to the zebra finch genome to build large contigs and characterize vast numbers of informative SNPs in nonmodel passerine bird species in a very efficient and cost-effective way.

opencc-zeroDec 2014View details →

Data from: Extending RAD tag analysis to microbial ecology: a comparison between multi locus sequence typing (MLST) and 2b-RAD to investigate Listeria monocytogenes genetic structure

The advent of next-generation sequencing (NGS) has dramatically changed bacterial typing technologies, increasing our ability to differentiate bacterial isolates. Despite it is now possible to sequence a bacterial genome in a few days and at reasonable costs, most genetic analyses do not require whole-genome sequencing, which also remains impractical for large population samples due to the cost of individual library preparation and bioinformatics. More traditional sequencing approaches, however, such as MultiLocus Sequence Typing (mlst) are quite laborious and time-consuming, especially for large-scale analyses. In this study, a genotyping approach based on restriction site-associated (RAD) tag sequencing, 2b-RAD, was applied to characterize Listeria monocytogenes strains. To verify the feasibility of the method, an in silico analysis was performed on 30 available complete genomes. For the same set of strains, in silico mlst analysis was conducted as well. Subsequently, 2b-RAD and mlst analyses were experimentally carried out on 58 isolates collected from food samples or food-processing sites. The obtained results demonstrate that 2b-RAD predicts mlst types and often provides more detailed information on population structure than mlst. Moreover, the majority of variants differentiating identical sequence type isolates mapped against accessory fragments, thus providing additional information to characterize strains. Although mlst still represents a reliable typing method, large-scale studies on molecular epidemiology and public health, as well as bacterial phylogenetics, population genetics and biosafety could benefit of a low cost and fast turnaround time approach such as the 2b-RAD analysis proposed here.

opencc-zeroFeb 2024View details →

Obovaria olivaria maf filtered vcf file from: RAD-tag and mitochondrial DNA sequencing reveal the genetic structure of a widespread and regionally imperiled freshwater mussel, Obovaria olivaria (Bivalvia: Unionidae)

Obovaria olivaria is a species of freshwater mussel native to the Mississippi River and Laurentian Great Lakes-St. Lawrence River drainages of North America. This mussel has experienced population declines across large parts of its distribution and is imperiled in many jurisdictions. Obovaria olivaria uses the similarly imperiled Acipenser fulvescens (Lake Sturgeon) as a host for its glochidia. We employed mitochondrial DNA sequencing and Restriction-site Associated DNA sequencing (RAD-seq) to assess patterns of genetic diversity and population structure of O. olivaria from 19 collection locations including the St. Lawrence River drainage, the Great Lakes drainage, the Upper Mississippi River drainage, the Ohioan River drainage and the Mississippi Embayment. Heterozygosity was highest in Upper Mississippi and Great Lakes populations, followed by a reduction in diversity and relative effective population size in the St. Lawrence populations. Pairwise FST ranged from 0.00 to 0.20, and analyses of genetic structure revealed two major ancestral populations, one including all St. Lawrence River/Ottawa River sites and the other including remaining sites; however, significant admixture and isolation by river distance across the range were evident. The genetic diversity and structure of O. olivaria is consistent with the existing literature on Acipenser fulvescens and suggest that, although northern and southern O. olivaria populations are genetically distinct, genetic structure in O. olivaria is largely clinal rather than discrete across its range. Conservation and restoration efforts of O. olivaria should prioritize the maintenance and restoration of locations where O. olivaria remain, especially in northern rivers, and to ensure connectivity that will facilitate dispersal of Acipenser fulvescens and movement of encysted glochidia.

Data from: Demographic inference from whole-genome and RAD sequencing data suggests alternating human impacts on goose populations since the last ice age

We investigated how population changes and fluctuations in the pink-footed goose might have been affected by climatic and anthropogenic factors. First, genomic data confirmed the existence of two separate populations: western (Iceland) and eastern (Svalbard/Denmark). Second, emographic inference suggests that the species survived the last glacial period as a single ancestral population with a low population size (100-1,000 individuals) that split into the current populations at the end of the Last Glacial Maximum with Iceland being the most plausible glacial refuge. While population changes during the last glaciation were clearly environmental, we hypothesize that more recent demographic changes are human-related: (1) the inferred population increase in the Neolithic is due to deforestation to establish new lands for agriculture, increasing available habitat for pink-footed geese (2) the decline inferred during the Middle Ages is due to human persecution and (3) improved protection explains the increasing demographic trends during the 20th century. Our results suggest both environmental (during glacial cycles) and anthropogenic effects (more recent) can be a threat to species survival.

opencc-by-4.0Nov 2024View details →

zenodo32/100

RaD_SoilRespiration_2020-2023_CodeRepository

Open the record for dataset details and reuse information.

opencc-zeroDec 2017View details →

Data from: Exploring the phylogeography of a hexaploid freshwater fish by RAD sequencing

The KwaZulu-Natal yellowfish (Labeobarbus natalensis) is an abundant cyprinid, endemic to KwaZulu-Natal Province, South Africa. In this study we developed a Single Nucleotide Polymorphism (SNP) dataset from double-digest Restriction-site Associated DNA (ddRAD) sequencing of samples across the distribution. We addressed several hidden challenges, primarily focussing on proper filtering of RAD data and selecting optimal parameters for data processing in polyploid lineages. We used the resulting high-quality SNP dataset to investigate the population genetic structure of L. natalensis. A small number of mitochondrial markers present in these data had disproportionate influence on the recovered genetic structure. The presence of singleton SNPs also confounded genetic structure. We found a well-supported division into northern and southern lineages, with further subdivision into five populations, one of which reflects north-south admixture. Approximate Bayesian Computation scenario testing supported a scenario where an ancestral population diverged into northern and southern lineages, which then diverged to yield the current five populations. All river systems showed similar levels of genetic diversity, which appears unrelated to drainage system size. Nucleotide diversity was highest in the smallest river system, the Mbokodweni, which, together with adjacent small coastal systems, should be considered as a key catchment for conservation.

opencc-zeroAug 2022View details →

Phylogenomics and biogeography of Castanea (chestnut) and Hamamelis (witch-hazel): Choosing between RAD-seq and Hyb-Seq approaches

Target enrichment and RAD-seq are well-established high throughput sequencing technologies that have been increasingly used for phylogenomic studies. Each method has its own pros and cons. The choice between them is a practical issue for plant systematists studying the evolutionary histories of biodiversity of relatively recent origins. However, few studies have compared the congruence and conflict between results from the two methods within the same group of organisms in plants. In this study, we employed RAD-seq and Hyb-Seq of Angiosperm 353 genes in phylogenomic and biogeographic studies of Hamamelis (the witch-hazels) and Castanea (chestnuts), two classic examples exhibiting the well-known eastern Asian (EA)-eastern North American (ENA) disjunct distribution, and compared them side by side. Our results showed congruences in phylogenetic inference and divergence time dating between the two data sets obtained through our customized procedures of library preparation and sequence trimming, although they differed in the number of loci and informative sites, the amount of missing data, and sampling within species. We suggest the selection of the two methods based on fund availability and sampling scale. Our phylogenetic analyses of RAD-seq and Hyb-Seq data resulted in well-resolved species relationships, and ancient introgressions were revealed in both genera by D-statistic test and PhyloNet. Biogeographic analyses including fossil data using total evidence-based dated tree and DEC model, applying specific inter-area dispersal probabilities, revealed a complicated history for each genus, indicating multiple intercontinental dispersals and local extinctions in areas outside of the taxa's modern ranges in both the Paleogene and Neogene. The study demonstrates the importance of including fossil taxa for a more accurate reconstruction of biogeographic histories of taxa to understand the EA and ENA floristic disjunction. Our results support an "out of western North America" migration of Castanea but an "out of Asia" migration of Hamamelis during their initial diversification, and the origins of the EA-ENA disjunction in both genera were results of vicariance.

opencc-zeroSep 2022View details →

RAD-seq of Ash individuals

Ash (Fraxinus spp.) is one of the most widely distributed tree genera in North America. Populations of ash in the United States and Canada have been decimated by the introduced pest Agrilus planipennis (Coleoptera: Buprestidae; emerald ash borer), having negative impacts on both forest ecosystems and economic interests. The majority of trees succumb to attack by A. planipennis, but some trees have been found to be tolerant to infestation despite years of exposure. Restriction site-associated DNA (RAD) sequencing was used to sequence ash individuals, both tolerant and susceptible to A. planipennis attack, in order to identify single nucleotide polymorphism (SNP) patterns related to tolerance and health declines. Understanding the mechanisms of host tolerance through genome-wide association has the potential to restore populations with cultivars that are able to withstand A. planipennis infestation. We were successful in using RAD-sequencing in order to identify SNPs that could contribute to tolerance of A. planipennis. This was a first step toward uncovering the genetic basis for host tolerance to A. planipennis.

opencc-zeroDec 2014View details →

Data from: "Identification of SNP markers for the endangered Ugandan red colobus (Procolobus rufomitratus tephrosceles) using RAD sequencing" in Genomic Resources Notes accepted 1 December 2014 to 31 January 2015

Despite dramatic growth in the field of primate genomics over the past decade, studies of primate population and conservation genomics in the wild have been hampered due to the difficulties inherent in studying non-model organisms and endangered species, such as lack of a reference genome and current challenges in de novo primate genome assembly. Here, we used Restriction-site Associated DNA (RAD) sequencing to develop a population-based SNP panel for the Ugandan red colobus (P. rufomitratus tephrosceles), which is a highly threatened monkey due to habitat loss. We analyzed blood samples from 24 individuals from Kibale National Park (Uganda) using single-end RAD sequencing. We obtained 70,773,857 reads, of which 58,814,906 passed the filtering steps. Using the program STACKS v. 1.11 we identified 113,376 loci, of which 50,558 were polymorphic and had a mean observed heterozygosity of 0.25. These data will be used to study the effects of habitat fragmentation on genomic diversity, dispersal, and disease transmission in this species. Our approach provides a good example of the potential of RAD sequencing in studies of wild primate populations.

opencc-zeroDec 2011View details →

Data from: Genome-wide RAD sequence data provide unprecedented resolution of species boundaries and relationships in the Lake Victoria cichlid adaptive radiation

Although population genomic studies using next generation sequencing (NGS) data are becoming increasingly common, studies focusing on phylogenetic inference using these data are in their infancy. Here, we use NGS data generated from reduced representation genomic libraries of restriction-site-associated DNA (RAD) markers to infer phylogenetic relationships among 16 species of cichlid fishes from a single rocky island community within Lake Victoria's cichlid adaptive radiation. Previous attempts at sequence-based phylogenetic analyses in Victoria cichlids have shown extensive sharing of genetic variation among species and no resolution of species or higher-level relationships. These patterns have generally been attributed to the very recent origin (<15 000 years) of the radiation, and ongoing hybridization between species. We show that as we increase the amount of sequence data used in phylogenetic analyses, we produce phylogenetic trees with unprecedented resolution for this group. In trees derived from our largest data supermatrices (3 to >5.8 million base pairs in width), species are reciprocally monophyletic with high bootstrap support, and the majority of internal branches on the tree have high support. Given the difficulty of the phylogenetic problem that the Lake Victoria cichlid adaptive radiation represents, these results are striking. The strict interpretation of the topologies we present here warrants caution because many questions remain about phylogenetic inference with very large genomic data set and because we can with the current analysis not distinguish between effects of shared ancestry and post-speciation gene flow. However, these results provide the first conclusive evidence for the monophyly of species in the Lake Victoria cichlid radiation and demonstrate the power that NGS data sets hold to resolve even the most difficult of phylogenetic challenges.

opencc-zeroDec 2013View details →

Data from: Use of RAD sequencing for delimiting species

RAD-tag sequencing is a promising method for conducting genome-wide evolutionary studies. However, to date, only a handful of studies empirically tested its applicability above the species level. In this communication, we use RAD tags to contribute to the delimitation of species within a diverse genus of deep-sea octocorals, Chrysogorgia, for which few classical genetic markers have proved informative. Previous studies have hypothesized that single mitochondrial haplotypes can be used to delimit Chrysogorgia species. On the basis of two lanes of Illumina sequencing, we inferred phylogenetic relationships among 12 putative species that were delimited using mitochondrial data, comparing two RAD analysis pipelines (Stacks and PyRAD). The number of homologous RAD loci decreased dramatically with increasing divergence, as >70% of loci are lost when comparing specimens separated by two mutations on the 700-nt long mitochondrial phylogeny. Species delimitation hypotheses based on the mitochondrial mtMutS gene are largely supported, as six out of nine putative species represented by more than one colony were recovered as discrete, well-supported clades. Significant genetic structure (correlating with geography) was detected within one putative species, suggesting that individuals characterized by the same mtMutS haplotype may belong to distinct species. Conversely, three mtMutS haplotypes formed one well-supported clade within which no population structure was detected, also suggesting that intraspecific variation exists at mtMutS in Chrysogorgia. Despite an impressive decrease in the number of homologous loci across clades, RAD data helped us to fine-tune our interpretations of classical mitochondrial markers used in octocoral species delimitation, and discover previously undetected diversity.

Data from: HyRAD-X, a versatile method combining exome capture and RAD sequencing to extract genomic information from ancient DNA

Over the last decade, protocols aimed at reproducibly sequencing reduced-genome subsets in non-model organisms have been widely developed. Their use is however limited to DNA of relatively high molecular weight. During the last year, several methods exploiting hybridization capture using probes based on RAD-sequencing loci have circumvented this limitation and opened avenues to the study of samples characterized by degraded DNA, such as historical specimens. Here, we present a major update to those methods, namely Hybridization capture from RAD-derived probes obtained from a reduced eXome template (hyRAD-X), a technique applying RAD-sequencing to messenger RNA from one or few fresh specimens to elaborate bench-top produced probes, i.e., a reduced representation of the exome, further used to capture homologous DNA from a samples set. In contrast to previous hybridization-capture methods, the reference catalog on which reads are aligned does not rely on de novo assembly of anonymous RAD-sequencing loci, but on an assembled transcriptome obtained from RNAseq data, thus increasing the accuracy of loci definition and Single-Nucleotide-Polmorphisms (SNP) call, and targeting, specifically, expressed genes. Finally, the capture step of hyRAD-X relies on RNA probes, increasing stringency of hybridization, making it well suited for low-content DNA samples. As a proof of concept, we applied hyRAD-X to subfossil needles from the coniferous tree Abies alba, collected in lake sediments (Origlio, Switzerland) and dating back from 7200-5800 years before present (BP). More specifically we investigated genetic variation before, during, and after an anthropogenic perturbation that caused an abrupt decrease in Abies alba population size, 6500-6200 years BP. HyRAD-X produced a matrix encompassing 524 exome-derived SNPs. Despite a lower observed heterozygosity was observed during the 6.500-6.200 years BP time slice, genetic composition was nearly identical before and after the perturbation, indicating that re-expansion of the population after the decline was driven by autochthonous specimens. To the best of our knowledge, this is the first time a population genomic study incorporating ancient DNA samples of tree subfossils is conducted at a moderate cost using reproducible exome-reduced complexity.

Data from: Revisiting comparisons of genetic diversity in stable and declining species: assessing genome-wide polymorphism in North American bumble bees using RAD sequencing

Genetic variation is of key importance for a species' evolutionary potential, and its estimation is a major component of conservation studies. New DNA sequencing technologies have enabled the analysis of large portions of the genome in nonmodel species, promising highly accurate estimates of such population genetic parameters. Restriction site-associated DNA sequencing (RADseq) is used to analyse thousands of variants in the bumble bee species Bombus impatiens, which is common, and Bombus pensylvanicus, which is in decline. Previous microsatellite-based analyses have shown that gene diversity is lower in the declining B. pensylvanicus than in B. impatiens. RADseq nucleotide diversities appear much more similar in the two species. Both species exhibit allele frequencies consistent with historical population expansions. Differences in diversity observed at microsatellites thus do not appear to have arisen from long-term differences in population size and are either recent in origin or may result from mutational processes. Additional research is needed to explain these discrepancies and to investigate the best ways to integrate next-generation sequencing data and more traditional molecular markers in studies of genetic diversity.