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.

Plants respond to pathogen infection by enhancing the antifungal gene expression of root-associated bacteria.

Plant health and fitness widely depend on interactions with soil microorganisms. Some bacteria such as pseudomonads can inhibit pathogens by producing antibiotics, and controlling these bacteria could help improve plant fitness. In the present study, we tested whether plants induce changes in the antifungal activity of root-associated bacteria as a response to root pathogens. We grew barley plants in a split-root system with one side of the root system challenged by the pathogen Pythium ultimum and the other side inoculated with the biocontrol strain Pseudomonas fluorescens CHA0. We used reporter genes to follow the expression of ribosomal RNA indicative of the metabolic state and of the gene phlA, required for production of 2,4-diacetylphloroglucinol, a key component of antifungal activity. Infection increased the expression of the antifungal gene phlA. No contact with the pathogen was required, indicating that barley influenced gene expression by the bacteria in a systemic way. This effect relied on increased exudation of diffusible molecules increasing phlA expression, suggesting that communication with rhizosphere bacteria is part of the pathogen response of plants. Tripartite interactions among plants, pathogens, and bacteria appear as a novel determinant of plant response to root pathogens.

aIr-1, a newly found locus on mouse chromosome 16 encoding a trans-acting activator factor for MHC class II gene expression.

RJ 2.2.5 is a human B cell line that has lost the capacity to express MHC class II genes. The human class II-positive phenotype is restored in somatic cell hybrids between RJ 2.2.5 and mouse spleen cells. By karyotype and molecular studies of an informative family of hybrids we have now shown that the reexpression of human class II gene products, as well as the maintenance of the mouse class II-positive phenotype, correlates with the presence of mouse chromosome 16. Thus, the existence on this mouse chromosome of a newly found locus, designated by us aIr-1, that determines a trans-acting activator function for class II gene expression, is established. Possible implications of this finding are discussed.

Differences in glucose transporter gene expression between rat pancreatic alpha- and beta-cells are correlated to differences in glucose transport but not in glucose utilization.

Glucose exerts inverse effects upon the secretory function of islet alpha- and beta-cells, suppressing glucagon release and increasing insulin release. This diverse action may result from differences in glucose transport and metabolism between the two cell types. The present study compares glucose transport in rat alpha- and beta-cells. beta-Cells transcribed GLUT2 and, to a lesser extent, GLUT 1; alpha-cells contained GLUT1 but no GLUT2 mRNA. No other GLUT-like sequences were found among cDNAs from alpha- or beta-cells. Both cell types expressed 43-kDa GLUT1 protein which was enhanced by culture. The 62-kDa beta-cell GLUT2 protein was converted to a 58-kDa protein after trypsin treatment of the cells without detectable consequences upon glucose transport kinetics. In beta-cells, the rates of glucose transport were 10-fold higher than in alpha-cells. In both cell types, glucose uptake exceeded the rates of glucose utilization by a factor of 10 or more. Glycolytic flux, measured as D-[5(3)H]glucose utilization, was comparable in alpha- and beta-cells between 1 and 10 mmol/liter substrate. In conclusion, differences in glucose transporter gene expression between alpha- and beta-cells can be correlated with differences in glucose transport kinetics but not with different glucose utilization rates.

Critical role for Ets, AP-1 and GATA-like transcription factors in regulating mouse Toll-like receptor 4 (Tlr4) gene expression.

TLR4 (Toll-like receptor 4) is essential for sensing the endotoxin of Gram-negative bacteria. Mutations or deletion of the TLR4 gene in humans or mice have been associated with altered predisposition to or outcome of Gram-negative sepsis. In the present work, we studied the expression and regulation of the Tlr4 gene of mouse. In vivo, TLR4 levels were higher in macrophages compared with B, T or natural killer cells. High basal TLR4 promoter activity was observed in RAW 264.7, J774 and P388D1 macrophages transfected with a TLR4 promoter reporter vector. Analysis of truncated and mutated promoter constructs identified several positive [two Ets (E twenty-six) and one AP-1 (activator protein-1) sites] and negative (a GATA-like site and an octamer site) regulatory elements within 350 bp upstream of the transcriptional start site. The myeloid and B-cell-specific transcription factor PU.1 bound to the proximal Ets site. In contrast, none among PU.1, Ets-1, Ets-2 and Elk-1, but possibly one member of the ESE (epithelium-specific Ets) subfamily of Ets transcription factors, bound to the distal Ets site, which was indispensable for Tlr4 gene transcription. Endotoxin did not affect macrophage TLR4 promoter activity, but it decreased TLR4 steady-state mRNA levels by increasing the turnover of TLR4 transcripts. TLR4 expression was modestly altered by other pro- and anti-inflammatory stimuli, except for PMA plus ionomycin which strongly increased promoter activity and TLR4 mRNA levels. The mouse and human TLR4 genes were highly conserved. Yet, notable differences exist with respect to the elements implicated in gene regulation, which may account for species differences in terms of tissue expression and modulation by microbial and inflammatory stimuli.

