Neurite, a finite difference large scale parallel program for the simulation of the electrical signal propagation in neurites under mechanical loading

With the growing body of research on traumatic brain injury and spinal cord injury, computational neuroscience has recently focused its modeling efforts on neuronal functional deficits following mechanical loading. However, in most of these efforts, cell damage is generally only characterized by purely mechanistic criteria, function of quantities such as stress, strain or their corresponding rates. The modeling of functional deficits in neurites as a consequence of macroscopic mechanical insults has been rarely explored. In particular, a quantitative mechanically based model of electrophysiological impairment in neuronal cells has only very recently been proposed (Jerusalem et al., 2013). In this paper, we present the implementation details of Neurite: the finite difference parallel program used in this reference. Following the application of a macroscopic strain at a given strain rate produced by a mechanical insult, Neurite is able to simulate the resulting neuronal electrical signal propagation, and thus the corresponding functional deficits. The simulation of the coupled mechanical and electrophysiological behaviors requires computational expensive calculations that increase in complexity as the network of the simulated cells grows. The solvers implemented in Neurite-explicit and implicit-were therefore parallelized using graphics processing units in order to reduce the burden of the simulation costs of large scale scenarios. Cable Theory and Hodgkin-Huxley models were implemented to account for the electrophysiological passive and active regions of a neurite, respectively, whereas a coupled mechanical model accounting for the neurite mechanical behavior within its surrounding medium was adopted as a link between lectrophysiology and mechanics (Jerusalem et al., 2013). This paper provides the details of the parallel implementation of Neurite, along with three different application examples: a long myelinated axon, a segmented dendritic tree, and a damaged axon. The capabilities of the program to deal with large scale scenarios, segmented neuronal structures, and functional deficits under mechanical loading are specifically highlighted.

Large-scale identification of gibberellin-related transcription factor defines group VII ERFs as functional of DELLA partners

DELLA proteins are the master negative regulators in gibberellin (GA) signaling acting in the nucleus as transcriptional regulators. The current view of DELLA action indicates that their activity relies on the physical interaction with transcription factors (TFs). Therefore, the identification of TFs through which DELLAs regulate GA responses is key to understanding these responses from a mechanistic point of view. Here, we have determined the TF interactome of the Arabidopsis (Arabidopsis thaliana) DELLA protein GIBBERELLIN INSENSITIVE and screened a collection of conditional TF overexpressors in search of those that alter GA sensitivity. As a result, we have found RELATED TO APETALA2.3, an ethylene-induced TF belonging to the group VII ETHYLENE RESPONSE FACTOR of the APETALA2/ethylene responsive element binding protein superfamily, as a DELLA interactor with physiological relevance in the context of apical hook development. The combination of transactivation assays and chromatin immunoprecipitation indicates that the interaction with GIBBERELLIN INSENSITIVE impairs the activity of RELATED TO APETALA2.3 on the target promoters. This mechanism represents a unique node in the cross regulation between the GA and ethylene signaling pathways controlling differential growth during apical hook development.

Standing together for reproducibility in large-scale computing: report on reproducibility@XSEDE

This is the final report on reproducibility@xsede, a one-day workshop held in conjunction with XSEDE14, the annual conference of the Extreme Science and Engineering Discovery Environment (XSEDE). The workshop's discussion-oriented agenda focused on reproducibility in large-scale computational research. Two important themes capture the spirit of the workshop submissions and discussions: (1) organizational stakeholders, especially supercomputer centers, are in a unique position to promote, enable, and support reproducible research; and (2) individual researchers should conduct each experiment as though someone will replicate that experiment. Participants documented numerous issues, questions, technologies, practices, and potentially promising initiatives emerging from the discussion, but also highlighted four areas of particular interest to XSEDE: (1) documentation and training that promotes reproducible research; (2) system-level tools that provide build- and run-time information at the level of the individual job; (3) the need to model best practices in research collaborations involving XSEDE staff; and (4) continued work on gateways and related technologies. In addition, an intriguing question emerged from the day's interactions: would there be value in establishing an annual award for excellence in reproducible research? Overview

The impact of large-scale circulation patterns on summer crop in Iberian Peninsula

Large-scale circulations patterns (ENSO, NAO) have been shown to have a significant impact on seasonal weather, and therefore on crop yield over many parts of the world(Garnett and Khandekar, 1992; Aasa et al., 2004; Rozas and Garcia-Gonzalez, 2012). In this study, we analyze the influence of large-scale circulation patterns and regional climate on the principal components of maize yield variability in Iberian Peninsula (IP) using reanalysis datasets. Additionally, we investigate the modulation of these relationships by multidecadal patterns. This study is performed analyzing long time series of maize yield, only climate dependent, computed with the crop model CERES-maize (Jones and Kiniry, 1986) included in Decision Support System for Agrotechnology Transfer (DSSAT v.4.5).

p-medicine: a medical informatics platform for integrated large scale heterogeneous patient data

Secure access to patient data is becoming of increasing importance, as medical informatics grows in significance, to both assist with population health studies, and patient specific medicine in support of treatment. However, assembling the many different types of data emanating from the clinic is in itself a difficulty, and doing so across national borders compounds the problem. In this paper we present our solution: an easy to use distributed informatics platform embedding a state of the art data warehouse incorporating a secure pseudonymisation system protecting access to personal healthcare data. Using this system, a whole range of patient derived data, from genomics to imaging to clinical records, can be assembled and linked, and then connected with analytics tools that help us to understand the data. Research performed in this environment will have immediate clinical impact for personalised patient healthcare.

Heterogeneous architecture for the optimization of large-scale graph processing data centres

Con el auge del Cloud Computing, las aplicaciones de proceso de datos han sufrido un incremento de demanda, y por ello ha cobrado importancia lograr m�ás eficiencia en los Centros de Proceso de datos. El objetivo de este trabajo es la obtenci�ón de herramientas que permitan analizar la viabilidad y rentabilidad de diseñar Centros de Datos especializados para procesamiento de datos, con una arquitectura, sistemas de refrigeraci�ón, etc. adaptados. Algunas aplicaciones de procesamiento de datos se benefician de las arquitecturas software, mientras que en otras puede ser m�ás eficiente un procesamiento con arquitectura hardware. Debido a que ya hay software con muy buenos resultados en el procesamiento de grafos, como el sistema XPregel, en este proyecto se realizará una arquitectura hardware en VHDL, implementando el algoritmo PageRank de Google de forma escalable. Se ha escogido este algoritmo ya que podr��á ser m�ás eficiente en arquitectura hardware, debido a sus características concretas que se indicaráan m�ás adelante. PageRank sirve para ordenar las p�áginas por su relevancia en la web, utilizando para ello la teorí��a de grafos, siendo cada página web un vértice de un grafo; y los enlaces entre páginas, las aristas del citado grafo. En este proyecto, primero se realizará un an�álisis del estado de la técnica. Se supone que la implementaci�ón en XPregel, un sistema de procesamiento de grafos, es una de las m�ás eficientes. Por ello se estudiará esta �ultima implementaci�ón. Sin embargo, debido a que Xpregel procesa, en general, algoritmos que trabajan con grafos; no tiene en cuenta ciertas caracterí��sticas del algoritmo PageRank, por lo que la implementaci�on no es �optima. Esto es debido a que en PageRank, almacenar todos los datos que manda un mismo v�értice es un gasto innecesario de memoria ya que todos los mensajes que manda un vértice son iguales entre sí e iguales a su PageRank. Se realizará el diseño en VHDL teniendo en cuenta esta caracter��ística del citado algoritmo,evitando almacenar varias veces los mensajes que son iguales. Se ha elegido implementar PageRank en VHDL porque actualmente las arquitecturas de los sistemas operativos no escalan adecuadamente. Se busca evaluar si con otra arquitectura se obtienen mejores resultados. Se realizará un diseño partiendo de cero, utilizando la memoria ROM de IPcore de Xillinx (Software de desarrollo en VHDL), generada autom�áticamente. Se considera hacer cuatro tipos de módulos para que as�� el procesamiento se pueda hacer en paralelo. Se simplificar�á la estructura de XPregel con el fin de intentar aprovechar la particularidad de PageRank mencionada, que hace que XPregel no le saque el m�aximo partido. Despu�és se escribirá el c�ódigo, realizando una estructura escalable, ya que en la computación intervienen millones de páginas web. A continuación, se sintetizar�á y se probará el código en una FPGA. El �ultimo paso será una evaluaci�ón de la implementaci�ón, y de posibles mejoras en cuanto al consumo.

