AlleleShift: an R package to predict and visualize population-level changes in allele frequencies in response to climate change

At any particular location, frequencies of alleles that are associated with adaptive traits are expected to change in future climates through local adaption and migration, including assisted migration (human-implemented when climate change is more rapid than natural migration rates). Making the assumption that the baseline frequencies of alleles across environmental gradients can act as a predictor of patterns in changed climates (typically future but possibly paleo-climates), a methodology is provided by AlleleShift of predicting changes in allele frequencies at the population level.

Methods

The prediction procedure involves a first calibration and prediction step through redundancy analysis (RDA), and a second calibration and prediction step through a generalized additive model (GAM) with a binomial family. As such, the procedure is fundamentally different to an alternative approach recently proposed to predict changes in allele frequencies from canonical correspondence analysis (CCA). The RDA step is based on the Euclidean distance that is also the typical distance used in Analysis of Molecular Variance (AMOVA). Because the RDA step or CCA approach sometimes predict negative allele frequencies, the GAM step ensures that allele frequencies are in the range of 0 to 1.

Results

AlleleShift provides data sets with predicted frequencies and several visualization methods to depict the predicted shifts in allele frequencies from baseline to changed climates. These visualizations include ‘dot plot’ graphics (function shift.dot.ggplot ), pie diagrams ( shift.pie.ggplot ), moon diagrams ( shift.moon.ggplot ), ‘waffle’ diagrams ( shift.waffle.ggplot ) and smoothed surface diagrams of allele frequencies of baseline or future patterns in geographical space ( shift.surf.ggplot ). As these visualizations were generated through the ggplot2 package, methods of generating animations for a climate change time series are straightforward, as shown in the documentation of AlleleShift and in the supplemental videos.

Availability

AlleleShift is available as an open-source R package from https://cran.r-project.org/package=AlleleShift and https://github.com/RoelandKindt/AlleleShift . Genetic input data is expected to be in the adegenet::genpop format, which can be generated from the adegenet::genind format. Climate data is available from various resources such as WorldClim and Envirem .