低氧指令因子1(hypoxia inducible factor， HIF-1)算作参与缺氧反应的中枢因子凡俗存在于动物体内， 它不错调控多种靶基因的抒发。据报谈， 高档动物HIF能调控100多个基因， 这些基因功能各类， 如调控血管生成、机体氧均衡、能量代谢、细胞凋一火和B淋巴细胞发育等[1-3]。HIF-1是由HIF-1α和HIF-1β两个亚基组成的异源二聚体， 其中HIF-1α亚基对氧气明锐， 在常氧细胞中， 可通过脯氨酸羟化酶和泛素卵白酶体将其降解， 而在低氧环境下黄胖系列， 脯氨酸羟化酶活性被扼制， HIF-1α正经抒发， 进而插足细胞核内与HIF-1β勾搭成异源二聚体。二聚体复合物与低氧反应原件(hypoxia responsive elements， HREs)勾搭， 再与一些转录因子共同作用调控下贱靶基因， 引起机体的生理生化反应[4]。在东谈主(Homo sapiens)、牦牛(Bos grunniens)、高原鼢鼠(Eospalax fontanierii)和大鼠(Rattus norvegicus)等[5-8]高档动物中， 已有好多HIF的报谈。水产动物的计议研究主要聚拢在鱼类。如武昌鱼(Megalobrama amblycephala)[9]和胭脂鱼(Myxocyprinus asiaticus)[10]等处于缺氧的环境条目下时， 一些组织中HIF-1α基因的转录水平会权贵上调。河鲈(Perca fluviatilis)、黄颡鱼(Pelteobagrus fulvidraco)、半滑舌鳎(Cynoglossus semilaevis)也有计议报谈[11-13]。无脊椎动物HIF-1的计议研究较少， 海洋软体动物中惟一太平洋牡蛎(Crassostrea gigas)[14]和杂色鲍(Haliotis diversicolor)[15]等少数几个物种HIF-1基因的计议报谈。Kawabe等[14]通过Northern blot和Western blot分析讲解， 缺氧条目下太平洋牡蛎HIF-1α mRNA和卵白齐大意被周期性指令抒发。

低氧在海水繁衍环境中频频出现， 常给繁衍生物带来不利影响。诚然， 贝类稳妥缺氧环境的智商相对较强， 但低氧情况下相似会出现免疫下落、滋长受阻和DNA毁伤等知足， 从而导致滋长速率放慢， 生殖力下落， 以致物化[16]。研究显现， 贝类不错通过增强氧的运送智商、降呆板量糟塌、栽植能量供应和免疫力等多种生理反应来搪塞低氧环境[17]。魁蚶(Scapharca broughtonii)是中国首要埋栖型经济贝类， 主要散布在中国、朝鲜半岛及菲律宾沿海等， 生涯在3~50 m水深的软泥或泥沙质海底[18]。赵庆[19]研究标明， 魁蚶的低氧耐受智商较强， 用DO值为0.5 mg/L的条目威胁10 d后， 存活率仍高达15%， 而四角蛤蜊(Mactra quadrangularis)5 d后物化率即达到90%;魁蚶的鳃、外衣膜和血液中SOD酶活性和T-AOC在低氧威胁后也发生了较大变化， 但其基因层面的反应握法和调控机制有待真切研究。

为探究魁蚶HIF-1α的基因结构特征和对低氧的反应， 本研究克隆获取了魁蚶HIF-1α基因的cDNA全长， 并对其组织散布、低氧威胁下的抒发进行了研究， 将为揭示魁蚶搪塞低氧的稳妥机制提供研究尊府， 也为进一步丰富贝类HIF的计议研究提供参考。

1 材料与方法 1.1 试验材料与起原

试验用魁蚶来自于长岛海区， 登科壳长为30 mm傍边的健康个体， 低温运至试验室， 于18℃充气海水中暂养14 d， 每天换水2次， 每次换水量50%， 每4 h投喂单胞藻一次。

1.2 试验经管及试验组织获取

