查荒源 于光晴 陈小云 肖武汉* 刘 兴*

(1. 大连海洋大学水产与生命学院, 大连 116023; 2. 中国科学院水生生物研究所, 武汉 430072)

低氧诱导因子(Hypoxia inducible factor, HIF)介导的转录激活调控是细胞感受氧气变化和应答低氧胁迫的最为关键的信号传导途径[1,2]。低氧信号传导途径非常保守, 从线虫到人类, 低氧信号都启动相同或相似的同源基因的转录和调控, 从而调控类似的生理和生化反应[3]。HIF由HIF-α(包括HIF-1α和HIF-2α)和HIF-1β两个亚单位组成, 每个亚单位都包含基本的bHLH-PAS (helix-loop-helix-PAS)结构域, 该结构域介导异源二聚体的形成及与DNA的结合。HIF-1β的蛋白水平稳定, 因而HIF的活性主要由HIF-α来决定[4,5]。在常氧条件下, HIF-α能够被一类依赖于氧气的脯氨酸羟化酶(Prolyl hydroxylase domain enzymes, PHDs, 包括PHD1, PHD2 和PHD3)所识别, 进而对HIF-α(包括HIF-1α和HIF-2α)上特定的脯氨酸残基进行羟基化修饰, 其中,PHD2是起主要作用的脯氨酸羟化酶[6]。被羟基化修饰的HIF-α能被pVHL(Von Hippel-Lindau)识别,pVHL能够与ElonginB、ElonginC结合形成E3泛素连接酶复合体, 介导HIF-α通过蛋白酶体降解[7]。在低氧条件下, 由于氧气的缺乏, 使得脯氨酸羟化酶的活性受到抑制, HIF-α的羟基化修饰受到抑制, 进而导致HIF-α的蛋白酶体降解途径受到抑制, 从而导致HIF-α的蛋白水平得以积累, 积累的HIF-α进入细胞核, 与HIF-1β形成异源二聚体, 该二聚体与下游基因启动子上的低氧应答元件(Hypoxia response element, HRE)结合, 从而激活下游基因的表达, 发挥相应的生物学功能[1,4]。

蛋白质赖氨酸残基上的甲基化修饰主要是由一类含有SET 保守结构域的蛋白质完成的, 目前在哺乳动物中已鉴定到50多种蛋白质可以归入这一家族, 组蛋白甲基转移酶SET9就是该家族的成员[8,9]。SET9可以甲基化修饰组蛋白H3上第4位的赖氨酸残基, 也可以对多种非组蛋白进行甲基化修饰[10—12]。SET9可以甲基化修饰p53上第372位的赖氨酸残基,从而激活p53的转录活性[13]。SET9可以甲基化RelA第37位的赖氨酸残基, 在肿瘤坏死因子TNFα的刺激下增强RelA与下游基因启动子的结合[14]; SET9也可以甲基化RelA第314和315位的赖氨酸残基, 从而促进与下游基因启动子结合的RelA的降解, 抑制NF-κB信号通路的激活[15]。

前期的研究工作发现SET9可以甲基化HIF-1α上第32位的赖氨酸残基, 抑制HIF-1α与下游基因启动子的结合, 从而抑制低氧信号通路下游基因的表达[16]。但是, SET9在体内调控低氧信号通路和在低氧胁迫中的生物学功能及其机制还并不清楚。

研究基因在体生物学功能最好的手段是建立该基因敲除的动物模型。斑马鱼(Danio rerio)是重要的模式生物, 建立目的基因敲除的斑马鱼品系有助于研究该基因的在体生物学功能。前期的研究工作以斑马鱼为模式生物, 建立了foxo3b、tet1及fih敲除的斑马鱼品系, 用以研究这些基因在低氧耐受中的生物学功能。foxo3b可以诱导vhl的表达, 从而促进hif-α的降解, foxo3b的敲除会导致斑马鱼低氧耐受能力的下降[17]; tet1可以稳定hif-α的蛋白水平,tet1的敲除会降低斑马鱼的低氧耐受能力[18]; fih可以抑制hif-α的转录激活活性, fih的敲除可以增强斑马鱼的低氧耐受能力[19]。

本研究利用CRISPR/Cas9基因编辑技术, 获得了甲基转移酶set9敲除的斑马鱼品系, 初步分析了set9基因的敲除对斑马鱼低氧耐受能力的影响, 为研究鱼类甲基转移酶set9在低氧耐受中的生物学功能及耐低氧鱼类新品种的培育提供了理论依据。

