Comparative modular analysis of gene expression in vertebrate organs.

ABSTRACT: BACKGROUND: The degree of conservation of gene expression between homologous organs largely remains an open question. Several recent studies reported some evidence in favor of such conservation. Most studies compute organs' similarity across all orthologous genes, whereas the expression level of many genes are not informative about organ specificity. RESULTS: Here, we use a modularization algorithm to overcome this limitation through the identification of inter-species co-modules of organs and genes. We identify such co-modules using mouse and human microarray expression data. They are functionally coherent both in terms of genes and of organs from both organisms. We show that a large proportion of genes belonging to the same co-module are orthologous between mouse and human. Moreover, their zebrafish orthologs also tend to be expressed in the corresponding homologous organs. Notable exceptions to the general pattern of conservation are the testis and the olfactory bulb. Interestingly, some co-modules consist of single organs, while others combine several functionally related organs. For instance, amygdala, cerebral cortex, hypothalamus and spinal cord form a clearly discernible unit of expression, both in mouse and human. CONCLUSIONS: Our study provides a new framework for comparative analysis which will be applicable also to other sets of large-scale phenotypic data collected across different species.

MAR elements regulate the probability of epigenetic switching between active and inactive gene expression.

Gene expression often cycles between active and inactive states in eukaryotes, yielding variable or noisy gene expression in the short-term, while slow epigenetic changes may lead to silencing or variegated expression. Understanding how cells control these effects will be of paramount importance to construct biological systems with predictable behaviours. Here we find that a human matrix attachment region (MAR) genetic element controls the stability and heritability of gene expression in cell populations. Mathematical modeling indicated that the MAR controls the probability of long-term transitions between active and inactive expression, thus reducing silencing effects and increasing the reactivation of silent genes. Single-cell short-terms assays revealed persistent expression and reduced expression noise in MAR-driven genes, while stochastic burst of expression occurred without this genetic element. The MAR thus confers a more deterministic behavior to an otherwise stochastic process, providing a means towards more reliable expression of engineered genetic systems.

Anaerobically controlled expression system derived from the arcDABC operon of Pseudomonas aeruginosa: application to lipase production.

The anaerobically inducible arcDABC operon encodes the enzymes of the arginine deiminase pathway in Pseudomonas aeruginosa. Upon induction, the arcAB mRNAs and proteins reach high intracellular levels, because of a strong anaerobically controlled promoter and mRNA processing in arcD, leading to stable downstream transcripts. We explored the usefulness of this system for the construction of expression vectors. The lacZ gene of Escherichia coli was expressed to the highest levels when fused close to the arc promoter. Insertion of lacZ further downstream into arcA or arcB did not stabilize the intrinsically unstable lacZ mRNA. On the contrary, lacZ mRNA appeared to be a vulnerable endonuclease target destabilizing arcAB mRNAs in the 5'-to-3' direction in P. aeruginosa. The native arc promoter was modified for optional expression in the -10 sequence and in the -40 region, which is a binding site for the anaerobic regulator ANR. In P. aeruginosa grown either anaerobically or with oxygen limitation in unshaken cultures, this promoter was stronger than the induced tac promoter. The P. aeruginosa lipAH genes, which encode extracellular lipase and lipase foldase, respectively, were fused directly to the modified arc promoter in an IncQ vector plasmid. Semianaerobic static cultures of P. aeruginosa PAO1 carrying this recombinant plasmid overproduced extracellular lipase 30-fold during stationary phase compared with the production by strain PAO1 without the plasmid. Severe oxygen limitation, in contrast, resulted in poor lipase productivity despite effective induction of the ANR-dependent promoter, suggesting that secretion of active lipase is blocked by the absence of oxygen. In conclusion, the modified arc promoter is useful for driving the expression of cloned genes in P. aeruginosa during oxygen-limited growth and stationary phase.

Expression analyses of three members of the AtPHO1 family reveal differential interactions between signaling pathways involved in phosphate deficiency and the responses to auxin, cytokinin, and abscisic acid.

