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.

Ultracompact dwarf galaxies in the Fornax Cluster

By utilizing the large multiplexing advantage of the Two-degree Field spectrograph on the Anglo-Australian Telescope, we have been able to obtain a complete spectroscopic sample of all objects in a predefined magnitude range, 16.5 < b(j) < 19.7 regardless of morphology, in an area toward the center of the Fornax Cluster of galaxies. Among the unresolved or marginally resolved targets, we have found five objects that are actually at the redshift of the Fornax Cluster; i.e., they are extremely compact dwarf galaxies or extremely large star clusters. All five have absorption-line spectra. With intrinsic sizes of less than 1.1 HWHM (corresponding to approximately 100 pc at the distance of the cluster), they are more compact and significantly less luminous than other known compact dwarf galaxies, yet much brighter than any globular cluster. In this paper we present new ground-based optical observations of these enigmatic objects. In addition to having extremely high central surface brightnesses, these objects show no evidence of any surrounding low surface brightness envelopes down to much fainter limits than is the case for, e.g., nucleated dwarf elliptical galaxies. Thus, if they are not merely the stripped remains of some other type of galaxy, then they appear to have properties unlike any previously known type of stellar system.

Detecting environmental impacts on metapopulations of mound spring invertebrates - Assessing an incidence function model

We use a stochastic patch occupancy model of invertebrates in the Mound Springs ecosystem of South Australia to assess the ability of incidence function models to detect environmental impacts on metapopulations. We assume that the probability of colonisation decreases with increasing isolation and the probability of extinction is constant across spring vents. We run the models to quasi-equilibrium, and then impose an impact by increasing the local extinction probability. We sample the output at various times pre- and postimpact, and examine the probability of detecting a significant change in population parameters. The incidence function model approach turns out to have little power to detect environmental impacts on metapopulations with small numbers of patches. (C) 2001 Elsevier Science Ltd. All rights reserved.

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.

Population pharmacokinetics of lamotrigine

The present study estimated the population pharmacokinetics of lamotrigine in patients receiving oral lamotrigine therapy with drug concentration monitoring, and determined intersubject and intrasubject variability. A total of 129 patients were analyzed from two clinical sites. Of these, 124 patients provided spare data (198 concentration-time points); nine patients (four from a previous group plus five from the current group) provided rich data (431 points). The population analysis was conducted using P-PHARM (TM) (SIMED Scientific Software, Cedex, France), a nonlinear mixed-effect modeling program. A single exponential elimination model (first-order absorption) with heteroscedastic weighting was used. Apparent clearance (CL/F) and volume of distribution (V/F) were the pharmacokinetic parameters estimated. Covariate analysis was performed to determine which factors explained any of the variability associated with lamotrigine clearance. Population estimates of CL/F and V/F for lamotrigine generated in the final model were 2.14 +/- 0.81 L/h and 78.1 +/- 5.1 L/kg. Intersubject and intrasubject variability for clearance was 38% and 38%, respectively. The covariates of concomitant valproate and phenytoin therapy accounted for 42% of the intersubject variability of clearance. Age, gender, clinic site, and other concomitant antiepileptic drugs did not influence clearance. This study of the population pharmacokinetics of lamotrigine in patients using the drug clinically provides useful data and should lead to better dosage individualization for lamotrigine.

Population pharmacokinetics of tacrolimus in children who receive cut-down or full liver transplants

