-
鱼类骨骼肌生长发育是一个由多种转录因子参与调控的复杂过程，其中生肌调节因子（Myogenic Regulatory Factors，MRFs）发挥至关重要的调控作用[1]。MRFs基因家族主要包括MyoD、MyoG、Myf5和MRF4（Myf6）[2]，其均编码一个高度保守的碱性螺旋-环-螺旋结构域（basic helix-loop-helix，bHLH）。bHLH结构域中HelixⅠ和HelixⅡ与由E2A基因编码的E12或E47蛋白结合形成异二聚体复合物，而basic碱性区域识别E-box（一种包含CANNTG的DNA序列）中的启动子进而与其特异性结合，从而激活骨骼肌特异性的基因表达[3, 4]。因此，MRFs基因家族调控着动物体骨骼肌发育整个过程，包括肌细胞的增殖、分化及肌纤维的形成，以及个体出生后骨骼肌功能的完善[5, 6]。
Chen 等[7]在斑马鱼（Danio rerio）的研究中发现，Myf5 在准体节中胚层（presomitic mesoderm，PSM）检测到表达，同时在 PSM 中未检测到其他生肌调节因子的表达，表明 Myf5 是肌发生的关键调控因子。同样的结果也出现在鲤鱼（Cyprinus carpio）[8]、牙鲆（Paralichthys olivaceus）[9]中。研究显示，在虹鳟（Oncorhynchus mykiss）和大西洋鳕鱼（Gadus morhua）的胚胎发育过程中，MyoD 基因在体节出现到原肠胚中期均检测到表达，说明MyoD在肌肉发生过程中起着重要作用[10 − 11]。Ganassi等[12]研究表明，MyoG基因敲除可显著抑制斑马鱼胚胎期单核肌纤维生成及肌节发育，表明其在胚胎发生过程中对肌肉分化的关键调控作用。Hinits 等[13]在斑马鱼中分离到的 MRF4 基因亚型，并在终末分化的肌肉纤维中大量表达，表明 MRF4 参与肌纤维分化过程。Schnapp等[14]通过基因敲除实验发现，MRF4能够在 MyoD 和 Myf5 基因敲除的斑马鱼中激活 MyoD 促进骨骼肌的形成。目前，生肌调节因子已在部分鱼类中被克隆。例如，鳜鱼（siniperca chuatsi）和达氏鲟（Acipenser dabryanus）中发现4种MRFs因子[15 − 16]；MyoD 在牙鲆、黄尾鰤（Seriola lalandi）、金头鲷（Sparus aurata）和红鳍东方鲀（Takifugu rubripes）中被克隆[17 − 20]； MyoG分别在赤眼鳟（Squaliobarbus crriculus）、刀鲚（Coilia nasus）等中被克隆[21 − 22]；MRF4在尼罗罗非鱼（Oreochromis niloticus）中被克隆[23]。
鞍带石斑鱼（Epinephelus lanceolatu）俗称龙趸、龙胆石斑鱼，隶属鲈形目（Perciformes）、鮨科（Serranidae）、石斑鱼属（Epinephalus），最大体长可达2 m，体质量可达400 kg，为体型最大的珊瑚礁鱼类[24]。鞍带石斑鱼是我国海南省重要的海水养殖品种，由于其生长迅速的特点，常用于与其他石斑鱼如棕点石斑鱼[25]、云纹石斑鱼[26]、驼背鲈[27]等进行杂交。本研究通过克隆鞍带石斑鱼生肌调节因子，获得其基因序列，并对其氨基酸序列进行同源性对比以及进化树分析，为进一步探究生肌调节因子对骨骼肌生长发育的调控提供理论基础。
-
鞍带石斑鱼【（20.00 ± 0.07）g，3月龄】购自海南省万宁市一养殖场。在实验开始前，将鱼置于室温（28 ℃）下的循环水槽中驯化两周。随后，挑选4尾健康的鞍带石斑鱼，用间氨基苯甲酸乙酯甲磺酸（MS-222）对其进行麻醉，在无菌操作下解剖取出鞍带石斑鱼7种不同的组织，包括骨骼肌、肝脏、肾脏、肠道、胃、心脏和脂肪组织，液氮速冻后于−80 ℃冰箱中保存。
GoScript™ Reverse Transcription System试剂盒购自美国Promeg公司；FastPure® Gel DNA Extraction Mini Kit购自南京诺唯赞生物科技有限公司；DH5α感受态细胞购自上海唯地生物技术有限公司；Gold ViewI型核酸染色剂购自北京索莱宝科技有限公司；RNAiso Plus、GoScript™ Reverse Transcription System试剂盒、SMARTer® RACE5′/3′kit试剂盒、pMD19-T Vector、PrimeScript™ RT reagent kit with gDNA Eraser试剂盒、TB Green® Premix Ex Taq™ Ⅱ（Tli RNaseH Plus）均购自宝日医生物技术（北京）有限公司。
-
NCBI数据库获得斜带石斑鱼（Epinephelus coioides）、斑马鱼、鲤鱼生肌调节因子cDNA序列，根据序列保守区域设计简并引物来扩增鞍带石斑鱼MRFs基因家族中间片段。根据中间片段序列设计用于扩增3′RACE（3′-F1、3′-F2）和5′RACE（5′-R1、5′-R1）的特异性引物。将中间片段、3′RACE和5′RACE的cDNA序列拼接，设计用于验证鞍带石斑鱼生肌调节因子ORF框的特异性引物（表1），并根据克隆得到的鞍带石斑鱼5种生肌调节因子基因序列，设计用于组织分布的引物荧光定量PCR引物（表2）。引物由生工生物工程（上海）股份有限公司合成。
引物 Primer 序列 Sequence 中间片段引物序列（5′－3′）
Middle fragment primer sequences（5′－3′）MyoD1 F1 ATGCACTTCTTCGAGGAC MyoD1 R1 TGTTGTCGGTGGAGATTC MyoD2 F1 CTCTGCTGACGACCTCTA MyoD2 R1 TTACTGCTGCTGGAATCG MyoG F1 TACGACCAATCCACCTATCA MyoG R1 TCCGCCTGCTGTAGAGAT MRF4 F1 CGCTATCTGGAAGAAGGG MRF4 R1 CTCTCTGGTTTGTCATTGC Myf5 F1 ACGTCTTCTCACCATCCCAG Myf5 R1 AGAAGGTGGCCATGTCTCTC 3′RACE所用引物序列（5′－3′）
3′RACE primer sequences（5′－3′）MyoD1 3′－F1 GTCCTCCCTCCTGCACCTCCATCAC MyoD1 3′－F2 GAGCACTACAGCGGGGACTCAGACG MyoD2 3′－F1 CCCCGCTCCAACTGCTCCGATGG MyoD2 3′－F2 GCAGCCACACAGACGATTCCAGCAGC MyoG 3′－F1 CCCAACCCGACCCAACCCAGAGTGT MyoG 3′－F2 CATGCGTGCCCTGACCTCCATCGTG MRF4 3′－F1 CCAACCCCAACCAGAGGCTACCCAA MRF4 3′－F2 CGACCTCTGCCGACCATTCCACTGC Myf5 3′－F1 CGCAACGCCATCCAGTACATCGAGAGC Myf5 3′－F2 TCTGAATGCAAACTACAGCGGCGGAT 5′RACE所用引物序列（5′－3′）
5′RACE primer sequences（5′－3′）MyoD1 5′－R1 CGTGGTGATGGAGGTGCAGGAGGGA MyoD1 5′－R2
MyoD2 5′－R1ATGGAGGTGCAGGAGGGAGG
CTTACCAAGCCGCCGCCTTTCTCGCMyoD2 5′－R2
MyoG 5′－R1CTTTCTCGCATCGTCGCTGCTTTCC
CTTGGGCAGCCTCTGGTTTGGGTTCMyoG 5′－R2
MRF4 5′－R1GACGCCTCTTCTCCCTCATCGTGGC
TGGGCGCCGACTTCCTCTTGCAGATMRF4 5′－R2 CATGAGGCACTGGCCGTCGCAG Myf5 5′－R1 TGCAGGCTCTCGATGTACTGGATGGCG Myf5 5′－R2 CTGGCTGGGGTTGGCTGAGGTGC 全长验证引物序列（5′－3′）
Full-length validation primer sequences（5′－3′）MyoD1 Full－F ACGCTTCTGTCTCTGTCACTGGAT MyoD1 Full－R ATGGTAAAGTCCCAGGTGTGAAGG MyoD2 Full－F CCCCGCTCCAACTGCTCCGATGG MyoD2 Full－R GCAGCCACACAGACGATTCCAGCAGC MyoG Full－F CCGAACACCACATACAGTAGAGC MyoG Full－R TGGAGGGTCTACTTGGGAAT MRF4 Full－F GGGCGCAACAATATGATGGACCTTT MRF4 Full－R GAATATAGTGGAAATAATTTAGGGA Myf5 Full－F ATGGAYGTCTTYTCVCCVTCCCA Myf5 Full－R TGGCGTTCGGTAAAAGCAG 接头引物序列（5′-3′）
Adapter primer sequences（5′－3′）3′ CDS AAGCAGTGGTATCAACGCAGAGTACTT
TTTTTTTTTTTTTTTTTTTTTTTTTTTTVN5′ CDS TTTTTTTTTTTTTTTTTTTTTTTTTVN NUPA AAGCAGTGGTATCAACGCAGAGT UPMA－S CTAATACGACTCACTATAGGGC UPMA－L CTAATACGACTCACTATAGGGCAAGCAG
TGGTATCAACGCAGAGTTable 1. Primer sequences for gene cloning
引物 Primer 序列 Sequence MyoD1 RT－F GGCATGACGGATTTCAACGG MyoD1 RT－R CTGCTGTTGTCGGTGGAGAT MyoD2 RT－F ATCATCACCACCACCACGTC MyoD2 RT－R TTCAGGGTCTCAAAGGCGTC MyoG RT－F CACGATGAGGGAGAAGAGGC MyoG RT－R TGCTGGTTGAGGGAAGACAC MRF4 RT－F CAGTGTGGCCAGTCATGAGT MRF4 RT－R TCTTCGGTGGTGATGCTGTC Myf5 RT－F GAGGAGGACGAGCATGTGAG Myf5 RT－R CTCGATGTACTGGATGGCGT EF1-α RT－F AGGGATGGAAGATTGAGCGC EF1-α RT－R CGTACCGGGCTTCAGGATAC Table 2. Nucleotide sequences of the primers used for Real-time PCR
-
使用Trizol试剂（Invitrogen，美国）提取总RNA，1.0%的琼脂糖凝胶电泳检测RNA质量，并在NanoDrop ND-
1000 分光光度计（NanoDrop，威尔明顿，特拉华州）上进行产率测定。RNA经脱氧核糖核酸去除酶（Takara，日本）去除DNA后，中间片段cDNA模板制备和3′端cDNA模板制备通过GoScript™ Reverse Transcription System试剂盒（Promega，美国）进行反转录；5′端cDNA模板制备采用SMARTer® RACE5′/3′kit（Takara，日本）试剂盒反转录合成；不同组织分布cDNA模板制备采用PrimeScript™ RT reagent kit with gDNA Eraser（Perfect Real Time）（Takara，日本）试剂盒进行反转录合成。 -
以所获得中间片段cDNA为模板，以引物 F1 和 R1 进行第1轮 PCR 扩增，PCR反应条件如下：94 ℃ 5 min；94 ℃ 30 s，TM值30 s，72 ℃ 60 s，35循环；72 ℃ 5 min；4 ℃ 保持。
-
以3' 端cDNA为模板，分别以鞍带石斑鱼MRFs基因家族基因正义引物3'-F1和反义接头引物UPMA（UPMA-L:UPMA-S为1﹕5混合使用）进行第1轮PCR扩增；将第1轮PCR产物稀释100倍，并取2 μl作为第2轮PCR扩增模板，以正义引物3'-F2和反义接头引物NUPA进行第2轮嵌套PCR扩增，第1轮PCR反应条件：94 ℃ 5 min；94 ℃ 30 s，TM值+5 ℃（每循环降落0.5 ℃） 30 s，72 ℃ 60 s，20循环；94 ℃ 30 s，TM值 30 s，72 ℃ 60 s，20循环；72 ℃ 7 min，4 ℃ 保持。第2轮PCR：94 ℃ 5 min；94 ℃ 30 s，TM值30 s，72 ℃ 60 s，40循环；72 ℃ 10 min，4 ℃ 保持。
-
以5'端cDNA为模板，分别以鞍带石斑鱼MRFs基因家族基因正义引物5'-R1和反义接头引物UPMA进行第1轮PCR扩增；将第1轮PCR产物稀释10倍，并分别取2 μL作为第2轮PCR扩增模板，分别以鞍带石斑鱼MRFs基因家族的正义引物3'-R2和反义接头引物NUPA进行第2轮嵌套PCR扩增，PCR反应条件同1.4.2。
-
PCR产物用1%琼脂糖凝胶电泳检测，并根据FastPure® Gel DNA Extraction Mini Kit（诺唯赞，中国）说明书切取含有目的片段凝胶，胶回收纯化PCR产物。
-
回收产物连接于pMD19-T Vector，转化到大肠杆菌（DH5α）感受态细胞中，涂在含氨苄青霉素LA平板上，挑取白色单菌斑培养，进行PCR阳性克隆检测，分别取100 μL菌液送至广州天一辉远公司测序。
-
根据测序结果，将MRFs基因家族cDNA序列拼接以获得基因全长序列，利用在线ORF Finder对5个基因ORF进行预测，并在其序列两端非编码区设计特异性引物（MyoD1 Full-F，MyoD1 Full-R；MyoD2 Full-F，MyoD2 Full-R；MyoG Full-F，MyoG Full-R；MRF4 Full-F，MRF4 Full-R；Myf5 Full-F，Myf5 Full-R）进行PCR全长扩增反应，以进行全长验证。
-
鞍带石斑鱼MyoD1、MyoD2、MyoG、MRF4、Myf5 5个基因cDNA序列使用ORF Finder（
http://www.bioinformatics.org/sms2/orf_find.html ）预测基因编码区；核苷酸和氨基酸序列同源性比对使用NCBI（https://www.ncbi.nlm.nih.gov/ ）blast程序进行；核酸序列翻译成基酸序列通过ExPASy（https://web.expasy.org/translate/ ）完成；采用SWISS-MODEL（https://swissmodel.expasy.org/ ）对蛋白质三级结构进行预测分析；多重氨基酸序列对比图谱使用 Clustal X和 ESPript 3.0构建；在NCBI数据上进行比对，获得MRFs基因同源序列。利用MEGA7.0软件中的ClustalW算法对MRFs基因的氨基酸序列进行多重序列比对，并采用Neighbor-Joining（NJ）法构建系统进化树，通过自引导检验（bootstrap）获得系统分支的置信度（重复次数为1 000）。 -
鞍带石斑鱼骨骼肌、脑、心脏、脂肪、胃等组织cDNA进行实时荧光定量 PCR 反应，根据TB Green® Premix Ex TagTM II（Takara，日本）试剂盒说明书进行操作。实时荧光定量在定量热循环仪（Roche Lightcycler® 480,瑞士）中进行。扩增体系为20 μL，其中包含10 μL 2×Taq PCR MasterMix II反应混合液，0.5 μL的正反引物，7 μL无核酶水和2 μLcDNA模板，反应程序为95 ℃ 30 s预变性，95 ℃ 5 s，55 ℃ 30 s，72 ℃ 30 s，共40循环，95 ℃ 5 s，65 ℃ 60 s绘制溶解曲线; 50 ℃ 30 s冷却。以EF1-α 作为内参基因。用5种浓度稀释不同倍数的cDNA样品制作标准曲线（每个cDNA样品设置3次重复试验），根据公式E=10（−1/slope）−1分析扩增效率。目的基因的相对表达水平采用R=2−ΔΔCt法进行计算[28]。
-
本研究将克隆得到的基因序列上传已上传至 GenBank，登录号分别为MW762676、MW762677、MW762678、MW762679、MW762680。
-
实验数据均以平均值±标准差表示，并进行了方差齐性检验和正态性检验。使用 SPSS22.0 软件对数据进行单因素方差分析（one-way ANOVA），使用 Duncan’s 进行多重比较分析，以 P < 0.05作为差异显著性判断标准。
-
克隆获得鞍带石斑鱼MyoD1、MyoD2、MyoG、MRF4和Myf5 cDNA序列长度分别为1 962 、1 466 、926 、1 052和1 133 bp，开放阅读框（ORF）分别为891、810、753、714 和726 bp，分别编码297、270、251、238和242个氨基酸。GenBank登录号分别为MW762676、MW762677、MW762678、MW762679和MW762680。结构分析表明，鞍带石斑鱼MRFs家族5个基因均由编码区（CDS）和非翻译区（UTR）组成（图1）。采用Garnier-Gibrat-Robson方法预测其2级结构，结果显示MyoD1、MyoD2、MyoG、MRF4和Myf5均包含多个α-螺旋和β-折叠结构域（图2）。3级结构预测进一步表明其均呈现典型的螺旋-环-螺旋（HLH）结构（图3）。
Figure 1. Gene structure of the Epinephelus lanceolatus MRFs family genes
Figure 2. Prediction of the secondary structure of MRFs protein in Epinephelus lanceolatus .
Figure 3. Tertiary structure of MRFs family proteins in Epinephelus lanceolatus
-
采用DNAMAN软件对鞍带石斑鱼5种生肌调节因子氨基酸序列进行比对，结果表明其具有60个左右氨基酸组成的保守碱性螺旋-环-螺旋（basic helix-loop-helix，bHLH）结构域（图4）。
Figure 4. Multiple alignments of the five MRFs amino acid sequences.
-
多重序列比对结果表明，鞍带石斑鱼MyoD1与金头鲷同源性最高，为92.93%，其次是草鱼（68.54%），与哺乳类同源性在54.49%～56.12%，与鸟类为64.47%（图5）
Figure 5. Multiple sequence alignment of the MyoD1 amino acid sequences between Epinephelus lanceolatus and other species
-
多重序列比对结果表明，鞍带石斑鱼MyoD2与斜带石斑鱼同源性最高（98.51%），其次是金头鲷为85.14%，与红鳍东方鲀（70.00%）最低（图6）。
Figure 6. Multiple sequence alignment of the MyoD2 amino acid sequences between Epinephelus lanceolatus andf other species
-
多重序列比对结果表明，鞍带石斑鱼MyoG与斜带石斑鱼和金头鲷同源性最高，同源性分别为98.80%、94.40%。与哺乳类（50.20%～51.39%）和鸟类（56.18%）同源性较低（图7）。
Figure 7. Multiple sequence alignment of the MyoG amino acid sequences between Epinephelus lanceolatus and other species
-
多重序列比对结果表明，鞍带石斑鱼MRF4与斜带石斑鱼同源性最高，为97.47%，其次是金头鲷（ 82.92%）和红鳍东方鲀（81.40%）。与哺乳类（50.20%～51.39%）和鸟类（56.18%）同源性最低（图8）。
Figure 8. Multiple sequence alignment of the MRF4 amino acid sequences between Epinephelus lanceolatus and other species
-
多重序列比对结果表明，鞍带石斑鱼Myf5与斜带石斑鱼同源性最高，为99.59%，其次是金头鲷（91.29%），与红鳍东方鲀、虹鳟及草鱼同源性分别为87.97%、81.82%和75.21%，与哺乳类（56.20%～56.59%）和鸟类（57.09%）同源性较低（图9）。
Figure 9. Multiple sequence alignment of the Myf5 amino acid sequences between Epinephelus lanceolatus and other species
-
采用MEGA7.0软件基于Neighbor-Joining（N-J）法构建了包括鞍带石斑鱼、人、小鼠、红原鸡、草鱼、斑马鱼等11种物种的系统进化树。结果显示，5种生肌调节因子（MRFs）明显分为两支，其中MyoD1、MyoD2与Myf5聚为一支，MRF4与MyoG聚为另一支。在系统进化过程中，每个基因均明显划分为鱼类、哺乳类和鸟类3大分支。在鱼类分支中，鞍带石斑鱼MRFs与斜带石斑鱼亲缘关系最近，二者聚成簇后与金头鲷聚在一起，再与红鳍东方鲀形成一小支，且各节点具有较高的自展值支持；与草鱼、斑马鱼、虹鳟的进化关系较远（图10）。
Figure 10. Phylogenetic tree based on MRFs amino acid sequences of different species（Neighbor-Joining method）
-
如图11-A所示，MyoD1基因在肠道、胃、肝脏和骨骼肌中均有表达，其中骨骼肌表达量最高，且尾肌表达量略高于肋肌，但两者无显著性差异（P > 0.05）；在心脏和脂肪组织中几乎检测不到MyoD1表达。MyoD2基因在以上6种组织中均有表达，骨骼肌中表达量最高，且尾肌表达量显著高于肋肌（P < 0.05），心脏中表达量最低（图11-B）。MyoG基因在6种组织中均有表达，骨骼肌中表达量最高，肋肌表达量略高于尾肌，但两者无显著性差异（P > 0.05）（图11-C）。MRF4基因在肠道、胃、心脏、骨骼肌和脂肪中均有表达，骨骼肌中表达量最高，且肋肌表达量显著高于尾肌（P < 0.05），而在肝脏中几乎检测不到MRF4表达（图11-D）。Myf5基因在除肝脏外的5种组织中均有表达，骨骼肌中表达量最高，且肋肌与尾肌表达差异与MyoD2结果一致（图11-E）。总体而言，5种生肌调节因子mRNA在骨骼肌表达量均显著高于其他组织（P < 0.000 1）。
Figure 11. Relative expression levels of 5 MRFs in different tissues of Epinephelus lanceolatus
-
本研究成功克隆了鞍带石斑鱼MyoD1、MyoD2、MyoG、MRF4和Myf5 cDNA序列，其长度分别为1 962、1 466 、926 、1 052 和
1133 bp，分别编码297、270、251、238和242个氨基酸。蛋白结构预测分析表明，鞍带石斑鱼5种MRFs均具有高度保守的碱性螺旋-环-螺旋结构域（bHLH），该区域由60个左右的氨基酸残基组成，其中碱性区域（basic）识别并结合E-box启动子，而HLH区域与E蛋白形成异二聚体复合物，从而激活肌肉特异性基因转录[3]。蛋白结构预测还显示，Myf5、MyoD1和MyoD2均属于Myf5基因超家族，且Myf5与MyoD在功能上具有相似性，均参与肌卫星细胞增殖及肌纤维增生调控。Kassar-Duchossoy等[29]研究发现，Myf5或MyoD单基因敲除小鼠的骨骼肌发育基本正常，但双基因敲除小鼠则完全丧失成肌能力，表明二者在肌生成中具有功能互补性。Relaix等[30]证实，Pax3和Pax7通过结合Myf5和MyoD启动子中的保守E-box和Pax结合位点，直接激活其表达。MyoD作为“主调控因子”的功能依赖于其启动子中的E-box元件（CANNTG），这一特征在Myf5中同样保守。Tapscott等[31]首次提出，E-box是肌源性分化的核心调控模块，且Myf5与MyoD共享这一元件的结合特性。这些研究表明，Myf5与MyoD在肌卫星细胞增殖中的冗余功能及调控机制的保守性，进一步支持了二者由同一个祖先基因复制进化而来的假说。同源性分析结果表明，鞍带石斑鱼MRFs氨基酸序列与其他鱼类的同源性较高，为68.54%～99.59%。特别是与鲈形目的同源性最高，为82.92%～99.59%；而与鸟类和哺乳类的同源性相对较低，为50.20%-64.47%。证明本研究所得，鞍带石斑鱼MRFs序列所编码的蛋白质与各物种的基因编码的蛋白质具有相似的功能，在鱼类中的高保守性同时也揭示了该蛋白质功能的高度保守性。结合系统进化树结果，鞍带石斑鱼与同属石斑鱼属的斜带石斑鱼聚在同一分支，并与鲈形目鱼类形成一簇，说明亲缘关系最近；而与哺乳类、鸟类的距离较远，说明与陆生动物亲缘性较远。与陈敦学基于线粒体构建的进化树结果相符[32]。鱼类MRFs与鸟类和哺乳动物的同源性存在较大差异，推测与物种进化、亲缘关系以及所含的氨基酸残基数目不同有关。
根据结构和功能的不同，MRFs基因家族可分为初级MRFs（Myf5和MyoD）和次级MRFs（MyoG和Myf6）。初级MRFs同源性较高，功能互补，主要调控体节细胞向生肌细胞系的定向分化；次级MRFs则促进成肌细胞增殖并分化为肌细胞，最终形成肌纤维[33]。在鱼类胚胎发育中，MyoD和Myf5表达通常早于MyoG和MRF4。例如，在牙鲆胚胎发育过程中，Myf5和MyoD在受精26 h后即可检测到，而MyoG的表达则首次出现在受精32.5 h后[34]。类似结果在鲤鱼[8]和虹鳟[35]中也有报道。然而，不同鱼类间MRFs表达模式也存在差异。例如，牙鲆MyoG在胚胎发育形成3个体节时首次表达，其表达水平持续上调至30个体节后逐渐下调[34]；而大西洋鲱（Clupea harengus）MyoG在胚胎发育形成34个体节前仅短暂表达，并在MyoD开始表达前从体节中消失[36]。这些研究表明，MRFs基因表达时序在鱼类中具有保守性，但也存在物种特异性差异。而进化分析结果同样显示，鞍带石斑鱼MRFs 基因家族明显分为两支，其中MyoD与Myf5聚成一支，而MyoG与MRF4聚成另一支，说明MyoD与Myf5同源性较高，可能与其在肌肉生长发育调控中的功能一致。
本研究发现，鞍带石斑鱼MRFs基因在骨骼肌中呈现高表达，这与在其他鱼类中的研究结果一致。例如，MyoD在斑马鱼胚胎体节及成体骨骼肌卫星细胞中特异性高表达，其敲除可导致体节肌细胞分化完全停滞，证实MyoD是斑马鱼骨骼肌发育的核心调控因子[37]。在鲶鱼（Silurus asotus）中，MyoG在骨骼肌快肌纤维分化阶段特异性高表达，通过激活肌球蛋白重链（MyHC）及肌钙蛋白编码基因，直接驱动肌纤维增粗[38]。松江鲈（Trachidermus fasciatus）的MRF4在骨骼肌中高表达，MRF4不仅维持成熟肌纤维功能稳态，还在肌肉再生过程中被激活，参与修复调控[39]。此外，虹鳟的Myf5主要在骨骼肌中表达，Myf5在肌肉前体细胞分化和胚胎肌肉发育中发挥主导作用[40]。这些结果表明，MRFs家族基因在鱼类骨骼肌中的高表达与其在肌肉生成和分化中的核心功能密切相关。此外，鞍带石斑鱼MRFs基因在其他组织中也检测到微量表达，推测其可能不仅参与肌肉分化，还在其他组织器官中发挥潜在作用。类似现象在其他鱼类中也有报道，如虹鳟Myf5在肝脏、肠、心脏、鳃和脑中有少量表达[40]，达氏鲟的MRFs基因在心脏、脑、肝、脾和肾中也有表达[16]。这些发现进一步支持了MRFs基因在鱼类多种组织中的多功能性。研究表明，鱼类早期生长发育的各个阶段均受到MRFs基因表达的影响。例如，在牙鲆中，MyoD在胚胎发生的初始阶段即在前体肌肉细胞中检测到表达[19]；在鳜鱼中，Myf5基因在原肠胚期这一早期发育阶段即开始表达[41]；此外，Johansen等[42]研究发现，虹鳟在不同生长阶段中，骨骼肌内的MyoD、Myf5、MyoG和MRF4表达水平存在显著差异。基于上述研究结果可以推断，MRFs在不同发育阶段的差异性表达可能与其调控肌肉生长的分子机制密切相关。
综上所述，本研究成功克隆得到鞍带石斑鱼5种生肌调节因子MyoD（MyoD1和MyoD2）、MyoG、MRF4、Myf5 cDNA 全长序列，并通过生物信息学软件对5种生肌调节因子 CDS区、氨基酸序列、蛋白结构、进化地位、亲缘关系以及不同组织中相对表达等进行了分析。鞍带石斑鱼5种生肌调节因子都具有高度保守的碱性螺旋-环-螺旋结构域（basic helix-loop-helix，bHLH），同时每个基因均与鲈形目鱼类分为一支，且5种生肌调节因子皆在骨骼肌中表达丰度较高，而在肝脏、心脏、肠道等其他组织中表达量较低。研究结果为进一步探究生肌调节因子在鞍带石斑鱼骨骼肌生长发育的作用机理奠定了基础。
Gene cloning and expression analysis of MRFs gene cDNA in the tissues of giant grouper Epinephelus lanceolatus
doi: 10.15886/j.cnki.rdswxb.20250030
- Received Date: 2025-02-12
- Accepted Date: 2025-03-21
- Rev Recd Date: 2025-03-14
-
Key words:
- Epinephelus lanceolatus /
- myogenic regulatory factors /
- cDNA cloning /
- gene expression
Abstract: Myogenic regulatory factors (MRFs) are critical regulators of skeletal muscle growth and development in animals. The full-length cDNA sequences of five MRFs from the giant grouper (Epinephelus lanceolatus) were cloned using RT-PCR and RAC, including MyoD1 (1,962 bp), MyoD2 (1,466 bp), MyoG (926 bp), MRF4 (1,052 bp), and Myf5 (1,133 bp), encoding 297, 270, 251, 238, and 242 amino acids, respectively. Amino acid sequence alignment revealed that all five MRFs contain a conserved basic helix-loop-helix (bHLH) domain consisting of 60 amino acids. Phylogenetic analysis demonstrated that the MRF genes are clustered into two distinct branches: MyoD (MyoD1 and MyoD2) groups with Myf5 in one clade, while MRF4 and MyoG form the other. All genes showed closest evolutionary relationships with Perciformes homologs. Real-time quantitative PCR analysis revealed predominant expression of these MRFs in skeletal muscle, with significantly lower expression levels in the liver, heart, and intestine. These findings provide foundation for elucidating the regulatory mechanisms of MRFs in skeletal muscle development of the giant grouper.
| Citation: | ZHANG Longting, ZHENG Zhi, WANG Yao, GAO Yujie. Gene cloning and expression analysis of MRFs gene cDNA in the tissues of giant grouper Epinephelus lanceolatus[J]. Journal of Tropical Biology. doi: 10.15886/j.cnki.rdswxb.20250030
|
DownLoad:
