Investigating the effects of predator removal and habitat management on nest success and breeding population size of a farmland passerine: a case study

Passerines are especially vulnerable to predation at the pre-independence stage. Although the role of nest success in British farmland passerine declines is contentious, improvement in nest success through sympathetic management could play a role in their reversal. Because habitat is known to interact with predation, management options for mitigation will need to consider effects of nest predation. We present results from an observational study of a population of Common Blackbird Turdus merula on a farm which has experienced a range of agri-environment and game-management options, including a period with nest predator control, as a case study to address some of these issues. We used an information theoretic model comparison procedure to look for evidence of interactions between habitat and nest predation, and then asked whether habitat management and nest predator abundances could explain population trends at the site through their effects on nest success. Interactions were detected between measures of predator abundance and habitat variables, and these varied with nest stage - habitat within the vicinity of the nest appeared to be important at the egg stage, and nest-placement characteristics were important at the nestling stage. Although predator control appeared to have a positive influence on Blackbird breeding population size, the non-experimental set-up meant we could not eliminate other potential explanations. Variation in breeding population size did not appear to be influenced by variation in nest success alone. Our study demonstrates that observational data can only go so far in detection of such effects, and we discuss how it might be taken further. Agri-environment and game-management techniques are likely to influence nest predation pressure on farmland passerines, but the patterns, mechanisms and importance to population processes remain not wholly understood.

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.

Estimating the hidden number of scrapie affected holdings in great Britain using a simple, truncated count model allowing for heterogeneity

None of the current surveillance streams monitoring the presence of scrapie in Great Britain provide a comprehensive and unbiased estimate of the prevalence of the disease at the holding level. Previous work to estimate the under-ascertainment adjusted prevalence of scrapie in Great Britain applied multiple-list capture-recapture methods. The enforcement of new control measures on scrapie-affected holdings in 2004 has stopped the overlapping between surveillance sources and, hence, the application of multiple-list capture-recapture models. Alternative methods, still under the capture-recapture methodology, relying on repeated entries in one single list have been suggested in these situations. In this article, we apply one-list capture-recapture approaches to data held on the Scrapie Notifications Database to estimate the undetected population of scrapie-affected holdings with clinical disease in Great Britain for the years 2002, 2003, and 2004. For doing so, we develop a new diagnostic tool for indication of heterogeneity as well as a new understanding of the Zelterman and Chao's lower bound estimators to account for potential unobserved heterogeneity. We demonstrate that the Zelterman estimator can be viewed as a maximum likelihood estimator for a special, locally truncated Poisson likelihood equivalent to a binomial likelihood. This understanding allows the extension of the Zelterman approach by means of logistic regression to include observed heterogeneity in the form of covariates-in case studied here, the holding size and country of origin. Our results confirm the presence of substantial unobserved heterogeneity supporting the application of our two estimators. The total scrapie-affected holding population in Great Britain is around 300 holdings per year. None of the covariates appear to inform the model significantly.

A simple variance formula for population size estimators by conditioning

This note considers the variance estimation for population size estimators based on capture–recapture experiments. Whereas a diversity of estimators of the population size has been suggested, the question of estimating the associated variances is less frequently addressed. This note points out that the technique of conditioning can be applied here successfully which also allows us to identify sources of variation: the variance due to estimation of the model parameters and the binomial variance due to sampling n units from a population of size N. It is applied to estimators typically used in capture–recapture experiments in continuous time including the estimators of Zelterman and Chao and improves upon previously used variance estimators. In addition, knowledge of the variances associated with the estimators by Zelterman and Chao allows the suggestion of a new estimator as the weighted sum of the two. The decomposition of the variance into the two sources allows also a new understanding of how resampling techniques like the Bootstrap could be used appropriately. Finally, the sample size question for capture–recapture experiments is addressed. Since the variance of population size estimators increases with the sample size, it is suggested to use relative measures such as the observed-to-hidden ratio or the completeness of identification proportion for approaching the question of sample size choice.

The extrapolation problem and how population modeling can help

