A straightforward multiallelic significance test for the Hardy-Weinberg equilibrium law

GENOMICS AND BIOINFORMATICS

RESEARCH ARTICLE

A straightforward multiallelic significance test for the Hardy-Weinberg equilibrium law

Marcelo S. Lauretto^I; Fabio Nakano^II; Silvio R. Faria Jr^III; Carlos A.B. Pereira^III; Julio M. Stern^III

^IEscola de Artes, Ciências e Humanidades, Universidade de São Paulo, São Paulo, Brazil

^IIInstituto de Química, Universidade de São Paulo, São Paulo, Brazil

^IIIInstituto de Matemática e Estatística Universidade de São Paulo, São Paulo, Brazil

^{Send correspondence to}

ABSTRACT

Much forensic inference based upon DNA evidence is made assuming Hardy-Weinberg Equilibrium (HWE) for the genetic loci being used. Several statistical tests to detect and measure deviation from HWE have been devised, and their limitations become more obvious when testing for deviation within multiallelic DNA loci. The most popular methods-Chi-square and Likelihood-ratio tests-are based on asymptotic results and cannot guarantee a good performance in the presence of low frequency genotypes. Since the parameter space dimension increases at a quadratic rate on the number of alleles, some authors suggest applying sequential methods, where the multiallelic case is reformulated as a sequence of "biallelic" tests. However, in this approach it is not obvious how to assess the general evidence of the original hypothesis; nor is it clear how to establish the significance level for its acceptance/rejection. In this work, we introduce a straightforward method for the multiallelic HWE test, which overcomes the aforementioned issues of sequential methods. The core theory for the proposed method is given by the Full Bayesian Significance Test (FBST), an intuitive Bayesian approach which does not assign positive probabilities to zero measure sets when testing sharp hypotheses. We compare FBST performance to Chi-square, Likelihood-ratio and Markov chain tests, in three numerical experiments. The results suggest that FBST is a robust and high performance method for the HWE test, even in the presence of several alleles and small sample sizes.

Key words: Hardy-Weinberg equilibrium, significance tests, FBST.

Introduction

The Hardy-Weinberg law is one of most important principles in population genetics, and establishes a direct relationship between allele and genotypic proportions in a population. This law states that in a large population of panmictic dioecious organisms with non-overlapping generations, the allelic and genotypic frequencies at a locus will stay unchanged, provided that migration, mutation, and natural selection do not affect that locus. When these conditions hold, it is said that the locus is under Hardy-Weinberg Equilibrium (HWE).

This principle was discussed independently by Yule, Pearson and Castle (between 1902 and 1904), for some particular allele frequencies (see references in Crow and Kimura, 1972). In 1908, Godfrey Hardy presented the general principle for two alleles (Hardy, 1908). This principle was called Hardy's law for 35 years, until Stern (1943) called attention to an article of Weinberg (1908) showing the same principle at the same time and demonstrating its validity for multiple alleles (Crow, 1999).

Since its postulation, several results in population genetics and much forensic inference based upon DNA evidence have been based on the assumption that HWE is valid for the genetic loci of interest. Some statistical tests to detect and measure deviation from HWE have been devised, and their limitations have become more obvious when testing for deviation within multiallelic DNA loci. The most common approach consists of goodness-of-fits tests, like Chi-square and Likelihood-ratio, which are heavily based on asymptotic results, and can sometimes lead to false rejection or acceptance of HWE when the sample sizes are small and/or some genotype sample frequencies are very small (Emigh, 1980). Another approach involves exact tests, but is restricted to small dimensions and allele numbers.

A Bayesian sequential method for multiallelic HWE test was proposed by Pereira et al. (2006), who suggested reformulating the multiallelic case as a sequence of "biallelic" tests. In that work, the central component is the Full Bayesian Significance Test (FBST), an intuitive Bayesian approach which does not assign positive probabilities to zero measure sets when testing sharp hypotheses (Pereira and Stern, 1999). Although the sequential method avoids the quadratic increase of parameter space dimension with respect to the number of alleles, it is not obvious how to assess the general evidence of the original hypothesis; nor is it clear how to establish the significance level for its acceptance or rejection (see DeGroot, 1970).

In this work, we propose a method for the multiallelic HWE test, based on the FBST. FBST has many theoretic and practical advantages over other approaches, and it has shown to be robust in several high-dimensional problems (Lauretto et al., 2008).

Background

In this section we introduce some notations, and the Hardy-Weinberg Equilibrium (HWE) formula. Let us consider k alleles A₁, A₂, ..., A_k in a locus. The main interest is to assess the population relative frequencies of the genotypes A_iA_j (i, j = 1, 2, ..., k) which we denote by p_ij. As usual in the literature (see Hardy, 1908), we assume that the allele frequencies do not depend on sex and thus are symmetric, that is, A_iA_j is equivalent to A_jA_i and p_ij = p_ji. Therefore, the parameter of interest is the (lower triangular) matrix of genotype proportions:

θ = (θ_ij), with θ_ii = p_ii, θ_ij = 2p_ij for 1 < j < i < k.

We denote by p₁, p₂, ..., p_k the (unknown) population frequencies of alleles A₁, A₂, ..., A_k, with p_i> 0 and p_i = 1. When the locus is under HWE, the genotype proportions are as follows:

In order to test the HWE in a locus, one considers a sample of n individuals drawn randomly from the population. Such a sample can be presented as the array x, whose elements 1 < j < i < k, are counts of genotypes A_iA_j. The sample size n is n = x_ij, and the sample frequency of allele A_i is n_i = x_ij + x_ji. Note that = 2n, since all loci have two alleles and this sum is the number of alleles in the whole sample. The sample proportion of allele A_i, _i, is given by

Assuming that each individual genotype does not depend on remaining individuals in the same generation, we can consider that the genotype frequencies x_ij follow a multinomial distribution with unknown parameter θ,

Testing Procedures

In this section we present three tests used in our comparative study. These and other approaches are described by Emigh (1980), Guo and Thompson (1992), Hernández and Weir (1989) and Montoya-Delgado et al. (2001).

Chi-square goodness-of-fit test

This test involves calculating the sample chi square value,

with

Under HWE, this quantity has a chi-square distribution with k(k – 1)/2 degrees of freedom.

The Chi-square goodness-of-fit test with continuity correction involves calculating the previous statistics, with the subtraction in each term of a correction constant c:

Usually c = 0.5 is the value chosen.

Likelihood-ratio tests

The likelihood function, given a sample, follows directly from the multinomial distribution presented in Eq. (3). A Likelihood-ratio test is constructed by comparing the likelihood maximized under the hypothesis, L₀, with the maximum likelihood, L₁, not constrained by the hypothesis. For HWE we have

with the sample allelic frequencies, , given by Eq. (2).

The test statistic

is asymptotically distributed as a chi-square distribution with k(k – 1)/2 degrees of freedom.

Markov Chain Monte Carlo (MCMC) method

Proposed by Guo and Thompson (1992), the method consists of an adaptation of the Metropolis algorithm, with the construction of a Markov chain with equilibrium distribution matching the genotype probabilities under HWE of samples that have the same allelic counts as the observed data.

Under HWE and conditional on sample allele counts, n₁, n₂, ..., n_k, the probability of obtaining the sample x is (see Levene, 1949):

Given the data x, the test evaluates

where = {y:Pr(y) < Pr(x), x ∈ Γ₀} and

The MCMC algorithm is performed in order to estimate the probability P in E. (9). Rejection or acceptance of the null hypothesis depends on whether P is smaller than a pre-specified significance level α.

Methods

The Full Bayesian Significance Test (FBST)

The Full Bayesian Significance Test (FBST) was proposed by Pereira and Stern (1999) as a coherent and intuitive Bayesian test. It assumes that the hypothesis, H, is defined as a subset defined by inequality and equality constraints:

For simplicity, we often use H for Θ_H. FBST is particularly focused on precise hypotheses: i.e., dim(Θ_H) < dim(Θ). In this work, f_x(θ) denotes the posterior probability density function, given the observation x. Bold 0 and 1 denote vectors of appropriate dimensions.

For the HWE test, the parametric space consists of all arrays of genotype proportions

The space on hypothesis is

As previously stated, we consider that the genotype frequencies x follow a multinomial distribution, given by Eq. (3). Taking as a priori the Dirichlet distribution with parameters (1,1...1), i.e., a uniform distribution, then the a posteriori is a Dirichlet distribution with parameters (x₁₁ + 1, x₂₁ + 1, x₂₂ + 1, ..., x_kk + 1) which is proportional to the likelihood function (DeGroot, 1970):

The computation of the evidence measure used on the FBST is performed in two steps:

1. The optimization step consists of finding the maximum (supremum) of the posterior under the null hypothesis, θ* = arg sup_hf_x (θ), f^* = f_x )θ^*)

2. The integration step consists of integrating the posterior density over the Tangential Set, , where the posterior is higher than anywhere in the hypothesis, i.e.,

(H) is the evidence against H, and EV(H) = 1 – (H) is the evidence supporting (or in favor of) H. For a better understanding of this evidence measure, Figure 1 illustrates two examples in the biallelic case, showing the null and tangential sets (Θ_H and ). Since θ₂₁ = 1 – θ₁₁ – θ₂₂, the parametric space is fully defined by homozygote proportions, θ₁₁ and θ₂₂. The parameter space corresponds to the area inside the triangle. Sample genotype counts for A₁₁, A₂₁, A₂₂ and (H) are also shown in each graph. Marker '*' represents the point of maximum a posteriori density in the constrained space Θ_H, and the level curve tangent to θ^* corresponds to frontier. Intuitively, if the hypothesis set is in a region of "low" posterior density (as in the example 1), then is "heavy" and therefore (H) is "large" (≅ 0.91), meaning "strong" evidence against H. On the other hand, as illustrated by the example 2, if hypothesis set is in a region of "high" posterior density, then is a "small" set, and hence (H) is "small" (≅ 0.36), meaning "weak" evidence against H.

For HWE test, the point θ^* = arg sup_Hf_x(θ) is given as follows.

Rewriting the posterior pdf under HWE, we have:

Taking its logarithm,

where n_i = x_ij + x_ji, and h is the sum of sample heterozygote frequencies, h = .

By the constraint p_i = 1, we have:

The gradients of l_x(θ) are given by:

Hence, the optimal point under HWE is given by the vector p^* = (, ..., ) which satisfies:

Summing over all constraints and after some algebra, we obtain the following solution:

The computation of θ^* from p^*follows from Eq. (1).

The integration step may be performed by generating a set of M points {θ⁽¹⁾, θ⁽²⁾, ..., θ^(M) with a Dirichlet distribution with parameter (x₁₁ + 1, x₂₁ + 1, x₂₂ + 1, ..., x_kk + 1) and computing the percentage of points with posterior density greater than f ^*:

where I(stat) = 1 if stat is true and 0 otherwise. A more precise and efficient Monte Carlo method for the integration step is presented by Lauretto et al. (2003).

As with any significance test, this procedure requires the choice of a threshold level, τ, for acceptance/rejection of the hypothesis at a significance level α. Several alternative methods have been developed for establishing this threshold:

• An empirical power analysis, developed by Stern and Zacks (2002) and Lauretto et al. (2003), provides critical levels that are consistent and also effective for small samples.

• A threshold based on reference sensitivity analysis and paraconsistent logic is given by Stern (2004).

• Pereira et al. (2008) relates the e-value threshold to standard p-value thresholds.

• Madruga et al. (2001) proves the existence of a loss function that renders the FBST a true Bayesian decision-theoretic procedure.

• An asymptotically consistent threshold for a given confidence level was given by Stern (2007), and Borges and Stern (2007), which we adopt in this work.

Let us consider the cumulative distribution of the evidence value against the hypothesis, (τ), Pr((H) < τ) given θ⁰, the true value of the parameter. Under appropriate regularity conditions, for increasing sample size, we can state the following:

If H is false, θ⁰∉ H, then (H) converges (in probability) to one, that is, (τ) → δ(1).

If H is true, θ⁰∈ H, then (τ), the confidence level, is approximated by the function Q(t, h, τ) = Chi2(t - h Chi2^-1 (t, τ)), where t = dim(Θ), h = dim(H) and Chi2(df, x) denotes the cumulative chi-square distribution with df degrees of freedom.

Hence, to reject H with a significance level α, we can set τ = Q^-1 (t, h, 1-α), i.e., set τ such that Q(t,h, τ) = 1 - α.

Results and Discussion

The numerical experiments used in the performance analysis are based on three typical datasets. Two examples consist of simulated data used in the literature as benchmarks for comparing the performance of competing methods, while the third example is from a real dataset. These examples are presented in Figure 2, as lower triangular matrices containing genotype frequencies. The first example is taken from Louis and Dempster (1987) and consists of a sample of size 45 of genotype frequencies distributed in four alleles. The second example is given by Guo and Thompson (1992) and consists of a sample of size 30 of simulated genotype frequencies simulated under HWE with underlying allele frequencies (0.2, 0.2, 0.2, 0.2, 0.05, 0.05, 0.05, 0.05). The third example is from a rheumatoid arthritis (RA) study performed by Wordsworth et al. (1992), where two hundred and thirty RA patients were genotyped for the HLA-DR locus. The DR4 allele was subdivided into Dw4, Dw14 and other subtypes. DRX represents all non-DR1, non-Dw4, non-Dw14 alleles. This example is used by Chen and Thomson (1999) as a benchmark.

Our main interest is to compare the error convergence of FBST and other methods presented in this work (MCMC, Chi-square and Likelihood-ratio) for increasing sample sizes. For each sample size n ∈ {30, 50, 100, 200}, we simulated two collections of 100 samples. The first collection consists of samples drawn under HWE, i.e., each sample is drawn with a multinomial distribution with parameters (n, θ^*), with θ^* = arg max_θ∈Hf_x (θ). The second collection consists of samples drawn with a multinomial distribution with parameters (n, θ^(h)), where θ^(h) is drawn under the posterior distribution. That is, each sampling iteration is performed in two steps:

a) draw θ^(h) ~ Dirichlet (x₁₁ + 1, x₂₁ + 1, ..., x_kk + 1), where x_ij (1 < j < i < k) are the frequencies in the original dataset; and

b) draw X^(h) ~Multinomial (n, θ^(h)).

The Type I error (rejection rate of a true hypothesis) is estimated by the proportion of samples in collection 1 such that HWE is rejected, and the Type II error (acceptance rate of a false hypothesis) is estimated by the proportion of samples in collection 2 such that HWE is accepted. The performance criterion used in this work is the average error, i.e., the average of Types I and II error rates. Two standard significance levels, α ∈ {0.01, 0.05}, were used to calibrate the asymptotic acceptance/rejection threshold of each method.

A variability measure for the errors was obtained by performing 10 batches of simulations and computing the mean and standard deviation of average errors across the batches.

Figure 3 presents the average errors for FBST, MCMC, Chi-square (with continuity correction) and Likelihood ratio for simulated data based on examples 1, 2 and 3. The bar height represents the mean of average errors, and the vertical line on the top of each bar is the error bar, representing the mean ± one standard deviation of average errors.

For simulated data based on examples 1 and 3, the best competitors are FBST and Likelihood-ratio test, while for simulated data based on example 2, the best competitors are FBST and MCMC. In every case, we notice that FBST is always the best competitor (especially for small sample sizes, n < 100) or is very close to it.

These numerical results suggest that FBST is more stable than the competitors discussed in this paper, in the sense that it has good comparative performance for different datasets and allele numbers.

Final Remarks

We have introduced a simple and straightforward procedure for testing deviance from Hardy-Weinberg Equilibrium (HWE) in the presence of several alleles. This procedure was implemented in C language, and integrated into a system for parentage testing developed with FAPESP support, where it is applied in the selection of loci to be used for parentage testing. Further details of this project can be found at http://watson.fapesp.br/PIPEM/Pipe13/genet1. htm. Currently, the routine is available by request to the corresponding author.

