Genome-wide location analysis and expression studies reveal a role for p110 CUX1 in the activation of DNA replication genes.

Proteolytic processing of the CUX1 transcription factor generates an isoform, p110 that accelerates entry into S phase. To identify targets of p110 CUX1 that are involved in cell cycle progression, we performed genome-wide location analysis using a promoter microarray. Since there are no antibodies that specifically recognize p110, but not the full-length protein, we expressed physiological levels of a p110 isoform with two tags and purified chromatin by tandem affinity purification (ChAP). Conventional ChIP performed on synchronized populations of cells confirmed that p110 CUX1 is recruited to the promoter of cell cycle-related targets preferentially during S phase. Multiple approaches including silencing RNA (siRNA), transient infection with retroviral vectors, constitutive expression and reporter assays demonstrated that most cell cycle targets are activated whereas a few are repressed or not affected by p110 CUX1. Functional classes that were over-represented among targets included DNA replication initiation. Consistent with this finding, constitutive expression of p110 CUX1 led to a premature and more robust induction of replication genes during cell cycle progression, and stimulated the long-term replication of a plasmid bearing the oriP replicator of Epstein Barr virus (EBV).

Genome wide search for genetic determinants of habitual alcohol, tobacco and coffee use, obesity-related traits, response to mental and physical stress and hemodynamic traits

Les habitudes de consommation de substances psychoactives, le stress, l’obésité et les traits cardiovasculaires associés seraient en partie reliés aux mêmes facteurs génétiques. Afin d’explorer cette hypothèse, nous avons effectué, chez 119 familles multi-générationnelles québécoises de la région du Saguenay-Lac-St-Jean, des études d’association et de liaison pangénomiques pour les composantes génétiques : de la consommation usuelle d’alcool, de tabac et de café, de la réponse au stress physique et psychologique, des traits anthropométriques reliés à l’obésité, ainsi que des mesures du rythme cardiaque (RC) et de la pression artérielle (PA). 58000 SNPs et 437 marqueurs microsatellites ont été utilisés et l’annotation fonctionnelle des gènes candidats identifiés a ensuite été réalisée. Nous avons détecté des corrélations phénotypiques significatives entre les substances psychoactives, le stress, l’obésité et les traits hémodynamiques. Par exemple, les consommateurs d’alcool et de tabac ont montré un RC significativement diminué en réponse au stress psychologique. De plus, les consommateurs de tabac avaient des PA plus basses que les non-consommateurs. Aussi, les hypertendus présentaient des RC et PA systoliques accrus en réponse au stress psychologique et un indice de masse corporelle (IMC) élevé, comparativement aux normotendus. D’autre part, l’utilisation de tabac augmenterait les taux corporels d’épinéphrine, et des niveaux élevés d’épinéphrine ont été associés à des IMC diminués. Ainsi, en accord avec les corrélations inter-phénotypiques, nous avons identifié plusieurs gènes associés/liés à la consommation de substances psychoactives, à la réponse au stress physique et psychologique, aux traits reliés à l’obésité et aux traits hémodynamiques incluant CAMK4, CNTN4, DLG2, DAG1, FHIT, GRID2, ITPR2, NOVA1, NRG3 et PRKCE. Ces gènes codent pour des protéines constituant un réseau d’interactions, impliquées dans la plasticité synaptique, et hautement exprimées dans le cerveau et ses tissus associés. De plus, l’analyse des sentiers de signalisation pour les gènes identifiés (P = 0,03) a révélé une induction de mécanismes de Potentialisation à Long Terme. Les variations des traits étudiés seraient en grande partie liées au sexe et au statut d’hypertension. Pour la consommation de tabac, nous avons noté que le degré et le sens des corrélations avec l’obésité, les traits hémodynamiques et le stress sont spécifiques au sexe et à la pression artérielle. Par exemple, si des variations ont été détectées entre les hommes fumeurs et non-fumeurs (anciens et jamais), aucune différence n’a été observée chez les femmes. Nous avons aussi identifié de nombreux traits reliés à l’obésité dont la corrélation avec la consommation de tabac apparaît essentiellement plus liée à des facteurs génétiques qu’au fait de fumer en lui-même. Pour le sexe et l’hypertension, des différences dans l’héritabilité de nombreux traits ont également été observées. En effet, des analyses génétiques sur des sous-groupes spécifiques ont révélé des gènes additionnels partageant des fonctions synaptiques : CAMK4, CNTN5, DNM3, KCNAB1 (spécifique à l’hypertension), CNTN4, DNM3, FHIT, ITPR1 and NRXN3 (spécifique au sexe). Ces gènes codent pour des protéines interagissant avec les protéines de gènes détectés dans l’analyse générale. De plus, pour les gènes des sous-groupes, les résultats des analyses des sentiers de signalisation et des profils d’expression des gènes ont montré des caractéristiques similaires à celles de l’analyse générale. La convergence substantielle entre les déterminants génétiques des substances psychoactives, du stress, de l’obésité et des traits hémodynamiques soutiennent la notion selon laquelle les variations génétiques des voies de plasticité synaptique constitueraient une interface commune avec les différences génétiques liées au sexe et à l’hypertension. Nous pensons, également, que la plasticité synaptique interviendrait dans de nombreux phénotypes complexes influencés par le mode de vie. En définitive, ces résultats indiquent que des approches basées sur des sous-groupes et des réseaux amélioreraient la compréhension de la nature polygénique des phénotypes complexes, et des processus moléculaires communs qui les définissent.