We argue that population modeling can add value to ecological risk assessment by reducing uncertainty when extrapolating from ecotoxicological observations to relevant ecological effects. We review other methods of extrapolation, ranging from application factors to species sensitivity distributions to suborganismal (biomarker and "-omics'') responses to quantitative structure activity relationships and model ecosystems, drawing attention to the limitations of each. We suggest a simple classification of population models and critically examine each model in an extrapolation context. We conclude that population models have the potential for adding value to ecological risk assessment by incorporating better understanding of the links between individual responses and population size and structure and by incorporating greater levels of ecological complexity. A number of issues, however, need to be addressed before such models are likely to become more widely used. In a science context, these involve challenges in parameterization, questions about appropriate levels of complexity, issues concerning how specific or general the models need to be, and the extent to which interactions through competition and trophic relationships can be easily incorporated.

Extending Zelterman’s approach for robust estimation of population size to zero-truncated clustered data

Estimation of population size with missing zero-class is an important problem that is encountered in epidemiological assessment studies. Fitting a Poisson model to the observed data by the method of maximum likelihood and estimation of the population size based on this fit is an approach that has been widely used for this purpose. In practice, however, the Poisson assumption is seldom satisfied. Zelterman (1988) has proposed a robust estimator for unclustered data that works well in a wide class of distributions applicable for count data. In the work presented here, we extend this estimator to clustered data. The estimator requires fitting a zero-truncated homogeneous Poisson model by maximum likelihood and thereby using a Horvitz-Thompson estimator of population size. This was found to work well, when the data follow the hypothesized homogeneous Poisson model. However, when the true distribution deviates from the hypothesized model, the population size was found to be underestimated. In the search of a more robust estimator, we focused on three models that use all clusters with exactly one case, those clusters with exactly two cases and those with exactly three cases to estimate the probability of the zero-class and thereby use data collected on all the clusters in the Horvitz-Thompson estimator of population size. Loss in efficiency associated with gain in robustness was examined based on a simulation study. As a trade-off between gain in robustness and loss in efficiency, the model that uses data collected on clusters with at most three cases to estimate the probability of the zero-class was found to be preferred in general. In applications, we recommend obtaining estimates from all three models and making a choice considering the estimates from the three models, robustness and the loss in efficiency. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)

Population genetic models can be used to study the evolution of the interacting behaviours of parents and their progeny

An example of the evolution of the interacting behaviours of parents and progeny is studied using iterative equations linking the frequencies of the gametes produced by the progeny to the frequencies of the gametes in the parental generation. This population genetics approach shows that a model in which both behaviours are determined by a single locus can lead to a stable equilibrium in which the two behaviours continue to segregate. A model in which the behaviours are determined by genes at two separate loci leads eventually to fixation of the alleles at both loci but this can take many generations of selection. Models of the type described in this paper will be needed to understand the evolution of complex behaviour when genomic or experimental information is available about the genetic determinants of behaviour and the selective values of different genomes. (c) 2007 Elsevier Inc. All rights reserved.

Chronic toxicity of ibuprofen to Daphnia magna: Effects on life history traits and population dynamics

The non-steroidal anti-inflammatory drug (NSAID) ibuprofen (IB) is a widely used pharmaceutical that can be found in several freshwater ecosystems. Acute toxicity studies with Daphnia magna suggest that the 48 h EC50 (immobilisation) is 10-100 mg IB l(-1). However, there are currently no chronic IB toxicity dataon arthropod populations, and the aquatic life impacts of such analgesic drugs are still undefined. We performed a 14-day exposure of D. magna to IB as a model compound (concentration range: 0, 20, 40 and 80 mg IB l(-1)) measuring chronic effects on life history traits and population performance. Population growth rate was significantly reduced at all IB concentrations, although survival was only affected at 80 mg IB l(-1). Reproduction, however, was affected at lower concentrations of IB (14-day EC50 of 13.4 mg IB l(-1)), and was completely inhibited at the highest test concentration. The results from this study indicate that the long-term crustacean population consequences of a chronic IB exposure at environmentally realistic concentrations (ng l(-1) to mu g l(-1)) would most likely be of minor importance. We discuss our results in relation to recent genomic studies, which suggest that the potential mechanism of toxicity in Daphnia is similar to the mode of action in mammals, where IB inhibits eicosanoid biosynthesis. (C) 2007 Elsevier Ireland Ltd. All rights reserved.

Population growth rate and carrying capacity for springtails Folsomia candida exposed to ivermectin

Forecasting the effects of stressors on the dynamics of natural populations requires assessment of the joint effects of a stressor and population density on the population response. The effects can be depicted as a contour map in which the population response, here assessed by Population growth rate, varies with stress and density in the same way that the height of land above sea level varies with latitude and longitude. We present the first complete map of this type using as our model Folsomia candida exposed to five different concentrations of the widespread anthelmintic veterinary medicine ivermectin in replicated microcosm experiments lasting 49 days. The concentrations of ivermectin in yeast were 0.0, 6.8 28.83 66.4 and 210.0 mg/L wet weight. Increasing density and chemical concentration both significantly reduced the population growth rate of Folsomia candida, in part through effects on food consumption and fecundity. The interaction between density and ivermectin concentration was "less-than-additive," implying that at high density populations were able to compensate for the effects of the chemical. This result demonstrates that regulatory protocols carried out at low density (as in most past experiments) may seriously overestimate effects in the field, where densities are locally high and populations are resource limited (e.g., in feces of livestock treated with ivermectin).

Population-level assessment of risks of pesticides to birds and mammals in the UK

It is generally acknowledged that population-level assessments provide,I better measure of response to toxicants than assessments of individual-level effects. population-level assessments generally require the use of models to integrate potentially complex data about the effects of toxicants on life-history traits, and to provide a relevant measure of ecological impact. Building on excellent earlier reviews we here briefly outline the modelling options in population-level risk assessment. Modelling is used to calculate population endpoints from available data, which is often about Individual life histories, the ways that individuals interact with each other, the environment and other species, and the ways individuals are affected by pesticides. As population endpoints, we recommend the use of population abundance, population growth rate, and the chance of population persistence. We recommend two types of model: simple life-history models distinguishing two life-history stages, juveniles and adults; and spatially-explicit individual-based landscape models. Life-history models are very quick to set up and run, and they provide a great deal or insight. At the other extreme, individual-based landscape models provide the greatest verisimilitude, albeit at the cost of greatly increased complexity. We conclude with a discussion of the cations of the severe problems of parameterising models.

Estimation of population growth or decline in genetically monitored populations

This article introduces a new general method for genealogical inference that samples independent genealogical histories using importance sampling (IS) and then samples other parameters with Markov chain Monte Carlo (MCMC). It is then possible to more easily utilize the advantages of importance sampling in a fully Bayesian framework. The method is applied to the problem of estimating recent changes in effective population size from temporally spaced gene frequency data. The method gives the posterior distribution of effective population size at the time of the oldest sample and at the time of the most recent sample, assuming a model of exponential growth or decline during the interval. The effect of changes in number of alleles, number of loci, and sample size on the accuracy of the method is described using test simulations, and it is concluded that these have an approximately equivalent effect. The method is used on three example data sets and problems in interpreting the posterior densities are highlighted and discussed.

Declining survival rates in a reintroduced population of the Mauritius kestrel: evidence for non-linear density dependence and environmental stochasticity

1. We studied a reintroduced population of the formerly critically endangered Mauritius kestrel Falco punctatus Temmink from its inception in 1987 until 2002, by which time the population had attained carrying capacity for the study area. Post-1994 the population received minimal management other than the provision of nestboxes. 2. We analysed data collected on survival (1987-2002) using program MARK to explore the influence of density-dependent and independent processes on survival over the course of the population's development. 3.We found evidence for non-linear, threshold density dependence in juvenile survival rates. Juvenile survival was also strongly influenced by climate, with the temporal distribution of rainfall during the cyclone season being the most influential climatic variable. Adult survival remained constant throughout. 4. Our most parsimonious capture-mark-recapture statistical model, which was constrained by density and climate, explained 75.4% of the temporal variation exhibited in juvenile survival rates over the course of the population's development. 5. This study is an example of how data collected as part of a threatened species recovery programme can be used to explore the role and functional form of natural population regulatory processes. With the improvements in conservation management techniques and the resulting success stories, formerly threatened species offer unique opportunities to further our understanding of the fundamental principles of population ecology.

Population subdivision and molecular sequence variation: Theory and analysis of Drosophila ananassae data

