Genome-wide association studies for fatty acid metabolic traits in five divergent pig populations

Multiple omics analysis reveals the regulation of SIRT5 on mitochondrial function and lipid metabolism during the differentiation of bovine preadipocytes

Preadipocyte differentiation represents a critical stage in adipogenesis, with mitochondria playing an undeniable pivotal role. Given the intricate interplay between transcription and metabolic signaling during adipogenesis, the regulation of sirtuin 5 (SIRT5) on mitochondrial function and lipid metabolism was revealed via multiple omics analysis. The findings suggest that SIRT5 plays a crucial role in promoting mitochondrial biosynthesis and maintaining mitochondrial function during preadipocyte differentiation. Moreover, SIRT5 modulates the metabolic levels of numerous bioactive substances by extensively regulating genes expression associated with differentiation, energy metabolism, lipid synthesis, and mitochondrial function. Finally, SIRT5 was found to suppress triacylglycerols (TAG) accumulation while enhancing the proportion and diversity of unsaturated fatty acids, and providing conditions for the expansion and stability of membrane structure during mitochondrial biosynthesis through numerous gene regulations. Our findings provide a foundation for the identification of crucial functional genes, signaling pathways, and metabolic substances associated with adipose tissue differentiation and metabolism.

Effects of Adipose Tissue-Specific Knockout of Delta-like Non-Canonical Notch Ligand 1 on Lipid Metabolism in Mice

Delta-like non-canonical Notch ligand 1 (DLK1), which inhibits the differentiation of precursor adipocytes, is a recognized marker gene for precursor adipocytes. Lipids play a crucial role in energy storage and metabolism as a vital determinant of beef quality. In this study, we investigated the mechanism of the DLK1 gene in lipid metabolism by constructing adipose tissue-specific knockout mice. We examined some phenotypic traits, including body weight, liver coefficient, fat index, the content of triglyceride (TG) and cholesterol (CHOL) in abdominal white adipose tissue (WAT) and blood. Subsequently, the fatty acid content and genes related to lipid metabolism expression were detected in DLK1-/- and wild-type mice via GC-MS/MS analysis and quantitative real-time PCR (qRT-PCR), respectively. The results illustrated that DLK1-/- mice exhibited significant abdominal fat deposition compared to wild-type mice. HE staining and immunohistochemistry (IHC) results showed that the white adipocytes of DLK1-/- mice were larger, and the protein expression level of DLK1-/- was significantly lower. Regarding the blood biochemical parameters of female mice, DLK1-/- mice had a strikingly higher triglyceride content (p < 0.001). The fatty acid content in DLK1-/- mice was generally reduced. There was a significant reduction in the expression levels of the majority of genes that play a crucial role in lipid metabolism. This study reveals the molecular regulatory mechanism of fat metabolism in mice and provides a molecular basis and reference for the future application of the DLK1 gene in the breeding of beef cattle with an excellent meat quality traits. It also provides a molecular basis for unravelling the complex and subtle relationship between adipose tissue and health.

Comparative transcriptome and proteome analyses of the longissimus dorsi muscle for explaining the difference between donkey meat and other meats

Domestic donkeys (Equus asinus) have been maintained mainly for service purposes in the past. Nowadays, there is an increasing interest in donkey milk and meat production in several countries, including China. Donkey meat is highly consumed because of its nutritional value and unique flavor. However, genomic studies on donkey meat are limited. Therefore, in this study, we aimed to examine the molecular difference of longissimus dorsi muscles of donkey, cow, and goat. RNA sequencing and Proteome sequencing technology were used to analyze the transcriptome and proteome of the longissimus dorsi muscle of donkey, cow, and goat. A total of 1338 and 1780 differentially expressed genes (DEGs) were identified in donkey meat compared with that in cow and goat meat, respectively. Most of the DEGs were involved in biological processes, including small GTPase-mediated signal transduction, protein ubiquitination, protein glycosylation, and MAP kinase tyrosine/serine/threonine phosphatase activity. Additionally, 764 and 1024 differentially expressed proteins (DEPs) were identified in cow vs. donkey, and goat vs. donkey, respectively; these DEPs were mainly involved in metabolism. Genetic variation and regulatory factors can combine as a database to provide more valuable molecular information for further analysis.

Establishing the relations of characteristic aroma precursors and volatile compounds for authenticating Tibetan pork

Tibetan pork has been favored for its unique aromas, which originate from chemical reactions between characteristic precursors in cooking. The precursors (e.g., fatty acids, free amino acids, reducing sugars, and thiamine) of Tibetan pork ((semi-) free range) from different regions in China, comprising Tibet, Sichuan, Qinghai, and Yunnan, and commercial pork (indoor reared) were compared in this study. Tibetan pork was characterized by higher ω-3 polyunsaturated fatty acids (i.e., C18:3n3), higher essential (i.e., valine, leucine, and isoleucine), aromatic (i.e., phenylalanine), and sulfur-containing (i.e., methionine and cysteine) free amino acids, higher thiamine, and lower reducing sugars. Boiled Tibetan pork exhibited higher heptanal, 4-heptenal, and 4-pentylbenzaldehyde compared with commercial pork. The results from multivariate statistical analysis revealed that precursors combined with volatiles exhibited discriminating capability for characterizing Tibetan pork. The precursors in Tibetan pork exerted a certain effect on characteristic aroma generation, probably arising from promoting chemical reactions in cooking.

Integration of non-additive genome-wide association study with a multi-tissue transcriptome analysis of growth and carcass traits in Duroc pigs

Growth and carcass traits are of economic importance in the pig production, which affect pork quality and profitability of finishing pig production. This study used whole-genome and transcriptome sequencing technologies to identify potential candidate genes affecting growth and carcass traits in Duroc pigs. The medium (50-60 k) single nucleotide polymorphism (SNP) arrays of 4 154 Duroc pigs from three populations were imputed to whole-genome sequence data, yielding 10 463 227 markers on 18 autosomes. The dominance heritabilities estimated for growth and carcass traits ranged from 0.000 ± 0.041 to 0.161 ± 0.054. Using non-additive genome-wide association study (GWAS), we identified 80 dominance quantitative trait loci for growth and carcass traits at genome-wide significance (false discovery rate < 5%), 15 of which were also detected in our additive GWAS. After fine mapping, 31 candidate genes for dominance GWAS were annotated, and 8 of them were highlighted that have been previously reported to be associated with growth and development (e.g. SNX14, RELN and ENPP2), autosomal recessive diseases (e.g. AMPH, SNX14, RELN and CACNB4) and immune response (e.g. UNC93B1 and PPM1D). By integrating the lead SNPs with RNA-seq data of 34 pig tissues from the Pig Genotype-Tissue Expression project (https://piggtex.farmgtex.org/), we found that the rs691128548, rs333063869, and rs1110730611 have significantly dominant effects for the expression of SNX14, AMPH and UNC93B1 genes in tissues related to growth and development for pig, respectively. Finally, the identified candidate genes were significantly enriched for biological processes involved in the cell and organ development, lipids catabolic process and phosphatidylinositol 3-kinase signalling (P < 0.05). These results provide new molecular markers for meat production and quality selection of pig as well as basis for deciphering the genetic mechanisms of growth and carcass traits.

RNA-seq transcriptome profiling of pigs' liver in response to diet with different sources of fatty acids

Pigs (Sus scrofa) are an animal model for metabolic diseases in humans. Pork is an important source of fatty acids (FAs) in the human diet, as it is one of the most consumed meats worldwide. The effects of dietary inclusion of oils such as canola, fish, and soybean oils on pig gene expression are mostly unknown. Our objective was to evaluate FA composition, identify changes in gene expression in the liver of male pigs fed diets enriched with different FA profiles, and identify impacted metabolic pathways and gene networks to enlighten the biological mechanisms' variation. Large White male pigs were randomly allocated to one of three diets with 18 pigs in each; all diets comprised a base of corn and soybean meal to which either 3% of soybean oil (SOY), 3% canola oil (CO), or 3% fish oil (FO) was added for a 98-day trial during the growing and finishing phases. RNA sequencing was performed on the liver samples of each animal by Illumina technology for differential gene expression analyses, using the R package DESeq2. The diets modified the FA profile, mainly in relation to polyunsaturated and saturated FAs. Comparing SOY vs. FO, 143 differentially expressed genes (DEGs) were identified as being associated with metabolism, metabolic and neurodegenerative disease pathways, inflammatory processes, and immune response networks. Comparing CO vs. SOY, 148 DEGs were identified, with pathways related to FA oxidation, regulation of lipid metabolism, and metabolic and neurodegenerative diseases. Our results help explain the behavior of genes with differential expression in metabolic pathways resulting from feeding different types of oils in pig diets.

Transcriptomics and Lipid Metabolomics Analysis of Subcutaneous, Visceral, and Abdominal Adipose Tissues of Beef Cattle

Fat deposition traits are influenced by genetics and environment, which affect meat quality, growth rate, and energy metabolism of domestic animals. However, at present, the molecular mechanism of fat deposition is not entirely understood in beef cattle. Therefore, the current study conducted transcriptomics and lipid metabolomics analysis of subcutaneous, visceral, and abdominal adipose tissue (SAT, VAT, and AAT) of Huaxi cattle to investigate the differences among these adipose tissues and systematically explore how candidate genes interact with metabolites to affect fat deposition. These results demonstrated that compared with SAT, the gene expression patterns and metabolite contents of VAT and AAT were more consistent. Particularly, SCD expression, monounsaturated fatty acid (MUFA) and triglyceride (TG) content were higher in SAT, whereas PCK1 expression and the contents of saturated fatty acid (SFA), diacylglycerol (DG), and lysoglycerophosphocholine (LPC) were higher in VAT. Notably, in contrast to PCK1, 10 candidates including SCD, ELOVL6, ACACA, and FABP7 were identified to affect fat deposition through positively regulating MUFA and TG, and negatively regulating SFA, DG, and LPC. These findings uncovered novel gene resources and offered a theoretical basis for future investigation of fat deposition in beef cattle.

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.

Association of IFNA16 and TNFRSF19 Polymorphisms with Intramuscular Fat Content and Fatty Acid Composition in Pigs

Interferon-alpha-16 (IFNA16) and tumor necrosis factor receptor superfamily member 19 (TNFRSF19) are cytokines that may play a role in adipogenesis and fatness. Single nucleotide polymorphisms (SNPs) of the porcine IFNA16 and TNFRSF19 genes were verified and their association with intramuscular fat (IMF) content and fatty acid (FA) composition were evaluated in commercial crossbred pigs. Two non-synonymous SNPs of the porcine IFNA16 c.413G > A and TNFRSF19 c.860G > C loci were detected in commercial crossbred pigs. The porcine IFNA16 c.413G >A polymorphism was significantly associated with stearic acid, total saturated FAs (SFAs), and the ratio of monounsaturated FAs (MUFAs) to SFAs (p < 0.05). Furthermore, the porcine TNFRSF19 c.860G > C polymorphism was found to be significantly associated with IMF content and arachidic acid levels (p < 0.05). The results revealed that porcine IFNA16 and TNFRSF19 polymorphisms are related to IMF content and/or FA composition and affirmed the importance of these cytokine genes as potential candidate genes for lipid deposition and FA composition in the muscle tissue of pigs.

Genetic architecture and genomic selection of fatty acid composition predicted by Raman spectroscopy in rainbow trout

Background

In response to major challenges regarding the supply and sustainability of marine ingredients in aquafeeds, the aquaculture industry has made a large-scale shift toward plant-based substitutions for fish oil and fish meal. But, this also led to lower levels of healthful n−3 long-chain polyunsaturated fatty acids (PUFAs)—especially eicosapentaenoic (EPA) and docosahexaenoic (DHA) acids—in flesh. One potential solution is to select fish with better abilities to retain or synthesise PUFAs, to increase the efficiency of aquaculture and promote the production of healthier fish products. To this end, we aimed i) to estimate the genetic variability in fatty acid (FA) composition in visceral fat quantified by Raman spectroscopy, with respect to both individual FAs and groups under a feeding regime with limited n-3 PUFAs; ii) to study the genetic and phenotypic correlations between FAs and processing yields- and fat-related traits; iii) to detect QTLs associated with FA composition and identify candidate genes; and iv) to assess the efficiency of genomic selection compared to pedigree-based BLUP selection.

