Predicting how many animals will be where: how to build, calibrate and evaluate individual-based models

Individual-based models (IBMs) can simulate the actions of individual animals as they interact with one another and the landscape in which they live. When used in spatially-explicit landscapes IBMs can show how populations change over time in response to management actions. For instance, IBMs are being used to design strategies of conservation and of the exploitation of fisheries, and for assessing the effects on populations of major construction projects and of novel agricultural chemicals. In such real world contexts, it becomes especially important to build IBMs in a principled fashion, and to approach calibration and evaluation systematically. We argue that insights from physiological and behavioural ecology offer a recipe for building realistic models, and that Approximate Bayesian Computation (ABC) is a promising technique for the calibration and evaluation of IBMs. IBMs are constructed primarily from knowledge about individuals. In ecological applications the relevant knowledge is found in physiological and behavioural ecology, and we approach these from an evolutionary perspective by taking into account how physiological and behavioural processes contribute to life histories, and how those life histories evolve. Evolutionary life history theory shows that, other things being equal, organisms should grow to sexual maturity as fast as possible, and then reproduce as fast as possible, while minimising per capita death rate. Physiological and behavioural ecology are largely built on these principles together with the laws of conservation of matter and energy. To complete construction of an IBM information is also needed on the effects of competitors, conspecifics and food scarcity; the maximum rates of ingestion, growth and reproduction, and life-history parameters. Using this knowledge about physiological and behavioural processes provides a principled way to build IBMs, but model parameters vary between species and are often difficult to measure. A common solution is to manually compare model outputs with observations from real landscapes and so to obtain parameters which produce acceptable fits of model to data. However, this procedure can be convoluted and lead to over-calibrated and thus inflexible models. Many formal statistical techniques are unsuitable for use with IBMs, but we argue that ABC offers a potential way forward. It can be used to calibrate and compare complex stochastic models and to assess the uncertainty in their predictions. We describe methods used to implement ABC in an accessible way and illustrate them with examples and discussion of recent studies. Although much progress has been made, theoretical issues remain, and some of these are outlined and discussed.

Stepwise construction of [Os6(CO)18(μ6-P)]-, a hexanuclear osmium cluster monoanion with trigonal-prismatic co-ordination for phosphorus: X-ray crystal structures of [Os6(μ-H)2(CO)20(MeCN)(μ3-PH)] and [PPh3Me][Os6(CO)18(μ6-P)]

Treatment of the labile cluster [Os3(CO)11(MeCN)] with PH3 affords the substituted product [Os3(CO)11(PH3)](1) in high yield. Subsequent reaction of (1) with Na2CO3 in MeOH, followed by acidification, gives the hydrido phosphido cluster [Os3(µ-H)(CO)10(µ-PH2)](2). When (2) is heated to 45–60 °C in the presence of [Os3(CO)11(MeCN)] a hexanuclear complex with the formulation [Os6(µ-H)2(CO)21(µ3-PH)](3) is obtained. If this reaction is repeated using [Os3(CO)10(MeCN)2] instead of [Os3(CO)11(MeCN)], an acetonitrile-containing product, [Os6(µ-H)2(CO)20(MeCN)(µ3-PH)](4), is obtained. An X-ray analysis of (4) shows that two Os3 triangular units are linked by a µ3-phosphinidene ligand, which symmetrically bridges an Os–Os edge of one triangle and is terminally co-ordinated to one Os atom of the second triangle. When (3) is treated with a weak base, such as [N(PPh3)2]Cl or [PPh3Me] Br, deprotonation to the corresponding cluster monoanion [Os6(µ-H)(CO)21(µ3-PH)]–(5) occurs. Treatment of (5) with a weak acid regenerates (3) in quantitative yield. Thermolysis of (3) leads to a closing up of the metal framework, affording the cluster [Os6(µ-H)(CO)18(µ6-P)], which readily deprotonates to give the anion [Os6(CO)18(µ6-P)]–(7) in the presence of [N(PPh3)2] Cl or [PPh3Me]Br. The same anion (7) may also be obtained by direct thermolysis of (5). An X-ray analysis of the [PPh3Me]+ salt of (7) confirms that the phosphorus occupies an interstitial site in a trigonal-prismatic hexaosmium framework, and co-ordinates to all six metal atoms with an average Os–P distance of 2.31 (1)Å. Proton and 31P n.m.r. data on all the new clusters are presented, and the position of the phosphorus resonance in the 31P n.m.r. spectrum is related to the changes in the environment of the phosphorus atom.

