尖毛虫属Actin Ⅰ、α-TBP和DNA pol α乱序基因的模式研究

杨 然 陈天兵 黄 俊 石文俊 区淑青 伊珍珍

(华南师范大学生命科学学院, 广州市亚热带生物多样性与环境生物监测重点实验室, 广州 510631)

尖毛虫属Actin Ⅰ、α-TBP和DNA pol α乱序基因的模式研究

杨 然 陈天兵 黄 俊 石文俊 区淑青 伊珍珍

(华南师范大学生命科学学院, 广州市亚热带生物多样性与环境生物监测重点实验室, 广州 510631)

研究旨在对尖毛虫属内现有物种的3种乱序小核基因结构进行比较, 探讨其乱序模式。于湛江湖光红树林水域中采集到一个尖毛虫属物种Oxytricha sp. (ZJ), 成功扩增了其肌动蛋白Ⅰ(ActinⅠ)、端粒结合蛋白(α-TBP)、DNA聚合酶α (DNA pol α)3个乱序基因的完整大核基因序列和完整/部分小核基因序列, 并结合已有资料对比研究了尖毛虫属这3个乱序基因的进化。结果表明: (1) Oxytricha sp. (ZJ)与O. nova的小核Actin Ⅰ基因具有相同的乱序模式, 区别于其余的尖毛虫属物种; 在增加尖毛虫属物种的基础上, 对前人推测提出了质疑,我们认为MDS-IES接合处移动现象在乱序MDSs之间并非比非乱序MDSs之间更保守。(2) Oxytricha sp. (ZJ)与O. nova的小核α-TBP基因具有相同乱序模式和相似长度的IESs。(3) Oxytricha sp. (ZJ)的小核DNA pol α基因乱序模式, 区别于任一已报道物种, 与属内O. trifallax最为相近。基于序列分析, 在DNA pol α基因中发现了一例IES转换为MDS的痕迹, 以及由此导致原先MDS的丢失。研究发现在编码区内IES向MDS的转变, 使得本应删除的序列成为基因组永久保留的一部分。

尖毛虫; 乱序基因; 大核命运序列; 内部删除序列; 肌动蛋白Ⅰ型; 端粒结合蛋白; DNA聚合酶α