Results

Proportions of the various FAs in fish were indirectly estimated using Raman scattering spectroscopy. Fish were genotyped using the 57 K SNP Axiom™ Trout Genotyping Array. Following quality control, the final analysis contained 29,652 SNPs from 1382 fish. Heritability estimates for traits ranged from 0.03 ± 0.03 (n-3 PUFAs) to 0.24 ± 0.05 (n-6 PUFAs), confirming the potential for genomic selection. n-3 PUFAs are positively correlated to a decrease in fat deposition in the fillet and in the viscera but negatively correlated to body weight. This highlights the potential interest to combine selection on FA and against fat deposition to improve nutritional merit of aquaculture products. Several QTLs were identified for FA composition, containing multiple candidate genes with indirect links to FA metabolism. In particular, one region on Omy1 was associated with n-6 PUFAs, monounsaturated FAs, linoleic acid, and EPA, while a region on Omy7 had effects on n-6 PUFAs, EPA, and linoleic acid. When we compared the effectiveness of breeding programmes based on genomic selection (using a reference population of 1000 individuals related to selection candidates) or on pedigree-based selection, we found that the former yielded increases in selection accuracy of 12 to 120% depending on the FA trait.

Conclusion

This study reveals the polygenic genetic architecture for FA composition in rainbow trout and confirms that genomic selection has potential to improve EPA and DHA proportions in aquaculture species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08062-7.

Genome-wide analysis of expression QTL (eQTL) and allele-specific expression (ASE) in pig muscle identifies candidate genes for meat quality traits

BACKGROUND

Genetic analysis of gene expression level is a promising approach for characterizing candidate genes that are involved in complex economic traits such as meat quality. In the present study, we conducted expression quantitative trait loci (eQTL) and allele-specific expression (ASE) analyses based on RNA-sequencing (RNAseq) data from the longissimus muscle of 189 Duroc × Luchuan crossed pigs in order to identify some candidate genes for meat quality traits.

RESULTS

Using a genome-wide association study based on a mixed linear model, we identified 7192 cis-eQTL corresponding to 2098 cis-genes (p ≤ 1.33e-3, FDR ≤ 0.05) and 6400 trans-eQTL corresponding to 863 trans-genes (p ≤ 1.13e-6, FDR ≤ 0.05). ASE analysis using RNAseq SNPs identified 9815 significant ASE-SNPs in 2253 unique genes. Integrative analysis between the cis-eQTL and ASE target genes identified 540 common genes, including 33 genes with expression levels that were correlated with at least one meat quality trait. Among these 540 common genes, 63 have been reported previously as candidate genes for meat quality traits, such as PHKG1 (q-value = 1.67e-6 for the leading SNP in the cis-eQTL analysis), NUDT7 (q-value = 5.67e-13), FADS2 (q-value = 8.44e-5), and DGAT2 (q-value = 1.24e-3).

CONCLUSIONS

The present study confirmed several previously published candidate genes and identified some novel candidate genes for meat quality traits via eQTL and ASE analyses, which will be useful to prioritize candidate genes in further studies.

Identification of strong candidate genes for backfat and intramuscular fatty acid composition in three crosses based on the Iberian pig

Meat quality has an important genetic component and can be modified by the fatty acid (FA) composition and the amount of fat contained in adipose tissue and muscle. The present study aimed to find genomic regions associated with the FA composition in backfat and muscle (longissimus dorsi) in 439 pigs with three different genetic backgrounds but having the Iberian breed in common. Genome-wide association studies (GWAS) were performed between 38,424 single-nucleotide polymorphisms (SNPs) covering the pig genome and 60 phenotypic traits related to backfat and muscle FA composition. Nine significant associated regions were found in backfat on the Sus scrofa chromosomes (SSC): SSC1, SSC2, SSC4, SSC6, SSC8, SSC10, SSC12, and SSC16. For the intramuscular fat, six significant associated regions were identified on SSC4, SSC13, SSC14, and SSC17. A total of 52 candidate genes were proposed to explain the variation in backfat and muscle FA composition traits. GWAS were also reanalysed including SNPs on five candidate genes (ELOVL6, ELOVL7, FADS2, FASN, and SCD). Regions and molecular markers described in our study may be useful for meat quality selection of commercial pig breeds, although several polymorphisms were breed-specific, and further analysis would be needed to evaluate possible causal mutations.

About lipid metabolism in Hermetia illucens (L. 1758): on the origin of fatty acids in prepupae

Although increasingly targeted in animal nutrition, black soldier fly larvae or prepupae (BSF, Hermetia illucens L. 1758) require the characterization and modulation of their fatty acid profile to become fully integrated within the feed sector. This improvement will only be possible by the understanding of underlaying biochemical pathways of fatty acid synthesis in BSF. In this study, we hypothesized a labelling of de novo synthesized fatty acids in BSF by the incorporation of deuterated water (D₂O) in their feed. Three batches of fifty larvae were reared on two diets with different polyunsaturated fatty acid profiles moistened with 40% of H₂O or D₂O: chicken feed or 40% of chicken feed and 60% of flax cake. Although the occurrence of D₂O in insect feed increased the larval development time and decreased prepupal weight, it was possible to track the biosynthesis of fatty acids through deuterium labelling. Some fatty acids (decanoic, lauric or myristic acid) were exclusively present in their deuterated form while others (palmitic, palmitoleic or oleic acid) were found in two forms (deuterated or not) indicating that BSF can partially produce these fatty acids via biosynthesis pathways and not only by bioaccumulation from the diet. These results suggest the importance of carbohydrates as a source of acetyl-CoA in the constitution of the BSF fatty acid profile but also the potential importance of specific enzymes (e.g. thioesterase II or Δ12 fat2 desaturase) in BSF fatty acid metabolism. Finally, nearly no deuterated polyunsaturated fatty acids were found in BSF fed with deuterium confirming that BSF is not able to produce these types of fatty acids. Despite the high levels of linolenic acid in flax-enriched diets, BSF will simply bioaccumulate around 13% of this fatty acid and will metabolize approximately two-thirds of it into saturated fatty acids as lauric or myristic acid.

