- Research article
- Open Access

# Selective Phenotyping, Entropy Reduction, and the Mastermind game

- Julien Gagneur
^{1}Email author, - Markus C Elze
^{2}and - Achim Tresch
^{3}

**12**:406

https://doi.org/10.1186/1471-2105-12-406

© Gagneur et al; licensee BioMed Central Ltd. 2011

**Received:**1 April 2011**Accepted:**20 October 2011**Published:**20 October 2011

## Abstract

### Background

With the advance of genome sequencing technologies, phenotyping, rather than genotyping, is becoming the most expensive task when mapping genetic traits. The need for efficient selective phenotyping strategies, *i*.*e*. methods to select a subset of genotyped individuals for phenotyping, therefore increases. Current methods have focused either on improving the detection of causative genetic variants or their precise genomic location separately.

### Results

Here we recognize selective phenotyping as a Bayesian model discrimination problem and introduce SPARE (Selective Phenotyping Approach by Reduction of Entropy). Unlike previous methods, SPARE can integrate the information of previously phenotyped individuals, thereby enabling an efficient incremental strategy. The effective performance of SPARE is demonstrated on simulated data as well as on an experimental yeast dataset.

### Conclusions

Using entropy reduction as an objective criterion gives a natural way to tackle both issues of detection and localization simultaneously and to integrate intermediate phenotypic data. We foresee entropy-based strategies as a fruitful research direction for selective phenotyping.

## Keywords

- Quantitative Trait Locus
- Recombinant Inbred Line
- Quantitative Trait Locus Mapping
- Quantitative Trait Locus Effect
- Quantitative Trait Locus Interval

## Background

Quantitative trait mapping aims at finding associations between genetic loci and phenotypes in a population of individuals. As genotyping costs are rapidly decreasing, often phenotyping becomes the limiting factor for mapping genetic traits. Recent techniques of genotyping have enabled access to complete or quasi-complete genotypes, from model organisms like yeast to human individuals [1–3]. Moreover, the number of genotyped individuals has largely increased with genome-wide association studies covering tens of thousands of samples [4]. Yet many phenotypes cannot be measured at this scale, including behavioral traits or susceptibility to drug treatments. The need is therefore increasing for selective phenotyping, *i.e*. the selection of a limited number of samples among a larger library of genotyped samples in order to optimally map a quantitative trait locus (QTL) [5].

Existing methods for selective phenotying methods can be classified into 2 categories (See [6] for a more complete review). The first category aims to maximize the resolution of the mapped loci, typically by maximizing the number of recombination events in the selected set (e.g., [7, 8]). The second category has focused on improving the detection power of QTL effects, by selecting a balanced distribution of alleles (e.g., [5], also discussed in [9]). However, no method has been proposed to date that combines both objectives [6].

Here we present SPARE (Selective Phenotyping Approach by Reduction of Entropy), a method that tackles both the issue of QTL effect detection power as well as the issue of locus resolution in a single framework. The problem of selective phenotyping is seen as an experimental design question in which the choice of samples must be made so that the outcome of their phenotyping allows, in expectation, optimal discrimination between the candidate genetic models. Here we follow a Bayesian mapping approach and the objective criterion is the minimization of the entropy of the a posteriori probability of the candidate models. We describe our method for single-causative locus models and fully genotyped populations with only two segregating genotypes, such as backcross, double-haploid, or recombinant inbred lines. An implementation of SPARE in the statistical programming language R [10] is available in Additional File 1 and at http://www.lmb.uni-muenchen.de/tresch/mastermind.html.

## Methods

### Mapping Strategy

In our setting, we have genotyped a set $\mathcal{J}=\{1,\mathrm{...},N\}$ of individuals (samples) at a set $\mathcal{L}=\{1,\mathrm{...},M\}$ of loci. The genotype of a sample $j\in \mathcal{J}$ at a locus $\lambda \in \mathcal{L}$ is denoted by ${x}_{j,\lambda}\in \mathcal{G}$, where $\mathcal{G}$ is the set of possible genotypes. The complete genotype information is thus stored in the matrix $x={({x}_{j,\lambda})}_{j\in \mathcal{J},\lambda \mathcal{L}\in ,}$, which is considered as a constant. We assume only two segregating genotypes, such as backcross, double-haploid or recombinant inbred lines and thus $\mathcal{G}=\{\pm 1\}$. Phenotyping a subset of samples $S\subset \mathcal{J}$ produces a phenotype vector ${y}^{S}\in {\mathcal{P}}^{s}$, where $\mathcal{P}$ denotes the set of possible phenotypes (typically, phenotypes are considered either binary or continuous). The genotype submatrix of *x* corresponding to the samples *S* is denoted by ${x}^{S}={\left({x}_{s,\lambda}\right)}_{s\in S,\lambda \in \mathcal{L}}$. The vector of phenotype measurements *y*^{
S
}is the realization of a random variable *Y*^{
S
}. If *S* consists of only one individual *s*, we shall write by abuse of notation *y*^{
S
}instead of *y*^{{s}}and *x*^{
s
}instead of *x*^{{s}}.

*y*

^{ s }of an individual with genotype

*x*

^{s}, given the QTL locus is λ,

*θ*is a suitable set of parameters for the respective model. Assuming observational independence given the genotype, the QTL and the parameters, the likelihood for the locus λ is obtained by summing over the parameters:

Specific parameterizations of (1) as well as the choice of appropriate priors Pr(*θ* | λ,*x*^{
S
}) in (2) will be discussed later in this section. Note that the individual observations in *y*^{
S
}, conditioned only on λ and *x*^{
S
}, are not independent in general.

The locus prior *π*(λ) can for instance be chosen to be uniform.

### Cost function

*H*(Λ) as a measure of the uncertainty about the causative QTL. After having phenotyped a sample subset

*S*, the remaining uncertainty about the QTL is:

*S*without knowledge of the phenotypic outcome

*y*

^{ S }, so the best we can do is to minimize the expectation of

*H*(Λ|

*Y*

^{ S }=

*y*

^{ S },

*x*

^{ S }) with respect to

*Y*

^{ S }|

*x*

^{ S }. Given a maximum number

*n*of individuals to be phenotyped, the task is to select a sample subset $S\subset \mathcal{J}$ of size |

*S*| =

*n*which minimizes

### Modeling of the genotype-phenotype relation

*P*(ϵ = ±1) =

*w*

_{±1}. The typical choices are either

*w*

_{1}= 1,

*w*

_{-1}= 0 for an a priori known genotype-phenotype association, as in typical linkage analysis studies, or ${w}_{1}={w}_{-1}=\frac{1}{2}$ for a completely unknown association. If the phenotype is fully inheritable and perfectly measured, the binary phenotypes

*Y*

^{ s }deterministically depends upon ${x}_{\lambda}^{s}$, given the set of parameters

*θ*which consists only of the variable ϵ:

*θ*= (

*q*

_{±1},ϵ). Equation (2) has to be evaluated for all loci λ. Therefore we make use of conjugate priors for

*θ*to obtain a closed form solution of the integral. The conjugate prior is a Beta distribution

*α*

_{±1},

*β*

_{±1}. This fully determines the prior if we furthermore assume parameter independence,

*i.e*.

*θ*= (

*μ*

_{±1}, ϵ). We choose conjugate Gaussian priors

The parameters *σ*^{2}, η_{±1}, and *ν*^{2} are considered given and fixed. Again this completely specifies the prior under the assumption of parameter independence, *Pr*(*θ*) = *Pr*(ϵ)*Pr*(*μ*_{1})*Pr*(*μ*_{-1}). For all three cases, closed form solutions of equation (2) are derived in the supplements (Additional File 2).

### Incremental selection

*S*for which the phenotypes are already known. Equations (4), (5) and (6) hold with

*S*∪

*T*in place of

*S*. The cost function is now defined given

