Massively parallel sequencing has revolutionized many areas of biology, but sequencing large amounts of DNA in many individuals is cost-prohibitive and unnecessary for many studies. Genomic complexity reduction techniques such as sequence capture and restriction enzyme-based methods enable the analysis of many more individuals per unit cost. Despite their utility, current complexity reduction methods have limitations, especially when large numbers of individuals are analyzed. Here we develop a much improved restriction site-associated DNA (RAD) sequencing protocol and a new method called Rapture (RAD capture). The new RAD protocol improves versatility by separating RAD tag isolation and sequencing library preparation into two distinct steps. This protocol also recovers more unique (nonclonal) RAD fragments, which improves both standard RAD and Rapture analysis. Rapture then uses an in-solution capture of chosen RAD tags to target sequencing reads to desired loci. Rapture combines the benefits of both RAD and sequence capture, i.e., very inexpensive and rapid library preparation for many individuals as well as high specificity in the number and location of genomic loci analyzed. Our results demonstrate that Rapture is a rapid and flexible technology capable of analyzing a very large number of individuals with minimal sequencing and library preparation cost. The methods presented here should improve the efficiency of genetic analysis for many aspects of agricultural, environmental, and biomedical science.

Accounting for specific alleles can be necessary to prevent the loss of significant biodiversity and ecosystem services.

The delineation of conservation units (CUs) is a challenging issue that has profound implications for minimizing the loss of biodiversity and ecosystem services. CU delineation typically seeks to prioritize evolutionary significance and genetic methods play a pivotal role in the delineation process by quantifying overall differentiation between populations. While CUs that primarily reflect overall genetic differentiation do protect adaptive differences between distant populations, they do not necessarily protect adaptive variation within highly connected populations. Advances in genomic methodology facilitate the characterization of adaptive genetic variation, but the potential utility of this information for CU delineation is unclear. Here we use genomic methods to investigate the evolutionary basis of premature migration in Pacific salmon, a complex behavioral and physiological adaptation that exists within highly-connected populations and has experienced severe declines. Strikingly, we find that premature migration is associated with the same single locus across multiple populations in each of two different species. Patterns of variation at this locus suggest that the premature migration alleles arose from a single evolutionary event within each species and were subsequently spread to distant populations through straying and positive selection. Our results reveal that complex adaptive variation can depend on rare mutational events at a single locus, demonstrate

Massively parallel sequencing has revolutionized many areas of biology but sequencing large amounts of DNA in many individuals is cost-prohibitive and unnecessary for many studies.Genomic complexity reduction techniques such as sequence capture and restriction enzyme- (which was not peer-reviewed) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity.

Aim Climate warming is causing extensive loss of glaciers in mountainous regions, yet our understanding of how glacial recession influences evolutionary processes and genetic diversity is limited. Linking genetic structure with the influences shaping it can improve understanding of how species respond to environmental change. Here, we used genome‐scale data and demographic modelling to resolve the evolutionary history of Lednia tumana, a rare, aquatic insect endemic to alpine streams. We also employed a range of widely used data filtering approaches to quantify how they influenced population structure results. Location Alpine streams in the Rocky Mountains of Glacier National Park, Montana, USA. Taxon Lednia tumana, a stonefly (Order Plecoptera) in the family Nemouridae. Methods We generated single nucleotide polymorphism data through restriction‐site associated DNA sequencing to assess contemporary patterns of genetic structure for 11 L. tumana populations. Using identified clusters, we assessed demographic history through model selection and parameter estimation in a coalescent framework. During population structure analyses, we filtered our data to assess the influence of singletons, missing data and total number of markers on results. Results Contemporary patterns of population structure indicate that L. tumana exhibits a pattern of isolation‐by‐distance among populations within three genetic clusters that align with geography. Mean pairwise genetic differentiation (FST) among populations was 0.033. Coalescent‐based demographic modelling supported divergence with gene flow among genetic clusters since the end of the Pleistocene (~13‐17 kya), likely reflecting the south‐to‐north recession of ice sheets that accumulated during the Wisconsin glaciation. Main conclusions We identified a link between glacial retreat, evolutionary history and patterns of genetic diversity for a range‐restricted stonefly imperiled by climate change. This finding included a history of divergence with gene flow, an unexpected conclusion for a mountaintop species. Beyond L. tumana, this study demonstrates the complexity of assessing genetic structure for weakly differentiated species, shows the degree to which rare alleles and missing data may influence results, and highlights the usefulness of genome‐scale data to extend population genetic inquiry in non‐model species.