A compositional approach to stable isotope data analysis

Isotopic data are currently becoming an important source of information regardingsources, evolution and mixing processes of water in hydrogeologic systems. However, itis not clear how to treat with statistics the geochemical data and the isotopic datatogether. We propose to introduce the isotopic information as new parts, and applycompositional data analysis with the resulting increased composition. Results areequivalent to downscale the classical isotopic delta variables, because they are alreadyrelative (as needed in the compositional framework) and isotopic variations are almostalways very small. This methodology is illustrated and tested with the study of theLlobregat River Basin (Barcelona, NE Spain), where it is shown that, though verysmall, isotopic variations comp lement geochemical principal components, and help inthe better identification of pollution sources

CODAPACK3D. A new version of Compositional Data Package

The statistical analysis of compositional data should be treated using logratios of parts,which are difficult to use correctly in standard statistical packages. For this reason afreeware package, named CoDaPack was created. This software implements most of thebasic statistical methods suitable for compositional data.In this paper we describe the new version of the package that now is calledCoDaPack3D. It is developed in Visual Basic for applications (associated with Excel©),Visual Basic and Open GL, and it is oriented towards users with a minimum knowledgeof computers with the aim at being simple and easy to use.This new version includes new graphical output in 2D and 3D. These outputs could bezoomed and, in 3D, rotated. Also a customization menu is included and outputs couldbe saved in jpeg format. Also this new version includes an interactive help and alldialog windows have been improved in order to facilitate its use.To use CoDaPack one has to access Excel© and introduce the data in a standardspreadsheet. These should be organized as a matrix where Excel© rows correspond tothe observations and columns to the parts. The user executes macros that returnnumerical or graphical results. There are two kinds of numerical results: new variablesand descriptive statistics, and both appear on the same sheet. Graphical output appearsin independent windows. In the present version there are 8 menus, with a total of 38submenus which, after some dialogue, directly call the corresponding macro. Thedialogues ask the user to input variables and further parameters needed, as well aswhere to put these results. The web site http://ima.udg.es/CoDaPack contains thisfreeware package and only Microsoft Excel© under Microsoft Windows© is required torun the software.Kew words: Compositional data Analysis, Software

Multivariate ARIMA Compositional Time Series Analysis

A compositional time series is obtained when a compositional data vector is observed atdifferent points in time. Inherently, then, a compositional time series is a multivariatetime series with important constraints on the variables observed at any instance in time.Although this type of data frequently occurs in situations of real practical interest, atrawl through the statistical literature reveals that research in the field is very much in itsinfancy and that many theoretical and empirical issues still remain to be addressed. Anyappropriate statistical methodology for the analysis of compositional time series musttake into account the constraints which are not allowed for by the usual statisticaltechniques available for analysing multivariate time series. One general approach toanalyzing compositional time series consists in the application of an initial transform tobreak the positive and unit sum constraints, followed by the analysis of the transformedtime series using multivariate ARIMA models. In this paper we discuss the use of theadditive log-ratio, centred log-ratio and isometric log-ratio transforms. We also presentresults from an empirical study designed to explore how the selection of the initialtransform affects subsequent multivariate ARIMA modelling as well as the quality ofthe forecasts

Compositional analysis of bivariate discrete probabilities

A joint distribution of two discrete random variables with finite support can be displayed as a two way table of probabilities adding to one. Assume that this table hasn rows and m columns and all probabilities are non-null. This kind of table can beseen as an element in the simplex of n · m parts. In this context, the marginals areidentified as compositional amalgams, conditionals (rows or columns) as subcompositions. Also, simplicial perturbation appears as Bayes theorem. However, the Euclideanelements of the Aitchison geometry of the simplex can also be translated into the tableof probabilities: subspaces, orthogonal projections, distances.Two important questions are addressed: a) given a table of probabilities, which isthe nearest independent table to the initial one? b) which is the largest orthogonalprojection of a row onto a column? or, equivalently, which is the information in arow explained by a column, thus explaining the interaction? To answer these questionsthree orthogonal decompositions are presented: (1) by columns and a row-wise geometric marginal, (2) by rows and a columnwise geometric marginal, (3) by independenttwo-way tables and fully dependent tables representing row-column interaction. Animportant result is that the nearest independent table is the product of the two (rowand column)-wise geometric marginal tables. A corollary is that, in an independenttable, the geometric marginals conform with the traditional (arithmetic) marginals.These decompositions can be compared with standard log-linear models.Key words: balance, compositional data, simplex, Aitchison geometry, composition,orthonormal basis, arithmetic and geometric marginals, amalgam, dependence measure,contingency table

