使用MatrixEQTL进行cis/trans-eQTL分析

最新推荐文章于 2023-03-30 10:20:43 发布
VIP文章 最新推荐文章于 2023-03-30 10:20:43 发布
阅读量3.8k 收藏 14
点赞数 2
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/weixin_43569478/article/details/108079736
版权

欢迎关注”生信修炼手册”!

Matrix是一款经典的eQTL分析软件，可以支持cis和trans-eQTL的分析，官网如下

http://www.bios.unc.edu/research/genomic_software/Matrix_eQTL/

其运行速度快，官方给出的和其他工具的运行速度比较结果如下

支持线性回归和ANOVA等模型，同时支持协变量的矫正。该软件有对应的R包，使用简便，基本用法如下

1. 读取snp相关文件

对于SNP而言，有以下两种文件

  1. SNP分型结果

  2. SNP染色体位置

SNP分型结果文件内容示意如下

该文件读取的代码如下

SNP染色体位置文件的内容示意如下

该文件读取的代码如下

2.读取gene相关文件

对于gene而言，有以下两种文件

  1. 基因表达量结果

  2. 基因染色体位置

基因表达量结果文件内容示意如下

该文件读取的代码如下

基因染色体位置结果文件内容示意如下

该文件读取的代码如下

3. 读取协变量

协变量是可选的，内容示意如下

该文件读取的代码如下

4. 进行eQTL分析

代码如下

运行完成之后，查看对应的输出结果即可，cis和trans的输出结果是类似的，示意如下

cis和trans本质上就是通过SNP位点和gene之间的距离进行区分，首先对所有的SNP-gene对进行计算，然后根据距离阈值将结果进行拆分，小于阈值的SNP-gene对归于cis-eQTL, 大于阈值的归于trans-eQTL。

·end·

—如果喜欢，快分享给你的朋友们吧—

往期精彩

  GWAS meta分析

  基因型填充

