The effects of subdivision on the genetic divergence of populations and species

Abstract. An island model of migration is used to study the effects of subdivision within populations and species on sample genealogies and on between-population or between-species measures of genetic variation. The model assumes that the number of demes within each population or species is large. When populations (or species), connected either by gene flow or historical association, are themselves subdivided into demes, changes in the migration rate among demes alter both the structure of genealogies and the time scale of the coalescent process. The time scale of the coalescent is related to the effective size of the population, which depends on the migration rate among demes. When the migration rate among demes within populations is low, isolation (or speciation) events seem more recent and migration rates among populations seem higher because the effective size of each population is increased. This affects the probability of reciprocal monophyly of two samples, the chance that a gene tree of a sample matches the species tree, and relative likelihoods of different types of polymorphic sites. It can also have a profound effect on the estimation of divergence times.     

TESTING MODELS OF MIGRATION AND ISOLATION AMONG POPULATIONS OF CHINOOK SALMON (ONCORHYNCHUS TSCHAWYTSCHA)

The chinook salmon (Oncorhynchus tschawytscha) is a behaviorally, morphologically, and ecologically variable species distributed over a large geographic range. Although previous genetic surveys have revealed considerable genetic differences among populations with different life history types and from different major river drainages, it is not clear to what degree these genetically distinct populations are connected by low levels of gene flow. The work described in this paper addresses this question by surveying DNA restriction site variation at six nuclear genes from nine populations encompassing most of the species's North American range, and then attempting to fit the patterns of variation observed at these genes to five different evolutionary models using computer simulations of the coalescent process. Two commonly used constant population size models, one hypothesizing no migration among populations and one hypothesizing equal rates of migration among populations, were found to be statistically inconsistent with the observed patterns of variation. The other three models, which involved either recent divergence among populations coupled with large changes in populations size, unequal migration rates among populations, or selection, were all found to be consistent with the observed patterns of variation. Assuming selective neutrality, these results suggest that either the populations have all descended from a common ancestral population within the last ~50,000 years and have all suffered large declines in effective population size since that time, or that they have a more ancient divergence time but are connected by low levels of gene flow. These conclusions rest on the assumption of selective neutrality. With the methods employed, it was not possible to simultaneously test hypotheses of both selective neutrality and population structure.     

ECOLOGICAL AND HISTORICAL ASSOCIATIONS OF GENE FLOW IN DARTERS (TELEOSTEI: PERCIDAE)

Life history should relate to gene flow (Nm) through its influence on dispersal and effective population size. Comparative studies designed to elucidate this relationship must contend with historical events that can yield misleading estimates of gene flow and statistical problems associated with inclusion of life-history traits correlated with phylogeny. We studied the relationships of life-history characters and gene flow in 15 species of darters, a monophyletic group of stream fishes. Populations of coexisting species were sampled in three geographic regions with different Pleistocene glaciation histories. Gene flow was estimated indirectly from allozymes using two methods, 8 and private alleles. Isolation-by-distance was also tested using regression of pairwise estimates of gene flow (M?) on distance. Theta and private-alleles methods produced congruent estimates of Nm, except in a study region hypothesized to have been historically fragmented and then united following Pleistocene glaciation. A relatively weak association between life-history traits and Nm (based on θ) was observed when species from the historically fragmented region were included in stepwise regression analysis, because Nm was low despite life-history differences among taxa in this region. Excluding observations from this region produced stronger associations between clutch size and Nm (r² = 0.57), and between female size, egg size, and Nm (r² = 0.95). Additional analyses that corrected for female body size and phylogenetic nonindependence agreed that darters with high fecundity and small eggs exhibited high gene flow, whereas darters with small clutches and large eggs had low gene flow. The latter combination of life-history traits primarily is exhibited in species from headwater habitats where parental investment presumably confers survivorship on offspring. Reduced gene flow and genetic divergence among demes appear to be evolutionary consequences of this strategy.     

INTRAPOPULATIONAL VARIATION IN POLLEN-MEDIATED GENE FLOW IN PLANTAGO LANCEOLATA L.

The extent of genetic variation in pollen-mediated gene flow distance was measured for a population of the wind-pollinated herb, Plantago lanceolata. Twenty-nine genotypes were collected from a natural population, clonally replicated, and grown to reproductive maturity in a greenhouse. Relative gene flow distances were measured for each replicate and genotype in a wind tunnel. Approximately five percent of the total variation in gene flow distance was attributable to variation among genotypes. Most of the remaining variation was attributable to differences between the early part of the growing season, when most of the flowering occurs, and the remainder of the season, when flowering is sparse. The rankings of the genotypes' gene flow distances showed significant concordance between the early-season measures and measures from later in the summer. There was no correlation between the average inflorescence height for a plant and the average gene flow distance for that plant. For analysis of pollen characteristics suspected as causal factors of variation in gene flow distance, the variation in gene flow distance was divided into two components: that due to environmental (seasonal) differences, and that due to the differences among genotypes within runs. Variation in both the buoyant properties of the pollen grains and their adhesiveness was significantly partially correlated with environmental variation in gene flow distance, while only the buoyant properties of the pollen grains were significantly partially correlated with among-genotype variation in gene flow distances.     

GENETIC DIFFERENTIATION AMONG CHEWING LOUSE POPULATIONS (MALLOPHAGA: TRICHODECTIDAE) IN A POCKET GOPHER CONTACT ZONE (RODENTIA: GEOMYIDAE)

Genetic variation among populations of chewing lice (Geomydoecus actuosi) was examined in relation to chromosomal and electrophoretic variation among populations of their hosts (Thomomys bottae) at a contact zone. Louse demes were characterized by low levels of genetic heterozygosity (H? = 0.039) that may result from founder effects during primary infestation of hosts, compounded by seasonal reductions in louse population size. Louse populations sampled from different hosts showed high levels of genetic structuring both within and among host localities. Microgeographic differentiation of louse populations is high (mean F_ST = 0.092) suggesting that properties of this host–parasite system promote differentiation of louse populations living on different individual hosts. Among-population differentiation in lice (F_ST = 0.240) was similar to that measured among host populations (F_ST = 0.236), suggesting a close association between gene flow in pocket gophers and gene flow in their lice.     

PHASE III OF WRIGHT'S SHIFTING BALANCE PROCESS AND THE VARIANCE AMONG DEMES IN MIGRATION RATE

Interdemic selection by the differential migration of individuals out from demes of high fitness and into demes of low fitness (Phase III) is one of the most controversial aspects of Wright's Shifting Balance Theory. I derive a relationship between Phase III migration and the interdemic selection differential, S, and show its potential effect on F_ST. The relationship reveals a diversifying effect of interdemic selection by Phase III migration on the genetic structure of a metapopulation. Using experimental metapopulations, I explored the effect of Phase III migration on F_ST by comparing the genetic variance among demes for two different patterns of migration: (1) island model migration and (2) Wright's Phase III migration. Although mean migration rates were the same, I found that the variance among demes in migration rate was significantly higher with Phase III than with island model migration. As a result, F_ST for the frequency of a neutral marker locus was higher with Phase III than it was with island model migration. By increasing F_ST, Phase III enhanced the genetic differentiation among demes for traits not subject to interdemic selection. This feature makes Wright's process different from individual selection which, by reducing effective population size, decreases the genetic variance within demes for all other traits. I discussed this finding in relation to the efficacy of Phase III and random migration for effecting peak shifts, and the contribution of genes with indirect effects to among‐deme variation.     

Colonization, dispersal, and hybridization influence phylogeography of North Atlantic sea urchins (Strongylocentrotus droebachiensis)

We used frequency-based and coalescent-based phylogeographic analysis of sea urchin (Strongylocentrotus droebachiensis) mitochondrial DNA (mtDNA) sequences and previously published microsatellite data to understand the relative influence of colonization and gene flow from older (north Pacific) and younger (northeast Atlantic) sea urchin populations on genetic variation in the northwest Atlantic. We found strong evidence of survival of northwestern Atlantic populations in local Pleistocene glacial refugia: most haplotypes were the same as or closely related to Pacific haplotypes, with deep gene genealogies that reflect divergence times within the northwestern Atlantic that are much older than the last glacial maximum. We detected gene flow across the North Atlantic in the form of haplotypes shared with or recently descended from European populations. We also found evidence of significant introgression of haplotypes from a closely related species (S. pallidus). The relative magnitude of gene flow estimated by coalescent methods (and the effective population size differences among oceanic regions) depended on the genetic marker used. In general, we found very small effective population size in the northeastern Atlantic and high trans-Arctic gene flow between the Pacific and northwestern Atlantic. Both analyses suggested significant back-migration to the Pacific. However, microsatellites more strongly reflected older Pacific migration (with similar effective population sizes across the Arctic), whereas mtDNA sequences appeared to be more sensitive to recent trans- Atlantic dispersal (with larger differences in effective population size). These differences across marker types might have several biological or methodological causes, and they suggest caution in interpretation of the results from a single locus or class of markers.     

Consequences of fragmentation for the ability to adapt to novel environments in experimental <Emphasis Type="Italic">Drosophila</Emphasis> metapopulations

We used experimental populations of Drosophila melanogaster, which had either been subdivided (metapopulations) or kept undivided for 40 generations, to study the consequences of population subdivision for the tolerance and adaptive response after six generations of exposure to novel environmental factors (high temperature, medium with ethanol or salt added) for traits with different genetic architectures. In this setup, we attempted to separate the effects of the loss of fitness due to inbreeding (i.e., the survival upon first exposure to stress) from the loss of adaptive potential due to the lack of genetic variation. To place our experimental results in a more general perspective, we used individual-based simulations combining different options of levels of gene flow, intensity of selection and genetic architecture to derive quantitative hypotheses of the effects of these factors on the adaptive response to stress. We observed that population subdivision resulted in substantial inter-deme variation in tolerance due to redistribution of genetic variation from within demes to among demes. In line with the simulation results, the adaptive response was generally lower in the subdivided than in the undivided populations, particularly so for high temperature. We observed pronounced differences between stress factors that are likely related to the different genetic architectures involved in resistance to these factors. From a conservation genetics viewpoint, our results have two important implications: (i) Long-term fragmentation in combination with restricted gene flow will limit the adaptive potential of individual subpopulations because adaptive variation will become distributed among populations rather than within populations. (ii) The genetic architecture of the trait(s) under selection is of great significance to understand the possible responses to novel stresses that may be expected.     

Geographic and host‐mediated population genetic structure in a cestode parasite of the three‐spined stickleback

Comparative studies of genetic diversity and population structure can shed light on the ecological and evolutionary factors that influence host–parasite interactions. Here we examined whether geography, time and genetic variation in Alaskan three‐spined stickleback (Gasterosteus aculeatus Linneaus) hosts shape the population genetic structure of the diphyllobothridean cestode parasite Schistocephalus solidus (Müller, 1776). Host lineages and haplotypes were identified by sequencing the mitochondrial cytochrome b gene, and host population structure was assessed by Bayesian clustering analysis of allelic variation at 11 microsatellite loci. Parasite population structure was characterized according to allelic variation at eight microsatellite loci. Mantel tests and canonical redundancy analysis were conducted to evaluate the proportion of parasite genetic variation attributable to time and geography vs. host lineage, haplotype, and genotypic cluster. Host and parasite population structure were largely discordant across the study area, probably reflecting differences in gene flow, environmental influences external to the host, and genomic admixture among host lineages. We found that geography explained the greatest proportion of parasite genetic variation, but that variation also reflects time, host lineage, and host haplotype. Associations with host haplotypes suggest that one parasite genotypic cluster exhibits a narrower host range, predominantly infecting the most common host haplotypes, whereas the other parasite cluster infects all haplotypes equally, including rare haplotypes. Although experimental infection trials might prove otherwise, distributional differences in hosts preferentially infected by S. solidus could underlie the observed pattern of population structure.     

Indo-West Pacific patterns of genetic differentiation in the high-dispersal starfish Linckia laevigata

Genetic variation in four natural populations of the starfish Linckia laevigata from the Indo-West Pacific was examined using restriction fragment analysis of a portion of the mtDNA including the control region. Digestion with seven restriction enzymes identified 47 haplotypes in a sample of 326 individuals. Samples collected from reef sites within each location were not significantly differentiated based on Φ_ST or spatial distribution of haplotypes, indicating that dispersal is high over short to moderate distances. Evidence of gene flow is further supported by the low divergence among haplotypes and the lack of any clear geographical structuring among different haplotypes in the gene phylogeny. However, analysis of molecular variance ( AMOVA ), Φ_ST and contingency χ² analyses of the spatial distribution of haplotypes demonstrate the presence of significant broad scale population genetic structure among the four widespread locations examined. RFLP data are consistent with high gene flow between the Philippines and Western Australia and moderate gene flow between the Great Barrier Reef (GBR) and Fiji, but only limited gene flow between either the Philippines or Western Australia and either the GBR or Fiji. The presence of mtDNA structure contrasts with previous allozyme data which suggest that dispersal among widely separated locations is equivalent to dispersal among populations within the highly connected GBR studies. This discordance between patterns of gene flow inferred from these two markers cannot be fully accounted for by differences in effective population size for mtDNA. This might suggest that while mtDNA variation may represent contemporary patterns of gene flow, allozyme variation among populations is yet to reach equilibrium between drift and migration over the range surveyed.     

Microsatellite Analyses of Blacktip Reef Sharks (Carcharhinus melanopterus) in a Fragmented Environment Show Structured Clusters

The population dynamics of shark species are generally poorly described because highly mobile marine life is challenging to investigate. Here we investigate the genetic population structure of the blacktip reef shark (Carcharhinus melanopterus) in French Polynesia. Five demes were sampled from five islands with different inter-island distances (50–1500 km). Whether dispersal occurs between islands frequently enough to prevent moderate genetic structure is unknown. We used 11 microsatellites loci from 165 individuals and a strong genetic structure was found among demes with both F-statistics and Bayesian approaches. This differentiation is correlated with the geographic distance between islands. It is likely that the genetic structure seen is the result of all or some combination of the following: low gene flow, time since divergence, small effective population sizes, and the standard issues with the extent to which mutation models actually fit reality. We suggest low levels of gene flow as at least a partial explanation of the level of genetic structure seen among the sampled blacktip demes. This explanation is consistent with the ecological traits of blacktip reef sharks, and that the suitable habitat for blacktips in French Polynesia is highly fragmented. Evidence for spatial genetic structure of the blacktip demes we studied highlights that similar species may have populations with as yet undetected or underestimated structure. Shark biology and the market for their fins make them highly vulnerable and many species are in rapid decline. Our results add weight to the case that total bans on shark fishing are a better conservation approach for sharks than marine protected area networks.     

THE PROBABILITY OF FIXATION OF A NEW KARYOTYPE IN A CONTINUOUS POPULATION

We investigate the probability of fixation of a chromosome rearrangement in a subdivided population, concentrating on the limit where migration is so large relative to selection (m ? s) that the population can be thought of as being continuously distributed. We study two demes, and one- and two-dimensional populations. For two demes, the probability of fixation in the limit of high migration approximates that of a population with twice the size of a single deme: migration therefore greatly reduces the fixation probability. However, this behavior does not extend to a large array of demes. Then, the fixation probability depends primarily on neighborhood size (Nb), and may be appreciable even with strong selection and free gene flow (≈exp(-B ≈ Nb√s) in one dimension, ≈exp(-B ≈ Nb) in two dimensions). Our results are close to those for the more tractable case of a polygenic character under disruptive selection.     

Local extinction and recolonization, species effective population size, and modern human origins

A primary objection from a population genetics perspective to a multiregional model of modern human origins is that the model posits a large census size, whereas genetic data suggest a small effective population size. The relationship between census size and effective size is complex, but arguments based on an island model of migration show that if the effective population size reflects the number of breeding individuals and the effects of population subdivision, then an effective population size of 10,000 is inconsistent with the census size of 500,000 to 1,000,000 that has been suggested by archeological evidence. However, these models have ignored the effects of population extinction and recolonization, which increase the expected variance among demes and reduce the inbreeding effective population size. Using models developed for population extinction and recolonization, we show that a large census size consistent with the multiregional model can be reconciled with an effective population size of 10,000, but genetic variation among demes must be high, reflecting low interdeme migration rates and a colonization process that involves a small number of colonists or kin-structured colonization. Ethnographic and archeological evidence is insufficient to determine whether such demographic conditions existed among Pleistocene human populations, and further work needs to be done. More realistic models that incorporate isolation by distance and heterogeneity in extinction rates and effective deme sizes also need to be developed. However, if true, a process of population extinction and recolonization has interesting implications for human demographic history.     

Contemporary factors influencing genetic diversity in the Alaska humpback whitefish Coregonus clupeaformis complex

Thirteen microsatellite loci were used to address three hypotheses regarding genetic diversity in the humpback whitefish Coregonus clupeaformis complex in Alaska. The test results provided further insight into the factors influencing C. clupeaformis complex population structure and level of genetic variation. First, the microsatellite data did not provide evidence of two spatially distinct Beringian and Eurasian refugial groups as revealed in previous phylogeographic analyses of mitochondrial DNA variation. Rather, the population structure inferred from the microsatellite variation appears to reveal the influence of factors on a more recent time scale, including gene flow among the refugial groups and isolation of some anadromous and freshwater resident populations. Second, anadromous C. clupeaformis complex collections exhibited higher intra‐population genetic diversity than freshwater resident collections. This outcome is consistent with previous meta analyses suggesting that freshwater resident populations probably have smaller historical effective population sizes and less conspecific gene flow because the habitat tends to be smaller and supports fewer and smaller populations. Finally, the analysis of contemporary immigration rates was consistent with, but did not provide statistical support for, the hypothesis that gene flow among anadromous C. clupeaformis complex populations along coastal Alaska is influenced by the Alaska Coastal Current. Further studies are needed to evaluate gene flow among coastal Alaska C. clupeaformis complex populations.     

GENETIC VARIANCE FOR RATE OF POPULATION INCREASE IN NATURAL POPULATIONS OF FLOUR BEETLES,TRIBOLIUM SPP.

The genetic and ecological effects of population subdividsion were investigated for two wild strains of Tribolium castaneum and two wild strains of T. confusum and compared with the effects of population subdivision on the synthetic laboratory strain of T. castaneum (c-SM), used extensively in earlier experiments. For the c-SM strain, it has been shown repeatedly, for a variety of different population structures (different combinations of effective numbers, N_e, and migration rates, m), that large heritable differences in population growth rate arise among demes during 10 to 15 generations of population subdivision. Because this laboratory strain was synthesized by mass mating several “inbred” strains in 1973 (80 to 100 generations ago), it is possible that it has genetic variation for fitness (measured as the heritable variance among demes in the rate of population increase) unusually large compared to natural populations of flour beetles. In this paper, I report that natural populations of flour beetle exhibit as much or more phenotypic and genetic variation in the effects of population structure on fitness than the laboratory strain, c-SM. The observation of substantial heritable variation for fitness in natural populations is unexpected under additive theory and may be indicative of nonadditive genetic variance.     

The geographic structure of morphological variation in eight species of fiddler crabs (Ocypodidae: genus Uca) from the eastern United States and Mexico

Species with larger geographic distributions are more likely to encounter a greater variety of environmental conditions and barriers to gene flow than geographically‐restricted species. Thus, even closely‐related species with similar life‐history strategies might vary in degree and geographic structure of variation if they differ in geographic range size. In the present study, we investigated this using samples collected across the geographic ranges of eight species of fiddler crabs (Crustacea: Uca) from the Atlantic and Gulf coasts of North America. Morphological variation in the carapace was assessed using geometric morphometric analysis of 945 specimens. Although the eight Uca species exhibit different degrees of intraspecific variation, widespread species do not necessarily exhibit more intraspecific or geographic variation in carapace morphology. Instead, species with more intraspecific variation show stronger morphological divergence among populations. This morphological divergence is partly a result of allometric growth coupled with differences in maximum body size among populations. On average, 10% of total within‐species variation is attributable to allometry. Possible drivers of the remaining morphological differences among populations include gene flow mediated by ocean currents and plastic responses to various environmental stimuli, with isolation‐by‐distance playing a less important role. The results obtained indicate that morphological divergence among populations can occur over shorter distances than expected based on dispersal potential. © 2010 The Linnean Society of London, Biological Journal of the Linnean Society, 2010, 100 , 248–270.     

Evidence of divergent growth rates among populations of the lizard Anolis carolinensis based on experimental manipulations of egg size

Geographic variation in body size is of special interest because it affects nearly all aspects of an organism’s life. I examined whether differences in body size among four populations of the green anole lizard, Anolis carolinensis, were attributable to maternal investment in egg size and/or growth rates of embryos and juveniles. Larger body size and larger egg size relative to female size in the northern part of the range have been documented in this species, and suggested to be adaptive responses to more extreme winters. The current study confirmed the trends in adult size and egg size in the north, but rejected the trend of larger egg size relative to body size in the south. To control for differences in maternal investment in egg size among populations, I performed yolk removals on eggs from two northern populations to produce comparably sized eggs relative to one southern population. This manipulation was designed to minimize the confounding effect of maternal investment in yolk, the primary energy reserves for eggs, so that intrinsic differences in embryonic growth due to metabolism could be investigated. I found that differences in juvenile and, potentially, embryonic growth rates existed among populations of A. carolinensis, both due to and independent of differences in egg size. Juveniles from the northernmost population were bigger not only due to larger egg size, but also due to faster juvenile growth and possibly differences in developmental stage of oviposition or conversion of egg mass to hatchling mass. Larger body size may hold a number of advantages in northern populations of this species, including starvation resistance through winters and better competitive access to food resources and warmer microhabitats.     

Contemporary temperature-driven divergence in a Nordic freshwater fish under conditions commonly thought to hinder adaptation

Background  

Evaluating the limits of adaptation to temperature is important given the IPCC-predicted rise in global temperatures. The rate and scope of evolutionary adaptation can be limited by low genetic diversity, gene flow, and costs associated with adaptive change. Freshwater organisms are physically confined to lakes and rivers, and must therefore deal directly with climate variation and change. In this study, we take advantage of a system characterised by low genetic variation, small population size, gene flow and between-trait trade-offs to study how such conditions affect the ability of a freshwater fish to adapt to climate change. We test for genetically-based differences in developmental traits indicating local adaptation, by conducting a common-garden experiment using embryos and larvae from replicate pairs of sympatric grayling demes that spawn and develop in natural cold and warm water, respectively. These demes have common ancestors from a colonization event 22 generations ago. Consequently, we explore if diversification may occur under severely constraining conditions.     

Tests of two methods for identifying founder effects in metapopulations reveal substantial type II error

Genetic analysis has been promoted as a way to reconstruct recent historical dynamics (“historical demography”) by screening for signatures of events, such as bottlenecks, that disrupt equilibrium patterns of variation. Such analyses might also identify “metapopulation” processes like extinction and recolonization or source-sink dynamics, but this potential remains largely unrealized. Here we use simulations to test the ability of two currently used strategies to distinguish between a set of interconnected subpopulations (demes) that have undergone bottlenecks or extinction and recolonization events (metapopulation dynamics) from a set of static demes. The first strategy, decomposed pairwise regression, provides a holistic test for heterogeneity among demes in their patterns of isolation-by-distance. This method suffered from a type II error rate of 59–100 %, depending on parameter conditions. The second strategy tests for deviations from mutation-drift equilibrium on a deme-by-deme basis to identify sites likely to have experienced recent bottlenecks or founder effects. Although bottleneck tests have good statistical power for single populations with recent population declines, their validity in structured populations has been called into question, and they have not been tested in a metapopulation context with immigration (or colonization) and population recovery. Our simulations of hypothetical metapopulations show that population recovery can rapidly eliminate the statistical signature of a bottleneck, and that moderate levels of gene flow can generate a false signal of recent population growth for demes in equilibrium. Although we did not cover all possible metapopulation scenarios, the performance of the tests was disappointing. Our results indicate that these methods might often fail to identify population bottlenecks and founder effects if population recovery and/or gene flow are influential demographic features of the study system.     

Breeding system and population genetic structure in philomycid slugs (Mollusca: Pulmonata)

Four species of philomycid slugs were collected at 13 sites in Tennessee and one site in Virginia. A total of 361 individuals were examined electrophoretically at 11 allozyme loci. Outcrossing was the primary breeding in all populations of all four species. Genetic similarities among demes, among species and among genera are similar to values obtained in studies of other animals with similar vagility and demography and are in good agreement with systematics based on morphological characters. In the species examined in greatest detail, Philomycus carolinianus , there was substantial genetic differentiation among demes. This differentiation appears to result from genetic drift and restricted gene flow. Genetic similarity between P. carolinianus demes did not correlate with the geographic distance between demes and hierarchical F-statistic analysis showed more than half of the differentiation among demes to result from differentiation at the smallest geographic scale. The Tennessee River system also appears to be a major barrier to gene flow among the P. carolinianus demes sampled.