Email updates

Keep up to date with the latest news and content from Genetics Selection Evolution and BioMed Central.

Open Access Research

Genome-wide mapping of Quantitative Trait Loci for fatness, fat cell characteristics and fat metabolism in three porcine F2 crosses

Hermann Geldermann1*, Stanislav Čepica2, Antonin Stratil2, Heinz Bartenschlager3 and Siegfried Preuss3

Author Affiliations

1 Animal Breeding and Biotechnology, University of Hohenheim, Stuttgart, Germany

2 Institute of Animal Physiology and Genetics, Academy of Sciences of the Czech Republic, Liběchov, Czech Republic

3 Department of Animal Breeding and Biotechnology, University of Hohenheim, Stuttgart, Germany

For all author emails, please log on.

Genetics Selection Evolution 2010, 42:31  doi:10.1186/1297-9686-42-31

The electronic version of this article is the complete one and can be found online at:

Received:26 January 2010
Accepted:28 July 2010
Published:28 July 2010

© 2010 Geldermann et al; licensee 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.



QTL affecting fat deposition related performance traits have been considered in several studies and mapped on numerous porcine chromosomes. However, activity of specific enzymes, protein content and cell structure in fat tissue probably depend on a smaller number of genes than traits related to fat content in carcass. Thus, in this work traits related to metabolic and cytological features of back fat tissue and fat related performance traits were investigated in a genome-wide QTL analysis. QTL similarities and differences were examined between three F2 crosses, and between male and female animals.


A total of 966 F2 animals originating from crosses between Meishan (M), Pietrain (P) and European wild boar (W) were analysed for traits related to fat performance (11), enzymatic activity (9) and number and volume of fat cells (20). Per cross, 216 (M × P), 169 (W × P) and 195 (W × M) genome-wide distributed marker loci were genotyped. QTL mapping was performed separately for each cross in steps of 1 cM and steps were reduced when the distance between loci was shorter. The additive and dominant components of QTL positions were detected stepwise by using a multiple position model.


A total of 147 genome-wide significant QTL (76 at P < 0.05 and 71 at P < 0.01) were detected for the three crosses. Most of the QTL were identified on SSC1 (between 76-78 and 87-90 cM), SSC7 (predominantly in the MHC region) and SSCX (in the vicinity of the gene CAPN6). Additional genome-wide significant QTL were found on SSC8, 12, 13, 14, 16, and 18. In many cases, the QTL are mainly additive and differ between F2 crosses. Many of the QTL profiles possess multiple peaks especially in regions with a high marker density. Sex specific analyses, performed for example on SSC6, SSC7 and SSCX, show that for some traits the positions differ between male and female animals. For the selected traits, the additive and dominant components that were analysed for QTL positions on different chromosomes, explain in combination up to 23% of the total trait variance.


Our results reveal specific and partly new QTL positions across genetically diverse pig crosses. For some of the traits associated with specific enzymes, protein content and cell structure in fat tissue, it is the first time that they are included in a QTL analysis. They provide large-scale information to analyse causative genes and useful data for the pig industry.


Reduced fatness improves carcass value, and therefore numerous studies on QTL mapping in pig concern fat deposition related traits (see reviews [1,2]). More recently, the results have been compiled in the database PigQTLdb ([3,4]; webcite). As shown in several studies, QTL profiles depend largely on genetic resources, trait definition and statistical models. Taken together, these studies have detected major QTL affecting fat traits on porcine chromosomes SSC1, 2, 4, 6, 7 and X.

Traits like volume of adipose tissue and fat metabolism are influenced by lipogenesis and lipolysis rates, relationship between lipogenesis and lipolysis, energy intake and adipocyte differentiation. In pig, fat accretion is related to the activity of NADPH-generating enzymes in adipose tissue [5]. Strutz [6] has reported genetic correlations of about -0.4 to -0.6 between carcass fat content and activity of NADPH-generating enzymes. The content of soluble proteins in porcine fat tissue is an indicator of metabolic activity and has been reported to be genetically correlated (about -0.5) with fat content in carcass [7]. Furthermore, data on the diameter and number of porcine fat cells and on cell size differences between lean and obese pigs have been reported [8,9].

Activity of specific enzymes, protein content and cell structure in fat tissue probably depend on a smaller number of genes than production traits related to fat content in carcass. Thus, we have measured metabolic and cytological features for back fat tissue together with performance traits related to carcass fat deposition and used these traits in a genome-wide QTL analysis.

The positions of the QTL were compared among three F2 porcine crosses as well as between male and female animals. For some traits, we analysed the combined influence of QTL positioned on different chromosomes on the trait variance. We detected a total of 76 QTL (P < 0.05) and 71 QTL (P < 0.01) with genome-wide significant effects for the three crosses, but numerous QTL were observed only in one or two of the crosses.



A total of 966 F2 pigs were generated with founder animals from the Meishan and Pietrain breeds and the European wild boar (Table 1). All pigs were maintained under standardized housing in one experimental station. Generation of animals for the three F2 crosses and conditions of feeding are described elsewhere [2,10].

Table 1. Pedigrees of the three F2 crosses with animal numbers used in the calculations


Blood samples were collected from founders, F1 and F2 animals. Blood was taken from the v. jugularis of living animals or during stunning and separated into plasma, erythrocytes and leucocytes. DNA was isolated from the leukocyte fraction by chloroform-phenol extraction according to standard protocols.

Adipose tissue of the back fat area between the skin and m. longissimus dorsi at 13th/14th rib was collected from the F2 animals directly after stunning. For each animal, a piece of back fat tissue was sampled and stored immediately in liquid nitrogen. After thawing the subcutaneous adipose tissue at the connective tissue border was separated into an inner and outer layer sample. For both samples, connective tissue and blood vessels were removed and the adipose tissue used immediately.

Trait measurements

As shown in Table 2, 40 traits were recorded, including 11 performance traits associated with fat deposition (Table 2a). Six other traits related to enzyme activities and three to protein content were measured in fat tissue (Table 2b). The relative numbers or volumes of fat cells were determined using different parameters defining 20 traits (Table 2c). Traits related to protein content, enzyme activities and fat cells are described in the following sections.

Table 2. Definition of traitsa

Soluble proteins and enzymes

Each fat tissue sample was cut into small pieces (about 1 mm thick) and then homogenized at 0°C in a 0.15 M KCl solution. The homogenate was centrifuged (20 min, 20000 g, +4°C) and the supernatant filtered (Filter No. 11303, pore diameter 1.2 μm, Sartorius, Göttingen, Germany). The filtrate was kept at +4°C and immediately used to measure protein content and enzyme activities. Protein contents were estimated according to [11]. For each fat tissue sample, protein content was measured three times and averaged. To measure each enzyme activity, 0.1 mL of the filtrate was mixed:

- for isocitrate dehydrogenase (ICDH): with 1.0 mL of 0.075 M glycyl-glycine buffer (pH 7.4), 0.1 mL of 0.05 M MnCl2-4H2O, 0.2 mL of 0.002 M NADP, 1.5 mL H2O, and 0.1 mL of 0.06 M 1.5 DL-isocitrate;

- for malate dehydrogenase (MDH): with 2.0 mL of 0.3 M Tris/HCl buffer (pH 8.5), 0.6 mL of 0.01 M MnSO4-H2O, 0.6 mL of 0.002 M NADP, 2.1 mL H2O, and 0.6 mL of 1 M malate;

- for 6-phosphogluconate dehydrogenase (6PGDH) and glucose-6-phosphate dehydrogenase (G6PDH): with 0.5 mL of 0.25 M glycyl-glycine buffer (pH 8.0), 0.5 mL of 0.2 M MgCl2-6H2O, 0.2 mL of 0.0075 M NADP, 0.8 mL H2O, 0.3 mL of 0.01 M 6-phosphogluconate (6PG), and 0.01 M glucose-6-phosphate (G6P).

The mixtures were incubated for 3 min at 30°C, and the absorbance was measured at 340 nm with a photometer (Perkin Elmer, Wellesley, MA, USA) for 5 min. The activity was calculated in IU per g of tissue. For each fat tissue sample, enzyme activities were measured twice and averaged. For further details on protein and enzyme traits see Table 2b.

Fat cell traits

According to the methods described in [12-14], each fat tissue sample was cut up with minimal pressure into slices about 1 mm thick. One g of tissue was suspended in 3 mL KRB buffer (Krebs-Ringer bicarbonate buffer with 5 mM glucose and 25 mM HEPES, pH 7.4) containing 3 mg/mL collagenase and slowly stirred at 37°C for 1 h.

The prepared cell suspension was filtered (PP filter, 1000 μm, Sartorius, Göttingen, Germany), collected in 3 mL KRB buffer, sedimented and again suspended in 3 mL KRB buffer. Then, 500 μL cell suspension were incubated with 5 mL collidine-HCl buffer (1 M 2,4,6-trimethylpyridine, 0.1 M HCl, 0.26 M NaCl, pH 7.4) and 3 mL OsO4 solution (3% w/v OsO4 in collidine-HCl buffer) for 24 h at room temperature. The number of suspended cells was measured with a Coulter-Counter (Model TA II, Beckman, Krefeld, Germany) in different size fractions. In practise, the particle counter measured the changes of resistance caused by individual particles passing the opening of a capillary wall with electrodes on both sides. Using an automatic coincidence correction guarantied that particles passing simultaneously were counted separately. Assuming spherical particles, the particle numbers and volumes were calculated for size classes with cell diameters of 23, 29, 36, 41, 57, 73, 92, 114, 146, 183, and 228 μm.

Marker loci and genotyping

Marker loci were selected to be informative, evenly distributed over the chromosomes, and nearly the same for the three crosses. Only when the information content of a selected locus within a cross was low, was an alternative flanking locus chosen for that cross. For regions with previously detected QTL for performance traits [2] on SSC2, SSC4 and SSCX, high marker density maps were built. Per cross, 216 (M × P), 169 (W × P) and 195 (W × M) polymorphic markers were genotyped (Table 3). Marker loci parameters (map position, number of alleles, observed informative meioses etc.) and polymorphism types are provided in Additional file 1.

Table 3. Overview of marker loci and chromosomesa

Additional file 1. Markers used for linkage and QTL analysis. The used marker loci are shown together with literature references and positions on the USDA MARC map. Moreover, the map positions, numbers of alleles and numbers of informative meioses are listed for each of the three crosses.

Format: DOCX Size: 77KB Download fileOpen Data

Statistical analyses

Linkage mapping of marker loci and calculation of trait values

Linkage mapping was performed using the CriMap software, version 2.4 [15,16]. The information content of each locus for mapping was assessed by the number of informative meioses (Additional file 1). The number of informative meioses averaged across all loci was 558 (702) for the M × P cross, 520 (722) for the W × P cross and 623 (732) for the W × M cross, the number in brackets being the maximum number of informative meioses for a locus. The frequencies of the observed informative meioses per cross were 0.79 (M × P), 0.72 (W × P) and 0.85 (W × M).

Additional file 2 contains the numbers of observations, phenotypic means, standard deviations and determination coefficients of the traits for the F2 animals of each cross.

Additional file 2. Parameters of the traits. Numbers of observations, phenotypic means, standard deviations and determination coefficients are given for the considered traits of the F2 animals and each cross.

Format: DOCX Size: 20KB Download fileOpen Data

QTL analysis

The least square method was applied for QTL mapping [17] and was performed separately for each of the three crosses in steps of 1 cM; the steps were reduced when the distance between marker loci was shorter. As described for the autosomes in [3] and for chromosome X in [18], the conditional probabilities for the transfer of an allele from the founder to the F2 individual were calculated for any position of the linkage array by considering all marker loci of a linkage group simultaneously and stored as additive and dominant components. From these linear components, the additive and dominant effects were calculated for each trait in a generalized linear model procedure (GLM) including the continuous (age at slaughter) and discontinuous (two-month classes of seasonal influence, sex, litter number) independent variables. Only 91 W × M F2 animals were measured for fat cell traits, which were not adjusted for the effects of season and litter number in our models because of insufficient connectedness of these independent variables. The mean square estimates of the additive and dominant components in relation to the error variance was calculated from the complete model, and the position on a chromosome with the highest F ratio value was considered as the most likely QTL position. Genome-wide (P < 0.05) significant QTL maxima (major peaks) were determined for all traits (Table 4). Additional genome-wide significant minor peaks were registered per trait and chromosome with P < 0.01 for performance traits (Table 2a) and P < 0.05 for the other traits (Table 2b and 2c) when they were more than 20 cM away from the major peak and from the already considered minor peaks.

Table 4. Genome-wide significant QTL for fat related traits identified in the three Hohenheim crosses

For chromosomes SSC6, 7 and X, we performed separate calculations for female and male animals in order to test sex-specific differences in QTL positions and genetic effects. The model for these data sets includes all independent variables, with the exception of sex.

Threshold values of the test statistic were derived by permutation tests [19], using 1000 repetitions. All permutations were calculated for different traits in data sets for crosses and chromosomes separately. Applying a Bonferroni correction [20], the P < 0.01 and P < 0.05 genome-wide thresholds were calculated for chromosomes 7, 16 and × and then averaged across the chromosomes and crosses, since the thresholds between the crosses and traits showed only slight differences (Additional file 3).

Additional file 3. Genome-wide threshold values. The threshold values, which were calculated according to [19] and with 1000 permutations, are listed for the P < 0.05 and P < 0.01 significance levels.

Format: DOCX Size: 12KB Download fileOpen Data

Testing multifactorial influences on selected traits, the additive and dominant components of significant QTL positions detected across all the chromosomes were included stepwise by using a multiple position model which included the environmental variables. Components with a significant proportion of the explained variance remained in the final model (see results in Table 5).

Table 5. Combined analysis of significant QTL positions a

Results and Discussion

Genome-wide distribution of QTL

