Capture–recapture analysis with a latent class model allowing for local dependence and observed heterogeneity

In epidemiology, capture–recapture models are commonly used to estimate the size of an unknown population based on several incomplete lists of individuals. The method operates under two main assumptions: independence between the lists (local independence) and homogeneity of capture probabilities of individuals. In practice, these assumptions are rarely satisfied. We introduce a multinomial latent class model that can account for both list dependence and heterogeneity. Parameter estimation is performed by maximizing the conditional likelihood function with the use of the EM algorithm. In addition, a new approach for evaluating the standard errors of the parameter estimates is discussed, which considerably reduces the computational burden associated with the evaluation of the variance of the population size estimate.     

Time to the MRCA of a sample in a Wright-Fisher model with variable population size

Determining the expected distribution of the time to the most recent common ancestor of a sample of individuals may deliver important information about the genetic markers and evolution of the population. In this paper, we introduce a new recursive algorithm to calculate the distribution of the time to the most recent common ancestor of the sample from a population evolved by any conditional multinomial sampling model. The most important advantage of our method is that it can be applied to a sample of any size drawn from a population regardless of its size growth pattern. We also present a very efficient method to implement and store the genealogy tree of the population evolved by the Galton–Watson process. In the final section we present results applied to a simulated population with a single bottleneck event and to real populations of known size histories.     

Transition Densities and Sample Frequency Spectra of Diffusion Processes with Selection and Variable Population Size

Advances in empirical population genetics have made apparent the need for models that simultaneously account for selection and demography. To address this need, we here study the Wright–Fisher diffusion under selection and variable effective population size. In the case of genic selection and piecewise-constant effective population sizes, we obtain the transition density by extending a recently developed method for computing an accurate spectral representation for a constant population size. Utilizing this extension, we show how to compute the sample frequency spectrum in the presence of genic selection and an arbitrary number of instantaneous changes in the effective population size. We also develop an alternate, efficient algorithm for computing the sample frequency spectrum using a moment-based approach. We apply these methods to answer the following questions: If neutrality is incorrectly assumed when there is selection, what effects does it have on demographic parameter estimation? Can the impact of negative selection be observed in populations that undergo strong exponential growth?     

Pest control by the release of insects carrying a female-killing allele on multiple loci

With recent advances in genetics, many new strategies for pest control have become feasible. This is the second article in which we model new techniques for pest control based on the mass release of genetically modified insects. In this article we model the release of insects carrying a dominant and redundant female killing or sterilizing (FK) allele on multiple genetic loci. If such insects are released into a target population, the FK allele can become widely spread in the population through the males while reducing the population each generation by killing females. We allow the number of loci used to vary from 1 to 20. We also allow the FK allele to carry a fitness cost in males due to the gene insertions. Using a model, we explore the effectiveness and optimal strategies for such releases. In the most ideal circumstances (no density-dependence and released insects equal in fitness to wild ones), FK releases are several orders of magnitude more effective than equal sized sterile male releases. For example, a single release of 19 FK-bearing males for every two wild males, with the released males carrying the FK allele on 10 loci, reduces the target population to 0.002% of no-release size. An equal sized sterile release reduces the target population to 5% of no-release size. We also show how the effectiveness of the technique decreases as the fitness cost of the FK alleles in males increases. For example, the above mentioned release reduces the target population to 0.7% of no-release size if each FK allele carries a fitness cost in males of 5%. Adding a simple model for density-dependence and assuming that each of the released males carries the FK allele on six loci, we show that the release size necessary to reduce the target population to 1/100 of no-release size in 10 generations of releases varies from 0.44:1 to 4:1 (depending on parameter values). We also calculate the optimal number of loci on which to put the FK allele under various circumstances.     

A discrete,stochastic model of colonial phytoplankton population size structure: Development and application to in situ imaging-in-flow cytometer observations of Dinobryon