*T*and

*y*

^{ T }:

By means of cost (14), |*S*| samples are selected and afterwards phenotyped. These selection/phenotyping rounds are iterated until the desired number of individuals is reached. In our applications, we have implemented the case where selection is incremental thus the integration above is done over a single dimension.

#### Remark

It was a discussion about an optimal strategy in the famous recreational Mastermind game that led us to formula (14) and initiated our research. Indeed if we choose loci, genotypes, phenotypes and genotype-phenotype modeling in an abusive way, identifying the hidden target sequence of colors is tantamount to bringing the cost function down to zero. We refer the reader to the supplements for a precise description of the algorithm. There, we implement a SPARE variant that plays the Mastermind game (and does it considerably better than most of us). The Mastermind game has been recently proposed as an educational tool to train biology students to experimental design [12]. Multiple algorithms have been proposed to play Mastermind ([13, 14] among many others), including entropy-based approaches [15].

### Simulations

Simulations for two scenarios were done: the *strong linkage* scenario and the *unbalanced allele* scenario. For both scenarios, linkage mapping for a backcross experiment was simulated with three biological replicates of each of the two parental F0 haploid strains and 100 F1 haploid trains. The 100 offspring genotypes were simulated by drawing a number of recombination events given the chromosome length in cM. The positions of the recombination events were drawn uniformly randomly along the chromosome. One unique QTL was drawn from one of the marker positions. The Gaussian model was applied. The fixed parameters were set to *σ*^{2} = 1,*η*_{±1} = 3, *ν*^{2} = 4, *w*_{±1} = 0.5. The Gaussian means *μ*_{±1} were drawn from the distributions (13), and the phenotypic values for the parental and the F1 strains were drawn according to (12). In the *strong linkage* scenario, we assumed a single chromosome genotyped at 100 markers equally spaced by 1 cM. No particular filtering is applied on the progeny so that at every marker, the two parental alleles are equally likely to occur. In the *unbalanced allele* scenario, we assumed a single chromosome genotyped at 100 markers equally spaced by 50 cM. Two alleles have low frequencies (20% in expectation) in the population: the causative allele and another randomly drawn one. This is obtained by first generating a much larger F1 population and then taking a random subset of 100 individuals weighting the probability to choose each individual according to its genotype for the two low frequency alleles. This procedure generates a population in which linkage is respected and the two low frequency alleles have in expectation a 20% frequency.

On these datasets, SPARE was initialized given the observation of the parental phenotypes. This setting matches the typical situation in which a cross is used to map a phenotypic difference observed between parental strains. SPARE, however, does not require initialization with parental phenotypes. We iterated the selection step until 30 out of the 100 F1 strains were included. For each selection size, we selected the samples by optimizing the MMA criterion using a two-step implementation [5] (a greedy search followed by greedy swaps). Optimization of the SBBL criterion was done in the same fashion. For each of the three methods, one obtains a series of 30 selection sets of size 1 to 30. The selection of SPARE consists of an increasing chain of sets, whereas the result of the other methods does not necessarily do so. We let MMA and SSBL run in a non-incremental fashion as these methods were not designed for this approach and, in contrast to SPARE, do not take the phenotypic information into account.

## Results

### Performance for detection and for resolution of QTLs

We compared our method to two methods representative of the two current selective phenotyping strategy families, *i.e*. one optimizing the resolution power and one the detection power. On the one side, the SSBL method [8] minimizes the Sum of the Squared Bin Lengths, where a bin is defined as a genomic interval within which there are no crossovers in any individual and that is bounded on either side by a crossover in at least one individual [16]. Minimizing the SSBL favors selections with evenly distributed crossovers along the genome. It therefore facilitates discriminating between markers in linkage and thus narrowing down QTL intervals. Methods focusing on resolution power are particularly advantageous when linkage is strong between alleles. On the other side, MMA, for Minimum Moment of Aberrations [5], maximizes the genotype dissimilarities at every marker in the selection thereby increasing the detection power. For two segregating genotypes, this favors balanced representation of the two alleles at each marker. This can be particularly advantageous for example in situation where the causative allele has a low frequency in the population, because it will enrich for carriers of the allele in the selection and facilitate the detection of its effect.

Simulations for two scenarios were done. In the *strong linkage* scenario, alleles at each marker are equally distributed but linkage is strong; therefore, resolution power is the main issue. In the *unbalanced allele* scenario, linkage is weak but the causative allele has a low frequency; thus, detection power is the main challenge (see Methods for details). Unbalanced allele frequencies are common in outbred populations. As in our simulation, unbalanced allele frequency can also occur in recombinant inbred line panels if combinations of alleles are lethal or lead to phenotypes incompatible with the assay of interest. The whole procedure was repeated 100 times for each scenario. In our formulation, performing QTL mapping amounts computing the posterior distribution over the candidate loci. We used the Kullback-Leibler distance to estimate how close the posterior distribution for a particular selection is to the posterior obtained when phenotypes of the whole population are taken into account. This measure allows us to compare the performance of different methods based on the selection they propose.

### Advantage of incremental phenotyping

### Application to a yeast dataset

We then applied SPARE to a linkage analysis dataset in yeast [17]. This dataset contains 184 meiotic recombinants genotyped at high resolution [1] and were phenotyped for growth in presence of arsenate (a binary trait). Using the phenotypic information of all 184 segregants, the original study identified and validated a single locus as determinant, the gene *PHO84*. SPARE, SSBL and MMA were run on this dataset, restricting for computational efficiency reasons to 314 markers regularly spaced at 100 cM distance. Different prior parameters were tried. The optimal prior, which takes as pseudo-counts of the Beta distributions the counts observed in the total population of 184 segregants: *q*_{1} ~ *Beta*(76, 17) and *q*_{-1} ~ *Beta*(12, 79), the uniform prior, and a strong prior but with an association to the locus closer to random (60%): *q*_{1} ~ *Beta*(60, 40) and *q*_{-1} ~ *Beta*(40, 60). The genotype-phenotype association is known here from the parental phenotypes and we therefore set *w*_{1} = 1 and *w*_{-1} = 0. As for the simulations, the Kullback-Leibler distance to the posterior estimated on the whole population was taken as a measure of the accuracy of the mapping for any selection.

## Discussion

We introduced SPARE, a selective phenotyping method based on Bayesian model discrimination that tackles both resolution and detection issues simultaneously. The objective function we minimize is an expected entropy which favors both a balanced selection of genotypes at all loci as well as evenly distributed crossovers across the genome. This formulation provides a natural way to deal with the issue of resolution: loci in linkage are seen as alternative competing models. Interestingly, linkage between neighboring loci does not have to be explicitly modeled. Indeed, two loci in linkage lead to similar likelihood. Moreover, our approach is able to make use of prior phenotype information and thereby proposes selections that enable the rapid reduction of QTL intervals. Simulations showed that SPARE can provide, on average, a more efficient selection strategy than previous methods. Application to a linkage analysis dataset in yeast showed that SPARE can propose, with appropriate priors, selections that enable a faster mapping of the QTL.

Even though SPARE is a Bayesian method, one should note that just like other selective phenotyping procedure, its application is independent of the method employed to perform QTL mapping, which can be frequentist. It is worth noting however that comparing QTL-mapping results such as LOD score curves for different selections is a not a well-defined task. Indeed, the objective of QTL mapping is not only to detect the causative locus, for example with a significant LOD score, but also to have a narrow confidence interval on its location. One contribution of the present study is to propose the entropy of the posterior to quantify how certain, *i.e*. how 'high and well-peaked', the mapping is. It relies on a Bayesian framework.

SPARE requires a prior distribution of the trait given the causative locus. The application to the yeast dataset showed that sensitivity to this prior distribution mattered. SPARE could only give the best results when applied with the prior distribution close to the optimal one. This is currently a weak point of the method.

Two extensions could broaden the applicability of the current framework: First, the restriction to the case of simple traits with a single causative locus may be dropped. One would like to detect complex traits, *i.e*. those involving multiple causative loci. Second, a selection step that chooses more than one sample at a time is desirable as phenotyping individuals in small batches might be more cost-effective than for one individual at a time. However, both improvements are computationally prohibitive if our cost function is used without modification, because then the task turns into a combinatorial optimization problem. Instead, efficient approximations must be developed to cope with the explosion in the number of hypotheses (Equation 5), and to solve the integration step in high dimensional space (Equation 8 and 14). The analogy to the Mastermind game suggests directions to a possible solution. For instance Cotta and colleagues [18] do not consider all possible hypotheses in the Mastermind game; they evaluate a representative subset of it, and this subset is evolved during the course of the game. Alternatively, the integration required in Equation 8 might be replaced by a sampling method. Another alternative is to change the cost function. It has been shown that Mastermind strategies that simply maximize the number of possible answers do surprisingly well [19]. An analogous phenotyping strategy would be promising even if more than one locus was associated with the phenotype. However, it was unclear to us how to transfer that approach to continuous phenotypes.

Our simulations also showed that random selections are actually performing reasonably well compared to standard methods in a linkage analysis approach, with balanced allele frequencies. While unexpected, this observation can be explained. Indeed, optimal selections are selections that are simultaneously balanced at each polymorphism and with evenly distributed break points along the genome. The natural process of recombination actually promotes this situation. Random draws thus typically show these desired features.

Recently, the use of Fisher's information to evaluate experimental designs in the context of QTL mapping, including selective phenotyping has been discussed [9]. In contrast to our method, Sen and colleagues assume the QTL locations given and use criteria based on Fisher's information to compare experimental designs. Hence, while the method of Sen and collaborators allows dealing with multiple-loci models, it does not address the question of resolution. The relationships between Fisher's information and Shannon entropy criteria in this context constitute an interesting direction of research.

## Conclusion

SPARE is a Bayesian model discrimination approach to perform selective phenotyping. Using entropy reduction as an objective criterion enabled tackling both issues of detection and localization jointly for the first time and integrating intermediate phenotypic data. Altogether, entropy-based strategies appear as a promising family of approaches for the question of selective phenotyping.

## Declarations

### Acknowledgements

The authors would like to thank the lab of Lars Steinmetz for financial support and Raeka Aiyar for careful reading of the manuscript. The authors are grateful to a referee for drawing their attention to the publication of Cotta et al. AT was supported by the Deutsche Forschungsgemeinschaft, SFB646, and by an LMU_{excellent} guest professorship to Patrick Cramer.

## Authors’ Affiliations

## References

