Polymorphism in the ELOVL6 gene is associated with a major QTL effect on fatty acid composition in pigs

Whole genome sequencing revealed genetic structure, domestication, and selection of Chinese indigenous ducks

The genetic architecture and domestication history of Chinese indigenous ducks, particularly those with distinct traits like the Longsheng duck (LSD), are not well characterized. This study used whole-genome resequencing data from 540 ducks across 30 populations to explore the genetic diversity and structure of these ducks. Our findings suggest extensive interspecific hybridization between mallard and spot-billed duck. Comparing with other Chinese ducks, LSD is a distinct breed with moderate genetic diversity. Selective signal analysis identified several key genes impacting neural development, fat metabolism, immunity, and circadian rhythms like SLC25A20 and PER2. These genes showed strong selective pressures that parallelled other domestication processes. Additionally, EDNRB2 was identified as a potential gene influencing the unique coloration of LSD's plumage, bill, and webbed feet, associated with distinct mutations in non-coding regions. Comparative analysis with other indigenous breeds further pinpointed genes associated with LSD-specific traits, including plumage color, reproductive capabilities, and fat deposition, such as MITF, SPATA2, EIF2S2, PLIN3, ATP1B1, and CCDC80. Our findings clarify the population genetics of Chinese indigenous ducks. They also highlight key genes and mutations that shape the unique characteristics of LSD. These insights pave the way for further research into the genetic resources and domestication patterns of Chinese ducks.

Whole-genome sequencing reveals patterns of runs of homozygosity underlying genetic diversity and selection in domestic rabbits

BACKGROUND

Runs of homozygosity (ROH) are continuous segments of homozygous genotypes inherited from both parental lineages. These segments arise due to the transmission of identical haplotypes. The genome-wide patterns and hotspot regions of ROH provide valuable insights into genetic diversity, demographic history, and selection trends. In this study, we analyzed whole-genome resequencing data from 117 rabbits to identify ROH patterns and inbreeding level across eleven rabbit breeds, including seven Chinese indigenous breeds and four exotic breeds, and to uncover selective signatures based on ROH islands.

RESULTS

We detected a total of 31,429 ROHs across the autosomes of all breeds, with the number of ROHs (NROH) per breed ranging from 1316 to 7476. The mean sum of ROHs length (SROH) per individual was 493.84 Mb, covering approximately 22.79% of the rabbit autosomal genome. The majority of the detected ROHs ranged from 1 to 2 Mb in length, with an average ROH length (LROH) of 1.84 Mb. ROHs longer than 6 Mb constituted only 0.83% of the detected ROHs. The average inbreeding coefficient derived from ROHs (FROH) was 0.23, with FROH values ranging from 0.14 to 0.38 across breeds. Among Chinese indigenous breeds, the Jiuyishan rabbit exhibited the highest values of NROH, SROH, LROH, and FROH, whereas the Fujian Yellow rabbit had the lowest FROH values. In exotic rabbit breeds, the Japanese White rabbit displayed the highest values for NROH, SROH, LROH, and FROH, while the Flemish Giant rabbit had the lowest values for these metrics. Additionally, we identified 17 ROH islands in Chinese indigenous breeds and 22 ROH islands in exotic rabbit breeds, encompassing 124 and 186 genes, respectively. In Chinese indigenous breeds, we identified prominent genes associated with reproduction, including CFAP206, RNF133, CPNE4, ASTE1, and ATP2C1, as well as genes related to adaptation, namely CADPS2, FEZF1, and EPHA7. In contrast, the exotic breeds exhibited a prevalence of genes associated with fat deposition, such as ELOVL3 and NPM3, as well as growth and body weight related genes, including FAM184B, NSMCE2, and TWNK.

CONCLUSIONS

This study enhances our understanding of genetic diversity and selection pressures in domestic rabbits, offering valuable implications for breeding management and conservation strategies.

Gene co-expression network analysis for porcine intramuscular fatty acid composition

In pigs, meat quality depends markedly on the fatty acid (FA) content and composition of the intramuscular fat, which is partly determined by the gene expression in this tissue. The aim of this work was to identify the link between muscle gene expression and its FA composition. In an (Iberian × Duroc) × Duroc backcrossed pig population, we identified modules of co-expressed genes, and correlation analyses were performed for each of them versus the phenotypes, finding four relevant modules. Two of the modules were positively correlated with saturated FAs (SFAs) and monounsaturated FAs (MUFAs), while negatively correlated with polyunsaturated FAs (PUFAs) and the omega-6/omega-3 ratio. The gene-enrichment analysis showed that these modules had over-representation of pathways related with the biosynthesis of unsaturated FAs, the Peroxisome proliferator-activated receptor signalling pathway and FA elongation. The two other relevant modules were positively correlated with PUFA and the n-6/n-3 ratio, but negatively correlated with SFA and MUFA. In this case, they had an over-representation of pathways related with fatty and amino acid degradation, and with oxidative phosphorylation. Using a graphical Gaussian model, we inferred a network of connections between the genes within each module. The first module had 52 genes with 87 connections, and the most connected genes were ADIPOQ, which is related with FA oxidation, and ELOVL6 and FABP4, both involved in FA metabolism. The second module showed 196 genes connected by 263 edges, being FN1 and MAP3K11 the most connected genes. On the other hand, the third module had 161 genes connected by 251 edges and ATG13 was the top neighbouring gene, while the fourth module had 224 genes and 655 connections, and its most connected genes were related with mitochondrial pathways. Overall, this work successfully identified relevant muscle gene networks and modules linked with FA composition, providing further insights on how the physiology of the pigs influences FA composition.

Identification of candidate regulatory genes for intramuscular fatty acid composition in pigs by transcriptome analysis

BACKGROUND

Intramuscular fat (IMF) content and its fatty acid (FA) composition are typically controlled by several genes, each with a small effect. In the current study, to pinpoint candidate genes and putative regulators involved in FA composition, we performed a multivariate integrative analysis between intramuscular FA and transcriptome profiles of porcine longissimus dorsi (LD) muscle. We also carried out a combination of network, regulatory impact factor (RIF), in silico prediction of putative target genes, and functional analyses to better support the biological relevance of our findings.

RESULTS

For this purpose, we used LD RNA-Seq and intramuscular FA composition profiles of 129 Iberian × Duroc backcrossed pigs. We identified 378 correlated variables (13 FA and 365 genes), including six FA (C20:4n-6, C18:2n-6, C20:3n-6, C18:1n-9, C18:0, and C16:1n-7) that were among the most interconnected variables in the predicted network. The detected FA-correlated genes include genes involved in lipid and/or carbohydrate metabolism or in regulation of IMF deposition (e.g., ADIPOQ, CHUK, CYCS, CYP4B1, DLD, ELOVL6, FBP1, G0S2, GCLC, HMGCR, IDH3A, LEP, LGALS12, LPIN1, PLIN1, PNPLA8, PPP1R1B, SDR16C5, SFRP5, SOD3, SNW1, and TFRC), meat quality (GALNT15, GOT1, MDH1, NEU3, PDHA1, SDHD, and UNC93A), and transport (e.g., EXOC7 and SLC44A2). Functional analysis highlighted 54 over-represented gene ontology terms, including well-known biological processes and pathways that regulate lipid and carbohydrate metabolism. RIF analysis suggested a pivotal role for six transcription factors (CARHSP1, LBX1, MAFA, PAX7, SIX5, and TADA2A) as putative regulators of gene expression and intramuscular FA composition. Based on in silico prediction, we identified putative target genes for these six regulators. Among these, TADA2A and CARHSP1 had extreme RIF scores and present novel regulators in pigs. In addition, the expression of TADA2A correlated (either positively or negatively) with C20:4n-6, C18:2n-6, C20:3n-6, C18:1n-9, and that of CARHSP1 correlated (positively) with the C16:1n-7 lipokine. We also found that these two transcription factors share target genes that are involved in lipid metabolism (e.g., GOT1, PLIN1, and TFRC).

CONCLUSIONS

This integrative analysis of muscle transcriptome and intramuscular FA profile revealed valuable information about key candidate genes and potential regulators for FA and lipid metabolism in pigs, among which some transcription factors are proposed to control gene expression and modulate FA composition differences.

Genome-wide association study for stayability at different calvings in Nellore beef cattle

BACKGROUNDING

Stayability, which may be defined as the probability of a cow remaining in the herd until a reference age or at a specific number of calvings, is usually measured late in the animal's life. Thus, if used as selection criteria, it will increase the generation interval and consequently might decrease the annual genetic gain. Measuring stayability at an earlier age could be a reasonable strategy to avoid this problem. In this sense, a better understanding of the genetic architecture of this trait at different ages and/or at different calvings is important. This study was conducted to identify possible regions with major effects on stayability measured considering different numbers of calvings in Nellore cattle as well as pathways that can be involved in its expression throughout the female's productive life.

RESULTS

The top 10 most important SNP windows explained, on average, 17.60% of the genetic additive variance for stayability, varying between 13.70% (at the eighth calving) and 21% (at the fifth calving). These SNP windows were located on 17 chromosomes (1, 2, 4, 6, 7, 8, 9, 10, 11, 12, 13, 14, 18, 19, 20, 27, and 28), and they harbored a total of 176 annotated genes. The functional analyses of these genes, in general, indicate that the expression of stayability from the second to the sixth calving is mainly affected by genetic factors related to reproductive performance, and nervous and immune systems. At the seventh and eighth calvings, genes and pathways related to animal health, such as density bone and cancer, might be more relevant.

CONCLUSION

Our results indicate that part of the target genomic regions in selecting for stayability at earlier ages (from the 2th to the 6th calving) would be different than selecting for this trait at later ages (7th and 8th calvings). While the expression of stayability at earlier ages appeared to be more influenced by genetic factors linked to reproductive performance together with an overall health/immunity, at later ages genetic factors related to an overall animal health gain relevance. These results support that selecting for stayability at earlier ages (perhaps at the second calving) could be applied, having practical implications in breeding programs since it could drastically reduce the generation interval, accelerating the genetic progress.

The Identification of Functional Genes Affecting Fat-Related Meat Traits in Meat-Type Pigeons Using Double-Digest Restriction-Associated DNA Sequencing and Molecular Docking Analysis

The Chinese indigenous Shiqi (SQ) pigeon and the imported White King (WK) pigeon are two meat-type pigeon breeds of economical and nutritional importance in China. They displayed significant differences in such meat quality traits as intramuscular fat (IMF) content and fatty acid (FA) compositions in the breast muscles. In this study, we aimed to screen candidate genes that could affect fat-related meat quality traits in meat-type pigeons. We investigated the polymorphic variations at the genomic level using double-digest restriction-associated DNA (ddRAD) sequencing in 12 squabs of SQ and WK pigeons that exhibited significant inter-breed differences in IMF content as well as FA and amino acid compositions in the breast muscles, and screened candidate genes influencing fat-related traits in squabs through gene ontology analysis and pathway analysis. By focusing on 6019 SNPs, which were located in genes with correct annotations and had the potential to induce changes in the encoded proteins, we identified 19 genes (ACAA1, ACAA2, ACACB, ACADS, ACAT1, ACOX3, ACSBG1, ACSBG2, ACSL1, ACSL4, ELOVL6, FADS1, FADS2, HACD4, HADH, HADHA, HADHB, MECR, OXSM) as candidate genes that could affect fat-related traits in squabs. They were significantly enriched in the pathways of FA metabolism, degradation, and biosynthesis (p < 0.05). Results from molecular docking analysis further revealed that three non-synonymous amino acid alterations, ACAA1(S357N), ACAA2(T234I), and ACACB(H1418N), could alter the non-bonding interactions between the enzymatic proteins and their substrates. Since ACAA1, ACAA2, and ACACB encode rate-limiting enzymes in FA synthesis and degradation, alterations in the enzyme-substrate binding affinity may subsequently affect the catalytic efficiency of enzymes. We suggested that SNPs in these three genes were worthy of further investigation into their roles in explaining the disparities in fat-related traits in squabs.

Genome-wide association study reveals candidate genes for traits related to meat quality in Colombian Creole hair sheep

Genome-wide association studies (GWAS) allow identifying genomic regions related to traits of economic importance in animals of zootechnical interest. The objective of this research was to conduct a genome-wide association study on meat quality traits using the Illumina OvineSNPs50 BeadChip array. The animals were sampled in the departments of Córdoba, Cesar, and Valle del Cauca. The genotypes obtained with the Illumina OvineSNP50 BeadChip microarray were analyzed SNP (single-nucleotide polymorphism) data to conduct a GWAS for pH and water-holding capacity (WHC) traits measured after 7 days of maturation, in the Longissimus dorsi (LD) muscle, in 167 Creole hair sheep of 12 months old belonging to Pelibuey (CHSP, n = 60), Ethiopian (CHSE, n = 44), and Sudan (CHSS, n = 63) breeds. The GWAS was done using a mixed linear model (MLMA) and based on the Ovis aries v3.1 genome. The CHSE showed the lowest meat juice release and, consequently, the highest water-holding capacity (WHC = 30.6 ± 0.1), suggesting that this breed has better performance in the meat industry compared with CHSS (WHC = 41.7 ± 0.1) and CHSP (WHC = 36.8 ± 0.1), since there is a relationship between WHC and juiciness. For the character pH, it was not possible to annotate genes related to meat quality, while, for the WHC, they have obtained 11 candidate genes associated (ELOVL2, ARAP2, LOC101102527, SHOC2, AIPL1, CSRNP3, IFRD, KDM8, NANS, DAPK1, IBN2, TPM2). Particularly, ELOVL2, ARAP2, IBN2, and TPM2 genes are involved in muscle contraction and fatty acid composition in sheep. In this study, we generated a baseline for GWAS related to meat quality traits in Colombian Creole hair sheep that can be used for future genomic selection plans.

Transcriptomic Profiling of Subcutaneous Backfat in Castrated and Intact Alentejano Pigs Finished Outdoors with Commercial and Fiber-Rich Diets

In this work, we studied the backfat transcriptome of surgically castrated (C), intact (I) and intact fed an experimental diet (IE) outdoor-reared male Alentejano (AL) pigs. The experimental diet was a high-fiber diet with locally produced legumes and by-products associated with a boar taint reduction effect. At slaughter (~160 kg), backfat samples were collected for total RNA sequencing. Intact pigs presented leaner carcasses, more total collagen, and more unsaturated intramuscular fat content than C animals. A total of 2726 differentially expressed genes (DEGs, |log2 FC|> 0.58, q < 0.05) were identified between C and I with overexpressed genes related to muscular activity (MYH1, ACTA1) or collagen metabolism (COL1A1, COL1A2) in I pigs. Between C and IE, 1639 DEGs of genes involved in lipidic metabolism (LEP, ME1, FABP4, ELOVL6) were overexpressed in C. Finally, only 28 DEGs were determined between I and IE. Clustering results indicated a drastic influence of the testis in the transcriptome of subcutaneous fat of AL pigs, while the diet had a marginal effect. Diet can reduce stress by increasing satiety in animals, and could have induced an increase of skatole degradation due to the higher expression of the CYP2A19 gene in the IE group.

Chromosome-level genome and population genomics provide novel insights into adaptive divergence in allopatric Eleutheronema tetradactylum

Understanding the adaptive ecological divergence provides important information for revealing biodiversity generation and maintenance. Adaptive ecology divergence in populations occurs in various environments and locations, but its genetic underpinnings remain elusive. We generated a chromosome-level genome of Eleutheronema tetradactylum (~582 Mb) and re-sequenced 50 allopatric E. tetradactylum in two independent environmental axes in China and Thailand Coastal waters as well as 11 cultured relatives. A low level of whole genome-wide diversity explained their decreased adaptive potential in the wild environment. Demographic analysis showed evidence of historically high abundance followed by a continuous distinct decline, plus signs of recent inbreeding and accumulation of deleterious mutations. Extensive signals of selective sweeps with signs of local adaptation to environmental differentiation between China and Thailand at genes related to thermal and salinity adaptation were discovered, which might be the driving factors of the geographical divergence of E. tetradactylum. Many genes and pathways subjected to strong selection under artificial breeding were associated with fatty acids and immunity (ELOVL6L, MAPK, p53/NF-kB), likely contributing to the eventual adaptation of artificial selective breeding. Our comprehensive study provided crucial genetic information for E. tetradactylum, with implications for the further conservation efforts of this threatened and ecologically valuable fish.

Transcriptome analysis of perirenal fat from Spanish Assaf suckling lamb carcasses showing different levels of kidney knob and channel fat

Introduction

Suckling lamb meat is highly appreciated in European Mediterranean countries because of its mild flavor and soft texture. In suckling lamb carcasses, perirenal and pelvic fat depots account for a large fraction of carcass fat accumulation, and their proportions are used as an indicator of carcass quality.

Material and Methods

This study aimed to characterize the genetic mechanisms that regulate fat deposition in suckling lambs by evaluating the transcriptomic differences between Spanish Assaf lambs with significantly different proportions of kidney knob and channel fat (KKCF) depots in their carcasses (4 High-KKCF lambs vs. 4 Low-KKCF lambs).

Results

The analyzed fat tissue showed overall dominant expression of white adipose tissue gene markers, although due to the young age of the animals (17-36 days), the expression of some brown adipose tissue gene markers (e.g., UCP1, CIDEA) was still identified. The transcriptomic comparison between the High-KKCF and Low-KKCF groups revealed a total of 80 differentially expressed genes (DEGs). The enrichment analysis of the 49 DEGs with increased expression levels in the Low-KKCF lambs identified significant terms linked to the biosynthesis of lipids and thermogenesis, which may be related to the higher expression of the UCP1 gene in this group. In contrast, the enrichment analysis of the 31 DEGs with increased expression in the High-KKCF lambs highlighted angiogenesis as a key biological process supported by the higher expression of some genes, such as VEGF-A and THBS1, which encode a major angiogenic factor and a large adhesive extracellular matrix glycoprotein, respectively.

Discussion

The increased expression of sestrins, which are negative regulators of the mTOR complex, suggests that the preadipocyte differentiation stage is being inhibited in the High-KKCF group in favor of adipose tissue expansion, in which vasculogenesis is an essential process. All of these results suggest that the fat depots of the High-KKCF animals are in a later stage of development than those of the Low-KKCF lambs. Further genomic studies based on larger sample sizes and complementary analyses, such as the identification of polymorphisms in the DEGs, should be designed to confirm these results and achieve a deeper understanding of the genetic mechanisms underlying fat deposition in suckling lambs.

Global analysis of the association between pig muscle fatty acid composition and gene expression using RNA-Seq

Fatty acids (FAs) play an essential role as mediators of cell signaling and signal transduction, affecting metabolic homeostasis and determining meat quality in pigs. However, FAs are transformed by the action of several genes, such as those encoding desaturases and elongases of FAs in lipogenic tissues. The aim of the current work was to identify candidate genes, biological processes, and pathways involved in the modulation of intramuscular FA profile from longissimus dorsi muscle. FA profile by gas chromatography of methyl esters and gene expression by RNA-Seq were determined in 129 Iberian × Duroc backcrossed pigs. An association analysis between the muscle transcriptome and its FA profile was performed, followed by a concordance and functional analysis. Overall, a list of well-known (e.g., PLIN1, LEP, ELOVL6, SC5D, NCOA2, ACSL1, MDH1, LPL, LGALS12, TFRC, GOT1, and FBP1) and novel (e.g., TRARG1, TANK, ENSSSCG00000011196, and ENSSSCG00000038429) candidate genes was identified, either in association with specific or several FA traits. Likewise, several of these genes belong to biological processes and pathways linked to energy, lipid, and carbohydrate metabolism, which seem determinants in the modulation of FA compositions. This study can contribute to elucidate the complex relationship between gene expression and FA profile in pig muscle.

ELOVL gene family plays a virtual role in response to breeding selection and lipid deposition in different tissues in chicken (Gallus gallus)

Background

Elongases of very long chain fatty acids (ELOVLs), a family of first rate-limiting enzymes in the synthesis of long-chain fatty acids, play an essential role in the biosynthesis of complex lipids. Disrupting any of ELOVLs affects normal growth and development in mammals. Genetic variations in ELOVLs are associated with backfat or intramuscular fatty acid composition in livestock. However, the effects of ELOVL gene family on breeding selection and lipid deposition in different tissues are still unknown in chickens.

Results

Genetic variation patterns and genetic associations analysis showed that the genetic variations of ELOVL genes were contributed to breeding selection of commercial varieties in chicken, and 14 SNPs in ELOVL2-6 were associated with body weight, carcass or fat deposition traits. Especially, one SNP rs17631638T > C in the promoter of ELOVL3 was associated with intramuscular fat content (IMF), and its allele frequency was significantly higher in native and layer breeds compared to that in commercial broiler breeds. Quantitative real-time PCR (qRT-PCR) determined that the ELOVL3 expressions in pectoralis were affected by the genotypes of rs17631638T > C. In addition, the transcription levels of ELOVL genes except ELOVL5 were regulated by estrogen in chicken liver and hypothalamus with different regulatory pathways. The expression levels of ELOVL1-6 in hypothalamus, liver, abdominal fat and pectoralis were correlated with abdominal fat weight, abdominal fat percentage, liver lipid content and IMF. Noteworthily, expression of ELOVL3 in pectoralis was highly positively correlated with IMF and glycerophospholipid molecules, including phosphatidyl choline, phosphatidyl ethanolamine, phosphatidyl glycerol and phospholipids inositol, rich in ω-3 and ω-6 long-chain unsaturated fatty acids, suggesting ELOVL3 could contribute to intramuscular fat deposition by increasing the proportion of long-chain unsaturated glycerophospholipid molecules in pectoralis.

Conclusions

In summary, we demonstrated the genetic contribution of ELOVL gene family to breeding selection for specialized varieties, and revealed the expression regulation of ELOVL genes and their potential roles in regulating lipid deposition in different tissues. This study provides new insights into understanding the functions of ELOVL family on avian growth and lipid deposition in different tissues and the genetic variation in ELOVL3 may aid the marker-assisted selection of meat quality in chicken.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08932-8.

Early fish domestication affects methylation of key genes involved in the rapid onset of the farmed phenotype

Animal domestication is a process of environmental modulation and artificial selection leading to permanent phenotypic modifications. Recent studies showed that phenotypic changes occur very early in domestication, i.e., within the first generation in captivity, which raises the hypothesis that epigenetic mechanisms may play a critical role on the early onset of the domestic phenotype. In this context, we applied reduced representation bisulphite sequencing to compare methylation profiles between wild Nile tilapia females and their offspring reared under farmed conditions. Approximately 700 differentially methylated CpG sites were found, many of them associated not only with genes involved in muscle growth, immunity, autophagy and diet response but also related to epigenetic mechanisms, such as RNA methylation and histone modifications. This bottom-up approach showed that the phenotypic traits often related to domestic animals (e.g., higher growth rate and different immune status) may be regulated epigenetically and prior to artificial selection on gene sequences. Moreover, it revealed the importance of diet in this process, as reflected by differential methylation patterns in genes critical to fat metabolism. Finally, our study highlighted that the TGF-β1 signalling pathway may regulate and be regulated by several differentially methylated CpG-associated genes. This could be an important and multifunctional component in promoting adaptation of fish to a domestic environment while modulating growth and immunity-related traits.

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.

Transcriptome study digs out BMP2 involved in adipogenesis in sheep tails

Background

Hu sheep and Tibetan sheep in China are characterized by fat tails and thin tails, respectively. Several transcriptomes have been conducted in different sheep breeds to identify the differentially expressed genes (DEGs) underlying this trait. However, these studies identified different DEGs in different sheep breeds.

Results

Hence, RNA sequencing was performed on Hu sheep and Tibetan sheep. We obtained a total of 45.57 and 43.82 million sequencing reads, respectively. Two libraries mapped reads from 36.93 and 38.55 million reads after alignment to the reference sequences. 2108 DEGs were identified, including 1247 downregulated and 861 upregulated DEGs. GO and KEGG analyses of all DEGs demonstrated that pathways were enriched in the regulation of lipolysis in adipocytes and terms related to the chemokine signalling pathway, lysosomes, and glycosaminoglycan degradation. Eight genes were selected for validation by RT–qPCR. In addition, the transfection of BMP2 overexpression into preadipocytes resulted in increased PPAR-γ expression and expression. BMP2 potentially induces adipogenesis through LOX in preadipocytes. The number of lipid drops in BMP2 overexpression detected by oil red O staining was also greater than that in the negative control.

Conclusion

In summary, these results showed that significant genes (BMP2, HOXA11, PPP1CC and LPIN1) are involved in the regulation of adipogenesis metabolism and suggested novel insights into metabolic molecules in sheep fat tails.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08657-8.

Development and Validation of KASP Assays for the Genotyping of Racing Performance-Associated Single Nucleotide Polymorphisms in Pigeons

Pigeon racing's recent upturn in popularity can be attributed in part to the huge prize money involved in these competitions. As such, methods to select pigeons with desirable genetic characteristics for racing or for selective breeding have also been gaining more interest. Polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) for genotyping-specific genes is one of the most commonly used molecular techniques, which can be costly, laborious and time consuming. The present study reports the development of an alternative genotyping method that employs Kompetitive Allele Specific Polymerase Chain Reaction (KASP) technology with specifically designed primers to detect previously reported racing performance-associated polymorphisms within the LDHA, MTYCB, and DRD4 genes. To validate, KASP assays and PCR-RFLP assays results from 107 samples genotyped for each of the genes were compared and the results showed perfect (100%) agreement of both methods. The developed KASP assays present an alternative rapid, reliable, and cost-effective method to identify polymorphisms in pigeons.

GWAS and Post-GWAS High-Resolution Mapping Analyses Identify Strong Novel Candidate Genes Influencing the Fatty Acid Composition of the Longissimus dorsi Muscle in Pigs

Fatty acid (FA) composition is one of the most important parameters for the assessment of meat quality in pigs. The FA composition in pork can also affect human health. Our aim was to identify quantitative trait loci (QTLs) and positional candidate genes affecting the FA profile of the longissimus dorsi muscle in a large F2 intercross between Landrace and Korean native pigs comprising 1105 F2 progeny by genome-wide association studies (GWAS) and post-GWAS high-resolution mapping analyses. We performed GWAS using the PorcineSNP60K BeadChip and a linear mixed model. Four genome-wide significant QTL regions in SSC8, SSC12, SSC14, and SSC16 were detected (p < 2.53 × 10-7). Several co-localizations of QTLs in SSC12 for oleic acid, linoleic acid, arachidonic acid, monounsaturated FAs, polyunsaturated FAs, and the polyunsaturated/saturated FA ratio were observed. To refine the QTL region in SSC12, a linkage and linkage disequilibrium analysis was applied and could narrow down the critical region to a 0.749 Mb region. Of the genes in this region, GAS7, MYH2, and MYH3 were identified as strong novel candidate genes based on further conditional association analyses. These findings provide a novel insight into the genetic basis of FA composition in pork and could contribute to the improvement of pork quality.

Transcriptomic diversity in longissimus thoracis muscles of Barbari and Changthangi goat breeds of India

The present study is an attempt to examine the differential expression of genes in longissimus thoracis muscles between meat and wool type Indian goat breeds. Barbari goat is considered the best meat breed while Changthangi is famous for its fine fibre quality. RNA sequencing data was generated from four biological replicates of longissimus thoracis muscles of Barbari and Changthangi goats. A clear demarcation could be observed between the breeds in terms of expression of genes associated with lipid metabolism (FASN, SCD, THRSP, DGAT2 and FABP3). Most significant genes with high connectivity identified by gene co-expression network analysis were associated with triacylglycerol biosynthesis pathway in Barbari goat. Highly interactive genes identified in Changthangi goat were mainly associated with muscle fibre type. This study provides an insight into the differential expression of genes in longissimus thoracis muscles between Barbari and Changthangi goats that are adapted to and reared in different agro-climatic regions.

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.

