A general model of error-prone PCR

In this paper, we generalise a previously-described model of the error-prone polymerase chain reaction (PCR) reaction to conditions of arbitrarily variable amplification efficiency and initial population size. Generalisation of the model to these conditions improves the correspondence to observed and expected behaviours of PCR, and restricts the extent to which the model may explore sequence space for a prescribed set of parameters. Error-prone PCR in realistic reaction conditions is predicted to be less effective at generating grossly divergent sequences than the original model. The estimate of mutation rate per cycle by sampling sequences from an in vitro PCR experiment is correspondingly affected by the choice of model and parameters. (c) 2005 Elsevier Ltd. All rights reserved.

Linking the coevolutionary and population dynamics of host-parasitoid interactions

The interplay between coevolutionary and population or community dynamics is currently the focus of much empirical and theoretical consideration. Here, we develop a simulation model to study the coevolutionary and population dynamics of a hypothetical host-parasitoid interaction. In the model, host resistance and parasitoid virulence are allowed to coevolve. We investigate how trade-offs associated with these traits modify the system's coevolutionary and population dynamics. The most important influence on these dynamics comes from the incorporation of density-dependent costs of resistance ability. We find three main outcomes. First, if the costs of resistance are high, then one or both of the players go extinct. Second, when the costs of resistance are intermediate to low, cycling population and coevolutionary dynamics are found, with slower evolutionary changes observed when the costs of virulence are also low. Third, when the costs associated with resistance and virulence are both high, the hosts trade-off resistance against fecundity and invest little in resistance. However, the parasitoids continue to invest in virulence, leading to stable host and parasitoid population sizes. These results support the hypothesis that costs associated with resistance and virulence will maintain the heritable variation in these traits found in natural populations and that the nature of these trade-offs will greatly influence the population dynamics of the interacting species.

CAMCR: Computer-Assisted Mixture model analysis for Capture–Recapture count data

Population size estimation with discrete or nonparametric mixture models is considered, and reliable ways of construction of the nonparametric mixture model estimator are reviewed and set into perspective. Construction of the maximum likelihood estimator of the mixing distribution is done for any number of components up to the global nonparametric maximum likelihood bound using the EM algorithm. In addition, the estimators of Chao and Zelterman are considered with some generalisations of Zelterman’s estimator. All computations are done with CAMCR, a special software developed for population size estimation with mixture models. Several examples and data sets are discussed and the estimators illustrated. Problems using the mixture model-based estimators are highlighted.

Sensitivity analysis of a cyanobacterial growth and movement model under two different flow regimes

Bloom-forming and toxin-producing cyanobacteria remain a persistent nuisance across the world. Modelling cyanobacterial behaviour in freshwaters is an important tool for understanding their population dynamics and predicting the location and timing of the bloom events in lakes, reservoirs and rivers. A new deterministic–mathematical model was developed, which simulates the growth and movement of cyanobacterial blooms in river systems. The model focuses on the mathematical description of the bloom formation, vertical migration and lateral transport of colonies within river environments by taking into account the major factors that affect the cyanobacterial bloom formation in rivers including light, nutrients and temperature. A parameter sensitivity analysis using a one-at-a-time approach was carried out. There were two objectives of the sensitivity analysis presented in this paper: to identify the key parameters controlling the growth and movement patterns of cyanobacteria and to provide a means for model validation. The result of the analysis suggested that maximum growth rate and day length period were the most significant parameters in determining the population growth and colony depth, respectively.

Limitations of the MRL mouse as a model for cardiac regeneration