Social bees and food plant associations in the Nilgiri Biosphere Reserve, India

The diversity of social bees was assessed at 15 sites across five locations of the Nilgiri Biosphere Reserve, Western Ghats, India, from January to December 2007. We also conducted floristic analyses of local vegetation in each site using one-hectare sample plots. All woody species with a dbh (diameter at breast height) : 30 cm were recorded within the plots. A total area of 9.72 ha was assessed for floristic composition. Similarity of floristic composition between sites was determined using the Jaccard's distance measure and a dendrogram constructed based on the hierarchical clustering of floristic dissimilarities between sites. A Bee Importance Index (BII) was developed to give a measure of the bee diversity at each site. This index was a sum of the species richness of bee species in a site and their visitation frequencies to flowers, calculated as mean flower visits hour 1 within 2 focal patches within one hectare plots. The visits of bee species to flowers were also recorded. The Jaccard distance measure indicated that the montane sites were quite dissimilar to the low elevation sites in floristic diversity. The BII was 7-9 for the wet forest sites and ranged from 4-6 for drier forest sites. Seventy three plant species were identified as social bee plants and of them 45% were visited by one species of bee, 37% by two bee species and 18% by more than two bee species, indicating a certain degree of floral specialization among bees.

The anti-proliferative Effects of enterolactone in prostate cancer cells: evidence for the role of DNA licencing genes, mi-R106b cluster expression, and PTEN dosage

The mammalian lignan, enterolactone, has been shown to reduce the proliferation of the earlier stages of prostate cancer at physiological concentrations in vitro. However, efficacy in the later stages of the disease occurs at concentrations difficult to achieve through dietary modification. We have therefore investigated what concentration(s) of enterolactone can restrict proliferation in multiple stages of prostate cancer using an in vitro model system of prostate disease. We determined that enterolactone at 20 μM significantly restricted the proliferation of mid and late stage models of prostate disease. These effects were strongly associated with changes in the expression of the DNA licencing genes (GMNN, CDT1, MCM2 and 7), in reduced expression of the miR-106b cluster (miR-106b, miR-93, and miR-25), and in increased expression of the PTEN tumour suppressor gene. We have shown anti-proliferative effects of enterolactone in earlier stages of prostate disease than previously reported and that these effects are mediated, in part, by microRNA-mediated regulation.

Comparative Cluster/Double Star observations of the high and low latitude dayside magnetopause

The launch of the Double Star mission has provided the opportunity to monitor events at distinct locations on the dayside magnetopause, in coordination with the quartet of Cluster spacecraft. We present results of two such coordinated studies. In the first, 6 April 2004, both Cluster and the Double Star TC-1 spacecraft were on outbound transits through the dawn-side magnetosphere. Cluster observed northward moving FTEs with +/- polarity, whereas TC-1 saw -/+ polarity FTEs. The strength, motion and occurrence of the FTE signatures changes somewhat according to changes in IMF clock angle. These observations are consistent with ongoing reconnection on the dayside magnetopause, resulting in a series of flux transfer events (FTEs) seen both at Cluster and TC-1. The observed polarity and motion of each FTE signature advocates the existence of an active reconnection region consistently located between the positions of Cluster and TC-1, lying north and south of the reconnection line, respectively. This scenario is supported by the application of a model, designed to track flux tube motion, to conditions appropriate for the prevailing interplanetary conditions. The results from the model confirm the observational evidence that the low-latitude FTE dynamics is sensitive to changes in convected upstream conditions. In particular, changing the interplanetary magnetic field (IMF) clock angle in the model predicts that TC-1 should miss the resulting FTEs more often than Cluster, as is observed. For the second conjunction, on the 4 Jan 2005, the Cluster and TC-1 spacecraft all exited the dusk-side magnetosphere almost simultaneously, with TC-1 lying almost equatorial and Cluster at northern latitudes at about 4 RE from TC-1. The spacecraft traverse the magnetopause during a strong reversal in the IMF from northward to southward and a number of magnetosheath FTE signatures are subsequently observed. One coordinated FTE, studied in detail by Pu et al, [this issue], carries and inflowing energetic electron population and shows a motion and orientation which is similar at all spacecraft and consistent with the predictions of the model for the flux tube dynamics, given a near sub-solar reconnection line. This event can be interpreted either as the passage of two parallel flux tubes arising from adjacent x-line positions, or as a crossing of a single flux tube at different positions.

An integrated model of soil-canopy spectral radiances, photosynthesis, fluorescence, temperature and energy balance

This paper presents the model SCOPE (Soil Canopy Observation, Photochemistry and Energy fluxes), which is a vertical (1-D) integrated radiative transfer and energy balance model. The model links visible to thermal infrared radiance spectra (0.4 to 50 μm) as observed above the canopy to the fluxes of water, heat and carbon dioxide, as a function of vegetation structure, and the vertical profiles of temperature. Output of the model is the spectrum of outgoing radiation in the viewing direction and the turbulent heat fluxes, photosynthesis and chlorophyll fluorescence. A special routine is dedicated to the calculation of photosynthesis rate and chlorophyll fluorescence at the leaf level as a function of net radiation and leaf temperature. The fluorescence contributions from individual leaves are integrated over the canopy layer to calculate top-of-canopy fluorescence. The calculation of radiative transfer and the energy balance is fully integrated, allowing for feedback between leaf temperatures, leaf chlorophyll fluorescence and radiative fluxes. Leaf temperatures are calculated on the basis of energy balance closure. Model simulations were evaluated against observations reported in the literature and against data collected during field campaigns. These evaluations showed that SCOPE is able to reproduce realistic radiance spectra, directional radiance and energy balance fluxes. The model may be applied for the design of algorithms for the retrieval of evapotranspiration from optical and thermal earth observation data, for validation of existing methods to monitor vegetation functioning, to help interpret canopy fluorescence measurements, and to study the relationships between synoptic observations with diurnally integrated quantities. The model has been implemented in Matlab and has a modular design, thus allowing for great flexibility and scalability.

Transpolar voltage and polar cap flux during the substorm cycle and steady convection events

Transpolar voltages observed during traversals of the polar cap by the Defense Meteorological Satellite Program (DMSP) F-13 spacecraft during 2001 are analyzed using the expanding-contracting polar cap model of ionospheric convection. Each of the 10,216 passes is classified by its substorm phase or as a steady convection event (SCE) by inspection of the AE indices. For all phases, we detect a contribution to the transpolar voltage by reconnection in both the dayside magnetopause and in the crosstail current sheet. Detection of the IMF influence is 97% certain during quiet intervals and >99% certain during substorm/SCE growth phases but falls to 75% in substorm expansion phases: It is only 27% during SCEs. Detection of the influence of the nightside voltage is only 19% certain during growth phases, rising during expansion phases to a peak of 96% in recovery phases: During SCEs, it is >99%. The voltage during SCEs is dominated by the nightside, not the dayside, reconnection. On average, substorm expansion phases halt the growth phase rise in polar cap flux rather than reversing it. The main destruction of the excess open flux takes place during the 6- to 10-hour interval after the recovery phase (as seen in AE) and at a rate which is relatively independent of polar cap flux because the NENL has by then retreated to the far tail. The best estimate of the voltage associated with viscous-like transfer of closed field lines into the tail is around 10 kV.

An exploration of the drivers to bio-security collective action among a sample of UK cattle and sheep farmers

