# Number and mode of inheritance of QTL influencing backfat thickness on SSC2p in Sino-European pig pedigrees

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

► To cite this version:

Flavie Tortereau, Hélène Gilbert, Henri Cm H. C. Heuven, Jean Pierre Bidanel, Martien Am M. A. Groenen, et al.. Number and mode of inheritance of QTL influencing backfat thickness on SSC2p in Sino-European pig pedigrees. Genetics Selection Evolution, 2011, 43 (march), Non paginé. 10.1186/1297-9686-43-11. hal-01000887

HAL Id: hal-01000887

https://hal.science/hal-01000887

Submitted on 29 May 2020

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published 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
Number and mode of inheritance of QTL
influencing backfat thickness on SSC2p
in Sino-European pig pedigrees

Flavie Tortereau<fnref n="1" />,³*, Hélène Gilbert², Henri CM Heuven³, Jean-Pierre Bidanel², Martien AM Groenen³ and Juliette Riquet<fnref n="1" />

**Abstract**

**Background:** In the pig, multiple QTL associated with growth and fatness traits have been mapped to chromosome 2 (SSC2) and among these, at least one shows paternal expression due to the IGF2-intron3-G3072A substitution. Previously published results on the position and imprinting status of this QTL disagree between analyses from French and Dutch F2 crossbred pig populations obtained with the same breeds (Meishan crossed with Large White or Landrace).

**Methods:** To study the role of paternal and maternal alleles at the IGF2 locus and to test the hypothesis of a second QTL affecting backfat thickness on the short arm of SSC2 (SSC2p), a QTL mapping analysis was carried out on a combined pedigree including both the French and Dutch F2 populations, on the progeny of F1 males that were heterozygous (A/G) and homozygous (G/G) at the IGF2 locus. Simulations were performed to clarify the relations between the two QTL and to understand to what extent they can explain the discrepancies previously reported.

**Results:** The QTL analyses showed the segregation of at least two QTL on chromosome 2 in both pedigrees, i.e. the IGF2 locus and a second QTL segregating at least in the G/G F1 males and located between positions 30 and 51 cM. Statistical analyses highlighted that the maternally inherited allele at the IGF2 locus had a significant effect but simulation studies showed that this is probably a spurious effect due to the segregation of the second QTL.

**Conclusions:** Our results show that two QTL on SSC2p affect backfat thickness. Differences in the pedigree structures and in the number of heterozygous females at the IGF2 locus result in different imprinting statuses in the two pedigrees studied. The spurious effect observed when a maternally allele is present at the IGF2 locus, is in fact due to the presence of a second closely located QTL. This work confirms that pig chromosome 2 is a major region associated with fattening traits.

**Introduction**

Many QTL associated with economically important
traits like growth, fatness and meat quality have been
detected since the 2000 s, as reviewed by Bidanel and
Rotschild in 2002 [1]. However, even for those that have
been fine-mapped, successful identification of the causal
mutation is rare. In 1999, a paternally expressed QTL
affecting backfat thickness (BFT) and muscle mass was
identified on the short arm of SSC2 close to the IGF2

gene in crosses between Large White (LW) and
European Wild Boar [2] and between LW and Pietrain
[3]. In 2003, Van Laere et al. [4] reported that the IGF2-
intron3-G3072A substitution is the causal mutation.
This mutation affects the binding site of a repressor and
up-regulates IGF2 expression in skeletal muscles and
heart, inducing major maternally imprinted effects on
muscle growth, heart size and fat deposition. Therefore,
selection for animals carrying allele A at this locus is a
major issue in pig production. Analysis of the frequency
and effects of this mutation in pig populations of differ-
ent genetic origins showed that both wild (G) and

<footnote>* Correspondence: flavie.tortereau@toulouse.inra.fr</footnote>

<footnote n="1">¹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>
mutant (A) alleles still segregate in modern populations (LW × Pietrain cross) [5,6], and that allele A is very rare or even nonexistent in local breeds and Wild Boars [7]. The strong favourable effect of allele A was confirmed in both Spanish [8] and Polish [9] LW and Landrace (LR) breeds. In 2004, Jungerius et al. [5] demonstrated that the mutation also explains the major imprinted QTL for backfat thickness in a cross between Meishan (MS) and European White pigs (LW and LR). Yet, although significant effects of the IGF2 mutation were revealed both by ultrasonic and carcass BFT measurements, the presence of a second QTL at a position near 40 cM, as previously described in this population by de Koning [10], cannot be excluded [5]. In the French LW × MS cross, QTL affecting loin weight and BFT on carcass have also been detected near the IGF2 locus [11]. However, surprisingly, no imprinting effect could be detected [12], although the breeds involved are similar (European White breeds and MS) in the Dutch and French studies, and the MS animals in both crosses are related. It has been shown that spurious imprinting effects can exist because of maternal effects [13] or because of linkage disequilibrium [14]. The aims of the present work were to estimate more precisely the IGF2 substitution effect by combining the two MS × European intercrosses, and to investigate further the genetic determinism of the SSC2p chromosomal region by testing the hypothesis of an additional QTL segregating on SSC2 in these populations. In addition, simulation studies were conducted to investigate how the presence of two QTL could affect the apparent mode of inheritance of IGF2 alleles.

