A Method for Detecting Population Genetic Structure in Diverse, High Gene-Flow Species

Detecting small amounts of genetic subdivision across geographic space remains a persistent challenge. Often a failure to detect genetic structure is mistaken for evidence of panmixia, when more powerful statistical tests may uncover evidence for subtle geographic differentiation. Such slight subdivision can be demographically and evolutionarily important as well as being critical for management decisions. We introduce here a method, called spatial analysis of shared alleles (SAShA), that detects geographically restricted alleles by comparing the spatial arrangement of allelic co-occurrences with the expectation under panmixia. The approach is allele-based and spatially explicit, eliminating the loss of statistical power that can occur with user-defined populations and statistical averaging within populations. Using simulated data sets generated under a stepping-stone model of gene flow, we show that this method outperforms spatial autocorrelation (SA) and UST under common real-world conditions: at relatively high migration rates when diversity is moderate or high, especially when sampling is poor. We then use this method to show clear differences in the genetic patterns of 2 nearshore Pacific mollusks, Tegula funebralis (5 Chlorostoma funebralis) and Katharina tunicata, whose overall patterns of within-species differentiation are similar according to traditional population genetics analyses. SAShA meaningfully complements UST/FST, SA, and other existing geographic genetic analyses and is especially appropriate for evaluating species with high gene flow and subtle genetic differentiation.

Spatial and temporal variation in population genetic structure of wild Nile tilapia (Oreochromis niloticus) across Africa

Background: Reconstructing the evolutionary history of a species is challenging. It often depends not only on the past biogeographic and climatic events but also the contemporary and ecological factors, such as current connectivity and habitat heterogeneity. In fact, these factors might interact with each other and shape the current species distribution. However, to what extent the current population genetic structure reflects the past and the contemporary factors is largely unknown. Here we investigated spatio-temporal genetic structures of Nile tilapia (Oreochromis niloticus) populations, across their natural distribution in Africa. While its large biogeographic distribution can cause genetic differentiation at the paleo-biogeographic scales, its restricted dispersal capacity might induce a strong genetic structure at micro-geographic scales. Results: Using nine microsatellite loci and 350 samples from ten natural populations, we found the highest genetic differentiation among the three ichthyofaunal provinces and regions (Ethiopian, Nilotic and Sudano-Sahelian) (R(ST) = 0.38 - 0.69). This result suggests the predominant effect of paleo-geographic events at macro-geographic scale. In addition, intermediate divergences were found between rivers and lakes within the regions, presumably reflecting relatively recent interruptions of gene flow between hydrographic basins (R(ST) = 0.24 - 0.32). The lowest differentiations were observed among connected populations within a basin (R(ST) = 0.015 in the Volta basin). Comparison of temporal sample series revealed subtle changes in the gene pools in a few generations (F = 0 - 0.053). The estimated effective population sizes were 23 - 143 and the estimated migration rate was moderate (m similar to 0.094 - 0.097) in the Volta populations. Conclusions: This study revealed clear hierarchical patterns of the population genetic structuring of O. niloticus in Africa. The effects of paleo-geographic and climatic events were predominant at macro-geographic scale, and the significant effect of geographic connectivity was detected at micro-geographic scale. The estimated effective population size, the moderate level of dispersal and the rapid temporal change in genetic composition might reflect a potential effect of life history strategy on population dynamics. This hypothesis deserves further investigation. The dynamic pattern revealed at micro-geographic and temporal scales appears important from a genetic resource management as well as from a biodiversity conservation point of view.

Rapid parallel adaptive radiations from a single hybridogenic ancestral population

The Alpine lake whitefish (Coregonus lavaretus) species complex is a classic example of a recent radiation, associated with colonization of the Alpine lakes following the glacial retreat (less than 15 kyr BP). They have formed a unique array of endemic lake flocks, each with one to six described sympatric species differing in morphology, diet and reproductive ecology. Here, we present a genomic investigation of the relationships between and within the lake flocks. Comparing the signal between over 1000 AFLP loci and mitochondrial control region sequence data, we use phylogenetic tree-based and population genetic methods to reconstruct the phylogenetic history of the group and to delineate the principal centres of genetic diversity within the radiation. We find significant cytonuclear discordance showing that the genomically monophyletic Alpine whitefish clade arose from a hybrid swarm of at least two glacial refugial lineages. Within this radiation, we find seven extant genetic clusters centred on seven lake systems. Most interestingly, we find evidence of sympatric speciation within and parallel evolution of equivalent phenotypes among these lake systems. However, we also find the genetic signature of human-mediated gene flow and diversity loss within many lakes, highlighting the fragility of recent radiations.

Hidden Biodiversity in an Ecologically Important Freshwater Amphipod: Differences in Genetic Structure between Two Cryptic Species

Cryptic species, i.e. species that are morphologically hard to distinguish, have been detected repeatedly in various taxa and ecosystems. In order to evaluate the importance of this finding, we have to know in how far cryptic species differ in various aspects of their biology. The amphipod Gammarus fossarum is a key invertebrate in freshwater streams and contains several cryptic species. We examined the population genetic structure, genetic diversity and demographic history of two of them (type A and type B) using microsatellite markers and asked whether they show significant differences. We present results of population genetic analyses based on a total of 37 populations from the headwaters of two major European drainages, Rhine and Rhone. We found that, in both species, genetic diversity was geographically structured among and within drainages. For type A in the Rhine and type B in the Rhone, we detected significant patterns of isolation by distance. The increase of genetic differentiation with geographical distance, however, was much higher in type A than in type B. This result indicates substantial interspecific differences in population history and/or the extent of current gene flow between populations. In the Rhine, type B does not show evidence of isolation by distance, and population differentiation is relatively low across hundreds of kilometres. The majority of these populations also show signatures of recent bottlenecks. These patterns are consistent with a recent expansion of type B into the Rhine drainage. In summary, our results suggest considerable and previously unrecognized interspecific differences in the genetic structure of these cryptic keystone species.

Genetic diversity, phenotypic variation and local adaptation in the alpine landscape: case studies with alpine plant species

Plant survival in alpine landscapes is constantly challenged by the harsh and often unpredictable environmental conditions. Steep environmental gradients and patchy distribution of habitats lead to small size and spatial isolation of populations and restrict gene flow. Agricultural land use has further increased the diversity of habitats below and above the treeline. We studied the consequences of the highly structured alpine landscape for evolutionary processes in four study plants: Epilobium fleischeri, Geum reptans, Campanula thyrsoides and Poa alpina. The main questions were: (1) How is genetic diversity distributed within and among populations and is it affected by altitude, population size or land use? (2) Do reproductive traits such as allocation to sexual or vegetative reproduction vary with altitude or land use? Furthermore, we studied if seed weight increases with altitude. Within-population genetic diversity of the four species was high and mostly not related to altitude and population size. Nevertheless, genetic differentiation among populations was pronounced and strongly increasing with distance. In Poa alpina genetic diversity was affected by land use. Results suggest considerable genetic drift among populations of alpine plants. Reproductive allocation was affected by altitude and land use in Poa alpina and by succession in Geum reptans. Seed weight was usually higher in alpine species than in related lowland species. We conclude that the evolutionary potential to respond to global change is mostly intact in alpine plants, even at high altitude. Phenotypic variability is shaped by adaptive as well as by random evolutionary processes; moreover plastic responses to growth conditions seem to be crucial for survival of plants in the alpine landscape.

Global diversity, population stratification and selection of human copy number variation

In order to explore the diversity and selective signatures of duplication and deletion human copy number variants (CNVs), we sequenced 236 individuals from 125 distinct human populations. We observed that duplications exhibit fundamentally different population genetic and selective signatures than deletions and are more likely to be stratified between human populations. Through reconstruction of the ancestral human genome, we identify megabases of DNA lost in different human lineages and pinpoint large duplications that introgressed from the extinct Denisova lineage now found at high frequency exclusively in Oceanic populations. We find that the proportion of CNV base pairs to single nucleotide variant base pairs is greater among non-Africans than it is among African populations, but we conclude that this difference is likely due to unique aspects of non-African population history as opposed to differences in CNV load.

Genetic composition, genetic diversity, and small-scale environmental variation matter for the experimental reintroduction of a rare plant

Aims Reintroduction has become an important tool for the management of endangered plant species. We tested the little-explored effects of small-scale environmental variation, genotypic composition (i.e. identity of genotypes), and genotypic diversity on the population survival of the regionally rare clonal plant Ranunculus reptans. For this species of periodically inundated lakeshores genetic differentiation had been reported between populations and between short-flooded and long-flooded microsites within populations.Methods We established 306 experimental test populations at a previously unoccupied lake shore, comprising either monocultures of 32 genotypes, mixtures of genotypes within populations or mixtures of genotypes between populations. In 2000, three years after planting out at the experimental site, a long-lasting flood caused the death of half of the experimental populations. In 2003, an extreme drought resulted in the lowest summer water levels ever measured.Important findings Despite these climatic extremes, 27 of the established populations survived until the end of the experiment in December 2003. The success of experimental populations largely differed between microsites. Moreover, the success of genotype monocultures depended on genotype and source population. Genetic differentiation between microsites played a minor role for the success of reintroduction. After the flood, populations planted with genotypes from different source populations increased in abundance, whereas populations with genotypes from single source populations and genotype monocultures decreased. We conclude that sources for reintroductions need to be selected carefully. Moreover, mixtures of plants from different populations appear to be the best choice for successful reintroduction, at least in unpredictably varying environments.

Low but contrasting neutral genetic differentiation shaped by winter temperature in European great tits

Gene flow is usually thought to reduce genetic divergence and impede local adaptation by homogenising gene pools between populations. However, evidence for local adaptation and phenotypic differentiation in highly mobile species, experiencing high levels of gene flow, is emerging. Assessing population genetic structure at different spatial scales is thus a crucial step towards understanding mechanisms underlying intraspecific differentiation and diversification. Here, we studied the population genetic structure of a highly mobile species – the great tit Parus major – at different spatial scales. We analysed 884 individuals from 30 sites across Europe including 10 close-by sites (< 50 km), using 22 microsatellite markers. Overall we found a low but significant genetic differentiation among sites (FST = 0.008). Genetic differentiation was higher, and genetic diversity lower, in south-western Europe. These regional differences were statistically best explained by winter temperature. Overall, our results suggest that great tits form a single patchy metapopulation across Europe, in which genetic differentiation is independent of geographical distance and gene flow may be regulated by environmental factors via movements related to winter severity. This might have important implications for the evolutionary trajectories of sub-populations, especially in the context of climate change, and calls for future investigations of local differences in costs and benefits of philopatry at large scales.

Hybridization between distant lineages increases adaptive variation during a biological invasion: stickleback in Switzerland

The three-spined stickleback is a widespread Holarctic species complex that radiated from the sea into freshwaters after the retreat of the Pleistocene ice sheets. In Switzerland, sticklebacks were absent with the exception of the far northwest, but different introduced populations have expanded to occupy a wide range of habitats since the late 19th century. A well-studied adaptive phenotypic trait in sticklebacks is the number of lateral plates. With few exceptions, freshwater and marine populations in Europe are fixed for either the low plated phenotype or the fully plated phenotype, respectively. Switzerland, in contrast, harbours in close proximity the full range of phenotypic variation known from across the continent. We addressed the phylogeographic origins of Swiss sticklebacks using mitochondrial partial cytochrome b and control region sequences. We found only five different haplotypes but these originated from three distinct European regions, fixed for different plate phenotypes. These lineages occur largely in isolation at opposite ends of Switzerland, but co-occur in a large central part. Across the country, we found a strong correlation between a microsatellite linked to the high plate ectodysplasin allele and the mitochondrial haplotype from a region where the fully plated phenotype is fixed. Phylogenomic and population genomic analysis of 481 polymorphic amplified fragment length polymorphism loci indicate genetic admixture in the central part of the country. The same part of the country also carries elevated within-population phenotypic variation. We conclude that during the recent invasive range expansion of sticklebacks in Switzerland, adaptive and neutral between-population genetic variation was converted into within-population variation, raising the possibility that hybridization between colonizing lineages contributed to the ecological success of sticklebacks in Switzerland.

Mutational spectrum and linkage disequilibrium patterns at the ornithine transcarbamylase gene (OTC)

Ornithine transcarbamylase (OTC; EC 2.1.3.3) is a hepatic enzyme involved in ammonia elimination via the urea cycle. Since the sequence of the OTC gene was reported many types of mutations continue to be found in OTC deficiency patients, continuing to increase the already wide mutational spectrum known for this gene. In this study we present the clinical, biochemical and molecular features of thirteen late-onset OTC deficiency patients. Mutations were identified in all these patients, among which six were novel point substitutions (L59R, A137P, L148S, Y176L, L186P, and K210N) and one was a 2-bp deletion at exon 4 (341-342delAA). In addition, a de novo genomic deletion of maternal origin encompassing exons 1 to 5 was also identified by the analysis of LD patterns using intragenic polymorphic markers. This work exemplifies the potential value of population genetic studies for the detection of large deletions.