MAXTEMP: A method to maximize precision of the temporal method for estimating Ne in genetic monitoring programs

Abstract

We introduce a new software program, MAXTEMP, that maximizes precision of the temporal method for estimating effective population size (Ne) in genetic monitoring programs, which are increasingly used to systematically track changes in global biodiversity. Scientists and managers are typically most interested in Ne for individual generations, either to match with single-generation estimates of census size (N) or to evaluate consequences of specific management actions or environmental events. Systematically sampling every generation produces a time series of single-generation estimates of temporal [Formula], which can then be used to estimate Ne; however, these estimates have relatively low precision because each reflects just a single episode of genetic drift. Systematic sampling also produces an array of multigenerational temporal estimates that collectively contain a great deal of information about genetic drift that, however, can be difficult to interpret. Here we show how additional information contained in multigenerational temporal estimates can be leveraged to increase precision of [Formula] for individual generations. Using information from one additional generation before and after a target generation can reduce the standard deviation of [Formula] by up to 50%, which not only tightens confidence intervals around [Formula] but also reduces the incidence of extreme estimates. Practical application of MAXTEMP is illustrated with data for a long-term genetic monitoring program for California delta smelt. A second feature of MAXTEMP, which allows one to estimate Ne in an unsampled generation using a combination of temporal and single-sample estimates of Ne from sampled generations, is also described and evaluated.