低氧威胁经管参照赵庆试验方法加以矫正[19]:行使蜕变氮气和空气的充气速率抵制DO， 实时监测并限定蜕变以保证氧浓度妥贴试验要求。试验共设4个DO浓度梯度， 分袂为0.5 mg/L、2.5 mg/L、4.5 mg/L和7.5 mg/L(对照组)， 经管0 h、4 h、8 h、16 h、24 h、36 h、48 h和64 h后， 在每组中立时取3个魁蚶的血淋巴、鳃、外衣膜、斧足、闭壳肌和肝胰腺6个组织置于液氮中， 用于RNA索要。基因克隆和组织散布研究材料来自于对照组。

1.3 基因克隆与序列分析 1.3.1 总RNA索要与cDNA合成

RNA索要取舍TRIzol法[20]， 稍加转变。主要圭臬为:将组织样品在液氮中研磨， TRIzol经管10 min， 加入氯仿抽提卵白质， 然后再行使异丙醇千里淀RNA， 用75%酒精连气儿洗净异丙醇后室温干燥， DEPC水熔解。临了测定RNA浓度和行使琼脂糖凝胶电泳检测RNA的质料、纯度和完好意思性。行使PrimeScriptTM Ⅱ 1st Strand cDNA Synthesis Kit试剂盒(TaKaRa)按照阐述书操作圭臬合成cDNA第一链， PrimeScriptTM RT reagent Kit with gDNA Eraser (Perfect Real Time)获取用于qRT-PCR试验的cDNA第一链， 置于‒20℃中备用。

1.3.2 SbHIF-1α基因克隆

凭据试验室已得转录组数据库中的序列， 行使Primer 5分袂遐想基因的两个特异性引物(HIF-1α-F、HIF-1α-R)扩增SbHIF-1α cDNA的中间片断。凭据扩增出的中间片断再分袂遐想出两对5'和3'RACE引物(RHIF- 1α-R1、RHIF-1α-R2)和(RHIF-1α-F1、RHIF-1α-F2)使用RACE试剂盒(Clontech)分袂进行巢式扩增获取基因cDNA全长。第一轮使用试剂盒中的UPM分袂和(RHIF-1α-R1、RHIF-1α-F1)进行四个体系扩增， 再以第一轮扩增家具稀释50倍为模板， 以NUP和(RHIF-1α-R2、RHIF-1α-F2)进行第二轮扩增。第一轮PCR反应圭臬为94℃ 30 s， 72℃ 1 min， 5个轮回; 94℃ 30 s， 70℃ 30 s， 72℃ 1 min， 5个轮回; 94℃ 30 s， 68℃ 30 s， 72℃ 1 min， 20个轮回; 4℃保存。第二轮PCR反应圭臬为94℃ 30 s， 68℃ 30 s， 72℃ 2 min， 20个轮回; 72℃ 2 min; 4℃保存。所用引物见表 1。

表 1 基因克隆和荧光定量PCR所用引物序列 Tab.1 The primers used in gene clone and qRT-PCR experiments

用1.2%的琼脂糖胶对PCR家具进行电泳检测， 对贪图条带进行切胶回收， 贯穿T载体并升沉至克隆菌， 再挑取阳性单克隆菌株， 测序考证。

1.3.3 序列分析

使用SeqMan对测序效果进行拼接， 获取完好意思cDNA序列。通过EditSeq软件翻译序列的ORF， 展望卵白分子量大小及等电点; 行使Interproscan和Smart在线软件(， -heidelberg.de/)展望编码卵白的功能域; 通过同源建模处事器(SWISS-MODEL)对卵白三级结构进行展望; 经Blast()查找同源序列， 行使ClustalX和DnaMan软件进行多序列同源比对， 行使Mega 6.0构建系统进化树。

1.4 荧光定量PCR(qRT-PCR)分析

用上述试验所得cDNA稀释15倍后算作qRT-PCR反应模板， 遐想特异性引物(QHIF- 1α-F、QHIF-1α-R)， 采选β-actin算作内参基因(表 1)。反应体系为TB Green Premix Ex Taq Ⅱ 10 μL正反向引物各0.8 μL， ROX ReferenceDye Ⅱ(50×) 0.4 μL， 模板2 μL， ddH2O 6 μL。反应圭臬为: 95℃ 30 s; 95℃ 5 s， 60℃ 34 s， 40个轮回; 95℃ 15 s， 60℃ 1 min， 95℃ 15 s。应用Applied Biosystems 7500 Fast Real-Time PCR System进行qRT-PCR试验， 每个样品作念3个重迭， 数据效果取舍2‒ΔΔCt法揣度基因的相对抒发量。应用SPSS 17.0软件进行单成分方差分析， P < 0.05时为权贵性相反， P < 0.01时为极权贵相反。取舍GraphPad prism 6.0软件进行作图分析。

