畜牧兽医学报  2021, Vol. 52 Issue (9): 2439-2451. DOI: 10.11843/j.issn.0366-6964.2021.09.007    PDF    
5-Aza-dC对牛成肌细胞MyoD1启动子甲基化及mRNA表达的影响
李雄1,2, 田念念1,2, 宋林锦1,2, 陈晨1,2, 许厚强1,2     
1. 贵州大学 高原山地动物遗传育种与繁殖教育部重点实验室/贵州省动物遗传育种与繁殖重点实验室, 贵阳 550025;
2. 贵州大学动物科学学院, 贵阳 550025
摘要:旨在探究5-氮杂-2-脱氧胞苷(5-Aza-2’-deoxycytidine,5-Aza-dC)对牛成肌细胞的增殖和肌分化因子(myogenic differentiation 1,MyoD1)启动子甲基化以及mRNA表达的影响。本研究复苏了前期冻存的关岭牛成肌细胞,并进行培养,待其生长到对数生长期,再通过不同浓度的5-Aza-dC对牛成肌细胞进行处理,利用流式细胞仪检测细胞凋亡和周期;结合高通量检测方法检测MyoD1启动子甲基化水平,并利用qRT-PCR检测MyoD1及相关基因的表达水平。研究发现,0.1 μmol·L-15-Aza-dC为最适浓度,该浓度下细胞抗凋亡因子Bcl的表达极显著降低(P < 0.01),促凋亡因子Bax的表达极显著提高(P < 0.01),促凋亡因子Caspase-9的表达显著提高(P < 0.05);周期因子Cyclin A2的表达显著提高(P < 0.05),而细胞因子Cyclin B1、Cyclin D的表达无显著变化;检测空白组和试验组MyoD1启动子甲基化水平发现,试验组甲基化水平极显著低于空白组(P < 0.01);而mRNA的表达水平极显著高于空白组(P < 0.01)。5-Aza-dC能够通过改变Caspase-9、Bcl、Bax、Cyclin A2等基因的表达来调节关岭牛成肌细胞的增殖和凋亡,并能够有效降低MyoD1基因启动子的甲基化水平(P < 0.01);极显著提高其mRNA的表达量(P < 0.01)。低浓度的5-Aza-dC能通过促进细胞凋亡及调控关岭牛MyoD1的甲基化水平来调控MyoD1的表达;同时推测,MyoD1基因启动子的甲基化水平能够影响关岭牛成肌细胞的生长发育,可为遗传标记辅助关岭牛的改良提供理论参考。
关键词关岭牛    5-氮杂-2-脱氧胞苷    DNA甲基化    细胞周期    细胞凋亡    
Effects of 5-Aza-dC on MyoD1 Promoter Methylation and mRNA Expression in Bovine Myoblasts
LI Xiong1,2, TIAN Niannian1,2, SONG Linjin1,2, CHEN Chen1,2, XU Houqiang1,2     
1. Key Laboratory of Animal Genetics, Breeding and Reproduction in Guizhou/Key Laboratory of Animal Genetics, Breeding and Reproduction in the Plateau Mountains Region of Ministry of Education, Guizhou University, Guiyang 550025, China;
2. College of Animal Sciences, Guizhou University, Guiyang 550025, China
Abstract: The aim of this study was to investigate the effects of 5-Aza-2'-deoxycytidine (5-Aza-dC) on the proliferation of bovine myoblasts, promoter methylation and mRNA expression of myogenic differentiation 1 (MyoD1). In this study, the Guanling cattle myoblasts cryopreserved in the early stage were recovered and cultured. After they had grown to the logarithmic growth phase, the bovine myoblasts were treated by different concentrations of 5-Aza-dC. The flow cytometry was used to detect cell apoptosis and cycle. The methylation level of MyoD1 promoter was detected by high-throughput detection methods, the expression level of MyoD1 and related genes were detected by qRT-PCR. The results showed that 0.1 μmol·L-1 5-Aza-dC was the optimal concentration, and the expression of anti-apoptotic factor Bcl was significantly decreased (P < 0.01), the expression of pro-apoptotic factor Bax was significantly increased (P < 0.01), and the expression of pro-apoptotic factor Caspase-9 was significantly increased (P < 0.05); The expression of Cyclin A2 was significantly increased (P < 0.05), but the expression of Cyclin B1 and Cyclin D were not significantly affected. The MyoD1 promoter methylation level in blank group and experimental group was detected, and the MyoD1 promoter methylation level in experimental group was significantly lower than that in blank group (P < 0.01), while the mRNA expression level was significantly higher in experimental group than that in blank group (P < 0.01). In summary, 5-Aza-dC can regulate the proliferation and apoptosis of Guanling cattle myoblasts by changing the expressions of Caspase-9, Bcl, Bax, Cyclin A2 and other genes, and can effectively reduce the methylation level of MyoD1 gene promoter (P < 0.01), and significantly increase its mRNA expression (P < 0.01). Low concentration of 5-Aza-dC can promote cell apoptosis and regulate MyoD1 methylation level to regulate the expression of MyoD1 in Guanling cattle. At the same time, it is speculated that the methylation level of the MyoD1 gene promoter can affect the growth and development of Guanling cattle myoblasts, which can provide a theoretical reference for screening genetic markers to assist the improvement of Guanling cattle.
Key words: Guanling cattle    5-Aza-2'-deoxycytidine    DNA methylation    cell cycle    apoptosis    

表观遗传调控多发生在相互作用水平上,在该水平上的研究可深入了解表观遗传学对生物学过程的促进作用。大量研究表明,表观遗传学是指基因的核苷酸序列不发生改变,而基因功能发生稳定的、可遗传的变化,主要包括非编码RNA调控[1]、DNA甲基化[2]、染色质重塑[3]、组蛋白修饰[4]等。DNA甲基化(DNA methylation)是指在DNA甲基化转移酶的作用下,活性甲基化合物(如S-腺苷基甲硫氨酸)将甲基催化转移到基因组CpG二核苷酸的胞嘧啶5’碳位共价键结合一个甲基基团,导致染色体结构、DNA构象、DNA稳定性等发生变化,从而控制基因表达[5]。对DNA甲基化研究主要集中于主动甲基化和去甲基化方面,DNA高甲基化能够减弱或者沉默基因的表达,而去甲基化时则能使基因活化或重新表达[6]。5-Aza-dC是目前用得较多的去甲基化方法之一,其对研究DNA甲基化对表观遗传的影响及特异性基因的表达和调控机理有着重要的作用。5-Aza-dC是一种胞唾陡类似物,常用名为地西他滨甲基化酶抑制剂,是目前已知的较为有效的一种去甲基化制剂[7]。5-Aza-dC大多应用于医学方面,如Greville等[8]发现,5-Aza-dC处理后,从化疗敏感/非转移细胞系中分泌的n-聚糖的分支和唾液化增加。这些变化与卵巢癌细胞系中MGAT5和ST3GAL4转录本mRNA表达水平升高有关。Xie等[9]研究发现,5-Aza-dC能有效抑制DNMTs的表达和活性,逆转MnSODGSTT1的表达。所以,通过药物、试剂等对相关基因去甲基化探究其表达影响,也是从侧面研究甲基化机制的一种方法。

近年来,大量研究表明,DNA甲基化所引起的表观遗传对于动物的生长发育有重要影响,并筛选与生长相关的甲基化位点作为分子标记[10]。肌肉组织是动物机体的重要组成部分,而DNA甲基化也能调控肌肉的发育过程[11]。Li等[12]研究了3个品种猪不同部位骨骼肌,发现品种间、性别间存在DNA甲基化的相似性和特异性,并鉴定了差异甲基化区域。李秀金等[13]通过对五指山和长白猪DNA甲基化和转录组数据的联合分析证实,DNA去甲基化酶Tet1在猪胚胎成肌细胞和C2C12细胞中通过Myogenin启动子去甲基化,上调Myogenin表达,从而促进成肌细胞的分化。万火福等[14]用比色法测定不同品种肌肉组织DNA甲基化水平,发现DNA甲基化水平在不同品种及不同组织中会呈现一种动态变化,且高水平的DNA甲基化可能促进霞烟鸡腿肌组织的发育。综上表明,不同畜禽肌肉发育过程中的甲基化模式具有一定特异性,不同家畜中甲基化状态会产生特定变化,这种甲基化状态变化能够通过调控不同阶段相关基因的表达,从而保证肌肉发育正常进行。

动物肌肉发育是一个极其复杂的过程,动物体内相关基因的甲基化状态也呈现出动态变化趋势[15]。肌分化因子(myogenic differentiation 1,MyoD1)是一种骨骼肌特异性转录因子,能够识别、启动骨骼肌特异基因的转录,具有诱导成纤维细胞和其他非骨骼肌细胞分化为骨骼肌细胞的能力[16]。近年来,MyoD1甲基化在猪[17]、羊[18]、鸡[19]等家畜上主要通过对肌肉进行全基因组甲基化检测其甲基化水平;沈兰等[20]发现,MyoD1的甲基化谱可能有助于预测接受标准化学放射治疗的浸润性宫颈癌患者的反应。目前,MyoD1甲基化对牛骨骼肌细胞发育影响的研究报道较少。本试验以关岭牛为研究对象,利用甲基化酶抑制剂(5-Aza-dC)处理牛成肌细胞,研究MyoD1启动子区甲基化水平及其对成肌细胞生长发育和增殖的影响,为关岭牛的改良提供理论参考。

1 材料与方法 1.1 试验材料

牛成肌细胞来源于贵州大学动物科学学院动物种质资源库。

1.2 主要试剂

DMEM/F12培养基、胎牛血清、TRIzol试剂购于Gibco公司;5-Aza-dC (5-氮杂-2′-脱氧胞苷)、细胞凋亡检测试剂盒、细胞周期试剂盒均购自美国Sigma公司;2×Es Taq Master Mix、DM 2000 Marker、丙三醇、无水乙醇、DEPC水购自康为世纪生物科技有限公司;SYBR-Green染料(SYBR Green PCR Master Mix)购自Bio-Rad公司;二甲基亚砜购于Hyclone公司;逆转录试剂盒和CCK8试剂盒购自MCE公司;一抗(Desmin)购自鼎国生物技术有限责任公司;二抗购自Jackson immuno research。

1.3 试验方法

1.3.1 细胞复苏及鉴定   取出冻存的细胞,复苏后转入含15%胎牛血清的DMEM完全培养基中。置于CO2培养箱(培养条件5% CO2、饱和湿度、37 ℃) 中继续培养,让其长到对数生长期。将培养的细胞接种于6孔细胞培养板上,待其生长至85%左右,进行细胞鉴定;试验参考间接免疫荧光(indirect immunofluorescence,IIF)的方法执行。

