biopipen 0.21.1__tar.gz → 0.22.0__tar.gz

{biopipen-0.21.1 → biopipen-0.22.0}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: biopipen
-Version: 0.21.1
+Version: 0.22.0
 Summary: Bioinformatics processes/pipelines that can be run from `pipen run`
 License: MIT
 Author: pwwang

biopipen-0.22.0/biopipen/__init__.py

@@ -0,0 +1 @@
+__version__ = "0.22.0"

{biopipen-0.21.1 → biopipen-0.22.0}/biopipen/core/filters.py

@@ -7,6 +7,7 @@ from typing import Any, List, Mapping
 
 from argx import Namespace
 from liquid.filters.manager import FilterManager
+from pipen_report.filters import register_component, render_ui, _tag
 
 filtermanager = FilterManager()
 
@@ -206,3 +207,144 @@ def r(
         return r(vars(obj), ignoreintkey, todot, sortkeys, skip, _i)
 
     return repr(obj)
+
+
+@register_component("fgsea")
+def _render_fgsea(
+    cont: Mapping[str, Any],
+    job: Mapping[str, Any],
+    level: int,
+    na_arg: str = "10",
+) -> str:
+    """Render fgsea report"""
+    # cont["dir"] is required
+    n_pathways = int(na_arg)
+    pathways = []
+    with Path(cont["dir"]).joinpath("fgsea.txt").open() as f:
+        next(f)  # skip header
+        for line in f:
+            pathway, _ = line.split("\t", 1)
+            pathways.append(pathway)
+            if len(pathways) >= n_pathways:
+                break
+
+    components = [
+        # Summary
+        {
+            "title": "Enrichment Analysis Summary",
+            "ui": "tabs",
+            "contents": [
+                {
+                    "title": "Plot",
+                    "ui": "flat",
+                    "contents": [
+                        {
+                            "kind": "image",
+                            "src": str(Path(cont["dir"]).joinpath("gsea_table.png")),
+                        }
+                    ],
+                },
+                {
+                    "title": "Table",
+                    "ui": "flat",
+                    "contents": [
+                        {
+                            "kind": "table",
+                            "src": str(Path(cont["dir"]).joinpath("fgsea.txt")),
+                        }
+                    ],
+                },
+            ]
+        },
+        # Pathways
+        {
+            "title": f"Enriched Pathways (Top {n_pathways})",
+            "ui": "table_of_images",
+            "contents": [
+                {
+                    "src": str(Path(cont["dir"]) / f"fgsea_{pw.replace('/', '-')}.png"),
+                    "title": pw,
+                }
+                for pw in pathways
+            ]
+        },
+    ]
+
+    return render_ui(components, "accordion", job, level)
+
+
+@register_component("pdf")
+def _render_pdf(
+    cont: Mapping[str, Any],
+    job: Mapping[str, Any],
+    level: int,
+) -> str:
+    """Render pdf report"""
+    # cont["src"] is required
+    height = cont.get("height", "600")
+    return _tag(
+        "embed",
+        src=str(cont["src"]),
+        type="application/pdf",
+        width="100%",
+        height=height,
+    )
+
+
+@register_component("gsea")
+def _render_gsea(
+    cont: Mapping[str, Any],
+    job: Mapping[str, Any],
+    level: int,
+) -> str:
+    """Render gsea report"""
+    # cont["dir"] is required
+    raise NotImplementedError()
+
+
+@register_component("enrichr")
+def _render_enrichr(
+    cont: Mapping[str, Any],
+    job: Mapping[str, Any],
+    level: int,
+) -> str:
+    """Render enrichr report"""
+    # cont["dir"] is required
+    dbs = [sumfile.stem[8:] for sumfile in Path(cont["dir"]).glob("Enrichr-*.txt")]
+    components = []
+
+    for db in dbs:
+        components.append(
+            {
+                "title": db,
+                "ui": "tabs",
+                "contents": [
+                    {
+                        "title": "Plot",
+                        "ui": "flat",
+                        "contents": [
+                            {
+                                "kind": "image",
+                                "src": str(
+                                    Path(cont["dir"]).joinpath(f"Enrichr-{db}.png")
+                                ),
+                            }
+                        ],
+                    },
+                    {
+                        "title": "Table",
+                        "ui": "flat",
+                        "contents": [
+                            {
+                                "kind": "table",
+                                "src": str(
+                                    Path(cont["dir"]).joinpath(f"Enrichr-{db}.txt")
+                                ),
+                            }
+                        ],
+                    },
+                ],
+            }
+        )
+
+    return render_ui(components, "accordion", job, level)

{biopipen-0.21.1 → biopipen-0.22.0}/biopipen/ns/scrna.py

@@ -257,6 +257,16 @@ class SeuratClustering(Proc):
             `object` is specified internally, and `-` in the key will be replaced with `.`.
             - resolution (type=float): The resolution of the clustering
             - <more>: See <https://satijalab.org/seurat/reference/findclusters>
+        cache (type=auto): Whether to cache the seurat object with cluster information.
+            If `True`, the seurat object will be cached in the job output directory, which will be not cleaned up when job is rerunning.
+            The cached seurat object will be saved as `<signature>.cached.RDS` file, where `<signature>` is the signature determined by
+            the input and envs of the process.
+            See <https://github.com/satijalab/seurat/issues/7849>, <https://github.com/satijalab/seurat/issues/5358> and
+            <https://github.com/satijalab/seurat/issues/6748> for more details.
+            To not use the cached seurat object, you can either set `cache` to `False` or delete the cached file at
+            `.pipen/<Pipeline>/SeuratClustering/0/output/<signature>.cached.RDS`.
+            You can also specify the directory to save the cached seurat object by setting `cache` to the directory path.
+
 
     Requires:
         r-seurat:
@@ -286,6 +296,7 @@ class SeuratClustering(Proc):
         "RunUMAP": {"reduction": "pca", "dims": 30},
         "FindNeighbors": {},
         "FindClusters": {"resolution": 0.8},
+        "cache": True,
     }
     script = "file://../scripts/scrna/SeuratClustering.R"
 
