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.

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

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.

Suppression of long-branch attraction artefacts in the animal phylogeny using a site-heterogeneous model

Affiliation: Département de Biochimie, Université de Montréal

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

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

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 sur le rôle des déséquilibres génomiques dans le Syndrome d’Impatiences Musculaires de l’Éveil

Le Syndrome d’Impatiences Musculaires de l’Éveil (SIME) est une maladie neurologique caractérisée par un besoin urgent de bouger les jambes. C’est également l’une des causes les plus fréquentes d’insomnie. C’est une maladie très répandue, avec une prévalence de presque 15 % dans la population générale. Les maladies multifactorielles comme le SIME sont souvent le résultat de l’évolution d’une composante génétique et d’une composante environnementale. Dans le cadre du SIME, les études d’association génomique ont permis l’identification de 4 variants à effet modéré ou faible. Cependant, ces quatre variants n’expliquent qu’une faible partie de la composante génétique de la maladie, ce qui confirme que plusieurs nouveaux variants sont encore à identifier. Le rôle des déséquilibres génomiques (Copy Number Variations ou CNVs) dans le mécanisme génétique du SIME est à ce jour inconnu. Cependant, les CNVs se sont récemment positionnés comme une source d’intérêt majeur de variation génétique potentiellement responsable des phénotypes. En collaboration avec une équipe de Munich, nous avons réalisé deux études CNVs à échelle génomique (biopuces à SNP et hybridation génomique comparée (CGH)) sur des patients SIME d’ascendance germanique. À l’aide d’une étude cas-contrôle, nous avons pu identifier des régions avec une occurrence de CNVs différentes pour les patients SIME, comparés à différents groupes contrôles. L’une de ces régions est particulièrement intéressante, car elle est concordante à la fois avec des précédentes études familiales ainsi qu’avec les récentes études d’associations génomiques.

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

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.

Effect of fatty acids on hyphal growth in the pathogenic yeast Candida albicans

Candida albicans est une levure pathogène qui, à l’état commensal, colonise les muqueuses de la cavité orale et du tractus gastro-intestinal. De nature opportuniste, C. albicans cause de nombreuses infections, allant des candidoses superficielles (muguet buccal, vulvo-vaginite) aux candidoses systémiques sévères. C. albicans a la capacité de se développer sous diverses morphologies, telles que les formes levures, pseudohyphes et hyphes. Des stimuli environnementaux mimant les conditions retrouvées chez l’hôte (température de 37°C, pH neutre, présence de sérum) induisent la transition levure-à-hyphe (i.e. morphogenèse ou filamentation). Cette transition morphologique contribue à la pathogénicité de C. albicans, du fait que des souches présentant un défaut de filamentation sont avirulentes. Non seulement la morphogenèse est un facteur de virulence, mais elle constituerait aussi une cible pour le développement d’antifongiques. En effet, il a déjà été démontré que l’inhibition de la transition levure-à-hyphe atténuait la virulence de C. albicans lors d’infections systémiques. Par ailleurs, des études ont démontré que de nombreuses molécules pouvaient moduler la morphogenèse. Parmi ces molécules, certains acides gras, dont l’acide linoléique conjugué (CLA), inhibent la formation d’hyphes. Ainsi, le CLA posséderait des propriétés thérapeutiques, du fait qu’il interfère avec un déterminant de pathogénicité de C. albicans. Par contre, avant d’évaluer son potentiel thérapeutique dans un contexte clinique, il est essentiel d’étudier son mode d’action. Ce projet vise à caractériser l’activité anti-filamentation des acides gras et du CLA et à déterminer le mécanisme par lequel ces molécules inhibent la morphogenèse chez C. albicans. Des analyses transcriptomiques globales ont été effectuées afin d’obtenir le profil transcriptionnel de la réponse de C. albicans au CLA. L’acide gras a entraîné une baisse des niveaux d’expression de gènes encodant des protéines hyphes-spécifiques et des régulateurs de morphogenèse, dont RAS1. Ce gène code pour la GTPase Ras1p, une protéine membranaire de signalisation qui joue un rôle important dans la transition levure-à-hyphe. Des analyses de PCR quantitatif ont confirmé que le CLA inhibait l’induction de RAS1. De plus, le CLA a non seulement causé une baisse des niveaux cellulaires de Ras1p, mais a aussi entraîné sa délocalisation de la membrane plasmique. En affectant les niveaux et la localisation cellulaire de Ras1p, le CLA nuit à l’activation de la voie de signalisation Ras1p-dépendante, inhibant ainsi la morphogenèse. Il est possible que le CLA altère la structure de la membrane plasmique et affecte indirectement la localisation membranaire de Ras1p. Ces travaux ont permis de mettre en évidence le mode d’action du CLA. Le potentiel thérapeutique du CLA pourrait maintenant être évalué dans un contexte d’infection, permettant ainsi de vérifier qu’une telle approche constitue véritablement une stratégie pour le traitement des candidoses.

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