Large discrepancies in differentiation of allozymes, nuclear and mitochondrial DNA loci in recently founded Pacific populations of the pearl oyster Pinctada margaritifera

This study presents a comparative analysis of population structure applied to the pearl oyster (Pinctada margaritifera) from the Central Pacific islands using three classes of molecular markers: two mitochondrial genes (mtDNA), five anonymous nuclear loci (anDNA), and eight polymorphic allozymes. Very low levels of haplotype diversity and nucleotidic divergence detected for mtDNA validate the hypothesis of a recent (re)colonization of Polynesian lagoons after their exondation during the last glaciations. Some nuclear loci, however, showed highly significant FST values, indicating a reduced amount of larval exchange between archipelagos at present. A large interlocus variance of FST was nevertheless observed. We discuss whether this pattern is inherent to the stochasticity of the drift process since recolonization, or if it could result from balancing selection acting on certain loci. This study illustrates once more the need to combine the analysis of several kinds of loci when unrelated phenomena are likely to leave their footprints on genetic structure.     

Mitochondrial DNA under siege in avian phylogeography

Mitochondrial DNA (mtDNA) has been the workhorse of research in phylogeography for almost two decades. However, concerns with basing evolutionary interpretations on mtDNA results alone have been voiced since the inception of such studies. Recently, some authors have suggested that the potential problems with mtDNA are so great that inferences about population structure and species limits are unwarranted unless corroborated by other evidence, usually in the form of nuclear gene data. Here we review the relative merits of mitochondrial and nuclear phylogeographical studies, using birds as an exemplar class of organisms. A review of population demographic and genetic theory indicates that mitochondrial and nuclear phylogeographical results ought to concur for both geographically unstructured populations and for populations that have long histories of isolation. However, a relatively common occurrence will be shallow, but geographically structured mtDNA trees--without nuclear gene corroboration--for populations with relatively shorter periods of isolation. This is expected because of the longer coalescence times of nuclear genes (approximately four times that of mtDNA); such cases do not contradict the mtDNA inference of recent isolation and evolutionary divergence. Rather, the nuclear markers are more lagging indicators of changes in population structure. A review of the recent literature on birds reveals the existence of relatively few cases in which nuclear markers contradict mitochondrial markers in a fashion not consistent with coalescent theory. Preliminary information from nuclear genes suggests that mtDNA patterns will prove to be robust indicators of patterns of population history and species limits. At equilibrium, mitochondrial loci are generally a more sensitive indicator of population structure than are nuclear loci, and mitochondrial estimates of F(ST)-like statistics are generally expected to exceed nuclear ones. Hence, invoking behavioural or ecological explanations of such differences is not parsimonious. Nuclear genes will prove important for quantitative estimates of the depths of haplotype trees, rates of population growth and values of gene flow.     

Variation in mitochondrial DNA reveals multiple distant glacial refugia in black spruce (Picea mariana), a transcontinental North American conifer

Range-wide genetic variation of black spruce (Picea mariana) was studied using polymerase chain reaction-random fragment length polymorphism markers of the mitochondrial genome. Four polymorphic mitochondrial DNA (mtDNA) loci were surveyed and two or three alleles were detected at each locus, resulting in 10 multilocus mtDNA types or mitotypes. A significant subdivision of population genetic diversity was detected (GST = 0.671; NST = 0.726), suggesting low levels of gene flow among populations. The distribution of mitotypes was not random (NST > GST; P < 0.05) and revealed four partially overlapping zones, presumably representative of different glacial populations. Comparison of the genetic structure derived from mtDNA markers and the colonization paths previously deduced from the fossil and pollen records allow us to infer at least three southern and one northeastern glacial populations for black spruce. The patterns revealed in this study suggest that black spruce shares its biogeographical history with other forest-associated North American species.     

Genetic similarity of disjunct populations of the giant sea bass Stereolepis gigas

Sequence variation in nuclear and mitochondrial genes of the giant sea bass Stereolepis gigas collected from the Pacific coast and the northern Sea of Cortez was examined. Restriction fragment length polymorphism analysis and direct sequencing showed extremely low mtDNA sequence diversity (13 closely related haplotypes with no evidence of geographical population subdivision). The mitochondrial haplotype mismatch distribution is consistent with a population expansion following the Last Pleistocene glaciation. Differences in single nucleotide polymorphism frequencies between Pacific and Sea of Cortez populations were detected at two of four nuclear loci, which may reflect natural selection or genetic drift in populations with low effective numbers of males. Although Pacific coast and Sea of Cortez populations of giant sea bass do not exhibit the mitochondrial phylogenetic break characteristic of many species with disjunct Pacific and Gulf populations, the possibility of genetic differentiation at nuclear loci suggests that a cautious approach to broodstock selection for captive breeding and restoration programmes be exercised.     

Microsatellite analysis supports mitochondrial phylogeography of the hellbender (Cryptobranchus alleganiensis)

We investigated genetic diversity of the hellbender (Cryptobranchus alleganiensis) throughout its range in the eastern US using nuclear markers and compared our results to a previously published mitochondrial analysis. A variety of nuclear markers, including protein-coding gene introns and microsatellites were tested but only microsatellites were variable enough for population level analysis. Microsatellite loci showed moderate among population sharing of alleles, in contrast to the reciprocal monophyly exhibited by mitochondrial DNA. However, analyses using F-statistics and Bayesian clustering algorithms showed considerable population subdivision and clustered hellbender populations into the same major groups as the mtDNA. The microsatellites combined with the mtDNA data suggest that gene flow is severely restricted or non-existent among eight major groups, and potentially among populations (rivers) within groups. The combined mtDNA and microsatellite data suggest that the currently recognized hellbender subspecies are paraphyletic. We suggest that the eight independent groups identified in our study should be managed as such, rather than basing conservation decisions on the two named subspecies of hellbender.     

Assessing genetic structure with multiple classes of molecular markers: a case study involving the introduced fire ant Solenopsis invicta.

We used 30 genetic markers of 6 different classes to describe hierarchical genetic structure in introduced populations of the fire ant Solenopsis invicta. These included four classes of presumably neutral nuclear loci (allozymes, codominant random amplified polymorphic DNAs (RAPDs), microsatellites, and dominant RAPDs), a class comprising two linked protein-coding nuclear loci under selection, and a marker of the mitochondrial DNA (mtDNA). Patterns of structure revealed by F statistics and exact tests of differentiation were highly concordant among the four classes of neutral nuclear markers, although the microsatellites were the most effective markers for detecting structure. The results from the mtDNA complemented those from the neutral nuclear markers by revealing that strong limitations to female-mediated gene flow were the cause of the local structure registered by the nuclear markers. The pattern of structure inferred from the selected nuclear loci was markedly different from the patterns derived from the other sets of markers but was predictable on the basis of the presumed mode of selection acting on these loci. In general, the results for all six classes of markers can be explained by known features of the social and reproductive biology of fire ants. Thus, the results from these diverse sets of markers, combined with detailed natural history data, provide an unusually complete picture of how the fundamental evolutionary forces of gene flow, drift, and selection govern the distribution of genetic variation within and between fire ant populations.     

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.     

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.     

Hierarchical Analysis of Genetic Structure in Native Fire Ant Populations: Results from Three Classes of Molecular Markers

We describe genetic structure at various scales in native populations of the fire ant Solenopsis invicta using two classes of nuclear markers, allozymes and microsatellites, and markers of the mitochondrial genome. Strong structure was found at the nest level in both the monogyne (single queen) and polygyne (multiple queen) social forms using allozymes. Weak but significant microgeographic structure was detected above the nest level in polygyne populations but not in monogyne populations using both classes of nuclear markers. Pronounced mitochondrial DNA (mtDNA) differentiation was evident also at this level in the polygyne form only. These microgeographic patterns are expected because polygyny in ants is associated with restricted local gene flow due mainly to limited vagility of queens. Weak but significant nuclear differentiation was detected between sympatric social forms, and strong mtDNA differentiation also was found at this level. Thus, queens of each form seem unable to establish themselves in nests of the alternate type, and some degree of assortative mating by form may exist as well. Strong differentiation was found between the two study regions using all three sets of markers. Phylogeographic analyses of the mtDNA suggest that recent limitations on gene flow rather than longstanding barriers to dispersal are responsible for this large-scale structure.     

Population structure and genetic diversity of metamorphic and paedomorphic populations of the tiger salamander,Ambystoma tigrinum

Genetic relationships, population subdivision and genetic diversity were estimated from mtDNA and allozyme data for two subspecies of tiger salamander, one of which is obligately metamorphic and the other polymorphic for paedomorphosis (larval reproduction). Far greater genetic differentiation exists between subspecies than within subspecies, suggesting that the subspecies have evolved in allopatry. Values of F_st calculated from both mtDNA and allozymes were greater than 0.400 for each subspecies. Significant population subdivision was detected even on a microgeographic scale. This extensive population subdivision indicates that populations can respond to extremely localized selection pressures. In the case of paedomorphosis, populations in permanent water should evolve paedomorphosis as long as the appropriate genes exist. For both mtDNA and allozymes, comparisons of population structure within the polymorphic subspecies and between polymorphic and metamorphic subspecies reveal no discernible effects of paedomorphosis. However, a comparison of paedomorphic and metamorphic populations of the polymorphic subspecies showed significantly higher mtDNA diversity in paedomorphic populations. The discrepancy between the allozyme and mtDNA results may be due to the lower effective population size of mtDNA compared to autosomal genes.     

Concordant evolution of mitochondrial and nuclear genomes in the wheat pathogen Phaeosphaeria nodorum

We compared patterns of mitochondrial restriction fragment length polymorphism (RFLP) diversity with patterns of nuclear RFLP diversity to investigate the effects of selection, gene flow, and sexual reproduction on the population genetic structure and evolutionary history of the wheat pathogen Phaeosphaeria nodorum. A total of 315 fungal isolates from Texas, Oregon, and Switzerland were analyzed using seven nuclear RFLP probes that hybridized to discrete loci and purified mitochondrial DNA that hybridized to the entire mtDNA genome. Forty-two different mitochondrial haplotypes and 298 different nuclear haplotypes were detected. The two most frequent mtDNA haplotypes were present in every population and represented 32% of all isolates. High levels of gene flow, low levels of population subdivision, no evidence for either host specificity or cyto-nuclear disequilibrium were inferred from the analysis of both genomes. The concordance in estimates of these population genetic parameters from both genomes suggests that the two genomes experienced similar degrees of migration, genetic drift and selection.     

Population structure of the African savannah elephant inferred from mitochondrial control region sequences and nuclear microsatellite loci

Two hundred and thirty-six mitochondrial DNA nucleotide sequences were used in combination with polymorphism at four nuclear microsatellite loci to assess the amount and distribution of genetic variation within and between African savannah elephants. They were sampled from 11 localities in eastern, western and southern Africa. In the total sample, 43 haplotypes were identified and an overall nucleotide diversity of 2.0% was observed. High levels of polymorphism were also observed at the microsatellite loci both at the level of number of alleles and gene diversity. Nine to 14 alleles per locus across populations and 44 alleles in the total sample were found. The gene diversity ranged from 0.51 to 0.72 in the localities studied. An analysis of molecular variance showed significant genetic differentiation between populations within regions and also between regions. The extent of subdivision between populations at the mtDNA control region was approximately twice as high as shown by the microsatellite loci (mtDNA F(ST) = 0.59; microsatellite R(ST) = 0.31). We discuss our results in the light of Pleistocene refugia and attribute the observed pattern to population divergence in allopatry accompanied by a recent population admixture following a recent population expansion.     

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.     

Nuclear DNA microsatellite analysis of genetic diversity and gene flow in the Scandinavian brown bear (Ursus arctos)

In the 1930s, the Scandinavian brown bear was close to extinction due to vigorous extermination programmes in Norway and Sweden. Increased protection of the brown bear in Scandinavia has resulted in the recovery of four subpopulations, which currently contain close to 1000 individuals. Effective conservation and management of the Scandinavian brown bear requires knowledge of the current levels of genetic diversity and gene flow among the four subpopulations. Earlier studies of mitochondrial DNA (mtDNA) diversity revealed extremely low levels of genetic variation, and population structure that grouped the three northern subpopulations in one genetic clade and the southernmost subpopulation in a second highly divergent clade. In this study, we extended the analysis of genetic diversity and gene flow in the Scandinavian brown bear using data from 19 nuclear DNA microsatellite loci. Results from the nuclear loci were strikingly different than the mtDNA results. Genetic diversity levels in the four subpopulations were equivalent to diversity levels in nonbottlenecked populations from North America, and significantly higher than levels in other bottlenecked and isolated brown bear populations. Gene flow levels between subpopulations ranged from low to moderate and were correlated with geographical distance. The substantial difference in results obtained using mtDNA and nuclear DNA markers stresses the importance of collecting data from both types of genetic markers before interpreting data and making recommendations for the conservation and management of natural populations. Based on the results from the mtDNA and nuclear DNA data sets, we propose one evolutionarily significant unit and four management units for the brown bear in Scandinavia.     

New evidence from mitochondrial DNA of a progenitor-derivative species relationship between black spruce and red spruce (Pinaceae)

Mitochondrial DNA (mtDNA) markers were used to assess the genetic diversity in allopatric populations of black spruce (Picea mariana [Mill.] BSP) and red spruce (P. rubens Sarg.). Patterns of mitochondrial haplotypes (mitotypes) were strikingly different between the two species. All mtDNA markers surveyed were polymorphic in black spruce, revealing four different mitotypes and high levels of mtDNA diversity (P(p) = 100%, A = 2.0, H = 0.496). In contrast, populations of red spruce had only two mitotypes and harbored low levels of ggenetic diversity (P(p) = 13.2%, A = 1.1, H = 0.120). When the southernmost allopatric populations of red spruce were considered, only one mitotype was detected. As previously reported for nuclear gene loci, the diversity observed for mtDNA in red spruce was a subset of that found in black spruce. Comparison of present and previously published data supports the hypothesis of a recent progenitor-derivative relationship between these species, red spruce presumably being derived by allopatric speciation of an isolated population of black spruce during the Pleistocene.     

Relative information content of polymorphic microsatellites and mitochondrial DNA for inferring dispersal and population genetic structure in the olive sea snake, Aipysurus laevis

Polymorphic microsatellites are widely considered more powerful for resolving population structure than mitochondrial DNA (mtDNA) markers, particularly for recently diverged lineages or geographically proximate populations. Weaker population subdivision for biparentally inherited nuclear markers than maternally inherited mtDNA may signal male-biased dispersal but can also be attributed to marker-specific evolutionary characteristics and sampling properties. We discriminated between these competing explanations with a population genetic study on olive sea snakes, Aipysurus laevis. A previous mtDNA study revealed strong regional population structure for A. laevis around northern Australia, where Pleistocene sea-level fluctuations have influenced the genetic signatures of shallow-water marine species. Divergences among phylogroups dated to the Late Pleistocene, suggesting recent range expansions by previously isolated matrilines. Fine-scale population structure within regions was, however, poorly resolved for mtDNA. In order to improve estimates of fine-scale genetic divergence and to compare population structure between nuclear and mtDNA, 354 olive sea snakes (previously sequenced for mtDNA) were genotyped for five microsatellite loci. F statistics and Bayesian multilocus genotype clustering analyses found similar regional population structure as mtDNA and, after standardizing microsatellite F statistics for high heterozygosities, regional divergence estimates were quantitatively congruent between marker classes. Over small spatial scales, however, microsatellites recovered almost no genetic structure and standardized F statistics were orders of magnitude smaller than for mtDNA. Three tests for male-biased dispersal were not significant, suggesting that recent demographic expansions to the typically large population sizes of A. laevis have prevented microsatellites from reaching mutation-drift equilibrium and local populations may still be diverging.     

Microsatellite Analysis of Genetic Variation and Population Subdivision for the Black Muntjac, Muntiacus crinifrons

The black muntjac (Muntiacus crinifrons) is a rare deer found only in a restricted region in east China. Recent studies of mitochondrial DNA diversity have shown a markedly low level of nucleotide diversity for the species, and the Suichang population was genetically differentiated from the two other populations, in Huangshan and Tianmushan mountains. In this study, we extended the analysis of genetic diversity and population subdivision for the black muntjac using data from 11 highly polymorphic nuclear DNA microsatellite loci. Contrary to the results based on mtDNA data, the microsatellite loci revealed that the black muntjac retained a rather high nuclear genetic diversity (overall average H (E) = 0.78). Nevertheless, both types of markers supported the idea that the extant black muntjac population is genetically disrupted (overall phi (ST) = 0.16 for mtDNA and overall F (ST) = 0.053 for microsatellite, both P < 0.001). The correlation between genetic differentiation and geographic distance was not significant (Mantel test; P > 0.05), implying that the patterns of genetic differentiation observed in this study might result from recent habitat fragmentation or loss. Based on the results from the mtDNA and nuclear DNA data sets, two management units were defined for the species, Huangshan/Tianmushan and Suichang. We also recommend that a new captive population be established with individuals from the Suichang region as a founder source.     

Phylogeography recapitulates topography: very fine-scale local endemism of a saproxylic 'giant' springtail at Tallaganda in the Great Dividing Range of south-east Australia

Comparative phylogeography can reveal processes and historical events that shape the biodiversity of species and communities. As part of a comparative research program, the phylogeography of a new, endemic Australian genus and species of log-dependent (saproxylic) collembola was investigated using mitochondrial sequences, allozymes and anonymous single-copy nuclear markers. We found the genetic structure of the species corresponds with five a priori microbiogeographical regions, with population subdivision at various depths owing to palaeoclimatic influences. Closely related mtDNA haplotypes are codistributed within a single region or occur in adjacent regions, nuclear allele frequencies are more similar among more proximate populations, and interpopulation migration is rare. Based on mtDNA divergence, a late Miocene-late Pliocene coalescence is likely. The present-day distribution of genetic diversity seems to have been impacted by three major climatic events: Pliocene cooling and drying (2.5-7 million years before present, Mybp), early Pleistocene wet-dry oscillations (c. 1.2 Mybp) and the more recent glacial-interglacial cycles that have characterized the latter part of the Quaternary (<0.4 Mybp).     

Nuclear and cytoplasmic variation within and between Eurasian Larix (Pinaceae) species

The genetic variation in nuclear and cytoplasmic markers was investigated in 28 populations belonging to Eurasian Larix species (L. decidua, L. sibirica, L. gmelinii, L. olgensis, and L. kaempferi). Nuclear genetic variation was assessed at 214 AFLP loci, and both PCR-RFLP and four microsatellite loci were used to estimate variation of the chloroplast DNA. Variation of the mitochondrial genome was measured using RFLPs. Although population differentiation at both nuclear and chloroplast markers was much weaker than at mitochondrial DNA, it nonetheless corroborated the grouping observed with mitochondrial DNA. The AFLPs led to the same population grouping as mtDNA. Notably, the presence of two ancient western and eastern groups within L. sibirica was confirmed and possible postglacial routes inferred. The genetic composition of the northernmost L. sibirica population in our sample established that it is located at the confluence of the eastern and western recolonization routes. The joint use of the three markers also indicated that populations around Lake Baikal are hybrids between L. gmelinii and L. sibirica, with L. gmelinii primarily acting as the pollinator. Finally, AFLP-based estimates of nucleotide variation were an order of magnitude larger than the strikingly low estimates of nucleotide variation recently reported in Pinus sylvestris.     

The Atlantic-Mediterranean transition: discordant genetic patterns in two seabream species, Diplodus puntazzo (Cetti) and Diplodus sargus (L.)

Sparids are a group of demersal perciform fish of high commercial value, which have experienced an extensive radiation, particularly in the Mediterranean, where they occupy a variety of different niches. The present study focuses on two species: Diplodus sargus and D. puntazzo, presenting a wide distribution from the Mediterranean to the eastern Atlantic coasts. They display similar ecological behaviour and are evolutionary closely related. Both are highly appreciated in fisheries and D. puntazzo is currently under domestication process. However, little is know on their population structure and it is an open question whether any genetic differentiation exists at the geographic level. To address this issue we examined sequence variation of a portion of the mitochondrial DNA (mtDNA) control region in population samples of each of the two species collected over a wide geographic range. In addition to the mtDNA, analysis of nuclear loci (allozymes) was included in the study to compare patterns revealed by nuclear and mitochondrial markers. The studied samples covered an area from the eastern Mediterranean to the Portuguese coasts immediately outside the Gibraltar Strait. The two species revealed a level of sequence polymorphism remarkably different for the control region with the D. puntazzo and D. sargus showing 111 and 28 haplotypes, respectively. Such a difference was not detected with allozyme markers. The two species also showed large differences in their population structure. While D. puntazzo presented a marked genetic divergence between the Atlantic and Mediterranean samples, D. sargus showed little intraspecific differentiation. These results were supported using both mtDNA and allozyme markers, and were interpreted as the consequence of differences in the history of the two species such as fluctuations in the effective population size due to bottlenecks and expansions, possibly combined with present-day differences in levels of gene flow.