Genetic clusters and sex-biased gene flow in a unicolonial Formica ant.

BACKGROUND: Animal societies are diverse, ranging from small family-based groups to extraordinarily large social networks in which many unrelated individuals interact. At the extreme of this continuum, some ant species form unicolonial populations in which workers and queens can move among multiple interconnected nests without eliciting aggression. Although unicoloniality has been mostly studied in invasive ants, it also occurs in some native non-invasive species. Unicoloniality is commonly associated with very high queen number, which may result in levels of relatedness among nestmates being so low as to raise the question of the maintenance of altruism by kin selection in such systems. However, the actual relatedness among cooperating individuals critically depends on effective dispersal and the ensuing pattern of genetic structuring. In order to better understand the evolution of unicoloniality in native non-invasive ants, we investigated the fine-scale population genetic structure and gene flow in three unicolonial populations of the wood ant F. paralugubris. RESULTS: The analysis of geo-referenced microsatellite genotypes and mitochondrial haplotypes revealed the presence of cryptic clusters of genetically-differentiated nests in the three populations of F. paralugubris. Because of this spatial genetic heterogeneity, members of the same clusters were moderately but significantly related. The comparison of nuclear (microsatellite) and mitochondrial differentiation indicated that effective gene flow was male-biased in all populations. CONCLUSION: The three unicolonial populations exhibited male-biased and mostly local gene flow. The high number of queens per nest, exchanges among neighbouring nests and restricted long-distance gene flow resulted in large clusters of genetically similar nests. The positive relatedness among clustermates suggests that kin selection may still contribute to the maintenance of altruism in unicolonial populations if competition occurs among clusters.

Microsatellite markers reveal shallow genetic differentiation between cohorts of the common sea urchin Paracentrotus lividus (Lamarck) in northwest Mediterranean

Temporal variability was studied in the common sea urchin Paracentrotus lividus through the analysis of the genetic composition of three yearly cohorts sampled over two consecutive springs in a locality in northwestern Mediterranean. Individuals were aged using growth ring patterns observed in tests and samples were genotyped for five microsatellite loci. No reduction of genetic diversity was observed relative to a sample of the adult population from the same location or within cohorts across years. FST and amova results indicated that the differentiation between cohorts is rather shallow and not significant, as most variability is found within cohorts and within individuals. This mild differentiation translated into estimates of effective population size of 90100 individuals. When the observed excess of homozygotes was taken into account, the estimate of the average number of breeders increased to c. 300 individuals. Given our restricted sampling area and the known small-scale heterogeneity in recruitment in this species, our results suggest that at stretches of a few kilometres of shoreline, large numbers of progenitors are likely to contribute to the larval pool at each reproduction event. Intercohort variation in our samples is six times smaller than spatial variation between adults of four localities in the western Mediterranean. Our results indicate that, notwithstanding the stochastic events that take place during the long planktonic phase and during the settlement and recruitment processes, reproductive success in this species is high enough to produce cohorts genetically diverse and with little differentiation between them. Further research is needed before the link between genetic structure and underlying physical and biological processes can be well established.

Genetic and karyotypic structure in the shrews of the Sorex araneus group: are they independent?

