畜牧兽医学报  2021, Vol. 52 Issue (2): 322-330. DOI: 10.11843/j.issn.0366-6964.2021.02.005    PDF    
CA5B基因的序列特征和表达分析
郭晋1, 范新浩2, 杨亚岚2, 梁国明1, 唐中林1,2     
1. 中国农业科学院北京畜牧兽医研究所, 北京 100193;
2. 中国农业科学院农业基因研究所, 广东岭南现代农业科学与技术实验室, 深圳 518120
摘要:旨在通过分析猪CA5B(线粒体碳酸酐酶VB,carbonic anhydrase 5b,mitochondria;CA5B/CAVB/Car5b)基因的序列特征和时空表达,解析CA5B对猪组织器官发育和代谢的影响。本研究采集6头240日龄贵州猪(公、母各3头)的心、肝、脾、肺、肾、皮下脂肪、背最长肌、睾丸、卵巢,采集通城猪和长白猪27个生长发育时间点(出生前33、40、45、50、55、60、65、70、75、80、85、90、95、100、105 d和出生后0、9、20、30、40、60、80、100、120、140、160、180 d)(每个品种3头)的骨骼肌,提取上述组织器官RNA进行转录组测序,分析CA5B基因的组织表达谱以及骨骼肌27个生长发育时间点的时序表达;基于NCBI数据库中CA5B蛋白序列,进行物种间的保守性分析,构建系统进化树,分析蛋白互作网络;通过Targetscan、PicTar和miRanda等软件预测与CA5B结合的潜在miRNAs,并在猪RNA编辑数据库中分析CA5B基因潜在的RNA编辑位点。结果发现,CA5B基因在猪、人、小鼠、猴、牛、马、山羊、绵羊、狗和猫中高度保守,相对于小鼠而言,猪CA5B蛋白优先与人聚为一类。CA5B蛋白与线粒体跨膜运输蛋白、类肌球蛋白卷曲螺旋蛋白、细胞色素家族等互作;CA5B在猪的脂肪、睾丸以及卵巢中高表达,在骨骼肌生长发育过程中胚胎期的表达高于出生后阶段;靶标预测结果表明,CA5B可能受ssc-miR-22-5p调控(P < 0.05);此外,CA5B还存在3个RNA编辑位点。结果提示,CA5B可能在猪骨骼肌生长发育和能量储存利用中起重要作用。
关键词    CA5B    时空表达    序列特征    miRNA    RNA编辑    
Sequence Characteristics and Expression Analysis of CA5B Gene in Pigs
GUO Jin1, FAN Xinhao2, YANG Yalan2, LIANG Guoming1, TANG Zhonglin1,2     
1. Institute of Animal Science, Chinese Academy of Agricultural Sciences, Beijing 100193, China;
2. Guangdong Laboratory of Lingnan Modern Agriculture, Agricultural Genomics Institute at Shenzhen, Chinese Academy of Agricultural Sciences, Shenzhen 518120, China
Abstract: This study aimed to analyze the effects of CA5B (carbonic anhydrase 5b, mitochondria; CA5B/CAVB/Car5b) on the organ/tissue development and metabolism in pigs through investigating its sequence characteristics and spatial expression patterns. The heart, liver, spleen, lung, kidney, subcutaneous fat, longissimus dorsi, testis and ovary tissues of 3 boars and 3 sows of 240-day-old Guizhou pigs were collected, and the skeletal muscle tissues at 27 growth and development time points (33, 40, 45, 50, 55, 60, 65, 70, 75, 80, 85, 90, 95, 100, 105 days before birth and 0, 9, 20, 30, 40, 60, 80, 100, 120, 140, 160, 180 days after birth) of 3 Tongcheng pigs and 3 Landrace pigs were also collected. We performed transcription sequencing using RNA extracted from the above tissues, and analyzed the expression profile of CA5B gene in 9 tissues and it's temporal expression in skeletal muscle at 27 development time points. Then the conservative analysis among species, construction of phylogenetic tree and protein interaction network analysis of CA5B gene were conducted based on its protein sequence extracted from NCBI database. Finally, three softwares (Targetscan, PicTar, miRanda) were used to predict the potential miRNAs binded to CA5B gene. The RNA editing sites located in CA5B gene were also detected using the RNA editing database. The results showed that CA5B gene was highly conservative expressed in pig, human, mouse, monkey, cattle, horse, goat, sheep, dog and cat, moreover, the pig and human had higher homology than mouse based on the protein sequence of CA5B. Interaction network showed that CA5B protein mainly interacted with mitochondrial transmembrane transport protein, myosin-like coiled-coil protein, cytochrome family and other proteins. The CA5B was highly expressed in fat, testis and ovary, and the expression level in embryonic period was higher than that in postnatal period in skeletal muscle. In addition, target prediction results suggested that CA5B might be regulated by ssc-miR-22-5p (P < 0.05), and there were 3 RNA editing sites in CA5B. The results suggested that CA5B potentially play an important role in skeletal muscle growth, development, energy storage and utilization in pig.
Key words: pig    CA5B    spatiotemporal expression    sequence characteristics    miRNA    RNA editing    

