Associations between gene polymorphisms in two crucial metabolic pathways and growth traits in pigs

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.

Transcriptome analysis reveals the potential roles of long non-coding RNAs in feed efficiency of chicken

Feed efficiency is an important economic trait and reduces the production costs per unit of animal product. Up to now, few studies have conducted transcriptome profiling of liver tissue in feed efficiency-divergent chickens (Ross vs native breeds). Also, molecular mechanisms contributing to differences in feed efficiency are not fully understood, especially in terms of long non-coding RNAs (lncRNAs). Hence, transcriptome profiles of liver tissue in commercial and native chicken breeds were analyzed. RNA-Seq data along with bioinformatics approaches were applied and a series of lncRNAs and target genes were identified. Furthermore, protein-protein interaction network construction, co-expression analysis, co-localization analysis of QTLs and functional enrichment analysis were used to functionally annotate the identified lncRNAs. In total, 2,290 lncRNAs were found (including 1,110 annotated, 593 known and 587 novel), of which 53 (including 39 known and 14 novel), were identified as differentially expressed genes between two breeds. The expression profile of lncRNAs was validated by RT-qPCR. The identified novel lncRNAs showed a number of characteristics similar to those of known lncRNAs. Target prediction analysis showed that these lncRNAs have the potential to act in cis or trans mode. Functional enrichment analysis of the predicted target genes revealed that they might affect the differences in feed efficiency of chicken by modulating genes associated with lipid metabolism, carbohydrate metabolism, growth, energy homeostasis and glucose metabolism. Some gene members of significant modules in the constructed co-expression networks were reported as important genes related to feed efficiency. Co-localization analysis of QTLs related to feed efficiency and the identified lncRNAs suggested several candidates to be involved in residual feed intake. The findings of this study provided valuable resources to further clarify the genetic basis of regulation of feed efficiency in chicken from the perspective of lncRNAs.

Local adaptations of Mediterranean sheep and goats through an integrative approach

Small ruminants are suited to a wide variety of habitats and thus represent promising study models for identifying genes underlying adaptations. Here, we considered local Mediterranean breeds of goats (n = 17) and sheep (n = 25) from Italy, France and Spain. Based on historical archives, we selected the breeds potentially most linked to a territory and defined their original cradle (i.e., the geographical area in which the breed has emerged), including transhumant pastoral areas. We then used the programs PCAdapt and LFMM to identify signatures of artificial and environmental selection. Considering cradles instead of current GPS coordinates resulted in a greater number of signatures identified by the LFMM analysis. The results, combined with a systematic literature review, revealed a set of genes with potentially key adaptive roles in relation to the gradient of aridity and altitude. Some of these genes have been previously implicated in lipid metabolism (SUCLG2, BMP2), hypoxia stress/lung function (BMPR2), seasonal patterns (SOX2, DPH6) or neuronal function (TRPC4, TRPC6). Selection signatures involving the PCDH9 and KLH1 genes, as well as NBEA/NBEAL1, were identified in both species and thus could play an important adaptive role.

Genetic variations (eQTLs) in muscle transcriptome and mitochondrial genes, and trans-eQTL molecular pathways in feed efficiency from Danish breeding pigs

Feed efficiency (FE) is a key trait in pig production, as improvement in FE has positive economic and environmental impact. FE is a complex phenotype and testing animals for FE is costly. Therefore, there has been a desire to find functionally relevant single nucleotide polymorphisms (SNPs) as biomarkers, to improve our biological understanding of FE as well as accuracy of genomic prediction for FE. We have performed a cis- and trans- eQTL (expression quantitative trait loci) analysis, in a population of Danbred Durocs (N = 11) and Danbred Landrace (N = 27) using both a linear and ANOVA model based on muscle tissue RNA-seq. We analyzed a total of 1425x19179 or 2.7x10⁷ Gene-SNP combinations in eQTL detection models for FE. The 1425 genes were from RNA-Seq based differential gene expression analyses using 25880 genes related to FE and additionally combined with mitochondrial genes. The 19179 SNPs were from applying stringent quality control and linkage disequilibrium filtering on genotype data using a GGP Porcine HD 70k SNP array. We applied 1000 fold bootstrapping and enrichment analysis to further validate and analyze our detected eQTLs. We identified 13 eQTLs with FDR < 0.1, affecting several genes found in previous studies of commercial pig breeds. Examples include MYO19, CPT1B, ACSL1, IER5L, CPT1A, SUCLA2, CSRNP1, PARK7 and MFF. The bootstrapping results showed statistically significant enrichment (p-value<2.2x10^-16) of eQTLs with p-value < 0.01 in both cis and trans-eQTLs. Enrichment analysis of top trans-eQTLs revealed high enrichment for gene categories and gene ontologies associated with genomic context and expression regulation. This included transcription factors (p-value = 1.0x10^-13), DNA-binding (GO:0003677, p-value = 8.9x10^-14), DNA-binding transcription factor activity (GO:0003700,) nucleus gene (GO:0005634, p-value<2.2x10^-16), negative regulation of expression (GO:0010629, p-value<2.2x10^-16). These results would be useful for future genome assisted breeding of pigs to improve FE, and in the improved understanding of the functional mechanism of trans eQTLs.

