畜牧兽医学报  2018, Vol. 49 Issue (7): 1366-1376. DOI: 10.11843/j.issn.0366-6964.2018.07.005    PDF    
牦牛Smad 4基因3'UTR区双荧光素酶载体构建及与bta-miR-146a的靶向验证
牛家强1,2, 王玉恒1,2, 索朗斯珠1,2, 强巴央宗1, 徐业芬1,2, 郭敏3, 程玲华3, 杨士承1,2     
1. 西藏农牧学院动物科学学院, 林芝 860000;
2. 西藏农牧学院西藏高原动物疫病研究 自治区高校重点实验室, 林芝 860000;
3. 中国农业大学动物科技学院, 北京 100193
摘要:旨在鉴定牦牛Smad 4基因与bta-miR-146a的靶向关系,为探究它们在牦牛卵巢卵泡发育过程中的可能分子调控机制提供理论基础。利用miRBase数据库和MEGA 6.0软件分析miR-146a在不同物种中的保守性,同时利用TargetScan软件预测其潜在靶基因,通过构建psiCHECK2双荧光素酶报告基因载体进行靶向验证。生物信息学分析结果显示,miR-146a成熟序列在脊椎动物中具有高度的保守性,并预测到178个与卵巢卵泡发育相关的潜在靶基因,从中选取Smad 4基因分析发现其3'非编码区域(3'UTR)与bta-miR-146a存在结合位点;进一步成功构建牦牛Smad 4基因3'UTR野生型及突变型psiCHECK2双荧光素酶报告质粒,荧光素酶活性检测结果显示,bta-miR-146a mimic对牦牛Smad 4基因3'UTR野生型报告质粒荧光活性有明显的下调作用,分别与突变型及NC对照组差异显著(P < 0.05)。本研究初步证实,牦牛Smad 4基因是bta-miR-146a的靶基因,从而提示bta-miR-146a可能通过调控Smad 4基因在牦牛卵巢卵泡发育过程中的表达发挥作用。
关键词牦牛    bta-miR-146a    Smad 4    荧光素酶活性    
Construction of Yak Smad 4 Gene 3'UTR Dual-luciferase Reporter Vector and Its Targeting Validation to bta-miR-146a
NIU Jia-qiang1,2, WANG Yu-heng1,2, SUOLANG Si-zhu1,2, CHANGBA Yang-zong1, XU Ye-fen1,2, GUO Min3, CHENG Ling-hua3, YANG Shi-cheng1,2     
1. Department of Animal Science, Tibet Agricultural and Animal Husbandry College, Linzhi 860000, China;
2. Provincial Key Laboratory of Tibet Plateau Animal Epidemic Disease Research, Tibet Agricultural and Animal Husbandry College, Linzhi 860000, China;
3. College of Animal Science and Technology, China Agricultural University, Beijing 100193, China
Abstract: The aim of this study was to identify the targeting relationship between Smad 4 gene and bta-miR-146a, which would provide a theoretical basis for their regulating mechanism in follicular development of yak ovary. The conservation of miR-146a in different species was analyzed by miRBase database and MEGA 6.0 software. At the same time, TargetScan software was used to predict the potential target genes of miR-146a, and the targeting validation was performed by construction of psiCHECK2 dual-luciferase reporter vector. The results of bioinformatics analysis showed that the miR-146a mature sequence was highly conserved in vertebrates, and a total of 178 potential target genes related to ovarian follicular development were predicted, then, among these genes, Smad 4 gene was selected to be further analyzed, showing that its 3'UTR had the binding site to bta-miR-146a. The wild and mutant type psiCHECK2 dual-luciferase reporter vectors of 3'UTR of yak Smad 4 gene were successfully constructed, and the dual-luciferase assay result showed that the luciferase activity of wild vector of 3'UTR of yak Smad 4 gene was distinctly down-regulated by the bta-miR-146a mimics, and was significantly lower than that of mutant vector and NC group(P < 0.05), respectively. In conclusion, this research preliminarily confirmed that yak Smad 4 gene was the target gene of bta-miR-146a, suggesting that bta-miR-146a might play roles in the development of yak follicles by regulating Smad 4 gene.
Key words: yak     bta-miR-146a     Smad 4     luciferase activity    

