Homolonto: generating homology relationships by pairwise alignment of ontologies and application to vertebrate anatomy.

MOTIVATION: The anatomy of model species is described in ontologies, which are used to standardize the annotations of experimental data, such as gene expression patterns. To compare such data between species, we need to establish relations between ontologies describing different species. RESULTS: We present a new algorithm, and its implementation in the software Homolonto, to create new relationships between anatomical ontologies, based on the homology concept. Homolonto uses a supervised ontology alignment approach. Several alignments can be merged, forming homology groups. We also present an algorithm to generate relationships between these homology groups. This has been used to build a multi-species ontology, for the database of gene expression evolution Bgee. AVAILABILITY: download section of the Bgee website http://bgee.unil.ch/

Systems biology of plants : from intraspecific genome variation to computational morphodynamics

Recently, the introduction of second generation sequencing and further advance-ments in confocal microscopy have enabled system-level studies for the functional characterization of genes. The degree of complexity intrinsic to these approaches needs the development of bioinformatics methodologies and computational models for extracting meaningful biological knowledge from the enormous amount of experi¬mental data which is continuously generated. This PhD thesis presents several novel bioinformatics methods and computational models to address specific biological questions in Plant Biology by using the plant Arabidopsis thaliana as a model system. First, a spatio-temporal qualitative analysis of quantitative transcript and protein profiles is applied to show the role of the BREVIS RADIX (BRX) protein in the auxin- cytokinin crosstalk for root meristem growth. Core of this PhD work is the functional characterization of the interplay between the BRX protein and the plant hormone auxin in the root meristem by using a computational model based on experimental evidence. Hyphotesis generated by the modelled to the discovery of a differential endocytosis pattern in the root meristem that splits the auxin transcriptional response via the plasma membrane to nucleus partitioning of BRX. This positional information system creates an auxin transcriptional pattern that deviates from the canonical auxin response and is necessary to sustain the expression of a subset of BRX-dependent auxin-responsive genes to drive root meristem growth. In the second part of this PhD thesis, we characterized the genome-wide impact of large scale deletions on four divergent Arabidopsis natural strains, through the integration of Ultra-High Throughput Sequencing data with data from genomic hybridizations on tiling arrays. Analysis of the identified deletions revealed a considerable portion of protein coding genes affected and supported a history of genomic rearrangements shaped by evolution. In the last part of the thesis, we showed that VIP3 gene in Arabidopsis has an evo-lutionary conserved role in the 3' to 5' mRNA degradation machinery, by applying a novel approach for the analysis of mRNA-Seq data from random-primed mRNA. Altogether, this PhD research contains major advancements in the study of natural genomic variation in plants and in the application of computational morphodynamics models for the functional characterization of biological pathways essential for the plant. - Récemment, l'introduction du séquençage de seconde génération et les avancées dans la microscopie confocale ont permis des études à l'échelle des différents systèmes cellulaires pour la caractérisation fonctionnelle de gènes. Le degrés de complexité intrinsèque à ces approches ont requis le développement de méthodologies bioinformatiques et de modèles mathématiques afin d'extraire de la masse de données expérimentale générée, des information biologiques significatives. Ce doctorat présente à la fois des méthodes bioinformatiques originales et des modèles mathématiques pour répondre à certaines questions spécifiques de Biologie Végétale en utilisant la plante Arabidopsis thaliana comme modèle. Premièrement, une analyse qualitative spatio-temporelle de profiles quantitatifs de transcripts et de protéines est utilisée pour montrer le rôle de la protéine BREVIS RADIX (BRX) dans le dialogue entre l'auxine et les cytokinines, des phytohormones, dans la croissance du méristème racinaire. Le noyau de ce travail de thèse est la caractérisation fonctionnelle de l'interaction entre la protéine BRX et la phytohormone auxine dans le méristème de la racine en utilisant des modèles informatiques basés sur des preuves expérimentales. Les hypothèses produites par le modèle ont mené à la découverte d'un schéma différentiel d'endocytose dans le méristème racinaire qui divise la réponse transcriptionnelle à l'auxine par le partitionnement de BRX de la membrane plasmique au noyau de la cellule. Cette information positionnelle crée une réponse transcriptionnelle à l'auxine qui dévie de la réponse canonique à l'auxine et est nécessaire pour soutenir l'expression d'un sous ensemble de gènes répondant à l'auxine et dépendant de BRX pour conduire la croissance du méristème. Dans la seconde partie de cette thèse de doctorat, nous avons caractérisé l'impact sur l'ensemble du génome des délétions à grande échelle sur quatre souches divergentes naturelles d'Arabidopsis, à travers l'intégration du séquençage à ultra-haut-débit avec l'hybridation génomique sur puces ADN. L'analyse des délétions identifiées a révélé qu'une proportion considérable de gènes codant était affectée, supportant l'idée d'un historique de réarrangement génomique modelé durant l'évolution. Dans la dernière partie de cette thèse, nous avons montré que le gène VÏP3 dans Arabidopsis a conservé un rôle évolutif dans la machinerie de dégradation des ARNm dans le sens 3' à 5', en appliquant une nouvelle approche pour l'analyse des données de séquençage d'ARNm issue de transcripts amplifiés aléatoirement. Dans son ensemble, cette recherche de doctorat contient des avancées majeures dans l'étude des variations génomiques naturelles des plantes et dans l'application de modèles morphodynamiques informatiques pour la caractérisation de réseaux biologiques essentiels à la plante. - Le développement des plantes est écrit dans leurs codes génétiques. Pour comprendre comment les plantes sont capables de s'adapter aux changements environnementaux, il est essentiel d'étudier comment leurs gènes gouvernent leur formation. Plus nous essayons de comprendre le fonctionnement d'une plante, plus nous réalisons la complexité des mécanismes biologiques, à tel point que l'utilisation d'outils et de modèles mathématiques devient indispensable. Dans ce travail, avec l'utilisation de la plante modèle Arabidopsis thalicinci nous avons résolu des problèmes biologiques spécifiques à travers le développement et l'application de méthodes informatiques concrètes. Dans un premier temps, nous avons investigué comment le gène BREVIS RADIX (BRX) régule le développement de la racine en contrôlant la réponse à deux hormones : l'auxine et la cytokinine. Nous avons employé une analyse statistique sur des mesures quantitatives de transcripts et de produits de gènes afin de démontrer que BRX joue un rôle antagonisant dans le dialogue entre ces deux hormones. Lorsque ce-dialogue moléculaire est perturbé, la racine primaire voit sa longueur dramatiquement réduite. Pour comprendre comment BRX répond à l'auxine, nous avons développé un modèle informatique basé sur des résultats expérimentaux. Les simulations successives ont mené à la découverte d'un signal positionnel qui contrôle la réponse de la racine à l'auxine par la régulation du mouvement intracellulaire de BRX. Dans la seconde partie de cette thèse, nous avons analysé le génome entier de quatre souches naturelles d'Arabidopsis et nous avons trouvé qu'une grande partie de leurs gènes étaient manquant par rapport à la souche de référence. Ce résultat indique que l'historique des modifications génomiques conduites par l'évolution détermine une disponibilité différentielle des gènes fonctionnels dans ces plantes. Dans la dernière partie de ce travail, nous avons analysé les données du transcriptome de la plante où le gène VIP3 était non fonctionnel. Ceci nous a permis de découvrir le rôle double de VIP3 dans la régulation de l'initiation de la transcription et dans la dégradation des transcripts. Ce rôle double n'avait jusqu'alors été démontrée que chez l'homme. Ce travail de doctorat supporte le développement et l'application de méthodologies informatiques comme outils inestimables pour résoudre la complexité des problèmes biologiques dans la recherche végétale. L'intégration de la biologie végétale et l'informatique est devenue de plus en plus importante pour l'avancée de nos connaissances sur le fonctionnement et le développement des plantes.

