Genome-Wide Analysis of lncRNA and mRNA Expression During Differentiation of Abdominal Preadipocytes in the Chicken

Metabolome and Transcriptome Profiling Reveals Age-Associated Variations in Meat Quality and Molecular Mechanisms of Taihe Black-Bone Silky Fowls

To explore the changes in meat quality and molecular mechanisms during the growth and development of Taihe black-bone silky fowl, this study employed liquid chromatography-mass spectrometry (LC-MS/MS) metabolomics to elucidate the dynamic changes of key differential metabolites (DMs) affecting meat quality, indicating that chicken at D120 had higher levels of ω-3 polyunsaturated fatty acids (PUFAs), creatine, anserine, and homocarnosine, while D150 had the most stachydrine and D210 had the most acylcarnitines. Additionally, D120 and D180 had more umami and sweet compounds. Furthermore, key metabolic pathways influenced by age included purine metabolism, the pentose phosphate pathway, nicotinate and nicotinamide metabolism, and taurine and hypotaurine metabolism. Transcriptomic identified differential expression genes (DEGs) are predominantly enriched in focal adhesion, the TGF-β signaling pathway, and the MAPK signaling pathway. Integrated metabolomics and transcriptomics revealed complex regulatory networks of DEGs and DMs in key metabolic pathways. This research enhanced our understanding of the biology of Taihe black-bone silky fowl meat quality, revealing possible biomarkers.

Long Noncoding RNA 6302 Regulates Chicken Preadipocyte Differentiation by Targeting SLC22A16

The excessive deposition of abdominal adipocytes in chickens is detrimental to poultry production. However, the regulatory factors that affect abdominal adipogenesis in chickens are still poorly understood. SLC22A16 is differentially expressed in abdominal preadipocytes and 10-day differentiated adipocytes in chickens, but its role in regulating chicken adipogenesis has not been reported. In this study, the function of SLC22A16 in chicken abdominal preadipocytes was investigated. SLC22A16 is significantly upregulated during abdominal adipocyte differentiation. The overexpression of SLC2A16 upregulated the expression of adipogenic marker genes and proliferation-related genes, and promoted the proliferation of adipocytes and the accumulation of triglycerides. The knockdown of SLC22A16 downregulated the expression of adipogenic marker genes and proliferation-related genes, inhibited the proliferation of adipocytes, and impaired the accumulation of triglycerides in adipocytes. In addition, LNC6302 was differentially expressed in abdominal preadipocytes and mature adipocytes, and was significantly positively correlated with the expression of SLC22A16. Interference with LNC6302 inhibits the expression of adipogenic marker genes and proliferation-related genes. The data supported the notion that LNC6302 promotes the differentiation of chicken abdominal adipocytes by cis-regulating the expression of SLC22A16. This study identified the role of SLC22A16 in the differentiation and proliferation of chicken adipocytes, providing a potential target for improving abdominal adipogenesis in chickens.

Identification of central regulators related to abdominal fat deposition in chickens based on weighted gene co-expression network analysis

Abdominal fat (AF) is one of the most important economic traits in chickens. Excessive AF in chickens will reduce feed utilization efficiency and negatively affect reproductive performance and disease resistance. However, the regulatory network of AF deposition needs to be further elucidated. In the present study, 300 one-day-old female Wannan chickens were reared to 17 wk of age, and 200 Wannan hens were selected to determine the abdominal fat percentage (AFP). Twenty AF tissue samples with the lowest AFP were selected as the low abdominal fat group (L-AFG), and 20 AF tissue samples with the highest AFP were selected as the high abdominal fat group (H-AFG). Eleven samples from L-AFG and 14 samples from H-AFG were selected for RNA-seq and used for weighted gene co-expression network analysis (WGCNA). Among the 25 RNA-seq samples, 5 samples with the lowest and highest AFP values were selected for differential expression gene analysis. Compared with the L-AFG, 225 and 101 genes were upregulated and downregulated in the H-AFG, respectively. A total of 20,503 genes were used to construct the WGCNA, and 44 co-expression gene modules were identified. Among these modules, 3 modules including turquoise, darkorange2, and floralwhite were identified as significantly associated with AFP traits. Furthermore, several genes including acyl-CoA oxidase 1 (ACOX1), stearoyl-CoA desaturase (SCD), aldehyde dehydrogenase 6 family member A1 (ALDH6A1), jun proto-oncogene, AP-1 transcription factor subunit (JUN), and fos proto-oncogene, AP-1 transcription factor subunit (FOS) involved in the "PPAR signaling pathway," "fatty acid metabolism," and "MAPK signaling pathway" were identified as central regulators that contribute to AF deposition. These results provide valuable information for further understanding of the gene expression and regulation of AF traits and contribute to future molecular breeding for AF in chickens.

A review: analysis of technical challenges in cultured meat production and its commercialization

The cultured meat technology has developed rapidly in recent years, but there are still many technical challenges that hinder the large-scale production and commercialization of cultured meat. Firstly, it is necessary to lay the foundation for cultured meat production by obtaining seed cells and maintaining stable cell functions. Next, technologies such as bioreactors are used to expand the scale of cell culture, and three-dimensional culture technologies such as scaffold culture or 3D printing are used to construct the three-dimensional structure of cultured meat. At the same time, it can reduce production costs by developing serum-free medium suitable for cultured meat. Finally, the edible quality of cultured meat is improved by evaluating food safety and sensory flavor, and combining ethical and consumer acceptability issues. Therefore, this review fully demonstrates the current development status and existing technical challenges of the cultured meat production technology with regard to the key points described above, in order to provide research ideas for the industrial production of cultured meat.

Weighted single-step GWAS identified candidate genes associated with carcass traits in a Chinese yellow-feathered chicken population

Carcass traits in broiler chickens are complex traits that are influenced by multiple genes. To gain deeper insights into the genetic mechanisms underlying carcass traits, here we conducted a weighted single-step genome-wide association study (wssGWAS) in a population of Chinese yellow-feathered chicken. The objective was to identify genomic regions and candidate genes associated with carcass weight (CW), eviscerated weight with giblets (EWG), eviscerated weight (EW), breast muscle weight (BMW), drumstick weight (DW), abdominal fat weight (AFW), abdominal fat percentage (AFP), gizzard weight (GW), and intestine length (IL). A total of 1,338 broiler chickens with phenotypic and pedigree information were included in this study. Of these, 435 chickens were genotyped using a 600K single nucleotide polymorphism chip for association analysis. The results indicate that the most significant regions for 9 traits explained 2.38% to 5.09% of the phenotypic variation, from which the region of 194.53 to 194.63Mb on chromosome 1 with the gene RELT and FAM168A identified on it was significantly associated with CW, EWG, EW, BMW, and DW. Meanwhile, the 5 traits have a strong genetic correlation, indicating that the region and the genes can be used for further research. In addition, some candidate genes associated with skeletal muscle development, fat deposition regulation, intestinal repair, and protection were identified. Gene ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analyses suggested that the genes are involved in processes such as vascular development (CD34, FGF7, FGFR3, ITGB1BP1, SEMA5A, LOXL2), bone formation (FGFR3, MATN1, MEF2D, DHRS3, SKI, STC1, HOXB1, HOXB3, TIPARP), and anatomical size regulation (ADD2, AKT1, CFTR, EDN3, FLII, HCLS1, ITGB1BP1, SEMA5A, SHC1, ULK1, DSTN, GSK3B, BORCS8, GRIP2). In conclusion, the integration of phenotype, genotype, and pedigree information without creating pseudo-phenotype will facilitate the genetic improvement of carcass traits in chickens, providing valuable insights into the genetic architecture and potential candidate genes underlying carcass traits, enriching our understanding and contributing to the breeding of high-quality broiler chickens.

RNA-seq analysis reveals the critical role of the novel lncRNA BIANCR in intramuscular adipogenesis through the ERK1/2 signaling pathway

BACKGROUND

Long non-coding RNAs (lncRNAs) regulate numerous biological processes, including adipogenesis. Research on adipogenesis will assist in the treatment of human metabolic diseases and improve meat quality in livestock, such as the content of intramuscular fat (IMF). However, the significance of lncRNAs in intramuscular adipogenesis remains unclear. This research aimed to reveal the lncRNAs transcriptomic profiles in the process of bovine intramuscular adipogenesis and to identify the lncRNAs involved in the adipogenesis of bovine intramuscular adipocytes.

RESULTS

In this research, a landscape of lncRNAs was identified with RNA-seq in bovine intramuscular adipocytes at four adipogenesis stages (0 d, 3 d, 6 d, and 9 d after differentiation). A total of 7035 lncRNAs were detected, including 3396 novel lncRNAs. Based on the results of differential analysis, co-expression analysis, and functional prediction, we focused on the bovine intramuscular adipogenesis-associated long non-coding RNA (BIANCR), a novel lncRNA that may have an important regulatory function. The knockdown of BIANCR inhibited proliferation and promoted apoptosis of intramuscular preadipocytes. Moreover, BIANCR knockdown inhibited intramuscular adipogenesis by regulating the ERK1/2 signaling pathway.

CONCLUSION

This study obtained the landscape of lncRNAs during adipogenesis in bovine intramuscular adipocytes. BIANCR plays a crucial role in adipogenesis through the ERK1/2 signaling pathway. The results are noteworthy for improving beef meat quality, molecular breeding, and metabolic disease research.

Identification of Key Modules and Candidate Genes for Powdery Mildew Resistance of Wheat-Agropyron cristatum Translocation Line WAT-2020-17-6 by WGCNA

