中华绒螯蟹（Eriocheir sinensis）Peroxiredoxin 6和Thioredoxin1基因的克隆及表达

中华绒螯蟹是我国重要的水产经济动物，近年来养殖规模不断扩大，产量持续增加。但是，伴随着养殖规模的扩大，养殖环境也日益恶化并导致了大量疾病的发生，严重制约了中华绒螯蟹养殖业的健康发展。因此，疾病预防和控制对中华绒螯蟹养殖业的可持续发展具有举足轻重的作用。与其他无脊椎动物一样，中华绒螯蟹的免疫系统没有免疫球蛋白和淋巴细胞，而是依靠由细胞免疫和体液免疫构成的固有免疫系统来对病原进行识别和清除。中华绒螯蟹的固有免疫机制的研究有助于推动中华绒螯蟹病害防治工作的开展。 本研究采用大规模EST测序方法，结合末端快速扩增（rapid amplification of cDNA ends，RACE）技术从中华绒螯蟹血淋巴中克隆到了过氧化物还原酶（peroxiredoxin，EsPrx6）和硫氧还蛋白（thioredoxin，EsTrx1）基因的cDNA 全长序列；采用实时荧光定量PCR 技术检测了这两个基因在健康个体中表达的组织分布情况以及鳗弧菌刺激后血淋巴细胞中的时序表达规律；同时，将这两个基因的编码区克隆到pET 系列载体，并在大肠杆菌中实现了重组表达，并进行了体外活性检测。 过氧化物还原酶是一个抗氧化蛋白超家族，在保护机体免受活性氧（reactive oxygen species，ROS）的伤害中发挥着重要作用。中华绒螯蟹Prx6（EsPrx6） 基因的cDNA 全长为1076 bp，5` UTR（untranslated region，UTR） 为69 bp，3` UTR 为347 bp，开放阅读框（open reading frame，ORF）为660 bp，编码219 个氨基酸的蛋白。mRNA 3`-端具有多聚腺苷酸加尾信号（polyadenylation signal）AATAAA 和polyA 尾巴。EsPrx6 的预测分子量为 24 kDa，理论等电点为6.21，具有一个保守的Prx 结构域、一个AhpC 结构域和过氧化物酶催化活性中心PVCTTE，表明EsPrx6 属于1-Cys 型Prx。在所检测的组织中均有EsPrx6 的表达，其中以肝胰腺表达量最高，为血淋巴细胞中表达量的17.4 倍。鳗弧菌刺激后，血淋巴细胞中EsPrx6 的表达下降，到12 h 时，实验组显著低于对照组（P<0.05）；随时间推移，表达水平逐渐回升，但在整个实验期间，都没有恢复到起始水平。将EsPrx6 进行体外重组并在大肠杆菌E. coli BL21（DE3）中实现表达，重组EsPrx6 具有预期的抗氧化活性和过氧化物酶活性，其中抗氧化活力为14.69 U/mg 蛋白，高于相同条件下GSH 的抗氧化力（P<0.05），过氧化物酶活力为23.46 U/mg 蛋白。结果表明，EsPrx6 作为一种重要的抗氧化剂，在中华绒螯蟹抵御ROS 可能引起的氧化损伤方面具有重要作用。 硫氧还蛋白是广泛存在于生物体内的一种具有硫醇依赖性的具有还原活性的蛋白。中华绒螯蟹Trx1（EsTrx1）基因的cDNA 全长为641 bp，5` UTR 为17 bp，3` UTR 为306 bp，开放阅读框为318 bp，编码105 个氨基酸。EsTrx1 的预测分子量为12.2 kDa，理论等电点为4.8。EsTrx1 不含信号肽，其氨基酸序列与其他动物的Trx1s 具有高度相似性，如与地中海黄蝎的Trx1 相似度达到73%；而与其他物种Trx2 的同源性很低，相似度仅为14.3-22.8%，表明EsTrx1 属于Trx1 亚族。实时荧光定量PCR 检测发现，EsTrx1 在鳃、性腺、肝胰腺、肌肉、心脏和血淋巴细胞中都有表达。血淋巴细胞中EsTrx1 mRNA 的表达量在菌刺激后上升，刺激后6 h，实验组表达量显著高于对照组和空白组（P<0.05），然后逐渐恢复到刺激前水平。为进一步探讨其生物学功能，将EsTrx1 进行体外重组并在大肠杆菌E. coli BL21（DE3）得到表达，重组EsTrx1 具有预期的氧化还原调节活性，抗氧化活力为3.06 U/mg，且抗氧化活力高于GSH（P<0.05）。rEsTrx1 的二硫键还原活力为5.03，低于凡纳滨对虾的二硫键还原活力（10.44），接近于大肠杆菌（4.93），小牛胸腺（6.50）和小牛肝脏（5.09），而高于鲍鱼Trx2（1.83）活力。结果表明，EsTrx1 在生理条件下能够作为一种重要的抗氧化剂，参与对细菌感染的免疫应答反应。

栉孔扇贝与海湾扇贝免疫学比较研究

扇贝是我国重要的海水养殖品种，但自1997年以来，养殖扇贝陆续爆发的大规模死亡，不但造成了巨大的经济损失，而且直接威胁到现有产业的生存和发展。引起扇贝大规模死亡原因是多方面的，其主要原因是病原滋生、养殖环境恶化和种质衰退。因此，深入研究扇贝免疫防御机制，探讨提高机体抗病力的有效途径和方法，改良种质和培育抗病品系，促进我国扇贝养殖业的健康、可持续发展。本文通过比较栉孔扇贝和海湾扇贝正常状态和重金属污染胁迫下相关免疫指标的变化及HSP70和HSP90在重金属胁迫后表达规律，以期更好的了解贝类的防御机制，为扇贝病害防治提供资料。 本研究采用流式细胞仪技术对栉孔扇贝和海湾扇贝血细胞死亡率、细胞吞噬率和呼吸爆发进行了测定，发现正常状态下两种扇贝的细胞死亡率相差不大，呼吸爆发的基础值相差也不大，而海湾扇贝细胞吞噬率显著高于栉孔扇贝。在血淋巴和肝胰腺中，海湾扇贝SOD、ACP酶活性均显著高于栉孔扇贝；海湾扇贝的MDA含量也高于栉孔扇贝，但差异不显著。鳗弧菌感染实验发现，海湾扇贝累积死亡率远低于栉孔扇贝。上述结果表明，海湾扇贝对细菌等异物的吞噬和杀灭能力以及机体自身的抗氧化能力高于栉孔扇贝，这为海湾扇贝比栉孔扇贝具有更高的抗逆性提供了证据。 不同浓度Pb2+刺激后，海湾扇贝血细胞的死亡率明显低于栉孔扇贝。海湾扇贝细胞吞噬率和呼吸爆发均高于栉孔扇贝。对两种扇贝体液免疫指标的测定发现，海湾扇贝的SOD、ACP活性和MDA含量均高于栉孔扇贝。结果表明，在Pb2+刺激下两种扇贝都处于重金属氧化胁迫下，其免疫系统均受到了影响，但相同剂量的Pb2+对两种扇贝的毒害程度不同，海湾扇贝受到的影响要低于栉孔扇贝。 采用同源克隆和cDNA 末端快速扩增（RACE）技术，从海湾扇贝血淋巴中克隆到了热休克蛋白90（AiHSP90）基因。AiHSP90的cDNA序列全长为2669bp，其中5' 非编码区（UTR）为90bp，3' UTR为44bp，开放阅读框（Open Reading Frame, ORF）包括2175bp，2524-2528位AATAA为推测的AiHSP90基因的 mRNA 多聚腺苷酸加尾信号（Polyadenylation Signal）。开放阅读框编码724个氨基酸残基，推测的理论分子量为83.08 kDa，等电点为4.81。 利用荧光实时定量PCR 技术检测了AiHSP90基因在革兰氏阴性菌鳗弧菌和革兰氏阳性菌藤黄微球菌感染后不同时间点相对于β-actin基因的表达。在鳗弧菌和藤黄微球菌刺激后，除72h外AiHSP90均上调表达，在9h其mRNA表达达到最大值，分别约为空白组的6.8倍和3倍；随后随着时间的推移，其表达下降，在48h恢复到接近空白组的水平；到72h时降到最低点。鳗弧菌对AiHSP90的诱导比藤黄微球菌要高，这也反映出鳗弧菌对海湾扇贝有更强的毒力，诱导的HSP90表达也更强。本研究首次揭示出在软体动物中HSP90在细菌感染条件下的表达情况，结果表明海湾扇贝AiHSP90参与了机体对细菌感染的应答。 采用荧光实时定量PCR技术，对栉孔扇贝和海湾扇贝在不同浓度、不同种类重金属胁迫后HSP70和HSP90基因表达规律进行了比较。结果显示，镉和铅对栉孔扇贝处理10天后即诱导HSP70以较高水平表达，20天后镉处理组HSP70恢复到与对照接近的水平，而在铅处理组也呈现出下降的趋势。而海湾扇贝HSP70基因表达与栉孔扇贝恰好相反，处理20天后其HSP70表达才出现较高水平。铜离子处理对诱导两种扇贝HSP70的影响相似，对栉孔扇贝HSP70的诱导强度更大。对于同一种重金属，不同的浓度诱导的HSP70表达也各不相同。而对于HSP90，海湾扇贝在受到三种重金属刺激后HSP90的表达量均高于栉孔扇贝，而且长时间的刺激（20天）诱导了更高的HSP90基因表达。不同的浓度的同一种重金属离子诱导的HSP90表达也各不相同。这些结果表明，在重金属胁迫下HSP70以及HSP90在栉孔扇贝和海湾扇贝中均表现出不同的应答方式。结合前面海湾扇贝抗逆性较强的综合结果，两种扇贝热休克蛋白基因表达的比较提示我们栉孔扇贝HSP70的激烈反应预示着抗逆能力的不足，而海湾扇贝HSP90的高表达却预示着更强的抗逆能力，有待于更多应激实验如热刺激、病原刺激等的进一步验证。

扇贝大防御素和G型溶菌酶的基因克隆与重组表达

扇贝养殖是我国传统的海水养殖产业，但自1997 年以来，养殖扇贝陆续爆发的大规模死亡，不但造成了巨大的经济损失，而且严重影响了该产业的健康发展。扇贝病害的不断爆发以及病因的多样性迫切要求制定新的疾病防治措施和开发新型的抗菌物质。 从扇贝自身的免疫防御因子入手，筛选和克隆参与免疫防御的功能基因，一方面可以研究抗病功能基因在病原感染或环境胁迫条件下的表达规律，深入探讨扇贝的免疫防御机制，并可作为抗病良种选育的分子标记，指导扇贝的遗传改良和抗病品系的培育；另一方面，可对抗菌效应物实现重组表达，开发新型的病害预防治疗制剂，取代目前普遍使用的抗生素和化学药物。抗菌效应物是机体在免疫应答过程中产生的多肽类物质，对侵入生物体内的细菌、病毒具有很强的免疫杀灭作用，对抗菌效应物的研究有助于深入了解机体先天性免疫防御的机制。 本研究采用大规模EST测序方法，结合cDNA末端快速扩增（RACE）技术，从海湾扇贝血淋巴中克隆到了大防御素基因（big defensin, AiBD）的全长cDNA序列，该cDNA全长为531 bp，其中5' 非编码区（UTR）为24 bp，开放阅读框（Open Reading Frame, ORF）含有369 bp，编码122 个氨基酸残基；随后为138-bp 的3' UTR，包括一个多聚腺苷酸信号序列（AATAAA）和ploy A尾巴。分析表明，海湾扇贝大防御素是以前体的形式合成，前体分子包括信号肽、前域和成熟肽三部分。采用Northern blot方法，以DIG标记的DNA探针检测了 AiBD mRNA在不同组织中的表达。结果发现，AiBD 基因的转录本主要在血淋巴中表达，在鳃中也有微量的表达，而在外套膜、闭壳肌、性腺及肝胰腺中检测不到杂交信号。采用QRT-PCR（quantitative real time PCR）对鳗弧菌感染后海湾扇贝血淋巴中AiBD mRNA 的表达量进行了检测，结果发现在感染后8 h 内， AiBD mRNA 的相对表达量平缓升高；随着刺激时间的增长，AiBD基因的mRNA表达量急剧增加，在刺激后16 h 和32 h 分别达到了空白组的72.3 倍和131.1 倍。为了研究海湾扇贝大防御素的抗菌活性，将其成熟肽编码区克隆到毕赤酵母表达载体pPIC9K并实现了重组表达。抑菌实验表明，重组AiBD具有广谱的抗菌活性，其对供试的三株革兰氏阳性菌（藤黄微球菌、溶壁微球菌、金黄色葡萄球菌）都表现出显著的抗菌活性，而对革兰氏阴性菌（鳗弧菌、亮弧菌）的抑菌活性则相对较弱；此外，重组AiBD对表达宿主也表现出杀菌活性，证明其具有抗真菌活性。 根据栉孔扇贝G 型溶菌酶基因的cDNA序列，利用构建的Genome Walking 文库获得了栉孔扇贝G 型溶菌酶基因的全长序列，该基因序列全长为8131 bp，由六个外显子和五个内含子组成。六个外显子长度分别为55 bp，60 bp，90 bp，113 bp，145 bp 和140 bp；五个内含子的长度分别为1126 bp，2161 bp，2744 bp，750 bp和592 bp；内含子的两侧都具有RNA正确剪接所必需的识别位点（GT/AG）。利用TRANSFAC 软件对栉孔扇贝G 型溶菌酶基因的5' 侧翼序列分析发现，该基因的5' 侧翼具有 TATA box 和 CAAT box 的共有序列；此外，在该基因的5' 侧翼发现了C/EBP、NF-κB、OCT-1 和 NF-IL6 等参与免疫基因激活的转录因子潜在结合位点。采用Northern blot方法，以生物素标记的RNA 探针检测了栉孔扇贝G 型溶菌酶基因在不同组织中的表达。结果发现，该基因的转录本主要在鳃、性腺及肝胰腺中表达，在血细胞和外套膜中也有微量的的表达，而在闭壳肌中检测不到杂交信号，这表明栉孔扇贝G 型溶菌酶可能兼备参与机体免疫防御和消化的功能。为了研究栉孔扇贝G 型溶菌酶的抗菌活性，将其成熟肽编码区克隆到毕赤酵母表达载体pPIC9K并实现了重组表达。抑菌实验表明，重组产物具有显著的抗阳性菌活性，其对供试的藤黄微球菌、溶壁微球菌表现出明显的抑制作用，对金黄色葡萄球菌未检测到抑制活性；而对革兰氏阴性菌仅表现出微弱的抑菌活性（亮弧菌和鳗弧菌），对大肠杆菌则基本无抑制活性。

皱纹盘鲍生长和抗逆相关新基因的初步研究

在刘晓研究员已构建的皱纹盘鲍全长cDNA文库及克隆随机测序的结果基础上，分析发现了皱纹盘鲍胸腺肽β（Tβ）cDNA全长，该序列全长487bp，生物信息学分析显示其含有一个135bp的开放阅读框（ORF），编码一个含44个氨基酸的蛋白质分子，另外其5'和3'非翻译区长度分别是69bp和283bp，3'端非翻译区具有AATAAA mRNA加尾信号及PolyA寡聚腺苷酸尾巴。采用荧光定量PCR的方法研究皱纹盘鲍Tβ的组织特异性发现，Tβ在各组织中均有分布，但在外套膜中高表达，推测和其组织结构中富含参与免疫分泌的颗粒有关。利用荧光定量PCR检测感染鳗弧菌后皱纹盘鲍腹足肌中Tβ在不同时间点的转录表达发现，Tβ基因表现出具有显著差异的时间依赖性，12h后Tβ基因的表达量明显升高，且在处理后24h达到最高值，然后表达量开始回落，至处理后72h，Tβ表达量基本降至初始值。推测皱纹盘鲍在受到病原菌侵染后Tβ基因表达量上调可能是由于Tβ基因参与了皱纹盘鲍的免疫反应，它可能在皱纹盘鲍的抗逆机制中发挥着重要作用。 发现了一个在皱纹盘鲍大小个体中差异表达的长1555bp含有完整3’端的新基因片断，暂命名为TH。利用荧光定量PCR检测皱纹盘鲍RwRh、J1Rh、J1Jm、RwJm 4个家系大小个体中TH基因的差异表达，发现其在大个体中的表达量显著高于小个体，推测其可能与皱纹盘鲍生长速度有某种关联。荧光定量PCR对皱纹盘鲍TH基因时空差异表达的研究发现，其在发育过程中的表达呈现出一定的规律性，即在围口壳幼虫期前始终保持较低水平的表达，当围口壳形成后表达量急剧增加，至幼鲍时达到最高值，在成鲍中也维持较高水平的表达，提示该基因可能调控壳的发育。同时，皱纹盘鲍TH基因具有组织特异性表达，在外套膜中的表达量显著高于其他组织，推测和其组织结构中富含参与壳形成的颗粒有关。另外还比较了皱纹盘鲍RwRh、J1Rh、J1Jm、RwJm 4个家系中TH基因的差异表达，结果发现在这4个家系中的相对表达呈现出一定的规律，即J1Rh＞RwRh＞RwJm＞J1Jm，可能有助于解释家系之间生物学性状上的差异。

栉孔扇贝CfLITAF基因的克隆与表达分析

关于无脊椎动物中是否存在细胞因子一直都存在着异议。从细胞因子本身入手，运用细胞生物学和免疫学手段，已经在软体动物、节肢动物和棘皮动物等无脊椎动物中证实了高等动物细胞因子样活性物质的存在，然而利用分子生物学手段至今没有得到任何核苷酸或蛋白序列信息，而从与细胞因子相关的分子如调控其表达的转录因子入手来研究无脊椎动物中的细胞因子样物质的存在，或许会得到意想不到的惊喜。 最近从人的单核细胞系THP-1中分离纯化出一种新型的转录因子，它可以被LPS诱导表达，产生的蛋白在TNF-α的表达过程中起着非常重要的作用，因而命名为LPS-induced TNF-α factor（LITAF）。随后该基因在小鼠和鸡中也被克隆分离出来，其编码的蛋白已经证实在TNF-α的表达过程中起着不可或缺的作用。迄今为止，在无脊椎动物中还没有相关同源基因的报道。 本研究采用大规模EST测序方法，结合锚定PCR技术从栉孔扇贝体内克隆到了高等动物LITAF基因的同源基因CfLITAF的全长cDNA序列。该cDNA全长为1240bp，其中5' 非编码区（UTR）为112 bp；开放阅读框（Open Reading Frame, ORF）包括450 bp，编码149个氨基酸残基，该多肽的理论分子量为16.08 kDa，等电点为6.77；3' UTR为678bp，包含一个poly A尾巴。SMART结构域预测发现CfLITAF蛋白含有一个保守的LITAF结构域。实时定量PCR检测CfLITAF基因在栉孔扇贝主要免疫器官中的表达分布情况，发现在所检测的六种组织血细胞、外套膜、肌肉、心、腮、性腺中均能检测到CfLITAF基因转录本的不同丰度的表达。血淋巴和肌肉中的表达量相差不大，而在心、外套膜、腮和性腺中的表达量要高于血淋巴和肌肉中的表达量，分别是血中表达量的2.45、2.82、9.23和24.93倍。 为了验证其表达量是否会受到LPS的诱导，利用QRT-PCR（quantitative real time PCR）检测了CfLITAF基因在受到LPS和PGN刺激后在血细胞中的表达规律。结果如下：在LPS刺激后3h，CfLITAF的表达量显著上调，达到了原初水平的1.5倍，随着时间的延长，其表达量逐渐恢复到刺激前水平；在用PGN刺激血细胞后的9个小时内均没有检测到CfLITAF基因mRNA表达量的显著变化。实验中所用血细胞为同一批原代培养的栉孔扇贝血淋巴细胞，细胞活力通过台盼蓝拒染实验测定。 CfLITAF基因的克隆与表达分析，不仅为进一步认识栉孔扇贝的免疫防御机制提供了实验参考，更重要的是为无脊椎动物中细胞因子样物质的研究提供了一个分子水平上的线索，为无脊椎动物中细胞因子样物质的研究奠定了理论基础。

栉孔扇贝抗氧化酶基因的克隆与表达分析

扇贝是我国海水养殖的重要品种，但自1994年以来，养殖扇贝陆续爆发的大规模死亡，不但造成了巨大的经济损失，而且直接威胁到现有产业的生存和发展。扇贝病害的不断爆发以及病因的多样性迫切要求制定新的疾病防治措施和开发新型的抗菌物质。因此，深入研究扇贝免疫防御机制，探讨提高机体抗病力的有效途径和方法，改良种质和培育抗病品系，无疑是解决目前困扰扇贝养殖业健康可持续发展的必经之路。 抗氧化酶可以清除活性氧，是维持机体内氧环境平衡，抵抗外界环境影响的重要免疫因子。本研究采用大规模 EST 测序方法和同源克隆的方法，结合 cDNA 末端快速扩增（RACE）技术，从栉孔扇贝中克隆到了超氧化物岐化酶(Superoxide dismutase，SOD)、过氧化氢酶(Catalase，CAT)、谷胱甘肽过氧化物酶（Glutathione peroxidase，GPX）等抗氧化酶基因的全长 cDNA序列，并对其基因结构进行了分析。同时，用实时定量PCR方法对这三个基因在健康扇贝血淋巴细胞、肾、鳃、肌肉、性腺等组织和在分别用鳗弧菌，溶壁微球菌和假丝酵母处理扇贝后不同时间段的表达差异情况进行了研究。 超氧化物歧化酶基因CfSOD的cDNA 全长为1022 bp，其中开放阅读框（Open Reading Frame, ORF）含有 459 bp，编码 153个氨基酸残基，无信号肽，为胞内蛋白。经BLASTP分析发现，CfSOD与其它动物具有较高的同源性。CfSOD中存在两个Cu/Zn-SOD的签名序列；另外Cu结合必须氨基酸(His-45，-47，-62 和-119)和Zn结合必须氨基酸(His-62，-70，-79和Asp-82)在CfSOD中保守。实时定量PCR 检测发现，CfSOD在鳃、血细胞和肾中有较高的表达。在鳗弧菌和溶壁微球菌刺激后，CfSOD的相对表达量逐渐下降，然后分别在32小时和16小时的时候恢复到刺激前的表达水平。在假丝酵母刺激后，CfSOD的mRNA表达没有显著差异。 栉孔扇贝过氧化氢酶基因CfCAT的cDNA全长为3146 bp，其中开放阅读框含有1521bp，编码507个氨基酸残基，无信号肽，为胞内蛋白。经BLASTP分析发现，CfCAT与其它动物具有较高的同源性。 CfCAT中存在过氧化氢酶近端活性位点和过氧化氢酶近端血红素配体签名序列，另外存在两个糖基化位点 NFS和 NFT，同时在CfCAT 的C末端存在过氧化物酶体定位信号AQL，为典型过氧化氢酶。实时定量PCR 检测发现，健康的扇贝中CfCAT在鳃和性腺中有较高的表达。CfCAT基因在鳗弧菌刺后表达升高，在4小时达到最高，约是刺激前表达量的6.8倍（P<0.05），后随着时间的变化而逐渐下降。在8小时表达量达到为刺激前表达量的1.3倍（P<0.05），在16和32小时略高于刺激前的水平。在溶壁微球菌刺激后CfCAT基因表达量也呈上升趋势，在刺激后4小时达到刺激前表达量的约2倍，然后有所下降，在16 小时又上升到刺激前表达量的2.9倍。CfCAT基因在假丝酵母刺激后的表达略有升高，4小时约是刺激前的1.2倍（P<0.05），在其他时间段变化不明显。 栉孔扇贝谷胱甘肽过氧化物酶基因CfGPX的cDNA 全长为1290 bp，其中开放阅读框含有705bp，编码235个氨基酸残基，有一个24核苷酸的信号肽序列。经BLASTP 分析发现，CfSOD与其它动物具有较高的同源性。CfGPX中发现谷胱甘肽过氧化物酶活性位点的签名序列， 另外发现硒半胱氨酸和硒半胱氨酸插入序列，为含硒型谷胱甘肽过氧化物酶。实时定量PCR 检测发现，未经处理的扇贝中CfGPX在性腺、肌肉、血和肾中有较高的表达。CfGPX基因在鳗弧菌刺后表达量快速上升，在6 小时的时候表达量达到最高，为刺激前的4.0倍（P>0.05），后随着时间的变化而逐渐下降。在溶壁微球菌刺激后CfGPX在前6小时表达略有降低，在6小时的时候表达量为刺激前的0.5倍（P<0.05），后随着时间的变化而逐渐升高。在16小时的时候表达量为刺激前的2.1倍（P<0.05）。在假丝酵母刺激后， CfGPX的表达量略有下降在8小时的时候表达量为刺激前的0.8倍（P<0.05）。 实验证明栉孔扇贝的超氧化物歧化酶基因CfSOD，过氧化氢酶基因CfCAT，谷胱甘肽过氧化物酶基因CfGPX基因在机体抵抗外界微生物刺激中起到了重要的作用。

扇贝丝氨酸蛋白酶抑制剂基因CfKZSPI及Aikunitz的克隆与重组表达

扇贝养殖是我国重要的海水养殖产业，然而自1997 年以来，养殖扇贝陆续爆发的大规模死亡，不但造成了巨大的经济损失，而且严重影响了该产业的健康发展。丝氨酸蛋白酶抑制因子及丝氨酸蛋白酶在无脊椎动物的免疫应答中起着核心作用，它们的协同作用直接导致外界病源入侵的信号转导和级联放大，并进一步激活一系列防御体系，如黑化反应、血液凝结和抗菌肽的合成等。因此，克隆扇贝参与免疫防御的丝氨酸蛋白酶抑制剂基因并对其功能进行研究，将有助于进一步研究扇贝的免疫防御机制，丰富和发展无脊椎动物免疫学的内容。 运用大规模EST技术和RACE技术从栉孔扇贝中克隆出一个Kazal型丝氨酸蛋白酶抑制剂基因，定名为CfKZSPI。该基因cDNA序列全长1788bp，其中5' 非编码区（Untranslated Region, UTR）为97 bp，3' UTR161 bp，有一个典型的多聚腺苷酸信号序列（AATAAA）和一个ploy A 尾巴，开放阅读框（Open Reading Frame, ORF）含有1530 bp，编码509 个氨基酸残基。对其推测氨基酸序列进行分析，发现其中包括22个氨基酸残基组成的信号肽序列和12个Kazal型丝氨酸蛋白酶抑制剂结构域。采用QRT-PCR（quantitative real time PCR）对鳗弧菌浸泡刺激后栉孔扇贝血淋巴中CfKZSPI 的 mRNA表达量进行了检测，发现其mRNA 的表达量在鳗弧菌刺激后3h明显上升，达到空白组的43.6倍；然后在6h时有所下降，为空白组的15.0倍；随着菌刺激时间的增长，CfKZSPI基因的 mRNA 表达量急剧增加，在刺激后8h，12h，24h分别达到空白组的174.1，207.8，675.4倍。统计分析发现3h（P=0.019<0.05）和12h（P=0.020<0.05）时，CfKZSPI基因mRNA表达量与空白组差异均显著。为了研究栉孔扇贝CfKZSPI的蛋白活性，将其第十二个结构域克隆到pET-32a（+）载体中，转化大肠杆菌Rosetta-gami（DE3）表达菌株，获得可溶性表达的蛋白rCfKZSPI-12，对其进行抑制蛋白酶活性的分析，发现其对胰蛋白酶有很强的抑制活性，而对凝血酶没有抑制活性。当rCfKZSPI-12与胰蛋白酶分子比率为1：1时，约90%的蛋白酶活性被抑制。运用狄更斯作图法研究rCfKZSPI-12对胰蛋白酶的抑制能力，结果发现其对胰蛋白酶的抑制常数为173 nmol L-1。 采用同样方法从海湾扇贝cDNA文库中克隆出一个Kunitz型丝氨酸蛋白酶抑制剂基因，定名为Aikunitz。该基因全长632 bp，其中5' UTR 为105 bp，3' UTR 为 245 bp，有一个典型的多聚腺苷酸信号序列（AATAAA）和一个ploy A 尾巴，ORF 含有282 bp，编码93 个氨基酸残基。推测的氨基酸序列N末端有一个20个氨基酸残基组成的信号肽序列，成熟蛋白包括一个Kunitz型丝氨酸蛋白酶抑制剂结构域。采用QRT-PCR对鳗弧菌和藤黄微球菌感染后海湾扇贝血淋巴中Aikunitz 的mRNA的表达量进行了检测，结果发现其在鳗弧菌刺激后3h到9h持续上升，9h时表达量为PBS对照组的4.49倍（P=0.008<0.05），然后开始下降，在72h时表达量为对照组的0.24倍（P=0.021<0.05）；而在藤黄微球菌刺激后3h到12h其表达量上升，其中6h时为空白组的5.95倍（P=0.0004<0.01）；12h以后迅速下降，其中24h的表达量为对照组的0.38倍（P=0.028<0.05）。将Aikunitz基因编码的成熟蛋白按照重组CfKZSPI-12的方法进行重组表达，并对重组蛋白进行抑制蛋白酶和抑菌活性分析。结果发现其对胰蛋白酶和弹性蛋白酶两种丝氨酸蛋白酶都没有抑制作用。抑菌实验同样发现，重组Aikunitz 对供试的革兰氏阳性菌藤黄微球菌和革兰氏阴性菌鳗弧菌和大肠杆菌都不显示明显抑菌活性。

迟缓爱德华氏菌基因组fosmid文库的构建、opp基因簇的克隆及aroA基因缺失突变株的构建

本研究构建了迟缓爱德华氏菌(Edwardsiella tarda)LSE40 基因组fosmid文库，该文库共包含2 500个克隆，插入片段平均大小为33.6kb，文库总容量约84Mb，覆盖E.tarda LSE40基因组（按5Mb计算）超过16倍。随机挑取1 000个fosmid克隆进行双末端测序共得到1 741条高质量的序列，序列平均长度546 bp，全长949 997bp，约为E.tarda基因组的19％。将这些序列提交到KEGG自动注释服务器KAAS对所得序列进行代谢途径分析，得到E.tarda LSE40的KO (KEGG Orthology) 注释。分析结果表明，与代谢途径相关的基因有932条序列，与环境信息处理相关基因283条，与遗传信息处理相关基因220条，与细胞进程和人类疾病相关的基因分别为64条和16条。同时将序列进行BlastX，按照微生物致病性共同主题找到61个毒力相关基因。Fosmid文库的建立和部分基因组序列的生物信息学分析为进一步研究E.tarda LSE40的致病机制、代谢机制和生理生态机制提供了丰富的物质基础。 通过比较基因组的方法，从E.tarda LSE40 fosmid文库克隆到编码寡肽透过酶的opp基因簇，该基因簇全长6 741bp，含有5个ORF，依次编码OppA-B-C-D-F 5个蛋白；位于oppA和oppB的间隔区和oppF之后的非编码区各有一个茎环结构，推测分别为oppA和opp基因簇的转录终止子。以细菌OppA的保守结构域SBP_bac_5构建系统发生树，结果显示E.tarda LSE40与同属细菌E.ictaluri的亲缘关系最近，与肠杆菌科细菌的亲缘关系较近，与革兰氏阳性细菌的亲缘关系较远，表明OppA的SBP_bac_5结构域可作为细菌分类鉴定的依据。 从E.tarda LSE40 fosmid文库克隆aroA基因全序列，该序列全长1 287bp，编码428个氨基酸，与鲶鱼爱德华氏菌(E. ictaluri)氨基酸相似性在94%，与其他肠杆菌科菌如Escherichia coli和Yersinia enterocolitica相似性在73%-74%。通过In-frame deletion构建了E.tarda LSE40 aroA缺失突变株。与野生型相比，aroA突变株的半数致死量LD50提高了62倍。在牙鲆接种~106cfu/ml的E.tarda细菌时，接种野生型细菌的牙鲆在6天内全部死亡，濒死鱼的细菌数达7.97×108cfu/ 100mg；而接种aroA突变株的牙鲆没有出现死亡，28天后检测不到细菌的存在。实验结果为进一步评价aroA突变株作为减毒活疫苗打下了基础。

大熊猫RPS15 cDNA的克隆及序列分析

本研究运用RT—PCR技术，首次从大熊描Ailuropoda melanoleucu的肌肉组织总RNA中成功克隆了核糖体蛋白S15（RPS15）基因的表达序列．并对其进行了初步分析。结果表明：大熊猫RPS15基因的表达序列全长为442bp，开放阅读框（ORF）为438bp，编码145个氨基酸，该蛋白的分子量为17．0401KDa，等电点为10．3，含有2个依赖于cAMP和cGMP的蛋白激酶磷酸化位点，5个蛋白激酶C磷酸化位点，4个N-酰基化位点及1个RPS19蛋白signature位点。进一步分析发现，大熊猫RPS15基因的表达序列及其编码的氨基酸序列与已报道的部分哺乳动物具有很高的相似性。

Isolation and characterization of a novel variant of HMW glutenin subunit gene from the St genome of Pseudoroegneria stipifolia

The x- and y-type high molecular weight (HMW) glutenin subunits are conserved seed storage proteins in wheat and related species. Here we describe investigations on the HMW glutenin subunits from several Pseudoroegneria accessions. The electrophoretic mobilities of the HMW glutenin subunits from Pd. stipifolia, Pd tauri and Pd strigosa were much faster than those of orthologous wheat subunits, indicating that their protein size may be smaller than that of wheat subunits. The coding sequence of the Glu-1St1 subunit (encoded by the Pseudoroegneria stipifolia accession PI325181) was isolated, and found to represent the native open reading frame (ORF) by in vitro expression. The deduced amino acid sequence of Glu-1St1 matched with that determined from the native subunit by mass spectrometric analysis. The domain organization in Glu-1St1 showed high similarity with that of typical HMW glutenin subunits. However, Glu-1St1 exhibited several distinct characteristics. First, the length of its repetitive domain was substantially smaller than that of conventional subunits, which explains its much faster electrophoretic mobility in SDS-PAGE. Second, although the N-terminal domain of Glu-1St1 resembled that of y-type subunit, its C-terminal domain was more similar to that of x-type subunit. Third, the N- and C-terminat domains of Glu-1St1 shared conserved features with those of barley D-hordein, but the repeat motifs and the organization of its repetitive domain were more similar to those of HMW glutenin subunits than to D-hordein. We conclude that Glu-1St1 is a novel variant of HMW glutenin subunits. The analysis of Glu-1St1 may provide new insight into the evolution of HMW glutenin subunits in Triticeae species. (C) 2007 Elsevier Ltd. All rights reserved.