A role for altered TLR gene expression in association with increased expression of CD200R in the induction of mucosal tissue CD4(+) Treg in aged mice following gavage with a liver extract along with intramuscular monophosphoryl lipid A (MPLA) injection.

Previous studies showed a fetal sheep liver extract (FSLE), in association with LPS, injected into aged (>20 months) mice reversed the altered polarization (increased IL-4 and IL-10 with decreased IL-2 and IFN-gamma) in cytokine production seen from ConA stimulated lymphoid cells of those mice. Aged mice show a >60% decline in numbers and suppressive function of both CD4(+)CD25(+)Foxp3(+)Treg and so-called Tr3 (CD4(+)TGFbeta(+)). Their number/function is restored to levels seen in control (8-week-old) mice by FSLE. We have reported at length on the ability of a novel pair of immunoregulatory molecules, members of the TREM family, namely CD200:CD200R, to control development of dendritic cells (DCs) which themselves regulate production of Foxp3(+) Treg. The latter express a distinct subset of TLRs which control their function. We report that a feature of the altered Treg expression following combined treatment with FSLE and monophosphoryl lipid A, MPLA (a bioactive component of lipid A of LPS) is the altered gene expression both of distinct subsets of TLRs and of CD200Rs. We speculate that this may represent one of the mechanisms by which FSLE and MPLA alter immunity in aged mice.

Na(+)-K(+)-ATPase gene expression during in vitro development of rat fetal forebrain.

The existence of at least three isoforms of Na(+)-K(+)-ATPase in adult brain tissues [alpha 1, kidney type; alpha 2 [or alpha(+)]; alpha 3] suggests that these genes might be regulated in a cell-specific and time-dependent manner during development. We have studied this question in serum-free aggregating cell cultures of mechanically dissociated rat fetal telencephalon. At the protein level, the relative rate of synthesis of the pool of alpha 1-, alpha 2-, and alpha 3-subunits increased approximately twofold over 15 days of culture, leading to a marked increase in the immunochemical pool of alpha-subunits as measured by a panspecific polyclonal antibody. Concomitantly, Na(+)-K(+)-ATPase enzyme-specific activity increased three- (lower forebrain) to sixfold (upper forebrain). The transcripts of all three alpha-isoforms and beta-subunit were detected in vitro in similar proportion to the level observed in vivo. alpha 3-mRNA (3.7 kb) was more abundant than alpha 1 (3.7 kb) or alpha 2 (5.3 and 3.4 kb). Cytosine arabinoside (0.4 microM) and cholera toxin (0.1 microM) were used to selectively eliminate glial cells or neurons, respectively. It was found that alpha 2-mRNA is predominantly transcribed in glial cell cultures, whereas alpha 3- and beta 1-mRNA (2.7, 2.3, and 1.8 kb) are predominant in neuronal cultures.

" Arabidopsis Thaliana " response to insect feeding : new components controlling defense gene expression and plant resistance