2 效果与分析 2.1 基因cDNA序列分析

将测序效果进行拼接比对， 最终得到HIF-1α基因的cDNA全长序列(定名为SbHIF-1α， GenBank登录号为MH936548)(图 1)。

图 1 SbHIF-1α核苷酸序列及编码氨基酸序列 方框内分袂代表肇端密码子(ATG)和停止密码子(TGA)， 单下划线的为加尾信号(aataaa)， 双下划线的氨基酸为糖基化位点， 加粗下划线的氨基酸为脯氨酸羟化酶勾搭位点， poly(A)用海浪线默示. Fig.1 Nucleotide and deduced amino acid sequence of SbHIF-1α The start codon (ATG) and stop codon (TAG) are boxed; polyadenylation signals (aataaa) are underlined with underscore; N-glycosylation sites are double underlined; the proline hydroxylase binding site is overstriking underlined; poly (A) tail are marked by wavy lines.

SbHIF-1α基因cDNA全长为2741 bp， 包括194 bp的5'-UTR， 411 bp的3'-UTR和2136 bp的ORF， ORF编码711个氨基酸残基， 展望卵白分子量为80.8 kDa， 表面等电点为5.57。3'-UTR近poly A隔邻具有加尾信号序列(AATAAA); 软件展望SbHIF-1α含有两个糖基化位点， 一个脯氨酸羟化酶勾搭位点LxxLAP。

2.2 SbHIF-1α的同源性分析

在线比对效果显现， SbHIF-1α编码的氨基酸序列与其他物种HIF-1α编码的氨基酸序列具有较高相似度， 与杂色鲍、好意思洲牡蛎(C. virginica) HIF-1α的相似度分袂为93%和92%， 与东谈主HIF-1α的相似度也高达77%。序列保守性相对较强， 氨基酸对比信息见图 2。

图 2 不同物种HIF-1α氨基酸序列比对 深色袒护默示保守的氨基酸残基， 淡色袒护默示相似的氨基酸残基. 分析所用的氨基酸序列GenBank登录号如下:西伯利亚鲟， KY174955.1;好意思洲牡蛎， HM441076.1;草鱼， AY450269.2;半滑舌鳎， NM_001294187.1;皱纹盘鲍， MH135278.1;东谈主， NM_181054.2;鲢， HM146310.1;团头鲂， GU363498.1;虾夷扇贝， XM_021497943.1;小家鼠， BC026139.1;麦穗鱼， FJ794604.1;非洲爪蟾， NM_001086980.1. Fig.2 Multiple sequence alignment of SbHIF-1α with other hemoglobin amino acid sequences Fully conserved residues are highlighted in dark. Strongly conserved residues among sequences are shaded in light colour. The GenBank accession numbers of protein sequences used for analysis are as follows: Acipenser baerii， KY174955.1; Crassostrea virginica， HM441076.1; Ctenopharyngodon Idella， AY450269.2; Cynoglossus semilaevis， NM_001294187.1; Haliotis discus hannai， MH135278.1; Homo sapiens， NM_181054.2; Hypophthalmichthys molitrix， HM146310.1; Megalobrama amblycephala， GU363498.1; Mizuhopecten yessoensis， XM_021497943.1; Mus musculus， BC026139.1; Pseudorasbora parva， FJ794604.1; Xenopus laevis， NM_001086980.1.