The PHO1 protein is involved in loading inorganic phosphate (Pi) to the root xylem. Ten genes homologous to AtPHO1 are present in the Arabidopsis thaliana (L.) Heyn genome. From this gene family, transcript levels of only AtPHO1, AtPHO1;H1 and AtPHO1;H10 were increased by Pi-deficiency. While the up-regulation of AtPHO1;H1 and AtPHO1;H10 by Pi deficiency followed the same rapid kinetics and was dependent on the PHR1 transcription factor, phosphite only strongly suppressed the expression of AtPHO1;H1 and had a minor effect on AtPHO1;H10. Addition of sucrose was found to increase transcript levels of both AtPHO1 and AtPHO1;H1 in Pi-sufficient or Pi-deficient plants, but to suppress AtPHO1:H10 under the same conditions. Treatments of plants with auxin or cytokinin had contrasting effect depending on the gene and on the Pi status of the plants. Thus, while both hormones down-regulated expression of AtPHO1 independently of the plant Pi status, auxin and cytokinin up-regulated AtPHO1;H1 and AtPHO1;H10 expression in Pi-sufficient plants and down-regulated expression in Pi-deficient plants. Treatments with abscisic acid inhibited AtPHO1 and AtPHO1;H1 expression in both Pi-sufficient and Pi-deficient plants, but increased AtPHO1;H10 expression under the same conditions. The inhibition of expression by abscisic acid of AtPHO1 and AtPHO1;H1, and of the Pi-starvation responsive genes AtPHT1;1 and AtIPS1, was dependant on the ABI1 type 2C protein phosphatase. These results reveal that various levels of cross talk between the signal transduction pathways to Pi, sucrose and phytohormones are involved in the regulation of expression of the three AtPHO1 homologues.

Chemoendocrine compared with endocrine adjuvant therapies for node-negative breast cancer: predictive value of centrally reviewed expression of estrogen and progesterone receptors--International Breast Cancer Study Group.

PURPOSE: To centrally assess estrogen receptor (ER) and progesterone receptor (PgR) levels by immunohistochemistry and investigate their predictive value for benefit of chemo-endocrine compared with endocrine adjuvant therapy alone in two randomized clinical trials for node-negative breast cancer. PATIENTS AND METHODS: International Breast Cancer Study Group Trial VIII compared cyclophosphamide, methotrexate, and fluorouracil (CMF) chemotherapy for 6 cycles followed by endocrine therapy with goserelin with either modality alone in pre- and perimenopausal patients. Trial IX compared three cycles of CMF followed by tamoxifen for 5 years versus tamoxifen alone in postmenopausal patients. Central Pathology Office reviewed 883 (83%) of 1,063 patients on Trial VIII and 1,365 (82%) of 1,669 on Trial IX and determined ER and PgR by immunohistochemistry. Disease-free survival (DFS) was compared across the spectrum of expression of each receptor using the Subpopulation Treatment Effect Pattern Plot methodology. RESULTS: Both receptors displayed a bimodal distribution, with substantial proportions showing no staining (receptor absent) and most of the remainder showing a high percentage of stained cells. Chemo-endocrine therapy yielded DFS superior to endocrine therapy alone for patients with receptor-absent tumors, and in some cases also for those with low levels of receptor expression. Among patients with ER-expressing tumors, additional prediction of benefit was suggested in absent or low PgR in Trial VIII but not in Trial IX. CONCLUSION: Low levels of ER and PgR are predictive of the benefit of adding chemotherapy to endocrine therapy. Low PgR may add further prediction among pre- and perimenopausal but not postmenopausal patients whose tumors express ER.

Evolution of mRNA expression levels of autosomal and sex-linked genes in amniotes

