Ecological niche modeling, cytogenetics and phylogeography of the genera Geomalacus and Letourneuxia (Gastropoda, Pulmonata) from the Iberian-Moroccan region

First described more that 150 years ago, the systematics of the genera Geomalacus and Letourneuxia (Arionidae, Gastropoda, Pulmonata) is still challenging. The taxonomic classification of arionid species is based on extremely labile characters such as body size or color that depends both on diet and environment, as well as age. Moreover, there is little information on the genetic diversity and population structure of the Iberian slugs that could provide extra clues to disentangle their problematic classification. The present work uses different analytical tools such as habitat suitability (Ecological Niche Modeling - ENM), cytogenetic analysis and phylogeography to establish the geographical distribution and evolutionary history of these pulmonate slugs. The potential distribution of the four Geomalacus species was modeled using ENM, which allowed the identification of new locations for G. malagensis, including a first report in Portugal. Also, it was predicted a much wider distribution for G. malagensis and G. oliveirae than previously known. Classical cytogenetic analyses were assayed with reproductive and a novel use of somatic tissues (mouth and tentacles) returning the number of chromosomes for the four Geomalacus species and L. numidica (n = 31, 2n = 62) and the respective karyotypes. G. malagensis and L. numidica present similar chromosome morphologies and karyotypic formulae, being more similar to each other than the Geomalacus among themselves. We further reconstructed the phylogeny of the genera Geomalacus and Letourneuxia using partial sequences of the mitochondrial cytochrome oxidase subunit I (COI) and the nuclear ribosomal small subunit (18S rRNA), and applied an independent evolutionary rate method, the indicator vectors correlation, to evaluate the existence of cryptic diversity within species. The five nominal species of Geomalacus and Letourneuxia comprise 14 well-supported cryptic lineages. Letourneuxia numidica was retrieved as a sister group of G. malagensis. G. oliveirae is paraphyletic with respect to G. anguiformis. According to our dating estimates, the most recent common ancestor of Geomalacus dates back to the Middle Miocene (end of the Serravallian stage). The major lineage splitting events within Geomalacus occurred during the dry periods of the Zanclean stage (5.3-3.6 million years) and some lineages were confined to more humid mountain areas of the Iberian Peninsula, which lead to a highly geographically structured mitochondrial genetic diversity. The major findings of this are the following: (1) provides updated species distribution maps for the Iberian Geomalacus expanding the known geographic distribution of the concerned species, (2) unravels the cryptic diversity within the genera Geomalacus and Letourneuxia, (3) Geomalacus oliveirae is paraphyletic with G. anguiformis and (4) Letourneuxia numidica is sister group of G. malagensis.

How historical and present climat conditions affected the distribution of the mediterranean whter shrew?:a phylogeographical and ecological approach

Tese de doutoramento, Biologia (Biologia Evolutiva), Universidade de Lisboa, Faculdade de Ciências, 2014

The identification and characterization of inter- and intra-species genetic diversity derived from retrotransposons in humans

Retrotransposons, which used to be considered as “junk DNA”, have begun to reveal their immense value to genome evolution and human biology due to recent studies. They consist of at least ~45% of the human genome and are more or less the same in other mammalian genomes. Retrotransposon elements (REs) are known to affect the human genome through many different mechanisms, such as generating insertion mutations, genomic instability, and alteration in gene expression. Previous studies have suggested several RE subfamilies, such as Alu, L1, SVA and LTR, are currently active in the human genome, and they are an important source of genetic diversity between human and other primates, as well as among humans. Although several groups had used Retrotransposon Insertion Polymorphisms (RIPs) as markers in studying primate evolutionary history, no study specifically focused on identifying Human-Specific Retrotransposon Element (HS-RE) and their roles in human genome evolution. In this study, by computationally comparing the human genome to 4 primate genomes, we identified a total of 18,860 HS-REs, among which are 11,664 Alus, 4,887 L1s, 1,526 SVAs and 783 LTRs (222 full length entries), representing the largest and most comprehensive list of HS-REs generated to date. Together, these HS-REs contributed a total of 14.2Mb sequence increase from the inserted REs and Target Site Duplications (TSDs), 71.6Kb increase from transductions, and 268.2 Kb sequence deletion of from insertion-mediated deletion, leading to a net increase of ~14 Mb sequences to the human genome. Furthermore, we observed for the first time that Y chromosome might be a hot target for new retrotransposon insertions in general and particularly for LTRs. The data also allowed for the first time the survey of frequency of TE insertions inside other TEs in comparison with TE insertion into none-TE regions. In summary, our data suggest that retrotransposon elements have played a significant role in the evolution of Homo sapiens.

