MODELING HETEROTACHY IN PHYLOGENETICS

Il a été démontré que l’hétérotachie, variation du taux de substitutions au cours du temps et entre les sites, est un phénomène fréquent au sein de données réelles. Échouer à modéliser l’hétérotachie peut potentiellement causer des artéfacts phylogénétiques. Actuellement, plusieurs modèles traitent l’hétérotachie : le modèle à mélange des longueurs de branche (MLB) ainsi que diverses formes du modèle covarion. Dans ce projet, notre but est de trouver un modèle qui prenne efficacement en compte les signaux hétérotaches présents dans les données, et ainsi améliorer l’inférence phylogénétique. Pour parvenir à nos fins, deux études ont été réalisées. Dans la première, nous comparons le modèle MLB avec le modèle covarion et le modèle homogène grâce aux test AIC et BIC, ainsi que par validation croisée. A partir de nos résultats, nous pouvons conclure que le modèle MLB n’est pas nécessaire pour les sites dont les longueurs de branche diffèrent sur l’ensemble de l’arbre, car, dans les données réelles, le signaux hétérotaches qui interfèrent avec l’inférence phylogénétique sont généralement concentrés dans une zone limitée de l’arbre. Dans la seconde étude, nous relaxons l’hypothèse que le modèle covarion est homogène entre les sites, et développons un modèle à mélanges basé sur un processus de Dirichlet. Afin d’évaluer différents modèles hétérogènes, nous définissons plusieurs tests de non-conformité par échantillonnage postérieur prédictif pour étudier divers aspects de l’évolution moléculaire à partir de cartographies stochastiques. Ces tests montrent que le modèle à mélanges covarion utilisé avec une loi gamma est capable de refléter adéquatement les variations de substitutions tant à l’intérieur d’un site qu’entre les sites. Notre recherche permet de décrire de façon détaillée l’hétérotachie dans des données réelles et donne des pistes à suivre pour de futurs modèles hétérotaches. Les tests de non conformité par échantillonnage postérieur prédictif fournissent des outils de diagnostic pour évaluer les modèles en détails. De plus, nos deux études révèlent la non spécificité des modèles hétérogènes et, en conséquence, la présence d’interactions entre différents modèles hétérogènes. Nos études suggèrent fortement que les données contiennent différents caractères hétérogènes qui devraient être pris en compte simultanément dans les analyses phylogénétiques.

Surveillance de maladies chroniques à l'aide des données administratives : cas de l'asthme au Québec

Thèse numérisée par la Division de la gestion de documents et des archives de l'Université de Montréal.

Inférence robuste à la présence des valeurs aberrantes dans les enquêtes

