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.

MOLECULAR REORGANIZATION FOLLOWING SENSORY STIMULATION IN THE ADULT MOUSE BARREL CORTEX

Sensory information is an important factor in shaping neuronal circuits during development and adulthood. In the barrel cortex of adult rodents, cells from layer IV are able to adapt their functional state to an increased flow of sensory information from the mystacial whisker follicles. Previous studies in our group have shown that whisker stimulation induces the formation of inhibitory synapses in the corresponding barrel (Knott et al., 2002) and decreases neuronal responses toward the deflection of the stimulated whisker (Quairiaux et al., 2007). Together these observations have turned the barrel cortex into a model to study homeostatic plasticity. At the cellular level, neuronal activity triggers intracellular signaling cascades leading to a transcriptional response. To further characterize the molecular pathways involved in the synaptic changes after whisker stimulation in the adult mouse, a previous doctoral student in our group performed a microarray analysis on laser-dissected barrels in sections through layer IV. This study identified the regulation (up and down) of a series of genes in the stimulated barrels (thesis of Johnston-Wenger, 2010). We here focused on ten genes that presented the highest fold change according to the microarray analysis. Out of these genes, 7 are known as neuronal activity-dependent genes (Tnncl, Nptx2, Sorcs3, Ptgs2, Nr4a2, Npas4 and Adcyapl) whereas three have so far not been related to neuronal plasticity (Scn7a, Pcdhl5 and Cede3). The study aimed at confirming the results of the microarray analysis and localizing molecular modifications in the stimulated barrel column at the cellular level. In situ hybridization for Pcdhl5 after different periods of whisker stimulation (3, 6, 9, 15, 24 hrs) allowed us to confirm that the 1.25 fold change used for the microarray analysis is an appropriate threshold for considering a regulation significant after sensory-stimulation. Moreover, we confirmed with in situ hybridization a significant upregulation of the genes of interest in the stimulated barrels. In situ hybridization and immunohistochemistry allowed us to observe the distribution of the genes of interest and the corresponding protein products at the cellular level. Three observations were made: 1) alterations of the expression was restricted to the stimulated barrels for all genes tested; 2) within a barrel column not all cells responded to whisker stimulation with an altered gene expression; 3) in the stimulated barrels, two different patterns of mRNA and protein expression can be distinguished. We hypothesize that this segregation of the activity-induced gene expression reflects the segregation of the two principal thalamocortical pathways conveying the sensory information to the barrel cortex. Moreover, only neurons reaching the critical threshold will modify their gene expression program resulting in structural as well as physiological modifications that prevent the subsequent propagation of the excess of excitation to the postsynaptic targets. The activity-induced gene expression is therefore adapted in a cell-type-specific manner to induce a homeostatic response to the entire neuronal network involved in the integration of the sensory information. This to our knowledge the first study showing the distinct, but complementary contribution of the two thalamocortical pathways in experience-dependent plasticity in the adult mouse barrel cortex. -- L'information sensorielle nous permet de continuellement façonner nos circuits neuronaux autant durant le développement qu'à l'âge adulte. Chez le rongeur l'information sensorielle perçue par les vibrisses est intégrée au niveau du cortex somatosensoriel primaire (appelé en anglais « barrel cortex ») dont les cellules de la couche IV sont capables d'adapter leur état fonctionnel en réponse à une augmentation d'activité neuronale. Ce modèle expérimental a permis à notre groupe de recherche d'observer des changements rapides du circuit neuronal en fonction de l'activité sensorielle. En effet, la stimulation continue d'une vibrisse d'une souris adulte pendant 24 heures induit non seulement un remaniement synaptique (Knott et al., 2002), mais également des changements physiologiques au niveau des neurones du tonneau correspondant (Quairiaux et al., 2007). Ces observations nous permettent d'affirmer que le « barrel cortex » est un modèle approprié pour y étudier la plasticité synaptique. Au niveau cellulaire, l'activité neuronale déclenche des cascades de signalisation intracellulaire résultant en une réponse transcriptionnelle. Afin de caractériser les voies moléculaires impliquées dans la plasticité synaptique, une puce à ARN nous a permis de comparer l'expression de gènes entre un tonneau correspondant à une vibrisse stimulée et un tonneau d'une vibrisse non-stimulée (Nathalie). Cette analyse a révélé un certain nombre de gènes régulés de manière positive ou négative par l'augmentation de l'activité neuronale. Nous nous sommes concentrés sur 10 gènes dont l'expression est fortement régulée. L'expression de sept d'entre eux a déjà été démontrée comme dépendante de l'activité neuronale (Tnncl, Nptx2, Sorcs3, Ptgs2, Nr4a2, Npas4 otAdcyapl) alors que l'expression des trois autres (Scn7a, Pcdhl5 et Cedei) n'a pour le moment pas encore été liée à la plasticité neuronale. Le but de cette thèse est de confirmer les résultats de la puce à ARN et de déterminer dans quel type cellulaire ces gènes sont exprimés. L'hybridation in situ pour le gène Pcdhl5, après différentes périodes de stimulation des vibrisses (3, 6, 9, 15 et 24 heures), nous a permis de confirmer que le seuil de 1.25x utilisé dans l'analyse de la puce à ARN est approprié pour considérer qu'un gène est régulé de manière significative par la stimulation sensorielle. Nous avons également pu confirmer à l'aide de cette technique que la stimulation sensorielle augmente significativement l'expression de ces dix gènes. L'expression de ces gènes au niveau cellulaire a été observée à l'aide des techniques d'hybridation in situ et d'immunohistochimie. Trois observations ont été faites : 1) la régulation de ces gènes est restreinte aux tonneaux correspondants aux vibrisses stimulées ; 2) au niveau d'une colonne corticale correspondant aux vibrisses stimulées, seules certaines cellules présentent une altération de leur expression génique ; 3) au niveau des tonneaux stimulés, deux profils d'expression d'ARNm et de protéines sont observés. Notre hypothèse est que cette distribution pourrait correspondre à la terminaison ségrégée des deux voies thalamocortical qui amènent l'information sensorielle dans le cortex cérébral. De plus, seul les neurones atteignant le seuil critique d'activation modifient leur expression génique en réponse à la stimulation sensorielle. Ces changements d'expression géniques vont permettre à la cellule de modifier ses propriétés structurales et physiologiques de manière a prevenir la propagation d'un excès d'activité neuronale au niveau de ses cibles postsynaptics. L'activité neuronale agit donc spécifiquement sur certains types cellulaires de maniere a induire une réponse homéostatique au niveau du réseau neuronal impliqué dans l'integration de l'information sensorielle. Nos travaux démontrent pour une première fois que les deux voies sensorielles contribuent d'une manière distincte et complémentaire à la plasticité corticale induite par un changement de l'activité sensorielle chez la souris adulte.