微小RNA(microRNA, miRNA)是一类内源性非蛋白质编码的RNA分子,约含有18~24个核苷酸,其一般通过作用于靶基因的3′UTR引起靶基因降解或者翻译受阻,从而调控基因表达[1-2]。1993年Lee等[3]在秀丽隐杆线虫中首次发现微小RNA-Lin-4,它能通过与其3′UTR结合,调控线虫细胞的发育。从此越来越多的学者展开对miRNA的相关研究。迄今为止,已通过克隆和测序技术分析出人[4]、小鼠[5-6]、猪[7]、牛[8]、羊[9]等多个物种的卵巢组织miRNA的表达谱,而且对颗粒细胞[10]、卵母细胞[11]、卵泡液[12]和黄体[13]等特定卵巢组织的miRNA的研究也有重大发现。Otsuka等[14-16]在小鼠卵巢中敲除miRNA生成过程的关键酶Dicerl(dicerl hypomorphic)后,导致雌鼠不育;Medeiros等[17]研究发现,miR-290-295缺陷会引起原始生殖细胞的丢失,从而导致雌性小鼠不育;Yin等[18]研究发现,miR-383通过靶定RNA结合基序单链结合蛋白1(RNA binding motif, single stranded interacting protein 1, RBMS 1)基因能促进小鼠颗粒细胞增殖和激素分泌;Pangas等[19-20]研究发现,miR-224通过调控Smad 4基因表达,从而促进小鼠颗粒细胞增殖和芳香酶表达;Yan等[21]研究发现,miR-145通过靶定细胞周期素D2(cyclinB 1, Ccnd 2)基因抑制小鼠颗粒细胞的增殖;Lin等[22]在猪卵巢卵泡中发现,miR-26b通过靶向DNA损伤相关基因的铵转运蛋白(ammonium transporter, AMT)基因,从而调控颗粒细胞凋亡;Ma和Xu等[13, 23]在猪和牛颗粒细胞中均发现miR-378直接靶向芳香酶,从而调控雌二醇的分泌,这些研究进一步表明,miRNA在动物卵巢卵泡中发挥着不可替代的作用[24-26]。前期的研究发现,bta-miR-146a在牦牛心、肝、脾、肺、肾、淋巴结、卵巢等组织中均有表达,其中在卵巢组织中的表达量较高[27]。然而关于bta-miR-146a靶基因的研究鲜见报道。

因此,本试验通过生物信息学软件分析bta-miR-146a的靶基因并从中选取卵巢卵泡发育相关的靶基因进行靶基因验证,为bta-miR-146a在牦牛卵巢卵泡发育中可能发挥的功能机制提供理论基础。

1 材料与方法 1.1 主要的试剂与仪器

293T细胞株(广州吉妮欧生物科技有限公司);限制性内切酶(Thermo Fisher);无内毒素质粒提取试剂盒(D6915,Omega);连接酶(EL0011,Thermo Fisher);PCR产物回收试剂盒(D6492,Omega);KOD酶(KOD-FX,Toyobo);DMEM(11965-092,Gibico);Fetal Bovine Serum(10099-141,Gibico);Trypsin-EDTA (25300-062,Gibico);Lipofectamine 2000 Reagent(11668-019,Thermo Fisher);Dual-Luciferase Reporter Assay System(E1910,Promega)。

荧光倒置显微镜(1X71,OLYMPUS);离心机(3K15,Sigma);电泳槽(JK-SPFT,北京君意东方电泳设备有限公司);电泳仪(JY300C,北京君意东方电泳设备有限公司);二氧化碳培养箱(NCO15AC,Sanyo);凝胶成像系统(JY04S-3C, 北京君意东方电泳设备有限公司);其他试剂均为国产分析纯。

1.2 生物信息软件预测及保守性分析 1.2.1 miR-146a在不同物种中的保守性分析

