胚胎发育是哺乳动物种群繁衍的关键环节。在植入后胚胎发育的原肠运动期间,上胚层细胞发生迁移和重排,形成具有不同发育潜能的外、中、内三胚层,为胚胎轴建立和器官发育奠定细胞基础
[1]。临床研究表明,孕期流产高发于胚胎发育第5~12周,而5周前由于胚胎较小缺乏临床数据,但预估流产率可达50%
[2]。胚胎发育异常是导致孕期流产的重要因素之一。小鼠三胚层分化缺陷引起胚胎发育阻滞,最终导致发育致死
[3]。探究三胚层分化的调控机制对阐明孕期流产风险具有重要意义。
转录因子在胚胎发育和胚层分化中扮演着关键角色。其通过识别并结合特定的脱氧核糖核酸(deoxyribonucleic acid,DNA)序列,调节下游靶基因的转录活性,进而指导细胞分化和组织形成调控胚胎发育进程。例如,转录因子配对盒基因6(paired box 6,PAX6)是神经系统发育的关键调控因子,其缺失导致神经系统发育异常,影响神经前体细胞的生成和分化
[4]。Brachyury(T)是中胚层形成必不可少的转录因子,它的表达下降导致中胚层结构缺失,如脊索、尿囊和原条缺失
[5]。SRY相关HMG盒基因17(SRY-related HMG-box gene 17,SOX17)敲除导致内胚层发育严重异常,无法形成前肠和中肠结构,胚胎在受精后第8.5天至第9.5天发育致死
[6]。由此可见,转录因子通过参与胚层分化,对维持胚胎发育至关重要。
转录因子SP/KLF家族在果蝇、大鼠、家兔、人和斑马鱼等物种中广泛表达,参与细胞增殖、凋亡、分化等生物学过程
[7-8]。SP/KLF家族的特异性蛋白1(specificity protein 1,SP1)是第1个被鉴定和纯化的转录因子
[9],与复发性流产、胎儿生长受限等疾病密切相关
[10-11]。果蝇中SP1/EGR样锌指蛋白HUCKEBEIN缺失,导致内胚层原基发育受损,外胚层和中胚层原基异常扩增
[12]。小鼠胚胎中
Sp1缺失导致胚胎发育严重迟缓,引起体轴不对称、翻转不完全以及卵黄囊异常生长等多种发育异常,最终致死于胚胎发育第11天
[13]。尽管SP1高表达于第8.5天胚胎的滋养巨细胞、神经与心脏组织
[14]。但其是否在三胚层分化时期表达,以及如何调控三胚层分化尚不明确。
由于细胞谱系复杂、细胞数量稀少,以及缺乏相应的小量化染色质免疫共沉淀测序(chromatin immunoprecipitation-sequencing,ChIP-seq)技术,小鼠早期胚胎发育研究一直是领域内的重点与难点。为此,本研究优化了ChIP-seq技术,并利用其检测转录因子SP1在三胚层中的基因组分布状态,结合Sp1敲除的小鼠胚胎模型,探究了SP1在三胚层分化过程中所发挥的作用。
1 材料与方法
1.1 实验动物
C57BL/6N雌鼠(4~6周龄)购自北京维通利华实验动物技术有限公司。假孕雌鼠与结扎雄鼠均为昆明鼠(6~10周龄),购自重庆医科大学实验动物中心。小鼠均饲养于本校实验动物中心独立通风笼盒(individual ventilated cage,IVC)环境,自由饮食取水。动物饲养和实验程序严格遵守重庆医科大学实验动物福利伦理相关规定。
1.2 主要试剂及仪器
本研究涉及的试剂与仪器信息如下:小鼠胚胎M2操作液(M7167,Sigma-Aldrich)、KSOM培养液(MR-101-D,Millipore)、SP1抗体(21962-1-AP,Proteintech)、ChIP-seq文库制备试剂盒(E7645S,NEB)、Smart-seq2文库制备试剂盒(TD502,Vazyme)、荧光定量仪Qubit4(Invitrogen)、电穿孔仪(Bio-Rad)、37 ℃加热台(Jeio Tech)。
1.3 Sp1 单向导RNA(single guide RNA,sgRNA)及验证引物
根据小鼠Sp1基因序列,利用IDT网站的CRISPR靶点在线设计工具设计4条sgRNA:5'-CCTAGAACTGGTGAACCACC-3'、5'-AGCTTAATGGTGCCCAAGAC-3'、5'-GTTTGTGGTACCAAGACACT-3'、5'-CTGACCCCTTGGTTTGCAAT-3'。sgRNA交由金斯瑞生物科技有限公司合成。为验证Sp1基因是否成功敲除,通过NCBI Primer designing tool网站设计两对引物,用于聚合酶链式反应(polymerase chain reaction,PCR)扩增。外侧引物(上游:5'-CCAAACCTGGGTTCTAATCCCT-3'、下游:5'-TTTCTCATATGGTGCTGAACATGG-3'),内侧引物(上游:5'-CCAAACCTGGGTTCTAATCCCT-3'、下游:5'-CGCTACCCCCATTATTGCCA-3')。引物交由北京擎科生物科技股份有限公司合成。
1.4 Sp1敲除胚胎构建及基因型鉴定
对C57BL/6N雌鼠进行超促排卵。腹腔注射10 IU孕马血清促性腺激素(pregnant mare serum gonadotropin,PMSG)以促进卵泡发育,46~48 h后注射10 IU人绒毛膜促性腺激素(human chorionic gonadotropin,HCG)诱发排卵。与C57BL/6N雄鼠合笼,次日检查阴道栓确认交配成功。对孕鼠进行颈椎脱臼处死,并取出受精卵。通过CRISPR-Cas9结合受精卵电转核酸酶(zygote electroporation of nucleases,ZEN)
[15]技术导入sgRNA,敲除
Sp1的第2和第3号外显子。将受精卵体外培养至2-细胞期后,移植至第0.5天假孕雌鼠输卵管中。在胚胎发育至第7.5天处死代孕雌鼠,分离胚外脏内胚层,加入10 μL DNA模板制备溶液(H
2O 39.47 μL、5×Lysis buffer 10 μL、蛋白酶K 0.53 μL)提取DNA。在55 ℃孵育60 min后,95 ℃孵育15 min灭活蛋白酶K。
接着进行目的片段扩增,反应体系为:2×Phanta Flash Master Mix 5 μL、H2O 3 μL、DNA模板 1 μL、上游引物0.5 μL、下游引物0.5 μL。反应条件:98 ℃ 30 s;(98 ℃ 15 s,62 ℃ 15 s,72 ℃ 45 s)×32个循环;72 ℃ 5 min。取10 μL PCR扩增产物加入琼脂糖凝胶(1.5%)上样孔中,经150 V电泳30 min后,使用凝胶成像系统成像。
1.5 第7.5天胚胎胚层分离
取孕7.5 d雌鼠,颈椎脱臼处死,从子宫剥离胚胎。按照《小鼠胚胎操作实验手册》第3版所述的胚层分离方法,分离第7.5天胚胎
[16]。首先去除胚胎周围的Reichert's膜,接着切除外胎盘锥,用酶消化胚区和胚外区,使用口吸管分离胚外外胚层(extraembryonic ectoderm,ExE)、胚外脏内胚层(extraembryonic visceral endoderm,ExVE)、内胚层(endoderm,End),再持玻璃针剥离中胚层(mesoderm,Mes),切下原条(primitive streak,PS),胚胎剩余部分为外胚层(ectoderm,Ect)。
1.6 ChIP-seq文库构建和测序
实验方法参照2016年STAR ChIP-seq并做修改
[17]。首先,收集胚胎样本,加入细胞裂解液于冰上孵育30 min以破碎细胞释放染色质,经微球菌核酸酶处理5 min。终止酶解后,用SP1抗体进行孵育,免疫沉淀DNA。之后按照NEB试剂盒E7645S说明书进行文库构建,使用Qubit 4荧光定量仪测定DNA浓度。将构建好的DNA文库送至北京诺禾致源公司测序。
1.7 RNA-seq文库构建和测序
实验方法参考Smart-seq2
[18]。首先,收集样本并加入0.2%Triton X-100提取反转录信使RNA(messenger RNA,mRNA),通过oligo(dT)反转录mRNA成补脱氧核糖核酸(complementary deoxyribonucleic acid,cDNA)。之后按照Vazyme试剂盒TD502说明书进行文库构建;使用Qubit 4荧光定量仪测定DNA浓度。将构建好的DNA文库送至北京诺禾致源公司进行测序。
1.8 ChIP-seq数据处理
首先将测序所得读段(Reads)通过Bowtie2 v2.2.2映射到小鼠mm9基因组上,参数设置为-t-q-N1-L25。去掉不能映射至基因组的读段、映射到多个位置的读段以及PCR导致的重复读段。在验证生物学重复之间的可重复性后,将不同重复的数据整合分析。之后将全基因组分为100 bp长度序列框,统计每100 bp区域内的信号富集,并使用每千碱基每百万测序读数的读段数(reads per kilobase per million of sequenced reads,RPKM)均一化信号。数据可视化采用加利福尼亚大学圣克鲁斯分校(University of California,Santa Cruz,UCSC)提供的在线软件(
http://genome.ucsc.edu/)。
1.9 RNA-seq数据处理
首先将测序所得读段通过TopHat v2.0.11映射到小鼠mm9基因组上。去掉不能映射到基因组的读段。之后基于refFlat提供的基因组注释信息,使用cufflinks v2.0.2计算每千碱基每百万测序片段的片段数(fragments per kilobase per million of sequenced reads,FPKM)。数据的可视化采用UCSC大学提供的在线软件(
http://genome.ucsc.edu/)。
2 结 果
2.1 SP1在小鼠第7.5天胚胎三胚层中的分布特征
首先,建立小鼠胚胎胚层分离技术并验证分离纯度。采用酶消化联合机械分离的方式获取小鼠胚胎第7.5天的胚层组织(
图1A),各设置两个生物学重复。利用RNA-seq检测各胚层基因表达水平。结果显示各胚层均特异性表达该组织的标记基因
[16,19],表明胚层分离获取的纯度较高(
图1B)。
之后,利用STAR ChIP-seq检测SP1在各胚层中的分布情况。SP1的信号富集模式在2个生物学样品之间高度一致,且信噪比良好(
图2A)。两生物学重复的相关系数约为0.7~0.8(
图2B),证明样本之间重复性较好。在后续分析中,本研究将不同生物学重复得到的测序数据进行合并,以提高SP1在基因组上的检测效率。
结果显示,SP1大量富集在基因启动子区域,在基因主体与基因间区分布较少。在Ect、Mes、End、PS与ExE中启动子上的分布比例分别为57.29%(4 596/8 022)、47.48%(3 574/7 527)、44.21%(1 699/3 843)、66.24%(2 943/4 443)与42.49%(7 677/18 067)。尽管ExE中SP1在启动子上的分布比例较低,但其信号峰数量相对其他胚层更多(
图3)。
2.2 Sp1敲除导致发育基因表达下降
为探究SP1在胚层发育中的功能,利用CRISPR-Cas9结合ZEN技术在受精卵时期敲除
Sp1(
图4),转移至代孕小鼠输卵管中,并于第7.5天收获对照组与基因组胚胎。
将胚胎分别进行基因分型后,首先在Ect中验证了
Sp1敲除后其信号富集变化。结果显示SP1信号峰的富集强度大量下降(
图5)。之后,采用RNA-seq检测
Sp1缺失后的基因表达变化。结果表明,相较于其他胚层,在Mes与ExE中
Sp1缺失对基因表达的影响最为明显,差异表达基因数量更多(
图6)。
由于SP1是转录激活因子,本研究着重关注RNA表达水平受其影响而下降的基因。基因本体(gene ontology,GO)分析结果显示,在不同胚层组织中,
Sp1敲除引起表达下调的基因通常参与细胞分化、器官发育或组织形态发生等生物学进程。例如,Ect中表达下调的基因与神经分化、细胞发育相关;Mes中下调基因参与心脏发育、细胞命运决定;End中下调基因调节多细胞器官发育、肺形态发育;PS中下调基因参与与前后轴形态决定、体节发生;ExE中下调基因与细胞迁移、血管生成等生物学进程相关(
图7)。
2.3 SP1调控胚胎发育的潜在靶基因
结合ChIP-seq与RNA-seq数据,本研究进一步探究了SP1调控胚层分化的下游靶基因。将符合以下条件的基因鉴定为靶基因:①
Sp1敲除后基因表达下降;②基因启动子区域具有SP1结合位点。结果显示,SP1在ExE中的靶基因数目较多(
图8A)。
Gbx2可能是Ect中受SP1直接调控的下游基因,而
Vim、
Enkur、
Tbx6、
Id3等基因则分别是Mes、End、PS与ExE的潜在靶基因(
表1)。UCSC截图的结果验证了这些基因启动子区域存在SP1结合信号,同时基因表达在敲除后下降(
图8B)。
3 讨 论
转录因子对维持哺乳动物胚胎发育与指导细胞命运决定至关重要。SP1作为广泛表达的转录因子,在细胞周期、DNA修复、组织发育及多种肿瘤进程中均发挥重要功能。然而,现有研究大多集中于成体组织或特定细胞系,SP1在早期胚胎及胚层分化中的作用尚未明确。本研究基于小量化ChIP-seq技术与胚层分离策略,率先解析了SP1在小鼠第7.5天胚胎基因组中的结合模式。通过对基因组信号分布特征的分析,发现SP1在不同胚层中均富集于基因启动子区域。敲除Sp1基因后,参与器官发育、形态发生等生物学过程的发育基因,其RNA表达水平明显下降,提示SP1作为启动子特异性转录激活因子,可能在胚层分化过程中发挥关键作用。进一步整合ChIP-seq和RNA-seq数据,鉴定了多个潜在受SP1直接调控的靶基因,如Gbx2、Tbx6和Vim等,为解析SP1在早期胚胎发育中的转录调控功能提供了新的证据。
在Ect中SP1结合并调控
Gbx2的表达。
Gbx2作为神经外胚层标志物,在中枢神经系统发育中扮演关键角色,其缺失可导致神经嵴细胞的迁移障碍及分化异常。斑马鱼中Sp5l是类似SP1的转录因子,
Sp5l敲低同样导致GBX2的表达活性降低
[20]。在PS中,SP1靶向调控
Tbx6基因,该基因在左右轴建立及体节发育过程中发挥核心作用,其突变小鼠表现为体节消失及异位神经管形成。
SP1在ExE中结合峰的数量最多,且
Sp1缺失后该组织中RNA表达下调的基因数目高于其他胚层。ExE是母体与胚胎相互作用的关键组织,其功能异常可能导致胎盘发育受损及妊娠失败。先前研究报道胎盘中SP1表达下调与人类复发性流产相关
[21],本研究进一步鉴定了SP1在ExE中的潜在靶基因,如
Id3和
Ddx3x等。在胎盘发育过程中,ID3介导表皮生长因子受体(epidermal growth factor receptor,EGFR)信号途径调控细胞增殖、血管生成和细胞因子产生
[22],而
Ddx3x缺陷胚胎呈现胎盘发育异常,导致胚胎死亡
[23],提示SP1可能通过调控血管生成、细胞迁移相关靶基因,影响胚胎着床及胎盘发育。综上所述,这些结果为明确SP1在三胚层分化及胚胎发育中的调控作用提供了新的支撑。
然而,SP1的转录调控网络可能涉及多个协同调控因子,其相互作用关系尚不清楚。已有研究表明,SP1可通过与Wnt/β-catenin信号通路协同调控发育相关基因,同时其活性可能受到蛋白激酶C及组蛋白修饰(如H3K27me3)的影响
[24-26]。因此,需进一步解析SP1如何与其他转录因子或表观遗传因子协同作用,从而调控胚层特异性基因的表达。此外,虽然本研究鉴定出多个潜在受SP1调控的发育基因,但其表达水平是否直接受SP1调控,以及其在胚胎发育中的生物学功能仍需进一步验证。
重庆市教育委员会重大资助项目(KJZD-M202300402)