La phylogénie moléculaire du genre nord-américain Eurybia (Asteraceae : Astereae) et ses proches parents (Oreostemma, Herrickia, Triniteurybia)

Eurybia et ses proches parents Oreostemma, Herrickia et Triniteurybia sont appelés le grade des eurybioïdes. Comprenant 31 espèces vivaces, ce grade appartient au clade Nord-américain de la tribu des Astereae. Les analyses moléculaires antérieures ont montré que ce groupe est à la fois paraphylétique aux Machaerantherinae et un groupe frère aux Symphyotrichinae. Les relations infragénériques partiellement résolues et faiblement supportées empêchent d’approfondir l'histoire évolutive des groupes et ce, particulièrement dans le genre principal Eurybia. Le but de cette étude est de reconstruire les relations phylogénétiques au sein des eurybioïdes autant par l'inclusion de toutes les espèces du grade que par l’utilisation de différents types de régions et de méthodes d'inférence phylogénétique. Cette étude présente des phylogénies basées sur l'ADN ribosomal nucléaire (ITS, ETS), de l'ADN chloroplastique (trnL-F, trnS-G, trnC-ycf6) et d’un locus du génome nucléaire à faible nombre de copie (CNGC4). Les données sont analysées séparément et combinées à l’aide des approches de parcimonie, bayesienne et de maximum de vraisemblance. Les données ADNnr n’ont pas permis de résoudre les relations entre les espèces polyploïdes des Eurybia. Les analyses combinées avec des loci d’ADNnr et d’ADNnr+cp ont donc été limitées à des diploïdes. Les analyses combinées ont montré une meilleure résolution et un meilleur support que les analyses séparées. La topologie de l’ADNnr+cp était la mieux résolue et supportée. La relation phylogénétique de genres appartenant au grade des eurybioïdes est comme suit : Oreostemma (Herrickia s.str. (Herrickia kingii (Eurybia (Triniteurybia - Machaerantherinae)))). Basé sur la topologie combinée de l’ADNnr+cp, nous avons effectué des analyses de biogéographie à l’aide des logiciels DIVA et LaGrange. Ces analyses ont révélé une première radiation des eurybioïdes dans l’Ouest de l’Amérique du Nord, suivi de deux migrations indépendantes dans l’Est de l’Amérique du Nord chez les Eurybia. Due au relatif manque de variabilité de l’ADNnr, l’ADNcp et CNGC4, où le triage de lignés incomplet était dominant, l'origine du grade est interprétée comme récente, possiblement du Pliocène. La diversification du groupe a été probablement favorisée par les glaciations Pléistocènes.

Approches algorithmiques pour l’inférence d’histoires de duplication en tandem avec inversions et délétions pour des familles multigéniques

[Français] Une fraction importante des génomes eucaryotes est constituée de Gènes Répétés en Tandem (GRT). Un mécanisme fondamental dans l’évolution des GRT est la recombinaison inégale durant la méiose, entrainant la duplication locale (en tandem) de segments chromosomiques contenant un ou plusieurs gènes adjacents. Différents algorithmes ont été proposés pour inférer une histoire de duplication en tandem pour un cluster de GRT. Cependant, leur utilisation est limitée dans la pratique, car ils ne tiennent pas compte d’autres événements évolutifs pourtant fréquents, comme les inversions, les duplications inversées et les délétions. Cette thèse propose différentes approches algorithmiques permettant d’intégrer ces événements dans le modèle de duplication en tandem classique. Nos contributions sont les suivantes: • Intégrer les inversions dans un modèle de duplication en tandem simple (duplication d’un gène à la fois) et proposer un algorithme exact permettant de calculer le nombre minimal d’inversions s’étant produites dans l’évolution d’un cluster de GRT. • Généraliser ce modèle pour l’étude d’un ensemble de clusters orthologues dans plusieurs espèces. • Proposer un algorithme permettant d’inférer l’histoire évolutive d’un cluster de GRT en tenant compte des duplications en tandem, duplications inversées, inversions et délétions de segments chromosomiques contenant un ou plusieurs gènes adjacents.

Phylogéographie comparée d’un système multitrophique : les parasitoïdes du genre Horismenus spp. ont-ils échappé au processus de domestication du haricot au Mexique?

