Classes Alléliques d’Haplotypes et Sélection Positive dans le Génome Humain

L'identification de régions génomiques cibles de la sélection naturelle positive permet de mieux comprendre notre passé évolutif et de trouver des variants génétiques fonctionnels importants. Puisque la fréquence des allèles sélectionnés augmente dans la population, la sélection laisse des traces sur les séquences d'ADN et ces empreintes sont détectées lorsque la variabilité génétique d'une région est différente de celle attendue sous neutralité sélective. On propose une nouvelle approche pour analyser les données de polymorphismes : le calcul des classes alléliques d’haplotypes (HAC), permettant d'évaluer la diversité globale des haplotypes en étudiant leur composition allélique. L'idée de l'approche est de déterminer si un site est sous sélection positive récente en comparant les distributions des HAC obtenues pour les deux allèles de ce site. Grâce à l'utilisation de données simulées, nous avons étudié ces distributions sous neutralité et sous sélection en testant l'effet de différents paramètres populationnels. Pour tester notre approche empiriquement, nous avons analysé la variation génétique au niveau du gène de lactase dans les trois populations inclues dans le projet HapMap.

Étude du polymorphisme du gène majeur d’histocompatibilité de classe IIb (MHIIb) chez l’omble de fontaine (Salvelinus fontinalis)

