Increase of a Caribbean leatherback turtle Dermochelys coriacea nesting population linked to long-term nest protection

The leatherback turtle Dermochelys coriacea is considered to be at serious risk of global extinction, despite ongoing conservation efforts. Intensive long-term monitoring of a leatherback nesting population on Sandy Point (St. Croix, US Virgin Islands) offers a unique opportunity to quantify basic population parameters and evaluate effectiveness of nesting beach conservation practices. We report a significant increase in the number of females nesting annually from ca. 18-30 in the 1980s to 186 in 2001, with a corresponding increase in annual hatchling production from ca. 2000 to over 49,000. We then analyzed resighting data from 1991 to 2001 with an open robust-design capture-mark-recapture model to estimate annual nester survival and adult abundance for this population. The expected annual survival probability was estimated at ca. 0.893 (95% CL 0.87-0.92) and the population was estimated to be increasing ca. 13% pa since the early 1990s. Taken together with DNA fingerprinting that identify mother-daughter relations, our findings suggest that the increase in the size of the nesting population since 1991 was probably due to an aggressive program of beach protection and egg relocation initiated more than 20 years ago. Beach protection and egg relocation provide a simple and effective conservation strategy for this Northern Caribbean nesting population as long as adult survival at sea remains relatively high. (c) 2005 Elsevier Ltd. All rights reserved.

Optimal design for model discrimination and parameter estimation for itraconazole population pharmacokinetics in cystic fibrosis patients

Optimal sampling times are found for a study in which one of the primary purposes is to develop a model of the pharmacokinetics of itraconazole in patients with cystic fibrosis for both capsule and solution doses. The optimal design is expected to produce reliable estimates of population parameters for two different structural PK models. Data collected at these sampling times are also expected to provide the researchers with sufficient information to reasonably discriminate between the two competing structural models.

Analysis of population pharmacokinetic data using NONMEM and WinBUGS

The aim of this report is to describe the use of WinBUGS for two datasets that arise from typical population pharmacokinetic studies. The first dataset relates to gentamicin concentration-time data that arose as part of routine clinical care of 55 neonates. The second dataset incorporated data from 96 patients receiving enoxaparin. Both datasets were originally analyzed by using NONMEM. In the first instance, although NONMEM provided reasonable estimates of the fixed effects parameters it was unable to provide satisfactory estimates of the between-subject variance. In the second instance, the use of NONMEM resulted in the development of a successful model, albeit with limited available information on the between-subject variability of the pharmacokinetic parameters. WinBUGS was used to develop a model for both of these datasets. Model comparison for the enoxaparin dataset was performed by using the posterior distribution of the log-likelihood and a posterior predictive check. The use of WinBUGS supported the same structural models tried in NONMEM. For the gentamicin dataset a one-compartment model with intravenous infusion was developed, and the population parameters including the full between-subject variance-covariance matrix were available. Analysis of the enoxaparin dataset supported a two compartment model as superior to the one-compartment model, based on the posterior predictive check. Again, the full between-subject variance-covariance matrix parameters were available. Fully Bayesian approaches using MCMC methods, via WinBUGS, can offer added value for analysis of population pharmacokinetic data.

Protocols for listing threatened species can forecast extinction

Risk-ranking protocols are used widely to classify the conservation status of the world's species. Here we report on the first empirical assessment of their reliability by using a retrospective study of 18 pairs of bird and mammal species (one species extinct and the other extant) with eight different assessors. The performance of individual assessors varied substantially, but performance was improved by incorporating uncertainty in parameter estimates and consensus among the assessors. When this was done, the ranks from the protocols were consistent with the extinction outcome in 70-80% of pairs and there were mismatches in only 10-20% of cases. This performance was similar to the subjective judgements of the assessors after they had estimated the range and population parameters required by the protocols, and better than any single parameter. When used to inform subjective judgement, the protocols therefore offer a means of reducing unpredictable biases that may be associated with expert input and have the advantage of making the logic behind assessments explicit. We conclude that the protocols are useful for forecasting extinctions, although they are prone to some errors that have implications for conservation. Some level of error is to be expected, however, given the influence of chance on extinction. The performance of risk assessment protocols may be improved by providing training in the application of the protocols, incorporating uncertainty in parameter estimates and using consensus among multiple assessors, including some who are experts in the application of the protocols. Continued testing and refinement of the protocols may help to provide better absolute estimates of risk, particularly by re-evaluating how the protocols accommodate missing data.