Biomarkers of human gastrointestinal tract regions.

Dysregulation of intestinal epithelial cell performance is associated with an array of pathologies whose onset mechanisms are incompletely understood. While whole-genomics approaches have been valuable for studying the molecular basis of several intestinal diseases, a thorough analysis of gene expression along the healthy gastrointestinal tract is still lacking. The aim of this study was to map gene expression in gastrointestinal regions of healthy human adults and to implement a procedure for microarray data analysis that would allow its use as a reference when screening for pathological deviations. We analyzed the gene expression signature of antrum, duodenum, jejunum, ileum, and transverse colon biopsies using a biostatistical method based on a multivariate and univariate approach to identify region-selective genes. One hundred sixty-six genes were found responsible for distinguishing the five regions considered. Nineteen had never been described in the GI tract, including a semaphorin probably implicated in pathogen invasion and six novel genes. Moreover, by crossing these genes with those retrieved from an existing data set of gene expression in the intestine of ulcerative colitis and Crohn's disease patients, we identified genes that might be biomarkers of Crohn's and/or ulcerative colitis in ileum and/or colon. These include CLCA4 and SLC26A2, both implicated in ion transport. This study furnishes the first map of gene expression along the healthy human gastrointestinal tract. Furthermore, the approach implemented here, and validated by retrieving known gene profiles, allowed the identification of promising new leads in both healthy and disease states.

Schistosome vaccines: a critical appraisal

