  Assessing the statistical power to detect susceptibility variants plays a
critical role in GWA studies both from the prospective and retrospective points
of view. Power is empirically estimated by simulating phenotypes under a
disease model H1. For this purpose, the "gold" standard consists in simulating
genotypes given the phenotypes (e.g. Hapgen). We introduce here an alternative
approach for simulating phenotypes under H1 that does not require generating
new genotypes for each simulation. In order to simulate phenotypes with a fixed
total number of cases and under a given disease model, we suggest three
algorithms: i) a simple rejection algorithm; ii) a numerical Markov Chain
Monte-Carlo (MCMC) approach; iii) and an exact and efficient backward sampling
algorithm. In our study, we validated the three algorithms both on a
toy-dataset and by comparing them with Hapgen on a more realistic dataset. As
an application, we then conducted a simulation study on a 1000 Genomes Project
dataset consisting of 629 individuals (314 cases) and 8,048 SNPs from
Chromosome X. We arbitrarily defined an additive disease model with two
susceptibility SNPs and an epistatic effect. The three algorithms are
consistent, but backward sampling is dramatically faster than the other two.
Our approach also gives consistent results with Hapgen. Using our application
data, we showed that our limited design requires a biological a priori to limit
the investigated region. We also proved that epistatic effects can play a
significant role even when simple marker statistics (e.g. trend) are used. We
finally showed that the overall performance of a GWA study strongly depends on
the prevalence of the disease: the larger the prevalence, the better the power.
