HIERFSTAT, a package for R to compute and test hierarchical F-statistics

The package HIERFSTAT for the statistical software R, created by the R Development Core Team, allows the estimate of hierarchical F-statistics from a hierarchy with any numbers of levels. In addition, it allows testing the statistical significance of population differentiation for these different levels, using a generalized likelihood-ratio test. The package HIERFSTAT is available at http://www.unil.ch/popgen/softwares/hierfstat.htm.

Physiologically based pharmacokinetic modeling of arsenic in the mouse

A remarkable feature of the carcinogenicity of inorganic arsenic is that while human exposures to high concentrations of inorganic arsenic in drinking water are associated with increases in skin, lung, and bladder cancer, inorganic arsenic has not typically caused tumors in standard laboratory animal test protocols. Inorganic arsenic administered for periods of up to 2 yr to various strains of laboratory mice, including the Swiss CD-1, Swiss CR:NIH(S), C57Bl/6p53(+/-), and C57Bl/6p53(+/+), has not resulted in significant increases in tumor incidence. However, Ng et al. (1999) have reported a 40% tumor incidence in C57Bl/6J mice exposed to arsenic in their drinking water throughout their lifetime, with no tumors reported in controls. In order to investigate the potential role of tissue dosimetry in differential susceptibility to arsenic carcinogenicity, a physiologically based pharmacokinetic (PBPK) model for inorganic arsenic in the rat, hamster, monkey, and human (Mann et al., 1996a, 1996b) was extended to describe the kinetics in the mouse. The PBPK model was parameterized in the mouse using published data from acute exposures of B6C3F1 mice to arsenate, arsenite, monomethylarsonic acid (MMA), and dimethylarsinic acid (DMA) and validated using data from acute exposures of C57Black mice. Predictions of the acute model were then compared with data from chronic exposures. There was no evidence of changes in the apparent volume of distribution or in the tissue-plasma concentration ratios between acute and chronic exposure that might support the possibility of inducible arsenite efflux. The PBPK model was also used to project tissue dosimetry in the C57Bl/6J study, in comparison with tissue levels in studies having shorter duration but higher arsenic treatment concentrations. The model evaluation indicates that pharmacokinetic factors do not provide an explanation for the difference in outcomes across the various mouse bioassays. Other possible explanations may relate to strain-specific differences, or to the different durations of dosing in each of the mouse studies, given the evidence that inorganic arsenic is likely to be active in the later stages of the carcinogenic process. [Authors]

Altered vision during motion: an unusual symptom of cerebellar dysfunction, quantifiable by a simple clinical test.

PURPOSE: To report a series of patients with cerebellar dysfunction and altered vision during motion, and to quantify their visual impairment in motion with a simple clinical test. METHODS: Twenty consecutive patients suffering from cerebellar dysfunction and altered vision during motion were examined between 1994 and 2007. A control group consisted of 20 age- and sex-matched healthy people. All patients had a full neuro-ophthalmic examination. Near visual acuity (NVA) was measured at rest (static NVA) and during chair rotation (dynamic NVA). Distance visual acuity (DVA) was measured at rest (static DVA) and during rotation of the patient's head (dynamic DVA). RESULTS: Only four of the 20 patients reported altered vision during motion spontaneously. The remaining 16 patients admitted this unusual visual disturbance only when asked specifically. All patients exhibited abnormal eye movements, including saccadic smooth pursuit (20/20), dysmetric saccades (15/20), nystagmus (19/20) and impaired suppression of vestibulo-ocular reflex (VOR) (20/20). During rotation of the examination chair (dynamic NVA), the drop in NVA averaged 5.6 lines (range 1-10 lines). During rotation of the patient's head (dynamic DVA), the drop in DVA averaged only 2.5 lines (range 0-10 lines). For the control group, there was no significant drop in NVA under dynamic conditions. CONCLUSION: Patients with cerebellar dysfunction rarely complain spontaneously of altered vision during motion. However, specific questioning may bring up this unusual symptom. The use of a simple clinical test, consisting of NVA measurement during rotation of the examination chair (dynamic NVA), allows practitioners to quantify the level of visual impairment in patients presenting altered VOR modulation.

Dosimetric aspects of a national survey of diagnostic and interventional radiology in Switzerland.

The effective dose delivered to the patient was determined, by modeling, for 257 types of examinations covering the different modalities of diagnostic and interventional radiology. The basic operational dosimetric quantities considered were obtained from the parameters of the examinations on the basis of dosimetric models. These models required a precise characterization of each examination. The operational dosimetric quantities were converted into doses to organs and effective doses using appropriate conversion factors. The determination of the collective effective dose to the Swiss population requires a number of corrections to account for the variability of several parameters: sensitivity of the detection system, age, gender, and build of the patient. The use of various dosimetric models is illustrated in this paper for a limited number of examination types covering the different radiological modalities, for which the established typical effective doses are given. With regard to individual doses, the study indicated that the average effective doses per type of examination can be classified into three levels: (a) the weakly irradiating examinations (less than 0.1 mSv), which represent 78% of the examinations and 4% of the collective dose, (b) the moderately irradiating examinations (between 0.1 mSv and 10 mSv), which represent 21% of the examinations and 72% of the collective dose, (c) the strongly irradiating examinations (more than 10 mSv), which represent 1% of the examinations and 24% of the collective dose.

Determinants of brain cell metabolic phenotypes and energy substrate utilization unraveled with a modeling approach.

Although all brain cells bear in principle a comparable potential in terms of energetics, in reality they exhibit different metabolic profiles. The specific biochemical characteristics explaining such disparities and their relative importance are largely unknown. Using a modeling approach, we show that modifying the kinetic parameters of pyruvate dehydrogenase and mitochondrial NADH shuttling within a realistic interval can yield a striking switch in lactate flux direction. In this context, cells having essentially an oxidative profile exhibit pronounced extracellular lactate uptake and consumption. However, they can be turned into cells with prominent aerobic glycolysis by selectively reducing the aforementioned parameters. In the case of primarily oxidative cells, we also examined the role of glycolysis and lactate transport in providing pyruvate to mitochondria in order to sustain oxidative phosphorylation. The results show that changes in lactate transport capacity and extracellular lactate concentration within the range described experimentally can sustain enhanced oxidative metabolism upon activation. Such a demonstration provides key elements to understand why certain brain cell types constitutively adopt a particular metabolic profile and how specific features can be altered under different physiological and pathological conditions in order to face evolving energy demands.

Pharmacogenetic Study on Risperidone Long-Acting Injection: Influence of Cytochrome P450 2D6 and Pregnane X Receptor on Risperidone Exposure and Drug-Induced Side-Effects.

Risperidone is metabolized by polymorphic enzymes, and a large variability in plasma concentration and therapeutic response is observed. Risperidone long-acting injection (RLAI) avoids the first-pass effect, and little is known about the influence of gene polymorphisms involved in its pharmacokinetics. The influence on plasma concentrations of risperidone (RIS), its metabolite 9-hydroxy-risperidone, and on adverse effects were investigated for polymorphisms of cytochrome P450 2D6 (CYP2D6) (*3, *4, *5, *6), CYP3A (CYP3A4*1B, CYP3A4 rs4646437, CYP3A5*3, CYP3A7*1C), ABCB1 (1236C>T, 2677G>T, 3435C>T), NR1/2 coding for pregnane X receptor (rs1523130, rs2472677, rs7643645), and for CYP3A activity measured by a phenotyping test. Forty-two patients with at least 4 consecutive unchanged doses of RLAI were included in a multicenter cross-sectional study. A 55% lower dose-adjusted plasma levels of RIS were observed for CYP2D6 ultrarapid metabolizers (n = 5) as compared with CYP2D6 intermediate metabolizers (P < 0.007). NR1/2 polymorphism (rs7643645A>G) influenced RIS exposure with a 2.8-fold lower active moiety (P = 0.031) in GG compared with the AA genotype. This was confirmed in a second independent cohort (n = 16). Furthermore, high-density lipoprotein cholesterol was positively correlated with CYP3A activity (P = 0.01), and the NR1/2 (rs2472677) polymorphism was associated with different adverse effects including prolactin plasma levels adjusted for age and sex. In conclusion, our results confirmed the influence of CYP2D6 genotype on plasma levels of RIS. This is the first report on the influence of NR1/2 polymorphisms on RLAI exposure and on drug-induced adverse effects. These results should be validated in larger cohorts.

Planar homography: accuracy analysis and applications

Projective homography sits at the heart of many problems in image registration. In addition to many methods for estimating the homography parameters (R.I. Hartley and A. Zisserman, 2000), analytical expressions to assess the accuracy of the transformation parameters have been proposed (A. Criminisi et al., 1999). We show that these expressions provide less accurate bounds than those based on the earlier results of Weng et al. (1989). The discrepancy becomes more critical in applications involving the integration of frame-to-frame homographies and their uncertainties, as in the reconstruction of terrain mosaics and the camera trajectory from flyover imagery. We demonstrate these issues through selected examples

Assessment of changes in body water by bioimpedance in acutely ill surgical patients.

OBJECTIVE: To evaluate the relationship between changes in body bioelectrical impedance (BI) at 0.5, 50 and kHz and the changes in body weight, as an index of total body water changes, in acutely ill surgical patients during the rapid infusion of isotonic saline solution. DESIGN: Prospective clinical study. SETTING: Multidisciplinary surgical ICU in a university hospital. PATIENTS: Twelve male patients treated for acute surgical illness (multiple trauma n = 5, major surgery n = 7). Selection criteria: stable cardiovascular parameters, normal cardiac function, signs of hypovolemia (CVP < or = 5 mmHg, urine output < 1 ml/kg x h). INTERVENTIONS: After baseline measurements, a 60 min fluid challenge test was performed with normal saline solution, 0.25 ml/kg/min [corrected]. MEASUREMENTS AND RESULTS: Body weight (platform digital scale), total body impedance (four-surface electrode technique; measurements at 0.5, 50 and 100 kHz) and urine output. Fluid retention induced a progressive decrease in BI at 0.5, 50 and 100 kHz, but the changes were significant for BI 0.5 and BI 100 only, from 40 min after the beginning of the fluid therapy onwards. There was a significant negative correlation between changes in water retention and BI 0.5, with individual correlation coefficients ranging from -0.72 to 0.95 (p < 0.01-0.0001). The slopes of the regression lines indicated that for each kg of water change, there was a mean decrease in BI of 18 ohm, but a substantial inter-individual variability was noted. CONCLUSION: BI measured at low frequency can represent a valuable index of acute changes in body water in a group of surgical patients but not in a given individual.

Plasma membrane H+-ATPase regulation is required for auxin gradient formation preceding phototropic growth.

Phototropism is a growth response allowing plants to align their photosynthetic organs toward incoming light and thereby to optimize photosynthetic activity. Formation of a lateral gradient of the phytohormone auxin is a key step to trigger asymmetric growth of the shoot leading to phototropic reorientation. To identify important regulators of auxin gradient formation, we developed an auxin flux model that enabled us to test in silico the impact of different morphological and biophysical parameters on gradient formation, including the contribution of the extracellular space (cell wall) or apoplast. Our model indicates that cell size, cell distributions, and apoplast thickness are all important factors affecting gradient formation. Among all tested variables, regulation of apoplastic pH was the most important to enable the formation of a lateral auxin gradient. To test this prediction, we interfered with the activity of plasma membrane H(+)-ATPases that are required to control apoplastic pH. Our results show that H(+)-ATPases are indeed important for the establishment of a lateral auxin gradient and phototropism. Moreover, we show that during phototropism, H(+)-ATPase activity is regulated by the phototropin photoreceptors, providing a mechanism by which light influences apoplastic pH.

Molecular modeling of peptide-MHC class I complexes : applications to peptide based immunotherapy

Summary The specific CD8+ T cell immune response against tumors relies on the recognition by the T cell receptor (TCR) on cytotoxic T lymphocytes (CTL) of antigenic peptides bound to the class I major histocompatibility complex (MHC) molecule. Such tumor associated antigenic peptides are the focus of tumor immunotherapy with peptide vaccines. The strategy for obtaining an improved immune response often involves the design of modified tumor associated antigenic peptides. Such modifications aim at creating higher affinity and/or degradation resistant peptides and require precise structures of the peptide-MHC class I complex. In addition, the modified peptide must be cross-recognized by CTLs specific for the parental peptide, i.e. preserve the structure of the epitope. Detailed structural information on the modified peptide in complex with MHC is necessary for such predictions. In this thesis, the main focus is the development of theoretical in silico methods for prediction of both structure and cross-reactivity of peptide-MHC class I complexes. Applications of these methods in the context of immunotherapy are also presented. First, a theoretical method for structure prediction of peptide-MHC class I complexes is developed and validated. The approach is based on a molecular dynamics protocol to sample the conformational space of the peptide in its MHC environment. The sampled conformers are evaluated using conformational free energy calculations. The method, which is evaluated for its ability to reproduce 41 X-ray crystallographic structures of different peptide-MHC class I complexes, shows an overall prediction success of 83%. Importantly, in the clinically highly relevant subset of peptide-HLAA*0201 complexes, the prediction success is 100%. Based on these structure predictions, a theoretical approach for prediction of cross-reactivity is developed and validated. This method involves the generation of quantitative structure-activity relationships using three-dimensional molecular descriptors and a genetic neural network. The generated relationships are highly predictive as proved by high cross-validated correlation coefficients (0.78-0.79). Together, the here developed theoretical methods open the door for efficient rational design of improved peptides to be used in immunotherapy. Résumé La réponse immunitaire spécifique contre des tumeurs dépend de la reconnaissance par les récepteurs des cellules T CD8+ de peptides antigéniques présentés par les complexes majeurs d'histocompatibilité (CMH) de classe I. Ces peptides sont utilisés comme cible dans l'immunothérapie par vaccins peptidiques. Afin d'augmenter la réponse immunitaire, les peptides sont modifiés de façon à améliorer l'affinité et/ou la résistance à la dégradation. Ceci nécessite de connaître la structure tridimensionnelle des complexes peptide-CMH. De plus, les peptides modifiés doivent être reconnus par des cellules T spécifiques du peptide natif. La structure de l'épitope doit donc être préservée et des structures détaillées des complexes peptide-CMH sont nécessaires. Dans cette thèse, le thème central est le développement des méthodes computationnelles de prédiction des structures des complexes peptide-CMH classe I et de la reconnaissance croisée. Des applications de ces méthodes de prédiction à l'immunothérapie sont également présentées. Premièrement, une méthode théorique de prédiction des structures des complexes peptide-CMH classe I est développée et validée. Cette méthode est basée sur un échantillonnage de l'espace conformationnel du peptide dans le contexte du récepteur CMH classe I par dynamique moléculaire. Les conformations sont évaluées par leurs énergies libres conformationnelles. La méthode est validée par sa capacité à reproduire 41 structures des complexes peptide-CMH classe I obtenues par cristallographie aux rayons X. Le succès prédictif général est de 83%. Pour le sous-groupe HLA-A*0201 de complexes de grande importance pour l'immunothérapie, ce succès est de 100%. Deuxièmement, à partir de ces structures prédites in silico, une méthode théorique de prédiction de la reconnaissance croisée est développée et validée. Celle-ci consiste à générer des relations structure-activité quantitatives en utilisant des descripteurs moléculaires tridimensionnels et un réseau de neurones couplé à un algorithme génétique. Les relations générées montrent une capacité de prédiction remarquable avec des valeurs de coefficients de corrélation de validation croisée élevées (0.78-0.79). Les méthodes théoriques développées dans le cadre de cette thèse ouvrent la voie du design de vaccins peptidiques améliorés.

Measuring conduct and cost parameters in the Spanish air transport market

This paper estimates a model of airline competition for the Spanish air transport market. I test the explanatory power of alternative oligopoly models with capacity constraints. In addition, I analyse the degree of density economies. Results show that Spanish airlines conduct follows a price-leadership scheme so that it is less competitive than the Cournot solution. I also find evidence that thin routes can be considered as natural monopolies

Modeling water movement in horizontal columns using fractal theory

