About the Authors:
Simona D. Frederiksen
Contributed equally to this work with: Simona D. Frederiksen, Peter Karlskov-Mortensen
Affiliation: Department of Veterinary and Animal Sciences, Faculty of Health and Medical Sciences, University of Copenhagen, Copenhagen, Denmark
Peter Karlskov-Mortensen
Contributed equally to this work with: Simona D. Frederiksen, Peter Karlskov-Mortensen
Affiliation: Department of Veterinary and Animal Sciences, Faculty of Health and Medical Sciences, University of Copenhagen, Copenhagen, Denmark
Sameer D. Pant
Affiliation: School of Animal and Veterinary Sciences, Charles Sturt University, Wagga Wagga, Australia
Maryse Guerin
Affiliation: INSERM UMR_S1166, Integrative Biology of Atherosclerosis Team, Paris, France
Philippe Lesnik
Affiliation: INSERM UMR_S1166, Integrative Biology of Atherosclerosis Team, Paris, France
Claus B. Jørgensen
Affiliation: Department of Veterinary and Animal Sciences, Faculty of Health and Medical Sciences, University of Copenhagen, Copenhagen, Denmark
Susanna Cirera
Affiliation: Department of Veterinary and Animal Sciences, Faculty of Health and Medical Sciences, University of Copenhagen, Copenhagen, Denmark
Camilla S. Bruun
Affiliation: Department of Veterinary and Animal Sciences, Faculty of Health and Medical Sciences, University of Copenhagen, Copenhagen, Denmark
Thomas Mark
Affiliation: Novo Nordisk, Scandinavia AB, Region Denmark, Maaloev, Denmark
Merete Fredholm
* E-mail: [email protected]
Affiliation: Department of Veterinary and Animal Sciences, Faculty of Health and Medical Sciences, University of Copenhagen, Copenhagen, Denmark
ORCID http://orcid.org/0000-0002-3563-7648
Introduction
Obesity is defined as excessive accumulation of fat in the body to the extent that it may have a negative effect on health. According to the World Health Organization (WHO), in 2016 obesity had more than doubled worldwide since 1980. Obesity can be socially stigmatizing, however, in itself it is not the primary health problem. Rather, it is the co-morbidities such as cardiovascular disease (CVD) and type 2 diabetes (T2D) which pose the major health problems. Common to many of the co-morbidities is an unhealthy metabolic profile with insulin resistance and dyslipidemia characterized by elevated triglycerides (TG) and low-density lipoprotein cholesterol (LDL-C) levels and by decreased high-density lipoprotein cholesterol (HDL-C) concentrations [1]. Approximately 25% of obese humans do, however, not present the metabolic complications and do not suffer from an increased susceptibility towards obesity related diseases [2, 3]. These individuals are characterized as 'metabolically healthy obese' (MHO) [4].
To understand why MHO individuals have a better prognosis than the metabolically unhealthy individuals it is relevant to elucidate the genetic mechanisms underlying regulation of blood lipids such as cholesterol and TG.
Cholesterol is a ubiquitous steroid and a vital component of cellular membranes in vertebrates. But at the same time, subendothelial accumulation of cholesterol is the cause of atherosclerotic lesions leading to vascular diseases, heart attacks, aortic aneurysms and peripheral vascular diseases [5] which, altogether, represent the most frequent causes of demise in the industrialized world [6].
Cholesterol is transported in plasma in lipoprotein particles. These consist of a major structural apolipoprotein, peripheral apolipoproteins, structural lipids (phospholipids and cholesterol) and a cargo of TG and steryl esters. Dietary lipids are absorbed via the small intestine and packaged as TG into large particles known as chylomicrons (CM) which distribute lipids throughout the body for direct use and storage. CM remnants return to the liver where the remaining lipid cargo enters hepatic lipid pathways and a second class of lipoprotein particles is produced; the very low density lipoproteins (VLDL). Like CM, VLDL distributes TG to the periphery. As VLDL delivers TG they shrink and shed peripheral apolipoproteins, hence the proportion of cholesterol increases and VLDL becomes low-density lipoproteins (LDL). Both VLDL and LDL are characterized by the major structural apolipoprotein APOB. LDL can bind to cells expressing the LDL-receptor (LDLR) and as such, LDL acts as an efficient cholesterol delivery system into cells [5].
The reverse transport of cholesterol from the periphery to the liver is mediated by high-density lipoproteins (HDL) containing the major structural apolipoprotein APOA1. Nascent HDL particles contain few lipids but collect free cholesterol and phospholipids from peripheral tissues. Different mechanisms enrich the HDL particles with lipids, which ultimately are delivered to the liver via Scavenger Receptor B1 (SCARB1) for subsequent excretion with the bile [5].
Heritability for LDL-C and HDL-C levels in humans are estimated to be around 70%.[7] Large-scale genome wide association studies (GWAS) in humans have discovered over 150 common genetic variants associated with plasma lipids [8]. However, these loci only explain a small fraction of the total variance in blood lipids [9] and the bulk of genetic factors for dyslipidemia are still unaccounted for.
A way to elucidate novel mechanisms involved in blood lipid regulation is to study animal models. Animal models provide the benefits of a strictly controlled diet and environment which is impossible in large human studies. Consequently, random noise is greatly reduced and a corresponding increase in power is presumed when animal models are used to study human conditions. Also, the individual breeds are genetically more homogeneous than humans due to domestication and artificial selection [e.g. ref. 10]. GWAS [11] and quantitative trait locus (QTL) mapping [12] have been performed in mice to identify loci in the genome with effect on blood lipids. But more frequently, spontaneous dyslipidemic and genetically engineered mouse models have been used to study the effect of specific dyslipidemia associated genes first identified in human GWAS studies [13]. However, market differences in metabolism and adipose tissue biology exists between rodents and humans [e.g. ref. 14].
The pig is an animal model with a close similarity to humans in body size, physiology, organ development and disease progression. It has been widely used as a model for cardiovascular and metabolic diseases [15]. The presence of atherosclerotic lesions in aorta was described in pigs as early as 1954 [16]. The close similarity to vascular lesions in human atherosclerosis have later been confirmed [17, 18], and the pig has proved its value as a model for this disease [19].
In the present study, a GWAS aimed at identifying loci with effect on blood lipid levels is reported. The GWAS was performed in two pig crosses established by using Göttingen Minipig as the parental boar line in both crosses whereas Duroc and Yorkshire production pigs were used as parental sow lines in the two crosses respectively. The Göttingen Minipig is an obesity prone pig breed [20] often used in studies of obesity, diabetes and metabolic syndrome [21–24], whereas Duroc and Yorkshire pigs have been bred for leanness for decades. Haplotypes were defined around the most significantly associated single nucleotide polymorphisms (SNPs), and the effect of the haplotypes with the highest additive effects were studied further.
Materials and methods
Experimental animal model, sample collection and blood lipid analysis
Two populations were produced as F2 crosses using purebred Göttingen Minipig (M), Duroc (D) and Yorkshire (Y) as parental lines (Ellegaard Göttingen Minipigs A/S, Dalmose, Denmark; DanBred International, Herlev, Denmark). The Minipig-Duroc (MD) crossbred F1 animals were founded by seven M boars and seven D sows and 28 F1 MD gilts and 16 F1 MD boars were used to produce 285 F2 animals. Similarly, the Minipig-Yorkshire (MY) crossbred F1 animals were founded by seven M boars and seven Y sows and 279 F2 animals were produced by 26 MY gilts and 13 MY boars (for further information see [25, 26]).
All pigs were raised under controlled environmental conditions and fed the same diet ad libitum. The project was approved by the Danish Animal Experiments Inspectorate. Animal care and maintenance were conducted according to the Danish “Act about Animal Husbandry” and “Animal Protection Act” (Act 432, July 9, 2004; Act 1150, Sep. 12, 2015). All pigs were housed at a regular pig farm, and slaughtered at a commercial slaughterhouse by stunning and bleeding under veterinary supervision.
Blood samples for blood lipid analysis were drawn from the jugular vein at about two month of age (63 ± 10 days, abbreviated Age1) and at slaughter (242 ± 48 days, abbreviated Age2). Plasma lipid levels were assayed by standardized techniques using a Konelab 20 Clinical Chemistry Analyzer (Thermo Scientific, Sweden) and commercial reagent kits from Roche Diagnostics for Total Cholesterol (CT) and from ThermoElectron for TG and HDL-C levels (direct method). LDL-C levels were calculated using the Friedewald formula. Observations >5 SD from the mean were considered outliers and excluded. Box-Cox transformation was used to adjust for non-normality. After transformation, skewness and kurtosis were calculated for each phenotypic distribution and Q-Q plots were made to evaluate normality.
Estimation of genetic parameters
Genetic parameters were estimated using Best Linear Unbiased Prediction (BLUP) based on Average-Restricted Maximum Likelihood (AI-REML) using DMU version 6, release 5.2 [27]. Variance components were estimated using two different univariate models depending on whether the phenotypes were measured at Age1 (Model 1) or Age2 (Model 2):where y is the phenotype for animal i, μ is the population mean, SEX and CROSS are fixed effects for sex and cross (MD or MY) for animal i, AGE for animal i is a covariate with regression coefficient γ and A is the additive genetic effect for animal i with the regression coefficient ξ. εi is the residual error. A and ε are assumed to be independent and normally distributed with variances AσANIMAL2 and Iσε2, respectively. A is the additive genetic relationship matrix based on pedigree information, σANIMAL2 the additive genetic variance, I an identity matrix of appropriate size and σε2 the residual error variance.
A bivariate linear mixed model was fitted to estimate phenotypic and genetic correlations between pairs of traits (Trait 1 and Trait 2):where y1 and y2 are vectors of the transformed phenotypic measurements of Trait 1 and Trait 2 respectively. X1 and X2 are design matrices connecting the measurements to b1 and b2 which are vectors of the environmental fixed effects, while Z1 and Z2 are design matrices linking the measurements to a1 and a2 which are vectors for the random genetic effects for each trait. Finally, e1 and e2 are vectors of the random environmental effects. The effects fitted in the bivariate linear mixed models are the same as the ones fitted in the univariate linear mixed models. When fitting the bivariate model, the non-diagonals of the additive genetic variance and covariance matrix G for animal effects and the variance and covariance matrix R for random residual effects were included as follows:where covA12 = covA21 is the additive genetic covariance and covE12 = covE21 the random residual covariance between Trait 1 and 2. The phenotypic covariance covP12 = covP21 between Trait 1 and 2 can be calculated as: covP12 = covA12 + covE12 [28].
Genotyping and Quality Control
Each animal was genotyped using the 60k porcine Illumina SNPchip. Sex chromosomes were excluded from the single-marker analyses. The Sus Scrofa 10.2 pig genome assembly was used to derive map positions for all SNPs. After genotyping, quality control (QC) was conducted using GenABEL version 1.8–0 [29] excluding: animals and SNPs with more than 5% missing genotypes; SNPs with a minor allele frequency (MAF) less than 5%; SNPs that significantly deviated from Hardy-Weinberg equilibrium (HWE) at a significance threshold of p < 1E-05. After QC, 549 pigs and 44,554 SNPs remained.
GWAS using single-marker test and estimation of SNP effects
GWAS was carried out for each phenotype using a single-marker test as implemented in GenABEL [30]. The analyses were performed in two steps: In step 1, the polygenic linear mixed models were defined as Model 1 (for Age1) and Model 2 (for Age2) where Ai is the random additive polygenic effect for animal i based on identity by state (IBS). Ai is assumed to be normally distributed with a (co)variance of GσANIMAL2 where G is the genomic relationship matrix and σANIMAL2 the additive polygenic variance. Estimation of effects was performed by means of maximum likelihood (ML) [31].
Step 2 use estimated residuals from step 1 (contains part of the QTL variance) to estimate SNP effects [30]:where is the estimated residual error from step 1, μ is the mean, SNPm is assumed a fixed effect with 3 genotype scores 0, 1 and 2 referring to the number of minor allele copies for the m’th SNP and εm is the residual error. Each single SNP was modeled independently. The allele substitution effect of the m’th SNP was calculated as the average phenotypic change when replacing a major allele with the minor allele [32].
SimpleM [33–35] was used to correct for multiple testing by calculating the effective number of independent tests Meff = 12916. Missing genotypic values were replaced with the common allele genotype in the genotype matrix prior to the analysis. It resulted in a cut off value of–log10(0.05/12916) = 5.41 corresponding to a nominal significance level (α) of 3.87*10−6.
For each significant SNP, the phenotypic variance explained by the SNP was calculated as [36]:where Vm is the phenotypic variance explained by the m’th SNP, qm is MAF for the m’th SNP and is the estimated effect for the m’th SNP. The variance explained by the SNP is assumed to be due to additive genetic variance. The proportion of additive genetic variance explained by each SNP was calculated as VSNP / σA2. The additive genetic effect (α) for each top-SNP was calculated using BLUP as implemented in DMU where the effect of having 0, 1 or 2 copies of a specific SNP was fitted in the models described above (Model 1 and Model 2).
Linkage disequilibrium and haplotype analysis
The method described by Gabriel et al.(2002) [37] as implemented in Haploview [38] was used to define haploblocks around the top SNPs identified in the GWAS. Default parameters were used except for “Fraction of strong LD in informative comparisons” which was set to >0.80 (default >0.95). Phased haplotypes for each haploblock for each individual were estimated using the—hap-phase option in Plink [39]. The additive genetic effect of each haplotype was estimated using BLUP by fitting the haplotype as a covariate in the models described above (Model 1 and Model 2). Effect were fitted for one haplotype variant at a time by computing the effect of having 0, 1 or 2 copies of that variant. Effects were only calculated for haplotype variants observed in the founder animals and not for variants found in less than 20 animals in the F2 population.
To evaluate haplotype effects in individuals with high- and low-BMI respectively, all F2 animals were sorted according to BMI and the one-third with highest BMI (high BMI) and lowest (low BMI) respectively were selected. Each of the groups were subdivided according to haplotypes in the HDL-C/CT and the LDL-C QTL regions and mean and standard deviation for BMI_Age2, CT_Age2, TG_Age2, HDL-C/CT_Age1, HDL-C_Age2, HDL-C/CT_Age2, LDL-C_Age2 and LDL-C/HDL-C at Age2 were calculated separately using Students-t test for animals which were homozygous for either the GM or the Yorkshire/Duroc haplotypes in the HDL-C/CT and LDL-C QTL regions.
Results
A moderate genetic correlation was found between Age1 and Age2 both in regard to LDL-C and HDL-C/CT, whereas, for TG there was no correlation between measurements at the two ages (see Table 1).
[Figure omitted. See PDF.]
Table 1. Genetic (rA) and phenotypic (rP) correlation between phenotypes measured at different age.
https://doi.org/10.1371/journal.pone.0178828.t001
The GWA study resulted in the identification of several QTL regions below the nominal significance level on sus scrofa chromosome 3 (SSC3) (see Table 2). A single SNP on SSC3 at 124.7 Mb is associated with TG_Age2, whereas, there is no indication for an association with TG_Age1 at this position (data not shown). Conversely, the same loci on SSC3 are associated with LDL-C at Age 1 and Age 2 and the same is true for HDL-C/CT. The most significant associations for LDL-C are found around 125.6 Mb and the most significant SNPs for HDL-C/CT are found at 122.8 Mb.
[Figure omitted. See PDF.]
Table 2. QTL regions and SNP effects.
https://doi.org/10.1371/journal.pone.0178828.t002
LD analysis (Fig 1) revealed that the two top SNPs for HDL-C/CT_Age1 and HDL-C/CT_Age2 are located close together and belong to the same LD-block defined by six SNPs from 122.7–122.9 Mb (designated '1' in Fig 1). The top SNP’s for TG_Age2 and LDL-C_Age1 and Age2 show minimal LD to surrounding SNPs and are not included in any LD blocks when default parameters in Haploview are used (results not shown). Relaxing the parameters for block definition, as described above, includes the top SNP for LDL-C_Age2 into an LD block defined by 10 SNPs spanning a region from 125.6–125.9 Mb (designated '3' Fig 1), but the top SNP for LDL-C_Age1 is still outside the LD block. In the same way the relaxed parameters include the top SNP for TG_Age2 into a four SNP haploblock from124.67–124.75 Mb (designated '2' in Fig 1) clearly separated from the HDL-C/CT and the LDL-C loci. The SNPs located in the haploblocks associated with the lipid levels are listed in SI 1.
[Figure omitted. See PDF.]
Fig 1. LD in a 2.8 Mb region (122.4–126.1 Mb) on Sus Scrofa chromosome 3.
Frames mark top SNPs for HDL-C/CT_Age1, HDL-C/CT_Age2, TG_Age2, LDL-C_Age1 and LDL-C_Age2, in this order from left to right. Triangles indicate haploblock structure in the QTL region. Block 1: QTL for HDL-C/CT; Block 2: QTL for TG; Block 3: QTL for LDL-C.
https://doi.org/10.1371/journal.pone.0178828.g001
The additive genetic effect of haplotypes present in more than 20 animals is shown in Table 3. The largest additive effects appear for the HDL-C/CT and LDL-C haplotypes for which there is also a clear breed specific segregation. For the HDL-C/CT associated haplotypes one specific haplotype, AAAGGG, impose a higher HDL-C/CT level. Twelve out of fourteen Göttingen Minipig founders are homozygous for this haplotype and the remaining two has one copy of AAAGGG. The haplotype GGGAAA, for which all Yorkshire founders are homozygous, has an equivalent negative effect, that is, animals with this haplotype have lower HDL-C/CT levels. The same trend is true for the haplotypes segregating from the Duroc breed, and among the Duroc haplotypes, AGAAGA, has an even stronger negative effect compared to the Yorkshire haplotype on the HDL-C/CT level. The effect of the individual haplotypes is similar at Age 1 and Age 2 (data not shown).
[Figure omitted. See PDF.]
Table 3. Haplotype effects—SSC3.
https://doi.org/10.1371/journal.pone.0178828.t003
For the LDL-C associated LD block, the haplotype AGGAGCCAAA originating from the Göttingen Minipig founders, decrease the LDL-C level. The Yorkshire founder specific haplotypes, GGAAGAAGAG and GGAAGAAGCG, increase LDL-C level; however, there is also one Yorkshire haplotype, GGAAGAAACG, which decreases LDL-C level. Two of the Duroc haplotypes increase the LDL-C level slightly. The effect of each haplotype is the same at Age 1 and Age 2 (data not shown) even though the blocks are only loosely in LD according to the LD analysis.
To further evaluate the phenotypic consequences of the haplotypes segregating from Göttingen Minipig and Yorkshire/Duroc respectively we compared lipid and TG levels in animals homozygous for either Göttingen Minipig haplotypes or the Duroc/Yorkshire haplotypes. The comparison was performed in the one-third of the F2 animals with the highest and one-third with lowest BMI, respectively. The mean BMI in the combined F2 population (564 animals) is 127±20 while in the selected high BMI group (136 animals) the mean BMI is 154±13, and in the selected low BMI group (132 animals) the mean BMI is 111±9. As seen in Table 4, within the low BMI group of animals there are no differences in TG and lipid values between the genetically diverse animals. However, when comparing animals homozygous for Göttingen Minipig haplotypes with animals homozygous for Yorksire/Duroc haplotypes within the high BMI group of animals both TG and lipid values (except HDL-C) are significantly higher in animals with the Yorkshire/Duroc haplotypes.
[Figure omitted. See PDF.]
Table 4. Effects of HDL-C/CT and LDL-C haplotypes with high additive effect in high- and low-BMI individuals.
https://doi.org/10.1371/journal.pone.0178828.t004
Discussion
In agreement with previously performed GWA studies in pigs [26, 40, 41], we have identified a QTL region influencing lipoprotein traits on SSC3. The homologous region on human chromosome 2 has also been implicated in the regulation of blood lipids in human GWA studies [42]. Our study shows that three separate haplopblocks influence HDL-C/CT -, TG -, and LDL-C levels respectively (see Fig 1). The HDL-C/CT and HDL-C QTL loci flank the TG_Age2 locus and thus, this 2.8 Mb region on SSC3 seems to encompass several loci with a regulatory effect on different plasma lipids. This is underpinned by the observation that the SSC3 locus affects HDL-C and LDL-C levels independent of age but for TG-levels there is only an association to TG_Age2 and no evidence suggesting an association to TG_Age1. This indicates that different genes or regulatory mechanisms in the region are responsible for effects on cholesterol and triglyceride levels, respectively. This is also in agreement with the observed correlations between phenotypes at different ages (Table 1) and the lack of a correlation between TG at Age1 and Age2.
There are no annotated known genes in the HDL-C/CT and TG associated SSC3 haploblocks in the Sscrofa 10.2 assembly but the gene encoding ras homolog family member B (RHOB) is located in the LDL-C associated haploblock. RHOB is a member of the RHO GTP-binding protein family which regulates expression of CD36 [43]. CD36 is a scavenger protein with high affinity for plasma lipoproteins including LDL and a high expression in tissues such as skeletal muscle, heart, mammary, epithelium and adipose tissue, with a very active fatty acid metabolism [44]. The gene encoding apolipoprotein B (APOB) is located at 125.23–125.35 Mb. It is an obvious candidate gene for lipid related traits but in the present study it is not included in associated LD-blocks even though it was covered by four SNPs with a MAF above the exclusion threshold. The gene encoding lipid droplet associated hydrolase (LDAH) is located at 125.42–125.53 Mb which is also outside the identified LD-blocks. Further studies including sequencing and genotyping of additional markers are warranted to ascertain LD structure and to identify genetic variation. It is very likely that the observed differences in blood lipids are caused by variation in regulatory components within the identified QTL regions in particular since a high number of different genes involved in regulation of lipoproteins are located within this region of the genome.
Investigation of the additive genetic effect of the most common haplotypes segregating within the F2 population indicates that the haplotypes originating from the Göttingen Minipig breed provide a healthier lipid profile compared to the haplotypes segregating from the Yorkshire and Duroc breeds. I.e., Göttingen Minipig haplotypes in haploblock 1 appears to increase HDL-C/CT compared to the Yorkshire and Duroc haplotypes, which have the opposite effect. Conversely, in haploblock 3 the Göttingen Minipig haplotype lowers the level of LDL-C and the Yorkshire and Duroc haplotype increases LDL-C level with the exception of the GGAAGAAACG haplotype segregating from the Yorkshire breed. This effect is confirmed by comparing animals homozygous for the Göttingen Minipig and Duroc/Yorkshire haplotypes respectively within groups of animals with high and low BMI respectively. Within the high BMI group animals with Göttingen Minipig haplotypes have a significantly lower TG level, lower levels of LDL-C, and LDL-C/HDL-C ratio and higher ratio of HDL-C to CT (Table 4). Within the low BMI group there is no difference in the phenotypes between the two genetic variants. Thus, overall Göttingen Minipig seems to have a genotype that supports a more healthy blood lipid profile in spite of the fact that they are prone to obesity. Or conversely, pigs without the Göttingen Minipig haplotypes develop a more unhealthy, dyslipidemic profile together with obesity, compared to pigs with Göttingen Minipig haplotypes which uphold a healthy lipid profile despite development of severe obesity.
The Gottingen Minipig breed was developed in the 1960's using Minnesota Minipigs, Vietnamese Potbelly Pigs, and German Landrace as founders (Simianer and Köhn, 2010). A likely explanation for the MHO profile in Göttingen Minipigs is that the obesity prone minipig founders have been adapted to overcome obesity by natural selection. Thus, unexpectedly, the results presented here indicate that the Göttingen Minipig breed is not well suited for studies of the obesity related co-morbidities but may be a valuable model to advance understanding of the MHO phenotype in humans.
The MHO profile identified in our F2 population is comparable to the MHO profile in humans which appear to be protected against obesity related metabolic complications. However, although MHO is an important, emerging phenotype in humans no universally accepted definition has been established for this phenotype yet [45]. It is also debated to what extend MHO individuals will remain healthy [46, 47]. On the other hand, studying MHO subjects may lead to better intervention strategies for metabolically unhealthy obese people, and elucidate if they by lifestyle changes or by the use of medicine can switch to a better metabolic profile.
In conclusion, we have substantiated that different genetic loci have an effect on TG early and late in life. We have also substantiated that the genomic region close to genes implicated in lipoprotein metabolism (RHOB, APOB, LDAH) comprise regulatory elements of importance to the regulation of lipid metabolism. Interestingly, the largest additive genetic effects of the haplotypes identified in this study show that haplotypes segregating from the obesity prone Göttingen Minipig breed is able to uphold a healthy lipid profile despite development of obesity. Thus, the genetic profile of the Göttinging Minipig breed seems to support a phenotype comparable to the MHO phenotype in humans promoting this pig breed as a model for further studies of this particular phenotype. Our future studies will be directed at identification the genetic variation in the regulatory components involved in lipid metabolism and further genetic characterization of the healthy metabolic phenotype.
Acknowledgments
We thank Tina Neergaard Mahler, Minna Baron Jakobsen, and Charlotte Bjørner Larsen for excellent technical assistance. We acknowledge the financial support by The Danish Independent Research Council Grant Number DFF 1335–00127.
Author Contributions
1. Conceptualization: MF PKM CBJ SC CSB.
2. Data curation: SDF PKM SDP.
3. Formal analysis: SDF PKM SDP.
4. Funding acquisition: MF.
5. Investigation: SDF PKM SDP TM MF.
6. Methodology: SDF PKM SDP MG PL.
7. Project administration: MF.
8. Resources: MF.
9. Supervision: MF SDP TM.
10. Validation: MF TM.
11. Visualization: PKM.
12. Writing – original draft: SDF PKM MF.
13. Writing – review & editing: SDF PKM SDP MG PL CBJ SC CSB TM MF.
Citation: Frederiksen SD, Karlskov-Mortensen P, Pant SD, Guerin M, Lesnik P, Jørgensen CB, et al. (2017) Haplotypes on pig chromosome 3 distinguish metabolically healthy from unhealthy obese individuals. PLoS ONE 12(6): e0178828. https://doi.org/10.1371/journal.pone.0178828
1. Grundy SM. Dyslipidaemia in 2015: Advances in treatment of dyslipidaemia. Nature reviews Cardiology. 2016;13(2):74–5. Epub 2016/01/15. pmid:26763537
2. Navarro E, Funtikova AN, Fito M, Schroder H. Can metabolically healthy obesity be explained by diet, genetics, and inflammation? Molecular nutrition & food research. 2015;59(1):75–93. Epub 2014/11/25.
3. Bosello O, Donataccio MP, Cuzzolaro M. Obesity or obesities? Controversies on the association between body mass index and premature mortality. Eating and weight disorders: EWD. 2016. Epub 2016/04/05.
4. Blüher M. The distinction of metabolically ‘healthy’ from ‘unhealthy’ obese individuals. Current Opinion in Lipidology. 2010;21(1):38–43. pmid:19915462
5. Vance JE, Vance DE. Biochemistry of Lipids, Lipoproteins and Membranes. 5 ed. Amsterdam: Elsevier Science; 2008.
6. Go AS, Mozaffarian D, Roger VL, Benjamin EJ, Berry JD, Blaha MJ, et al. Heart Disease and Stroke Statistics—2014 Update: A report from the American Heart Association. Circulation. 2014;129(3):e28–e292. pmid:24352519
7. Weiss LA, Pan L, Abney M, Ober C. The sex-specific genetic architecture of quantitative traits in humans. Nat Genet. 2006;38(2):218–22. pmid:16429159
8. Global Lipids Genetics Consortium. Discovery and refinement of loci associated with lipid levels. Nat Genet. 2013;45(11):1274–83. pmid:24097068
9. Teslovich TM, Musunuru K, Smith AV, Edmondson AC, Stylianou IM, Koseki M, et al. Biological, clinical and population relevance of 95 loci for blood lipids. Nature. 2010;466(7307):707–13. pmid:20686565
10. Groenen MA. A decade of pig genome sequencing: a window on pig domestication and evolution. Genetics, selection, evolution: GSE. 2016;48:23. Epub 2016/03/31. pmid:27025270
11. Pletcher MT, McClurg P, Batalov S, Su AI, Barnes SW, Lagler E, et al. Use of a dense single nucleotide polymorphism map for in silico mapping in the mouse. PLoS biology. 2004;2(12):e393. Epub 2004/11/10. pmid:15534693
12. Wang X, Paigen B. Genetics of Variation in HDL Cholesterol in Humans and Mice. Circulation Research. 2005;96(1):27–42. pmid:15637305
13. Kuivenhoven JA, Hegele RA. Mining the genome for lipid genes. Biochimica et Biophysica Acta (BBA)—Molecular Basis of Disease. 2014;1842(10):1993–2009.
14. Arner P. Resistin: yet another adipokine tells us that men are not mice. Diabetologia. 2005;48(11):2203–5. pmid:16193286
15. Lunney JK. Advances in Swine Biomedical Model Genomics. International journal of biological sciences. 2007;3(3):179–84. pmid:17384736
16. Gottlieb H, Lalich JJ. The Occurrence of Arteriosclerosis in the Aorta of Swine. The American Journal of Pathology. 1954;30(4):851–5. pmid:13180692
17. Rapacz J, Hasler-Rapacz J, Taylor K, Checovich W, Attie A. Lipoprotein mutations in pigs are associated with elevated plasma cholesterol and atherosclerosis. Science. 1986;234(4783):1573–7. pmid:3787263
18. Prescott MF, McBride CH, Hasler-Rapacz J, Von Linden J, Rapacz J. Development of complex atherosclerotic lesions in pigs with inherited hyper-LDL cholesterolemia bearing mutant alleles for apolipoprotein B. The American Journal of Pathology. 1991;139(1):139–47. pmid:1853929
19. Hamamdzic D, Wilensky RL. Porcine models of accelerated coronary atherosclerosis: role of diabetes mellitus and hypercholesterolemia. Journal of diabetes research. 2013;2013:761415. Epub 2013/07/12. pmid:23844374
20. Bollen PJ, Madsen LW, Meyer O, Ritskes-Hoitinga J. Growth differences of male and female Gottingen minipigs during ad libitum feeding: a pilot study. Lab Anim. 2005;39(1):80–93. Epub 2005/02/11. pmid:15703128
21. Johansen T, Hansen HS, Richelsen B, Malmlof R. The obese Gottingen minipig as a model of the metabolic syndrome: dietary effects on obesity, insulin sensitivity, and growth hormone profile. Comparative medicine. 2001;51(2):150–5. Epub 2002/04/02. pmid:11922179
22. Christoffersen B, Golozoubova V, Pacini G, Svendsen O, Raun K. The young Gottingen minipig as a model of childhood and adolescent obesity: influence of diet and gender. Obesity. 2013;21(1):149–58. Epub 2013/03/19. pmid:23505180
23. Cirera S, Jensen MS, Elbrond VS, Moesgaard SG, Christoffersen BO, Kadarmideen HN, et al. Expression studies of six human obesity-related genes in seven tissues from divergent pig breeds. Anim Genet. 2014;45(1):59–66. pmid:24033492
24. Christoffersen B, Ribel U, Raun K, Golozoubova V, Pacini G. Evaluation of different methods for assessment of insulin sensitivity in Gottingen minipigs: introduction of a new, simpler method. American journal of physiology Regulatory, integrative and comparative physiology. 2009;297(4):R1195–201. Epub 2009/08/28.
25. Kogelman LJ, Kadarmideen HN, Mark T, Karlskov-Mortensen P, Bruun CS, Cirera S, et al. 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. 2013;4:29. Epub 2013/03/22. pmid:23515185
26. Pant SD, Karlskov-Mortensen P, Jacobsen MJ, Cirera S, Kogelman LJ, Bruun CS, et al. Comparative Analyses of QTLs Influencing Obesity and Metabolic Phenotypes in Pigs and Humans. PLoS One. 2015;10(9):e0137356. PubMed Central PMCID: PMCPMC4562524. pmid:26348622
27. Madsen P, Jensen J, Labouriau R, Christensen O, Sahana G. DMU—A Package for Analyzing Multivariate Mixed Models in quantitative Genetics and Genomics. URL: pureaudk/portal/files/82154310/Paperpdf. 2014.
28. Wilson AJ, Reale D, Clements MN, Morrissey MM, Postma E, Walling CA, et al. An ecologist's guide to the animal model. The Journal of animal ecology. 2010;79(1):13–26. Epub 2010/04/23. pmid:20409158
29. Aulchenko YS, Ripke S, Isaacs A, van Duijn CM. GenABEL: an R library for genome-wide association analysis. Bioinformatics (Oxford, England). 2007;23(10):1294–6. Epub 2007/03/27.
30. Amin N, van Duijn CM, Aulchenko YS. A genomic background based method for association analysis in related individuals. PLoS One. 2007;2(12):e1274. Epub 2007/12/07. pmid:18060068
31. Aulchenko YS, de Koning DJ, Haley C. Genomewide rapid association using mixed model and regression: a fast and simple method for genomewide pedigree-based quantitative trait loci association analysis. Genetics. 2007;177(1):577–85. Epub 2007/07/31. pmid:17660554
32. Chen W-M, Abecasis GR. Family-Based Association Tests for Genomewide Association Scans. The American Journal of Human Genetics. 2007;81(5):913–26. pmid:17924335
33. Gao X. Multiple testing corrections for imputed SNPs. Genetic Epidemiology. 2011;35(3):154–8. pmid:21254223
34. Gao X, Becker LC, Becker DM, Starmer JD, Province MA. Avoiding the high Bonferroni penalty in genome-wide association studies. Genetic Epidemiology. 2010;34(1):100–5. pmid:19434714
35. Gao X, Starmer J, Martin ER. A multiple testing correction method for genetic association studies using correlated single nucleotide polymorphisms. Genet Epidemiol. 2008;32(4):361–9. Epub 2008/02/14. pmid:18271029
36. Falconer D, Mackay T. Introduction to Quantitative Genetics. 4 ed. England: Longman Group Limited; 1996.
37. Gabriel SB, Schaffner SF, Nguyen H, Moore JM, Roy J, Blumenstiel B, et al. The Structure of Haplotype Blocks in the Human Genome. Science. 2002;296(5576):2225–9. pmid:12029063
38. Barrett JC, Fry B, Maller J, Daly MJ. Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics (Oxford, England). 2005;21(2):263–5. Epub 2004/08/07.
39. Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, Bender D, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. American journal of human genetics. 2007;81(3):559–75. Epub 2007/08/19. pmid:17701901
40. Manunza A, Casellas J, Quintanilla R, Gonzalez-Prendes R, Pena RN, Tibau J, et al. A genome-wide association analysis for porcine serum lipid traits reveals the existence of age-specific genetic determinants. BMC genomics. 2014;15:758. Epub 2014/09/06. pmid:25189197
41. Yang H, Huang X, Zeng Z, Zhang W, Liu C, Fang S, et al. Genome-Wide Association Analysis for Blood Lipid Traits Measured in Three Pig Populations Reveals a Substantial Level of Genetic Heterogeneity. PLoS One. 2015;10(6):e0131667. Epub 2015/06/30. pmid:26121138
42. Helgadottir A, Gretarsdottir S, Thorleifsson G, Hjartarson E, Sigurdsson A, Magnusdottir A, et al. Variants with large effects on blood lipids and the role of cholesterol and triglycerides in coronary disease. Nat Genet. 2016;48(6):634–9. pmid:27135400
43. Ruiz-Velasco N, Domínguez A, Vega MA. Statins upregulate CD36 expression in human monocytes, an effect strengthened when combined with PPAR-γ ligands Putative contribution of Rho GTPases in statin-induced CD36 expression. Biochemical Pharmacology. 2004;67(2):303–13. pmid:14698043
44. Calvo D, Gomez-Coronado D, Suarez Y, Lasuncion MA, Vega MA. Human CD36 is a high affinity receptor for the native lipoproteins HDL, LDL, and VLDL. Journal of lipid research. 1998;39(4):777–88. Epub 1998/04/29. pmid:9555943
45. Roberson LL, Aneni EC, Maziak W, Agatston A, Feldman T, Rouseff M, et al. Beyond BMI: The "Metabolically healthy obese" phenotype & its association with clinical/subclinical cardiovascular disease and all-cause mortality—a systematic review. BMC Public Health. 2014;14:14. pmid:24400816
46. Hwang LC, Bai CH, Sun CA, Chen CJ. Prevalence of metabolically healthy obesity and its impacts on incidences of hypertension, diabetes and the metabolic syndrome in Taiwan. Asia Pacific journal of clinical nutrition. 2012;21(2):227–33. pmid:22507609
47. Hamer M, Stamatakis E. Metabolically healthy obesity and risk of all-cause and cardiovascular disease mortality. The Journal of clinical endocrinology and metabolism. 2012;97(7):2482–8. pmid:22508708
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
© 2017 Frederiksen et al. This is an open access article distributed under the terms of the Creative Commons Attribution License: http://creativecommons.org/licenses/by/4.0/ (the “License”), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. Notwithstanding the ProQuest Terms and Conditions, you may use this content in accordance with the terms of the License.
Abstract
We have established a pig resource population specifically designed to elucidate the genetics involved in development of obesity and obesity related co-morbidities by crossing the obesity prone Göttingen Minipig breed with two lean production pig breeds. In this study we have performed genome wide association (GWA) to identify loci with effect on blood lipid levels. The most significantly associated single nucleotide polymorphisms (SNPs) were used for linkage disequilibrium (LD) and haplotype analyses. Three separate haploblocks which influence the ratio between high density lipoprotein cholesterol and total cholesterol (HDL-C/CT), triglycerides (TG) and low density lipoprotein cholesterol (LDL-C) levels respectively were identified on Sus Scrofa chromosome 3 (SSC3). Large additive genetic effects were found for the HDL-C/CT and LDL-C haplotypes. Haplotypes segregating from Göttingen Minipigs were shown to impose a positive effect on blood lipid levels. Thus, the genetic profile of the Göttingen Minipig breed seems to support a phenotype comparable to the metabolic healthy obese (MHO) phenotype in humans.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer