Vntr Allele Frequency Distributions under the Stepwise Mutation Model: A Computer Simulation Approach

Variable numbers of tandem repeats (VNTRs) are a class of highly informative and widely dispersed genetic markers. Despite their wide application in biological science, little is known about their mutational mechanisms or population dynamics. The objective of this work was to investigate four summary measures of VNTR allele frequency distributions: number of alleles, number of modes, range in allele size and heterozygosity, using computer simulations of the one-step stepwise mutation model (SMM). We estimated these measures and their probability distributions for a wide range of mutation rates and compared the simulation results with predictions from analytical formulations of the one-step SMM. The average heterozygosity from the simulations agreed with the analytical expectation under the SMM. The average number of alleles, however, was larger in the simulations than the analytical expectation of the SMM. We then compared our simulation expectations with actual data reported in the literature. We used the sample size and observed heterozygosity to determine the expected value, 5th and 95th percentiles for the other three summary measures, allelic size range, number of modes and number of alleles. The loci analyzed were classified into three groups based on the size of the repeat unit: microsatellites (1-2 base pair (bp) repeat unit), short tandem repeats [(STR) 3-5 bp repeat unit], and minisatellites (15-70 bp repeat unit). In general, STR loci were most similar to the simulation results under the SMM for the three summary measures (number of alleles, number of modes and range in allele size), followed by the microsatellite loci and then by the minisatellite loci, which showed deviations in the direction of the infinite allele model (IAM). Based on these differences, we hypothesize that these three classes of loci are subject to different mutational forces.     

A Recent Shift from Polygyny to Monogamy in Humans Is Suggested by the Analysis of Worldwide Y-Chromosome Diversity

Molecular genetic data contain information on the history of populations. Evidence of prehistoric demographic expansions has been detected in the mitochondrial diversity of most human populations and in a Y-chromosome STR analysis, but not in a previous study of 11 Y-chromosome SNPs in Europeans. In this paper, we show that mismatch distributions and tests of mutation/drift equilibrium based on up to 166 Y-chromosome SNPs, in 46 samples from all continents, also fail to support an increase of the male effective population size. Computer simulations show that the low nuclear versus mitochondrial mutation rates cannot explain these results. However, ascertainment bias, i.e., when only highly variable SNP sites are typed, may be concealing any Y SNPs evidence for a recent, but not an ancient, increase in male effective population sizes. The results of our SNP analyses can be reconciled with the expansion of male effective population sizes inferred from STR loci, and with mitochondrial evidence, by admitting that humans were essentially polygynous during much of their history. As a consequence, until recently only a few men may have contributed a large fraction of the Y-chromosome pool at every generation. The number of breeding males may have increased, and the variance of their reproductive success may have decreased, through a recent shift from polygyny to monogamy, which is supported by ethnological data and possibly accompanied the shift from mobile to sedentary communities.     

Towards Improvements in the Estimation of the Coalescent: Implications for the Most Effective Use of Y Chromosome Short Tandem Repeat Mutation Rates

Over the past two decades, many short tandem repeat (STR) microsatellite loci on the human Y chromosome have been identified together with mutation rate estimates for the individual loci. These have been used to estimate the coalescent age, or the time to the most recent common ancestor (TMRCA) expressed in generations, in conjunction with the average square difference measure (ASD), an unbiased point estimator of TMRCA based upon the average within-locus allele variance between haplotypes. The ASD estimator, in turn, depends on accurate mutation rate estimates to be able to produce good approximations of the coalescent age of a sample. Here, a comparison is made between three published sets of per locus mutation rate estimates as they are applied to the calculation of the coalescent age for real and simulated population samples. A novel evaluation method is developed for estimating the degree of conformity of any Y chromosome STR locus of interest to the strict stepwise mutation model and specific recommendations are made regarding the suitability of thirty-two commonly used Y-STR loci for the purpose of estimating the coalescent. The use of the geometric mean for averaging ASD and across loci is shown to improve the consistency of the resulting estimates, with decreased sensitivity to outliers and to the number of STR loci compared or the particular set of mutation rates selected.     

Allele Frequencies at Microsatellite Loci: The Stepwise Mutation Model Revisited

We summarize available data on the frequencies of alleles at microsatellite loci in human populations and compare observed distributions of allele frequencies to those generated by a simulation of the stepwise mutation model. We show that observed frequency distributions at 108 loci are consistent with the results of the model under the assumption that mutations cause an increase or decrease in repeat number by one and under the condition that the product Nu, where N is the effective population size and u is the mutation rate, is larger than one. We show that the variance of the distribution of allele sizes is a useful estimator of Nu and performs much better than previously suggested estimators for the stepwise mutation model. In the data, there is no correlation between the mean and variance in allele size at a locus or between the number of alleles and mean allele size, which suggests that the mutation rate at these loci is independent of allele size.     

Genetic variation at five trimeric and tetrameric tandem repeat loci in four human population groups.

Trimeric and tetrameric short tandem repeats (STRs) represent a rich source of highly polymorphic markers in the human genome that may be studied with the polymerase chain reaction (PCR). We report the analysis of a multilocus genotype survey of 97-380 chromosomes in U.S. Black, White, Mexican-American, and Asian populations at five STR loci located on chromosomes 1, 4, 11, and X. The heterozygote frequencies of the loci ranged from 0.36 to 0.91 and the number of alleles from 6 to 20 for the 20 population and locus combinations. Relative allele frequencies exhibited differences between populations and unimodal, bimodal, and complex distributions. Although deviations were noted at some locus-population test combinations, genotype data from the loci were consistent overall with Hardy-Weinberg equilibrium by three tests. Population subheterogeneity within each ethnic group was not detected by two additional tests. No mutations were detected in a total of 860 meioses for two loci studied in the CEPH kindreds and five loci studied in other families. An indirect estimate of the mutation rates gave values from 2.3 x 10(-5) to 15.9 x 10(-5) for the five loci. Higher mutation rates appear to be associated with greater numbers of tandem repeats in the core motif. The most frequent genotype for all five loci combined appears to have a frequency of 7.59 x 10(-4). Together, these results suggest that trimeric and tetrameric STR loci are useful markers for the study of new mutations and genetic linkage analysis and for application to personal identification in the medical and forensic sciences.     

Population genetics of short tandem repeat (STR) loci

To investigate the population genetics of short tandem repeat (STR) polymorphisms in human populations, we have studied the allele frequency distributions of four STR loci (HUMTH01, HUMVWA31, HUMF13A1 and HUMFES) in 16 different population surveys which can be categorised within three broadly defined ethnic groups: Caucasian, Asian (Indian subcontinent), and African (Afro-Caribbean and US black). We have observed that allele frequency distributions of populations within ethnic groups are similar; consequently, genetic distances are an order of magnitude lower than between ethnic groups. Inbreeding coefficients (F-statistics) and calculations of the number of mean heterozygous loci per individual, along with estimates of variance, did not suggest that the populations were substructured. This included a study of an immigrant Asian population known to comprise at least three different sub-groups. Finally, an indication of the discriminating power is given by calculation of likelihood ratios (LR) of each individual tested across all four loci. Approximately 70% of Caucasians give an LR of greater than 10,000; the test is even more discriminating in Afro-Caribbeans-approximately 90% of tests are greater than 10,000.Editor's commentsThe authors present data generated by the move from VNTR to STR loci for human identification. The data they present for samples within major racial groupings address some of the concerns about population substructuring discussed by Balding and Nichols in this volume.     

广西融水苗族3个STR基因座的群体遗传学研究

为了解广西融水苗族人群无关个体的3个短串联重复序列(short tandem repeat,STR):HUMCSF1PO,HUMTPOX,HUMTH01遗传多态性分布情况,本文用枸橼酸钠抗凝法采集血样,酚—氯仿抽提法提取DNA,应用复合扩增技术对血样DNA的3个STR基因座进行扩增和检测。结果显示:在三个STR位点共检测出19种等位基因,48种基因型,频率分布分别在0.0024—0.4663和0.0048—0.3173之间;基因型的分布符合Hardy-Weinberg平衡定律(P>0.05)。计算种族、民族之间的遗传距离,并对之进行比较得出:广西融水苗族与美国高加索人及美国非洲人存在显著差异,且与美国非洲人之间的差异大于与美国高加索人之间的差异;广西融水苗族与广西侗族的关系近于与其他少数民族的关系。     

Analysis of allele frequency of seven microsatellite loci of Y chromosome in three Tuva populations

The allele frequency distribution of seven microsatellite loci of the nonrecombining region of the Y chromosome (Y-STRs) was analyzed in three geographically distant indigenous populations of the Tuva Republic. The populations did not differ in allele frequency distribution of the seven Y-STRs. The Y-chromosome microsatellite loci in Tuvinians showed a high diversity (H = 0.575) that was nearly identical in all three populations. The genetic distance Ddm between the three populations was low, suggesting no subdivision of the modern male population of Tuva. Estimates of the period of linear changes in Ddm showed that Y-chromosome microsatellites can be used to reconstruct evolutionary events dating back no more than 40,000-50,000 years. The problems of human population phylogeny are discussed on the basis of data on Y-chromosome STRs.     

Detecting population growth, selection and inherited fertility from haplotypic data in humans

The frequency of a rare mutant allele and the level of allelic association between this allele and one or several closely linked markers are frequently measured in genetic epidemiology. Both quantities are related to the time elapsed since the appearance of the mutation in the population and the intrinsic growth rate of the mutation (which may be different from the average population growth rate). Here, we develop a method that uses these two kinds of genetic data to perform a joint estimation of the age of the mutation and the minimum growth rate that is compatible with its present frequency. In absence of demographic data, it provides a useful estimate of population growth rate. When such data are available, contrasts among estimates from several loci allow demographic processes, affecting all loci similarly, to be distinguished from selection, affecting loci differently. Testing these estimates on populations for which data are available for several disorders shows good congruence with demographic data in some cases whereas in others higher growth rates are obtained, which may be the result of selection or hidden demographic processes.     

Mutation analysis of 28 autosomal short tandem repeats in the Chinese Han population

Short tandem repeats (STRs) have been extensively used in forensic genetics. However, according to previous studies, the mutation rates of STRs are relatively high and are affected by many factors. Therefore, it is important to analyze STR mutations and determine the influence of underlying factors on STR mutation rates. Mutation rates of 28 autosomal STRs were determined from 8708 paternity testing cases in the Chinese Han population, and the relationships between STR mutation rates and population, sex, age, allele length and heterozygosity were investigated. A total of 279 mutations were observed at 27 loci in a total of 233,530 meiosis cases, including 273 (97.8%) one-step, 5 (1.8%) two-step and 1 (0.4%) three-step mutations. The overall average mutation rate was 1.19?×?10^–3 (95% CI 1.06?×?10^–3???1.34?×?10^–3) ranging from 0 (TPOX) to 2.79?×?10^–3 (D13S325). Mutation rate comparisons revealed statistically significant differences at several STRs among populations. Paternal mutations occurred more frequently than maternal mutations, at a ratio of 6.04:1, and the mutation rate tended to increase with paternal age. Moreover, our study revealed a bias towards contraction mutations for long alleles and expansion mutations for short alleles. No obvious bias was observed in the overall mutation direction. In addition, STR loci with higher expected heterozygosity (H_exp) tended to have higher mutation rates. This work revealed the relationships between STR mutation rates and several influencing factors, providing useful data and information for further research on STR mutations in forensic genetics.

     

Comparative study of genetic variation at 15 STR loci in three isolated populations of the Bosnian mountain area

Fifteen autosomal STR loci (D3S1358, TH01, D21S11, D18S51, Penta E, D5S818, D13S317, D7S820, D16S539, CSF1PO, Penta D, VWA, D8S1179, TPOX, and FGA) were studied in three geographically close but isolated populations from the Bosnian mountain area. The three villages are Bobovica, Dejcici, and Lukomir. DNA was obtained from 83 individuals, and the allele frequencies and genetic diversity among the three sample groups were compared. In addition, seven of the STR loci (CSF1PO, D13S317, D3S1358, D5S818, D7S820, FGA, TH01) were used in a comparative population analysis of the Bjelasnica-Treskavica region and the Adriatic islands of Brac, Hvar, and Korcula. Although the sample sizes are relatively small, the observed variation within any of the small isolated populations is high and comparable to less isolated groups. In addition, even though the populations are geographically isolated, the STR data are similar among the populations. The most significant frequency differences were observed at the TH01 locus. Although the specific allele distributions in any untyped population cannot be determined a priori, we find support for a high degree of diversity for the STR loci in most populations. In addition, the multiple locus profile is highly informative not only for various population studies but also for forensic studies, even when specific population data are not available.     

Genetic variation of new 21 autosomal short tandem repeat loci in a Chinese Salar ethnic group

In the present study, we reported the allele frequencies for new 21 autosomal short tandem repeat (STR) loci, including D6S474, D12ATA63, D22S1045, D10S1248, D1S1677, D11S4463, D1S1627, D3S4529, D2S441, D6S1017, D4S2408, D19S433, D17S1301, D1GATA113, D18S853, D20S482, D14S1434, D9S1122, D2S1776, D10S1435 and D5S2500 loci. Forensic statistical parameters were estimated from a sample set of 120 unrelated healthy individuals from the Salar ethnic group in Xunhua Salar Autonomous County of Qinghai province, China. A total of 151 alleles were observed at 21 STR loci in the population, and their allele frequencies were in the range of 0.004–0.554. All STR loci showed a high degree of genetic polymorphisms, and the combined probability of exclusion, combined power of discrimination and combined probability of matching for all 21 STR loci were 0.9999993134, 0.99999999999999999991739 and 8.2607 × 10⁻²⁰, respectively. For all the 21 STR loci in the Salar ethnic group, the observed genotypic data showed no significant deviation from those expected under the Hardy-Weinberg equilibrium. The allele frequency distributions for the 21 autosomal STR loci were compared between the Salar group and its neighboring populations and significant differences were detected among these populations at D1S1677, D2S441, D3S4529, D4S2408, D6S1017, D11S4463, D12ATA63, D14S14343, D18S853, D19S433 and D22S1045 loci.     

Precision and accuracy of divergence time estimates from STR and SNPSTR variation

Inference of intraspecific population divergence patterns typically requires genetic data for molecular markers with relatively high mutation rates. Microsatellites, or short tandem repeat (STR) polymorphisms, have proven informative in many such investigations. These markers are characterized, however, by high levels of homoplasy and varying mutational properties, often leading to inaccurate inference of population divergence. A SNPSTR is a genetic system that consists of an STR polymorphism closely linked (typically < 500 bp) to one or more single-nucleotide polymorphisms (SNPs). SNPSTR systems are characterized by lower levels of homoplasy than are STR loci. Divergence time estimates based on STR variation (on the derived SNP allele background) should, therefore, be more accurate and precise. We use coalescent-based simulations in the context of several models of demographic history to compare divergence time estimates based on SNPSTR haplotype frequencies and STR allele frequencies. We demonstrate that estimates of divergence time based on STR variation on the background of a derived SNP allele are more accurate (3% to 7% bias for SNPSTR versus 11% to 20% bias for STR) and more precise than STR-based estimates, conditional on a recent SNP mutation. These results hold even for models involving complex demographic scenarios with gene flow, population expansion, and population bottlenecks. Varying the timing of the mutation event generating the SNP revealed that estimates of divergence time are sensitive to SNP age, with more recent SNPs giving more accurate and precise estimates of divergence time. However, varying both mutational properties of STR loci and SNP age demonstrated that multiple independent SNPSTR systems provide less biased estimates of divergence time. Furthermore, the combination of estimates based separately on STR and SNPSTR variation provides insight into the age of the derived SNP alleles. In light of our simulations, we interpret estimates from data for human populations.     

Efficient DNA fingerprinting method for the identification of cross-culture contamination of cell lines

In order to identify cross-culture contamination of cell lines, we applied DNA fingerprinting using variable number of tandem repeat (VNTR) loci and short tandem repeat (STR) loci amplified by polymerase chain reaction (PCR) instead of a radioisotope labeled multilocus probe. Eleven cell lines were used for the Apo B and D1S80 loci detection, and twelve cell lines were examined in the Y-chromosome analysis. The data obtained from the sister cell lines NALM-6 and B85, two MOLM-1 cultures from two cryopreserved tubes, and four subclones of BALM-9 and its sister cell line BALM-10, displayed clear and distinct bands of each PCR product for both Apo B and D1S80. Detection of a Y-chromosome DNA sequence is another very informative marker for the identification of cell lines, if the Y-chromosome is present. We examined eight cell lines for the expression of four STR loci; the data thus generated were compared with the results previously reported from other laboratories. The resulting electrophoretic banding patterns showed that our "home-made" STR detection system is a useful and efficient tool for the authentication of cell lines. PCR detection of VNTR and STR loci represents a simple, rapid and powerful DNA fingerprinting technique to authenticate human cell lines and to detect cross-culture contamination. This PCR technique may be used in lieu of the more time-consuming, labor-intensive and radioactive Southern blot multilocus method.     

CPI distribution and cutoff values for duo kinship testing

DNA-based tests commonly use 13 STR (short tandem repeat) loci in human identification and paternity testing--the Combined DNA Index System or CODIS. Its average degree of accuracy of paternity identification is greater than 0.9999 under the circumstance of a mother, a child and a putative father. However, the possibility of false inclusions increases under circumstances such as [1] only two members of a family group are available--a duo case during determination of paternity or [2] identification of human remains while only one living relative is present. In Taiwan, the National Unidentified Human Remains Database uses the CODIS 13 STR for the identification of family members. Two or more reference samples in the DNA database have been found to share one allele at all loci tested. Then the Combined Paternity Index (CPI) is used to determine and provide an estimate of kinship in such cases. Combining 499,500 sets of DNA data for the 13 STR CODIS loci, totally 431 (0.086%) cases are false inclusions where all 13 loci shared at least one allele. Simulated partial DNA profiles (not all 13 loci yielded results) were created to mimic the mutation and degradation process. All 431 real duo cases were analyzed to evaluate sensitivity and specificity. This report provided four kinship-matching situations with CPI cutoff values when the number of allele-sharing loci exceeded 11. CPI values greater or lesser than the suggested cutoff point will provide a greater degree of confidence in determining whether two samples are derived from first-degree relatives.     

Out-of-Africa, the peopling of continents and islands: tracing uniparental gene trees across the map

Genetic relationships between human groups were first studied by comparisons of relative allele frequency at multiple loci. Geographical study of detailed, highly resolved trees of single, non-recombining uniparental loci (mitochondrial DNA: mtDNA and Y chromosome/non-recombining Y: NRY), following specific lineages rather than populations, then revolutionized knowledge of the peopling of the world, although, curiously, the use of geographically highly specific mutations that protect against malaria, found on individual autosomal globin genes, were first in single-locus phylogeography. mtDNA, with its high single nucleotide polymorphism (SNP) mutation rates and relative ease of dating, led the way and gave stronger proof of the recent near replacement of all human species by anatomically modern humans (AMH). AMH left Africa via a single southern exit about 70 000 years ago and rapidly spread around the Indian Ocean towards the Antipodes, long before a small branch left a South Asian colony, earlier on the trail, to populate Europe. The worldwide skeleton phylogeny of mtDNA is fully resolved, but a regional analysis will continue to illuminate subsequent migrations. NRY with a lower SNP mutation rate still has a dating problem relating to use the of single tandem repeats (STRs), but has validated mtDNA results and with more geographical specificity and genomic size, as with the autosomal human genome, has much more detail to offer for the future.     

