# Combining two Meishan F2 crosses improves the detection of QTL on pig chromosomes 2, 4 and 6

Flavie Tortereau, Hélène Gilbert, Henri Cm Heuven, Jean Pierre Bidanel, Martien Am Groenen, Juliette Riquet

► To cite this version:

Flavie Tortereau, Hélène Gilbert, Henri Cm Heuven, Jean Pierre Bidanel, Martien Am Groenen, et al.. Combining two Meishan F2 crosses improves the detection of QTL on pig chromosomes 2, 4 and 6. Genetics Selection Evolution, 2010, 42, online (november), Non paginé. 10.1186/1297-9686-42-42. hal-01193594

HAL Id: hal-01193594

https://hal.science/hal-01193594

Submitted on 30 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_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Combining two Meishan F2 crosses improves the detection of QTL on pig chromosomes 2, 4 and 6

Flavie Tortereau¹,³, Hélène Gilbert², Henri CM Heuven³, Jean-Pierre Bidanel², Martien AM Groenen³, Juliette Riquet¹*

**Abstract**

**Background:** In pig, a number of experiments have been set up to identify QTL and a multitude of chromosomal regions harbouring genes influencing traits of interest have been identified. However, the mapping resolution remains limited in most cases and the detected QTL are rather inaccurately located. Mapping accuracy can be improved by increasing the number of phenotyped and genotyped individuals and/or the number of informative markers. An alternative approach to overcome the limited power of individual studies is to combine data from two or more independent designs.

**Methods:** In the present study we report a combined analysis of two independent design (a French and a Dutch F2 experimental designs), with 2000 F2 individuals. The purpose was to further map QTL for growth and fatness on pig chromosomes 2, 4 and 6. Using QTL-map software, uni- and multiple-QTL detection analyses were applied separately on the two pedigrees and then on the combination of the two pedigrees.

**Results:** Joint analyses of the combined pedigree provided (1) greater significance of shared QTL, (2) exclusion of false suggestive QTL and (3) greater mapping precision for shared QTL.

**Conclusions:** Combining two Meishan x European breeds F2 pedigrees improved the mapping of QTL compared to analysing pedigrees separately. Our work was facilitated by the access to raw phenotypic data and DNA of animals from both pedigrees and the combination of the two designs with the addition of new markers allowed us to fine map QTL without phenotyping additional animals.

**Background**

Over the past fifteen years, the construction of genetic maps in livestock species has enhanced efforts to dissect the molecular basis of the genetic variation of agriculturally important traits. In pig, a number of experiments have been set up to identify QTL and many chromosomal regions harbouring genes influencing traits of interest have been identified [1] and reported in QTLdb http://www.genome.iastate.edu/cgi-bin/QTLdb/index[2]. However, in most cases mapping resolution remains limited and the QTL detected are rather inaccurately located. Mapping accuracy can be improved by increasing the number of phenotyped and genotyped individuals and/or the number of informative markers. However, collecting this additional information is often

time-consuming and/or expensive. An alternative approach to overcome the limited power of individual studies is to combine data from two or more independent designs. Combining several pedigrees together increases the number of animals without additional phenotyping or genotyping costs. Without access to raw data, meta-analysis of published results can be an informative approach to increase precision. Allison and Heo [3] have proposed meta-analytical techniques that can be used under difficult conditions. However, these analyses are complicated by the differences among testing methods and experimental designs and finally, the gain in accuracy of QTL mapping is limited. Availability of the raw data to analyse jointly independent data sets is probably a better way to combine different QTL mapping designs. In pig, some studies aiming at combining pedigrees in order to increase the power of QTL detection have already been carried out. Walling et al. [4] have combined French, British, Dutch, American,

<footnote>* Correspondence: juliette.riquet@toulouse.inra.fr
¹INRA, UMR 0444 Génétique Cellulaire, F-31326 Castanet-Tolosan, France
Full list of author information is available at the end of the article</footnote>
Swedish and German studies to detect QTL on pig chromosome 4 or SSC4 (for Sus scrofa chromosome 4) and Perez-Enciso et al. [5] have combined pedigrees from Spanish, French and German designs to refine the location of a QTL for growth and fatness traits on SSCX. However, these analyses are complicated by the lack of common markers and often by slight differences in trait definitions and measurements. In addition, parental populations are usually different, and the QTL segregating in the various designs are not necessarily the same. Under these conditions, combining data sets from different origins may not be optimal to improve estimations of QTL localisations and effects.

Here we report an analysis of QTL located on SSC2, SSC4 and SSC6, that combines French and Dutch F2 pedigrees involving Meishan (MS), Large White (LW) and Landrace (LR) breeds. The analysis was focused on these three pig chromosomes because previous detection analyses [6-12] have shown that the QTL identified on these chromosomes contribute less to the global variance of the traits than QTL detected for example on SSC7 or SSCX. To optimize this joint analysis and reconstruct a unique genetic map from the 2000 F2 offspring of this combined design, additional microsatellite markers were included in either one or both pedigrees. Using single- and multiple-QTL mapping analyses on each pedigree and on the combined pedigree, we investigated the benefits of combining pedigrees (i.e. doubling the pedigree size) to refine the location of QTL for growth and fatness on SSC2, 4 and 6.

# Methods

## Pedigrees and phenotypic data

QTL mapping data from two experimental F2 crosses between European pig breeds x Meishan were used: (1) the French PORQTL pedigree produced at INRA [6], by mating six Large White sires and six Meishan dams and then six F1 sires and 20 F1 dams to produce 1052 F2 animals; all pigs were born and raised at the INRA GEPA experimental unit (Poitou-Charentes) and (2) a Dutch pedigree, obtained at the University of Wageningen (WU) [13,14] by mating 19 Meishan sires and 126 Large White or Landrace dams and then 39 F1 sires and 265 F1 dams to produce 1212 F2 offspring; this Dutch design was conducted in five different breeding companies. Among the 39 Dutch half-sib families, we selected the 24 largest families, amounting finally to 1919 French and Dutch F2 animals.

Details on the phenotypic data have been reported respectively for the French pedigree in [6] for growth and fatness traits, [15] for teat number, [8] for carcass composition traits and [7] for IntraMuscular Fat (IMF) and for the Dutch pedigree in [10-12] for growth, fatness and meat quality traits [10-12] and in [9] for teat

<table><caption>Table 1 Studied traits in the French and Dutch pedigrees</caption><thead><tr><th>Pedigree</th><th>Trait</th><th>N</th><th>Mean</th><th>SD</th></tr></thead><tbody><tr><td>French</td><td>birth weight (kg)</td><td>1052</td><td>1.25</td><td>0.27</td></tr><tr><td>French</td><td>weaning weight (kg)</td><td>1052</td><td>5.46</td><td>1.13</td></tr><tr><td>French</td><td>teat number</td><td>1046</td><td>14.82</td><td>1.56</td></tr><tr><td>French</td><td>carcass weight (kg)</td><td>529</td><td>59.11</td><td>10.48</td></tr><tr><td>French</td><td>BFT (mm)</td><td>521</td><td>17.00</td><td>5.00</td></tr><tr><td>French</td><td>LD depth (mm)</td><td>521</td><td>35.00</td><td>9.00</td></tr><tr><td>French</td><td>meat percentage (%)</td><td>521</td><td>50.44</td><td>3.66</td></tr><tr><td>French</td><td>LWG (kg/day)</td><td>960</td><td>0.49</td><td>0.10</td></tr><tr><td>French</td><td>IMF (%)</td><td>248</td><td>1.69</td><td>0.54</td></tr><tr><td>French</td><td>X2 (mm)</td><td>521</td><td>17.80</td><td>5.22</td></tr><tr><td>French</td><td>shoulder weight (kg)</td><td>489</td><td>4.71</td><td>0.84</td></tr><tr><td>French</td><td>midriff weight (kg)</td><td>489</td><td>1.16</td><td>0.33</td></tr><tr><td>French</td><td>ham weight (kg)</td><td>489</td><td>5.79</td><td>0.92</td></tr><tr><td>French</td><td>loin weight (kg)</td><td>489</td><td>8.07</td><td>1.46</td></tr><tr><td>French</td><td>leaf fat weight (kg)</td><td>489</td><td>0.42</td><td>0.23</td></tr><tr><td>French</td><td>foot weight (kg)</td><td>489</td><td>1.03</td><td>0.22</td></tr><tr><td>French</td><td>belly weight (kg)</td><td>489</td><td>2.96</td><td>0.65</td></tr><tr><td>French</td><td>kidney weight (kg)</td><td>487</td><td>0.11</td><td>0.02</td></tr><tr><td>French</td><td>head weight (kg)</td><td>481</td><td>4.81</td><td>0.82</td></tr><tr><td>Dutch</td><td>birth weight (kg)</td><td>867</td><td>1.22</td><td>0.22</td></tr><tr><td>Dutch</td><td>weaning weight (kg)</td><td>864</td><td>8.23</td><td>2.00</td></tr><tr><td>Dutch</td><td>teat number</td><td>869</td><td>15.42</td><td>1.20</td></tr><tr><td>Dutch</td><td>carcass weight (kg)</td><td>548</td><td>70.24</td><td>9.59</td></tr><tr><td>Dutch</td><td>BFT (mm)</td><td>565</td><td>22.19</td><td>5.70</td></tr><tr><td>Dutch</td><td>LD depth (mm)</td><td>563</td><td>40.82</td><td>6.72</td></tr><tr><td>Dutch</td><td>meat percentage (%)</td><td>565</td><td>48.53</td><td>4.25</td></tr><tr><td>Dutch</td><td>LWG (kg/day)</td><td>551</td><td>0.53</td><td>0.08</td></tr><tr><td>Dutch</td><td>IMF (%)</td><td>557</td><td>1.87</td><td>0.88</td></tr><tr><td>Dutch</td><td>pH_24 (LD)</td><td>565</td><td>5.67</td><td>0.27</td></tr><tr><td>Dutch</td><td>L*</td><td>563</td><td>53.85</td><td>4.73</td></tr><tr><td>Dutch</td><td>a*</td><td>565</td><td>17.13</td><td>1.83</td></tr><tr><td>Dutch</td><td>b*</td><td>565</td><td>9.52</td><td>1.91</td></tr><tr><td>Dutch</td><td>pH_24 (ham)</td><td>565</td><td>5.83</td><td>0.32</td></tr><tr><td>Dutch</td><td>driploss (%)</td><td>563</td><td>2.68</td><td>1.52</td></tr><tr><td>Dutch</td><td>cookloss (%)</td><td>564</td><td>26.35</td><td>3.46</td></tr><tr><td>Dutch</td><td>shear force (N)</td><td>564</td><td>39.20</td><td>10.75</td></tr></tbody></table>

Traits indicated in bold are common to the two independent pedigrees; N = number of analysed F2 for the corresponding trait; SD = Standard Deviation; BFT = Back-Fat Thickness (referred as X4 in [8]); LD = *Longissimus Dorsi*; LWG = Life Weight Gain; IMF = IntraMuscular Fat content; X2 is another measurement of back-fat thickness; pH_24 = pH measured 24 h after slaughter; a* = meat redness; b* = meat yellowness; L* = meat lightness

number. Nine traits related to growth, fatness and teat numbers (Table 1) were included in a joint analysis of the pedigrees. Seven of these nine traits i.e. birth weight, weaning weight, carcass weight, teat number, IMF, Back-Fat Thickness (BFT) between the 3rd and 4th ribs of a carcass at 6 cm from the spine and *Longissimus Dorsi* (LD) depth were chosen because they had been recorded in both designs with the same conditions. The two remaining common traits i.e. Life Weight Gain
(LWG) and meat percentage were already available for the Dutch pedigree and had to be computed for the French pedigree. Meat percentage was computed in the Dutch pedigree with the Hennessy Grading Probe formula taking BFT and muscle depth into account. For the French pedigree, we applied a similar formula as that used in France at the time of the experiment [16] which is also based on back fat thickness and muscle depth (meat percentage = 55.698 - 0.710 × BFT + 0.198 × LD). LWG is an average daily gain estimated throughout the entire animals' life and is calculated as weight/age.

Some additional economically important traits that had been recorded only in one of the two pedigrees but for which a significant QTL had been previously detected, were analysed only in the corresponding pedigree. These traits were related to additional fatness (X2, measured between the 3rd and the 4th lumbar vertebrae at 8 cm from the spine), cut weights (shoulder, midriff, ham, loin, leaf fat, foot, belly, kidney and head) for the French pedigree and meat quality (pH in m. *Longissimus Dorsi* and in m. *Semimembranosus* taken 24 h after slaughter, L*, a* and b* colour values of m. *Longissimus Dorsi*, driploss, cookloss and shear force) for the Dutch pedigree (Table 1). These traits are not shared by both designs, but are economically important and thus were re-analysed in this study with additional microsatellite markers.

## Genotyping

In order to compare QTL detection results among the French, Dutch and combined pedigrees, a consensus linkage map based on genotyping data from the two Meishan x European breeds F2 populations was generated. The aim was to have a density of one marker every 10 cM within the previously described QTL regions and every 20 cM on the rest of the chromosomes. QTL regions extended from the telomere of the p arm to microsatellite SW240 on SSC2, between microsatellites S0301 and S0214 on SSC4 and along two regions on SSC6 (between SW2535 and SW1057 and between S0059 and SW607). Initially, French and Dutch pedigree were genotyped over these three chromosomes with 22 and 29 microsatellite markers respectively [6,11]. Five microsatellite markers on SSC2, five on SSC4 and six on SSC6 were common to both designs. Additional informative microsatellite markers were included for one or both pedigree(s), to obtain a unique set of common markers. Among the markers genotyped on both pedigrees on SSC6, microsatellite 261M17 was specifically designed from the BAC end-sequence bT261M17SP6 with primers 5'-CTCTTCCATTCCCT-GATTGC-3' and 5'-CCCCTTCCTCACCTCTTTTCT-3' to fill the gap between S0121 (122 cM) and SW322

(152 cM). On the common map, this new microsatellite is located 12.8 cM from S0121 and 17.4 cM from SW322. Finally, for SSC2 four additional microsatellites were analysed in the French pedigree and two in the Dutch population, for SSC4, two in the French pedigree and one in the Dutch population and for SSC6, four in the French pedigree and three in the Dutch population. New genotyping data were obtained at INRA as previously described [6]. All the genotypes were validated and stored in the Gemma database [17]. Only common markers were kept in the analysis, except for S0217 and SW2466 that were used only on the Dutch animals and SW1089 and SW607 only on the French animals. Microsatellites S0217 and SW1089 and microsatellites SW2466 and SW607 which mapped to the same position respectively on SSC4 and SSC6 were considered as unique markers in the combined analysis. Marker maps were ordered using CriMap package [18], considering all the F2 animals of the two designs. The sex-averaged maps are presented in Figure 1.

## Statistical analyses

Before QTL detection, phenotypes were corrected for the usual fixed effects using a linear model (GLM procedure, SAS® 9.1, SAS Institute, Inc.). For traits previously analysed, published fixed effects and covariates were used and for the other traits, fixed effects and covariates that were significant at the 5% level in the variance analysis were kept in the final model. Corrected data showed similar variances for the traits common in both designs but recorded independently so that no standardisation was applied.

QTL analysis using these corrected data was performed with the QTLMap software developed at INRA [19,20] based on interval mapping without any hypothesis on the number of QTL alleles present in the Meishan and European breeds.

The test statistic was computed as the ratio of likelihoods under the hypothesis of one (H1) vs. no (H0) QTL linked to the set of markers considered. Under hypothesis H1, a QTL effect for each sire and each dam (only dams with more than 10 offspring were taken into account) was fitted to the data. All sufficiently probable (above 0.10) dam phases were considered, so that the likelihood A could not be entirely linearised. For every cM along a linkage group, the likelihood A could be written as:

$$ \Lambda^x = \prod_{i,j} \sum_{hd_{ij}} p(hd_{ij} | M_i, \hat{h}s_i) \prod_{ijk} f(\tilde{y}_{ijk} | \hat{h}s_i, hd_{ij}, M_i) $$

where: $\prod_{i,j}$ is a product over full-sib families of sire i and dam ij, $\sum_{hd_{ij}}$ is a summation over dam phases hd$_{ij}$ with a
<figure><img src="image_2.png" /><figcaption>**Figure 1** Linkage maps of SSC2, SSC4 and SSC6 for the combined pedigree. Microsatellite markers in bold were added for the analysis of the French animals and markers in italics for the Dutch animals; *: on SSC4, microsatellite S0217 is genotyped only on Dutch animals and SW1089 on French animals; on SSC6, SW2466 is genotyped only on Dutch animals and SW607 on French animals</figcaption></figure>

probability greater than 0.10, $\hat{h}_{s_i} = \arg\max_{hs_i} p(hs_i|M_i)$, $p(hs_i|M_i)$ = linkage phase $hs_i$ probability for sire $i$ given marker information $M_i$, $p(hd_{ij}|M_i, \hat{h}_{s_i})$ = linkage phase for dam $ij$ given marker information $M_i$ and sire linkage phase, $f(\tilde{y}_{ijk}|\hat{h}_{s_i}, hd_{ij}, M_i)$ = density function of the adjusted phenotype $\tilde{y}_{ijk}$ of the offspring $ijk$ of the $ij$th dam and the $ith$ sire, conditional on the chromosome segments transmitted by the sire ($q_s$) and the dam ($q_d$). $\tilde{y}_{ijk}$ is supposed to be normally distributed with a mean $\sum_{q_s=1} \sum_{q_d=1} p(d_{ijk}^x = (q_s, q_d)|hs_i, hd_{ij}, M_i)(\mu_i^{xq_s} + \mu_{ij}^{xq_d})$ and a variance

$\sigma_i^2$ estimated within each sire family, where $p(d_{ijk}^x = (q_s, q_d)|hs_i, hd_{ij}, M_i)$ is the transmission probability from parents $i$ and $ij$ to offspring $ijk$, and $\mu_i^{xq_s}$ and $\mu_{ij}^{xq_d}$ can be parameterised as $\mu_{i(j)}^{x1} = \mu_{i(j)} + \alpha_{i(j)}^x / 2$ and $\mu_{i(j)}^{x2} = \mu_{i(j)} - \alpha_{i(j)}^x / 2, \alpha_{i(j)}^x$ and $\alpha_{ij}^x$ being the within-half-sib and within-full-sib average QTL substitution effects and $\mu_{i(j)}$ being the family mean for parent $i(j)$. Average substitution effects, which in the present case are equivalent to additive values ($a$), were hence estimated within each sire family as $\mu_i^{x1} - \mu_i^{x2}$ and within each dam family as $\mu_{ij}^{x1} - \mu_{ij}^{x2}$, and
averaged over families to estimate the average QTL effect in
the population.

To guarantee an accurate estimation of the sire QTL
effects, only sire families with more than 30 progeny
were retained in the analysis, thus 15 sire families were
omitted from the Dutch pedigree. Due to number of
progeny per dam, dam effects were estimated for all
dams in the French pedigree, whereas none was esti-
mated in the Dutch pedigree.

The maximum LRT along the linkage group indicated
the most likely position for a QTL. Significance thresh-
olds were empirically computed using 1000 simulations
under the null hypothesis, assuming an infinitesimal
polygenic model (i.e. the trait is controlled by an infinite
number of additive loci, each with infinitesimal effect
and is thus not influenced by a major QTL) and a nor-
mal distribution of performance traits [21]. In practice,
for progeny p, simulated phenotypes yp were sampled as
the sum of a polygenic part us and an environmental
part ep with normal distributions of mean = 0 and var-
iances depending on the heritability of the trait, the
total phenotypic variance being 1. The polygenic parts
were sampled on the F1 sires and dams (us and ud,
respectively) and the transmission to the progeny was
simulated as up = 0.5 us + 0.5 ud, resulting inyp = 0.5
us + 0.5 ud + ep for progeny p. Simulations were pre-
ferred to permutations because of the family structures
[22]. With these structures, permutations have to be
performed within full-sib families to respect data varia-
bility in the different families. In our case, data number
within each family was not sufficient to achieve an
extensive description of the distribution of the test sta-
tistic under the null hypothesis. QTL were considered
significant if the chromosome-wide significance thresh-
old exceeded 5% and suggestive if it exceeded 10%.
Chromosome-wide significance thresholds were pre-
ferred to genome-wide significance thresholds since only
three chromosomes were included in the analyses. Esti-
mated significant thresholds (at the 5% chromosome-
wide level) varied with traits and pedigree ranging
between 45 and 50 for each independent pedigree and
80 and 85 for the combined pedigree. Confidence inter-
vals around a QTL position were empirically determined
by the “2-LOD drop-off” method [23].

For each chromosomal region, QTL detection analyses
were applied separately on the French and Dutch pedi-
grees, and then on the combination of both pedigrees
thereafter named “combined pedigree”.

Additional analyses were carried out with QTLMap to
test the segregation of two linked QTL in a linkage
group [24] and revealed two situations: (1) when a sig-
nificant QTL had been previously detected (H0 versus
H1), the null hypothesis was “one QTL segregating at
the maximum likelihood position estimated under H1”

and (2) when no QTL had been previously detected (H0
versus H1), the null hypothesis was “no QTL”. In both
cases, the alternative hypothesis (H2) was “two linked
QTL segregating on the linkage group”. The LRT were
computed following a grid-search strategy, using first
5 cM steps along the chromosome to identify significant
regions in which then finer steps (1 cM) were applied.
Significance thresholds were empirically estimated by a
thousand simulations under the null hypothesis as
described by Gilbert and Le Roy [25]. When H0 was
“no QTL”, thresholds were the same as those computed
previously for single QTL tests. When H0 was “one
QTL segregating at the maximum likelihood position
x_max estimated under H1”, simulations were done
assuming that the trait was controlled by a QTL located
at x_max and having the effect estimated for the maxi-
mum likelihood at position x_max in the single QTL
analysis, all F1 being considered as heterozygous for the
QTL. Estimated significant thresholds (at the 5% chro-
mosome-wide level) varied with traits and pedigree, ran-
ging between 85 and 90 for each independent pedigree
and between 140 and 150 for the combined pedigree.

Finally, QTL detection was also carried out on the
adjusted data using a Line-Cross model (LC) with the
online version of QTLExpress [26]. In this report, the
method is only briefly described since results are not
shown in detail. The LC model assumed that Meishan
and European breeds were fixed for alternate QTL
alleles. With the LC model, the adjusted phenotypes
were fitted to a linear model including additive and
dominant components [27] and the chromosome-wide
significance thresholds were determined by permuta-
tions of data as described by Churchill and Doerge [28].

Results

Linkage maps

For chromosomes SCC2, SSC4 and SSC6, the estimated
marker orders of their linkage maps were consistent
with those of the published USDA-MARC linkage maps
http://www.marc.usda.gov/genome/swine/swine.html
and their sex-averaged lengths were 149 cM, 116 cM
and 166 cM, respectively (Figure 1).

QTL detection

Table 2 shows the QTL detection results for each chro-
mosome separately and for each independent pedigree
and the combined pedigree. These results were obtained
by using the half-full sib model with the QTLMap soft-
ware. Additional analyses were done with the Line Cross
model (which assumes that parental breeds are fixed for
alternate QTL alleles) with the QTLExpress online Soft-
ware [26]. The same QTL were detected with the HFS
and LC models (data not shown), except for the QTL
underlying birth weight that was only described with the
<table><caption>Table 2 QTL detected in the two independent and the combined pedigrees by the single QTL detection analysis with significance level <10%</caption><thead><tr><th>Pedigree</th><th>Trait</th><th>SSC</th><th>Position (cM)</th><th>LRT</th><th>Threshold</th><th>Number of segregating sites</th><th>QTL effect</th><th>References</th></tr></thead><tbody><tr><td>French</td><td>BFT</td><td>2</td><td>0 [0-8]</td><td>47.2</td><td>+</td><td>2</td><td>+ 0.25</td><td>[8]</td></tr><tr><td>French</td><td>loin</td><td>2</td><td>0 [0-3]</td><td>58.4</td><td>**</td><td>3</td><td>- 0.30</td><td>[8]</td></tr><tr><td>French</td><td>IMF</td><td>4</td><td>0 [0-12]</td><td>37.0</td><td>+</td><td>3</td><td>- 0.009</td><td>/</td></tr><tr><td>French</td><td>BFT</td><td>4</td><td>61 [55-66]</td><td>45.5</td><td>+</td><td>4</td><td>+ 0.19</td><td>[8]</td></tr><tr><td>French</td><td>LWG</td><td>4</td><td>66 [56-71]</td><td>69.6</td><td>**</td><td>2</td><td>- 0.18</td><td>[6]</td></tr><tr><td>French</td><td>x2</td><td>4</td><td>14 [6-23]</td><td>47.0</td><td>+</td><td>2</td><td>- 0.006</td><td>/</td></tr><tr><td>French</td><td>head weight</td><td>4</td><td>42 [17-49]</td><td>44.6</td><td>+</td><td>3</td><td>- 0.25</td><td>/</td></tr><tr><td>French</td><td>birth weight</td><td>4</td><td>53 [45-58]</td><td>61.3</td><td>**</td><td>4</td><td>- 0.07</td><td>[6]</td></tr><tr><td>French</td><td>belly weight</td><td>4</td><td>55 [48-61]</td><td>56.9</td><td>**</td><td>2</td><td>- 0.04</td><td>[8]</td></tr><tr><td>French</td><td>teat number</td><td>6</td><td>110 [98-130]</td><td>46.0</td><td>+</td><td>3</td><td>+ 0.03</td><td>[15]</td></tr><tr><td>French</td><td>loin weight</td><td>6</td><td>99 [92-108]</td><td>46.6</td><td>+</td><td>3</td><td>+ 0.22</td><td>/</td></tr><tr><td>French</td><td>midriff weight</td><td>6</td><td>144 [137-150]</td><td>59.1</td><td>**</td><td>4</td><td>- 0.15</td><td>/</td></tr><tr><td>Dutch</td><td>weaning weight</td><td>2</td><td>21 [17-24]</td><td>46.2</td><td>+</td><td>5</td><td>+ 0.03</td><td>/</td></tr><tr><td>Dutch</td><td>BFT</td><td>2</td><td>28 [18-46]</td><td>60.4</td><td>**</td><td>9</td><td>+ 0.22</td><td>[11]</td></tr><tr><td>Dutch</td><td>meat</td><td>2</td><td>30 [19-47]</td><td>63.8</td><td>**</td><td>8</td><td>- 0.31</td><td>/</td></tr><tr><td>Dutch</td><td>a*</td><td>2</td><td>26 [19-36]</td><td>51.2</td><td>*</td><td>10</td><td>+ 0.02</td><td>[10]</td></tr><tr><td>Dutch</td><td>IMF</td><td>4</td><td>0 [0-6]</td><td>47.0</td><td>+</td><td>9</td><td>+ 0.14</td><td>[11]</td></tr><tr><td>Dutch</td><td>birth weight</td><td>6</td><td>136 [130-144]</td><td>47.4</td><td>+</td><td>8</td><td>- 0.13</td><td>/</td></tr><tr><td>Dutch</td><td>teat number</td><td>6</td><td>154 [146-161]</td><td>47.5</td><td>+</td><td>8</td><td>- 0.05</td><td>/</td></tr><tr><td>Dutch</td><td>L*</td><td>6</td><td>103 [91-118]</td><td>61.3</td><td>**</td><td>7</td><td>+ 0.20</td><td>/</td></tr><tr><td>Combined</td><td>BFT</td><td>2</td><td>1 [0-5]</td><td>95.6</td><td>*</td><td>12</td><td>+ 0.34</td><td></td></tr><tr><td>Combined</td><td>meat</td><td>2</td><td>32 [25-46]</td><td>97.8</td><td>**</td><td>10</td><td>- 0.32</td><td></td></tr><tr><td>Combined</td><td>IMF</td><td>4</td><td>0 [0-5]</td><td>84.0</td><td>*</td><td>12</td><td>+ 0.11</td><td></td></tr><tr><td>Combined</td><td>teat number</td><td>4</td><td>46 [37-51]</td><td>75.7</td><td>+</td><td>10</td><td>- 0.04</td><td></td></tr><tr><td>Combined</td><td>birth weight</td><td>4</td><td>53 [49-57]</td><td>82.1</td><td>*</td><td>9</td><td>- 0.10</td><td></td></tr><tr><td>Combined</td><td>LWG</td><td>4</td><td>83 [59-90]</td><td>98.6</td><td>**</td><td>11</td><td>- 0.30</td><td></td></tr><tr><td>Combined</td><td>teat number</td><td>6</td><td>104 [94-113]</td><td>84.8</td><td>*</td><td>11</td><td>- 0.01</td><td></td></tr></tbody></table>

Traits indicated in bold are common to the two independent pedigrees; LRT is the maximum Likelihood Ratio Test; the QTL effect is expressed in standard deviation units; the effect is given as Meishan - European alleles; +, *, ** are the 10%, 5% and 1% chromosome-wide significance levels respectively; the references of previously published data are given, /: when no previous result was available for this trait on the chromosome or when this QTL was not detected in the previous analysis.

HFS model in the combined pedigree. Moreover, it is worth noting that with LC model, most QTL were significant at a 1% chromosome-wide threshold whereas with HFS model significant QTL were detected at different thresholds (Table 2).

## SSC2

Single-QTL detection analyses identified a chromosome-wide significant QTL affecting loin weight with the French design and a suggestive QTL affecting back fat thickness. These two QTL are located on the telomeric end of pig chromosome 2 (SCC2p) where the IGF2 gene is located and have already been described [6,8]. With the Dutch design, three chromosome-wide significant QTL affecting meat percentage, BFT and a* colour value were detected around position 25 cM, the latter two QTL already described in [11]. An additional suggestive QTL affecting the weight at weaning was also detected in the same region with the Dutch design. Combining

the French and Dutch pedigrees, only two QTL reached the 5% chromosome-wide significance threshold: significant LRT values were obtained for a QTL influencing back fat thickness at position 1 cM and for a QTL affecting meat percentage at position 32 cM. Using multiple-QTL analyses no additional QTL was detected on SSC2. Additional analyses related to parent-of-origin effect were computed but are not reported in the present study.

## SSC4

With the French design on SSC4, two significant QTL affecting birth weight and belly weight around position 55 cM, and one QTL affecting life weight gain around position 66 cM were detected, all three QTL already described in [6,8]. Two additional suggestive QTL also previously described in [7] were identified affecting intra-muscular fat content at position 0 cM and back fat thickness at 61 cM [8]. Two new suggestive QTL were
detected, affecting X2 at 14 cM, and head weight at 42 cM. Using multiple-QTL detection analyses, a pair of QTL localised at positions 30 cM and 74 cM was detected for teat number. With the Dutch design, only one QTL, previously described in [11] and affecting intra-muscular fat content at 0 cM reached the 10% chromosome-wide suggestive. When combining the French and Dutch pedigrees this QTL reached the 5% chromosome-wide significance threshold (Figure 2). A QTL affecting birth weight around position 55 cM, and a QTL affecting life weight gain at 83 cM, detected only with the French design, were also confirmed in the combined analysis. Using multiple-QTL tests, the hypothesis of two QTL affecting this trait was more likely than a single-QTL hypothesis: the test for loci at 59 and 90 cM reached the 5% chromosome-wide significance threshold (Figure 3). Additionally, a suggestive QTL influencing teat number was detected in the combined analysis at 46 cM using the single-QTL analysis. The two-QTL model retained in the analysis of the French pedigree for this trait was not significant with the combined analysis.

SSC6

With the French design on SSC6, a significant QTL affecting midriff weight at 144 cM and two other suggestive QTL one influencing loin weight (99 cM) and one affecting teat number (110 cM) were detected. With the Dutch design, a significant QTL influencing L* parameter at 103 cM and two other suggestive QTL one affecting birth weight (136 cM) and one affecting teat number (154 cM) were identified. For this last trait, a two-QTL model was significant at the 5% chromosome-wide significance threshold for two loci localized at 50 and 155 cM. When combining the French and Dutch pedigrees only one significant QTL affecting teat number at 104 cM was detected at the 5% chromosome-wide threshold.

Discussion

The aim of this study was to combine two F2 designs produced independently in France and in the Netherlands to detect QTL influencing economically important traits. These two designs were selected on the following criteria: (1) comparable founder breeds (Meishan and

<figure><img src="image_3.png" /><figcaption>**Figure 2** QTL underlying Intra Muscular Fat content on SSC4 for the three studied pedigrees. The solid line gives the result for the combined pedigree, the circled line for the French pedigree and the crossed line for the Dutch pedigree; for each analysis, the LRT is presented in proportion to the 5% threshold on the chromosome.</figcaption></figure>
<figure><img src="image_4.png" /><figcaption>**Figure 3** Two-QTL analysis results for Life Weight Gain on SSC4 with the combined pedigree. The z axis gives the value of the LRT divided by the 5% threshold at the chromosome-wide level, the surface shown in gray corresponding to a ratio higher than 1.</figcaption></figure>

Large White and/or Landrace European breeds) and (2) the quasi-equal number of offspring produced in both protocols. Furthermore, although the European breeds were not identical, the Meishan sires used in the Dutch pedigree are related to the French Meishan dams. This supports the assumption that common Meishan QTL alleles segregated in both designs. European and Meishan breeds should have contrasted haplotypes (haplotypes being highly similar in both Meishan pedigrees) and QTL should segregate for the same loci. However, the two populations differ with respect to the number of families (six F1 sire families in the French design versus 39 F1 sires families in the Dutch design) and the reciprocal cross used to produce the F1 animals. To combine these designs the six French families and 24 of the 39 Dutch families, composed of more than 30 offspring, were retained. Our study focused on three pig chromosomes, SSC2, SSC4 and SSC6, for which QTL had already been detected. Despite the lack of overlap between some of the identified QTL from the pedigrees

analysed separately, joint analyses of the combined pedigree should provide (1) greater significance of shared QTL, (2) exclusion of false suggestive QTL and (3) greater mapping precision for shared QTL. First, we investigated how the addition of new genotypes influenced the two designs. QTL detections were performed independently for each pedigree and for all the traits to be compared to the results previously published. We were then able to estimate the advantage between a combined analysis and independent analyses.

## QTL detected with the French design

In the analysis with the French design, the results were consistent with those previously reported although some differences were observed. In 2001, Bidanel et al. [6] have reported a highly significant QTL underlying BFT on SSC4 between two markers located at positions 43 and 83 cM. Addition of microsatellite SW35 at 52.7 cM in the present study resulted in a loss of significance for this trait (10% chromosome-wide threshold). On
the opposite, we detected in the same region a sugges-
tive QTL affecting head weight which had not been pre-
viously identified in this pedigree [8]. This shows that
the addition of a single highly informative marker in a
region with a low marker density can change results to
a great extent. Cepica et al. [29] have reported a QTL
affecting head weight in the same region on this chro-
mosome. On SSC6, Milan et al. [8] have described a
suggestive QTL for belly weight between positions 2
and 32 cM In our analysis, this QTL was neither signifi-
cant, nor suggestive. In this case also, a microsatellite,
SW2535, was also added above SW2406, which allowed
a better coverage of the telomeric part of the chromo-
some. Due to the addition of this marker, we can con-
clude that the previously suggested QTL is probably a
false-positive. Finally, with the French pedigree, three
new QTL on SSC6: two QTL, one affecting loin weight
(99 cM) and one affecting teat number (110 cM), and
one significant QTL influencing midriff weight (144
cM). Previous studies had revealed QTL affecting loin
weight on SSC6 at 83 cM in crosses involving Pietrain,
Large-White and Leicoma animals [30], and between
122 and 149 cM with a Duroc x Pietrain design [31].

QTL detected with the Dutch design

With the Dutch pedigree alone, it was possible that our
results differed from those previously reported because
of (1) the addition of new markers (as for the French
pedigree), (2) the selection of 24 families among the 38
that were used by de Koning et al. [11] and (3) the
model used (mixture of full and half-sib families vs. line
cross model). Initially, de Koning et al. have detected
many QTL regions using a line-cross model and our
results are closer to those obtained using a half-sib
model [10-12]. Among the previously QTL detected, we
did not confirm two suggestive QTL influencing the b*
colour value [10] and LWG on SSC6 [12]. These differ-
ences could be due to the addition of new markers
showing that these QTL are false-positive ones. How-
ever, we cannot exclude that they were lost because of
the reduction of the number of families. It is possible
that these QTL were segregating in some of the 14
excluded families, which results in loss of power. How-
ever, with the addition of new markers in the Dutch
pedigree, we detected a new QTL underlying the L* col-
our value on SSC6 at 103 cM. No QTL influencing this
meat quality trait has been reported in this region. In a
Landrace x Iberian design, Ovilo et al. [32] have
described one QTL at 40 cM on SSC6 influencing meat
lightness in interaction with another locus on SSC18.
We also detected a QTL underlying birth weight at 136
cM while Yue et al. [33] had reported a QTL for the
same trait in a cross between Meishan and Pietrain
breeds at 98.2 cM.

Our results illustrate the influence of marker density
to exclude false-positive QTL and to detect new ones.
In 2008, Hu and Xu [34] have demonstrated using simu-
lations that when the interval between two adjacent
markers decreases, the power of QTL detection
increases. Here we confirm this statement with real
datasets.

QTL detected combining the French and Dutch designs

The size of a QTL design has a major influence on the
power of QTL detection. A QTL can be detected thanks
to the recombination events occurring during the
gamete production in F1 animals. In an F2 protocol, the
number of recombination events is limited. Therefore,
to confirm or refine the position of QTL initially
detected with an F2 protocol, the number of crossing-
over events cannot be increased without additional ani-
mals. Producing additional animals is expensive and
time-consuming. Thus, a joint analysis of independent
pedigrees is an easy alternative to enhance the number
of F2 animals in the design. One suggestive QTL affect-
ing intramuscular fat content on SSC4 had been
detected in both the French and Dutch pedigrees inde-
pendently. By combining the two designs, this QTL
reached a 5% chromosome-wide significant threshold.
As shown on Figure 2, the two pedigrees contributed
equally to the LRT. The proportion of segregating
fathers per pedigree and/or the QTL effect were prob-
ably too small to detect this QTL in each pedigree inde-
pendently but the combined pedigree resulted in more
power and the QTL reached the significance threshold.
For SSC2, independent analyses of each pedigree indi-
cated that a single QTL influencing BFT and meat per-
centage was segregating on the telomeric end of the
short arm of this chromosome but not at the same posi-
tion. With the French pedigree, LRT value was maximal
in the IGF2 region (0 cM) whereas with the Dutch
design it was around 28 cM. Analysis of the combined
pedigree confirmed that in fact two linked loci are seg-
regating in this region, i.e. the IGF2 gene and a locus
located around 30 cM (Figure 4). QTL influencing meat
percentage have been previously detected between 40
and 60 cM with crosses involving either Meishan, Pie-
train and Wild Boar breeds [35], or Large White and
Meishan animals [36].

Concerning the two-QTL detection analysis, our study
provides evidence that combining two pedigrees and
adding new markers increases the power of QTL detec-
tion. With a single-QTL analysis, we can conclude that
there is a single QTL at 80 cM on SSC4 influencing life
weight gain whereas in fact this maximum statistic is
probably due to the combined actions of two QTL. The
combination of two designs indicates that the presence
of two different QTL located at 60 and 90 cM and
<figure><img src="image_5.png" /></figure>

common to both populations is more likely than the segregation of a unique QTL at 80 cM (Figure 3).

## Limits of the combined analysis

The benefit of combined analyses is essentially obtained when the QTL is segregating in both pedigrees. If QTL were segregating in only one of the two designs, detection power and estimated QTL effect could be reduced. We observed this situation for some QTL initially detected independently in one of both pedigrees. For example, a significant QTL affecting birth weight was detected only with the French design at 53 cM on SSC4 and was confirmed with the combined pedigree but it was not a major segregating QTL in the Dutch pedigree, the significance of this LRT was reduced in the combined analysis. Thus, this QTL is specific to the French pedigree. At the extreme, a QTL previously detected in one pedigree may be undetectable in the combined analysis. This situation was observed only for suggestive QTL such as that affecting birth weight on SSC6 detected in the Dutch pedigree only. Thus these QTL are also specific to the pedigree in which they were detected. The same situation has been reported by Walling et al. [4] who have shown that among seven different pedigrees a QTL affecting birth weight detected on SSC4 segregated in the French pedigree only but after

adding the six other pedigrees, although the QTL was still detected, its significance was lower and its effect was divided by two.

## Influence of the breeds used

QTL segregating in several independent designs will be largely influenced by the breeds used. On the one hand, by combining two pedigrees (involving only three breeds), Kim et al. in 2005 [37] have detected 10 new QTL undetected in either of the two pedigrees. In this case, the combination of pedigrees increased the number of interesting regions. On the other hand, Pérez-Enciso et al. [5] in 2005 have demonstrated that by combining pedigrees, the possibility of detecting new QTL is sometimes reduced. This analysis was performed using five independent crosses involving six pig breeds (Iberian, Landrace, Large White, Meishan, Pietrain and Wild Boar). If a QTL allele is fixed in one breed involved in one cross only, the addition of other pedigrees that do not involve these breeds can reduce the effect of this QTL and therefore make it less detectable. This is also supported by the report of Guo et al. [38] in 2008 who have shown that if a QTL is population-dependant it is highly probable that combining pedigrees will provide no benefit. To avoid this drawback of a joint analysis, Guo et al. have combined two Meishan
x Large White pedigrees with F0 animals coming from the same populations. Combining pedigrees seems more advantageous if pedigrees involve similar breeds like in our combined study. In our case, the breeds involved in both pedigrees were genetically very similar, which led to the detection of common QTL (BFT on SSC2 or IMF on SSC4). In some cases, the QTL detected with the combined pedigree had only been described in one of the two designs like birth weight on SSC4. This result can be due to either a population specific QTL, or to differences in the proportion of families contributing to the statistic signal.

