EVOLUTION OF THE SCRAMBLED PATTERN OF THE ACTIN Ⅰ, Α-TBP AND DNA POL Α GENE WITHIN THE GENUS OXYTRICHA (PROTOZOA, CILIATES)
-
摘要: 研究旨在对尖毛虫属内现有物种的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的转变，使得本应删除的序列成为基因组永久保留的一部分。Abstract: 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.
-
纤毛虫原生动物具有二型核, 小核基因通常被内部删除序列(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 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 Ⅰ gene
MDSs are represented in grey, IESs are represented in black; MDS2 is inverted; white box represents unavailable gene regim; the same applies below
![]()
图 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
![]()
图 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
![]()
图 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)
![]()
图 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
![]()
![]()
![]()
图 8 Oxytricha sp. (ZJ) 进化过程中IES向MDS的转变
Figure 8. IES to MDS transformation during Oxytricha sp. (ZJ) evolution
表 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 
下载: 导出CSV
表 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
下载: 导出CSV
表 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
下载: 导出CSV
表 4 Oxytricha sp. (ZJ) 小核Actin Ⅰ基因的MDSs、IESs及指针序列
Table 4 Characteristics of MDSs, IESs, and pointer sequences in Oxytricha sp. (ZJ)Actin Ⅰ gene
MDSn 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
下载: 导出CSV
表 5 Oxytricha sp. (ZJ) 小核α-TBP基因的MDSsIESs及指针序列
Table 5 Characteristics of MDSs, IESs, and pointer sequences in Oxytricha sp. (ZJ)α-TBP gene
MDSn 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 clones
下载: 导出CSV
表 6 Oxytricha sp. (ZJ) 小核DNApol α基因片段MDSs, IESs及指针序列
Table 6 Characteristics of MDSs, IESs, and pointer sequences in Oxytricha sp. (ZJ) DNApol α gene segment
MDSn 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 information
下载: 导出CSV
-
[1] Prescott D M, Greslin A F. Scrambled Actin Ⅰ gene in the micronucleus of Oxytricha nova[J]. Developmental Genetics, 1992, 13(1):66-74 doi: 10.1002/(ISSN)1520-6408
[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 http://adsabs.harvard.edu/abs/1999NYASA.870..349C
[3] Nowacki M, Landweber L F. Epigenetic inheritance in ciliates[J]. Current Opinion in Microbiology, 2009, 12(6):638-643 doi: 10.1016/j.mib.2009.09.012
[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 http://www.pnas.org/content/86/16/6264.full.pdf
[5] Prescott D M. Genome gymnastics:unique modes of DNA evolution and processing in ciliates[J]. Nature Reviews Genetics, 2000, 1(3):191-198 doi: 10.1038/35042057
[6] Katz L A, KovnerA M. Alternative processing of scrambled genes generates protein diversity in the ciliate Chilodonella uncinata[J]. Journal of Experimental Zoo-logy Part B Molecular and Developmental Evolution, 2010, 314(6):480-488 doi: 10.1002/jez.b.21354
[7] Lauth M R, Spear B B, Heumann J, et al. DNA of cilia-ted protozoa:DNA sequence diminution during macronuclear development of Oxytricha[J]. Cell, 1976, 7(1):67-74 doi: 10.1016/0092-8674(76)90256-7
[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 http://genesdev.cshlp.org/content/6/5/788.full.pdf
[9] Hoffman D C, Prescott D M. The germline gene enco-ding DNA polymerase α in the hypotrichous ciliate Oxytricha nova is extremely scrambled[J]. Nucleic Acids Research, 1996, 24(17):3337-3340 doi: 10.1093/nar/24.17.3337
[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 Ame-rica, 1995, 92(9):3888-3892 doi: 10.1073/pnas.92.9.3888
[11] DuBois M L, Prescott D M. Volatility of internal elimi-nated segments in germ line genes of hypotrichous cilia-tes[J]. Molecular and Cellular Biology, 1997, 17(1):326-337 doi: 10.1128/MCB.17.1.326
[12] Hogan D J, Hewitt E A, Orr K E, et al. Evolution of IESs and scrambling in the Actin Ⅰ gene in hypotrichous cilia-tes[J]. Proceedings of the National Academy of Scien-ces of the United States of America, 2001, 98(26):15101-15106 doi: 10.1073/pnas.011578598
[13] Dalby A B, Prescott D M. The scrambled Actin Ⅰ gene in Uroleptus pisces[J]. Chromosoma, 2004, 112(5):247-254 doi: 10.1007/s00412-003-0270-4
[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 doi: 10.1007/s00239-005-0166-4
[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 doi: 10.1007/s004120050311
[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 doi: 10.1093/molbev/msj089
[17] Hoffman D C, Prescott D M. Evolution of internal elimina-ted segments and scrambling in the micronuclear gene encoding DNA polymerase alpha in two Oxytricha species[J]. Nucleic Acids Research, 1997, 25(10):1883-1889 doi: 10.1093/nar/25.10.1883
[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 doi: 10.1073/pnas.97.7.3298
[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 doi: 10.1078/143446104774199628
[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 doi: 10.1073/pnas.0507682102
[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 doi: 10.1371/journal.pbio.1001473
[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 doi: 10.1093/nar/27.5.1243
[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 phylogenetie inferenee under mixed models[J]. Bioinformatics, 2003, 19(12):1572-1574 doi: 10.1093/bioinformatics/btg180
[27] Page R M. TreeView:An application to view phylogene-tic trees on personal computers[J]. Computer Applications in the Biosciences, 1996, 12(4):357-358 http://bioinformatics.oxfordjournals.org/content/12/4/357.full.pdf
[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 doi: 10.1093/molbev/msm092
[29] Yi Z Z, Huang L J, Yang R, et al. Actin evolution in cilia-tes (Protist, Alveolata) is characterized by high diversity and three duplication events[J]. Molecular Phylogenetics and Evolution, 2016, 96:45-54 doi: 10.1016/j.ympev.2015.11.024
[30] Duharcourt S, Lepère G, Meyer E. Developmental geno-me rearrangements in ciliates:a natural genomic subtraction mediated by non-coding transcripts[J]. Trends in Genetics, 2009, 25(8):344-350 doi: 10.1016/j.tig.2009.05.007
[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 doi: 10.1146/annurev-genet-110410-132432
-
期刊类型引用(1)
1. 何燕，尹家奇，生欣. 艾美游仆虫大核基因组与转录组测序及结构特征. 水生生物学报. 2021(05): 1014-1023 . 
本站查看
其他类型引用(0)
计量
- 文章访问数: 1418
- HTML全文浏览量: 532
- PDF下载量: 12
- 被引次数: 1