Population structure, history and gene flow in a group of closely related land snails: genetic variation in Partula from the Society Islands of the Pacific

Previous studies of Partula land snails from the Society Islands, French Polynesia, have shown that populations within species are highly differentiated in terms of their morphology, behaviour, ecology and molecular genetic variation. Despite this level of variability, differences between species are sometimes small, possibly reflecting the fact that reproductive isolation is not always complete and there exists the opportunity for genetic exchange between taxa through hybridization. The present study uses sequence data from a mitochondrial gene to further investigate genetic variation in Society Island Partula. Most populations are found in this study to be highly differentiated, but within individual species there seems to be no simple relationship either between genetic distance and geographical proximity, or between variation in mitochondria and that in allozymes or morphological characteristics. Among species there appears to be no simple correlation between degrees of reproductive isolation and genetic relatedness according to mitochondrial DNA. The results suggest that past events as well as ongoing drift and selection may have been important in affecting patterns of variation. Similarities among species at specific localities suggest that there must have been some genetic exchange in the past, although this may not necessarily reflect ongoing rates of hybridization. The discrepancy between results for different markers probably reflects the differential effects of drift and selection on mitochondrial and nuclear genes.     

Change and stability in a steep morph‐frequency cline in the snail Cepaea nemoralis (L.) over 43 years

Populations of the polymorphic land snail Cepaea nemoralis (L.) from Deepdale, Derbyshire, UK, sampled in 1965–67, showed a pattern of area effects, with steep clines among groups of populations differing in shell colour and banding morph frequencies. In 2010, most of these populations were resampled. In particular, a continuous transect made in 1967 of 42 quadrats (18.34 × 18.34 m) across a steep cline in several morph frequencies was completely resampled. In the dale as a whole, yellow shells had increased in frequency. In the transect, the frequencies of banding morphs showed no significant changes, although colour morphs showed some changes. Pink shells had increased in frequency in a section in which scrub had developed, and brown shells had increased in frequency in the area in which they had originally been at the highest frequency. In each case, the selection coefficients were of the order of 4%. Yellow had increased elsewhere. Nevertheless, both in the dale as a whole and in the transect, the pattern of geographical change in morph frequencies had remained essentially the same. Estimates of migration based on previous studies of marked snails and on modelling of the effect of drift and migration suggest that, regardless of whether the cline is a product of differential selection or of the gradual merging of previously separate founding populations, it has been in existence for a long time, and that migration occurs over greater distances than estimated from direct observation on marked snails. Although we can demonstrate that selection has occurred, the origin and maintenance of the cline and others similar to it remain in doubt; the development and maintenance of polymorphism in this species may require consideration of several processes operating on different time scales. © 2012 The Linnean Society of London     

Polymorphism in relation to habitat in the snail Cepaea hortensis in Iceland.

Polymorphism for shell colour and banding pattern in Cepaea hortensis was studied in a confined area in south Iceland. Morph freqquencies can be related to habitat. Yellow unbanded snails are more frequent in grassland and herb meadows thatn in "darker" habitats such as in ddense Angelica . Fused banding is relatively more frequent in "daarker" habitats. Predation by birds is not known to occur and rodent predation in winter or genetic drift cannot explain the observed correlations. Habitats differ in their microclimate and it is suggested that climatic selection is important. Differences in morph frequencies between juvenile and adult snails support this view.     

Gene flow and melanism in Lake Erie garter snake populations

Melanistic garter snakes ( Thamnophis sirtalis ) are unusually common near Lake Erie, apparently because selection for thermoregulatory ability in cool lake-shore habitats (which favours melanistic morphs) outweighs selection for crypsis (which favours striped morphs). However, morph frequencies are highly variable among sites, suggesting that random genetic drift also influences colour pattern. In an effort to better understand the evolutionary processes influencing garter snake colour patterns, we estimated Fx and Nm (the number of migrants per generation) among island and mainland populations from patterns of allozymic variation detected using electrophoresis. Estimates of Nm were high, ranging from 2.7 to 37.6 between pairs of study sites and making it unlikely that differences in morph frequencies among sites were solely the result of random genetic drift. Furthermore, differences in F _st estimates between colour pattern (a one-locus two-allele trait) and allozyme loci suggest that colour pattern alleles are not in Hardy-Weinberg equilibrium, most likely as a result of natural selection. Comparison of allozymic data from Lake Erie with those from more distant sites suggests that gene flow occurs over long distances in T. sirtalis.     

Population genetics study of hemoglobinopathies in Uzbekistan. II. Population dynamics of hemoglobinopathies

It was shown that on comparing variability of selective neutral genetic marker systems with that of the beta-thalassemia system for the populations of different hierarchical level, the relative importance of selection and genetic drift could be evaluated. The genetic differentiation of the beta-thalassemia gene frequencies in elementary populations (villages) could be solely explained by genetic drift. On the other hand, the differentiation of district populations (the sizes of the populations being 10(6] for beta-thalassemia gene frequencies could be explained by selection forces. This is supported by the fact that the genetic distances and FST values are only significant for the beta-thalassemia gene and not for the neutral genetic systems, when the district populations are compared.     

Discrepancies in the estimation of gene flow in partula

Methods for estimating gene flow (Nm) from genetic data should provide important insights into the dynamics of natural populations. If they are to be used with confidence, however, the methods must be shown to produce valid results. Estimates of Nm have been obtained for the snails Partula taeniata and Partula suturalis, based on F(ST) and on the frequencies of private alleles, p(1). Jackknifing was used to reduce the bias of estimates and to obtain confidence limits. The estimates derived from F(ST) are consistent with the low vagility of snails, and with direct field studies of gene flow in P. taeniata. In contrast, the estimates derived from p(1) were up to seven times as large, less precise and less consistent. Although the underlying causes of these discrepancies are not clear, the results suggest that F(ST) is the more reliable indirect estimator of gene flow, at least for Partula.     

Gene flow and melanism in garter snakes revisited: a comparison of molecular markers and island vs. coalescent models

Within populations, the stochastic effect of genetic drift and deterministic effect of natural selection are potentially weakened or altered by gene flow among populations. The influence of gene flow on Lake Erie populations of the common garter snake has been of particular interest because of a discontinuous colour pattern polymorphism (striped vs. melanistic) that is a target of natural selection. We reassessed the relative contributions of gene flow and genetic drift using genetic data and population size estimates. We compared all combinations of two marker systems and two analytical approaches to the estimation of gene flow rates: allozymes (data previously published), microsatellite DNA (new data), the island model ( F _ST-based approach), and a coalescence-based approach. For the coalescence approach, mutation rates and sampling effects were also investigated. While the two markers produced similar results, gene flow based on F _ST was considerably higher (Nm > 4) than that from the coalescence-based method (Nm < 1). Estimates of gene flow are likely to be inflated by lack of migration-drift equilibrium and changing population size. Potentially low rates of gene flow (Nm < 1), small population size at some sites, and positive correlations of number of microsatellite DNA alleles and island size and between M , mean ratio of number of alleles to range in allele size, and island size suggest that in addition to selection, random genetic drift may influence colour pattern frequencies. © 2003 The Linnean Society of London, Biological Journal of the Linnean Society , 2003, 79, 389–399.     

Altitude,language, and class I HLA allele frequencies in Papua New Guinea

Class I HLA gene frequencies show considerable variation over short geographical distances in Papua New Guinea. Hypotheses to account for this invoke natural selection, population structure, the pattern of population movement, or past demographic changes. To determine the role of the various factors in shaping this distribution, we have studied correlations between HLA-based genetic distances, geographical distances, altitude, and linguistic differences in Papua New Guinea. Linguistic differences at the family or stock level within the Trans-New Guinea Phylum generally correspond to genetic differences. However, on the basis of their HLA gene frequencies, speakers of Austronesian (AN) languages do not form a distinct group of populations. Linguistic variation and spatial autocorrelation do not fully account for the altitudinal cline differences noted in gene frequencies, particularly at the HLA-A locus. We propose that the distribution of HLA gene frequencies in Papua New Guinea is partially under the control of selection operating differentially along the altitude gradient. © 1994 Wiley-Liss, Inc.     

A new method to uncover signatures of divergent and stabilizing selection in quantitative traits

While it is well understood that the pace of evolution depends on the interplay between natural selection, random genetic drift, mutation, and gene flow, it is not always easy to disentangle the relative roles of these factors with data from natural populations. One popular approach to infer whether the observed degree of population differentiation has been influenced by local adaptation is the comparison of neutral marker gene differentiation (as reflected in FST) and quantitative trait divergence (as reflected in QST). However, this method may lead to compromised statistical power, because FST and QST are summary statistics which neglect information on specific pairs of populations, and because current multivariate tests of neutrality involve an averaging procedure over the traits. Further, most FST-QST comparisons actually replace QST by its expectation over the evolutionary process and are thus theoretically flawed. To overcome these caveats, we derived the statistical distribution of population means generated by random genetic drift and used the probability density of this distribution to test whether the observed pattern could be generated by drift alone. We show that our method can differentiate between genetic drift and selection as a cause of population differentiation even in cases with FST=QST and demonstrate with simulated data that it disentangles drift from selection more accurately than conventional FST-QST tests especially when data sets are small.     

Contrasting patterns of body shape and neutral genetic divergence in marine and lake populations of threespine sticklebacks

Comparisons of neutral marker and quantitative trait divergence can provide important insights into the relative roles of natural selection and neutral genetic drift in population differentiation. We investigated phenotypic and genetic differentiation among Fennoscandian threespine stickleback (Gasterosteus aculeatus) populations, and found that the highest degree of differentiation occurred between sea and freshwater habitats. Within habitats, morphological divergence was highest among the different freshwater populations. Pairwise phenotypic and neutral genetic distances among populations were positively correlated, suggesting that genetic drift may have contributed to the morphological differentiation among habitats. On the other hand, the degree of phenotypic differentiation (PST) clearly surpassed the neutral expectation set by FST, suggesting a predominant role for natural selection over genetic drift as an explanation for the observed differentiation. However, separate PST/FST comparisons by habitats revealed that body shape divergence between lake and marine populations, and even among marine populations, can be strongly influenced by natural selection. On the other hand, genetic drift can play an important role in the differentiation among lake populations.     

Negative frequency‐dependent selection maintains shell banding polymorphisms in two marine snails (Littorina fabalis and Littorina saxatilis)

The presence of shell bands is common in gastropods. Both the marine snails Littorina fabalis and Lttorina saxatilis are polymorphic for this trait. Such polymorphism would be expected to be lost by the action of genetic drift or directional selection, but it appears to be widespread at relatively constant frequencies. This suggests it is maintained by balancing selection on the trait or on a genetically linked trait. Using long time series of empirical data, we compared potential effects of genetic drift and negative frequency‐dependent selection (NFDS) in the two species. The contribution of genetic drift to changes in the frequency of bands in L. fabalis was estimated using the effective population size estimated from microsatellite data, while the effect of genetic drift in L. saxatilis was derived from previously published study. Frequency‐dependent selection was assessed by comparing the cross‐product estimator of fitness with the frequency of the polymorphism across years using a regression analysis. Both studied species showed patterns of NFDS. In addition, in L. fabalis, contributions from genetic drift could explain some of the changes in banding frequency. Overdominance and heterogeneous selection did not fit well to our data. The possible biological explanations resulting in the maintenance of the banding polymorphism are discussed.     

Evidence of a reproductive barrier between two forms of the marine periwinkle Littorina fabalis (Gastropoda)

Studies of allozyme variation may reveal unexpected patterns of genetic variation which challenge earlier conclusions of species delimitations based on morphological data. However, allozyme variation alone may not be sufficient to resolve this kind of problem. For example, populations of the marine intertidal snail Littorina fabalis (=Littorina mariae) from wave exposed parts end from protected parts of the same shores are distinguished by different alleles of arginine kinase (Ark) while indifferent, or very nearly so, in another 29 loci. Intermediate populations have large deficiencies of exposed/sheltered heterozygote classes of Ark and we have earlier suggested habitat-related selection in this locus as the explanation. In this study we estimated growth rate of individual snails of different Ark-genotypes in three different habitats (exposed, sheltered and intermediate). In all habitats the snails homozygous for alleles of ‘exposed’ type grew faster and matured at a larger size than did snails homozygous for alleles of ‘sheltered’ types. This relationship was indirectly confirmed in three additional sites of intermediate exposure where exposed AA-genotypes dominated among large (>8 mm) snails while the sheltered genotypes dominated among small (<5 mm) snails of truly sympatric samples. We furthermore found small differences in allele frequencies of two other loci (Pgi and Pgm-2) and in shell colour frequencies, comparing sympatric snails of exposed and sheltered Ark-homozygotes. Although we found no signs of habitat-related selection among snails of different Ark-genotype, or selection against heterozygotes, we cannot reject selection in Ark, as our experiments only covered one island, one season and grown-up snails. The coupling between allozyme and phenotypic characters in strictly sympatric samples of snails suggests the presence of two gene pools. Perhaps the large and small forms of L. fabalis represent very closely related cryptic taxa. However, introgression between them seems a possible explanation for the striking similarities in the vast majority of morphological and allozyme characters.     

Geographic patterns of genetic differentiation and plumage colour variation are different in the pied flycatcher (Ficedula hypoleuca)

The pied flycatcher is one of the most phenotypically variable bird species in Europe. The geographic variation in phenotypes has often been attributed to spatial variation in selection regimes that is associated with the presence or absence of the congeneric collared flycatcher. Spatial variation in phenotypes could however also be generated by spatially restricted gene flow and genetic drift. We examined the genetic population structure of pied flycatchers across the breeding range and applied the phenotypic Q _ST ( P _ST)– F _ST approach to detect indirect signals of divergent selection on dorsal plumage colouration in pied flycatcher males. Allelic frequencies at neutral markers were found to significantly differ among populations breeding in central and southern Europe whereas northerly breeding pied flycatchers were found to be one apparently panmictic group of individuals. Pairwise differences between phenotypic ( P _ST) and neutral genetic distances ( F _ST) were positively correlated after removing the most differentiated Spanish and Swiss populations from the analysis, suggesting that genetic drift may have contributed to the observed phenotypic differentiation in some parts of the pied flycatcher breeding range. Differentiation in dorsal plumage colouration however greatly exceeded that observed at neutral genetic markers, which indicates that the observed pattern of phenotypic differentiation is unlikely to be solely maintained by restricted gene flow and genetic drift.     

History or current selection? A molecular analysis of 'area effects' in the land snail Cepaea nemoralis.

We have used molecular variation in microsatellite and mitochondrial DNA to throw light on the origins of enigmatic geographical patterns, known as 'area effects', in the shell polymorphisms of the land snail Cepaea nemoralis. Our aim was to assess the relative importance of recent selection and historical events in the formation of these patterns. On the Marlborough Downs in Wiltshire, England, the 'type locality' for area effects, the frequencies of microsatellites are significantly associated with the frequencies of alleles for shell banding. A less clear association is found between microsatellites and shell colour. Mitochondrial haplotypes show no significant relationships. Although the correlated geographical patterns could be the results of random genetic drift from an initially uniform array of populations, the magnitudes of the patterns, and of the correlations between them, seem too strong to have arisen by drift since the last glaciation. Our results suggest that invasions from refugia have been the most important factors in forming area effects.     

Prehistoric inter-archipelago trading of Polynesian tree snails leaves a conservation legacy

Inter-archipelago exchange networks were an important aspect of prehistoric Polynesian societies. We report here a novel genetic characterization of a prehistoric exchange network involving an endemic Pacific island tree snail, Partula hyalina. It occurs in the Society (Tahiti only), Austral and Southern Cook Islands. Our genetic data, based on museum, captive and wild-caught samples, establish Tahiti as the source island. The source lineage is polymorphic in shell coloration and contains a second nominal species, the dark-shelled Partula clara, in addition to the white-shelled P. hyalina. Prehistoric inter-island introductions were non-random: they involved white-shelled snails only and were exclusively inter-archipelago in scope. Partulid shells were commonly used in regional Polynesian jewellery, and we propose that the white-shelled P. hyalina, originally restricted to Tahiti, had aesthetic value throughout these archipelagoes. Demand within the Society Islands could be best met by trading dead shells, but a low rate of inter-archipelago exchange may have prompted the establishment of multiple founder populations in the Australs and Southern Cooks. The alien carnivorous land snail Euglandina rosea has recently devastated populations of all 61 endemic species of Society Island partulid snails. Southern Cooks and Australs P. hyalina now represent the only unscathed wild populations remaining of this once spectacular land snail radiation.     

Distribution patterns of morph frequencies in the snail Cepaea hortensis in Iceland

In Iceland Cepaea hortensis Müll. reaches its northern limit of distribution occurring in restricted areas along the south coast. Samples of snails collected along steep cliffs in two such restricted areas were effectively monomorphic for shell colour (yellow), but were polymorphic for banding pattern and revealed strong area effects. In both localities samples from sites facing east contained varying frequencies of unbanded snails whereas in ecologically similar sites facing south unbanded were absent. In an experimental population released at a south-facing site unbanded snails failed to become established whereas banded individuals survived and bred. It is suggested that the area effects are maintained through non-visual (most likely climatic) selection though other possibilities cannot be ruled out.     

Plant species diversity and composition of experimental grasslands affect genetic differentiation of Lolium perenne populations

Contrasting hypotheses exist about the relationship between plant species diversity and genetic diversity. However, experimental data of species diversity effects on genetic differentiation among populations are lacking. To address this, Lolium perenne was sown with an equal number of seeds in 78 experimental grasslands (Jena Experiment) varying in species richness (1, 2, 4, 8 to 16) and functional group richness and composition (1-4; grasses, legumes, small herbs, tall herbs). Population sizes were determined 4years after sowing, and single-nucleotide polymorphism (SNP) DNA markers based on bulk samples of up to 100 individuals per population were applied. Genetic distances between the field populations and the initially sown seed population increased with sown species richness. The degree of genetic differentiation from the original seed population was largely explained by actual population sizes, which suggests that genetic drift was the main driver of differentiation. Weak relationships among relative allele frequencies and species diversity or actual population sizes, and a positive correlation between actual population sizes and expected heterozygosity also supported the role of genetic drift. Functional composition had additional effects on genetic differentiation of L. perenne populations, indicating a selection because of genotype-specific interactions with other species. Our study supports that genetic diversity is likely to be lower in plant communities with a higher number of interspecific competitors. Negative effects of species richness on population sizes may increase the probability of genetic drift, and selection because of genotype-specific interactions depending on species and genotypic community composition may modulate this relationship.     

Population biology of Avena

Summary Pollen and seed dispersal patterns were analyzed in both natural and experimental populations of Avena barbata. Localized estimates of gene flow rates and plant densities gave estimates of neighborhood size in the range of 40 to 400 plants; the estimates of mean rate and distance of gene flow seemed to vary widely due to variable wind direction, rodent activity, microsite heterogeneity, etc. The relative sizes of neighborhoods in several populations were correlated with the patchy distribution of different genotypes (scored for lemma color and leaf sheath hairiness) within short distances, but patch sizes had a wide range among different sites. Highly localized gene flow patterns seemed to account for the observed pattern of highly patchy variation even when the dispersal curves for both pollen and seed were platykurtic in many cases. Measures of the stability of patches in terms of their size, dispersion in space and genetic structure in time are needed in order to sort out the relative roles of founder effects, random drift (due to small neighborhood size), and highly localized selection. However, our observations suggest that many variables and stochastic processes are involved in such studies so as to allow only weak inference about the underlying role of natural selection, drift and factors of population regulatien.     

Isozyme Variation and Genetic Diversity of the European Barley Powdery Mildew Population

Isozyme and virulence analyses of Erysiphe graminis bordei were performed with samples collected from different sites from nearly all over Europe. Isozymes and unspecific proteins extracted from conidia were separated by starch gel electrophoresis and isoelectric focusing, respectively, and the resulting isozyme banding patterns were compared with the corresponding virulence data. One isozyme phenotype prevailed in all samples. Only 7.9% of 280 isolates showed divergent banding patterns. Expected frequencies of isolates with divergent banding patterns were calculated for each subsample. In the Italian subsample, isolates with divergent banding patterns were significantly more frequent than expected. At the same time, isolates from Italy had significantly fewer virulence factors than those from N.W. Europe, indicating weaker selection by host resistance genes. It is suggested that isozyme uniformity in the homogeneous north-western European barley powdery mildew population has arisen from strong selection pressures for specific virulence genes. The direction of this selection, acting upon a mainly asexually reproducing population, has changed over space and time due to the introduction of new resistance genes, forcing local populations through bottlenecks. This may have led to random loss of genetic variation (genetic drift) in the barley powdery mildew gene pool.     

Genetic variability and drift load in populations of an aquatic snail

Abstract Population genetic theory predicts that in small populations, random genetic drift will fix and accumulate slightly deleterious mutations, resulting in reduced reproductive output. This genetic load due to random drift (i.e., drift load) can increase the extinction risk of small populations. We studied the relationship between genetic variability (indicator of past population size) and reproductive output in eight isolated, natural populations of the hermaphroditic snail Lymnaea stagnalis . In a common laboratory environment, snails from populations with the lowest genetic variability mature slower and have lower fecundity than snails from genetically more variable populations. This result suggests that past small population size has resulted in increased drift load, as predicted. The relationship between genetic variability and reproductive output is independent of the amount of nonrandom mating within populations. However, reproductive output and the current density of snails in the populations were not correlated. Instead, data from the natural populations suggest that trematode parasites may determine, at least in part, population densities of the snails.