Gene prediction and comparison

使用基于基因组序列的从头预测方法(Ab initio methods同时分别使用头预测软件 GENSCAN AUGUSTUS预测exon和intron的剪切位点和基于证据支持的基因预测(Evidence-based methods与dog和human比较得到不同基因座数目。最后做了一个panda的参考基因集。

 

measure the quality of gene prediction):通过比较不同基因组(代表性动物)之间的key参数,得知panda与人是相似的。(evaluate the rate of missing exons in the predicted genes拿dog和panda与人比gene coverage(为什么dog也要比一下?因为dog与panda亲缘关系很近,所以以此作为可靠性依据),得知确实有oxen missing(因为1.末端2.因为小同时被长inron分离,所以难找);功能基因情况与dog相似。结论:预测效果好(因为可与好注释结果相比较)

 

identify gene orthologues among panda, dog, mouse and human):使用InParanoid and Multiparanoid被至少两个物种共享的同源基因,四个物种共享的同源基因,panda特有基因是狗特有基因(可能是因为dog genome多基于证据,而不是Ab initio methods预测)的两倍。

obtain greater insight into the evolutionary dynamics 进化动力学of the genesmaximum-likelihood比较每个物种(4个:dog,human,muse,panda)和祖先的区别(同源基因扩张或缩小),得知缺失的功能在进化上有很大作用(由图可知),同时发现参与受体活动的基因是显著差异基因。

maximum-likelihood:最大似然估计:将每个位置所有可能出现的残基替换概率(每个物种和祖先)进行累加,产生特定位点的似然值

looked for signatures of positive selection对这些基因做似然比(branch-site)测试,建立一个模型:input:1.人和panda的同源基因;2.panda和rat, mouse and dog中的一个此数据严格过滤),使用三种H1假设,以panda,Dog和5个物种至少一个作为假设主体,以确定panda,dog或五个物种中的一个是否受到正选择。

 

功能相同点:Fisher’s exact and Mann–Whitney U tests证明:Dog和 panda的共同的正向选择基因是与免疫和防卫有关基因与原先在哺乳动物中选择的正选择基因一致。GO分析得到panda的正选择基因在免疫系统中。

Fisher’s exact and Mann–Whitney U tests:H0非参数实验

GO分析:功能分析

功能不同点:从GO分析来说,熊猫(和其他6个哺乳动物基因相比)特有两个基因,在dog中两个基因被过表达。

(gain insight into some of the traits unique to the panda)

 

(diet):

吃竹子并不是熊猫基因决定(因为有编码消化系统基因,但未发现消化纤维的基因),而是消化系统微生物决定的。

介绍了编码五种味觉的基因,bitterness与dog有相似基因,Umami中T1R1是假基因(因为T1R1发生了移码错误),同时发现panda的该基因的ds/dn比dog低(因为panda该基因缺失)。因为该基因的缺失,使得panda虽然分类为食肉动物,但是实际上食草(该基因控制的一个受体可以感知食物中的某一个蛋白质,以此感知Umami,因为该基因缺失了,所以panda不具有该功能,所以无法感受到食物中的“鲜”)。

同义置换和非同义置换率(dS/dN ratio)的估计

fecundity):我们发现panda genome 中有与其他哺乳动物相似的性发育和性功能的基因。值得再讨论。

