A note on the estimation of allele frequencies

A new method of obtaining allele frequency estimates is described. The method may prove useful when maximum likelihood estimates are not available. No assumptions regarding the absence of alleles are required. Errors of estimates have not been obtained, but the process seems to converge to maximum likelihood.     

Estimation of allele frequencies for VNTR loci

VNTR loci provide valuable information for a number of fields of study involving human genetics, ranging from forensics (DNA fingerprinting and paternity testing) to linkage analysis and population genetics. Alleles of a VNTR locus are simply fragments obtained from a particular portion of the DNA molecule and are defined in terms of their length. The essential element of a VNTR fragment is the repeat, which is a short sequence of basepairs. The core of the fragment is composed of a variable number of identical repeats that are linked in tandem. A sample of fragments from a population of individuals exhibits substantial variation in length because of variation in the number of repeats. Each distinct fragment length defines an allele, but any given fragment is measured with error. Therefore the observed distribution of fragment lengths is not discrete but is continuous, and determination of distinct allele classes is not straightforward. A mixture model is the natural statistical method for estimating the allele frequencies of VNTR loci. In this article we develop nonparametric methods for obtaining the distribution of allele sizes and estimates of their frequencies. Methods for obtaining maximum-likelihood estimates are developed. In addition, we suggest an empirical Bayes method to improve the maximum-likelihood estimates of the gene frequencies; the empirical Bayes procedure effects a local smoothing. The latter method works particularly well when measurement error is large relative to the repeat size, because the estimated distribution of allele frequencies when maximum likelihood is used is unreliable because of an alternating pattern of over- and underestimation. We define alleles and estimate the allele frequencies for two VNTR loci from the human genome (D17S79 and D2S44), from data obtained from Lifecodes, Inc.     

Systematic detection of epistatic interactions based on allele pair frequencies

Epistatic genetic interactions are key for understanding the genetic contribution to complex traits. Epistasis is always defined with respect to some trait such as growth rate or fitness. Whereas most existing epistasis screens explicitly test for a trait, it is also possible to implicitly test for fitness traits by searching for the over- or under-representation of allele pairs in a given population. Such analysis of imbalanced allele pair frequencies of distant loci has not been exploited yet on a genome-wide scale, mostly due to statistical difficulties such as the multiple testing problem. We propose a new approach called Imbalanced Allele Pair frequencies (ImAP) for inferring epistatic interactions that is exclusively based on DNA sequence information. Our approach is based on genome-wide SNP data sampled from a population with known family structure. We make use of genotype information of parent-child trios and inspect 3×3 contingency tables for detecting pairs of alleles from different genomic positions that are over- or under-represented in the population. We also developed a simulation setup which mimics the pedigree structure by simultaneously assuming independence of the markers. When applied to mouse SNP data, our method detected 168 imbalanced allele pairs, which is substantially more than in simulations assuming no interactions. We could validate a significant number of the interactions with external data, and we found that interacting loci are enriched for genes involved in developmental processes.     

HLA Class I allele frequencies in the Lebanese population

The highly polymorphic Human Leukocyte Antigen system encompasses different loci that have been studied in transplantation as well as diseases and population associated research. This study is the first and largest of its kind to describe the distribution of HLA-A, -B and -C alleles in Lebanon. Respectively, 1994, 1309 and 1163 Lebanese individuals referred for HLA typing and possible bone marrow/kidney donation were tested for HLA-A, HLA-B and HLA-C alleles using the polymerase chain reaction/Sequence specific priming (PCR-SSP) method. Our data were compared to that of several populations with interesting and common findings shared with the Moroccan, Jordanian, Tunisian, Omani, Korean, Chinese, Japanese, Peruan, Bulgarian, Irish, Polish, Spanish, Swiss, American, African and Brazilian populations. The following data concerning the Lebanese population will help future investigators to study the relation of HLA-A, -B and -C alleles with common diseases in Lebanon and will add to the available international literature. This new data will serve as a major reference report in the region.     

Determining SNP allele frequencies in DNA pools

Single nucleotide polymorphisms (SNPs) are among the most common types of polymorphism used for genetic association studies. A method to allow the accurate quantitation of their allele frequencies from DNA pools would both increase throughput and decrease costs for large-scale genotyping. However, to date, most DNA pooling studies have concentrated on the use of microsatellite polymorphisms. In the case of SNPs that are restriction fragment length polymorphisms (RFLPs), studies have tended to use methods for the quantitation of allele frequency from pools that rely on densitometric evaluation of bands on an autoradiograph. Radiation-based methods have well-known drawbacks, and we present two alternative methods for the determination of SNP allele frequencies. For RFLPs, we used agarose gel electrophoresis of digested PCR products with ethidium bromide staining combined with densitometric analysis of gel images on a PC. For all types of SNP, we used allele-specific fluorescent probes in the Taqman assay to determine the relative frequencies of two different alleles. Both methods gave accurate and reproducible results, suggesting they are suitable for use in DNA pooling experiments.     

HLA class II allele frequencies in the Lebanese population

Aims

Being one of the most polymorphic genetic systems , the Human Leukocyte Antigen system is divided into class I (HLA-A, HLA-B and HLA-C) and class II (HLA-DP, -DQ and -DR). This study is the first and largest of its kind to describe the distribution of HLA-DQB1 and HLA-DRB1 alleles in Lebanon and the region.

Methods

Respectively, 560 and 563 Lebanese individuals referred for HLA typing and possible bone marrow/kidney donation were tested for HLA-DQB1 and HLA-DRB1 alleles using the polymerase chain reaction/sequence specific priming (PCR-SSP) method.

Results

Our data were compared to that of several populations with interesting common findings between the Lebanese, Jordanian, Bahraini, Saudi, Kuwaiti, Tunisian, Korean, Japanese, Thai, Irish, Bulgarian and Polish populations.

Conclusion

These data about the Lebanese population are going to aid future researchers to study the relation of HLA-DQB1 and HLA-DRB1 alleles with major and common diseases in the Lebanese population and will add to the available international literature associated with these loci. In addition it will serve as a reference for the future national bone marrow registry program in our country. We also reviewed the literature for the described association between HLA-DRB1 and -DQB1 loci and different disease entities.     

Inferring the genetic basis of inbreeding depression in plants.

Recent progress in the genetic analysis of inbreeding depression in plants is reviewed. While the debate over the importance of genes of dominance versus overdominance effect continues, the scope of inferences has widened and now includes such facets as the interactions between genes, the relative abundance of major versus minor genes, life cycle stage expression, and mutation rates. The types of inferences are classified into the genomic, where many genes are characterized as an average, and the genic, where individual genes are characterized. Genomic inferences can be based upon natural levels of inbreeding depression, purging experiments, the comparison of individuals of differing F (e.g., prior inbreeding), and various crossing designs. Genic inferences mainly involve mapping and characterizing loci with genetic markers, involving either a single cross or, ideally, several crosses. Alternative statistical models for analyzing polymorphic loci causing inbreeding depression should be a fruitful problem for geneticists to pursue. Key words : inbreeding depression, genetic load, self-fertilization, QTL mapping.     

DAT1 VNTR allele frequencies in the Omani population

To determine the genetic variability at the human DAT1 VNTR locus, we screened 110 healthy Omani blood donors using the polymerase chain reaction (PCR) and agarose gel electrophoresis. Two common alleles (DAT1*9 and DAT1*10) were observed with a frequency of 0.332 and 0.609, respectively. There were also five rare alleles (*3, *6, *7, *8, and *11). The frequency of the observed genotypes was not significantly different from the expected Hardy-Weinberg distribution. We compared the Omani DAT1 alleles with similar data from other populations.     

Estimation of allele frequencies in ethnically heterogeneous populations

