This is a test version of Biostars. For the public version, visit https://www.biostars.org.
LncRNA pipeline

HI ALL, I am working with lncrna, my pipeline is mapping with star and then sorting, folloed by stringtie and from the output using awk command, extract lncrna: -

awk '($3 == "transcript" && $5 - $4 > 200) { exon_count = gsub(/exon_number "[0-9]+"/, "&"); } (exon_count > 2) { print }'
lncrna rna-seq ngs
• 2,121 views • link
updated by GenoMax GenoMax • written by Researcher • Researcher

2 answers

Subhikhsa, subtle changes

=200

and exon count could be equivalent or more than 3 to avoid intergenic boundaries and capture lincRNAs

Can you please send the screenshot?

Prash

PS: Please don't type in caps. It means SHOUTING in net jargon

• 0 views • link
written by Prash • Prash

okay thank you sir.

1   StringTie   transcript  14529   29358   1000    -   .   gene_id "STRG.1"; transcript_id "STRG.1.1"; cov "9.250028"; FPKM "1.439988"; TPM "3.138408";
1   StringTie   transcript  14529   24901   1000    -   .   gene_id "STRG.1"; transcript_id "STRG.1.2"; cov "15.766621"; FPKM "2.454452"; TPM "5.349400";
1   StringTie   transcript  135141  135895  1000    -   .   gene_id "STRG.2"; transcript_id "STRG.2.1"; reference_id "ENST00000494149"; ref_gene_id "ENSG00000268903"; cov "1.669579"; FPKM "0.259910"; TPM "0.566465";
1   StringTie   transcript  137682  137965  1000    -   .   gene_id "STRG.3"; transcript_id "STRG.3.1"; reference_id "ENST00000595919"; ref_gene_id "ENSG00000269981"; cov "1.830701"; FPKM "0.284992"; TPM "0.621132";
1   StringTie   transcript  164282  168955  1000    -   .   gene_id "STRG.4"; transcript_id "STRG.4.1"; cov "1.205568"; FPKM "0.187676"; TPM "0.409033";
1   StringTie   transcript  184957  197026  1000    -   .   gene_id "STRG.5"; transcript_id "STRG.5.1"; cov "3.289948"; FPKM "0.512159"; TPM "1.116234";
1   StringTie   transcript  185217  195411  1000    -   .   gene_id "STRG.5"; transcript_id "STRG.5.2"; reference_id "ENST00000623083"; ref_gene_id "ENSG00000279457"; ref_gene_name "WASH9P"; cov "8.402778"; FPKM "1.308094"; TPM "2.850948";
1   StringTie   transcript  185217  187958  1000    -   .   gene_id "STRG.5"; transcript_id "STRG.5.3"; cov "4.316052"; FPKM "0.671897"; TPM "1.464378";
1   StringTie   transcript  185217  199874  1000    -   .   gene_id "STRG.5"; transcript_id "STRG.5.4"; cov "2.504329"; FPKM "0.389859"; TPM "0.849685";
1   StringTie   transcript  185217  199874  1000    -   .   gene_id "STRG.5"; transcript_id "STRG.5.5"; cov "4.394231"; FPKM "0.684067"; TPM "1.490903";
1   StringTie   transcript  185217  197026  1000    -   .   gene_id "STRG.5"; transcript_id "STRG.5.6"; cov "1.085559"; FPKM "0.168993"; TPM "0.368315";

these are few lines from the lncrna extracted gtf file, that i obtained after using the above command .

• 0 views • link
updated by GenoMax GenoMax • written by Researcher • Researcher

I also have one more doubt. so if this gtf file is right, to find the expression of lncrna, ill have to run htseq for the mapped bam file along with this lncrna gtf file to get lncrna expression count, which i can further use for dseq2. is that right?

• 0 views • link
written by Researcher • Researcher

This is the pipeline I suggest where we benchmarked . However, it all depends on the questions you are asking! Is it towards the regulatory potential of lncRNAs or lncRNA x gene interatcions etc.

• 0 views • link
written by Prash • Prash

Thank you for sharing the pipeline sir. My end goal is lncrna and gene interaction

• 0 views • link
written by Researcher • Researcher

Log in to answer this question.

More posts like this