Ab initio methods|Evidence-based methods|maximum-likelihood|branch-site|H1|H0|GO|dS/dN ratio的更多相关文章

  1. 最大似然估计实例 | Fitting a Model by Maximum Likelihood (MLE)

    参考:Fitting a Model by Maximum Likelihood 最大似然估计是用于估计模型参数的,首先我们必须选定一个模型,然后比对有给定的数据集,然后构建一个联合概率函数,因为给定 ...

  2. Maximum Likelihood及Maximum Likelihood Estimation

    1.What is Maximum Likelihood? 极大似然是一种找到最可能解释一组观测数据的函数的方法. Maximum Likelihood is a way to find the mo ...

  3. 似然函数 | 最大似然估计 | likelihood | maximum likelihood estimation | R代码

    学贝叶斯方法时绕不过去的一个问题,现在系统地总结一下. 之前过于纠结字眼,似然和概率到底有什么区别?以及这一个奇妙的对等关系(其实连续才是f,离散就是p). 似然函数 | 似然值 wiki:在数理统计 ...

  4. Don’t Use Accessor Methods in Initializer Methods and dealloc 【初始化和dealloc方法中不要调用属性的存取方法,而要直接调用 _实例变量】

    1.问题:    在dealloc方法中使用[self.xxx release]和[xxx release]的区别? 用Xcode的Analyze分析我的Project,会列出一堆如下的提示:Inco ...

  5. 最大似然估计(Maximum Likelihood,ML)

    先不要想其他的,首先要在大脑里形成概念! 最大似然估计是什么意思?呵呵,完全不懂字面意思,似然是个啥啊?其实似然是likelihood的文言翻译,就是可能性的意思,所以Maximum Likeliho ...

  6. Linear Regression and Maximum Likelihood Estimation

    Imagination is an outcome of what you learned. If you can imagine the world, that means you have lea ...

  7. Maximum Likelihood Method最大似然法

    最大似然法,英文名称是Maximum Likelihood Method,在统计中应用很广.这个方法的思想最早由高斯提出来,后来由菲舍加以推广并命名. 最大似然法是要解决这样一个问题:给定一组数据和一 ...

  8. Maximum likelihood from incomplete data via the EM algorithm (1977)

    Maximum likelihood from incomplete data via the EM algorithm (1977)  

  9. [Bayes] Maximum Likelihood estimates for text classification

    Naïve Bayes Classifier. We will use, specifically, the Bernoulli-Dirichlet model for text classifica ...

随机推荐

  1. 点击实现CSS样式切换

    如图所示 代码如下图: 特别要注意的是:a标签不会继承上级的color,所以要单独为其设置 参看代码(并非上图代码)如下: <!DOCTYPE html> <html> < ...

  2. 通过node.js读取html页面及其页面中引入的css样式

    Node 是什么? Node 是一个基于Chrome V8 引擎的JavaScript 运行环境. Node 不是一种独立的语言.Node不是 JavaScript 框架, Node是一个除了浏览器之 ...

  3. java基础笔记(三)——main方法

    1.解析public static void main(String[] args)方法 JVM在运行程序时,会首先查找main()方法作为入口,main是JVM识别的特殊方法名. public是权限 ...

  4. 使用CSS3的translate和transition功能,控制一个两个div块的联动

    之前的工作有接触到一些css3的新特性.div块的移动和回到初始位置,可以利用在开发中的很多地方.这里记录下来,下次就不用辛苦的灾区百度了. <html> <head> < ...

  5. Nginx托管.Net Core应用程序

    Nginx托管.Net Core应用程序 一.安装.Net Core 参考官方文档:https://www.microsoft.com/net/core#linuxcentos 1.添加dotnet产 ...

  6. 061 Rotate List 旋转链表

    给定一个链表,将链表向右旋转 k 个位置,其中 k 是非负数.示例:给定 1->2->3->4->5->NULL 且 k = 2,返回 4->5->1-> ...

  7. 大都市 meg

    Description 在经济全球化浪潮的影响下,习惯于漫步在清晨的乡间小路的邮递员Blue Mary也开始骑着摩托车传递邮件了.不过,她经常回忆起以前在乡间漫步的情景. 昔日,乡下有依次编号为1.. ...

  8. iOS 面试常问之多线程

    本片围绕多线程全面展开叙述. 1.为什么要有多线程/多线程是用来干什么的? 2.多线程是什么? 3.如何创建多线程? 4.多线程在哪些情况下会使用/多线程使用场景? 5.三种多线程的优缺点? 6.线程 ...

  9. JAVA基础之转换流和缓冲流

    个人理解: 在理解的字符流和字节流的区别后.要是想读取指定的编码格式的文件时,特别是不是默认的格式时,就需要转换流了,需要注意的是字符流是需要清除缓冲区的:当需要快速的进行读取时,则需要缓冲流.存在即 ...

  10. Spring的七种事务传播机制

    概述 当我们调用一个基于Spring的Service接口方法(如UserService#addUser())时,它将运行于Spring管理的事务环境中,Service接口方法可能会在内部调用其它的Se ...