miR-148a-3p靶向PPARγ基因抑制鹅颗粒细胞孕酮的合成

2021-06-30解广娟胡深强王继文

畜牧兽医学报 2021年6期

邓艳，解广娟，胡深强，李亮，王继文

(四川农业大学畜禽遗传资源发掘与创新利用四川省重点实验室，成都 613111)

microRNAs(miRNAs)是一类由约22个核苷酸组成的短链非编码RNA分子，其通过自身种子序列与靶基因特定位点进行碱基配对结合[1-2]，广泛参与动物卵巢类固醇激素合成与分泌的调控过程[3-5]。

近年来，大量研究已报道，miRNAs能调节卵巢颗粒细胞类固醇激素的生物合成[6-8]。研究发现，miR-15a参与了人卵巢颗粒细胞中黄体酮、雄激素和雌激素的释放，并且在过表达miR-15a后可以促进孕酮和睾酮生成[9]。Zhang等[10]发现，在FSH处理后的牛颗粒细胞中，miR-31和miR-143通过靶向FSHR基因抑制孕酮的分泌。在猪颗粒细胞中，miR-378通过靶向芳香化酶的3′UTR来调节卵巢中雌二醇的分泌[4]。而在小鼠颗粒细胞中，miR-383可以靶向RBMS1促进雌二醇的代谢[5]。这些结果表明，miRNAs在哺乳动物卵巢颗粒细胞的类固醇激素合成与分泌方面有重要调控作用。相比miRNAs在哺乳动物中的研究，其在禽类卵巢中的研究大多集中在对卵巢细胞增殖和凋亡的调控。比如，Kang等[11]研究表明，miR-26a-5p通过靶向TNRC6A基因促进鸡卵巢卵泡膜细胞的增殖。miR-26b-5p受到转录因子SMAD4的调控，靶向INHBB介导TGF-β/SMADs信号通路的表达，从而促进鹅卵泡颗粒细胞的增殖[12]。此外，鹅颗粒细胞的增殖和凋亡也受miR-181a-5p的调控，其通过靶向SIRT1参与FOXO信号通路[13]。然而，miRNAs 对禽类卵泡颗粒细胞的类固醇激素的合成与分泌的调控机制还未见报道。

过氧化物酶体增殖物激活型受体(peroxisome proliferator activated receptors，PPARs)是核受体超家族的配体依赖型转录因子，主要包括α、β、γ三种亚型[14-15]。其中，PPARγ参与调节卵巢的卵泡发育、排卵及卵子的质量和类固醇激素(雌激素和孕激素)的分泌[16-17]。Kim等[16]研究表明，PPARγ通过调控孕激素受体参与小鼠颗粒细胞类固醇激素的合成。然而，PPARγ基因是否影响鹅卵巢颗粒细胞的生物学功能目前尚不清楚。

本课题组前期转录组测序发现，miR-148a-3p在鹅4～6 mm、8～10 mm和F5卵泡颗粒层中的表达逐渐升高。因此，本研究首先检测miR-148a-3p在鹅不同阶段卵泡颗粒层中的表达，然后在体外培养的鹅卵泡颗粒细胞中通过转染mimics和inhibitor模拟miR-148a-3p的过表达和抑制，利用qPCR法检测与类固醇激素合成分泌相关基因的表达变化，接着通过双荧光素酶报告系统验证靶标关系，最后干扰靶基因检测miR-148a-3p对孕酮合成变化的影响，为进一步阐明miR-148a-3p在禽类颗粒细胞中的功能提供依据。

1 材料与方法

1.1 试验动物和样品采集

本试验共选用12只天府肉鹅母系母鹅，其中，6只进行组织样品采集，6只进行颗粒细胞的分离培养。所选用的试验鹅均由四川农业大学家禽育种场提供，试验鹅在相同的条件下饲养，开产时间和体重基本一致。鹅麻醉致死后，迅速取出整个卵巢，用游标卡尺测量不同直径大小的卵泡(2～4 mm、4～6 mm、 6～8 mm、8～10 mm、F5、F4、F3、F2、F1)，根据甘翔等[18]的方法进一步分离不同阶段卵泡颗粒层，PBS漂洗3～4次，将颗粒层附着的卵黄清洗干净后，剪碎并置于-80 ℃冻存，用于RNA提取。

