Ecology and evolution of social-organization: insights from fire ants and other highly eusocial insects

Social organisms exhibit conspicuous intraspecific variation in all facets of their social organization. A prominent example of such variation in the highly eusocial Hymenoptera is differences in the number of reproductive queens per colony, Differences in queen number in ants are associated with differences in a host of reproductive and social traits, including queen phenotype and breeding strategy, mode of colony reproduction, and pattern of sex allocation. We examine the causes and consequences of changes in colony queen number and associated traits using the fire ant Solenopsis invicta as a principal model. Ecological constraints on mode of colony founding may act as important selective forces causing the evolution of queen number in this and many other ants, with social organization generally perpetuated across generations by means of the social environment molding appropriate queen phenotypes and reproductive strategies. Shifts in colony queen number have profound effects on genetic structure within nests and may also influence genetic structure at higher levels (aggregations of nests or local demes) because of the association of queen number with particular mating and dispersal habits. Divergence of breeding habits between populations with different social organizations has the potential to promote genetic differentiation between these social variants. Thus, evolution of social organization can be important in generating intrinsic selective regimes that channel subsequent social evolution and in initiating the development of significant population genetic structure, including barriers to gene flow important in cladogenesis.

Chromosome localization of microsatellite markers in the shrews of the Sorex araneus group.

The extremely high rate of karyotypic evolution that characterizes the shrews of the Sorex araneus group makes this group an exceptionally interesting model for population genetics and evolutionary studies. Here, we attempted to map 46 microsatellite markers at the chromosome arm level using flow-sorted chromosomes from three karyotypically different taxa of the Sorex araneus group (S. granarius and the chromosome races Cordon and Novosibirsk of S. araneus). The most likely localizations were provided for 35 markers, among which 25 were each unambiguously mapped to a single locus on the corresponding chromosomes in the three taxa, covering the three sexual chromosomes (XY1Y2) and nine of the 18 autosomal arms of the S. araneus group. The results provide further evidence for a high degree of conservation in genome organization in the S. araneus group despite the presence of numerous Robertsonian rearrangements. These markers can therefore be used to compare the genetic structure among taxa of the S. araneus group at the chromosome level and to study the role of chromosomal rearrangements in the genetic diversification and speciation process of this group.

Isolation and characterization of highly polymorphic microsatellite loci in the bladder campion, Silene vulgaris (Caryophyllaceae)

This study reports the isolation and characterization of seven highly polymorphic microsatellite loci in Silene vulgaris (Caryophyllaceae). The loci were isolated from two libraries constructed from genomic DNA enriched for CA and GA repeats. These markers yielded nine to 40 alleles per locus (mean 22.1) in a survey of 45 individuals from a single population located in the western Swiss Alps. Average observed heterozygosity ranged from 16.2 to 77.4%. These microsatellite loci should be valuable tools for studying fine-scale genetic structure.

Origin and expansion of the allotetraploid Aegilops geniculata, a wild relative of wheat.

*This study reconstructs the phylogeography of Aegilops geniculata, an allotetraploid relative of wheat, to discuss the impact of past climate changes and recent human activities (e.g. the early expansion of agriculture) on the genetic diversity of ruderal plant species. *We combined chloroplast DNA (cpDNA) sequencing, analysed using statistical parsimony network, with nonhierarchical K-means clustering of amplified fragment length polymorphism (AFLP) genotyping, to unravel patterns of genetic structure across the native range of Ae. geniculata. The AFLP dataset was further explored by measurement of the regional genetic diversity and the detection of isolation by distance patterns. *Both cpDNA and AFLP suggest an eastern Mediterranean origin of Ae. geniculata. Two lineages have spread independently over northern and southern Mediterranean areas. Northern populations show low genetic diversity but strong phylogeographical structure among the main peninsulas, indicating a major influence of glacial cycles. By contrast, low genetic structuring and a high genetic diversity are detected in southern Mediterranean populations. Finally, we highlight human-mediated dispersal resulting in substantial introgression between resident and migrant populations. *We have shown that the evolutionary trajectories of ruderal plants can be similar to those of wild species, but are interfered by human activities, promoting range expansions through increased long-distance dispersal and the creation of suitable habitats.