Objective Myocardial repair following injury in mammals is restricted such that damaged areas are replaced by scar tissue, impairing cardiac function. MRL mice exhibit exceptional regenerative healing in an ear punch wound model. Some myocardial repair with restoration of heart function has also been reported following cryoinjury. Increased cardiomyocyte proliferation and a foetal liver stem cell population were implicated. We investigated molecular mechanisms facilitating myocardial repair in MRL mice to identify potential therapeutic targets in non-regenerative species. Methods Expressions of specific cell-cycle regulators that might account for regeneration (CDKs 1, 2, 4 and 6; cyclins A, E, D1 and B1; p21, p27 and E2F5) were compared by immunoblotting in MRL and control C57BL/6 ventricles during development. Flow cytometry was used to investigate stem cell populations in livers from foetal mice, and infarct sizes were compared in coronary artery-ligated and sham-treated MRL and C57BL/6 adult mice. Key findings No differences in the expressions of cell cycle regulators were observed between the two strains. Expressions of CD34+Sca1+ckit-, CD34+Sca1+ckit+ and CD34+Sca1-ckit+ increased in livers from C57BL/6 vs MRL mice. No differences were observed in infarct sizes, levels of fibrosis, Ki67 staining or cardiac function between MRL and C57BL/6 mice. Conclusions No intrinsic differences were observed in cell cycle control molecules or stem cell populations between MRL and control C57BL mouse hearts. Pathophysiologically relevant ischaemic injury is not repaired more efficiently in MRL myocardium, questioning the use of the MRL mouse as a reliable model for cardiac regeneration in response to pathophysiologically relevant forms of injury.

How predation and landscape fragmentation affect vole population dynamics

Background: Microtine species in Fennoscandia display a distinct north-south gradient from regular cycles to stable populations. The gradient has often been attributed to changes in the interactions between microtines and their predators. Although the spatial structure of the environment is known to influence predator-prey dynamics of a wide range of species, it has scarcely been considered in relation to the Fennoscandian gradient. Furthermore, the length of microtine breeding season also displays a north-south gradient. However, little consideration has been given to its role in shaping or generating population cycles. Because these factors covary along the gradient it is difficult to distinguish their effects experimentally in the field. The distinction is here attempted using realistic agent-based modelling. Methodology/Principal Findings: By using a spatially explicit computer simulation model based on behavioural and ecological data from the field vole (Microtus agrestis), we generated a number of repeated time series of vole densities whose mean population size and amplitude were measured. Subsequently, these time series were subjected to statistical autoregressive modelling, to investigate the effects on vole population dynamics of making predators more specialised, of altering the breeding season, and increasing the level of habitat fragmentation. We found that fragmentation as well as the presence of specialist predators are necessary for the occurrence of population cycles. Habitat fragmentation and predator assembly jointly determined cycle length and amplitude. Length of vole breeding season had little impact on the oscillations. Significance: There is good agreement between our results and the experimental work from Fennoscandia, but our results allow distinction of causation that is hard to unravel in field experiments. We hope our results will help understand the reasons for cycle gradients observed in other areas. Our results clearly demonstrate the importance of landscape fragmentation for population cycling and we recommend that the degree of fragmentation be more fully considered in future analyses of vole dynamics.

SATgene dietary model to implement diets of differing fat composition in prospectively genotyped groups (apoE) using commercially available foods

esponse to dietary fat manipulation is highly heterogeneous, yet generic population-based recommendations aimed at reducing the burden of CVD are given. The APOE epsilon genotype has been proposed to be an important determinant of this response. The present study reports on the dietary strategy employed in the SATgenɛ (SATurated fat and gene APOE) study, to assess the impact of altered fat content and composition on the blood lipid profile according to the APOE genotype. A flexible dietary exchange model was developed to implement three isoenergetic diets: a low-fat (LF) diet (target composition: 24 % of energy (%E) as fat, 8 %E SFA and 59 %E carbohydrate), a high-saturated fat (HSF) diet (38 %E fat, 18 %E SFA and 45 %E carbohydrate) and a HSF-DHA diet (HSF diet with 3 g DHA/d). Free-living participants (n 88; n 44 E3/E3 and n 44 E3/E4) followed the diets in a sequential design for 8 weeks, each using commercially available spreads, oils and snacks with specific fatty acid profiles. Dietary compositional targets were broadly met with significantly higher total fat (42·8 %E and 41·0 %E v. 25·1 %E, P ≤ 0·0011) and SFA (19·3 %E and 18·6 %E v. 8·33 %E, P ≤ 0·0011) intakes during the HSF and HSF-DHA diets compared with the LF diet, in addition to significantly higher DHA intake during the HSF-DHA diet (P ≤ 0·0011). Plasma phospholipid fatty acid analysis revealed a 2-fold increase in the proportion of DHA after consumption of the HSF-DHA diet for 8 weeks, which was independent of the APOE genotype. In summary, the dietary strategy was successfully implemented in a free-living population resulting in well-tolerated diets which broadly met the dietary targets set.

