I want to run HTSeq on count number of mapped mRNA reads on Human genes (GENCode GTF file). However, I have to modify this gtf file because I want to allow a search region of an offset of -200 bp upstream of each gene. I think I need to add exon of length 200bp to each gene, so, HTSeq can count the reads that overlap such regions. Do you know if there is some script (code) which can do that?
-
-
Just extend the bounds of each exon #1. I'm not aware of a premade script to do this, but you could write one in python or R (import the GTF as a GRanges object, split by transcript, and apply an appropriate function) easily enough. -
Thank you very much @dpryan. I think I will write this code. I guess it is simple but needs more careful analysis. This is because we need to check for each gene (or transcript) if the new upstream region will not hit any upstream gene. If such case happens, then we need to shrink the upstream region accordingly. Is this right? or we let genes overlap with each other (then HTSeq will identify reads as ambiguous)Comment
-
There's no reason to shrink regions if they overlap another gene. If you're going to grow a genes extent and that causes an overlap then so be it. Having a few more ambiguous alignments in such cases is a non-problem.Comment
-
Thanks @dpryan. I agree with you that we do not need to shrink, however, I suggested this to avoid eliminating reads from being counted by HTSeq count. Actually, I do not like the ambiguity approach in HTSeq count because it artificially changes the way the expression count is defined. I just posted a question about this.Comment
-
There's nothing artificial about how htseq-count handles that. It's the correct solution for the downstream analyses for which it's intended.Comment
-
I agree with you if you mean that HTSeq-count is intended to generate expression matrix for detection of differentially expressed genes. Unfortunately, this is just one single applications, while, many applications need expression matrix that reflect the actual expression mechanism. Please correct me if you think this is not true and/or not accurate.Comment
-
Correct and there's no need to bother modifying htseq-count when there are so many existing tools that are intended for those purposes.Comment
-
Thanks very much @dpryan for your replies and the suggestionsComment
-
CRISPR/Cas9 sparked the gene editing revolution for both research and therapeutics.1 But this system still showed severe issues that limited its applications. The most prominent were the heavy reliance on PAM sequences, delivery limitations, double-stranded breaks that prompt unintended edits and cell death, and editing inefficiency (both in targeting and in knock-in reliability).
Despite this, “CRISPR helped turn genome editing from a specialized technique into...07-31-2026, 11:01 AM -
Proteomics platforms are evolving rapidly, with advances in mass spectrometry and affinity-based approaches expanding what researchers can detect and at what scale. As the field moves toward deeper proteome coverage and clinical applications, scientists face an increasingly complex landscape of tools. This article will explore how researchers are navigating these choices to find the right platform for their work.
The systematic characterization of the human proteome has...07-20-2026, 11:48 AM -
Genomics studies in neuroscience face a special challenge due to the brain’s complexity and scarcity of samples. Mapping changes in cell type and state using conventional next-generation sequencing methods remains challenging. Advances in technologies like single-cell sequencing, spatial transcriptomics, and long-read sequencing have opened the door to deeper studies of the brain and diseases like Alzheimer’s, amyotrophic lateral sclerosis (ALS), and schizophrenia.
...07-09-2026, 11:10 AM
Comment