As one of the serious diseases of wheat, powdery mildew (Blumeria graminis f. sp. tritici) is a long-term threat to wheat production. Therefore, it is of great significance to explore new powdery mildew-resistant genes for breeding. The wild relative species of wheat provide gene resources for resistance to powdery mildew breeding. Agropyron cristatum (2n = 4x = 28, genomes PPPP) is an important wild relative of wheat, carrying excellent genes for high yield, disease resistance, and stress resistance, which can be used for wheat improvement. To understand the molecular mechanism of powdery mildew resistance in the wheat-A. cristatum translocation line WAT2020-17-6, transcriptome sequencing was performed, and the resistance genes were analyzed by weighted gene co-expression network analysis (WGCNA). In the results, 42,845 differentially expressed genes were identified and divided into 18 modules, of which six modules were highly correlated with powdery mildew resistance. Gene ontology (GO) enrichment analysis showed that the six interested modules related to powdery mildew resistance were significantly enriched in N-methyltransferase activity, autophagy, mRNA splicing via spliceosome, chloroplast envelope, and AMP binding. The candidate hub genes of the interested modules were further identified, and their regulatory relationships were analyzed based on co-expression data. The temporal expression pattern of the 12 hub genes was verified within 96 h after powdery mildew inoculation by RT-PCR assay. In this study, we preliminarily explained the resistance mechanism of the wheat-A. cristatum translocation lines and obtained the hub candidate genes, which laid a foundation in the exploration of resistance genes in A. cristatum for powdery mildew-resistant breeding in wheat.

Effects of Kadsura coccinea L. Fruit Extract on Growth Performance, Meat Quality, Immunity, Antioxidant, Intestinal Morphology and Flora of White-Feathered Broilers

This study aimed to determine whether adding Kadsura coccinea fruit extract to the diet of broilers could replace antibiotics. For this study, 300 one-day-old AA white feathered broilers were divided into five groups (no sex separated), with six repetitions per group (n = 10), as follows: blank control group (basal feed, CK group), positive drug (basal feed + 300 mg/kg aureomycin, PD group), and Kadsura coccinea low-dose, medium-dose, and high-dose groups (basal feed + 100 mg/kg, 200 mg/kg, and 300 mg/kg of Kadsura coccinea fruit extract, LD group, MD group and HD group). The experiment period was divided into early (1−21 days) and late (22−42 days) stage. We found that supplementation with Kadsura coccinea fruit extract in the diet significantly improved the growth performance of broilers (p < 0.05), reduced the feed to meat ratio (p < 0.05), reduced the fat percentage (p < 0.05), while had no significant effect on meat quality (p > 0.05) and Kadsura coccinea fruit extract could promote the development of immune organs to different extents, enhance antioxidant capacity, the contents of SOD and GSH-Px in serum were significantly increased (p < 0.05), improve the ratio of villus height to crypt depth. Finally, Kadsura coccinea fruit extract increased the relative abundance of probiotics and beneficial bacteria (Bacteroidales, NK4A214, Subdoligranulum and Eubacterium hallii) (p < 0.05) and reduced the relative abundance of harmful bacteria (Erysipelatoclostridium) (p < 0.05) in the gut of broilers. Compared with positive drug group, most of the indexes in the medium-dose group were better or had similar effects. We believe that Kadsura coccinea fruit extract can be used as a potential natural antibiotic substitute in livestock and poultry breeding programs.

CHD6 promotes broad nucleosome eviction for transcriptional activation in prostate cancer cells

Despite being a member of the chromodomain helicase DNA-binding protein family, little is known about the exact role of CHD6 in chromatin remodeling or cancer disease. Here we show that CHD6 binds to chromatin to promote broad nucleosome eviction for transcriptional activation of many cancer pathways. By integrating multiple patient cohorts for bioinformatics analysis of over a thousand prostate cancer datasets, we found CHD6 expression elevated in prostate cancer and associated with poor prognosis. Further comprehensive experiments demonstrated that CHD6 regulates oncogenicity of prostate cancer cells and tumor development in a murine xenograft model. ChIP-Seq for CHD6, along with MNase-Seq and RNA-Seq, revealed that CHD6 binds on chromatin to evict nucleosomes from promoters and gene bodies for transcriptional activation of oncogenic pathways. These results demonstrated a key function of CHD6 in evicting nucleosomes from chromatin for transcriptional activation of prostate cancer pathways.

WGCNA Analysis of Important Modules and Hub Genes of Compound Probiotics Regulating Lipid Metabolism in Heat-Stressed Broilers

This study aimed to study compound probiotics’ (Lactobacillus casei, Lactobacillus acidophilus and Bifidobacterium) effects on production performance, lipid metabolism and meat quality in heat-stressed broilers. A total of 400 one-day-old AA broilers were randomly divided into four groups, each containing the same five replicates, with 20 broilers in each replicate. The control (21 °C) and experiment 2 were fed a basic corn−soybean meal diet. Experiment 1 (21 °C) and experiment 3 were fed a basic corn−soybean meal diet with 10 g/kg compound probiotics on days 7 and 28, respectively. The ambient temperature of experiment 2 and experiment 3 was increased to 30−32 °C (9:00−17:00) for 28−42 days, while the temperature for the other time was kept at 21 °C. The results showed that, compared with the control, the production performance and the content of high-density lipoprotein cholesterol in experiment 1 and triglyceride (TG) in experiment 2 increased (p < 0.05). Compared with experiment 2, TG decreased and the production performance increased in experiment 3 (p < 0.05). However, there was no significant change in meat quality indicators. Weighted gene co-expression network analysis (WGCNA) was used to analyze the intramuscular fat, abdominal fat and five blood lipid indicators. We found five related modules. Fatty acid biosynthesis, glycerolipid metabolism, and fat digestion and absorption were the pathways for KEGG enrichment. Additionally, NKX2-1, TAS2R40, PTH, CPB1, SLCO1B3, GNB3 and AQP7 may be the hub genes of compound probiotics regulating lipid metabolism in heat-stressed broilers. In conclusion, this study identified the key genes of compound probiotics regulating lipid metabolism and provided a theoretical basis for the poultry breeding industry to alleviate heat stress.

Resveratrol Inhibits Proliferation and Differentiation of Porcine Preadipocytes by a Novel LincRNA-ROFM/miR-133b/AdipoQ Pathway

Resveratrol (RES) has a wide range of biological and pharmacological activities with various health benefits for humans as a food additive. In animal production, RES has been considered a potential functional feed additive for producing high-quality pork. Long noncoding RNAs (lncRNAs) have emerged as essential regulators of fat metabolism, and phytochemicals can regulate fat metabolism through lncRNA. However, it is unclear whether RES can improve back-fat thickness by regulating lncRNA. In this study, we identified a novel lncRNA, which was named a long intergenic non-protein coding RNA, a regulator of fat metabolism (LincRNA-ROFM), from our previous lncRNA sequencing data. LincRNA-ROFM can inhibit adipocyte proliferation and differentiation. In-depth analyses showed that LincRNA-ROFM acts as a molecular sponge for miR-133b, and adiponectin (AdipoQ) is a direct target of miR-133b in porcine preadipocytes. In addition, the expression of LincRNA-ROFM was positively correlated with AdipoQ. RES can promote the expression of LincRNA-ROFM by PPARα and C/EBPα. Altogether, our research showed that LincRNA-ROFM acts as a ceRNA to sequester miR-133b and is upregulated by RES, leading to heightened AdipoQ expression, and thus decreased adipocyte proliferation and differentiation, which reduces back-fat thickness of pigs. Taken together, the RES/LincRNA-ROFM/miR-133b/AdipoQ regulatory network preliminarily explains the mechanism of action of RES in inhibiting fat deposition, which provides new insight into the downstream mechanism of RES inhibition of fat deposits by regulating the lncRNA.

Genomic Analyses for Selective Signatures and Genes Involved in Hot Adaptation Among Indigenous Chickens From Different Tropical Climate Regions

Climate change, especially weather extremes like extreme cold or extreme hot, is a major challenge for global livestock. One of the animal breeding goals for sustainable livestock production should be to breed animals with excellent climate adaptability. Indigenous livestock and poultry are well adapted to the local climate, and they are good resources to study the genetic footprints and mechanism of the resilience to weather extremes. In order to identify selection signatures and genes that might be involved in hot adaptation in indigenous chickens from different tropical climates, we conducted a genomic analysis of 65 indigenous chickens that inhabit different climates. Several important unique positively selected genes (PSGs) were identified for each local chicken group by the cross-population extended haplotype homozygosity (XP-EHH). These PSGs, verified by composite likelihood ratio, genetic differentiation index, nucleotide diversity, Tajima’s D, and decorrelated composite of multiple signals, are related to nerve regulation, vascular function, immune function, lipid metabolism, kidney development, and function, which are involved in thermoregulation and hot adaptation. However, one common PSG was detected for all three tropical groups of chickens via XP-EHH but was not confirmed by other five types of selective sweep analyses. These results suggest that the hot adaptability of indigenous chickens from different tropical climate regions has evolved in parallel by taking different pathways with different sets of genes. The results from our study have provided reasonable explanations and insights for the rapid adaptation of chickens to diverse tropical climates and provide practical values for poultry breeding.

Transcriptome Sequencing to Identify Important Genes and lncRNAs Regulating Abdominal Fat Deposition in Ducks

Simple Summary

Abdominal fat deposition affects the quality of duck meat and the feed conversion ratio. Here, we performed transcriptome sequencing of the abdominal adipose tissue of ducks with high and low abdominal fat rate by RNA sequencing, exploring the key regulatory genes and lncRNAs related to abdominal fat deposition. As a result, several candidate genes, lncRNAs, and pathways related to abdominal fat deposition in ducks were retrieved. This study lays the foundations for exploring molecular mechanisms underlying the regulation of abdominal fat deposition in ducks, providing a theoretical reference for breeding high-quality meat-producing ducks.

Abstract

Abdominal fat deposition is an important trait in meat-producing ducks. F2 generations of 304 Cherry Valley and Runzhou Crested White ducks were studied to identify genes and lncRNAs affecting abdominal fat deposition. RNA sequencing was used to study abdominal fat tissue of four ducks each with high or low abdominal fat rates. In all, 336 upregulated and 297 downregulated mRNAs, and 95 upregulated and 119 downregulated lncRNAs were identified. Target gene prediction of differentially expressed lncRNAs identified 602 genes that were further subjected to Gene Ontology and KEGG pathway analysis. The target genes were enriched in pathways associated with fat synthesis and metabolism and participated in biological processes, including Linoleic acid metabolism, lipid storage, and fat cell differentiation, indicating that these lncRNAs play an important role in abdominal fat deposition. PPAPA, FOXO3, FASN, PNPLA2, FKBP5, TCF7L2, BMP2, FGF2, LIFR, ZBTB16, SIRT, GYG2, NCOR1, and NR3C1 were involved in the regulation of abdominal fat deposition. PNPLA2, TCF7L2, FGF2, LIFR, BMP2, FKBP5, GYG2, and ZBTB16 were regulated by the lncRNAs TCONS_00038080, TCONS_0033547, TCONS_00066773, XR_001190174.3, XR_003492471.1, XR_003493494.1, XR_001192142.3, XR_002405656.2, XR_002401822.2, XR_003497063.1, and so on. This study lays foundations for exploring molecular mechanisms underlying the regulation of abdominal fat deposition in ducks and provides a theoretical basis for breeding high-quality meat-producing ducks.

Identification of a Long Noncoding RNA (lncPRDM16) Inhibiting Preadipocyte Proliferation in the Chicken

Long noncoding RNAs are vital to a variety of biological and physiological processes through multiple modes of functional interaction with DNA, RNA, and proteins. In chickens, numerous lncRNAs were discovered to be important to growth or disease progression. However, the detailed molecular function and role of lncRNAs remain less explored. Here, we performed lncRNA sequencing on abdominal adipose tissues from broiler lines divergently selected for abdominal fat content, and significantly differentially expressed lncRNAs were found, including lncPRDM16, a divergently transcribed and conserved lncRNA near PRDM16. Full lengths of two transcripts of lncPRDM16 were obtained, and their genomic structures were compared. Expression dynamics of lncPRDM16 in different tissues and during preadipocyte proliferation and differentiation were profiled. Moreover, a 250-nucleotide sequence at 5'-end was found to be inevitable to the function of lncPRDM16 in inhibiting preadipocyte proliferation and regulating the promoter activities of both lncPRDM16 and PRDM16. Taken together, we identified the 5'-end functional elements of lncPRDM16 and their potential importance in inhibiting preadipocyte proliferation. Our findings provide the foundation for further exploration of lncPRDM16 function and potential improvement of chicken muscle quality.

Differentially Expressed lncRNAs Related to the Development of Abdominal Fat in Gushi Chickens and Their Interaction Regulatory Network

Chickens are one of the most important sources of meat worldwide, and the growth status of abdominal fat is closely related to production efficiency. Long noncoding RNAs (lncRNAs) play an important role in lipid metabolism and deposition regulation. However, research on the expression profile of lncRNAs related to the development of abdominal fat in chickens after hatching and their interaction regulatory networks is still lacking. To characterize the lncRNA expression profile during the development of chicken abdominal fat, abdominal adipose tissues from 6-, 14-, 22-, and 30-week-old Chinese Gushi chickens were herein used to construct 12 cDNA libraries, and a total of 3,827 new lncRNAs and 5,466 previously annotated lncRNAs were revealed. At the same time, based on the comparative analysis of five combinations, 276 differentially expressed lncRNAs (DE-lncRNAs) were screened. Functional enrichment analysis showed that the predicted target genes of these DE-lncRNAs were significantly enriched in pathways related to the posttranscriptional regulation of gene expression, negative regulation of cell proliferation, cell adhesion and other biological processes, glycosphingolipid biosynthesis, PPAR signaling, fatty acid degradation, fatty acid synthesis and others. In addition, association analysis of the lncRNA transcriptome profile was performed, and DE-lncRNA-related lncRNA-mRNA, lncRNA-miRNA and lncRNA-miRNA-mRNA interaction regulatory networks were constructed. The results showed that DE-lncRNA formed a complex network with PPAR pathway components, including PPARD, ACOX1, ADIPOQ, CPT1A, FABP5, ASBG2, LPL, PLIN2 and related miRNAs, including mir-200b-3p, mir-130b-3p, mir-215-5p, mir-122-5p, mir-223 and mir-125b-5p, and played an important regulatory role in biological processes such as lipid metabolism, adipocyte proliferation and differentiation. This study described the dynamic expression profile of lncRNAs in the abdominal fat of Gushi chickens for the first time and constructed the DE-lncRNA interaction regulatory network. The results expand the number of known lncRNAs in chicken abdominal fat and provide valuable resources for further elucidating the posttranscriptional regulatory mechanism of chicken abdominal fat development or deposition.

The emerging roles and mechanisms of exosomal non-coding RNAs in the mutual regulation between adipose tissue and other related tissues in obesity and metabolic diseases

Exosomes (EXs) are the major types of extracellular vesicles (EVs) of 30-100 nm diameter that can be secreted by most cells to the extracellular environment. EXs transport endogenous cargoes (proteins, lipids, RNAs, etc.) to target cells and thereby triggers the release of these bioactive components, which then play important roles in regulating numerous biological processes under both physiological and pathological conditions. Throughout the studies in recent years, growing evidences have shown that EXs-derived non-coding RNAs (EXs-ncRNAs) are emerging as key players in cell-to-cell communication between adipose tissue and other related tissues in obesity and metabolic diseases. In this review, we will summarize the recent findings about EXs-ncRNAs, especially focus on the following aspects: 1) the biogenesis of EXs and emerging roles of EXs-ncRNAs, 2) the role of EXs-ncRNAs (EXs-miRNAs, EXs-lncRNAs, EXs-circRNAs, etc.) that were secreted by adipose-related tissues in promoting the differentiation of preadipocytes into mature and fully functional adipocytes, and 3) the crosstalk between the adipose tissue derived EXs-ncRNAs and the development of insulin resistance, obesity and various cancers. This review aims to reveal the emerging roles and mechanisms of EXs-ncRNAs in the mutual regulation of adipose tissue and its related tissues in obesity and metabolic diseases, so as to provide references for elucidating the etiology of obesity and related metabolic diseases and screening novel therapeutic targets.

Effects of Feeding Fermented Mulberry Leaf Powder on Growth Performance, Slaughter Performance, and Meat Quality in Chicken Broilers

Simple Summary

Mulberry leaf is widely used in ruminants feeding, such as sheep, beef cattle, and dairy calves. Due to the high content of crude fiber in mature mulberry leaves and branches and the presence of anti-nutritional factors such as tannin, excessive addition will affect the production performance and health of livestock and poultry, and limit its large-scale application in animal production to a certain extent. The disadvantages of woody plants can be improved by microbial fermentation, which can reduce the content of anti-nutritional factors, and increase the content of peptides and amino acids, probiotics, and bioactive components. In this study, Lactobacillus, Saccharomycetes, and Bacillus subtilis were used to make mixed strains to ferment mulberry leaf powder, and different proportions were added to the diet of yellow feathered chicken broilers. The results showed that the addition of fermented mulberry leaf in the diet could improve the digestion and absorption of nutrients, and then improve its growth performance, and increase the contents of inosine monophosphate (IMP), total amino acids, essential amino acids, and delicious amino acids in breast and thigh muscle, and improved polyunsaturated fatty acids and essential fatty acids in breast muscle; this also has a positive effect on improving meat quality.

Abstract

This study was conducted to investigate the effects of feeding fermented mulberry leaf powder (FMLP) on growth performance, slaughter performance, and meat quality of broilers. A total of 360 1-day-old chickens were randomly divided into 5 groups. The control group was fed basal diet (CON), 3% FMLP, 6% FMLP, 9% FMLP, and 3% unfermented mulberry leaf powder. The (MLP) group was fed basal diet supplemented with 3%, 6%, 9% fermented mulberry leaf powder, and 3% MLP, respectively. The experiment lasted for 56 days, with 1–28 days as the starter phase and 29–56 days as the grower phase. The results on the growth performance showed that diets supplemented with 3% FMLP significantly increased the ratio of villus height to crypt depth in the duodenum, jejunum, and ileum of broilers, enhanced the activity of intestinal amylase and digestibility of dry matter and crude protein, improved the average daily gain (ADG), and decreased the feed to gain ratio (F/G) (p < 0.05). Compared with the control group diet, the 3% FMLP group diet significantly increased the breast muscle yield (p < 0.05), reduced the abdominal fat ratio (0.1 < p < 0.05), and improved the slaughter performance of broilers. The 3% MLP group diet increased the shear force of breast muscle (p < 0.05) and thigh muscle of broilers compared to the control group, and adding FMLP could reverse the above results. Additionally, relative to the control group, FMLP supplementation improved the contents of inosine monophosphate (IMP), total amino acids (TAA), essential amino acids (EAA), and delicious amino acids (DAA) in breast and thigh muscle, and improved polyunsaturated fatty acids (PUFA) and essential fatty acids (EFA) in breast muscle; the 6% and 9% FMLP groups showed preferably such effects (p < 0.05). In conclusion, dietary supplementation of FMLP can improve the digestion and absorption of nutrients, and then improve the growth performance of broilers; it also has a positive effect on improving slaughter performance and meat quality.

