Bioinformática: interfaz web para estudiar el efecto de diferentes condiciones sobre la expresión de los genes

El trabajo realizado se divide en dos bloques bien diferenciados, ambos relacionados con el análisis de microarrays. El primer bloque consiste en agrupar las condiciones muestrales de todos los genes en grupos o clústers. Estas agrupaciones se obtienen al aplicar directamente sobre la microarray los siguientes algoritmos de agrupación: SOM,PAM,SOTA,HC y al aplicar sobre la microarray escalada con PC y MDS los siguientes algoritmos: SOM,PAM,SOTA,HC y K-MEANS. El segundo bloque consiste en realizar una búsqueda de genes basada en los intervalos de confianza de cada clúster de la agrupación activa. Las condiciones de búsqueda ajustadas por el usuario se validan para cada clúster respecto el valor basal 0 y respecto el resto de clústers, para estas validaciones se usan los intervalos de confianza. Estos dos bloques se integran en una aplicación web ya existente, el applet PCOPGene, alojada en el servidor: http://revolutionresearch.uab.es.

Selective and powerful stress gene expression in Arabidopsis in response to malondialdehyde.

The provenance, half-life and biological activity of malondialdehyde (MDA) were investigated in Arabidopsis thaliana. We provide genetic confirmation of the hypothesis that MDA originates from fatty acids containing more than two methylene-linked double bonds, showing that tri-unsaturated fatty acids are the in vivo source of up to 75% of MDA. The abundance of the combined pool of free and reversibly bound MDA did not change dramatically in stress, although a significant increase in the free MDA pool under oxidative conditions was observed. The half-life of infiltrated MDA indicated rapid metabolic turnover/sequestration. Exposure of plants to low levels of MDA using a recently developed protocol powerfully upregulated many genes on a cDNA microarray with a bias towards those implicated in abiotic/environmental stress (e.g. ROF1 and XERO2). Remarkably, and in contrast to the activities of other reactive electrophile species (i.e. small vinyl ketones), none of the pathogenesis-related (PR) genes tested responded to MDA. The use of structural mimics of MDA isomers suggested that the propensity of the molecule to act as a cross-linking/modifying reagent might contribute to the activation of gene expression. Changes in the concentration/localisation of unbound MDA in vivo could strongly affect stress-related transcription.

Análisis de la expresión de las moléculas relacionadas con el complejo principal de histocompatibilidad de clase I y de NKG2D como respuesta al daño oxidativo del DNA en el cáncer urotelial de vejiga. Relación con las características clínico patológicas y supervivencia

Es tracta d'un estudi retrospectiu de casos de 280 pacients diagnosticats de tumor vesical primari amb un seguiment mínim de 8 anys. S'ha construït un Tissue microarray i mitjançant mètodes semiquantitatius d’inmunohistoquímica es determinarà l'expressió de les molècules MICA (MHC class I chain-related gene A) i del seu receptor NKG2D (Natural-Killer group 2-member D) a nivell tissular, relacionant-lo amb variables anatomopatològiques segons els grups de risc, hàbit tabàquic i sexe. Finalment valorarem l'expressió de MICA/NKG2D com a factor independent de recidiva / progressió tumoral. En la literatura només existeixen 2 treballs que relacionin MICA amb el càncer vesical.

Gene expression mapping in neuropathic pain : molecular plasticity in nociceptors and superficial dorsal horn