Fractal mathematics has been used to characterize water and solute transport in porous media and also to characterize and simulate porous media properties. The objective of this study was to evaluate the correlation between the soil infiltration parameters sorptivity (S) and time exponent (n) and the parameters dimension (D) and the Hurst exponent (H). For this purpose, ten horizontal columns with pure (either clay or loam) and heterogeneous porous media (clay and loam distributed in layers in the column) were simulated following the distribution of a deterministic Cantor Bar with fractal dimension H" 0.63. Horizontal water infiltration experiments were then simulated using Hydrus 2D software. The sorptivity (S) and time exponent (n) parameters of the Philip equation were estimated for each simulation, using the nonlinear regression procedure of the statistical software package SAS®. Sorptivity increased in the columns with the loam content, which was attributed to the relation of S with the capillary radius. The time exponent estimated by nonlinear regression was found to be less than the traditional value of 0.5. The fractal dimension estimated from the Hurst exponent was 17.5 % lower than the fractal dimension of the Cantor Bar used to generate the columns.

Modeling changes in organic carbon stocks for distinct soils in southeastern brazil after four eucalyptus rotations using the century model

Soil organic matter (SOM) plays an important role in carbon (C) cycle and soil quality. Considering the complexity of factors that control SOM cycling and the long time it usually takes to observe changes in SOM stocks, modeling constitutes a very important tool to understand SOM cycling in forest soils. The following hypotheses were tested: (i) soil organic carbon (SOC) stocks would be higher after several rotations of eucalyptus than in low-productivity pastures; (ii) SOC values simulated by the Century model would describe the data better than the mean of observations. So, the aims of the current study were: (i) to evaluate the SOM dynamics using the Century model to simulate the changes of C stocks for two eucalyptus chronosequences in the Rio Doce Valley, Minas Gerais State, Brazil; and (ii) to compare the C stocks simulated by Century with the C stocks measured in soils of different Orders and regions of the Rio Doce Valley growing eucalyptus. In Belo Oriente (BO), short-rotation eucalyptus plantations had been cultivated for 4.0; 13.0, 22.0, 32.0 and 34.0 years, at a lower elevation and in a warmer climate, while in Virginópolis (VG), these time periods were 8.0, 19.0 and 33.0 years, at a higher elevation and in a milder climate. Soil samples were collected from the 0-20 cm layer to estimate C stocks. Results indicate that the C stocks simulated by the Century model decreased after 37 years of poorly managed pastures in areas previously covered by native forest in the regions of BO and VG. The substitution of poorly managed pastures by eucalyptus in the early 1970´s led to an average increase of C of 0.28 and 0.42 t ha-1 year-1 in BO and VG, respectively. The measured C stocks under eucalyptus in distinct soil Orders and independent regions with variable edapho-climate conditions were not far from the values estimated by the Century model (root mean square error - RMSE = 20.9; model efficiency - EF = 0.29) despite the opposite result obtained with the statistical procedure to test the identity of analytical methods. Only for lower soil C stocks, the model over-estimated the C stock in the 0-20 cm layer. Thus, the Century model is highly promising to detect changes in C stocks in distinct soil orders under eucalyptus, as well as to indicate the impact of harvest residue management on SOM in future rotations.

Determination of hydraulic properties of a tropical soil of Hawaii using column experiments and inverse modeling

A method for determining soil hydraulic properties of a weathered tropical soil (Oxisol) using a medium-sized column with undisturbed soil is presented. The method was used to determine fitting parameters of the water retention curve and hydraulic conductivity functions of a soil column in support of a pesticide leaching study. The soil column was extracted from a continuously-used research plot in Central Oahu (Hawaii, USA) and its internal structure was examined by computed tomography. The experiment was based on tension infiltration into the soil column with free outflow at the lower end. Water flow through the soil core was mathematically modeled using a computer code that numerically solves the one-dimensional Richards equation. Measured soil hydraulic parameters were used for direct simulation, and the retention and soil hydraulic parameters were estimated by inverse modeling. The inverse modeling produced very good agreement between model outputs and measured flux and pressure head data for the relatively homogeneous column. The moisture content at a given pressure from the retention curve measured directly in small soil samples was lower than that obtained through parameter optimization based on experiments using a medium-sized undisturbed soil column.