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.     

THE QUANTITATIVE AND MOLECULAR GENETIC ARCHITECTURE OF A SUBDIVIDED SPECIES

In an effort to elucidate the evolutionary mechanisms that determine the genetic architecture of a species, we have analyzed 17 populations of the microcrustacean Daphnia pulex for levels of genetic variation at the level of life-history characters and molecular markers in the nuclear and mitochondrial genomes. This species is highly subdivided, with approximately 30% of the variation for nuclear molecular markers and 50% of the variation for mitochondrial markers being distributed among populations. The average level of genetic subdivision for quantitative traits is essentially the same as that for nuclear markers, which superficially suggests that the life-history characters are diverging at the neutral rate. However, the existence of a strong correlation between the levels of population subdivision and broadsense heritabilities of individual traits argues against this interpretation, suggesting instead that the among-population divergence of some quantitative traits (most notably body size) is being driven by local adaptation to different environments. The fact that the mean phenotypes of the individual populations are also strongly correlated with local levels of homozygosity indicates that variation in local inbreeding plays a role in population differentiation. Rather than being a passive consequence of local founder effects, levels of homozygosity may be selected for directly for their effects on the phenotype (adaptive inbreeding depression). There is no relationship between the levels of variation within populations for molecular markers and quantitative characters, and this is explained by the fact that the average standing genetic variation for life-history characters in this species is equivalent to only 33 generations of variation generated by mutation.     

Phenotypic and genetic divergence among island populations of sika deer (<Emphasis Type="Italic">Cervus nippon</Emphasis>) in southern Japan: a test of the local adaptation hypothesis

Transplant and common garden experiments have been used in studies on local adaptation, but are difficult to be conducted for large animals with long life span. A previous study on the southern Japanese islands demonstrated that relative limb lengths of sika deer (Cervus nippon) were short on islands with steep slopes. We hypothesized that this morphological variation was evidence for local adaptation, and tested this hypothesis by comparing phenotypic divergence with neutral genetic divergence among eight populations of the sika deer in the southern Japanese islands. Divergence patterns differed between the phenotypic and neutral genetic features. Genetic similarity was high among individuals on Kyushu (OI, KGS, and KGK). Individuals on Tanegashima (TN) and Yakushima (YK) also constituted a group, whereas individuals on Tsushima (TS), Wakamatsujima (WM), and Kuchinoerabujima (KE) formed a genetically distinct group. Phenotypic data indicated that individuals from TS, OI, KGS, and KGK exhibited similarity, whereas individuals on YK formed an isolated group that was separated from the other populations. The degree of phenotypic divergence was larger than that of neutral genetic divergence between TN and YK. These results suggest that divergent selection worked between two of the eight island populations (TN and YK). The morphological trait of captive-bred individuals from TN and YK, which had never experienced their original environments, retained their original morphological features. By combining the results of multiple analyses, we found that the difference in relative limb length between the two populations was consistent with local adaptation hypothesis, although conclusive results were not obtained for the other populations.     

Population Structure Limits Parallel Evolution in Sticklebacks

Population genetic theory predicts that small effective population sizes (N_e) and restricted gene flow limit the potential for local adaptation. In particular, the probability of evolving similar phenotypes based on shared genetic mechanisms (i.e., parallel evolution), is expected to be reduced. We tested these predictions in a comparative genomic study of two ecologically similar and geographically codistributed stickleback species (viz. Gasterosteus aculeatus and Pungitius pungitius). We found that P. pungitius harbors less genetic diversity and exhibits higher levels of genetic differentiation and isolation-by-distance than G. aculeatus. Conversely, G. aculeatus exhibits a stronger degree of genetic parallelism across freshwater populations than P. pungitius: 2,996 versus 379 single nucleotide polymorphisms located within 26 versus 9 genomic regions show evidence of selection in multiple freshwater populations of G. aculeatus and P. pungitius, respectively. Most regions involved in parallel evolution in G. aculeatus showed increased levels of divergence, suggestive of selection on ancient haplotypes. In contrast, haplotypes involved in freshwater adaptation in P. pungitius were younger. In accordance with theory, the results suggest that connectivity and genetic drift play crucial roles in determining the levels and geographic distribution of standing genetic variation, providing evidence that population subdivision limits local adaptation and therefore also the likelihood of parallel evolution.     

Relative influences of historical and contemporary forces shaping the distribution of genetic variation in the Atlantic killifish, Fundulus heteroclitus

A major goal of population genetics research is to identify the relative influences of historical and contemporary processes that serve to structure genetic variation. Most population genetic models assume that populations exist in a state of migration-drift equilibrium. However, in the past this assumption has rarely been verified, and is likely rarely achieved in natural populations. We assessed the equilibrium status at both local and regional scales of the Atlantic killifish, Fundulus heteroclitus . This species is a model organism for the study of adaptive clinal variation, but has also experienced a complicated history of range expansion and secondary contact following allopatric divergence, potentially obscuring the influence of contemporary evolutionary processes. Presumptively neutral genetic markers (microsatellites) demonstrated zones of secondary intergradation among coastal populations centred around northern New Jersey and the Chesapeake Bay region. Analysis of genetic variation indicated isolation by distance among some populations and provided supporting evidence that the Delaware Bay, but not the Chesapeake Bay, has acted as a barrier to dispersal among coastal populations. Bayesian estimates indicated large effective population sizes and low migration rates, and were in good agreement with empirically derived estimates of population and neighbourhood size from mark–recapture studies. These data indicate that populations are not in migration-drift equilibrium at a regional scale, and suggest that contributing factors include large population size combined with relatively low migration rates. These conditions should be considered when interpreting the evolutionary significance of the distribution of genetic variation among F. heteroclitus populations.     

Plasticity matches phenotype to local conditions despite genetic homogeneity across 13 snake populations

In a widespread species, a matching of phenotypic traits to local environmental optima is generally attributed to site-specific adaptation. However, the same matching can occur via adaptive plasticity, without requiring genetic differences among populations. Adult sea kraits (Laticauda saintgironsi) are highly philopatric to small islands, but the entire population within the Neo-Caledonian Lagoon is genetically homogeneous because females migrate to the mainland to lay their eggs at communal sites; recruits disperse before settling, mixing up alleles. Consequently, any matching between local environments (e.g. prey sizes) and snake phenotypes (e.g. body sizes and relative jaw sizes (RJSs)) must be achieved via phenotypic plasticity rather than spatial heterogeneity in gene frequencies. We sampled 13 snake colonies spread along an approximately 200 km northwest–southeast gradient (n > 4500 individuals) to measure two morphological features that affect maximum ingestible prey size in gape-limited predators: body size and RJS. As proxies of habitat quality (HQ), we used protection status, fishing pressure and lagoon characteristics (lagoon width and distance of islands to the barrier reef). In both sexes, spatial variation in body sizes and RJSs was linked to HQ; albeit in different ways, consistent with sex-based divergences in foraging ecology. Strong spatial divergence in morphology among snake colonies, despite genetic homogeneity, supports the idea that phenotypic plasticity can facilitate speciation by creating multiple phenotypically distinct subpopulations shaped by their environment.     

Rate of Adaptation in Sexuals and Asexuals: A Solvable Model of the Fisher–Muller Effect

The adaptation of large asexual populations is hampered by the competition between independently arising beneficial mutations in different individuals, which is known as clonal interference. In classic work, Fisher and Muller proposed that recombination provides an evolutionary advantage in large populations by alleviating this competition. Based on recent progress in quantifying the speed of adaptation in asexual populations undergoing clonal interference, we present a detailed analysis of the Fisher–Muller mechanism for a model genome consisting of two loci with an infinite number of beneficial alleles each and multiplicative (nonepistatic) fitness effects. We solve the deterministic, infinite population dynamics exactly and show that, for a particular, natural mutation scheme, the speed of adaptation in sexuals is twice as large as in asexuals. This result is argued to hold for any nonzero value of the rate of recombination. Guided by the infinite population result and by previous work on asexual adaptation, we postulate an expression for the speed of adaptation in finite sexual populations that agrees with numerical simulations over a wide range of population sizes and recombination rates. The ratio of the sexual to asexual adaptation speed is a function of population size that increases in the clonal interference regime and approaches 2 for extremely large populations. The simulations also show that the imbalance between the numbers of accumulated mutations at the two loci is strongly suppressed even by a small amount of recombination. The generalization of the model to an arbitrary number L of loci is briefly discussed. If each offspring samples the alleles at each locus from the gene pool of the whole population rather than from two parents, the ratio of the sexual to asexual adaptation speed is approximately equal to L in large populations. A possible realization of this scenario is the reassortment of genetic material in RNA viruses with L genomic segments.     

Genetic Diversity,Heteroplasmy, and Recombination in Mitochondrial Genomes of Daphnia pulex,Daphnia pulicaria,and Daphnia obtusa

Genetic variants of mitochondrial DNA at the individual (heteroplasmy) and population (polymorphism) levels provide insight into their roles in multiple cellular and evolutionary processes. However, owing to the paucity of genome-wide data at the within-individual and population levels, the broad patterns of these two forms of variation remain poorly understood. Here, we analyze 1,804 complete mitochondrial genome sequences from Daphnia pulex, Daphnia pulicaria, and Daphnia obtusa. Extensive heteroplasmy is observed in D. obtusa, where the high level of intraclonal divergence must have resulted from a biparental-inheritance event, and recombination in the mitochondrial genome is apparent, although perhaps not widespread. Global samples of D. pulex reveal remarkably low mitochondrial effective population sizes, <3% of those for the nuclear genome. In addition, levels of population diversity in mitochondrial and nuclear genomes are uncorrelated across populations, suggesting an idiosyncratic evolutionary history of mitochondria in D. pulex. These population-genetic features appear to be a consequence of background selection associated with highly deleterious mutations arising in the strongly linked mitochondrial genome, which is consistent with polymorphism and divergence data suggesting a predominance of strong purifying selection. Nonetheless, the fixation of mildly deleterious mutations in the mitochondrial genome also appears to be driving positive selection on genes encoded in the nuclear genome whose products are deployed in the mitochondrion.     

The genetic architecture of adaptation under migration-selection balance

Many ecologically important traits have a complex genetic basis, with the potential for mutations at many different genes to shape the phenotype. Even so, studies of local adaptation in heterogeneous environments sometimes find that just a few quantitative trait loci (QTL) of large effect can explain a large percentage of observed differences between phenotypically divergent populations. As high levels of gene flow can swamp divergence at weakly selected alleles, migration-selection-drift balance may play an important role in shaping the genetic architecture of local adaptation. Here, we use analytical approximations and individual-based simulations to explore how genetic architecture evolves when two populations connected by migration experience stabilizing selection toward different optima. In contrast to the exponential distribution of allele effect sizes expected under adaptation without migration (Orr 1998), we find that adaptation with migration tends to result in concentrated genetic architectures with fewer, larger, and more tightly linked divergent alleles. Even if many small alleles contribute to adaptation at the outset, they tend to be replaced by a few large alleles under prolonged bouts of stabilizing selection with migration. All else being equal, we also find that stronger selection can maintain linked clusters of locally adapted alleles over much greater map distances than weaker selection. The common empirical finding of QTL of large effect is shown to be expected with migration in a heterogeneous landscape, and these QTL may often be composed of several tightly linked alleles of smaller effect.     

Environmental heterogeneity enhances clonal interference

Clonal interference (CI) is a phenomenon that may be important in several asexual microbes. It occurs when population sizes are large and mutation rates to new beneficial alleles are of significant magnitude. Here we explore the role of gene flow and spatial heterogeneity in selection strength in the adaptation of asexuals. We consider a subdivided population of individuals that are adapting, through new beneficial mutations, and that migrate between different patches. The fitness effect of each mutation depends on the patch and all mutations considered are assumed to be unconditionally beneficial. We find that spatial variation in selection pressure affects the rate of adaptive evolution and its qualitative effects depend on the level of gene flow. In particular, we find that both low migration and high levels of heterogeneity lead to enhanced CI. In contrast, for high levels of migration the rate of fixation of adaptive mutations is higher when environmental heterogeneity is present. In addition, we observe that the level of fitness variation is higher and simultaneous fixation of multiple mutations tends to occur in the regime of low migration rates and high heterogeneity.     

Genetic and phylogenetic consequences of island biogeography

Abstract.— Island biogeography theory predicts that the number of species on an island should increase with island size and decrease with island distance to the mainland. These predictions are generally well supported in comparative and experimental studies. These ecological, equilibrium predictions arise as a result of colonization and extinction processes. Because colonization and extinction are also important processes in evolution, we develop methods to test evolutionary predictions of island biogeography. We derive a population genetic model of island biogeography that incorporates island colonization, migration of individuals from the mainland, and extinction of island populations. The model provides a means of estimating the rates of migration and extinction from population genetic data. This model predicts that within an island population the distribution of genetic divergences with respect to the mainland source population should be bimodal, with much of the divergence dating to the colonization event. Across islands, this model predicts that populations on large islands should be on average more genetically divergent from mainland source populations than those on small islands. Likewise, populations on distant islands should be more divergent than those on close islands. Published observations of a larger proportion of endemic species on large and distant islands support these predictions.     

Population history in Arabidopsis halleri using multilocus analysis

A. halleri is a psuedometallophyte with a patchy distribution in Europe and is often spread by human activity. To determine the population history and whether this history is consistent with potential human effects, we surveyed nucleotide variation using 24 loci from 12 individuals in a large A. halleri population. The means of total and silent nucleotide variation (θ_W) are within the range expected for the species. The population genetic neutrality tests Tajima’s D and Wall’s B had significant composite results rejecting panmixia, and Approximate Bayesian Computation analysis revealed that a subdivision model better explained the variation than the standard neutral model, refugia (or admixture), bottleneck or change of population size models. A categorical regression analysis further supports the subdivision model, and under the subdivision model, the neutrality tests are no longer significant. The best support was for two source populations, a situation consistent with the mixing of two populations possibly mediated by human activity. This scenario might limit the genetic diversity and adaptive potential of the population. The non‐neutral population variation described here should be considered in bioinformatic searches for adaptation.     

Do highly divergent loci reside in genomic regions affecting reproductive isolation? A test using next-generation sequence data in Timema stick insects

ABSTRACT: BACKGROUND: Genetic divergence during speciation with gene flow is heterogeneous across the genome, with some regions exhibiting stronger differentiation than others. Exceptionally differentiated regions are often assumed to experience reduced introgression, i.e., reduced flow of alleles from one population into another because such regions are affected by divergent selection or cause reproductive isolation. In contrast, the remainder of the genome can be homogenized by high introgression. Although many studies have documented variation across the genome in genetic differentiation, there are few tests of this hypothesis that explicitly quantify introgression. Here, we provide such a test using 38,304 SNPs in populations of Timema cristinae stick insects. We quantify whether loci that are highly divergent between geographically separated ('allopatric') populations exhibit unusual patterns of introgression in admixed populations. To the extent this is true, highly divergent loci between allopatric populations contribute to reproductive isolation in admixed populations. RESULTS: As predicted, we find a substantial association between locus-specific divergence between allopatric populations and locus-specific introgression in admixed populations. However, many loci depart from this relationship, sometimes strongly so. We also report evidence for selection against foreign alleles due to local adaptation. CONCLUSIONS: Loci that are strongly differentiated between allopatric populations sometimes contribute to reproductive isolation in admixed populations. However, geographic variation in selection and local adaptation, in aspects of genetic architecture (such as organization of genes, recombination rate variation, number and effect size of variants contributing to adaptation, etc.), and in stochastic evolutionary processes such as drift can cause strong differentiation of loci that do not always contribute to reproductive isolation. The results have implications for the theory of 'genomic islands of speciation'.     

Indo-Pacific population structure and evolutionary history of the coconut crab Birgus latro

Mitochondrial DNA variation was used to examine population structure in a widespread, marine-dispersed species, Birgus latro . Crabs were collected from eight locations throughout the species' Indo-Pacific distribution. Purified mtDNA from 160 individuals was cut with five restriction enzymes, revealing high haplotype diversity (0.96) and moderate nucleotide diversity (0.75%). Island populations from the Indian Ocean (Christmas I.) and Pacific Ocean were significantly different ( G _ST= 0.37) and had distinct mtDNA lineages with a net sequence divergence of 1.4%. Pacific island populations had diverged in a manner consistent with isolation by distance, with only the most peripheral populations being significantly different. The results for mtDNA are largely concordant with those from allozymes, although estimates of gene flow between the Indian and Pacific Oceans were much lower when based on mtDNA. The mtDNA phylogeny also permitted a deeper examination of the evolutionary and demographic history of Birgus latro . Long-term separation of populations is evident in the complete phylogenetic subdivision of mtDNA lineages between the Indian and Pacific Ocean populations sampled. The starlike phylogeny of alleles from the Pacific suggests a rapid population expansion in the Pacific during the Pleistocene. Including information about allele phylogeny, as well as distribution and frequency, obscured contemporary population structure, but provided unique insights into the evolutionary history of the species.     

Population amalgamation and genetic variation: observations on artificially agglomerated tribal populations of Central and South America.

The interpretation of data on genetic variation with regard to the relative roles of different evolutionary factors that produce and maintain genetic variation depends critically on our assumptions concerning effective population size and the level of migration between neighboring populations. In humans, recent population growth and movements of specific ethnic groups across wide geographic areas mean that any theory based on assumptions of constant population size and absence of substructure is generally untenable. We examine the effects of population subdivision on the pattern of protein genetic variation in a total sample drawn from an artificial agglomerate of 12 tribal populations of Central and South America, analyzing the pooled sample as though it were a single population. Several striking findings emerge. (1) Mean heterozygosity is not sensitive to agglomeration, but the number of different alleles (allele count) is inflated, relative to neutral mutation/drift/equilibrium expectation. (2) The inflation is most serious for rare alleles, especially those which originally occurred as tribally restricted "private" polymorphisms. (3) The degree of inflation is an increasing function of both the number of populations encompassed by the sample and of the genetic divergence among them. (4) Treating an agglomerated population as though it were a panmictic unit of long standing can lead to serious biases in estimates of mutation rates, selection pressures, and effective population sizes. Current DNA studies indicate the presence of numerous genetic variants in human populations. The findings and conclusions of this paper are all fully applicable to the study of genetic variation at the DNA level as well.     

Patterns of parapatric speciation

Abstract. Geographic variation may ultimately lead to the splitting of a subdivided population into reproductively isolated units in spite of migration. Here, we consider how the waiting time until the first split and its location depend on different evolutionary factors including mutation, migration, random genetic drift, genetic architecture, and the geometric structure of the habitat. We perform large-scale, individual-based simulations using a simple model of reproductive isolation based on a classical view that reproductive isolation evolves as a by-product of genetic divergence. We show that rapid parapatric speciation on the time scale of a few hundred to a few thousand generations is plausible even when neighboring subpopulations exchange several individuals each generation. Divergent selection for local adaptation is not required for rapid speciation. Our results substantiates the claims that species with smaller range sizes (which are characterized by smaller local densities and reduced dispersal ability) should have higher speciation rates. If mutation rate is small, local abundances are low, or substantial genetic changes are required for reproductive isolation, then central populations should be the place where most splits take place. With high mutation rates, high local densities, or with moderate genetic changes sufficient for reproductive isolation, speciation events are expected to involve mainly peripheral populations.     

Evolutionary and statistical properties of three genetic distances

Many genetic distances have been developed to summarize allele frequency differences between populations. I review the evolutionary and statistical properties of three popular genetic distances: DS, DA, and theta;, using computer simulation of two simple evolutionary histories: an isolation model of population divergence and an equilibrium migration model. The effect of effective population size, mutation rate, and mutation mechanism upon the parametric value between pairs of populations in these models explored, and the unique properties of each distance are described. The effect of these evolutionary parameters on study design is also investigated and similar results are found for each genetic distance in each model of evolution: large sample sizes are warranted when populations are relatively genetically similar; and loci with more alleles produce better estimates of genetic distance.     

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.     

Are habitat fragmentation,local adaptation and isolation‐by‐distance driving population divergence in wild rice Oryza rufipogon?

Habitat fragmentation weakens the connection between populations and is accompanied with isolation by distance (IBD) and local adaptation (isolation by adaptation, IBA), both leading to genetic divergence between populations. To understand the evolutionary potential of a population and to formulate proper conservation strategies, information on the roles of IBD and IBA in driving population divergence is critical. The putative ancestor of Asian cultivated rice (Oryza sativa) is endangered in China due to habitat loss and fragmentation. We investigated the genetic variation in 11 Chinese Oryza rufipogon populations using 79 microsatellite loci to infer the effects of habitat fragmentation, IBD and IBA on genetic structure. Historical and current gene flows were found to be rare (m_h = 0.0002–0.0013, m_c = 0.007–0.029), indicating IBD and resulting in a high level of population divergence (F_ST = 0.343). High within‐population genetic variation (H_E = 0.377–0.515), relatively large effective population sizes (N_e = 96–158), absence of bottlenecks and limited gene flow were found, demonstrating little impact of recent habitat fragmentation on these populations. Eleven gene‐linked microsatellite loci were identified as outliers, indicating local adaptation. Hierarchical AMOVA and partial Mantel tests indicated that population divergence of Chinese O. rufipogon was significantly correlated with environmental factors, especially habitat temperature. Common garden trials detected a significant adaptive population divergence associated with latitude. Collectively, these findings imply that IBD due to historical rather than recent fragmentation, followed by local adaptation, has driven population divergence in O. rufipogon.