Genome-wide association study for carcass quality traits and growth in purebred and crossbred pigs1

Carcass quality traits such as back fat (BF), loin depth (LD), and ADG are of extreme economic importance for the swine industry. This study aimed to (i) estimate the genetic parameters for such traits and (ii) conduct a single-step genome-wide association study (ssGWAS) to identify genomic regions that affect carcass quality and growth traits in purebred (PB) and three-way crossbred (CB) pigs. A total of 28,497 PBs and 135,768 CBs pigs were phenotyped for BF, LD, and ADG. Of these, 4,857 and 3,532 were genotyped using the Illumina PorcineSNP60K Beadchip. After quality control, 36,328 SNPs were available and were used to perform an ssGWAS. A bootstrap analysis (n = 1,000) and a signal enrichment analysis were performed to declare SNP significance. Genome regions were based on the variance explained by significant 10-SNP sliding windows. Estimates of PB heritability (SE) were 0.42 (0.019) for BF, 0.39 (0.020) for LD, and 0.35 (0.021) for ADG. Estimates of CB heritability were 0.49 (0.042) for BF, 0.27 (0.029) for LD, and 0.12 (0.021) for ADG. Genetic correlations (SE) across the two populations were 0.81 (0.02), 0.79 (0.04), and 0.56 (0.05), for BF, LD, and ADG, respectively. The variance explained by significant regions for each trait in PBs ranged from 1.51% to 1.35% for BF, from 4.02% to 3.18% for LD, and from 2.26% to 1.45% for ADG. In CBs, the variance explained by significant regions ranged from 1.88% to 1.37% for BF, from 1.29% to 1.23% for LD, and from 1.54% to 1.32% for ADG. In this study, we have described regions of the genome that determine carcass quality and growth traits of PB and CB pigs. These results provide evidence that there are overlapping and nonoverlapping regions in the genome influencing carcass quality and growth traits in PBs and three-way CB pigs.

Health effects of kiwi wine on rats: an untargeted metabolic fingerprint study based on GC-MS/TOF

Kiwi wine is a popular fermentation product of kiwi fruit in Asian countries. To better understand the potential health effects of kiwi wine, an untargeted gas chromatography-mass spectrometer (GC-MS) approach was taken to assess the metabolic fingerprint of rats after dietary ingestion of kiwi wine. 7 differentially expressed endogenous metabolites from serum and 8 from urine were enriched in carbohydrate metabolism, amino acid metabolism pathway, fat metabolism and other metabolisms and selected from the KEGG. The above results showed that kiwi wine mainly led to a pronounced perturbation of energy metabolism (especially carbohydrate metabolism) during the consumption period. After stopping the supply of kiwi wine 30 days later, 6 and 3 endogenous metabolites from serum and urine respectively were screened and involved in a small part of carbohydrate related amino acid metabolism and fat metabolism, which indicated that the effect of kiwi wine sustained a lasting effect on energy metabolism, amino acid metabolism and lipid metabolism after stopping the supply. Thus, kiwi wine might have a positive function on health associated with the metabolism of its constituents. To the best of our knowledge, this study provides a nutrition field view for the development of the kiwi wine agricultural industry via an untargeted GC-MS metabolomic approach.

Influence of Delipation on the Energy Metabolism in Pig Parthenogenetically Activated Embryos

This study was designed not only to measure the effect of delipation on the developmental viability of pig parthenogenetically activated (PA) embryos, but also to evaluate the changes of mitochondria DNA (mtDNA), reactive oxygen species (ROS) level, adenosine triphosphate (ATP) content and gene (Acsl3, Acadsb, Acaa2, Glut1) expression level at different stages after delipation. Results showed that no effect was observed on the cleavage ability, but significant lower blastocyst rate was obtained in delipated embryos. Copy number of mtDNA decreased gradually from MII to four-cell stages and subsequently kept consistent with blastocyst stage both in delipated and control embryos, but the copy number of mtDNA in delipated embryos was similar to that in the control groups no matter at which developmental stage was observed. Both in delipated and control embryos, ATP content progressive decreased from one-cell to blastocyst stages, while just at one-cell stage, a significant decrease of ATP level was observed in delipated embryos compared with that of control. The level of ROS increased obviously after delipation at cleavage stage, but no difference was seen at blastocyst stage. Finally, the expression level of genes related to fatty acids beta-oxidation (Acadsb and Acaa2) was decreased, while the expression level of genes related to glucose metabolism (Glut 1) was upregulated after delipation. In conclusion, the reduction of lipids in pig oocytes will affect the developmental competence of pig PA embryos by disturbed energy metabolism and ROS stress.