Sex-specific estimates of dispersal show female philopatry and male dispersal in a promiscuous amphibian, the alpine salamander (Salamandra atra).

Amphibians display wide variations in life-history traits and life cycles that should prove useful to explore the evolution of sex-biased dispersal, but quantitative data on sex-specific dispersal patterns are scarce. Here, we focused on Salamandra atra, an endemic alpine species showing peculiar life-history traits. Strictly terrestrial and viviparous, the species has a promiscuous mating system, and females reproduce only every 3 to 4 years. In the present study, we provide quantitative estimates of asymmetries in male vs. female dispersal using both field-based (mark-recapture) and genetic approaches (detection of sex-biased dispersal and estimates of migration rates based on the contrast in genetic structure across sexes and age classes). Our results revealed a high level of gene flow among populations, which stems exclusively from male dispersal. We hypothesize that philopatric females benefit from being familiar with their natal area for the acquisition and defence of an appropriate shelter, while male dispersal has been secondarily favoured by inbreeding avoidance. Together with other studies on amphibians, our results indicate that a species' mating system alone is a poor predictor of sex-linked differences in dispersal, in particular for promiscuous species. Further studies should focus more directly on the proximate forces that favour or limit dispersal to refine our understanding of the evolution of sex-biased dispersal in animals.

Mating frequency of ant queens with alternative dispersal strategies, as revealed by microsatellite analysis of sperm

In social Hymenoptera (ants, bees, and wasps), the number of males that mate with the same queen affects social and genetic organization of the colony. However, the selective forces leading to single mating in certain conditions and multiple mating in others remain enigmatic. In this study, I investigated whether queens of the wood ant Formica paralugubris adopting different dispersal strategies varied in their mating frequency (the number of males with whom they mated). The frequency of multiple mating was determined by using microsatellite markers to genotype the sperm stored in the spermatheca of queens, and the validity of this method was confirmed by analysing mother-offspring combinations obtained from experimental single-queen colonies. Dispersing queens, which may found new colonies, did not mate with more males than queens that stayed within polygynous colonies, where the presence of numerous reproductive individuals ensured a high level of genetic diversity. Hence, this study provides no support to the hypotheses that multiple mating is beneficial because it increases genetic variability within colonies. Most of the F. paralugubris queens mated with a single male, whatever their dispersal strategy and life history. Moreover, multiple mating had little effect on colony genetic structure: the effective mating frequency was 1.11 when calculated from within-brood relatedness, and 1.13 when calculated from the number of mates detected in the sperm. Hence, occasional multiple mating by F. paralugubris queens may have no adaptive significance.

Restricted gene flow at specific parts of the shrew genome in chromosomal hybrid zones.

The species and races of the shrews of the Sorex araneus group exhibit a broad range of chromosomal polymorphisms. European taxa of this group are parapatric and form contact or hybrid zones that span an extraordinary variety of situations, ranging from absolute genetic isolation to almost free gene flow. This variety seems to depend for a large part on the chromosome composition of populations, which are primarily differentiated by various Robertsonian fusions of a subset of acrocentric chromosomes. Previous studies suggested that chromosomal rearrangements play a causative role in the speciation process. In such models, gene flow should be more restricted for markers on chromosomes involved in rearrangements than on chromosomes common in both parent species. In the present study, we address the possibility of such differential gene flow in the context of two genetically very similar but karyotypically different hybrid zones between species of the S. araneus group using microsatellite loci mapped to the chromosome arm level. Interspecific genetic structure across rearranged chromosomes was in general larger than across common chromosomes. However, the difference between the two classes of chromosomes was only significant in the hybrid zone where the complexity of hybrids is expected to be larger. These differences did not distinguish populations within species. Therefore, the rearranged chromosomes appear to affect the reproductive barrier between karyotypic species, although the strength of this effect depends on the complexity of the hybrids produced.

Predicting the deleterious effects of mutation load in fragmented populations.