RNA-Seq based genetic variant discovery provides new insights into controlling fat deposition in the tail of sheep

Genetic basis of fat deposition in sheep tail have not been completely elucidated yet. Understanding the genetic mechanisms controlling fat-tail size can improve breeding strategies to modulate fat deposition. RNA sequencing has made it possible to discover genetic variants that may underlie various phenotypic differences. Hence, to identify genetic variants that are important for describing different fat-tail phenotypes in sheep, RNA sequencing was used for single nucleotide polymorphism (SNP) calling in two Iranian sheep breeds (Lori-Bakhtiari, fat-tailed; n = 4, vs Zel, thin-tailed; n = 4). Using a stringent pipeline, a total of 112,344 known SNPs were genotyped, of which 30,550 and 42,906 SNPs were shared by at least two Lori-Bakhtiari and Zel, respectively. Comparing these SNPs showed 2,774 (including 209 missense and 25 deleterious SNPs) and 10,470 (including 1,054 missense and 116 deleterious SNPs) breed-specific SNPs in Lori-Bakhtiari and Zel sheep, respectively. Potential breed-specific SNPs were detected by considering those located in QTL regions associated with fatness or reported as important candidates in previous similar studies. Of the breed-specific SNPs, 724 and 2,905 were located in the QTL regions. Functional enrichment analysis of the affected genes revealed several enriched gene ontologies and KEGG pathways related to fat metabolism. Based on the results, several affected genes were proposed to be strongly linked with fat deposition such as DGAT2, ACSL1, ACACA, ADIPOQ, ACLY, FASN, CPT2, SCD, ADCY6, PER3, CSF1R, SLC22A4, GFPT1, CDS2, BMP6, ACSS2, ELOVL6, HOXA10 and FABP4. Moreover, several SNPs were found in the candidate genes related to fatty acid oxidation introducing them as promising candidates responsible for lower fat content in tail of Zel. Our findings provided new insights into the genetic mechanisms of fat deposition in sheep, which can serve to designing appropriate breeding programs.

Identification of eQTLs associated with lipid metabolism in Longissimus dorsi muscle of pigs with different genetic backgrounds

Intramuscular fat content and its fatty acid composition affect porcine meat quality and its nutritional value. The present work aimed to identify genomic variants regulating the expression in the porcine muscle (Longissimus dorsi) of 45 candidate genes for lipid metabolism and fatty acid composition in three experimental backcrosses based on the Iberian breed. Expression genome-wide association studies (eGWAS) were performed between the muscle gene expression values, measured by real-time quantitative PCR, and the genotypes of 38,426 SNPs distributed along all chromosomes. The eGWAS identified 186 eSNPs located in ten Sus scrofa regions and associated with the expression of ACSM5, ACSS2, ATF3, DGAT2, FOS and IGF2 (FDR < 0.05) genes. Two expression quantitative trait loci (eQTLs) for IGF2 and ACSM5 were classified as cis-acting eQTLs, suggesting a mutation in the same gene affecting its expression. Conversely, ten eQTLs showed trans-regulatory effects on gene expression. When the eGWAS was performed for each backcross independently, only three common trans-eQTL regions were observed, indicating different regulatory mechanisms or allelic frequencies among the breeds. In addition, hotspot regions regulating the expression of several genes were detected. Our results provide new data to better understand the functional regulatory mechanisms of lipid metabolism genes in muscle.

Genome-wide detection of copy number variants in European autochthonous and commercial pig breeds by whole-genome sequencing of DNA pools identified breed-characterising copy number states

In this study, we identified copy number variants (CNVs) in 19 European autochthonous pig breeds and in two commercial breeds (Italian Large White and Italian Duroc) that represent important genetic resources for this species. The genome of 725 pigs was sequenced using a breed-specific DNA pooling approach (30-35 animals per pool) obtaining an average depth per pool of 42×. This approach maximised CNV discovery as well as the related copy number states characterising, on average, the analysed breeds. By mining more than 17.5 billion reads, we identified a total of 9592 CNVs (~683 CNVs per breed) and 3710 CNV regions (CNVRs; 1.15% of the reference pig genome), with an average of 77 CNVRs per breed that were considered as private. A few CNVRs were analysed in more detail, together with other information derived from sequencing data. For example, the CNVR encompassing the KIT gene was associated with coat colour phenotypes in the analysed breeds, confirming the role of the multiple copies in determining breed-specific coat colours. The CNVR covering the MSRB3 gene was associated with ear size in most breeds. The CNVRs affecting the ELOVL6 and ZNF622 genes were private features observed in the Lithuanian Indigenous Wattle and in the Turopolje pig breeds respectively. Overall, the genome variability unravelled here can explain part of the genetic diversity among breeds and might contribute to explain their origin, history and adaptation to a variety of production systems.

Genome-wide association study identifies quantitative trait loci regions involved in muscle acidic profile in Large White heavy pigs

The widespread use of genome-wide association studies resulted in the discovery of genomic regions associated with fatty acid (FA) composition in different porcine tissues, but little information exists about the genes involved in FA composition of meat obtained from heavy pigs selected for the production of Italian dry-cured hams. To this objective, we genotyped with a single nucleotide polymorphism (SNP) panel 795 Italian Large White heavy pigs to identify the markers and genomic regions associated with Semimembranosus muscle FA profile. Heritability estimates for intramuscular fat FA profile were of low-to-moderate magnitude, suggesting that these traits may be improved with genomic selection. On the whole, 45 SNPs were significantly associated with 14 FAs, and 4 of them (ALGA008109, ALGA0081097, CASI0010164 and SIRI0000267) were associated with more than 1 FA. The palmitoleic : palmitic and oleic : stearic ratios displayed the highest number of significant markers and the most significant associations (Bonferroni adjusted P < 5.00E-07). Of particular interest, the palmitoleic : palmitic ratio was strongly associated with markers located at 111 to 114 Mb on chromosome 14, in the same chromosomal region where Stearoyl-CoA desaturase Δ9 (SCD) gene is located. Several significant chromosomal regions were found; some of them harbour key genes playing pivotal roles in FA desaturation and elongation, such as SCD and some members of the Elongation of Very Long-Chain FA (ELOVL) gene family. The results suggest that the identification of causal mutations in these regions may provide a set of markers useful for selection schemes aimed at improving FA composition in pork products.

Effects of long-term feeding of rapeseed meal on skeletal muscle transcriptome, production efficiency and meat quality traits in Norwegian Landrace growing-finishing pigs

This study was performed to investigate the effects of dietary inclusion of 20% rapeseed meal (RSM) as an alternative to soybean meal (SBM) in a three-month feeding experiment with growing finishing pigs. Dietary alteration affected growth performance, several carcass traits and transcriptional responses in the skeletal muscle, but did not affect measured meat quality traits. In general, pigs fed the RSM test diet exhibited reduced growth performance compared to pigs on SBM control diet. Significant transcriptional changes in the skeletal muscle of growing pigs fed RSM diet were likely the consequence of an increased amount of fiber and higher polyunsaturated fatty acids, and presence of bioactive phytochemicals, such as glucosinolates. RNAseq pipeline using Tophat2-Cuffdiff identified 57 upregulated and 63 downregulated genes in RSM compared to SBM pigs. Significantly enriched among downregulated pathways was p53-mediated signalling involved in cellular proliferation, while activation of negative growth regulators (IER5, KLF10, BTG2, KLF11, RETREG1, PRUNE2) in RSM fed pigs provided further evidence for reduced proliferation and increased cellular death, in accordance with the observed reduction in performance traits. Upregulation of well-known metabolic controllers (PDK4, UCP3, ESRRG and ESRRB), involved in energy homeostasis (glucose and lipid metabolism, and mitochondrial function), suggested less available energy and nutrients in RSM pigs. Furthermore, several genes supported more pronounced proteolysis (ABTB1, OTUD1, PADI2, SPP1) and reduced protein synthesis (THBS1, HSF4, AP1S2) in RSM muscle tissue. In parallel, higher levels of NR4A3, PDK4 and FGF21, and a drop in adropin, ELOVL6 and CIDEC/FSP27 indicated increased lipolysis and fatty acid oxidation, reflective of lower dressing percentage. Finally, pigs exposed to RSM showed greater expression level of genes responsive to oxidative stress, indicated by upregulation of GPX1, GPX2, and TXNIP.

Deep transcriptome analysis using RNA-Seq suggests novel insights into molecular aspects of fat-tail metabolism in sheep

Fat-tail content of sheep breeds is varied and the molecular mechanisms regulating fat-tail development have not been well characterized. Aiming at better identifying the important candidate genes and their functional pathways contributing to fat deposition in the tail, a comparative transcriptome analysis was performed between fat- (Lori-Bakhtiari) and thin-tailed (Zel) Iranian sheep breeds using RNA-seq. The experiment was conducted on six male lambs (three lambs per each breed) at seven months of age. Four different combinations of aligners and statistical methods including Hisat2 + edgeR, Hisat2 + DESeq2, STAR + edgeR and STAR + DESeq2 were used to identify the differentially expressed genes (DEGs). The DEGs were selected for functional enrichment analysis and protein-protein interaction (PPI) network construction. Module analysis was also conducted to mine the functional sub-networks from the PPI network. In total, 264 genes including 80 up- and 184 down-regulated genes were identified as DEGs. The RNA-Seq results were validated by Q-RT-PCR. Functional analysis of DEGs and the module analysis of PPI network demonstrated that in addition to pathways affecting lipid metabolism, a series of enriched functional terms related to "response to interleukin", "MAPK signaling pathways", "Wnt signaling pathway", "ECM-receptor interaction", "regulation of actin cytoskeleton", and "response to cAMP" might contribute to the deposition of fat in tails of sheep. Overall results using RNA-Seq analysis characterized important candidate genes involved in the fatty acid metabolism and regulation of fat deposition, suggesting novel insights into molecular aspects of fat-tail metabolism in sheep. Selected DEGs should be further investigated as potential markers associated with the fat-tail development in sheep breeds.

Uncovering the mechanism whereby dietary nicotinic acid increases the intramuscular fat content in finishing steers by RNA sequencing analysis

In our previous study, we found that a higher dosage of nicotinic acid (NA) in the diet dramatically increases intramuscular fat (IMF) content and improves meat quality in finishing steers. We hypothesised that increased IMF results from the regulation of genes associated with adipogenesis. To address this hypothesis, RNA-seq was used to investigate gene-expression profiles of longissimus muscles from the same 16 cattle that were also used in our previous study and treated with or without dietary NA. Four cDNA libraries were constructed and sequenced. The repeatability and reproducibility of RNA-seq data were confirmed by quantitative reverse-transcription polymerase-chain reaction. In total, 123 differentially expressed genes (DEGs) were identified between longissimus muscles treated and those not treated with dietary NA. Of the 123 DEGs, 117 genes were upregulated by the NA treatment. These DEGs were enriched in 21 pathways, including the extracellular matrix (ECM) –receptor interaction, PPAR signalling pathway, adipocytokine signalling pathway and transforming growth factor-β signalling pathway, all of which are associated with lipid metabolism. Furthermore, candidate genes related to adipocyte differentiation and adipogenesis (PLIN1, PLIN2, ADPN, LEP, LCN2 and SOCS3), lipid metabolism (FABP4, RBP4, GAL, ANXA1, ANXA2 and PTX3) and fatty acid synthesis and esterification (ELOVL6, ACSM1, SOT1 and PTGIS) were upregulated in the NA group. Three genes involved in glucose metabolism (PGAM1, UGDH and GLUT3) were also transcriptionally upregulated. However, MYH4 that encodes glycolytic Type IIb muscle fibres was downregulated by dietary NA. These gene expression results indicated a confirmation of our hypothesis that dietary NA increases the IMF content of longissimus muscle through upregulating the expression of the genes related to adipocyte differentiation, adipogenesis and lipid and glucose metabolism.

Genetic variants of fatty acid elongase 6 in Chinese Holstein cow

In the dairy industry, genetic variants have contributed to the improvement of milk production traits. Fatty acid elongase 6 (ELOVL6), which elongates saturated fatty acids (SFAs) and monounsaturated fatty acids (MUFAs), plays a distinct role in the balance of long-chain fatty acids composition in animals. ELOVL6 catalyzes the elongation of palmitic acids (C16:0) which is the most common saturated fatty acid found in animals and also an essential precursor to synthesize other long-chain fatty acids. However, the genetic variant research of bovine ELOVL6 on milk quality traits is still unknown. Therefore, our study aimed to detect the single nucleotide polymorphism (SNPs) of bovine ELOVL6 and explored the relationship between SNPs and milk quality traits including milk yield, fat content, protein content and somatic cell score. In this study, three SNPs, including SNP1 (g16379651A>G), SNP2 (g16458976A>G) and SNP3 (g16511290A>G), have been identified in intron 3 and 3'UTR regions of ELOVL6 in Chinese Holstein (CH) cows. Besides, the results of genetic diversity analysis, linkage disequilibrium and haplotype analysis indicated that these SNPs presented moderate polymorphisms which reflected relatively high genetic diversity. No strong linkage among these SNPs were detected in sampled population of cows. Moreover, the results of correlation analyses demonstrated that these SNPs of bovine ELOVL6 were significantly related to milk yield (P < 0.05). The SNP1 was also correlated with somatic cell score, whereas the SNP3 was associated with fat content. The 21 combined genotypes (diplotypes) were highly significantly correlated (P < 0.01) with milk yield. These results revealed that the genetic variants of bovine ELOVL6 influenced the milk production of CH cows. Hence, the three SNPs could be regarded as molecular markers in marker-assisted selection (MAS) of the dairy cow breeding.

Detection of genetic variants between different Polish Landrace and Puławska pigs by means of RNA-seq analysis

Variant calling analysis based on RNA sequencing data provides information about gene variants. RNA-seq is cheaper and faster than is DNA sequencing. However, it requires individual hard filters during data processing due to post-transcriptional modifications such as splicing and RNA editing. In the present study, RNA-seq transcriptome data on two Polish pig breeds (Puławska, PUL, n = 8, and Polish Landrace, PL, n = 8) were included. The pig breeds are significantly different with regard to meat qualities such as texture, water exudation, growth traits and fat content in carcasses. A total of 2451 significant mutations were identified by a chi square tests, and functional analysis was carried out using Panther, KEGG and Kobas. Interesting missense gene variants and mutations located in regulatory regions were found in a few genes related to fatty acid metabolism and lipid storage such as ACSL5, ALDH3A2, FADS1, SCD, PLA2G12A and ATGL. A validation of mutational influences on pig traits was performed for ALDH3A2, ATGL, PLA2G12A and MYOM1 variants using association analysis including 215 pigs of the PL and PUL breeds. The ALDH3A2ENSSSCT00000019636.2:c.470T>C polymorphism was found to affect the weight of the ham and loin eye area. In turn, an ENSSSCT00000004091.2:c.2836G>A MYOM1 mutation, which could be implicated in myofibrillar network organisation, had an effect on meatiness and loin texture parameters. The study aimed to estimate the usefulness of RNA-seq results for a purpose other than differentially expressed gene analysis. The analysis performed indicated interesting gene variants that could be used in the future as markers during selection.

Association study highlights the influence of ELOVL fatty acid elongase 6 gene region on backfat fatty acid composition in Large White pig breed

Dietary fatty acid (FA) composition has an impact on human health. There is an increasing request from consumers for healthier food and pork industry must respond to it without worsening performance and the technological properties of pork products. The inclusion of genetic markers for carcass FA composition in pig selection schemes could be a useful tool to reach the right balance between unsaturated and saturated FAs to satisfy market demands. With the aim of finding genomic regions associated with porcine backfat FA composition, a genome-wide association study was performed on 798 Italian Large White pigs genotyped using Illumina PorcineSNP60 k. The strongest associations with backfat contents of palmitic, palmitoleic, oleic, medium-chain and long-chain FAs were found for the Sus scrofa chromosome (SSC) 8 region located at 119 to 122 Mb, where the gene ELOVL FA elongase 6 is mapped. Palmitic, palmitoleic, stearic and oleic acid contents were also found associated with SSC14, in particular with the genomic region at 121 to 124 Mb, where stearoyl-CoA desaturase Δ9 gene lies. On the other hand, the genomic regions associated with backfat contents of arachidic, arachidonic, n-6 and n-3 FAs showed to harbour mainly genes involved in dietary lipids and carbohydrates digestion, absorption and utilisation. To our knowledge, this is the first study performed in Large White pigs identifying markers and genomic regions associated with backfat FA composition. The results validate in Large White some associations previously detected in other pig breeds and indicate the involvement of distinct metabolic pathways in the deposition pattern of essential and non-essential FAs.

Genome-wide identification and analysis of elongase of very long chain fatty acid genes in the silkworm, Bombyx mori

Very long chain fatty acids (VLCFAs), such as sphingolipids, are components of cellular lipids, which are essential for cell proliferation. Mutations in the genes that encode proteins participating in VLCFA biosynthesis may cause inherited diseases, such as macular degeneration. Elongases of very long chain fatty acid (ELOVL) are enzymes that are involved in the biosynthesis of VLCFAs. Here, a total of 13 ELOVL genes, distributed across three chromosomes, were identified in the silkworm genome; all the ELOVL members contain a distinct ELO domain and a conserved HXXHH motif. Phylogenetic reconstruction was performed to analyze the evolutionary relationships among different species and to predict gene functions. The 13 ELOVL genes were assigned to the ELOVL3/6, ELOVL1/7, and ELOVL4 clades. Microarray and semiquantitative PCR analyses indicated that these genes are differentially expressed among various tissues, in turn suggesting functional divergence in the growth and development of each tissue. Further investigation showed that the expression level of the BGIBMGA000424 gene is significantly negatively correlated with the cocoon-shell weight among different silkworm strains. Taken together, the present study is the first comprehensive analysis of ELOVL genes in silkworm, and the results may serve as a foundation for further analysis of the physiological functions of ELOVL genes in silkworm.

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.

The effect of lipid metabolism-related genes on intramuscular fat content and fatty acid composition in multiple muscles