Cette thèse comporte trois articles dont un est publié et deux en préparation. Le sujet central de la thèse porte sur le traitement des valeurs aberrantes représentatives dans deux aspects importants des enquêtes que sont : l’estimation des petits domaines et l’imputation en présence de non-réponse partielle. En ce qui concerne les petits domaines, les estimateurs robustes dans le cadre des modèles au niveau des unités ont été étudiés. Sinha & Rao (2009) proposent une version robuste du meilleur prédicteur linéaire sans biais empirique pour la moyenne des petits domaines. Leur estimateur robuste est de type «plugin», et à la lumière des travaux de Chambers (1986), cet estimateur peut être biaisé dans certaines situations. Chambers et al. (2014) proposent un estimateur corrigé du biais. En outre, un estimateur de l’erreur quadratique moyenne a été associé à ces estimateurs ponctuels. Sinha & Rao (2009) proposent une procédure bootstrap paramétrique pour estimer l’erreur quadratique moyenne. Des méthodes analytiques sont proposées dans Chambers et al. (2014). Cependant, leur validité théorique n’a pas été établie et leurs performances empiriques ne sont pas pleinement satisfaisantes. Ici, nous examinons deux nouvelles approches pour obtenir une version robuste du meilleur prédicteur linéaire sans biais empirique : la première est fondée sur les travaux de Chambers (1986), et la deuxième est basée sur le concept de biais conditionnel comme mesure de l’influence d’une unité de la population. Ces deux classes d’estimateurs robustes des petits domaines incluent également un terme de correction pour le biais. Cependant, ils utilisent tous les deux l’information disponible dans tous les domaines contrairement à celui de Chambers et al. (2014) qui utilise uniquement l’information disponible dans le domaine d’intérêt. Dans certaines situations, un biais non négligeable est possible pour l’estimateur de Sinha & Rao (2009), alors que les estimateurs proposés exhibent un faible biais pour un choix approprié de la fonction d’influence et de la constante de robustesse. Les simulations Monte Carlo sont effectuées, et les comparaisons sont faites entre les estimateurs proposés et ceux de Sinha & Rao (2009) et de Chambers et al. (2014). Les résultats montrent que les estimateurs de Sinha & Rao (2009) et de Chambers et al. (2014) peuvent avoir un biais important, alors que les estimateurs proposés ont une meilleure performance en termes de biais et d’erreur quadratique moyenne. En outre, nous proposons une nouvelle procédure bootstrap pour l’estimation de l’erreur quadratique moyenne des estimateurs robustes des petits domaines. Contrairement aux procédures existantes, nous montrons formellement la validité asymptotique de la méthode bootstrap proposée. Par ailleurs, la méthode proposée est semi-paramétrique, c’est-à-dire, elle n’est pas assujettie à une hypothèse sur les distributions des erreurs ou des effets aléatoires. Ainsi, elle est particulièrement attrayante et plus largement applicable. Nous examinons les performances de notre procédure bootstrap avec les simulations Monte Carlo. Les résultats montrent que notre procédure performe bien et surtout performe mieux que tous les compétiteurs étudiés. Une application de la méthode proposée est illustrée en analysant les données réelles contenant des valeurs aberrantes de Battese, Harter & Fuller (1988). S’agissant de l’imputation en présence de non-réponse partielle, certaines formes d’imputation simple ont été étudiées. L’imputation par la régression déterministe entre les classes, qui inclut l’imputation par le ratio et l’imputation par la moyenne sont souvent utilisées dans les enquêtes. Ces méthodes d’imputation peuvent conduire à des estimateurs imputés biaisés si le modèle d’imputation ou le modèle de non-réponse n’est pas correctement spécifié. Des estimateurs doublement robustes ont été développés dans les années récentes. Ces estimateurs sont sans biais si l’un au moins des modèles d’imputation ou de non-réponse est bien spécifié. Cependant, en présence des valeurs aberrantes, les estimateurs imputés doublement robustes peuvent être très instables. En utilisant le concept de biais conditionnel, nous proposons une version robuste aux valeurs aberrantes de l’estimateur doublement robuste. Les résultats des études par simulations montrent que l’estimateur proposé performe bien pour un choix approprié de la constante de robustesse.

Une nouvelle méthode de détection du glaucome par la mesure de l'asymétrie interoculaire : l’asymétrie du rapport de la surface neurorétinienne sur la surface du disque optique ou rim area to disc area asymmetry ratio (RADAAR)

