This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Chicken embryonic heart RNA Seq reads are not mapping to exons- need some advice.

Hi Biostars community, I am making my first foray into RNA Sequencing on chicken tissues and I have run into a problem. We have conducted an experiment on isolated AV Cushion cells and performed an RNA Seq to gather some exploratory data. We prepped rna using the RNEasy plus kit, which has a genomic dna removal column, and at the end we judged the RNA Quality to be good (RINs all were >9.7). As far as I know, library preparation and sequencing went off went off without a hitch (my collaborators performed a polyA enrichment and then sent the samples for sequencing on an Illumin Novaseq 6000 platform). Alignment of my reads to the Galgal6 genome using STAR worked ok (~63% of reads were unique paired aligned reads, and about 37% were not mapped). However, when I looked at the coverage breakdown, only 3% of the reads mapped fully within an exon (about 18% mapped partly within an exon, 43.15% mapped fully within an intron, and 33.65% was mapping to intergenic DNA). Does anyone have an idea why we are getting so few of our reads mapping to exons? Can anyone who has done chicken RNASeq before share their thoughts/experience? Potentially a problem with polyA enrichment maybe? I am new to the game and need help!

heart chicken rna-seq

Do you have DNA contamination in this data by chance?

We shouldn't, we prepared our RNA using the RNEasy Plus kit which has a genomic DNA eliminator column. That said, I am keeping an open mind- how could I tell if I did?

Since you said

43.15% mapped fully within an intron, and 33.65% was mapping to intergenic DNA

Did you browse the alignments and see if the reads are scattered all over. Ideally in RNAseq they should be piling up under exons.

0 answers

No answers yet.

Log in to answer this question.