Slowing down a pine invasion despite uncertainty in demography and dispersal

1. Management decisions regarding invasive plants often have to be made quickly and in the face of fragmentary knowledge of their population dynamics. However, recommendations are commonly made on the basis of only a restricted set of parameters. Without addressing uncertainty and variability in model parameters we risk ineffective management, resulting in wasted resources and an escalating problem if early chances to control spread are missed. 2. Using available data for Pinus nigra in ungrazed and grazed grassland and shrubland in New Zealand, we parameterized a stage-structured spread model to calculate invasion wave speed, population growth rate and their sensitivities and elasticities to population parameters. Uncertainty distributions of parameters were used with the model to generate confidence intervals (CI) about the model predictions. 3. Ungrazed grassland environments were most vulnerable to invasion and the highest elasticities and sensitivities of invasion speed were to long-distance dispersal parameters. However, there was overlap between the elasticity and sensitivity CI on juvenile survival, seedling establishment and long-distance dispersal parameters, indicating overlap in their effects on invasion speed. 4. While elasticity of invasion speed to long-distance dispersal was highest in shrubland environments, there was overlap with the CI of elasticity to juvenile survival. In shrubland invasion speed was most sensitive to the probability of establishment, especially when establishment was low. In the grazed environment elasticity and sensitivity of invasion speed to the severity of grazing were consistently highest. Management recommendations based on elasticities and sensitivities depend on the vulnerability of the habitat. 5. Synthesis and applications. Despite considerable uncertainty in demography and dispersal, robust management recommendations emerged from the model. Proportional or absolute reductions in long-distance dispersal, juvenile survival and seedling establishment parameters have the potential to reduce wave speed substantially. Plantations of wind-dispersed invasive conifers should not be sited on exposed sites vulnerable to long-distance dispersal events, and trees in these sites should be removed. Invasion speed can also be reduced by removing seedlings, establishing competitive shrubs and grazing. Incorporating uncertainty into the modelling process increases our confidence in the wide applicability of the management strategies recommended here.

Responses in carcass lean pH at 24 hours post-mortem in lines of Large White pigs selected for growth rate on a fixed ration over a set time

A restricted maximum likelihood analysis applied to an animal model showed no significant differences (P > 0.05) in pH value of the longissimus dorsi measured at 24 h post-mortem (pH24) between high and low lines of Large White pigs selected over 4 years for post-weaning growth rate on restricted feeding. Genetic and phenotypic correlations between pH24 and production and carcass traits were estimated using all performance testing records combined with the pH24 measurements (5.05–7.02) on slaughtered animals. The estimate of heritability for pH24 was moderate (0.29 ± 0.18). Genetic correlations between pH24 and production or carcass composition traits, except for ultrasonic backfat (UBF), were not significantly different from zero. UBF had a moderate, positive genetic correlation with pH24 (0.24 ± 0.33). These estimates of genetic correlations affirmed that selection for increased growth rate on restricted feeding is likely to result in limited changes in pH24 and pork quality since the selection does not put a high emphasis on reduced fatness.

Toward better outcomes with tacrolimus therapy: Population pharmacokinetics and individualized dosage prediction in adult liver transplantation

