Skip to main content

Non-additive effects of RBP4, ESR1 and IGF2 polymorphisms on litter size at different parities in a Chinese-European porcine line



The aim of this work was to study the effects on litter size of variants of the porcine genes RBP4, ESR1 and IGF2, currently used in genetic tests for different purposes. Moreover, we investigated a possible effect of the interaction between RBP4-Msp I and ESR1-Pvu II polymorphisms. The IGF2-intron3-G3072A polymorphism is actually used to select lean growth, but other possible effects of this polymorphism on reproductive traits need to be evaluated.


Detection of polymorphisms in the genomic and cDNA sequences of RBP4 gene was carried out. RBP4-Msp I and IGF2-intron3-G3072A were genotyped in a hyperprolific Chinese-European line (Tai-Zumu) and three new RBP4 polymorphisms were genotyped in different pig breeds. A bivariate animal model was implemented in association analyses considering the number of piglets born alive at early (NBA12) and later parities (NBA3+ ) as different traits. A joint analysis of RBP4-Msp I and ESR1-Pvu II was performed to test their possible interaction. In the IGF2 analysis, paternal or maternal imprinting effects were also considered.


Four different RBP4 haplotypes were detected (TGAC, GGAG, GAAG and GATG) in different pig breeds and wild boars. A significant interaction effect between RBP4-Msp I and ESR1-Pvu II polymorphisms of 0.61 ± 0.29 piglets was detected on NBA3+. The IGF2 analysis revealed a significant increase on NBA3+ of 0.74 ± 0.37 piglets for the paternally inherited allele A.


All the analyzed pig and wild boar populations shared one of the four detected RBP4 haplotypes. This suggests an ancestral origin of the quoted haplotype. The joint use of RBP4-Msp I and ESR1-Pvu II polymorphisms could be implemented to select for higher prolificacy in the Tai-Zumu line. In this population, the paternal allele IGF2-intron3-3072A increased litter size from the third parity. The non-additive effects on litter size reported here should be tested before implementation in other pig breeding schemes.


The use of molecular information in pig breeding programs may enhance genetic gains by increasing the accuracy of genetic evaluation and decreasing generation intervals [1]. More than twelve single nucleotide polymorphisms (SNP) on candidate porcine genes have been associated with litter size or with its main components [2] and some genetic tests have been developed and implemented by breeding companies. For example, variants of the genes ESR1, PRLR, RBP4 and FSHB have been shown to have effects on litter size ranging from 0.25 to over 1 piglet per litter [3].

The retinol binding protein 4 (RBP4) gene codes for a member of the RBP protein family present in the uterus and in embryos during the early stages of gestation [4]. These proteins bind retinol, the bound retinol is then internalized by the cells and triggers embryogenesis [5]. Messer et al. [6] have proposed RBP4 as a possible candidate gene associated with litter size. Subsequently, Rothschild et al. [7], have carried out a study on animals from six commercial lines and reported a significant effect of an intronic polymorphism, the RBP4-Msp I, on the total number of born piglets. Many other studies have shown the existence of a relationship between this polymorphism and litter size [812].

The protein coded by the estrogen receptor 1 (ESR1) gene promotes the expression of different transcription factors involved in the reproductive function of female tissues (ovaries, cervix, uterus...). The ESR1-PvuII polymorphism has been studied previously in the Tai-Zumu line by our group but no significant effect on litter size was observed [13]. Recently, Gonçalves et al. [14] have performed an interesting study in a commercial population that revealed a significant interaction on litter size between RBP4-Msp I and ESR1-PvuII polymorphisms.

A polymorphism detected in the porcine insulin-like growth factor 2 (IGF2) gene, the IGF2-intron3-G3072A SNP [15], has been described as the causal factor of the SSC2 imprinted QTL, which affects fat deposition and muscle growth [16, 17]. Pigs inheriting the paternal allele A have lower backfat thickness and higher lean growth. These effects have been confirmed in different experimental crosses and commercial populations [1820]. Thus, it is likely that allele A has been favored in populations where artificial selection has focused on decreasing fat deposition and increasing lean content. IGF2 is a peptide hormone that participates in the IGF axis, which plays an important role in the promotion of cell proliferation and in the inhibition of apoptosis [21]. Some authors have demonstrated a direct participation of IGF2 in the reproductive function in mouse and farm animals [22, 23]. In addition, selection on lean growth and consequent decrease of fat percentage could reduce prolificacy since larger litter sizes impose greater demand on the sow's energy reserves [24]. Therefore, selecting the paternal inherited allele A could have undesired effects on litter size, which should be evaluated [3].

Estimates of the genetic parameters of litter size in pigs are usually obtained using repeatability models where different parities are considered as different records of the same trait. However, various results support the hypothesis that early and later parities may be partially controlled by different genes and should be considered as different traits. Therefore the use of multitrait models would be more appropriate [2527].

