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.

Identification of reference genes for expression analysis by real-time quantitative PCR in drought-stressed soybean

The objective of this work was to validate, by quantitative PCR in real time (RT-qPCR), genes to be used as reference in studies of gene expression in soybean in drought-stressed trials. Four genes commonly used in soybean were evaluated: Gmβ-actin, GmGAPDH, GmLectin and GmRNAr18S. Total RNA was extracted from six samples: three from roots in a hydroponic system with different drought intensities (0, 25, 50, 75 and 100 minutes of water stress), and three from leaves of plants grown in sand with different soil moistures (15, 5 and 2.5% gravimetric humidity). The raw cycle threshold (Ct) data were analyzed, and the efficiency of each primer was calculated for an overall analysis of the Ct range among the different samples. The GeNorm application was used to evaluate the best reference gene, according to its stability. The GmGAPDH was the least stable gene, with the highest mean values of expression stability (M), and the most stable genes, with the lowest M values, were the Gmβ-actin and GmRNAr18S, when both root and leaves samples were tested. These genes can be used in RT-qPCR as reference gene for expression analysis.

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/.

Conserved developmental expression of Fezf in chordates and Drosophila and the origin of the Zona Limitans Intrathalamica (ZLI) brain organizer

Background: The zona limitans intrathalamica (ZLI) and the isthmus organizer (IsO) are two major secondary organizers of vertebrate brain development. These organizers are located at the interface of the expression domains of key patterning genes (Fezf-Irx and Otx-Gbx, respectively). To gain insights into the evolutionary origin of the ZLI, we studied Fezf in bilaterians. Results: In this paper, we identified a conserved sequence motif (Fezf box) in all bilaterians. We report the expression pattern of Fezf in amphioxus and Drosophila and compare it with those of Gbx, Otx and Irx. We found that the relative expression patterns of these genes in vertebrates are fully conserved in amphioxus and flies, indicating that the genetic subdivisions defining the location of both secondary organizers in early vertebrate brain development were probably present in the last common ancestor of extant bilaterians. However, in contrast to vertebrates, we found that Irx-defective flies do not show an affected Fezf expression pattern. Conclusions: The absence of expression of the corresponding morphogens from cells at these conserved genetic boundaries in invertebrates suggests that the organizing properties might have evolved specifically in the vertebrate lineage by the recruitment of key morphogens to these conserved genetic locations.

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.

Conserved developmental expression of Fezf in chordates and Drosophila and the origin of the Zona Limitans Intrathalamica (ZLI) brain organizer

Background: The zona limitans intrathalamica (ZLI) and the isthmus organizer (IsO) are two major secondary organizers of vertebrate brain development. These organizers are located at the interface of the expression domains of key patterning genes (Fezf-Irx and Otx-Gbx, respectively). To gain insights into the evolutionary origin of the ZLI, we studied Fezf in bilaterians. Results: In this paper, we identified a conserved sequence motif (Fezf box) in all bilaterians. We report the expression pattern of Fezf in amphioxus and Drosophila and compare it with those of Gbx, Otx and Irx. We found that the relative expression patterns of these genes in vertebrates are fully conserved in amphioxus and flies, indicating that the genetic subdivisions defining the location of both secondary organizers in early vertebrate brain development were probably present in the last common ancestor of extant bilaterians. However, in contrast to vertebrates, we found that Irx-defective flies do not show an affected Fezf expression pattern. Conclusions: The absence of expression of the corresponding morphogens from cells at these conserved genetic boundaries in invertebrates suggests that the organizing properties might have evolved specifically in the vertebrate lineage by the recruitment of key morphogens to these conserved genetic locations.

How to introduce connectance in a frame of an expression for diversity

Systems made of parts that are totally connected do not work, neither ecosys- tems nor artifacts. Relative connectance is inversely related to diversity, and both magnitudes can find a common frame of expression, in which some constant expressing the constraints of any organization might be embodied. If S is Simp- son's index, the expression (1 - S)IS as a measure of diversity offers some advantages or, at least, helps further reasoning. Such expression is the ratio between total interspecific possible interactions and possible intraspecific inter- actions.

Effect of iodide on Fas, Fas-ligand and Bcl-w mRNA expression in thyroid of NOD mice pretreated with methimazole