1.3.2 CCK8检测不同浓度5-Aza-dC对成肌细胞增殖凋亡的影响   细胞复苏后,当细胞长到85%以上进行分瓶,当关岭牛成肌细胞处于对数生长期时,用0.25%胰蛋白酶消化约4 min,并进行计数,以每孔5×103个接种于96孔板内(每孔加入100 μL)。置于CO2培养箱内培养12 h;12 h后加入添加了不同浓度5-Aza-dC的培养基,浓度分别为0.05、0.1、0.2、0.3 μmol·L-1以及空白组(0 μmol·L-1),每组设3个重复,分别处理0、24和48 h。在每个孔内加入10 μL CKK-8试剂,酶标仪(450 nm)检测吸光密度(OD);确定在本实验室条件下对细胞生长影响最小的5-Aza-dC最适浓度,并以此浓度进行后续试验。

1.4 流式细胞仪检测细胞凋亡

将对数生长期细胞接种于6孔板内,标记处理组和空白组置于CO2培养箱内培养12 h,用含有5-Aza-dC的DMEM/F12(“1.3.2”获得5-Aza-dC最适浓度)培养基换液,培养48 h,用3 mL PBS溶液轻轻润洗培养板内细胞;加入1 mL 0.25%不含EDTA的胰酶消化;将细胞轻轻重悬于所保留的培养基中,密度大约为1×106个·mL-1,1 000 r·min-1离心6 min,弃去上清液,PBS洗涤2次,加入冰预冷的70%的乙醇固定,4 ℃,3 h;离心弃去固定液,3.5 mL PBS重悬5 min;将细胞悬液转移到干净的离心管内,加入2 μL PI,5 μL Annexin V-FITC,室温避光反应15 min;室温1 000 r·min-1离心5 min,去除上清;将细胞用0.5 mL Propidium Iodide预冷缓冲液,用流式细胞仪检测分析。

1.5 流式细胞仪检测细胞周期

5-Aza-dC处理细胞48 h后,用预冷PBS洗2次,加入1 mL 0.25%不含EDTA的胰酶消化,1 300 r·min-1离心6 min,获取细胞沉淀,弃上清,重悬细胞,迅速打入3.5 mL 70%预冷的乙醇中,吹打均匀,4 ℃储存过夜;离心乙醇固定过的细胞,弃上清;洗涤细胞3次以去除残留的乙醇;加入100 μL RNase A重悬细胞,37 ℃水浴30 min,加入500 μL PI染色液,4 ℃避光反应30 min;应用流式细胞仪检测。

1.6 qRT-PCR检测MyoD1、DNMT1及相关基因mRNA表达

1.6.1 RNA提取及第一条cDNA的合成   5-Aza-dC处理细胞48 h后,处理组和空白组的细胞用Trizol提取总RNA。将RNA定量至2 ng,根据逆转录试剂盒进行反应合成cDNA,反应产物置于-20 ℃保存备用。

1.6.2 引物的设计与合成   根据NCBI公布的牛MyoD1(NM_001040478.2)、DNMT1(NM_182651.2)、Caspase-9 (NM_001205504.2)、Bax(NM_173894.1)、Bcl(CX950904.1)、Cyclin A2(NM_001075123.1)、Cyclin B1(BC116011.1)、Cyclin D(XM_024991034.1)以及内参基因GAPDH(NM_001034034.2)。利用Primer Premier 5.0软件进行定量引物设计,最后由重庆擎科生物有限公司合成,引物序列见表 1

表 1 引物序列信息 Table 1 Primer sequence information

1.6.3 qRT-PCR检测MyoD1、DNMT1及相关基因表达   以处理组和空白组cDNA为模板分别扩增MyoD1以及相关基因,采用SYBR Green荧光染料法在实时荧光定量PCR仪上进行定量反应。根据伯乐公司Bio-rad iTaq Univer SYBR (1725121)荧光定量qPCR mix试剂盒制作上机PCR混池,混池体系为10 μL:ULtraSYBR Mixture(High ROX) 5 μL,上、下游引物各1 μL,cDNA模板0.7 μL,ddH2O 2.3 μL。反应条件:95 ℃ 1 min;95 ℃ 10 s;58 ℃ 30 s,40个循环后,进行熔解曲线分析;以每5 s上升0.5 ℃的速率从65 ℃升高到95 ℃,荧光信号在循环结束时检测,每个样品做3个生物学重复。

1.7 高通量检测MyoD1启动子区甲基化

依据原始提供的牛MyoD1启动子区区域进行引物设计(表 2),选取检测位置区域(NC_037342.1)将样品送往重庆擎科生物科技有限公司进行高通量测序。

表 2 甲基化测序引物序列 Table 2 Primer sequences for methylation sequencing
1.8 统计分析

运用Excel 2010软件对各组织中目的基因和内参基因mRNA表达量进行处理,运用2-ΔΔCt法计算。采用独立样本t检验和单因素方差分析进行数据处理(SPSS 21.0),多重比较采用Duncan’s多重极差检验法,运用Graphpad6作图,P < 0.05表示差异显著,P < 0.01表示差异极显著。

2 结果 2.1 间接免疫荧光细胞鉴定

