When several queens reign: how polygyny affects social life in ants

Astract: The aim of this thesis was to investigate how the presence of multiple queens (polygyny) affects social organization in colonies of the ant Formica exsecta. This is important because polygyny results in reduced relatedness among colony members and therefore reflects a potential paradox for altruistic cooperation being explained by inclusive fitness theory. The reason for this is that workers in polygynous colonies rear no longer only their siblings (high inclusive fitness gain) but also more distantly ox even unrelated brood (low or no inclusive fitness gain). All research projects conducted in this thesis are novel and significant contributions to the understanding of the social evolution of insect societies. We used a mixture of experimental and observational methodologies in laboratory and field colonies of F. exsecta to examine four important aspects of social life that are impacted by polygyny. First, we investigated the influence of queen number on colony sex allocation and found that the number of queens present in a colony significantly affects colony sex ratio investment. The data were consistent with the queen-replenishment hypothesis, which is based on the observation that newly mated queens are often recruited back to their parental nest. According to this theory, colonies containing many queens should only produce males due to local resource competition (i.e. related queens compete for common resources), whereas colonies hosting few queens benefit most from producing new queens to ensure colony survival. Second, we examined how reproduction is partitioned among nestmate queens. We detected a novel pattern of reproductive partitioning whereby a high proportion of queens were completely specialized in the production of only a subset of offspring classes produced within a colony, which might translate into great differences in reproductive success between queens. Third, we could demonstrate that F. exsecta workers indiscriminately reared highly related and unrelated brood although such nepotistic behaviour (preferential rearing of relatives) would be predicted by inclusive fitness theory. The absence of nepotism is probably best explained by its negative effects on overall colony efficiency. Finally, we conducted a detailed population genetic analysis, which revealed that the genetic population structure is different for queens and workers. Our data were best explained with queens forming family-based groups (multicolonial population structure), whereas workers from several nests seemed to be grouped into larger unites (unicolonial population structure) with workers moving freely between neighbouring nests. Altogether, the presented work significantly increased our understanding of the complex organization of polygynous social insect colonies and shows how an important life history trait such as queen number affects social organization at various levels. Résumé: Le but de cette thèse était d'étudier comment la présence de plusieurs reines par colonie (polygynie) influence la vie sociale chez la fourmi Formica exsecta. Ce sujet est important parce que la polygynie chez les insectes sociaux présente un passible paradoxe au niveau de la théorie du "fitness inclusive". Ce paradoxe est basé sur le fait que les ouvrières n'élèvent plus uniquement leurs frères et soeurs (gain de "fitness inclusive" maximale), mais également des individus moins ou pas du tout apparentés (gain de "fitness inclusive" réduit ou absent). Tous les projets de recherche présentés au cours de cette thèse apportent une meilleure compréhension et connaissance au niveau de l'organisation des colonies chez les insectes sociaux. Nous avons employé des méthodes d'observation et de laboratoire afin de mettre en évidence des aspects importants de la vie sociale chez les fourmis influencés par la polygynie. Quatre aspects ont été caractérisés : (1) l'influence du nombre de reines sur le sexe ratio produit par la colonie. Nous avons démontré que les colonies contenant beaucoup de reines produisaient rarement des reines tandis que les colonies contenant peu de reines souvent investissaient beaucoup de ressources dans la production des reines. Ces résultats sont en accord avec la "queen-replenishment hypothesis" qui est basé sur l'observation que les nouvelles reines sont recrutées dans la colonie où elles étaient nées. Cette hypothèse postule que la production des reines est défavorable dans les colonies contenant beaucoup de reines, parce que ces reines apparentées, rentrent en compétition pour des ressources communes. Au contraire, la production des reines est favorable dans des colonies contenant peu de reines afin d'assurer la survie de la colonie ; (2) comment les reines dans une colonie répartissent leur reproduction. Nous avons mis en évidence un nouveau pattern de cette répartition où une grande proportion de reines est complètement spécialisée dans la production d'un seul type de couvain ce qui probablement aboutit à des différences significatives entre reines dans le succès reproducteur ; (3) la capacité des ouvrières à discriminer un couvain de soeur d'un couvain non apparenté. Les résultats ont montré que les ouvrières ne font pas de discrimination entre le couvain de soeur et le couvain non apparenté ce qui n'est pas en accord avec la théorie de la "fitness inclusive". Cette absence de discrimination est probablement due à des effets négatifs comme par exemple la diminution de la production du couvain; (4) la structure génétique d'une population de F. exsecta. Nous avons mis en évidence que la structure génétique entre des groupes de reines est significativement différente de la structure génétique entre des groupes d'ouvrières. Les données suggèrent que les reines forment des groupes basés sur une structure familiale tandis que les ouvrières sont groupées dans des unités plus grandes.