碳酸酐酶(carbonic anhydrase,CAs)是一类可催化二氧化碳和碳酸氢盐的可逆水合金属酶[1],在不同物种的组织和细胞中表达水平不同[2]。由于CAs抑制剂与糖异生、尿素生成、脂肪生成、肿瘤发生等多种生理反应有关[3-4],CAs抑制剂现已成为肥胖和肿瘤等疾病治疗的药物[5-6]。线粒体(mitochrodria, MT)通过它自身表达的DNA、细胞质中核基因编码的蛋白质和酶完成细胞内的生物氧化过程[7]。研究发现,线粒体CA调节葡萄糖的氧化代谢、ROS的产生和氧化应激、线粒体的生物发生[8]。线粒体碳酸酐酶两种同工酶CA5A(carbonic anhydrase,mitochondrial,5A;CA5A)和CA5B(carbonic anhydrase,mitochondrial,5B;CA5B)(以前被称为CAVA和CAVB)在线粒体丙酮酸代谢中发挥关键作用,其通过可逆水合二氧化碳产生HCO3-:CO2+ H2O $\overline {\leftharpoondown}$H2CO3 $\overline {\leftharpoondown}$H++HCO3-[9],为丙酮酸转化为草酰乙酸提供HCO3-,使草酰乙酸进入三羧酸循环途径,是葡萄糖有氧代谢过程的关键底物。

人源碳酸酐酶(hCAs)参与大量的生理和病理过程,包括糖异生、脂肪生成、尿素生成、致瘤性以及各种病原体的生长和毒性[6, 10]CA5(以前被称为CAV)在糖异生、脂肪生成[6]、代谢调控[8, 11]和疾病发生[6, 8, 11-12]中发挥关键作用。目前,这些hCAs已作为药物靶点被深入研究,几种hCA抑制剂已在临床应用[6]。小鼠线粒体碳酸酐酶CA5B与CA5A具有同源性[13],小鼠CA5B基因在成年皮下脂肪和生殖器脂肪中表达量最高,小鼠胰岛细胞中CA5在功能上可能与调节胰岛素分泌有关[14],此外,小鼠CA5参与调节脂肪细胞的丙酮酸羧化[15]和丙酮酸糖异生[16]。双敲除(double knockout,DKO)CA5A和CA5B使小鼠高氨血症增强和空腹血糖降低,表明CA5A和CA5B有助于尿素生成和调节糖异生[17]

猪是人类最重要的蛋白来源和生物医学模式动物之一。碳水化合物和脂肪是参与能量代谢的主要物质,因此,研究分析猪代谢相关基因对于肉品质的改良有重要作用。目前,猪CA5B基因还没有注释,其表达及调控机制尚不清楚。本研究利用猪转录组数据分析CA5B基因在发育过程中的表达特征以及组织中的表达情况,为研究CA5B的功能提供依据。

1 材料与方法 1.1 试验材料

