I am working with a few novel transcripts of genes- before I confirm their existence experimentally, I would like to perform some bioinformatic analysis. I have already considered coding potential, protein domain prediction, transcription factor binding sites, sequence homology, and RNA secondary structure (still a little unsure how to use this one). These transcripts were discovered using RNA-Seq. Are there any other elements of genes/confirmed transcripts that I should look for in the sequence of my transcripts and corresponding softwares? (I can find the software myself if necessary, but I have run out of characteristics to search for). I would like to characterize these transcripts structurally and functionally as completely as possible, including potential protein function, mRNA degradation, etc - some new features for what to look into would be appreciated.
What are some (bioinformatic) methods to characterize potentially novel gene transcripts?
$\begingroup$ The current version of your question is unclear. The RNA-Seq reads you are analyzing clearly exist (otherwise there would be nothing to analyze) yet your question makes it sound like the transcript does not exist. If the transcript does not exist then what is your alternative hypothesis? Where could the reads come from? $\endgroup$– mdperryJul 8, 2015 at 22:16
$\begingroup$ Thanks for your response! Let me clarify: these transcripts are created by alternative splicing, and are unique in that they contain retained introns. Their expression levels are quite low, and I have been told that there is significant likelihood that they are merely artifacts ( I am not too familiar with rna-seq to understand where this may arise) -"novel" transcripts with intron retention apparently do not often hold up to experimental validation. I want to gain evidence that they exist in silico before I move on to qpcr. $\endgroup$– pMarkovJul 8, 2015 at 22:55
$\begingroup$ Furthermore, I would like to characterize these transcripts as much as possible - gain insight into whether they would just be noncoding transcripts or be of some interest- I attempted to look into what types of protein domains may arise, whether the transcript retained introns have coding potential, and new tfbs. Essentially I want to gain as much insight into the structure and potential function of this transcript as possible, using only the sequence. I am merely a student, so please ask for further clarification! *I edited to question to match my responses here. $\endgroup$– pMarkovJul 8, 2015 at 22:58
$\begingroup$ @mdperry- As an example, I am currently looking into potential miRNA binding site introduced to the transcript through intron retention - I am open to investigation of any interesting features that these transcripts may bring! $\endgroup$– pMarkovJul 8, 2015 at 23:06
$\begingroup$ You might want to check the genomic sequences and how the transcript is spliced so that RNA-seq data could be confirmed. You may also find alternative splicing isoforms. You could check poly A signals in the genomic sequence. $\endgroup$– 243Jul 9, 2015 at 1:58
It sounds like you have considered most of the obvious alternatives (and thank you for clarifying the question). I suppose the first question in an alternatively spliced transcript with a retained intron is whether the open reading frame of the protein is maintained. If there is a termination codon that now becomes in-frame due to the intron then the protein would ordinarily be truncated. The mRNA could also get targeted for degradation due to the nonsense-mediated decay pathway (i.e., from introducing an early stop-codon).
If the alternatively-spliced transcript lacks any significant ORFs, or if there is an ORF but no suitable translation initiation codon then you have entered the realm of non-coding RNAs, of which there are two loose categories, microRNAs (there are several classes) or long non-coding RNAs. most miRNAs have recognizable sequence motifs, may be contained in mir-Base, and may be complementary (in part) to a regulatory target elsewhere in the genome. lncRNAs are less well-defined. Perhaps the best criteria are transcripts that are stable enough to be detected, but lack noticeable protein coding features. I don't even think there is a consensus on the minimum length of a transcript to be considered an lncRNA.
What sort of databases have you searched? For example the ancient dbEST contains all kinds of short cDNA reads from all over the planet, dating back to 1992 (or so). Evolutionary conservation of exons and splice-sites, or reading frames can all be used to support the hypothesis that a transcribed region has a biological function.
For a really deep fishing expedition you can use tblastn, etc. to take a potential translated ORF and search all 6 conceptually translated reading frames of all the sequences in a database--it takes longer, and there can be many spurious matches that you have to search through, but if you are studying something that has never been annotated it might be worth considering.
RE: folding RNA (secondary structure) every sequence can be folded into some kind of structure (try it and see), but structurally conserved compact loops typically fall into one of a small class of selected families (like GNRY). Stems are easy to find, but if the single-stranded loop is huge, how likely is that structure to fold up in real time?
RE: TFBS not clear to me at all what you mean by this. Transcription factors bind dsDNA not ssRNA. There are families of protein domains that bind RNA but not typically with such sequence-specificity. So unless this alternatively spliced transcript is being transcribed from some type of alternative promoter, I am not sure where you would even be looking for TFBS. It's true that some large introns can indeed contain regulatory regions for the gene as a whole (cf. immunoglobulin heavy chain enhancer is conveniently stored in an intron), but it is not clear to me if that would have an affect when transcribed into RNA.
I do like your idea that a retained intron could now contain a binding site for a regulatory miRNA, but if I understand the theory this would lead to the transcript being degraded--so it would be more difficult to detect.
$\begingroup$ Thanks! Again, I believe I have not shared all necessary information; in some of the transcripts, the alternative spicing generates an alternative 5' UTR, in which I am searching for TFBS (these alternate 5' UTR's appear to be retained introns). The miRNA idea comes from that these transcripts are actually differentially expressed in the disease phenotype I am analyzing, so I was looking into introduced miRNA binding sites as a potential reason. The ORF analysis is interesting-I had not considered the meaning behind a lack of orfs within the retained intron. Thanks for your time! $\endgroup$– pMarkovJul 12, 2015 at 0:57
For identifying function do a homology search. There is little functional annotation of lncRNAs. So homology based information can be obtained only for protein sequences. So you can try these:
- Check the coding potential. Find ORFs (perhaps set a minimum length cutoff). To be stringent you can also check for Kozak consensus sequences (for eukaryotes) in these novel transcripts. This program apparently includes the Kozak rule in start codon prediction.
- Translate the RNA and run psi-BLAST. Psi-BLAST is better than normal BLASTp in identifying distant homologs.
- You can obtain the GO (functional annotation) data for these homologs and relate them to your novel transcript.
- For lncRNAs, you can find if they overlap with any other known transcript. Many lncRNAs seem to have overlapping loci with their target genes (both sense/antisense).
- Search for RNA motifs. This article reports an extensive study towards identification of different RNA motifs (RBP binding sites) in humans and Drosophila. Other interesting motifs include quadruplex forming motifs and ARE.
- If you have the genome sequence then you can map your transcripts to the genome and identify the introns. Actually your RNAseq assembly, if reference guided, would produce a GTF file. You can extract introns using the GTF. You may check for potential miRNAs in these introns, based on stem loop predictions. However it is better to have a small-RNAseq data for predicting miRNAs reliably.
- As you had already considered, you can check for miRNA binding sites. For this purpose miRanda is the best for novel transcripts (TargetScan is a pain). You can use RNAhybrid too and take a consensus of miRanda and RNAhybrid results.
- You can do other kinds of analysis based on just the sequence of the transcripts. These include checking for codon usage (for potential protein coding RNAs), GC richness (comparison with the known transcripts) etc.