Upstream regulation of neurotransmitter signaling genes in metastatic colorectal cancer (mCRC) patients: Data from CALGB/SWOG 80405 trial (Alliance) and MONSTAR-SCREEN-2.
Authors: Millstein Joshua , Battaglin Francesca , Yang Yan , Soni Shivani , Algaze Sandra , Ashouri Karam , Mittal Pooja , Torres-Gonzalez Lesly , Zhang Wu , Ou Fang-Shu , Hochster Howard S. , Venook Alan P. , Innocenti Federico , Mumenthaler Shannon M. , Hashimoto Tadayoshi , Fujisawa Takao , Nakamura Yoshiaki , Bando Hideaki , Yoshino Takayuki , Lenz Heinz-Josef 
Published: June 2024

3523

Background: Evidence suggests neurotransmitter signaling may play a crucial role in tumorigenesis of gastrointestinal (GI) cancers, including migration, proliferation, and apoptosis, as well as angiogenesis and the microenvironment. We developed and applied a strategy to identify upstream regulators of neurotransmitters and their associations with patient outcomes. Methods: Prior literature yielded 85 neurotransmitter signaling genes mapped to RNA-seq gene expression in tumor tissue from 545 patients with mCRC enrolled in the phase III CALGB/SWOG 80405 trial. These samples were also sequenced by FoundationOne CDx for DNA damage. Using Cox proportional hazards regression, we identified a subset of CDx genes associated with overall survival (OS) or progression-free survival (PFS). We extended the Causal Inference Test (CIT) for time-to-event data and identified a subset of the 85 signaling genes consistent with mediation of DNA mutation effects on OS and PFS. Another CIT analysis identified potential upstream neurotransmitter regulators, gene expression mediators of mutation effects on the signaling genes. Top genes were then assessed in whole-transcriptome sequencing (WTS) expression data with Caris MI profile collected from 398 patients with advanced CRC enrolled in the MONSTAR-SCREEN-2 trial. Results: At the 0.05 alpha level, we identified 20 signaling genes with gene expression associated with PFS and/or OS, and 68 genes with evidence of DNA mutations associated with PFS and/or OS. False discovery rate (FDR) < 0.05 was applied to filter the 56,674 expression features on PFS and OS associations, yielding 173 unique genes as possible mediators. The CIT yielded 53 putative mediators at FDR < 0.05, suggesting upstream regulation of neurotransmitter signaling. Five genes were identified as targeting at least 8 signaling genes each, PIGU, ANO9, LY6G6D, and GPR160, and FSCN1, where the first 4 are favorable and only FSCN1 is a poorer prognostic. Following up on these 5 genes in MONSTAR-SCREEN-2 by dichotomizing at the median value of WTS expression, 3 were nominally associated with OS, PIGU ( P = 0.050; HR, 0.59; 95%CI, 0.34-1.00), ANO9( P = 0.029; HR, 0.55; 95%CI, 0.32-0.95), LY6G6D( P = 0.0024; HR, 0.44; 95%CI, 0.25-0.76), and FSCN1 was suggestive ( P = 0.079; HR, 1.62; 95%CI, 0.94-2.77), with direction of effects in agreement with 80405 results for all four. Conclusions: Our novel approach identified previously unknown upstream regulators of neurotransmitter signaling associated with treatment outcomes in mCRC. PIGU, FSCN1, and LY6G6D have been suggested in prior literature as targets for CRC and ANO9 has been suggested as a target for GI cancer. These results bolster the strategy of using mediation analysis with large-scale molecular data to identify new therapeutic targets in CRC. NCT00265850.