The scientific community lacks models for the dynamic changes in population size structure that occur in colonial phytoplankton. This is surprising, as size is a key trait affecting many aspects of phytoplankton ecology, and colonial forms are very common. We aim to fill this gap with a new discrete, stochastic model of dynamic changes in phytoplankton colonies' population size structure. We use the colonial phytoplankton Dinobryon as a proof-of-concept organism. The model includes four stochastic functions—division, stomatocyst production, colony breakage, and colony loss—to determine Dinobryon population size structure and populations counts. Although the functions presented here are tailored to Dinobryon, the model is readily adaptable to represent other colonial taxa. We demonstrate how fitting our model to in situ observations of colony population size structure can provide a powerful approach to explore colony size dynamics. Here, we have (1) collected high-frequency in situ observations of Dinobryon in Lac (Lake) Montjoie (Quebec, Canada) in 2013 with a moored Imaging FlowCytobot (IFCB) and (2) fit the model to those observations with a genetic algorithm solver that extracts parameter estimates for each of the four stochastic functions. As an example of the power of this model-data integration, we also highlight ecological insights into Dinobryon colony size and stomatocyst production. The Dinobryon population was enriched in larger, flagellate-rich colonies near bloom initiation and shifted to smaller and emptier colonies toward bloom decline.     

The structured ancestral selection graph and the many-demes limit

We show that the unstructured ancestral selection graph applies to part of the history of a sample from a population structured by restricted migration among subpopulations, or demes. The result holds in the limit as the number of demes tends to infinity with proportionately weak selection, and we have also made the assumptions of island-type migration and that demes are equivalent in size. After an instantaneous sample-size adjustment, this structured ancestral selection graph converges to an unstructured ancestral selection graph with a mutation parameter that depends inversely on the migration rate. In contrast, the selection parameter for the population is independent of the migration rate and is identical to the selection parameter in an unstructured population. We show analytically that estimators of the migration rate, based on pairwise sequence differences, derived under the assumption of neutrality should perform equally well in the presence of weak selection. We also modify an algorithm for simulating genealogies conditional on the frequencies of two selected alleles in a sample. This permits efficient simulation of stronger selection than was previously possible. Using this new algorithm, we simulate gene genealogies under the many-demes ancestral selection graph and identify some situations in which migration has a strong effect on the time to the most recent common ancestor of the sample. We find that a similar effect also increases the sensitivity of the genealogy to selection.     

Effects of sample size on the performance of species distribution models

A wide range of modelling algorithms is used by ecologists, conservation practitioners, and others to predict species ranges from point locality data. Unfortunately, the amount of data available is limited for many taxa and regions, making it essential to quantify the sensitivity of these algorithms to sample size. This is the first study to address this need by rigorously evaluating a broad suite of algorithms with independent presence–absence data from multiple species and regions. We evaluated predictions from 12 algorithms for 46 species (from six different regions of the world) at three sample sizes (100, 30, and 10 records). We used data from natural history collections to run the models, and evaluated the quality of model predictions with area under the receiver operating characteristic curve (AUC). With decreasing sample size, model accuracy decreased and variability increased across species and between models. Novel modelling methods that incorporate both interactions between predictor variables and complex response shapes (i.e. GBM, MARS-INT, BRUTO) performed better than most methods at large sample sizes but not at the smallest sample sizes. Other algorithms were much less sensitive to sample size, including an algorithm based on maximum entropy (MAXENT) that had among the best predictive power across all sample sizes. Relative to other algorithms, a distance metric algorithm (DOMAIN) and a genetic algorithm (OM-GARP) had intermediate performance at the largest sample size and among the best performance at the lowest sample size. No algorithm predicted consistently well with small sample size ( n  < 30) and this should encourage highly conservative use of predictions based on small sample size and restrict their use to exploratory modelling.     

Simulating plasticity as a framework for understanding habitat selection and its role in adaptive capacity and extinction risk through an expansion of CDMetaPOP

