Population structure within and between subspecies of the Mediterranean triplefin fish Tripterygion delaisi revealed by highly polymorphic microsatellite loci

Although F(ST) values are widely used to elucidate population relationships, in some cases, when employing highly polymorphic loci, they should be regarded with caution, particularly when subspecies are under consideration. Tripterygion delaisi presents two subspecies that were investigated here, using 10 microsatellite loci. A Bayesian approach allowed us to clearly identify both subspecies as two different evolutionary significant units. However, low F(ST) values were found between subspecies as a consequence of the large number of alleles per locus, while homoplasy could be disregarded as indicated by the standardized genetic distance G'(ST). Heterozygosity saturation was observed in highly polymorphic loci containing more than 15 alleles, and this threshold was used to define two loci pools. The less variable loci pool revealed higher genetic variance between subspecies, while the more variable pool showed higher genetic variance between populations. Furthermore, higher differentiation was also observed between populations using G'(ST) with the more variable loci. Nonetheless, a more reliable population structure within subspecies was obtained when all loci were included in the analyses. In T. d. xanthosoma, isolation by distance was detected between the eight analysed populations, and six genetically homogeneous clusters were inferred by Bayesian analyses that are in accordance with F(ST) values. The neighbourhood-size method also indicated rather small dispersal capabilities. In conclusion, in fish with limited adult and larval dispersal capabilities, continuous rocky habitat seems to allow contact between populations and prevent genetic differentiation, while large discontinuities of sand or deep-water channels seems to reduce gene flow.     

Patterns of differentiation in a colour polymorphism and in neutral markers reveal rapid genetic changes in natural damselfly populations

The existence and mode of selection operating on heritable adaptive traits can be inferred by comparing population differentiation in neutral genetic variation between populations (often using F(ST) values) with the corresponding estimates for adaptive traits. Such comparisons indicate if selection acts in a diversifying way between populations, in which case differentiation in selected traits is expected to exceed differentiation in neutral markers [F(ST )(selected) > F(ST )(neutral)], or if negative frequency-dependent selection maintains genetic polymorphisms and pulls populations towards a common stable equilibrium [F(ST) (selected) < F(ST) (neutral)]. Here, we compared F(ST) values for putatively neutral data (obtained using amplified fragment length polymorphism) with estimates of differentiation in morph frequencies in the colour-polymorphic damselfly Ischnura elegans. We found that in the first year (2000), population differentiation in morph frequencies was significantly greater than differentiation in neutral loci, while in 2002 (only 2 years and 2 generations later), population differentiation in morph frequencies had decreased to a level significantly lower than differentiation in neutral loci. Genetic drift as an explanation for population differentiation in morph frequencies could thus be rejected in both years. These results indicate that the type and/or strength of selection on morph frequencies in this system can change substantially between years. We suggest that an approach to a common equilibrium morph frequency across all populations, driven by negative frequency-dependent selection, is the cause of these temporal changes. We conclude that inferences about selection obtained by comparing F(ST) values from neutral and adaptive genetic variation are most useful when spatial and temporal data are available from several populations and time points and when such information is combined with other ecological sources of data.     

Microsatellite allele sizes: a simple test to assess their significance on genetic differentiation

The mutation process at microsatellite loci typically occurs at high rates and with stepwise changes in allele sizes, features that may introduce bias when using classical measures of population differentiation based on allele identity (e.g., F(ST), Nei's Ds genetic distance). Allele size-based measures of differentiation, assuming a stepwise mutation process [e.g., Slatkin's R(ST), Goldstein et al.'s (deltamu)(2)], may better reflect differentiation at microsatellite loci, but they suffer high sampling variance. The relative efficiency of allele size- vs. allele identity-based statistics depends on the relative contributions of mutations vs. drift to population differentiation. We present a simple test based on a randomization procedure of allele sizes to determine whether stepwise-like mutations contributed to genetic differentiation. This test can be applied to any microsatellite data set designed to assess population differentiation and can be interpreted as testing whether F(ST) = R(ST). Computer simulations show that the test efficiently identifies which of F(ST) or R(ST) estimates has the lowest mean square error. A significant test, implying that R(ST) performs better than F(ST), is obtained when the mutation rate, mu, for a stepwise mutation process is (a) >/= m in an island model (m being the migration rate among populations) or (b) >/= 1/t in the case of isolated populations (t being the number of generations since population divergence). The test also informs on the efficiency of other statistics used in phylogenetical reconstruction [e.g., Ds and (deltamu)(2)], a nonsignificant test meaning that allele identity-based statistics perform better than allele size-based ones. This test can also provide insights into the evolutionary history of populations, revealing, for example, phylogeographic patterns, as illustrated by applying it on three published data sets.     

Allozyme and microsatellite loci provide discordant estimates of population differentiation in the endangered dusky grouper (Epinephelus marginatus) within the Mediterranean Sea

The dusky grouper, Epinephelus marginatus, inhabits coastal reefs in the Mediterranean Sea and Atlantic Ocean. A decline in the abundance of this long-lived protogynous hermaphrodite has led to its listing as an endangered species in the Mediterranean, and heightened management concerns regarding its genetic variability and population substructure. To address these concerns, we analysed genetic variation at seven microsatellite and 28 allozyme loci in dusky groupers sampled from seven areas (for microsatellites) and three areas (for allozymes) in the west-central Mediterranean. Levels of genetic variability were higher for microsatellites than for allozymes (mean H(E) = 0.78 and 0.07, respectively), but similar to those observed in other marine fishes with comparable markers. Both microsatellites and allozymes revealed significant genetic differentiation among all areas analysed with each class of marker, but the magnitude of differentiation revealed by allozymes over three locales (F(ST) = 0.214) was greater than that detected with microsatellites over seven areas, or over the three areas shared with the allozyme analysis (F(ST) = 0.018 and approximately 0, respectively). A large proportion of the allozyme differentiation was due to a single locus (ADA*) possibly influenced by selection, but allozyme differentiation over the three areas was still highly significant (F(ST) = 0.06, P < 0.0001), and the 95% confidence intervals for allozyme and microsatellite F(ST) did not overlap when this locus was excluded. There was no evidence of isolation by distance with either class of markers. Our results lead us to conclude that dusky groupers are not panmictic in the Mediterranean Sea and suggest that they should be managed on a local basis. However, more work is needed to elucidate genetic relationships among populations.     

Genome scan to detect genetic structure and adaptive genes of natural populations of Cryptomeria japonica

We investigated 29 natural populations of Cryptomeria japonica using 148 cleaved amplified polymorphic sequence markers to elucidate their genetic structure and identify candidate adaptive genes of this species. In accordance with the inferred evolutionary history of the species during and after the last glacial episode, the genetic diversity was higher in western populations than in northern populations. The results of phylogenetic and genetic structure analyses suggest that populations of the two main varieties of the species have clearly diverged from each other and that two of the examined loci are strongly associated with the differentiation between the two varieties. Using a coalescent simulation based on F(ST) and H(e) values, we detected five genes that had higher, and two that had lower, values than the respective 99% confidence intervals (C.I.s) that are theoretically expected intervals under a neutral infinite-island model. We also detected 13 outlier loci using a coalescent simulation based on the assumption that the 2 varieties originated from the splitting of an ancestral population. Four of these loci were detected by both methods, two of which were detected in a genetic structure analysis as loci associated with differentiation between the two varieties of the species, and are strong candidates for genes that have been subject to selection.     

Efficiency of model-based Bayesian methods for detecting hybrid individuals under different hybridization scenarios and with different numbers of loci

Accurate detection of offspring resulting from hybridization between individuals of distinct populations has a range of applications in conservation and population genetics. We assessed the hybrid identification efficiency of two methods (implemented in the STRUCTURE and NEWHYBRIDS programs) which are tailored to identifying hybrid individuals but use different approaches. Simulated first- and second-generation hybrids were used to assess the performance of these two methods in detecting recent hybridization under scenarios with different levels of genetic divergence and varying numbers of loci. Despite the different approaches of the methods, the hybrid detection efficiency was generally similar and neither of the two methods outperformed the other in all scenarios assessed. Interestingly, hybrid detection efficiency was only minimally affected by whether reference population allele frequency information was included or not. In terms of genotyping effort, efficient detection of F1 hybrid individuals requires the use of 12 or 24 loci with pairwise F(ST) between hybridizing parental populations of 0.21 or 0.12, respectively. While achievable, these locus numbers are nevertheless higher than the number of loci currently commonly applied in population genetic studies. The method of STRUCTURE seemed to be less sensitive to the proportion of hybrids included in the sample, while NEWHYBRIDS seemed to perform slightly better when individuals from both backcross and F1 hybrid classes were present in the sample. However, separating backcrosses from purebred parental individuals requires a considerable genotyping effort (at least 48 loci), even when divergence between parental populations is high.     

Measuring population differentiation using GST or D? A simulation study with microsatellite DNA markers under a finite island model and nonequilibrium conditions

The genetic differentiation of populations is a key parameter in population genetic investigations. Wright's F(ST) (and its relatives such as G(ST) ) has been a standard measure of differentiation. However, the deficiencies of these indexes have been increasingly realized in recent years, leading to some new measures being proposed, such as Jost's D (Molecular Ecology, 2008; 17, 4015). The existence of these new metrics has stimulated considerable debate and induced some confusion on which statistics should be used for estimating population differentiation. Here, we report a simulation study with neutral microsatellite DNA loci under a finite island model to compare the performance of G(ST) and D, particularly under nonequilibrium conditions. Our results suggest that there exist fundamental differences between the two statistics, and neither G(ST) nor D operates satisfactorily in all situations for quantifying differentiation. D is very sensitive to mutation models but G(ST) noticeably less so, which limits D's utility in population parameter estimation and comparisons across genetic markers. Also, the initial heterozygosity of the starting populations has some important effects on both the individual behaviours of G(ST) and D and their relative behaviours in early differentiation, and this effect is much greater for D than G(ST) . In the early stages of differentiation, when initial heterozygosity is relatively low (<0.5, if the number of subpopulations is large), G(ST) increases faster than D; the opposite is true when initial heterozygosity is high. Therefore, the state of the ancestral population appears to have some lasting impacts on population differentiation. In general, G(ST) can measure differentiation fairly well when heterozygosity is low whatever the causes; however, when heterozygosity is high (e.g. as a result of either high mutation rate or high initial heterozygosity) and gene flow is moderate to strong, G(ST) fails to measure differentiation. Interestingly, when population size is not very small (e.g. N ≥ 1000), G(ST) measures differentiation quite linearly with time over a long duration when gene flow is absent or very weak even if mutation rate is not low (e.g. μ = 0.001). In contrast, D, as a differentiation measure, performs rather robustly in all these situations. In practice, both indexes should be calculated and the relative levels of heterozygosities (especially H(S) ) and gene flow taken into account. We suggest that a comparison of the two indexes can generate useful insights into the evolutionary processes that influence population differentiation.     

Isolation and characterization of 16 novel microsatellite loci in two inbred strains of the Chinese hamster (Cricetulus griseus)

The Shanyi inbred A and E strains of the Chinese hamster are widely used in biomedical research, but detailed genetic characterization has been lacking. We developed microsatellite markers that could be used for genetic diversity analysis and linkage map construction. We isolated and characterized 16 novel microsatellite loci from a microsatellite-enriched genomic DNA library. These loci were genotyped in 48 animals from the two strains, and the polymorphic information content was determined. In the Shanyi A and E populations, 14 and 15 loci were found to be polymorphic, respectively, with polymorphic information content ranging from 0.1393 to 0.8082 and from 0.1109 to 0.7397, respectively. A total of 115 alleles were found for the 16 microsatellite loci in the two populations; the mean observed heterozygosity (H(O)) was 0.5191 and 0.4333 for the A and E populations, respectively, indicating marked genetic variation within the two populations. Correspondingly, the F(ST) values ranged from 0.002 to 0.9253, with an overall mean of 0.1935, indicating significant genetic difference between the two strains. The population differentiation levels were substantiated by Nei's genetic distance and full Bayesian analyses computed with STRUCTURE. Despite the genetic diversity and differentiation within and between the two inbred populations, the 48 individuals were correctly allocated into their original populations with high statistical confidence based on these 16 microsatellite loci. These novel microsatellite loci should be useful genetic markers for these two Chinese hamster inbred strains.     

Joint estimation of contemporary seed and pollen dispersal rates among plant populations

There are few statistical methods for estimating contemporary dispersal among plant populations. A maximum-likelihood procedure is introduced here that uses pre- and post-dispersal population samples of biparentally inherited genetic markers to jointly estimate contemporary seed and pollen immigration rates from a set of discrete external sources into a target population. Monte Carlo simulations indicate that accurate estimates and reliable confidence intervals can be obtained using this method for both pollen and seed migration rates at modest sample sizes (100 parents/population and 100 offspring) when population differentiation is moderate (F(ST) ≥ 0.1), or by increasing pre-dispersal samples (to about 500 parents/population) when genetic divergence is weak (F(ST) = 0.01). The method exhibited low sensitivity to the number of source populations and achieved good accuracy at affordable genetic resolution (10 loci with 10 equifrequent alleles each). Unsampled source populations introduced positive biases in migration rate estimates from sampled sources, although they were minor when the proportion of immigration from the latter was comparatively low. A practical application of the method to a metapopulation of the Australian resprouter shrub Banksia attenuata revealed comparable levels of directional seed and pollen migration among dune groups, and the estimate of seed dispersal was higher than a previous estimate based on conservative assignment tests. The method should be of interest to researchers and managers assessing broad-scale nonequilibrium seed and pollen gene flow dynamics in plants.     

Molecular and quantitative genetic differentiation across Europe in yellow dung flies

Relating geographic variation in quantitative traits to underlying population structure is crucial for understanding processes driving population differentiation, isolation and ultimately speciation. Our study represents a comprehensive population genetic survey of the yellow dung fly Scathophaga stercoraria, an important model organism for evolutionary and ecological studies, over a broad geographic scale across Europe (10 populations from the Swiss Alps to Iceland). We simultaneously assessed differentiation in five quantitative traits (body size, development time, growth rate, proportion of diapausing individuals and duration of diapause), to compare differentiation in neutral marker loci (F(ST)) to that of quantitative traits (Q(ST)). Despite long distances and uninhabitable areas between sampled populations, population structuring was very low but significant (F(ST) = 0.007, 13 microsatellite markers; F(ST) = 0.012, three allozyme markers; F(ST) = 0.007, markers combined). However, only two populations (Iceland and Sweden) showed significant allelic differentiation to all other populations. We estimated high levels of gene flow [effective number of migrants (Nm) = 6.2], there was no isolation by distance, and no indication of past genetic bottlenecks (i.e. founder events) and associated loss of genetic diversity in any northern or island population. In contrast to the low population structure, quantitative traits were strongly genetically differentiated among populations, following latitudinal clines, suggesting that selection is responsible for life history differentiation in yellow dung flies across Europe.     

Population genetic structure of the endangered tropical tree species Caryocar brasiliense, based on variability at microsatellite loci

We report the population genetic structure of the endangered tropical tree species Caryocar brasiliense, based on variability at 10 microsatellite loci. Additionally, we compare heterozygosity and inbreeding estimates for continuous and fragmented populations and discuss the consequences for conservation. For a total of 314 individuals over 10 populations, the number of alleles per locus ranged from 20 to 27 and expected and observed heterozygosity varied from 0.129 to 0.924 and 0.067 to 1.000, respectively. Significant values of theta and R(ST) showed important genetic differentiation among populations. theta was much lower than R(ST), suggesting that identity by state and identity by descent have diverged in these populations. Although a significant amount of inbreeding was found under the identity by descent model (f = 0.11), an estimate of inbreeding for microsatellite markers based on a more adequate stepwise mutation model showed no evidence of nonrandom mating (R(IS) = 0.04). Differentiation (pairwise F(ST)) was positively correlated with geographical distance, as expected under the isolation by distance model. No effect of fragmentation on heterozygosity or inbreeding could be detected. This is most likely due to the fact that Cerrado fragmentation is a relatively recent event (approximately 60 years) compared to the species life cycle. Also, the populations surveyed from both fragmented and disturbed areas were composed mainly of adult individuals, already present prior to ecosystem fragmentation. Adequate hypothesis testing of the effect of habitat fragmentation will require the recurrent analysis of juveniles across generations in both fragmented and nonfragmented areas.     

Power for detecting genetic divergence: differences between statistical methods and marker loci

Information on statistical power is critical when planning investigations and evaluating empirical data, but actual power estimates are rarely presented in population genetic studies. We used computer simulations to assess and evaluate power when testing for genetic differentiation at multiple loci through combining test statistics or P values obtained by four different statistical approaches, viz. Pearson's chi-square, the log-likelihood ratio G-test, Fisher's exact test, and an F(ST)-based permutation test. Factors considered in the comparisons include the number of samples, their size, and the number and type of genetic marker loci. It is shown that power for detecting divergence may be substantial for frequently used sample sizes and sets of markers, also at quite low levels of differentiation. The choice of statistical method may be critical, though. For multi-allelic loci such as microsatellites, combining exact P values using Fisher's method is robust and generally provides a high resolving power. In contrast, for few-allele loci (e.g. allozymes and single nucleotide polymorphisms) and when making pairwise sample comparisons, this approach may yield a remarkably low power. In such situations chi-square typically represents a better alternative. The G-test without Williams's correction frequently tends to provide an unduly high proportion of false significances, and results from this test should be interpreted with great care. Our results are not confined to population genetic analyses but applicable to contingency testing in general.     

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

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

Genetic polymorphism of six DNA loci in six population groups of India

The genetic profile based on autosomal markers, four microsatellite DNA markers (D8S315, FES, D8S592, and D2S1328) and two minisatellite DNA markers (TPMT and PDGFA), were analyzed in six endogamous populations to examine the effect of geographic and linguistic affiliation on the genetic affinities among the groups. The six populations are from three different states of India and are linguistically different. Marathas from western India speak Marathi, an Indo-European language. Arayas, Muslims, Ezhavas, and Nairs from Kerala state of South India speak Malayalam, and Iyers from Tamil Nadu state speak Tamil. Genomic DNA was extracted from peripheral blood samples of random, normal, healthy individuals. Locus-specific PCR amplification was carried out, followed by electrophoresis of the amplicons and genotyping. All the loci were highly polymorphic and followed Hardy-Weinberg equilibrium, except for loci D8S315 and PDGFA in Iyers and Marathas, respectively. All six loci had high heterozygosity (average heterozygosity ranged from 0.73 to 0.76) and high polymorphism information content (0.57-0.90). The extent of gene differentiation among the six populations (G(ST) = 0.030) was greater than that for four Kerala populations (G(ST) = 0.011), suggesting proximity between the four Kerala populations. This result conforms with the cultural and linguistic background of the populations. The extent of diversity found among the populations probably resulted from the strict endogamous practices that they follow.     

Allelic and population variation of microsatellite loci in aspen (Populus tremuloides)

To develop a robust basis for inferences about population genetics and evolution, this work assayed 192 aspens (Populus tremuloides) from 11 sites in Wisconsin, USA, for allelic and population variation at 16 microsatellite loci distributed across the Populus genome. Frequency distributions of fluorochrome-labeled alleles resolved by capillary electrophoresis were analyzed for relationships to repeat size and number. Population-level statistics were compared with those of other studies, especially in Populus. All loci were polymorphic, varying widely in the number of alleles per locus (mean = 8.25, range 2-20). Expected and observed heterozygosities were high (0.45 and 0.41, respectively), with little differentiation among populations (F(ST) = 0.006-0.045) and a moderate level of inbreeding (F(IS) = 0.09), intermediate among levels reported in studies based on isozymes. Contrary to several other reports, allele frequencies clustered tightly around the modal frequency, and the genetic diversity (measured as alleles per locus or as expected heterozygosity) was not related to either the repeat unit size or to the number of repeats.     

Variation for neutral markers is correlated with variation for quantitative traits in the plant pathogenic fungus Mycosphaerella graminicola

We compared genetic variation and population differentiation at RFLP marker loci with seven quantitative characters including fungicide resistance, temperature sensitivity, pycnidial size, pycnidial density, colony size, percentage of leaves covered by pycnidia (PLACP) and percentage of leaves covered by lesions (PLACL) in Mycosphaerella graminicola populations sampled from four regions. Wide variation in population differentiation was found across the quantitative traits assayed. Fungicide resistance, temperature sensitivity, and PLACP displayed a significantly higher Q(ST) than G(ST), consistent with selection for local adaptation, while pycnidial size, pycnidial density and colony size displayed a lower or significantly lower Q(ST) than G(ST), consistent with constraining selection. There was not a statistical difference between Q(ST) and G(ST) in PLACL. We also found a positive and significant correlation between genetic variation in molecular marker loci and quantitative traits at the multitrait scale, suggesting that estimates of overall genetic variation for quantitative traits in M. graminicola could be derived from analysis of the molecular genetic markers.     

Microsatellite variation in cassava (Manihot esculenta, Euphorbiaceae) and its wild relatives: further evidence for a southern Amazonian origin of domestication

Genetic variation at five microsatellite loci was used to investigate the evolutionary and geographical origins of cassava (Manihot esculenta subsp. esculenta) and the population structure of cassava's wild relatives. Two hundred and twelve individuals were sampled, representing 20 crop accessions, 27 populations of cassava's closest wild relative (M. esculenta subsp. flabellifolia), and six populations of a potentially hybridizing species (M. pruinosa). Seventy-three alleles were observed across all loci and populations. These data indicate the following on cassava's origin: (1) genetic variation in the crop is a subset of that found in the wild M. esculenta subspecies, suggesting that cassava is derived solely from its conspecific wild relative. (2) Phenetic analyses group cassava with wild populations from the southern border of the Amazon basin, indicating this region as the likely site of domestication. (3) Manihot pruinosa, while closely related to M. esculenta (and possibly hybridizing with it where sympatric), is probably not a progenitor of the crop. Genetic differentiation among the wild populations is moderately high (F:(ST) = 0.42, rho(ST) = 0.54). This differentiation has probably arisen primarily through random genetic drift (rather than mutation) following recent population divergence.     

Genetic variation and structure of house sparrow populations: is there an island effect?

Population genetic structure and intrapopulation levels of genetic variation have important implications for population dynamics and evolutionary processes. Habitat fragmentation is one of the major threats to biodiversity. It leads to smaller population sizes and reduced gene flow between populations and will thus also affect genetic structure. We use a natural system of island and mainland populations of house sparrows along the coast of Norway to characterize the different population genetic properties of fragmented populations. We genotyped 636 individuals distributed across 14 populations at 15 microsatellite loci. The level of genetic differentiation was estimated using F‐statistics and specially designed Mantel tests were conducted to study the influence of population type (i.e. mainland or island) and geographic distance on the genetic population structure. Furthermore, the effects of population type, population size and latitude on the level of genetic variation within populations were examined. Our results suggest that genetic processes on islands and mainland differed in two important ways. First, the intrapopulation level of genetic variation tended to be lower and the occurrence of population bottlenecks more frequent on islands than the mainland. Second, although the general level of genetic differentiation was low to moderate, it was higher between island populations than between mainland populations. However, differentiation increased in mainland populations somewhat faster with geographical distance. These results suggest that population bottleneck events and genetic drift have been more important in shaping the genetic composition of island populations compared with populations on the mainland. Such knowledge is relevant for a better understanding of evolutionary processes and conservation of threatened populations.     

Genetic diversity and differentiation of dromedarian camel of India

Estimation of genetic variability and relationship among different livestock breeds is important for management of genetic resources for their sustainable utilization and conservation. This is more important when the livestock species, like camel, have shown a sharp decline in head count during the last decade. In the present study we estimated genetic variability and relationship among four camel breeds of India using 23 microsatellite loci. A total of 252 alleles were observed across all the four populations with mean number of alleles per locus as 8.04, 7.30, 6.39, and 7.43 for Bikaneri, Jaisalmeri, Kutchi, and Mewari breeds, respectively. The mean observed heterozygosity of the four breeds were 0.58, 0.57, 0.56, and 0.60 for Bikaneri, Jaisalmeri, Kutchi, and Mewari breeds, respectively and were lower than expected heterozygosity values. The mean estimates of F statistics were 0.227+/-0.044 (F(IT)), 0.157+/-0.038 (F(IS)), and 0.082+/-0.019 (F(ST)). The values were significantly different from zero for all the three measures and point towards the existence of population structure and moderate differentiation in four camel breeds. The exact test also indicated significant population differentiation (P < 0.001). The analysis of molecular variance revealed 12% of the variation attributed to among populations and 88% within populations. Sixty-nine percent of the individuals could be correctly assigned using "leave one out" procedure. All the individuals of Mewari and 42 out of 44 Jaisalmeri were correctly assigned. The existence of strong population structure in Jaisalmeri and Mewari camel was further substantiated by Nei's standard genetic distance as well as interindividual allele sharing distance. Thus these two breeds owing to selection for specific traits are distinct from other camel breeds.     

Genetic status of Asiatic black bear (Ursus thibetanus) reintroduced into South Korea based on mitochondrial DNA and microsatellite loci analysis

The Asiatic black bear is one of the most endangered mammals in South Korea owing to population declines resulting from human exploitation and habitat fragmentation. To restore the black bear population in South Korea, 27 bear cubs from North Korea and Russian Far East (Primorsky Krai) were imported and released into Jirisan National Park, a reservoir of the largest wild population in South Korea, in 2004. To monitor the success of this reintroduction, the genetic diversity and population structure of the reintroduced black bears were measured using both mitochondrial and nuclear DNA markers. Mitochondrial D-loop region DNA sequences (615 bp) of 43 Japanese black bears from previous study and 14 Southeast Asian black bears in this study were employed to obtain phylogenetic inference of the reintroduced black bears. The mitochondrial phylogeny indicated Asiatic black bear populations from Russian Far East and North Korea form a single evolutionary unit distinct from populations from Japan and Southeast Asia. Mean expected heterozygosity (H(E)) across 16 microsatellite loci was 0.648 for Russian and 0.676 for North Korean populations. There was a moderate but significant level of microsatellite differentiation (F(ST) = 0.063) between black bears from the 2 source areas. In addition, genetic evidences revealed that 2 populations are represented as diverging groups, with lingering genetic admixture among individuals of 2 source populations. Relatedness analysis based on genetic markers indicated several discrepancies with the pedigree records. Implication of the phylogenetic and genetic evidences on long-term management of Asiatic black bears in South Korea is discussed.