Synchronous versus asynchronous modeling of gene regulatory networks.

MOTIVATION: In silico modeling of gene regulatory networks has gained some momentum recently due to increased interest in analyzing the dynamics of biological systems. This has been further facilitated by the increasing availability of experimental data on gene-gene, protein-protein and gene-protein interactions. The two dynamical properties that are often experimentally testable are perturbations and stable steady states. Although a lot of work has been done on the identification of steady states, not much work has been reported on in silico modeling of cellular differentiation processes. RESULTS: In this manuscript, we provide algorithms based on reduced ordered binary decision diagrams (ROBDDs) for Boolean modeling of gene regulatory networks. Algorithms for synchronous and asynchronous transition models have been proposed and their corresponding computational properties have been analyzed. These algorithms allow users to compute cyclic attractors of large networks that are currently not feasible using existing software. Hereby we provide a framework to analyze the effect of multiple gene perturbation protocols, and their effect on cell differentiation processes. These algorithms were validated on the T-helper model showing the correct steady state identification and Th1-Th2 cellular differentiation process. AVAILABILITY: The software binaries for Windows and Linux platforms can be downloaded from http://si2.epfl.ch/~garg/genysis.html.

Changes on Schistosoma mansoni (Digenea: Schistosomatidae) worm load in Nectomys squamipes (Rodentia: Sigmodontinae) concurrently infected with Echinostoma paraensei (Digenea: Echinostomatidae)

