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.

Expression of calbindin immunoreactivity by subpopulations of primary sensory neurons in chick embryo dorsal root ganglion cells grown in coculture or conditioned medium.

Primary sensory neurons which innervate neuromuscular spindles in the chicken are calbindin-immunoreactive. The influence exerted by developing skeletal muscle on the expression of calbindin immunoreactivity by subpopulations of dorsal root ganglion (DRG) cells in the chick embryo was tested in vitro in coculture with myoblasts, in conditioned medium (CM) prepared from myoblasts and in control cultures of DRG cells alone. Control cultures of DRG cells grown at the 6th embryonic day (E6) did not show any calbindin-immunostained ganglion cell. In coculture of myoblasts previously grown for 14 days, about 3% of calbindin-immunoreactive ganglion cells were detected while about 1% were observed in some cultures grown in CM. Fibroblasts from various sources were devoid of effect. Skin or kidney cells were more active than myoblasts to initiate calbindin expression by subpopulations of DRG cells in coculture or, to a lesser degree, in CM. The results suggest that cellular factors would rather induce calbindin expression in certain sensory neurons than ensure a selective neuronal survival.

MHC class I expression dependent on bacterial infection and parental factors in whitefish embryos (Salmonidae).

Ecological conditions can influence not only the expression of a phenotype, but also the heritability of a trait. As such, heritable variation for a trait needs to be studied across environments. We have investigated how pathogen challenge affects the expression of MHC genes in embryos of the lake whitefish Coregonus palaea. In order to experimentally separate paternal (i.e. genetic) from maternal and environmental effects, and determine whether and how stress affects the heritable variation for MHC expression, embryos were produced in full-factorial in vitro fertilizations, reared singly, and exposed at 208 degree days (late-eyed stage) to either one of two strains of Pseudomonas fluorescens that differ in their virulence characteristics (one increased mortality, while both delayed hatching time). Gene expression was assessed 48 h postinoculation, and virulence effects of the bacterial infection were monitored until hatching. We found no evidence of MHC class II expression at this stage of development. MHC class I expression was markedly down-regulated in reaction to both pseudomonads. While MHC expression could not be linked to embryo survival, the less the gene was expressed, the earlier the embryos hatched within each treatment group, possibly due to trade-offs between immune function and developmental rate or further factors that affect both hatching timing and MHC expression. We found significant additive genetic variance for MHC class I expression in some treatments. That is, changes in pathogen pressures could induce rapid evolution in MHC class I expression. However, we found no additive genetic variance in reaction norms in our study population.

The evolution of lncRNA repertoires and expression patterns in tetrapods.

Only a very small fraction of long noncoding RNAs (lncRNAs) are well characterized. The evolutionary history of lncRNAs can provide insights into their functionality, but the absence of lncRNA annotations in non-model organisms has precluded comparative analyses. Here we present a large-scale evolutionary study of lncRNA repertoires and expression patterns, in 11 tetrapod species. We identify approximately 11,000 primate-specific lncRNAs and 2,500 highly conserved lncRNAs, including approximately 400 genes that are likely to have originated more than 300 million years ago. We find that lncRNAs, in particular ancient ones, are in general actively regulated and may function predominantly in embryonic development. Most lncRNAs evolve rapidly in terms of sequence and expression levels, but tissue specificities are often conserved. We compared expression patterns of homologous lncRNA and protein-coding families across tetrapods to reconstruct an evolutionarily conserved co-expression network. This network suggests potential functions for lncRNAs in fundamental processes such as spermatogenesis and synaptic transmission, but also in more specific mechanisms such as placenta development through microRNA production.

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.

Influence of environmental factors on the expression of phenotypic characters by chicken dorsal root ganglion cells in culture.

Primary sensory neurons were grown under four conditions of culture. The influence of nonneuronal cells, horse serum or both was studied on the phenotypic expression of certain neuronal subpopulations. The number of neurons expressing acetylcholinesterase, alpha-bungarotoxin-binding sites or a high uptake capacity for glutamine was enhanced by nonneuronal cells. The horse serum increases the neuronal subpopulation exhibiting a carbonic anhydrase activity. Certain phenotypic changes fit conditions consistent with an epigenetic induction rather than a cell selection.

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.

