Population-based continuous optimization, probabilistic modelling and mean shift

Evolutionary algorithms perform optimization using a population of sample solution points. An interesting development has been to view population-based optimization as the process of evolving an explicit, probabilistic model of the search space. This paper investigates a formal basis for continuous, population-based optimization in terms of a stochastic gradient descent on the Kullback-Leibler divergence between the model probability density and the objective function, represented as an unknown density of assumed form. This leads to an update rule that is related and compared with previous theoretical work, a continuous version of the population-based incremental learning algorithm, and the generalized mean shift clustering framework. Experimental results are presented that demonstrate the dynamics of the new algorithm on a set of simple test problems.

A stochastic metapopulation model accounting for habitat dynamics

A stochastic metapopulation model accounting for habitat dynamics is presented. This is the stochastic SIS logistic model with the novel aspect that it incorporates varying carrying capacity. We present results of Kurtz and Barbour, that provide deterministic and diffusion approximations for a wide class of stochastic models, in a form that most easily allows their direct application to population models. These results are used to show that a suitably scaled version of the metapopulation model converges, uniformly in probability over finite time intervals, to a deterministic model previously studied in the ecological literature. Additionally, they allow us to establish a bivariate normal approximation to the quasi-stationary distribution of the process. This allows us to consider the effects of habitat dynamics on metapopulation modelling through a comparison with the stochastic SIS logistic model and provides an effective means for modelling metapopulations inhabiting dynamic landscapes.

A modelling approach to estimate the effect of exotic pollinators on exotic weed population dynamics: bumblebees and broom in Australia

The role of mutualisms in contributing to species invasions is rarely considered, inhibiting effective risk analysis and management options. Potential ecological consequences of invasion of non-native pollinators include increased pollination and seed set of invasive plants, with subsequent impacts on population growth rates and rates of spread. We outline a quantitative approach for evaluating the impact of a proposed introduction of an invasive pollinator on existing weed population dynamics and demonstrate the use of this approach on a relatively data-rich case study: the impacts on Cytisus scoparius (Scotch broom) from proposed introduction of Bombus terrestris. Three models have been used to assess population growth (matrix model), spread speed (integrodifference equation), and equilibrium occupancy (lattice model) for C. scoparius. We use available demographic data for an Australian population to parameterize two of these models. Increased seed set due to more efficient pollination resulted in a higher population growth rate in the density-independent matrix model, whereas simulations of enhanced pollination scenarios had a negligible effect on equilibrium weed occupancy in the lattice model. This is attributed to strong microsite limitation of recruitment in invasive C. scoparius populations observed in Australia and incorporated in the lattice model. A lack of information regarding secondary ant dispersal of C. scoparius prevents us from parameterizing the integrodifference equation model for Australia, but studies of invasive populations in California suggest that spread speed will also increase with higher seed set. For microsite-limited C. scoparius populations, increased seed set has minimal effects on equilibrium site occupancy. However, for density-independent rapidly invading populations, increased seed set is likely to lead to higher growth rates and spread speeds. The impacts of introduced pollinators on native flora and fauna and the potential for promoting range expansion in pollinator-limited 'sleeper weeds' also remain substantial risks.

Collective coherent population trapping in a thermal field

We analyze the efficiency of coherent population trapping (CPT) in a superposition of the ground states of three-level atoms under the influence of the decoherence process induced by a broadband thermal field. We show that in a single atom there is no perfect CPT when the atomic transitions are affected by the thermal field. The perfect CPT may occur when only one of the two atomic transitions is affected by the thermal field. In the case when both atomic transitions are affected by the thermal field, we demonstrate that regardless of the intensity of the thermal field the destructive effect on the CPT can be circumvented by the collective behavior of the atoms. An analytic expression was obtained for the populations of the upper atomic levels which can be considered as a measure of the level of thermal decoherence. The results show that the collective interaction between the atoms can significantly enhance the population trapping in that the population of the upper state decreases with an increased number of atoms. The physical origin of this feature is explained by the semiclassical dressed-atom model of the system. We introduce the concept of multiatom collective coherent population trapping by demonstrating the existence of collective (entangled) states whose storage capacity is larger than that of the equivalent states of independent atoms.

Genetic population structure and call variation in a passerine bird, the satin bowerbird, Ptilonorhynchus violaceus

Geographic variation in vocalizations is widespread in passerine birds, but its origins and maintenance remain unclear. One hypothesis to explain this variation is that it is associated with geographic isolation among populations and therefore should follow a vicariant pattern similar to that typically found in neutral genetic markers. Alternatively, if environmental selection strongly influences vocalizations, then genetic divergence and vocal divergence may be disassociated. This study compared genetic divergence derived from 11 microsatellite markers with a metric of phenotypic divergence derived from male bower advertisement calls. Data were obtained from 16 populations throughout the entire distribution of the satin bowerbird, an Australian wet-forest-restricted passerine. There was no relationship between call divergence and genetic divergence, similar to most other studies on birds with learned vocalizations. Genetic divergence followed a vicariant model of evolution, with the differentiation of isolated populations and isolation-by-distance among continuous populations. Previous work on Ptilonorhynchus violaceus has shown that advertisement call structure is strongly influenced by the acoustic environment of different habitats. Divergence in vocalizations among genetically related populations in different habitats indicates that satin bowerbirds match their vocalizations to the environment in which they live, despite the homogenizing influence of gene flow. In combination with convergence of vocalizations among genetically divergent populations occurring in the same habitat, this shows the overriding importance that habitat-related selection can have on the establishment and maintenance of variation in vocalizations.

Kidney side population reveals multilineage potential and renal functional capacity but also cellular heterogeneity

Side population (SP) cells in the adult kidney are proposed to represent a progenitor population. However, the size, origin, phenotype, and potential of the kidney SP has been controversial. In this study, the SP fraction of embryonic and adult kidneys represented 0.1 to 0.2% of the total viable cell population. The immunophenotype and the expression profile of kidney SP cells was distinct from that of bone marrow SP cells, suggesting that they are a resident nonhematopoietic cell population. Affymetrix expression profiling implicated a role for Notch signaling in kidney SP cells and was used to identify markers of kidney SP. Localization by in situ hybridization confirmed a primarily proximal tubule location, supporting the existence of a tubular niche, but also revealed considerable heterogeneity, including the presence of renal macrophages. Adult kidney SP cells demonstrated multilineage differentiation in vitro, whereas microinjection into mouse metanephroi showed that SP cells had a 3.5- to 13-fold greater potential to contribute to developing kidney than non-SP main population cells. However, although reintroduction of SP cells into an Adriamycin-nephropathy model reduced albuminuria:creatinine ratios, this was without significant tubular integration, suggesting a humoral role for SP cells in renal repair. The heterogeneity of the renal SP highlights the need for further fractionation to distinguish the cellular subpopulations that are responsible for the observed multilineage capacity and transdifferentiative and humoral activities.

Predicting population dynamics of weed biological control agents: science or gazing into crystal balls?

Various factors can influence the population dynamics of phytophages post introduction, of which climate is fundamental. Here we present an approach, using a mechanistic modelling package (CLIMEX), that at least enables one to make predictions of likely dynamics based on climate alone. As biological control programs will have minimal funding for basic work (particularly on population dynamics), we show how predictions can be made using a species geographical distribution, relative abundance across its range, seasonal phenology and laboratory rearing data. Many of these data sets are more likely to be available than long-term population data, and some can be incorporated into the exploratory phase of a biocontrol program. Although models are likely to be more robust the more information is available, useful models can be developed using information on species distribution alone. The fitted model estimates a species average response to climate, and can be used to predict likely geographical distribution if introduced, where the agent is likely to be more abundant (i.e. good locations) and more importantly for interpretation of release success, the likely variation in abundance over time due to intra- and inter-year climate variability. The latter will be useful in predicting both the seasonal and long-term impacts of the potential biocontrol agent on the target weed. We believe this tool may not only aid in the agent selection process, but also in the design of release strategies, and for interpretation of post-introduction dynamics and impacts. More importantly we are making testable predictions. If biological control is to become more of a science making and testing such hypothesis will be a key component.

Imitating the neighbours: vocal dialect matching in a mimic-model system

Vocal mimicry provides a unique system for investigating song learning and cultural evolution in birds. Male lyrebirds produce complex vocal displays that include extensive and accurate mimicry of many other bird species. We recorded and analysed the songs of the Albert's lyrebird (Menura alberti) and its most commonly imitated model species, the satin bowerbird (Ptilonorhynchus violaceus), at six sites in southeast Queensland, Australia. We show that each population of lyrebirds faithfully reproduces the song of the local population of bowerbirds. Within a population, lyrebirds show less variation in song structure than the available variation in the songs of the models. These results provide the first quantitative evidence for dialect matching in the songs of two species that have no direct ecological relationship.

Population pharmacokinetics of itraconazole and its active metabolite hydroxy-itraconazole in paediatric cystic fibrosis and bone marrow transplant patients

Objective: The objective of the study was to characterise the population pharmacokinetic properties of itraconazole and its active metabolite hydroxyitraconazole in a representative paediatric population of cystic fibrosis and bone marrow transplant (BMT) patients and to identify patient characteristics influencing the pharmacokinetics of itraconazole. The ultimate goals were to determine the relative bioavailability between the two oral formulations (capsules vs oral solution) and to optimise dosing regimens in these patients. Methods: All paediatric patients with cystic fibrosis or patients undergoing BMT at The Royal Children's Hospital, Brisbane, QLD, Australia, who were prescribed oral itraconazole for the treatment of allergic bronchopulmonary aspergillosis (cystic fibrosis patients) or for prophylaxis of any fungal infection (BMT patients) were eligible for the study. Blood samples were taken from the recruited patients as per an empirical sampling design either during hospitalisation or during outpatient clinic visits. ltraconazole and hydroxy-itraconazole plasma concentrations were determined by a validated high-performance liquid chromatography assay with fluorometric detection. A nonlinear mixed-effect modelling approach using the NONMEM software to simultaneously describe the pharmacokinetics of itraconazole and its metabolite. Results: A one-compartment model with first-order absorption described the itraconazole data, and the metabolism of the parent drug to hydroxy-itraconazole was described by a first-order rate constant. The metabolite data also showed one-compartment characteristics with linear elimination. For itraconazole the apparent clearance (CLitraconazole) was 35.5 L/hour, the apparent volume of distribution (V-d(itraconazole)) was 672L, the absorption rate constant for the capsule formulation was 0.0901 h(-1) and for the oral solution formulation was 0.96 h-1. The lag time was estimated to be 19.1 minutes and the relative bioavailability between capsules and oral solution (F-rel) was 0.55. For the metabolite, volume of distribution, V-m/(F (.) f(m)), and clearance, CL/(F (.) fm), were 10.6L and 5.28 L/h, respectively. The influence of total bodyweight was significant, added as a covariate on CLitraconazoie/F and V-d(itraconazole)/F (standardised to a 70kg person) using allometric three-quarter power scaling on CLitraconazole/F, which therefore reflected adult values. The unexplained between-subject variability (coefficient of variation %) was 68.7%, 75.8%, 73.4% and 61.1% for CLitraconazoie/F, Vd(itraconazole)/F, CLm/(F (.) fm) and F-rel, respectively. The correlation between random effects of CLitraconazole and Vd((itraconazole)) was 0.69. Conclusion: The developed population pharmacokinetic model adequately described the pharmacokinetics of itraconazole and its active metabolite, hydroxy-itraconazole, in paediatric patients with either cystic fibrosis or undergoing BMT. More appropriate dosing schedules have been developed for the oral solution and the capsules to secure a minimum therapeutic trough plasma concentration of 0.5 mg/L for these patients.

Migraine With Aura and Migraine Without Aura Are Not Distinct Entities: Further Evidence From a Large Dutch Population Study

It is often debated whether migraine with aura (MA) and migraine without aura (MO) are etiologically distinct disorders. A previous study using latent class analysis (LCA) in Australian twins showed no evidence for separate subtypes of MO and MA. The aim of the present study was to replicate these results in a population of Dutch twins and their parents, siblings and partners (N = 10,144). Latent class analysis of International Headache Society (IHS)-based migraine symptoms resulted in the identification of 4 classes: a class of unaffected subjects (class 0), a mild form of nonmigrainous headache (class 1), a moderately severe type of migraine (class 2), typically without neurological symptoms or aura (8% reporting aura symptoms), and a severe type of migraine (class 3), typically with neurological symptoms, and aura symptoms in approximately half of the cases. Given the overlap of neurological symptoms and nonmutual exclusivity of aura symptoms, these results do not support the MO and MA subtypes as being etiologically distinct. The heritability in female twins of migraine based on LCA classification was estimated at .50 (95% confidence intervals [0CI} .27 -.59), similar to IHS-based migraine diagnosis (h(2) = .49, 95% Cl .19-.57). However, using a dichotomous classification (affected-unaffected) decreased heritability for the IHS-based classification (h(2) = .33, 95% Cl .00-.60), but not the LCA-based classification (h(2) = .51, 95% Cl. 23-.61). Importantly, use of the LCA-based classification increased the number of subjects classified as affected. The heritability of the screening question was similar to more detailed LCA and IHS classifications, suggesting that the screening procedure is an important determining factor in genetic studies of migraine.

On parameter estimation in population models

We describe methods for estimating the parameters of Markovian population processes in continuous time, thus increasing their utility in modelling real biological systems. A general approach, applicable to any finite-state continuous-time Markovian model, is presented, and this is specialised to a computationally more efficient method applicable to a class of models called density-dependent Markov population processes. We illustrate the versatility of both approaches by estimating the parameters of the stochastic SIS logistic model from simulated data. This model is also fitted to data from a population of Bay checkerspot butterfly (Euphydryas editha bayensis), allowing us to assess the viability of this population. (c) 2006 Elsevier Inc. All rights reserved.

Accounting for management costs in sensitivity analyses of matrix population models

Traditional sensitivity and elasticity analyses of matrix population models have been used to p inform management decisions, but they ignore the economic costs of manipulating vital rates. For exam le, the growth rate of a population is often most sensitive to changes in adult survival rate, but this does not mean that increasing that rate is the best option for managing the population because it may be much more expensive than other options. To explore how managers should optimize their manipulation of vital rates, we incorporated the cost of changing those rates into matrix population models. We derived analytic expressions for locations in parameter space where managers should shift between management of fecundity and survival, for the balance between fecundity and survival management at those boundaries, and for the allocation of management resources to sustain that optimal balance. For simple matrices, the optimal budget allocation can often be expressed as simple functions of vital rates and the relative costs of changing them. We applied our method to management of the Helmeted Honeyeater (Lichenostomus melanops cassidix; an endangered Australian bird) and the koala (Phascolarctos cinereus) as examples. Our method showed that cost-efficient management of the Helmeted Honeyeater should focus on increasing fecundity via nest protection, whereas optimal koala management should focus on manipulating both fecundity and survival simultaneously, These findings are contrary to the cost-negligent recommendations of elasticity analysis, which would suggest focusing on managing survival in both cases. A further investigation of Helmeted Honeyeater management options, based on an individual-based model incorporating density dependence, spatial structure, and environmental stochasticity, confirmed that fecundity management was the most cost-effective strategy. Our results demonstrate that decisions that ignore economic factors will reduce management efficiency.