1 材料与方法

1.1 主要试剂及仪器

主要试剂 DNA胶回收试剂盒购自Bioflux公司, 质粒小提试剂盒购自天根公司, 限制性内切酶、RNAiso plus购自TaKaRa, All-in-One cDNA Synthesis SuperMix kit试剂盒、SYBR Green Fast quantitative PCR master mix试剂盒购自Biotool公司, T4 DNA连接酶、DNase Ⅰ购自Fermentas公司,T7体外转录试剂盒购自Ambion公司, Agarose、2×Sosoo Mix重组酶购自Tsingke公司, DNA Marker、Trans10感受态细胞购自TransGen公司, 多聚甲醛、DEPC购自Amresco公司。

主要仪器 显微注射仪PLI-100(Harvard Aparatus), 斑马鱼生化培养箱SHP-150(上海精宏),低氧工作站Invivo2(Ruskinn), 凝胶成像系统(上海培清科技有限公司), 离心机5702、5430、5425D、5415R(Eppendorf), 荧光定量PCR仪Step One Real-Time PCR System(ABI), 恒温水浴锅(上海精宏),Nanodrop 2000浓度定量仪(Thermo), -80℃超低温冰箱(Thermo)。

1.2 斑马鱼品系

AB品系的斑马鱼购自国家斑马鱼资源中心。

1.3 序列分析

使用ClustalX2软件对人SET9、小鼠Set9和斑马鱼set9的氨基酸序列进行比对, 并根据已知的人SET9的结构域标注3个物种的SET结构域及其甲基转移酶的酶活性位点。

1.4 斑马鱼set9基因的敲除与鉴定

利用在线设计工具 (http://crispr.mit.edu) 设计Cas9作用的靶位点; 将体外合成的Cas9 mRNA和gRNA(guide RNA)按照400 ng/μL﹕60 ng/μL混合后,注射到一细胞期胚胎动物极; 24h后, 利用异双链迁移率测定法(Heteroduplex mobility assay, HMA)筛选有敲除效果的F0代嵌合体; F0代性成熟后与AB品系斑马鱼交配, 筛选到能稳定遗传的F1代杂合子;F1代性成熟后自交可获得对应的野生型和纯合子突变体, 利用PCR扩增基因组中包含Cas9作用靶位点的DNA片段, 经测序确定突变的碱基信息。

1.5 实时荧光定量PCR(qRT-PCR)

取受精后3d的野生型及set9基因敲除的斑马鱼幼鱼, 使用TaKaRa公司的RNAiso plus进行RNA的提取, 使用Biotool公司的All-in-One cDNA Synthesis SuperMix kit试剂盒进行cDNA的合成, 使用Biotool公司的SYBR Green Fast quantitative PCR master mix试剂盒进行qRT-PCR检测, 使用5′-GCATGT GCTGGGTCTACTATC-3′和5′-GGCCATCAGGG TACACATAAG-3′扩增set9的mRNA,β-actin作为内参, 其扩展引物为5′-TACAATGAGCTCCGTGTT GC-3′和5′-ACATACAATGGCAGGGGTGTT-3′。

1.6 斑马鱼低氧处理

使用Ruskinn Invivo2 I-400低氧工作站进行斑马鱼的低氧胁迫处理。对于受精后3d的斑马鱼幼鱼, 低氧胁迫时氧气浓度调整到2%, 将出膜且发育程度一致的野生型及基因敲除的斑马鱼置于6 cm培养皿中, 每皿50条, 加入5 mL的胚胎培养液（Egg water）, 放置于低氧工作站中, 每12h统计1次鱼的存活情况, 利用GraphPad 7.0绘制存活曲线; 对于3月龄的斑马鱼成鱼, 低氧胁迫时氧气浓度调整到5%, 将个体大小相近(个体重量±0.02 g)的野生型及基因敲除的斑马鱼放置于250 mL锥形瓶中, 加入200 mL的曝气水, 放置于低氧工作站中, 并进行跟踪拍照, 记录鱼的运动及存活情况。

1.7 TUNEL染色

将个体大小相近(个体重量±0.02 g)的3月龄野生型及基因敲除的斑马鱼放置于低氧工作站, 并在低氧胁迫6h后将对应的野生型及敲除品系的斑马鱼进行脑部的取材和固定, 利用TUNEL染色试剂盒进行染色并进行拍照, 比较斑马鱼脑部细胞凋亡的情况。

1.8 统计学处理

所有相关实验至少重复3次。使用GraphPad 7.0软件展示统计的数据, 结果表示为平均值±SEM。对于所有比对,P＜0.05被视为具有统计学意义。

2 结果

2.1 斑马鱼甲基转移酶set9的序列特征

斑马鱼set9基因的全长CDS长度为1977 bp, 编码373个氨基酸组成的蛋白质。利用ClustalX2软件对人SET9、小鼠Set9和斑马鱼set9的氨基酸序列进行比对, 发现序列的保守性很高。作为一个含SET结构域的赖氨酸甲基转移酶, 斑马鱼set9也含有一个保守的SET结构域(图 1中虚线框标注)并且其发挥甲基转移酶活性的酶活性位点在三个物种间都是保守的(图 1中实线框标注)。

图1 人SET9、小鼠Set9、斑马鱼set9的氨基酸序列对比Fig. 1 Alignment of amino acid sequences of human SET9, mouse Set9 and zebrafish set9

2.2 建立set9基因敲除的斑马鱼品系

根据CRISPR/Cas9基因编辑技术的原理, 在斑马鱼set9基因的第一个外显子上, 设计了基因编辑的靶位点, 其序列为5′-GGACAGCGATGATGAC AACA-3′(图 2A)。提取基因编辑品系鱼尾基因组,使用引物5′-AGTGCAGGGCTGATGATG-3′和引物5′-CTAATCCTCGACGGAATGAGAC-3′扩增包含靶位点的DNA片段, 经异双链迁移率测定法筛选获得有效果的F0代(图 2B)。

F0发育至性成熟后与AB品系交配, 筛选到能稳定遗传的F1代, 获得一个缺失8个碱基的品系。由于在翻译起始密码子之后, 缺失8个碱基, 造成移码突变, 终止密码子提前出现, 使得突变的品系只表达不包含SET结构域的14个氨基酸的多肽(图 2A)。由于翻译的提前终止, 会导致mRNA的不稳定性,进而导致mRNA的快速降解, 因此, 进一步利用qRT-PCR技术分析野生型和突变品系中set9的mRNA水平, 结果显示突变品系中set9的mRNA水平显著下降(图 2C)。

图2 利用CRISPR/Cas9技术敲除斑马鱼set9基因Fig. 2 Targeting strategy for generating zebrafish set9 mutant by CRISPR/Cas9

以上结果表明, 利用CRISPR/Cas9基因编辑技术, 获得了set9基因敲除的斑马鱼品系。在实验室正常饲养条件下,set9基因敲除的斑马鱼与其同胞野生型相比, 发育、生长和繁殖等并无肉眼可见的差异。

2.3 set9的敲除能增强斑马鱼低氧耐受能力

利用获得的set9基因敲除的斑马鱼品系, 分析set9基因在斑马鱼低氧耐受中的功能。首先, 选取出生后3d的斑马鱼幼鱼进行低氧处理。该时期的斑马鱼幼鱼刚孵化出膜, 可以方便地获取足够数量的样本, 且个体差异较小, 因而可以很好地评估斑马鱼低氧耐受能力的差异。将野生型及set9基因敲除的斑马鱼幼鱼放置于氧气含量为2%的低氧工作站中, 随着时间的推移, 记录斑马鱼幼鱼的存活时间, 制作存活曲线。结果显示, 在低氧胁迫处理24h后, 野生型的斑马鱼已经开始死亡, 而set9基因敲除的斑马鱼在处理36h后才出现明显的死亡; 在低氧胁迫处理60h后, 野生型的斑马鱼均已死亡, 而set9敲除的斑马鱼还有部分存活(图 3)。由此可见,set9基因的敲除可以增强斑马鱼幼鱼的低氧耐受能力。

图3 set9基因的敲除能增强斑马鱼幼鱼的低氧耐受能力Fig. 3 Deletion of set9 enhances hypoxia tolerance in zebrafish larvae

进一步, 选取3月龄的斑马鱼成鱼进行低氧耐受能力的比较。将野生型及set9基因敲除的斑马鱼成鱼放置于氧气含量为5%的低氧工作站中, 录制视频记录斑马鱼的存活及运动的情况。结果显示,在低氧胁迫处理4h后, 野生型的斑马鱼即开始出现死亡; 在低氧胁迫处理8h后, 野生型的斑马鱼均已死亡, 而set9敲除的斑马鱼则刚开始有鱼出现濒死的迹象(图 4)。以上结果显示,set9基因的敲除可以增强斑马鱼成鱼的低氧耐受能力。

图4 set9基因的敲除能增强斑马鱼成鱼的低氧耐受能力Fig. 4 Deletion of set9 enhances hypoxia tolerance in zebrafish larvae

2.4 set9的敲除减轻低氧导致的脑细胞凋亡

为阐明set9敲除的斑马鱼高低氧耐受力的机制, 利用TUNEL染色检测低氧胁迫处理后野生型及set9敲除的斑马鱼脑中细胞凋亡的情况。5%O2处理6h后, 选择存活的野生型及set9敲除的斑马鱼, 对其脑组织取材、固定、切片并进行TUNEL染色。结果表明, 野生型斑马鱼脑中凋亡细胞数量明显高于突变体(图 5)。

图5 set9基因的敲除能减少低氧胁迫下的斑马鱼脑部的细胞凋亡Fig. 5 Deletion of set9 diminishes hypoxia-induced apoptosis in adult zebrafish brains

3 讨论

SET9是一种重要的赖氨酸甲基转移酶, SET9最开始报道为通过甲基化组蛋白H3的第4位赖氨酸残基来实现对基因表达的调控[11]。此外, 多种非组蛋白也被发现是SET9的底物, 包括p53、TAF10、E2F1、AR、RelA和FOXO3a等[12]。SET9通过调控这些非组蛋白的甲基化修饰参与了细胞对多种环境胁迫的响应: SET9对FOXO3a的甲基化修饰可以调控氧化损伤诱导的神经细胞死亡[20,21]; SET9通过甲基化修饰p53和E2F1等转录因子的活性来响应DNA损伤等[13,22,23]。最近, 首次在细胞水平报道了SET9通过直接甲基化低氧信号通路中的关键转录因子HIF-1α第32位的赖氨酸残基, 进而调控HIF-1α的转录活性, 从而在低氧信号转导及细胞低氧胁迫应答中发挥重要的生物学功能[16]。然而, 对于SET9在体内的生物学功能还并不清楚。

水体溶氧量随季节、地理位置、水流状态和水体环境等变化而发生很大的变化, 鱼类在长期演化历程中发展出适应不同溶氧水环境的策略[3]。而另一方面, 随着“工厂化”养殖等集约化养殖模式的发展, 单位水体养殖密度和产量的提高, 也迫切地要求养殖品种的耐低氧能力有相应的提高。对于鱼类低氧适应机制的研究, 不仅对于认识鱼类物种的分化及低氧信号传导途径的系统演化等方面具有十分重要的理论意义, 而且对于培育耐低氧鱼类新品种具有十分重要的实践意义。低氧信号传导途径非常保守, 近年来的研究表明其在鱼类低氧适应及低氧耐受中具有重要功能[24—26]。鉴于低氧诱导因子及其调控的分子网络在低氧信号传导中的关键作用, 挖掘低氧诱导因子的调控因子并解析其调控模式的分子机制是阐明鱼类低氧耐受分子机制的主要手段, 可以为培育耐低氧鱼类新品种提供候选靶标。

在哺乳动物中, SET9通过甲基化HIF-1α第32位的赖氨酸残基, 从而实现对低氧信号通路的调控。SET9在不同物种间具有极高的保守性, 其发挥甲基转移酶活性的SET结构域和酶活性位点, 在不同物种间都是非常保守的。然而, 斑马鱼中HIF-1α第32位的赖氨酸残基并不保守, 并不具备SET9作用的经典底物所具备的“R/K S K”3个氨基酸组成的基序[13,16]。因此, 斑马鱼set9调控低氧信号传导途径及影响鱼类低氧耐受的分子机制还有待进一步的探讨, 很可能与哺乳类SET9调控低氧信号传导的机制有所不同。

本研究利用CRISPR/Cas9基因编辑技术, 获得了甲基转移酶set9敲除的斑马鱼品系。斑马鱼set9基因的敲除可以增强斑马鱼的低氧耐受能力,减少低氧胁迫导致脑组织中的细胞凋亡水平, 且在正常饲养条件下, 其发育、生长和繁殖与野生型斑马鱼相比, 并无差异。本研究为深入了解甲基转移酶set9的在体生物学功能及其机制提供了线索, 也为鱼类耐低氧新品种的培育提供了候选的靶标。