采集240日龄贵州小型猪(公、母各3头)的心、肝、脾、肺、肾、皮下脂肪、背最长肌、睾丸、卵巢组织。此外,采集通城猪和长白猪各27个生长发育时间点(出生前33、40、45、50、55、60、65、70、75、80、85、90、95、100、105 d和出生后0、9、20、30、40、60、80、100、120、140、160、180 d)的骨骼肌,每个品种每个时间点采集3头猪。上述组织样品立即现场液氮快速冻存,然后转移至-80 ℃冰箱保存备用。

1.2 试验设备及试剂

台式低温离心机(Ependorff)、样品研磨仪(Tekmar)、制冰机(SANYO)、1.5 mL离心管(AXY -GEN)、总RNA提取试剂盒(天根)。

1.3 试验方法

1.3.1 RNA提取   取低温保存的猪上述各种组织样本提取总RNA,提取流程参照天根试剂盒说明书,提取的总RNA保存在-80 ℃冰箱备用。

1.3.2 RNA转录组测序   取冻存的RNA,根据Illumina® TruSeqTM RNA样品制备流程,构建cDNA文库,利用Illumina HiSeqTM2500,采用双端测序方法进行转录组测序,测序数据已上传至NCBI数据库(GSE73763和PRJNA488311)。

1.3.3 RNA测序数据处理   用TopHat2[18]将得到的RNA测序数据比对到猪的参考基因组上,经过注释和计算,得到CA5B基因在各个样本中RPKM(每万reads中来自某基因每千碱基长度的reads数)的表达情况。

1.3.4 CA5B功能预测   从NCBI数据库下载猪、人、小鼠、猴、牛、马、山羊、绵羊、狗和猫的CA5B基因的蛋白质序列,利用EBI的多序列比对软件Clustal Omega(https://www.ebi.ac.uk/Tools/msa/clustalo/)对不同物种的CA5B进行多序列比对分析,利用STRING数据库(https://string-db.org)分析CA5B与其它蛋白的相互作用,利用DREP[19]数据库(http://www.rnanet.org/editing/)分析CA5B潜在的RNA编辑位点,利用Targetscan、PicTar和miRanda工具预测CA5B基因的3′UTR区域与miRNA结合位点,用R语言的cor. test()函数计算CA5B基因和miRNA表达量之间的相关系数。

2 结果 2.1 CA5B的序列保守性分析

由于小鼠和人CA5B基因功能已有报道研究,为了探究猪CA5B基因功能,首先分析了猪、小鼠与人CA5B基因的保守性。从NCBI数据库下载猪(XP_020935797.1)、人(EAW98898.1)、小鼠(NP_851832.2)、猴、牛、马、山羊、绵羊、狗和猫的CA5B蛋白序列,利用EBI的多序列比对软件Clustal Omega (https://www.ebi.ac.uk/Tools/msa/clustalo/) 进行多序列比对分析。图 1为10个物种蛋白编码序列的比对结果,在人、猪、小鼠、猴、牛、马、山羊、绵羊、狗和猫中,除小鼠外,其余序列相似性均在90%左右;就猪、人和小鼠的比对结果来看,猪和小鼠CA5B序列一致性为86.4%,而猪和人同源性为90.5%。说明猪CA5B蛋白序列与人相似度更近,与小鼠次之。

图 1 不同物种CA5B基因的氨基酸序列比对结果 Fig. 1 Alignment of CA5B amino acid sequence among different species
2.2 CA5B的系统发育树

为研究CA5B基因在不同物种中的进化关系,利用MEG5软件,以猪、人、小鼠、猴、牛、马、山羊、绵羊、狗和猫的CA5B的蛋白序列构建系统发育树。结果发现,猪、狗与猫CA5B聚为一支,牛、绵羊与山羊CA5B聚为一支,这两支属于同一大分支,与人、猴CA5B具有相同的来源。这8个物种与马聚为一大支,最后与小鼠CA5B聚为一支(图 2)。表明,与小鼠相比较,猪与人CA5B在进化上更近。

图 2 CA5B序列的系统发育树 Fig. 2 Phylogenetic tree for CA5B sequence
2.3 CA5B的蛋白互作网络分析

利用STRING数据库(https://string-db.org)分析CA5B与其它蛋白的相互作用。结果发现,与CA5B互作的蛋白包括TMEM35、CYP24A1和XPNPEP2等基因编码的蛋白质(图 3)。通过在NCBI数据库(https://www.ncbi.nlm.nih.gov)和Uniprot数据库(https://www.uniprot.org/uniprot/)中检索发现,与CA5B蛋白互作的这些基因在猪中的注释较少,在人和小鼠中有较多注释。

图 3 CA5B的蛋白互作网络 Fig. 3 Interaction network for CA5B protein

CA5B互作蛋白功能:跨膜蛋白(transmembrane protein 35A,TMEM35)在细胞的过氧化物酶体中释放的可溶性肽可调节神经元轴突的生长;雄性至死因子3(male-specific lethal 3, MSL3)具有组蛋白乙酰转移酶活性及与甲基化组蛋白结合功能,在染色质重塑和转录调节中发挥作用;溶质载体家族成员(solute carrier family 25 member 30, SLC25A30)为在小鼠和人心肌、骨骼肌中特异性表达的线粒体转运蛋白,完成线粒体内外的物质运输;富含脯氨酸的Gla 1(proline rich Gla 1, PRRG1)可编码跨膜蛋白TMEM171,可与钙离子结合;细胞色素P450(cytochrome P450, family 24, subfamily a, polypeptide 1, CYP24A1)基因编码的蛋白质作用于线粒体,其将骨化三醇降解为骨化醇,在维生素D分解和钙稳态中起关键作用;TXLNG(taxilin gamma, taxilinγ)基因可编码类肌球蛋白卷曲螺旋蛋白,在人和小鼠的细胞内囊泡运输中发挥作用,可能参与调控细胞周期和调节骨密度,该基因的表达可能引起脂多糖的上调;肽酶M20结构域(peptidase M20 domain containing 2, PM20D2)编码氨基酰化酶和氨基水解酶,催化蛋白质水解;X-脯氨酰氨肽酶2(X-prolyl aminopeptidase 2, XPNPEP2)基因编码氨基肽酶,催化邻近脯氨酸残基的N端裂解。

根据CA5B蛋白互作网络关系中互作蛋白的功能,推测CA5B可能参与调节生长和发育,线粒体内物质转运和线粒体外跨膜运输,并且可能与骨骼肌钙稳态调节、细胞蛋白质分解、脂多糖代谢调控有关。

2.4 CA5B基因的组织表达谱

为了研究CA5B在各组织器官中的表达特点,本研究利用已有的贵州小型猪心、肝、脾、肺、肾、皮下脂肪、背最长肌、睾丸、卵巢9种组织的RNA数据[20],分析CA5B在各种组织器官中的表达水平。结果发现,CA5B在成年猪的脂肪、睾丸、卵巢中高表达,在骨骼肌、心肌和肾中表达水平较低(图 4)。

图 4 CA5B基因在猪中的组织表达谱 Fig. 4 Tissue expression patterns of CA5B gene in pigs
2.5 CA5B在骨骼肌生长发育中的时序表达

为了了解CA5B在骨骼肌生长发育过程中的作用,本研究利用已有的高通量测序数据,首次分析了该基因在通城猪和长白猪27个生长发育时间点(E33-D180)中骨骼肌的动态表达。发现,CA5B的表达水平在通城猪和长白猪的胚胎期高于出生后时期,在胚胎期逐渐升高,出生后降低;并且在胚胎发育后期80 d(E80)到出生0 d(D0)表达水平较高,通城猪在胚胎期80 d(E80)时表达水平最高,而长白猪在胚胎期85 d(E85)时表达水平最高(图 5)。可以推断,在通城猪和长白猪骨骼肌生长发育过程中,CA5B基因在胚胎期有重要作用,并且在胚胎骨骼肌发育后期至出生前发挥重要功能。

图 5 CA5B在通城猪和长白猪27个生长发育时间点骨骼肌中的动态表达 Fig. 5 The dynamic expression of CA5B gene across 27 growth and development time points in Tongcheng and Landrace pigs
2.6 CA5B的RNA编辑

为探究CA5B的RNA编辑位点,在DREP[19]数据库(http://www.rnanet.org/editing/)中分析CA5B潜在的RNA编辑位点,结果在CA5B基因的内含子区域和3′UTR区域分别找到了1个和2个A-to-G类型的RNA编辑位点,这3个位点均位于重复元件Pre0_SS区域,属于SINE/tRNA家族(表 1)。

表 1 CA5B的RNA编辑位点信息 Table 1 The RNA editing sites for CA5B
2.7 CA5B与miRNA结合的生物信息分析

为找到可能调节CA5B的miRNA,本研究利用Targetscan、PicTar和miRanda工具预测CA5B基因的3′UTR区域与miRNA的结合位点,结果发现,ssc-miR-22-5p可以与CA5B基因3′UTR靶向结合(图 6)。通过计算CA5B基因和miRNA表达量之间的相关系数和P值,发现CA5B与ssc-miR-22-5p在猪骨骼肌生长发育过程中的表达呈显著负相关(R=-0.31, P=0.022,图 7)。表明CA5B的表达潜在受ssc-miR-22-5p调控。

图 6 CA5B与miR-22-5p靶向结合序列 Fig. 6 Targeting sequence of CA5B combined with miR-22-5p
图 7 CA5B与miR-22-5p表达相关性 Fig. 7 Expression correlation of CA5B with miR-22-5p
3 讨论

线粒体在细胞和组织代谢中发挥重要作用,CA在线粒体内不可或缺而被广泛研究。CA5B和CA5A是线粒体丙酮酸代谢的重要底物,在葡萄糖有氧代谢过程中发挥关键作用。从组织表达情况来看,小鼠CA5B和CA5A的表达呈现出不同的组织特异性,在大多数组织中都检测到CA5B,而CA5A的表达仅限于肝、骨骼肌和肾[13]。从序列上来看,CA5A基因序列在小鼠和人之间有78%相似性,CA5B基因序列在小鼠和人之间有95%相似性,CA5B的保守性更高[13]。人和小鼠CA同工酶序列之间的系统发育树表明,CA5B和CA5A都是从人和小鼠的祖先基因进化而来,并已经进化出了不同的生理作用[13]。本研究将小鼠CA5B的蛋白结构域序列与人和猪报道的CA5B结构域序列进行比较发现,CA5B结构域在猪和人之间的保守性(90.5%一致性)高于猪和小鼠CA5B保守性(86.4%一致性)。系统发育分析显示,人和猪CA5B在进化树上优先聚为一类,最后与小鼠聚在一起。因此,可以推断CA5B基因在物种间的保守性更高。

根据分析结果发现,CA5B蛋白互作网络中的基因编码肌球蛋白卷曲螺旋蛋白、钙稳态调节蛋白、线粒体跨膜蛋白、线粒体转运蛋白以及氨基水解蛋白。据文献报道,离子通道和代谢底物转运中,人线粒体溶质载体SLC家族与CA5B蛋白表达相关[21],推测CA5B在猪中可能也参与调控线粒体内物质转运过程,在线粒体内能量代谢中发挥重要作用。根据CA5B高表达于猪的脂肪、睾丸和卵巢组织,这些都是ATP合成和利用的重要组织器官。CA5B基因在小鼠成年皮下脂肪和生殖器脂肪表达量最高,其参与调节脂肪丙酮酸羧化[15]、糖异生和尿素生成[17]。因此,推测CA5B在猪中可能也参与调控糖和脂肪的代谢。在猪骨骼肌27个生长发育时间点中,CA5B基因胚胎期表达水平高于出生后时期。推测猪的CA5B与骨骼肌生长发育以及能量储存和利用有密切关系。与此同时,还发现CA5B可能受ssc-miR-22-5p调控。miRNA对机体发育和代谢具有重要的调节作用[22-23],是靶基因表达的关键调节因子[24],参与调节骨骼肌细胞增殖[25-26]、分化[27-28]和再生[29],一些miRNA调节成脂基因的表达、脂肪的形成[30]和分化[31]、白色脂肪棕色化[32]以及棕色脂肪扩张[33]。据研究报道,miR-22-5p是脂肪来源间充质干细胞的的重要参考基因[34],还在骨溶解中发挥重要作用[35]。因而,进一步推测ssc-miR-22-5p可能通过靶向CA5B参与调节猪的体内发育和代谢。最后,本研究分析发现了CA5B基因的3个A-to-G编辑位点。RNA编辑位点可能参与调节基因表达和细胞功能[36]、骨骼肌发育[19]、疾病发生[37]和肿瘤的发展[38],可能对猪基因的表达进行调控。

猪是重要的家畜,为我们的生活提供大量的蛋白质和脂肪。因此,研究分析参与调节猪糖代谢和脂肪生成的CA5B基因对于肉品质的改良有一定的参考作用。但是,现在对猪CA5B的研究较少,因此需要对该基因进行深入的挖掘以及功能验证。

4 结论

本研究初步探讨了CA5B基因在10个物种间的保守性和进化特征,CA5B蛋白可能参与细胞内线粒体内外ATP的合成、物质转运以及体内氨基酸的分解代谢。ssc-miR-22-5p和3个RNA编辑位点潜在参与调控CA5B基因表达。CA5B可能在猪胚胎期的发育过程中发挥重要功能,在猪出生后能量储存和利用的组织中有关键作用。以上结果为更深入研究猪CA5B的功能提供了重要的信息和奠定了基础。

致谢: 感谢北京畜牧兽医研究所猪基因工程与种质创新团队提供的科研平台,感谢实验室师兄、师姐、师弟、师妹在样品采集和数据分析时提供的帮助和支持,感谢团队所有人员的照顾和帮助。

参考文献
[1] HASSAN M I, SHAJEE B, WAHEED A, et al. Structure, function and applications of carbonic anhy-drase isozymes[J]. Bioorg Med Chem, 2013, 21(6): 1570–1582. DOI: 10.1016/j.bmc.2012.04.044
[2] MCKENNA R, FROST S C. Overview of the carbonic anhy-drase family[J]. Subcell Biochem, 2014, 75: 3–5.
[3] DODGSON S J, FORSTER Ⅱ R E. Carbonic anhydrase: inhibition results in decreased urea pro-duction by hepatocytes[J]. J Appl Physiol (1985), 1986, 60(2): 646–652. DOI: 10.1152/jappl.1986.60.2.646
[4] RASMUSSEN J K, BOEDTKJER E. Carbonic anhydrase inhibitors modify intracellular pH tran-sients and contractions of rat middle cerebral arteries during CO2/HCO3- fluctuations[J]. J Cereb Blood Flow Metab, 2018, 38(3): 492–505. DOI: 10.1177/0271678X17699224
[5] RAHMAN S, BIBI S, JAVED T, et al. Review: therapeutic potential of carbonic anhydrase inhibitors[J]. Pak J Pharm Sci, 2019, 32(2): 709–720.
[6] COSTA G, CARTA F, AMBROSIO F A, et al. A computer-assisted discovery of novel potential anti-obesity compounds as selective carbonic anhydrase VA inhibitors[J]. Eur J Med Chem, 2019, 181: 111565. DOI: 10.1016/j.ejmech.2019.111565
[7] VAN DER BLIEK A M, SEDENSKY M M, MORGAN P G. Cell biology of the mitochondrion[J]. Genetics, 2017, 207(3): 843–871. DOI: 10.1534/genetics.117.300262
[8] BERNARDINO R L, DIAS T R, MOREIRA B P, et al. Carbonic anhydrases are involved in mitochondrial biogenesis and control the production of lactate by human Sertoli cells[J]. FEBS J, 2019, 286(7): 1393–1406. DOI: 10.1111/febs.14779
[9] PRICE T O, ERANKI V, BANKS W A, et al. Topiramate treatment protects blood-brain barrier pericytes from hyperglycemia-induced oxidative damage in diabetic mice[J]. Endocrinology, 2012, 153(1): 362–372. DOI: 10.1210/en.2011-1638
[10] LI Z, JIANG L, CHEW S H, et al. Carbonic anhydrase 9 confers resistance to ferroptosis/apoptosis in malignant mesothelioma under hypoxia[J]. Redox Biol, 2019, 26: 101297. DOI: 10.1016/j.redox.2019.101297
[11] VAN KARNEBEEK C D, SLY W S, ROSS C J, et al. Mitochondrial carbonic anhydrase VA deficiency resulting from CA5A alterations presents with hyperammonemia in early childhood[J]. Am J Hum Genet, 2014, 94(3): 453–461. DOI: 10.1016/j.ajhg.2014.01.006
[12] SOLESIO M E, PEIXOTO P M, DEBURE L, et al. Carbonic anhydrase inhibition selectively prevents amyloid β neurovascular mitochondrial toxicity[J]. Aging Cell, 2018, 17(4): e12787. DOI: 10.1111/acel.12787
[13] SHAH G N, HEWETT-EMMETT D, GRUBB J H, et al. Mitochondrial carbonic anhydrase CA VB: differences in tissue distribution and pattern of evolution from those of CA VA suggest distinct physiological roles[J]. Proc Natl Acad Sci U S A, 2000, 97(4): 1677–1682. DOI: 10.1073/pnas.97.4.1677
[14] PARKKILA A K, SCARIM A L, PARKKILA S, et al. Expression of carbonic anhydrase V in pancreatic beta cells suggests role for mitochondrial carbonic anhydrase in insulin secretion[J]. J Biol Chem, 1998, 273(38): 24620–24623. DOI: 10.1074/jbc.273.38.24620
[15] HAZEN S A, WAHEED A, SLY W S, et al. Differentiation-dependent expression of CA V and the role of carbonic anhydrase isozymes in pyruvate carboxylation in adipocytes[J]. FASEB J, 1996, 10(4): 481–490. DOI: 10.1096/fasebj.10.4.8647347
[16] DODGSON S J, FORSTER R E. Inhibition of CA V decreases glucose synthesis from pyruvate[J]. Arch Biochem Biophys, 1986, 251(1): 198–204. DOI: 10.1016/0003-9861(86)90066-4
[17] SHAH G N, RUBBELKE T S, HENDIN J, et al. Targeted mutagenesis of mitochondrial carbonic anhydrases VA and VB implicates both enzymes in ammonia detoxification and glucose metabolism[J]. Proc Natl Acad Sci U S A, 2013, 110(18): 7423–7428. DOI: 10.1073/pnas.1305805110
[18] KIM D, PERTEA G, TRAPNELL C, et al. TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions[J]. Genome Biol, 2013, 14(4): R36. DOI: 10.1186/gb-2013-14-4-r36
[19] YANG Y L, ZHU M, FAN X H, et al. Developmental atlas of the RNA editome in Sus scrofa skeletal muscle[J]. DNA Res, 2019, 26(3): 261–272. DOI: 10.1093/dnares/dsz006
[20] LIANG G M, YANG Y L, NIU G L, et al. Genome-wide profiling of Sus scrofa circular RNAs across nine organs and three developmental stages[J]. DNA Res, 2017, 24(5): 523–535. DOI: 10.1093/dnares/dsx022
[21] HAMURO J, DEGUCHI H, FUJITA T, et al. Polarized expression of ion channels and solute carrier family transporters on heterogeneous cultured human corneal endothelial cells[J]. Invest Ophthalmol Vis Sci, 2020, 61(5): 47. DOI: 10.1167/iovs.61.5.47
[22] HOU X H, TANG Z L, LIU H L, et al. Discovery of MicroRNAs associated with myogenesis by deep sequencing of serial developmental skeletal muscles in pigs[J]. PLoS One, 2012, 7(12): e52123. DOI: 10.1371/journal.pone.0052123
[23] YANG Y L, LI Y, LIANG R Y, et al. Dynamic expression of MicroRNA-127 during porcine prenatal and postnatal skeletal muscle development[J]. Journal of Integr Agric, 2014, 13(6): 1331–1339. DOI: 10.1016/S2095-3119(13)60419-0
[24] LU T X, ROTHENBERG M E. MicroRNA[J]. J Allergy Clin Immunol, 2018, 141(4): 1202–1207. DOI: 10.1016/j.jaci.2017.08.034
[25] WU N Z, GU T T, LU L, et al. Roles of miRNA-1 and miRNA-133 in the proliferation and differentiation of myoblasts in duck skeletal muscle[J]. J Cell Physiol, 2019, 234(4): 3490–3499. DOI: 10.1002/jcp.26857
[26] JIANG A W, DONG C, LI B J, et al. MicroRNA-206 regulates cell proliferation by targeting G6PD in skeletal muscle[J]. FASEB J, 2019, 33(12): 14083–14094. DOI: 10.1096/fj.201900502RRRR
[27] WANG J, TAN J Y, QI Q, et al. miR-487b-3p suppresses the proliferation and differentiation of myoblasts by targeting IRS1 in skeletal muscle myogenesis[J]. Int J Biol Sci, 2018, 14(7): 760–774. DOI: 10.7150/ijbs.25052
[28] KONG D L, HE M, YANG L, et al. MiR-17 and miR-19 cooperatively promote skeletal muscle cell differentiation[J]. Cell Mol Life Sci, 2019, 76(24): 5041–5054. DOI: 10.1007/s00018-019-03165-7
[29] ZENG P, HAN W H, LI C Y, et al. miR-378 attenuates muscle regeneration by delaying satellite cell activation and differentiation in mice[J]. Acta Biochim Biophys Sin (Shanghai), 2016, 48(9): 833–839. DOI: 10.1093/abbs/gmw077
[30] GERIN I, BOMMER G T, MCCOIN C S, et al. Roles for miRNA-378/378* in adipocyte gene expression and lipogenesis[J]. Am J Physiol Endocrinol Metab, 2010, 299(2): E198–E206. DOI: 10.1152/ajpendo.00179.2010
[31] ESAU C, KANG X L, PERALTA E, et al. MicroRNA-143 regulates adipocyte differentiation[J]. J Biol Chem, 2004, 279(50): 52361–52365. DOI: 10.1074/jbc.C400438200
[32] LIU W Y, BI P P, SHAN T Z, et al. miR-133a regulates adipocyte browning in vivo[J]. PLoS Genet, 2013, 9(7): e1003626. DOI: 10.1371/journal.pgen.1003626
[33] PAN D N, MAO C X, QUATTROCHI B, et al. MicroRNA-378 controls classical brown fat expansion to counteract obesity[J]. Nat Commun, 2014, 5: 4725. DOI: 10.1038/ncomms5725
[34] RAGNI E, ORFEI C P, DE LUCA P, et al. miR-22-5p and miR-29a-5p are reliable reference genes for analyzing extracellular vesicle-associated miRNAs in adipose-derived mesenchymal stem cells and are stable under inflammatory priming mimicking osteoarthritis condition[J]. Stem Cell Rev Rep, 2019, 15(5): 743–754. DOI: 10.1007/s12015-019-09899-y
[35] YANG X C, ZHANG Y Y, LI Y S, et al. MALAT1 enhanced the proliferation of human osteoblasts treated with ultra-high molecular weight polyethylene by targeting VEGF via miR-22-5p[J]. Int J Mol Med, 2018, 41(3): 1536–1546.
[36] GATSIOU A, VLACHOGIANNIS N, LUNELLA F F, et al. Adenosine-to-Inosine RNA editing in health and disease[J]. Antioxid Redox Signal, 2018, 29(9): 846–863. DOI: 10.1089/ars.2017.7295
[37] STELLOS K, GATSIOU A, STAMATELOPOULOS K, et al. Adenosine-to-inosine RNA editing controls cathepsin S expression in atherosclerosis by enabling HuR-mediated post-transcriptional regulation[J]. Nat Med, 2016, 22(10): 1140–1150. DOI: 10.1038/nm.4172
[38] VELAZQUEZ-TORRES G, SHOSHAN E, IVAN C, et al. A-to-I miR-378a-3p editing can prevent melanoma progression via regulation of PARVA expression[J]. Nat Commun, 2018, 9(1): 461. DOI: 10.1038/s41467-018-02851-7