对魁蚶HIF-1α卵白功保守结构域进行分析， 效果见图 3， 魁蚶HIF-1α卵白的N端量对保守， 共有4个保守结构域: 1个HLH结构域， 由第16位氨基酸到第69位氨基酸之间的54个氨基酸残基组成; 2个PAS结构域， 其中PAS-A由第79位氨基酸到第145位氨基酸之间的67个氨基酸残基组成; PAS-B由第214位氨基酸到第280位氨基酸之间的67个氨基酸残基组成; 1个PAC结构域， 由第286位氨基酸到第329位氨基酸之间的44个氨基酸残基组成， 这些保守结构域与图 2深色保守区域相对应。

图 3 魁蚶HIF-1α保守结构域 Fig.3 Domains of hypoxia inducible factor-1α of Scapharca broughtonii

对魁蚶bHIF-1α卵白三级结构进行展望， 如图 4所示。其中一致性最高的区域为22~326位氨基酸残基， 正好是两个图 3展望4个保守结构域场合区域。

图 4 魁蚶HIF-1α三级结构 Fig.4 The tertiary structure of hypoxia inducible factor-1α of Scapharca broughtonii

行使Mega 6.0软件的Maximum likelihood模子构建的进化树如图 5。本研究中的魁蚶先与虾夷扇贝(Mizuhopecten yessoensis)聚为一支， 再与杂色鲍和好意思洲牡蛎汇聚， 它们的亲缘关系最近， 形成一个相对寥寂的大分支， 接着与脊索动物聚为一支， 而与属节肢动物门的棉铃、拟穴青蟹和中华绒螯蟹的亲缘关系则相对较远。这显现魁蚶HIF-1α基因的分子进化地位与其生物学分类地位基本一致。

图 5 Maximum likelihood模子构建的SbHIF-1α系统进化树 进化树中所用卵白的氨基酸序列的GenBank登录号:拟穴青蟹， KU644140.1;中华绒螯蟹， KF825558.1;棉铃虫， KU552106.1;好意思洲牡蛎， HM441076.1;日本鹌鹑， XM_015865743.1;皱纹盘鲍， MH135278.1;东谈主， NM_181054.2;鲢， HM146310.1;团头鲂， GU363498.1;虾夷扇贝， XM_021497943.1;小家鼠， BC026139.1;麦穗鱼， FJ794604.1;非洲爪蟾， NM_001086980.1. Fig.5 The phylogenetic tree based on the sequences of different SbHIF-1α family members The GenBank accession numbers of protein sequences used for phylogenetic analysis are as follows: Scylla paramamosain， KU644140.1; Eriocheir sinensis， KF825558.1; Helicoverpa armigera， KU552106.1; Crassostrea virginica， HM441076.1; Coturnix japonica， XM_015865743.1; Haliotis discus hannai ， MH135278.1; Homo sapiens， NM_181054.2; Hypophthalmichthys molitrix， HM146310.1; Megalobrama amblycephala， GU363498.1; Mizuhopecten yessoensis， XM_021497943.1; Mus musculus， BC026139.1; Pseudorasbora parva， FJ794604.1; Xenopus laevis， NM_001086980.1. 2.3 SbHIF-1α组织抒发分析

qRT-PCR检测 SbHIF-1α基因在血淋巴、鳃、外衣膜、斧足、闭壳肌和肝胰腺6个组织中的抒发效果如图 6所示。从图不错看出， 在所检测的6个组织中均能检测到SbHIF-1α的转录本， 血淋巴中抒发量最高， 鳃组织次之， 闭壳肌中最低。权贵性分析效果显现， SbHIF-1α在外衣膜、闭壳肌、肝胰腺、斧足四个组织中的抒发量莫得权贵相反(P > 0.05)， 相对于闭壳肌， 血淋巴中抒发量是闭壳肌的11.45倍， 达到权贵相反(P < 0.05)， 鳃中的抒发量也达到相抗拒平(P < 0.05);同期， SbHIF-1α在血淋巴和鳃中的抒发量也有权贵相反。

图 6 SbHIF-1α在魁蚶不同组织中的抒发散布 a:血淋巴， b:鳃， c:外衣膜， d:斧足， e:闭壳肌， f:肝胰腺. 不同字母间默示相反权贵(P < 0.05). Fig.6 Distribution of SbHIF-1α gene in different tissues of Scapharca broughtonii a: Haemocyte， b: Gill， c: Mantle， d: Foot， e: Adductor muscle， f: Hepatopancreas. Different letters indicate significant difference (P < 0.05). 2.4 SbHIF-1α在低氧威胁后各组织的抒发分析

