SIMO Multiomics Integration Agent
The SIMO Multiomics Integration Agent performs spatial integration of multi-omics datasets through probabilistic alignment. Unlike previous tools limited to transcriptomics, SIMO integrates spatial transcriptomics with single-cell RNA-seq and expands to chromatin accessibility, DNA methylation, and proteomics data.
When to Use This Skill
- When integrating spatial transcriptomics with single-cell multi-omics data.
- For constructing comprehensive tissue atlases with spatial context.
- To map epigenomic states (ATAC-seq, methylation) onto spatial coordinates.
- When analyzing multi-modal cellular phenotypes in tissue architecture.
- For spatial deconvolution combining multiple modalities.
Core Capabilities
Spatial-scRNA Integration: Probabilistically align single-cell RNA-seq to spatial coordinates.
Chromatin Accessibility Mapping: Project scATAC-seq profiles onto spatial tissue locations.
DNA Methylation Spatial Mapping: Integrate single-cell methylation data with spatial context.
Multi-Modal Fusion: Combine transcriptomic, epigenomic, and proteomic layers.
Probabilistic Cell-Type Assignment: Assign cell types to spatial spots with uncertainty quantification.
Spatial Niche Identification: Discover cellular niches defined by multi-omic signatures.
Supported Modalities
| Modality |
Input Format |
Spatial Reference |
| scRNA-seq |
AnnData, Seurat |
Visium, MERFISH, Xenium |
| scATAC-seq |
SnapATAC2, ArchR |
Visium, Slide-seq |
| scMethyl |
Bismark, allcools |
Any spatial modality |
| CITE-seq (protein) |
AnnData |
Spatial proteomics |
| Multi-ome (RNA+ATAC) |
Muon, SnapATAC2 |
All platforms |
Integration Algorithm
| Step |
Method |
Purpose |
| Feature Selection |
HVG + marker genes |
Reduce dimensionality |
| Embedding |
Variational autoencoder |
Shared latent space |
| Alignment |
Optimal transport |
Probabilistic matching |
| Spatial Mapping |
Gaussian processes |
Smooth spatial predictions |
| Uncertainty |
Posterior sampling |
Confidence intervals |
Workflow
Input: Spatial transcriptomics (Visium/MERFISH/Xenium), reference single-cell multi-omics.
Preprocessing: Normalize, select features, QC both datasets.
Embedding: Learn joint latent representation across modalities.
Probabilistic Alignment: Compute cell-to-spot assignment probabilities.
Spatial Imputation: Transfer modalities to spatial coordinates.
Niche Analysis: Identify spatial domains by multi-omic signatures.
Output: Integrated spatial multi-omics object, niche assignments, visualizations.
Example Usage
User: "Integrate our scRNA-seq and scATAC-seq data with the spatial transcriptomics to understand chromatin states in different tissue regions."
Agent Action:
python3 Skills/Genomics/SIMO_Multiomics_Integration_Agent/simo_integration.py \
--spatial_data visium_data.h5ad \
--scrna_ref scrna_atlas.h5ad \
--scatac_ref scatac_atlas.h5ad \
--modalities rna,atac \
--n_spots_per_cell 5 \
--uncertainty_quantification true \
--output integrated_spatial_multiome.h5ad
Output Components
| Output |
Description |
Format |
| Integrated Object |
Multi-modal spatial data |
AnnData/Muon |
| Cell Type Map |
Spatial cell type assignments |
GeoTIFF, CSV |
| Chromatin Accessibility Map |
Spatial ATAC patterns |
BigWig, CSV |
| Niche Assignments |
Spatial domain labels |
CSV, Zarr |
| Uncertainty Maps |
Per-spot confidence |
GeoTIFF |
| Gene Activity Scores |
ATAC-derived gene activity |
AnnData layer |
Spatial Platforms Supported
| Platform |
Resolution |
Spots/Cells |
Genes |
| 10x Visium |
55 μm |
~5,000 |
Whole transcriptome |
| 10x Visium HD |
8 μm |
~300,000 |
Whole transcriptome |
| 10x Xenium |
Subcellular |
>100,000 |
300-5,000 panel |
| MERFISH |
Subcellular |
>1M |
100-10,000 panel |
| Slide-seq |
10 μm |
~60,000 |
Whole transcriptome |
| CosMx |
Subcellular |
>1M |
1,000-6,000 panel |
AI/ML Components
Variational Integration:
- Multi-modal VAE for joint embeddings
- Contrastive learning for modality alignment
- Batch correction across datasets
Probabilistic Mapping:
- Optimal transport with entropic regularization
- Gaussian process spatial smoothing
- Bayesian uncertainty estimation
Niche Discovery:
- Multi-view clustering
- Spatial autocorrelation (Moran's I)
- Graph neural networks for niche boundaries
Prerequisites
- Python 3.10+
- Scanpy, Squidpy, Muon
- scvi-tools, SnapATAC2
- POT (Python Optimal Transport)
- PyTorch, GPyTorch
Related Skills
- scGPT_Agent - For foundation model embeddings
- Spatial_Epigenomics_Agent - For spatial epigenomics analysis
- Cell_Cell_Communication - For ligand-receptor analysis
- Nicheformer_Spatial_Agent - For spatial niche modeling
Special Considerations
- Batch Effects: Pre-align datasets from different protocols
- Spot Deconvolution: Lower resolution platforms need deconvolution
- Sparsity: scATAC data requires aggregation strategies
- Compute: Multi-modal integration is memory-intensive
- Validation: Verify spatial patterns with known marker distributions
Applications
| Application |
Use Case |
| Tumor Microenvironment |
Map chromatin states of immune infiltrates |
| Development |
Track lineage chromatin dynamics spatially |
| Neurodegeneration |
Spatial mapping of epigenetic changes |
| Fibrosis |
Understand spatial activation programs |
Author
AI Group - Biomedical AI Platform
1---2name: simo-multiomics-integration-agent3description: AI-powered spatial integration of multi-omics datasets using probabilistic alignment for comprehensive tissue atlas construction and cellular state mapping.4---56<!--7# COPYRIGHT NOTICE8# This file is part of the "Universal Biomedical Skills" project.9# Copyright (c) 2026 MD BABU MIA, PhD <md.babu.mia@mssm.edu>10# All Rights Reserved.11#12# This code is proprietary and confidential.13# Unauthorized copying of this file, via any medium is strictly prohibited.14#15# Provenance: Authenticated by MD BABU MIA1617-->18192021# SIMO Multiomics Integration Agent2223The **SIMO Multiomics Integration Agent** performs spatial integration of multi-omics datasets through probabilistic alignment. Unlike previous tools limited to transcriptomics, SIMO integrates spatial transcriptomics with single-cell RNA-seq and expands to chromatin accessibility, DNA methylation, and proteomics data.2425## When to Use This Skill2627* When integrating spatial transcriptomics with single-cell multi-omics data.28* For constructing comprehensive tissue atlases with spatial context.29* To map epigenomic states (ATAC-seq, methylation) onto spatial coordinates.30* When analyzing multi-modal cellular phenotypes in tissue architecture.31* For spatial deconvolution combining multiple modalities.3233## Core Capabilities34351. **Spatial-scRNA Integration**: Probabilistically align single-cell RNA-seq to spatial coordinates.36372. **Chromatin Accessibility Mapping**: Project scATAC-seq profiles onto spatial tissue locations.38393. **DNA Methylation Spatial Mapping**: Integrate single-cell methylation data with spatial context.40414. **Multi-Modal Fusion**: Combine transcriptomic, epigenomic, and proteomic layers.42435. **Probabilistic Cell-Type Assignment**: Assign cell types to spatial spots with uncertainty quantification.44456. **Spatial Niche Identification**: Discover cellular niches defined by multi-omic signatures.4647## Supported Modalities4849| Modality | Input Format | Spatial Reference |50|----------|--------------|-------------------|51| scRNA-seq | AnnData, Seurat | Visium, MERFISH, Xenium |52| scATAC-seq | SnapATAC2, ArchR | Visium, Slide-seq |53| scMethyl | Bismark, allcools | Any spatial modality |54| CITE-seq (protein) | AnnData | Spatial proteomics |55| Multi-ome (RNA+ATAC) | Muon, SnapATAC2 | All platforms |5657## Integration Algorithm5859| Step | Method | Purpose |60|------|--------|---------|61| Feature Selection | HVG + marker genes | Reduce dimensionality |62| Embedding | Variational autoencoder | Shared latent space |63| Alignment | Optimal transport | Probabilistic matching |64| Spatial Mapping | Gaussian processes | Smooth spatial predictions |65| Uncertainty | Posterior sampling | Confidence intervals |6667## Workflow68691. **Input**: Spatial transcriptomics (Visium/MERFISH/Xenium), reference single-cell multi-omics.70712. **Preprocessing**: Normalize, select features, QC both datasets.72733. **Embedding**: Learn joint latent representation across modalities.74754. **Probabilistic Alignment**: Compute cell-to-spot assignment probabilities.76775. **Spatial Imputation**: Transfer modalities to spatial coordinates.78796. **Niche Analysis**: Identify spatial domains by multi-omic signatures.80817. **Output**: Integrated spatial multi-omics object, niche assignments, visualizations.8283## Example Usage8485**User**: "Integrate our scRNA-seq and scATAC-seq data with the spatial transcriptomics to understand chromatin states in different tissue regions."8687**Agent Action**:88```bash89python3 Skills/Genomics/SIMO_Multiomics_Integration_Agent/simo_integration.py \90 --spatial_data visium_data.h5ad \91 --scrna_ref scrna_atlas.h5ad \92 --scatac_ref scatac_atlas.h5ad \93 --modalities rna,atac \94 --n_spots_per_cell 5 \95 --uncertainty_quantification true \96 --output integrated_spatial_multiome.h5ad97```9899## Output Components100101| Output | Description | Format |102|--------|-------------|--------|103| Integrated Object | Multi-modal spatial data | AnnData/Muon |104| Cell Type Map | Spatial cell type assignments | GeoTIFF, CSV |105| Chromatin Accessibility Map | Spatial ATAC patterns | BigWig, CSV |106| Niche Assignments | Spatial domain labels | CSV, Zarr |107| Uncertainty Maps | Per-spot confidence | GeoTIFF |108| Gene Activity Scores | ATAC-derived gene activity | AnnData layer |109110## Spatial Platforms Supported111112| Platform | Resolution | Spots/Cells | Genes |113|----------|------------|-------------|-------|114| 10x Visium | 55 μm | ~5,000 | Whole transcriptome |115| 10x Visium HD | 8 μm | ~300,000 | Whole transcriptome |116| 10x Xenium | Subcellular | >100,000 | 300-5,000 panel |117| MERFISH | Subcellular | >1M | 100-10,000 panel |118| Slide-seq | 10 μm | ~60,000 | Whole transcriptome |119| CosMx | Subcellular | >1M | 1,000-6,000 panel |120121## AI/ML Components122123**Variational Integration**:124- Multi-modal VAE for joint embeddings125- Contrastive learning for modality alignment126- Batch correction across datasets127128**Probabilistic Mapping**:129- Optimal transport with entropic regularization130- Gaussian process spatial smoothing131- Bayesian uncertainty estimation132133**Niche Discovery**:134- Multi-view clustering135- Spatial autocorrelation (Moran's I)136- Graph neural networks for niche boundaries137138## Prerequisites139140* Python 3.10+141* Scanpy, Squidpy, Muon142* scvi-tools, SnapATAC2143* POT (Python Optimal Transport)144* PyTorch, GPyTorch145146## Related Skills147148* scGPT_Agent - For foundation model embeddings149* Spatial_Epigenomics_Agent - For spatial epigenomics analysis150* Cell_Cell_Communication - For ligand-receptor analysis151* Nicheformer_Spatial_Agent - For spatial niche modeling152153## Special Considerations1541551. **Batch Effects**: Pre-align datasets from different protocols1562. **Spot Deconvolution**: Lower resolution platforms need deconvolution1573. **Sparsity**: scATAC data requires aggregation strategies1584. **Compute**: Multi-modal integration is memory-intensive1595. **Validation**: Verify spatial patterns with known marker distributions160161## Applications162163| Application | Use Case |164|-------------|----------|165| Tumor Microenvironment | Map chromatin states of immune infiltrates |166| Development | Track lineage chromatin dynamics spatially |167| Neurodegeneration | Spatial mapping of epigenetic changes |168| Fibrosis | Understand spatial activation programs |169170## Author171172AI Group - Biomedical AI Platform173174175<!-- AUTHOR_SIGNATURE: 9a7f3c2e-MD-BABU-MIA-2026-MSSM-SECURE -->