Digital spiral analysis for objective assessment of fine motor timing variability in Parkinson's disease

OBJECTIVES: To develop a method for objective assessment of fine motor timing variability in Parkinson’s disease (PD) patients, using digital spiral data gathered by a touch screen device. BACKGROUND: A retrospective analysis was conducted on data from 105 subjects including65 patients with advanced PD (group A), 15 intermediate patients experiencing motor fluctuations (group I), 15 early stage patients (group S), and 10 healthy elderly subjects (HE) were examined. The subjects were asked to perform repeated upper limb motor tasks by tracing a pre-drawn Archimedes spiral as shown on the screen of the device. The spiral tracing test was performed using an ergonomic pen stylus, using dominant hand. The test was repeated three times per test occasion and the subjects were instructed to complete it within 10 seconds. Digital spiral data including stylus position (x-ycoordinates) and timestamps (milliseconds) were collected and used in subsequent analysis. The total number of observations with the test battery were as follows: Swedish group (n=10079), Italian I group (n=822), Italian S group (n = 811), and HE (n=299). METHODS: The raw spiral data were processed with three data processing methods. To quantify motor timing variability during spiral drawing tasks Approximate Entropy (APEN) method was applied on digitized spiral data. APEN is designed to capture the amount of irregularity or complexity in time series. APEN requires determination of two parameters, namely, the window size and similarity measure. In our work and after experimentation, window size was set to 4 and similarity measure to 0.2 (20% of the standard deviation of the time series). The final score obtained by APEN was normalized by total drawing completion time and used in subsequent analysis. The score generated by this method is hence on denoted APEN. In addition, two more methods were applied on digital spiral data and their scores were used in subsequent analysis. The first method was based on Digital Wavelet Transform and Principal Component Analysis and generated a score representing spiral drawing impairment. The score generated by this method is hence on denoted WAV. The second method was based on standard deviation of frequency filtered drawing velocity. The score generated by this method is hence on denoted SDDV. Linear mixed-effects (LME) models were used to evaluate mean differences of the spiral scores of the three methods across the four subject groups. Test-retest reliability of the three scores was assessed after taking mean of the three possible correlations (Spearman’s rank coefficients) between the three test trials. Internal consistency of the methods was assessed by calculating correlations between their scores. RESULTS: When comparing mean spiral scores between the four subject groups, the APEN scores were different between HE subjects and three patient groups (P=0.626 for S group with 9.9% mean value difference, P=0.089 for I group with 30.2%, and P=0.0019 for A group with 44.1%). However, there were no significant differences in mean scores of the other two methods, except for the WAV between the HE and A groups (P<0.001). WAV and SDDV were highly and significantly correlated to each other with a coefficient of 0.69. However, APEN was not correlated to neither WAV nor SDDV with coefficients of 0.11 and 0.12, respectively. Test-retest reliability coefficients of the three scores were as follows: APEN (0.9), WAV(0.83) and SD-DV (0.55). CONCLUSIONS: The results show that the digital spiral analysis-based objective APEN measure is able to significantly differentiate the healthy subjects from patients at advanced level. In contrast to the other two methods (WAV and SDDV) that are designed to quantify dyskinesias (over-medications), this method can be useful for characterizing Off symptoms in PD. The APEN was not correlated to none of the other two methods indicating that it measures a different construct of upper limb motor function in PD patients than WAV and SDDV. The APEN also had a better test-retest reliability indicating that it is more stable and consistent over time than WAV and SDDV.

Automatic spiral analysis for objective assessment of motor symptoms in Parkinson's disease