Abstract: Plants cannot run away to escape attacking herbivores, but they defend themselves by producing anti-digestive proteins and toxic compounds (for example glucosinolates). The first goal of this thesis was to study changes in gene expression after insect attack using microarrays. The responses of Arabidopsis thaliana to feeding by the specialist Pieris rapae and the generalist Spodoptera liffora is were compared. We found that the transcript profiles after feeding by the two chewing insects were remarkably similar, although the generalist induced a slightly stronger response. The second goal was to evaluate the implication of the four signals jasmonic acid (JA), salicylic acid (SA), ethylene (ET), and abscisic acid (ABA) in the control of insect-regulated gene expression. Using signaling mutants, we observed that JA was the predominant signal and that ABA modulated defense gene expression. In contrast, SA and ET appeared to control slightly gene expression, but only after feeding by S. litforalis. The third goal was to establish whether plant responses are really effective against insects. In accordance with the transcript profile, both insects were affected by the JA-dependent defenses, as they performed better on the JA-insensitive mutant. S. littoralis also performed better on ABA-deficient mutants, providing evidence for the role of ABA in defense against insects. When testing indole or aliphatic glucosinolate deficient mutants, we found that they were also more susceptible to insect feeding, providing some of the first genetic evidence for the defensive role of glucosinolates in planta. Finally, a glutathione-deficient mutant, pad2-1, was also more susceptible to insect feeding and we could attribute this phenotype to a lowered accumulation of the major indole glucosinolate. In this thesis, we provide a comprehensive list of insect-regulated genes, including many transcription factors that constitute interesting candidate genes for the further study of insect-induced expression changes. Understanding how the plant responses to insects are regulated will provide tools for a better management of insect pest in the field. Résumé: Les plantes ne peuvent s'échapper pour fuir les insectes qui les attaquent, mais elles se défendent en produisant des protéines anti-digestives et des composés toxiques (par exemple des glucosinolates). Le premier but de cette thèse était d'étudier les changements de l'expression génétique lors d'attaque par des insectes en utilisant des puces à ADN. Nous avons comparé la réponse d'Arabidopsis thaliana à deux espèces d'insectes avec des habitudes alimentaires différentes : le spécialiste Pieris rapae et le généraliste Spodoptera littoralis. Nous avons trouvé que les profils de transcription après l'attaque par les deux insectes sont remarquablement similaires, bien que le généraliste induise une réponse légèrement plus forte. Le deuxième but était de déterminer l'implication de quatre signaux dans le contrôle de la réponse :l'acide jasmonique (JA), l'acide salicylique (SA), l'éthylène (ET), et l'acide abscissique (ABA). En utilisant de mutants de signalisation, nous avons montré que l'acide jasmonique était le signal prédominant et que l'acide abscissique modulait l'expression génétique. D'autre part, l'acide salicylique et l'éthylène contrôlent à un degré moindre l'expression génétique, mais seulement après l'attaque par S. littoralís. Le troisième but était d'établir si les réponses des plantes sont efficaces contre les insectes. En accord avec le profil de transcription, les deux espèces d'insectes se sont mieux développées sur un mutant insensible au JA, indiquant que les défenses contrôlées par ce signal sont cruciales pour la plante. De plus, les larves de S. littorales se sont mieux développées sur des mutants déficients en ABA, ce qui fournit une preuve du rôle de l'acide abscissique dans la défense contre les insectes. En testant des mutants déficients en glucosinolates de type indole ou aliphatique, nous avons trouvé qu'ils étaient plus sensibles aux insectes, démontrant ainsi le rôle défensif des glucosinolates in planta. Finalement, le mutant déficient en glutathion pad2-1 était aussi plus sensible à l'attaque des insectes, et nous avons pu attribuer ce phénotype à une plus faible augmentation d'un indole glucosinolate dans ce mutant. Dans cette thèse, nous avons mis en évidence un nombre important de gènes contrôlés par les insectes, comprenant de nombreux facteurs de transcription qui constituent des candidats intéressants pour`étudier plus en détail les changements d'expression génétique induits par les insectes. Une meilleure compréhension de la réponse des plantes contre l'attaque des insectes devrait nous permettre de développer de nouvelles stratégies pour mieux gérer les ravageurs des cultures.

The mammalian circadian timing system: from gene expression to physiology.

Many physiological processes in organisms from bacteria to man are rhythmic, and some of these are controlled by self-sustained oscillators that persist in the absence of external time cues. Circadian clocks are perhaps the best characterized biological oscillators and they exist in virtually all light-sensitive organisms. In mammals, they influence nearly all aspects of physiology and behavior, including sleep-wake cycles, cardiovascular activity, endocrinology, body temperature, renal activity, physiology of the gastro-intestinal tract, and hepatic metabolism. The master pacemaker is located in the suprachiasmatic nuclei, two small groups of neurons in the ventral part of the hypothalamus. However, most peripheral body cells contain self-sustained circadian oscillators with a molecular makeup similar to that of SCN (suprachiasmatic nucleus) neurons. This organization implies that the SCN must synchronize countless subsidiary oscillators in peripheral tissues, in order to coordinate cyclic physiology. In this review, we will discuss some recent studies on the structure and putative functions of the mammalian circadian timing system, but we will also point out some apparent inconsistencies in the currently publicized model for rhythm generation.

A comprehensive analysis of human gene expression profiles identifies stromal immunoglobulin κ C as a compatible prognostic marker in human solid tumors.

PURPOSE: Although the central role of the immune system for tumor prognosis is generally accepted, a single robust marker is not yet available. EXPERIMENTAL DESIGN: On the basis of receiver operating characteristic analyses, robust markers were identified from a 60-gene B cell-derived metagene and analyzed in gene expression profiles of 1,810 breast cancer; 1,056 non-small cell lung carcinoma (NSCLC); 513 colorectal; and 426 ovarian cancer patients. Protein and RNA levels were examined in paraffin-embedded tissue of 330 breast cancer patients. The cell types were identified with immunohistochemical costaining and confocal fluorescence microscopy. RESULTS: We identified immunoglobulin κ C (IGKC) which as a single marker is similarly predictive and prognostic as the entire B-cell metagene. IGKC was consistently associated with metastasis-free survival across different molecular subtypes in node-negative breast cancer (n = 965) and predicted response to anthracycline-based neoadjuvant chemotherapy (n = 845; P < 0.001). In addition, IGKC gene expression was prognostic in NSCLC and colorectal cancer. No association was observed in ovarian cancer. IGKC protein expression was significantly associated with survival in paraffin-embedded tissues of 330 breast cancer patients. Tumor-infiltrating plasma cells were identified as the source of IGKC expression. CONCLUSION: Our findings provide IGKC as a novel diagnostic marker for risk stratification in human cancer and support concepts to exploit the humoral immune response for anticancer therapy. It could be validated in several independent cohorts and carried out similarly well in RNA from fresh frozen as well as from paraffin tissue and on protein level by immunostaining.

Concordance among digital gene expression, microarrays, and qPCR when measuring differential expression of microRNAs.

Profiling microRNA (miRNA) expression is of widespread interest given the critical role of miRNAs in many cellular functions. Profiling can be achieved via hybridization-based (microarrays), sequencing-based, or amplification-based (quantitative reverse transcription-PCR, qPCR) technologies. Among these, microarrays face the significant challenge of accurately distinguishing between mature and immature miRNA forms, and different vendors have developed different methods to meet this challenge. Here we measure differential miRNA expression using the Affymetrix, Agilent, and Illumina microarray platforms, as well as qPCR (Applied Biosystems) and ultra high-throughput sequencing (Illumina). We show that the differential expression measurements are more divergent when the three types of microarrays are compared than when the Agilent microarray, qPCR, and sequencing technology measurements are compared, which exhibit a good overall concordance.

Evolution of gene expression in fire ants: the effects of developmental stage, caste, and species.

Ants provide remarkable examples of equivalent genotypes developing into divergent and discrete phenotypes. Diploid eggs can develop either into queens, which specialize in reproduction, or workers, which participate in cooperative tasks such as building the nest, collecting food, and rearing the young. In contrast, the differentiation between males and females generally depends upon whether eggs are fertilized, with fertilized (diploid) eggs giving rise to females and unfertilized (haploid) eggs giving rise to males. To obtain a comprehensive picture of the relative contributions of gender (sex), caste, developmental stage, and species divergence to gene expression evolution, we investigated gene expression patterns in pupal and adult queens, workers, and males of two species of fire ants, Solenopsis invicta and S. richteri. Microarray hybridizations revealed that variation in gene expression profiles is influenced more by developmental stage than by caste membership, sex, or species identity. The second major contributor to variation in gene expression was the combination of sex and caste. Although workers and queens share equivalent diploid nuclear genomes, they have highly distinctive patterns of gene expression in both the pupal and the adult stages, as might be expected given their extraordinary level of phenotypic differentiation. Overall, the difference in the proportion of differentially expressed genes was greater between workers and males than between workers and queens or queens and males, consistent with the fact that workers and males share neither gender nor reproductive capability. Moreover, between-species comparisons revealed that the greatest difference in gene expression patterns occurred in adult workers, a finding consistent with the fact that adult workers most directly experience the distinct external environments characterizing the different habitats occupied by the two species. Thus, much of the evolution of gene expression in ants may occur in the worker caste, despite the fact that these individuals are largely or completely sterile. Analyses of gene expression evolution revealed a combination of positive selection and relaxation of stabilizing selection as important factors driving the evolution of such genes.

Detection of plant-modulated alterations in antifungal gene expression in Pseudomonas fluorescens CHA0 on roots by flow cytometry.

The biocontrol activity of the root-colonizing Pseudomonas fluorescens strain CHA0 is largely determined by the production of antifungal metabolites, especially 2,4-diacetylphloroglucinol. The expression of these metabolites depends on abiotic and biotic environmental factors, in particular, elements present in the rhizosphere. In this study, we have developed a new method for the in situ analysis of antifungal gene expression using flow cytometry combined with green fluorescent protein (GFP)-based reporter fusions to the phlA and prnA genes essential for the production of the antifungal compounds 2,4-diacetylphloroglucinol and pyrrolnitrin, respectively, in strain CHA0. Expression of phlA-gfp and prnA-gfp in CHA0 cells harvested from the rhizosphere of a set of plant species as well as from the roots of healthy, leaf pathogen-attacked, and physically stressed plants were analyzed using a FACSCalibur. After subtraction of background fluorescence emitted by plant-derived particles and CHA0 cells not carrying the gfp reporters, the average gene expression per bacterial cell could be calculated. Levels of phlA and prnA expression varied significantly in the rhizospheres of different plant species. Physical stress and leaf pathogen infection lowered phlA expression levels in the rhizosphere of cucumber. Our results demonstrate that the newly developed approach is suitable to monitor differences in levels of antifungal gene expression in response to various plant-derived factors. An advantage of the method is that it allows quantification of bacterial gene expression in rhizosphere populations at a single-cell level. To our best knowledge, this is the first study using flow cytometry for the in situ analysis of biocontrol gene expression in a plant-beneficial bacterium in the rhizosphere.

Differential regulation of Na-K-ATPase isoform gene expression by T3 during rat brain development.

A fetal rat telencephalon organotypic cell culture system was found to reproduce the developmental pattern of Na-K-adenosinetriphosphatase (ATPase) gene expression observed in vivo [Am. J. Physiol. 258 (Cell Physiol. 27): C1062-C1069, 1990]. We have used this culture system to study the effects of triiodothyronine (T3; 0.003-30 nM) on mRNA abundance and basal transcription rates of Na-K-ATPase isoforms. Steady-state mRNA levels were low at culture day 6 (corresponding to the day of birth) but distinct for each isoform alpha 3 much greater than beta 1 = beta 2 greater than alpha 2 greater than alpha 1. At culture day 6, T3 did not modify mRNA abundance of any isoform. At culture day 12 (corresponding to day 7 postnatal), T3 increased the mRNA level of alpha 2 (4- to 7-fold), beta 2 (4- to 5-fold), alpha 1 (3- to 6-fold), and beta 1 (1.5-fold), whereas alpha 3 mRNA levels remained unchanged. Interestingly, the basal transcription rate for each isoform differed strikingly (alpha 2 greater than alpha 1 much greater than beta 1 = beta 2 greater than alpha 3) but remained stable throughout 12 days of culture and was not regulated by T3. Thus we observed an inverse relationship between rate of transcription and rate of mRNA accumulation for each alpha-isoform, suggesting that alpha 1- and alpha 2-mRNA are turning over rapidly whereas alpha 3-mRNA is turning over slowly. Our data indicate that one of the mechanisms by which T3 selectively controls Na-K-ATPase gene expression during brain development in vitro occurs at the posttranscriptional level.

Comparative modular analysis of gene expression in vertebrate development

The focus of my PhD research was the concept of modularity. In the last 15 years, modularity has become a classic term in different fields of biology. On the conceptual level, a module is a set of interacting elements that remain mostly independent from the elements outside of the module. I used modular analysis techniques to study gene expression evolution in vertebrates. In particular, I identified ``natural'' modules of gene expression in mouse and human, and I showed that expression of organ-specific and system-specific genes tends to be conserved between such distance vertebrates as mammals and fishes. Also with a modular approach, I studied patterns of developmental constraints on transcriptome evolution. I showed that none of the two commonly accepted models of the evolution of embryonic development (``evo-devo'') are exclusively valid. In particular, I found that the conservation of the sequences of regulatory regions is highest during mid-development of zebrafish, and thus it supports the ``hourglass model''. In contrast, events of gene duplication and new gene introduction are most rare in early development, which supports the ``early conservation model''. In addition to the biological insights on transcriptome evolution, I have also discussed in detail the advantages of modular approaches in large-scale data analysis. Moreover, I re-analyzed several studies (published in high-ranking journals), and showed that their conclusions do not hold out under a detailed analysis. This demonstrates that complex analysis of high-throughput data requires a co-operation between biologists, bioinformaticians, and statisticians.