0151d00a4a1d73a78c35f6158c6c936ff338faeb max Fri Apr 24 10:37:34 2026 -0700 NMD Escape: MANE subtrack, Rule 1 bug fix, transcript filter. refs #33737 - Add nmdEscMane subtrack (MANE Select Plus Clinical 1.5), built from /gbdb/hg38/mane/mane.bb. Reuses nmdEscTranscripts.html. - Fix Rule 1: measure 50 bp upstream of the transcript's last splice junction (including 3'UTR introns) rather than stripping 3'UTR from the exon list first. The old logic painted the entire last CDS exon as NMD-escape whenever the transcript had only one CDS exon, even when a 3'UTR intron sat far past the stop codon (e.g. NBDY: 207 bp of CDS over-painted for a junction 2.6 kb past the stop). - Add --rule1-mode {cds,mrna} (default cds): cds counts only CDS bp on the walk-back (paints up to 50 bp of CDS matching the rule label literally); mrna counts mRNA bp and clips to CDS (tracks the 55 bp rule literature). Documented in makeDoc. - Rule 4: when a 3'UTR intron exists, the last CDS-containing exon has a downstream EJC and is now eligible for the long-exon rule. - Mouseover lists contributing transcript accessions when 1-3 items collapse into a region; falls back to a count above that. - Add filterText/filterType/filterLabel on all three escape subtracks so a user can narrow the display to one transcript. - genePredNmdEsc: --gene-sym-field (default 17 for Gencode; pass 18 for MANE, whose HGNC symbol lives in bigGenePred geneName2). - Add findShortTxLongUtrIntron.py helper for finding MANE transcripts with long UTR introns (used to pick NMD edge-case test cases). Post-fix collapsed-region counts (--rule1-mode=cds): MANE 1.5: 67,752 Gencode V49: 233,375 RefSeq Curated: 112,356 diff --git src/hg/makeDb/trackDb/human/hg38/nmdEscTranscripts.html src/hg/makeDb/trackDb/human/hg38/nmdEscTranscripts.html index b9207fc0d6f..cc74d2318c6 100644 --- src/hg/makeDb/trackDb/human/hg38/nmdEscTranscripts.html +++ src/hg/makeDb/trackDb/human/hg38/nmdEscTranscripts.html @@ -1,37 +1,43 @@

Description

The NMD escape ruleset tracks show predicted regions where a premature termination codon (PTC) or frameshift variant is likely to cause the transcript to escape nonsense-mediated decay (NMD), leading to the production of an aberrant truncated protein rather than degradation of the mRNA.

The following rules were applied to transcript annotations to define predicted NMD escape regions (Nagy et al, Trends Biochem Sci 1998 and Lindeboom et al, Nat Genet 2016):

    -
  1. 50 bp rule: The entire last coding exon plus the last 50 bp of - the penultimate coding exon. A PTC here has no downstream exon-exon - junction (or is too close to the last one) for NMD to be triggered. - Non-protein-coding 3' exons are not counted when identifying the last - coding junction. Note: when the penultimate coding exon is shorter than - 50 bp, the annotated region extends only to the upstream junction of - that exon and does not walk further upstream. A small number of - transcripts with unusually short penultimate coding exons are affected.
    2. +
  2. 50 bp rule: Coding positions within 50 bp (mRNA distance) + upstream of the transcript's last splice junction, plus any coding + sequence downstream of that junction. A PTC in this window has no + downstream exon-exon junction (or is too close to the last one) for + NMD to be triggered. The last junction is determined from all exons + of the transcript, including 3′UTR introns, since those introns + deposit EJCs that can trigger NMD. For transcripts with no 3′UTR + intron (the common case), this reduces to the entire last coding exon + plus the last 50 bp of the penultimate coding exon. For transcripts + with a 3′UTR intron (~4.5% of MANE transcripts), the last + junction sits downstream of the stop codon; the escape region is only + the stretch of CDS within 50 bp (mRNA distance) of that junction, so + if the junction is more than 50 bp past the stop codon no CDS position + escapes via this rule.
  3. No downstream EJC rule: Transcripts with a single coding exon and no 3′UTR intron. No exon-exon junction exists downstream of the stop codon, so no EJC is deposited that could trigger NMD at a PTC. This covers truly intronless transcripts as well as transcripts whose only introns are in the 5′UTR (where EJCs are cleared by the scanning 40S ribosomal subunit or sit upstream of the stop and are never encountered by the terminating ribosome). Transcripts with a single coding exon but a 3′UTR intron are excluded, because that intron deposits an EJC downstream of the stop codon that can trigger NMD.
  4. Start-proximal region: The first 100 bp of coding nucleotides. PTCs in this region do not lead to NMD, a phenomenon known as start-proximal NMD insensitivity. One proposed mechanism, supported by experimental evidence, is re-initiation of translation at a downstream AUG codon.
  5. Long exon rule: Coding exons longer than 400 bp (excluding the last coding exon, which is already covered by the 50 bp rule). Lindeboom et al. @@ -94,33 +100,34 @@ downstream EJC reduces UPF1-EJC contact.

    6. Display Conventions and Configuration

    Regions from overlapping transcripts with the same coordinates are collapsed into a single item. The gene symbol is shown as the item name. Mouseover displays the NMD escape rule and the number of transcripts. The details page lists all contributing transcript IDs.

    Items are colored by the NMD escape rule that applies:

    Data Access

    The data underlying this track can be explored interactively with the