Allele Specific Expression Analysis
2
0
Entering edit mode
6.2 years ago

Hey All,

Does anyone have any experience in Allele specific gene expression analysis? Any recommended tools and pipelines? And any idea on how to retain only RNA-seq reads overlapping known SNP sites that differ between the two distinct parental genomes? I have fair knowledge of RNA Seq and SNP analysis. But ASE is nowhere explained properly. Anything would be helpful

EDIT: I have created two pseudo parental genomes and mapped my samples onto both of them using STAR and with options --outFilterMismatchNmax 0 ( Allowing for 0 mismatches ) and --outFilterMultimapNmax 1 (Allow no multi-mapped reads ). This ensures allele-specific mapping of most SNP-containing reads to only one strain-specific genome. I am interested only in the X chromosome ( X linked genes ). I don't know what to do after this.

Susmita

SNP Allele specific RNA-Seq X chromosome STAR • 4.8k views
ADD COMMENT
0
Entering edit mode

Hi, After alignment you can use bedtools to intersect your bam or bed on SNP coordinates and extract only those reads which overlap SNPs.

ADD REPLY
0
Entering edit mode

Hi Susmita, Have you by any chance managed to analyse your allele-specific RNA-seq data? best TA

ADD REPLY
1
Entering edit mode

Yes I did managed to do the analysis and some more. But its not allele specific RNA-Seq data. I do allele specific on RNA-Seq data. Any questions?

ADD REPLY
0
Entering edit mode

yes that is what i meant .. Do you mind sharing the workflow you have used with me?

Tamer

ADD REPLY
1
Entering edit mode
6.2 years ago
Tm ★ 1.1k

You can have a look at ASEReadCounter of GATK or SCALE.

ADD COMMENT
0
Entering edit mode

I do not have single cell RNA Seq data. But I'll look into the ASEReadCounter of GATK for sure!

ADD REPLY
0
Entering edit mode

You can also refer this article where they have used STAR for mapping purpose followed by ASEReadCounter tool for ASE analysis. Here they have used human lymphoblastoid cell line RNAseq data and not single cell RNAseq data.

ADD REPLY
0
Entering edit mode

does ASEReadCounter work for INDELs? as I tried using it and did not get any results

ADD REPLY
0
Entering edit mode

Try HaplotypeCaller from GATK or VarScan2.

ADD REPLY
1
Entering edit mode
6.2 years ago

Also take a look at PHASER software https://www.nature.com/articles/ncomms12817. The authors are famous in the field.

ADD COMMENT

Login before adding your answer.

Traffic: 2178 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6