RNA-Seq [Inkit wiki]

用户工具

站点工具


stat:rnaseq

RNA-Seq

教程

流程

https://bioconductor.org/packages/release/BiocViews.html#___GeneExpressionWorkflow
Raw Data: Fastq (.gz) →

  1. QC: FastQC
  2. Filter
  3. Alignment: STAR
  4. Count
  5. Cluster
  6. Heatmap
  7. Differential

1 QC

FastQC

apt install fastqc

fastqc --noextract RawData/I409/I409_1.fq.gz -o results/1_initial_qc/ 

2 Alignment

基因组注释数据

Alignment Indexing

对齐工具

HISAT2

https://daehwankimlab.github.io/hisat2/download/ * https://notebook.community/ssjunnebo/pathogen-informatics-training/Notebooks/RNA-Seq/genome-mapping

hisat2-build -p 32 fasta/Rattus_norvegicus.mRatBN7.2.dna.toplevel.fa hisat_index
hisat2 -x hisat_index/hisat_index -1 I409/I409_1.fq.gz -2 I409/I409_2.fq.gz -S I409.sam -p 32

STAR

2.7.10b

apt install rna-star

# 创建索引,索引文件创建一次即可. 需要从Ensembl下载对应物种的Fasta文件和GTF文件。
STAR --runMode genomeGenerate --genomeDir star_index --genomeFastaFiles fasta/* --sjdbGTFfile gtf/* --runThreadN 14

# 运行分析
STAR --genomeDir star_index --readFilesIn filtered/sample_filtered.fq --outSAMtype BAM SortedByCoordinate --quantMode GeneCounts --runThreadN 14

STAR --genomeDir star_index --readFilesIn rna4/RawData/I409/I409_1.fq --outSAMtype BAM SortedByCoordinate --quantMode GeneCounts --runThreadN 14

kallisto

# bat
kallisto bus [arguments] FASTQ-files

kallisto quant -i rattus_index_ki/transcriptome.idx -o reads.kallisto_quant -t 64 --fusion --pseudobam --genomebam --gtf gtf\Rattus_norvegicus.mRatBN7.2.109.gtf rna4\RawData\I409\I409_1.fq.gz rna4\RawData\I409\I409_2.fq.gz

3 Count

featureCounts -p -M -O -T 32 -a gtf/Rattus_norvegicus.mRatBN7.2.109.gtf -o output.txt data.sam [bam]
stat/rnaseq.txt · 最后更改: 2024/04/04 23:32 由 inkit