Conserving marine biodiversity: insights from life-history trait candidate genes in Atlantic cod (Gadus morhua)

Recent technological developments have facilitated an increased focus on identifying genomic regions underlying adaptive trait variation in natural populations, and it has been advocated that this information should be important for designating population units for conservation. In marine fishes, phenotypic studies have suggested adaptation through divergence of life-history traits among natural populations, but the distribution of adaptive genetic variation in these species is still relatively poorly known. In this study, we extract information about the geographical distribution of genetic variation for 33 single nucleotide polymorphisms (SNPs) associated with life-history trait candidate genes, and compare this to variation in 70 putatively neutral SNPs in Atlantic cod (Gadus morhua). We analyse samples covering the major population complexes in the eastern Atlantic and find strong evidence for non-neutral levels and patterns of population structuring for several of the candidate gene-associated markers, including two SNPs in the growth hormone 1 gene. Thus, this study aligns with findings from phenotypic studies, providing molecular data strongly suggesting that these or closely linked genes are under selection in natural populations of Atlantic cod. Furthermore, we find that patterns of variation in outlier markers do not align with those observed at selectively neutral markers, and that outlier markers identify conservation units on finer geographical scales than those revealed when analysing only neutral markers. Accordingly, results also suggest that information about adaptive genetic variation will be useful for targeted conservation and management in this and other marine species.     

Imprints from genetic drift and mutation imply relative divergence times across marine transition zones in a pan-European small pelagic fish (Sprattus sprattus)

Geographic distributions of most temperate marine fishes are affected by postglacial recolonisation events, which have left complex genetic imprints on populations of marine species. This study investigated population structure and demographic history of European sprat (Sprattus sprattus L.) by combining inference from both mtDNA and microsatellite genetic markers throughout the species' distribution. We compared effects from genetic drift and mutation for both genetic markers in shaping genetic differentiation across four transition zones. Microsatellite markers revealed significant isolation by distance and a complex population structure across the species' distribution (overall θ(ST)=0.038, P<0.01). Across transition zones markers indicated larger effects of genetic drift over mutations in the northern distribution of sprat contrasting a stronger relative impact of mutation in the species' southern distribution in the Mediterranean region. These results were interpreted to reflect more recent divergence times between northern populations in accordance with previous findings. This study demonstrates the usefulness of comparing inference from different markers and estimators of divergence for phylogeographic and population genetic studies in species with weak genetic structure, as is the case in many marine species.     

A comparison of biallelic markers and microsatellites for the estimation of population and conservation genetic parameters in Atlantic salmon (Salmo salar)

Biallelic markers such as single nucleotide polymorphisms (SNPs) and insertion/deletion polymorphisms have become increasingly popular markers for various population genetics applications. However, the effort required to develop biallelic markers in nonmodel organisms is still substantial. In this study, we compared the estimation of various population genetic parameters (genetic divergence and structuring, isolation-by-distance, genetic diversity) using a limited number of biallelic markers (in total 7 loci) to those estimated with 14 microsatellite loci in 21 Atlantic salmon (Salmo salar) populations from northern Europe. Pairwise FST values were significantly correlated between biallelic loci and microsatellite datasets, as was overall heterozygosity when both anadromous and nonanadromous populations were analyzed together. However, when the anadromous and nonanadromous samples were analyzed separately, only genetic divergence correlations remained significant. Biallelic markers alone were not sufficient for reliable neighbor-joining clustering of populations but gave highly similar isolation-by-distance signals when compared with microsatellites. Finally, although several population prioritization measures for conservation exhibited significant correlation between different marker types, the specific populations highlighted as being most valuable for conservation purposes varied depending on the marker type and conservation criteria applied. This study demonstrates that a relatively small set of biallelic markers can be sufficient for obtaining concordant results in most of the analyses compared with microsatellites, although estimates of genetic distance are generally more concordant than estimates of genetic diversity. This suggests that a relatively small number of biallelic markers can provide useful information for various population genetic applications. However, we emphasize that the use of much higher number of loci is preferable, especially when the genetic differences between populations are subtle or individual multilocus genotype-based analyses are to be performed.     

Genetic diversity in European chestnut populations by means of genomic and genic microsatellite markers

Microsatellite or simple sequence repeats (SSRs) are one of the most used markers in population genetic studies. SSR markers developed from expressed sequence tags (EST) have proved useful to examine functional diversity in relation to adaptive variation. The information provided by both genomic and genic microsatellite markers could offer more accurate indication on the distribution of the genetic diversity among and within populations assuming different evolutionary drivers. This is the first study on chestnut (Castanea sativa Mill.) in which the genetic diversity was evaluated by means of genomic (SSRs) and genic (EST-SSRs) microsatellite markers. We genotyped nine natural European chestnut populations distributed throughout representative areas of contrasting climatic conditions in the Mediterranean basin. Genomic SSRs showed significantly higher levels of diversity in terms of number of alleles, effective number of alleles, expected heterozygosity and level of polymorphism. Furthermore, there were significant differences in the level of differentiation among populations. The UPGMA analysis revealed different clustering pattern between populations, being the grouping according to geographic distances in the case of genomic SSRs and two differentiated groups based on the northern–southern distribution of the populations for EST-SSRs. Furthermore, the EST-SSR transferability among related Castanea and Quercus species was stated. Our results confirm that combining genomic SSRs and EST-SSRs is a useful tool to give complementary information to explain the genetic and adaptive diversity in chestnut.     

Genetic characteristics of 207 microsatellite markers in the Korean population and in other Asian populations

Microsatellites, short tandem repeats, are useful markers for genetic analysis because of their high frequency of occurrence over the genome, high information content due to variable repeat lengths, and ease of typing. To establish a panel of microsatellite markers useful for genetic studies of the Korean population, the allele frequencies and heterozygosities of 207 microsatellite markers in 119 unrelated Korean, Indian and Pakistani individuals were compared. The average heterozygosity of the Korean population was 0.71, similar to that of the Indian and Pakistani populations. More than 80% of the markers showed heterozygosity of over 0.6 and were valuable as genetic markers for genome-wide screening for disease susceptibility loci in these populations. To identify the allelic distributions of the multilocus genetic data from these microsatellite markers, the population structures were assessed by clustering. These markers supported, with the most probability, three clustering groups corresponding to the three geographical populations. When we assumed only two hypothetical clusters (K), the Korean population was separate from the others, suggesting a relatively deep divergence of the Korean population. The present 207 microsatellite markers appear to reflect the historical and geographical origins of the different populations as well as displaying a similar degree of variation to that seen in previously published genetic data. Thus, these markers will be useful as a reference for human genetic studies on Asians.     

X-Chromosome STR markers data in a Cabo Verde immigrant population of Lisboa

Population genetic data of 12 X chromosomal short tandem repeats markers (DXS10074, DXS10079, DXS10101, DXS10103, DXS10134, DXS10135, DXS10146, DXS10148, DXS7132, DXS7423, DXS8378 and HPRTB) were analysed in 54 females and 95 males of an immigrant population from Cabo Verde living in Lisboa. The obtained results for forensic statistical parameters such as observed heterozigosity, polymorphism information content, power of discrimination and mean exclusion chance, based on single allele frequencies, reveal that this multiplex system is highly informative and can represent an important tool for genetic identification purposes in the immigrant population of Cabo Verde. Since the studied short tandem repeats genetic markers are distributed on four linkage groups, that can provide independent genotype information, we studied those groups as haploytes. The forensic efficiency parameters for the linked groups were all higher than 0.97, with linkage group I being the most polymorphic and linkage group III the less informative.     

The mode of evolution of molecular markers in populations of house mice under artificial selection for locomotor behavior

A complete understanding of the mode of evolution of molecular markers is important for making inferences about different population genetic parameters, especially because a number of studies have reported patterns of allelic variation at molecular markers that are not in agreement with neutral evolutionary expectations. In the present study, house mice (Mus domesticus) from the fourteenth generation of a selection experiment for increased voluntary wheel-running activity were used to test how selection on a complex behavior affects the distribution of allelic variation by examining patterns of variation at six microsatellite and four allozyme loci. This population had a hierarchical structure that allowed for simultaneous testing of the effects of selection and genetic drift on the distribution of allelic variation by comparing observed patterns of allele frequencies and estimates of genetic divergence at multiple hierarchical levels to expectations under models of neutral evolution. The levels of genetic divergence among replicate lines and between selection groups, estimated from microsatellite data or pooled microsatellite and allozyme data, were not significantly different from expectations under neutral evolution. Furthermore, the pattern of change of allele frequencies between the base population and generation 14 was largely in agreement with expectations under neutral evolution (although the PGM locus exhibited a pattern of change within populations that was difficult to explain under neutral evolution). Overall the results generally provide support for the neutral evolution of molecular markers.     