Objective: To develop a method for objective quantification of PD motor symptoms related to Off episodes and peak dose dyskinesias, using spiral data gathered by using a touch screen telemetry device. The aim was to objectively characterize predominant motor phenotypes (bradykinesia and dyskinesia), to help in automating the process of visual interpretation of movement anomalies in spirals as rated by movement disorder specialists. Background: A retrospective analysis was conducted on recordings from 65 patients with advanced idiopathic PD from nine different clinics in Sweden, recruited from January 2006 until August 2010. In addition to the patient group, 10 healthy elderly subjects were recruited. Upper limb movement data were collected using a touch screen telemetry device from home environments of the subjects. Measurements with the device were performed four times per day during week-long test periods. On each test occasion, the subjects were asked to trace pre-drawn Archimedean spirals, using the dominant hand. The pre-drawn spiral was shown on the screen of the device. The spiral test was repeated three times per test occasion and they were instructed to complete it within 10 seconds. The device had a sampling rate of 10Hz and measured both position and time-stamps (in milliseconds) of the pen tip. Methods: Four independent raters (FB, DH, AJ and DN) used a web interface that animated the spiral drawings and allowed them to observe different kinematic features during the drawing process and to rate task performance. Initially, a number of kinematic features were assessed including ‘impairment’, ‘speed’, ‘irregularity’ and ‘hesitation’ followed by marking the predominant motor phenotype on a 3-category scale: tremor, bradykinesia and/or choreatic dyskinesia. There were only 2 test occasions for which all the four raters either classified them as tremor or could not identify the motor phenotype. Therefore, the two main motor phenotype categories were bradykinesia and dyskinesia. ‘Impairment’ was rated on a scale from 0 (no impairment) to 10 (extremely severe) whereas ‘speed’, ‘irregularity’ and ‘hesitation’ were rated on a scale from 0 (normal) to 4 (extremely severe). The proposed data-driven method consisted of the following steps. Initially, 28 spatiotemporal features were extracted from the time series signals before being presented to a Multilayer Perceptron (MLP) classifier. The features were based on different kinematic quantities of spirals including radius, angle, speed and velocity with the aim of measuring the severity of involuntary symptoms and discriminate between PD-specific (bradykinesia) and/or treatment-induced symptoms (dyskinesia). A Principal Component Analysis was applied on the features to reduce their dimensions where 4 relevant principal components (PCs) were retained and used as inputs to the MLP classifier. Finally, the MLP classifier mapped these components to the corresponding visually assessed motor phenotype scores for automating the process of scoring the bradykinesia and dyskinesia in PD patients whilst they draw spirals using the touch screen device. For motor phenotype (bradykinesia vs. dyskinesia) classification, the stratified 10-fold cross validation technique was employed. Results: There were good agreements between the four raters when rating the individual kinematic features with intra-class correlation coefficient (ICC) of 0.88 for ‘impairment’, 0.74 for ‘speed’, 0.70 for ‘irregularity’, and moderate agreements when rating ‘hesitation’ with an ICC of 0.49. When assessing the two main motor phenotype categories (bradykinesia or dyskinesia) in animated spirals the agreements between the four raters ranged from fair to moderate. There were good correlations between mean ratings of the four raters on individual kinematic features and computed scores. The MLP classifier classified the motor phenotype that is bradykinesia or dyskinesia with an accuracy of 85% in relation to visual classifications of the four movement disorder specialists. The test-retest reliability of the four PCs across the three spiral test trials was good with Cronbach’s Alpha coefficients of 0.80, 0.82, 0.54 and 0.49, respectively. These results indicate that the computed scores are stable and consistent over time. Significant differences were found between the two groups (patients and healthy elderly subjects) in all the PCs, except for the PC3. Conclusions: The proposed method automatically assessed the severity of unwanted symptoms and could reasonably well discriminate between PD-specific and/or treatment-induced motor symptoms, in relation to visual assessments of movement disorder specialists. The objective assessments could provide a time-effect summary score that could be useful for improving decision-making during symptom evaluation of individualized treatment when the goal is to maximize functional On time for patients while minimizing their Off episodes and troublesome dyskinesias.

Ecotoxicological assessment of sediments from the Santos and São Vicente estuarine system- Brazil

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Discrimination of geomorphic surfaces with multivariate analysis of soil attributes in sandstone - basalt lithosequence

