| 查看: 3879 | 回复: 16 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
用Trinity进行de novo拼装
|
|||
| 最近在公司做了转录组测序，植物样本，illumina Hiseq2000，三个样本共100M的clean reads，一起拼接。该公司用Trinity软件进行de novo拼接，拼出来contigs达40多万条，后来又用CD-HIT进行聚类，结果仍有30多万条，其中包含很多转录本的存在。从文献中看一般植物转录组de novo拼接也就几万条序列，太多的转录本会影响到后面表达量的比较。不知道是由于Trinity这个拼接软件的原因还是其他什么原因。 |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有4人回复
-
售SCI文章，我:8O.5.5.1O.54,科目全,可十急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有6人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有6人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有8人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有10人回复
-
2026国自然函评费到账
已经有11人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有3人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- Trinity拼装contig数量
已经有0人回复
- 陈连福NGS生物信息分析 2015暑期培训班招生简章
已经有326人回复
- 模具咬花数字化：更简单更快地实现丰富的表面质感
已经有0人回复
- 单词演义（五）
已经有40人回复
- 关于de novo拼接结果的疑问
已经有1人回复
- 【迅雷】1024分辨率《加勒比海盗4》BD中英双字无水印
已经有6人回复
- 【迅雷】2011动作科幻大片《加勒比海盗4》BD国语配音中字1024高清
已经有5人回复
- 【迅雷】2011最新动作冒险大片《加勒比海盗4：惊涛怪浪》DVD中英双字
已经有65人回复
» 抢金币啦！回帖就可以得到:
-
告别流水线爽文！108万字完结仙侠《女娲石之灵石奇缘（仙界篇）》值得细品
+1/9020
-
散金祈福
+5/4015
-
南方医科大学|生物材料与医学工程团队诚聘科研助理
+2/270
-
散金求好运
+1/192
-
诚征另一半
+1/172
-
那年春暖花开，今朝如梦初醒 —— 北漂十五载，相亲十年，归乡孑然一身
+1/98
-
山东征女友，坐标济南
+1/89
-
焦虑
+1/66
-
上海交通大学张航课题组招聘博士后（电化学能量存储与转换）
+1/43
-
瑞典-博后-固态电池-粘结剂-锂离子电池方向
+1/37
-
南京理工大学优青团队催化化学方向招收推免研究生
+1/26
-
***大叔找对象（沪深等地）
+1/23
-
pcold 原核表达载体
+1/9
-
祈祷老公中国基
+1/7
-
南方医科大学深圳医院陈烨/余涛教授课题组联合培养博士后招聘
+1/6
-
2026年国自然
+1/3
-
找女友
+1/2
-
新审查指南施行半年，AI方向的专利申请是不是更看"技术味"了？
+1/2
-
出全新蔡司 Plan?APOCHROMAT 20×/0.8 平场复消色差物镜
+1/2
-
找女友
+1/1
★ ★ ★ ★ ★
小木虫: 金币+0.5, 给个红包，谢谢回帖
zhaohq1209: 金币+4, 信息学的GG就是不一般 2012-04-13 11:28:43
小木虫: 金币+0.5, 给个红包，谢谢回帖
zhaohq1209: 金币+4, 信息学的GG就是不一般 2012-04-13 11:28:43
|
我看到你的测序结果后感到很奇怪，有几个问题想问你： （1）你的三个样本的reads长度是多少啊？如果reads太短，也会导致拼接的结果变化。还有就是你用哪种方式建库，single-end，mated-paired 和 paired-end？一般，如果用paired-end技术，可能拼接结果会更好点。 （2）你的三个样本是混在一起建库后测序的吗？还是分别建成三个库，分别加三个不同barcode的加以区别吗？如果是后者，我有一个问题，一般Hiseq 2000的一个lane可以获得5-30G左右的数据，而且一个lane里面最多也就可以加到24样品啊，一般的公司在一个lane最多也就加到8个样品，所以你的一个样品获得数据5G/24=210Mb左右才对啊，三个样品应该是至少也应该是600Mbp左右啊！如果你的reads总共才100Mbp的话，拼接成这样，应该是因为你的测序depth太低（因为植物的基因很大，除了拟南芥小点（125Mb）外，一般都大于400Mbp）。你的植物的物种是木本，还是草本植物；如果是木本植物，它的基因组可能会更大点，也会影响你的测序depth。 （3）在使用Tiniity拼接时，你输出的最小的contig长度是多少呢（即“--min_contig_length”参数设置为多少，程序默认200bp）？ （4）在使用Tiniity拼接时，他使用哪个方法拼接的：Inchworm、 Chrysalis 和Butterfly （A）Inchworm assembles the RNA-seq data into the unique sequences of transcripts, often generating full-length transcripts for a dominant isoform, but then reports just the unique portions of alternatively spliced transcripts. （B） Chrysalis clusters the Inchworm contigs into clusters and constructs complete de Bruijn graphs for each cluster. Each cluster represents the full transcriptonal complexity for a given gene (or sets of genes that share sequences in common). Chrysalis then partitions the full read set among these disjoint graphs. （C）Butterfly then processes the individual graphs in parallel, tracing the paths that reads and pairs of reads take within the graph, ultimately reporting full-length transcripts for alternatively spliced isoforms, and teasing apart transcripts that corresponds to paralogous genes. 三种方法得到的结果也是有所不同的。 |
8楼2012-04-13 09:38:17
2楼2012-04-11 16:50:53
3楼2012-04-11 16:57:55
4楼2012-04-12 08:26:10
