SEQanswers

SEQanswers (http://seqanswers.com/forums/index.php)
-   General (http://seqanswers.com/forums/forumdisplay.php?f=16)
-   -   how to visualize locations of lncRNA close to my gene of interest (http://seqanswers.com/forums/showthread.php?t=48619)

angel-sakura 11-30-2014 03:38 AM

how to visualize locations of lncRNA close to my gene of interest
 
Hi, everyone, I am looking for the lncRNAs of some specific genes.
Q1: Now I have extracted all predicted lncRNAs (using PLEK) that close to my target genes (500 KB upstream or downstream). I can do blastn to get the location of those lncRNAs on gene scaffold I extracted from genomic sequence (my target gene +/-500KB). But I don't know how to visualize them. Could someone please show me the software or protocol?
Q2: Is there any lncRNA database or lncRNA database with consensus sequences that I can blast with?
Q3: Is there any software or command that I can predict structure/function/lncRNA-protein, lncRNA-DNA,lncRNA-microRNA interaction of lncRNAs (with multiple sequences as input)?

Thank you very much!
Y

dpryan 12-01-2014 12:58 AM

1. I'd use IGV and make either a BED or GTF file describing the lncRNA locations in the scaffold. If you add your gene of interest to it too then you can easily visualize where things are relative to each other (you can also do this in R with packages like Gviz, but I think just using IGV will be a bit easier for your needs).
2. RNAcentral has the lncRNA sequences from Ensembl and lincRNAdb, so have a look at it. It's interface isn't great at the moment, so you'll have to download and parse the fasta file (I actually have some tools to do this that I need to post on github).
3. No clue, hopefully someone else knows.

angel-sakura 12-01-2014 02:22 AM

Thank you very much for your kind reply!
Could you also recommend some commands or softwares that I can use to transfer lncRNAs (fasta format) into BED or GTF format?

dpryan 12-01-2014 02:26 AM

The simplest thing is to probably use blast to align them to you either the genome or the targeted region and then use some of the replies in this thread over at Biostars to convert that to GFF/GTF. I should note that I've not personally done this, so I can't say how easy/difficult it is...but it should at least work well in theory (famous last words!).

angel-sakura 12-02-2014 04:47 AM

lol, I haven't done this before either, but it is worth trying it.
Thank you !

docnj 12-28-2014 06:35 AM

Hi angel-sakura, You said that you have extracted all predicted lncRNAs that close to the target genes by using PLEK. I'm confronted with the same problem myself, too. I have one interested target gene, but don't know how to find lncRNAs close to it. May I ask how you conducted the prediction through using PLEK? If you can also provide procedures in details, that would be appreciated very much!! Thank you!

angel-sakura 01-02-2015 09:24 PM

Quote:

Originally Posted by docnj (Post 157194)
Hi angel-sakura, You said that you have extracted all predicted lncRNAs that close to the target genes by using PLEK. I'm confronted with the same problem myself, too. I have one interested target gene, but don't know how to find lncRNAs close to it. May I ask how you conducted the prediction through using PLEK? If you can also provide procedures in details, that would be appreciated very much!! Thank you!

Hi, sorry for the late reply, I was on holidays. :)
Check this link:http://sourceforge.net/projects/plek/files/
You can downloand the PLEK and its manual from above link. The PLEK is straight forward, you just input the fasta file, then it will produce an output containing coding and non-coding genes.
As for finding lncRNAs close to your gene of interet. The steps are: 1. extract the upstream and downstream sequences (such as 500kb) of your target gene. 2. Blast RNAseq data to fragment extracted from step 1. 3. predict the blast result with PLEK.

Cheers


All times are GMT -8. The time now is 12:48 PM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2021, vBulletin Solutions, Inc.