A method for reconstructing allele frequencies characteristic of an original ethnically homogeneous population before the start of migration processes is described. Information on both the ethnic group studied and offspring of interethnic marriages is used to estimate the allele frequencies. This makes it possible to increase the informativeness of the sample, which, in the case of ethnic heterogeneity, depends not only on allele frequencies and the total sample size, but also on the ethnic structure of the sample. The problem of estimating allele frequency in an ethnically heterogeneous sample has been solved analytically for diallelic loci. It has been demonstrated that, if offspring of interethnic marriages with the same degree of outbreeding is added to a sample of the ethnic group studied, the sample informativeness does not change. To utilize the information contained in the phenotypes of the offspring of interethnic marriages, representatives of the population from which migration occurs should be included into the sample. The size of the sample ensuring the preassigned accuracy of estimation is minimized at a certain ratio between the numbers of the offspring of interethnic marriages and the "immigrants." To analyze polyallelic loci, a software package has been developed that allows estimating allele frequencies, determining the errors of these estimates, and planning the sample ensuring the preassigned accuracy of estimation. The package is available free at http://mga.bionet.bionet.nsc.ru/PopMixed/PopMixed.html.     

Microsatellite allele frequencies in humans and chimpanzees, with implications for constraints on allele size

The distributions of allele sizes at eight simple-sequence repeat (SSR) or microsatellite loci in chimpanzees are found and compared with the distributions previously obtained from several human populations. At several loci, the differences in average allele size between chimpanzees and humans are sufficiently small that there might be a constraint on the evolution of average allele size. Furthermore, a model that allows for a bias in the mutation process shows that for some loci a weak bias can account for the observations. Several alleles at one of the loci (Mfd 59) were sequenced. Differences between alleles of different lengths were found to be more complex than previously assumed. An 8-base-pair deletion was present in the nonvariable region of the chimpanzee locus. This locus contains a previously unrecognized repeated region, which is imperfect in humans and perfect in chimpanzees. The apparently greater opportunity for mutation conferred by the two perfect repeat regions in chimpanzees is reflected in the higher variance in repeat number at Mfd 59 in chimpanzees than in humans. These data indicate that interspecific differences in allele length are not always attributable to simple changes in the number of repeats.      

Mutant allele frequencies in domestic cats of Taiwan

Preliminary data for Taiwanese cats generally agree with previous findings in far eastern populations, especially Vladivostok. However, surveys are still too few in number to achieve any detailed perspective for this vast region.     

On the distribution of allele frequencies in a diffusion model

An expression is derived and values tabulated for the expected allele frequencies and their variances, arranged in decreasing order in a population, from the finite and infinite alleles diffusion model in Watterson (1976). The neutral model and also a model with heterozygote selection are considered. Some observed ABO blood group allele frequencies are compared with the tabulated expected frequencies in the neutral three allele model. This extends the results of Watterson and Guess (1977) who tabulate the expected value of the most common allele. One test of neutrality previously advocated is to consider the distribution of F, the population homozygosity, conditional on G, the product of allele frequencies. However it is shown here that for a large number of alleles, F and G are asymptotically independent, the test would not be a good one in this case. A limit theorem is derived for the distribution of allele frequencies in the neutral model when the mutation rate is large. In this case F is shown to be asymptotically normal. An inequality is derived for the probability that the oldest allele in a population is amongst the r most frequent types. An inequality is also found for the probability that a sample will only contain representatives of the r most frequent allele types in the population.     

A method for adjusting allele frequencies in the case of microsatellite allele drop‐out

A new method is proposed to adjust allele frequencies when allelic drop‐out is common. This method assumes Hardy–Weinberg equilibrium (HWE), and treats the problematic alleles as a one‐locus two‐allele system with dominance. By assuming that the homozygote frequency of the ‘recessive’ allele is measured correctly, we can back calculate the allele frequency of the ‘dominant’ allele, and adjust the heterozygote frequency accordingly. The drawback is that multilocus genotypes cannot be constructed and tests that use deviations from Hardy–Weinberg such as tests for bottlenecks become impossible. An example is given where a large homozygote excess (F_IS = 0.44) is adjusted to a reasonable level (F_IS = 0.046). The effect of scoring error was set in relation to sampling error and while F_IS values can be seriously biased, F_ST values are not necessarily so, if scoring error and sample size are both low. As sample size increases, the effect of scoring error increases.