QTL定位:MapQTL6

MapQTL是一款用于分析二倍体物种实验群体中数量性状位点(QTL)的Windows计算机程序,对于用户比较友好,更加便于使用。通过QTL分析,可以检测基因组上负责所研究数量性状表型变异的区域,估计相关的基因型效应,并以表格和(可调的)图表的形式呈现分析结果。

官网:https://www.kyazma.nl/index.php/MapQTL/

MapQTL支持的群体类型非常多,包括:BC1,F2,RIL,DH,CP,BCpxFy,IMxFy,DH2等。分析方法包含IM,MQM(=CIM)等,可以根据自己的数据需求进行选择。

没有付费的版本只可以分析两个连锁群、两个性状,也不可以导出结果,不可以复制。

1.准备输入数据

在官方的安装包里,提供了各种群体的输入文件示例,本篇文章以CP群体为例。必要的数据文件有三个,标记数据(.loc),图谱(.map)和表型(.qua)。MapQTL和joinmap是同一家开发的,文件格式通用,可以直接将joinmap中的.loc和.map带入到软件中进行计算,省去了输入文件整理的繁琐。

1.1标记数据

前四行为数据的整体介绍,“=”前是固定内容,后为自己填写的内容。正文每一行为一个标记,每一列为一个个体。


CP.loc

1.2图谱.map

.map文件中为连锁群划分的信息,按连锁群排列,第一列为标记,第二列为position。标记顺序要与.loc中一致。


CP.map

1.3表型.qua

前三行为基本介绍,同.loc,个体数目,表型数目,缺失值表示方法,每一列为一个表型,trait_1,trait_2,trait_3……,样本的顺序要与.loc一致。


CP.qua

要注意.loc和.qua中的个体排列顺序要一致,文件名前缀要一致,否则会报错。

2.导入输入数据

也不知道最多能导入多少标记,我自己的数据1w多标记,在导入最后一步.map的时候软件会崩,所以最后把标记分成了两部分,分开来跑。

①新建项目:File——New Project


New Project

②导入标记:File——Load Data,选择.loc文件


Load Data

③导入表型数据:再次File——Load Data,选择.qua文件
④导入图谱数据:第三次File——Load Data,选择.map文件

导入成功后,可以看到所有文件的基本信息:


3.QTL定位

①右键点击Populations下的性状,选择分析所有性状,此时性状下会有红色高亮。




②鼠标右键单机Maps下的连锁群信息,选中后红色高亮



③菜单下选择QTL定位方法,此处以MQM为例

④在工具栏中单击计算按钮,进行QTL定位

4.结果导出

运算完成后,点击左侧性状名称,在右侧Results下可以查看定位结果,包括连锁群、位置、标记名称、LOD值、贡献率等。点击LOD可进行升序或降序的排列,从而筛选QTL。



付费版可直接将结果导出,工具栏,导出按钮。



结果生成的图形,也可以直接进行查看,勾选想要看到的信息,然后再Charts中查看图形。

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

©著作权归作者所有,转载或内容合作请联系作者
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。

推荐阅读更多精彩内容

  • 基因定位最有效且最常用的方法就是构建遗传连锁图谱进行基因定位,该方法对于数量性状和质量性状的基因定位都适用。今天,...
    组学大讲堂阅读 7,878评论 0 24
  • Rqtl2是Rqtl的加强更新版,可以更好地处理高维数据和复杂的群体设计,对于图形输出也有了较强的改善。 0.功能...
    Wei_Sun阅读 1,541评论 0 4
  • 在上篇帖子中介绍了如何用qtl这个R包进行单QTL定位,但有很多复杂性状都是由多位点共同控制的。有些位点之间可能存...
    Wei_Sun阅读 2,184评论 0 3
  • 如果分子标记覆盖整个基因组,当标记与控制特定性状基因连锁时,不同标记基因型的表型值存在显著差异,通过数量性状观测值...
    WooWoods阅读 13,647评论 1 10
  • 作物遗传育种需要发掘对农艺性状有重要贡献的QTL。某一重要农艺性状往往是由多个QTL调控。对这些QTL进行精...
    陈牛牛阅读 5,005评论 0 29