Column Generation Algorithms for Nonlinear Optimization II: Numerical Investigations

García et al. present a class of column generation (CG) algorithms for nonlinear programs. Its main motivation from a theoretical viewpoint is that under some circumstances, finite convergence can be achieved, in much the same way as for the classic simplicial decomposition method; the main practical motivation is that within the class there are certain nonlinear column generation problems that can accelerate the convergence of a solution approach which generates a sequence of feasible points. This algorithm can, for example, accelerate simplicial decomposition schemes by making the subproblems nonlinear. This paper complements the theoretical study on the asymptotic and finite convergence of these methods given in [1] with an experimental study focused on their computational efficiency. Three types of numerical experiments are conducted. The first group of test problems has been designed to study the parameters involved in these methods. The second group has been designed to investigate the role and the computation of the prolongation of the generated columns to the relative boundary. The last one has been designed to carry out a more complete investigation of the difference in computational efficiency between linear and nonlinear column generation approaches. In order to carry out this investigation, we consider two types of test problems: the first one is the nonlinear, capacitated single-commodity network flow problem of which several large-scale instances with varied degrees of nonlinearity and total capacity are constructed and investigated, and the second one is a combined traffic assignment model

MBIS: Multivariate Bayesian Image Segmentation Tool

We present MBIS (Multivariate Bayesian Image Segmentation tool), a clustering tool based on the mixture of multivariate normal distributions model. MBIS supports multi-channel bias field correction based on a B-spline model. A second methodological novelty is the inclusion of graph-cuts optimization for the stationary anisotropic hidden Markov random field model. Along with MBIS, we release an evaluation framework that contains three different experiments on multi-site data. We first validate the accuracy of segmentation and the estimated bias field for each channel. MBIS outperforms a widely used segmentation tool in a cross-comparison evaluation. The second experiment demonstrates the robustness of results on atlas-free segmentation of two image sets from scan-rescan protocols on 21 healthy subjects. Multivariate segmentation is more replicable than the monospectral counterpart on T1-weighted images. Finally, we provide a third experiment to illustrate how MBIS can be used in a large-scale study of tissue volume change with increasing age in 584 healthy subjects. This last result is meaningful as multivariate segmentation performs robustly without the need for prior knowledge.

Risk of Leaching in Soils Amended by Compost and digestate from Municipal Solid Waste

New European directives have proposed the direct application of compost and digestate produced from municipal solid wastes as organic matter sources in agricultural soils. Therefore information about phosphorus leaching from these residues when they are applied to the soil is increasingly mportant. Leaching experiments were conducted to determine the P mobility in compost and digestate mixtures, supplying equivalent amounts to 100 kg P ha?1 to three different types of soils. The tests were performed in accordance with CEN/TS 14405:2004 analyzing the maximum dissolved reactive P and the kinetic rate in the leachate. P biowaste fractionation indicated that digestate has a higher level of available P than compost has. In contrast, P losses in leaching experiments with soil-compost mixtureswere higher than in soil-digestate mixtures. For bothwastes, therewas no correlation between disolved reactive P lost and the water soluble P.The interaction between soil and waste, the long experimentation time, and the volume of leachate obtained caused the waste?s wettability to become an influential parameter in P leaching behavior. The overall conclusion is that kinetic data analysis provides valuable information concerning the sorption mechanism that can be used for predicting the large-scale behavior of soil systems.

Functional response of Ulmus minor Mill. to drought, flooding and dutch elm disease

Ulmus minor es una especie arbórea originaria de Europa cuyas poblaciones han sido diezmadas por el hongo patógeno causante de la enfermedad de la grafiosis. La conservación de los olmos exige plantearse su propagación a través de plantaciones y conocer mejor su ecología y biología. Ulmus minor es un árbol de ribera, pero frecuentemente se encuentra alejado del cauce de arroyos y ríos, donde la capa freática sufre fuertes oscilaciones. Por ello, nuestra hipótesis general es que esta especie es moderadamente resistente tanto a la inundación como a la sequía. El principal objetivo de esta tesis doctoral es entender desde un punto de vista funcional la respuesta de U. minor a la inundación, la sequía y la infección por O. novo-ulmi; los factores que posiblemente más influyen en la distribución actual de U. minor. Con este objetivo se persigue dar continuidad a los esfuerzos de conservación de esta especie que desde hace años se dedican en varios centros de investigación a nivel mundial, ya que, entender mejor los mecanismos que contribuyen a la resistencia de U. minor ante la inoculación con O. novo-ulmi y factores de estrés abiótico ayudará en la selección y propagación de genotipos resistentes a la grafiosis. Se han planteado tres experimentos en este sentido. Primero, se ha comparado la tolerancia de brinzales de U. minor y U. laevis – otro olmo ibérico – a una inmersión controlada con el fin de evaluar su tolerancia a la inundación y comprender los mecanismos de aclimatación. Segundo, se ha comparado la tolerancia de brinzales de U. minor y Quercus ilex – una especie típica de ambientes Mediterránea secos – a la falta de agua en el suelo con el fin de evaluar el grado de tolerancia y los mecanismos de aclimatación a la sequía. El hecho de comparar dos especies contrastadas responde al interés en entender mejor cuales son los procesos que conducen a la muerte de una planta en condiciones de sequía – asunto sobre el que hay una interesante discusión desde hace algunos años. En tercer lugar, con el fin de entender mejor la resistencia de algunos genotipos de U. minor a la grafiosis, se han estudiado las diferencias fisiológicas y químicas constitutivas e inducidas por O. novo-ulmi entre clones de U. minor seleccionados a priori por su variable grado de resistencia a esta enfermedad. En el primer experimento se observó que los brinzales de U. minor sobrevivieron 60 días inmersos en una piscina con agua no estancada hasta una altura de 2-3 cm por encima del cuello de la raíz. A los 60 días, los brinzales de U. laevis se sacaron de la piscina y, a lo largo de las siguientes semanas, fueron capaces de recuperar las funciones fisiológicas que habían sido alteradas anteriormente. La conductividad hidráulica de las raíces y la tasa de asimilación de CO2 neta disminuyeron en ambas especies. Por el contrario, la tasa de respiración de hojas, tallos y raíces aumentó en las primeras semanas de la inundación, posiblemente en relación al aumento de energía necesario para desarrollar mecanismos de aclimatación a la inundación, como la hipertrofia de las lenticelas que se observó en ambas especies. Por ello, el desequilibrio del balance de carbono de la planta podría ser un factor relevante en la mortalidad de las plantas ante inundaciones prolongadas. Las plantas de U. minor (cultivadas en envases de 16 litros a media sombra) sobrevivieron por un prolongado periodo de tiempo en verano sin riego; la mitad de las plantas murieron tras 90 días sin riego. El cierre de los estomas y la pérdida de hojas contribuyeron a ralentizar las pérdidas de agua y tolerar la sequía en U. minor. Las obvias diferencias en tolerancia a la sequía con respecto a Q. ilex se reflejaron en la distinta capacidad para ralentizar la aparición del estrés hídrico tras dejar de regar y para transportar agua en condiciones de elevada tensión en el xilema. Más relevante es que las plantas con evidentes síntomas de decaimiento previo a su muerte exhibieron pérdidas de conductividad hidráulica en las raíces del 80% en ambas especies, mientras que las reservas de carbohidratos apenas variaron y lo hicieron de forma desigual en ambas especies. Árboles de U. minor de 5 y 6 años de edad (plantados en eras con riego mantenido) exhibieron una respuesta a la inoculación con O. novo-ulmi consistente con ensayos previos de resistencia. La conductividad hidráulica del tallo, el potencial hídrico foliar y la tasa de asimilación de CO2 neta disminuyeron significativamente en relación a árboles inoculados con agua, pero solo en los clones susceptibles. Este hecho enlaza con el perfil químico “más defensivo” de los clones resistentes, es decir, con los mayores niveles de suberina, ácidos grasos y compuestos fenólicos en estos clones que en los susceptibles. Ello podría restringir la propagación del hongo en el árbol y preservar el comportamiento fisiológico de los clones resistentes al inocularlos con el patógeno. Los datos indican una respuesta fisiológica común de U. minor a la inundación, la sequía y la infección por O. novo-ulmi: pérdida de conductividad hidráulica, estrés hídrico y pérdida de ganancia neta de carbono. Pese a ello, U. minor desarrolla varios mecanismos que le confieren una capacidad moderada para vivir en suelos temporalmente anegados o secos. Por otro lado, el perfil químico es un factor relevante en la resistencia de ciertos genotipos a la grafiosis. Futuros estudios deberían examinar como este perfil químico y la resistencia a la grafiosis se ven alteradas por el estrés abiótico. ABSTRACT Ulmus minor is a native European elm species whose populations have been decimated by the Dutch elm disease (DED). An active conservation of this species requires large-scale plantations and a better understanding of its biology and ecology. U. minor generally grows close to water channels. However, of the Iberian riparian tree species, U. minor is the one that spread farther away from rivers and streams. For these reasons, we hypothesize that this species is moderately tolerant to both flooding and drought stresses. The main aim of the present PhD thesis is to better understand the functional response of U. minor to the abiotic stresses – flooding and drought – and the biotic stress – DED – that can be most influential on its distribution. The overarching goal is to aid in the conservation of this emblematic species through a better understanding of the mechanisms that contribute to resistance to abiotic and biotic stresses; an information that can help in the selection of resistant genotypes and their expansion in large-scale plantations. To this end, three experiments were set up. First, we compared the tolerance to experimental immersion between seedlings of U. minor and U. laevis – another European riparian elm species – in order to assess their degree of tolerance and understand the mechanisms of acclimation to this stress. Second, we investigated the tolerance to drought of U. minor seedlings in comparison with Quercus ilex (an oak species typical of dry Mediterranean habitats). Besides assessing and understanding U. minor tolerance to drought at the seedling stage, the aim was to shed light into the functional alterations that trigger drought-induced plant mortality – a matter of controversy in the last years. Third, we studied constitutive and induced physiological and biochemical differences among clones of variable DED resistance, before and following inoculation with Ophiostoma novo-ulmi. The goal is to shed light into the factors of DED resistance that is evident in some genotypes of U. minor, but not others. Potted seedlings of U. minor survived for 60 days immersed in a pool with running water to approximately 2-3 cm above the stem collar. By this time, U. minor seedlings died, whereas U. laevis seedlings moved out of the pool were able to recover most physiological functions that had been altered by flooding. For example, root hydraulic conductivity and leaf photosynthetic CO2 uptake decreased in both species; while respiration initially increased with flooding in leaves, stems and roots possibly to respond to energy demands associated to mechanisms of acclimation to soil oxygen deficiency; as example, a remarkable hypertrophy of lenticels was soon observed in flooded seedlings of both species. Therefore, the inability to maintain a positive carbon balance somehow compromises seedling survival under flooding, earlier in U. minor than U. laevis, partly explaining their differential habitats. Potted seedlings of U. minor survived for a remarkable long time without irrigation – half of plants dying only after 90 days of no irrigation in conditions of high vapour pressure deficit typical of summer. Some mechanisms that contributed to tolerate drought were leaf shedding and stomata closure, which reduced water loss and the risk of xylem cavitation. Obviously, U. minor was less tolerant to drought than Q. ilex, differences in drought tolerance resulting mostly from the distinct capacity to postpone water stress and conduct water under high xylem tension among species. More relevant was that plants of both species exhibited similar symptoms of root hydraulic failure (i.e. approximately 80% loss of hydraulic conductivity), but a slight and variable depletion of non-structural carbohydrate reserves preceding dieback. Five- and six-year-old trees of U. minor (planted in the field with supplementary watering) belonging to clones of contrasted susceptibility to DED exhibited a different physiological response to inoculation with O. novo-ulmi. Stem hydraulic conductivity, leaf water potential and photosynthetic CO2 uptake decreased significantly relative to control trees inoculated with water only in DED susceptible clones. This is consistent with the “more defensive” chemical profile observed in resistant clones, i.e. with higher levels of saturated hydrocarbons (suberin and fatty acids) and phenolic compounds than in susceptible clones. These compounds could restrict the spread of O. novo-ulmi and contribute to preserving the near-normal physiological function of resistant trees when exposed to the pathogen. These results evidence common physiological responses of U. minor to flooding, drought and pathogen infection leading to xylem water disruption, leaf water stress and reduced net carbon gain. Still, seedlings of U. minor develop various mechanisms of acclimation to abiotic stresses that can play a role in surviving moderate periods of flood and drought. The chemical profile appears to be an important factor for the resistance of some genotypes of U. minor to DED. How abiotic stresses such as flooding and drought affect the capacity of resistant U. minor clones to face O. novo-ulmi is a key question that must be contemplated in future research.