Compositional data and Simpson's paradox

Simpson's paradox, also known as amalgamation or aggregation paradox, appears whendealing with proportions. Proportions are by construction parts of a whole, which canbe interpreted as compositions assuming they only carry relative information. TheAitchison inner product space structure of the simplex, the sample space of compositions, explains the appearance of the paradox, given that amalgamation is a nonlinearoperation within that structure. Here we propose to use balances, which are specificelements of this structure, to analyse situations where the paradox might appear. Withthe proposed approach we obtain that the centre of the tables analysed is a naturalway to compare them, which avoids by construction the possibility of a paradox.Key words: Aitchison geometry, geometric mean, orthogonal projection

Compositional Time Series: An Application

The composition of the labour force is an important economic factor for a country.Often the changes in proportions of different groups are of interest.I this paper we study a monthly compositional time series from the Swedish LabourForce Survey from 1994 to 2005. Three models are studied: the ILR-transformed series,the ILR-transformation of the compositional differenced series of order 1, and the ILRtransformationof the compositional differenced series of order 12. For each of thethree models a VAR-model is fitted based on the data 1994-2003. We predict the timeseries 15 steps ahead and calculate 95 % prediction regions. The predictions of thethree models are compared with actual values using MAD and MSE and the predictionregions are compared graphically in a ternary time series plot.We conclude that the first, and simplest, model possesses the best predictive power ofthe three models

Compositional evolution with mass transfer in closed systems

Evolution of compositions in time, space, temperature or other covariates is frequentin practice. For instance, the radioactive decomposition of a sample changes its composition with time. Some of the involved isotopes decompose into other isotopes of thesample, thus producing a transfer of mass from some components to other ones, butpreserving the total mass present in the system. This evolution is traditionally modelledas a system of ordinary di erential equations of the mass of each component. However,this kind of evolution can be decomposed into a compositional change, expressed interms of simplicial derivatives, and a mass evolution (constant in this example). A rst result is that the simplicial system of di erential equations is non-linear, despiteof some subcompositions behaving linearly.The goal is to study the characteristics of such simplicial systems of di erential equa-tions such as linearity and stability. This is performed extracting the compositional differential equations from the mass equations. Then, simplicial derivatives are expressedin coordinates of the simplex, thus reducing the problem to the standard theory ofsystems of di erential equations, including stability. The characterisation of stabilityof these non-linear systems relays on the linearisation of the system of di erential equations at the stationary point, if any. The eigenvelues of the linearised matrix and theassociated behaviour of the orbits are the main tools. For a three component system,these orbits can be plotted both in coordinates of the simplex or in a ternary diagram.A characterisation of processes with transfer of mass in closed systems in terms of stability is thus concluded. Two examples are presented for illustration, one of them is aradioactive decay

Modulation of expression and activity of cytochrome P450s and alteration of praziquantel kinetics during murine schistosomiasis

In this study, we investigated the expression and activity of liver cytochrome P450s (CYPs) and praziquantel (PZQ) kinetics in mice infected with Schistosoma mansoni. Swiss Webster (SW) mice of both genders were infected (100 cercariae) on postnatal day 10 and killed on post-infection days (PIDs) 30 or 55. Non-infected mice of the same age and sex served as controls. Regardless of mouse sex, infection depressed the activities of CYP1A [ethoxy/methoxy-resorufin-O-dealkylases (EROD/MROD)], 2B9/10 [pentoxy/benzyloxy-resorufin-O-dealkylases (PROD, BROD)], 2E1 [p-nitrophenol-hydroxylase (PNPH)] and 3A11 [erythromycin N-demethylase (END)] on PID 55 but not on PID 30. On PID 55, infection decreased liver CYP mRNA levels (real-time reverse transcription-polymerase chain reaction). On PID 30, whereas mRNA levels remained unaltered in males, they were depressed in females. Plasma PZQ (200 and 400 mg/kg body weight intraperitoneally) levels were measured (high-performance liquid chromatography) at different post-treatment intervals. In males and females, infection delayed the PZQ clearance on PID 55, but not on PID 30. Therefore, it can be concluded that schistosomiasis down-modulated CYP expression and activity and delayed PZQ clearance on PID 55, when a great number of parasite eggs were lodged in the liver. On PID 30, when egg-laying was initiated by the worms, no change of CYP expression and activity was found, except for a depression of CYP1A2 and 3A11 mRNAs in female mice.

Modulation of phosphorylation of neuronal cytoskeletal proteins by neuronal depolarization.

1. The neuronal cytoskeletal protein tau and the carboxy tails of cytoskeletal proteins neurofilament-M (NF-M) and neurofilament-H (NF-H) are phosphorylated on serine residues by the cyclin-dependent kinase cdk-5. 2. In aggregating neuronal-glial cultures we show that veratridine-mediated cation influx causes dephosphorylation of tau, NF-M and NF-H. Dephosphorylation was blocked specifically by cyclosporine A but not by okadiac acid at concentrations up to 200 nM. 3. These results suggest that veratridine-triggered cation influx causes activation of PP-2B (calcineurin) leading to dephosphorylation of these cytoskeletal proteins.

Modulation of MDR1 and MRP3 Gene Expression in Lung Cancer Cells after Paclitaxel and Carboplatin Exposure

Carboplatin-paclitaxel is a reference regimen in the treatment of locally advanced or disseminated non-small cell lung cancer (NSCLC). This paper discusses the multidrug resistance developed with this drug combination, which is one of the major obstacles to successful treatment. In order to understand and overcome the drug resistance pattern of NSCLC after carboplatin plus paclitaxel exposure, levels of mRNA expression of multidrug resistance 1 (MDR1) and multidrug resistance-associated protein 3 (MRP3) were investigated in primary NSCLC cell lines (A-549 and A-427) and a metastasis-derived NSCLC cell line (NODO). Our results showed that exposure of the three NSCLC lines to plasma concentrations of paclitaxel (5 μM) produced an increase in MDR1 expression, while MRP3 showed no alteration in expression. By contrast, the same cells exposed to carboplatin plasma concentrations (30 μM) showed overexpression of MRP3. In these cells, MDR1 showed no expression changes. Interestingly, the combination of both paclitaxel and carboplatin caused increased expression of the MDR1 drug resistance gene rather than the individual treatments. These results suggest that carboplatin and paclitaxel may induce drug resistance mediated by MDR1 and MRP3, which may be enhanced by the simultaneous use of both drugs.

MicroRNA-200 family modulation in distinct breast cancer phenotypes.

The epithelial to mesenchymal transition (EMT) contributes to tumor invasion and metastasis in a variety of cancer types. In human breast cancer, gene expression studies have determined that basal-B/claudin-low and metaplastic cancers exhibit EMT-related characteristics, but the molecular mechanisms underlying this observation are unknown. As the family of miR-200 microRNAs has been shown to regulate EMT in normal tissues and cancer, here we evaluated whether the expression of the miR-200 family (miR-200f) and their epigenetic state correlate with EMT features in human breast carcinomas. We analyzed by qRT-PCR the expression of miR-200f members and various EMT-transcriptional inducers in a series of 70 breast cancers comprising an array of phenotypic subtypes: estrogen receptor positive (ER+), HER2 positive (HER2+), and triple negative (TN), including a subset of metaplastic breast carcinomas (MBCs) with sarcomatous (homologous or heterologous) differentiation. No MBCs with squamous differentiation were included. The DNA methylation status of miR-200f loci in tumor samples were inspected using Sequenom MassArray® MALDI-TOF platform. We also used two non-tumorigenic breast basal cell lines that spontaneously undergo EMT to study the modulation of miR-200f expression during EMT in vitro. We demonstrate that miR-200f is strongly decreased in MBCs compared with other cancer types. TN and HER2+ breast cancers also exhibited lower miR-200f expression than ER+ tumors. Significantly, the decreased miR-200f expression found in MBCs is accompanied by an increase in the expression levels of EMT-transcriptional inducers, and hypermethylation of the miR-200c-141 locus. Similar to tumor samples, we demonstrated that downregulation of miR-200f and hypermethylation of the miR-200c-141 locus, together with upregulation of EMT-transcriptional inducers also occur in an in vitro cellular model of spontaneous EMT. Thus, the expression and methylation status of miR-200f could be used as hypothetical biomarkers to assess the occurrence of EMT in breast cancer.

Gut microbiota in children with type 1 diabetes differs from that in healthy children: a case-control study.

BACKGROUND A recent study using a rat model found significant differences at the time of diabetes onset in the bacterial communities responsible for type 1 diabetes modulation. We hypothesized that type 1 diabetes in humans could also be linked to a specific gut microbiota. Our aim was to quantify and evaluate the difference in the composition of gut microbiota between children with type 1 diabetes and healthy children and to determine the possible relationship of the gut microbiota of children with type 1 diabetes with the glycemic level. METHODS A case-control study was carried out with 16 children with type 1 diabetes and 16 healthy children. The fecal bacteria composition was investigated by polymerase chain reaction-denaturing gradient gel electrophoresis and real-time quantitative polymerase chain reaction. RESULTS The mean similarity index was 47.39% for the healthy children and 37.56% for the children with diabetes, whereas the intergroup similarity index was 26.69%. In the children with diabetes, the bacterial number of Actinobacteria and Firmicutes, and the Firmicutes to Bacteroidetes ratio were all significantly decreased, with the quantity of Bacteroidetes significantly increased with respect to healthy children. At the genus level, we found a significant increase in the number of Clostridium, Bacteroides and Veillonella and a significant decrease in the number of Lactobacillus, Bifidobacterium, Blautia coccoides/Eubacterium rectale group and Prevotella in the children with diabetes. We also found that the number of Bifidobacterium and Lactobacillus, and the Firmicutes to Bacteroidetes ratio correlated negatively and significantly with the plasma glucose level while the quantity of Clostridium correlated positively and significantly with the plasma glucose level in the diabetes group. CONCLUSIONS This is the first study showing that type 1 diabetes is associated with compositional changes in gut microbiota. The significant differences in the number of Bifidobacterium, Lactobacillus and Clostridium and in the Firmicutes to Bacteroidetes ratio observed between the two groups could be related to the glycemic level in the group with diabetes. Moreover, the quantity of bacteria essential to maintain gut integrity was significantly lower in the children with diabetes than the healthy children. These findings could be useful for developing strategies to control the development of type 1 diabetes by modifying the gut microbiota.

Modulation of the Endocannabinoids N-Arachidonoylethanolamine (AEA) and 2-Arachidonoylglycerol (2-AG) on Executive Functions in Humans.

Animal studies point to an implication of the endocannabinoid system on executive functions. In humans, several studies have suggested an association between acute or chronic use of exogenous cannabinoids (Δ9-tetrahydrocannabinol) and executive impairments. However, to date, no published reports establish the relationship between endocannabinoids, as biomarkers of the cannabinoid neurotransmission system, and executive functioning in humans. The aim of the present study was to explore the association between circulating levels of plasma endocannabinoids N-arachidonoylethanolamine (AEA) and 2-Arachidonoylglycerol (2-AG) and executive functions (decision making, response inhibition and cognitive flexibility) in healthy subjects. One hundred and fifty seven subjects were included and assessed with the Wisconsin Card Sorting Test; Stroop Color and Word Test; and Iowa Gambling Task. All participants were female, aged between 18 and 60 years and spoke Spanish as their first language. Results showed a negative correlation between 2-AG and cognitive flexibility performance (r = -.37; p<.05). A positive correlation was found between AEA concentrations and both cognitive flexibility (r = .59; p<.05) and decision making performance (r = .23; P<.05). There was no significant correlation between either 2-AG (r = -.17) or AEA (r = -.08) concentrations and inhibition response. These results show, in humans, a relevant modulation of the endocannabinoid system on prefrontal-dependent cognitive functioning. The present study might have significant implications for the underlying executive alterations described in some psychiatric disorders currently associated with endocannabinoids deregulation (namely drug abuse/dependence, depression, obesity and eating disorders). Understanding the neurobiology of their dysexecutive profile might certainly contribute to the development of new treatments and pharmacological approaches.

