From 0a39340ab9bdbeccb902aad1ce5e45c55bf70e34 Mon Sep 17 00:00:00 2001
From: Ahmet Nihat Simsek <ahmet.n.simsek@gmail.com>
Date: Wed, 5 Mar 2025 09:24:42 +0100
Subject: [PATCH] gene exp reproducibility e2e test

---
 e2e/features/molecular/test_genes.py | 17 +++++++++++++++++
 1 file changed, 17 insertions(+)

diff --git a/e2e/features/molecular/test_genes.py b/e2e/features/molecular/test_genes.py
index 6f4ad0fdd..f5eb31b92 100644
--- a/e2e/features/molecular/test_genes.py
+++ b/e2e/features/molecular/test_genes.py
@@ -85,3 +85,20 @@ def test_no_probes_found_in_concept():
         gene=siibra.vocabularies.GENE_NAMES.G0S2,
     )
     assert features == []
+
+
+@xfail_if_allen_api_unavailable
+def test_reproducing_query_results():
+    region = siibra.get_region("julich 2.9", "V1")
+    features = siibra.features.get(
+        region, siibra.features.molecular.GeneExpressions,
+        gene=siibra.vocabularies.GENE_NAMES.GABARAPL2
+    )
+    assert len(features) == 1
+    gene_exp = features[0]
+    assert set(gene_exp.data['gene']) == {'GABARAPL2'}
+    assert gene_exp.data.iloc[0]['level'] == 10.1143
+    assert gene_exp.data.iloc[2]['zscore'] == -0.6653
+    assert gene_exp.data.iloc[3]['mni_xyz'] == (-6.0, -86.0, 10.0)
+    assert gene_exp.data.iloc[4]['probe_id'] == 1046316
+    assert gene_exp.data.iloc[4]['donor_name'] == "H0351.1015"