Next Article in Journal
The Role of Cardiolipin in Mitochondrial Function and Neurodegenerative Diseases
Next Article in Special Issue
Exploring DIX-DIX Homo- and Hetero-Oligomers in Wnt Signaling with AlphaFold2
Previous Article in Journal
Cooperation of Various Cytoskeletal Components Orchestrates Intercellular Spread of Mitochondria between B-Lymphoma Cells through Tunnelling Nanotubes
Previous Article in Special Issue
When You Come to a Fork in the Road, Take It: Wnt Signaling Activates Multiple Pathways through the APC/Axin/GSK-3 Complex
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Synthetic Lethality between Cohesin and WNT Signaling Pathways in Diverse Cancer Contexts

by
Maria Michela Pallotta
,
Maddalena Di Nardo
and
Antonio Musio
*
Institute for Biomedical Technologies (ITB), National Research Council (CNR), 56124 Pisa, Italy
*
Author to whom correspondence should be addressed.
Cells 2024, 13(7), 608; https://doi.org/10.3390/cells13070608
Submission received: 27 February 2024 / Revised: 25 March 2024 / Accepted: 29 March 2024 / Published: 30 March 2024

Abstract

:
Cohesin is a highly conserved ring-shaped complex involved in topologically embracing chromatids, gene expression regulation, genome compartmentalization, and genome stability maintenance. Genomic analyses have detected mutations in the cohesin complex in a wide array of human tumors. These findings have led to increased interest in cohesin as a potential target in cancer therapy. Synthetic lethality has been suggested as an approach to exploit genetic differences in cancer cells to influence their selective killing. In this study, we show that mutations in ESCO1, NIPBL, PDS5B, RAD21, SMC1A, SMC3, STAG2, and WAPL genes are synthetically lethal with stimulation of WNT signaling obtained following LY2090314 treatment, a GSK3 inhibitor, in several cancer cell lines. Moreover, treatment led to the stabilization of β-catenin and affected the expression of c-MYC, probably due to the occupancy decrease in cohesin at the c-MYC promoter. Finally, LY2090314 caused gene expression dysregulation mainly involving pathways related to transcription regulation, cell proliferation, and chromatin remodeling. For the first time, our work provides the underlying molecular basis for synthetic lethality due to cohesin mutations and suggests that targeting the WNT may be a promising therapeutic approach for tumors carrying mutated cohesin.

1. Introduction

Cohesin is a member of the larger family of SMC (structural maintenance of chromosomes) complexes. The core of cohesin is composed of four subunits: SMC1A, SMC3, RAD21, and STAG1/2. Its activity is closely regulated by many accessory factors. Cohesin loading onto chromatin, mediated by the complex NIPBL–MAU2, occurs in G1 in yeast or at the end of the telophase of the previous cell cycle in mammalian cells. ESCO1/2 factors allow cohesin establishment in the S phase whereas PDS5 proteins ensure its maintenance. Once loaded onto the cohesin ring, the interaction of SMC3 head domain with DNA stimulates SMC3–SMC1A interface ATPase activity, responsible for opening of the cohesin ring. This allows WAPL to open the SMC3–RAD21 interface, totally releasing the DNA. Cohesin is phosphorylated at RAD21 and STAG subunits by polo-like kinase 1 during prophase and prometaphase; this process permits the removal of cohesin from chromosome arms. A modest amount of cohesin remains on chromosomes, especially at centromeres. Centromeric cohesin is certainly preserved by SGO1–PP2A complex, which stores cohesin in a hypophosphorylated state and maintains centromeric cohesion. At the metaphase–anaphase transition, the remaining cohesin is dissolved by the endopeptidase separase (ESPL1), which cleaves the cohesin’s RAD21 subunit. This cleavage permits opening of the cohesin ring, causing it to dissociate from chromosomes [1,2].
Ample experimental evidence highlights the role of cohesin as a pivotal player in various aspects of chromosome biology: correct chromosome segregation, DNA repair, genome stability, chromatin structure, 3D genome organization, and gene expression regulation. Considering its roles, it is not surprising that cancer genome and exome sequencing have revealed that cohesin subunits undergo a wide spectrum of somatic mutations in cancer. In fact, cohesin mutations have been identified in breast cancer [3,4], urothelial bladder carcinoma [5,6,7,8,9], lung carcinoma [10], Ewing’s sarcoma [11,12,13], colorectal carcinoma [14,15,16,17,18], glioblastoma [19,20], melanoma [21], and myeloid neoplasms [22,23,24,25]. According to the COSMIC (Catalogue of Somatic Mutations in Cancer, URL https://cancer.sanger.ac.uk/cosmic accessed on 01 August 2023), STAG1 (5%), NIPBL (4.9%), STAG2 (3.4%), and PDS5B (3.4%) are the most frequently mutated in cancer [26]. The multifaceted roles of cohesin may aid in hampering the growth of cohesin-mutant cancer cells by interfering with pathways that depend on correct cohesin function. Synthetic lethality indicates a genetic interaction where the concomitant dysfunction of two genes causes cell death, whereas with the alteration of a single gene, the cell remains viable. For instance, STAG2 shows a strong synthetic lethal interaction with its paralog STAG1 [27,28] while the inhibition of poly-ADP ribose polymerase (PARP) displays synthetic lethality with cohesin variants [29,30]. Recently, we showed that cohesin-mutant cancer cells were highly sensitive to treatment with LY2090314, a GSK3 inhibitor that acts as an agonist of the WNT signaling pathway. In particular, our investigation revealed that the stabilization of β-catenin in response to mutated cohesin probably plays a significant role in the heightened sensitivity of WNT target genes [31]. However, an outstanding question remains regarding the precise mechanisms underlying β-catenin stabilization and WNT activation in cohesin-mutant cells. Furthermore, to date, only a limited number of cohesin mutations have been tested for synthetic lethality. A larger number of genetic interactions can identify key processes that result in synthetic lethality when altered, providing a therapeutic strategy for cohesin-mutant cancer. In this study we first performed whole exome sequencing (WES) in cells belonging to the NCI-60 cancer cell line panel in order to identify mutations in cohesin genes. Via this approach, we detected fifteen mutations on both core cohesin and regulatory cohesin genes, namely ESCO1, ESPL1, NIPBL, PDS5B, RAD21, SMC1A, SMC3, STAG2, and WAPL genes. Next, cell lines harboring cohesin gene variants were treated with LY2090314. We found that treatment led to a strong decrease in cell viability and β-catenin stabilization in all cell lines except for that carrying ESPL1 mutation. Of note, the combination of LY2090314 treatment and cohesin deficiency caused a decrease in c-MYC (from hereon MYC) expression, likely to have been triggered by the occupancy decrease in both cohesin and β-catenin on the MYC promotor. Finally, RNA-seq analysis showed that LY2090314 treatment caused gene expression dysregulation affecting pathways related to chromatin remodeling, regulation of transcription, regulation of cell proliferation, and intracellular signal transduction. All together, these results indicate that all mutated cohesin genes, except for ESPL1, are synthetically lethal with stimulation of WNT signaling and this may offer a potential therapeutic target in tumors where cohesin is mutated.

2. Materials and Methods

2.1. Cell Culture and Treatment

The NCI-60 cancer cell line panel was purchased from the National Cancer Institute (USA). Cells were grown in RPMI supplemented with 10% FBS, 100 U/mL penicillin, 0.1 mg/mL streptomycin, and 1% L-glutamine in a humidified 5% CO2 atmosphere at 37 °C. Cells were treated with LY2090314 (Merck, Darmstadt, Germany) 100 nM for 24 h.

2.2. Whole Exome Sequencing (WES)

Genomic DNA was isolated from a diverse panel of NCI-60 cancer cell lines utilizing the PureDireX-Genomin DNA Isolation Kit (Bio-Helix, New Taipei City, Taiwan) employing a nonorganic extraction procedure. Subsequently, library preparation and exome enrichment were performed using the Nextera Flex for Enrichment solution (Illumina, San Diego, CA, USA) in conjunction with ‘SureSelect Human All Exon V7’ probes (Agilent Technologies, Santa Clara, CA, USA), targeting 35 Mb of human exonic content. Quantification and quality assessment of the samples were conducted utilizing the Qubit 2.0 Fluorometer (Thermo Fisher Scientific, Waltham, MA, USA) and Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA). The constructed libraries were subjected to sequencing on the NovaSeq 6000 platform (Illumina, San Diego, CA, USA) in 150 pair-end mode. Raw sequencing data underwent initial processing, including format conversion and de-multiplexing, using the Bcl2Fastq 2.0.21 version of the Illumina pipeline. Adapter sequences were subsequently masked from the raw fastq data utilizing Cutadapt v1.112 (accessed on 01 September 2023), employing parameters such as --anywhere (on both adapter sequences), --overlap 5, --times 2, --minimum-length 35, and --mask-adapter.
Further analysis was carried out using the Illumina DRAGEN Somatic Pipelines 3.5.7, (https://emea.illumina.com/products/by-type/informatics-products/basespacesequencehub/apps/edico-genome-inc-dragen-somatic-pipeline.html), which facilitated read mapping to the GRCh38/hg38 assembly and variant identification. Variants were functionally annotated using Annovar [32] (version Build #21596031) and SnpEff [33] (version 5.2). SNP and indel summary reports were generated, containing comprehensive information such as variant coordinates, base pair changes, amino acid change annotation, and functional annotations, including clinical significance, population frequencies, and various scores (SIFT, PolyPhen, LRT, MutationTaster), as well as HPO and other pertinent information essential for variant prioritization.

2.3. MTT Assay

HCT116, HCT15, A549, CCRF CEM, SF539, RFX393, UO31, IGROV1, HCC2998 cancer cell lines, along with MCF7 as the control cell line, were seeded onto 96-well plates at a concentration of 2 × 103 cells/well. Subsequently, MTT (3-(4,5-dimethyl-thiazol-2-yl)-2,5-diphenyltetrazolium bromide) solution at a concentration of 5 mg/mL was added to each well, and the plates were incubated at 37 °C for 4 h to allow the formation of formazan crystals. Following the incubation period, the formazan crystals were dissolved by adding acidified isopropanol (in 0.01 M HCl) to each well. The optical density of the resulting solution was measured using a microplate reader at a wavelength of 595 nm, allowing the quantification of cell viability.

2.4. Co-Immunoprecipitation (Co-IP)

For co-immunoprecipitation (co-IP) experiments, 800 μg of total protein extract from MCF7 cells was dissolved in 1 mL of incubation buffer. The solution underwent preclearing by incubation with 20μL of Dynabeads protein A (Invitrogen) for 1 h to reduce non-specific binding. Following preclearing, the supernatants were then incubated overnight at 4 °C with 7 μg of β-catenin antibody coupled to 40μL of Dynabeads protein A, facilitating the specific immunoprecipitation of target proteins. Subsequently, the immunocomplexes were boiled in sample buffer to dissociate the bound proteins and proteins were separated by SDS-PAGE. The Western blotting procedure involved probing the membrane with specific antibodies against NIPBL, RAD21, SMC1A, andSMC3.

2.5. Western Blotting

The harvested cell pellets were resuspended in lysis buffer containing tris HCl pH 8.0 (25 μM), NaCl (55 μM), EDTA (1 μM), and Protease Inhibitor Cocktail (Merck, Darmstadt, Germany). The protein concentration in the lysates was determined using the Bradford Protein Assay (Thermo Fisher Scientific, Waltham, MA, USA). Subsequently, proteins were loaded at a concentration of 5μg per lane and separated by SDS-PAGE. Following electrophoresis, proteins were transferred onto nitrocellulose membranes (Amersham, Marlborough, MA, USA) for subsequent immunoblotting. Non-specific binding was blocked with freshly prepared 5% BSA for 1 h on a shaking platform at room temperature. The membranes were then incubated with primary antibodies at a dilution of 1:1000, unless otherwise indicated, allowing specific detection of target proteins. After washing to remove unbound primary antibodies, the membranes were incubated with secondary antibody–peroxidase conjugate (Merck, Darmstadt, Germany) to facilitate detection. Chemiluminescence detection (Amersham, Marlborough, MA, USA) was employed to visualize the protein bands on the membranes, and images were captured using a Chemidoc imaging system (Bio-Rad, Hercules, CA, USA). To ensure equal loading of protein samples, anti-tubulin antibody was used as a loading control.

2.6. Antibodies

Primary antibodies were as follows: anti-γ-tubulin (T5326, 1:5000, Sigma-Aldrich), anti-β-catenin (#9562, CST), anti-GSK3 beta antibody [3D10] (ab93926, 1:1000, Abcam, Cambridge, UK), anti-SMC1A (A700-018, 1:1000, Fortis Life Sciences, Waltham, MA, USA), anti-NIPBL (A301-779A, 1:1000, Fortis Life Sciences, Waltham, MA, USA), anti-RAD21 (A300-080A, 1:1000, Fortis Life Sciences), and anti-SMC3 (A300-060A, 1:1000, Fortis Life Sciences, Waltham, MA, USA).

2.7. Chromatin Immunoprecipitation (ChIP)

To investigate the localization of β-catenin, NIPBL, RAD21, SMC1A, and SMC3 at the MYC locus, chromatin immunoprecipitation (ChIP) assays were conducted following previously established protocols with minor adjustments [34,35]. Initially, cells were crosslinked with 1% formaldehyde for 15 min, and the crosslinking reaction was quenched with the addition of 125 mM glycine. The crosslinked cells were then lysed successively with lysis buffer 1 (50 mM HEPES-KOH, pH 7.5, 140 mM NaCl, 1 mM EDTA, 10% glycerol, 0.5% NP-40, and 0.25% Triton X-100) and lysis buffer 2 (10 mM Tris-HCl, 100 mM NaCl, and 1 mM EDTA), and sonicated in lysis buffer 3 (10 mM Tris-HCl, 100 mM NaCl, 1 mM EDTA, 0.1% sodium deoxycholate, and 0.5% sarkosyl) to obtain chromatin fragments ranging from 500 to 200 base pairs. The sheared chromatin samples were then immunoprecipitated using Dynabeads protein A (Invitrogen) pre-bound with 10 μg of antibodies against β-catenin, NIPBL, RAD21, SMC1A, and SMC3. Mouse IgG (Sigma) served as a negative control. Following immunoprecipitation, the beads were washed successively with low-salt buffer (20 mM tris–HCl, pH 8, 150 mM NaCl, 0.5 mM EDTA, 0.1% SDS, 1% Triton X-100), high-salt buffer (20 mM tris–HCl, pH 8, 500 mM NaCl, 0.5 mM EDTA, 0.1% SDS, 1% Triton X-100), and RIPA buffer (50 mM HEPES-KOH, pH 7.5, 0.5 mM EDTA, 10% NP-40, 0.7% sodium deoxycolate and 0.5 M LiCl) to remove non-specifically bound proteins. The immunoprecipitated chromatin was eluted from the beads overnight at 65 °C and subsequently treated with proteinase K. DNA was then extracted using phenol-chloroform extraction, followed by purification with the QIAquick Purification Kit (Qiagen, Hilden, Germany). For each cell line, three independent ChIP assays were performed and analyzed via quantitative PCR (qPCR). Each sample was run in duplicate and repeated three times to ensure reproducibility. Specific primers targeting either the promoter or the first exon of MYC gene were utilized for qPCR analysis (Supplementary Table S1).
The results were expressed using the percent input method, where signals obtained from the ChIP were normalized to signals obtained from an input sample, representing the amount of chromatin used in the ChIP.

2.8. RNA-Sequencing (RNA-seq)

For RNA-seq analysis, two samples from each cell line, both before and after LY2090314 treatment, underwent individual processing as previously described [35,36]. Library preparation was performed using the Universal Plus mRNA-Seq kit (Tecan Genomics, Redwood City, CA, USA) according to the manufacturer’s instructions (library type: fr-secondstrand). RNA samples were quantified and subjected to quality assessment using either the Agilent 2100 Bioanalyzer RNA assay (Agilent Technologies, Santa Clara, CA, USA) or the Caliper LabChip GX (PerkinElmer, Waltham, MA, USA). Final libraries were evaluated for quality using both the Qubit 2.0 Fluorometer (Thermo Fisher Scientific, Waltham, MA, USA) and either the Agilent Bioanalyzer DNA assay or the Caliper LabChip GX (PerkinElmer, Waltham, MA, USA). Subsequently, the libraries were prepared for sequencing and sequenced in single-end 75 bp mode on the NextSeq 500 (Illumina, San Diego, CA, USA)/paired-end 150 bp mode on the NovaSeq 6000 (Illumina, San Diego, CA, USA). Raw sequencing data underwent initial processing using the Illumina BCL Convert v3.9.31 pipeline for format conversion and de-multiplexing. The reads were then aligned to the reference hg38-iGenomes genome using STAR, a splice junction mapper for RNA-seq reads [37]. This alignment method facilitated the identification of splice junctions between exons. Preprocessing of RNA-seq data for differential expression analysis was performed using Htseq-count to count the overlap of reads with genes. DESeq2 [38,39] was subsequently employed to compare expression levels of genes and transcripts by fitting a generalized linear model (GLM) for each gene. Normalization was conducted using the median-of-ratios method [38], and statistical significance was determined using a Wald test.
Only protein-coding genes were considered for analysis, and gene-level expression values were determined using fragments per kilobase million (FPKM) mapped. Genes with FPKM > 1 were considered expressed and analyzed, with a false discovery rate (FDR) threshold of <0.001 for statistical significance.

2.9. Pathway Analysis and Function

To gain insights into the biological significance of the differentially expressed genes, functional analysis was conducted using the Database for Annotation, Visualization, and Integrated Discovery (DAVID) v2023q2 [40] (URL “https://david.ncifcrf.gov” accessed on 09 October 2023). This online tool enables the annotation and visualization of gene lists, providing valuable information on enriched biological processes. For each term identified by DAVID, a p-value was calculated, and terms with a significance level of p < 0.05 were considered enriched.

2.10. cDNA Synthesis and Quantitative Real-Time PCR (qPCR)

Total RNA extraction was carried out using the RNAeasy Mini-kit (Qiagen), following the manufacturer’s protocol. Subsequently, cDNA synthesis was performed using SuperScript™ II reverse transcriptase in conjunction with oligo-dT primers (Invitrogen). PCR analyses were conducted using the Rotor Gene 3000 (Corbett, Sydney, Australia). qPCR reactions were run in triplicate to ensure reproducibility and reliability of the results. To normalize the qPCR data, expression levels were referenced to the housekeeping gene HPRT. Specific primer pairs utilized for the validation of RNA-seq data and MYC expression are listed in Supplementary Table S2.

2.11. Statistical Analyses

Statistical analyses were conducted using the IBM SPSS statistical package (version 28.0) for Windows. Results were reported as mean ± standard deviation (SD). Differences between continuous variables were assessed using Student’s t-test. A p-value < 0.05 was considered statistically significant.

3. Results

3.1. Identification of Mutated Cohesin Genes in Cancer Cell Lines

The list of cell lines in the NCI-60 panel and their tissue origins are given in Supplementary Table S3. The panel consists of fifty-five cell lines: two from the prostate, four from blood, six each from the colorectal system, breast, ovary, and central nervous system, eight each from the kidney and the lung and nine from melanoma. The identification of mutations in cohesin genes was performed via WES. Although a limitation of cell line sequencing is the lack of available normal-matched tissue for comparison, the NCI-60 panel does allow comparisons between cell lines from nine distinct tissues of origin. We detected fifteen mutations in nine cell lines each carrying at least one cohesin mutation. In particular, one was derived from each of lung, ovary, blood, and central nervous system, two from kidney, and three from colon cancer (Table 1). Identified variants were mapped on both core cohesin and regulatory cohesin genes, involving ESCO1, ESPL1, NIPBL, PDS5B, RAD21, SMC1A, SMC3, STAG2, and WAPL genes. Most mutations (80%, 12 out of 15) were nonsense or frameshift leading to premature stop codons. The remaining three mutations were missense, the first two (c.3288G > T, c.5038A > G) affecting ESPL1 and the last (c.2028G > T) involving SMC1A (Table 1). Overall, WES analysis showed that 16% (9 out of 55) of cell lines belonging to the NCI-60 panel carry at least one mutation in cohesin or regulatory cohesin genes.

3.2. Effect of LY2090314 Treatment on Cancer Cell Lines Carrying Mutations in Cohesin Genes

Since the availability of cancer cell lines harboring mutations in cohesin genes is very rare, this finding allowed us to investigate whether all identified cohesin mutations were synthetically lethal with stimulation of WNT signaling and whether targeting the WNT pathway may prove to be a novel therapeutic strategy for cohesin-mutant cancers. To this aim, all cohesin-mutated cancer cells were treated with LY2090314 (100 nM for 24 h), a GSK3 inhibitor and stimulator of the WNT signaling pathway. The activation of WNT signaling via ligand binding or GSK3 inhibition releases β-catenin, thus allowing it to accumulate in the nucleus and activate transcription of WNT target gene/s [41,42].
After treatment, we analyzed GSK3 and β-catenin protein levels via Western blotting. We found that there was no difference in GSK3 levels between cohesin-mutated and control cells after treatment of cells with LY2090314 (Figure 1A and Figure S1A). This was somewhat expected since, upon administration, LY2090314 binds to and inhibits GSK-3 in an ATP-competitive manner but treatment has no effect on GSK3 protein level, as previously shown [43]. However, this result contrasts with what has been previously described [31]. It is possible that this discrepancy arises from the combination of several factors. In fact, in our work, the treatment with LY2090314 was applied for 24 h and the total proteins were subsequently extracted, while in the previous work, the treatment lasted 6 h and then the proteins were extracted from the cytoplasmic or membrane fraction [31].
Instead, LY2090314 treatment markedly increased β-catenin in all mutated cells compared with controls (Figure 1B and Figure S1B). Next, we investigated the cell viability after WNT stimulation by using MTT assay. Interestingly, we found that LY2090314 inhibited the growth of all cohesin-mutated cell lines apart from the MCF7 control cell line and IGROV1 carrying an ESPL1 mutation. It is interesting to note that cells carrying multiple cohesin mutations exhibited reduced viability compared with those with a single mutation. For instance, the A549 cell line, which harbors mutations in both WAPL and NIPBL genes, showed diminished viability compared with the CCRFCEM cell line, which only has a mutation in NIPBL, highlighting the contribution of WAPL to synthetic lethality. In addition, the MTT assay allowed discrimination of the role of the SMC1A mutation. Indeed, the IGROV1 cell line, which carries the mutation in ESPL1, showed no differences in viability compared with the control line. In contrast, the HCT15 cell line, which harbors mutations in both ESPL1 and SMC1A, exhibited growth delay after GSK3 inhibitor treatment (Figure 1C). In order to investigate a possible interaction between cohesin and β-catenin, we performed co-IP between cohesin members and β-catenin. We found that NIPBL, RAD21, SMC1A, and SMC3 were detected in IP–β-catenin precipitates, while no signal was detected in the control Western blotting using IgG-coated beads (Figure 1D and Figure S2).
These results suggest that β-catenin is stabilized following WNT activation when cohesin activity is perturbed. Furthermore, for the first time, we have demonstrated that β-catenin interacts with cohesin factors.

3.3. The Interplay between β-Catenin, Cohesin and MYC

It is well known that the accumulation of β-catenin promotes the transcription of many oncogenes such as MYC [44]. In addition, previous studies have shown that MYC transcription is positively regulated by cohesin in Drosophila, zebrafish, mouse, and human cells [45,46,47,48]. Therefore, we analyzed the expression of MYC in cohesin-mutated cancer cells after LY2090314 treatment.
We found that treatment caused a strong decrease in MYC expression level in all cell lines, except for IGROV1 and MCF7 cells, which on the contrary displayed a significant increase. In particular, SF539, HCC2998, RFX393, CCRFCEM, and UO31 showed a decrease of < 80%, A549 and HCT15 almost 50%, while IGROV1 and MCF7 displayed increases of 45% and 68%, respectively (Figure 2A).
WNT signalling was previously shown to induce enhancer-promoter looping at the MYC gene in cancer cells [49]. Therefore, we investigated both cohesin and β-catenin binding at MYC promoters via ChIP-qPCR analysis. In order to investigate the effect of gene mutation on the binding of β-catenin, we focused on cell lines harboring a single cohesin mutation. The MYC locus contains three exons and four distinct promoters. P0 transcripts begin at many initiation sites, while both P1 and P2 are placed at the beginning of the first exon. Finally, P3 maps between the first and second exon. Exon 1, a large noncoding exon, contains an enhancer element and a region that oversees the elongation of RNA transcripts [50,51]. We found that β-catenin did not bind to P0 promoter while its recruitment at the first exon showed a significant decrease in cell lines with a single cohesin mutation. In contrast, β-catenin binding displayed a significant increase in MCF7 and IGROV1 cell lines (Figure 2B). Next, we also investigated the recruitment of core cohesin members and regulatory factors at MYC promoters. ChIP data showed that binding of NIPBL, RAD21, SMC1A, and SMC3 decreased significantly at both P0 and P1/P2 promoters in CCRFCEM, RFX393, UO31, and HCT116 cell lines, respectively, whereas ESPL1 showed a significant increase in the IGROV1 cell line (p < 0.05) (Figure 2C,D and Figure S3).
All together, these results suggest that LY2090314 treatment leads to decreased MYC expression, and this probably depends on the reduced binding of both cohesin and β-catenin at the first exon of MYC.

3.4. LY2090314 Treatment Causes Gene Expression Dysregulation

In order to investigate the effect of LY2090314 treatment on gene expression against a cohesin-mutated background, we performed RNA-seq analysis. Treatment caused the dysregulated expression of many genes ranging from 6162 in RFX393 to 3253 in HCT116 cell lines (FDR < 0.001, Table 2).
In the RFX393, SF539, and UO31 cell lines carrying mutations in the core cohesin genes, RAD21, SMC1A, and STAG2 showed the highest number of dysregulated genes. The complete list of dysregulated genes is reported in DataSet. Cell lines shared 133 dysregulated genes (Supplementary Table S4). The classification of dysregulated genes according to molecular function and biological process via the DAVID tool (URL “https://david.ncifcrf.gov” accessed on 09 October 2023) showed that identified pathways are related to positive regulation of transcription (GO:0045893), positive regulation of cell proliferation (GO:0008284), chromatin remodeling (GO:0006338), and intracellular signal transduction (GO:0035556). The first ten pathways for a number of genes are shown in Figure 3, while the remaining pathways are listed in Supplementary Table S5.
To find a comprehensive relationship between the 133 dysregulated proteins, the STRING network analysis tool (URL “https://string-db.org” accessed on 10 October 2023) was applied. The protein interactions were filtered, so that unconnected root nodes and indirect interactions were discarded. This analysis predicts that the dysregulated proteins are related to two clusters, centered on KDM5B, ATP5A1, and YAP1 and ERBB4, respectively, whose members are involved in proteosome and cancer progression (Figure 3. The heat map related to the core of the two clusters (see Figure 4) showed that 97% of genes maintained the same trend, though to a different extent (Supplementary Figure S4, Supplementary Table S6). RNA-seq data were validated in five genes, namely PSMD7, PSMD4, PSMB9, SFHM1, and PAX6, via quantitative RT-PCR experiments (Supplementary Figure S5). These genes were chosen because most of them are involved in cancer development [52,53,54,55] and their differential expression could explain the decreased viability after LY2090314 treatment.
Collectively, following treatment with LY2090314, the cell lines shared 133 dysregulated genes involved in fundamental biological processes whose changes in gene expression could underlie synthetic lethality in cells with mutated cohesin genes.

4. Discussion

Advances in cancer genomics have allowed the identification of potential therapeutic targets in relation to genomic defects that were formerly disregarded due to their non-driver mutation classification [56,57]. The cohesin complex is frequently mutated in cancer and it is thought that cohesin dysregulation plays a significant role in cancer development and progression [26]. These findings have attracted increased interest supporting the idea of cohesin as a possible therapeutic target [17,58]. Recently, we showed that cohesin subunits SMC3, RAD21, and STAG2 are more responsive to treatment with LY2090314, a GSK3 inhibitor and an agonist of WNT signaling [31]. This phenomenon is called synthetic lethality, where the occurrence of a single genetic event is tolerable for cell survival, whereas the co-occurrence of a second event results in cell death. To date, only a limited number of cohesin gene mutations have been tested for synthetic lethality and the mechanism underlying synthetic lethality is still under debate.
Since the availability of cell lines carrying cohesin mutations is limited, we first performed WES for an NCI-60 cancer panel. We identified at least one cohesin mutation in nine cell lines derived from colon, kidney, lung, ovary, blood, and central nervous system cancers. Core cohesin subunits and regulatory factors were both represented. In fact, we detected missense, frameshift, or nonsense mutations in ESCO1, ESPL1, NIPBL, PDS5B, RAD21, SMC1A, SMC3, STAG2, and WAPL genes.
Next, these cell lines were treated with LY2090314 compound to test its effect on cell viability in cells carrying cohesin mutation deriving from different cancer contexts. We found no difference in GSK3 levels between MCF7 cells, a control cell line without cohesin mutations, and mutant cohesin cells. In contrast, we found that β-catenin was stabilized in all mutant cohesin cells following LY2090314 treatment. It is worth noting that HCT116 and HCT15 cell lines are constitutively active for WNT signaling. In fact, HCT116 cells carry a three-base deletion in the β-catenin encoding gene causing an in-frame deletion at Ser45 while HCT15 harbors a truncated mutation in the APC gene [59,60,61,62]. The migration of β-catenin into the nucleus represents a pivotal event in the canonical WNT pathway [63,64]. This process hinges on the destabilization of the β-catenin destruction complex, a crucial assembly comprising CK1, GSK3, Axin, and APC, culminating in the potentiation of WNT signaling [65,66,67,68]. These findings indicate that mutations in APC or β-catenin thwart the destruction complex’s ability to bind to β-catenin, leading to aberrant accumulation of β-catenin within the cell. In addition, our results suggest that β-catenin normally targeted for degradation is stabilized in cohesin-deficient cells and is available in the active form following inhibition of GSK3, which in turn enhances WNT signaling. The WNT pathway and β-catenin are involved in establishing the body axis during embryonic growth [69,70,71]. In adults, they contribute to tissue stem cell maintenance, homeostasis, apoptosis, and cell renewal, proliferation, and repair [72,73]. The aberrant expression of WNT signaling/β-catenin has been widely considered to play a pivotal role in tumorigenesis, cancer progression, and metastasis (reviewed in [74]). Furthermore, we found that all mutated cells showed growth delay after treatment, with the only exception being a cell line harboring ESPL1 mutation. This result could reflect the different roles of cohesin proteins under investigation. RAD21, SMC1A, SMC3 belong to the core cohesin complex, while NIPBL and WAPL have an antagonist role. NIPBL is involved in the loading of cohesin onto chromosomes [75]. Once loaded, cohesin can translocate along the chromosome or be removed by the cohesin-releasing factor WAPL [76,77,78]. Instead, ESPL1, although essential for sister chromatid cohesion dissolution, is not a member of the cohesin complex [79]. Therefore, unlike other proteins, ESPL1 intervenes only in the final phase of the cohesin life cycle and, thus, its mutations have no effect on sensitivity to treatment with LY2090314.
It is well known that the accumulation of β-catenin promotes the transcription of many oncogenes, such as MYC [44,80]. MYC is a transcription factor that promotes the rapid down- and upregulation of about 15% of human genes and, as consequence, it is thought to coordinate the plethora of transcriptional changes that foster cell proliferation [81,82].
Here, we found that activation of the WNT pathway caused a considerable decrease in both cell viability and MYC expression level in all cell lines, apart from that carrying ESPL1 mutation, as previously stated. This data indicates that LY2090314 exerts its effects through the downregulation of MYC, due to β-catenin stabilization and enhancement of WNT signaling, consequently affecting cell growth (Figure 5). This study suggests that WNT enhancement may be effective for the treatment of tumors containing cohesin mutations.
Based on the results presented in this study, we show that SMC1A, RAD21, and SMC3 bind all promoter regions of MYC. This result indicates that MYC is positively regulated by cohesin, further supporting this notion as previously observed in Drosophila, zebrafish, and mouse cells [45,48]. The recruitment of cohesin on MYC promoter regions showed a significant decrease following LY2090314 treatment. Interestingly, we found that β-catenin recruitment also decreases at the MYC locus and β-catenin interacts with cohesin complex. This suggests that the regulation of MYC expression is finely regulated, involving both cohesin and β-catenin. As an effect of MYC dysregulation, we found the alteration of the expression of hundreds of genes. Cell lines carrying mutations in RAD21, SMC1A, and STAG2 displayed the highest numbers of dysregulated genes. This observation further confirms the role of cohesin in regulating gene expression, suggesting that altered cohesin activity leads to aberrant transcription [83]. We hypothesize that this points to changes in gene expression as a cancer-promoting mechanism. This notion is supported by the finding that RAD21, SMC1A, and STAG2 have been thought to play a pivotal role in tumorigenesis [84,85,86]. All cell lines carrying cohesin mutations shared 133 genes. These genes were involved in several biological processes, such as regulation of transcription, intracellular signal transduction, negative regulation of apoptotic process, positive regulation of cell proliferation, and chromatin remodeling. In addition, STRING analysis showed that dysregulated proteins belong to two clusters, centered on KDM5B, ATP5A1 and YAP1 and ERBB4. which are involved in proteosomes and tumorigenesis [87,88,89,90,91,92]. The altered expression of these proteins could underlie the synthetic lethal phenomenon exerted by LY2090314 treatment in a genetic environment with mutated cohesin.

5. Conclusions

In summary, we identified the molecular mechanism underlying synthetic lethality due to cohesin mutations. Moreover, synthetic lethality offers a promising approach for developing valuable therapeutic interventions in cancer when direct targeting of driver genes is impractical. Thus, our findings indicate that targeting the WNT pathway may be effective in treating cancers that harbor cohesin mutations.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/cells13070608/s1, Table S1: Primers used to investigate cohesin loading at the MYC locus. Table S2: Primers used to validate RNA-seq data and MYC expression by RT-PCR. Table S3: NCI-60 cancer panel cell used in this work. Table S4: Dysregulated genes in common among cohesion-mutated cancer cells. Table S5: Classification of dysregulated genes by molecular function and biological process through the DAVID tool. Table S6: Fold change of the core genes as identified through STRING analysis. Figure S1: Effects of LY2090314 (LY) treatment in cancer cell lines carrying mutation in cohesin genes. Figure S2: β-catenin and cohesin interaction. Figure S3. Effects of LY2090314 (LY) treatment on the MYC locus. Figure S4: Heatmap of dysregulated genes which were held in common among cancer cell lines carrying mutation in cohesin genes. Figure S5: RNA-seq data was validated by RT-PCR.

Author Contributions

M.M.P. and M.D.N. performed WES, in vitro cell biology experiments, and RNA-seq analysis. M.M.P. and A.M. conceived the study, designed experiments, and interpreted the data. A.M. wrote the manuscript with input from M.M.P. and M.D.N. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by a grant from Fondazione Associazione Italiana per la Ricerca sul Cancro (IG 23284) to A.M. and by a seed grant to M.M.P.

Data Availability Statement

The relevant data supporting the findings of this study are available in this article and its Supplementary Information files. All NGS raw files have been deposited in the NCBI Sequence Read Archive under accession number PRJNA1019503.

Conflicts of Interest

The authors declare no conflicts of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Hoencamp, C.; Rowland, B.D. Genome control by SMC complexes. Nat. Rev. Mol. Cell Biol. 2023, 24, 633–650. [Google Scholar] [CrossRef] [PubMed]
  2. Yatskevich, S.; Rhodes, J.; Nasmyth, K. Organization of Chromosomal DNA by SMC Complexes. Annu. Rev. Genet. 2019, 53, 445–482. [Google Scholar] [CrossRef] [PubMed]
  3. Koedoot, E.; van Steijn, E.; Vermeer, M.; Gonzalez-Prieto, R.; Vertegaal, A.C.O.; Martens, J.W.M.; Le Devedec, S.E.; van de Water, B. Splicing factors control triple-negative breast cancer cell mitosis through SUN2 interaction and sororin intron retention. J. Exp. Clin. Cancer Res. 2021, 40, 82. [Google Scholar] [CrossRef] [PubMed]
  4. Oishi, Y.; Nagasaki, K.; Miyata, S.; Matsuura, M.; Nishimura, S.I.; Akiyama, F.; Iwai, T.; Miki, Y. Functional pathway characterized by gene expression analysis of supraclavicular lymph node metastasis-positive breast cancer. J. Hum. Genet. 2007, 52, 271–279. [Google Scholar] [CrossRef] [PubMed]
  5. Balbas-Martinez, C.; Sagrera, A.; Carrillo-de-Santa-Pau, E.; Earl, J.; Marquez, M.; Vazquez, M.; Lapi, E.; Castro-Giner, F.; Beltran, S.; Bayes, M.; et al. Recurrent inactivation of STAG2 in bladder cancer is not associated with aneuploidy. Nat. Genet. 2013, 45, 1464–1469. [Google Scholar] [CrossRef] [PubMed]
  6. Solomon, D.A.; Kim, J.S.; Bondaruk, J.; Shariat, S.F.; Wang, Z.F.; Elkahloun, A.G.; Ozawa, T.; Gerard, J.; Zhuang, D.; Zhang, S.; et al. Frequent truncating mutations of STAG2 in bladder cancer. Nat. Genet. 2013, 45, 1428–1430. [Google Scholar] [CrossRef] [PubMed]
  7. Guo, G.; Sun, X.; Chen, C.; Wu, S.; Huang, P.; Li, Z.; Dean, M.; Huang, Y.; Jia, W.; Zhou, Q.; et al. Whole-genome and whole-exome sequencing of bladder cancer identifies frequent alterations in genes involved in sister chromatid cohesion and segregation. Nat. Genet. 2013, 45, 1459–1463. [Google Scholar] [CrossRef] [PubMed]
  8. Taylor, C.F.; Platt, F.M.; Hurst, C.D.; Thygesen, H.H.; Knowles, M.A. Frequent inactivating mutations of STAG2 in bladder cancer are associated with low tumour grade and stage and inversely related to chromosomal copy number changes. Hum. Mol. Genet. 2014, 23, 1964–1974. [Google Scholar] [CrossRef]
  9. Cancer Genome Atlas Research, N. Comprehensive molecular characterization of urothelial bladder carcinoma. Nature 2014, 507, 315–322. [Google Scholar] [CrossRef]
  10. Zhu, H.E.; Li, T.; Shi, S.; Chen, D.X.; Chen, W.; Chen, H. ESCO2 promotes lung adenocarcinoma progression by regulating hnRNPA1 acetylation. J. Exp. Clin. Cancer Res. 2021, 40, 64. [Google Scholar] [CrossRef]
  11. Crompton, B.D.; Stewart, C.; Taylor-Weiner, A.; Alexe, G.; Kurek, K.C.; Calicchio, M.L.; Kiezun, A.; Carter, S.L.; Shukla, S.A.; Mehta, S.S.; et al. The genomic landscape of pediatric Ewing sarcoma. Cancer Discov. 2014, 4, 1326–1341. [Google Scholar] [CrossRef]
  12. Brohl, A.S.; Solomon, D.A.; Chang, W.; Wang, J.; Song, Y.; Sindiri, S.; Patidar, R.; Hurd, L.; Chen, L.; Shern, J.F.; et al. The genomic landscape of the Ewing Sarcoma family of tumors reveals recurrent STAG2 mutation. PLoS Genet. 2014, 10, e1004475. [Google Scholar] [CrossRef] [PubMed]
  13. Tirode, F.; Surdez, D.; Ma, X.; Parker, M.; Le Deley, M.C.; Bahrami, A.; Zhang, Z.; Lapouble, E.; Grossetete-Lalami, S.; Rusch, M.; et al. Genomic landscape of Ewing sarcoma defines an aggressive subtype with co-association of STAG2 and TP53 mutations. Cancer Discov. 2014, 4, 1342–1353. [Google Scholar] [CrossRef] [PubMed]
  14. Di Nardo, M.; Astigiano, S.; Baldari, S.; Pallotta, M.M.; Porta, G.; Pigozzi, S.; Antonini, A.; Emionite, L.; Frattini, A.; Valli, R.; et al. The synergism of SMC1A cohesin gene silencing and bevacizumab against colorectal cancer. J. Exp. Clin. Cancer Res. 2024, 43, 49. [Google Scholar] [CrossRef] [PubMed]
  15. Barber, T.D.; McManus, K.; Yuen, K.W.; Reis, M.; Parmigiani, G.; Shen, D.; Barrett, I.; Nouhi, Y.; Spencer, F.; Markowitz, S.; et al. Chromatid cohesion defects may underlie chromosome instability in human colorectal cancers. Proc. Natl. Acad. Sci. USA 2008, 105, 3443–3448. [Google Scholar] [CrossRef] [PubMed]
  16. Cucco, F.; Servadio, A.; Gatti, V.; Bianchi, P.; Mannini, L.; Prodosmo, A.; De Vitis, E.; Basso, G.; Friuli, A.; Laghi, L.; et al. Mutant cohesin drives chromosomal instability in early colorectal adenomas. Hum. Mol. Genet. 2014, 23, 6. [Google Scholar] [CrossRef] [PubMed]
  17. Sarogni, P.; Palumbo, O.; Servadio, A.; Astigiano, S.; D’Alessio, B.; Gatti, V.; Cukrov, D.; Baldari, S.; Pallotta, M.M.; Aretini, P.; et al. Overexpression of the cohesin-core subunit SMC1A contributes to colorectal cancer development. J. Exp. Clin. Cancer Res. 2019, 38, 108. [Google Scholar] [CrossRef] [PubMed]
  18. Wang, J.; Yu, S.; Cui, L.; Wang, W.; Li, J.; Wang, K.; Lao, X. Role of SMC1A overexpression as a predictor of poor prognosis in late stage colorectal cancer. BMC Cancer 2015, 15, 90. [Google Scholar] [CrossRef] [PubMed]
  19. Brennan, C.W.; Verhaak, R.G.; McKenna, A.; Campos, B.; Noushmehr, H.; Salama, S.R.; Zheng, S.; Chakravarty, D.; Sanborn, J.Z.; Berman, S.H.; et al. The somatic genomic landscape of glioblastoma. Cell 2013, 155, 462–477. [Google Scholar] [CrossRef]
  20. Bailey, M.L.; O’Neil, N.J.; van Pel, D.M.; Solomon, D.A.; Waldman, T.; Hieter, P. Glioblastoma cells containing mutations in the cohesin component STAG2 are sensitive to PARP inhibition. Mol. Cancer Ther. 2014, 13, 724–732. [Google Scholar] [CrossRef]
  21. Ryu, B.; Kim, D.S.; Deluca, A.M.; Alani, R.M. Comprehensive expression profiling of tumor cell lines identifies molecular signatures of melanoma progression. PLoS ONE 2007, 2, e594. [Google Scholar] [CrossRef] [PubMed]
  22. Kon, A.; Shih, L.Y.; Minamino, M.; Sanada, M.; Shiraishi, Y.; Nagata, Y.; Yoshida, K.; Okuno, Y.; Bando, M.; Nakato, R.; et al. Recurrent mutations in multiple components of the cohesin complex in myeloid neoplasms. Nat. Genet. 2013, 45, 1232–1237. [Google Scholar] [CrossRef] [PubMed]
  23. Thota, S.; Viny, A.D.; Makishima, H.; Spitzer, B.; Radivoyevitch, T.; Przychodzen, B.; Sekeres, M.A.; Levine, R.L.; Maciejewski, J.P. Genetic alterations of the cohesin complex genes in myeloid malignancies. Blood 2014, 124, 1790–1798. [Google Scholar] [CrossRef]
  24. Thol, F.; Bollin, R.; Gehlhaar, M.; Walter, C.; Dugas, M.; Suchanek, K.J.; Kirchner, A.; Huang, L.; Chaturvedi, A.; Wichmann, M.; et al. Mutations in the cohesin complex in acute myeloid leukemia: Clinical and prognostic implications. Blood 2014, 123, 914–920. [Google Scholar] [CrossRef] [PubMed]
  25. Papaemmanuil, E.; Gerstung, M.; Bullinger, L.; Gaidzik, V.I.; Paschka, P.; Roberts, N.D.; Potter, N.E.; Heuser, M.; Thol, F.; Bolli, N.; et al. Genomic Classification and Prognosis in Acute Myeloid Leukemia. N. Engl. J. Med. 2016, 374, 2209–2221. [Google Scholar] [CrossRef] [PubMed]
  26. Di Nardo, M.; Pallotta, M.M.; Musio, A. The multifaceted roles of cohesin in cancer. J. Exp. Clin. Cancer Res. 2022, 41, 96. [Google Scholar] [CrossRef] [PubMed]
  27. van der Lelij, P.; Lieb, S.; Jude, J.; Wutz, G.; Santos, C.P.; Falkenberg, K.; Schlattl, A.; Ban, J.; Schwentner, R.; Hoffmann, T.; et al. Synthetic lethality between the cohesin subunits STAG1 and STAG2 in diverse cancer contexts. eLife 2017, 6, e26980. [Google Scholar] [CrossRef] [PubMed]
  28. Benedetti, L.; Cereda, M.; Monteverde, L.; Desai, N.; Ciccarelli, F.D. Synthetic lethal interaction between the tumour suppressor STAG2 and its paralog STAG1. Oncotarget 2017, 8, 37619–37632. [Google Scholar] [CrossRef] [PubMed]
  29. Liu, Y.; Xu, H.; Van der Jeught, K.; Li, Y.; Liu, S.; Zhang, L.; Fang, Y.; Zhang, X.; Radovich, M.; Schneider, B.P.; et al. Somatic mutation of the cohesin complex subunit confers therapeutic vulnerabilities in cancer. J. Clin. Investig. 2018, 128, 2951–2965. [Google Scholar] [CrossRef]
  30. Mondal, G.; Stevers, M.; Goode, B.; Ashworth, A.; Solomon, D.A. A requirement for STAG2 in replication fork progression creates a targetable synthetic lethality in cohesin-mutant cancers. Nat. Commun. 2019, 10, 1686. [Google Scholar] [CrossRef]
  31. Chin, C.V.; Antony, J.; Ketharnathan, S.; Labudina, A.; Gimenez, G.; Parsons, K.M.; He, J.; George, A.J.; Pallotta, M.M.; Musio, A.; et al. Cohesin mutations are synthetic lethal with stimulation of WNT signaling. eLife 2020, 9, e61405. [Google Scholar] [CrossRef]
  32. Wang, K.; Li, M.; Hakonarson, H. ANNOVAR: Functional annotation of genetic variants from high-throughput sequencing data. Nucleic Acids Res. 2010, 38, e164. [Google Scholar] [CrossRef]
  33. Cingolani, P.; Platts, A.; Wang, L.L.; Coon, M.; Nguyen, T.; Wang, L.; Land, S.J.; Lu, X.; Ruden, D.M. A program for annotating and predicting the effects of single nucleotide polymorphisms, SnpEff: SNPs in the genome of Drosophila melanogaster strain w1118; iso-2; iso-3. Fly 2012, 6, 80–92. [Google Scholar] [CrossRef]
  34. Cucco, F.; Palumbo, E.; Camerini, S.; D’Alessio, B.; Quarantotti, V.; Casella, M.L.; Rizzo, I.M.; Cukrov, D.; Delia, D.; Russo, A.; et al. Separase prevents genomic instability by controlling replication fork speed. Nucleic Acids Res. 2018, 46, 267–278. [Google Scholar] [CrossRef]
  35. Mannini, L.; Lamaze, F.C.; Cucco, F.; Amato, C.; Quarantotti, V.; Rizzo, I.M.; Krantz, I.D.; Bilodeau, S.; Musio, A. Mutant cohesin affects RNA polymerase II regulation in Cornelia de Lange syndrome. Sci. Rep. 2015, 5, 16803. [Google Scholar] [CrossRef] [PubMed]
  36. Mannini, L.; Cucco, F.; Quarantotti, V.; Amato, C.; Tinti, M.; Tana, L.; Frattini, A.; Delia, D.; Krantz, I.D.; Jessberger, R.; et al. SMC1B is present in mammalian somatic cells and interacts with mitotic cohesin proteins. Sci. Rep. 2015, 5, 18472. [Google Scholar] [CrossRef]
  37. Dobin, A.; Davis, C.A.; Schlesinger, F.; Drenkow, J.; Zaleski, C.; Jha, S.; Batut, P.; Chaisson, M.; Gingeras, T.R. STAR: Ultrafast universal RNA-seq aligner. Bioinformatics 2013, 29, 15–21. [Google Scholar] [CrossRef] [PubMed]
  38. Anders, S.; Huber, W. Differential expression analysis for sequence count data. Genome Biol. 2010, 11, R106. [Google Scholar] [CrossRef] [PubMed]
  39. Love, M.I.; Huber, W.; Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014, 15, 550. [Google Scholar] [CrossRef]
  40. Sherman, B.T.; Hao, M.; Qiu, J.; Jiao, X.; Baseler, M.W.; Lane, H.C.; Imamichi, T.; Chang, W. DAVID: A web server for functional enrichment analysis and functional annotation of gene lists (2021 update). Nucleic Acids Res. 2022, 50, W216–W221. [Google Scholar] [CrossRef]
  41. Logan, C.Y.; Nusse, R. The Wnt signaling pathway in development and disease. Annu. Rev. Cell Dev. Biol. 2004, 20, 781–810. [Google Scholar] [CrossRef] [PubMed]
  42. Clevers, H.; Loh, K.M.; Nusse, R. Stem cell signaling. An integral program for tissue renewal and regeneration: Wnt signaling and stem cell control. Science 2014, 346, 1248012. [Google Scholar] [CrossRef] [PubMed]
  43. Kunnimalaiyaan, S.; Schwartz, V.K.; Jackson, I.A.; Clark Gamblin, T.; Kunnimalaiyaan, M. Antiproliferative and apoptotic effect of LY2090314, a GSK-3 inhibitor, in neuroblastoma in vitro. BMC Cancer 2018, 18, 560. [Google Scholar] [CrossRef] [PubMed]
  44. Higgs, M.R.; Lerat, H.; Pawlotsky, J.M. Hepatitis C virus-induced activation of beta-catenin promotes c-Myc expression and a cascade of pro-carcinogenetic events. Oncogene 2013, 32, 4683–4693. [Google Scholar] [CrossRef] [PubMed]
  45. Rhodes, J.M.; Bentley, F.K.; Print, C.G.; Dorsett, D.; Misulovin, Z.; Dickinson, E.J.; Crosier, K.E.; Crosier, P.S.; Horsfield, J.A. Positive regulation of c-Myc by cohesin is direct, and evolutionarily conserved. Dev. Biol. 2010, 344, 637–649. [Google Scholar] [CrossRef]
  46. Schaaf, C.A.; Misulovin, Z.; Sahota, G.; Siddiqui, A.M.; Schwartz, Y.B.; Kahn, T.G.; Pirrotta, V.; Gause, M.; Dorsett, D. Regulation of the Drosophila Enhancer of split and invected-engrailed gene complexes by sister chromatid cohesion proteins. PLoS ONE 2009, 4, e6202. [Google Scholar] [CrossRef] [PubMed]
  47. Kawauchi, S.; Calof, A.L.; Santos, R.; Lopez-Burks, M.E.; Young, C.M.; Hoang, M.P.; Chua, A.; Lao, T.; Lechner, M.S.; Daniel, J.A.; et al. Multiple organ system defects and transcriptional dysregulation in the Nipbl(+/−) mouse, a model of Cornelia de Lange Syndrome. PLoS Genet. 2009, 5, e1000650. [Google Scholar] [CrossRef]
  48. Pallotta, M.M.; Di Nardo, M.; Sarogni, P.; Krantz, I.D.; Musio, A. Disease-associated c-MYC downregulation in human disorders of transcriptional regulation. Hum. Mol. Genet. 2022, 31, 1599–1609. [Google Scholar] [CrossRef] [PubMed]
  49. Liu, J.; Xiao, Q.; Xiao, J.; Niu, C.; Li, Y.; Zhang, X.; Zhou, Z.; Shu, G.; Yin, G. Wnt/beta-catenin signalling: Function, biological mechanisms, and therapeutic opportunities. Signal Transduct. Target. Ther. 2022, 7, 3. [Google Scholar] [CrossRef]
  50. Bentley, D.L.; Groudine, M. Novel promoter upstream of the human c-myc gene and regulation of c-myc expression in B-cell lymphomas. Mol. Cell Biol. 1986, 6, 3481–3489. [Google Scholar] [CrossRef]
  51. Yang, J.Q.; Remmers, E.F.; Marcu, K.B. The first exon of the c-myc proto-oncogene contains a novel positive control element. EMBO J. 1986, 5, 3553–3562. [Google Scholar] [CrossRef] [PubMed]
  52. Wang, J.; Wang, T.; Feng, Y.K.; Liu, Y.; Fu, B.; Liu, X.T.; Wu, Q.Z. Deubiquitinating enzyme PSMD7 promotes bladder cancer development: Involvement of RAB1A stabilization. Cell Signal 2024, 114, 110996. [Google Scholar] [CrossRef] [PubMed]
  53. Teng, X.; Yang, T.; Yuan, B.; Yang, Y.; Liu, J.; Wang, X.; Wang, Y.; Ma, T.; Yin, X.; Yu, H.; et al. Prognostic analysis of patients with breast cancer based on tumor mutational burden and DNA damage repair genes. Front. Oncol. 2023, 13, 1177133. [Google Scholar] [CrossRef] [PubMed]
  54. Tamilzhalagan, S.; Muthuswami, M.; Periasamy, J.; Lee, M.H.; Rha, S.Y.; Tan, P.; Ganesan, K. Upregulated, 7q21-22 amplicon candidate gene SHFM1 confers oncogenic advantage by suppressing p53 function in gastric cancer. Cell Signal 2015, 27, 1075–1086. [Google Scholar] [CrossRef] [PubMed]
  55. Kiselev, Y.; Andersen, S.; Johannessen, C.; Fjukstad, B.; Standahl Olsen, K.; Stenvold, H.; Al-Saad, S.; Donnem, T.; Richardsen, E.; Bremnes, R.M.; et al. Transcription factor PAX6 as a novel prognostic factor and putative tumour suppressor in non-small cell lung cancer. Sci. Rep. 2018, 8, 5059. [Google Scholar] [CrossRef] [PubMed]
  56. Liu, Y.; Hu, X.; Han, C.; Wang, L.; Zhang, X.; He, X.; Lu, X. Targeting tumor suppressor genes for cancer therapy. Bioessays 2015, 37, 1277–1286. [Google Scholar] [CrossRef] [PubMed]
  57. Muller, F.L.; Aquilanti, E.A.; DePinho, R.A. Collateral Lethality: A new therapeutic strategy in oncology. Trends Cancer 2015, 1, 161–173. [Google Scholar] [CrossRef] [PubMed]
  58. Antony, J.; Chin, C.V.; Horsfield, J.A. Cohesin Mutations in Cancer: Emerging Therapeutic Targets. Int. J. Mol. Sci. 2021, 22, 6788. [Google Scholar] [CrossRef]
  59. Ilyas, M.; Tomlinson, I.P.; Rowan, A.; Pignatelli, M.; Bodmer, W.F. Beta-catenin mutations in cell lines established from human colorectal cancers. Proc. Natl. Acad. Sci. USA 1997, 94, 10330–10334. [Google Scholar] [CrossRef]
  60. Zhang, Q.; Yang, X.; Wu, J.; Ye, S.; Gong, J.; Cheng, W.M.; Luo, Z.; Yu, J.; Liu, Y.; Zeng, W.; et al. Reprogramming of palmitic acid induced by dephosphorylation of ACOX1 promotes beta-catenin palmitoylation to drive colorectal cancer progression. Cell Discov. 2023, 9, 26. [Google Scholar] [CrossRef]
  61. Kawasaki, Y.; Sato, R.; Akiyama, T. Mutated APC and Asef are involved in the migration of colorectal tumour cells. Nat. Cell Biol. 2003, 5, 211–215. [Google Scholar] [CrossRef] [PubMed]
  62. Rowan, A.J.; Lamlum, H.; Ilyas, M.; Wheeler, J.; Straub, J.; Papadopoulou, A.; Bicknell, D.; Bodmer, W.F.; Tomlinson, I.P. APC mutations in sporadic colorectal tumors: A mutational “hotspot” and interdependence of the “two hits”. Proc. Natl. Acad. Sci. USA 2000, 97, 3352–3357. [Google Scholar] [CrossRef]
  63. Basu, S.; Haase, G.; Ben-Ze’ev, A. Wnt signaling in cancer stem cells and colon cancer metastasis. F1000Res 2016, 5, 699. [Google Scholar] [CrossRef] [PubMed]
  64. Zhan, T.; Rindtorff, N.; Boutros, M. Wnt signaling in cancer. Oncogene 2017, 36, 1461–1473. [Google Scholar] [CrossRef] [PubMed]
  65. Amit, S.; Hatzubai, A.; Birman, Y.; Andersen, J.S.; Ben-Shushan, E.; Mann, M.; Ben-Neriah, Y.; Alkalay, I. Axin-mediated CKI phosphorylation of beta-catenin at Ser 45: A molecular switch for the Wnt pathway. Genes. Dev. 2002, 16, 1066–1076. [Google Scholar] [CrossRef] [PubMed]
  66. Behrens, J.; Jerchow, B.A.; Wurtele, M.; Grimm, J.; Asbrand, C.; Wirtz, R.; Kuhl, M.; Wedlich, D.; Birchmeier, W. Functional interaction of an axin homolog, conductin, with beta-catenin, APC, and GSK3beta. Science 1998, 280, 596–599. [Google Scholar] [CrossRef] [PubMed]
  67. Kishida, S.; Yamamoto, H.; Ikeda, S.; Kishida, M.; Sakamoto, I.; Koyama, S.; Kikuchi, A. Axin, a negative regulator of the wnt signaling pathway, directly interacts with adenomatous polyposis coli and regulates the stabilization of beta-catenin. J. Biol. Chem. 1998, 273, 10823–10826. [Google Scholar] [CrossRef] [PubMed]
  68. Lee, E.; Salic, A.; Kruger, R.; Heinrich, R.; Kirschner, M.W. The roles of APC and Axin derived from experimental and theoretical analysis of the Wnt pathway. PLoS Biol. 2003, 1, E10. [Google Scholar] [CrossRef] [PubMed]
  69. Clark, E.G.; Petersen, C.P. BMP suppresses WNT to integrate patterning of orthogonal body axes in adult planarians. PLoS Genet. 2023, 19, e1010608. [Google Scholar] [CrossRef]
  70. Bennett, J.T.; Stickney, H.L.; Choi, W.Y.; Ciruna, B.; Talbot, W.S.; Schier, A.F. Maternal nodal and zebrafish embryogenesis. Nature 2007, 450, E1–E2. [Google Scholar] [CrossRef]
  71. Arnold, C.P.; Benham-Pyle, B.W.; Lange, J.J.; Wood, C.J.; Sanchez Alvarado, A. Wnt and TGFbeta coordinate growth and patterning to regulate size-dependent behaviour. Nature 2019, 572, 655–659. [Google Scholar] [CrossRef] [PubMed]
  72. Steinhart, Z.; Angers, S. Wnt signaling in development and tissue homeostasis. Development 2018, 145, dev146589. [Google Scholar] [CrossRef]
  73. Bastakoty, D.; Young, P.P. Wnt/beta-catenin pathway in tissue injury: Roles in pathology and therapeutic opportunities for regeneration. FASEB J. 2016, 30, 3271–3284. [Google Scholar] [CrossRef]
  74. Bugter, J.M.; Fenderico, N.; Maurice, M.M. Mutations and mechanisms of WNT pathway tumour suppressors in cancer. Nat. Rev. Cancer 2021, 21, 5–21. [Google Scholar] [CrossRef] [PubMed]
  75. Ciosk, R.; Shirayama, M.; Shevchenko, A.; Tanaka, T.; Toth, A.; Nasmyth, K. Cohesin’s binding to chromosomes depends on a separate complex consisting of Scc2 and Scc4 proteins. Mol. Cell 2000, 5, 243–254. [Google Scholar] [CrossRef] [PubMed]
  76. Lengronne, A.; Katou, Y.; Mori, S.; Yokobayashi, S.; Kelly, G.P.; Itoh, T.; Watanabe, Y.; Shirahige, K.; Uhlmann, F. Cohesin relocation from sites of chromosomal loading to places of convergent transcription. Nature 2004, 430, 573–578. [Google Scholar] [CrossRef] [PubMed]
  77. Gandhi, R.; Gillespie, P.J.; Hirano, T. Human Wapl is a cohesin-binding protein that promotes sister-chromatid resolution in mitotic prophase. Curr. Biol. 2006, 16, 2406–2417. [Google Scholar] [CrossRef] [PubMed]
  78. Kueng, S.; Hegemann, B.; Peters, B.H.; Lipp, J.J.; Schleiffer, A.; Mechtler, K.; Peters, J.M. Wapl controls the dynamic association of cohesin with chromatin. Cell 2006, 127, 955–967. [Google Scholar] [CrossRef] [PubMed]
  79. Hauf, S.; Waizenegger, I.C.; Peters, J.M. Cohesin cleavage by separase required for anaphase and cytokinesis in human cells. Science 2001, 293, 1320–1323. [Google Scholar] [CrossRef]
  80. Shang, S.; Hua, F.; Hu, Z.W. The regulation of beta-catenin activity and function in cancer: Therapeutic opportunities. Oncotarget 2017, 8, 33972–33989. [Google Scholar] [CrossRef]
  81. Zeller, K.I.; Zhao, X.; Lee, C.W.; Chiu, K.P.; Yao, F.; Yustein, J.T.; Ooi, H.S.; Orlov, Y.L.; Shahab, A.; Yong, H.C.; et al. Global mapping of c-Myc binding sites and target gene networks in human B cells. Proc. Natl. Acad. Sci. USA 2006, 103, 17834–17839. [Google Scholar] [CrossRef]
  82. Fernandez, P.C.; Frank, S.R.; Wang, L.; Schroeder, M.; Liu, S.; Greene, J.; Cocito, A.; Amati, B. Genomic targets of the human c-Myc protein. Genes. Dev. 2003, 17, 1115–1129. [Google Scholar] [CrossRef]
  83. Horsfield, J.A. Full circle: A brief history of cohesin and the regulation of gene expression. FEBS J. 2023, 290, 1670–1687. [Google Scholar] [CrossRef]
  84. Musio, A. The multiple facets of the SMC1A gene. Gene 2020, 743, 144612. [Google Scholar] [CrossRef]
  85. Cheng, H.; Zhang, N.; Pati, D. Cohesin subunit RAD21: From biology to disease. Gene 2020, 758, 144966. [Google Scholar] [CrossRef] [PubMed]
  86. Romero-Perez, L.; Surdez, D.; Brunet, E.; Delattre, O.; Grunewald, T.G.P. STAG Mutations in Cancer. Trends Cancer 2019, 5, 506–520. [Google Scholar] [CrossRef] [PubMed]
  87. Metzler, V.M.; de Brot, S.; Haigh, D.B.; Woodcock, C.L.; Lothion-Roy, J.; Harris, A.E.; Nilsson, E.M.; Ntekim, A.; Persson, J.L.; Robinson, B.D.; et al. The KDM5B and KDM1A lysine demethylases cooperate in regulating androgen receptor expression and signalling in prostate cancer. Front. Cell Dev. Biol. 2023, 11, 1116424. [Google Scholar] [CrossRef] [PubMed]
  88. He, R.; Xhabija, B.; Gopi, L.K.; Kurup, J.T.; Xu, Z.; Liu, Z.; Kidder, B.L. H3K4 demethylase KDM5B regulates cancer cell identity and epigenetic plasticity. Oncogene 2022, 41, 2958–2972. [Google Scholar] [CrossRef] [PubMed]
  89. Mohamed, T.I.A.; Ezugwu, A.E.; Fonou-Dombeu, J.V.; Mohammed, M.; Greeff, J.; Elbashir, M.K. A novel feature selection algorithm for identifying hub genes in lung cancer. Sci. Rep. 2023, 13, 21671. [Google Scholar] [CrossRef]
  90. Thilakasiri, P.; O’Keefe, R.N.; To, S.Q.; Chisanga, D.; Eissmann, M.F.; Carli, A.L.; Duscio, B.; Baloyan, D.; Dmello, R.S.; Williams, D.; et al. Mechanisms of cellular crosstalk in the gastric tumor microenvironment are mediated by YAP1 and STAT3. Life Sci. Alliance 2024, 7, e202302411. [Google Scholar] [CrossRef]
  91. Hu, R.; Hou, H.; Li, Y.; Zhang, M.; Li, X.; Chen, Y.; Guo, Y.; Sun, H.; Zhao, S.; Liao, M.; et al. Combined BET and MEK Inhibition synergistically suppresses melanoma by targeting YAP1. Theranostics 2024, 14, 593–607. [Google Scholar] [CrossRef] [PubMed]
  92. Wills, C.; He, Y.; Summers, M.G.; Lin, Y.; Phipps, A.I.; Watts, K.; Law, P.J.; Al-Tassan, N.A.; Maughan, T.S.; Kaplan, R.; et al. A genome-wide search for determinants of survival in 1926 patients with advanced colorectal cancer with follow-up in over 22,000 patients. Eur. J. Cancer 2021, 159, 247–258. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Effects of LY2090314 (LY) treatment in cancer cell lines carrying mutations in cohesin genes. (A) Western blotting showing that GSK3 levels did not change in MCF7 (1, untreated; 2, plus LY); HCT116 (3, untreated; 4, plus LY2090314); HCT15 (5, untreated; 6, plus LY), and IGROV1 (7, untreated; 8, plus LY4) cell lines. (B) LY2090314 treatment caused an increase in the expression of β-catenin in MCF7 (1, untreated; 2, plus LY); HCT116 (3, untreated; 4, plus LY); HCT15 (5, untreated; 6, plus LY), and IGROV1 (7, untreated; 8, plus LY) cell lines. (C) MTT assay shows that LY2090314 treatment caused the reduction of viability in all cell lines except for MCF7, the control cell lines without mutation in cohesin genes, and IGROV1 carrying mutation in the ESPL1 gene. (D) β-catenin was found to be co-precipitated with SMC1A, whereas no SMC1A signal was detected in the IPs using IgG-coated beads. * p < 0.05.
Figure 1. Effects of LY2090314 (LY) treatment in cancer cell lines carrying mutations in cohesin genes. (A) Western blotting showing that GSK3 levels did not change in MCF7 (1, untreated; 2, plus LY); HCT116 (3, untreated; 4, plus LY2090314); HCT15 (5, untreated; 6, plus LY), and IGROV1 (7, untreated; 8, plus LY4) cell lines. (B) LY2090314 treatment caused an increase in the expression of β-catenin in MCF7 (1, untreated; 2, plus LY); HCT116 (3, untreated; 4, plus LY); HCT15 (5, untreated; 6, plus LY), and IGROV1 (7, untreated; 8, plus LY) cell lines. (C) MTT assay shows that LY2090314 treatment caused the reduction of viability in all cell lines except for MCF7, the control cell lines without mutation in cohesin genes, and IGROV1 carrying mutation in the ESPL1 gene. (D) β-catenin was found to be co-precipitated with SMC1A, whereas no SMC1A signal was detected in the IPs using IgG-coated beads. * p < 0.05.
Cells 13 00608 g001
Figure 2. Effects of LY2090314 (LY) treatment on MYC locus. (A) MYC was downregulated in all cell lines carrying mutation in cohesin genes, apart from IGROV1 harboring ESPL1 mutation and MCF7, the control cell line. (B) LY treatment caused decreased binding of β-catenin at P1/P2 promoters of MYC locus in all cell lines carrying a single cohesin mutation, again except for IGROV1 and MCF7 cell lines. (C,D) ChIP experiments with antibodies against SMC1A or SMC3 showed decreased cohesin binding at both P0 and P1/2 promoters of the MYC gene following LY treatment in UO31 and HCT116, respectively. ** p < 0.01, * p < 0.05.
Figure 2. Effects of LY2090314 (LY) treatment on MYC locus. (A) MYC was downregulated in all cell lines carrying mutation in cohesin genes, apart from IGROV1 harboring ESPL1 mutation and MCF7, the control cell line. (B) LY treatment caused decreased binding of β-catenin at P1/P2 promoters of MYC locus in all cell lines carrying a single cohesin mutation, again except for IGROV1 and MCF7 cell lines. (C,D) ChIP experiments with antibodies against SMC1A or SMC3 showed decreased cohesin binding at both P0 and P1/2 promoters of the MYC gene following LY treatment in UO31 and HCT116, respectively. ** p < 0.01, * p < 0.05.
Cells 13 00608 g002
Figure 3. Effects of LY2090314 treatment on gene expression. GO term enrichment analysis of biological process that were significantly over-represented when considering differentially expressed genes following LY2090314 treatment.
Figure 3. Effects of LY2090314 treatment on gene expression. GO term enrichment analysis of biological process that were significantly over-represented when considering differentially expressed genes following LY2090314 treatment.
Cells 13 00608 g003
Figure 4. Effects of LY2090314 treatment on gene expression. STRING network analysis of the 133 dysregulated proteins. They were grouped in two different clusters.
Figure 4. Effects of LY2090314 treatment on gene expression. STRING network analysis of the 133 dysregulated proteins. They were grouped in two different clusters.
Cells 13 00608 g004
Figure 5. Model for the synthetic lethal interaction between LY2090314 treatment and mutated cohesin. In cancer cells with and without cohesin gene mutations, LY2090314 inhibits GSK3, which acts as an agonist of the WNT signaling pathway. In cells without cohesin mutation, this leads to an increase in both cohesin and β-catenin binding at MYC promoters, stimulating cell growth. Instead, in the presence of mutated cohesin, the recruitment of both β-catenin and cohesin at the MYC locus is decreased and MYC expression is downregulated, leading to decreased viability.
Figure 5. Model for the synthetic lethal interaction between LY2090314 treatment and mutated cohesin. In cancer cells with and without cohesin gene mutations, LY2090314 inhibits GSK3, which acts as an agonist of the WNT signaling pathway. In cells without cohesin mutation, this leads to an increase in both cohesin and β-catenin binding at MYC promoters, stimulating cell growth. Instead, in the presence of mutated cohesin, the recruitment of both β-catenin and cohesin at the MYC locus is decreased and MYC expression is downregulated, leading to decreased viability.
Cells 13 00608 g005
Table 1. Cohesin gene mutations identified in the NCI-60 cancer cell line panel.
Table 1. Cohesin gene mutations identified in the NCI-60 cancer cell line panel.
Cell LineTissueGeneDNA MutationProtein Mutation
HCT116Colon SMC3c.2425C > Tp.Gln809 *
HCT15Colon ESPL1c.3288G > Tp.Glu1096Asp
SMC1Ac.2028G > Tp.Glu676Asp
A549LungNIPBLc.7201G > Tp.Gly2401 *
WAPLc.916C > Tp.Gln306 *
CCRFCEMBloodNIPBLc.8212C > Tp.Gln2738 *
SF539Central nervous systemSTAG2c.775C > Tp.Arg259 *
RFX393KidneyRAD21L1c.31C > Tp.Arg11 *
UO31KidneySMC1Ac.855_856insp.Glu286fs
IGROV1OvaryESPL1c.5038A > Gp.Arg1680Gly
HCC2998Colon SMC3c.2830C > Tp.Arg944 *
ESCO1c.1369G > Tp.Glu457 *
ESCO1c.925G > Tp.Glu309 *
STAG2c.1147G > Tp.Glu383 *
PDS5Bc.2740G > Tp.Glu914 *
* premature stop codon.
Table 2. Numbers of dysregulated genes in cancer cells carrying mutations in cohesin genes.
Table 2. Numbers of dysregulated genes in cancer cells carrying mutations in cohesin genes.
Cell LineDysregulated GenesUpregulated GenesDownregulated Genes
A549480822402568
CCRFCEM540331242279
HCC2998424221772065
HCT15369320001693
HCT116325316041649
IGROV1274713331414
RFX393616228263336
SF539585529212934
UO31593626313305
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Pallotta, M.M.; Di Nardo, M.; Musio, A. Synthetic Lethality between Cohesin and WNT Signaling Pathways in Diverse Cancer Contexts. Cells 2024, 13, 608. https://doi.org/10.3390/cells13070608

AMA Style

Pallotta MM, Di Nardo M, Musio A. Synthetic Lethality between Cohesin and WNT Signaling Pathways in Diverse Cancer Contexts. Cells. 2024; 13(7):608. https://doi.org/10.3390/cells13070608

Chicago/Turabian Style

Pallotta, Maria Michela, Maddalena Di Nardo, and Antonio Musio. 2024. "Synthetic Lethality between Cohesin and WNT Signaling Pathways in Diverse Cancer Contexts" Cells 13, no. 7: 608. https://doi.org/10.3390/cells13070608

APA Style

Pallotta, M. M., Di Nardo, M., & Musio, A. (2024). Synthetic Lethality between Cohesin and WNT Signaling Pathways in Diverse Cancer Contexts. Cells, 13(7), 608. https://doi.org/10.3390/cells13070608

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop