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 clinique et génétique d’une nouvelle forme d’ataxie spinocérébelleuse pure associée à l’Érythrokératodermie

Nous présentons ici la description clinique et génétique d’un syndrome neurocutané unique. Le laboratoire du Dr Cossette a entrepris la caractérisation clinique et génétique d'une famille canadienne-française qui a été identifiée par les Drs Giroux et Barbeau en 1972 et qui comprend plus de 100 personnes sur six générations. Les membres atteints de cette famille présentent des lésions typiques d'érythrokératodermie (EK) (OMIM 133190, EKV1 et EKV2), associées à une ataxie spinocérébelleuse pure. Dans cette famille, l'ataxie est caractérisée par des troubles de la coordination et de la démarche causés par une dégénérescence du cervelet et de la moelle épinière. Cette ataxie est transmise selon un mode autosomique dominant. Une étude antérieure de cette variante d'EK avec ataxie avait suggéré une liaison sur le chromosome 1p34-p35, soit la même région que les formes EKV de type 1 et 2, causées respectivement par des mutations dans les gènes connexin-31 (GJB3; OMIM 603324) et connexin-30.3 (GJB4; OMIM 605425). Cependant, aucune mutation n'a été retrouvée dans ces gènes pour la famille canadienne-française. Nous avons récemment recontacté la famille et effectué des examens détaillés, incluant une imagerie par résonance magnétique (IRM) et un électromyogramme (EMG). Les manifestations neurologiques des individus atteints sont compatibles avec une nouvelle forme d’ataxie cérébelleuse pure à transmission autosomique dominante (ADCA de type III dans la classification de Harding) que nous avons appelée SCA34. Une cartographie complète du génome nous a permis de localiser le gène SCA34 sur le chromosome 6p12.3-q16.2. Également, en collaboration avec les Drs Alexis Brice (Hôpital Pitié-La Salpêtrière, Paris) et Alfredo Brusco (Hôpital San Giovanni Battista di Torino, Italie), nous avons confirmé que trois autres familles européennes avec SCA inexpliquée étaient également liées au locus SCA34. Notre laboratoire a récemment entrepris la recherche des mutations responsables de SCA34. Les résultats de ce criblage de gènes candidats sont présentés dans le chapitre 3 de cette thèse.

Exclusion de liaison génétique au locus SPAX2 de cas canadiens-français d’ataxie spastique

Les ataxies héréditaires sont des désordres neuro-dégénératifs qui causent une ataxie comme symptôme primaire; soit une perte de coordination des mouvements volontaires, un sens de l’équilibre déficient et un trouble à la motricité. Elles forment un groupe cliniquement et génétiquement hétérogène. De ce fait, de nombreuses classifications existent basées sur différents critères. Cependant, le consensus actuel veut que le mode de transmission soit le critère premier de classement. On estime la prévalence mondiale des ataxies héréditaires à 6/100 000 bien que ce nombre diffère entre régions. C’est le cas du Québec où la structuration historique du bassin génétique canadien-français a menée à des effets fondateurs régionaux, ce qui a eu comme conséquence de hausser la prévalence régionale de certaines maladies. L’Acadie est également une région canadienne-française avec des effets fondateurs où le taux de prévalence de certaines ataxies héréditaires est plus élevé. Nous avons recruté huit familles canadiennes-françaises provenant de diverses régions du Québec, ayant un lien génétique plus ou moins rapproché avec l’Acadie, dans lesquelles nous avons observé dix cas d’une forme d’ataxie spastique autosomique récessive relativement légère qui a résistée à l’analyse des gènes d’ataxies connues. Nous avons émis l’hypothèse d’être en présence d’une nouvelle forme d’ataxie à effet fondateur pour la population canadienne-française. Afin d’identifier le gène muté responsable de cette ataxie, un criblage génomique des marqueurs SNP pour les individus recrutés fut effectué. Puis, par cartographie de l’homozygotie, une région de 2,5 Mb fut identifiée sur le chromosome 17p13 dans une famille. Une revue de la littérature nous a permis de constater, qu’en 2007, quatre familles nord-africaines atteintes d’une ataxie dénommée SPAX2 qui présentaient des manifestations cliniques semblables avaient déjà été liées au même locus sur le chromosome 17. Afin de supporter notre hypothèse que les malades étaient porteurs de deux copies de la même mutation fondatrice et de cartographier plus finement notre région d’intérêt, les haplotypes de tous les atteints de nos huit familles furent étudiés. Nous avons établie qu’un intervalle de 200 kb (70 SNP), soit du marqueur rs9900036 à rs7222052, était partagé par tous nos participants. Les deux gènes les plus prometteurs des 18 se trouvant dans la région furent séquencés. Aucune mutation ne fut trouvée dans les gènes SLC25A11 et KIF1C. Par la suite, une analyse de liaison génétique stricte avec calcul de LOD score nous a permis d’exclure ce locus de 200 kb comme étant celui porteur du gène muté causant l’ataxie dans la majorité de nos familles. Nous avons donc conclus que malgré qu’une famille soit homozygote pour une grande région du chromosome 17, l’absence d’Informativité des marqueurs SNP dans la région de 200 kb fut responsable de l’apparent partage d’haplotype homozygote. Le travail reste donc entier afin d’identifier les mutations géniques responsables de la présentation ataxique chez nos participants de souche acadienne.