Molecular determinants of congenital hypothyroidism due to thyroid dysgenesis

L’hypothyroïdie congénitale par dysgénésie thyroïdienne (HCDT) est la condition endocrinienne néonatale la plus fréquemment rencontrée, avec une incidence d’un cas sur 4000 naissances vivantes. L’HCDT comprend toutes les anomalies du développement de la thyroïde. Parmi ces anomalies, le diagnostic le plus fréquent est l’ectopie thyroïdienne (~ 50% des cas). L’HCDT est fréquemment associée à un déficit sévère en hormones thyroïdiennes (hypothyroïdisme) pouvant conduire à un retard mental sévère si non traitée. Le programme de dépistage néonatal assure un diagnostic et un traitement précoce par hormones thyroïdiennes. Cependant, même avec un traitement précoce (en moyenne à 9 jours de vie), un retard de développement est toujours observé, surtout dans les cas les plus sévères (c.-à-d., perte de 10 points de QI). Bien que des cas familiaux soient rapportés (2% des cas), l’HCTD est essentiellement considérée comme une entité sporadique. De plus, plus de 92% des jumeaux monozygotiques sont discordants pour les dysgénésies thyroïdiennes et une prédominance féminine est rapportée (spécialement dans le cas d’ectopies thyroïdiennes), ces deux observations étant clairement incompatible avec un mode de transmission héréditaire mendélien. Il est donc cohérent de constater que des mutations germinales dans les facteurs de transcription thyroïdiens connus (NKX2.1, PAX8, FOXE1, and NKX2.5) ont été identifiées dans seulement 3% des cas sporadiques testés et furent, de plus, exclues lors d’analyse d’association dans certaines familles multiplex. Collectivement, ces données suggèrent que des mécanismes non mendéliens sont à l’origine de la majorité des cas de dysgénésie thyroïdienne. Parmi ces mécanismes, nous devons considérer des modifications épigénétiques, des mutations somatiques précoces (au stade du bourgeon thyroïdien lors des premiers stades de l’embryogenèse) ou des défauts développementaux stochastiques (c.-à-d., accumulation aléatoire de mutations germinales ou somatiques). Voilà pourquoi nous proposons un modèle «2 hits » combinant des mutations (épi)génétiques germinales et somatiques; ce modèle étant compatible avec le manque de transmission familial observé dans la majorité des cas d’HCDT. Dans cette thèse, nous avons déterminé si des variations somatiques (épi)génétiques sont associées à l’HCTD via une approche génomique et une approche gène candidat. Notre approche génomique a révélé que les thyroïdes ectopiques ont un profil d’expression différent des thyroïdes eutopiques (contrôles) et que ce profil d’expression est enrichi en gènes de la voie de signalisation Wnt. La voie des Wnt est cruciale pour la migration cellulaire et pour le développement de plusieurs organes dérivés de l’endoderme (p.ex. le pancréas). De plus, le rôle de la voie des Wnt dans la morphogénèse thyroïdienne est supporté par de récentes études sur le poisson-zèbre qui montrent des anomalies du développement thyroïdien lors de la perturbation de la voie des Wnt durant différentes étapes de l’organogénèse. Par conséquent, l’implication de la voie des Wnt dans l’étiologie de la dysgénésie thyroïdienne est biologiquement plausible. Une trouvaille inattendue de notre approche génomique fut de constater que la calcitonine était exprimée autant dans les thyroïdes ectopiques que dans les thyroïdes eutopiques (contrôles). Cette trouvaille remet en doute un dogme de l’embryologie de la thyroïde voulant que les cellules sécrétant la calcitonine (cellules C) proviennent exclusivement d’une structure extrathyroïdienne (les corps ultimobranchiaux) fusionnant seulement avec la thyroïde en fin de développement, lorsque la thyroïde a atteint son emplacement anatomique définitif. Notre approche gène candidat ne démontra aucune différence épigénétique (c.-à-d. de profil de méthylation) entre thyroïdes ectopiques et eutopiques, mais elle révéla la présence d’une région différentiellement méthylée (RDM) entre thyroïdes et leucocytes dans le promoteur de FOXE1. Le rôle crucial de FOXE1 dans la migration thyroïdienne lors du développement est connu et démontré dans le modèle murin. Nous avons démontré in vivo et in vitro que le statut de méthylation de cette RDM est corrélé avec l’expression de FOXE1 dans les tissus non tumoraux (c.-à-d., thyroïdes et leucocytes). Fort de ces résultats et sachant que les RDMs sont de potentiels points chauds de variations (épi)génétiques, nous avons lancé une étude cas-contrôles afin de déterminer si des variants génétiques rares localisés dans cette RDM sont associés à la dysgénésie thyroïdienne. Tous ces résultats générés lors de mes études doctorales ont dévoilé de nouveaux mécanismes pouvant expliquer la pathogenèse de la dysgénésie thyroïdienne, condition dont l’étiologie reste toujours une énigme. Ces résultats ouvrent aussi plusieurs champs de recherche prometteurs et vont aider à mieux comprendre tant les causes des dysgénésies thyroïdiennes que le développement embryonnaire normal de la thyroïde chez l’homme.

Évaluation neurobiologique des souris spontanément hypertendues : Du vieillissement à la génomique comparative

Le but de cette thèse est premièrement d’évaluer l’effet du vieillissement sur les fonctions psychomotrices des souches de souris sélectionnées génétiquement en fonction de leur tension artérielle (TA); deuxièmement, de localiser les déterminants génétiques des phénotypes psychophysiologiques à partir de souches recombinantes congéniques (RCS). Ces travaux ont mené à la publication de 4 articles. Le premier article décrit l’évaluation des fonctions psychomotrices des souches avec une tension artérielle élevée (HBP), basse (LBP) et normale (NBP). La performance aux épreuves d’exploration, d’habiletés motrices et d’apprentissage spatial, a été mesurée sur deux cohortes âgées respectivement de 12 mois et de trois mois. Indépendamment de l’âge, les HBPs sont hyperactives dans l’open-field (OF), mais pas dans le test d’exploration de trous. Inversement, les LBP explorent moins d’espaces que les NBP et, à trois mois seulement, sont hypoactives dans l’OF. Par ailleurs, les HBPs et les LBP présentent des déficits précoces de coordination motrice et des fonctions visuo-motrices. Le second article concerne l’évaluation longitudinale de la coordination motrice, de l’anxiété et de l’apprentissage spatial des souches HBP, LBP et NBP, à l’âge de deux mois et de 12 mois. Le vieillissement accentue l’hyperactivité des HBPs dans l’OF. Par contre, l’hypoactivité des souris LBP est détectable seulement à l’âge de deux mois. Indépendamment de l’âge, les souris HBP et LBP montrent une perception réduite du danger dans l’épreuve d’anxiété et des dysfonctions visuo-motrices au labyrinthe aquatique. Enfin, des déficits précoces de coordination motrice se manifestent seulement chez les HBPs. Il reste à déterminer si les déficits observés sont liés à des déterminants génétiques indépendants ou secondaires aux altérations de la tension artérielle. Le troisième article présente la comparaison entre les souches consanguines A/J et C57Bl/6J (B6) aux épreuves de l’OF, de la planche à trous, du labyrinthe aquatique et du cintre (coordination motrice). Les B6 explore d’avantage l’OF et la planche à trous. Les B6 sont moins rapides sur le cintre, mais supérieurs aux A/J dans le labyrinthe aquatique, avec une plate-forme invisible ou visible. Ces résultats démontrent l’implication de déterminants génétiques. Cette thèse se termine par un quatrième article sur la localisation des déterminants génétiques de la susceptibilité au stress dans les RCS, dérivées de A/J et B6, et présentant un agencement spécifique de 12.5% du génome. La réactivité émotionnelle est évaluée dans l’OF et le plus-maze; la réponse de stress est mesurée par radio télémétrie de la température interne pendant le stress d’immobilisation (SI) sous diète régulière et riche en sel; l’excrétion des électrolytes urinaires est dosée après 24 heures de diète salée. Les loci les plus significatifs sont situés dans les régions suivantes: de l’émotionalité dans l’OF (Emo1) sur le chr. 1 (LOD=4.6) correspondant à la région homologue impliquée dans la cohorte d’hypertension familiale du Saguenay; de la dopa décarboxylase (ddc) sur le chr. 11 pour l’émergence du plus-maze (LOD=4.7); de la protéine liant l’endotoxine (lbp) sur le chr. 2 pour l’hypothermie initiale en réponse au SI (LOD=4); et de HSP90 sur le chr. 12 pour l’excrétion de Ca++ (LOD=4.6). Des banques de données sont ensuite interrogées pour recenser les polymorphismes des régions régulatrices ou codantes des gènes candidats chez les souches ancestrales A/J et B6, dont les séquences sont disponibles pour le génome entier. Des utilitaires web permettent de dévoiler les changements dans la structure secondaire de l’ARNm, l’interférence avec des microARN ou avec d’autres motifs de liaison. Plusieurs SNPs fonctionnels ont été identifiés pour le QTL du chr. 1, particulièrement dans les éléments de régulation; ceux-ci impliquant des gènes reliés avec les réponses inflammatoire/immunitaire ou avec le système cardiovasculaire. La quantification par la PCR confirme une régulation à la baisse d’atp1a2 dans le cœur et le cerveau des souches susceptibles à l’anxiété. Ces résultats confirment l’intrication des altérations de la susceptibilité au stress et de la régulation de la TA.

Caractérisation clinique et génétique d’une nouvelle forme d’ataxie autosomique récessive dans la population québécoise

Les ataxies autosomiques récessives sont un groupe de troubles neurologiques hétérogènes caractérisés par une incoordination brute des mouvements musculaires impliquant le dysfonctionnement nerveux du cervelet qui coordonne le mouvement. Plusieurs formes héréditaires ont été décrites dont la plus connue : l’ataxie de Friedriech. Dans cette thèse nous rapportons l'identification et la caractérisation d’une nouvelle forme dans la population québécoise. L’ataxie récessive spastique avec leucoencéphalopathie (ARSAL; aussi connue comme l’ataxie autosomique récessive spastique de type 3 (SPAX3); OMIM 611390) est la deuxième ataxie spastique décrite dans la population canadienne française. En effet, près de 50 % de nos cas sont originaires de la région de Portneuf. En 2006, nous avons décrit les caractéristiques cliniques de cette nouvelle forme d’ataxie. Un premier criblage du génome entier, constitué de plus de 500 marqueurs microsatellites, a permis la localisation du locus sur le chromosome 2q33-34. Suite au séquençage de plus de 37 gènes candidats et afin de rétrécir cet intervalle candidat, nous avons utilisé une micro-puce d’ADN constituée de marqueurs SNP «single nucleotide polymorphism» et nous avons identifié un deuxième intervalle candidat de 0.658Mb au locus 2q33 dans lequel se trouvent moins de 9 gènes. L’identification et la caractérisation de ces mutations a nécessité l’utilisation de diverses technologies de pointe. Trois mutations (une délétion et deux réarrangements complexes) dans le gène mitochondrial tRNA-synthetase (MARS2) ont été identifiées dans notre cohorte. Nous émettons l’hypothèse que la nature des mutations complexes est responsable d’un dérèglement de la transcription du gène, ce qui a un impact néfaste sur la fonction mitochondriale et le tissu neuronal.

Recherche des facteurs génétiques à l’origine de la maladie de Parkinson dans la population canadienne-française du Québec

La maladie de Parkinson (MP) est une affection neurodégénérative invalidante et incurable. Il est maintenant clairement établi que d’importants déterminants génétiques prédisposent à son apparition. La recherche génétique sur des formes familiales de la MP a mené à la découverte d’un minimum de six gènes causatifs (SNCA, LRRK2, Parkin, PINK1, DJ-1 and GBA) et certains, par exemple LRRK2, contiennent des variations génétiques qui prédisposent également aux formes sporadiques. La caractérisation des protéines codées par ces gènes a mené à une meilleure compréhension des mécanismes moléculaires sousjacents. Toutefois, en dépit de ces efforts, les causes menant à l’apparition de la MP restent inconnues pour la majorité des patients. L’objectif général des présents travaux était d’identifier des mutations prédisposant à la MP dans la population canadienne-française du Québec à partir d’une cohorte composée principalement de patients sporadiques. Le premier volet de ce projet consistait à déterminer la présence de mutations de LRRK2 dans notre cohorte en séquençant directement les exons contenant la majorité des mutations pathogéniques et en effectuant une étude d’association. Nous n’avons identifié aucune mutation et l’étude d’association s’est avérée négative, suggérant ainsi que LRRK2 n’est pas une cause significative de la MP dans la population canadienne-française. La deuxième partie du projet avait pour objectif d’identifier de nouveaux gènes causatifs en séquençant directement des gènes candidats choisis à cause de leurs implications dans différents mécanismes moléculaires sous-tendant la MP. Notre hypothèse de recherche était basée sur l’idée que la MP est principalement due à des mutations individuellement rares dans un grand nombre de gènes différents. Nous avons identifié des mutations rares dans les gènes PICK1 et MFN1. Le premier code pour une protéine impliquée dans la régulation de la transmission du glutamate tandis que le second est un des acteurs-clés du processus de fusion mitochondriale. Nos résultats, qui devront être répliqués, suggèrent que le séquençage à grande échelle pourrait être une méthode prometteuse d’élucidation des facteurs de prédisposition génétiques à la MP ; ils soulignent l’intérêt d’utiliser une population fondatrice comme les canadiens-français pour ce type d’étude et devraient permettre d’approfondir les connaissances sur la pathogénèse moléculaire de la MP.

Génomique fonctionnelle des cellules corticotropes hypophysaires : contrôle génétique de la gestion systémique des stress

L'axe hypothalamo-hypophyso-surrénalien (HPA) permet de maintenir l'homéostasie de l'organisme face à divers stress. Qu'ils soient de nature psychologique, physique ou inflammatoire/infectieux, les stress provoquent la synthèse et la libération de CRH par l'hypothalamus. Les cellules corticotropes hypophysaires perçoivent ce signal et en réaction, produisent et sécrètent l'ACTH. Ceci induit la synthèse des glucocorticoïdes (Gc) par le cortex surrénalien; ces stéroïdes mettent le système métabolique en état d’alerte pour la réponse au stress et à l’agression. Les Gc ont le rôle essentiel de contrôler les défenses de l'organisme, en plus d'exercer une rétro-inhibition sur l'axe HPA. L'ACTH est une petite hormone peptidique produite par le clivage d'un précurseur: la pro-opiomélanocortine (POMC). À cause de sa position critique dans la normalisation de l'homéostasie, le contrôle transcriptionnel du gène Pomc a fait l'objet d'études approfondies au cours des dernières décennies. Nous savons maintenant que la région promotrice du gène Pomc permet une expression ciblée dans les cellules POMC hypophysaires. L'étude du locus Pomc par des technologies génomiques m'a permis de découvrir un nouvel élément de régulation qui est conservé à travers l'évolution des mammifères. La caractérisation de cet enhancer a démontré qu'il dirige une expression restreinte à l'hypophyse, et plus particulièrement dans les cellules corticotropes. De façon intéressante, l'activité de cet élément dépend d'un nouveau site de liaison recrutant un homodimère du facteur de transcription Tpit, dont l'expression est également limitée aux cellules POMC de l'hypophyse. La découverte de cet enhancer ajoute une toute nouvelle dimension à la régulation de l'expression de POMC. Les cytokines pro-inflammatoires IL6/LIF et les Gc sont connus pour leur antagonisme sur la réaction inflammatoire et sur le promoteur Pomc via l'action des facteurs de transcription Stat3 et GR respectivement. L'analyse génomique des sites liés ii par ces deux facteurs nous a révélé une interrelation complexe et a permis de définir un code transcriptionnel entre ces voies de signalisation. En plus de leur action par interaction directe avec l’ADN au niveau des séquences régulatrices, ces facteurs interagissent directement entre eux avec des résultats transcriptionnels différents. Ainsi, le recrutement de GR par contact protéine:protéine (tethering) sur Stat3 étant lié à l'ADN provoque un antagonisme transcriptionnel. Inversement, le tethering de Stat3 sur GR supporte une action synergique, tout comme leur co-recrutement à l'ADN sur des sites contigus ou composites. Lors d'une activation soutenue, ce synergisme entre les voies IL6/LIF et Gc induit une réponse innée de défense cellulaire. Ainsi lors d'un stress majeur, ce mécanisme de défense est mis en branle dans toutes les cellules et tissus. En somme, les travaux présentés dans cette thèse définissent les mécanismes transcriptionnels engagés dans le combat de l'organisme contre les stress. Plus particulièrement, ces mécanismes ont été décrits au niveau de la réponse globale des corticotropes et du gène Pomc. Il est essentiel pour l'organisme d'induire adéquatement ces mécanismes afin de faire face aux stress et d'éviter des dérèglements comme les maladies inflammatoires et métaboliques.

