| 查看: 1672 | 回复: 7 | |||
[交流]
长生的仙草--PacBio第三代测序技术Nature最新文章解读
|
|||
|
Oropetium，属于禾本科，是目前已知最小基因组的草类，基因组大小为245M，基因组虽小却变化多端，具有28,466个编码蛋白基因和43%的重复序列，还有30%以上更紧凑的常染色体区域。使用常规NGS测序数据进行拼接，简直就是不可能的任务。 幸运的是，来自唐纳德植物科学中心的 Robert VanBuren, Doug Bryant以及Todd Mockler以及来自于加州大学伯克利分校等单位的研究者，2014年成为了由PacBio公司赞助的“世界上最有趣基因组项目““Most Interesting Genome in the World” SMRT Grant Program的获奖者，该项目由Pacific Biosciences公司和Sage Science共同发起。Pacific Biosciences公司承担了其中的第三代测序和基因图谱绘制工作。 重大突破： 2015年11月11日，Nature杂志在线发表了他们的文章：Single-molecule sequencing of the desiccationtolerant grass Oropetium thomaeum，基于Pacific Biosciences公司的第三代单分子测序技术独立测序，并结合Bionano单分子光学图谱技术，获得几乎完整的耐旱草Oropetium thomaeum基因组。 本文完全采用第三代单分子测序技术，基于P6-C4试剂，用32个 SMRT cells在RS II系统上运行，测序时间不到一周，产生了覆盖72×Oropetium基因组的序列。获得了Read N50>16kb的长度长，拼接得到Oropetium 基因组的99%（244Mb），共组装成625个contigs，Contig N50达到2.4Mb（而目前已发表的所有植物基因组Contig N50平均大小仅为50kb）。其中35个最大的Contigs覆盖了一半的基因组，107个Contig就已包含了90%的基因组序列。 作者轻松组装得到了Oropetium“几乎完整”的基因组草图。包括无Gap覆盖基因区域和基因间序列，比如着丝粒，端粒，转座子及rRNA簇等区域，这些区域在以往的草图中都是典型的无法拼接的区域。仅仅72x，测序的准确率就高达99.99995%，和经典的Sanger法相似，高于大部分基于NGS数据拼接的准确性。 其中全长为135,324bp的叶绿体基因组被拼接为一个单独的Contig，其中还包含两个25kb长的反向重复区域，这在以往的NGS测序拼接中都被误认为是一个拷贝。 Oropetium基因组的组装还获得了所有的18个端粒区序列。这一区域有40-900个不等的重复序列，9个着丝粒区域中的3个也被完整拼接到400kb的超长反向重复序列中，这其中包含155bp的重复单元单体，或2倍，3倍或4倍的重复单元。 Oropetium基因组的完成充分证明，采用PacBioSMRT 单分子实时测序技术和BioNano单分子基因组图谱绘制系统，可以完成前所未有的高质量基因组。 作者在文中指出，该项目从DNA提取到获得完整基因组组装只用了不到一个月，PacBio成本与基于Illumina的基因组组装也是接近的。PacBio SMRT测序可以获得全新的高质量基因组组装，并用于基因间序列的完整ENCODE类型分析，这对于其他基于NGS的方法是不可能的。 应用价值： 该物种完整基因组的完成，同时也为比较基因组学研究提供了非常有价值的资源。 禾本科（学名：Poaceae）是最主要的植物大科之一，有660属近10000种，是种子植物中最有经济价值的大科，与人们的生活息息相关，比如：小麦、稻米、玉米、大麦、高粱等，是人类粮食和牲畜饲料的主要来源。在生态上，禾本科植物适应性广，分布遍及全球，从热带至寒带，从平原到高山，湖泊、沼泽、沙漠地区均有它们的踪迹，居高等植物之冠。 上图就是文章中对禾本科的多个属，进行进化树分析，包括高粱属，稻属，短柄草属进行了比对，Oropetium基因组更为紧凑，重复区域更少，以基因组数据比对进行的种属分类就更为准确，Oropetium的优良抗逆性能也可以在跨物种作物育种中得到利用。 测序常识： 1.什么是Reads? 高通量测序平台产生的原始序列就称为reads。 2. 什么是Contig？ 拼接软件基于reads之间的overlap重叠区，拼接获得的序列称为Contig（重叠群）。Contig就是完全用序列拼接出来的，中间没有缺失，但不代表没有数量和位置上的错误。基于二代测序技术测序片段较短，很难拼出长序列，为了克服这个问题，需要使用3k/ 5k /10k /20k不同长度的文库，构建大片段，再把大片段两端100bp序列进行测序定位，这样把Contig给拼接起来。 3.何谓N50？ N50就是50%的序列中最短的那条序列长度，把序列按长到短排序，累计相加，当累计长度达到拼接基因组大小50%的时候，这一条序列是多长，那就是N50。这个对于Contig 和Scaffold适用，所以经常看到Contig N50 和Scaffold N50，这两个指标越大表示基因组组装的越好，其间的Gap越少。 4.什么是Scaffold？ 基因组de novo从头测序，通过reads拼接获得Contigs后，往往还需要构建更大的文库，以获得一定大小片段（如3Kb、6Kb、10Kb、20Kb）两端的序列。基于这些序列，可以确定一些Contig之间的顺序关系，这些先后顺序已知的Contigs组成Scaffold。对于动物来说，一般可以达到2M级别，植物很多只有几十K左右。 5.PacBio三代测序解决问题的意义？ PacBio SMRT三代单分子测序技术平均读长能达到16Kb，对于如简单重复序列，回文序列，反向重复，转座子等NGS技术发现不了或者不能定位的基因结构变异区域，直接测通，减少了拼接过程带来的错误，而且可以覆盖如高G/C含量等NGS技术的盲区，所以Contig N50就能达到之前NGS技术的Scaffold效果，甚至更好。所以如这篇文章报道的，其Contig N50就达到了2.4M。 PacBioSMRT第三代测序技术如此强悍，相信不久的未来，大量的基于PacBio三代技术的高质量基因组数据会不断涌现。 Pacbio在中国 天津生物芯片于2013年引进国内首批PacBio RS II测序仪系统，将PacBio RS II测序平台与二代测序平台相结合，目前可以达到平均8K以上的读长，它应用于小基因组完成图、动植物复杂基因组组装、表观基因组学、全长转录组测序等方面。2014年，率先升级至P6-C4试剂盒，读长高达42Kb，reads N50高达18Kb，每个SMRT Cell的数据量高达1.2Gb，碱基准确度>99.999%（QV50）。在这两年时间里，天津生物芯片通过上百个大小项目的实战经验积累，进一步优化和完善了自己基于第三代测序的数据分析平台，成为第三代测序分析服务领域领跑者。某二倍体基因组大小为350Mb左右，且杂合度比较高。天津生物芯片使用P6-C4试剂组合，测其全基因组DNA。最终测得subreads的总数据量为30.4Gb，经过过滤的subreads长度的N5O值为15450bp，subreads的平均长度为11054bp。采用TBC研发的拼接策略仅用Pacbio进行denovo拼接，最终获得了321Mb基因组序列，最终拼接结果的Contig N50值达1.35mb。 如有转载，请注明出处。 |
回复此楼» 猜你喜欢
新审查指南施行半年，AI方向的专利申请是不是更看"技术味"了？
已经有22人回复
-
专利优先审查新办法出了，哪些成果值得走快速预审通道？
已经有33人回复
-
植物学论文润色/翻译怎么收费?
已经有165人回复
-
新审查指南堵死了"一案双申"的老路，实用新型和发明只能二选一了？
已经有27人回复
-
课题快结题了才发现没东西可申请专利？四个维度帮你挖一挖
已经有22人回复
-
论文投出去之前到底要不要先申请专利？被审稿周期坑过一次后的复盘
已经有8人回复
-
横向课题结题时企业突然问专利归谁，合同里没约定的话默认算谁的？
已经有22人回复
-
教师节聊点实在的：很多老师的专利正在悄悄失效
已经有21人回复
-
INVESTIGATION OF LIFE SCIENCE 征稿（双刊号，版面费 200
已经有0人回复
-
第一次独立答复专利审查意见，才发现这活跟毕业答辩是一个套路
已经有16人回复
-
使用QE质谱仪检测杂质含量
已经有0人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 抢金币啦！回帖就可以得到:
-
项目结题验收，希望一切顺利
+1/978
-
诚征女友
+1/471
-
浙大物院表面物理课题组长期招收普博及直博学生
+2/462
-
华南理工大学 光电材料方向 王竣可教授课题组招聘博士后、博士生及科研助理
+5/430
-
坐标北京，诚征女友，奔结婚，非诚勿扰！
+5/225
-
美国肯塔基大学Dr. Sheng Tong博士后招聘（空间可控的基因组编辑与肿瘤免疫治疗）
+1/96
-
中国科学院上海有机化学研究所左智伟课题组招收2027年入学博士研究生
+1/86
-
各位大佬帮忙分析一下面上意见，是否上会？几个A、B、C？
+1/74
-
江南大学（211）招收2027年低空技术与工程、土木工程专业推荐免试研究生
+2/68
-
招收 哈工大深圳﹣大湾区大学联合培养博士生1名，2027年入学
+1/37
-
211诚招材化相关专业硕士、博士、博后（长期有效）推免-保研-调剂-申博
+1/31
-
浙江师范大学引进人才被解聘，拒绝支付50万绩效并被索赔211万
+1/26
-
深圳大学应用技术学院招聘凝聚态物理博士后
+1/16
-
深圳大学2027级光电信息工程/物理学研究生推免
+1/16
-
以色列理工-生物质、塑料等催化转化2-3名---全奖博士研究生（随时入学）和科研助理
+1/16
-
北理工集成电路杰青团队 | 诚招科助理
+1/4
-
口腔特刊客座招募客编Biofunctional Materials牙科生物功能材料专题，Scopus期刊
+1/4
-
北京理工大学-集成电路与电子学院杰青团队-招博士后
+1/4
-
北京理工大学-集成电路与电子学院杰青团队-招博士后
+1/4
-
第一次独立答复专利审查意见，才发现这活跟毕业答辩是一个套路
+1/4
8楼2015-11-24 13:31:18
简单回复
dmbb2楼
2015-11-24 10:58
回复
想飞的fish(金币+1): 谢谢参与
2015-11-24 11:03
回复
想飞的fish(金币+1): 谢谢参与
假大空4楼
2015-11-24 11:09
回复
想飞的fish(金币+1): 谢谢参与
skyish5楼
2015-11-24 11:10
回复
想飞的fish(金币+1): 谢谢参与
nono20096楼
2015-11-24 11:42
回复
想飞的fish(金币+1): 谢谢参与
·
楚云猛7楼
2015-11-24 13:13
回复