LncRNAs in domesticated animals: from dog to livestock species

Animal genomes are pervasively transcribed into multiple RNA molecules, of which many will not be translated into proteins. One major component of this transcribed non-coding genome is the long non-coding RNAs (lncRNAs), which are defined as transcripts longer than 200 nucleotides with low coding-potential capabilities. Domestic animals constitute a unique resource for studying the genetic and epigenetic basis of phenotypic variations involving protein-coding and non-coding RNAs, such as lncRNAs. This review presents the current knowledge regarding transcriptome-based catalogues of lncRNAs in major domesticated animals (pets and livestock species), covering a broad phylogenetic scale (from dogs to chicken), and in comparison with human and mouse lncRNA catalogues. Furthermore, we describe different methods to extract known or discover novel lncRNAs and explore comparative genomics approaches to strengthen the annotation of lncRNAs. We then detail different strategies contributing to a better understanding of lncRNA functions, from genetic studies such as GWAS to molecular biology experiments and give some case examples in domestic animals. Finally, we discuss the limitations of current lncRNA annotations and suggest research directions to improve them and their functional characterisation.

Construction of miRNA-mRNA network in the differentiation of chicken preadipocytes

1. MicroRNAs (miRNAs) play key roles in regulating lipid metabolism, adipogenesis and fat deposition in chicken. To date, there are only a few miRNAs that had been confirmed to be involved in chicken adipogenesis. The detailed mechanisms by which miRNAs regulate chicken adipogenesis remain largely unknown.

2. To identify candidate miRNAs involved in chicken preadipocyte differentiation and explore potential mechanisms behind their functions, the following study analysed and identified miRNA and mRNA expression levels in undifferentiated and differentiated preadipocytes. Hub miRNA-mRNA interactions were identified, and the degree of connectivity of DE miRNAs in the network was established.

3. A total of 145 DE miRNAs and 660 DE mRNAs were identified between undifferentiated and differentiated preadipocytes. An miRNA-mRNA network was constructed, including 29 DE miRNAs and 155 DE mRNAs, forming 470 miRNA-mRNA interactions. Functional enrichment analysis showed that DE mRNAs in the network were significantly enriched in 712 biological processes and 13 KEGG pathways. Based on the connectivity degree, five DE miRNAs with higher degrees miR-195-x, gga-miR-200a-3p, gga-miR-135a-5p, novel-m0067-5p and novel-m0270-5p were identified as hub miRNAs. Fifty-eight DE mRNAs interacted with these five hub miRNAs and formed 70 miRNA-mRNA interactions.

4. This study constructed a miRNA-mRNA network associated with chicken preadipocyte differentiation and identified five hub miRNAs in the network. The findings identified the number of chicken adipogenic miRNAs and laid the foundation for elucidating the miRNA-mediated regulatory mechanism in chicken adipogenesis.

Novel 61-bp Indel of RIN2 Is Associated With Fat and Hatching Weight Traits in Chickens

The Ras and Rab interactor 2 (RIN2) gene, which encodes RAS and Rab interacting protein 2, can interact with GTP-bound Rab5 and participate in early endocytosis. This study found a 61-bp insertion/deletion (indel) in the RIN2 intron region, and 3 genotypes II, ID, and DD were observed. Genotype analysis of mutation sites was performed on 665 individuals from F₂ population and 8 chicken breeds. It was found that the indel existed in each breed and that yellow feathered chickens were mainly of the DD genotype. Correlation analysis of growth and carcass traits in the F₂ population of Xinghua and White Recessive Rock chickens showed that the 61-bp indel was significantly correlated with abdominal fat weight, abdominal fat rate, fat width, and hatching weight (P < 0.05). RIN2 mRNA was expressed in all the tested tissues, and its expression in abdominal fat was higher than that in other tissues. In addition, the expression of the RIN2 mRNA in the abdominal fat of the DD genotype was significantly higher than that of the II genotype (P < 0.05). The transcriptional activity results showed that the luciferase activity of the pGL3-DD vector was significantly higher than that of the pGL3-II vector (P < 0.01). Moreover, the results indicate that the polymorphisms in transcription factor binding sites (TFBSs) of 61-bp indel may affect the transcriptional activity of RIN2, and thus alter fat traits in chicken. The results of this study showed that the 61-bp indel was closely related to abdominal fat-related and hatching weight traits of chickens, which may have reference value for molecular marker-assisted selection of chickens.

Whole Transcriptome Analysis Reveals a Potential Regulatory Mechanism of LncRNA-FNIP2/miR-24-3p/FNIP2 Axis in Chicken Adipogenesis

Lipid biosynthesis is a complex process, which is regulated by multiple factors including lncRNA. However, the role of lncRNA in chicken abdominal fat accumulation is still unclear. In this research, we collected liver tissues from six high abdominal fat rate Sanhuang broilers and six low abdominal fat rate Sanhuang broilers to perform lncRNA sequencing and small RNA sequencing. A total of 2,265 lncRNAs, 245 miRNAs, and 5,315 mRNAs were differently expressed. Among of them, 1,136 differently expressed genes were enriched in the metabolic process. A total of 36 differently expressed genes, which were considered as differently expressed lncRNAs' targets, were enriched in the metabolic process. In addition, we also found out that eight differently expressed miRNAs could target 19 differently expressed genes. FNIP2 and PEX5L were shared in a cis-regulatory network and a differently expressed miRNA target relationship network. LncRNA-FNIP2/miR-24-3p/FNIP2 axis was considered as a potential candidate that may participate in lipid synthesis. Experimentally, the objective reality of lncRNA-FNIP2/miR-24-3p/FNIP2 axis was clarified and the regulation effect of lncRNA-FNIP2/miR-24-3p/FNIP2 axis on synthesis was validated. In brief, our study reveals a potential novel regulatory mechanism that lncRNA-FNIP2/miR-24-3p/FNIP2 axis was considered as being involved in lipid synthesis during chicken adipogenesis in liver.

Integrative Analyses of mRNA Expression Profile Reveal SOCS2 and CISH Play Important Roles in GHR Mutation-Induced Excessive Abdominal Fat Deposition in the Sex-Linked Dwarf Chicken

Sex-linked dwarf (SLD) chicken, which is caused by a recessive mutation of the growth hormone receptor (GHR), has been widely used in the Chinese broiler industry. However, it has been found that the SLD chicken has more abdominal fat deposition than normal chicken. Excessive fat deposition not only reduced the carcass quality of the broilers but also reduced the immunity of broilers to diseases. To find out the key genes and the precise regulatory pathways that were involved in the GHR mutation-induced excessive fat deposition, we used high-fat diet (HFD) and normal diet to feed the SLD chicken and normal chicken and analyzed the differentially expressed genes (DEGs) among the four groups. Results showed that the SLD chicken had more abdominal fat deposition and larger adipocytes size than normal chicken and HFD can promote abdominal fat deposition and induce adipocyte hypertrophy. RNA sequencing results of the livers and abdominal fats from the above chickens revealed that many DEGs between the SLD and normal chickens were enriched in fat metabolic pathways, such as peroxisome proliferator-activated receptor (PPAR) signaling, extracellular matrix (ECM)-receptor pathway, and fatty acid metabolism. Importantly, by constructing and analyzing the GHR-downstream regulatory network, we found that suppressor of cytokine signaling 2 (SOCS2) and cytokine-inducible SH2-containing protein (CISH) may involve in the GHR mutation-induced abdominal fat deposition in chicken. The ectopic expression of SOCS2 and CISH in liver-related cell line leghorn strain M chicken hepatoma (LMH) cell and immortalized chicken preadipocytes (ICP) revealed that these two genes can regulate fatty acid metabolism, adipocyte differentiation, and lipid droplet accumulation. Notably, overexpression of SOCS2 and CISH can rescue the hyperactive lipid metabolism and excessive lipid droplet accumulation of primary liver cell and preadipocytes that were isolated from the SLD chicken. This study found some genes and pathways involved in abdominal fat deposition of the SLD chicken and reveals that SOCS2 and CISH are two key genes involved in the GHR mutation-induced excessive fat deposition of the SLD chicken.

Effect of Vitamin E Supplementation on Deposition and Gene Expression Profiling of Abdominal Fat in Broiler Chickens

The aim of this study was to study the regulation of abdominal fat deposition by DL-α-tocopherol acetate (vitamin E) in broilers. Diets supplemented with 50 IU vitamin E significantly diminished abdominal fat deposition in broilers at day 35. Transcriptome sequencing results for abdominal fat tissues of the control (FC) and 50 IU vitamin E-supplemented (FT) groups identified 602 differentially expressed genes (DEGs), which were enriched in cellular process, cell and cell part, and binding Gene Ontology terms. Pathway functional analysis revealed that the DEGs were enriched in 42 metabolic pathways. Notably, the most enriched pathway, fatty acid biosynthesis, was found to play a key role in lipid metabolism. Further, the key regulators of lipid metabolism, including fatty acid synthase, acetyl-CoA carboxylase alpha, and acyl-CoA synthetase long-chain family member 1, demonstrated decreased expression following vitamin E supplementation. Herein, we have identified pathways and genes regulated by vitamin E, thereby providing novel insights into the nutrients regulating abdominal fat deposition in broilers.