The water rat, Nectomys squamipes, closely involved in schistosomiasis transmission in Brazil, has been found naturally infected simultaneously by Schistosoma mansoni and Echinostoma paraensei. Laboratory experiments were conducted to verify parasitic interaction in concurrent infection. It was replicated four times with a total of 42 water rats and essayed two times with 90 mice pre-infected with E. paraensei. Rodents were divided into three groups in each replication. A wild strain recently isolated from Sumidouro, RJ, and a laboratory strain of S. mansoni from Belo Horizonte (BH) was used. Rats infected with E. paraensei were challenged 4 weeks later with S. mansoni and mice 2 or 6 weeks after the infection with S. mansoni. Necropsy took place 8 weeks following S. mansoni infection. The N. squamipes treatment groups challenged with S. mansoni RJ strain showed a significant decrease (80 and 65%) in the S. mansoni parasite load when compared with their respective control groups. There was a significant change or no change in the hosts challenged with the BH strain. The persistence time of E. paraensei within host was extended in relation to control groups, with a consequent enhancement of the number of recovered worm. An E. paraensei strain-specific influence on S. mansoni parasitism is reported. This paper presents some experimental data about this interaction in N. squamipes and Mus musculus.

Assessing key safety concerns of a Wolbachia-based strategy to control dengue transmission by Aedes mosquitoes

Mosquito-borne diseases such as dengue fever, chikungunya or malaria affect millions of people each year and control solutions are urgently needed. An international research program is currently being developed that relies on the introduction of the bacterial endosymbiont Wolbachia pipientis into Aedes aegypti to control dengue transmission. In order to prepare for open-field testing releases of Wolbachia-infected mosquitoes, an intensive social research and community engagement program was undertaken in Cairns, Northern Australia. The most common concern expressed by the diverse range of community members and stakeholders surveyed was the necessity of assuring the safety of the proposed approach for humans, animals and the environment. To address these concerns a series of safety experiments were undertaken. We report in this paper on the experimental data obtained, discuss the limitations of experimental risk assessment and focus on the necessity of including community concerns in scientific research.

Calcium, vitamin D and cardiovascular disease.

The relationship between calcium and cardiovascular diseases (CVD) has been explored for a long time. Studies exploring the effect of calcium intake or calcium supplementation on cardiovascular risk suggest that systolic blood pressure increases under low calcium intake and decreases with calcium supplementation. A lower calcium intake has been associated with an increased risk of stroke. However, the impact of calcium supplementation on stroke risk remains unclear. Calcium supplementation may increase the risk of myocardial infarction. The relationship between vitamin D and CVD has been explored more recently. Negative correlations between vitamin D levels and the risk of hypertension, myocardial infarction, and stroke have been reported in several observational studies. The effect of vitamin D supplementation on blood pressure is still unclear and no effect of vitamin D supplementation on coronary heart disease or stroke has been clearly demonstrated. There is a lack of randomized clinical trials primarily addressing the effect of these parameters on CVD. Therefore, the real impact of calcium and vitamin D on cardiovascular outcomes remains to be documented by appropriate experimental data.