Identification de modifications post-traductionnelles de Staufen1 et étude de leur fonction régulatrice

La régulation post-transcriptionnelle joue un rôle de premier plan dans le contrôle fin de l’expression génique en permettant une modulation de la synthèse de protéines dans le temps et l’espace, en fonction des besoins de la cellule. Ainsi, des protéines reconnaissant des éléments d’ARN présents sur des transcrits peuvent influencer toutes les étapes de leur existence, soit leur épissage, leur export nucléaire, leur localisation subcellulaire, leur traduction et leur dégradation. Staufen1 (Stau1) est un membre de la famille des protéines liant l’ARN double-brin qui contribue à la régulation post-transcriptionnelle par son implication dans des mécanismes qui vont promouvoir l’épissage alternatif, le transport, la dé-répression de la traduction et l’induction de la dégradation d’ARN messagers (ARNm) spécifiques. L’identité des cibles potentielles de Stau1 est maintenant connue puisqu’une étude à l’échelle du génome a montré que la protéine s’associe à près de 7% du transcriptome des cellules HEK293T. Ces ARNm se classent dans un large éventail de catégories fonctionnelles, mais il est tout de même intéressant de noter qu’une grande proportion d’entre eux code pour des protéines reliées au métabolisme cellulaire et à la régulation de processus cellulaires. En considérant toutes ces informations, nous avons émis l’hypothèse que les différentes activités de Stau1 puissent être modulées afin de contrôler adéquatement l’expression des transcrits liés par la protéine. Dans la mesure où certains ARNm faisant partie des complexes définis par la présence de Stau1 codent pour des régulateurs clés de la prolifération cellulaire, nous avons voulu examiner si l’expression de la protéine varie au cours du cycle de division cellulaire. Nous avons montré que l’abondance de Stau1 est maximale en début de mitose et qu’elle diminue ensuite lorsque les cellules complètent la division cellulaire. Nous avons ensuite découvert que cette baisse d’expression de Stau1 en sortie de mitose dépend du complexe promoteur d’anaphase/cyclosome (APC/C). En soutien à l’idée que Stau1 soit une cible de cette ubiquitine ligase de type E3, nous avons de plus démontré que Stau1 est ubiquitiné et dégradé par le protéasome. Ce contrôle des niveaux de Stau1 semble important puisque la surexpression de la protéine retarde la sortie de mitose et entraîne une diminution importante de la prolifération cellulaire. Par ailleurs, nous avons supposé que les différentes fonctions de Stau1 puissent également être sujettes à une régulation. Compte tenu que les activités de nombreuses protéines liant l’ARN peuvent être contrôlées par des modifications post-traductionnelles telles que la phosphorylation, nous avons voulu tester la possibilité que Stau1 soit phosphorylé. L’immunopurification de Stau1 et son analyse par spectrométrie de masse nous a permis d’identifier trois phosphosites dans la protéine. L’évaluation du rôle de ces événements de phosphorylation à l’aide de mutants phoshomimétiques ou non-phoshorylables a révélé que la modification de Stau1 pourrait compromettre son association à la protéine UPF1. Comme cette interaction est nécessaire pour déstabiliser les transcrits liés par Stau1, nos résultats suggèrent fortement que la fonction de Stau1 dans la dégradation d’ARNm est régulée négativement par sa phosphorylation. Toutes ces données mettent en lumière l’importance des modifications post-traductionnelles telles que l’ubiquitination et la phosphorylation dans la modulation de l’expression et des fonctions de Stau 1. Somme toute, il est vraisemblable que ces mécanismes de contrôle puissent avoir un impact significatif sur le destin des ARNm liés par Stau1, particulièrement dans un contexte de progression dans le cycle cellulaire.

Functional Analysis of Nuclear beta-Adrenergic Receptors in the Myocardium

