A genetic demographic analysis of Lake Malawi rock‐dwelling cichlids using spatio‐temporal sampling

We estimated the effective population sizes (N_e) and tested for short‐term temporal demographic stability of populations of two Lake Malawi cichlids: Maylandia benetos, a micro‐endemic, and Maylandia zebra, a widespread species found across the lake. We sampled a total of 351 individuals, genotyped them at 13 microsatellite loci and sequenced their mitochondrial D‐loop to estimate genetic diversity, population structure, demographic history and effective population sizes. At the microsatellite loci, genetic diversity was high in all populations. Yet, genetic diversity was relatively low for the sequence data. Microsatellites yielded mean N_e estimates of 481 individuals (±99 SD) for M. benetos and between 597 (±106.3 SD) and 1524 (±483.9 SD) individuals for local populations of M. zebra. The microsatellite data indicated no deviations from mutation–drift equilibrium. Maylandia zebra was further found to be in migration–drift equilibrium. Temporal fluctuations in allele frequencies were limited across the sampling period for both species. Bayesian Skyline analyses suggested a recent expansion of M. zebra populations in line with lake‐level fluctuations, whereas the demographic history of M. benetos could only be estimated for the very recent past. Divergence time estimates placed the origin of M. benetos within the last 100 ka after the refilling of the lake and suggested that it split off the sympatric M. zebra population. Overall, our data indicate that micro‐endemics and populations in less favourable habitats have smaller N_e, indicating that drift may play an important role driving their divergence. Yet, despite small population sizes, high genetic variation can be maintained.     

Decomposing demographic contributions to the effective population size with moose as a case study

Levels of random genetic drift are influenced by demographic factors, such as mating system, sex ratio and age structure. The effective population size (N_e) is a useful measure for quantifying genetic drift. Evaluating relative contributions of different demographic factors to N_e is therefore important to identify what makes a population vulnerable to loss of genetic variation. Until recently, models for estimating N_e have required many simplifying assumptions, making them unsuitable for this task. Here, using data from a small, harvested moose population, we demonstrate the use of a stochastic demographic framework allowing for fluctuations in both population size and age distribution to estimate and decompose the total demographic variance and hence the ratio of effective to total population size (N_e/N) into components originating from sex, age, survival and reproduction. We not only show which components contribute most to N_e/N currently, but also which components have the greatest potential for changing N_e/N. In this relatively long‐lived polygynous system we show that N_e/N is most sensitive to the demographic variance of older males, and that both reproductive autocorrelations (i.e., a tendency for the same individuals to be successful several years in a row) and covariance between survival and reproduction contribute to decreasing N_e/N (increasing genetic drift). These conditions are common in nature and can be caused by common hunting strategies. Thus, the framework presented here has great potential to increase our understanding of the demographic processes that contribute to genetic drift and viability of populations, and to inform management decisions.     

EFFECTIVE POPULATION SIZE AND GENETIC DRIFT IN TRISTYLOUS EICHHORNIA PANICULATA (PONTEDERIACEAE)

Populations of the tristylous, annual Eichhornia paniculata are markedly differentiated with respect to frequency of mating types. This variation is associated with evolutionary changes in mating system, from predominant outcrossing to high self-fertilization. To assess the potential influence of genetic drift acting on this variation, we estimated effective population size in 10 populations from northeastern Brazil using genetic and demographic methods. Effective size (N_e) was inferred from temporal changes in allele frequency at two to eight isozyme loci and also calculated using five demographic variables: 1) the number of flowering individuals (N); 2) temporal fluctuations in N; 3) variance in flower number; 4) frequency of mating types; and 5) selfing rate. Average N_e based on isozyme data was 15.8, range 3.4–70.6, and represented a fraction (mean N_e/N = 0.106) of the census number of individuals (mean N = 762.8; range: 30.5–5,040). Temporal variation in N and variance in flower number each reduced N_e to about a half of N whereas mating type frequencies and selfing rate caused only small reductions in N_e relative to N. All estimates of N_e based on demographic variables were considerably larger than those obtained from genetic data. The two kinds of estimates were in general agreement, however, when all demographic variables were combined into a single measure. Monte Carlo simulations indicated that effective size must be fewer than about 40 for drift to overcome the frequency-dependent selection that maintains the polymorphism for mating type. Applying the average N_e/N value to 167 populations censused in northeastern Brazil indicated that 72% had effective sizes below this number. This suggests that genetic drift is likely to play a dominant role in natural populations of E. paniculata.     

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.     