Modelització de la producció de melanina a partir de polifenol oxidasa d'Agaricus bisporus

Les melanines són un grup heterogeni de polímers producte de reaccions enzimàtiques en els teixits vegetals que contenen compostos fenòlics o polifenòlics. Estudis recents han descobert algunes propietats benèfiques de les melanines sobre la salut, tals com antioxidants, antiinflamatòries, immunològiques i propietats anti-tumorals. Així, no només la seva eliminació ha de ser examinada, sinó que també podria considerar-se la seva addició a aliments funcionals de nova creació. D’aquesta manera, es requereix conèixer el mecanisme cinètic de la lanogènesi abans de la seva possible utilització industrial. S’ha desenvolupat un model cinètic per explicar la formació de melanina a partir de L-tirosina utilitzant polifenol oxidasa d’Agaricus bisporus i monitoritzant l'absorbància de la solució. Aquesta expressió permet descriure la formació de melanina en funció del temps de reacció i obtenir alguns paràmetres importants que defineixen el producte, com el coeficient d'extinció. L’absorbància comença a créixer després d'un període de latència en què es produeixen productes intermedis incolors. El coeficient d'extinció dels productes resultants no és un valor constant, perquè depèn de les condicions de cada experiment. La tirosinasa tingué un menor efecte catalitzador sobre la L-tirosina (primera reacció que catalitza), que sobre L-DOPA (segona reacció).

Insulin and NPY pathways and the control of GnRH function and puberty onset.

Energy balance exerts a critical influence on reproductive function. Leptin and insulin are among the metabolic factors signaling the nutritional status of an individual to the hypothalamus, and their role in the overall modulation of the activity of GnRH neurons is increasingly recognized. As such, they participate to a more generalized phenomenon: the signaling of peripheral metabolic changes to the central nervous system. The physiological importance that the interactions occurring between peripheral metabolic factors and the central nervous system bear for the control of food intake is increasingly recognized. The central mechanisms implicated are the focus of attention of very many research groups worldwide. We review here the experimental data that suggest that similar mechanisms are at play for the metabolic control of the neuroendocrine reproductive function. It is appearing that metabolic signals are integrated at the levels of first-order neurons equipped with the proper receptors, ant that these neurons send their signals towards hypothalamic GnRH neurons which constitute the integrative element of this network.

THCCOOH concentrations in whole blood: Are they useful in discriminating occasional from heavy smokers?

Some forensic and clinical circumstances require knowledge of the frequency of drug use. Care of the patient, administrative, and legal consequences will be different if the subject is a regular or an occasional cannabis smoker. To this end, 11-nor-9-carboxy-Δ9-tetrahydrocannabinol (THCCOOH) has been proposed as a criterion to help to distinguish between these two groups of users. However, to date this indicator has not been adequately assessed under experimental conditions. We carried out a controlled administration study of smoked cannabis with a placebo. Cannabinoid levels were determined in whole blood using tandem mass spectrometry. Significantly high differences in THCCOOH concentrations were found between the two groups when measured during the screening visit, prior to the smoking session, and throughout the day of the experiment. Receiver operating characteristic (ROC) curves were determined and two threshold criteria were proposed in order to distinguish between these groups: a free THCCOOH concentration below 3 µg/L suggested an occasional consumption (≤ 1 joint/week) while a concentration higher than 40 µg/L corresponded to a heavy use (≥ 10 joints/month). These thresholds were tested and found to be consistent with previously published experimental data. The decision threshold of 40 µg/L could be a cut-off for possible disqualification for driving while under the influence of cannabis. A further medical assessment and follow-up would be necessary for the reissuing of a driving license once abstinence from cannabis has been demonstrated. A THCCOOH level below 3 µg/L would indicate that no medical assessment is required. Copyright © 2013 John Wiley & Sons, Ltd.

