Using whole-exome sequencing data in an exome-wide association study approach to identify genetic risk factors influencing acute lymphoblastic leukemia response : a focus on asparaginase complications & vincristine-induced peripheral neuropathy

Abstract

Le traitement de la leucémie lymphoblastique aiguë (LLA) de l’enfant, une affection d'origine maligne des cellules progénitrices lymphoïdes, s’est considérablement amélioré au cours des dernières décennies. En effet, le taux de succès du traitement a dépassé 90% dans des conditions favorables. Cependant, des toxicités liées au traitement peuvent être fatales et entrainer l’interruption ou la cessation du traitement. L'allergie, la pancréatite et la thrombose sont des complications fréquentes du traitement de la LLA et sont associées à l'utilisation de l'asparaginase (ASNase), tandis qu’une toxicité fréquente due à la vincristine (VCR) induit la neuropathie périphérique (VIPN). Étant donné que l’ajustement du schéma posologique afin d’augmenter l'efficacité et diminuer la toxicité est un processus sensible, ceci demeure un défi majeur dans plusieurs protocoles de traitement. La pharmacogénétique étudie comment des altérations de la composante génétique peuvent influer sur la variabilité interindividuelle observée dans la réponse au traitement. Une meilleure compréhension de la base moléculaire de cette variabilité pourrait améliorer considérablement les résultats du traitement, en permettant la personnalisation de ce dernier en fonction du profil génétique du patient. Des études récentes suggèrent l’avantage d’appliquer l’analyse de l’exome à la découverte de variants associés à des traits humains complexes ainsi qu’à des phénotypes de réactions médicamenteuses. L'objectif de notre travail était d'utiliser les données de séquençage pour réaliser des études d'association à l'échelle de l'exome, y compris des étapes de filtrage et de validation, afin d'identifier de nouveaux variants génétiques susceptibles de moduler le risque de développer des complications associées à ASNase et à VIPN. Douze SNP étaient associés à des complications due à l’ASNase dans la cohorte initiale, dont 3 étaient associés à une allergie, 3 à une pancréatite et 6 à une thrombose. Parmi ceux-ci, les variants rs3809849, rs11556218 et rs34708521 des gènes MYBBP1A, IL16 et SPEF2 respectivement ont été associés à des complications multiples et leur association à une pancréatite a été répliquée dans une cohorte de validation indépendante. En ce qui concerne la VCR, trois variantes ont été associées à la modulation du risque de VIPN: rs2781377 dans SYNE2, rs10513762 dans MRPL47 et rs3803357 dans BAHD1. Nous démontrons également le puissant effet combiné de la présence de plusieurs variants de risque pour chacune des toxicités étudiées et fournissons des modèles de prédiction du risque pour la pancréatite et le VIPN basés sur la méthode d’évaluation du risque génétique pondérée et qui ont été validés à l’interne. De plus, étant donné une association du polymorphisme du gène MYBBP1A avec de multiples issus de traitement, nous avons cherché à comprendre comment cette altération génétique se traduit par des variabilités de réponse aux traitements à l’ASNase. En utilisant la technique CRISPR-CAS9 pour induire l'inactivation de gènes dans des lignées cellulaires cancéreuses PANC1 (pancréatiques) nous avons testé la différence de viabilité entre les cellules inactivées et les cellules du type sauvage à la suite de la suppression du gène et du traitement par ASNase. Nos résultats suggèrent un rôle fonctionnel de ce gène dans la modulation de la viabilité, de la capacité de prolifération et de la morphologie des cellules knock-out, ainsi que dans leur sensibilité à l'ASNase, et plaident en outre pour que le gène influence l’issus du traitement de la LLA par ASNase. Le présent travail démontre que l’utilisation de l’approche de séquençage de l’exome entier dans le contexte d’une étude d’association à l’échelle de l’exome est une stratégie valide « sans hypothèse » pour identifier de nouveaux marqueurs génétiques modulant l’effet du traitement de la LLA de l’enfant, et souligne l’importance de l'effet synergique de la combinaison des locus à risque.

Treatment of childhood acute lymphoblastic leukemia (ALL), a malignant disorder of lymphoid progenitor cells has improved significantly over the past decades and treatment success rates have surpassed 90% in favorable settings. However, treatment-related toxicities can be life-threatening and cause treatment interruption or cessation. Allergy, pancreatitis and thrombosis are common complications of ALL treatment associated with the use of asparaginase (ASNase), while vincristine-induced peripheral neuropathy (VIPN) is a frequent toxicity of vincristine (VCR). It is a sensitive process and a constant struggle to adjust the dosing regimen to ensure maximum efficacy and minimum toxicity. Pharmacogenetics studies show alterations in the genetic component between individuals can influence the observed variability in treatment response. A better understanding of the molecular basis of this variability in drug effect could significantly improve treatment outcome by allowing the personalization of ALL treatment based on the genetic profile of the patient. Emerging reports suggest the benefit of applying exome analysis to uncover variants associated with complex human traits as well as drug response phenotypes. Our objective in this work was to use available whole-exome sequencing data to perform exome-wide association studies followed by stepwise filtering and validation processes to identify novel variants with a potential to modulate the risk of developing ASNase complications and VIPN. Twelve SNPs were associated with ASNase complications in the discovery cohort including 3 associated with allergy, 3 with pancreatitis and 6 with thrombosis. Of those, rs3809849 in MYBBP1A, rs11556218 in IL16 and rs34708521 in SPEF2 genes were associated with multiple complications and their association with pancreatitis was replicated in an independent validation cohort. As for VCR, three variants were associated with modulating the risk of VIPN: rs2781377 in SYNE2, rs10513762 in MRPL47 and rs3803357 in BAHD1. We also demonstrate a strong combined effect of harbouring multiple risk variants for each of the studied toxicities, and provide internally-validated risk-prediction models based on the weighted genetic risk score method for pancreatitis and VIPN. Furthermore, given the association of the polymorphism in MYBBP1A gene with multiple treatment outcomes, we aimed at understanding how this genetic alteration translates into differences in ASNase treatment response through cell-based functional analysis. Using CRISPR-CAS9 technology we produced gene knockout of PANC1 (pancreatic) cancer cell-lines and tested the difference in viability between the knockouts and wild-type cells following gene deletion and ASNase treatment. Our results suggest a functional role of this gene in modulating the viability, proliferation capacity and the morphology of the knockout cells as well as their sensitivity to ASNase and further advocates the implication of the gene in influencing the outcome of ALL treatment with ASNase. The present work demonstrates that using whole-exome sequencing data in the context of exome-wide association study is a successful “hypothesis-free” strategy for identifying novel genetic markers modulating the effect of childhood ALL treatment and highlights the importance of the synergistic effect of combining risk loci.