Application of compositional data analysis to geochemical data of marine sediments

In an earlier investigation (Burger et al., 2000) five sediment cores near the RodriguesTriple Junction in the Indian Ocean were studied applying classical statistical methods(fuzzy c-means clustering, linear mixing model, principal component analysis) for theextraction of endmembers and evaluating the spatial and temporal variation ofgeochemical signals. Three main factors of sedimentation were expected by the marinegeologists: a volcano-genetic, a hydro-hydrothermal and an ultra-basic factor. Thedisplay of fuzzy membership values and/or factor scores versus depth providedconsistent results for two factors only; the ultra-basic component could not beidentified. The reason for this may be that only traditional statistical methods wereapplied, i.e. the untransformed components were used and the cosine-theta coefficient assimilarity measure.During the last decade considerable progress in compositional data analysis was madeand many case studies were published using new tools for exploratory analysis of thesedata. Therefore it makes sense to check if the application of suitable data transformations,reduction of the D-part simplex to two or three factors and visualinterpretation of the factor scores would lead to a revision of earlier results and toanswers to open questions . In this paper we follow the lines of a paper of R. Tolosana-Delgado et al. (2005) starting with a problem-oriented interpretation of the biplotscattergram, extracting compositional factors, ilr-transformation of the components andvisualization of the factor scores in a spatial context: The compositional factors will beplotted versus depth (time) of the core samples in order to facilitate the identification ofthe expected sources of the sedimentary process.Kew words: compositional data analysis, biplot, deep sea sediments

Dynamic graphics of parametrically linked multivariate methods used in compositional data analysis

Many multivariate methods that are apparently distinct can be linked by introducing oneor more parameters in their definition. Methods that can be linked in this way arecorrespondence analysis, unweighted or weighted logratio analysis (the latter alsoknown as "spectral mapping"), nonsymmetric correspondence analysis, principalcomponent analysis (with and without logarithmic transformation of the data) andmultidimensional scaling. In this presentation I will show how several of thesemethods, which are frequently used in compositional data analysis, may be linkedthrough parametrizations such as power transformations, linear transformations andconvex linear combinations. Since the methods of interest here all lead to visual mapsof data, a "movie" can be made where where the linking parameter is allowed to vary insmall steps: the results are recalculated "frame by frame" and one can see the smoothchange from one method to another. Several of these "movies" will be shown, giving adeeper insight into the similarities and differences between these methods

Spatial Analysis of Cell Composition Data

A version of Matheron’s discrete Gaussian model is applied to cell composition data.The examples are for map patterns of felsic metavolcanics in two different areas. Q-Qplots of the model for cell values representing proportion of 10 km x 10 km cell areaunderlain by this rock type are approximately linear, and the line of best fit can be usedto estimate the parameters of the model. It is also shown that felsic metavolcanics in theAbitibi area of the Canadian Shield can be modeled as a fractal

New insights on river water chemistry by using non-centred simplicial principal component analysis: a case study

The use of perturbation and power transformation operations permits the investigation of linear processes in the simplex as in a vectorial space. When the investigated geochemical processes can be constrained by the use of well-known starting point, the eigenvectors of the covariance matrix of a non-centred principalcomponent analysis allow to model compositional changes compared with a reference point.The results obtained for the chemistry of water collected in River Arno (central-northern Italy) have open new perspectives for considering relative changes of the analysed variables and to hypothesise the relative effect of different acting physical-chemical processes, thus posing the basis for a quantitative modelling

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

Clinical performance of cervical restorations: a meta-analysis.

