Genomic architecture of sickle cell disease clinical variation in children from West Africa : a case-control study design

Contexte : L’anémie falciforme ou drépanocytose est un problème de santé important, particulièrement pour les patients d’origine africaine. La variation phénotypique de l’anémie falciforme est problématique pour le suivi et le traitement des patients. L’architecture génomique responsable de cette variabilité est peu connue. Principe : Mieux saisir la contribution génétique de la variation clinique de cette maladie facilitera l’identification des patients à risque de développer des phénotypes sévères, ainsi que l’adaptation des soins. Objectifs : L’objectif général de cette thèse est de combler les lacunes relatives aux connaissances sur l’épidémiologie génomique de l’anémie falciforme à l’aide d’une cohorte issue au Bénin. Les objectifs spécifiques sont les suivants : 1) caractériser les profils d’expressions génomiques associés à la sévérité de l’anémie falciforme ; 2) identifier des biomarqueurs de la sévérité de l’anémie falciforme ; 3) identifier la régulation génétique des variations transcriptionelles ; 4) identifier des interactions statistiques entre le génotype et le niveau de sévérité associé à l’expression ; 5) identifier des cibles de médicaments pour améliorer l’état des patients atteints d’anémie falciforme. Méthode : Une étude cas-témoins de 250 patients et 61 frères et soeurs non-atteints a été menée au Centre de Prise en charge Médical Intégré du Nourrisson et de la Femme Enceinte atteints de Drépanocytose, au Bénin entre février et décembre 2010. Résultats : Notre analyse a montré que des profils d’expressions sont associés avec la sévérité de l’anémie falciforme. Ces profils sont enrichis de génes des voies biologiques qui contribuent à la progression de la maladie : l’activation plaquettaire, les lymphocytes B, le stress, l’inflammation et la prolifération cellulaire. Des biomarqueurs transcriptionnels ont permis de distinguer les patients ayant des niveaux de sévérité clinique différents. La régulation génétique de la variation de l’expression des gènes a été démontrée et des interactions ont été identifiées. Sur la base de ces résultats génétiques, des cibles de médicaments sont proposées. Conclusion: Ce travail de thèse permet de mieux comprendre l’impact de la génomique sur la sévérité de l’anémie falciforme et ouvre des perspectives de développement de traitements ciblés pour améliorer les soins offerts aux patients.

Genetic determinants of clinical heterogeneity in sickle cell disease

L’anémie falciforme est une maladie monogénique causée par une mutation dans le locus de la β-globine. Malgré le fait que l’anémie falciforme soit une maladie monogénique, cette maladie présente une grande hétérogénéité clinique. On présume que des facteurs environnementaux et génétiques contribuent à cette hétérogénéité. Il a été observé qu’un haut taux d’hémoglobine fœtale (HbF) diminuait la sévérité et la mortalité des patients atteints de l’anémie falciforme. Le but de mon projet était d’identifier des variations génétiques modifiant la sévérité clinique de l’anémie falciforme. Dans un premier temps, nous avons effectué la cartographie-fine de trois régions précédemment associées avec le taux d’hémoglobine fœtale. Nous avons ensuite effectué des études d’association pan-génomiques avec deux complications cliniques de l’anémie falciforme ainsi qu’avec le taux d’hémoglobine fœtale. Hormis les régions déjà identifiées comme étant associées au taux d’hémoglobine fœtale, aucun locus n’a atteint le niveau significatif de la puce de génotypage. Pour identifier des groupes de gènes modérément associés au taux d’hémoglobine fœtale qui seraient impliqués dans de mêmes voies biologiques, nous avons effectué une étude des processus biologiques. Finalement, nous avons effectué l’analyse de 19 exomes de patients Jamaïcains ayant des complications cliniques mineures de l’anémie falciforme. Compte tenu de la taille des cohortes de réplication disponibles, nous n’avons pas les moyens de valider statistiquement les variations identifiées par notre étude. Cependant, nos résultats fournissent de bons gènes candidats pour des études fonctionnelles et pour les réplications futures. Nos résultats suggèrent aussi que le β-hydroxybutyrate en concentration endogène pourraient influencer le taux d’hémoglobine fœtale. De plus, nous montrons que la cartographie-fine des régions associées par des études pan-génomiques peut identifier des signaux d’association additionnels et augmenter la variation héritable expliquée par cette région.

