Query & integrate data¶

import lamindb as ln
import bionty as bt

ln.track("wukchS8V976U0000")

→ connected lamindb: testuser1/test-facs

→ created Transform('wukchS8V'), started new Run('iPAtW0uV') at 2024-11-21 06:55:40 UTC

→ notebook imports: bionty==0.53.1 lamindb==0.76.16

Inspect the CellMarker registry ¶

Inspect your aggregated cell marker registry as a DataFrame:

bt.CellMarker.df().head()

Show code cell output Hide code cell output

	uid	name	synonyms	description	gene_symbol	ncbi_gene_id	uniprotkb_id	source_id	organism_id	run_id	created_at	created_by_id
id
41	3ZFziy5ims8J	CD14/19	None	None	None	None	None	NaN	1	2	2024-11-21 06:55:33.764984+00:00	1
40	31nZfqQo8yZg	CD103		None	ITGAE	3682	P38570	28.0	1	2	2024-11-21 06:55:33.754551+00:00	1
39	1iLDs6cZIpxj	CD69		None	CD69	969	Q07108	28.0	1	2	2024-11-21 06:55:33.754519+00:00	1
38	525YfNUB967z	CD49B		None	ITGA2	3673	P17301	28.0	1	2	2024-11-21 06:55:33.754486+00:00	1
37	3IPMBjs68Vy1	CXCR4		None	CXCR4	7852	P61073	28.0	1	2	2024-11-21 06:55:33.754447+00:00	1

Search for a marker (synonyms aware):

bt.CellMarker.search("PD-1").df().head(2)

Show code cell output Hide code cell output

	uid	name	synonyms	description	gene_symbol	ncbi_gene_id	uniprotkb_id	source_id	organism_id	run_id	created_at	created_by_id
id
29	33vFR1q26vnM	PD1	PID1\|PD-1\|PD 1	None	PDCD1	5133	A0A0M3M0G7	28	1	1	2024-11-21 06:55:22.086901+00:00	1

Look up markers with auto-complete:

markers = bt.CellMarker.lookup()
markers.cd8

Query panels and collections based on markers, e.g., which collections have 'CD8' in the flow panel:

panels_with_cd8 = ln.FeatureSet.filter(cell_markers=markers.cd8).all()

ln.Artifact.filter(feature_sets__in=panels_with_cd8).df()

Show code cell output Hide code cell output

	uid	version	is_latest	description	key	suffix	type	size	hash	n_objects	n_observations	_hash_type	_accessor	visibility	_key_is_virtual	storage_id	transform_id	run_id	created_at	created_by_id
id
1	HaHgs30aOKQLB9fw0000	None	True	Alpert19	None	.h5ad	dataset	33374864	QNP1c3p6scaAwPo9AW8fLw	None	166537	md5	AnnData	1	True	1	1	1	2024-11-21 06:55:24.902221+00:00	1
2	pNVChGDo9gJecbF00000	None	True	Oetjen18_t1	None	.h5ad	dataset	46506448	WbPHGIMM_5GT68rC8ZydHA	None	241552	md5	AnnData	1	True	1	2	2	2024-11-21 06:55:34.191247+00:00	1

Access registries:

features = ln.Feature.lookup()

Find shared cell markers between two files:

artifacts = ln.Artifact.filter(feature_sets__in=panels_with_cd8).list()

shared_markers = artifacts[0].features["var"] & artifacts[1].features["var"]
shared_markers.list("name")