1.2 miR-148a-3p在不同物种中的保守性分析

在miRBase(http://www．Mirbase.org/)数据库中在线检索人(Homosapiens, GenBank_No.：NC_000003.12)、小鼠(Musmusculus, GenBank_No.：NC_000072.6)、鸡(Gallusgallus, GenBank_No.：NC_006099.5)、猪(Susscrofa, GenBank_No.：NC_010455.5)、牛(BosTaurus, GenBank_No.：NC_037349.1)、恒河猴(Macacamulatta, GenBank_No.：NC_041755.1)、绿蜥蜴(Anoliscarolinensis, GenBank_No.：NC_014777.1)、大西洋鲑鱼(Salmosalar, GenBank_No.：NC_027314.1)8个物种的miR-148a-3p成熟体序列，鹅的miR-148a-3p成熟体序列由课题组前期miRNA测序所得，采用MEGA 7进行保守性分析。

1.3 鹅等级卵泡颗粒细胞分离培养与转染

参考Deng等[19]的方法进行鹅等级卵泡颗粒细胞的体外培养，待细胞汇合度为70%～80%时，根据Lipofectamine®3000(Invitrogen，USA)试剂使用说明，先将lipofectamine 3000分别和miR-148a-3p mimics、mimics NC、inhibitor、inhibitor NC、Control在DMEM/F12(HyClone，USA)培养基内形成脂质体复合物，然后共转染到颗粒细胞中，在含5% CO2的37 ℃培养箱中培养。转染后6 h对细胞状态进行一次观察，转染24 h后收取上清液和细胞用于后续的检测。miR-148a-3p mimics、mimics NC、inhibitor、inhibitor NC由上海吉玛制药技术有限公司合成(表1)。

表1 miR-148a-3p mimics和inhibitor合成序列

1.4 鹅卵泡颗粒层组织和颗粒细胞RNA的提取与反转录

利用RNAiso Plus(TaKaRa, China)提取鹅卵泡颗粒层组织和细胞的总RNA，利用Nano Drop 2000紫外分光光度计检测总RNA的浓度和纯度，使用Prime ScripTMRT reagent Kit(TaKaRa, China)进行mRNA反转录，按照上海吉玛制药技术有限公司的miRNA反转录试剂盒反转录miRNA，反应结束后置于-20 ℃保存备用。

1.5 qPCR检测miR-148a-3p及关键基因的表达

利用课题组前期对鹅颗粒细胞miRNA测序得到的miR-148a-3p序列，由上海吉玛制药技术有限公司合成miR-148a-3p定量检测试剂盒，其他mRNAs 定量检测引物由华大基因科技有限公司合成(表2)。使用SYBR®Premix Ex TaqTM II(TaKaRa, China)检测mRNAs表达量，mRNAs反应体系为：SYBR®Premix Ex TaqTM II 6.25 μL，上、下游引物各0.5 μL，模板1 μL，无菌水4.25 μL，总反应体系为12.5 μL。反应条件：95 ℃ 10 min；95 ℃ 30 s，退火温度30 s，72 ℃延伸30 s， 72 ℃ 10 min, 总共40个循环；设置熔解曲线为65～95 ℃，每5 s增加0.5 ℃，以检测引物特异性。使用miR-148a-3p定量检测试剂盒检测miR-148a-3p的表达量，miR-148a-3p的反应体系为：2×Real-time PCR Master Mix (SYBR) 10 μL，上、下游引物各0.4 μL，TaqDNA polymerase 0.2 μL，模板2 μL，无菌水7 μL，总反应体系为20 μL。反应条件为：95 ℃ 3 min；95 ℃ 12 s，退火温度40 s，共40个循环。每个样本设3个技术重复，miRNA表达量计算以U6为内参，mRNAs表达量计算以GAPDH和β-actin为内参。

表2 用于qPCR的引物

1.6 突变型载体和野生型载体构建

通过miRecords_version4、miRTarBase、Tarbase 7.0等数据库中收录的经试验验证的miR-148a-3p靶基因，综合TragetScan 7.1、miRanda、miRWalk 2.0、DIANA-microT等在线软件预测的结果发现，PPARγ-3′UTR与miR-148a-3p种子序列存在结合位点(图1)。根据鹅PPARγ-3′UTR 序列设计含结合位点在内的扩增引物(表3)，将所得到的PPARγ-3′UTR目的片段和pmir GLO空载质粒(Promega，USA)分别用SacⅠ和XbaⅠ进行双酶切，酶切反应完毕后进行纯化回收。将经过切胶纯化后的目的片段和已切开的pmir GLO(Promega，USA)空载质粒进行连接，构建野生型靶标验证载体PGL/PPARγ-WT。将野生型载体的潜在结合位点突变为其互补序列(图1)，从而得到突变型重组载体即PGL/PPARγ-Mut，突变载体由上海吉玛制药技术有限公司完成。

黑色加粗字体为种子序列和突变序列

表3 PPARγ-3′UTR序列扩增引物

1.7 双荧光素酶报告试验验证miR-148a-3p的靶基因

使用中国仓鼠卵巢细胞系(CHO，中国科学院昆明动物研究所)用于双荧光素酶报告系统检测。CHO细胞系复苏后接种于培养皿中，待细胞长满后用胰酶消化，然后用含10% FBS(Gibco, USA)的DMEM/F12(HyClone, USA)培养基重悬细胞，以4×106个·mL-1的密度接种于96孔板。待细胞密度达到约60%～70%时，按照Lipofectamine®3000(Invitrogen，USA)说明书将miR-148a-3p mimics分别与PGL/PPARγ-WT或PGL/PPARγ-Mut两两进行共转染，每个处理设置3孔重复。转染后24 h使用Dual-Luciferase®Reporter试剂盒(Promega, China)检测细胞荧光素酶活性。以萤火虫荧光素酶作为试验报告基因，海肾荧光素酶作为对照报告基因，二者的比值即为相对荧光素活性。

1.8 miR-148a-3p靶基因的干扰

从NCBI获取鹅PPARγ CDS区预测序列，由上海吉玛制药技术有限公司构建si-PPARγ，序列见表4。结合Lipofectamine®3000(Invitrogen，USA)试剂使用说明，使用12孔板转染siRNA及其NC，转染24 h后收集样品进行关键基因表达量的检测及激素含量的检测。

表4 PPARγ干扰序列

1.9 ELISA检测孕酮含量

细胞转染24 h后收取细胞上清液，根据ELISA试剂盒(上海恒远科技有限公司，中国)说明书在板条标准品孔内加50 μL不同浓度标准品，样本孔加10 μL待测样本和40 μL稀释液，空白孔不加，其余孔中加入100 μL辣根过氧化物酶(HRP)标记的抗体，封板膜密封，37 ℃恒温箱温育60 min；弃去液体，纸上拍干，每孔加满洗涤液，静置1 min，甩去洗涤液，纸上拍干，重复洗板5次；每孔各加入50 μL底物A、B，37 ℃避光孵育15 min；每孔加入50 μL终止液混匀，15 min内在450 nm波长处检测各孔OD值。在Excel中，以标准品浓度为横坐标，OD值为纵坐标，绘制出标准品线性回归曲线，按曲线方程计算各样本浓度值，从而得到不同处理组中孕酮含量。

1.10 数据的统计与分析

采用SPSS 19.0软件对试验数据进行单因素方差分析，组间差异显著性检验采用独立样本t检验进行分析，利用公式2-△△Ct计算所选目标基因的相对表达量[20]。P<0.05 表示具有统计学意义，数据用“平均值±标准误”表示。

2 结果

2.1 miR-148a-3p在不同物种中的保守性分析

通过miRBase数据库下载不同物种的miR-148a-3p的成熟体序列，比对不同物种的成熟体序列，结果发现，miR-148a-3p的成熟体序列在各个物种中高度保守，其核心的种子序列基本一致(表5)。

表5 miR-148a-3p序列保守性分析

2.2 miR-148a-3p在鹅不同阶段卵泡颗粒层中的表达分析

结果如图2所示，随着卵泡直径的增加，miR-148a-3p的表达量呈现先升高后下降的趋势，且在等级卵泡F4颗粒层中表达最高。此外，miR-148a-3p在等级阶段卵泡颗粒层的表达也高于等级前卵泡(图2)，因此，选择等级卵泡颗粒细胞用于后续的细胞试验。

不同小写字母表示差异显著(P<0.05)，下同

2.3 转染miR-148a-3p mimics和inhibitor对鹅卵泡颗粒细胞功能的影响

为了探究miR-148a-3p对鹅卵泡颗粒细胞功能的影响，在鹅等级卵泡颗粒细胞中转染miR-148a-3p mimics和inhibitor。定量结果显示，与mimics NC组相比，miR-148a-3p mimics能显著下调3β-HSD表达(图3A,P<0.05)，而与inhibitor NC组相比，miR-148a-3p inhibitor能显著上调3β-HSD表达(图3C,P<0.05)。同时，ELISA结果显示，与mimics NC组相比，miR-148a-3p mimics能显著抑制孕酮的合成(图3B,P<0.05)，而与inhibitor NC组相比，miR-148a-3p inhibitor能显著促进孕酮的合成(图3D,P<0.05)。由此可知，miR-148a-3p能抑制鹅等级卵泡颗粒细胞中孕酮的合成。

A. miR-148a-3p mimics对3β-HSD的mRNA表达的影响；B. miR-148a-3p mimics对孕酮合成的影响；C. miR-148a-3p inhibitor对3β-HSD的mRNA表达的影响；D. miR-148a-3p inhibitor对孕酮合成的影响。NC. 阴性对照组；mimics. miR-148a-3p mimics；inhibitor. miR-148a-3p inhibitor；Control. 空白对照组

2.4 miR-148a-3p靶基因的确定

为了进一步探究miR-148a-3p抑制鹅卵泡颗粒细胞孕酮合成的分子机制，在CHO细胞系中共转染miR-148a-3p mimics与野生型载体PGL/PPARγ-WT或突变型载体PGL/PPARγ-Mut，采用双荧光素酶报告系统检测双荧光素酶的活性。结果显示，当miR-148a-3p与PPARγ的结合位点突变后，双荧光素酶活性显著升高(图4，P<0.05)，证实PPARγ是miR-148a-3p的靶基因。

图4 miR-148a-3p与PPARγ靶标关系的双荧光素检测

2.5 miR-148a-3p mimics和inhibitor对PPARγ表达的影响

为了探究miR-148a-3p对靶基因的转录调控的影响，在鹅卵泡颗粒细胞中转染miR-148a-3p mimics和miR-148-3p inhibitor后检测靶基因PPARγ的mRNA表达变化。结果如图5所示，与mimics NC组相比，miR-148a-3p mimics能显著降低鹅颗粒细胞中PPARγ的表达(图5A,P<0.05)。与inhibitor NC组相比，miR-148a-3p inhibitor能显著上调PPARγ的表达(图5B,P<0.05)。由此可知，miR-148a-3p能调控靶基因PPARγ的转录水平。

A. miR-148a-3p mimics对PPARγ的mRNA表达的影响；B. miR-148a-3p inhibitor对PPARγ的mRNA表达的影响

2.6 干扰PPARγ对鹅卵泡颗粒细胞孕酮分泌的影响

为了探究miR-148a-3p是否会通过其靶基因的转录调控影响鹅卵泡颗粒细胞的激素合成，在颗粒细胞中干扰PPARγ后检测关键基因的表达变化及孕酮含量。结果如图6所示，与对照组和si-PPARγ NC组相比，干扰PPARγ基因后，3β-HSD的表达量显著降低(图6A,P<0.05)，孕酮的含量也降低，但差异不显著(图6B,P>0.05)。由此可知，鹅颗粒细胞中孕酮的含量受PPARγ的转录调控。

A.干扰PPARγ对3β-HSD的mRNA表达的影响；B.干扰PPARγ对孕酮合成的影响

3 讨论

miR-148a-3p的序列已经在多个物种中被鉴定[21]，而在鹅上，miR-148a-3p的序列是通过前期转录组测序获得。本研究通过对多个物种的序列进行保守性比对分析发现，miR-148a-3p在不同物种中高度保守，推测miR-148a-3p在不同的物种中可能具有类似的功能，但关于miR-148a-3p调控鹅卵泡颗粒细胞类固醇激素分泌的相关研究未见报道。

目前的研究发现，miR-148a-3p在许多生理过程中均发挥重要作用，比如免疫、肿瘤和癌症[22-23]，骨骼肌细胞的分化、增殖和凋亡以及脂肪细胞的增殖分化[24-26]。本研究中，miR-148a-3p在鹅不同发育阶段的卵泡颗粒层中均有表达，且随着卵泡直径的增加，其表达量呈现先升高后下降的趋势，在F4、F5阶段卵泡达到顶峰。F4、F5阶段卵泡为等级前卵泡发育到等级卵泡的过渡时期[27]，推测miR-148a-3p可能在卵泡发育或卵泡选择过程中有重要作用，这也表明，miR-148a-3p也参与禽类的繁殖生理过程。进一步体外培养鹅颗粒细胞，转染miR-148a-3p mimics能显著下调3β-HSD表达，抑制孕酮的合成；而转染miR-148a-3p inhibitor能显著上调3β-HSD表达，促进孕酮的合成。3β-HSD是孕酮合成的关键基因[28]，禽类卵泡颗粒层中3β-HSD的表达量随着卵泡的发育逐渐升高[29-30]。这些结果表明，miR-148a-3p可以抑制3β-HSD的表达降低孕酮的合成，进而影响鹅卵泡的选择和发育。miRNAs主要通过与靶基因的3′-UTR结合调控靶基因表达继而影响靶细胞的生物学功能[1]。研究表明，在食管癌中，miR-148a-3p通过靶向DNA甲基转移酶1(DNA (cytosine-5)-methyltransferase 1,DNMT1)调控癌细胞的增殖和入侵[31]。在人卵巢癌细胞中，miR-148a-3p通过靶向络氨酸蛋白激酶Met(tyrosine-protein kinase Met,c-Met)抑制癌细胞的增殖[32]。而在正常的细胞中，如C2C12肌细胞系和原代老鼠细胞中，miR-148a能直接通过靶向含蛋白激酶1的Rho相关螺旋线圈(Rho-associated coiled-coil containing protein kinase 1,ROCK1)的3′UTR来促进成肌细胞分化[33]。在牛的骨骼肌细胞中，miR-148a-3p通过转录后水平下调类Krüppel因子6(Krüppel-like factor 6,KLF6)的表达抑制成肌细胞的增殖，促进其凋亡[34]。而在鸡的骨骼肌卫星细胞中，miR-148a-3p通过直接靶向间质同源框2(mesenchyme homeobox 2,Meox2)来促进细胞的分化，抑制凋亡，但不影响增殖，同时能增加PI3K/AKT信号通路[35]。本课题组前期通过靶基因在线预测软件发现，miR-148a-3p与PPARγ3′-UTR存在潜在结合位点，双荧光素酶报告载体试验结果显示，miR-148a-3p可以直接靶向结合PPARγ调控鹅颗粒细胞的生理功能。进一步检测miR-148a-3p对PPARγ表达量的影响时发现，转染miR-148a-3p mimics能显著降低鹅颗粒细胞中PPARγ表达，而转染miR-148a-3p inhibitor后，PPARγ的表达显著升高。研究发现，PPARγ可以影响大鼠[36]、猪[37]、绵羊[38]、牛[39]卵巢颗粒细胞产生类固醇激素。Komar[40]研究表明，PPARγ可以引起人孕酮分泌量的改变。本研究中，干扰PPARγ能显著降低3β-HSD表达水平和孕酮的含量。这些结果就表明，在鹅卵泡颗粒细胞中，miR-148a-3p通过靶向PPARγ来抑制3β-HSD的表达，进而降低孕酮的合成。