How can I get counts from transgenic expression in RNA-seq data?
1
0
Entering edit mode
2.1 years ago

Hello all,

This is probably a pretty simple question, but I've never done an analysis with transgenic expression involved before and I want to make sure I do it correctly. Basically, I have RNA-seq data, a genome and annotation files for the host, and the sequence of a transgene and its expression vector. The transgene is exogenous and therefore not in the host genome. Mapping the reads to the host transcriptome is straightforward enough, but what would be the best method by which to get counts data for the exogenous transgene (along with any portions of the vector itself which may be present)?

Thanks,

Erik

RNA-seq • 1.2k views
ADD COMMENT
3
Entering edit mode
2.1 years ago

You should add the transgene as a chromosome to the FASTA you used to generate your reference index. You could also add it to your GTF but I think it would be considerably more work. Hopefully you are mapping with something like STAR which should accommodate this well.

ADD COMMENT
1
Entering edit mode

Alternatively, you can take the unmapped reads from the host genome and map it to a reference genome of your expression vector.

ADD REPLY
0
Entering edit mode

Great, thank you! I'm using STAR so this should work just fine. Again, thanks so much!!

ADD REPLY
0
Entering edit mode

Oh, and one more question about this. I have to add it to the GFF file if I want STAR to map it, don't I?

ADD REPLY
1
Entering edit mode

Not necessarily unless you are using STAR to do your gene/exon counts. Otherwise, you could just count the alignment to the specific chromosome.

ADD REPLY

Login before adding your answer.

Traffic: 2521 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6