畜牧兽医学报  2022, Vol. 53 Issue (6): 1712-1722. DOI: 10.11843/j.issn.0366-6964.2022.06.006    PDF    
lncRNA-6617调控猪肌内前体脂肪细胞分化的筛选与功能研究
孟珊, 杨阳, 李睿霄, 姬梦婷, 张娜, 路畅, 蔡春波, 高鹏飞, 郭晓红, 曹果清, 李步高     
山西农业大学动物科学学院,太谷 030801
摘要:旨在探究lncRNA-6617的生物学特性及其对猪肌内前体脂肪细胞分化的影响,以及初步探究其上游调节机制,为后续深入研究猪肌内脂肪沉积调控机制提供参考依据。本研究以饲养在相同条件下的1、90、180日龄健康马身公猪(n=3)和大白公猪(n=3)为研究对象,利用RNA-seq技术筛选脂肪型猪(马身猪)和瘦肉型猪(大白猪)肌肉特异表达的lncRNAs,并通过生物信息学分析、RT-PCR及Sanger测序鉴定lncRNA-6617生物学特性及其猪组织中的表达特征。在猪肌内前体脂肪细胞中干扰lncRNA-6617,检测lncRNA-6617对肌内脂肪细胞分化能力的影响。随后,检测m6A修饰相关酶在不同品种猪中的表达差异,研究lncRNA-6617的上游调控机制。结果显示,lncRNA-6617位于猪18号染色体,RAMP3基因反义链的第3内含子区域,无蛋白编码能力,主要分布于细胞核; lncRNA-6617在猪所检测各组织中均有表达,背部皮下脂肪中表达量最高,且在马身猪肌肉组织中的表达极显著高于大白猪(P < 0.01);在猪肌内前体脂肪细胞分化过程中,lncRNA-6617表达量呈上升趋势,且在分化第5天表达量达到最高; 干扰lncRNA-6617诱导分化7 d后,分化的成熟脂肪细胞明显减少,成脂关键基因CEBPαPPARγaP2的表达均极显著下调(P < 0.01)。m6A修饰的去甲基化酶FTO在马身猪背最长肌和股二头肌中的表达极显著高于大白猪(P < 0.01),在腰大肌中的表达显著高于大白猪(P < 0.05),上游调节因子ZNF217在马身猪肌肉组织中的表达极显著高于大白猪(P < 0.01),与lncRNA-6617表达模式一致。干扰FTO后,lncRNA-6617的表达极显著下调(P < 0.01)。本研究筛选了在马身猪肌肉组织中显著高表达的lncRNA-6617,并深入研究了lncRNA-6617促进猪肌内前体脂肪细胞分化的功能及其上游调节机制,丰富了脂肪细胞生成的表观调控网络。
关键词马身猪    lncRNA-6617    m6A    肌内脂肪细胞    成脂分化    
Screening and Functional Study of lncRNA-6617 Regulating Porcine Intramuscular Preadipocytes Differentiation
MENG Shan, YANG Yang, LI Ruixiao, JI Mengting, ZHANG Na, LU Chang, CAI Chunbo, GAO Pengfei, GUO Xiaohong, CAO Guoqing, LI Bugao     
College of Animal Science, Shanxi Agricultural University, Taigu 030801, China
Abstract: The aim of this study was to investigate the biological characteristics of lncRNA-6617 and its effect on the differentiation of porcine intramuscular preadipocytes, preliminarily explore its upstream regulatory mechanism, and to provide the base for further clarifying the regulation mechanism of intramuscular fat deposition in pigs. In this study, healthy Mashen boars (n=3) and Large White boars (n=3) of 1, 90, and 180 days of age were raised under the same conditions. RNA-seq technology was used to screen muscle-specific lncRNAs in fatty-type pigs (Mashen pigs) and lean-type pigs (Large White pigs). Bioinformatics analysis, RT-PCR and Sanger sequencing were used to identify lncRNA-6617 biological characteristics and expression characteristics in pig tissues. The effect of lncRNA-6617 on the differentiation ability of intramuscular adipocytes were detected after interfering with lncRNA-6617 in porcine intramuscular preadipocytes. Subsequently, the expression differences of m6A modification-related enzymes in different breeds of pigs were detected, and the upstream regulatory mechanism of lncRNA-6617 was studied. The results showed that lncRNA-6617 was located on pig chromosome 18, the third intron region of the RAMP3 gene antisense strand, had no protein coding ability, and mainly distributed in the nucleus. lncRNA-6617 was expressed in all detected pig tissues, with the highest expression in back subcutaneous fat, and its expression was significantly higher in MS pigs muscle tissue than that in LW pigs(P < 0.01). In the differentiation process of porcine intramuscular preadipocytes, the expression level of lncRNA-6617 showed an upward trend and reached the peak on day 5. After interfering with lncRNA-6617, the number of differentiated mature adipocytes were significantly reduced, and the expression of key adipogenic genes CEBPα, PPARγ and aP2 were significantly down-regulated (P < 0.01) on day 7 of differentiation. The expression of m6A-modified demethylase FTO were significantly higher in the longissimus dorsi and biceps femoris muscle of MS pigs than that of LW pigs (P < 0.01), and were significantly higher than that of LW pigs in the psoas muscle (P < 0.05). And the expression of upstream regulator ZNF217 were significantly higher in the muscle tissue of MS pigs than that in LW pigs (P < 0.01), which was consistent with the expression pattern of lncRNA-6617. After interfering with FTO, the expression of lncRNA-6617 was significantly down-regulated (P < 0.01). In conclusion, lncRNA-6617 highly expressed in the muscle tissue of Mashen pig was screened and identified. The function of lncRNA-6617 in promoting the differentiation of porcine intramuscular preadipocytes and its upstream regulation mechanism were further studied, which riched the epigenetic regulation network of adipogenesis.
Key words: Mashen pig    lncRNA-6617    m6A    intramuscular adipocytes    adipogenic differentiation    

脂肪组织传统上被认为是机体内主要的能量储存器官[1]。随着分子研究层面的不断深入,研究发现脂肪组织也是机体最大的内分泌器官之一,通过分泌大量的脂肪细胞因子影响全身新陈代谢[2]。脂肪组织分布在动物机体内的多个部位,如皮下脂肪和内脏脂肪等,其中存在于肌纤维之间和肌纤维束之间的肌内脂肪,是决定动物肉品质的重要因素之一,影响着肉的口感和风味[3-4]。研究显示,肌内脂肪细胞来源于肌肉血管基质部分(sromal vascular fraction, SVF)中的肌内前体脂肪细胞,逐步分化形成成熟的脂肪细胞[5-6]。近年来,随着功能基因组学研究和转录组测序技术的应用,越来越多的调控因子被报道参与脂肪生成的调控。长链非编码RNA(long noncoding RNA, lncRNA)是一类转录本长度大于200 nt的RNA分子,通过影响下游关键因子的表达在脂肪细胞生成过程中发挥着重要的调控作用[7-10]。但关于lncRNA调控肌内脂肪生成的研究还较少,其调控网络尚不够完善。目前,关于lncRNA调控脂肪细胞生成的研究主要集中于下游的分子机制,lncRNA本身是如何被调控?其调控的分子机制还尚待研究。N6-甲基腺苷(N6-methyladenosine,m6A)是一类在真核生物RNA上广泛存在的表观遗传修饰[11],研究发现,lncRNA本身存在着m6A修饰,m6A修饰水平的变化影响着lncRNA功能的发挥,从而参与细胞增殖、纤维化、癌症等的调控[12-14]。2014年,首次发现了去甲基化酶FTO通过m6A修饰的靶mRNA的选择性剪切,影响脂肪生成[15]。但是m6A修饰是否可以通过介导lncRNA的作用参与成脂分化的调控还需进一步揭示。

马身猪是山西省地方品种,表现出肉质优良、繁殖力较高、抗逆性强[16],且背最长肌中肌内脂肪含量极显著高于大白猪[17]。研究发现,在马身猪肌肉组织中,脂肪细胞生成的关键调控因子Zfp423、PPARγCEBPα表达量显著升高[18],说明马身猪的肌内脂肪细胞生成相较于大白猪显著增加。因此,本研究通过对肌内脂肪含量高的马身猪和肌内脂肪含量较低的大白猪背最长肌组织进行lncRNA测序[19],从而筛选到一个在马身猪中特异高表达的lncRNA-6617。通过生物信息学软件分析其生物学特性,比较分析lncRNA-6617在马身猪和大白猪肌肉组织间的表达差异。利用体外分离培养的猪肌内前体脂肪细胞为模型,干扰lncRNA-6617,探究其对猪肌内前体脂肪细胞成脂分化的影响。为了寻找lncRNA-6617上游作用机制的调控方式,检测马身猪和大白猪肌肉组织中m6A修饰相关酶的表达差异,为后续分子作用机制的深入研究提供理论基础。

1 材料与方法 1.1 试验材料

选取饲养在相同条件下的1、90、180日龄健康马身公猪和大白公猪各3头,分别屠宰采集其心、肝、脾、肺、背部和腹部皮下脂肪、背最长肌、腰大肌、股二头肌等组织,在液氮中速冻,随后于-80 ℃冰箱保存备用。

1.2 主要试剂与仪器

1.2.1 主要试剂   胎牛血清和Opti-MEM购自美国Gibco公司; 胰酶、青链霉素、高糖DMEM和胰岛素购自美国Sigma公司; Dex、IBMX、4%多聚甲醛和饱和油红O染液购自北京Solarbio公司; 吲哚美辛购自阿拉丁; Trizol Reagen和PARISTM Kit购自美国Life Technologies公司; Prime Script RT reagent Kit with gDNA Eraser和SYBR Premix Ex Taq II等购自北京全式金; Lipofectamine 2000 Reagent购自Invitrogen公司; 引物由上海生工合成; siRNA委托上海吉玛公司设计并合成。

1.2.2 主要仪器   细胞成像系统(EVOS FL Auto,Life Technologies,美国),Nanodrop 2000核酸蛋白测定仪(Thermo,美国),普通PCR仪(Bio-RAD,美国),实时荧光定量PCR仪(Bio-RAD,美国)。

1.3 试验方法

1.3.1 RNA提取及cDNA文库的构建和测序   按照Trizol说明书方法提取1、90和180日龄马身猪和大白猪背最长肌组织RNA,经质量检测合格后,进行Illumina Hiseq上机测序。

1.3.2 猪lncRNA-6617的筛选及生物信息学分析   通过edgeR[20]和cuffdiff[21]软件对1日龄马身猪与大白猪背最长肌lncRNA测序结果分析,选定马身猪中特异高表达的lncRNA-6617为目标lncRNA。利用UCSC-Blast和Ensemble-Blast明确lncRNA-6617的基因组定位; 通过CPC[22]、CPAT[23]在线网站预测lncRNA-6617的蛋白编码能力。

1.3.3 lncRNA-6617的鉴定   根据测序得到的lncRNA-6617序列信息,用Primer Premier 6.0软件进行引物设计,以马身猪背最长肌组织cDNA为模板,通过RT-PCR对lncRNA进行扩增,随后用1%琼脂糖凝胶电泳检测扩增产物,对目的片段进行胶回收,连接转化,挑取阳性单菌落进行PCR鉴定,将筛选出的阳性菌送至华大基因测序。

1.3.4 猪lncRNA-6617的核质定位   试验选用未分化的猪肌内前体脂肪细胞,按照试剂盒说明书进行细胞核和细胞质RNA分离,把RNA反转录成cDNA,通过qRT-PCR检测lncRNA-6617在猪肌内前体脂肪细胞的细胞核和细胞质中的表达量。

1.3.5 猪肌内前体脂肪细胞的成脂分化   马身猪肌内前体脂肪细胞由实验室前期分离培养和保存,将传代细胞接种于6孔板,采用“激素鸡尾酒”法[24]诱导成脂。待细胞完全汇合后,更换为诱导培养基(完全培养基+10 μg·mL-1 Ins+0.5 mmoL·L-1 IBMX+1 μmoL·L-1 DEX+100 μmoL·L-1 IND),4 d后更换为维持培养基(完全培养基+10 μg·mL-1 Ins),每2 d换1次培养基,直至脂滴融合成大脂滴。

1.3.6 细胞转染及油红O染色   将肌内前体脂肪细胞接种于6孔板,密度达到60%左右时进行转染。转染前更换为新鲜培养基,siRNA试剂按照Lipofectamine 2000说明书步骤转染,每组设置3个重复孔。转染48 h后收集细胞,检测干扰效率。同时收集转染后诱导分化7 d的细胞,检测分化中后期成脂关键基因的表达。待细胞脂滴明显形成后,用预冷的PBS清洗两遍,4%多聚甲醛固定30 min,弃去甲醛,60%异丙醇浸洗1 min,弃去异丙醇。待细胞完全干燥,油红O染液(饱和油红O∶蒸馏水= 3∶2)染色10 min,PBS洗去多余的油红O染液,显微镜下观察拍照。

1.3.7 总RNA的提取与实时定量PCR   按照Trizol说明书方法提取RNA。按照全式金反转录试剂盒说明书合成第一链cDNA,稀释为20倍。然后以该cDNA作为模板,通过qRT-PCR检测基因表达情况,β-actin为内参基因,引物序列见表 1。反应体系为10 μL:SYBR为5 μL,cDNA为4.4 μL,上、下游引物各为0.3 μL。反应程序:95 ℃预变性3 min; 95 ℃变性10 s,60 ℃退火20 s,72 ℃延伸30 s,40个循环。基因相对表达量根据2-ΔΔCt法计算。

