Nondetectability of restriction fragments and independence of DNA fragment sizes within and between loci in RFLP typing of DNA.

We provide experimental evidence showing that, during the restriction-enzyme digestion of DNA samples, some of the HaeIII-digested DNA fragments are small enough to prevent their reliable sizing on a Southern gel. As a result of such nondetectability of DNA fragments, individuals who show a single-band DNA profile at a VNTR locus may not necessarily be true homozygotes. In a population database, when the presence of such nondetectable alleles is ignored, we show that a pseudodependence of alleles within as well as across loci may occur. Using a known statistical method, under the hypothesis of independence of alleles within loci, we derive an efficient estimate of null allele frequency, which may be subsequently used for testing allelic independence within and across loci. The estimates of null allele frequencies, thus derived, are shown to agree with direct experimental data on the frequencies of HaeIII-null alleles. Incorporation of null alleles into the analysis of the forensic VNTR database suggests that the assumptions of allelic independence within and between loci are appropriate. In contrast, a failure to incorporate the occurrence of null alleles would provide a wrong inference regarding the independence of alleles within and between loci.     

Measures of association between disease and genotype.

This paper is concerned with the statistical aspects of the phenomenon of disease occurring more frequently in individuals with some genotypes than in individuals with others. A correlation coefficient is defined to quantify association between disease and genotype. A distinction is made between the concepts of independence of allele and disease and independence of genotype and disease. This distinction is used to define two components of association which describe separate aspects of association of disease with genotype. One component is a measure of the association of disease with allele; the other a measure of the effect of allele interaction on association of disease and genotype. One aspect of the usefulness of the partition into components which is discussed is in expressing the recurrence risk of disease for a relative of an affected individual. A chi-squared analysis is provided to test hypotheses about the components of association and other hypotheses of genetic interest. This analysis is illustrated using a study done to determine the effect of the sex-linked dwarfing gene in male chickens on resistance to E. coli infection. This analysis shows a significant allele interaction effect on resistance to disease but no association of disease with alleles. In conclusion, some extensions and limitations of the proposed concepts and procedures are discussed.     

Correlation-based inference for linkage disequilibrium with multiple alleles

The correlation between alleles at a pair of genetic loci is a measure of linkage disequilibrium. The square of the sample correlation multiplied by sample size provides the usual test statistic for the hypothesis of no disequilibrium for loci with two alleles and this relation has proved useful for study design and marker selection. Nevertheless, this relation holds only in a diallelic case, and an extension to multiple alleles has not been made. Here we introduce a similar statistic, R(2), which leads to a correlation-based test for loci with multiple alleles: for a pair of loci with k and m alleles, and a sample of n individuals, the approximate distribution of n(k - 1)(m - 1)/(km)R(2) under independence between loci is chi((k-1)(m-1))(2). One advantage of this statistic is that it can be interpreted as the total correlation between a pair of loci. When the phase of two-locus genotypes is known, the approach is equivalent to a test for the overall correlation between rows and columns in a contingency table. In the phase-known case, R(2) is the sum of the squared sample correlations for all km 2 x 2 subtables formed by collapsing to one allele vs. the rest at each locus. We examine the approximate distribution under the null of independence for R(2) and report its close agreement with the exact distribution obtained by permutation. The test for independence using R(2) is a strong competitor to approaches such as Pearson's chi square, Fisher's exact test, and a test based on Cressie and Read's power divergence statistic. We combine this approach with our previous composite-disequilibrium measures to address the case when the genotypic phase is unknown. Calculation of the new multiallele test statistic and its P-value is very simple and utilizes the approximate distribution of R(2). We provide a computer program that evaluates approximate as well as "exact" permutational P-values.     

The Maintenance of Single-Locus Polymorphism. IV. Models with Mutation from Existing Alleles

The ability of viability selection to maintain allelic polymorphism is investigated using a constructionist approach. In extensions to the models we have previously proposed, a population is bombarded with a series of mutations whose fitnesses in conjunction with other alleles are functions of the corresponding fitnesses with a particular allele, the parent allele, already in the population. Allele frequencies are iterated simultaneously, thus allowing alleles to be driven to extinction by selection. Such models allow very high levels of polymorphism to evolve: up to 38 alleles in one case. Alleles that are lethal as homozygotes can evolve to surprisingly high frequencies. The joint evolution of allele frequencies and viabilities highlights the necessity to consider more than the current morphology of a population. Comparisons are made with the neutral theory of evolution and it is suggested that failure to reject neutrality using the Ewens-Watterson test cannot be regarded as evidence for the neutral theory.     

Maintenance of Genetic Variability under the Joint Effect of Mutation, Selection and Random Drift

Formulae are developed for the distribution of allele frequencies (the frequency spectrum), the mean number of alleles in a sample, and the mean and variance of heterozygosity under mutation pressure and under either genic or recessive selection. Numerical computations are carried out by using these formulae and Watterson's (1977) formula for the distribution of allele frequencies under overdominant selection. The following properties are observed: (1) The effect of selection on the distribution of allele frequencies is slight when 4Ns 相似文献   

Analysis of diversity of Russian and Ukrainian bread wheat (Triticum aestivum L.) cultivars for high-molecular-weight glutenin subunits

The allelic diversity of high-moleculat-weght glutenin subunits (H WIGS) in Russian and Ukrainian bread wheat cultivars was analyzed. The diversity of spring wheat cultivars for alleles of the Glu-1 loci is characterized by medium values of the polymorphism index (polymorphism information content, PlC), and in winter wheats it varies from high at the Glu-A1 locus to low at the Glu-D1 locus. The spring and winter cultivars differ significantly in the frequencies of alleles of the glutenin loci. The combination of the Glu-A1b, Glu-B1c, and Glu-D1a alleles prevails among the spring cultivars, and the combination of the Glu-A1a, Glu-B1c, and Glu-D1d alleles prevails among the winter cultivars. The distribution of the Glu-1 alleles significantly depends on the moisture and heat supply in the region of origin of the cultivars. Drought resistance is associated with the Glu-D1a allele in the spring wheat and with the Glu-B1b allele in the winter wheat. The sources of the Glu-1 alleles were identified in the spring and wheat cultivars. The analysis of independence of the distribution of the spring and winter cultivars by the market classes and by the alleles of the HMWGS loci showed a highly significant association of the alleles of three Glu-1 loci with the market classes in foreign cultivars and independence or a weak association in the Russian and Ukrainian cultivars. This seems to be due to the absence of a statistically substantiated system of classification of the domestic cultivars on the basis of their quality.     

Evaluation of ancestral inbreeding coefficients: Ballou’s formula versus gene dropping

Estimation of the purging of detrimental effects through inbreeding and selection is an important issue in conservation genetics opening new perspectives for the management of small populations. In 1997 Ballou proposed the ancestral inbreeding coefficient, which is calculated recursively via pedigree inbreeding coefficients, as a tool for evaluating the purging of deleterious alleles in zoo populations. The formula of Ballou assumes independence of inbreeding and ancestral inbreeding coefficients at any stage of the recursion. This study investigates the consequences of this inaccuracy on the estimation of true ancestral inbreeding, i.e. the proportion of alleles within a genome that has undergone inbreeding in the past. As an alternative we propose the estimation of ancestral inbreeding by the method of gene dropping. The methods are compared by stochastic simulation for various models with respect to mode of inheritance (neutral, detrimental and lethal alleles) and different settings for population size and initial allele frequencies. In all scenarios the proportion of alleles within a genome that has undergone inbreeding in the past was overestimated by Ballou’s formula. The overestimation was more pronounced in smaller populations but was not affected by genetic model or initial allele frequency. In contrast, the ancestral inbreeding coefficient calculated by gene dropping provided a robust estimate of ancestral inbreeding in most models and settings. A marginal overestimation was observed only in models with lethal alleles. Therefore, we recommend applying the gene dropping approach to estimate ancestral inbreeding coefficients.     

