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('n9wwSEOV') at 2024-12-20 15:06:49 UTC

→ notebook imports: bionty==0.53.2 lamindb==0.77.3

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-12-20 15:06:43.586359+00:00	1
40	31nZfqQo8yZg	CD103		None	ITGAE	3682	P38570	28.0	1	2	2024-12-20 15:06:43.576793+00:00	1
39	1iLDs6cZIpxj	CD69		None	CD69	969	Q07108	28.0	1	2	2024-12-20 15:06:43.576761+00:00	1
38	525YfNUB967z	CD49B		None	ITGA2	3673	P17301	28.0	1	2	2024-12-20 15:06:43.576729+00:00	1
37	3IPMBjs68Vy1	CXCR4		None	CXCR4	7852	P61073	28.0	1	2	2024-12-20 15:06:43.576689+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-12-20 15:06:31.623098+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	key	description	suffix	type	size	hash	n_objects	n_observations	_hash_type	_accessor	visibility	_key_is_virtual	storage_id	transform_id	version	is_latest	run_id	created_at	created_by_id
id
1	l2w0Itc5044I3LuP0000	None	Alpert19	.h5ad	dataset	33374864	QNP1c3p6scaAwPo9AW8fLw	None	166537	md5	AnnData	1	True	1	1	None	True	1	2024-12-20 15:06:34.437856+00:00	1
2	0HzoLuLmCFhJiQBE0000	None	Oetjen18_t1	.h5ad	dataset	46506448	WbPHGIMM_5GT68rC8ZydHA	None	241552	md5	AnnData	1	True	1	2	None	True	2	2024-12-20 15:06:43.965714+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")