Gene expression variation and expression quantitative trait mapping of human chromosome 21 genes.

Inter-individual differences in gene expression are likely to account for an important fraction of phenotypic differences, including susceptibility to common disorders. Recent studies have shown extensive variation in gene expression levels in humans and other organisms, and that a fraction of this variation is under genetic control. We investigated the patterns of gene expression variation in a 25 Mb region of human chromosome 21, which has been associated with many Down syndrome (DS) phenotypes. Taqman real-time PCR was used to measure expression variation of 41 genes in lymphoblastoid cells of 40 unrelated individuals. For 25 genes found to be differentially expressed, additional analysis was performed in 10 CEPH families to determine heritabilities and map loci harboring regulatory variation. Seventy-six percent of the differentially expressed genes had significant heritabilities, and genomewide linkage analysis led to the identification of significant eQTLs for nine genes. Most eQTLs were in trans, with the best result (P=7.46 x 10(-8)) obtained for TMEM1 on chromosome 12q24.33. A cis-eQTL identified for CCT8 was validated by performing an association study in 60 individuals from the HapMap project. SNP rs965951 located within CCT8 was found to be significantly associated with its expression levels (P=2.5 x 10(-5)) confirming cis-regulatory variation. The results of our study provide a representative view of expression variation of chromosome 21 genes, identify loci involved in their regulation and suggest that genes, for which expression differences are significantly larger than 1.5-fold in control samples, are unlikely to be involved in DS-phenotypes present in all affected individuals.

Tandem repeat sequence variation as causative cis-eQTLs for protein-coding gene expression variation: the case of CSTB.

Association studies have revealed expression quantitative trait loci (eQTLs) for a large number of genes. However, the causative variants that regulate gene expression levels are generally unknown. We hypothesized that copy-number variation of sequence repeats contribute to the expression variation of some genes. Our laboratory has previously identified that the rare expansion of a repeat c.-174CGGGGCGGGGCG in the promoter region of the CSTB gene causes a silencing of the gene, resulting in progressive myoclonus epilepsy. Here, we genotyped the repeat length and quantified CSTB expression by quantitative real-time polymerase chain reaction in 173 lymphoblastoid cell lines (LCLs) and fibroblast samples from the GenCord collection. The majority of alleles contain either two or three copies of this repeat. Independent analysis revealed that the c.-174CGGGGCGGGGCG repeat length is strongly associated with CSTB expression (P = 3.14 × 10(-11)) in LCLs only. Examination of both genotyped and imputed single-nucleotide polymorphisms (SNPs) within 2 Mb of CSTB revealed that the dodecamer repeat represents the strongest cis-eQTL for CSTB in LCLs. We conclude that the common two or three copy variation is likely the causative cis-eQTL for CSTB expression variation. More broadly, we propose that polymorphic tandem repeats may represent the causative variation of a fraction of cis-eQTLs in the genome.

Natural gene-expression variation in Down syndrome modulates the outcome of gene-dosage imbalance.

Down syndrome (DS) is characterized by extensive phenotypic variability, with most traits occurring in only a fraction of affected individuals. Substantial gene-expression variation is present among unaffected individuals, and this variation has a strong genetic component. Since DS is caused by genomic-dosage imbalance, we hypothesize that gene-expression variation of human chromosome 21 (HSA21) genes in individuals with DS has an impact on the phenotypic variability among affected individuals. We studied gene-expression variation in 14 lymphoblastoid and 17 fibroblast cell lines from individuals with DS and an equal number of controls. Gene expression was assayed using quantitative real-time polymerase chain reaction on 100 and 106 HSA21 genes and 23 and 26 non-HSA21 genes in lymphoblastoid and fibroblast cell lines, respectively. Surprisingly, only 39% and 62% of HSA21 genes in lymphoblastoid and fibroblast cells, respectively, showed a statistically significant difference between DS and normal samples, although the average up-regulation of HSA21 genes was close to the expected 1.5-fold in both cell types. Gene-expression variation in DS and normal samples was evaluated using the Kolmogorov-Smirnov test. According to the degree of overlap in expression levels, we classified all genes into 3 groups: (A) nonoverlapping, (B) partially overlapping, and (C) extensively overlapping expression distributions between normal and DS samples. We hypothesize that, in each cell type, group A genes are the most dosage sensitive and are most likely involved in the constant DS traits, group B genes might be involved in variable DS traits, and group C genes are not dosage sensitive and are least likely to participate in DS pathological phenotypes. This study provides the first extensive data set on HSA21 gene-expression variation in DS and underscores its role in modulating the outcome of gene-dosage imbalance.