Adaptive capacity can present challenges for modelling as it encompasses multiple ecological and evolutionary processes such as natural selection, genetic drift, gene flow and phenotypic plasticity. Spatially explicit, individual-based models provide an outlet for simulating these complex interacting eco-evolutionary processes. We expanded the existing Cost-Distance Meta-POPulation (CDMetaPOP) framework with inducible plasticity modelled as a habitat selection behaviour, using temperature or habitat quality variables, with a genetically based selection threshold conditioned on past individual experience. To demonstrate expected results in the new module, we simulated hypothetical populations and then evaluated model performance in populations of redband trout (Oncorhynchus mykiss gairdneri) across three watersheds where temperatures induce physiological stress in parts of the stream network. We ran simulations using projected warming stream temperature data under four scenarios for alleles that: (1) confer thermal tolerance, (2) bestow plastic habitat selection, (3) give both thermal tolerance and habitat selection preference and (4) do not provide either thermal tolerance or habitat selection. Inclusion of an adaptive allele decreased declines in population sizes, but this impact was greatly reduced in the relatively cool stream networks. As anticipated with the new module, high-temperature patches remained unoccupied by individuals with the allele operating plastically after exposure to warm temperatures. Using complete habitat avoidance above the stressful temperature threshold, habitat selection reduced the overall population size due to the opportunity cost of avoiding areas with increased, but not guaranteed, mortality. Inclusion of plasticity within CDMetaPOP will provide the potential for genetic or plastic traits and ‘rescue’ to affect eco-evolutionary dynamics for research questions and conservation applications.     

Dynamic population artificial bee colony algorithm for multi-objective optimal power flow

This paper proposes a novel artificial bee colony algorithm with dynamic population (ABC-DP), which synergizes the idea of extended life-cycle evolving model to balance the exploration and exploitation tradeoff. The proposed ABC-DP is a more bee-colony-realistic model that the bee can reproduce and die dynamically throughout the foraging process and population size varies as the algorithm runs. ABC-DP is then used for solving the optimal power flow (OPF) problem in power systems that considers the cost, loss, and emission impacts as the objective functions. The 30-bus IEEE test system is presented to illustrate the application of the proposed algorithm. The simulation results, which are also compared to nondominated sorting genetic algorithm II (NSGAII) and multi-objective ABC (MOABC), are presented to illustrate the effectiveness and robustness of the proposed method.     

Effects of size‐ and sex‐selective harvesting: An integral projection model approach

Harvesting is often size‐selective, and in species with sexual size dimorphism, it may also be sex‐selective. A powerful approach to investigate potential consequences of size‐ and/or sex‐selective harvesting is to simulate it in a demographic population model. We developed a population‐based integral projection model for a size‐ and sex‐structured species, the commonly exploited pike (Esox lucius). The model allows reproductive success to be proportional to body size and potentially limited by both sexes. We ran all harvest simulations with both lower size limits and slot limits, and to quantify the effects of selective harvesting, we calculated sex ratios and the long‐term population growth rate (λ). In addition, we quantified to what degree purely size‐selective harvesting was sex‐selective, and determined when λ shifted from being female to male limited under size‐ and sex‐selective harvesting. We found that purely size‐selective harvest can be sex‐selective, and that it depends on the harvest limits and the size distributions of the sexes. For the size‐ and sex‐selective harvest simulations, λ increased with harvest intensity up to a threshold as females limited reproduction. Beyond this threshold, males became the limiting sex, and λ decreased as more males were harvested. The peak in λ, and the corresponding sex ratio in harvest, varied with both the selectivity and the intensity of the harvest simulation. Our model represents a useful extension of size‐structured population models as it includes both sexes, relaxes the assumption of female dominance, and accounts for size‐dependent fecundity. The consequences of selective harvesting presented here are especially relevant for size‐ and sex‐structured exploited species, such as commercial fisheries. Thus, our model provides a useful contribution toward the development of more sustainable harvesting regimes.     

