In silico analysis of mitochondrial proteins

Le rôle important joué par la mitochondrie dans la cellule eucaryote est admis depuis longtemps. Cependant, la composition exacte des mitochondries, ainsi que les processus biologiques qui sy déroulent restent encore largement inconnus. Deux facteurs principaux permettent dexpliquer pourquoi létude des mitochondries progresse si lentement : le manque defficacité des méthodes didentification des protéines mitochondriales et le manque de précision dans lannotation de ces protéines. En conséquence, nous avons développé un nouvel outil informatique, YimLoc, qui permet de prédire avec succès les protéines mitochondriales à partir des séquences génomiques. Cet outil intègre plusieurs indicateurs existants, et sa performance est supérieure à celle des indicateurs considérés individuellement. Nous avons analysé environ 60 génomes fongiques avec YimLoc afin de lever la controverse concernant la localisation de la bêta-oxydation dans ces organismes. Contrairement à ce qui était généralement admis, nos résultats montrent que la plupart des groupes de Fungi possèdent une bêta-oxydation mitochondriale. Ce travail met également en évidence la diversité des processus de bêta-oxydation chez les champignons, en corrélation avec leur utilisation des acides gras comme source dénergie et de carbone. De plus, nous avons étudié le composant clef de la voie de bêta-oxydation mitochondriale, lacyl-CoA déshydrogénase (ACAD), dans 250 espèces, couvrant les 3 domaines de la vie, en combinant la prédiction de la localisation subcellulaire avec la classification en sous-familles et linférence phylogénétique. Notre étude suggère que les gènes ACAD font partie dune ancienne famille qui a adopté des stratégies évolutionnaires innovatrices afin de générer un large ensemble denzymes susceptibles dutiliser la plupart des acides gras et des acides aminés. Finalement, afin de permettre la prédiction de protéines mitochondriales à partir de données autres que les séquences génomiques, nous avons développé le logiciel TESTLoc qui utilise comme données des Expressed Sequence Tags (ESTs). La performance de TESTLoc est significativement supérieure à celle de tout autre outil de prédiction connu. En plus de fournir deux nouveaux outils de prédiction de la localisation subcellulaire utilisant différents types de données, nos travaux démontrent comment lassociation de la prédiction de la localisation subcellulaire à dautres méthodes danalyse in silico permet daméliorer la connaissance des protéines mitochondriales. De plus, ces travaux proposent des hypothèses claires et faciles à vérifier par des expériences, ce qui présente un grand potentiel pour faire progresser nos connaissances des métabolismes mitochondriaux.

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.

Cell lines and animal model in the analysis of pharmacogenomics markers in childhood acute lymphoblastic leukemia

