Cohesin mutations are synthetic lethal with stimulation of WNT signaling
Abstract
Mutations in genes encoding subunits of the cohesin complex are common in several cancers, but may also expose druggable vulnerabilities. We generated isogenic MCF10A cell lines with deletion mutations of genes encoding cohesin subunits SMC3, RAD21, and STAG2 and screened for synthetic lethality with 3009 FDA-approved compounds. The screen identified several compounds that interfere with transcription, DNA damage repair and the cell cycle. Unexpectedly, one of the top ‘hits’ was a GSK3 inhibitor, an agonist of Wnt signaling. We show that sensitivity to GSK3 inhibition is likely due to stabilization of β-catenin in cohesin-mutant cells, and that Wnt-responsive gene expression is highly sensitized in STAG2-mutant CMK leukemia cells. Moreover, Wnt activity is enhanced in zebrafish mutant for cohesin subunits stag2b and rad21. Our results suggest that cohesin mutations could progress oncogenesis by enhancing Wnt signaling, and that targeting the Wnt pathway may represent a novel therapeutic strategy for cohesin-mutant cancers.
Introduction
The cohesin complex is essential for sister chromatid cohesion, DNA replication, DNA repair, and genome organization. Three subunits, SMC1A, SMC3, and RAD21, form the core ring-shaped structure of human cohesin (Dorsett and Ström, 2012; Horsfield et al., 2012). A fourth subunit of either STAG1 or STAG2 binds to cohesin by contacting RAD21 and SMC subunits (Shi et al., 2020), and is required for the association of cohesin with DNA (Dorsett and Ström, 2012; Horsfield et al., 2012; Shi et al., 2020). The STAG subunits of cohesin are also capable of binding RNA in the nucleus (Pan et al., 2020). Cohesin associates with DNA by interaction with loading factors NIPBL and MAU2 (Wendt, 2017), its stability on DNA is regulated by the acetyltransferases ESCO1 (Wutz et al., 2020) and ESCO2 (van der Lelij et al., 2009), and its removal is facilitated by PDS5 and WAPL (Wutz et al., 2017; Shintomi and Hirano, 2009). The cohesin ring itself acts as a molecular motor to extrude DNA loops, and this activity is thought to underlie its ability to organize the genome (Vian et al., 2018; Mayerova et al., 2020). Cohesin works together with CCCTC-binding factor (CTCF) to mediate three-dimensional genome structures, including enhancer-promoter loops that instruct gene accessibility and expression (Hansen, 2020; Rowley and Corces, 2018). The consequences of cohesin mutation therefore manifest as chromosome segregation errors, DNA damage, and deficiencies in genome organization leading to gene expression changes.
All cohesin subunits are essential to life: homozygous mutations in genes encoding complex members are embryonic lethal (Horsfield et al., 2012). However, haploinsufficient germline mutations in NIPBL, ESCO2, and in cohesin genes, cause human developmental syndromes known as the ‘cohesinopathies’ (Horsfield et al., 2012). Somatic mutations in cohesin genes are prevalent in several different types of cancer, including bladder cancer (15–40%), endometrial cancer (19%), glioblastoma (7%), Ewing’s sarcoma (16–22%) and myeloid leukemias (5–53%) (De Koninck and Losada, 2016; Hill et al., 2016; Waldman, 2020). The prevalence of cohesin gene mutations in myeloid malignancies (Kon et al., 2013; Papaemmanuil et al., 2016; Thol et al., 2014; Thota et al., 2014; Yoshida et al., 2013) reflects cohesin’s role in determining lineage identity and differentiation of hematopoietic cells (Galeev et al., 2016; Mazumdar et al., 2015; Mullenders et al., 2015; Viny et al., 2015). Of the cohesin genes, STAG2 is the most frequently mutated, with about half of cohesin mutations in cancer involving STAG2 (Waldman, 2020).
While cancer-associated mutations in genes encoding RAD21, SMC3, and STAG1 are always heterozygous (Thota et al., 2014; Kon et al., 2013; Tsai et al., 2017), mutations in the X chromosome-located genes SMC1A and STAG2 can result in complete loss of function due to hemizygosity (males), or silencing of the wild type during X-inactivation (females). STAG2 and STAG1 have redundant roles in cell division, therefore complete loss of STAG2 is tolerated due to partial compensation by STAG1. Loss of both STAG2 and STAG1 leads to lethality (Benedetti et al., 2017; van der Lelij et al., 2017). STAG1 inhibition in cancer cells with STAG2 mutation causes chromosome segregation defects and subsequent lethality (Liu et al., 2018). Therefore, although partial depletion of cohesin can confer a selective advantage to cancer cells, a complete block of cohesin function will cause cell death. The multiple roles of cohesin provide an opportunity to inhibit the growth of cohesin-mutant cancer cells via chemical interference with pathways that depend on normal cohesin function. For example, poly ADP-ribose polymerase (PARP) inhibitors were previously shown to exhibit synthetic lethality with cohesin mutations (Waldman, 2020; Liu et al., 2018; Mondal et al., 2019; McLellan et al., 2012; O'Neil et al., 2013). PARP inhibitors prevent DNA double-strand break repair (Zaremba and Curtin, 2007), a process that also relies on cohesin function.
To date, only a limited number of compounds have been identified as inhibitors of cohesin-mutant cells (Waldman, 2020). Here, we sought to identify additional compounds of interest by screening libraries of FDA-approved molecules against isogenic MCF10A cells with deficiencies in RAD21, SMC3, or STAG2. Unexpectedly, our screen identified a novel sensitivity of cohesin-deficient cells to a GSK3 inhibitor that acts as an agonist of the Wnt signaling pathway. We found that β-catenin stabilization upon cohesin deficiency likely contributes to an acute sensitivity of Wnt target genes. The results raise the possibility that sensitization to Wnt signaling in cohesin-mutant cells may participate in oncogenesis, and suggest that Wnt agonism could be therapeutically useful for cohesin-mutant cancers.
Results
Cohesin gene deletion in MCF10A cells results in minor cell cycle defects
To avoid any complications with pre-existing oncogenic mutations, we chose the relatively ‘normal’ MCF10A line for creation and screening of isogenic deletion clones of cohesin genes SMC3, RAD21, and STAG2. MCF10A is a near-diploid, immortalized, breast epithelial cell line that exhibits normal epithelial characteristics (Tait et al., 1990) and has been successfully used for siRNA and small molecule screening (Telford et al., 2015). Two sgRNAs per gene were designed targeting the 5' and 3' UTR regions, respectively, of RAD21, SMC3, and STAG2 genes. Single cells were isolated and grown into clones that were genotyped for complete gene deletions (Figure 1, Supplementary file 1). We isolated several RAD21 and SMC3 deletion clones, and selected single clones for further characterization that grew normally and were essentially heterozygous. In the selected RAD21 deletion clone, one of three RAD21 alleles (on chromosome 8, triploid in MCF10A) was confirmed deleted, with one wild type and one undetermined allele. In the selected SMC3 deletion clone, one of two SMC3 alleles (on chromosome 10) was deleted. In the selected STAG2 deletion clone, homozygous loss of STAG2 was determined by the absence of STAG2 mRNA and protein. For convenience here, we named the three cohesin-mutant clones RAD21+/-, SMC3+/-, and STAG2-/-.

Creation of MCF10A isogenic cell lines with cohesin gene deletions.
(A) Top, schematic diagram shows the deletion strategy for genes encoding cohesin subunits RAD21, SMC3, and STAG2 using two sgRNAs targeting the 5'UTR and the 3'UTR of each gene. Bottom, heterozygous clones were identified by PCR using specific primer pairs flanking the deletion region. Representative DNA gel shows the PCR products yielded using specific primer pairs for MCF10A parental and RAD21+/- deletion clone. M, ladder marker. (B,C,D) Schematic deletion strategy and summary of the allele sequences for the STAG2 homozygous deletion clone, and the RAD21 and SMC3 heterozygous deletion clones. (E) RNA levels of the targeted genes in MCF10A cohesin-deficient clones. (F) Representative immunoblot and (G) quantification of cohesin protein levels. γ-tubulin was used as loading control. n = 3 independent experiments, mean ±s.d., one-tailed student t test: **p≤0.01; ****p≤0.0001. Guide RNAs and PCR primers can be found in Supplementary file 1.
The RAD21+/-, SMC3+/-, and STAG2-/- clones had essentially normal cell cycle progression when compared to parental cells, although the RAD21+/- and SMC3+/- clones proliferated more slowly than the others (Figure 2A,B). Chromosome spreads revealed that only the STAG2-/- clone had noteworthy chromosome cohesion defects characterized by partial or complete loss of chromosome cohesion and gain or loss of more than one chromosome (Figure 2C,D; Figure 2—figure supplement 1A). Remarkably, SMC3+/- cells had noticeably larger nuclei that appeared to be less dense (Figure 2E,F; Figure 2—figure supplement 1B), possibly owing to decompaction of DNA in this clone. RAD21+/- and SMC3+/- clones exhibited occasional lagging chromosomes and micronuclei, while the STAG2-/- clone did not (Figure 2G; Figure 2—figure supplement 1C–E). Cell growth and morphology in monolayer culture was essentially normal in all three cohesin-mutant clones (Figure 2—figure supplement 2).

minor deficiencies in cell cycle progression, chromosome segregation and nuclear morphology in cohesin-deficient isogenic cell lines.
(A) Proliferation curves of MCF10A parental and cohesin-deficient clones. n = 3 independent experiments, mean ± sd, two-way ANOVA: **p≤0.01; ****p≤0.0001. Doubling times in hours of MCF10A parental cells are 15.5 ± 0.1; RAD21+/-, 17.1 ± 0.2; SMC3+/-, 21.4 ± 0.9; STAG2-/-, 16.4 ± 0.2 respectively. (B) Flow cytometry analysis of cell cycle progression. (C) Representative metaphase spread images of cohesin-deficient cells. Black arrow indicates a normal chromosome; Red arrow indicates premature chromatid separation (PCS). (D) Quantification of chromosome cohesion defects. A minimum of 20 metaphase spreads were examined per individual experiment, n = 2 independent experiments, mean ±s.d. (E) Representative confocal images of nuclear morphology. Scale bar, 15 μM. Yellow arrow indicates a micronucleus. (F) Quantification of nuclear area. (G) Quantification of micronuclei (MN). A minimum of 1000 cells was examined per individual experiment. n = 5 independent experiments (square symbols), mean ±s.d., one-way ANOVA: *p≤0.05; ****p≤0.0001. Source data is available for A,D,F,G in Figure 2—source data 1.
-
Figure 2—source data 1
Raw data for Figure 2a,d,f,g.
- https://cdn.elifesciences.org/articles/61405/elife-61405-fig2-data1-v2.xlsx
Overall, the cohesin deletion clones appear to have infrequent but specific cell cycle anomalies that are shared between some, but not all clones. Anomalies include loss of chromosome cohesion, lagging chromosomes, or micronuclei, but these features do not appear to majorly impact on cell cycle progression or morphology.
Cohesin-depleted cells have altered nucleolar morphology and are sensitive to DNA damaging agents
Cohesin-deficient cells have been demonstrated to display compromised nucleolar morphology and ribosome biogenesis (Bose et al., 2012; Harris et al., 2014), as well as sensitivity to DNA damaging agents (Bailey et al., 2014; Mondal et al., 2019). Analysis of our RAD21+/-, SMC3+/-, and STAG2-/- clones confirmed these findings. Cohesin-deleted cells in steady-state growth had abnormal nucleolar morphology, as revealed by fibrillarin and nucleolin staining (Figure 3—figure supplement 1). Furthermore, we found that treatment with DNA intercalator/transcription inhibitor Actinomycin D caused marked fragmentation of nucleoli in all three cohesin-deficient cell lines as determined by fibrillarin staining (Figure 3A,B). Actinomycin D treatment increased γ-H2AX and TP53 in the nuclei of RAD21+/- and SMC3+/- cells, in particular, indicating that these cells are compromised for DNA damage repair relative to the parental MCF10A cells (Figure 3C–F). In contrast, immunostaining for γ-H2AX and TP53 levels were comparable at baseline between cohesin-deficient clones and parental cells. Interestingly, the STAG2-/- deletion clone was much more resistant to DNA damage caused by Actinomycin D (Figure 3C–F), even though nucleoli are abnormal in this clone (Figure 3—figure supplement 1).

Cohesin-deficient cells have increased sensitivity to nucleolar stress and DNA damaging agents.
(A) Representative image and (B) quantification of nucleolar dispersal observed in parental (WT) MCF10A cells and cohesin-deficient clones after exposure to a DNA damaging agent, Actinomycin D (ActD) 8 nM. Fibrillarin (FBL) staining was used as a marker for nucleoli. (C) Representative image and (D) quantification of DNA damage foci observed in parental (WT) MCF10A cells and cohesin-deficient clones after exposure to ActD. An antibody detecting γH2AX was used to visualize foci of DNA double-strand breaks. (E) Representative image and (F) quantification of nuclear p53 in parental (WT) MCF10A cells and cohesin-deficient clones after exposure to ActD. A minimum of 500 cells was examined per individual experiment. n = 3 independent experiments, mean ± s.d., one-way ANOVA: *p≤0.05; **p≤0.01; ***p≤0.0005. Scale bar, 15 μM. Source data is available for Figure 2B,D,F in Figure 3—source data 1.
-
Figure 3—source data 1
Raw data for Figure 3.
- https://cdn.elifesciences.org/articles/61405/elife-61405-fig3-data1-v2.xlsx
Overall, increased susceptibility of RAD21+/- and SMC3+/- clones to DNA damage is consistent with cohesin’s role in DNA double-strand break repair (Sjögren and Ström, 2010), and highlights the different requirements for RAD21 and SMC3 versus STAG2 in this process. Cohesin mutations were previously shown to sensitize cells to PARP inhibitor, Olaparib (Mondal et al., 2019; Matto et al., 2015). We confirmed a mild to moderate sensitivity to Olaparib in our cohesin-deficient MCF10A clones relative to parental cells (Figure 3—figure supplement 2).
Collectively, characterization of our cohesin-deficient clones provided confidence that they represent suitable models for synthetic lethal screening. To confirm that the phenotypes of our chosen clones are representative, we selected a further two clones with heterozygous deletions in SMC3 and RAD21, and monitored their growth, morphology and chromosome cohesion (Figure 3—figure supplement 3). These analyses showed that the two additional deletion clones were similar to those already characterized (Figures 1–3), providing surety that our cohesin-deficient cell lines have properly representative phenotypes. Furthermore, none of the cohesin-deficient clones exhibited enhanced cell death compared with parental MCF10A cells (Figure 2B; Figure 3—figure supplement 3).
A synthetic lethal screen of FDA-approved compounds identifies common sensitivity of cohesin-mutant cells to WNT activation and BET inhibition
To identify additional compounds that inhibit the growth of cohesin-deficient cells, we screened the cohesin-deficient MCF10A cell lines with five dose concentrations (1–5000 nM) of 3009 compounds, including FDA-approved compounds (2399/3009), kinase inhibitors (429/3009), and epigenetic modifiers (181/3009) (Figure 4A). DMSO and Camptothecin were included as negative and positive viability controls, respectively (Figure 4—figure supplement 1). We assayed cell viability after 48 hr of compound treatment.

A synthetic lethal screen identifies common sensitivity of cohesin-deficient cells to WNT activation and BET inhibition.
(A) Schematic overview of the synthetic lethal screen. (B,C,D) Overview of the differential area over the curve (AOC) activity of all compounds tested in cohesin-deficient cell lines relative to parental MCF10A cells in the primary screen. A threshold of differential AOC ≥ 0.15 (red dashed lines) was used to filter candidate compounds of interest. (E) Venn diagram showing the number of common and unique compounds that inhibited RAD21+/-, SMC3+/-, and STAG2-/- in the primary screen. (F,G) Dose-response curves of I-BET-762 and LY2090314. Source data is available for B–E in Figure 4—source data 1.
-
Figure 4—source data 1
Raw cell counts for SL screen compound treatments.
- https://cdn.elifesciences.org/articles/61405/elife-61405-fig4-data1-v2.xlsx
-
Figure 4—source data 2
Table and AOC measurements of all hit compounds from the screen.
- https://cdn.elifesciences.org/articles/61405/elife-61405-fig4-data2-v2.xlsx
Synthetic lethal candidate compounds were ranked based on the differential area over curve (AOC) values that are derived from a growth rate-based metric (GR) (Figure 4B–D; Figure 4—source data 1). The GR metric takes into account the varying growth rate of dividing cells to mitigate inconsistent comparison of compound effects across cohesin-deficient cell lines (Hafner et al., 2016). Compounds that caused ≥30% growth inhibition in cohesin-deficient clones compared with parental MCF10A cells were selected for further analysis. The screen identified candidate 206 synthetic lethal compounds, of which 18 inhibited all three cohesin-deficient cell lines ≥ 30% more than the parental MCF10A cells (Figure 4E; Figure 4—source data 2; Table 1; Table 1—source data 1).
Significant inhibitors of all three cohesin-deficient clones.
Compound | Rank 1° screen | Rank 2° screen | Target | Pathway | 1° differential AOC activity | ||
---|---|---|---|---|---|---|---|
RAD21+/- | SMC3+/- | STAG2-/- | |||||
WAY-600 | 1 | 3 | mTORC1/2 | PI3K/AKT/mTOR | 0.30 | 0.37 | 0.32 |
I-BET-762 | 2 | 8 | BET proteins | Epigenetics | 0.25 | 0.27 | 0.29 |
LY2090314 | 3 | 6 | GSK3 | WNT | 0.25 | 0.22 | 0.25 |
Vistusertib (AZD2014) | 4 | 1 | mTORC1/2 | PI3K/AKT/mTOR | 0.18 | 0.33 | 0.26 |
P276-00 | 5 | 17 | CDK1/4/9 | Cell Cycle | 0.17 | 0.29 | 0.28 |
MK-8745 | 6 | 22 | Aurora A | Cell Cycle | 0.18 | 0.30 | 0.26 |
Ethacridine lactate | 7 | 28 | Anti-infection | Microbiology | 0.18 | 0.25 | 0.27 |
CUDC-101 | 8 | 36 | EGFR, HDAC, HER2 | Epigenetics | 0.16 | 0.25 | 0.24 |
Dabrafenib (GSK2118436) | 9 | 9 | BRAF | MAPK | 0.18 | 0.27 | 0.18 |
SAR131675 | 10 | 34 | VEGFR3 | Protein Tyrosine Kinase | 0.15 | 0.17 | 0.33 |
ZM 447439 | 11 | 41 | Aurora A/B | Cell Cycle | 0.16 | 0.22 | 0.24 |
Gitoxigenin Diacetate | 12 | 32 | NA | Other | 0.17 | 0.24 | 0.20 |
UNC669 | 13 | 62 | Epigenetic Reader Domain | Epigenetics | 0.18 | 0.23 | 0.20 |
4-Phenylbutyric Acid | 14 | 29 | Endoplasmic reticulum stress | Other | 0.20 | 0.17 | 0.20 |
Ipatasertib (GDC-0068) | 15 | 12 | AKT | PI3K/AKT/mTOR | 0.21 | 0.19 | 0.16 |
VX-702 | 16 | 43 | P38 MAPK | MAPK | 0.15 | 0.19 | 0.21 |
RVX-208 | 17 | 58 | BET proteins | Epigenetics | 0.16 | 0.20 | 0.17 |
Dihydroergotamine mesylate | 18 | 49 | NA | Other | 0.16 | 0.19 | 0.16 |
Olaparib | 351 | PARP1/2 | DNA Damage | 0.11 | 0.13 | 0.22 |
-
Table 1—source data 1
Compounds with growth inhibitory activity (AOC) ranked to produce Table 1.
- https://cdn.elifesciences.org/articles/61405/elife-61405-table1-data1-v2.xlsx
-
Table 1—source data 2
Compounds effective in the secondary screen ranked to produce Table 1.
- https://cdn.elifesciences.org/articles/61405/elife-61405-table1-data2-v2.xlsx
Most of the 206 compounds inhibited at least two cohesin-deficient cell lines and were classed in similar categories of inhibitor (Figure 4—figure supplement 2A–D; Figure 4—source data 2; Table 1). Of the 206 primary screen hits, 85 (including the 18 that inhibited all three cohesin-deficient clones, plus the most effective candidates from each inhibitor category) were subjected to secondary screening in an 11-point dose curve ranging from 0.5 nM to 10 μM. Notable sensitive pathways included: DNA damage repair, the PI3K/AKT/mTOR pathway, epigenetic control of transcription, and stimulation of the Wnt signaling pathway (Table 1; Table 1—source data 2; Figure 4—source data 2).
The identification of DNA damage repair inhibitors in our screen agrees with previous studies showing synthetic lethality of PARP inhibition with cohesin mutation (McLellan et al., 2012). Differential sensitivity of the cohesin deletion cell lines to PI3K/AKT/mTOR inhibitors is consistent with the observed nucleolar deficiencies in cohesin-mutant cells (Figure 3—figure supplement 1). The PI3K/AKT/mTOR pathway stimulates ribosome biogenesis and rDNA transcription; because rDNA is contained in nucleoli, it is likely that rDNA transcription and ribosome production is already compromised in the cohesin gene deletion cell lines. We had previously shown that BET inhibition is effective in blocking precocious gene expression in the chronic myelogenous leukemia cell line K562 containing a STAG2 mutation (Antony et al., 2020). Growth inhibition of cohesin-deficient MCF10A cells by I-BET-762 (Figure 4F; Figure 4—figure supplement 2E) reinforces the idea that targeting BET could be therapeutically effective in cohesin-mutant cancers.
Interestingly, we found that LY2090314, a GSK3 inhibitor and stimulator of the Wnt signaling pathway, inhibits all three cohesin deletion lines (Table 1, Figure 4G, Figure 4—figure supplement 2F). LY2090314 also inhibited the growth of K562 STAG2 R614* mutant leukemia cells that we had previously characterized (Antony et al., 2020; Figure 4—figure supplement 2G) as well as the two additional MCF10A RAD21- and SMC3-deficient clones (Figure 3—figure supplement 3G). Wnt signaling appears to act upstream of cohesin (Xu et al., 2014; Estarás et al., 2015), and also to be primarily downregulated downstream of cohesin mutation (Mills et al., 2018; Schuster et al., 2015; Avagliano et al., 2017). Therefore, we were prompted to further investigate why Wnt stimulation via GSK3 inhibition might cause lethality in cohesin-deficient cells.
Stabilization of β-catenin in cohesin-deficient MCF10A cells
In ‘off’ state of canonical Wnt signaling, β-catenin forms a complex including Axin, APC, GSK3, and CK1 proteins. β-catenin phosphorylated by CK1 and GSK3 is recognized by the E3 ubiquitin ligase β-Trcp, which targets β-catenin for proteasomal degradation. Activation of Wnt signaling by ligand binding, or by GSK3 inhibition, releases β-catenin, allowing it to accumulate in the nucleus where it binds TCF to activate Wnt target gene transcription (Logan and Nusse, 2004; Clevers et al., 2014). We found that there was no difference in GSK3 levels between parental MCF10A cells and the cohesin-deficient clones, and upon treatment of cells with LY2090314, levels of GSK3 were reduced in all cells as expected (Figure 5—figure supplement 1). In contrast, in untreated cells we found that β-catenin is stabilized in all three cohesin-deficient clones compared with parental MCF10A cells. When Wnt signaling is in the ‘off’ state, phospho-Ser33/37/Thr41 marks β-catenin targeted for degradation. Membrane-associated phospho-Ser33/37/Thr41-β-catenin was increased in cohesin-deficient clones, and disappeared upon treatment with LY2090314 (Figure 5A). Phospho-Ser675 β-catenin, the active form of β-catenin that is induced upon Wnt signaling, was noticeably increased in the cytoplasm of cohesin-deficient cells following treatment with LY2090314 (Figure 5B,C). The results imply that inactive β-catenin that would normally be targeted for degradation is instead stabilized in cohesin-deficient cells, and is available to be converted into the active form following inhibition of GSK3.

LY2090314-mediated WNT stimulation leads to β-catenin stabilization in cohesin-deficient MCF10A cells.
(A) Immunoblot of the membrane fraction of parental (WT) and cohesin-deficient MCF10A cells shows increased basal level of β-catenin phosphorylation at Ser33/37/Thr41. (B) Immunoblot of the cytoplasmic fraction shows increased level of both total and phosphorylated β-catenin at Ser675 after parental (WT) and cohesin-deficient MCF10A cells were treated with LY2090314 at 100 nM for 24 hr. (C) Quantification of protein levels for total and phosphorylated β-catenin at Ser675. n = 3 independent experiments, mean ± s.d., one-way ANOVA: *p≤0.05, **p≤0.01. (D) Immunofluorescence images show cytosolic accumulation of active β-catenin in cohesin-deficient MCF10A cells treated with LY2090314 100 nM for 24 hr, relative to parental (WT) MCF10A cells. White arrows indicate puncta of β-catenin (pSer675). Scale bar = 25 μM. Full length blots and molecular size markers are available for A,B in Figure 5—source data 1. Source quantification data is available for C in Figure 5—source data 2.
-
Figure 5—source data 1
Untrimmed blots for Figure 5A, B.
- https://cdn.elifesciences.org/articles/61405/elife-61405-fig5-data1-v2.pdf
-
Figure 5—source data 2
Quantitation of blots in Figure 5A, B.
- https://cdn.elifesciences.org/articles/61405/elife-61405-fig5-data2-v2.xlsx
-
Figure 5—source data 3
Untrimmed blots for Figure 5A, B.
- https://cdn.elifesciences.org/articles/61405/elife-61405-fig5-data3-v2.pdf
-
Figure 5—source data 4
Untrimmed blots for Figure 5A, B.
- https://cdn.elifesciences.org/articles/61405/elife-61405-fig5-data4-v2.pdf
To determine if stabilization of β-catenin is conserved in a second model of cohesin-mutant cancer, we performed an identical LY2090314 treatment on HCT116 cells that were stably transfected to express two SMC1A mutations identified in human colorectal carcinomas, c.2027A > G (leading to p.E676G change near the hinge domain) and c.2479 C > T (leading to a truncated protein, p.Q827X) (Cucco et al., 2014; Sarogni et al., 2019). The limitation is that these cells are not a model of cohesin deficiency, but rather, one in which normal cohesin function is perturbed by expression of a mutant version of SMC1A (Sarogni et al., 2019). We observed an increased basal level of phosphorylated β-catenin at Ser675 in cells that express either of these SMC1A mutants. Furthermore, LY2090314 treatment markedly increased total β-catenin in both the SMC1A mutants compared with HCT116 wild type controls (Figure 5—figure supplement 2A,B). The results indicate that abnormally high levels of active β-catenin are also present following Wnt stimulation when a fourth subunit of cohesin, SMC1A, is perturbed.
Immunofluorescence labeling of MCF10A cells showed that in LY2090314-treated cells, the stabilized active phospho-Ser675 β-catenin (and total β-catenin, Figure 5—figure supplement 3B) mainly locates to puncta in the cytoplasm (Figure 5D). In contrast, no β-catenin accumulation was observed in DMSO-treated cells (Figure 5—figure supplement 3A). β-catenin-containing puncta were also observed upon WNT3A treatment of cohesin-deficient clones (Figure 5—figure supplement 3C). This indicates that Wnt stimulation rather than a secondary effect of LY2090314 is responsible for β-catenin accumulation. We could not reliably detect an increase of β-catenin in the nucleus of cohesin-deficient MCF10A clones by immunofluorescence (Figure 5D), therefore we decided to use transcriptional response to determine the consequences of β-catenin stabilization for Wnt signaling.
Cohesin-deficient leukemia cells are hypersensitive to Wnt signaling
RNA-sequencing analysis of the cohesin gene deletion clones compared with parental MCF10A cells revealed that gene expression changes strongly track with the identity of the deleted cohesin gene (Figure 6—figure supplement 1A). However, expressed transcripts encoding Wnt signaling pathway components did not cluster differently in cohesin-deficient clones compared with parental MCF10A cells, and cohesin mutation-based clustering remained dominant (Figure 6—figure supplement 1B). We reasoned that stimulation of Wnt signaling in a more responsive cell type might be necessary to determine how stabilized β-catenin in cohesin-deficient cells affects transcription.
Myeloid leukemias are frequently characterized by cohesin mutations, particularly in STAG2. Furthermore, activation of Wnt signaling is associated with transformation in AML (Wang et al., 2010; Beghini et al., 2012; Kang et al., 2020) and AML patients were identified with mutations in AXIN1 and APC that lead to stabilization of β-catenin (Erbilgin et al., 2012). CMK is a Down Syndrome-derived megakaryoblastic cell line that typifies myeloid leukemias that are particularly prone to cohesin mutation (Yoshida et al., 2013). We edited CMK to contain the STAG2-AML associated mutation R614* (Antony et al., 2020) and selected single clones with complete loss of the STAG2 protein to create the cell line CMK-STAG2-/- (Figure 6—figure supplement 2A–C).
Immunofluorescence showed that β-catenin is increased by 20% in the nuclei of CMK-STAG2-/- cells relative to parental cells upon WNT3A stimulation (Figure 6A,B). To determine immediate early transcriptional responses to Wnt signaling, RNA-sequencing was performed on CMK-STAG2-/- and parental CMK cells at baseline and after stimulation with WNT3A for 4 hr. Around 76% more genes were upregulated in response to WNT3A in CMK-STAG2-/- compared with CMK parental cells (616 in CMK-STAG2-/- vs 350 in CMK parental, FDR ≤ 0.05, Figure 6C), while around the same number were downregulated. About one quarter of differentially expressed genes overlapped between CMK-STAG2-/- and CMK parental cells.

Cohesin-STAG2 mutant CMK cells show increased sensitivity to Wnt signaling.
(A) Immunofluorescence images showing slightly increased nuclear accumulation of β-catenin in STAG2-CMK cells (STAG2-/-) compared to parental (WT) following treatment with LY2090314 at 100 nM for 24 hr. (B) Quantification of nuclear total β-catenin in parental (WT) and STAG2-CMK cells (STAG2-/-) cells. Fluorescence of nuclear total β-catenin was determined relative to the nuclear area. Image J was used quantify cells from 10 different confocal fields. The graph depicts s.e.m. from analyses of 170–188 cell nuclei, and the p value was calculated using a student’s t test. (C) Histogram showing the number of genes upregulated or downregulated at FDR ≤ 0.05 upon WNT3A treatment in parental (WT) and STAG2-/- CMK cells. WNT3A stimulation was performed on three biological replicates, from three independent experiments. (D) Overlap of genes significantly upregulated or downregulated (FDR ≤ 0.05) upon WNT3A treatment between parental (WT) and STAG2-/- CMK cells. (E) Top enriched pathways (ranked by gene count) from the significantly downregulated and upregulated genes (FDR ≤ 0.05) following WNT3A treatment in STAG2-/- and parental (WT) CMK cells using the ClusterProfiler R package on the Gene Ontology Biological Process dataset modeling for both cell types (WT or STAG2-/-) and regulation pattern (up- or downregulation). Source data is available for C,D in Figure 6—source data 1.
-
Figure 6—source data 1
Gene expression data for Figure 6.
- https://cdn.elifesciences.org/articles/61405/elife-61405-fig6-data1-v2.xlsx
Strikingly, 402 genes that were not Wnt-responsive in CMK parental cells became Wnt-responsive upon introduction of the STAG2 R614* mutation (Figure 6D). Genes upregulated in CMK-STAG2-/- markedly increased the number of Wnt-sensitive biological pathways following WNT3A treatment compared with parental CMK cells (Figure 6E). A strongly upregulated cluster of 244 transcripts in CMK-STAG2-/- included genes encoding signaling molecules (JAG2, IL6ST, SEMA3G) and transcription factors (SP7, KLF3, SMAD3, EP300, and EPHA4) (Figure 6—figure supplement 2D). Genes in this cluster were enriched for LEF1 and TCF7 binding sites, indicating their potential to be directly regulated by β-catenin. Enriched biological pathways included Wnt signaling, cell cycle and metabolism (Figure 6—figure supplement 2E). Pathway analyses of genes significantly downregulated only in CMK-STAG2-/- also showed enrichment for metabolism (Figure 6—figure supplement 2F,G).
Overall, our results show that CMK-STAG2-/- cells are exquisitely sensitive to Wnt signaling. Introduction of the STAG2 R614* mutation amplified expression of Wnt-responsive genes and sensitized genes and pathways that are not normally Wnt-responsive in CMK. This sensitivity could be due at least in part to stabilized β-catenin.
Conservation of WNT sensitivity in cohesin-deficient zebrafish
To determine if enhanced Wnt sensitivity is conserved in a cohesin-deficient animal model, we used two previously described zebrafish cohesin mutants: stag2bnz207 (Ketharnathan et al., 2020), which has a seven base-pair deletion in stag2b leading to a prematurely truncated Stag2b protein, and rad21nz171 (Horsfield et al., 2007), which has a nonsense point mutation in the rad21 gene that eliminates Rad21 protein. To provide a readout for Wnt signaling, we used transgenic zebrafish carrying a TCF/β-catenin reporter in which exogenous Wnt stimulation induces nuclear mCherry: Tg(7xTCF-Xla.Siam:nslmCherry)ia5 (Moro et al., 2012). The Wnt reporter construct, which drives nuclear mCherry red fluorescence in Wnt-responsive cells, was introduced into zebrafish carrying the stag2bnz207 and rad21nz171 mutation by crossing.
Zebrafish embryos heterozygous for Tg(7xTCF-Xla.Siam:nslmCherry)ia5, and either homozygous for stag2bnz207, rad21nz171, or wild type, were exposed to 0.15 M LiCl, an agonist of the Wnt signaling pathway (Figure 7). Expression of mCherry in the midbrain of embryos was detected at 20 hr post-fertilization (hpf) by epifluorescence and confocal imaging. A constitutive low level of mCherry is present in the developing midbrain of both untreated wild type (Figure 7A,B,I,J), and rad21nz171 embryos (Figure 7M,N). On the other hand, untreated stag2bnz207 embryos exhibited noticeably higher basal mCherry levels than wild type siblings (Figure 7E,F compared with A,B). This observation indicates that the Wnt pathway is more intrinsically active in these embryos. However, the addition of LiCl did not result in much additional fluorescence owing to stag2bnz207 mutation (Figure 7G,H compared with C,D). While mCherry expression in rad21nz171 mutants resembled that in wild type embryos at baseline, LiCl treatment dramatically increased the existing midbrain mCherry expression in rad21nz171 mutants compared with wild type embryos (Figure 7O,P compared with K,L). This observation indicates that rad21nz171 mutants are more sensitive to Wnt stimulation than wild type.

Zebrafish stag2b and rad21 cohesin mutants show increased sensitivity to Wnt signaling.
Wnt reporter Tg(7xTCF-Xla.Siam:nlsmCherry)ia5 control embryos, Tg(7xTCF-Xla.Siam:nlsmCherry)ia5;stag2bnz207 and Tg(7xTCF-Xla.Siam:nlsmCherry)ia5;rad21nz171 cohesin-mutant embryos were treated with 0.15 M LiCl from 4 hpf to 20 hpf. (A–H) max projections of 4 (10 μm) optical sections. (A,C,E,G) TD (transmitted light detector) images merged with confocal images. (B,D,F,H) confocal images alone. (I,K,M,O) Brightfield/fluorescent and (J,N,L,P) confocal images of the same embryos in I,K,M,O. (A–D) and (I–J) Tg(7xTCF-Xla.Siam:nlsmCherry)ia5 control embryos (WT) have low level fluorescence (Wnt reporter activity) in the midbrain that is increased following treatment. (E–H) Tg(7xTCF-Xla.Siam:nlsmCherry)ia5; stag2bnz207 embryos have elevated baseline levels of fluorescence (Wnt reporter activity) relative to controls (yellow arrow) with not much further increase upon LiCl treatment. (M–P) Tg(7xTCF-Xla.Siam:nlsmCherry)ia5; rad21nz171 cohesin-mutant embryos show enhanced Wnt reporter activity in the midbrain (green arrow) upon LiCl treatment, relative to controls. hpf, hours post-fertilization. mb, midbrain. Scale bar, 50 µm.
Overall, the results show that Wnt signaling is sensitized in cohesin-mutant zebrafish embryos, similar to our observations in MCF10A and CMK cell lines. Our findings agree with previous work showing that canonical Wnt signaling is hyperactivated in cohesin-loader nipblb-loss-of-function zebrafish embryos (Mazzola et al., 2019). Altogether, the results suggest that hyperactivation of Wnt signaling is a conserved feature of cohesin-deficient cells, and that enhanced sensitivity to Wnt is at least in part due to stabilization of β-catenin.
Discussion
The cohesin complex and its regulators are encoded by several different loci, and genetic alterations in any one of them may occur in up to 26% of patients included in The Cancer Genome Atlas (TCGA) studies (Romero-Pérez et al., 2019). Therefore, we were motivated to identify compounds that interfere with cell viability in more than just one type of cohesin mutant. The generation of RAD21, SMC3, and STAG2 cohesin mutations in the breast epithelial cell line MCF10A resulted in mild cell cycle and nucleolar phenotypes that are consistent with the many cellular roles of cohesin. Differences between RAD21 and SMC3 heterozygotes vs STAG2 homozygotes could be explained by the fact that RAD21 and SMC3 are obligate members of the cohesin ring, whereas STAG2 can be compensated by STAG1.
Synthetic lethal sensitivities of cohesin-mutant cells that emerged from our screen mostly related to the phenotypes of cohesin-mutant cells, and/or their previously identified vulnerabilities. For example, cohesin-mutant cells were sensitive to inhibitors of the PI3K/Akt/mTOR pathway that feeds into ribosome biogenesis, epigenetic inhibitors that could interfere with cohesin’s genome organization and gene expression roles, and a limited sensitivity to PARP inhibitors. The observed sensitivity to PI3K/AKT/mTOR inhibitors is consistent with the nucleolar disruption present in cohesin-deficient cell lines, and with cohesin’s known involvement in rDNA transcription and ribosome biogenesis (Bose et al., 2012). We have also previously described evidence for sensitivity of cohesin-deficient cell lines to bromodomain (BET) inhibition (Antony et al., 2020). However, the Wnt agonist LY2090314, which mimics Wnt activation by inhibiting GSK3-β (Atkinson et al., 2015), emerged as a novel class of compound that inhibited the growth of all three cohesin mutants tested. Because there is pathway convergence between Wnt signaling and the PI3K/AKT/mTOR pathway (Prossomariti et al., 2020) it is possible that these pathways represent a common avenue of sensitivity. The identification of PI3K/AKT/mTOR inhibitors in the screen supports this notion.
Wnt signals are important for stem cell maintenance and renewal in multiple mammalian tissues (Clevers et al., 2014). Canonical Wnt signaling is required for self-renewal of leukemia stem cells (LSCs) (Kang et al., 2020), and is reactivated in more differentiated granulocyte-macrophage progenitors (GMPs) when they give rise to LSCs (Wang et al., 2010). Wnt signaling is an important regulator of hematopoietic stem cell (HSC) self-renewal as well (Rattis et al., 2004). Our results suggest that cohesin-deficient cells are sensitive to Wnt agonism because they already have stabilization of β-catenin. Interestingly, multiple experimental systems showed that dose-dependent reduction in cohesin function causes HSC expansion accompanied by a block in differentiation (Viny and Levine, 2018; Mazumdar and Majeti, 2017). It is possible that the effects of cohesin deficiency on HSC development could be in part due to enhanced Wnt signaling.
Cohesin mutations are particularly frequent (53%) in Down Syndrome-associated Acute Megakaryoblastic Leukemia (DS-AMKL) patients (Yoshida et al., 2013). In the edited DS-AMKL cell line STAG2-CMK, there was a dramatically enhanced immediate early transcriptional response to WNT3A, including induction of hundreds of genes that are not normally responsive to WNT3A in this cell line. It is possible that STAG2 deficiency in CMK leads to an altered chromatin structure that sensitizes genes to Wnt signaling. However, we observed increased nuclear localization of β-catenin in STAG2-CMK, increasing the likelihood that stabilization of β-catenin in STAG2-CMK cells accounts at least in part for amplified gene expression in response to WNT3A. Interestingly, DS-AMKL leukemias are often associated with amplified Wnt signaling caused by Trisomy 21 (Emmrich et al., 2014) suggesting a possible synergy between cohesin mutations and Wnt signaling dysregulation in DS-AMKL.
We observed an enhanced Wnt reporter response to Wnt activation in rad21- and stag2b-mutant zebrafish, arguing that Wnt sensitivity upon cohesin mutation is a conserved phenomenon. Previous research shows that cohesin genes are at once both targets (Xu et al., 2014; Ghiselli et al., 2003) and upstream regulators (Estarás et al., 2015; Schuster et al., 2015; Avagliano et al., 2017; Mazzola et al., 2019) of Wnt signaling. For example, depletion of cohesin-loader nipbl in zebrafish embryos was reported to downregulate Wnt signaling at 24 hpf (Pistocchi et al., 2013), but to upregulate it at 48 hpf (Mazzola et al., 2019). What determines the directionality of cohesin’s effect on Wnt signaling is unclear. It is possible that feedback loops operate differently in cell- and signal-dependent contexts (MacDonald et al., 2009) in cohesin-deficient backgrounds.
Is Wnt pathway activation a conserved feature of cohesin-mutant cancers? Using publicly available data at TCGA (Grossman et al., 2016), we correlated expression of genes in the Wnt pathway (Hallmark WNT β-catenin signaling) with nonsense mutation of STAG2, the most common of the cohesin mutations, in the four most represented cancers (bladder cancer, endometrial carcinoma, glioblastoma multiforme and cervical kidney renal papillary cell carcinoma). In comparative ranking, the Wnt pathway ranks 8th for enrichment in upregulated genes (after DNA repair/G2 checkpoint, MYC/E2F targets, oxidative phosphorylation and the unfolded protein response; Supplementary file 2). While our analysis could indicate conserved upregulation of the Wnt pathway in STAG2 mutant cancers, a caveat is that the true situation is likely to be more complex. If β-catenin is stabilized, upstream and in-parallel effectors such as Wnt ligands, receptors, signaling intermediates might trend to downregulation owing to negative feedback loop regulation of the Wnt pathway for example via DKK1 (Niida et al., 2004). Furthermore, because cancer collections are unlikely to be under Wnt stimulation at the time of RNA collection, it is possible that their true Wnt sensitivity will be missed in a steady-state transcription analysis.
A remaining question is exactly what causes β-catenin stabilization and Wnt hyperactivation in cohesin-mutant cells. One possibility is that stabilization of β-catenin upon cohesin deficiency is linked to energy metabolism. For example, in developing amniote embryos, glycolysis in actively growing cells of the embryonic tail bud raises the intracellular pH, which in turn causes β-catenin stabilization (Oginuma et al., 2020; Oginuma et al., 2017). This mode of glycolysis is similar to that observed in cancer cells and is known as the Warburg effect (Parks et al., 2013; Wang et al., 2016). Metabolic and oxidative stress disturbances in cohesin-mutant cells observed in this study and others (Bose et al., 2012; Cukrov et al., 2018) might similarly influence β-catenin stability and in turn, sensitivity to incoming Wnt signals. Further research will be needed to determine how β-catenin is stabilized in cohesin-deficient cells, whether β-catenin stabilization is part of a transition to malignancy, and if the associated Wnt sensitivity represents a therapeutic opportunity in cohesin-mutant cancers.
Materials and methods
Cell culture
Request a detailed protocolMCF10A (a spontaneously immortalized breast epithelial cell line) was purchased from Sigma (product #: CRL 10317). Sigma use STR analysis to verify the line. A PCR test in our laboratory showed that all cells were mycoplasma free. CMK (acute megakaryocytic leukemia associated with Down Syndrome) was a gift from Dr Motomi Osato, National University of Singapore. These cells tested negative with MycoAlert Mycoplasma Detection Kit in our laboratory. K562 (chronic myelogenous leukemia) were purchased from ATCC, CCL-243. These cells tested mycoplasma free using B3903-Mycoplasma Detection Kit - QuickTest-com from Biotool. The HTC116 cell line was purchased form ATCC. STR analysis is used by ATCC to verify all cell lines.
MCF10A and its cohesin-deficient derivatives were maintained in Dulbecco’s modified Eagle medium (DMEM) (Life Technologies) supplemented with 5% horse serum (Life Technologies), 20 ng/mL of epidermal growth factor (EGF) (Peprotech), 0.5 mg/mL of hydrocortisone, 100 ng/mL of cholera toxin and 10 μg/mL of insulin (local pharmacy). All supplements were purchased from Sigma-Aldrich unless otherwise stated. K562 cells were maintained in Isocove’s Modified Dulbecco’s Media (IMDM) (Life Technologies) containing 10% fetal bovine serum. CMK cells were maintained in RPMI 1640 media containing 10% fetal bovine serum. CMK cells and adherent K562-STAG2 null line were detached for subculture using trypsin-EDTA (0.005% final concentration, Life Technologies). For WNT stimulation, recombinant Human WNT-3A Protein (5036-WN, R and D systems) was used at 200 ng/mL for 4 hr. Human colorectal carcinoma HCT116 cells were grown in DMEM with 10% fetal calf serum and antibiotics. Stable transfection of HCT116 cells with SMC1A mutations c.2027A > G and c.2479 C > T was described previously (Sarogni et al., 2019). All cells were cultured at 37°C in 5% CO2.
Generation of isogenic cell lines using CRISPR-CAS9 editing
Request a detailed protocolWe used CRISPR-CAS9 sgRNAs targeting the 5' and 3' UTR regions of RAD21, SMC3 and STAG2 gene to create MCF10A RAD21+/-, MCF10A SMC3+/-, and MCF10A STAG2-/-. Specific sgRNAs were cloned into the px458 plasmid (Ran et al., 2013) and transfected into MCF10A cells. Single GFP-positive cells were isolated into 96-well plates using a FACSAriaII (Becton Dickinson) and clonally expanded. PCR screening and Sanger sequencing identified cells with heterozygous deletion of RAD21 or SMC3 and homozygous deletion of STAG2. Primer and guide sequences are provided in Supplementary file 1. Editing of K562 to create STAG2 R614* mutation has been described previously (Antony et al., 2020). The CMK line with the STAG2 R614* mutation was generated using the same sgRNA and method.
Quantitative PCR (RT-qPCR)
Request a detailed protocolTotal RNA was extracted using NucleoSpin RNA kit (Machery-Nagel). cDNA was synthesized using qScript cDNA SuperMix (Quanta Biosciences). RT-qPCR was performed on a LightCycler 480 II (Roche Life Science) using SYBR Premix Ex Taq (Takara). Expression values relative to reference genes cyclophilin and glyceraldehyde 3-phosphate dehydrogenase (GAPDH) were derived using qBase Plus (Biogazelle). Primer sequences are provided in Supplementary file 1.
Antibodies
Primary antibodies used are as follows: anti-RAD21 (ab992), anti-SMC3 (#5696, CST), anti-STAG2 (ab4463), anti-γ-Tubulin (T5326; Sigma-Aldrich) in 1:5000, anti-fibrillarin (ab5821), anti-nucleolin (ab13541), anti-gamma H2AX (ab26350), anti-TP53 (ab131442), anti-β-catenin (#9562, CST), anti-phospho-β-catenin (Ser675) (#9567, CST), anti-phospho- phospho-β-catenin (Ser33/37/Thr41) (#9561, CST). Secondary antibodies used for immunofluorescence were anti-mouse Alexa Fluor 488 (1:2000, Life Technologies), anti-rabbit Alexa Fluor 488 (1:2000, Life Technologies), anti-rabbit Alexa Fluor 568 (1:2000, Life Technologies). All antibodies were used in 1:1000 for immunoblotting or immunofluorescence unless otherwise stated.
Immunoblotting
Request a detailed protocolImmunoblotting was performed as described previously (Antony et al., 2015; Dasgupta et al., 2016). Primary antibodies were detected using IRDye 800CW Donkey anti-Goat IgG and IRDye 680CW Goat anti-mouse IgG, IRDye 800CW Goat anti-rabbit IgG (LICOR). LI-COR Odyssey and LI-COR Image Studio software was used to image and quantify blots.
Proliferation assays
Request a detailed protocolMCF10A parental and isogenic cohesin-deficient cell lines were seeded in 96-well plates at 2000 cells per well. Cell confluence was monitored by time-lapse microscopy using IncuCyte FLR with a 10X objective for 5 days.
Cell cycle analysis
Request a detailed protocolCells were synchronized using double thymidine block as described previously (Dasgupta et al., 2016). Samples were harvested, fixed, and stained with 10 μg/mL propidium iodide (PI) and 250 μg/mL RNase A, 37°C for 30 min. Cells were then analyzed using a Beckman Coulter Gallios Flow Cytometer.
Immunofluorescence and imaging
Request a detailed protocolCells stained for FBL, gH2AX, and TP53 were imaged using an Opera Phenix high content screening system, with 63x water objectives in confocal mode. Spot enumeration and signal intensities were analyzed using Harmony software (PerkinElmer). CMK or MCF10A cells were fixed with 4% (v/v) paraformaldehyde in PBS for 10 min, then permeablized with 0.1% Triton in PBS. CMK cells were spun onto slides using SHANDON cytospin prior to fixation. Cells were blocked with 2% (w/v) bovine serum albumin in PBS, then incubated with primary antibody overnight at 4°C. The next day, cells were washed with PBS and incubated with secondary antibody and Hoechst 33342 (1 µg/mL) at room temperature for 1 hr. Cells were then washed and mounted on slides using ProLong Gold Anti-fade (Life Technologies) or DAKO mounting medium. Imaging was done using a Nikon C2 confocal microscope with NIS Elements software and images were processed and quantified using Image J or FIJI software.
High-throughput compound screen
Request a detailed protocol3009 compounds from an FDA-approved, kinase inhibitors and epigenetic libraries (Compounds Australia, Griffith University, Australia) were screened against MCF10A parental and MCF10A cohesin-deficient clones. Compounds are stored by Compounds Australia under robust environmental conditions and supplied in assay ready plate format. MCF10A parental cells and cohesin-deficient derivatives were screened as single biological replicates with two technical replicates for each drug concentration. MCF10A parental cells were seeded at 600 cells/well, RAD21+/- at 800 cells/well, SMC3+/- at 1200 cells/well and STAG2-/- at 700 cells/well, into 384-well CellCarrier-384 Ultra microplates (PerkinElmer) with EL406 microplate washer dispenser (BioTek). 24 hr later, growth media was aspirated from the plates and 35 µL of fresh MCF10A complete growth media was added into each well using BioTek EL406TM Microplate washer dispenser. 5 µL of 8x concentrated compound solution (diluted in MCF10A complete growth media) was added into each well with a JANUS automated liquid handling system (PerkinElmer). At the time of drug treatment, one untreated plate was retained to determine the cell number at t = 0. For the rest of the assay plate, compounds were added in 40 µL of complete medium per well using an Echo 550 Liquid Handler (Labcyte). Positive (Campthothecin, 40 nM) and negative controls (DMSO 0.1%) were added to each assay plate for quality control. Duplicated control plates of each cell line were also prepared for cell count quantification prior to the addition of drugs, to be used in the growth rate inhibition (GR) metrics. After 48 hr, control and treated cells in were fixed and stained simultaneously using 4% PFA/0.5 µg/mL DAPI/0.1% Triton X-100 solution. Cells were imaged using an Opera Phenix high content screening system and analyzed using Columbus software (PerkinElmer) using 20x water objectives. To determine compound activity, dose-response curves for each compound was plotted and using an R package, GRmetrics (Clark et al., 2017). Synthetic lethal candidate compounds were selected and ranked based on the differential area over curve (AOC) metrics (Hafner et al., 2016) derived from dose-response curves of MCF10A parental and cohesin-deficient cell lines. Compounds that caused ≥30% growth inhibition (AOC ≥0.15) in cohesin-deficient clones compared with parental MCF10A cells were selected as hits. A secondary screen was performed with 85 candidate hits identified from primary screen (two technical replicates per clone per concentration). The secondary screen was added to independently validate primary screen data. Activity of the selected compounds was tested in eleven concentrations (0.5 nM to 10 µM). Secondary screen hit compounds were selected based on the same threshold used in primary screen.
Cell viability assays to validate compounds identified from the screen
Request a detailed protocolIndividual compounds were purchased from Sigma-Aldrich, Selleck Chem or MedChemExpress and dissolved in DMSO at recommended concentrations. MCF10A parental cells were seeded in 96-well plates at 3000 cells/well, RAD21+/- at 4500 cells/well, SMC3+/- at 5000 cells/well and STAG2-/- at 3500 cells/well. Cells were incubated for 24 hr then treated with compounds for 48 hr. DAPI-stained cells were counted using a Lionheart FX automated microscope (BioTek). K562 parental and STAG2-/- cells were seeded at 2000 cells per well in 96 well plates, incubated with the compound for 48 hr after which viability was measured using 3-(4,5-dimethylthiazol-2-yl)−2,5-diphenyltetrazolium bromide or MTT.
RNA-sequencing and analyses
Request a detailed protocolTotal RNA was extracted using NucleoSpin RNA kit (Machery-Nagel). MCF10A libraries from three biological replicates of each cell type were prepared using NEBNext Ultra RNA Library Prep Kit (Illumina) and sequenced on HiSeq X by Annoroad Gene Technology Ltd. (Beijing, China), contracted through Custom Science (NZ). CMK lines were treated with 200 ng/mL of recombinant human WNT3A (R and D systems) for 4 hr. Libraries were prepared from baseline (non-treated) and WNT3A treated CMK cell lines using Illumina TruSeq stranded mRNA library and sequenced on the HiSeq 2500 V4 at the Otago Genomics Facility (Dunedin, New Zealand). RNA-sequencing reads were first trimmed for sequencing adapters and low quality (q < 20). Cleaned reads were then aligned to the human genome GRCh37 (hg19) using HISAT2 version 2.0.5 with gene annotation from Ensembl version 75. Read counts were retrieved by exon and summarized by gene using featureCount (Liao et al., 2014) version v1.5.3. Differentially expressed genes in the STAG2-/- mutants versus wild type were identified using DESeq2 (Love et al., 2014). P-values were adjusted for multi-test following Benjamini-Hochberg methodology. Pathways analyses were performed using Molecular Signature Databases ReactomePA (Yu and He, 2016) and clusterProfiler (Yu et al., 2012) on differentially expressed genes.
Zebrafish methods and imaging
Request a detailed protocolWild type (WIK), TCF reporter line Tg(7xTCF-Xla.Siam:nslmCherry)ia5 (Moro et al., 2012), stag2bnz207 (Ketharnathan et al., 2020) and rad21nz171 (Horsfield et al., 2007) mutant fish lines were maintained according to established husbandry methods (Westerfield, 1995). We crossed heterozygous rad21nz171 or homozygous stag2bnz207 mutants to the TCF reporter line then crossed these fish (rad21nz171/+;TCF/+) to rad21nz171/+, or (stag2bnz207/+;TCF/+) to stag2bnz207/+. Progeny carrying the TCF reporter were incubated with 0.15 M lithium chloride (LiCl) from 4 hpf for 16 hr. LiCl salt was directly dissolved in embryo water and added to embryos in six well plates. 50 embryos were used per treatment group. At 20 hpf, embryos were washed, anesthetized with MS-222 (200 mg/L) and mounted in low melting agarose (0.6%) or 3% methyl cellulose for imaging. Z-stacks were acquired using a Nikon C2 confocal microscope. Embryos were genotyped post-imaging.
Statistical analyses
Request a detailed protocolAll statistical analyses were carried out using R Studio or Prism version eight software (GraphPad).
Data availability
All RNA sequencing data has been deposited at the GEO database under accession codes GSE154086. All data generated or analysed during this study are included in the manuscript and supporting files. Source data files have been provided for Figures 1-5 and Table 1.
-
NCBI Gene Expression OmnibusID GSE154086. Expression profiling in cohesin mutant MCF10A epithelial and CMK leukaemia cells.
References
-
Cohesin modulates transcription of estrogen-responsive genesBiochimica Et Biophysica Acta (BBA) - Gene Regulatory Mechanisms 1849:257–269.https://doi.org/10.1016/j.bbagrm.2014.12.011
-
BET inhibition prevents aberrant RUNX1 and ERG transcription in STAG2 mutant leukaemia cellsJournal of Molecular Cell Biology 12:397–399.https://doi.org/10.1093/jmcb/mjz114
-
Glioblastoma cells containing mutations in the cohesin component STAG2 are sensitive to PARP inhibitionMolecular Cancer Therapeutics 13:724–732.https://doi.org/10.1158/1535-7163.MCT-13-0749
-
Mutant cohesin drives chromosomal instability in early colorectal adenomasHuman Molecular Genetics 23:6773–6778.https://doi.org/10.1093/hmg/ddu394
-
HDAC8 inhibition blocks SMC3 deacetylation and delays cell cycle progression without affecting Cohesin-dependent transcription in MCF7 Cancer cellsJournal of Biological Chemistry 291:12761–12770.https://doi.org/10.1074/jbc.M115.704627
-
Cohesin mutations in cancerCold Spring Harbor Perspectives in Medicine 6:a026476.https://doi.org/10.1101/cshperspect.a026476
-
Genetic alterations in members of the wnt pathway in acute leukemiaLeukemia & Lymphoma 53:508–510.https://doi.org/10.3109/10428194.2011.613133
-
The cohesin SMC3 is a target the for β-Catenin/TCF4 transactivation pathwayJournal of Biological Chemistry 278:20259–20267.https://doi.org/10.1074/jbc.M209511200
-
Toward a shared vision for Cancer genomic dataNew England Journal of Medicine 375:1109–1112.https://doi.org/10.1056/NEJMp1607591
-
Cohesion promotes nucleolar structure and functionMolecular Biology of the Cell 25:337–346.https://doi.org/10.1091/mbc.e13-07-0377
-
Cohesin mutations in human CancerBiochimica et Biophysica Acta 1866:1–11.https://doi.org/10.1016/j.bbcan.2016.05.002
-
Cohesin-dependent regulation of runx genesDevelopment 134:2639–2649.https://doi.org/10.1242/dev.002485
-
Deregulated notch and wnt signaling activates early-stage myeloid regeneration pathways in leukemiaJournal of Experimental Medicine 217:20190787.https://doi.org/10.1084/jem.20190787
-
Cohesin components Stag1 and Stag2 differentially influence haematopoietic mesoderm development in zebrafish embryosFrontiers in Cell and Developmental Biology 8:617545.https://doi.org/10.3389/fcell.2020.617545
-
Somatic mutation of the cohesin complex subunit confers therapeutic vulnerabilities in CancerJournal of Clinical Investigation 128:2951–2965.https://doi.org/10.1172/JCI98727
-
The wnt signaling pathway in development and diseaseAnnual Review of Cell and Developmental Biology 20:781–810.https://doi.org/10.1146/annurev.cellbio.20.010403.113126
-
Cohesin biology: from passive rings to molecular motorsTrends in Genetics 36:387–389.https://doi.org/10.1016/j.tig.2020.03.001
-
The role of mutations in the cohesin complex in acute myeloid leukemiaInternational Journal of Hematology 105:31–36.https://doi.org/10.1007/s12185-016-2119-7
-
NIPBL: a new player in myeloid cell differentiationHaematologica 104:1332–1341.https://doi.org/10.3324/haematol.2018.200899
-
Cohesin loss alters adult hematopoietic stem cell homeostasis, leading to myeloproliferative neoplasmsJournal of Experimental Medicine 212:1833–1850.https://doi.org/10.1084/jem.20151323
-
Synthetic lethality and Cancer: cohesin and PARP at the replication forkTrends in Genetics 29:290–297.https://doi.org/10.1016/j.tig.2012.12.004
-
Cohesin SA1 and SA2 are RNA binding proteins that localize to RNA containing regions on DNANucleic Acids Research 48:5639–5655.https://doi.org/10.1093/nar/gkaa284
-
Genomic classification and prognosis in acute myeloid leukemiaNew England Journal of Medicine 374:2209–2221.https://doi.org/10.1056/NEJMoa1516192
-
Disrupting proton dynamics and energy metabolism for Cancer therapyNature Reviews Cancer 13:611–623.https://doi.org/10.1038/nrc3579
-
Are wnt/β-Catenin and PI3K/AKT/mTORC1 distinct pathways in colorectal Cancer?Cellular and Molecular Gastroenterology and Hepatology 10:491–506.https://doi.org/10.1016/j.jcmgh.2020.04.007
-
Genome engineering using the CRISPR-Cas9 systemNature Protocols 8:2281–2308.https://doi.org/10.1038/nprot.2013.143
-
Wnt signaling in the stem cell nicheCurrent Opinion in Hematology 11:88–94.https://doi.org/10.1097/01.moh.0000133649.61121.ec
-
Organizational principles of 3D genome architectureNature Reviews Genetics 19:789–800.https://doi.org/10.1038/s41576-018-0060-8
-
Overexpression of the cohesin-core subunit SMC1A contributes to colorectal Cancer developmentJournal of Experimental & Clinical Cancer Research 38:108.https://doi.org/10.1186/s13046-019-1116-0
-
A neural crest origin for cohesinopathy heart defectsHuman Molecular Genetics 24:7005–7016.https://doi.org/10.1093/hmg/ddv402
-
S-phase and DNA damage activated establishment of sister chromatid cohesion--importance for DNA repairExperimental Cell Research 316:1445–1453.https://doi.org/10.1016/j.yexcr.2009.12.018
-
Ultrastructural and immunocytochemical characterization of an immortalized human breast epithelial cell line, MCF-10Cancer Research 50:6087–6094.
-
Dose-dependent role of the cohesin complex in normal and malignant hematopoiesisJournal of Experimental Medicine 212:1819–1832.https://doi.org/10.1084/jem.20151317
-
Cohesin mutations in myeloid malignancies made simpleCurrent Opinion in Hematology 25:61–66.https://doi.org/10.1097/MOH.0000000000000405
-
Emerging themes in cohesin Cancer biologyNature Reviews Cancer 20:504–515.https://doi.org/10.1038/s41568-020-0270-1
-
Resolving the genomic localization of the kollerin Cohesin-Loader complexMethods in Molecular Biology 1515:115–123.https://doi.org/10.1007/978-1-4939-6545-8_7
-
BookThe Zebrafish Book. a Guide for the Laboratory Use of Zebrafish (Brachydanio rerio)University of Oregon Press.
-
The landscape of somatic mutations in down syndrome–related myeloid disordersNature Genetics 45:1293–1299.https://doi.org/10.1038/ng.2759
-
clusterProfiler: an R package for comparing biological themes among gene clustersOMICS: A Journal of Integrative Biology 16:284–287.https://doi.org/10.1089/omi.2011.0118
-
ReactomePA: an R/Bioconductor package for reactome pathway analysis and visualizationMolecular BioSystems 12:477–479.https://doi.org/10.1039/C5MB00663E
-
PARP inhibitor development for systemic Cancer targetingAnti-Cancer Agents in Medicinal Chemistry 7:515–523.https://doi.org/10.2174/187152007781668715
Decision letter
-
Richard M WhiteSenior Editor; Memorial Sloan Kettering Cancer Center, United States
-
Ravi MajetiReviewing Editor; Stanford University, United States
-
Aaron VinyReviewer
-
Massa ValentinaReviewer; Università degli Studi di Milano, Italy
In the interests of transparency, eLife publishes the most substantive revision requests and the accompanying author responses.
Acceptance summary:
The manuscript describes an interesting finding that stimulation of Wnt signaling can lead to lethality in cohesin mutant cells. This suggests a potentially novel therapeutic strategy for cohesion-mutant cancers.
Decision letter after peer review:
Thank you for submitting your article "Cohesin mutations are synthetic lethal with stimulation of WNT signaling" for consideration by eLife. Your article has been reviewed by three peer reviewers, one of whom is a member of our Board of Reviewing Editors, and the evaluation has been overseen by Richard White as the Senior Editor. The following individuals involved in review of your submission have agreed to reveal their identity: Aaron Viny (Reviewer #2); Massa Valentina (Reviewer #3).
The reviewers have discussed the reviews with one another and the Reviewing Editor has drafted this decision to help you prepare a revised submission.
We would like to draw your attention to changes in our revision policy that we have made in response to COVID-19 (https://elifesciences.org/articles/57162). Specifically, when editors judge that a submitted work as a whole belongs in eLife but that some conclusions require a modest amount of additional new data, as they do with your paper, we are asking that the manuscript be revised to either limit claims to those supported by data in hand, or to explicitly state that the relevant conclusions require additional supporting data.
Our expectation is that the authors will eventually carry out the additional experiments and report on how they affect the relevant conclusions either in a preprint on bioRxiv or medRxiv, or if appropriate, as a Research Advance in eLife, either of which would be linked to the original paper.
Summary:
In this manuscript, Chin et al., investigate synthetic lethal candidates in conjunction with cohesin mutations in cancer. The authors engineer mutant cell line variants for several members of the cohesin complex and demonstrate that they exhibit established features of cohesin-mutant cells. Then they conduct a small molecule screen to identify compounds and pathways with increased sensitivity in these cohesin-mutant cells. They identify a GSK3 inhibitor suggesting that Wnt pathway activation is synthetic lethal with cohesin mutants. They show evidence suggesting this occurs though β-catenin stabilization and activation of Wnt target genes. Finally, they use a genetic zebrafish model to show hyperactivation of Wnt pathway in Rad21-deficient embryos. In general, the link between cohesin mutants and Wnt pathway activation is of interest, but there are several issues with the manuscript.
Essential revisions:
1) While there is evidence presented here that cohesin-mutant cells exhibit increased Wnt pathway activation, there is no data demonstrating that it is key to the cancer phenotypes of cohesin-mutant cancers. Such data would enhance the significance of the findings here.
2) Analysis of the cohesin-mutant MCF10A clones relied on investigation of a single clone per gene. This is not sufficient to draw convincing conclusions, particularly with CRIPSR/Cas9 targeting and potential off target effects. Several clones per gene should be investigated to account for potential clonal variability, particularly given the differences in phenotypes observed between these clones.
3) The chemical screen identified a number of cohesin mutant sensitivities, and the common hits are shown, however, given the phenotypic differences between the RAD21/SMC3 clone and the STAG2 null clones in micronuclei and lagging strands, did the 76 STAG2 unique chemical compounds share a functional class? Why do the authors think SMC3 and STAG2 had far more overlap that each of these did with RAD21? Why did RAD21 have so few hits overall by comparison?
4) Western blotting shows cohesin mutant clones have increased total and membrane β-catenin that lead to transcriptional responses. Is there evidence of Wnt pathway/β-catenin transcriptional activation in publicly available datasets of cohesin perturbation or cohesin-mutant cancers?
5) The zebrafish experiments are fundamental and elegantly done. However, I think that, given the role of RAD21 and of the other cohesin subunits, it would be pivotal to include similar experiments with zebrafish modelling other genes of the complex.
6) Cohesin subunit deficient clones of MCF10A were generated using CRISPR-Cas9 editing. While this is clearly described in the Materials and methods section, it would be helpful to state the model system as well in results.
7) For cell cycle analysis (Results section) it would be extremely important, also considering the literature, to assess cell death rate.
8) In Figure 1F/1G/Figure 1—figure supplement 1B what do the square symbols represent? Please include in the legend.
9) In Figure 4A/B/Figure 4—figure supplement 1A/1B please include size markers. Also Figure 4—figure supplement 1 needs a loading control such as tubulin.
https://doi.org/10.7554/eLife.61405.sa1Author response
Essential revisions:
1) While there is evidence presented here that cohesin-mutant cells exhibit increased Wnt pathway activation, there is no data demonstrating that it is key to the cancer phenotypes of cohesin-mutant cancers. Such data would enhance the significance of the findings here.
The reviewers make a good point. It would be ideal to test if cohesin mutant primary cancer cells are Wnt-sensitivity or hyperactivity as part of their phenotype, but actually testing this in vivo is quite a difficult task. Primary cells with and without cohesin mutations would need to be compared in a Wnt-stimulated scenario, because there may be no apparent phenotype to measure in cells at baseline or endpoint conditions where Wnt signalling is not active.
While we unfortunately do not have access to primary cells, in a revision experiment for this paper we managed to show that Wnt sensitivity is extendable to another cancer cell type (colorectal cancer cell line HCT116), and a different cohesin subunit mutation: SMC1A.
To provide the new data, we collaborated with our Italian colleague Antonio Musio, whose lab has the appropriate models. Antonio’s group have characterised two cohesin SMC1A mutations identified in human colorectal carcinomas, c.2027A>G and c.2479 C>T (Cucco et al., 2014; Sarogni et al., 2019). Using the same treatment conditions as we did for MCF10A cells, Antonio’s group showed increased b-catenin protein levels, particularly in response to LY2090314, in HCT116 cells that stably overexpress SMC1A containing these two different mutations. Importantly, SMC1A is a cohesin gene that wasn’t tested in our screen or subjected to downstream analysis in our lab, and yet these experiments confirm our independent findings.
We have added a new figure supplement with these data (Figure 5—figure supplement 2). Antonio and colleague Michela Pallotta (who performed the work) have been included as co-authors on the manuscript. The new text in the results (which also explains the model’s limitation) is reproduced below:
“To determine if stabilization of β-catenin is conserved in a second model of cohesin mutant cancer, we performed an identical LY2090314 treatment on HCT116 cells that were stably transfected to express two SMC1A mutations identified in human colorectal carcinomas, c.2027A>G (leading to p.E676G change near the hinge domain) and c.2479 C>T (leading to a truncated protein, p.Q827X) [53, 54]. The limitation is that these cells are not a model of cohesin deficiency, but rather, one in which normal cohesin function is perturbed by expression of a mutant version of SMC1A [54]. We observed an increased basal level of phosphorylated β-catenin at Ser675 in cells that express either of these SMC1A mutants. Furthermore, LY2090314 treatment markedly increased total β-catenin in both the SMC1A mutants compared with HCT116 wild type controls (Figure 5—figure supplement 2A,B). The results indicate that abnormally high levels of active β-catenin are also present following Wnt stimulation when a fourth subunit of cohesin, SMC1A, is perturbed.”
Interestingly, the mutant SMC1A-expressing HCT116 cells are more aggressive in tumour formation (Sarogni et al., 2019). Future experiments beyond the scope of this study could follow up to determine if Wnt signalling is important for this phenotype.
Below in response to point #4, we describe an additional analysis of TCGA data, which supports the idea that cohesin mutant cancers have over-active Wnt signalling.
2) Analysis of the cohesin-mutant MCF10A clones relied on investigation of a single clone per gene. This is not sufficient to draw convincing conclusions, particularly with CRIPSR/Cas9 targeting and potential off target effects. Several clones per gene should be investigated to account for potential clonal variability, particularly given the differences in phenotypes observed between these clones.
Owing to the expense of high throughput screening, it was too costly to screen more than one clone per gene in the screen. However, even though there’s only a single clone representing each cohesin subunit in the SL screen, we believe the screen was actually quite conservative and robust, because the top hits we followed up were only selected if they affected all three cohesin mutant clones.
However, in response to this valid comment, we have now analysed one additional RAD21 and one addition SMC3 clone that were produced, but not used in the screen. These clones (RAD21+/- clone 2 and SMC3 +/- clone 2) behave similarly to the original clones that were used in the screen. Data showing allele details, cell morphology, mRNA levels, chromosome cohesion defects, cell viability, growth curves and response to LY2090314 are shown in Figure 3—figure supplement 3. The cells grow slightly more slowly than wild type and are differentially sensitive to LY2090314. These results indicate that the clones used for the screen were likely to be typical and any variability was not likely to be due to off-target effects. The new data are described in the Results section:
“Collectively, characterization of our cohesin-deficient clones provided confidence that they represent suitable models for synthetic lethal screening. To confirm that the phenotypes of our chosen clones are representative, we selected a further two clones with heterozygous deletions in SMC3 and RAD21, and monitored their growth, morphology and chromosome cohesion (Figure 3—figure supplement 3). These analyses showed that the two additional deletion clones were similar to those already characterized (Figure 1, Figure 2, Figure 3), providing surety that our cohesin-deficient cell lines have properly representative phenotypes. Furthermore, none of the cohesin-deficient clones exhibited enhanced cell death compared with parental MCF10A cells (Figure 2B; Figure 3—figure supplement 3).”
Regarding the subtly different phenotypes of the cohesin-deficient clones, it is perhaps not surprising that the clones do not have identical phenotypes, particularly for the STAG2 homozygous deletion.
First, STAG2 can be substituted by STAG1 and therefore a STAG2 mutation can be homozygous, unlike SMC3 and RAD21. This means that in the STAG2 mutant, STAG1 cohesin is predominant, whereas in RAD21 and SMC3 mutants, there would be predicted to be a reduction in dose of STAG1- and STAG2-containing cohesin.
Second, differences in chromosome segregation phenotypes between the different cohesin-deficient clones might be due to STAG2 having a role in centromere chromosome cohesion (Canudas and Smith, 2009) such that its homozygous loss results in “cleaner” mis-segregation events and more frequent complete loss of chromosome cohesion, as shown in Figure 2D (please note that the key of Figure 2D has been corrected to show that the black section of the bar is 'normal' and the 'green' sections correspond to loss of cohesion!).
The functional differences between STAG2 and STAG1 cohesin are described nicely in a recent review from Ana Losada and colleagues (Cuadrado and Losada, 2020).
In the screen there was quite a high degree of overlap of hits for the clones we generated, particularly SMC3 and RAD21, where only 5 of 27 compounds affecting RAD21 mutants were not also a hit in SMC3 mutants, and only 1 hit that was exclusive to RAD21. This suggests that despite potential differences between the clones, there were a number of common sensitivities. All clones have been well characterised, and findings extended to other cell lines and zebrafish.
3) The chemical screen identified a number of cohesin mutant sensitivities, and the common hits are shown, however, given the phenotypic differences between the RAD21/SMC3 clone and the STAG2 null clones in micronuclei and lagging strands, did the 76 STAG2 unique chemical compounds share a functional class?
The unique hits for each clone are given in Figure 4—figure supplement 2. This figure shows that functional classes are largely shared among the cohesin subunits. From this figure, classes of compound that affected only cells containing the STAG2 mutation can be identified. These include “cytoskeletal”, “GPCR and G protein”, ‘Immunology and inflammation’ and “transmembrane transporters”. These extra categories specific to STAG2 mutation look like they reflect a response to particular signalling pathways, and could result from the more specialised role of STAG2-cohesin in tissue-specific transcription.
Why do the authors think SMC3 and STAG2 had far more overlap that each of these did with RAD21?
Figure 4E shows that a number of compounds affecting RAD21 mutant cells also inhibited SMC3 mutants (n=4) and STAG2 mutants (n=4), and there were 18 hits in common with all 3 mutants. In addition, RAD21 mutant cells were sensitive to only one compound that didn’t also affect the other mutants, and this particular compound is involved in DNA damage repair. It’s just that there were fewer compounds overall that affected RAD21 mutants. That’s probably why there wasn’t as much overlap as observed with SMC3 and STAG2.
Why did RAD21 have so few hits overall by comparison?
That’s a very good question. We don’t have an answer but can speculate that additional compound sensitivities in SMC3 and STAG2 mutants might derive from important dose-dependent functions that they have in the nucleus, which are not equally sensitised by limiting the amount of RAD21-containing cohesin. We wonder if SMC3 depletion exposes an important role in chromosome maintenance because of our observation of increased nuclear decompaction and micronuclei the SMC3 mutant clone (Figure 2E-G). The specialised role of STAG2 in gene expression (see answer to Q2 above) might give the STAG2 mutant an increased drug sensitivity profile that includes a wider range of pathways.
4) Western blotting shows cohesin mutant clones have increased total and membrane β-catenin that lead to transcriptional responses. Is there evidence of Wnt pathway/β-catenin transcriptional activation in publicly available datasets of cohesin perturbation or cohesin-mutant cancers?
This is an excellent suggestion, we have followed this up bioinformatically using publicly available TCGA data. We could only really look at STAG2-mutant cancers with high sample numbers in TCGA, because we needed a decent number of samples that had a cohesin mutation AND gene expression data. Here’s a summary of what we did:
STAG2 nonsense mutations (frameshift or stop-gained events) with gene expression were retrieved from TCGA by cancer type. Only cancer types having more than two STAG2 nonsense mutation samples were considered. After this filtering, four cancer types are remaining: Bladder Urothelial Carcinoma (BLCA), Uterine Corpus Endometrial Carcinoma (UCEC), Glioblastoma multiforme (GBM) and Cervical Kidney renal papillary cell carcinoma (KIRP).
For each cancer type, an equal number of control samples were also retrieved. Samples were considered as control if no STAG2 nonsense mutation was present.
When the 150 samples (all cancer type) were ranked by the Wald statistic, STAG2 was the first ranked gene. The comparison of this ranking to the gene set HALLMARK_WNT_Β_CATENIN ranked this pathway as 8th most enriched, overall showing an enrichment in up-regulated genes in the Wnt pathway. The results can be found in Supplementary file 2.
So, in the analysis we did, it looks like Wnt hyperactivation IS a feature of STAG2 mutant cancers. However, we treat the results quite cautiously, and only in the Discussion section, because of the caveats with this approach:
“Is Wnt pathway activation a conserved feature of cohesin-mutant cancers? Using publicly available data at TCGA [74], we correlated expression of genes in the Wnt pathway (Hallmark WNT β-catenin signaling) with nonsense mutation of STAG2, the most common of the cohesin mutations, in the four most represented cancers (bladder cancer, endometrial carcinoma, glioblastoma multiforme and cervical kidney renal papillary cell carcinoma). In comparative ranking, the Wnt pathway ranks 8th for enrichment in upregulated genes (after DNA repair/G2 checkpoint, MYC/E2F targets, oxidative phosphorylation and the unfolded protein response; Supplementary File 2). While our analysis could indicate conserved upregulation of the Wnt pathway in STAG2 mutant cancers, a caveat is that the true situation is likely to be more complex. If β-catenin is stabilized, upstream and in-parallel effectors such as Wnt ligands, receptors, signalling intermediates might trend to downregulation owing to negative feedback loop regulation of the Wnt pathway for example via DKK1 [75]. Furthermore, because cancer collections are unlikely to be under Wnt stimulation at the time of RNA collection, it is possible that their true Wnt sensitivity will be missed in a steady-state transcription analysis.”
5) The zebrafish experiments are fundamental and elegantly done. However, I think that, given the role of RAD21 and of the other cohesin subunits, it would be pivotal to include similar experiments with zebrafish modelling other genes of the complex.
We are happy to include timely new data that addresses this concern.
We recently generated zebrafish carrying mutations for the stag genes, stag1a, stag1b and stag2b. Characterisation of the Stag paralogues in zebrafish indicates that Stag2b is the most abundant Stag2, and likely to be most closely related to human STAG2. A description of the mutants and their early haematopoietic phenotypes has recently been accepted for publication (preprint here: https://www.biorxiv.org/content/10.1101/2020.10.19.346122v1).
The stag2b mutant is homozygous viable. We crossed this mutant to the TCF reporter line then crossed these fish (stag2b/+; TCF/+) to stag2b/+. Progeny carrying the TCF reporter were stimulated (or not) with LiCl, imaged, and genotyped post-imaging. Data for (stag2b/stag2b; TCF/+) and (+/+; TCF/+) genotypes were collected for a revision to Figure 6, now Figure 7 in the manuscript.
The results show higher baseline activity of the TCF Wnt reporter in stag2b homozygotes, but when stimulated with LiCl, we could not reliably determine extra expression above WT. Therefore, the stag2b mutant zebrafish also show enhanced Wnt signalling, but unlike for rad21-/- embryos, the enhancement is at baseline, pre-stimulation. Overall, the zebrafish experiments are consistent with over-active Wnt signalling being associated with cohesin mutation during embryogenesis. This information has been added to the Results:
“Zebrafish embryos heterozygous for Tg(7xTCF-Xla.Siam:nslmCherry)ia5, and either homozygous for stag2bnz207, rad21nz171 or wild type, were exposed to 0.15 M LiCl, an agonist of the Wnt signaling pathway (Figure 7). Expression of mCherry in the midbrain of embryos was detected at 20 hours post-fertilization (hpf) by epifluorescence and confocal imaging. A constitutive low level of mCherry is present in the developing midbrain of both untreated wild type (Figure 7A,B,I,J), and rad21nz171 embryos (Figure 7M,N). On the other hand, untreated stag2bnz207 embryos exhibited noticeably higher basal mCherry levels than wild type siblings (Figure 7E,F compared with A,B). This observation indicates that the Wnt pathway is more intrinsically active in these embryos. However, the addition of LiCl did not result in much additional fluorescence owing to stag2bnz207 mutation (Figure 7G,H compared with C,D). While mCherry expression in rad21nz171 mutants resembled that in wild type embryos at baseline, LiCl treatment dramatically increased the existing midbrain mCherry expression in rad21nz171 mutants compared with wild type embryos (Figure 7O,P compared with K,L). This observation indicates that rad21nz171 mutants are more sensitive to Wnt stimulation than wild type.”
Anastasia Labudina, who performed these experiments, has been added as an author on the manuscript.
6) Cohesin subunit deficient clones of MCF10A were generated using CRISPR-Cas9 editing. While this is clearly described in the Materials and methods section, it would be helpful to state the model system as well in results.
A brief description has been added to subsection “Cohesin gene deletion in MCF10A cells results in minor cell cycle defects”, and Supplementary file 1 has been redesignated as a Figure (Figure 1) describing generation of the mutant clones.
7) For cell cycle analysis (Results section) it would be extremely important, also considering the literature, to assess cell death rate.
We have included trypan blue measurement of cell survival in Figure 3—figure supplement 3E, described in the text: “Furthermore, none of the cohesin-deficient clones exhibited enhanced cell death compared with parental MCF10A cells (Figure 2B; Figure 3—figure supplement 3E).”
None of the cohesin-deficient clones have particularly enhanced cell death. Cell survival information is already included in the paper: the cell cycle flow cytometry analyses in Figure 2B do not show a sub-G1 population, which would be indicative of cell death.
8) In Figure 1F/1G/Figure 1—figure supplement 1B what do the square symbols represent? Please include in the legend.
This Figure is now Figure 2. The square symbols refer to the biological replicates for each experiment. There were actually 5 biological replicates per experiment, and this has now been corrected in the figure legend for both figures.
9) In Figure 4A/B/Figure 4—figure supplement 1A/1B please include size markers. Also Figure 4—figure supplement 1 needs a loading control such as tubulin.
This Figure is now Figure 5. For Figure 5A,B, we have provided a file showing the original blots with molecular size markers as Figure 5—source data 1. A loading control and molecular sizing has been added to Figure 5—figure supplement 1A/B. Original blots with molecular sizes have been provided for Figure 5—figure supplement 1 (Figure 5—source data 3) and Figure 5—figure supplement 2 (Figure 5—source data 4).
https://doi.org/10.7554/eLife.61405.sa2Article and author information
Author details
Funding
Health Research Council of New Zealand (15/229)
- Julia A Horsfield
Health Research Council of New Zealand (19/415)
- Ross D Hannan
- Julia A Horsfield
Associazione Italiana per la Ricerca sul Cancro (IG23284)
- Antonio Musio
The Maurice Wilkins centre for Molecular Biodiscovery (3705733)
- Jisha Antony
- Julia A Horsfield
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
The authors are grateful to the ANU Centre for Therapeutic Discovery at the John Curtin School of Medical Research at Australian National University for assistance with screening. We acknowledge Compounds Australia (https://www.compoundsaustralia.com) for their provision of specialized compound management and logistics research services to the project.
Ethics
Animal experimentation: Work with zebrafish was approved by the University of Otago (Dunedin) Animal Ethics Committee (AUP19/17) and conducted using approved institutional animal care standard operating procedures.
Senior Editor
- Richard M White, Memorial Sloan Kettering Cancer Center, United States
Reviewing Editor
- Ravi Majeti, Stanford University, United States
Reviewers
- Aaron Viny
- Massa Valentina, Università degli Studi di Milano, Italy
Version history
- Received: July 24, 2020
- Accepted: December 4, 2020
- Accepted Manuscript published: December 7, 2020 (version 1)
- Version of Record published: December 17, 2020 (version 2)
Copyright
© 2020, Chin et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 2,639
- Page views
-
- 368
- Downloads
-
- 15
- Citations
Article citation count generated by polling the highest count across the following sources: Crossref, PubMed Central, Scopus.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Cancer Biology
- Computational and Systems Biology
Pancreatic cancer is one of the deadliest cancer types with poor treatment options. Better detection of early symptoms and relevant disease correlations could improve pancreatic cancer prognosis. In this retrospective study, we used symptom and disease codes (ICD-10) from the Danish National Patient Registry (NPR) encompassing 6.9 million patients from 1994 to 2018,, of whom 23,592 were diagnosed with pancreatic cancer. The Danish cancer registry included 18,523 of these patients. To complement and compare the registry diagnosis codes with deeper clinical data, we used a text mining approach to extract symptoms from free text clinical notes in electronic health records (3078 pancreatic cancer patients and 30,780 controls). We used both data sources to generate and compare symptom disease trajectories to uncover temporal patterns of symptoms prior to pancreatic cancer diagnosis for the same patients. We show that the text mining of the clinical notes was able to complement the registry-based symptoms by capturing more symptoms prior to pancreatic cancer diagnosis. For example, ‘Blood pressure reading without diagnosis’, ‘Abnormalities of heartbeat’, and ‘Intestinal obstruction’ were not found for the registry-based analysis. Chaining symptoms together in trajectories identified two groups of patients with lower median survival (<90 days) following the trajectories ‘Cough→Jaundice→Intestinal obstruction’ and ‘Pain→Jaundice→Abnormal results of function studies’. These results provide a comprehensive comparison of the two types of pancreatic cancer symptom trajectories, which in combination can leverage the full potential of the health data and ultimately provide a fuller picture for detection of early risk factors for pancreatic cancer.
-
- Cancer Biology
Cancer stem cells (CSCs) undergo epithelial-mesenchymal transition (EMT) to drive metastatic dissemination in experimental cancer models. However, tumour cells undergoing EMT have not been observed disseminating into the tissue surrounding human tumour specimens, leaving the relevance to human cancer uncertain. We have previously identified both EpCAM and CD24 as CSC markers that, alongside the mesenchymal marker Vimentin, identify EMT CSCs in human oral cancer cell lines. This afforded the opportunity to investigate whether the combination of these three markers can identify disseminating EMT CSCs in actual human tumours. Examining disseminating tumour cells in over 12,000 imaging fields from 74 human oral tumours, we see a significant enrichment of EpCAM, CD24 and Vimentin co-stained cells disseminating beyond the tumour body in metastatic specimens. Through training an artificial neural network, these predict metastasis with high accuracy (cross-validated accuracy of 87-89%). In this study, we have observed single disseminating EMT CSCs in human oral cancer specimens, and these are highly predictive of metastatic disease.