Differential-Evolution Control Parameter Optimization for Unmanned Aerial Vehicle Path Planning

The differential evolution algorithm has been widely applied on unmanned aerial vehicle (UAV) path planning. At present, four random tuning parameters exist for differential evolution algorithm, namely, population size, differential weight, crossover, and generation number. These tuning parameters are required, together with user setting on path and computational cost weightage. However, the optimum settings of these tuning parameters vary according to application. Instead of trial and error, this paper presents an optimization method of differential evolution algorithm for tuning the parameters of UAV path planning. The parameters that this research focuses on are population size, differential weight, crossover, and generation number. The developed algorithm enables the user to simply define the weightage desired between the path and computational cost to converge with the minimum generation required based on user requirement. In conclusion, the proposed optimization of tuning parameters in differential evolution algorithm for UAV path planning expedites and improves the final output path and computational cost.     

Comparison of photo‐matching algorithms commonly used for photographic capture–recapture studies

Photographic capture–recapture is a valuable tool for obtaining demographic information on wildlife populations due to its noninvasive nature and cost‐effectiveness. Recently, several computer‐aided photo‐matching algorithms have been developed to more efficiently match images of unique individuals in databases with thousands of images. However, the identification accuracy of these algorithms can severely bias estimates of vital rates and population size. Therefore, it is important to understand the performance and limitations of state‐of‐the‐art photo‐matching algorithms prior to implementation in capture–recapture studies involving possibly thousands of images. Here, we compared the performance of four photo‐matching algorithms; Wild‐ID, I3S Pattern+, APHIS, and AmphIdent using multiple amphibian databases of varying image quality. We measured the performance of each algorithm and evaluated the performance in relation to database size and the number of matching images in the database. We found that algorithm performance differed greatly by algorithm and image database, with recognition rates ranging from 100% to 22.6% when limiting the review to the 10 highest ranking images. We found that recognition rate degraded marginally with increased database size and could be improved considerably with a higher number of matching images in the database. In our study, the pixel‐based algorithm of AmphIdent exhibited superior recognition rates compared to the other approaches. We recommend carefully evaluating algorithm performance prior to using it to match a complete database. By choosing a suitable matching algorithm, databases of sizes that are unfeasible to match “by eye” can be easily translated to accurate individual capture histories necessary for robust demographic estimates.     

Increasing sequence correlation limits the efficiency of recombination in a multisite evolution model

The accumulation of preexisting beneficial alleles in a haploid population, under selection and infrequent recombination, and in the absence of new mutation events is studied numerically by means of detailed Monte Carlo simulations. On the one hand, we confirm our previous work, in that the accumulation rate follows modified single-site kinetics, with a timescale set by an effective selection coefficient s(eff) as shown in a previous work, and we confirm the qualitative features of the dependence of s(eff) on the population size and the recombination rate reported therein. In particular, we confirm the existence of a threshold population size below which evolution stops before the emergence of best-fit individuals. On the other hand, our simulations reveal that the population dynamics is essentially shaped by the steady accumulation of pairwise sequence correlation, causing sequence congruence in excess of what one would expect from a uniformly random distribution of alleles. By sequence congruence, we understand here the opposite of genetic distance, that is, the fraction of monomorphic sites of specified allele type in a pair of genomes (individual sequences). The effective selection coefficient changes more rapidly with the recombination rate and has a higher threshold in this parameter than found in the previous work, which neglected correlation effects. We examine this phenomenon by monitoring the time dependence of sequence correlation based on a set of sequence congruence measures and verify that it is not associated with the development of linkage disequilibrium. We also discuss applications to HIV evolution in infected individuals and potential implications for drug therapy.     

Modelling horse management in the Australian Alps