Genome-Wide Association Study of Body Weight Traits in Chinese Fine-Wool Sheep

Body weight is an important economic trait for sheep and it is vital for their successful production and breeding. Therefore, identifying the genomic regions and biological pathways that contribute to understanding variability in body weight traits is significant for selection purposes. In this study, the genome-wide associations of birth, weaning, yearling, and adult weights of 460 fine-wool sheep were determined using resequencing technology. The results showed that 113 single nucleotide polymorphisms (SNPs) reached the genome-wide significance levels for the four body weight traits and 30 genes were annotated effectively, including AADACL3, VGF, NPC1, and SERPINA12. The genes annotated by these SNPs significantly enriched 78 gene ontology terms and 25 signaling pathways, and were found to mainly participate in skeletal muscle development and lipid metabolism. These genes can be used as candidate genes for body weight in sheep, and provide useful information for the production and genomic selection of Chinese fine-wool sheep.

Effect of ELOVL6 on the lipid metabolism of bovine adipocytes

This study investigated the effect of ELOVL6 (elongation of very long chain fatty acids protein 6) and its underlying mechanism on lipid metabolism in bovine adipocytes. The ELOVL6 gene was overexpressed in bovine adipocytes by adenoviruses, and RNA sequencing was performed. Overexpression of ELOVL6 showed reduced proportions of C14:0 (Myristic) and C16:0 (palmitate) fatty acids and increased proportions of C18.0 (stearate) and C20:4n6 (arachidonic) fatty acids in adipocytes. In addition, a total of 2170 differentially expressed genes (DEGs) were found, containing 1802 up-regulated and 368 down-regulated genes. KEGG pathway analysis revealed that the down-regulated genes were linked with the regulation of lipolysis and the Wnt signaling pathway. The up-regulated genes were mainly involved in the FoxO signaling pathway; the PI3K-Akt signaling pathway; and the cAMP signaling pathway. In conclusion, our results suggest that ELOVL6 could affect the fatty acid composition in bovine adipocytes. We identified numerous related DEGs and pathways, which may provide a basis for studying the function and molecular mechanism of the ELOVL6 gene in regulating lipid metabolism.

Advances of Molecular Markers and Their Application for Body Variables and Carcass Traits in Qinchuan Cattle

This review considers the unique characteristics of Chinese cattle and intramuscular fat content (IMF) as factors influencing meat quality, including tenderness, flavor, and juiciness of meat. Due to its nutritional qualities, meat contributes to a healthy and balanced diet. The intramuscular fat content and eating quality of beef are influenced by many factors, which can generally be divided into on-farm and pre-slaughter factors (breed, sex of cattle, age at slaughter, housing system, diet, and pre-slaughter handling) and postmortem factors (post-slaughter processing, chilling temperature, and packaging). Meat quality traits can also be influenced by the individual genetic background of the animal. Worldwide, the function of genes and genetic polymorphisms that have potential effects on fattening of cattle and beef quality have been investigated. The use of DNA markers is recognized as a powerful and efficient approach to achieve genetic gain for desirable phenotypic characteristics, which is helpful for economic growth. The polymorphisms of the SIRT4, SIRT6, SIRT7, CRTC3, ABHD5, KLF6, H-FABP, and ELOVL6 genes for body and growth characteristics of cattle, and also for beef quality, are considered with the aim of highlighting the significance of beef intramuscular fat content, and that growth, body, and meat quality characteristics are polygenically regulated.

About the existence of common determinants of gene expression in the porcine liver and skeletal muscle

BACKGROUND

The comparison of expression QTL (eQTL) maps obtained in different tissues is an essential step to understand how gene expression is genetically regulated in a context-dependent manner. In the current work, we have compared the transcriptomic and eQTL profiles of two porcine tissues (skeletal muscle and liver) which typically show highly divergent expression profiles, in 103 Duroc pigs genotyped with the Porcine SNP60 BeadChip (Illumina) and with available microarray-based measurements of hepatic and muscle mRNA levels. Since structural variation could have effects on gene expression, we have also investigated the co-localization of cis-eQTLs with copy number variant regions (CNVR) segregating in this Duroc population.

RESULTS

The analysis of differential expresssion revealed the existence of 1204 and 1490 probes that were overexpressed and underexpressed in the gluteus medius muscle when compared to liver, respectively (|fold-change| > 1.5, q-value < 0.05). By performing genome scans in 103 Duroc pigs with available expression and genotypic data, we identified 76 and 28 genome-wide significant cis-eQTLs regulating gene expression in the gluteus medius muscle and liver, respectively. Twelve of these cis-eQTLs were shared by both tissues (i.e. 42.8% of the cis-eQTLs identified in the liver were replicated in the gluteus medius muscle). These results are consistent with previous studies performed in humans, where 50% of eQTLs were shared across tissues. Moreover, we have identified 41 CNVRs in a set of 350 pigs from the same Duroc population, which had been genotyped with the Porcine SNP60 BeadChip by using the PennCNV and GADA softwares, but only a small proportion of these CNVRs co-localized with the cis-eQTL signals.