Acknowledgments

The authors are grateful for the support of EACH-USP, IQ-USP, IME-USP, Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES), Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) and Fundação de Apoio à Pesquisa do Estado de São Paulo (FAPESP).

Received: November 4, 2008; Accepted: May 25, 2009.

Guest Editor: José Carlos Merino Mombach

• Borges W and Stern JM (2007) The rules of logic composition for the Bayesian epistemic e-values. Logic J IGPL 15:401-420.
• Chen JJ and Thomson G (1999) The variance for the disequilibrium coefficient in the individual Hardy-Weinberg test. Biometrics 55:1269-1272.
• Crow JF and Kimura M (1972) An Introduction to Population Genetics Theory. Harper & Row, New York, 591 pp.
• Crow JF (1999) Hardy, Weinberg and language impediments. Genetics 152:821-825.
• DeGroot MH (1970) Optimal Statistical Decisions. Wiley, New York, 489 pp.
• Emigh TH (1980) A comparison of tests for Hardy-Weinberg equilibrium. Biometrics 36:627-642.
• Guo SW and Thompson EA (1992) Performing the exact test of Hardy-Weinberg proportion for multiple alleles. Biometrics 48:361-372.
• Hardy GH (1908) Mendelian proportions in a mixed population. Science 28:49-50.
• Hernández JL and Weir BS (1989) A disequilibrium coefficient approach to Hardy-Weinberg testing. Biometrics 45:53-70.
• Lauretto MS, Pereira CAB, Stern JM and Zacks S (2003) Full Bayesian significance test applied to multivariate normal structure models. Braz J Probabil Stat 7:147-168.
• Lauretto MS, Pereira CAB and Stern JM (2008) The full Bayesian significance test for mixture models: Results in gene expression clustering. Genet Mol Res 7:883-897.
• Levene H (1949) On a matching problem arising in genetics. Ann Math Stat 20:91-94.
• Louis EJ and Dempster ER (1987) An exact test for Hardy-Weinberg and multiple alleles. Biometrics 43:805-811.
• Madruga M, Esteves LG and Wechsler S (2001) On the Bayesianity of Pereira-Stern tests. Test 10:291-299.
• Montoya-Delgado LE, Irony TZ, Pereira CAB and Whittle MR (2001) An unconditional exact test for the Hardy-Weinberg equilibrium law: Sample-space ordering using the Bayes factor. Genetics 158:875-883.
• Pereira CAB, Nakano F, Stern JM and Whittle MR (2006) Genuine Bayesian multiallelic significance test for the Hardy-Weinberg equilibrium law. Genet Mol Res 5:619-631.
• Pereira CAB and Stern JM (1999) Evidence and credibility: Full Bayesian significance test for precise hypotheses. Entropy J 1:69-80.
• Pereira CAB, Stern JM and Wechsler S (2008) Can a significance test be genuinely Bayesian? Bayesian Anal 3:79-100.
• Stern C (1943) The Hardy-Weinberg law. Science 97:137-138.
• Stern JM and Zacks S (2002) Testing the independence of Poisson variates under the Holgate bivariate distribution. The power of a new evidence test. Stat Probabil Lett 60:313-320.
• Stern JM (2004) Inconsistency analysis for statistical tests of hypothesis. In: López-Díaz M and Gil MA (eds) Soft Methodology and Random Information Systems. Springer, New York, pp 567-574.
• Stern JM (2007) Cognitive constructivism, eigen-solutions and sharp statistical hypotheses. Cybernet Hum Knowing 14:9-36.
• Weinberg W (1908) Über den Nachweis der Vererbung beim Menschen. Jahresh Wuertt Verh Vaterl Naturkd 64:369-382.
• Wordsworth P, Pile KD, Buckley JD, Lanchbury JSS, Ollier B, Lathrop M and Bell JI (1992) HLA heterozygosity contributes to susceptibility to rheumatoid arthritis. Am J Hum Genet 51:585-591.

Publication Dates

• Publication in this collection
  10 Aug 2009
• Date of issue
  2009

History

• Accepted
  25 May 2009
• Received
  04 Nov 2008