图 1所示,Desmin抗体发生特异性结合,呈阳性反应显红色,细胞核在DAPI的染色下呈蓝色,Merge后能完全重合,表明培养的细胞为成肌细胞。

A. DAPI核染;B. Desmin染色;C. 合图Merge;D. DAPI核染对照;E.不加Desmin一抗对照染色 A. DAPI nuclear staining; B. Desmin staining; C.Merge; D. Nuclear staining control; E. Rabbit-derived Desmin antibody was not added as control 图 1 成肌细胞鉴定(100×) Fig. 1 Myoblast identification(100×)
2.2 CCK8检测不同浓度5-Aza-dC对成肌细胞增殖的影响

本试验在成肌细胞中添加不同浓度的甲基化转移酶抑制剂5-Aza-dC后,用CCK8检测0、24、48 h成肌细胞增殖情况。表 3可以看出,与对照组相比,成肌细胞中添加0.05或0.1 μmol·L-1 5-Aza-dC时细胞存活率呈上升趋势,但是随着浓度增加细胞存活率逐渐下降,其中0.3 μmol·L-1 5-Aza-dC浓度下细胞存活率极显著下降(P < 0.01);处理24 h后,0.1 μmol·L-15-Aza-dC处理细胞的存活率呈显著上升(P < 0.05),在0.1 μmol·L-1存活率最高。处理48 h后,细胞存活趋势与24 h相同,但是同24 h相比提高不明显。表明0.1 μmol·L-1 5-Aza-dC处理组适合用于后续试验。

表 3 不同浓度5-Aza-dC对细胞存活率的影响 Table 3 Effects of different concentrations of 5-Aza-dC on cell survival rate
2.3 流式细胞仪检测成肌细胞凋亡

图 2可知,0.1 μmol·L-1 5-Aza-dC处理的细胞凋亡率较空白组极显著增加(P < 0.01)。为进一步检测甲基化转移酶抑制剂5-Aza-dC是否与BclBaxCaspase-9等细胞凋亡因子有关,通过qRT-PCR检测,0.1 μmol·L-1 5-Aza-dC处理的细胞能够极显著的降低细胞抗凋亡因子Bcl的表达(P < 0.01),极显著的提高促凋亡因子Bax表达(P < 0.01),显著提高促凋亡因子Caspase-9的表达(P < 0.05)。

A.空白组细胞凋亡;B.试验组细胞凋亡;C.细胞凋亡率;D.细胞凋亡因子mRNA的表达。Q1.机械性损伤细胞;Q2.晚期凋亡细胞;Q3.正常细胞;Q4.早期凋亡细胞 A. Apoptosis in blank group; B. Apoptosis in experimental group; C. Apoptosis rate; D. Expression of apoptosis factor mRNA. Q1. Cells with mechanical injury; Q2. Non-viable apoptotic cell; Q3. Normal cells; Q4. Viable apoptotic cell 图 2 5-Aza-dC对关岭牛成肌细胞凋亡及相关因子表达的影响 Fig. 2 Effects of 5-Aza-dC on apoptosis and expression of related factors in Guanling cattle myoblasts
2.4 流式细胞仪检测成肌细胞周期

通过流式细胞仪检测0.1 μmol·L-1 5-Aza-dC处理的细胞及空白组细胞,结果显示(图 3),0.1 μmol·L-15-Aza-dC处理组G0/G1期细胞比例较高于空白组,但差异不显著(P>0.05),S期及G2/M期细胞低于空白组,差异也不显著(P>0.05)。为进一步检测甲基化转移酶抑制剂5-Aza-dC对细胞周期相关因子Cyclin B1、Cyclin A2、Cyclin D表达的影响,通过qRT-PCR检测相关细胞因子的表达,发现5-Aza-dC显著的提高了周期因子Cyclin A2的表达(P < 0.05),对细胞因子Cyclin B1、Cyclin D的表达影响不显著(P>0.05)。

A.空白组细胞周期;B.试验组细胞周期;C.细胞周期不同时期比例;D细胞周期因子mRNA的表达 A. Cell cycle in blank group; B. Cell cycle in experimental group; C. The ratio of different stages of the cell cycle; D. Expression of Cyclin mRNA 图 3 5-Aza-dC对关岭牛成肌细胞周期及相关因子表达的影响 Fig. 3 Effects of 5-Aza-dC on myoblasts cell cycle and expression of related factors in Guanling cattle
2.5 qRT-PCR检测MyoD1、DNMT1 mRNA相对表达量

2.5.1 PCR扩增结果   由图 4可知,PCR扩增产物条带清晰,没有拖带和二聚体,片段大小符合试验预期。

A.关岭牛MyoD1基因荧光引物电泳图;B.关岭牛DNMT1基因荧光引物电泳图。M. DNA相对分子质量标准 A. Fluorescence primer electrophoresis image of MyoD1 gene in Guanling cattle; B. Fluorescence primer electrophoresis image of DNMT1 gene in Guanling cattle. M. DL2000 marker 图 4 荧光引物PCR产物扩增结果 Fig. 4 PCR product amplification results of fluorescent primers

2.5.2 Real-time PCR检测结果扩增曲线和峰图   由图 5图 6图 7可知设计的引物通过实时荧光定量PCR检测,扩增曲线拐点清晰,熔解曲线均为单峰,无杂峰,引物特异性好。