Inter‐oceanic variation in patterns of host‐associated divergence in a seabird ectoparasite

Aim Parasites with global distributions and wide host spectra provide excellent models for exploring the factors that drive parasite diversification. Here, we tested the relative force of host and geography in shaping population structure of a widely distributed and common ectoparasite of colonial seabirds, the tick Ixodes uriae. Location Two natural geographic replicates of the system: numerous seabird colonies of the North Pacific and North Atlantic Ocean basins. Methods Using eight microsatellite markers and tick samples from a suite of multi‐specific seabird colonies, we examined tick population structure in the North Pacific and compare patterns of diversity and structure to those in the Atlantic basin. Analyses included population genetic estimations of diversity and population differentiation, exploratory multivariate analyses, and Bayesian clustering approaches. These different analyses explicitly took into account both the geographic distance among colonies and host use by the tick. Results Overall, little geographic structure was observed among Pacific tick populations. However, host‐related genetic differentiation was evident, but was variable among host types and lower than in the North Atlantic. Main conclusions Tick population structure is concordant with the genetic structure observed in seabird host species within each ocean basin, where seabird populations tend to be less structured in the North Pacific than in the North Atlantic. Reduced tick genetic structure in the North Pacific suggests that host movement among colonies, and thus tick dispersal, is higher in this region. In addition to information on parasite diversity and gene flow, our findings raise interesting questions about the subtle ways that host behaviour, distribution and phylogeographic history shape the genetics of associated parasites across geographic landscapes.     

Estimation of multilocus linkage disequilibria in diploid populations with dominant markers

Analysis of population structure and organization with DNA-based markers can provide important information regarding the history and evolution of a species. Linkage disequilibrium (LD) analysis based on allelic associations between different loci is emerging as a viable tool to unravel the genetic basis of population differentiation. In this article, we derive the EM algorithm to obtain the maximum-likelihood estimates of the linkage disequilibria between dominant markers, to study the patterns of genetic diversity for a diploid species. The algorithm was expanded to estimate and test linkage disequilibria of different orders among three dominant markers and can be technically extended to manipulate an arbitrary number of dominant markers. The feasibility of the proposed algorithm is validated by an example of population genetic studies of hickory trees, native to southeastern China, using dominant random amplified polymorphic DNA markers. Extensive simulation studies were performed to investigate the statistical properties of this algorithm. The precision of the estimates of linkage disequilibrium between dominant markers was compared with that between codominant markers. Results from simulation studies suggest that three-locus LD analysis displays increased power of LD detection relative to two-locus LD analysis. This algorithm is useful for studying the pattern and amount of genetic variation within and among populations.     

Genetic diversity and bottleneck analysis of Indian bellary sheep by microsatellite markers

Bellary sheep population variability and structure was investigated genetically utilizing FAO recommended microsatellite markers. Genetic variation at 20 microsatellite loci, population structure, and genetic bottleneck hypothesis were examined. Estimates of genetic variability such as effective number of alleles and gene diversities revealed substantial genetic variation frequently displayed by microsatellite markers. A total of 133 alleles were detected. Average polymorphism across the studied loci and expected gene diversity in the population were 1.419 +/- 0.405 and 0.684 +/- 0.140, respectively. No significant genotypic linkage disequilibrium was detected across population, suggesting no evidence of linkage between loci. The population was observed to be significantly differentiated into different groups, showed fairly high level of inbreeding (f = 0.253 +/- 0.050) and global heterozygote deficit. Population structure analysis indicated the intermixing/introduction of unique/rare alleles in these migrating flocks. A normal L-shaped distribution of mode-shift test, non-significant heterozygosity excess on the basis of different models, as revealed from Sign, Standardized differences and Wilcoxon sign rank tests suggested that there was no recent bottleneck. The study revealed that even breed with increasing population trend needs genetic management for the conservation and improvement.     

A multi-marker test based on family data in genome-wide association study

Background

Requirements for successful implementation of multivariate animal threshold models including phenotypic and genotypic information are not known yet. Here simulated horse data were used to investigate the properties of multivariate estimators of genetic parameters for categorical, continuous and molecular genetic data in the context of important radiological health traits using mixed linear-threshold animal models via Gibbs sampling. The simulated pedigree comprised 7 generations and 40000 animals per generation. Additive genetic values, residuals and fixed effects for one continuous trait and liabilities of four binary traits were simulated, resembling situations encountered in the Warmblood horse. Quantitative trait locus (QTL) effects and genetic marker information were simulated for one of the liabilities. Different scenarios with respect to recombination rate between genetic markers and QTL and polymorphism information content of genetic markers were studied. For each scenario ten replicates were sampled from the simulated population, and within each replicate six different datasets differing in number and distribution of animals with trait records and availability of genetic marker information were generated. (Co)Variance components were estimated using a Bayesian mixed linear-threshold animal model via Gibbs sampling. Residual variances were fixed to zero and a proper prior was used for the genetic covariance matrix.

Results

Effective sample sizes (ESS) and biases of genetic parameters differed significantly between datasets. Bias of heritability estimates was -6% to +6% for the continuous trait, -6% to +10% for the binary traits of moderate heritability, and -21% to +25% for the binary traits of low heritability. Additive genetic correlations were mostly underestimated between the continuous trait and binary traits of low heritability, under- or overestimated between the continuous trait and binary traits of moderate heritability, and overestimated between two binary traits. Use of trait information on two subsequent generations of animals increased ESS and reduced bias of parameter estimates more than mere increase of the number of informative animals from one generation. Consideration of genotype information as a fixed effect in the model resulted in overestimation of polygenic heritability of the QTL trait, but increased accuracy of estimated additive genetic correlations of the QTL trait.

Conclusion

Combined use of phenotype and genotype information on parents and offspring will help to identify agonistic and antagonistic genetic correlations between traits of interests, facilitating design of effective multiple trait selection schemes.     

Population structure, habitat features and genetic structure of managed red deer populations

Management of game ungulates alters population structure and habitat features, with potential effects on genetic structure. Here, we study 26 red deer (Cervus elaphus) populations in Spain. We used census data and habitat features as well as genetic information at 11 microsatellite markers from 717 individuals. We found that metapopulations presented a distribution associated with forest interruptions. Within metapopulations, fences did not have a significant effect on red deer genetic structure. The metapopulations we studied presented similar population structure, but they differed in habitat features and genetic structure. The metapopulation with higher resource availability showed a genetic structure pattern in which genetic relatedness between geographically close individuals was high while relatedness between geographically distant individuals was low. Contrarily, the metapopulation with lower resource availability presented a genetic structure pattern in which the genetic relatedness between individuals of different populations was independent of the geographic distance. We discuss the possible connection between resource availability and genetic structure. Finally, we did not find any population or environmental variable related to genetic differentiation within metapopulations.     

Genetic diversity and phylogenetic analysis of blackbuck (Antilope cervicapra) in southern India

Blackbuck (Antilope cervicapra) is a threatened species endemic to the Indian subcontinent. Many populations of blackbuck are found in southern India. Populations of blackbuck are negatively affected in many places for various reasons, such as habitat destruction and poaching. Their range decreased sharply during the 20th century. There is very limited information available on the population dynamics of blackbuck in southern India. For the phylogenetic and genetic diversity analyses of blackbuck populations among different distribution ranges in southern India, we sequenced mt DNA of cytochrome b (Cyt b) for 120, cytochrome c oxidase subunit-1 (COI) for 137 and the control region (CR) for 137 fecal pellets from eleven different locations in southern India. We analyzed the genetic structure of three mitochondrial markers, the CR, Cyt b and the COI region, separately and in a combined dataset. The haplotype diversity and nucleotide diversity of CR were 0.969 and 0.047, respectively, and were higher than those of Cyt b and COI. A Bayesian phylogeny and an MJ network based on the CR and combined dataset (105 sequences) signified several distinct haplotype clusters within blackbuck, whereas no clusters were identified with the Cyt b and COI phylogenetic analyses. The analysis of molecular variance of the combined data set revealed 52.46% genetic variation within the population. Mismatch distribution analysis revealed that blackbuck populations underwent complex changes with analysis of the combined dataset in each population and analysis of each marker separately in the overall population. The results provide evidence that blackbuck in different geographic locations has a distinct population structure due to habitat fragmentation after the formation of the Western and Eastern Ghats.

     

省沽油群体遗传多样性的AFLP分析

