Jesse N. Weber

Double Digest RADseq: An Inexpensive Method for De Novo SNP Discovery and Genotyping in Model and Non-Model Species

The ability to efficiently and accurately determine genotypes is a keystone technology in modern genetics, crucial to studies ranging from clinical diagnostics, to genotype-phenotype association, to reconstruction of ancestry and the detection of selection. To date, high capacity, low cost genotyping has been largely achieved via “SNP chip” microarray-based platforms which require substantial prior knowledge of both genome sequence and variability, and once designed are suitable only for those targeted variable nucleotide sites. This method introduces substantial ascertainment bias and inherently precludes detection of rare or population-specific variants, a major source of information for both population history and genotype-phenotype association. Recent developments in reduced-representation genome sequencing experiments on massively parallel sequencers (commonly referred to as RAD-tag or RADseq) have brought direct sequencing to the problem of population genotyping, but increased cost and procedural and analytical complexity have limited their widespread adoption. Here, we describe a complete laboratory protocol, including a custom combinatorial indexing method, and accompanying software tools to facilitate genotyping across large numbers (hundreds or more) of individuals for a range of markers (hundreds to hundreds of thousands). Our method requires no prior genomic knowledge and achieves per-site and per-individual costs below that of current SNP chip technology, while requiring similar hands-on time investment, comparable amounts of input DNA, and downstream analysis times on the order of hours. Finally, we provide empirical results from the application of this method to both genotyping in a laboratory cross and in wild populations. Because of its flexibility, this modified RADseq approach promises to be applicable to a diversity of biological questions in a wide range of organisms.

Adaptive Variation in Beach Mice Produced by Two Interacting Pigmentation Genes

Little is known about the genetic basis of ecologically important morphological variation such as the diverse color patterns of mammals. Here we identify genetic changes contributing to an adaptive difference in color pattern between two subspecies of oldfield mice (Peromyscus polionotus). One mainland subspecies has a cryptic dark brown dorsal coat, while a younger beach-dwelling subspecies has a lighter coat produced by natural selection for camouflage on pale coastal sand dunes. Using genome-wide linkage mapping, we identified three chromosomal regions (two of major and one of minor effect) associated with differences in pigmentation traits. Two candidate genes, the melanocortin-1 receptor (Mc1r) and its antagonist, the Agouti signaling protein (Agouti), map to independent regions that together are responsible for most of the difference in pigmentation between subspecies. A derived mutation in the coding region of Mc1r, rather than change in its expression level, contributes to light pigmentation. Conversely, beach mice have a derived increase in Agouti mRNA expression but no changes in protein sequence. These two genes also interact epistatically: the phenotypic effects of Mc1r are visible only in genetic backgrounds containing the derived Agouti allele. These results demonstrate that cryptic coloration can be based largely on a few interacting genes of major effect.

Discrete genetic modules are responsible for complex burrow evolution in Peromyscus mice

Relative to morphological traits, we know little about how genetics influence the evolution of complex behavioural differences in nature. It is unclear how the environment influences natural variation in heritable behaviour, and whether complex behavioural differences evolve through few genetic changes, each affecting many aspects of behaviour, or through the accumulation of several genetic changes that, when combined, give rise to behavioural complexity. Here we show that in nature, oldfield mice (Peromyscus polionotus) build complex burrows with long entrance and escape tunnels, and that burrow length is consistent across populations, although burrow depth varies with soil composition. This burrow architecture is in contrast with the small, simple burrows of its sister species, deer mice (P. maniculatus). When investigated under laboratory conditions, both species recapitulate their natural burrowing behaviour. Genetic crosses between the two species reveal that the derived burrows of oldfield mice are dominant and evolved through the addition of multiple genetic changes. In burrows built by first-generation backcross mice, entrance-tunnel length and the presence of an escape tunnel can be uncoupled, suggesting that these traits are modular. Quantitative trait locus analysis also indicates that tunnel length segregates as a complex trait, affected by at least three independent genetic regions, whereas the presence of an escape tunnel is associated with only a single locus. Together, these results suggest that complex behaviours—in this case, a classic ‘extended phenotype’—can evolve through multiple genetic changes each affecting distinct behaviour modules.

Demystifying the RAD fad

We are writing in response to the population and phylogenomics meeting review by Andrews & Luikart ( ) entitled ‘Recent novel approaches for population genomics data analysis’. Restriction‐site‐associated DNA (RAD) sequencing has become a powerful and useful approach in molecular ecology, with several different published methods now available to molecular ecologists, none of which can be considered the best option in all situations. A&L report that the original RAD protocol of Miller et al. ( ) and Baird et al. ( ) is superior to all other RAD variants because putative PCR duplicates can be identified (see Baxter et al. ), thereby reducing the impact of PCR artefacts on allele frequency estimates (Andrews & Luikart ). In response, we (i) challenge the assertion that the original RAD protocol minimizes the impact of PCR artefacts relative to that of other RAD protocols, (ii) present additional biases in RADseq that are at least as important as PCR artefacts in selecting a RAD protocol and (iii) highlight the strengths and weaknesses of four different approaches to RADseq which are a representative sample of all RAD variants: the original RAD protocol (mbRAD, Miller et al. ; Baird et al. ), double digest RAD (ddRAD, Peterson et al. ), ezRAD (Toonen et al. ) and 2bRAD (Wang et al. ). With an understanding of the strengths and weaknesses of different RAD protocols, researchers can make a more informed decision when selecting a RAD protocol.

The evolution of burrowing behaviour in deer mice (genus Peromyscus)

The evolutionary history of most behaviours remains unknown. Here, we assay burrowing behaviour of seven species of deer mice in standardized environments to determine how burrowing evolved in this genus (Peromyscus). We found that several, but not all, species burrow even after many generations of captive breeding. Specifically, there were significant and repeatable differences in both the frequency of burrowing and burrow shape between species. Moreover, these observed species-specific behaviours resemble those reported in wild mice. These results suggest that there is probably a strong genetic component to burrowing in deer mice. We also generated a phylogeny for these seven species using characters from four mtDNA and two autosomal loci. Mapping burrowing behaviour onto this phylogeny suggests a sequence for how complex burrowing evolves: from small, simple burrows to long, multitunnel burrows with defined entrance and escape tunnels. In particular, the most ‘complex’ burrows of P. polionotus appear to be derived. These behavioural data, when examined in a phylogenetic context, show that even closely related species differ in their burrowing behaviours and that the most complex burrows probably evolved by the gradual accumulation of genetic change over time.

Contrasting effects of environment and genetics generate a continuum of parallel evolution

Parallel evolution of similar traits by independent populations in similar environments is considered strong evidence for adaptation by natural selection. Often, however, replicate populations in similar environments do not all evolve in the same way, thus deviating from any single, predominant outcome of evolution. This variation might arise from non-adaptive, population-specific effects of genetic drift, gene flow or limited genetic variation. Alternatively, these deviations from parallel evolution might also reflect predictable adaptation to cryptic environmental heterogeneity within discrete habitat categories. Here, we show that deviations from parallel evolution are the consequence of environmental variation within habitats combined with variation in gene flow. Threespine stickleback (Gasterosteus aculeatus) in adjoining lake and stream habitats (a lake–stream ‘pair’) diverge phenotypically, yet the direction and magnitude of this divergence is not always fully parallel among 16 replicate pairs. We found that the multivariate direction of lake–stream morphological divergence was less parallel between pairs whose environmental differences were less parallel. Thus, environmental heterogeneity among lake–stream pairs contributes to deviations from parallel evolution. Additionally, likely genomic targets of selection were more parallel between environmentally more similar pairs. In contrast, variation in the magnitude of lake–stream divergence (independent of direction) was better explained by differences in lake–stream gene flow; pairs with greater lake–stream gene flow were less morphologically diverged. Thus, both adaptive and non-adaptive processes work concurrently to generate a continuum of parallel evolution across lake–stream stickleback population pairs.

Evaluation of TagSeq, a reliable low‐cost alternative for RNAseq