At present, collective action regarding bio-security among UK cattle and sheep farmers is rare. Despite the occurrence of catastrophic livestock diseases such as bovine spongiform encephalopathy (BSE) and foot and mouth disease (FMD), within recent decades, there are few national or local farmer-led animal health schemes. To explore the reasons for this apparent lack of interest, we utilised a socio-psychological approach to disaggregate the cognitive, emotive and contextual factors driving bio-security behaviour among cattle and sheep farmers in the United Kingdom (UK). In total, we interviewed 121 farmers in South-West England and Wales. The main analytical tools included a content, cluster and logistic regression analysis. The results of the content analysis illustrated apparent 'dissonance' between bio-security attitudes and behaviour.(1) Despite the heavy toll animal disease has taken on the agricultural economy, most study participants were dismissive of the many measures associated with bio-security. Justification for this lack of interest was largely framed in relation to the collective attribution or blame for the disease threats themselves. Indeed, epidemic diseases were largely related to external actors and agents. Reasons for outbreaks included inadequate border control, in tandem with ineffective policies and regulations. Conversely, endemic livestock disease was viewed as a problem for 'bad' farmers and not an issue for those individuals who managed their stock well. As such, there was little utility in forming groups to address what was largely perceived as an individual problem. Further, we found that attitudes toward bio-security did not appear to be influenced by any particular source of information per se. While strong negative attitudes were found toward specific sources of bio-security information, e.g. government leaflets, these appear to simply reflect widely held beliefs. In relation to actual bio-security behaviours, the logistic regression analysis revealed no significant difference between in-scheme and out of scheme farmers. We concluded that in order to support collective action with regard to bio-security, messages need to be reframed and delivered from a neutral source. Efforts to support group formation must also recognise and address the issues relating to perceptions of social connectedness among the communities involved. (c) 2008 Elsevier B.V. All rights reserved.

Managing threatened species – the ecological toolbox, evolutionary theory and the declining-population paradigm

1. The management of threatened species is an important practical way in which conservationists can intervene in the extinction process and reduce the loss of biodiversity. Understanding the causes of population declines (past, present and future) is pivotal to designing effective practical management. This is the declining-population paradigm identified by Caughley. 2. There are three broad classes of ecological tool used by conservationists to guide management decisions for threatened species: statistical models of habitat use, demographic models and behaviour-based models. Each of these is described here, illustrated with a case study and evaluated critically in terms of its practical application. 3. These tools are fundamentally different. Statistical models of habitat use and demographic models both use descriptions of patterns in abundance and demography, in relation to a range of factors, to inform management decisions. In contrast, behaviourbased models describe the evolutionary processes underlying these patterns, and derive such patterns from the strategies employed by individuals when competing for resources under a specific set of environmental conditions. 4. Statistical models of habitat use and demographic models have been used successfully to make management recommendations for declining populations. To do this, assumptions are made about population growth or vital rates that will apply when environmental conditions are restored, based on either past data collected under favourable environmental conditions or estimates of these parameters when the agent of decline is removed. As a result, they can only be used to make reliable quantitative predictions about future environments when a comparable environment has been experienced by the population of interest in the past. 5. Many future changes in the environment driven by management will not have been experienced by a population in the past. Under these circumstances, vital rates and their relationship with population density will change in the future in a way that is not predictable from past patterns. Reliable quantitative predictions about population-level responses then need to be based on an explicit consideration of the evolutionary processes operating at the individual level. 6. Synthesis and applications. It is argued that evolutionary theory underpins Caughley’s declining-population paradigm, and that it needs to become much more widely used within mainstream conservation biology. This will help conservationists examine critically the reliability of the tools they have traditionally used to aid management decision-making. It will also give them access to alternative tools, particularly when predictions are required for changes in the environment that have not been experienced by a population in the past.

Donepezil for the symptomatic treatment of patients with mild to moderate Alzheimer's disease: a meta-analysis of individual patient data from randomised controlled trials

Background: The objective was to evaluate the efficacy and tolerability of donepezil (5 and 10 mg/day) compared with placebo in alleviating manifestations of mild to moderate Alzheimer's disease (AD). Method: A systematic review of individual patient data from Phase II and III double-blind, randomised, placebo-controlled studies of up to 24 weeks and completed by 20 December 1999. The main outcome measures were the ADAS-cog, the CIBIC-plus, and reports of adverse events. Results: A total of 2376 patients from ten trials were randomised to either donepezil 5 mg/day (n = 821), 10 mg/day (n = 662) or placebo (n = 893). Cognitive performance was better in patients receiving donepezil than in patients receiving placebo. At 12 weeks the differences in ADAS-cog scores were 5 mg/day-placebo: - 2.1 [95% confidence interval (CI), - 2.6 to - 1.6; p < 0.001], 10 mg/day-placebo: - 2.5 ( - 3.1 to - 2.0; p < 0.001). The corresponding results at 24 weeks were - 2.0 ( - 2.7 to - 1.3; p < 0.001) and - 3.1 ( - 3.9 to - 2.4; p < 0.001). The difference between the 5 and 10 mg/day doses was significant at 24 weeks (p = 0.005). The odds ratios (OR) of improvement on the CIBIC-plus at 12 weeks were: 5 mg/day-placebo 1.8 (1.5 to 2.1; p < 0.001), 10 mg/day-placebo 1.9 (1.5 to 2.4; p < 0.001). The corresponding values at 24 weeks were 1.9 (1.5 to 2.4; p = 0.001) and 2.1 (1.6 to 2.8; p < 0.001). Donepezil was well tolerated; adverse events were cholinergic in nature and generally of mild severity and brief in duration. Conclusion: Donepezil (5 and 10 mg/day) provides meaningful benefits in alleviating deficits in cognitive and clinician-rated global function in AD patients relative to placebo. Increased improvements in cognition were indicated for the higher dose. Copyright © 2004 John Wiley & Sons, Ltd.

Bone morphogenetic proteins (BMP)-4,-6, and-7 potently suppress basal and luteinizing hormone-induced androgen production by bovine theca interna cells in primary culture: Could ovarian hyperandrogenic dysfunction be caused by a defect in thecal BMP signaling?

We reported recently that bovine theca interna cells in primary culture express several type-I and type-II receptors for bone morphogenetic proteins (BMPs). The same cells express at least two potential ligands for these receptors (BMP-4 and - 7), whereas bovine granulosa cells and oocytes express BMP-6. Therefore, BMPs of intrafollicular origin may exert autocrine/paracrine actions to modulate theca cell function. Here we report that BMP-4, - 6, and - 7 potently suppress both basal ( P < 0.0001; respective IC50 values, 0.78, 0.30, and 1.50 ng/ml) and LH-induced ( P < 0.0001; respective IC50 values, 5.00, 0.55, and 4.55 ng/ml) androgen production by bovine theca cells while having only a moderate effect on progesterone production and cell number. Semiquantitative RT-PCR showed that all three BMPs markedly reduced steady-state levels of mRNA for P450c17. Levels of mRNA encoding steroidogenic acute regulatory protein, P450scc, and 3 beta-hydroxysteroid dehydrogenase were also reduced but to a much lesser extent. Immunocytochemistry confirmed a marked reduction in cellular content of P450c17 protein after BMP treatment ( P < 0.001). Exposure to BMPs led to cellular accumulation of phosphorylated Smad1, but not Smad2, confirming that the receptors signal via a Smad1 pathway. The specificity of the BMP response was further explored by coincubating cells with BMPs and several potential BMP antagonists, chordin, gremlin, and follistatin. Gremlin and chordin were found to be effective antagonists of BMP-4 and - 7, respectively, and the observation that both antagonists enhanced ( P < 0.01) androgen production in the absence of exogenous BMP suggests an autocrine/paracrine role for theca-derived BMP- 4 and - 7 in modulating androgen production. Collectively, these data indicate that an intrafollicular BMP signaling pathway contributes to the negative regulation of thecal androgen production and that ovarian hyperandrogenic dysfunction could be a result of a defective autoregulatory pathway involving thecal BMP signaling.

