Functional analysis of MSTN promoter in scallop (Chalmys nobilis)
-
摘要:
为了解肌肉生长抑制素myostatin (MSTN)在华贵栉孔扇贝(Chalmys nobilis)闭壳肌肌肉生长和发育过程中所起的负调控作用，对华贵栉孔扇贝的MSTN启动子序列进行了生物信息学分析。结果显示，MSTN启动子序列长1 358 bp，有4个转录起始位点。核心启动子区为–100~–51 bp。有1个TATA-box (–92~–86 bp)和2个E-box等顺式作用元件；潜在的转录因子结合位点有MEF2、MEF3、FoxO、MTBF和MyoD等；启动子区域无CpG岛。成功构建了6个MSTN启动子不同长度片段的荧光素酶表达载体，瞬时转染到293T细胞并进行双荧光素酶报告基因活性检测，表明6个启动子片段均有转录活性，PGL-534的活性最高，其次为PGL-274、PGL-22和PGL-102，最低的为PGL-995。–216~–364区域可能存在负调控基因表达的转录因子结合位点，–364~–825区域可能存在正调控基因表达的转录因子结合位点。
Abstract:In order to investigate the negative regulater role of myostatin (MSTN) in growth and development of adductor of Chalmys nobilis, we analyzed the promoter sequence of MSTN of C. nobilis. The results show that the promoter sequence was 1 358 bp in length, including four transcription start sites (TSS). The core promoter region was located from –100 bp to –51 bp. Two kinds of cis-regulatory element, a TATA-box (located from –92 bp to –86 bp) and two E-boxes, were detected in promoter. Potential transcription factor binding sites including MEF2, MEF3, FoxO, MTBF, MyoD and so on were found in promoter. No CpG island was found. Six luciferase expression vectors with different lengths of MSTN promoter were successfully constructed and transiently transfected into 293T cells for an analysis of the activity of dual luciferase reporter gene. It is shown that all the six promoter sequences had transcriptional activity, with PGL-534 the highest, followed by PGL-274, PGL-22 and PGL-102, and PGL-995 the lowest. There might be some binding sites of potential negative transcription factors from –216 bp to –364 bp, and potential positive transcription factors from –364 bp to –825 bp.
-
Keywords:
- Chalmys nobilis /
- MSTN /
- promoter /
- transient transfection /
- functional analysis
-
华贵栉孔扇贝(Chalmys nobilis)具有生长快、产量高、养殖周期短等特点，是中国华南沿海地区主要的海水养殖贝种之一，具有重要的经济价值。闭壳肌是华贵栉孔扇贝主要的食用部分。目前对扇贝闭壳肌的研究较少，主要有闭壳肌与其他性状的相关性研究，如墨西哥湾扇贝(Argopecten irradians concentricus)体质量和壳长与闭壳肌质量紧密相关[1]；华贵栉孔扇贝软体部质量对闭壳肌质量的直接影响最大，其次是壳宽[2]；影响虾夷扇贝(C. farreri)闭壳肌质量的主要因素是壳宽[3]；另外闭壳肌颜色和蛋白分布也有涉及[4-5]。但很少见到闭壳肌生长和发育的相关研究。
肌肉生长抑制素myostatin (MSTN)是转化生长因子-β (transforming growth factor-beta)超家族的成员，在动物肌肉生长和发育过程中起负调控作用[6]。Myostatin蛋白包括1个分泌信号序列、蛋白水解处理位点、含有9个半胱氨酸残基的保守型羧基末端区[7-8]。Myostatin的前体蛋白由信号序列、N-末端前肽区域和1个活性配基的C-末端区域组成[7-8]。前体蛋白经过2次蛋白质水解切割后，Myostatin即可活化[8-9]。成熟的Myostatin蛋白是一个C-末端由二硫键连接的二聚体，与MSTN基因受体结合，从而发挥其生物活性[8-9]。
目前在海湾扇贝(A. irradians)、栉孔扇贝 (C. farreri)、华贵栉孔扇贝、贻贝(Mytilus chilensis)、小狮爪海扇蛤(Nodipecten subnodosus)和竹蛏(Sinonovacula constricta)等贝类中有myostatin基因的研究报道[10-16]。Morelos等[14]研究1龄小狮爪海扇蛤一年中myostatin在闭壳肌中的表达情况、闭壳肌质量和肌纤维数量和大小，证实myostatin表达和闭壳肌性状有关联。在一些贝类中发现MSTN有与生长性状相关联的SNP位点[15,17-18]。Hu等[11]预测了栉孔扇贝MSTN启动子中有MEF2、COMP、MTBF、E-box等调控元件。
笔者已获得了华贵栉孔扇贝MSTN启动子序列[16]。本研究拟对启动子序列进行生物信息学分析。通过构建不同缺失长度的报告基因系统，检测启动子的转录活性，为今后进一步研究MSTN基因的调控机制奠定基础。
1. 材料与方法
1.1 实验材料
华贵栉孔扇贝采购于广东深圳南澳，剪取其闭壳肌组织，液氮冷冻后放置于– 80 ℃冰箱保存。
1.2 实验方法
1.2.1 DNA提取
用海洋动物DNA提取试剂盒(天根)提取DNA，终质量浓度为50 ng·mL–1，– 20 ℃保存备用。
1.2.2 myostatin启动子的生物信息学分析
笔者已获得myostatin基因启动子序列[16]。用在线软件BDGP (http://www.fruitfly.org/seq_tools/promoter.html)预测核心启动子区域和转录起始位点(transcription start site，TSS)。用MatInspector软件预测潜在的转录因子结合位点。用MethPrimer软件(http://www.urogene.org/methprimer2/)预测CpG岛。
1.2.3 不同长度片段表达载体构建
根据潜在的转录因子结合位点分析结果，设计6个正向引物和1个反向引物。上、下游引物分别添加KpnⅠ和HindⅢ酶切位点和保护碱基(表1)。以PGL-22为正向引物，P1322A为反向引物，扩增出最长的启动子片段，胶回收后，连接到pMD18-T载体，转化到DH5α大肠杆菌中，挑单克隆菌送上海生工生物工程有限公司测序。测序正确后，用质粒提取试剂盒(生工)提取质粒，以质粒DNA为模板，分别以PGL-102、PGL-274、PGL-534、PGL-995和PGL-1143为正向引物，P1322A为反向引物，通过PCR扩增不同片段长度的MSTN基因启动子片段，随后连接转化，测序验证序列正确后提取质粒。
表 1 实验所用引物Table 1. Primers used in this experiment引物
primer序列 (5'−3')
sequencePGL-22 CGGggtaccACACGGCGAAAAAATGGAGC PGL-102 CGGggtaccGGATGCCATAAATCAAAACCACAAC PGL-274 CGGggtaccAAAACGCCGCCAAACG PGL-534 CGGggtaccACTTCAGGCTGTATCGCAAAT PGL-995 CGGggtaccACCCGTTGGCAGCGTTCA PGL-1143 CGGggtaccTCTAAATGCTAACCCTTGTGCTG P1322A CTAaagcttTATAGCGGTTACGTTACAGATGGTT 注：小写字体为酶切位点 Note: Lowercase letters represent restriction enzyme loci. 用KpnⅠ和HindⅢ2种内切酶同时在37 ℃水浴酶切含有不同长度片段的pMD18-T重组质粒和pGL3-basic质粒，分别纯化后，用Solution I连接酶(TaKaRa) 16 ℃连接过夜，转化到DH5α大肠杆菌中，用去内毒素质粒提取试剂盒(OMEGA)提取重组质粒，测序以验证序列是否正确。
1.2.4 细胞培养
在24孔培养板中接种HEK293T细胞，使转染时细胞密度在70%~80%。所用培养基为DMEM高糖培养基(Hyclone)，加入10% FBS (Gibco)。细胞放置于37 ℃、5% CO2培养箱中培养约20 h。
1.2.5 瞬时转染
内参质粒为pRL-TK (Promega)，分别将各pGL3-basic重组质粒和内参质粒以20∶1的比例共转染到细胞。转染对照组为pGL3-basic载体和pRL-TK载体。每组做5个平行。瞬时转染采用转染试剂盒Lipofectamine 2000 Transfection Reagent (Invitrogen)，根据试剂盒的说明进行实验。
1.2.6 双荧光素酶活性检测
用萤光素酶检测试剂盒Dual-Luciferase®Reporter Assay System(Promega)分析报告基因Luciferase活性。转染24 h后，吸去旧的培养基，加入1×PBS清洗1~2次，吸出PBS后，每孔中加入细胞裂解液PLB (passive lysis buffer，Promega) 200 μL，室温孵育10 min；将裂解液10 000 r·min–1离心5 min，取上清作为待测液；取100 μL裂解液上清加入96孔发光板中，加入100 μL萤火虫荧光素酶检测工作液(碧云天)，吹吸混匀；上机测发光值，积分时间5 s；加入海肾荧光素酶检测工作液(碧云天) 100 μL，吹吸混匀，用化学发光分析仪BK-L96C (中生北控)监测发光值。用SPSS 19.0软件中的单因素方差分析检测各样品间差异的显著性。
2. 结果
2.1 华贵栉孔扇贝MSTN基因启动子的生物信息学分析
华贵栉孔扇贝MSTN基因5'调控区序列长1 358 bp (GenBank登录号：KY888888)。BDGP预测有4个TSS，分别为位于ATG上游第70个碱基“A”，第233个碱基“G”，第1 059个碱基“A”和第1 299个碱基“T”(图1)。核心启动子区在–100~–51 bp。MatInspector在线软件分析潜在的转录因子结合位点包括MEF2、MEF3、FoxO、MTBF和MyoD等，存在1个TATA-box (位于–92~–86bp)和2个E-box顺式作用元件 (图1)。MethPrimer软件预测无CpG岛。
![]()
图 1 MSTN启动子序列分析下划直线表示转录因子结合位点，下划波浪线为顺式作用元件，斜体表示核心启动子区域，阴影部分表示转录起始位点Figure 1. Analysis of MSTN gene promoter sequenceThe potential transcription factor binding sites are straightly underlined. The cis-regulatory elements are wavily underlined. Core promoter region is in italic. Nucletide of each transcription start sites are in shadow.2.2 MSTN基因启动子不同长度片段的扩增和表达载体构建
6个正向引物(PGL-22、PGL-102、PGL-274、PGL-534、PGL-995和PGL-1143)和1个下游引物(P1322A)分别扩增出6个长度不同的启动子片段，长度依次为1 301 bp、1 221 bp、1 049 bp、789 bp、328 bp和180 bp (图2)。构建的pGL3-basic重组质粒名字与正向引物名字一一对应，即PGL-22、PGL-102、PGL-274、PGL-534、PGL-995和PGL-1143，经测序确定序列正确。
![]()
图 2 MSTN基因启动子不同长度片段的扩增电泳图电泳顺序：PGL-22、PGL-102、PGL-274、PGL-534、PGL-995、PGL-1143、DL2 000Figure 2. Electrophoretograms of MSTN gene promoter with different lengthsElectrophoresis order: PGL-22, PGL-102, PGL-274, PGL-534, PGL-995, PGL-1143, DL2 0002.3 MSTN启动子活性分析
将构建的6个不同长度缺失片段pGL3-basic重组质粒分别与pRL-TK质粒共转染到HEK293T细胞。pGL3-basic空载体作为阴性对照。荧光素酶活性检测结果见图3，6个启动子片段的pGL3-basic重组载体与阴性对照pGL-basic空载体相比，荧光强度均有显著性差异(P<0.05)。除PGL-22和PGL-274外，其他启动子片段的荧光强度彼此间均有显著性差异(P<0.05)。PGL-534的转录活性最高，其次为PGL-274、PGL-22和PGL-102，最低为PGL-995。
![]()
图 3 MSTN启动子不同长度片段的相对活性相同字母表示差异不显著，不同字母表示差异显著(P<0.05)Figure 3. Relative luciferase activity of MSTN gene with different lengthsThe same letters indicate no significant difference, while different letters indicate significant difference (P<0.05).3. 讨论
启动子是基因的“开关”，决定着基因的活动。通过生物信息学预测启动子上的顺式作用元件和反式作用因子，是研究启动子结构和功能的基础。
本研究对华贵栉孔扇贝MSTN启动子序列进行生物信息学分析，结果显示核心启动子区为–100~–51 bp。华贵栉孔扇贝MSTN启动子包含4个TSS，分别位于起始密码ATG上游的第70、第233、第1 059和第1 299个碱基。栉孔扇贝和人的MSTN启动子有3个TSS[11,19]，但在牛中只有1个[20]。TATA-box是真核生物启动子中重要组成部分，在猪[21]、大黄鱼 (Larimichthys crocea)[22]、栉孔扇贝[11]和华贵栉孔扇贝MSTN启动子中均发现有1个TATA-box，金头鲷 (Sparus aurata) MSTN启动子有2个TATA-box[23]。在大黄鱼和金头鲷MSTN启动子中，还检测到CAAT-box，但在扇贝中均未发现。E-box (核心序列为5'-CANNTG-3')在以上物种MSTN启动子中均检测到，说明E-box在不同物种MSTN启动子区域具有保守性。E-box突变会降低启动子活性[24]。E-box能被生肌调控家族因子(myogenic regular factor，MRFs)中的碱基螺旋-环-螺旋(bHLH)结构特异性识别，从而调控MSTN的转录[25]。
华贵栉孔扇贝MSTN启动子中存在的转录因子结合位点有MEF2、MEF3、FoxO、MTBF和MyoD。在其他物种中也发现有一样的转录因子结合位点，如在栉孔扇贝中发现MEF2和MTBF[11]，在金头鲷中发现MEF2[23]，在猪中发现MEF2和MyoD[21]，在牛中发现MEF2、MyoD和FoxO1[26]。均为肌肉特异性相关的转录因子结合位点，参与MSTN基因的转录和表达[27-28]。MEF2广泛存在于肌肉细胞中，可与许多肌肉特异性基因的启动子结合[29]。MEF2有MEF2A、MEF2B、MEF2C、MEF2D等4种亚型。
CpG岛主要存在于启动子或第一外显子中，在基因表达调控中起负调控作用。本研究未发现华贵栉孔扇贝MSTN启动子存在CpG岛，这与栉孔扇贝MSTN启动子的分析结果相同。
报告基因是检测启动子活性的常用方法。本研究构建了6个不同缺失片段的启动子报告基因，发现它们都有转录活性。这些报告基因均包含核心启动子区域，说明预测的核心启动子区域结果可靠。6个不同缺失片段启动子的转录活性与启动子长度无关。如PGL-534的活性最高，但该序列并不是最长。这说明MSTN启动子序列中有的序列是促进该基因表达的调控元件，有的则是抑制该基因表达的调控元件。这些表达元件在一起产生作用，产生不同的转录活性。如PGL-1143转录活性比PGL-995的强，说明–216~–364 bp区域可能存在负调控基因表达的转录因子结合位点；PGL-995的转录活性远低于PGL-534，说明–364~–825 bp区域可能存在正调控基因表达的转录因子结合位点。这些推论还需要进一步研究证实。
-
![]()
图 1 MSTN启动子序列分析
下划直线表示转录因子结合位点，下划波浪线为顺式作用元件，斜体表示核心启动子区域，阴影部分表示转录起始位点
Figure 1. Analysis of MSTN gene promoter sequence
The potential transcription factor binding sites are straightly underlined. The cis-regulatory elements are wavily underlined. Core promoter region is in italic. Nucletide of each transcription start sites are in shadow.
![]()
图 2 MSTN基因启动子不同长度片段的扩增电泳图
电泳顺序：PGL-22、PGL-102、PGL-274、PGL-534、PGL-995、PGL-1143、DL2 000
Figure 2. Electrophoretograms of MSTN gene promoter with different lengths
Electrophoresis order: PGL-22, PGL-102, PGL-274, PGL-534, PGL-995, PGL-1143, DL2 000
![]()
图 3 MSTN启动子不同长度片段的相对活性
相同字母表示差异不显著，不同字母表示差异显著(P<0.05)
Figure 3. Relative luciferase activity of MSTN gene with different lengths
The same letters indicate no significant difference, while different letters indicate significant difference (P<0.05).
表 1 实验所用引物
Table 1 Primers used in this experiment
引物
primer序列 (5'−3')
sequencePGL-22 CGGggtaccACACGGCGAAAAAATGGAGC PGL-102 CGGggtaccGGATGCCATAAATCAAAACCACAAC PGL-274 CGGggtaccAAAACGCCGCCAAACG PGL-534 CGGggtaccACTTCAGGCTGTATCGCAAAT PGL-995 CGGggtaccACCCGTTGGCAGCGTTCA PGL-1143 CGGggtaccTCTAAATGCTAACCCTTGTGCTG P1322A CTAaagcttTATAGCGGTTACGTTACAGATGGTT 注：小写字体为酶切位点 Note: Lowercase letters represent restriction enzyme loci. 
下载: 导出CSV
-
[1] 杨世平, 刘慧玲, 劳赞, 等. 墨西哥湾扇贝经济性状对闭壳肌重决定效应分析[J]. 渔业现代化, 2013, 40(3): 41-45. doi: 10.3969/j.issn.1007-9580.2013.03.009 [2] 刘志刚, 章启忠, 王辉. 华贵栉孔扇贝主要经济性状对闭壳肌重的影响效果分析[J]. 热带海洋学报, 2009, 28(1): 61-66. doi: 10.3969/j.issn.1009-5470.2009.01.010 [3] 张存善, 常亚青, 曹学彬, 等. 虾夷扇贝体形性状对软体重和闭壳肌重的影响效果分析[J]. 水产学报, 2009, 33(1): 87-94. [4] 刘志刚, 王辉, 吕文刚, 等. 华贵栉孔扇贝(Chlamys nobilis)壳色与闭壳肌颜色遗传规律的研究[J]. 海洋与湖沼, 2012, 43(2): 237-243. [5] 宋扬, 张晴, 周晏琳, 等. 虾夷扇贝横纹肌和平滑肌的蛋白分布及理化性质[J]. 水产学报, 2017, 41(11): 1798-1805. [6] MCPHERRON A C, LEE S J. Double muscling in cattle due to mutations in the myostatin gene[J]. Proc Natl Acad Sci USA, 1997, 94(23): 12457-12461. doi: 10.1073/pnas.94.23.12457
[7] DUSHYANTH K, BHATTACHARYA T K, SHUKLA R, et al. Gene expression and polymorphism of myostatin gene and its association with growth traits in chicken[J]. Anim Biotechnol, 2016, 27(4): 269-277. doi: 10.1080/10495398.2016.1182541
[8] QIAN Z, MI X, WANG X, et al. cDNA cloning and expression analysis of myostatin/GDF11 in shrimp, Litopenaeus vannamei[J]. Comp Biochem Physiol A, 2013, 165(1): 30-39. doi: 10.1016/j.cbpa.2013.02.001
[9] DALL'OLIO S, FONTANESI L, NANNI COSTA L, et al. Analysis of horse myostatin gene and identification of single nucleotide polymorphisms in breeds of different morphological types[J]. J Biomed Biotechnol, 2010, 2010: 542945.
[10] NÚÑEZ-ACUÑA G, GALLARDO-ESCÁRATE C. The myostatin gene of Mytilus chilensis evidences a high level of polymorphism and ubiquitous transcript expression[J]. Gene, 2014, 536(1): 207-212. doi: 10.1016/j.gene.2013.11.067
[11] HU X, GUO H, HE Y, et al. Molecular characterization of myostatin gene from Zhikong scallop Chlamys farreri (Jones et Preston 1904)[J]. Genes Genet Syst, 2010, 85(3): 207-218. doi: 10.1266/ggs.85.207
[12] KIM H W, MYKLES D L, GOETZ F W, et al. Characterization of a myostatin-like gene from the bay scallop, Argopecten irradians[J]. Biochim Biophys Acta, 2004, 1679(2): 174-179. doi: 10.1016/j.bbaexp.2004.06.005
[13] GUO L, LI L, ZHANG S D, et al. Molecular cloning and characterization of the myostatin gene in a cultivated variety of bay scallop, Argopecten irradians[J]. Aquaculture, 2012, 350(2): 192-199.
[14] MORELOS R M, RAMÍREZ J L, GARCÍA-GASCA A, et al. Expression of the myostatin gene in the adductor muscle of the Pacific lion-paw scallop Nodipecten subnodosus in association with growth and environmental conditions[J]. J Exp Zool A, 2015, 323(4): 239-255. doi: 10.1002/jez.1914
[15] NIU D, WANG L, BAI Z, et al. Identification and expression characterization of the myostatin (MSTN) gene and association analysis with growth traits in the razor clam Sinonovacula constricta[J]. Gene, 2015, 555(2): 297-304. doi: 10.1016/j.gene.2014.11.020
[16] FAN S, XU Y, LIU B, et al. Molecular characterization and expression analysis of the myostatin gene and its association with growth traits in Noble scallop (Chlamys nobilis)[J]. Comp Biochem Physiol B, 2017, 212: 24-31. doi: 10.1016/j.cbpb.2017.07.004
[17] WANG X, MENG X, SONG B, et al. SNPs in the myostatin gene of the mollusk Chlamys farreri: association with growth traits[J]. Comp Biochem Physiol B, 2010, 155(3): 327-330. doi: 10.1016/j.cbpb.2009.12.001
[18] GUO L, LI L, ZHANG S, et al. Novel polymorphisms in the myostatin gene and their association with growth traits in a variety of bay scallop, Argopecten irradians[J]. Anim Genet, 2011, 42(3): 339-340.
[19] GONZALEZ-CADAVID N F, TAYLOR W E, YARASHESKI K, et al. Organization of the human myostatin gene and expression in healthy men and HIV-infected men with muscle wasting[J]. Proc Natl Acad Sci USA, 1998, 95(25): 14938-14943. doi: 10.1073/pnas.95.25.14938
[20] JEANPLONG F, SHARMA M, SOMERS W G, et al. Genomic organization and neonatal expression of the bovine myostatin gene[J]. Mol Cell Biochem, 2001, 220(1/2): 31-37. doi: 10.1023/A:1010801511963
[21] BONGIORNI S, TILESI F, BICORGNA S, et al. Promoter polymorphisms in genes involved in porcine myogenesis influence their transcriptional activity[J]. BMC Genet, 2014, 15: 119.
[22] XUE L, DONG X, ZHANG X, et al. Organization and functional analysis of the 5' flanking regions of myostatin-1 and 2 genes from Larimichthys crocea[J]. DNA Cell Biol, 2012, 31(5): 845-855. doi: 10.1089/dna.2011.1263
[23] FUNKENSTEIN B, BALAS V, REBHAN Y, et al. Characterization and functional analysis of the 5' flanking region of Sparus aurata myostatin-1 gene[J]. Comp Biochem Physiol A, 2009, 153(1): 55-62. doi: 10.1016/j.cbpa.2008.09.031
[24] 杜荣, 安晓荣, 陈永福, 等. 绵羊Myostatin基因启动子的功能分析[J]. 中国科学C辑, 2007, 37(5): 551-557. doi: 10.3969/j.issn.1674-7232.2007.05.009 [25] SPILLER M P, KAMBADUR R, JEANPLONG F, et al. The myostatin gene is a downstream target gene of basic helix-loop-helix transcription factor MyoD[J]. Mol Cell Biol, 2002, 22(20): 7066-7082. doi: 10.1128/MCB.22.20.7066-7082.2002
[26] 夏丹. 关岭牛MSTN基因启动子相关转录因子的筛选、克隆及表达纯化[D]. 贵阳: 贵州大学, 2016: 23. [27] LI M, CHAN K, CAI D, et al. Identification and purification of an intrinsic human muscle myogenic factor that enhances muscle repair and regeneration[J]. Arch Biochem Biophys, 2000, 384(2): 263-268. doi: 10.1006/abbi.2000.2100
[28] 李佳, 邓捷, 张军林, 等. 肌肉增强子因子2对猪肌肉生长抑制素启动子活性的调节[J]. 生物工程学报, 2012, 28(8): 918-926. [29] OLSON E N, PERRY M, SCHULZ R A. Regulation of muscle differentiation by the MEF2 family of MADS box transcription factors[J]. Dev Biol, 1995, 172(1): 2-14. doi: 10.1006/dbio.1995.0002
-
期刊类型引用(1)
1. 赵德风，陈然，肖国强，滕爽爽. 泥蚶转录因子c-Myc与下游ATP结合盒转运蛋白基因ABCA3启动子的结合鉴定. 浙江农业科学. 2023(06): 1317-1322 . 
百度学术
其他类型引用(1)
计量
- 文章访问数: 4374
- HTML全文浏览量: 1947
- PDF下载量: 38
- 被引次数: 2
粤公网安备 44010502001741号