@@ -418,6 +429,7 @@ class SeuratClusterStats(Proc):
             - group-by: Same as `ident`. How the points are colored.
             - split-by: The column name in metadata to split the cells into different plots.
             - shape-by: The column name in metadata to use as the shape.
+            - subset: An expression to subset the cells, will be passed to `tidyrseurat::filter()`.
             - devpars (ns): The device parameters for the plots.
                 - res (type=int): The resolution of the plots.
                 - height (type=int): The height of the plots.
@@ -482,6 +494,7 @@ class SeuratClusterStats(Proc):
             "group-by": None,
             "split-by": None,
             "shape-by": None,
+            "subset": None,
             "reduction": "dim",
             "devpars": {"res": 100, "height": 800, "width": 1000},
         },
@@ -642,6 +655,7 @@ class CellsDistribution(Proc):
             Ignored if `cells_order` is specified.
         subset: An expression to subset the cells, will be passed to `dplyr::filter()` on metadata.
             This will be applied prior to `each`.
+        descr: The description of the case, will be shown in the report.
         devpars (ns): The device parameters for the plots.
             - res (type=int): The resolution of the plots
             - height (type=int): The height of the plots
@@ -676,6 +690,7 @@ class CellsDistribution(Proc):
         "cells_orderby": None,
         "cells_n": 10,
         "subset": None,
+        "descr": None,
         "devpars": {},
         "each": None,
         "section": "DEFAULT",
@@ -1669,7 +1684,7 @@ class RadarPlots(Proc):
         "cluster_order": [],
         "breaks": [],
         "direction": "intra-cluster",