表 1 引物信息 Table 1 Primer information

1.3.8 数据统计与分析   试验均设置3个生物学重复,所有数据结果使用SPSS Statistics 22.0软件one-way ANOVA和独立样本T检验进行显著性分析,采用Duncan′s法进行多重比较,P < 0.05表示差异显著,P < 0.01表示差异极显著。

2 结果 2.1 lncRNA-6617的筛选

利用edgeR和cuffdiff软件分析不同发育阶段马身猪和大白猪背最长肌组织中差异表达的lncRNA,结果显示,1日龄马身猪与90日龄马身猪、1日龄马身猪与180日龄马身猪、90日龄马身猪与180日龄马身猪相比,分别有88、60、18个差异表达的lncRNAs(图 1A1B),1日龄马身猪与1日龄大白猪、90日龄马身猪与90日龄大白猪、180日龄马身猪与180日龄大白猪相比,分别有87、21、16个差异表达的lncRNAs(图 1C1D)。

MS. 马身猪; LW. 大白猪。下同 MS. Mashen pig; LW. Large White pig. The same as below 图 1 马身猪和大白猪不同阶段差异表达的lncRNA数量 Fig. 1 The number of differentially expressed lncRNAs at different stages of MS and LW pigs

通过比较分析1日龄马身猪与大白猪背最长肌差异lncRNA的表达量,筛选到在马身猪中特异高表达的lncRNA-6617(图 2)。

图 2 1日龄马身猪和大白猪差异lncRNA热图(A)和火山图(B) Fig. 2 The difference lncRNAs heat map (A) and volcano map (B) of 1-day-old MS and LW pigs
2.2 lncRNA-6617的生物学鉴定

利用测序得到的lncRNA-6617序列放入UCSC和Ensemble数据库进行比对,发现lncRNA-6617位于猪第18号染色体上,由RAMP3基因所在DNA链的反义链转录而来,与RAMP3基因内含子3具有278 bp的碱基互补区域(图 3A),提示该lncRNA为反义型lncRNA。

A. 基因组位置; B. 蛋白编码能力预测:a为CPAT结果,b为CPC结果; C. 亚细胞定位; D. PCR扩增; E. 测序结果与已知序列比对 A. Genomic location; B. Protein coding ability prediction: a is the CPAT result, b is the CPC result. C. Subcellular localization; D. PCR amplification; E. Alignment of sequencing results with known sequences 图 3 lncRNA-6617的生物学特性 Fig. 3 Biological characteristics of lncRNA-6617

通过CPAT及CPC在线网站对lncRNA-6617进行蛋白编码能力预测分析,选择已报道的非编码基因NEAT1为参考,结果表明,lncRNA-6617序列不具有编码蛋白质的能力(图 3B)。利用未分化的猪肌内前体脂肪细胞进行核质分离试验,确定lncRNA-6617的亚细胞定位,结果显示,lncRNA-6617主要分布于细胞核(图 3C)。

为验证序列的存在,以转录组测序获得的序列为模板设计引物,通过RT-PCR扩增lncRNA-6617。琼脂糖凝胶电泳结果显示,扩增出长278 bp的目的条带(图 3D),条带清晰且单一,Sanger测序结果表明条带序列比对一致(图 3E)。

2.3 lncRNA-6617的表达模式分析

为明确lncRNA-6617组织表达特性,通过qRT-PCR检测了lncRNA-6617在马身猪心、肝、脾、肺、背部和腹部皮下脂肪、背最长肌、腰大肌、股二头肌等9个组织中的表达情况。结果显示,lncRNA-6617在各组织中均有表达(图 4A),背部皮下脂肪中表达量最高,其次是腹部皮下脂肪等组织。进一步发现,lncRNA-6617的表达量在肌内脂肪含量高的马身猪肌肉组织中极显著高于肌内脂肪含量较低的大白猪(P < 0.01,图 4B)。

A. 组织表达谱,不同小写字母表示差异显著(P < 0.05);B. 马身猪和大白猪肌肉组织间的比较; C、D. 诱导肌内前体脂肪细胞分化过程中PPARγ和lncRNA-6617的时序表达,**. P < 0.01,*. P < 0.05,下同 A. Tissue expression profile, different lowercase letters indicate significant difference(P < 0.05); B. Comparison of muscle tissues between Mashen pigs and large White pigs; C, D. Temporal expression during induction of intramuscular preadipocytes differentiation, **. P < 0.01, *. P < 0.05, the same as below 图 4 lncRNA-6617表达模式 Fig. 4 lncRNA-6617 expression patterns

对猪肌内前体脂肪细胞进行成脂诱导分化,分别收集诱导0、1、3、5、7 d的细胞,qRT-PCR检测PPARγ和lncRNA-6617在成脂诱导分化中的表达情况。结果表明,随着成脂分化进程的进行,lncRNA-6617表达量持续升高,与成脂核心调控因子PPARγ的趋势一致(图 4C4D),相较于未分化阶段,该lncRNA主要集中在分化中后期表达。

2.4 干扰lncRNA-6617抑制猪肌内前体脂肪细胞成脂分化

将siRNA和对照组silnc-NC分别转染至猪肌内前体脂肪细胞中,48 h后收集细胞检测干扰效率。结果显示,与silnc-NC组相比,silnc-6617组极显著降低了lncRNA-6617的表达(P < 0.01,图 5A)。干扰lncRNA-6617后,对猪肌内前体脂肪细胞进行诱导分化7 d后,分化的成熟脂肪细胞明显减少(图 5C),成脂关键基因CEBPαPPARγaP2的表达均极显著下调(P < 0.01,图 5B),表明干扰lncRNA-6617抑制猪肌内前体脂肪细胞分化。

A. 干扰效率; B. 干扰lncRNA-6617后分化7 d成脂相关基因的表达变化; C. 油红O染色(100×) A. Interference efficiency; B. The expression changes of adipogenesis-related genes after interfering with lncRNA-6617 on day 7 of diferentiation; C. Oil Red O staining (100×) 图 5 干扰lncRNA-6617抑制猪肌内前体脂肪细胞成脂分化 Fig. 5 Interfering with lncRNA-6617 inhibits adipogenic differentiation of pig intramuscular preadipocytes
2.5 m6A去甲基化酶FTO调控lncRNA-6617的表达

为了探究lncRNA-6617上游的调控机制,本研究分析了马身猪和大白猪肌肉组织m6A修饰水平的差异,结果显示,m6A去甲基化酶FTO的表达水平在马身猪背最长肌和股二头肌中极显著高于大白猪(P < 0.01),在马身猪腰大肌中显著高于大白猪(P < 0.05),上游调节因子ZNF217的表达水平在马身猪各个肌肉组织中均极显著高于大白猪(P < 0.01,图 6)。

图 6 马身猪和大白猪肌肉组织中m6A修饰相关酶的表达 Fig. 6 The expression of m6A modification related enzymes in muscle tissues of MS and LW pigs

基于SRAMP数据库,利用lncRNA-6617序列预测m6A修饰位点,发现在lncRNA-6617上存在m6A潜在的修饰位点GGACU/A(图 7A)。在猪肌内前体脂肪细胞中干扰FTO,发现干扰FTO极显著降低了lncRNA-6617的表达(P < 0.01,图 7B),说明m6A修饰可能介导lncRNA-6617的表达。

A. lncRNA-6617上m6A修饰位点预测; B. FTO干扰效率以及干扰FTO后lncRNA-6617的表达变化 A. Prediction of m6A modification site on lncRNA-6617; B. The expression changes of FTO and lncRNA-6617 after interfering with FTO 图 7 FTO可能调控lncRNA-6617的表达 Fig. 7 FTO may regulate the expression of lncRNA-6617
3 讨论

lncRNA是一类转录本长度大于200 nt的RNA分子,通常不编码蛋白质,而是以RNA的形式在表观遗传修饰、转录调控以及转录后修饰等生物学事件中参与蛋白编码基因调控[25-27]。随着高通量测序的发展,一系列调控脂肪细胞生成的lncRNA相继被鉴定,其发挥作用的机制也不断被揭示。lncRNA通过多种作用方式影响关键转录因子CEBPαPPARγ等调控脂肪细胞生成[28-30]。最近研究也发现,lncIMF4[31]和IMFlnc1[32]参与调控猪肌内脂肪细胞的生成。但是,目前在猪肌内脂肪细胞生成的研究中,更多侧重于测序数据的分析,对于lncRNA的调控成脂分化功能研究尚不深入。