- Mancera E, Bourgon R, Brozzi A, Huber W: High-resolution mapping of meiotic crossovers and non-crossovers in yeast. Nature. 2008, 454 (7203): 479-485. 10.1038/nature07135.PubMed CentralView ArticlePubMedGoogle Scholar
- Levy S, Sutton G, Ng PC, Feuk L, Halpern AL, Walenz BP, Axelrod N, Huang J, Kirkness EF, Denisov G, Lin Y, Macdonald JR, Pang AWC, Shago M, Stockwell TB, Tsiamouri A, Bafna V, Bansal V, Kravitz SA, Busam DA, Beeson KY, Mcintosh TC, Remington KA, Abril JF, Gill J, Borman J, Rogers YH, Frazier ME, Scherer SW, Strausberg RL, Venter JC: The Diploid Genome Sequence of an Individual Human. Plos Biol. 2007, 5 (10): e254-10.1371/journal.pbio.0050254.PubMed CentralView ArticlePubMedGoogle Scholar
- Wheeler DA, Srinivasan M, Egholm M, Shen Y, Chen L, Mcguire A, He W, Chen YJ, Makhijani V, Roth GT, Gomes X, Tartaro K, Niazi F, Turcotte CL, Irzyk GP, Lupski JR, Chinault C, Song XZ, Liu Y, Yuan Y, Nazareth L, Qin X, Muzny DM, Margulies M, Weinstock GM, Gibbs RA, Rothberg JM: The complete genome of an individual by massively parallel DNA sequencing. Nature. 2008, 452 (7189): 872-10.1038/nature06884.View ArticlePubMedGoogle Scholar
- Altshuler D, Daly MJ, Lander ES: Genetic mapping in human disease. Science. 2008, 322 (5903): 881-8. 10.1126/science.1156409.PubMed CentralView ArticlePubMedGoogle Scholar
- Jin C, Lan H, Attie AD, Churchill GA, Bulutuglo D, Yandell BS: Selective Phenotyping for Increased Efficiency in Genetic Mapping Studies. Genetics. 2004, 168 (4): 2285-2293. 10.1534/genetics.104.027524.PubMed CentralView ArticlePubMedGoogle Scholar
- Rosa GJM, Leon ND, Rosa AJM: Review of microarray experimental design strategies for genetical genomics studies. Physiological Genomics. 2006, 28: 15-23. 10.1152/physiolgenomics.00106.2006.View ArticlePubMedGoogle Scholar
- Jannink JL: Selective Phenotyping to Accurately Map Quantitative Trait Loci. Crop Science. 2005, 45 (3): 901-908. 10.2135/cropsci2004.0278.View ArticleGoogle Scholar
- Xu Z, Zou F, Vision TJ: Improving Quantitative Trait Loci Mapping Resolution in Experimental Crosses by the Use of Genotypically Selected Samples. Genetics. 2005, 170: 401-408. 10.1534/genetics.104.033746.PubMed CentralView ArticlePubMedGoogle Scholar
- Sen S, Johannes F, Broman K: Selective Genotyping and Phenotyping Strategies in a Complex Trait Context. Genetics. 2009, 181 (4): 1613-10.1534/genetics.108.094607.PubMed CentralView ArticlePubMedGoogle Scholar
- R Development Core Team: R: A Language and Environment for Statistical Computing. 2011, R Foundation for Statistical Computing, Vienna, Austria, [ISBN 3-900051-07-0], [http://www.R-project.org]Google Scholar
- Box GEP, Hill WJ: Discrimination among Mechanistic Models. Technometrics. 1967, 9: 57-71. 10.2307/1266318.View ArticleGoogle Scholar
- Strom A, Barolo S: Using the Game of Mastermind to Teach, Practice, and Discuss Scientific Reasoning Skills. PLoS Biology. 2011, 9: e1000578-10.1371/journal.pbio.1000578.PubMed CentralView ArticlePubMedGoogle Scholar
- Knuth D: The Computer as a Master Mind. Journal of Recreational Mathematics. 1976, 9: 1-6.Google Scholar
- Koyama Lai: An Optimal Mastermind Strategy. Journal of Recreational Mathematics. 1993, 25: 251-256.Google Scholar
- Neuwirth E: Some strategies for Mastermind. Mathematical Methods of Operations Research. 1982, 26: 257-278.View ArticleGoogle Scholar
- Vision T, Brown D, Shmoys D, Durrett R, Tanksley S: Selective Mapping: A Strategy for Optimizing the Construction of High-Density Linkage Maps. Genetics. 2000, 155: 407-PubMed CentralPubMedGoogle Scholar
- Gagneur J, Sinha H, Perocchi F, Bourgon R, Huber W, Steinmetz LM: Genome-wide allele- and strand-specific expression profiling. Mol Syst Biol. 2009, 5: 274-PubMed CentralView ArticlePubMedGoogle Scholar
- Cotta C, Guervos JJM, Garcia AM, Runarsson TP: Entropy-Driven Evolutionary Approaches to the Mastermind Problem. PPSN (2)'10. 2010, Springer, 421-431.Google Scholar
- Kooi B: Yet another mastermind strategy. ICGA Journal. 2005, 28: 13-20.Google Scholar

## Copyright

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.