La leucémie aiguë lymphoblastique (LAL) est le cancer pédiatrique le plus fréquent. Elle est la cause principale de mortalité liée au cancer chez les enfants due à un groupe de patient ne répondant pas au traitement. Les patients peuvent aussi souffrir de plusieurs toxicités associées à un traitement intensif de chimiothérapie. Les études en pharmacogénétique de notre groupe ont montré une corrélation tant individuelle que combinée entre les variants génétiques particuliers d’enzymes dépendantes du folate, particulièrement la dihydrofolate réductase (DHFR) ainsi que la thymidylate synthase (TS), principales cibles du méthotrexate (MTX) et le risque élevé de rechute chez les patients atteints de la LAL. En outre, des variations dans le gène ATF5 impliqué dans la régulation de l’asparagine synthetase (ASNS) sont associées à un risque plus élevé de rechute ou à une toxicité ASNase dépendante chez les patients ayant reçu de l’asparaginase d’E.coli (ASNase). Le but principal de mon projet de thèse est de comprendre davantage d’un point de vue fonctionnel, le rôle de variations génétiques dans la réponse thérapeutique chez les patients atteints de la LAL, en se concentrant sur deux composants majeurs du traitement de la LAL soit le MTX ainsi que l’ASNase. Mon objectif spécifique était d’analyser une association trouvée dans des paramètres cliniques par le biais d’essais de prolifération cellulaire de lignées cellulaires lymphoblastoïdes (LCLs, n=93) et d’un modèle murin de xénogreffe de la LAL. Une variation génétique dans le polymorphisme TS (homozygosité de l’allèle de la répétition triple 3R) ainsi que l’haplotype *1b de DHFR (défini par une combinaison particulière d’allèle dérivé de six sites polymorphiques dans le promoteur majeur et mineur de DHFR) et de leurs effets sur la sensibilité au MTX ont été évalués par le biais d’essais de prolifération cellulaire. Des essais in vitro similaires sur la réponse à l’ASNase de E. Coli ont permis d’évaluer l’effet de la variation T1562C de la région 5’UTR de ATF5 ainsi que des haplotypes particuliers du gène ASNS (définis par deux variations génétiques et arbitrairement appelés haplotype *1). Le modèle murin de xénogreffe ont été utilisé pour évaluer l’effet du génotype 3R3R du gène TS. L’analyse de polymorphismes additionnels dans le gène ASNS a révélé une diversification de l’haplotype *1 en 5 sous-types définis par deux polymorphismes (rs10486009 et rs6971012,) et corrélé avec la sensibilité in vitro à l’ASNase et l’un d’eux (rs10486009) semble particulièrement important dans la réduction de la sensibilité in vitro à l’ASNase, pouvant expliquer une sensibilité réduite de l’haplotype *1 dans des paramètres cliniques. Aucune association entre ATF5 T1562C et des essais de prolifération cellulaire en réponse à ASNase de E.Coli n’a été détectée. Nous n’avons pas détecté une association liée au génotype lors d’analyse in vitro de sensibilité au MTX. Par contre, des résultats in vivo issus de modèle murin de xénogreffe ont montré une relation entre le génotype TS 3R/3R et la résistance de manière dose-dépendante au traitement par MTX. Les résultats obtenus ont permis de fournir une explication concernant un haut risque significatif de rechute rencontré chez les patients au génotype TS 3R/3R et suggèrent que ces patients pourraient recevoir une augmentation de leur dose de MTX. À travers ces expériences, nous avons aussi démontré que les modèles murins de xénogreffe peuvent servir comme outil préclinique afin d’explorer l’option d’un traitement individualisé. En conclusion, la connaissance acquise à travers mon projet de thèse a permis de confirmer et/ou d’identifier quelques variants dans la voix d’action du MTX et de l’ASNase qui pourraient faciliter la mise en place de stratégies d’individualisation de la dose, permettant la sélection d’un traitement optimum ou moduler la thérapie basé sur la génétique individuelle.

Anatomic and molecular characterization of the 1 endocrine pancreas of a teleostean fish: Atlantic wolffish (Anarhichas lupus)

Background: The biologic attributes of the endocrine pancreas and the comparative endocrinology of islet amyloid polypeptide (IAPP) of fish are not well described in the literature. This study describes the endocrine pancreas of one teleostean fish. Ten captive Atlantic wolffish (Anarhichas lupus) from the Montreal Biodome were submitted for necropsy and their pancreata were collected. Results: Grossly, all the fish pancreata examined contained 1-3 nodules of variable diameter (1-8 mm). Microscopically, the nodules were uniform, highly cellular, and composed of polygonal to elongated cells. Immunofluorescence for pancreatic hormones was performed. The nodules were immunoreactive for insulin most prominent centrally, but with IAPP and glucagon only in the periphery of the nodules. Exocrine pancreas was positive for chromogranin A. Not previously recognized in fish, IAPP immunoreactivity occurred in α, glucagon-containing, cells and did not co-localize with insulin in β cells. The islet tissues were devoid of amyloid deposits. IAPP DNA sequencing was performed to compare the sequence among teleost fish and the potency to form amyloid fibrils. In silico analysis of the amino acid sequences 19-34 revealed that it was not amyloidogenic. Conclusions: Amyloidosis of pancreatic islets would not be expected as a spontaneous disease in the Atlantic wolffish. Our study underlines that this teleost fish is a potential candidate for pancreatic xenograft research.

Identification and characterization of new biomarkers in aggressive subtypes of breast cancer