Population subdivision complicates analysis of molecular variation. Even if neutrality is assumed, three evolutionary forces need to be considered: migration, mutation, and drift. Simplification can be achieved by assuming that the process of migration among and drift within subpopulations is occurring fast compared to Mutation and drift in the entire population. This allows a two-step approach in the analysis: (i) analysis of population subdivision and (ii) analysis of molecular variation in the migrant pool. We model population subdivision using an infinite island model, where we allow the migration/drift parameter Theta to vary among populations. Thus, central and peripheral populations can be differentiated. For inference of Theta, we use a coalescence approach, implemented via a Markov chain Monte Carlo (MCMC) integration method that allows estimation of allele frequencies in the migrant pool. The second step of this approach (analysis of molecular variation in the migrant pool) uses the estimated allele frequencies in the migrant pool for the study of molecular variation. We apply this method to a Drosophila ananassae sequence data set. We find little indication of isolation by distance, but large differences in the migration parameter among populations. The population as a whole seems to be expanding. A population from Bogor (Java, Indonesia) shows the highest variation and seems closest to the species center.

Effect of fecal water on an in vitro model of colonic mucosal barrier function

Fecal water (FW) has been shown to exert, in cultured cells, cytotoxic and genotoxic effects that have implications for colorectal cancer (CRC) risk. We have investigated a further biological activity of FW, namely, the ability to affect gap junctions in CACO2 cell monolayers as an index of mucosal barrier function, which is known to be disrupted in cancer. FW samples fi-om healthy, free-living, European subjects that were divided into two broad age groups, adult (40 +/- 9.7 yr; n = 53) and elderly (76 +/- 7.5 yr; n = 55) were tested for effects on gap junction using the transepithelial resistance (TER) assay. Overall, treatment of CACO2 cells with FW samples fi-om adults increased TER (+ 4 %), whereas FW from elderly subjects decreased TER (-5%); the difference between the two groups was significant (P < 0.05). We also measured several components of FW potentially associated with modulation of TER, namely, short-chain fatty acid (SCFA) and ammonia. SCFAs (propionic, acetic, and n-butyric) were significantly lower in the elderly population (-30%, -35%, and -21%, respectively, all P pound 0.01). We consider that FW modulation of in vitro epithelial barrier function is a potentially useful noninvasive biomarker, but it requires further validation to establish its relationship to CRC risk.

Revised QUICKI provides a strong surrogate estimate of insulin sensitivity when compared with the minimal model

OBJECTIVE: To compare insulin sensitivity (Si) from a frequently sampled intravenous glucose tolerance test (FSIGT) and subsequent minimal model analyses with surrogate measures of insulin sensitivity and resistance and to compare features of the metabolic syndrome between Caucasians and Indian Asians living in the UK. SUBJECTS: In all, 27 healthy male volunteers (14 UK Caucasians and 13 UK Indian Asians), with a mean age of 51.2 +/- 1.5 y, BMI of 25.8 +/- 0.6 kg/m(2) and Si of 2.85 +/- 0.37. MEASUREMENTS: Si was determined from an FSIGT with subsequent minimal model analysis. The concentrations of insulin, glucose and nonesterified fatty acids (NEFA) were analysed in fasting plasma and used to calculate surrogate measure of insulin sensitivity (quantitative insulin sensitivity check index (QUICKI), revised QUICKI) and resistance (homeostasis for insulin resistance (HOMA IR), fasting insulin resistance index (FIRI), Bennetts index, fasting insulin, insulin-to-glucose ratio). Plasma concentrations of triacylglycerol (TAG), total cholesterol, high density cholesterol, (HDL-C) and low density cholesterol, (LDL-C) were also measured in the fasted state. Anthropometric measurements were conducted to determine body-fat distribution. RESULTS: Correlation analysis identified the strongest relationship between Si and the revised QUICKI (r = 0.67; P = 0.000). Significant associations were also observed between Si and QUICKI (r = 0.51; P = 0.007), HOMA IR (r = -0.50; P = 0.009), FIRI and fasting insulin. The Indian Asian group had lower HDL-C (P = 0.001), a higher waist-hip ratio (P = 0.01) and were significantly less insulin sensitive (Si) than the Caucasian group (P = 0.02). CONCLUSION: The revised QUICKI demonstrated a statistically strong relationship with the minimal model. However, it was unable to differentiate between insulin-sensitive and -resistant groups in this study. Future larger studies in population groups with varying degrees of insulin sensitivity are recommended to investigate the general applicability of the revised QUICKI surrogate technique.