In addition to differences in protein-coding gene sequences, changes in expression resulting from mutations in regulatory sequences have long been hypothesized to be responsible for phenotypic differences between species. However, unlike comparison of genome sequences, few studies, generally restricted to pairwise comparisons of closely related mammalian species, have assessed between-species differences at the transcriptome level. They reported that gene expression evolves at different rates in various organs and in a pattern that is overall consistent with neutral models of evolution. In the first part of my thesis, I investigated the evolution of gene expression in therian mammals (i.e.7 placental and marsupials), based on microarray data from human, mouse and the gray short-tailed opossum (Monodelphis domestica). In addition to autosomal genes, a special focus was given to the evolution of X-linked genes. The therian X chromosome was recently shown to be younger than previously thought and to harbor a specific gene content (e.g., genes involved in brain or reproductive functions) that is thought to have been shaped by specific sex-related evolutionary forces. Sex chromosomes derive from ordinary autosomes and their differentiation led to the degeneration of the Y chromosome (in mammals) or W chromosome (in birds). Consequently, X- or Z-linked genes differ in gene dose between males and females such that the heterogametic sex has half the X/Z gene dose compared to the ancestral state. To cope with this dosage imbalance, mammals have been reported to have evolved mechanisms of dosage compensation.¦In the first project, I could first show that transcriptomes evolve at different rates in different organs. Out of the five tissues I investigated, the testis is the most rapidly evolving organ at the gene expression level while the brain has the most conserved transcriptome. Second, my analyses revealed that mammalian gene expression evolution is compatible with a neutral model, where the rates of change in gene expression levels is linked to the efficiency of purifying selection in a given lineage, which, in turn, is determined by the long-term effective population size in that lineage. Thus, the rate of DNA sequence evolution, which could be expected to determine the rate of regulatory sequence change, does not seem to be a major determinant of the rate of gene expression evolution. Thus, most gene expression changes seem to be (slightly) deleterious. Finally, X-linked genes seem to have experienced elevated rates of gene expression change during the early stage of X evolution. To further investigate the evolution of mammalian gene expression, we generated an extensive RNA-Seq gene expression dataset for nine mammalian species and a bird. The analyses of this dataset confirmed the patterns previously observed with microarrays and helped to significantly deepen our view on gene expression evolution.¦In a specific project based on these data, I sought to assess in detail patterns of evolution of dosage compensation in amniotes. My analyses revealed the absence of male to female dosage compensation in monotremes and its presence in marsupials and, in addition, confirmed patterns previously described for placental mammals and birds. I then assessed the global level of expression of X/Z chromosomes and contrasted this with its ancestral gene expression levels estimated from orthologous autosomal genes in species with non-homologous sex chromosomes. This analysis revealed a lack of up-regulation for placental mammals, the level of expression of X-linked genes being proportional to gene dose. Interestingly, the ancestral gene expression level was at least partially restored in marsupials as well as in the heterogametic sex of monotremes and birds. Finally, I investigated alternative mechanisms of dosage compensation and found that gene duplication did not seem to be a widespread mechanism to restore the ancestral gene dose. However, I could show that placental mammals have preferentially down-regulated autosomal genes interacting with X-linked genes which underwent gene expression decrease, and thus identified a novel alternative mechanism of dosage compensation.

Submicroscopic deletion in patients with Williams-Beuren syndrome influences expression levels of the nonhemizygous flanking genes.

Genomic imbalance is a common cause of phenotypic abnormalities. We measured the relative expression level of genes that map within the microdeletion that causes Williams-Beuren syndrome and within its flanking regions. We found, unexpectedly, that not only hemizygous genes but also normal-copy neighboring genes show decreased relative levels of expression. Our results suggest that not only the aneuploid genes but also the flanking genes that map several megabases away from a genomic rearrangement should be considered possible contributors to the phenotypic variation in genomic disorders.

Expression hierarchy of T cell epitopes from melanoma differentiation antigens: unexpected high level presentation of tyrosinase-HLA-A2 Complexes revealed by peptide-specific, MHC-restricted, TCR-like antibodies.

Peptide Ags presented by class I MHC molecules on human melanomas and that are recognized by CD8(+) T cells are the subjects of many studies of antitumor immunity and represent attractive candidates for therapeutic approaches. However, no direct quantitative measurements exist to reveal their expression hierarchy on the cell surface. Using novel recombinant Abs which bind these Ags with a peptide-specific, MHC-restricted manner, we demonstrate a defined pattern of expression hierarchy of peptide-HLA-A2 complexes derived from three major differentiation Ags: gp100, Melan-A/Mart-1, and tyrosinase. Studying melanoma cell lines derived from multiple patients, we reveal a surprisingly high level of presentation of tyrosinase-derived complexes and moderate to very low expression of complexes derived from other Ags. No correlation between Ag presentation and mRNA expression was found; however, protein stability may play a major role. These results provide new insights into the characteristics of Ag presentation and are particularly important when such targets are being considered for immunotherapy. These results may shed new light on relationships between Ag presentation and immune response to cancer Ags.

Analysis of gene expression patterns in animals

During my PhD, my aim was to provide new tools to increase our capacity to analyse gene expression patterns, and to study on a large-scale basis the evolution of gene expression in animals. Gene expression patterns (when and where a gene is expressed) are a key feature in understanding gene function, notably in development. It appears clear now that the evolution of developmental processes and of phenotypes is shaped both by evolution at the coding sequence level, and at the gene expression level.Studying gene expression evolution in animals, with complex expression patterns over tissues and developmental time, is still challenging. No tools are available to routinely compare expression patterns between different species, with precision, and on a large-scale basis. Studies on gene expression evolution are therefore performed only on small genes datasets, or using imprecise descriptions of expression patterns.The aim of my PhD was thus to develop and use novel bioinformatics resources, to study the evolution of gene expression. To this end, I developed the database Bgee (Base for Gene Expression Evolution). The approach of Bgee is to transform heterogeneous expression data (ESTs, microarrays, and in-situ hybridizations) into present/absent calls, and to annotate them to standard representations of anatomy and development of different species (anatomical ontologies). An extensive mapping between anatomies of species is then developed based on hypothesis of homology. These precise annotations to anatomies, and this extensive mapping between species, are the major assets of Bgee, and have required the involvement of many co-workers over the years. My main personal contribution is the development and the management of both the Bgee database and the web-application.Bgee is now on its ninth release, and includes an important gene expression dataset for 5 species (human, mouse, drosophila, zebrafish, Xenopus), with the most data from mouse, human and zebrafish. Using these three species, I have conducted an analysis of gene expression evolution after duplication in vertebrates.Gene duplication is thought to be a major source of novelty in evolution, and to participate to speciation. It has been suggested that the evolution of gene expression patterns might participate in the retention of duplicate genes. I performed a large-scale comparison of expression patterns of hundreds of duplicated genes to their singleton ortholog in an outgroup, including both small and large-scale duplicates, in three vertebrate species (human, mouse and zebrafish), and using highly accurate descriptions of expression patterns. My results showed unexpectedly high rates of de novo acquisition of expression domains after duplication (neofunctionalization), at least as high or higher than rates of partitioning of expression domains (subfunctionalization). I found differences in the evolution of expression of small- and large-scale duplicates, with small-scale duplicates more prone to neofunctionalization. Duplicates with neofunctionalization seemed to evolve under more relaxed selective pressure on the coding sequence. Finally, even with abundant and precise expression data, the majority fate I recovered was neither neo- nor subfunctionalization of expression domains, suggesting a major role for other mechanisms in duplicate gene retention.

Nitric Oxide Induces the Expression of the Monocarboxylate Transporter MCT4 in Cultured Astrocytes by a cGMP-Independent Transcriptional Activation

