Do habitat suitability models reliably predict the recovery areas of threatened species?

1. Identifying those areas suitable for recolonization by threatened species is essential to support efficient conservation policies. Habitat suitability models (HSM) predict species' potential distributions, but the quality of their predictions should be carefully assessed when the species-environment equilibrium assumption is violated.2. We studied the Eurasian otter Lutra lutra, whose numbers are recovering in southern Italy. To produce widely applicable results, we chose standard HSM procedures and looked for the models' capacities in predicting the suitability of a recolonization area. We used two fieldwork datasets: presence-only data, used in the Ecological Niche Factor Analyses (ENFA), and presence-absence data, used in a Generalized Linear Model (GLM). In addition to cross-validation, we independently evaluated the models with data from a recolonization event, providing presences on a previously unoccupied river.3. Three of the models successfully predicted the suitability of the recolonization area, but the GLM built with data before the recolonization disagreed with these predictions, missing the recolonized river's suitability and badly describing the otter's niche. Our results highlighted three points of relevance to modelling practices: (1) absences may prevent the models from correctly identifying areas suitable for a species spread; (2) the selection of variables may lead to randomness in the predictions; and (3) the Area Under Curve (AUC), a commonly used validation index, was not well suited to the evaluation of model quality, whereas the Boyce Index (CBI), based on presence data only, better highlighted the models' fit to the recolonization observations.4. For species with unstable spatial distributions, presence-only models may work better than presence-absence methods in making reliable predictions of suitable areas for expansion. An iterative modelling process, using new occurrences from each step of the species spread, may also help in progressively reducing errors.5. Synthesis and applications. Conservation plans depend on reliable models of the species' suitable habitats. In non-equilibrium situations, such as the case for threatened or invasive species, models could be affected negatively by the inclusion of absence data when predicting the areas of potential expansion. Presence-only methods will here provide a better basis for productive conservation management practices.

Influence of the morphology of the dural sac on surgical decision making in lumbar spinal stenosis

Introduction: Surgical decision making in lumbar spinal stenosis (LSS) takes into account primarily clinical symptoms as well as concordant radiological findings. We hypothesized that a wide variation of operative threshold would be found in particular as far as judgment of severity of radiological stenosis is concerned. Patients and methods: The number of surgeons who would proceed to decompression was studied relative to the perceived severity of radiological stenosis based either on measurements of dural sac cross sectional area (DSCA) or on the recently described morphological grading as seen on axial T2 MRI images. A link to an electronic survey page with a set of ten axial T2 MRI images taken from ten patients with either low back pain or LSS were sent to members of three national or international spine societies. Those 10 images were randomly presented initially and re-shuffled on a second page including this time DSCA measurements in mm2, ranging from 14 to 226 mm2, giving a total of 20 images to appraise. Morphological grades were ranging from grade A to D. Surgeons were asked if they would consider decompression given the radiological appearance of stenosis and that symptoms of neurological claudication were severe in patients who were otherwise fit for surgery. Fisher's exact test was performed following dichotomization of data when appropriate. Results: A total of 142 spine surgeons (113 orthopedic spine surgeons, 29 neurosurgeons) responded from 25 countries. A substantial agreement was observed in operating patients with severe (grade C) or extreme (grade D) stenosis as defined by the morphological grade compared to lesser stenosis (A&B) grades (p<0.0001). Decision to operate was not dependent on number of years in practice, medical density in practicing country or specialty although more neurosurgeons would operate on grade C stenosis (p<0.005). Disclosing the DSCA measurement did not alter the decision to operate. Although 20 surgeons only had prior knowledge of the description of the morphological grading, their responses showed no statistically significant difference with those of the remaining 122 physicians. Conclusions: This study showed that surgeons across borders are less influenced by DSCA in their decision making than by the morphological appearance of the dural sac. Classifying LSS according to morphology rather than surface measurements appears to be consistent with current clinical practice.

Management and prognosis of status epilepticus according to hospital setting: a prospective study.

