@sjcrh/proteinpaint-client 2.211.0 → 2.211.1-0

package/dist/chunk-7J45QGK7.js

@@ -0,0 +1,217 @@
+import {
+  DATermTypes,
+  DMRCATE_DEFAULTS
+} from "./chunk-UODENFH4.js";
+import {
+  dofetch3
+} from "./chunk-RMUK3TLD.js";
+import {
+  DMR_SCAN_ELEMENT_TYPE
+} from "./chunk-57Z4VYLM.js";
+import {
+  rgb
+} from "./chunk-Q5RDQNIT.js";
+
+// plots/volcano/colors.ts
+function getGroupColors(config) {
+  const groups = config?.samplelst?.groups;
+  const termValues = config?.tw?.term?.values;
+  const rawDown = termValues?.[groups?.[0]?.name]?.color || "red";
+  const rawUp = termValues?.[groups?.[1]?.name]?.color || "blue";
+  return {
+    controlColor: toHex(rawDown, "red"),
+    caseColor: toHex(rawUp, "blue")
+  };
+}
+function toHex(color, fallback) {
+  const c = rgb(color || fallback);
+  return c.displayable() ? c.formatHex() : rgb(fallback).formatHex();
+}
+
+// plots/volcano/model/VolcanoModel.ts
+var VolcanoModel = class {
+  /** TODO: This model is used in both the volcano and gsea.
+   * In the future, create base model in DA and use specific
+   * classes for the volcano and gsea. */
+  constructor(plot, termType) {
+    this.plot = plot;
+    this.app = plot.app;
+    this.termType = termType;
+  }
+  /** May use mapper instead as more termTypes are added */
+  async getData(config, settings) {
+    this.config = config;
+    this.settings = settings;
+    if (this.termType === DATermTypes.GENE_EXPRESSION) {
+      const body = await this.getGERequestBody();
+      const response = await dofetch3("termdb/DE", { body, signal: this.plot.api?.getAbortSignal() });
+      if (response && !response.error) response.daRequest = body;
+      return response;
+    }
+    if (this.termType === DATermTypes.DNA_METHYLATION) {
+      const body = await this.getDMRequestBody();
+      const response = await dofetch3("termdb/diffMeth", { body, signal: this.plot.api?.getAbortSignal() });
+      if (response && !response.error) response.daRequest = body;
+      return response;
+    }
+    if (this.termType === DATermTypes.SINGLECELL_CELLTYPE) {
+      const body = await this.getSCCTRequestBody();
+      return await dofetch3("termdb/singlecellDEgenes", { body, signal: this.plot.api?.getAbortSignal() });
+    }
+    if (this.termType === DATermTypes.PROTEOME_DAP) {
+      const body = this.getDapRequestBody();
+      return await dofetch3("termdb/dapVolcano", { body, signal: this.plot.api?.getAbortSignal() });
+    }
+    if (this.termType === DATermTypes.SINGLECELL_GENE_EXPRESSION) {
+    }
+    throw new Error(`Volcano plot does not support route for termType='${this.termType}'`);
+  }
+  //Gene expression
+  async getGERequestBody() {
+    await this.getOtherSamples(this.config.samplelst);
+    const state = this.app.getState();
+    const body = {
+      kind: "DE",
+      genome: this.app.vocabApi.vocab.genome,
+      dslabel: this.app.vocabApi.vocab.dslabel,
+      method: this.settings.method,
+      min_count: this.settings.minCount,
+      min_total_count: this.settings.minTotalCount,
+      samplelst: this.config.samplelst,
+      filter: state.termfilter.filter,
+      filter0: state.termfilter.filter0,
+      cpm_cutoff: this.settings.cpmCutoff,
+      volcanoRender: this.getVolcanoRender()
+    };
+    const pseudobulk = this.config.tw?.pseudobulk;
+    if (pseudobulk) body.pseudobulk = pseudobulk;
+    this.addConfounderTw(body);
+    return body;
+  }
+  //DNA methylation
+  async getDMRequestBody() {
+    await this.getOtherSamples(this.config.samplelst);
+    const state = this.app.getState();
+    const body = {
+      kind: "DM",
+      genome: this.app.vocabApi.vocab.genome,
+      dslabel: this.app.vocabApi.vocab.dslabel,
+      samplelst: this.config.samplelst,
+      filter: state.termfilter.filter,
+      filter0: state.termfilter.filter0,
+      min_samples_per_group: this.settings.minSamplesPerGroup,
+      exclude_sex_chr: this.settings.excludeSexChr,
+      /* Omitted rather than sent as 'promoter' when it is the default, so a request
+      from a promoter-only dataset is byte-identical to what this client sent before
+      the element picker existed. The server resolves an absent element_type to
+      'promoter'. This does NOT preserve cache keys -- the key object gained the
+      field server-side, so every pre-existing dm/ entry is orphaned on deploy
+      regardless of what the client sends. */
+      ...this.settings.elementType && this.settings.elementType != "promoter" ? { element_type: this.settings.elementType } : {},
+      // the scan's own knobs; the server ignores them for any other element type
+      ...this.settings.elementType == DMR_SCAN_ELEMENT_TYPE ? {
+        scan: {
+          ...this.settings.scanChromosome ? { chromosome: this.settings.scanChromosome } : {},
+          backgroundCorrection: !!this.settings.backgroundCorrection,
+          minCpgs: this.settings.minCpgs,
+          profileBinBp: this.settings.profileBinBp,
+          // only when changed: an explicit default would orphan every cached scan
+          ...this.settings.lambda != DMRCATE_DEFAULTS.lambda ? { lambda: this.settings.lambda } : {},
+          ...this.settings.C != DMRCATE_DEFAULTS.C ? { C: this.settings.C } : {},
+          ...this.settings.fdrCutoff != DMRCATE_DEFAULTS.fdrCutoff ? { fdrCutoff: this.settings.fdrCutoff } : {}
+        }
+      } : {},
+      volcanoRender: this.getVolcanoRender()
+    };
+    if (this.settings.elementType != DMR_SCAN_ELEMENT_TYPE) this.addConfounderTw(body);
+    return body;
+  }
+  /** Parameters telling the server to run the `volcano` Rust renderer and return a
+   * volcano PNG + top-significant rows instead of the full dot list. */
+  getVolcanoRender() {
+    const dotRadius = Math.max(this.settings.width, this.settings.height) / 80;
+    const { caseColor, controlColor } = getGroupColors(this.config);
+    const useDeltaBeta = this.termType === DATermTypes.DNA_METHYLATION && this.settings.xAxis == "delta_beta";
+    return {
+      significanceThresholds: {
+        pValueCutoff: this.settings.pValue,
+        pValueType: this.settings.pValueType,
+        foldChangeCutoff: useDeltaBeta ? this.settings.deltaBetaCutoff : this.settings.foldChangeCutoff
+      },
+      ...useDeltaBeta ? { xField: "delta_beta" } : {},
+      // Tied to the delta-beta axis: the control that sets it is only offered for
+      // methylation, and centering a log2 fold-change axis is a different conversation.
+      ...useDeltaBeta && this.settings.centerDeltaBeta ? { centerX: true } : {},
+      pixelWidth: this.settings.width,
+      pixelHeight: this.settings.height,
+      colorSignificant: toHex(this.settings.defaultSignColor, "red"),
+      colorSignificantUp: caseColor,
+      colorSignificantDown: controlColor,
+      colorNonsignificant: toHex(this.settings.defaultNonSignColor, "black"),
+      dotRadius,
+      maxInteractiveDots: this.settings.maxInteractiveDots,
+      // Render the PNG at device-pixel resolution so it stays sharp on
+      // retina screens. The server reports the plot extent in CSS-space,
+      // so SVG overlay coords are unaffected.
+      //
+      // Oversample by 2× so the PNG also stays sharp when the user
+      // *zooms in after* the initial render (the captured DPR is frozen
+      // at fetch time — bigger headroom = more tolerable post-render
+      // zoom before pixelation appears). The server clamp keeps the
+      // bitmap memory bounded.
+      devicePixelRatio: (typeof window !== "undefined" ? window.devicePixelRatio : 1) * 2
+    };
+  }
+  //This is a workaround until the server can accept an arr of confounder tws
+  addConfounderTw(body) {
+    const confounders = this.config?.confounderTws;
+    if (confounders?.length) {
+      body.tw = this.config.confounderTws[0];
+      if (confounders.length > 1) body.tw2 = this.config.confounderTws[1];
+    }
+  }
+  //Single cell cell type
+  getSCCTRequestBody() {
+    const body = {
+      genome: this.app.vocabApi.vocab.genome,
+      dslabel: this.app.vocabApi.vocab.dslabel,
+      sample: this.config.sample,
+      termId: this.config.termId,
+      categoryName: this.config.categoryName,
+      volcanoRender: this.getVolcanoRender()
+    };
+    return body;
+  }
+  getDapRequestBody() {
+    const { organism, assay, cohort } = this.config.proteomeDetails;
+    return {
+      genome: this.app.vocabApi.vocab.genome,
+      dslabel: this.app.vocabApi.vocab.dslabel,
+      organism,
+      assay,
+      cohort,
+      volcanoRender: this.getVolcanoRender()
+    };
+  }
+  /** retrieve the sampleId/sampleName for samples in
+   * the "others" group instead of using {in: false} */
+  async getOtherSamples(samplelst) {
+    const othersSamplesGroup = samplelst.groups.find((g) => !g.in);
+    if (!othersSamplesGroup) return;
+    const state = this.app.getState();
+    const samplesGroup = samplelst.groups.find((g) => g.in);
+    othersSamplesGroup.values = [];
+    for (const s of await this.plot.vocabApi.getFilteredSampleList(state.termfilter.filter)) {
+      if (!samplesGroup.values.some((i) => i.sampleId == s.id)) {
+        othersSamplesGroup.values.push({ sampleId: s.id, sample: s.name });
+      }
+    }
+    othersSamplesGroup.in = true;
+  }
+};
+
+export {
+  getGroupColors,
+  VolcanoModel
+};
+//# sourceMappingURL=chunk-7J45QGK7.js.map

package/dist/chunk-7L3N6RGH.js

@@ -0,0 +1,134 @@
+import {
+  addGeneSearchbox,
+  isoformSelect,
+  pickCollectionFraction,
+  sayerror
+} from "./chunk-UODENFH4.js";
+import {
+  Menu
+} from "./chunk-7XZA2XR2.js";
+import {
+  dofetch3
+} from "./chunk-RMUK3TLD.js";
+import {
+  ISOFORM_EXPRESSION,
+  getColors
+} from "./chunk-57Z4VYLM.js";
+
+// termdb/handlers/isoformExpression.ts
+var SearchHandler = class {
+  constructor() {
+    this.currentGene = null;
+  }
+  init(opts) {
+    this.callback = opts.callback;
+    this.app = opts.app;
+    this.termCollectionSelectionMode = opts.termCollectionSelectionMode;
+    const holder = opts.holder.append("div").style("padding", "10px 0px");
+    this.dom = {
+      errDiv: holder.append("div").style("margin", "5px 0px").style("display", "none")
+    };
+    const geneSearch = addGeneSearchbox({
+      tip: new Menu({ padding: "0px" }),
+      genome: opts.genomeObj,
+      row: holder,
+      searchOnly: "gene",
+      callback: async () => {
+        try {
+          this.dom.errDiv.style("display", "none");
+          if (!geneSearch.geneSymbol) throw new Error("No gene selected");
+          if (geneSearch.geneSymbol === this.currentGene) return;
+          this.currentGene = geneSearch.geneSymbol;
+          if (this.dom.isoformDiv) this.dom.isoformDiv.remove();
+          this.dom.isoformDiv = holder.append("div");
+          await this.showIsoforms(geneSearch.geneSymbol, opts.genomeObj);
+        } catch (e) {
+          this.dom.errDiv.style("display", "block");
+          sayerror(this.dom.errDiv, "Error: " + (e instanceof Error ? e.message : String(e)));
+        }
+      }
+    });
+  }
+  async showIsoforms(gene, genomeObj) {
+    if (!gene) throw new Error("No gene selected");
+    const data = await dofetch3("genelookup", { body: { genome: genomeObj.name, input: gene, deep: 1 } });
+    if (!data.gmlst?.length) throw new Error(`No isoforms found for ${gene}`);
+    const enstCandidates = data.gmlst.filter((gm) => gm.isoform?.startsWith("ENST"));
+    if (enstCandidates.length === 0) throw new Error(`No Ensembl transcript isoforms found for ${gene}`);
+    const { available } = await dofetch3("termdb/isoformAvailability", {
+      body: {
+        genome: genomeObj.name,
+        dslabel: this.app.vocabApi.vocab.dslabel,
+        isoforms: enstCandidates.map((gm) => gm.isoform)
+      }
+    });
+    const availableSet = new Set(available || []);
+    const enstModels = enstCandidates.filter((gm) => availableSet.has(gm.isoform));
+    if (enstModels.length === 0) throw new Error(`No isoforms with data found for ${gene}`);
+    if (gene !== this.currentGene) return;
+    const div = this.dom.isoformDiv;
+    div.append("div").style("margin-bottom", "8px").style("opacity", 0.65).text(`${gene} \u2014 select isoform(s):`);
+    isoformSelect({
+      holder: div,
+      allgm: enstModels,
+      multiSelect: true,
+      // a single checked isoform yields an individual term, 2+ yield a collection
+      getSubmitLabel: (selectedCount) => selectedCount === 1 ? "Select One Isoform" : `Create Collection (${selectedCount})`,
+      onMultiSelect: (selected) => {
+        if (selected.length === 1) {
+          this.selectIsoform(selected[0].isoform, gene);
+        } else {
+          this.selectCollection(selected, gene);
+        }
+      }
+    });
+  }
+  getUnit() {
+    return this.app.vocabApi.termdbConfig.queries.isoformExpression?.unit || "TPM";
+  }
+  selectIsoform(isoform, gene) {
+    const name = `${isoform} ${this.getUnit()}`;
+    this.callback({ isoform, gene, name, type: ISOFORM_EXPRESSION });
+  }
+  selectCollection(gms, gene) {
+    const unit = this.getUnit();
+    const termlst = gms.map((gm) => ({
+      id: gm.isoform,
+      name: gm.isoform,
+      type: ISOFORM_EXPRESSION,
+      isoform: gm.isoform
+    }));
+    const colorScale = getColors(termlst.length);
+    const term = {
+      type: "termCollection",
+      isCustom: true,
+      memberType: "numeric",
+      name: `${gene} Isoforms (${unit})`,
+      termlst,
+      propsByTermId: Object.fromEntries(termlst.map((term2) => [term2.id, { color: colorScale(term2.id) }])),
+      isleaf: true
+    };
+    if (this.termCollectionSelectionMode === "fraction") {
+      if (!this.dom?.isoformDiv) throw new Error("isoform result holder is missing");
+      this.dom.fractionDiv?.remove();
+      this.dom.fractionDiv = this.dom.isoformDiv.append("div").style("margin-top", "10px");
+      pickCollectionFraction({
+        holder: this.dom.fractionDiv,
+        term,
+        callback: (tw) => this.callback(tw)
+      });
+      return;
+    }
+    this.callback(term);
+  }
+};
+function filterIsoforms(gmlst, availableItems) {
+  const itemSet = new Set(availableItems);
+  return gmlst.filter((gm) => gm.isoform?.startsWith("ENST") && (itemSet.size === 0 || itemSet.has(gm.isoform)));
+}
+
+export {
+  SearchHandler,
+  filterIsoforms
+};
+//# sourceMappingURL=chunk-7L3N6RGH.js.map