The monocarboxylate transporter MCT4 is a proton-linked carrier particularly important for lactate release from highly glycolytic cells. In the central nervous system, MCT4 is exclusively expressed by astrocytes. Surprisingly, MCT4 expression in primary cultures of mouse cortical astrocytes is conspicuously low, suggesting that an external, nonastrocytic signal is necessary to obtain the observed pattern of expression in vivo. Here, we demonstrate that nitric oxide (NO), delivered by various NO donors, time- and dose-dependently induces MCT4 expression in cultured cortical astrocytes both at the mRNA and protein levels. In contrast, NO does not enhance the expression of MCT1, the other astrocytic monocarboxylate transporter. The transcriptional effect of NO is not mediated by a cGMP-dependent mechanism as shown by the absence of effect of a cGMP analog or of a selective guanylate cyclase inhibitor. NO causes an increase in astrocytic lactate transport capacity which requires the enhancement of MCT4 expression as both are prevented by the use of a specific siRNA against MCT4. In addition, cumulated lactate release by astrocytes over a period of 24 h was also enhanced by NO treatment. Our data suggest that NO represents a putative intercellular signal to control MCT4 expression in astrocytes and in doing so, to facilitate lactate transfer to other surrounding cell types in the central nervous system. (C) 2011 Wiley-Liss, Inc.

Optimization Of Promoter Sequences In Order To Mimic Physiological Pde6b Expression

Purpose: Many retinal degenerations result from defective retina-specific gene expressions. Thus, it is important to understand how the expression of a photoreceptor-specific gene is regulated in vivo in order to achieve successful gene therapy. The present study aims to design an AAV2/8 vector that can regulate the transcript level in a physiological manner to replace missing PDE6b in Rd1 and Rd10 mice. In previous studies (Ogieta, et al., 2000), the short 5' flanking sequence of the human PDE6b gene (350 bp) was shown to be photoreceptor-specific in transgenic mice. However, the efficiency and specificity of the 5' flanking region of the human PDE6b was not investigated in the context of gene therapy during retinal degeneration. In this study, two different sequences of the 5' flanking region of the human PDE6b gene were studied as promoter elements and their expression will be tested in wild type and diseased retinas (Rd 10 mice).Methods: Two 5' flanking fragments of the human PDE6b gene: (-93 to +53 (150 bp) and -297 to +53 (350 bp)) were cloned in different plasmids in order to check their expression in vitro and in vivo by constructing an AAV2/8 vector. These elements drove the activity of either luciferase (pGL3 plasmids) or EGFP. jetPEI transfection in Y 79 cells was used to evaluate gene expression through luciferase activity. Constructs encoding EGFP under the control of the two promoters were performed in AAV2.1-93 (or 297)-EGFP plasmids to produce AAV2/8 vectors.Results: When pGL3-93 (150 bp) or pGL3-297 (350 bp) were transfected in the Y-79 cells, the smaller fragment (150 bp) showed higher gene expression compared to the 350 bp element and to the SV40 control, as previously reported. The 350 bp drove similar levels of expression when compared to the SV40 promoter. In view of these results, the fragments (150 bp or 350 bp) were integrated into the AAV2.1-EGFP plasmid to produce AAV2/8 vector, and we are currently evaluating the efficiency and specificity of the produced constructs in vivo in normal and diseased retinas.Conclusions: Comparisons of these vectors with vectors bearing ubiquitous promoters should reveal which construct is the most suitable to drive efficient and specific gene expression in diseased retinas in order to restore a normal function on the long term.

Essential role of Smad3 in the inhibition of inflammation-induced PPARbeta/delta expression.

Wound healing proceeds by the concerted action of a variety of signals that have been well identified. However, the mechanisms integrating them and coordinating their effects are poorly known. Herein, we reveal how PPARbeta/delta (PPAR: peroxisome proliferator-activated receptor) follows a balanced pattern of expression controlled by a crosstalk between inflammatory cytokines and TGF-beta1. Whereas conditions that mimic the initial inflammatory events stimulate PPARbeta/delta expression, TGF-beta1/Smad3 suppresses this inflammation-induced PPARbeta/delta transcription, as seen in the late re-epithelialization/remodeling events. This TGF-beta1/Smad3 action involves an inhibitory effect on AP-1 activity and DNA binding that results in an inhibition of the AP-1-driven induction of the PPARbeta/delta promoter. As expected from these observations, wound biopsies from Smad3-null mice showed sustained PPARbeta expression as compared to those of their wild-type littermates. Together, these findings suggest a mechanism for setting the necessary balance between inflammatory signals, which trigger PPARbeta/delta expression, and TGF-beta1/Smad3 that governs the timely decrease of this expression as wound healing proceeds to completion.

