Unconfigured Ad

Collapse
X
 
  • Time
  • Show
Clear All
new posts
  • qoiopipq
    Junior Member
    • Apr 2016
    • 7

    #1

    Tissue specific 3' UTRs using DEXSeq

    I used IsoSCM to identify alternative 3' UTRs with different lengths in different tissue types(output files are GTFs containing tissue specific 3' UTRs), then I want to quantify the differential usage of these alternative 3' UTRs across different tissue types. I know that neither Ensembl nor UCSC GTF reference file contains the tissue specific 3' UTRs, so I tried to re-annotate IsoSCM output GTF files with reference GTF. The idea is that I want to get a GTF file containing all tissue specific 3' UTRs, which can be used for the following dexseq_count step.

    But, if multiple 3' UTRs/terminal exons (same start site, different stop sites) are annotated with same transcript ID, dexseq_prepare_annotation will only include the longest one in the GFF output.
    Here is an example of two possible 3' UTRs of gene Fubp1 in thymus GTF file:
    Code:
    chr3	sol	exon	152232395	152233253	.	+	.	gene_id "Fubp1"; gene_name "Fubp1"; p_id "P24872"; transcript_id "NM_057172"; tss_id "TSS5999";
    chr3	sol	exon	152232395	152232485	.	+	.	gene_id "Fubp1"; gene_name "Fubp1"; p_id "P24872"; transcript_id "NM_057172"; tss_id "TSS5999";
    Two different 3' UTRs of gene Fubp1 in liver GTF file:
    Code:
    chr3	sol	exon	152232395	152233965	.	+	.	gene_id "Fubp1"; gene_name "Fubp1"; p_id "P24872"; transcript_id "NM_057172"; tss_id "TSS5999";
    chr3	sol	exon	152232395	152233885	.	+	.	gene_id "Fubp1"; gene_name "Fubp1"; p_id "P24872"; transcript_id "NM_057172"; tss_id "TSS5999";
    One more 3' UTR of gene Fubp1 in spleen GTF file:
    Code:
    chr3	sol	exon	152232395	152232532	.	+	.	gene_id "Fubp1"; gene_name "Fubp1"; p_id "P24872"; transcript_id "NM_057172"; tss_id "TSS5999";
    After merging them and put through dexseq_prepare_annotation step, I only got one exon to represent 3' UTR:
    Code:
    chr3	dexseq_prepare_annotation.py	exonic_part	152232395	152233965	.	+	.	transcripts "NM_057172"; exonic_part_number "018"; gene_id "Fubp1"
    So, I am wondering how can I include all alternative 3' UTRs in the GFF file, then I can use it to count for the usage of these alternative 3' UTRs ??

    Cheers

Latest Articles

Collapse

  • kinpict.com
    Reply to 16S/ITS Profiling vs. Shotgun Metagenomics: What Changes in Your Workflow?
    by kinpict.com
    The part about this topic feels surprisingly grounded, and the details on the page is the bit that made me slow down. I’d keep <a href="https://kinpict.com/">ai family photo generator</a> bookmarked for later.
    Today, 06:37 PM
  • SEQadmin2
    16S/ITS Profiling vs. Shotgun Metagenomics: What Changes in Your Workflow?
    by SEQadmin2





    NGS has revolutionized clinical microbiology by offering new ways to identify and surveil pathogens without prior knowledge of the target. One popular approach is 16S rRNA and ITS amplicon sequencing, which targets specific microbial markers. But studies are increasingly incorporating shotgun metagenomics, a technique that sequences all DNA from a sample without targeting specific microbial markers. This approach allows more detailed taxonomic classification,...
    10-30-2026, 06:46 AM

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by SEQadmin2, 11-12-2026, 07:59 AM
0 responses
27 views
0 reactions
Last Post SEQadmin2  
Started by SEQadmin2, 11-05-2026, 09:25 AM
0 responses
35 views
0 reactions
Last Post SEQadmin2  
Started by SEQadmin2, 11-03-2026, 11:54 AM
2 responses
56 views
0 reactions
Last Post kinpict.com  
Started by SEQadmin2, 10-28-2026, 09:57 AM
0 responses
175 views
0 reactions
Last Post SEQadmin2  
Working...