The aim of this research was to study the possible effects of porcine RBP4, ESR1 and IGF2 polymorphisms on the prolificacy of a hyperprolific Chinese-European composite pig line. For this purpose, the detection of new polymorphisms in the RBP4 gene and analysis of their possible effects on litter size were carried out. In addition, the interaction between RBP4 and ESR1 polymorphisms was investigated on our material. The IGF2-intron3-G3072A polymorphism, already used in selection to increase lean growth, was analyzed in order to check if selection on the paternal allele A could affect litter size. All the analyses were carried out using a bivariate model to discriminate the genetic effects on early and later parities.



Research protocols followed the guidelines stated in the Guide for the Care and Use of Agricultural Animals in Agricultural Research and Teaching (FASS, 1999). Data from a Chinese-European composite dam line (Tai-Zumu) were provided by GENE+. This line was developed from Meishan and Jiaxing sows inseminated by hyperprolific French Large White boars, and it was selected for lean growth during nine generations [28]. The pedigree available for this composite line contained 2973 animals of which 2570 sows had 6472 litter size records distributed among 59 farm-year-season classes. The number of litters per parity class is reported in Table 1. Different subsets of genotyped sows were used for the different association analyses carried out.

Table 1 Estimates of heritabilities and genetic correlations for litter size at different parities

Sequencing of the porcine RBP4 gene

Genomic DNA was isolated from blood samples according to a standard protocol [29]. Total RNA was extracted with Tri-Reagent (Sigma-Aldrich Chemie, Madrid, Spain) from liver samples. First strand cDNA was synthesized using 5 μg of total RNA, Superscript™ II Reverse Transcriptase (Invitrogen, Life Technologies, Barcelona) and random hexamers following the supplier's instructions.

The PCR reactions were performed in a 25 μL final volume containing standard buffer (75 mM Tris-HCl pH 9.0, 50 mM KCl, 20 mM (NH4)2SO4), MgCl2 concentrations optimized for each amplified fragment (Additional file 1, Table S1), 200 μM dNTP, 0.5 μM of each primer, 0.5 U of Tth polymerase (Biotools, Madrid, Spain) and 70 ng of genomic DNA or 2 μL of cDNA. Thermocycling conditions were as follows: 94°C (5 min), 40 cycles at 94°C (30 s), the specific annealing temperature (Additional File 1, Table S1) for each primer pair (45 s) and 72°C (45 s), with a final extension step at 72°C (10 min). The amplified products were sequenced using BigDye-Terminator Cycle Sequencing 3.0 in an ABI 3730 automatic sequencer (Applied Biosystems, Warrington, UK). The sequences were edited and aligned using Winstar software.

A 565 bp fragment spanning from exon 2 to 4 of the RBP4 gene was amplified from genomic DNA samples of three Tai-Zumu individuals using the PCR protocol published by Rothschild et al [7]. These authors reported an RBP4-MspI polymorphism but the exact information about its location was not available. The final sequence was submitted to GenBank (accession number: GU932906). Moreover, two overlapping RBP4 cDNA fragments spanning from exon 2 to 6 and covering the complete coding sequence (CDS) were amplified from Tai-Zumu individuals. The primer pairs (RBP4F1-RBP4R1 and RBP4F2-RBP4R2, Additional File 1, Table S1) were designed from the available porcine RBP4 mRNA sequence (GenBank accession number: NM_214057).

SNP genotyping

Five intronic and one exonic SNP were detected in the RBP4 sequences obtained. One of the intronic SNP (c.249-63G>C) was identified as the RBP4-MspI polymorphism previously reported by Rothschild et al [7]. This SNP was genotyped on genomic DNA samples using the published PCR-RFLP protocol. Allele G named as restriction pattern 1 corresponds to three main bands of 190/157/134 bp and allele C named as restriction pattern 2 corresponds to four main bands of 190/134/112/45 bp [9]. A pyrosequencing protocol that allowed simultaneous genotyping of three intronic SNP (c.248+15T>G, c.248+16G>A and c.248+27A>T) was developed using primers RBP4F3-RBP4R3-RBP4Pyr3 (Additional File 1, Table S1). In addition to Tai-Zumu individuals, samples from wild boars as well as Iberian, Landrace, Duroc, Large-White and Meishan breeds were also analyzed. RBP4 haplotypes were determined using the PHASE software.

The ESR1-Pvu II genotyping data were taken from Muñoz et al [13] and the IGF2-intron3-G3072A polymorphism was genotyped by pyrosequencing as described by Van Laere et al. [15] in a PSQ HS 96 system (Pyrosequencing AB, Uppsala, Sweden).

Statistical analysis

A multitrait animal model was used to estimate genetic parameters. Under this approach, the numbers of piglets born alive at each one of the six parity classes (1 to 5 and ≥ 6) were treated as different traits.

where y 1 to y ≥ 6 represent litter size records (NBA) at each parity class, β 1 to β ≥ 6 are the vectors of fixed effects for the six different traits considered, which include the genetic line of the litter's sire (Tai-Zumu or Landrace), parity order and farm-year-season, u 1 to u ≥ 6 and e 1 to e ≥ 6 are vectors of random additive genetics and residual effects for each trait, respectively. Matrices X 1 to X ≥ 6 and Z 1 to Z ≥ 6 are incidence matrixes that associate respectively elements of β 1 to β ≥ 6 and u 1 to u ≥ 6 with the records in y 1 to y ≥ 6 .p ≥ 6 is the vector of permanent environmental effects for each sow with records in the last parity class being W the incidence matrix relating the elements of p ≥ 6 with the records in y ≥6 . The expectation of y i (i = 1 to 5 and ≥ 6) is X i β i and the variance-covariance structure of random effects was assumed to be:

where and are the direct additive genetic and residual variances for trait i, respectively, is the direct genetic covariance between trait i and j (j = 1 to 5 and ≥ 6) and their residual covariance.

A preliminary analysis of the whole data set was performed using this multitrait model. Then, a bivariate model was used to carry out a subsequent analysis of litter size data. In this model, the number of piglets born alive at the first and second parity (NBA12) and the number of piglets born alive at the third and subsequent parities (NBA3+ ) were considered as two different traits. The reduced model can be written as:

Finally, three specific bivariate models were used for the different association analyses, depending on indicator variable values included in X matrices:

  1. i.

    Mendelian inheritance: used in the analysis of the effect of RBP4, ESR1 and IGF2 polymorphisms. It includes additive (α) and dominant (δ) effects. The value of α for each sow depends on her genotype (α = -1, 0, 1) and δ assumes a zero value for homozygote individuals and 1 for the heterozygotes.

  2. ii.

    Mendelian inheritance with epistasis effects: used in the joint analysis of RBP4 and ESR1 polymorphisms. Besides α and δ values, additive x additive interaction (Ψ) effects are also included. Ψ are equal to -1, 0 or 1 depending on the genotypic combination of the analyzed polymorphisms (AA11 = -1; AA12 = 0; AA22 = 1; AB-- = 0; BB11 = -1; BB12 = 0 and BB22 = 1)

  3. iii.

    Paternal or maternal imprinting: used in the analysis of IGF2 SNP. Two association analyses were performed fitting the paternal imprinting effects. Additive and dominant effects were included in the first analysis but not in the second one. In the first analysis, imprinting effects are included (λ) for the heterozygote sows: on the one hand, λ = -1/2 or λ = 1/2 if they have inherited respectively allele G or allele A from the father and on the other hand, λ = 0 for homozygote individuals. In the second analysis, the sows that have received the paternal allele G (GG or GA) have λ = -1/2 and those that have received the paternal allele A (AA or AG) have λ = 1/2. A similar parameterization was used for maternal imprinting.

The statistical significance of each effect was tested comparing the full and reduced models by the χ2 approach to the distribution of the log-likelihood ratios. Variance components and parameter estimates were obtained using VCE-5 program [30] and association analysis were performed using Qxpack package [31]


Variance ratios

Estimated values of heritability (h2 = σ2u2y) for NBA at each parity class and estimated genetic correlations between parities are shown in Table 1. Heritability values for the last two parity classes clearly exceed those of the first four classes. Genetic correlations are greater between adjacent parities, but their values tend to decrease as the number of interspersed parities increases. Although different parities should be considered as different traits, the lower number of genotyped dams compared to the total number of sows requires the use of simpler models to perform the association analyses of this study. According to the structure of the genetic correlations, the records of the first and second parities were grouped in one trait (NBA12) and the remaining in another one (NBA3+ ). Parameter estimates for both traits obtained from the whole data set are shown in Table 2. On the one hand, estimates of parity order effects on NBA12 were expressed as deviation from the first parity (-0.08 ± 0.36) and on NBA3+ as deviation from the third parity (4-3 = 0.25 ± 0.18; 5-3 = -0.04 ± 0.43 & ≥ 6-3 = 0.07 ± 0.69). On the other hand, the estimated effect of genetic line of the litter's sire was not statistically significant, i.e. 0.25 ± 0.18 for NBA12 and -0.07 ± 0.25 for NBA3+.

Table 2 Genealogical data and estimates of heritabilities, permanent environmental effects and correlations between NBA12 and NBA3+

RBP4 and ESR1

After sequencing and aligning the 565 bp genomic fragment of the RBP4 gene, five intronic SNP were detected: c.111+47T>C, located in intron 2 and c.248+15G>T, c.248+16A>G, c.248+27A>T and c.249-63G>C located in intron 3. The c.249-63G>C SNP was identified as the polymorphism RBP4-MspI [7] and corresponds to the second position of a recognition site of the Msp I restriction enzyme (CC GG). Moreover, a silent SNP, c.156G>A was detected on exon 3. In addition, two overlapping cDNA fragments of 485 and 479 bp, respectively, were amplified and sequenced. The assembled fragments form an 861 bp sequence that covers the complete CDS. As a result, the SNP c.156G>A was confirmed but no other exonic polymorphism could be detected. From the comparison of the different sequences, SNP c.111+47T>C, c.156G>A, c.248+15G>T and c.248+16A>G seem to be cosegregating in the sequenced Tai-Zumu individuals. In order to check their segregation pattern, SNP c.248+15G>T, c.248+16A>G, c.248+27A>T and c.249-63G>C (RBP4-Msp I) were genotyped on different domestic pig populations (Tai-Zumu, Duroc, Landrace, Large-White, Meishan and Iberian) and wild boars. The results distinguished four different haplotypes for the quoted positions: TGAC, GGAG, GAAG and GATG. Their respective frequencies in the different populations are shown in Table 3.

