Interpretable transcriptome-to-phenotype modeling of cell-painting nuclear morphology features from RNA-seq under low-dose radiation exposure

With rapid advancements in high-throughput multi-modal profiling techniques across molecular, cellular to tissue scales, translating such multi-modal data into knowledge discovery for foundational understanding of cellular mechanisms is central in modern biomedical sciences. In this study, we focus on understanding how low-dose radiation exposure perturbs cellular morphology by linking high-dimensional transcriptomic responses to quantitative cell-phenotype readouts over time. We present a time-resolved inverse modeling framework that associates gene expression changes via RNA-sequencing with nuclear morphology features obtained from cell-painting imaging. Morphology responses were defined as treated-control differences for multiple nuclear features including size, shape, intensity, and textures, indexed by radiation dose and week. To capture time-dependent associations while maintaining interpretability, both RNA-sequencing and cell-painting data were stratified into four temporal phases (weeks 1-2, 3-4, 5-6, 7-9) and phase-dependent effects are encoded via gene-phase interaction predictors. To reduce confounding by dose trends and to evaluate generalization across time, we used a two-stage leave-one-week-out procedure: (i) a dose-only baseline model produced out-of-week residuals for each morphology feature, and (ii) elastic-net regression on phase-aware predictors modeled residual variation not explained by dose. Hyperparameters were selected via an exhaustive grid search scored by the correlation between observed residuals and out-of-week residual predictions, with additional sparsity diagnostics based on nonzero coefficient counts per fold. Stable predictors were identified by selection frequency and sign consistency across folds, then pruned further for multicollinearity and parsimony. Final reduced models were fit using ordinary least squares with heteroskedasticity-consistent standard errors to report effect estimates robust to non-constant variance. This workflow yields a transparent, time-stratified set of transcriptomic predictors associated with longitudinal nuclear morphology changes and provides a reproducible foundation for downstream biological interpretation and validation.