Cette thèse constitue à la fois un apport de nature clinique et technologique dans l’approche diagnostique du glaucome. Plus précisément, nous nous proposons d’étudier une nouvelle façon de détecter le glaucome par la mesure de l’asymétrie du rapport de la surface de l’anneau neurorétinien et de la surface de la papille ou du disque optique ou rim to disc area asymmetry ratio (RADAAR). Pour atteindre cet objectif, nous avons recours à une base de données composée d’une population subdivisée en 4 différents groupes de diagnostic (normal, glaucome possible, glaucome probable et glaucome définitif). Les mesures du RADAAR sont calculées de différentes façons à partir des paramètres stéréométriques de la tête du nerf optique des sujets, produits par la microscopie confocale à balayage laser (Heidelberg Retina Tomograph (HRT) (Heidelberg Engineering, Germany)). Nous procédons à une analyse de données grâce au logiciel SPSS où nous mettons en exergue la distribution du RADAAR dans les différentes populations, sa validité et son utilité dans le dépistage du glaucome. Nous enrôlons donc 523 sujets dans cette étude avec 82 sujets atteints de glaucome définitif. La moyenne d’âge est de 62 ans. Il y a plus de femmes que d’hommes et plus de Caucasiens que d’Africains Caribéens. Nous trouvons que la distribution de la mesure du RADAAR est différente des sujets d’un groupe de diagnostic à l’autre. En termes de performance, la sensibilité de la mesure du RADAAR est très basse c'est-à-dire que sa capacité de détecter la maladie est basse. En revanche la mesure du RADAAR est plus spécifique c'est-à-dire que sa capacité d’identifier les sujets exempts de la maladie est plus grande. Elle tendrait à être aussi plus performante chez les Africains Caribéens que chez les Caucasiens. De même, elle serait plus sensible chez les hommes que chez les femmes. La mesure du RADAAR est utile si on l’associe à une autre méthode de diagnostic comme l’analyse de Régression de Moorfields (MRA) incluse dans le logiciel du HRT3 spécialement lors de la détection du glaucome dans la population à haut risque. En définitive, nous déterminons que la mesure du RADAAR se veut un outil d’aide au diagnostic. Elle est particulièrement intéressante dans le contexte de dépistage de glaucome.

Impact des polluants agricoles sur la génétique des populations d’une espèce sentinelle : le ouaouaron (Rana catesbeiana)

L’expansion agricole ne cesse d’agir sur la perte d’habitats essentiels et nécessaires au développement des espèces. Bien que plusieurs espèces réussissent à survivre dans ces habitats peu adéquats, la persistance et la santé de plusieurs populations semblent compromises par l’utilisation souvent intensive de polluants chimiques agricoles et de fertilisants. Cette étude a pour but de déterminer l’impact des contaminants et de l’écologie du paysage sur la diversité génétique des populations de ouaouarons retrouvées en milieu agricole. Notre hypothèse de départ stipule qu’une exposition chronique aux polluants agricoles induira des différences génétiques au niveau des populations exposées. Le bassin versant de la rivière Yamaska a été désigné comme site d’étude puisqu’il fait partie de la région agricole la plus importante du Québec et parce qu’on y retrouve un gradient d’utilisation des terres pour l’agriculture (faible, moyen, élevé). Le ouaouaron a été choisi à titre de modèle biologique puisque ses caractéristiques physiologiques et écologiques en font une espèce sentinelle capable de rendre compte de l’état de santé global des écosystèmes. La caractérisation génétique des populations a été effectuée à partir de marqueurs d’AFLP (Amplified Fragment Length Polymorphism). Les résultats montrent que la diversité génétique est liée à la colonisation à partir de l’embouchure de la rivière Yamaska et que quelques populations sont génétiquement différenciées. De plus, nous avons démontré une relation positive entre le nombre de locus polymorphes et l’atrazine, l’indice de contamination et le métolachlore et la concentration en azote ainsi qu’entre l’hétérozygotie attendue et la concentration en phosphate.

Le peuplement fondateur de la région de Lotbinière et ses conséquences démogénétiques

