Cloning, expression analysis of TRIAP1 and its related miRNAs screening in Penaeus monodon
-
摘要:
凋亡抑制因子1 (tumor protein p53 regulated inhibitor of apoptosis 1, TRIAP1)可以通过抑制细胞凋亡参与生物体内应激反应。为探究TRIAP1基因及调控其表达的miRNAs在斑节对虾 (Penaeus monodon) 应对哈氏弧菌(Vibrio harveyi)刺激时的表达调控情况，该研究通过RACE技术获得了斑节对虾TRIAP1基因(PmTRIAP1)cDNA全长序列，并利用实时荧光定量PCR (qRT-PCR)技术，对PmTRIAP1在斑节对虾不同组织中的表达分布模式，以及哈氏弧菌刺激下PmTRIAP1与相关miRNAs的表达关联情况进行了探究。结果显示，PmTRIAP1 cDNA全长2 522 bp，包括11 bp的5'非编码区(UTR)、2 289 bp的3'UTR和222 bp的开放阅读框(ORF)，共编码73个氨基酸。PmTRIAP1基因在各组织中均有表达，其中在血细胞中表达量最高；哈氏弧菌刺激下，PmTRIAP1表达量显著降低，而Pm-miR-145-3p和Pm-miR-454-3p表达量显著升高，表明PmTRIAP1和Pm-miR-145-3p、Pm-miR-454-3p的表达量呈负相关关系；双荧光素酶报告实验结果显示，Pm-miR-145-3p和Pm-miR-454-3p可通过结合PmTRIAP1 3'UTR降低荧光素酶活性，表明Pm-miR-145-3p和Pm-miR-454-3p均可靶向调控PmTRIAP1的表达。
-
关键词:
- 斑节对虾 /
- TRIAP1 /
- Pm-miR-145-3p /
- Pm-miR-454-3p /
- 哈氏弧菌
Abstract:TRIAP1 can participate in vivo responses by inhibiting apoptosis. In order to analyze the expression regulation of TRIAP1 and miRNAs that regulate TRIAP1's expression under Vibrio harveyi challenge in Penaeus monodon, we obtained the P. monodon TRIAP1 gene (PmTRIAP1) cDNA by RACE, then detected the tissue distribution of PmTRIAP1 by quantitative real-time PCR (qRT-PCR), and investigated the expression association between PmTRIAP1 and its related miRNAs under V. harveyi challenge. The results show that the full-length PmTRIAP1 cDNA sequence was 2 522 bp, containing an 11 bp 5' non-coding region (UTR), a 2 289 bp 3'-UTR and a 222 bp open reading frame encoding 73 amino acids. The qRT-PCR analysis shows that PmTRIAP1 was ubiquitously expressed in all the tested tissues and the highest expression level was obtained in hemocyte. The expression of PmTRIAP1 decreased significantly; however, the expressions of Pm-miR-145-3p and Pm-miR-454-3p increased significantly under V. harveyi challenge, which suggests that the expression levels of PmTRIAP1 and Pm-miR-145-3p, Pm-miR-454-3p were negatively correlated. Dual luciferase reporter assay shows that Pm-miR-145-3p and Pm-miR-454-3p can bind PmTRIAP1 3'UTR to reduce luciferase activity. It is suggested that both Pm-miR-145-3p and Pm-miR-454-3p can target regulate the expression of PmTRIAP1.
-
Keywords:
- Penaeus monodon /
- TRIAP1 /
- Pm-miR-145-3p /
- Pm-miR-454-3p /
- Vibrio harveyi
-
肿瘤蛋白p53调控的凋亡抑制因子1 (tumor protein p53 regulated inhibitor of apoptosis 1, TRIAP1)是一类受p53靶向调控的功能型蛋白因子，在细胞凋亡等程序性细胞死亡过程中起重要作用[1]。研究表明，TRIAP1主要通过与热休克蛋白70 (HSP70)相互作用阻止凋亡酶激活因子(APAF1)和细胞色素C结合，并抑制含半胱氨酸的天冬氨酸蛋白水解酶9 (Caspase9)的激活，从而参与调控细胞凋亡过程[1]。例如，低水平的基因毒性物质刺激下，生物机体内的p53可特异性结合TRIAP1基因第二外显子上的p53结合位点诱导TRIAP1产生，进而抑制细胞凋亡并启动受损基因修复机制[1]；Fook-Alves等[2]对经特异性沉默TRIAP1基因后的多发性骨髓瘤细胞系进行研究时发现，Caspase9和APAF1基因的表达量呈显著下降趋势，晚期细胞凋亡率显著升高；Potting等[3]在Hela细胞中也获得了类似的研究结果，降低TRIAP1的表达会阻碍心磷脂(CL)累积、增加细胞色素C释放，进而引发细胞凋亡的发生。近年来，TRIAP1在高等动物凋亡相关的生理活动中的研究屡屡见诸报端，但在低等水产养殖动物中却鲜有报道。
MicroRNA (miRNA)是一类长约18~25个核苷酸的内源性非编码小RNA分子，它可以依赖于自身与靶基因2种不同机制的互补性(完全匹配或者种子序列匹配)以类似RNA干扰的途径调控目的基因翻译或降解靶mRNA和新合成的多肽链[4-6]。成熟的miRNA与靶基因的mRNA非翻译区(UTR)结合的特异性主要由miRNA 5'末端被称为种子序列的2~8个碱基决定[7-8]。现已证实，无论是像甲壳动物、棘皮动物、软体动物等的低等无脊椎水产养殖动物，还是像鱼类一样的高等脊椎水产养殖动物，miRNA作为转录后调控的关键因子，参与调控机体内多种重要的生命活动过程，如生长发育、生殖和免疫等[9-10]。近年来，诸如miR-15、miR-16和miR-29等大批与细胞凋亡发生密切相关的miRNA被陆续鉴定出来[11-12]，并已初步证实上述miRNA主要通过靶向调控细胞凋亡相关基因的表达介导细胞凋亡的发生。例如，Zhang等[13]在日本比目鱼(Paralichthys olivaceus)的研究中得出miR-731可以抑制p53的表达从而抑制细胞凋亡；Yi等[14]在日本囊对虾(Penaeus japonicus)的研究中发现过表达miR-1000可抑制p53的表达从而抑制细胞凋亡发生，而沉默miR-1000则显著提高p53的表达从而促进细胞凋亡的发生。
斑节对虾(P. monodon)等无脊椎动物因缺乏脊椎动物特有的获得性免疫系统，主要依赖先天性免疫系统抵御外界不良刺激。细胞凋亡属于对虾先天性免疫系统的血细胞防御系统中的重要一环，在斑节对虾抵御外界不良刺激过程中起重要作用[15-16]。基于TRIAP1参与调控细胞凋亡的发生，本研究克隆了斑节对虾TRIAP1基因(PmTRIAP1)，通过实时定量PCR (qRT-PCR)技术研究了其在不同组织中的表达分布规律以及哈氏弧菌(Vibrio harveyi)刺激下其在血细胞和肝胰腺中mRNA水平的变化情况；并对Pm-miR-145-3p、Pm-miR-454-3p和PmTRIAP1的关联关系进行了探究。旨在通过对PmTRIAP1及调控其表达的miRNAs表达调控的研究，为阐明PmTRIAP1在斑节对虾中抵御病菌入侵的分子机制奠定基础。
1. 材料与方法
1.1 实验材料
斑节对虾来自中国水产科学研究院南海水产研究所珠海斗门实验基地，体质量为(15±3) g，在(25±1) ℃、盐度3的海水中暂养3 d。每天更换养殖池内2/3的海水并按时投喂常规商品化对虾饲料，实验前24 h停止喂食。
哈氏弧菌由中国水产科学研究院南海水产研究所渔业生物病害防治研究室馈赠。双荧光素酶报告实验质粒(pMIR-REPORT)由和元生物技术有限公司(上海)构建。
人胚肾293T (HEK-293T)细胞株来源于中科院细胞库，细胞培养条件为DMEM+10% FBS、37 ℃、5% 二氧化碳 (CO2)。
1.2 实验方法
1.2.1 总RNA提取和cDNA合成
随机选取3尾健康斑节对虾迅速解剖，取心脏、淋巴、血细胞、肌肉、鳃、卵巢、脑、肠、胸神经、胃、肝胰脏11种组织速冻于液氮中。其中血细胞获取方法是：使用1 mL注射器抽取血液，转移至2 mL全新EP管中，立即4 ℃、4 000 g离心5 min，倒掉上清，取其沉淀。按照Trizol (Invitrogen，美国)说明书提取上述样品的总RNA，溶于DEPC水中。用NanoDrop-2000/2000c分光光度计检测浓度，通过1.2%琼脂糖凝胶电泳检测质量。使用PrimeScript TM逆转录试剂盒(TaKaRa)将RNA逆转录为cDNA。保存于−80 ℃。将cDNA样品工作液稀释至约40 ng·μL−1，并在−20 ℃下储存以备后用。
1.2.2 PmTRIAP1全长扩增
从实验室的斑节对虾转录组数据库中获得PmTRIAP1的部分cDNA序列。设计并合成特异性引物TRIAP1-3GSP1/3GSP2 (表1)，使用SMARTer RACE 5'/3'试剂盒 (TaKaRa)通过降落PCR和巢式PCR技术扩增目的基因3'末端。使用引物TRIAP1-F/R (表1)进行常规PCR验证PmTRIAP1基因cDNA的开放阅读框 (ORF) 的正确性。
表 1 实验所用引物Table 1. Sequences of primers used in this study引物名称
primer name序列 (5'−3')
sequence用途
applicationTRIAP1-F ATGAACAGTGTGAGCAAGGATT 开放阅读框验证 TRIAP1-R CTATTCCTTCGGCTTTTCTTCCG 开放阅读框验证 TRIAP1-3GSP1 GGAGGAAGAAAACGCTGGCAGTGG 3' RACE 扩增 TRIAP1-3GSP2 TGCAGCTGAGCTTGAGACCCAAAA 3' RACE 扩增 qTRIAP1-F TGAACAGTGTGAGCAAGGATT PmTRIAP1 RT-PCR qTRIAP1-R AGTGCTTCTTGGACACAGCTTTGG PmTRIAP1 RT-PCR EF-1α-F AAGCCAGGTATGGTTGTCAACTTT RT-PCR 参照基因 EF-1α-R CGTGGTGCATCTCCACAGACT RT-PCR 参照基因 qmiR-145-F GGATTCCTGGAAATACTG miR-145 RT-PCR qmiR-454-F TAGTGCAATATTGGCTA miR-454 RT-PCR qmiR-R CAGTGCGTGTCGTGGAGT miRNAs RT-PCR U6-S CTCGCTTCGGCAGCACA RT-PCR 内参基因 U6-A AACGCTTCACGAATTTGCGT RT-PCR 内参基因 1.2.3 PmTRIAP1的生物信息学分析
使用NCBI (http://blast.ncbi.nlm.nih.gov/Blast)的BLAST程序分析PmTRIAP1的全长cDNA序列。用ExPaSy计算工具 (http://web.expasy.org/protparam/)计算多肽链的分子量、氨基酸含量和等电点。用Clustal X 2.1进行多重序列比对(http://www.clustal.org/)。使用MEGA 7.0的邻接 (neighbor-joining, NJ)方法构建进化树。所有TRIAP1比对序列信息均从GenBank数据库获得。
1.2.4 PmTRIAP1组织表达分析
使用qRT-PCR检测1.2.1所述的斑节对虾11个组织中PmTRIAP1的表达。PCR反应总体系为12.5 μL，其中SYBR Premix Ex Taq (TaKaRa) 6.25 μL、cDNA模板1 μL、上游引物qTRAIP1-F和下游引物qTRIAP1-R各0.5 μL、双蒸水补足至12.5 μL。反应条件为94 ℃预变性30 s；94 ℃变性15 s，59 ℃退火延伸30 s，40个循环；溶解温度从65 ℃升至95 ℃；40 ℃冷却5 min。使用Roche Light-Cycler 480 Ⅱ实时定量PCR仪进行PCR，用2−ΔΔCt法计算基因相对表达量。结果用SPSS 22.0软件分析，表示为“平均值±标准差(
$\overline X \pm {\rm SD}$ )”，并进行单因素方差(One-Way ANOVA)分析，P<0.05为差异显著。PmTRIAP1基因特异性引物(qTRIAP1-F, qTRIAP1-R)、内参基因EF-1a引物(GenBank：DQ021452.1)见表1。1.2.5 预测与PmTRIAP1-3'UTR结合的miRNAs
结合实验室前期建立的斑节对虾小RNA文库(尚未发表)，用软件Mireap (http://www.miRNA.org/miRNA/)、MiRanda (http://www.miRNA/miRNA.org)和Targetscan (http://genes.mit.edu/targetscan)从中筛选出与PmTRIAP1-3'UTR结合的miRNAs。3个软件预测出来的共有miRNAs为最终结果。Mireap比较发夹前体结构和限制性位点的保守性、计算前体的折叠自由能来预测miRNAs；MiRanda评估2个序列(靶基因和miRNA)的互补性，通过计算自由能和进化保守度来过滤信息从而预测miRNAs；TargetScan根据靶基因与种子序列的最佳互补性评出最终得分预测miRNAs。Mireap和Targetscan中的参数是默认的，MiRanda中的参数修改为S≥140和E≤−14.0。
1.2.6 斑节对虾哈氏弧菌刺激实验
使用哈氏弧菌对斑节对虾进行刺激，并检测PmTRIAP1和预测的miRNAs的表达变化情况。在2216E培养基中培养哈氏弧菌，持续摇动(180 r·min−1) 12 h。用磷酸盐缓冲液(PBS，pH 7.4)制备得到1.0×108 CFU·mL−1的菌悬液。
本实验设置1个对照组和1个实验组，每组随机选取25尾健康的斑节对虾，分别饲养。稳定3 d后腹腔注射。对照组注射100 μL PBS溶液，实验组注射100 μL 1×108 CFU·mL−1的哈氏弧菌悬液。注射后分别于第0、第6、第24、第48、第72小时取样，5组重复。取其肝胰腺和血细胞储存于−80 ℃备用。
按照1.2.1的方法提取2种组织的总RNA并合成cDNA，PmTRIAP1荧光定量的表达分析、反应体系、反应程序、内参基因和数据处理方法同1.2.4。
miRNAs荧光定量实验反应总体系为25 μL，其中ddH2O 9 μL、SYBR Advantage Premix (2×) (TaKaRa) 12.5 μL、ROX Dye (50×) 0.5 μL、cDNA模板2 μL、上游引物qmiR-145或454和下游引物qmiR各0.5 μL。加好试剂后，将荧光定量PCR孔板放入Roche Light-Cycler 480 Ⅱ中，进行以下程序：94 ℃变性 10 s；94 ℃ 5 s，60 ℃ 20 s，40个循环；溶解曲线：94 ℃ 60 s，55 ℃ 30 s；95 ℃ 30 s。
所用引物qmiR-145、qmiR-454、qmiR、内参基因引物U6-S、U6-A见表1。数据处理方法同1.2.4。
1.2.7 双荧光素酶报告基因检测miRNAs
Pm-miR-454-3p和Pm-miR-145-3p类似物(agomir-454: 5'-TAGTGCAATATTGCTAATAGGCT-3'和agomir-145: 5'-GGATTCCTGGAAATACTGTTCT-3')及阴性对照agomir NC (agomir NC: 5'-TGATTTTTCAGCGCCTTGAAAG-3')委托上海和元生物技术有限公司合成。将Pm-miR-145-3p对应的靶基因PmTRIAP1 3' UTR (1 816~1 837 bp，5'-CAACAAAATCAACAGGAATCTG-3')和3' UTR突变序列 (5'-CAACAAAATCAACGATGCGATG-3') 以及Pm-miR-454-3p对应的靶基因PmTRIAP1 3'UTR (1 224~1 248 bp，5'-CAAGGGTATACCATTTATTGCACTT-3') 和3' UTR突变序列 (5'-CAAGGGTATACCATTTAGCTAGACT-3') 分别连接到pmir-report荧光质粒，从而构建含有PmTRIAP1 3'-UTR的嗜凤梨果蝇野生型质粒(miR-454-PmTRIAP1 3'UTR-WT和miR-145-PmTRIAP1 3'UTR-WT)和突变型质粒(miR-454-PmTRIAP1 3'UTR-MUT和miR-145-PmTRIAP1 3'UTR-MUT)。
质粒转染细胞。将293T细胞以70%的汇合度接种到96孔板中，24 h后转染萤光素酶报告基因质粒和RNA，每个样品设置6个复孔。96孔板转染质粒时每孔比例为萤火虫萤光素酶载体(Firefly)∶海肾萤光素酶载体 (Renilla)∶转染试剂=0.2 μg∶0.01 μg∶0.25 μL。配置miRNA和转染试剂稀释液，最终miRNA浓度为100 nmol·L−1，转染试剂每孔0.25 μL，常温孵育5 min，然后将稀释好的DNA和miRNA分别与转染试剂混匀，常温孵育20 min。每孔弃掉50 μL培养基，将25 μL DNA转染混合液和25 μL miRNA转染混合液分别添加到每孔细胞样品中。转染6 h后更换新鲜完全培养基。
双报告基因检测。质粒共转染48 h后，弃去培养基，用100 μL 1×100 PBS洗涤96孔板，吸干剩余的PBS。5×PLB (passive lysis buffer)用去离子水稀释成1×PLB，使用前置于常温。每孔加50 μL稀释好的1×PLB，摇床常温条件下摇15 min，进行裂解。白色不透光的96孔酶标板中每孔加入上述的10 μL上清液。在每个孔中加入100 μL预混合的LARⅡ(用Luciferase Assay Buffer溶解冻干粉Luciferase Assay substrate存于−20 ℃)，2 s后测数据。每孔加入100 μL预混合的Stop Glo & Reagent，静止2 s后测数据。
统计学分析。通过GraphPad Prism (Version 5，GraphPad Software)作图，并进行双侧t检验。
2. 结果
2.1 PmTRIAP1基因的分子特征和系统发育分析
通过SMART RACE技术，克隆获得斑节对虾PmTRIAP1基因cDNA全长(GenBank accession No. MK314973，图1)。PmTRIAP1 cDNA全长2 522 bp，含有11 bp的5'UTR、2 289 bp的3'UTR和222 bp的ORF，编码71个氨基酸，理论分子量为8.6 kD，理论等电点(PI)为4.9。
![]()
图 1 PmTRIAP1的核酸和氨基酸序列上方为核苷酸序列，下方为氨基酸序列；起始密码子(ATG)和终止密码子(TAG)用方框标出；倾斜加粗部分为加尾信号(AATAAA)Figure 1. Nucleotide and deduced amino acid sequences of PmTRIAP1The nucleotide sequence is listed above, and the amino acid sequence is listed below. The initiation code (ATG) and the termination code (TAG) are boxed. The polyadenylation signal sequence (AATAAA) is in italic and bold.多重序列比对结果显示，PmTRIAP1氨基酸序列与其他物种中的TRIAP1氨基酸序列相比有着较高的保守性。其中，PmTRIAP1与柑橘木虱(Diaphorina citri) TRIAP1的相似度最高(61%)，其他相似性较高的物种有墨西哥脂鲤(Astyanax mexicanus) 56%和家鸡(Gallus gallus) 51% (图2)。
![]()
图 2 PmTRIAP1和其他物种TRIAP1氨基酸序列多重比对Figure 2. Multiple alignment of deduced amino acid sequences of TRIAP1 from P. monodon and other speciesNCBI下载13个不同物种的TRIAP1序列与PmTRIAP1构建进化树，结果显示斑节对虾单独聚为一支 (图3)。
![]()
图 3 基于PmTRIAP1氨基酸序列的NJ系统进化树Figure 3. NJ phylogenetic tree based on PmTRIAP1 amino acid sequences2.2 PmTRIAP1的组织表达分析
qRT-PCR测定PmTRIAP1在11种组织中的表达情况(图4)。结果显示PmTRIAP1在11种组织中均有表达。血细胞中表达量最高，肝胰腺和肠中表达量较高，胸神经中表达量最低。
![]()
图 4 PmTRIAP1在各个组织中的相对表达量图中值为平均值±标准差(n=3)；不同字母表示差异显著(P<0.05)Figure 4. Relative expression of PmTRIAP1 in different tissuesValues are shown as $ \overline X $±SD (n=3). Different lowercase letters indicate significant difference (P<0.05).2.3 哈氏弧菌刺激下PmTRIAP1的表达分析
在血细胞中，哈氏弧菌刺激后，PmTRIAP1表达量在前6 h无显著变化，第24小时显著降低，第48小时降至最低水平，第72小时有所回升但仍显著低于对照组水平(图5-A)。在肝胰腺中，哈氏弧菌刺激后PmTRIAP1的表达水平一直低于磷酸缓冲盐组，第6小时表达无显著变化，第24小时表达显著降低，第48小时降到最低，第72小时明显回升，与第0小时相比无显著差异 (图5-B)。
![]()
图 5 哈氏弧菌注射后PmTRIAP1在血细胞 (A) 和肝胰腺 (B) 中的相对表达量图中值为平均值±标准差(n=3)；*. 显著性差异(P<0.05)；图6同此Figure 5. Relative expression levels of PmTRIAP1 in hemocyte (A) and hepatopancreas (B) after V. harveyi challengeValues are shown as $ \overline X $±SD (n=3); *. significant difference (P<0.05); the same case in Fig.6.2.4 初步筛选出靶向调节PmTRIAP1的miRNAs
Targetscan、miRanda和Mireap分别筛选靶向调控PmTRIAP1的miRNAs。经3个软件联合预测、筛选，获得2个潜在调控PmTRIAP1的miRNAs (表2)。
表 2 Targetscan、miRanda和Mireap筛选出2个miRNAsTable 2. Two miRNAs screened by Targetscan，miRanda and MireapmiRNA名称 miRNA name 序列 (5'−3') sequence Pm-miR-145-3p GGATTCCTGGAAATACTGTTCT Pm-miR-454-3p TAGTGCAATATTGCTAATAGGCT 通过1.2.6的哈氏弧菌刺激实验检测了血细胞中感染哈氏弧菌后2种miRNAs的表达模式。qRT-PCR结果显示Pm-miR-454-3p和Pm-miR-145-3p的表达水平在第6、第24、第48、第72小时均显著上调 (图6)。
![]()
图 6 哈氏弧菌注射后Pm-miR-145-3p (A) 和Pm-miR-454-3p (B) 在血细胞中的相对表达量Figure 6. Relative expression level of Pm-miR-145-3p (A) and Pm-miR-454-3p (B) in hemocyte after V. harveyi challenge2.5 双荧光素酶报告基因检测
为进一步验证Pm-miR-454-3p和Pm-miR-145-3p对PmTRIAP1的靶向调控作用，分别根据预测的Pm-miR-454-3p (图7-A)和Pm-miR-145-3p (图7-B)的调控位点开展双荧光素报告实验。
![]()
图 7 Pm-miR-454-3p和Pm-miR-145-3p与靶基因PmTRIAP1 3'-UTR结合位点序列信息A. PmTRIAP1 3'-UTR和Pm-miR-454-3p；B. PmTRIAP1 3'-UTR和Pm-miR-145-3pFigure 7. Sequence information of binding and mutant sites of each miRNA at 3'-UTR of target genesA. 3'-UTR of PmTRIAP1 and Pm-miR-454-3p; B. 3'-UTR of PmTRIAP1 and Pm-miR-145-3p双荧光素酶报告实验结果表明，构建的野生型报告质粒miR-145-PmTRIAP1 3'UTR-WT的荧光素酶活性受到Pm-miR-145-3p类似物(agomir-145)的抑制作用，活性大约为对照组miR-145-PmTRIAP1 3'UTR-MUT+agomir NC的76% (P=0.004 7)，实验组miR-145-PmTRIAP1 3'UTR-WT+agomir NC和miR-145-PmTRIAP1 3'UTR-MUT+agomir-145均未检测到明显活性变化(图8-A)。构建的野生型报告质粒miR-454-PmTRIAP1 3'UTR-WT的荧光素酶活性受到Pm-miR-454-3p类似物(agomir-454)的抑制作用，活性大约为对照组miR-454-PmTRIAP1 3'UTR-MUT+agomir NC的43% (P=0.042 6)，实验组miR-454-PmTRIAP1 3'UTR-WT+agomir NC和miR-454-PmTRIAP1 3'UTR-MUT+agomir-454均未检测到明显的活性变化(图8-B)。
![]()
图 8 miRNAs与靶基因PmTRIAP1间互作关系的双荧光素报告验证实验A. PmTRIAP1和Pm-miR-145-3p之间互作关系；B. PmTRIAP1和Pm-miR-454-3p之间互作关系Figure 8. Interaction between miRNAs and PmTRIAP1 validated by dual-luciferase reporter assaysA. interaction of PmTRIAP1 and Pm-miR-145-3p; B. interaction of PmTRIAP1 and Pm-miR-454-3p3. 讨论
斑节对虾是华南和东南亚最具商业价值的水产养殖物种之一。目前，弧菌等病害频发给对虾养殖产业造成了巨大的经济损失，与抵御弧菌入侵相关的分子调控机制也尚未清晰。本研究拟通过对斑节对虾TRIAP1基因与弧菌刺激之间的关联性研究，为解析斑节对虾抵御弧菌入侵的分子机制提供一定理论基础。
本研究成功克隆获得斑节对虾PmTRIAP1 cDNA全长序列，并对PmTRIAP1在斑节对虾各组织中的表达分布规律进行了探究。由于PmTRIAP1在与斑节对虾生殖调控密切相关的卵巢组织、与免疫相关的血淋巴组织、与生长相关的肌肉组织等各待测组织中均有广泛表达，预示该基因在斑节对虾体内可能参与调控多种不同的生理过程。PmTRIAP1的表达量在血细胞中最高，肝胰腺中次之。已有研究表明血细胞在免疫中扮演重要角色，肝胰腺是对虾非特异性免疫组织，也是重要的解毒器官，在应激反应中发挥重要[17-18]，进而初步说明PmTRIAP1可能参与了斑节对虾的免疫调控过程[19]。有鉴于此，本实验选取了肝胰腺和血细胞作为研究对象，对哈氏弧菌刺激下PmTRIAP1以及可能调控PmTRIAP1的miRNAs的表达变化情况进行了探究。
哈氏弧菌是水产养殖中十分常见的条件致病菌，在海洋环境中广泛存在，能够引发水生动物大规模的疾病爆发[20]。研究表明哈氏弧菌感染能引发细胞凋亡，例如，Li等[21]发现哈氏弧菌可导致日本牙鲆(Paralichthys olivaceus)复发性疾病爆发并诱导宿主细胞凋亡。TRIAP1作为一种细胞凋亡抑制因子，TRIAP1基因的表达水平对于细胞凋亡的调控至关重要。例如，Potting等[3]发现TRIAP1蛋白阻止细胞凋亡的发生；Fook-Alves等[2]发现在鼻咽喉癌细胞中TRIAP1基因的沉默导致了细胞凋亡增加。本研究发现，在哈氏弧菌刺激下，肝胰腺中PmTRIAP1的表达量在前6 h内无明显变化，第24小时表达水平显著降低，而在第72小时PmTRIAP1表达量与对照组却没有明显差异(图5-B)；在血细胞中，PmTRIAP1表达规律与在肝胰腺中的基本相似，不同的是，在第72小时PmTRIAP1的表达量仍然显著低于对照组(图5-A)。本研究结果进一步证明了PmTRIAP1基因的表达水平与哈氏弧菌感染之间存在显著的关联性，并推断哈氏弧菌感染后，斑节对虾能通过降低自身PmTRIAP1的表达来促进细胞凋亡的发生，进而实现对哈氏弧菌的抵御功能。
生物体内，miRNA通过碱基互补配对的方式识别靶mRNA，并根据互补程度的不同，指导沉默复合体降解靶mRNA或者阻遏靶mRNA的翻译。例如，Zhang等[22]研究发现在灿烂弧菌(V. splendidus)感染仿刺参时，miR-137呈现下调趋势，而其靶基因AjBHMT呈现上调趋势；Huang和Zhang[23]研究发现，在对虾白斑综合症病毒 (white spot syndrome virus, WSSV)侵染对虾时，机体内miR-7显著上调，而WSSV病毒的早期基因wsv-477显著下调；Li等[24]发现在鼻咽喉癌细胞中，miR-320b的下调可介导TRIAP1过表达。本实验利用生物信息学软件分析筛选出了可能靶向调控PmTRIAP1的2个miRNAs，并研究了哈氏弧菌刺激后2个miRNAs的表达变化情况。结果表明，在哈氏弧菌刺激24 h后，Pm-miR-145-3p和Pm-miR-454-3p表达量均显著升高(图6)，而PmTRIAP1表达量显著降低，表明PmTRIAP1可能是Pm-miR-145-3p和Pm-miR-454-3p的靶基因(图5)。为了进一步证明Pm-miR-145-3p与Pm-miR-454-3p可以靶向调控PmTRIAP1的表达，对其进行了双荧光素酶报告检测实验。结果显示，Pm-miR-145-3p显著降低了含有PmTRIAP1 3'UTR的荧光素酶活性，表明Pm-miR-145-3p可以靶向调控PmTRIAP1。
综上，本研究克隆了PmTRIAP1基因cDNA全长，探究了在哈氏弧菌刺激下PmTRIAP1 mRNA和相关miRNAs的表达变化情况。在哈氏弧菌刺激下，PmTRIAP1 mRNA表达量显著降低，Pm-miR-145-3p和Pm-miR-454-3p表达量显著升高，说明了PmTRIAP1基因、Pm-miR-145-3p和Pm-miR-454-3p都参与了斑节对虾抵御哈氏弧菌的应激反应。研究发现，PmTRIAP1和2个miRNAs的表达变化规律呈负相关，结合双荧光素酶报告实验，进一步推测出Pm-miR-145-3p和Pm-miR-454-3p靶向调控PmTRIAP1，进而对哈氏弧菌刺激做出应答。本实验为深入探究斑节对虾在抵御细菌感染时，先天免疫调节机制的激活提供了参考。
-
![]()
图 1 PmTRIAP1的核酸和氨基酸序列
上方为核苷酸序列，下方为氨基酸序列；起始密码子(ATG)和终止密码子(TAG)用方框标出；倾斜加粗部分为加尾信号(AATAAA)
Figure 1. Nucleotide and deduced amino acid sequences of PmTRIAP1
The nucleotide sequence is listed above, and the amino acid sequence is listed below. The initiation code (ATG) and the termination code (TAG) are boxed. The polyadenylation signal sequence (AATAAA) is in italic and bold.
![]()
图 2 PmTRIAP1和其他物种TRIAP1氨基酸序列多重比对
Figure 2. Multiple alignment of deduced amino acid sequences of TRIAP1 from P. monodon and other species
![]()
图 3 基于PmTRIAP1氨基酸序列的NJ系统进化树
Figure 3. NJ phylogenetic tree based on PmTRIAP1 amino acid sequences
![]()
图 4 PmTRIAP1在各个组织中的相对表达量
图中值为平均值±标准差(n=3)；不同字母表示差异显著(P<0.05)
Figure 4. Relative expression of PmTRIAP1 in different tissues
Values are shown as $ \overline X $±SD (n=3). Different lowercase letters indicate significant difference (P<0.05).
![]()
图 5 哈氏弧菌注射后PmTRIAP1在血细胞 (A) 和肝胰腺 (B) 中的相对表达量
图中值为平均值±标准差(n=3)；*. 显著性差异(P<0.05)；图6同此
Figure 5. Relative expression levels of PmTRIAP1 in hemocyte (A) and hepatopancreas (B) after V. harveyi challenge
Values are shown as $ \overline X $±SD (n=3); *. significant difference (P<0.05); the same case in Fig.6.
![]()
图 6 哈氏弧菌注射后Pm-miR-145-3p (A) 和Pm-miR-454-3p (B) 在血细胞中的相对表达量
Figure 6. Relative expression level of Pm-miR-145-3p (A) and Pm-miR-454-3p (B) in hemocyte after V. harveyi challenge
![]()
图 7 Pm-miR-454-3p和Pm-miR-145-3p与靶基因PmTRIAP1 3'-UTR结合位点序列信息
A. PmTRIAP1 3'-UTR和Pm-miR-454-3p；B. PmTRIAP1 3'-UTR和Pm-miR-145-3p
Figure 7. Sequence information of binding and mutant sites of each miRNA at 3'-UTR of target genes
A. 3'-UTR of PmTRIAP1 and Pm-miR-454-3p; B. 3'-UTR of PmTRIAP1 and Pm-miR-145-3p
![]()
图 8 miRNAs与靶基因PmTRIAP1间互作关系的双荧光素报告验证实验
A. PmTRIAP1和Pm-miR-145-3p之间互作关系；B. PmTRIAP1和Pm-miR-454-3p之间互作关系
Figure 8. Interaction between miRNAs and PmTRIAP1 validated by dual-luciferase reporter assays
A. interaction of PmTRIAP1 and Pm-miR-145-3p; B. interaction of PmTRIAP1 and Pm-miR-454-3p
表 1 实验所用引物
Table 1 Sequences of primers used in this study
引物名称
primer name序列 (5'−3')
sequence用途
applicationTRIAP1-F ATGAACAGTGTGAGCAAGGATT 开放阅读框验证 TRIAP1-R CTATTCCTTCGGCTTTTCTTCCG 开放阅读框验证 TRIAP1-3GSP1 GGAGGAAGAAAACGCTGGCAGTGG 3' RACE 扩增 TRIAP1-3GSP2 TGCAGCTGAGCTTGAGACCCAAAA 3' RACE 扩增 qTRIAP1-F TGAACAGTGTGAGCAAGGATT PmTRIAP1 RT-PCR qTRIAP1-R AGTGCTTCTTGGACACAGCTTTGG PmTRIAP1 RT-PCR EF-1α-F AAGCCAGGTATGGTTGTCAACTTT RT-PCR 参照基因 EF-1α-R CGTGGTGCATCTCCACAGACT RT-PCR 参照基因 qmiR-145-F GGATTCCTGGAAATACTG miR-145 RT-PCR qmiR-454-F TAGTGCAATATTGGCTA miR-454 RT-PCR qmiR-R CAGTGCGTGTCGTGGAGT miRNAs RT-PCR U6-S CTCGCTTCGGCAGCACA RT-PCR 内参基因 U6-A AACGCTTCACGAATTTGCGT RT-PCR 内参基因 
下载: 导出CSV
表 2 Targetscan、miRanda和Mireap筛选出2个miRNAs
Table 2 Two miRNAs screened by Targetscan，miRanda and Mireap
miRNA名称 miRNA name 序列 (5'−3') sequence Pm-miR-145-3p GGATTCCTGGAAATACTGTTCT Pm-miR-454-3p TAGTGCAATATTGCTAATAGGCT
下载: 导出CSV
-
[1] PARK W R, NAKAMURA Y. p53CSV, a novel p53-inducible gene involved in the p53-dependent cell-survival pathway[J]. Cancer Res, 2005, 65(4): 1197-1206. doi: 10.1158/0008-5472.CAN-04-3339
[2] FOOK-ALVES V L, de OLIVEIRA M B, ZANATTA D B, et al. TP53 regulated inhibitor of apoptosis 1 (TRIAP1) stable silencing increases late apoptosis by upregulation of caspase9 and APAF1 in RPMI8226 multiple myeloma cell line[J]. Biochim Biophys Acta, 2016, 1862(6): 1105-1110. doi: 10.1016/j.bbadis.2016.03.011
[3] POTTING C, TATSUTA T, KÖNIG T, et al. TRIAP1/PRELI complexes prevent apoptosis by mediating intramitochondrial transport of phosphatidic acid[J]. Cell Metab, 2013, 18(2): 287-295. doi: 10.1016/j.cmet.2013.07.008
[4] OLIVIERI F, RIPPO M R, MONSURRÒ V, et al. MicroRNAs linking inflamm-aging, cellular senescence and cancer[J]. Ageing Res Rev, 2013, 12(4): 1056-1068. doi: 10.1016/j.arr.2013.05.001
[5] BEITZINGER M, PETERS L, ZHU J Y, et al. Identification of human microRNA targets from isolated argonaute protein complexes[J]. RNA Biol, 2014, 4(2): 76-84.
[6] OU J, LI Y, DING Z, et al. Transcriptome-wide identification and characterization of the Procambarus clarkii microRNAs potentially related to immunity against Spiroplasma eriocheiris infection[J]. Fish Shellfish Immunol, 2013, 35(2): 607-617. doi: 10.1016/j.fsi.2013.05.013
[7] KUHN D E, MARTIN M M, FELDMAN D S, et al. Experimental validation of miRNA targets[J]. Methods, 2008, 44(1): 47-54. doi: 10.1016/j.ymeth.2007.09.005
[8] TAN T T, CHEN M, HARIKRISHNA J A, et al. Deep parallel sequencing reveals conserved and novel miRNAs in gill and hepatopancreas of giant freshwater prawn[J]. Fish Shellfish Immunol, 2013, 35(4): 1061-1069. doi: 10.1016/j.fsi.2013.06.017
[9] HE L, WANG Y L, LI Q, et al. Profiling microRNAs in the testis during sexual maturation stages in Eriocheir sinensis[J]. Anim Reprod Sci, 2015, 162: 52-61. doi: 10.1016/j.anireprosci.2015.09.008
[10] LI S, ZHU S, LI C, et al. Characterization of microRNAs in mud crab Scylla paramamosain under Vibrio parahaemolyticus infection[J]. PLoS One, 2013, 8(8): e73392.
[11] CHEN Y, FU L, WEN X, et al. Oncogenic and tumor suppressive roles of microRNAs in apoptosis and autophagy[J]. Apoptosis, 2014, 19(8): 1177-1189. doi: 10.1007/s10495-014-0999-7
[12] AN F, GONG B, WANG H, et al. miR-15b and miR-16 regulate TNF mediated hepatocyte apoptosis via BCL2 in acute liver failure[J]. Apoptosis, 2012, 17(7): 702-716. doi: 10.1007/s10495-012-0704-7
[13] ZHANG B C, ZHOU Z J, SUN L. pol-miR-731, a teleost miRNA upregulated by megalocytivirus, negatively regulates virus-induced type I interferon response, apoptosis, and cell cycle arrest[J]. Sci Rep, 2016, 6: 28354. doi: 10.1038/srep28354
[14] GONG Y, JU C, ZHANG X. The miR-1000-p53 pathway regulates apoptosis and virus infection in shrimp[J]. Fish Shellfish Immunol, 2015, 46(2): 516-522. doi: 10.1016/j.fsi.2015.07.022
[15] IWANAGA S, LEE B L. Recent advances in the innate immunity of invertebrate animals[J]. BMB Rep, 2005, 38(2): 128-150. doi: 10.5483/BMBRep.2005.38.2.128
[16] LITTLE T J, HULTMARK D, READ A F. Invertebrate immunity and the limits of mechanistic immunology[J]. Nat Immunol, 2005, 6(7): 651. doi: 10.1038/ni1219
[17] 戴文婷, 傅明骏, 赵超, 等. 斑节对虾CDK2基因全长cDNA克隆及表达分析[J]. 南方水产科学, 2015, 11(2): 1-11. doi: 10.3969/j.issn.2095-0780.2015.02.001 [18] 史进选, 傅明骏, 赵超, 等. 斑节对虾GRP94基因的克隆及其在不同应激条件下的表达与分析[J]. 南方水产科学, 2016, 12(5): 61-70. doi: 10.3969/j.issn.2095-0780.2016.05.008 [19] 傅明骏, 赵超, 杨其彬, 等. 斑节对虾过氧化氢酶基因全长cDNA克隆及表达分析[J]. 南方水产科学, 2015, 11(6): 107-113. doi: 10.3969/j.issn.2095-0780.2015.06.015 [20] 王瑞旋, 耿玉静, 冯娟, 等. 杂色鲍哈维弧菌耐药质粒的鉴定和分析[J]. 南方水产科学, 2012, 8(2): 1-6. doi: 10.3969/j.issn.2095-0780.2012.02.001 [21] LI M F, WANG C L, SUN L. A pathogenic Vibrio harveyi lineage causes recurrent disease outbreaks in cultured Japanese flounder (Paralichthys olivaceus) and induces apoptosis in host cells[J]. Aquaculture, 2011, 319(1/2): 30-36.
[22] ZHANG P, LI C, ZHANG R, et al. The roles of two miRNAs in regulating the immune response of sea cucumber[J]. Genetics, 2015, 201(4): 1397-1410. doi: 10.1534/genetics.115.178871
[23] HUANG T, ZHANG X. Functional analysis of a crustacean microRNA in host-virus interaction[J]. J Virol, 2012, 86(23): 12997-13004. doi: 10.1128/JVI.01702-12
[24] LI Y, TANG X, HE Q, et al. Overexpression of mitochondria mediator gene TRIAP1 by miR-320b loss is associated with progression in nasopharyngeal carcinoma[J]. PLoS Genet, 2016, 12(7): e1006183. doi: 10.1371/journal.pgen.1006183
计量
- 文章访问数: 4662
- HTML全文浏览量: 1978
- PDF下载量: 42
粤公网安备 44010502001741号
