Scale-dependence in geographic variation in a freshwater gastropod across the Palearctic

ABSTRACT

The reality of spatial clinal variation in morphological traits of freshwater pulmonate snails (Gastropoda: Pulmonata) has repeatedly been questioned or totally discounted. There is a lack of sound statistical evidence in the articles hitherto published on this subject supporting these claims. Here, by means of different analytical methods (analysis of spatial autocorrelation, linear regression analysis, canonical correlation analysis and others), we demonstrate that shell variation in the dwarf pond snail, Galba truncatula, is patterned in space throughout the northern and central Palearctic, with latitudinally-oriented clines in body size and in some shell proportions. Shell size in G. truncatula decreases with latitude and temperature, representing a special case of converse Bergmann cline. However, the temperature itself is hardly the main driver of shell size variation. It is argued that the shorter growing seasons at high latitudes may represent a better explanation for the observed trend. Shell proportions in the dwarf pond snails vary weakly at the macrogeographic scale, being spatially patterned at lower (mesogeographic) scales around 1200–1500?km. In general, spatial variation in G. truncatula shell size is decoupled from variation in shell shape, demonstrating clear scale-dependence similar to that found in different species of terrestrial (non-aquatic) pulmonate snails.     

Evaluating the contributions of change in investment and change in efficiency to age‐related declines in male and female reproduction

It is commonly observed that reproduction decreases with age, often at a different rate in males and females. This phenomenon is generally interpreted as senescence. Such reproductive declines may stem from at least two sources: a change in resource allocation and a decline in the ability to convert resources into offspring. This distinction is important because a shift in resource allocation may be favoured by selection, while reduced efficiency is purely deleterious. We propose a way to distinguish whether a decline in reproduction is purely deleterious based on estimating reproductive investment, output, and their ratio, efficiency. We apply this approach to the hermaphroditic snail Physa acuta and demonstrate that both male and female functions decline with age. The male decline largely stems from reduced investment into male activity while female decline is due to increased reproductive inefficiency. This shows that age‐related declines in reproduction can occur for a number of different reasons, a distinction that is usually masked by the general term ‘senescence’. This approach could be applied to any species to evaluate age‐related reproductive decline. We advocate that future studies measure age trajectories of reproductive investment and output to explore the potential processes hidden behind the observation that reproduction declines with age.     

Body size-dependent gender role in a simultaneous hermaphrodite freshwater snail, Physa acuta

We examined whether gender role in the simultaneous hermaphroditefreshwater snail, Physa acuta, is determined by relative bodysize in a manner predicted by the size-advantage model. We observedthe body-size combinations of pairs in the laboratory by usingfield-collected populations. Smaller individuals tended to playthe "male" role (sperm donor), and larger snails the "female"(sperm recipient). Next, we analyzed the mating behaviors involvedin gender-role decision in snail pairs of three different body-sizecombinations, using "large" and "small" snails. Smaller snailswere more likely to approach the partner as a male in different-sizecombination (large/small), whereas frequent initial approachesas a male and rejection behavior as a female were observed inthe large/large combination. Third, we examined the body sizepreference when a snail can freely choose the partner from twoother individuals of different body sizes (large/large/smallor large/small/small). Small individuals had a significant tendencyto act as the male and positively selected large snails as thefemale partner in both triple combinations. However, the largeindividual acted as both the male and the female with nearlyequal frequency. In the size-differing pairings, copulationsoccurred after fewer male approaches and fewer rejections thanin pairings involving two large snails, suggesting that bodysize difference is one of the behavioral solutions in genderconflict. Clear gender-role switching associated with body sizewas not seen. Smaller snails thus have a tendency to play themale role more frequently but adopt both gender roles when theirbody size is sufficiently large.     

Small effective population size in the long-toed salamander

The effective population sizes (Ne) of six populations of the long-toed salamander (Ambystoma macrodactylum) from Montana and Idaho, USA were estimated from allozyme data from samples collected in 1978, 1996 and 1997 using the temporal allele frequency method. Five of the six estimates ranged from 23 to 207 (mean = 123 +/- 79); one estimate was indistinguishable from infinity. In order to infer the actual Ne of salamander populations, we compared the frequency distribution of our observed Ne estimates with distributions obtained from simulated populations of known Ne. Our observed Ne estimate distribution was consistent with distributions from simulated populations with Ne values of 10, 25, and 50, suggesting an actual Ne for each of the six salamander populations of less than 100. This Ne estimate agrees with most other Ne estimates for amphibians. We conclude by discussing the conservation implications of small Ne values in amphibians in the context of increasing isolation of populations due to habitat fragmentation.     