The species of the common shrew (Sorex araneus) group are morphologically very similar but exhibit high levels of karyotypic variation. Here we used genetic variation at 10 microsatellite markers in a data set of 212 individuals mostly sampled in the western Alps and composed of five karyotypic taxa (Sorex coronatus, Sorex antinorii and the S. araneus chromosome races Cordon, Bretolet and Vaud) to investigate the concordance between genetic and karyotypic structure. Bayesian analysis confirmed the taxonomic status of the three sampled species since individuals consistently grouped according to their taxonomical status. However, introgression can still be detected between S. antinorii and the race Cordon of S. araneus. This observation is consistent with the expected low karyotypic complexity of hybrids between these two taxa. Geographically based cryptic substructure was discovered within S. antinorii, a pattern consistent with the different postglaciation recolonization routes of this species. Additionally, we detected two genetic groups within S. araneus notwithstanding the presence of three chromosome races. This pattern can be explained by the probable hybrid status of the Bretolet race but also suggests a relatively low impact of chromosomal differences on genetic structure compared to historical factors. Finally, we propose that the current data set (available at http://www.unil.ch/dee/page7010_en.html#1) could be used as a reference by those wanting to identify Sorex individuals sampled in the western Alps.

Scale-specific sex-biased dispersal in the Valais shrew unveiled by genetic variation on the Y chromosome, autosomes, and mitochondrial DNA.

We investigated sex specificities in the evolutionary processes shaping Y chromosome, autosomes, and mitochondrial DNA patterns of genetic structure in the Valais shrew (Sorex antinorii), a mountain dwelling species with a hierarchical distribution. Both hierarchical analyses of variance and isolation-by-distance analyses revealed patterns of population structure that were not consistent across maternal, paternal, and biparentally inherited markers. Differentiation on a Y microsatellite was lower than expected from the comparison with autosomal microsatellites and mtDNA, and it was mostly due to genetic variance among populations within valleys, whereas the opposite was observed on other markers. In addition, there was no pattern of isolation by distance for the Y, whereas there was strong isolation by distance on mtDNA and autosomes. We use a hierarchical island model of coancestry dynamics to discuss the relative roles of the microevolutionary forces that may induce such patterns. We conclude that sex-biased dispersal is the most important driver of the observed genetic structure, but with an intriguing twist: it seems that dispersal is strongly male biased at large spatial scale, whereas it is mildly biased in favor of females at local scale. These results add to recent reports of scale-specific sex-biased dispersal patterns, and emphasize the usefulness of the Y chromosome in conjunction with mtDNA and autosomes to infer sex specificities.

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.

Linking genetic and species diversities in a community of freshwater gastropods

Summary Biodiversity is usually studied through species or genetic diversities. To date, these two levels of diversity have remained the independent .fields of investigations of community ecologists and population geneticists. However, recent joint analyses of species and genetic diversities have suggested that common processes may underlie the two levels. Positive correlations between species diversity and genetic diversity may arise when the effects of drift and migration overwhelm selective effects. The first goal of this thesis was to make a joint investigation of the patterns of species and genetic diversity in a community of freshwater gastropods living in a floodplain habitat. The second goal was to determine, as far as possible, the relative influences of the processes underlying the patterns observed at each level. In chapter 2 we investigate the relative influences of different evolutionary forces in shaping the genetic structure of Radix balthica populations. Results revealed that the structure inferred using quantitative traits was lower or equal to the one inferred using neutral molecular markers. Consequently, the pattern of structure observed could be only due to random drift, possibly to uniform selection, but definitely not to selection for local optima. In chapter 3, we analyze the temporal variation of species and genetic diversities in five localities. An extended period of drought occurred at the end of the study period leading to decay of both species and genetic diversities. This parallel loss of diversity following a natural perturbation highlighted the role sometimes predominant of random drift over selection on patterns of biodiversity in a floodplain habitat. In chapter 4, we compare the spatial genetic structures of two sympatric species: Radix balthica and Planorbis carinatus. We found that R. balthica populations are weakly structured and have moderate to high values of gene diversity. In contrast, P. carinatus populations are highly structured and poorly diverse. Then we measured correlations between various indices of species and genetic diversity using genetic data .from the two species. We found only one significant correlation: between species richness and gene diversity of P. carinatus. This result highlights the .need to use genetic date from more than one species to infer correlations between species and genetic diversities. Overall, this thesis provided new insights into the common processes underlying patterns of species and genetic diversity. Résumé La biodiversité est généralement étudiée au niveau de la diversité génétique ou spécifique. Ces deux niveaux sont restés jusqu'à maintenant les domaines d'investigation séparés des généticiens des populations et des écologistes des communautés. Cependant, des analyses conjointes des diversités génétique et spécifique ont récemment suggéré que des processus similaires pouvaient influencer ces deux niveaux. Des corrélations positives entre les diversités génétique et spécifique pourraient être dues aux effets de migration et de dérive qui dominent les effets sélectifs. Le premier but de cette thèse était de faire une étude conjointe des diversités génétique et spécifique dans une communauté de gastéropodes d'eau douce. Le second objectif était de déterminer les influences relatives des différents processus liés à chaque niveau de diversité. Dans le chapitre 2 nous cherchons à déterminer quelles forces évolutives influencent la structure génétique de quatre populations de Radix balthica. La structure mesurée sur des traits quantitatifs s'est révélée être plus faible ou égale à celle mesurée avec des marqueurs moléculaires neutres. La structure observée pourrait ainsi être due uniquement à la dérive génétique, potentiellement à la sélection uniforme, mais en aucun cas à la sélection locale pour différents optima. Dans le chapitre 3 nous analysons la variation temporelle des diversités génétique et spécifique dans cinq localités. Une récente période de sécheresse a causé une diminution parallèle des deux niveaux de diversité. Cette perturbation à mis en évidence le rôle parfois prépondérant de la dérive par rapport à celui de la sélection dans le déterminisme de la biodiversité dans un écosytème alluvial. Dans le chapitre 4, nous comparons la structure génétique spatiale de deux espèces vivant en sympatrie : Radix balthica et Planorbis carinatus. Les populations de R. balthica sont peu structurées et présentent un niveau de diversité relativement élevé alors que celles de P. carinatus sont fortement structurées et peu diversifiées. Nous avons ensuite mesuré différentes corrélations entre les diversités génétique et spécifique, mais la seule relation significative a été trouvée entre la richesse spécifique et la diversité génétique de P. carinatus. Ainsi, cette thèse a permis de découvrir de nouveaux aspects des processus qui influencent en parallèle la diversité aux niveaux génétique et spécifique.

Geographic patterns of genetic variation in a broadly distributed marine vertebrate: new insights into loggerhead turtle stock structure from expanded mitochondrial DNA sequences

Previous genetic studies have demonstrated that natal homing shapes the stock structure of marine turtle nesting populations. However, widespread sharing of common haplotypes based on short segments of the mitochondrial control region often limits resolution of the demographic connectivity of populations. Recent studies employing longer control region sequences to resolve haplotype sharing have focused on regional assessments of genetic structure and phylogeography. Here we synthesize available control region sequences for loggerhead turtles from the Mediterranean Sea, Atlantic, and western Indian Ocean basins. These data represent six of the nine globally significant regional management units (RMUs) for the species and include novel sequence data from Brazil, Cape Verde, South Africa and Oman. Genetic tests of differentiation among 42 rookeries represented by short sequences (380 bp haplotypes from 3,486 samples) and 40 rookeries represented by long sequences (~800 bp haplotypes from 3,434 samples) supported the distinction of the six RMUs analyzed as well as recognition of at least 18 demographically independent management units (MUs) with respect to female natal homing. A total of 59 haplotypes were resolved. These haplotypes belonged to two highly divergent global lineages, with haplogroup I represented primarily by CC-A1, CC-A4, and CC-A11 variants and haplogroup II represented by CC-A2 and derived variants. Geographic distribution patterns of haplogroup II haplotypes and the nested position of CC-A11.6 from Oman among the Atlantic haplotypes invoke recent colonization of the Indian Ocean from the Atlantic for both global lineages. The haplotypes we confirmed for western Indian Ocean RMUs allow reinterpretation of previous mixed stock analysis and further suggest that contemporary migratory connectivity between the Indian and Atlantic Oceans occurs on a broader scale than previously hypothesized. This study represents a valuable model for conducting comprehensive international cooperative data management and research in marine ecology.

Population structure in a highly pelagic seabird, the Cory's shearwater (Calonectris diomedea): an examination of genetics, morphology and ecology

Increasing evidence suggests oceanic traits may play a key role in the genetic structuring of marine organisms. Whereas genetic breaks in the open ocean are well known in fishes and marine invertebrates, the importance of marine habitat characteristics in seabirds remains less certain. We investigated the role of oceanic transitions versus population genetic processes in driving population differentiation in a highly vagile seabird, the Cory"s shearwater, combining molecular, morphological and ecological data from 27 breeding colonies distributed across the Mediterranean (Calonectris diomedea diomedea) and the Atlantic (C. d. borealis). Genetic and biometric analyses showed a clear differentiation between Atlantic and Mediterranean Cory"s shearwaters. Ringing-recovery data indicated high site fidelity of the species, but we found some cases of dispersal among neighbouring breeding sites (<300 km) and a few long distance movements (>1000 km) within and between each basin. In agreement with this, comparison of phenotypic and genetic data revealed both current and historical dispersal events. Within each region, we did not detect any genetic substructure among archipelagos in the Atlantic, but we found a slight genetic differentiation between western and eastern breeding colonies in the Mediterranean. Accordingly, gene flow estimates suggested substantial dispersal among colonies within basins. Overall, genetic structure of the Cory"s shearwater matches main oceanographic breaks (Almería-Oran Oceanic Front and Siculo-Tunisian Strait), but spatial analyses suggest that patterns of genetic differentiation are better explained by geographic rather than oceanographic distances. In line with previous studies, genetic, phenotypic and ecological evidence supported the separation of Atlantic and Mediterranean forms, suggesting the 2 taxa should be regarded as different species.

Effect of race, genetic population structure, and genetic models in two-locus association studies: clustering of functional renin-angiotensin system gene variants in hypertension association studies

Previous genetic association studies have overlooked the potential for biased results when analyzing different population structures in ethnically diverse populations. The purpose of the present study was to quantify this bias in two-locus association studies conducted on an admixtured urban population. We studied the genetic structure distribution of angiotensin-converting enzyme insertion/deletion (ACE I/D) and angiotensinogen methionine/threonine (M/T) polymorphisms in 382 subjects from three subgroups in a highly admixtured urban population. Group I included 150 white subjects; group II, 142 mulatto subjects, and group III, 90 black subjects. We conducted sample size simulation studies using these data in different genetic models of gene action and interaction and used genetic distance calculation algorithms to help determine the population structure for the studied loci. Our results showed a statistically different population structure distribution of both ACE I/D (P = 0.02, OR = 1.56, 95% CI = 1.05-2.33 for the D allele, white versus black subgroup) and angiotensinogen M/T polymorphism (P = 0.007, OR = 1.71, 95% CI = 1.14-2.58 for the T allele, white versus black subgroup). Different sample sizes are predicted to be determinant of the power to detect a given genotypic association with a particular phenotype when conducting two-locus association studies in admixtured populations. In addition, the postulated genetic model is also a major determinant of the power to detect any association in a given sample size. The present simulation study helped to demonstrate the complex interrelation among ethnicity, power of the association, and the postulated genetic model of action of a particular allele in the context of clustering studies. This information is essential for the correct planning and interpretation of future association studies conducted on this population.

Genetic variation and population structure of the endangered Hyacinth Macaw (Anodorhynchus hyacinthinus): implications for conservation

The Hyacinth Macaw (Anodorhynchus hyacinthinus) is one of 14 endangered species in the family Psittacidae occurring in Brazil, with an estimated total population of 6,500 specimens. We used nuclear molecular markers (single locus minisatellites and microsatellites) and 472 bp of the mitochondrial DNA control region to characterize levels of genetic variability in this species and to assess the degree of gene flow among three nesting sites in Brazil (Pantanal do Abobral, Pantanal de Miranda and Piaui). The origin of five apprehended specimens was also investigated. The results suggest that, in comparison to other species of parrots, Hyacinth Macaws possess relatively lower genetic variation and that individuals from two different localities within the Pantanal (Abobral and Miranda) belong to a unique interbreeding population and are genetically distinct at nuclear level from birds from the state of Piaui. The analyses of the five apprehended birds suggest that the Pantanal is not the source of birds for illegal trade, but their precise origin could not be assigned. The low genetic variability detected in the Hyacinth Macaw does not seem to pose a threat to the survival of this species. Nevertheless, habitat destruction and nest poaching are the most important factors negatively affecting their populations in the wild. The observed genetic structure emphasizes the need of protection of Hyacinth Macaws from different regions in order to maintain the genetic diversity of this species.

Genetic diversity and population structure of pumas (Puma concolor) in southeastern Brazil: implications for conservation in a human-dominated landscape

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Phylogeography and genetic population structure of Caribbean sharpnose shark Rhizoprionodon porosus

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Genetic diversity and population structure of Musa accessions in ex situ conservation

Abstract Background Banana cultivars are mostly derived from hybridization between wild diploid subspecies of Musa acuminata (A genome) and M. balbisiana (B genome), and they exhibit various levels of ploidy and genomic constitution. The Embrapa ex situ Musa collection contains over 220 accessions, of which only a few have been genetically characterized. Knowledge regarding the genetic relationships and diversity between modern cultivars and wild relatives would assist in conservation and breeding strategies. Our objectives were to determine the genomic constitution based on Internal Transcribed Spacer (ITS) regions polymorphism and the ploidy of all accessions by flow cytometry and to investigate the population structure of the collection using Simple Sequence Repeat (SSR) loci as co-dominant markers based on Structure software, not previously performed in Musa. Results From the 221 accessions analyzed by flow cytometry, the correct ploidy was confirmed or established for 212 (95.9%), whereas digestion of the ITS region confirmed the genomic constitution of 209 (94.6%). Neighbor-joining clustering analysis derived from SSR binary data allowed the detection of two major groups, essentially distinguished by the presence or absence of the B genome, while subgroups were formed according to the genomic composition and commercial classification. The co-dominant nature of SSR was explored to analyze the structure of the population based on a Bayesian approach, detecting 21 subpopulations. Most of the subpopulations were in agreement with the clustering analysis. Conclusions The data generated by flow cytometry, ITS and SSR supported the hypothesis about the occurrence of homeologue recombination between A and B genomes, leading to discrepancies in the number of sets or portions from each parental genome. These phenomenons have been largely disregarded in the evolution of banana, as the “single-step domestication” hypothesis had long predominated. These findings will have an impact in future breeding approaches. Structure analysis enabled the efficient detection of ancestry of recently developed tetraploid hybrids by breeding programs, and for some triploids. However, for the main commercial subgroups, Structure appeared to be less efficient to detect the ancestry in diploid groups, possibly due to sampling restrictions. The possibility of inferring the membership among accessions to correct the effects of genetic structure opens possibilities for its use in marker-assisted selection by association mapping.

Novel tools for conservation genetics in marine fish: population structure and evolution of the European hake (Merluccius merluccius) inferred by SNP variation and applications to traceability

The research presented in my PhD thesis is part of a wider European project, FishPopTrace, focused on traceability of fish populations and products. My work was aimed at developing and analyzing novel genetic tools for a widely distributed marine fish species, the European hake (Merluccius merluccius), in order to investigate population genetic structure and explore potential applications to traceability scenarios. A total of 395 SNPs (Single Nucleotide Polymorphisms) were discovered from a massive collection of Expressed Sequence Tags, obtained by high-throughput sequencing, and validated on 19 geographic samples from Atlantic and Mediterranean. Genome-scan approaches were applied to identify polymorphisms on genes potentially under divergent selection (outlier SNPs), showing higher genetic differentiation among populations respect to the average observed across loci. Comparative analysis on population structure were carried out on putative neutral and outlier loci at wide (Atlantic and Mediterranean samples) and regional (samples within each basin) spatial scales, to disentangle the effects of demographic and adaptive evolutionary forces on European hake populations genetic structure. Results demonstrated the potential of outlier loci to unveil fine scale genetic structure, possibly identifying locally adapted populations, despite the weak signal showed from putative neutral SNPs. The application of outlier SNPs within the framework of fishery resources management was also explored. A minimum panel of SNP markers showing maximum discriminatory power was selected and applied to a traceability scenario aiming at identifying the basin (and hence the stock) of origin, Atlantic or Mediterranean, of individual fish. This case study illustrates how molecular analytical technologies have operational potential in real-world contexts, and more specifically, potential to support fisheries control and enforcement and fish and fish product traceability.