Consequences of multiple mating‐system shifts for population and range‐wide genetic structure in a coastal dune plant

Evolutionary transitions from outcrossing to selfing can strongly affect the genetic diversity and structure of species at multiple spatial scales. We investigated the genetic consequences of mating‐system shifts in the North American, Pacific coast dune endemic plant Camissoniopsis cheiranthifolia (Onagraceae) by assaying variation at 13 nuclear (n) and six chloroplast (cp) microsatellite (SSR) loci for 38 populations across the species range. As predicted from the expected reduction in effective population size (N_e) caused by selfing, small‐flowered, predominantly selfing (SF) populations had much lower nSSR diversity (but not cpSSR) than large‐flowered, predominantly outcrossing (LF) populations. The reduction in nSSR diversity was greater than expected from the effects of selfing on N_e alone, but could not be accounted for by indirect effects of selfing on population density. Although selfing should reduce gene flow, SF populations were not more genetically differentiated than LF populations. We detected five clusters of nSSR genotypes and three groups of cpSSR haplotypes across the species range consisting of parapatric groups of populations that usually (but not always) differed in mating system, suggesting that selfing may often initiate ecogeographic isolation. However, lineage‐wide genetic variation was not lower for selfing clusters, failing to support the hypothesis that selection for reproductive assurance spurred the evolution of selfing in this species. Within three populations where LF and SF plants coexist, we detected genetic differentiation among diverged floral phenotypes suggesting that reproductive isolation (probably postzygotic) may help maintain the striking mating‐system differentiation observed across the range of this species.     

Unexpected high genetic diversity in small populations suggests maintenance by associative overdominance

The effective population size (N_e) is a central factor in determining maintenance of genetic variation. The neutral theory predicts that loss of variation depends on N_e, with less genetic drift in larger populations. We monitored genetic drift in 42 Drosophila melanogaster populations of different adult census population sizes (10, 50 or 500) using pooled RAD sequencing. In small populations, variation was lost at a substantially lower rate than expected. This observation was consistent across two ecological relevant thermal regimes, one stable and one with a stressful increase in temperature across generations. Estimated ratios between N_e and adult census size were consistently higher in small than in larger populations. The finding provides evidence for a slower than expected loss of genetic diversity and consequently a higher than expected long‐term evolutionary potential in small fragmented populations. More genetic diversity was retained in areas of low recombination, suggesting that associative overdominance, driven by disfavoured homozygosity of recessive deleterious alleles, is responsible for the maintenance of genetic diversity in smaller populations. Consistent with this hypothesis, the X‐chromosome, which is largely free of recessive deleterious alleles due to hemizygosity in males, fits neutral expectations even in small populations. Our experiments provide experimental answers to a range of unexpected patterns in natural populations, ranging from variable diversity on X‐chromosomes and autosomes to surprisingly high levels of nucleotide diversity in small populations.     

Effective size of density‐dependent populations in fluctuating environments

Reliable estimates of effective population size are of central importance in population genetics and evolutionary biology. For populations that fluctuate in size, harmonic mean population size is commonly used as a proxy for (multi‐) generational effective size. This assumes no effects of density dependence on the ratio between effective and actual population size, which limits its potential application. Here, we introduce density dependence on vital rates in a demographic model of variance effective size. We derive an expression for the ratio in a density‐regulated population in a fluctuating environment. We show by simulations that yearly genetic drift is accurately predicted by our model, and not proportional to as assumed by the harmonic mean model, where N is the total population size of mature individuals. We find a negative relationship between and N. For a given N, the ratio depends on variance in reproductive success and the degree of resource limitation acting on the population growth rate. Finally, our model indicate that environmental stochasticity may affect not only through fluctuations in N, but also for a given N at a given time. Our results show that estimates of effective population size must include effects of density dependence and environmental stochasticity.     

Life‐history evolution under fluctuating density‐dependent selection and the adaptive alignment of pace‐of‐life syndromes

We present a novel perspective on life‐history evolution that combines recent theoretical advances in fluctuating density‐dependent selection with the notion of pace‐of‐life syndromes (POLSs) in behavioural ecology. These ideas posit phenotypic co‐variation in life‐history, physiological, morphological and behavioural traits as a continuum from the highly fecund, short‐lived, bold, aggressive and highly dispersive ‘fast’ types at one end of the POLS to the less fecund, long‐lived, cautious, shy, plastic and socially responsive ‘slow’ types at the other. We propose that such variation in life histories and the associated individual differences in behaviour can be explained through their eco‐evolutionary dynamics with population density – a single and ubiquitous selective factor that is present in all biological systems. Contrasting regimes of environmental stochasticity are expected to affect population density in time and space and create differing patterns of fluctuating density‐dependent selection, which generates variation in fast versus slow life histories within and among populations. We therefore predict that a major axis of phenotypic co‐variation in life‐history, physiological, morphological and behavioural traits (i.e. the POLS) should align with these stochastic fluctuations in the multivariate fitness landscape created by variation in density‐dependent selection. Phenotypic plasticity and/or genetic (co‐)variation oriented along this major POLS axis are thus expected to facilitate rapid and adaptively integrated changes in various aspects of life histories within and among populations and/or species. The fluctuating density‐dependent selection POLS framework presented here therefore provides a series of clear testable predictions, the investigation of which should further our fundamental understanding of life‐history evolution and thus our ability to predict natural population dynamics.     

DOES MATE LIMITATION IN SELF‐INCOMPATIBLE SPECIES PROMOTE THE EVOLUTION OF SELFING? THE CASE OF LEAVENWORTHIA ALABAMICA

Genetic diversity at the S‐locus controlling self‐incompatibility (SI) is often high because of negative frequency‐dependent selection. In species with highly patchy spatial distributions, genetic drift can overwhelm balancing selection and cause stochastic loss of S‐alleles. Natural selection may favor the breakdown of SI in populations with few S‐alleles because low S‐allele diversity constrains the seed production of self‐incompatible plants. We estimated S‐allele diversity, effective population sizes, and migration rates in Leavenworthia alabamica, a self‐incompatible mustard species restricted to discrete habitat patches in rocky glades. Patterns of polymorphism were investigated at the S‐locus and 15 neutral microsatellites in three large and three small populations with 100‐fold variation in glade size. Populations on larger glades maintained more S‐alleles, but all populations were estimated to harbor at least 20 S‐alleles, and mate availabilities typically exceeded 0.80, which is consistent with little mate limitation in nature. Estimates of the effective size (N_e) in each population ranged from 600 to 1600, and estimated rates of migration (m) ranged from 3 × 10⁻⁴ to nearly 1 × 10⁻³. According to theoretical models, there is limited opportunity for genetic drift to reduce S‐allele diversity in populations with these attributes. Although pollinators or resources limit seed production in small glades, limited S‐allele diversity does not appear to be a factor promoting the incipient breakdown of SI in populations of this species that were studied.     

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.     

Trapped within the city: integrating demography,time since isolation and population‐specific traits to assess the genetic effects of urbanization

Urbanization is a severe form of habitat fragmentation that can cause many species to be locally extirpated and many others to become trapped and isolated within an urban matrix. The role of drift in reducing genetic diversity and increasing genetic differentiation is well recognized in urban populations. However, explicit incorporation and analysis of the demographic and temporal factors promoting drift in urban environments are poorly studied. Here, we genotyped 15 microsatellites in 320 fire salamanders from the historical city of Oviedo (Est. 8th century) to assess the effects of time since isolation, demographic history (historical effective population size; N_e) and patch size on genetic diversity, population structure and contemporary N_e. Our results indicate that urban populations of fire salamanders are highly differentiated, most likely due to the recent N_e declines, as calculated in coalescence analyses, concomitant with the urban development of Oviedo. However, urbanization only caused a small loss of genetic diversity. Regression modelling showed that patch size was positively associated with contemporary N_e, while we found only moderate support for the effects of demographic history when excluding populations with unresolved history. This highlights the interplay between different factors in determining current genetic diversity and structure. Overall, the results of our study on urban populations of fire salamanders provide some of the very first insights into the mechanisms affecting changes in genetic diversity and population differentiation via drift in urban environments, a crucial subject in a world where increasing urbanization is forecasted.     

Variance in Female Reproductive Success Differentially Impacts Effective Population Size in the Short-Nosed Fruit Bat, <Emphasis Type="Italic">Cynopterus sphinx</Emphasis>

Effective population size (N _e) quantifies the effects of micro-evolutionary processes and the rate of loss of genetic diversity in a population. Several demographic and mating parameters reduce N _e. Theoretical studies elucidate the impacts of various demographic and mating system parameters on N _e, while empirical studies illustrate realized N _e for species with differing life histories and mating systems. However, effect of intra-specific variation in mating system on effective size remains largely unexplored. In this paper we investigated the effect of promiscuous and polygynous mating on N _e in two wild populations of the short-nosed fruit bat, Cynopterus sphinx. N _e/N (ratio of effective population size to census size) was lower than unity in both populations, and much lower for the polygynous population compared to promiscuous population. Elasticity analyses reveal that N _e/N was sensitive to deviations in the sex ratio. Variance in female reproductive success had a higher impact on N _e compared to variance in male reproductive success in the promiscuous population. However, for the polygynous population, impact of variance in male reproductive success on N _e was higher than that of variance in female reproductive success. Our results suggest that depending on mating system, different populations of the same species could have alternate evolutionary trajectories. The rate of loss of genetic diversity would be lower for the promiscuous population compared to the polygynous population. Our study is the first to highlight which parameters would most significantly impact population specific N _e under different mating systems.     

Coupling Demographic and Genetic Variability from Archived Collections of European Anchovy (Engraulis encrasicolus)

It is well known that temporal fluctuations in small populations deeply influence evolutionary potential. Less well known is whether fluctuations can influence the evolutionary potentials of species with large census sizes. Here, we estimated genetic population parameters from as survey of polymorphic microsatellite DNA loci in archived otoliths from Adriatic European anchovy (Engraulis encrasicolus), a fish with large census sizes that supports numerous local fisheries. Stocks have fluctuated greatly over the past few decades, and the Adriatic fishery collapsed in 1987. Our results show a significant reduction of mean genetic parameters as a consequence of the population collapse. In addition, estimates of effective population size (N_e) are much smaller than those expected in a fishes with large population census sizes (N_c). Estimates of N_e indicate low effective population sizes, even before the population collapse. The ratio N_e/N_e ranged between 10⁻⁶ and 10⁻⁸, indicating a large discrepancy between the anchovy gene pool and population census size. Therefore, anchovy populations may be more vulnerable to fishery effort and environmental change than previously thought.     

Adaptive divergence despite strong genetic drift: genomic analysis of the evolutionary mechanisms causing genetic differentiation in the island fox (Urocyon littoralis)

The evolutionary mechanisms generating the tremendous biodiversity of islands have long fascinated evolutionary biologists. Genetic drift and divergent selection are predicted to be strong on islands and both could drive population divergence and speciation. Alternatively, strong genetic drift may preclude adaptation. We conducted a genomic analysis to test the roles of genetic drift and divergent selection in causing genetic differentiation among populations of the island fox (Urocyon littoralis). This species consists of six subspecies, each of which occupies a different California Channel Island. Analysis of 5293 SNP loci generated using Restriction‐site Associated DNA (RAD) sequencing found support for genetic drift as the dominant evolutionary mechanism driving population divergence among island fox populations. In particular, populations had exceptionally low genetic variation, small N_e (range = 2.1–89.7; median = 19.4), and significant genetic signatures of bottlenecks. Moreover, islands with the lowest genetic variation (and, by inference, the strongest historical genetic drift) were most genetically differentiated from mainland grey foxes, and vice versa, indicating genetic drift drives genome‐wide divergence. Nonetheless, outlier tests identified 3.6–6.6% of loci as high F_ST outliers, suggesting that despite strong genetic drift, divergent selection contributes to population divergence. Patterns of similarity among populations based on high F_ST outliers mirrored patterns based on morphology, providing additional evidence that outliers reflect adaptive divergence. Extremely low genetic variation and small N_e in some island fox populations, particularly on San Nicolas Island, suggest that they may be vulnerable to fixation of deleterious alleles, decreased fitness and reduced adaptive potential.     