-        "section": None,
+        "section": "DEFAULT",
         "devpars": {
             "res": 100,
             "width": 1200,
@@ -1719,6 +1734,8 @@ class MetaMarkers(Proc):
         dbs (list): The dbs to do enrichment analysis for significant
             markers See below for all libraries.
             <https://maayanlab.cloud/Enrichr/#libraries>
+        subset: The subset of the cells to do the analysis.
+            An expression passed to `dplyr::filter()`.
         p_adjust (choice): The method to adjust the p values, which can be used to filter the significant markers.
             See also <https://rdrr.io/r/stats/p.adjust.html>
             - holm: Holm-Bonferroni method
@@ -1759,6 +1776,7 @@ class MetaMarkers(Proc):
         "group-by": None,
         "idents": None,
         "each": None,
+        "subset": None,
         "prefix_each": True,
         "p_adjust": "BH",
         "dbs": [

{biopipen-0.21.1 → biopipen-0.22.0}/biopipen/ns/tcr.py

@@ -180,6 +180,7 @@ class Immunarch(Proc):
     - The dynamics of repertoires across time points/samples, provided by [`immunarch::trackClonotypes`](https://immunarch.com/reference/trackClonotypes.html)
     - The spectratype of clonotypes, provided by [`immunarch::spectratype`](https://immunarch.com/reference/spectratype.html)
     - The distributions of kmers and sequence profiles, provided by [`immunarch::getKmers`](https://immunarch.com/reference/getKmers.html)
+    - The V-J junction circos plots, implemented within the script of this process.
 
     Environment Variable Design:
         With different sets of arguments, a single function of the above can perform different tasks.
@@ -657,6 +658,24 @@ class Immunarch(Proc):
                 The values will be passed to the corresponding arguments above.
                 If any of these arguments are not specified, the default case will be added, with the name `DEFAULT` and the
                 values of `envs.kmers.k`, `envs.kmers.head`, `envs.kmers.vis_args` and `envs.kmers.devpars`.
+        vj_junc (ns): Arguments for VJ junction circos plots.
+            This analysis is not included in `immunarch`. It is a separate implementation using [`circlize`](https://github.com/jokergoo/circlize).
+            - by: Groupings to show VJ usages. Typically, this is the `Sample` column, so that the VJ usages are shown for each sample.
+                But you can also use other columns, such as `Subject` to show the VJ usages for each subject.
+                Multiple columns should be separated by `,`.
+            - by_clones (flag): If True, the VJ usages will be calculated based on the distinct clonotypes, instead of the individual cells.
+            - subset: Subset the data before plotting VJ usages.
+                The whole data will be expanded to cell level, and then subsetted.
+                Clone sizes will be re-calculated based on the subsetted data, which will affect the VJ usages at cell level (by_clones=False).
+            - devpars (ns): The parameters for the plotting device.
+                - width (type=int): The width of the plot.
+                - height (type=int): The height of the plot.
+                - res (type=int): The resolution of the plot.
+            - cases (type=json;order=9): If you have multiple cases, you can use this argument to specify them.
+                The keys will be used as the names of the cases. The values will be passed to the corresponding arguments above.
+                If any of these arguments are not specified, the values in `envs.vj_junc` will be used.
+                If NO cases are specified, the default case will be added, with the name `DEFAULT` and the
+                values of `envs.vj_junc.by`, `envs.vj_junc.by_clones` `envs.vj_junc.subset` and `envs.vj_junc.devpars`.
     """  # noqa: E501
     input = "immdata:file,metafile:file"
     output = "outdir:dir:{{in.immdata | stem}}.immunarch"
@@ -801,6 +820,14 @@ class Immunarch(Proc):
             },
             "cases": {},
         },
+        # VJ junction
+        "vj_junc": {
+            "by": "Sample",
+            "by_clones": True,
+            "devpars": {"width": 800, "height": 800, "res": 100},
+            "subset": None,
+            "cases": {},
+        },
     }
     script = "file://../scripts/tcr/Immunarch.R"
     plugin_opts = {

biopipen-0.22.0/biopipen/reports/delim/SampleInfo.svelte

@@ -0,0 +1,16 @@
+{% from "utils/misc.liq" import report_jobs, table_of_images -%}
+
+<script>
+    import { Image, DataTable, Descr } from "$libs";
+</script>
+
+{%- macro report_job(job, h=1) -%}
+    {{ job | render_job: h=h }}
+{%- endmacro -%}
+
+{%- macro head_job(job) -%}
+    <h1>{{job.in.infile | stem | escape }}</h1>
+{%- endmacro -%}
+
+{{ report_jobs(jobs, head_job, report_job) }}
+

biopipen-0.22.0/biopipen/reports/scrna/CellsDistribution.svelte

@@ -0,0 +1,15 @@
+{% from "utils/misc.liq" import report_jobs, table_of_images -%}
+<script>
+    import { Image, DataTable, Descr } from "$libs";
+    import { Tabs, Tab, TabContent } from "$ccs";
+</script>
+
+{%- macro report_job(job, h=1) -%}
+    {{ job | render_job: h=h }}
+{%- endmacro -%}
+
+{%- macro head_job(job) -%}
+    <h1>{{job.in.srtobj | stem0 | escape}}</h1>
+{%- endmacro -%}
+
+{{ report_jobs(jobs, head_job, report_job) }}

biopipen-0.22.0/biopipen/reports/scrna/MarkersFinder.svelte

@@ -0,0 +1,18 @@
+{% from "utils/misc.liq" import report_jobs -%}
+
+<script>
+    import { Image, DataTable, Descr } from "$libs";
+    import { Tabs, Tab, TabContent, InlineNotification, Accordion, AccordionItem } from "$ccs";
+</script>
+
+
+{%- macro report_job(job, h=1) -%}
+    {{ job | render_job: h=h }}
+{%- endmacro -%}
+
+
+{%- macro head_job(job) -%}
+    <h1>{{job.in.srtobj | stem | escape}}</h1>
+{%- endmacro -%}
+
+{{ report_jobs(jobs, head_job, report_job) }}

biopipen-0.22.0/biopipen/reports/scrna/MetaMarkers.svelte

@@ -0,0 +1,18 @@
+{% from "utils/misc.liq" import report_jobs, table_of_images -%}
+{% from "utils/gsea.liq" import enrichr_report -%}
+<script>
+    import { Image, DataTable, Descr } from "$libs";
+    import { Tabs, Tab, TabContent, InlineNotification, Accordion, AccordionItem } from "$ccs";
+</script>
+
+
+{%- macro report_job(job, h=1) -%}
+    {{ job | render_job: h=h }}
+{%- endmacro -%}
+
+
+{%- macro head_job(job) -%}
+  <h1>{{job.in.srtobj | stem | escape}}</h1>
+{%- endmacro -%}
+
+{{ report_jobs(jobs, head_job, report_job) }}

biopipen-0.22.0/biopipen/reports/scrna/RadarPlots.svelte

@@ -0,0 +1,15 @@
+{% from "utils/misc.liq" import report_jobs, table_of_images -%}
+<script>
+    import { Image, DataTable } from "$libs";
+    import { Tabs, Tab, TabContent } from "$ccs";
+</script>
+
+{%- macro report_job(job, h=1) -%}
+    {{ job | render_job: h=h }}
+{%- endmacro -%}
+
+{%- macro head_job(job) -%}
+    <h1>{{job.in.srtobj | stem}}</h1>
+{%- endmacro -%}
+
+{{ report_jobs(jobs, head_job, report_job) }}

biopipen-0.22.0/biopipen/reports/scrna/ScFGSEA.svelte

@@ -0,0 +1,16 @@
+{% from "utils/gsea.liq" import fgsea_report -%}
+{% from "utils/misc.liq" import report_jobs -%}
+<script>
+    import { Image, DataTable, Descr } from "$libs";
+    import { Accordion, AccordionItem, Tabs, Tab, TabContent } from "$ccs";
+</script>
+
+{%- macro report_job(job, h=1) -%}
+    {{ job | render_job: h=h }}
+{%- endmacro -%}
+
+{%- macro head_job(job) -%}
+    <h1>{{job.in.srtobj | stem0 | escape}}</h1>
+{%- endmacro -%}
+
+{{ report_jobs(jobs, head_job, report_job) }}

biopipen-0.22.0/biopipen/reports/scrna/SeuratClusterStats.svelte

@@ -0,0 +1,16 @@
+{% from "utils/misc.liq" import report_jobs, table_of_images -%}
+{% from_ os import path %}
+<script>
+    import { DataTable, Image, Descr } from "$libs";
+    import { Tabs, Tab, TabContent } from "$ccs";
+</script>
+
+{%- macro report_job(job, h=1) -%}
+    {{ job | render_job: h=h }}
+{%- endmacro -%}
+
+{%- macro head_job(job) -%}
+    <h1>{{job.in.srtobj | stem}}</h1>
+{%- endmacro -%}
+
+{{ report_jobs(jobs, head_job, report_job) }}

biopipen-0.22.0/biopipen/reports/scrna/SeuratPreparing.svelte

@@ -0,0 +1,15 @@
+{% from "utils/misc.liq" import report_jobs -%}
+{% from_ os import path %}
+<script>
+    import { Image, DataTable, Descr } from "$libs";
+</script>
+
+{%- macro report_job(job, h=1) -%}
+    {{ job | render_job: h=h }}
+{%- endmacro -%}
+
+{%- macro head_job(job) -%}
+    <h1>{{job.in.metafile | stem}}</h1>
+{%- endmacro -%}
+
+{{ report_jobs(jobs, head_job, report_job) }}

biopipen-0.22.0/biopipen/reports/scrna/TopExpressingGenes.svelte

@@ -0,0 +1,17 @@
+{% from "utils/misc.liq" import report_jobs -%}
+{% from "utils/gsea.liq" import enrichr_report -%}
+<script>
+    import { Image, DataTable, Descr } from "$libs";
+    import { Accordion, AccordionItem, Tabs, Tab, TabContent } from "$ccs";
+</script>
+
+
+{%- macro report_job(job, h=1) -%}
+    {{ job | render_job: h=h }}
+{%- endmacro -%}
+
+{%- macro head_job(job) -%}
+  <h1>{{job.in.srtobj | stem0 | escape}}</h1>
+{%- endmacro -%}
+
+{{ report_jobs(jobs, head_job, report_job) }}

biopipen-0.22.0/biopipen/reports/scrna_metabolic_landscape/MetabolicFeatures.svelte

@@ -0,0 +1,32 @@
+{% from "utils/misc.liq" import report_jobs, table_of_images -%}
+{% from "utils/gsea.liq" import fgsea_report_script, fgsea_report, gsea_report -%}
+
+<script>
+  import { Image, DataTable, Descr } from "$libs";
+  import { Tabs, Tab, TabContent, Accordion, AccordionItem } from "$ccs";
+</script>
+
+{%- macro report_job(job, h=2) -%}
+  {% if envs.fgsea %}
+    {{ job | render_job: h=h }}
+  {% else %}
+    {%- for ssdir in job.out.outdir | glob: "*" -%}
+      {%- if basename(ssdir) == "ALL" -%}
+        {%- set h = 1 -%}
+      {%- else -%}
+        <h{{h}}>{{ ssdir | stem }}</h{{h}}>
+      {%- endif -%}
+
+      {% for cldir in ssdir | glob: '*' %}
+        <h{{h+1}}>{{ cldir | basename }}</h{{h+1}}>
+        {{ gsea_report(cldir, h+2, envs, envs.top) }}
+      {% endfor %}
+    {%- endfor -%}
+  {% endif %}
+{%- endmacro -%}
+
+{%- macro head_job(job) -%}
+  <h1>{{job.in.sobjfile | stem | escape}}</h1>
+{%- endmacro -%}
+
+{{ report_jobs(jobs, head_job, report_job) }}

biopipen-0.22.0/biopipen/reports/scrna_metabolic_landscape/MetabolicFeaturesIntraSubset.svelte

@@ -0,0 +1,28 @@
+{% from "utils/misc.liq" import report_jobs, table_of_images -%}
+{% from "utils/gsea.liq" import fgsea_report, gsea_report -%}
+
+<script>
+  import { Image, DataTable, Descr } from "$libs";
+  import { Tabs, Tab, TabContent, Accordion, AccordionItem } from "$ccs";
+</script>
+
+{%- macro report_job(job, h=2) -%}
+  {% if envs.fgsea %}
+    {{ job | render_job: h=h }}
+  {% else %}
+    {%  for groupdir in job.out.outdir | glob: "*" %}
+      <h{{h}}>{{groupdir | basename}}</h{{h}}>
+      {%- set dsdirs = groupdir | glob: "*" -%}
+      {% for dsdir in groupdir | glob: "*" %}
+          <h{{h+1}}>{{ dsdir | basename }}</h{{h+1}}>
+          {{ gsea_report(dsdir, h+2, envs, envs.top) }}
+      {% endfor %}
+    {% endfor %}
+  {% endif %}
+{%- endmacro -%}
+
+{%- macro head_job(job) -%}
+  <h1>{{job.in.sobjfile | stem | escape}}</h1>
+{%- endmacro -%}
+
+{{ report_jobs(jobs, head_job, report_job) }}

biopipen-0.22.0/biopipen/reports/scrna_metabolic_landscape/MetabolicPathwayActivity.svelte

@@ -0,0 +1,89 @@
+{% from "utils/misc.liq" import report_jobs -%}
+
+<script>
+    import { Image, Descr } from "$libs";
+    import { ListItem, UnorderedList } from "$ccs";
+</script>
+
+<h1>Introduction</h1>
+
+<Descr>
+    Metabolic landscape of single cells in the tumor microenvironment.
+</Descr>
+
+<h2>Workflow of the original analysis</h2>
+<Image src="https://github.com/LocasaleLab/Single-Cell-Metabolic-Landscape/raw/master/pipeline.png" />
+
+<h2>Reference</h2>
+<UnorderedList>
+    <ListItem><a href="https://www.nature.com/articles/s41467-019-11738-0" target="_blank">
+        Zhengtao, Ziwei Dai, and Jason W. Locasale.
+        "Metabolic landscape of the tumor microenvironment at single cell resolution."
+        Nature communications 10.1 (2019): 1-12.
+    </a></ListItem>
+    <ListItem><a href="https://github.com/LocasaleLab/Single-Cell-Metabolic-Landscape" target="_blank">
+        Orginal pipeline
+    </a></ListItem>
+</UnorderedList>
+
+<h2>Analyses with this pipeline</h2>
+
+<Descr>
+The cells are grouped at 2 dimensions: `grouping`, usually the cell types, and `subsetting`, usually
+the groups that bring biological meaning (i.e. different timepoints or sample types (tumor/normal)).
+</Descr>
+
+<UnorderedList>
+<ListItem>
+    MetabolicPathwayActivity (this page)
+    <p>Investigating the metabolic pathways of the cells in different groups and subsets.</p>
+    <p>The cells are first grouped by subsets and then the metabolic activities are examined for each groups in different subsets.</p>
+</ListItem>
+<ListItem>
+    <a href="../MetabolicPathwayHeterogeneity/index.html">MetabolicPathwayHeterogeneity</a>
+    <p>Showing metabolic pathways enriched in genes with highest contribution to the metabolic heterogeneities</p>
+</ListItem>
+<ListItem>
+    <a href="../MetabolicFeatures/index.html">MetabolicFeatures</a>
+    <p>Gene set enrichment analysis against the metabolic pathways for groups in different subsets.</p>
+</ListItem>
+<ListItem>
+    <a href="../MetabolicFeaturesIntraSubsets/index.html">MetabolicFeaturesIntraSubsets</a>
+    <p>Gene set enrichment analysis against the metabolic pathways for subsets in different groups.</p>
+</ListItem>
+</UnorderedList>
+
+
+{%- macro report_job(job, h=2) -%}
+  {%- for ssdir in job.out.outdir | glob: "*" -%}
+  {%- if not isdir(ssdir) -%}
+    {%- continue -%}
+  {%- endif -%}
+  <h{{h}}>{{ ssdir | stem }}</h{{h}}>
+
+  <h{{ h+1 }}>Metabolic pathway activities by <code>{{envs.grouping}}</code></h{{ h+1 }}>
+  <Image src="{{ssdir | joinpaths: 'KEGGpathway_activity_heatmap.png'}}" />
+
+  <h{{ h+1 }}>Distributions of pathway activities by <code>{{envs.grouping}}</code></h{{ h+1 }}>
+  <Image src="{{ssdir | joinpaths: 'pathway_activity_violinplot.png'}}" />
+  {%- endfor -%}
+
+  {% if job.out.outdir | glob: "*.group-*.png" -%}
+  <h{{h}}>Merged heatmaps</h{{h}}>
+  {% for group_hm in job.out.outdir | glob: "*.group-*.png" -%}
+    {%- if group_hm.endswith(".group-unclustered.png") -%}
+        <h{{h+1}}>{{group_hm | stem | replace: ".group-unclustered", " (Group Unclustered)"}}</h{{h+1}}>
+        <Image src="{{group_hm}}" />
+    {%- else -%}
+        <h{{h+1}}>{{group_hm | stem | replace: ".group-clustered", " (Group Clustered)"}}</h{{h+1}}>
+        <Image src="{{group_hm}}" />
+    {%- endif -%}
+  {%- endfor -%}
+  {%- endif -%}
+{%- endmacro -%}
+
+{%- macro head_job(job) -%}
+  <h1>{{job.in.sobjfile | stem | escape}}</h1>
+{%- endmacro -%}
+
+{{ report_jobs(jobs, head_job, report_job) }}

biopipen-0.22.0/biopipen/reports/scrna_metabolic_landscape/MetabolicPathwayHeterogeneity.svelte

@@ -0,0 +1,15 @@
+{% from "utils/misc.liq" import report_jobs, table_of_images -%}
+
+<script>
+    import { Image, Descr } from "$libs";
+</script>
+
+{%- macro report_job(job, h=2) -%}
+  {{ job | render_job: h=h }}
+{%- endmacro -%}
+
+{%- macro head_job(job) -%}
+  <h1>{{job.in.sobjfile | stem | escape}}</h1>
+{%- endmacro -%}
+
+{{ report_jobs(jobs, head_job, report_job) }}

biopipen-0.22.0/biopipen/reports/tcr/CDR3AAPhyschem.svelte

@@ -0,0 +1,32 @@
+{% from "utils/misc.liq" import report_jobs -%}
+<script>
+    import { Image, DataTable, Descr } from "$libs";
+    import { Tile, UnorderedList, ListItem, Link, Tabs, Tab, TabContent } from "$ccs";
+</script>
+
+<h1>Summary</h1>
+<Descr>
+    The idea is to perform a regression between two groups of cells
+    (e.g. Treg vs Tconv) at different length of CDR3 AA sequences.
+    The regression will be performed for each physicochemical feature of the
+    AA (hydrophobicity, volume and isolectric point).
+
+    <h2>Reference</h2>
+    <UnorderedList>
+        <ListItem><Link href="https://www.nature.com/articles/ni.3491" target="_blank">https://www.nature.com/articles/ni.3491</Link></ListItem>
+        <ListItem><Link href="https://www.nature.com/articles/s41590-022-01129-x" target="_blank">https://www.nature.com/articles/s41590-022-01129-x</Link></ListItem>
+        <ListItem>Wimley, W. C. &amp; White, S. H. Experimentally determined hydrophobicity scale for proteins at membrane - interfaces. Nat. Struct. Biol. 3, 842-848 (1996).</ListItem>
+        <ListItem>Hdbk of chemistry &amp; physics 72nd edition. (CRC Press, 1991).</ListItem>
+        <ListItem>Zamyatnin, A. A. Protein volume in solution. Prog. Biophys. Mol. Biol. 24, 107-123 (1972).</ListItem>
+    </UnorderedList>
+</Descr>
+
+{%- macro report_job(job, h=1) -%}
+    {{ job | render_job: h=h }}
+{%- endmacro -%}
+
+{%- macro head_job(job) -%}
+    <h1>{{job.out.outdir | stem | replace: ".immunarch", ""}}</h1>
+{%- endmacro -%}
+
+{{ report_jobs(jobs, head_job, report_job) }}