Cette étude vise à comparer l’histoire évolutive des parasitoïdes du genre Horismenus (Hymenoptera: Eulophidae) à celle de leurs hôtes bruches (Coleoptera: Bruchidae) et plante hôte (Phaseolus vulgaris L.) cultivée dans le contexte d’agriculture traditionnelle, au sein de son centre de domestication Mésoaméricain. Nous avons analysé la structure génétique de 23 populations de quatre espèces de parasitoïdes au Mexique, en utilisant un fragment du gène mitochondrial COI afin de les comparer aux structures précédemment publiées des hôtes bruches et du haricot commun. Nous avons prédit que les structures génétiques des populations d’hôtes (bruches et plante) et de parasitoïdes seraient similaires puisque également influencées par la migration entremise par l’humain (HMM) étant donnée que les parasitoïdes se développent telles que les bruches à l’intérieur des haricots. Compte tenu des stratégies de manipulation reproductive utilisées par l’alpha-protéobactérie endosymbionte Wolbachia spp. pour assurer sa transmission, la structure génétique des populations de parasitoïdes inférée à partir du génome mitochondrial devrait être altérée conséquemment à la transmission conjointe des mitochondries et des bactéries lors de la propagation de l’infection dans les populations de parasitoïdes. Les populations du parasitoïde H. missouriensis sont infectées par Wolbachia spp. Tel que prédit, ces populations ne sont pas différenciées (FST = 0,06), ce qui nous empêche d’inférer sur une histoire évolutive parallèle. Contrairement aux bruches, Acanthoscelides obtectus et A. ovelatus, la HMM n'est pas un processus contemporain qui influence la structure génétique des populations du parasitoïde H. depressus, étant donné la forte différenciation (FST = 0,34) qui existe entre ses populations. La structure génétique observée chez H. depressus est similaire à celle de sa plante hôte (i.e. dispersion aléatoire historique à partir d'un pool génique ancestral très diversifié) et est probablement le résultat d’un flux génique important en provenance des populations de parasitoïdes associées aux haricots spontanées à proximité des champs cultivés. L’étude de l’histoire évolutive intégrant plusieurs niveaux trophiques s’est avérée fructueuse dans la détection des différentes réponses évolutives entre les membres du module trophique face aux interactions humaines et parasitaires, et montre la pertinence d’analyser les systèmes écologiques dans leur ensemble.

Analyse de la corrélation conditionnelle dérivée de la coévolution d’un système de trois gènes par un modèle du maximum de vraisemblance

Les gènes codant pour des protéines peuvent souvent être regroupés et intégrés en modules fonctionnels par rapport à un organelle. Ces modules peuvent avoir des composantes qui suivent une évolution corrélée pouvant être conditionnelle à un phénotype donné. Les gènes liés à la motilité possèdent cette caractéristique, car ils se suivent en cascade en réponse à des stimuli extérieurs. L’hyperthermophilie, d’autre part, est interreliée à la reverse gyrase, cependant aucun autre élément qui pourrait y être associé avec certitude n’est connu. Ceci peut être dû à un déplacement de gènes non orthologues encore non résolu. En utilisant une approche bio-informatique, une modélisation mathématique d’évolution conditionnelle corrélée pour trois gènes a été développée et appliquée sur des profils phylétiques d’archaea. Ceci a permis d’établir des théories quant à la fonction potentielle du gène du flagelle FlaD/E ainsi que l’histoire évolutive des gènes lui étant liés et ayant contribué à sa formation. De plus, une histoire évolutive théorique a été établie pour une ligase liée à l’hyperthermophilie.

Étude de l’évolution des génomes par duplications, pertes et réarrangements

La duplication est un des évènements évolutifs les plus importants, car elle peut mener à la création de nouvelles fonctions géniques. Durant leur évolution, les génomes sont aussi affectés par des inversions, des translocations (incluant des fusions et fissions de chromosomes), des transpositions et des délétions. L'étude de l'évolution des génomes est importante, notamment pour mieux comprendre les mécanismes biologiques impliqués, les types d'évènements qui sont les plus fréquents et quels étaient les contenus en gènes des espèces ancestrales. Afin d'analyser ces différents aspects de l'évolution des génomes, des algorithmes efficaces doivent être créés pour inférer des génomes ancestraux, des histoires évolutives, des relations d'homologies et pour calculer les distances entre les génomes. Dans cette thèse, quatre projets reliés à l'étude et à l'analyse de l'évolution des génomes sont présentés : 1) Nous proposons deux algorithmes pour résoudre des problèmes reliés à la duplication de génome entier : un qui généralise le problème du genome halving aux pertes de gènes et un qui permet de calculer la double distance avec pertes. 2) Nous présentons une nouvelle méthode pour l'inférence d'histoires évolutives de groupes de gènes orthologues répétés en tandem. 3) Nous proposons une nouvelle approche basée sur la théorie des graphes pour inférer des gènes in-paralogues qui considère simultanément l'information provenant de différentes espèces afin de faire de meilleures prédictions. 4) Nous présentons une étude de l'histoire évolutive des gènes d'ARN de transfert chez 50 souches de Bacillus.