Récemment plusieurs récepteurs couplés aux protéines G (RCPGs) ont été caractérisés au niveau des membranes intracellulaires, dont la membrane nucléaire. Notre objectif était de déterminer si les sous-types de récepteurs β-adrénergiques (βAR) et leurs machineries de signalisation étaient fonctionnels et localisés à la membrane nucléaire des cardiomyocytes. Nous avons démontré la présence des β1AR et β3AR, mais pas du β2AR à la membrane nucléaire de myocytes ventriculaires adultes par immunobuvardage, par microscopie confocale, et par des essais fonctionnels. De plus, certains partenaires de signalisation comme les protéines GαS, Gαi, l’adénylate cyclase II, et V/VI y étaient également localisés. Les sous-types de βAR nucléaires étaient fonctionnels puisqu'ils pouvaient lier leurs ligands et activer leurs effecteurs. En utilisant des noyaux isolés, nous avons observé que l'agoniste non-sélectif isoprotérénol (ISO), et que le BRL37344, un ligand sélectif du β3AR, stimulaient l'initiation de la synthèse de l’ARN, contrairement à l'agoniste sélectif du β1AR, le xamotérol. Cette synthèse était abolie par la toxine pertussique (PTX). Cependant, la stimulation des récepteurs nucléaires de type B de l’endothéline (ETB) causaient une réduction de l'initiation de la synthèse d’ARN. Les voies de signalisations impliquées dans la régulation de la synthèse d’ARN par les RCPGs ont ensuite été étudiées en utilisant des noyaux isolés stimulés par des agonistes en présence ou absence de différents inhibiteurs des voies MAP Kinases (proteines kinases activées par mitogènes) et de la voie PI3K/PKB. Les protéines impliquées dans les voies de signalisation de p38, JNK, ERK MAP Kinase et PKB étaient présents dans les noyaux isolés. L'inhibition de PKB par la triciribine, inhibait la synthèse d’ARN. Nous avons ensuite pu mettre en évidence par qPCR que la stimulation par l’ISO entrainait une augmentation du niveau d'ARNr 18S ainsi qu’une diminution de l'expression d’ARNm de NFκB. En contraste, l’ET-1 n’avait aucun effet sur le niveau d’expression de l’ARNr 18S. Nous avons ensuite montré que la stimulation par l’ISO réduisait l’expression de plusieurs gènes impliqués dans l'activation de NFκB, tandis que l’inhibition de ERK1/2 et PKB renversait cet effet. Un microarray global nous a ensuite permis de démontrer que les βARs et les ETRs nucléaires régulaient un grand nombre de gènes distincts. Finalement, les βARs et ETRs nucléaires augmentaient aussi une production de NO de noyaux isolés, ce qui pouvait être inhibée par le LNAME. Ces résultats ont été confirmés dans des cardiomyocytes intacts en utilisant des analogues cagés et perméables d’ISO et de l'ET-1: l'augmentation de NO nucléaire détectée par DAF2-DA, causée par l'ET-1 et l'ISO, pouvait être prévenue par le LNAME. Finalement, l’augmentation de l’initiation de la transcription induite par l'ISO était aussi bloquée par le L-NAME ou par un inbitheur de PKG, le KT5823, suggérant que la voie NO-GC-PKG est impliquée dans la régulation de la transcription par les βAR. En conclusion, les βARs et les ETRs nucléaires utilisent des voies de signalisation différentes et exercent ainsi des effets distincts sur l’expression des gènes cardiaques. Ils représentent donc une avenue intéressante pour le développement de drogues pharmacologiques.

Phosphorylation du CTD de l'ARN polymérase II et impact de l'histone H2A.Z sur le positionnement des nucléosomes chez S. cerevisiae

La phosphorylation du domaine C-terminal de l’ARN polymérase II permet à ce complexe protéique d’exécuter la transcription des gènes, en plus de coupler à la transcription des événements moléculaires comme la maturation des ARNm. Mes résultats montrent que même si cette phosphorylation suit un patron similaire à l’ensemble des gènes, il existe des exceptions pouvant être dues à des mécanismes alternatifs de phosphorylation du CTD. Le présent ouvrage s’intéresse également au rôle qu’occupe la variante d’histone H2A.Z dans l’organisation de la chromatine. Des études précédentes on montré que le positionnement de certains nucléosomes le long de l’ADN serait influencé par H2A.Z et aurait une influence sur la capacité de transcrire les gènes. Par une approche génomique utilisant les puces à ADN, j’ai cartographié l’impact de la délétion de H2A.Z sur la structure des nucléosomes. Enfin, des résultats intéressants sur la dynamique d’incorporation de H2A.Z à la chromatine ont été obtenus.

Identification des réseaux transcriptionnnels de résistance aux antifongiques chez Candida albicans