Gene co-expression network analysis to identify critical modules and candidate genes of drought-resistance in wheat

AIM

To establish a gene co-expression network for identifying principal modules and hub genes that are associated with drought resistance mechanisms, analyzing their mechanisms, and exploring candidate genes.

METHODS AND FINDINGS

42 data sets including PRJNA380841 and PRJNA369686 were used to construct the co-expression network through weighted gene co-expression network analysis (WGCNA). A total of 1,896,897,901 (284.30 Gb) clean reads and 35,021 differentially expressed genes (DEGs) were obtained from 42 samples. Functional enrichment analysis indicated that photosynthesis, DNA replication, glycolysis/gluconeogenesis, starch and sucrose metabolism, arginine and proline metabolism, and cell cycle were significantly influenced by drought stress. Furthermore, the DEGs with similar expression patterns, detected by K-means clustering, were grouped into 29 clusters. Genes involved in the modules, such as dark turquoise, yellow, and brown, were found to be appreciably linked with drought resistance. Twelve central, greatly correlated genes in stage-specific modules were subsequently confirmed and validated at the transcription levels, including TraesCS7D01G417600.1 (PP2C), TraesCS5B01G565300.1 (ERF), TraesCS4A01G068200.1 (HSP), TraesCS2D01G033200.1 (HSP90), TraesCS6B01G425300.1 (RBD), TraesCS7A01G499200.1 (P450), TraesCS4A01G118400.1 (MYB), TraesCS2B01G415500.1 (STK), TraesCS1A01G129300.1 (MYB), TraesCS2D01G326900.1 (ALDH), TraesCS3D01G227400.1 (WRKY), and TraesCS3B01G144800.1 (GT).

CONCLUSIONS

Analyzing the response of wheat to drought stress during different growth stages, we have detected three modules and 12 hub genes that are associated with drought resistance mechanisms, and five of those genes are newly identified for drought resistance. The references provided by these modules will promote the understanding of the drought-resistance mechanism. In addition, the candidate genes can be used as a basis of transgenic or molecular marker-assisted selection for improving the drought resistance and increasing the yields of wheat.

Review: Long non-coding RNA in livestock

Less than 2% of mammalian genomes code for proteins, but 'the majority of its bases can be found in primary transcripts' - a phenomenon termed the pervasive transcription, which was first reported in 2007. Even though most of the transcripts do not code for proteins, they play a variety of biological functions, with regulation of gene expression appearing as the most common one. Those transcripts are divided into two groups based on their length: small non-coding RNAs, which are maximally 200 bp long, and long non-coding RNAs (lncRNAs), which are longer than 200 nucleotides. The advances in next-generation sequencing methods provided a new possibility of investigating the full set of RNA molecules in the cell. In this review, we summarized the current state of knowledge on lncRNAs in three major livestock species - Sus scrofa, Bos taurus and Gallus gallus, based on the literature and the content of biological databases. In the NONCODE database, the largest number of identified lncRNA transcripts is available for pigs, but cattle have the largest number of lncRNA genes. Poultry is represented by less than a half of records. Genomic annotation of lncRNAs showed that the majority of them are assigned to introns (pig, poultry) or intergenic (cattle). The comparison with well-annotated human and mouse genomes indicates that such annotation is a result of lack of proper lncRNA annotation data. Since lncRNAs play an important role in genomic studies, their characterization in farm animals' genomes is critical in bridging the gap between genotype and phenotype.

MiRNAs and mRNAs Analysis during Abdominal Preadipocyte Differentiation in Chickens

Simple Summary

We sequenced the miRNAs and mRNAs of preabdominal fat cells and differentiated adipocytes, and target genes of miRNA combined with mRNA transcriptome data jointly. We found that the MAPK signal pathway, insulin signal pathway, fatty acid metabolism, ECM( extracellular matrix)–receptor interaction, and other signal pathways were involved in the differentiation of preabdominal fat cells. In addition, we found that some miRNAs–mRNAs combinations were strongly related to the differentiation of fat cells (miR-214−ACSBG2, NFKB2, CAMK2A, ACLY, CCND3, PLK3, ITGB2; miR-148a-5p−ROCK2; miR-10a-5p−ELOVL5; miR-146b-5p−LAMA4; miR-6615-5p−FLNB; miR-1774−COL6A1). Our findings provide important resources for the study of adipocyte differentiation.

Abstract

The excessive deposition of abdominal fat has become an important factor in restricting the production efficiency of chickens, so reducing abdominal fat deposition is important for improving growth rate. It has been proven that miRNAs play an important role in regulating many physiological processes of organisms. In this study, we constructed a model of adipogenesis by isolating preadipocytes (Ab-Pre) derived from abdominal adipose tissue and differentiated adipocytes (Ab-Ad) in vitro. Deep sequencing of miRNAs and mRNAs expressed in Ab-Pre and Ab-Ad groups was conducted to explore the effect of miRNAs and mRNAs on fat deposition. We identified 80 differentially expressed miRNAs (DEMs) candidates, 58 of which were up-regulated and 22 down-regulated. Furthermore, six miRNAs and six mRNAs were verified by qRT-PCR, and the results showed that the expression of the DEMs and differentially expressed genes (DEGs) in the two groups was consistent with our sequencing results. When target genes of miRNA were combined with mRNA transcriptome data, a total of 891 intersection genes were obtained, we predicted the signal pathways of cross genes enrichment to the MAPK signal pathway, insulin signal pathway, fatty acid metabolism, and ECM–receptor interaction. Meanwhile, we constructed miRNA and negatively correlated mRNA target networks, including 12 miRNA–mRNAs pairs, which showed a strong association with the abdominal adipocyte differentiation (miR-214−ACSBG2, NFKB2, CAMK2A, ACLY, CCND3, PLK3, ITGB2; miR-148a-5p−ROCK2; miR-10a-5p−ELOVL5; miR-146b-5p−LAMA4; miR-6615-5p−FLNB; miR-1774−COL6A1). Overall, these findings provide a background for further research on lipid metabolism. Thus, we can better understand the molecular genetic mechanism of chicken abdominal fat deposition.

Outcomes of a low birth weight phenotype on piglet gut microbial composition and intestinal transcriptomic profile

Decades of selection for increased litter size has caused a proportion of sows to consistently produce low birth weight (LBW) litters resulting in economic loss for producers due to reduced piglet survivability and growth. We hypothesized that piglets from LBW litters would have altered gut microbial composition, intestinal architecture, and intestinal transcriptomic profiles compared with piglets from high birth weight (HBW) litters. Sows were designated LBW (n = 45) or HBW (n = 46) based on litter birth weights of three successive parities. LBW piglets were 22% lighter (P < 0.001) at birth; however, no longer differed (P > 0.05) in weight at weaning compared with HBW piglets. LBW piglets had reduced (P < 0.05) fecal microbial diversity with a 114% increase in fecal Enterobacteriaceae (P < 0.05), as well as reduced (P < 0.05) abundance of cecal Roseburia and Faecalibacterium, fiber-degrading butyrate producers. Several genes associated with metabolic (PER2, CES1, KLHL38, and HK2) and immune pathways (IL-1B, IRF8, and TNIP3) were differentially expressed, suggesting altered metabolic and immune function in LBW piglets. In conclusion, LBW piglets had potentially unfavorable shifts in microbial structure in comparison to HBW piglets accompanied with alterations in metabolic and immune gene expression. Results indicate some biological consequences linking LBW phenotype to changes in production efficiency later in life.

Genetic regulatory networks for salt-alkali stress in Gossypium hirsutum with differing morphological characteristics

Background

Cotton grows in altering environments that are often unfavorable or stressful for its growth and development. Consequently, the plant must cope with abiotic stresses such as soil salinity, drought, and excessive temperatures. Alkali-salt stress response remains a cumbersome biological process and is regulated via a multifaceted transcriptional regulatory network in cotton.

Results

To discover the molecular mechanisms of alkali-salt stress response in cotton, a comprehensive transcriptome analysis was carried out after alkali-salt stress treatment in three accessions of Gossypium hirsutum with contrasting phenotype. Expression level analysis proved that alkali-salt stress response presented significant stage-specific and tissue-specific. GO enrichment analysis typically suggested that signal transduction process involved in salt-alkali stress response at SS3 and SS12 stages in leaf; carbohydrate metabolic process and oxidation-reduction process involved in SS48 stages in leaf; the oxidation-reduction process involved at all three phases in the root. The Co-expression analysis suggested a potential GhSOS3/GhCBL10-SOS2 network was involved in salt-alkali stress response. Furthermore, Salt-alkali sensitivity was increased in GhSOS3 and GhCBL10 Virus-induced Gene Silencing (VIGS) plants.

Conclusion

The findings may facilitate to elucidate the underlying mechanisms of alkali-salt stress response and provide an available resource to scrutinize the role of candidate genes and signaling pathway governing alkali-salt stress response.

Screening and Identification of Differentially Expressed and Adipose Growth-Related Protein-Coding Genes During the Deposition of Perirenal Adipose Tissue in Rabbits

BACKGROUND