Identification des gènes responsables des hyperplasies surrénaliennes macronodulaires bilatérales familiales avec récepteurs aberrants

La majorité des hyperplasies macronodulaires bilatérales des surrénales avec syndrome de Cushing ACTH-indépendant (AIMAH) est due à l’expression aberrante de divers récepteurs hormonaux au niveau du cortex surrénalien. Les gènes responsables des AIMAH familiales avec récepteurs aberrants n’ont pas été identifiés. Le but de ce projet est de les identifier. Une étude de liaison, visant à identifier la ou les régions du génome comprenant le ou les gènes pouvant être en cause dans les AIMAH familiales, a été réalisée en utilisant l’ADN des membres d’une famille (10 malades et 7 sains) originaire du Québec, atteinte d’AIMAH et syndrome de Cushing et caractérisée par l’expression des récepteurs β-adrénergique et V1-vasopressine. Diverses régions chromosomiques entre les personnes atteintes et non-atteintes de la famille ont été soulignées. Un total de 707453 SNPs a été obtenu, et après analyse statistique, 159 SNPs significatifs, pouvant être associés au phénotype, ont été mis en évidence entre les deux groupes. Il a été constaté que la majorité de ces SNPs se situaient sur les régions chromosomiques 1q32.1 et 16q12.2. Une étude du transcriptome a aussi été réalisée en utilisant l’ADN des tumeurs de deux patients de la famille, ainsi que l’ADN d'autres tumeurs surrénaliennes. Les analyses statistiques ont permis d’identifier 15 gènes susceptibles d’être reliés à la maladie (11 surexprimés et 4 sous-exprimés). En utilisant les données de ces deux études, nous avons ciblé six gènes du chromosome 1 (ATP2B4, PPP1R12B, SOX13, CACNA1S, ADORA1et PHLDA3), un du chromosome 16 (CHD9) et un du chromosome 13 (SPRY2), afin de rechercher la présence de mutations. Le séquençage n’a révélé aucun changement de nucléotide dans les gènes PPP1R12B et SOX13. Dans les gènes ATP2B4, CACNA1S, ADORA1et PHLDA3, le séquençage a révélé des changements de nucléotides n’entrainant soit pas de changement d’acide aminé soit un changement d’acide aminé jugé « non pertinent », du fait qu’il ne permettait pas de différencier les sujets sains des sujets atteints. Pour ce qui est de CHD9 et SPRY2, le séquençage a permis d’identifier des changements de nucléotides entrainant des changements d’acides aminés de façon plus fréquente chez les sujets atteints par rapport aux sujets sains. En conclusion, nos travaux nous ont donc permis d’identifier, par étude de liaison et par analyse du transcriptome, des gènes candidats qui pourraient être responsables de cette pathologie. Le séquençage de ces gènes candidats a révélé des mutations de CHD9 et SPRY2. Ces résultats s’avèrent prometteurs puisque ces deux gènes produisent des protéines impliquées dans le remodelage de la chromatine et dans la régulation de la signalisation des protéines kinases. Le phénotypage et le génotypage des patients atteints doivent être poursuivis pour vérification.

Identification des bases génétiques des myopathies à multi-minicores avec ou sans cardiomyopathie

Thèse réalisée en cotutelle avec l'Université Pierre et Marie Curie, Paris 6(UPMC, Paris, France).

Cartographie génétique d’une forme familiale autosomale dominante d’anévrysmes du septum inter-ventriculaire et de communications inter-ventriculaires au chromosome 10p15

Les malformations cardiaques congénitales (CHM) représentent 28 % de toutes les malformations congénitales majeures et touchent 8 pour 1000 naissances à terme. Elles sont la cause de mortalité et de morbidité non infectieuse la plus fréquente chez les enfants de moins d’une année de vie. Les communications interventriculaires (VSD) forment le sous-type de CHM le plus fréquent et l’aggrégation familiale est extrêmement rare. Le but de cette étude était d’identifier les facteurs génétiques et les régions chromosomiques contribuant aux VSD. Une grande famille ségréguant diverses formes de pathologies septales, incluant des VSD, des anévrysmes du septum interventriculaire (VSA) et des communications interauriculaires (ASD), a été examinées et caractérisées cliniquement et génétiquement. Dix-huit membres de la famille, sur trois générations, ont pu être étudiés. (10 affectés : 4 VSD, 3 VSA, 2 ASD et une tétralogie de Fallot). L’analyse de liaison multipoint paramétrique démontre un logarithme des probabilités maximal (LOD) de 3.29 liant significativement le chromosome 10p15.3-10p15.2 aux traits observés dans cette famille. Le pointage LOD oriente vers une région pauvre en gènes qui a déjà été associée aux malformations du septum interventriculaire, mais qui est distincte de la région du syndrome de DiGeorge de type 2 sur le chromosome 10p. De plus, plusieurs scénarios d’analyse de liaison suggèrent que la tétralogie de Fallot est une phénocopie et qu’elle est donc génétiquement différente des autres pathologies du septum observées dans cette famille. En bref, cette étude associe une forme rare de VSD/VSA au chromosome 10p15 et permet d’étendre le spectre de l’hétérogénéité des pathologies septales. Mots-clés : Malformations cardiaques congénitales, malformations du septum, tétralogie de Fallot, analyse de liaison, chromosome 10p15, génétique moléculaire

Caractérisation génétique d’une forme d’ataxie tardive.

Les ataxies forment un groupe de maladies neurodégénératives qui sont caractérisées par un manque de coordination des mouvements volontaires. Mes travaux ont porté sur une forme d'ataxie à début tardif (LOCA), après l’âge de 50 ans. Les principales caractéristiques cliniques sont: atrophie cérébelleuse à l’IRM (88%), dysarthrie (81%), atrophie du lobe frontal (50%) et nystagmus (52%). La ségrégation dans les familles de cette ataxie est en faveur d’une transmission récessive. Afin d'identifier le gène responsable de LOCA, nous avons recruté 38 patients affectés d'une forme tardive d'ataxie, issus du SLSJ, des Cantons de l’Est ou d’autres régions du Québec. Un premier criblage du génome a été effectué avec des marqueurs microsatellites sur une famille clé. Une analyse de liaison paramétrique nous a suggéré une liaison au chromosome 13 (4.4Mb). Une recherche d’un haplotype partagé entre 17 familles LOCA a diminué la taille de l'intervalle candidat à 1.6Mb, mais l’haplotype s’est avéré fréquent dans la population canadienne-française. Un second criblage du génome avec des marqueurs SNP nous a permis d’évaluer par cartographie d’homozygotie la possibilité qu’une mutation fondatrice partagée dans des sous-groupes de malades. Plusieurs stratégies d'analyse ont été effectuées, entre autre par regroupement régional. Aucun loci candidats ne fut identifié avec confiance. Nous avons donc combiné les données de génotypage avec le séquençage exomique afin d'identifier le gène responsable. L'analyse de six individus atteints nous a permis d'obtenir une liste de variants rare contenant quatre gènes potentiels. Cette analyse doit se poursuivre pour identifier le gène responsable de LOCA.

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.

Exploration génétique de l'hypothyroïdie congénitale par dysgénésie thyroïdienne