RESUME : La douleur neuropathique est le résultat d'une lésion ou d'un dysfonctionnement du système nerveux. Les symptômes qui suivent la douleur neuropathique sont sévères et leur traitement inefficace. Une meilleure approche thérapeutique peut être proposée en se basant sur les mécanismes pathologiques de la douleur neuropathique. Lors d'une lésion périphérique une douleur neuropathique peut se développer et affecter le territoire des nerfs lésés mais aussi les territoires adjacents des nerfs non-lésés. Une hyperexcitabilité des neurones apparaît au niveau des ganglions spinaux (DRG) et de la corne dorsale (DH) de la moelle épinière. Le but de ce travail consiste à mettre en évidence les modifications moléculaires associées aux nocicepteurs lésés et non-lésés au niveau des DRG et des laminae I et II de la corne dorsale, là où l'information nociceptive est intégrée. Pour étudier les changements moléculaires liés à la douleur neuropathique nous utilisons le modèle animal d'épargne du nerf sural (spared nerve injury model, SNI) une semaine après la lésion. Pour la sélection du tissu d'intérêt nous avons employé la technique de la microdissection au laser, afin de sélectionner une sous-population spécifique de cellules (notamment les nocicepteurs lésés ou non-lésés) mais également de prélever le tissu correspondant dans les laminae superficielles. Ce travail est couplé à l'analyse à large spectre du transcriptome par puce ADN (microarray). Par ailleurs, nous avons étudié les courants électriques et les propriétés biophysiques des canaux sodiques (Na,,ls) dans les neurones lésés et non-lésés des DRG. Aussi bien dans le système nerveux périphérique, entre les neurones lésés et non-lésés, qu'au niveau central avec les aires recevant les projections des nocicepteurs lésés ou non-lésés, l'analyse du transcriptome montre des différences de profil d'expression. En effet, nous avons constaté des changements transcriptionnels importants dans les nocicepteurs lésés (1561 gènes, > 1.5x et pairwise comparaison > 77%) ainsi que dans les laminae correspondantes (618 gènes), alors que ces modifications transcriptionelles sont mineures au niveau des nocicepteurs non-lésés (60 gènes), mais important dans leurs laminae de projection (459 gènes). Au niveau des nocicepteurs, en utilisant la classification par groupes fonctionnels (Gene Ontology), nous avons observé que plusieurs processus biologiques sont modifiés. Ainsi des fonctions telles que la traduction des signaux cellulaires, l'organisation du cytosquelette ainsi que les mécanismes de réponse au stress sont affectés. Par contre dans les neurones non-lésés seuls les processus biologiques liés au métabolisme et au développement sont modifiés. Au niveau de la corne dorsale de la moelle, nous avons observé des modifications importantes des processus immuno-inflammatoires dans l'aire affectée par les nerfs lésés et des changements associés à l'organisation et la transmission synaptique au niveau de l'aire des nerfs non-lésés. L'analyse approfondie des canaux sodiques a démontré plusieurs changements d'expression, principalement dans les neurones lésés. Les analyses fonctionnelles n'indiquent aucune différence entre les densités de courant tétrodotoxine-sensible (TTX-S) dans les neurones lésés et non-lésés même si les niveaux d'expression des ARNm des sous-unités TTX-S sont modifiés dans les neurones lésés. L'inactivation basale dépendante du voltage des canaux tétrodotoxine-insensible (TTX-R) est déplacée vers des potentiels positifs dans les cellules lésées et non-lésées. En revanche la vitesse de récupération des courants TTX-S et TTX-R après inactivation est accélérée dans les neurones lésés. Ces changements pourraient être à l'origine de l'altération de l'activité électrique des neurones sensoriels dans le contexte des douleurs neuropathiques. En résumé, ces résultats suggèrent l'existence de mécanismes différenciés affectant les neurones lésés et les neurones adjacents non-lésés lors de la mise en place la douleur neuropathique. De plus, les changements centraux au niveau de la moelle épinière qui surviennent après lésion sont probablement intégrés différemment selon la perception de signaux des neurones périphériques lésés ou non-lésés. En conclusion, ces modulations complexes et distinctes sont probablement des acteurs essentiels impliqués dans la genèse et la persistance des douleurs neuropathiques. ABSTRACT : Neuropathic pain (NP) results from damage or dysfunction of the peripheral or central nervous system. Symptoms associated with NP are severe and difficult to treat. Targeting NP mechanisms and their translation into symptoms may offer a better therapeutic approach.Hyperexcitability of the peripheral and central nervous system occurs in the dorsal root ganglia (DRG) and the dorsal horn (DH) of the spinal cord. We aimed to identify transcriptional variations in injured and in adjacent non-injured nociceptors as well as in corresponding laminae I and II of DH receiving their inputs.We investigated changes one week after the injury induced by the spared nerve injury model of NP. We employed the laser capture microdissection (LCM) for the procurement of specific cell-types (enrichment in nociceptors of injured/non-injured neurons) and laminae in combination with transcriptional analysis by microarray. In addition, we studied functionál properties and currents of sodium channels (Nav1s) in injured and neighboring non-injured DRG neurons.Microarray analysis at the periphery between injured and non-injured DRG neurons and centrally between the area of central projections from injured and non-injured neurons show significant and differential expression patterns. We reported changes in injured nociceptors (1561 genes, > 1.5 fold, >77% pairwise comparison) and in corresponding DH laminae (618 genes), while less modifications occurred in non-injured nociceptors (60 genes) and in corresponding DH laminae (459 genes). At the periphery, we observed by Gene Ontology the involvement of multiple biological processes in injured neurons such as signal transduction, cytoskeleton organization or stress responses. On contrast, functional overrepresentations in non-injured neurons were noted only in metabolic or developmentally related mechanisms. At the level of superficial laminae of the dorsal horn, we reported changes of immune and inflammatory processes in injured-related DH and changes associated with synaptic organization and transmission in DH corresponding to non-injured neurons. Further transcriptional analysis of Nav1s indicated several changes in injured neurons. Functional analyses of Nav1s have established no difference in tetrodotoxin-sensitive (TTX-S) current densities in both injured and non-injured neurons, despite changes in TTX-S Nav1s subunit mRNA levels. The tetrodotoxin-resistant (TTX-R) voltage dependence of steady state inactivation was shifted to more positive potentials in both injured and non-injured neurons, and the rate of recovery from inactivation of TTX-S and TTX-R currents was accelerated in injured neurons. These changes may lead to alterations in neuronal electrogenesis. Taken together, these findings suggest different mechanisms occurring in the injured neurons and the adjacent non-injured ones. Moreover, central changes after injury are probably driven in a different manner if they receive inputs from injured or non-injured neurons. Together, these distinct and complex modulations may contribute to NP.

