The development of high‐throughput sequencing technologies is dramatically increasing the use of single nucleotide polymorphisms (
Molecular ecologists seek to genotype hundreds to thousands of loci from hundreds to thousands of individuals at minimal cost per sample. Current methods, such as restriction‐site‐associated
- NSF-PAR ID:
- 10243960
- Publisher / Repository:
- Wiley-Blackwell
- Date Published:
- Journal Name:
- Molecular Ecology Resources
- Volume:
- 16
- Issue:
- 5
- ISSN:
- 1755-098X
- Page Range / eLocation ID:
- p. 1264-1278
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract SNP s) across the field of genetics, but most parentage studies of wild populations still rely on microsatellites. We developed a bioinformatic pipeline for identifyingSNP panels that are informative for parentage analysis from restriction site‐associatedDNA sequencing (RAD seq) data. This pipeline includes options for analysis with or without a reference genome, and provides methods to maximize genotyping accuracy and select sets of unlinked loci that have high statistical power. We test this pipeline on small populations of Mexican gray wolf and bighorn sheep, for which parentage analyses are expected to be challenging due to low genetic diversity and the presence of many closely related individuals. We compare the results of parentage analysis acrossSNP panels generated with or without the use of a reference genome, and betweenSNP s and microsatellites. For Mexican gray wolf, we conducted parentage analyses for 30 pups from a single cohort where samples were available from 64% of possible mothers and 53% of possible fathers, and the accuracy of parentage assignments could be estimated because true identities of parents were known a priori based on field data. For bighorn sheep, we conducted maternity analyses for 39 lambs from five cohorts where 77% of possible mothers were sampled, but true identities of parents were unknown. Analyses with and without a reference genome producedSNP panels with ≥95% parentage assignment accuracy for Mexican gray wolf, outperforming microsatellites at 78% accuracy. Maternity assignments were completely consistent across allSNP panels for the bighorn sheep, and were 74.4% consistent with assignments from microsatellites. Accuracy and consistency of parentage analysis were not reduced when using as few as 284SNP s for Mexican gray wolf and 142SNP s for bighorn sheep, indicating our pipeline can be used to developSNP genotyping assays for parentage analysis with relatively small numbers of loci. -
more » « less
Abstract Next‐generation sequencing technologies now allow researchers of non‐model systems to perform genome‐based studies without the requirement of a (often unavailable) closely related genomic reference. We evaluated the role of restriction endonuclease (
RE ) selection in double‐digest restriction‐site‐associatedDNA sequencing (ddRAD seq) by generating reduced representation genome‐wide data using four differentRE combinations. Our expectation was thatRE selections targeting longer, more complex restriction sites would recover fewer loci thanRE with shorter, less complex sites. We sequenced a diverse sample of non‐model arachnids, including five congeneric pairs of harvestmen (Opiliones) and four pairs of spiders (Araneae). Sample pairs consisted of either conspecifics or closely related congeneric taxa, and in total 26 sample pair analyses were tested. Sequence demultiplexing, read clustering and variant calling were performed in thepy program. The 6‐base pair cutterRAD Eco combined with methylated site‐specific 4‐base pair cutterRI MspI produced, on average, the greatest numbers of intra‐individual loci and shared loci per sample pair. As expected, the number of shared loci recovered for a sample pair covaried with the degree of genetic divergence, estimated with cytochrome oxidase I sequences, although this relationship was non‐linear. Our comparative results will prove useful in guiding protocol selection for ddRAD seq experiments on many arachnid taxa where reference genomes, even from closely related species, are unavailable. -
more » « less
Premise The ability to sequence genome‐scale data from herbarium specimens would allow for the economical development of data sets with broad taxonomic and geographic sampling that would otherwise not be possible. Here, we evaluate the utility of a basic double‐digest restriction site–associated
DNA sequencing (ddRAD seq) protocol usingDNA s from four genera extracted from both silica‐dried and herbarium tissue.Methods DNA s fromDraba ,Boechera ,Solidago , andIlex were processed with a ddRAD seq protocol. The effects ofDNA degradation, taxon, and specimen age were assessed.Results Although taxon, preservation method, and specimen age affected data recovery, large phylogenetically informative data sets were obtained from the majority of samples.
Discussion These results suggest that herbarium samples can be incorporated into dd
RAD seq project designs, and that specimen age can be used as a rapid on‐site guide for sample choice. The detailed protocol we provide will allow users to pursue herbarium‐based ddRAD seq projects that minimize the expenses associated with fieldwork and sample evaluation. -
more » « less
Abstract Reproductive isolation can be initiated by changes in one or a few key traits that prevent random mating among individuals in a population. During the early stages of speciation, when isolation is often incomplete, there will be a heterogeneous pattern of differentiation across regions of the genome between diverging populations, with loci controlling these key traits appearing the most distinct as a result of strong diversifying selection. In this study, we used Illumina‐sequenced dd
RAD tags to identify genomewide patterns of differentiation in three recently diverged island populations of theMonarcha castaneiventris flycatcher of the Solomon Islands. Populations of this species have diverged in plumage colour, and these differences in plumage colour, in turn, are used in conspecific recognition and likely important in reproductive isolation. Previous candidate gene sequencing identified point mutations inMC 1RASIP SNP s are extremely low between populations of distinct plumage colour, with no loci reaching the level of differentiation found in either candidate gene. Further, we found that a phylogenetic analysis based on theseSNP s produced a taxonomy wherein the two melanic populations appear to have evolved convergently, rather than from a single common ancestor, in contrast to their original classification as a single subspecies. Finally, we found evidence that the pattern of low genomic differentiation is the result of both incomplete lineage sorting and gene flow between populations. -
more » « less
Abstract Identification of genes underlying genomic signatures of natural selection is key to understanding adaptation to local conditions. We used targeted resequencing to identify
SNP markers in 5321 candidate adaptive genes associated with known immunological, metabolic and growth functions in ovids and other ungulates. We selectively targeted 8161 exons in protein‐coding and nearby 5′ and 3′ untranslated regions of chosen candidate genes. Targeted sequences were taken from bighorn sheep (Ovis canadensis ) exon capture data and directly from the domestic sheep genome (Ovis aries v. 3; oviAri3). The bighorn sheep sequences used in the Dall's sheep (Ovis dalli dalli ) exon capture aligned to 2350 genes on the oviAri3 genome with an average of 2 exons each. We developed a microfluidic qPCR‐basedSNP chip to genotype 476 Dall's sheep from locations across their range and test for patterns of selection. Using multiple corroborating approaches (lositan andbayescan ), we detected 28SNP loci potentially under selection. We additionally identified candidate loci significantly associated with latitude, longitude, precipitation and temperature, suggesting local environmental adaptation. The three methods demonstrated consistent support for natural selection on nine genes with immune and disease‐regulating functions (e.g. Ovar‐DRA ,APC ,BATF 2,MAGEB 18), cell regulation signalling pathways (e.g.KRIT 1,PI 3K,ORRC 3), and respiratory health (CYSLTR 1). Characterizing adaptive allele distributions from novel genetic techniques will facilitate investigation of the influence of environmental variation on local adaptation of a northern alpine ungulate throughout its range. This research demonstrated the utility of exon capture for gene‐targetedSNP discovery and subsequentSNP chip genotyping using low‐quality samples in a nonmodel species.