Postnatal expression pattern of calcium-binding proteins in organotypic thalamic cultures and in the dorsal thalamus in vivo.

The present study describes the postnatal expression of calbindin, calretinin and parvalbumin and glutamic acid decarboxylase (GAD) and microtubule-associated protein 2 (MAP2) in organotypic monocultures of rat dorsal thalamus compared to the thalamus in vivo. Cultures were maintained for up to 7 weeks. Cortex-conditioned medium improved the survival of thalamic cultures. MAP2-immunoreactive material was present in somata and dendrites of small and large-sized neurons throughout the cultures. Parvalbumin immunoreactivity was present in larger multipolar or bitufted neurons along the edge of a culture. These neurons also displayed strong parvalbumin mRNA and GAD mRNA expression, and GABA immunoreactivity. They likely corresponded to cells of the nucleus reticularis thalami. Parvalbumin mRNA, but neither parvalbumin protein nor GAD mRNA, was expressed in neurons with large somata within the explant. They likely represented relay cells. GAD mRNA, but not parvalbumin mRNA, was expressed in small neurons within the explants. Small neurons also displayed calbindin- and calretinin-immunoreactivity. The small neurons likely represented local circuit neurons. The time course of expression of the calcium-binding proteins revealed that all were present at birth with the predicted molecular weights. A low, but constant parvalbumin expression was observed in vitro without the developmental increase seen in vivo, which most likely represented parvalbumin from afferent sources. In contrast, the explantation transiently downregulated the calretinin and calbindin expression, but the neurons recovered the expression after 14 and 21 days, respectively. In conclusion, thalamic monocultures older than three weeks represent a stable neuronal network containing well differentiated neurons of the nucleus reticularis thalami, relay cells and local circuit neurons.

CleanEx: a database of heterogeneous gene expression data based on a consistent gene nomenclature and linked to an improved annotation system

Résumé: L'automatisation du séquençage et de l'annotation des génomes, ainsi que l'application à large échelle de méthodes de mesure de l'expression génique, génèrent une quantité phénoménale de données pour des organismes modèles tels que l'homme ou la souris. Dans ce déluge de données, il devient très difficile d'obtenir des informations spécifiques à un organisme ou à un gène, et une telle recherche aboutit fréquemment à des réponses fragmentées, voir incomplètes. La création d'une base de données capable de gérer et d'intégrer aussi bien les données génomiques que les données transcriptomiques peut grandement améliorer la vitesse de recherche ainsi que la qualité des résultats obtenus, en permettant une comparaison directe de mesures d'expression des gènes provenant d'expériences réalisées grâce à des techniques différentes. L'objectif principal de ce projet, appelé CleanEx, est de fournir un accès direct aux données d'expression publiques par le biais de noms de gènes officiels, et de représenter des données d'expression produites selon des protocoles différents de manière à faciliter une analyse générale et une comparaison entre plusieurs jeux de données. Une mise à jour cohérente et régulière de la nomenclature des gènes est assurée en associant chaque expérience d'expression de gène à un identificateur permanent de la séquence-cible, donnant une description physique de la population d'ARN visée par l'expérience. Ces identificateurs sont ensuite associés à intervalles réguliers aux catalogues, en constante évolution, des gènes d'organismes modèles. Cette procédure automatique de traçage se fonde en partie sur des ressources externes d'information génomique, telles que UniGene et RefSeq. La partie centrale de CleanEx consiste en un index de gènes établi de manière hebdomadaire et qui contient les liens à toutes les données publiques d'expression déjà incorporées au système. En outre, la base de données des séquences-cible fournit un lien sur le gène correspondant ainsi qu'un contrôle de qualité de ce lien pour différents types de ressources expérimentales, telles que des clones ou des sondes Affymetrix. Le système de recherche en ligne de CleanEx offre un accès aux entrées individuelles ainsi qu'à des outils d'analyse croisée de jeux de donnnées. Ces outils se sont avérés très efficaces dans le cadre de la comparaison de l'expression de gènes, ainsi que, dans une certaine mesure, dans la détection d'une variation de cette expression liée au phénomène d'épissage alternatif. Les fichiers et les outils de CleanEx sont accessibles en ligne (http://www.cleanex.isb-sib.ch/). Abstract: The automatic genome sequencing and annotation, as well as the large-scale gene expression measurements methods, generate a massive amount of data for model organisms. Searching for genespecific or organism-specific information througout all the different databases has become a very difficult task, and often results in fragmented and unrelated answers. The generation of a database which will federate and integrate genomic and transcriptomic data together will greatly improve the search speed as well as the quality of the results by allowing a direct comparison of expression results obtained by different techniques. The main goal of this project, called the CleanEx database, is thus to provide access to public gene expression data via unique gene names and to represent heterogeneous expression data produced by different technologies in a way that facilitates joint analysis and crossdataset comparisons. A consistent and uptodate gene nomenclature is achieved by associating each single gene expression experiment with a permanent target identifier consisting of a physical description of the targeted RNA population or the hybridization reagent used. These targets are then mapped at regular intervals to the growing and evolving catalogues of genes from model organisms, such as human and mouse. The completely automatic mapping procedure relies partly on external genome information resources such as UniGene and RefSeq. The central part of CleanEx is a weekly built gene index containing crossreferences to all public expression data already incorporated into the system. In addition, the expression target database of CleanEx provides gene mapping and quality control information for various types of experimental resources, such as cDNA clones or Affymetrix probe sets. The Affymetrix mapping files are accessible as text files, for further use in external applications, and as individual entries, via the webbased interfaces . The CleanEx webbased query interfaces offer access to individual entries via text string searches or quantitative expression criteria, as well as crossdataset analysis tools, and crosschip gene comparison. These tools have proven to be very efficient in expression data comparison and even, to a certain extent, in detection of differentially expressed splice variants. The CleanEx flat files and tools are available online at: http://www.cleanex.isbsib. ch/.