Table 3 Haplotypic frequencies of RBP4 gene in different porcine populations

In a first step, an association analysis of the RBP4-MspI SNP was performed in 534 sows with 957 litter size records for NBA12 and 1043 for NBA3+. Allele 1 (frequency = 0.51) was significantly associated with a higher number of piglets born alive in the two first parities (NBA12), but not in the third and subsequent parities (NBA3+ ). The estimated additive effect on NBA12 was 0.42 piglets per litter (P≤0.016), and no dominance effects were observed (Table 4). A separate analysis of ESR1-PvuII SNP was carried out on 403 sows (56 AA, 180 AB and 167 BB), with 733 litter size records for NBA12 and 934 for NBA3+. No significant effect on litter size was evidenced. In addition, a joint analysis between RBP4-MspI and ESR1-PvuII polymorphisms was performed using data from 375 sows with 679 litter size records for NBA12 and 874 for NBA3+. The number of sows for each one of the nine genotypic combinations ranged from 12 (ESR1-PvuII AA/RBP4-MspI 22) to 81 (ESR1-PvuII AB/RBP4-MspI 12). The additive effect of RBP4-MspI on NBA12 was confirmed and a significant interaction effect was detected on NBA3+ (Table 4). The genotypes of the largest litter sizes corresponded to the combinations (ESR1 AA/RBP4 11) and (ESR1 BB /RBP4 22) and the least prolific to the alternative combination (ESR1 BB/RBP4 11) and (ESR1 AA/RBP4 22) (Figure 1). The estimated differences for NBA3+ between both groups of sows are 1.09 ± 0.54 piglets (P < 0.046).

Table 4 Individual and joint analysis of RBP4-Msp I and ESR1-Pvu I effects on NBA12 and NBA3+
Figure 1
figure 1

Interaction effects between genotypes RBP4 -MspI and ESR1 -PvuII on NBA 12 and NBA 3 + .


Results obtained in the different association analyses fitting IGF2 SNP effects are shown in Table 5.A Mendelian inheritance analysis was performed on 550 genotyped sows (192 GG, 264 GA and 94 AA), with 985 records for NBA12 and 1057 records for NBA3+ , but no significant result was obtained. Otherwise, to implement a model of imprinting inheritance requires that the paternal or maternal inheritance of the alleles can be determined in the heterozygote sows. This was possible for 56 of the 264 total heterozygotes: 31 with the paternal allele G and 25 with the paternal allele A. The analysis was performed on 342 sows with 613 records for NBA12 and 710 records for NBA3+. When additive and dominant effects were taken into account, a suggestive additive effect of the paternal allele A was detected on NBA3+ (0.36 ± 0.21, P < 0.052). If only paternal imprinting effects are considered, a significant increase produced by the paternal allele A of the number of piglets alive was detected on NBA3+. Maternal imprinting effects were not evidenced in a complementary analysis (Table 5).

Table 5 Results of association analysis of IGF2-intron3-G3072A SNP with litter size at different parities


If most of the genes affecting NBA at different parities were the same, homogenous heritability estimates and high values of genetic correlations would be expected. However, as shown in Table 1, heterogeneous values of heritability and genetic correlation were found. These results, as others previously obtained from different pig breeds, suggest that different genes or combinations of genes may affect litter size in each one of the parities [2527]. Thus, multitrait models instead of the repeatability model should be used to analyse porcine litter size data, although simpler bivariate models distinguishing early and later parities may be adequate for reduced data sets.

Porcine RBP4 studies performed so far have mainly focused on association analyses between the RBP4-MspI polymorphism and litter size. The current study allowed us to detect four RBP4 haplotypes in six different pig breeds and European wild boars. TGAC is the only haplotype shared by all the populations analyzed and hence it is probably the ancestral haplotype. GGAG was exclusively detected in Meishan and GATG in Meishan and Tai-Zumu. The other haplotype (GAAG) was detected in all the pig breeds and wild boars analyzed except Iberian pigs that only displayed the TGAC haplotype. Some authors have reported introgression of Asian alleles in many European breeds, but not in Iberian pigs [3234]. Our results confirm this and suggest an Asian origin for haplotypes GGAG, GAAG and GATG. The low frequency (0.009) of haplotype GAAG in wild boars can be explained by the existence of uncontrolled mating between wild boars and domestic pigs in a region where wild boars coexist with open air pig production. Another aspect to consider is that the number of detected haplotypes is higher for Meishan individuals than for those from European breeds. This is consistent with Amaral et al. [35] who reported a higher haplotypic diversity and lower proportion of fixed markers in Chinese breeds. Similar situations have already been reported for other genes (PRLR, BMPR1B, ESR1) related to reproductive traits [13, 36].