CONCLUSION

Despite the fact that there are considerable differences in the gene expression patterns of the porcine liver and skeletal muscle, we have identified a substantial proportion of common cis-eQTLs regulating gene expression in both tissues. Several of these cis-eQTLs influence the mRNA levels of genes with important roles in meat (CTSF) and carcass quality (TAPT1), lipid metabolism (TMEM97) and obesity (MARC2), thus evidencing the practical importance of dissecting the genetic mechanisms involved in their expression.

Comparing the mRNA expression profile and the genetic determinism of intramuscular fat traits in the porcine gluteus medius and longissimus dorsi muscles

BACKGROUND

Intramuscular fat (IMF) content and composition have a strong impact on the nutritional and organoleptic properties of porcine meat. The goal of the current work was to compare the patterns of gene expression and the genetic determinism of IMF traits in the porcine gluteus medius (GM) and longissimus dorsi (LD) muscles.

RESULTS

A comparative analysis of the mRNA expression profiles of the pig GM and LD muscles in 16 Duroc pigs with available microarray mRNA expression measurements revealed the existence of 106 differentially expressed probes (fold-change > 1.5 and q-value < 0.05). Amongst the genes displaying the most significant differential expression, several loci belonging to the Hox transcription factor family were either upregulated (HOXA9, HOXA10, HOXB6, HOXB7 and TBX1) or downregulated (ARX) in the GM muscle. Differences in the expression of genes with key roles in carbohydrate and lipid metabolism (e.g. FABP3, ORMDL1 and SLC37A1) were also detected. By performing a GWAS for IMF content and composition traits recorded in the LD and GM muscles of 350 Duroc pigs, we identified the existence of one region on SSC14 (110-114 Mb) displaying significant associations with C18:0, C18:1(n-7), saturated and unsaturated fatty acid contents in both GM and LD muscles. Moreover, we detected several genome-wide significant associations that were not consistently found in both muscles. Further studies should be performed to confirm whether these associations are muscle-specific. Finally, the performance of an eQTL scan for 74 genes, located within GM QTL regions and with available microarray measurements of gene expression, made possible to identify 14 cis-eQTL regulating the expression of 14 loci, and six of them were confirmed by RNA-Seq.

CONCLUSIONS

We have detected significant differences in the mRNA expression patterns of the porcine LD and GM muscles, evidencing that the transcriptomic profile of the skeletal muscle tissue is affected by anatomical, metabolic and functional factors. A highly significant association with IMF composition on SSC14 was replicated in both muscles, highlighting the existence of a common genetic determinism, but we also observed the existence of a few associations whose magnitude and significance varied between LD and GM muscles.

A genome-wide detection of selection signatures in conserved and commercial pig breeds maintained in Poland

BACKGROUND

Identification of selection signatures can provide a direct insight into the mechanism of artificial selection and allow further disclosure of the candidate genes related to the animals' phenotypic variation. Domestication and subsequent long-time selection have resulted in extensive phenotypic changes in domestic pigs, involving a number of traits, like behavior, body composition, disease resistance, reproduction and coat color. In this study, based on genotypes obtained from PorcineSNP60 Illumina assay we attempt to detect both diversifying and within-breed selection signatures in 530 pigs belonging to four breeds: Polish Landrace, Puławska, Złotnicka White and Złotnicka Spotted, of which the last three are a subject of conservative breeding and substantially represent the native populations.

RESULTS

A two largely complementary statistical methods were used for signatures detection, including: pairwise F_ST and relative extended haplotype homozygosity (REHH) test. Breed-specific diversifying selection signals included several genes involved in processes connected with fertility, growth and metabolism which are potentially responsible for different phenotypes of the studied breeds. The diversifying selection signals also comprised PPARD gene that was previously found to have a large effect on the shape of the external ear in pigs or two genes encoding neuropeptide Y receptors (Y2 and Y5) involved in fat deposition and stress response which are important features differentiating the studied breeds. REHH statistics allowed detecting several within-breed selection signatures overlapping with genes connected with a range of functions including, among others: metabolic pathways, immune system response or implantation and development of the embryo.

CONCLUSIONS

The study provides many potential candidate genes with implication for traits selected in the individual breeds and gives strong basis for further studies aiming at identification of sources of variation among the studied pig breeds.

A polymorphism in the fatty acid desaturase-2 gene is associated with the arachidonic acid metabolism in pigs

Arachidonic acid (C20:4) is related to a wide range of biological effects including lipid homeostasis. The fatty acid desaturase-2 (FADS2) gene encodes for the delta-6-desaturase, which is involved in the biosynthesis of C20:4 from linoleic acid (C18:2). The purpose of this study was to characterise mutations in the promoter of the porcine FADS2, evaluating in particular the effect of one haplotype tagging polymorphism (rs321384923A > G) on the biosynthesis pathway of C20:4. A total of 1,192 Duroc barrows with records on fatty acid composition in muscle and subcutaneous fat were genotyped. Pigs carrying the A allele showed, irrespective of fat content, both enhanced FADS2 expression and higher C20:4 in muscle and exhibited increased ratios of C20:4 to C18:2 and of C20:4 to eicosadienoic acid (C20:2) in both muscle and adipose tissue. Despite the inverse relationship observed between C20:4 and fat content, the rs321384923 polymorphism had no impact on lean weight. It is concluded that the haplotype encompassing the rs321384923 polymorphism at the porcine FADS2 affects the n-6 fatty acid profile by specifically modifying the desaturation efficiency of C18:2 to C20:4 rather than by concomitant variations in C18:2 following changes in fat content.

A comparative analysis of label-free liquid chromatography-mass spectrometry liver proteomic profiles highlights metabolic differences between pig breeds

The liver is a complex organ governing several physiological processes that define biological mechanisms affecting growth, feed efficiency and performance traits in all livestock species, including pig. Proteomics may contribute to a better understanding of the relationship between liver functions and complex production traits in pigs and to characterize this species as biomedical model. This study applied, for the first time, a label‐free liquid chromatography-mass spectrometry (LC‐MS) proteomic approach to compare the liver proteome profiles of two important heavy pig breeds, Italian Duroc and Italian Large White. Liver specimens were collected (after slaughtering) from performance tested pigs of these two breeds, raised in standard conditions. The label‐free LC‐MS method captured a total of 501 proteins of which 200 were subsequently considered in the between breeds comparison. A statistical pipeline based on the sparse Partial Least Squares Discriminant Analysis (sPLS-DA), coupled with stability and significance tests, was applied for the identification of up or down regulated proteins between breeds. This analysis revealed a total of 25 proteins clearly separating Italian Duroc and Italian Large White pigs. Among the top proteins differentiating the two breeds, 3-ketoacyl-CoA thiolase, mitochondrial (ACAA2) and histone H2B type 2-F (HIST2H2BF) were up-regulated in Italian Duroc pigs and carboxylesterase 3 (CES3) and ketohexokinase (KHK) were up-regulated in Italian Large White pigs. Fatty acid synthase (FASN), involved in fatty acid metabolism and encoded by a gene located in a QTL region for fatty acid composition, was up-regulated in Italian Large White pigs. The in silico protein interaction analysis showed that 16 of these proteins were connected in one big module. Bioinformatic functional analysis indicated that differentially expressed proteins were involved in several biological processes related to the metabolism of lipids, amino-acids, carbohydrates, cofactors and antibiotics/drugs, suggesting that these functions might distinguish Italian Duroc and Italian Large White pigs. This pilot comparative proteomic analysis of the porcine liver highlighted several biological factors that could determine the peculiar production potentials of these two heavy pig breeds, derived by their different genetic backgrounds.

Expression analysis of candidate genes for fatty acid composition in adipose tissue and identification of regulatory regions

The aim of this work was to study the genetic basis of the backfat expression of lipid-related genes associated with meat quality traits in pigs. We performed a genome-wide association study with the backfat gene expression measured in 44 genes by qPCR and the PorcineSNP60 BeadChip genotypes in 115 Iberian x Landrace backcross animals. A total of 193 expression-associated SNPs located in 19 chromosomal regions were associated with expression levels of ACSM5, ELOVL6, FABP4, FADS2, and SLC27A4 genes. Three expression quantitative trail loci (eQTLs) corresponding to ACSM5, FABP4, and FADS2 were classified as cis-acting eQTLs, whereas the remaining 16 eQTLs have trans-regulatory effects. Remarkably, a SNP in the ACSM5 promoter region and a SNP in the 3′UTR region of FABP4 were the most associated polymorphisms with the ACSM5 and FABP4 expression levels, respectively. Moreover, relevant lipid-related genes mapped in the trans-eQTLs regions associated with the ACSM5, FABP4, FADS2, and SLC27A4 genes. Interestingly, a trans-eQTL hotspot on SSC13 regulating the gene expression of ELOVL6, ELOLV5, and SCD, three important genes implicated in the elongation and desaturation of fatty acids, was identified. These findings provide new data to further understand the functional regulatory mechanisms implicated in the variation of fatty acid composition in pigs.

TRIENNIAL GROWTH AND DEVELOPMENT SYMPOSIUM: Genetics and breeding for intramuscular fat and oleic acid content in pigs

The intramuscular fat (IMF) and oleic acid (OL) content have been favorably related to pork quality and human health. This influences the purchasing behavior of consumers and, therefore, also shifts the attention of breeding companies toward whether these traits are included into the breeding goal of the lines producing for high-valued markets. Because IMF and OL are unfavorably associated with lean content, a key economic trait, the real challenge for the industry is not simply to increase IMF and OL, but rather to come up with the right trade-off between them and lean content. In this paper we review the efforts performed to genetically improve IMF and OL, with particular reference to the research we conducted in a Duroc line aimed at producing high quality fresh and dry-cured pork products. Based on this research, we conclude that there are selection strategies that lead to response scenarios where IMF, OL, and lean content can be simultaneously improved. Such scenarios involve regular recording of IMF and OL, so that developing a cost-efficient phenotyping system for these traits is paramount. With the economic benefits of genomic selection needing further assessment in pigs, selection on a combination of pedigree-connected phenotypes and genotypes from a panel of selected genetic markers is presented as a suitable alternative. Evidence is provided supporting that at least a polymorphism in the leptin receptor and another in the stearoyl-CoA desaturase genes should be in that panel. Selection for IMF and OL results in an opportunity cost on lean growth. The extent to which it is affordable relies on the consumers' willingness to pay for premium products and on the cost to benefit ratio of alternative management strategies, such as specific dietary manipulations. How the genotype can influence the effect of the diet on IMF and OL remains a topic for further research.

Genome wide association study and genomic prediction for fatty acid composition in Chinese Simmental beef cattle using high density SNP array

BACKGROUND

Fatty acid composition of muscle is an important trait contributing to meat quality. Recently, genome-wide association study (GWAS) has been extensively used to explore the molecular mechanism underlying important traits in cattle. In this study, we performed GWAS using high density SNP array to analyze the association between SNPs and fatty acids and evaluated the accuracy of genomic prediction for fatty acids in Chinese Simmental cattle.

RESULTS

Using the BayesB method, we identified 35 and 7 regions in Chinese Simmental cattle that displayed significant associations with individual fatty acids and fatty acid groups, respectively. We further obtained several candidate genes which may be involved in fatty acid biosynthesis including elongation of very long chain fatty acids protein 5 (ELOVL5), fatty acid synthase (FASN), caspase 2 (CASP2) and thyroglobulin (TG). Specifically, we obtained strong evidence of association signals for one SNP located at 51.3 Mb for FASN using Genome-wide Rapid Association Mixed Model and Regression-Genomic Control (GRAMMAR-GC) approaches. Also, region-based association test identified multiple SNPs within FASN and ELOVL5 for C14:0. In addition, our result revealed that the effectiveness of genomic prediction for fatty acid composition using BayesB was slightly superior over GBLUP in Chinese Simmental cattle.

CONCLUSIONS

We identified several significantly associated regions and loci which can be considered as potential candidate markers for genomics-assisted breeding programs. Using multiple methods, our results revealed that FASN and ELOVL5 are associated with fatty acids with strong evidence. Our finding also suggested that it is feasible to perform genomic selection for fatty acids in Chinese Simmental cattle.

DNA methylation landscape of fat deposits and fatty acid composition in obese and lean pigs

Obese and lean type pig breeds exhibit differences in their fat deposits and fatty acid composition. Here, we compared the effect of genome-wide DNA methylation on fatty acid metabolism between Landrace pigs (LP, leaner) and Rongchang pigs (RP, fatty). We found that LP backfat (LBF) had a higher polyunsaturated fatty acid content but a lower adipocyte volume than RP backfat (RBF). LBF exhibited higher global DNA methylation levels at the genome level than RBF. A total of 483 differentially methylated regions (DMRs) were located in promoter regions, mainly affecting olfactory and sensory activity and lipid metabolism. In LBF, the promoters of genes related to ATPase activity had significantly stronger methylation. This fact may suggest lower energy metabolism levels, which may result in less efficient lipid synthesis in LBF. Furthermore, we identified a DMR in the miR-4335 and miR-378 promoters and validated their methylation status by bisulfite sequencing PCR. The hypermethylation of the promoters of miR-4335 and miR-378 in LBF and the resulting silencing of the target genes may result in LBF's low content in saturated fatty acids and fat deposition capacity. This study provides a solid basis for exploring the epigenetic mechanisms affecting fat deposition and fatty acid composition.