Un effet fondateur survient lorsqu’un petit nombre d’immigrants forment une nouvelle population et qu’ainsi les descendants ont une majorité de gènes provenant de ces quelques ancêtres. L’effet fondateur québécois, qui résulte de l’établissement de quelques milliers d’immigrants français aux XVIIe et XVIIIe siècles, est bien documenté. Mais des effets fondateurs régionaux ont aussi été identifiés. Ce mémoire de maîtrise vise à déterminer si un effet fondateur régional est à l’oeuvre dans la région de Lotbinière (Chaudière-Appalaches), dont le peuplement initial remonte à la fin du XVIIe siècle. Le fichier BALSAC et le Registre de la population du Québec ancien ont permis de constituer deux groupes de descendants, 715 individus mariés à la fin du XVIIIe siècle, et 60 autres mariés à la fin du XXe siècle. Par généalogies ascendantes et descendantes, les fondateurs immigrants et régionaux de la région ont par la suite été identifiés. Les résultats indiquent que l’effet fondateur régional avait encore une forte empreinte chez le groupe de descendants du XVIIIe siècle, mais que l’impact s’atténue en ce qui concerne les descendants contemporains. L’homogénéité démontrée par les coefficients d’apparentement et l’indice de contribution génétique uniforme, le petit nombre de fondateurs régionaux et le fait que 65 % des gènes contemporains étaient déjà introduits en 1800 sont des signes qui pointent vers un effet fondateur régional. Par contre, le nonisolement de la région, la proportion modérée de gènes contemporains introduits par les premiers fondateurs régionaux et les niveaux de consanguinité semblables aux autres régions du centre du Québec, incitent à nuancer cette conclusion. En fait, il y a possiblement deux Lotbinière : le Lotbinière ancien, sur la rive et le Lotbinière nouveau, dans les terres; chacun ayant son pool génique et son historique de peuplement propre.

Réfugiés climatiques : statut et traitement

L’état actuel des travaux ne rend pas compte de l’ampleur des questions philosophiques et des enjeux moraux suscités par l’apparition sur la scène internationale des réfugiés climatiques. Il est pourtant urgent de leur accorder une protection à travers des accords internationaux. Les philosophes qui se sont penchés sur le sujet ont été induits en erreur tant par la multiplicité des termes employés que leur absence de définitions. Ce travail critique la tendance actuelle des militants écologistes à vouloir englober des populations aux problèmes divers sous le terme de réfugié. Banaliser l’emploi du terme de réfugié n’est pas seulement fallacieux mais également dangereux. A terme, les militants se tourneront vers la Convention de Genève pour revendiquer que les populations déplacées soient considérées comme des réfugiés. Or la Convention de Genève n’est pas un outil adéquat pour remédier au sort de ces populations. De plus, on ne peut élargir le statut de réfugié pour inclure ces populations sans risquer de perdre sa crédibilité et son efficience. Suivre la pente qu’emprunte les militants nous mènerait à accorder le même traitement aux réfugiés climatiques et aux réfugiés politiques, ce qui est une erreur. Notre hypothèse est que les habitants des petits pays insulaires à l’inverse des autres populations ont besoin d’un élargissement de la Convention de Genève. Nous arguerons que nous avons des devoirs et des responsabilités envers eux que nous n’avons pas envers les réfugiés politiques. Pour défendre ce point de vue, il faut définir clairement ce qu’est un réfugié climatique et justifier cette appellation. Nous devrons donc confronter la notion de réfugié climatique à d’autres notions concurrentes. Une fois les termes définis, nous envisagerons les enjeux éthiques à travers le prisme des questions de justice globale. Nous verrons que pour déterminer qui devrait remédier au sort des réfugiés climatique, il ne suffit pas de se référer à la responsabilité causale. Cela nous mènera à arguer que bien que séduisant, le principe pollueur-payeur n’est pas un outil adéquat pour guider la réflexion. Nous serons également amenés à nous interroger sur la pertinence d’une institution environnementale globale.

Clustering algorithms and shape factor methods to discriminate among small GTPase phenotypes using DIC image analysis.

