SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Reply
 
Thread Tools
Old 02-08-2021, 06:41 AM   #1
icd20
Junior Member
 
Location: London

Join Date: Dec 2020
Posts: 1
Default Help with ERV analysis

Hi,
Very new to this! Apologies if not enough information given.
I'm trying to count the change in ERV expression before and after treatment of human cell lines with a drug. I have been using the gEVE database, which provides a GTF annotation for human genomes (http://geve.med.u-tokai.ac.jp/download/).
My pipeline is: generate genome indices using STAR --> trim reads (fasta files) using trimmomatic --> sort, remove blacklisted regions --> align trimmed reads to genome using STAR again --> count reads using featureCounts.

All seems fine except counting reads: I get a very low percentage of assigned reads (usually 0.1%) and these are mainly due to "Unassigned_NoFeatures". This is my input for featureCounts:
fc_ERV <- featureCounts(files=filenames, annot.ext="Hsap38.geve.v1.gtf", isGTFAnnotationFile=TRUE, GTF.featureType="CDS", countMultiMappingReads=TRUE, genome='Homo_sapiens.GRCh38.dna.primary_assembly.fa', isPairedEnd=TRUE,nthread=20)
I have used "GTF.featureType="CDS"" because the GTF from gEVE does not have an exon column, only CDS.
Any ideas on what I'm doing wrong? I am new to Bioinformatics, any help would be much appreciated.

Last edited by icd20; 02-08-2021 at 06:53 AM.
icd20 is offline   Reply With Quote
Reply

Tags
erv, featurecounts, gtf annotation file, rsubread

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 11:18 AM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2021, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO