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.     

Tracing the Origin and Evolutionary Fate of Recent Gene Retrocopies in Natural Populations of the House Mouse

Although the contribution of retrogenes to the evolution of genes and genomes has long been recognized, the evolutionary patterns of very recently derived retrocopies that are still polymorphic within natural populations have not been much studied so far. We use here a set of 2,025 such retrocopies in nine house mouse populations from three subspecies (Mus musculus domesticus, M. m. musculus, and M. m. castaneus) to trace their origin and evolutionary fate. We find that ancient house-keeping genes are significantly more likely to generate retrocopies than younger genes and that the propensity to generate a retrocopy depends on its level of expression in the germline. Although most retrocopies are detrimental and quickly purged, we focus here on the subset that appears to be neutral or even adaptive. We show that retrocopies from X-chromosomal parental genes have a higher likelihood to reach elevated frequencies in the populations, confirming the notion of adaptive effects for “out-of-X” retrogenes. Also, retrocopies in intergenic regions are more likely to reach higher population frequencies than those in introns of genes, implying a more detrimental effect when they land within transcribed regions. For a small subset of retrocopies, we find signatures of positive selection, indicating they were involved in a recent adaptation process. We show that the population-specific distribution pattern of retrocopies is phylogenetically informative and can be used to infer population history with a better resolution than with SNP markers.     

Ecological opportunities and specializations shaped genetic divergence in a highly mobile marine top predator

Environmental conditions can shape genetic and morphological divergence. Release of new habitats during historical environmental changes was a major driver of evolutionary diversification. Here, forces shaping population structure and ecotype differentiation (‘pelagic’ and ‘coastal’) of bottlenose dolphins in the North-east Atlantic were investigated using complementary evolutionary and ecological approaches. Inference of population demographic history using approximate Bayesian computation indicated that coastal populations were likely founded by the Atlantic pelagic population after the Last Glacial Maxima probably as a result of newly available coastal ecological niches. Pelagic dolphins from the Atlantic and the Mediterranean Sea likely diverged during a period of high productivity in the Mediterranean Sea. Genetic differentiation between coastal and pelagic ecotypes may be maintained by niche specializations, as indicated by stable isotope and stomach content analyses, and social behaviour. The two ecotypes were only weakly morphologically segregated in contrast to other parts of the World Ocean. This may be linked to weak contrasts between coastal and pelagic habitats and/or a relatively recent divergence. We suggest that ecological opportunity to specialize is a major driver of genetic and morphological divergence. Combining genetic, ecological and morphological approaches is essential to understanding the population structure of mobile and cryptic species.     

Patterns of genomic and allochronic strain divergence in the fall armyworm,Spodoptera frugiperda (J.E. Smith)

Speciation is the process through which reproductive isolation develops between distinct populations. Because this process takes time, speciation studies often necessarily examine populations within a species that are at various stages of divergence. The fall armyworm, Spodoptera frugiperda (J.E. Smith), is comprised of two strains (R = Rice & C = Corn) that serve as a novel system to explore population divergence in sympatry. Here, we use ddRADSeq data to show that fall armyworm strains in the field are largely genetically distinct, but some interstrain hybridization occurs. Although we detected F1 hybrids of both R‐ and C‐strain maternal origin, only hybrids with R‐strain mtDNA were found to contribute to subsequent generations, possibly indicating a unidirectional barrier to gene flow. Although these strains have been previously defined as “host plant‐associated,” we recovered an equal proportion of R‐ and C‐strain moths in fields dominated by C‐strain host plants. As an alternative to host‐associated divergence, we tested the hypothesis that differences in nightly activity patterns could account for reproductive isolation by genotyping temporally collected moths. Our data indicates that strains exhibit a significant shift in the timing of their nightly activities in the field. This divergence in phenology creates a prezygotic reproductive barrier that likely maintains the genetic isolation between strains. Thus, we conclude that it may be ecologically inaccurate to refer to the C‐ and R‐ strain as “host‐associated” and they should more appropriately be considered “allochronic strains.”     

Genome‐wide differentiation in closely related populations: the roles of selection and geographic isolation

Population divergence in geographic isolation is due to a combination of factors. Natural and sexual selection may be important in shaping patterns of population differentiation, a pattern referred to as ‘isolation by adaptation’ (IBA). IBA can be complementary to the well‐known pattern of ‘isolation by distance’ (IBD), in which the divergence of closely related populations (via any evolutionary process) is associated with geographic isolation. The barn swallow Hirundo rustica complex comprises six closely related subspecies, where divergent sexual selection is associated with phenotypic differentiation among allopatric populations. To investigate the relative contributions of selection and geographic distance to genome‐wide differentiation, we compared genotypic and phenotypic variation from 350 barn swallows sampled across eight populations (28 pairwise comparisons) from four different subspecies. We report a draft whole‐genome sequence for H. rustica, to which we aligned a set of 9493 single nucleotide polymorphisms (SNPs). Using statistical approaches to control for spatial autocorrelation of phenotypic variables and geographic distance, we find that divergence in traits related to migratory behaviour and sexual signalling, as well as geographic distance, together explain over 70% of genome‐wide divergence among populations. Controlling for IBD, we find 42% of genomewide divergence is attributable to IBA through pairwise differences in traits related to migratory behaviour and sexual signalling alone. By (i) combining these results with prior studies of how selection shapes morphological differentiation and (ii) accounting for spatial autocorrelation, we infer that morphological adaptation plays a large role in shaping population‐level differentiation in this group of closely related populations.     

Historical connectivity,contemporary isolation and local adaptation in a widespread but discontinuously distributed species endemic to Taiwan,Rhododendron oldhamii (Ericaceae)

Elucidation of the evolutionary processes that constrain or facilitate adaptive divergence is a central goal in evolutionary biology, especially in non-model organisms. We tested whether changes in dynamics of gene flow (historical vs contemporary) caused population isolation and examined local adaptation in response to environmental selective forces in fragmented Rhododendron oldhamii populations. Variation in 26 expressed sequence tag-simple sequence repeat loci from 18 populations in Taiwan was investigated by examining patterns of genetic diversity, inbreeding, geographic structure, recent bottlenecks, and historical and contemporary gene flow. Selection associated with environmental variables was also examined. Bayesian clustering analysis revealed four regional population groups of north, central, south and southeast with significant genetic differentiation. Historical bottlenecks beginning 9168–13,092 years ago and ending 1584–3504 years ago were revealed by estimates using approximate Bayesian computation for all four regional samples analyzed. Recent migration within and across geographic regions was limited. However, major dispersal sources were found within geographic regions. Altitudinal clines of allelic frequencies of environmentally associated positively selected outliers were found, indicating adaptive divergence. Our results point to a transition from historical population connectivity toward contemporary population isolation and divergence on a regional scale. Spatial and temporal dispersal differences may have resulted in regional population divergence and local adaptation associated with environmental variables, which may have played roles as selective forces at a regional scale.     

Natural Selection and Neutral Evolution Jointly Drive Population Divergence between Alpine and Lowland Ecotypes of the Allopolyploid Plant Anemone multifida (Ranunculaceae)

Population differentiation can be driven in large part by natural selection, but selectively neutral evolution can play a prominent role in shaping patters of population divergence. The decomposition of the evolutionary history of populations into the relative effects of natural selection and selectively neutral evolution enables an understanding of the causes of population divergence and adaptation. In this study, we examined heterogeneous genomic divergence between alpine and lowland ecotypes of the allopolyploid plant, Anemone multifida. Using peak height and dominant AFLP data, we quantified population differentiation at non-outlier (neutral) and outlier loci to determine the potential contribution of natural selection and selectively neutral evolution to population divergence. We found 13 candidate loci, corresponding to 2.7% of loci, with signatures of divergent natural selection between alpine and lowland populations and between alpine populations (Fst  = 0.074–0.445 at outlier loci), but neutral population differentiation was also evident between alpine populations (F_ST  = 0.041–0.095 at neutral loci). By examining population structure at both neutral and outlier loci, we determined that the combined effects of selection and neutral evolution are associated with the divergence of alpine populations, which may be linked to extreme abiotic conditions and isolation between alpine sites. The presence of outlier levels of genetic variation in structured populations underscores the importance of separately analyzing neutral and outlier loci to infer the relative role of divergent natural selection and neutral evolution in population divergence.     

Effects of insularity on genetic diversity within and among natural populations

We conducted a quantitative literature review of genetic diversity (GD) within and among populations in relation to categorical population size and isolation (together referred to as “insularity”). Using populations from within the same studies, we were able to control for between‐study variation in methodology, as well as demographic and life histories of focal species. Contrary to typical expectations, insularity had relatively minor effects on GD within and among populations, which points to the more important role of other factors in shaping evolutionary processes. Such effects of insularity were sometimes seen—particularly in study systems where GD was already high overall. That is, insularity influenced GD in a study system when GD was high even in non‐insular populations of the same study system—suggesting an important role for the “scope” of influences on GD. These conclusions were more robust for within population GD versus among population GD, although several biases might underlie this difference. Overall, our findings indicate that population‐level genetic assumptions need to be tested rather than assumed in nature, particularly for topics underlying current conservation management practices.     

Mantel test in population genetics

The comparison of genetic divergence or genetic distances, estimated by pairwise F_ST and related statistics, with geographical distances by Mantel test is one of the most popular approaches to evaluate spatial processes driving population structure. There have been, however, recent criticisms and discussions on the statistical performance of the Mantel test. Simultaneously, alternative frameworks for data analyses are being proposed. Here, we review the Mantel test and its variations, including Mantel correlograms and partial correlations and regressions. For illustrative purposes, we studied spatial genetic divergence among 25 populations of Dipteryx alata (“Baru”), a tree species endemic to the Cerrado, the Brazilian savannas, based on 8 microsatellite loci. We also applied alternative methods to analyze spatial patterns in this dataset, especially a multivariate generalization of Spatial Eigenfunction Analysis based on redundancy analysis. The different approaches resulted in similar estimates of the magnitude of spatial structure in the genetic data. Furthermore, the results were expected based on previous knowledge of the ecological and evolutionary processes underlying genetic variation in this species. Our review shows that a careful application and interpretation of Mantel tests, especially Mantel correlograms, can overcome some potential statistical problems and provide a simple and useful tool for multivariate analysis of spatial patterns of genetic divergence.     

Contemporary evolutionary divergence for a protected species following assisted colonization

Background

Contemporary evolution following assisted colonization may increase the probability of persistence for refuge populations established as a bet-hedge for protected species. Such refuge populations are considered “genetic replicates” that might be used for future re-colonization in the event of a catastrophe in the native site. Although maladaptive evolutionary divergence of captive populations is well recognized, evolutionary divergence of wild refuge populations may also occur on contemporary time scales. Thus, refuge populations may lose their “value” as true genetic replicates of the native population. Here, we show contemporary evolutionary divergence in body shape in an approximately 30-year old refuge population of the protected White Sands pupfish (Cyprinodon tularosa) resulting in a body-shape mismatch with its native environment.

Methodology/Principal Findings

Geometric morphometic data were collected from C. tularosa cultures raised in experimental mesocosms. Cultures were initiated with fish from the two native populations, plus hybrids, in high or low salinity treatments representing the salinities of the two native habitats. We found that body shape was heritable and that shape variation due to phenotypic plasticity was small compared to shape variation due to population source. C. tularosa from the high salinity population retained slender body shapes and fish from the low salinity population retained deep body shapes, irrespective of mesocosm salinity. These data suggest that the observed divergence of a recently established pupfish population was not explained by plasticity. An analysis of microsatellite variation indicated that no significant genetic drift occurred in the refuge population, further supporting the adaptive nature of changes in body shape. These lines of evidence suggest that body shape divergence of the refuge population reflects a case of contemporary evolution (over a 30-year period).

Conclusions/Significance

These results suggest assisted colonization can introduce novel, and/or relaxed selection, and lead to unintended evolutionary divergence.     

Measuring genetic distances between breeds: use of some distances in various short term evolution models

Many works demonstrate the benefits of using highly polymorphic markers such as microsatellites in order to measure the genetic diversity between closely related breeds. But it is sometimes difficult to decide which genetic distance should be used. In this paper we review the behaviour of the main distances encountered in the literature in various divergence models. In the first part, we consider that breeds are populations in which the assumption of equilibrium between drift and mutation is verified. In this case some interesting distances can be expressed as a function of divergence time, t, and therefore can be used to construct phylogenies. Distances based on allele size distribution (such as (δμ)² and derived distances), taking a mutation model of microsatellites, the Stepwise Mutation Model, specifically into account, exhibit large variance and therefore should not be used to accurately infer phylogeny of closely related breeds. In the last section, we will consider that breeds are small populations and that the divergence times between them are too small to consider that the observed diversity is due to mutations: divergence is mainly due to genetic drift. Expectation and variance of distances were calculated as a function of the Wright-Malécot inbreeding coefficient, F. Computer simulations performed under this divergence model show that the Reynolds distance [57]is the best method for very closely related breeds.     

Reconstruction of Microbial Haplotypes by Integration of Statistical and Physical Linkage in Scaffolding

DNA sequencing technologies provide unprecedented opportunities to analyze within-host evolution of microorganism populations. Often, within-host populations are analyzed via pooled sequencing of the population, which contains multiple individuals or “haplotypes.” However, current next-generation sequencing instruments, in conjunction with single-molecule barcoded linked-reads, cannot distinguish long haplotypes directly. Computational reconstruction of haplotypes from pooled sequencing has been attempted in virology, bacterial genomics, metagenomics, and human genetics, using algorithms based on either cross-host genetic sharing or within-host genomic reads. Here, we describe PoolHapX, a flexible computational approach that integrates information from both genetic sharing and genomic sequencing. We demonstrated that PoolHapX outperforms state-of-the-art tools tailored to specific organismal systems, and is robust to within-host evolution. Importantly, together with barcoded linked-reads, PoolHapX can infer whole-chromosome-scale haplotypes from 50 pools each containing 12 different haplotypes. By analyzing real data, we uncovered dynamic variations in the evolutionary processes of within-patient HIV populations previously unobserved in single position-based analysis.     

Genetic structure of chimpanzee populations

Little is known about the history and population structure of our closest living relatives, the chimpanzees, in part because of an extremely poor fossil record. To address this, we report the largest genetic study of the chimpanzees to date, examining 310 microsatellites in 84 common chimpanzees and bonobos. We infer three common chimpanzee populations, which correspond to the previously defined labels of “western,” “central,” and “eastern,” and find little evidence of gene flow between them. There is tentative evidence for structure within western chimpanzees, but we do not detect distinct additional populations. The data also provide historical insights, demonstrating that the western chimpanzee population diverged first, and that the eastern and central populations are more closely related in time.     

Divergence time of the two regional medaka populations in Japan as a new time scale for comparative genomics of vertebrates

The southern and northern Japanese populations of the medaka fish provide useful tools to gain insights into the comparative genomics and speciation of vertebrates, because they can breed to produce healthy and fertile offspring despite their highly divergent genetic backgrounds compared with those of human–chimpanzee. Comparative genomics analysis has suggested that such large genetic differences between the two populations are caused by higher molecular evolutionary rates among the medakas than those of the hominids. The argument, however, was based on the assumption that the two Japanese populations diverged approximately at the same time (4.0–4.7 Myr ago) as the human–chimpanzee lineage (5.0–6.0 Myr ago). This can be misleading, because the divergence time of the two populations was calculated based on estimated, extremely higher molecular evolutionary rates of other fishes with an implicit assumption of a global molecular clock. Here we show that our estimate, based on a Bayesian relaxed molecular-clock analysis of whole mitogenome sequences from 72 ray-finned fishes (including 14 medakas), is about four times older than that of the previous study (18 Myr). This remarkably older estimate can be reconciled with the vicariant events of the Japanese archipelago, and the resulting rates of molecular evolution are almost identical between the medaka and hominid lineages. Our results further highlight the fact that reproductive isolation may not evolve despite a long period of geographical isolation.     

The Effect of Linkage on Establishment and Survival of Locally Beneficial Mutations

We study invasion and survival of weakly beneficial mutations arising in linkage to an established migration–selection polymorphism. Our focus is on a continent–island model of migration, with selection at two biallelic loci for adaptation to the island environment. Combining branching and diffusion processes, we provide the theoretical basis for understanding the evolution of islands of divergence, the genetic architecture of locally adaptive traits, and the importance of so-called “divergence hitchhiking” relative to other mechanisms, such as “genomic hitchhiking”, chromosomal inversions, or translocations. We derive approximations to the invasion probability and the extinction time of a de novo mutation. Interestingly, the invasion probability is maximized at a nonzero recombination rate if the focal mutation is sufficiently beneficial. If a proportion of migrants carries a beneficial background allele, the mutation is less likely to become established. Linked selection may increase the survival time by several orders of magnitude. By altering the timescale of stochastic loss, it can therefore affect the dynamics at the focal site to an extent that is of evolutionary importance, especially in small populations. We derive an effective migration rate experienced by the weakly beneficial mutation, which accounts for the reduction in gene flow imposed by linked selection. Using the concept of the effective migration rate, we also quantify the long-term effects on neutral variation embedded in a genome with arbitrarily many sites under selection. Patterns of neutral diversity change qualitatively and quantitatively as the position of the neutral locus is moved along the chromosome. This will be useful for population-genomic inference. Our results strengthen the emerging view that physically linked selection is biologically relevant if linkage is tight or if selection at the background locus is strong.     

Characterizing genic and nongenic molecular markers: comparison of microsatellites and SNPs

The implications of transitioning to single nucleotide polymorphism (SNPs) from microsatellite markers (MSs) have been investigated in a number of population genetics studies, but the effect of genomic location on the amount of information each type of marker reveals has not been explored in detail. We developed novel SNP markers flanking 1 kb regions of 13 genic (within gene or <1 kb away from gene) and 13 nongenic (>10 kb from annotated gene) MSs in the threespine stickleback genome to obtain comparable data for both types of markers. We analysed patterns of genetic diversity and divergence on various geographic scales after converting the SNP loci within each genomic region into haplotypes. Marker type (SNP haplotype or MS) and location (genic or nongenic) significantly affected most estimates of population diversity and divergence. Between‐lineage divergence was significantly higher in SNP haplotypes (genic and nongenic), however, within‐lineage divergence was similar between marker types. Most divergence and diversity measures were uncorrelated between markers, except for population differentiation which was correlated between MSs and SNP haplotypes (both genic and nongenic). Broad‐scale population structure and assignment were similarly resolved by both marker types, however, only the MSs were able to delimit fine‐scale population structuring, particularly when genic and nongenic markers were combined. These results demonstrate that estimates of genetic variability and differentiation among populations can be strongly influenced by marker type, their genomic location in relation to genes and by the interaction of these two factors. This highlights the importance of having an awareness of the inherent strengths and limitations associated with different molecular tools to select the most appropriate methods for accurately addressing various ecological and evolutionary questions.     

The computer program STRUCTURE does not reliably identify the main genetic clusters within species: simulations and implications for human population structure

One of the primary goals of population genetics is to succinctly describe genetic relationships among populations, and the computer program STRUCTURE is one of the most frequently used tools for doing so. The mathematical model used by STRUCTURE was designed to sort individuals into Hardy–Weinberg populations, but the program is also frequently used to group individuals from a large number of populations into a small number of clusters that are supposed to represent the main genetic divisions within species. In this study, I used computer simulations to examine how well STRUCTURE accomplishes this latter task. Simulations of populations that had a simple hierarchical history of fragmentation showed that when there were relatively long divergence times within evolutionary lineages, the clusters created by STRUCTURE were frequently not consistent with the evolutionary history of the populations. These difficulties can be attributed to forcing STRUCTURE to place individuals into too few clusters. Simulations also showed that the clusters produced by STRUCTURE can be strongly influenced by variation in sample size. In some circumstances, STRUCTURE simply put all of the individuals from the largest sample in the same cluster. A reanalysis of human population structure suggests that the problems I identified with STRUCTURE in simulations may have obscured relationships among human populations—particularly genetic similarity between Europeans and some African populations.     

Population genetic structure,differentiation, and diversity in Tetrix subulata pygmy grasshoppers: roles of population size and immigration

Genetic diversity within and among populations and species is influenced by complex demographic and evolutionary processes. Despite extensive research, there is no consensus regarding how landscape structure, spatial distribution, gene flow, and population dynamics impact genetic composition of natural populations. Here, we used amplified fragment length polymorphisms (AFLPs) to investigate effects of population size, geographic isolation, immigration, and gene flow on genetic structure, divergence, and diversity in populations of Tetrix subulata pygmy grasshoppers (Orthoptera: Tetrigidae) from 20 sampling locations in southern Sweden. Analyses of 1564 AFLP markers revealed low to moderate levels of genetic diversity (PPL = 59.5–90.1; Hj = 0.23–0.32) within and significant divergence among sampling localities. This suggests that evolution of functional traits in response to divergent selection is possible and that gene flow is restricted. Genetic diversity increased with population size and with increasing proportion of long‐winged phenotypes (a proxy of recent immigration) across populations on the island of Öland, but not on the mainland. Our data further suggested that the open water separating Öland from the mainland acts as a dispersal barrier that restricts migration and leads to genetic divergence among regions. Isolation by distance was evident for short interpopulation distances on the mainland, but gradually disappeared as populations separated by longer distances were included. Results illustrate that integrating ecological and molecular data is key to identifying drivers of population genetic structure in natural populations. Our findings also underscore the importance of landscape structure and spatial sampling scheme for conclusions regarding the role of gene flow and isolation by distance.     

Estimating a geographically explicit model of population divergence

Patterns of genetic variation can provide valuable insights for deciphering the relative roles of different evolutionary processes in species differentiation. However, population-genetic models for studying divergence in geographically structured species are generally lacking. Since these are the biogeographic settings where genetic drift is expected to predominate, not only are population-genetic tests of hypotheses in geographically structured species constrained, but generalizations about the evolutionary processes that promote species divergence may also be potentially biased. Here we estimate a population-divergence model in montane grasshoppers from the sky islands of the Rocky Mountains. Because this region was directly impacted by Pleistocene glaciation, both the displacement into glacial refugia and recolonization of montane habitats may contribute to differentiation. Building on the tradition of using information from the genealogical relationships of alleles to infer the geography of divergence, here the additional consideration of the process of gene-lineage sorting is used to obtain a quantitative estimate of population relationships and historical associations (i.e., a population tree) from the gene trees of five anonymous nuclear loci and one mitochondrial locus in the broadly distributed species Melanoplus oregonensis. Three different approaches are used to estimate a model of population divergence; this comparison allows us to evaluate specific methodological assumptions that influence the estimated history of divergence. A model of population divergence was identified that significantly fits the data better compared to the other approaches, based on per-site likelihood scores of the multiple loci, and that provides clues about how divergence proceeded in M. oregonensis during the dynamic Pleistocene. Unlike the approaches that either considered only the most recent coalescence (i.e., information from a single individual per population) or did not consider the pattern of coalescence in the gene genealogies, the population-divergence model that best fits the data was estimated by considering the pattern of gene lineage coalescence across multiple individuals, as well as loci. These results indicate that sampling of multiple individuals per population is critical to obtaining an accurate estimate of the history of divergence so that the signal of common ancestry can be separated from the confounding influence of gene flow-even though estimates suggest that gene flow is not a predominant factor structuring patterns of genetic variation across these sky island populations. They also suggest that the gene genealogies contain information about population relationships, despite the lack of complete sorting of gene lineages. What emerges from the analyses is a model of population divergence that incorporates both contemporary distributions and historical associations, and shows a latitudinal and regional structuring of populations reminiscent of population displacements into multiple glacial refugia. Because the population-divergence model itself is built upon the specific events shaping the history of M. oregonensis, it provides a framework for estimating additional population-genetic parameters relevant to understanding the processes governing differentiation in geographically structured species and avoids the problems of relying on overly simplified and inaccurate divergence models. The utility of these approaches, as well as the caveats and future improvements, for estimating population relationships and historical associations relevant to genetic analyses of geographically structured species are discussed.     

Population differentiation decreases with depth in deep-sea bivalves

The deep sea is the largest ecosystem on Earth. Recent exploration has revealed that it supports a highly diverse and endemic benthic invertebrate fauna, yet the evolutionary processes that generate this remarkable species richness are virtually unknown. Environmental heterogeneity, topographic complexity, and morphological divergence all tend to decrease with depth, suggesting that the potential for population differentiation may decrease with depth. To test this hypothesis, we use mitochondrial DNA (16S rRNA gene) to examine patterns of population differentiation in four species of protobranch bivalves (Nuculoma similis, Deminucula atacellana, Malletia abyssorum, and Ledella ultima) distributed along a depth gradient in the western North Atlantic. We sequenced 268 individuals from formalin-fixed samples and found 45 haplotypes. The level of sequence divergence among haplotypes within species was similar, but shifted from between populations at bathyal depths to within populations at abyssal depths. Levels of population structure as measured by phiST were considerably greater in the upper bathyal species (N. similis = 0.755 and D. atacellana = 0.931; 530-3834 m) than in the lower bathyal/abyssal species (M. abyssorum = 0.071 and L. ultima = 0.045; 2864-4970 m). Pairwise genetic distances among the samples within each species also decreased with depth. Population trees (UPGMA) based on modified coancestry coefficients and nested clade analysis both indicated strong population-level divergence in the two upper bathyal species but little for the deeper species. The population genetic structure in these protobranch bivalves parallels depth-related morphological divergence observed in deep-sea gastropods. The higher level of genetic and morphological divergence, coupled with the strong biotic and abiotic heterogeneity at bathyal depths, suggests this region may be an active area of species formation. We suggest that the steep, topographically complex, and dynamic bathyal zone, which stretches as a narrow band along continental margins, plays a more important role in the evolutionary radiation of the deep-sea fauna than the much more extensive abyss.