Managing threatened species: the ecological toolbox, evolutionary theory and declining-population paradigm

1. The management of threatened species is an important practical way in which conservationists can intervene in the extinction process and reduce the loss of biodiversity. Understanding the causes of population declines (past, present and future) is pivotal to designing effective practical management. This is the declining-population paradigm identified by Caughley. 2. There are three broad classes of ecological tool used by conservationists to guide management decisions for threatened species: statistical models of habitat use, demographic models and behaviour-based models. Each of these is described here, illustrated with a case study and evaluated critically in terms of its practical application. 3. These tools are fundamentally different. Statistical models of habitat use and demographic models both use descriptions of patterns in abundance and demography, in relation to a range of factors, to inform management decisions. In contrast, behaviour-based models describe the evolutionary processes underlying these patterns, and derive such patterns from the strategies employed by individuals when competing for resources under a specific set of environmental conditions. 4. Statistical models of habitat use and demographic models have been used successfully to make management recommendations for declining populations. To do this, assumptions are made about population growth or vital rates that will apply when environmental conditions are restored, based on either past data collected under favourable environmental conditions or estimates of these parameters when the agent of decline is removed. As a result, they can only be used to make reliable quantitative predictions about future environments when a comparable environment has been experienced by the population of interest in the past. 5. Many future changes in the environment driven by management will not have been experienced by a population in the past. Under these circumstances, vital rates and their relationship with population density will change in the future in a way that is not predictable from past patterns. Reliable quantitative predictions about population-level responses then need to be based on an explicit consideration of the evolutionary processes operating at the individual level. 6. Synthesis and applications. It is argued that evolutionary theory underpins Caughley's declining-population paradigm, and that it needs to become much more widely used within mainstream conservation biology. This will help conservationists examine critically the reliability of the tools they have traditionally used to aid management decision-making. It will also give them access to alternative tools, particularly when predictions are required for changes in the environment that have not been experienced by a population in the past.

Analysis of the sulfate-reducing bacterial and methanogenic archaeal populations in contrasting Antarctic sediments

The distribution and activity of communities of sulfate-reducing bacteria (SRB) and methanogenic archaea in two contrasting Antarctic sediments were investigated. Methanogenesis dominated in freshwater Lake Heywood, while sulfate reduction dominated in marine Shallow Bay. Slurry experiments indicated that 90% of the methanogenesis in Lake Heywood was acetoclastic. This finding was supported by the limited diversity of clones detected in a Lake Heywood archaeal clone library, in which most clones were closely related to the obligate acetate-utilizing Methanosaeta concilii. The Shallow Bay archaeal clone library contained clones related to the C-1-utilizing Methanolobus and Methanococcoides and the H-2-utilizing Methanogenium. Oligonucleotide probing of RNA extracted directly from sediment indicated that archaea represented 34% of the total prokaryotic signal in Lake Heywood and that Methanosaeta was a major component (13.2%) of this signal. Archaea represented only 0.2% of the total prokaryotic signal in RNA extracted from Shallow Bay sediments. In the Shallow Bay bacterial clone library, 10.3% of the clones were SRB-like, related to Desulfotalea/Desulforhopalus, Desulfofaba, Desulfosarcina, and Desulfobacter as well as to the sulfur and metal oxidizers comprising the Desulfuromonas cluster. Oligonucleotide probes for specific SRB clusters indicated that SRB represented 14.7% of the total prokaryotic signal, with Desulfotalea/Desulforhopalus being the dominant SRB group (10.7% of the total prokaryotic signal) in the Shallow Bay sediments; these results support previous results obtained for Arctic sediments. Methanosaeta and Desulfotalea/Desulforhopalus appear to be important in Lake Heywood and Shallow Bay, respectively, and may be globally important in permanently low-temperature sediments.