DOC_ID : T11-0001
Doc_ID: A08-0001GRCm39en104Star275a
Editor: Mira
Reviewer:
Description
RSEM is a software package for estimating gene and isoform expression levels from RNA-Seq data. The RSEM package provides an user-friendly interface, supports threads for parallel computation of the EM algorithm, single-end and paired-end read data, quality scores, variable-length reads and RSPD estimation. In addition, it provides posterior mean and 95% credibility interval estimates for expression levels.
Build RSEM references using RefSeq, Ensembl, or GENCODE annotations
RefSeq and Ensembl are two frequently used annotations. For human and mouse, GENCODE annotaions are also available. Here, we show how to build RSEM references using Ensembl annotation. It is important to use every genome version with it’s compatible gtf file.
Source
- URL :
- DNA Sequence :http://ftp.ensembl.org/pub/release-104/fasta/mus_musculus/dna/=>Mus_musculus.GRCm39.dna.primary_assembly.fa.gz
- GTF : http://ftp.ensembl.org/pub/release-104/gtf/mus_musculus/=>Mus_musculus.GRCm39.104.gtf.gz
- File size :
- Mus_musculus.GRCm39.dna.primary_assembly.fa.gz 806.4MB (解壓縮後2.71GB)
- Mus_musculus.GRCm39.104.gtf.gz 31.3MB (解壓縮後0.83GB)
- Mus_musculus.GRCm39.dna.primary_assembly.fa.gz 806.4MB (解壓縮後2.71GB)
Genome assemble version : GRCm39 Release 104
Detail information :
使用STAR version 2.7.5a及RSEM version 1.3.3 來製作index,請確認已完成GApp standard analysis environment的安裝
#Activate standard analysis environment
conda activate GApp
#移動到Ref資料夾
cd ~/GA_bundle/Ref/
#創建放置fasta及GTF的資料夾
mkdir -p Mus_musculus/GRCm39en104Star2.7.5a/
#移動至資料夾
cd Mus_musculus/GRCm39en104Star2.7.5a/
#下載fasta
wget ftp://ftp.ensembl.org/pub/release-104/fasta/mus_musculus/dna/Mus_musculus.GRCm39.dna.primary_assembly.fa.gz
#解壓縮檔案
gunzip Mus_musculus.GRCm39.dna.primary_assembly.fa.gz
#下載GTF
wget ftp://ftp.ensembl.org/pub/release-104/gtf/mus_musculus/Mus_musculus.GRCm39.104.gtf.gz
#解壓縮檔案
gunzip Mus_musculus.GRCm39.104.gtf.gz
#用RSEM及STAR prepare index
rsem-prepare-reference \
–gtf ~/GA_bundle/Ref/Mus_musculus/GRCm39en104Star2.7.5a/Mus_musculus.GRCm39.104.gtf \
–star \
-p 20 \
~/GA_bundle/Ref/Mus_musculus/GRCm39en104Star2.7.5a/Mus_musculus.GRCm39.dna.primary_assembly.fa \
~/GA_bundle/Ref/Mus_musculus/GRCm39en104Star2.7.5a/GRCm39.104.genome
使用gtfToGenePred及samtools version 1.10製作用來提供基因位置及rRNA位置的訊息的ref_flat file及 ribosomal intervals file
#利用gtfToGenePred將gtf轉為ref_flat檔
gtfToGenePred -genePredExt -geneNameAsName2 -ignoreGroupsWithoutExons Mus_musculus.GRCm39.104.gtf /dev/stdout |awk ‘BEGIN { OFS=”\t”} {print $12, $1, $2, $3, $4, $5, $6, $7, $8, $9, $10}’ > Mus_musculus.GRCm39.104.gtf.refflat
#由reference genome抽取genome大小資訊
#Step1
samtools faidx Mus_musculus.GRCm39.dna.primary_assembly.fa
#Step2
cut -f1,2 Mus_musculus.GRCm39.dna.primary_assembly.fa.fai > sizes.genome
#Step3
perl -lane ‘print “\@SQ\tSN:$F[0]\tLN:$F[1]\tAS:GRCm39″‘ sizes.genome |grep -v _ >> Mus_musculus.GRCm39.104.gtf.rRNA.refflat
#合併由GTF抽出的rRNA資訊
grep ‘gene_biotype “rRNA”‘ Mus_musculus.GRCm39.104.gtf |awk ‘$3 == “gene”‘ |cut -f1,4,5,7,9 |perl -lane ‘/gene_id “([^”]+)”/ or die “no gene_id on $.”;print join “\t”, (@F[0,1,2,3], $1)’ |sort -k1V -k2n -k3n >> Mus_musculus.GRCm39.104.gtf.rRNA.refflat
Statistics
Summary
Assembly | GRCm39 (Genome Reference Consortium Mouse Reference 39), INSDC Assembly GCA_000001635.9, Jun 2020 |
Base Pairs | 2,728,222,451 |
Golden Path Length | 2,728,222,451 |
Assembly provider | Genome Reference Consortium |
Annotation provider | Ensembl |
Annotation method | Full genebuild |
Genebuild started | Aug 2020 |
Genebuild released | |
Genebuild last updated/patched | Mar 2021 |
Database version | 104.39 |
Gencode version | GENCODE M27 |
Gene counts (Primary assembly)
Coding genes | 22,468 (incl 253 readthrough) |
Non coding genes | 16,060 |
Small non coding genes | 5,526 |
Long non coding genes | 9,972 (incl 76 readthrough) |
Misc non coding genes | 562 |
Pseudogenes | 13,650 (incl 1 readthrough) |
Gene transcripts | 142,434 |
Other
Short Variants | 83,318,087 |
Structural variants | 791,878 |
Index and modification
Index
Index software | File list |
STAR rsem-prepare-reference | chrLength.txt chrName.txt chrNameLength.txt chrStart.txt exonGeTrInfo.tab exonInfo.tab geneInfo.tab Genome genomeParameters.txt GRCm39.104.genome.chrlist GRCm39.104.genome.grp GRCm39.104.genome.idx.fa GRCm39.104.genome.n2g.idx.fa GRCm39.104.genome.seq GRCm39.104.genome.ti GRCm39.104.genome.transcripts.fa Log.out SA SAindex sjdbInfo.txt sjdbList.fromGTF.out.tab sjdbList.out.tab transcriptInfo.tab |
gtfToGenePred samtools perl | Mus_musculus.GRCm39.104.gtf.gtf.refflat Mus_musculus.GRCm39.104.gtf.rRNA.refflat Mus_musculus.GRCm39.dna.primary_assembly.fa.fai sizes.genome |
Bundle files
Type | File list |
NA |