OBJECTIVES: To carry out a meta-analysis in order to assess the influencing factors on retention loss and marginal discoloration of cervical restorations made of composites and glass ionomer (derivates). METHODS: The literature was searched for prospective clinical studies on cervical restorations with an observation period of at least 18 months. RESULTS: Fifty clinical studies involving 40 adhesive systems matched the inclusion criteria. On average, 10% of the cervical fillings were lost and 24% exhibited marginal discoloration after 3 years. The variability ranged from 0% to 50% for retention loss and from 0% to 74% for marginal discoloration. Hardly any secondary caries was detected. When linear mixed models with a study and experiment effect were used, the analysis revealed that the adhesive/restorative class had the most significant influence, with 2-step self-etching adhesive systems performing best and 1-step self-etching adhesive systems performing worst; 3-step etch-and-rinse systems, glass ionomers/resin-modified glass ionomers, 2-step etch-and-rinse systems and polyacid-modified resin composites were ranked in between. Restorations placed in teeth whose dentin/enamel had been prepared/roughened showed a statistically significant higher retention rate than those placed in teeth with unprepared dentin (p<0.05). Beveling of the enamel and the type of isolation used (rubberdam/cotton rolls) had no significant influence. SIGNIFICANCE: The clinical performance of cervical restorations is significantly influenced by the type of adhesive system used and/or the adhesive class to which the system belonged and whether the dentin/enamel is prepared or not. 2-Step self-etching- and 3-step etch&rinse systems shall be chosen over 1-step self-etching systems and glass ionomer derivates. The dentin (and enamel) surface shall be roughened before placement of the restoration.

Factors associated with use of community mental health services by schizophrenia patients using multilevel analysis

BACKGROUND Persons with schizophrenia and related disorders may be particularly sensitive to a number of determinants of service use, including those related with illness, socio-demographic characteristics and organizational factors. The objective of this study is to identify factors associated with outpatient contacts at community mental health services of patients with schizophrenia or related disorders. METHODS This cross-sectional study analyzed 1097 patients. The main outcome measure was the total number of outpatient consultations during one year. Independent variables were related to socio-demographic, clinical and use of service factors. Data were collected from clinical records. RESULTS The multilevel linear regression model explained 46.35% of the variance. Patients with significantly more contacts with ambulatory services were not working and were receiving welfare benefits (p = 0.02), had no formal education (p = 0.02), had a global level of severity of two or three (four being the most severe) (p < 0.001), with one or more inpatient admissions (p < 0.001), and in contact with both types of professional (nurses and psychiatrists) (p < 0.001). The patients with the fewest ambulatory contacts were those with diagnoses of persistent delusional disorders (p = 0.04) and those who were attended by four of the 13 psychiatrists (p < 0.001). CONCLUSIONS As expected, the variables that explained the use of community service could be viewed as proxies for severity of illness. The most surprising finding, however, was that a group of four psychiatrists was also independently associated with use of ambulatory services by patients with schizophrenia or related disorders. More research is needed to carefully examine how professional support networks interact to affect use of mental health.

Inflammation in disseminated lesions: an analysis of CD4+, CD20+, CD68+, CD31+ and vW+ cells in non-ulcerated lesions of disseminated leishmaniasis

Disseminated leishmaniasis (DL) differs from other clinical forms of the disease due to the presence of many non-ulcerated lesions (papules and nodules) in non-contiguous areas of the body. We describe the histopathology of DL non-ulcerated lesions and the presence of CD4-, CD20-, CD68-, CD31- and von Willebrand factor (vW)-positive cells in the inflamed area. We analysed eighteen biopsies from non-ulcerated lesions and quantified the inflamed areas and the expression of CD4, CD20, CD68, CD31 and vW using Image-Pro software (Media Cybernetics). Diffuse lymphoplasmacytic perivascular infiltrates were found in dermal skin. Inflammation was observed in 3-73% of the total biopsy area and showed a significant linear correlation with the number of vW+ vessels. The most common cells were CD68+ macrophages, CD20+ B-cells and CD4+ T-cells. A significant linear correlation between CD4+ and CD20+ cells and the size of the inflamed area was also found. Our findings show chronic inflammation in all DL non-ulcerated lesions predominantly formed by macrophages, plasmacytes and T and B-cells. As the inflamed area expanded, the number of granulomas and extent of the vascular framework increased. Thus, we demonstrate that vessels may have an important role in the clinical evolution of DL lesions.

