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

PLA2G12A protects against diet-induced obesity and insulin resistance by enhancing energy expenditure and clearance of circulating triglycerides

Secreted phospholipase A2s are involved in the development of obesity, type 2 diabetes mellitus (T2DM) and cardiovascular disease, which have become serious and growing health concerns worldwide. Integration of genome-wide association study and gene co-expression networks analysis showed that the secreted phospholipase A2 group XIIA (PLA2G12A) may participate in hepatic lipids metabolism. Nevertheless, the role of PLA2G12A in lipid metabolism and its potential mechanism remain elusive. Here, we used AAV9 vector carrying human PLA2G12A gene to exogenously express hPLA2G12A in the liver of mice. We demonstrated that the overexpression of hPLA2G12A resulted in a significant decrease in serum lipid levels in wild-type mice fed with chow diet or high-fat diet (HFD). Moreover, hPLA2G12A treatment protected against diet-induced obesity and insulin resistance in mice fed a HFD. Notably, we found that hPLA2G12A treatment confers protection against obesity and hyperlipidemia independent of its enzymatic activity, but rather by increasing physical activity and energy expenditure. Furthermore, we demonstrated that hPLA2G12A treatment induced upregulation of ApoC2 and Cd36 and downregulation of Angptl8, which contributed to the increase in clearance of circulating triglycerides and hepatic uptake of fatty acids without affecting hepatic de novo lipogenesis, very low-density lipoprotein secretion, or intestinal lipid absorption. Our study highlights the potential of PLA2G12A gene therapy as a promising approach for treating obesity, insulin resistance and T2DM.

Characterization and Function Analysis of miRNA Editing during Fat Deposition in Chinese Indigenous Ningxiang Pigs

This study aimed to identify active miRNA editing sites during adipose development in Ningxiang pigs and analyze their characteristics and functions. Based on small RNA-seq data from the subcutaneous adipose tissues of Ningxiang pigs at four stages-30 days (piglet), 90 days (nursery), 150 days (early fattening), and 210 days (late fattening)-we constructed a developmental map of miRNA editing in the adipose tissues of Ningxiang pigs. A total of 505 miRNA editing sites were identified using the revised pipeline, with C-to-U editing types being the most prevalent, followed by U-to-C, A-to-G, and G-to-U. Importantly, these four types of miRNA editing exhibited base preferences. The number of editing sites showed obvious differences among age groups, with the highest occurrence of miRNA editing events observed at 90 days of age and the lowest at 150 days of age. A total of nine miRNA editing sites were identified in the miRNA seed region, with significant differences in editing levels (p < 0.05) located in ssc-miR-23a, ssc-miR-27a, ssc-miR-30b-5p, ssc-miR-15a, ssc-miR-497, ssc-miR-15b, and ssc-miR-425-5p, respectively. Target gene prediction and KEGG enrichment analyses indicated that the editing of miR-497 might potentially regulate fat deposition by inhibiting adipose synthesis via influencing target binding. These results provide new insights into the regulatory mechanism of pig fat deposition.

Identification of genomic regions, genetic variants and gene networks regulating candidate genes for lipid metabolism in pig muscle

The intramuscular fat content and fatty acid composition of porcine meat have a significant impact on its quality and nutritional value. This research aimed to investigate the expression of 45 genes involved in lipid metabolism in the longissimus dorsi muscle of three experimental pig backcrosses, with a 25% of Iberian background. To achieve this objective, we conducted an expression Genome-Wide Association Study (eGWAS) using gene expression levels in muscle measured by high-throughput real-time qPCR for 45 target genes and genotypes from the PorcineSNP60 BeadChip or Axiom Porcine Genotyping Array and 65 single nucleotide polymorphisms (SNPs) located in 20 genes genotyped by a custom-designed Taqman OpenArray in a cohort of 354 animals. The eGWAS analysis identified 301 eSNPs associated with 18 candidate genes (ANK2, APOE, ARNT, CIITA, CPT1A, EGF, ELOVL6, ELOVL7, FADS3, FASN, GPAT3, NR1D2, NR1H2, PLIN1, PPAP2A, RORA, RXRA and UCP3). Three cis-eQTL (expression quantitative trait loci) were identified for GPAT3, RXRA, and UCP3 genes, which indicates that a genetic polymorphism proximal to the same gene is affecting its expression. Furthermore, 24 trans-eQTLs were detected, and eight candidate regulatory genes were located in these genomic regions. Additionally, two trans-regulatory hotspots in Sus scrofa chromosomes 13 and 15 were identified. Moreover, a co-expression analysis performed on 89 candidate genes and the fatty acid composition revealed the regulatory role of four genes (FABP5, PPARG, SCD, and SREBF1). These genes modulate the levels of α-linolenic, arachidonic, and oleic acids, as well as regulating the expression of other candidate genes associated with lipid metabolism. The findings of this study offer novel insights into the functional regulatory mechanism of genes involved in lipid metabolism, thereby enhancing our understanding of this complex biological process.

Tandem mass tag-based quantitative proteomics analysis reveals the effects of the α-lactalbumin peptides GINY and DQW on lipid deposition and oxidative stress in HepG2 cells

The objective of this study was to investigate the mechanism by which the α-lactalbumin peptides Gly-Ile-Asn-Tyr (GINY) and Asp-Gln-Trp (DQW) ameliorate free fatty acid-induced lipid deposition in HepG2 cells. The results show that GINY and DQW reduced triglyceride, total cholesterol, and free fatty acid levels significantly in free fatty acid-treated HepG2 cells. Based on proteomic analysis, GINY and DQW alleviated lipid deposition and oxidative stress mainly through the peroxisome proliferator-activated receptor (PPAR) pathway, fatty acid metabolism, oxidative phosphorylation, and response to oxidative stress. In vitro experiments confirmed that GINY and DQW upregulated the mRNA and protein expression of fatty acid β-oxidation-related and oxidative stress-related genes, and downregulated the mRNA and protein expression of lipogenesis-related genes by activating peroxisome proliferator-activated receptor α (PPARα). Meanwhile, GINY and DQW reduced free fatty acid-induced lipid droplet accumulation and reactive oxygen species generation, and enhanced the mitochondrial membrane potential and ATP levels. Furthermore, GINY and DQW enhanced carnitine palmitoyl-transferase 1a (CPT-1a) and superoxide dismutase activities, and diminished acetyl-coenzyme A carboxylase 1 (ACC1) and fatty acid synthase (FASN) activities in a PPARα-dependent manner. Interestingly, GW6471 (a PPARα inhibitor) weakened the effects of GINY and DQW on the PPARα pathway. Hence, our findings suggest that GINY and DQW have the potential to alleviate nonalcoholic fatty liver disease by activating the PPARα pathway.

Identification of transcriptional regulatory variants in pig duodenum, liver, and muscle tissues

BACKGROUND

In humans and livestock species, genome-wide association studies (GWAS) have been applied to study the association between variants distributed across the genome and a phenotype of interest. To discover genetic polymorphisms affecting the duodenum, liver, and muscle transcriptomes of 300 pigs from 3 different breeds (Duroc, Landrace, and Large White), we performed expression GWAS between 25,315,878 polymorphisms and the expression of 13,891 genes in duodenum, 12,748 genes in liver, and 11,617 genes in muscle.

RESULTS

More than 9.68 × 1011 association tests were performed, yielding 14,096,080 significantly associated variants, which were grouped in 26,414 expression quantitative trait locus (eQTL) regions. Over 56% of the variants were within 1 Mb of their associated gene. In addition to the 100-kb region upstream of the transcription start site, we identified the importance of the 100-kb region downstream of the 3'UTR for gene regulation, as most of the cis-regulatory variants were located within these 2 regions. We also observed 39,874 hotspot regulatory polymorphisms associated with the expression of 10 or more genes that could modify the protein structure or the expression of a regulator gene. In addition, 2 motifs (5'-GATCCNGYGTTGCYG-3' and a poly(A) sequence) were enriched across the 3 tissues within the neighboring sequences of the most significant single-nucleotide polymorphisms in each cis-eQTL region.

CONCLUSIONS

The 14 million significant associations obtained in this study are publicly available and have enabled the identification of expression-associated cis-, trans-, and hotspot regulatory variants within and across tissues, thus shedding light on the molecular mechanisms of regulatory variations that shape end-trait phenotypes.

Multi-omics analysis dissects the genetic architecture of seed coat content in Brassica napus

Background

Brassica napus is an important vegetable oil source worldwide. Seed coat content is a complex quantitative trait that negatively correlates with the seed oil content in B. napus.

Results