Naïvement perçu, le processus d’évolution est une succession d’événements de duplication et de mutations graduelles dans le génome qui mènent à des changements dans les fonctions et les interactions du protéome. La famille des hydrolases de guanosine triphosphate (GTPases) similaire à Ras constitue un bon modèle de travail afin de comprendre ce phénomène fondamental, car cette famille de protéines contient un nombre limité d’éléments qui diffèrent en fonctionnalité et en interactions. Globalement, nous désirons comprendre comment les mutations singulières au niveau des GTPases affectent la morphologie des cellules ainsi que leur degré d’impact sur les populations asynchrones. Mon travail de maîtrise vise à classifier de manière significative différents phénotypes de la levure Saccaromyces cerevisiae via l’analyse de plusieurs critères morphologiques de souches exprimant des GTPases mutées et natives. Notre approche à base de microscopie et d’analyses bioinformatique des images DIC (microscopie d’interférence différentielle de contraste) permet de distinguer les phénotypes propres aux cellules natives et aux mutants. L’emploi de cette méthode a permis une détection automatisée et une caractérisation des phénotypes mutants associés à la sur-expression de GTPases constitutivement actives. Les mutants de GTPases constitutivement actifs Cdc42 Q61L, Rho5 Q91H, Ras1 Q68L et Rsr1 G12V ont été analysés avec succès. En effet, l’implémentation de différents algorithmes de partitionnement, permet d’analyser des données qui combinent les mesures morphologiques de population native et mutantes. Nos résultats démontrent que l’algorithme Fuzzy C-Means performe un partitionnement efficace des cellules natives ou mutantes, où les différents types de cellules sont classifiés en fonction de plusieurs facteurs de formes cellulaires obtenus à partir des images DIC. Cette analyse démontre que les mutations Cdc42 Q61L, Rho5 Q91H, Ras1 Q68L et Rsr1 G12V induisent respectivement des phénotypes amorphe, allongé, rond et large qui sont représentés par des vecteurs de facteurs de forme distincts. Ces distinctions sont observées avec différentes proportions (morphologie mutante / morphologie native) dans les populations de mutants. Le développement de nouvelles méthodes automatisées d’analyse morphologique des cellules natives et mutantes s’avère extrêmement utile pour l’étude de la famille des GTPases ainsi que des résidus spécifiques qui dictent leurs fonctions et réseau d’interaction. Nous pouvons maintenant envisager de produire des mutants de GTPases qui inversent leur fonction en ciblant des résidus divergents. La substitution fonctionnelle est ensuite détectée au niveau morphologique grâce à notre nouvelle stratégie quantitative. Ce type d’analyse peut également être transposé à d’autres familles de protéines et contribuer de manière significative au domaine de la biologie évolutive.

Living and fishing in a marine protected area : balancing traditional fisheries with conservation in Karimunjawa National Park, Indonesia

Cette recherche a porté sur quelques enjeux importants liés à la gestion des aires marines protégées (AMP) en Indonésie en examinant comment celles-ci sont en mesure d'adapter leurs politiques afin de mieux répondre à l'évolution des conditions socioéconomiques et écologiques, quels ont été les impacts socioéconomiques de ces aires, et quelles sont les préoccupations environnementales des acteurs locaux dont les moyens de subsistance dépendent des ressources règlementées. Le «livelihoods framework » a servi de guide pour notre analyse des changements socioéconomiques dans la région, tandis que la notion d’« environmentality » d’Agrawal a fourni les bases théoriques pour l'examen de la formation de sujets environnementaux au parc national de Karimunjawa. Cette étude a montré que les changements de politique apportés au plan de la gestion du parc sont un pas dans la bonne direction, mais que les objectifs importants liés sa cogestion n'ont jamais été entièrement réalisés dans la pratique. Les résultats montrent également que d'importants changements socioéconomiques surviennent dans le parc, de nombreux pêcheurs se tournent vers des moyens de subsistance alternatifs, afin de compenser la baisse des prises de poissons. Enfin, cette étude a révélé que d'importants changements positifs dans les préoccupations environnementales sont survenus depuis la modification du zonage du parc, mais que ceux-ci ne se sont pas entièrement traduits en conformité avec les règles et règlements de l'AMP.

Immigration et alimentation à Montréal aux XVIIe et XVIIIe siècles : essai d’interprétation à partir d’analyses isotopiques sur des populations archéologiques

