按照T-RFLP操作流程得到数据后,需要对进行分析处理。T-RFLP数据分析软件有很多,GeneMarker是其中一个。分析T-RFLP数据的主要目的就是看有几个峰,每个峰所占的比例是多少,与Quantity One分析DGGE结果有点类似。使用GeneMarker分析T-RFLP结果的大致方法如下,比较粗略,具体各个参数的意义,请参考GeneMarker手册或相关教程:
1. 打开GeneMarker, 点”Open Data” -”Add”将待分析的T-RFLP数据打开。
2. 点击”Run”,会跳出一个对话框,里面有很多参数需要设置。其中的Size Standard一定要与T-RFLP实验时所使用的Size Standard一致,其它参数可以根据情况适当调整。官方网站上给出了下面一些参数设置建议:
Suggested Analysis Parameters
I. Template Selection
1. Panel: None
2. Analysis Type: AFLP
II. Data Process
1. Raw Data Analysis: select Smooth
2. Allele Call: deselect Auto Range and set start to 60, end to 500
3. Peak Detection Threshold: Intensity > 40
4. Stutter Peak Filter: deselect or set left and right values to 0
III. Additional settings
1. Peak Score: Reject < 0 check 1< Pass
2. AFLP-Unconfidence at Rightside Score: 1
3. 运行完毕之后就打开了GeneMarker分析窗口,可以看到样品的Gel image和electropherogram,点击左侧的样品名称可以逐个查看样品的electropherogram,通过上方的工具条可以进行放大、缩小,查看列表等一系列操作。
4. 保存数据。首先设定Report Settings,其中Report style选择Bin Table (AFLP/MLPA), Option里面可以选择Show intensity或Show Peak Area,我一般使用Show Peak Area。然后点击旁边的Save Report即可保存为excel格式文件。
5. 后续处理。 首先归一化处理(以峰面积之和中等大小的样品为基准),然后去除样品一些小的峰,多小算小?这个要根据实际情况自己设定一个值,然后合并小于1bp的相邻的峰。最后计算每个峰所占的百分比。
十一月 3rd, 2010 at 10:20
您好!
我刚学T-RFLP和DGGE,看到你的博客实在很钦佩!我实验中遇到一些问题,周围的同志没有做这方面的,所以想问问你,可以吗?
1.用GENEMARKER分析AFLP时,按照你的步骤,后来出来的表格由几种染色(blue\green\red…),结果不太一样,我该看哪种颜色的呢?
2.几种样品比较,除了1、0数据外的“?”怀疑数据,后期怎么处理?
我对软件方面的知识完全不懂,问题也可能比较幼稚,可周围又找不到相关的同学来问这些问题,呵呵
回复
Ye Lin Reply:
十一月 3rd, 2010 at 20:37
你的引物是用什么颜色标记的就看什么颜色,一般情况红色是Marker吧!
不知道你的实验是怎么做的,具体数据怎么分析,没法帮你啊
回复
十一月 3rd, 2010 at 20:46
呵呵,谢谢你!我现在大概知道颜色是怎么回事了。我想比较不同底泥样品中的微生物群落多样性。genemarker给出的bin table里,在不同片段位置处有的样品是”1″或“0”,是不是想说明这个长度的片段有的样品“有”或“没有”?那“?”呢,我怎么处理?呵呵,这个问题是不是很弱
回复
Ye Lin Reply:
十一月 3rd, 2010 at 21:04
我们一般都用峰面积,选项里面有个“Show Peak Area”
回复
十一月 21st, 2010 at 19:52
很有内涵的博客,收藏了,以后有问题请教你
回复
十二月 29th, 2010 at 15:19
想请教下,合并小于1BP的片段,您是用4舍5入得方法吗?比如121.4这类片段怎么处理,我的做法是找121 和122的,两者得到的菌我都取,请问我这做种做法可取吗?
回复
Ye Lin Reply:
十二月 29th, 2010 at 16:15
这种细节问题没有明确的规定!你采用的做法能找到参考文献就行。
回复
五月 17th, 2011 at 18:10
你好,我现在刚开始做t-rflp和dgge,很是菜希望以后多多指教能告诉我你的联系方式,qq号么
回复
叶林 Reply:
五月 17th, 2011 at 21:59
有问题在这里给我留言吧,好不好,QQ交流效率太低了。
回复
五月 31st, 2011 at 11:56
你好!能否请教一个问题,我是用Peak Scanner Software 分析数据的,但是又一次样品的size standards 的片段比实际的多一些,而且都是一些很大的数据,比如我的marker中最大的片段为625 bp,但是分析结果显示有900多 bp的片段,这是怎么回事呢?如何分析呢?还有就是同样大小的片段通过软件分析可能与不同的菌种相匹配,这怎么取舍数据呢?谢谢!
回复
叶林 Reply:
五月 31st, 2011 at 20:51
这个软件我没用过,检查一下size standards相关设置有没有错,没错的话,可能是marker被污染了吧。
同样大小的片段用这种方法根本区分不开,我觉得没法与不同的菌种建立联系。
回复
五月 31st, 2011 at 23:11
谢谢!不知道你是否通过这个网站:https://secure.limnology.wisc.edu/trflp/index.jsp
分析过T-RFLP的片段呢?
回复
叶林 Reply:
五月 31st, 2011 at 23:43
没用过,这个网站是用T-RFLP确定序列的分类单元吗?这几乎不可能的吧
回复
隋心 Reply:
七月 29th, 2011 at 21:08
你好,我现在也在分析T-RFLP数据,可是用https://secure.limnology.wisc.edu/trflp/index.jsp
老是说我数据格式错误,我想询问一下,该怎么使用这个网站,可以指导一下么?谢谢你了。期待你的回复
回复
叶林 Reply:
八月 1st, 2011 at 20:31
我看那个网站上有联系方式,你发Email问一下吧。
回复
六月 1st, 2011 at 02:41
是的,这个网站就是根据tRF,判断分类单元的,的但是往往比对出来同一个片段会有不同的结果,我现在还不知道是否把多个酶的结果同时上传会有不同。
回复
叶林 Reply:
六月 2nd, 2011 at 10:22
这样判断分类单元,发论文,应该不会被接受吧!把序列全部测出来,进行分类都可能有问题。这样只根据序列上的一个或几个酶切位点分类,可信度太差了。
回复
六月 1st, 2011 at 02:43
对了GeneMarker 这个软件可以免费下载使用吗?我刚才找了一下,没找到呀。
回复
六月 1st, 2011 at 12:46
你好!我刚刚下载了一个试用版的,参照你的说明基本能使用了,但是还有三个问题想请教:一是我目前的软件不能保存结果,所有和保存有关的键都是灰色的,什么原因呢?二是分析后给出的片段也都相差很小,只有零点几,我看到上面说合并,是指的把这些相差不多的,在1以内的合并成一个片段吗?三是数据表格中的0$0, 1$一个数值,这代表的什么意思?0$0是不是指样品中没有这个峰,而有一的表示有这个峰,峰面积是多少。对吗!
非常感谢你的帮助!
回复
叶林 Reply:
六月 2nd, 2011 at 10:19
这个软件不是免费软件,试用版不能保存结果很正常。我们实验室买了一个License。
这个方法我好久没用过了,以后也不打算再用了,一些具体的细节问题不记得了,你研究一下软件的说明书吧
片段长度在1以内的合并,这个我记得,我们当初就是这样做的。
回复
六月 2nd, 2011 at 22:28
还想请教:根据你的经验,你觉得T-RFLP最大的问题在哪呢?谢谢!
回复
叶林 Reply:
六月 2nd, 2011 at 22:43
有些不同的细菌,在酶切的时候,要么都有相同的切点,要么都没有切点,找不到一种合适的酶将它们区分开。
回复
六月 2nd, 2011 at 22:44
你们现在主要用什么方法呢?454吗?
非常感谢你的帮助!
回复
叶林 Reply:
六月 2nd, 2011 at 22:49
嗯,454和illumina的高通量测序
不用客气,呵呵
回复
六月 2nd, 2011 at 22:47
如果采用多个酶分别进行,然后一起分析呢?你觉得这样是否会好一些呢?
回复
叶林 Reply:
六月 2nd, 2011 at 22:53
会有一点点改善,但不能从根本上解决问题,我以前专门写了个程序试过上百个酶,也没找到办法把我样品中的不同序列分开。
还是测序最直接,序列测出了,什么都清楚了。
回复
六月 2nd, 2011 at 22:59
是呀,还是测序最直接,但是如果没有较多的资助,那还是有困难的,真是不能两全呀!
真的非常高兴能搜到你的博客,你是做环境微生物的吗?我昨天找到你的一篇文章,还没有仔细看看呢,一定好好学习一下,有什么问题再请教。
genemarker 这个软件有破解版吗?
回复
叶林 Reply:
六月 3rd, 2011 at 22:39
我不知道有没有破解版
回复
七月 11th, 2011 at 17:03
对于第五步“然后合并小于1bp的相邻的峰”。
请问合并后的峰面积是取和还是取平均?
例如:
菌种 点位(Hha I) 峰面积(HhaI) 点位(HaeIII) 峰面积(Hae III)
Craurococcus(D85828) 59 10667 71.2 1869
Craurococcus(D85828) 59 10667 72.9 1924
用双酶切来分析得出在Hha I的同一个位点,Hae III的不同位点对应有一种菌,如果这两行结果要合并,那么HaeIII的峰面积是取1869和1924之和还是它们的平均呢?理由是?
回复
八月 15th, 2011 at 18:46
我是南大的硕士,现在手里有一批trflp数据(由size,pk high,pk area 等字段组成的excel.我第一步要怎么走啊?
回复