Long-term Effective Population Size Research Articles

Population structure of the malaria vector Anopheles moucheti in the equatorial forest region of Africa

BackgroundAnopheles moucheti is a major malaria vector in forested areas of Africa. However, despite its important epidemiological role, it remains poorly known and insufficiently studied. Here, levels of genetic differentiation were estimated between different A. moucheti populations sampled throughout its distribution range in Central Africa.MethodsPolymorphism at ten microsatellite markers was compared in mosquitoes sampled in Cameroon, the Democratic Republic of Congo and an island on Lake Victoria in Uganda. Microsatellite data were used to estimate genetic diversity within populations, their relative long-term effective population size, and the level of genetic differentiation between them.ResultsAll specimens collected in Tsakalakuku (Democratic Republic of Congo) were identified as A. m. bervoetsi while other samples consisted of A. m. moucheti. Successful amplification was obtained at all microsatellite loci within all A. m. moucheti samples while only six loci amplified in A. m. bervoetsi. Allelic richness and heterozygosity were high for all populations except the island population of Uganda and A. m. bervoetsi. High levels of genetic differentiation were recorded between A. m. bervoetsi and each A. m. moucheti sample as well as between the island population of A. m. moucheti and mainland populations. Significant isolation by distance was evidenced between mainland populations.ConclusionHigh levels of genetic differentiation supports complete speciation of A. m. bervoetsi which should henceforth be recognized as a full species and named A. bervoetsi. Isolation by distance is the main force driving differentiation between mainland populations of A. m. moucheti. Genetically and geographically isolated populations exist on Lake Victoria islands, which might serve as relevant field sites for evaluation of innovative vector control strategies.

Population size and genome size in fishes: a closer look

The several thousand-fold range in genome size among animals has remained a subject of active research and debate for more than half a century, but no satisfactory explanation has yet been provided. Many one-dimensional models have been postulated, but so far none has been successful in accounting for observed patterns in genome size diversity. The recent model based on differences in effective population size appeared to gain empirical support with a study of genome size and inferred effective population size in fishes, but there were several questionable aspects of the analysis. First, it was based on an assumption that microsatellite heterozygosity indicates long-term effective population size, whereas in actuality these markers evolve quickly and are sensitive to demographic events. Second, it included both ancient polyploids and non-polyploids, the former of which did not gain their current genome sizes through the accumulation of slightly deleterious mutations as required in the model. Third, the analysis neglected the tremendous influence that Pleistocene glaciation bottlenecks had on heterozygosities in freshwater (and far less so, marine) fishes. In sum, it is apparent that genomes reached their current sizes in most fishes long before contemporary microsatellite heterozygosities were shaped, and that ancient polyploidy rather than the accumulation of mildly deleterious transposon insertions in small populations is the dominant factor that has influenced the large end of the range of genome sizes among fishes.

More genes underwent positive selection in chimpanzee evolution than in human evolution

Observations of numerous dramatic and presumably adaptive phenotypic modifications during human evolution prompt the common belief that more genes have undergone positive Darwinian selection in the human lineage than in the chimpanzee lineage since their evolutionary divergence 6-7 million years ago. Here, we test this hypothesis by analyzing nearly 14,000 genes of humans and chimps. To ensure an accurate and unbiased comparison, we select a proper outgroup, avoid sequencing errors, and verify statistical methods. Our results show that the number of positively selected genes is substantially smaller in humans than in chimps, despite a generally higher nonsynonymous substitution rate in humans. These observations are explainable by the reduced efficacy of natural selection in humans because of their smaller long-term effective population size but refute the anthropocentric view that a grand enhancement in Darwinian selection underlies human origins. Although human and chimp positively selected genes have different molecular functions and participate in different biological processes, the differences do not ostensibly correspond to the widely assumed adaptations of these species, suggesting how little is currently known about which traits have been under positive selection. Our analysis of the identified positively selected genes lends support to the association between human Mendelian diseases and past adaptations but provides no evidence for either the chromosomal speciation hypothesis or the widespread brain-gene acceleration hypothesis of human origins.

Isolation and interchange among insular spinner dolphin communities in the South Pacific revealed by individual identification and genetic diversity

Gray's spinner dolphins Stenella longirostris longirostris are found in apparently rela- tively small and discrete communities around many islands throughout the Pacific. However, the boundaries of these communities, on the scale of a dolphin's lifespan or across generations, are unknown. Here we report a combined demographic and genetic approach to describing the isolation and interchange of insular spinner dolphins among island communities of the Society Archipelago, French Polynesia. Dorsal fin photographs for individual identification and biospy samples for genetic analyses (n = 154) were collected from 6 island communities during 189 small-boat surveys over 3 yr. Capture-recapture analyses at Moorea (our primary study site), based on long-term observations of distinctively marked individuals and microsatellite genotypes (12 loci), indicated a local community of about 150 dolphins. This community appeared relatively closed on an intra-generational scale, as confirmed by resightings of individuals across 15 yr. Surveys around neighbouring islands indicated the presence of similar distinct communities, likely to follow demographic patterns similar to Moorea, with relatively low levels of interchange between communities. Overall, significant differentiation at both mitochondrial and nuclear levels indicated restricted gene flow among neighbouring communi- ties, although some individual movement was documented. High levels of insular mitochondrial DNA (mDNA) genetic diversity (female long-term effective population size (Nef) ~ 100 000) contrasted with demographic characteristics. No evidence of bottlenecks was found in microsatellite allele frequen- cies or mtDNA haplotypes, discounting the possibility of a recent founder effect. Instead, we suggest that this genetic pattern is the result of metapopulation structure, based on numerous insular commu- nities evolutionarily connected through male and female gene flow.

Coding sequence polymorphism in avian mitochondrial genomes reflects population histories

Nucleotide sequence diversity at mitochondrial protein-coding loci from 72 species of birds from different geographical regions was analysed in order to test the hypothesis that temperate zone species show population genetic effects of past glaciation. Temperate zone species showed reduced nucleotide diversity in comparison to tropical mainland species, suggesting that the latter have long-term effective population sizes due to population bottleneck effects during the most recent glaciation. This hypothesis was further supported by evidence of an unusually high estimated rate of population growth in species breeding in North America and wintering in the New World tropics (Nearctic migrants), consistent with population expansion after a bottleneck. Nearctic migrants also showed evidence of an abundance of rare nonsynonymous (amino acid-altering) polymorphisms, a pattern suggesting that slightly deleterious polymorphisms drifted to high frequencies during a bottleneck and are now being eliminated by selection. Because the shape of the North American land mass limited the area available for refugia during glaciation, the bottleneck effects are predicted to have been particularly strong in Nearctic migrants, and this prediction was supported. The reduced genetic diversity of Nearctic migrants provides an additional basis for concern for the survival of these species, which are threatened by loss of habitat in the winter range and by introduced disease.

Contrasting patterns of population subdivision and historical demography in three western Mediterranean lizard species inferred from mitochondrial DNA variation

Pleistocene climatic oscillations were a major force shaping genetic variability in many taxa. We analyse the relative effects of the ice ages across a latitudinal gradient in the Western Mediterranean region, testing two main predictions: (i) species with historical distributions in northern latitudes should have experienced greater loss of suitable habitat, resulting in higher extinction of historical lineages than species distributed in southern latitudes, where the effects of the ice ages were not as drastic. This would be reflected in the observation of lower diversity and number of differentiated lineages in northern areas. (ii) a signature of demographic expansion following the climate amelioration should be obvious in northern species, whereas in the south evidence of long-term effective population size stability should be observed. We used as models three species of wall lizards (Podarcis bocagei, Podarcis carbonelli and Podarcis vaucheri) that replace each other along the study area. We investigated the patterns of mitochondrial DNA diversity and subdivision and obtained demographic parameter estimates for each species. Our results suggest that P. bocagei, the northernmost species, bears low genetic diversity, a shallow coalescent history and marks of a demographic expansion. In contrast, P. vaucheri, the species with a southernmost distribution, shows deeper coalescence events, complex geographical substructure and no evidence for population growth. The species with an intermediate distribution, P. carbonelli, shows average levels of diversity, substructure and population growth. Taken together, these results conform to our main predictions and are explained by a differential influence of the ice ages on distinct latitudes.

Inbreeding levels and prey abundance interact to determine fecundity in natural populations of two species of wolf spider

Long-term effective population size is expected, and has been shown, to correlate positively with various measures of population fitness. Here we examine the interacting effects of population size (as a surrogate for genetic factors) and prey consumption rates (as a surrogate for environmental quality) on fecundity in two sympatric species of wolf spider, Rabidosa punctulata and Rabidosa rabida. Population size was estimated in each of seven genetically isolated populations in each of 3 years using mark-recapture methods. Fecundity was estimated as the mean number of live offspring produced by ∼15 females sampled from each population of each species each year for 3 years. Prey consumption rates were estimated by sampling ∼300 spiders per population per year and assaying the proportion of spiders with prey. Larger populations have higher fecundity and more genetic diversity than smaller populations. Variation among populations in fecundity for a given year could be attributed most strongly to differences in population size, with variation in prey consumption rates and the interaction between population size and prey consumption playing smaller but still important roles. During the most stressful environmental conditions, the smallest populations of both species experienced disproportionately low-fecundity rates, more than doubling the estimated number of lethal equivalents during those years. The evidence presented in this paper for inbreeding-environment interactions at the population level and further evidence for a log-linear relationship between population size and fitness have important implications for conservation.

The genetic effective and adult census size of an Australian population of tiger prawns (Penaeus esculentus)

This study compares estimates of the census size of the spawning population with genetic estimates of effective current and long-term population size for an abundant and commercially important marine invertebrate, the brown tiger prawn (Penaeus esculentus). Our aim was to focus on the relationship between genetic effective and census size that may provide a source of information for viability analyses of naturally occurring populations. Samples were taken in 2001, 2002 and 2003 from a population on the east coast of Australia and temporal allelic variation was measured at eight polymorphic microsatellite loci. Moments-based and maximum-likelihood estimates of current genetic effective population size ranged from 797 to 1304. The mean long-term genetic effective population size was 9968. Although small for a large population, the effective population size estimates were above the threshold where genetic diversity is lost at neutral alleles through drift or inbreeding. Simulation studies correctly predicted that under these experimental conditions the genetic estimates would have non-infinite upper confidence limits and revealed they might be overestimates of the true size. We also show that estimates of mortality and variance in family size may be derived from data on average fecundity, current genetic effective and census spawning population size, assuming effective population size is equivalent to the number of breeders. This work confirms that it is feasible to obtain accurate estimates of current genetic effective population size for abundant Type III species using existing genetic marker technology.

Nucleotide Variation and Haplotype Diversity in a 10-kb Noncoding Region in Three Continental Human Populations

Noncoding regions are usually less subject to natural selection than coding regions and so may be more useful for studying human evolution. The recent surveys of worldwide DNA variation in four 10-kb noncoding regions revealed many interesting but also some incongruent patterns. Here we studied another 10-kb noncoding region, which is in 6p22. Sixty-six single-nucleotide polymorphisms were found among the 122 worldwide human sequences, resulting in 46 genotypes, from which 48 haplotypes were inferred. The distribution patterns of DNA variation, genotypes, and haplotypes suggest rapid population expansion in relatively recent times. The levels of polymorphism within human populations and divergence between humans and chimpanzees at this locus were generally similar to those for the other four noncoding regions. Fu and Li's tests rejected the neutrality assumption in the total sample and in the African sample but Tajima's test did not reject neutrality. A detailed examination of the contributions of various types of mutations to the parameters used in the neutrality tests clarified the discrepancy between these test results. The age estimates suggest a relatively young history in this region. Combining three autosomal noncoding regions, we estimated the long-term effective population size of humans to be 11,000 +/- 2800 using Tajima's estimator and 17,600 +/- 4700 using Watterson's estimator and the age of the most recent common ancestor to be 860,000 +/- 258,000 years ago.

Genetic population structure as indirect measure of dispersal ability in a Lake Tanganyika cichlid

Diversification and speciation processes are influenced by intrinsic (ecological specialization, dispersal) and extrinsic (habitat structure and instability) factors, but the effect of ecological characteristics on dispersal is difficult to assess. This study uses mitochondrial control region sequences to investigate the population structure and demographic history of the endemic Lake Tanganyika cichlid Neolamprologus caudopunctatus with a preference for the rock-sand interface along two stretches of continuous, rocky shoreline, and across a sandy bay representing a potential dispersal barrier. Populations along uninterrupted habitat were not differentiated; whereas, the sandy bay separated two reciprocally monophyletic clades. The split between the two clades between 170,000 and 260,000 years BP coincides with a period of rising water level following a major lowstand, and indicates that clades remained isolated throughout subsequent lake level fluctuations. Low long-term effective population sizes were inferred from modest genetic diversity estimates, and may be due to recent population expansions starting from small population sizes 45,000-60,000 years BP. Comparisons with available data from specialized rock-dwelling species of the same area suggest that habitat structure and lake level fluctuations determine phylogeographic patterns on large scales, while fine-scale population structure and demography are modulated by species-specific ecologies.

Different Post-Pleistocene Histories of Eurasian Parids

Previous phylogeographic studies of the great tit (Parus major) and the willow tit (Parus montanus) found a general absence of phylogeographic structure for both species and suggested that each species underwent range contraction during the last Ice Age and survived in relatively low numbers, P. major in southern Europe and P. montanus in southeastern Asia. However, prior studies did not sample the entire range of either species. We analyzed sequence data for the complete mitochondrial ND2 gene from 87 P. major and 139 P. montanus from 15 new Eurasian localities, both to test prior conclusions and to provide better coverage of each species' range. Our analyses confirmed the absence of phylogeographic structure in P. major and P. montanus and supported the prior refuge hypothesis for P. major. For P. montanus, we concluded that besides surviving the Ice Age in southeastern Asia, as previously hypothesized, it apparently sustained a relatively large population in northern Eurasian riverine thickets and then expanded eastward. Genetic diversity was low in P. major (pi = 0.0012, h = 0.64) and moderate in P. montanus (pi = 0.0021, h = 0.88), suggesting higher long-term effective population sizes and the older ages of populations in P. montanus. If molecular substitution rates are similar, P. montanus colonized its current Eurasian range earlier than P. major. Differences between prior studies and ours likely result from sampling gaps in earlier studies.

DNA microsatellite characterization of the jaguar (Panthera onca) in Colombia

The Colombian jaguar population is thought to contain two different subspecies, Panthera onca centralis and Panthera onca onca. The genetic structure of this population was evaluated using 12 microsatellite loci (n = 62 samples). In addition, 22 jaguar DNA samples from Guatemala, Paraguay, Perú, Bolivia, Venezuela and Brazil were analyzed for these microsatellite loci (n = 84 samples). The results of this study indicate six primary themes. First, the levels of gene diversity were very high. Second, the majority of the loci analyzed showed an absence of Hardy-Weinberg equilibrium, probably due to the Wahlund effect (= population subdivision). Third, several microsatellite loci showed significant heterogeneity between the two supposed subspecies in the country. Nevertheless, gene flow was present between them, and heterogeneity was relatively low, although the assignment analyses showed good classification of the jaguars studied into their respective subspecies. Fourth, the long-term historical effective population sizes were calculated through a maximum likelihood procedure for single and multi-step mutation models. Fifth, seven out of twelve DNA microsatellites studied significantly deviated from a single-step mutation model. However, the overall mean multi-step mutation percentage for these 12 DNA microsatellites was only 6%. Therefore, 94% of mutations were uni-step. Sixth, no bottleneck events were detected in the Colombian jaguar population overall.

Conflicting patterns of mitochondrial and nuclear DNA diversity in Phylloscopus warblers

Molecular variation is often used to infer the demographic history of species, but sometimes the complexity of species history can make such inference difficult. The willow warbler, Phylloscopus trochilus, shows substantially less geographical variation than the chiffchaff, Phylloscopus collybita, both in morphology and in mitochondrial DNA (mtDNA) divergence. We therefore predicted that the willow warbler should harbour less nuclear DNA diversity than the chiffchaff. We analysed sequence data obtained from multiple samples of willow warblers and chiffchaffs for the mtDNA cytochrome b gene and four nuclear genes. We confirmed that the mtDNA diversity among willow warblers is low (pi = 0.0021). Sequence data from three nuclear genes (CHD-Z, AFLP-WW1 and MC1R) not linked to the mitochondria demonstrated unexpectedly high nucleotide diversity (pi values of 0.0172, 0.0141 and 0.0038) in the willow warbler, on average higher than the nucleotide diversity for the chiffchaff (pi values of 0.0025, 0.0017 and 0.0139). In willow warblers, Tajima's D analyses showed that the mtDNA diversity, but not the nuclear DNA diversity, has been reduced relative to the neutral expectation of molecular evolution, suggesting the action of a selective sweep affecting the maternally inherited genes. The large nuclear diversity seen within willow warblers is not compatible with processes of neutral evolution occurring in a population with a constant population size, unless the long-term effective population size has been very large (N(e) > 10(6)). We suggest that the contrasting patterns of genetic diversity in the willow warbler may reflect a more complex evolutionary history, possibly including historical demographic fluctuations or historical male-biased introgression of nuclear genes from a differentiated population of Phylloscopus warblers.