Management of the feral Horse (Equus caballus) in the Australian Alps bioregion is a difficult and emotive issue, with interested parties working from vastly differing perspectives. Compounding this, information regarding ecology and distribution of horses, and the cost and effectiveness of management strategies is often unknown or uncertain. Resolving these issues requires an objective approach with the flexibility to incorporate different potential scenarios. We used a spatially explicit population model to compare the potential effects of two different management strategies on populations of horses in the Australian Alps bioregion: culling from helicopters versus trapping and mustering. We populated the model using the results of population surveys conducted in 2014, vegetation data and cost estimates. We then provided an estimate of the effect of each strategy on population size across the Alps, and their corresponding costs, compared to no management. To account for uncertainties, we simulated different scenarios for horse population densities, dispersal rates and population growth rates. Management using aerial culling was more effective than mustering in every scenario modelled, and three to six times cheaper. Aerial culling was only slightly more effective within its control region. However, because mustering is necessarily restricted by road access, this translated to a substantial improvement in population control – up to 2000 horses where growth and dispersal rates were high. Our results unequivocally suggest aerial culling as the only strategy that could effectively control horses within the modelled range of scenarios; this result stands in addition to its other potential benefits of lower cost, animal stress and landscape disturbance. A major advantage of this modelling approach is that we can easily update it with new data, test different measures of effectiveness and add new scenarios to adapt to the rapidly changing situation on the ground, both in terms of the ecology and the political climate.     

Proceedings of the SMBE Tri-National Young Investigators' Workshop 2005. Mutation rate and the cost of complexity

Two recent theoretical studies of adaptation suggest that more complex organisms tend to adapt more slowly. Specifically, in Fisher's "geometric" model of a finite population where multiple traits are under optimizing selection, the average progress ensuing from a single mutation decreases as the number of traits increases--the "cost of complexity." Here, I draw on molecular and histological data to assess the extent to which on a large phylogenetic scale, this predicted decrease in the rate of adaptation per mutation is mitigated by an increase in the number of mutations per generation as complexity increases. As an index of complexity for multicellular organisms, I use the number of visibly distinct types of cell in the body. Mutation rate is the product of mutational target size and population mutation rate per unit target. Despite much scatter, genome size appears to be positively correlated with complexity (as indexed by cell-type number), which along with other considerations suggests that mutational target size tends to increase with complexity. In contrast, effective population mutation rate per unit target appears to be negatively correlated with complexity. The net result is that mutation rate probably does tend to increase with complexity, although probably not fast enough to eliminate the cost of complexity.     

Inferring landscape resistance to gene flow when genetic drift is spatially heterogeneous

In connectivity models, land cover types are assigned cost values characterizing their resistance to species movements. Landscape genetic methods infer these values from the relationship between genetic differentiation and cost distances. The spatial heterogeneity of population sizes, and consequently genetic drift, is rarely included in this inference although it influences genetic differentiation. Similarly, migration rates and population spatial distributions potentially influence this inference. Here, we assessed the reliability of cost value inference under several migration rates, population spatial patterns and degrees of population size heterogeneity. Additionally, we assessed whether considering intra-population variables, here using gravity models, improved the inference when drift is spatially heterogeneous. We simulated several gene flow intensities between populations with varying local sizes and spatial distributions. We then fit gravity models of genetic distances as a function of (i) the ‘true’ cost distances driving simulations or alternative cost distances, and (ii) intra-population variables (population sizes, patch areas). We determined the conditions making the identification of the ‘true’ costs possible and assessed the contribution of intra-population variables to this objective. Overall, the inference ranked cost scenarios reliably in terms of similarity with the ‘true’ scenario (cost distance Mantel correlations), but this ‘true’ scenario rarely provided the best model goodness of fit. Ranking inaccuracies and failures to identify the ‘true’ scenario were more pronounced when migration was very restricted (<4 dispersal events/generation), population sizes were most heterogeneous and some populations were spatially aggregated. In these situations, considering intra-population variables helps identify cost scenarios reliably, thereby improving cost value inference from genetic data.     

