# QTL fine mapping with Bayesian $C(\pi)$: a simulation study

Irene van den Berg van den Berg, Sébastien S. Fritz, Didier Boichard

<figure><img src="image_1.png" /></figure>

## To cite this version:

Irene van den Berg van den Berg, Sébastien S. Fritz, Didier Boichard. QTL fine mapping with Bayesian $C(\pi)$: a simulation study. Genetics Selection Evolution, 2013, 45 (june), Non paginé. 10.1186/1297-9686-45-19. [hal-01000784](https://hal.science/hal-01000784)

HAL Id: hal-01000784

https://hal.science/hal-01000784

Submitted on 29 May 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# QTL fine mapping with Bayes C(π): a simulation study

Irene van den Berg<fnref n="1" />,<fnref n="2" />,³*, Sébastien Fritz⁴ and Didier Boichard<fnref n="1" />,<fnref n="2" />

**Abstract**

**Background:** Accurate QTL mapping is a prerequisite in the search for causative mutations. Bayesian genomic selection models that analyse many markers simultaneously should provide more accurate QTL detection results than single-marker models. Our objectives were to (a) evaluate by simulation the influence of heritability, number of QTL and number of records on the accuracy of QTL mapping with Bayes Cπ and Bayes C; (b) estimate the QTL status (homozygous vs. heterozygous) of the individuals analysed. This study focussed on the ten largest detected QTL, assuming they are candidates for further characterization.

**Methods:** Our simulations were based on a true dairy cattle population genotyped for 38 277 phased markers. Some of these markers were considered biallelic QTL and used to generate corresponding phenotypes. Different numbers of records (4387 and 1500), heritability values (0.1, 0.4 and 0.7) and numbers of QTL (10, 100 and 1000) were studied. QTL detection was based on the posterior inclusion probability for individual markers, or on the sum of the posterior inclusion probabilities for consecutive markers, estimated using Bayes C or Bayes Cπ. The QTL status of the individuals was derived from the contrast between the sums of the SNP allelic effects of their chromosomal segments.

**Results:** The proportion of markers with null effect (π) frequently did not reach convergence, leading to poor results for Bayes Cπ in QTL detection. Fixing π led to better results. Detection of the largest QTL was most accurate for medium to high heritability, for low to moderate numbers of QTL, and with a large number of records. The QTL status was accurately inferred when the distribution of the contrast between chromosomal segment effects was bimodal.

**Conclusions:** QTL detection is feasible with Bayes C. For QTL detection, it is recommended to use a large dataset and to focus on highly heritable traits and on the largest QTL. QTL statuses were inferred based on the distribution of the contrast between chromosomal segment effects.

**Background**

The first step in the identification of causative mutations underlying quantitative traits is to search for genomic regions associated with these traits, called quantitative trait loci (QTL). With the availability of genetic markers, a large number of QTL has been reported for various species and traits [1-3]. However, taking the final step from QTL to causative mutation has proven more difficult. In farm animals, most economically important

traits are complex and involve a large number of genes, each with a small effect. Furthermore, due to long range linkage disequilibrium, as is the case in dairy cattle [4], several markers can be in linkage disequilibrium with the causative mutation rather than only the markers closest to the mutation. As a consequence, reported QTL often have large confidence intervals [5], making the search for causative mutations difficult, and mutations underlying such traits have been identified only in a few cases [6]. Accurate and precise QTL mapping is thus a prerequisite to the search for causative mutations.

A large variety of models for QTL detection exists, including single-marker models, interval mapping [7] and Bayesian models [8]. Since quantitative traits are likely to be influenced by a large number of QTL, models that

<footnote>* Correspondence: irene.vanderberg@jouy.inra.fr</footnote>

<footnote n="1">¹INRA, UMR1313 Génétique animale et biologie intégrative, Domaine de Vilvert, 78350 Jouy-en-Josas, France</footnote>

<footnote n="2">²AGROPARISTECH, UMR1313 Génétique animale et biologie intégrative, 16 rue Claude Bernard, 75231 Paris 05, France</footnote>

<footnote>Full list of author information is available at the end of the article</footnote>
analyse all markers simultaneously should give more ac-
curate results than models that analyse only one or a
few markers at a time. With the Bayesian models devel-
oped for genomic selection, a large number of markers
can be analysed simultaneously [9]. In a simulation study
by Sahana et al. [10], Bayesian QTL mapping resulted in
higher power and more precise QTL locations than
single-marker, haplotype-based or mixed-model ap-
proaches. When dense single nucleotide polymorphism
(SNP) chips are used, loci influencing a trait are assumed
to be in linkage disequilibrium with one or more
markers. Genomic selection models exploit this linkage
disequilibrium to estimate the effects of the markers and
genomic breeding values [9]. Although the focus in the
development of these models has been to estimate
breeding values rather than to identify QTL, the esti-
mated marker effects show potential use for QTL map-
ping [11]. For the original models proposed by
Meuwissen et al. [9], either all markers were included in
the model (Bayes A), or the proportion of markers in-
cluded in the model was fixed at a certain value (Bayes B).
To overcome the statistical problems associated with these
models reported by Gianola et al. [12], Habier et al. [13]
proposed several modifications of the original models, in-
cluding Bayes Cπ and Bayes C. In Bayes Cπ, the propor-
tion of markers included in the model is assigned a prior
distribution and estimated during the analysis, while in
Bayes C, the proportion is given a fixed value. In a
simulation by Sun et al. [14], Bayes Cπ was successful
in identifying large QTL. Bayes C has been used to
identify QTL for various traits in beef cattle [15,16],
the horse [17] and pigs [18,19].

The performance of the Bayesian genomic selection
models is known to be influenced by the genetic archi-
tecture of the trait and the number of records used to
estimate the marker effects; the accuracy of predicted
breeding values decreases with decreasing heritability
[20,21] and an increasing number of QTL [21,22] and a
large number of records is needed to obtain high accur-
acy [4]. Although several studies have evaluated the im-
pact of these factors on the accuracy of breeding value
estimation, their impact on the use of these models for
QTL mapping is unknown.

The regions detected by QTL mapping generally con-
tain a large number of polymorphisms that could be the
causative mutation. However, assuming only one true
causative mutation in a QTL, an individual's genotype
for the causative mutation should be in concordance
with its status for the QTL: if the individual is homozy-
gous for the QTL, it should be homozygous for the
causative mutation, and if the individual is heterozygous
for the QTL, it should also be heterozygous for the
causative mutation [23]. By comparing the status of the
QTL and polymorphisms in the QTL region, all

polymorphisms with statuses that are not in concord-
ance with the statuses of the QTL can be eliminated as
potential causative mutations. Using this approach,
Karlsson et al. [24] were able to eliminate more than half
of the candidate polymorphisms for a causative mutation
involved in coat colour in dogs. However, if non-
concordant polymorphisms are eliminated, it is important
that prediction of the QTL status is correct. One
possibility to determine the QTL status is to compare
the estimated marker effects of the two haplotypes that
an individual carries in the QTL region: if the effects of
the haplotypes are very different, the individual is likely
heterozygous for the QTL, while if the difference is close
to zero, it is likely to be homozygous.

The first objective of this study was to evaluate by simu-
lation the influence of heritability, number of QTL and
number of records on the accuracy of QTL mapping with
genomic prediction methods Bayes Cπ and Bayes C. Our
second objective was to estimate the accuracy of the pre-
diction of the status of the causative mutation using the
estimated marker effects from Bayes C. The study fo-
cussed on the largest QTL because they are the natural
candidates for molecular identification of a causative mu-
tation. It aimed at answering the following questions: are
we able to detect the largest QTL? Are they true QTL?
Are we able to infer correctly the QTL status of the
individuals?

Methods

Simulation

QTL and phenotypes were simulated for different herita-
bilities and numbers of QTL. An existing real dataset
was used as a base for the simulation. It consisted of
4387 genotyped French Holstein bulls, with a pedigree
of 12 142 individuals. The genotypes were obtained with
the Illumina Bovine SNP50 Bead Chip® by Labogena.
After removal of SNPs with a minor allele frequency
below 5%, 38 277 SNPs were retained for analysis. Simu-
lations were carried out in nine scenarios, with heritabil-
ity (h²) of 0.1, 0.4 or 0.7 and 10, 100 or 1000 biallelic
QTL, in both the full dataset and a subset consisting of
1500 randomly selected individuals from the full dataset.
Table 1 summarizes the nine simulation scenarios. For
the scenarios with 100 QTL and a heritability of 0.4 or
0.7, the simulation was run 11 times. It was run only
once for the others.

The simulated QTL were divided into three groups: large, medium and small QTL. In the scenarios with 10 QTL, there were one large, two medium and seven small QTL, in the scenarios with 100 QTL, five large, 20 medium and 75 small QTL, and in the scenarios with 1000 QTL, 50 large, 200 medium and 750 small QTL. With the high QTL density in the scenarios with 1000 QTL, the goal was to assess whether the model could
<table><caption>**Table 1 Simulated scenarios**</caption><thead><tr><th>Number of QTL</th><th>10</th><th>100</th><th>1000</th></tr></thead><tbody><tr><td>Number of large QTL</td><td>1</td><td>5</td><td>50</td></tr><tr><td>Number of medium QTL</td><td>2</td><td>20</td><td>200</td></tr><tr><td>Number of small QTL</td><td>7</td><td>75</td><td>750</td></tr><tr><td><strong>Heritability</strong></td><td><strong>0.1</strong></td><td><strong>0.4</strong></td><td><strong>0.7</strong></td><td><strong>0.1</strong></td><td><strong>0.4</strong></td><td><strong>0.7</strong></td><td><strong>0.1</strong></td><td><strong>0.4</strong></td><td><strong>0.7</strong></td></tr><tr><td>Variance of large QTL</td><td>2.86</td><td>11.45</td><td>20.05</td><td>0.35</td><td>1.40</td><td>2.45</td><td>0.04</td><td>0.14</td><td>0.25</td></tr><tr><td>Variance of medium QTL</td><td>0.95</td><td>3.82</td><td>6.68</td><td>0.12</td><td>0.47</td><td>0.82</td><td>0.01</td><td>0.05</td><td>0.08</td></tr><tr><td>Variance of small QTL</td><td>0.32</td><td>1.27</td><td>2.23</td><td>0.04</td><td>0.16</td><td>0.27</td><td>0.00</td><td>0.02</td><td>0.03</td></tr></tbody></table>

Number of large, medium, and small QTL, with the variance explained by each QTL; the phenotypic variance was 100 and, together, the simulated QTL explained 70% of the genetic variance.

distinguish the larger QTL from the smaller ones. Together, the QTL explained 70% of the genetic variance. Large and medium QTL explained respectively 9 and 3 times the variance of a small QTL. The locations of the QTL were determined by random sampling of the SNP number. The QTL were subsequently discarded from the genotype file used in later analyses to estimate the SNP effects. The phenotypic variance $\sigma_p^2$ was set at 100.

Allele substitution effects ($a_k$) were scaled as follows. The frequency $p_k$ of one of the alleles of the $k^{th}$ QTL was computed in order to adjust QTL effect $a_k$ such that $2p_k(1-p_k)a_k^2$ equalled the part of the variance explained by the QTL. To generate residual polygenic effects, breeding values of founder animals were drawn from a normal distribution $N(0, \sigma_u^2)$ with polygenic variance $\sigma_u^2$ set to 30% of the total genetic variance. Breeding values of individuals with known parents were drawn from a normal distribution $N(0.5(u_s + u_d), 0.5\sigma_u^2)$ where $u_s$ and $u_d$ equal the breeding values of an individual's sire and dam. Breeding values of individuals with a known sire but unknown dam were drawn from a normal distribution $N(0.5u_s, 0.75\sigma_u^2)$. The simulated performance of an individual was calculated as the sum of its polygenic value, the sum of its QTL effects and an environmental effect drawn from a normal distribution $N(0, (1-h^2)\sigma_p^2)$.

## QTL detection

The effect of each SNP was estimated using methods Bayes C$\pi$ and Bayes C [13] as implemented in the GS3 software [25]. For both methods, the statistical model was:

$$y_i = \mu + u_i + \sum_{k=1}^{K} z_{ik} + e_i,$$

where $y_i$ is the phenotype for individual $i$, $\mu$ the overall mean, $K$ the number of markers, $z_{ik}$ the genotype of individual $i$ for marker $k$ coded 0, 1 or 2 depending on the

number of copies of a given marker allele the individual carried, $a_k$ the additive effect of marker $k$, and $e_i$ the random residual for individual $i$.

All unknown parameters were assigned prior distributions and sampled with a Monte Carlo Markov chain (MCMC) using Gibbs sampling. The MCMC was run for 180 000 iterations, with a burn-in of 20 000 iterations and thin interval of 50. The prior used for $a_k$ was a mixture distribution that equalled:

$$a_k | \pi, \sigma_a^2 \sim \begin{cases} 0 & \text{with probability } \pi, \\ N(0, \sigma_a^2) & \text{with probability } (1-\pi), \end{cases}$$

Where $\sigma_a^2$ is the common marker effect variance and the hyper parameter $\pi$ is the prior probability that the effect of marker $k$ is 0. Variances $\sigma_u^2$, $\sigma_a^2$ and $\sigma_e^2$ were assigned inverted chi-square distributions with $\nu = 4.2$ degrees of freedom and scale parameter $S^2 = \frac{\hat{\sigma}^2(\nu-2)}{\nu}$ where $\hat{\sigma}^2$ equals the prior value of $\sigma_u^2$, $\sigma_a^2$ or $\sigma_e^2$. When Bayes C$\pi$ was used, $\pi$ was assigned a uniform distribution and sampled with the MCMC. For Bayes C, $\pi$ was fixed at 0.99, based on the estimates obtained with Bayes C$\pi$ in the scenarios for which $\pi$ converged.

As QTL detection criterion, the posterior inclusion probability (PIP) was used, which is the proportion of iterations that included a specific marker in the model. Detection was based on either the PIP of individual markers, or the PIP summed over intervals of adjacent markers. Different interval sizes were tested, consisting of 20, 40, 60 or 80 markers, corresponding on average across the genome respectively to 1.3, 2.6, 3.9 and 5.3 Mb. Because our goal was to identify large QTL rather than all QTL, the 10 markers with the highest PIP that were more than 4 Mb apart or the 10 non-overlapping intervals with the highest PIP were declared QTL. In the simulation, a detected QTL was denoted false positive if there was no true QTL within a distance of 2 Mb on either side of the marker or within the interval, depending on the criterion that was used.
## Prediction of the status at the causative mutation

Prior to analysis, genotypes were phased in order to define haplotypes. Phasing was done using DagPhase [26], accounting for both family structure and population linkage disequilibrium. To predict the status of an individual for a simulated causative mutation, first the estimated effects of all markers in an identified QTL interval were summed up for each haplotype that the individual carried. Second, the absolute value of the difference between the estimated effects of an individual's two haplotypes was used to group individuals by partitioning around medoids (PAM) [27], as implemented in the fpc R-package [28], where the number of clusters (n) is estimated on the basis of optimum average silhouette [29]. The PAM algorithm was as follows:

1. *n* medoids were randomly selected from the data.

2. All non-medoids were assigned to the closest medoid and the costs of configurations, when the medoids and data points were switched, were calculated.

3. The configuration with the lowest cost was selected.

4. Steps 2 and 3 were repeated until there was no longer any change in the medoids.

A range of 2 to 4 was used to estimate *n*. The status of animals in the cluster with the lowest haplotype difference was denoted homozygous, and that of animals in the cluster with the highest difference was denoted heterozygous. If more than two clusters were present, the status of animals in the other clusters was denoted unknown. To evaluate the accuracy of status prediction, predicted statuses were compared with the true statuses in order to compute the proportion of correct, incorrect and unknown statuses.

## Results

### Posterior distribution of π

The posterior distribution of the proportion of markers included in the Bayes Cπ model (1-π) for the full dataset is shown in Figure 1. Only in the scenarios with 10 or 100 QTL and a heritability of 0.4 or 0.7, did the posterior distribution of π show a strong peak. In the other scenarios, values taken by π ranged from 0 to 1. In the subset data, a clear peak in the posterior distribution of π was only observed with a heritability of 0.7 and 10 or 100 QTL. For scenarios without a clear peak, the average PIP was much higher than for scenarios with a strong peak. As shown in Figure 2, for a scenario without a strong peak, the average PIP was high and the peaks were much smaller and broader for Bayes Cπ than for Bayes C. Table 2 shows the average number of markers included in the model in the full dataset. For all

scenarios, the number of included markers was much higher than the number of simulated QTL.

### Influence of heritability and number of records

The accuracy of QTL detection was influenced by both heritability and the number of simulated QTL. Table 3 shows the numbers of correctly identified QTL and false positives for the nine scenarios, both when intervals consisting of 60 subsequent markers were selected on the basis of the sum of the PIP of all markers in the interval, and when markers were individually selected according to their PIP. Both for selection of intervals and of markers, the number of correctly identified QTL and the number of false positives were highly affected by heritability and number of QTL. The number of false positives decreased with increasing heritability. Furthermore, the number of false positives was higher with 100 QTL than with 10 QTL. With 1000 QTL, due to the large number of QTL, there was on average 1.6 QTL present per haplotype region and thus hardly any false positives were found. Most detected QTL were, however, small, meaning that most large QTL were not in the 10 selected regions.

In most cases, selection of intervals resulted in fewer or equal false positives than selection of markers, except in the scenario with 10 QTL and a heritability of 0.7, which had two and one false positives for selection based on intervals and markers, respectively.

Generally, the proportion of large QTL detected was greater than the proportions of medium and small QTL detected. However, with 1000 QTL, hardly any large QTL were detected when intervals were selected but larger QTL were detected when selecting individual markers. For example, with 1000 QTL and a heritability of 0.7, selecting intervals led to the detection of only one large QTL, while selection of markers detected five large QTL.

For the scenarios with 100 QTL and a heritability of 0.4 and 0.7, the simulation was run ten times. The results of these repeated simulations are shown in Table 4. The mean number of false positives over ten runs was 0.40 and 0.10 for heritabilities of 0.4 and 0.7, respectively.

### Number of records, interval size and number of detected QTL

The number of records had a strong influence on the accuracy of QTL detection. Figure 3 shows the number of false positives for both the full dataset and the subset. The number of false positives was higher for the subset than for the full dataset. The difference was largest in the scenarios with 10 QTL and a heritability of 0.4, which had only one false positive with the full dataset but five false positives with the subset.
<figure><img src="image_3.png" /><figcaption>**Figure 1** Posterior distribution of the proportion of markers included in Bayes $C\pi$ (1-$\pi$) for simulated scenarios. The MCMC chain was run for 180 000 iterations; scenarios differed in the number of QTL (10, 100, or 1000) and heritability ($h^2$ = 0.1, 0.4, or 0.7).</figcaption></figure>

Figure 4 shows the effect of interval size on the number of false positives. In the scenarios with 1000 QTL, the genome was saturated with QTL, and there were hardly any intervals without a QTL. As we were interested in identifying the largest QTL, we redefined the number of false positives as the number of detected intervals without QTL or with only small QTL for these scenarios. The number of false positives decreased with increasing interval size. This decrease was stronger when the number of QTL was greater, and was most pronounced with a heritability of 0.4.

Figure 5 shows the number of false positives when more than 10 intervals were selected. The number of false positives in the scenarios with 1000 QTL was defined as previously described. In all scenarios, the

number of false positives increased when more intervals were selected. This increase was lower with a higher heritability, especially with 100 QTL.

### Prediction of status at the causative mutation

Figure 6 shows the proportion of correct, incorrect and undefined statuses for the scenarios with 100 QTL and a heritability of 0.4 and 0.7. On average, QTL status prediction was correct for approximately 50% of the individuals but this percentage differed strongly between QTL, ranging from 20.4% to 97.7%. Size of the QTL did not affect the percentage of correctly predicted statuses. The shape of the distributions of the haplotype differences varied largely between QTL. As shown in Figure 7, for a QTL with status predicted with almost 100%
<figure><img src="image_4.png" /><figcaption>**Figure 2** Posterior inclusion probability (p) across the genome for Bayes Cπ and Bayes C. Scenario with 100 simulated QTL and heritability equal to 0.1; the red, orange and green dots on the x-axes show the true positions of the large, medium, and small simulated QTL, respectively.</figcaption></figure>

accuracy, the homozygous and heterozygous individuals formed two clearly separated groups, while for the QTL that gave the poorest status predictions, there was complete overlap between homozygous and heterozygous individuals. The accuracy of status prediction was affected by the number of QTL in the neighbourhood of the causative mutation. Figure 8 shows that the proportion of correctly predicted statuses was highest when there was no other QTL within 10 Mb of the causative mutation, although this effect was larger with a heritability of 0.4 than with a heritability of 0.7. Furthermore, when there was only one QTL present, clustering almost always resulted in two groups and no unknowns, while with more QTL close to one another, the individuals were divided into three or four groups, so that status was undefined for some animals.

## Discussion

The accuracy of QTL mapping was affected by heritability, number of QTL, and number of records. In the subset data, the number of false positives was very high in all scenarios, demonstrating the need for a large dataset. With a low heritability, the QTL effects were smaller than with a higher heritability, making it more difficult to accurately locate the QTL and thus resulting in a greater number of false positives. Therefore, better

results will be obtained when progeny-based evaluations are used rather than individual phenotypes, since the heritability of progeny means is higher. This strategy based on daughter-yield-deviations or deregressed estimated breeding values is quite frequently used in dairy cattle [30].

When the number of QTL to explain a given amount of genetic variance increases, the individual QTL effects become smaller and the QTL thus become more difficult to detect and locate. Consequently, the number of false

<table><caption>**Table 3** Correctly identified QTL and false positives using Bayes C</caption><thead><tr><th rowspan="3">Size and<br>number of<br>simulated QTL</th><th colspan="6">Heritability</th></tr><tr><th colspan="3">Selection of intervals</th><th colspan="3">Selection of markers</th></tr><tr><th>0.1</th><th>0.4</th><th>0.7</th><th>0.1</th><th>0.4</th><th>0.7</th></tr></thead><tbody><tr><td>Large</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td><td>1</td></tr><tr><td>Medium</td><td>2</td><td>0</td><td>2</td><td>2</td><td>0</td><td>1</td></tr><tr><td>Small</td><td>7</td><td>4</td><td>4</td><td>5</td><td>3</td><td>5</td></tr><tr><td>False positives</td><td></td><td>5</td><td>3</td><td>2</td><td>6</td><td>3</td></tr><tr><td>Large</td><td>5</td><td>1</td><td>2</td><td>4</td><td>1</td><td>2</td></tr><tr><td>Medium</td><td>20</td><td>1</td><td>6</td><td>6</td><td>2</td><td>2</td></tr><tr><td>Small</td><td>75</td><td>1</td><td>2</td><td>2</td><td>1</td><td>1</td></tr><tr><td>False positives</td><td></td><td>7</td><td>3</td><td>2</td><td>7</td><td>5</td></tr><tr><td>Large</td><td>50</td><td>0</td><td>1</td><td>1</td><td>1</td><td>4</td></tr><tr><td>Medium</td><td>200</td><td>7</td><td>9</td><td>8</td><td>5</td><td>7</td></tr><tr><td>Small</td><td>750</td><td>8</td><td>16</td><td>16</td><td>14</td><td>16</td></tr><tr><td>False positives</td><td></td><td>1</td><td>0</td><td>0</td><td>1</td><td>0</td></tr></tbody></table>

Number of correctly identified QTL categorised by size of the QTL variance (large, medium or small), and number of false positives, with QTL detection based on the posterior inclusion probability summed up over intervals of 60 adjacent markers or the posterior inclusion probability of individual markers; the 10 best intervals or markers were denoted as QTL; more than one QTL can be present in an interval, so that the total number of QTL identified can be larger than 10.

<table><caption>**Table 2** Posterior means of the number of markers included in Bayes Cπ</caption><thead><tr><th rowspan="3">Number<br>of QTL<br>simulated</th><th colspan="3">Heritability</th><th colspan="3">Subset</th></tr><tr><th colspan="3">Full data</th><th colspan="3">Subset</th></tr><tr><th>0.1</th><th>0.4</th><th>0.7</th><th>0.1</th><th>0.4</th><th>0.7</th></tr></thead><tbody><tr><td>10</td><td>25 173</td><td>503</td><td>43</td><td>28 633</td><td>24 782</td><td>1374</td></tr><tr><td>100</td><td>29 504</td><td>743</td><td>418</td><td>17 663</td><td>27 072</td><td>11 906</td></tr><tr><td>1000</td><td>20 698</td><td>26 612</td><td>25 249</td><td>19 681</td><td>21 766</td><td>26 895</td></tr></tbody></table>
<table><caption>**Table 4** Correctly identified QTL and false positives using Bayes C over 10 replicates</caption><thead><tr><th rowspan="2"></th><th colspan="4">Heritability = 0.4</th><th colspan="4">Heritability = 0.7</th></tr><tr><th>Large</th><th>Medium</th><th>Small</th><th>False</th><th>Large</th><th>Medium</th><th>Small</th><th>False</th></tr></thead><tbody><tr><th>Mean</th><td>3.70</td><td>4.5</td><td>3.60</td><td>0.40</td><td>2</td><td>5.4</td><td>8.5</td><td>0.1</td></tr><tr><th>Min</th><td>3</td><td>3</td><td>2</td><td>0</td><td>0</td><td>3</td><td>6</td><td>0</td></tr><tr><th>Max</th><td>4</td><td>8</td><td>6</td><td>2</td><td>3</td><td>8</td><td>12</td><td>1</td></tr><tr><th>SD</th><td>0.48</td><td>1.78</td><td>1.51</td><td>0.70</td><td>0.94</td><td>1.90</td><td>2.55</td><td>0.32</td></tr></tbody></table>

Mean, minimum, maximum and standard deviation of the number of correctly identified QTL categorised by effect size (large, medium or small), and number of false positives, based on 10 replicates; QTL detection was based on the posterior inclusion probability summed up over intervals of 60 adjacent markers; the best 10 intervals were denoted as QTL; 100 QTL were simulated, of which 5 were large, 20 medium and 75 small.

positives increased when the number of QTL increased from 10 to 100. However, with a further increase to 1000 QTL, the number of false positives decreased to 1 or 0 because with 1000 QTL on a total of 37 227 markers, an interval of 60 markers includes on average 1.6 QTL. Randomly selecting intervals would give the same result. The motivation for QTL mapping is the identification of large rather than small QTL. In the scenarios with 1000 QTL, mainly small QTL were found. It is, however, important to realise that large QTL in the scenarios with 1000 QTL had about the same size as the small QTL in the scenarios with 100 QTL and were much smaller than the small QTL in the scenarios with 10 QTL (Table 1). For example, with a heritability of 0.4, a large QTL in the scenario with 1000 QTL explains 0.14% of the total phenotypic variance, while small QTL in the scenarios with 100 and 10 QTL explain 0.16% and 1.27%, respectively. The small size of the large QTL in the scenarios with 1000 QTL clearly makes it difficult to detect these QTL. Furthermore, with the large number of small and medium QTL in these scenarios, several small and medium QTL together can have a larger effect than a single large QTL. Based on our results, the

variance explained by the QTL in an interval should roughly explain at least 1% of the phenotypic variance in order for the QTL to be detected.

QTL detection can be based on several criteria, including marker effect, variance and PIP, either for each marker individually or summed up over an interval of adjacent markers. In this study, PIP was used as the detection criterion, for both markers and intervals. If several markers are in linkage disequilibrium with a QTL, the QTL effect can spread over several markers. In that case, selection based on intervals can give more accurate results than selection based on individual markers [10]. In most scenarios evaluated here, fewer false positives were found when intervals were selected than when markers were selected. However, as shown in Table 3, selection of intervals resulted in the detection of mainly small QTL and hardly any large QTL in the scenarios with 1000 QTL, while there were 50 large QTL simulated and only the best 10 intervals were selected. In contrast, when detection was based on individual markers, a much larger proportion of the detected QTL were large QTL. The intervals that were selected in this case were generally located in the proximity of several medium or large QTL. Selection of intervals can thus result in the identification of regions that are close to the largest number of QTL rather than regions that include the largest QTL. Therefore, for polygenic traits, selection of markers can give better results, in terms of identifying large QTL, than selection of intervals. Since the number of QTL influencing a trait is generally unknown, it is difficult to decide which approach to use, and the best option could be to combine both approaches and declare only QTL that have a high ranking for both criteria.

When the interval size increased, the number of false positives decreased, as shown in Figure 4. However, increasing interval size makes the subsequent search for causative mutations more difficult, since the search area

<figure><img src="image_5.png" /></figure>
<figure><img src="image_6.png" /><figcaption>**Figure 4** Number of false positives as a function of interval size. The ten intervals with the highest summed posterior inclusion probabilities were selected.</figcaption></figure>

increases as well. Based on our results, interval sizes containing 40 to 60 markers and corresponding to an average 2.6 to 3.9 Mb, seem to be a reasonable compromise for the usual designs currently available in animal breeding.

In results reported here, the 10 best markers or intervals were selected and denoted as QTL, which is only a small fraction of the total number of QTL in the scenarios with 100 and 1000 QTL. However, when more QTL are selected, the number of false positives increases rapidly. This is in agreement with other studies, in which Bayes A and B [9] and Bayes Cπ [14] were successful in locating the large QTL but missed the smaller ones. The Bayesian genomic selection models are thus suitable for QTL mapping when the aim is to find the major QTL.

In addition to previously discussed factors, linkage disequilibrium between markers and QTL is required for the markers to pick up the effect of a QTL. The accuracy of genomic prediction with Bayesian models is lower in populations with a large effective population size and thus a low amount of linkage disequilibrium [4]. In this study, data were simulated using genotypes of French Holstein cattle. When applied to populations with a different effective population size, the accuracy of QTL mapping is expected

to be different than observed here. In contrast to genomic prediction, the motivation for QTL mapping is not to capture the maximum amount of information on the trait but rather to capture information on markers close to the causative mutation. In a population in which linkage disequilibrium is present only over short distances, only markers near the causative mutation will have an effect, and thus a more accurate estimate of the QTL location is expected than with a population with long range linkage disequilibrium.

The posterior distribution of π was affected by heritability, number of QTL and number of records. In scenarios for which the posterior distribution of π lacked a clear peak were the scenarios with a low heritability or a large number of QTL. In these scenarios, the variance explained by an individual QTL was small, increasing the difficulty of locating the QTL. For scenarios with a clear peak, the number of SNPs included in the model was much larger than the number of simulated QTL. This agrees with results of Habier et al. [13] who reported that the number of SNPs with effect estimated on the basis of π was overestimated. In their results, the number of markers included in the model increased with decreasing heritability. The accuracy of breeding value

<figure><img src="image_7.png" /><figcaption>**Figure 5** Number of false positives according to the number of selected intervals. Intervals consisted of 60 adjacent markers.</figcaption></figure>
<figure><img src="image_8.png" /><figcaption>**Figure 6** Proportion of correctly and incorrectly predicted and undefined statuses at the detected QTL. For scenarios with heritabilities (h²) of 0.4 and 0.7, 100 QTL were simulated 10 times, of which five were large (L), 20 medium (M) and 75 small (S).</figcaption></figure>

predictions from Bayesian genomic selection models de-
creases and becomes more similar to the accuracy of GBLUP
with decreasing heritability, increasing number of QTL
and decreasing number of records, [4,20-23]. The GBLUP
model uses all markers to estimate breeding values on
the basis of the genomic relationship between animals
rather than based on individual marker effects [31].
In scenarios that did not show a clear peak in the pos-
terior distribution of π, individual marker effects were
small and thus more difficult to estimate. A possible
explanation for the lack of a clear peak could be that in
these cases, Bayes Cπ acts similar to the GBLUP model
and relies more on the genomic relationship between ani-
mals than on QTL effects, and therefore includes a large
and variable proportion of markers in the model to cap-
ture relationships between animals. As a result, the aver-
age PIP of markers was high, and the peaks of PIP were
relatively small deviations from this mean. In addition, a
strong negative correlation was observed between π and
σ². When π was fixed, the deviations from the mean were

larger and sharper, as in Figure 2, making it easier to
locate QTL.

The accuracy of the predicted statuses of the causative
mutations differed strongly between QTL. For some
QTL, the accuracy was almost 100%, while for others the
status of almost none of the individuals was correctly
predicted. The accuracy of status predictions was affected
by the presence of other QTL nearby. When the effect of
the haplotypes is influenced by several QTL rather than by
only one, the status of the QTL can be predicted correctly
based on the difference between haplotype effects only if
the statuses of all QTL influencing the haplotype effect are
in concordance with each other. When this is not the case,
the differences between haplotype effects of homozygous
and heterozygous individuals overlap and it is not possible
to cluster the individuals in two groups. Although the pre-
diction of QTL status was influenced by the presence of
other QTL nearby, their presence seemed to explain only
part of the variation in prediction accuracy between QTL.
It is not clear what caused the rest of the variation.

<figure><img src="image_9.png" /><figcaption>**Figure 7** Distribution of absolute difference between estimated haplotype effects in an identified QTL region. Haplotype effects were estimated by the sum of the estimated marker effects for all markers in the haplotype; distribution for the QTL with the highest (left) and lowest (right) status prediction accuracy.</figcaption></figure>
<figure><img src="image_10.png" /><figcaption>**Figure 8** Effect of the number of QTL on accuracy of status prediction for the causative mutation. Box plots showing the proportion of correctly and incorrectly predicted and undefined statuses, for scenarios with heritabilities (h²) of 0.4 and 0.7; in both scenarios, 100 QTL were simulated 10 times, of which five were large (L), 20 medium (M) and 75 small (S); the x-axis shows the number of QTL within a distance of 10 Mb of one another.</figcaption></figure>

In order to identify the causative mutation, it is important that the predicted QTL status of individuals is correct. Otherwise the true polymorphism could be eliminated because of lack of concordance with the QTL status. Although the status prediction was far from perfect, histograms of the haplotype differences clearly demonstrated whether the clustering was successful or not. As shown in Figure 7, when two clearly separated groups were visible, almost 100% of the statuses were correctly predicted, while histograms did not show separate groups for QTL with low prediction accuracy.

## Conclusions

The performance of Bayes Cπ and Bayes C was highly affected by heritability, number of QTL and number of records. The number of false positives was high with low heritability, and decreased with increasing heritability. In polygenic scenarios, larger QTL were detected when detection was based on the PIP of individual markers than when based on the PIP summed up over a group of adjacent markers. Thus, a first conclusion is that Bayes C is a suitable method for QTL mapping when it is applied to traits with moderate to high heritability and when a large number of records are available.

For Bayes Cπ, there was no strong peak in the posterior distribution of π when heritability was low or the number of QTL was large, and a large number of SNPs were selected. In such cases, localisation of the QTL was more accurate with Bayes C. Furthermore, with moderate heritability, a large number of records was needed for π to show a clear peak in the posterior distribution. Therefore, a second conclusion is to recommend a high fixed value for π.

The accuracy of the predicted status of the causative mutation is affected by the size of the QTL, as well as by the presence of other QTL nearby. Although the accuracy in general is not very high, a histogram of the haplotype differences can be used to decide whether the clustering is correct or not and thus usable for further analysis.

## Competing interests

The authors declare that they have no competing interests.

## Authors' contributions

IB and DB designed and carried out the study, and drafted the manuscript. SF generated and provided phased data. All authors read and approved the final manuscript.
## Acknowledgements

IB benefited from an Erasmus-Mundus fellowship and a grant by Apisgene, within the framework of the European Graduate School in Animal Breeding and Genetics. This work was part of the "Rules and Tools" project, financed by the French National Research Agency (ANR-09-GENM-002-01). Genotype data originated from the Cartofine project funded by ANR and Apisgene.

## Author details

¹INRA, UMR1313 Génétique animale et biologie intégrative, Domaine de Vilvert, 78350 Jouy-en-Josas, France. ²AGROPARISTECH, UMR1313 Génétique animale et biologie intégrative, 16 rue Claude Bernard, 75231 Paris 05, France. ³Department of Molecular Biology and Genetics, Faculty of Science and Technology, Aarhus University, DK-8830 Tjele, Denmark. ⁴UNCEIA, 149 rue de Bercy, 75012 Paris, France.

Received: 9 November 2012 Accepted: 7 June 2013

Published: 19 June 2013

## References

1. Khatkar MS, Thomson PC, Tammes I, Raadsma HW: *Quantitative trait loci mapping in dairy cattle: review and meta-analysis*. *Genet Sel Evol* 2004, **36**:163-190.

2. Abasht B, Dekkers JCM, Lamont SJ: *Review of quantitative trait loci identified in the chicken*. *Poultry Sci* 2006, **85**:2079-2096.

3. Rothschild MF, Hu ZL, Jiang Z: *Advances in QTL mapping in pigs*. *Int J Biol Sci* 2007, **3**:192-197.

4. Goddard ME, Hayes BJ: *Mapping genes for complex traits in domestic animals and their use in breeding programmes*. *Nat Rev Genet* 2009, **10**:381-391.

5. Boichard D, Grohs C, Bourgeois F, Cerqueira F, Faugeras R, Neau A, Rupp R, Amigues Y, Boscher MY, Levéziel H: *Detection of genes influencing economic traits in three French dairy cattle breeds*. *Genet Sel Evol* 2003, **35**:77-101.

6. Braunschweig MH: *Mutations in the bovine ABCG2 and the ovine MSTN gene added to the few quantitative trait nucleotides identified in farm animals: a mini-review*. *J Appl Genet* 2010, **51**:289-297.

7. Jansen RC, Stam P: *High resolution of quantitative traits into multiple loci via interval mapping*. *Genetics* 1994, **136**:1447-1455.

8. Yi N, Shriner D: *Advances in Bayesian multiple quantitative trait loci mapping in experimental crosses*. *Heredity* 2008, **100**:240-252.

9. Meuwissen THE, Hayes BJ, Goddard ME: *Prediction of total genetic value using genome-wide dense marker maps*. *Genetics* 2001, **157**:1819-1829.

10. Sahana G, Guldbrandtsen B, Janss L, Lund MS: *Comparison of association mapping methods in a complex pedigreed population*. *Genet Epidemiol* 2010, **34**:455-462.

11. Hayes BJ, Pryce J, Chamberlain AJ, Bowman PJ, Goddard ME: *Genetic architecture of complex traits and accuracy of genomic prediction: coat colour, milk-fat percentage and type in Holstein cattle as contrasting model traits*. *PLoS Genet* 2010, **6**:e1001139.

12. Gianola D, Delos Campos G, Hill WG, Manfredi E, Fernando R: *Additive genetic variability and the Bayesian alphabet*. *Genetics* 2009, **183**:347-363.

13. Habier D, Fernando RL, Kizilkaya K, Garrick DJ: *Extension of the Bayesian alphabet for genomic selection*. *BMC Bioinformatics* 2011, **12**:186.

14. Sun X, Habier D, Fernando RL, Garrick DJ, Dekkers JCM: *Genomic breeding value prediction and QTL mapping of QTLMAS 2010 data using Bayesian methods*. *BMC Proc* 2011, **5**:S13.

15. Peters SO, Kizilkaya K, Garrick DJ, Fernando RL, Reecy JM, Weaber RL, Silver GA, Thomas MG: *Bayesian genome-wide association analysis of growth and yearling ultrasound measures of carcass traits in Brangus heifers*. *J Anim Sci* 2012, **90**:3398-3409.

16. Peters SO, Kizilkaya K, Garrick DJ, Fernando RL, Reecy JM, Weaber RL, Silver GA, Thomas MG: *Heritability and Bayesian genome-wide association study of first service conception and pregnancy in Brangus heifers*. *J Anim Sci* 2013, **91**:605-612.

17. Schurink A, Wolc A, Ducro BJ, Frankena K, Garrick DJ, Dekkers JCM, van Arendonk JAM: *Genome-wide association study of insect bite hypersensitivity in two horse populations in the Netherlands*. *Genet Sel Evol* 2012, **44**:31.

18. Onteru SK, Fan B, Nikkilä MT, Garrick DJ, Stalder KJ, Rothschild MF: *Whole-genome association analyses for lifetime reproductive traits in the pig*. *J Anim Sci* 2011, **89**:988-995.

19. Fan B, Onteru SK, Du ZQ, Garrick DJ, Stalder KJ, Rothschild MF: *Genome-wide association study identifies loci for body composition and structural soundness traits in pigs*. *PLoS One* 2011, **6**:e14726.

20. Calus MPL, Veerkamp RF: *Accuracy of breeding values when using and ignoring the polygenic effect in genomic breeding value estimation with a marker density of one SNP per cM*. *J Anim Breed Genet* 2007, **124**:362-368.

21. Daetwyler HD, Pong-Wong R, Villanueva, Woolliams JA: *The impact of genetic architecture on genome-wide evaluation methods*. *Genetics* 2010, **185**:1021-1031.

22. Clark SA, Hickey JM, van der Werf JH: *Different models of genetic variation and their effect on genomic evaluation*. *Genet Sel Evol* 2011, **43**:18.

23. Ron M, Weller JL: *From QTL to QTN identification in livestock – winning by points rather than knock-out: a review*. *Anim Genet* 2007, **38**:429-439.

24. Karlsson EK, Baranowska I, Wade CM, Salmon, Hillbertz NH, Zody MC, Anderson N, Biagi TM, Patterson N, Pielberg GR, Kulbokas EJ 3rd, Comstock KE, Keller ET, Mesirov JP, Von Euler H, Kämpe O, Hedhammar A, Lander ES, Andersson G, Andersson L, Lindblad-Toh K: *Efficient mapping of Mendelian traits in dogs through genome-wide association*. *Nat Genet* 2007, **39**:1321-1328.

25. Legarra A, Ricard A, Filangi O: *GS3*. User Manual; 2012 [http://snp.toulouse.inra.fr/~alegarra/manualgs3_last.pdf]

26. Druet T, Georges M: *A hidden markov model combining linkage and linkage disequilibrium information for haplotype reconstruction and quantitative trait locus fine mapping*. *Genetics* 2010, **184**:789-798.

27. Kaufman L, Rousseeuw PJ: *Finding Groups in Data: An Introduction to Cluster Analysis*. Hoboken: John Wiley & Sons; 2005.

28. Hennig C: *fpc: Flexible procedures for clustering*. R package version 2.0-3; 2010. http://cran.rproject.org/web/packages/fpc/.

29. Rousseeuw PJ: *Silhouettes: a graphical aid to the interpretation and validation of cluster analysis*. *J Comput Appl Math* 1987, **20**:53-65.

30. Mrode RA, Swanson GJT: *Calculating cow and daughter yield deviations and partitioning of genetic evaluations under a random regression model*. *Livest Prod Sci* 2004, **86**:253-260.

31. Hayes BJ, Visscher PM, Goddard ME: *Increased accuracy of artificial selection by using the realized relationship matrix*. *Genet Res* 2009, **91**:47-60.

doi:10.1186/1297-9686-45-19

Cite this article as: van den Berg et al.: QTL fine mapping with Bayes C (π): a simulation study. *Genetics Selection Evolution* 2013 **45**:19.

Submit your next manuscript to BioMed Central and take full advantage of:

* Convenient online submission

* Thorough peer review

* No space constraints or color figure charges

* Immediate publication on acceptance

* Inclusion in PubMed, CAS, Scopus and Google Scholar

* Research which is freely available for redistribution

Submit your manuscript at
www.biomedcentral.com/submit

BioMed Central