Estimation of Long-Term Effective Population Sizes Through the History of Durum Wheat Using Microsatellite Data

Estimation of long-term effective population size (N(e)) from polymorphism data alone requires an independent knowledge of mutation rate. Microsatellites provide the opportunity to estimate N(e) because their high mutation rate can be estimated from observed mutations. We used this property to estimate N(e) in allotetraploid wheat Triticum turgidum at four stages of its history since its domestication. We estimated the mutation rate of 30 microsatellite loci. Allele-specific mutation rates mu were predicted from the number of repeats of the alleles. Effective population sizes were calculated from the diversity parameter theta = 4N(e)mu. We demonstrated from simulations that the unbiased estimator of theta based on Nei's heterozygosity is the most appropriate for estimating N(e) because of a small variance and a relative robustness to variations in the mutation model compared to other estimators. We found a N(e) of 32,500 individuals with a 95% confidence interval of [20,739; 45,991] in the wild ancestor of wheat, 12,000 ([5790; 19,300]) in the domesticated form, 6000 ([2831; 9556]) in landraces, and 1300 ([689; 2031]) in recent improved varieties. This decrease illustrates the successive bottlenecks in durum wheat. No selective effect was detected on our loci, despite a complete loss of polymorphism for two of them.

A separation-of-timescales approach to the coalescent in a continuous population.

This article presents an analysis of a model of isolation by distance in a continuous, two-dimensional habitat. An approximate expression is derived for the distribution of coalescence times for a pair of sequences sampled from specific locations in a rectangular habitat. Results are qualitatively similar to previous analyses of isolation by distance, but account explicitly for the location of samples relative to the habitat boundaries. A separation-of-timescales approach takes advantage of the fact that the sampling locations affect only the recent coalescent behavior. When the population size is larger than the number of generations required for a lineage to cross the habitat range, the long-term genealogical process is reasonably well described by Kingman's coalescent with time rescaled by the effective population size. This long-term effective population size is affected by the local dispersal behavior as well as the geometry of the habitat. When the population size is smaller than the time required to cross the habitat, deep branches in the genealogy are longer than would be expected under the standard neutral coalescent, similar to the pattern expected for a panmictic population whose population size was larger in the past.

Colonization history of the high‐arctic pink‐footed gooseAnser brachyrhynchus

Population structure and phylogeography of the pink-footed goose, Anser brachyrhynchus Baillon 1833, was studied using mtDNA control region sequences (221 bp) from 142 individuals. Present breeding areas of the species in Greenland, Iceland, and Svalbard were largely covered by ice during the late Pleistocene. In pairwise comparisons phiST estimates showed significant differentiation among eastern and western populations, whereas sampling localities within both areas were not differentiated. The mtDNA data indicate that the populations have separated recently (less than 10 000 years ago) and present breeding areas were colonized from one refugial population. The levels of haplotype and nucleotide diversity were approximately five times higher for the eastern population compared to the western population and suggest that the latter was colonized by a subset of eastern birds. Time to the most recent common ancestor of the species is 32 000-46 000 years, i.e. the present mtDNA variation of the pink-footed goose has accumulated during the last 0.1 My. Estimates of the long-term female effective population size (5400-7700 for the eastern population) imply that the refugial population of the pink-footed goose has been large. Tundra habitats were more extensive in cold periods of the late Pleistocene than today and may have sustained population sizes that allowed the accumulation of extant genetic polymorphism. It is not probable that the postulated small refugial areas in the high latitudes had a significant role in maintaining this diversity.

Conservation genetics in aquatic species: General approaches and case studies in fishes and springsnails of arid lands

Many aquatic species have become endangered because of the elimination and fragmentation of their habitat, human alteration of the environment, and introduced species. New applications and insights from genetic studies of endangered species are being used to formulate management plans to prevent extinction of endangered taxa. Here we discuss three applications of genetic techniques to the conservation of endangered aquatic species in the desert Southwest, USA. First, we discuss estimation of long-term effective population size in three endangered big-river fishes, bonytail chub (Gila elegans), humpback chub (G. cypha) and razorback suckers (Xyrauchen texanus), of the lower Colorado River. For all three species, the current census number is much smaller than the estimated effective population size in which these species evolved. Second, we discuss the determination of species, ESUs, and MUs in the endangered Sonoran topminnow. Molecular genetic data show that the Gila and Yaqui topminnows are different species. Experimental examination of pre-mating and post-mating reproductive isolation demonstrates that some reproductive isolation has already developed between these two species. Finally, springsnails (Pyrgulopsis) of the lower Colorado River basin exist only in remote and isolated springs. Examination of molecular variation generally supports the morphological designation of 16 different species, although a few species have several ESUs. There does not appear to be a correlation of genetic and geographic distance between species, suggesting that they are quite old. These genetic studies provide insight into the conservation of these rare aquatic species. Although mitigating ecological factors may be most significant in preventing extinction, genetic studies can provide the evolutionary context for endangered species and identify what units are most in need of conservation.