Plusieurs souches cliniques de Candida albicans résistantes aux médicaments antifongiques azolés surexpriment des gènes encodant des effecteurs de la résistance appartenant à deux classes fonctionnelles : i) des transporteurs expulsant les azoles, CDR1, CDR2 et MDR1 et ii) la cible des azoles 14-lanostérol déméthylase encodée par ERG11. La surexpression de ces gènes est due à la sélection de mutations activatrices dans des facteurs de transcription à doigts de zinc de la famille zinc cluster (Zn2Cys6) qui contrôlent leur expression : Tac1p (Transcriptional activator of CDR genes 1) contrôlant l’expression de CDR1 et CDR2, Mrr1p (Multidrug resistance regulator 1), régulant celle de MDR1 et Upc2p (Uptake control 2), contrôlant celle d’ERG11. Un autre effecteur de la résistance clinique aux azoles est PDR16, encodant une transférase de phospholipides, dont la surexpression accompagne souvent celle de CDR1 et CDR2, suggérant que les trois gènes appartiennent au même régulon, potentiellement celui de Tac1p. De plus, la régulation transcriptionnelle du gène MDR1 ne dépend pas seulement de Mrr1p, mais aussi du facteur de transcription de la famille basic-leucine zipper Cap1p (Candida activator protein 1), un régulateur majeur de la réponse au stress oxydatif chez C. albicans qui, lorsque muté, induit une surexpression constitutive de MDR1 conférant la résistance aux azoles. Ces observations suggèrent qu’un réseau de régulation transcriptionnelle complexe contrôle le processus de résistance aux antifongiques azolés chez C. albicans. L’objectif de mon projet au doctorat était d’identifier les cibles transcriptionnelles directes des facteurs de transcription Tac1p, Upc2p et Cap1p, en me servant d’approches génétiques et de génomique fonctionnelle, afin de i) caractériser leur réseau transcriptionnel et les modules transcriptionnels qui sont sous leur contrôle direct, et ii) d’inférer leurs fonctions biologiques et ainsi mieux comprendre leur rôle dans la résistance aux azoles. Dans un premier volet, j’ai démontré, par des expériences de génétique, que Tac1p contrôle non seulement la surexpression de CDR1 et CDR2 mais aussi celle de PDR16. Mes résultats ont identifié une nouvelle mutation activatrice de Tac1p (N972D) et ont révélé la participation d’un autre régulateur dans le contrôle transcriptionnel de CDR1 et PDR16 dont l’identité est encore inconnue. Une combinaison d’expériences de transcriptomique et d’immunoprécipitation de la chromatine couplée à l’hybridation sur des biopuces à ADN (ChIP-chip) m’a permis d’identifier plusieurs gènes dont l’expression est contrôlée in vivo et directement par Tac1p (PDR16, CDR1, CDR2, ERG2, autres), Upc2p (ERG11, ERG2, MDR1, CDR1, autres) et Cap1p (MDR1, GCY1, GLR1, autres). Ces expériences ont révélé qu’Upc2p ne contrôle pas seulement l’expression d’ERG11, mais aussi celle de MDR1 et CDR1. Plusieurs nouvelles propriétés fonctionnelles de ces régulateurs ont été caractérisées, notamment la liaison in vivo de Tac1p aux promoteurs de ses cibles de façon constitutive et indépendamment de son état d’activation, et la liaison de Cap1p non seulement à la région du promoteur de ses cibles, mais aussi celle couvrant le cadre de lecture ouvert et le terminateur transcriptionnel putatif, suggérant une interaction physique avec la machinerie de la transcription. La caractérisation du réseau transcriptionnel a révélé une interaction fonctionnnelle entre ces différents facteurs, notamment Cap1p et Mrr1p, et a permis d’inférer des fonctions biologiques potentielles pour Tac1p (trafic et la mobilisation des lipides, réponse au stress oxydatif et osmotique) et confirmer ou proposer d’autres fonctions pour Upc2p (métabolisme des stérols) et Cap1p (réponse au stress oxydatif, métabolisme des sources d’azote, transport des phospholipides). Mes études suggèrent que la résistance aux antifongiques azolés chez C. albicans est intimement liée au métabolisme des lipides membranaires et à la réponse au stress oxydatif.

Caractérisation clinique et génétique d’une famille canadienne-française atteinte de la neuropathie héréditaire sensitive avec rétinite pigmentaire et ataxie

La complexité de l’étude des neuropathies héréditaires provient de leur hétérogénéité clinique et génétique et de la diversité des fibres composant les nerfs périphériques. Cette complexité se reflète dans les nombreuses classifications différentes. Les neuropathies héréditaires se classifient entre autres selon leur mode de transmission et leur atteinte sensitive, autonomique et motrice. Les neuropathies héréditaires sensitives et autonomiques (NHSA) se présentent avec une perte de la sensation distale aux membres, accompagnée d’autres manifestations selon le type de NHSA. L’étude des NHSA est facilitée lorsqu’il existe des grappes de familles originaires de régions du Québec où des effets fondateurs pour des maladies récessives ont déjà été identifiés. Nous avons recruté une grande famille canadienne-française originaire de Paspébiac dans la Baie-des-Chaleurs dans laquelle nous avons identifié quatre cas atteints d’une neuropathie héréditaire sensitive avec rétinite pigmentaire et ataxie (NHSRPA). Nous avons émis l’hypothèse que nous étions en présence d’une nouvelle forme de neuropathie héréditaire sensitive récessive à effet fondateur. Afin d’identifier la position chromosomique du gène muté responsable de la NHSRPA, nous avons tout d’abord complété un criblage du génome en génotypant des marqueurs microsatellites «single tandem repeat» (STR) sur des individus clés et nous avons ensuite procédé à une analyse de liaison génétique paramétrique. Ces études nous ont permis de lier cette famille au chromosome 1 et de définir un premier intervalle candidat de 6,7 Mb. Grâce à un génotypage de marqueurs «single nucleotide polymorphism» (SNP), nous avons réduit l’intervalle candidat à 5,3 Mb au locus 1q32,2-q32,3. Cette région contient 44 gènes candidats. Une revue plus fine de la littérature a fait ressortir qu’une famille espagnole et une américaine de souche hollandaise souffrant de la même maladie avaient déjà été liées au même locus. L’origine possiblement basque de notre famille gaspésienne nous a poussé à comparer l’haplotype porteur avec celui de la famille espagnole qui, quoi que gitane, provient du pays basque espagnol. Ces travaux ont démontré le partage d’une région de 203 kb. Afin de rétrécir davantage notre intervalle candidat, nous avons comparé les haplotypes des cas entre les deux familles et nous avons identifié un dernier intervalle candidat de 60 SNP au locus 1q32,3. Cette région ne contient que quatre gènes candidats dont le plus intéressant est le gène «activating transcription factor» (ATF3). À ce jour, aucune mutation n’a été trouvée dans le gène ATF3 et les gènes FAM71A, BATF3 et NSL1. Des expériences supplémentaires sont nécessaires afin d’identifier le gène muté responsable de la NHSRPA.

