影响藏猪与大约克夏猪肌肉品质的lncRNA的筛选及功能分析

doi:10.11843/j.issn.0366-6964.2022.01.006

摘要/Abstract

摘要： 旨在分析藏猪和大约克夏猪（180日龄）背最长肌组织中与猪肉质性状相关的长链非编码RNA （lncRNAs）和基因，探讨lncRNAs在猪肉品质中的分子调控作用，为猪肉质性状的改良提供理论依据。本研究根据先前公布的数据，以饲养在相同条件下的180日龄健康藏公猪（n=3）和大约克夏公猪（n=3）为研究对象，收集背最长肌肉品质表型数据，利用相同个体背最长肌转录组数据筛选差异表达lncRNAs和差异表达基因（DEGs）。根据差异表达lncRNAs和DEGs位置关系以及表达的相关性预测lncRNAs顺式作用（cis）和反式作用（trans）的靶DEGs，并对trans靶DEGs进行GO生物学过程和KEGG通路富集分析。进一步筛选潜在影响肌肉生长发育和脂肪沉积的候选差异表达lncRNAs和靶DEGs，构建lncRNAs-DEGs-生物学过程调控网络。最后将筛选的候选靶基因与藏猪和大约克夏猪的肉质性状（滴水损失、pH、肉色、大理石花纹、肌内脂肪含量、肌纤维横截面积、眼肌面积和维生素B₁含量）进行关联分析。结果显示，在藏猪和大约克夏猪背最长肌组织中共鉴定出1 546个lncRNAs和20 219个mRNAs，其中包括49个差异表达lncRNAs和154个DEGs，差异表达lncRNAs分别顺式和反式调控7个和138个靶DEGs。功能富集分析发现，大量靶DEGs显著富集到与肌肉生长发育和脂肪沉积相关的生物学过程，包括磷脂酰肌醇3-激酶信号的正调控、肌肉结构发育、脂肪细胞分化、JAK-STAT级联的正调控、MAPK级联的正调控等。调控网络显示，MSTRG.34836.10、MSTRG.151.1调控的靶DEGs最多，其次是MSTRG.21406.1、MSTRG.26709.1。关联分析发现，这些lncRNAs调控的候选靶DEGs有11个与肉质性状显著相关。综上表明，本研究鉴定的部分差异表达lncRNAs与肌肉生长发育和脂肪沉积有关，分析发现这些差异表达lncRNAs可能通过调控HES1、LEP、TGFB2和PER2等基因的表达影响猪肉品质，为进一步解析猪肉质性状的潜在分子机制奠定了基础。

关键词: 藏猪, 大约克夏猪, RNA-seq, lncRNAs, 肉质性状

Abstract: This study aimed to identify the long non-coding RNAs (lncRNAs) and genes related to pork quality traits of the longissimus dorsi in Tibetan pigs and Yorkshire pigs (180 days old), and to further explore their molecular regulatory role, which could provide theoretical basis for the improvement breeding of pork quality traits. According to previously published data, 180-day-old healthy Tibetan boars (n=3) and Yorkshire boars (n=3) raised under the same conditions were used as the research objects to collect phenotype data of the quality of the longissimus dorsi and the transcriptome data of the longissimus dorsi of the same individual were used to screen differentially expressed lncRNAs and differentially expressed genes (DEGs). Cis- and trans-target DEGs were predicted based on the positional relationship and expression correlation of differentially expressed lncRNAs and DEGs, then trans-target DEGs were enriched by GO biological processes and KEGG pathway analysis. Subsequently, lncRNAs-DEGs-biological process regulatory network related to muscle growth and fat deposition was constructed. Then, association events between candidate target genes and meat quality traits were analyzed, which comprised of drip loss, pH, meat color, marbling, intramuscular fat content, muscle fiber cross-sectional area, loin muscle area and vitamin B₁ content. A total of 1 546 lncRNAs and 20 219 mRNAs were identified in the longissimus dorsi of Tibetan pigs and Yorkshire pigs,there were 49 differentially expressed lncRNAs and 154 DEGs of them based on statistical analysis. Of these, 7 cis and 138 trans regulatory DEGs targeted by differentially expressed lncRNAs were predicted. Functional enrichment analysis showed that a large number of candidate DEGs were significantly enriched in biological processes of muscle growth and fat deposition, including positive regulation of phosphatidylinositol 3-kinase signaling, muscle structure development, fat cell differentiation, positive regulation of JAK-STAT cascade, and positive regulation of MAPK cascade, etc. Evidence from regulatory network showed that MSTRG.34836.10 and MSTRG.151.1 had the most target DEGs, followed by MSTRG.21406.1, MSTRG.26709.1. Association analysis indicated that 11 candidate target DEGs were significantly related to the meat quality traits. Taken together, some differentially expressed lncRNAs identified in this study were related to muscle growth and fat deposition. It was also suggested that these differentially expressed lncRNAs had an effect on pork quality by regulating the expression of HES1, LEP, TGFB2, and PER2. This study laid the foundation for further analysis of the potential molecular mechanism of pork quality traits.

Key words: Tibetan pig, Yorkshire pig, RNA-seq, lncRNAs, meat quality traits

中图分类号:

S828.2

宋明坤, 薛明明, 张力戈, 颜铎, 夏宁, 商鹏, 李新建, 韩雪蕾, 乔瑞敏, 李秀领, 李明, 王克君. 影响藏猪与大约克夏猪肌肉品质的lncRNA的筛选及功能分析[J]. 畜牧兽医学报, 2022, 53(1): 53-65.

SONG Mingkun, XUE Mingming, ZHANG Lige, YAN Duo, XIA Ning, SHANG Peng, LI Xinjian, HAN Xuelei, QIAO Ruimin, LI Xiuling, LI Ming, WANG Kejun. Screening and Functional Analysis of lncRNA Affecting Muscle Quality in Tibetan and Yorkshire Pigs[J]. Acta Veterinaria et Zootechnica Sinica, 2022, 53(1): 53-65.

参考文献 45

[1]	ZAPPATERRA M,GIOIOSA S,CHILLEMI G,et al.Dissecting the gene expression networks associated with variations in the major components of the fatty acid Semimembranosus muscle profile in Large White heavy pigs[J]. Animals (Basel),2021,11(3):628.
[2]	ZHANG J,CUI L,MA J,et al.Transcriptome analyses reveal genes and pathways associated with fatty acid composition traits in pigs[J].Anim Genet,2017,48(6):645-652.
[3]	LI M Z,TIAN S L,JIN L,et al.Genomic analyses identify distinct patterns of selection in domesticated pigs and Tibetan wild boars[J].Nat Genet,2013,45(12):1431-1438.
[4]	SHANG P,LI W T,LIU G,et al.Identification of lncRNAs and genes responsible for fatness and fatty acid composition traits between the Tibetan and Yorkshire pigs[J].Int J Genomics,2019,2019:5070975.
[5]	SHANG P,LI W T,TAN Z K,et al.Population genetic analysis of ten geographically isolated Tibetan pig populations[J].Animals (Basel),2020,10(8):1297.
[6]	MA Y F,HAN X M,HUANG C P,et al.Population genomics analysis revealed origin and high-altitude adaptation of Tibetan pigs[J].Sci Rep,2019,9(1):11463.
[7]	李梦柔,韦明邦,张健,等.基于iTRAQ技术挖掘猪肾组织高原低氧适应性关键蛋白[J].畜牧兽医学报, 2021, 52(5):1238-1246.LI M R,WEI M B,ZHANG J,et al.Mining key proteins of high altitude hypoxic adaptability in porcine kidney tissue based on iTRAQ technology[J].Acta Veterinaria et Zootechnica Sinica,2021,52(5):1238-1246.(in Chinese)
[8]	MI S,LI X,ZHANG C H,et al.Characterization and discrimination of Tibetan and Duroc×(Landrace×Yorkshire) pork using label-free quantitative proteomics analysis[J].Food Res Int,2019,119:426-435.
[9]	GAN M L,SHEN L Y,FAN Y,et al.High altitude adaptability and meat quality in Tibetan pigs:a reference for local pork processing and genetic improvement[J].Animals (Basel),2019,9(12):1080.
[10]	STATELLO L,GUO C J,CHEN L L,et al.Gene regulation by long non-coding RNAs and its biological functions[J].Nat Rev Mol Cell Biol,2021,22(2):96-118.
[11]	YANG S M,LIM K H,KIM S H,et al.Molecular landscape of long noncoding RNAs in brain disorders[J].Mol Psychiatry,2021,26(4):1060-1074.
[12]	WANG L J,YANG X,ZHU Y H,et al.Genome-wide identification and characterization of long noncoding RNAs of brown to white adipose tissue transformation in goats[J].Cells,2019,8(8):904.
[13]	CHEN R,LEI S,JIANG T,et al.Roles of lncRNAs and circRNAs in regulating skeletal muscle development[J]. Acta Physiol (Oxf),2020,228(2):e13356.
[14]	WANG J,CHEN M Y,CHEN J F,et al.LncRNA IMFlnc1 promotes porcine intramuscular adipocyte adipogenesis by sponging miR-199a-5p to up-regulate CAV-1[J].BMC Mol Cell Biol,2020,21(1):77.
[15]	LI Q Q,HUANG Z Y,ZHAO W J,et al.Transcriptome analysis reveals long intergenic non-coding RNAs contributed to intramuscular fat content differences between Yorkshire and Wei Pigs[J].Int J Mol Sci,2020, 21(5):1732.
[16]	TAO X,LIANG Y,YANG X M,et al.Transcriptomic profiling in muscle and adipose tissue identifies genes related to growth and lipid deposition[J].PLoS One,2017,12(9):e0184120.
[17]	SUN L,LUO H T,BU D C,et al.Utilizing sequence intrinsic composition to classify protein-coding and long non-coding transcripts[J].Nucleic Acids Res,2013,41(17):e166.
[18]	KONG L,ZHANG Y,YE Z Q,et al.CPC:assess the protein-coding potential of transcripts using sequence features and support vector machine[J].Nucleic Acids Res,2007,35(Web Server issue):W345-W349.
[19]	ANDERS S,PYL P T,HUBER W.HTSeq—a Python framework to work with high-throughput sequencing data[J].Bioinformatics,2015,31(2):166-169.
[20]	LOVE M I,HUBER W,ANDERS S.Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2[J].Genome Biol,2014,15(12):550.
[21]	GUIL S,ESTELLER M.Cis-acting noncoding RNAs:friends and foes[J].Nat Struct Mol Biol,2012,19(11): 1068-1075.
[22]	QUINLAN A R.BEDTools:the Swiss-army tool for genome feature analysis[J].Curr Protoc Bioinformatics, 2014, 47: 11.12.1-11.12.34.
[23]	ZHOU Y Y,ZHOU B,PACHE L,et al.Metascape provides a biologist-oriented resource for the analysis of systems-level datasets[J].Nat Commun,2019,10(1):1523.
[24]	SAITO R,SMOOT M E,ONO K,et al.A travel guide to Cytoscape plugins[J].Nat Methods,2012,9(11): 1069-1076.
[25]	KAZIMIERCZYK M,KASPROWICZ M K,KASPRZYK M E,et al.Human long noncoding RNA interactome:detection,characterization and function[J].Int J Mol Sci,2020,21(3):1027.
[26]	JOSHI M,RAJENDER S.Long non-coding RNAs (lncRNAs) in spermatogenesis and male infertility[J].Reprod Biol Endocrinol,2020,18(1):103.
[27]	SHANG P,WANG Z X,CHAMBA Y,et al.A comparison of prenatal muscle transcriptome and proteome profiles between pigs with divergent growth phenotypes[J].J Cell Biochem,2019,120(4):5277-5286.
[28]	TSAI M C,MANOR O,WAN Y,et al.Long noncoding RNA as modular scaffold of histone modification complexes[J].Science,2010,329(5992):689-693.
[29]	CHEN G T,CHENG X F,SHI G L,et al.Transcriptome analysis reveals the effect of long intergenic noncoding RNAs on pig muscle growth and fat deposition[J].Biomed Res Int,2019,2019:2951427.
[30]	DENG S,BOTHE I,BAYLIES M K.The formin diaphanous regulates myoblast fusion through actin polymerization and Arp2/3 regulation[J].PLoS Genet,2015,11(8):e1005381.
[31]	LI R Y,LI B J,JIANG A W,et al.Exploring the lncRNAs related to skeletal muscle fiber types and meat quality traits in pigs[J].Genes (Basel),2020,11(8):883.
[32]	GOPINATH S D,WEBB A E,BRUNET A,et al.FOXO3 promotes quiescence in adult muscle stem cells during the process of self-renewal[J].Stem Cell Reports,2014,2(4):414-426.
[33]	SHAN T Z,LIU J Q,WU W C,et al.Roles of Notch signaling in adipocyte progenitor cells and mature adipocytes[J].J Cell Physiol,2017,232(6):1258-1261.
[34]	BERETTA M,BAUER M,HIRSCH E.PI3K signaling in the pathogenesis of obesity:the cause and the cure[J].Adv Biol Regul,2015,58:1-15.
[35]	YOSHIDA T,DELAFONTAINE P.Mechanisms of IGF-1-mediated regulation of skeletal muscle hypertrophy and atrophy[J].Cells,2020,9(9):1970.
[36]	SIN T K,ZHANG G H,ZHANG Z C,et al.Cancer-induced muscle wasting requires p38β MAPK activation of p300[J].Cancer Res,2021,81(4):885-897.
[37]	DA SILVA C,DURANDT C,KALLMEYER K,et al.The role of Pref-1 during adipogenic differentiation:an overview of suggested mechanisms[J].Int J Mol Sci,2020,21(11):4104.
[38]	COLLOTTA D,HULL W,MASTROCOLA R,et al.Baricitinib counteracts metaflammation,thus protecting against diet-induced metabolic abnormalities in mice[J].Mol Metab,2020,39:101009.
[39]	ADDINSALL A B,CACCIANI N,AKKAD H,et al.JAK/STAT inhibition augments soleus muscle function in a rat model of critical illness myopathy via regulation of complement C3/3R[J].J Physiol,2021,599(11):2869-2886.
[40]	YIN H D,HE H R,SHEN X X,et al.MicroRNA profiling reveals an abundant miR-200a-3p promotes skeletal muscle satellite cell development by targeting TGF-β2 and regulating the TGF-β2/SMAD signaling pathway[J].Int J Mol Sci,2020,21(9):3274.
[41]	LAHMANN I,BRÖHL D,ZYRIANOVA T,et al.Oscillations of MyoD and Hes1 proteins regulate the maintenance of activated muscle stem cells[J].Genes Dev,2019,33(9-10):524-535.
[42]	PÉREZ-MONTARELO D,RODRÍGUEZ M C,FERNÁNDEZ A,et al.Haplotypic diversity of porcine LEP and LEPR genes involved in growth and fatness regulation[J].J Appl Genet,2015,56(4):525-533.
[43]	HAJIHOSSEINLO A,JAFARI S,AJDARY M.The relationship of GH and LEP gene polymorphisms with fat-tail measurements (fat-tail dimensions) in fat-tailed Makooei breed of Iranian sheep[J].Adv Biomed Res,2015,4:172.
[44]	DARDENO T A,CHOU S H,MOON H S,et al.Leptin in human physiology and therapeutics[J].Front Neuroendocrinol,2010,31(3):377-393.
[45]	GRIMALDI B,BELLET M M,KATADA S,et al.PER2 controls lipid metabolism by direct regulation of PPARγ[J].Cell Metab,2010,12(5):509-520.