不同壁材微胶囊饲料对黄姑鱼稚鱼生长和消化酶活力的影响
Effects of different wall materials of the microencapsulated diets on the growth and digestive enzyme activity of Nibea albiflora Larvae
-
摘要: 研究采用湿法制粒流化床包衣工艺, 分别以明胶、乙基纤维素、玉米醇溶蛋白为壁材制备微胶囊饲料, 比较其对黄姑鱼稚鱼生长和消化酶活力的影响. 粒径(178-590) m的3种微胶囊饲料质量均大于50%. 扫描电镜观察微胶囊饲料的表面均有一层较为致密的包衣薄膜. 壁材明胶、乙基纤维素、玉米醇溶蛋白微胶囊饲料的包含率分别为95.4%、95.6%和95.8%; 脂类包埋率分别为72.6%、76.5%和64.3%; 氮保留率分别为53.5%、62.3%和54.6%. 将3种微胶囊饲料分别饲喂15日龄黄姑鱼稚鱼30d. 明胶组和玉米醇溶蛋白组稚鱼的体重、全长均显著高于乙基纤维素组(P0.05), 但成活率差异不显著(P0.05); 明胶组稚鱼的体重、全长和成活率均高于玉米醇溶蛋白组, 但差异均不显著(P0.05). 明胶组稚鱼的胰蛋白酶活力显著高于乙基纤维素组和玉米醇溶蛋白组(P0.05), 但淀粉酶和碱性磷酸酶活力的差异均不显著(P0.05). 与乙基纤维素、玉米醇溶蛋白相比, 明胶更适合作为黄姑鱼稚鱼微胶囊饲料壁材.Abstract: In this study, we evaluated the effects of different wall materials of microencapsulated diets on the growth and the digestive enzyme activity of N. albiflora larvae. Three groups of diets were microencapsulated with gelatin, ethyl celluse, and zein. The microencapsulated diets were prepared with wet granulation and fluidized bed coating process. More than 50% of the microencapsulated diets had diameters between 250 m and 590 m. Scanning electron microscopy microphotographs showed that the surface of the microencapsulated diet was covered by a dense film. The inclusion efficiency of the gelatin, ethyl celluse, and zein groups was 95.4%, 95.6% and 95.8%, respectively; the lipid encapsulation efficiency was 72.6%, 76.5% and 64.3%, respectively; the nitrogen retention efficiency was 53.5%, 62.3% and 54.6%, respectively. A 30-day feeding experiment was conducted with the larvae at DAH 15. At DAH 20, the larvae were weaned from copepods to the experimental diets. The wet weight and total length of larvae from the gelatin group and the zein group were significantly larger than those from the ethyl cellulose group (P0.05). There was no significant difference in the survival rate between the groups. The trypsin activity of the larvae from the gelatin group was significantly higher (P0.05) than the ethyl cellulose group and the zein group. There were no significant differences in the amylases activity and the alkaline phosphatases activity between the groups. Therefore compared to ethyl cellulose and zein, gelatin should be more suitable for the wall material of the microencapsulated diets for N. albiflora larvae.
-
Keywords:
- Nibea albiflora /
- Larvae /
- Microencapsulated diet /
- Growth /
- Digestive enzyme activity
-
纤毛虫原生动物具有二型核, 小核基因通常被内部删除序列(Internal eliminated sequences, IESs) 分隔成多个大核命运序列(Macronucleus destined sequences, MDSs) 片段。在接合生殖过程中, 通过小核IESs的删除和MDSs的重组连接、端粒序列的生成和基因的高度扩增等过程形成具功能的大核基因[1-5]。染色体重排广泛存在于3种纤毛虫类群中: Armophorea、旋唇纲Spirotrichea、层咽纲Phyllo-pharyngea[6], 许多小核基因MDSs的排列顺序会发生变化, 形成乱序基因。其中在旋唇纲腹毛类纤毛虫中该现象尤为明显, 例如尖毛虫(Oxytricha) 和棘尾虫(Stylonychia) 中分别有96%和98%的小核基因组序列在小核发育成大核的过程中被删除[7]。
在腹毛类纤毛虫中, 有3种乱序基因为人们熟知, 分别为肌动蛋白Ⅰ型(Actin Ⅰ)[1]、端粒结合蛋白(α-TBP)[8]及DNA聚合酶α (DNA polymerase α)[9]。目前已报道10个物种12个种群的小核Actin Ⅰ基因序列[4, 1, 10-14]、6个物种的小核α-TBP基因[8, 15, 16]、7个物种的小核DNApol α基因[9, 17-20]。Hogan等[12]推测了腹毛类纤毛虫小核Actin Ⅰ基因的进化路线, 从共同的原始祖先到3个MDSs的Urostyla grandis, 进化为拥有4个MDSs的Engelmanniella mobilis, MDS2发生倒置, 再进化为拥有8个MDSs的Stylony-chia lemnae, 然后分化为两个分支, 一支为3个尖毛虫类群, 一支为S. pustulata和Oxytricha sp. (Misty)。之后, Chen等[21]在成功扩增2种尾柱目纤毛虫大小核Actin Ⅰ基因的基础上, 将腹毛类小核Actin Ⅰ进化路线进一步修正, 提出腹毛类不同目之间可能存在不同的进化路线。Wong等[16]推测了α-TBP基因可能的进化路线, 认为祖先类群为非乱序物种(如Holosticha sp.), 历经乱序排列的Uroleptus sp.、Pa-raurostyla weissei, 进化为中间MDS发生易位的O. trifallax, 再经历两个MDSs发生融合(如O. nova/S. myti-lus)。DNApol α的进化路线尚未有人给出。综上, 已报道的腹毛类3种乱序基因进化路线仅为较为粗略的大致框架, 细节部分依然亟待修订和补充。
本研究在湛江湖光红树林水域中采集得到一个腹毛类尖毛虫物种Oxytricha sp. (ZJ)。通过测定O. sp. (ZJ) 大、小核的Actin Ⅰ、α-TBP和DNApol α基因序列, 并对尖毛虫属内现有物种3个乱序的小核基因结构进行比较, 探讨尖毛虫属基因乱序模式, 增加人们对相关基因进化的认识。
1. 材料与方法
1.1 样品与DNA提取
Oxytricha sp. (ZJ) 于2010年11月采自湛江湖光红树林水域(东经109°40′-110°35′, 北纬20°14′-21°35′), 镜检后分离[22]。于5‰人工海水中投入米粒滋生细菌, 进行单克隆培养。至大量繁殖阶段, 筛绢过滤除去大的杂质后, 5000 r/min, 3min离心富集虫体。DNA提取按照伊珍珍等[23]改良的Qiagen Dneasy Blood & Tissue Kit试剂盒使用方法操作。电泳检测完整性, 微量核酸测量仪检测其浓度。
1.2 PCR扩增与测定
大核Actin Ⅰ基因序列的获得 利用引物Actin800F和Actin800R1[18](表 1) 获得约800 bp的大核Actin Ⅰ基因片段。PCR反应总体系为50 μL, 包括去离子水36 μL, 10×Ex buffer 5 μL, dNTP mix (2.5 mmol/L) 5 μL, 模板2 μL, 正向反向引物(25 mmol/L) 各0.8 μL, Ex Taq (5 U/μL) 0.4 μL。PCR扩增程序如下: 95℃预变性5min; 95℃变性30s, 52℃退火30s, 72℃延伸1min, 扩增35个循环; 72℃终延伸10min。经1.0%琼脂糖凝胶电泳检测后, 切胶, 用Universal DNA Purification Kit试剂盒回收纯化, 购自天根生化科技(北京) 有限公司, 连接至pMD18-T载体(TaKaRa, 日本), 转化大肠杆菌DH5α感受态细胞, 挑选5个阳性克隆后, 送美吉公司双向测序, 测序引物为M13F-47与M13R-48。
表 1 大、小核Actin Ⅰ基因所用引物名称及序列Table 1. The name and sequence of primers for macro-and micronuclear Actin Ⅰ gene引物名称Name 引物序列Sequence (5′-3′) Actin800F AACTGGGAYGAYATGGARAAGAT Actin800R1 ATCCACATKSHGGCGAAGGT J-F (1203) GGTTATTGCCAGCCCAGACAGAA J-M3F1 TTTCTCAGGTGAGGATGCTCC J-M3F2 TCCCATCAATCGTCGGTAGACCC AP12 GTAATACGACTCTATAGGGCACGCGTGGTCGACGGCCCGGGCTGGTCCCCAAAACCCCAAAACCCCAAAA AP1 GTAATACGACTCACTATAGGGC AP2 ACTATAGGGCACGCGTGGT J-F (673) TGTGATGCTGGTGATGGTGTTAC J-(133RC) CCCTACTAATCTCATTCTTGGCAT J-M9F TGCCCTATCAAAGTATTTCCTA J-F (42) CCATCAGTAAATTTGACATTACTAT J-M8R CTCATTCAACCTGTCCGTAAC J-I1F CATTATGGATTTATGTTTATTGGTG J-I2F CCCAAGTCAATCAACAAAAGTATC J-R (1550) TTGAAACTAATCTGAATAAAAATG 基于上述获得的序列设计引物J-F (1203)、J-M3F1、J-M3F2 (表 1)。利用Chang等[19]建立的端粒抑制PCR (TSP-PCR) 扩增。分别以AP12 (表 1)/Actin800R1、Actin800F/AP12、J-M3F1/AP12为上下游引物进行第一轮PCR, 扩增条件为: 94℃ 2min, 72℃ 4min, 重复7个循环, 再94℃ 2min, 67℃ 4min, 重复32个循环, 然后67℃ 4min。再分别以AP1 (表 1)/Actin800R1、Actin800F/AP1、J-M3F1/AP1为上下游引物做第二轮巢式PCR, AP2 (表 1)/Actin800R1、J-F (1203)/AP2、J-M3F2/AP2为上下游引物做第三轮巢式PCR。扩增条件为: 94℃ 25s, 72℃ 4min, 重复5个循环, 再94℃ 25s, 67℃ 4min, 重复20个循环, 然后67℃ 4min。连续3轮PCR, 获得基因左端序列片段AP2/Actin800R1和右端序列J-F (1203)/AP2以及J-M3F2/AP2。扩增产物经切胶回收、转化、克隆测序, 拼接得到大核基因序列全长, 并提交至GenBank, 登录号为KX602204。
小核Actin Ⅰ基因序列的获得 根据已报道Oxytricha属的小核Actin Ⅰ基因乱序结构, 在推测出的乱序区MDS3/4和MDS1上分别设计上、下游引物, J-F (673)(表 2) 和J-(133RC)(表 2), 使用巢式PCR和降落PCR扩增(95℃预变性5min, 95℃变性30s, 67℃退火30s, 之后退火温度依次降低1℃, 72℃延伸2min, 扩增15个循环; 95℃变性30s, 52℃ 30s, 72℃延伸2min, 扩增30个循环; 72℃终延伸10min)。得到目的条带。切胶回收、转化、克隆测序, 获得小核部分基因片段。
表 2 大、小核α-TBP基因所用引物名称及序列Table 2. The name and sequence of primers for macro-and micronuclear α-TBP gene引物名称Name 引物序列Sequence (5′-3′) αTBP-R GTGTCCTTRATGARRTAGAAWCC J-R2 TTGAAGAAGTTAGCACCGAGACCGT J-R1 GTCGGAGGGTTCAATCTTGGTTACG J-M5F1 TCTACAACAGTTCATGGGCTCTC J-M5F2 CACTTGAGGAAATCAGTGGTGCTG J-M2R TGAGTCACCAGTTCCCTTTTAGC J-M1F1 TTGAGATATGTGGCTGGATTTG J-M1F2 GGCTGGATTTGAAAATACGAGTAG J-I4R GTAGGTAAATTAAAGATTTGGC J-M9F TTAAGGATGGCTCAGGTCAAGT J-M8F1 TTGTTGAGCCCACTTTCTGATG J-M8F2 CGAAAACTTGGAGGATCTTAGCA J-M11F GGTGAGGTAGTAAGAATCAGATC J-M10R GTGGGAATTTAAGTTTAAGAGC J-I9F GAGAAGTCTATAAGTAAAAGTGCC J-R (2133) GATTATAAAATGATTAAAGATCATG J-I13R2 CTTATCATAAACATTATTCAC 在已获得的小核基因序列的基础上设计多条引物J-M9F、J-F (42)、J-M8R、J-I1F、J-I2F、J-R (1550)(表 1)。再分别以J-M9F/J-M8R、J-M3F1/J-(133RC)、J-I1F/J-M8R为上下游引物进行首轮PCR, J-F (42)/J-M8R、J-M3F2/J-R (1550)、J-I2F/J-M8R进行巢式PCR, 扩增得到片段J-F (42)/J-M8R、J-M3F2/J-R (1550)、J-I2F/J-M8R。覆盖小核基因全长。拼接获得小核基因序列全长(图 1), 并提交至GenBank, 登录号为KX602205。
![]()
图 1 Oxytricha sp. (ZJ) 大、小核Actin Ⅰ基因结构、引物位置及获得的PCR产物示意图灰色方块代表MDSs, 黑色方块代表IESs; 其中MDS2倒置; 白色方块代表未得到的基因区域, 下同Figure 1. Schematic illustration of macro-and micronuclear structures, PCR primer locations and products obtained in Oxytricha sp. (ZJ)Actin Ⅰ geneMDSs are represented in grey, IESs are represented in black; MDS2 is inverted; white box represents unavailable gene regim; the same applies below大、小核α-TBP基因序列的获得 根据上文描述的方法测定Oxytricha sp. (ZJ) 的α-TBP基因的完整大核序列。大核扩增所用引物为αTBP-F、αTBP-R、J-R2、J-R1、J-M5F1、J-M5F2 (表 2)。小核所用引物为J-M2R、J-M1F1、J-M1F2、J-I4R、J-M9F、J-M8F1、J-M8F2、J-M11F、J-M10R、J-I9F、J-R (2133)、J-I13R2 (表 2、图 2), 并提交至Gen-Bank, 大核基因登录号为KX602208, 小核基因为KX602209。
![]()
图 2 Oxytricha sp. (ZJ) 大、小核α-TBP基因结构、引物位置及获得的PCR产物示意图Figure 2. Schematic illustration of macro-and micronuclear structures, PCR primer locations and products obtained in Oxytricha sp. (ZJ)α-TBP gene大、小核DNApol α基因序列的获得 基于上文描述的方式获得Oxytricha sp. (ZJ) 的DNApol α基因完整的大核序列。引物为Pol α-F、Pol α-R、J-R2 (DPα)、J-R1 (DPα)、J-F3 (DPα)、J-F4 (DPα) (表 3), 得到DNApol α基因完整的大核序列。
表 3 大、小核DNApol α基因所用引物名称及序列Table 3. The name and sequence of primers for macro-and micronuclear DNApol α gene引物名称Name 引物序列Sequence (5′-3′) Pol α-F CTAYTGGATHGATGCWCATG Pol α-R CACCAVTCTCTWCKNACCATATC J-R2(DPα) CTAGTAACTGGCTTGCATCTCC J-R1(DPα) AGTGCCTCTAGCTTTGCCTTTC J-F3(DPα) CGATGGAGTGTTTAAGAGTTTAC J-F4(DPα) ATCACCAGATGCTAAAGTTGTGC J-M32F CGTTTCTGAACTAAATCTCT J-M1R1 CCTCTCTTCCTGACTATGTTGTA J-M1R2 TTTAGTTCCATCTTTACGCCCAG J-L1F TCCACTCATCATTTTTTTACTGC J-R5R CTAGTAACTGGCTTGCATCTCC J-R6R AATGTGTTTCCAGTGACWGTGCT J-R9R TCTTTTCTCTGTCTGAGTTGCTA J-5000R CATGCTCACTCTCTTAATACCTC J-5100R GCTAAACCCCATGTAGAAAAGTTG 因在DNApol α基因小核序列扩增时, 扩增结果均为大、小核基因序列的嵌合体, 故只获得部分序列。所用引物为J-M32F、J-M1R1、J-M1R2。J-L1F、J-R5R、J-R6R、J-R9R、J-5000R、J-5100R (表 3、图 3), 并提交至GenBank, 大核基因登录号为KX602206, 小核基因为KX602207。
![]()
图 3 Oxytricha sp. (ZJ) 大、小核DNApol α基因结构、引物位置及获得的PCR产物示意图Figure 3. Schematic illustration of macro-and micronuclear structures, PCR primer locations and products obtained in Oxytricha sp. (ZJ) DNApol α gene大核Actin Ⅰ基因系统树的构建 因为已完成大核α-TBP基因和DNApol α基因测序的尖毛虫物种数目太少(仅3种), 故而只依据Actin Ⅰ基因序列构建系统树。除本工作测得序列(序列号为: KX-602204) 外, 另外7条序列源于Genbank数据库(序列号为: U63567、AF508049、AF508050、AY044839、U18940、U63566、M22480)。使用Bioedit7.1.9.0[24]内置Clustal进行多序列比对, 序列末端截平, 并手动调整删除空格或比对模糊序列。本工作使用Mr-Modeltest 2[25]进行核酸模型筛选GTR+G+I, 使用MrBayes3.1.2[26]构建了贝叶斯树(Bayesian infe-rence, BI)。用马尔科夫链的蒙特卡洛方法(MCMC) 设置为4条链, 运行1000000代, 每100代进行抽样, 以后验概率来表示各分支的可信性。TreeView v5.1.2[27]、MEGA7[28]等软件用于树形图的美化。
2. 结果
2.1 Oxytricha sp. (ZJ) 的大、小核Actin Ⅰ基因结构特征
Oxytricha sp. (ZJ) 的大核Actin Ⅰ基因全长为1546 bp (含端粒序列)。含有编码375个氨基酸的开放阅读框(ORF), 其AT含量为53.81%。5′和3′非编码区的序列长度分别为198 bp (含20 bp端粒序列) 和226 bp (含20 bp端粒序列), AT含量分别为73.71%和62.87%。5个阳性菌斑测序累计检测到59个碱基变异位点, 其中ORF含有38个碱基变异位点, 涉及12个氨基酸的差异, 未引起移码突变或无义突变。剩余的碱基变异位点, 有5个位于5′非编码区, 16个位于3′非编码区。推测多克隆间序列差异是由于基因重复现象引起[29]。
Oxytricha sp. (ZJ) 的小核Actin Ⅰ基因序列长度为1775 bp。通过与大核序列比对确定该小核基因由8个IESs和9个MDSs组成。IESs序列长度为11-79 bp, MDSs序列长度为19-583 bp。MDS/IES交接处有5-13 bp的指针重复序列。8个IESs均为富含AT的序列(68.00%-100.00%)(表 4)。MDSs按照3-4-6-5-7-9-2-1-8的顺序排列, 其中MDS2为倒置排列(图 1)。
表 4 Oxytricha sp. (ZJ) 小核Actin Ⅰ基因的MDSs、IESs及指针序列Table 4. Characteristics of MDSs, IESs, and pointer sequences in Oxytricha sp. (ZJ)Actin Ⅰ geneMDSn MDSn的长度Length of MDSn MDSn和(n+1) 间的指针序列Pointer sequence between MDSn and (n+1) MDSi和MDSj间IES的长度Length of IESs between MDSi and MDSj IESs AT含量IESs AT content (%) 1 196 ACAACAAACTTGC 3-4 79 79.75 2 35 ? …GGCTGG 4-6 30 73.33 3 338 TCTCT 6-5 68 83.82 4 583 CCAAGTCAATCAA 5-7 30 80.00 5 19 TATTRCCAGC 7-9 11 100.00 6 139 CTGAGCAATCA 9-2 21 90.48 7 87 TTAACWATTTR 2-1 26 80.77 8 ? AGGTGG…? 1-8 25 68.00 9 99 6-8 注: MDSs和IESs的长度都包括了一拷贝的指针序列; 不同克隆间指针序列的差异用兼并碱基表示; 问号代表数值或序列无法确定
Note: Length of MDSs include pointers on both ends; IESs exclude pointer sequences; Use degenerate symbols reflects differences between clones; Question marks indicate undetermined information基于尖毛虫属内8个物种的大核Actin Ⅰ基因序列建的贝叶斯树(图 4a) 中, Oxytricha trifallax、O. fallax与O. granulifera以最高置信值聚为一支, 其余5个物种O. longa、O. sp. (Misty)、O. nova、O. sp. (Aspen)、O. sp. (ZJ) 聚为一支。
![]()
图 4 基于Actin Ⅰ基因序列构建的贝叶斯树(a) 与Actin Ⅰ基因结构进化路线图(b)Figure 4. Baysian trees based on the sequences of Actin Ⅰ gene (a), and evolutionary route map of Actin Ⅰ gene structure (b)2.2 Oxytricha sp. (ZJ) 的大、小核α-TBP基因结构特征
Oxytricha sp. (ZJ) 的大核α-TBP基因全长为2140 bp。含有编码495个氨基酸的ORF, 其AT含量为53.56%。5′和3′非编码区的序列长度分别为170 bp (含20 bp端粒序列) 和438 bp (含20 bp端粒序列), AT含量分别为74.83%和66.50%。通过与前人数据[22]比较, 找到一个长度为53 bp、富含AT (77.36%) 的内含子序列。5个阳性克隆测序累计检测到107个碱基变异位点, 其中ORF含有49个碱基变异位点, 涉及19个氨基酸的差异, 未引起移码突变或无义突变。剩余的碱基变异位点, 有2个位于5′非编码区, 53个位于3′非编码区, 3个位于内含子。
Oxytricha sp. (ZJ) 小核α-TBP基因序列长度为2562 bp。通过与大核序列比对确定该小核基因由13个IESs和14个MDSs组成。IESs序列长度为16-71 bp, MDSs序列长度为13-731 bp。在MDS/IES交接处有5-11 bp的指针重复序列。13个IESs中, 除MDS1/MDS3间的AT含量较低外(52.38%), 其他12个IESs均为富含AT序列(62.50%-87.32%)(表 5)。MDSs按照1-3-5-7-9-11-2-4-6-8-10-12-13-14的顺序排列, 标号为奇数的MDSs与标号为偶数的MDSs各自簇集在一起(图 2)。
表 5 Oxytricha sp. (ZJ) 小核α-TBP基因的MDSsIESs及指针序列Table 5. Characteristics of MDSs, IESs, and pointer sequences in Oxytricha sp. (ZJ)α-TBP geneMDSn MDSn的长度Length of MDSn MDSn和(n+1) 间的指针序列Pointer sequence between MDSn and (n+1) MDSi和MDSj间IES的长度Length of IESs between MDSi and MDSj IESs AT含量IESs AT content (%) 1 223 GGTGCAC 1-3 21 52.38 2 282 GTCTYGTA 3-5 33 72.73 3 35 ATTCACA 5-7 29 68.97 4 22 ATTAGAGT 7-9 42 78.57 5 158 TCGCWCATTC 9-11 16 62.50 6 27 AGAAGTC 11-2 50 86.00 7 13 TGCCACT 2-4 20 65.00 8 173 TCAAGCTTAA 4-6 45 73.33 9 21 AGTTTT 6-8 19 68.42 10 48 CTGGTGA 8-10 30 76.67 11 26 ATCAGCTACTT 10-12 25 76.00 12 28 AAAGT 12-13 34 85.29 13 731 TAAGTRAT 13-14 71 87.32 14 312 注: MDSs和IESs的长度都包括了一拷贝的指针序列; 不同克隆间指针序列的差异用兼并碱基表示
Note: Length of MDSs include pointers on both ends; IESs exclude pointer sequences; Use degenerate symbols reflects differences between clones2.3 Oxytricha sp. (ZJ) 的大、小核DNApol α基因结构特征
Oxytricha sp. (ZJ) 大核DNApol α基因全长为5118 bp。含有编码1534个氨基酸的ORF, 其AT含量约为62.15%。含两个长度分别为96 bp (intron-1) 和49 bp (intron-2) 的内含子, AT含量分别为71.88%和69.39%。5′和3′非编码区的序列长度分别为158 bp (含20 bp端粒序列) 和216 bp (含20 bp端粒序列), AT含量分别为71.85%和74.09%。5个阳性克隆测序累计检测到205个碱基变异位点, 其中ORF含有184个碱基变异位点, 涉及39个氨基酸的差异, 未引起移码突变或无义突变。剩余的碱基变异位点, 有5个位于5′非编码区, 5个位于3′非编码区, 8个位于intron-1, 3个位于intron-2。
已获得Oxytricha sp. (ZJ) 的小核DNApol α基因片段长度为2068 bp。包括15个富含AT序列(62.50%-90.91%) 的IESs, 14个包含两端指针序列的完整MDSs (5-7-9-11-13-15-17-19-21-23-25-27-29-30)(表 6), 以及2个MDSs (2和32) 的部分区域。并由此确定了另13个MDSs (1-6-8-10-12-14-16-18-20-22-24-26-28) 的位置(不包含指针序列), 指针序列只能在相连的3个MDSs (29-30-31) 间被确定出来(图 3)。
表 6 Oxytricha sp. (ZJ) 小核DNApol α基因片段MDSs, IESs及指针序列Table 6. Characteristics of MDSs, IESs, and pointer sequences in Oxytricha sp. (ZJ) DNApol α gene segmentMDSn MDSn的长度Length of MDSn MDSn和(n+1) 间的指针序列Pointer sequence between MDSn and (n+1) MDSi和MDSj间IES的长度Length of IESs between MDSi and MDSj IESs AT含量IESs AT content (%) 2 ? ? 2-5 24 62.50 5 39 ? 5-7 22 77.27 7 85 ? 7-9 11 90.91 9 70 ? 9-11 17 82.35 11 39 ? 11-13 20 80.00 13 55 ? 13-15 26 84.62 15 32 ? 15-17 5 80.00 17 101 ? 17-19 42 78.57 19 86 ? 19-21 17 88.24 21 100 ? 21-23 14 85.71 23 82 ? 23-25 31 83.87 25 40 ? 25-27 14 85.71 27 97 ? 27-29 50 74.00 29 418 AACAT 29-30 43 79.09 30 160 TGAAGAA 30-31 85 82.35 31 ? 注: MDSs的长度包括了两端的指针序列; IESs的长度不含指针序列; 问号代表数值无法确定
Note: Length of MDSs include pointers on both ends; IESs exclude pointer sequences; Question marks indicate undetermined information3. 讨论
3.1 尖毛虫属小核Actin Ⅰ基因的乱序模式
Actin Ⅰ基因是腹毛类纤毛虫中最早发现的乱序基因[4]。目前尖毛虫属已知小核Actin Ⅰ基因结构的物种共包括5个:Oxytricha trifallax、O. nova、O. fallax、O. sp. (Aspen)、O. sp. (Misty)。根据Hogan等[12]推测的腹毛类小核Actin Ⅰ基因结构进化路线, 尖毛虫属内小核Actin Ⅰ基因乱序模式分为两种情况:O. sp. (Misty) 与Stylonychia pustulata拥有较为相似的模式(模式1);O. trifallax、O. nova、O. sp. (Aspen) 模式(模式2) 相近。而之后报道的O. fallax与我们新测序的O. sp. (ZJ) 模式可归入模式2。
尖毛虫属内小核Actin Ⅰ基因乱序模式2所涉及的5个物种可以分为2种情况: (1)Oxytricha fal-lax与O. trifallax两个种群Actin Ⅰ的小核基因乱序模式均为3-4-6-5-7-9-10-2-1-8 (模式2.1)。DuBois和Prescott等[12]将三者视为O. trifallax-O. fallax群组, 代表了一个单一种。(2)O. sp. (ZJ) 小核Actin Ⅰ基因由8个IESs和9个MDSs组成, 同属内的O. nova在组成数目上完全相同, 且2个物种的MDSs排列顺序均为3-4-6-5-7-9-2-1-8 (模式2.2)[1], 比模式2.1少了长的IES6与非常短的MDS10。O. sp. (Aspen)[11]仅有部分小核Actin Ⅰ基因完成测序, Hogan等[12]推测其乱序模式与O. nova (模式2.2) 相似, 鉴于数据的不完整性, 在此我们不做详细讨论。
对比基于大核Actin Ⅰ基因序列构建的系统树(图 4a) 与小核Actin Ⅰ基因结构进化路线(图 4b) 可以看出, 二者存在一定的相似之处, 但却不完全相同。在图 4a中作为姐妹枝出现的Oxytricha fallax与O. trifallax的小核Actin Ⅰ乱序模式均为2.1, 拥有10个MDSs (图 4b)。但是, 在大核Actin Ⅰ基因系统树中(图 4a), 同属于小核Actin Ⅰ乱序模式2.2 (图 4b) 的3个物种O. nova、O. sp. (Aspen)、O. sp. (ZJ) 并未聚在一起。这表明, 大小核Actin Ⅰ基因进化路线并不完全一致。
模式2.1、2.2之间, 在IESs的长度方面均具有较大差异。在模式2.1中, Oxytricha trifallax的IESs长度为19-276 bp[10]; 而在模式2.2中, O. sp. (ZJ) 和O. nova的IESs长度分别是11-79 bp (本研究)、11-109 bp[4]。这符合前人的推测:由于IESs的快速进化与易变性, 在进化过程中其长度也发生了改变, 而小核Actin Ⅰ基因进化过程中, IESs的长度比其序列更为保守[15]。凭IESs的序列很难推断进化的图式。
在模式2.1中, Oxytricha trifallax的指针序列长度为2-14 bp[11]; 而在模式2.2中, O. nova和O. sp. (ZJ) 的指针序列长度范围较为接近, 分别为5-13 bp[4]和4-13 bp (表 4)。但相同乱序模式的物种指针序列并未发现具有更多序列相似性。
3.2 尖毛虫属小核α-TBP基因的乱序模式
目前尖毛虫属已知小核α-TBP基因结构的物种共包括2个:Oxytricha trifallax、O. nova。我们的研究表明, O. sp. (ZJ) 小核α-TBP基因由13个IESs和14个MDSs组成, 与同属的O. nova一致[8], 说明在这两个种分化过程中该小核基因没有发生新的MDS片段的重组或易位。与之相对比, O. trifallax多了位于非乱序区的3个MDSs和3个IESs[15](图 2)。猜测这两种不同的模式是由于IES的插入/丢失造成。O. sp. (ZJ) 与O. nova的MDSs排列模式为1-3-5-7-9-11-2-4-6-8-10-12-13-14[8]; 而O. trifallax排列模式为1-3-5-7-10-12-2-4-6-8-9-11-13-14-15-16-17[15]。O. sp. (ZJ)、O. nova中的MDS8、MDS13-14分别与O. trifallax中的MDS8-9、MDS14-17同源(图 6)。
![]()
图 5 Oxytricha sp. (ZJ)、O. nova和O. trifallax大、小核Actin Ⅰ基因结构模式图(参照DuBois等, 1995[10])梯形代表端粒, 黑色圆点代表起始密码子与终止密码子Figure 5. Schematic illustration of macro-and micronuclear Actin Ⅰ gene structures of Oxytricha sp. (ZJ), O. nova and O. trifallax (reference DuBois, et al., 1995[10])Ladder represent telomeres, black dots mark the pointer of start and stop codons![]()
Oxytricha sp. (ZJ)、O. nova的IESs数目均为13个, 二者的长度差别不大, 分别为16-71 bp、11-55 bp; 而O. trifallax的IESs数目为16个, 长度为5-108 bp[15]。O. sp. (ZJ) 和O. nova中的IESs10, 11, 12与O. trifallax中的IESs11, 12, 13相对应。与小核Actin Ⅰ的结果相似, O. sp. (ZJ) 与O. nova具有同样α-TBP基因乱序模式, 且IESs的长度较为一致, O. sp. (ZJ)、O. nova、O. trifallax的IES平均长度分别为33、26和35 bp。O. sp. (ZJ)(16 bp) 和O. nova (8 bp) 最短的IES是IES5, 而O. trifallax最短的为IES9 (5 bp)[8, 15]。由此可见, IESs的长度、数目具有物种特异性。
模式2.1的Oxytricha trifallax指针序列长度为7-17 bp[15], 模式2.2的O. nova和O. sp. (ZJ) 的指针序列长度分别为为3-19 bp[8]、5-11 bp (表 5), 长度范围差异大。相同乱序模式的物种指针序列并未发现具有更多序列相似性。
3.3 尖毛虫属小核DNApol α基因的乱序模式
与小核Actin Ⅰ和α-TBP基因相比, DNApol α基因乱序程度更为复杂, 基因序列更长, 乱序模式为非随机排列[9]。目前, 尖毛虫属已知小核DNApol α基因结构的物种共包括2个:Oxytricha trifallax、O. nova。我们的研究未能获取O. sp. (ZJ) 的小核DNApol α基因全长。根据扩增出的O. sp. (ZJ) 小核DNApol α基因部分序列, 我们确定出14个完整MDSs (排列模式为5-7-9-11-13-15-17-19-21-23-25-27-29-30), 并推测出未扩增成功的16个MDSs (MDSs1-2、偶数MDSs6-28、MDSs31-32) 的位置, 同时推测MDS3与MDS4可能丢失。O. trifallax包含至少51个MDSs和50个IESs[17], O. nova包含至少45个MDSs和44个IESs[9]。通过比对分析(图 7), 我们发现O. sp. (ZJ) 含有3个非乱序排列的MDSs, 分别是MDS29、MDS30、MDS31, 而在O. trifallax和O. nova中分别是9个和7个非乱序排列的MDSs[17]。
![]()
模式2.1的Oxytricha trifallax指针序列长度为2-15 bp[17], 模式2.2的O. nova指针序列长度为2-16 bp[9]。本研究中未获得O. sp. (ZJ) 的DNApol α基因序列全长, 很多指针序列无法推测(表 6), 故而本研究未对该基因不同物种间指针序列进行对比。
我们确定出的Oxytricha sp. (ZJ) 小核DNApol α基因的14个完整MDSs位于倒置区域。O. sp. (ZJ) 的MDS7对应O. nova的MDS6, 与O. trifallax的MDS7-8-9同源, 并且MDS8在O. trifallax的小核基因主体部分中是缺失的[17]。据此合理推断是: (1) 此区域由3个MDSs组成, 且中间MDS的缺失代表了祖先种的基因结构特征。O. trifallax从祖先种分离时保留了该结构特征。而在O. sp. (ZJ) 和O. nova分化前的进化过程中将3个MDSs融合成了一个MDS。(2) 为第一种可能的反向进程, 即祖先种在此区域为一个MDS结构, 而在O. trifallax从祖先种分化时, 独立的发生了一次重组将中间的部分转离了小核基因的主体部分, 类似移位现象, 使得一个MDS变成了3个MDSs。之前Wong等[16]的研究表明, 在种系分离自共同祖先后, 独立发生编码区域的重新划分和融合。
因Oxytricha sp. (ZJ) 的MDS7与O. trifallax的MDS7-8-9同源, 且O. trifallax在此处的IES与O. sp. (ZJ) 的MDS7对应区域的序列完全一致(图 8), 表明在O. sp. (ZJ) 的进化过程中可能将该IES同化成了MDS的一部分, 即发生IES向MDS的转换, 伴随此种转换, 使得原先3个MDSs融合成了一个MDS, 而原先位于中间的MDS被丢弃。这个转换过程呈现出腹毛类纤毛虫小核基因在进化过程中方式的多样性。这支持了上述第一种解释, 即在O. sp. (ZJ) 和O. nova分化前的进化过程中将3个MDSs融合成了一个MDS。Duharcourt等[30]在Paramecium中与Chalker与Yao[31]在Tetrahymena中都发现了旧大核的IESs在新大核中的保留现象。而Mollenbeck等[19]通过比较Stylonychia lemnae两种群的小核Actin Ⅰ基因发现, 发现IES向MDS转变的迹象, 但这种转变发生在非编码区, 未对基因的开放阅读框造成影响。本实验则发现编码区内IES向MDS的转变, 使得本应删除的序列成为基因组永久保留的一部分。
![]()
图 8 Oxytricha sp. (ZJ) 进化过程中IES向MDS的转变Figure 8. IES to MDS transformation during Oxytricha sp. (ZJ) evolution -
[1] Rnnestad I, Yfera M, Ueberschr B, et al. Feeding behaviour and digestive physiology in larval fish: current knowledge, and gaps and bottlenecks in research [J]. Reviews in Aquaculture, 2013, 5(Suppl. 1): s59-s98
[2] Yu H R, Mai K S, Ma H M, et al. Evaluation of microdiets versus frozen copepods on growth, survival and selected chemical composition of large yellow croaker Pseudosciaena crocea larvae [J]. Acta Hydrobiologica Sinica, 2012, 36(1): 49-56 [于海瑞, 麦康森, 马洪明, 等. 微颗粒饲料与冰冻桡足类对大黄鱼稚鱼生长、存活和体成分的影响. 水生生物学报, 2012, 36(1): 49-56]
[3] Gouin S. Microencapsulation: industrial appraisal of existing technologies and trends [J]. Trends in Food Science Technology, 2004, 15: 330-347
[4] Xie Z, Wang F, Lou B, et al. Retention efficiency and release of nutrients in the digestive tract of larval shrimp (Penaeus japonicus Bate) using different microencapsulated diets [J]. Aquaculture International, 2013, 21: 667-678
[5] Lou B, Shi H L, Mao G M, et al. Studies on techniques of the artificial breeding and large size gingerling cultivating of Spotted Maigre (Nibea albiflora) [J]. Modern Fisheries Information, 2011, 26(3): 20-23 [楼宝, 史会来, 毛国民, 等. 黄姑鱼全人工繁育及大规格苗种培育技术研究. 现代渔业信息, 2011, 26(3): 20-23]
[6] Xie Z, Wang F, Liu H, et al. Gelatin-walled microencapsulated diet for larval shrimp (Penaeus japonicus Bate) manufactured using fluidized bed coating process [J]. Aquaculture Research, 2010, 42: 65-73
[7] Chong A, Hashim R, Chow-Yang, et al. Partial characterization and activities of proteases from the digestive tract of discus fish (Symphysodon aequifasciata) [J]. Aquaculture, 2002, 203: 321-333
[8] Natalia Y, Hashim R, Ali A, et al. Characterization of digestive enzymes in a carnivorous ornamental fish, the Asian bony tongue Scleropages formosus (Osteoglossidae) [J]. Aquaculture, 2004, 233: 305-320
[9] Bessey O A, Lowry O H, Brock M J. A method for the rapid determination of alkaline phosphatase with five cubic millimeters of serum [J]. The Journal of Biological Chemistry, 1946, 164: 321-329
[10] Bradford M M. A rapid and sensitive method for the quantization of microgram quantities of protein utilizing the principle of protein dye binding [J]. Analytical Biochemistry, 1976, 72: 248-254
[11] Agrawal R, Naveen Y. Pharmaceutical processing-a review on wet granulation technology [J]. International Journal of Pharmaceutical Frontier Research, 2011, 1(1): 65-83
[12] Gibbs B F, Kermasha S, Alli I, et al. Encapsulation in the food industry: a review [J]. International Journal of Food Sciences and Nutrition, 1999, 50: 213-224
[13] Rosas C, Tut J, Baeza J, et al. Effect of type of binder on growth, digestibility, and energetic balance of Octopus maya [J]. Aquaculture, 2008, 275: 291-297
[14] Chen X B. The related research on microdiets processing technology of flouder fish [D]. Thesis for Master of Science. Ocean University of China, Qingdao, 2011 [陈笑冰. 鲆鲽鱼类微颗粒饲料加工工艺的相关研究. 硕士学位论文, 中国海洋大学, 青岛, 2011]
[15] Niu H X, Zhu A X, Chang J, et al. Study on preparation process and release mechanism of microencapsulation lysine [J]. Science and Technology of Food Industry, 2012, 125(8): 16-21 [牛化欣, 祝爱侠, 常杰, 等. 微胶囊赖氨酸制备工艺及其缓释效果研究. 食品工业科技, 2012, 125(8): 16-21]
[16] Anas A, Philip R, Singh I S B. Chitosan as a wall material for a microencapsulated delivery system for Macrobrachium rosenbergii (de Man) larvae [J]. Aquaculture Research, 2008, 39: 885-890
[17] Partridge G J, Southgate P C. The effect of binder composition on ingestion and assimilation of microbound diets (MBD) by barramundi Lates calcarifer Bloch larvae [J]. Aquacultutre Research, 1999, 30: 879-886
[18] Fernndez-Daz C, Yfera M. Capacity of gilthead seabream, Sparus aurata L., larvae to break down dietary microcapsules [J]. Aquaculture, 1995, 134: 269-278
[19] Xie Z, Wang F, Zhu A, et al. Effects of diets microencapsulated with different wall materials on growth and digestive enzymes of the larvae of Penaeus japonicus Bate [J]. Journal of Shellfish Research, 2011, 30: 133-138
[20] Tesser M B, Portella M C. Degradation analysis of microencapsulated diet in pacu (Piaractus mesopotamicus Holmberg, 1887) larvae intestine through scanning electron microscopy (SEM) [J]. Acta Scientiarum Animal Sciences, 2003, 1: 49-52
-
期刊类型引用(1)
1. 何燕，尹家奇，生欣. 艾美游仆虫大核基因组与转录组测序及结构特征. 水生生物学报. 2021(05): 1014-1023 . 
本站查看
其他类型引用(0)
下载:
计量
- 文章访问数: 1301
- HTML全文浏览量: 0
- PDF下载量: 718
- 被引次数: 1