Learning from weak, proxy, or relative supervision is common when ground-truth labels are unavailable, but robustness under distribution shift remains poorly understood because the supervision mechanism itself may change across environments. We formalize this phenomenon as supervision drift, defined as changes in P(y∣x,c) across contexts, and study it in CRISPR-Cas13d transcriptomic perturbation experiments where guide efficacy is inferred indirectly from RNA-seq responses. Using publicly available data spanning two human cell lines and multiple post-induction timepoints, we construct a controlled non-IID benchmark with explicit domain (cell line) and temporal shifts, while reusing a fixed weak-label construction across all contexts to avoid changing targets. Across linear and tree-based models, weak supervision supports meaningful learning in-domain (ridge R2=0.356, Spearman ρ=0.442) and partial cross-cell-line transfer (ρ≈0.40). In contrast, temporal transfer collapses across all model classes considered, yielding negative R2 and weak or near-zero ρ (ridge R2=−0.145, ρ=0.008; XGBoost R2=−0.155, ρ=0.056; random forest R2=−0.322, ρ=0.139). Additional robustness analyses using externally recomputed weak labels, shift-score quantification, and simple mitigation baselines preserve the same qualitative pattern. Feature-label association and feature-importance analyses remain relatively stable across cell lines but change sharply over time, indicating that failures arise from supervision drift rather than model capacity or simple covariate shift. These results show that strong in-domain performance under weak supervision can be misleading and motivate feature stability as a lightweight diagnostic for non-transferability before deployment.
Figures & tables
Figure 1: Overview of the experimental framework for studying generalization under weak supervision. Transcript-level features are extracted from multiple biological contexts (K562 Day 2, HEK293FT Day 2, and HEK293FT Day 7). A single fixed weak supervision signal y is constructed once from the HEK293FT Day 7 minus Day 2 transcriptomic response and reused across all experiments, ensuring that the target remains unchanged across contexts. Models are trained within a single source context (e.g., HEK293FT Day 2 or K562 Day 2) to predict this fixed weak label, without access to other contexts during training. Generalization is then evaluated under three settings: (i) in-domain (HEK293FT Day 2 → HEK293FT Day 2), (ii) cross-domain (K562 Day 2 → HEK293FT Day 2), and (iii) temporal (HEK293FT Day 2 → HEK293FT Day 7). By holding the supervision signal fixed and varying only the feature context, this design isolates changes in the feature–supervision relationship. While partial transfer is observed across cell lines, temporal generalization fails, revealing supervision drift, where the conditional relationship P(y∣x,c) changes across contexts.
Setting
Training features
Test features
Target
In-domain
HEK293FT Day 2
HEK293FT Day 2
fixed y
Cross-domain
K562 Day 2
HEK293FT Day 2
fixed y
Temporal
HEK293FT Day 2
HEK293FT Day 7
fixed y
Table 1: Data assignments for the three evaluation settings. The weak supervision target y is the same fixed transcript-level vector in all cases.
Figure 2: In-domain predictability under weak supervision (HEK293FT, Day 2). Performance is reported using R2 , mean absolute error (MAE; lower is better), and Spearman rank correlation ρ . Ridge regression and a depth-capped random forest achieve comparable in-domain performance, indicating that transcript-level features contain recoverable relative signal under indirect supervision.
Figure 3: Cross-domain generalization from K562 to HEK293FT at Day 2. Performance is reported using R2 (solid bars) and Spearman rank correlation ρ (hatched bars). Both ridge regression ( R2=0.331 , ρ=0.399 ) and a depth-capped random forest ( R2=0.310 , ρ=0.364 ) retain a substantial portion of their in-domain predictive signal, indicating partial cross-cell-line transferability under weak supervision.
Figure 4: Temporal generalization heatmap from Day 2 to Day 7 in HEK293FT. Performance is reported using R2 and Spearman rank correlation ρ (higher is better). All model families considered, including ridge regression, random forest (depth-capped), and XGBoost, exhibit severe performance degradation, indicating consistent failure to transfer predictive signals across time under weak supervision.
Context
Feature
Spearman Corr.
n
K562_D2
logTPM
−0.282
251,955
K562_D2
rank_pct_within_sample
−0.282
251,955
HEK293FT_D2
logTPM
−0.468
251,955
HEK293FT_D2
rank_pct_within_sample
−0.468
251,955
HEK293FT_D7
logTPM
−0.001
251,955
HEK293FT_D7
rank_pct_within_sample
−0.001
251,955
Table 2: Feature–label Spearman correlations across biological contexts. Correlations are computed between transcript-level features (logTPM and within-sample rank percentiles) and weak supervision signals in each context. Feature relevance is consistent across cell lines at Day 2 but collapses at Day 7, indicating a breakdown in the feature–supervision relationship under temporal shift.
Feature
Context Pair
∣Δ Corr. ∣
logTPM
K562_D2 vs HEK293FT_D2
0.186
logTPM
HEK293FT_D2 vs HEK293FT_D7
0.468
rank_pct_within_sample
K562_D2 vs HEK293FT_D2
0.186
rank_pct_within_sample
HEK293FT_D2 vs HEK293FT_D7
0.468
Table 3: Feature stability across structured distribution shifts, measured as the absolute difference in feature–label Spearman correlation between context pairs. Feature relevance remains relatively stable under cross-domain shift (K562 Day 2 vs. HEK293FT Day 2) but exhibits substantially larger instability under temporal shift (HEK293FT Day 2 vs. Day 7), suggesting supervision drift rather than feature-specific effects.
Setting
Model
R2
Spearman ρ
In-domain (HEK D2)
Ridge
0.051
0.205
In-domain (HEK D2)
RandomForest
0.081
0.131
Cross-domain (K562 → HEK)
Ridge
-0.227
-0.205
Cross-domain (K562 → HEK)
RandomForest
-0.363
-0.033
Temporal (D2 → D7)
Ridge
0.006
0.064
Temporal (D2 → D7)
RandomForest
-0.118
-0.029
Table 4: Performance using externally recomputed weak labels. The qualitative behavior remains unchanged: positive in-domain performance, substantial degradation under cross-domain transfer, and collapse under temporal transfer.
Model
ρ(ShiftScore,Performance)
Context Pairs
Ridge
-0.239
6
RandomForest
-0.478
6
Table 5: Correlation between supervision drift magnitude and predictive performance. Larger shifts in feature–label association correspond to worse generalization.
Setting
Model Variant
R2
Spearman ρ
In-domain
Plain
0.359
0.444
In-domain
Context-aware
0.359
0.444
Cross-domain
Plain
0.331
0.399
Cross-domain
Context-aware
0.331
0.399
Temporal
Plain
-0.145
0.008
Temporal
Context-aware
-0.145
0.008
Table 6: Context-aware modeling does not improve temporal generalization.
Figure 5: Feature ablation study across evaluation scenarios. Performance is shown for expression-based features (logTPM only), rank-based features (rank only), and their combination (both), evaluated under in-domain, cross-domain, and temporal settings. Left: coefficient of determination ( R2 ). Middle: mean absolute error (MAE; lower is better). Right: Spearman rank correlation ( ρ ). For comparisons between logTPM-only and rank-only feature sets, the most informative metrics are R2 and MAE, since Spearman correlation is invariant to monotone transformations within a context. Across all feature representations, temporal generalization fails, indicating that supervision instability, not feature choice alone, limits transfer over time.
Appendix figures & tables6 assets
Supplementary material from the paper’s appendix.
Appendix
Evaluation Setting
Features
R2
MAE
Spearman ρ
In-domain (HEK293FT D2)
logTPM
0.377
0.171
0.386
Cross-domain (K562 → HEK293FT)
logTPM
0.340
0.170
0.393
Temporal (HEK293FT D2 → D7)
logTPM
-0.155
0.225
0.056
Appendix
Table 7: XGBoost baseline results under weak supervision. Despite improved in-domain and cross-domain performance relative to linear models, temporal generalization remains poor, indicating that increased model capacity does not resolve temporal distribution shift.
Setting
Model
R2
Spearman ρ
In-domain (HEK D2)
Ridge
0.051
0.205
In-domain (HEK D2)
RandomForest
0.081
0.131
Cross-domain (K562 → HEK)
Ridge
-0.227
-0.205
Cross-domain (K562 → HEK)
RandomForest
-0.363
-0.033
Temporal (D2 → D7)
Ridge
0.006
0.064
Temporal (D2 → D7)
RandomForest
-0.118
-0.029
Appendix
Table 8: Performance under externally recomputed weak labels. The same qualitative failure pattern persists.
Model
ρ(ShiftScore,Performance)
Context Pairs
Ridge
-0.239
6
RandomForest
-0.478
6
Appendix
Table 9: Correlation between shift score and predictive performance. Larger supervision drift corresponds to worse generalization.
Setting
Model Variant
R2
Spearman ρ
In-domain
Plain
0.359
0.444
In-domain
Context-aware
0.359
0.444
Cross-domain
Plain
0.331
0.399
Cross-domain
Context-aware
0.331
0.399
Temporal
Plain
-0.145
0.008
Temporal
Context-aware
-0.145
0.008
Appendix
Table 10: Context-aware modeling does not improve temporal generalization.
Feature Set
Context Pair
SHAP Rank Corr.
Note
logTPM
K562 D2 → HEK293FT D2
–
Single feature
logTPM
HEK293FT D2 → HEK293FT D7
–
Single feature
Rank pct.
K562 D2 → HEK293FT D2
–
Single feature
Rank pct.
HEK293FT D2 → HEK293FT D7
–
Single feature
Both
K562 D2 → HEK293FT D2
High
Stable ordering
Both
HEK293FT D2 → HEK293FT D7
Moderate
Temporal instability
Appendix
Table 11: SHAP rank correlation diagnostics across contexts. Rank-based agreement is undefined for single-feature models. For the combined feature representation (logTPM + within-sample rank percentile), SHAP enables rank-based comparison, revealing stable importance ordering under cross-domain shift and partial instability under temporal shift.
Sample
Percent Mapped (%)
Mapped Reads
Processed Reads
Notes
HEK293FT D2
87.18
431M
495M
High-quality mapping
K562 D2
86.41
376M
435M
High-quality mapping
HEK293FT D7
71.66
52.6M
73.4M
Moderate mapping
Long-term (SRX26635950)
0.02
39K
169M
Failed mapping
Long-term (SRX27451361)
0.00
0
34M
Failed mapping
Appendix
Table 12: Salmon mapping summary across samples. Long-term samples exhibit near-zero mapping rates and were excluded from downstream modeling, while short-term samples demonstrate consistently high alignment quality.
Engineering Research Center of the Ministry of Education for Intelligent Control System and Intelligent Equipment, Yanshan University, Qinhuangdao, China