SEQanswers

Go Back   SEQanswers > General



Similar Threads
Thread Thread Starter Forum Replies Last Post
RNA-Seq: SOAPsplice: Genome-Wide ab initio Detection of Splice Junctions from RNA-Seq Newsbot! Literature Watch 0 02-04-2012 02:00 AM
RNA-Seq: RNA-Seq Analysis of Gene Expression and Alternative Splicing by Double-Rando Newsbot! Literature Watch 0 03-03-2011 02:00 AM
RNA-Seq: MapSplice: Accurate mapping of RNA-seq reads for splice junction discovery. Newsbot! Literature Watch 2 10-14-2010 08:35 AM
Where is the mRNA-seq reads and KSHV virus or cell line BC3's alternative splice? jiexu Bioinformatics 1 10-12-2010 09:39 AM
Alternative Hits in BWA generated .sam jtladner Bioinformatics 1 06-09-2010 02:07 PM

Reply
 
Thread Tools
Old 07-17-2011, 07:31 AM   #1
Louis_Lemire
Junior Member
 
Location: Ft. White, FL

Join Date: Jun 2011
Posts: 5
Question Alternative splice or RNA-seq generated error?

I am about as newbie as a newbie can get so to be honest I'm a little reluctant to post - however, I have discovered a result looking at the raw Illumina reads that is not readily answerable this early in my RNA-seq workflow, and, with my limited knowledge at this point, thought I might ask to see if I could get an opinion on the below results.

I am particularly interested in an hypothetical unannotated paralog/alt splice that may not align with my genome so I am spending quite a bit of time perusing the raw illumina reads in order to take a close look at some of the more conserved regions of my research proteins looking for paralogs, etc. (as well as to get a 'feel' for the raw reads, how they behave, etc, on manual queries). After I finish with this preliminary analysis I have a few hundred hours of learning before I can comment with any confidence on sequence assembly matters - I enjoy computers but I am far removed from a Linux wizard.

Below is a result I found generated from high quality reads (>Q30) which suggests an alternative splice. In the code box below there are three lines:

Line 1: partial exon 2 of one of my research proteins
Line 2: raw Illumina reads linked by grep query, all have >Q30, and all cross the putative splice site
Line 3: partial exon 7 of the same protein in Line 1

The <....> bracket indicates the beginning of an intronic sequence at the end of exon 2.

Code:

                              ********** ***::****:*
e2                         ...RGHTGLFAGG<ASTYQVGLELC...>
     ...GHALLFRTSVMAKVEIQAVSTCRGHTGLFAGG<ASTFHVGLEAC...>
e7   ...GHALLYRTTVMAKLEIQAVSTCR...      <--- intron --->
        *****:**:****:*********
The above result appears to be an alternative splice. However, I was wondering if it may be an error generated by RNA-seq preparation of exp material, i.e., two pieces of DNA randomly cut and joined. There were about 10 copies of the middle region above all yielding high quality reads and all crossing an apparent splice site.

Q: What is the likelyhood that the above is real and not a machine artifact?

Last edited by Louis_Lemire; 07-17-2011 at 08:36 AM. Reason: grammer
Louis_Lemire is offline   Reply With Quote
Old 07-22-2011, 02:38 AM   #2
Louis_Lemire
Junior Member
 
Location: Ft. White, FL

Join Date: Jun 2011
Posts: 5
Default

I may have found an explanation for this strange exon7-exon2 splice. Li et al. (2008) discuss a statistical approach towards identifying the degree of alternative splicing in a differential gene expression paper. [1] Li points out that the splicing of exons in reverse order is 'impossible' [2,3] and uses these rare events as measures of alternative splice false discovery events. Li found that these type of splicing events amount to roughly 1% of mapped junctions.

Presumably during the workup of the cDNA for Illumina reading there is a low probability that the DNA can form a hair-pin turn back on itself and purportedly recombine - a rare event. That this happened with my research protein was coincidental.

[1] Hairi Li et al (2008) Determination of tag density required for digital transcriptome analysis: Application to an androgen-sensitive prostate cancer model. PNAS 105(52):20179-20184
[2] D. L. Black et al. (2003) Mechanisms of alternative pre-messenger RNA splicing. Annu Rev Biochem 72:291-336.
[3] J. M. Johnson et al. (2003) Genome-wide survey of human alternative pre--mRNA splicing with exon junction microarrays. Science 302:2141-2144
Louis_Lemire is offline   Reply With Quote
Reply

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 04:10 AM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2018, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO