# CTAS: a network control theory-based approach to identify key regulatory TFs of AS events during epithelial–mesenchymal transition

## TL;DR

This paper introduces CTAS, a new method to identify key transcription factors that control alternative splicing during epithelial–mesenchymal transition.

## Contribution

CTAS uses network control theory to uncover multi-layered regulatory logic from bulk data, enabling the identification of TFs controlling EMT-related AS events.

## Key findings

- CTAS reconstructs EMT trajectories with high accuracy (Spearman’s ρ = 0.99946) and infers directed networks with 89.9% ROC AUC.
- In TCGA BRCA data, CTAS identifies HOXA3, PRDM8, and TWIST2 as top TF controllers of AS events during EMT.
- Dynamic shifts in nine AS events were detected, with ZNF521 and HIC1 highlighted as candidate regulators in a CD44 subnetwork.

## Abstract

Alternative splicing (AS) is a key driver of transcriptomic diversity and plays a pivotal role in epithelial–mesenchymal transition (EMT). During EMT, dynamic splicing changes contribute to cell plasticity and metastasis, yet the upstream regulatory logic remains unclear. Although transcription factors (TFs) are thought to influence AS programs, they typically act through RNA-binding proteins (RBPs), forming a hierarchical TF\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{upgreek}
\usepackage{mathrsfs}
\setlength{\oddsidemargin}{-69pt}
\begin{document}
$\rightarrow $\end{document}RBP\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{upgreek}
\usepackage{mathrsfs}
\setlength{\oddsidemargin}{-69pt}
\begin{document}
$\rightarrow $\end{document}AS cascade. Current computational strategies struggle to recover such multi-layered regulation from bulk cross-sectional data, limiting our ability to identify TFs that ultimately control EMT-related AS events. To address this gap, we developed CTAS, a network control theory-based approach to identify key regulatory TFs of AS events during EMT. CTAS integrates pseudotime ordering, trend analysis, sparse directed network inference, and control-theoretic screening into a unified framework. In simulations, CTAS reconstructs EMT trajectories with Spearman’s \documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{upgreek}
\usepackage{mathrsfs}
\setlength{\oddsidemargin}{-69pt}
\begin{document}
$\rho = 0.99946$\end{document} and directed networks with ROC AUC = 89.9%, and remains robust under noise. Applied to a TCGA BRCA cohort, CTAS builds a directed TF\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{upgreek}
\usepackage{mathrsfs}
\setlength{\oddsidemargin}{-69pt}
\begin{document}
$\to $\end{document}RBP\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{upgreek}
\usepackage{mathrsfs}
\setlength{\oddsidemargin}{-69pt}
\begin{document}
$\to $\end{document}AS network and identifies HOXA3 (1.00), PRDM8 (0.86), and TWIST2 (0.83) as top TF controllers, alongside significant dynamic shifts in nine AS events detected by Wilcoxon test (\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{upgreek}
\usepackage{mathrsfs}
\setlength{\oddsidemargin}{-69pt}
\begin{document}
$P <.05$\end{document}). A focused CD44 subnetwork further highlights ZNF521 (0.86) and HIC1 (0.65) as candidate regulators. These findings demonstrate that CTAS transforms cross-sectional data into dynamic regulatory insights and yields experimentally testable TFs that control AS during EMT.