BACKGROUND: The treatment of status epilepticus (SE) is based on relatively little evidence although several guidelines have been published. A recent study reported a worse SE prognosis in a large urban setting as compared to a peripheral hospital, postulating better management in the latter. The aim of this study was to analyse SE episodes occurring in different settings and address possible explanatory variables regarding outcome, including treatment quality. METHODS: Over six months we prospectively recorded consecutive adults with SE (fit lasting five or more minutes) at the Centre Hospitalier Universitaire Vaudois (CHUV) and in six peripheral hospitals (PH) in the same region. Demographical, historical and clinical variables were collected, including SE severity estimation (STESS score) and adherence to Swiss SE treatment guidelines. Outcome at discharge was categorised as "good" (return to baseline), or "poor" (persistent neurological sequelae or death). RESULTS: Of 54 patients (CHUV: 36; PH 18), 33% had a poor outcome. Whilst age, SE severity, percentage of SE episodes lasting less than 30 minutes and total SE duration were similar, fewer patients had a good outcome at the CHUV (61% vs 83%; OR 3.57; 95% CI 0.8-22.1). Mortality was 14% at the CHUV and 5% at the PH. Most treatments were in agreement with national guidelines, although less often in PH (78% vs 97%, P = 0.04). CONCLUSION: Although not statistically significant, we observed a slightly worse SE prognosis in a large academic centre as compared to smaller hospitals. Since SE severity was similar in the two settings but adherence to national treatment guidelines was higher in the academic centre, further investigation on the prognostic role of SE treatment and outcome determinants is required.

Modelling the influence of change in fire regime on the local distribution of a Mediterranean pyrophytic plant species (Cistus salviifolius) at its northern range limit

Aims: To assess the potential distribution of an obligate seeder and active pyrophyte, Cistus salviifolius, a vulnerable species in the Swiss Red List; to derive scenarios by changing the fire return interval; and to discuss the results from a conservation perspective. A more general aim is to assess the impact of fire as a natural factor influencing the vegetation of the southern slopes of the Alps. Locations: Alps, southern Switzerland. Methods: Presence-absence data to fit the model were obtained from the most recent field mapping of C. salviifolius. The quantitative environmental predictors used in this study include topographic, climatic and disturbance (fire) predictors. Models were fitted by logistic regression and evaluated by jackknife and bootstrap approaches. Changes in fire regime were simulated by increasing the time-return interval of fire (simulating longer periods without fire). Two scenarios were considered: no fire in the past 15 years; or in the past 35 years. Results: Rock cover, slope, topographic position, potential evapotranspiration and time elapsed since the last fire were selected in the final model. The Nagelkerke R-2 of the model for C. salviifolius was 0.57 and the Jackknife area under the curve evaluation was 0.89. The bootstrap evaluation revealed model robustness. By increasing the return interval of fire by either up to 15 years, or 35 years, the modelled C. salviifolius population declined by 30-40%, respectively. Main conclusions: Although fire plays a significant role, topography and rock cover appear to be the most important predictors, suggesting that the distribution of C. salviifolius in the southern Swiss Alps is closely related to the availability of supposedly competition-free sites, such as emerging bedrock, ridge locations or steep slopes. Fire is more likely to play a secondary role in allowing C. salviifolius to extend its occurrence temporarily, by increasing germination rates and reducing the competition from surrounding vegetation. To maintain a viable dormant seed bank for C. salviifolius, conservation managers should consider carrying out vegetation clearing and managing wild fire propagation to reduce competition and ensure sufficient recruitment for this species.

Lack of Enantioselectivity in the SULT1A3-catalyzed Sulfoconjugation of Normetanephrine Enantiomers: An In Vitro and Computational Study.

(1R)-Normetanephrine is the natural stereoisomeric substrate for sulfotransferase 1A3 (SULT1A3)-catalyzed sulfonation. Nothing appears known on the enantioselectivity of the reaction despite its potential significance in the metabolism of adrenergic amines and in clinical biochemistry. We confronted the kinetic parameters of the sulfoconjugation of synthetic (1R)-normetanephrine and (1S)-normetanephrine by recombinant human SULT1A3 to a docking model of each normetanephrine enantiomer with SULT1A3 and the 3'-phosphoadenosine-5'-phosphosulfate cofactor on the basis of molecular modeling and molecular dynamics simulations of the stability of the complexes. The K(M) , V(max) , and k(cat) values for the sulfonation of (1R)-normetanephrine, (1S)-normetanephrine, and racemic normetanephrine were similar. In silico models were consistent with these findings as they showed that the binding modes of the two enantiomers were almost identical. In conclusion, SULT1A3 is not substrate-enantioselective toward normetanephrine, an unexpected finding explainable by a mutual adaptability between the ligands and SULT1A3 through an "induced-fit model" in the catalytic pocket. Chirality, 00:000-000, 2012.© 2012 Wiley Periodicals, Inc.

Long-Term Safety of Canakinumab in Patients with Gouty Arthritis.

Background/Purpose: Gouty arthritis (GA) is a chronic inflammatory disease. Targeting the inflammatory pathway through IL-1_ inhibition with canakinumab (CAN) may provide significant long-term benefits. CAN safety versus triamcinolone acetonide (TA) over initial 24 weeks (blinded study) for patients (pts) with history of frequent attacks (_3 in year before baseline) was reported earlier from core (_-RELIEVED [_-REL] and _-REL-II) and first extension (E1) studies1. Herein we present full 18-month long-term CAN safety data, including open-label second extension (E2) studies. Methods: GA pts completing _-REL E1 and _-REL-II E1 studies1 were enrolled in these 1-year, open-label, E2 studies. All pts entering E2, whether randomized to CAN or TA, received CAN 150 mg sc on demand upon new attack. Data are presented only for pts randomized to CAN, and are reported cumulatively, i.e. including corresponding data from previously reported core and E1 studies. Long-term safety outcomes and safety upon re-treatment are presented as incidence rate per 100 patient-years (pyr) of study participation for AEs and SAEs. Deaths are reported for all pts (randomized to CAN or TA). Selected predefined notable laboratory abnormalities are shown (neutrophils, platelets, liver and renal function tests). Long-term attack rate per year is also provided. Results: In total, 69/115 (60%) and 72/112 (64.3%) of the pts randomized to CAN in the two core studies entered the two E2 studies, of which 68 and 64 pts, respectively completed the E2 studies. The 2 study populations had differing baseline comorbidity and geographic origin. Lab data (not time adjusted) for neutropenia appears worse after retreatment in _-REL E2, and deterioration of creatinine clearance appears worse after retreatment (Table 1). The time-adjusted incidence rates for AEs were 302.4/100 pyr and 360/100 pyr, and for SAEs were 27.9/100 pyr and 13.9/100 pyr in _-REL E2 and _-REL-II E2 respectively (Table 1). The time-adjusted incidence rates of any AEs, infection AEs, any SAEs, and selected SAEs before and after re-treatment are presented in Table 1. Incidence rates for AEs and SAEs declined after re-treatment, with the exception of SAEs in _-REL-II E2, which increased from 2.9/100 pyr to 10.9/100 pyr (no infection SAEs after retreatment in _-REL-II E2, and other SAEs fit no special pattern). In the total safety population (N_454, core and all extensions), there were 4 deaths, 2 in the core studies previously reported1 and 2 during the _-REL E2 study (one patient in the CAN group died from pneumonia; one patient in the TA group who never received CAN died of pneumococcal sepsis). None of the deaths was suspected by investigators to be study drug related. The mean rates of new attacks per year on CAN were 1.21 and 1.18 in _-REL E2 and in _-REL-II E2. Conclusion: The clinical safety profile of CAN upon re-treatment was maintained long-term with no new infection concerns

Economics of dialysis dependence following renal replacement therapy for critically ill acute kidney injury patients.

BACKGROUND: The obective of this study was to perform a cost-effectiveness analysis comparing intermittent with continuous renal replacement therapy (IRRT versus CRRT) as initial therapy for acute kidney injury (AKI) in the intensive care unit (ICU). METHODS: Assuming some patients would potentially be eligible for either modality, we modeled life year gained, the quality-adjusted life years (QALYs) and healthcare costs for a cohort of 1000 IRRT patients and a cohort of 1000 CRRT patients. We used a 1-year, 5-year and a lifetime horizon. A Markov model with two health states for AKI survivors was designed: dialysis dependence and dialysis independence. We applied Weibull regression from published estimates to fit survival curves for CRRT and IRRT patients and to fit the proportion of dialysis dependence among CRRT and IRRT survivors. We then applied a risk ratio reported in a large retrospective cohort study to the fitted CRRT estimates in order to determine the proportion of dialysis dependence for IRRT survivors. We conducted sensitivity analyses based on a range of differences for daily implementation cost between CRRT and IRRT (base case: CRRT day $632 more expensive than IRRT day; range from $200 to $1000) and a range of risk ratios for dialysis dependence for CRRT as compared with IRRT (from 0.65 to 0.95; base case: 0.80). RESULTS: Continuous renal replacement therapy was associated with a marginally greater gain in QALY as compared with IRRT (1.093 versus 1.078). Despite higher upfront costs for CRRT in the ICU ($4046 for CRRT versus $1423 for IRRT in average), the 5-year total cost including the cost of dialysis dependence was lower for CRRT ($37 780 for CRRT versus $39 448 for IRRT on average). The base case incremental cost-effectiveness analysis showed that CRRT dominated IRRT. This dominance was confirmed by extensive sensitivity analysis. CONCLUSIONS: Initial CRRT is cost-effective compared with initial IRRT by reducing the rate of long-term dialysis dependence among critically ill AKI survivors.

Comparative performance of supertree algorithms in large data sets using the soapberry family (Sapindaceae) as a case study.

For the last 2 decades, supertree reconstruction has been an active field of research and has seen the development of a large number of major algorithms. Because of the growing popularity of the supertree methods, it has become necessary to evaluate the performance of these algorithms to determine which are the best options (especially with regard to the supermatrix approach that is widely used). In this study, seven of the most commonly used supertree methods are investigated by using a large empirical data set (in terms of number of taxa and molecular markers) from the worldwide flowering plant family Sapindaceae. Supertree methods were evaluated using several criteria: similarity of the supertrees with the input trees, similarity between the supertrees and the total evidence tree, level of resolution of the supertree and computational time required by the algorithm. Additional analyses were also conducted on a reduced data set to test if the performance levels were affected by the heuristic searches rather than the algorithms themselves. Based on our results, two main groups of supertree methods were identified: on one hand, the matrix representation with parsimony (MRP), MinFlip, and MinCut methods performed well according to our criteria, whereas the average consensus, split fit, and most similar supertree methods showed a poorer performance or at least did not behave the same way as the total evidence tree. Results for the super distance matrix, that is, the most recent approach tested here, were promising with at least one derived method performing as well as MRP, MinFlip, and MinCut. The output of each method was only slightly improved when applied to the reduced data set, suggesting a correct behavior of the heuristic searches and a relatively low sensitivity of the algorithms to data set sizes and missing data. Results also showed that the MRP analyses could reach a high level of quality even when using a simple heuristic search strategy, with the exception of MRP with Purvis coding scheme and reversible parsimony. The future of supertrees lies in the implementation of a standardized heuristic search for all methods and the increase in computing power to handle large data sets. The latter would prove to be particularly useful for promising approaches such as the maximum quartet fit method that yet requires substantial computing power.

Mechanisms involved in the maintenance of inbreeding depression in gynodioecious "Silene vulgaris" (Caryophyllaceae): an experimental investigation

