Hi all,
I am looking for your help about tools/pipeline for allele specific expression analysis. From various papers, I have come up with this pipeline:
1) Call variant with GATK using WGS samples
2) Create 2 transcriptomes using these variants and my annotation: so I have one transcriptome for allele 1 and one transcriptome for allele 2
3) Map my RNA-seq data against this pool transcriptome using bowtie2
4) and get count with ASE-TIGAR (http://nagasakilab.csml.org/ase-tigar/) per gene per haplotype
5) Perform downstream analysis with DEseq2.
It works and I got some results but I am wondering if you know any other tools than ASE-TIGAR to get the counts just for comparison/validation purpose.
Also, if you have other pipeline recommendation, I would be happy to hear ! Thanks for your help.
0 answers
No answers yet.
Log in to answer this question.
-
What are the best tools for quantifying allele-specific expression from bulk RNA-seq data these day…
written by Paulo •
-
Extract genes from GTF with variants from VCF
written by Picasa
-
Parameter to filter a VCF of SNPs?
written by Picasa
-
Differences in RNAseq Variant Calling and Allele Specific Expression
written by serpalma.v •
-
allelic specific expression with Kallisto
written by Assa Yeroslaviz
-
How to get counts data for allele specific expression analysis (ASE) from aligned bam files (RNA-se…
written by Naresh D J •
-
Input file MULTIPOOL - reference allele counts from mpileup/bam
written by pedroamandrade •
-
Feasible Annotation of 4 million unique SNPs
written by serpalma.v •
-
CNV analysis on Illumina NGS data with PLINK or other toolsets
written by plague •
-
Allele specific expression (ASE) using unfiltered reads from RNA-seq data?
written by mbk0asis
For the part of DESeq2, I would recommend having a look at this page http://rpubs.com/mikelove/ase
Thanks for your link. Do you have experience with ASE study ?
a bit, though the field itself is just emerging with no well-accepted standards.
You could take a look at this quantification tool