Human-induced habitat fragmentation constitutes a major threat to biodiversity. Both genetic and demographic factors combine to drive small and isolated populations into extinction vortices. Nevertheless, the deleterious effects of inbreeding and drift load may depend on population structure, migration patterns, and mating systems and are difficult to predict in the absence of crossing experiments. We performed stochastic individual-based simulations aimed at predicting the effects of deleterious mutations on population fitness (offspring viability and median time to extinction) under a variety of settings (landscape configurations, migration models, and mating systems) on the basis of easy-to-collect demographic and genetic information. Pooling all simulations, a large part (70%) of variance in offspring viability was explained by a combination of genetic structure (F(ST)) and within-deme heterozygosity (H(S)). A similar part of variance in median time to extinction was explained by a combination of local population size (N) and heterozygosity (H(S)). In both cases the predictive power increased above 80% when information on mating systems was available. These results provide robust predictive models to evaluate the viability prospects of fragmented populations.

Characterization of thirteen microsatellite markers in river and brook lampreys (Lampetra fluviatilis and L-planeri)

We describe the development based on 454 pyrosequencing technology of thirteen microsatellite markers for two closely related species of lamprey: Lampetra fluviatilis and L. planeri. The number of alleles per locus ranged from 2 to 5 in L. fluviatilis and from 2 to 6 in L. planeri. Gene diversity ranged from 0.062 to 0.718 in L. fluviatilis and from 0.322 to 0.677 in L. planeri. These markers will be helpful to study population genetic structure of both species and resolve their taxonomic status as separate species or ecotypes of a single species.

Inferring recent migration rates from individual genotypes.

We present a novel and straightforward method for estimating recent migration rates between discrete populations using multilocus genotype data. The approach builds upon a two-step sampling design, where individual genotypes are sampled before and after dispersal. We develop a model that estimates all pairwise backwards migration rates (m(ij), the probability that an individual sampled in population i is a migrant from population j) between a set of populations. The method is validated with simulated data and compared with the methods of BayesAss and Structure. First, we use data for an island model and then we consider more realistic data simulations for a metapopulation of the greater white-toothed shrew (Crocidura russula). We show that the precision and bias of estimates primarily depend upon the proportion of individuals sampled in each population. Weak sampling designs may particularly affect the quality of the coverage provided by 95% highest posterior density intervals. We further show that it is relatively insensitive to the number of loci sampled and the overall strength of genetic structure. The method can easily be extended and makes fewer assumptions about the underlying demographic and genetic processes than currently available methods. It allows backwards migration rates to be estimated across a wide range of realistic conditions.

Mites as biological tags of their hosts.

Movements and spatial distribution of host populations are expected to shape the genetic structure of their parasite populations. Comparing the genetic patterns of both interacting species may improve our understanding of their evolutionary history. Moreover, genetic analyses of parasites with horizontal transmission may serve as indicators of historical events or current demographic processes that are not apparent in the genetic signature of their hosts. Here, we compared mitochondrial variation in populations of the ectoparasitic mite Spinturnix myoti with the genetic pattern of its host, the Maghrebian bat Myotis punicus in North Africa and in the islands of Corsica and Sardinia. Mite mitochondrial differentiation among populations was correlated with both host mitochondrial and nuclear differentiation, suggesting spatial co-differentiation of the lineages of the two interacting species. Therefore our results suggest that parasite dispersal is exclusively mediated by host movements, with open water between landmasses as a main barrier for host and parasite dispersal. Surprisingly the unique presence of a continental European mite lineage in Corsica was inconsistent with host phylogeographical history and strongly suggests the former presence of European mouse-eared bats on this island. Parasites may thus act as biological tags to reveal the presence of their now locally extinct host.

Phylogeography and recolonization of the Swiss Alps by the Valais shrew (Sorex antinorii), inferred with autosomal and sex-specific markers.