纤毛虫原生动物具有二型核, 小核基因通常被内部删除序列(Internal eliminated sequences, IESs)分隔成多个大核命运序列(Macronucleus destined sequences, MDSs)片段。在接合生殖过程中, 通过小核IESs的删除和MDSs的重组连接、端粒序列的生成和基因的高度扩增等过程形成具功能的大核基因[1—5]。染色体重排广泛存在于3种纤毛虫类群中: Armophorea、旋唇纲Spirotrichea、层咽纲Phyllopharyngea[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个物种的小核DNA pol α基因[9,17—20]。Hogan等[12]推测了腹毛类纤毛虫小核Actin Ⅰ基因的进化路线,从共同的原始祖先到3个MDSs的Urostyla grandis,进化为拥有4个MDSs的Engelmanniella mobilis, MDS2发生倒置, 再进化为拥有8个MDSs的Stylonychia lemnae, 然后分化为两个分支, 一支为3个尖毛虫类群, 一支为S. pustulata和Oxytricha sp. (Misty)。之后, Chen等[21]在成功扩增2种尾柱目纤毛虫大小核Actin Ⅰ基因的基础上, 将腹毛类小核Actin Ⅰ进化路线进一步修正, 提出腹毛类不同目之间可能存在不同的进化路线。Wong等[16]推测了α-TBP基因可能的进化路线, 认为祖先类群为非乱序物种(Holosticha sp.), 历经乱序排列的Uroleptus sp.、Paraurostyla weissei, 进化为中间MDS发生易位的O. trifallax, 再经历两个MDSs发生融合(如O. nova/S. mytilus)。DNA pol α的进化路线尚未有人给出。综上,已报道的腹毛类3种乱序基因进化路线仅为较为粗略的大致框架, 细节部分依然亟待修订和补充。

本研究在湛江湖光红树林水域中采集得到一个腹毛类尖毛虫物种Oxytricha sp. (ZJ)。通过测定O. sp. (ZJ)大、小核的Actin Ⅰ、α-TBP和DNA pol α基因序列, 并对尖毛虫属内现有物种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 Ⅰ 基因所用引物名称及序列Tab. 1 The name and sequence of primers for macro- and micronuclear Actin Ⅰ gene

基于上述获得的序列设计引物J-F (1203)、JM3F1、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为上下游引物做第二轮巢式P C R, A P 2 (表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。扩增产物经切胶回收、转化、克隆测序, 拼接得到大核基因序列全长, 并提交至G e n B a n k, 登录号为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)。得到目的条带。切胶回收、转化、克隆测序, 获得小核部分基因片段。

在已获得的小核基因序列的基础上设计多条引物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/JM8R进行巢式PCR, 扩增得到片段J-F (42)/J-M8R、J-M3F2/J-R (1550)、J-I2F/J-M8R。覆盖小核基因全长。拼接获得小核基因序列全长(图 1), 并提交至GenBank, 登录号为KX602205。

表 2 大、小核α-TBP基因所用引物名称及序列Tab. 2 The name and sequence of primers for macro- and micronuclear α-TBP gene

大、小核α-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、JM9F、J-M8F1、J-M8F2、J-M11F、J-M10R、J-I9F、J-R (2133)、J-I13R2 (表 2、图 2), 并提交至Gen-Bank, 大核基因登录号为KX602208, 小核基因为KX602209。

大、小核DNA pol α基因序列的获得基于上文描述的方式获得Oxytricha sp. (ZJ)的DNA pol α基因完整的大核序列。引物为Pol α-F、Pol α-R、 J-R2 (DPα)、J-R1 (DPα)、J-F3 (DPα)、J-F4 (DPα) (表 3), 得到DNA pol α基因完整的大核序列。

因在DNA pol α基因小核序列扩增时, 扩增结果均为大、小核基因序列的嵌合体, 故只获得部分序列。所用引物为J-M32F、J-M1R1、J-M1R2。JL1F、J-R5R、J-R6R、J-R9R、J-5000R、J-5100R (表 3、图 3), 并提交至GenBank, 大核基因登录号为KX602206, 小核基因为KX602207。

大核Actin Ⅰ基因系统树的构建因为已完成大核α-TBP基因和DNA pol α基因测序的尖毛虫物种数目太少(仅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 inference, 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]。

图 1 Oxytricha sp. (ZJ)大、小核Actin Ⅰ基因结构、引物位置及获得的PCR产物示意图Fig. 1 Schematic illustration of macro- and micronuclear structures, PCR primer locations and products obtained in Oxytricha sp. (ZJ) Actin Ⅰ gene灰色方块代表MDSs, 黑色方块代表IESs; 其中MDS2倒置; 白色方块代表未得到的基因区域, 下同MDSs are represented in grey, IESs are represented in black; MDS2 is inverted; white box represents unavailable gene regim; the same applies below

表 3 大、小核DNA pol α基因所用引物名称及序列Tab. 3 The name and sequence of primers for macro- and micronuclear DNA pol α gene

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)。

基于尖毛虫属内8个物种的大核Actin Ⅰ基因序列建的贝叶斯树(图 4a)中, Oxytricha trifallax、O. fallax与O. granulifera以最高置信值聚为一支,其余5个物种O. longa、O. sp. (Misty)、O. nova、O. sp. (Aspen)、O. sp. (ZJ)聚为一支。

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个位于内含子。

图 2 Oxytricha sp. (ZJ)大、小核α-TBP基因结构、引物位置及获得的PCR产物示意图Fig. 2 Schematic illustration of macro- and micronuclear structures, PCR primer locations and products obtained in Oxytricha sp. (ZJ) α-TBP gene

