STRAND: Sequence‑conditioned transport for single‑cell perturbations

Boyang Fu1, Sameer Gabbita2, George Dasoulas3, Xiang Lin1, Shanghua Gao1, Xiaorui Su1, Soumya Ghosh3, Marinka Zitnik1,4,5

1Harvard Medical School  2Johns Hopkins University  3Merck & Co., Inc.  4Kempner Institute, Harvard University  5Broad Institute of MIT and Harvard
NeurIPS 2026

Perturbations of the same gene can produce different transcriptional responses depending on the targeted locus, such as an alternative transcription start site or a distal enhancer, yet gene-level models represent these interventions with a single gene identifier. STRAND encodes the regulatory DNA sequence at the target together with a mask over the perturbed positions, and uses this encoding to transport observed control cells to their predicted perturbed state, at 128 bp model resolution within a 524 kb window.

The animation shows the four steps of one prediction: STRAND reads a 524,288 bp window of reference DNA in 128 bp bins, masks the bins at the target locus, encodes the masked sequence as a perturbation condition, and transports a control cell to its predicted perturbed state. In the expression readout, green marks an increase and red a decrease relative to the control cell.

Predicted responses differ by perturbed locus and by cell line

Each panel below reports in silico perturbations, in which STRAND masks a window of the reference sequence and predicts the resulting single-cell state.

STRAND attains the best average rank on few-shot and unseen-gene benchmarks

The benchmarks compare STRAND with gene-level perturbation models on CRISPRi Perturb-seq screens in K562, Jurkat and RPE1 cells. The three tabs correspond to few-shot prediction with five cells of each test gene in training, prediction for genes held out from training, and transfer to a cell line held out from training.

How does STRAND compare with gene-level baselines?

STRAND predicts perturbed expression by transporting control cells, conditioned on the DNA sequence at the target

STRAND combines two frozen pretrained encoders with three trainable modules. The Borzoi DNA encoder2, in its enhanced Flashzoi version3, embeds the 524 kb sequence window in 128 bp bins, and the perturbation module combines these embeddings with the mask into a perturbation condition. The SCimilarity RNA encoder4 places each control cell in a latent space of cell states, where an I²SB bridge5, a learned stochastic path from control to perturbed states, transports the cell according to the perturbation condition. A decoder then predicts the expression of each gene from the transported state.

How does STRAND work?

What does each component contribute?

Each variant is retrained with one module replaced and all other components and losses held fixed. In the first variant, the perturbation module is replaced by DNA embeddings that are averaged over the window and kept fixed; in the second, the latent transport module is replaced by a deterministic MLP predictor that takes the control and DNA embeddings as input. Values are means over test perturbations, and average rank is the mean rank of each variant across the seven metrics (Table S10, p.27 of the preprint).

Future Direction: STRAND-guided Prioritization of Candidate Perturbation Sites

We envision combining STRAND's predicted perturbation effects with published experimental findings and genomic annotations to prioritize candidate regions for follow-up experiments.

  1. Define the queryThe user specifies the target gene, the cell line and the genomic region to scan.
  2. Run a STRAND scanSTRAND predicts the perturbation effect of each fixed-width window as the window moves across the selected region.
  3. Shortlist candidate regionsThe scan groups neighboring high-scoring windows into regions and keeps those with the largest predicted perturbation effects.
  4. Retrieve evidenceFor each candidate region, experimental studies, regulatory annotations and results from earlier perturbation screens are collected.
  5. Prioritize candidatesAn expert evaluates STRAND predictions alongside the available biological evidence and ranks candidate regions for further investigation.

Paper and code

Paper

STRAND: Sequence-Conditioned Transport for Single-Cell Perturbations. NeurIPS 2026.

Preprint, arXiv:2602.10156

Code

The repository holds the code to train STRAND, reproduce the benchmarks, and scan a locus.

GitHub, mims-harvard/STRAND

Data

STRAND is trained and evaluated on single-cell CRISPRi Perturb-seq screens, curated in PerturbQA: K562 and RPE1 from Replogle et al. and Jurkat from Nadig et al.

Replogle et al., Cell 2022
Nadig et al. 2025, GEO GSE264667

Citation

@inproceedings{fu2026strand,
  title     = {{STRAND}: Sequence-Conditioned Transport for Single-Cell Perturbations},
  author    = {Fu, Boyang and Gabbita, Sameer and Dasoulas, George and Lin, Xiang and
               Gao, Shanghua and Su, Xiaorui and Ghosh, Soumya and Zitnik, Marinka},
  booktitle = {Advances in Neural Information Processing Systems},
  year      = {2026}
}

Correspondence: Marinka Zitnik, marinka@hms.harvard.edu

References

  1. Chiabrando et al. The mitochondrial heme exporter FLVCR1b mediates erythroid differentiation. J Clin Invest 122, 4569–4579 (2012). doi:10.1172/JCI62422; PubMed 23187127
  2. Linder et al. Predicting RNA-seq coverage from DNA sequence as a unifying model of gene regulation. Nat Genet 57, 949–961 (2025). doi:10.1038/s41588-024-02053-6
  3. Hingerl, Karollus and Gagneur. Flashzoi: an enhanced Borzoi for accelerated genomic analysis. Bioinformatics 41, btaf467 (2025). doi:10.1093/bioinformatics/btaf467
  4. Heimberg et al. A cell atlas foundation model for scalable search of similar human cells. Nature 638, 1085–1094 (2025). doi:10.1038/s41586-024-08411-y
  5. Liu et al. I²SB: Image-to-Image Schrödinger Bridge. ICML (2023). arXiv:2302.05872