Nonobese diabetic (NOD) mice and a derived strain, NOD.H.2h4, have been used as a model for experimental spontaneous thyroiditis and thyroiditis induced by iodide excess after a goiter-inducing period. Some authors have proposed that iodide, given after methimazole or propylthiouracil, is capable of inducing apoptosis in thyroid cells and that anti-thyroid drugs can modulate the expression of apoptosis components such as Fas and its ligand (Fas-L). Here we evaluated the effect of potassium iodide (20 µg/animal for 4 days, ip) given to NOD mice at the 10th week of life after exposure to methimazole (1 mg/ml) in drinking water from the 4th to the 10th week of life. Fas, Fas-L and Bcl-w expression were analyzed semiquantitatively by RT-PCR immediately after potassium iodide administration (group MI44D) or at week 32 (MI32S). Control groups were added at 10 (C10) and 32 weeks (C32), as well as a group that received only methimazole (CM10). An increase in the expression of Fas-L and Bcl-w (P<0.01, ANOVA) was observed in animals of group MI44D, while Fas was expressed at higher levels (P = 0.02) in group C32 (72.89 ± 47.09 arbitrary units) when compared to group C10 (10.8 ± 8.55 arbitrary units). Thus, the analysis of Fas-L and Bcl-w expression in the MI44D group and Fas in group C32 allowed us to detect two different patterns of expression of these apoptosis components in thyroid tissue of NOD mice.

Expression of neuronal nitric oxide synthase in the developing superficial layers of the rat superior colliculus

We investigated the level of expression of neuronal nitric oxide synthase (nNOS) in the retinorecipient layers of the rat superior colliculus during early postnatal development. Male and female Lister rats ranging in age between the day of birth (P0) and the fourth postnatal week were used in the present study. Two biochemical methods were used, i.e., in vitro measurement of NOS specific activity by the conversion of [³H]-arginine to [³H]-citrulline, and analysis of Western blotting immunoreactive bands from superior colliculus homogenates. As revealed by Western blotting, very weak immunoreactive bands were observed as early as P0-2, and their intensity increased progressively at least until P21. The analysis of specific activity of NOS showed similar results. There was a progressive increase in enzymatic activity until near the end of the second postnatal week, and a nonsignificant tendency to an increase until the end of the third week was also observed. Thus, these results indicated an increase in the amount of nNOS during the first weeks after birth. Our results confirm and extend previous reports using histochemistry for NADPH-diaphorase and immunocytochemistry for nNOS, which showed a progressive increase in the number of stained cells in the superficial layers during the first two postnatal weeks, reaching an adult pattern at the end of the third week. Furthermore, our results suggested that nNOS is present in an active form in the rat superior colliculus during the period of refinement of the retinocollicular pathway.

A role for histone-like protein H1 (H-NS) in the regulation of hemolysin expression by Serratia marcescens

The histone-like protein H1 (H-NS) is an abundant structural component of the bacterial nucleoid and influences many cellular processes including recombination, transcription and transposition. Mutations in the hns gene encoding H-NS are highly pleiotropic, affecting the expression of many unrelated genes. We have studied the role of H-NS on the regulation of hemolysin gene expression in Serratia marcescens. The Escherichia coli hns mutant carrying S. marcescens hemolysin genes on a plasmid constructed by ligation of the 3.2-kb HindIII-SacI fragment of pR02 into pBluescriptIIKS, showed a high level of expression of this hemolytic factor. To determine the osmoregulation of wild-type and hns defective mutants the cells were grown to mid-logarithmic phase in LB medium with 0.06 or 0.3 M NaCl containing ampicillin and kanamycin, whereas to analyze the effect of pH on hemolysin expression, the cells were grown to late-logarithmic phase in LB medium buffered with 0.1 M Tris-HCl, pH 4.5 to 8.0. To assay growth phase-related hemolysin production, bacterial cells were grown in LB medium supplemented with ampicillin and kanamycin. The expression of S. marcescens hemolysin genes in wild-type E. coli and in an hns-defective derivative at different pH and during different growth phases indicated that, in the absence of H-NS, the expression of hemolysin did not vary with pH changes or growth phases. Furthermore, the data suggest that H-NS may play an important role in the regulation of hemolysin expression in S. marcescens and its effect may be due to changes in DNA topology influencing transcription and thus the amount of hemolysin expression. Implications for the mechanism by which H-NS influences gene expression are discussed.

Modulation of epithelial sodium channel (ENaC) expression in mouse lung infected with Pseudomonas aeruginosa

Affiliation: André Dagenais: Centre hospitalier de l'Université de Montréal/ Hôtel-Dieu, Département de médecine, Université de Montréal. Yves Berthiaume: Médecine et spécialités médicales, Faculté de médecine