Genetic diversity of Brazilian Pantaneiro horse and relationships among horse breeds

The objective of this work was to evaluate the genetic diversity of Brazilian Pantaneiro horse by microsatellite markers, investigate the effect of genetic bottlenecks and estimate genetic differentiation among four horse breeds. Genetic variation was estimated through allele frequencies and mean breed heterozygosity. Nei's genetic distances among the breeds Pantaneiro, Thoroughbred, Arabian, Spanish Pure Breed (Andalusian), and Uruguay Creole were calculated, and it was used to construct an UPGMA dendrogram. Clustering at different K values was calculated to infer population structure and assign individuals to populations. Nei's distances showed a minimum distance between Pantaneiro horse and Spanish Pure Breed (0.228), and similar distances from Spanish Pure Breed to Thoroughbred and to Arabian (0.355 and 0.332). It was observed a great level of diversity, clear distance from Pantaneiro horse to other breeds, and genetic uniformity within breed. It was verified a certain level of substructure of Pantaneiro horse showing no influences from the other studied breeds.

Genetic diversity in caribou linked to past and future climate change

Climate-driven range fluctuations during the Pleistocene have continuously reshaped species distribution leading to populations of contrasting genetic diversity. Contemporary climate change is similarly influencing species distribution and population structure, with important consequences for patterns of genetic diversity and species' evolutionary potential1. Yet few studies assess the impacts of global climatic changes on intraspecific genetic variation2, 3, 4, 5. Here, combining analyses of molecular data with time series of predicted species distributions and a model of diffusion through time over the past 21 kyr, we unravel caribou response to past and future climate changes across its entire Holarctic distribution. We found that genetic diversity is geographically structured with two main caribou lineages, one originating from and confined to Northeastern America, the other originating from Euro-Beringia but also currently distributed in western North America. Regions that remained climatically stable over the past 21 kyr maintained a high genetic diversity and are also predicted to experience higher climatic stability under future climate change scenarios. Our interdisciplinary approach, combining genetic data and spatial analyses of climatic stability (applicable to virtually any taxon), represents a significant advance in inferring how climate shapes genetic diversity and impacts genetic structure.

Reduced genetic diversity, increased isolation and multiple introductions of invasive giant hogweed in the western Swiss Alps.

Abstract The giant hogweed (Heracleum mantegazzianum) has successfully invaded 19 European countries as well as parts of North America. It has become a problematic species due to its ability to displace native flora and to cause public health hazards. Applying population genetics to species invasion can help reconstruct invasion history and may promote more efficient management practice. We thus analysed levels of genetic variation and population genetic structure of H. mantegazzianum in an invaded area of the western Swiss Alps as well as in its native range (the Caucasus), using eight nuclear microsatellite loci together with plastid DNA markers and sequences. On both nuclear and plastid genomes, native populations exhibited significantly higher levels of genetic diversity compared to invasive populations, confirming an important founder event during the invasion process. Invasive populations were also significantly more differentiated than native populations. Bayesian clustering analysis identified five clusters in the native range that corresponded to geographically and ecologically separated groups. In the invaded range, 10 clusters occurred. Unlike native populations, invasive clusters were characterized by a mosaic pattern in the landscape, possibly caused by anthropogenic dispersal of the species via roads and direct collection for ornamental purposes. Lastly, our analyses revealed four main divergent groups in the western Swiss Alps, likely as a consequence of multiple independent establishments of H. mantegazzianum.

Influence of biocontrol strain Pseudomonas fluorescens CHA0 and its antibiotic overproducing derivative on the diversity of root colonizing pseudomonads

Non-target effects of biocontrol strains of Pseudomonas on the population of resident pseudomonads should be assessed prior to their large scale application in the environment. The rifampicin resistant bacterium P. fluorescens CHA0-Rif and its antibiotic overproducing derivative CHA0-Rif/pME3424 were introduced into soil microcosms and the population of resident pseudomonads colonizing cucumber roots was investigated after 10 and 52 days. Both CHA0-Rif and CHA0-Rif/pME3424 displaced a part of the resident pseudomonad population after 10 days. To investigate the population structure, utilization of 10 carbon sources and production of two exoenzymes was assessed for 5600 individual pseudomonad isolates and 1700 isolates were subjected to amplified ribosomal DNA restriction analysis of the spacer region (spacer-ARDRA). After 10 days, only the proportion of pseudomonads able to degrade -tryptophan was reduced in treatments inoculated with either biocontrol strain. In parallel the phenotypic diversity was reduced. These effects were only observed 10 days after inoculation, and they were similar for inoculation with CHA0-Rif and CHA0-Rif/pME3424. Changes in the population structure of resident pseudomonads on cucumber roots during plant growth were more pronounced than changes due to the inoculants. The inoculants did not affect the genotypic diversity detected with spacer-ARDRA, but the genotypic fingerprints corresponded only partially to the phenotypic profiles. Overall CHA0-Rif had a small and transient impact on the population of resident pseudomonads and the effect was essentially the same for the genetically engineered derivative CHA0-

Estimating gene flow across hybrid zones: how reliable are microsatellites?

Due to their high polymorphism, microsatellites have become one of the most valued genetic markers in population biology. We review the first two published studies on hybrid zones of the common shrew based on microsatellites. Both reveal surprisingly high interracial gene flow. It can be shown that these are overestimates. Indeed, in classical population genetics models as F-statistics, mutation is neglected. This constitutes an acceptable assumption as long as migration is higher than mutation. However, in hybrid zones where genetic exchanges can be rare, neglecting mutation will lead to strong overestimates of migration when working with microsatellites which display mutation rates up to 10(-3). As there seems to be no straightforward way to correct for this bias, interracial gene flow estimates based on microsatellites should be taken with caution. This problem should however not conceal the enormous potential of microsatellites to unravel the genetics of hybrid zones.

Inbreeding and selection on sex ratio in the bark beetle Xylosandrus germanus.

ABSTRACT: BACKGROUND: Local Mate Competition (LMC) theory predicts a female should produce a more female-biased sex ratio if her sons compete with each other for mates. Because it provides quantitative predictions that can be experimentally tested, LMC is a textbook example of the predictive power of evolutionary theory. A limitation of many earlier studies in the field is that the population structure and mating system of the studied species are often estimated only indirectly. Here we use microsatellites to characterize the levels of inbreeding of the bark beetle Xylosandrus germanus, a species where the level of LMC is expected to be high. RESULTS: For three populations studied, genetic variation for our genetic markers was very low, indicative of an extremely high level of inbreeding (FIS = 0.88). There was also strong linkage disequilibrium between microsatellite loci and a very strong genetic differentiation between populations. The data suggest that matings among non-siblings are very rare (3%), although sex ratios from X. germanus in both the field and the laboratory have suggested more matings between non-sibs, and so less intense LMC. CONCLUSIONS: Our results confirm that caution is needed when inferring mating systems from sex ratio data, especially when a lack of biological detail means the use of overly simple forms of the model of interest.

Mapping genetic variants associated with beta-adrenergic responses in inbred mice.