Within each cross, we identified QTL which explain more than about 4.3% of the error variance (VF2) with a P < 0.05 genome-wide significance level (threshold with F ratio > 8.5). As shown in Table 4, a total of 147 genome-wide QTL were found (76 at P < 0.05, and 71 at P < 0.01) for the three crosses. The numbers of significant QTL were 30 at P < 0.05 and 33 at P < 0.01 for M × P, 22 at P < 0.05 and 25 at P < 0.01 for W × P, and 24 at P < 0.05 and 13 at P < 0.01 for W × M. However, since we tested three populations and 40 traits in 120 genome scans, about six false positive QTL may occur at a genome-wide 5% significance.

The numbers of QTL detected per trait were about three times higher for the performance traits (Table 2a) than for the other groups of traits (protein, enzyme, fat cell traits, Table 2b and 2c). This finding can be explained by the fact that performance traits are likely to be influenced by a higher number of genes than protein, enzyme and fat cell traits.

In Table 4, the QTL positions and the flanking marker loci for the Hohenheim maps are indicated together with the corresponding USDA MARC map positions. Figure 1 shows the genome-wide QTL distribution for the three crosses. For performance traits, if only the major QTL and adjusted positions on USDA MARC map are considered, the following results can be emphasized:

thumbnailFigure 1. Genomic distribution of QTL. The distribution of the QTL detected in the Hohenheim crosses (M: Meishan; P: Pietrain; W: European wild boar) and with F ratio values above the genome-wide thresholds P = 0.05 is shown on the pig chromosomes (SSC); for each cross, the sex-averaged map in Kosambi morgan (M) is adjusted to the length calculated for the Hohenheim M × P cross; results for SSCX were obtained from female animals; the different symbols for the three trait groups represent major QTL peaks (black) and minor QTL peaks (red) that show distances > 20 cM to the major peak and to other minor peak observed for the same trait.

An accumulation of QTL for fat deposition traits (performance traits) was observed on SSC1. For the W × P cross, QTL were mainly located at positions 76-78 and 87-90 cM. QTL at positions 89-91 and 105-108 cM were detected in the W × M cross, besides two other QTL at positions 57 cM and 113 cM. QTL at 114 and 136 cM were observed in the M × P cross. A QTL for enzyme activity was found with a 5% significance level in the W × P cross, and several QTL were detected in W × P and W × M crosses for fat cell parameters at about 91, 104 and 111-113 cM, three of them near SW705, where [21] has detected QTL for fat cell traits.

On SSC2, only QTL related to performance traits were found in the W × P cross (at 57 cM and 73 cM) in spite of the fact that in the Pietrain breed, the allele IGF2-intron3-3072A responsible for a paternally expressed QTL at the proximal end (0.6 cM) of SSC2 affecting muscle growth and fat deposition is nearly fixed, while in wild boar and the Meishan breed only the wild allele IGF2-intron3-3072G is detected [22]. Therefore, F1 males from W × P and M × P crosses should be IGF2 heterozygous and about half of the F2 animals should possess the allele IGF2-intron3-3072A. The IGF2-intron3-3072 locus was not tested in the crosses as no suitable assay was available. However, its location corresponds to the interval between the markers SW2443/SWC9 and S0141, in which no QTL for performance traits was observed in this study.

Two QTL (P < 0.01) were detected on SSC4, one related to performance traits (37 cM, M × P cross) and one to fat cell traits (74 cM, W × P cross). Another QTL for fat cell traits was found at position 53-55 cM (W × P cross).

Several QTL for performance traits were also found on SSC6 in the M × P cross between the markers TGFB1 and NPPB at around 76 cM. The QTL for both traits on SSC6 in the W × P cross were located in the same interval. Whereas Bidanel et al. [23] have confirmed this QTL position, other authors [24,25] have mapped a QTL for back fat thickness on SSC6 in the vicinity of SW1881 corresponding to position 121 cM.

All 20 QTL (P < 0.01) on SSC7 were found in the major histocompatibility complex (MHC), of which 19 were located approximately 10 cM around the genes TNFA and TNFB. These 19 QTL seem to be distributed in three clusters, one slightly proximal to marker KE6, one slightly distal to TNFA/TNFB and one about 6 cM distal to TNFA/TNFB. The remaining QTL (performance trait AFW, M × P cross) was detected about 9 cM distal to TNFA/TNFB. A total of 18 QTL was observed in the M × P cross for performance (4), enzyme activity (5) and fat cell traits (9), and only two QTL were detected in the W × M cross (one for performance and one for enzyme activity traits). These differences of QTL between crosses might be affected by the information content of marker loci. The QTL for back fat thickness located near TNFA/TNFB have also been reported by [26-29] and Miller [21] has located QTL for fat cell traits at the same position.

On SSC8, 12, 13, 14, 16 and 18, several QTL for traits related to protein content and fat cells were observed, three of them with P < 0.01. Amongst these, two concerning fat cell traits were found on SSC8 for the W × M cross at 108 cM (calculated from 91 observations only), and one QTL detected on SSC14 for another fat cell trait was located between the markers SW210 and SW2488 in the W × P cross.