Cellular variability of RpoS expression underlies subpopulation activation of an integrative and conjugative element.

Conjugative transfer of the integrative and conjugative element ICEclc in the bacterium Pseudomonas knackmussii is the consequence of a bistable decision taken in some 3% of cells in a population during stationary phase. Here we study the possible control exerted by the stationary phase sigma factor RpoS on the bistability decision. The gene for RpoS in P. knackmussii B13 was characterized, and a loss-of-function mutant was produced and complemented. We found that, in absence of RpoS, ICEclc transfer rates and activation of two key ICEclc promoters (P(int) and P(inR)) decrease significantly in cells during stationary phase. Microarray and gene reporter analysis indicated that the most direct effect of RpoS is on P(inR), whereas one of the gene products from the P(inR)-controlled operon (InrR) transmits activation to P(int) and other ICEclc core genes. Addition of a second rpoS copy under control of its native promoter resulted in an increase of the proportion of cells expressing the P(int) and P(inR) promoters to 18%. Strains in which rpoS was replaced by an rpoS-mcherry fusion showed high mCherry fluorescence of individual cells that had activated P(int) and P(inR), whereas a double-copy rpoS-mcherry-containing strain displayed twice as much mCherry fluorescence. This suggested that high RpoS levels are a prerequisite for an individual cell to activate P(inR) and thus ICEclc transfer. Double promoter-reporter fusions confirmed that expression of P(inR) is dominated by extrinsic noise, such as being the result of cellular variability in RpoS. In contrast, expression from P(int) is dominated by intrinsic noise, indicating it is specific to the ICEclc transmission cascade. Our results demonstrate how stochastic noise levels of global transcription factors can be transduced to a precise signaling cascade in a subpopulation of cells leading to ICE activation.

MIMAS 3.0 is a Multiomics Information Management and Annotation System.

