Extracting promoter sequence with GenomicFeatures ?
Entering edit mode
2.6 years ago
lokraj2003 ▴ 90

I am trying to extract promoter sequence of a non-model organism (Orf virus). I used flank function from GenomicFeatures to extract 50 base pairs upstream of start codon.

promoters <- flank(gr, 100)

Then I used those coordinates on the costum BSgenome Orf virus genome package that I made by myself to extract DNA sequences of the 50 bp of each gene.

promoter_seq <- getSeq(orf_genome, start=promoters$start,end=promoters$end)

Since, this function takes strand into consideration, it is taking 50 bp from the upstream of the positive sense gene and 50 bp downstream of negative sense gene (50 bp after start codon). Is this approach correct to extract promoters ? I don't know the exact location of promoters, so, I am just extracting 50 bp now. Also, I guess the promoter sequence extracted from the negative strand has to be reverse complemented. Is there any way to do something like "if the strand is negative , reverse complement the DNA sequence" ?

GenomicFeatures Bioconductor • 761 views

Login before adding your answer.

Traffic: 1737 users visited in the last hour
Help About
Access RSS

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

Powered by the version 2.3.6