RNA-seq生信分析流程
RNA-seq生信分析流程
RNA-seq是近些年发展起来的针对转录组的测序技术,其能够获得mRNA、smallRNA以及各种非编码RNA的序列。
在不同细胞或者在相同细胞的不同发育阶段细胞中这些RNA的表达水平是不同的,依赖RNA-seq测序技术对这些不同细胞进行差异表达分析,可以得出转录组层面上的表达模式,当前已广泛应用于基础研究、临床诊断和药物研发等领域。
文章目录
- 1.数据质控
- 2.mapping
- 3.确定RNA丰度
- 4.差异基因表达分析
- 叮!
1.数据质控
由于当前的测序方法需要再测序样品的两端加上接头,以便于进行大量非特异性扩增以及高通量测序,因此测序结果首先要进行去接头和质控处理。数据质控即包括去接头和质控处理两部分,可以进行数据质控的工具有FastQC、cutadaptor、Trimmomatic、Trim galore等,其中较常见的是Trim galore。
2.mapping
进行数据质控后需要将获得的序列mapping回基因组,用于mapping的工具也有很多如tophat、Bowtie、hisat2、BWA等。几个工具正hisat2和Bowtie使用较多,Bowtie使用时较为节省内存并且超级快速,其并不只是一个简单的拼接工具,其适合的工作是将小序列mapping到大基因组上。hisat2速度也很快,也较为常用。
3.确定RNA丰度
确定各种RNA的丰度是非常重要的。不同样品间可以使用TMM值、upper quartile值来计算,同种样品则可以使用RPKM、quantile值进行比较。可以使用的工具有htseq-count或stringtie、Bedtools、GFold等。
4.差异基因表达分析
差异表达分析用于确定各RNA表达差异是否有显著性,差异分析的工具有CuffDiff、DEseq、EdgeR等,每种工具都有一定的特点如CuffDiff更适合于不同样本间的比较,而另外两种更适合于相同样本的比较。获取到差异分析结果后,还可以选择使用GO分析或者KEGG通路分析进行后续分析,用于确定这些差异表达代表的生物学意义。
叮!
当前除了RNA-seq外还有sRNA-seq、miRNA-seq、chip-seq、甲基化测序等测序技术,但是这些测序技术后续的生信处理的原理基本相似。本文仅非常简要地介绍了一下RNA-seq的大体步骤,很多细节信息未涉及,技术实现也未涉及到,因为现在还没做过这方面的课题,仅仅是记录一下之前课上学到的内容。以后若是有做到相关内容再更新吧!!!
RNA-seq生信分析流程相关推荐
- 生信分析流程构建的几大流派
导言 构建生信分析流程是生物信息学从业人员必备的技能之一,对该项能力的评估常常是各大公司招录人员的参考项目之一. 在进行 ngsjs 项目时,我做了一张示意图来表示一些高通量测序数据分析项目重现性的要 ...
- RNAseq生信分析流程简介
文章目录 RNAseq简介 RNAseq分析流程 1. 实验设计 1.1 目标RNA提取策略和文库类型选择 1.2 测序深度或文库大小 1.3 实验重复 2. RNAseq信息分析 2.1 质控 和 ...
- 生信分析流程框架概述
目录 1. 工作流管理系统 Scripts Make nextflow 2. 工作流引擎 Cromwell Argo 3. 工作流语言 Common Workflow Language (CWL) 1 ...
- 生物信息学分析服务器搭建教程,Snakemake搭建生信分析流程-步骤
什么是snakemake 官网地址:https://snakemake.readthedocs.io/en/stable/ 使用snakemake的好处 process data with the s ...
- 图形化开放式生信分析系统开发 - 9 Illumina测序仪测序数据自动拆分
前文链接: 图形化开放式生信分析系统开发 - 1 需求分析及技术实现 图形化开放式生信分析系统开发 - 2 样本信息处理 图形化开放式生信分析系统开发 - 3 生信分析流程的进化 图形化开放式生信分析 ...
- 图形化开放式生信分析系统开发- 1基本需求分析及技术实现
图形化开放式生信分析系统开发- 1基本需求分析及技术实现 起因/背景 软件获取:到官网sliverworkspace.com免费下载个人版,最新版本 2.0.277363 几张图片 下面进入正题,以具 ...
- 图形化开放式生信分析系统开发 - 7 分析报告的模板定制与自动生成
前文链接: 图形化开放式生信分析云平台产品开发 - 1 需求分析及技术实现 图形化开放式生信分析云平台产品开发 - 2 样本信息处理 图形化开放式生信分析云平台产品开发 - 3 生信分析流程的进化 图 ...
- 2020.09.30【RNA-seq流程】丨转录组生信分析全流程
RNA-Seq生信分析全流程 摘要 第一部分 step.1 下载数据 step.2 数据质控 第二部分 step.3序列比对 step.4 计算基因表达量 step.5 插入片段长度检验 step.6 ...
- 在B站学习大名鼎鼎的StatQuest 系列统计和生信分析视频(中文字幕)- 也见证助理教授到创业者的华丽转身...
生物信息学习的正确姿势 NGS系列文章包括NGS基础.在线绘图.转录组分析 (Nature重磅综述|关于RNA-seq你想知道的全在这).ChIP-seq分析 (ChIP-seq基本分析流程).单细胞 ...
最新文章
- 如何更优雅地对接第三方API
- window下python扩展包大全
- 目前有量子计算机,中国“祖冲之”刚刚成为当前最强大的量子计算机
- Spring AOP源码解析(一)——核心概念
- 手把手教你学Dapr - 4. 服务调用
- 火狐配置java_java selenium+firefox环境搭建
- 工作242:关于第二个git仓库提交代码
- 如何建立好的索引.--针对Distinct
- 安装Ubunutu音频视频库
- 【项目管理】项目管理计划
- ini文件中文乱码 python_python编程之ini文件处理-configparser模块应用
- Kaggle比赛(二)House Prices: Advanced Regression Techniques
- oracle期末试题及答案,oracle期末考试试题及答案
- 微型计算机相关的英文文献,微型计算机控制系统(单片机控制系统)外文文献翻译.doc...
- selenium自动化入门之实现163邮箱发送邮件
- 互联网巨头开放给创业者的新机会!
- Cluster - 基本概念
- 【语音智能平台】京东Alpha VS 讯飞开放平台 (一)
- 玩抖音是浪费时间吗?
- 数字信号处理之信号的抽取和内插
热门文章
- Oracle字符集讨论(转)
- 计算机等级考试二级 Python 语言程序设计考试大纲(2022 版)
- PCB电路板元件目标检测数据集(YOLO标签,1400张,9类目标)
- 从时间管理到管理情绪,这些自我管理的技巧你知道几个?
- word中怎样输入标准破折号“——”
- 如何为餐饮商家打造代运营解决方案?
- SpringBoot2.x(三)热部署devtool和配置文件自动注入实战
- 这是我见过描写天津女孩中最真实的
- 台式计算机32位和64位的区别,电脑操作系统32位和64位的区别
- html页面高度和宽度,html和body元素的高度和宽度