Inbreeding depression and drift load in small populations at demographic disequilibrium

Inbreeding depression is a major driver of mating system evolution and has critical implications for population viability. Theoretical and empirical attention has been paid to predicting how inbreeding depression varies with population size. Lower inbreeding depression is predicted in small populations at equilibrium, primarily due to higher inbreeding rates facilitating purging and/or fixation of deleterious alleles (drift load), but predictions at demographic and genetic disequilibrium are less clear. In this study, we experimentally evaluate how lifetime inbreeding depression and drift load, estimated by heterosis, vary with census (N_c) and effective (estimated as genetic diversity, H_e) population size across six populations of the biennial Sabatia angularis as well as present novel models of inbreeding depression and heterosis under varying demographic scenarios at disequilibrium (fragmentation, bottlenecks, disturbances). Our experimental study reveals high average inbreeding depression and heterosis across populations. Across our small sample, heterosis declined with H_e, as predicted, whereas inbreeding depression did not vary with H_e and actually decreased with N_c. Our theoretical results demonstrate that inbreeding depression and heterosis levels can vary widely across populations at disequilibrium despite similar H_e and highlight that joint demographic and genetic dynamics are key to predicting patterns of genetic load in nonequilibrium systems.     

Genomewide patterns of variation in genetic diversity are shared among populations,species and higher‐order taxa

Genomewide screens of genetic variation within and between populations can reveal signatures of selection implicated in adaptation and speciation. Genomic regions with low genetic diversity and elevated differentiation reflective of locally reduced effective population sizes (N_e) are candidates for barrier loci contributing to population divergence. Yet, such candidate genomic regions need not arise as a result of selection promoting adaptation or advancing reproductive isolation. Linked selection unrelated to lineage‐specific adaptation or population divergence can generate comparable signatures. It is challenging to distinguish between these processes, particularly when diverging populations share ancestral genetic variation. In this study, we took a comparative approach using population assemblages from distant clades assessing genomic parallelism of variation in N_e. Utilizing population‐level polymorphism data from 444 resequenced genomes of three avian clades spanning 50 million years of evolution, we tested whether population genetic summary statistics reflecting genomewide variation in N_e would covary among populations within clades, and importantly, also among clades where lineage sorting has been completed. All statistics including population‐scaled recombination rate (ρ), nucleotide diversity (π) and measures of genetic differentiation between populations (F_ST, PBS, d_xy) were significantly correlated across all phylogenetic distances. Moreover, genomic regions with elevated levels of genetic differentiation were associated with inferred pericentromeric and subtelomeric regions. The phylogenetic stability of diversity landscapes and stable association with genomic features support a role of linked selection not necessarily associated with adaptation and speciation in shaping patterns of genomewide heterogeneity in genetic diversity.     

PSMC analysis of effective population sizes in molecular ecology and its application to black‐and‐white Ficedula flycatchers

Climatic fluctuations during the Quaternary period governed the demography of species and contributed to population differentiation and ultimately speciation. Studies of these past processes have previously been hindered by a lack of means and genetic data to model changes in effective population size (N_e) through time. However, based on diploid genome sequences of high quality, the recently developed pairwise sequentially Markovian coalescent (PSMC) can estimate trajectories of changes in N_e over considerable time periods. We applied this approach to resequencing data from nearly 200 genomes of four species and several populations of the Ficedula species complex of black‐and‐white flycatchers. N_e curves of Atlas, collared, pied and semicollared flycatcher converged 1–2 million years ago (Ma) at an N_e of ≈ 200 000, likely reflecting the time when all four species last shared a common ancestor. Subsequent separate N_e trajectories are consistent with lineage splitting and speciation. All species showed evidence of population growth up until 100–200 thousand years ago (kya), followed by decline and then start of a new phase of population expansion. However, timing and amplitude of changes in N_e differed among species, and for pied flycatcher, the temporal dynamics of N_e differed between Spanish birds and central/northern European populations. This cautions against extrapolation of demographic inference between lineages and calls for adequate sampling to provide representative pictures of the coalescence process in different species or populations. We also empirically evaluate criteria for proper inference of demographic histories using PSMC and arrive at recommendations of using sequencing data with a mean genome coverage of ≥18X, a per‐site filter of ≥10 reads and no more than 25% of missing data.     