图 3 Oxytricha sp. (ZJ)大、小核DNA pol α基因结构、引物位置及获得的PCR产物示意图Fig. 3 Schematic illustration of macro- and micronuclear structures, PCR primer locations and products obtained in Oxytricha sp. (ZJ) DNA pol α gene

表 4 Oxytricha sp. (ZJ)小核Actin Ⅰ基因的MDSs、IESs及指针序列Tab. 4 Characteristics of MDSs, IESs, and pointer sequences in Oxytricha sp. (ZJ) Actin Ⅰ gene

图 4 基于Actin Ⅰ基因序列构建的贝叶斯树(a)与Actin Ⅰ基因结构进化路线图(b)Fig. 4 Baysian trees based on the sequences of Actin Ⅰ gene (a), and evolutionary route map of Actin Ⅰ gene structure (b)

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)。

2.3 Oxytricha sp. (ZJ)的大、小核DNA pol α基因结构特征

Oxytricha sp. (ZJ)大核DNA pol α基因全长为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)的小核DNA pol α基因片段长度为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)。

3 讨论

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 pustu-lata拥有较为相似的模式(模式1); O. trifallax、O. nova、O. sp. (Aspen)模式(模式2)相近。而之后报道的O. fallax与我们新测序的O. sp. (ZJ)模式可归入模式2。

表 5 Oxytricha sp. (ZJ)小核α-TBP基因的MDSsIESs及指针序列Tab. 5 Characteristics of MDSs, IESs, and pointer sequences in Oxytricha sp. (ZJ) α-TBP gene

表 6 Oxytricha sp. (ZJ)小核DNA pol α 基因片段MDSs, IESs及指针序列Tab. 6 Characteristics of MDSs, IESs, and pointer sequences in Oxytricha sp. (ZJ) DNA pol α gene segment

尖毛虫属内小核Actin Ⅰ基因乱序模式2所涉及的5个物种可以分为2种情况: (1) Oxytricha fallax与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)。但相同乱序模式的物种指针序列并未发现具有更多序列相似性。

MDSs在进化的过程中存在MIJ (MDS-IES junctions)移动现象[10,15,16,20], 这种移动通常在有限范围之内。DuBois等[10]提出, 乱序MDSs之间的MIJ移动比非乱序MDSs之间的更保守, 移动距离更短。在本研究中(图 5), 我们得出不同结论。对比模式2.2 O. sp. (ZJ)与O. nova发生的MIJ移动, 乱序MDSs之间移动最大达79 bp (MDS7-MDS9), 唯一的非乱序MDSs (MDS3-MDS4)之间为161 bp, 符合前人推测, 乱序MDSs之间的MIJ移动更为保守。而对比模式2.1的O. trifallax与模式2.2的O. sp. (ZJ), 乱序MDSs之间MIJ的移动范围为12—66 bp, 非乱序MDS3/4的对应边缘却相对保守, 仅相差23 bp, 由此可见MIJ移动距离与MDSs是否乱序并无直接关系,预示着DuBois等[10]的观点并不具有通用性。有意思的是, 非乱序MDS3/MDS4之间的MIJ移动距离,同为模式2.2的O. sp. (ZJ)与O. nova之间移动距离为161 bp, 不同模式的O. trifallax与O. nova之间移动距离为138 bp[15], 而不同模式的O. trifallax与O. sp. (ZJ)之间移动距离为23 bp。我们推测祖先种的小核Actin Ⅰ基因在MDS3-4对应区域具有两个IESs和3个MDSs, 在O. nova (模式2.2)进化过程中丢失了前一个IES, 在O. trifallax (模式2.1)和O. sp. (ZJ)(模式2.2)分别的进化过程中丢失了后一个IES, 最终3个物种均在该区域保留了2个MDSs和1个IES, 造成了上述MIJ移动距离的差距。因此, 我们赞成DuBois等[10]的观点, 认为MDS3/4之间MIJ移动可能是由于IES的插入/缺失造成的。由于MIJ移动距离较大, 我们不赞成是指针序列的突变[15,16]造成此MIJ移动。同时我们注意到属于同一乱序模式2.2的两个物种[O. nova和O. sp. (ZJ)]并未按照同样模式保留或缺失IESs, 该现象有待我们进一步研究。

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)。

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), 长度范围差异大。相同乱序模式的物种指针序列并未发现具有更多序列相似性。

图 5 Oxytricha sp. (ZJ)、O. nova和O. trifallax大、小核Actin Ⅰ基因结构模式图(参照DuBois 等, 1995[10])Fig. 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

图 6 Oxytricha sp. (ZJ)、O. nova和O. trifallax大、小核α-TBP基因结构模式图(参照Prescott 等, 1998[15])Fig. 6 Schematic illustration of macro- and micronuclear α-TBP gene structures of Oxytricha sp. (ZJ), O. nova and O. trifallax (reference Prescott, et al., 1998[15])

3.3 尖毛虫属小核DNA pol α基因的乱序模式

与小核Actin Ⅰ和α-TBP基因相比, DNA pol α基因乱序程度更为复杂, 基因序列更长, 乱序模式为非随机排列[9]。目前, 尖毛虫属已知小核DNA pol α基因结构的物种共包括2个: Oxytricha trifallax、O. nova。我们的研究未能获取O. sp. (ZJ)的小核DNA pol α基因全长。根据扩增出的O. sp. (ZJ)小核DNA pol α基因部分序列, 我们确定出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)的DNA pol α基因序列全长, 很多指针序列无法推测(表 6), 故而本研究未对该基因不同物种间指针序列进行对比。

我们确定出的Oxytricha sp. (ZJ)小核DNA pol α基因的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的转变, 使得本应删除的序列成为基因组永久保留的一部分。

图 7 Oxytricha sp. (ZJ)、O. nova和O. trifallax大、小核DNA pol α基因片段结构示意图(参照Hoffman 等, 1997[17])Fig. 7 Schematic illustration of macro- and micronuclear DNA pol α gene structures of Oxytricha sp. (ZJ), O. nova and O. trifallax (reference Hoffman, et al., 1997[17])

图 8 Oxytricha sp. (ZJ)进化过程中IES向MDS的转变Fig. 8 IES to MDS transformation during Oxytricha sp. (ZJ) evolution

[1]Prescott D M, Greslin A F. Scrambled Actin Ⅰ gene in the micronucleus of Oxytricha nova [J]. Developmental Genetics, 1992, 13(1): 66—74

[2]Curtis E A, Landweber L F. Evolution of gene scrambling in ciliate micronuclear genes [J]. Annals of the New York Academy of Sciences, 1999, 870: 349—350

[3]Nowacki M, Landweber L F. Epigenetic inheritance in ciliates [J]. Current Opinion in Microbiology, 2009, 12(6): 638—643

[4]Greslin A F, Prescott D M, Oka Y, et al. Reordering of nine exons is necessary to form a functional actin gene in Oxytricha nova [J]. Genetics, 1989, 86(16): 6264—6268

[5]Prescott D M. Genome gymnastics: unique modes of DNA evolution and processing in ciliates [J]. Nature Reviews Genetics, 2000, 1(3): 191—198

[6]Katz L A, KovnerA M. Alternative processing of scrambled genes generates protein diversity in the ciliate Chilodonella uncinata [J]. Journal of Experimental Zoology Part B Molecular and Developmental Evolution, 2010, 314(6): 480—488

[7]Lauth M R, Spear B B, Heumann J, et al. DNA of ciliated protozoa: DNA sequence diminution during macronuclear development of Oxytricha [J]. Cell, 1976, 7(1): 67—74

[8]Mitcham J L, Lynn A J, Prescott D M. Analysis of a scrambled gene: the gene encoding α-telomere-binding protein in Oxytricha nova [J]. Genes & Development, 1992, 6(5): 788—800