En 2015, la récidive tumorale et les métastases du cancer du sein demeurent une cause importante de décès à travers le monde. Toutefois, ces cancers sont souvent hétérogènes car en dépit d’un phénotype similaire, l’évolution clinique et la réponse au traitement peuvent varier considérablement. Il y a donc un intérêt évident à identifier et à caractériser de nouveaux biomarqueurs pour permettre classer les tumeurs mammaires dans des sous-groupes plus homogènes. Notre hypothèse est que chaque cancer mammaire possède des caractéristiques distinctes au plan des altérations du génome et des profils d’expression géniques et que ces changements se traduisent cliniquement par une prédisposition à former des métastases ou à répondre ou non à la chimiothérapie et aux thérapies ciblées. Dans le cadre de nos travaux, nous nous sommes intéressés aux sous-types agressifs de tumeurs mammaires et notamment les cancers de type triple négatif. Nous avons aussi tenté d’identifier des marqueurs capables de distinguer l’une de l’autre les tumeurs de type luminal A et luminal B. Pour ce faire, nous avons d’abord utilisé une stratégie in silico à partir de données publiques (micro-puces d’ADN et séquençage de l’ARN). Nous avons ensuite construit sept micro-matrices tissulaires (TMA) provenant de tissus mammaires normaux et tumoraux fixés à la formaline et enrobés en paraffine. Ces outils nous ont permis d’évaluer par immunohistochimie les niveaux d’expression différentielle des marqueurs suivants : ANXA1, MMP-9, DP103 et MCM2. Ceux-ci ont été comparés aux marqueurs usuels du cancer du sein (ER, PR, HER2, CK5/6 et FOXA1) et corrélés aux données cliniques (survie globale et métastase). Nos résultats indiquent que ces nouveaux marqueurs jouent un rôle important dans l’évolution clinique défavorable des tumeurs de haut grade. Dans un premier article nous avons montré que l’expression d’ANXA1 est dérégulée dans les cancers de type triple-négatif et aussi, dans une certaine mesure, dans les tumeurs HER2+. Nous croyons qu’ANXA1 permet de mieux comprendre le processus d’hétérogénéité tumorale et facilite l’identification des tumeurs de haut grade. Nous proposons également qu’ d’ANXA1 stimule la transition épithélio-mésenchymateuse (EMT) et la formation des métastases. Dans un second temps, nous avons montré que les niveaux d’expression de MMP-9 reflètent la différenciation cellulaire et corrèlent avec les sous-types de cancers mammaires ayant un mauvais pronostic. Nous estimons que MMP-9 permet de mieux comprendre et d’identifier les tumeurs mammaires à haut risque. De fait, la surexpression de MMP-9 est associée à une augmentation des métastases, une récidive précoce et une diminution de la survie globale. Dans le cadre d’un troisième article, nous avons montré que la surexpression du marqueur de prolifération MCM2 s’observe dans les cancers triple-négatifs, HER2+ et Luminal B par comparaison aux cancers luminal A (p< 0.0001). Nos résultats suggèrent qu’en utilisant un seuil de 40% de noyaux marqués, nous pourrions distinguer l’une de l’autre les tumeurs de type luminal A et luminal B. Cela dit, avant de pouvoir envisager l’utilisation de ce marqueur en clinique, une étude de validation sur une nouvelle cohorte de patientes s’impose. En somme, les résultats de nos travaux suggèrent qu’ANXA1, MMP-9 et MCM2 sont des marqueurs intéressants pour mieux comprendre les mécanismes physiopathologiques impliqués dans la progression tumorale et le développement des métastases. À terme, ces nouveaux marqueurs pourraient être utilisés seuls ou en combinaison avec d’autres gènes candidats pour permettre le développement de trousses « multigènes » ou d’essais protéomiques multiplex pour prédire l’évolution clinique des cancers mammaires.

Structural aspects of the ribosome evolution and function

Les résultats ont été obtenus avec le logiciel "Insight-2" de Accelris (San Diego, CA)

Modulation de l’expression du gène CFTR par le produit du gène FIC1 responsable de la cholestase familiale intra-hépatique progressive de type 1 : Identification des mécanismes moléculaires impliqués

Réalisé en cotutelle avec l'Université de Cergy-Pontoise

Annotation des ARN non codants du génome de Candida albicans par méthode bioinformatique

