Concordance of genetic divergence among sockeye salmon populations at allozyme, nuclear DNA, and mitochondrial DNA markers

Abstract.— We examined genetic variation at 21 polymorphic allozyme loci, 15 nuclear DNA loci, and mitochondrial DNA in four spawning populations of sockeye salmon ( Oncorhynchus nerka ) from Cook Inlet, Alaska, to test for differences in the patterns of divergence among different types of markers. We were specifically interested in testing the suggestion that natural selection at allozyme loci compromises the effectiveness of these markers for describing the amount and patterns of gene flow among populations. We found concordance among markers in the amount of genetic variation within and among populations, with the striking exception of one allozyme locus ( sAH ), which exhibited more than three times the amount of among-population differentiation as other loci. A consideration of reports of discordance between allozymes and other loci indicates that these differences usually result from one or two exceptional loci. We conclude that it is important to examine many loci when estimating genetic differentiation to infer historical amounts of gene flow and patterns of genetic exchange among populations. It is less important whether those loci are allozymes or nuclear DNA markers.     

Concordance of allozyme and microsatellite differentiation in a marine fish, but evidence of selection at a microsatellite locus

Previous studies have reported higher levels of divergence for microsatellites than for allozymes in several species, suggested to reflect stabilizing selection on the allozymes. We compared the differentiation patterns of 11 allozyme and nine microsatellite loci using 679 spawning Atlantic herring (Clupea harengus) collected in the Baltic and North Seas to test for differential natural selection on these markers. Observed distributions of F statistics for the two types of markers are conspicuously dissimilar, but we show that these differences can largely be explained by sampling phenomena caused by different allele frequency distributions and degrees of variability. The results show consistently low levels of differentiation for both marker types, with the exception of one outlier microsatellite locus with a notably high F(ST). The aberrant pattern at this locus is primarily due to two alleles occurring at markedly high frequencies in the Baltic, suggesting selection at this locus, or a closely linked one. When excluding this locus, the two marker types show similar, weak differentiation patterns with F(ST) values between the Baltic and the North Seas of 0.001 and 0.002 for allozymes and microsatellites, respectively. This small heterogeneity, and weak isolation by distance, is easier to distinguish statistically with microsatellites than with allozymes that have fewer alleles and skewed frequency distributions. The allozymes, however, also detect surprisingly low levels of divergence. Our results support suggestions that previously described differences between marker types are primarily caused by a small number of outlier loci.     

Allozyme and microsatellite loci provide discordant estimates of population differentiation in the endangered dusky grouper (Epinephelus marginatus) within the Mediterranean Sea

The dusky grouper, Epinephelus marginatus, inhabits coastal reefs in the Mediterranean Sea and Atlantic Ocean. A decline in the abundance of this long-lived protogynous hermaphrodite has led to its listing as an endangered species in the Mediterranean, and heightened management concerns regarding its genetic variability and population substructure. To address these concerns, we analysed genetic variation at seven microsatellite and 28 allozyme loci in dusky groupers sampled from seven areas (for microsatellites) and three areas (for allozymes) in the west-central Mediterranean. Levels of genetic variability were higher for microsatellites than for allozymes (mean H(E) = 0.78 and 0.07, respectively), but similar to those observed in other marine fishes with comparable markers. Both microsatellites and allozymes revealed significant genetic differentiation among all areas analysed with each class of marker, but the magnitude of differentiation revealed by allozymes over three locales (F(ST) = 0.214) was greater than that detected with microsatellites over seven areas, or over the three areas shared with the allozyme analysis (F(ST) = 0.018 and approximately 0, respectively). A large proportion of the allozyme differentiation was due to a single locus (ADA*) possibly influenced by selection, but allozyme differentiation over the three areas was still highly significant (F(ST) = 0.06, P < 0.0001), and the 95% confidence intervals for allozyme and microsatellite F(ST) did not overlap when this locus was excluded. There was no evidence of isolation by distance with either class of markers. Our results lead us to conclude that dusky groupers are not panmictic in the Mediterranean Sea and suggest that they should be managed on a local basis. However, more work is needed to elucidate genetic relationships among populations.     

Genetic differentiation and history of populations of the Italian treefrog Hyla intermedia: lack of concordance between mitochondrial and nuclear markers

The genetic differentiation among 33 populations of the Italian treefrog, Hyla intermedia (Anura: Hylidae), was investigated using both biparentally (23 allozyme loci) and maternally (partial mitochondrial cytochrome b gene) inherited markers. Two main population groups were evidenced by both markers, located north and south of the northern Apennines. However, the pattern of differentiation between these two groups was much less pronounced at allozymes than at mtDNA, leading to gene flow estimates that were 25 times lower at mitochondrial than at nuclear level. Also, the mtDNA divergence between the two groups was particularly marked for two cospecific lineages of anuran amphibians (the P-distance being on average 9.04%), while their average genetic distance at allozymes was comparatively low (D _NEI = 0.07). This contrasting pattern of nuclear versus mitochondrial genetic variation is discussed in the context of: (1) marker specific selection, (2) secondary contact and sex-biased gene flow and (3) ancestral polymorphism and colonization from north to south. Finally we emphasize how, for population genetic studies, the use of multiple markers having distinct evolutionary properties can help unravel the existence of more complex evolutionary histories.     

Genetic variation at allozyme and RAPD loci in sessile oak Quercus petraea (Matt.) Liebl.: the role of history and geography

The nuclear genetic variation within and among 21 populations of sessile oak was estimated at 31 RAPD loci in conjunction with previous estimates of variation at eight allozyme loci. The aim of the study was to assess the relative role of isolation-by-distance and postglacial history on patterns of nuclear variation. Because of its small effective population size and maternal transmission, the chloroplast genome is a good marker of population history. Both kinds of nuclear variation (RAPD and allozyme) were therefore compared, first, to the geographical distances among populations and, secondly, to chloroplast DNA restriction polymorphism in the same populations. Multiple Mantel tests were used for this purpose. Although RAPDs revealed less genetic diversity than allozymes, levels of genetic differentiation ( G _ST) were identical. The standard genetic distance calculated at all RAPD loci was correlated with geographical distances but not with the genetic distance calculated from chloroplast DNA data. Conversely, allozyme variation was correlated with chloroplast DNA variation, but not with geography. Possibly, divergent selection at two allozyme loci during the glacial period could explain this pattern. Because of its greater number of loci assayed, RAPDs probably provided a less biased picture of the relative role of geography and history.     

Allozyme electrophoresis still represents a powerful technique in the management of coral reefs

Understanding genetic variability and gene flow between populations of scleractinian corals separated by one to several hundred kilometers is crucially important as we head into a century of climate change in which an understanding of the connectivity of populations is a critically important question in management. Genetic methods that directly use molecular variance in the DNA should offer greater precision in detecting differences among individuals and populations than the more traditional allozyme electrophoresis. However, this paper highlights the point that the limited number of DNA markers that have been identified for scleractinian coral genetic studies do not necessarily offer greater precision than that offered by allozymes. In fact, at present allozyme electrophoresis yields greater information than the eight different DNA markers used in this study. Given the relative ease of use of allozymes and the wealth of comparable data sets from numerous previously published studies, allozyme electrophoresis should not be dismissed for population structure and connectivity studies on coral reefs. While continued effort should be placed into searching for new DNA markers, until a more sensitive DNA marker becomes available for scleractinian corals, allozyme electrophoresis remains a powerful and relevant technique for understanding the connectivity of coral population studies.     

Amplified fragment length polymorphism (AFLP) analysis of the genetic structure of the zebra mussel, Dreissena polymorpha, in the Mississippi River

1. We predicted that zebra mussel, Dreissena polymorpha (Pallas), genetic structure in the Mississippi River would follow a model of invasive species genetics, which predicts low genetic structure among populations of recently established species. This prediction was upheld in our previous genetic study using allozymes, however, one locus yielded anomalous results. 2. We employed amplified fragment length polymorphism (AFLP) analysis as a neutral marker to assess the amount of genetic structure within and among populations, and as a test of expected population structure from both invasion genetic theory, and the results from our previous study. 3. There was greater spatial differentiation, as measured by F_st, observed using AFLP's than for allozymes (P < 0.001). There was no evidence that AFLP variation conformed to an isolation by distance model, and genetic relationships of populations, as measured by AFLP markers, were not similar to those detected in our allozyme survey. 4. The lack of concordance between these two genetic marker systems probably reflects their differential responses to drift, migration, and selection occurring during this rapid invasion. Strong population structure is counter to predictions that populations of invasive species will not be differentiated, as with observations based on allozyme markers. Therefore, newly established species may require genetic surveys using multiple marker systems to evaluate population structure.     

Introgression between Two Cutthroat Trout Subspecies with Substantial Karyotypic, Nuclear and Mitochondrial Genomic Divergence

The authors used allozymes encoded by nuclear genes and restriction enzyme analysis of mitochondrial DNA (mtDNA) to study secondary contact between westslope (Salmo clarki lewisi) and Yellowstone cutthroat trout (Salmo clarki bouvieri) in Forest Lake, Montana. Eleven diagnostic allozyme loci identified this as a random-mating hybrid swarm. No parental, first-generation hybrid or backcross genotypes were detected in the sample (N = 33), and genotype distributions at all the variable loci conform to binomial expectations. There is little linkage disequilibrium between the diagnostic loci, indicating that the nuclear genomes of the two subspecies are largely randomly associated. The allozymes and mtDNA give identical estimates of the proportional genetic contribution of each subspecies. Thus, males and females from both subspecies have contributed equally to this hybrid swarm. Although these subspecies have accumulated substantial genetic divergence between their nuclear (Nei's D = 0.34) and mitochondrial (2% sequence divergence) genomes, this has not resulted in a genetic barrier to exchange between them.     

Comparative allozyme and microsatellite population structure in a narrow endemic plant species, Centaurea corymbosa Pourret (Asteraceae)

Centaurea corymbosa Pourret (Asteraceae) is a narrow endemic species known only from six populations located in a 3-km2 area in the south of France. Earlier field experiments have suggested that pollen and seed dispersal were highly restricted within and among populations. Consistent with the field results, populations were highly differentiated for five allozyme loci and among-population variation fitted an isolation-by-distance model. In the present study, we investigated the genetic structure of C. corymbosa using six microsatellite loci. As with allozymes, microsatellites revealed no within-population structure and a large differentiation among populations. However, allozyme loci were less powerful than microsatellites in detecting the extent of gene flow assessed by assignment tests. The patterns of structuration greatly varied among loci for both types of marker; we suggest that differences in single-locus pattern could mainly be an effect of stochastic variation for allozymes and an effect of variation in mutation rate for microsatellites. In contrast to the multilocus results, the two most polymorphic microsatellite loci did not show any isolation-by-distance pattern. Our results suggest that highly variable loci might not always be the best suited markers to quantify levels of gene flow among populations.     

Distinguishing adaptive from nonadaptive genetic differentiation: comparison of Q(ST) and F(ST) at two spatial scales

Genetic differentiation in 20 hierarchically sampled populations of wild barley was analyzed with quantitative traits, allozymes and Random Amplified Polymorphic DNAs (RAPDs), and compared for three marker types at two hierarchical levels. Regional subdivision for both molecular markers was much lower than for quantitative traits. For both allozymes and RAPDs, most loci exhibited minor or no regional differentiation, and the relatively high overall estimates of the latter were due to several loci with exceptionally high regional differentiation. The allozyme- and RAPD-specific patterns of differentiation were concordant in general with one another, but not with quantitative trait differentiation. Divergent selection on quantitative traits inferred from very high regional Q(ST) was in full agreement with our previous results obtained from a test of local adaptation and multilevel selection analysis. In contrast, most variation in allozyme and RAPD variation was neutral, although several allozyme loci and RAPD markers were exceptional in their levels of regional differentiation. However, it is not possible to answer the question whether these exceptional loci are directly involved in the response to selection pressure or merely linked to the selected loci. The fact that Q(ST) and F(ST) did not differ at the population scale, that is, within regions, but differed at the regional scale, for which local adaptation has been previously shown, implies that comparison of the level of subdivision in quantitative traits, as compared with molecular markers, is indicative of adaptive population differentiation only when sampling is carried out at the appropriate scale.     

Do discrepancies between microsatellite and allozyme variation reveal differential selection between sea and lagoon in the sea bass (Dicentrarchus labrax)?

In the present study the genetic structure of Dicentrarchus labrax (14 samples from the Mediterranean) was analysed at six microsatellite loci, in order to test the hypothesis that some enzymatic loci undergo selection between marine and lagoon habitat. Eight of the 14 samples were analysed at both microsatellite and allozyme markers. The analysis of the genetic variation among the Mediterranean samples showed that (i) &Fcirc;ST values obtained with the six microsatellite loci were much smaller than those obtained with the 28 allozymes and (ii) microsatellite loci seemed to reflect more the geographical proximity than an ecological one. Thirteen enzymatic loci exhibited moderate to high values compared with microsatellites. This was interpreted as evidence that these allozymes are non-neutral. However, only six loci seemed to be implicated in differentiation between marine and lagoon samples, the causes of selection being unknown for the others. A possible scenario of population dynamics of the sea bass between marine and lagoon habitat is suggested.     

Genetic differentiation and local adaptation in the salt-marsh beetle Pogonus chalceus: a comparison between allozyme and microsatellite loci

The genetic structure of Pogonus chalceus from 11 Atlantic and seven Mediterranean Western European populations was analysed. Results from seven allozyme loci were compared to five microsatellites to test the hypothesis that some enzymatic loci undergo selection. Both allozyme and microsatellite results showed that Mediterranean beetle populations are genetically distinct from Atlantic populations. The analysis of the genetic structure showed that FST values derived from all microsatellite loci were much smaller than those obtained from allozymes. The enzymatic locus Idh-1 exhibited a high value compared to the other loci, suggesting that it is non-neutral. The same Idh-1 locus was implicated in differentiation between temporary and stable populations, as followed also from a highly significant relationship between the allele two of this Idh-1 locus and dispersal power population estimates. The 'parallel evolution' model may account for the diversification of locally adapted Pogonus chalceus populations between different microhabitats.     

Allozymes and nDNA markers show different levels of population differentiation in the mussel Mytilus?edulis on British coasts

In order to test the hypothesis that natural selection can be operating on allozymes, six nuclear DNA markers were used to study population structure in six Mytilus edulis (Linnaeus, 1758) populations around the British Isles and the results compared with previous allozyme data from the same (two sites) or very close localities (the other four sites). Allozymes and nDNA markers show significant differentiation between populations at the P < 0.01 level, with the exception of the nDNA marker PLIIa. The results also indicate significant differences (modified version of the Lewontin–Krakauer test) between the two sets of markers which are about one order of magnitude higher for nDNA than for the allozyme loci. The bootstrap analysis was carried out to correct for the possibility that this difference is a sample size artefact. These results are consistent with the operation of balancing selection on allozymes or on loci in linkage disequilibrium with them. Handling editor: C. Sturmbauer     

Reconciling patterns of inter-ocean molecular variance from four classes of molecular markers in blue marlin (Makaira nigricans)

Different classes of molecular markers occasionally yield discordant views of population structure within a species. Here, we examine the distribution of molecular variance from 14 polymorphic loci comprising four classes of molecular markers within approximately 400 blue marlin individuals (Makaira nigricans). Samples were collected from the Atlantic and Pacific Oceans over 5 years. Data from five hypervariable tetranucleotide microsatellite loci and restriction fragment length polymorphism (RFLP) analysis of whole molecule mitochondrial DNA (mtDNA) were reported and compared with previous analyses of allozyme and single-copy nuclear DNA (scnDNA) loci. Temporal variance in allele frequencies was nonsignificant in nearly all cases. Mitochondrial and microsatellite loci revealed striking phylogeographic partitioning among Atlantic and Pacific Ocean samples. A large cluster of alleles was present almost exclusively in Atlantic individuals at one microsatellite locus and for mtDNA, suggesting that, if gene flow occurs, it is likely to be unidirectional from Pacific to Atlantic oceans. Mitochondrial DNA inter-ocean divergence (FST) was almost four times greater than microsatellite or combined nuclear divergences including allozyme and scnDNA markers. Estimates of Neu varied by five orders of magnitude among marker classes. Using mathematical and computer simulation approaches, we show that substantially different distributions of FST are expected from marker classes that differ in mode of inheritance and rate of mutation, without influence of natural selection or sex-biased dispersal. Furthermore, divergent FST values can be reconciled by quantifying the balance between genetic drift, mutation and migration. These results illustrate the usefulness of a mitochondrial analysis of population history, and relative precision of nuclear estimates of gene flow based on a mean of several loci.     

Selection on outlier loci and their association with adaptive phenotypes in Littorina saxatilis contact zones

A fundamental issue in speciation research is to evaluate phenotypic variation and the genomics driving the evolution of reproductive isolation between sister taxa. Above all, hybrid zones are excellent study systems for researchers to examine the association of genetic differentiation, phenotypic variation and the strength of selection. We investigated two contact zones in the marine gastropod Littorina saxatilis and utilized landmark‐based geometric morphometric analysis together with amplified fragment length polymorphism (AFLP) markers to assess phenotypic and genomic divergence between ecotypes under divergent selection. From genetic markers, we calculated the cline width, linkage disequilibrium and the average effective selection on a locus. Additionally, we conducted an association analysis linking the outlier loci and phenotypic variation between ecotypes and show that a proportion of outlier loci are associated with key adaptive phenotypic traits.     

Genotypic selection in Daphnia populations consisting of inbred sibships

The genetic basis of fitness reduction associated with inbreeding is still poorly understood. Here we use associations between allozyme genotypes and fitness to investigate the genetic basis of inbreeding depression in experimental outdoor populations of the water flea, Daphnia magna. In Daphnia, a phase of clonal reproduction follows hatching from sexually produced resting eggs, and changes in genotype frequencies during the clonal phase can be used to estimate fitness. Our experiment resembles natural colonization of ponds in that single clones colonize an empty pool, expand asexually and produce sexual offspring by selfing (sisters mate with their clonal brothers). These offspring diapause and form populations consisting of selfed sibships in the following spring. In 12 of 13 experimental populations, genotypes of selfed hatchlings after diapause conformed to Mendelian expectations. During the subsequent ca. 10 asexual generations, however, genotype frequencies changed significantly at 19 of 27 single loci studied within populations, mostly in favour of heterozygotes, with heterozygosity at multiple loci affecting the change in genotype frequency multiplicatively. Because variance in heterozygosity among siblings at a given marker reflects only heterozygosity in the chromosomal region around this marker, our results suggest that selection at fitness-associated loci in the chromosomal regions near the markers were responsible for these changes. The genotype frequency changes were more consistent with selection acting on linked loci than on the allozymes themselves. Taken together, the evidence for abundant selection in the chromosomal regions of the markers and the fact that changes in genotype frequencies became apparent only after several generations of clonal selection, point to a genetic load consisting of many alleles of small or intermediate effects, which is consistent with the strong genetic differentiation and repeated genetic bottlenecks in the metapopulation from which the animals for this study were obtained.     

GEOGRAPHIC DISTRIBUTION OF MOLECULAR VARIANCE WITHIN THE BLUE MARLIN (MAKAIRA NIGRICANS): A HIERARCHICAL ANALYSIS OF ALLOZYME,SINGLE-COPY NUCLEAR DNA,AND MITOCHONDRIAL DNA MARKERS

This study presents a comparative hierarchical analysis of variance applied to three classes of molecular markers within the blue marlin (Makaira nigricans). Results are reported from analyses of four polymorphic allozyme loci, four polymorphic anonymously chosen single-copy nuclear DNA (scnDNA) loci, and previously reported restriction fragment length polymorphisms (RFLPs) of mitochondrial DNA (mtDNA). Samples were collected within and among the Atlantic and Pacific Oceans over a period of several years. Although moderate levels of genetic variation were detected at both polymorphic allozyme (H = 0.30) and scnDNA loci (H = 0.37), mtDNA markers were much more diverse (h = 0.85). Allele frequencies were significantly different between Atlantic and Pacific Ocean samples at three of four allozyme loci and three of four scnDNA loci. Estimates of allozyme genetic differentiation (θ_O) ranged from 0.00 to 0.15, with a mean of 0.08. The θ_O values for scnDNA loci were similar to those of allozymes, ranging from 0.00 to 0.12 with a mean of 0.09. MtDNA RFLP divergence between oceans (θ_O = 0.39) was significantly greater than divergence detected at nuclear loci (95% nuclear confidence interval = 0.04–0.11). The fourfold smaller effective population size of mtDNA and male-mediated gene flow may account for the difference observed between nuclear and mitochondrial divergence estimates.     

Lack of geographic variation in anonymous nuclear polymorphisms in the American oyster, Crassostrea virginica

Comparing geographic variation of noncoding nuclear DNA polymorphisms, which presumably are neutral to natural selection, with geographic variation of allozymes is potentially a good way to detect the effects of selection on allozyme polymorphisms. A previous study of four anonymous nuclear markers in the American oyster, Crassostrea virginica, found dramatic differences in allele frequency between the Gulf of Mexico and the Atlantic Ocean. In contrast, 14 allozyme polymorphisms were fairly uniform in frequency between the two areas. This led to the conclusion that all of the allozyme polymorphisms were kept uniform in frequency by balancing selection. To test the robustness of this pattern, six additional anonymous nuclear DNA polymorphisms were surveyed in oysters from Panacea, Fla, and Charleston, S.C. on the Gulf and Atlantic coasts, respectively. Unlike the previously studied DNA markers, the six DNA polymorphisms examined here show geographic variation that is not significantly greater than that of allozymes. The reason for the discrepancy between the two sets of DNA polymorphisms is unclear.      

Genetic architecture of sexual and asexual populations of the aphid Rhopalosiphum padi based on allozyme and microsatellite markers

Cyclical parthenogens, including aphids, are attractive models for comparing the genetic outcomes of sexual and asexual reproduction, which determine their respective evolutionary advantages. In this study, we examined how reproductive mode shapes genetic structure of sexual (cyclically parthenogenetic) and asexual (obligately parthenogenetic) populations of the aphid Rhopalosiphum padi by comparing microsatellite and allozyme data sets. Allozymes showed little polymorphism, confirming earlier studies with these markers. In contrast, microsatellite loci were highly polymorphic and showed patterns very discordant from allozyme loci. In particular, microsatellites revealed strong heterozygote excess in asexual populations, whereas allozymes showed heterozygote deficits. Various hypotheses are explored that could account for the conflicting results of these two types of genetic markers. A strong differentiation between reproductive modes was found with both types of markers. Microsatellites indicated that sexual populations have high allelic polymorphism and heterozygote deficits (possibly because of population subdivision, inbreeding or selection). Little geographical differentiation was found among sexual populations confirming the large dispersal ability of this aphid. In contrast, asexual populations showed less allelic polymorphism but high heterozygosity at most loci. Two alternative hypotheses are proposed to explain this heterozygosity excess: allele sequence divergence during long-term asexuality or hybrid origin of asexual lineages. Clonal diversity of asexual lineages of R. padi was substantial suggesting that they could have frozen genetic diversity from the pool of sexual lineages. Several widespread asexual genotypes were found to persist through time, as already seen in other aphid species, a feature seemingly consistent with the general-purpose genotype hypothesis.     

Parallel microgeographic patterns of genetic diversity and divergence revealed by allozyme, RAPD, and microsatellites in Triticum dicoccoides at Ammiad, Israel

The levels of genetic diversity were compared by means of 35 allozyme, 60 RAPD, and 25 microsatellite (SSR) markers for 75–175 individuals of tetraploid wild emmer wheat (Triticum dicoccoides) collected in 1993 from a microgeographic microsite, Ammiad, north of the Sea of Galilee, Israel. This microsite included four major habitats, which showed highly significant differentiation in ecological factors, in particular with respect to rock cover, proximity and height, and surface soil moisture in the early growing season of T. dicoccoides. Higher within-subpopulation genetic diversity was found in the primarily non-coding DNA regions (RAPD and SSR) rather than in the protein-coding (allozymes) regions. However, much larger gene differentiation (G _ST) among the subpopulations was observed in the protein-coding allozymes than in the RAPDs and SSRs. Larger genetic distance was found at SSR loci, followed by allozyme and RAPD loci. The subpopulations in drier habitats tend to have higher allozyme, RAPD and SSR diversities (He), the relatively wet Karst subpopulation showed only about half He of the other relatively drier habitats. The subpopulations with larger difference of soil moisture between habitats tend to show larger genetic distances at allozyme, RAPD and SSR loci. These results suggest that climatic selection through aridity stress may be an important factor acting on both structural protein-coding and presumably partly regulatory non-coding DNA regions, resulting in microscale adaptive patterns, although hitchhiking and random drift may also intervene. These results have profound implications for genetic conservation both in situ and ex situ.