Genetic assignment of individuals to source populations using network estimation tools

Abstract

<div lang="en" class="abs"><p class="abs_header">Abstract</p><p><em>1.</em> Dispersal, the movement of individuals between populations, is crucial in many ecological and genetic processes. However, direct identification of dispersing individuals is difficult or impossible in natural populations. By using genetic assignment methods, individuals with unknown genetic origin can be assigned to source populations. This knowledge is necessary in studying many key questions in ecology, evolution and conservation.</p><p><em>2.</em> We introduce a network-based tool BONE (Baseline Oriented Network Estimation) for genetic population assignment, which borrows concepts from undirected graph inference. In particular, we use sparse multinomial Least Absolute Shrinkage and Selection Operator (LASSO) regression to estimate probability of the origin of all mixture individuals and their mixture proportions without tedious selection of the LASSO tuning parameter. We compare BONE with three genetic assignment methods implemented in R packages radmixture, assignPOP and RUBIAS.</p><p><em>3.</em> Probability of the origin and mixture proportion estimates of both simulated and real data (an insular house sparrow metapopulation and Chinook salmon populations) given by BONE are competitive or superior compared to other assignment methods. Our examples illustrate how the network estimation method adapts to population assignment, combining the efficiency and attractive properties of sparse network representation and model selection properties of the <em>L</em>₁ regularization. As far as we know, this is the first approach showing how one can use network tools for genetic identification of individuals’ source populations.</p><p><em>4.</em> BONE is aimed at any researcher performing genetic assignment and trying to infer the genetic population structure. Compared to other methods, our approach also identifies outlying mixture individuals that could originate outside of the baseline populations. BONE is a freely available R package under the GPL licence and can be downloaded at GitHub. In addition to the R package, a tutorial for BONE is available at <a href="https://github.com/markkukuismin/BONE/">https://github.com/markkukuismin/BONE/</a>.</p></div>