Characterization of differentially expressed and lipid metabolism-related lncRNA-mRNA interaction networks during the growth of liver tissue through rabbit models

Background Characterization the long non-coding RNAs (lncRNAs) and their regulated mRNAs involved in lipid metabolism during liver growth and development is of great value for discovering new genomic biomarkers and therapeutic targets for fatty liver and metabolic syndrome. Materials and methods Liver samples from sixteen rabbit models during the four growth stages (birth, weaning, sexual maturity, and somatic maturity) were used for RNA-seq and subsequent bioinformatics analyses. Differentially expressed (DE) lncRNAs and mRNAs were screened, and the cis/trans-regulation target mRNAs of DE lncRNAs were predicted. Then the function enrichment analyses of target mRNAs were performed through Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, respectively. The target protein interaction (PPI) and lncRNA-mRNA co-expression networks were constructed using string version 11.0 platform and R Stats. Finally, six lncRNAs and six mRNAs were verified taking RT-qPCR. Results Liver Oil Red O detection found that the liver showed time-dependent accumulation of lipid droplets. 41,095 lncRNAs, 30,744 mRNAs, and amount to 3,384 DE lncRNAs and 2980 DE mRNAs were identified from 16 cDNA sequencing libraries during the growth of liver. 689 out of all DE lncRNAs corresponded to 440 DE mRNAs by cis-regulation and all DE mRNAs could be regulated by DE lncRNAs by trans-regulation. GO enrichment analysis showed significant enrichment of 892 GO terms, such as protein binding, cytosol, extracellular exsome, nucleoplasm, and oxidation-reduction process. Besides, 52 KEGG pathways were significantly enriched, including 11 pathways of lipid metabolism were found, like Arachidonic acid metabolism, PPAR signaling pathway and Biosynthesis of unsaturated fatty acids. After the low expression DE mRNAs and lncRNAs were excluded, we further obtained the 54 mRNAs were regulated by 249 lncRNAs. 351 interaction pairs were produced among 38 mRNAs and 215 lncRNAs through the co-expression analysis. The PPI network analysis found that 10 mRNAs such as 3β-Hydroxysteroid-Δ24 Reductase (DHCR24), lathosterol 5-desaturase (SC5D), and acetyl-CoA synthetase 2 (ACSS2) were highly interconnected hub protein-coding genes. Except for MSTRG.43041.1, the expression levels of the 11 genes by RT-qPCR were the similar trends to the RNA-seq results. Conclusion The study revealed lncRNA-mRNA interation networks that regulate lipid metabolism during liver growth, providing potential research targets for the prophylaxis and treatment of related diseases caused by liver lipid metabolism disorders.

[1]  Liqin Yang,et al.  Decoding the role of long noncoding RNAs in the healthy aging of centenarians , 2021, Briefings Bioinform..

[2]  Stephen T. C. Wong,et al.  Acetyl-CoA Synthetase 2: A Critical Linkage in Obesity-Induced Tumorigenesis in Myeloma. , 2021, Cell metabolism.

[3]  Xianbo Jia,et al.  Screening and Identification of Differentially Expressed and Adipose Growth-Related Protein-Coding Genes During the Deposition of Perirenal Adipose Tissue in Rabbits , 2020, Diabetes, metabolic syndrome and obesity : targets and therapy.

[4]  R. Báez.  Advances in Lipid Metabolism , 2020 .

[5]  Runping Liu,et al.  Long non-coding RNA H19 in the liver-gut axis: A diagnostic marker and therapeutic target for liver diseases. , 2020, Experimental and molecular pathology.

[6]  Weifeng He,et al.  The role of long noncoding RNAs in hepatocellular carcinoma , 2020, Molecular Cancer.

[7]  J. Rabinowitz,et al.  Dietary fructose feeds hepatic lipogenesis via microbiota-derived acetate , 2020, Nature.

[8]  Xiaolan Fan,et al.  The Roles and Mechanisms of lncRNAs in Liver Fibrosis , 2020, International journal of molecular sciences.

[9]  Min Li,et al.  Predicting Human lncRNA-Disease Associations Based on Geometric Matrix Completion , 2019, IEEE Journal of Biomedical and Health Informatics.

[10]  C. Yuan,et al.  Pathophysiological Functions of the lncRNA TUG1. , 2019, Current pharmaceutical design.

[11]  M. Shiomi.  The History of the WHHL Rabbit, an Animal Model of Familial Hypercholesterolemia (I) , 2019, Journal of atherosclerosis and thrombosis.

[12]  Xianbo Jia,et al.  Screening and identification of MicroRNAs expressed in perirenal adipose tissue during rabbit growth , 2019, bioRxiv.

[13]  S. Chong,et al.  Roles and Regulation of Long non-coding RNAs in Hepatocellular Carcinoma. , 2019, Cancer research.

[14]  Chuanzheng Sun,et al.  LncRNA MEG3 functions as a ceRNA in regulating hepatic lipogenesis by competitively binding to miR-21 with LRP6. , 2019, Metabolism: clinical and experimental.

[15]  Shan Jiang,et al.  Cholesterol depletion induced by RNA interference targeting DHCR24 protects cells from liposome-induced cytotoxicity , 2019, Preparative biochemistry & biotechnology.