Heterosis or Neutrality?

Various statistics have been proposed on an ad hoc basis to test whether alleles at a locus are selectively neutral. By considering population models in which selection operates, this paper shows that the population homozygosity is a powerful test statistic for testing departures from neutrality, in the direction of heterozygote advantage or disadvantage. The sample homozygosity plays a similar role when only sample data are available. Some numerical examples are included, showing the application of the test.—An analysis is made of the effect of heterosis on such quantities as the expected number of alleles in the population or sample, the effective number of alleles, the expected homozygosity, and on the population and sample allele frequency distributions generally.     

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.     

Polymorphism for Dimerizing Ability at the Esterase-5 Locus in DROSOPHILA PSEUDOOBSCURA

Three new alleles are reported at the esterase-5 locus of Drosophila pseudoobscura. All three of these alleles are different from those previously reported in their ability to dimerize. One allele will not form heterodimers or homodimers and exists only as a monomer. A second allele does form heterodimers but will not form homodimers. The third allele forms both hetero- and homodimers as well as forming monomers. Estimates of the frequency of these alleles in a natural population are given. The existence of these is discussed with respect to recently proposed models for a molecular mechanism for heterosis.     

Parental origin of germ-line and somatic mutations in the retinoblastoma gene

Segregation analysis of polymorphic sites within the retinoblastoma (RB) gene and on chromosome 13, as well as the parental origin of the lost allele in the tumor, were analyzed in 24 families with RB patients. Four mutant alleles transmitted through the germ-line and seven de novo germ-line mutant alleles were identified in 11 patients with hereditary RB. Segregation analysis within the RB gene and on chromosome 13 was useful for DNA diagnosis of susceptibility to RB in relatives of hereditary patients, even if mutations were not identified. All seven de novo germ-line mutant alleles were paternally derived. The bias toward the paternal allele for de novo germ-line mutations of the RB gene was statistically significant. Seven paternal alleles and six maternal alleles were lost in 13 non-hereditary RB tumors with no bias in the parental origin of the somatic allele loss. These results suggest that the physical environment or a deficiency in DNA repair during spermatogenesis may be associated with significant risk factors for de novo germ-line mutations.     

On averaging power for genetic association and linkage studies

A power calculation is crucial in planning genetic studies. In genetic association studies, the power is often calculated using the expected number of individuals with each genotype calculated from an assumed allele frequency under Hardy-Weinberg equilibrium. Since the allele frequency is often unknown, the number of individuals with each genotype is random and so a power calculation assuming a known allele frequency may be incorrect. Ambrosius et al. recently showed that the power ignoring this randomness may lead to studies with insufficient power and proposed averaging the power due to the randomness. We extend the method of averaging power in two directions. First, for testing association in case-control studies, we use the Cochran-Armitage trend test and find that the time needed for calculating the averaged power is much reduced compared to the chi-square test with two degrees of freedom studied by Ambrosius et al. A real study is used for illustration of the method. Second, we extend the method to linkage analysis, where the number of identical-by-descent alleles shared by siblings is random. The distribution of identical-by-descent numbers depends on the underlying genetic model rather than the allele frequency. The robust test for linkage analysis is also examined using the averaged powers. We also recommend a sensitivity analysis when the true allele frequency or the number of identical-by-descent alleles is unknown.     

Spontaneous mutability in UV-sensitive excision-defective strains of Saccharomyces.

The genes RAD1, RAD2, RAD3 and RAD4 encode enzymes in the pathway leading to excision repair of UV-induced DNA damage in Saccharomyces cerevisiae. Four mutant alleles of these loci (rad1-1, rad2-2, rad3-12, and rad4-3) were studied for their effect on spontaneous reversion rate to lysine and histidine independence, by means of the 1000-compartment fluctuation test of von Borstel, Cain and Steinberg. Of these four excision-defective alleles, only rad3-12 was found to substantially increase the spontaneous reversion rate of the nonsense-suppressible lys1-1 allele, both through locus reversion as well as by forward mutation at one of eight suppressor loci. Similarly, only rad3-12 conferred a considerable increase in the reversion frequency of the missense his1-7 mutant. As the RAD3 gene product is believed to mediate the first step in the excision-repair pathway, it is assumed that spontaneous lesions in the rad3 strain are channelled into a mutagenic repair pathway, thus accounting for the enhanced spontaneous mutation rate.     

Methods for estimating gene frequencies and detecting selection in bacterial populations

Recent breakthroughs in molecular technology, most significantly the polymerase chain reaction (PCR) and in situ hybridization, have allowed the detection of genetic variation in bacterial communities without prior cultivation. These methods often produce data in the form of the presence or absence of alleles or genotypes, however, rather than counts of alleles. Using relative allele frequencies from presence-absence data as estimates of population allele frequencies tends to underestimate the frequencies of common alleles and overestimate those of rare ones, potentially biasing the results of a test of neutrality in favor of balancing selection. In this study, a maximum-likelihood estimator (MLE) of bacterial allele frequencies designed for use with presence-absence data is derived using an explicit stochastic model of the host infection (or bacterial sampling) process. The performance of the MLE is evaluated using computer simulation and a method is presented for evaluating the fit of estimated allele frequencies to the neutral infinite alleles model (IAM). The methods are applied to estimate allele frequencies at two outer surface protein loci (ospA and ospC) of the Lyme disease spirochete, Borrelia burgdorferi, infecting local populations of deer ticks (Ixodes scapularis) and to test the fit to a neutral IAM.     

Allelic Negative Complementation at the Abruptex Locus of DROSOPHILA MELANOGASTER

The mutations of the Abruptex locus in Drosophila melanogaster fall into three categories. There are recessive lethal alleles and viable alleles. The latter can be divided into suppressors and nonsuppressors of Notch mutations. The recessive lethals are lethal in heterozygous combination with Notch. As a rule the recessive lethals are lethal also in heterozygous combination with the viable alleles. Heterozygous combinations of certain viable alleles are also lethal. In such heterozygotes, one heteroallele is a suppressor of Notch and the other is a nonsuppressor. Other heterozygous combinations of viable alleles are viable and have an Abruptex phenotype. The insertion of the wild allele of the Abruptex locus as an extra dose (carried by a duplication) into the chromosomal complement of the fly fully restores the viability of the otherwise lethal heterozygotes if two viable alleles are involved. The extra wild allele also restores the viability of heterozygotes in which a lethal and a suppressor allele are present. If, however, a lethal and a nonsuppressor are involved, the wild allele only partly restores the viability, and the effect of the wild allele is weakest if two lethal alleles are involved. It seems likely that of the viable alleles the suppressors of Notch are hypermorphic and the nonsuppressors are hypomorphic. The lethal alleles share properties of both types, and are possibly antimorphic mutations. It is suggested that the locus is responsible for a single function which, however, consists of two components. The hypermorphic mutations are defects of the one component and the hypomorphic mutations of the other. In heterozygotes their cumulative action leads to decreased viability. The lethal alleles are supposed to be defects of the function as a whole. The function controlled by the locus might be a regulative function.     

A fast, unbiased and exact allelic test for case-control association studies

Association studies are traditionally performed in the case-control framework. As a first step in the analysis process, comparing allele frequencies using the Pearson's chi-square statistic is often invoked. However such an approach assumes the independence of alleles under the hypothesis of no association, which may not always be the case. Consequently this method introduces a bias that deviates the expected type I error-rate. In this article we first propose an unbiased and exact test as an alternative to the biased allelic test. Available data require to perform thousands of such tests so we focused on its fast execution. Since the biased allelic test is still widely used in the community, we illustrate its pitfalls in the context of genome-wide association studies and particularly in the case of low-level tests. Finally, we compare the unbiased and exact test with the Cochran-Armitage test for trend and show it perfoms similarly in terms of power. The fast, unbiased and exact allelic test code is available in R, C++ and Perl at: http://stat.genopole.cnrs.fr/software/fueatest.     

Marker selection for the transmission/disequilibrium test, in recently admixed populations.

Recent admixture between genetically differentiated populations can result in high levels of association between alleles at loci that are <=10 cM apart. The transmission/disequilibrium test (TDT) proposed by Spielman et al. (1993) can be a powerful test of linkage between disease and marker loci in the presence of association and therefore could be a useful test of linkage in admixed populations. The degree of association between alleles at two loci depends on the differences in allele frequencies, at the two loci, in the founding populations; therefore, the choice of marker is important. For a multiallelic marker, one strategy that may improve the power of the TDT is to group marker alleles within a locus, on the basis of information about the founding populations and the admixed population, thereby collapsing the marker into one with fewer alleles. We have examined the consequences of collapsing a microsatellite into a two-allele marker, when two founding populations are assumed for the admixed population, and have found that if there is random mating in the admixed population, then typically there is a collapsing for which the power of the TDT is greater than that for the original microsatellite marker. A method is presented for finding the optimal collapsing that has minimal dependence on the disease and that uses estimates either of marker allele frequencies in the two founding populations or of marker allele frequencies in the current, admixed population and in one of the founding populations. Furthermore, this optimal collapsing is not always the collapsing with the largest difference in allele frequencies in the founding populations. To demonstrate this strategy, we considered a recent data set, published previously, that provides frequency estimates for 30 microsatellites in 13 populations.     

Gene flow's effect on the genetic architecture of a local adaptation and its consequences for QTL analyses

This paper uses computer simulations to determine how gene flow between populations affects (1) the genetic architecture of a local adaptation and (2) properties of alleles segregating in quantitative trait locus (QTL) mapping populations. Results suggest that the average magnitude of an allele that causes a phenotypic difference between populations declines as the migration rate increases, but with an increase in migration, alleles of larger magnitude cause proportionally more of the phenotypic difference between populations. Gene flow between populations that are used in a QTL study tends to cause the average magnitude and percent variance explained (PVE) of an allele in a mapping population to increase. Thus, although the average magnitude of an allele causing a difference declines with migration the average magnitude or PVE of an allele in a QTL mapping population may increase. The reason is that the probability an allele is sampled for a QTL mapping population is in direct proportion to its frequency and alleles of larger magnitude tend to segregate at relatively higher frequencies than alleles of smaller effect with an increased migration. As the rate of gene flow between populations increases, the proportion of the phenotypic difference explained by alleles that are segregating in a QTL mapping population (and therefore potentially detected) decreases. Lastly, results suggest QTL alleles of large effect (>20% PVE) should be commonly found, provided the divergence time between populations is not too long or optima of populations are not too far apart.     

Genetic quality of the Miyaluo captive forest musk deer (Moschus berezovskii) population as assessed by microsatellite loci

The Miyaluo captive forest musk deer population (Sichuan Province, China) is one of the largest captive breeding populations in the world. In order to evaluate the genetic quality and provide available genetic management strategy, seven polymorphism microsatellite loci were applied to assess the genetic variation of the Miyaluo forest musk deer. The results indicated that a total of 168 alleles were detected from these seven microsatellite loci in 361 individuals, and the number of the alleles per locus ranged from 12 to 41 with a mean of 24. The average observed heterozygosity, expected heterozygosity, and PIC were 0.782, 0.854, and 0.837, respectively. Considering the results of the loci Hardy–Weinberg equilibrium test, the comparison of the common allele frequency as well as the private allele between the adults and juveniles, we concluded that the heterozygosity and the genetic diversity of the Miyaluo captive breeding population are increasing due to the input of new individuals from other populations. However, the frequency of some alleles declined sharply, and some were even lost indicating that there is a risk for diversity loss. Thus, we proposed an improved management and breeding strategy for the captive breeding population of the forest musk deer.