RATES OF FITNESS DECLINE AND REBOUND SUGGEST PERVASIVE EPISTASIS

Unraveling the factors that determine the rate of adaptation is a major question in evolutionary biology. One key parameter is the effect of a new mutation on fitness, which invariably depends on the environment and genetic background. The fate of a mutation also depends on population size, which determines the amount of drift it will experience. Here, we manipulate both population size and genotype composition and follow adaptation of 23 distinct Escherichia coli genotypes. These have previously accumulated mutations under intense genetic drift and encompass a substantial fitness variation. A simple rule is uncovered: the net fitness change is negatively correlated with the fitness of the genotype in which new mutations appear—a signature of epistasis. We find that Fisher's geometrical model can account for the observed patterns of fitness change and infer the parameters of this model that best fit the data, using Approximate Bayesian Computation. We estimate a genomic mutation rate of 0.01 per generation for fitness altering mutations, albeit with a large confidence interval, a mean fitness effect of mutations of ?0.01, and an effective number of traits nine in mutS^? E. coli. This framework can be extended to confront a broader range of models with data and test different classes of fitness landscape models.     

Population-Based Reversible Jump Markov Chain Monte Carlo

We present an extension of population-based Markov chain MonteCarlo to the transdimensional case. A major challenge is thatof simulating from high- and transdimensional target measures.In such cases, Markov chain Monte Carlo methods may not adequatelytraverse the support of the target; the simulation results willbe unreliable. We develop population methods to deal with suchproblems, and give a result proving the uniform ergodicity ofthese population algorithms, under mild assumptions. This resultis used to demonstrate the superiority, in terms of convergencerate, of a population transition kernel over a reversible jumpsampler for a Bayesian variable selection problem. We also givean example of a population algorithm for a Bayesian multivariatemixture model with an unknown number of components. This isapplied to gene expression data of 1000 data points in six dimensionsand it is demonstrated that our algorithm outperforms some competingMarkov chain samplers. In this example, we show how to combinethe methods of parallel chains (Geyer, 1991), tempering (Geyer& Thompson, 1995), snooker algorithms (Gilks et al., 1994),constrained sampling and delayed rejection (Green & Mira,2001).     

Time to the MRCA of a sample in a Wright–Fisher model with variable population size

Determining the expected distribution of the time to the most recent common ancestor of a sample of individuals may deliver important information about the genetic markers and evolution of the population. In this paper, we introduce a new recursive algorithm to calculate the distribution of the time to the most recent common ancestor of the sample from a population evolved by any conditional multinomial sampling model. The most important advantage of our method is that it can be applied to a sample of any size drawn from a population regardless of its size growth pattern. We also present a very efficient method to implement and store the genealogy tree of the population evolved by the Galton–Watson process. In the final section we present results applied to a simulated population with a single bottleneck event and to real populations of known size histories.     

Presence-only data and the em algorithm

Summary .  In ecological modeling of the habitat of a species, it can be prohibitively expensive to determine species absence. Presence-only data consist of a sample of locations with observed presences and a separate group of locations sampled from the full landscape, with unknown presences. We propose an expectation–maximization algorithm to estimate the underlying presence–absence logistic model for presence-only data. This algorithm can be used with any off-the-shelf logistic model. For models with stepwise fitting procedures, such as boosted trees, the fitting process can be accelerated by interleaving expectation steps within the procedure. Preliminary analyses based on sampling from presence–absence records of fish in New Zealand rivers illustrate that this new procedure can reduce both deviance and the shrinkage of marginal effect estimates that occur in the naive model often used in practice. Finally, it is shown that the population prevalence of a species is only identifiable when there is some unrealistic constraint on the structure of the logistic model. In practice, it is strongly recommended that an estimate of population prevalence be provided.