Note: This response was posted by the corresponding author to Review Commons. The content has not been altered except for formatting.
Learn more at Review Commons
Reply to the reviewers
We thank the three reviewers for their thoughtful and constructive comments which help us to improve the manuscript. Please find our responses below. * *
Reviewer #1 (Evidence, reproducibility and clarity (Required)):
Summary This study investigates how altered expression of cleavage and polyadenylation (CPA) factors affects alternative polyadenylation (APA), transcription termination and cellular phenotypes in colorectal cancer (CRC) cell lines. The authors combine genetic perturbations of CPA factors with chemical inhibition of CPSF73 and assess effects on clonogenic potential, transcription-replication conflicts, APA profiles, and transcription termination-associated RNAPII phosphorylation patterns. The main comparisons are performed between healthy (1CT), primary tumour (SW480, HCT116) and metastatic (SW620) cell lines, which are reported to contain altered expression levels of CPA factors. The data suggest differential dependence on CPA factors between primary tumour-derived and metastatic CRC cell lines, as well as changes in transcription termination patterns. The data are overall well-presented with clear figures. However, in several cases the strength of the conclusions appears to exceed the support provided by the data, and alternative interpretations should be considered.
Major comments 1. Clonogenic sensitivity to CPA factor perturbation and comparability of clonogenic assays between cell lines: -The data indicate that clonogenic potential in SW480 is strongly dependent on CPSF73 and PCF11, whereas SW620 appear less sensitive. However, the interpretation is complicated by differences in depletion efficiencies. In SW620 cells, PCF11 depletion appears inefficient, and protein levels remain higher than in siLUC-treated SW480 cells (Fig. 1D and S1C; also in comparison to 1CT by inference of Fig. 1C). Thus, the apparent resistance of SW620 cells could reflect insufficient depletion rather than true biological tolerance. The effectiveness of siCPSF73 treatments is difficult to assess from the presented data. Quantification of protein knockdown levels should be provided and incorporated into the interpretation.
-In Fig. 1D, 1E, and S1D, colony formation of DMSO- or siLUC-treated SW620 and SW480 cells differs markedly in absolute terms. However, the graphs are normalized separately for each cell line, which obscures this difference. This raises two concerns: First, the baseline clonogenic capacity differs between the lines and should be discussed. Second, it is unclear whether direct comparisons between cell lines are valid when normalization is performed independently. For example, in absolute terms, 1 µM JTE-607 appears to have a similar effect in SW620 cells as 5 µM in SW480 cells, which would contradict the conclusion that metastatic cells are more tolerant to CPA perturbations. This issue should be explicitly addressed.
We thank the reviewer for those thoughtful comments.
a) Assessing the biological meaning of differences in PCF11 depletion efficiency between SW480 and SW620 cell lines is inherently tricky, because the two cell lines differ 3-fold in their baseline PCF11 level (Fig. 1C). Even with equal efficiencies of knock-down, the number of PCF11 molecules per cell left after the treatment will differ. We haven't mentioned this in our original manuscript but will highlight this issue in the revised version - as we agree it is an important consideration for the interpretation of the results.
b) As requested, we will add quantification of western blots from 3 biological replicates to the revised manuscript, to demonstrate the depletion efficiencies. We agree that the single western blot presented by itself was not sufficient; the efficiency of SW620 knock-down is not lower compared to SW480.
c) The baseline clonogenic capacity of SW480 and SW620 has been previously calculated and compared in two publications (PMID: 31961892 and 29796953). In both cases, the SW620 cells showed higher clonogenic potential than SW480, which was calculated based on the number of clones containing more than 50 cells.
d) The reason behind normalization of our data to a control sample is the difference in cell size between the cell lines, which prohibits their direct comparison.
For the colony formation assays, we seeded the same number of cells and cultured them for the same amount of time. However, the difference in cell size, leads to a huge difference in colony sizes (Figure 1D), therefore it was not possible to set the same parameters for counting colonies of SW480 and SW620 cells. Therefore, we decided to use an approach frequently used in high profile cancer studies (e.g. Li at al., 2023, PMID: 37620362, Waterhouse et al., 2025, PMID: 40328966, Yang et al., 2026, PMID: 41484364) and normalize each biological replicate to the control sample to analyze the response to the treatment only.
e) During revision, we might additionally perform CellTiter 96® Non-Radioactive Cell Proliferation Assay (MTT) to test how another cancerous characteristic of SW480 and SW620 cells are affected by JTE-607.
f) We will also perform colony formation and/or MTT assays for 3 additional cell lines: HCT116 (primary tumor-derived) and T87 and COLO-205 (metastasis-derived, which we are currently in the process of obtaining) to assess their sensitivity to JTE-607.
g) The result of higher sensitivity of SW620 cells compared to SW480 cells has been obtained not only for PCF11 knock-down, where inter-cell line differences of baseline protein level make interpretations more difficult, but also for CPSF73 knock-down (Fig. 1D), which baseline level was similar and knock-down was equally efficient in both cell lines, and for CPSF73 inhibition (Fig. 1E); with the use of normalization procedures used frequently in literature (see point d).
Therefore, we argue that our conclusion that SW480 cells are more sensitive than SW620 to the abrogation of 3' pre-mRNA cleavage and transcription termination is valid. However, we are willing to weaken our conclusion if the reviewer does not agree with our point of view.
For the additional cancer-specific experiments proposed above, we suggest the usage of JTE-607 as drug treatment is more robust, reproducible, and medically relevant compared to knock-down experiments.
- Interpretation of transcription termination markers: -The study uses RNAPII T4ph as a marker of transcription termination, which is well justified based on the ref. [30], but still the mechanistic basis of this modification is not fully understood. Changes in T4ph localization are interpreted as consequences of CPA activity, but possible differences in kinase or phosphatase activities between cell lines are not considered that could affect the T4ph levels or localization. Therefore, conclusions based solely on T4ph redistribution should be presented with greater caution, and alternative explanations should be acknowledged.
While in our experience RNAPII T4ph is the most sensitive and useful termination marker, we agree with the referee that its metabolism and function is insufficiently understood - this is an important and interesting direction for future investigation.
In order to increase the robustness of our study, during revision we will additionally perform nascent transcriptomics on SW480 and SW620 using a different method, POINT-seq. POINT-seq in contrast to T4ph mNET-seq relies neither on RNAPII modification status nor is affected by pausing. We will also probe global T4ph-RNAPII levels in our cellular model by western blot. We will then adjust our manuscript accordingly.
-Line 240 states that premature termination is increased in primary tumour cells. However, the data show increased T4ph signal (Fig. 4B) but no change in total RNAPII occupancy in gene bodies (Fig. 4A). This does not directly demonstrate increased termination. Additional evidence or a more cautious interpretation would be appropriate.
The reviewer is right in pointing out the difference between the Total-RNAPII and T4ph-RNAPII signals across the gene body. We will provide a clearer description and explanation in the revised manuscript.
T4ph-RNAPII is present at low levels in human cells. S2ph and S5ph are the dominant modifications, accounting for ~75% of phospho-counts, whereas T4ph has a relative abundance of ~15% (PMID: 26799765). In addition, T4ph is concentrated at gene ends and typically very low in the gene body (PMID: 28017589, 30819644, doi: 10.1101/2025.07.14.664659). Consequently, it is very easy to spot its gene-body increase in metagene analysis (Figure 4B), even when it happens only on a subset of genes in cancer samples (e.g. Fig. 4D).
Total-RNAPII signal in the gene body largely reflects S2ph-modified RNAPII levels so its metagene analysis is not sufficiently sensitive to detect differences in gene-body T4ph-RNAPII.
Consequently, RNAPII-T4ph and RNAPII-total mNET-seq show distinct metagene patterns and different responses to termination changes. RNAPII-T4ph mNET-seq is a sensitive method to detect changes in termination patterns, while total-RNAPII is much less specific and sensitive with respect to transcription termination.
- Cleavage-termination distance as a predictor of transcript levels: -Figure 5A presents median distances across all genes. It would be informative to perform a gene-wise comparison between cell lines (difference in cleavage-T4ph distance for the same gene, e.g. in 1CT vs. HCT116, individual differences plotted across all genes). This analysis could help clarify how frequently individual genes experience the effect (shortening of the cleavage-T4ph distance between 1CT and tumour cells) that is observed globally.
Thank you for this valuable suggestion. We have performed the gene-wise comparison which is indeed very informative. Firstly, we observed the same trend as for all active protein-coding genes - shorter distance in all CRC cell lines compared to 1CT cells with the lowest values of the cleavage-termination distance in the primary tumor cells. Secondly, and even more importantly, this analysis additionally shows that the shortening effect is global - only a small percentage of genes do not undergo shortening of the cleavage-T4ph distance between 1CT and tumor cells.
We will incorporate the results of this analysis into the figures of the revised manuscript.
-The manuscript claims that proximity between pre-mRNA 3′-end cleavage and transcription termination predicts increased nuclear transcript levels. However, the correlation coefficients are small (Spearman r ~ -0.2 at most), indicating weak predictive power. Therefore, the use of the term "predicts," especially in the manuscript title, appears to overstate the strength of the relationship. The authors should either moderate this claim or provide additional analysis to support stronger predictive value.
We agree with the reviewer that the term "predicts" is not ideal in this context and are happy to substitute "is associated with". The title would then read: "Proximity of pre-mRNA 3′ end processing and transcription termination is associated with enhanced gene expression".
Minor comments -Figures 1B and S1A: The discontinuous y-axis makes it difficult to assess relative protein level differences between normal and cancer samples. Statistical testing should be included to evaluate significance.
We had decided against statistical testing due to the problems with biological interpretation of such analyses and its limitations for proteins present in the cell at low levels and/or highly variable between samples. PCF11 is such protein. It is an order of magnitude less abundant compared to other RNA 3' processing factors, and its levels are variable as shown in our Fig. 1B (re-analyzed proteomics data from Wiśniewski et al., 2015). Therefore, the increase in PCF11 levels in this dataset is not statistically significant in Mann-Whitney test, while it is significant for CPSF73.
The variability of PCF11 levels can be also observed in the Clinical Proteomic Tumor Analysis Consortium (CPTAC) data in the Human Protein Atlas (while no absolute quantification was performed there).
In two independently obtained proteomics patient datasets (Wiśniewski et al., 2015; CPTAC), as well as in Western blot assays from our cell culture model, an increase in PCF11 protein abundance is observed in cancer cells. This consistency across different datasets and our model holds greater biological relevance than the statistical analysis of highly varied samples. Nevertheless, if the reviewer requires statistics, we will include them in the revised manuscript.
The discontinuous y-axis was applied due to broad range of protein molecules. Presentation of data with linear continuous scale did not allow to present the difference between normal and cancer samples for all the proteins on the same graph.
Alternatively, if the reviewer and editor prefer, we are happy to present the data with log10 transformed scale. The disadvantage of log-scale is that the differences between normal and cancer samples are less obvious to the eye, the advantage is a continuous y-axis.
-Lines 217-218: The text should emphasize that nuclear RNA abundance may not reflect cytoplasmic mRNA levels, particularly when APA alters 3′UTRs and may affect mRNA stability.
We agree and will incorporate the reviewer's suggestion in our revised manuscript.
-Lines 261-264: The cleavage-termination distance metric should be more clearly defined as the distance between the polyadenylation site and the T4ph signal peak.
We plan to incorporate a drawing into the figure, to better explain our cleavage-termination definition.
We also performed the cleavage site to T4ph signal peak (highest signal in the termination window) distance calculations, and they show the same trend as our original method (Figure 5A), with no changes to the conclusions we made. We will incorporate this additional analysis into a supplementary figure.
**Referees cross-commenting**
Reviewer #3:
On the contrary to implied in the reviewer report, this manuscript does not report the effects of CPSF73 inhibitor JTE-607 on APA. On this note, as the authors discuss uncoupling of cleavage and transcription termination, they could consider (this is not a request) testing how the cleavage inhibitor JTE-607 impacts the distribution of transcription termination marker T4ph, and whether the effects would be different in different cell lines where the coupling appears to be different. This could give mechanistic insights into the sources of the differences between cell lines.
In order to get a mechanistic idea why shorter cleavage-termination distance is associated with higher gene expression, we plan to test the cleavage efficiency on genes, which show differences in cleavage-termination distance and expression levels, between SW480 and SW620 cell lines. To this end, we will perform POINT-seq, checking differences between those cell lines in control conditions and with JTE-607. We believe that this new experimental approach will provide a deeper mechanistic insight, compared to performing further correlation analyses repeating the same experiment types.
Reviewer #1 (Significance (Required)):
This study addresses an important question in RNA biology and cancer research: how altered expression or pharmacological targeting of CPA factors affects alternative polyadenylation, transcription termination, and cellular phenotypes in CRC models. This topic is timely, as CPSF73 has been proposed as a therapeutic target, making it important to understand the molecular and cellular consequences of modulating CPA factor activities. A key strength and robust finding of the study is the identification of unexpected relationships between pre-mRNA 3′-end processing and transcription termination during CRC progression. Notably, the authors report that changes in alternative polyadenylation and transcription termination appear to be uncoupled and may even occur in opposite directions. This challenges simplified models in which these processes are tightly coordinated and suggests that their (mis)regulation in cancer cells may be more complex than previously appreciated. Secondly, the study provides an interesting observation that gene-specific changes in cleavage-T4ph distance correlate negatively with changes in nuclear levels of processed transcripts. This suggests a potential relationship between the spatial coupling of 3′-end processing and transcription termination and transcript abundance. If validated mechanistically, this could represent a conceptual advance in understanding how transcription termination dynamics influence gene expression outputs. However, the observed correlations are relatively weak, and the mechanistic basis of this relationship remains unclear. As such, this advance is primarily descriptive at this stage.
As indicated in response to the cross-commenting point above, one possible mechanistic explanation why shorter cleavage-termination distance could be associated with higher gene expression, is increased cleavage efficiency when the cleavage-termination distance is short. To test this hypothesis, we will perform POINT-seq on SW480 and SW620 cell lines, in control and CPSF73 inhibition conditions. We have previously demonstrated that POINT-seq technique allows calculation of cleavage efficiencies, and its alterations (doi: 10.1101/2025.07.14.664659).
So far, our data (Fig. 5F, G) indicates that PCF11 is involved in this process since PCF11 downregulation resulted in lengthening the distance between 3′-end cleavage and RNAPII terminal pausing. This lengthening was in parallel correlated with the decrease of the nuclear RNA levels. However, PCF11 participates in multiple steps of gene expression - pre-mRNA cleavage, alternative polyadenylation, RNAPII pausing, and mRNA export - making the underlying mechanism difficult to pinpoint without additional experiments.
Importantly, our work provides the first clear evidence that changes in cleavage site usage and termination region usage can become uncoupled. We hope that continued tool development, together with studies like ours, will ultimately enable a full mechanistic understanding.
Several interpretations of experimental data would benefit from more cautious framing or additional analysis. In particular, the relationship between changes in CPA factor expression levels and sensitivity to the CPSF73 inhibitor JTE-607 across CRC cell lines remains unclear from the presented data.
During the revision we will explain more clearly the rationale for our interpretation of the data. In cases where more cautious framing would still be needed, we will include alternative interpretations.
This work will be of interest primarily to basic researchers in RNA processing and transcription regulation, gene expression control, cancer cell biology and pharmacological targeting of RNA-processing machineries.
Reviewer field of expertise: My expertise is in RNA processing and gene regulation. I do not have specific expertise clinical oncology or cancer biology.
Reviewer #2 (Evidence, reproducibility and clarity (Required)):
Factors involved in pre-mRNA cleavage and polyadenylation (CPA) are upregulated in many cancers and have been found to be associated with poor prognosis. In their manuscript "Proximity of pre-mRNA 3′ end processing and transcription termination predicts enhanced gene expression", Stepien et al. use colorectal cancer (CRC)-derived cell lines as a model of CPA overexpression to study its biological consequences. To this end, the authors initially confirm increased expression of CPA factors in these cell lines and demonstrate that their knock-down strongly decreases the colony-forming ability of primary tumour-derived CRC cells. They further assess various phenotypes that are expected to depend on CPA activity based on the current knowledge in the field, including poly-A site selection, occurrence of transcription-replication conflicts, and the site of transcription termination. Contrary to expectations, they find a proximal shift in transcription termination to be the most prominent change in CRC ell lines with high CPA levels, despite no clear preference for proximal poly-A site usage in these cells, suggesting an uncoupling of both processes. The authors combine their 3'-end mapping data and T4P-mNET-seq data mapping terminating RNAPII to score cleavage-termination distance at individual genes and find shorter distances to correlate with increased gene expression in the different cell lines. Overall, this is a carefully conducted study, and the claims and conclusions are well supported by the data.
I have some minor comments: 1. PLA assay to quantify transcription-replication conflicts (Figure 2). The quantified data looks very convincing and is also in good agreement with the proximal shift in transcription termination that is demonstrated later in the paper. However, the PLA channel signal in the microscopy image examples shown in panel A looks very blurry, and it is hard to imagine that one would be able to count # foci based on this. This may just be an issue with the resolution of the image provided. Apparently, there are much less foci in the treated samples shown in panel B - maybe microscopy images for these could be provided as well? Also, since none of the treatments impact the # of TRCs, it would have been nice to include a positive control known to induce TRCs to demonstrate that the assay works (if such a control is known) - this is optional, and I would not ask to repeat the entire experiment just for this additional control (but maybe the authors have done it and the data is already available?).
We apologize for the low resolution of the picture presented in Figure 2. We were unable to upload high resolution picture file during the first submission, for technical reasons. We will improve it in the revised manuscript.
The difference in baseline PLA foci between Fig. 2A and 2B reflects a known sensitivity of the PLA assay to cell confluency. As these two experiments were performed at different confluences, direct cross-panel comparison is not appropriate. For this reason, all quantitative comparisons in the manuscript are made strictly within the same plate, the same PLA reaction, and between wells with comparable confluency, which avoids introducing bias from these technical variables. For clarity, we plan to incorporate the above information into the Methods section. To validate assay specificity within each experiment, we confirmed that EdU-positive cells consistently showed higher PLA foci counts than EdU-negative cells from the same wells, demonstrating that quantification reflects genuine PCNA-associated signal above background. With this internal validation in place, each panel's comparisons remain valid and interpretable on their own terms.
No classical positive control exists for a PolII-pThr4/PCNA PLA interaction, as this is a relatively unexplored proximity event with no established positive control condition. We used single-antibody negative controls to establish assay specificity, although we didn't quantify and show it. We also used EdU-negative cells within the same wells as an internal background baseline, ensuring that measured foci reflect genuine signal above background. As a proxy for positive controls, we relied on the detection of changes in PLA foci number between the tested conditions, such as the effect of 4h XRN2 degradation. Also, the consistency of biological replicates and the differences between cell lines made us quite secure we were detecting reproducible and biologically relevant differences.
- Figure 2A-C: please include information on number of cells quantified
We will incorporate this information into the revised manuscript.
- Figure 2C: In the label, please include degron, e.g. HCT116 CPSF73-AID rather than just HCT116
We will modify the label according to the reviewer's suggestion.
- Figure 5C: When quantifying nascent txn based on mNET-seq, to which extent would one expect terminally paused RNAPII along the gene body (premature termination events) to contribute to the increased signal? That is, could an increase in stalling be mistaken for an increase in transcription? Based on the metagene plot in Fig 2A it doesn't look like it, but the authors may be able to estimate the effect (if any) from their data.
We thank the reviewer for pointing this out.
As reviewer #1 observed, and we comment above (Rev.1 point 2b), the increase of premature termination events in cancer cells, which can be readily detected by RNAPII T4ph mNET-seq increase in the gene body, does not globally perturb total RNAPII mNET-seq profiles (see metagenes in figure 4A and 4B).
Nevertheless, mNET-seq method does indeed detect both nascent transcription levels and RNAPII pausing, which is particularly relevant when wanting to make conclusions on a single gene level. In order to increase the robustness of our study and make stronger conclusions about nascent transcription rates, independent of stalling, during revision we will perform POINT-seq experiments in SW480 and SW620 cells. That method, in contrast to mNET-seq, is not pausing sensitive.
Reviewer #2 (Significance (Required)):
The observed uncoupling of poly-A site selection and size of termination window is unexpected and raises important questions on how these coupled processes can be regulated independently.
Strengths of the study: i) Parallel assessment of different CRC-based cell lines provides evidence of phenotype stability across patients. ii) Brings together strong technical expertise combining different state-of-the-art methodologies to map and correlate poly-A site usage, site of transcription termination, and levels of nascent transcription within the same cell lines under the same conditions, providing a comprehensive dataset.
Limitations: i) For the time being, observation limited to CRC cell lines.
While this is the first time that we are able to show the pre-mRNA 3' cleavage and transcription termination uncoupling so clearly, we have previously reported findings in other cell types which pointed to this direction. We found in HeLa cells (PMID: 30819644) that genes preferentially using distal polyadenylation sites exhibit more proximal RNAPII terminal pausing compared to genes that predominantly use proximal polyadenylation sites. Recently, we also found in U2OS cells after SETD2 KO and renal cell carcinoma cell lines with SETD2 mutation, that readthrough transcription occurs independently of APA (doi: 10.1101/2025.07.14.664659). This phenomenon could be frequent, but it has not been investigated until now, as cleavage and termination were usually studied separately.
In terms of the correlation between cleavage-termination distance and expression levels, in our study so far, we found it in CRC (HCT116, SW480, SW620) and cervical cancer (HeLa) cell lines. During revision we plan to test it additionally in pancreatic cell lines, with high sensitivity to JTE-607 treatment (BxPC3), medium (Panc1), and low sensitivity (MiaPaCa2).
ii) Mechanism behind proximal shift of termination to be determined.
We agree with the reviewer that the mechanism underlying the proximal transcription termination is missing. Our unpublished data show correlation between RNAPII pausing and transcription termination factors occupancy on chromatin. However, since more factors are involved, such as elongation speed and chromatin architecture, resolving the mechanism requires further extensive studies.
I expect this work to be of interest to an audience interested in transcription and regulation of gene expression more broadly, with potential translational relevance for cancer therapy.
Reviewer #3 (Evidence, reproducibility and clarity (Required)):
The manuscript by Stępień et al. aims to investigate the roles of pre-mRNA 3′ end processing and transcription termination using colorectal cancer (CRC) cell lines (normal colon epithelial cells: 1CT; primary tumours: HCT116 and SW480; metastatic tumour: SW620). By using publicly available proteomic datasets and their cellular models, the authors first demonstrate elevated expression of several cleavage and polyadenylation (CPA) and termination factors, including CPSF73 and PCF11, in CRC cells. They further assess the functional relevance of CPSF73 and PCF11, showing that siRNA-mediated knockdown of these factors reduces colony formation, particularly in primary cancer cell lines. However, they do not observe a clear association between CPA/termination factors and transcription-replication collisions (TRCs), suggesting that TRCs may not underlie the altered colony formation phenotype.
The authors also examine alternative polyadenylation (APA) using the CPSF73 inhibitor JTE-607 and report complex APA patterns: primary tumour cell lines display a bias toward distal APA site usage, whereas the metastatic SW620 line preferentially uses proximal sites. They further evaluate transcription termination and observe a proximal shift (early termination) in primary CRC cells, and to a lesser extent in SW620 cells. Noting the apparent discrepancy between APA site shifting and proximal termination, the authors introduce a new metric termed cleavage-termination distance, defined as the distance between the coordinates of the major PAS and RNAPII termination site. They report an association between shortening of cleavage-termination distance and increased gene expression, which may contribute to the upregulation of cancer-related genes.
Overall, this is a well-written manuscript that highlights potential roles of pre-mRNA processing and transcription termination in gene expression control, with implications for cancer biology. Nevertheless, several issues should be addressed to strengthen the study:
Overall, this is a well-written manuscript that reveals potential roles of pre-mRNA processing and transcription termination in gene expression control, with implications for cancer biology. Nevertheless, I have a few comments that may help strengthen the study.
Major comments: 1. The study includes two primary tumour cell lines but only one metastatic cell line (SW620), which is derived from the same patient as SW480. It remains unclear whether the observed effects represent general characteristics of metastatic tumour cells or are specific to this particular cell line.
Our primary workhorse in this study are the cell lines SW480 and SW620, which are derived from the same patient, to avoid the confounding variable of genetic diversity between cell lines. Unfortunately, these are the only paired CRC cell lines currently available in cell banks.
We would not want to perform further (expensive and time-consuming) genomic assays on additional CRC metastatic cell lines since the cell lines available were isolated from other types of metastasis (liver or lung, while SW620 comes from lymph node) and other patients - which would make interpreting any results obtained with them difficult. However, we plan to check the sensitivity of one more primary (HCT116) and two more metastatic (T87 and COLO-205) cancer cell lines to JTE-607 treatment in colony formation or MTT assay to find out whether the differences in CRC cell sensitivity are more cancer-stage or patient specific.
Further on, we plan to check whether our finding of alterations in cleavage-termination distance might have clinically relevant prognostic value, even outside of the context of CRC. To this end, we will test the hypothesis that a short cleavage-termination distance could be a prognostic marker for sensitivity of cells to JTE-607 treatment. It has been previously demonstrated that pancreatic cancer (PC) cell lines differ in sensitivity to JTE-607 (PMID: 38191171). We will perform T4ph-mNET-seq and nuclear 3'mRNA-seq experiments on PC cell lines to check the cleavage-termination distance in JTE-607-sensitive (BxPC3), medium sensitive (Panc1) and least JTE-607 sensitive (MiaPaCa2) cells, and for presence or absence of correlation of this distance with the cell sensitivity to JTE-607.
The rationale for focusing on colorectal cancer in this study requires further clarification. Although the Introduction provides a comprehensive review of CPSF73 and PCF11 in other cancer types, evidence specific to colorectal cancer is limited. Are these factors known to be mutated or dysregulated in CRC? Is their expression associated with patient survival? The authors could strengthen their rationale by performing a basic analysis using publicly available datasets (e.g., TCGA), such as evaluating expression levels in tumour versus normal tissue and generating Kaplan-Meier survival curves.
We will respond to these questions in the revision.
- In Figure 5 and Supplementary Figure 5, the authors analyse cleavage-termination distance across oncogenes and tumour suppressor genes and observe a negative correlation between cleavage-termination distance and gene expression level. This is an interesting finding and suggests a possible mechanism for enhancing expression of cancer-related genes. It would be valuable to extend this analysis more systematically-for example, by stratifying genes based on cleavage-termination distance and performing gene ontology enrichment analysis / GSEA to identify functional categories enriched among genes with shorter or longer distances. The authors could further relate these gene sets to, for example, distinct phenotypes between primary vs metastatic tumours.
This is an excellent suggestion. We will perform the above analyses carefully during the revision. Our initial analysis done upon receiving the reviews suggests that the genes, whose cleavage-termination distance decreases during tumorigenesis, while gene expression increases, are enriched for RNA processing, DNA damage response, chromatin organization and ribosome biogenesis factors. On the other hand, increased cleavage-termination distance and decreased gene expression are mostly associated with organelle assembly and protein localization. We will deepen this analysis and discuss the implication to cancer biology in our revised manuscript.
Minor comments: 4. In Figure 2A, the number of RNAPII-PCNA PLA foci appear comparable between SW480 and SW620, whereas in Figure 2B this seems to be much lower in SW620 compared to SW480. Could the authors clarify this discrepancy?
The difference in baseline PLA foci between Fig. 2A and 2B reflects a known sensitivity of the PLA assay to cell confluency. As these two experiments were performed at different confluencies, direct cross-panel comparison is not appropriate. For this reason, all quantitative comparisons in the manuscript are made strictly within the same plate, the same PLA reaction, and between wells with comparable confluency, which avoids introducing bias from these technical variables. For clarity, we plan to incorporate the above information into the Methods section. To validate assay specificity within each experiment, we confirmed that EdU-positive cells consistently showed higher PLA foci counts than EdU-negative cells from the same wells, demonstrating that quantification reflects genuine PCNA-associated signal above background. With this internal validation in place, each panel's comparisons remain valid and interpretable on their own terms.
- Is the cleavage-termination distance metric influenced by gene length? If so, should this parameter be normalised to gene length to avoid potential bias?
No, gene length is not a bias in the cleavage-termination distance.
- We performed correlation analysis and there is no significant correlation between the cleavage-termination distance and gene length, in any of cell line pairs in our model: HCT116 vs 1CT (spearman r=0.001, p=0.945); SW480 vs 1CT (spearman r=0.036, p=0.0654); SW620 vs 1CT (spearman r=-0.018, p=0.325).
- Additionally, we quantified the decrease in cleavage-termination distance on the very same gene, just in different cell lines. We will incorporate this result into the manuscript.
- The data and analysis scripts generated in this study have not yet been made publicly available and therefore cannot be fully evaluated.
We apologize for this omission. The revised manuscript will contain the link to our publicly available scripts in GitHub and the GEO access.
**Referees cross-commenting**
I agree with the reports from both Reviewer #1 and Reviewer #2.
I would like to thank Reviewer #1 for pointing out my mistaken. The authors did not use JTE-607 to study APA; rather, they studied the differences in APA between cell lines. I apologise for the confusion.
Reviewer #3 (Significance (Required)):
General assessment: This study investigates the contribution of pre-mRNA 3′ end processing and transcription termination to colorectal cancer (CRC) biology using a combination of cell line comparisons (primary versus metastatic tumours), chemical, and RNAi perturbations, and bioinformatic analyses.
The major strengths of the work include: • The use of CRC cell lines representing normal, primary, and metastatic states, including matched primary and metastatic lines derived from the same patient. • A systematic analysis of alternative polyadenylation (APA) and transcription termination, revealing a potential uncoupling between these two closely related processes. • The introduction of a novel quantitative metric-cleavage-termination distance-to examine the relationship between PAS usage and RNAPII termination. • The identification of a negative association between cleavage-termination distance and gene expression, suggesting an additional regulatory layer influencing gene expression.
However, certain limitations should be considered: • The generalisability of conclusions regarding metastatic CRC is limited by reliance on a single metastatic cell line.
We believe that the experiments we outlined above in response to Reviewer #3 point 1 will allow us to extend the generalizability of conclusion.
• The translational relevance of the findings could be further strengthened through patient-level or clinical data analysis.
We agree with the reviewer. Due to technical limitations, it is not possible to perform nascent transcriptomic experiments on patient material at this time. However, we will attempt to strengthen the translational relevance by additional experiments and analysis as indicated in response to Reviewer #3 points 1-3.
Advance: The study proposes potentially novel roles for 3′ end cleavage and transcription termination in regulating gene expression in colorectal cancer. In particular, the conceptual distinction between APA site shifting and transcription termination, together with the introduction of the cleavage-termination distance metric, represents a conceptual advance.
Audience: The work is primarily positioned within basic research. With additional translational context, it may also attract interest from a broader audience.
Field of expertise: transcriptional regulation and bioinformatics