A multinomial model of applying recognition to judge between multiple alternatives

Proponents of the “fast and frugal” approach to decision-making suggest that inferential judgments are best made on the basis of limited information. For example, if only one of two cities is recognized and the task is to judge which city has the larger population, the recognition heuristic states that the recognized city should be selected. In preference choices with >2 options, it is also standard to assume that a “consideration set”, based upon some simple criterion, is established to reduce the options available. A multinomial processing tree model is outlined which provides the basis for estimating the extent to which recognition is used as a criterion in establishing a consideration set for inferential judgments.

Improvement of 3D protein models using multiple templates guided by single-template model quality assessment

Motivation: Modelling the 3D structures of proteins can often be enhanced if more than one fold template is used during the modelling process. However, in many cases, this may also result in poorer model quality for a given target or alignment method. There is a need for modelling protocols that can both consistently and significantly improve 3D models and provide an indication of when models might not benefit from the use of multiple target-template alignments. Here, we investigate the use of both global and local model quality prediction scores produced by ModFOLDclust2, to improve the selection of target-template alignments for the construction of multiple-template models. Additionally, we evaluate clustering the resulting population of multi- and single-template models for the improvement of our IntFOLD-TS tertiary structure prediction method. Results: We find that using accurate local model quality scores to guide alignment selection is the most consistent way to significantly improve models for each of the sequence to structure alignment methods tested. In addition, using accurate global model quality for re-ranking alignments, prior to selection, further improves the majority of multi-template modelling methods tested. Furthermore, subsequent clustering of the resulting population of multiple-template models significantly improves the quality of selected models compared with the previous version of our tertiary structure prediction method, IntFOLD-TS.

Likelihood-free estimation of model evidence

Statistical methods of inference typically require the likelihood function to be computable in a reasonable amount of time. The class of “likelihood-free” methods termed Approximate Bayesian Computation (ABC) is able to eliminate this requirement, replacing the evaluation of the likelihood with simulation from it. Likelihood-free methods have gained in efficiency and popularity in the past few years, following their integration with Markov Chain Monte Carlo (MCMC) and Sequential Monte Carlo (SMC) in order to better explore the parameter space. They have been applied primarily to estimating the parameters of a given model, but can also be used to compare models. Here we present novel likelihood-free approaches to model comparison, based upon the independent estimation of the evidence of each model under study. Key advantages of these approaches over previous techniques are that they allow the exploitation of MCMC or SMC algorithms for exploring the parameter space, and that they do not require a sampler able to mix between models. We validate the proposed methods using a simple exponential family problem before providing a realistic problem from human population genetics: the comparison of different demographic models based upon genetic data from the Y chromosome.

A framework for adapting survey design through time for wildlife population assessment

Sampling strategies for monitoring the status and trends in wildlife populations are often determined before the first survey is undertaken. However, there may be little information about the distribution of the population and so the sample design may be inefficient. Through time, as data are collected, more information about the distribution of animals in the survey region is obtained but it can be difficult to incorporate this information in the survey design. This paper introduces a framework for monitoring motile wildlife populations within which the design of future surveys can be adapted using data from past surveys whilst ensuring consistency in design-based estimates of status and trends through time. In each survey, part of the sample is selected from the previous survey sample using simple random sampling. The rest is selected with inclusion probability proportional to predicted abundance. Abundance is predicted using a model constructed from previous survey data and covariates for the whole survey region. Unbiased design-based estimators of status and trends and their variances are derived from two-phase sampling theory. Simulations over the short and long-term indicate that in general more precise estimates of status and trends are obtained using this mixed strategy than a strategy in which all of the sample is retained or all selected with probability proportional to predicted abundance. Furthermore the mixed strategy is robust to poor predictions of abundance. Estimates of status are more precise than those obtained from a rotating panel design.

Linking pesticide exposure and spatial dynamics: An individual-based model of wood mouse (Apodemus sylvaticus) populations in agricultural landscapes