The GATG haplotype showed a low frequency in the Tai-Zumu population (Table 3) and thus performing an association analysis with one of the SNP instead of the haplotypes seemed more suitable. The SNP chosen was RBP4-Msp I because it presents intermediate allelic frequencies in the population. Given the distribution of haplotypes observed in the Tai-Zumu population, the analysis carried out with the RBP4-MspI SNP would be equivalent to comparing haplotype TGAC to haplotypes GAAG and GATG. Individual and joint association analyses of RBP4-Msp I with NBA12 and NBA3+ revealed a favourable additive effect of allele 1 on NBA12. This result is in accordance with that detected by Rothschild et al. [7]. They have reported a 0.23 piglet/litter effect of the RBP4-MspI allele on the total number of piglets born in six lines from different genetic origins. Also, Spöter et al. [37] have detected both additive and dominant effects of 0.24 and 0.31 piglet/litter on NBA, in the German Landrace breed but not in the German Large-White breed. Similar negative results were obtained by other authors in a Duroc x Large White synthetic line and in a Polish breed [8, 11]. Experiments where frequencies of RBP4-Msp I alleles were compared in control and selected lines for increased litter size did not reveal any significant result [9, 10].

These diverse results indicate that the causal mutation could be in linkage disequilibrium with the porcine RBP4-MspI SNP. Besides, a possible dependence on the genetic background should be taken into account, because epistatic effects could be affecting pig prolificacy as recently reported [38, 39]. Gonçalves et al. [14] have pointed out that effects of the RBP4-MspI polymorphism on litter size depend on the genotype of the ESR1-Pvu II allele in a comparison between sows from three genotypic classes. The litter size for second and later parities of sows carrying either ESR1 allele A/RBP4 genotype 11 or ESR1 allele B/RBP4 genotype 22 was greater than that of sows grouped in the third class (ESR1 AA/RBP4 22 and ESR1 BB/RBP4 11). The results of our joint association analysis allow us to corroborate more precisely the results obtained by Gonçalves et al. [14] i.e., sows with genotypic combinations ESR1 AA/RBP4 11 and ESR1 BB/RBP4 22 were the most prolific for NBA3+. These findings may reflect a physiological interaction between estrogens and RBP4 proteins. Once, the first secretion of RBP has occurred in the embryo, embryonic estrogens are secreted in the maternal uterus where they induce an increase of expression and secretion of RBP proteins. These proteins enter the embryo cells rising the RBP receptors density and allowing the embryo development to continue [40]. Therefore the joint selection of RBP4-MspI and ESR1-PvuII could be implemented to improve prolificacy in Tai-Zumu pigs, although its use in other commercial populations requires confirmation of the observed interaction.

Implementation of molecular markers in selection requires exhaustive verification in order to ensure that no undesirable effect arises in other economically important traits. So far, some studies have been developed to check the effect of IGF2-intron3-G3072A on prolificacy, with uneven results in different populations, although the methodology used and the available information varied among the studies. Using a Mendelian inheritance model, Horak et al. and Katska-Kiazkiewicz et al. [11, 41] have detected significant effects of different IGF2 polymorphisms on litter size in Czech and Polish pigs, respectively. In addition, Rempel et al. [42] have not detected any significant effect of IGF2-intron3-G3072A in a composite pig line. Assuming an imprinting inheritance model, Buys et al. [43] have detected an increase on litter size due to the paternal inherited allele G in dam lines based on Large-White and Landrace breeds. However, in other studies an increase in prolificacy was detected on the heterozygote individuals who inherited the paternal allele A [44, 45].

In the current study, both types of inheritance were taken into account. A significant effect was only detected under the inheritance model of paternal imprinting, i.e. an increase of 0.74 piglet on NBA3+. Hence, it is clear that the results depend on the model employed. Note that imprinting phenomena could arise from CpG island methylation events that trigger the silencing of the genes on a chromosomal region [46, 47]. Indeed, the IGF2-intron3-G3072A mutation is located in a CpG island and its causality on pig lean growth has been well confirmed [48]. Although more studies are required to explain the effects on prolificacy, selection of the paternal IGF2-intron3-G3072A mutation could be implemented in the Tai-Zumu population due to its beneficial effects both on lean growth and litter size in third and subsequent parities.


A multitrait model is recommended to analyze the effects of various polymorphisms on litter size since early and later parities can be partially controlled by different genes.

Analysis of the RBP4 gene in wild boars and six porcine populations allowed to detect four haplotypes. Only one of the four detected haplotypes was shared by all the analyzed pig and wild boar populations indicating an ancestral origin of the quoted haplotype. Otherwise, RBP4-MspI does not seem to be the causative mutation associated with an increase in litter size. However, an interaction effect between RBP4-Msp I and ESR1-Pvu II on NBA3+ was detected in the Tai-Zumu population. According to this, the joint use of the most favorable genotypic combination could be implemented in order to select for higher litter size.

Selecting the paternally inherited IGF2-intron3-3072A allele in Tai Zumu increases litter size from the third parity. The causative mutation could be situated either in the IGF2 gene or very close to this gene.


  1. Dekkers JC: Commercial application of marker- and gene-assisted selection in livestock: strategies and lessons. J Anim Sci. 2004, 82 (E-Suppl): E313-328.

    PubMed  Google Scholar 

  2. Distl O: Mechanisms of regulation of litter size in pigs on the genome level. Reprod Domest Anim. 2007, 42 (Suppl 2): 10-16. 10.1111/j.1439-0531.2007.00887.x.

    Article  PubMed  Google Scholar 

  3. Rothschild MF: Porcine genomics delivers new tools and results: this little piggy did more than just go to market. Genet Res. 2004, 83: 1-6. 10.1017/S0016672303006621.

    CAS  Article  PubMed  Google Scholar 

  4. Trout WE, Hall JA, Stallings-Mann ML, Galvin JM, Anthony RV, Roberts RM: Steroid regulation of the synthesis and secretion of retinol-binding protein by the uterus of the pig. Endocrinology. 1992, 130: 2557-2564. 10.1210/en.130.5.2557.

    CAS  PubMed  Google Scholar 

  5. Yelich JV, Pomp D, Geisert RD: Detection of transcripts for retinoic acid receptors, retinol-binding protein, and transforming growth factors during rapid trophoblastic elongation in the porcine conceptus. Biol Reprod. 1997, 57: 286-294. 10.1095/biolreprod57.2.286.

    CAS  Article  PubMed  Google Scholar 

  6. Messer LA, Wang L, Yelich J, Pomp D, Geisert RD, Rothschild MF: Linkage mapping of the retinol-binding protein 4 (RBP4) gene to porcine chromosome 14. Mamm Genome. 1996, 7: 396-10.1007/s003359900117.

    CAS  Article  PubMed  Google Scholar 

  7. Rothschild MF, Messer L, Day A, Wales R, Short T, Southwood O, Plastow G: Investigation of the retinol-binding protein 4 (RBP4) gene as a candidate gene for increased litter size in pigs. Mamm Genome. 2000, 11: 75-77. 10.1007/s003350010015.

    CAS  Article  PubMed  Google Scholar 

  8. Drogemuller C, Hamann H, Distl O: Candidate gene markers for litter size in different German pig lines. J Anim Sci. 2001, 79: 2565-2570.

    CAS  PubMed  Google Scholar 

  9. Linville RC, Pomp D, Johnson RK, Rothschild MF: Candidate gene analysis for loci affecting litter size and ovulation rate in swine. J Anim Sci. 2001, 79: 60-67.

    CAS  PubMed  Google Scholar 

  10. Blowe CD, Boyette KE, Ashwell MS, Eisen EJ, Robison OW, Cassady JP: Characterization of a line of pigs previously selected for increased litter size for RBP4 and follistatin. J Anim Breed Genet. 2006, 123: 389-395. 10.1111/j.1439-0388.2006.00620.x.

    CAS  Article  PubMed  Google Scholar 

  11. Katska-Ksiazkiewicz L, Lechniak-Cieslak D, Korwin-Kossakowska A, Alm H, Rynska B, Warzych E, Sosnowski J, Sender G: Genetical and biotechnological methods of utilization of female reproductive potential in mammals. Reprod Biol. 2006, 6 (Suppl 1): 21-36.

    PubMed  Google Scholar 

  12. Spotter A, Distl O: Genetic approaches to the improvement of fertility traits in the pig. Vet J. 2006, 172: 234-247. 10.1016/j.tvjl.2005.11.013.

    CAS  Article  PubMed  Google Scholar 

  13. Muñoz G, Ovilo C, Estelle J, Silio L, Fernandez A, Rodriguez C: Association with litter size of new polymorphisms on ESR1 and ESR2 genes in a Chinese-European pig line. Genet Sel Evol. 2007, 39: 195-206. 10.1186/1297-9686-39-2-195.

    PubMed Central  Article  PubMed  Google Scholar 

  14. Gonçalves IDV, Gonçalves PBD, J.C. S, Portela VV, Borges LFK, Oliveira JFC, Lovatto PA: Interaction between estrogen receptor and retinol-binding protein polymorphisms as a tool for the selection of prolific pigs. Genet Mol Biol. 2008, 31: 481-486. 10.1590/S1415-47572008000300014.

    Article  Google Scholar 

  15. 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. 10.1038/nature02064.

    CAS  Article  PubMed  Google Scholar 

  16. 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. 10.1038/5938.

    CAS  Article  PubMed  Google Scholar 

  17. 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. 10.1038/5935.

    CAS  Article  PubMed  Google Scholar 

  18. 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 x European white pig intercross. Genet Res. 2004, 84: 95-101. 10.1017/S0016672304007098.

    CAS  Article  PubMed  Google Scholar 

  19. 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 x Landrace cross. J Anim Sci. 2005, 83: 2723-2728.

    CAS  PubMed  Google Scholar 

  20. Heuven HCM, Bovenhuis H: Effect of IGF2 on growth characteristics of F2 Meishan x White. Proceedings of 56th Annual Meeting of the European Association for Animal Production: 5-8 June 2005; Upssala. 2005

    Google Scholar 

  21. Oksbjerg N, Gondret F, Vestergaard M: Basic principles of muscle development and growth in meat-producing mammals as affected by the insulin-like growth factor (IGF) system. Domest Anim Endocrinol. 2004, 27: 219-240. 10.1016/j.domaniend.2004.06.007.

    CAS  Article  PubMed  Google Scholar 

  22. Badinga L, Song S, Simmen RC, Clarke JB, Clemmons DR, Simmen FA: Complex mediation of uterine endometrial epithelial cell growth by insulin-like growth factor-II (IGF-II) and IGF-binding protein-2. J Mol Endocrinol. 1999, 23: 277-285. 10.1677/jme.0.0230277.

    CAS  Article  PubMed  Google Scholar 

  23. Schams D, Berisha B, Kosmann M, Einspanier R, Amselgruber WM: Possible role of growth hormone, IGFs, and IGF-binding proteins in the regulation of ovarian function in large farm animals. Domest Anim Endocrinol. 1999, 17: 279-285. 10.1016/S0739-7240(99)00044-2.

    CAS  Article  PubMed  Google Scholar 

  24. Mathur P, Liu Y: Marker assisted selection for the Canadian Swine industry. Proceedings of the 28th Annual National Swine Improvement Federation; 4-5 December 2003; Des Moines. 2003, 146-149.

    Google Scholar 

  25. Hanenberg EHAT, Knol EF, Mercks JWM: Estimates of genetic parameters for reproduction traits at different parities in Dutchs Landrace pigs. Livest Prod Sci. 2001, 69: 179-186. 10.1016/S0301-6226(00)00258-X.

    Article  Google Scholar 

  26. Serenius T, Sevón-Aimonen ML, Mäntyasaari EA: Effect of service sire and validity of repeatibility model in litter size and farrowing interval of Finnish Landrace and Large White populations. Livest Prod Sci. 2003, 81: 213-222. 10.1016/S0301-6226(02)00300-7.

    Article  Google Scholar 

  27. Fernández A, Rodrigañez J, Zuzúarregui J, Rodriguez MC, Silio L: Genetic parameters for litter size and weight at different parities in Iberian pigs. Span Journ Agr Res. 2008, 6: 98-106.

    Article  Google Scholar 

  28. Milan D, Demeure O, Laval G, Iannuccelli N, Genet C, Bonnet M, Burgaud G, Riquet J, Gasnier C, Bidanel JP: Identification de régions dugénome répondant à la sélection dans une lignée porcine sino-européenne: la Tai- zumu. Journ Recher Porc. 2003, 35: 309-316.

    Google Scholar 

  29. Sambrook J, Fritsh EF, Maniatis T: Molecular Cloning: A Laboratory Manual. 1998, Cold Spring Harbor: Cold Spring Harbor Laboratory Press

    Google Scholar 

  30. Kovac M, Groenveld E, L.A GC: VCE-5, a package for the estimation of dispersion parameters. Proceedings of the 7th World Congress on Genetics Applied to Livestock Production: 19-23 August 2002; Montpellier. 2002, 741-742.

    Google Scholar 

  31. Perez-Enciso M, Misztal I: Qxpak: a versatile mixed model application for genetical genomics and QTL analyses. Bioinformatics. 2004, 20: 2792-2798. 10.1093/bioinformatics/bth331.

    CAS  Article  PubMed  Google Scholar 

  32. Alves E, Ovilo C, Rodriguez MC, Silio L: Mitochondrial DNA sequence variation and phylogenetic relationships among Iberian pigs and other domestic and wild pig populations. Anim Genet. 2003, 34: 319-324. 10.1046/j.1365-2052.2003.01010.x.

    CAS  Article  PubMed  Google Scholar 

  33. Clop A, Amills M, Noguera JL, Fernandez A, Capote J, Ramon MM, Kelly L, Kijas JM, Andersson L, Sanchez A: Estimating the frequency of Asian cytochrome B haplotypes in standard European and local Spanish pig breeds. Genet Sel Evol. 2004, 36: 97-104. 10.1186/1297-9686-36-1-97.

    PubMed Central  Article  PubMed  Google Scholar 

  34. 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. 10.1534/genetics.107.084269.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  35. Amaral AJ, Megens HJ, Crooijmans RP, Heuven HC, Groenen MA: Linkage disequilibrium decay and haplotype block structure in the pig. Genetics. 2008, 179: 569-579. 10.1534/genetics.107.084277.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  36. Tomás A: Analysis of QTL and candidate genes related to prolificaciy and piglet survivability in an Iberian x Meishan F2 population. Universitat Autònoma de Barcelona, Departament de Ciència Animal i dels Aliments. 2006

    Google Scholar 

  37. Spotter A, Muller S, Hamann H, Distl O: Effect of polymorphisms in the genes for LIF and RBP4 on litter size in two German pig lines. Reprod Domest Anim. 2009, 44: 100-105. 10.1111/j.1439-0531.2007.01004.x.

    CAS  Article  PubMed  Google Scholar 

  38. Noguera JL, Rodriguez C, Varona L, Tomas A, Munoz G, Ramirez O, Barragan C, Arque M, Bidanel JP, Amills M, Ovilo C, Sanchez A: A bi-dimensional genome scan for prolificacy traits in pigs shows the existence of multiple epistatic QTL. BMC Genomics. 2009, 10: 636-10.1186/1471-2164-10-636.

    PubMed Central  Article  PubMed  Google Scholar 

  39. Fernandez-Rodriguez A, Rodriguez C, Varona L, Balcells I, Noguera JL, Ovilo C, Fernandez AI: Analysis of candidate genes underlying two epistatic quantitative trait loci on SSC12 affecting litter size in pig. Anim Genet. 2010, 41: 73-80. 10.1111/j.1365-2052.2009.01966.x.

    CAS  Article  PubMed  Google Scholar 

  40. Harney JP, Ott TL, Geisert RD, Bazer FW: Retinol-binding protein gene expression in cyclic and pregnant endometrium of pigs, sheep, and cattle. Biol Reprod. 1993, 49: 1066-1073. 10.1095/biolreprod49.5.1066.

    CAS  Article  PubMed  Google Scholar 

  41. Horák P, Miková G, Urban T, Putnová L, Knoll A, Dvorak J: Association of polymorphism in the IGF2 gene with litter size in Black Pied Prestice pigs. Czech J Anim Sci. 2001, 46 (11): 505-508.

    Google Scholar 

  42. Rempel LA, Nonneman DJ, Wise TH, Erkens T, Peelman LJ, Rohrer GA: Association analyses of candidate SNP on reproductive traits in swine. J Anim Sci. 2010, 88 (1): 1-15. 10.2527/jas.2009-1985.

    CAS  Article  PubMed  Google Scholar 

  43. Buys N, Van den Abeele A, Stinckens A, Deley J, Georges M: Effect of the IGF2-intron3-G3072A mutation on prolificacy in sows. Proceedings of 8th Congress on Genetic Applied to Livestock Production: 13-18 August 2006; Belo Horizonte. 2006, CD

    Google Scholar 

  44. Heuven HCM, Coster A, Madsen O, Bovenhuis H: Bi-polar imprinting of IGF2 affects litter size in Meishan-F2 crossbred sows. Proceedings of 60th Annual Meeting of the European Association for Animal Production; 24-27 August 2009; Barcelona. 2009, Wageningen Academic Publishers

    Google Scholar 

  45. Jafarikia M, Maignel L, Wyss S, Berkel WV, Sullivan B: Effect of IGF2 gene on sow productivity traits. Proceedings of 60th Annual Meeting of the European Association for Animal Production; -27 August 2009; Barcelona. 2009, Wageningen Academic Publishers

    Google Scholar 

  46. 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. 10.1073/pnas.140216397.

    PubMed Central  CAS  Article  PubMed  Google Scholar 

  47. Hore TA, Rapkins RW, Graves JA: Construction and evolution of imprinted loci in mammals. Trends Genet. 2007, 23: 440-448. 10.1016/j.tig.2007.07.003.

    CAS  Article  PubMed  Google Scholar 

  48. Markljung E, Jiang L, Jaffe JD, Mikkelsen TS, Wallerman O, Larhammar M, Zhang X, Wang L, Saenz-Vash V, Gnirke A, Lindroth AM, Barrés R, Yan J, Strömberg S, De S, Pontén F, Lander ES, Carr SA, Zierath JR, Kullander K, Wadelius C, Lindblad-Toh K, Andersson G, Hjälm G, Andersson L: ZEBD6 a novel transcription factor derived from a domesticated DNA transposon regulates IGF2 expression and muscle growth. Plos Biology. 7 (12): e1000256-10.1371/journal.pbio.1000256.