Intramuscular fat content (IMF) and fatty acid composition are two important factors that have a significant effect on meat quality. Previous studies about lipid deposition mainly focussed on breed effects, but the regulation mechanism of lipid metabolism among multiple muscles is not clear. Here, we hypothesised that there are correlations between lipid metabolism-related genes and muscle fibre types composition and lipid deposition in multiple muscles. We analysed the relationship between the expression of 18 lipid metabolism-related genes and muscle fibre types composition, and their relation with IMF and fatty acid composition in 14 different muscles. The IMF content and fatty acid composition were significantly different among the muscle tissues (P < 0.01). IMF was significantly higher in the trapezius and semitendinosus muscles compared with the others (P < 0.05); the content was about four times higher than that of the peroneal longus. Moreover, the trapezius and masseter had a higher monounsaturated fatty acid (MUFA) : saturated fatty acid (SFA) (>1.48) (P < 0.05) and polyunsaturated fatty acid (PUFA) : SFA (>0.45) (P < 0.05), which are more conducive to human health. Second, the expression levels of seven genes, AdPLA (r = 0.605, P < 0.05), DGAT2 (r = 0.553, P < 0.05), FABP4 (r = 0.637, P < 0.05), ELOVL6 (r = 0.57, P < 0.05), FASN (r = 0.556, P < 0.05), PPARγ (r = –0.51, P < 0.05) and SCD (r = 0.579, P < 0.05) were associated with IMF. In addition, the expression of FASN was positively correlated with MUFA (r = 0.556, P < 0.05) and total fatty acids (r = 0.547, P < 0.05). Finally, the MyHC IIa content (IIa) was positively correlated with MUFA (r = 0.56, P < 0.05) but negatively correlated with SFA (r = –0.553, P < 0.05) in different muscles. Besides, MUFA : SFA was positively correlated with IIa (r = 0.609, P < 0.05), although negatively correlated with MyHC IIx content (IIx) (r = –0.566, P < 0.05). These results added beneficial information to improve meat quality, and our understanding of the mechanism of fat deposition and fatty acid composition regulation. They also provide potential genetic markers for the study of muscular fatty acid composition.

Genetic mapping using 1.4M SNP array refined loci for fatty acid composition traits in Chinese Erhualian and Bamaxiang pigs

Chinese indigenous pigs display marked genetic and phenotypic differences compared with western commercial pigs. In this study, we tested the association between 660K SNPs and longissimus muscle fatty acid composition traits in Chinese Erhualian (n = 331) and Bamaxiang (n = 315) pigs based on a customized 1.4 million SNP array. We identified a total of 64 significant associations for 20 fatty acid composition traits at the p-value threshold of 1 × 10^-6 among which 42 associations in low linkage disequilibrium (r²  < .2) with previously reported loci were considered novel. We substantially improved the strength and precision of the associations at four previously detected loci near FADS2, ELOVL7, ELOVL6 and FASN genes, facilitating follow-up candidate gene studies. Moreover, we also identified loci near ABCD2, ACSBG1, ELOVL5, HPGDS, DAGT2, ACAD10 and ACSL1 genes with function relevant to metabolism of fatty acids. In this study, valuable genetic variants and candidate genes associated with fatty acid composition traits were identified in Erhualian and Bamaxiang pigs. Some identified loci could be used to improve pork nutrition in pig breeding practice. Using the SNP array with higher marker density and less ascertainment bias improved QTL detection power and precision in Chinese indigenous pigs.

Genome-wide association study confirm major QTL for backfat fatty acid composition on SSC14 in Duroc pigs

Background

Fatty acid composition contributes importantly to meat quality and is essential to the nutritional value of the meat. Identification of genetic factors underlying levels of fatty acids can be used to breed for pigs with healthier meat. The aim of this study was to conduct genome-wide association studies (GWAS) to identify QTL regions affecting fatty acid composition in backfat from the pig breeds Duroc and Landrace.

Results

Using data from the Axiom porcine 660 K array, we performed GWAS on 454 Duroc and 659 Landrace boars for fatty acid phenotypes measured by near-infrared spectroscopy (NIRS) technology (C16:0, C16:1n-7, C18:0, C18:1n-9, C18:2n-6, C18:3n-3, total saturated fatty acids, monounsaturated fatty acids and polyunsaturated fatty acids). Two QTL regions on SSC4 and SSC14 were identified in Duroc for the de novo synthesized fatty acids traits, whereas one QTL on SSC8 was detected in Landrace for C16:1n-7. The QTL region on SSC14 has been reported in previous studies and a putative causative mutation has been suggested in the promoter region of the SCD gene. Whole genome re-sequencing data was used for genotype imputation and to fine map the SSC14 QTL region in Norwegian Duroc. This effort confirms the location of the QTL on this chromosome as well as suggesting other putative candidate genes in the region. The most significant single nucleotide polymorphisms (SNPs) located on SSC14 explain between 55 and 76% of the genetic variance and between 27 and 54% of the phenotypic variance for the de novo synthesized fatty acid traits in Norwegian Duroc. For the QTL region on SSC8 in Landrace, the most significant SNP explained 19% of the genetic variance and 5% of the phenotypic variance for C16:1n-7.

Conclusions

This study confirms a major QTL affecting fatty acid composition on SSC14 in Duroc, which can be used in genetic selection to increase the level of fatty acid desaturation. The SSC14 QTL was not segregating in the Landrace population, but another QTL on SSC8 affecting C16:1n-7 was identified and might be used to increase the level of desaturation in meat products from this breed.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3752-0) contains supplementary material, which is available to authorized users.

Integration of liver gene co-expression networks and eGWAs analyses highlighted candidate regulators implicated in lipid metabolism in pigs

In the present study, liver co-expression networks and expression Genome Wide Association Study (eGWAS) were performed to identify DNA variants and molecular pathways implicated in the functional regulatory mechanisms of meat quality traits in pigs. With this purpose, the liver mRNA expression of 44 candidates genes related with lipid metabolism was analysed in 111 Iberian x Landrace backcross animals. The eGWAS identified 92 eSNPs located in seven chromosomal regions and associated with eight genes: CROT, CYP2U1, DGAT1, EGF, FABP1, FABP5, PLA2G12A, and PPARA. Remarkably, cis-eSNPs associated with FABP1 gene expression which may be determining the C18:2(n-6)/C18:3(n-3) ratio in backfat through the multiple interaction of DNA variants and genes were identified. Furthermore, a hotspot on SSC8 associated with the gene expression of eight genes was identified and the TBCK gene was pointed out as candidate gene regulating it. Our results also suggested that the PI3K-Akt-mTOR pathway plays an important role in the control of the analysed genes highlighting nuclear receptors as the NR3C1 or PPARA. Finally, sex-dimorphism associated with hepatic lipid metabolism was identified with over-representation of female-biased genes. These results increase our knowledge of the genetic architecture underlying fat composition traits.

Genome scan linkage analysis identifies a major quantitative trait loci for fatty acid composition in longissimus dorsi muscle in an F2 intercross between Landrace and Korean native pigs

OBJECTIVE

This study was conducted to locate quantitative trait loci (QTL) influencing fatty acid (FA) composition in a large F₂ intercross between Landrace and Korean native pigs.

METHODS

Eighteen FA composition traits were measured in more than 960 F₂ progeny. All experimental animals were genotyped with 165 microsatellite markers located throughout the pig autosomes.

RESULTS

We detected 112 QTLs for the FA composition; Forty seven QTLs reached the genome-wide significant threshold. In particular, we identified a cluster of highly significant QTLs for FA composition on SSC12. QTL for polyunsaturated fatty acid on pig chromosome 12 (F-value = 97.2 under additive and dominance model, nominal p-value 3.6×10^-39) accounted for 16.9% of phenotypic variance. In addition, four more QTLs for C18:1, C18:2, C20:4, and monounsaturated fatty acids on the similar position explained more than 10% of phenotypic variance.

CONCLUSION

Our findings of a major QTL for FA composition presented here could provide helpful information to locate causative variants to improve meat quality traits in pigs.

Genetic Marker Discovery in Complex Traits: A Field Example on Fat Content and Composition in Pigs

Among the large number of attributes that define pork quality, fat content and composition have attracted the attention of breeders in the recent years due to their interaction with human health and technological and sensorial properties of meat. In livestock species, fat accumulates in different depots following a temporal pattern that is also recognized in humans. Intramuscular fat deposition rate and fatty acid composition change with life. Despite indication that it might be possible to select for intramuscular fat without affecting other fat depots, to date only one depot-specific genetic marker (PCK1 c.2456C>A) has been reported. In contrast, identification of polymorphisms related to fat composition has been more successful. For instance, our group has described a variant in the stearoyl-coA desaturase (SCD) gene that improves the desaturation index of fat without affecting overall fatness or growth. Identification of mutations in candidate genes can be a tedious and costly process. Genome-wide association studies can help in narrowing down the number of candidate genes by highlighting those which contribute most to the genetic variation of the trait. Results from our group and others indicate that fat content and composition are highly polygenic and that very few genes explain more than 5% of the variance of the trait. Moreover, as the complexity of the genome emerges, the role of non-coding genes and regulatory elements cannot be disregarded. Prediction of breeding values from genomic data is discussed in comparison with conventional best linear predictors of breeding values. An example based on real data is given, and the implications in phenotype prediction are discussed in detail. The benefits and limitations of using large SNP sets versus a few very informative markers as predictors of genetic merit of breeding candidates are evaluated using field data as an example.

Deciphering the regulation of porcine genes influencing growth, fatness and yield-related traits through genetical genomics

Genetical genomics approaches aim at identifying quantitative trait loci for molecular traits, also known as intermediate phenotypes, such as gene expression, that could link variation in genetic information to physiological traits. In the current study, an expression GWAS has been carried out on an experimental Iberian × Landrace backcross in order to identify the genomic regions regulating the gene expression of those genes whose expression is correlated with growth, fat deposition, and premium cut yield measures in pig. The analyses were conducted exploiting Porcine 60K SNP BeadChip genotypes and Porcine Expression Microarray data hybridized on mRNA from Longissimus dorsi muscle. In order to focus the analysis on productive traits and reduce the number of analyses, only those probesets whose expression showed significant correlation with at least one of the seven phenotypes of interest were selected for the eGWAS. A total of 63 eQTL regions were identified with effects on 36 different transcripts. Those eQTLs overlapping with phenotypic QTLs on SSC4, SSC9, SSC13, and SSC17 chromosomes previously detected in the same animal material were further analyzed. Moreover, candidate genes and SNPs were analyzed. Among the most promising results, a long non-coding RNA, ALDBSSCG0000001928, was identified, whose expression is correlated with premium cut yield. Association analysis and in silico sequence domain annotation support TXNRD3 polymorphisms as candidate to regulate ALDBSSCG0000001928 expression, which can be involved in the transcriptional regulation of surrounding genes, affecting productive and meat quality traits.

Genetic architecture of fatty acid composition in the longissimus dorsi muscle revealed by genome-wide association studies on diverse pig populations

Background

Fatty acid composition in muscle is an important factor that affects the nutritive value and taste of pork. To investigate the genetic architecture of fatty acid composition of pork, we measured fatty acid contents in longissimus dorsi muscle of 1244 pigs from three divergent populations and conducted genome-wide association studies (GWAS) for fatty acid contents.

Results

We detected 26 genome-wide significant quantitative trait loci (QTL) on eight chromosomes (SSC for Sus scrofa) for eight fatty acids. These loci not only replicated previously reported QTL for C18:0 on SSC14 and C20:0 on SSC16, but also included several novel QTL such as those for C20:1 on SSC7, C14:0 on SSC9, and C14:0, C16:0 and C16:1 on SSC12. Furthermore, we performed a meta-analysis of GWAS on five populations, including the three populations that were investigated in this study and two additional populations that we had previously examined. This enhanced the strength of the associations detected between fatty acid composition and several marker loci, especially for those for C18:0 on SSC14 and C20:0 on SSC16. The genes ELOVL5, ELOVL6, ELOVL7, FASN, SCD and THRSP, which have functions that are directly relevant to fatty acid metabolism, are proximal to the top associated markers at six significant QTL.

Conclusions

The findings improve our understanding of the genetic architecture of fatty acid composition in pork and contribute to further fine-map and characterize genes that influence fatty acid composition.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-016-0184-2) contains supplementary material, which is available to authorized users.

Genome-wide analysis of microRNAs identifies the lipid metabolism pathway to be a defining factor in adipose tissue from different sheep

MicroRNAs are short (17–24 nt) non-coding RNAs that are involved in post-transcriptional regulation of gene expression in multicellular organisms by affecting both the stability and translation of mRNAs. In recent years, deep sequencing of the transcriptome is increasingly being utilized with the promise of higher sensitivity for the identification of differential expression patterns as well as the opportunity to discover new transcripts, including new alternative isoforms and miRNAs. Here, we utilized RNA-seq technology to perform a genome-wide analysis of miRNAs from the adipose tissue of the two species of sheep to look for clues that might explain the fat deposition differences between the sheep. The RNA-seq analysis detected 3132 miRNAs from the adipose tissue of the Small-tail Han and Dorset sheep, of which 2893 were defined as potential new miRNAs. In addition, 54 miRNAs were differentially expressed between the two breeds of sheep. Gene ontology and pathway analyses of the predicted target genes that negatively associated with the differentially expressed miRNAs revealed that there was less active lipid metabolism in the adipose tissue of Small Tail Han sheep. This study can help understand the underling mechanisms responsible for the morphological differences related to fat deposition between two breeds of sheep.

Integration of transcriptome and whole genomic resequencing data to identify key genes affecting swine fat deposition

Fat deposition is highly correlated with the growth, meat quality, reproductive performance and immunity of pigs. Fatty acid synthesis takes place mainly in the adipose tissue of pigs; therefore, in this study, a high-throughput massively parallel sequencing approach was used to generate adipose tissue transcriptomes from two groups of Songliao black pigs that had opposite backfat thickness phenotypes. The total number of paired-end reads produced for each sample was in the range of 39.29-49.36 millions. Approximately 188 genes were differentially expressed in adipose tissue and were enriched for metabolic processes, such as fatty acid biosynthesis, lipid synthesis, metabolism of fatty acids, etinol, caffeine and arachidonic acid and immunity. Additionally, many genetic variations were detected between the two groups through pooled whole-genome resequencing. Integration of transcriptome and whole-genome resequencing data revealed important genomic variations among the differentially expressed genes for fat deposition, for example, the lipogenic genes. Further studies are required to investigate the roles of candidate genes in fat deposition to improve pig breeding programs.

Epigenetic regulation of the ELOVL6 gene is associated with a major QTL effect on fatty acid composition in pigs

Background

In previous studies on an Iberian x Landrace cross, we have provided evidence that supported the porcine ELOVL6 gene as the major causative gene of the QTL on pig chromosome 8 for palmitic and palmitoleic acid contents in muscle and backfat. The single nucleotide polymorphism (SNP) ELOVL6:c.-533C > T located in the promoter region of ELOVL6 was found to be highly associated with ELOVL6 expression and, accordingly, with the percentages of palmitic and palmitoleic acids in longissimus dorsi and adipose tissue. The main goal of the current work was to further study the role of ELOVL6 on these traits by analyzing the regulation of the expression of ELOVL6 and the implication of ELOVL6 polymorphisms on meat quality traits in pigs.

Results

High-throughput sequencing of BAC clones that contain the porcine ELOVL6 gene coupled to RNAseq data re-analysis showed that two isoforms of this gene are expressed in liver and adipose tissue and that they differ in number of exons and 3’UTR length. Although several SNPs in the 3’UTR of ELOVL6 were associated with palmitic and palmitoleic acid contents, this association was lower than that previously observed with SNP ELOVL6:c.-533C > T. This SNP is in full linkage disequilibrium with SNP ELOVL6:c.-394G > A that was identified in the binding site for estrogen receptor alpha (ERα). Interestingly, the ELOVL6:c.-394G allele is associated with an increase in methylation levels of the ELOVL6 promoter and with a decrease of ELOVL6 expression. Therefore, ERα is clearly a good candidate to explain the regulation of ELOVL6 expression through dynamic epigenetic changes in the binding site of known regulators of ELOVL6 gene, such as SREBF1 and SP1.

Conclusions

Our results strongly suggest the ELOVL6:c.-394G > A polymorphism as the causal mutation for the QTL on pig chromosome 8 that affects fatty acid composition in pigs.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-015-0111-y) contains supplementary material, which is available to authorized users.

A splice mutation in the PHKG1 gene causes high glycogen content and low meat quality in pig skeletal muscle

Glycolytic potential (GP) in skeletal muscle is economically important in the pig industry because of its effect on pork processing yield. We have previously mapped a major quantitative trait loci (QTL) for GP on chromosome 3 in a White Duroc × Erhualian F₂ intercross. We herein performed a systems genetic analysis to identify the causal variant underlying the phenotype QTL (pQTL). We first conducted genome-wide association analyses in the F₂ intercross and an F₁₉ Sutai pig population. The QTL was then refined to an 180-kb interval based on the 2-LOD drop method. We then performed expression QTL (eQTL) mapping using muscle transcriptome data from 497 F₂ animals. Within the QTL interval, only one gene (PHKG1) has a cis-eQTL that was colocolizated with pQTL peaked at the same SNP. The PHKG1 gene encodes a catalytic subunit of the phosphorylase kinase (PhK), which functions in the cascade activation of glycogen breakdown. Deep sequencing of PHKG1 revealed a point mutation (C>A) in a splice acceptor site of intron 9, resulting in a 32-bp deletion in the open reading frame and generating a premature stop codon. The aberrant transcript induces nonsense-mediated decay, leading to lower protein level and weaker enzymatic activity in affected animals. The mutation causes an increase of 43% in GP and a decrease of>20% in water-holding capacity of pork. These effects were consistent across the F₂ and Sutai populations, as well as Duroc × (Landrace × Yorkshire) hybrid pigs. The unfavorable allele exists predominantly in Duroc-derived pigs. The findings provide new insights into understanding risk factors affecting glucose metabolism, and would greatly contribute to the genetic improvement of meat quality in Duroc related pigs.

Glycogen storage diseases (GSD) are a group of inherited disorders characterized by storage of excess glycogen, which are mainly caused by the abnormality of a particular enzyme essential for releasing glucose from glycogen. GSD-like conditions have been described in a wide variety of species. Pigs are a valuable model for the study of human GSD. Moreover, pigs affected by GSD usually produce inferior pork with a lower ultimate pH (so-called “acid meat”) and less processing yield due to post-mortem degradation of the excess glycogen. So far, only one causal variant, PRKAG3 R225Q, has been identified for GSD in pigs. Here we reported a loss-of-function mutation in the PHKG1 gene that causes the deficiency of the glycogen breakdown, consequently leading to GSD and acid meat in Duroc-sired pigs. Eliminating the undesirable mutation from the breeding stock by a diagnostic DNA test will greatly reduce the incidence of GSD and significantly improve pork quality and productivity in the pig.