Cross-validation of bioelectrical impedance analysis for the assessment of body composition in a representative sample of 6- to 13-year-old children.

BACKGROUND/OBJECTIVES: (1) To cross-validate tetra- (4-BIA) and octopolar (8-BIA) bioelectrical impedance analysis vs dual-energy X-ray absorptiometry (DXA) for the assessment of total and appendicular body composition and (2) to evaluate the accuracy of external 4-BIA algorithms for the prediction of total body composition, in a representative sample of Swiss children. SUBJECTS/METHODS: A representative sample of 333 Swiss children aged 6-13 years from the Kinder-Sportstudie (KISS) (ISRCTN15360785). Whole-body fat-free mass (FFM) and appendicular lean tissue mass were measured with DXA. Body resistance (R) was measured at 50 kHz with 4-BIA and segmental body resistance at 5, 50, 250 and 500 kHz with 8-BIA. The resistance index (RI) was calculated as height(2)/R. Selection of predictors (gender, age, weight, RI4 and RI8) for BIA algorithms was performed using bootstrapped stepwise linear regression on 1000 samples. We calculated 95% confidence intervals (CI) of regression coefficients and measures of model fit using bootstrap analysis. Limits of agreement were used as measures of interchangeability of BIA with DXA. RESULTS: 8-BIA was more accurate than 4-BIA for the assessment of FFM (root mean square error (RMSE)=0.90 (95% CI 0.82-0.98) vs 1.12 kg (1.01-1.24); limits of agreement 1.80 to -1.80 kg vs 2.24 to -2.24 kg). 8-BIA also gave accurate estimates of appendicular body composition, with RMSE < or = 0.10 kg for arms and < or = 0.24 kg for legs. All external 4-BIA algorithms performed poorly with substantial negative proportional bias (r> or = 0.48, P<0.001). CONCLUSIONS: In a representative sample of young Swiss children (1) 8-BIA was superior to 4-BIA for the prediction of FFM, (2) external 4-BIA algorithms gave biased predictions of FFM and (3) 8-BIA was an accurate predictor of segmental body composition.

Effectiveness of holistic neuropsychological rehabilitation for Spanish population with acquired brain injury measured using Rasch analysis.

INTRODUCTION The Rasch model is increasingly used in the field of rehabilitation because it improves the accuracy of measurements of patient status and their changes after therapy. OBJECTIVE To determine the long-term effectiveness of a holistic neuropsychological rehabilitation program for Spanish outpatients with acquired brain injury (ABI) using Rasch analysis. METHODS Eighteen patients (ten with long evolution - patients who started the program > 6 months after ABI- and eight with short evolution) and their relatives attended the program for 6 months. Patients' and relatives' answers to the European Brain Injury Questionnaire and the Frontal Systems Behavior Scale at 3 time points (pre-intervention. post-intervention and 12 month follow-up) were transformed into linear measures called logits. RESULTS The linear measures revealed significant improvements with large effects at the follow-up assessment on cognitive and executive functioning, social and emotional self-regulation, apathy and mood. At follow-up, the short evolution group achieved greater improvements in mood and cognitive functioning than the long evolution patients. CONCLUSIONS The program showed long-term effectiveness for most of the variables, and it was more effective for mood and cognitive functioning when patients were treated early. Relatives played a key role in the effectiveness of the rehabilitation program.

An immuno-electron microscopical analysis of transcribing multinucleosomal templates: what happens to the histones?