Évolution moléculaire : un modèle Markov-modulé pour les processus de substitution

Les processus Markoviens continus en temps sont largement utilisés pour tenter d’expliquer l’évolution des séquences protéiques et nucléotidiques le long des phylogénies. Des modèles probabilistes reposant sur de telles hypothèses sont conçus pour satisfaire la non-homogénéité spatiale des contraintes fonctionnelles et environnementales agissant sur celles-ci. Récemment, des modèles Markov-modulés ont été introduits pour décrire les changements temporels dans les taux d’évolution site-spécifiques (hétérotachie). Des études ont d’autre part démontré que non seulement la force mais également la nature de la contrainte sélective agissant sur un site peut varier à travers le temps. Ici nous proposons de prendre en charge cette réalité évolutive avec un modèle Markov-modulé pour les protéines sous lequel les sites sont autorisés à modifier leurs préférences en acides aminés au cours du temps. L’estimation a posteriori des différents paramètres modulants du noyau stochastique avec les méthodes de Monte Carlo est un défi de taille que nous avons su relever partiellement grâce à la programmation parallèle. Des réglages computationnels sont par ailleurs envisagés pour accélérer la convergence vers l’optimum global de ce paysage multidimensionnel relativement complexe. Qualitativement, notre modèle semble être capable de saisir des signaux d’hétérogénéité temporelle à partir d’un jeu de données dont l’histoire évolutive est reconnue pour être riche en changements de régimes substitutionnels. Des tests de performance suggèrent de plus qu’il serait mieux ajusté aux données qu’un modèle équivalent homogène en temps. Néanmoins, les histoires substitutionnelles tirées de la distribution postérieure sont bruitées et restent difficilement interprétables du point de vue biologique.

Analysing artefacts dependencies to evolving software systems

Les logiciels sont en constante évolution, nécessitant une maintenance et un développement continus. Ils subissent des changements tout au long de leur vie, que ce soit pendant l'ajout de nouvelles fonctionnalités ou la correction de bogues. Lorsque les logiciels évoluent, leurs architectures ont tendance à se dégrader et deviennent moins adaptables aux nouvelles spécifications des utilisateurs. En effet, les architectures de ces logiciels deviennent plus complexes et plus difficiles à maintenir à cause des nombreuses dépendances entre les artefacts. Par conséquent, les développeurs doivent comprendre les dépendances entre les artefacts des logiciels pour prendre des mesures proactives qui facilitent les futurs changements et ralentissent la dégradation des architectures des logiciels. D'une part, le maintien d'un logiciel sans la compréhension des les dépendances entre ses artefacts peut conduire à l'introduction de défauts. D'autre part, lorsque les développeurs manquent de connaissances sur l'impact de leurs activités de maintenance, ils peuvent introduire des défauts de conception, qui ont un impact négatif sur l'évolution du logiciel. Ainsi, les développeurs ont besoin de mécanismes pour comprendre comment le changement d'un artefact impacte le reste du logiciel. Dans cette thèse, nous proposons trois contributions principales : La spécification de deux nouveaux patrons de changement et leurs utilisations pour fournir aux développeurs des informations utiles concernant les dépendances de co-changement. La spécification de la relation entre les patrons d'évolutions des artefacts et les fautes. La découverte de la relation entre les dépendances des anti-patrons et la prédisposition des différentes composantes d'un logiciel aux fautes.

Ramifications génétiques et démographiques de l'effet fondateur québécois

Les événements fondateurs et les expansions territoriales peuvent promouvoir une cascade de changements génétiques et ont ainsi pu jouer un rôle important au cours de l’histoire évolutive de l’Homme moderne. Or, chez les populations humaines, les conséquences évolutives et la dynamique démographique des processus de colonisation demeurent largement méconnues et difficiles à étudier. Dans cette thèse, nous avons utilisé les généalogies de la population fondatrice canadienne-française ainsi que des données génomiques pour étudier ces questions. Les analyses génomiques et généalogiques, remarquablement concordantes, ont dévoilé un nouveau portrait détaillé de la structure de la population du Québec, incluant un continuum de diversité génétique dans l’axe ouest/est et des sous-populations significativement différenciées. L’analyse de l’immigration fondatrice a montré que virtuellement tous les Canadiens français sont métissés. Allant à l’encontre d’une prétendue homogénéité génétique de la population, nos résultats démontrent que le peuplement des régions a engendré une rapide différentiation génétique et expliquent certaines signatures régionales de l’effet fondateur. De plus, en suivant les changements évolutifs dans les généalogies, nous avons montré que les caractéristiques des peuplements fondateurs peuvent affecter les traits liés à la fécondité et au succès reproducteur. Cette thèse offre une meilleure compréhension du patrimoine génétique du Québec et apporte des éléments de réponse sur les conséquences évolutives des événements fondateurs.

Phylogénie et biogéographie du genre Bauhinia s.l. (Leguminosae)

Bauhinia s.l. est le plus vaste genre de la tribu des Cercideae (Ceasalpinioideae, Leguminoseae), avec plus de 300 espèces. Il présente une distribution pantropicale et une grande variabilité morphologique. Ces deux caractéristiques ont limité les études taxonomiques sur le genre complet, résultant en plusieurs études taxonomiques de certains groupes seulement. En 1987, Wunderlin et al. proposent une vaste révision taxonomique de la tribu des Cercideae, basée sur des données morphologiques, et divisent le genre Bauhinia en quatre sous-genres. En 2005, Lewis et Forest publient une nouvelle classification préliminaire basée sur des données moléculaires, mais sur un échantillonnage taxonomique restreint. Leurs conclusions remettent en question le monophylétisme du genre Bauhinia et suggèrent plutôt la reconnaissance de huit genres au sein du grade Bauhinia s.l. Afin de vérifier les hypothèses de Lewis et Forest, et obtenir une vision plus claire de l’histroire de Bauhinia s.l., nous avons séquencé deux régions chloroplastiques (trnL-trnF et matK-trnK) et deux régions nucléaires (Leafy et Legcyc) pour un vaste échantillonnage représentatif des Cercideae. Une première phylogénie de la tribu a tout d’abord été réalisée à partir des séquences de trnL-trnF seulement et a confirmé le non-monoplylétisme de Bauhinia s.l., avec l’inclusion du genre Brenierea, traditionnellement reconnu comme genre frère de Bauhinia s.l. Afin de ne pas limiter notre vision de l’histoire évolutive des Cercideae à un seul type de données moléculaires et à une seule région, une nouvelle série d’analyse a été effectuée, incluant toutes les séquences chloroplastiques et nucléaires. Une phylogénie individuelle a été reconstruite pour chacune des régions du génome, et un arbre d’espèce ainsi qu’un arbre de supermatrice ont été reconstruits. Bien que certaines contradictions apparaissent entre les phylogénies, les grandes lignes de l’histoire des Cercideae ont été résolues. Bauhinia s.l. est divisée en deux lignées : les groupes Phanera et Bauhinia. Le groupe Bauhinia est constitué des genres Bauhinia s.s., Piliostigma et Brenierea. Le groupe Phanera est constitué des genres Gigasiphon, Tylosema, Lysiphyllum, Barklya, Phanera et Schnella. Les genres Cercis, Adenolobus et Griffonia sont les groupes-frères du clade Bauhinia s.l. Au minimum un événement de duplication de Legcyc a été mis en évidence pour la totalité de la tribu des Cercideae, excepté Cercis, mais plusieurs évènements sont suggérés à la fois par Legcyc et Leafy. Finalement, la datation et la reconstruction des aires ancestrales de la tribu ont été effectuées. La tribu est datée de 49,7 Ma et est originaire des régions tempérées de l’hémisphère nord, probablement autour de la mer de Thétys. La tribu s’est ensuite dispersée vers les régions tropicales sèches de l’Afrique, où la séparation des groupes Bauhinia et Phanera a eu lieu. Ces deux groupes se sont ensuite dispersés en parallèle vers l’Asie du sud-est au début du Miocène. À la même période, une dispersion depuis l’Afrique de Bauhinia s.s. a permis la diversification des espèces américaines de ce genre, alors que le genre Schnella (seul genre américain du groupe Phanera) est passé par l’Australie afin de rejoindre le continent américain. Cette dispersion vers l’Australie sera également à l’origine des genres Lysiphyllum et Barklya

Genomic architecture of adaptive color pattern divergence and convergence in Heliconius butterflies

Identifying the genetic changes driving adaptive variation in natural populations is key to understanding the origins of biodiversity. The mosaic of mimetic wing patterns in Heliconius butterflies makes an excellent system for exploring adaptive variation using next-generation sequencing. In this study, we use a combination of techniques to annotate the genomic interval modulating red color pattern variation, identify a narrow region responsible for adaptive divergence and convergence in Heliconius wing color patterns, and explore the evolutionary history of these adaptive alleles. We use whole genome resequencing from four hybrid zones between divergent color pattern races of Heliconius erato and two hybrid zones of the co-mimic Heliconius melpomene to examine genetic variation across 2.2 Mb of a partial reference sequence. In the intergenic region near optix, the gene previously shown to be responsible for the complex red pattern variation in Heliconius, population genetic analyses identify a shared 65-kb region of divergence that includes several sites perfectly associated with phenotype within each species. This region likely contains multiple cis-regulatory elements that control discrete expression domains of optix. The parallel signatures of genetic differentiation in H. erato and H. melpomene support a shared genetic architecture between the two distantly related co-mimics; however, phylogenetic analysis suggests mimetic patterns in each species evolved independently. Using a combination of next-generation sequencing analyses, we have refined our understanding of the genetic architecture of wing pattern variation in Heliconius and gained important insights into the evolution of novel adaptive phenotypes in natural populations.

Aplicació de mètodes genètics en l'avaluació de la gestió per a la conservació de poblacions de truita comuna, Salmo trutta

Durant anys, el principal mètode de gestió de les poblacions de truita comuna (Salmo trutta L.) ha estat la repoblació amb exemplars exògens. El seguiment genètic de les poblacions de truita comuna dels Pirineus orientals, realitzat en aquest tesi, indica que els al.lels procedents d'aquestes repoblacions estan conduint a una homogeneització de les poblacions naturals i a la pèrdua de la seva història evolutiva. D'aquí la importància de la detecció de la introgressió en el desenvolupament de noves estratègies de gestió i conservació de les poblacions d'aquesta espècie. En aquest treball, s'ha avaluat l'eficàcia de diferents marcadors i mètodes que ens ofereix la genètica de poblacions en la detecció de la introgressió present a les poblacions naturals. Alhora que s'ha analizat la influència que han tingut les reserves genètiques, aplicades amb posterioritat a les repoblacions, i que intenten equilibrar l'explotació i la conservació dels recursos genètics de les poblacions natives.

Overview and response to reviewers of The Singing Neanderthals

Why are humans musical? Why do people in all cultures sing or play instruments? Why do we appear to have specialized neurological apparatus for hearing and interpreting music as distinct from other sounds? And how does our musicality relate to language and to our evolutionary history? Anthropologists and archaeologists have paid little attention to the origin of music and musicality — far less than for either language or ‘art’. While art has been seen as an index of cognitive complexity and language as an essential tool of communication, music has suffered from our perception that it is an epiphenomenal ‘leisure activity’, and archaeologically inaccessible to boot. Nothing could be further from the truth, according to Steven Mithen; music is integral to human social life, he argues, and we can investigate its ancestry with the same rich range of analyses — neurological, physiological, ethnographic, linguistic, ethological and even archaeological — which have been deployed to study language. In The Singing Neanderthals Steven Mithen poses these questions and proposes a bold hypothesis to answer them. Mithen argues that musicality is a fundamental part of being human, that this capacity is of great antiquity, and that a holistic protolanguage of musical emotive expression predates language and was an essential precursor to it. This is an argument with implications which extend far beyond the mere origins of music itself into the very motives of human origins. Any argument of such range is bound to attract discussion and critique; we here present commentaries by archaeologists Clive Gamble and Iain Morley and linguists Alison Wray and Maggie Tallerman, along with Mithen's response to them. Whether right or wrong, Mithen has raised fascinating and important issues. And it adds a great deal of charm to the time-honoured, perhaps shopworn image of the Neanderthals shambling ineffectively through the pages of Pleistocene prehistory to imagine them humming, crooning or belting out a cappella harmonies as they went.