Original Research ARTICLE
An F2 pig resource population as a model for genetic studies of obesity and obesity-related diseases in humans: design and genetic parameters
- Department of Veterinary Clinical and Animal Sciences, Faculty of Health and Medical Sciences, University of Copenhagen, Copenhagen, Denmark
Obesity is a rising worldwide public health problem. Difficulties to precisely measure various obesity traits and the genetic heterogeneity in human have been major impediments to completely disentangle genetic factors causing obesity. The pig is a relevant model for studying human obesity and obesity-related (OOR) traits. Using founder breeds divergent with respect to obesity traits we have created an F2 pig resource population (454 pigs), which has been intensively phenotyped for 36 OOR traits. The main rationale for our study is to characterize the genetic architecture of OOR traits in the F2 pig design, by estimating heritabilities, genetic, and phenotypic correlations using mixed- and multi-trait BLUP animal models. Our analyses revealed high coefficients of variation (15–42%) and moderate to high heritabilities (0.22–0.81) in fatness traits, showing large phenotypic and genetic variation in the F2 population, respectively. This fulfills the purpose of creating a resource population divergent for OOR traits. Strong genetic correlations were found between weight and lean mass at dual-energy x-ray absorptiometry scanning (0.56–0.97). Weight and conformation also showed strong genetic correlations with slaughter traits (e.g., rg between abdominal circumference and leaf fat at slaughtering: 0.66). Genetic correlations between fat-related traits and the glucose level vary between 0.35 and 0.74 and show a strong correlation between adipose tissue and impaired glucose metabolism. Our power calculations showed a minimum of 80% power for QTL detection for all phenotypes. We revealed genetic correlations at population level, for the first time, for several difficult to measure and novel OOR traits and diseases. The results underpin the potential of the established F2 pig resource population for further genomic, systems genetics, and functional investigations to unravel the genetic background of OOR traits.
Obesity is a growing worldwide public health problem with a prevalence that has almost doubled between 1980 and 2008 (World Health Organization, 2012), resulting in one of the most significant public health challenges of the twenty-first century. In 2008, more than 1.4 billion adults were overweight and more than 0.5 billion were obese. Obesity has been defined by O’Rahilly (2009) as “a state in which the total amount of triglyceride stored in adipose tissue is abnormally increased.” As a result, several other organs and important metabolic systems are dysregulated in obese individuals (Goossens, 2008). Obesity is associated with several severe diseases such as type 2 diabetes and cardiovascular diseases, which are all closely related to each other and resulting from the interaction between genetic and environmental factors (Bener et al., 2005). The recent increase in the worldwide prevalence of obesity has increased the need for knowledge about the genetic and biological background of obesity and related metabolic diseases. A deeper insight into the genetic background and, in turn, novel metabolic molecular pathways will help to explain why some people are susceptible to obesity whereas others remain lean. Consequently, this could lead to improvement of prevention, therapy, and treatment of obesity and obesity-related (OOR) diseases (Farooqi and O’Rahilly, 2007).
Many years of research has revealed the difficulty in pinpointing the genetic background of OOR traits and diseases. Several genome-wide association studies (GWAS) have been performed to identify common genetic variants contributing to OOR traits and diseases. Using a GWAS-approach, 6–11% of the genetic variance could be explained in body mass index (BMI) by the identification of 32 loci, while heritability was estimated at 40–70% (Speliotes et al., 2010) and 25–30% of the genetic variance could be explained in blood lipids by the identification of 95 loci (Teslovich et al., 2010). For type 2 diabetes at least 10 loci were identified using GWAS strategies (Scott et al., 2007; Sladek et al., 2007), but their contribution to individual risk prediction of type 2 diabetes is relatively small. In human studies of these traits, it has been shown to be extremely difficult to pinpoint the most important biological pathways involved in complex traits due to confounding between genetics and environment and due to genetic heterogeneity. Moreover, precise measurements on various OOR traits in human research are often difficult, very expensive, and in some cases even impossible, resulting in the inability to detect the total genetic variance. A complex trait such as obesity is a combination of numerous phenotypic conditions. By dissecting them into endophenotypes, the phenotype can be more accurately defined which will improve the genetic detection of diseases (Bougnères, 2002). An animal model will not only help reduce costs substantially but also helps control the environment. It has the potential to create a large dataset of various unique phenotypic measurements for genetic research to increase knowledge about the genetic and biological background of OOR traits.
Expression studies of several genes involved in glucose metabolism in adipose tissues have produced markedly different results in humans and rodents hindering the translation from findings in rodent models to human biology (Arner, 2005). The pig resembles humans more than rodents, as the pig shares metabolic and digestive features as well as cardiovascular systems and other anatomical and physiological features with humans (Spurlock and Gabler, 2008). The recent annotation of the porcine genome (Groenen et al., 2012) as well as QTL and selective sweep analysis (Rubin et al., 2012) has demonstrated that homologous variants are present in human and pigs. Moreover, it has been shown that human obesity genes, such as FTO and MC4R, have homologous variants in pigs (Kim et al., 2000; Fontanesi et al., 2009). As fatness in pigs is an economically important trait in pig production, several studies have focused on the genetics of fatness in commercial pig lines. A large study on fatness traits in pigs uncovered 663 genes with expression profiles that were significantly correlated with the trait “fat area” (in cm2, measured between the 13th and 14th rib) as well as biologically relevant pathways and networks (Ponsuksili et al., 2011). Several QTL and GWAS studies have been conducted using pigs from the same breeds used as founders in the studied population. In the Yorkshire pigs, e.g., several micro-satellite markers were significantly associated with average daily gain and backfat thickness (Kim et al., 2006). In a Duroc × Pietrain F2 population, several QTLs were significantly associated with serum lipids, including a large QTL on SSC1 (Uddin et al., 2011). This has also been reported to be in association with average daily gain and carcass composition in the same population (Liu et al., 2007). Recently, a region on SSC6 for backfat thickness was detected in Duroc pigs (Okumura et al., 2012) and in commercial pigs several QTL regions were detected in association with body composition traits in commercial pigs (Large White and Large White × Landrace crosses) (Fan et al., 2011). As per the review of Bellinger et al. (2006) pigs are a useful model for type 2 diabetes, as they show type 2 diabetes related characteristics. From this review, it can be concluded that both the Göttingen minipig and the Yorkshire pigs are valuable sources to investigate the genetic and biological background of type 2 diabetes.
Another great advantage of using the pig as a model in genetic research is the possibility to create a population segregating for the traits of interest. That is, by mating genetically and phenotypically divergent individuals and subsequently intercrossing the offspring, an informative resource population can be established for genetic studies. An important prerequisite for a successful outcome of exploiting a resource population is the ability to precisely define and measure phenotypes, endophenotypes, and intermediate phenotypes (Kadarmideen et al., 2006) (e.g., transcriptome or proteome) and control the environment and diet in an experimental setting. Such resource populations allow not only genetic (GWAS) studies but also systems genetics and biology studies where genetic variations in intermediate phenotypes can be studied by integrative approaches (Kadarmideen et al., 2011; Kogelman et al., 2011).
An F2 pig resource population was established to increase knowledge about the genetic and biological background of OOR traits because of the above-mentioned arguments. It is well accepted that before conducting quantitative trait loci (QTL) mapping or GWAS in resource populations, a thorough knowledge on (additive) genetic variance of each trait and (genetic) relationships amongst them are investigated, particularly if such parameters for novel traits have not yet been reported elsewhere. For instance, genetic parameters for 16 different trypanotolerance traits were estimated (Van Der Waaij et al., 2003) for detection of QTLs for these trypanotolerance traits in a F2 design (Hanotte et al., 2003). In this paper, we report the establishment of the F2 pig resource population and quantify its genetic make-up by estimating heritabilities of 36 different (and partly novel) OOR traits and genetic and phenotypic correlations amongst them.
Materials and Methods
An F2 population was created using Danish production pig breeds, i.e., purebred Yorkshire (YY) and Duroc (DD) sows from a DanBred breeding herd and Göttingen minipig (MM) boars from Ellegaard A/S in the parental generation. The production pig breeds have been selected for leanness and growth among other traits during the last 60 years, whereas Göttingen minipigs are prone to obesity and bred principally for their small size and ease of handling. They are also known to share the metabolic impairments seen in obese humans (Johansen et al., 2001). For the parental animals two different maternal lines were used to ensure a higher likelihood of detecting narrow QTL regions, as it is expected that many of the same genes are fixed in the different breeds, but the LD structure will differ between them.
Seven female Duroc pigs and seven female Yorkshire pigs were artificially inseminated with 14 male Göttingen minipigs (28 F0 animals), resulting in 14 full-sib F1 families. In total, 127 F1 animals were produced: 60 males and 67 females, from which 81 F1 animals were randomly selected to produce the F2 generation. A total of 29 Minipig-Yorkshire (MY) and 24 Minipig-Duroc (MD) gilts were mated to respectively 15 MY and 13 MD non-litter mate boars (each boar used once or twice) to produce the F2 generation, so all F2 animals were a combination of Minipig × Yorkshire or Minipig × Duroc (Figure 1).
Figure 1. Overview of the experimental design. The resource population consisting of 454 F2 pigs was created and extensively phenotyped using a cross between lean production pigs (Yorkshire/Duroc) × Göttingen minipig.
The F2 animals were produced at the research farm, University of Copenhagen Tåstrup, Denmark. Animal care and maintenance have been conducted according to the Danish “Animal Maintenance Act” (Act 432 dated 09/06/2004). A total of 474 piglets were born in three different groups, specified by season of birth (SOB), whereby the first SOB consisted of approximately 250 piglets born in June or July, the second SOB consisted of approximately 125 piglets born in August, September, or October, and the third SOB consisted of approximately 100 piglets born in December or January. Litters in the F2 generation consisted on average of eight piglets per litter and ranged from 2 to 12 born piglets. Piglets were weaned at 6–8 weeks of age. In total, the F2 generation consisted of 474 born piglets, 469 of which were still alive at 5 weeks of age and used in further analyses. The complete pedigree consists of 563 animals.
Collection of Phenotypes
Intensive phenotypic recording was performed on all pigs from birth to slaughter. Body weight was measured at birth and at the following approximate time points: 2 weeks (13–16 days), 5 weeks (34–38 days), 2 months (64 ± 11 days), 125 days (115 ± 27 days), and 7 months (220 ± 45 days). After weaning, at approximately 2 months of age (64 ± 11 days), the piglets were anesthetized with 1 ml/10 kg of a mixture containing zolazepam and tiletamin (Zoletil® 50 Vet., ChemVet, Denmark), ketamin (Ketaminol® Vet. 100 mg/ml, Intervet, Denmark), detomidin (Domitor 1 mg/ml, Pfizer Animal Health, Denmark), and buthorphanol (Torbugesic®, 10 mg/ml, Scanvet, Denmark). This full sedation is approved and granted by the Experimental Animal Inspectorate on July 1 2011 (approval number2007/561-1434). Body composition was determined using dual-energy x-ray absorptiometry (DXA) scanning (Hologic Explorer, Santax Medico, Aarhus, Denmark). The scanning was performed with the pigs in ventral recumbency with the front and hind limbs stretched backwards (Figure 2). Using the differential attenuation of low (38 keV) and high-energy (70-keV) x-rays the body composition of fat and other soft tissues was estimated. The analysis of body composition was performed using the scanner software package, providing measurements of fat mass, lean mass, total mass, fat percentage, bone mineral content (BMC), and bone mineral density (BMD). Previous research indicated that those DXA measurements are very accurate in comparison with chemical analysis (Mitchell et al., 1998). Subsequently, blood was collected from the jugular vein and the following body conformation traits were measured (Figure 3): abdominal circumference (Figure 3A), length of the body from rump to the front of the shoulder blade (Figure 3B), shoulder height (Figure 3C), and thorax height (Figure 3D).
Figure 2. A dual-energy x-ray absorptiometry (DXA) scan of one of the pigs at approximately 2.5 months old.
Figure 3. Body conformation measurements in the pig. Body conformation measurements in the pig including abdominal circumference at the navel (A), the length of the body from rump to the front of the shoulder blade (B), shoulder height (C), and thorax height (D).
At approximately 2.5 months of age, pigs were transported to a production pig farm (Ringsted, Denmark) and placed in pens (20 m2) with an average of 15 piglets of the same sex per pen. Pigs were fed ad libitum with standard production pig feed (Classic, dlg.dk) and with free access to water. All pigs were weighed two to four times within the period from 2.5 months to slaughter at approximately 80–100 kg. At the time of weighing, body conformation measurements were recorded (Figure 3). Fasting glucose was measured on a subset of the pigs (146 animals) fasted for 24 h at the age of 6–8 months. The measurements were performed with an OneTouch UltraSmart Blood Glucose Monitoring System (Lifescan, Inc., Ireland) on blood collected from the jugular vein. The pigs were slaughtered at a commercial slaughterhouse. The pigs were euthanized by bleeding after electrical stunning and the following phenotypes were recorded: weight of leaf fat, weight of the great omentum (272 animals), and weight of mesenteric fat as measured from an 8 cm diameter large section of the omentum in the triangle between ileum and cecum (234 animals). A Fat-o-Meat’er (Carometec A/S, Herlev, Denmark) was used to measure thickness of back fat and meat in millimeters. Measurements were taken between the third and fourth lumbar vertebra (position 1) and between the third and fourth last rib (position 2), both at 8 cm from midline. Meat percentage was calculated as:
where is the meat in percentage, X1 is the backfat measurement at position 1, X2 is the backfat measurement at position 2 and X3 is the muscle thickness at position 2. The Fat-o-Meat’er detects changes in light reflectance to determine shifts between fat and muscle. However, it is known that fat percentage in pigs is positively correlated to the reflectance (Hodgson et al., 1991) resulting in inaccurate measurements in very obese animals. Therefore, Fat-o-Meat’er measurements were set as missing values when reflectance was above 65 (equal to 2.5σ above the mean), thereby excluding inaccurate measurements. Furthermore, carcass weight (excl. skin head and leaf fat) was also recorded.
For subsequent data analysis some variables were calculated using the recorded measurements. Average daily gain was estimated over the whole growth period (from birth to approximately 7 months of age), and as daily gain pre- and post-weaning. In human studies, BMI is a commonly used indicator of obesity (Saunders et al., 2007). It was calculated as weight/(height)2 and therefore in pigs it has been calculated as weight/(length)2, whereby length is the length of the body from rump to the front of the shoulder blade (Figure 3B). Bergman et al. (2011) defined a new parameter to estimate percentage of adiposity more accurately, the body adiposity index (BAI). BAI in humans is proportional to the hip circumference/(height)1.5. Therefore, in this project it has been calculated as abdominal circumference/(length)1.5, where the abdominal circumference is illustrated in Figure 3A.
In total 15F2 animals were removed from the dataset, due to illness or death, which created a F2 resource population of 454 animals.
Statistical Analysis and Power Calculations
Data exploration of all variables was performed to detect outliers, distributions, and relations among variables. Exploration of distribution showed that DXA traits (DXA fat, DXA lean, DXA total, DXA%fat, BMC, and BMD) were not normally distributed and therefore, the first four were transformed using log10 transformation and the last two were transformed using a square-root transformation, respectively, to approach the normal distribution. Quantile–quantile plots indicated that the residuals of transformed variables were approximately normally distributed. All fixed effects were initially analyzed using the GLM procedure of SAS (SAS Version 9.3, 2012) to determine significant fixed effects with respect to all analyzed variables. According to these preliminary analyses, sex, SOB, and age (at time of measurement) were taken into account as fixed effects. For the variables DXA fat, DXA lean, DXA total, SLfat, SLfat_om, and SLfat_int, the length of the animal at the relevant age was also taken into account as fixed effect. The length of the animal is taken into account, because of the difference in conformation of the production pig (long and lean) and the Göttingen minipig (short and thick). Table A1 in Appendix presents the solutions for all fixed effects for the different variables. In addition, in order to detect whether there was sufficient power to detect QTL for the phenotypes under study, we have conducted statistical power/sample size calculations at very conservative assumptions, using qtlDesign R software (Sen et al., 2007). We assumed QTL explains a small (0.5%) to large (15%) proportion of the phenotypic variance at moderate and slightly high heritabilities (of 0.25 and 0.50) with an average distance between SNP markers less than 1 cm. The empirical methods of calculating power for QTL mapping for categorical traits (Kadarmideen et al., 2000) was also adapted for our design to show power on underlying liability scale.
Variance Components Estimation
An univariate animal model was fitted to the data using ASReml (Gilmour et al., 2009) to estimate variance components and heritabilities for each trait considered. The model used for variance component estimations was:
where GENDER is the fixed effect of the sex i of the animal (two classes), SOB is the fixed effect of the SOB j when the pigs were born (three classes), AGE is the age k of the pig in days at the time of measurement as covariate, and bk is the associated fixed regression coefficient. LENGTH is the length of the animal in cm at the time of measuring as a covariate and bl is the associated fixed regression coefficient. The latter effect was only included for fat mass, lean mass, total mass (at DXA scanning) and weight of leaf fat at slaughtering. Random effects included were ANIMAL, the additive genetic effect of the animal l, and e, the error term. The following (co)variances among random effects were fitted:
where A is the additive genetic relationship matrix and I is an identity matrix. The animal and residual terms were assumed to be where is the additive genetic variance and the error variance.
Bivariate Models for Estimation of Genetic and Phenotypic Correlation
In order to estimate heritabilities and genetic and phenotypic correlations, we fitted a series of bivariate animal models for all combinations of traits. The bivariate animal model can be represented as follows:
where y1 and y2 are vectors of observations for the concerning traits, b1 and b2 are vectors of fixed effects for the concerning traits, a1 and a2 are vectors of random animal effects for the concerning traits, X1 and X2 the incidence matrices relating records of the concerning traits to fixed effects, and Z1 and Z2 are incidence matrices relating observations with random animal effects. The expectations and variance are:
where and are additive genetic variances for traits 1 and 2, is the additive genetic covariance between traits 1 and 2, are residual variances for traits 1 and 2, and is the residual covariance between traits 1 and 2. A and I are as described above, while a and e were assumed to be normally distributed with a mean of zero and (co)variances as specified above.
All analyses were implemented using ASReml (Gilmour et al., 2009). Heritabilities, genetic correlations, and phenotypic correlations were calculated from the estimated (co-) variance components. Estimation of heritability with the lowest standard error (univariate or bivariate) was used as heritability for the concerning trait.
Results and Discussion
Descriptive Statistics of the F2 Pig Resource Population
We have created an F2 pig resource population with an extensive phenotypic variation for OOR traits using Göttingen minipigs (genetically predisposed to obesity) and Duroc or Yorkshire pigs (genetically selected for leanness) as founder breeds. This resource population has been comprehensively phenotyped for several OOR traits with the aim to conduct, among other things, QTL mapping, and GWAS. Table 1 shows descriptive statistics of all analyzed traits directly related to OOR, while Table 2 shows descriptive statistics of all analyzed OOR traits indirectly related to OOR in the F2 resource population. Both tables show the mean, SD, range, and coefficient of variation (CV). It is clear from the descriptive statistics that there is substantial phenotypic variation of the traits measured within the pedigree, as illustrated by the high CV values (from 7.5 to 42.4%). In particular, CV values for the fatness traits are very high (f.e., 42.2% for SLfat), confirming that the design of the pedigree has fulfilled the purpose of creating a resource population with extensive phenotypic variation. Bahelka et al. (2007) reported CV’s for different slaughter traits which were all lower than those reported in our F2 pig resource population, e.g., 8.53% for carcass weight, 19.08% for backfat, and 23.05% for fat (kg). It is also important to note that fasting glucose levels (FGLs) ranged from 1.8 to 11.15 mmol/L with a CV of 22.5%, indicating that metabolic parameters involved in, for instance, type 2 diabetes show a wide phenotypic variation in the resource population. Hence, we believe the established F2 resource population is an excellent model suitable for further genetic, genomic, systems genetics, and functional investigations of OOR traits.
Table 1. Summary of all traits directly associated with obesity and obesity-related traits (mean, range, and SD) with their abbreviations.
Table 2. Summary of all traits indirectly associated with obesity and obesity-related traits (mean, range, and SD) with their abbreviations.
Dual-energy x-ray absorptiometry measurements of fat percentage in the F2 resource population showed an average of 18.6% (SE = 2.9). Previous research showed that DXA is an accurate method for estimating body fat, lean mass, and BMC in pork carcasses (Suster et al., 2004). Mitchell et al. (1998) reported fat percentage of 11.43(±3.35) kg, which is slightly lower than the estimate of fat percentage in our population (18.64 ± 2.87 kg). The BMC was higher (338 ± 17) than in our research population (301 ± 94).
Heritabilities are presented in Table 3 and all estimates were statistically highly significant. In general, most estimated heritabilities were high, and higher than previous reports for some traits studied in field populations.
Table 3. Heritabilities and standard errors estimated with the bivariate model using all 563 animals (454 with records).
There was not sufficient information in the data to analyze weight at birth, at 2 weeks of age and at 5 weeks of age. Heritability was high for weight at 2 months (0.78, SE = 0.14) and at 125 days (0.67, SE = 0.17), and slightly lower at 7 months (0.39, SE = 0.12). Huisman et al. (2002) used two different random regression models to estimate heritability of field data of production pigs using three time points to measure: approximately 73, 135, and 190 days. This resulted in heritabilities of 0.14–0.53, 0.17–0.61, and 0.23–0.56, respectively. Johnson and Nugent (2003) estimated heritabilities for weight slightly lower, i.e., between 0.12 and 0.27 at 100 days of age for different commercial production breeds (Landrace, Yorkshire, Duroc, and Hampshire), and between 0.15 and 0.27 at 177 days of age. Heritabilities for body length were estimated for the same breeds between 0.16 and 0.32, while in our population heritability varied between 0.52 ± 0.14 (2 months of age) and 0.26 ± 0.11 (7 months of age). Heritabilities for conformation measurements tended to be lower for measurements at 7 months of age than at 2 months of age, which may be due to a change in the ratio between phenotypic and genetic variance.
The heritability of average daily gain was 0.54 (SE = 0.14), while the heritability of pre-weaning daily gain was 0.81 (SE = 0.14) and post-weaning gain was 0.60 (SE = 0.14). Estimation of heritabilities for daily gain for different production breeds and crosses ranged between 0.20 and 0.31 (Brandt and Täubert, 1998). Kadarmideen et al. (2004) estimated heritability for post-weaning gain in the range of 0.28 for Large White and Landrace breeds, whereas Hofer et al. (1992) estimated heritabilities for total average daily gain ranging between 0.12 and 0.25. Heritabilities for DXA traits were moderate to high, with for example, a heritability of 0.57 (0.14) for the percentage of fat. Estimates of slaughter traits varied between 0.08 (SLfat_om) and 0.54 (SLcw) with high standard errors (0.07–0.15), which is due to the relatively small number of animals. Previous studies indicate high heritabilities for slaughter characteristics. For instance, Kadarmideen et al. (2004) reported a heritability of 0.60 for premium meat cuts and 0.66 for intramuscular fat, while Singh et al. (2001) reported a heritability of 0.59 for carcass weight and 0.54 for backfat thickness.
We found most estimated heritabilities to be higher than previous reports, particularly for traits studied in field populations. One of the reasons for this is that our F2 pig resource population was created with the aim to generate large genetic variability in selected traits, which will consequently increase the heritability (as well as power to detect QTLs). In addition, our results cannot be directly compared with results based on field data and production pigs.
Genetic Predictions for Fatness and Fasting Glucose Level
The CVs show a wide range of phenotypic variation: 15–45% for traits directly related to OOR traits. Heritabilities (e.g., 0.57 for DXA%fat) showed that most of this overall variation is due to additive genetic effects, as expected from the F2 intercross. This was further investigated by plotting the distributions of estimated breeding values (EBV’s) for DXA%fat and FGL, as presented in Figure 4. It shows the wide and normally distributed genetic values within those OOR traits, resulting in some animals which are highly genetically predisposed for high or low body fatness and FGLs.
Figure 4. Distributions of estimated breeding values for percentage of fat (DXA fat) and fasting glucose level (GLU). Distributions of estimated breeding values (EBV’s) for (A) percentage of fat at DXA scanning and (B) fasting glucose level.
Phenotypic and Genetic Correlations
A multi-trait model that analyzes several traits simultaneously (Kadarmideen et al., 2004) in a one-step estimation procedure is desirable due to the high accuracy in prediction. However this was not feasible due to singularity problems and/or illogical estimates of covariance components so we fitted a series of two-trait models. A series of two-trait or bivariate animal models yielded estimates of phenotypic and genetic correlations and reported in Tables 4 and 5. Only statistically significant genetic and phenotypic correlations are reported and discussed here onward. However, we have presented phenotypic and genetic correlations (with standard errors) between all traits (regardless of statistical significance) for the complete dataset in Table S1 in Supplementary Material.
Table 4. Genetic (above the diagonal) and phenotypic (below the diagonal) correlations between DXA and weight and conformation traits with standard errors, estimated by individual animal model using all 563 animals.
Table 5. Genetic correlations of slaughter traits with weight and conformation with standard errors, estimated by individual animal model using all 563 animals.
Weight measurements show a high genetic correlation with conformation measurements (0.36–0.99) and DXA scanning results, as expected (Table 4). Previous research also indicated high correlations (0.51–0.75) between length and weight, for different production breeds (Johnson and Nugent, 2003).
In human research, BMI is popular due to simplicity of measurement and proven association with obesity. However, it is also known that BMI has many shortcomings (Heber, 2010). For example, ethnic differences and the influence of muscle mass (athletes) on weight result in misdiagnosis and misunderstanding of obesity. BAI is proposed as a better index for body adiposity (Bergman et al., 2011). BAI is expected to be a better indicator of per cent adiposity than measurements of body weight. However, critics raise the question whether this measure solves above-mentioned problems with BMI (Heymsfield and Shen, 2011). Heritability estimates for BMI and BAI are comparable, 0.56 and 0.58 respectively. To evaluate both measurements in order to determine which measurement would explain most about the measurements at slaughter, genetic correlations between BMI/BAI and slaughter traits were assessed. As expected, BAI in general gives a more negative genetic correlation with slaughter traits than BMI, e.g., 0.72 vs. −0.49 for meat percentage and −0.79 vs. 0.54 for backfat (mm). This stronger relationship indicates that BAI is a better indicator for adiposity than BMI. These results are novel, as to date no other studies have so far evaluated correlations between these traits.
Previous research has indicated a negative correlation between extreme obesity and BMD as well as between extreme obesity and BMC (Weiler et al., 2000). This correlation is mainly seen in postmenopausal women, resulting from decreased production of ovarian estrogens (Nunez et al., 2007), and is less strong in men than in women (Felson et al., 1993), resulting in a higher BMD for men (Arabi et al., 2004; Brown et al., 2004). In this F2 resource population, no significant difference between males and females was found (P = 0.09); however, females have a slightly higher BMC and BMD than males, which is not in line with human research. Conversely, a study in a Duroc × Pietrain crossbred population shows no gender effect for BMD and BMC (Laenoi et al., 2011). Strong positive genetic correlations were found between BMC and fat mass indicated by DXA scanning (0.85, SE = 0.06) and the same was evident for genetic correlations with backfat (mm) at slaughtering (0.89, SE = 0.15). Since obesity causes a lower BMD and BMC in humans, these results are not in line with the expected negative correlation. A possible explanation for these findings is the fact that measurements in our study are made on pigs at a very young age and thus, the results are not comparable to human studies.
It is recognized that adipose tissue is a highly active metabolic and endocrine organ, and dysfunction of adipose tissue is associated with insulin resistance (Goossens, 2008). This resource population has a wide variation in adiposity, and consequently it can be expected that there is also a wide variation in the level of metabolic parameters affected by obesity. As mentioned previously, pigs are a good model for human type 2 diabetes (Bellinger et al., 2006). In humans normal fasting plasma glucose levels are below 5.6 mmol/L, impaired plasma glucose levels are between 5.6 and 7 mmol/L, while type 2 diabetes has plasma glucose levels above 7 mmol/L. Previous research indicated lower fasting plasma glucose levels for the minipig (i.e., 3.6 mmol/L Larsen et al., 2003), whereas the domestic pig is estimated at the same level as humans (Larsen and Rolin, 2004). In our resource population, FGLs range between 1.80 and 11.15, with an average of 4.50 (SE = 2.18). The high heritability of FGL in this population (0.49) shows that there is a clear genetic background for FGLs in this population. There is a scarcity of research published on the genetics of FGLs in pigs, but our estimate indicates that glucose metabolism is heritable in pig populations. Studies of heritability for FGLs in humans show heritabilities between 0.38 and 0.51 for twin studies in Northern Europe, 0.08–0.72 for family studies and around 0.25 after adjustments in a genetically isolated population (Santos et al., 2006). Genetic parameter estimations show a positive correlation with weight (0.61), length (0.57), abdominal circumference (0.41), and fat mass at DXA scanning (0.59). The genetic correlations with the fat characteristics at slaughtering show a strong positive correlation (0.90), but a strong negative correlation with the meat characteristics (−0.56 to −0.76). As expected, these results show an elevated plasma glucose level with higher fat mass in the F2 population, which may indicate impaired glucose metabolism in some of the animals in this population (Goossens, 2008).
At approximately 100 kg, the pigs were slaughtered and dismembered. Weight characteristics show a positive genetic correlation with fat characteristics (between 0.40 and 0.76), and a negative genetic correlation with meat characteristics (between −0.45 and −0.76), which shows that an increase in weight of the pigs mainly results in an increase of fat (Table 5). The same was the case for conformation measurements: an increase in length mainly results in increase of backfat (0.56, SE = 0.25), and an increase of abdominal circumference results in increase of fat mass (0.66, SE = 0.23). Previous research indicated a positive correlation between length and backfat, between 0.10 and 0.41 for different production breeds (Johnson and Nugent, 2003), which is lower than estimates in this study (0.56, SE = 0.25). As mentioned earlier, the advantage of this F2 resource population is the possibility of measuring traits which are complicated or expensive to measure in human studies. One of these traits is the average daily gain. In this study strong genetic correlations are found between ADG and fat characteristics at slaughtering, e.g., between ADG and SLfat 0.74 (0.14). These results show that a strong growth in this F2 resource population mainly results in the accumulation of fat mass, as the genetic correlation with lean meat is negative (e.g., between ADGII and SL%meat −0.33, SE = 0.29). As expected on the basis of earlier mentioned correlations, a strong positive genetic correlation is found between fat characteristics at slaughtering and BMC (0.75, SE = 0.20) and the FGL (0.39, SE = 0.32).
In the current paper we describe a unique resource population consisting of 450 F2 pigs created by using divergent breeds (Yorkshire/Duroc × Göttingen Minipig) as founders. The population was designed with the purpose of establishing a pig model in which the genetic background of OOR and associated traits can be studied. The statistical analyses made on the phenotypic traits recorded in the population clearly demonstrate that potential QTL’s causing variation in the traits of interest are likely segregating in the F2 pig resource population. This study considered a total of 563 animals (454 with phenotypes and others are F0 and F1 parental animals) for estimation of genetic parameters. While this may be argued as low number of animals for parameter estimation, several other studies used fewer number of animals (fewer than 563 animals used in this study) to estimate genetic parameters, e.g., by Van Der Waaij et al. (2003). The established F2 population was measured for phenotypes with expected heritabilities in the range of around 0.25 to 0.50. From power calculations based on this assumption of heritabilities, we have concluded that around 450 F2 animals would be sufficient to detect even small QTLs explaining only around 2% of genetic variance at 80% power. Thus, genes involved in the phenotypes measured within the population can be identified in future genomic and systems genetic studies. This is also supported by previous QTL mapping studies using far fewer number of F2 animals (e.g., Ovilo et al., 2002; Hanotte et al., 2003). The whole resource population has been genotyped using the 60 K SNP Chip, which will be used in genome-wide association and fine QTL mapping studies as parts of the BioChild (www.biochild.ku.dk) and UNIK projects. Moreover, several tissues were collected at slaughtering, of which some will be used for RNA-sequencing to study transcriptomic profiles in extremely obese vs. lean pigs, as identified by the quantitative genetic analysis results reported here. In conclusion, the creation of this pig resource population resulted in a great potential to study the genetic architecture of OOR diseases.
Quantitative trait loci mapping or GWAS requires a thorough knowledge on genetic variances and any (genetic) relationships amongst traits that are investigated, particularly if such parameters for novel traits have not yet been reported elsewhere. This study has contributed to creating such knowledge for the established F2 pig resource population and provided a strong support for genomic investigations of OOR traits via GWAS and/or QTL mapping. Specifically, high heritabilities and strong correlations between weight, conformation, DXA, and slaughter traits between the different traits were reported here. This study, to the best of our knowledge, is the first experimental study in pigs in which a wide range of obesity traits, including novel DXA scan traits, information on different fatness traits at slaughter, and FGL have been measured in a large number of experimental animals. Consequently, it is the first study to report quantitative genetics characteristics on these novel OOR traits (in terms of estimated heritabilities and genetic correlations). In conclusion, the established F2 population gives a good rationale for further genome-wide and systems genetics approaches for detecting DNA variants, pathways, and gene networks affecting human obesity.
Merete Fredholm was the experimental project leader who originally conceived the F2 experiment and collected/supervised the measurements on all pigs. Haja N. Kadarmideen was the statistical genetics analysis project leader who originally conceived the quantitative genetics design, models, and analysis part and supervised Lisette J. A. Kogelman in the analyses. Lisette J. A. Kogelman collated, edited, and analyzed the data. Lisette J. A. Kogelman wrote the first draft of the manuscript with contributions from Haja N. Kadarmideen, Thomas Mark, and Merete Fredholm. Peter Karlskov-Mortensen, Camilla S. Bruun, Susanna Cirera, Thomas Mark, Mette J. Jacobsen, and Claus B. Jørgensen contributed to pedigree development and phenotypic measurements of the pigs. All authors read and approved the final version of the manuscript.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The authors would like to thank the Ministry of Science and Technology for the grant to the “UNIK Project for Food Fitness and Pharma for Health” and the University of Copenhagen for a PhD stipend awarded to Lisette Kogelman for financial support. They would also like to thank all the technical staff involved in breeding, managing, and slaughtering the animals as well as all the veterinary students helping with phenotypic measurements.
The Supplementary Material for this article can be found online at http://www.frontiersin.org/Livestock_Genomics/10.3389/fgene.2013.00029/abstract
Table S1. Phenotypic and genetic correlations of all analyzed variables. Excel file including all phenotypic and genetic correlations, with their standard errors.
Arabi, A., Nabulsi, M., Maalouf, J., Choucair, M., Khalife, H., Vieth, R., et al. (2004). Bone mineral density by age, gender, pubertal stages, and socioeconomic status in healthy Lebanese children and adolescents. Bone 35, 1169–1179.
Bahelka, I., Hanusová, E., Peškovicová, D., and Demo, P. (2007). The effect of sex and slaughter weight on intramuscular fat content and its relationship to carcass traits of pigs. Czech J. Anim. Sci. 52, 122–129.
Brown, L. B., Streeten, E. A., Shuldiner, A. R., Almasy, L. A., Peyser, P. A., and Mitchell, B. D. (2004). Assessment of sex-specific genetic and environmental effects on bone mineral density. Genet. Epidemiol. 27, 153–161.
Fan, B., Onteru, S. K., Du, Z.-Q., Garrick, D. J., Stalder, K. J., and Rothschild, M. F. (2011). Genome-wide association study identifies loci for body composition and structural soundness traits in pigs. PLoS ONE 6:e14726. doi:10.1371/journal.pone.0014726
Felson, D. T., Zhang, Y., Hannan, M. T., and Anderson, J. J. (1993). Effects of weight and body mass index on bone mineral density in men and women: the framingham study. J. Bone Miner. Res. 8, 567–573.
Fontanesi, L., Scotti, E., Buttazzoni, L., Davoli, R., and Russo, V. (2009). The porcine fat mass and obesity associated (FTO) gene is associated with fat deposition in Italian duroc pigs. Anim. Genet. 40, 90–93.
Groenen, M. A. M., Archibald, A. L., Uenishi, H., Tuggle, C. K., Takeuchi, Y., Rothschild, M. F., et al. (2012). Analyses of pig genomes provide insight into porcine demography and evolution. Nature 491, 393–398.
Hanotte, O., Ronin, Y., Agaba, M., Nilsson, P., Gelhaus, A., Horstmann, R., et al. (2003). Mapping of quantitative trait loci controlling trypanotolerance in a cross of tolerant West African N’Dama and susceptible East African Boran cattle. Proc. Natl. Acad. Sci. U.S.A. 100, 7443–7448.
Hodgson, R. R., Davis, G. W., Smith, G. C., Savell, J. W., and Cross, H. R. (1991). Relationships between pork loin palatability traits and physical characteristics of cooked chops. J. Anim. Sci. 69, 4858–4865.
Hofer, A., Hagger, C., and Künzi, N. (1992). Genetic evaluation of on-farm tested pigs using an animal model I. Estimation of variance components with restricted maximum likelihood. Livest. Prod. Sci. 30, 69–82.
Johansen, T., Hansen, H. S., Richelsen, B., and Malmlöf, K. (2001). The obese gottingen minipig as a model of the metabolic syndrome: dietary effects on obesity, insulin sensitivity, and growth hormone profile. Comp. Med. 51, 150–155.
Johnson, Z. B., and Nugent, R. A. (2003). Heritability of body length and measures of body density and their relationship to backfat thickness and loin muscle area in swine. J. Anim. Sci. 81, 1943–1949.
Kadarmideen, H. N., Janss, L. L. G., and Dekkers, J. C. M. (2000). Power of quantitative trait locus mapping for polygenic binary traits using generalized and regression interval mapping in multi-family half-sib designs. Genet. Res. 76, 305–317.
Kadarmideen, H. N., Schwörer, D., Ilahi, H., Malek, M., and Hofer, M. (2004). Genetics of osteochondral disease and its relationship with meat quality and quantity, growth, and feed conversion traits in pigs. J. Anim. Sci. 82, 3118–3127.
Kim, K. S., Larsen, N., Short, T., Plastow, G., and Rothschild, M. F. (2000). A missense variant of the porcine melanocortin-4 receptor (MC4R) gene is associated with fatness, growth, and feed intake traits. Mamm. Genome 11, 131–135.
Kogelman, L. J. A., Byrne, K., Vuocolo, T., Watson-Haigh, N., Kadarmideen, H. N., Kijas, J., et al. (2011). Genetic architecture of gene expression in ovine skeletal muscle. BMC Genomics 12:607. doi:10.1186/1471-2164-12-607
Laenoi, W., Uddin, M., Cinar, M., Grosze-Brinkhaus, C., Tesfaye, D., Jonas, E., et al. (2011). Quantitative trait loci analysis for leg weakness-related traits in a Duroc × Pietrain crossbred population. Genet. Sel. Evol. 43, 13.
Larsen, M., Gotfredsen, C., Wilken, M., Carr, R., Pørksen, N., and Rolin, B. (2003). Loss of beta-cell mass leads to a reduction of pulse mass with normal periodicity, regularity and entrainment of pulsatile insulin secretion in Göttingen minipigs. Diabetologia 46, 195–202.
Liu, G., Jennen, D. G. J., Tholen, E., Juengst, H., Kleinwächter, T., Hölker, M., et al. (2007). A genome scan reveals QTL for growth, fatness, leanness and meat quality in a Duroc-Pietrain resource population. Anim. Genet. 38, 241–252.
Nunez, N. P., Carpenter, C. L., Perkins, S. N., Berrigan, D., Jaque, S. V., Ingles, S. A., et al. (2007). Extreme obesity reduces bone mineral density: complementary evidence from mice and women. Obesity 15, 1980–1987.
Okumura, N., Matsumoto, T., Hayashi, T., Hirose, K., Fukawa, K., Itou, T., et al. (2012). Genomic regions affecting backfat thickness and cannon bone circumference identified by genome-wide association study in a Duroc pig population. Anim. Genet. 1365–2052.
Ovilo, C., Clop, A., Noguera, J. L., Oliver, M. A., Barragán, C., Rodríguez, C., et al. (2002). Quantitative trait locus mapping for meat quality traits in an Iberian × Landrace F2 pig population. J. Anim. Sci. 80, 2801–2808.
Ponsuksili, S., Murani, E., Brand, B., Schwerin, M., and Wimmers, K. (2011). Integrating expression profiling and whole-genome association for dissection of fat traits in a porcine model. J. Lipid Res. 52, 668–678.
Rubin, C.-J., Megens, H.-J., Barrio, A. M., Maqbool, K., Sayyab, S., Schwochow, D., et al. (2012). Strong signatures of selection in the domestic pig genome. Proc. Natl. Acad. Sci. U.S.A. 109, 19529–19536.
Santos, R. L. P., Zillikens, M. C., Rivadeneira, F. R., Pols, H. A. P., Oostra, B. A., Duijn, C. M., et al. (2006). Heritability of fasting glucose levels in a young genetically isolated population. Diabetologia 49, 667–672.
Scott, L. J., Mohlke, K. L., Bonnycastle, L. L., Willer, C. J., Li, Y., Duren, W. L., et al. (2007). A genome-wide association study of type 2 diabetes in finns detects multiple susceptibility variants. Science 316, 1341–1345.
Speliotes, E. K., Willer, C. J., Berndt, S. I., Monda, K. L., Thorleifsson, G., Jackson, A. U., et al. (2010). Association analyses of 249,796 individuals reveal 18 new loci associated with body mass index. Nat. Genet. 42, 937–948.
Suster, D., Leury, B. J., Hofmeyr, C. D., D’souza, D. N., and Dunshea, F. R. (2004). The accuracy of dual energy X-ray absorptiometry (DXA), weight, and P2 back fat to predict half-carcass and primal-cut composiion in pigs within and across research experiments. Aust. J. Agric. Res. 55, 973–982.
Teslovich, T. M., Musunuru, K., Smith, A. V., Edmondson, A. C., Stylianou, I. M., Koseki, M., et al. (2010). Biological, clinical and population relevance of 95 loci for blood lipids. Nature 466, 707–713.
Uddin, M., Duy, D., Cinar, M., Tesfaye, D., Tholen, E., Juengst, H., et al. (2011). Detection of quantitative trait loci affecting serum cholesterol, LDL, HDL, and triglyceride in pigs. BMC Genet. 12:62. doi:10.1186/1471-2156-12-62
Van Der Waaij, E. H., Hanotte, O., Van Arendonk, J. A. M., Kemp, S. J., Kennedy, D., Gibson, J. P., et al. (2003). Population parameters for traits defining trypanotolerance in an F2 cross of N’Dama and Boran cattle. Livest. Prod. Sci. 84, 219–230.
World Health Organization, (2012). Obesity and overweight, Fact sheet N°311, updated March 2013. Available at: http://www.who.int/mediacentre/factsheets/fs311/en/
Keywords: animal model, F2 design, obesity, diabetes, heritabilities, genetic correlations, genetic predictions
Citation: Kogelman LJA, Kadarmideen HN, Mark T, Karlskov-Mortensen P, Bruun CS, Cirera S, Jacobsen MJ, Jørgensen CB and Fredholm M (2013) An F2 pig resource population as a model for genetic studies of obesity and obesity-related diseases in humans: design and genetic parameters. Front. Genet. 4:29. doi: 10.3389/fgene.2013.00029
Received: 07 December 2012; Accepted: 22 February 2013;
Published online: 18 March 2013.
Edited by:Tad S. Sonstegard, United States Department of Agriculture, USA
Reviewed by:Luca Fontanesi, University of Bologna, Italy
Jeffrey O’Connell, University of Maryland School of Medicine, USA
Copyright: © 2013 Kogelman, Kadarmideen, Mark, Karlskov-Mortensen, Bruun, Cirera, Jacobsen, Jørgensen and Fredholm. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in other forums, provided the original authors and source are credited and subject to any copyright notices concerning any third-party graphics etc.
*Correspondence: Merete Fredholm, Faculty of Health and Medical Sciences, University of Copenhagen, Grønnegårdsvej 3, 1870 Frederiksberg C, Copenhagen, Denmark. e-mail: firstname.lastname@example.org