Epigenetic and Transcriptomic Characterization of Pure Adipocyte Fractions From Obese Pigs Identifies Candidate Pathways Controlling Metabolism

Identification of potential obese-specific biomarkers and pathways associated with abdominal subcutaneous fat deposition in pig using a comprehensive bioinformatics strategy

Abdominal subcutaneous fat deposition (ASFD) is not only related to meat quality in the pig industry but also to human health in medicine. It is of great value to elucidate the potential molecular mechanisms of ASFD. The present study aims to identify obese-specific biomarkers and key pathways correlated with ASFD in pigs. The ASF-related mRNA expression dataset GSE136754 was retrieved from the Gene Expression Omnibus (GEO) database and systematically analyzed using a comprehensive bioinformatics method. A total of 565 differentially expressed genes (DEGs) were identified between three obese and three lean pigs, and these DEGs were mainly involved in the p53 signaling pathway, MAPK signaling pathway and fatty acid metabolism. A protein-protein interaction (PPI) network, consisting of 540 nodes and 1,065 edges, was constructed, and the top ten genes with the highest degree scores-ABL1, HDAC1, CDC42, HDAC2, MRPS5, MRPS10, MDM2, JUP, RPL7L1 and UQCRFS1-were identified as hub genes in the whole PPI network. Especially HDAC1, MDM2, MRPS10 and RPL7L1 were identified as potential robust obese-specific biomarkers due to their significant differences in single gene expression levels and high ROC area; this was further verified by quantitative real-time PCR (qRT-PCR) on abdominal subcutaneous fat samples from obese-type (Saba) and lean-type (Large White) pigs. Additionally, a mRNA-miRNA-lncRNA ceRNA network consisting of four potential biomarkers, 15 miRNAs and 51 lncRNAs was established, and two targeted lncRNAs with more connections, XIST and NEAT1, were identified as potentially important regulatory factors. The findings of this study may provide novel insights into the molecular mechanism involved in ASFD.

Elucidating breed-specific variants of native pigs in Korea: insights into pig breeds' genomic characteristics

Although conserving native pig breeds is important in Korea, research on the genomic aspects to identify breed-specific variations in native pig breeds is uncommon. Single nucleotide polymorphisms (SNPs) can be a powerful source for identifying breed-specific variants. We used whole genome sequencing data, including Jeju Native Pig (JNP), Korean Native Pig (KNP), Korean Wild Boar (KWB), and other western commercial pig breeds to determine native pig breed-specific SNPs. Furthermore, the goal was not only to determine the genomic specificity of native pig breeds but also to identify SNPs that carry breed-specific information (breed-informative SNPs) that can be related to breed characteristics. The representative characteristics of native pigs are their unique meat quality and disease resistance. We surveyed the gene ontology (GO) of native pigs with breed-specific SNPs. Examining the genes associated with GO may contribute to revealing the reasons for the unique characteristics of native pig breeds. The enriched GOs terms were neuron projection development, cell surface receptor signaling pathway, ion homeostasis in JNP, cell adhesion and wound healing in KNP, and DNA repair and reproduction in KWB. We expect that this study of breed-specific SNPs will enable us to gain a deeper understanding of native pigs in Korea.

SNP discovery and association study for growth, fatness and meat quality traits in Iberian crossbred pigs