Mature male parr contribution to the effective size of an anadromous Atlantic salmon (Salmo salar) population over 30 years

We describe temporal changes in the genetic composition of a small anadromous Atlantic salmon (Salmo salar) population from South Newfoundland, an area where salmon populations are considered threatened (COSEWIC 2010). We examined the genetic variability (13 microsatellite loci) in 869 out‐migrating smolt and post‐spawning kelt samples, collected from 1985 to 2011 for a total of 22 annual collections and a 30 year span of assigned cohorts. We estimated the annual effective number of breeders (N_b) and the generational effective population size (N_e) through genetic methods and demographically using the adult sex ratio. Comparisons between genetic and demographic estimates show that the adult spawners inadequately explain the observed N_e estimates, suggesting that mature male parr are significantly increasing N_b and N_e over the study period. Spawning as parr appears to be a viable and important strategy in the near absence of adult males.     

Female polyandry and size‐assortative mating in isolated local populations of the Japanese common toad Bufo japonicus

In anurans, female polyandry under male harassment is distributed across taxa because of external aquatic fertilization. According to the sexual selection theory, male–male competition for access to females is affected by the operational sex ratio (OSR) and population density. The Japanese common toad, Bufo japonicus, is widespread in mainland Japan, and like the European common toad, B. bufo, it engages in explosive breeding. In this study, we observed the breeding behaviour of B. japonicus in isolated local populations for over four years in two breeding ponds with different population sizes and densities: large‐low (L‐pond) and small‐high (S‐pond). We analysed the relative polyandry ratio in egg clutches laid by females and estimated the size‐assortative mating pattern to be an indicator of male–male competition in the two ponds. Both ponds tended to exhibit a size‐assortative mating pattern; however, the frequency of polyandry was different in the two ponds (L‐pond = 20% and S‐pond = 90%). Our results showed that polyandry could occur without multiple amplexus with a high population density, i.e. eggs were often fertilized by free‐swimming sperm in the small shallow pond. We propose that high female polyandry ratios without continuous male harassment are generated because of a male‐biased OSR and a high population density in the small pond. © 2014 The Linnean Society of London, Biological Journal of the Linnean Society, 2014, 113 , 236–242.     

Naturally rare versus newly rare: demographic inferences on two timescales inform conservation of Galápagos giant tortoises

Long‐term population history can influence the genetic effects of recent bottlenecks. Therefore, for threatened or endangered species, an understanding of the past is relevant when formulating conservation strategies. Levels of variation at neutral markers have been useful for estimating local effective population sizes (N_e) and inferring whether population sizes increased or decreased over time. Furthermore, analyses of genotypic, allelic frequency, and phylogenetic information can potentially be used to separate historical from recent demographic changes. For 15 populations of Galápagos giant tortoises (Chelonoidis sp.), we used 12 microsatellite loci and DNA sequences from the mitochondrial control region and a nuclear intron, to reconstruct demographic history on shallow (past ~100 generations, ~2500 years) and deep (pre‐Holocene, >10 thousand years ago) timescales. At the deep timescale, three populations showed strong signals of growth, but with different magnitudes and timing, indicating different underlying causes. Furthermore, estimated historical N_e of populations across the archipelago showed no correlation with island age or size, underscoring the complexity of predicting demographic history a priori. At the shallow timescale, all populations carried some signature of a genetic bottleneck, and for 12 populations, point estimates of contemporary N_e were very small (i.e., < 50). On the basis of the comparison of these genetic estimates with published census size data, N_e generally represented ~0.16 of the census size. However, the variance in this ratio across populations was considerable. Overall, our data suggest that idiosyncratic and geographically localized forces shaped the demographic history of tortoise populations. Furthermore, from a conservation perspective, the separation of demographic events occurring on shallow versus deep timescales permits the identification of naturally rare versus newly rare populations; this distinction should facilitate prioritization of management action.