SbHIF-1α在低氧威胁后各组织中的抒发变化如图 7所示。算作对照组， DO值为7.5 mg/L的对照组中SbHIF-1α基因在各组织的抒发量低且抒发正经。全体来看， 血淋巴(图 7A)和鳃(图 7B)中SbHIF-1α比较其他组织反应较为积极， 抒发量变化幅度也相对高于其他组织， 最高值分袂为对照组的519.43倍和124倍， 且抒发量随威胁时候加多而总体升高。血淋巴中， 0.5 mg/L、2.5 mg/L的低氧威胁下， SbHIF-1α抒发随时候延伸而迟缓加多， 4.5 mg/L组在36 h急剧高潮达到最高值后又快速回调。鳃中SbHIF-1α的抒发量趋势与血淋巴相似， 但威胁组在16 h时有个小的回调， 36 h运转权贵上调， 64 h达到最高值。而SbHIF-1α这种升高后缩小再升高的变化趋势也体现外衣膜(图 7C)、斧足(图 7D)、闭壳肌(图 7E)和肝胰腺(图 7F)中， 不外这4个组织中SbHIF-1α的变化幅度相对于血淋巴和鳃小好多， 如0.5 mg/L试验组外衣膜、闭壳肌和肝胰腺中SbHIF-1α的最大值分袂是对照组的95.98倍、42.35倍和52.31倍。另外， 6个组织中SbHIF-1α对不同DO经管反应进度具有一定的相似性， 随DO值缩小即威胁进度的增强， SbHIF-1α反应进度愈加浓烈。如血淋巴中， 0.5 mg/L的低氧经管组相较2.5 mg/L和4.5 mg/L这两个低氧组反应更积极， 相对抒发变化更大， 从威胁运转便积极反应， 高潮趋势较2.5 mg/L和4.5 mg/L这两个低氧组更为赫然， 权贵水平在4 h时就达到极权贵水平(P < 0.01)。

图 7 各组织SbHIF-1α在低氧威胁激后mRNA抒发水平变化 合并试验组不同字母默示相反性极权贵(P < 0.01). Fig.7 SbHIF-1α mRNA expression level after hypoxia stress in various organizations Different letters in the same experimental group indicate highly significant difference (P < 0.01). 3 商榷

高档动物耐低氧智商差， 任何万古候的缺氧齐会导致机体缺少以致物化[21]。但是， 某些低等动物如海洋无脊椎动物依然进化到在一段时候内行使极少的氧气也不错存活[22]。为了搪塞缺氧取环境， 生物机体形成了一系列的蜕变机制， 其中HIF-1等于一种最主要的应答细胞内氧气浓度缩小的转录因子成员之一， 不错对多种基因进行调控[23]。海洋软体动物中， 依然有太平洋牡蛎[14]和杂色鲍[24]等HIF-1α基因被进行克隆和抒发研究。本研究通过PCR和RACE时间克隆获取了一种HIF-1α基因cDNA全长序列。序列和结构分析标明， SbHIF-1α具有与其他物种HIF-1α相似性较高的基因序列， 且领有HIF-1α基因保守的HLH、bHLH、PAS-A、PAS-B和PAC结构域， 这阐述SbHIF-1α基因在进化上保留了比较保守的结构域和功能位点[25]。不外， SbHIF-1α与其他物种的HIF-1α比较也存在一定的相反。比如， 东谈主等高档动物以及节肢动物门中小长臂虾的HIF-1α中齐具有两个保守的脯氨酸羟化酶勾搭位点在[26-27]， 而SbHIF-1α只含有一个脯氨酸羟化酶勾搭位点LxxLAP， 但这与太平洋牡蛎[14]和杂色鲍[24]等近缘物种一致， 这可能是软体动物在HIF基因进化上形成的私有特色。魁蚶和栉孔扇贝HIF-1α两个基因的氨基酸序列相似度为90%， 序列相反不赫然， 但是两种贝类的耐低氧智商相反赫然[28]。据报谈， 从线虫到东谈主类， 低氧信号齐启动结构相通或相似的同源基因转录和调控， 并激活雷同的生理、生化反应[29]。是以形成不同贝类对氧气不同稳妥智商的原因可能是发生在HIF-1信号通路效应基因的永别上， 这一效果与鱼类不同物种HIF-1基因相似而耐低氧智商差距大的效果相似[30]。

本研究中， 在血淋巴、外衣膜、鳃、斧足、肝胰腺和闭壳肌共6个组织中均检测到了SbHIF-1α基因转录本， 阐述了HIF散布的凡俗性， 这一效果也与HIF-1α在好意思洲牡蛎和蓝蟹(Callinectes sapidus)等水生无脊椎动物中各组织的抒发散布一致[31-32]。SbHIF-1α在闭壳肌中抒发量最低， SbHIF-1α在血淋巴中抒发量最高， 是闭壳肌的11.45倍， 其他组织的相对抒发量也齐比较低。Cai等[33]发现杂色鲍HIF-1α基因亦然在血淋巴和鳃中抒发量较高， 其他组织相对较低， 与本研究效果较为一致。不外， 好意思洲牡蛎、凡纳滨对虾(Litopenaeus vannamei)以及杂色鲍等其他水生无脊椎动物中HIF-1α均是在鳃组织中抒发量最高[31， 33-34]， 推测原因可能是魁蚶血淋巴比鳃在对外界刺激的明慧方面阐述更首要的作用， 赵庆[19]和黄永欢等[35]研究也显现魁蚶的血淋巴在搪塞外界刺激应答中阐述较积极的作用。

伦理片a在线线2

本研究丰富了不同浓度低氧威胁下， 各个组织HIF-1α基因的抒发谱。魁蚶受到低氧威胁时， HIF-1α在各组织中均作念出实时反应。血淋巴受较强低氧威胁后SbHIF-1α的抒发量变化总体趋势是迟缓升高， 且溶氧越低的试验组对威胁反应更积极， 高潮趋势较更为赫然， 0.5 mg/L试验组在4 h时就达到极权贵水平(P < 0.01)。这阐述魁蚶HIF-1α在低氧威胁初期到低氧威胁的通盘这个词进程齐在阐述首要的作用， 这一丝与Cai等[33]研究效果一致， 而本研究更全面地分析了不同低氧条目下及更长经管时候下HIF-1α的抒发握法。本研究中鳃、外衣膜、斧足、闭壳肌和肝胰腺组织低氧威胁使各经管组 SbHIF-1α的抒发量变化的总体趋势基本是升高后缩小再高潮的趋势， 这一效果与依然研究的太平洋牡蛎、凡纳滨对以及杂色鲍等水活泼物的计议研究效果有一定永别[31， 33]， 这阐述与其他其他海洋软体动物比较， 魁蚶可能具有一套尽头的耐低氧稳妥机制， 推测与其含有红细胞计议， 计议机理还需进一步研究。导致这一变化趋势的原因可能是机体受到低氧刺激后， SbHIF-1α积极反应后与HIF-1β基因相勾搭进一门径控下贱基因， 从而抒发量先减少。但是由于更高强度的低氧威胁可能使机体需要更多的稳妥调控来搪塞， 是以SbHIF-1α基因在各组织中的抒发量会继续高潮以对机体进行更多方面的稳妥性调控。上述变化趋势与魁蚶铁卵白基因在受到外界刺激后的变化趋势一致[36]， 推测这两个基因受到外界环境威胁后的蜕变机制相似。本研究开展了SbHIF-1α在外衣膜、斧足、闭壳肌和肝胰腺组织在低氧威胁下转录水平的研究， 卵白质水平的反应握法还有待进一步进行。

现在， 对于HIF-1的研究多聚拢在哺乳类上， 而贝类HIF-1的研究相对较少。而HIF-1算作机体蜕变氧均衡的主要转录因子， 在呼吸蜕变和搪塞外界刺激的免疫明慧方面阐述作用， 仍有许多后续真切研究亟待开展。本研究初步揭示了魁蚶HIF-1α基因的结构特征和抒发时势黄胖系列， 为魁蚶HIF-1信号通路和低氧稳妥机制研究提供了参考。