Variance effective population size is affected by census size in sub-structured populations

Measurement of allele frequency shifts between temporally spaced samples has long been used for assessment of effective population size (N_e), and this ‘temporal method’ provides estimates of N_e referred to as variance effective size (N_eV). We show that N_eV of a local population that belongs to a sub-structured population (a metapopulation) is determined not only by genetic drift and migration rate (m), but also by the census size (N_c). The realized N_eV of a local population can either increase or decrease with increasing m, depending on the relationship between N_e and N_c in isolation. This is shown by explicit mathematical expressions for the factors affecting N_eV derived for an island model of migration. We verify analytical results using high-resolution computer simulations, and show that the phenomenon is not restricted to the island model migration pattern. The effect of N_c on the realized N_eV of a local subpopulation is most pronounced at high migration rates. We show that N_c only affects local N_eV, whereas N_eV for the metapopulation as a whole, inbreeding (N_eI), and linkage disequilibrium (N_eLD) effective size are all independent of N_c. Our results provide a possible explanation to the large variation of N_e/N_c ratios reported in the literature, where N_e is frequently estimated by N_eV. They are also important for the interpretation of empirical N_e estimates in genetic management where local N_eV is often used as a substitute for inbreeding effective size, and we suggest an increased focus on metapopulation N_eV as a proxy for N_eI.     

Effective population size as a driver for divergence of an antimicrobial peptide (Hymenoptaecin) in two common European bumblebee species

Social insects are the target of numerous pathogens. This is because the high density of closely‐related individuals frequently interacting with each other enhances the transmission and establishment of pathogens. This high selective pressure results in the rapid evolution of immune genes, which might be counteracted by a reduced effective population size (N_e) lowering the effectiveness of selection. We tested the effect of N_e on the evolutionary rate of an important immune gene for the antimicrobial peptide Hymenoptaecin in two common central European bumblebee species: Bombus terrestris and Bombus lapidarius. Both species are similar in their biology and are expected to be under similar selective pressures because pathogen prevalence does not differ between species. However, previous studies indicated a higher N_e in B. terrestris compared to B. lapidarius. We found high intraspecific variability in the coding sequence but low variability for silent polymorphisms in B. lapidarius. Estimates of long‐ and short‐term N_e were three‐ to four‐fold higher N_e in B. terrestris, although the species did not differ in census population sizes. The difference in N_e might result in less efficient selection and suboptimal adaptation of immune genes (e.g. hymenoptaecin) in B. lapidarius, and thus this species might become less resistant and more tolerant, turning into a superspreader of diseases.     

Estimating genetic drift and effective population size from temporal shifts in dominant gene marker frequencies

Measurement of temporal change in allele frequencies represents an indirect method for estimating the genetically effective size of populations. When allele frequencies are estimated for gene markers that display dominant gene expression, such as, e.g. random amplified polymorphic DNA (RAPD) and amplified fragment length polymorphism (AFLP) markers, the estimates can be seriously biased. We quantify bias for previous allele frequency estimators and present a new expression that is generally less biased and provides a more precise assessment of temporal allele frequency change. We further develop an estimator for effective population size that is appropriate when dealing with dominant gene markers. Comparison with estimates based on codominantly expressed genes, such as allozymes or microsatellites, indicates that about twice as many loci or sampled individuals are required when using dominant markers to achieve the same precision.     

Effects of population characteristics and structure on estimates of effective population size in a house sparrow metapopulation

Effective population size (N_e) is a key parameter to understand evolutionary processes and the viability of endangered populations as it determines the rate of genetic drift and inbreeding. Low N_e can lead to inbreeding depression and reduced population adaptability. In this study, we estimated contemporary N_e using genetic estimators (LDNE, ONeSAMP, MLNE and CoNe) as well as a demographic estimator in a natural insular house sparrow metapopulation. We investigated whether population characteristics (population size, sex ratio, immigration rate, variance in population size and population growth rate) explained variation within and among populations in the ratio of effective to census population size (N_e/N_c). In general, N_e/N_c ratios increased with immigration rates. Genetic N_e was much larger than demographic N_e, probably due to a greater effect of immigration on genetic than demographic processes in local populations. Moreover, although estimates of genetic N_e seemed to track N_c quite well, the genetic N_e‐estimates were often larger than N_c within populations. Estimates of genetic N_e for the metapopulation were however within the expected range (<N_c). Our results suggest that in fragmented populations, even low levels of gene flow may have important consequences for the interpretation of genetic estimates of N_e. Consequently, further studies are needed to understand how N_e estimated in local populations or the total metapopulation relates to actual rates of genetic drift and inbreeding.     

Gene genealogies when the sample size exceeds the effective size of the population

We study the properties of gene genealogies for large samples using a continuous approximation introduced by R. A. Fisher. We show that the major effect of large sample size, relative to the effective size of the population, is to increase the proportion of polymorphisms at which the mutant type is found in a single copy in the sample. We derive analytical expressions for the expected number of these singleton polymorphisms and for the total number of polymorphic, or segregating, sites that are valid even when the sample size is much greater than the effective size of the population. We use simulations to assess the accuracy of these predictions and to investigate other aspects of large-sample genealogies. Lastly, we apply our results to some data from Pacific oysters sampled from British Columbia. This illustrates that, when large samples are available, it is possible to estimate the mutation rate and the effective population size separately, in contrast to the case of small samples in which only the product of the mutation rate and the effective population size can be estimated.     

Generation time and effective population size in Polar Eskimos

North Greenland Polar Eskimos are the only hunter-gatherer population, to our knowledge, who can offer precise genealogical records spanning several generations. This is the first report from Eskimos on two key parameters in population genetics, namely, generation time (T) and effective population size (Ne). The average mother-daughter and father-son intervals were 27 and 32 years, respectively, roughly similar to the previously published generation times obtained from recent agricultural societies across the world. To gain an insight for the generation time in our distant ancestors, we calculated maternal generation time for two wild chimpanzee populations. We also provide the first comparison among three distinct approaches (genealogy, variance and life table methods) for calculating Ne, which resulted in slightly differing values for the Eskimos. The ratio of the effective to the census population size is estimated as 0.6-0.7 for autosomal and X-chromosomal DNA, 0.7-0.9 for mitochondrial DNA and 0.5 for Y-chromosomal DNA. A simulation of alleles along the genealogy suggested that Y-chromosomal DNA may drift a little faster than mitochondrial DNA in this population, in contrast to agricultural Icelanders. Our values will be useful not only in prehistoric population inference but also in understanding the shaping of our genome today.     

Inbreeding effective population size under some artificial selection schemes

Summary It is well known that truncation selection is the most efficient form of directional selection in terms of changing gene frequency. In this paper we show circumstances where truncation selection followed by a balanced mating generates inbreeding effective population size smaller than that generated by a selection that assigns mating frequencies to individuals according to their breeding values, where both selection schemes give the same expected performance of selected individuals (selection differential). Breeding values of selected individuals and the weight used to determine mating frequencies are assumed to be linearly distributed on a performance scales, x. To assign mating frequencies to the individuals in the weighting system, the selected individuals are grouped using a constant , and ith group in the interval x_i, x_i + . With small number of groups, say 2 or 3, the weighting system in general generates inbreeding effective population size that is larger than that generated by a truncation selection. As the number of the groups increases, truncation selection generates larger effective numbers.     

Sensitivity analysis of effective population size to demographic parameters in house sparrow populations

The ratio between the effective and the census population size, , is an important measure of the long‐term viability and sustainability of a population. Understanding which demographic processes that affect most will improve our understanding of how genetic drift and the probability of fixation of alleles is affected by demography. This knowledge may also be of vital importance in management of endangered populations and species. Here, we use data from 13 natural populations of house sparrow (Passer domesticus) in Norway to calculate the demographic parameters that determine . Using the global variance‐based Sobol’ method for the sensitivity analyses, we found that was most sensitive to demographic variance, especially among older individuals. Furthermore, the individual reproductive values (that determine the demographic variance) were most sensitive to variation in fecundity. Our results draw attention to the applicability of sensitivity analyses in population management and conservation. For population management aiming to reduce the loss of genetic variation, a sensitivity analysis may indicate the demographic parameters towards which resources should be focused. The result of such an analysis may depend on the life history and mating system of the population or species under consideration, because the vital rates and sex–age classes that is most sensitive to may change accordingly.     

Estimating effective population size from linkage disequilibrium: severe bias in small samples