Here we provide insights into the genetic basis of natural variation of seed coat content by transcriptome-wide association studies (TWAS) and genome-wide association studies (GWAS) using 382 B. napus accessions. By population transcriptomic analysis, we identify more than 700 genes and four gene modules that are significantly associated with seed coat content. We also characterize three reliable quantitative trait loci (QTLs) controlling seed coat content by GWAS. Combining TWAS and correlation networks of seed coat content-related gene modules, we find that BnaC07.CCR-LIKE (CCRL) and BnaTT8s play key roles in the determination of the trait by modulating lignin biosynthesis. By expression GWAS analysis, we identify a regulatory hotspot on chromosome A09, which is involved in controlling seed coat content through BnaC07.CCRL and BnaTT8s. We then predict the downstream genes regulated by BnaTT8s using multi-omics datasets. We further experimentally validate that BnaCCRL and BnaTT8 positively regulate seed coat content and lignin content. BnaCCRL represents a novel identified gene involved in seed coat development. Furthermore, we also predict the key genes regulating carbon allocation between phenylpropane compounds and oil during seed development in B. napus.

Conclusions

This study helps us to better understand the complex machinery of seed coat development and provides a genetic resource for genetic improvement of seed coat content in B. napus breeding.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02647-5.

Genomic loci mispositioning in Tmem120a knockout mice yields latent lipodystrophy

Little is known about how the observed fat-specific pattern of 3D-spatial genome organisation is established. Here we report that adipocyte-specific knockout of the gene encoding nuclear envelope transmembrane protein Tmem120a disrupts fat genome organisation, thus causing a lipodystrophy syndrome. Tmem120a deficiency broadly suppresses lipid metabolism pathway gene expression and induces myogenic gene expression by repositioning genes, enhancers and miRNA-encoding loci between the nuclear periphery and interior. Tmem120a-/- mice, particularly females, exhibit a lipodystrophy syndrome similar to human familial partial lipodystrophy FPLD2, with profound insulin resistance and metabolic defects that manifest upon exposure to an obesogenic diet. Interestingly, similar genome organisation defects occurred in cells from FPLD2 patients that harbour nuclear envelope protein encoding LMNA mutations. Our data indicate TMEM120A genome organisation functions affect many adipose functions and its loss may yield adiposity spectrum disorders, including a miRNA-based mechanism that could explain muscle hypertrophy in human lipodystrophy.

Variability in porcine microRNA genes and its association with mRNA expression and lipid phenotypes

BACKGROUND

Mature microRNAs (miRNAs) play an important role in repressing the expression of a wide range of mRNAs. The presence of polymorphic sites in miRNA genes and their corresponding 3'UTR binding sites can disrupt canonical conserved miRNA-mRNA pairings, and thus modify gene expression patterns. However, to date such polymorphic sites in miRNA genes and their association with gene expression phenotypes and complex traits are poorly characterized in pigs.

RESULTS

By analyzing whole-genome sequences from 120 pigs and wild boars from Europe and Asia, we identified 285 single nucleotide polymorphisms (SNPs) that map to miRNA loci, and 109,724 SNPs that are located in predicted 7mer-m8 miRNA binding sites within porcine 3'UTR. In porcine miRNA genes, SNP density is reduced compared with their flanking non-miRNA regions. By sequencing the genomes of five Duroc boars, we identified 12 miRNA SNPs that were subsequently genotyped in their offspring (N = 345, Lipgen population). Association analyses of miRNA SNPs with 38 lipid-related traits and hepatic and muscle microarray expression phenotypes recorded in the Lipgen population were performed. The most relevant detected association was between the genotype of the rs319154814 (G/A) SNP located in the apical loop of the ssc-miR-326 hairpin precursor and PPP1CC mRNA levels in the liver (q-value = 0.058). This result was subsequently confirmed by qPCR (P-value = 0.027). The rs319154814 (G/A) genotype was also associated with several fatty acid composition traits.

CONCLUSIONS

Our findings show a reduced variability of porcine miRNA genes, which is consistent with strong purifying selection, particularly in the seed region that plays a critical role in miRNA binding. Although it is generally assumed that SNPs mapping to the seed region are those with the most pronounced consequences on mRNA expression, we show that a SNP mapping to the apical region of ssc-miR-326 is significantly associated with hepatic mRNA levels of the PPP1CC gene, one of its predicted targets. Although experimental confirmation of such an interaction is reported in humans but not in pigs, this result highlights the need to further investigate the functional effects of miRNA polymorphisms that are located outside the seed region on gene expression in pigs.

Expression analysis of porcine miR-33a/b in liver, adipose tissue and muscle and its potential role in fatty acid metabolism

mir-33a and mir-33b are co-transcribed with the SREBF2 and SREBF1 transcription factors, respectively. The main role of SREBF1 is the regulation of genes involved in fatty acid metabolism, while SREBF2 regulates genes participating in cholesterol biosynthesis and uptake. Our objective was to study the expression of both miR-33a and miR-33b, together with their host SREBF genes, in liver, adipose tissue and muscle to better understand the role of miR-33a/b in the lipid metabolism of pigs. In our study, the expression of miR-33a, miR-33b and SREBF2 in liver, adipose tissue, and muscle was studied in 42 BC1_LD (25% Iberian x 75% Landrace backcross) pigs by RT-qPCR. In addition, the expression of in-silico predicted target genes and fatty acid composition traits were correlated with the miR-33a/b expression. We observed different tissue expression patterns for both miRNAs. In adipose tissue and muscle a high correlation between miR-33a and miR-33b expression was found, whereas a lower correlation was observed in liver. The expression analysis of in-silico predicted target-lipid related genes showed negative correlations between miR-33b and CPT1A expression in liver. Conversely, positive correlations between miR-33a and PPARGC1A and USF1 gene expression in liver were observed. Lastly, positive and negative correlations between miR-33a/b expression and saturated fatty acid (SFA) and polyunsaturated fatty acid (PUFA) content, respectively, were identified. Overall, our results suggested that both miRNAs are differentially regulated and have distinct functions in liver, in contrast to muscle and adipose tissue. Furthermore, the correlations between miR-33a/b expression both with the expression of in-silico predicted target-lipid related genes and with fatty acid composition, opens new avenues to explore the role of miR33a/b in the regulation of lipid metabolism.

Novel Genomic Regions Associated with Intramuscular Fatty Acid Composition in Rabbits

Intramuscular fat (IMF) content and its composition affect the quality of meat. Selection for IMF generated a correlated response on its fatty acid composition. The increase of IMF content is associated with an increase of its saturated (SFA) and monounsaturated (MUFA) fatty acids, and consequently a decrease of polyunsaturated fatty acids (PUFA). We carried out a genome wide association study (GWAS) for IMF composition on two rabbit lines divergently selected for IMF content, using a Bayes B procedure. Association analyses were performed using 475 individuals and 90,235 Single Nucleotide Polymorphisms (SNPs). The main objectives were to identify genomic regions associated with the IMF composition and to generate a list of candidate genes. Genomic regions associated with the intramuscular fatty acid composition were spread across different rabbit chromosomes (OCU). An important region at 34.0-37.9 Mb on OCU1 was associated with C14:0, C16:0, SFA, and C18:2n6, explaining 3.5%, 11.2%, 11.3%, and 3.2% of the genomic variance, respectively. Another relevant genomic region was found to be associated at 46.0-48.9 Mb on OCU18, explaining up to 8% of the genomic variance of MUFA/SFA. The associated regions harbor several genes related to lipid metabolism, such as SCD, PLIN2, and ERLIN1. The main genomic regions associated with the fatty acids were not previously associated with IMF content in rabbits. Nonetheless, MTMR2 is the only gene that was associated with both the IMF content and composition in rabbits. Our study highlighted the polygenic nature of the fatty acids in rabbits and elucidated its genetic background.

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.

Dietary intake of bioactive ingredients impacts liver and adipose tissue transcriptomes in a porcine model of prepubertal early obesity

Global prevalence of obesity has increased to epidemic proportions over the past 40 years, with childhood obesity reaching alarming rates. In this study, we determined changes in liver and adipose tissue transcriptomes of a porcine model for prepubertal early obesity induced by a high-calorie diet and supplemented with bioactive ingredients. A total of 43 nine-weeks-old animals distributed in four pens were fed with four different dietary treatments for 10 weeks: a conventional diet; a western-type diet; and a western-type diet with Bifidobacterium breve and rice hydrolysate, either adding or not omega-3 fatty acids. Animals fed a western-type diet increased body weight and total fat content and exhibited elevated serum concentrations of cholesterol, whereas animals supplemented with bioactive ingredients showed lower body weight gain and tended to accumulate less fat. An RNA-seq experiment was performed with a total of 20 animals (five per group). Differential expression analyses revealed an increase in lipogenesis, cholesterogenesis and inflammatory processes in animals on the western-type diet while the supplementation with bioactive ingredients induced fatty acid oxidation and cholesterol catabolism, and decreased adipogenesis and inflammation. These results reveal molecular mechanisms underlying the beneficial effects of bioactive ingredient supplementation in an obese pig model.