Les molécules classiques du complexe majeur d’histocompatibilité de classe II (CMHII) sont des glycoprotéines de surface spécialisées dans la présentation de peptides, principalement dérivés de pathogènes extracellulaires, aux récepteurs des lymphocytes T CD4+ afin d’initier la réponse immunitaire adaptative. Elles sont encodées, avec celles du CMH de classe I, par les gènes les plus polymorphiques identifiés jusqu’à maintenant, avec plusieurs loci et une grande diversité allélique à chacun d’eux. De plus, le polymorphisme des gènes du CMHII n’est pas limité qu’aux séquences codantes. Il est également observé dans les promoteurs où on a démontré ses effets sur le niveau d’expression des gènes. La variation de la régulation d’un gène est considérée comme un facteur important et pour laquelle des modifications morphologiques, physiologiques et comportementales sont observées chez tous les organismes. Des séquences d’ADN répétées impliquées dans cette régulation ont été identifiées dans les régions non-codantes des génomes. D’un autre côté, la sélection par les pathogènes permettrait l’évolution et le maintien du polymorphisme des gènes du CMH chez les vertébrés. À ce sujet, plusieurs études ont montré l’implication de différents allèles du CMH dans la résistance ou la susceptibilité aux maladies. Cette étude avait pour objectifs de caractériser le polymorphisme du gène MHIIb chez l’omble de fontaine (Salvelinus fontinalis) et de documenter ses effets au niveau de la survie conférée par des allèles et/ou génotypes particuliers lors d’une infection, ainsi que sur la variation du niveau d’expression du gène dans différentes conditions. Dans une première partie, nous avons identifié un total de 6 allèles du gène MHIIb, désignés Safo-DAB*0101 à Safo-DAB*0601, qui montrent une grande similarité avec les séquences codantes provenant de poissons téléostéens et de l’humain. L’analyse des séquences du domaine b1 a permis de détecter l’effet d’une pression sélective positive pour maintenir le polymorphisme dans cette région de la molécule. Quatre de ces allèles ont été testés lors d’une expérience d’infection avec le pathogène Aeromonas salmonicida afin d’évaluer l’effet qu’ils pouvaient avoir sur la survie des poissons. Nous avons trouvé que l’allèle DAB*0101 était significativement associé à la résistance à la furonculose. En plus d’avoir été identifié chez les individus homozygotes pour cet allèle, l’effet a également été remarqué au niveau de la survie les poissons de génotype DAB*0101/*0201. À l’opposé, les facteurs de risque élevé obtenus pour les génotypes DAB*0201/*0301 et DAB*0301/*0401 suggèrent plutôt une association à la susceptibilité. Étant donné la faible fréquence à laquelle l’allèle DAB*0101 a été retrouvé dans la population, le modèle de la sélection dépendante de la fréquence pourrait expliquer l’avantage conféré par ce dernier et souligne l’importance de ce mécanisme pour le maintien du polymorphisme du gène MHIIb chez l’omble de fontaine. Dans une seconde partie, nous avons rapporté la présence d’un minisatellite polymorphique formé d’un motif de 32 nucléotides dans le second intron du gène MHIIb, et pour lequel un nombre exclusif de répétitions du motif a été associé à chaque allèle (69, 27, 20, 40, 19 et 25 répétitions pour les allèles DAB*0101 à DAB*0601 respectivement). L’expression relative de quatre allèles a été évaluée dans des poissons hétérozygotes aux températures de 6 ºC et 18 ºC. Les résultats indiquent que les allèles possédant un long minisatellite montrent une réduction de l’expression du gène d’un facteur 1,67 à 2,56 par rapport aux allèles qui en contiennent un court. De même, des allèles qui incluent des minisatellites de tailles similaires n’affichent pas de différence significative au niveau de l’abondance du transcrit aux deux températures. De plus, l’effet répressif associé aux longs minisatellites est amplifié à la température de 18 ºC dans des poissons de trois génotypes différents. Nous avons finalement observé une augmentation significative par un facteur 2,08 de l’expression totale du gène MHIIb à la température de 6 ºC. Ces résultats appuient l’implication des séquences d’ADN répétées dans la régulation de l’activité transcriptionnelle d’un gène et suggèrent qu’un minisatellite sensible aux différences de températures pourrait être soumis aux forces sélectives et jouer un rôle important dans l’expression de gènes et l’évolution des organismes poïkilothermes.

Développement de procédés efficaces pour la construction et la production de vecteurs adénoviraux

L’adénovirus possède plusieurs caractéristiques faisant de ce virus un candidat de choix pour la construction de vecteurs utiles dans les études de génomique fonctionnelle. Dans la majorité de ces applications, on a recours à un vecteur adénoviral de première génération délété de sa région E1. L’utilisation de vecteurs adénoviraux comprend deux maillons faibles : la construction du vecteur et la production subséquente de ce dernier. Le développement de méthodes alternatives est donc nécessaire pour renforcer ces deux maillons, permettant ainsi une utilisation étendue de ces vecteurs. Ce développement va s’articuler sur deux axes : l’ingénierie du vecteur de transfert pour la construction de l’adénovirus recombinant et l’ingénierie d’une lignée cellulaire pour la production du vecteur. En utilisant un vecteur de transfert adénoviral co-exprimant, à partir d’un promoteur régulable à la tétracycline, la protéase de l’adénovirus et une protéine de fluorescence verte (GFP) par l’intermédiaire d’un site d’entrée ribosomal interne (IRES), notre groupe a établi que la sélection positive, via l’expression ectopique de la protéase, est un processus efficace pour la création de librairie d’adénovirus recombinants. Par contre, la diversité atteinte dans ce premier système est relativement faible, environ 1 adénovirus recombinant par 1 000 cellules. Le travail effectué dans le cadre de cette thèse vise à construire un nouveau transfert de vecteur dans lequel l’expression de la protéase sera indépendante de celle du transgène permettant ainsi d’optimiser l’expression de la protéase. Ce travail d’optimisation a permis de réduire le phénomène de transcomplémentation du virus parental ce qui a fait grimper la diversité à 1 virus recombinant par 75 cellules. Ce système a été mis à l’épreuve en générerant une librairie adénovirale antisens dirigée contre la GFP. La diversité de cette librairie a été suffisante pour sélectionner un antisens réduisant de 75% l’expression de la GFP. L’amplification de ce vecteur adénoviral de première génération doit se faire dans une lignée cellulaire exprimant la région E1 telle que les cellules 293. Par contre, un adénovirus de première génération se répliquant dans les cellules 293 peut échanger, par recombinaison homologue, son transgène avec la région E1 de la cellule créant ainsi un adénovirus recombinant réplicatif (RCA), compromettant ainsi la pureté des stocks. Notre groupe a déjà breveté une lignée cellulaire A549 (BMAdE1) exprimant la région E1, mais qui ne peut pas recombiner avec le transgène du virus. Par contre, le niveau de réplication de l’adénovirus dans les BMAdE1 est sous-optimal, à peine 15-30% du niveau obtenu dans les cellules 293. Le travail fait dans le cadre de cette thèse a permis de mettre en évidence qu’une expression insuffisante d’E1B-55K était responsable de la mauvaise réplication du virus dans les BMAdE1. Nous avons produit de nouveaux clones à partir de la lignée parentale via une transduction avec un vecteur lentiviral exprimant E1B-55K. Nous avons confirmé que certains clones exprimaient une plus grande quantité d’E1B-55K et que ces clones amplifiaient de manière plus efficace un vecteur adénoviral de première génération. Ce clone a par la suite été adapté à la culture en suspension sans sérum.

Estimation des corrélations phylogénétiques entre paramètres d'évolution moléculaire et Traits d'histoire de vie

Depuis quelques années, l'évolution moléculaire cherche à caractériser les variations et l'intensité de la sélection grâce au rapport entre taux de substitution synonyme et taux de substitution non-synonyme (dN/dS). Cette mesure, dN/dS, a permis d'étudier l'histoire de la variation de l'intensité de la sélection au cours du temps ou de détecter des épisodes de la sélection positive. Les liens entre sélection et variation de taille efficace interfèrent cependant dans ces mesures. Les méthodes comparatives, quant a elle, permettent de mesurer les corrélations entre caractères quantitatifs le long d'une phylogénie. Elles sont également utilisées pour tester des hypothèses sur l'évolution corrélée des traits d'histoire de vie, mais pour être employées pour étudier les corrélations entre traits d'histoire de vie, masse, taux de substitution ou dN/dS. Nous proposons ici une approche combinant une méthode comparative basée sur le principe des contrastes indépendants et un modèle d'évolution moléculaire, dans un cadre probabiliste Bayésien. Intégrant, le long d'une phylogénie, sur les reconstructions ancestrales des traits et et de dN/dS nous estimons les covariances entre traits ainsi qu'entre traits et paramètres du modèle d'évolution moléculaire. Un modèle hiérarchique, a été implémenté dans le cadre du logiciel coevol, publié au cours de cette maitrise. Ce modèle permet l'analyse simultané de plusieurs gènes sans perdre la puissance donnée par l'ensemble de séquences. Un travail deparallélisation des calculs donne la liberté d'augmenter la taille du modèle jusqu'à l'échelle du génome. Nous étudions ici les placentaires, pour lesquels beaucoup de génomes complets et de mesures phénotypiques sont disponibles. À la lumière des théories sur les traits d'histoire de vie, notre méthode devrait permettre de caractériser l'implication de groupes de gènes dans les processus biologique liés aux phénotypes étudiés.

Le rôle de la MAPK non-conventionnelle ERK3 dans le développement thymique

Les premières cellules progénitrices lympho-myéloïdes (LMPP) entrent dans le thymus et commencent leur processus de différenciation au stade double négatif (DN, CD4-CD8-). Après un réarrangement fonctionnel de la chaine bêta de leur récepteur des cellules T (RCT), les cellules reçoivent des signaux de différenciation, de prolifération, de survie et de sélection et passent au stade double positif (DP, CD4+CD8+). Ensuite, la chaine alpha du RCT est réarrangée et testée via les sélections positive et négative. Si le RCT reçoit un signal ni trop fort, ni trop faible, les cellules passent au stade simple positif où elles exprimeront soit la molécule CD4 ou CD8. ERK3, une MAPK non-conventionnelle, joue un rôle important dans le développement thymique. Des études précédentes ont démontré qu’une déficience en ERK3 diminue de 50 % la cellularité thymique et de 75% le nombre de cellules simples positives CD4 (CD4SP). Nous avons posé comme hypothèses qu’il y a une augmentation de l’apoptose chez les thymocytes DP de souris Erk3-/- et que cette déficience chez les thymocytes DP affecterait la sélection positive des cellules DP, réduisant ainsi le nombre de thymocytes CD4SP. Afin de vérifier la première hypothèse, nous avons regardé les niveaux d’apoptose grâce à la cytométrie en flux et par immunohistochimie. Dans les deux cas, nous étions incapables de détecter une différence du niveau d’apoptose chez les thymocytes DP entre les souris Erk3+/+ et Erk3-/-. Ensuite, nous nous sommes posés la question suivante : La demi-vie des thymocytes DP Erk3-/- était-elle plus courte que le type sauvage? La demi-vie des thymocytes DP a été mesurée à l’aide de l’étude des réarrangements secondaires de la chaine alpha du RCT par PCR semi-quantitatif et à l’aide de cultures de thymus fœtaux. En effet, ERK3 semble important pour prolonger la demi-vie des thymocytes DP. Ensuite, nous avons utilisé des marqueurs cellulaires différentiels (CD69, CD5 et RCT) pour regarder si les thymocytes DP sont capables de passer la sélection positive. En effet, il y a moins de thymocytes DP Erk3-/- qui sont CD69fort, CD5fort et RCTfort. Finalement, nous voulons savoir si les fonctions de ERK3 passent par MK5, son seul partenaire d’interaction connu à ce jour. Après la caractérisation du thymus de la souris Mk5-/-, nous observons que seulement la réduction du nombre de thymocytes CD4SP est identique à celle des thymocytes CD4SP de la souris Erk3-/-. En conclusion, ces résultats révèlent des fonctions importantes pour la molécule ERK3 lors du processus de sélection positive, le maintient de la demi-vie des thymocytes DP et lors de la régulation de développement thymique de manière MK5-dépendante et -indépendante.

Rôle de la tolérance centrale et périphérique des lymphocytes T autoréactifs dans deux nouveaux modèles murins double transgéniques

Les maladies autoimmunes sont des affections chroniques, le plus souvent invalidantes, qui touchent plus de 5% de la population dans les pays développés. L’autoimmunité résulte de la rupture des mécanismes de tolérance du système immunitaire vis-à-vis des autoantigènes exprimés par les tissus de l’organisme, entraînant la destruction d’un ou de plusieurs organes-cibles par les lymphocytes T et/ou B. L’hépatite autoimmune et le diabète autoimmun se caractérisent par la destruction sélective des hépatocytes et des cellules beta pancréatiques, respectivement. De plus en plus d’arguments suggèrent une implication des lymphocytes T CD8+ dans le déclenchement, la progression et la régulation des réponses associées à plusieurs maladies autoimmunes. Dans ce projet, nous avons suivi l’évolution de clones de lymphocytes T CD8+ spécifiques à un antigène particulier dont le site d’expression différait. Pour ce faire, nous avons développé deux nouveaux modèles murins double transgéniques par croisement entre une lignée de souris exprimant un TCR transgénique spécifique à la nucléoprotéine (NP) du virus de la chorioméningite lymphocytaire (LCMV), et une souris exprimant cette NP-LCMV : 1) uniquement dans les hépatocytes (modèle d’hépatite autoimmune), ou 2) simultanément dans le thymus et le pancréas (modèle de diabète autoimmun). L’avidité fonctionnelle des lymphocytes T CD8+ spécifiques à la NP chez les souris TCR transgéniques était inversement proportionnelle au niveau d’expression du TCR. Le répertoire lymphocytaire dans le thymus, la rate, les ganglions et le sang périphérique a été caractérisé pour chacune des lignées de souris double transgéniques, de même que la capacité fonctionnelle et le phénotype (marqueurs d’activation/mémoire) des lymphocytes T CD8+ autoréactifs. Chacun des deux nouveaux modèles présentés dans cette étude ont montré que les lymphocytes T CD8+ spécifiques à la NP sont aptes à briser la tolérance centrale et périphérique et à provoquer une réaction d’autoimmunité spontanée. Dans le modèle d’hépatite autoimmune, où l’expression de l’autoantigène était restreinte au foie, la surexpression du TCR transgénique a entraîné une délétion thymique quasi-totale des lymphocytes T CD8+ spécifiques à la NP prévenant le développement d’une hépatite spontanée. alors qu’un niveau de TCR comparable à celui d’une souris de type sauvage a permis une sélection positive des lymphocytes autoréactifs qui se sont accumulés dans le foie où ils se sont activés pour provoquer une hépatite autoimmune spontanée. Dans le modèle de diabète autoimmun, où l’autoantigène était exprimé dans le pancréas et le thymus, les souris des deux lignées double transgéniques ont montré une délétion thymique partielle, peu importe le niveau d’expression du TCR. Seuls les mâles adultes développaient un diabète spontané et une partie de leurs lymphocytes T CD8+ exprimaient une combinaison particulière de marqueurs d’activation/mémoire (CD44, CD122, PD-1). Cette population lymphocytaire était absente chez les souris femelles et les mâles sains. L’étude de la tolérance des lymphocytes T CD8+ autoréactifs dans nos deux nouveaux modèles murins double transgéniques a permis d’identifier des mécanismes alternatifs possiblement impliqués dans la tolérance et l’activation, et de mieux comprendre le rôle des lymphocytes T CD8+ autoréactifs dans le processus autoimmun menant à l’hépatite autoimmune et au diabète autoimmun. Ces découvertes seront utiles pour développer de nouvelles approches thérapeutiques ciblant les lymphocytes T CD8+ autoréactifs.