Fas and Fas ligand in embryos and adult mice: ligand expression in several immune-privileged tissues and coexpression in adult tissues characterized by apoptotic cell turnover.

The cell surface receptor Fas (FasR, Apo-1, CD95) and its ligand (FasL) are mediators of apoptosis that have been shown to be implicated in the peripheral deletion of autoimmune cells, activation-induced T cell death, and one of the two major cytolytic pathways mediated by CD8+ cytolytic T cells. To gain further understanding of the Fas system., we have analyzed Fas and FasL expression during mouse development and in adult tissues. In developing mouse embryos, from 16.5 d onwards, Fas mRNA is detectable in distinct cell types of the developing sinus, thymus, lung, and liver, whereas FasL expression is restricted to submaxillary gland epithelial cells and the developing nervous system. Significant Fas and FasL expression were observed in several nonlymphoid cell types during embryogenesis, and generally Fas and FasL expression were not localized to characteristic sites of programmed cell death. In the adult mouse, RNase protection analysis revealed very wide expression of both Fas and FasL. Several tissues, including the thymus, lung, spleen, small intestine, large intestine, seminal vesicle, prostate, and uterus, clearly coexpress the two genes. Most tissues constitutively coexpressing Fas and FasL in the adult mouse are characterized by apoptotic cell turnover, and many of those expressing FasL are known to be immune privileged. It may be, therefore, that the Fas system is implicated in both the regulation of physiological cell turnover and the protection of particular tissues against potential lymphocyte-mediated damage.

Expression of genes encoding the pre-TCR and CD3 complex during thymus development.

The mature TCR is composed of a clonotypic heterodimer (alpha beta or gamma delta) associated with the invariant CD3 components (gamma, delta, epsilon and zeta). There is now considerable evidence that more immature forms of the TCR-CD3 complex (consisting of either CD3 alone or CD3 associated with a heterodimer of TCR beta and pre-T alpha) can be expressed at the cell surface on early thymocytes. These pre-TCR complexes are believed to be necessary for the ordered progression of early T cell development. We have analyzed in detail the expression of both the pre-TCR and CD3 complex at various stages of adult thymus development. Our data indicate that all CD3 components are already expressed at the mRNA level by the earliest identifiable (CD4lo) thymic precursor. In contrast, genes encoding the pre-TCR complex (pre-T alpha and fully rearranged TCR beta) are first expressed at the CD44loCD25+CD4-CD8- stage. Detectable surface expression of both CD3 and TCR beta are delayed relative to expression of the corresponding genes, suggesting the existence of other (as yet unidentified) components of the pre-TCR complex.

T cell receptor delta gene rearrangement and T early alpha (TEA) expression in immature alpha beta lineage thymocytes: implications for alpha beta/gamma delta lineage commitment.

Mature T cells comprise two mutually exclusive lineages expressing heterodimeric alpha beta or gamma delta antigen receptors. During development, beta, gamma, and delta genes rearrange before alpha, and mature gamma delta cells arise in the thymus prior to alpha beta cells. The mechanism underlying commitment of immature T cells to the alpha beta or gamma delta lineage is controversial. Since the delta locus is located within the alpha locus, rearrangement of alpha genes leads to deletion of delta. We have examined the rearrangement status of the delta locus immediately prior to alpha rearrangement. We find that many thymic precursors of alpha beta cells undergo VDJ delta rearrangements. Furthermore, the same cells frequently coexpress sterile T early alpha (TEA) transcripts originating 3' of C delta and 5' of the most upstream J alpha, thus implying that individual alpha beta lineage cells undergo sequential VDJ delta and VJ alpha rearrangements. Finally, VDJ delta rearrangements in immature alpha beta cells appear to be random, supporting models in which alpha beta lineage commitment is determined independently of the rearrangement status at the TCR delta locus.

Visualizing olfactory receptor expression and localization in Drosophila.