Rabbit is a good model for genetic and medical studies in other livestock species. The rabbit shows low adipose tissue deposition, and the phenomena indicates that there is some specificity of adipose deposition during the rabbit growth. However, little is known about genes that regulate the growth of adipose tissue in rabbits.

MATERIALS AND METHODS

Deep RNA-seq and comprehensive bioinformatics analyses were used to characterize the genes of rabbit visceral adipose tissue (VAT) at 35, 85 and 120 days after birth. Differentially expressed genes (DEGs) were identified at the three growth stages by DESeq. To explore the function of the candidate genes, Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were performed. Six DEGs were randomly selected, and their expression profiles were validated by q-PCR.

RESULTS

A total of 20,303 known transcripts and 99,199 new transcripts from 8 RNA sequencing libraries were identified, and 34 differentially expressed genes (DEGs) were screened. GO enrichment and KEGG pathway analyses revealed that the DEGs were mainly involved in lipid metabolism regulation including acylglycerol metabolic process and mobilization, and decomposition of lipids to generate ATP in adipocytes and fatty acid metabolism, included LOC100342322 and LOC100342572. In addition, 133 protein-coding genes that play a role in adipose growth and development were screened, including acyl-CoA synthetase long-chain family member 5 (ACSL5) and fatty acid-binding protein 2 (FABP2). The validation results of six DEGs by q-PCR showed similar trends with the results of RNA-seq.

CONCLUSION

In summary, this study provides the first report of the coding genes profiles of rabbit adipose tissue during different growth stages. These data allow for the identification of candidate genes for subsequent studies on rabbit genetics and regulation of adipose cells, and provide an animal model for studying obesity in humans.

Combining information from genome-wide association and multi-tissue gene expression studies to elucidate factors underlying genetic variation for residual feed intake in Australian Angus cattle

BACKGROUND

Genome-wide association studies (GWAS) are extensively used to identify single nucleotide polymorphisms (SNP) underlying the genetic variation of complex traits. However, much uncertainly often still exists about the causal variants and genes at quantitative trait loci (QTL). The aim of this study was to identify QTL associated with residual feed intake (RFI) and genes in these regions whose expression is also associated with this trait. Angus cattle (2190 steers) with RFI records were genotyped and imputed to high density arrays (770 K) and used for a GWAS approach to identify QTL associated with RFI. RNA sequences from 126 Angus divergently selected for RFI were analyzed to identify the genes whose expression was significantly associated this trait with special attention to those genes residing in the QTL regions.

RESULTS

The heritability for RFI estimated for this Angus population was 0.3. In a GWAS, we identified 78 SNPs associated with RFI on six QTL (on BTA1, BTA6, BTA14, BTA17, BTA20 and BTA26). The most significant SNP was found on chromosome BTA20 (rs42662073) and explained 4% of the genetic variance. The minor allele frequencies of significant SNPs ranged from 0.05 to 0.49. All regions, except on BTA17, showed a significant dominance effect. In 1 Mb windows surrounding the six significant QTL, we found 149 genes from which OAS2, STC2, SHOX, XKR4, and SGMS1 were the closest to the most significant QTL on BTA17, BTA20, BTA1, BTA14, and BTA26, respectively. In a 2 Mb windows around the six significant QTL, we identified 15 genes whose expression was significantly associated with RFI: BTA20) NEURL1B and CPEB4; BTA17) RITA1, CCDC42B, OAS2, RPL6, and ERP29; BTA26) A1CF, SGMS1, PAPSS2, and PTEN; BTA1) MFSD1 and RARRES1; BTA14) ATP6V1H and MRPL15.

CONCLUSIONS

Our results showed six QTL regions associated with RFI in a beef Angus population where five of these QTL contained genes that have expression associated with this trait. Therefore, here we show that integrating information from gene expression and GWAS studies can help to better understand the genetic mechanisms that determine variation in complex traits.

Genome-wide identification and characterization of long non-coding RNAs during differentiation of visceral preadipocytes in rabbit

Emerging evidence suggests that long non-coding RNAs (lncRNAs) are critical regulators of diverse biological processes, including adipogenesis. Despite being considered an ideal animal model for studying adipogenesis, little is known about the roles of lncRNAs in the regulation of rabbit preadipocyte differentiation. In the present study, visceral preadipocytes isolated from newborn rabbits were cultured in vitro and induced for differentiation, and global lncRNA expression profiles of adipocytes collected at days 0, 3, and 9 of differentiation were analyzed by RNA-seq. A total of 2066 lncRNAs were identified from nine RNA-seq libraries. Compared to protein-coding transcripts, lncRNA transcripts exhibited characteristics of a longer length and lower expression level. Furthermore, 486 and 357 differentially expressed (DE) lncRNAs were identified when comparing day 3 vs. day 0 and day 9 vs. day 3, respectively. Target genes of DE lncRNAs were predicted by the cis-regulating approach. Prediction of functions revealed that DE lncRNAs when comparing day 3 vs. day 0 were involved in gene ontology (GO) terms of developmental growth, growth, developmental cell growth, and stem cell proliferation, and involved in Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways of PI3K-Akt signaling pathway, fatty acid biosynthesis, and the insulin signaling pathway. The DE lncRNAs when comparing day 9 vs. day 3 were involved in GO terms that associated with epigenetic modification and were involved in the KEGG pathway of cAMP signaling pathway. This study provides further insight into the regulatory function of lncRNAs in rabbit visceral adipose and facilitates a better understanding of different stages of preadipocyte differentiation.

High-Throughput RNA Sequencing Reveals NDUFC2-AS lncRNA Promotes Adipogenic Differentiation in Chinese Buffalo (Bubalus bubalis L)

The buffalo (Bubalus bubalis L.) is prevalent in China and the increasing demand for meat production has changed its role from being a beast of burden to a meat source. The low fat deposition level has become one of the main barriers for its use in meat production. It is urgent to reveal factors involved in fat deposition in buffalo. This study performed RNA sequencing to investigate both long noncoding RNAs (lncRNAs) and mRNAs of adipose tissues in young and adult buffalos. A total of 124 lncRNAs and 2008 mRNAs showed differential expression patterns between young and adult samples. Coexpression analysis and functional enrichment revealed 585 mRNA–lncRNA pairs with potential function in fat deposition. After validation by qRT-PCR, we focused on a lncRNA transcribed from the ubiquinone oxidoreductase subunit C2 (NDUFC2) antisense (AS) strand which showed high correlation with thyroid hormone responsive protein (THRSP). NDUFC2-AS lncRNA is highly expressed in adipose tissue and maturation adipocytes and mainly exists in the nucleus. Functional assays demonstrated that NDUFC2-AS lncRNA promotes adipogenic differentiation by upregulating the expression levels of THRSP and CCAAT enhancer binding protein alpha (C/EBPα) in buffalo. These results indicate that NDUFC2-AS lncRNA promotes fat deposition in buffalo.

Dynamics of transcriptome changes during subcutaneous preadipocyte differentiation in ducks

BACKGROUND

Pekin duck is an important animal model for its ability for fat synthesis and deposition. However, transcriptional dynamic regulation of adipose differentiation driven by complex signal cascades remains largely unexplored in this model. This study aimed to explore adipogenic transcriptional dynamics before (proliferation) and after (differentiation) initial preadipocyte differentiation in ducks.

RESULTS

Exogenous oleic acid alone successfully induced duck subcutaneous preadipocyte differentiation. We explored 36 mRNA-seq libraries in order to study transcriptome dynamics during proliferation and differentiation processes at 6 time points. Using robust statistical analysis, we identified 845, 652, 359, 2401 and 1933 genes differentially expressed between -48 h and 0 h, 0 h and 12 h, 12 h and 24 h, 24 h and 48 h, 48 h and 72 h, respectively (FDR < 0.05, FC > 1.5). At the proliferation stage, proliferation related pathways and basic cellular and metabolic processes were inhibited, while regulatory factors that initiate differentiation enter the ready-to-activate state, which provides a precondition for initiating adipose differentiation. According to weighted gene co-expression network analysis, pathways positively related to adipogenic differentiation are significantly activated at the differentiation stage, while WNT, FOXO and other pathways that inhibit preadipocyte differentiation are negatively regulated. Moreover, we identified and classified more than 100 transcription factors that showed significant changes during differentiation, and found novel transcription factors that were not reported to be related to preadipoctye differentiation. Finally, we manually assembled a proposed regulation network model of subcutaneous preadipocyte differentiation base on the expression data, and suggested that E2F1 may serve as an important link between the processes of duck subcutaneous preadipocyte proliferation and differentiation.

CONCLUSIONS

For the first time we comprehensively analyzed the transcriptome dynamics of duck subcutaneous preadipocyte proliferation and differentiation. The current study provides a solid basis for understanding the synthesis and deposition of subcutaneous fat in ducks. Furthermore, the information generated will allow future investigations of specific genes involved in particular stages of duck adipogenesis.

Genome-Wide Identification of Long Non-Coding RNAs and Their Regulatory Networks Involved in Apis mellifera ligustica Response to Nosema ceranae Infection