利用4对AFLP引物对来自安徽石台、湖北大悟和河南桐柏的3个省沽油群体进行遗传多样性分析,共扩增出489条带,其中多态性条带460条,多态性百分比为93.99%。不同群体Nei’s基因多样性指数(H)和Shannon’s信息指数(I)变化范围分别为0.1920～0.2046和0.2937～0.3151,其中湖北大悟群体遗传多样性最高。物种水平和种群水平的H分别为0.2190和0.1964,群体内变异占总变异的89.68%,表明省沽油遗传变异主要存在于各群体内部。3个群体的平均遗传距离为0.0292,UPGMA聚类分析结果说明省沽油各群体间亲缘关系较近并和地域具有相关性。建议省沽油的遗传资源保护应以种内遗传多样性的保护为主。     

Limited Reintroduction Does Not Always Lead to Rapid Loss of Genetic Diversity: An Example from the American Chestnut (Castanea dentata; Fagaceae)

In restoring species, reasons for introducing limited numbers of individuals at different locations include costs of introduction and maintenance, limited founder supply, and risk “bet hedging.” However, populations initiated from few founders may experience increased genetic drift, inbreeding, and diversity loss. We examined the genetic diversity of an isolated stand of more than 5,000 American chestnut trees relative to that of the 9 surviving stand founders (out of 10 total) planted in the 1880s. We used minisatellite DNA probes to reveal 84 genetic markers (circa 24 loci) among the nine founders, and their genetic diversity was compared with three separate plots of descendant trees, as well as with two natural stands. The descendants were circa 7.3% more heterozygous than the founders (mean estimated H= 0.556 vs. 0.518, respectively; p < 0.0001). Genetic differentiation was not pronounced (F_ST < 0.031), and no markers, including those at low frequency among the founders, were lost in the descendants. The founders and natural transects were not significantly different in H or similarity (mean proportion of bands shared). Special planting or mating protocols for establishment of a vigorous American chestnut population from a low number of founders may not be required to avoid strong effects of genetic drift and inbreeding. These results demonstrate that loss of genetic diversity following reintroduction of a limited number of founders is not always inevitable, such as this case where the species is highly outcrossing, expression of heterozygous advantage may occur, the original founders remain as gene contributors over generations, and the establishing population expands constantly and rapidly.     

Geographical variation of genetic and phenotypic traits in the Mexican sailfin mollies, Poecilia velifera and P. petenensis

Comparing the patterns of population divergence using both neutral genetic and phenotypic traits provides an opportunity to examine the relative importance of evolutionary mechanisms in shaping population differences. We used microsatellite markers to examine population genetic structure in the Mexican sailfin mollies Poecilia velifera and P. petenensis. We compared patterns of genetic structure and divergence to that in two types of phenotypic traits: morphological characters and mating behaviours. Populations within each species were genetically distinct, and conformed to a model of isolation by distance, with populations within different geographical regions being more genetically similar to one another than were populations from different regions. Bayesian clustering and barrier analyses provided additional support for population separation, especially between geographical regions. In contrast, none of the phenotypic traits showed any type of geographical pattern, and population divergence in these traits was uncorrelated with that found in neutral markers. There was also a weaker pattern of regional differences among geographical regions compared to neutral genetic divergence. These results suggest that while divergence in neutral traits is likely a product of population history and genetic drift, phenotypic divergence is governed by different mechanisms, such as natural and sexual selection, and arises at spatial scales independent from those of neutral markers.     

Social and genetic population structure of free-ranging cheetah in Botswana: implications for conservation

Once widely distributed throughout Africa, cheetahs (Acinonyx jubatus) occur today within fragmented populations and are listed as vulnerable by the IUCN. Botswana currently hosts the second largest cheetah population throughout the species’ range. This study initiated a molecular genetic survey of wild Botswana cheetah populations. It focused on the relatedness within presumed social groups using 14 microsatellite markers and revealed a higher proportion of unrelated male coalitions than was expected. Based on the unrelated cheetahs only, the estimation of the genetic variation corresponded with results from recent studies on different African populations. The analysis of unrelated individuals indicated limited genetic differentiation between cheetahs from different regions of Botswana. This suggests that the Botswana cheetah population might represent a unique panmictic population as long as sufficient levels of gene flow are maintained within the distribution range. This baseline information will now be incorporated to develop management strategies and set priorities for cheetah conservation in Botswana.     

Bayesian prediction of breeding values for multivariate binary and continuous traits in simulated horse populations using threshold-linear models with Gibbs sampling

Simulated data were used to determine the properties of multivariate prediction of breeding values for categorical and continuous traits using phenotypic, molecular genetic and pedigree information by mixed linear-threshold animal models via Gibbs sampling. Simulation parameters were chosen such that the data resembled situations encountered in Warmblood horse populations. Genetic evaluation was performed in the context of the radiographic findings in the equine limbs. The simulated pedigree comprised seven generations and 40 000 animals per generation. The simulated data included additive genetic values, residuals and fixed effects for one continuous trait and liabilities of four binary traits. For one of the binary traits, quantitative trait locus (QTL) effects and genetic markers were simulated, with three different scenarios with respect to recombination rate (r) between genetic markers and QTL and polymorphism information content (PIC) of genetic markers being studied: r = 0.00 and PIC = 0.90 (r0p9), r = 0.01 and PIC = 0.90 (r1p9), and r = 0.00 and PIC = 0.70 (r0p7). For each scenario, 10 replicates were sampled from the simulated horse population, and six different data sets were generated per replicate. Data sets differed in number and distribution of animals with trait records and the availability of genetic marker information. Breeding values were predicted via Gibbs sampling using a Bayesian mixed linear-threshold animal model with residual covariances fixed to zero and a proper prior for the genetic covariance matrix. Relative breeding values were used to investigate expected response to multi- and single-trait selection. In the sires with 10 or more offspring with trait information, correlations between true and predicted breeding values ranged between 0.89 and 0.94 for the continuous traits and between 0.39 and 0.77 for the binary traits. Proportions of successful identification of sires of average, favourable and unfavourable genetic value were 81% to 86% for the continuous trait and 57% to 74% for the binary traits in these sires. Expected decrease of prevalence of the QTL trait was 3% to 12% after multi-trait selection for all binary traits and 9% to 17% after single-trait selection for the QTL trait. The combined use of phenotype and genotype data was superior to the use of phenotype data alone. It was concluded that information on phenotypes and highly informative genetic markers should be used for prediction of breeding values in mixed linear-threshold animal models via Gibbs sampling to achieve maximum reduction in prevalences of binary traits.     

Genetic structure of Araucaria angustifolia (Araucariaceae) populations in Brazil: implications for the in situ conservation of genetic resources

The distribution of the genetic variation within and among natural populations of A. ANGUSTIFOLIA growing in different regions in Brazil was assessed at microsatellite and AFLP markers. Both markers revealed high gene diversity ( H = 0.65; AR = 9.1 for microsatellites and H = 0.27; P = 77.8 % for AFLPs), moderate overall differentiation ( RST = 0.13 for microsatellites and FST = 0.10 for AFLPs), but high divergence of the northernmost, geographically isolated population. In a Bayesian analysis, microsatellite data suggested population structure at two levels: at K = 2 and at K = 3 in agreement to the geographical distribution of populations. This result was confirmed by the UPGMA dendrogram based on microsatellite data (bootstrap support > 95 %). Non-hierarchical AMOVA revealed high variation among populations from different A POSTERIORI defined geographical groups. The genetic distance between sample locations increased with geographical distance for microsatellites ( R = 0.62; P = 0.003) and AFLPs ( R = 0.32; P = 0.09). This pattern of population differentiation may be correlated with population history such as geographical isolation and postglacial colonization of highlands. Implications of the population genetic structure for the conservation of genetic resources are discussed.     

High genetic differentiation among remnant populations of the endangered Caesalpinia echinata Lam. (Leguminosae–Caesalpinioideae)

Forest fragments along the Atlantic coastland of Brazil have been highly impacted by extensive human activities for the last 400 years. Caesalpinia echinata (Leguminosae– Caesalpinioideae), brazilwood, was overexploited during this period due to its economical importance as a dye. As a result, the species has become endangered and today its total population size is very restricted. We have assessed the distribution of genetic variation between five natural populations of brazilwood by means of RAPD (random amplified polymorphic DNA) markers. Of the total genetic variability, 28.5% was attributable to differences between two geographical groups, 29.6% to population differences within groups and 42.0% to individual differences within populations. The high level of population differentiation observed is in contrast to that expected for a primarily outcrossed woody perennial plant, and suggests that there may be a degree of inbreeding. Our results are in agreement with previous studies which postulated that C. echinata has always occurred in clumps, being common in some places but rare in between. From a conservation point of view, different populations representing different regions should be protected and, yet, plants with different origins should not be synthesized into populations in a recovery process at the risk of loss and dilution of genetic information. This study demonstrates that RAPD markers were effective in establishing a clear correlation between genetic and geographical distance and in identifying areas of maximum diversity, and may be used as an initial approach to assess the partitioning of genetic variation in this endangered species.