在miRBase数据库(http://www.mirbase.org)中查询并下载牛(Bos taurus)、人(Homo saspiens)、小鼠(Mus musculus)、大鼠(Rattus norvegicus)、黑猩猩(Gorilla gorilla)、猕猴(Macaca mulatta)、狗(dog)、鸡(chicken)等8个物种的miR-146a成熟序列,再运用MEGA6.0软件分析其保守性。

1.2.2 miR-146a靶基因预测

由于牦牛miRNA尚无合适的数据库,因此选择在牛miRNA数据库中进行靶基因预测。利用生物信息学软件TargetScan(http://www.targetscan.org)在线预测,以miR-146a为miRNA,预测其靶基因,再从中挑选可能与哺乳动物卵巢卵泡生长发育相关的靶基因为研究对象。

1.3 牦牛Smad 4基因3′UTR野生型载体构建 1.3.1 引物设计和PCR扩增

根据GenBank ID: ENST00000398417.2确定的牛Smad 4基因3′UTR区,再利用BLAST对基因序列及载体序列比对分析后, 发现Xho Ⅰ与Not Ⅰ两个限制性内切酶可以将目的片段克隆到载体中,因此选择psiCHECK2载体(图 1)。设计合成特异引物(F:5′-CCG CTCGAGGGTCTTTTAACTTAACGTTG-3′;R:5′-ATAAGAATGCGGCCGCCGG TAAGATATTAGATGGAG -3′),其中导入了XhoⅠ与NotⅠ内切酶酶切位点(引物斜体下划线部分),用于PCR扩增牦牛Smad 4基因的3′UTR,预期扩增目的片段总长度为1 123 bp。

图 1 载体图谱 Figure 1 Vector map

以牦牛卵巢组织基因组DNA为模板,进行PCR扩增,反应体系为50 μL:2× KOD-FX Buffer 25 μL,dNTP 10 μL,上下游引物(10 μmol·L-1)各1.5 μL,KOD-FX DNA聚合酶1 μL,DNA模板2 μL,ddH2O 9 μL。反应条件:95 ℃预变性5 min;95 ℃变性20 s, 60 ℃退火10 s,72 ℃延伸30 s,进行40个循环;72 ℃延伸7 min;4 ℃保存。PCR产物进行1%琼脂糖电泳检测确定特异性产物。

1.3.2 扩增产物回收

按照Omega回收试剂盒说明书操作:200 μL PCR反应体积加入800 μL Buffer CP到1.5 mL离心管,剧烈震荡,短暂离心;将上述混合物转入吸附柱,13 000 r·min-1离心1 min,弃滤液;加入700 μL洗脱液,13 000 r·min-1离心1 min,弃滤液;加入500 μL洗脱液,13 000 r·min-1离心1 min,弃滤液;13 000 r·min-1离心2 min,甩去吸附柱上的残留液体;把吸附柱转移到一个新的1.5 mL离心管,在吸附柱中加入30 μL ddH2O,室温放置1 min;13 000 r·min-1离心2 min,剩余滤液即是回收的DNA。

1.3.3 酶切和产物纯化

XhoⅠ与NotⅠ双酶切Smad 4启动子基因的PCR产物和载体,酶切反应体系:10×FastDigest Green Buffer 2 μL,vector 10 μL,内切酶各1 μL,ddH2O补足到20 μL。反应条件:37 ℃ 4 h。酶切后按照上述1.3.2方法回收并纯化酶切产物。

1.3.4 连接、转化和菌落鉴定

取纯化酶切产物进行连接反应,反应总体系为20 μL:10×Reaction Buffer 2 μL,Target 10 μL,Vector 5 μL,T4 DNA Ligase 1 μL,ddH2O补足到20 μL。反应条件:16 ℃ 30 min。将连接产物加到100 μL Top10感受态细胞中混匀,冰浴30 min后置于42 ℃水浴90 s,取出后立即于冰浴中放置2 min, 加入800 μL 37 ℃预热的LB液体培养基(不含抗生素),37 ℃振荡培养45 min复苏;8 000 r·min-1离心1 min,弃上清液,将100 μL混匀菌液加到含抗生素LB固体琼脂培养基上,倒置于37 ℃培养箱中,培养12~16 h长出菌落。

挑取上述平板上的单菌落,溶到3 μL水中,取0.5 μL做模板,按照上述步骤进行PCR扩增。试验共挑取8个菌落,用载体上游和下游引物进行PCR反应,再用1%琼脂糖电泳检测扩增产物,然后将阳性克隆产物送北京信诺金达生物科技有限公司进行测序鉴定。

1.4 牦牛Smad 4基因3′-UTR突变型载体构建

利用突变PCR的方法, 将牦牛Smad 4基因3′UTR靶序列AGTTCTC突变为CTGGAGA。突变引物F:5′-TGGTTAATAGTACATACTGTGAGAGCAAATGAA-3′;R:5′-CTCACAGTATGTACTATTAACCAGCTGCACAAC-3′。Smad 4基因3′UTR双荧光素酶报告基因突变型载体(以下简称“Smad4-3′UTR-MT载体”)构建的基本过程:以上述野生型载体(以下简称“Smad4-3′UTR-WT载体”)为模板, 分别扩增含突变序列的两段序列, 并进行拼接, 之后进行酶切、连接和转化,最后突变载体同样进行测序分析和扩增。

1.5 细胞培养与转染

复苏293T细胞,用细胞培养液(DMEM+10% FBS)培养于37 ℃、5%CO2、95%相对湿度的培养箱中,经2~3次传代,细胞生长状态较好并呈对数生长期时,开始准备转染。在96孔板中,每孔接种1×104个细胞,每孔加100 μL细胞培养液,放入37 ℃、5% CO2、95%相对湿度的培养箱中培养24 h。然后按照Lipofectamine 2000转染试剂说明书进行转染:10 μL OPTI-MEM培养基稀释mimic或mimic NC,15 μL OPTI-MEM培养基稀释载体,25 μL OPTI-MEM培养基稀释lipo2000转染试剂,室温放置5 min。将三者轻轻吹打混合,静置20 min后,加入至细胞内。转染6 h后更换为新鲜细胞培养液。bta-miR-146a mimics/NC均由广州市锐博生物科技有限公司合成。mimics转染浓度为20 μmol·L-1,质粒浓度为100 ng。总共设4个试验组,分别为Smad4-3′UTR-WT与mimic共转染组、Smad4-3′UTR-WT与mimic NC共转染组、Smad4-3′UTR-MT与mimic共转染组、Smad4-3′UTR-MT与mimic NC共转染组,每组4个重复。

1.6 荧光素酶活性测定

转染48 h后,各组细胞间用细胞裂解液裂解细胞,取5 μL细胞裂解液与萤火虫荧光素酶缓冲液及底物5 μL,混合测定荧光强度,然后加入海肾荧光素酶反应缓冲液及腔肠素底物5 μL,混匀再次测定海肾荧光素酶活性。每个样品处理一致,比较海肾荧光素酶活性。

1.7 统计学方法

用SPSS19.0统计软件中单因素方差分析对试验组荧光素酶活性进行差异显著分析,试验统计数据均以“平均值±标准差(x±Sd)”表示。

2 结果 2.1 bta-miR-146a保守性分析

利用miRBase数据库(http://www.mirbase.org)中下载的牛(Bos taurus)、人(Homo saspiens)、小鼠(Mus musculus)、大鼠(Rattus norvegicus)、黑猩猩(Gorilla gorilla)、猕猴(Macaca mulatta)、狗(dog)、鸡(chicken)8个物种的miR-146a成熟序列,再运用MEGA6.0软件分析其保守性(表 1),结果发现,miR-146a的成熟序列在8个物种均为22个碱基,碱基基本一致,进一步表明miR-146a成熟序列在脊椎动物中具有高度的保守性。

表 1 miR-146a在不同物种中的成熟序列比对分析 Table 1 The alignment of mature miR-146a sequences from 8 sepcies
2.2 bta-miR-146a靶基因预测

采用TargetScan软件,以牛数据库为基础,搜索miR-146a的靶基因,结果发现,bta-miR-146a存在248个潜在靶基因,其中有178个潜在靶基因在卵巢组织中表达(表 2),进一步分析发现,Smad 4基因也是bta-miR-146a的一个潜在靶基因。前期研究发现,牦牛Smad 4基因是一个表达谱很广的基因,尤其在卵巢、输卵管、子宫等组织中均有mRNA表达,而且牛Smad 4基因3′UTR与bta-miR-146a在397~404 bp位置存在互补结合(图 2),因此预测牦牛Smad 4基因可能是bta-miR-146a的靶基因之一。

表 2 miR-146a靶基因相关信息 Table 2 Related informations of bta-miR-146a target genes
图 2 Smad 4基因3′UTR和bta-miR-146a相互结合靶位点预测 Figure 2 Target point prediction for bovine Smad 4 gene 3′UTR to bta-miR-146a
2.3 牦牛Smad 4-3′UTR-WT载体构建

利用设计的特异性引物,以牦牛卵巢组织基因组DNA为模板PCR扩增Smad 4基因3′UTR,并用1%琼脂糖凝胶电泳检测扩增产物,然后对经过酶切、连接和转化后的Smad 4基因3′UTR载体菌落进行PCR扩增,用1%琼脂糖检测扩增产物。由图 3图 4可知,PCR扩增产物在1 000~2 000 bp之间出现特异性条带,并且其大小与预期片段大小一致。将阳性克隆产物送公司测序,测序结果(图 5)证实牦牛Smad 4基因3′UTR靶序列AGTTCTC成功地插入双荧光素酶报告载体psiCHECK2中,此Smad 4-3′UTR-WT载体可以用于荧光素酶活性的检测。

M. DNA相对分子质量标准;1~2. PCR扩增产物 N. DNA marker DL2000;1-2. PCR products of yak Smad 4 gene 3′UTR 图 3 牦牛Smad 4基因3′UTR区RT-PCR电泳图 Figure 3 RT-PCR electrophoresis of yak Smad 4 gene 3′UTR
M. DNA相对分子质量标准;1~8.菌落扩增产物 N. DNA marker DL2000;1-8. Colonies amplification products 图 4 牦牛Smad 4基因载体菌落RT-PCR电泳图 Figure 4 RT-PCR electrophoresis of yak Smad 4 vector colonies
A. Smad 4-3′UTR-WT载体序列;B. Smad 4-3′UTR-MT载体序列 A. Sequence of Smad 4-3′UTR-WT; B. Sequence of Smad 4-3′UTR-MT 图 5 牦牛Smad 4基因3′UTR构建载体序列测序图 Figure 5 Yak Smad 4 gene 3′UTR vector sequencing map
2.4 牦牛Smad 4-3′UTR-MT载体构建

以牦牛Smad 4-3′UTR-WT载体为模板,进行PCR反应,扩增含突变序列的两段序列并进行拼接,之后进行酶切、连接、转化、提取质粒,并进行测序分析,测序结果证实牦牛Smad 4基因3′UTR靶序列AGTTCTC突变为CTGGAGA(图 5),此Smad 4-3′UTR-MT载体构建成功,可以用于后续荧光素酶活性的检测。

2.5 双荧光素酶活性检测结果

将牦牛Smad 4基因3′UTR野生型/突变型载体与bta-miR-146a mimics/NC共转染293T细胞后,检测双荧光素酶活性。由图 6可知,Smad 4-3′UTR-WT+ mimic组的荧光素酶活性显著低于Smad 4-3′UTR-WT+ mimic NC组(P < 0.05),而Smad 4-3′UTR-MT+ mimic组的荧光素酶活性与Smad 4-3′UTR-MT+ mimic NC组差异不显著(P>0.05)。此结果表明,bta-miR-146a对Smad 4-3′UTR-WT的荧光素酶活性有显著抑制作用,而bta-miR-146a对Smad 4-3′UTR-MT的荧光素酶活性没有显著影响。

3 讨论

miRNA是一类内源性的非编码RNA,大小为18~24个核苷酸,主要通过碱基互补配对的方式识别靶基因,并根据互补程度的不同指导沉默复合体降解靶基因或者阻遏靶基因的翻译[25]。有研究表明,miRNA对卵泡和卵母细胞的生长、成熟有不可替代的作用[26]。在小鼠的研究中发现,miR-383能促进颗粒细胞增殖、激素分泌[18];miR-224能促进颗粒细胞增殖[28-29]、激素分泌以及芳香酶表达[30],延迟卵丘卵母细胞成熟[31];miR-145能抑制颗粒细胞的增殖[21];miR-146能经过视黄酸信号途径调控精原细胞分化[32]。在猪的研究中发现,miR-26b抑制芳香酶表达和激素分泌,并且促进颗粒细胞的凋亡[22]。我们先前通过PCR半定量方法检测了Smad 4基因和miRNA在牦牛各组织中的表达情况,结果发现,Smad 4在牦牛卵巢、输卵管和子宫等生殖组织中均有表达,而且bta-miR-146a在牦牛卵巢组织中表达量较高,推测Samd 4基因是bta-miR-146a的靶基因[27]。因此本试验首先利用miRBase数据库和MEGA 6.0分析了miR-146a成熟序列在牛、人、小鼠、大鼠、黑猩猩、猕猴、狗和鸡等8个物种中的保守性,结果发现,miR-146a成熟序列在脊椎动物中具有高度的保守性,暗示miR-146a在基本生命过程中发挥着一定的生物学功能。然后采用TargetScan软件,以牛数据库为基础,搜索miR-146a的靶基因并且分析bta-miR-146a与靶基因的结合位点,结果发现,Smad 4是bta-miR-146a的248个潜在靶基因之一,而且Smad 4基因3′UTR区AGUUCUC序列与bta-miR-146a成熟序列UCAAGAG存在保守的互补结合位点,更进一步提示Smad 4基因可能是bta-miR-146a的靶基因之一。

各组间,不同字母表示差异显著(P < 0.05),相同字母表示差异不显著(P>0.05) Different letters indicate significant differences(P < 0.05), the same letter indicate no significant difference(P>0.05) among different groups 图 6 双荧光素酶活性检测结果 Figure 6 The results of dual-luciferase activity assay

在哺乳动物中,Smads蛋白是将转化生长因子-β(transforming growth factor-β,TGF-β)信号传导至细胞核的介导蛋白,在TGF-β信号传导过程中起着关键性作用[33-34]。目前已发现9种哺乳动物的Smads蛋白,根据功能不同分为受体激活型Smad蛋白(recepter activated smad,R-Smad)、抑制型Smad蛋白(inhibitory smads,I-Smads)和通用调节型Smad蛋白(common mediator smads,Co-Smads),其中Smad 4为唯一Co-Smads蛋白[35-36]。在TGF-β信号传导过程中,R-Smad首先被激活的受体磷酸化,然后与Smad 4结合形成异源性聚合体,进入细胞核,从而调控目的基因的表达[33, 37]。当敲除Smad 4基因时会阻断信号传导通路,使上游信号分子所携带的信息无法正常传递至细胞核,从而导致细胞损伤[19, 38-39]。故Smad 4在TGF-β信号传导过程中起着中枢性作用[3, 33, 37]。而且现研究也已证明,Smad 4基因在小鼠多种组织中发挥着重要作用,例如,Smad 4基因在卵巢中表达广泛,其中在卵泡中表达更加明显,而且随着卵巢的成熟,表达量逐渐上升[40];当敲除心内膜中Smad 4基因会导致上皮细胞转化功能受阻、心内膜细胞增殖停止和心室心内膜结构异型等症状,从而影响心正常生长发育[41];敲除胚胎平滑肌细胞Smad 4基因会对平滑肌细胞分化、增殖、迁移等产生影响,从而导致胚胎心血管机能不全[42];当敲除表皮Smad 4基因会使毛囊干细胞持续性激活消耗,从而导致皮肤肿瘤形成[43];当特异性敲除Smad 4基因后小鼠排卵数、产仔数等生殖能力随时间而逐渐下降,并且有腔卵泡明显减少,闭锁卵泡大大增加,表明Smad 4基因在小鼠卵巢卵泡发育过程中确实发挥着关键性作用[19, 39]

为了验证牦牛Smad 4基因与bta-miR-146a的靶向关系,本研究以牦牛卵巢基因组DNA为模板,克隆Smad 4基因3′UTR区,并将结合序列AGTTCTC和突变序列CTGGAGA插入双荧光素酶报告载体psiCHECK2中,从而构建Smad4-3′UTR-WT和Smad4-3′UTR-MT载体,再通过载体与miRNA共转染细胞后,检测双荧光素酶活性,结果发现,bta-miR-146a可以直接与Smad 4基因3′UTR区的部分碱基互补配对结合,表明Smad 4基因是bta-miR-146a的靶基因。这在理论上直接证明bta-miR-146a可以抑制Samd 4基因表达,此结果为进一步探索bta-miR-146a在牦牛卵巢卵泡发育中的功能奠定了理论基础。

4 结论

本试验先利用miRBase数据库和MEGA 6.0软件分析发现miR-146a具有高度的保守性,再利用TargetScan软件成功预测到Smad 4基因是bta-miR-146a的潜在靶基因之一,而且其3′UTR区部分序列与bta-miR-146a存在结合位点。然后通过成功构建Smad 4基因3′UTR区双荧光素酶报告基因载体,初步证实Smad 4基因是bta-miR-146a的靶基因,从而提示bta-miR-146a在牦牛卵泡发育过程中可能通过调控Smad 4基因发挥作用,但有待深入研究。

参考文献
[1] MOHR A M, MOTT J L. Overview of microRNA biology[J]. Semin Liver Dis, 2015, 35(1): 3–11. DOI: 10.1055/s-00000069
[2] DONG H F, LEI J P, DING L, et al. microRNA:Function, detection, and bioanalysis[J]. Chem Rev, 2013, 113(8): 6207–6233. DOI: 10.1021/cr300362f
[3] LEE R C, FEINBAUM R L, AMBROS V. The C.elegans heterochronic gene lin-4 encodes small RNAs with antisense complementarity to lin-14[J]. Cell, 1993, 75(5): 843–854. DOI: 10.1016/0092-8674(93)90529-Y
[4] LANDGRAF P, RUSU M, SHERIDAN R, et al. A mammalian microRNA expression atlas based on small RNA library sequencing[J]. Cell, 2007, 129(7): 1401–1414. DOI: 10.1016/j.cell.2007.04.040
[5] AHN H W, MORIN R D, ZHAO H, et al. microRNA transcriptome in the newborn mouse ovaries determined by massive parallel sequencing[J]. Mol Hum Reprod, 2010, 16(7): 463–471. DOI: 10.1093/molehr/gaq017
[6] RO S, SONG R, PARK C, et al. Cloning and expression profiling of small RNAs expressed in the mouse ovary[J]. RNA, 2007, 13(12): 2366–2380. DOI: 10.1261/rna.754207
[7] LI M Z, LIU Y K, WANG T, et al. Repertoire of porcine microRNAs in adult ovary and testis by deep sequencing[J]. Int J Biol Sci, 2011, 7(7): 1045–1055. DOI: 10.7150/ijbs.7.1045
[8] HOSSAIN M M, GHANEM N, HOELKER M, et al. Identification and characterization of miRNAs expressed in the bovine ovary[J]. BMC Genomics, 2009, 10: 443. DOI: 10.1186/1471-2164-10-443
[9] MCBRIDE D, CARRÉ W, SONTAKKE S D, et al. Identification of miRNAs associated with the follicular-luteal transition in the ruminant ovary[J]. Reproduction, 2012, 144(2): 221–233. DOI: 10.1530/REP-12-0025
[10] FIEDLER S D, CARLETTI M Z, HONG X M, et al. Hormonal regulation of MicroRNA expression in periovulatory mouse mural granulosa cells[J]. Biol Reprod, 2008, 79(6): 1030–1037. DOI: 10.1095/biolreprod.108.069690
[11] PAN B, TOMS D, SHEN W, et al. microRNA-378 regulates oocyte maturation via the suppression of aromatase in porcine cumulus cells[J]. Am J Physiol Endocrinol Metab, 2015, 308(6): E525–E534. DOI: 10.1152/ajpendo.00480.2014
[12] DA SILVEIRA J C, VEERAMACHANENI D N R, WINGER Q A, et al. Cell-secreted vesicles in equine ovarian follicular fluid contain miRNAs and proteins:A possible new form of cell communication within the ovarian follicle[J]. Biol Reprod, 2012, 86(3): 71.
[13] MA T H, JIANG H, GAO Y, et al. Microarray analysis of differentially expressed microRNAs in non-regressed and regressed bovine corpus luteum tissue; microRNA-378 may suppress luteal cell apoptosis by targeting the interferon gamma receptor 1 gene[J]. J Appl Genet, 2011, 52(4): 481–486. DOI: 10.1007/s13353-011-0055-z
[14] OTSUKA M, JING Q, GEORGEL P, et al. Hypersusceptibility to vesicular stomatitis virus infection in Dicer1-deficient mice is due to impaired miR24 and miR93 expression[J]. Immunity, 2007, 27(1): 123–134. DOI: 10.1016/j.immuni.2007.05.014
[15] OTSUKA M, ZHENG M, HAYASHI M, et al. Impaired microRNA processing causes corpus luteum insufficiency and infertility in mice[J]. J Clin Invest, 2008, 118(5): 1944–1954. DOI: 10.1172/JCI33680
[16] HONG X M, LUENSE L J, MCGINNIS L K, et al. Dicer1 is essential for female fertility and normal development of the female reproductive system[J]. Endocrinology, 2008, 149(12): 6207–6212. DOI: 10.1210/en.2008-0294
[17] MEDEIROS L A, DENNIS L M, GILL M E, et al. Mir-290-295 deficiency in mice results in partially penetrant embryonic lethality and germ cell defects[J]. Proc Natl Acad Sci U S A, 2011, 108(34): 14163–14168. DOI: 10.1073/pnas.1111241108
[18] YIN M M, LV M R, YAO G D, et al. Transactivation of microRNA-383 by steroidogenic factor-1 promotes estradiol release from mouse ovarian granulosa cells by targeting RBMS1[J]. Mol Endocrinol, 2012, 26(7): 1129–1143. DOI: 10.1210/me.2011-1341
[19] PANGAS S A, LI X H, ROBERTSON E J, et al. Premature luteinization and cumulus cell defects in ovarian-specific Smad4 knockout mice[J]. Mol Endocrinol, 2006, 20(6): 1406–1422. DOI: 10.1210/me.2005-0462
[20] YAO G D, YIN M M, LIAN J, et al. microRNA-224 is involved in transforming growth factor-β-mediated mouse granulosa cell proliferation and granulosa cell function by targeting Smad4[J]. Mol Endocrinol, 2010, 24(3): 540–551. DOI: 10.1210/me.2009-0432
[21] YAN G J, ZHANG L X, FANG T, et al. microRNA-145 suppresses mouse granulosa cell proliferation by targeting activin receptor IB[J]. FEBS Lett, 2012, 586(19): 3263–3270. DOI: 10.1016/j.febslet.2012.06.048
[22] LIN F, LI R, PAN Z X, et al. miR-26b promotes granulosa cell apoptosis by targeting ATM during follicular atresia in porcine ovary[J]. PLoS One, 2012, 7(6): e38640. DOI: 10.1371/journal.pone.0038640
[23] XU S Y, LINHER-MELVILLE K, YANG B B, et al. Micro-RNA378(miR-378) regulates ovarian estradiol production by targeting aromatase[J]. Endocrinology, 2011, 152(10): 3941–3951. DOI: 10.1210/en.2011-1147
[24] BARTEL D P. microRNAs:Genomics, biogenesis, mechanism, and function[J]. Cell, 2004, 116(2): 281–297. DOI: 10.1016/S0092-8674(04)00045-5
[25] AMBROS V. The functions of animal microRNAs[J]. Nature, 2004, 431(7006): 350–355. DOI: 10.1038/nature02871
[26] 徐源, 孙铁成, 张爱玲, 等. 卵巢miRNA研究进展[J]. 畜牧兽医学报, 2014, 45(4): 509–516.
XU Y, SUN T C, ZHANG A L, et al. Progress on the research of ovarian miRNA[J]. Acta Veterinaria et Zootechnica Sinica, 2014, 45(4): 509–516. (in Chinese)
[27] 王玉恒, 徐业芬, 牛家强, 等. 牦牛Smad 4基因mRNA组织表达谱及其miRNAs初步研究[J]. 畜牧兽医学报, 2017, 48(10): 1882–1891.
WANG Y H, XU Y F, NIU J Q, et al. Preliminary studies on tissue expression profiles of Smad 4 mRNA and its miRNAs in yak[J]. Acta Veterinaria et Zootechnica Sinica, 2017, 48(10): 1882–1891. (in Chinese)
[28] LIANG M, YAO G D, YIN M M, et al. Transcriptional cooperation between p53 and NF-κB p65 regulates microRNA-224 transcription in mouse ovarian granulosa cells[J]. Mol Cell Endocrinol, 2013, 370(1-2): 119–129. DOI: 10.1016/j.mce.2013.02.014
[29] YAO G D, LIANG M, LIANG N, et al. microRNA-224 is involved in the regulation of mouse cumulus expansion by targeting Ptx3[J]. Mol Cell Endocrinol, 2014, 382(1): 244–253. DOI: 10.1016/j.mce.2013.10.014
[30] WANG Y, REN J W, GAO Y, et al. microRNA-224 targets SMAD family member 4 to promote cell proliferation and negatively influence patient survival[J]. PLoS One, 2013, 8(7): e68744. DOI: 10.1371/journal.pone.0068744
[31] LI X F, WANG H D, SHENG Y, et al. microRNA-224 delays oocyte maturation through targeting Ptx3 in cumulus cells[J]. Mech Dev, 2017, 143: 20–25. DOI: 10.1016/j.mod.2016.12.004
[32] HUSZAR J M, PAYNE C J. MicroRNA 146(Mir146) modulates spermatogonial differentiation by retinoic acid in mice[J]. Biol Reprod, 2013, 88(1): 15.
[33] HELDIN C H, MIYAZONO K, TEN DIJKE P. TGF-β signalling from cell membrane to nucleus through SMAD proteins[J]. Nature, 1997, 390(6659): 465–471. DOI: 10.1038/37284
[34] LEE K B, ZHANG K, FOLGER J K, et al. Evidence supporting a functional requirement of SMAD4 for bovine preimplantation embryonic development:A potential link to embryotrophic actions of follistatin[J]. Biol Reprod, 2014, 91(3): 62.
[35] LUUKKO K, YLIKORKALA A, MÄKELÄ T P. Developmentally regulated expression of Smad3, Smad4, Smad6, and Smad7 involved in TGF-beta signaling[J]. Mech Dev, 2001, 101(1-2): 209–212. DOI: 10.1016/S0925-4773(00)00556-6
[36] WONG C, ROUGIER-CHAPMAN E M, FREDERICK J P, et al. Smad3-Smad4 and AP-1 complexes synergize in transcriptional activation of the c-Jun promoter by transforming growth factor β[J]. Mol Cell Biol, 1999, 19(3): 1821–1830. DOI: 10.1128/MCB.19.3.1821
[37] ITOH S, TEN DIJKE P. Negative regulation of TGF-β receptor/Smad signal transduction[J]. Curr Opin Cell Biol, 2007, 19(2): 176–184. DOI: 10.1016/j.ceb.2007.02.015
[38] LI J, QIN Y, ZHAO F K, et al. Anterior segment dysgenesis correlation with epithelial-mesenchymal transition in Smad4 knockout mice[J]. Int J Ophthalmol, 2016, 9(7): 943–947.
[39] YU C, ZHANG Y L, FAN H Y. Selective Smad4 knockout in ovarian preovulatory follicles results in multiple defects in ovulation[J]. Mol Endocrinol, 2013, 27(6): 966–978. DOI: 10.1210/me.2012-1364
[40] HUANG Z, YUAN X, WANG M, et al. Molecular cloning of the SMAD4 gene and its mRNA expression analysis in ovarian follicles of the Yangzhou goose (Anser cygnoides)[J]. Br Poult Sci, 2016, 57(4): 515–521. DOI: 10.1080/00071668.2016.1180670
[41] MOSKOWITZ I P, WANG J, PETERSON M A, et al. Transcription factor genes Smad4 and Gata4 cooperatively regulate cardiac valve development.[corrected][J]. Proc Natl Acad Sci U S A, 2011, 108(10): 4006–4011. DOI: 10.1073/pnas.1019025108
[42] MAO X, DEBENEDITTIS P, SUN Y, et al. Vascular smooth muscle cell Smad4 gene is important for mouse vascular development[J]. Arterioscler Thromb Vasc Biol, 2012, 32(9): 2171–2177. DOI: 10.1161/ATVBAHA.112.253872
[43] YANG L L, MAO C M, TENG Y, et al. Targeted disruption of Smad4 in mouse epidermis results in failure of hair follicle cycling and formation of skin tumors[J]. Cancer Res, 2005, 65(19): 8671–8678. DOI: 10.1158/0008-5472.CAN-05-0800