2020.11.3【读书笔记】丨生物信息学与功能基因组学(第二章 序列数据的获取 下)
- 2.7 基因组浏览器
- 三个主流基因组浏览器:Ensemble、UCSC和NCBI
- 基因组组装
- 定义:
- 所获得的一个物种DNA序列按照染色体的形式进行的一种组装。
- 内容:
- 对基因组的注释,如起止位点、外显子、DNA重复元件或其他基因组特征。
- 基因组参照联盟(GRC)
- 维护人类、小鼠和斑马鱼的参考基因组。
- 对基因组组装需要明确几点:
- 每个染色体的起止位置是什么
- 基因组序列中包括多少gap区域,能否被弥补?
- 基因组上的结构变异位点如何被展现?可否被发现并更正?
- 一个基因组组装中有多少碱基是错误的?
- UCSC Genome Browser
- UCSC浏览器目前支持36个脊椎和非脊椎动物基因组的分析,是目前针对人类和其他重要物种应用最广泛的基因组浏览器。UCSC Genome Browser 提供了不同分辨率的染色体位置图形视图。
- 图示说明如何使用UCSC
- Ensembl Genome Browser
- 提供各种真核生物的一系列综合网站。
- 该网站目的是自动分析和注释基因组数据并能以浏览器的形式来战士基因组数据。
- Ensemble固定标识符
- NCBI Map Viewer
- Map Viewer包括了后生动物、真菌和植物等多个物种的染色体图谱。
- Map Viewer允许基于文本和序列的搜索
- 提供基因组四个层次的详细信息
- 物种主页
- 基因组视图
- 图谱视图
- 序列视图
- 定义:
- 2.8 获取序列数据的案例:单个基因、蛋白质
- 以人类组蛋白和HIV-1pol蛋白为例
- 组蛋白
- NCBI Protein收录了470000个组蛋白条目
- 获取步骤:
- 1.在NCBI Protein确定Homo sapiens物种
- 2.检索“txid9606[Organism;exp]histone”,可以得到8000个人类组蛋白,超过2000个有RefSeq索引编号;去掉组蛋白去乙酰没和组蛋白乙酰转移酶可以得到含有1700多个RefSeq索引编号的蛋白质
- 3.通过RefSeq了解组蛋白概述
- 4.随意挑选一个组蛋白进行研究
- 5.在组蛋白质序列数据库(Marinno-Ramiirez等,2011)中可发现人类人类基因组中有113个组蛋白基因,包含一个基因簇
- 6.通过其他蛋白质家族数据库(Pfam和InterPro)可以提供关于蛋白质和基因家族的简介描述以及其中代表性的蛋白质和基因。
- HIV-1 pol
- HIV-1的RNA依赖性DNA聚合酶是一个反转录酶,其编码基因成为pol(polymerase)。如何获取反转录酶的DNA和蛋白质序列?
- 获取步骤:
- 1.在NCBI Nucleotide输入hiv-1,得到3000+个RefSeq
- 2.确定物种后,得到一个RefSeq条目(NC_001802.1)
- 注意:NCBI其他数据库(UniGene、OMIM等)并不适合查找病毒反转录酶的序列
- 访问数据集:区域和特征的大规模查询
- 把一个基因和多个基因进行对比
- 提问:
- 完整的人类球蛋白基因都包括那些基因?
- 这些基因都位于哪条染色体上?
- 11号染色体上有多少外显子?每个外线组内部有多少重复元件?
- 重点使用资源:Ensemble(BioMart项目)和UCSC(Table Browser)
- 提问:
- BioMart项目
- 提供了对多个数据库中的海量数据的便捷访问
- 该数据库基于以下两个原则建立:
- 数据不可知性的建模
- 从第三方数据库引入数据集合后后使用一种关联模式来访问这些数据
- 数据联合的形式
- 将许多不同的数据库联合成一个单一、完整、虚拟的数据库。
- 数据不可知性的建模
- UCSC Table Browser
- UCSC Table Browser可以在UCSC Genome Browser可视化的同事可以以精确和完整的表格形式来进行呈现。这些表格可被浏览,查阅和下载
- 以人类基因组为例
- 以人类基因组为例
- UCSC Table Browser可以在UCSC Genome Browser可视化的同事可以以精确和完整的表格形式来进行呈现。这些表格可被浏览,查阅和下载
- 自定义数据轨:BED文件的多用途性
- 对浏览器信息进行自定义的原因
- 要获取给定类别的信息
- 要上传我们感兴趣的信息
- 可使用BEDtools软件来分析BED文件
- BED文件格式必含内容
- 染色体号、染色体起始位置和染色体终止位置
- BED文件格式可选内容
- 第四列:名字
- 第五列:分数/正负链
- 第七、八列:加粗起始/终止位置
- 第九列:RGB颜色值
- 第十至十二列:区块(外显子)数目、大小、起始。
- 可支持的自定义文件格式
- 表格
- 注意事项
- 染色体可否使用数字表示
- 位置计数是0-based还是1-based
- 不同数据来源的计数表示
- 不同数据来源的计数表示
- 表格
- 对浏览器信息进行自定义的原因
- Galaxy:可重复的、基于网页的高通量研究
- 简介
- Galaxu是基于网页的分析平台,可接受包括如BioMart和UCSC Table Browser 等多个来源的输入文件。
- 优势
- 提供大量整合的工具可用于多种类型数据的输入和分析
- 基于网页并提供了大量在其他平台上需要通过命令行来执行的软件包
- 每一个分析步骤会被记录、储存并与他人共享,促进可重复性研究。
- 使用场景
- 获取蛋白质序列并进行序列比对
- 基因组DNA序列比对
- 在对染色体进行探索时,可获取人类基因组是的微卫星序列;我们可以创制表格,对未知排序找出最长的一个。
- 分析二代测序数据时,可以输入FASTQ文件,执行序列比对,并分析BAM和VCF文件
- 由于含有RNA-seq成套软件被大家经常使用,比对工具Bowtie和BWA都可以在Galaxy上使用
- 简介
- 把一个基因和多个基因进行对比
- 2.9 生物医学文献获取
- 美国国家医学图书馆(NLM)创立数目数据库-MEDLINE(医学文献分析和联机检索系统)可通过NCBI旗下,PubMed进行免费访问。
- PubMed搜索
- 使用过滤器可以显示可免费获取的文章
- 截断单词使用通配符(例:therap*)可以获得单词的不同变体。
- PubMed搜索
- 美国国家医学图书馆(NLM)创立数目数据库-MEDLINE(医学文献分析和联机检索系统)可通过NCBI旗下,PubMed进行免费访问。
- 2.10 展望
- 三大数据库:GenBank、EMBL-Bank、DDBJ
- 没有一个唯一正确的方式来获取信息,有许多方法都是可行的
- 数据库资源都是紧密相关并相互提供访问链接
- 2.11 常见问题
- 明确目标(基因名称、需要DNA还是蛋白质序列)
- 添加适当的限制条件进行特定搜索
- 人类基因组组织(HUGO)和基因命名委员会(HGNC)可以展示人类基因的官方基因符合并附有数据资源的相关链接。
- 2.12 给学生的建议
- 深入了解每个生物信息学网站(EBI、NCBI、Ensemble、UCSC)
- 尝试比对一次研究一个基因和一组基因
2020.11.3【读书笔记】丨生物信息学与功能基因组学(第二章 序列数据的获取 下)相关推荐
- 【读书笔记】《Effective Java》第二章 第2条:遇到多个构造器参数时要考虑使用Builder
一.前言 <Effective Java>读书笔记系列 第二章 第1条:创建和销毁对象 第二章 第2条:遇到多个构造器参数时要考虑使用Builder 二.介绍 我们开发中偶尔会遇到一些需要 ...
- 2021.11.22【读书笔记】丨生物信息学与功能基因组学(第五章第五节 用类似于BLAST的比对工具快速搜索基因组DNA)
5.5 用类似于BLAST的比对工具快速搜索基因组DNA 需求:随着基因组DNA数据库数量增长,对比对工具要求越来越高 能在基因组DNA中找到外显子 比对时考虑基因组DNA包含的测序错误 有相应的算法 ...
- 2021.11.21【读书笔记】丨生物信息学与功能基因组学(第五章 高级数据库搜索 中 )
5.3 寻找远缘相关蛋白质:位置特异性迭代BLAST(PSI-BLAST)和DELTA-BLAST PAM250矩阵给探测远缘相关蛋白质提供了一个更好的打分系统,可以改变打分矩阵来检测远缘蛋白质,但仍 ...
- 2021.12.19【读书笔记】丨生物信息学与功能基因组学(第五章 高级数据库搜索 下)
5.5 用类似于BLAST的比对工具快速搜索基因组DNA 需求:随着基因组DNA数据库数量增长,对比对工具要求越来越高 能在基因组DNA中找到外显子 比对时考虑基因组DNA包含的测序错误 有相应的算法 ...
- 2022.01.23【读书笔记】丨生物信息学与功能基因组学(第六章 多重序列比对 上)
学习目标 理解使用ClustalW进行多重序列比对(MSA)的三个主要阶段: 描述几种其他的多重序列比对(MSA)程序,了解他们的工作原理,比对它们与ClustalW的异同: 理解进行基准研究的重要性 ...
- 【读书笔记】Linux内核完全注释第二章:微型计算机组成结构
非常主观的笔记,只是摘抄一下自己觉得比较重要的地方,并不会概述所有内容(概括能力太差是这样的) 2.微型计算机组成结构 2.1 微型计算机组成原理 除了 CPU 以外,现代 PC 机主板主要使用 2 ...
- 【读书笔记】周志华 机器学习 第二章 模型评估和选择
第二章 模型评估和选择 1 欠拟合和过拟合 偏差和方差 1.1 欠拟合和过拟合 1.2 偏差和方差 2 评估方法 性能度量 2.1 评估方法 2.1.1 留出法 2.2.2 交叉验证法 2.2.3 自 ...
- CSAPP读书笔记与习题作业练习-第3章
CSAPP读书笔记与习题作业练习-第3章 疑问 疑问一 练习题3.26中为什么通过循环+异或可以判断二进制数中1的个数的奇偶性(已解决) 习题 练习题3.1 练习题3.2 练习题3.3 练习题3.4 ...
- 读书笔记|《金字塔原理》_第一章
读书笔记|<金字塔原理>_第一章 [章节]--第一章为什么要用金字塔结构 [讲了什么] 1.首先是为什么要用金字塔结构? 2.如何将思想组织成金字塔结构? 2.1 归类分组,将思想组织成金 ...
最新文章
- CUBLAS_STATUS_ALLOC_FAILED CUDA_ERROR_OUT_OF_MEMORY
- Python 输入一些数,统计最大值及其出现的频率,求一个数的全部质因数
- Java类加载器( 死磕9)
- bzoj 3190 赛车 半平面交
- ERROR 1044 (42000): Access denied for user ''@'...
- 神经网络入门(详细 )
- python小工具—图片拼接PDF
- Linux学习总结(1)——Linux命令大全完整版
- DXP2004生成PCB不显示连线
- 还在纠结程序员能不能干一辈子吗?有答案了!!
- Java Web程序设计第12章课后第2题
- Android 应用程序签名、发布
- a标签去掉下划线_HTML 所有标签列表
- 建立视图实发工资mysql_MYSQL人事工资管理系统--设计视图(四)
- MKD3用法之一构造虚假WIFI热点
- 趣拍云短视频SDK全面升级,简单易用引开发者点赞
- 英语口语259之每日十句口语
- 多对多关系的多表关联查询
- 海康、大华、星邦网络摄像头的 RTSP协议 地址与格式
- PSD文件模板更改照片和文字
热门文章
- 首批 5G 网络就位之后,下一步在何方?
- 02 FizzBuzz小游戏
- 转 PCB各层的含义 (solder paste 区别)
- 微软ASP.NET 电商网站开发实战 MVC6 +HTML5 +WCF+WebAPI+NoSQL+mongoDB+Redis+Core视频 代码 面试题...
- 安装蓝点Linux2.0 要小心(转)
- [ 物联网篇 ] 41 - 嵌入式系统中的epoll_wait 的异常处理
- 汇佳学校体育|昌熙桓:84场60次夺冠,高尔夫是我终身的热爱
- java 泛化_(转)Java中的泛化,依赖,关联,聚合,合成
- caj文献放进zotero
- pharo-smalltalk语言