QTL for protein content were detected on SSCX for the W × M cross at 80-81 cM in the immediate vicinity of CAPN6. QTL related to enzyme activities were found on SSCX in the M × P cross at positions 29, 57 and 112-114 cM. Another QTL for fat cell traits was found at about 56 cM, at the same position where [30] described a QTL for backfat thickness.

Effects of F2 crosses on QTL profiles

As shown in Figure 1 and Table 4, most of the QTL were observed within a few chromosome regions only, and the QTL were often specific to one or two of the three F2 crosses. For example, QTL on SSCX occur mainly in the crosses M × P and W × M and with a cross-specific distribution. The QTL detected in similar chromosomal intervals in two of the three crosses indicate that alleles transmitted from one of the resource groups are different from the alleles in the two other resources.

High allelic effects caused by a distinct founder breed were observed, for example, on SSC4 (near ATP1A2), SSC6 (near RYR1) and SSC7 (between TNFA and S0102). The relevant SSC7 interval includes the MHC region where Meishan cryptic alleles are responsible for a decrease in fat deposition and enzyme activity traits and an increase in the proportion of small fat cells' numbers and volumes (observed in the F2 M × P and W × M crosses). The same effects of Meishan alleles on SSC7 have been reported for fat deposition as well as for numbers and volumes of adipocytes in a Large White × Meishan backcross [31]. On the contrary, Meishan alleles that increase fat deposition were located in the M × P and W × M crosses on SSC1 between TGFBR1 and SW705. Moreover, Pietrain alleles in the crosses with Meishan as well as with wild boar on SSC6 at TGFB1/A1BG had negative effects on obesity. None of the regions with significant effects on fat deposition traits was common to all three crosses, except the one for fat cell traits between TNFB and PSMA4 on SSC7 at about 55 to 90 cM referring to the USDA MARC map.

Figure 2 demonstrates the cross-specific QTL profiles for SSC1, SSC7 and SSCX. The QTL for protein content on SSCX at CAPN6 (mapped at 81 cM on USDA MARC map, [18,32]) was observed only in the W × M cross. Numerous QTL profiles on SSC1 and SSC7 were similar between the M × P and W × M crosses indicating that allele effects in Meishan were highly different to those in Pietrain and wild boar. However, SSC7 QTL were similar among all three crosses for an interval between about 50 and 100 cM (which contains the MHC, see Figure 2), revealing that major QTL effects are caused by alleles that segregate in all the crosses.

thumbnailFigure 2. Examples of F ratio profiles in the different Hohenheim crosses shown for chromosomes 1, 7 and X. The solid line indicates the P = 0.01, the dashed line the P = 0.05 genome-wide thresholds, and the dotted line the P = 0.05 chromosome-wide threshold for F ratio values; traits are defined in Table 2; positions of markers are given in Kosambi centimorgan (cM) on the linkage maps of Hohenheim crosses; results for SSCX were obtained from female animals; markers are described in Additional file 1; data sets for the two sexes are shown with the averaged linkage map distances whereas for SSCX the female map distances are used; M: Meishan; P: Pietrain; W: European wild boar.

Several of the F ratio profiles reveal multiple peaks per chromosome (Figure 2, Table 4). This might be caused by pleiotropic effects of the involved genes. However, multiple peaks in the F ratio profile of a trait per chromosome may also result from linkage disequilibria among alleles of linked loci in F2 animals, whereby the linkage disequilibrium increases while the distances between the considered loci decrease. Significances of QTL peaks can be influenced by different information contents of the marker loci used in the flanking regions of a QTL. Thus, more markers and multipoint regression analyses may help to determine the contribution of single QTL peaks to the total genetic variance of the trait considered.

Examples of multiple and cross-specific QTL peaks per chromosome are also shown in Figure 2 for SSC7. In the M × P and W × M crosses, the major QTL profiles on SSC7 span from about 55 to 90 cM (including the genes CYP21A2, KE6, TNFA, TNFB), and in the W × P cross the major QTL were found at about 105 cM (between S0066 and S0115). The 30 cM interval covering the largest QTL on SSC7 contains the MHC known to include numerous functional genes in man and mouse. In this interval, genome-wide significant QTL were detected especially in both Meishan crosses. Concerning fat deposition traits, this could be due mainly to a smaller difference between the purebred estimates for wild boar and Pietrain compared to that between these two breeds and the Meishan breed [10]. For instance, the difference in average back fat depth (ABFD) between Pietrain and wild boar was 2.13 mm, whereas it was 7.77 mm between Pietrain and Meishan and 9.90 mm between wild boar and Meishan. A further example of effects of crosses on the patterns and positions of QTL was observed for SSC6 in the region of the loci LIPE, TGFB1, A1BG, EAH and NPPB (USDA MARC map 75 to 80 cM, Table 4). Important QTL were detected in this region for both M × P and W × P crosses. The additive effects for the grand-paternal inheritance indicate a negative influence of distinct Pietrain founder alleles on performance traits associated with fatness.

Differences of QTL profiles calculated separately for female and male F2 offspring

QTL analyses for female and male F2 offspring are shown for example on SSC6, SSC7 and SSCX and use averaged linkage map distances for the autosomes and the female map distances for SSCX. Figure 3 shows QTL effects for the traits FVCM, FMR and PCOI, which differ between female and male F2 animals. For example, the trait FVCM in females of the M × P cross are highly influenced by a QTL at 75 cM on SSC7 whereas the F ratio value for males shows non-significance at that position. Males of the W × P cross show a QTL for the trait FMR on SSC6 at a position near 125 cM, which is located about 40 cM distal to the position found for all (female and male) animals. The trait PCOI represents an example of sex specific QTL positions on SSCX (QTL at 94 and 104 cM for females and males, respectively, W × M cross).

thumbnailFigure 3. Examples of F ratio profiles calculated for all (females and males), female or male F2 animals. For further explanations see Figure 2.

Sex specific QTL positions have also been reported on SSCX for muscle, fatness and growth related traits in the W × M cross [18]. Sex specific and fat related QTL have been described on chromosome 5 in chicken [33] and on several chromosomes in mouse [34]. Gene expression studies in male and female F2 mice have shown a large degree of sexually dimorphic gene expression in several tissues [35,36]. An expression QTL (eQTL) study [37] has shown that most of the eQTL were cis eQTL (mapping to the location of the gene) and sex-shared. Genetic mechanisms possibly underlying sex-specific expression, like sex linkage, sex-specific allelic effects or genomic imprinting, are discussed in [38].

Combined analysis of significant QTL positions across chromosomes

Across all the pig chromosomes and for selected traits, we have carried out a combined analysis of the additive and dominant components of significant QTL positions. Taking each trait separately, the components of those positions were included step by step in a multiple position model. In the final model, only components with significant variance proportions were included. Examples of the results are shown in Table 5 and elucidate why the explained phenotypic variance in the F2 generations increased markedly up to about 23%, and the determination coefficients (r2) of the initial model (analyses without genetic independent variables) were more or less doubled. For each trait, several QTL positions, partially located on the same chromosome, remained significant in the combined analysis. This means that the combined analysis indicates a predominant contribution of a few QTL regions to the genetic variance of a trait. Therefore, multiple testing elucidates chromosome intervals which can be significant for breeding programmes.


As demonstrated in this report, in pig, fat related traits correspond to numerous specific QTL positions across the genome. For some of the traits associated with specific enzymes, protein content and cell structure in fat tissue, it is the first time that they are included in a QTL analysis. We have found that QTL positions differ between F2 crosses, and differ partly for their additive and dominant effects. Some of these QTL show sex specific effects. Many of the QTL profiles possess multiple peaks especially in regions with a high marker density, and confidence intervals mostly exceed 10 cM [39]. Therefore, QTL intervals are rarely narrowed down to a sufficiently small number of candidate loci to be able to suggest one as the most probable causative gene.

Nevertheless, porcine chromosome regions, which contain QTL, can be aligned with loci of expressed genes, as well as with orthologous genes in man and mouse using data from PigQTLdb ([4]; webcite). Today, QTL intervals can be compared with the pig genome sequence data (Sscrofa9, Wellcome Trust Sanger Institute 2009, webcite) to investigate the action of single genes and their variants. The selection of putative causative genes may consider groups of genes that are regulated in parallel and are members of the same metabolic pathway. Thus, the results of genome-wide QTL mapping are important for subsequent analyses of specific genes as well as for selecting DNA markers for breeding purposes.

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

HG is responsible for most of the concept and design, for finding funding, and for drafting the tables and manuscript. SC, AS and SP have carried out the genotyping of marker loci and revised the manuscript. HB performed the statistical analysis, created the figures and helped to draft the manuscript.

All authors have read and approved the final manuscript.


The investigation was supported by the German Research Foundation (DFG, grant nos. Mu616/6 and Ge291/20), the EC programmes BRIDGE and INCO-Copernicus (Contract no. ERBIC15CT960902), the Czech Science Foundation (Grant No. 523/07/0353 and 523/06/1302), and the Institutional Research Plan of the IAPG AS CR (AV0Z50450515). The Meishan pigs used in the experiments originated from a population provided by the Wageningen Agricultural University and Euribrid, BV Boxmeer, The Netherlands.


  1. Bidanel JP, Rothschild M: Current status of quantitative trait locus mapping in pigs.

    PigNews Inform 2002, 23:39N-54N. OpenURL

  2. Geldermann H, Müller E, Moser G, Reiner G, Bartenschlager H, Cepica S, Stratil A, Kuryl J, Moran C, Davoli R: Genome-wide linkage and QTL mapping in porcine F2 families generated from Pietrain, Meishan and Wild Boar crosses.

    J Anim Breed Genet 2003, 120:363-393. Publisher Full Text OpenURL

  3. Hu ZL, Dracheva S, Jang W, Maglott D, Bastiaansen J, Rothschild MF, Reecy JM: A QTL resource and comparison tool for pigs: PigQTLDB.

    Mamm Genome 2005, 16:792-800. PubMed Abstract | Publisher Full Text OpenURL

  4. Hu ZL, Fritz ER, Reecy JM: AnimalQTLdb: a livestock QTL database tool set for positional QTL information mining and beyond.

    Nucleic Acids Res 2007, 35:D604-D609. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  5. Müller E, Rogdakis E: Genetische Regulation des Fettstoffwechsels beim Schwein.

    Hohenheimer Arbeiten, Reihe Tierische Produktion 1985, 131:7-28. OpenURL

  6. Strutz C: Phenotypic and genetic parameters of NADPH-generating enzymes in porcine adipose tissue. In Ph D thesis. University of Hohenheim; 1977. OpenURL

  7. Rogdakis E, Strutz C: Phänotypische und genetische Parameter für den löslichen Proteingehalt im Schweinegewebe.

    Züchtungskunde 1978, 50:312-316. OpenURL

  8. Etherton TD: Subcutaneous adipose tissue cellularity of swine with different propensities for adipose tissue growth.

    Growth 1980, 44:182-191. PubMed Abstract OpenURL

  9. Fiedler I, Wiesemuller W, Michelchen G, Kuhn G: Fat clogging, fat cell count and fat cell size in swine in relation to age and feeding intensity.

    Arch Tierernahr 1990, 40:681-687. PubMed Abstract OpenURL

  10. Müller E, Moser G, Bartenschlager H, Geldermann H: Trait values of growth, carcass and meat quality in Wild Boar, Meishan and Pietrain pigs as well as their crossbred generations.

    J Anim Breed Genet 2000, 117:189-202. Publisher Full Text OpenURL

  11. Lowry OH, Rosebrough NJ, Farr AL, Randall RJ: Protein measurement with the Folin phenol reagent.

    J Biol Chem 1951, 193:265-275. PubMed Abstract | Publisher Full Text OpenURL

  12. Hirsch J, Gallian E: Methods for the determination of adipose cell size in man and animals.

    J Lipid Res 1968, 9:110-119. PubMed Abstract | Publisher Full Text OpenURL

  13. Etherton TD, Thompson EH, Allen CE: Improved techniques for studies of adipocyte cellularity and metabolism.

    Lipid Res 1977, 18:552-557. OpenURL

  14. Kalbitz S: Untersuchungen der Insulinrezeptoren und des Glukosemetabolismus im Rückenspeck von Schweinen verschiedener Zuchtlinien. In Ph D thesis. University of Hohenheim; 1986. OpenURL

  15. Green P, Falls K, Crooks S: Documentation for CRI-MAP, Version 2.4. Washington University School of Medicine, St. Louis, MO; 1990. OpenURL

  16. Keats BJ, Sherman SL, Morton NE, Robson EB, Buetow KH, Cartwright PE, Chakravarti A, Francke U, Green PP, Ott J: Guidelines for human linkage maps: an international system for human linkage maps (ISLM, 1990).

    Genomics 1991, 9:557-560. PubMed Abstract | Publisher Full Text OpenURL

  17. Haley CS, Knott SA, Elsen JM: Mapping quantitative trait loci in crosses between outbred lines using least squares.

    Genetics 1994, 136:1195-1207. PubMed Abstract | PubMed Central Full Text OpenURL

  18. Cepica S, Bartenschlager H, Geldermann H: Mapping of QTL on chromosome X for fat deposition, muscling and growth traits in a wild boar × Meishan F2 family using a high-density gene map.

    Anim Genet 2007, 38:634-638. PubMed Abstract | Publisher Full Text OpenURL

  19. Churchill GA, Doerge RW: Empirical threshold values for quantitative trait mapping.

    Genetics 1994, 138:963-971. PubMed Abstract | PubMed Central Full Text OpenURL

  20. Miller RG: Simultaneous Statistical Interference. 2nd edition. Springer Verlag, Berlin; 1981:6-8. OpenURL

  21. Liu LT, Ren J, Li L, Li WB, Yang B, Zhu WC, Guo BL, Ma JW, Huang LS: A genome scan for loci affecting adipocyte size and number in abdominal fat in a White Duroc × Erhualian F2 resource population.

    Animal 2009, 3:632-636. Publisher Full Text OpenURL

  22. Van Laere A-S, 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. PubMed Abstract | Publisher Full Text OpenURL

  23. 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. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  24. Ovilo C, Perez-Enciso M, Barragan C, Clop A, Rodriguez C, Oliver MA, Toro MA, Noguera JL: A QTL for intramuscular fat and backfat thickness is located on porcine chromosome 6.

    Mamm Genome 2000, 11:344-346. PubMed Abstract | Publisher Full Text OpenURL

  25. Edwards DB, Ernst CW, Tempelman RJ, Rosa GJ, Raney NE, Hoge MD, Bates RO: Quantitative trait loci mapping in an F2 Duroc × Pietrain resource population: I. Growth traits.

    J Anim Sci 2008, 86:241-253. PubMed Abstract | Publisher Full Text OpenURL

  26. 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. PubMed Abstract | PubMed Central Full Text OpenURL

  27. Rohrer GA: Identification of quantitative trait loci affecting birth characters and accumulation of backfat and weight in a Meishan-White Composite resource population.

    J Anim Sci 2000, 78:2547-2553. PubMed Abstract | Publisher Full Text OpenURL

  28. Malek M, Dekkers JC, Lee HK, Baas TJ, Rothschild MF: A molecular genome scan analysis to identify chromosomal regions influencing economic traits in the pig. I. Growth and body composition.

    Mamm Genome 2001, 12:630-636. PubMed Abstract | Publisher Full Text OpenURL

  29. Kim KS, Thomsen H, Bastiaansen J, Nguyen NT, Dekkers JC, Plastow GS, Rothschild MF: Investigation of obesity candidate genes on porcine fat deposition quantitative trait loci regions.

    Obes Res 2004, 12:1981-1994. PubMed Abstract | Publisher Full Text OpenURL

  30. Harlizius B, Rattink AP, De Koning DJ, Faivre M, Joosten RG, Van Arendonk JA, Groenen MA: The X chromosome harbors quantitative trait loci for backfat thickness and intramuscular fat content in pigs.

    Mamm Genome 2000, 11:800-802. PubMed Abstract | Publisher Full Text OpenURL

  31. Demars J, Riquet J, Sanchez MP, Billon Y, Hocquette JF, Lebret B, Iannuccelli N, Bidanel JP, Milan D, Gondret F: Metabolic and histochemical characteristics of fat and muscle tissues in homozygous or heterozygous pigs for the body composition QTL located on chromosome 7.

    Physiol Genomics 2007, 30:232-241. PubMed Abstract | Publisher Full Text OpenURL

  32. Cepica S, Masopust M, Knoll A, Bartenschlager H, Yerle M, Rohrer GA, Geldermann H: Linkage and RH mapping of 10 genes to a QTL region for fatness and muscling traits on pig chromosome X.

    Anim Genet 2006, 37:603-604. PubMed Abstract | Publisher Full Text OpenURL

  33. Abasht B, Pitel F, Lagarrigue S, Le Bihan-Duval E, Le Roy P, Demeure O, Vignoles F, Simon J, Cogburn L, Aggrey S, Vignal A, Douaire M: Fatness QTL on chicken chromosome 5 and interaction with sex.

    Genet Sel Evol 2006, 38:297-311. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  34. Taylor BA, Tarantino LM, Phillips SJ: Gender-influenced obesity QTLs identified in a cross involving the KK type II diabetes-prone mouse strain.

    Mamm Genome 1999, 10:963-968. PubMed Abstract | Publisher Full Text OpenURL

  35. Wang S, Yehya N, Schadt EE, Wang H, Drake TA, Lusis AJ: Genetic and genomic analysis of a fat mass trait with complex inheritance reveals marked sex specificity.

    PLoS Genet 2006, 2:e15. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  36. Yang X, Schadt EE, Wang S, Wang H, Arnold AP, Ingram-Drake L, Drake TA, Lusis AJ: Tissue-specific expression and regulation of sexually dimorphic genes in mice.

    Genome Res 2006, 16:995-1004. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  37. Bhasin JM, Chakrabarti E, Peng DQ, Kulkarni A, Chen X, Smith JD: Sex specific gene regulation and expression QTLs in mouse macrophages from a strain intercross.

    PLoS One 2008, 3:e1435. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  38. Bonduriansky R, Chenoweth SF: Intralocus sexual conflict.

    Trends Ecol Evol 2009, 24:280-288. PubMed Abstract | Publisher Full Text OpenURL

  39. Darvasi A, Soller M: A simple method to calculate resolving power and confidence interval of QTL map location.

    Behav Genet 1997, 27:125-132. PubMed Abstract | Publisher Full Text OpenURL

  40. Herbst K: Entwicklung, Stand und Perspektiven der Schweineproduktion in der Bundesrepublik Deutschland.

    Züchtungskunde 1980, 52:304-323. OpenURL