888e7470c14eeecdca310ed36bb45c3c00ae8052 lrnassar Tue Apr 21 15:14:04 2026 -0700 QA fixes for MPRA superTrack. refs #37359 Fix broken mpraVarDb bigDataUrl — pointed at /gbdb/hg38/mpra/mpravardb.bb but the file is at /gbdb/hg38/mpra/mpravardb/mpravardb.bb, causing hgTrackDb -strict to silently drop the subtrack. Rebuild mpravardb.bb after two fixes in mpravardbToBed.py: sanitize UTF-8 in user-visible string fields (curly quotes, primes, NBSP mojibake) that the browser does not transcode, eliminating ~246k non-ASCII occurrences across 42% of rows; and change safe_float / pval_to_score to write NaN and return score 0 for NA / out-of-range p-values instead of 0.0 and score 1000 (previously inflated untested variants to the top of score-sorted views). trackDb stanza cleanup: shorten mpraVarDb longLabel, drop superfluous type bed 4 from superTrack, make bigBed 9+13 explicit, remove redundant mouseOverField, align parent mpra on, add filterValues for cell_line/assay/cellLine and filterByRange sliders for percentile_rank / fdr / log2FC, add labelFields and maxWindowToDraw. Description pages: add cross-species disclosure (mouse reporter cells used to assay human sequences), update mpraVarDb header to post-liftOver count 239,028 with Studies-table footnote, fix mpraVarDb.html download-server paths, soften imprecise "51 MPRA experiments" claim in mpra.html and mprabase.html. relatedTracks.ra: reciprocal mpra <-> wgEncodeReg4 and mpra <-> cCREs. Expand mpra.txt makedoc with upstream provenance and QA-rebuild log. diff --git src/hg/makeDb/trackDb/human/hg38/mpraVarDb.html src/hg/makeDb/trackDb/human/hg38/mpraVarDb.html index 99e1158d586..7048cc83bf5 100644 --- src/hg/makeDb/trackDb/human/hg38/mpraVarDb.html +++ src/hg/makeDb/trackDb/human/hg38/mpraVarDb.html @@ -1,25 +1,31 @@
-The MPRAVarDB track shows 242,818 variants from 18 MPRA studies compiled -in the MPRAVarDB database +The MPRAVarDB track shows 239,028 variants successfully mapped to hg38 +(from 242,818 total) across 18 MPRA studies compiled in the MPRAVarDB database (Wang et al., 2024). Each variant was experimentally tested in an MPRA experiment to evaluate whether it affects transcriptional regulatory activity. The database covers over 30 cell lines and 30 human diseases and traits, including neurodegenerative diseases, immune disorders, melanoma, multiple myeloma, and autoimmune diseases.
++Note on cell lines: The cell line shown for each variant is the reporter +cell line in which the human regulatory element was assayed. Several studies +used mouse cell lines (e.g. Neuro-2a, N2A, NIH/3T3, MIN6) as reporter systems +for human sequences; these variants retain human (hg38) coordinates. +
Items are colored by statistical significance:
Each item shows the variant name (rsID when available, otherwise chr:pos:ref>alt), the reference and alternate alleles, the associated disease or trait, cell line, log2 fold change, p-value, and FDR.
@@ -153,62 +159,66 @@+Variant counts above are from the source publications (pre-liftOver totals). +Of 242,818 total source variants, 239,028 lifted successfully to hg38; see Methods. +
Data was downloaded from the
MPRAVarDB web server.
Variants originally mapped to hg19 (213,689 of 242,818) were lifted to hg38
using liftOver. 114 variants could not be mapped and were excluded.
The remaining variants were merged with the 29,129 natively hg38-mapped variants
to produce a total of 239,028 hg38 records.
The data can be explored interactively in table format with the Table Browser or the Data Integrator and exported from there to spreadsheet or tab-sep tables. From scripts, the data can be accessed through our API, track=mpraVarDb.
For automated download and analysis, the genome annotation is stored in a bigBed file that can be downloaded from -our download server. +our download server. The file for this track is called mpravardb.bb. Individual regions or the whole genome annotation can be obtained using our tool bigBedToBed, which can be compiled from the source code or downloaded as a precompiled binary for your system. Instructions for downloading source code and binaries can be found here. The tool can also be used to obtain features within a given range, e.g. -bigBedToBed http://hgdownload.soe.ucsc.edu/gbdb/hg38/mpra/mpravardb.bb -chrom=chr21 -start=0 -end=100000000 stdout +bigBedToBed http://hgdownload.soe.ucsc.edu/gbdb/hg38/mpra/mpravardb/mpravardb.bb -chrom=chr21 -start=0 -end=100000000 stdout
The original annotation source data can be downloaded from the MPRAVarDB web server.
Thanks to Tao Wang and colleagues at the University of Florida for creating and maintaining the MPRAVarDB database.
Wang T, Matreyek KA, Yang X.