Sex‐chromosome differentiation parallels postglacial range expansion in European tree frogs (Hyla arborea)

Occasional XY recombination is a proposed explanation for the sex‐chromosome homomorphy in European tree frogs. Numerous laboratory crosses, however, failed to detect any event of male recombination, and a detailed survey of NW‐European Hyla arborea populations identified male‐specific alleles at sex‐linked loci, pointing to the absence of XY recombination in their recent history. Here, we address this paradox in a phylogeographic framework by genotyping sex‐linked microsatellite markers in populations and sibships from the entire species range. Contrasting with postglacial populations of NW Europe, which display complete absence of XY recombination and strong sex‐chromosome differentiation, refugial populations of the southern Balkans and Adriatic coast show limited XY recombination and large overlaps in allele frequencies. Geographically and historically intermediate populations of the Pannonian Basin show intermediate patterns of XY differentiation. Even in populations where X and Y occasionally recombine, the genetic diversity of Y haplotypes is reduced below the levels expected from the fourfold drop in copy numbers. This study is the first in which X and Y haplotypes could be phased over the distribution range in a species with homomorphic sex chromosomes; it shows that XY‐recombination patterns may differ strikingly between conspecific populations, and that recombination arrest may evolve rapidly (<5000 generations).     

No evidence that Y‐chromosome differentiation affects male fitness in a Swiss population of common frogs

The canonical model of sex‐chromosome evolution assigns a key role to sexually antagonistic (SA) genes on the arrest of recombination and ensuing degeneration of Y chromosomes. This assumption cannot be tested in organisms with highly differentiated sex chromosomes, such as mammals or birds, owing to the lack of polymorphism. Fixation of SA alleles, furthermore, might be the consequence rather than the cause of recombination arrest. Here we focus on a population of common frogs (Rana temporaria) where XY males with genetically differentiated Y chromosomes (nonrecombinant Y haplotypes) coexist with both XY° males with proto‐Y chromosomes (only differentiated from X chromosomes in the immediate vicinity of the candidate sex‐determining locus Dmrt1) and XX males with undifferentiated sex chromosomes (genetically identical to XX females). Our study finds no effect of sex‐chromosome differentiation on male phenotype, mating success or fathering success. Our conclusions rejoin genomic studies that found no differences in gene expression between XY, XY° and XX males. Sexual dimorphism in common frogs might result more from the differential expression of autosomal genes than from sex‐linked SA genes. Among‐male variance in sex‐chromosome differentiation seems better explained by a polymorphism in the penetrance of alleles at the sex locus, resulting in variable levels of sex reversal (and thus of X‐Y recombination in XY females), independent of sex‐linked SA genes.     

Dmrt1 polymorphism and sex‐chromosome differentiation in Rana temporaria

Sex‐determination mechanisms vary both within and among populations of common frogs, opening opportunities to investigate the molecular pathways and ultimate causes shaping their evolution. We investigated the association between sex‐chromosome differentiation (as assayed from microsatellites) and polymorphism at the candidate sex‐determining gene Dmrt1 in two Alpine populations. Both populations harboured a diversity of X‐linked and Y‐linked Dmrt1 haplotypes. Some males had fixed male‐specific alleles at all markers (“differentiated” Y chromosomes), others only at Dmrt1 (“proto‐” Y chromosomes), while still others were genetically indistinguishable from females (undifferentiated X chromosomes). Besides these XX males, we also found rare XY females. The several Dmrt1 Y haplotypes differed in the probability of association with a differentiated Y chromosome, which we interpret as a result of differences in the masculinizing effects of alleles at the sex‐determining locus. From our results, the polymorphism in sex‐chromosome differentiation and its association with Dmrt1, previously inferred from Swedish populations, are not just idiosyncratic features of peripheral populations, but also characterize highly diverged populations in the central range. This implies that an apparently unstable pattern has been maintained over long evolutionary times.     

A cryptic heterogametic transition revealed by sex-linked DNA markers in Palearctic green toads

In sharp contrast to birds and mammals, most cold‐blooded vertebrates have homomorphic (morphologically undifferentiated) sex chromosomes. This might result either from recurrent X‐Y recombination (occurring e.g. during occasional events of sex reversal) or from frequent turnovers (during which sex‐determining genes are overthrown by new autosomal mutations). Evidence for turnovers is indeed mounting in fish, but very few have so far been documented in amphibians, possibly because of practical difficulties in identifying sex chromosomes. Female heterogamety (ZW) has long been established in Bufo bufo, based on sex reversal and crossing experiments. Here, we investigate a sex‐linked marker identified from a laboratory cross between Palearctic green toads (Bufo viridis subgroup). The F₁ offspring produced by a female Bufo balearicus and a male Bufo siculus were phenotypically sexed, displaying an even sex ratio. A sex‐specific marker detected in highly reproducible AFLP genotypes was cloned. Sequencing revealed a noncoding, microsatellite‐containing fragment. Reamplification and genotyping of families of this and a reciprocal cross showed B. siculus to be male heterogametic (XY) and suggested the same system for B. balearicus. Our results thus reveal a cryptic heterogametic transition within bufonid frogs and help explain patterns of hybrid fitness within the B. viridis subgroup. Turnovers of genetic sex‐determination systems may be more frequent in amphibians than previously thought and thus contribute to the prevalence of homomorphic sex chromosomes in this group.     

Identification and characterization of nuclear microsatellite loci in the aquatic moss Platyhypnidium riparioides (Brachytheciaceae)

Eight microsatellite loci from the aquatic moss Platyhypnidium riparioides were identified using the method of microsatellite‐enriched libraries. Polymorphism was assessed in a sample of four populations of 20 individuals each from four streams of the Meuse hydrographic basin in southern Belgium. The markers amplified three to seven alleles per locus. Comparison of observed and expected heterozygosities as well as F‐statistics (F_ST = 0.62) reveals a significant genetic differentiation among populations. These markers will be useful for further investigation of population genetic structure and diversity at different nested spatial scales.     

Comparison of quantitative and molecular variation in agroforestry populations of the shea tree (Vitellaria paradoxa C.F. Gaertn) in Mali

In this study we investigated the within- and between-population genetic variation using microsatellite markers and quantitative traits of the shea tree, Vitellaria paradoxa, an important agroforestry tree species of the Sudano–Sahelian region in Africa. Eleven populations were sampled across Mali and in northern Côte d’Ivoire. Leaf size and form and growth traits were measured in a progeny test at the nursery stage. Eight microsatellites were used to assess neutral genetic variation. Low levels of heterozygosity were recorded (1.6–3.0 alleles/locus; H_E = 0.25–0.42) and the fixation index (F_IS = −0.227–0.186) was not significantly different from zero suggesting that Hardy–Weinberg equilibrium is encountered in all populations sampled. Quantitative traits exhibited a strong genetic variation between populations and between families within populations. The degree of population differentiation of the quantitative traits (Q_ST = 0.055–0.283, Q_STmean = 0.189) strongly exceeds that in eight microsatellite loci (F_ST = −0.011–0.142, F_STmean = 0.047). Global and pairwise F_ST values were very low and not significantly different from zero suggesting agroforestry practices are amplifying gene flow (Nm = 5.07). The population means for quantitative traits and the rainfall variable were not correlated, showing variation was not linked with this climatic cline. It is suggested that this marked differentiation for quantitative traits, independent of environmental clines and despite a high gene flow, is a result of local adaptation and human selection of shea trees. This process has induced high linkage disequilibrium between underlying loci of polygenic characters.     

Within‐population polymorphism of sex‐determination systems in the common frog (Rana temporaria)

In sharp contrast with birds and mammals, the sex chromosomes of ectothermic vertebrates are often undifferentiated, for reasons that remain debated. A linkage map was recently published for Rana temporaria (Linnaeus, 1758) from Fennoscandia (Eastern European lineage), with a proposed sex‐determining role for linkage group 2 (LG2). We analysed linkage patterns in lowland and highland populations from Switzerland (Western European lineage), with special focus on LG2. Sibship analyses showed large differences from the Fennoscandian map in terms of recombination rates and loci order, pointing to large‐scale inversions or translocations. All linkage groups displayed extreme heterochiasmy (total map length was 12.2 cM in males, versus 869.8 cM in females). Sex determination was polymorphic within populations: a majority of families (with equal sex ratios) showed a strong correlation between offspring phenotypic sex and LG2 paternal haplotypes, whereas other families (some of which with female‐biased sex ratios) did not show any correlation. The factors determining sex in the latter could not be identified. This coexistence of several sex‐determination systems should induce frequent recombination of X and Y haplotypes, even in the absence of male recombination. Accordingly, we found no sex differences in allelic frequencies on LG2 markers among wild‐caught male and female adults, except in one high‐altitude population, where nonrecombinant Y haplotypes suggest sex to be entirely determined by LG2. Multifactorial sex determination certainly contributes to the lack of sex‐chromosome differentiation in amphibians.     

Novel Y‐chromosome short tandem repeats in Sus scrofa and their variation in European wild boar and domestic pig populations

Y‐chromosome markers are important tools for studying male‐specific gene flow within and between populations, hybridization patterns and kinship. However, their use in non‐human mammals is often hampered by the lack of Y‐specific polymorphic markers. We identified new male‐specific short tandem repeats (STRs) in Sus scrofa using the available genome sequence. We selected four polymorphic loci (5–10 alleles per locus), falling in one duplicated and two single‐copy regions. A total of 32 haplotypes were found by screening 211 individuals from eight wild boar populations across Europe and five domestic pig populations. European wild boar were characterized by significantly higher levels of haplotype diversity compared to European domestic pigs (H_D = 0.904 ± 0.011 and H_D = 0.491 ± 0.077 respectively). Relationships among STR haplotypes were investigated by combining them with single nucleotide polymorphisms at two linked genes (AMELY and UTY) in a network analysis. A differentiation between wild and domestic populations was observed (F_ST = 0.229), with commercial breeds sharing no Y haplotype with the sampled wild boar. Similarly, a certain degree of geographic differentiation was observed across Europe, with a number of local private haplotypes and high diversity in northern populations. The described Y‐chromosome markers can be useful to track male inheritance and gene flow in wild and domestic populations, promising to provide insights into evolutionary and population genetics in Sus scrofa.     

Cryptic recombination in the ever-young sex chromosomes of Hylid frogs

Sex chromosomes are expected to evolve suppressed recombination, which leads to degeneration of the Y and heteromorphism between the X and Y. Some sex chromosomes remain homomorphic, however, and the factors that prevent degeneration of the Y in these cases are not well understood. The homomorphic sex chromosomes of the European tree frogs (Hyla spp.) present an interesting paradox. Recombination in males has never been observed in crossing experiments, but molecular data are suggestive of occasional recombination between the X and Y. The hypothesis that these sex chromosomes recombine has not been tested statistically, however, nor has the X‐Y recombination rate been estimated. Here, we use approximate Bayesian computation coupled with coalescent simulations of sex chromosomes to quantify X‐Y recombination rate from existent data. We find that microsatellite data from H. arborea, H. intermedia and H. molleri support a recombination rate between X and Y that is significantly different from zero. We estimate that rate to be approximately 10⁵ times smaller than that between X chromosomes. Our findings support the notion that very low recombination rate may be sufficient to maintain homomorphism in sex chromosomes.     

A genetic linkage map of red drum,Sciaenops ocellatus

Second‐generation, sex‐specific genetic linkage maps were generated for the economically important estuarine‐dependent marine fish Sciaenops ocellatus (red drum). The maps were based on F₁ progeny from each of two single‐pair mating families. A total of 237 nuclear‐encoded microsatellite markers were mapped to 25 linkage groups. The female map contained 226 markers, with a total length of 1270.9 centiMorgans (cM) and an average inter‐marker interval of 6.53 cM; the male map contained 201 markers, with a total length of 1122.9 cM and an average inter‐marker interval of 6.03 cM. The overall recombination rate was approximately equal in the two sexes (♀:♂ = 1.03:1). Recombination rates in a number of linkage intervals, however, differed significantly between the same sex in both families and between sexes within families. The former occurred in 2.4% of mapped intervals, while the latter occurred in 51.2% of mapped intervals. Sex‐specific recombination rates varied within chromosomes, with regions of both female‐biased and male‐biased recombination. Original clones from which the microsatellite markers were generated were compared with genome sequence data for the spotted green puffer, Tetraodon nigroviridis; a total of 43 matches were located in 17 of 21 chromosomes of T. nigroviridis, while seven matches were in unknown portions of the T. nigroviridis genome. The map for red drum provides a new, useful tool for aquaculture, population genetics, and comparative genomics of this economically important marine species.     

Sequence Polymorphism and Geographical Variation at a Positively Selected MHC-DRB Gene in the Finless Porpoise (Neophocaena phocaenoides): Implication for Recent Differentiation of the Yangtze Finless Porpoise?

Sequence polymorphism at the MHC class II DRB locus was investigated in three finless porpoise (Neophocaena phocaenoides) populations in Chinese waters. Intragenic recombination and strong positive selection were the main forces in generating sequence diversity in the DRB gene. MHC sequence diversity changed significantly along the study period. Significant decrease in heterozygosity and lost alleles have been detected in the Yangtze River population and South China Sea population since 1990. Furthermore, there is a trend of increasing population differentiation over time. Especially, the genetic differentiation between the Yangtze River population and the Yellow Sea population was very low prior to 1990 (F _ST = 0.036, P = 0.009), but became very significant after 1990 (F _ST = 0.134, P < 0.001), suggesting a recent augmentation of genetic differentiation between both populations probably in a relatively short-term period. Porpoises from the Yangtze River displayed divergent frequencies of shared and private alleles from those displayed by two marine populations, which suggest that the former riverine population has been under a different selection regime (characteristic of a fresh water environment) than that of its marine counterparts.     

The Pillars of Hercules as a bathymetric barrier to gene flow promoting isolation in a global deep‐sea shark (Centroscymnus coelolepis)

Knowledge of the mechanisms limiting connectivity and gene flow in deep‐sea ecosystems is scarce, especially for deep‐sea sharks. The Portuguese dogfish (Centroscymnus coelolepis) is a globally distributed and near threatened deep‐sea shark. C. coelolepis population structure was studied using 11 nuclear microsatellite markers and a 497‐bp fragment from the mtDNA control region. High levels of genetic homogeneity across the Atlantic (Φ_ST = ?0.0091, F_ST = 0.0024, P > 0.05) were found suggesting one large population unit at this basin. The low levels of genetic divergence between Atlantic and Australia (Φ_ST = 0.0744, P < 0.01; F_ST = 0.0015, P > 0.05) further suggested that this species may be able to maintain some degree of genetic connectivity even across ocean basins. In contrast, sharks from the Mediterranean Sea exhibited marked genetic differentiation from all other localities studied (Φ_ST = 0.3808, F_ST = 0.1149, P < 0.001). This finding suggests that the shallow depth of the Strait of Gibraltar acts as a barrier to dispersal and that isolation and genetic drift may have had an important role shaping the Mediterranean shark population over time. Analyses of life history traits allowed the direct comparison among regions providing a complete characterization of this shark's populations. Sharks from the Mediterranean had markedly smaller adult body size and size at maturity compared to Atlantic and Pacific individuals. Together, these results suggest the existence of an isolated and unique population of C. coelolepis inhabiting the Mediterranean that most likely became separated from the Atlantic in the late Pleistocene.     

Genetically distinct populations of northern shrimp,Pandalus borealis,in the North Atlantic: adaptation to different temperatures as an isolation factor

The large‐scale population genetic structure of northern shrimp, Pandalus borealis, was investigated over the species’ range in the North Atlantic, identifying multiple genetically distinct groups. Genetic divergence among sample localities varied among 10 microsatellite loci (range: F_ST = ?0.0002 to 0.0475) with a highly significant average (F_ST = 0.0149; P < 0.0001). In contrast, little or no genetic differences were observed among temporal replicates from the same localities (F_ST = 0.0004; P = 0.33). Spatial genetic patterns were compared to geographic distances, patterns of larval drift obtained through oceanographic modelling, and temperature differences, within a multiple linear regression framework. The best‐fit model included all three factors and explained approximately 29% of all spatial genetic divergence. However, geographic distance and larval drift alone had only minor effects (2.5–4.7%) on large‐scale genetic differentiation patterns, whereas bottom temperature differences explained most (26%). Larval drift was found to promote genetic homogeneity in parts of the study area with strong currents, but appeared ineffective across large temperature gradients. These findings highlight the breakdown of gene flow in a species with a long pelagic larval phase (up to 3 months) and indicate a role for local adaptation to temperature conditions in promoting evolutionary diversification and speciation in the marine environment.     

Phylogeography,more than elevation,accounts for sex chromosome differentiation in Swiss populations of the common frog (Rana temporaria)*

Sex chromosomes in vertebrates range from highly heteromorphic (as in most birds and mammals) to strictly homomorphic (as in many fishes, amphibians, and nonavian reptiles). Reasons for these contrasted evolutionary trajectories remain unclear, but species such as common frogs with polymorphism in the extent of sex chromosome differentiation may potentially deliver important clues. By investigating 92 common frog populations from a wide range of elevations throughout Switzerland, we show that sex chromosome differentiation strongly correlates with alleles at the candidate sex-determining gene Dmrt1. Y-specific Dmrt1 haplotypes cluster into two main haplogroups, Y_A and Y_B, with a phylogeographic signal that parallels mtDNA haplotypes: Y_A populations, with mostly well-differentiated sex chromosomes, occur primarily south of the main alpine ridge that bisects Switzerland, whereas Y_B populations, with mostly undifferentiated (proto-)sex chromosomes, occur north of this ridge. Elevation has only a marginal effect, opposing previous suggestions of a major role for climate on sex chromosome differentiation. The Y-haplotype effect might result from differences in the penetrance of alleles at the sex-determining locus (such that sex reversal and ensuing X-Y recombination are more frequent in Y_B populations), and/or fixation of an inversion on Y_A (as supported by the empirical observation that Y_A haplotypes might not recombine in XY_A females).     

Development of polymorphic microsatellite loci in the hermaphroditic freshwater snails Drepanotrema surinamense and Drepanotrema depressissimum

We characterized new variable microsatellites in two congeneric species of hermaphroditic freshwater snails (Drepanotrema depressissimum and D. surinamense), as well as conditions for multiplexing and simultaneously genotyping sets of loci. D. depressissimum had high mean gene diversity (> 0.8) and large number of alleles (= 10.9) per population. Most loci and populations were at Hardy–Weinberg equilibrium. The F_ST estimates were low among lesser Antilles populations and larger with a Venezuelan population. Far less diversity was found in D. surinamense with mean number of alleles and gene diversity per population of 2.8 and 0.34, respectively. Very few heterozygous individuals were observed. The most likely explanation is a high selfing rate (> 0.825) in this species. Unsurprisingly, the F_ST estimates among populations were much higher than in D. depressissimum. Cross‐species amplification in three congeneric species was on the whole unsuccessful.     

Mitochondrial DNA phylogeography of Spodoptera exigua across a broad geographic area in China

The beet armyworm, Spodoptera exigua (Lepidoptera: Noctuidae), is an important agriculture pest in China that causes serious economic losses in some of the main crop‐producing areas. To monitor and manage this pest effectively, it is necessary to investigate its phylogeographic patterns in China. In this study, we used a partial sequence of mitochondrial DNA Cytb gene consisting of 724 bps to investigate the genetic diversity of the beet armyworm. A total of 765 individuals from 47 populations across the main distribution range of the species were collected, and 112 haplotypes were identified. Moderate‐to‐high levels of genetic diversity (Hd = 0.672 ± 0.017, Pi = 0.00268 ± 0.00021) for the total populations were obtained. Phylogenetic and median‐joining network analyses indicated there was no distinct geographic distribution pattern among haplotypes. Overall, the study also revealed significant differentiation among some populations (P < 0.05). The F_ST values of Shenyang population (SY2012–SY 2014), as well as Baoding (BD), Taian (TA), Lucheng (LC), Zhengzhou (ZZ) and Wuhan (WH), were significantly different from those of the populations in most other locations. Hierarchical AMOVA showed there was no significantly genetic structure between populations located in seven geographic regions and four main bioclimatic zones. Finally, unimodal mismatch distribution combined with negative Tajima's D (D = ?2.696, P < 0.001) and Fu's F_S (F_S = ?207.228, P > 0.05) indicated recent population expansion of S. exigua at large spatial scales in China.     

Size differentiation in Finnish house sparrows follows Bergmann's rule with evidence of local adaptation

Bergmann's rule predicts that individuals are larger in more poleward populations and that this size gradient has an adaptive basis. Hence, phenotypic divergence in size traits between populations (P_ST) is expected to exceed the level of divergence by drift alone (F_ST). We measured 16 skeletal traits, body mass and wing length in 409 male and 296 female house sparrows Passer domesticus sampled in 12 populations throughout Finland, where the species has its northernmost European distributional margin. Morphometric differentiation across populations (P_ST) was compared with differentiation in 13 microsatellites (F_ST). We find that twelve traits phenotypically diverged more than F_ST in both sexes, and an additional two traits diverged in males. The phenotypic divergence exceeded F_ST in several traits to such a degree that findings were robust also to strong between‐population environmental effects. Divergence was particularly strong in dimensions of the bill, making it a strong candidate for the study of adaptive molecular genetic divergence. Divergent traits increased in size in more northern populations. We conclude that house sparrows show evidence of an adaptive latitudinal size gradient consistent with Bergmann's rule on the modest spatial scale of ca. 600 km.     

