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:
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!
Here is my best guess workflow:
- 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.)
- 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.
- 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.)
- 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.
- (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!
Comment