O conceito de superfície geomórfica permite uma interligação entre os diferentes ramos da ciência do solo, tais como geologia, geomorfologia e pedologia. Esta associação favorece a compreensão da distribuição espacial dos solos na paisagem, e torna possível compreender o comportamento dos atributos do solo, que estão principalmente relacionadas com a estratigrafia e formas do relevo. Assim, este estudo visa à aplicação da estatística multivariada para categorizar superfícies geomórficas em uma litossequência arenito-basalto, de modo a fornecer uma base para a avaliação do solo em áreas afins. A área de estudo está localizada no município de Pereira Barreto, São Paulo, Brasil. A área escolhida possui 530 hectares, onde foram localizadas e mapeadas três superfícies geomórficas (I, II e III). Na área, 134 amostras foram coletadas nas profundidades de 0,0-0,2 m e 0,8-1,0 m, foram determinados os conteúdos de areia, silte e argila, pH em CaCl2, conteúdo de MO, P, Ca, Mg, K, Al e H+Al. Com base nos resultados, foram realizadas a análise univariada e multivariada de variância, clusters e principal componente, a fim de comparar as três superfícies geomórficas. A análise estatística univariada dos atributos do solo não foi eficiente na identificação das três superfícies geomórficas. Utilizando-se os atributos físicos e químicos do solo, as técnicas estatísticas multivariada permitiram à separação dos três grupos de corpos naturais do solo que foram equivalentes as três superfícies geomórficas mapeadas. Estes resultados são interessantes, pois demonstram a viabilidade da utilização de classificação numérica das superfícies geomórficas para ajudar no mapeamento de solo.

Application of chemometric methods in the evaluation of chemical and spectroscopic data on organic matter from Oxisols in sewage sludge applications

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Structural damage detection by fuzzy clustering

The development of strategies for structural health monitoring (SHM) has become increasingly important because of the necessity of preventing undesirable damage. This paper describes an approach to this problem using vibration data. It involves a three-stage process: reduction of the time-series data using principle component analysis (PCA), the development of a data-based model using an auto-regressive moving average (ARMA) model using data from an undamaged structure, and the classification of whether or not the structure is damaged using a fuzzy clustering approach. The approach is applied to data from a benchmark structure from Los Alamos National Laboratory, USA. Two fuzzy clustering algorithms are compared: fuzzy c-means (FCM) and Gustafson-Kessel (GK) algorithms. It is shown that while both fuzzy clustering algorithms are effective, the GK algorithm marginally outperforms the FCM algorithm. (C) 2008 Elsevier Ltd. All rights reserved.

Hidrocarbonetos policíclicos aromáticos no meio ambiente: diferenciação de fontes em sedimentos e metabólitos em bile de peixes

Many studies on environmental ecosystems quality related to polycyclic aromatic hydrocarbons (PAH) have been carried out routinely due to their ubiquotus presence worldwide and to their potential toxicity after its biotransformation. PAH may be introduced into the environmet by natural and anthropogenic processes from direct runoff and discharges and indirect atmospheric deposition. Sources of naturally occurring PAHs include natural fires, natural oil seepage and recent biological or diagenetic processes. Anthropogenic sources of PAHs, acute or chronic, are combustion of organic matter (petroleum, coal, wood), waste and releases/spills of petroleum and derivatives (river runoff, sewage outfalls, maritime transport, pipelines). Besides the co-existence of multiples sources of PAH in the environmental samples, these compounds are subject to many processes that lead to geochemical fates (physical-chemical transformation, biodegradation and photo-oxidation), which leads to an alteration of their composition. All these facts make the identification of the hydrocarbons sources, if petrogenic, pyrolytic or natural, a challenge. One of the objectives of this study is to establish tools to identify the origin of hydrocarbons in environmental samples. PAH diagnostic ratios and PAH principal component analysis were tested on a critical area: Guanabara Bay sediments. Guanabara Bay is located in a complex urban area of Rio de Janeiro with a high anthropogenic influence, being an endpoint of chronic pollution from the Greater Rio and it was the scenario of an acute event of oil release in January 2000. It were quantified 38 compounds, parental and alkylated PAH, in 21 sediment samples collected in two surveys: 2000 and 2003. The PAH levels varied from 400 to 58439 ng g-1. Both tested techniques for origin identification of hydrocarbons have shown their applicability, being able to discriminate the PAH sources for the majority of the samples analysed. The bay sediments were separated into two big clusters: sediments with a clear pattern of petrogenic introduction of hydrocarbons (from intertidal area) and sediments with combustion characteristics (from subtidal region). Only a minority of the samples could not display a clear contribution of petrogenic or pyrolytic input. The diagnostic ratios that have exhibited high ability to distinguish combustion- and petroleum-derived PAH inputs for Guanabara Bay sediments were Phenanthrene+Anthracene/(Phenanthrene+Anthracene+C1Phenanthrene); Fluorantene/(Fluorantene+Pyrene); Σ (other 3-6 ring PAHs)/ Σ (5 alkylated PAH series). The PCA results prooved to be a useful tool for PAH source identification in the environment, corroborating the diagnostic indexes. In relation to the temporal evaluation carried out in this study, it was not verified significant changes on the class of predominant source of the samples. This result indicates that the hydrocarbons present in the Guanabara Bay sediments are mainly related to the long-term anthropogenic input and not directly related to acute events such as the oil spill of January 2000. This findings were similar to various international estuarine sites. Finally, this work had a complementary objective of evaluating the level of hydrocarbons exposure of the aquatic organisms of Guanabara Bay. It was a preliminary study in which a quantification of 12 individual biliar metabolites of PAH was performed in four demersal fish representing three different families. The analysed metabolites were 1-hydroxynaphtalene, 2-hidroxinaphtalene, 1hydroxyphenanthrene, 9-hydroxyphenanthrene, 2-hydroxyphenanthrene, 1hydroxypyrene, 3-hidroxibiphenil, 3- hydroxyphenanthrene, 1-hydroxychrysene, 9hydroxyfluorene, 4-hydroxyphenanthrene, 3-hydroxybenz(a)pyrene. The metabolites concentrations were found to be high, ranging from 13 to 177 µg g-1, however they were similar to worldwide regions under high anthropogenic input. Besides the metabolites established by the used protocol, it was possible to verified high concentrations of three other compounds not yet reported in the literature. They were related to pyrolytic PAH contribution to Guanabara Bay aquatic biota: 1-hydroxypyrine and 3-hydroxybenz(a)pyrine isomers