Étude de facteurs génétiques prédictifs dans le neuroblastome, en particulier les anomalies du chromosmoe 14q

Le neuroblastome (NB) représente 8% de tous les cancers pédiatriques et est caractérisé par sa grande hétérogénéité clinique. Afin d’évaluer son pronostic, plusieurs facteurs génétiques sont utilisés : amplification de MYCN, délétion 1p, gain 11q et gain 17q. Les buts de notre travail étaient d’abord de vérifier si l’hybridation in situ en fluorescence (FISH) permet une analyse complète de ces anomalies et ensuite, en utilisant une analyse globale du génome telle le polymorphisme nucléotidique simple (SNP), de vérifier la concordance avec les résultats de la FISH et le pronostic potentiel des anomalies du 14q, en particulier du gène AKT. Nous avons donc établi un panel de sondes pour la FISH qui a été appliqué sur 16 tumeurs non-fixées. Après isolation de l’ADN de 36 tumeurs, nous avons effectué une analyse génotypique par SNP utilisant les puces « Affymetrix Genome-Wide Human SNP Array 6.0 » contenant 945,826 sondes non polymorphiques et 906,000 sondes polymorphiques. Nos résultats ont démontré que la FISH permet l’évaluation complète des anomalies génétiques importantes du NB et que les anomalies déséquilibrées sont détectées très précisément par SNP. Les anomalies du 14q tendent à être associées avec des facteurs cliniques comme le grade et l’évolution, contrairement aux anomalies d’AKT. L’analyse du 14q a révélé trois gènes d’intérêt, MAX, BCL11B et GPHN, qui devraient être analysés sur un plus grand échantillon. Ainsi, l’étude par FISH semble adaptée pour détecter les anomalies génétiques classiques du NB, alors que celles retrouvées en 14q représentent de potentielles cibles thérapeutiques pour cette tumeur.

Étude génétique et fonctionnelle de variantes de la région chromosomique 3p21 associée aux maladies inflammatoires de l'intestin

Des études de liaison et d’association génétiques ont permis d’identifier certains des facteurs de risque génétiques aux maladies inflammatoires de l’intestin (MII) dans la région chromosomique 3p21. Dans cette région, le polymorphisme nucléotidique simple (SNP) codant non-synonyme du gène MST1, rs3197999, encodant pour la mutation R689C, a été associé et répliqué à la fois à la colite ulcéreuse (CU) et à la maladie de Crohn (MC). Un autre SNP, corrélé à des SNP codants non-synonymes du gène MST1R, a également été associé à la MC. Afin de déterminer si d’autres variantes des gènes MST1 et MST1R sont associés à la CU, nous avons testé pour association des SNP de ces gènes. Seul un proxy de R689C a montré un signal d’association significatif aux MII, ce qui suggère que R689C est la variante causale aux MII dans le gène MST1. En cherchant à déterminer si la région 3p21 contenait plusieurs signaux d’association mutuellement indépendants, trois SNP ont été identifiés comme possible facteurs de risque indépendants, et ont été génotypés dans des cas de CU et de MC et des témoins, puis nos résultats d’association ont été combinés à ceux provenant de trois autres cohortes indépendantes. Les trois SNP, R689C (MST1), rs6802890 et rs7629936 (CDHR4), sont associés aux MII, mais une étude d’association conditionnelle suggère qu’il existe en fait deux signaux d’association mutuellement indépendants dans la région 3p21. Le signal principal provient de R689C, une mutation de la protéine MSP. Cette protéine a un rôle dans l’inflammation chez les macrophages murins, et la migration, la cicatrisation et la survie chez les cellules épithéliales. Dans cette étude, le rôle de la MSP a été investigué dans des modèles de macrophages humains et de cellules épithéliales de côlon, et seule la phosphorylation d’AKT, un acteur dans la voie de signalisation de la survie cellulaire, a été modulée par la MSP dans nos modèles. Ce projet a donc permis d’apporter des connaissances sur les facteurs de risques génétiques aux MII dans la région 3p21, en identifiant 2 signaux d’association indépendants, et en nous informant sur le rôle de MST1, duquel provient le signal d’association principal, chez les cellules humaines.

Développement d’outils pour l’analyse de données de ChIP-seq et l’identification des facteurs de transcription

La méthode ChIP-seq est une technologie combinant la technique de chromatine immunoprecipitation avec le séquençage haut-débit et permettant l’analyse in vivo des facteurs de transcription à grande échelle. Le traitement des grandes quantités de données ainsi générées nécessite des moyens informatiques performants et de nombreux outils ont vu le jour récemment. Reste cependant que cette multiplication des logiciels réalisant chacun une étape de l’analyse engendre des problèmes de compatibilité et complique les analyses. Il existe ainsi un besoin important pour une suite de logiciels performante et flexible permettant l’identification des motifs. Nous proposons ici un ensemble complet d’analyse de données ChIP-seq disponible librement dans R et composé de trois modules PICS, rGADEM et MotIV. A travers l’analyse de quatre jeux de données des facteurs de transcription CTCF, STAT1, FOXA1 et ER nous avons démontré l’efficacité de notre ensemble d’analyse et mis en avant les fonctionnalités novatrices de celui-ci, notamment concernant le traitement des résultats par MotIV conduisant à la découverte de motifs non détectés par les autres algorithmes.