Variability of fifteen autosomal microsatellite DNA loci in five populations of aboriginal South Siberians

The allele distributions for 15 STR loci included in the AmpFISTR SGM Plus and AmpFISTR Profiler Plus kits ("Applied Biosystems", USA) were determined in 261 healthy unrelated individuals belonging to five indigenous populations of South Siberia: in Buryats, Altaians, Tofalars, Sojots and Khakassians. No significant differences in allele frequencies were found between populations studied. Combined power of discrimination (PD) for the STR loci investigated were estimated for the populations under study.     

Genetic analysis of the Utah population: a comparison of STR and VNTR loci

Genetic data are reported for nine short tandem repeat (STR) loci (D3S1358, vWA, FGA, D8S1179, D21S11, D18S51, D5S818, D13S317, and D7S820) and six variable number of tandem repeat (VNTR) loci (D2S44, D10S28, D4S139, D1S7, D5S110, and D17S79) in samples of Utah African Americans, European Americans, and Hispanics. Little evidence of departures from Hardy-Weinberg equilibrium or gametic equilibrium was found in these populations. Because of their relatively higher mutation rates, the VNTR loci exhibited higher average heterozygosity and lower FST levels than did the STR loci. Genetic distance analysis showed congruence between the two types of systems, and a genetic distance analysis of the STR data showed that the three Utah populations are genetically similar to the same ethnic groups in other parts of the United States. In addition, this analysis showed that the African American population is the most genetically divergent, with greater similarity between the Hispanic and European American populations. This analysis demonstrates a high degree of consistency for population designations commonly used in forensic analysis.     

Uniparental genetic markers in South Amerindians

A comprehensive review of uniparental systems in South Amerindians was undertaken. Variability in the Y-chromosome haplogroups were assessed in 68 populations and 1,814 individuals whereas that of Y-STR markers was assessed in 29 populations and 590 subjects. Variability in the mitochondrial DNA (mtDNA) haplogroup was examined in 108 populations and 6,697 persons, and sequencing studies used either the complete mtDNA genome or the highly variable segments 1 and 2. The diversity of the markers made it difficult to establish a general picture of Y-chromosome variability in the populations studied. However, haplogroup Q1a3a* was almost always the most prevalent whereas Q1a3* occurred equally in all regions, which suggested its prevalence among the early colonizers. The STR allele frequencies were used to derive a possible ancient Native American Q-clade chromosome haplotype and five of six STR loci showed significant geographic variation. Geographic and linguistic factors moderately influenced the mtDNA distributions (6% and 7%, respectively) and mtDNA haplogroups A and D correlated positively and negatively, respectively, with latitude. The data analyzed here provide rich material for understanding the biological history of South Amerindians and can serve as a basis for comparative studies involving other types of data, such as cultural data.     

Mutation rate in human microsatellites: influence of the structure and length of the tandem repeat.

In 10,844 parent/child allelic transfers at nine short-tandem-repeat (STR) loci, 23 isolated STR mismatches were observed. The parenthood in each of these cases was highly validated (probability >99.97%). The event was always repeat related, owing to either a single-step mutation (n=22) or a double-step mutation (n=1). The mutation rate was between 0 and 7 x 10(-3) per locus per gamete per generation. No mutations were observed in three of the nine loci. Mutation events in the male germ line were five to six times more frequent than in the female germ line. A positive exponential correlation between the geometric mean of the number of uninterrupted repeats and the mutation rate was observed. Our data demonstrate that mutation rates of different loci can differ by several orders of magnitude and that different alleles at one locus exhibit different mutation rates.