[16]  C. Vantaggiato,et al.  Antipsychotics Promote Metabolic Disorders Disrupting Cellular Lipid Metabolism and Trafficking , 2019, Trends in Endocrinology & Metabolism.

[17]  C. Mantzoros,et al.  Obesity and nonalcoholic fatty liver disease: From pathophysiology to therapeutics. , 2019, Metabolism: clinical and experimental.

[18]  H. Cao,et al.  Therapeutic Effect of Chitooligosaccharide Tablets on Lipids in High-Fat Diets Induced Hyperlipidemic Rats , 2019, Molecules.

[19]  Jessica Yun,et al.  Role of non-coding RNAs in liver disease progression to hepatocellular carcinoma , 2019, Archives of Pharmacal Research.

[20]  M. Brunelli,et al.  Oil Red O Is a Useful Tool to Assess Donor Liver Steatosis on Frozen Sections During Transplantation. , 2018, Transplantation proceedings.

[21]  B. Tu,et al.  ACSS2 promotes systemic fat storage and utilization through selective regulation of genes involved in lipid metabolism , 2018, Proceedings of the National Academy of Sciences.

[22]  G. Rankin,et al.  Effects of cytochrome P450 single nucleotide polymorphisms on methadone metabolism and pharmacodynamics , 2018, Biochemical pharmacology.

[23]  Mohd Effendy ABD WAHID,et al.  Evaluation of Herb-Drug Interaction of Synacinn™ and Individual Biomarker through Cytochrome 450 Inhibition Assay. , 2018, Drug metabolism letters.

[24]  R. Guigó,et al.  The reference epigenome and regulatory chromatin landscape of chronic lymphocytic leukemia , 2018, Nature Medicine.

[25]  Byung-Hoon Lee,et al.  Duodenal Niemann-Pick C1-like 1 expression was negatively correlated with liver X receptor expression in nonalcoholic fatty liver disease , 2018, The Korean journal of internal medicine.

[26]  Xiao Zhu,et al.  Long Noncoding RNAs: Advances in Lipid Metabolism. , 2018, Advances in clinical chemistry.

[27]  Meiyi Song,et al.  Circular RNA in Liver: Health and Diseases. , 2018, Advances in experimental medicine and biology.

[28]  R. Goyal,et al.  LncRNA as a Therapeutic Target for Angiogenesis , 2017, Current topics in medicinal chemistry.

[29]  Timofei S Zatsepin,et al.  lncRNA in the liver: Prospects for fundamental research and therapy by RNA interference. , 2016, Biochimie.

[30]  E. Lander,et al.  Local regulation of gene expression by lncRNA promoters, transcription and splicing , 2016, Nature.

[31]  E. Tsochatzis,et al.  The multiple-hit pathogenesis of non-alcoholic fatty liver disease (NAFLD). , 2016, Metabolism: clinical and experimental.

[32]  Rohit Loomba,et al.  Recommendations for Diagnosis, Referral for Liver Biopsy, and Treatment of Nonalcoholic Fatty Liver Disease and Nonalcoholic Steatohepatitis. , 2015, Mayo Clinic proceedings.

[33]  Steven L Salzberg,et al.  HISAT: a fast spliced aligner with low memory requirements , 2015, Nature Methods.

[34]  S. Salzberg,et al.  StringTie enables improved reconstruction of a transcriptome from RNA-seq reads , 2015, Nature Biotechnology.

[35]  Davide Heller,et al.  STRING v10: protein–protein interaction networks, integrated over the tree of life , 2014, Nucleic Acids Res..

[36]  W. Huber,et al.  Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2 , 2014, Genome Biology.

[37]  Xiaobing Fu,et al.  lncRNAs: insights into their function and mechanics in underlying disorders. , 2014, Mutation research. Reviews in mutation research.

[38]  J. Murray,et al.  Polymorphisms of CYP51A1 from Cholesterol Synthesis: Associations with Birth Weight and Maternal Lipid Levels and Impact on CYP51 Protein Structure , 2013, PloS one.

[39]  Yi Zhao,et al.  Utilizing sequence intrinsic composition to classify protein-coding and long non-coding transcripts , 2013, Nucleic acids research.

[40]  Maite Huarte,et al.  Long non-coding RNAs: challenges for diagnosis and therapies. , 2013, Nucleic acid therapeutics.

[41]  Cole Trapnell,et al.  Integrative annotation of human large intergenic noncoding RNAs reveals global properties and specific subclasses. , 2011, Genes & development.

[42]  Vincent Wai-Sun Wong,et al.  Disease progression of non-alcoholic fatty liver disease: a prospective study with paired liver biopsies at 3 years , 2010, Gut.

[43]  Cole Trapnell,et al.  Role of Rodent Secondary Motor Cortex in Value-based Action Selection Nih Public Access Author Manuscript , 2006 .

[44]  Yoshihiro Yamanishi,et al.  KEGG for linking genomes to life and the environment , 2007, Nucleic Acids Res..

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

[46]  Juan Miguel García-Gómez,et al.  BIOINFORMATICS APPLICATIONS NOTE Sequence analysis Manipulation of FASTQ data with Galaxy , 2005 .

[47]  P. Shannon,et al.  Cytoscape: a software environment for integrated models of biomolecular interaction networks. , 2003, Genome research.

[48]  Thomas D. Schmittgen,et al.  Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. , 2001, Methods.