I got RNASeq data in several samples. I checked the FastQC, seems the read quality are good (Hiseq 2000). But the problem is many reads are mapped to intronic region, and the regions have no any reference exons there (Refseq, ensembl, gencode). We don't know what they are. We guess the problem happend in library preparation, the concentration was low. Now the data has come out and we can't re-sequencing, so we want to remove the reads mapped to intronic region, is there a method to do that? Or anyone have an idea about the intronic reads. Thanks.
•
written
by
camelbbs
2 answers
You can easily use BEDOPS to solve this problem quickly. It includes bedops and various conversion scripts for putting data into BED format, which bedops can process.
Assuming your reads are in BAM format:
$ bam2bed < reads.bam \
| bedops --not-element-of -1 - introns.bed \
> reads-not-in-introns.bed
The file reads-not-in-introns.bed is a sorted BED file containing all reads that do not overlap intronic elements.
You can then pass this result to bedmap to do counting of reads over other region sets (whole-genome or subsets).
Note that we assume your introns are in BED format and are sorted, e.g.:
$ sort-bed unsorted-introns.bed > introns.bed
Alternatively, if your introns are in some other format — say, GTF — then BEDOPS conversion scripts will losslessly turn them into sorted BED, e.g.:
$ gtf2bed < introns.gtf > introns.bed
•
written
by
Alex Reynolds
Log in to answer this question.
-
Low Mapping Rate with Kallisto on RNA-seq Data
written by meetmet •
-
Requesting Insights on Potential Complex Structural Event Indicated by IGV
written by pragnapcu •
-
Receiving more reads from exons than whole genes using featurecounts
written by till.schwaemmle •
-
Kmer content failed in FastQC analysis
written by josmantorres •
-
Starrseq data how can I characterise reads that maps out of the targeted regions.
written by morovatunc
-
Fastq reads mapped to the desired position on reference genome
written by Louis Kok •
-
shall we do PCA analysis in target region sequencing project
written by yangchaogab •
-
Quality Drop In The Reverse Read Of Paired End Library (Illumina Hiseq 1000)
written by ancient_learner
-
How To Calculate The Rpkm For Isoforms
written by camelbbs
-
Why Are There Many Rna-Seq Hits To Intronic Regions?
written by lh3
Ask yourself: "why should I remove intronic reads?" Do you want to remove outcome that you do not understand, until your experiment fits your expectations?
What does low concentration have to do with getting unwanted reads, what 'makes up sequences' that are not real in case of low concentration? See also Why are there many RNA-seq hits to intronic regions? Intronic sequences might be novel transcripts, remains of nascent RNA, lincRNA, antisense RNA, if close to exons, wrong exon boundaries in the annotation.