Odor detection and discrimination by olfactory systems in vertebrates and invertebrates depend both on the selective expression of individual olfactory receptor genes in subpopulations of olfactory sensory neurons, and on the targeting of the encoded proteins to the exposed, ciliated endings of sensory dendrites. Techniques to visualize the expression and localization of olfactory receptor gene products in vivo have been essential to reveal the molecular logic of peripheral odor coding and to permit investigation of the developmental and cellular neurobiology of this sensory system. Here, we describe methods for detection of olfactory receptor transcripts and proteins in the antennal olfactory organ of the fruit fly, Drosophila melanogaster, an important genetic model organism. We include protocols both for antennal cryosections and whole-mount antennae. These methods can be adapted for detection of receptor expression in other olfactory and gustatory tissues in Drosophila, as well as in the chemosensory systems of other insects.

Triiodothyronine and nerve growth factor are required to induce cytoplasmic dynein expression in rat dorsal root ganglion cultures.

Beside the several growth factors which play a crucial role in the development and regeneration of the nervous system, thyroid hormones also contribute to the normal development of the central and peripheral nervous system. In our previous work, we demonstrated that triiodothyronine (T3) in physiological concentration enhances neurite outgrowth of primary sensory neurons in cultures. Neurite outgrowth requires microtubules and microtubule associated proteins (MAPs). Therefore the effects of exogenous T3 or/and nerve growth factors (NGF) were tested on the expression of cytoskeletal proteins in primary sensory neurons. Dorsal root ganglia (DRG) from 19 day old rat embryos were cultured under four conditions: (1) control cultures in which explants were grown in the absence of T3 and NGF, (2) cultures grown in the presence of NGF alone, (3) in the presence of T3 alone or (4) in the presence of NGF and T3 together. Analysis of proteins by SDS-polyacrylamide gel electrophoresis revealed the presence of several proteins in the molecular weight region around 240 kDa. NGF and T3 together induced the expression of one protein, in particular, with a molecular weight above 240 kDa, which was identified by an antibody against MAP1c, a protein also known as cytoplasmic dynein. The immunocytochemical detection confirmed that this protein was expressed only in DRG explants grown in the presence of NGF and T3 together. Neither control explants nor explants treated with either NGF or T3 alone expressed dynein. In conclusion, a combination of nerve growth factor and thyroid hormone is necessary to regulate the expression of cytoplasmic dynein, a protein that is involved in retrograde axonal transport.

Constitutive expression of FGF4 disrupts the development of the eye and the anterior CNS during mouse embryogenesi, but does not influence the expressionof shh in these areas.

The biological consequences of constitutive fibroblast growth factor-4 (fgf4) expression have been analysed during anterior CNS development of mouse chimeric embryos. Severe mutant embryos exhibit exencephaly, absence of eye development and anomalous differentiation of neuropithelium. These embryos also show ectopic limb buds resembling the early phases of limb development. Because our results show that anterior CNS in those chimeric embrios does not express shh ectopically, we suggest that malformations may be due to interference between the ectopic expression of fgf4 in the cephalic area and the receptors for the members of the FGF family that regulate brain and eye development, namely fgf8. If this is correct, the results indirectly suport the crucial role of fgf8 in patterning the anterior CNS.

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.

Leishmania major: histone H1 gene expression from the sw3 locus.

Histone H1 in the parasitic protozoan Leishmania is a developmentally regulated protein encoded by the sw3 gene. Here we report that histone H1 variants exist in different Leishmania species and strains of L. major and that they are encoded by polymorphic genes. Amplification of the sw3 gene from the genome of three strains of L. major gave rise to different products in each strain, suggesting the presence of a multicopy gene family. In L. major, these genes were all restricted to a 50-kb Bg/II fragment found on a chromosomal band of 1.3 Mb (chromosome 27). The detection of RFLPs in this locus demonstrated its heterogeneity within several species and strains of Leishmania. Two different copies of sw3 (sw3.0 and sw3.1) were identified after screening a cosmid library containing L. major strain Friedlin genomic DNA. They were identical in their 5' UTRs and open reading frames, but differed in their 3' UTRs. With respect to the originally cloned copy of sw3 from L. major strain LV39, their open reading frames lacked a repeat unit of 9 amino acids. Immunoblots of L. guyanensis parasites transfected with these cosmids revealed that both copies could give rise to the histone H1 protein. The characterization of this locus will now make possible a detailed analysis of the function of histone H1 in Leishmania, as well as permit the dissection of the molecular mechanisms governing the developmental regulation of the sw3 gene.