太湖新银鱼卵巢发育不同阶段个体生化组成和能量密度研究
CHEMICAL COMPOSITIONS AND ENERGY DENSITIES OF FEMALE NEOSALANX TAIHUENSIS WITH OVARIES AT DIFFERENT DEVELOPMENTAL STAGES
-
摘要: 对丹江口水库太湖新银鱼雌鱼性腺发育不同阶段个体的生化组成和能量密度变化进行研究, 以分析性腺发育的物质和能量来源。卵巢8-11月间基本处于Ⅱ期, 之后迅速发育; 1-3月主要为Ⅳ期或Ⅴ期, 4月开始以Ⅵ期个体为主。卵巢发育的Ⅱ期和Ⅲ期, 太湖新银鱼全鱼中的生化组成和能量密度随卵巢发育而显著增高。当发育至Ⅳ期时, 全鱼和去性腺鱼体中的干物质、粗脂肪、粗灰分相对含量和能量密度仍明显增加, 而粗蛋白相对含量在全鱼当中有所减少。由Ⅳ期发育到Ⅴ期时, 全鱼中的各生化组成无明显变化, 能量密度显著增加; 去性腺鱼体中的粗蛋白相对含量显著增加, 粗脂肪显著减少, 其他生化组成和能量密度无明显变化。产后Ⅵ期, 全鱼中的粗脂肪、粗灰分相对含量和能量密度显著下降, 干物质和粗蛋白相对含量则显著上升。性腺发育成熟前(Ⅳ期前), 去性腺鱼体中的总能量和粗蛋白绝对含量明显增加, 而性腺中能量和粗蛋白绝对含量增加不明显, 粗脂肪含量的增加明显。由Ⅳ期发育到Ⅴ期时, 全鱼中的总能量、粗蛋白和粗脂肪绝对含量均略有增加, 主要是由于性腺中能量、粗蛋白和粗脂肪含量的增加, 而去性腺鱼体中能量、粗蛋白和粗脂肪绝对含量均明显减少。产后Ⅵ期全鱼中的总能量、粗蛋白和粗脂肪绝对含量均急剧降低。上述结果初步表明太湖新银鱼卵巢发育成熟前在鱼体中储存物质和能量, 在卵巢发育成熟时这些储存的物质和能量可能向卵巢转移, 成为卵巢发育成熟的主要物质和能量来源。Abstract: Chemical composition and energy density were analyzed for female Neosalanx taihuensis with ovaries at different developmental stages to determine the energy source of gonadal development. Ovaries of females were at stage Ⅱor Ⅲ from August to December, and developed quickly to stage Ⅳ or Ⅴ from January to March of the next year, then reached stage Ⅵ in April. Chemical composition and energy density in the whole fish increased significantly from stage Ⅱ to stage Ⅲ. Contents of dry matter, lipids and ash, and energy density increased significantly in the whole fish and fish both with and without ovaries from stage Ⅲ to Ⅳ, whereas relative protein content decreased in the whole fish. At mature stage (Ⅴ), energy density in the whole fish increased, although chemical compositions did not have significant change. In fish body without ovaries, relative contents of protein increased, while lipid decreased, and chemical compositions and energy density had little changes. At the spent stage (Ⅵ), relative contents of lipid and ash decreased in the whole fish, conversely, relative contents of dry matter and protein increased. Total energy and absolute protein content in fish body without ovaries increased from stage Ⅲ to Ⅳ, while there was little increase if ovaries were included in the analysis. Absolute lipid content in ovaries increased from stage Ⅲ to Ⅳ. At mature stage (Ⅴ), total energy, absolute protein content, and lipid content increased slightly in the whole fish. This change primarily reflected their increases in ovaries. In the fish without ovaries, they decreased apparently. At the spent stage (Ⅵ), total energy, absolute protein content and lipid content decreased sharply in the whole fish. Our preliminary results indicated that materials and energy were stored in fish body during the premature stages, and were then transferred to ovaries at the mature stage.
-
纤毛虫原生动物具有二型核, 小核基因通常被内部删除序列(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] Ni Y, Zhu C D. Fishes of the Taihu Lake[M]. Shanghai: Shanghai Scientific Technical Publishers. 2005, 217-221[倪勇, 朱成德. 太湖鱼类志. 上海:上海科学技术出版社. 2005, 217-221]
[2] Wang X L, Xiang J G, Liu J S, et al. The influences to transplantation of icefish (Neosalanx taihuensis chen) on early-life feeding and growth of Hemiculter leucisculus basilewsky. Acta Hydrobiologica Sinica, 2012, 36(2): 276-283[王星璐, 向建国, 刘家寿, 等. 太湖新银鱼移植对 早期摄食和生长的影响. 水生生物学报, 2012, 36(2): 276-283]
[3] Hu C L, Chen W X, Liu J S. Situation and strategy research of icefish transplant in China[J]. Reservoir Fisheries, 1998, (2): 3-7[胡传林, 陈文祥, 刘家寿. 我国银鱼移植增殖现状及对策分析. 水利渔业, 1998, (2): 3-7]
[4] Wu L, Liu J S, Wang X L, et al. Identification of indi born in different spawning seasons using otolith microstructure to reveal life history of Neosalanx taihuensis[J]. Fisheries Sciences, 2011, 77(3): 321-327
[5] Gong W B, Li H T, Xie S Q, et al. Two spawning stocks of icefish Neosalanx taihuensis revealed from annual reproductive cycle analysis[J]. Fisheries Sciences, 2009, 75(5): 1157-1165
[6] Gong W B, Wu L, Cheng F, et al. A comparative study on reproductive characteristics of the spring and autumn spawning stocks of Neosalanx taihuensis Chen in the Three Gorges Reservoir[J]. Acta Hydrobiology Sinica, 2010, 34(5): 1065-1068[龚望宝, 吴朗, 程飞, 等. 三峡水库太湖新银鱼春季和秋季繁殖群体的繁殖生物学特征比较. 水生生物学报, 2010, 34(5): 1065-1068]
[7] Otero J, Gonzlez A F, Sieiro M P, et al. Reproductive cycle and energy allocation of Octopus vulgaris in Galician waters, NE Atlantic[J]. Fisheries Research, 2007, 85(1-2): 122-129
[8] Raymond J, Himmelman J H, Guderley H E. Biochemical content, energy composition and reproductive effort in the broadcasting sea star Asterias vulgaris over the spawning period[J]. Journal of Experimental Marine Biology and Ecology, 2007, 341(1): 32-44
[9] Joaquim S, Matias D, Matias M, et al. Reproductive activity and biochemical composition of the pullet carpet shell Venerupis senegalensis (Gmelin, 1791) (Mollusca: Bivalvia) from Ria de Aveiro (northwestern coast of Portugal)[J]. Scientia Marina, 2011, 75(2): 217-226
[10] Rodrguez-Gonzlez H, Hernndez-Llamas A, Villarreal H, et al. Gonadal development and biochemical composition of female crayfish Cherax quadricarinatus (Decapoda: Parastacidae) in relation to the Gonadosomatic Index at first maturation[J]. Aquaculture, 2006, 254(1-4): 637-645
[11] Ojea J, Pazos A J, Martnez D, et al. Seasonal variation in weight and biochemical composition of the tissues of Ruditapes decussatus in relation to the gametogenic cycle[J]. Aquaculture, 2004, 238(1-4): 451-468
[12] Yin M C. Fish Ecology[M]. Beijing: Agricultural Publishers of China. 1995, 111-112[殷名称. 鱼类生态学. 北京: 中国农业出版社. 1995, 111-112]
[13] Ren J S, Marsden I D, Ross A H, et al. Seasonal variation in the reproductive activity and biochemical composition of the Pacific oyster (Crassostrea gigas) from the Marlborough Sounds, New Zealand[J]. New Zealand Journal of Marine and Freshwater Research, 2003, 37(1): 171-182
[14] Dridi S, Romdhane M S, Elcafsi M. Seasonal variation in weight and biochemical composition of the Pacific oyster, Crassostrea gigas in relation to the gametogenic cycle and environmental conditions of the Bizert lagoon[J]. Aquaculture, 2007, 263(1-4): 238-248
[15] Dima J B, Vido N A, Leal G A, et al. Fluctuations in the biochemical composition of the Patagonian stone crab Platyxanthus patagonicus A. Milne Edwards, 1879 (Platyxanthidae: Brachyura) throughout its reproductive cycle[J]. Scientia Marina, 2009, 73(3): 423-430
[16] Bayne B L. Aspects of reproduction in bivalve molluscs[A]. In: Wiley M L (Eds.), Estuarine Processes[C]. New York: Academic Press. 1976, 432-448
[17] Liu W G, Li Q, Shirasu K, et al. Reproductive cycle and biochemical composition of the Zhe oyster Crassostrea plicatula Gmelin in an eastern coastal bay of China[J]. Aquaculture, 2006, 261(2): 752-759
[18] Najmudeen T M. Variation in biochemical composition during gonad maturation of the tropical abalone Haliotis varia Linnaeus 1758 (Vetigastropoda: Haliotidae)[J]. Marine Biology Research, 2007, 3(6): 454-461
[19] Castille F L, Lawrence A L. Relationship between maturation and biochemical composition of the gonads and digestive glands of the shrimps Penaeus aztecus Ives and Penaeus setiferus (L.)[J]. Journal of Crustacean Biology, 1989, 9(2): 202-211
[20] Li Q, Yang L, Ke Q Z, et al. Gametogenic cycle and biochemical composition of the clam Mactra chinensis (Mollusca: Bivalvia): Implications for aquaculture and wild stock management.[J]. Marine Biology Research, 2011, 7(4): 407-415
[21] Han D J, Peng J H, Jian D, et al. The assessment of food organism resources and trophic states of Danjiangkou Reservoir[J]. Journal of Lake Sciences, 1997, 9(1): 57-62[韩德举, 彭建华, 简东, 等. 丹江口水库的饵料生物资源及水体营养状态评价. 湖泊科学, 1997, 9(1): 57-62]
[22] Gong S Y, Wu M S, Jin D X, et al. Primary study on feeding habits of Neosalanx taihuensis in Fushui Reservoir[J]. Reservoir Fisheries, 1999, 19(3): 7-9[龚世园, 吴敏生, 金顶玄, 等. 富水水库太湖新银鱼食性的初步研究. 水利渔业. 1999, 19(3): 7-9]
[23] Li Y X. Reproductive Ecology of Anadromous Coilia ectenes in the Yangtze River[D]. Thesis for Doctor of Science. Institute of Hydrobiology, Chinese Academy of Sciences, Wuhan. 2009[黎雨轩. 长江洄游性刀鲚的繁殖生态学研究. 博士学位论文, 中国科学院水生生物研究所, 武汉. 2009]
-
期刊类型引用(1)
1. 何燕，尹家奇，生欣. 艾美游仆虫大核基因组与转录组测序及结构特征. 水生生物学报. 2021(05): 1014-1023 . 
本站查看
其他类型引用(0)
下载:
计量
- 文章访问数: 1180
- HTML全文浏览量: 4
- PDF下载量: 622
- 被引次数: 1