[9]Hoffman D C, Prescott D M. The germline gene encoding DNA polymerase α in the hypotrichous ciliate Oxytricha nova is extremely scrambled [J]. Nucleic Acids Research, 1996, 24(17): 3337—3340

[10]DuBois M L, Prescott D M. Scrambling of the Actin Ⅰgene in two Oxytricha species [J]. Proceedings of the National Academy of Sciences of the United States of America, 1995, 92(9): 3888—3892

[11]DuBois M L, Prescott D M. Volatility of internal eliminated segments in germ line genes of hypotrichous ciliates [J]. Molecular and Cellular Biology, 1997, 17(1): 326—337

[12]Hogan D J, Hewitt E A, Orr K E, et al. Evolution of IESs and scrambling in the Actin Ⅰ gene in hypotrichous ciliates [J]. Proceedings of the National Academy of Sciences of the United States of America, 2001, 98(26): 15101—15106

[13]Dalby A B, Prescott D M. The scrambled Actin Ⅰ gene in Uroleptus pisces [J]. Chromosoma, 2004, 112(5): 247—254

[14]Möllenbeck M, Cavalcanti A R, Jönsson F, et al. Interconversion of germline-limited and somatic DNA in a scrambled gene [J]. Journal of Molecular Evolution, 2006, 63(1): 69—73

[15]Prescott J D, DuBois M L, Prescott D M. Evolution of the scrambled germline gene encoding a-telomere binding protein in three hypotrichous ciliates [J]. Chromosoma, 1998, 107(5): 293—303

[16]Wong L C, Landweber L F. Evolution of programmed DNA rearrangements in a scrambled gene [J]. Molecular Biology and Evolution, 2006, 23(4):756—763

[17]Hoffman D C, Prescott D M. Evolution of internal eliminated segments and scrambling in the micronuclear gene encoding DNA polymerase alpha in two Oxytricha species [J]. Nucleic Acids Research, 1997, 25(10): 1883—1889

[18]Landweber L F, Kuo T C, Curtis E A. Evolution and assembly of an extremely scrambled gene [J]. Proceedings of the National Academy of Sciences of the United States of America, 2000, 97(7): 3298—3303

[19]Chang W J, Stover N A, Addis V M, et al. A micronuclear locus containing three protein-coding genes remains linked during macronuclear development in the spirotrichous ciliate Holosticha [J]. Protist, 2004, 155(2): 245—255

[20]Chang W J, Bryson P D, Liang H, et al. The evolutionary origin of a complex scrambled gene [J]. Proceedings of the National Academy of Sciences of the United States of America, 2005, 102(42): 15149—15154

[21]Chen T B, Yi Z Z, Huang J, et al. Evolution of the germline Actin gene in hypotrichous ciliates: multiple nonscrambled IESs at extremely conserved locations in two Urostylids [J]. Journal of Eukaryotic Microbiology, 2014, 62(2): 188—195

[22]Swart E C, Bracht J R, Magrini V, et al. The Oxytricha trifallax macronuclear genome: a complex eukaryotic genome with 16,000 tiny chromosomes [J]. PLoS Biology, 2013, 11(1): e1001473

[23]Prescott D M. The evolutionary scrambling and developmental unscrambling of germline genes in hypotrichous ciliates [J]. Nucleic Acids Research, 1999, 27(5): 1243—1250

[24]Hall T A. BioEdit: A user-friendly biological sequence alignment editor and analysis program for Windows 95/98/ NT [J]. Nucleic Acids Symposium Series, 1999, 41: 95—98

[25]Nylander J A. MrModel test v2. Uppsala University. 2004

[26]Ronquist F, Huelsenbeck J P. MrBayes3: Bayesian phylo-genetie inferenee under mixed models [J]. Bioinformatics, 2003, 19(12): 1572—1574

[27]Page R M. TreeView: An application to view phylogenetic trees on personal computers [J]. Computer Applications in the Biosciences, 1996, 12(4): 357—358

[28]Tamura B K, Dudley J, Nei M, et al. MEGA4: Molecular evolutionary genetics analysis (MEGA) software version 4.0 [J]. Molecular Biology and Evolution, 2007, 24(8): 1596—1599

[29]Yi Z Z, Huang L J, Yang R, et al. Actin evolution in ciliates (Protist, Alveolata) is characterized by high diversity and three duplication events [J]. Molecular Phylogenetics and Evolution, 2016, 96: 45—54

[30]Duharcourt S, Lepère G, Meyer E. Developmental genome rearrangements in ciliates: a natural genomic subtraction mediated by non-coding transcripts [J]. Trends in Genetics, 2009, 25(8): 344—350

[31]Chalker D L, Yao M C. DNA elimination in ciliates: transposon domestication and genome surveillance [J]. Annual Review of Genetics, 2011, 45: 227—246

EVOLUTION OF THE SCRAMBLED PATTERN OF THE ACTIN Ⅰ, A-TBP AND DNA POL A GENE WITHIN THE GENUS OXYTRICHA (PROTOZOA, CILIATES)

YANG Ran, CHEN Tian-Bing, HUANG Jun, SHI Wen-Jun, OU Shu-Qing and YI Zhen-Zhen
(Guangzhou Key Laboratory of Subtropical Biodiversity and Biomonitoring, School of Life Science, South China Normal University, Guangzhou 510631, China)

Genes in the germline (micronuclear) genome of protozoan ciliates are interrupted by multiple, non-coding sequences called internal eliminated segments (IESs) and macronuclear destined segments (MDSs). In the micronucleus, the MDSs are not arranged sequentially but scrambled for some genes. Studying scrambled gene structures will provide a better understanding of the molecular mechanisms of its evolution. In this study, we compared the complete macronuclear gene sequences and complete/incomplete micronuclear gene sequences of Actin I, α-telomere-binding protein (α-TBP), and DNA polymerase α (DNA pol α) of Oxytricha species (Oxytricha sp. (ZJ)) collected from water samples of Huguang mangrove in Zhanjiang with those of other Oxytricha species, and had 3 major findings. The scrambled patterns of the micronuclear Actin I gene in O. sp. (ZJ) were similar with those of O. nova, but shown significant differences with other species in genus Oxytricha. We revealed that the shifts of MDS-IES junctions between scrambled MDSs are not more conservative than those between non-scrambled MDSs, which is inconsistent with a previous study. The scrambled pattern of micronuclear α-TBP gene in O. sp. (ZJ) is same with that of O. nova, and the lengths of their IESs of these two species are also similar. The scrambled pattern of micronuclear DNA pol α gene in O. sp. (ZJ) is different from any of the previously reported species, and it is only somewhat similar to O. trifallax. Moreover, an IES transited to MDS was observed in DNA pol α gene, which caused the missing of its original MDS.

Oxytricha; Scramble gene; Macronuclear destined sequences (MDSs); Internal eliminated sequences (IESs); ActinⅠ; α-telomere-binding protein (α-TBP); DNA polymerase α (DNA pol α)

Q344+.1

A

1000-3207(2017)02-0285-11

10.7541/2017.35

2016-04-18;

2016-07-20

国家自然科学基金面上项目(31471973); 广东省特支计划百千万工程青年拔尖人才项目; 2015年大学生创新创业训练计划华南师范大学校级项目资助 [Supported by the National Natural Science Foundation of China (31471973); Special Support Program of Guangdong Province; 2015 SCNU Undergraduate Training Program for Innovation and Entrepreneurship]

杨然(1991—), 女, 山东济南人; 硕士研究生; 主要从事水生生物分子系统与进化研究。E-mail: yyangran@163.com

伊珍珍(1981—), 女, 山东淄博人; 副研究员, 硕导; 主要研究方向为纤毛虫原生动物分子系统与进化。E-mail: zyi@scnu.edu.cn