由于lncRNA具有明显的时空表达特异性[33],使得研究肌内脂肪细胞生成需要进一步筛选鉴定特异表达的lncRNA。此前有研究者分别对绒山羊[34]和肉用山羊[35]肌内脂肪细胞成脂分化前后的细胞进行高通量测序,筛选与肌内沉积相关的功能候选基因。李嫒等[36]利用RNA-Seq技术对大白猪和莱芜猪肌内脂肪中circRNA进行鉴定,以筛选出调控猪脂肪沉积的circRNA,证明高通量测序是筛选功能基因的重要手段。本研究利用本课题组对肌内脂肪含量高的马身猪和肌内脂肪含量较低的大白猪背最长肌组织lncRNA测序结果,筛选到在马身猪中特异高表达的lncRNA-6617。进一步通过生物信息学分析和功能鉴定,发现lncRNA-6617在马身猪肌肉组织中显著高表达,说明lncRNA-6617可能是影响肌内脂肪沉积差异的一个关键因子。同时发现lncRNA-6617在猪肌内前体脂肪细胞成脂分化早期表达量逐渐上升,干扰lncRNA-6617后,显著抑制肌内脂肪细胞的生成,说明lncRNA-6617是调控猪肌内脂肪细胞生成的关键因子。lncRNA-6617功能的发现可以进一步在理论上完善肌内脂肪细胞生成的调控网络,为解析肌内脂肪组织发育机制和遗传选育筛选靶标的确立提供理论依据。

作为遗传信息的传递载体,RNA上的修饰对其功能的发挥起着至关重要的作用[37]。前期已有研究对金华猪、长白猪及长金猪的脂肪组织与mRNA m6A修饰水平进行相关性分析,发现m6A修饰酶FTO在肉脂型金华猪脂肪组织中表达量最高[38]。同样有研究发现,肌内脂肪含量高的金华猪背最长肌中m6A修饰水平显著低于含量较低的长白猪[39]。本研究中,FTO的表达在肌内脂肪含量高的马身猪肌肉组织中显著高于肌内脂肪含量较低的大白猪,且与lncRNA-6617的组织表达模式一致,提示马身猪和大白猪肌内脂肪含量的差异可能与m6A修饰相关。研究发现,在猪肌内前体脂肪细胞中,过表达甲基转移酶METTL3显著抑制了脂肪细胞的生成[40]。m6A修饰作为RNA表观修饰方式之一,同样影响着lncRNA功能的发挥。于是本研究在肌内前体脂肪细胞中干扰FTO,发现lncRNA-6617的表达显著下调,说明m6A修饰的去甲基化酶FTO可能介导lncRNA-6617功能的发挥,但其具体作用机制还需进一步研究。

4 结论

本研究鉴定了在马身猪肌肉组织中特异高表达的lncRNA-6617;干扰lncRNA-6617后,显著抑制猪肌内前体脂肪细胞分化; lncRNA-6617的表达也可能受到上游m6A去甲基化酶FTO的调控; 这些结果说明,m6A修饰可能介导lncRNA-6617参与肌内脂肪细胞生成的调控,为后续深入研究猪肌内前体脂肪细胞生成的分子机制提供了理论基础。

参考文献
[1]
LIANG W, YE D D. The potential of adipokines as biomarkers and therapeutic agents for vascular complications in type 2 diabetes mellitus[J]. Cytokine Growth Factor Rev, 2019, 48: 32-39. DOI:10.1016/j.cytogfr.2019.06.002
[2]
BRADLEY D, SHANTARAM D, SMITH A, et al. Adipose tissue T regulatory cells: implications for health and disease[J]. Adv Exp Med Biol, 2021, 1278: 125-139.
[3]
HOCQUETTE J F, GONDRET F, BAÉZA E, et al. Intramuscular fat content in meat-producing animals: development, genetic and nutritional control, and identification of putative markers[J]. Animal, 2010, 4(2): 303-319. DOI:10.1017/S1751731109991091
[4]
张宁芳, 成志敏, 乐宝玉, 等. 猪肌源性前体脂肪细胞的分离培养和成脂诱导分化研究[J]. 畜牧兽医学报, 2018, 49(12): 2612-2621.
ZHANG N F, CHENG Z M, LE B Y, et al. Isolation, culture and adipogenic differention of pig myogenic preadipocytes cell[J]. Acta Veterinaria et Zootechnica Sinica, 2018, 49(12): 2612-2621. DOI:10.11843/j.issn.0366-6964.2018.12.010 (in Chinese)
[5]
BERRY D C, STENESEN D, ZEVE D, et al. The developmental origins of adipose tissue[J]. Development, 2013, 140(19): 3939-3949. DOI:10.1242/dev.080549
[6]
RODRIGUEZ A M, ELABD C, DELTEIL F, et al. Adipocyte differentiation of multipotent cells established from human adipose tissue[J]. Biochem Biophys Res Commun, 2004, 315(2): 255-263. DOI:10.1016/j.bbrc.2004.01.053
[7]
YANG Q, WAN Q, ZHANG L T, et al. Analysis of LncRNA expression in cell differentiation[J]. RNA Biol, 2018, 15(3): 413-422. DOI:10.1080/15476286.2018.1441665
[8]
REY F, URRATA V, GILARDINI L, et al. Role of long non-coding RNAs in adipogenesis: state of the art and implications in obesity and obesity-associated diseases[J]. Obes Rev, 2021, 22(7): e13203.
[9]
SQUILLARO T, PELUSO G, GALDERISI U, et al. Long non-coding RNAs in regulation of adipogenesis and adipose tissue function[J]. Elife, 2020, 9: e59053. DOI:10.7554/eLife.59053
[10]
GHAFOURI-FARD S, TAHERI M. The expression profile and role of non-coding RNAs in obesity[J]. Eur J Pharmacol, 2021, 892: 173809. DOI:10.1016/j.ejphar.2020.173809
[11]
MEYER K D, SALETORE Y, ZUMBO P, et al. Comprehensive analysis of mRNA methylation reveals enrichment in 3'UTRs and near stop codons[J]. Cell, 2012, 149(7): 1635-1646. DOI:10.1016/j.cell.2012.05.003
[12]
WARDA A S, KRETSCHMER J, HACKERT P, et al. Human METTL16 is a N6-methyladenosine (m6A) methyltransferase that targets pre-mRNAs and various non-coding RNAs[J]. EMBO Rep, 2017, 18(11): 2004-2014. DOI:10.15252/embr.201744940
[13]
LIU H M, XU Y X, YAO B, et al. A novel N6-methyladenosine (m6A)-dependent fate decision for the lncRNA THOR[J]. Cell Death Dis, 2020, 11(8): 613. DOI:10.1038/s41419-020-02833-y
[14]
LIU P H, ZHANG B, CHEN Z, et al. m6A-induced lncRNA MALAT1 aggravates renal fibrogenesis in obstructive nephropathy through the miR-145/FAK pathway[J]. Aging (Albany NY), 2020, 12(6): 5280-5299.
[15]
ZHAO X, YANG Y, SUN B F, et al. FTO-dependent demethylation of N6-methyladenosine regulates mRNA splicing and is required for adipogenesis[J]. Cell Res, 2014, 24(12): 1403-1419. DOI:10.1038/cr.2014.151
[16]
GAO P F, CHENG Z M, LI M, et al. Selection of candidate genes affecting meat quality and preliminary exploration of related molecular mechanisms in the Mashen pig[J]. Asian-Australas J Anim Sci, 2019, 32(8): 1084-1094. DOI:10.5713/ajas.18.0718
[17]
GUO X H, QIN B Y, YANG X F, et al. Comparison of carcass traits, meat quality and expressions of MyHCs in muscles between Mashen and Large White pigs[J]. Ital J Anim Sci, 2019, 18(1): 1410-1418. DOI:10.1080/1828051X.2019.1674701
[18]
ZHAO J X, LI K, YANG Q Y, et al. Enhanced adipogenesis in Mashen pigs compared with Large White pigs[J]. Ital J Anim Sci, 2017, 16(2): 217-225. DOI:10.1080/1828051X.2017.1285682
[19]
GAO P F, GUO X H, DU M, et al. LncRNA profiling of skeletal muscles in Large White pigs and Mashen pigs during development,[J]. J Anim Sci, 2017, 95(10): 4239-4250. DOI:10.2527/jas2016.1297
[20]
ROBINSON M D, MCCARTHY D J, SMYTH G K. edgeR: a bioconductor package for differential expression analysis of digital gene expression data[J]. Bioinformatics, 2010, 26(1): 139-140. DOI:10.1093/bioinformatics/btp616
[21]
GHOSH S, CHAN C K K. Analysis of RNA-Seq data using TopHat and cufflinks[J]. Methods Mol Biol, 2016, 1374: 339-361.
[22]
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(S2): W345-W349.
[23]
WANG L G, PARK H J, DASARI S, et al. CPAT: coding-potential assessment tool using an alignment-free logistic regression model[J]. Nucleic Acids Res, 2013, 41(6): e74. DOI:10.1093/nar/gkt006
[24]
RUIZ-OJEDA F J, RUPÉREZ A I, GOMEZ-LLORENTE C, et al. Cell models and their application for studying adipogenic differentiation in relation to obesity: a review[J]. Int J Mol Sci, 2016, 17(7): 1040. DOI:10.3390/ijms17071040
[25]
SUN L, GOFF L A, TRAPNELL C, et al. Long noncoding RNAs regulate adipogenesis[J]. Proc Natl Acad Sci U S A, 2013, 110(9): 3387-3392. DOI:10.1073/pnas.1222643110
[26]
JANDURA A, KRAUSE H M. The new RNA world: growing evidence for long noncoding RNA functionality[J]. Trends Genet, 2017, 33(10): 665-676. DOI:10.1016/j.tig.2017.08.002
[27]
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. DOI:10.1038/s41580-020-00315-9
[28]
ZHU E D, ZHANG J J, LI Y C, et al. Long noncoding RNA Plnc1 controls adipocyte differentiation by regulating peroxisome proliferator-activated receptor γ[J]. FASEB J, 2019, 33(2): 2396-2408. DOI:10.1096/fj.201800739RRR
[29]
ZUO C Q, PAN Y Q, LENG D, et al. Transcriptome analysis of long non-coding RNAs reveals NR_015556 lncRNA is a novel regulator for adipocyte differentiation[J]. Biochem Biophys Res Commun, 2022, 601: 79-85. DOI:10.1016/j.bbrc.2022.02.069
[30]
CHEN C, CUI Q M, ZHANG X, et al. Long non-coding RNAs regulation in adipogenesis and lipid metabolism: emerging insights in obesity[J]. Cell Signal, 2018, 51: 47-58. DOI:10.1016/j.cellsig.2018.07.012
[31]
SUN Y M, CAI R, WANG Y Q, et al. A newly identified LncRNA LncIMF4 controls adipogenesis of porcine intramuscular preadipocyte through attenuating autophagy to inhibit lipolysis[J]. Animals (Basel), 2020, 10(6): 926.
[32]
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. DOI:10.1186/s12860-020-00324-8
[33]
刘怡冰, 敖红, 邢凯, 等. 脂肪相关lncRNAs在动物中的研究进展[J]. 中国畜牧兽医, 2019, 46(7): 2021-2029.
LIU Y B, AO H, XING K, et al. Research progress on fat-related lncRNAs in animals[J]. China Animal Husbandry & Veterinary Medicine, 2019, 46(7): 2021-2029. (in Chinese)
[34]
杜琛, 付绍印, 高鸿雁, 等. 绒山羊肌内脂肪细胞成熟前后比较转录组分析[J]. 畜牧兽医学报, 2014, 45(5): 714-721.
DU C, FU S Y, GAO H Y, et al. Transcriptome analysis of intramuscular preadipocytes and matureadipocyte in cashmere goats[J]. Acta Veterinaria et Zootechnica Sinica, 2014, 45(5): 714-721. (in Chinese)
[35]
谢光杰, 王永, 许晴, 等. 简州大耳羊肌内脂肪细胞成脂分化差异表达基因的筛选与鉴定[J]. 畜牧兽医学报, 2020, 51(7): 1525-1536.
XIE G J, WANG Y, XU Q, et al. Selection and validation of the differentially expressed genes during the adipogenic differentiation of Jianzhou Da'er goat intramuscular adipocytes[J]. Acta Veterinaria et Zootechnica Sinica, 2020, 51(7): 1525-1536. (in Chinese)
[36]
李嫒, 张秀秀, 黄万龙, 等. 大白猪和莱芜猪肌内脂肪组织circRNAs的鉴定与分析[J]. 畜牧兽医学报, 2018, 49(7): 1343-1353.
LI A, ZHANG X X, HUANG W L, et al. Identification and analysis of circRNAs in intramuscular adipose tissues between Large White and Laiwu pigs[J]. Acta Veterinaria et Zootechnica Sinica, 2018, 49(7): 1343-1353. (in Chinese)
[37]
CHI K R. The RNA code comes into focus[J]. Nature, 2017, 542(7642): 503-506. DOI:10.1038/542503a
[38]
朱琳娜. FTO、METTL3基因表达对猪脂肪细胞mRNA N6-甲基腺苷水平及脂肪沉积的影响研究[D]. 杭州: 浙江大学, 2014.
ZHU L N. Effect of FTO, METTL3 gene expression on mRNA m6A methylation and lipid metabolism in porcine subcutaneous adipocytes[D]. Hangzhou: Zhejiang University, 2014. (in Chinese)
[39]
江芹. mRNA m6A修饰关键基因对猪肌内脂肪沉积的影响及机制[D]. 杭州: 浙江大学, 2018.
JIANG Q. Effects and mechanisms of key mRNA m6A genes in porcine intramuscular fat deposition[D]. Hangzhou: Zhejiang University, 2018. (in Chinese)
[40]
WANG X X, ZHU L N, CHEN J Q, et al. mRNA m6A methylation downregulates adipogenesis in porcine adipocytes[J]. Biochem Biophys Res Commun, 2015, 459(2): 201-207. DOI:10.1016/j.bbrc.2015.02.048

(编辑   郭云雁)