Using one male-inherited, one female-inherited and eight biparentally inherited markers, we investigate the population genetic structure of the Valais shrew (Sorex antinorii) in the Swiss Alps. Bayesian analysis on autosomal microsatellites suggests a clear genetic differentiation between two groups of populations. This geographically based structure is consistent with two separate postglacial recolonization routes of the species into Switzerland from Italian refugia after the last Pleistocene glaciations. Sex-specific markers also confirm genetic structuring among western and eastern areas, since very few haplotypes for either Y chromosome or mtDNA genome are shared between the two regions. Overall, these results suggest that two already well-differentiated genetic lineages colonized the Swiss Alps and came into secondary contact in the Rhône Valley. Low level of admixture between the two lineages is likely explained by the mountainous landscape structure of lateral valleys orthogonal to the main Rhône valley.

Comparative phylogeography in a specific and obligate pollination antagonism.

In specific and obligate interactions the nature and abundance of a given species can have important effects on the survival and population dynamics of associated organisms. In a phylogeographic framework, we therefore expect that the fates of organisms interacting specifically are also tightly interrelated. Here we investigate such a scenario by analyzing the genetic structures of species interacting in an obligate plant-insect pollination lure-and-trap antagonism, involving Arum maculatum (Araceae) and its specific psychodid (Diptera) visitors Psychoda phalaenoides and Psycha grisescens. Because the interaction is asymmetric (i.e., only the plant depends on the insect), we expect the genetic structure of the plant to be related with the historical pollinator availability, yielding incongruent phylogeographic patterns between the interacting organisms.Using insect mtDNA sequences and plant AFLP genome fingerprinting, we inferred the large-scale phylogeographies of each species and the distribution of genetic diversities throughout the sampled range, and evaluated the congruence in their respective genetic structures using hierarchical analyses of molecular variances (AMOVA). Because the composition of pollinator species varies in Europe, we also examined its association with the spatial genetic structure of the plant.Our findings indicate that while the plant presents a spatially well-defined genetic structure, this is not the case in the insects. Patterns of genetic diversities also show dissimilar distributions among the three interacting species. Phylogeographic histories of the plant and its pollinating insects are thus not congruent, a result that would indicate that plant and insect lineages do not share the same glacial and postglacial histories. However, the genetic structure of the plant can, at least partially, be explained by the type of pollinators available at a regional scale. Differences in life-history traits of available pollinators might therefore have influenced the genetic structure of the plant, the dependent organism, in this antagonistic interaction.

Multilocus phylogeography of the common lizard Zootoca vivipara at the Ibero-Pyrenean suture zone reveals lowland barriers and high-elevation introgression.

BACKGROUND: The geographic distribution of evolutionary lineages and the patterns of gene flow upon secondary contact provide insight into the process of divergence and speciation. We explore the evolutionary history of the common lizard Zootoca vivipara (= Lacerta vivipara) in the Iberian Peninsula and test the role of the Pyrenees and the Cantabrian Mountains in restricting gene flow and driving lineage isolation and divergence. We also assess patterns of introgression among lineages upon secondary contact, and test for the role of high-elevation trans-mountain colonisations in explaining spatial patterns of genetic diversity. We use mtDNA sequence data and genome-wide AFLP loci to reconstruct phylogenetic relationships among lineages, and measure genetic structure RESULTS: The main genetic split in mtDNA corresponds generally to the French and Spanish sides of the Pyrenees as previously reported, in contrast to genome-wide AFLP data, which show a major division between NW Spain and the rest. Both types of markers support the existence of four distinct and geographically congruent genetic groups, which are consistent with major topographic barriers. Both datasets reveal the presence of three independent contact zones between lineages in the Pyrenean region, one in the Basque lowlands, one in the low-elevation mountains of the western Pyrenees, and one in the French side of the central Pyrenees. The latter shows genetic evidence of a recent, high-altitude trans-Pyrenean incursion from Spain into France. CONCLUSIONS: The distribution and age of major lineages is consistent with a Pleistocene origin and a role for both the Pyrenees and the Cantabrian Mountains in driving isolation and differentiation of Z. vivipara lineages at large geographic scales. However, mountain ranges are not always effective barriers to dispersal, and have not prevented a recent high-elevation trans-Pyrenean incursion that has led to asymmetrical introgression among divergent lineages. Cytonuclear discordance in patterns of genetic structure and introgression at contact zones suggests selection may be involved at various scales. Suture zones are important areas for the study of lineage formation and speciation, and our results show that biogeographic barriers can yield markedly different phylogeographic patterns in different vertebrate and invertebrate taxa.