A. MyoD1扩增曲线;B. MyoD1熔解峰 A. MyoD1 amplification curve; B. MyoD1 melt peak 图 5 MyoD1基因扩增曲线和峰图 Fig. 5 Amplification curve and melt peak of MyoD1 gene
A. DNMT1扩增曲线;B. DNMT1熔解峰 A. DNMT1 amplification curve; B. DNMT1 melt peak 图 6 DNMT1基因扩增曲线和峰图 Fig. 6 Amplification curve and melt peak of DNMT1 gene
A. GAPDH扩增曲线;B. GAPDH熔解峰 A. GAPDH amplification curve; B.GAPDH melt peak 图 7 GAPDH基因扩增曲线和峰图 Fig. 7 Amplification curve and melt peak of GAPDH gene

2.5.3 MyoD1、DNMT1基因表达分析   由图 8可知,空白组中MyoD1的mMRNA的表达量极显著低于试验组的表达量(P < 0.01),DNMT1的表达量也极显著低于试验组(P < 0.01)。

A. MyoD1相对表达量;B. DNMT1相对表达量 A. Relative expression level of MyoD1;B. Relative expression level of DNMT1 图 8 5-Aza-dC对MyoD1、DNMT1基因表达影响 Fig. 8 Effects of 5-Aza-dC on MyoD1 and DNMT1 gene expression
2.6 MyoD1启动子区甲基化检测

表 4可知,处理48 h后,空白组的MyoD1启动子区CpG甲基化率为28.0%,CHG甲基化率为2.0%,CHH甲基化率为1.9%,试验组的MyoD1启动子区CpG甲基化率为19.8%,CHG甲基化率为1.3%,CHH甲基化率为1.8%。并且0.1 μmol·L-15-Aza-dC处理的试验组极显著的降低了MyoD1启动子区的甲基化率(P < 0.01,图 9)。

表 4 5mC位点覆盖统计 Table 4 Coverage statistics of 5mC sites
图 9 MyoD1基因启动子甲基化率 Fig. 9 MyoD1 promoter methylation rate
3 讨论

DNA甲基化作为一种相对稳定的修饰状态,在DNA甲基转移酶的作用下,可随DNA的复制过程遗传给新生的子代DNA,是一种重要的表观遗传机制[21]。DNA甲基化通过DNA甲基转移酶(DNA methylthansferase,DNMT)来实现,并分为2类,即维持DNA甲基化转移酶(DNMT1)和从头甲基化酶[22-24]。随着对DNA甲基化的继续深入研究,发现一些可以抑制甲基化酶的药物,甲基化酶抑制剂是一种嘧啶核苷类似物,可以致使多种抑癌基因重新激活恢复功能,重新发挥抗肿瘤的作用,以及诱导细胞凋亡,抑制肿瘤细胞的生长,因而有其广阔的临床应用前景,应用较广泛的是甲基化抑制剂5-氮杂-2甲基脱氧胞苷(5-Aza-dC)。利用药物手段调节DNA甲基化水平可以帮助了解基因甲基化对其的影响。

细胞凋亡是细胞的一种基本生物学现象,在多细胞生物去除冗余或异常细胞中起着必要的作用,在生物体进化、内环境稳定以及多个系统发育中起着重要作用[25-27]。Kiianitsa等[28]发现,PARP1在5-Aza-dC处理的人成纤维细胞中形成共价DNA加合物并可以诱导凋亡,阻断DNA复制。Shire等[29]发现,5-Aza-dC恢复了SULF1阴性细胞系中SULF1 mRNA的表达,与此相关的是硫酸酯酶活性的增加和促进肝癌细胞凋亡。李瑾[30]研究发现,5-Aza-dC通过影响DNA甲基化水平,引起细胞周期阻滞,使细胞发生凋亡从而具有抗肿瘤作用。

吴元等[31]用5-Aza-dc处理T24细胞,结果发现,Bcl2蛋白表达降低,Bax、Cleaved Caspase-3表达上升,但不显著。Kiianitsa等[32]发现,在5-Aza-dC处理的人成纤维细胞中PARP1形成共价的DNA加合物并诱导细胞凋亡。而本试验在牛成肌细胞中添加不同浓度的甲基化转移酶抑制剂5-Aza-dC后,用CCK8检测0、24、48 h细胞增殖情况。结果显示,0.1 μmol·L-1 5-Aza-dC处理组是本试验组中最适浓度并用于后续试验,通过流式细胞仪对细胞凋亡和细胞周期进行检测,发现0.1 μmol·L-1 5-Aza-dC处理组能够极显著促进成肌细胞凋亡(P < 0.01);对细胞凋亡因子检测表明,该浓度极显著地降低细胞抗凋亡因子Bcl的表达(P < 0.01),极显著地提高了促凋亡因子Bax表达(P < 0.01),显著提高了促凋亡因子Caspase-9的表达(P < 0.05);这表明5-Aza-dC能够极显著促进关岭牛成肌细胞的凋亡。同时对细胞周期的检测发现,5-Aza-dC对细胞周期的影响差异不显著(P>0.05),为进一步验证5-Aza-dC对细胞周期的影响,通过qRT-PCR发现,周期因子Cyclin A2表达显著提高(P < 0.05),而细胞因子Cyclin B1、Cyclin D表达未发生显著变化(P>0.05)。李凤珍[33]研究发现,低浓度范围的5-Aza-dC(0.005 μmol·L-1)不会显著的影响细胞周期。李瑾[30]研究发现,10 μmol·L-1的5-Aza-dC处理人肺癌细胞A549 48 h后,人肺癌细胞A549细胞周期阻滞于G2/M期,且细胞内的活性氧水平显著升高(P < 0.05)。本试验发现,使用0.1 μmol·L-1 5-Aza-dC处理后,周期因子Cyclin A2的表达显著提高(P < 0.05),推测可能是由于5-Aza-dC浓度偏高。细胞凋亡和细胞周期都是维持体内细胞数量动态平衡的基本措施以及维持胚胎正常发育、机体健康的基本生物学现象[34-37]。本研究发现,5-Aza-dC能够促进细胞凋亡,但是对于细胞周期影响不明显,对于5-Aza-dC的最适浓度和时间还要进一步验证。

MyoD1基因在调节肌肉细胞的生长、代谢及促进肌肉细胞的增殖、分化等功能上均具有的重要作用[38]。Megiorni等[39]发现,5-Aza-dC处理后的DNA去甲基化能够上调miR-378a-3p水平,提高MyHC基因表达量。MyoD1同是重要的肌肉发育基因,通过0.1 μmol·L-1 5-Aza-dC处理细胞后发现,试验组极显著的降低MyoD1启动子区的甲基化率(P < 0.01)。而空白组中MyoD1的mRNA的表达量极显著低于试验组(P < 0.01),同时也发现,DNMT1的表达量极显著低于试验组(P < 0.01)。可能由于5-Aza-dC的作用主要是在DNA复制过程中能够与甲基转移酶(DNMTs) 结合形成共价复合物而达到去甲基化,能够明显抑制该酶的甲基转移活性但是并不影响DNMTs表达。本研究结果提示,5-Aza-dC通过调节Caspase-9、BclBaxCyclin A2 mRNA表达水平影响牛成肌细胞的增殖和凋亡,并且低浓度水平5-Aza-dC(0.1 μmol·L-1)极显著降低了MyoD1启动子区甲基化水平。

4 结论

本研究发现,5-Aza-dC能够通过改变Caspase-9、BclBaxCyclin A2等基因的表达来调节关岭牛成肌细胞的增殖、凋亡。检测发现,0.1 μmol·L-1 5-Aza-dC能极显著降低MyoD1基因启动子区甲基化水平(P < 0.01),极显著提高其mRNA的表达量(P < 0.01)。因此,低浓度水平5-Aza-dC(0.1 μmol·L-1)能有效降低成肌细胞中MyoD1启动子区甲基化水平进而提高MyoD1 mRNA的表达量。所以,低浓度的5-Aza-dC能通过促进细胞凋亡及调控关岭牛MyoD1的甲基化水平来调控MyoD1的表达;同时,MyoD1启动子甲基化可作为遗传标记,为后续关岭牛的品种改良提供理论依据。

参考文献
[1]
ZHAO J G, LU Z G, WANG L, et al. Plant responses to heat stress: physiology, transcription, noncoding RNA, and epigenetics[J]. Int J Mol Sci, 2020, 22(1): 117. DOI:10.3390/ijms22010117
[2]
SENGUPTA A, ROY S S, CHOWDHURY S. Non-duplex G-quadruplex DNA structure: a developing story from predicted sequences to DNA structure-dependent epigenetics and beyond[J]. Acc Chem Res, 2021, 54(1): 46-56. DOI:10.1021/acs.accounts.0c00431
[3]
MADHANI H D. Unbelievable but true: epigenetics and chromatin in fungi[J]. Trends Genet, 2021, 37(1): 12-20. DOI:10.1016/j.tig.2020.09.016
[4]
HARRISON R E S, WENG K G, WANG Y X, et al. Phase separation and histone epigenetics in genome regulation[J]. Curr Opin Solid State Mater Sci, 2021, 25(1): 100892. DOI:10.1016/j.cossms.2020.100892
[5]
ZHANG B, BAN D M, GOU X, et al. Genome-wide DNA methylation profiles in Tibetan and Yorkshire pigs under high-altitude hypoxia[J]. J Anim Sci Biotechnol, 2019, 10(1): 25. DOI:10.1186/s40104-019-0316-y
[6]
鲍正攀. 调控miR-148a基因表达及5-Aza-dc处理对水牛乳腺上皮细胞增殖及泌乳相关基因表达的影响[D]. 南宁: 广西大学, 2018.
BAO Z P. Effect of regulation of miR-148a gene expression and 5-Aza-dc on the proliferation of buffalo mammary epithelial cells and their expression of genes related to laction[D]. Nanning: Guangxi University, 2018. (in Chinese)
[7]
陈观娣. 5-Aza-dC和放疗对宫颈癌细胞系的增殖抑制作用与P16、MGMT基因表达的关系探讨[D]. 广州: 南方医科大学, 2018.
CHEN G D. The inhibition effect of 5-Aza-dC and radiotherapy on cervical cancer cell lines proliferation and the correlation with p16 and MGMT gene expression[D]. Guangzhou: Southern Medical University, 2018. (in Chinese)
[8]
GREVILLE G, LLOP E, HOWARD J, et al. 5-AZA-dC induces epigenetic changes associated with modified glycosylation of secreted glycoproteins and increased EMT and migration in chemo-sensitive cancer cells[J]. Clin Epigenetics, 2021, 13(1): 34. DOI:10.1186/s13148-021-01015-7
[9]
XIE M Y, YANG Y, LIU P, et al. 5-aza-2'-deoxycytidine in the regulation of antioxidant enzymes in retinal endothelial cells and rat diabetic retina[J]. Int J Ophthalmol, 2019, 12(1): 1-7.
[10]
于翔, 王晨柳, 文逸凡, 等. DNA甲基化的表观遗传机制及其在动物遗传育种中应用[J]. 中国牛业科学, 2020, 46(5): 59-62.
YU X, WANG C L, WEN Y F, et al. Epigenetic mechanism of DNA methylation and its application in animal genetics and breeding[J]. China Cattle Science, 2020, 46(5): 59-62. DOI:10.3969/j.issn.1001-9111.2020.05.016 (in Chinese)
[11]
高丽. Myostatin基因编辑牛肌肉卫星细胞成分化过程中DNA甲基化修饰的作用机制研究[D]. 呼和浩特: 内蒙古大学, 2019.
GAO L. The effect of DNA methylations on the muscle satellite cell myogenic differentiation of myostatin gene edited cattle[D]. Hohhot: Inner Mongolia University, 2019. (in Chinese)
[12]
LI M Z, WU H L, LUO Z G, et al. An atlas of DNA methylomes in porcine adipose and muscle tissues[J]. Nat Commun, 2012, 3: 850. DOI:10.1038/ncomms1854
[13]
李秀金, 王金辉, 黄运茂, 等. DNA甲基化影响不同猪种胚胎期肌肉发育差异的研究进展[J]. 中国畜牧杂志, 2020, 56(6): 24-29.
LI X J, WANG J H, HUANG Y M, et al. Research progresses on the role of DNA methylation in regulating distinct embryonic muscle development among pig breeds[J]. Chinese Journal of Animal Science, 2020, 56(6): 24-29. (in Chinese)
[14]
万火福, 韦凤英, 廖玉英, 等. 鸡肌肉组织基因组DNA甲基化与肉质性状相关性分析[J]. 中国家禽, 2018, 40(12): 5-9.
WAN H F, WEI F Y, LIAO Y Y, et al. Study on the associations of genomic DNA methylation in muscle tissue with meat quality of chicken[J]. China Poultry, 2018, 40(12): 5-9. (in Chinese)
[15]
WEBER M, HELLMANN I, STADLER M B, et al. Distribution, silencing potential and evolutionary impact of promoter DNA methylation in the human genome[J]. Nat Genet, 2007, 39(4): 457-466. DOI:10.1038/ng1990
[16]
WARDLE F C. Master control: transcriptional regulation of mammalian Myod[J]. J Muscle Res Cell Motil, 2019, 40(2): 211-226. DOI:10.1007/s10974-019-09538-6
[17]
刘会敬. 中外猪种胚胎不同发育时期骨骼肌转录组和全基因组甲基化分析[D]. 武汉: 华中农业大学, 2018.
LIU H J. Integration of DNA methylome and transcriptome analysis of prenatal skeletal muscles between Tongcheng and Yorkshire pigs[D]. Wuhan: Huazhong Agricultural University, 2018. (in Chinese)
[18]
曹阳. 杜寒杂交羊与小尾寒羊背最长肌全基因组甲基化与转录组比较分析[D]. 长春: 吉林大学, 2017.
CAO Y. Comparative analysis on genome-wide methylation and RNA-seq in longissimus dorsi muscle between dorper×small tailed Han crossbred and small tailed Han sheep[D]. Changchun: Jilin University, 2017. (in Chinese)
[19]
薛倩. 京海黄鸡腿肌全基因组甲基化和转录组分析研究[D]. 扬州: 扬州大学, 2017.
XUE Q. Supported by China agriculture research system (CARS-42) and priority academic program development of Jiangsu higher education institutions[D]. Yangzhou: Yangzhou University, 2017. (in Chinese)
[20]
沈兰, 陈春华, 叶冬云, 等. ESR1、MYOD1和hTERT基因启动子DNA甲基化与宫颈癌放化疗敏感性的关系[J]. 河北医药, 2020, 42(1): 33-37.
SHEN L, CHEN C H, YE D Y, et al. Relationship between DNA methylation of ESR1, MYOD1, hTERT gene promoter and sensitivity to radiotherapy and chemotherapy in patients with cervical cancer[J]. Hebei Medical Journal, 2020, 42(1): 33-37. (in Chinese)
[21]
BORCHIELLINI M, UMMARINO S, DI RUSCIO A. The bright and dark side of DNA methylation: a matter of balance[J]. Cells, 2019, 8(10): 1243. DOI:10.3390/cells8101243
[22]
刘胜辉. 黑色素瘤相关抗原MAGE-A家族在喉鳞状细胞癌中的表达及临床意义[D]. 石家庄: 河北医科大学, 2016.
LIU S H. Expressions and significance of melanoma antigen-A family in laryngeal squamous cell carcinoma[D]. Shijiazhuang: Hebei Medical University, 2017. (in Chinese)
[23]
VENTÓS-ALFONSO A, YLLA G, MONTAÑES J C, et al. DNMT1 promotes genome methylation and early embryo development in cockroaches[J]. iScience, 2020, 23(12): 101778. DOI:10.1016/j.isci.2020.101778
[24]
IGUCHI E, TAKAI A, TAKEDA H, et al. DNA methyltransferase 3B plays a protective role against hepatocarcinogenesis caused by chronic inflammation via maintaining mitochondrial homeostasis[J]. Sci Rep, 2020, 10(1): 21268. DOI:10.1038/s41598-020-78151-2
[25]
LONG Y, UBYCH K, JAGU E, et al. FRET-based method for direct, real-time measurement of DNA methyltransferase activity[J]. Bioconjug Chem, 2020, 32(1): 192-198.
[26]
D'ARCY M S. Cell death: a review of the major forms of apoptosis, necrosis and autophagy[J]. Cell Biol Int, 2019, 43(6): 582-592. DOI:10.1002/cbin.11137
[27]
MAJTNEROVÁ P, ROUŠAR T. An overview of apoptosis assays detecting DNA fragmentation[J]. Mol Biol Rep, 2018, 45(5): 1469-1478. DOI:10.1007/s11033-018-4258-9
[28]
KⅡANITSA K, ZHANG Y B, MAIZELS N. Treatment of human cells with 5-aza-dC induces formation of PARP1-DNA covalent adducts at genomic regions targeted by DNMT1[J]. DNA Repair (Amst), 2020, 96: 102977. DOI:10.1016/j.dnarep.2020.102977
[29]
SHIRE A, LOMBERK G, LAI J P, et al. Restoration of epigenetically silenced sulfatase 1 expression by 5-Aza-2'-deoxycytidine sensitizes hepatocellular carcinoma cells to chemotherapy-induced apoptosis[J]. Med Epigenet, 2015, 3(1): 1-18. DOI:10.1159/000375461
[30]
李瑾. 5-氮杂-2-脱氧胞苷对牛胎儿成纤维细胞和人肺癌细胞A549的影响及作用机理的研究[D]. 呼和浩特: 内蒙古大学, 2015.
LI J. The effect of 5-Aza-dC-2'-deoxycytidine on bovine fetal fibroblasts and human lung cancer A549 associated mechanism[D]. Hohhot: Inner Mongolia University, 2015. (in Chinese)
[31]
吴元, 张文涛, 马文超, 等. SOX7启动子甲基化对膀胱癌增殖及凋亡的影响[J]. 医学研究生学报, 2019, 32(12): 1285-1290.
WU Y, ZHANG W T, MA W C, et al. Effect of SOX7 promoter methylation on proliferation and apoptosis of bladder cancer[J]. Journal of Medical Postgraduates, 2019, 32(12): 1285-1290. (in Chinese)
[32]
KⅡANITSA K, ZHANG Y B, MAIZELS N. Treatment of human cells with 5-aza-dC induces formation of PARP1-DNA covalent adducts at genomic regions targeted by DNMT1[J]. DNA Repair (Amst), 2020, 96: 102977. DOI:10.1016/j.dnarep.2020.102977
[33]
李凤珍. 5-氮-2-脱氧胞苷对牛成纤维细胞周期、染色体和凋亡的影响[D]. 武汉: 华中农业大学, 2009.
LI F Z. 5-Aza-2-deoxyctidine on cell cycle, chromosome and apoptosis of fibroblast[D]. Wuhan: Huazhong Agricultural University, 2009. (in Chinese)
[34]
COLELLA F, SCILLITANI G, PIERRI C L. Sweet as honey, bitter as bile: mitochondriotoxic peptides and other therapeutic proteins isolated from animal tissues, for dealing with mitochondrial apoptosis[J]. Toxicology, 2021, 447: 152612. DOI:10.1016/j.tox.2020.152612
[35]
WHITE T L, DESHPANDE N, KUMAR V, et al. Cell cycle re-entry and arrest in G2/M phase induces senescence and fibrosis in Fuchs Endothelial Corneal Dystrophy[J]. Free Radic Biol Med, 2021, 164: 34-43. DOI:10.1016/j.freeradbiomed.2020.12.445
[36]
刘玉思. 神经管畸形动物模型中凋亡调控因子(miR-322、PCSK9)的表达规律及作用机制的研究[D]. 沈阳: 中国医科大学, 2020.
LIU Y S. The study on the expression and mechanism of apoptosis regulatory factors (miR-322, PCSK9) in animal model of neural tube defect[D]. Shenyang: China Medical University, 2020. (in Chinese)
[37]
刘丹. Rno_circ_0005139/miR-324-3p/Wnt5a调控大鼠肛门直肠畸形胚胎发育过程中细胞增殖、凋亡及作用机制的研究[D]. 沈阳: 中国医科大学, 2020.
LIU D. Rno_circ_0005139/miR-324-3p/Wnt5a regulates cell proliferation, apoptosis and the mechanism in the development of rat anorectal malformations[D]. Shenyang: China Medical University, 2020. (in Chinese)
[38]
DUMONT N A, RUDNICKI M A. Characterizing satellite cells and myogenic progenitors during skeletal muscle regeneration[J]. Methods Mol Biol, 2017, 1560: 179-188.
[39]
MEGIORNI F, CIALFI S, MCDOWELL H P, et al. Deep Sequencing the microRNA profile in rhabdomyosarcoma reveals down-regulation of miR-378 family members[J]. BMC Cancer, 2014, 14: 880. DOI:10.1186/1471-2407-14-880

(编辑   郭云雁)