Iberian pigs and its crosses are produced to obtain high-quality meat products. The objective of this work was to evaluate a wide panel of DNA markers, selected by biological and functional criteria, for association with traits related to muscle growth, fatness, meat quality and metabolism. We used 18 crossbred Iberian pigs with divergent postnatal growth patterns for whole genome sequencing and SNP discovery, with over 13 million variants being detected. We selected 1023 missense SNPs located on annotated genes and showing different allele frequencies between pigs with makerdly different growth patterns. We complemented this panel with 192 candidate SNPs obtained from literature mining and from muscle RNAseq data. The selected markers were genotyped in 480 Iberian × Duroc pigs from a commercial population, in which phenotypes were obtained, and an association study was performed for the 1005 successfully genotyped SNPs showing segregation. The results confirmed the effects of several known SNPs in candidate genes (such as LEPR, ACACA, FTO, LIPE or SCD on fatness, growth and fatty acid composition) and also disclosed interesting effects of new SNPs in less known genes such as LRIG3, DENND1B, SOWAHB, EPHX1 or NFE2L2 affecting body weight, average daily gain and adiposity at different ages, or KRT10, NLE1, KCNH2 or AHNAK affecting fatness and FA composition. The results provide a valuable basis for future implementation of marker-assisted selection strategies in swine and contribute to a better understanding of the genetic architecture of relevant traits.

Modeling microRNA-driven post-transcriptional regulation using exon-intron split analysis in pigs

The contribution of microRNAs (miRNAs) to mRNA post-transcriptional regulation has often been explored by the post hoc selection of downregulated genes and determining whether they harbor binding sites for miRNAs of interest. This approach, however, does not discriminate whether these mRNAs are also downregulated at the transcriptional level. Here, we have characterized the transcriptional and post-transcriptional changes in mRNA expression in two porcine tissues: gluteus medius muscle of fasted and fed Duroc gilts and adipose tissue of lean and obese Duroc-Göttingen minipigs. Exon-intron split analysis of RNA-seq data allowed us to identify downregulated mRNAs with high post-transcriptional signals in fed or obese states, and we assessed whether they harbor binding sites for upregulated miRNAs in any of these two physiological states. We found 26 downregulated mRNAs with high post-transcriptional signals in the muscle of fed gilts and 21 of these were predicted targets of miRNAs upregulated in fed pigs. For adipose tissue, 44 downregulated mRNAs in obese minipigs displayed high post-transcriptional signals, and 25 of these were predicted targets of miRNAs upregulated in the obese state. These results suggest that the contribution of miRNAs to mRNA repression is more prominent in the skeletal muscle system. Finally, we identified several genes that may play relevant roles in the energy homeostasis of the pig skeletal muscle (DKK2 and PDK4) and adipose (SESN3 and ESRRG) tissues. By differentiating transcriptional from post-transcriptional changes in mRNA expression, exon-intron split analysis provides a valuable view of the regulation of gene expression, complementary to canonical differential expression analyses.

Acetate rescues defective brain-adipose metabolic network in obese Wistar rats by modulation of peroxisome proliferator-activated receptor-γ

We investigated the hypothesis that acetate ameliorates brain-adipose metabolic dysfunction (BAMED) in high fat diet (HFD)-induced obesity, possibly by modulation of peroxisome proliferator-activated receptor-γ (PPAR-γ). Ten-week-old male Wistar rats were randomly assigned into four groups (n = 6/group): Control, acetate and obese with or without acetate groups received vehicle (distilled water; po), acetate (200 mg/kg, po) and 40% HFD with or without acetate respectively. The treatments lasted for 12 weeks. Obese animals showed increase in body weight, visceral fat mass, insulin and triglyceride-glucose index and a reduction in insulin sensitivity. In addition, obese animals also showed increase in plasma/hypothalamic and adipose pyruvate dehydrogenase kinase-4, lactate-pyruvate ratio, malondialdehyde, γ-glutamyl transferase, and a decrease in glucose-6-phosphate dehydrogenase, glutathione, nitric oxide and PPAR-γ. HFD also elevated plasma/hypothalamic lipid and decreased adipose lipid profile, increased hypothalamic and adipose tumor necrosis factor-α, interleukin-6 and histone deacetylase (HDAC), and elevated plasma/adipose leptin. These alterations were reversed by concomitant administration of acetate. The present results demonstrate that obesity is characterized by BAMED, which is accompanied by altered HDAC/PPAR-γ. The results in addition suggest that acetate, an HDAC inhibitor rescues BAMED with consequent normalization of body weight and visceral fat mass by modulation of PPAR-γ and suppression of oxidative stress.

Large-scale transcriptome sequencing in broiler chickens to identify candidate genes for breast muscle weight and intramuscular fat content

BACKGROUND

In broiler production, breast muscle weight and intramuscular fat (IMF) content are important economic traits. Understanding the genetic mechanisms that underlie these traits is essential to implement effective genetic improvement programs. To date, genome-wide association studies (GWAS) and gene expression analyses have been performed to identify candidate genes for these traits. However, GWAS mainly detect associations at the DNA level, while differential expression analyses usually have low power because they are typically based on small sample sizes. To detect candidate genes for breast muscle weight and IMF contents (intramuscular fat percentage and relative content of triglycerides, cholesterol, and phospholipids), we performed association analyses based on breast muscle transcriptomic data on approximately 400 Tiannong partridge chickens at slaughter age.

RESULTS

First, by performing an extensive simulation study, we evaluated the statistical properties of association analyses of gene expression levels and traits based on the linear mixed model (LMM) and three regularized linear regression models, i.e., least absolute shrinkage and selection operator (LASSO), ridge regression (RR), and elastic net (EN). The results show that LMM, LASSO and EN with tuning parameters that are determined based on the one standard error rule exhibited the lowest type I error rates. Using results from all three models, we detected 43 candidate genes with expression levels that were associated with breast muscle weight. In addition, candidate genes were detected for intramuscular fat percentage (1), triglyceride content (2), cholesterol content (1), and phospholipid content (1). Many of the identified genes have been demonstrated to play roles in the development and metabolism of skeletal muscle or adipocyte. Moreover, weighted gene co-expression network analyses revealed that many candidate genes were harbored by gene co-expression modules, which were also significantly correlated with the traits of interest. The results of Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analyses indicated that these modules are involved in muscle development and contraction, and in lipid metabolism.

CONCLUSIONS

Our study provides valuable insight into the transcriptomic bases of breast muscle weight and IMF contents in Chinese indigenous yellow broilers. Our findings could be useful for the genetic improvement of these traits in broiler chickens.

Diet-Dependent Changes of the DNA Methylome Using a Göttingen Minipig Model for Obesity

Objective: Environmental factors can influence obesity by epigenetic mechanisms. The aim of this study was to investigate obesity-related epigenetic changes and the potential for reversal of these changes in the liver of Göttingen minipigs subjected to diet interventions. Methods: High-throughput liquid hybridization capture-based bisulfite sequencing (LHC-BS) was used to quantify the methylation status of gene promotor regions in liver tissue in three groups of male castrated Göttingen minipigs: a standard chow group (SD, N = 7); a group fed high fat/fructose/cholesterol diet (FFC, N = 10) and a group fed high fat/fructose/cholesterol diet during 7 months and reversed to standard diet for 6 months (FFC/SD, N = 12). Expression profiling by qPCR of selected metabolically relevant genes was performed in liver tissue from all pigs. Results: The pigs in the FFC diet group became morbidly obese. The FFC/SD diet did not result in a complete reversal of the body weight to the same weight as in the SD group, but it resulted in reversal of all lipid related metabolic parameters. Here we identified widespread differences in the patterning of cytosine methylation of promoters between the different feeding groups. By combining detection of differentially methylated genes with a rank-based hypergeometric overlap algorithm, we identified 160 genes showing differential methylation in corresponding promoter regions in the FFC diet group when comparing with both the SD and FFC/SD groups. As expected, this differential methylation under FFC diet intervention induced de-regulation of several metabolically-related genes involved in lipid/cholesterol metabolism, inflammatory response and fibrosis generation. Moreover, five genes, of which one is a fibrosis-related gene (MMP9), were still perturbed after diet reversion. Conclusion: Our findings highlight the potential of exploring diet-epigenome interactions for treatment of obesity.