Patient outcomes in transplantation would improve if dosing of immunosuppressive agents was individualized. The aim of this study is to develop a population pharmacokinetic model of tacrolimus in adult liver transplant recipients and test this model in individualizing therapy. Population analysis was performed on data from 68 patients. Estimates were sought for apparent clearance (CL/F) and apparent volume of distribution (V/F) using the nonlinear mixed effects model program (NONMEM). Factors screened for influence on these parameters were weight, age, sex, transplant type, biliary reconstructive procedure, postoperative day, days of therapy, liver function test results, creatinine clearance, hematocrit, corticosteroid dose, and interacting drugs. The predictive performance of the developed model was evaluated through Bayesian forecasting in an independent cohort of 36 patients. No linear correlation existed between tacrolimus dosage and trough concentration (r(2) = 0.005). Mean individual Bayesian estimates for CL/F and V/F were 26.5 8.2 (SD) L/hr and 399 +/- 185 L, respectively. CL/F was greater in patients with normal liver function. V/F increased with patient weight. CL/F decreased with increasing hematocrit. Based on the derived model, a 70-kg patient with an aspartate aminotransferase (AST) level less than 70 U/L would require a tacrolimus dose of 4.7 mg twice daily to achieve a steady-state trough concentration of 10 ng/mL. A 50-kg patient with an AST level greater than 70 U/L would require a dose of 2.6 mg. Marked interindividual variability (43% to 93%) and residual random error (3.3 ng/mL) were observed. Predictions made using the final model were reasonably nonbiased (0.56 ng/mL), but imprecise (4.8 ng/mL). Pharmacokinetic information obtained will assist in tacrolimus dosing; however, further investigation into reasons for the pharmacokinetic variability of tacrolimus is required.

Of sheep and rain: large-scale population dynamics of the red kangaroo

1. We analysed time-series data from populations of red kangaroos (Macropus rufus, Desmarest) inhabiting four areas in the pastoral zone of South Australia. We formulated a set of a priori models to disentangle the relative effects of the covariates: rainfall, harvesting, intraspecific competition, and domestic herbivores, on kangaroo population-growth rate. 2. The statistical framework allowed for spatial variation in the growth-rate parameters, response to covariates, and environmental variability, as well as spatially correlated error terms due to shared environment. 3. The most parsimonious model included all covariates but no area-specific parameter values, suggesting that kangaroo densities respond in the same way to the covariates across the areas. 4. The temporal dynamics were spatially correlated, even after taking into account the potentially synchronizing effect of rainfall, harvesting and domestic herbivores. 5. Counter-intuitively, we found a positive rather than negative effect of domestic herbivore density on the population-growth rate of kangaroos. We hypothesize that this effect is caused by sheep and cattle acting as a surrogate for resource availability beyond rainfall. 6. Even though our system is well studied, we must conclude that approximating resources by surrogates such as rainfall is more difficult than previously thought. This is an important message for studies of consumer-resource systems and highlights the need to be explicit about population processes when analysing population patterns.

Some considerations on the design of population pharmacokinetic studies

The goal of this manuscript is to introduce a framework for consideration of designs for population pharmacokinetic orpharmacokinetic-pharmacodynamic studies. A standard one compartment pharmacokinetic model with first-order input and elimination is considered. A series of theoretical designs are considered that explore the influence of optimizing the allocation of sampling times, allocating patients to elementary designs, consideration of sparse sampling and unbalanced designs and also the influence of single vs. multiple dose designs. It was found that what appears to be relatively sparse sampling (less blood samples per patient than the number of fixed effects parameters to estimate) can also be highly informative. Overall, it is evident that exploring the population design space can yield many parsimonious designs that are efficient for parameter estimation and that may not otherwise have been considered without the aid of optimal design theory.

The population pharmacokinetics of citalopram after deliberate self-poisoning: A Bayesian approach

Defining the pharmacokinetics of drugs in overdose is complicated. Deliberate self-poisoning is generally impulsive and associated with poor accuracy in dose history. In addition, early blood samples are rarely collected to characterize the whole plasma-concentration time profile and the effect of decontamination on the pharmacokinetics is uncertain. The aim of this study was to explore a fully Bayesian methodology for population pharmacokinetic analysis of data that arose from deliberate self-poisoning with citalopram. Prior information on the pharmacokinetic parameters was elicited from 14 published studies on citalopram when taken in therapeutic doses. The data set included concentration-time data from 53 patients studied after 63 citalopram overdose events (dose range: 20-1700 mg). Activated charcoal was administered between 0.5 and 4 h after 17 overdose events. The clinical investigator graded the veracity of the patients' dosing history on a 5-point ordinal scale. Inclusion of informative priors stabilised the pharmacokinetic model and the population mean values could be estimated well. There were no indications of non-linear clearance after excessive doses. The final model included an estimated uncertainty of the dose amount which in a simulation study was shown to not affect the model's ability to characterise the effects of activated charcoal. The effect of activated charcoal on clearance and bioavailability was pronounced and resulted in a 72% increase and 22% decrease, respectively. These findings suggest charcoal administration is potentially beneficial after citalopram overdose. The methodology explored seems promising for exploring the dose-exposure relationship in the toxicological settings.

The effect of temperature on the size and population density of dinoflagellates in larvae of the reef coral Porites astreoides

Pre-settlement events play an important role in determining larval success in marine invertebrates with bentho-pelagic life histories, yet the consequences of these events typically are not well understood. The purpose of this study was to examine the pre-settlement impacts of different seawater temperatures on the size and population density of dinoflagellate symbionts in brooded larvae of the Caribbean coral Porites astreoides. Larvae were collected from P. astreoides at 14-20 m depth on Conch Reef (Florida) in June 2002, and incubated for 24 h at 15 temperatures spanning the range 25.1 degrees-30.0 degrees C in mean increments of 0.4 +/- 0.1 degrees C (+/- SD). The most striking feature of the larval responses was the magnitude of change in both parameters across this 5 degrees C temperature range within 24 h. In general, larvae were largest and had the highest population densities of Symbiodinium sp. between 26.4 degrees-27.7 degrees C, and were smallest and had the lowest population densities at 25.8 degrees C and 28.8 degrees C. Larval size and symbiont population density were elevated slightly (relative to the minimal values) at the temperature extremes of 25.1 degrees C and 30 degrees C. These data demonstrate that coral larvae are highly sensitive to seawater temperature during their pelagic phase, and respond through changes in size and the population densities of Symbiodinium sp. to ecologically relevant temperature signals within 24 h. The extent to which these changes are biologically meaningful will depend on the duration and frequency of exposure of coral larvae to spatio-temporal variability in seawater temperature, and whether the responses have cascading effects on larval success and their entry to the post-settlement and recruitment phase.

Population-based detection of systolic and diastolic dysfunction with amino-terminal pro-B-type natriuretic peptide

Background There is limited information regarding the clinical utility of amino-terminal pro-B-type natriuretic pepticle (NT-proBNP) for the detection of left ventricular (LV) dysfunction in the community. We evaluated predictors of circulating NT-proBNP levels and determined the utility of NT-proBNP to detect systolic and diastolic LV dysfunction in older adults. Methods. A population-based sample of 1229 older adults (mean age 69.4 years, 50.1% women) underwent echocardiographic assessment of cardiac structure and function and measurement of circulating NT-proBNP levels. Results Predictors of NT-proBNP included age, female sex, body mass index, and cardiorenal parameters (diastolic dysfunction [DID] severity; LV mass and left atrial volume; right ventricular overload; decreasing ejection fraction [EF] and creatinine clearance). The performance of NT-proBNP to detect any degree of LV dysfunction, including mild DID, was poor (area under the curve 0.56-0.66). In contrast, the performance of NT-proBNP for the detection of EF 0.90 regardless of age and sex; history of hypertension, diabetes, coronary artery disease; or body mass category. The ability of NT-proBNP to detect EF

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.