Afin d’étudier l’influence de la migration sur l’alimentation à Montréal aux XVIIe et XVIIIe siècles, 64 individus de la collection du cimetière Notre-Dame, daté de 1691 à 1796, ont fait l’objet d’analyses ostéologiques et isotopiques. Les analyses isotopiques ont portées sur le carbone (d13C) et l’azote (d15N) du collagène des os, ainsi que sur le d13C et l’oxygène (d18O) du carbonate de l’apatite des os et des dents (prémolaires et troisièmes molaires). Le d18O des dents a permis de définir approximativement trois régions d’origine (région de Montréal, région enrichie en 18O (i.e. Acadie, Louisiane, Nouvelle-Angleterre, France, Antilles et Afrique) et région appauvrie en 18O (intérieur des terres et plus au nord) pour 58 individus, et sept possibles parcours migratoire (N=27). Plus de la moitié de l’échantillon est composé d’individus possiblement natifs de Montréal (55 %). De plus, les résultats indiquent que les gens étaient peu mobiles avant l’âge de 16 ans. Toutefois, 12 individus ont entrepris des déplacements entre 7 et 16 ans, majoritairement d’un environnement enrichi vers Montréal (N=5) ou de Montréal vers une région appauvrie (N=5). L’âge de recrutement des mousses sur les navires, la traite de la fourrure, la coupe du bois et possiblement aussi l’esclavage pourraient expliquer cette « jeune » migration. Sur le plan alimentaire, les végétaux de type C3, la viande nourrie aux ressources C3 et le poisson faisaient partie du menu montréalais. Les plantes C4 (majoritairement maïs mais aussi sucre de canne [rhum]) étaient consommées en quantité variable. La question de l’influence de la migration sur l’alimentation n’a pu être explorée en profondeur en raison de contraintes liées à la contamination du d18O du carbonate des os. La combinaison des données ostéologiques et isotopiques à la distribution spatiale des sépultures, a permis d’étudier un aspect de l’archéologie funéraire à l’échelle individuelle (identité possible), sans toutefois fournir de résultats probants, à l’échelle du cimetière et de son organisation globale.

Relevé spectroscopique et étude des propriétés physiques des étoiles naines blanches à moins de 40 parsecs du Soleil

Les étoiles naines blanches représentent la fin de l’évolution de 97% des étoiles de notre galaxie, dont notre Soleil. L’étude des propriétés globales de ces étoiles (distribution en température, distribution de masse, fonction de luminosité, etc.) requiert l’élaboration d’ensembles statistiquement complets et bien définis. Bien que plusieurs relevés d’étoiles naines blanches existent dans la littérature, la plupart de ceux-ci souffrent de biais statistiques importants pour ce genre d’analyse. L’échantillon le plus représentatif de la population d’étoiles naines blanches demeure à ce jour celui défini dans un volume complet, restreint à l’environnement immédiat du Soleil, soit à une distance de 20 pc (∼ 65 années-lumière) de celui-ci. Malheureusement, comme les naines blanches sont des étoiles intrinsèquement peu lumineuses, cet échantillon ne contient que ∼ 130 objets, compromettant ainsi toute étude statistique significative. Le but de notre étude est de recenser la population d’étoiles naines blanches dans le voisinage solaire a une distance de 40 pc, soit un volume huit fois plus grand. Nous avons ainsi entrepris de répertorier toutes les étoiles naines blanches à moins de 40 pc du Soleil à partir de SUPERBLINK, un vaste catalogue contenant le mouvement propre et les données photométriques de plus de 2 millions d’étoiles. Notre approche est basée sur la méthode des mouvements propres réduits qui permet d’isoler les étoiles naines blanches des autres populations stellaires. Les distances de toutes les candidates naines blanches sont estimées à l’aide de relations couleur-magnitude théoriques afin d’identifier les objets se situant à moins de 40 pc du Soleil, dans l’hémisphère nord. La confirmation spectroscopique du statut de naine blanche de nos ∼ 1100 candidates a ensuite requis 15 missions d’observations astronomiques sur trois grands télescopes à Kitt Peak en Arizona, ainsi qu’une soixantaine d’heures allouées sur les télescopes de 8 m des observatoires Gemini Nord et Sud. Nous avons ainsi découvert 322 nouvelles étoiles naines blanches de plusieurs types spectraux différents, dont 173 sont à moins de 40 pc, soit une augmentation de 40% du nombre de naines blanches connues à l’intérieur de ce volume. Parmi ces nouvelles naines blanches, 4 se trouvent probablement à moins de 20 pc du Soleil. De plus, nous démontrons que notre technique est très efficace pour identifier les étoiles naines blanches dans la région peuplée du plan de la Galaxie. Nous présentons ensuite une analyse spectroscopique et photométrique détaillée de notre échantillon à l’aide de modèles d’atmosphère afin de déterminer les propriétés physiques de ces étoiles, notamment la température, la gravité de surface et la composition chimique. Notre analyse statistique de ces propriétés, basée sur un échantillon presque trois fois plus grand que celui à 20 pc, révèle que nous avons identifié avec succès les étoiles les plus massives, et donc les moins lumineuses, de cette population qui sont souvent absentes de la plupart des relevés publiés. Nous avons également identifié plusieurs naines blanches très froides, et donc potentiellement très vieilles, qui nous permettent de mieux définir le côté froid de la fonction de luminosité, et éventuellement l’âge du disque de la Galaxie. Finalement, nous avons aussi découvert plusieurs objets d’intérêt astrophysique, dont deux nouvelles étoiles naines blanches variables de type ZZ Ceti, plusieurs naines blanches magnétiques, ainsi que de nombreux systèmes binaires non résolus.

Variations dans la réponse de la diversité génétique de populations de couleuvres insulaires faisant face à la perte d’habitat

Projet de recherche réalisé avec Bernard Angers comme directeur de maîtrise, Denis Réale en tant que co-directeur et grâce à la collaboration active d'Emmanuel Milot.

Étude exploratoire sur les origines du peuplement de l’île de Madagascar : une approche de craniométrie comparative appliquée sur des populations modernes

Les origines du peuplement de l’île de Madagascar ne sont encore que partiellement explorées à l’heure actuelle. Différentes populations ont contribué au peuplement de l’île, de nombreuses théories sur les origines de ce peuplement ont émergé et varient grandement selon les sources consultées. Selon l’archéologie et l’anthropologie culturelle, l’arrivée des premiers peuples remonterait à deux millénaires avant notre ère et plusieurs strates de vagues migratoires venues d’Afrique et d’Asie se sont succédées. Pour une vision complète du peuplement de toute l’île, ce sont les études en linguistique et en génétique qui ont donné les meilleures pistes en s’orientant vers une origine à prédominance indonésienne plutôt qu’africaine. Il reste cependant à confronter ces données diverses à celles issues de l’approche phénotypique, qui est peu utilisée. Mon objectif est donc d’explorer cette hypothèse à partir des données craniométriques, et ainsi de tester les modèles de peuplement proposés grâce à d’autres approches. Cet échantillon malgache (N=207) a été subdivisé sur la base de diverses données (géographie, ethnies et affiliations linguistiques). Après des analyses intra-groupe et intergroupes, ce dernier a été comparé à d’autres données craniométriques personnelles et publiées (N=1184). Deux types d’approches statistiques (multivariées classiques et issues de l’approche de la génétique des populations ou RMET) ont été utilisées afin d’obtenir des paramètres diversifiés et complémentaires. Les résultats issus des deux approches tendent vers une origine mixte (Afrique et Asie), dont la prépondérance varie en fonction de la région et du sexe. En effet, les hommes malgaches ont une origine triple (sud de l’Asie du Sud-est, sud de l’Afrique et côtes sud-est africaines), alors que les femmes ont plutôt une origine double (Afrique et Asie) selon l’approche multivariée classique. D’après les analyses RMET, on note que les individus des régions du nord et de l’est de l’île se rapprochent des populations de Tanzanie et les Malgaches présentent des similarités avec les populations indiennes. De plus, on remarque que les Malgaches du groupe nord présentent par rapport aux autres groupes un degré d’hétérogénéité plus élevé (Fst). Ce phénomène est dû probablement à des apports de populations plus diverses dès le début du peuplement de l’île dans cette région. Cette étude, basée sur un petit échantillon, confirme néanmoins les thèses antérieures sur la diversité du peuplement malgache et de plus elle démontre que les composantes prédominantes (Afrique ou Asie) varient selon les régions et le sexe.