Long non-coding RNAs (lncRNAs) are a diverse class of transcripts that structurally resemble mRNAs but do not encode proteins, and lncRNAs have been proven to play pivotal roles in a wide range of biological processes in animals and plants. However, knowledge of expression patterns and potential roles of honeybee lncRNA response to Nosema ceranae infection is completely unknown. Here, we performed whole transcriptome strand-specific RNA sequencing of normal midguts of Apis mellifera ligustica workers (Am7CK, Am10CK) and N. ceranae-inoculated midguts (Am7T, Am10T), followed by comprehensive analyses using bioinformatic and molecular approaches. A total of 6353 A. m. ligustica lncRNAs were identified, including 4749 conserved lncRNAs and 1604 novel lncRNAs. These lncRNAs had minimal sequence similarities with other known lncRNAs in other species; however, their structural features were similar to counterparts in mammals and plants, including shorter exon and intron length, lower exon number, and lower expression level, compared with protein-coding transcripts. Further, 111 and 146 N. ceranae-responsive lncRNAs were identified from midguts at 7-days post-inoculation (dpi) and 10 dpi compared with control midguts. Twelve differentially expressed lncRNAs (DElncRNAs) were shared by Am7CK vs. Am7T and Am10CK vs. Am10T comparison groups, while the numbers of unique DElncRNAs were 99 and 134, respectively. Functional annotation and pathway analysis showed that the DElncRNAs may regulate the expression of neighboring genes by acting in cis and trans fashion. Moreover, we discovered 27 lncRNAs harboring eight known miRNA precursors and 513 lncRNAs harboring 2257 novel miRNA precursors. Additionally, hundreds of DElncRNAs and their target miRNAs were found to form complex competitive endogenous RNA (ceRNA) networks, suggesting that these DElncRNAs may act as miRNA sponges. Furthermore, DElncRNA-miRNA-mRNA networks were constructed and investigated, the results demonstrated that a portion of the DElncRNAs were likely to participate in regulating the host material and energy metabolism as well as cellular and humoral immune host responses to N. ceranae invasion. Our findings revealed here offer not only a rich genetic resource for further investigation of the functional roles of lncRNAs involved in the A. m. ligustica response to N. ceranae infection, but also a novel insight into understanding the host-pathogen interaction during honeybee microsporidiosis.

Comprehensive Analysis of lncRNAs and circRNAs Reveals the Metabolic Specialization in Oxidative and Glycolytic Skeletal Muscles

The biochemical and functional differences between oxidative and glycolytic muscles could affect human muscle health and animal meat quality. However, present understanding of the epigenetic regulation with respect to lncRNAs and circRNAs is rudimentary. Here, porcine oxidative and glycolytic skeletal muscles, which were at the growth curve inflection point, were sampled to survey variant global expression of lncRNAs and circRNAs using RNA-seq. A total of 4046 lncRNAs were identified, including 911 differentially expressed lncRNAs (p < 0.05). The cis-regulatory analysis identified target genes that were enriched for specific GO terms and pathways (p < 0.05), including the oxidation-reduction process, glycolytic process, and fatty acid metabolic. All these were closely related to different phenotypes between oxidative and glycolytic muscles. Additionally, 810 circRNAs were identified, of which 137 were differentially expressed (p < 0.05). Interestingly, some circRNA-miRNA-mRNA networks were found, which were closely linked to muscle fiber-type switching and mitochondria biogenesis in muscles. Furthermore, 44.69%, 39.19%, and 54.01% of differentially expressed mRNAs, lncRNAs, and circRNAs respectively were significantly enriched in pig quantitative trait loci (QTL) regions for growth and meat quality traits. This study reveals a mass of candidate lncRNAs and circRNAs involved in muscle physiological functions, which may improve understanding of muscle metabolism and development from an epigenetic perspective.

GROWTH AND DEVELOPMENT SYMPOSIUM: STEM AND PROGENITOR CELLS IN ANIMAL GROWTH: Long noncoding RNAs in adipogenesis and adipose development of meat animals12

Sequencing technology, especially next-generation RNA sequencing, has greatly facilitated the identification and annotation of long noncoding RNAs (lncRNAs). In mammals, a large number of lncRNAs have been identified, which regulate various biological processes. An increasing number of lncRNAs have been identified which could function as key regulators of adipogenesis (adipocyte formation), a key step of the development of adipose tissue. Because proper adipose tissue development is a key factor affecting animal growth efficiency, lean/fat ratio, and meat quality, summarizing the roles and recent advances of lncRNAs in adipogenesis is needed in order to develop strategies to effectively manage fat deposition. In this review, we updated lncRNAs contributed to the regulation of adipogenesis, focusing on their roles in fat development of farm animals.

Integrated analysis of lncRNA and mRNA repertoires in Marek's disease infected spleens identifies genes relevant to resistance

BACKGROUND

Marek's disease virus (MDV) is an oncogenic herpesvirus that can cause T-cell lymphomas in chicken. Long noncoding RNA (lncRNA) is strongly associated with various cancers and many other diseases. In chickens, lncRNAs have not been comprehensively identified. Here, we profiled mRNA and lncRNA repertoires in three groups of spleens from MDV-infected and non-infected chickens, including seven tumorous spleens (TS) from MDV-infected chickens, five spleens from the survivors (SS) without lesions after MDV infection, and five spleens from noninfected chickens (NS), to explore the underlying mechanism of host resistance in Marek's disease (MD).

RESULTS

By using a precise lncRNA identification pipeline, we identified 1315 putative lncRNAs and 1166 known lncRNAs in spleen tissue. Genomic features of putative lncRNAs were characterized. Differentially expressed (DE) mRNAs, putative lncRNAs, and known lncRNAs were profiled among three groups. We found that several specific intergroup differentially expressed genes were involved in important biological processes and pathways, including B cell activation and the Wnt signaling pathway; some of these genes were also found to be the hub genes in the co-expression network analyzed by WGCNA. Network analysis depicted both intergenic correlation and correlation between genes and MD traits. Five DE lncRNAs including MSTRG.360.1, MSTRG.6725.1, MSTRG.6754.1, MSTRG.15539.1, and MSTRG.7747.5 strongly correlated with MD-resistant candidate genes, such as IGF-I, CTLA4, HDAC9, SWAP70, CD72, JCHAIN, CXCL12, and CD8B, suggesting that lncRNAs may affect MD resistance and tumorigenesis in chicken spleens through their target genes.

CONCLUSIONS

Our results provide both transcriptomic and epigenetic insights on MD resistance and its pathological mechanism. The comprehensive lncRNA and mRNA transcriptomes in MDV-infected chicken spleens were profiled. Co-expression analysis identified integrated lncRNA-mRNA and gene-gene interaction networks, implying that hub genes or lncRNAs exert critical influence on MD resistance and tumorigenesis.

Identification and analysis of long non-coding RNAs in response to H5N1 influenza viruses in duck (Anas platyrhynchos)

BACKGROUND

Long non-coding RNAs (lncRNAs) are important component of mammalian genomes, where their numbers are even larger than that of protein-coding genes. For example, human (Homo sapiens) (96,308 vs. 20,376) and mouse (Mus musculus) (87,774 vs. 22,630) have more lncRNA genes than protein-coding genes in the NONCODEv5 database. Recently, mammalian lncRNAs were reported to play critical roles in immune response to influenza A virus infections. Such observation inspired us to identify lncRNAs related to immune response to influenza A virus in duck, which is the most important natural host of influenza A viruses.

RESULTS

We explored features of 62,447 lncRNAs from human, mouse, chicken, zebrafish and elegans, and developed a pipeline to identify lncRNAs using the identified features with transcriptomic data. We then collected 151,970 assembled transcripts from RNA-Seq data of 21 individuals from three tissues and annotated 4094 duck lncRNAs. Comparing to duck protein-coding transcripts, we found that 4094 lncRNAs had smaller number of exons (2.4 vs. 10.2) and longer length of transcripts (1903.0 bp vs. 1686.9 bp) on average. Among them, 3586 (87.6%) lncRNAs located in intergenic regions and 619 lncRNAs showed differential expression in ducks infected by H5N1 virus when compared to control individuals. 58 lncRNAs were involved into two co-expressional modules related to anti-influenza A virus immune response. Moreover, we confirmed that eight lncRNAs showed remarkably differential expression both in vivo (duck individuals) and in vitro (duck embryo fibroblast cells, DEF cells) after infected with H5N1 viruses, implying they might play important roles in response to influenza A virus infection.

CONCLUSIONS

This study presented an example to annotate lncRNA in new species based on model species using transcriptome data. These data and analysis provide information for duck lncRNAs' function in immune response to influenza A virus.

Tissue Expression Difference between mRNAs and lncRNAs

Messenger RNA (mRNA) and long noncoding RNA (lncRNA) are two main subgroups of RNAs participating in transcription regulation. With the development of next generation sequencing, increasing lncRNAs are identified. Many hidden functions of lncRNAs are also revealed. However, the differences in lncRNAs and mRNAs are still unclear. For example, we need to determine whether lncRNAs have stronger tissue specificity than mRNAs and which tissues have more lncRNAs expressed. To investigate such tissue expression difference between mRNAs and lncRNAs, we encoded 9339 lncRNAs and 14,294 mRNAs with 71 expression features, including 69 maximum expression features for 69 types of cells, one feature for the maximum expression in all cells, and one expression specificity feature that was measured as Chao-Shen-corrected Shannon's entropy. With advanced feature selection methods, such as maximum relevance minimum redundancy, incremental feature selection methods, and random forest algorithm, 13 features presented the dissimilarity of lncRNAs and mRNAs. The 11 cell subtype features indicated which cell types of the lncRNAs and mRNAs had the largest expression difference. Such cell subtypes may be the potential cell models for lncRNA identification and function investigation. The expression specificity feature suggested that the cell types to express mRNAs and lncRNAs were different. The maximum expression feature suggested that the maximum expression levels of mRNAs and lncRNAs were different. In addition, the rule learning algorithm, repeated incremental pruning to produce error reduction algorithm, was also employed to produce effective classification rules for classifying lncRNAs and mRNAs, which gave competitive results compared with random forest and could give a clearer picture of different expression patterns between lncRNAs and mRNAs. Results not only revealed the heterogeneous expression pattern of lncRNA and mRNA, but also gave rise to the development of a new tool to identify the potential biological functions of such RNA subgroups.

