Genetic variability in the European bison (Bison bonasus) population from Białowieża forest over 50 years

The aim of this study was to assess potential post-bottleneck temporal genetic differentiation following the reintroduction of the species into the Białowieża Forest. Variability of 12 polymorphic microsatellite markers was analysed for 178 individuals born between 1955 and 2005, divided by birth year into five temporal groups. Low overall allelic richness (AR) per locus (AR = 2.0) and a low overall expected heterozygosity (H_E = 0.28) were observed. The overall F _IS was not significantly different from zero. The mean F _IS values were, however, significantly different from zero for individuals born between 1955 and 1965 ( F _IS = 0.19). A Bayesian computation was used to estimate effective population size (N_e) for each temporal group. We observed relatively small values of N_e ranging from 7.9 to 28.4. The low N_e values confirm that, despite a rapid growth of the bison population following the founder event, N_e increased only slowly.  © 2009 The Linnean Society of London, Biological Journal of the Linnean Society , 2009, 97 , 801–809.     

Structure of allozyme variation in Nordic Silene nutans (Caryophyllaceae): population size, geographical position and immigration history

We investigated allozyme variation in 34 populations of the perennial herb Silene nutans from Sweden and northern Finland, areas that were ice-covered during the last (Weichselian) glaciation. The present geographical structure of genetic variation in S. nutans in Sweden and northern Finland appears to have been mainly shaped by ancient historical processes. Patterns of variation in allele frequencies suggest two major postglacial immigration routes into Sweden, with populations entering the area from both the south and the east and forming a contact zone with admixed populations in central Sweden. While estimates of within-population genetic diversity and allelic richness are significantly correlated with present population size and geographical position (latitude), population size is not correlated with latitude. Low genetic diversity in the northern populations is more likely to have resulted from ancient stochastic events during the process of immigration than from recent population fragmentation. F _IS values are high and increase with latitude. Evidence of recent bottlenecks was detected in several southern Swedish populations: these can be interpreted in terms of population fragmentation as a result of anthropogenic disturbance. Soil pH is uncorrelated with population size and position.  © 2004 The Linnean Society of London, Biological Journal of the Linnean Society, 2004, 81 , 357–371.     

Stability of genetic structure and effective population size inferred from temporal changes of microsatellite DNA polymorphisms in the land snail Helix aspersa (Gastropoda: Helicidae)

Temporal evolution of genetic variability may have far-reaching consequences for a diverse array of evolutionary processes. Within the polders of the Bay of Mont-Saint-Michel (France), populations of the land snail Helix aspersa are characterized by a metapopulation structure with occasional extinction processes resulting from farming practices. A temporal survey of genetic structure in H . aspersa was carried out using variability at four microsatellite loci, in ten populations sampled two years apart. Levels of within-population genetic variation, as measured by allelic richness, H _e or F _is , did not change over time and similar levels of population differentiation were demonstrated for both sampling years. The extent of genetic differentiation between temporal samples of the same population established (i) a stable structure for six populations, and (ii) substantial genetic changes for four populations. Using classical F -statistics and a maximum likelihood method, estimates of the effective population size ( N _e) illustrated a mixture of stable populations with high N _e, and unstable populations characterized by very small N _e estimates (of 5–11 individuals). Owing to human disturbances, intermittent gene flow and genetic drift are likely to be the predominant evolutionary processes shaping the observed genetic structure. However, the practice of multiple matings and sperm storage is likely to provide a reservoir of variability, minimizing the eroding genetic effects of population size reduction and increasing the effective population size.  © 2004 The Linnean Society of London, Biological Journal of the Linnean Society , 2004, 82 , 89–102.     

Population genetic structure of the butterfly Melitaea didyma (Nymphalidae) along a northern distribution range border

The population genetic structure of the butterfly Melitaea didyma was studied along the northern distribution range border in Central Germany by means of allozyme electrophoresis. Individuals were sampled from a total of 21 habitat patches from four regions, and two provinces. Sampling was designed to estimate local vs. regional differentiation. High levels of variability were found, H _e= 0.14–0.21. The mean expected sample heterozygosity from one region, Mosel, was significantly lower than from the Hammelburg region, H _e= 0.17 and 0.19, respectively. Two hierarchical levels of genetic differentiation were found. Within regions individuals sampled from different patches behaved as belonging to one population with high levels of gene flow (Hammelburg F _ST= 0.015, Mosel F _ST= 0.044), though local isolation barriers did create a substructuring of these populations. The inbreeding coefficients, F _IS, were constant over all sample levels, suggesting a similar distribution of habitat patches within regions. Between regions gene flow was limited. An isolation by distance analysis indicated that the hierarchical structure, at the provincial level, may be breaking down due to isolation of regional populations. A more general observation was that the sampling design may greatly have influenced the estimation of genetic differentiation. Depending on which samples were included, overall F _ST estimates ranged from 0.059–0.090.     

Widespread gene flow and high genetic variability in populations of water voles Arvicola terrestris in patchy habitats

Theory predicts that the impact of gene flow on the genetic structure of populations in patchy habitats depends on its scale and the demographic attributes of demes (e.g. local colony sizes and timing of reproduction), but empirical evidence is scarce. We inferred the impact of gene flow on genetic structure among populations of water voles Arvicola terrestris that differed in average colony sizes, population turnover and degree of patchiness. Colonies typically consisted of few reproducing adults and several juveniles. Twelve polymorphic microsatellite DNA loci were examined. Levels of individual genetic variability in all areas were high ( H _O= 0.69–0.78). Assignments of juveniles to parents revealed frequent dispersal over long distances. The populations showed negative F _IS values among juveniles, F _IS values around zero among adults, high F _ST values among colonies for juveniles, and moderate, often insignificant, F _ST values for parents. We inferred that excess heterozygosity within colonies reflected the few individuals dispersing from a large area to form discrete breeding colonies. Thus pre-breeding dispersal followed by rapid reproduction results in a seasonal increase in differentiation due to local family groups. Genetic variation was as high in low-density populations in patchy habitats as in populations in continuous habitats used for comparison. In contrast to most theoretical predictions, we found that populations living in patchy habitats can maintain high levels of genetic variability when only a few adults contribute to breeding in each colony, when the variance of reproductive success among colonies is likely to be low, and when dispersal between colonies exceeds nearest-neighbour distances.     

Genetic differentiation in pikeperch (Sander lucioperca): the relative importance of gene flow, drift and common history

To investigate the relative importance of homogenizing factors, such as gene flow, and diversifying factors, such as drift, genetic variation in pikeperch ( Sander lucioperca ) in two Fennoscandian regions (North and South) was analysed with microsatellites. Allelic richness and the degree of differentiation were significantly higher in the North ( F _ST= 0·20) than in the South ( F _ST= 0·064). In northern areas, assignments of genotypes were almost exclusively to the population of origin, but in southern areas, the proportion of correct assignments was significantly lower. Most samples exhibited significant heterozygote deficits, and the level of relatedness was higher than expected from randomness. These combined results suggest that there has been more gene flow between populations in southern areas than in northern areas, where the importance of genetic drift has been greater. Effective population sizes were small ( c. 100) and did not differ between areas. The effect of a common history appears minor, and thus processes such as genetic drift and gene flow have been more influential in shaping the patterns of genetic diversity in this species.     

Genetic structure of populations of the threatened eastern massasauga rattlesnake, Sistrurus c. catenatus: evidence from microsatellite DNA markers

Throughout its distribution in North America, the threatened eastern massasauga rattlesnake ( Sistrurus c. catenatus ) persists in a series of habitat-isolated disjunct populations of varying size. Here, we use six microsatellite DNA loci to generate information on the degree of genetic differentiation between, and the levels of inbreeding within populations to understand how evolutionary processes operate in these populations and aid the development of conservation plans for this species. Samples were collected from 199 individuals from five populations in Ontario, New York and Ohio. Our results show that all sampled populations: (i) differ significantly in allele frequencies even though some populations are < 50 km apart, and may contain genetically distinct subpopulations < 2 km apart; (ii) have an average of 23% of alleles that are population specific; and (iii) have significant F _IS values (mean overall F _IS= 0.194) probably due to a combination of Wahlund effects resulting from fine-scale genetic differentiation within populations and the presence of null alleles. Our results imply that massasauga populations may be genetically structured on an extremely fine scale even within continuous populations, possibly due to limited dispersal. Additional information is needed to determine if dispersal and mating behaviour within populations can account for this structure and whether the observed differentiation is due to random processes such as drift or to local adaptation. From a conservation perspective, our results imply that these massasauga populations should be managed as demographically independent units and that each has high conservation value in terms of containing unique genetic variation.     

Comparison of quantitative and molecular genetic variation of native vs. invasive populations of purple loosestrife (Lythrum salicaria L., Lythraceae)

Study of adaptive evolutionary changes in populations of invasive species can be advanced through the joint application of quantitative and population genetic methods. Using purple loosestrife as a model system, we investigated the relative roles of natural selection, genetic drift and gene flow in the invasive process by contrasting phenotypical and neutral genetic differentiation among native European and invasive North American populations ( Q _ST −  F _ST analysis). Our results indicate that invasive and native populations harbour comparable levels of amplified fragment length polymorphism variation, a pattern consistent with multiple independent introductions from a diverse European gene pool. However, it was observed that the genetic variation reduced during subsequent invasion, perhaps by founder effects and genetic drift. Comparison of genetically based quantitative trait differentiation ( Q _ST) with its expectation under neutrality ( F _ST) revealed no evidence of disruptive selection ( Q _ST >  F _ST) or stabilizing selection ( Q _ST <  F _ST). One exception was found for only one trait (the number of stems) showing significant sign of stabilizing selection across all populations. This suggests that there are difficulties in distinguishing the effects of nonadaptive population processes and natural selection. Multiple introductions of purple loosestrife may have created a genetic mixture from diverse source populations and increased population genetic diversity, but its link to the adaptive differentiation of invasive North American populations needs further research.     

Genetic diversity and genetic structure of populations of Ranunculus japonicus Thunb. (Ranunculaceae)

Abstract In order to clarify the genetic diversity and population structure of Ranunculus japonicus , allozymic analysis was conducted on 60 populations in southwestern Japan. Considerable genetic variati ons were detected among the populations of R. japonicus . The genetic diversities within species ( H _es = 0.215) and within populations ( H _ep = 0.172) were slightly higher than those of other perennial herbs with widespread distribution and outcrossing plants. Significantly higher values of fixation index were detected in some populations, which might have arisen from restricted mating partners. The majority of genetic variation (approx. 80%) resided within a population and a moderate level of genetic differentiation ( G _ST = 0.203) was observed among populations. The F _ST value (0.203) suggests the existence of a substantial population structure in this species. The highly significant correlation between geographic distance and F _ST values indicates that isolation by distance has played an important role in the construction of the genetic structure of this species.     

Population divergence in the amphicarpic species Amphicarpaea edgeworthii Benth. (Fabaceae): microsatellite markers and leaf morphology

Comparative analyses of the genetic differentiation in microsatellite markers ( F _ST) and leaf morphology characters ( Q _ST) of Amphicarpaea edgeworthii Benth. were conducted to gain insight into the roles of random processes and natural selection in the population divergence. Simple sequence repeat analyses on 498 individuals of 19 natural populations demonstrate that a significant genetic differentiation occurs among populations (mean F _ST = 0.578), and A. edgeworthii is a highly self-fertilized species (mean selfing rate s  = 0.989). The distribution pattern of genetic diversity in this species shows that central populations possess high genetic diversity (e.g. population WL with H _E = 0.673 and population JG with H _E = 0.663), whereas peripheral ones have a low H _E as in population JD (0.011). The morphological divergence of leaf shape was estimated by the elliptical Fourier analysis on the data from 11 natural and four common garden populations. Leaf morphology analyses indicate the morphological divergence does not show strong correlation with the genetic differentiation ( R  = 0.260, P  = 0.069). By comparing the 95% confidence interval of Q _ST with that of F _ST, Q _ST values for five out of 12 quantitative traits are significantly higher than the average F _ST value over eight microsatellite loci. The comparison of F _ST and Q _ST suggests that two kinds of traits can be driven by different evolutionary forces, and the population divergence in leaf morphology is shaped by local selections.  © 2009 The Linnean Society of London, Biological Journal of the Linnean Society , 2009, 96 , 505–516.     

Genetic diversity of a successful colonizer: isolated populations of Metrioptera roeselii regain variation at an unusually rapid rate

Newly founded isolated populations need to overcome detrimental effects of low genetic diversity. The establishment success of a population may therefore depend on various mechanisms such as assortative mating, purging of deleterious alleles, creation of new mutations and/or repeated inflow of new genotypes to reduce the effects of inbreeding and further loss of genetic variation. We compared the level of genetic variation in introduced populations of an insect species (Metrioptera roeselii) far beyond its natural distribution with levels found in their respective founder populations and coupled the data with timing since establishment. This allowed us to analyze if the introduced populations showed signs of temporal changes in genetic variation and have made it possible to evaluate underlying mechanisms. For this, we used neutral genetic markers, seven microsatellite loci and a 676–bp‐long sequence of the mtDNA COI gene. All tested indices (allelic richness, unbiased expected heterozygosity, effective size, haplotype diversity, and nucleotide diversity) except inbreeding coefficient had significantly higher values in populations within the founding populations inside the continuous area of the species distribution compared with the introduced populations. A logarithmic model showed a significant correlation of both allelic richness and unbiased expected heterozygosity with age of the isolated populations. Considering the species' inferred colonization history and likely introduction pathways, we suggest that multiple introductions are the main mechanism behind the temporal pattern observed. However, we argue that influences of assortative mating, directional selection, and effects of an exceptional high intrapopulation mutation rate may have impacts. The ability to regain genetic diversity at this level may be one of the main reasons why M. roeselii successfully continue to colonize northern Europe.     

High genetic diversity in the remnant island population of hihi and the genetic consequences of re-introduction

The maintenance of genetic diversity is thought to be fundamental for the conservation of threatened species. It is therefore important to understand how genetic diversity is affected by the re-introduction of threatened species. We use establishment history and genetic data from the remnant and re-introduced populations of a New Zealand endemic bird, the hihi Notiomystis cincta, to understand genetic diversity loss and quantify the genetic effects of re-introduction. Our data do not support any recent bottleneck events in the remnant population. Furthermore, all genetic diversity measures indicate the remnant hihi population has retained high levels of genetic diversity relative to other New Zealand avifauna with similar histories of decline. Genetic diversity (N(A) , alleles per locus, allelic richness, F(IS) and H(S) ) did not significantly decrease in new hihi populations founded through re-introduction when compared to their source populations, except in the Kapiti Island population (allelic richness and H(S) ) which had very slow post-re-introduction population growth. The N(e) /N(c) ratio in the remnant population was high, but decreased in first-level re-introductions, which together with significant genetic differentiation between populations (F(ST) & Fisher's exact tests) suggest that extant populations are diverging as a result of founder effects and drift. Importantly, simulations of future allele loss predict that the number of alleles lost will be higher in populations with a slow population growth, fewer founding individuals and with nonrandom mating. Interestingly, this species has very high levels of extra-pair paternity which may reduce reproductive variance by allowing social and floater males to reproduce a life history trait that together with a large remnant population size may help maintain higher levels of genetic diversity than expected.     

Understanding population structure and historical demography in a conservation context: population genetics of an endangered fern

The construction of the world's largest hydroelectric scheme across the Yangtze River, the Three Gorges Dams (TGD), in the centre of a southern-central Chinese biodiversity hot spot, the Three Gorges Reservoir Area (TGRA), has attracted international concern and conservation action. To examine whether landscape changes to date have impacted regional flora, and to establish long-term monitoring baselines, we assessed the distribution and dynamics of an endangered and TGRA endemic fern, Adiantum reniforme var. sinense . For eight nuclear microsatellites, high levels of genetic diversity ( H _E = 0.653–0.781) and slightly elevated inbreeding ( F _IS = 0.077–0.197) were found across 13 surveyed populations. The population history of this fern is characterized by a balance of gene flow and genetic drift, where historical dispersal, inferred from coalescent ( F =  0.129) and genetic differentiation ( F _ST = 0.094 and R _ST = 0.180) approaches, is moderate, reflecting an isolation by distance relationship. Importantly, most populations exhibited mutation-drift disequilibrium, suggesting a recent population decline, which is congruent with the known demographic history of the species following dam-related activities. Based on these results, populations of A. reniforme var. sinense are expected to lose genetic diversity and increase genetic structure as dam-related activities decrease size and increase genetic isolation of remnants.     

High Genetic Diversity and Differentiation in Relict Lowland Populations of Gentianella austriaca (A. and J. Kern.) Holub (Gentianaceae)

Abstract: The biennial Gentianella austriaca (A. & J. Kern.) Holub, representing a nutrient-poor grassland taxon of low competition power, is becoming rare in the lowlands of eastern Austria due to changes in land use. To estimate effects of isolation and decreasing population sizes, as well as evolutionary relationships, we investigated variation in isozymes and morphological characters within and between seven populations from the mountains, foothills, and lowlands. Additionally, data on reproduction, habitat, germination and population sizes were collected to examine possible causes of variation and differentiation. We found highest genetic diversity (v_a, v_go) in the lowland and foothill populations, and highest genetic differentiation (D_ja, D_jgo) (i.e., lowest genetic identity: Nei's I) in the lowland populations. The low diversity of the mountain populations might indicate that they are derived from lowland populations. Surprisingly, highest genetic diversity was found in the smallest population. This indicates that in small remnant populations of taxa with a mixed mating system, genetic diversity might be maintained even after many generations after reducing population size dramatically. We found some relationship between genetic diversity and high fitness (germination success) and (inversely) with seed size. Plant size and reproductive success are negatively correlated with altitude, whereas flower size and seed size seem to be subject to other forces of selection. Combining all morphometric, reproductive and genetic traits, the lowland populations are most strongly differentiated and therefore of highest conservation priority.     

Assortative mating among rock-dwelling cichlid fishes supports high estimates of species richness from Lake Malawi

It has been estimated that Lake Malawi, Africa, contains 500–650 endemic species of cichlid fishes, the largest number of vertebrate species endemic to any comparable sized area on the planet. As many of these putative species cannot be distinguished anatomically, these estimates of species richness depend to a great extent on the assumption that sympatrically occurring male colour morphs represent biological species. We have tested this assumption using a combination of behavioural observations of courtship and microsatellite DNA analysis for six putative species of the Pseudotropheus ( Tropheops ) complex and three of the Pseudotropheus ( Maylandia ) complex occurring sympatrically at Nkhata Bay. We were unable to demonstrate assortative courtship for the species pairs Pseudotropheus ( Maylandia ) zebra / P . 'gold zebra' or P. ( Tropheops ) 'band'/ P. ( T .) 'rust' because we were unable to distinguish between the females of these taxa. All other taxa showed clear assortative courtship, except for P. ( T .) 'deep', a deep-water species which was rarely observed. Fixation indices (θ_ST for the infinite allele model, and R _ST for the stepwise mutation model) calculated from six microsatellite DNA loci demonstrated significant deviations from panmixia in all pairwise comparisons of putative species, indicating little or no gene flow between populations. All taxa showed high levels of allelic diversity providing evidence that genetic bottlenecking may have been of limited importance in the speciation process. Assortative mating among taxa differing only in male colouration is supportive of theories that speciation in these fishes has been driven by sexual selection by female choice.     

Extremely low levels of genetic diversity in the terrestrial orchid Epipactis thunbergii (Orchidaceae) in South Korea: implications for conservation

Levels of allozyme variation, population genetic structure, and fine-scale genetic structure (FSGS) of the rare, both sexually and clonally reproducing terrestrial orchid Epipactis thunbergii were examined for eight ( N  = 734) populations in a 20 × 20-km area in South Korea. Twenty-three putative allozyme loci resolved from 15 enzyme systems were used. Extremely low levels of allozyme variation were found within populations: the mean frequency of polymorphic loci was 3.8% [isocitrate dehydrogenase ( Idh-2 ) with two alleles was polymorphic across populations], the mean number of alleles per locus was 1.04, and the mean expected heterozygosity was 0.013. The overall fixation index was not significantly different from zero ( F _IS = 0.069), although the species is self-compatible. However, a significantly high degree of population differentiation was found between populations at Idh-2 ( F _ST = 0.388) in the studied area. Furthermore, spatial autocorrelation analyses revealed a significant FSGS (up to 3 m) within populations. These observations suggest that the main explanatory factors for the extremely low levels of genetic diversity and the shaping of the population genetic structure of E. thunbergii are genetic drift as a result of a small effective population size, a restricted gene flow, and the isolation of populations. Considering the current genetic structure of E. thunbergii , three guidelines are suggested for the development of conservation strategies for the species in South Korea: (1) protection of habitats of standing populations; (2) prohibition by law of any collection of E. thunbergii ; and (3) protection of nearby pollinator populations, given the fact that fruit set in natural habitats is very low.  © 2007 The Linnean Society of London, Botanical Journal of the Linnean Society , 2007, 155 , 161–169.     

Microsatellite analysis reveals high genetic diversity but low genetic structure in Ethiopian indigenous cattle populations

Ethiopian cattle are under threat from uncontrolled mating practices and are at high risk of becoming genetically homogeneous. Therefore, to evaluate genetic diversity, population structure and degree of admixture, 30 microsatellite markers were genotyped using 351 DNA samples from 10 Ethiopian cattle populations and the Holstein breed. The mean number of alleles per cattle population ranged from 6.93 ± 2.12 in Sheko to 7.50 ± 2.35 in Adwa. The mean observed and expected heterozygosities were 0.674 ± 0.015 and 0.726 ± 0.019 respectively. Ethiopian cattle populations have maintained a high level of within-population genetic differentiation (98.7%), the remainder being accounted for by differentiation among populations (1.3%). A highly significant deficiency in heterozygotes was detected within populations ( F _IS = 0.071; P  <   0.001) and total inbreeding ( F _IT = 0.083; P  <   0.001). The study populations were highly admixed but distinct from pure Bos taurus and Bos indicus breeds. The various levels of admixture and high genetic diversity make Ethiopian cattle populations suitable for future genetic improvement and utilization under a wide range of agro-ecologies in Ethiopia.     

Effects of a founder event and supplementary introductions on genetic variation in a captive breeding population of the endangered Spanish killifish

Twelve polymorphic allozyme loci were employed to assess the genetic change in a captive breeding population of the endangered killifish Aphanius baeticus in the Doñana National Park, south‐western Spain. The initial founder event did not significantly reduce the allelic richness or the expected heterozygosity. No genetic bottleneck signature was detected by tests for deviation from mutation‐drift equilibrium. The F _ST between the wild source and captive population, however, was relatively high (0·053 or 0·122 when excluding or including the locus IDHP‐1 * respectively), after just two to three generations in captivity. Two generations after the incorporation of 68 new wild specimens (greater than five generations after founding) decreased the genetic differences and the F _ST(0·041 excluding IDHP‐1 *). The restoration efforts appeared to be helpful and the study of 12 polymorphic loci and a sensitive parameter such as F _ST were useful for monitoring genetic changes in captivity. Nonetheless, future monitoring should include additional highly polymorphic loci (microsatellites) to achieve higher power to detect genetic change. Such restoration and monitoring efforts should help to avoid rapid inbreeding, adaptation to captivity, and to maintain the long‐term evolutionary potential in small isolated populations.     

Long-term stability and effective population size in North Sea and Baltic Sea cod (Gadus morhua)

DNA from archived otoliths was used to explore the temporal stability of the genetic composition of two cod populations, the Moray Firth (North Sea) sampled in 1965 and 2002, and the Bornholm Basin (Baltic Sea) sampled in 1928 and 1997. We found no significant changes in the allele frequencies for the Moray Firth population, while subtle but significant genetic changes over time were detected for the Bornholm Basin population. Estimates of the effective population size ( N _e ) generally exceeded 500 for both populations when employing a number of varieties of the temporal genetic method. However, confidence intervals were very wide and N _e 's most likely range in the thousands. There was no apparent loss of genetic variability and no evidence of a genetic bottleneck for either of the populations. Calculations of the expected levels of genetic variability under different scenarios of N _e showed that the number of alleles commonly reported at microsatellite loci in Atlantic cod is best explained by N _e 's exceeding thousand. Recent fishery-induced bottlenecks can, however, not be ruled out as an explanation for the apparent discrepancy between high levels of variability and recently reported estimates of N _e  << 1000. From life history traits and estimates of survival rates in the wild, we evaluate the compatibility of the species' biology and extremely low N _e / N ratios. Our data suggest that very small N _e 's are not likely to be of general concern for cod populations and, accordingly, most populations do not face any severe threat of losing evolutionary potential due to genetic drift.     

COLOR ASSORTATIVE MATING CONTRIBUTES TO SYMPATRIC DIVERGENCE OF NEOTROPICAL CICHLID FISH

It is still debated vigorously whether sexual selection can result in speciation without physical barriers to gene flow. In this study, we used field data and molecular methods to investigate the gold–normal color polymorphism in two endemic cichlid fish species of crater lake Xiloá, Nicaragua. We found significant assortative mating by color in both Amphilophus xiloaensis and A. sagittae . Focusing on A. xiloaensis , microsatellite allele frequencies, an assignment test, and model-based cluster analysis demonstrates significant and clear genetic differentiation ( F _ST= 0.03) between gold and normal individuals in sympatry. In addition, we find genetic differentiation between all three sympatric and ecologically distinct Midas cichlid species of Lake Xiloá, A. amarillo , A. sagittae , and A. xiloaensis ( F _ST= 0.03 – 0.19), and clear genetic isolation of these species from their closest relative ( A. citrinellus ) in the neighboring great lake Managua. The A. xiloaensis gold morph is genetically more distinct from the lake's other two Midas cichlid species than is A. xiloaensis -normal. Thus, we have identified sexual isolation based on color that is evident in population genetics and mate choice. Our results suggest that sexual selection through color assortative mating may play an important role in incipient sympatric speciation in Midas cichlids of Nicaragua.