7e87cadb469b4e0eb4fb7f973154357cfe7fc345 lrnassar Mon Sep 21 15:56:18 2026 -0700 QA fixes for the mei (Mobile Insertions) track collection. refs #37524 Fix two data bugs found during QA and rebuild the affected bigBeds. meiEul1dbToBed.py looked up samples and individuals by name, but euL1db joins on 1-based row numbers, so neither join ever matched and the individual count, tissues, clinical conditions and populations were empty on all 8,991 insertions while the contributing-samples table printed row numbers. Both loaders now key on the row number, the table prints the sample name, and the adjacent population filter is case-insensitive so it actually drops "unknown". meiHgsvc3CsvToBed.py took alt[1:] on every record, which dropped the first base of the element on the 96 GRCh38 and 111 T2T-CHM13 records where PALMER2 is the only caller and ALT carries no anchor base; it now prefers INFO SEQ, which always matches SVLEN. Correct seven statements on the description pages against their sources: the HGSVC3 single-caller split was attributed to PALMER rather than L1ME-AID, its orthogonal concordance was 90.8% rather than 92.5%, euL1db was credited with aligning the L1HS consensus when the paper says it was processed from our RepeatMasker track, DeepMEI's network was described as a classifier rather than a genotyper and given the wrong training set, euL1db listed two detection methods absent from the data, and HMEID contradicted itself on the MELT ASSESS cutoff. Also: the SweGen bigDataUrl now points at _swegen.bb so the restricted callset is kept off the download server; the container page no longer claims the whole collection is long-read, lists the two euL1db subtracks, scopes its display conventions to the subtracks they describe, and cites all six papers; dead and wrong track links are repointed and pinned to a db; $db replaces hardcoded hg38 in paths on pages that serve three assemblies; the euL1db labels no longer carry hg38 counts and a lift note that made no sense on hg19; all six subtracks gain a dataVersion; the euL1db filter ranges match the data; and five autoSql field descriptions match what the files contain. Document the gbdb symlinks and the QA changes in doc/hg38/mei.txt, correct the HMEID bedToBigBed type there, and add an hg19.txt pointer since hg19 carries the two euL1db subtracks. diff --git src/hg/makeDb/trackDb/human/meiSwegen.html src/hg/makeDb/trackDb/human/meiSwegen.html index 306855a3594..0c700f6d14c 100644 --- src/hg/makeDb/trackDb/human/meiSwegen.html +++ src/hg/makeDb/trackDb/human/meiSwegen.html @@ -1,36 +1,39 @@ <h2>Description</h2> <p> This track shows <b>mobile element insertions (MEIs)</b> identified by <a href="https://melt.igs.umaryland.edu/" target="_blank">MELT</a> on the <a href="https://swefreq.nbis.se/dataset/SweGen" target="_blank">SweGen</a> cohort of 1,000 Swedish whole-genome samples (Ameur et al. 2017). Each site is an insertion of an Alu, L1 (LINE-1), SVA or HERV-K mobile element relative to the reference. The SweGen short-variant frequency data for the same cohort is shown in the -<a href="hgTrackUi?g=swegen">SweGen variant frequencies</a> subtrack of +<a href="hgTrackUi?db=hg38&g=swefreq">SweGen variant frequencies</a> subtrack of the SNV Frequencies collection. </p> <table class="stdTbl"> -<tr><th>Class</th><th>MEIs</th></tr> -<tr><td>Alu</td><td>14,467</td></tr> -<tr><td>L1</td><td>2,429</td></tr> -<tr><td>SVA</td><td>1,131</td></tr> -<tr><td>HERVK</td><td>73</td></tr> -<tr><th>Total (GRCh37)</th><th>18,100</th></tr> -<tr><th>Total (after liftOver to hg38)</th><th>18,090</th></tr> +<tr> + <th>Class</th> + <th>Original GRCh37 callset</th> + <th>Shown in this track (hg38)</th> +</tr> +<tr><td>Alu</td><td>14,467</td><td>14,461</td></tr> +<tr><td>L1</td><td>2,429</td><td>2,426</td></tr> +<tr><td>SVA</td><td>1,131</td><td>1,131</td></tr> +<tr><td>HERVK</td><td>73</td><td>72</td></tr> +<tr><th>Total</th><th>18,100</th><th>18,090</th></tr> </table> <p> For each MEI, the track reports the mobile element class, the insertion length, the MELT subfamily call (e.g. AluYa5, L1Ta), the target-site duplication sequence, the MELT ASSESS quality score, nearby gene context if the insertion lies in or close to a gene, the allele count (MELT_AN; despite the name this is the number of allele observations, not the allele number), the alt-allele frequency, and the MELT FILTER status. </p> <h2>Display Conventions and Configuration</h2> <p> An insertion has zero length on the reference: it attaches between @@ -93,31 +96,32 @@ requires a brief approval). The VCF uses GRCh37 contigs without a "chr" prefix; the conversion adds the prefix, drops the VCF header, maps SVTYPE codes (<tt>ALU</tt>, <tt>LINE1</tt>, <tt>SVA</tt>, <tt>HERVK</tt>) to the element class names used here, copies INFO fields through to the BED, and writes a bed9+9 file with 1-bp anchor intervals. The hg19 BED was then lifted to hg38 with UCSC <tt>liftOver</tt> (<tt>-tab -bedPlus=9</tt>), which mapped 18,090 of 18,100 records; 10 records fell into hg38-deleted regions and were dropped. The lifted BED was sorted and converted to bigBed using the <a href="https://github.com/ucscGenomeBrowser/kent/blob/master/src/hg/makeDb/scripts/mei/meiSwegen.as" target="_blank">meiSwegen.as</a> schema. Conversion and lift steps are documented in the <a href="https://github.com/ucscGenomeBrowser/kent/blob/master/src/hg/makeDb/doc/hg38/mei.txt" target="_blank">makeDoc file</a>; the scripts live in <a href="https://github.com/ucscGenomeBrowser/kent/tree/master/src/hg/makeDb/scripts/mei" -target="_blank">src/hg/makeDb/scripts/mei</a>. +target="_blank">src/hg/makeDb/scripts/mei</a>, +and the track configuration is in <a href="https://github.com/ucscGenomeBrowser/kent/blob/master/src/hg/makeDb/trackDb/human/mei.ra" target="_blank">trackDb/human/mei.ra</a>. </p> <h3>Why the original GRCh37 MELT VCF rather than the GRCh38 SVDB files</h3> <p> The SweGen download portal also distributes a hg38 variant set (<tt>SweGen38_{ALU,L1,SVA,HERV}.vcf</tt>) for the same 1,001 samples, produced with SVDB after re-running on GRCh38. We chose to lift the original GRCh37 MELT VCF instead because the hg38 SVDB files contain 138,853 records (about 7.7× the MELT site count), and roughly 60% of those records are singletons (<tt>OCC=1</tt>) without any quality filter. They also drop most of the per-site annotation: no MELT subfamily call (e.g. AluYa5, L1Ta), no insertion length (<tt>SVLEN=0</tt> everywhere), no target-site duplication, no MELT ASSESS quality score, no gene context and no FILTER stratification (every site is marked <tt>PASS</tt>). The GRCh37 MELT VCF, lifted to @@ -135,31 +139,31 @@ download portal</a>. </p> <h2>Credits</h2> <p> Thanks to Adam Ameur, Diana Ekman, Jesper Eisfeldt, Daniel Nilsson and the SweGen consortium for generating and releasing the MELT MEI callset, and to SciLifeLab for producing the underlying SweGen WGS data. </p> <h2>References</h2> <p> -Ameur A, Dahlberg J, Olason P, Vezzi F, Karlsson R, Martin M, Viklund J, Kähäri AK, Lundin P, Che H +Ameur A, Dahlberg J, Olason P, Vezzi F, Karlsson R, Martin M, Viklund J, Kähäri AK, Lundin P, Che H <em>et al</em>. <a href="https://doi.org/10.1038/ejhg.2017.130" target="_blank"> SweGen: a whole-genome data resource of genetic variability in a cross-section of the Swedish population</a>. <em>Eur J Hum Genet</em>. 2017 Nov;25(11):1253-1260. PMID: <a href="https://www.ncbi.nlm.nih.gov/pubmed/28832569" target="_blank">28832569</a>; PMC: <a href="https://www.ncbi.nlm.nih.gov/pmc/articles/PMC5765326/" target="_blank">PMC5765326</a> </p> <p> Gardner EJ, Lam VK, Harris DN, Chuang NT, Scott EC, Pittard WS, Mills RE, 1000 Genomes Project Consortium, Devine SE. <a href="http://genome.cshlp.org/lookup/pmidlookup?view=long&pmid=28855259" target="_blank"> The Mobile Element Locator Tool (MELT): population-scale mobile element discovery and biology</a>. <em>Genome Res</em>. 2017 Nov;27(11):1916-1929.