Effective population size (N _e) is a central concept in evolutionary biology and conservation genetics. It predicts rates of loss of neutral genetic variation, fixation of deleterious and favourable alleles, and the increase of inbreeding experienced by a population. A method exists for the estimation of N _e from the observed linkage disequilibrium between unlinked loci in a population sample. While an increasing number of studies have applied this method in natural and managed populations, its reliability has not yet been evaluated. We developed a computer program to calculate this estimator of N _e using the most widely used linkage disequilibrium algorithm and used simulations to show that this estimator is strongly biased when the sample size is small (<‰100) and below the true N _e. This is probably due to the linkage disequilibrium generated by the sampling process itself and the inadequate correction for this phenomenon in the method. Results suggest that N _e estimates derived using this method should be regarded with caution in many cases. To improve the method’s reliability and usefulness we propose a way to determine whether a given sample size exceeds the population N _e and can therefore be used for the computation of an unbiased estimate.     

Estimation of effective population size and detection of a recent population decline coinciding with habitat fragmentation in a ground beetle

We assess the impact of habitat fragmentation on the effective size (N(e)) of local populations of the flightless ground beetle Carabus violaceus in a small (<25 ha) and a large (>80 ha) forest fragment separated by a highway. N(e) was estimated based on the temporal variation of allele frequencies at 13 microsatellite loci using two different methods. In the smaller fragment, N(e) estimates ranged between 59 and a few hundred, whereas values between 190 and positive infinity were estimated for the larger fragment. In both samples, we detected a signal of population decline, which was stronger in the small fragment. The estimated time of onset of this N(e) reduction was consistent with the hypothesis that recent road constructions have divided a continuous population into several isolated subpopulations. In the small fragment, N(e) of the local population may be so small that its long-term persistence is endangered.     

An empirical test of the estimation of historical effective population size using Drosophila melanogaster

The availability of a large number of high-density markers (SNPs) allows the estimation of historical effective population size (N_e) from linkage disequilibrium between loci. A recent refinement of methods to estimate historical N_e from the recent past has been shown to be rather accurate with simulation data. The method has also been applied to real data for numerous species. However, the simulation data cannot encompass all the complexities of real genomes, and the performance of any estimation method with real data is always uncertain, as the true demography of the populations is not known. Here, we carried out an experimental design with Drosophila melanogaster to test the method with real data following a known demographic history. We used a population maintained in the laboratory with a constant census size of about 2800 individuals and subjected the population to a drastic decline to a size of 100 individuals. After a few generations, the population was expanded back to the previous size and after a few further generations again expanded to twice the initial size. Estimates of historical N_e were obtained with the software GONE both for autosomal and X chromosomes from samples of 17 individuals sequenced for the whole genome. Estimates of the historical effective size were able to infer the patterns of changes that occurred in the populations showing generally good performance of the method. We discuss the limitations of the method and the application of the software carried out so far.     

Contemporary effective population and metapopulation size (Ne and meta‐Ne): comparison among three salmonids inhabiting a fragmented system and differing in gene flow and its asymmetries

We estimated local and metapopulation effective sizes ( and meta‐) for three coexisting salmonid species (Salmo salar, Salvelinus fontinalis, Salvelinus alpinus) inhabiting a freshwater system comprising seven interconnected lakes. First, we hypothesized that might be inversely related to within‐species population divergence as reported in an earlier study (i.e., F_ST: S. salar> S. fontinalis> S. alpinus). Using the approximate Bayesian computation method implemented in ONeSAMP, we found significant differences in () between species, consistent with a hierarchy of adult population sizes (). Using another method based on a measure of linkage disequilibrium (LDNE: ), we found more finite values for S. salar than for the other two salmonids, in line with the results above that indicate that S. salar exhibits the lowest among the three species. Considering subpopulations as open to migration (i.e., removing putative immigrants) led to only marginal and non‐significant changes in , suggesting that migration may be at equilibrium between genetically similar sources. Second, we hypothesized that meta‐ might be significantly smaller than the sum of local s (null model) if gene flow is asymmetric, varies among subpopulations, and is driven by common landscape features such as waterfalls. One ‘bottom‐up’ or numerical approach that explicitly incorporates variable and asymmetric migration rates showed this very pattern, while a number of analytical models provided meta‐ estimates that were not significantly different from the null model or from each other. Our study of three species inhabiting a shared environment highlights the importance and utility of differentiating species‐specific and landscape effects, not only on dispersal but also in the demography of wild populations as assessed through local s and meta‐s and their relevance in ecology, evolution and conservation.     

Monitoring the effective population size of a brown bear (Ursus arctos) population using new single-sample approaches

The effective population size (N(e) ) could be the ideal parameter for monitoring populations of conservation concern as it conveniently summarizes both the evolutionary potential of the population and its sensitivity to genetic stochasticity. However, tracing its change through time is difficult in natural populations. We applied four new methods for estimating N(e) from a single sample of genotypes to trace temporal change in N(e) for bears in the Northern Dinaric Mountains. We genotyped 510 bears using 20 microsatellite loci and determined their age. The samples were organized into cohorts with regard to the year when the animals were born and yearly samples with age categories for every year when they were alive. We used the Estimator by Parentage Assignment (EPA) to directly estimate both N(e) and generation interval for each yearly sample. For cohorts, we estimated the effective number of breeders (N(b) ) using linkage disequilibrium, sibship assignment and approximate Bayesian computation methods and extrapolated these estimates to N(e) using the generation interval. The N(e) estimate by EPA is 276 (183-350 95% CI), meeting the inbreeding-avoidance criterion of N(e) > 50 but short of the long-term minimum viable population goal of N(e) > 500. The results obtained by the other methods are highly consistent with this result, and all indicate a rapid increase in N(e) probably in the late 1990s and early 2000s. The new single-sample approaches to the estimation of N(e) provide efficient means for including N(e) in monitoring frameworks and will be of great importance for future management and conservation.     

Census vs. effective population size in chinook salmon: large- and small-scale environmental perturbation effects

Population viability has often been assessed by census of reproducing adults. Recently this method has been called into question and estimation of the effective population size (Ne) proposed as a complementary method to determine population health. We examined genetic diversity in five populations of chinook salmon (Oncorhynchus tshawytscha) from the upper Fraser River watershed (British Columbia, Canada) at 11 microsatellite loci over 20 years using DNA extracted from archived scale samples. We tested for changes in genetic diversity, calculated the ratio of the number of alleles to the range in allele size to give the statistic M, calculated Ne from the temporal change in allele frequency, used the maximum likelihood method to calculate effective population size (NeM), calculated the harmonic mean of population size, and compared these statistics to annual census estimates. Over the last two decades population size has increased in all five populations of chinook examined; however, Ne calculated for each population was low (81-691) and decreasing over the time interval measured. Values of NeM were low, but substantially higher than Ne calculated using the temporal method. The calculated values for M were generally low (M < 0.70), indicating recent population reductions for all five populations. Large-scale historic barriers to migration and development activities do not appear to account for the low values of Ne; however, available spawning area is positively correlated with Ne. Both Ne and M estimates indicate that these populations are potentially susceptible to inbreeding effects and may lack the ability to respond adaptively to stochastic events. Our findings question the practice of relying exclusively on census estimates for interpreting population health and show the importance of determining genetic diversity within populations.     

Post‐fragmentation population structure in a cooperative breeding Afrotropical cloud forest bird: emergence of a source‐sink population network

The impact of demographic parameters on the genetic population structure and viability of organisms is a long‐standing issue in the study of fragmented populations. Demographic and genetic tools are now readily available to estimate census and effective population sizes and migration and gene flow rates with increasing precision. Here we analysed the demography and genetic population structure over a recent 15‐year time span in five remnant populations of Cabanis's greenbul (Phyllastrephus cabanisi), a cooperative breeding bird in a severely fragmented cloud forest habitat. Contrary to our expectation, genetic admixture and effective population sizes slightly increased, rather than decreased between our two sampling periods. In spite of small effective population sizes in tiny forest remnants, none of the populations showed evidence of a recent population bottleneck. Approximate Bayesian modelling, however, suggested that differentiation of the populations coincided at least partially with an episode of habitat fragmentation. The ratio of meta‐N_e to meta‐N_c was relatively low for birds, which is expected for cooperative breeding species, while N_e/N_c ratios strongly varied among local populations. While the overall trend of increasing population sizes and genetic admixture may suggest that Cabanis's greenbuls increasingly cope with fragmentation, the time period over which these trends were documented is rather short relative to the average longevity of tropical species. Furthermore, the critically low N_c in the small forest remnants keep the species prone to demographic and environmental stochasticity, and it remains open if, and to what extent, its cooperative breeding behaviour helps to buffer such effects.