Identification of long non-coding RNAs involved in fat 1 deposition in pigs using two high-throughput sequencing 2 methods 3

: 23 Background: Adipose is an important body tissue in pigs, and fatty traits are critical in pig 24 production. The function of long non-coding RNA (lncRNA) in fat deposition and metabolism has 25 been proven in previous studies. In this study, we focused on lncRNAs associated with fattening 26 traits in pigs. The adipose tissue of six Landrace pigs with either extremely-high or -low backfat 27 thickness ( n high = 3, n low = 3) was collected, after which we performed strand-specific RNA 28 sequencing using biological replicates and pooling methods. 29 Results: A total of 19,631 genes and 2,013 lncRNAs were identified using the coding potential 30 calculator, coding-non-coding index, and Pfam database, including 334 known transcripts and 31 1,679 novel transcripts. Using edgeR, we determined that 220 lncRNAs and 1,512 genes were 32 differentially expressed (|Fold Change| > 2 and false discovery rate < 0.05) between the two 33 groups in biological replicate RNA sequencing (RNA-seq), and 127 lncRNAs and 2,240 genes 34 were differently expressed in pooling RNA-seq. Further Kyoto Encyclopedia of Genes and 35 Genomes and Gene Ontology enrichment analysis of the differentially expressed genes found that 36 some of the genes were involved in several key pathways related to fat development. After 37 targeting gene prediction, we determined that some cis-target genes of the differentially expressed 38 lncRNAs play an important role in fat deposition. For example, ACSL3 is cis-targeted by lncRNA 39 TCONS-00052400, and it can activate the conversion of long-chain fatty acids. In addition, 40 lncRNA TCONS_00041740 was up-regulated in the high backfat thickness group, and its 41 cis-target gene ACACB was also up-regulated in this group. It has been reported that ACACB is the 42 rate-limiting enzyme in fatty acid oxidation. Conclusions: Since these genes have necessary functions in fat metabolism, the results imply that the lncRNAs detected in our study may affect fat deposition in pigs through regulation of their 45 target genes. In summary, our study explored the regulation of lncRNA and their target genes on 46 fat deposition in pigs and provided new insights for further investigation of the biological 47 functions of lncRNA. In addition, the gene PLA2G12B in the BH group. PLA2G12B A2 group is encoded by this gene and belongs to the phospholipase A2 ( PLA2 ) group of enzymes, which plays a role in lipid hydrolysis by releasing free fatty acids and lysophospholipids [52] . Studies have shown that a reduction of PLA2G12B decreases the amount of serum triglyceride (TG)-rich VLDL particles secreted by the liver, resulting in a reduction in TG content [53] . PLA2G12B can also participate in the pathogenesis of idiopathic membranous nephropathy (iMN) by regulating lipid metabolism [54] . In a previous study by Guan, PLA2G12B -null mice had obvious accumulation of large lipid droplets in the liver, displaying the fatty liver phenotype [55] . These results indicate that the genes exclusively expressed in the high or low backfat groups may also have a certain regulatory effect on lipid metabolism. to and TCONS-00052400 to ACSL3 These results can provide useful information for understanding the regulation of by lncRNA in pigs. However, further genetic experiments are still needed to validate the association of the lncRNA and mRNA functions presented in this study.

[1]  Hemin Ni,et al.  Transcriptome analysis of miRNA and mRNA in the livers of pigs with highly diverged backfat thickness , 2019, Scientific Reports.

[2]  Xihui Sheng,et al.  Comparative adipose transcriptome analysis digs out genes related to fat deposition in two pig breeds , 2019, Scientific Reports.

[3]  C. Zou,et al.  Transcriptome Analysis Reveals the Effect of Long Intergenic Noncoding RNAs on Pig Muscle Growth and Fat Deposition , 2019, BioMed research international.

[4]  Yanfang Jiang,et al.  Lipid metabolism participates in human membranous nephropathy identified by whole-genome gene expression profiling. , 2019, Clinical science.

