Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • DrOM
    replied
    Hallo!

    have you advanced with the genomic coordinates for TargetScan6.2? I think you can skip step #2, since the target file also provides the position in the ungapped UTRs of that species. Therefore it should be ok to strip the gaps of the alignment file, get the genomic coordinates for the UTRs and proceed like you outlined above.

    It would be great if you could share your results!!
    Cheers

    Pablo

    Leave a comment:


  • jmw86069
    started a topic TargetScan v6.1 genome coordinates

    TargetScan v6.1 genome coordinates

    I'm trying to find or create genome coordinates for TargetScan v6.1, released March 2012. UCSC provides a track from v5.1, but nothing newer. This thread (link UCSC mailing list) suggests it's trickier than one would think. So before embarking on that quest, I wanted to poll this board to see if anyone already has created genome coordinates from TargetScan v6.1?

    Here is my best guess workflow:
    1. Download UTR_Sequences.txt from TargetScan for the given species, grep out only the taxon ID needed (10090 for mouse from the TargetScanMouse data; 9606 for human from the TargetScanHuman data; etc.)
    2. Create gapped-to-ungapped mapping of the UTR_Sequences.txt data (which is in gapped alignment format) so I can convert the reported TargetScan gapped coordinates to un-gapped UTR coordinates.
    3. Strip out the gaps from each sequence to create an un-gapped UTR FASTA file, then BLAT them to the appropriate genome (mm9, hg18/hg19, etc.), in order to convert UTR ungapped coordinates to genome coordinates. This step seems necessary because the TargetScan UTRs were created using some logic, and aren't directly from any available UCSC track (that I can tell anyway.)
    4. Use the Conserved_Family_Info.txt file for the given organism, which provides gapped UTR coordinates for each predicted miRNA site. Convert each coordinate using the logic above.
    5. (gapped UTR coordinate) --> (ungapped UTR coordinate) --> (genome coordinate)


    I don't know how much I lose by using BLAT to re-align UTRs back to the genome. If anyone has a suggested alternative (e.g. MUMmer?), I'd greatly welcome it! I'm really hoping for well-behaved alignments!

Latest Articles

Collapse

  • seqadmin
    Exploring the Dynamics of the Tumor Microenvironment
    by seqadmin




    The complexity of cancer is clearly demonstrated in the diverse ecosystem of the tumor microenvironment (TME). The TME is made up of numerous cell types and its development begins with the changes that happen during oncogenesis. “Genomic mutations, copy number changes, epigenetic alterations, and alternative gene expression occur to varying degrees within the affected tumor cells,” explained Andrea O’Hara, Ph.D., Strategic Technical Specialist at Azenta. “As...
    07-08-2024, 03:19 PM

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by seqadmin, Yesterday, 06:46 AM
0 responses
9 views
0 likes
Last Post seqadmin  
Started by seqadmin, 07-24-2024, 11:09 AM
0 responses
26 views
0 likes
Last Post seqadmin  
Started by seqadmin, 07-19-2024, 07:20 AM
0 responses
160 views
0 likes
Last Post seqadmin  
Started by seqadmin, 07-16-2024, 05:49 AM
0 responses
127 views
0 likes
Last Post seqadmin  
Working...
X