Method using a co-culture system with high-purity intramuscular preadipocytes and satellite cells from chicken pectoralis major muscle

Intramuscular fat is important in improving meat quality; however, the lack of high-purity intramuscular preadipocyte (IMP) in vitro has severely limited the in-depth research on the mutual regulation of myocytes and adipocytes in chicken. In this study, we establish a new method by combining the mature adipocyte ceiling method and the transwell co-culture system. Mature intramuscular adipocyte (MIA) and muscle satellite cell (MSC) were obtained from digested pectoralis major, and MIAs were transformed into IMPs by dedifferentiation with ceiling culture. MSCs were then purified by differential adhesion for 2 h. The results by inverted-microscope observation, MTT assay, Oil Red O staining, and q-PCR revealed that the de-differentiated cells from MIA were identified as the IMPs, and had the same the cellular morphology, the capacity on differentiation, proliferation and passage with the abdominal fat preadipocytes (P > 0.05). The applicability of the obtained IMPs in co-cultured experiment with the MSC revealed that it could meet the requirements of the experimental study. Finally, a co-culture system of IMPs and MSCs was established using a transwell chamber. The co-cultured results indicated that MSCs in the proliferative stage tend to accelerate the differentiation of IMPs to induce more fat content in co-cultured IMPs than in single-culture IMPs (P < 0.05), in the non-proliferative stage, the results tend to show the opposite (P < 0.05). The mRNA levels of related genes significantly changed in accordance with the fat content in cells. The results strongly supported the view that the established co-culture system was effective and feasible. In summary, we successfully found a new method to explore the interaction between myocytes and adipocytes of chicken. Our findings can deepen research on the regulation of chicken myocytes and adipocytes.

Identification and characterization of long non-coding RNAs in porcine granulosa cells exposed to 2,3,7,8-tetrachlorodibenzo-p-dioxin

Background

Long non-coding RNAs (lncRNAs) may regulate gene expression in numerous biological processes including cellular response to xenobiotics. The exposure of living organisms to 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD), a persistent environmental contaminant, results in reproductive defects in many species including pigs. The aims of the study were to identify and characterize lncRNAs in porcine granulosa cells as well as to examine the effects of TCDD on the lncRNA expression profile in the cells.

Results

One thousand six hundred sixty-six lncRNAs were identified and characterized in porcine granulosa cells. The identified lncRNAs were found to be shorter than mRNAs. In addition, the number of exons was lower in lncRNAs than in mRNAs and their exons were longer. TCDD affected the expression of 22 lncRNAs (differentially expressed lncRNAs [DELs]; log₂ fold change  ≥ 1, P-adjusted < 0.05) in the examined cells. Potential functions of DELs were indirectly predicted via searching their target cis- and trans-regulated protein-coding genes. The co-expression analysis revealed that DELs may influence the expression of numerous genes, including those involved in cellular response to xenobiotics, dioxin metabolism, endoplasmic reticulum stress and cell proliferation. Aryl hydrocarbon receptor (AhR) and cytochrome P450 1A1 (CYP1A1) were found among the trans-regulated genes.

Conclusions

These findings indicate that the identified lncRNAs may constitute a part of the regulatory mechanism of TCDD action in granulosa cells. To our knowledge, this is the first study describing lncRNAs in porcine granulosa cells as well as TCDD effects on the lncRNA expression profile. These results may trigger new research directions leading to better understanding of molecular processes induced by xenobiotics in the ovary.

Genome-wide identification of tissue-specific long non-coding RNA in three farm animal species

Background

Numerous long non-coding RNAs (lncRNAs) have been identified and their roles in gene regulation in humans, mice, and other model organisms studied; however, far less research has been focused on lncRNAs in farm animal species. While previous studies in chickens, cattle, and pigs identified lncRNAs in specific developmental stages or differentially expressed under specific conditions in a limited number of tissues, more comprehensive identification of lncRNAs in these species is needed. The goal of the FAANG Consortium (Functional Annotation of Animal Genomes) is to functionally annotate animal genomes, including the annotation of lncRNAs. As one of the FAANG pilot projects, lncRNAs were identified across eight tissues in two adult male biological replicates from chickens, cattle, and pigs.

Results

Comprehensive lncRNA annotations for the chicken, cattle, and pig genomes were generated by utilizing RNA-seq from eight tissue types from two biological replicates per species at the adult developmental stage. A total of 9393 lncRNAs in chickens, 7235 lncRNAs in cattle, and 14,429 lncRNAs in pigs were identified. Including novel isoforms and lncRNAs from novel loci, 5288 novel lncRNAs were identified in chickens, 3732 in cattle, and 4870 in pigs. These transcripts match previously known patterns of lncRNAs, such as generally lower expression levels than mRNAs and higher tissue specificity. An analysis of lncRNA conservation across species identified a set of conserved lncRNAs with potential functions associated with chromatin structure and gene regulation. Tissue-specific lncRNAs were identified. Genes proximal to tissue-specific lncRNAs were enriched for GO terms associated with the tissue of origin, such as leukocyte activation in spleen.

Conclusions

LncRNAs were identified in three important farm animal species using eight tissues from adult individuals. About half of the identified lncRNAs were not previously reported in the NCBI annotations for these species. While lncRNAs are less conserved than protein-coding genes, a set of positionally conserved lncRNAs were identified among chickens, cattle, and pigs with potential functions related to chromatin structure and gene regulation. Tissue-specific lncRNAs have potential regulatory functions on genes enriched for tissue-specific GO terms. Future work will include epigenetic data from ChIP-seq experiments to further refine these annotations.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5037-7) contains supplementary material, which is available to authorized users.

Genomic Insights Into the Multiple Factors Controlling Abdominal Fat Deposition in a Chicken Model

Genetic selection for an increased growth rate in meat-type chickens has been accompanied by excessive fat accumulation particularly in abdominal cavity. These progressed to indirect and often unhealthy effects on meat quality properties and increased feed cost. Advances in genomics technology over recent years have led to the surprising discoveries that the genome is more complex than previously thought. Studies have identified multiple-genetic factors associated with abdominal fat deposition. Meanwhile, the obesity epidemic has focused attention on adipose tissue and the development of adipocytes. The aim of this review is to summarize the current understanding of genetic/epigenetic factors associated with abdominal fat deposition, or as it relates to the proliferation and differentiation of preadipocytes in chicken. The results discussed here have been identified by different genomic approaches, such as QTL-based studies, the candidate gene approach, epistatic interaction, copy number variation, single-nucleotide polymorphism screening, selection signature analysis, genome-wide association studies, RNA sequencing, and bisulfite sequencing. The studies mentioned in this review have described multiple-genetic factors involved in an abdominal fat deposition. Therefore, it is inevitable to further study the multiple-genetic factors in-depth to develop novel molecular markers or potential targets, which will provide promising applications for reducing abdominal fat deposition in meat-type chicken.

Sequencing and characterization of lncRNAs in the breast muscle of Gushi and Arbor Acres chickens

Chicken muscle quality is one of the most important factors determining the economic value of poultry, and muscle development and growth are affected by genetics, environment, and nutrition. However, little is known about the molecular regulatory mechanisms of long non-coding RNAs (lncRNAs) in chicken skeletal muscle development. Our study aimed to better understand muscle development in chickens and thereby improve meat quality. In this study, Ribo-Zero RNA-Seq was used to investigate differences in the expression profiles of muscle development related genes and associated pathways between Gushi (GS) and Arbor Acres (AA) chickens. We identified two muscle tissue specific expression lncRNAs. In addition, the target genes of these lncRNAs were significantly enriched in certain biological processes and molecular functions, as demonstrated by Gene Ontology (GO) analysis, and these target genes participate in five signaling pathway, as revealed by an analysis of the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. Taken together, these data suggest that different lncRNAs might be involved in regulating chicken muscle development and growth and provide new insight into the molecular mechanisms of lncRNAs.

Differential regulation of mRNAs and lncRNAs related to lipid metabolism in two pig breeds

Long non-coding RNAs (lncRNAs) can regulate lipid metabolism and adipogenesis. However, there is little research on the role of lncRNAs in fat deposition in pig. In this study, RNA-seq technology was used to analyze the gene expression profiles of subcutaneous adipose tissue in Laiwu (LW) and Large White (LY) pigs. Then, key lncRNAs and genes associated with lipid metabolism and adipogenic differentiation were identified. Fifty four lncRNAs and 482 known mRNAs were differentially expressed in the two pig breeds. GO (Gene Ontology) and KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway analyses revealed that differentially expressed genes and the target genes of differentially expressed lncRNAs were significantly enriched in PPAR signaling pathway and biological processes including fat cell differentiation and fatty acid metabolism. Key lncRNAs might regulate adipogenic differentiation and fatty acid metabolism by regulating genes involved in above signaling pathway and biological processes. Specifically, XLOC_014379, XLOC_011279, XLOC_064871, XLOC_019518 and XLOC_013639 might target SCD, LPIN1, TRIB3, EGR2 and FABP3, respectively, and then play critical regulatory role. These results are useful for understanding fat deposition in pig, breeding livestock with high quality meat, and preventing and treating lipid metabolic disease.