Syndrome du QT long: étude clinique à l’Institut Cardiologique de Montréal et recherche de nouvelles variantes causales par séquençage à haut débit

Le syndrome du QT long congénital (LQTS) est une canalopathie génétique, à l’origine de syncopes et mort subite. Le dépistage génétique identifie des variantes génétiques dans 50-70% des cas, suggérant l’implication d’autres gènes. Nous avons recueilli les caractéristiques des patients avec un LQTS à l’ICM, et recruté 12 patients avec un génotype négatif pour le LQTS pour un séquençage à haut débit des exons afin d’identifier de nouvelles variantes causales. Nous avons développé une approche analytique par étapes : (1) les gènes connus du, (2) les gènes dans des loci identifiés par des études d’association sur le QT, et (3) les gènes montrant la même variante chez plusieurs patients. L’analyse génétique a identifié de nouvelles variantes dans: (1) KCNJ2, ANK2 et AKAP9, et (2) dans NOS1AP. (3) Deux patientes avec des phénotypes semblables présentent la même variante homozygote dans TECRL, un nouveau gène candidat dont le rôle est inconnu.

Searching for novel gene functions in yeast : identification of thousands of novel molecular interactions by protein-fragment complementation assay followed by automated gene function prediction and high-throughput lipidomics

La compréhension de processus biologiques complexes requiert des approches expérimentales et informatiques sophistiquées. Les récents progrès dans le domaine des stratégies génomiques fonctionnelles mettent dorénavant à notre disposition de puissants outils de collecte de données sur l’interconnectivité des gènes, des protéines et des petites molécules, dans le but d’étudier les principes organisationnels de leurs réseaux cellulaires. L’intégration de ces connaissances au sein d’un cadre de référence en biologie systémique permettrait la prédiction de nouvelles fonctions de gènes qui demeurent non caractérisées à ce jour. Afin de réaliser de telles prédictions à l’échelle génomique chez la levure Saccharomyces cerevisiae, nous avons développé une stratégie innovatrice qui combine le criblage interactomique à haut débit des interactions protéines-protéines, la prédiction de la fonction des gènes in silico ainsi que la validation de ces prédictions avec la lipidomique à haut débit. D’abord, nous avons exécuté un dépistage à grande échelle des interactions protéines-protéines à l’aide de la complémentation de fragments protéiques. Cette méthode a permis de déceler des interactions in vivo entre les protéines exprimées par leurs promoteurs naturels. De plus, aucun biais lié aux interactions des membranes n’a pu être mis en évidence avec cette méthode, comparativement aux autres techniques existantes qui décèlent les interactions protéines-protéines. Conséquemment, nous avons découvert plusieurs nouvelles interactions et nous avons augmenté la couverture d’un interactome d’homéostasie lipidique dont la compréhension demeure encore incomplète à ce jour. Par la suite, nous avons appliqué un algorithme d’apprentissage afin d’identifier huit gènes non caractérisés ayant un rôle potentiel dans le métabolisme des lipides. Finalement, nous avons étudié si ces gènes et un groupe de régulateurs transcriptionnels distincts, non préalablement impliqués avec les lipides, avaient un rôle dans l’homéostasie des lipides. Dans ce but, nous avons analysé les lipidomes des délétions mutantes de gènes sélectionnés. Afin d’examiner une grande quantité de souches, nous avons développé une plateforme à haut débit pour le criblage lipidomique à contenu élevé des bibliothèques de levures mutantes. Cette plateforme consiste en la spectrométrie de masse à haute resolution Orbitrap et en un cadre de traitement des données dédié et supportant le phénotypage des lipides de centaines de mutations de Saccharomyces cerevisiae. Les méthodes expérimentales en lipidomiques ont confirmé les prédictions fonctionnelles en démontrant certaines différences au sein des phénotypes métaboliques lipidiques des délétions mutantes ayant une absence des gènes YBR141C et YJR015W, connus pour leur implication dans le métabolisme des lipides. Une altération du phénotype lipidique a également été observé pour une délétion mutante du facteur de transcription KAR4 qui n’avait pas été auparavant lié au métabolisme lipidique. Tous ces résultats démontrent qu’un processus qui intègre l’acquisition de nouvelles interactions moléculaires, la prédiction informatique des fonctions des gènes et une plateforme lipidomique innovatrice à haut débit , constitue un ajout important aux méthodologies existantes en biologie systémique. Les développements en méthodologies génomiques fonctionnelles et en technologies lipidomiques fournissent donc de nouveaux moyens pour étudier les réseaux biologiques des eucaryotes supérieurs, incluant les mammifères. Par conséquent, le stratégie présenté ici détient un potentiel d’application au sein d’organismes plus complexes.