Phylogéographie des principales vipères européennes (Vipera ammodytes, V. aspis et V. berus), structuration génétique et multipaternité chez Vipera berus

Résumé: Les vipères du genre Vipera sont des serpents venimeux distribués dans la totalité du Paléarctique. Malgré cette répartition considérable, elles sont extrêmement menacées, leur déclin étant principalement dû à la destruction et à la fragmentation de leur habitat ainsi qu'à la persécution humaine. Afin d'apporter de nouveaux éléments dans le contexte de la protection de ce groupe de reptiles, nous avons utilisé durant ce travail de thèse différents marqueurs moléculaires pour étudier la structuration génétique à petite et à large échelle chez trois espèces appartenant au genre Vipera. La première étude, une phylogéographie moléculaire de la vipère ammodytes (Vipera ammodytes), a montré dans l'ensemble de l'aire de répartition une forte structuration génétique provenant d'isolements antérieures au Pléistocène. La présence d'un nombre important de clades dans le centre des Balkans suggère que cette région a fourni de nombreux refuges isolés durant les glaciations. Ces dernières ont également eu un impact considérable sur la diversité génétique au sein de la majorité des clades, suite à d'importants goulots d'étranglement durant le Pléistocène. L'étude de la phylogéographie de la vipère aspic (Vipera aspis) a montré une différenciation génétique entre les populations présentes de chaque côté des Alpes, mais également une forte structuration interne avec la mise en évidence d'un refuge en France. Cette étude est la première à établir clairement l'utilisation d'un refuge français pour un vertébré terrestre. La troisième partie de cette thèse a étudié la phylogéographie de la vipère péliade (Vipera berus), espèce cible de ce travail. En plus de la mise en évidence d'un groupe génétique inattendu (localisé dans le nord de l'Italie, le sud de l'Autriche, le nord de la Slovénie et l'extrême sud-est de la Suisse), la variabilité génétique au sein du groupe nordique (comprenant les animaux de l'entier de l'aire de répartition de l'espèce à l'exception des individus du groupe italien et les animaux provenant des Balkans) est suffisamment importante pour conclure à l'utilisation de refuges glaciaires nordiques durant les dernières glaciations, en complément des refuges habituellement décrits pour la majorité des espèces animales (soit les péninsules ibérique, italienne et balakanique). Ces résultats nous ont conduit à effectuer une étude morphologique (quatrième partie) comparant les vipères péliades du "clade italien" et du "clade nordique" décrits ci-dessus. Seules de petites différences morphologiques ont pu être mises en évidence, malgré une séparation de ces groupes estimée à plus d'un million d'années. Une étude à plus petite échelle, centrée sur le Massif jurassien et certaines populations alpines et françaises, a été entreprise afin d'estimer leur diversité génétique et d'évaluer la structuration génétique entre les populations à l'aide de marqueurs microsatellites (cinquième partie). Une importante structuration a été observée entre les populations distantes de plus de 3 kilomètres, la structuration entre les populations plus proches étant plus limitée. De plus, une diversité génétique plus faible dans les populations jurassiennes et alpines comparativement aux populations du massif central et de la côte atlantique a été constatée, probablement due à une perte de diversité génétique lors de la recolonisation post-glaciaire. La sixième étude s'est intéressée au succès reproducteur des mâles de vipères péliades en conditions naturelles. Une corrélation entre la taille des mâles et leur succès reproducteur a été relevée, les individus de plus grande taille ayant un succès reproducteur plus élevé. Le taux de multipaternité a aussi été investigué, démontrant que la proportion de pontes issues de plusieurs pères est élevée (69%) malgré la faible densité de vipères observée sur le site étudié. Finalement, aucun lien entre le nombre de pères au sein d'une ponte et la mortalité des jeunes à la naissance n'a pu être mis en évidence, contrastant avec des travaux précédents. En conclusion, l'observation de la structuration très marquée chez les vipères péliades devrait permettre d'affiner les méthodes de protection de l'espèce dans le massif jurassien. A plus large échelle, l'importante structuration génétique observée chez les vipères ammodytes, aspic et péliade résultant de l'utilisation de nombreux refuges glaciaires, complémentaires aux refuges habituellement utilisés par les espèces animales, démontre l'intérêt de l'analyse phylogéographique des reptiles pour la compréhension des phénomènes de colonisation et d' extinction des populations durant la fin du Tertiaire et le Quaternaire. La mise en évidence chez les différentes espèces de vipères étudiées de nombreux groupes génétiques distincts (ESUs) devrait conduire à des modifications de la taxonomie ainsi qu'au statut de protection de ces espèces. Abstract: The vipers of the genus Vipera are venomous snakes widespread throughout the Palaearctic regions. Despite a large distribution area, several species are extremely threatened, especially due to the destruction and fragmentation of their habitats, as well as by human persecution. In order to increase the knowledge on these species and to improve their protection, several molecular markers have been used to investigate the genetic structure on small and large scales, within three species of the genus Vipera. The first study, a molecular phylogeography of the nose-horned viper (Vipera ammodytes), showed a considerable structuring throughout the distribution area, due to isolation into refugia before the Pleistocene. A high number of clades in the centre of the Balkans suggests that this region harboured numerous isolated glacial refugia during the last glaciation. Moreover, low genetic diversity within several clades implies that most populations of nose-horned vipers have suffered bottlenecks during the Pleistocene. The study of the phylogeography of the asp viper (Vipera aspis) showed genetic differentiation between populations on each side of the Alps, as well as considerable internal genetic structure, suggesting the use of a glacial refugium in France. This study is the first to establish firmly the occurrence of a French refugia for a terrestrial vertebrate. The third part of this work involved a phylogeographic study of the adder (Vipera berus), the target species of this thesis. Three clades were revealed: a Balkan clade (corresponding to the subspecies V. b. sachalinensis), an unexpected Italian clade (limited to northern Italy, southern Austria, northern Slovenia and southeasternmost corner of Switzerland) and a Northern clade clade (including adders of the whole distribution area excepted animals from the Balkan and the Italian clades). The genetic variability within the Northern clade is sufficiently high to conclude that a northern glacial refugia during the last glaciation, in addition to those refugia already described for the main species (Iberian, Italian and Balkan peninsula). These results motivated a morphological study (part four) comparing the adders from the Italian and the Northern clades describe above. Only small morphological differences have been found, despite the split between these two clades have taken place more than 1 million years ago. A study on a local scale, focused on the Jura Mountains, on a few populations in the Alps and France was, performed to estimate the genetic diversity and the genetic structure between populations using microsatellite markers (part five). Considerable structure was observed between populations separated by more than 3 kilometres, whereas the structure between closer populations is less marked. Moreover, lower genetic diversity in the populations from Jura Mountains and Alps was noticed compared to populations from Massif Central of Atlantic coast. Such loss of genetic variation probably followed post-glacial recolonisation. The sixth study focused on the reproductive success of male adders in the wild. A positive correlation between body length and reproductive success was observed. Multiple paternity was also observed in most of clutches (69%) despite the low density of adders in the study area. Finally, no relationship was found between the number of fathers in a clutch and the survival of offspring at birth, contradicting previous studies. To conclude, the observation of a significant genetic structure in Vipera berus will enable recommendations to be made to improve protection of this species in the Jura Mountain. On a larger scale, the considerable genetic structure found within Vipera ammdoytes, V. aspis and V. berus, resulting from isolation in additional glacial refugia to those already described for other species, demonstrates the relevance of phylogeographic studies of reptiles to better understand the colonisation and disappearance during the last Tertiary and the Quaternary. The observation of several groups of evolutionary significant units (ESUs) within the three studied species might lead to a revision of the taxonomy, as well as their conservation status.