Functional impact of genetic variation on gene expression

Numerous links between genetic variants and phenotypes are known and genome-wide association studies dramatically increased the number of genetic variants associated with traits during the last decade. However, how changes in the DNA perturb the molecular mechanisms and impact on the phenotype of an organism remains elusive. Studies suggest that many traitassociated variants are in the non-coding region of the genome and probably act through regulation of gene expression. During my thesis I investigated how genetic variants affect gene expression through gene regulatory mechanisms. The first chapter was a collaborative project with a pharmaceutical company, where we investigated genome-wide copy number variation (CNVs) among Cynomolgus monkeys (Macaca fascicularis) used in pharmaceutical studies, and associated them to changes in gene expression. We found substantial copy number variation and identified CNVs linked to tissue-specific expression changes of proximal genes. The second and third chapters focus on genetic variation in humans and its effects on gene regulatory mechanisms and gene expression. The second chapter studies two human trios, where the allelic effects of genetic variation on genome-wide gene expression, protein-DNA binding and chromatin modifications were investigated. We found abundant allele specific activity across all measured molecular phenotypes and show extended coordinated behavior among them. In the third chapter, we investigated the impact of genetic variation on these phenotypes in 47 unrelated individuals. We found that chromatin phenotypes are organized into local variable modules, often linked to genetic variation and gene expression. Our results suggest that chromatin variation emerges as a result of perturbations of cis-regulatory elements by genetic variants, leading to gene expression changes. The work of this thesis provides novel insights into how genetic variation impacts gene expression by perturbing regulatory mechanisms. -- De nombreux liens entre variations génétiques et phénotypes sont connus. Les études d'association pangénomique ont considérablement permis d'augmenter le nombre de variations génétiques associées à des phénotypes au cours de la dernière décennie. Cependant, comprendre comment ces changements perturbent les mécanismes moléculaires et affectent le phénotype d'un organisme nous échappe encore. Des études suggèrent que de nombreuses variations, associées à des phénotypes, sont situées dans les régions non codantes du génome et sont susceptibles d'agir en modifiant la régulation d'expression des gènes. Au cours de ma thèse, j'ai étudié comment les variations génétiques affectent les niveaux d'expression des gènes en perturbant les mécanismes de régulation de leur expression. Le travail présenté dans le premier chapitre est un projet en collaboration avec une société pharmaceutique. Nous avons étudié les variations en nombre de copies (CNV) présentes chez le macaque crabier (Macaca fascicularis) qui est utilisé dans les études pharmaceutiques, et nous les avons associées avec des changements d'expression des gènes. Nous avons découvert qu'il existe une variabilité substantielle du nombre de copies et nous avons identifié des CNVs liées aux changements d'expression des gènes situés dans leur voisinage. Ces associations sont présentes ou absentes de manière spécifique dans certains tissus. Les deuxième et troisième chapitres se concentrent sur les variations génétiques dans les populations humaines et leurs effets sur les mécanismes de régulation des gènes et leur expression. Le premier se penche sur deux trios humains, père, mère, enfant, au sein duquel nous avons étudié les effets alléliques des variations génétiques sur l'expression des gènes, les liaisons protéine-ADN et les modifications de la chromatine. Nous avons découvert que l'activité spécifique des allèles est abondante abonde dans tous ces phénotypes moléculaires et nous avons démontré que ces derniers ont un comportement coordonné entre eux. Dans le second, nous avons examiné l'impact des variations génétiques de ces phénotypes moléculaires chez 47 individus, sans lien de parenté. Nous avons observé que les phénotypes de la chromatine sont organisés en modules locaux, qui sont liés aux variations génétiques et à l'expression des gènes. Nos résultats suggèrent que la variabilité de la chromatine est due à des variations génétiques qui perturbent des éléments cis-régulateurs, et peut conduire à des changements dans l'expression des gènes. Le travail présenté dans cette thèse fournit de nouvelles pistes pour comprendre l'impact des différentes variations génétiques sur l'expression des gènes à travers les mécanismes de régulation.

