QTL上位性分析:QTLNetwork2.0

官网:http://ibi.zju.edu.cn/index.html/BCL/software/qtlnetwork/download.html

QTLNetwork是一款用于QTL上位性分析,以及QTL与环境互作(QE)效应分析的Windows版软件,适用于DH,RI,BC,F2,IF2和BxFy等群体,同时对结果进行可视化。

1.输入数据

QTLNetwork2.0输入数据包括两部分:标记连锁图谱文件(map file)和数据文件(data file)。图谱文件包含所有连锁群上的标记,以及他们的顺序和遗传距离;数据文件包含每个个体的标记,和表型值。软件下载成功后,随软件带有示例数据文件 (\SampleData) ,可以参考示例文件的格式进行整理。

1.1标记连锁图谱文件(map file)格式

在文件的最开始是一些一般说明,最前面都以下划线“_”开始:


General Description
  • _DistanceUnit:遗传距离单位, “cM” stands for centi-Morgan and “M” stands for Morgan
  • _MapFunction:图谱的创建方式,“K” is for Kosambi function and “H” for Haldane function
  • _Chromosomes:图谱中染色体或连锁群的总数目
  • _MarkerNumbers:每条染色体或连锁群上的标记数目,与连锁群一一对应

一般说明的下面是图谱文件的主体内容,以* MapBegin * 开始,* MapEnd *结束,中间是标记的顺序和位置,这里的位置是两个标记之间的距离。


Map Body

1.2数据文件(data file)文件格式

数据文件主要由三部分内容组成:一般说明,标记数据和表型数据。

a.一般说明
General description
  • _Population:群体类型,常见的有RI,BC,F2,IF2,BxFy等
  • _Genotypes:群体中的基因型总数
  • _Observations:每个性状的观测值总数目。例如200个个体,两次重复,此处填400,没有重复则为200,与表型的个数无关。
  • _Environments:试验设计中是否有多个环境,yes或者no
  • _Replications:试验设计中是否有重复或区组,yes或者no
  • _TraitNumber:性状的总数目
  • _TotalMarker:标记数目的总数,这个数值必须与map file中 _MarkerNumbers 的总数一致
  • _MarkerCode:定义一个标记编码方案。P1:标记与P1一致,P2:标记与P2一致,F1:标记与F1一致,F1P1:标记不与P2一致,F1P2:标记不与P1一致。在文件中每个标记用一个字母或一个数字表示,缺失值用“.”表示。
b.标记数据

这部分以* MarkerBegin * 开始,* MarkerEnd *结束,标记顺序要与map data中一致,支持两种格式。

第一种每一行为一个个体,每一列为一个标记,每行的末尾以“;”结束:


Type I

第二种每列行为一个个体,每一行为一个标记,每行的末尾以“;”结束:


Type II

标记的名字和标记顺序必须与map data中一致,标记名字中不可以有分隔符。

c.表型数据

这部分以* TraitBegin * 开始,* TraitEnd *结束,如果存在多环境和重复,需要将不同环境、不同重复下的表型值均展示出来,以下图为例,第一列为环境,第二列重复,第三列基因型,当没有多环境多重复时,这几列可以省略不写。第四列开始为不同性状的观测值,每一行以“;”结尾:


Trait data body

2.创建新项目

2.1输入数据

在完成输入数据准备后,准备创建新的项目。
Project——New(QTLNetwork Data Format):在对应位置选择map文件和data文件。


image.png

输入数据导入成功:


2.2导入为其他数据格式

QTLNetwork支持将输入成功的数据转出为QTLNetwork或者QTL Cartographer数据格式,Project——Export Source Data:


2.3开始计算

在数据导入成功,检查无误后,开始计算。单击Run,设置计算参数:


  • Map Epistasis:只显示单个QTL的效应和上位性
  • High Order Epistasis :显示AD, DA, DD等高阶上位性,这个选项只适用于F2和IF2群体。
  • Do 2D Genome Scan:选择这个选项来绘制有或没有单位点效应的上位性QTL。否则,程序将只检测具有单位点效应的QTL间的上位相互作用。
  • Permutation:设置置换检验次数来控制一类错误。
  • MCMC:通过MCMC( Monte Carlo Markov Chain method)方法来估计QTL效应,否则使用的是MIM(mixed linear model)方法。
  • Superior Genotype Prediction:根据QTL效应估计优势基因型。
基因组扫描配置设置

包含窗口大小、步长等。


显著性配置设置

运行时,会有计算进度条:



计算结束后自动跳转计算结果:



图中的颜色以及图形解释说明书中有详细介绍:

在工具栏setting中,可以设置图形参数:

3.保存图表及计算结果

完成所有计算后,可以导出图形和报告,工具栏中Edit——save picture和Edit——Save report(.pre)。



报告中分为部分:
第一部分是变异组成,基因型变异,环境变异,基因型环境互作效应,重复、残差、表型变异等;
第二部分是QTL效应,包含QTL的位置、加性效应、显性效应等;
第三部分是QTL的遗传力,加性效应遗传力、显性效应遗传力等;
第四部分是上位性检测,两个QTL的位置,以及相互作用;
第五部分是上位性遗传力;
第六部分是基因型值,包括双亲,F1和估计优势基因型的基因效应;
第一部分是估计优势基因型。

引用转载请注明出处,如有错误敬请指出。

最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 193,812评论 5 459
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 81,626评论 2 371
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 141,144评论 0 319
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 52,052评论 1 263
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 60,925评论 4 355
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 46,035评论 1 272
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 36,461评论 3 381
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 35,150评论 0 253
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 39,413评论 1 290
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 34,501评论 2 307
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 36,277评论 1 325
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 32,159评论 3 312
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 37,528评论 3 298
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 28,868评论 0 17
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 30,143评论 1 250
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 41,407评论 2 341
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 40,615评论 2 335

推荐阅读更多精彩内容