Immuno-electron microscopy was used to visualize the structure of reconstituted chromatin after in vitro transcription by purified T7 RNA polymerase. T7 RNA polymerase disrupts the nucleosomal structure in the transcribed region. This disruption is not influenced by the template, linear or supercoiled, and the presence or absence of nucleosomal positioning sequences in the transcribed region. In this study, we used monoclonal autoantibodies reacting with the nucleosome core particles and epitopes within several regions of the four different core histones. Some of the residues recognized by the autoantibodies are accessible on the surface of the nucleosomes and some are more internal and therefore less exposed at the surface. We show that the loss of the nucleosomal configuration during transcription is due to the loss of histone/DNA binding and that at least part of the histones are transferred to the nascent RNA chains. Consequently, after in vitro transcription by T7 RNA polymerase, the nucleosomal template does not conserve its original configuration, and no interaction of antigen/antibodies is observed anymore in the region that has been transcribed. Therefore, we conclude that in our in vitro transcription assay, nucleosomes are detached from the template, and not simply unfolded with histones remaining attached to the DNA.

Characterisation of "fou2", a constitutive wound-activated mutant and analysis of the proteome and leaf physiology in the region proximal to wounds in Arabidopsis

Résumé : Les jasmonates (JA), une famille d'hor1none végétale, jouent un rôle central dans la réponse à la blessure, et aux attaques d'insectes et de pathogènes. Les JA sont principalement dérivés d'un acide gras, l'acide linolénique. L'addition par une lipoxygénase d'une molécule d'oxygène à l'acide linolénique initie la synthèse de JA. Cependant les mécanismes régulant l'activation de la biosynthèse de JA ne sont pas encore connus. C'est pour cette raison que dans ce travail, nous avons caractérisé chez Arabidopsis thaliana (l'Arabette des Dames) un mutant fou2 dont l'activité lipoxygénase est plus élevée que celle d'une plante sauvage. Les niveaux de JA sont constitutivement plus élevés et l'activation de la synthèse de JA après blessure est fortement plus induite chez fou2 que chez le type sauvage. En outre, fou2 est plus résistant au pathogène Botrytis cinerea et à la chenille Spodoptera littoralis. Afin de comprendre quel mécanisme chez fou2 génére ce phénotype, nous avons cloné le gène responsable du phénotype de fou2. Le mutant fou2 porte une mutation dans le gène d'un canal à deux pores transportant probablement du potassium, du lumen de la vacuole végétale vers le compartiment cytosolique. L'analyse du protéome de fou2 a permis d'identifier une expression plus élevée de sept protéines régulées par les JA ou le stress. La découverte de l'implication d'un canal dans le phénotype de fou2 renforce l'hypothèse que les flux de cations pourraient être impliqués dans les étapes précoces de la synthèse des JA. Nous avons également étudié le protéome et la physiologie d'une feuille blessée, Pour évaluer les changements d'expression protéique en réponse à la blessure et contrôlés par les JA, nous avons quantifié l'expression de 5937 protéines chez une plante d'Arabidopsis sauvage et chez un mutant incapable de synthétiser des JA. Parmi ces 5937 protéines, nous avons identifié 99 protéines régulées par la blessure chez le type sauvage. Nous avons observé pour 65% des protéines dont l'expression protéique changeait après blessure une bonne corrélation entre la quantité de transcrits et de protéines. Plusieurs enzymes de la voie des chorismates impliquées dans la biosynthèse des acides aminés phénoliques étaient induites par les JA après blessure. Une quantification des acides aminés a montré que les niveaux d'acides aminés phénoliques augmentaient significativement après blessure. La blessure induisait aussi des changements dans l'expression de protéines impliquées dans la réponse au stress et particulièrement au stress oxydatif. Nous avons quantifié l'état réduit et oxydé du glutathion, un tripeptide qui, sous sa forme réduite, est l'antioxydant majeur des cellules. Nous avons trouvé une quantité significativement plus élevée de glutathion oxydé chez le type sauvage blessé que chez la plante aus blessée. Ce résultat suggère que la génération d'un stress oxydatif et la proportion relative de glutathions réduits et oxydés sont contrôlés par les JA après blessure. Abstract : Plants possess a family of potent fatty acid-derived wound-response and developmental regulators: the jasmonates. These compounds are derived from the tri?unsaturated fatty acid a-linolenic-acid (18:3). Addition of an oxygen molecule to 18:3 by 13-lipoxygenases (13-LOX) initiates JA biosynthesis. Actually components regulating the activation of JA biosynthesis are poorly defined. Therefore we characterized in Arabidopsis thaliana the fatty acid Qxygenation upregulated 2 (fou2) mutant, which was previously isolated in a screen for mutants with an enhanced 13-LOX activity. As a consequence of this increased 13-LOX activity, JA levels in fou2 are higher than in wild type (WT) and wounding strongly increased JA biosynthesis compared to WT. fou2 was more resistant to the fungus Botrytis cinerea and the generalist caterpillar Spodaptera littomlis, The fou2 mutant carries a missense mutation in the Two Pore Channel 1 gene (TPCJ), which encodes a vacuolar cation channel transporting probably K* into the cytosol. Patchclamp analysis of fou2 vacuolar membranes showed faster time-dependent conductivity and activation of the mutated channel at lower membrane potentials than wild-type. Proteomic analysis of fou2 leaves identified increased levels of seven biotic stress- and JA- inducible proteins. The discovery of the implication of a channel in the fou2 phenotype strenghtens the hypothesis that cation fluxes might be implicated in early steps of JA synthesis. We further concentrated on the proteome and leaf physiology in the region proximal to wounds in Arabidopsis using the WT and the aos JA-biosynthesis deficient mutant in order to find JA- induced proteins changes. We used two successive proteomic methods to assess protein changes in response to wounding Arabidopsis leaves, two dimensional electrophoresis (2DE) and linear trap quadrupole ion-trap mass spectrometry. In total 5937 proteins were quantified. We identified 99 wound-regulated proteins in the WT. Most these proteins were also wound-regulated at the transcript level showing a good correlation between transcript and protein abundance. We identified several wound-regulated enzymes involved in amino acid biosynthesis and confirmed this result by amino acid quantification. Proteins involved in stress reponses were upregulated, particularly in redox species regulation. We found a significantly higher quantity of oxidized glutathione in wounded WT relative to wounded aos leaves. This result suggests that levels of reduced glutathione are controlled by JA after wounding.

Trabecular bone score (TBS) the new parameter of 2D texture analysis for the evaluation of 3D bone micro architecture status

X-ray is a technology that is used for numerous applications in the medical field. The process of X-ray projection gives a 2-dimension (2D) grey-level texture from a 3- dimension (3D) object. Until now no clear demonstration or correlation has positioned the 2D texture analysis as a valid indirect evaluation of the 3D microarchitecture. TBS is a new texture parameter based on the measure of the experimental variogram. TBS evaluates the variation between 2D image grey-levels. The aim of this study was to evaluate existing correlations between 3D bone microarchitecture parameters - evaluated from μCT reconstructions - and the TBS value, calculated on 2D projected images. 30 dried human cadaveric vertebrae were acquired on a micro-scanner (eXplorer Locus, GE) at isotropic resolution of 93 μm. 3D vertebral body models were used. The following 3D microarchitecture parameters were used: Bone volume fraction (BV/TV), Trabecular thickness (TbTh), trabecular space (TbSp), trabecular number (TbN) and connectivity density (ConnD). 3D/2D projections has been done by taking into account the Beer-Lambert Law at X-ray energy of 50, 100, 150 KeV. TBS was assessed on 2D projected images. Correlations between TBS and the 3D microarchitecture parameters were evaluated using a linear regression analysis. Paired T-test is used to assess the X-ray energy effects on TBS. Multiple linear regressions (backward) were used to evaluate relationships between TBS and 3D microarchitecture parameters using a bootstrap process. BV/TV of the sample ranged from 18.5 to 37.6% with an average value at 28.8%. Correlations' analysis showedthat TBSwere strongly correlatedwith ConnD(0.856≤r≤0.862; p<0.001),with TbN (0.805≤r≤0.810; p<0.001) and negatively with TbSp (−0.714≤r≤−0.726; p<0.001), regardless X-ray energy. Results show that lower TBS values are related to "degraded" microarchitecture, with low ConnD, low TbN and a high TbSp. The opposite is also true. X-ray energy has no effect onTBS neither on the correlations betweenTBS and the 3Dmicroarchitecture parameters. In this study, we demonstrated that TBS was significantly correlated with 3D microarchitecture parameters ConnD and TbN, and negatively with TbSp, no matter what X-ray energy has been used. This article is part of a Special Issue entitled ECTS 2011. Disclosure of interest: None declared.

Prediction of fat-free mass using bioelectrical impedance analysis in young adults from five populations of African origin.

Background/objectives:Bioelectrical impedance analysis (BIA) is used in population and clinical studies as a technique for estimating body composition. Because of significant under-representation in existing literature, we sought to develop and validate predictive equation(s) for BIA for studies in populations of African origin.Subjects/methods:Among five cohorts of the Modeling the Epidemiologic Transition Study, height, weight, waist circumference and body composition, using isotope dilution, were measured in 362 adults, ages 25-45 with mean body mass indexes ranging from 24 to 32. BIA measures of resistance and reactance were measured using tetrapolar placement of electrodes and the same model of analyzer across sites (BIA 101Q, RJL Systems). Multiple linear regression analysis was used to develop equations for predicting fat-free mass (FFM), as measured by isotope dilution; covariates included sex, age, waist, reactance and height(2)/resistance, along with dummy variables for each site. Developed equations were then tested in a validation sample; FFM predicted by previously published equations were tested in the total sample.Results:A site-combined equation and site-specific equations were developed. The mean differences between FFM (reference) and FFM predicted by the study-derived equations were between 0.4 and 0.6âeuro0/00kg (that is, 1% difference between the actual and predicted FFM), and the measured and predicted values were highly correlated. The site-combined equation performed slightly better than the site-specific equations and the previously published equations.Conclusions:Relatively small differences exist between BIA equations to estimate FFM, whether study-derived or published equations, although the site-combined equation performed slightly better than others. The study-derived equations provide an important tool for research in these understudied populations.

Quantitative sequence analysis of FBN1 premature termination codons provides evidence for incomplete NMD in leukocytes.

We improved, evaluated, and used Sanger sequencing for quantification of single nucleotide polymorphism (SNP) variants in transcripts and gDNA samples. This improved assay resulted in highly reproducible relative allele frequencies (e.g., for a heterozygous gDNA 50.0+/-1.4%, and for a missense mutation-bearing transcript 46.9+/-3.7%) with a lower detection limit of 3-9%. It provided excellent accuracy and linear correlation between expected and observed relative allele frequencies. This sequencing assay, which can also be used for the quantification of copy number variations (CNVs), methylations, mosaicisms, and DNA pools, enabled us to analyze transcripts of the FBN1 gene in fibroblasts and blood samples of patients with suspected Marfan syndrome not only qualitatively but also quantitatively. We report a total of 18 novel and 19 known FBN1 sequence variants leading to a premature termination codon (PTC), 26 of which we analyzed by quantitative sequencing both at gDNA and cDNA levels. The relative amounts of PTC-containing FBN1 transcripts in fresh and PAXgene-stabilized blood samples were significantly higher (33.0+/-3.9% to 80.0+/-7.2%) than those detected in affected fibroblasts with inhibition of nonsense-mediated mRNA decay (NMD) (11.0+/-2.1% to 25.0+/-1.8%), whereas in fibroblasts without NMD inhibition no mutant alleles could be detected. These results provide evidence for incomplete NMD in leukocytes and have particular importance for RNA-based analyses not only in FBN1 but also in other genes.