Learning-induced gene expression in the heads of two Nasonia species that differ in long-term memory formation.

BACKGROUND: Cellular processes underlying memory formation are evolutionary conserved, but natural variation in memory dynamics between animal species or populations is common. The genetic basis of this fascinating phenomenon is poorly understood. Closely related species of Nasonia parasitic wasps differ in long-term memory (LTM) formation: N. vitripennis will form transcription-dependent LTM after a single conditioning trial, whereas the closely-related species N. giraulti will not. Genes that were differentially expressed (DE) after conditioning in N. vitripennis, but not in N. giraulti, were identified as candidate genes that may regulate LTM formation. RESULTS: RNA was collected from heads of both species before and immediately, 4 or 24 hours after conditioning, with 3 replicates per time point. It was sequenced strand-specifically, which allows distinguishing sense from antisense transcripts and improves the quality of expression analyses. We determined conditioning-induced DE compared to naïve controls for both species. These expression patterns were then analysed with GO enrichment analyses for each species and time point, which demonstrated an enrichment of signalling-related genes immediately after conditioning in N. vitripennis only. Analyses of known LTM genes and genes with an opposing expression pattern between the two species revealed additional candidate genes for the difference in LTM formation. These include genes from various signalling cascades, including several members of the Ras and PI3 kinase signalling pathways, and glutamate receptors. Interestingly, several other known LTM genes were exclusively differentially expressed in N. giraulti, which may indicate an LTM-inhibitory mechanism. Among the DE transcripts were also antisense transcripts. Furthermore, antisense transcripts aligning to a number of known memory genes were detected, which may have a role in regulating these genes. CONCLUSION: This study is the first to describe and compare expression patterns of both protein-coding and antisense transcripts, at different time points after conditioning, of two closely related animal species that differ in LTM formation. Several candidate genes that may regulate differences in LTM have been identified. This transcriptome analysis is a valuable resource for future in-depth studies to elucidate the role of candidate genes and antisense transcription in natural variation in LTM formation.