Fatores abióticos que afetam a distribuição do gênero Astyanax baird & Girard, 1854 em riachos de cabeceiras de Chapada dos Guimarães, bacia do Rio Cuiabá, Mato Grosso

The main factors that influence the species distribution are related mainly with the dispersion capacity, behavior, presence of other species and to the factors physicist-chemistries. Had been evaluated ecological factors that affect the distribution of the genus Astyanax in Chapada dos Guimaraes headstreams, Cuiaba River basin, Mato Grosso, Brazil. Samples in 27 streams had been made in Chapada dos Guimaraes, and the fish had been collected through the method of electric fish, and the ambient data had been analyzed with one technique of called ordinance of Principal Components Analysis (PCA). 250 individuals had been collect, being 95 A. asuncionensis, 62 A. abramis, 36 A. lineatus and 57 individuals had been captured identified like A. scabripinnis (Jenyns, 1842) sensu Eigenmann, 1927. In relation the ambient characteristics, the variables gotten allowed two axes of the PCA, who explain 53,114 % of the variance of the seven variables analyzed. Axle 1, Main Component 1, explained 31,147 % and Axle 2, Main Component 2, explained 21,967 % of the variation of the ambient data. The formation of axle 1 was influenced positively by the depth and negative by the conductivity and pH. The formation of axle 2 was influenced positively by the Turbidity and negative to the Dissolved Oxygen. The ordinance of the species throughout for the axes corroborated for the understanding of the distribution of specimens.

Variabilidade espacial e temporal de parâmetros físico-químicos nos rios Turvo, Preto e Grande no estado de São Paulo, Brasil

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

(1)H NMR Fingerprinting of Brazilian Commercial Gasoline: Pattern-Recognition Analyses for Origin Authentication Purposes

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Analysis of floristic composition and structure as an aid to monitoring protected areas of dense rain forest in southeastern Brazil

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Conformational analysis: A new approach by means of chemometrics

In conformational analysis, the systematic search method completely maps the space but suffers from the combinatorial explosion problem because the number of conformations increases exponentially with the number of free rotation angles. This study introduces a new methodology of conformational analysis that controls the combinatorial explosion. It is based on a dimensional reduction of the system through the use of principal component analysis. The results are exactly the same as those obtained for the complete search but, in this case, the number of conformations increases only quadratically with the number of free rotation angles. The method is applied to a series of three drugs: omeprazole. pantoprazole, lansoprazole-benzimidazoles that suppress gastric-acid secretion by means of H(+), K(+)-ATPase enzyme inhibition. (C) 2002 John Wiley Sons. Inc.