Using SNP Weights Derived From Gene Expression Modules to Improve GWAS Power for Feed Efficiency in Pigs

The "large p small n" problem has posed a significant challenge in the analysis and interpretation of genome-wide association studies (GWAS). The use of prior information to rank genomic regions and perform SNP selection could increase the power of GWAS. In this study, we propose the use of gene expression data from RNA-Seq of multiple tissues as prior information to assign weights to SNP, select SNP based on a weight threshold, and utilize weighted hypothesis testing to conduct a GWAS. RNA-Seq libraries from hypothalamus, duodenum, ileum, and jejunum tissue of 30 pigs with divergent feed efficiency phenotypes were sequenced, and a three-way gene x individual x tissue clustering analysis was performed, using constrained tensor decomposition, to obtain a total of 10 gene expression modules. Loading values from each gene module were used to assign weights to 49,691 commercial SNP markers, and SNP were selected using a weight threshold, resulting in 10 SNP sets ranging in size from 101 to 955 markers. Weighted GWAS for feed intake in 4,200 pigs was performed separately for each of the 10 SNP sets. A total of 36 unique significant SNP associations were identified across the ten gene modules (SNP sets). For comparison, a standard unweighted GWAS using all 49,691 SNP was performed, and only 2 SNP were significant. None of the SNP from the unweighted analysis resided in known QTL related to swine feed efficiency (feed intake, average daily gain, and feed conversion ratio) compared to 29 (80.6%) in the weighted analyses, with 9 SNP residing in feed intake QTL. These results suggest that the heritability of feed intake is driven by many SNP that individually do not attain genome-wide significance in GWAS. Hence, the proposed procedure for prioritizing SNP based on gene expression data across multiple tissues provides a promising approach for improving the power of GWAS.

A novel duplicated insertion/deletion (InDel) of the CPT1a gene and its effects on growth traits in goat

Carnitine palmitoyltransferase 1a (CPT1a) is a rate-limited enzyme in the mitochondrial fatty acid β-oxidation pathway. It acts as a bridge between PPARα and the fatty acid oxidation pathways and is closely related to ruminant growth and development. In this study, one 12 bp InDel polymorphism of the CPT1a gene was identified in 700 goats, and we designated these three genotypes II, ID, and DD. Association analysis showed that the InDel polymorphism was closely associated with trunk index (p = 0.008) and body length index (p = 0.034) in Hainan black goats, and body length (p = 0.010), chest circumference (p = 0.004), chest depth (p = 0.029), and huckle bone width (p = 0.002) in Nubian goats, as well as the chest circumference (p = 0.016) in the Fuqing goat breed. In both kids and adult goats, qRT-PCR results showed that the CPT1a gene was expressed in all tissues, showing the highest mRNA levels in the liver, lung, spleen, and kidney, followed by the adipose tissue and brain. This indicates an association between the InDel of the CPT1a gene and growth traits in selected goat breeds, which may facilitate marker-assisted selection in goat genetics and breeding.

Machine learning applied to transcriptomic data to identify genes associated with feed efficiency in pigs

Background

To date, the molecular mechanisms that underlie residual feed intake (RFI) in pigs are unknown. Results from different genome-wide association studies and gene expression analyses are not always consistent. The aim of this research was to use machine learning to identify genes associated with feed efficiency (FE) using transcriptomic (RNA-Seq) data from pigs that are phenotypically extreme for RFI.

Methods

RFI was computed by considering within-sex regression on mean metabolic body weight, average daily gain, and average backfat gain. RNA-Seq analyses were performed on liver and duodenum tissue from 32 high and 33 low RFI pigs collected at 153 d of age. Machine-learning algorithms were used to predict RFI class based on gene expression levels in liver and duodenum after adjusting for batch effects. Genes were ranked according to their contribution to the classification using the permutation accuracy importance score in an unbiased random forest (RF) algorithm based on conditional inference. Support vector machine, RF, elastic net (ENET) and nearest shrunken centroid algorithms were tested using different subsets of the top rank genes. Nested resampling for hyperparameter tuning was implemented with tenfold cross-validation in the outer and inner loops.

Results

The best classification was obtained with ENET using the expression of 200 genes in liver [area under the receiver operating characteristic curve (AUROC): 0.85; accuracy: 0.78] and 100 genes in duodenum (AUROC: 0.76; accuracy: 0.69). Canonical pathways and candidate genes that were previously reported as associated with FE in several species were identified. The most remarkable pathways and genes identified were NRF2-mediated oxidative stress response and aldosterone signalling in epithelial cells, the DNAJC6, DNAJC1, MAPK8, PRKD3 genes in duodenum, and melatonin degradation II, PPARα/RXRα activation, and GPCR-mediated nutrient sensing in enteroendocrine cells and SMOX, IL4I1, PRKAR2B, CLOCK and CCK genes in liver.

Conclusions

ML algorithms and RNA-Seq expression data were found to provide good performance for classifying pigs into high or low RFI groups. Classification was better with gene expression data from liver than from duodenum. Genes associated with FE in liver and duodenum tissue that can be used as predictive biomarkers for this trait were identified.

Electronic supplementary material

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

Genetic variants and haplotype combination in the bovine CRTC3 affected conformation traits in two Chinese native cattle breeds (Bos Taurus)

CREB-regulated transcription coactivator 3 (CRTC3) plays an extensive role in glucose and lipid metabolism. This study investigated the genetic variation and haplotype combination in CRTC3 and verified their contribution to bovine growth traits. Firstly, investigated the mRNA expression of CRTC3 in adult Qinchuan cattle and evaluated the effects that genetic variation of CRTC3 had on conformation and carcass traits in two Chinese cattle breeds (Qinchuan and Jiaxian). Four SNPs (single nucleotide polymorphisms) were identified including two in introns (SNP1: g.62652 A > G and SNP4: g.91297C > T) and two in exons (SNP2 g.62730C > T and SNP3: g.66478G > C). The association and haplotype combination results showed that there was an association with some growth and carcass traits(P < 0.05). Individuals with haplotype combination H1H1 (-AACCCCTT-) were associated with a conformation of a larger framed animal and an animal that produced a larger loin area. Variations in the CRTC3 genes and the haplotype combination H1H1 may be considered as molecular markers for carcass traits that are associated with more lean meat yield for use in cattle breeding programs in China.

Analysing the Expression of Eight Clock Genes in Five Tissues From Fasting and Fed Sows

In a previous study, we observed that circadian clock genes are differentially expressed in the skeletal muscle of fasting and fed sows. The goal of the current work was to investigate if these genes are also differentially expressed in tissues containing the central (hypothalamus) and peripheral (duodenum, dorsal fat, muscle, and liver) clocks. As animal material, we used 12 sows that fasted 12 h before slaughtering (T0) and 12 sows that were fed ad libitum 7 h prior slaughtering (T2). Tissue samples were collected immediately after slaughter and total RNA was subsequently extracted. The expression of the ARNTL, BHLHE40, CRY2, NPAS2, NR1D1, PER1, PER2, and SIK1 genes was measured by quantitative reverse transcription PCR. The numbers of clock genes showing differential expression before and after feeding varied depending on the tissue i.e., four in dorsal fat and duodenum, six in skeletal muscle, and seven in the liver. In contrast, none of the eight analysed genes displayed a significant differential expression in hypothalamus, the tissue where the central clock resides. This result supports that the differential expression of clock genes in the four tissues mentioned above is probably induced by nutrition and not by the central clock entrained by light. Moreover, we have observed that the NPAS2 and ARNTL genes display positive log2(FC) values in the five tissues under analysis, whilst the CRY2, PER1 (except dorsal fat) and PER2 (except hypothalamus) genes generally show negative log2(FC) values. Such result might be explained by the existence of a negative feedback loop between the ARNTL/NPAS2 and CRY/PER genes. Collectively, these results support that nutrition plays an important role in modulating the timing of porcine peripheral circadian clocks. Such regulation could be essential for coordinating the subsequent metabolic response to nutrient supply.

A genome-wide association study of growth and fatness traits in two pig populations with different genetic backgrounds

Improvement in growth and fatness traits are the main objectives in pig all breeding programs. Tenth rib backfat thickness (10RIBBFT) and days to 100 kg (D100), which are good predictors of carcass lean content and growth rate, respectively, are economically important traits and also main breeding target traits in pigs. To investigate the genetic mechanisms of 10RIBBFT and D100 of pigs, we sampled 1,137 and 888 pigs from 2 Yorkshire populations of American and British origin, respectively, and conducted genome-wide association study (GWAS) through combined analysis and meta-analysis, to identify SNPs associated with 10RIBBFT and D100. A total of 11 and 7 significant SNPs were identified by combined analysis for 10RIBBFT and D100, respectively. And in meta-analysis, 8 and 7 significant SNPs were identified for 10RIBBFT and D100, respectively. Among them, 6 and 5 common significant SNPs in two analysis results were, respectively, identified associated with 10RIBBFT and D100, and correspondingly explained 2.09% and 0.52% of the additive genetic variance of 10RIBBFT and D100. Further bioinformatics analysis revealed 10 genes harboring or close to these common significant SNPs, 5 for 10RIBBFT and 5 for D100. In particular, Gene Ontology analysis highlighted 6 genes, PCK1, ANGPTL3, EEF1A2, TNFAIP8L3, PITX2, and PLA2G12, as promising candidate genes relevant with backfat thickness and growth. PCK1, ANGPTL3, EEF1A2, and TNFAIP8L3 could influence backfat thickness through phospholipid transport, regulation of lipid metabolic process through the glycerophospholipid biosynthesis and metabolism pathway, the metabolism of lipids and lipoproteins pathway. PITX2 has a crucial role in skeletal muscle tissue development and animal organ morphogenesis, and PLA2G12A plays a role in the lipid catabolic and phospholipid catabolic processes, which both are involved in the body weight pathway. All these candidate genes could directly or indirectly influence fat production and growth in Yorkshire pigs. Our findings provide novel insights into the genetic basis of growth and fatness traits in pigs. The candidate genes for D100 and 10RIBBFT are worthy of further investigation.

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.

Integrative approach using liver and duodenum RNA-Seq data identifies candidate genes and pathways associated with feed efficiency in pigs

This study aims identifying candidate genes and pathways associated with feed efficiency (FE) in pigs. Liver and duodenum transcriptomes of 37 gilts showing high and low residual feed intake (RFI) were analysed by RNA-Seq. Gene expression data was explored through differential expression (DE) and weighted gene co-expression network analyses. DE analysis revealed 55 and 112 differentially regulated genes in liver and duodenum tissues, respectively. Clustering genes according to their connectivity resulted in 23 (liver) and 25 (duodenum) modules of genes with a co-expression pattern. Four modules, one in liver (with 444 co-expressed genes) and three in duodenum (gathering 37, 126 and 41 co-expressed genes), were significantly associated with FE indicators. Intra-module analyses revealed tissue-specific candidate genes; 12 of these genes were also identified as DE between individuals with high and low RFI. Pathways enriched by the list of genes showing DE and/or belonging to FE co-expressed modules included response to oxidative stress, inflammation, immune response, lipid metabolism and thermoregulation. Low overlapping between genes identified in duodenum and liver tissues was observed but heat shock proteins were associated to FE in both tissues. Our results suggest tissue-specific rather than common transcriptome regulatory processes associated with FE in pigs.

Association of genetic variants and expression levels of porcine FABP4 and FABP5 genes

The FABP4 and FABP5 genes, coding for fatty acid transport proteins, have long been studied as positional candidate genes for SSC4 QTL affecting fat deposition and composition traits in pigs. Polymorphisms in these genes, FABP4:g.2634_2635insC and FABP5:g.3000T>G, have previously been associated with fatness traits in an Iberian by Landrace cross (IBMAP). The aim of the present work was to evaluate the functional implication of these genetic variants. For this purpose, FABP4 and FABP5 mRNA expression levels in 114 BC1_LD animals (25% Iberian × 75% Landrace) were analyzed using real-time quantitative PCR in backfat and muscle. FABP4 gene expression in backfat, but not in muscle, was associated with FABP4:g.2634_2635insC. In contrast, FABP5:g.3000T>G was not associated with gene expression levels. An expression-based genome-wide association study highlighted the FABP4:g.2634_2635insC polymorphism as the polymorphism most associated with FABP4 gene expression in backfat. Furthermore, other genomic regions associated in trans with the mRNA expression of FABP4 in backfat and FABP5 in muscle were also identified. Finally, two putative transcription binding sites for PPARG and NR4A2 may be affected by the FABP4:g.2634_2635insC polymorphism, modifying FABP4 gene expression. Our results reinforce FABP4 as a candidate gene for fatness traits on SSC4.