New insight into the SSC8 genetic determination of fatty acid composition in pigs

Background

Fat content and fatty acid composition in swine are becoming increasingly studied because of their effect on sensory and nutritional quality of meat. A QTL (quantitative trait locus) for fatty acid composition in backfat was previously detected on porcine chromosome 8 (SSC8) in an Iberian x Landrace F₂ intercross. More recently, a genome-wide association study detected the same genomic region for muscle fatty acid composition in an Iberian x Landrace backcross population. ELOVL6, a strong positional candidate gene for this QTL, contains a polymorphism in its promoter region (ELOVL6:c.-533C < T), which is associated with percentage of palmitic and palmitoleic acids in muscle and adipose tissues. Here, a combination of single-marker association and the haplotype-based approach was used to analyze backfat fatty acid composition in 470 animals of an Iberian x Landrace F₂ intercross genotyped with 144 SNPs (single nucleotide polymorphisms) distributed along SSC8.

Results

Two trait-associated SNP regions were identified at 93 Mb and 119 Mb on SSC8. The strongest statistical signals of both regions were observed for palmitoleic acid (C16:1(n-7)) content and C18:0/C16:0 and C18:1(n-7)/C16:1(n-7) elongation ratios. MAML3 and SETD7 are positional candidate genes in the 93 Mb region and two novel microsatellites in MAML3 and nine SNPs in SETD7 were identified. No significant association for the MAML3 microsatellite genotypes was detected. The SETD7:c.700G > T SNP, although statistically significant, was not the strongest signal in this region. In addition, the expression of MAML3 and SETD7 in liver and adipose tissue varied among animals, but no association was detected with the polymorphisms in these genes. In the 119 Mb region, the ELOVL6:c.-533C > T polymorphism showed a strong association with percentage of palmitic and palmitoleic fatty acids and elongation ratios in backfat.

Conclusions

Our results suggest that the polymorphisms studied in MAML3 and SETD7 are not the causal mutations for the QTL in the 93 Mb region. However, the results for ELOVL6 support the hypothesis that the ELOVL6:c.-533C > T polymorphism has a pleiotropic effect on backfat and intramuscular fatty acid composition and that it has a role in the determination of the QTL in the 119 Mb region.

From SNP co-association to RNA co-expression: novel insights into gene networks for intramuscular fatty acid composition in porcine

Background

Fatty acids (FA) play a critical role in energy homeostasis and metabolic diseases; in the context of livestock species, their profile also impacts on meat quality for healthy human consumption. Molecular pathways controlling lipid metabolism are highly interconnected and are not fully understood. Elucidating these molecular processes will aid technological development towards improvement of pork meat quality and increased knowledge of FA metabolism, underpinning metabolic diseases in humans.

Results

The results from genome-wide association studies (GWAS) across 15 phenotypes were subjected to an Association Weight Matrix (AWM) approach to predict a network of 1,096 genes related to intramuscular FA composition in pigs. To identify the key regulators of FA metabolism, we focused on the minimal set of transcription factors (TF) that the explored the majority of the network topology. Pathway and network analyses pointed towards a trio of TF as key regulators of FA metabolism: NCOA2, FHL2 and EP300. Promoter sequence analyses confirmed that these TF have binding sites for some well-know regulators of lipid and carbohydrate metabolism. For the first time in a non-model species, some of the co-associations observed at the genetic level were validated through co-expression at the transcriptomic level based on real-time PCR of 40 genes in adipose tissue, and a further 55 genes in liver. In particular, liver expression of NCOA2 and EP300 differed between pig breeds (Iberian and Landrace) extreme in terms of fat deposition. Highly clustered co-expression networks in both liver and adipose tissues were observed. EP300 and NCOA2 showed centrality parameters above average in the both networks. Over all genes, co-expression analyses confirmed 28.9% of the AWM predicted gene-gene interactions in liver and 33.0% in adipose tissue. The magnitude of this validation varied across genes, with up to 60.8% of the connections of NCOA2 in adipose tissue being validated via co-expression.

Conclusions

Our results recapitulate the known transcriptional regulation of FA metabolism, predict gene interactions that can be experimentally validated, and suggest that genetic variants mapped to EP300, FHL2, and NCOA2 modulate lipid metabolism and control energy homeostasis in pigs.

Genome-wide analysis of porcine backfat and intramuscular fat fatty acid composition using high-density genotyping and expression data

Background

Porcine fatty acid composition is a key factor for quality and nutritive value of pork. Several QTLs for fatty acid composition have been reported in diverse fat tissues. The results obtained so far seem to point out different genetic control of fatty acid composition conditional on the fat deposits. Those studies have been conducted using simple approaches and most of them focused on one single tissue. The first objective of the present study was to identify tissue-specific and tissue-consistent QTLs for fatty acid composition in backfat and intramuscular fat, combining linkage mapping and GWAS approaches and conducted under single and multitrait models. A second aim was to identify powerful candidate genes for these tissue-consistent QTLs, using microarray gene expression data and following a targeted genetical genomics approach.

Results

The single model analyses, linkage and GWAS, revealed over 30 and 20 chromosomal regions, 24 of them identified here for the first time, specifically associated to the content of diverse fatty acids in BF and IMF, respectively. The analyses with multitrait models allowed identifying for the first time with a formal statistical approach seven different regions with pleiotropic effects on particular fatty acids in both fat deposits. The most relevant were found on SSC8 for C16:0 and C16:1(n-7) fatty acids, detected by both linkage and GWAS approaches. Other detected pleiotropic regions included one on SSC1 for C16:0, two on SSC4 for C16:0 and C18:2, one on SSC11 for C20:3 and the last one on SSC17 for C16:0. Finally, a targeted eQTL scan focused on regions showing tissue-consistent effects was conducted with Longissimus and fat gene expression data. Some powerful candidate genes and regions were identified such as the PBX1, RGS4, TRIB3 and a transcription regulatory element close to ELOVL6 gene to be further studied.

Conclusions

Complementary genome scans have confirmed several chromosome regions previously associated to fatty acid composition in backfat and intramuscular fat, but even more, to identify new ones. Although most of the detected regions were tissue-specific, supporting the hypothesis that the major part of genes affecting fatty acid composition differs among tissues, seven chromosomal regions showed tissue-consistent effects. Additional gene expression analyses have revealed powerful target regions to carry the mutation responsible for the pleiotropic effects.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-14-845) contains supplementary material, which is available to authorized users.

Analysis of porcine adipose tissue transcriptome reveals differences in de novo fatty acid synthesis in pigs with divergent muscle fatty acid composition

Background

In pigs, adipose tissue is one of the principal organs involved in the regulation of lipid metabolism. It is particularly involved in the overall fatty acid synthesis with consequences in other lipid-target organs such as muscles and the liver. With this in mind, we have used massive, parallel high-throughput sequencing technologies to characterize the porcine adipose tissue transcriptome architecture in six Iberian x Landrace crossbred pigs showing extreme phenotypes for intramuscular fatty acid composition (three per group).

Results

High-throughput RNA sequencing was used to generate a whole characterization of adipose tissue (backfat) transcriptome. A total of 4,130 putative unannotated protein-coding sequences were identified in the 20% of reads which mapped in intergenic regions. Furthermore, 36% of the unmapped reads were represented by interspersed repeats, SINEs being the most abundant elements. Differential expression analyses identified 396 candidate genes among divergent animals for intramuscular fatty acid composition. Sixty-two percent of these genes (247/396) presented higher expression in the group of pigs with higher content of intramuscular SFA and MUFA, while the remaining 149 showed higher expression in the group with higher content of PUFA. Pathway analysis related these genes to biological functions and canonical pathways controlling lipid and fatty acid metabolisms. In concordance with the phenotypic classification of animals, the major metabolic pathway differentially modulated between groups was de novo lipogenesis, the group with more PUFA being the one that showed lower expression of lipogenic genes.

Conclusions

These results will help in the identification of genetic variants at loci that affect fatty acid composition traits. The implications of these results range from the improvement of porcine meat quality traits to the application of the pig as an animal model of human metabolic diseases.

DAG expression: high-throughput gene expression analysis of real-time PCR data using standard curves for relative quantification

BACKGROUND

Real-time quantitative PCR (qPCR) is still the gold-standard technique for gene-expression quantification. Recent technological advances of this method allow for the high-throughput gene-expression analysis, without the limitations of sample space and reagent used. However, non-commercial and user-friendly software for the management and analysis of these data is not available.

RESULTS

The recently developed commercial microarrays allow for the drawing of standard curves of multiple assays using the same n-fold diluted samples. Data Analysis Gene (DAG) Expression software has been developed to perform high-throughput gene-expression data analysis using standard curves for relative quantification and one or multiple reference genes for sample normalization. We discuss the application of DAG Expression in the analysis of data from an experiment performed with Fluidigm technology, in which 48 genes and 115 samples were measured. Furthermore, the quality of our analysis was tested and compared with other available methods.

CONCLUSIONS

DAG Expression is a freely available software that permits the automated analysis and visualization of high-throughput qPCR. A detailed manual and a demo-experiment are provided within the DAG Expression software at http://www.dagexpression.com/dage.zip.