Genetic association study of plasma creatine kinase levels in the Montreal Heart Institute Hospital Cohort

Il a déjà été démontré que les statines (ou inhibiteurs de la HMG-CoA réductase) sont efficaces pour réduire le LDL-cholestérol et elles se sont depuis établies comme étant le pilier dans le traitement de la dyslipidémie. Toutefois, environ 10 pourcent des utilisateurs de statines souffrent d'effets indésirables, généralement sous forme de myopathie qui est souvent accompagnée d’un taux élevé de la créatine kinase (CK) plasmatique. Il est fréquent que les patients doivent arrêter les statines à cause d’un taux de CK dépassant un seuil de référence. Nous avons examiné le taux de CK de près de 6000 participants de la biobanque de l’ICM, qui ont récemment été génotypés à l'aide de la micropuce d'ADN ExomChip d'Illumina. Des études antérieures ont démontré une association significative entre le taux de CK plasmatique et des polymorphismes génétiques et nous avons cherché à répliquer ces résultats par association génétique et à l'aide du test SKAT pour les polymorphismes rares. Nous avons répliqué les résultats dans le gène CKM (rs11559024, p=1.59x10-23) et le gène LILRB5 (rs12975366, p=1.44x10-26) dans le chromosome 19. Nous espérons que ces résultats seront éventuellement utilisés en clinique pour la prédiction des taux de référence de CK personnalisés selon le profil génétique des patients utilisateurs de statines.

Impact d'un traumatisme crânio-cérébral léger sur l’architecture du sommeil et le transcriptome dans un modèle murin

Le traumatisme crânien léger (TCL) est l'un des troubles neurologiques les plus courants affectant la santé publique. Aussi, les troubles du sommeil sont fréquents chez les patients atteints de TCL. Les études chez les rongeurs montrent que certains marqueurs de plasticité synaptique diminuent après le TCL, ce qui pourrait nuire à la plasticité du cerveau. Nous suggérons que la perte de sommeil intensifie l'effet négatif de TCL, qui peut refléter les changements des marqueurs de plasticité synaptique ou des changements des voies physiologiques qui régulent le sommeil. En utilisant un modèle de traumatisme crânien sur crâne fermé (closed head injury), nous avons étudié la relation bidirectionnelle entre le TCL et le sommeil en évaluant les effets de TCL sur l’activité électrique du cerveau par électroencéphalographie (EEG), et ceux de la privation de sommeil (PS) sur l'expression génique post-TCL. Premièrement, l'activité EEG a été enregistrée pour voir si l'architecture du sommeil est altérée suite au TCL. Nous avons ensuite voulu tester si la PS suite TCL induit des changements dans l'expression des gènes : Arc, Homer1a, Hif1a, Bdnf, Fos et éphrines, qui ont été liés à la plasticité synaptique et à la régulation du sommeil. Nous avons également étudié l'effet de la PS post-TCL sur le génome complet dans les régions cibles (cortex et l'hippocampe). Les principaux résultats obtenus dans cette étude confirment que TCL modifie de manière significative l'activité spectrale pendant l'éveil, le sommeil Rapid Eye Movement (REM) et le sommeil non-REM dans le deuxième 24 heures post-TCL. Fait intéressant, la capacité de maintenir de longues périodes d'éveil a été altérée immédiatement après TCL (première 24h post-TCL). La dynamique de l'activité delta pendant l'éveil a été modifié par le TCL. Parallèlement à ces modifications, des changements dans l'expression des gènes ont été observés dans le cortex et l'hippocampe. Seulement Arc et EfnA3 ont montré une interaction TCL / PS et ce dans l’hippocampe, tandis que l'expression de tous les autres gènes semblait être affectée par la PS ou TCL indépendamment. Nos résultats montrent pour la première fois que le TCL induit l'expression de deux chimiokines (Ccl3 et Cxcl5) à la fois dans le cortex cérébral et l'hippocampe 2,5 jours post-TCL. Également, nous avons observé que le TCL induit une diminution de l'expression de Lgals3 et S100A8 dans le cortex, et une augmentation d’Olig2 dans l'hippocampe. Les résultats concernant les effets de la PS sur le génome complet du cortex et de l'hippocampe montrent des changements significatifs dans les gènes impliqués dans diverses fonctions physiologiques, telles que les rythmes circadiens, la réponse inflammatoire, ainsi que de l'activation des cellules gliales. En général, nos résultats précisent les changements dans la qualité de l’éveil ainsi que dans l'expression de divers gènes après TCL.