Background. The aim of this study was to investigate the population pharmacokinetics of tacrolimus in pediatric liver transplant recipients and to identify factors that may explain pharmacokinetic variability. Methods. Data were collected retrospectively from 35 children who received oral immunosuppressant therapy with tacrolimus. Maximum likelihood estimates were sought for the typical values of apparent clearance (CL/F) and apparent volume of distribution (V/F) with the program NONMEM. Factors screened for influence on the pharmacokinetic parameters were weight, age, gender, postoperative day, days since commencing tacrolimus therapy, transplant type (whole child liver or cut-down adult liver), liver function tests (bilirubin, alkaline phosphatase [ALP], aspartate aminotransferase [AST], gamma -glutamyl transferase [GGT], alanine aminotransferase [ALT]), creatinine clearance, hematocrit, corticosteroid dose, and concurrent therapy with metabolic inducers and inhibitors of tacrolimus. Results. No clear correlation existed between tacrolimus dosage and blood concentrations (r(2) =0.003). Transplant type, age, and liver function test values were the most important factors (P

Mark-recapture may reveal more about ecology than about population trends: Demography of a threatened ghost bat (Macroderma gigas) population

Since European settlement in Australia, the geographical range of ghost bats (Macroderma gigas) has contracted northwards. Ghost bats are thought to occur in disjunct populations with little interpopulation migration, raising concerns over the current status and future viability of the southernmost colony, which has also been threatened by mining activity. To address these concerns, demographic parameters of the southernmost colony were estimated from a mark-recapture study conducted during 1975-1981. Female bats gave birth to a single young in late spring, but only 40% (22-70%, 95% CI) of females bred in their second year, increasing to 93% (87-97%, 95% CI) for females greater than or equal to 2 years old. Sixty-five percent of juveniles caught were female. Annual adult survival ranged between 0.57-0.77 for females and 0.43-0.66 for males, and was lowest over winter-spring and greatest in autumn-winter. Juvenile survival for the first year ranged between 0.35-0.46 for females and 0.29-0.42 for males. Adult survival varied among seasons, was negatively associated with rainfall, but was not associated with temperature beyond being lower in late winter. Poor survival may result from the inferior daytime roosts that bats must use if water seepage forces them to leave their normal roosts. Although these age-specific rates of fecundity and survival suggested a declining population, mark-recapture estimates of the population trend indicated stability over the study period. Counts at daytime roosts also suggested a population decline, but were considered unreliable because of an increasing tendency of bats to avoid detection. It is therefore likely that some assumptions in estimating survival were violated. These results provide a caution against the uncritical use of population projections derived from mark-recapture estimates of demographic parameters, and the use of untested indices as the basis for conservation decisions.

Inferring population history from microsatellite and enzyme data in serially introduced cane toads, Bufo marinus

Much progress has been made on inferring population history from molecular data. However, complex demographic scenarios have been considered rarely or have proved intractable. The serial introduction of the South-Central American cane Load Bufo marinas in various Caribbean and Pacific islands involves four major phases: a possible genetic admixture during the first introduction, a bottleneck associated with founding, a transitory, population boom, and finally, a demographic stabilization. A large amount of historical and demographic information is available for those introductions and can be combined profitably with molecular data. We used a Bayesian approach to combine this information With microsatellite (10 loci) and enzyme (22 loci) data and used a rejection algorithm to simultaneously estimate the demographic parameters describing the four major phases of the introduction history,. The general historical trends supported by microsatellites and enzymes were similar. However, there was a stronger support for a larger bottleneck at introductions for microsatellites than enzymes and for a more balanced genetic admixture for enzymes than for microsatellites. Verb, little information was obtained from either marker about the transitory population boom observed after each introduction. Possible explanations for differences in resolution of demographic events and discrepancies between results obtained with microsatellites and enzymes were explored. Limits Of Our model and method for the analysis of nonequilibrium populations were discussed.

Comparison of two population pharmacokinetic programs, NONMEM and P-PHARM, for tacrolimus

Objectives: To compare the population modelling programs NONMEM and P-PHARM during investigation of the pharmacokinetics of tacrolimus in paediatric liver-transplant recipients. Methods: Population pharmacokinetic analysis was performed using NONMEM and P-PHARM on retrospective data from 35 paediatric liver-transplant patients receiving tacrolimus therapy. The same data were presented to both programs. Maximum likelihood estimates were sought for apparent clearance (CL/F) and apparent volume of distribution (V/F). Covariates screened for influence on these parameters were weight, age, gender, post-operative day, days of tacrolimus therapy, transplant type, biliary reconstructive procedure, liver function tests, creatinine clearance, haematocrit, corticosteroid dose, and potential interacting drugs. Results: A satisfactory model was developed in both programs with a single categorical covariate - transplant type - providing stable parameter estimates and small, normally distributed (weighted) residuals. In NONMEM, the continuous covariates - age and liver function tests - improved modelling further. Mean parameter estimates were CL/F (whole liver) = 16.3 1/h, CL/F (cut-down liver) = 8.5 1/h and V/F = 565 1 in NONMEM, and CL/F = 8.3 1/h and V/F = 155 1 in P-PHARM. Individual Bayesian parameter estimates were CL/F (whole liver) = 17.9 +/- 8.8 1/h, CL/F (cutdown liver) = 11.6 +/- 18.8 1/h and V/F = 712 792 1 in NONMEM, and CL/F (whole liver) = 12.8 +/- 3.5 1/h, CL/F (cut-down liver) = 8.2 +/- 3.4 1/h and V/F = 221 1641 in P-PHARM. Marked interindividual kinetic variability (38-108%) and residual random error (approximately 3 ng/ml) were observed. P-PHARM was more user friendly and readily provided informative graphical presentation of results. NONMEM allowed a wider choice of errors for statistical modelling and coped better with complex covariate data sets. Conclusion: Results from parametric modelling programs can vary due to different algorithms employed to estimate parameters, alternative methods of covariate analysis and variations and limitations in the software itself.

Population pharmacokinetics of tacrolimus in adult kidney transplant recipients

Objectives: The aims of this study were to investigate the population pharmacokinetics of tacrolimus in adult kidney transplant recipients and to identify factors that explain variability. Methods: Population analysis was performed on retrospective data from 70 patients who received oral tacrolimus twice daily. Morning blood trough concentrations were measured by liquid chromatography-tandem mass spectrometry. Maximum likelihood estimates were sought for apparent clearance (CL/F) and apparent volume of distribution (V/F), with the use of NONMEM (GloboMax LLC, Hanover, Md). Factors screened for influence on these parameters were weight, age, gender, postoperative day, days of tacrolimus therapy, liver function tests, creatinine clearance, hematocrit fraction, corticosteroid dose, and potential interacting drugs. Results. CL/F was greater in patients with abnormally low hematocrit fraction (data from 21 patients only), and it decreased with increasing days of therapy and AST concentrations (P