Dynamic analysis of gastruloids on microraft array (DAGOMR)
Data files
Nov 11, 2026 version files 20.77 GB
-
DAGOMR_Images.zip
20.74 GB
-
DAGOMR.zip
36.59 MB
-
README.md
7.20 KB
Abstract
Two-dimensional human gastruloids are a reproducible and scalable model system to study aspects of early embryonic development. However, the heterogeneity of single gastruloids, which consist of multiple cells that self-organize into different tissue layers, cultured under homogenous conditions have yet to be quantified. The microraft array platform was used to capture time-series imaging of SOX2 dynamics in developing gastruloids and to isolate individual gastruloids for downstream RNA-sequencing. Low-dimensional latent representations were constructed from the imaging time-series of each gastruloid to quantify the emergence of various patterning behaviors across and within experimental batches. The heterogenous gene expression profiles of sequenced gastruloids were also dissected by experimental batch. Imaging and sequencing appeared to provide orthogonal information, thus more robust computational approaches are likely required to decipher the complex relationship between the data modalities. The heterogeneity of gastruloids, in addition to other multicellular models, calls for careful experimental design considerations.
Dataset DOI: 10.5061/dryad.1vhhmgr80
Description of the data and file structure
This repository hosts datasets describing metrics, time-series imaging-based latent representations, and gene expression derived from two-dimensional human gastruloids for the following manuscript in APL Bioengineering: Heterogeneous Patterning Behavior Emerges from Human Stem Cell-Based Embryo Model.
Files and variables
File: DAGOMR.zip
Endpoint Gastruloid Analysis
Features manually extracted from high-quality gastruloids at the final time point.
Batch: experimental batch ID
Raft: microraft index
Channel: full (entire gastruloid) or SOX2 (SOX2+ region)
Feature: radius (radius in µm) or eccen (eccentricity)
Value: value of corresponding feature
Batch: experimental batch ID
Channel: full (entire gastruloid) or SOX2 (SOX2+ region)
Feature: radius (radius in µm) or eccen (eccentricity)
Mean: mean value of corresponding feature
Stdev: standard deviation of corresponding feature
Endpoint Gastruloid Grading
Quality grading of 50 fully developed gastruloids at the final time point.
Batch: experimental batch ID
Raft: microraft index
Grade: microraft grade [centered (centered gastruloid on microraft), defect (presence of microraft defect), aggregate (presence of adjacent aggregate), border (gastruloid contacting microraft border)]
SOX2_Center: whether gastruloid exhibits central circular SOX2+ region (yes or no)
Gastruloid Metadata
Metadata of gastruloid samples isolated for sequencing.
Sample_Name: sequencing sample ID
Plate_Date: date of collection plate for sequencing
Well: well address
Array_Date: date of array (initial culture)
Batch: experimental batch ID
Raft: microraft index
Collect: whether gastruloid was successfully collected (0/1 for failure/success)
RNA: whether gastruloid was successfully sequenced (0/1 for failure/success)
Weight to Image Loss
16-dimension (16D) latent representations created from individual YX slices from time-series imaging of SOX2-mCitrine reporter in developing gastruloids while varying weight to image loss.
dim_{n}: latent feature n value
timepoint: time point [0, 8]
sample_id: {array date}_array1_{raft number}
split: training (train), validation (val), or test (test) dataset
sample_id: {array date}_array1_{raft number}
timepoint: time point [0, 8]
PC{n}: principal component n value
Latent Representations
16-dimensional (16D) latent representations created from individual YX slices (static) or TYX stacks (dynamic) from time-series imaging of SOX2-mCitrine reporter in developing gastruloids.
emb_{n}: latent feature n value
sample_id: {array date}_array1_{raft number}
split: training (train), validation (val), or test (test) dataset
batch: experimental batch ID
raft: microraft index
umap_1: UMAP1 values
umap_2: UMAP2 values
emb_{n}: latent feature n value
sample_id: {array date}_array1_{raft number}
timepoint: time point [0, 8]
split: training (train), validation (val), or test (test) dataset
batch: experimental batch ID
raft: microraft index
timepoint: time point [0, 8]
umap_1: UMAP1 values
umap_2: UMAP2 values
Model Loss as Function of Latent Dimensions
Model loss (image mean squared error and time point cross entropy) as a function of number of latent dimensions.
latent_dim: number of latent dimensions
image_loss_weight: weight to image loss [0, 1]
timepoint_loss_weight: weight to time point loss [0, 1]
best_test_image_loss: lowest image loss
best_test_image_epoch: minimum epochs to achieve lowest image loss
best_test_timepoint_loss: lowest time point loss
best_test_timepoint_epoch: minimum epochs to achieve lowest time point loss
Model Replicates
16-dimensional (16D) latent representations created from individual YX slices (static) from time-series imaging of SOX2-mCitrine reporter in developing gastruloids across multiple model runs/seeds.
seed: model run seed [0, 4]
sample_id: {array date}_array1_{raft number}
timepoint: time point [0, 8]
dim_n: latent feature n
dim_n: latent feature n
timepoint: time point [0, 8]
sample_id: {array date}_array1_{raft number}
split: training (train), validation (val), or test (test) dataset
sample_id: {array date}_array1_{raft number}
timepoint: time point [0, 8]
umap_1: UMAP1 values
umap_2: UMAP2 values
Bulk RNA-seq
Sequencing data from single euploid and mosaic gastruloids.
Gene: gene name
IJ.{plate date}.{well address}: gene expression count for sequencing sample ID
User_ID: gene name
ensembl_ID: gene name
symbol: gene name
search_label: gene name
IJ{plate date}{well address}: normalized gene expression count for sequencing sample ID
Cluster: gene expression cluster ID
Sample: sequencing sample ID
User_ID: gene name
ensembl_ID: ENSEMBL gene ID
symbol: gene name
search_label: {gene name} | {ENSEMBL gene ID}
IJ{plate date}{well address}: gene expression count normalized by standard deviation for sequencing sample ID
Sample_ID: sample ID
PC{n}: principal component n value
Batch: experimental batch ID
Raft: microraft index
Experimental batch ID corresponding to each sequencing sample ID.
File: DAGOMR_Images.zip
Processed
Processed imaging time-series slices (YX) [{array date}_green_array1_{microraft index}_{time point}.png] of SOX2-mCitrine reporter in living gastruloids for each time point. Subdirectory names correspond to experimental batch ID.
Raw
Raw imaging time-series stacks (TCYX) [array1_raft{microraft index}_tl.tif}] of living gastruloids. Subdirectory “0623_Masks” contains manually drawn masks for the whole gastruloid (roi) and SOX2+ region (sox). Subdirectory names correspond to experimental batch ID.
BF: brightfield (transmitted light)
SOX2: SOX2-mCitrine
Code/software
For code used, please refer to doi.org/10.5281/zenodo.17073719.
