RNA-seq分析(Fastqc+Trimmomatic+STAR+HTseq-count+DESeq2)

2024-03-28 21:18

本文主要是介绍RNA-seq分析(Fastqc+Trimmomatic+STAR+HTseq-count+DESeq2),希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

最近做RNA-seq,正好把流程整理下,也希望分享和相互学习。
具体将以Fastqc + Trimmomatic + STAR + HTseq-count + DEseq2的流程来进行。

查看数据完整性

for dir in `ls`; do cd $dir; md5sum -c MD5_*txt; cd ..; done

预处理

FastQC + Trimmomatic

fastqc -t 5 sample_R1.fq.gz
fastqc -t 5 sample_R2.fq.gz
java -jar ~/tools/Trimmomatic/Trimmomatic-0.36/trimmomatic-0.36.jar PE -threads 20 sample_R1.fq.gz sample_R2.fq.gz -baseout sample_filtered.fq.gz ILLUMINACLIP:~/tools/Trimmomatic/Trimmomatic-0.36/adapters/TruSeq3-PE.fa:2:30:10 LEADING:3 TRAILING:3 SLIDINGWINDOW:4:15 HEADCROP:8 MINLEN:36 HEADCROP:15

fastqc后发现有些样本per tile sequence content 1,Per base sequence content,Adapter Content,Kmer Content没有通过。主要问题是去除些质量差的reads;根据前15个左右碱基比不均一,用HEADCROP去掉。用的是TruSeq的adapter,故而加上,用Trimmomatic。
Trimmomatic相关学习内容,见2,3.4.

STAR

make index
人和小鼠的基因组和参考注释用Tophat的igenomes下:

STAR --runThreadN 30 --runMode genomeGenerate --genomeDir STARINDEX_20180118/ --genomeFastaFiles WholeGenomeFasta/genome.fa --sjdbGTFfile ../Annotation/Genes/genes.gtf --sjdbOverhang 134

do the alignment.
可以基于第一次比对的结果,用SJ.out.tab于重新Genome的Index,然后再比对(在用找SNP和Indel时尤其推荐)。7

STAR --runThreadN 30 --genomeDir ~/Ref/UCSC_hg19/Homo_sapiens/UCSC/hg19/Sequence/STARIndex_20180118 --readFilesIn sample_filtered_1P.fq.gz sample_filtered_2P.fq.gz --outFileNamePrefix ./Hs_treat3/Hs_treat3 --readFilesCommand zcat

参考内容:5, 6,
Trim reads map to multiple regions.

samtools view -bS -F 4 Hs_treat3Aligned.out.sam > Hs_treat3_mapped.bam
samtools sort -n Hs_treat3_mapped.bam Hs_treat3_sort

HTSeq

用htseq-count计算read counts。8,9

htseq-count -f bam -s no Hs_treat3_sort.bam ~/Ref/UCSC_hg19/Homo_sapiens/UCSC/hg19/Annotation/Genes/genes.gtf > sample.count

DESeq2差异分析

library(DESeq2)
condition <- factor(c("A","A","B","B"))
dds <- DESeqDataSetFromMatrix(hs, DataFrame(condition), ~ condition)
dds <- dds[ rowSums(counts(dds)) > 1, ]   #过滤low count数据
nrow(dds)
dds <- DESeq(dds)     #差异分析
res <- results(dds)   #用result()函数获取结果
summary(res)  #summary()函数统计结果
count_r <- counts(dds, normalized=T)  #提取normalized count matrix

10

这篇关于RNA-seq分析(Fastqc+Trimmomatic+STAR+HTseq-count+DESeq2)的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/856718

相关文章

Android kotlin中 Channel 和 Flow 的区别和选择使用场景分析

《Androidkotlin中Channel和Flow的区别和选择使用场景分析》Kotlin协程中,Flow是冷数据流,按需触发,适合响应式数据处理;Channel是热数据流,持续发送,支持... 目录一、基本概念界定FlowChannel二、核心特性对比数据生产触发条件生产与消费的关系背压处理机制生命周期

怎样通过分析GC日志来定位Java进程的内存问题

《怎样通过分析GC日志来定位Java进程的内存问题》:本文主要介绍怎样通过分析GC日志来定位Java进程的内存问题,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教... 目录一、GC 日志基础配置1. 启用详细 GC 日志2. 不同收集器的日志格式二、关键指标与分析维度1.

MySQL count()聚合函数详解

《MySQLcount()聚合函数详解》MySQL中的COUNT()函数,它是SQL中最常用的聚合函数之一,用于计算表中符合特定条件的行数,本文给大家介绍MySQLcount()聚合函数,感兴趣的朋... 目录核心功能语法形式重要特性与行为如何选择使用哪种形式?总结深入剖析一下 mysql 中的 COUNT

MySQL中的表连接原理分析

《MySQL中的表连接原理分析》:本文主要介绍MySQL中的表连接原理分析,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教... 目录1、背景2、环境3、表连接原理【1】驱动表和被驱动表【2】内连接【3】外连接【4编程】嵌套循环连接【5】join buffer4、总结1、背景

python中Hash使用场景分析

《python中Hash使用场景分析》Python的hash()函数用于获取对象哈希值,常用于字典和集合,不可变类型可哈希,可变类型不可,常见算法包括除法、乘法、平方取中和随机数哈希,各有优缺点,需根... 目录python中的 Hash除法哈希算法乘法哈希算法平方取中法随机数哈希算法小结在Python中,

Java Stream的distinct去重原理分析

《JavaStream的distinct去重原理分析》Javastream中的distinct方法用于去除流中的重复元素,它返回一个包含过滤后唯一元素的新流,该方法会根据元素的hashcode和eq... 目录一、distinct 的基础用法与核心特性二、distinct 的底层实现原理1. 顺序流中的去重

关于MyISAM和InnoDB对比分析

《关于MyISAM和InnoDB对比分析》:本文主要介绍关于MyISAM和InnoDB对比分析,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教... 目录开篇:从交通规则看存储引擎选择理解存储引擎的基本概念技术原理对比1. 事务支持:ACID的守护者2. 锁机制:并发控制的艺

MyBatis Plus 中 update_time 字段自动填充失效的原因分析及解决方案(最新整理)

《MyBatisPlus中update_time字段自动填充失效的原因分析及解决方案(最新整理)》在使用MyBatisPlus时,通常我们会在数据库表中设置create_time和update... 目录前言一、问题现象二、原因分析三、总结:常见原因与解决方法对照表四、推荐写法前言在使用 MyBATis

Python主动抛出异常的各种用法和场景分析

《Python主动抛出异常的各种用法和场景分析》在Python中,我们不仅可以捕获和处理异常,还可以主动抛出异常,也就是以类的方式自定义错误的类型和提示信息,这在编程中非常有用,下面我将详细解释主动抛... 目录一、为什么要主动抛出异常?二、基本语法:raise关键字基本示例三、raise的多种用法1. 抛

github打不开的问题分析及解决

《github打不开的问题分析及解决》:本文主要介绍github打不开的问题分析及解决,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教... 目录一、找到github.com域名解析的ip地址二、找到github.global.ssl.fastly.net网址解析的ip地址三