Probing the different chaperone activities of the bacterial HSP70-HSP40 system using a thermolabile luciferase substrate.

During mild heat-stress, a native thermolabile polypeptide may partially unfold and transiently expose water-avoiding hydrophobic segments that readily tend to associate into a stable misfolded species, rich in intra-molecular non-native beta-sheet structures. When the concentration of the heat-unfolded intermediates is elevated, the exposed hydrophobic segments tend to associate with other molecules into large stable insoluble complexes, also called "aggregates." In mammalian cells, stress- and mutation-induced protein misfolding and aggregation may cause degenerative diseases and aging. Young cells, however, effectively counteract toxic protein misfolding with a potent network of molecular chaperones that bind hydrophobic surfaces and actively unfold otherwise stable misfolded and aggregated polypeptides. Here, we followed the behavior of a purified, initially mostly native thermolabile luciferase mutant, in the presence or absence of the Escherichia coli DnaK-DnaJ-GrpE chaperones and/or of ATP, at 22 °C or under mild heat-stress. We concomitantly measured luciferase enzymatic activity, Thioflavin-T fluorescence, and light-scattering to assess the effects of temperature and chaperones on the formation, respectively, of native, unfolded, misfolded, and/or of aggregated species. During mild heat-denaturation, DnaK-DnaJ-GrpE+ATP best maintained, although transiently, high luciferase activity and best prevented heat-induced misfolding and aggregation. In contrast, the ATP-less DnaK and DnaJ did not maintain optimal luciferase activity and were less effective at preventing luciferase misfolding and aggregation. We present a model accounting for the experimental data, where native, unfolded, misfolded, and aggregated species spontaneously inter-convert, and in which DnaK-DnaJ-GrpE+ATP specifically convert stable misfolded species into unstable unfolded intermediates.

Analytical expressions for the flame front speed in the downward combustion of thin solid fuels and comparison to experiments

We derive analytical expressions for the propagation speed of downward combustion fronts of thin solid fuels with a background flow initially at rest. The classical combustion model for thin solid fuels that consists of five coupled reaction-convection-diffusion equations is here reduced into a single equation with the gas temperature as the single variable. For doing so we apply a two-zone combustion model that divides the system into a preheating region and a pyrolyzing region. The speed of the combustion front is obtained after matching the temperature and its derivative at the location that separates both regions.We also derive a simplified version of this analytical expression expected to be valid for a wide range of cases. Flame front velocities predicted by our analyticalexpressions agree well with experimental data found in the literature for a large variety of cases and substantially improve the results obtained from a previous well-known analytical expression

Reaction-diffusion wave fronts: Multigeneration biological species under climate change

A generalization of reaction-diffusion models to multigeneration biological species is presented. It is based on more complex random walks than those in previous approaches. The new model is developed analytically up to infinite order. Our predictions for the speed agree to experimental data for several butterfly species better than existing models. The predicted dependence for the speed on the number of generations per year allows us to explain the change in speed observed for a specific invasion

Metabolic Flux and Compartmentation Analysis in the Brain In vivo.

Through significant developments and progresses in the last two decades, in vivo localized nuclear magnetic resonance spectroscopy (MRS) became a method of choice to probe brain metabolic pathways in a non-invasive way. Beside the measurement of the total concentration of more than 20 metabolites, (1)H MRS can be used to quantify the dynamics of substrate transport across the blood-brain barrier by varying the plasma substrate level. On the other hand, (13)C MRS with the infusion of (13)C-enriched substrates enables the characterization of brain oxidative metabolism and neurotransmission by incorporation of (13)C in the different carbon positions of amino acid neurotransmitters. The quantitative determination of the biochemical reactions involved in these processes requires the use of appropriate metabolic models, whose level of details is strongly related to the amount of data accessible with in vivo MRS. In the present work, we present the different steps involved in the elaboration of a mathematical model of a given brain metabolic process and its application to the experimental data in order to extract quantitative brain metabolic rates. We review the recent advances in the localized measurement of brain glucose transport and compartmentalized brain energy metabolism, and how these reveal mechanistic details on glial support to glutamatergic and GABAergic neurons.

