李 晶,张梦瑶,刘开东,柳 楠,贺建宁
(1.青岛农业大学 动物科技学院,山东 青岛 266109;2.曲阜市畜牧兽医技术服务中心,山东 曲阜 273100;3.青岛畜牧兽医研究所,山东 青岛 266121)
敖汉细毛羊是我国很早培育出的细毛羊品种,它主要分布在东北部地区,是我国最具有代表性的细毛羊之一[1],其羊毛是纺织工业的重要原料。敖汉细毛羊毛囊的结构和特性直接决定毛发品质和商业价值[2],毛囊拥有着相对复杂的结构,它的形成是由多种细胞来协同完成的。毛囊通过一系列的分化生长最后成为羊毛,其生长发育也有些复杂,它是由多种因子来相互调节,因此,通过研究相关因子来探讨毛囊的生长发育显得尤为重要。
Hox家族在生物体内是一类与发育调控相关十分重要的基因。Hox家族这类转录因子都是高度保守的[3],对毛囊细胞增殖分化也起着十分重要的作用[4-5]。研究发现,在皮肤和毛囊生长发育过程中Hox家族基因起到了主导作用[6],Hoxa5主导了皮肤细胞增殖、分化与凋亡[5],对毛囊退行期进行了调控,并诱导了细胞分化[7]。动物形态的进化和多样性与Hox基因数目、序列及其表达方式密切相关,已有研究结果表明,当Hox家族基因发生变化时,动物在被毛生长时会表现出缺陷[8]。BMP6基因是TGF-β超家族中的一员[9-10],TGF通路中的许多转录因子都受BMP6调控,均表现出上调和下调趋势。杨燕燕等[11]在生殖发育上对BMP6进行了研究,它控制着羊的多胎,尤其是在发情期基因表达相当明显。也有研究表明,BMPs通路中的转录组因子对毛囊发育主要表现出抑制作用,而它的下游靶标BMPRs主要表现出促进的作用[12-13]。程旭等[14]通过对铁调素在BMP6-HJV-SMAD信号通路的作用,验证了BMP6的重要性。由此可见,Hoxa5和BMP6基因都能抑制毛囊生长,但二者之间的关系鲜有报道。
对于家畜皮肤毛囊的一些研究大多集中在单个基因的功能验证上,还有一些是与骨的调节、繁殖等与皮肤毛囊无关的研究,本试验主要的侧重点是利用同源重组技术构建表达载体通过与单、共转染比较基因表达量的变化来研究Hoxa5和BMP6之间的作用关系,旨在为更深层次上研究毛囊功能基因奠定基础。
健康的30日龄胎羊由青岛畜牧所奥特种羊场提供;TRIzol、KpnⅠ、倒置显微镜、荧光定量PCR等均购自Bio-Rad公司[9]。
1.2.1 引物设计 选取绵羊Hoxa5基因(GenBank登录号为NM-001009431.1)CDS区 和BMP6基因(GenBank登录号为NM-001110277.1)CDS区,软件导入pcDNA3.1表达载体的序列,通过Snap Gene设计Hoxa5基因和BMP6基因上下游引物,序列分别如下:
F(Hoxa5):5′-ATGCGTTAACTCAGGAATCGGATGCATCGGTCAAGCTGGCAT-3′;R(Hoxa5):5′-TGGCCTTAAAATCTCCGAGGAAATACTCGCTCACGCGGTAG-3′。F(BMP6):5′-TTACGGGCATCGGTAGGGAAATTCGTATTGCATGCAAAGCT-3′;R(BMP6):5′-TGTCCGAGACCTGCACATCCACGGCGGCCCTTCAGCATACCCCTA-3′。扩增体系:20 μL,cDNA 1 μL,上下游引物各1 μL,Mix 17 μL。
1.2.2 pcDNA3.1与Hoxa5、BMP6基因同源重组 对Hoxa5、BMP6进行切胶,用KpnⅠ对pcDNA3.1进行单酶切,体系:KpnⅠ 1 μL,pcDNA3.1 5 μL,Buffer 5 μL,ddH2O 39 μL。通过电泳对BMP6、Hoxa5基因片段及pcDNA3.1进行胶回收,SoSo连接体系:BMP6、Hoxa5基因上、下游引物各1 μL,pcDNA3.1 1 μL,SoSo 5 μL,ddH2O 3 μL,转化DH5α,过夜对菌液进行振荡培养,将菌液生长良好的送往公司进行测序[12-13]。
1.2.3 胎儿成纤维细胞培养 在细胞室将30日龄胎羊从母羊子宫内取出,先用75%的酒精对胎儿进行消毒,再用PBS将胎儿表面的酒精冲洗掉。在培养皿中用手术刀将胎儿的头部和四肢、内脏等全部剔除,只用胎羊的躯干,将其剪碎,分布在培养皿中,在每一个肌肉块上面加2滴已经预热的胎牛血清,放进CO2培养箱(37 ℃,5.0% CO2),24 h对其进行观察并且换液[15-16]。
1.2.4 细胞转染 细胞进行传代以后,等密度再次达到95%左右时,借助转染试剂对重组的载体进行瞬时转染。取脂质体20 μL,PBS 480 μL,加入到无菌的离心管中并混匀,静置5 min。进一步进行转染,转染时分为2组,分别是共转染组和单转染组,共转染组中加入Hoxa5、BMP6质粒各20 μL,单转染组加入Hoxa5质粒40 μL或者BMP6质粒40 μL,2组均加入PBS 460 μL,具体详细操作步骤参照Lipofectamine2000说明书[17]。
1.2.5 细胞转染后的qPCR 重组载体转染成纤维细胞后需要进一步的增殖培养,当密度达到95%左右时,对成纤维细胞的RNA进行提取,并将其反转录成cDNA[18]。使用Snap Gene软件对Hoxa5、BMP6和GAPDH基因进行qPCR引物设计,序列见表1。
表1 引物信息Tab.1 Primers information
qPCR反应条件:95 ℃反应7 min;95 ℃ 40 s,60 ℃退火40 s,35个循环;72 ℃ 30 s,4 ℃保存。每组3个生物学重复,对数据利用Ct(2-ΔΔCt)公式进行计算,得出Hoxa5、BMP6的相对表达量。借助SPSS 20.0软件对数据进行差异显著性的分析,以P<0.01作为差异显著性判断标准[17]。
1.2.6 Hoxa5、BMP6蛋白表达的检测 提取成纤维细胞中的组蛋白,进行Western Blot,首先进行转膜,后进行封闭,时长2 h;进行一抗、二抗的孵育,均对抗体进行2 000倍的稀释。孵育完成后利用1∶1的曝光液进行曝光拍照[17]。
以cDNA为模板,添加引物进行PCR 扩增获取Hoxa5、BMP6基因片段。1%琼脂糖凝胶电泳结果表明(图1),Hoxa5基因条带在图中804 bp处,BMP6基因条带在图中537 bp处,都与目的基因CDS区大小吻合,可用于后期的试验。
M.2000 Marker;1.Hoxa5基因扩增片段(804 bp);2.BMP6基因扩增片段(537 bp)。M.2000 Marker;1.Hoxa5 product(804 bp);2.BMP6 product(537 bp).
使用KpnⅠ将pcDNA3.1表达载体从环状酶切成链状,用1%琼脂糖凝胶电泳检测是否酶切成功(图2)。1,2是对pcDNA3.1的KpnⅠ单酶切,pcDNA3.1长度为5 428 bp,图中条带所在位置约是5 428 bp,结果显示酶切完成,回收条带进行重组。
图2 pcDNA3.1载体单酶切鉴定Fig.2 Single enzyme digestion identification of pcDNA3.1 vector
将pcDNA3.1酶切成链状后,用SoSo试剂盒将BMP6、Hoxa5分别与pcDNA3.1相连接,转化,振荡培养,进行菌液PCR验证如图3所示,在约804 bp处有明亮电泳条带,是Hoxa5基因;537 bp处的是BMP6基因。菌液测序结果用Blast进行对比,结果如图4,说明BMP6、Hoxa5没有碱基的突变,已成功的连接到pcDNA3.1载体上。将连接好的重组载体从菌液中提取出来,如图5所示,1-2为pcDNA3.1-BMP6重组质粒,大小为5 965 bp,3-4为pcDNA3.1-Hoxa5重组质粒,大小为6 232 bp,与图中条带相符。
M.2000 Marker;1,4.BMP6菌液PCR;2,3.Hoxa5菌液PCR。M.2000 Marker;1,4.BMP6 bacterial PCR;2,3. Hoxa5 bacterial PCR.
对胎儿成纤维细胞的生长情况进行观察记录(图6),24 h能够观察到细胞已经开始贴到培养皿底部,周围也有一些漂浮的组织。培养至96 h后,密度再次达到95%左右,此时即可对原代细胞进行传代[18]。
对引物进行验证,1%琼脂糖凝胶电泳结果(图7)表明,引物扩增效果良好,没有出现二聚体等杂物,在128 bp 处的是内参基因GAPDH,在186 bp处的是Hoxa5基因,在116 bp处的是BMP6基因,综上可得引物均可用于实时荧光定量PCR试验。
图4 BMP6(上)和Hoxa5(下)表达载体菌液的测序比对Fig.4 Sequencing comparison of BMP6 (up) and Hoxa5 (down) expression vector bacilli
M.10000 Marker;1-2.pcDNA3.1-BMPR6;3-4.pcDNA3.1-Hoxa5。M. 10000 Marker;1-2. pcDNA3.1-BMPR6 recombinant plasmids;3-4. pcDNA3.1-Hoxa5 recombinant plasmids
图 8,9分别为Hoxa5、BMP6和GAPDH的qPCR过程中的扩增曲线、熔解曲线,两线效果显示均良好,在qPCR过程中没有出现非特异性产物。对2组的mRNA表达进行差异显著性分析(图10),分析显示Hoxa5与BMP6共转染后的Hoxa5基因相对表达量极显著低于单转染(P<0.01),而共转染后的BMP6基因相对表达量在共转染却极显著高于单转染(P<0.01)。
图6 成纤维细胞原代(左)、传代(右)培养(×100)Fig.6 Primary (left) and passage (right) culture of fibroblasts (×100)
A:M.500 Marker;1-2.内参基因GAPDH;B:M.500 Marker;1-2.目的基因BMP6;C:M.500 Marker;1-2.目的基因Hoxa5。A:M.500 Marker;1-2 .Reference gene GAPDH;B:M. 500 Marker;1-2.Target gene BMP6;C:M.500 Marker;1-2. Target gene Hoxa5.
图8 Hoxa5、BMP6和GAPDH 熔解曲线Fig.8 Hoxa5,BMP6 and GAPDH melting curves
图9 Hoxa5、BMP6和GAPDH扩增曲线 Fig.9 Hoxa5, BMP6 and GAPDH amplification curve
对Hoxa5、BMP6基因的共转染、单转染蛋白的表达(图11)进行比较分析,从灰暗度以及条带的大小,使用ImageJ和SPSS软件进行分析,Hoxa5与BMP6共转染后的Hoxa5蛋白表达量极显著低于单转染过程(P<0.01),而共转染的BMP6蛋白表达量却极显著高于单转染过程(P<0.01)(图12),这与mRNA相对表达量的结果一致。
不同大写字母表示差异极显著(P<0.01)。图12同。Different capital letters indicate extremely significant difference(P<0.01).The same as Fig.12.
图11 Hoxa5、BMP6蛋白表达 Fig.11 Hoxa5 and BMP6 expression protein
图12 Hoxa5、BMP6蛋白相对表达量Fig.12 Relative expression of Hoxa5 and BMP6 protein
许多研究表明,毛囊的发育依靠一系列的信号分子,这些信号分子在真皮细胞和上皮细胞间穿梭着,介导了这2个细胞群体间的相互作用,毛囊能够稳定的进行发育全靠这两者的介导。毛囊的形成过程中离不开所谓的原始信号,此信号来源于真皮下[18],并且在真皮下去激活Wnt信号[19-20]。目前,许多信号通路都参与毛囊的发育,最常见的有Wnt通路、FGF家族、BMP家族、TGF家族[21]等。Botchkarev等[22]利用小鼠来研究BMPs在皮肤中的表达位置以及表达时期,表明BMPs家族基因对毛囊的发育存在着一定的抑制作用。BMP6基因和BMP4一样同属于BMP家族,是一种抑制毛囊发育的信号分子,Hoxa5基因是TNF家族中尤为重要的基因,同样也是一种抑制毛囊发育的信号分子,两基因对毛囊的发育均有抑制作用[23]。在信号通路中,BMP家族中的信号分子对TNF通路中信号分子有干涉作用,有些是辅助基因的表达,有些是抑制基因的表达。现在对单个基因功能的验证普遍存在,多个基因相互作用还是尤为少见,本研究通过利用同源重组手段来验证多基因相互之间的作用,帮助筛选相关基因,有助于后期育种和羊毛生产。
为探讨Hoxa5、BMP6基因在细胞水平上是否存在互作效应,本研究利用同源重组技术,快速构建了pc-DNA3.1-Hoxa5、pc-DNA3.1-BMP6表达载体,并通过共转染、单转染后基因mRNA、蛋白表达量的比较分析,发现Hoxa5基因有促进BMP6基因表达,BMP6基因有抑制Hoxa5基因表达的现象。本研究在构建表达载体过程中所用同源重组技术较传统的T克隆技术节约了时间,减少了碱基突变;同时由分子水平到细胞水平得到了Hoxa5和BMP6基因的互作关系,一方面验证了与羊毛生长有关的基因,为信号通路的解析添加了论证,另一方面为进一步在个体水平上验证其功能奠定基础。