Optimizing Presence–Absence Surveys For Detecting Population Trends

Presence-absence surveys are a commonly used method for monitoring broad-scale changes in wildlife distributions. However, the lack of power of these surveys for detecting population trends is problematic for their application in wildlife management. Options for improving power include increasing the sampling effort or arbitrarily relaxing the type I error rate. We present an alternative, whereby targeted sampling of particular habitats in the landscape using information from a habitat model increases power. The advantage of this approach is that it does not require a trade-off with either cost or the Pr(type I error) to achieve greater power. We use a demographic model of koala (Phascolarctos cinereus) population dynamics and simulations of the monitoring process to estimate the power to detect a trend in occupancy for a range of strategies, thereby demonstrating that targeting particular habitat qualities can improve power substantially. If the objective is to detect a decline in occupancy, the optimal strategy is to sample high-quality habitats. Alternatively, if the objective is to detect an increase in occupancy, the optimal strategy is to sample intermediate-quality habitats. The strategies with the highest power remained the same under a range of parameter assumptions, although observation error had a strong influence on the optimal strategy. Our approach specifically applies to monitoring for detecting long-term trends in occupancy or abundance. This is a common and important monitoring objective for wildlife managers, and we provide guidelines for more effectively achieving it.

Investigating the use of sweet sorghum as a model for sugar accumulation in sugarcane

This paper outlines a current investigation of sugar accumulation in sweet sorghum to assist in understanding and simplifying this complex trait in sugarcane. A recombinant inbred line (RIL) sorghum population, between a sweet and a grain sorghum, has been developed and phenotyped for various morphological and agronomic traits related to grain yield, biomass and stem sugar content. A genetic linkage map will be constructed for the sweet sorghum population with the objective of identifying genomic regions associated with sucrose accumulation in sweet sorghum. This will lead to further work, including comparative mapping in sugarcane, to identify the extent to which sweet sorghum can be used as a model for investigating sugar accumulation in sugarcane.

Bayesian population pharmacokinetic analysis of sirolimus

Objective: It is usual that data collected from routine clinical care is sparse and unable to support the more complex pharmacokinetic (PK) models that may have been reported in previous rich data studies. Informative priors may be a pre-requisite for model development. The aim of this study was to estimate the population PK parameters of sirolimus using a fully Bayesian approach with informative priors. Methods: Informative priors including prior mean and precision of the prior mean were elicited from previous published studies using a meta-analytic technique. Precision of between-subject variability was determined by simulations from a Wishart distribution using MATLAB (version 6.5). Concentration-time data of sirolimus retrospectively collected from kidney transplant patients were analysed using WinBUGS (version 1.3). The candidate models were either one- or two-compartment with first order absorption and first order elimination. Model discrimination was based on computation of the posterior odds supporting the model. Results: A total of 315 concentration-time points were obtained from 25 patients. Most data were clustered at trough concentrations with range of 1.6 to 77 hours post-dose. Using informative priors, either a one- or two-compartment model could be used to describe the data. When a one-compartment model was applied, information was gained from the data for the value of apparent clearance (CL/F = 18.5 L/h), and apparent volume of distribution (V/F = 1406 L) but no information was gained about the absorption rate constant (ka). When a two-compartment model was fitted to the data, the data were informative about CL/F, apparent inter-compartmental clearance, and apparent volume of distribution of the peripheral compartment (13.2 L/h, 20.8 L/h, and 579 L, respectively). The posterior distribution of the volume distribution of central compartment and ka were the same as priors. The posterior odds for the two-compartment model was 8.1, indicating the data supported the two-compartment model. Conclusion: The use of informative priors supported the choice of a more complex and informative model that would otherwise have not been supported by the sparse data.