BACKGROUND: DNA sequence integrity, mRNA concentrations and protein-DNA interactions have been subject to genome-wide analyses based on microarrays with ever increasing efficiency and reliability over the past fifteen years. However, very recently novel technologies for Ultra High-Throughput DNA Sequencing (UHTS) have been harnessed to study these phenomena with unprecedented precision. As a consequence, the extensive bioinformatics environment available for array data management, analysis, interpretation and publication must be extended to include these novel sequencing data types. DESCRIPTION: MIMAS was originally conceived as a simple, convenient and local Microarray Information Management and Annotation System focused on GeneChips for expression profiling studies. MIMAS 3.0 enables users to manage data from high-density oligonucleotide SNP Chips, expression arrays (both 3'UTR and tiling) and promoter arrays, BeadArrays as well as UHTS data using MIAME-compliant standardized vocabulary. Importantly, researchers can export data in MAGE-TAB format and upload them to the EBI's ArrayExpress certified data repository using a one-step procedure. CONCLUSION: We have vastly extended the capability of the system such that it processes the data output of six types of GeneChips (Affymetrix), two different BeadArrays for mRNA and miRNA (Illumina) and the Genome Analyzer (a popular Ultra-High Throughput DNA Sequencer, Illumina), without compromising on its flexibility and user-friendliness. MIMAS, appropriately renamed into Multiomics Information Management and Annotation System, is currently used by scientists working in approximately 50 academic laboratories and genomics platforms in Switzerland and France. MIMAS 3.0 is freely available via http://multiomics.sourceforge.net/.

Interfaz web para mostrar los genes según las relaciones que mantienen sus expresiones

La investigació entre les relacions dels nivells d’expressió dels gens aporta molta informació sobre els processos biològics i patològics. Mitjançant la tècnica de les microarrays es possibilita la investigació de les relacions d’expressió de milers de gens a la vegada. La finalitat d’aquest projecte es fent ús de l’aplicatiu web PCOPGene-Net, permetre la identificació dels gens per les relacions d’expressió no lineals que tenen amb la resta de gens i permetre també la identificació de les relacions d’expressió no lineals entre els gens d’una microarray.

Comparación de datos de expresión génica del servidor local con datos de una base de datos remota

En la presente memoria se detalla con precisión las diversas fases del trabajo para construir una aplicación web en el servidor http://revolutionresearch.uab.es que permite enriquecer los clusters de la microarray del usuario con información biomédica de una base de datos remota. Los clusters de origen estadístico (o no) de la microarray del usuario se enriquecen a partir de cruzar sus genes marcadores con la base de datos de genes marcadores de microarrays (base de datos remota) con clusters basados en información biomédica. La base de datos de genes marcadores de microarrays ha sido obtenida a partir de la base de datos de GEO Profiles del NCBI.

Gene expression databases and data mining.

The DNA microarray technology has arguably caught the attention of the worldwide life science community and is now systematically supporting major discoveries in many fields of study. The majority of the initial technical challenges of conducting experiments are being resolved, only to be replaced with new informatics hurdles, including statistical analysis, data visualization, interpretation, and storage. Two systems of databases, one containing expression data and one containing annotation data are quickly becoming essential knowledge repositories of the research community. This present paper surveys several databases, which are considered "pillars" of research and important nodes in the network. This paper focuses on a generalized workflow scheme typical for microarray experiments using two examples related to cancer research. The workflow is used to reference appropriate databases and tools for each step in the process of array experimentation. Additionally, benefits and drawbacks of current array databases are addressed, and suggestions are made for their improvement.

Messenger RNA expression of transporter and ion channel genes in undifferentiated and differentiated Caco-2 cells compared to human intestines.

PURPOSE: The purpose of this work was to study the influence of cell differentiation on the mRNA expression of transporters and channels in Caco-2 cells and to assess Caco-2 cells as a model for carrier-mediated drug transport in the intestines. METHOD: Gene mRNA expression was measured using a custom-designed microarray chip with 750 deoxyoligonucleotide probes (70mers). Each oligomer was printed four times on poly-lysine-coated glass slides. Expression profiles were expressed as ratio values between fluorescence intensities of Cy3 and Cy5 dye-labeled cDNA derived from poly(A) + RNA samples of Caco-2 cells and total RNA of human intestines. RESULTS: Significant differences in the mRNA expression profile of transporters and channels were observed upon differentiation of Caco-2 cells from 5 days to 2 weeks in culture, including changes for MAT8, S-protein, and Nramp2. Comparing Caco-2 cells of different passage number revealed few changes in mRNAs except for GLUT3, which was down-regulated 2.4-fold within 13 passage numbers. Caco-2 cells had a similar expression profile when either cultured in flasks or on filters but differed more strongly from human small and large intestine, regardless of the differentiation state of Caco-2 cells. Expression of several genes highly transcribed in small or large intestines differed fourfold or more in Caco-2 cells. CONCLUSIONS: Although Caco-2 cells have proven a suitable model for studying carrier-mediated transport in human intestines, the expression of specific transporter and ion channel genes may differ substantially.

Differential gene expression in response to mechanical wounding and insect feeding in Arabidopsis.

Wounding in multicellular eukaryotes results in marked changes in gene expression that contribute to tissue defense and repair. Using a cDNA microarray technique, we analyzed the timing, dynamics, and regulation of the expression of 150 genes in mechanically wounded leaves of Arabidopsis. Temporal accumulation of a group of transcripts was correlated with the appearance of oxylipin signals of the jasmonate family. Analysis of the coronatine-insensitive coi1-1 Arabidopsis mutant that is also insensitive to jasmonate allowed us to identify a large number of COI1-dependent and COI1-independent wound-inducible genes. Water stress was found to contribute to the regulation of an unexpectedly large fraction of these genes. Comparing the results of mechanical wounding with damage by feeding larvae of the cabbage butterfly (Pieris rapae) resulted in very different transcript profiles. One gene was specifically induced by insect feeding but not by wounding; moreover, there was a relative lack of water stress-induced gene expression during insect feeding. These results help reveal a feeding strategy of P. rapae that may minimize the activation of a subset of water stress-inducible, defense-related genes.

Modular analysis of gene expression data with R.

SUMMARY: Large sets of data, such as expression profiles from many samples, require analytic tools to reduce their complexity. The Iterative Signature Algorithm (ISA) is a biclustering algorithm. It was designed to decompose a large set of data into so-called 'modules'. In the context of gene expression data, these modules consist of subsets of genes that exhibit a coherent expression profile only over a subset of microarray experiments. Genes and arrays may be attributed to multiple modules and the level of required coherence can be varied resulting in different 'resolutions' of the modular mapping. In this short note, we introduce two BioConductor software packages written in GNU R: The isa2 package includes an optimized implementation of the ISA and the eisa package provides a convenient interface to run the ISA, visualize its output and put the biclusters into biological context. Potential users of these packages are all R and BioConductor users dealing with tabular (e.g. gene expression) data. AVAILABILITY: http://www.unil.ch/cbg/ISA CONTACT: sven.bergmann@unil.ch

Developmental and Anatomical Constraints on Vertebrate Genome Evolution

Pendant ma thèse de doctorat, j'ai utilisé des espèces modèles, comme la souris et le poisson-zèbre, pour étudier les facteurs qui affectent l'évolution des gènes et leur expression. Plus précisément, j'ai montré que l'anatomie et le développement sont des facteurs clés à prendre en compte, car ils influencent la vitesse d'évolution de la séquence des gènes, l'impact sur eux de mutations (i.e. la délétion du gène est-elle létale ?), et leur tendance à se dupliquer. Où et quand il est exprimé impose à un gène certaines contraintes ou au contraire lui donne des opportunités d'évoluer. J'ai pu comparer ces tendances aux modèles classiques d'évolution de la morphologie, que l'on pensait auparavant refléter directement les contraintes s'appliquant sur le génome. Nous avons montré que les contraintes entre ces deux niveaux d'organisation ne peuvent pas être transférées simplement : il n'y a pas de lien direct entre la conservation du génotype et celle de phénotypes comme la morphologie. Ce travail a été possible grâce au développement d'outils bioinformatiques. Notamment, j'ai travaillé sur le développement de la base de données Bgee, qui a pour but de comparer l'expression des gènes entre différentes espèces de manière automatique et à large échelle. Cela implique une formalisation de l'anatomie, du développement et de concepts liés à l'homologie grâce à l'utilisation d'ontologies. Une intégration cohérente de données d'expression hétérogènes (puces à ADN, marqueurs de séquence exprimée, hybridations in situ) a aussi été nécessaire. Cette base de données est mise à jour régulièrement et disponible librement. Elle devrait contribuer à étendre les possibilités de comparaison de l'expression des gènes entre espèces pour des études d'évo-devo (évolution du développement) et de génomique. During my PhD, I used model species of vertebrates, such as mouse and zebrafish, to study factors affecting the evolution of genes and their expression. More precisely I have shown that anatomy and development are key factors to take into account, influencing the rate of gene sequence evolution, the impact of mutations (i.e. is the deletion of a gene lethal?), and the propensity of a gene to duplicate. Where and when genes are expressed imposes constraints, or on the contrary leaves them some opportunity to evolve. We analyzed these patterns in relation to classical models of morphological evolution in vertebrates, which were previously thought to directly reflect constraints on the genomes. We showed that the patterns of evolution at these two levels of organization do not translate smoothly: there is no direct link between the conservation of genotype and phenotypes such as morphology. This work was made possible by the development of bioinformatics tools. Notably, I worked on the development of the database Bgee, which aims at comparing gene expression between different species in an automated and large-scale way. This involves the formalization of anatomy, development, and concepts related to homology, through the use of ontologies. A coherent integration of heterogeneous expression data (microarray, expressed sequence tags, in situ hybridizations) is also required. This database is regularly updated and freely available. It should contribute to extend the possibilities for comparison of gene expression between species in evo-devo and genomics studies.

Gout in the spotlight.

Understanding how uric acid crystals provoke inflammation is crucial to improving our management of acute gout. It is well known that urate crystals stimulate monocytes and macrophages to elaborate inflammatory cytokines, but the tissue response of the synovium is less well understood. Microarray analysis of mRNA expression by these lining cells may help to delineate the genes that are modulated. Employing a murine air-pouch model, a number of genes expressed by innate immune cells were found to be rapidly upregulated by monosodium urate crystals. These findings provide new research avenues to investigate the physiopathology of gouty inflammation, and may eventually lead to new therapeutic targets in acute gout.

Interpretation of array comparative genome hybridization data: a major challenge.

The advent and application of high-resolution array-based comparative genome hybridization (array CGH) has led to the detection of large numbers of copy number variants (CNVs) in patients with developmental delay and/or multiple congenital anomalies as well as in healthy individuals. The notion that CNVs are also abundantly present in the normal population challenges the interpretation of the clinical significance of detected CNVs in patients. In this review we will illustrate a general clinical workflow based on our own experience that can be used in routine diagnostics for the interpretation of CNVs.