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 @@ <h2>Description</h2> <p> -The <b>MPRAVarDB</b> track shows 242,818 variants from 18 MPRA studies compiled -in the MPRAVarDB database +The <b>MPRAVarDB</b> track shows 239,028 variants successfully mapped to hg38 +(from 242,818 total) across 18 MPRA studies compiled in the MPRAVarDB database (<a href="https://pubmed.ncbi.nlm.nih.gov/38617248/" target="_blank">Wang et al., 2024</a>). 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. </p> +<p> +<b>Note on cell lines:</b> 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. +</p> <h2>Display Conventions</h2> <p> Items are colored by statistical significance: <ul> <li><b><span style="color: #C80000;">Dark red</span></b>: FDR < 0.05 (significant after multiple testing correction) — 22,465 variants (9.3%)</li> <li><b><span style="color: #FFA500;">Orange</span></b>: nominal p-value < 0.05 but FDR ≥ 0.05 — 17,780 variants (7.3%)</li> <li><b><span style="color: #BEBEBE;">Grey</span></b>: not significant (p-value ≥ 0.05) — 202,573 variants (83.4%)</li> </ul> </p> <p> 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. </p> @@ -153,62 +159,66 @@ <tr> <td><a href="https://pubmed.ncbi.nlm.nih.gov/31477794/" target="_blank">Rao et al., 2021</a></td> <td>88</td> <td>Alcohol use disorder</td> <td>BLA, CE, NAC, SFC</td> <td>SNPs in 3'UTR of 88 genes from allele-specific expression analysis (30 AUD subjects vs 30 controls)</td> </tr> <tr> <td><a href="https://pubmed.ncbi.nlm.nih.gov/27259154/" target="_blank">Ulirsch et al., 2016</a></td> <td>62</td> <td>Red blood cell traits</td> <td>K562, K562+GATA1</td> <td>2,756 variants in strong LD with 75 sentinel variants associated with RBC traits</td> </tr> </table> +<p> +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. +</p> <h2>Methods</h2> <p> Data was downloaded from the <a href="https://mpravardb.rc.ufl.edu/" target="_blank">MPRAVarDB web server</a>. Variants originally mapped to hg19 (213,689 of 242,818) were lifted to hg38 using <code>liftOver</code>. 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. </p> <h2>Data Access</h2> <p> The data can be explored interactively in table format with the <a href="../cgi-bin/hgTables">Table Browser</a> or the <a href="../cgi-bin/hgIntegrator">Data Integrator</a> and exported from there to spreadsheet or tab-sep tables. From scripts, the data can be accessed through our <a href="https://api.genome.ucsc.edu" target="_blank">API</a>, track=<i>mpraVarDb</i>. </p> <p> For automated download and analysis, the genome annotation is stored in a bigBed file that can be downloaded from -<a href="http://hgdownload.soe.ucsc.edu/gbdb/hg38/mpra" target="_blank">our download server</a>. +<a href="http://hgdownload.soe.ucsc.edu/gbdb/hg38/mpra/mpravardb" target="_blank">our download server</a>. The file for this track is called <tt>mpravardb.bb</tt>. Individual regions or the whole genome annotation can be obtained using our tool <tt>bigBedToBed</tt>, 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 <a href="http://hgdownload.soe.ucsc.edu/downloads.html#utilities_downloads" target="_blank">here</a>. The tool can also be used to obtain features within a given range, e.g. -<tt>bigBedToBed http://hgdownload.soe.ucsc.edu/gbdb/hg38/mpra/mpravardb.bb -chrom=chr21 -start=0 -end=100000000 stdout</tt> +<tt>bigBedToBed http://hgdownload.soe.ucsc.edu/gbdb/hg38/mpra/mpravardb/mpravardb.bb -chrom=chr21 -start=0 -end=100000000 stdout</tt> </p> <p> The original annotation source data can be downloaded from the <a href="https://mpravardb.rc.ufl.edu/" target="_blank">MPRAVarDB web server</a>. </p> <h2>Credits</h2> <p> Thanks to Tao Wang and colleagues at the University of Florida for creating and maintaining the MPRAVarDB database. </p> <h2>References</h2> <p> Wang T, Matreyek KA, Yang X.