Twelve polymorphic expressed sequence tags‐derived markers for Plasmopara halstedii,the causal agent of sunflower downy mildew

Twelve expressed sequence tags‐derived markers were isolated from Plasmopara halstedii (Oomycetes), the causal agent of sunflower downy mildew. A total of 25 single nucleotide polymorphisms and five indels were detected by single‐strand conformation polymorphism analysis and developed for high‐throughput genotyping of 32 isolates. There was a high level of genetic diversity (H_E = 0.484). Observed heterozygosity ranged from 0 to 0.143 indicating that P. halstedii is probably a selfing species. These markers were also useful in detecting significant genetic variations among French populations (F_ST = 0.193) and between French and Russian populations (F_ST = 0.23). Cross‐amplification tests on three closely related species indicated that no loci amplified in other Oomycete species.     

Generalization of the QST framework in hierarchically structured populations: Impacts of inbreeding and dominance

Q_ST is a differentiation parameter based on the decomposition of the genetic variance of a trait. In the case of additive inheritance and absence of selection, it is analogous to the genic differentiation measured on individual loci, F_ST. Thus, Q_ST?F_ST comparison is used to infer selection: selective divergence when Q_ST > F_ST, or convergence when Q_ST < F_ST. The definition of Q‐statistics was extended to two‐level hierarchical population structures with Hardy–Weinberg equilibrium. Here, we generalize the Q‐statistics framework to any hierarchical population structure. First, we developed the analytical definition of hierarchical Q‐statistics for populations not at Hardy–Weinberg equilibrium. We show that the Q‐statistics values obtained with the Hardy–Weinberg definition are lower than their corresponding F‐statistics when F_IS > 0 (higher when F_IS < 0). Then, we used an island model simulation approach to investigate the impact of inbreeding and dominance on the Q_ST?F_ST framework in a hierarchical population structure. We show that, while differentiation at the lower hierarchical level (Q_SR) is a monotonic function of migration, differentiation at the upper level (Q_RT) is not. In the case of additive inheritance, we show that inbreeding inflates the variance of Q_RT, which can increase the frequency of Q_RT > F_RT cases. We also show that dominance drastically reduces Q‐statistics below F‐statistics for any level of the hierarchy. Therefore, high values of Q‐statistics are good indicators of selection, but low values are not in the case of dominance.     

Relative performance of Bayesian clustering software for inferring population substructure and individual assignment at low levels of population differentiation

Traditional methods for characterizing genetic differentiation among populations rely on a priori grouping of individuals. Bayesian clustering methods avoid this limitation by using linkage and Hardy–Weinberg disequilibrium to decompose a sample of individuals into genetically distinct groups. There are several software programs available for Bayesian clustering analyses, all of which describe a decrease in the ability to detect distinct clusters as levels of genetic differentiation among populations decrease. However, no study has yet compared the performance of such methods at low levels of population differentiation, which may be common in species where populations have experienced recent separation or high levels of gene flow. We used simulated data to evaluate the performance of three Bayesian clustering software programs, PARTITION, STRUCTURE, and BAPS, at levels of population differentiation below F _ST=0.1. PARTITION was unable to correctly identify the number of subpopulations until levels of F _ST reached around 0.09. Both STRUCTURE and BAPS performed very well at low levels of population differentiation, and were able to correctly identify the number of subpopulations at F _ST around 0.03. The average proportion of an individual’s genome assigned to its true population of origin increased with increasing F _ST for both programs, reaching over 92% at an F _ST of 0.05. The average number of misassignments (assignments to the incorrect subpopulation) continued to decrease as F _ST increased, and when F _ST was 0.05, fewer than 3% of individuals were misassigned using either program. Both STRUCTURE and BAPS worked extremely well for inferring the number of clusters when clusters were not well-differentiated (F _ST=0.02–0.03), but our results suggest that F _ST must be at least 0.05 to reach an assignment accuracy of greater than 97%.