24小时热门版块排行榜    

查看: 1207  |  回复: 0

jaxb

银虫 (正式写手)

[交流] 利用基于半导体的第二代测序技术高通量的发掘兔基因组SNP

利用基于半导体的第二代测序技术高通量的发掘兔基因组SNP
意大利博洛尼亚大学F. Bertolini……L. Fontanesi*
虽然兔基因组测序和组装(oryCun2.0)已经完成,而且也有一些关于兔微卫星、SNP和CNV的研究,但是兔基因组大规模多态性发掘还没有报道。因此本研究利用离子流个人基因组测序仪(PGM)测序了HaeIII
and RsaI两个限制性酶的简化代表库(RRLs),识别了10只家兔(1只法国香槟兔,3只巨型格仔兔,1只勃艮第兔,2只维兰特兔和3只商业肉兔)基因组的SNPs。经过ION TORRENT suite v.2.2软件的过滤,两个库分别获得了2,917,879和4,046,871个reads,平均长度分别为96bp和103bp,总共长度分别为280.51Mb(248.49Mb的质量>20)和417.28Mb(360.89Mb的质量>20)。其中分别有2,627,434(90%)和3,685,226(91%)个reads比对到了兔基因组,占兔基因组的7.01%和9.06%(合并后为15.82%)。序列数据已经上传到欧洲核苷酸数据库(EMBL),序列号为ERP002438。实验要求至少有4个reads以上覆盖的SNP为高可信度SNP,因此只能在31,389,157bp和58,542,265bp兔基因组长度上寻找,其SNP的覆盖率分别是5.38和6.68。两个库的数据合并覆盖了兔基因组上107,922,811bp长度(<4reads的也算),深度为5.7X。本研究对SNP进行过滤,总共推断出62,491个SNPs,平均每1727bp一个SNP(62,491/107,922,811)。其中8号染色体和2号染色体SNP最多,分别为6,289和4,928个。SNP的密度最高的也是18号染色体,为11.1kb每SNP,最低在14号染色体,为50.8kb。研究又用桑格测序进行了验证,检测到了22个SNPs(OCU7上9个,OCU18上12个,OCU19上1个),其中16个SNPs在库中,包括一个reads小于4X的。也就是6/22的假阴性率和0/22的假阳性率。也就是说68%的高可信度SNPs(15个)检出率和100%假阳性率。大部分的SNPs都存在于基因间(67.3%)和内含子区域(32.5%)。位于编码区域的SNPs中有538个是同义突变,有479个是错义突变,有16个是无义突变。本研究结果为以后该物种的应用和基础研究提供了非常有用的信息。

交流与问题:
1,为什么选择reads的覆盖率为4X?从文中可以看到假阴性的原因主要就是SNP对reads的限制太严格。作为仅是对SNP的发掘,应该是可以有适当的假阳性,尽量减少假阴性。
2,桑格测序进行验证的时候,区域是如何选取的?可以这样理解吗,两个库中在这些区域上高可信度的SNPs为15个,且全部被桑格测序所发现,还有桑格测序发现的6个没有存在于库中?
3,三种突变合起来的SNP是1033个,占所有发现的1.65%。这与前面基因间和内含子区域的SNPs占99.8%有点不符合,是什么原因?
4,总得SNP密度为1727bp一个SNP,怎么密度最高的8号染色体密度反而是11.1kb一个SNP呢?是不是小数点打错了啊?
5,基因间区域和内含子区域怎么区别?基因间区域( intergenic)就是不编码区域吧,不属于内含子区域?还是说内含子只是指基因内外显子之间的区域?

来源:http://www.ncbi.nlm.nih.gov/pubmed/24444082-Bertolini F, Schiavo G, Scotti E, Ribani A, Martelli PL, Casadio R, Fontanesi L. High-throughput SNP discovery in the rabbit (Oryctolagus cuniculus) genome by next-generation semiconductor-based sequencing. Anim Genet, 2014, 45(2):304-7.
回复此楼
丰富你的动物食品!
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

相关版块跳转 我要订阅楼主 jaxb 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 时间戳变了,能看出什么问题? +18 基诺咪客 2026-08-17 23/1150 2026-08-20 17:19 by Godzela
[基金申请] 只有每年这种时候来逛逛小木虫 +21 yaoyewhu2008 2026-08-20 22/1100 2026-08-20 17:03 by yaoyewhu2008
[基金申请] 今天系统多次维护,明天很可能放榜! +9 zju2000 2026-08-16 10/500 2026-08-20 15:56 by Leogzhya
[基金申请] 放榜前的不淡定 40+3 snowwithsea 2026-08-19 12/600 2026-08-20 15:56 by lfy8008
[基金申请] 言之凿凿,明天官网会通知放榜? +14 医学老男孩 2026-08-20 16/800 2026-08-20 15:47 by LNP@mRNA
[基金申请] 重要消息,中午系统在维护 +11 yuleib84 2026-08-18 12/600 2026-08-20 11:09 by xskun
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +16 Ldrop2023 2026-08-13 20/1000 2026-08-20 11:07 by xskun
[基金申请] 今天基金会出结果吗?20260819 +14 kkkl_v 2026-08-19 15/750 2026-08-20 08:21 by sunzitan
[基金申请] filecode,4个jtjc了 +11 ziyangfang 2026-08-19 13/650 2026-08-19 21:28 by aasahr
[基金申请] 今天放榜没戏了吧 +4 yuleib84 2026-08-19 5/250 2026-08-19 19:52 by hhs666
[基金申请] 今天放榜吗? +14 布布和一二 2026-08-19 15/750 2026-08-19 18:07 by gltch
[基金申请] 2027广东省杰青 +4 奶牛小黑 2026-08-15 10/500 2026-08-19 11:02 by wanfengnew
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +6 Tide man 2026-08-14 7/350 2026-08-19 09:56 by ZJTJZ
[基金申请] 明天放榜? +5 Shxjjxjkx 2026-08-18 5/250 2026-08-18 18:14 by -大大大大大-
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 6/300 2026-08-18 15:36 by wwm09
[基金申请] 时间戳又变了8-15 +14 archvillain 2026-08-15 26/1300 2026-08-18 13:37 by phantomgost
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
信息提示
请填处理意见