## Methods

### Animals and phenotypic data

The French and Dutch F2 MS × European breeds crosses and the recorded phenotypes have been described previously [15-17]. Briefly, the French INRA-PORQTL pedigree consisted of 12 F0 (six LW sires and six MS dams), 26 F1 (six sires and 20 dams) and 521 castrated male F2 pigs. All animals were born and raised at the INRA GEPA experimental research unit (Surgères, Charentes). The Dutch pedigree, obtained from the University of Wageningen (WU), was initiated by mating 19 MS sires to 126 LW and LR dams, resulting in an F1 population of 39 sires and 265 F1 dams, which produced a total of 1212 F2 offspring. The Dutch pedigree was bred in five different breeding companies. Among the 39 Dutch half-sib families, only the 24 largest (more than 30 progeny) were retained in the present analysis in order to homogenize the family structure of the two pedigrees.

Among the traits recorded in the two populations, BFT measured between the third and the fourth rib of

carcass at 6 cm from the spine [10,11] was considered here as the main common trait shared in both designs affected by the QTL under study. This trait was recorded on 565 Dutch pigs (castrated males and females) and on 521 French pigs (castrated males only).

Phenotypic data were first adjusted for fixed effects and covariates with the GLM procedure in SAS® (SAS® 9.1, SAS® Institute, Inc.). The models used to adjust the data included the effects of batch, slaughter day and carcass weight for the INRA pedigree and breeding company, sex, slaughter day and carcass weight for the Dutch pedigree.

### Genetic data

Animals from both pedigrees were genotyped for 11 microsatellites evenly spaced on chromosome 2 (SW2443 (0 cM); SWC9 (2 cM); SW2623 (11 cM); SW256 (23 cM); S0141 (37 cM); SW240 (51 cM); S0091 (76 cM); S0010 (90 cM); S0368 (96 cM); S0378 (108 cM) and S0036 (149 cM)), as previously reported [18].

Genotyping of the IGF2-intron3-G3072A substitution was performed on some of the F0 and F1 animals of both pedigrees. Previously, F1 boars and their parents [5] from the Dutch pedigree had been genotyped by the pyrosequencing technique (Pyrosequencing AB) described in [4]. In the French design, all F0 and F1 animals were genotyped by PCR-RFLP using primers 5'-GGACCGAGCCAGGGACGAGCCT-3' and 5'-GGGA GGTCCCAGAAAAAGTC-3'. The polymerase chain reaction was carried out at 57°C using the GC-RICH PCR system (ROCHE), in presence of 1 M GC-RICH Resolution solution, and 1.5 mM of MgCl2. PCR-RFLP with the restriction enzyme ApeK1 was used to detect the mutation according to the manufacturer's recommendations for time, temperature and buffer conditions. Then, genotypes of all F2 animals at the mutation were inferred for non recombinant haplotypes inherited from F1 individuals, using information from the pedigree and from the transmission of parental haplotypes for surrounding markers (SW2443 and SWC9). No genotype was assigned for recombinant F2 piglets with a heterozygous A/G parent or if the mother had not been genotyped for the mutation. The parental origin of the allele inherited at the A/G substitution was also inferred when possible according to the phase they inherited from their parents.

### QTL analyses

QTL detection was performed on the adjusted data using the QTLMap software [19,20] as explained in [18]. Parameter estimates were obtained by maximization of the likelihood with a Newton-Raphson algorithm, and a Likelihood Ratio Test (LRT) was computed at each cM along SSC2. The maximum LRT along the linkage
group indicated the most likely position for a QTL. For each sire, the substitution effect corresponds to the difference between the Meishan and the European alleles, a positive effect indicating an increased value of the trait due to Meishan alleles. The average QTL substitution effect was computed as the mean of the absolute values of the sire substitution effects. QTL significance thresholds were empirically computed using 1000 simulations under the null hypothesis, assuming an infinitesimal polygenic model for the trait, as described by Gilbert and Le Roy [21].

QTL detection analyses were carried out first for the French and Dutch pedigrees separately, and then for the combined pedigree. A potential second QTL segregating within these pedigrees was investigated with two different methods. First, the multi-QTL option of QTLMAP was used to detect two linked QTL on SSC2 for BFT. The alternative hypothesis (H1) of two QTL segregating was compared to the null hypothesis of one QTL segregating at the IGF2 locus. The LRT were computed following a grid-search strategy, using 5 cM steps along the chromosome. Significance thresholds were empirically estimated by 1000 simulations under the null hypothesis, as described by Gilbert and Le Roy [21]. In a second analysis, the segregation of a potential additional QTL was investigated: (1) by analysing the data from the progeny of sires homozygous at the IGF2 locus (G/ G) and (2) by performing a QTL detection analysis on the full combined pedigree with a model that included IGF2 as a fixed effect.

## Mode of inheritance of the QTL