RNAseq is a relatively new tool for ecological genetics that offers researchers insight into changes in gene expression in response to a myriad of natural or experimental conditions. However, standard RNAseq methods (e.g., Illumina TruSeq® or NEBNext®) can be cost prohibitive, especially when study designs require large sample sizes. Consequently, RNAseq is often underused as a method, or is applied to small sample sizes that confer poor statistical power. Low cost RNAseq methods could therefore enable far greater and more powerful applications of transcriptomics in ecological genetics and beyond. Standard mRNAseq is costly partly because one sequences portions of the full length of all transcripts. Such whole‐mRNA data are redundant for estimates of relative gene expression. TagSeq is an alternative method that focuses sequencing effort on mRNAs’ 3’ end, reducing the necessary sequencing depth per sample, and thus cost. We present a revised TagSeq library construction procedure, and compare its performance against NEBNext®, the ‘gold‐standard’ whole mRNAseq method. We built both TagSeq and NEBNext® libraries from the same biological samples, each spiked with control RNAs. We found that TagSeq measured the control RNA distribution more accurately than NEBNext®, for a fraction of the cost per sample (~10%). The higher accuracy of TagSeq was particularly apparent for transcripts of moderate to low abundance. Technical replicates of TagSeq libraries are highly correlated, and were correlated with NEBNext® results. Overall, we show that our modified TagSeq protocol is an efficient alternative to traditional whole mRNAseq, offering researchers comparable data at greatly reduced cost.

Five Hundred Microsatellite Loci for Peromyscus

Mice of the genus Peromyscus, including several endangered subspecies, occur throughout North America and have been important models for conservation research. We describe 526 primer pairs that amplify microsatellite DNA loci for Peromyscus maniculatus bairdii, 467 of which also amplify in Peromyscus polionotus subgriseus. For 12 of these loci, we report diversity data from a natural population. These markers will be an important resource for future genomic studies of Peromyscus evolution and mammalian conservation.

Partitioning the effects of isolation by distance, environment, and physical barriers on genomic divergence between parapatric threespine stickleback.

Genetic divergence between populations is shaped by a combination of drift, migration, and selection, yielding patterns of isolation‐by‐distance (IBD) and isolation‐by‐environment (IBE). Unfortunately, IBD and IBE may be confounded when comparing divergence across habitat boundaries. For instance, parapatric lake and stream threespine stickleback (Gasterosteus aculeatus) may have diverged due to selection against migrants (IBE), or mere spatial separation (IBD). To quantitatively partition the strength of IBE and IBD, we used recently developed population genetic software (BEDASSLE) to analyze partial genomic data from three lake‐stream clines on Vancouver Island. We find support for IBD within each of three outlet streams (unlike prior studies of lake‐stream stickleback). In addition, we find evidence for IBE (controlling for geographic distance): the genetic effect of habitat is equivalent to geographic separation of ∼1.9 km of IBD. Remarkably, of our three lake‐stream pairs, IBE is strongest where migration between habitats is easiest. Such microgeographic genetic divergence would require exceptionally strong divergent selection, which multiple experiments have failed to detect. Instead, we propose that nonrandom dispersal (e.g., habitat choice) contributes to IBE. Supporting this conclusion, we show that the few migrants between habitats are a nonrandom subset of the phenotype distribution of the source population.

Resist Globally, Infect Locally: A Transcontinental Test of Adaptation by Stickleback and Their Tapeworm Parasite

Parasite infections are a product of both ecological processes affecting host-parasite encounter rates and evolutionary dynamics affecting host susceptibility. However, few studies examine natural infection variation from both ecological and evolutionary perspectives. Here, we describe the ecological and evolutionary factors generating variation in infection rates by a tapeworm (Schistocephalus solidus) in a vertebrate host, the threespine stickleback (Gasterosteus aculeatus). To explore ecological aspects of infection, we measured tapeworm prevalence in Canadian stickleback inhabiting two distinct environments: marine and freshwater. Consistent with ecological control of infection, the tapeworm is very rare in marine environments, even though marine fish are highly susceptible. Conversely, commonly infected freshwater stickleback exhibit substantial resistance in controlled laboratory trials, suggesting that high exposure risk overwhelms their recently evolved resistance. We also tested for parasite adaptation to its host by performing transcontinental reciprocal infections, using stickleback and tapeworm populations from Europe and western Canada. More infections occurred in same-continent host-parasite combinations, indicating parasite “local” adaptation, at least on the scale of continents. However, the recently evolved immunity of freshwater hosts applies to both local and foreign parasites. The pattern of adaptation described here is not wholly compatible with either of the common models of host-parasite coevolution (i.e., matching infection or targeted recognition). Instead, we propose a hybrid, eco-evolutionary model to explain the remarkable pattern of global host resistance and local parasite infectivity.