The wood mouse is a common and abundant species in agricultural landscape and is a focal species in pesticide risk assessment. Empirical studies on the ecology of the wood mouse have provided sufficient information for the species to be modelled mechanistically. An individual-based model was constructed to explicitly represent the locations and movement patterns of individual mice. This together with the schedule of pesticide application allows prediction of the risk to the population from pesticide exposure. The model included life-history traits of wood mice as well as typical landscape dynamics in agricultural farmland in the UK. The model obtains a good fit to the available population data and is fit for risk assessment purposes. It can help identify spatio-temporal situations with the largest potential risk of exposure and enables extrapolation from individual-level endpoints to population-level effects. Largest risk of exposure to pesticides was found when good crop growth in the “sink” fields coincided with high “source” population densities in the hedgerows. Keywords: Population dynamics, Pesticides, Ecological risk assessment, Habitat choice, Agent-based model, NetLogo

Inferring the biggest and best: a measurement model for applying recognition to evoke consideration sets and judge between multiple alternatives

It has long been supposed that preference judgments between sets of to-be-considered possibilities are made by means of initially winnowing down the most promising-looking alternatives to form smaller “consideration sets” (Howard, 1963; Wright & Barbour, 1977). In preference choices with >2 options, it is standard to assume that a “consideration set”, based upon some simple criterion, is established to reduce the options available. Inferential judgments, in contrast, have more frequently been investigated in situations in which only two possibilities need to be considered (e.g., which of these two cities is the larger?) Proponents of the “fast and frugal” approach to decision-making suggest that such judgments are also made on the basis of limited, simple criteria. For example, if only one of two cities is recognized and the task is to judge which city has the larger population, the recognition heuristic states that the recognized city should be selected. A multinomial processing tree model is outlined which provides the basis for estimating the extent to which recognition is used as a criterion in establishing a consideration set for inferential judgments between three possible options.

Solar ultraviolet-B radiation and vitamin D: a cross-sectional population-based study using data from the 2007 to 2009 Canadian Health Measures Survey

Background: Exposure to solar ultraviolet-B (UV-B) radiation is a major source of vitamin D3. Chemistry climate models project decreases in ground-level solar erythemal UV over the current century. It is unclear what impact this will have on vitamin D status at the population level. The purpose of this study was to measure the association between ground-level solar UV-B and serum concentrations of 25-hydroxyvitamin D (25(OH)D) using a secondary analysis of the 2007 to 2009 Canadian Health Measures Survey (CHMS). Methods: Blood samples collected from individuals aged 12 to 79 years sampled across Canada were analyzed for 25(OH)D (n=4,398). Solar UV-B irradiance was calculated for the 15 CHMS collection sites using the Tropospheric Ultraviolet and Visible Radiation Model. Multivariable linear regression was used to evaluate the association between 25(OH)D and solar UV-B adjusted for other predictors and to explore effect modification. Results: Cumulative solar UV-B irradiance averaged over 91 days (91-day UV-B) prior to blood draw correlated significantly with 25(OH)D. Independent of other predictors, a 1 kJ/m 2 increase in 91-day UV-B was associated with a significant 0.5 nmol/L (95% CI 0.3-0.8) increase in mean 25(OH)D (P =0.0001). The relationship was stronger among younger individuals and those spending more time outdoors. Based on current projections of decreases in ground-level solar UV-B, we predict less than a 1 nmol/L decrease in mean 25(OH)D for the population. Conclusions: In Canada, cumulative exposure to ambient solar UV-B has a small but significant association with 25(OH)D concentrations. Public health messages to improve vitamin D status should target safe sun exposure with sunscreen use, and also enhanced dietary and supplemental intake and maintenance of a healthy body weight.

The determinants of local population growth: a study of Oxfordshire in the nineteenth century

This paper presents a new econometric model for analysing population growth at the village and town level. It develops and applies a theory of the equilibrium distribution of population over space. The theory emphasises geographical fundamentals, such as rivers as transport corridors, and soil types that govern agricultural specialisation; also institutional factors such as town government, market charters and the concentration of land ownership. Nineteenth century Oxfordshire is used as a case study, but the method can also be applied at a multi-county and national level. The results show that the development of railways in nineteenth-century Oxfordshire accelerated a long-term shake-out in the market system, whereby rural markets disappeared and urban markets grew. This shake-out had significant implications for population growth at the local level.