Comparative sediment quality assessment in different littoral ecosystems from Spain (Gulf of Cadiz) and Brazil (Santos and Sao Vicente estuarine system)

The goal of this work was to establish comparisons among environmental degradation in different areas from Southern Spain (Gulf of Cadiz) and Brazil (Santos and Sao Vicente estuary), by using principal component analyses (PCA) to integrate sediment toxicity (amphipods mortality) and chemical-physical data (Zn, Cd, Pb; Cu, Ni, Co, V, PCBs, PAHs concentrations, OC and fines contents). The results of PCA extraction of Spanish data showed that Bay of Cadiz, CA-1 did not present contamination or degradation; CA-2 exhibited contamination by PCBs, however it was not related to the amphipods mortality. Ria of Huelva was the most impacted site, showing contamination caused principally by hydrocarbons, in HV-1 and HV-2, but heavy metals were also important contaminants at HV-1, HV-2 and HV-3. Algeciras Bay was considered as not degraded in GR-3 and -4, but in GR-3' high contamination by PAHs was found. In the Brazilian area, the most degraded sediments were found in the stations situated at the inner parts of the estuary (SSV-2, SSV-3, and SSV-4), followed by SSV-6, which is close to the Submarine Sewage Outfall of Santos - SSOS. Sediments from SSV-1 and SSV-5 did not present chemical contamination, organic contamination or significant amphipod mortality. The results, of this investigation showed that both countries present environmental degradation related to PAHs: in Spain, at Ria of Huelva and Gudarranque river's estuary areas; and in Brasil, in the internal portion of the Santos and Sao Vicente estuary. The same situation is found for heavy metals, since all of the identified metals are related to toxicity in the studied areas, with few exceptions (V for both Brazil and Spain, and Cd and Co for Brazilian areas). The contamination by PCBs is more serious for Santos and Sao Vicente estuary than for the investigated areas in Gulf of Cadiz, where such compound did not relate to the toxicity. (c) 2006 Elsevier Ltd. All rights reserved.

A partial least squares and principal component regression study of quinone compounds with trypanocidal activity

A quantitative structure-activity relationship (QSAR) study of 19 quinone compounds with trypanocidal activity was performed by Partial Least Squares (PLS) and Principal Component Regression (PCR) methods with the use of leave-one-out crossvalidation procedure to build the regression models. The trypanocidal activity of the compounds is related to their first cathodic potential (Ep(c1)). The regression PLS and PCR models built in this study were also used to predict the Ep(c1) of six new quinone compounds. The PLS model was built with three principal components that described 96.50% of the total variance and present Q(2) = 0.83 and R-2 = 0.90. The results obtained with the PCR model were similar to those obtained with the PLS model. The PCR model was also built with three principal components that described 96.67% of the total variance with Q(2) = 0.83 and R-2 = 0.90. The most important descriptors for our PLS and PCR models were HOMO-1 (energy of the molecular orbital below HOMO), Q4 (atomic charge at position 4), MAXDN (maximal electrotopological negative difference), and HYF (hydrophilicity index).

Allometric analysis of the ontogenetic variation and evolution of the skull in Caiman spix, 1825 (Crocodylia: Alligatoridae)

We studied ontogenetic variation in the shape of the skull among species of Caiman using principal component analysis. Comparison of multivariate allometric coefficients and ontogenetic trends between size and shape reveals that C. sclerops and C. yacare have similar ontogenetic processes, and they are more related to each other than either is to C. latirostris. Allometric relationships of the characters measured are similar in all species studied. The greater differences were in the width measurements, with higher coefficients in shape (second principal component) for C. latirostris, and length measurements with higher coefficients in shape for C. yacare and C. sclerops. The ontogenetic process leading to change in skull shape in the group seems to be plesiomorphic for elongation and derived for broadening. Statistical comparison of the ontogenetic trends with models of allometric heterochrony suggests that C. latirostris has diverged from the other species by a neotenic process, and that C. sclerops is separated from C. yacare by ontogenetic scaling (progenesis).