[5]  G. Maulucci,et al.  Impact of apolipoprotein A1- or lecithin:cholesterol acyltransferase-deficiency on white adipose tissue metabolic activity and glucose homeostasis in mice. , 2019, Biochimica et biophysica acta. Molecular basis of disease.

[6]  A. Mizokami,et al.  Phospholipase C-related catalytically inactive protein: A novel signaling molecule for modulating fat metabolism and energy expenditure. , 2019, Journal of oral biosciences.

[7]  Panpeng Wei,et al.  Identification and comparison of long non-conding RNA in Jinhua and Landrace pigs. , 2018, Biochemical and biophysical research communications.

[8]  X. Ruan,et al.  Knocking down Stard3 decreases adipogenesis with decreased mitochondrial ROS in 3T3-L1 cells. , 2018, Biochemical and biophysical research communications.

[9]  H. Randeva,et al.  Differential expression of Lp-PLA2 in obesity and type 2 diabetes and the influence of lipids , 2018, Diabetologia.

[10]  Hongbo Zhang,et al.  Genome-Wide Analysis of Long Non-Coding RNAs in Potato and Their Potential Role in Tuber Sprouting Process , 2017, International journal of molecular sciences.

[11]  B. Yang,et al.  Transcriptome analyses reveal genes and pathways associated with fatty acid composition traits in pigs , 2017, Animal genetics.

[12]  Yong-tang Zheng,et al.  Annotation and cluster analysis of spatiotemporal- and sex-related lncRNA expression in rhesus macaque brain , 2017, Genome research.

[13]  Nathalie Spruyt,et al.  The long non-coding RNA 91H increases aggressive phenotype of breast cancer cells and up-regulates H19/IGF2 expression through epigenetic modifications. , 2017, Cancer letters.

[14]  H. Pan,et al.  Activation of farnesoid X receptor promotes triglycerides lowering by suppressing phospholipase A2 G12B expression , 2016, Molecular and Cellular Endocrinology.

[15]  Jiazhong Guo,et al.  Genome-wide identification and characterization of long non-coding RNAs in developmental skeletal muscle of fetal goat , 2016, BMC Genomics.

[16]  Jian Li,et al.  Effect of High Fat Dietary Intake during Maternal Gestation on Offspring Ovarian Health in a Pig Model , 2016, Nutrients.

[17]  Hong Chen,et al.  Long non-coding RNA ADNCR suppresses adipogenic differentiation by targeting miR-204. , 2016, Biochimica et biophysica acta.

[18]  F. Zhu,et al.  Identification of genes for controlling swine adipose deposition by integrating transcriptome, whole-genome resequencing, and quantitative trait loci data , 2016, Scientific Reports.

[19]  Robert D. Finn,et al.  The Pfam protein families database: towards a more sustainable future , 2015, Nucleic Acids Res..

[20]  R. Shiekhattar,et al.  Regulation of transcription by long noncoding RNAs. , 2014, Annual review of genetics.

[21]  Paul Theodor Pyl,et al.  HTSeq—a Python framework to work with high-throughput sequencing data , 2014, bioRxiv.

[22]  A. Atala Re: lncRNA-dependent mechanisms of androgen-receptor-regulated gene activation programs. , 2014, The Journal of urology.

[23]  Björn Usadel,et al.  Trimmomatic: a flexible trimmer for Illumina sequence data , 2014, Bioinform..

[24]  Hui Liu,et al.  Identification and Characterization of Long Non-Coding RNAs Related to Mouse Embryonic Brain Development from Available Transcriptomic Data , 2013, PloS one.

[25]  Cole Trapnell,et al.  TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions , 2013, Genome Biology.

[26]  Huanming Yang,et al.  The sequence and analysis of a Chinese pig genome , 2012, GigaScience.

[27]  Steven L Salzberg,et al.  Fast gapped-read alignment with Bowtie 2 , 2012, Nature Methods.

[28]  David R. Kelley,et al.  Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks , 2012, Nature Protocols.

[29]  Chuan-Yun Li,et al.  KOBAS 2.0: a web server for annotation and identification of enriched pathways and diseases , 2011, Nucleic Acids Res..

[30]  Robert C. Thompson,et al.  NGSQC: cross-platform quality analysis pipeline for deep sequencing data , 2010, BMC Genomics.

[31]  Davis J. McCarthy,et al.  edgeR: a Bioconductor package for differential expression analysis of digital gene expression data , 2009, Bioinform..

[32]  D. Mashek,et al.  Suppression of Long Chain Acyl-CoA Synthetase 3 Decreases Hepatic de Novo Fatty Acid Synthesis through Decreased Transcriptional Activity* , 2009, The Journal of Biological Chemistry.

[33]  F. Visseren,et al.  Adipose tissue dysfunction in obesity, diabetes, and vascular diseases. , 2008, European heart journal.

[34]  Jin-cui Ye,et al.  Long Chain Acyl-CoA Synthetase 3-mediated Phosphatidylcholine Synthesis Is Required for Assembly of Very Low Density Lipoproteins in Human Hepatoma Huh7 Cells* , 2008, Journal of Biological Chemistry.

[35]  R. Coleman,et al.  Long-chain acyl-CoA synthetases and fatty acid channeling. , 2007, Future lipidology.

[36]  Yong Zhang,et al.  CPC: assess the protein-coding potential of transcripts using sequence features and support vector machine , 2007, Nucleic Acids Res..

[37]  P. Stadler,et al.  RNA Maps Reveal New RNA Classes and a Possible Function for Pervasive Transcription , 2007, Science.

[38]  G. Plastow,et al.  Commercial pigs: an untapped resource for human obesity research? , 2006, Drug discovery today.

[39]  Jianglin Fan,et al.  A minipig model of high‐fat/high‐sucrose diet‐induced diabetes and atherosclerosis , 2004, International journal of experimental pathology.

[40]  S. Wakil,et al.  Acetyl-CoA carboxylase 2 mutant mice are protected against obesity and diabetes induced by high-fat/high-carbohydrate diets , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[41]  Brad T. Sherman,et al.  DAVID: Database for Annotation, Visualization, and Integrated Discovery , 2003, Genome Biology.

[42]  N. Lundeheim,et al.  Effect of birth litter size, birth parity number, growth rate, backfat thickness and age at first mating of gilts on their reproductive performance as sows. , 2001, Animal reproduction science.

[43]  Martin M. Matzuk,et al.  Continuous Fatty Acid Oxidation and Reduced Fat Storage in Mice Lacking Acetyl-CoA Carboxylase 2 , 2001, Science.

[44]  M. Ashburner,et al.  Gene Ontology: tool for the unification of biology , 2000, Nature Genetics.

[45]  T. Fujita,et al.  Pancreatic insulo‐acinar portal systems in humans, rats, and some other mammals: Scanning electron microscopy of vascular casts , 1997, Microscopy research and technique.

[46]  K. Z. Masoodi,et al.  Database for annotation, visualisation and integrated discovery , 2022, Bioinformatics for Everyone.

[47]  黄志伟 Hepatocyte nuclear factor-4 alpha regulates liver triglyceride metabolism in part through secreted phospholipase A2 GXIIB , 2011 .

[48]  Cedric E. Ginestet ggplot2: Elegant Graphics for Data Analysis , 2011 .

[49]  Christian Gieger,et al.  Six new loci associated with body mass index highlight a neuronal influence on body weight regulation , 2009, Nature Genetics.

[50]  D. Bellinger,et al.  Swine models of type 2 diabetes mellitus: insulin resistance, glucose tolerance, and cardiovascular complications. , 2006, ILAR journal.

[51]  C. Lottaz,et al.  BIOINFORMATICS APPLICATIONS NOTE , 2001 .

[52]  M. Ashburner,et al.  The Gene Ontology Consortium , 2000 .

[53]  Supplementary Supplementary Supplementary Supplementary Methods Methods Methods Methods Comparison Comparison Comparison Comparison of of of of Cnci Cnci Cnci Cnci Performance Performance Performance Performance with with with with Cpc Cpc Cpc Cpc and and and and Phylocsf Phylocsf Phylocsf Phylocsf , 2022 .