L'hypothyroïdie congénitale par dysgénésie thyroïdienne (HCDT, ectopie dans plus de 80 %) a une prévalence de 1 cas sur 4000 naissances vivantes. L’HCDT est la conséquence d'une défaillance de la thyroïde embryonnaire à se différencier, à se maintenir ou à migrer vers sa localisation anatomique (partie antérieure du cou), qui aboutit à une absence totale de la thyroïde (athyréose) ou à une ectopie thyroïdienne (linguale ou sublinguale). Les HCDT sont principalement non-syndromiques (soit 98% des cas sont non-familiale), ont un taux de discordance de 92% chez les jumeaux monozygotes, et ont une prédominance féminine et ethnique (i.e., Caucasienne). La majorité des cas d’HCDT n’a pas de cause connue, mais est associée à un déficit sévère en hormones thyroïdiennes (hypothyroïdie). Des mutations germinales dans les facteurs de transcription liés à la thyroïde (NKX2.1, FOXE1, PAX8, NKX2.5) ont été identifiées dans seulement 3% des patients atteints d’HCDT sporadiques et l’analyse de liaisons exclue ces gènes dans les rares familles multiplex avec HCDT. Nous supposons que le manque de transmission familiale claire d’HCDT peut résulter de la nécessité d’au moins deux « hits » génétiques différents dans des gènes importants pour le développement thyroïdien. Pour répondre au mieux nos questions de recherche, nous avons utilisé deux approches différentes: 1) une approche gène candidat, FOXE1, seul gène impliqué dans l’ectopie dans le modèle murin et 2) une approche en utilisant les techniques de séquençage de nouvelle génération (NGS) afin de trouver des variants génétiques pouvant expliquer cette pathologie au sein d’une cohorte de patients avec HCDT. Pour la première approche, une étude cas-contrôles a été réalisée sur le promoteur de FOXE1. Il a récemment été découvert qu’une région du promoteur de FOXE1 est différentiellement méthylée au niveau de deux dinucléotides CpG consécutifs, définissant une zone cruciale de contrôle de l’expression de FOXE1. L’analyse d’association basée sur les haplotypes a révélé qu’un haplotype (Hap1: ACCCCCCdel1C) est associé avec le HCDT chez les Caucasiens (p = 5x10-03). Une réduction significative de l’activité luciférase est observée pour Hap1 (réduction de 68%, p<0.001) comparé au promoteur WT de FOXE1. Une réduction de 50% de l’expression de FOXE1 dans une lignée de cellules thyroïdienne humaine est suffisante pour réduire significativement la migration cellulaire (réduction de 55%, p<0.05). Un autre haplotype (Hap2: ACCCCCCC) est observé moins fréquemment chez les Afro-Américain comparés aux Caucasiens (p = 1.7x10-03) et Hap2 diminue l’activité luciférase (réduction de 26%, p<0.001). Deux haplotypes distincts sont trouvés fréquemment dans les contrôles Africains (Black-African descents). Le premier haplotype (Hap3: GTCCCAAC) est fréquent (30.2%) chez les contrôles Afro-Américains comparés aux contrôles Caucasiens (6.3%; p = 2.59 x 10-9) tandis que le second haplotype (Hap4: GTCCGCAC) est trouvé exclusivement chez les contrôles Afro-Américains (9.4%) et est absent chez les contrôles Caucasiens (P = 2.59 x 10-6). Pour la deuxième approche, le séquençage de l’exome de l’ADN leucocytaire entre les jumeaux MZ discordants n’a révélé aucune différence. D'où l'intérêt du projet de séquençage de l’ADN et l’ARN de thyroïdes ectopiques et orthotopiques dans lesquelles de l'expression monoallélique aléatoire dans a été observée, ce qui pourrait expliquer comment une mutation monoallélique peut avoir des conséquences pathogéniques. Finalement, le séquençage de l’exome d’une cohorte de 36 cas atteints d’HCDT a permis d’identifier de nouveaux variants probablement pathogéniques dans les gènes récurrents RYR3, SSPO, IKBKE et TNXB. Ces quatre gènes sont impliqués dans l’adhésion focale (jouant un rôle dans la migration cellulaire), suggérant un rôle direct dans les défauts de migration de la thyroïde. Les essais de migration montrent une forte diminution (au moins 60% à 5h) de la migration des cellules thyroïdiennes infectées par shRNA comparés au shCtrl dans 2 de ces gènes. Des zebrafish KO (-/- et +/-) pour ces nouveaux gènes seront réalisés afin d’évaluer leur impact sur l’embryologie de la thyroïde.