La bio-informatique est un champ pluridisciplinaire qui utilise la biologie, l’informatique, la physique et les mathématiques pour résoudre des problèmes posés par la biologie. L’une des thématiques de la bio-informatique est l’analyse des séquences génomiques et la prédiction de gènes d’ARN non codants. Les ARN non codants sont des molécules d’ARN qui sont transcrites mais pas traduites en protéine et qui ont une fonction dans la cellule. Trouver des gènes d’ARN non codants par des techniques de biochimie et de biologie moléculaire est assez difficile et relativement coûteux. Ainsi, la prédiction des gènes d’ARNnc par des méthodes bio-informatiques est un enjeu important. Cette recherche décrit un travail d’analyse informatique pour chercher des nouveaux ARNnc chez le pathogène Candida albicans et d’une validation expérimentale. Nous avons utilisé comme stratégie une analyse informatique combinant plusieurs logiciels d’identification d’ARNnc. Nous avons validé un sous-ensemble des prédictions informatiques avec une expérience de puces à ADN couvrant 1979 régions du génome. Grace à cette expérience nous avons identifié 62 nouveaux transcrits chez Candida albicans. Ce travail aussi permit le développement d’une méthode d’analyse pour des puces à ADN de type tiling array. Ce travail présente également une tentation d’améliorer de la prédiction d’ARNnc avec une méthode se basant sur la recherche de motifs d’ARN dans les séquences.

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

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.

Understanding the transcriptional control of EIF4E and its dysregulation in acute myeloid leukemia: role of NF-κB

EIF4E, le facteur d’initiation de la traduction chez les eucaryotes est un oncogène puissant et qui se trouve induit dans plusieurs types de cancers, parmi lesquels les sous-types M4 et M5 de la leucémie aiguë myéloblastique (LAM). EIF4E est régulé à plusieurs niveaux cependant, la régulation transcriptionnelle de ce gène est peu connue. Mes résultats montrent que EIF4E est une cible transcriptionnelle directe du facteur nucléaire « kappa-light- chain- enhancer of activated B cells » (NF-κB).Dans les cellules hématopoïétiques primaires et les lignées cellulaires, les niveaux de EIF4E sont induits par des inducteurs de NF-κB. En effet, l’inactivation pharmaceutique ou génétique de NF-κB réprime l’activation de EIF4E. En effet, suite à l’activation de NF-κB chez l’humain, le promoteur endogène de EIF4E recrute p65 (RelA) et c-Rel aux sites évolutionnaires conservés κB in vitro et in vivo en même temps que p300 ainsi que la forme phosphorylée de Pol II. De plus, p65 est sélectivement associé au promoteur de EIF4E dans les sous-types LAM M4/M5 mais non pas dans les autres sous-types LAM ou dans les cellules hématopoïétiques primaires normales. Ceci indique que ce processus représente un facteur essentiel qui détermine l’expression différentielle de EIF4E dans la LAM. Les analyses de données d’expressions par séquençage de l’ARN provenant du « Cancer Genome Atlas » (TCGA) suggèrent que les niveaux d’ARNm de EIF4E et RELA se trouvent augmentés dans les cas LAM à pronostic intermédiaire ou faible mais non pas dans les groupes cytogénétiquement favorables. De plus, des niveaux élevés d’ARNm de EIF4E et RELA sont significativement associés avec un taux de survie relativement bas chez les patients. En effet, les sites uniques κB se trouvant dans le promoteur de EIF4E recrutent le régulateur de transcription NF-κB p65 dans 47 nouvelles cibles prévues. Finalement, 6 nouveaux facteurs de transcription potentiellement impliqués dans la régulation du gène EIF4E ont été prédits par des analyses de données ChIP-Seq provenant de l’encyclopédie des éléments d’ADN (ENCODE). Collectivement, ces résultats fournissent de nouveaux aperçus sur le control transcriptionnel de EIF4E et offrent une nouvelle base moléculaire pour sa dérégulation dans au moins un sous-groupe de spécimens de LAM. L’étude et la compréhension de ce niveau de régulation dans le contexte de spécimens de patients s’avère important pour le développement de nouvelles stratégies thérapeutiques ciblant l’expression du gène EIF4E moyennant des inhibiteurs de NF-κB en combinaison avec la ribavirine.