最低0.47元/天 解锁文章
生信修炼手册
关注
  • 2
    点赞
  • 踩
  • 14
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
《微信小程序进阶实战之分答应用开发（中级项目）》（完整版）
tt15050的博客
03-14 568
找到Matrix eQTL这个包，看下文章Matrix eQTL: ultra fast eQTL analysis via large matrix operations（https://doi.org/10.1093/bioinformatics/bts163 ） eQTL（表达数量性状位点）计算transcript-SNP 的关系，即分析SNP与基因的表达是否相关。由于计算数量巨大，很多人都...
eQTL资料阅读（一）
sinat_40922336的博客
05-24 3667
文章目录Expression quantitative trait loci: present and future1、简介2、Cis- trans-区别：eQTL原理及程序介绍ppt1、cis-和trans-区别2、ppt介绍内容What eQTL areHow linear models can be used to find themComplications that arise and possible approaches to deal with themHow to carry out la
Die Kristallstruktur von cis- und trans-N-iso-Propylamidodimethylindium, [(CH3)2In-N(H)iC3H7]2
06-29
Die Kristallstruktur von cis- und trans-N-iso-Propylamidodimethylindium, [(CH3)2In-N(H)iC3H7]2 Die Kristallstruktur von cis- und trans-N-iso-Propylamidodimethylindium, [(CH3)2In-N(H)iC3H7]2 F. Cordeddu, H.-D. Hausen† und J. Weidlein* Stuttgart, Institut für Anorganische Chemie der Universität Bei der Redaktion eingegangen am 26. November 2001. The Crystal Structure of cis- and trans-N-iso-Propylamidodimethyl Indium, [(CH3)2In- N(H)iC3H7]2 Abstract. According to the X-ray structure determin
eqtl-intro:eQTL分析入门课程
06-16
eQTL分析简介 该存储库包含介绍性 eQTL 分析课程的材料。 Docker 镜像 泊坞窗目录包含Dockerfile和数据需要建立一个包含练习的所有必要组件的码头工人形象，这门课程的组成部分。 生成的 Docker 映像可从。 docker pull humburg/eqtl-intro将拉取最新版本的镜像。 安装 Docker 使用这个镜像需要一个有效的 Docker 安装。 、 、 和许多安装说明可通过 Docker。 使用 Docker 镜像 docker 镜像包含和并允许使用 RStudio IDE 和 R 控制台。 要使用 IDE，请使用以下命令启动 Docker 容器 docker run -p 8787:8787 humburg/eqtl-intro 然后，通过访问端口 8787 上的 docker 主机，可以通过 Web 浏览器访问 RStudio。在 Linu
co-trans-ext 翻译工具-crx插件
04-04
语言:中文 (简体) 集搜狗翻译、百度翻译、有道翻译、金山词霸、谷歌翻译于一体的翻译工具 co-trans-ext 是一个集搜狗翻译、百度翻译、有道翻译、金山词霸、谷歌翻译于一体的翻译扩展。各平台可随时切换，取长补短，更适合于阅读各领域的文档。不止于翻译，解析详细，更利于对各门语言的深入学习与理解。支持划词翻译和网页翻译。
Light-Trans.rar_Trans-light
07-15
图形图像亮度变换，具有实际例子，使用matlab编写，不容错过哦！
eQTL hotspot : eQTL热点区域
庐州月光的博客
12-16 1541
欢迎关注”生信修炼手册”!众所周知，eQTL分为了cis-eQTL和trans-eQTL两种作用方式，cis模式下只能调控临近的基因，而trans模式突破了距离限制，在该模式下一个eQT...
值得借鉴的eQTL可视化形式
庐州月光的博客
12-15 2525
欢迎关注”生信修炼手册”!eQTL的分析结果本质就是一些调控基因表达的SNP位点，在结果展示时，最经典的可视化方式如下按照SNP分型结果对样本进行分组，用箱体图或者小提琴图的方式展示不同...
使用FastQTL进行cis-eQTL分析
庐州月光的博客
12-10 3171
欢迎关注”生信修炼手册”!FastQTL是一款专门用于cis-eQTL分析的软件，在GTEx项目中就是采用该软件进行cis-eQTL的分析，对应的文章发表在Bioinformatics杂...
QML 基础类型 vector 和 matrix
gooderfeng的专栏
05-28 1807
vector2D vector3D vector4D 和 matrix4x4 是QML中与opengl相关的一些基础类型。 vector2D vector3D vector4D几乎具有一致的借口和函数，主要用于保存一个点的信息。 matrix4x4 用来表示具有16个元素的数组，主要用来做opengl中的matrix计算，如修改坐标系，旋转，移动等操作。 vector*D 常用的函数
TRANS-JOUR.zip
02-09
TRANS-JOUR.zip
Rexx/Trans-开源
04-25
Rexx / Trans提供了启用Rexx的应用程序和Rexx解释器之间的抽象级别。 这允许使用SAA Rexx API编写应用程序并与Rexx / Trans链接。 在运行时，Rexx / Trans会为您找到合适的Rexx解释器
Science2014一篇关于eQTL的算法描述
08-25
关于eQTL的算法，进行线性回归.生物类，生物信息方向
GTEx:来自 GTEx 数据的 eQTL 检测
06-18
GTEx 祖松鹏通过贝叶斯非参数方法从 GTEx 数据中检测顺式 eQTL。 程序 给定距离截断值，获取每个基因周围的 SNP（文件：genelocsnp） 运行 R 脚本以获取 DECODE 结果。
解决R速度太慢问题
weixin_30439067的博客
08-15 2805
R的速度慢一直被人诟病，最近做一个比较大的dataset的分析，跑得实在太慢，发现症结是R的data frame的index太慢： 以下为测试： gene_list = 1:100000 eQTL_mat = matrix(nrow = length(gene_list), ncol = 7) # 创建一个matrix eQTL_df = as.data.frame(matrix(nr...
表达数量性状位点（eQTL）的概念及其相关分析原理
weixin_38097914的博客
03-30 1113
第三个是基因型数据，同样地，除开第一列，其它列都代表不同的样本，每一行代表的是不同的基因型（SNP），一般基因型数据用0，1，2这三个数字编码，代表的是效应等位基因剂量。这里gene1（因变量）一般就是一个基因的表达量，snp1（自变量）就是一个SNP的基因型，两者拟合，矫正相关干扰项（如sex和age等），error_term是指回归模型的误差项。第二个是基因表达量信息，同样地，除开第一列，其它列都代表不同的样本，每一行代表的是不同的基因（一般来说基因表达数据需要先进行标准化转换）。
eQTL
an05423833476591的博客
08-11 634
首先QTL是数量性状位点，比如身高是一个数量性状，其对应的控制基因的位点就是一个数量性状位点，而eQTL就是控制数量性状表达位点，即能控制数量性状基因（如身高基因）表达水平高低的那些基因的位点。 数量性状基因座：控制数量性状的基因在基因组中的位置称数量性状基因座。常利用DNA分子标记技术对这些区域进行定位，与连续变化的数量性状表型有密切关系 表达数量性状基因座（expressi...
ncRNA-eQTL:eQTL与ncRNA的碰撞
庐州月光的博客
12-11 581
欢迎关注”生信修炼手册”!ncRNA-eQTL数据库专注于研究不同肿瘤中调控ncRNA表达量的eQTL, 通过TCGA数据库获取不同肿瘤中的SNP分型信息，以及lncRNA和miRNA的...
全基因组关联分析学习资料（GWAS tutorial）20210313更新版
wendy_milk的博客
03-14 6115
时隔一年多，于今日（20210313）再次更新教程，后续如果新出教程的话会一直更新。 同时感谢多年陪伴我成长、一直默默关注我的你们~ by 橙子牛奶糖（陈文燕） 1.前言 很多人问我有没有关于全基因组关联分析（GWAS）原理的书籍或者文章推荐。 其实我个人觉得，做这个分析，先从跑流程开始，再去看原理。 为什么这么说呢，因为对于初学者来说，跑流程就像一个大黑洞，学习原理就像一个小黑洞。 很多人花了好几个月的时间在看原理，一旦丢给他数据去分析，依旧束手无策。 不会跑流程，内心依旧会很恐慌。就像从来没有入门一样。
feign.FeignException$InternalServerError: [500] during [GET] to [http://witmed-system-provider/trans/system/dictionary/department/clinic/list] [DepartmentTransport#getClinicDeptList()]: [{"timestamp":"2023-07-24T02:26:51.783+00:00","status":500,"error":"Internal Server Error","path":"/trans/system/dictionary/department/clinic/list"}]
最新发布
07-25
这个错误是由于在执行GET请求时，出现了一个内部服务器错误（Internal Server Error）。具体的错误信息是：{"timestamp":"2023-07-24T02:26:51.783+00:00","status":500,"error":"Internal Server Error","path":"/trans/system/dictionary/department/clinic/list"}。 这个错误可能是由服务器端的某些问题引起的，可能是由于代码逻辑错误、数据库连接问题或其他服务器端的异常情况导致的。如果你是系统的开发者，可以检查服务器端的日志或调试信息来查找具体的原因。如果你是系统的使用者，建议联系系统的开发团队或技术支持人员，向他们报告这个问题，并提供错误信息以帮助他们进行排查和修复。

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值