High connectivity on a global scale in the pelagic wahoo,Acanthocybium solandri (tuna family Scombridae)

The population genetic structure and phylogeography of wahoo, Acanthocybium solandri, were investigated on a global scale with intron six of lactate dehydrogenase‐A (ldhA6, 8 locations, N = 213) and mtDNA cytochrome b (Cytb, 10 locations, N = 322). Results show extensive sharing of haplotypes across the wahoo's entire global range, and analyses were unable to detect significant structure (nuclear F_ST = 0.0125, P = 0.106; mtDNA Φ_ST < 0.0001, P = 0.634). Power analyses indicated 95% confidence in detecting nuclear F_ST ≥ 0.0389 and mtDNA Φ_ST ≥ 0.0148. These findings appear unique, as most other tunas, billfishes, and oceanic sharks exhibit significant population structure on the scale of East–West Atlantic, Atlantic vs. Indian‐Pacific, or East–West Pacific. Overall nuclear heterozygosity (H = 0.714) and mtDNA haplotype diversity (h = 0.918) are both high in wahoo, while overall mtDNA nucleotide diversity (π= 0.006) and nuclear nucleotide diversity (π=0.004) are uniformly low, indicating a recent increase in population size. Coalescence analyses yield an estimate of effective female population size (N_eF) at ~816 000, and a population bottleneck ~690 000 years ago. However, conclusions about population history from our Cytb data set are not concordant with a control region survey, a finding that will require further investigation. This is the first example of a vertebrate with a single globally distributed population, a finding we attribute to extensive dispersal at all life stages. The indications of a worldwide stock for wahoo reinforce the mandate for international cooperation on fisheries issues.     

Genetic connectivity among swarming sites in the wide ranging and recently declining little brown bat (Myotis lucifugus)

Characterizing movement dynamics and spatial aspects of gene flow within a species permits inference on population structuring. As patterns of structuring are products of historical and current demographics and gene flow, assessment of structure through time can yield an understanding of evolutionary dynamics acting on populations that are necessary to inform management. Recent dramatic population declines in hibernating bats in eastern North America from white‐nose syndrome have prompted the need for information on movement dynamics for multiple bat species. We characterized population genetic structure of the little brown bat, Myotis lucifugus, at swarming sites in southeastern Canada using 9 nuclear microsatellites and a 292‐bp region of the mitochondrial genome. Analyses of F_ST, Φ_ST, and Bayesian clustering (STRUCTURE) found weak levels of genetic structure among swarming sites for the nuclear and mitochondrial genome (Global F_ST = 0.001, P < 0.05, Global Φ_ST = 0.045, P < 0.01, STRUCTURE K = 1) suggesting high contemporary gene flow. Hierarchical AMOVA also suggests little structuring at a regional (provincial) level. Metrics of nuclear genetic structure were not found to differ between males and females suggesting weak asymmetries in gene flow between the sexes. However, a greater degree of mitochondrial structuring does support male‐biased dispersal long term. Demographic analyses were consistent with past population growth and suggest a population expansion occurred from approximately 1250 to 12,500 BP, following Pleistocene deglaciation in the region. Our study suggests high gene flow and thus a high degree of connectivity among bats that visit swarming sites whereby mainland areas of the region may be best considered as one large gene pool for management and conservation.     

Isolation and characterization of polymorphic microsatellite loci for the study of paternity and population structure in the little penguin Eudyptula minor

We isolated and characterized eight novel microsatellite loci in the little penguin Eudyptula minor, using nonradioactive polymerase chain reaction‐based techniques to screen GA and GAAA repeats from enriched genomic DNA libraries. All eight loci were polymorphic and seven were variable in our main study population (mean H_E = 0.613, mean N_A = 7.14). Cross‐amplification using a microsatellite primer developed in Spheniscus demersus (African penguin) yielded one additional polymorphic locus. This locus combined with six of the little penguin loci is suitable for paternity assignment in little penguins (exclusion probability for seven unlinked loci = 0.993).     

Phylogeography of Pinus tabulaeformis Carr. (Pinaceae), a dominant species of coniferous forest in northern China

How coniferous trees in northern China changed their distribution ranges in response to Quaternary climatic oscillations remains largely unknown. Here we report a study of the phylogeography of Pinus tabulaeformis, an endemic and dominant species of coniferous forest in northern China. We examined sequence variation of maternally inherited, seed‐dispersed mitochondrial DNA (mtDNA) (nad5 intron 1 and nad4/3–4) and paternally inherited, pollen‐ and seed‐dispersed chloroplast DNA (cpDNA) (rpl16 and trnS‐trnG) within and among 30 natural populations across the entire range of the species. Six mitotypes and five chlorotypes were recovered among 291 trees surveyed. Population divergence was high for mtDNA variation (G_ST = 0.738, N_ST = 0.771) indicating low levels of seed‐based gene flow and significant phylogeographical structure (N_ST > G_ST, P < 0.05). The spatial distribution of mitotypes suggests that five distinct population groups exist in the species: one in the west comprising seven populations, a second with a north–central distribution comprising 15 populations, a third with a southern and easterly distribution comprising five populations, a fourth comprising one central and one western population, and a fifth comprising a single population located in the north‐central part of the species’ range. Each group apart from the fourth group is characterized by a distinct mitotype, with other mitotypes, if present, occurring at low frequency. It is suggested, therefore, that most members of each group apart from Group 4 are derived from ancestors that occupied different isolated refugia in a previous period of range fragmentation of the species, possibly at the time of the Last Glacial Maximum. Possible locations for these refugia are suggested. A comparison of mitotype diversity between northern and southern subgroups within the north‐central group of populations (Group 2) showed much greater uniformity in the northern part of the range both within and between populations. This could indicate a northward migration of the species from a southern refugium in this region during the postglacial period, although alternative explanations cannot be ruled out. Two chlorotypes were distributed across the geographical range of the species, resulting in lower levels of among‐population chlorotype variation. The geographical pattern of variation for all five chlorotypes provided some indication of the species surviving past glaciations in more than one refugium, although differentiation was much less marked, presumably due to the greater dispersal of cpDNA via pollen.     

Dissection of the genetic basis of heterosis in an elite maize hybrid by QTL mapping in an immortalized F2 population

The genetic basis of heterosis for grain yield and its components was investigated at the single- and two-locus levels using molecular markers with an immortalized F₂ (IF₂) population, which was developed by pair crosses among recombinant inbred lines (RILs) derived from the elite maize hybrid Yuyu22. Mid-parent heterosis of each cross in the IF₂ population was used to map heterotic quantitative trait loci. A total of 13 heterotic loci (HL) were detected. These included three HL for grain yield, seven for ear length, one for ear row number and two for 100-kernel weight. A total of 143 digenic interactions contributing to mid-parent heterosis were detected at the two-locus level involving all three types of interactions (additive × additive = AA, additive × dominance = AD or DA, dominance × dominance = DD). There were 25 digenic interactions for grain yield, 36 for ear length, 31 for ear row number and 51 for 100-kernel weight. Altogether, dominance effects of HL at the single-locus level as well as AA interactions played an important role in the genetic basis of heterosis for grain yield and its components in Yuyu22.     

genodive version 3.0: Easy‐to‐use software for the analysis of genetic data of diploids and polyploids

genodive version 3.0 is a user‐friendly program for the analysis of population genetic data. This version presents a major update from the previous version and now offers a wide spectrum of different types of analyses. genodive has an intuitive graphical user interface that allows direct manipulation of the data through transformation, imputation of missing data, and exclusion and inclusion of individuals, population and/or loci. Furthermore, genodive seamlessly supports 15 different file formats for importing or exporting data from or to other programs. One major feature of genodive is that it supports both diploid and polyploid data, up to octaploidy (2n = 8x) for some analyses, but up to hexadecaploidy (2n = 16x) for other analyses. The different types of analyses offered by genodive include multiple statistics for estimating population differentiation (φ_ST, F_ST, F?_ST, G_ST, G?_ST, G??_ST, D_est, R_ST, ρ), analysis of molecular variance‐based K‐means clustering, Hardy–Weinberg equilibrium, hybrid index, population assignment, clone assignment, Mantel test, Spatial Autocorrelation, 23 ways of calculating genetic distances, and both principal components and principal coordinates analyses. A unique feature of genodive is that it can also open data sets with nongenetic variables, for example environmental data or geographical coordinates that can be included in the analysis. In addition, genodive makes it possible to run several external programs (lfmm , structure , instruct and vegan ) directly from its own user interface, avoiding the need for data reformatting and use of the command line. genodive is available for computers running Mac OS X 10.7 or higher and can be downloaded freely from: http://www.patrickmeirmans.com/software .     

Isolation and characterization of microsatellite loci in the Cape fynbos heath Erica coccinea (Ericaceae)

Enriched genomic libraries were used to isolate and characterize dinucleotide microsatellite loci in Erica coccinea, a South African Cape fynbos heath species with distinct resprouter and seeder populations. Microsatellites were required to investigate the effect of the contrasting demographic pattern driven by these two post-fire responses in the population genetic structure of seeder and resprouter forms within this species. Eight microsatellite loci were characterised and amplified a total of 106 alleles in 2 samples each of 30 individuals from 1 resprouter and 1 seeder population. Mean allele numbers were 7.88 and 11.0 for the resprouter and seeder population, respectively. Both populations showed similar average observed and expected heterozygosity levels, H _O(resprouter) = 0.683, H _O(seeder) = 0.696; H _E(resprouter) = 0.726, H _E(seeder) = 0.756, and average positive inbreeding coefficients F _IS(resprouter) = 0.058, F _IS(seeder) = 0.080. This set of microsatellite loci will be used to conduct a population genetic survey of seeder and resprouter populations throughout the range of the species. Cross-species transferability was also assayed in four other South African and four European species of the genus Erica, supporting their potential use for population genetic analyses.     

Applicability of AFLP Fingerprinting Markers to Molecular Discrimination of the Three Main Fungal Species Associated with Grapevine Decline in Spain

Diplodia seriata, Phaeomoniella chlamydospora and Phaeoacremonium aleophilum are the three main species associated with grapevine decline in Spain. AFLP markers were developed to discriminate Spanish populations of these species. The markers were used to genotype isolates of D. seriata, P. chlamydospora and P. aleophilum. AFLP markers were valuable in performing population genetic studies as genetic variability (K_x) ranged from 0.07 in the P. chlamydospora population to 0.28 in the D. seriata population. Species‐specific markers obtained using only two AFLP combinations clearly discriminate D. seriata, P. chlamydospora and P. aleophilum and are a useful tool in simultaneous identification tests.     

Effects of elevated CO2 on life‐history traits of three successive generations of Frankliniella occidentalis and F. intonsa on kidney bean,Phaseolus vulgaris

The objective of this study was to determine how elevated CO₂ impacts on life‐history traits and life table parameters in three successive generations of invasive species Frankliniella occidentalis (Pergande) (Thysanoptera: Thripidae) and its related native species, Frankliniella intonsa (Pergande), fed with kidney bean leaves grown in ambient CO₂. The oviposition period, sex ratio, net reproductive rate (R₀), intrinsic rate of increase (r_m), and finite rate of increase (λ) of F. occidentalis increased in elevated CO₂, and larval duration, survival rate, mean generation time (T), and population doubling time (DT) decreased. For F. intonsa, larval duration, survival rate, oviposition period, longevity of female adults, R₀, r_m, and λ decreased in elevated CO₂, whereas sex ratio, T, and DT increased. These results indicated that the effects of elevated CO₂ would be beneficial to F. occidentalis, whereas it would be detrimental to F. intonsa. However, the effects of elevated CO₂ on F. occidentalis and F. intonsa differed over generations. In elevated CO₂, larval duration, survival rate, oviposition period, sex ratio, r_m, and λ of F. occidentalis increased linearly through successive generations, whereas T and DT decreased linearly, which suggested that the effects of elevated CO₂ on F. occidentalis would be slowly accentuated over time. For F. intonsa, larval duration, survival rate, oviposition period, r_m, and λ decreased linearly over generations, whereas sex ratio, T, and DT increased linearly. This indicated that the effects of elevated CO₂ on F. intonsa would slowly accentuate over time. We conclude that F. occidentalis would be more adapted to elevated CO₂ than F. intonsa.     

Effect of pH on survival,reproduction, egg viability and growth rate of five closely related rotifer species

In the laboratory we examined the effect of pH (5–10 with one interval) on survival, reproduction, egg viability and growth rate (intrinsic growth rate—r _m and population growth rate—r) of five Brachionus rotifer species (B. calyciflorus, B. quadridentatus, B. urceolaris, B. patulus and B. angularis). The pH was shown to exert a major influence on egg viability and growth rate (r _m and r) for each species. The age-specific survivorship curves within a species were not significantly different at pH 6–10. The optimal pH for each species is near-neutral pH (pH 6–8), and the fecundity decreased as the pH deviated from these values. For each Brachionus species, there was no significant difference between age-specific fecundity curves at pH 7 and 8. At acid pH (pH 5 or 6) higher egg mortality was observed for each species. The r _m and population r of the five Brachionus species incubated at different pHs were significantly influenced by pH. The pH supporting the highest r _m or r was obtained at pH 6–8, but varied due to species. In this study B. urceolaris and B. patulus could tolerate a broad range of pH, while the populations of B. calyciflorus, B. quadridentatus and B. angulari declined at acid conditions.     

Population genetic structure of island and mainland populations of the quokka, Setonix brachyurus (Macropodidae): a comparison of AFLP and microsatellite markers

Translocation and reintroduction are important tools for the conservation or recovery of species threatened with extinction in the wild. However, an understanding of the potential genetic consequences of mixing populations requires an understanding of the genetic variation within, and similarities among, donor and recipient populations. Genetic diversity was measured using two independent marker systems (microsatellites and AFLPs) for one island and four small remnant mainland populations of Setonix brachyurus, a threatened medium sized macropod restricted to fragmented habitat remnants and two off-shore islands in southwest Australia. Microsatellite diversity in the island population (R _s = 3.2, H _e = 71%) was similar to, or greater than, all mainland populations (R _s = 2.1–3.9, H _e = 34-71%). In contrast, AFLP diversity was significantly lower in the island population (PPL = 20.5; H _j = 0.118) compared to all mainland populations (mean PPL = 79.5–89.7; mean H _j = 0.23–0.29). Microsatellites differentiated all (mainland and island) populations from each other. However, AFLP only differentiated the island population from the mainland populations—all mainland populations were not significantly differentiated from each other for this marker. Given a known time since isolation of the island population from the mainland (6,000 years ago), and an overall more conservative rate of evolution of AFLP markers, our results are consistent with mainland populations fragmenting thousands of years ago (but <6,000 years), probably as a consequence of reduced rainfall and the constriction of the preferred mesic habitat of quokkas. Our results also support a recent history of severe population bottlenecks in mainland populations, and a long history of bottlenecks of the island population, but reflect a recent explosion in numbers since European occupation of the island. Our results indicate that translocation of island populations to supplement mainland populations would introduce genetically markedly differentiated, and possibly maladapted, individuals.     

Stability of the analytical solution of Penna model of biological aging

There are some analytical solutions of the Penna model of biological aging; here, we discuss the approach by Coe et al. (Phys. Rev. Lett. 89, 288103, 2002), based on the concept of self-consistent solution of a master equation representing the Penna model. The equation describes transition of the population distribution at time t to next time step (t + 1). For the steady state, the population n(a, l, t) at age a and for given genome length l becomes time-independent. In this paper we discuss the stability of the analytical solution at various ranges of the model parameters—the birth rate b or mutation rate m. The map for the transition from n(a, l, t) to the next time step population distribution n(a + 1, l, t + 1) is constructed. Then the fix point (the steady state solution) brings recovery of Coe et al. results. From the analysis of the stability matrix, the Lyapunov coefficients, indicative of the stability of the solutions, are extracted. The results lead to phase diagram of the stable solutions in the space of model parameters (b, m, h), where h is the hunt rate. With increasing birth rate b, we observe critical b ₀ below which population is extinct, followed by non-zero stable single solution. Further increase in b leads to typical series of bifurcations with the cycle doubling until the chaos is reached at some b _c. Limiting cases such as those leading to the logistic model are also discussed.     

Adaptability to elevated temperature and nitrogen addition is greater in a high-elevation population than in a low-elevation population of Hippophae rhamnoides

Hippophae rhamnoides L. is a broadleaf deciduous woody shrub occurring in southwest China, where it has been widely used in ecological restoration. In this study we investigated growth and physiological responses of 2-year-old healthy seedlings to elevated temperature, nitrogen (N) addition and their combination in two contrasting populations from high and low elevations. In closed-top chamber experiments, two populations were subjected to two temperature conditions (ambient temperature and temperature elevated by 2.2 ± 0.2°C) and two N levels (0 and 25 g N m⁻² a⁻¹). Compared with the control, increases in total leaf area (TLA), total chlorophyll content (TC), light-saturated photosynthetic rate (P _max), guaiacol peroxidase activity (POD), catalase activity (CAT) and carbon isotope composition (δ¹³C) were greater in the high-elevation population than in the low-elevation population under elevated temperature. On the other hand, decreases in root and shoot biomass ratio (RS), TC, P _max, light saturation point (L _SP), light compensation point (L _CP), superoxide dismutase (SOD), POD, CAT and δ¹³C were lower in the high-elevation population than in the low-elevation population under N addition. Moreover, the combination of elevated temperature and N addition decreased RS, P _max, apparent quantum efficiency (Φ), SOD, POD and δ¹³C significantly more in the low-elevation population than in the high-elevation population. These results demonstrated that there are different adaptive strategies among H. rhamnoides populations, the high-elevation population exhibiting higher adaptability to elevated temperature and N addition than the low-elevation population.     

Modeling the Dynamics of Viral Evolution Considering Competition Within Individual Hosts and at Population Level: The Effects of Treatment

We consider two viral strains competing against each other within individual hosts (at cellular level) and at population level (for infecting hosts) by studying two cases. In the first case, the strains do not mutate into each other. In this case, we found that each individual in the population can be infected by only one strain and that co-existence in the population is possible only when the strain that has the greater basic intracellular reproduction number, R _0c , has the smaller population number R _0p . Treatment against the one strain shifts the population equilibrium toward the other strain in a complicated way (see Appendix B). In the second case, we assume that the strain that has the greater intracellular number R _0c can mutate into the other strain. In this case, individual hosts can be simultaneously infected by both strains (co-existence within the host). Treatment shifts the prevalence of the two strains within the hosts, depending on the mortality induced by the treatment, which is, in turn, dependent upon the doses given to each individual. The relative proportions of the strains at the population level, under treatment, depend both on the relative proportions within the hosts (which is determined by the dosage of treatment) and on the number of individuals treated per unit time, that is, the rate of treatment. Implications for cases of real diseases are briefly discussed.     

Tandem repeat markers for population genetic studies of the protected California tiger salamander (<Emphasis Type="Italic">Ambystoma californiense</Emphasis>)

Habitat loss is the single greatest threat to persistence of the critically threatened California tiger salamander (Ambystoma californiense). To aid management plans that designate critical habitat for this species, I developed and characterized 21 tetranucleotide microsatellite markers using two native populations in Santa Barbara and Alameda Counties. Allelic variation and average heterozygosities were lower in the endangered Santa Barbara population (allele range 1–4, mean 2.4; H _O = 0.308 H _E = 0.288) compared with the threatened Alameda population (allele range 2–10, mean 6.7; H _O = 0.712, H _E = 0.722). In-depth population studies using these markers will provide vital information for plans to assign critical habitat that optimize gene flow among breeding populations, as well as for identifying non-native hybrid genotypes that threaten native A. californiense stocks. Beyond the conservation goals for A. californiense, the close phylogenetic relationships within the tiger salamander complex also suggest a broad utility for population studies using these markers.     

Development and characterisation of microsatellite loci in Farfugium japonicum (Asteraceae)

Eight microsatellite loci for the perennial herb Farfugium japonicum, including the rheophytic variety luchuense endemic to riparian areas of the Ryukyu Islands, Japan, were isolated and characterised. The number of alleles ranged from 5 to 14. The expected (H _E) and observed (H _O) heterozygosities were 0.344–0.885 and 0.121–0.754, respectively, from 69 individuals in one population. Six loci exhibited significantly fewer heterozygotes than expected under Hardy–Weinberg equilibrium (P < 0.05). The primers amplifying microsatellite sequences in F. japonicum may provide a population genetics tool useful in the establishment of a conservation strategy.     

Microsatellite loci for the critically endangered Santa Cruz long-toed salamander (Ambystoma macrodactylum croceum) and other Ambystoma taxa

The imperiled Santa Cruz long-toed salamander (Ambystoma macrodactylum croceum) is a geographically isolated lineage of the long-toed salamander species complex. Recovery of this taxon requires understanding its population genetic structure, which is necessary to identify critical habitat. To assess the population structure of A. m. croceum, I developed primers for 17 microsatellite loci and evaluated their performance. Polymorphism ranged from 3 to 10 alleles per locus, with an average of 6.3 per locus. Mean observed and expected heterozygosity were similar (H _O = 0.678; H _E = 0.708), with one locus demonstrating heterozygote deficiency. The performance of these and previously identified microsatellite loci were also characterized for two other long-toed salamander lineages, as well as two other species of Ambystoma. The utility of these markers has similar and overlapping applications for conservation and population genetics studies of other Ambystoma taxa, and will provide a basis for revising the recovery plan for A. m. croceum, particularly for designating corridors among interacting breeding populations and upland areas, and identifying areas where gene flow should be restored.     

Polymorphic microsatellite loci in the Coffee Berry Borer,Hypothenemus hampei (Coleoptera,Scolytidae)

Microsatellite loci were isolated from the coffee berry borer, Hypothenemus hampei, which is regarded as the major pest in coffee cultures. Seven polymorphic loci were obtained from an enriched genomic library. A low to moderate genetic diversity was observed per locus, with an observed number of alleles ranging from two to five in the 39 unrelated individuals sampled from Ethiopia. A clear deficit of heterozygotes within the population (mean heterozygosities, H_O = 0.10/H_E = 0.50) and an extreme inbreeding (F_IS, 0.70–1.00) were demonstrated. Cross‐species amplifications showed that some of the markers could be useful in two closely related Hypothenemus species.     

Naturalization of plant populations: the role of cultivation and population size and density

Field experimentation is required to assess the effects of environmental stochasticity on small immigrant plant populations—a widely understood but largely unexplored aspect of predicting any species’ likelihood of naturalization and potential invasion. Cultivation can mitigate this stochasticity, although the outcome for a population under cultivation nevertheless varies enormously from extinction to persistence. Using factorial experiments, we investigated the effects of population size, density, and cultivation (irrigation) on the fate of founder populations for four alien species with different life history characteristics (Echinochloa frumentacea, Fagopyrum esculentum, Helianthus annuus, and Trifolium incarnatum) in eastern Washington, USA. The fate of founder populations was highly variable within and among the 3 years of experimentation and illustrates the often precarious environment encountered by plant immigrants. Larger founder populations produced more seeds (P < 0.001); the role of founder population size, however, differed among years. Irrigation resulted in higher percent survival (P < 0.001) and correspondingly larger net reproductive rate (R ₀; P < 0.001). But the minimum level of irrigation for establishment, R ₀ > 1, differed among years and species. Sowing density did not affect the likelihood of establishment for any species. Our results underscore the importance of environmental stochasticity in determining the fate of founder populations and the potential of cultivation and large population size in countering the long odds against naturalization. Any implementation of often proposed post-immigration field trials to assess the risk of an alien species becoming naturalized, a requisite step toward invasion, will need to assess different sizes of founder populations and the extent and character of cultivation (intentional or unintentional) that the immigrants might receive.     

Mitochondrial DNA diversity of the vulnerable Chinese Egret (Egretta eulophotes) from China

The genetic diversity and population structure of the vulnerable Chinese Egret (Egretta eulophotes) were surveyed in the present study from three archipelagoes that cover the most southerly to the very northerly parts of the Chinese distribution range of this species, using a 433-bp fragment of the mitochondrial control region (CR). Among 90 individual samples, 31 different haplotypes were defined by 30 polymorphic sites. Overall haplotype diversity, nucleotide diversity and mean sequence divergence (p-distance) of this egret were 0.920, 0.0088 and 1.11%, respectively. NJ tree and parsimony network for the CR haplotypes of the Chinese Egret showed little genetic structure, and analysis of molecular variance indicated low but significant genetic differentiation (haplotype-based Φ_ST = 0.03267, P < 0.05 and distance-based Φ_ST = 0.04194, P < 0.05) among populations. The significant Fu’s F _S tests (Fu’s F _S  = −16.946, P < 0.01) and mismatch distribution analysis (τ = 4.463, SSD = 0.0081, P = 0.12) suggested that the low genetic differentiation and little geographical structure of the genetic differentiation might be explained by the population expansion. The Mantel test (haplotype-based F _ST, r = 0.639, P = 0.34 and distance-based F _ST, r = 0.947, P = 0.15) suggest that the significant genetic differentiation among populations was likely due to isolation by distance.