Relative mutation rates at di-, tri-, and tetranucleotide microsatellite loci

Data mining for simple sequence repeats in expressed sequence tags from barley, maize, rice, sorghum and wheat

Plant genomics projects involving model species and many agriculturally important crops are resulting in a rapidly increasing database of genomic and expressed DNA sequences. The publicly available collection of expressed sequence tags (ESTs) from several grass species can be used in the analysis of both structural and functional relationships in these genomes. We analyzed over 260000 EST sequences from five different cereals for their potential use in developing simple sequence repeat (SSR) markers. The frequency of SSR-containing ESTs (SSR-ESTs) in this collection varied from 1.5% for maize to 4.7% for rice. In addition, we identified several ESTs that are related to the SSR-ESTs by BLAST analysis. The SSR-ESTs and the related sequences were clustered within each species in order to reduce the redundancy and to produce a longer consensus sequence. The consensus and singleton sequences from each species were pooled and clustered to identify cross-species matches. Overall a reduction in the redundancy by 85% was observed when the resulting consensus and singleton sequences (3569) were compared to the total number of SSR-EST and related sequences analyzed (24 606). This information can be useful for the development of SSR markers that can amplify across the grass genera for comparative mapping and genetics. Functional analysis may reveal their role in plant metabolism and gene evolution.

Genomic microsatellites as evolutionary chronometers: a test in wild cats

Nuclear microsatellite loci (2- to 5-bp tandem repeats) would seem to be ideal markers for population genetic monitoring because of their abundant polymorphism, wide dispersal in vertebrate genomes, near selective neutrality, and ease of assessment; however, questions about their mode of generation, mutation rates and ascertainment bias have limited interpretation considerably. We have assessed the patterns of genomic diversity for ninety feline microsatellite loci among previously characterized populations of cheetahs, lions and pumas in recapitulating demographic history. The results imply that the microsatellite diversity measures (heterozygosity, allele reconstitution and microsatellite allele variance) offer proportionate indicators, albeit with large variance, of historic population bottlenecks and founder effects. The observed rate of reconstruction of new alleles plus the growth in the breadth of microsatellite allele size (variance) was used here to develop genomic estimates of time intervals following historic founder events in cheetahs (12,000 yr ago), in North American pumas (10,000-17,000 yr ago), and in Asiatic lions of the Gir Forest (1000-4000 yr ago).

Interchromosomal telomere length variation

Despite the recent discovery of interchromosomal telomere length variation, a role for heterogeneity in telomere maintenance has yet to be established. This study aimed to clarify relative telomere length differences between chromosomes. Combined chromosomal sorting and telomeric repeat content analysis in GM130B cells, the relative telomeric repeat content in each chromosome, were calculated. Each chromosome could be isolated except for chromosomes 1 and 2 and chromosomes 9 to 12, which were isolated in a group. Telomere length was correlated with the size of the corresponding chromosome. Concomitant relative telomeric repeat content analysis in each chromosome and terminal restriction fragment analysis using the whole genome revealed that the terminal restriction fragments of each chromosome were heterogenously distributed through the smear of the fragments. This is the first description of an association between telomere length and chromosome size.

Evidence of common ancestry for the maple syrup urine disease (MSUD) Y438N allele in non-Mennonite MSUD patients

Maple syrup urine disease (MSUD) is a rare (1/185,000) autosomal recessive inborn error of branched-chain amino acid metabolism characterized by increased plasma leucine, isoleucine, and valine levels. Though, genetically heterogeneous in the worldwide population, MSUD in Old Order Mennonites (1/150-176) is the result of a tyrosine to asparagine substitution (Y438N; previously Y393N) in the E1alpha subunit of the branched-chain alpha-keto acid dehydrogenase (BCKAD) complex. Due to endogamous practices, the presence of Y438N in all reported Mennonite MSUD patients has historically been attributed to a founder effect. However, we have also identified the Y438N defect in eight MSUD patients of non-Mennonite lineage. To evaluate the genetic origin of this defect in these non-Mennonite patients, we examined Mennonite MSUD families and non-Mennonite MSUD families using microsatellite markers located on chromosome 19q13.1-13.2 (location of E1alpha gene, BCKDHA). Haplotype analyses revealed a major and four minor haplotypes that cosegregate with the Y438N allele in the Old Order Mennonite MSUD patients and carrier relatives. Analyses of eight non-Mennonite MSUD patients reveal that three of the non-Mennonite MSUD patients shared common Mennonite Y438N haplotypes, strongly suggesting Mennonite ancestry. However, the remaining non-Mennonite patients carry Y438N haplotypes that are significantly different from the Mennonite Y438N haplotype, suggesting that the occurrence of the defect in these families is due to either pre-Mennonite or de novo events.

Non-fatherhood or mutation? A probabilistic approach to parental exclusion in paternity testing

The occurrence of germline mutations at microsatellite loci poses problems in ascertaining non-fatherhood status in paternity testing. We describe the appropriate probabilistic analysis for computing the likelihood ratio in favour of paternity while allowing for mutation, for all 18 relevant combinations of seemingly incompatible parental genotypes. We allow arbitrary and possibly different mutation rates in paternal and maternal germlines. We describe a stationary mutation model for expressing the required allele-specific transition mutation rates in terms of overall mutation rates, and compare the likelihood ratios calculated from this and from other mutation models suggested in the literature. We also show how to derive an upper bound on the likelihood ratio, depending only on the overall mutation rate.

DNA dinucleotide evolution in humans: fitting theory to facts

We examine length distributions of approximately 6000 human dinucleotide microsatellite loci, representing chromosomes 1-22, from the GDB database. Under the stepwise mutation model, results from theory and simulation are compared with the empirical data. In both constant and expanding population scenarios, a simple single-step model with parameters chosen to account for the observed variance of microsatellite lengths produces results inconsistent with the observed heterozygosity and the dispersion of length skewness. Complicating the model by allowing a variable mutation rate accounts for the homozygosity, and introducing a small probability of a large mutation step accounts for the dispersion in skewnesses. We discuss these results in light of the long-term evolution of microsatellites.

Characterization of AFLP markers in damselflies: prevalence of codominant markers and implications for population genetic applications

Amplified fragment length polymorphism (AFLP) analysis is becoming increasingly popular as a method for generating molecular markers for population genetic applications. For practical considerations, it is generally assumed in population studies that AFLPs segregate as dominant markers, i.e., that present and absent are the only possible states of a given locus. We tested the assumption of dominance in natural populations of the damselfly Nehalennia irene (Hagen) (Odonata: Coenagrionidae). Electro-blotted AFLP products from 21 samples were probed with individual markers. Eleven markers were analyzed, of which two were monomorphic and nine were polymorphic. Only two of the polymorphic markers behaved in a strictly dominant manner. The remaining seven polymorphic markers displayed various degrees of codominance, with 2-10 visible alleles in the sample. Of the three markers displaying the highest degree of variability, two contained microsatellite repeat tracts. Our results suggest that the assumption of dominance is unfounded. As a result, AFLP analysis may be unsuitable for estimating several important population genetic parameters, including genetic diversity.

Short tandem repeat (STR) haplotypes in HLA: an integrated 50-kb STR/linkage disequilibrium/gene map between the RING3 and HLA-B genes and identification of STR haplotype diversification in the class III region

We present a dense STR/linkage disequilibrium(LD)/gene map between the RING3 and HLA-B loci, reference allelic sizes on the most prevalent HLA haplotypes and their allelic frequencies in pedigree founders. This resource will facilitate LD, evolution and gene mapping studies, including comparisons of HLA and STR haplotypes and identification of HLA recombinants. The map was constructed by testing novel and previously reported STRs using a panel of 885 individuals in 211 families and 60 DNA samples from cell lines and bone marrow donors homozygous in the HLA-A, -B and -DR loci selected from over 15 000 entries into the registry of Swedish bone marrow donors. We have also analysed the variability of STR alleles/haplotypes on the most prevalent HLA haplotypes to identify STRs useful for fine mapping of disease genes in the region previously implicated in susceptibility to many disorders. The analysis of 40 HLA-A*01, B*0801, DRB1*03011, DQB1*0201 haplotypes in homozygous donors showed a surprising stability in 23 STRs between the class II recombination hot spot and HLA-B, with the average of 1.9% (16/838) variant alleles. However, 40% variant alleles were found at the D6S2670 locus in intron 19 of the tenascin-X gene both in the families and homozygous donors. The nucleotide sequence analysis of this STR showed a complex polymorphism consisting of tetra- (CTTT)(8-18) and penta-nucleotide (CTTTT)(1-2) repeats, separated by an intervening non-polymorphic sequence of 42 bp. The HLA-A1, B*0801, DRB1*03011, DQB1*0201 haplotypes had five (CTTT)(14-18)/(CTTTT)(2) variants with a predominant (CTTT)(16) allele, implicating the tetranucleotide component as the source of this ancestral haplotype diversification, which may be due to the location of D6S2670 in the region of the highest GC content in the human MHC.

Nuclear- and chloroplast-microsatellite variation in A-genome species of rice

Simple sequence length polymorphism analysis was carried out to reveal microsatellite variation and to clarify the phylogenetic relationships among A-genome species of rice. Total DNA from 29 cultivars (23 Oryza sativa and 6 O. glaberrima) and 30 accessions of wild A-genome species (12 O. rufipogon, 5 O. glumaepatula, 2 O. longistaminata, 6 O. meridionalis, and 5 O. barthii) was used as a template for PCR to detect 24 nuclear and 10 chloroplast microsatellite loci. Microsatellite allelic diversity was examined based on amplified banding patterns. Microsatellites amplified clearly in all 59 accessions, with an average of 18.4 alleles per locus. The polymorphism information content (PIC) value ranged from 0.85 to 0.94, with an average of 0.89. At the species level, high average PIC values were observed in O. sativa (0.79) and O. rufipogon (0.80). For chloroplast microsatellites, the average number of alleles per locus and the average PIC value were 2.9 and 0.38, respectively. While the magnitude of diversity was much greater for nuclear microsatellites than for chloroplast microsatellites, they showed parallel patterns of differentiation for each taxonomic group. Using the ratio of common alleles (estimated as size of amplified fragments) as a similarity index, the average percentages of common microsatellite alleles were calculated between taxa. For both nuclear and chloroplast microsatellites, O. sativa showed the highest similarity values to O. rufipogon, and O. glaberrima was most similar to O. barthii. These data support previous evidence that these cultivars originated from the corresponding wild ancestral species.Key words: simple sequence length polymorphism, SSLP, microsatellite marker, rice, Oryza sativa, allelic diversity, phylogenetics.

The forensic DNA implications of genetic differentiation between endogamous communities

In many indigenous minority populations, and among migrants from Asian and African populations now resident in western Europe, North America and Australia, there is a strong tradition of endogamy and a preference for consanguineous unions. These marriage practices can result in F(ST) values greatly in excess of the maximum value (0.01) currently recommended for forensic DNA purposes under guidelines established by the National Research Council (NRC) of the USA. To examine the possible extent of deviation from this accepted norm, three co-resident Pakistani communities were studied using 10 autosomal dinucleotide markers and six tetranucleotide markers on the Y-chromosome. The mean population subdivision coefficient (FST) value was 0.13 for the autosomal loci, and Y-chromosome loci exhibited even stronger differentiation with unique alleles identified in all three communities. The data indicate that even when sub-populations are virtually indistinguishable in terms of anthropology, geography, ethnicity or culture, they may still exhibit major genetic differentiation. Where significant population stratification is known to exist, more detailed genetic databases should be developed for forensic DNA purposes, based on reference data from each of the appropriate sub-populations and not on random or combined samples.

Microsatellite stability in human post-mortem tissues

Human identification and forensic criminal casework may involve DNA profiling of decomposed material. Somatic microsatellite (STR) instability may lead to false exclusions and theoretically to false inclusions, both in criminal cases and in human identification. Hence, the somatic and postmortal stability of the actual sequences is crucial to the reliability of such analyses. Somatic STR stability in human tissues has been documented in small series only and the effect of postmortal tissue decomposition on microsatellite stability remains to be elucidated. On this basis, we have systematically searched for somatic STR mutations in 26 deceased humans without signs of decomposition at autopsy and 25 autopsy cases with obvious signs of postmortal decomposition. A blood sample and six tissue samples were collected from each case. Seven STRs were chosen for study, the tetranucleotides HUMVWA31/A, HUMTH01, HUMF13A1, and HUMFES/FPS, and the hyperpolymorphic markers HUMAPOAI1, D11S554 and HUMACTBP2. Denaturing gel electrophoresis was performed on an ABD Prism 377 gene sequencer with Genescan 672 software (Applied Biosystems, Inc.). The bone DNA profile of each case was chosen as the standard DNA profile. All cases gave profiles from additional tissues. By intraindividual comparison of DNA profiles in the cases without signs of degradation we find that the short repetitive sequences under study are stable, that is without evidence of somatic mutations. The cases with varying degree of decomposition display postmortal microsatellite stability, we detect no somatic mutations or other possible postmortal changes that could lead to between-organ non-matches. In conclusion, PCR-based STR analyses are suitable in human identification and forensic casework dealing with different tissues, even when the substrate is heavily decomposed.

Differential distribution of simple sequence repeats in eukaryotic genome sequences

Complete chromosome/genome sequences available from humans, Drosophila melanogaster, Caenorhabditis elegans, Arabidopsis thaliana, and Saccharomyces cerevisiae were analyzed for the occurrence of mono-, di-, tri-, and tetranucleotide repeats. In all of the genomes studied, dinucleotide repeat stretches tended to be longer than other repeats. Additionally, tetranucleotide repeats in humans and trinucleotide repeats in Drosophila also seemed to be longer. Although the trends for different repeats are similar between different chromosomes within a genome, the density of repeats may vary between different chromosomes of the same species. The abundance or rarity of various di- and trinucleotide repeats in different genomes cannot be explained by nucleotide composition of a sequence or potential of repeated motifs to form alternative DNA structures. This suggests that in addition to nucleotide composition of repeat motifs, characteristic DNA replication/repair/recombination machinery might play an important role in the genesis of repeats. Moreover, analysis of complete genome coding DNA sequences of Drosophila, C. elegans, and yeast indicated that expansions of codon repeats corresponding to small hydrophilic amino acids are tolerated more, while strong selection pressures probably eliminate codon repeats encoding hydrophobic and basic amino acids. The locations and sequences of all of the repeat loci detected in genome sequences and coding DNA sequences are available at http://www.ncl-india.org/ssr and could be useful for further studies.

The use of intraallelic variability for testing neutrality and estimating population growth rate

To better understand the forces affecting individual alleles, we introduce a method for finding the joint distribution of the frequency of a neutral allele and the extent of variability at closely linked marker loci (the intraallelic variability). We model three types of intraallelic variability: (a) the number of nonrecombinants at a linked biallelic marker locus, (b) the length of a conserved haplotype, and (c) the number of mutations at a linked marker locus. If the population growth rate is known, the joint distribution provides the basis for a test of neutrality by testing whether the observed level of intraallelic variability is consistent with the observed allele frequency. If the population growth rate is unknown but neutrality can be assumed, the joint distribution provides the likelihood of the growth rate and leads to a maximum-likelihood estimate. We apply the method to data from published data sets for four loci in humans. We conclude that the Delta32 allele at CCR5 and a disease-associated allele at MLH1 arose recently and have been subject to strong selection. Alleles at PAH appear to be neutral and we estimate the recent growth rate of the European population to be approximately 0.027 per generation with a support interval of (0.017-0.037). Four of the relatively common alleles at CFTR also appear to be neutral but DeltaF508 appears to be significantly advantageous to heterozygous carriers.

Mutations at Y-STR loci: implications for paternity testing and forensic analysis

Knowledge about mutation rates and the mutational process of Y-chromosomal short-tandem-repeat (STR) or microsatellite loci used in paternity testing and forensic analysis is crucial for the correct interpretation of resulting genetic profiles. Therefore, we recently analysed a total of 4999 male germline transmissions from father/son pairs of confirmed paternity (probability > or = 99.9%) at 15 Y-STR loci which are commonly applied to forensics. We identified 14 mutations. Locus specific mutation rate estimates varied between 0 and 8.58 x 10(-3), and the overall average mutation rate estimate was 2.80 x 10(-3) (95% CIL 1.72 x 10(-3)-4.27 x 10(-3)). In two confirmed father/son pairs mutation at two Y-STRs were observed. The probability of two mutations occurring within the same single germline transmission was estimated to be statistically not unexpected. Additional alleles caused by insertion polymorphisms have been found at a number of Y-STRs and a frequency of 0.12% was estimated for DYS19. The observed mutational features for Y-STRs have important consequences for forensic applications such as the definition of criteria for exclusions in paternity testing and the interpretation of genetic profiles in stain analysis. In order to further enrich our knowledge of Y-STR mutations we suggest the establishment of a Y-STR mutation database and ask the forensic community for data contribution.

Microsatellite variation in solitary and social tuco-tucos: molecular properties and population dynamics

Variation at 15 microsatellite loci was characterized for a population of the solitary Patagonian tuco-tuco (Ctenomys haigi) and a population of group-living colonial tuco-tuco (C. sociabilis), both of which were located in the Limay River Valley of south western Argentina. All loci examined were characterized by uninterrupted di- or trinucleotide repeats in both species; seven of these loci had been isolated from C. haigi and eight had been isolated from C. sociabilis. Across all loci, there was a significant tendency for both number of alleles and heterozygosity to be greater in C. haigi than in C. sociabilis. Cloning and sequencing of multiple PCR products per locus per population revealed no significant biases in allele length, suggesting that this difference in variability was not due to ascertainment bias or to population-wide differences in rates of microsatellite evolution. Instead, differences in microsatellite variation between C. haigi and C. sociabilis were more consistent with current demographic and suspected historical differences between these populations. In particular, interpopulation differences in the distribution of microsatellite allele sizes and allele frequencies suggested that C. sociabilis has been more subject to historical reductions in population size (e.g. population bottlenecks) than has C. haigi. These findings indicate that analyses of microsatellite structure can yield important insights into the population-level phenomena that contribute to diversity at these markers, including differences in population history that continue to influence levels of genetic variability.

Microsatellite variation and fine-scale population structure in the wood frog (Rana sylvatica)

We investigated genetic population structure in wood frogs (Rana sylvatica) from a series of Prairie Pothole wetlands in the northern Great Plains. Amphibians are often thought to exist in demographic metapopulations, which require some movement between populations, yet genetic studies have revealed strong subdivision among populations, even at relatively fine scales (several km). Wood frogs are highly philopatric and studies of dispersal suggest that they may exhibit subdivision on a scale of approximately 1-2 km. We used microsatellites to examine population structure among 11 breeding assemblages separated by as little as 50 m up to approximately 5.5 km, plus one population separated from the others by 20 km. We found evidence for differentiation at the largest distances we examined and among a few neighbouring ponds, but most populations were strikingly similar in allele frequencies, suggesting high gene flow among all but the most distant populations. We hypothesize that the few significant differences among neighbouring populations at the finest scale may be a transient effect of extinction-recolonization founder events, driven by periodic drying of wetlands in this hydrologically dynamic landscape.

RAPD, RFLP and SSLP analyses of phylogenetic relationships between cultivated and wild species of rice

RAPD, RFLP, nuclear SSLP and chloroplast SSLP analyses were carried out to clarify the phylogenetic relationships among A-genome species of rice. In total, 12 cultivars of Oryza sativa (4 Japonica, 3 Javanica and 5 Indica), one cultivar of O. glaberrima, and 17 wild accessions (12 O. rufipogon, 2 O. glumaepatula, 1 O. longistaminata, 1 O. meridionalis and 1 O. barthii) were used. Their banding patterns were scored and compared to evaluate the similarity between accessions. Genetic differentiation within and between taxa was examined based on the average similarity indices. Except for chloroplast SSLP analysis, the average similarities were higher within O. sativa than within O. rufipogon, and O. sativa Indica had greater intrasubspecific variation than Japonica and Javanica. Comparisons between cultivated and wild species showed that O. sativa was closely related to O. rufipogon, while O. glaberrima was closely related to O. barthii. This indicated that two cultivated species, O. sativa and O. glaberrima, originated from O. rufipogon and O. barthii, respectively. Domestication of O. sativa seemed to be diphyletic, since strong similarity was observed between O. sativa Japonica-Javanica and O. rufipogon from China and between O. sativa Indica and O. rufipogon from tropical Asia. In addition, dendrograms for RAPD, RFLP, and nuclear and chloroplast SSLP analyses were constructed to reveal the overall genetic relationships among A-genome species. In all analyses, O. sativa and O. glaberrima formed groups with O. rufipogon and O. barthii, respectively. However, their manners of clustering with other wild species were not the same. The results of RAPD and RFLP analyses indicate that O. glumaepatula was relatively close to the groups of O. sativa and O. glaberrima whereas O. longistaminata and O. meridionalis were highly differentiated from other A-genome species. On the other hand, clear interspecific relationships were not obtained by nuclear or chloroplast SSLP analyses.

Haplotype evolution and linkage disequilibrium: A simulation study

We simulated the evolution of a three-site haplotype system, two restriction fragment length polymorphisms flanking one short tandem repeat polymorphism, under five different demographic scenarios, three with constant population size and two with population growth. The simulation was designed to observe the effects of population history, recombination fraction, and mutation rate on allele and haplotype frequencies, haplotype diversity, frequency of ancestral alleles, and linkage disequilibrium. The known ancestral haplotypes were often found at low frequencies and even became extinct after 5, 000 generations, especially with small effective population sizes. The original linkage disequilibrium was eroded and even reversed.

Worldwide genetic analysis of the CFTR region

Mutations at the cystic fibrosis transmembrane conductance regulator gene (CFTR) cause cystic fibrosis, the most prevalent severe genetic disorder in individuals of European descent. We have analyzed normal allele and haplotype variation at four short tandem repeat polymorphisms (STRPs) and two single-nucleotide polymorphisms (SNPs) in CFTR in 18 worldwide population samples, comprising a total of 1,944 chromosomes. The rooted phylogeny of the SNP haplotypes was established by typing ape samples. STRP variation within SNP haplotype backgrounds was highest in most ancestral haplotypes-although, when STRP allele sizes were taken into account, differences among haplotypes became smaller. Haplotype background determines STRP diversity to a greater extent than populations do, which indicates that haplotype backgrounds are older than populations. Heterogeneity among STRPs can be understood as the outcome of differences in mutation rate and pattern. STRP sites had higher heterozygosities in Africans, although, when whole haplotypes were considered, no significant differences remained. Linkage disequilibrium (LD) shows a complex pattern not easily related to physical distance. The analysis of the fraction of possible different haplotypes not found may circumvent some of the methodological difficulties of LD measure. LD analysis showed a positive correlation with locus polymorphism, which could partly explain the unusual pattern of similar LD between Africans and non-Africans. The low values found in non-Africans may imply that the size of the modern human population that emerged "Out of Africa" may be larger than what previous LD studies suggested.

Genetic variation at microsatellite loci in Spanish sheep

Genetic variation at 18 microsatellite loci was analysed in six indigenous Spanish sheep: Churra; Latxa; Manchega; Rasa-Aragonesa; Castellana and Merino. Merinos had frequently the highest number of alleles per locus, whereas Latxas showed the lowest one at many loci. Markers ordered decreasingly according to the number of variants differentiated in the whole population were: MAF70; TGLA13; CSSM66; BM143, BM6444; MAF36; MAF64; CSSM6; TGLA53; OarFCB11; MAF33; BM4621; MAF48; MAF65; BM1258; ILSTS002; ADCYC and OarCP34. Parameters of variability such as effective number of alleles and gene diversities corroborated the high level of variation frequently displayed by microsatellite markers. Comparison of allele distributions among populations and loci did not reveal consistent shapes. Distributions were centralised in some cases, whereas in others some kind of skewness was evident. Breed-specific alleles were detected at most loci, being frequent in Merinos and rare in Churras.

A power analysis of microsatellite-based statistics for inferring past population growth

We present results concerning the power to detect past population growth using three microsatellite-based statistics available in the current literature: (1) that based on between-locus variability, (2) that based on the shape of allele size distribution, and (3) that based on the imbalance between variance and heterozygosity at a locus. The analysis is based on the single-step stepwise mutation model. The power of the statistics is evaluated for constant, as well as variable, mutation rates across loci. The latter case is important, since it is a standard procedure to pool data collected at a number of loci, and mutation rates at microsatellite loci are known to be different. Our analysis indicates that the statistic based on the imbalance between allele size variance and heterozygosity at a locus has the highest power for detection of population growth, particularly when mutation rates vary across loci.

Molecular cloning and sequence analysis of the Danio rerio catalase gene

Catalase is an antioxidant enzyme that plays a central role in the protection against oxidative stress through the metabolism of hydrogen peroxide. Catalase has been well studied in plants, bacteria, and mammals, but little work has been done in other vertebrate species. We have cloned the zebrafish (Danio rerio) catalase cDNA containing the complete coding region and analyzed expression by both reverse transcription polymerase chain reaction and western blot. The deduced amino acid sequence predicts a protein of 526 amino acids with both the primary DNA and amino acid sequences highly conserved among vertebrate species. The major protein-heme contact points in the catalase enzyme complex are also well conserved, although several amino acids associated with the second and third levels of the major substrate channel are not, suggesting potential differences in substrate access or specificity. The 3' flanking region of the cDNA contains a dinucleotide repeat near the termination codon consisting of a near perfect CA array that is polymorphic. The rat and mouse catalase genes also contain a CA repeat sequence in the 3' untranslated region, which, along with an adjacent 5' stem-loop structure, has previously been shown to be a site for mRNA protein binding (Clerch, 1995, Arch. Biochem. Biophys. 317 (1995) 267-274). A stem-loop structure is also predicted adjacent to the zebrafish CA repeat, suggesting a similar role in catalase gene regulation.

High frequency of extrapair fertilization in a plural breeding bird, the Mexican jay, revealed by DNA microsatellites

We used tetra-nucleotide microsatellite DNA typing to estimate the frequency of extrapair fertilization (EPF) in a plural breeding species, the Mexican jay, Aphelocoma ultramarina, in Arizona. We found EPF in 32 of 51 complete broods (63%) and 55 of 139 nestlings (40%) for which the putative father had been identified (one of the highest rates of EPF known for birds). At least 96.1% of EPF fathers came from within the group. This is by far the highest known within-group EPF rate among socially monogamous, communally rearing species. Most (70%) males of breeding age (3+ years) had no genetic paternity in a given year. Social fathers (i.e. those with nests and mated females) rarely obtained EPFs; of 25 social fathers, 23 had young in only one nest and only two had young in two nests by virtue of EPF. Of the 27 males known to be EPF fathers without a nest of their own, none had young in more than one nest. Only 7% of EPF fathers had their own broods reaching banding age (day 14), compared with 29.7% of social fathers. The proportion of EPF young was significantly larger in smaller broods. Breeding females in all age classes were equally likely to have EPF young. Copyright 2000 The Association for the Study of Animal Behaviour.

Microsatellite mutations in the germline: implications for evolutionary inference

Microsatellite DNA sequences mutate at rates several orders of magnitude higher than that of the bulk of DNA. Such high rates mean that spontaneous mutations that form new-length variants can realistically be seen in pedigree analysis. Data on observed mutation events from various organisms are now accumulating, allowing inferences on DNA sequence evolution to be made through an unusually direct approach. Here I discuss and integrate microsatellite mutation data in an evolutionary context. A striking feature of the mutation process is that it seems highly heterogeneous, with distinct differences between species, repeat types, loci and alleles. Age and sex also affect the mutation rate. Within genomes at equilibrium, the microsatellite-length distribution is a delicate balance between biased mutation processes and point mutations acting towards the decay of repetitive DNA. Indeed, simple repeats do not evolve simply.

The genetic legacy of Paleolithic Homo sapiens sapiens in extant Europeans: a Y chromosome perspective

A genetic perspective of human history in Europe was derived from 22 binary markers of the nonrecombining Y chromosome (NRY). Ten lineages account for >95% of the 1007 European Y chromosomes studied. Geographic distribution and age estimates of alleles are compatible with two Paleolithic and one Neolithic migratory episode that have contributed to the modern European gene pool. A significant correlation between the NRY haplotype data and principal components based on 95 protein markers was observed, indicating the effectiveness of NRY binary polymorphisms in the characterization of human population composition and history.

Molecular cloning of a widely distributed microsatellite core sequence from the cultivated mushroom Agaricus bisporus

An Agaricus bisporus microsatellite with the tetranucleotide motif TATG tandemly repeated was isolated from an A. bisporus library enriched in repeated sequences. The use of the 16-mer oligonucleotide (TATG)4 indicates that many loci contain nearby copies of the microsatellite in opposite orientations. The wide distribution of the microsatellite in the A. bisporus genome was assessed (i) by polyacrylamide gel electrophoresis of the products generated by directed amplification of microsatellite-region DNA (DAMD) and (ii) by hybridization of these products with A. bisporus chromosomes separated by pulsed-field gel electrophoresis. This is, to our knowledge, the first microsatellite reported in the cultivated edible mushrooms. DAMD-PCR products were generated using DNA of three Pleurotus species (P. pulmonarius, P. sajor-caju, and P. florida), indicating that (TATG)4 repeats are also present in these cultivated species. The variability found within closely related strains indicates that such microsatellites are useful in fingerprinting and studying genetic variability in wild and commercial mushrooms.

Short tandem-repeat polymorphism/alu haplotype variation at the PLAT locus: implications for modern human origins

Two dinucleotide short tandem-repeat polymorphisms (STRPs) and a polymorphic Alu element spanning a 22-kb region of the PLAT locus on chromosome 8p12-q11.2 were typed in 1,287-1,420 individuals originating from 30 geographically diverse human populations, as well as in 29 great apes. These data were analyzed as haplotypes consisting of each of the dinucleotide repeats and the flanking Alu insertion/deletion polymorphism. The global pattern of STRP/Alu haplotype variation and linkage disequilibrium (LD) is informative for the reconstruction of human evolutionary history. Sub-Saharan African populations have high levels of haplotype diversity within and between populations, relative to non-Africans, and have highly divergent patterns of LD. Non-African populations have both a subset of the haplotype diversity present in Africa and a distinct pattern of LD. The pattern of haplotype variation and LD observed at the PLAT locus suggests a recent common ancestry of non-African populations, from a small population originating in eastern Africa. These data indicate that, throughout much of modern human history, sub-Saharan Africa has maintained both a large effective population size and a high level of population substructure. Additionally, Papua New Guinean and Micronesian populations have rare haplotypes observed otherwise only in African populations, suggesting ancient gene flow from Africa into Papua New Guinea, as well as gene flow between Melanesian and Micronesian populations.

Complex mutations in a high proportion of microsatellite loci from the protozoan parasite Plasmodium falciparum

Microsatellite loci are generally assumed to evolve via a stepwise mutational process and a battery of statistical techniques has been developed in recent years based on this or related mutation models. It is therefore important to investigate the appropriateness of these models in a wide variety of taxa. We used two approaches to examine mutation patterns in the malaria parasite Plasmodium falciparum: (i) we examined sequence variation at 12 tri-nucleotide repeat loci; and (ii) we analysed patterns of repeat structure and heterozygosity at 114 loci using data from 12 laboratory parasite lines. The sequencing study revealed complex patterns of mutation in five of the 12 loci studied. Alleles at two loci contain indels of 24 bp and 57 bp in flanking regions, while in the other three loci, blocks of imperfect microsatellites appear to be duplicated or inserted; these loci essentially consist of minisatellite repeats, with each repeat unit containing four to eight microsatellites. The survey of heterozygosity revealed a positive relationship between repeat number and microsatellite variability for both di- and trinucleotides, indicating a higher mutation rate in loci with longer repeat arrays. Comparisons of levels of variation in different repeat types indicate that the mutation rate of dinucleotide-bearing loci is 1.6-2.1 times faster than trinucleotides, consistent with the lower mean number of repeats in trinucleotide-bearing loci. However, despite the evidence that microsatellite arrays themselves are evolving in a manner consistent with stepwise mutation model in P. falciparum, the high frequency of complex mutations precludes the use of analytical tools based on this mutation model for many microsatellite-bearing loci in this protozoan. The results call into question the generality of models based on stepwise mutation for analysing microsatellite data, but also demonstrate the ease with which loci that violate model assumptions can be detected using minimal sequencing effort.

Instability of interstitial telomeric sequences in the human genome

The length variability of four human interstitial telomeric sequences (ITs) is described. Three of the ITs contain short telomeric stretches ranging between 53 and 84 bp and are localized in 21q22, 2q31, and 7q36; the fourth IT derives from the subtelomeric domain of chromosome 6p and contains a tract of a few hundred basepairs of exact and degenerate repeats. Using primers flanking the repeats, we amplified the genomic DNA from unrelated individuals and from family members, and we found that all the loci are polymorphic. At the 21q22 IT locus, two equally frequent alleles were found, while the number of alleles at the 2q31, 7q36, and 6pter IT loci was 8, 6, and 4, respectively. Sequence analysis revealed that in the three loci containing short ITs the alleles differ from one another for multiples of the hexanucleotide; it is likely that the mechanism leading to the polymorphism is DNA polymerase slippage. These loci were also unstable in gastric tumor cells characterized by microsatellite instability. At the 6pter IT locus, the four alleles range in length from about 500 to about 700 bp; this variability is probably due to unequal exchange or gene conversion. Our data indicate that stretches of exact internal telomeric repeats can be highly unstable, like microsatellites with shorter units, and that they can be useful polymorphic markers for linkage analysis, for forensic applications, and for the detection of genetic instability in tumors.

Microsatellite variability differs between dinucleotide repeat motifs-evidence from Drosophila melanogaster

Recently, the use of microsatellites as genetic markers has become very popular. While their evolutionary dynamics are not yet fully understood, the emerging picture is that several factors are influencing microsatellite mutation rates. Recent experiments demonstrated a significant effect of repeat motif length on microsatellite mutation rates. Here, we studied the influence of the base composition of the microsatellite. Forty-two microsatellite loci on the second chromosome with the three most abundant dinucleotide repeat motifs (TC/AG, AT/TA, GT/CA) were characterized for six different Drosophila melanogaster populations. Applying ANOVA to the variance in repeat number, we found a significant influence of repeat motif on microsatellite variability. Calculating relative mutation rates, GT/CA appears to have the highest mutation rate, and AT/TA appears to have the lowest. Similar differences in mutation rates were obtained by an alternative method which estimates microsatellite mutation rates from their genomic length distribution.

Evolution of an Alu DNA element of type Sx in the lineage of primates and the origin of an associated tetranucleotide microsatellite

A 394-bp DNA fragment, which in human is on chromosome 6 near the MOG (myelin oligodendrocyte glycoprotein) gene and encompasses an Alu element and an associated tetranucleotide microsatellite, was sequenced from a large range of primate species to follow its evolutionary divergence and to understand the origin of the microsatellite. This Alu element is found at the same orthologous position in all primates sequenced, but the tetranucleotide repeat is present only in Catarrhini between the 3'-oligo(dA) of the Alu element and the 3' flanking direct repeat. Little intraspecific variation was found. Sequence identity values for this orthologous primate Alu averaged 90% (82-99%) with transitions comprising between 70% and 100% of the observed nucleotide substitutions. Although the insertion of the Alu element predates the separation of these species, the original sequence of the site of integration can still be identified. This identification of the direct repeats suggests an active role of the oligo(dA) of the Alu element in the origin of the tetranucleotide repeats. The microsatellite probably appeared after the insertion of the Alu element, early in the lineage leading to the common ancestor of the hominoids and the Old World monkeys.

Distribution and abundance of microsatellites in the yeast genome can Be explained by a balance between slippage events and point mutations

We fit a Markov chain model of microsatellite evolution introduced by Kruglyak et al. to data on all di-, tri-, and tetranucleotide repeats in the yeast genome. Our results suggest that many features of the distribution of abundance and length of microsatellites can be explained by this simple model, which incorporates a competition between slippage events and base pair substitutions, with no need to invoke selection or constraints on the lengths. Our results provide some new information on slippage rates for individual repeat motifs, which suggest that AT-rich trinucleotide repeats have higher slippage rates. As our model predicts, we found that many repeats were adjacent to shorter repeats of the same motif. However, we also found a significant tendency of microsatellites of different motifs to cluster.

Microsatellite variation in Daphnia pulex from both sides of the Baltic Sea

Despite large genetic differentiation among neighbouring populations of many freshwater zooplankton species, a macrogeographical homogeneity of allozyme variation is generally observed. A study on breeding systems in Scandinavian populations of Daphnia pulex suggested a latitudinally related cline in breeding system with both diploid cyclic parthenogens and diploid obligate parthenogens at the latitude of 60-61 degrees N. Variation at neutral markers may be more affected by selection at linked loci in such species than in strictly sexual species. In this paper I present a study of variation at five microsatellite loci in a total of 34 populations from small ponds and rockpools on both sides of the Baltic Sea at 60-61 degrees N. Two major groups, which may represent different species of the D. pulex complex, are defined with the microsatellites. Neighbouring populations show both similar and well differentiated genetic composition. Populations separated by larger geographical distances show only a large differentiation and a macrogeographic pattern. The large differentiation observed at small distances can be explained with small effective population size: variation at the microsatellite loci has been shaped by population bottlenecks followed with expansion in size, and possibly by selection. No conclusive evidence is found for obligative parthenogenesis.

Mutational analyses of dinucleotide and tetranucleotide microsatellites in Escherichia coli: influence of sequence on expansion mutagenesis

Mutagenesis at [GT/CA](10), [TC/AG](11) and [TTCC/AAGG](9) microsatellite sequences inserted in the herpes simplex virus thymidine kinase (HSV-tk) gene was analyzed in isogenic mutL(+) and mutL(-) Escherichia coli. In both strains, significantly more expansion than deletion mutations were observed at the [TTCC/AAGG](9) motif relative to either dinucleotide motif. As the HSV-tk coding sequence contains an endogenous [G/C](7) mononucleotide repeat and approximately 1000 bp of unique sequence, we were able to compare mutagenesis among various sequence motifs. We observed that the relative risk of mutation in E.COLI: is: [TTCC/AAGG](9) > [GT/CA](10) approximately [TC/AG](11) > unique approximately [G/C](7). The mutation frequency varied 1400-fold in mutL(+) cells between the tetranucleotide motif and the mononucleotide motif, but only 50-fold in mutL(-) cells. The [G/C](7) sequence was destabilized the greatest and the tetranucleotide motif the least by loss of mismatch repair. These results demonstrate that the quantitative risk of mutation at various microsatellites greatly depends on the DNA sequence composition. We suggest alternative models for the production of expansion mutations during lagging strand replication of the [TTCC/AAGG](9) microsatellite.

Characteristics and frequency of germline mutations at microsatellite loci from the human Y chromosome, as revealed by direct observation in father/son pairs

A number of applications of analysis of human Y-chromosome microsatellite loci to human evolution and forensic science require reliable estimates of the mutation rate and knowledge of the mutational mechanism. We therefore screened a total of 4,999 meioses from father/son pairs with confirmed paternity (probability >/=99. 9%) at 15 Y-chromosomal microsatellite loci and identified 14 mutations. The locus-specific mutation-rate estimates were 0-8. 58x10-3, and the average mutation rate estimates were 3.17x10-3 (95% confidence interval [CI] 1.89-4.94x10-3) across 8 tetranucleotide microsatellites and 2.80x10-3 (95% CI 1.72-4.27x10-3) across all 15 Y-chromosomal microsatellites studied. Our data show a mutational bias toward length increase, on the basis of observation of more repeat gains than losses (10:4). The data are in almost complete agreement with the stepwise-mutation model, with 13 single-repeat changes and 1 double-repeat change. Sequence analysis revealed that all mutations occurred in uninterrupted homogenous arrays of >/=11 repeats. We conclude that mutation rates and characteristics of human Y-chromosomal microsatellites are consistent with those of autosomal microsatellites. This indicates that the general mutational mechanism of microsatellites is independent of recombination.

Human population expansion and microsatellite variation

Polymorphisms at di-, tri-, and tetranucleotide microsatellite loci have been analyzed in 14 worldwide populations. A statistical index of population expansion, denoted S(k), is introduced to detect historical changes in population size using the variation at the microsatellites. The index takes the value 0 at equilibrium with constant population size and is positive or negative according to whether the population is expanding or contracting, respectively. The use of S(k) requires estimation of properties of the mutation distribution for which we use both family data of Dib et al. for dinucleotide loci and our population data on tri- and tetranucleotide loci. Statistical estimates of the expansion index, as well as their confidence intervals from bootstrap resampling, are provided. In addition, a dynamical analysis of S(k) is presented under various assumptions on population growth or decline. The studied populations are classified as having high, intermediate, or low values of S(k) and genetic variation, and we use these to interpret the data in terms of possible population dynamics. Observed values of S(k) for samples of di-, tri-, and tetranucleotide data are compatible with population expansion earlier than 60,000 years ago in Africa, Asia, and Europe if the initial population size before the expansion was on the order of 500. Larger initial population sizes force the lower bound for the time since expansion to be much earlier. We find it unlikely that bottlenecks occurred in Central African, East Asian, or European populations, and the estimated expansion times are rather similar for all of these populations. This analysis presented here suggests that modern human populations departed from Africa long before they began to expand in size. Subsequently, the major groups (the African, East Asian, and European groups) started to grow at approximately same time. Populations of South America and Oceania show almost no growth. The Mbuti population from Zaire appears to have experienced a bottleneck during its expansion.

Heterogeneous mutation processes in human microsatellite DNA sequences

Although microsatellite polymorphisms are one of the most commonly used tools in genetic analyses, it remains to be understood how microsatellite DNA has evolved as a ubiquitous and highly abundant class of repetitive sequences in eukaryotic genomes. On the basis of analyses of spontaneous human microsatellite mutations of germline origin, I show here that different mutation biases underlie the evolution of microsatellite repeats. The within-locus mutation rate increases with allele length, but is not affected by the size difference between an individual's two alleles (allele span). Within loci, long alleles tend to mutate to shorter lengths, thereby acting to prevent infinite growth. Expansions are more common than contractions among dinucleotide repeats, whereas no such trend is evident among tetranucleotide repeats. This observation is consistent with the longer repeat lengths and higher frequency of di- compared with tetranucleotide repeats. An excess of paternally transmitted mutations (male-to-female ratio of 4.9) supports a male-biased mutation rate in the human genome.

Microsatellite mutations and inferences about human demography

Microsatellites have been widely used as tools for population studies. However, inference about population processes relies on the specification of mutation parameters that are largely unknown and likely to differ across loci. Here, we use data on somatic mutations to investigate the mutation process at 14 tetranucleotide repeats and carry out an advanced multilocus analysis of different demographic scenarios on worldwide population samples. We use a method based on less restrictive assumptions about the mutation process, which is more powerful to detect departures from the null hypothesis of constant population size than other methods previously applied to similar data sets. We detect a signal of population expansion in all samples examined, except for one African sample. As part of this analysis, we identify an "anomalous" locus whose extreme pattern of variation cannot be explained by variability in mutation size. Exaggerated mutation rate is proposed as a possible cause for its unusual variation pattern. We evaluate the effect of using it to infer population histories and show that inferences about demographic histories are markedly affected by its inclusion. In fact, exclusion of the anomalous locus reduces interlocus variability of statistics summarizing population variation and strengthens the evidence in favor of demographic growth.

The trimmed-haplotype test for linkage disequilibrium

Single-marker linkage-disequilibrium (LD) methods cannot fully describe disequilibrium in an entire chromosomal region surrounding a disease allele. With the advent of myriad tightly linked microsatellite markers, we have an opportunity to extend LD analysis from single markers to multiple-marker haplotypes. Haplotype analysis has increased statistical power to disclose the presence of a disease locus in situations where it correctly reflects the historical process involved. For maximum efficiency, evidence of LD ought to come not just from a single haplotype, which may well be rare, but in addition from many similar haplotypes that could have descended from the same ancestral founder but have been trimmed in succeeding generations. We present such an analysis, called the "trimmed-haplotype method." We focus on chromosomal regions that are small enough that disequilibrium in significant portions of them may have been preserved in some pedigrees and yet that contain enough markers to minimize coincidental occurrence of the haplotype in the absence of a disease allele: perhaps regions 1-2 cM in length. In general, we could have no idea what haplotype an ancestral founder carried generations ago, nor do we usually have a precise chromosomal location for the disease-susceptibility locus. Therefore, we must search through all possible haplotypes surrounding multiple locations. Since such repeated testing obliterates the sampling distribution of the test, we employ bootstrap methods to calculate significance levels. Trimmed-haplotype analysis is performed on family data in which genotypes have been assembled into haplotypes. It can be applied either to conventional parent-affected-offspring triads or to multiplex pedigrees. We present a method for summarizing the LD evidence, in any pedigree, that can be employed in trimmed-haplotype analysis as well as in other methods.

Integration of genetic maps by polynomial transformations

Currently available genetic maps differ in a variety of basic features; in particular, with respect to the total length of the genome. Consequently, the question arises as to the extent to which genetic maps are compatible to each other, as well as to the methods with which genetic maps can be transformed into one another. We propose a set of nonlinear, polynomial transformations that enable the integration of genetic maps at a sufficiently high overall precision. Our analysis of six major, publicly available maps, and iteratively optimized polynomials of up to degree 5, yielded differences of </= +/-0.8 cM between empirical and reconstructed marker locations for >90% of points. Similarly, we determined, at a slightly worse overall fit, those polynomials that enabled the reconstruction of sex-specific recombination estimates from sex-averaged data. Our results suggest that polynominal transformations may become a valuable extension of standard map construction methods due to a rapid integration of newly developed markers into existing maps. Am. J. Med. Genet. (Neuropsychiatr. Genet.) 96:108-113, 2000.

Microsatellite analysis of Pinus taeda L. in Zimbabwe

Deducing the origin of early 20th century introductions of Pinus taeda into Zimbabwe is possible given microsatellite markers and clear population differentiation in ancestral U.S. populations. This study was designed to determine whether P. taeda introductions into Zimbabwe came from one U.S. region or whether the present-day population is an admixture of introductions from east and west of the Mississippi River Valley. Principal components analysis, Cavalli-Sforza and Edwards' chord distances and presence of diagnostic alleles each indicate that the Zimbabwe population is an admixture. There were five novel alleles in the Zimbabwe population not represented in the indigenous U.S. populations, possibly because of de novo mutation, introgression with other introduced North American pines or sampling error.

Haplotype fine mapping by evolutionary trees

To refine the location of a disease gene within the bounds provided by linkage analysis, many scientists use the pattern of linkage disequilibrium between the disease allele and alleles at nearby markers. We describe a method that seeks to refine location by analysis of "disease" and "normal" haplotypes, thereby using multivariate information about linkage disequilibrium. Under the assumption that the disease mutation occurs in a specific gap between adjacent markers, the method first combines parsimony and likelihood to build an evolutionary tree of disease haplotypes, with each node (haplotype) separated, by a single mutational or recombinational step, from its parent. If required, latent nodes (unobserved haplotypes) are incorporated to complete the tree. Once the tree is built, its likelihood is computed from probabilities of mutation and recombination. When each gap between adjacent markers is evaluated in this fashion and these results are combined with prior information, they yield a posterior probability distribution to guide the search for the disease mutation. We show, by evolutionary simulations, that an implementation of these methods, called "FineMap," yields substantial refinement and excellent coverage for the true location of the disease mutation. Moreover, by analysis of hereditary hemochromatosis haplotypes, we show that FineMap can be robust to genetic heterogeneity.

Relative stabilities of dinucleotide and tetranucleotide repeats in cultured mammalian cells

The differences in rates of frameshift mutations between a dinucleotide repeat sequence [(CA)(17)] and a tetranucleotide repeat sequence [(GAAA)(17)] have been determined in immortalized, non-tumorigenic, mismatch repair-proficient mouse cells and in mismatch repair-defective human colorectal cancer cells. Clones with mutations were selected on the basis of restoration of activity of a bacterial neomycin resistance gene whose reading frame was disrupted by insertion of the microsatellite upstream of the translation initiation codon. This gene was introduced into the cells on a plasmid, which integrated into the genome of the host cells. Mutation rates of the tetra-nucleotide repeat were much lower than those of the dinucleotide repeat in both cell types. In addition, independent subclones of the colorectal cancer cell line were assayed by PCR for instability of endo-gen-ous tetranucleotide and dinucleotide repeat sequen-ces. In all cases, the mutation frequencies of the dinucleotide repeats were higher than those of the tetranucleotide repeats.

Trinucleotide repeats and other microsatellites in yeasts

Microsatellites are direct tandem DNA repeats found in all genomes. A particular class of microsatellites, called trinucleotide repeats, is responsible for a number of neurological disorders in humans. We review here our current state of knowledge on trinucleotide repeat instability, and discuss the molecular mechanisms that may be involved in trinucleotide repeat expansions leading to fatal diseases in humans. We also present original data on microsatellite distribution in several microbial genomes, and on the use of microsatellites as physical markers to accurately and easily genotype yeast strains.

Genetic diversity of two African and sixteen South American populations determined on the basis of six hypervariable loci

A total of 582 individuals (1,164 chromosomes) from two African, eight African-derived South American, five South American Amerindian, and three Brazilian urban populations were studied at four variable number of tandem repeat (VNTR) and two short tandem repeat (STR) hypervariable loci. These two sets of loci did not show distinct allele profiles, which might be expected if different processes promoted their molecular differentiation. The two African groups showed little difference between them, and their intrapopulational variation was similar to those obtained in the African-derived South American communities. The latter showed different degrees of interpopulation variability, despite the fact that they presented almost identical average degrees of non-African admixture. The F(ST) single locus estimates differed in the five sets of populations, probably due to genetic drift, indicating the need to consider population structure in the evaluation of their total variability. A high interpopulational diversity was found among Amerindian populations in relation to Brazilian African-derived isolated communities. This is probably a consequence of the differences in the patterns of gene flow and genetic drift that each of these semi-isolated groups experienced.