Gene expression profiling reveals distinct features of various porcine adipose tissues
Lipids in Health and Disease volume 12, Article number: 75 (2013)
The excessive accumulation of body fat is a major risk factor to develop a variety of metabolic diseases. To investigate the systematic association between the differences in gene expression profiling and adipose deposition, we used pig as a model, and measured the gene expression profiling of six variant adipose tissues in male and females from three pig breeds which display distinct fat level.
We identified various differential expressed genes among breeds, tissues and between sexes, and further used a clustering method to identify sets of functionally co-expression genes linked to different obesity-related phenotypes. Our results reveal that the subcutaneous adipose tissues mainly modulate metabolic indicators, nonetheless, the visceral adipose tissues as well as the intermuscular adipose tissue were mainly associated with the impaired inflammatory and immune response.
The present study provided the evidence of gene expression profiling that the subcutaneous adipose tissues are mainly affected the metabolism process, whereas the visceral and intermuscular adipose tissues should been term as the metabolic risk factors of obesity.
Obesity is becoming a major global health hazard mainly because of the increased ingestion of high energy food and the decreased amounts of physical activity in populations. A gain in the adipose mass in the body is generally accompanied by increased risk of metabolism-related diseases. However, it is not only the amount of body fat mass, but also the distribution of adipose tissue that is now recognised as an important predictor of metabolic abnormalities.
Based on their anatomical location, adipose tissue can be divided into two major types: visceral adipose tissue (VAT) and subcutaneous adipose tissue (SAT). The VAT is a higher risk factor for obesity than SAT. In particular, basal free fatty acid flux, lipolysis rate, and secretory protein expression are markedly higher in visceral adipocytes compared with subcutaneous adipocytes. Our previous reports have provided epigenetic evidence that VAT is mainly associated with impaired inflammatory and immune responses and is a metabolic risk factor of obesity, while SAT is mainly associated with metabolic processes[3, 4]. In addition, ethnic group and gender are important factors that affect adipose deposition[5, 6].
Pig (Sus scrofa) has been recognised as an attractive biomedical model for human because its physiology and anatomy are remarkably similar to humans and pigs are easy to breed and handle. Here, we present a comprehensive genome-wide comparison of gene expression profiling among six adipose tissues from different body sites in both sexes of three well-defined pig breeds displaying distinct adipose phenotypes. Using microarray technology, we identified various differentially expressed genes that may be related to the different adipose deposition phenotypes among the breeds, and between the sexes and anatomic locations. We also used a clustering method to identify sets of co-expressed genes that may be functionally linked to different obesity-related phenotypes. These results can be used to understand the different accumulation mechanisms of the distinct adipose tissue types.
Results and discussion
Our previous report based on same individuals have demonstrated that the adipose volume were significant difference among adipose tissues and breeds, and between sexes. In addition, the 24 kinds of reprehensive metabolism indicators in serum also revealed the same ranking from the leaner Landrace, the wild Tibetan and the fatty Rongchang pigs.
Identification of differentially expressed (DE) genes and relatedness of adipose samples
The number of DE genes (P < 0.05, three-way repeated-measures analysis of variance [ANOVA]) was lowest between male and female (1,225), higher among different adipose tissues (2,750), and highest among the different breeds (2,949), indicating significant biological differences between the groups in the latter two categories (Figure 1A). To further characterise variability in the gene expression profiling, we performed a principal components analysis (PCA) based on the DE genes in each of the compared groups. The six adipose tissues clustered into two groups, VATs and SATs, which reflected their distinct functional and metabolic features (Figure 1B). Interestingly, the intermuscular adipose (IAD) tissue that is deposited between muscle bundles, clustered with the VATs in terms of gene expression profiling, suggesting that IAD may be a potential new risk factor for obesity-related diseases, in agreement with our findings from a previous study on methylation.
Although the three pig breeds clustered into three groups, the Rongchang and Tibetan breeds were closer to each other than the Landrace breed (Figure 1C). These results not only reflect the marked phenotypic differences in adipose deposition among the three pig breeds which have undergone different/opposite breeding directions, but also highlight the deep phylogenetic split between the European (Landrace) and Asian pigs (Rongchang and Tibetan).
In addition, the male and female were also been roughly distinguish (Figure 1D), which consisted with the lowest variation of gene expression profiling between sexes, which only have less than half of the amount of DE genes than among tissue types and breeds.
DE genes between the VATs and SATs
The PCA classified six adipose tissues into two groups: the VAT group which included greater omentum adipose (GOM), mesenteric adipose (MAD), retroperitoneal adipose (RAD) and IAD; and the SAT group which included upper layer of backfat (ULB) and inner layer of backfat (ILB).
We performed a pairwise comparison analysis and identified 192 and 181 DE genes that were specifically up-regulated in the SAT and VAT groups, respectively (Figure 2 and Table 1). The DE genes which were specifically up-regulated in the VAT group were mainly enriched in the GO categories of the immune processes, including positive regulation of defence response (8 genes, P = 5.63 × 10-4), positive regulation of immune system process (10 genes, P = 1.11 × 10-4), adaptive immune response (6 genes, P = 3.54 × 10-4), innate immune response (7 genes, P = 7.86 × 10-4) and B cell mediated immunity (5 genes, P = 9.91 × 10-4) (Figure 3). Various genes known to be involved in the immune response were highly expressed in VATs compared with in SATs (P < 0.05). For example, C3 (complement 3), which was reported to play a central role in the activation of the complement system[8, 9] was up-regulated in the VATs. CLEC12B (C-type lectin domain family 12B) limits the activity of monocyte-derived immune cells after cell differentiation and possibly during inflammatory diseases (Figure 4). The action of these genes which were highly expressed in VATs may contribute to the obesity-induced chronic inflammation in adipose tissue that precedes the development of insulin resistance and type II diabetes.
The DE genes that were specifically up-regulated in the SAT group were mainly enriched in the GO categories of metabolic processes, including glucose metabolic process (6 genes, P = 1.16 × 10-2), hexose metabolic process (6 genes, P = 2.81 × 10-2), fatty acid metabolic process (6 genes, P = 3.15 × 10-2) and monosaccharide metabolic process (6 genes, P = 4.77 × 10-2) (Figure 3). This result is consistent with previous reports that, compared with the high-risk of visceral obesity caused by VATs, SATs are mainly associated with metabolic processes[4, 11].
DE genes among the three pig breeds
The number of DE genes (P < 0.05, three-way repeated-measures ANOVA) was highest between the fatty Rongchang and leaner Landrace pigs (1,855), lower between the wild Tibetan and Landrace pigs (1,440), and lowest between the Tibetan and Rongchang pigs (868) (Figure 2).
The highly expressed genes in the Landrace pigs were mainly involved in the GO and KEGG categories of skeletal system development processes and immune processes, including positive regulation of skeletal system development (7 genes, P = 2.22 × 10-4), bone development (8 genes, P = 4.16 × 10-2), positive regulation of immune system process (20 genes, P = 1.19 × 10-3) and adaptive immune response (8 genes, P = 2.39 × 10-2) (Figure 5). The highly expressed genes in the Rongchang pigs were mainly involved in the processes of lipid metabolism and immune processes, such as lipid oxidation (8 genes, P = 1.40 × 10-3), fatty acid oxidation (8 genes, P = 2.53 × 10-3), fatty acid metabolic process (21 genes, P = 4.71 × 10-4), activation of immune response (6 genes, P = 1.30 × 10-2) and positive regulation of immune response (8 genes, P = 5.93 × 10-3). This result reflects the phenotypic divergence of adipose deposition between these two pig breeds that have undergone artificial selection for different breeding targets. The Landrace breed is under strong artificial selection for more muscle production, while the Rongchang breed is selected for more adipose deposition.
DE genes between the sexes
Almost three-quarters of the DE genes between the sexes (879 of 1,225, 71.76%) were up-regulated in the males (Figure 2). These highly expressed genes were enriched mainly in the GO categories of the immune processes (Figure 6), including positive regulation of immune system process (28 genes, P = 4.49 × 10-6), adaptive immune response (11 genes, P = 1.71 × 10-3), humoral immune response (7 genes, P = 1.86 × 10-3), immunoglobulin mediated immune response (9 genes, P = 2.09 × 10-3), regulation of acute inflammatory response (5 genes, P = 1.15 × 10-2) and immune effector process (13 genes, P = 1.33 × 10-2). The DE genes (346 of 1,225, 28.24%), which were up-regulated in the females (Figure 2), were mainly involved in metabolic processes (Figure 6), including the coenzyme metabolic process (10 genes, P = 1.08 × 10-3), cofactor metabolic process (11 genes, P = 1.64 × 10-3), oxidoreduction coenzyme metabolic process (5 genes, P = 1.54 × 10-2) and organic acid biosynthetic process (8 genes, P = 1.10 × 10-2). This result reflects a sexual dimorphism in body fat distribution. Females deposit relatively more fat in subcutaneous/inguinal depots, which has beneficial effects on the control of body weight and metabolism; on the other hand, males tend to deposit more fat in intra-abdominal/gonadal depots, which has been found to be related to risk for a number of diseases.
Interestingly, we identified some genes that potentially contribute to the sexual differences in obesity development. Compared with the females, the males exhibited higher expression levels of RXFP2 (P =1.19 × 10-13) and CPB1 (P = 5.69 × 10-6) (Figure 7), which have been related to gubernaculum testis development and spermatogenesis[15, 16]. However, compared with the males, the females exhibited higher expression levels of PAEP (P = 3.71 × 10-5) and FABP7 (P = 8.42 × 10-9) (Figure 7), which have been associated mainly with the secretory endometrial glands and breast cancer[17, 18].
Validation of gene expression changes by Q-PCR
The expression changes of nine genes that are well-known to be related to adipose metabolism (ACADL, ADIPOQ, ADIPOR 1, CAV 3, MDH 1, MDH 2, ME 1, Insig 1 and SCAP) showed significant positive correlations between the Q-PCR results and their expressions determined by microarray (P < 0.05, Figure 8), which highlighted the high confidence of the results from the microarray approach.
Co-expressed gene sets linked to the obesity-related phenotypes
Generally, quantitative characters (such as adipose deposition) in mammals are determined by a set of genes that cooperate in certain biological processes. To further explore the relationships between the genes and obesity-related traits, we applied a clustering method to identify sets of functionally related genes linked to the phenotypic traits as reported previously.
We identified eight gene modules containing 947 genes for ULB and twelve gene modules containing 1,277 genes for ILB. In three typical VATs, we identified nine, nine and eleven gene modules containing 933, 1,058 and 1,409 genes for the MAD, GOM and RAD tissues, respectively. For the IAD, we identified eleven modules containing 1,423 genes (Figure 9).
Next, we explored the correlations between the modules of co-expressed genes and the various obesity phenotypes. Three gene modules in the SATs (i.e. ULB and ILB) were significantly correlated with adipocyte volume and the concentration of plasma lactate dehydrogenase in serum (false discovery rate adjusted permutation, P < 0.05) (Figure 10 and Table 2). The co-expressed genes within these modules were mainly involved in metabolic processes, such as fatty acid metabolic process (23 genes, P = 4.38 × 10-5), fatty acid beta-oxidation (3 genes, P = 2.65 × 10-3), coenzyme metabolic process (6 genes, P = 4.77 × 10-5) and carbohydrate metabolic process (12 genes, P = 2.94 × 10-3) (Table 3), in agreement with an earlier report that SATs contribute mainly to metabolic processes[20, 21].
The co-expressed genes in two of the VATs (i.e. GOM and MAD) mainly affected the concentrations of inflammatory and immune adipokines in serum, i.e. plasminogen activator inhibitor-1 (PAI-1), interleukin 6 (IL-6) and tumor necrosis factor-α (TNF-α) (Figure 10 and Table 2). The co-expressed genes within these modules were mainly involved in immune system process (21 genes, P = 6.66 × 10-4), immune response (8 genes, P = 6.84 × 10-3), cellular defence response (15 genes, P =1.25 × 10-10), cellular defence response (6 genes, P = 4.67 × 10-3) and B cell mediated immunity (9 genes, P = 1.05 × 10-7). This result suggested that the GOM and MAD tissues contributed mainly to the obesity-induced chronic inflammation in adipose tissue that precedes the development of insulin resistance and type II diabetes[22, 23]. The RAD tissue, on the other hand, was found to affect mainly the body density and the concentrations of three metabolic indicators in serum, namely, high density lipoprotein, insulin, and apolipoprotein A-1 (Apo-A1), which are mainly involved in metabolic processes such as coenzyme metabolic process (16 genes, P =3.02 × 10-9), primary metabolic process (241 genes, P = 9.31 × 10-9), and tricarboxylic acid cycle (7 genes, P = 5.49 × 10-6) (Table 3). Studies in human showed that GOM was strongly correlated to insulin, while no such correlations were found between lipid uptake and SAT or RAD.
Similar to the findings for the high-risk VATs, the co-expressed genes in IAD were also found to be significantly correlated with immune-related TNF-α (Figure 10 and Table 2), further suggesting that the IAD is an independent risk factor for metabolic diseases.
The present study describes a genome-wide analysis of gene expression profiling among various adipose depots between the sexes of three well-defined pig breeds displaying distinct adipose phenotypes. We present evidence that the SATs mainly affect metabolic processes, while the VATs and the IAD are associated mainly with immune and inflammation responses that are seen as the metabolic risk factors of obesity. This study also supports pig as an ideal model for studying human obesity.
Materials and methods
Three females and three males at 210-days-old for each of the leaner Landrace pigs (a Western breed), the wild Tibetan pigs (a feral, indigenous Chinese pig that has not undergone artificial selection) and the fatty Rongchang pigs (a Chinese breed) were used in this study as previously described. Animals were humanely sacrificed, in compliance with experimental animals established by the Ministry of Agriculture of China. Three VATs (i.e. GOM, MAD and RAD), two SATs (i.e. ULB and ILB) as well as IAD were rapidly separated from each carcass, immediately frozen in liquid nitrogen, and stored at −80°C until RNA extraction. For more information, please refer to Li et al..
Measurements of the obesity-related phenotypes
Measurements of the pig body density, the concentrations of 24 serum-circulating indicators of metabolism and the adipocyte volume are from our previous report based on same individuals. For more information, please refer to Li et al..
RNA isolation and microarray analysis
The total RNA (10 μg) of the 108 tissue samples was extracted with TRIzol (Invitrogen) and further purified using an Rneasy column (Qiagen). Hybridization, washing, and scanning were done according to standard Agilent protocols. The data analysis was performed with MultiExperiment Viewer (MeV). Array data have been uploaded to NCBI’s Gene Expression Omnibus (GEO) [accession number GSE30343]. For more information, please refer to Li et al..
Identification of differentially expressed genes
We used three-way repeated-measures ANOVA (n = 3 per breed per sex per tissue) to identified differentially expressed (DE) mRNA among six adipose tissues between sexes of three pig breeds (P < 0.05, corrected with adjusted Bonferroni method, FDR <0.05, 1,000 permutations), which were further analyzed using DAVID to examine whether these DE genes were enriched for specific functional catalogs (i.e. KEGG pathways and Gene Ontology-Biological Processes and -Molecular Functions) (P < 0.05, corrected with adjusted Bonferroni method).
Identification of co-expressed gene modules
To identify sets of functionally related genes in six adipose tissues, linked to the phenotypic traits, we applied a clustering method as described previously[4, 19] with some modifications. Spearman rank correlation coefficients were determined between all possible gene-pairs across 18 samples for each adipose tissue. The strongest correlated gene-pair was selected, and grouped together in a set that was assigned the average expression value of the two genes that constitute this set. After addition of this newly created set to the dataset, the two individual genes were removed from the data and the strongest correlation in the dataset was again selected. This resulted in either the expansion of a set already created or in the creation of a new set. We kept repeating this as an iterative process until the most significantly correlated pair was coefficient of correlation r < 0.80. Only the sets containing 50 or more genes are kept for further analysis. The gene sets - reflecting the average expression value of the genes constituting that set - were correlated with obesity phenotypic traits using a Spearman rank correlation coefficient.
Quantitative PCR (Q-PCR) was used to confirm the expression pattern observed in the microarray results. Q-PCR was performed on a CF96 Real-Time PCR Detection System (Bio-Rad) using SYBR® Green Real-time PCR Master Mix (TaKaRa, China). The PCR primer sequences are shown in Table 4. Porcine ACTB, TBP and TOP2B were simultaneously used as endogenous control genes. The 2-ΔΔCt method was used to determine the relative mRNA abundance for the surveyed samples.
Subcutaneous adipose tissue
Visceral adipose tissue
Upper layer of backfat
Inner layer of backfat
Greater omentum adipose
Wronska A, Kmiec Z: Structural and biochemical characteristics of various white adipose tissue depots. Acta Physiol (Oxf). 2012, 205: 194-208. 10.1111/j.1748-1716.2012.02409.x.
Matsuzawa Y: The metabolic syndrome and adipocytokines. FEBS Lett. 2006, 580: 2917-2921. 10.1016/j.febslet.2006.04.028
Li M, Wu H, Luo Z, Xia Y, Guan J, Wang T, Gu Y, Chen L, Zhang K, Ma J: An atlas of DNA methylomes in porcine adipose and muscle tissues. Nat Commun. 2012, 3: 850-861.
Li M, Wu H, Wang T, Xia Y, Jin L, Jiang A, Zhu L, Chen L, Li R, Li X: Co-methylated genes in different adipose depots of pig are associated with metabolic, inflammatory and immune processes. Int J Biol Sci. 2012, 8: 831-837.
Karastergiou K, Smith SR, Greenberg AS, Fried SK: Sex differences in human adipose tissues–the biology of pear shape. Biology of Sex Differences. 2012, 3: 13- 10.1186/2042-6410-3-13
Deurenberg P, Deurenberg-Yap M: Differences in body-composition assumptions across ethnic groups: practical consequences. Curr Opin Clin Nutr. 2001, 4: 377-383. 10.1097/00075197-200109000-00007.
Lunney JK: Advances in swine biomedical model genomics. Int J Biol Sci. 2007, 3: 179-184.
Mueller-Ortiz SL, Drouin SM, Wetsel RA: The alternative activation pathway and complement component C3 are critical for a protective immune response against Pseudomonas aeruginosa in a murine model of pneumonia. Infect Immun. 2004, 72: 2899-2906. 10.1128/IAI.72.5.2899-2906.2004
Abdelnoor AM, Kobeissy F, Farhat D, Hadi U: Circulating immune complexes and complement C3 and C4 levels in a selected group of patients with rhinitis in Lebanon. Clin Mol Allergy. 2004, 2: 6- 10.1186/1476-7961-2-6
Hoffmann SC, Schellack C, Textor S, Konold S, Schmitz D, Cerwenka A, Pflanz S, Watzl C: Identification of CLEC12B, an inhibitory receptor on myeloid cells. J Biol Chem. 2007, 282: 22370-22375. 10.1074/jbc.M704250200
Wolfs MG, Rensen SS, Bruin-Van Dijk EJ, Verdam FJ, Greve J-W, Sanjabi B, Bruinenberg M, Wijmenga C, van Haeften TW, Buurman WA: Co-expressed immune and metabolic genes in visceral and subcutaneous adipose tissue from severely obese individuals are associated with plasma HDL and glucose levels: a microarray study. BMC Med Genomics. 2010, 3: 34- 10.1186/1755-8794-3-34
Wajchenberg BL: Subcutaneous and visceral adipose tissue: their relation to the metabolic syndrome. Endocr Rev. 2000, 21: 697-738. 10.1210/er.21.6.697
Satoor SN, Puranik AS, Kumar S, Williams MD, Ghale M, Rahalkar A, Karandikar MS, Shouche Y, Patole M, Bhonde R: Location, location, location: Beneficial effects of autologous fat transplantation. Sci Rep. 2011, 1: 81-
Grove K, Fried S, Greenberg A, Xiao X, Clegg D: A microarray analysis of sexual dimorphism of adipose tissues in high-fat-diet-induced obese mice. Int J Obesity. 2010, 34: 989-1000. 10.1038/ijo.2010.12.
Yuan F, Li X, Lin J, Schwabe C, Büllesbach E, Rao C, Lei Z: The role of RXFP2 in mediating androgen-induced inguinoscrotal testis descent in LH receptor knockout mice. Reproduction. 2010, 139: 759-769. 10.1530/REP-09-0518
Luitjens C, Gallegos M, Kraemer B, Kimble J, Wickens M: CPEB proteins control two key steps in spermatogenesis in C. elegans. Genes Dev. 2000, 14: 2596-2609. 10.1101/gad.831700
Joshi SG, Bank JF, Henriques ES, Makarachi A, Matties G: Serum levels of a progestagen-associated endometrial protein during the menstrual cycle and pregnancy. J Clin Endocr Metab. 1982, 55: 642-648. 10.1210/jcem-55-4-642
Zhang H, Rakha E, Ball G, Spiteri I, Aleskandarany M, Paish E, Powe D, Macmillan R, Caldas C, Ellis I: The proteins FABP7 and OATP2 are associated with the basal phenotype and patient outcome in human breast cancer. Breast Cancer Res Treat. 2010, 121: 41-51. 10.1007/s10549-009-0450-x
Eisen MB, Spellman PT, Brown PO, Botstein D: Cluster analysis and display of genome-wide expression patterns. Proc Natl Acad Sci. 1998, 95: 14863-14868. 10.1073/pnas.95.25.14863
Smith SR, Lovejoy JC, Greenway F, Ryan D, deJonge L, de la Bretonne J, Volafova J, Bray GA: Contributions of total body fat, abdominal subcutaneous adipose tissue compartments, and visceral adipose tissue to the metabolic complications of obesity. Metabolism. 2001, 50: 425-435. 10.1053/meta.2001.21693
Tran TT, Yamamoto Y, Gesta S, Kahn CR: Beneficial effects of subcutaneous fat transplantation on metabolism. Cell Metab. 2008, 7: 410-420. 10.1016/j.cmet.2008.04.004
Iyer A, Fairlie DP, Prins JB, Hammock BD, Brown L: Inflammatory lipid mediators in adipocyte function and obesity. Nat Rev Endocrinol. 2010, 6: 71-82. 10.1038/nrendo.2009.264
Torres‒Leal FL, Fonseca‒Alaniz MH, Rogero MM, Tirapegui J: The role of inflamed adipose tissue in the insulin resistance. Cell Biochem Funct. 2010, 28: 623-631. 10.1002/cbf.1706
Maury E, Ehala-Aleksejev K, Guiot Y, Detry R, Vandenhooft A, Brichard SM: Adipokines oversecreted by omental adipose tissue in human obesity. Am J Physiol-Endoc M. 2007, 293: E656-E665.
Yim J, Heshka S, Albu J, Heymsfield S, Kuznia P, Harris T, Gallagher D: Intermuscular adipose tissue rivals visceral adipose tissue in independent associations with cardiovascular risk. Int J Obesity. 2007, 31: 1400-1405. 10.1038/sj.ijo.0803621.
Da Wei Huang BTS, Lempicki RA: Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2008, 4: 44-57. 10.1038/nprot.2008.211.
This work was supported by grants from the National High Technology Research and Development Program of China (863 Program) (2013AA102502),the Specialized Research Fund of Ministry of Agriculture of China (NYCYTX-009), the Project of Provincial Twelfth Five Years’ Animal Breeding of Sichuan Province (2011YZGG15) and the National Special Foundation for Transgenic Species of China (2011ZX08006-003).
The authors have declared that no competing interests exist.
MZL and XWL conceived and designed the experiment. CWZ and MZL performed the data analysis and drafted the manuscript. JDM and AAJ collected the samples and records. GQT, MMM, LZ and LB carried out the physiological, biochemical and histological studies, statistical analysis and prepared nucleic acids. CWZ and JZ performed gene expression microarray and data analysis. All authors read and approved the final manuscript.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Zhou, C., Zhang, J., Ma, J. et al. Gene expression profiling reveals distinct features of various porcine adipose tissues. Lipids Health Dis 12, 75 (2013). https://doi.org/10.1186/1476-511X-12-75