Analyses of variance (ANOVA) were carried out to infer the inheritance pattern of the SSC2 QTL, using data adjusted for the previously described fixed effects. Tests were applied to compare different effects $\alpha_i$ in the model $Y_{ij} = \mu + \alpha_i + \varepsilon_{ij}$, where $Y_{ij}$ is the adjusted performance of individual $j$ of genotype $i$ (see below), $\mu$ is the population mean, $\varepsilon_{ij}$ is the residual error of individual $j$ of genotype $i$, and $\alpha_i$ is the tested effect. Three different effects for $\alpha_i$ were built based on the following inheritance patterns:

- Only the paternally inherited allele at the mutation has an effect (model IGF2pat, $i = \{A,G\}$)

- Only the maternally inherited allele at the mutation has an effect (model IGF2mat, $i = \{A,G\}$)

- Both the paternally inherited allele and the maternally inherited allele at the mutation have effects (model IGF2patmat, $i = \{AA,AG,GA,GG\}$, the paternal allele being written first).

These analyses of variance were applied to all F2 individuals of the combined pedigree, of both pedigrees

separately, and to sub-groups of animals defined according to the genotype of the parents at the IGF2 mutation:

- F2 having A/G sires

- F2 having A/G dams

- F2 having A/G sires and G/G dams

- F2 having G/G sires and A/G dams

- F2 having A/G sires and A/G dams.

## Detection of spurious effects of the maternally inherited IGF2-allele

Simulation studies were performed with the QTLMAP software to evaluate the power of the inheritance pattern analyses and of the additional QTL studies proposed in this paper, in the presence of a major imprinted gene in the chromosomal region investigated. A QTL segregating at 44 cM affecting the trait was simulated while assuming a paternal effect of the IGF2 locus. Only the phenotypes were simulated; family structures and genotypes were obtained from real data from the two pedigrees. The effect of the IGF2pat model was set to 0.48 phenotypic standard deviations of the trait (as estimated in the data set, see below). The QTL was assumed to be bi-allelic, with the Q allele decreasing backfat level as compared to the q allele. The simulated QTL effect represented the substitution effect of allele q by allele Q. Simulations were then performed with the following parameter values:

- Frequency of the QTL alleles:

- in the F0 European breeds (French F0 males and Dutch F0 females) for allele Q: 0.25, 0.50, 0.75 or 1.00

- in the Meishan populations (French F0 females and Dutch F0 males) for allele q: 1.00. Fixation was assumed based on the small size of the original population and based on the fact that MS individuals were also homozygous for IGF2.

- Effect of the simulated QTL: 0.22 or 0.32 or 0.42 phenotypic standard deviations of the trait.

For each simulation, a QTL analysis was performed as described above and the value of the maximum LRT (LRTmax) and its position were recorded. Simulated phenotypes were exported to the SAS® software and analyses of variance were performed as previously described to determine which inheritance pattern was detected depending on the simulated parameters, applying either the IGF2pat or the IGF2mat models. For the analyses of the two pedigrees separately, 2000 replicate simulations were performed for each combination of frequency × effect parameters. For the combined pedigree, 2000 replicates were performed with an effect of
0.32 and a frequency of 0.50 for the QTL in both grand-
parental populations, as the two pedigrees were recipro-
cal. The percentage of replicates returning significant
results for each pattern of inheritance of IGF2 and
detection of the QTL were then computed from the
2000 replicates for each situation.

Results

Genotyping results for the IGF2 mutation

The IGF2-intron3-G3072A mutation was genotyped for
most of the F0 and F1 founders of both pedigrees
(Table 1). Presence of IGF2 wild type and mutant alleles
in the Dutch pedigree was reported previously [5]. To
summarize, all MS F0 sires were homozygous (G/G) for
the wild allele, and allelic heterogeneity was identified
for the LW F0 dams: in the two Dutch LR lines, all the
dams were homozygous (G/G), whereas in the three
other LW lines allele A was found with frequencies over
80%. Among the 24 sire families selected for our study,
12 F1 sires were homozygous (G/G), and 12 F1 sires
were heterozygous (A/G). These 24 F1 sire families
involved 65 heterozygous females (A/G) and 71 homo-
zygous females (G/G), while the genotype of 38 F1 dams
remained unknown.

All F0 and F1 animals were genotyped in the French
pedigree (Table 1). All MS F0 dams were homozygous
G/G. Among the six LW sires, five were heterozygous
(A/G) and one was homozygous for the mutant allele
(A/A). Among the six F1 sires, four were homozygous
(G/G) and two were heterozygous (A/G). These six F1
sire families involved 15 heterozygous females (A/G)
and five homozygous females (G/G).

The genotypes of the F2 pigs at the IGF2-intron3-
G3072A were inferred from the genotypes of their par-
ents at the mutation and the haplotypes inherited at the
surrounding SW2443 and SWC9 microsatellite markers.
A complete genotype at the IGF2 mutation could be
obtained for 90% of the F2 pigs. Analyses of variance to
test the inheritance pattern of the IGF2 mutation were
thus performed on 980 F2 animals (543 Dutch F2 and
437 French F2). For ANOVA studies with the IGF2pat-
mat model, the heterozygous (A/G) F2 pigs were split
into two groups depending on the parental origin of the

two alleles. For the combined pedigree, the total num-
bers of animals of each genotype at the mutation were
44 A/A, 568 G/G, 155 A/G and 213 G/A, with the first
allele identifying the paternal allele.

QTL detection

First, each pedigree was analysed independently. In the
French pedigree, the maximum of the test statistic was
obtained in the IGF2 region (0 cM) but was only sig-
nificant at the 10% threshold. Analysis of the Dutch
pedigree gave a significant result at the 5% threshold,
but the maximum of the test statistics was reached at
28 cM.

The QTL detection analysis was then performed on
the combined pedigree (Figure 1). The maximum LRT
value was obtained in the region surrounding the IGF2
position. However, between 13 and 51 cM, the values of
the test statistics were also higher than the 5%
threshold.

A multi-QTL analysis was then performed with the
combined pedigree but neither significant nor suggestive
results were obtained for the hypothesis of two QTL
segregating within both pedigrees.

The QTL detection analysis performed on the 14
families from sires heterozygous at the mutation revealed
a significant QTL close to the IGF2 locus (Figure 2). The
decrease of the test statistic values downstream from the
IGF2 gene was abrupt and no other region reached the
5% threshold. A complementary analysis was performed
on the 16 families originating from homozygous F1 sires
(G/G) and detected a significant QTL at 44 cM. The sub-
stitution effects estimated at this second QTL position
showed that, among the 16 sires analysed, three F1 sires
could not be validated as heterozygous for the QTL. The
13 remaining sires were heterozygous with MS alleles
associated with high BFT values in nine families and with
low BFT values in two families. For the two remaining
sires, the breed origin of the favourable allele could not
be determined. On average, the QTL effect was estimated
to be 0.32 s.d. of the trait. A similar result was obtained
with the combined pedigree using phenotypic data cor-
rected for the effect of the IGF2-intron3-G3072A geno-
type (data not shown). These results clearly indicate that

<table><caption>Table 1 Distribution of genotypes at IGF2-intron3-G3072-A substitution</caption>
  <thead>
    <tr>
      <th rowspan="2"></th>
      <th colspan="5">Dutch pedigree</th>
      <th colspan="5">French pedigree</th>
    </tr>
    <tr>
      <th>A/A</th>
      <th>A/G</th>
      <th>G/G</th>
      <th>unknown</th>
      <th>total</th>
      <th>A/A</th>
      <th>A/G</th>
      <th>G/G</th>
      <th>unknown</th>
      <th>total</th>
    </tr>
  </thead>
  <tbody>
    <tr>
      <th>F0 males</th>
      <td>0</td>
      <td>0</td>
      <td>19</td>
      <td>0</td>
      <td>19</td>
      <td>0</td>
      <td>5</td>
      <td>1</td>
      <td>0</td>
      <td>6</td>
    </tr>
    <tr>
      <th>F0 females</th>
      <td>30</td>
      <td>21</td>
      <td>22</td>
      <td>27</td>
      <td>100</td>
      <td>0</td>
      <td>0</td>
      <td>6</td>
      <td>0</td>
      <td>6</td>
    </tr>
    <tr>
      <th>F1 males</th>
      <td>0</td>
      <td>12</td>
      <td>12</td>
      <td>0</td>
      <td>24</td>
      <td>0</td>
      <td>2</td>
      <td>4</td>
      <td>0</td>
      <td>6</td>
    </tr>
    <tr>
      <th>F1 females</th>
      <td>0</td>
      <td>65</td>
      <td>71</td>
      <td>38</td>
      <td>174</td>
      <td>0</td>
      <td>15</td>
      <td>5</td>
      <td>0</td>
      <td>20</td>
    </tr>
    <tr>
      <th>F2</th>
      <td>342</td>
      <td>180</td>
      <td>20</td>
      <td>23</td>
      <td>565</td>
      <td>226</td>
      <td>188</td>
      <td>22</td>
      <td>85</td>
      <td>521</td>
    </tr>
  </tbody>
</table>

Numbers of heterozygous F2 are given regardless parental origin of alleles.
<figure><img src="image_2.png" /></figure>

**Figure 1** QTL detection on SSC2 in the three studied pedigrees. Solid, circled and crossed lines represent respectively the combined, French and Dutch pedigrees; for each analysis, the LRT is presented as a proportion of the 5% threshold on the chromosome.

a significant QTL affecting BFT is segregating around 40 cM on SSC2.

## Mode of inheritance of the QTL

Analyses of variance were performed on different sub-groups of animals to investigate the effect of the IGF2 mutation. Results obtained using the IGF2pat model confirmed the strong effect of the paternally inherited allele at the mutation (Table 2), since significant p-values were observed in all studied sub-groups of F2 pigs. The p-values obtained with the IGF2patmat model were always very similar to those obtained with the IGF2pat model (data not shown). To investigate the potential effect of the maternally inherited IGF2 allele, the IGF2mat model was also tested. When the analysis was performed on the F2 progeny of heterozygous dams, a significant p-value was obtained with the combined pedigree (p = 0.04). When the analysis was carried out on the progeny of heterozygous dams mated to

homozygous sires, a significant p-value was also observed (p = 0.01). Analysing each pedigree independently, results tended to be significant (p < 0.10) for these two progeny sub-groups in the French pedigree and for the F2 produced from A/G dams and G/G sires in the Dutch pedigree (Table 2).

## Detection of spurious effects of the maternally inherited IGF2-allele

The simulated QTL was detected in about 80% of replicates when its effect was at least 0.32 s.d. regardless of the frequency of allele Q in the European grand-parental population (Figure 3). When the simulated QTL had a small effect (0.22), the French pedigree tended to be more powerful than the Dutch pedigree to detect the QTL. With the Dutch design, the simulated QTL was detected in fewer than 50% of replicates. For the simulations performed with the combined pedigree, the QTL was detected in 88% of replicates.

ANOVA was first carried out with the IGF2pat model, using all families. For both pedigrees, the simulated effect of the paternally inherited allele at IGF2 was detected in most replicates (Figure 4). The Dutch pedigree gave more significant results than the French pedigree. When the frequency of the simulated Q allele increased in the European populations, the percentage of replicates resulting in a significant effect for the paternally inherited allele decreased. With the combined pedigree, 83% of replicates showed a significant effect of IGF2 on backfat thickness.

Using the model of maternal inheritance on simulated imprinted paternally expressed IGF2 effects, the proportion of results reaching significance for an effect of the maternal allele at IGF2 (IGF2mat) was expected to be low or null. Variance analyses were performed on the sub-group of progeny produced by heterozygous dams regardless of the genotypes of the sires. With the Dutch pedigree, few replicates led to validation of the maternal expression. In contrast, with the French pedigree, more significant results were obtained (Figure 4). When the simulated QTL had a large effect (0.42 s.d.) and a low frequency of the Q allele was simulated in the European F0 (0.25), up to 75% of the replicates gave a significant result for the IGF2mat model in the French pedigree.

With the combined pedigree, 6.6% of the simulations detected a significant maternally inherited allele effect. When only progeny from the G/G sires among the heterozygous dams' families were considered, the effect of the allele inherited from the mother at the IGF2 mutation was significant in 23% of replicates.

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

**Figure 2** QTL analyses on SSC2 on sub-groups of the combined pedigree. The solid line represents the QTL detection results from the segregating sire families (A/G sires) and the circled-marked line results from the no-segregating sire families (G/G sires) at the IGF2-intron3-G3072A mutation; for each analysis, the LRT is presented as a proportion of the 5% threshold on the chromosome.

## Discussion

The aim of this study was to confirm the existence of two QTL associated with BFT on SSC2p and to further
<table><caption>Table 2 Statistical analyses of inheritance of the IGF2 mutation effect</caption><thead><tr><th>F2 studied</th><th>Model</th><th>Combined pedigree</th><th>French pedigree</th><th>Dutch pedigree</th></tr></thead><tbody><tr><td>All</td><td>IGF2 pat</td><td>&lt;0.0001 ***</td><td>&lt;0.0001 ***</td><td>&lt;0.0001 ***</td></tr><tr><td>All</td><td>IGF2mat</td><td>0.30</td><td>0.15</td><td>0.95</td></tr><tr><td>From A/G F1 sires</td><td>IGF2pat</td><td>&lt;0.0001 ***</td><td>0.0007 ***</td><td>&lt;0.0001 ***</td></tr><tr><td>From A/G F1 sires</td><td>IGF2mat</td><td>0.94</td><td>0.91</td><td>0.69</td></tr><tr><td>From A/G F1 dams</td><td>IGF2pat</td><td>0.0001 ***</td><td>0.005 **</td><td>0.005 **</td></tr><tr><td>From A/G F1 dams</td><td>IGF2mat</td><td><strong>0.04</strong> *</td><td><strong>0.08</strong> †</td><td>0.26</td></tr><tr><td>From A/G sires and G/G dams</td><td>IGF2pat</td><td>&lt;0.0001 ***</td><td>0.01 *</td><td>0.001 ***</td></tr><tr><td>From A/G sires and G/G dams</td><td>IGF2mat</td><td>-</td><td>-</td><td>-</td></tr><tr><td>From G/G sires and A/G dams</td><td>IGF2pat</td><td>-</td><td>-</td><td>-</td></tr><tr><td>From G/G sires and A/G dams</td><td>IGF2mat</td><td><strong>0.01</strong> *</td><td><strong>0.07</strong> †</td><td><strong>0.08</strong> †</td></tr><tr><td>From A/G sires and A/G dams</td><td>IGF2pat</td><td>0.0004 ***</td><td>0.02 *</td><td>0.008 **</td></tr><tr><td>From A/G sires and A/G dams</td><td>IGF2mat</td><td>0.88</td><td>0.93</td><td>0.97</td></tr></tbody></table>