β-blockers and β-agonists are primarily used to treat cardiovascular diseases. Inter-individual variability in response to both drug classes is well recognized, yet the identity and relative contribution of the genetic players involved are poorly understood. This work is the first genome-wide association study (GWAS) addressing the values and susceptibility of cardiovascular-related traits to a selective β(1)-blocker, Atenolol (ate), and a β-agonist, Isoproterenol (iso). The phenotypic dataset consisted of 27 highly heritable traits, each measured across 22 inbred mouse strains and four pharmacological conditions. The genotypic panel comprised 79922 informative SNPs of the mouse HapMap resource. Associations were mapped by Efficient Mixed Model Association (EMMA), a method that corrects for the population structure and genetic relatedness of the various strains. A total of 205 separate genome-wide scans were analyzed. The most significant hits include three candidate loci related to cardiac and body weight, three loci for electrocardiographic (ECG) values, two loci for the susceptibility of atrial weight index to iso, four loci for the susceptibility of systolic blood pressure (SBP) to perturbations of the β-adrenergic system, and one locus for the responsiveness of QTc (p<10(-8)). An additional 60 loci were suggestive for one or the other of the 27 traits, while 46 others were suggestive for one or the other drug effects (p<10(-6)). Most hits tagged unexpected regions, yet at least two loci for the susceptibility of SBP to β-adrenergic drugs pointed at members of the hypothalamic-pituitary-thyroid axis. Loci for cardiac-related traits were preferentially enriched in genes expressed in the heart, while 23% of the testable loci were replicated with datasets of the Mouse Phenome Database (MPD). Altogether these data and validation tests indicate that the mapped loci are relevant to the traits and responses studied.

Genetic associations in the detection of QTLs for wheat spike-related traits

The objective of this work was to assess the genetic diversity and population structure of wheat genotypes, to detect significant and stable genetic associations, as well as to evaluate the efficiency of statistical models to identify chromosome regions responsible for the expression of spike-related traits. Eight important spike characteristics were measured during five growing seasons in Serbia. A set of 30 microsatellite markers positioned near important agronomic loci was used to evaluate genetic diversity, resulting in a total of 349 alleles. The marker-trait associations were analyzed using the general linear and mixed linear models. The results obtained for number of allelic variants per locus (11.5), average polymorphic information content value (0.68), and average gene diversity (0.722) showed that the exceptional level of polymorphism in the genotypes is the main requirement for association studies. The population structure estimated by model-based clustering distributed the genotypes into six subpopulations according to log probability of data. Significant and stable associations were detected on chromosomes 1B, 2A, 2B, 2D, and 6D, which explained from 4.7 to 40.7% of total phenotypic variations. The general linear model identified a significantly larger number of marker-trait associations (192) than the mixed linear model (76). The mixed linear model identified nine markers associated to six traits.

RiskDiff: A web tool for the analysis of the difference due to risk and demographic factors for incidence or mortality data.