An effective schistosome vaccine is a desirable control tool but progress towards that goal has been slow. Protective immunity has been difficult to demonstrate in humans, particularly children, so no routes to a vaccine have emerged from that source. The concept of concomitant immunity appeared to offer a paradigm for a vaccine operating against incoming larvae in the skin but did not yield the expected dividends. The mining of crude parasite extracts, the use of monoclonal antibodies and protein selection based on immunogenicity produced a panel of vaccine candidates, mostly of cytoplasmic origin. However, none of these performed well in independent rodent trials, but glutathione-S-transferease from Schistosoma haematobium is currently undergoing clinical trials as an anti-fecundity vaccine. The sequencing of the S. mansoni transcriptome and genome and the development of proteomic and microarray technologies has dramatically improved the possibilities for identifying novel vaccine candidates, particularly proteins secreted from or exposed at the surface of schistosomula and adult worms. These discoveries are leading to a new round of protein expression and protection experiments that will enable us to evaluate systematically all the major targets available for immune intervention. Only then will we know if schistosomes have an Achilles' heel.

Notch signaling is required for normal prostatic epithelial cell proliferation and differentiation.

Notch pathway is crucial for stem/progenitor cell maintenance, growth and differentiation in a variety of tissues. Using a transgenic cell ablation approach, we found in our previous study that cells expressing Notch1 are crucial for prostate early development and re-growth. Here, we further define the role of Notch signaling in regulating prostatic epithelial cell growth and differentiation using biochemical and genetic approaches in ex vivo or in vivo systems. Treatment of developing prostate grown in culture with inhibitors of gamma-secretase/presenilin, which is required for Notch cleavage and activation, caused a robust increase in proliferation of epithelial cells co-expressing cytokeratin 8 and 14, lack of luminal/basal layer segregation and dramatically reduced branching morphogenesis. Using conditional Notch1 gene deletion mouse models, we found that inactivation of Notch1 signaling resulted in profound prostatic alterations, including increased tufting, bridging and enhanced epithelial proliferation. Cells within these lesions co-expressed both luminal and basal cell markers, a feature of prostatic epithelial cells in predifferentiation developmental stages. Microarray analysis revealed that the gene expression in a number of genetic networks was altered following Notch1 gene deletion in prostate. Furthermore, expression of Notch1 and its effector Hey-1 gene in human prostate adenocarcinomas were found significantly down-regulated compared to normal control tissues. Taken together, these data suggest that Notch signaling is critical for normal cell proliferation and differentiation in the prostate, and deregulation of this pathway may facilitate prostatic tumorigenesis.

Sox10 promotes the formation and maintenance of giant congenital naevi and melanoma.

Giant congenital naevi are pigmented childhood lesions that frequently lead to melanoma, the most aggressive skin cancer. The mechanisms underlying this malignancy are largely unknown, and there are no effective therapies. Here we describe a mouse model for giant congenital naevi and show that naevi and melanoma prominently express Sox10, a transcription factor crucial for the formation of melanocytes from the neural crest. Strikingly, Sox10 haploinsufficiency counteracts Nras(Q61K)-driven congenital naevus and melanoma formation without affecting the physiological functions of neural crest derivatives in the skin. Moreover, Sox10 is also crucial for the maintenance of neoplastic cells in vivo. In human patients, virtually all congenital naevi and melanomas are SOX10 positive. Furthermore, SOX10 silencing in human melanoma cells suppresses neural crest stem cell properties, counteracts proliferation and cell survival, and completely abolishes in vivo tumour formation. Thus, SOX10 represents a promising target for the treatment of congenital naevi and melanoma in human patients.

Changes in the transcriptional profile of transporters in the intestine along the anterior-posterior and crypt-villus axes.

BACKGROUND: The purpose of this work was to characterize the expression of drug and nutrient carriers along the anterior-posterior and crypt-villus axes of the intestinal epithelium and to study the validity of utilizing whole gut tissue rather than purified epithelial cells to examine regional variations in gene expression. RESULTS: We have characterized the mRNA expression profiles of 76 % of all currently known transporters along the anterior-posterior axis of the gut. This is the first study to describe the expression profiles of the majority of all known transporters in the intestine. The expression profiles of transporters, as defined according to the Gene Ontology consortium, were measured in whole tissue of the murine duodenum, jejunum, ileum and colon using high-density microarrays. For nine transporters (Abca1, Abcc1, Abcc3, Abcg8, Slc10a2, Slc28a2, Slc2a1, Slc34a2 and Slc5a8), the mRNA profiles were further measured by RT-PCR in laser micro-dissected crypt and villus epithelial cells corresponding to the aforementioned intestinal regions. With respect to differentially regulated transporters, the colon had a distinct expression profile from small intestinal segments. The majority (59 % for p cutoff < or = 0.05) of transporter mRNA levels were constant across the intestinal sections studied. For the transporter subclass "carrier activity", which contains the majority of known carriers for biologically active compounds, a significant change (p < or = 0.05) along the anterior-posterior axis was observed. CONCLUSION: All nine transporters examined in laser-dissected material demonstrated good replication of the region-specific profiles revealed by microarray. Furthermore, we suggest that the distribution characteristics of Slc5a8 along the intestinal tract render it a suitable candidate carrier for monocarboxylate drugs in the posterior portion of the intestine. Our findings also predict that there is a significant difference in the absorption of carrier-mediated compounds in the different intestinal segments. The most pronounced differences can be expected between the adjoining segments ileum and colon, but the differences between the other adjoining segments are not negligible. Finally, for the examined genes, profiles measured in whole intestinal tissue extracts are representative of epithelial cell-only gene expression.

“Unmixing” Tissue Gene Expression Signatures from Tumor Biopsies

Emergent molecular measurement methods, such as DNA microarray, qRTPCR, andmany others, offer tremendous promise for the personalized treatment of cancer. Thesetechnologies measure the amount of specific proteins, RNA, DNA or other moleculartargets from tumor specimens with the goal of “fingerprinting” individual cancers. Tumorspecimens are heterogeneous; an individual specimen typically contains unknownamounts of multiple tissues types. Thus, the measured molecular concentrations resultfrom an unknown mixture of tissue types, and must be normalized to account for thecomposition of the mixture.For example, a breast tumor biopsy may contain normal, dysplastic and cancerousepithelial cells, as well as stromal components (fatty and connective tissue) and bloodand lymphatic vessels. Our diagnostic interest focuses solely on the dysplastic andcancerous epithelial cells. The remaining tissue components serve to “contaminate”the signal of interest. The proportion of each of the tissue components changes asa function of patient characteristics (e.g., age), and varies spatially across the tumorregion. Because each of the tissue components produces a different molecular signature,and the amount of each tissue type is specimen dependent, we must estimate the tissuecomposition of the specimen, and adjust the molecular signal for this composition.Using the idea of a chemical mass balance, we consider the total measured concentrationsto be a weighted sum of the individual tissue signatures, where weightsare determined by the relative amounts of the different tissue types. We develop acompositional source apportionment model to estimate the relative amounts of tissuecomponents in a tumor specimen. We then use these estimates to infer the tissuespecificconcentrations of key molecular targets for sub-typing individual tumors. Weanticipate these specific measurements will greatly improve our ability to discriminatebetween different classes of tumors, and allow more precise matching of each patient tothe appropriate treatment

Mutation screening of multiple genes in Spanish patients with autosomal recessive retinitis pigmentosa by targeted resequencing.

Retinitis Pigmentosa (RP) is a heterogeneous group of inherited retinal dystrophies characterised ultimately by the loss of photoreceptor cells. RP is the leading cause of visual loss in individuals younger than 60 years, with a prevalence of about 1 in 4000. The molecular genetic diagnosis of autosomal recessive RP (arRP) is challenging due to the large genetic and clinical heterogeneity. Traditional methods for sequencing arRP genes are often laborious and not easily available and a screening technique that enables the rapid detection of the genetic cause would be very helpful in the clinical practice. The goal of this study was to develop and apply microarray-based resequencing technology capable of detecting both known and novel mutations on a single high-throughput platform. Hence, the coding regions and exon/intron boundaries of 16 arRP genes were resequenced using microarrays in 102 Spanish patients with clinical diagnosis of arRP. All the detected variations were confirmed by direct sequencing and potential pathogenicity was assessed by functional predictions and frequency in controls. For validation purposes 4 positive controls for variants consisting of previously identified changes were hybridized on the array. As a result of the screening, we detected 44 variants, of which 15 are very likely pathogenic detected in 14 arRP families (14%). Finally, the design of this array can easily be transformed in an equivalent diagnostic system based on targeted enrichment followed by next generation sequencing.