Enabling Spatially Resolved Gene Expression Analysis in a Multilayer Paper-Based In Vitro Tumor Model
- 1Indian Institute of Science
Abstract
Spatial heterogeneity is a defining feature of solid tumors, yet molecular analysis of three-dimensional (3D) in vitro tumor models is often performed as a bulk measurement, potentially obscuring position-dependent cellular responses. Here, we establish an analytical workflow for spatially resolved gene-expression analysis in a multilayer paper-based 3D tumor model by addressing the challenge of RNA recovery from cell-laden paper-Matrigel scaffolds. Using MDA-MB-231 breast cancer cells, multiple RNA extraction strategies were systematically evaluated across conventional culture, Matrigel, paper, and composite paper-Matrigel environments. Conventional silica-column extraction yielded only 3.39 {+/-} 0.93 ng/L RNA from the composite scaffold, whereas TRIzol-chloroform extraction increased recovery to 201.13 {+/-} 38.32 ng/L. The optimized workflow subsequently enabled independent RNA recovery from all eight layers of a multi-layer paper-based in vitro tumor, with concentrations ranging from 113.06 {+/-} 5.00 to 245.26 {+/-} 15.76 ng/L and purity suitable for downstream qRT-PCR analysis. Layer-resolved profiling of untreated cultures revealed distinct spatial patterns in HIF1A, HIF2A, VEGFA, and GLUT4 expression across the multilayer construct. To evaluate the workflow under pharmacological perturbation, cultures were exposed to 5 M doxorubicin (DOX) for 48 h. Relative cell viability exhibited a non-monotonic spatial profile, whereas DOX-induced HIF1A and VEGFA expression generally became more pronounced toward the deeper layers when each treated layer was compared with its corresponding untreated control. Collectively, these results establish a workflow for converting individually addressable regions of a multilayer 3D tumor model into spatially resolved transcriptional measurements, enabling analysis of both baseline spatial heterogeneity and position-dependent molecular responses to pharmacological perturbation.
Related
- Spatially resolved gene expression analysis illuminates location-specific functions in the reef-building coral Pocillopora acuta
- Multilayer MEG source modelling enables depth-resolved laminar inference in humans
- BundleMorphometry: Spatially Resolved White Matter Bundle Shape Analysis Using Nonlinear Streamline-based Registration
- PathwayFinder enables linked exploration of cell types and gene modules in spatially resolved transcriptomics data
- Spatially resolved host-fungal transcriptomics profiling by fungal-DBiT
BibTeX
@misc{kumar2026enabling,
title = {{Enabling Spatially Resolved Gene Expression Analysis in a Multilayer Paper-Based In Vitro Tumor Model}},
author = {A. Kumar and B. J. Toley},
howpublished = {bioRxiv},
year = {2026},
month = oct,
doi = {10.64898/2026.10.18.759958},
url = {https://www.biorxiv.org/content/10.64898/2026.10.18.759958v1},
note = {CC BY-NC-ND},
}
RIS
TY - UNPB
TI - Enabling Spatially Resolved Gene Expression Analysis in a Multilayer Paper-Based In Vitro Tumor Model
AU - A. Kumar
AU - B. J. Toley
PY - 2026
DA - 2026/10/18
PB - bioRxiv
DO - 10.64898/2026.10.18.759958
UR - https://www.biorxiv.org/content/10.64898/2026.10.18.759958v1
AB - Spatial heterogeneity is a defining feature of solid tumors, yet molecular analysis of three-dimensional (3D) in vitro tumor models is often performed as a bulk measurement, potentially obscuring position-dependent cellular responses. Here, we establish an analytical workflow for spatially resolved gene-expression analysis in a multilayer paper-based 3D tumor model by addressing the challenge of RNA recovery from cell-laden paper-Matrigel scaffolds. Using MDA-MB-231 breast cancer cells, multiple RNA extraction strategies were systematically evaluated across conventional culture, Matrigel, paper, and composite paper-Matrigel environments. Conventional silica-column extraction yielded only 3.39 {+/-} 0.93 ng/L RNA from the composite scaffold, whereas TRIzol-chloroform extraction increased recovery to 201.13 {+/-} 38.32 ng/L. The optimized workflow subsequently enabled independent RNA recovery from all eight layers of a multi-layer paper-based in vitro tumor, with concentrations ranging from 113.06 {+/-} 5.00 to 245.26 {+/-} 15.76 ng/L and purity suitable for downstream qRT-PCR analysis. Layer-resolved profiling of untreated cultures revealed distinct spatial patterns in HIF1A, HIF2A, VEGFA, and GLUT4 expression across the multilayer construct. To evaluate the workflow under pharmacological perturbation, cultures were exposed to 5 M doxorubicin (DOX) for 48 h. Relative cell viability exhibited a non-monotonic spatial profile, whereas DOX-induced HIF1A and VEGFA expression generally became more pronounced toward the deeper layers when each treated layer was compared with its corresponding untreated control. Collectively, these results establish a workflow for converting individually addressable regions of a multilayer 3D tumor model into spatially resolved transcriptional measurements, enabling analysis of both baseline spatial heterogeneity and position-dependent molecular responses to pharmacological perturbation.
ER -