ANOVA were done on BFT standardized residuals with either the allele inherited at IGF2 from the sire (pat) or from the dam (mat); *** p-value < 0.001, ** p-value <0.01, * p-value < 0.05, † p-value < 0.10; in bold are indicated the p-values < 0.10 obtained with the maternal allelic effect model.

dissect the imprinting effect of the IGF2 gene, in order to resolve contradictory results published for the French and Dutch MS × European pedigrees [5,12]. The two designs were based on similar founder breeds (MS were crossed with LW and/or LR) and they contributed equally (considering the number of F2) to the so-called “combined pedigree”. However, two major differences should be noted. First, the two pedigrees are reciprocal: the MS breed was used as the sire breed in the Dutch pedigree and as the dam breed in the French one. The second discrepancy lies in the pedigree structure: a limited number of large F1 sire families were produced in the French pedigree, whereas the Dutch pedigree consisted in many small F1 sire families.

Our objectives were to (1) show that another QTL on SSC2p is associated with BFT in addition to the IGF2 gene and (2) determine the most likely hypothesis to explain the discrepancies regarding observed imprinting effects in the two pedigrees. Two hypotheses were proposed: a) the imprinting of IGF2 gene is not complete and b) the mendelian effect detected at 0 cM (at the position of the IGF2 locus) in the French pedigree [12] is in fact an artefact due to genetic linkage between IGF2 and a second QTL.

<figure><img src="image_4.png" /><figcaption>**Figure 3 Percentages of replicates for which the simulated QTL was detected, depending on the frequency of the Q allele simulated in the European grand-parental populations.** The QTL was considered detected when the 5% threshold was reached; analyses were performed for each pedigree independently (in grey for the French design, in black for the Dutch design) and for three values for the QTL effect (diamonds 0.22 s.d., squares 0.32 s.d. and triangles 0.42 s.d.).</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>**Figure 4 Percentage of replicates that resulted in a statistically significant effect for the allele inherited from the sire or from the dam at the IGF2 mutation.** The effects were considered significant when the p-value was lower than 5%; the dotted and solid lines represent, respectively, the results obtained for the allele inherited from the sire and from the dam at the IGF2 mutation; different frequencies of the Q allele in the European grand-parental populations were tested; analyses were performed for each pedigree independently (in grey for the French design, in black for the Dutch design) and for three values of the QTL effect (diamonds 0.22 s.d., squares 0.32 s.d. and triangles 0.42 s.d.).</figcaption></figure>
Detection of QTL underlying BFT on the short arm of SSC2

Analysis of the combined pedigree with the single-QTL model gave a most likely position of the QTL at 0 cM, but LRT values over the 5% threshold were obtained for all positions in the first 50 cM of SSC2 (Figure 1). When the two pedigrees were analysed independently, two different most likely QTL positions were obtained, in accordance with previously published results on the two pedigrees, i.e. at 0 cM in the French pedigree [11,17] and in an interval between positions 20 and 50 cM in the Dutch pedigree [22,23]. After genotyping the Dutch pedigree for the IGF2-intron3-G3072A mutation, it was concluded that the IGF2 QTN, localised at 0 cM, explained most of the observed paternally expressed QTL for BFT on SSC2 [5]. However, the presence of an additional QTL around 30 cM could not be excluded. In the present work, nearly every position on the short arm of SSC2 was significant, which is consistent with the large variability of positions found in the Dutch pedigree, and tends to confirm the hypothesis that more than one QTL associated with BFT is segregating in this region. However, this hypothesis had not been validated before for these two pedigrees.

With the combined pedigree, the multiple-QTL model gave no significant result, so that the hypothesis of two co-segregating QTL could not be validated. These results indicate that the F1 sires in general are not heterozygous for both loci and/or that the dataset does not provide enough information to validate the alternative hypothesis, which might be due to the proportion of homozygous sires for the IGF2 mutation in the combined pedigree (0.5).

The separate analyses of the A/G and G/G F1 fathers provided evidence for the segregation of two QTL, the IGF2 mutation and a second one most likely positioned at 44 cM, which is close to the position initially reported in the Dutch pedigree [10]. Segregation of a QTL affecting BFT around 40 cM was also reported by Lee et al. [24] in a Wild Boar × MS pedigree in which all founders were A/A for the IGF2-intron3-G3072A mutation. Therefore, our results confirm the hypothesis of a second QTL segregating in the Dutch LW × MS cross, as suggested by Jungerius et al. [5], and extend this observation to the French cross.

Is the IGF2 gene only paternally expressed?

The IGF2 gene has been studied in detail in several species and, in most cases, paternal expression has been described. Nevertheless, modifications of the imprinted status of genes have been reported in humans (most often associated with diseases) [25]. For the IGF2 gene, such pathological modifications have already been described [26,27]. However, a study on the loss of

imprinting of the IGF2 gene in colorectal cancers has also shown a loss of imprinting in normal mucosa and peripheral blood leukocytes [28]. Moreover, Sakatani et al. [29] have reported maternal expression of IGF2 in a healthy human population. Li et al. [30] have demonstrated that IGF2 P1 transcripts are bi-allelically expressed in all studied organs from adult healthy pigs. In addition, several studies have reported that the imprinting pattern of a locus can be variable over ontogenetic time [31] or under different environmental effects [32].

One present objective was to evaluate whether the differences of the imprinting status previously reported for the French [12] and the Dutch [5] pedigrees could be clarified. Significant effects of the allele inherited from the dam at the IGF2 mutation were obtained when analysing the progeny of A/G mothers, which is not in accordance with exclusive paternal expression of the IGF2 gene. Simulations showed that segregation of a second QTL at 40 cM from IGF2 can lead to the false detection of expression of the maternally inherited IGF2-allele in the French pedigree, whereas this was almost never observed in simulations of the Dutch or the combined pedigree. The discrepancies between simulation results can be explained by differences in the pedigree structures. The Dutch pedigree is based on 104 matings between F0 animals that produced many small half- and full-sib families (24 F1 sires and 174 dams), whereas the French pedigree is based on six matings only between the F0 animals that produced six large half- and full-sib families (six sires and 20 dams). Another difference between the two designs was the allele frequencies at the IGF2-intron3-G3072A mutation: in the French design, 75% of the F1 dams were heterozygous at the mutation whereas only 47% dams where heterozygous in the Dutch design. Since the simulations were performed using the real genotype data at the IGF2-intron3-G3072A mutation, we could not estimate the influence of the percentage of heterogeneous dams. The relative impact of the different family designs and allele frequencies at the mutation could not be differentiated with these simulations. Nevertheless, the ANOVA and simulation analyses suggest that the differences in the effect of the maternal allele at the IGF2 locus can be caused by the segregation of an additional QTL at 44 cM. The genetic linkage between IGF2 and a second QTL with an effect of 0.32 s.d. is high enough to create this artefactual maternal effect, even if the two loci are relatively distant.

Conclusion

Since 2003, several studies have reported the effect of the IGF2-intron3-G3072A mutation on BFT. Besides this QTN, several studies tend to show that additional
loci in the surrounding chromosomal region could influ-
ence the same trait. By combining two F2 designs, this
study demonstrates that a second significant QTL affect-
ing pig fatness is localised around 44 cM and that segre-
gation of this second locus can explain the maternal
effect that was observed in the French pedigree at the
IGF2 locus. Thus, selection schemes against BFT should
not only take the status at the IGF2 mutation into
account but also genotypes at other QTL in the region.

QTL for other economically important traits have
been described on SSC2p, including QTL affecting daily
feed intake [33] and teat number [34]. For these traits,
the influence of the IGF2-intron3-G3072A mutation via
pleiotropic effects has been excluded [5,33]. Therefore,
the short arm of SSC2 seems to be an important chro-
mosomal region for pig production. Thus, fine-mapping
the other QTL on this chromosome will be of major
interest. This task will require a cautious design of fine-
mapping experiments since the pedigree structures and
the variety of loci in the region can lead to false
conclusions.

**Acknowledgements**

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

**Author details**

**Authors' contributions**

FT carried out the QTL detection and statistical analyses. HG participated in those 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 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: 2 December 2010 Accepted: 6 March 2011

Published: 6 March 2011

**References**

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

2. Jeon JT, Carlborg O, Tornsten A, Giuffra E, Amarger V, Chardon P, Andersson-Eklund L, Andersson K, Hansson I, Lundstrom K, Andersson L: A paternally expressed QTL affecting skeletal and cardiac muscle mass in pigs maps to the IGF2 locus. *Nat Genet* 1999, **21**:157-158.

3. Nezer C, Moreau L, Brouwers B, Coppieters W, Detilleux J, Hanset R, Karim L, Kvasz A, Leroy P, Georges M: An imprinted QTL with major effect on muscle mass and fat deposition maps to the IGF2 locus in pigs. *Nat Genet* 1999, **21**:155-156.

4. Van Laere AS, Nguyen M, Braunschweig M, Nezer C, Collette C, Moreau L, Archibald AL, Haley CS, Buys N, Tally M, Andersson G, Georges M, Andersson L: A regulatory mutation in IGF2 causes a major QTL effect on muscle growth in the pig. *Nature* 2003, **425**:832-836.

5. Jungerius BJ, van Laere AS, Te Pas MF, van Oost BA, Andersson L, Groenen MA: The IGF2-intron3-G3072A substitution explains a major imprinted QTL effect on backfat thickness in a Meishan × European white pig intercross. *Genet Res* 2004, **84**:95-101.

6. Carrodeguas JA, Burgos C, Moreno C, Sanchez AC, Ventanas S, Tarrafeta L, Barcelona JA, Lopez MO, Oria R, Lopez-Buesa P: Incidence in diverse pig populations of an IGF2 mutation with potential influence on meat quality and quantity: An assay based on real time PCR (RT-PCR). *Meat Science* 2005, **71**:577-582.

7. Ojeda A, Huang LS, Ren J, Angiolillo A, Cho IC, Soto H, Lemus-Flores C, Makuza SM, Folch JM, Perez-Enciso M: Selection in the making: a worldwide survey of haplotypic diversity around a causative mutation in porcine IGF2. *Genetics* 2008, **178**:1639-1652.

8. Estelle J, Mercade A, Noguera JL, Perez-Enciso M, Ovilo C, Sanchez A, Folch JM: Effect of the porcine IGF2-intron3-G3072A substitution in an outbred Large White population and in an Iberian × Landrace cross. *J Anim Sci* 2005, **83**:2723-2728.

9. Oczkowicz M, Tyra M, Walinowicz K, Rozycki M, Rejduch B: Known mutation (A3072G) in intron 3 of the IGF2 gene is associated with growth and carcass composition in Polish pig breeds. *J Appl Genet* 2009, **50**:257-259.

10. 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.

11. 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.

12. Quintanilla R, Milan D, Bidanel JP: A further look at quantitative trait loci affecting growth and fatness in a cross between Meishan and Large White pig populations. *Genet Sel Evol* 2002, **34**:193-210.

13. Hager R, Cheverud JM, Wolf JB: Maternal effects as the cause of parent-of-origin effects that mimic genomic imprinting. *Genetics* 2008, **178**:1755-1762.

14. Sandor C, Georges M: On the detection of imprinted quantitative trait loci in line crosses: effect of linkage disequilibrium. *Genetics* 2008, **180**:1167-1175.

15. 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.

16. 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.

17. 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.

18. Tortereau F, Gilbert H, Heuven HC, Bidanel JP, Groenen MA, Riquet J: Combining two Meishan F2 crosses improves the detection of QTL on pig chromosomes 2, 4 and 6. *Genet Sel Evol* 2010, **42**:42.

19. 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.

20. 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.

21. 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.

22. de Koning DJ, Rattink AP, Harlizius B, van Arendonk JA, Brascamp EW, Groenen MA: Genome-wide scan for body composition in pigs reveals important role of imprinting. *Proc Natl Acad Sci USA* 2000, **97**:7947-7950.

23. Rattink AP, De Koning DJ, Faivre M, Harlizius B, van Arendonk JA, Groenen MA: Fine mapping and imprinting analysis for fatness trait QTLs in pigs. *Mamm Genome* 2000, **11**:656-661.

24. 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.

25. Falls JG, Pulford DJ, Wylie AA, Jirtle RL: Genomic imprinting: implications for human disease. *Am J Pathol* 1999, **154**:635-647.

<footnote n="1">¹INRA, UMR 0444 Génétique Cellulaire, F-31326 Castanet-Tolosan, France.</footnote>

<footnote n="2">²INRA, UMR 1313 Génétique Animale et Biologie Intégrative, F-78352 Jouy-en-Josas, France. ³Wageningen University, Animal Breeding and Genetics Group, 6700 AH Wageningen, The Netherlands.</footnote>
26. Cui H: Loss of imprinting of IGF2 as an epigenetic marker for the risk of human cancer. *Dis Markers* 2007, **23**:105-112.

27. Jirtle RL: IGF2 loss of imprinting: a potential heritable risk factor for colorectal cancer. *Gastroenterology* 2004, **126**:1190-1193.

28. Cui H, Horon IL, Ohlsson R, Hamilton SR, Feinberg AP: Loss of imprinting in normal tissue of colorectal cancer patients with microsatellite instability. *Nat Med* 1998, **4**:1276-1280.

29. Sakatani T, Wei M, Katoh M, Okita C, Wada D, Mitsuya K, Meguro M, Ikeguchi M, Ito H, Tycko B, Oshimura M: Epigenetic heterogeneity at imprinted loci in normal populations. *Biochem Biophys Res Commun* 2001, **283**:1124-1130.

30. Li C, Bin Y, Curchoe C, Yang L, Feng D, Jiang Q, O'Neill M, Tian XC, Zhang S: Genetic imprinting of H19 and IGF2 in domestic pigs (*Sus scrofa*). *Anim Biotechnol* 2008, **19**:22-27.

31. Wolf JB, Cheverud JM, Roseman C, Hager R: Genome-wide analysis reveals a complex pattern of genomic imprinting in mice. *PLoS Genet* 2008, **4**: e1000091.

32. Hager R, Cheverud JM, Wolf JB: Change in maternal environment induced by cross-fostering alters genetic and epigenetic effects on complex traits in mice. *Proc Biol Sci* 2009, **276**:2949-2954.

33. Houston RD, Haley CS, Archibald AL, Rance KA: A QTL affecting daily feed intake maps to Chromosome 2 in pigs. *Mamm Genome* 2005, **16**:464-470.

34. 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.

doi:10.1186/1297-9686-43-11

Cite this article as: Tortereau *et al*.: Number and mode of inheritance of QTL influencing backfat thickness on SSC2p in Sino-European pig pedigrees. *Genetics Selection Evolution* 2011 **43**:11.

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