Background Analysing the observed differences for incidence or mortality of a particular disease between two different situations (such as time points, geographical areas, gender or other social characteristics) can be useful both for scientific or administrative purposes. From an epidemiological and public health point of view, it is of great interest to assess the effect of demographic factors in these observed differences in order to elucidate the effect of the risk of developing a disease or dying from it. The method proposed by Bashir and Estève, which splits the observed variation into three components: risk, population structure and population size is a common choice at practice. Results A web-based application, called RiskDiff has been implemented (available at http://rht.iconcologia.net/riskdiff.htm webcite), to perform this kind of statistical analyses, providing text and graphical summaries. Code from the implemented functions in R is also provided. An application to cancer mortality data from Catalonia is used for illustration. Conclusions Combining epidemiological with demographical factors is crucial for analysing incidence or mortality from a disease, especially if the population pyramids show substantial differences. The tool implemented may serve to promote and divulgate the use of this method to give advice for epidemiologic interpretation and decision making in public health.

The odd couple: contrasting phylogeographic patterns in two sympatric sibling species of woodlouse-­‐hunter spiders in the Canary Islands

Comparative phylogeography seeks for commonalities in the spatial demographic history of sympatric organisms to characterize the mechanisms that shaped such patterns. The unveiling of incongruent phylogeographic patterns in co-occurring species, on the other hand, may hint to overlooked differences in their life histories or microhabitat preferences. The woodlouse-hunter spiders of the genus Dysdera have undergone a major diversi cation on the Canary Islands. The species pair Dysdera alegranzaensis and Dysdera nesiotes are endemic to the island of Lanzarote and nearby islets, where they co-occur at most of their known localities. The two species stand in sharp contrast to other sympatric endemic Dysdera in showing no evidence of somatic (non-genitalic) differentiation. Phylogenetic and population genetic analyses of mitochondrial cox1 sequences from an exhaustive sample of D. alegranzaensis and D. nesiotes specimens, and additional mitochondrial (16S, L1, nad1) and nuclear genes (28S, H3) were analysed to reveal their phylogeographic patterns and clarify their phylogenetic relationships. Relaxed molecular clock models using ve calibration points were further used to estimate divergence times between species and populations. Striking differences in phylogeography and population structure between the two species were observed. Dysdera nesiotes displayed a metapopulation-like structure, while D. alegranzaensis was characterized by a weaker geographical structure but greater genetic divergences among its main haplotype lineages, suggesting more complex population dynamics. Our study con rms that co-distributed sibling species may exhibit contrasting phylogeographic patterns in the absence of somatic differentiation. Further ecological studies, however, will be necessary to clarify whether the contrasting phylogeographies may hint at an overlooked niche partitioning between the two species. In addition, further comparisons with available phylogeographic data of other eastern Canarian Dysdera endemics con rm the key role of lava ows in structuring local populations in oceanic islands and identify localities that acted as refugia during volcanic eruptions

The Imprint of Geologic History on Within‐Island Diversification of Woodlouse-­Hunter Spiders (Araneae, Dysderidae) in the Canary Islands.

Geological processes and ecological adaptation are major drivers of diversification on oceanic islands. Although diversification in these islands is often interpreted as resulting from dispersal or island hopping rather than vicariance, this may not be the case in islands with complex geological histories. The island of Tenerife, in the Canary Islands, emerged in the late Miocene as 3 precursor islands that were subsequently connected and reisolated by volcanic cycles. The spider Dysdera verneaui is endemic to the island of Tenerife, where it is widely distributed throughout most island habitats, providing an excellent model to investigate the role of physical barriers and ecological adaptation in shaping within-island diversity. Here, we present evidence that the phylogeographic patterns of this species trace back to the independent emergence of the protoislands. Molecular markers (mitochondrial genes cox1, 16S, and nad1 and the nuclear genes ITS-2 and 28S) analyzed from 100 specimens (including a thorough sampling of D. verneaui populations and additional outgroups) identify 2 distinct evolutionary lineages that correspond to 2 precursor islands, each with diagnostic genital characters indicative of separate species status. Episodic introgression events between these 2 main evolutionary lineages explain the observed incongruence between mitochondrial and nuclear markers, probably as a result of the homogenization of their ITS-2 sequence types. The most widespread lineage exhibits a complex population structure, which is compatible with either secondary contact, following connection of deeply divergent lineages, or alternatively, a back colonization from 1 precursor island to another.

The Imprint of Geologic History on Within‐Island Diversification of Woodlouse-­Hunter Spiders (Araneae, Dysderidae) in the Canary Islands.

Geological processes and ecological adaptation are major drivers of diversification on oceanic islands. Although diversification in these islands is often interpreted as resulting from dispersal or island hopping rather than vicariance, this may not be the case in islands with complex geological histories. The island of Tenerife, in the Canary Islands, emerged in the late Miocene as 3 precursor islands that were subsequently connected and reisolated by volcanic cycles. The spider Dysdera verneaui is endemic to the island of Tenerife, where it is widely distributed throughout most island habitats, providing an excellent model to investigate the role of physical barriers and ecological adaptation in shaping within-island diversity. Here, we present evidence that the phylogeographic patterns of this species trace back to the independent emergence of the protoislands. Molecular markers (mitochondrial genes cox1, 16S, and nad1 and the nuclear genes ITS-2 and 28S) analyzed from 100 specimens (including a thorough sampling of D. verneaui populations and additional outgroups) identify 2 distinct evolutionary lineages that correspond to 2 precursor islands, each with diagnostic genital characters indicative of separate species status. Episodic introgression events between these 2 main evolutionary lineages explain the observed incongruence between mitochondrial and nuclear markers, probably as a result of the homogenization of their ITS-2 sequence types. The most widespread lineage exhibits a complex population structure, which is compatible with either secondary contact, following connection of deeply divergent lineages, or alternatively, a back colonization from 1 precursor island to another.

The Imprint of Geologic History on Within‐Island Diversification of Woodlouse-­Hunter Spiders (Araneae, Dysderidae) in the Canary Islands.

Geological processes and ecological adaptation are major drivers of diversification on oceanic islands. Although diversification in these islands is often interpreted as resulting from dispersal or island hopping rather than vicariance, this may not be the case in islands with complex geological histories. The island of Tenerife, in the Canary Islands, emerged in the late Miocene as 3 precursor islands that were subsequently connected and reisolated by volcanic cycles. The spider Dysdera verneaui is endemic to the island of Tenerife, where it is widely distributed throughout most island habitats, providing an excellent model to investigate the role of physical barriers and ecological adaptation in shaping within-island diversity. Here, we present evidence that the phylogeographic patterns of this species trace back to the independent emergence of the protoislands. Molecular markers (mitochondrial genes cox1, 16S, and nad1 and the nuclear genes ITS-2 and 28S) analyzed from 100 specimens (including a thorough sampling of D. verneaui populations and additional outgroups) identify 2 distinct evolutionary lineages that correspond to 2 precursor islands, each with diagnostic genital characters indicative of separate species status. Episodic introgression events between these 2 main evolutionary lineages explain the observed incongruence between mitochondrial and nuclear markers, probably as a result of the homogenization of their ITS-2 sequence types. The most widespread lineage exhibits a complex population structure, which is compatible with either secondary contact, following connection of deeply divergent lineages, or alternatively, a back colonization from 1 precursor island to another.

The interplay of spatial and climatic landscapes in the genetic distribution of a South American parrot

AimOur aim was to understand the interplay of heterogeneous climatic and spatial landscapes in shaping the distribution of nuclear microsatellite variation in burrowing parrots, Cyanoliseus patagonus. Given the marked phenotypic differences between populations of burrowing parrots we hypothesized an important role of geographical as well climatic heterogeneity in the population structure of this species. LocationSouthern South America. MethodsWe applied a landscape genetics approach to investigate the explicit patterns of genetic spatial autocorrelation based on both geography and climate using spatial principal component analysis (sPCA). This necessitated a novel statistical estimation of the species climatic landscape, considering temperature- and precipitation-based variables separately to evaluate their weight in shaping the distribution of genetic variation in our model system. ResultsGeographical and climatic heterogeneity successfully explained molecular variance in burrowing parrots. sPCA divided the species distribution into two main areas, Patagonia and the pre-Andes, which were connected by an area of geographical and climatic transition. Moreover, sPCA revealed cryptic and conservation-relevant genetic structure: the pre-Andean populations and the transition localities were each divided into two groups, each management units for conservation. Main conclusionssPCA, a method originally developed for spatial genetics, allowed us to unravel the genetic structure related to spatial and climatic landscapes and to visualize these patterns in landscape space. These novel climatic inferences underscore the importance of our modified sPCA approach in revealing how climatic variables can drive cryptic patterns of genetic structure, making the approach potentially useful in the study of any species distributed over a climatically heterogeneous landscape.