Effet du stress prolifératif sur la fonction des cellules souches hématopoïétiques : rôles des gènes Scl, E2A et Heb

Le système hématopoïétique est un tissu en constant renouvellement et les cellules souches hématopoïétiques (CSHs) sont indispensables pour soutenir la production des cellules matures du sang. Deux fonctions définissent les CSHs; la propriété d’auto-renouvellement, soit la capacité de préserver l’identité cellulaire suivant une division, et la multipotence, le potentiel de différenciation permettant de générer toutes les lignées hématopoïétiques. Chez l’adulte, la majorité des CSHs sont quiescentes et l’altération de cet état corrèle avec une diminution du potentiel de reconstitution des CSHs, suggérant que la quiescence protège les fonctions des CSHs. La quiescence est un état réversible et dynamique et les réseaux génétiques le contrôlant restent peu connus. Un nombre croissant d’évidences suggère que si à l’état d’homéostasie il y a une certaine redondance entre les gènes impliqués dans ces réseaux de contrôle, leurs rôles spécifiques sont révélés en situation de stress. La famille des bHLHs (basic helix-loop-helix) inclue différentes classes des protéines dont ceux qui sont tissu-spécifiques comme SCL, et les protéines E, comme E12/E47 et HEB. Certains bHLHs sont proposés êtres important pour la fonction des cellules souches, mais cela ne fait pas l’unanimité, car selon le contexte cellulaire, il y a redondance entre ces facteurs. La question reste donc entière, y a-t-il un rôle redondant entre les bHLHs d’une même classe pour la fonction à long-terme des CSHs? Les travaux présentés dans cette thèse visaient dans un premier temps à explorer le lien encore mal compris entre la quiescence et la fonction des CSHs en mesurant leurs facultés suite à un stress prolifératif intense et dans un deuxième temps, investiguer l’importance et la spécificité de trois gènes pour la fonction des CSHs adultes, soit Scl/Tal1, E2a/Tcf3 et Heb/Tcf12. Pour répondre à ces questions, une approche cellulaire (stress prolifératif) a été combinée avec une approche génétique (invalidation génique). Plus précisément, la résistance des CSHs au stress prolifératif a été étudiée en utilisant deux tests fonctionnels quantitatifs optimisés, soit un traitement basé sur le 5-fluorouracil, une drogue de chimiothérapie, et la transplantation sérielle en nombre limite. Dans la mesure où la fonction d’un réseau génique ne peut être révélée que par une perturbation intrinsèque, trois modèles de souris, i.e. Scl+/-, E2a+/- et Heb+/- ont été utilisés. Ceci a permis de révéler que l’adaptation des CSHs au stress prolifératif et le retour à l’équilibre est strictement contrôlé par les niveaux de Scl, lesquels règlent le métabolisme cellulaire des CSHs en maintenant l’expression de gènes ribosomaux à un niveau basal. D’autre part, bien que les composantes du réseau puissent paraître redondants à l’équilibre, mes travaux montrent qu’en situation de stress prolifératif, les niveaux de Heb restreignent la prolifération excessive des CSHs en induisant la sénescence et que cette fonction ne peut pas être compensée par E2a. En conclusion, les résultats présentés dans cette thèse montrent que les CSHs peuvent tolérer un stress prolifératif intense ainsi que des dommages à l’ADN non-réparés, tout en maintenant leur capacité de reconstituer l’hématopoïèse à long-terme. Cela implique cependant que leur métabolisme revienne au niveau de base, soit celui trouvé à l’état d’homéostasie. Par contre, avec l’augmentation du nombre de division cellulaire les CSHs atteignent éventuellement une limite d’expansion et entrent en sénescence.

Étude de la régulation de l'activité du ligand Delta dans le cadre de la signalisation Notch

La voie de signalisation Notch est conservée au cours de l'évolution. Elle joue un rôle clé dans le développement, et elle est impliquée dans de nombreuses décisions de destin cellulaire, dans le maintien des cellules souches, et dans le contrôle de la prolifération et de la différenciation cellulaires. Une dérégulation de la signalisation Notch est impliquée dans diverses maladies et cancers, y compris les tumeurs solides, comme les cancers du sein et du col de l'utérus, et les leucémies, comme la Leucémie Aiguë Lymphoblastique des cellules T (LAL-T). Notch est un récepteur transmembranaire activé par des ligands transmembranaires de la famille DSL (Delta/Serrate/Lag-2). Bien que plusieurs mutations oncogéniques ont été identifiées au niveau du récepteur Notch, de nombreux cancers modulés par Notch demeurent ligand-dépendants. Étonnamment, les mécanismes moléculaires régulant l'activation du ligand sont encore relativement peu caractérisés par rapport à ceux qui régissent le récepteur Notch lui-même. Utilisant un essai de co-culture avec un rapporteur luciférase de Notch, nous avons effectué le premier crible d'ARNi pan-génomique visant spécifiquement à identifier des régulateurs des ligands de Notch dans la cellule émettrice du signal. Nous avons ainsi pu découvrir de nouvelles classes de régulateurs communs pour les ligands Delta-like1 et 4. Ces régulateurs comprennent des inhibiteurs de protéases, des facteurs de transcription, et des gènes divers à fonction inconnue, tels que Tmem128 « Transmembrane protein 128 », ou à fonction préalablement caractérisée tels que la co-chaperonne moléculaire Cdc37 « Cell division cycle 37 homolog ». Par la suite, nous avons développé des cribles secondaires fonctionnels où nous avons démontré l'importance de ces régulateurs pour des événements Notch-dépendants, comme la différenciation des cellules T normales, et la survie des cellules souches pré-leucémiques isolées à partir d'un modèle murin de LAL-T. En outre, nous avons prouvé que les régulateurs les plus forts du crible de survie sont également nécessaires pour l'activité d'auto-renouvellement des cellules souches pré-leucémiques. Finalement, nous avons entamé une caractérisation moléculaire préliminaire de deux régulateurs nouvellement identifiés; Tmem128 et Cdc37 afin d'étudier leur mécanisme d'action sur les ligands. En conclusion, cette étude nous a permis d'identifier de nouveaux régulateurs de la voie Notch qui pourraient servir de cibles thérapeutiques potentielles dans les cancers; tel qu'illustré par le modèle LAL-T. La compréhension des détails moléculaires sous-jacents aux fonctions de ces régulateurs sera essentielle afin de développer des inhibiteurs pharmacologiques pour bloquer leur action et entraver la signalisation Notch dans le cancer.

Identification et étude de mécanismes régulant l’expression de MAPK

Les fichiers accompagnant le document sont en format Microsoft Excel 2010.

Régulation de la division asymétrique chez C. elegans

La division asymétrique est essentielle pour générer la diversité au cours du développement et permet aussi de réguler la balance entre renouvellement et différenciation des cellules souches chez l’adulte. Dans ces deux cas de figure, elle dépend respectivement d’une polarité intrinsèque ou d’une polarité extrinsèque. C. elegans est un excellent modèle pour étudier les mécanismes cellulaires et moléculaires de la division asymétrique in vivo. Chez l’embryon, le maintien d’un axe de polarité antéro-postérieur dépend des protéines PAR conservées et localisées de façon asymétrique en deux groupes mutuellement exclusifs; le groupe antérieur avec PAR-3, PAR-6, PKC-3 et le groupe postérieur avec PAR-2 et PAR-1. L’absence d’une protéine PAR entraine une perte de polarité et une létalité embryonnaire. Lors d’un crible par ARN interférence mené par Jean-Claude Labbé pour identifier les suppresseurs de la létalité associée à la perte de PAR-2, deux cyclines de type B, CYB-2.1 et CYB-2.2 ont été trouvées. J’ai déterminé que CYB-2.1 et CYB-2.2 interviennent dans la polarité sans perturber le cycle cellulaire et agissent vraisemblablement avec leur kinase associée, CDK-1, pour stabiliser les niveaux protéiques de PAR-6. Ces travaux permettent de mieux définir les liens étroits entre polarité et cycle cellulaire. La lignée germinale de C. elegans est un excellent modèle pour étudier les divisions des cellules souches germinales in vivo. Par contre, l’absence d’orientation préférentielle de ces divisions laisse envisager que la complexité morphologique de la niche pourrait engendrer une diversité d’axe possible. J’ai étudié la régulation morphologique de cette niche, une unique cellule somatique appelée distal tip cell (DTC), qui arborise de longues extensions au stade adulte. Mes résultats préliminaires favorisent un modèle dans lequel les cellules souches et progéniteurs germinaux (CSPG) supportent la formation de ces extensions. Enfin, j’ai obtenu des conditions favorables à l’étude de la division asymétrique extrinsèque dans ce modèle, en simplifiant l’architecture de la niche dans des conditions qui préservent les divisions cellulaires des cellules souches. Mes travaux ont permis de mieux comprendre les liens unissant les différents processus biologiques impliqués dans la division asymétrique, d’une part par l’étude du rôle qu’y jouent des régulateurs clés du cycle cellulaire au cours du développement et d’autre part par la caractérisation d’une communication bidirectionnelle entre la niche et les cellules souches chez l’adulte.

Identification of Genomic Variants Associated with Adolescent Idiopathic Scoliosis (AIS) in French-Canadian Population

La scoliose idiopathique est une déformation tridimensionnelle de la colonne vertébrale dont la pathogenèse reste obscure. Cette maladie affecte 2-4% des adolescents de 10-18 ans parmi les garçons et les filles. Il est à noter que les filles sont plus sévèrement affectées et ce en plus grand nombre que les garçons. Les études de jumeaux ont montré que les facteurs génétiques jouent un rôle important dans la scoliose idiopathique de l'adolescent (SIA). Depuis 2010, les études d'association pan génomiques ont été multipliées dans les recherches, visant à trouver des gènes candidats impliqués dans la SIA à travers des examens des polymorphismes nucléotidiques (SNPs). Un test génétique nommé "ScoliScore" a été publié pour essayer de prédire la progression de courbure dans la population caucasienne. Cependant, l'association n'a pas été reproduite dans une grande étude japonaise, soulignant l'importance d'une étude de réplication dans une population caucasienne indépendante. Dans ce contexte, mon projet de maîtrise a permis de génotyper plus de 1,4 millions de SNPs dans une cohorte canadienne-française dans le but: 1) de valider l'association de ScoliScoreTM; et 2) d’identifier les variants génomiques associées à la SIA dans la population québécoise. Notre étude a montré qu’aucun des variants constituant le test ScoliScoreTM n’était associé à la SIA. Ceci suggère que l'absence d'association dans une cohorte japonaise n'est pas due à l'appartenance ethnique. Aussi, nous avons identifié des variants génomiques associés significativement à l’initiation et/ou la progression de SIA dans la population québécoise, suggérant des gènes candidats impliqués dans la pathogenèse de SIA.

Développement d’une méthode bio-informatique permettant de relier les gènes aux métabolites

L’objectif de ce projet était de faire le lien entre gènes et métabolites afin d’éventuellement proposer des métabolites à mesurer en lien avec la fonction de gènes. Plus particulièrement, nous nous sommes intéressés aux gènes codant pour des protéines ayant un impact sur le métabolisme, soit les enzymes qui catalysent les réactions faisant partie intégrante des voies métaboliques. Afin de quantifier ce lien, nous avons développé une méthode bio-informatique permettant de calculer la distance qui est définie comme le nombre de réactions entre l’enzyme encodée par le gène et le métabolite dans la carte globale du métabolisme de la base de données Kyoto Encyclopedia of Genes and Genomes (KEGG). Notre hypothèse était que les métabolites d’intérêt sont des substrats/produits se trouvant à proximité des réactions catalysées par l’enzyme encodée par le gène. Afin de tester cette hypothèse et de valider la méthode, nous avons utilisé les études d’association pangénomique combinées à la métabolomique (mGWAS) car elles rapportent des associations entre variants génétiques, annotés en gènes, et métabolites mesurés. Plus précisément, la méthode a été appliquée à l’étude mGWAS par Shin et al. Bien que la couverture des associations de Shin et al. était limitée (24/299), nous avons pu valider de façon significative la proximité entre gènes et métabolites associés (P<0,01). En somme, cette méthode et ses développements futurs permettront d’interpréter de façon quantitative les associations mGWAS, de prédire quels métabolites mesurer en lien avec la fonction d’un gène et, plus généralement, de permettre une meilleure compréhension du contrôle génétique sur le métabolisme.