Download references


Financial support was provided by Spanish MCYT grant FIT01-0000-2001-027.We are grateful to Gene +, especially to Fernando Flamarique, Michel Sourdioux and Christian Gasnier for supplying data and blood samples. We gratefully acknowledge to María Angeles López and Rita Benitez for technical support and to Beatriz Villanueva for her valuables suggestions. M. Muñoz is funded by a PhD INIA grant.

Author information

Authors and Affiliations


Corresponding author

Correspondence to María Muñoz.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

MM carried out the polymorphism detection and the genotyping tasks in the RBP4 gene, drafted and finalized the manuscript. AIF carried out the genotyping of the IGF2-intron3-G3072A polymorphism. CO and GM carried out the genotyping task of the ESR1-Pvu II polymorphism. AF performed the statistical analysis and helped to revise the manuscript. EA participated in the design of the study of RBP4 gene, helped to draft, revise and complete the manuscript. LS and CR conceived, coordinated and led the project. Besides LS participated in revising and finishing the manuscript.

All authors read and approved the final manuscript.

Electronic supplementary material


Additional file 1: Table S1 - Primer sequences, annealing temperatures, MgCl 2 concentrations and amplicon sizes used for RBP4 sequencing and pyrosequencing. This table shows primers used for RBP4 sequencing and pyrosequencing. Annealing temperature, MgCl2 concentration and amplification size are indicated for each fragment. (PDF 72 KB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Rights and permissions

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

Reprints and Permissions

About this article

Cite this article

Muñoz, M., Fernández, A.I., Óvilo, C. et al. Non-additive effects of RBP4, ESR1 and IGF2 polymorphisms on litter size at different parities in a Chinese-European porcine line. Genet Sel Evol 42, 23 (2010).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


  • Wild Boar
  • Single Nucleotide Polymorphism
  • Litter Size
  • Paternal Allele
  • Wild Boar Population