Population genetics of the swimming crab Callinectes bellicosus (Brachyura: Portunidae) from the eastern Pacific Ocean

The population genetics and historical demography of the swimming crab Callinectes bellicosus from the eastern Pacific were assessed using mitochondrial DNA (mtDNA) sequences from portions of the cytochrome c oxidase subunit I (COI) and cytochrome b (Cytb) genes. Analysis of molecular variance of sequence data from crabs collected from nine localities, ranging from the upper to lower Gulf of California and the outer coast of the Baja California peninsula, revealed an absence of population structure, suggesting a high level of gene flow over a wide geographic area. Maximum-likelihood estimates of long-term effective population size obtained with the program FLUCTUATE, in addition to highly significant values obtained from neutrality tests (Tajima’s D, Fu and Li’s D, and Fu’s FS) and mismatch analysis, are consistent with a population expansion dating to the Pleistocene epoch. Phylogenetic analysis of C. bellicosus sequences using both neighbor-joining and Bayesian methods revealed a widely distributed subclade (clade II) cryptically embedded at low frequency in the main (clade I) population. Although sequence divergence between the two clades was low (1.1% COI; 0.6% Cytb), statistical support for the split was high. The Kimura-2-parameter genetic distance between C. bellicosus and the sympatric and morphologically similar C. arcuatus was high (d=0.17) and similar to the genetic distance between C. bellicosus and the allopatric C. sapidus from the western Atlantic (d=0.18), suggesting an ancient (Miocene) divergence of C. bellicosus and C. arcuatus.

Recent developments in conservation genetics

Because of the new molecular genetics information from various genome projects, new applications and insights for genetics studies in endangered species are forthcoming. Neutral variants are generally used for conservation applications and estimating evolutionary parameters, and with highly variable loci, many more markers, or extensive sequence data, these approaches should become much more informative. Detrimental and adaptive variation is of importance in conservation genetics but identification and characterization of such variation is more difficult. Neutral variants might be used to identify adaptive variants but the overlay of different mutational processes and selective regimes suggests that great caution should be used in making such predictions. Examples in endangered species discussed below are (1) estimation of long-term effective population size using neutral markers in three fishes of the lower Colorado River; (2) genetic restoration of the Florida panther by the introduction of Texas cougars; (3) impact of pathogens and adaptive variation on the winter-run chinook salmon; (4) examination of the adaptive variation for a major histocompatibility complex (MHC) gene in red wolves. Although these examples are all from animals, the techniques and approaches used should prove equally useful in endangered plant species.

Monitoring the genetic diversity of three Asian arowana ( Scleropages formosus) captive stocks using AFLP and microsatellites

Asian arowana ( Scleropages formosus) is listed by CITES as a highly endangered fish species. The genetic diversity and population structure of three Asian arowana captive stocks was analyzed using microsatellites and AFLP markers on 32 randomly collected individuals from each stock. Six AFLP primer pairs amplified a total of 324 bands across the three sample sets. The green stock showed the highest percentage of polymorphic bands (15.6%) and expected heterozygosity (0.26), followed by the red (13.0% and 0.24), then the golden one (12.7% and 0.22). Microsatellite analysis showed high allelic and gene diversity in all three varieties. The green stock showed higher allele number (100) and higher gene diversity (0.75) than the red (98 and 0.74) and golden ones (85 and 0.71), respectively. The estimates of long-term effective population size by two different methods ranged from 4741 to 7288 and from 25,056 to 64,641, respectively. Tests of heterozygosity and allele frequency distribution indicate that a recent bottleneck in any of the three captive populations is highly unlikely. Altogether, the data suggest that the captive breeding program for conservation and sustainable use of the endangered Asian arowana will be successful. In addition, we found that the differentiation among populations was intermediate ( F ST=0.047, R ST=0.103) and the genetic distance was the smallest between the green and red, whereas the largest between the red and golden varieties.