## Reduction of QTL interval

An important advantage of a combined analysis is to estimate more precisely mapping intervals. This parameter may be largely influenced by the total number of markers used and the number of common markers analysed on the different designs. In our study, new data were generated by genotyping animals of both pedigrees with additional microsatellite markers to maximise the power and the precision. Previous studies that combined two or more pedigrees only analysed pre-existing datasets. The originality of our study is the large number of common markers, i.e. 29 common markers evenly spaced along the three chromosomes. In 2003, Bennewitz et al. [39] have combined two bovine granddaugh ter designs and surprisingly, obtained larger confidence intervals than with single designs. Previously, using simulations they had estimated that an increase of population size led to a reduction of confidence intervals [40,41]. Thus, they have proposed that the increase in confidence intervals is probably due to the fact that the families were genotyped with different sets of markers and only a few belonging to both sets. Using a dataset of highly common markers makes it possible to increase the confidence interval of a QTL. In our analysis, we avoided this drawback by analysing a unique set of microsatellite markers on both pedigrees and on the combined one. With a single set of markers, we obtained narrower confidence intervals than previously reported estimated intervals calculated with the same drop-off method [6,8]. In the present study, confidence intervals varied from 3 to 30 cM, with an average size of 15 cM. In the initial analyses, the confidence intervals spanned on average 40 cM. Unfortunately, since in the first analyses reported by de Koning et al. [11,12] no confidence intervals were calculated, we could not carry out a comparison.

## Conclusion

Combining two Meishan x European breeds F2 pedigrees improved the mapping of QTL compared to separate analyses of the pedigrees. We detected new QTL,

confirmed some QTL which were previously described and excluded false positive QTL. Our work was facilitated by the access to raw phenotypic data and DNA of animals from both pedigrees. Analysis of a single set of markers proved more efficient to obtain smaller QTL intervals. The combination of two designs involving very similar breeds with the addition of new markers allowed us to fine-map QTL without phenotyping additional animals.

## Acknowledgements

The authors would like to thank Hélène Hayes for her precious help for the English revision of the manuscript and the reviewers for their useful comments.

## Author details

¹INRA, UMR 0444 Génétique Cellulaire, F-31326 Castanet-Tolosan, France.

²INRA, UMR 1313 Génétique Animale et Biologie Intégrative, F-78352 Jouyen-Josas, France. ³Wageningen University, Animal Breeding and Genetics Group, 6700 AH Wageningen, The Netherlands.

## Authors' contributions

FT carried out the genotyping and the QTL detection analyses. HG participated in the QTL detection analyses. HCMH provided the information and the phenotypic data sets of the Dutch pedigree. JPB provided the information and the phenotypic datasets of the French pedigree. MAMG and JR conceived of the study. JR participated in its design and coordination and helped to draft the manuscript. All authors contributed and approved the final manuscript.

## Competing interests

The authors declare that they have no competing interests.

Received: 4 May 2010 Accepted: 25 November 2010
Published: 25 November 2010

## References

1. Bidanel JP, Rothschild MF: Current status of quantitative trait locus mapping in pigs. *Pig News Inf* 2002, **23**:39N-53N.

2. Hu ZL, Recey JM: Animal QTLdb: beyond a repository. A public platform for QTL comparisons and integration with diverse types of structural genomic information. *Mamm Genome* 2007, **18**:1-4.

3. Allison DB, Heo M: Meta-analysis of linkage data under worst-case conditions: a demonstration using the human OB region. *Genetics* 1998, **148**:859-865.

4. Walling GA, Visscher PM, Andersson L, Rothschild MF, Wang L, Moser G, Groenen MA, Bidanel JP, Cepica S, Archibald AL, Geldermann H, de Koning DJ, Milan D, Haley CS: Combined analyses of data from quantitative trait loci mapping studies. *Chromosome 4* effects on porcine growth and fatness. *Genetics* 2000, **155**:1369-1378.

5. Perez-Enciso M, Mercade A, Bidanel JP, Geldermann H, Cepica S, Bartenschlager H, Varona L, Milan D, Folch JM: Large-scale, multibreed, multitrait analyses of quantitative trait loci experiments: the case of porcine X chromosome. *J Anim Sci* 2005, **83**:2289-2296.

6. Bidanel JP, Milan D, Iannuccelli N, Amigues Y, Boscher MY, Bourgeois F, Caritez JC, Gruand J, Le Roy P, Lagant H, Quintanilla R, Renard C, Gellin J, Ollivier L, Chevalet C: Detection of quantitative trait loci for growth and fatness in pigs. *Genet Sel Evol* 2001, **33**:289-309.

7. Larzul C, Bidanel JP, Iannuccelli N, Gruand J, Mourot J, Milan D: Détection de QTL pour la teneur en lipides intramusculaires et les activités enzymatiques lipogéniques dans un croisement entre porcs Meishan et Large White. *Journées Recherche Porcine* 2008, **40**:127-130.

8. Milan D, Bidanel JP, Iannuccelli N, Riquet J, Amigues Y, Gruand J, Le Roy P, Renard C, Chevalet C: Detection of quantitative trait loci for carcass composition traits in pigs. *Genet Sel Evol* 2002, **34**:705-728.

9. Hirooka H, de Koning DJ, Harlizius B, van Arendonk JA, Rattink AP, Groenen MA, Brascamp EW, Bovenhuis H: A whole-genome scan for
quantitative trait loci affecting teat number in pigs. *J Anim Sci* 2001, 79:2320-2326.

10. de Koning DJ, Harlizius B, Rattink AP, Groenen MA, Brascamp EW, van Arendonk JA: Detection and characterization of quantitative trait loci for meat quality traits in pigs. *J Anim Sci* 2001, 79:2812-2819.

11. de Koning DJ, Janss LL, Rattink AP, van Oers PA, de Vries BJ, Groenen MA, van der Poel JJ, de Groot PN, Brascamp EW, van Arendonk JA: Detection of quantitative trait loci for backfat thickness and intramuscular fat content in pigs (Sus scrofa). *Genetics* 1999, **152**:1679-1690.

12. de Koning DJ, Rattink AP, Harlizius B, Groenen MA, Brascamp EW, van Arendonk JA: Detection of characterization of quantitative trait loci for growth and reproduction traits in pigs. *Livest Prod Sci* 2001, **72**:185-198.

13. Janss LL, van Arendonk JA, Brascamp EW: Bayesian statistical analyses for presence of single genes affecting meat quality traits in a crossed pig population. *Genetics* 1997, **145**:395-408.

14. Janss LL, Van Arendonk JA, Brascamp EW: Segregation analyses for presence of major genes affecting growth, backfat, and litter size in Dutch Meishan crossbreds. *J Anim Sci* 1997, **75**:2864-2876.

15. Bidanel JP, Rosendo A, Iannuccelli N, Riquet J, Gilbert H, Caritez JC, Billon Y, Amigues Y, Prunier A, Milan D: Detection of quantitative trait loci for teat number and female reproductive traits in Meishan X Large White F2 pigs. *Animal* 2008, **2**:813-820.

16. Daumas G, Dhorne T: Nouvelles équations françaises de prédiction du taux de muscle des carcasses de porc. *Journées Recherche Porcine* 1994, **26**:151-156.

17. Iannuccelli E, Woloszyn N, Arhainx J, Gellin J, Milan D: GEMMA: a database to manage and automate microsatellite genotyping. *Anim Genet* 1996, **27**(Suppl 2):55.

18. Green P, Falls K, Crooks S: Documentation for CRIMAP version 2.4. Washington University School of Medicine, St Louis; 1990.

19. Mangin B, Goffinet B, Le Roy P, Boichard D, Elsen JM: Alternative models for QTL detection in livestock. II. Likelihood approximations and sire marker genotype estimations. *Genet Sel Evol* 1999, **31**:225-237.

20. Elsen JM, Mangin B, Goffinet B, Boichard D, Le Roy P: Alternative models for QTL detection in livestock. I. General introduction. *Genet Sel Evol* 1999, **31**:213-224.

21. Le Roy P, Elsen JM, Boichard D, Mangin B, Bidanel JP, Goffinet B: An algorithm for QTL detection in mixture of full and half sib families. *Proceedings of the 6th World Congress on Genetics Applied to Livestock Production* 1998, **26**:257-260.

22. Churchill GA, Doerge RW: Naive application of permutation testing leads to inflated type I error rates. *Genetics* 2008, **178**:609-610.

23. Lander E, Kruglyak L: Genetic dissection of complex traits: guidelines for interpreting and reporting linkage results. *Nat Genet* 1995, **11**:241-247.

24. Gilbert H, Le Roy P: Power of three multitrait methods for QTL detection in crossbred populations. *Genet Sel Evol* 2004, **36**:347-361.

25. Gilbert H, Le Roy P: Methods for the detection of multiple linked QTL applied to a mixture of full and half sib families. *Genet Sel Evol* 2007, **39**:139-158.

26. Seaton G, Haley CS, Knott SA, Kearsey M, Visscher PM: QTL Express: mapping quantitative trait loci in simple and complex pedigrees. *Bioinformatics* 2002, **18**:339-340.

27. Haley CS, Knott SA, Elsen JM: Mapping quantitative trait loci in crosses between outbred lines using least squares. *Genetics* 1994, **136**:1195-1207.

28. Churchill GA, Doerge RW: Empirical threshold values for quantitative trait mapping. *Genetics* 1994, **138**:963-971.

29. Cepica S, Stratil A, Kopecny M, Blazkova P, Schröffel J Jr, Davoli R, Fontanesi L, Reiner G, Bartenschlager H, Moser G: Linkage and QTL mapping for Sus scrofa chromosome 4. *J Anim Breed Genet* 2003, **120**(suppl 1):28-37.

30. Mohrmann M, Roehe R, Knap PW, Looft H, Plastow GS, Kalm E: Quantitative trait loci associated with AutoFOM grading characteristics, carcass cuts and chemical body composition during growth of Sus scrofa. *Anim Genet* 2006, **37**:435-443.

31. Edwards DB, Ernst CW, Raney NE, Doumit ME, Hoge MD, Bates RO: Quantitative trait locus mapping in an F2 Duroc x Pietrain resource population: II. Carcass and meat quality traits. *J Anim Sci* 2008, **86**:254-266.

32. Ovilo C, Clop A, Noguera JL, Oliver MA, Barragan C, Rodriguez C, Silio L, Toro MA, Coll A, Folch JM, Sanchez A, Babot D, Varona L, Perez-Enciso M:

Quantitative trait locus mapping for meat quality traits in an Iberian x Landrace F2 pig population. *J Anim Sci* 2002, **80**:2801-2808.

33. Yue G, Stratil A, Kopecny M, Schröffelova D, Schröffel J Jr, Hojny J, Cepica S, Davoli R, Zambonelli P, Brunsch C, Sternstein I, Moser G, Bartenschlager H, Reiner G, Geldermann H: Linkage and QTL mapping for Sus Scrofa chromosome 6. *J Anim Breed Genet* 2003, **120**(suppl 1):45-55.

34. Hu Z, Xu S: A simple method for calculating the statistical power for detecting a QTL located in a marker interval. *Heredity* 2008, **101**:48-52.

35. Lee SS, Chen Y, Moran C, Cepica S, Reiner G, Bartenschlager H, Moser G, Geldermann H: Linkage and QTL mapping for Sus scrofa chromosome 2. *J Anim Breed Genet* 2003, **120**(suppl 1):11-19.

36. Qu YC, Deng CY, Xiong YZ, Zheng R, Yu L, Su YH, Liu GL: The construction of the genetic map and QTL locating analysis on chromosome 2 in swine. *Yi Chuan Xue Bao* 2002, **29**:972-976.

37. Kim JJ, Rothschild MF, Beever J, Rodriguez-Zas S, Dekkers JC: Joint analysis of two breed cross populations in pigs to improve detection and characterization of quantitative trait loci. *J Anim Sci* 2005, **83**:1229-1240.

38. Guo YM, Lee GJ, Archibald AL, Haley CS: Quantitative trait loci for production traits in pigs: a combined analysis of two Meishan x Large White populations. *Anim Genet* 2008, **39**:486-495.

39. Bennewitz J, Reinsch N, Grohs C, Leveziel H, Malafosse A, Thomsen H, Xu N, Looft C, Kuhn C, Brockmann GA, Schwerin M, Weimann C, Hiendleder S, Erhardt G, Medjugorac I, Russ I, Forster M, Brenig B, Reinhardt F, Reents R, Averdunk G, Blumel J, Boichard D, Kalm E: Combined analysis of data from two granddaughter designs: A simple strategy for QTL confirmation and increasing experimental power in dairy cattle. *Genet Sel Evol* 2003, **35**:319-338.

40. Bennewitz J, Reinsch N, Kalm E: Improved confidence intervals in quantitative trait loci mapping by permutation bootstrapping. *Genetics* 2002, **160**:1673-1686.

41. Visscher PM, Thompson R, Haley CS: Confidence intervals in QTL mapping by bootstrapping. *Genetics* 1996, **143**:1013-1020.

doi:10.1186/1297-9686-42-42

Cite this article as: Tortereau *et al.*: Combining two Meishan F2 crosses improves the detection of QTL on pig chromosomes 2, 4 and 6. *Genetics Selection Evolution* 2010 **42**:42.

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