9541aea4204c8079ea2e464403f2225f6aa33cfd
max
  Mon Sep 14 07:59:00 2026 -0700
uniprot: show the splice variant track, and filter the CAT alignments

Two things that were quietly missing.

unipSplice has been built since the pipeline rewrite in #19351 and never had a
trackDb stanza - "git log -S unipSplice" on uniprot.ra returns nothing, so it was
never wired up rather than deliberately dropped. It holds UniProt's splice variant
features, 28962 of them on hg38 and 29195 on hs1, and has been invisible on every
assembly for years. Added to uniprot.ra and to the archive/contrib template, which
is where the hub-served assemblies get their trackDb. The file exists on the same
131 assemblies as unipDomain.bb, so no stanza points at anything missing.

The alignments on a CAT assembly were not being filtered with pslSelect. That
filter maps a UniProt accession to the transcripts UniProt cross-references for it,
and the README is blunt about how much it matters for protein families with nearly
identical transcripts. It understood two kinds of id, Ensembl and RefSeq, and CAT
names a transcript after its source gene, so hs1 matched neither and fell through
unfiltered. Every row of the CAT bigBed carries the Ensembl transcript it was lifted
from, so catSourceTransMap joins on that and UniProt's Ensembl cross-reference does
the rest.

The mapping had to become one-to-many for this: one Ensembl transcript can name
several of ours, both because paralogs are lifted from the same source and because
duplicate CAT names were given -dup suffixes earlier. It reproduces those suffixes
by walking the bigBed in the same order rather than attaching every paralog to every
source, which measured 234903 of 234903 hs1 transcripts mapped, no duplicates,
against 330749 for the loose version. For every other gene track the lists hold one
element and the behaviour is unchanged, which is checked: a single-valued entry
still writes exactly one pair line and still counts a version difference, a
multi-valued one writes a line per transcript, and an id we do not have is still
skipped so pslSelect -qPass passes it through.

refs #38300

diff --git src/hg/makeDb/trackDb/uniprot.ra src/hg/makeDb/trackDb/uniprot.ra
index 1d03dd1e992..02a8e9aef8c 100644
--- src/hg/makeDb/trackDb/uniprot.ra
+++ src/hg/makeDb/trackDb/uniprot.ra
@@ -221,22 +221,34 @@
         filterValues.status Manually reviewed (Swiss-Prot),Unreviewed (TrEMBL)
         mouseOver <b>UniProt record</b>: $uniProtId<br> <b>Position</b>: $position<br> <b>UniProt status</b>: $status
 
         track unipConflict
         shortLabel Seq. Conflicts
         longLabel UniProt Sequence Conflicts
         priority 13
         parent uniprot off
         visibility dense
         type bigBed 12 +
         bigDataUrl /gbdb/$D/uniprot/unipConflict.bb
         urls uniProtId="http://www.uniprot.org/uniprot/$$#Sequence_conflict_section" pmids="https://www.ncbi.nlm.nih.gov/pubmed/$$"
         filterValues.status Manually reviewed (Swiss-Prot),Unreviewed (TrEMBL)
         mouseOver <b>UniProt record</b>: $uniProtId<br> <b>Position</b>: $position<br> <b>UniProt status</b>: $status
 
+        track unipSplice
+        shortLabel Splice Variants
+        longLabel UniProt Splice Variants
+        priority 14
+        parent uniprot
+        visibility dense
+        type bigBed 12 +
+        bigDataUrl /gbdb/$D/uniprot/unipSplice.bb
+        urls uniProtId="http://www.uniprot.org/uniprot/$$#sequences" pmids="https://www.ncbi.nlm.nih.gov/pubmed/$$"
+        filterValues.status Manually reviewed (Swiss-Prot),Unreviewed (TrEMBL)
+        mouseOver <b>UniProt record</b>: $uniProtId<br> <b>Position</b>: $position<br> <b>UniProt status</b>: $status
+
 searchTable unipAliSwissprot
 searchType bigBed
 searchDescription UniProt/SwissProt manually curated Proteins mapped to the genome by UCSC
 
 searchTable unipAliTrembl
 searchType bigBed
 searchDescription UniProt/TrEMBL predicted Proteins mapped to the genome by UCSC