Summary Gynodioecy, the joint occurrence of females and hermaphrodites within natural populations, is a widely studied mating system ever since Darwin (1877). It is an exceptional mating system because continuous selection is necessary to maintain it. Since females only reproduce through ovules whereas hermaphrodites transmit genes through ovules and pollen, larger female fitness, in terms of seed output, is required to allow their maintenance. Two non-exclusive mechanisms can account for the maintenance of females. First, as females do not produce pollen they can reallocate their resources towards a higher ovule production. Second, hermaphrodites can self- and cross-fertilize whereas females are obligate outcrossers. Thus hermaphrodites should partly suffer from inbreeding depression (i.e.: the fitness decline of inbred relative to outbred individuals) and thereby produce less fit progeny than females. This thesis investigated the effects of self- and cross-fertilization of heimaphrodites over two consecutive generations. Inbreeding depression increased across the successive stages of the life- cycle (i.e.: from "seed traits" to "reproductive traits") displaying large inbreeding depression estimates (up to 0.76). This investigation not only detected large inbreeding depression estimates but also detected mechanisms involved in the maintenance of inbreeding depression. For instance cryptic self-incompatibility which is here a larger in vivo pollen performance of distant pollen compared to self-pollen; the expression of inbreeding depression especially in late life-cycle stages, and the appearance of females in the progeny of selfed hermaphrodites. The female biased sex ratio in the progeny of selfed hermaphrodites was a surprising result and could either come from the sex determining mechanisms (complex nucleo-cytoplasmic interaction(s)) and/or from inbreeding depression. Indeed, we not only got females and hermaphrodites but also partial male-sterile (PMS) individuals (i.e.: individuals with differing number of viable stamens). We detected that inbred pollen bearing plants (excluding females) have less viable stamens per flower than outbred plants. A positive correlation was detected between inbreeding depression for the number of viable stamens per flower and the difference in sex ratio between inbred and outbred individuals. A positive relationship was also detected between inbreeding depression for pollen viability and inbreeding depression for number of viable stamens per flower. Each correlation can either account for pleiotropic effects (a major gene acting on the two considered traits) or linkage disequilibrium between genes controlling each of the two related traits. If we hypothesize that these correlations are due to a major gene with pleiotropic effects, the positive relationship between inbreeding depression for number of viable stamens per flower and inbreeding depression for pollen viability showed that deleterious alleles present on a major gene coding for pollen production and viability depressed male fitness within inbred plants. The positive relationship between sex ratio difference between inbred and outbred individuals and inbreeding depression for number of viable stamens per flower indicates that (1) either number of viable stamens per flower is, in addition to inbreeding, also affected by the loci coding for sex determinism or, (2) the presence of females within the progeny of selfed hermaphrodites is a consequence of large inbreeding depression inhibiting pollen production, or (3) sex is here determined by a combination of loci coding for sex expression and inbreeding depression for male reproductive traits. In conclusion, Silene vulgaris has been shown to be a good model for understanding the evolution of mating systems that promote outbreeding. Résumé La gynodïoécie est définie comme étant la présence simultanée d'hermaphrodites et de femelles au sein de populations naturelles d'une même espèce. Ce système de reproduction a toujours fasciné le monde scientifique depuis Darwin, comme en témoigne ses écrits (1876, 1877) sur les systèmes de reproduction chez les plantes. Les femelles ne transmettent leurs gènes qu'à travers leurs ovules alors que les hermaphrodites transmettent leurs gènes à la fois par la voie mâle (le pollen) et la voie femelle (les ovules). La condition pour que la gynodïoécie se maintienne nécessite donc une fitness de la fonction femelle plus élevée chez les femelles que chez les hermaphrodites. Deux mécanismes mutuellement non exclusifs peuvent expliquer le maintien des femelles au sein de ces populations gynodioïques. D'une part, les femelles peuvent réallouer les ressources non utilisées pour la production de pollen et peuvent par conséquent produire plus d'ovules. D'autre part, la reproduction des femelles ne peut se faire que par allo-fécondation alors que les hermaphrodites, peuvent se reproduire à la fois par auto- et allo-fécondation. L'autofécondation s'accompagne en général d'une diminution de fitness de la descendance relativement à la progéniture issue d'allo-fécondation ; ce phénomène est connu sous le nom de dépression de consanguinité. Cette thèse avait pour but de mettre en évidence une éventuelle dépression de consanguinité chez Silene vulgaris, une espèce gynodioïque. Des hermaphrodites, issus de trois vallées alpines, ont été auto- et allo¬fécondés sur deux générations successives. La dépression de consanguinité pouvant s'exprimer à tous les stades de vie d'un individu, plusieurs traits de fitness, allant du nombre de graines par fruit à la production de gamètes ont été mesurés sur différents stades de vie successifs. L'estimation de la dépression de consanguinité totale atteignait des valeurs allant de 0.52 à 0.76 selon la vallée considérée, ce qui indiquerait que les hermaphrodites ont tout intérêt à limiter l'autofécondation et que les femelles ne devraient pas avoir de peine à subsister dans les vallées étudiées. Par la même occasion des mécanismes diminuant la purge potentielle du fardeau génétique, et permettant ainsi le maintien du « niveau » de dépression de consanguinité et par conséquence le maintien de la gynodïoécie ont été mis en évidence. En effet, nos résultats montrent que la dépression de consanguinité s'exprimait tard dans le cycle de vie permettant ainsi à un certain nombre individus consanguins de transmettre leurs allèles délétères à la génération suivante. D'autre part, la croissance in vivo des tubes polliniques d'auto-pollen était plus lente que celle de l'allo-pollen et donc en situation de compétition directe, les ovules devraient plutôt être issus d'allo-fécondation, diminuant ainsi les chances de purges d'allèles délétères. Enfin, l'apparition de femelles dans la progéniture d'hermaphrodites autofécondés diminue aussi les chances de purge d'allèles délétères. Il nous a été impossible de déterminer si l'apparition de femelles dans la descendance d'hermaphrodites autofécondés était due au déterminisme génétique du sexe ou si la différence de sexe ratio entre la descendance auto- et allo-fécondée était due à une éventuelle dépression de consanguinité inhibant la production de pollen. Nous avons observé que S. vulgaris ne présentaient pas uniquement des hermaphrodites et des femelles mais aussi toute sorte d'individus intermédiaires avec un nombre variable d'étamines viables. Nous avons pu mettre' en évidence des corrélations positives entre (1) la différence de sexe ratio (la proportion d'individus produisant du pollen) entre individus consanguins et non consanguins et une estimation de la dépression de consanguinité pour le nombre d'étamines viables d'individus produisant du pollen, ainsi qu'entre (2) la dépression de consanguinité pour le nombre d'étamines viables et celle estimée pour la viabilité du pollen. Chaque corrélation indique soit l'effet d'un (ou plusieurs) gène(s) pléiotropique(s), soit un déséquilibre de liaison entre les gènes. En considérant que ces corrélations sont le résultat d'effet pléiotropiques, la relation entre le nombre d'étamines viables par fleur et la viabilité du pollen, indiquerait un effet négatif de la consanguinité sur la production et la viabilité du pollen due partiellement à un gène majeur. La seconde corrélation indiquerait soit que les gènes responsables de la détermination du sexe agissent aussi sur l'expression de la fonction mâle soit que l'expression du sexe est sujette à la dépression de consanguinité, ou encore un mélange des deux. Aux regards de ces résultats, Silene vulgaris s'est avéré être un bon modèle de compréhension de l'évolution des systèmes de reproduction vers la séparation des sexes.

Benchmarking therapeutic drug monitoring software: A systematic evaluation of available computer tools

Introduction: Therapeutic drug monitoring (TDM) aims at optimizing treatment by individualizing dosage regimen based on measurement of blood concentrations. Maintaining concentrations within a target range requires pharmacokinetic and clinical capabilities. Bayesian calculation represents a gold standard in TDM approach but requires computing assistance. In the last decades computer programs have been developed to assist clinicians in this assignment. The aim of this benchmarking was to assess and compare computer tools designed to support TDM clinical activities.¦Method: Literature and Internet search was performed to identify software. All programs were tested on common personal computer. Each program was scored against a standardized grid covering pharmacokinetic relevance, user-friendliness, computing aspects, interfacing, and storage. A weighting factor was applied to each criterion of the grid to consider its relative importance. To assess the robustness of the software, six representative clinical vignettes were also processed through all of them.¦Results: 12 software tools were identified, tested and ranked. It represents a comprehensive review of the available software's characteristics. Numbers of drugs handled vary widely and 8 programs offer the ability to the user to add its own drug model. 10 computer programs are able to compute Bayesian dosage adaptation based on a blood concentration (a posteriori adjustment) while 9 are also able to suggest a priori dosage regimen (prior to any blood concentration measurement), based on individual patient covariates, such as age, gender, weight. Among those applying Bayesian analysis, one uses the non-parametric approach. The top 2 software emerging from this benchmark are MwPharm and TCIWorks. Other programs evaluated have also a good potential but are less sophisticated (e.g. in terms of storage or report generation) or less user-friendly.¦Conclusion: Whereas 2 integrated programs are at the top of the ranked listed, such complex tools would possibly not fit all institutions, and each software tool must be regarded with respect to individual needs of hospitals or clinicians. Interest in computing tool to support therapeutic monitoring is still growing. Although developers put efforts into it the last years, there is still room for improvement, especially in terms of institutional information system interfacing, user-friendliness, capacity of data storage and report generation.

Benchmarking therapeutic drug monitoring software: A systematic evaluation of available computer tools

Objectives: Therapeutic drug monitoring (TDM) aims at optimizing treatment by individualizing dosage regimen based on blood concentrations measurement. Maintaining concentrations within a target range requires pharmacokinetic (PK) and clinical capabilities. Bayesian calculation represents a gold standard in TDM approach but requires computing assistance. The aim of this benchmarking was to assess and compare computer tools designed to support TDM clinical activities.¦Methods: Literature and Internet were searched to identify software. Each program was scored against a standardized grid covering pharmacokinetic relevance, user-friendliness, computing aspects, interfacing, and storage. A weighting factor was applied to each criterion of the grid to consider its relative importance. To assess the robustness of the software, six representative clinical vignettes were also processed through all of them.¦Results: 12 software tools were identified, tested and ranked. It represents a comprehensive review of the available software characteristics. Numbers of drugs handled vary from 2 to more than 180, and integration of different population types is available for some programs. Nevertheless, 8 programs offer the ability to add new drug models based on population PK data. 10 computer tools incorporate Bayesian computation to predict dosage regimen (individual parameters are calculated based on population PK models). All of them are able to compute Bayesian a posteriori dosage adaptation based on a blood concentration while 9 are also able to suggest a priori dosage regimen, only based on individual patient covariates. Among those applying Bayesian analysis, MM-USC*PACK uses a non-parametric approach. The top 2 programs emerging from this benchmark are MwPharm and TCIWorks. Others programs evaluated have also a good potential but are less sophisticated or less user-friendly.¦Conclusions: Whereas 2 software packages are ranked at the top of the list, such complex tools would possibly not fit all institutions, and each program must be regarded with respect to individual needs of hospitals or clinicians. Programs should be easy and fast for routine activities, including for non-experienced users. Although interest in TDM tools is growing and efforts were put into it in the last years, there is still room for improvement, especially in terms of institutional information system interfacing, user-friendliness, capability of data storage and automated report generation.

Measurements of glial metabolic fluxes with C-11-acetate using positron emission and an adapted NMR-based metabolic modeling approach

Objectives: Acetate brain metabolism has the particularity to occur specifically in glial cells. Labeling studies, using acetate labeled either with 13C (NMR) or 11C (PET), are governed by the same biochemical reactions and thus follow the same mathematical principles. In this study, the objective was to adapt an NMR acetate brain metabolism model to analyse [1-11C]acetate infusion in rats. Methods: Brain acetate infusion experiments were modeled using a two-compartment model approach used in NMR.1-3 The [1-11C]acetate labeling study was done using a beta scintillator.4 The measured radioactive signal represents the time evolution of the sum of all labeled metabolites in the brain. Using a coincidence counter in parallel, an arterial input curve was measured. The 11C at position C-1 of acetate is metabolized in the first turn of the TCA cycle to the position 5 of glutamate (Figure 1A). Through the neurotransmission process, it is further transported to the position 5 of glutamine and the position 5 of neuronal glutamate. After the second turn of the TCA cycle, tracer from [1-11C]acetate (and also a part from glial [5-11C]glutamate) is transferred to glial [1-11C]glutamate and further to [1-11C]glutamine and neuronal glutamate through the neurotransmission cycle. Brain poster session: oxidative mechanisms S460 Journal of Cerebral Blood Flow & Metabolism (2009) 29, S455-S466 Results: The standard acetate two-pool PET model describes the system by a plasma pool and a tissue pool linked by rate constants. Experimental data are not fully described with only one tissue compartment (Figure 1B). The modified NMR model was fitted successfully to tissue time-activity curves from 6 single animals, by varying the glial mitochondrial fluxes and the neurotransmission flux Vnt. A glial composite rate constant Kgtg=Vgtg/[Ace]plasma was extracted. Considering an average acetate concentration in plasma of 1 mmol/g5 and the negligible additional amount injected, we found an average Vgtg = 0.08±0.02 (n = 6), in agreement with previous NMR measurements.1 The tissue time-activity curve is dominated by glial glutamate and later by glutamine (Figure 1B). Labeling of neuronal pools has a low influence, at least for the 20 mins of beta-probe acquisition. Based on the high diffusivity of CO2 across the blood-brain barrier; 11CO2 is not predominant in the total tissue curve, even if the brain CO2 pool is big compared with other metabolites, due to its strong dilution through unlabeled CO2 from neuronal metabolism and diffusion from plasma. Conclusion: The two-compartment model presented here is also able to fit data of positron emission experiments and to extract specific glial metabolic fluxes. 11C-labeled acetate presents an alternative for faster measurements of glial oxidative metabolism compared to NMR, potentially applicable to human PET imaging. However, to quantify the relative value of the TCA cycle flux compared to the transmitochondrial flux, the chemical sensitivity of NMR is required. PET and NMR are thus complementary.

Efficacy, tolerability and risk factors for virological failure of darunavir-based therapy for treatment-experienced HIV-infected patients: the Swiss HIV Cohort Study.

OBJECTIVES: Darunavir was designed for activity against HIV resistant to other protease inhibitors (PIs). We assessed the efficacy, tolerability and risk factors for virological failure of darunavir for treatment-experienced patients seen in clinical practice. METHODS: We included all patients in the Swiss HIV Cohort Study starting darunavir after recording a viral load above 1000 HIV-1 RNA copies/mL given prior exposure to both PIs and nonnucleoside reverse transcriptase inhibitors. We followed these patients for up to 72 weeks, assessed virological failure using different loss of virological response algorithms and evaluated risk factors for virological failure using a Bayesian method to fit discrete Cox proportional hazard models. RESULTS: Among 130 treatment-experienced patients starting darunavir, the median age was 47 years, the median duration of HIV infection was 16 years, and 82% received mono or dual antiretroviral therapy before starting highly active antiretroviral therapy. During a median patient follow-up period of 45 weeks, 17% of patients stopped taking darunavir after a median exposure of 20 weeks. In patients followed beyond 48 weeks, the rate of virological failure at 48 weeks was at most 20%. Virological failure was more likely where patients had previously failed on both amprenavir and saquinavir and as the number of previously failed PI regimens increased. CONCLUSIONS: As a component of therapy for treatment-experienced patients, darunavir can achieve a similar efficacy and tolerability in clinical practice to that seen in clinical trials. Clinicians should consider whether a patient has failed on both amprenavir and saquinavir and the number of failed PI regimens before prescribing darunavir.

Why is the cladoceran Simocephalus vetulus (Müller) not a 'bang-bang strategist '? A critique of the optimal-body-size model

Lynch's (1980a) optimal-body-size model is designed to explain some major trends in cladoceran life histories; in particular the fact that large and littoral species seem to be bang-bang strategists (they grow first and the reproduce) whereas smaller planktonic species seem to be intermediate strategists (they grow and reproduce simultaneously). Predation is assumed to be an important selective pressure for these trends. Simocephalus vetulus (Müller) does not fit this pattern; being a littoral and relatively large species but an intermediate strategist. As shown by computer simulations, this species would reduce its per capita rate of increase by adopting the strategy predicted by the optimal-body-size model. Two aspects of the model are criticized: (1) the optimization criterion is shown to be incorrect and (2) the prediction of an intermediate strategy is not justified. Structural constraints are suggested to be responsible for the intermediate strategy of S.vetulus. Biotic interactions seem to have little effect on the observed life-history patterns of this species.

Full-waveform inversion of cross-hole ground-penetrating radar data to characterize a gravel aquifer close to the Thur River, Switzerland

Cross-hole radar tomography is a useful tool for mapping shallow subsurface electrical properties viz. dielectric permittivity and electrical conductivity. Common practice is to invert cross-hole radar data with ray-based tomographic algorithms using first arrival traveltimes and first cycle amplitudes. However, the resolution of conventional standard ray-based inversion schemes for cross-hole ground-penetrating radar (GPR) is limited because only a fraction of the information contained in the radar data is used. The resolution can be improved significantly by using a full-waveform inversion that considers the entire waveform, or significant parts thereof. A recently developed 2D time-domain vectorial full-waveform crosshole radar inversion code has been modified in the present study by allowing optimized acquisition setups that reduce the acquisition time and computational costs significantly. This is achieved by minimizing the number of transmitter points and maximizing the number of receiver positions. The improved algorithm was employed to invert cross-hole GPR data acquired within a gravel aquifer (4-10 m depth) in the Thur valley, Switzerland. The simulated traces of the final model obtained by the full-waveform inversion fit the observed traces very well in the lower part of the section and reasonably well in the upper part of the section. Compared to the ray-based inversion, the results from the full-waveform inversion show significantly higher resolution images. At either side, 2.5 m distance away from the cross-hole plane, borehole logs were acquired. There is a good correspondence between the conductivity tomograms and the natural gamma logs at the boundary of the gravel layer and the underlying lacustrine clay deposits. Using existing petrophysical models, the inversion results and neutron-neutron logs are converted to porosity. Without any additional calibration, the values obtained for the converted neutron-neutron logs and permittivity results are very close and similar vertical variations can be observed. The full-waveform inversion provides in both cases additional information about the subsurface. Due to the presence of the water table and associated refracted/reflected waves, the upper traces are not well fitted and the upper 2 m in the permittivity and conductivity tomograms are not reliably reconstructed because the unsaturated zone is not incorporated into the inversion domain.