Ab initio benchmark study for the oxidative addition of CH4 to Pd: importance of basis-set flexibility and polarization

To obtain a state-of-the-art benchmark potential energy surface (PES) for the archetypal oxidative addition of the methane C-H bond to the palladium atom, we have explored this PES using a hierarchical series of ab initio methods (Hartree-Fock, second-order Møller-Plesset perturbation theory, fourth-order Møller-Plesset perturbation theory with single, double and quadruple excitations, coupled cluster theory with single and double excitations (CCSD), and with triple excitations treated perturbatively [CCSD(T)]) and hybrid density functional theory using the B3LYP functional, in combination with a hierarchical series of ten Gaussian-type basis sets, up to g polarization. Relativistic effects are taken into account either through a relativistic effective core potential for palladium or through a full four-component all-electron approach. Counterpoise corrected relative energies of stationary points are converged to within 0.1-0.2 kcal/mol as a function of the basis-set size. Our best estimate of kinetic and thermodynamic parameters is -8.1 (-8.3) kcal/mol for the formation of the reactant complex, 5.8 (3.1) kcal/mol for the activation energy relative to the separate reactants, and 0.8 (-1.2) kcal/mol for the reaction energy (zero-point vibrational energy-corrected values in parentheses). This agrees well with available experimental data. Our work highlights the importance of sufficient higher angular momentum polarization functions, f and g, for correctly describing metal-d-electron correlation and, thus, for obtaining reliable relative energies. We show that standard basis sets, such as LANL2DZ+ 1f for palladium, are not sufficiently polarized for this purpose and lead to erroneous CCSD(T) results. B3LYP is associated with smaller basis set superposition errors and shows faster convergence with basis-set size but yields relative energies (in particular, a reaction barrier) that are ca. 3.5 kcal/mol higher than the corresponding CCSD(T) values

Redox dysregulation in the pathophysiology of schizophrenia and bipolar disorder: insights from animal models.

Abstract Significance: Schizophrenia (SZ) and bipolar disorder (BD) are classified as two distinct diseases. However, accumulating evidence shows that both disorders share genetic, pathological, and epidemiological characteristics. Based on genetic and functional findings, redox dysregulation due to an imbalance between pro-oxidants and antioxidant defense mechanisms has been proposed as a risk factor contributing to their pathophysiology. Recent Advances: Altered antioxidant systems and signs of increased oxidative stress are observed in peripheral tissues and brains of SZ and BD patients, including abnormal prefrontal levels of glutathione (GSH), the major cellular redox regulator and antioxidant. Here we review experimental data from rodent models demonstrating that permanent as well as transient GSH deficit results in behavioral, morphological, electrophysiological, and neurochemical alterations analogous to pathologies observed in patients. Mice with GSH deficit display increased stress reactivity, altered social behavior, impaired prepulse inhibition, and exaggerated locomotor responses to psychostimulant injection. These behavioral changes are accompanied by N-methyl-D-aspartate receptor hypofunction, elevated glutamate levels, impairment of parvalbumin GABA interneurons, abnormal neuronal synchronization, altered dopamine neurotransmission, and deficient myelination. Critical Issues: Treatment with the GSH precursor and antioxidant N-acetylcysteine normalizes some of those deficits in mice, but also improves SZ and BD symptoms when given as adjunct to antipsychotic medication. Future Directions: These data demonstrate the usefulness of GSH-deficient rodent models to identify the mechanisms by which a redox imbalance could contribute to the development of SZ and BD pathophysiologies, and to develop novel therapeutic approaches based on antioxidant and redox regulator compounds. Antioxid. Redox Signal. 18, 1428-1443.