Acid-sensing ion channels : modulation by serine-proteases and involvement in acid-induced action potential generation in hippocampal neurons

SUMMARY Acid-sensing ion channels (ASICs) are non-voltage gated sodium channels. They are activated by rapid extracellular acidification and generate an inactivating inward current. Four ASIC genes have been cloned: ASIC1, 2, 3 and 4, with variants a and b for ASIC1and AS1C2. ASICs are expressed in neurons of the central (CNS) and peripheral nervous system (PNS). In the CNS, ASICs have a role in learning, memory, as well as in neuronal death in ischemia. In the PNS, ASICs are involved in the perception of acid-induced pain, as well as in mechanoperception. In one part of my thesis project, we addressed the question of the mechanism of regulation of ASIC1 a by the serine protease trypsin at the molecular level. Trypsin modifies the function of ASIC1 a but not of ASIC1b. In order to identify the channel region responsible for this effect, we created chimeras between ASIC1 a and 1b. Subsequently, to identify the exact trypsin target(s), we mutated predicted trypsin sites in the region identified by the chimera. In the second part of a project, we investigated the role of ASICs at the cellular level, in neuronal signaling. Using the whole-cell patch clamp in hippocampal neuronal culture, we studied the potential involvement of ASICs in action potential (AP) generation. In the first part of the thesis work, we showed that trypsin modifies ASIC1a function: it shifts the pH activation and the steady-state inactivation curve towards more acidic values and accelerates the time course of the channel recovery from inactivation. We also showed that trypsin cleaves ASIC1a and that the functional effect and a channel cleavage correlate. In the inactivated state, channels cannot be modified by trypsin. Cleavage occurs in a channel region that is also important for inactivation of all ASICs; a part of this region is critical for the inhibition of ASIC1 a by the spider toxin Psalmotoxin1. In the second part of the thesis work, we showed that ASIC activity can modulate AP generation. ASIC activity by itself can induce trains of APs. In situations in which this activity by itself is not sufficient to induce APs, it can contribute to AP generation. During high neuronal activity, ASIC activity can block already existing trains of APs. In conclusion, depending on the activity of neuron in a particular moment, ASICs can differently modulate AP generation; they can induce, facilitate or inhibit APs. We also showed that trypsin changes the capability of ASICs to modulate AP generation by shifting the pH dependence to more acidic values, which adapts channel gating to pH conditions which may occur in pathological conditions such as ischemia. Our finding that trypsin modifies ASIC1 a function identifies a novel pharmacological tool, and proposes a mechanism of ASIC1a regulation that may have a physiological importance. The identification of the exact site of trypsin action gives insight to the molecular mechanisms of ASIC regulation. This work proposes a role in modulation of AP generation for ASICs in the CNS. RESUME Les canaux ASIC sont les canaux ioniques activés par l'acidification rapide extracellulaire. Activés, ils génèrent un courant entrant qui inactive en présence de stimulus acide. Quatre gènes ASIC ont été clonés, ASIC1, 2, 3 et 4, avec les variants a et b pour ASIC1 et 2. Les ASICs sont exprimés dans les neurones du système nerveux central (SNC) et périphérique (SNP). Dans le SNC, les ASIC ont un rôle dans le mémoire, apprentissage et la mort neuronale dans t'ischémie. Dans le SNP, ils ont un rôle dans la perception de la douleur et méchanosensation. Dans une partie de mon projet de thèse, nous avons étudié les mécanismes de la régulation d'ASIC1a par la sérine-protéase trypsine au niveau moléculaire. La trypsine modifie la fonction d'ASIC1a et pas ASIC1b. Nous avons créé les chimères entre ASIC1 a et 1 b, afin d'identifier la région du canal responsable pour l'effet. Pour identifier le(s) site(s) exactes de l'action de la trypsine, nous avons muté les sites potentiels de la trypsine dans la région identifiée par les chimères. Dans la deuxième partie du projet, nous avons étudié le rôle des ASICs au niveau cellulaire. En utilisant la technique du patch clamp dans les cultures des neurones de l'hippocampe, nous avons étudié l'implication des ASICs dans la génération des potentiels d'action (PA). Nous avons montré que la trypsine agit sur le canal ASIC1a ; elle décale l'activation et « steady-state » inactivation vers les valeurs plus acides, et elle raccourcit le temps du « recovery » du canal. La trypsine coupe ASIC1a sur le résidu K145 et l'effet fonctionnel et la coupure corrèlent. Nous avons identifié la région du canal responsable pour l'inactivation de tous les ASICs ; une partie de cette région est responsable pour ['inhibition d'ASIC1 a par la Psalmotoxinel . Nous avons montré que les ASICs peuvent moduler la génération des PAs. L'activité des ASICs peut induire les trains des PAs. Quand l'activité des ASICs n'est pas suffisante pour induire le PA, elle peut contribuer à sa génération. Pendant l'activité neuronale forte, l'activité des ASICs peut bloquer les trains des PAs qui existent déjà. En conclusion, dépendant de l'activité neuronale, les ASICs peuvent moduler la génération des PAs différemment ; ils peuvent induire, faciliter ou inhiber les PAs. La trypsine change la capacité des ASICs de moduler les PAs. Après l'action de la trypsine, les ASICs peuvent moduler la génération des PAs dans les conditions légèrement acides, suivies par les fluctuations du pH acide, qui peuvent exister dans l'ischémie. Le fait que la trypsine agit sur ASIC1a définit l'outil pharmacologique et propose le mécanisme de la régulation d'ASICI a qui pourrait avoir l'importance physiologique. L'identification du site de l'action de la trypsine éclaircit les mécanismes moléculaires de la régulation des ASICs. Cette étude propose un rôle des ASICs dans la modulation de la génération des PAs. Résumé pour le public large Les neurones sont les cellules de système nerveux dont la fonction est la signalisation. Comme toutes les autres cellules, les neurones ont une membrane qui sépare l'intérieur du milieu extérieur. Cette membrane est imperméable pour des particules chargées (ions). Dans cette membrane existent les protéines spécifiques, « canaux », qui permettent le transport des ions d'un côté de la membrane à l'autre, comme réponse aux stimuli différents. Ce transport des ions à travers la membrane génère un courant, qu'on peut mesurer. Ce courant est la base de la communication entre les neurones, ou, ce qu'on appelle la signalisation neuronale. Quand ce courant est suffisamment grand, il permet la génération du potentiel d'action, qui est le message principal de communication neuronale. Les canaux ASIC (acid-sensing ion channel), que nous étudions dans le laboratoire, sont activés par les acides. Les acides sont relâchés dans beaucoup de situations dans le système nerveux. Les ASIC ont été découverts récemment (en 1996), et nous ne connaissons pas encore très bien toutes les fonctions de ces canaux. Nous savons qu'ils ont un rôle dans le mémoire, apprentissage, la sensation de la douleur et l'infarctus cérébral. Dans la première partie de ce projet de thèse, nous avons voulu mieux comprendre comment fonctionnent ces canaux. Pour faire ça, nous avons étudié la régulation des ASICs par une protéine, trypsine, qui coupe le canal ASIC. Nous avons étudié ou exactement la trypsine coupe le canal et quels effets ça produit sur la fonction du canal. Dans la deuxième partie du projet de thèse, nous avons voulu mieux connaître comment le canal fonctionne au niveau de la cellule, comment il interagit avec les autres canaux et si il a un rôle dans la génération des potentiels d'action. Nous avons pu montrer que la trypsine change la fonction du canal, ce qui lui permet de fonctionner différemment. Nous avons aussi déterminé ou exactement ta trypsine coupe le canal. Au niveau de la cellule, nous avons montré que les ASIC peuvent moduler la génération des potentiels d'action, étant, dépendant de l'activité du neurone, soit activateurs, soit inhibiteurs. La trypsine est une molécule qui peut être libérée dans le système nerveux pendant certaines conditions, comme l'infarctus cérébral. A cause de ça, les connaissances que la trypsine agit sur le anal ASIC pourraient être important physiologiquement. La connaissance de l'endroit exacte ou la trypsine coupe le canal nous aide à mieux comprendre la relation structure-fonction du canal. La modulation de la génération des potentiels d'actions par les ASIC indique que ces canaux peuvent avoir un rôle important dans la signalisation neuronale.