Classification of human astrocytic gliomas on the basis of gene expression: a correlated group of genes with angiogenic activity emerges as a strong predictor of subtypes.

The development of targeted treatment strategies adapted to individual patients requires identification of the different tumor classes according to their biology and prognosis. We focus here on the molecular aspects underlying these differences, in terms of sets of genes that control pathogenesis of the different subtypes of astrocytic glioma. By performing cDNA-array analysis of 53 patient biopsies, comprising low-grade astrocytoma, secondary glioblastoma (respective recurrent high-grade tumors), and newly diagnosed primary glioblastoma, we demonstrate that human gliomas can be differentiated according to their gene expression. We found that low-grade astrocytoma have the most specific and similar expression profiles, whereas primary glioblastoma exhibit much larger variation between tumors. Secondary glioblastoma display features of both other groups. We identified several sets of genes with relatively highly correlated expression within groups that: (a). can be associated with specific biological functions; and (b). effectively differentiate tumor class. One prominent gene cluster discriminating primary versus nonprimary glioblastoma comprises mostly genes involved in angiogenesis, including VEGF fms-related tyrosine kinase 1 but also IGFBP2, that has not yet been directly linked to angiogenesis. In situ hybridization demonstrating coexpression of IGFBP2 and VEGF in pseudopalisading cells surrounding tumor necrosis provided further evidence for a possible involvement of IGFBP2 in angiogenesis. The separating groups of genes were found by the unsupervised coupled two-way clustering method, and their classification power was validated by a supervised construction of a nearly perfect glioma classifier.

The effect of translocation-induced nuclear reorganization on gene expression.

Translocations are known to affect the expression of genes at the breakpoints and, in the case of unbalanced translocations, alter the gene copy number. However, a comprehensive understanding of the functional impact of this class of variation is lacking. Here, we have studied the effect of balanced chromosomal rearrangements on gene expression by comparing the transcriptomes of cell lines from controls and individuals with the t(11;22)(q23;q11) translocation. The number of differentially expressed transcripts between translocation-carrying and control cohorts is significantly higher than that observed between control samples alone, suggesting that balanced rearrangements have a greater effect on gene expression than normal variation. Many of the affected genes are located along the length of the derived chromosome 11. We show that this chromosome is concomitantly altered in its spatial organization, occupying a more central position in the nucleus than its nonrearranged counterpart. Derivative 22-mapping chromosome 22 genes, on the other hand, remain in their usual environment. Our results are consistent with recent studies that experimentally altered nuclear organization, and indicated that nuclear position plays a functional role in regulating the expression of some genes in mammalian cells. Our study suggests that chromosomal translocations can result in hitherto unforeseen, large-scale changes in gene expression that are the consequence of alterations in normal chromosome territory positioning. This has consequences for the patterns of gene expression change seen during tumorigenesis-associated genome instability and during the karyotype changes that lead to speciation.

Disruption of gene expression in hybrids of the fire ants Solenopsis invicta and Solenopsis richteri.

Transcriptome analysis is a powerful tool for unveiling the distribution and magnitude of genetic incompatibilities between hybridizing taxa. The nature of such incompatibilities is closely associated with the evolutionary histories of the parental species and may differ across tissues and between the sexes. In eusocial insects, the presence of castes that experience divergent selection regimes may result in additional distinct patterns of caste-specific hybrid incompatibilities. We analysed levels of expression of >14 000 genes in two life stages of each caste in the fire ants Solenopsis invicta and Solenopsis richteri and in their hybrids. We found strong contributions of both developmental stage and caste to gene expression patterns. In contrast, variability in expression was only weakly associated with taxonomic identity, with hybrid scores falling between those of the two parental species. Hybrid incompatibilities were surprisingly modest, with only 32 genes being mis-expressed, indicating low levels of disruption in gene regulation in hybrids; males and workers each mis-expressed at least seven times as many genes as queens. Interestingly, homologues of many of the mis-expressed genes have been implicated in behavioural variation in Drosophila melanogaster. General expression profiles of hybrids consistently were more similar to those of S. richteri than S. invicta, presumably because S. richteri trans-regulatory elements tend to be dominant and/or because there is an overall bias in the genetic composition of the hybrids towards S. richteri. Altogether, our results suggest that selection acting on each caste may contribute differently to interspecific divergence and speciation in this group of ants.

Evolution at Two Levels in Fire Ants: The Relationship between Patterns of Gene Expression and Protein Sequence Evolution.

Variation in protein sequence and gene expression each contribute to phenotypic diversity, and may be subject to similar selective pressures. Eusocial insects are particularly useful for investigating the evolutionary link between protein sequence and condition-dependent patterns of gene expression because gene expression plays a central role in determining differences between eusocial insect sexes and castes. We investigated the relationship between protein coding sequence evolution and gene expression patterns in the fire ants Solenopsis invicta, S. richteri, and their hybrids to gain greater insight into how selection jointly operates on gene expression and coding sequence. We found that genes with high expression variability within castes and sexes were frequently differentially expressed between castes and sexes, as well as between species and hybrids. These results indicate that genes showing high variation in expression in one context also tend to show high variation in expression in other contexts. Our analyses further revealed that variation in both intra- and interspecific gene expression was positively associated with rate of protein sequence evolution in Solenopsis. This suggests that selective constraints on a gene operate both at the level of protein sequence and at the level of gene expression regulation. Overall, our study provides one of the strongest demonstrations that selective constraints mediate both protein sequence evolution and gene expression variability across different biological contexts and timescales.

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.

Copy number variants, diseases and gene expression.

Copy number variation (CNV) has recently gained considerable interest as a source of genetic variation likely to play a role in phenotypic diversity and evolution. Much effort has been put into the identification and mapping of regions that vary in copy number among seemingly normal individuals in humans and a number of model organisms, using bioinformatics or hybridization-based methods. These have allowed uncovering associations between copy number changes and complex diseases in whole-genome association studies, as well as identify new genomic disorders. At the genome-wide scale, however, the functional impact of CNV remains poorly studied. Here we review the current catalogs of CNVs, their association with diseases and how they link genotype and phenotype. We describe initial evidence which revealed that genes in CNV regions are expressed at lower and more variable levels than genes mapping elsewhere, and also that CNV not only affects the expression of genes varying in copy number, but also have a global influence on the transcriptome. Further studies are warranted for complete cataloguing and fine mapping of CNVs, as well as to elucidate the different mechanisms by which they influence gene expression.

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

Trade-off between constitutive and inducible resistance against herbivores is only partially explained by gene expression and glucosinolate production.

The hypothesis that constitutive and inducible plant resistance against herbivores should trade-off because they use the same resources and impose costs to plant fitness has been postulated for a long time. Negative correlations between modes of deployment of resistance and defences have been observed across and within species in common garden experiments. It was therefore tested whether that pattern of resistance across genotypes follows a similar variation in patterns of gene expression and chemical defence production. Using the genetically tractable model Arabidopsis thaliana and different modes of induction, including the generalist herbivore Spodoptera littoralis, the specialist herbivore Pieris brassicae, and jasmonate application, constitutive and inducibility of resistance was measured across seven A. thaliana accessions that were previously selected based on constitutive levels of defence gene expression. According to theory, it was found that modes of resistance traded-off among accessions, particularly against S. littoralis, in which accessions investing in high constitutive resistance did not increase it substantially after attack and vice-versa. Accordingly, the average expression of eight genes involved in glucosinolate production negatively predicted larval growth across the seven accessions. Glucosinolate production and genes related to defence induction on healthy and herbivore-damaged plants were measured next. Surprisingly, only a partial correlation between glucosinolate production, gene expression, and the herbivore resistance results was found. These results suggest that the defence outcome of plants against herbivores goes beyond individual molecules or genes but stands on a complex network of interactions.

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.

Post-transcriptional regulation of circadian gene expression by miRNAs

La grande majorité des organismes vivants ont développé un système d'horloges biologiques internes, appelées aussi horloges circadiennes, contrôlant l'expression de gênes impliqués dans de nombreux processus moléculaires et comportementaux. Au cours de la dernière décennie, des analyses « microarray » et séquençages à haut débit sur divers tissus de mammifères, indiquent que jusqu'à 20% du transcriptome serait sous contrôle circadien. Il était jusqu'à présent admis que la majorité des ARNm ayant une accumulation rythmique était générée par une transcription qui était elle-même rythmique. Toutefois, de récentes études ont suggéré qu'une proportion considérable des ARNm cycliques serait en fait générée par des mécanismes post-transcriptionnelles, incluant une régulation par micro-ARN (miARN). Lorsque j'ai débuté mon travail de thèse, l'influence des miARN sur l'expression des gènes circadiens, au niveau pangénomique, était encore méconnue. Par l'utilisation d'un modèle murin, dont la biogenèse des miARN a été spécifiquement désactivée au niveau des cellules hépatiques (knockout conditionnel pour Dicer), je me suis donc intéressée au rôle que jouaient ces molécules régulatrices sur la rythmicité de l'expression génique dans le foie. Des séquençages sur l'ensemble du transcriptome révèlent que l'horloge interne du foie est étonnement résistante à la perte totale des miARN. Nous avons cependant trouvé que les miARN agissent de façon importante sur la régulation de l'expression des gènes contrôlés par l'horloge moléculaire. La corégulation par les miARN, affectant jusqu'à 30% des gènes transcrits de façon rythmiques, conduit ainsi à une modulation de phase et d'amplitude du rythme de l'abondance des ARNm. En revanche, seuls peu de transcrits dépendent uniquement des miARN pour la rythmicité de leur accumulation. Enfin, mon travail met en évidence plusieurs miARN spécifiques, qui semblent préférentiellement moduler l'expression des gènes cycliques et permet l'identification de voies hépatiques particulièrement sujettes à une double régulation par les miARN et l'horloge biologique interne. La première masse d'analyses a essentiellement porté sur le rôle que jouent les miARN au niveau de l'expression des gènes contrôlés par l'horloge interne. Dans deux études de suivi, je me suis penchée sur deux aspects supplémentaires et complémentaires de la manière dont les miARN et l'oscillation de l'expression des gènes interagissent. Dans les hépatocytes murins, spécifiquement privés de Dicer, je me suis demandée si un phénotype horloge avait pu être masqué, dû à un entraînement stable de l'horloge du foie par l'horloge maîtresse du cerveau. J'ai donc commencé une série d'expériences ambitieuses (impliquant la mesure de la rythmicité du foie in vivo, chez l'animal vivant) afin de déséquilibrer l'entrainement de l'horloge hépatique via l'utilisation d'un protocole nutritionnel spécifique. Les premiers résultats suggèrent que dans des conditions où l'animal subit une restriction alimentaire pendant la journée, les miARN sont importants dans la cinétique d'adaptation des organes périphériques à un nouvel horaire de sustentation. Dans une deuxième ligne de recherche, j'ai plus profondément étudié quels seraient les miARN responsables des rythmes post-transcriptionnels des ARNm, en utilisant le séquençage de « small » ARN sur 24h. L'analyse est en cours et se poursuivra après l'obtention de mon diplôme. De façon générale, mon travail révèle d'importants et nouveaux rôles des miARN dans la modulation de l'expression circadienne des gènes hépatiques. De plus, le set de données générées dans l'étude déjà publiée, peut dorénavant servir de ressource valable pour de prochaines investigations sur le rôle physiologique que les miARN jouent au niveau du foie. -- Most living organisms have developed internal timing systems, called circadian clocks, to drive the rhythmic expression of genes involved in many molecular and behavioral processes. Over the last decade, microarray analyses and high- throughput sequencing from various mammalian tissues have indicated that up to 20% of the transcriptome are under circadian control. It was generally assumed that the majority of rhythmic mRNA accumulation is generated by rhythmic transcription. However, recent studies have suggested that a considerable proportion of mRNA cycling may actually be generated by post-transcriptional mechanisms, including by microRNAs. When I started my thesis work, it was still unknown how miRNAs influence circadian gene expression in a genome-wide fashion. Using a mouse model in which miRNA biogenesis can be inactivated in hepatocytes (conditional Dicer knockout mouse), I have thus addressed the role that these regulatory molecules play in rhythmic gene expression in the liver. Whole transcriptome sequencing revealed that the hepatic core clock was surprisingly resilient to total miRNA loss. However, we found that miRNAs acted as important regulators of clock-controlled gene expression. Co- regulation by miRNAs, which affected up to 30% of rhythmically transcribed genes, thus led to the modulation of phases and amplitudes of mRNA abundance rhythms. By contrast, only very few transcripts were strictly dependent on miRNAs for their rhythmic accumulation. Finally, my work highlights several specific miRNAs that appear to preferentially modulate cyclic gene expression, and identifies pathways in the liver that are particularly prone to dual regulation through miRNAs and the clock. The first bulk of analyses mainly dealt with the role that miRNAs play at the level of rhythmic clock output gene expression. In two follow-up studies I further delved into two additional, complementary aspects of how miRNAs and gene expression oscillations interact. First, I addressed whether a core clock phenotype in the hepatocyte-specific Dicer knockout could have been masked due to the stable entrainment of the liver clock by the animals' master clock in the brain. I thus started a series of ambitious experiments (involving the in vivo recording of liver rhythms in live animals) to bring the stable entrainment of the liver clock out of equilibrium using specific feeding protocols. My first results suggest that under conditions when animals are challenged by food restriction to daytime, miRNAs are important for the kinetics of adapting to unusual mealtime in peripheral tissue. In a second line of research, I have more carefully investigated which miRNAs are responsible for post- transcriptional mRNA rhythms using small RNA sequencing around-the-clock. The analyses are ongoing and will be continued after my graduation. Overall, my work uncovered important and novel roles of miRNA activity in shaping hepatic circadian gene expression; moreover, the datasets collect in the published studies can serve as a valuable resource for further investigations into the physiological roles that miRNAs play in liver. -- L'alternance du jour et de la nuit dirige depuis longtemps la vie quotidienne des êtres humains et de la plupart des organismes sur terre. Ce cycle de 24 heures façonne beaucoup de changements comportementaux et physiologiques tels que la vigilance, la température corporelle et le sommeil. Les rythmes journaliers, appelés rythmes circadiens, sont dirigés par des horloges biologiques tournant dans presque chaque cellule du corps. Une structure dans le cerveau agit en tant qu'horloge maitresse pour synchroniser les horloges internes entre elles et en fonction des signaux de jour/nuit extérieurs. Dans les cellules "les gènes de l'horloge" sont activés et désactivés une fois par jour ce qui déclenche des cycles dans lesquels des protéines sont produites de manière circadienne. Ces rythmes protéiques sont spécialisés pour chaque tissu ou organe et peuvent les aider à réaliser leurs tâches quotidiennes. Les rythmes circadiens peuvent être générés d'autres manières n'impliquant pas directement les composants des gènes de l'horloge. Les ARN messagers (ARNm) sont des molécules intermédiaires dans la production de protéines à partir d'ADN. Dans le foie des souris jusqu'à 20% des molécules d'ARNm sont produites suivant des rythmes circadiens. Le foie réalise des tâches essentielles dans le contrôle du métabolisme incluant celui des hydrates de carbone, des graisses et du cholestérol. Un timing précis est important afin de traiter les substances nutritives correctement lors des repas il en résulte une variation des quantités de certains ARNm et protéines coïncidant avec les repas. Les microARNs constituent une autre classe de molécules ARN de très petite taille qui régulent l'efficacité de traduction des ARNm en protéines et la stabilité des ARNm. Lors de mon travail de thèse, j'ai exploré de manière approfondie l'influence de ces petits régulateurs sur les rythmes circadiens du foie de souris. Ces expériences qui impliquaient le "Knock-out" d'un gène essentiel à la production de microARNs montrent qu'au lieu de générer les rythmes des ARNm, les microARNs les ajustent pour répondre aux besoins spécifiques du foie comme assurer leur pic au bon moment de la journée. Le ciblage de microARNs spécifiques peut révéler de nouvelles stratégies pour rectifier ces rythmes lorsque par exemple les fonctions métaboliques ne fonctionnent plus normalement. -- The rising and setting of the sun have long driven the daily schedules of humans and most organisms on the earth. This 24-hr cycle shapes many behavioural and physiological changes, such as alertness, body temperature, and sleep. These daily rhythms, which are called circadian rhythms, are dictated by biological clocks that are ticking in almost every single cell of the body. A region in the brain acts as a master clock to synchronize the internal clocks with each other and with the outside light/dark cycles. In cells, "core clock genes" are turned on and off once per day, which triggers cycles that cause some proteins to be produced in a circadian manner. The protein rhythms are specialized to a particular tissue or organ, and may help them to carry out their designated daily tasks. However, circadian rhythms might also be produced by other ways that do not involve these core clock components. Messenger RNAs (mRNAs) are intermediate molecules in the production of proteins from DNA. In the mouse liver, up to 20% of mRNA molecules are produced in circadian cycles. The liver performs essential tasks that control metabolism-including that of carbohydrates, fats, and cholesterol. Precisely timing when certain mRNAs and proteins reach peaks and troughs in their activities to coincide with mealtimes is important for nutrients to be properly processed. Other RNA molecules called microRNAs, i.e. RNAs of very small size, regulate at which rate mRNA molecules are translated into proteins. In my thesis work, I have explored at the influence of these small regulators on circadian rhythms in the mouse liver in greater detail. These experiments, which involved "knocking out" a gene that is essential for the production of microRNAs, show that rather than generating the mRNA rhythms, the microRNAs appear to adjust them to meet the specific needs of the liver, such as ensuring that they peak at the right time-of-day. Targeting specific microRNA molecules may reveal new strategies to tweak these rhythms, which could help to improve conditions when metabolic functions go wrong.

Changes in reproductive roles are associated with changes in gene expression in fire ant queens.

Abstract In species with social hierarchies, the death of dominant individuals typically upheaves the social hierarchy and provides an opportunity for subordinate individuals to become reproductives. Such a phenomenon occurs in the monogyne form of the fire ant, Solenopsis invicta, where colonies typically contain a single wingless reproductive queen, thousands of workers and hundreds of winged nonreproductive virgin queens. Upon the death of the mother queen, many virgin queens shed their wings and initiate reproductive development instead of departing on a mating flight. Workers progressively execute almost all of them over the following weeks. To identify the molecular changes that occur in virgin queens as they perceive the loss of their mother queen and begin to compete for reproductive dominance, we collected virgin queens before the loss of their mother queen, 6 h after orphaning and 24 h after orphaning. Their RNA was extracted and hybridized against microarrays to examine the expression levels of approximately 10 000 genes. We identified 297 genes that were consistently differentially expressed after orphaning. These include genes that are putatively involved in the signalling and onset of reproductive development, as well as genes underlying major physiological changes in the young queens.