tooluniverse-expression-data-retrieval

Retrieve gene expression and omics datasets from ArrayExpress and BioStudies with gene disambiguation and quality assessment. Use for finding RNA-seq/microarray datasets by organism/tissue/condition, comparing across studies (case-control, time-series, dose-response), and assessing dataset suitabili

By mims-harvard · 391 installs

npx skills add mims-harvard/tooluniverse --skill tooluniverse-expression-data-retrieval

Source repository · Upstream listing

Gene Expression & Omics Data Retrieval Retrieve gene expression experiments and multi omics datasets with disambiguation and quality assessment. IMPORTANT : Always use English terms in tool calls. Respond in the user's language. LOOK UP DON'T GUESS : Never assume which datasets exist or their accessions. Always search to confirm. Domain Reasoning Before retrieving, determine: organism, tissue, experimental design (case control/time series/dose response). These affect which database to search and how to interpret results. RNA seq provides wider dynamic range; microarray has extensive legacy data. Prioritize experiments with =3 biological replicates, complete annotations, and both raw+processed data. Workflow Phase 0: Clarification (When Needed) Ask ONLY if: gene name ambiguous, tissue/condition unclear, organism not specified. Skip for: specific accessions (E MTAB , E GEOD , S BSST ), clear disease/tissue+organism, explicit platform requests. Phase 1: Query Disambiguation Resolve official gene symbol (HGNC for human, MGI for mouse). Note common aliases for search expansion. User Query Type Search Strategy Specific accession Direct retrieval Gene + condition "[gene] [condition]" + species filter Disease only "[disease]" + species filter Technology specific Add platform keywords Phase 2: Data Retrieval (Internal) Search silently. Do NOT narrate the process. Fallback Chains Primary Fallback ArrayExpress search BioStudies search arrayexpress get experiment biostudies get study arrayexpress get experiment files Note "Files unavailable" Phase 3: Report Dataset Profile Present as a Dataset Search Report . Hide search process. Include: 1. Search Summary : query, databases searched, result count 2. Top Experiments (per experiment): Accession, organism, type (RNA seq/microarray), platform, sample count, date Description, experimental design (conditions, replicates, tissue) Sample groups table, data files table Quality assessment (●●●/●●○/●○○) 3. Multi Omics Studies (from BioStudies): accession, type, data types included 4. Summary Table : all experiments ranked 5. Recommendations : best dataset for user's purpose, integration notes 6. Data Access : download links, database URLs Data Quality Tiers Tier Symbol Criteria High ●●● =3 bio replicates, complete metadata, processed data available Medium ●●○ 2 3 replicates OR some metadata gaps Low ●○○ No replicates, sparse metadata, or access issues Caution ○○○ Single sample, no replication, outdated platform Reasoning Framework Dataset quality : Prioritize =3 biological replicates, complete annotations, both raw+processed data. Single replicate experiments can inform but not be sole evidence. Platform comparison : RNA seq = wider dynamic range, novel transcripts. Microarray = probe limited but extensive legacy data. Cross platform combining requires batch correction. Metadata scoring : Rate 0 5 on: (1) sample annotations, (2) design documented, (3) pipeline described, (4) raw data deposited, (5) publication linked. Score <=2 warrants caution. GEO vs ArrayExpress : GEO has broader coverage (older studies); ArrayExpress enforces stricter metadata. BioStudies captures multi omics. Search both. Synthesis Questions 1. Does the dataset have sufficient replication and metadata for the intended analysis? 2. Are there batch effects or confounding variables? 3. Do multiple datasets show concordant patterns, and can they be integrated? Error Handling Error Response "No experiments found" Broaden keywords, remove species filter, try synonyms "Accession not found" Verify format, check if withdrawn "Files not available" Note: "Data files restricted by submitter" "API timeout" Retry once, note "(metadata retrieval incomplete)" Tool Reference ArrayExpress : arrayexpress search experiments (search), arrayexpress get experiment (metadata), arrayexpress get experiment files (downloads), arrayexpress get experiment samples (annotations) BioStudies : biostudies search (search), biostudies get study (metadata+sections), biostudies get study files (files) Additional Sources : GEO search rnaseq datasets / geo search datasets GEO (largest RNA seq repo) OmicsDI search datasets cross repository aggregation (GEO+ArrayExpress+PRIDE+MassIVE) GTEx get expression summary baseline tissue expression (54 normal tissues, param: gene symbol ) ENAPortal search studies sequencing studies (param: query with description="..." ) CxGDisc search datasets single cell datasets (needs exact disease ontology terms) GxA list experiments / GxA get experiment info / GxA get experiment expression EBI Gene Expression Atlas: curated baseline (tissue) and differential (condition comparison) bulk RNA seq/microarray experiments. Filter GxA list experiments by species and experiment type ("baseline" vs "differential") first; GxA get experiment expression needs a specific experiment accession (e.g. E MTAB 2836 ) plus gene id . ARCHS4 get gene expression / ARCHS4 get gene correlations ARCHS4: box plot expression statistics (min/Q1/median/Q3/max, log2) across human/mouse tissues and cell lines, and top co expressed genes by Pearson correlation across 300K+ uniformly reprocessed RNA seq samples. Good for a quick baseline expression sanity check or finding candidate co regulated genes without running a full differential expression pipeline. Bgee search genes / Bgee get gene expression / Bgee list species Bgee: curated baseline expression by tissue calls across 29+ animal species (vs. ARCHS4's human/mouse only reprocessed RNA seq, or GxA's curated but narrower experiment set). Bgee get gene expression needs BOTH an Ensembl gene ID and an NCBI taxonomy ID (e.g. gene id="ENSG00000141510" , species id="9606" for human TP53) resolve the gene ID with Bgee search genes first if you only have a symbol. Use Bgee list species to confirm a species is covered before assuming Bgee has it. Best for cross species comparison questions ("is gene X expressed in the same tissues in mouse and zebrafish?"), not for finding raw sequencing datasets to download. PubMed search articles dataset discovery via publications Search Parameters ArrayExpress : keywords (free text), species (scientific name), array (platform filter), limit BioStudies : query (free text), limit