写在前面

写这个贴子纯粹为了记录一下自己的双面打工人的生活,然后再总结一下自己的工作方法还有哪些地方存在问题,希望下次再遇到这些小问题的时候可以更快速地应对。

同源关系如何定义| How to define ortholog?

可以参考这篇文章中的记录:https://www.jianshu.com/p/82d4cf6c3eda

简单来说,就是speciation events得到的是一对或者一组基因,称为orthologs,而gene duplication事件形成的一对或者一组基因,称为paralogs。

而在实际的分析中,往往得到的是一组orthogroup(比如OrthoFinder 2.0的结果文件),其中不仅包含了ortholog还包含了paralog。

流程实现

注意

分析过程中需要注意的问题,

  • 从蛋白质序列的相似度、系统发育的角度,一个species内的paralog一般比两个species之间的ortholog要呈现出更高的相似度

    举个例子,比如水稻和拟南芥的某一个disease-resistant gene,Arab genome内部的paralog之间的序列相似度,相较于rice和Arab之间的相似度,前者的相似度更高,反映在最终的phylogeny上,也是Arab的paralog之间会更倾向于聚集在一个clade上,而不是单独与rice genome中的ortholog聚集在一个clade上

  • 先鉴定出一个species内部,一个gene的orthogroup,可能能够增加鉴定到另一个species中ortholog的概率

思路

分析的一个流程如下,

  • OrthoFinder 2.0
  • seqkit提取targeted gene
  • reciprocal BLAST
  • JCVI提取每个query的最优比对结果
  • MAFFT
  • IQtree 2.0

同源基因鉴定 | OrthoFinder 2.0 + MAFFT + IQtree相关推荐

  1. 测序结果峰图分析物种鉴定方法1.0

    # ITS ATCATTA| → → →← ← ← GACCTCA|AATCA# LSU TAACAAGG|TCCCC → → → ← ← ← GAAGAAGCC| 笨办法,仅供参考,后续会有更新!! ...

  2. 利用 OrthoFinder、IQtree、Notung、iTOL 绘制基因树

    本文为以 4 个兰科物种中 NB-ARC 结构域的序列信息作为输入,构建 有根基因树.步骤参考文献 Genome- Wide Analysis of the Nucleotide Binding Si ...

  3. 使用OrthoFinder进行基因家族分析

    谈论到直系同源基因分析的时候,大部分教程都是介绍OrthoMCL,这是2003年发表的一个工具,目前的引用次数已经达到了3000多,但这个软件似乎在2013年之后就不在更新,而且安装时还需要用到MyS ...

  4. iqtree+ggtree绘制进化树教程

    写在前面 最近,自己在绘制进化树的分析,那么自己也总结一下相关的教程吧.对于进化树绘制的教程全网依旧是很多的.自己作为初学者,也是根据这些教程来学习,来做分析.做物种或基因的进化树,软件很多,目前很多 ...

  5. 纯python统计基于orthofinder得到的系统发育树的关注点位置的树型结构

    对于某一个物种或某类物种在整个系统发育树中的位置存在一定争议的情况,使用直系同源基因构建单基因树,并对该物种或该类物种所在结构进行统计是可以对争议起到一定的解决作用的,在此留下全套流程和大家交流. 主 ...

  6. 估算带卷积核二分类0,3的网络的收敛时间和迭代次数

    制作一个网络分类minst的0和3求出这网络的迭代次数曲线表达式n(δ),和准确率表达式p-max(δ),用预期准确率去估算n,并推算需要的时间. 将minst的28*28的图片缩小到9*9,网络用一 ...

  7. 遗传:微生物组数据分析方法与应用

    本文版权归<遗传>杂志,已获授权,转载请联系杂志社 微生物组数据分析方法与应用 刘永鑫1,2,秦媛1,2,3,郭晓璇1,2,白洋1,2,3 1. 中国科学院遗传与发育生物学研究所,植物基因 ...

  8. 基于TBtools做基因家族分析教程 (全)

    基因家族分析笔记-全部开始记录 一. 写在前面 2023年4月中旬自己开始做基因家族的分析,对于这块自己没有接触过,因此也是一个挑战,没事!!!(安慰自己),对于基因家族的分析网上的教程很多,跟着步骤 ...

  9. MPB:上海巴斯德所崔杰组-RNA病毒组与生物信息学分析

    为进一步提高<微生物组实验手册>稿件质量,本项目新增大众评审环节.文章在通过同行评审后,采用公众号推送方式分享全文,任何人均可在线提交修改意见.公众号格式显示略有问题,建议电脑端点击文末阅 ...

最新文章

  1. android 实现全屏代码
  2. tomcat 随windows启动
  3. MATLAB求解常微分方程
  4. 1.17 StringBuffer类详解
  5. linux 内核升级 网络 不能上网,Deepin Linux 无法上网
  6. ELK技术栈—Logstash—基础介绍
  7. mysql中utf8_general_cs和utf8_general_ci和utf8_bin的一点区
  8. c库函数tmpfile
  9. keepalive之LVS-DR架构
  10. Invalid options object. Copy Plugin has been initialized using an options object that does not match
  11. 作者:​冯景华(1984-),男,国家超级计算天津中心主任助理、系统管理部部长。...
  12. 角点检测--基于梯度的方法(Moravec角点检测、Harris角点检测、Shi-Tomasi角点检测)
  13. 工业机器人实训耗材_工业机器人实训室
  14. LeetCode题解-6. ZigZag Conversion
  15. Java:打包成jar包
  16. Unity实战篇 | unity接入QQ登录 详细过程——Android篇
  17. web 前端后端分工
  18. 双目标定(三)标定流程(含矫正)
  19. 【UVM基础】UVM 的 build_phase 执行顺序
  20. uc浏览器手机端禁用复制

热门文章

  1. NIOS_II填坑之路——EPCS出现“Cannot open flash device”解决办法
  2. CSS-子元素选择器
  3. 在下列用户中什么管理计算机的权限最小,计算机管理员考试试题和答案
  4. 【整理】HLS视频协议第二弹--裁剪部分视频及m3u8文件,编写通用客户端以播放m3u8视频
  5. 在c++中使用easyx画一个实时走动的钟表(方法细节)
  6. SLAM技术学习笔记
  7. mac开发 oc 授权检测、申请、打开授权设置页面
  8. 宕机超 12 小时,损失过亿,唯品会基础平台负责人被“祭天”
  9. Web前端学习之正向代理与反向代理
  10. java如何输入2的31次方_续一: 如何优化Java程序:十进制转十六进制(2的31次方以内的正整数)...