Project flow¶
LaminDB allows tracking data lineage on the entire project level.
Here, we walk through exemplified app uploads, pipelines & notebooks following Schmidt et al., 2022.
A CRISPR screen reading out a phenotypic endpoint on T cells is paired with scRNA-seq to generate insights into IFN-γ production.
These insights get linked back to the original data through the steps taken in the project to provide context for interpretation & future decision making.
More specifically: Why should I care about data flow?
Data flow tracks data sources & transformations to trace biological insights, verify experimental outcomes, meet regulatory standards, increase the robustness of research and optimize the feedback loop of team-wide learning iterations.
While tracking data flow is easier when it’s governed by deterministic pipelines, it becomes hard when it’s governed by interactive human-driven analyses.
LaminDB interfaces workflow mangers for the former and embraces the latter.
# !pip install 'lamindb[jupyter,bionty,aws]'
!lamin init --storage ./mydata
Show code cell output
💡 connected lamindb: testuser1/mydata
Import lamindb:
import lamindb as ln
from IPython.display import Image, display
💡 connected lamindb: testuser1/mydata
Steps¶
In the following, we walk through exemplified steps covering different types of transforms (Transform).
Note
The full notebooks are in this repository.
App upload of phenotypic data 
¶
Register data through app upload from wetlab by testuser1:
# This function mimics the upload of artifacts via the UI
# In reality, you simply drag and drop files into the UI
def mock_upload_crispra_result_app():
ln.setup.login("testuser1")
transform = ln.Transform(name="Upload GWS CRISPRa result", type="upload")
ln.track(transform=transform)
output_path = ln.core.datasets.schmidt22_crispra_gws_IFNG(ln.settings.storage.root)
output_file = ln.Artifact(
output_path, description="Raw data of schmidt22 crispra GWS"
)
output_file.save()
mock_upload_crispra_result_app()
Show code cell output
💡 created Transform('8uwpz71JcxRcjhZ2') & created Run('2024-08-16 09:34:51.873982+00:00')
Hit identification in notebook 
¶
Access, transform & register data in drylab by testuser2 in notebook hit-identification.
Show code cell content
# the following mimics the integrated analysis notebook
# In reality, you would execute inside the notebook
import nbproject_test
from pathlib import Path
cwd = Path.cwd()
nbproject_test.execute_notebooks(cwd / "project-flow-scripts/hit-identification.ipynb", write=True)
Executing notebooks in /home/runner/work/lamin-usecases/lamin-usecases/docs/project-flow-scripts/hit-identification.ipynb
Scheduled: ['hit-identification']
hit-identification
✓ (4.229s)
Total time: 4.230s
Inspect data flow:
artifact = ln.Artifact.filter(description="hits from schmidt22 crispra GWS").one()
artifact.view_lineage()
Sequencer upload 
¶
Upload files from sequencer via script chromium_10x_upload.py:
!python project-flow-scripts/chromium_10x_upload.py
Show code cell output
💡 connected lamindb: testuser1/mydata
💡 created Transform('qCJPkOuZAi9q0000') & created Run('2024-08-16 09:34:59.331636+00:00')
scRNA-seq bioinformatics pipeline 
¶
Process uploaded files using a script or workflow manager: Pipelines – workflow managers and obtain 3 output files in a directory filtered_feature_bc_matrix/:
!python project-flow-scripts/cellranger.py
Show code cell output
💡 connected lamindb: testuser1/mydata
💡 created Transform('pNctUrJLlBe8Cio6') & created Run('2024-08-16 09:35:02.270144+00:00')
❗ this creates one artifact per file in the directory - consider ln.Artifact(dir_path) to get one artifact for the entire directory
!python project-flow-scripts/postprocess_cellranger.py
Show code cell output
💡 connected lamindb: testuser1/mydata
💡 created Transform('YqmbO6oMXjRj0000') & created Run('2024-08-16 09:35:04.101752+00:00')
Inspect data flow:
output_file = ln.Artifact.filter(description="perturbseq counts").one()
output_file.view_lineage()
Integrate scRNA-seq & phenotypic data ¶
Integrate data in notebook integrated-analysis.
Show code cell content
# the following mimics the integrated analysis notebook
# In reality, you would execute inside the notebook
nbproject_test.execute_notebooks(cwd / "project-flow-scripts/integrated-analysis.ipynb", write=True)
Executing notebooks in /home/runner/work/lamin-usecases/lamin-usecases/docs/project-flow-scripts/integrated-analysis.ipynb
Scheduled: ['integrated-analysis']
integrated-analysis
✓ (4.079s)
Total time: 4.081s
Review results¶
Let’s load one of the plots:
# track the current notebook as transform
ln.context.uid = "1LCd8kco9lZU0000"
ln.context.track()
💡 notebook imports: ipython==8.26.0 lamindb==0.76.0 nbproject_test==0.5.1
💡 created Transform('1LCd8kco9lZU0000') & created Run('2024-08-16 09:35:10.511787+00:00')
artifact = ln.Artifact.filter(key__contains="figures/matrixplot").one()
artifact.cache()
Show code cell output
PosixUPath('/home/runner/work/lamin-usecases/lamin-usecases/docs/mydata/.lamindb/WQ0o6uSSQQRpMq2afZfj.png')
display(Image(filename=artifact.path))
We see that the image artifact is tracked as an input of the current notebook. The input is highlighted, the notebook follows at the bottom:
artifact.view_lineage()
Alternatively, we can also look at the sequence of transforms:
transform = ln.Transform.search("Project flow").first()
transform.predecessors.df()
| uid | version | is_latest | name | key | description | type | source_code | hash | reference | reference_type | _source_code_artifact_id | created_by_id | updated_at | |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| id | ||||||||||||||
| 6 | lB3IyPLQSmvt0000 | None | True | Perform single cell analysis, integrate with C... | integrated-analysis.ipynb | None | notebook | None | None | None | None | None | 2 | 2024-08-16 09:35:08.863961+00:00 |
transform.view_lineage()
Understand runs¶
We tracked pipeline and notebook runs through run_context, which stores a Transform and a Run record as a global context.
Artifact objects are the inputs and outputs of runs.
What if I don’t want a global context?
Sometimes, we don’t want to create a global run context but manually pass a run when creating an artifact:
run = ln.Run(transform=transform)
ln.Artifact(filepath, run=run)
When does an artifact appear as a run input?
When accessing an artifact via cache(), load() or open(), two things happen:
The current run gets added to
artifact.input_ofThe transform of that artifact gets added as a parent of the current transform
You can then switch off auto-tracking of run inputs if you set ln.settings.track_run_inputs = False: Can I disable tracking run inputs?
You can also track run inputs on a case by case basis via is_run_input=True, e.g., here:
artifact.load(is_run_input=True)
Query by provenance¶
We can query or search for the notebook that created the artifact:
transform = ln.Transform.search("GWS CRIPSRa analysis").first()
And then find all the artifacts created by that notebook:
ln.Artifact.filter(transform=transform).df()
| uid | version | is_latest | description | key | suffix | type | _accessor | size | hash | _hash_type | n_objects | n_observations | visibility | _key_is_virtual | storage_id | transform_id | run_id | created_by_id | updated_at | |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| id | ||||||||||||||||||||
| 2 | LlmWdQBpAV8FyqOH1FP3 | None | True | hits from schmidt22 crispra GWS | None | .parquet | dataset | DataFrame | 18368 | ky0pLx8o9oFJnn0niPUHog | md5 | None | None | 1 | True | 1 | 2 | 2 | 2 | 2024-08-16 09:34:57.020416+00:00 |
Which transform ingested a given artifact?
artifact = ln.Artifact.filter().first()
artifact.transform
Transform(uid='8uwpz71JcxRcjhZ2', is_latest=True, name='Upload GWS CRISPRa result', type='upload', created_by_id=1, updated_at='2024-08-16 09:34:51 UTC')
And which user?
artifact.created_by
User(uid='DzTjkKse', handle='testuser1', name='Test User1', updated_at='2024-08-16 09:34:59 UTC')
Which transforms were created by a given user?
users = ln.User.lookup()
ln.Transform.filter(created_by=users.testuser1).df()
| uid | version | is_latest | name | key | description | type | source_code | hash | reference | reference_type | _source_code_artifact_id | created_by_id | updated_at | |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| id | ||||||||||||||
| 1 | 8uwpz71JcxRcjhZ2 | None | True | Upload GWS CRISPRa result | None | None | upload | None | None | None | None | NaN | 1 | 2024-08-16 09:34:51.871644+00:00 |
| 3 | qCJPkOuZAi9q0000 | None | True | chromium_10x_upload.py | chromium_10x_upload.py | None | script | None | None | None | None | 5.0 | 1 | 2024-08-16 09:34:59.710933+00:00 |
| 7 | 1LCd8kco9lZU0000 | None | True | Project flow | project-flow.ipynb | None | notebook | None | None | None | None | NaN | 1 | 2024-08-16 09:35:10.507303+00:00 |
Which notebooks were created by a given user?
ln.Transform.filter(created_by=users.testuser1, type="notebook").df()
| uid | version | is_latest | name | key | description | type | source_code | hash | reference | reference_type | _source_code_artifact_id | created_by_id | updated_at | |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| id | ||||||||||||||
| 7 | 1LCd8kco9lZU0000 | None | True | Project flow | project-flow.ipynb | None | notebook | None | None | None | None | None | 1 | 2024-08-16 09:35:10.507303+00:00 |
We can also view all recent additions to the entire database:
ln.view()
Show code cell output
Artifact
| uid | version | is_latest | description | key | suffix | type | _accessor | size | hash | _hash_type | n_objects | n_observations | visibility | _key_is_virtual | storage_id | transform_id | run_id | created_by_id | updated_at | |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| id | ||||||||||||||||||||
| 12 | WQ0o6uSSQQRpMq2afZfj | None | True | None | figures/matrixplot_fig2_score-wgs-hits-per-clu... | .png | None | None | 28814 | Uqpe3rI2qCa6KHvUIxadLw | md5 | None | None | 1 | True | 1 | 6 | 6 | 2 | 2024-08-16 09:35:09.719789+00:00 |
| 11 | T9o23F8CrsO5OERWNM3Z | None | True | None | figures/umap_fig1_score-wgs-hits.png | .png | None | None | 118999 | JMFvnvCcQzIdtXM8Y12MKg | md5 | None | None | 1 | True | 1 | 6 | 6 | 2 | 2024-08-16 09:35:09.502809+00:00 |
| 10 | g95uLDa2EAgq9aVuZRFn | None | True | perturbseq counts | schmidt22_perturbseq.h5ad | .h5ad | None | AnnData | 20659936 | la7EvqEUMDlug9-rpw-udA | md5 | None | None | 1 | False | 1 | 5 | 5 | 2 | 2024-08-16 09:35:05.784807+00:00 |
| 9 | 9CO1qIb3mXlXYlPfYYVe | None | True | None | perturbseq/filtered_feature_bc_matrix/barcodes... | .tsv.gz | None | None | 6 | S4dHcjum96m_9nLvClE0Fw | md5 | None | None | 1 | False | 1 | 4 | 4 | 2 | 2024-08-16 09:35:02.638021+00:00 |
| 8 | 8LWw1SucoJSPW2qxfr1C | None | True | None | perturbseq/filtered_feature_bc_matrix/matrix.m... | .mtx.gz | None | None | 6 | 9_IXVmuUm6zUsmO5rwUwCQ | md5 | None | None | 1 | False | 1 | 4 | 4 | 2 | 2024-08-16 09:35:02.637522+00:00 |
| 7 | QAfoO5RXrp6Qwld3irov | None | True | None | perturbseq/filtered_feature_bc_matrix/features... | .tsv.gz | None | None | 6 | Q4ym-bXrn5ewjfCS8YGYTw | md5 | None | None | 1 | False | 1 | 4 | 4 | 2 | 2024-08-16 09:35:02.636721+00:00 |
| 4 | 5vxv2eLvZJRkSTYdiq20 | None | True | None | fastq/perturbseq_R2_001.fastq.gz | .fastq.gz | None | None | 6 | XJ2KE5YTfP6hswybpXO1Vg | md5 | None | None | 1 | False | 1 | 3 | 3 | 1 | 2024-08-16 09:34:59.691350+00:00 |
Run
| uid | started_at | finished_at | is_consecutive | reference | reference_type | transform_id | report_id | environment_id | parent_id | created_by_id | |
|---|---|---|---|---|---|---|---|---|---|---|---|
| id | |||||||||||
| 1 | fGiJqifsTLhDBAuEENY4 | 2024-08-16 09:34:51.874128+00:00 | NaT | True | None | None | 1 | None | NaN | None | 1 |
| 2 | EZzWgz5dWMpQ8x8s8EmX | 2024-08-16 09:34:56.614499+00:00 | NaT | True | None | None | 2 | None | NaN | None | 2 |
| 3 | SpEHFbJMQ8GaVoATSURT | 2024-08-16 09:34:59.331804+00:00 | 2024-08-16 09:34:59.707660+00:00 | True | None | None | 3 | None | 6.0 | None | 1 |
| 4 | Pt3RzRKMlx15EqlmuIOy | 2024-08-16 09:35:02.270280+00:00 | NaT | None | None | None | 4 | None | NaN | None | 2 |
| 5 | B68ZWMHKyqriAzz1Wu0G | 2024-08-16 09:35:04.101900+00:00 | NaT | None | None | None | 5 | None | NaN | None | 2 |
| 6 | 0CQOeIp6hzh3A35ymgUl | 2024-08-16 09:35:08.868926+00:00 | NaT | True | None | None | 6 | None | NaN | None | 2 |
| 7 | G52qUVNHmY9cVSaOKZnU | 2024-08-16 09:35:10.511912+00:00 | NaT | True | None | None | 7 | None | NaN | None | 1 |
Storage
| uid | root | description | type | region | instance_uid | run_id | created_by_id | updated_at | |
|---|---|---|---|---|---|---|---|---|---|
| id | |||||||||
| 1 | YAgd422VemAH | /home/runner/work/lamin-usecases/lamin-usecase... | None | local | None | 54ZGqgkROOFf | None | 1 | 2024-08-16 09:34:50.040221+00:00 |
Transform
| uid | version | is_latest | name | key | description | type | source_code | hash | reference | reference_type | _source_code_artifact_id | created_by_id | updated_at | |
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| id | ||||||||||||||
| 7 | 1LCd8kco9lZU0000 | None | True | Project flow | project-flow.ipynb | None | notebook | None | None | None | None | NaN | 1 | 2024-08-16 09:35:10.507303+00:00 |
| 6 | lB3IyPLQSmvt0000 | None | True | Perform single cell analysis, integrate with C... | integrated-analysis.ipynb | None | notebook | None | None | None | None | NaN | 2 | 2024-08-16 09:35:08.863961+00:00 |
| 5 | YqmbO6oMXjRj0000 | None | True | postprocess_cellranger.py | postprocess_cellranger.py | None | script | None | None | None | None | NaN | 2 | 2024-08-16 09:35:04.099031+00:00 |
| 4 | pNctUrJLlBe8Cio6 | 7.2.0 | True | Cell Ranger | None | None | pipeline | None | None | https://www.10xgenomics.com/support/software/c... | None | NaN | 2 | 2024-08-16 09:35:02.267131+00:00 |
| 3 | qCJPkOuZAi9q0000 | None | True | chromium_10x_upload.py | chromium_10x_upload.py | None | script | None | None | None | None | 5.0 | 1 | 2024-08-16 09:34:59.710933+00:00 |
| 2 | T0T28btuB0PG0000 | None | True | GWS CRIPSRa analysis | hit-identification.ipynb | None | notebook | None | None | None | None | NaN | 2 | 2024-08-16 09:34:56.611231+00:00 |
| 1 | 8uwpz71JcxRcjhZ2 | None | True | Upload GWS CRISPRa result | None | None | upload | None | None | None | None | NaN | 1 | 2024-08-16 09:34:51.871644+00:00 |
User
| uid | handle | name | updated_at | |
|---|---|---|---|---|
| id | ||||
| 2 | bKeW4T6E | testuser2 | Test User2 | 2024-08-16 09:35:02.258320+00:00 |
| 1 | DzTjkKse | testuser1 | Test User1 | 2024-08-16 09:34:59.220626+00:00 |
Show code cell content
!lamin login testuser1
!lamin delete --force mydata
!rm -r ./mydata
✅ logged in with email testuser1@lamin.ai (uid: DzTjkKse)
Traceback (most recent call last):
File "/opt/hostedtoolcache/Python/3.10.14/x64/bin/lamin", line 8, in <module>
sys.exit(main())
File "/opt/hostedtoolcache/Python/3.10.14/x64/lib/python3.10/site-packages/rich_click/rich_command.py", line 367, in __call__
return super().__call__(*args, **kwargs)
File "/opt/hostedtoolcache/Python/3.10.14/x64/lib/python3.10/site-packages/click/core.py", line 1157, in __call__
return self.main(*args, **kwargs)
File "/opt/hostedtoolcache/Python/3.10.14/x64/lib/python3.10/site-packages/rich_click/rich_command.py", line 152, in main
rv = self.invoke(ctx)
File "/opt/hostedtoolcache/Python/3.10.14/x64/lib/python3.10/site-packages/click/core.py", line 1688, in invoke
return _process_result(sub_ctx.command.invoke(sub_ctx))
File "/opt/hostedtoolcache/Python/3.10.14/x64/lib/python3.10/site-packages/click/core.py", line 1434, in invoke
return ctx.invoke(self.callback, **ctx.params)
File "/opt/hostedtoolcache/Python/3.10.14/x64/lib/python3.10/site-packages/click/core.py", line 783, in invoke
return __callback(*args, **kwargs)
File "/opt/hostedtoolcache/Python/3.10.14/x64/lib/python3.10/site-packages/lamin_cli/__main__.py", line 171, in delete
return delete(instance, force=force)
File "/opt/hostedtoolcache/Python/3.10.14/x64/lib/python3.10/site-packages/lamindb_setup/_delete.py", line 98, in delete
n_objects = check_storage_is_empty(
File "/opt/hostedtoolcache/Python/3.10.14/x64/lib/python3.10/site-packages/lamindb_setup/core/upath.py", line 776, in check_storage_is_empty
raise InstanceNotEmpty(message)
lamindb_setup.core.upath.InstanceNotEmpty: Storage /home/runner/work/lamin-usecases/lamin-usecases/docs/mydata/.lamindb contains 5 objects ('_is_initialized' ignored) - delete them prior to deleting the instance