# BMAL1 Knockdown Leans Epithelial–Mesenchymal Balance toward Epithelial Properties and Decreases the Chemoresistance of Colon Carcinoma Cells

Yuan Zhang, Aurore Devocelle, Christophe Desterke, Lucas Eduardo Botelho de Souza, Éva Hadadi, Hervé Acloque, Adlen Foudi, Yao Xiang, Annabelle Ballesta, Yunhua Chang, et al.

<figure><img src="image_1.png" /></figure>

► To cite this version:

Yuan Zhang, Aurore Devocelle, Christophe Desterke, Lucas Eduardo Botelho de Souza, Éva Hadadi, et al.. BMAL1 Knockdown Leans Epithelial–Mesenchymal Balance toward Epithelial Properties and Decreases the Chemoresistance of Colon Carcinoma Cells. International Journal of Molecular Sciences, 2021, 22, 10.3390/ijms22105247. [https://doi.org/10.3390/ijms22105247](https://doi.org/10.3390/ijms22105247)

HAL Id: hal-04427569

https://hal.inrae.fr/hal-04427569

Submitted on 30 Jan 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>
<figure><img src="image_1.png" /></figure>

HAL
open science

# BMAL1 Knockdown Leans Epithelial–Mesenchymal Balance toward Epithelial Properties and Decreases the Chemoresistance of Colon Carcinoma Cells

Yuan Zhang, Aurore Devocelle, Christophe Desterke, Lucas Eduardo Botelho de Souza, Éva Hadadi, Hervé Acloque, Adlen Foudi, Yao Xiang, Annabelle Ballesta, Yunhua Chang, et al.

<figure><img src="image_2.png" /></figure>

## ► To cite this version:

Yuan Zhang, Aurore Devocelle, Christophe Desterke, Lucas Eduardo Botelho de Souza, Éva Hadadi, et al.. BMAL1 Knockdown Leans Epithelial–Mesenchymal Balance toward Epithelial Properties and Decreases the Chemoresistance of Colon Carcinoma Cells. International Journal of Molecular Sciences, 2021, 22, 10.3390/ijms22105247. [https://doi.org/10.3390/ijms22105247](https://doi.org/10.3390/ijms22105247)

HAL Id: hal-04427569

https://hal.inrae.fr/hal-04427569v1

Submitted on 30 Jan 2024

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_3.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

Article

BMAL1 Knockdown Leans Epithelial–Mesenchymal Balance
toward Epithelial Properties and Decreases the
Chemoresistance of Colon Carcinoma Cells

Yuan Zhang 1,2,3,†, Aurore Devocelle 2,4,†, Christophe Desterke 1,2, Lucas Eduardo Botelho de Souza 1,2,
Éva Hadadi 1,2, Hervé Acloque 1,2, Adlen Foudi 1,2 ID, Yao Xiang 5, Annabelle Ballesta 6, Yunhua Chang 1,2,5,‡ ID
and Julien Giron-Michel 2,4,*,‡ ID

<figure><img src="image_4.png" /></figure>

check for
updates

**Citation:** Zhang, Y.; Devocelle, A.; Desterke, C.; de Souza, L.E.B.; Hadadi, É.; Acloque, H.; Foudi, A.; Xiang, Y.; Ballesta, A.; Chang, Y.; et al. BMAL1 Knockdown Leans Epithelial–Mesenchymal Balance toward Epithelial Properties and Decreases the Chemoresistance of Colon Carcinoma Cells. *Int. J. Mol. Sci.* **2021**, *22*, 5247. https://doi.org/10.3390/ijms22105247

Academic Editor: Monica Fedele

Received: 13 April 2021

Accepted: 6 May 2021

Published: 16 May 2021

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_5.png" /></figure>

**Copyright:** © 2021 by the authors.
Licensee MDPI, Basel, Switzerland.
This article is an open access article
distributed under the terms and
conditions of the Creative Commons
Attribution (CC BY) license (https://
creativecommons.org/licenses/by/
4.0/).

**Abstract:** The circadian clock coordinates biological and physiological functions to day/night cycles. The perturbation of the circadian clock increases cancer risk and affects cancer progression. Here, we studied how BMAL1 knockdown (BMAL1-KD) by shRNA affects the epithelial–mesenchymal transition (EMT), a critical early event in the invasion and metastasis of colorectal carcinoma (CRC). In corresponding to a gene set enrichment analysis, which showed a significant enrichment of EMT and invasive signatures in BMAL1_high CRC patients as compared to BMAL1_low CRC patients, our results revealed that BMAL1 is implicated in keeping the epithelial–mesenchymal equilibrium of CRC cells and influences their capacity of adhesion, migration, invasion, and chemoresistance. Firstly, BMAL1-KD increased the expression of epithelial markers (E-cadherin, CK-20, and EpCAM) but decreased the expression of Twist and mesenchymal markers (N-cadherin and vimentin) in CRC cell lines. Finally, the molecular alterations after BMAL1-KD promoted mesenchymal-to-epithelial transition-like changes mostly appeared in two primary CRC cell lines (i.e., HCT116 and SW480) compared to the metastatic cell line SW620. As a consequence, migration/invasion and drug resistance capacities decreased in HCT116 and SW480 BMAL1-KD cells. Together, BMAL1-KD alerts the delicate equilibrium between epithelial and mesenchymal properties of CRC cell lines, which revealed the crucial role of BMAL1 in EMT-related CRC metastasis and chemoresistance.

**Keywords:** colorectal cancer; epithelial–mesenchymal transition (EMT); circadian clock; BMAL1; metastasis; chemoresistance

1. Introduction

Colorectal cancer (CRC) is one of the most common cancers and represents the fourth
leading cause of cancer death worldwide. In 2019, there were 145,600 new cases and
approximately 51,020 deaths in the United States [1]. Metastasis in liver and lung represents
one of the main causes of CRC-related mortality [2]. An insight into the molecular events
driving metastasis, particularly the enhanced invasiveness and therapeutic resistance, is
crucial for developing some novel treatment regimens to combat CRC. One such critical

<footnote n="1">1. INSERM UMR-S 935, CNRS Campus, 94801 Villejuif, France; zyljxmm88@gmail.com (Y.Z.); christophe.desterke@inserm.fr (C.D.); lucasebsouza@usp.br (L.E.B.d.S.); Eva.Hadadi@vub.be (É.H.); herve.acloque@inrae.fr (H.A.); adlen.foudi@inserm.fr (A.F.); yunhua.chang-marchand@inserm.fr (Y.C.)</footnote>

<footnote n="2">2 Orsay-Vallée Campus, Paris-Saclay University, 91190 Gif-sur-Yvette, France; aurore.devocelle@inserm.fr</footnote>

<footnote n="3">3 Institute of Life Sciences, Chongqing Medical University, Chongqing 400016, China</footnote>

<footnote n="4">4 INSERM UMR-S-MD 1197/Ministry of the Armed Forces, Biomedical Research Institute of the Armed Forces (IRBA), Paul-Brousse Hospital Villejuif and CTSA Clamart, 94807 Villejuif, France</footnote>

<footnote n="5">5 INSERM UMR-S 1151, Department of Immunology, Infectiology and Hematology, Institut Necker-Enfants Malades (INEM), Paris Descartes University, CNRS UMR 8253, 75730 Paris, France; yao.xiang@inserm.fr</footnote>

<footnote n="6">6 INSERM UMR-S 900, Institut Curie, MINES ParisTech CBIO, PSL Research University, 92210 Saint-Cloud, France; annabelle.ballesta@curie.fr</footnote>

<footnote>* Correspondence: julien.giron-michel@inserm.fr; Tel.: +33-1-45-59-53-27</footnote>

<footnote>† The two authors contribute equally to this work.</footnote>

<footnote>‡ The two authors contribute equally to this work.</footnote>
<figure><img src="image_4.png" /></figure>

Article

BMAL1 Knockdown Leans Epithelial–Mesenchymal Balance
toward Epithelial Properties and Decreases the
Chemoresistance of Colon Carcinoma Cells

Yuan Zhang 1,2,3,†, Aurore Devocelle 2,4,†, Christophe Desterke 1,2, Lucas Eduardo Botelho de Souza 1,2,
Éva Hadadi 1,2, Hervé Acloque 1,2, Adlen Foudi 1,2 ID, Yao Xiang 5, Annabelle Ballesta 6, Yunhua Chang 1,2,5,‡ ID
and Julien Giron-Michel 2,4,*,‡ ID

<figure><img src="image_5.png" /></figure>

check for
updates

**Citation:** Zhang, Y.; Devocelle, A.; Desterke, C.; de Souza, L.E.B.; Hadadi, É.; Acloque, H.; Foudi, A.; Xiang, Y.; Ballesta, A.; Chang, Y.; et al. BMAL1 Knockdown Leans Epithelial–Mesenchymal Balance toward Epithelial Properties and Decreases the Chemoresistance of Colon Carcinoma Cells. *Int. J. Mol. Sci.* **2021**, *22*, 5247. https://doi.org/10.3390/ijms22105247

Academic Editor: Monica Fedele

Received: 13 April 2021

Accepted: 6 May 2021

Published: 16 May 2021

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_6.png" /></figure>

Copyright: © 2021 by the authors.
Licensee MDPI, Basel, Switzerland.
This article is an open access article
distributed under the terms and
conditions of the Creative Commons
Attribution (CC BY) license (https://
creativecommons.org/licenses/by/
4.0/).

**Abstract:** The circadian clock coordinates biological and physiological functions to day/night cycles. The perturbation of the circadian clock increases cancer risk and affects cancer progression. Here, we studied how BMAL1 knockdown (BMAL1-KD) by shRNA affects the epithelial–mesenchymal transition (EMT), a critical early event in the invasion and metastasis of colorectal carcinoma (CRC). In corresponding to a gene set enrichment analysis, which showed a significant enrichment of EMT and invasive signatures in BMAL1_high CRC patients as compared to BMAL1_low CRC patients, our results revealed that BMAL1 is implicated in keeping the epithelial–mesenchymal equilibrium of CRC cells and influences their capacity of adhesion, migration, invasion, and chemoresistance. Firstly, BMAL1-KD increased the expression of epithelial markers (E-cadherin, CK-20, and EpCAM) but decreased the expression of Twist and mesenchymal markers (N-cadherin and vimentin) in CRC cell lines. Finally, the molecular alterations after BMAL1-KD promoted mesenchymal-to-epithelial transition-like changes mostly appeared in two primary CRC cell lines (i.e., HCT116 and SW480) compared to the metastatic cell line SW620. As a consequence, migration/invasion and drug resistance capacities decreased in HCT116 and SW480 BMAL1-KD cells. Together, BMAL1-KD alerts the delicate equilibrium between epithelial and mesenchymal properties of CRC cell lines, which revealed the crucial role of BMAL1 in EMT-related CRC metastasis and chemoresistance.

**Keywords:** colorectal cancer; epithelial–mesenchymal transition (EMT); circadian clock; BMAL1; metastasis; chemoresistance

**1. Introduction**

Colorectal cancer (CRC) is one of the most common cancers and represents the fourth
leading cause of cancer death worldwide. In 2019, there were 145,600 new cases and
approximately 51,020 deaths in the United States [1]. Metastasis in liver and lung represents
one of the main causes of CRC-related mortality [2]. An insight into the molecular events
driving metastasis, particularly the enhanced invasiveness and therapeutic resistance, is
crucial for developing some novel treatment regimens to combat CRC. One such critical

<footnote n="1">1. INSERM UMR-S 935, CNRS Campus, 94801 Villejuif, France; zyljxmm88@gmail.com (Y.Z.); christophe.desterke@inserm.fr (C.D.); lucasebsouza@usp.br (L.E.B.d.S.); Eva.Hadadi@vub.be (É.H.); herve.acloque@inrae.fr (H.A.); adlen.foudi@inserm.fr (A.F.); yunhua.chang-marchand@inserm.fr (Y.C.)</footnote>

<footnote n="2">2 Orsay-Vallée Campus, Paris-Saclay University, 91190 Gif-sur-Yvette, France; aurore.devocelle@inserm.fr</footnote>

<footnote n="3">3 Institute of Life Sciences, Chongqing Medical University, Chongqing 400016, China</footnote>

<footnote n="4">4 INSERM UMR-S-MD 1197/Ministry of the Armed Forces, Biomedical Research Institute of the Armed Forces (IRBA), Paul-Brousse Hospital Villejuif and CTSA Clamart, 94807 Villejuif, France</footnote>

<footnote n="5">5 INSERM UMR-S 1151, Department of Immunology, Infectiology and Hematology, Institut Necker-Enfants Malades (INEM), Paris Descartes University, CNRS UMR 8253, 75730 Paris, France; yao.xiang@inserm.fr</footnote>

<footnote n="6">6 INSERM UMR-S 900, Institut Curie, MINES ParisTech CBIO, PSL Research University, 92210 Saint-Cloud, France; annabelle.ballesta@curie.fr</footnote>

<footnote>* Correspondence: julien.giron-michel@inserm.fr; Tel.: +33-1-45-59-53-27</footnote>

<footnote>† The two authors contribute equally to this work.</footnote>

<footnote>‡ The two authors contribute equally to this work.</footnote>
event is the epithelial-to-mesenchymal transition (EMT), an early step of CRC metastasis, by which polarized epithelial cells lose their cell polarity and cell–cell adhesion but gain mesenchymal cell properties with enhanced migratory capacity, invasiveness, metastatic potential, and drug resistance [3]. Several deregulated signaling pathways, such as Sonic Hedgehog, transformation growth factor-β (TGF-β), and Wnt pathways, are involved in the initialization and maintenance of EMT [4]. During this process, the expression of E-cadherin, referred to as the “caretaker” of the epithelial phenotype, decreases and results in the disturbance of the E-cadherin/β-catenin complex [5], which is crucial to maintaining epithelial integrity [4]. Moreover, the disturbance of this complex, ultimately, will lead to the nuclear translocation of β-catenin and transcription of EMT-promoting genes through Wnt signaling activation.

The circadian timing system (CTS) coordinates organismal behavior, including physi-
ology and metabolism, over 24 h environmental cycles. In mammals, the suprachiasmatic
nuclei (SCN) of the hypothalamus allows the coordination of circadian rhythms via pe-
ipheral molecular clocks composed of about twenty genes located in each cell of the
organism. The circadian clock genes compose a complex network of auto-regulatory tran-
scription/translation feedback loops, such as positive (BMAL1 and CLOCK) or negative
(PER and CRY) loops and also feedback loops (REV-ERB α/β and ROR α/β/γ) to regulate
different biological rhythm progressions [6]. The circadian clock thus directly or indi-
rectly regulates the expression of thousands of genes in multiple cell types. Subsequently,
these rhythmic regulations control many cellular processes such as autophagy, nutrient
metabolism, redox regulation, cell cycle, DNA damage repair, protein folding, and cellular
secretion [6–9]. An intact circadian clock is therefore essential for proper functioning of
numerous physiological and behavioral processes, and its disruption may lead to disorders
including cancers [10]. Indeed, higher incidence rates of colorectal, breast, and endometrial
cancers have been reported in shift workers with presumed circadian disruption [11]. In
addition, circadian rhythms could significantly modify the efficacy and toxicity of more
than 50 anticancer drugs. Oxaliplatin, a first-line chemotherapeutic anti-colorectal cancer
drug is a good example. At first, its development was halted due to the excessive toxicities,
whereas further clinical chronotherapy trials in phases I, II, and III revealed its anticancer
efficacy by using the optimal drug delivery time which minimized adverse effects and
maximized therapeutic efficacy [12,13].

As a main positive loop element which initiates circadian oscillation, BMAL1 is con-
sidered an essential and non-redundant component of circadian clocks. Among canonical
clock genes, only BMAL1 knockout results in complete loss of circadian rhythm in both
SCN and peripheral tissues [14], which certificates the key role of BMAL1 in the circadian
system. Studies revealed that BMAL1 is involved in the pathogenesis of human cancers,
functioning either as tumor suppressor or oncogenic factor [15–18]. Although BMAL1
exhibits a globally repressive function in many tumors, a high level of BMAL1 and CLOCK
expressions are often associated with poorly differentiated or late-stage CRC as well as
liver metastasis [19,20]. In fact, the role of BMAL1 in CRC development remains poorly
understood, in particular its link to EMT, a key process in tumor progression. Here, we
evaluated in vitro how knockdown BMAL1 (BMAL1-KD) directly affects EMT as well as
migration/invasion and the drug-resistance capacities of two primary human CRC cell
lines (i.e., HCT116 and SW480) and a metastatic CRC cell line (i.e., SW620). Our results
revealed the importance of BMAL1 in modulating the epithelial–mesenchymal equilibrium
of CRC cells.

2. Results

2.1. *Transcriptome of CRC Patients Shows BMAL1 Correlation to EMT and Cancer Invasiveness*

Disruption of the circadian clock plays a crucial role in carcinogenesis, in particular
in human CRC [11–13]. Different from other circadian clock gene knockouts, BMAL1 is
the only gene whose deletion alone leads to complete loss of rhythmicity [14]. To evaluate
the biological functions of BMAL1 during colorectal tumorigenesis, we first investigated
event is the epithelial-to-mesenchymal transition (EMT), an early step of CRC metastasis, by which polarized epithelial cells lose their cell polarity and cell–cell adhesion but gain mesenchymal cell properties with enhanced migratory capacity, invasiveness, metastatic potential, and drug resistance [3]. Several deregulated signaling pathways, such as Sonic Hedgehog, transformation growth factor-β (TGF-β), and Wnt pathways, are involved in the initialization and maintenance of EMT [4]. During this process, the expression of E-cadherin, referred to as the “caretaker” of the epithelial phenotype, decreases and results in the disturbance of the E-cadherin/β-catenin complex [5], which is crucial to maintaining epithelial integrity [4]. Moreover, the disturbance of this complex, ultimately, will lead to the nuclear translocation of β-catenin and transcription of EMT-promoting genes through Wnt signaling activation.

The circadian timing system (CTS) coordinates organismal behavior, including physi-
ology and metabolism, over 24 h environmental cycles. In mammals, the suprachiasmatic
nuclei (SCN) of the hypothalamus allows the coordination of circadian rhythms via pe-
ipheral molecular clocks composed of about twenty genes located in each cell of the
organism. The circadian clock genes compose a complex network of auto-regulatory tran-
scription/translation feedback loops, such as positive (BMAL1 and CLOCK) or negative
(PER and CRY) loops and also feedback loops (REV-ERB α/β and ROR α/β/γ) to regulate
different biological rhythm progressions [6]. The circadian clock thus directly or indi-
rectly regulates the expression of thousands of genes in multiple cell types. Subsequently,
these rhythmic regulations control many cellular processes such as autophagy, nutrient
metabolism, redox regulation, cell cycle, DNA damage repair, protein folding, and cellular
secretion [6–9]. An intact circadian clock is therefore essential for proper functioning of
numerous physiological and behavioral processes, and its disruption may lead to disorders
including cancers [10]. Indeed, higher incidence rates of colorectal, breast, and endometrial
cancers have been reported in shift workers with presumed circadian disruption [11]. In
addition, circadian rhythms could significantly modify the efficacy and toxicity of more
than 50 anticancer drugs. Oxaliplatin, a first-line chemotherapeutic anti-colorectal cancer
drug is a good example. At first, its development was halted due to the excessive toxicities,
whereas further clinical chronotherapy trials in phases I, II, and III revealed its anticancer
efficacy by using the optimal drug delivery time which minimized adverse effects and
maximized therapeutic efficacy [12,13].

As a main positive loop element which initiates circadian oscillation, BMAL1 is con-
sidered an essential and non-redundant component of circadian clocks. Among canonical
clock genes, only BMAL1 knockout results in complete loss of circadian rhythm in both
SCN and peripheral tissues [14], which certificates the key role of BMAL1 in the circadian
system. Studies revealed that BMAL1 is involved in the pathogenesis of human cancers,
functioning either as tumor suppressor or oncogenic factor [15–18]. Although BMAL1
exhibits a globally repressive function in many tumors, a high level of BMAL1 and CLOCK
expressions are often associated with poorly differentiated or late-stage CRC as well as
liver metastasis [19,20]. In fact, the role of BMAL1 in CRC development remains poorly
understood, in particular its link to EMT, a key process in tumor progression. Here, we
evaluated in vitro how knockdown BMAL1 (BMAL1-KD) directly affects EMT as well as
migration/invasion and the drug-resistance capacities of two primary human CRC cell
lines (i.e., HCT116 and SW480) and a metastatic CRC cell line (i.e., SW620). Our results
revealed the importance of BMAL1 in modulating the epithelial–mesenchymal equilibrium
of CRC cells.

2. Results

2.1. *Transcriptome of CRC Patients Shows BMAL1 Correlation to EMT and Cancer Invasiveness*

Disruption of the circadian clock plays a crucial role in carcinogenesis, in particular
in human CRC [11–13]. Different from other circadian clock gene knockouts, BMAL1 is
the only gene whose deletion alone leads to complete loss of rhythmicity [14]. To evaluate
the biological functions of BMAL1 during colorectal tumorigenesis, we first investigated
the transcriptome of CRC patients by RNAseq in The Cancer Genome Atlas (TCGA) consortium. Starting with Z-score matrix quantification, patients were stratified with their extreme BMAL1 quantification (Z-scores over Gaussian distribution in absolute value): six patients were identified with low expression and five patients with a high expression. BMAL1 expression was found to be significantly different between these two groups of patients (p = 6.04 × 10⁻¹⁰, Figure 1A).

<figure><img src="image_6.png" /><figcaption>**Figure 1.** BMAL1 deregulation in colorectal carcinomas (CRC) is associated to epithelial–mesenchymal transition (EMT) and cancer invasiveness perturbations. (A) Boxplot of mRNA expression (Z-scores, RNAseq) in tumors of CRC with patients stratified on their BMAL1 quantification. (B) Enrichment for hallmark EMT gene sets in CRC tumors of patients with a high level of BMAL1 (NES: Normalized Enrichment Score). (C) Principal component analysis performed with EMT-enriched genes found deregulated in CRC tumor samples. (D) Expression heatmap of EMT genes connected to BMAL1 in CRC tumors. (E) Enrichment of multicancer invasiveness signature in CRC tumors of patients with a high level of BMAL1 (NES: Normalized Enrichment Score). (F) Principal component analysis performed with cancer invasiveness-enriched genes found deregulated in CRC tumor samples. (G) Expression heatmap of cancer invasiveness genes connected to BMAL1 in CRC tumors.</figcaption></figure>

Gene set enrichment analysis (GSEA) showed a significant enrichment of epithelial–mesenchymal transition signature in BMAL1_high CRC patients as compared to BMAL1_low CRC patients (Normalized Enrichment Score: NES = -1.99, *p* < 0.001, Figure 1B). Unsupervised principal component analysis with the main EMT-enriched genes (*FUCA1*, *SGCB*, *IL15*, *WIPF1*, *TGFB1*, *ECM2*, *ADM12*, *POSTN*, *PRRX1*, *VIM*, *BASP1*, and *EMP3*) allowed for the discrimination of the subgroup of BMAL1-dependent CRC patients (p = 3.23 × 10⁻⁵, Figure 1C). This stratification of patients with EMT genes was confirmed by expression heatmap with hierarchical clustering (Euclidean distances, Figure 1D). Similarly, GSEA analysis showed an association between BMAL1 expression in CRC patients with multi-
the transcriptome of CRC patients by RNAseq in The Cancer Genome Atlas (TCGA) consortium. Starting with Z-score matrix quantification, patients were stratified with their extreme BMAL1 quantification (Z-scores over Gaussian distribution in absolute value): six patients were identified with low expression and five patients with a high expression. BMAL1 expression was found to be significantly different between these two groups of patients (p = 6.04 × 10⁻¹⁰, Figure 1A).

<figure><img src="image_7.png" /><figcaption>**Figure 1.** BMAL1 deregulation in colorectal carcinomas (CRC) is associated to epithelial–mesenchymal transition (EMT) and cancer invasiveness perturbations. (A) Boxplot of mRNA expression (Z-scores, RNAseq) in tumors of CRC with patients stratified on their BMAL1 quantification. (B) Enrichment for hallmark EMT gene sets in CRC tumors of patients with a high level of BMAL1 (NES: Normalized Enrichment Score). (C) Principal component analysis performed with EMT-enriched genes found deregulated in CRC tumor samples. (D) Expression heatmap of EMT genes connected to BMAL1 in CRC tumors. (E) Enrichment of multicancer invasiveness signature in CRC tumors of patients with a high level of BMAL1 (NES: Normalized Enrichment Score). (F) Principal component analysis performed with cancer invasiveness-enriched genes found deregulated in CRC tumor samples. (G) Expression heatmap of cancer invasiveness genes connected to BMAL1 in CRC tumors.</figcaption></figure>

Gene set enrichment analysis (GSEA) showed a significant enrichment of epithelial–mesenchymal transition signature in BMAL1_high CRC patients as compared to BMAL1_low CRC patients (Normalized Enrichment Score: NES = -1.99, *p* < 0.001, Figure 1B). Unsupervised principal component analysis with the main EMT-enriched genes (*FUCA1*, *SGCB*, *IL15*, *WIPF1*, *TGFB1*, *ECM2*, *ADM12*, *POSTN*, *PRRX1*, *VIM*, *BASP1*, and *EMP3*) allowed for the discrimination of the subgroup of BMAL1-dependent CRC patients (p = 3.23 × 10⁻⁵, Figure 1C). This stratification of patients with EMT genes was confirmed by expression heatmap with hierarchical clustering (Euclidean distances, Figure 1D). Similarly, GSEA analysis showed an association between BMAL1 expression in CRC patients with multi-
cancer invasiveness signature (NES = -2.53, *p* < 0.001, Figure 1E). Unsupervised principal component analysis with the main cancer invasiveness-enriched genes (ASPN, GLT8D2, OLFML2B, CRISPLD2, ADAM12, POSTN, and PRRX1) allowed for the discrimination of subgroups of BMAL1-dependent CRC patients (p = 9.9 × 10⁻⁴, Figure 1F). This stratification of patients with cancer invasiveness genes was confirmed by expression heatmap classification (Euclidean distances, Figure 1G). These results confirmed that BMAL1 over-expression in tumors of CRC patients was associated to EMT and cancer invasiveness program deregulation.

## 2.2. BMAL1-KD Increases E-Cadherin Expression and E-Cadherin/β-Catenin Co-Localization at the Plasma Membrane of CRC Cells

To investigate BMAL1's involvement in the EMT process and invasiveness during colorectal tumorigenesis, two primary CRC cell lines (i.e., HCT116 and SW480) and a metastatic cell line, SW620, were transduced with lentiviruses encoding a scrambled shRNA (shScr) or a shRNA-targeting BMAL1 (shBMAL1). SW480 and SW620 cell lines are derived from the primary tumor and a metastasis site of the same patient, respectively. Previous results from our laboratory showed that BMAL1 mRNA and protein substantially decreased in the three stable BMAL1 knockdown (BMAL1-KD) cell lines compared to their proper control, thereby modulating different cellular activities such as proliferation, apoptosis, and senescence [21].

We here investigated whether BMAL1-KD could also be involved in the EMT process of colon cancer. As a main component of the cell-cell adhesion junctions, E-cadherin plays a central role in maintaining epithelium integrity [22–24]. Since loss of E-cadherin is a key feature of EMT in various cancers, the E-cadherin expression patterns in the three CRC cell lines were first analyzed to assess the impact of BMAL1-KD on the EMT process. Our results showed a significant increase in E-cadherin expression at mRNA (Figure 2A, RT-PCR) and protein (Figure 2B, Western blot) levels in the three BMAL1-KD cell lines compared to their respective control, except for mRNA expression in SW620 cells. This increase was particularly obvious in the BMAL1-KD SW480 cells. In addition, flow cytometry results also revealed a higher E-cadherin cell membrane expression on the BMAL1-KD CRC lines (Figure 2C).

Recruitment of β-catenin and others Catenins to E-cadherin at the cell membrane is essential in the adhesion complexes that play a crucial role for maintaining cell-cell integrity and the epithelial architecture. Disrupting E-cadherin/β-catenin complex, which conducts to a nuclear translocation of β-catenin, affects not only the cell adhesive repertoire, but also Wnt-signaling activation [25–27]. Thus, immunofluorescence (IF) analysis of E-cadherin (red staining) and β-catenin (blue staining) were realized to investigate if increased E-cadherin expression induced by BMAL1-KD altered membrane/nuclear localization of β-catenin in the different CRC cell lines. Consistent with previous flow cytometry analysis, a stronger staining of plasma membrane tethered E-cadherin, particularly at the interface between adjacent cells, was observed in the three BMAL1-KD CRC cell lines compared to their proper controls (Figure 2D). Moreover, an increased co-localization of β-catenin and E-cadherin at the plasma membrane was observed in HCT116 and SW480 BMAL1-KD cells. Indeed, HCT116 BMAL1-KD cells showed a predominantly cell membrane localization of β-catenin compared to its control in which β-catenin exhibited both a cell membrane and cytoplasmic localization. Interestingly, SW480 BMAL1-KD cells showed an evident plasma membrane β-catenin and E-cadherin co-localization compared to control cells in which no cortical staining of β-catenin was observed.

Western blot on cytoplasmic and nuclear extracts of BMAL1-KD and control CRC cell lines were applied to analyze the localization of β-catenin in the different cell fractions (Figure 2E). Consistently with IF experiments, a significant increase in β-catenin expression in the cytoplasmic compartments that contain cytosolic and membranous proteins was found in HCT116 and SW480 BMAL1-KD cells (Figure 2E, upper panels). In addition, both HCT116 and SW480 BMAL1-KD cells presented a significant decrease in β-catenin nuclear
cancer invasiveness signature (NES = -2.53, *p* < 0.001, Figure 1E). Unsupervised principal component analysis with the main cancer invasiveness-enriched genes (ASPN, GLT8D2, OLFML2B, CRISPLD2, ADAM12, POSTN, and PRRX1) allowed for the discrimination of subgroups of BMAL1-dependent CRC patients (p = 9.9 × 10⁻⁴, Figure 1F). This stratification of patients with cancer invasiveness genes was confirmed by expression heatmap classification (Euclidean distances, Figure 1G). These results confirmed that BMAL1 over-expression in tumors of CRC patients was associated to EMT and cancer invasiveness program deregulation.

## 2.2. BMAL1-KD Increases E-Cadherin Expression and E-Cadherin/β-Catenin Co-Localization at the Plasma Membrane of CRC Cells

To investigate BMAL1's involvement in the EMT process and invasiveness during colorectal tumorigenesis, two primary CRC cell lines (i.e., HCT116 and SW480) and a metastatic cell line, SW620, were transduced with lentiviruses encoding a scrambled shRNA (shScr) or a shRNA-targeting BMAL1 (shBMAL1). SW480 and SW620 cell lines are derived from the primary tumor and a metastasis site of the same patient, respectively. Previous results from our laboratory showed that BMAL1 mRNA and protein substantially decreased in the three stable BMAL1 knockdown (BMAL1-KD) cell lines compared to their proper control, thereby modulating different cellular activities such as proliferation, apoptosis, and senescence [21].

We here investigated whether BMAL1-KD could also be involved in the EMT process of colon cancer. As a main component of the cell-cell adhesion junctions, E-cadherin plays a central role in maintaining epithelium integrity [22–24]. Since loss of E-cadherin is a key feature of EMT in various cancers, the E-cadherin expression patterns in the three CRC cell lines were first analyzed to assess the impact of BMAL1-KD on the EMT process. Our results showed a significant increase in E-cadherin expression at mRNA (Figure 2A, RT-PCR) and protein (Figure 2B, Western blot) levels in the three BMAL1-KD cell lines compared to their respective control, except for mRNA expression in SW620 cells. This increase was particularly obvious in the BMAL1-KD SW480 cells. In addition, flow cytometry results also revealed a higher E-cadherin cell membrane expression on the BMAL1-KD CRC lines (Figure 2C).

Recruitment of β-catenin and others Catenins to E-cadherin at the cell membrane is essential in the adhesion complexes that play a crucial role for maintaining cell-cell integrity and the epithelial architecture. Disrupting E-cadherin/β-catenin complex, which conducts to a nuclear translocation of β-catenin, affects not only the cell adhesive repertoire, but also Wnt-signaling activation [25–27]. Thus, immunofluorescence (IF) analysis of E-cadherin (red staining) and β-catenin (blue staining) were realized to investigate if increased E-cadherin expression induced by BMAL1-KD altered membrane/nuclear localization of β-catenin in the different CRC cell lines. Consistent with previous flow cytometry analysis, a stronger staining of plasma membrane tethered E-cadherin, particularly at the interface between adjacent cells, was observed in the three BMAL1-KD CRC cell lines compared to their proper controls (Figure 2D). Moreover, an increased co-localization of β-catenin and E-cadherin at the plasma membrane was observed in HCT116 and SW480 BMAL1-KD cells. Indeed, HCT116 BMAL1-KD cells showed a predominantly cell membrane localization of β-catenin compared to its control in which β-catenin exhibited both a cell membrane and cytoplasmic localization. Interestingly, SW480 BMAL1-KD cells showed an evident plasma membrane β-catenin and E-cadherin co-localization compared to control cells in which no cortical staining of β-catenin was observed.

Western blot on cytoplasmic and nuclear extracts of BMAL1-KD and control CRC cell lines were applied to analyze the localization of β-catenin in the different cell fractions (Figure 2E). Consistently with IF experiments, a significant increase in β-catenin expression in the cytoplasmic compartments that contain cytosolic and membranous proteins was found in HCT116 and SW480 BMAL1-KD cells (Figure 2E, upper panels). In addition, both HCT116 and SW480 BMAL1-KD cells presented a significant decrease in β-catenin nuclear
expression (Figure 2E, lower panels). No evident difference of β-catenin localization was observed between SW620 BMAL1-KD and control cells, neither in IF, nor in Western blot.

<figure><img src="image_7.png" /></figure>

**Figure 2.** BMAL1-KD increased the plasma membrane colocation of E-cadherin and β-catenin but decreased β-catenin nuclear location. (A) Quantitative RT-PCR assay revealed an increased E-cadherin expression in two primary BMAL1 knockdown (BMAL1-KD) cell lines (n = 4, *p* < 0.05, **p** < 0.01) in comparison to control (Ctr) cells. (B) Increased E-cadherin expression in BMAL1-KD CRC cell lines was checked by Western blot. Left, Representative Western blots of 6 independent experiments are shown. Right, Bar charts represents E-cadherin expression normalized to HSC70 (n = 6, *p* < 0.05, **p** < 0.01). All data are shown as means ± SEM. (C) Enhanced E-cadherin surface expression on BMAL1-KD CRC cell lines was ascertained by flow cytometry analysis. Left, Representative staining of one experiment are shown. Cell staining with the anti-E-Cadherin antibody is represented by the dark-colored histograms and staining with the isotype-matched control antibody correspond to the gray histograms. At the right of each histogram are shown the mean fluorescence intensity values for each staining. Right, Graphs represent the relative expression of E-cadherin in three independent experiments (*p* < 0.05, ***p* < 0.001). (D) Confocal laser scanning microscopy showing stainings for E-cadherin (red) and β-catenin (blue). BMAL1-KD promotes elevated levels of E-cadherin and β-catenin at the plasma membrane of both HCT116 and SW480 cells, whereas BMAL1-KD only increased evidently the membrane localization of E-cadherin but not β-catenin in SW620 cells. Scale bars, 20 µm. These data are representative of three independent experiments. (E) Cytoplasmic (upper panels) and nuclear (lower panels) extracts of BMAL1-KD and control (Ctr) cell lines were analyzed by Western blot. Left, Representative stainings of one experiment are shown. Right, Graphs represent the quantification of Western blot signals, corresponding to β-catenin expression normalized to HSC70 (cytoplasmic extracts) or Lamin β1 (nuclear extracts). Data are represented as means ± SEM for 5 independent determinations (*p* < 0.05). A significant increase in β-catenin cytoplasmic localization was found in SW480 BMAL1-KD cells. HCT116 BMAL1-KD cells presented an increased tendency of cytoplasmic β-catenin (p = 0.056). A significant decrease in β-catenin nuclear localization (n = 5, *p* < 0.05) was found in the SW480 BMAL1-KD and HCT116 BMAL1-KD cell lines.
expression (Figure 2E, lower panels). No evident difference of β-catenin localization was observed between SW620 BMAL1-KD and control cells, neither in IF, nor in Western blot.

<figure><img src="image_8.png" /></figure>

**Figure 2.** BMAL1-KD increased the plasma membrane colocation of E-cadherin and β-catenin but decreased β-catenin nuclear location. (A) Quantitative RT-PCR assay revealed an increased E-cadherin expression in two primary BMAL1 knockdown (BMAL1-KD) cell lines (n = 4, *p* < 0.05, **p** < 0.01) in comparison to control (Ctr) cells. (B) Increased E-cadherin expression in BMAL1-KD CRC cell lines was checked by Western blot. Left, Representative Western blots of 6 independent experiments are shown. Right, Bar charts represents E-cadherin expression normalized to HSC70 (n = 6, *p* < 0.05, **p** < 0.01). All data are shown as means ± SEM. (C) Enhanced E-cadherin surface expression on BMAL1-KD CRC cell lines was ascertained by flow cytometry analysis. Left, Representative staining of one experiment are shown. Cell staining with the anti-E-Cadherin antibody is represented by the dark-colored histograms and staining with the isotype-matched control antibody correspond to the gray histograms. At the right of each histogram are shown the mean fluorescence intensity values for each staining. Right, Graphs represent the relative expression of E-cadherin in three independent experiments (*p* < 0.05, ***p* < 0.001). (D) Confocal laser scanning microscopy showing stainings for E-cadherin (red) and β-catenin (blue). BMAL1-KD promotes elevated levels of E-cadherin and β-catenin at the plasma membrane of both HCT116 and SW480 cells, whereas BMAL1-KD only increased evidently the membrane localization of E-cadherin but not β-catenin in SW620 cells. Scale bars, 20 µm. These data are representative of three independent experiments. (E) Cytoplasmic (upper panels) and nuclear (lower panels) extracts of BMAL1-KD and control (Ctr) cell lines were analyzed by Western blot. Left, Representative stainings of one experiment are shown. Right, Graphs represent the quantification of Western blot signals, corresponding to β-catenin expression normalized to HSC70 (cytoplasmic extracts) or Lamin β1 (nuclear extracts). Data are represented as means ± SEM for 5 independent determinations (*p* < 0.05). A significant increase in β-catenin cytoplasmic localization was found in SW480 BMAL1-KD cells. HCT116 BMAL1-KD cells presented an increased tendency of cytoplasmic β-catenin (p = 0.056). A significant decrease in β-catenin nuclear localization (n = 5, *p* < 0.05) was found in the SW480 BMAL1-KD and HCT116 BMAL1-KD cell lines.
### 2.3. BMAL1-KD Leans Epithelial–Mesenchymal Balance of CRC Cells toward Epithelial Properties

Downregulation of E-cadherin is considered as a major hallmark of EMT as restoring E-cadherin expression is sufficient for reversal of the transformed phenotype [3]. Increased E-cadherin expression in BMAL1-KD CRC cells prompted us to check if the expression of other epithelial and mesenchymal markers were also altered. Western blot (Figure 3A) and cell surface staining in flow cytometry (Figure 3B) results revealed higher levels of EpCAM, another key epithelial denominator, in BMAL1-KD CRC cell lines as compared to controls. Moreover, the total (Figure 3A) and intracellular (Figure 3B) expressions of the epithelial marker Cytokeratin 20 (CK-20) were enhanced while those of the mesenchymal marker vimentin were decreased in BMAL1-KD CRC cell lines.

<figure><img src="image_8.png" /><figcaption>**Figure 3.** BMAL1-KD leans epithelial–mesenchymal balance of CRC cells toward epithelial properties. (A) Western blot revealed that BMAL1-KD increased the expression of the epithelial markers EpCAM and cytokeratin-20 (CK-20) and decreased the mesenchymal marker (vimentin) in BMAL1-KD CRC cell lines. HSC70 was used as loading controls. (B) EpCAM surface expression and intracellular expression of vimentin and CK-20 were analyzed by flow cytometry. Left, Representative stainings of one experiment are shown. Staining of the cells with the antigen-specific antibody is represented by the dark-colored histograms and staining with the isotype-matched control antibody correspond to the gray histograms. At the right of each histogram are shown the mean fluorescence intensity values for each staining. Right, Graphs represent the relative expression of EpCAM, CK-20, and vimentin in three independent experiments (*p* < 0.05, **p** < 0.01, ***p** < 0.001). (C) Quantitative RT-PCR revealed that RNA expression level of mesenchymal hallmark N-cadherin was significantly decreased in SW480 BMAL1-KD and SW620 BMAL1-KD cell lines. The BMAL1-KD HCT116 cell line exhibited a decreased tendency of N-cadherin (p = 0.0762) (n = 5, **p** < 0.01). Data are shown as means ± SEM. (D) Quantitative RT-PCR analyzed the RNA expression level of three EMT-activating transcription factors Twist, Slug, and Snail in control and BMAL1-KD CRC cell lines. BMAL1-KD decreased significantly the expression of Twist but not Slug and Snail (n = 5, *p* < 0.05, **p** < 0.01). Data are shown as means ± SEM.</figcaption></figure>
### 2.3. BMAL1-KD Leans Epithelial–Mesenchymal Balance of CRC Cells toward Epithelial Properties

Downregulation of E-cadherin is considered as a major hallmark of EMT as restoring E-cadherin expression is sufficient for reversal of the transformed phenotype [3]. Increased E-cadherin expression in BMAL1-KD CRC cells prompted us to check if the expression of other epithelial and mesenchymal markers were also altered. Western blot (Figure 3A) and cell surface staining in flow cytometry (Figure 3B) results revealed higher levels of EpCAM, another key epithelial denominator, in BMAL1-KD CRC cell lines as compared to controls. Moreover, the total (Figure 3A) and intracellular (Figure 3B) expressions of the epithelial marker Cytokeratin 20 (CK-20) were enhanced while those of the mesenchymal marker vimentin were decreased in BMAL1-KD CRC cell lines.

<figure><img src="image_9.png" /><figcaption>**Figure 3.** BMAL1-KD leans epithelial–mesenchymal balance of CRC cells toward epithelial properties. (A) Western blot revealed that BMAL1-KD increased the expression of the epithelial markers EpCAM and cytokeratin-20 (CK-20) and decreased the mesenchymal marker (vimentin) in BMAL1-KD CRC cell lines. HSC70 was used as loading controls. (B) EpCAM surface expression and intracellular expression of vimentin and CK-20 were analyzed by flow cytometry. Left, Representative stainings of one experiment are shown. Staining of the cells with the antigen-specific antibody is represented by the dark-colored histograms and staining with the isotype-matched control antibody correspond to the gray histograms. At the right of each histogram are shown the mean fluorescence intensity values for each staining. Right, Graphs represent the relative expression of EpCAM, CK-20, and vimentin in three independent experiments (*p* < 0.05, **p** < 0.01, ***p** < 0.001). (C) Quantitative RT-PCR revealed that RNA expression level of mesenchymal hallmark N-cadherin was significantly decreased in SW480 BMAL1-KD and SW620 BMAL1-KD cell lines. The BMAL1-KD HCT116 cell line exhibited a decreased tendency of N-cadherin (p = 0.0762) (n = 5, **p** < 0.01). Data are shown as means ± SEM. (D) Quantitative RT-PCR analyzed the RNA expression level of three EMT-activating transcription factors Twist, Slug, and Snail in control and BMAL1-KD CRC cell lines. BMAL1-KD decreased significantly the expression of Twist but not Slug and Snail (n = 5, *p* < 0.05, **p** < 0.01). Data are shown as means ± SEM.</figcaption></figure>
Furthermore, quantitative RT-PCR revealed a significant decrease in the expression of the mesenchymal marker N-cadherin in SW480 and SW620 BMAL1-KD cells, this gene being expressed at very low levels in HCT116 cell lines (Figure 3C). Flow cytometry analysis confirmed a decreased surface expression of N-cadherin on the three BMAL1-KD cell lines (Figure S1). Finally, quantitative RT-PCR results showed that the expression of the EMT transcription factor (EMT-TF) Twist decreased in the three BMAL1-KD CRC cell lines compared to their respective control, whereas no significant changes in the mRNA levels of the two other key EMT-TFs, Slug and Snail, were observed (Figure 3D).

## 2.4. BMAL1-KD Induces Morphological Changes in CRC Cell Lines

Cells undergoing EMT or its reverse process MET, exhibited profound morphological changes. As BMAL1-KD leans the epithelial–mesenchymal balance of the colon carcinoma cell toward epithelial properties, we investigated the cell morphology of the three CRC cell lines by phase contrast microscopy (Figure 4A). HCT116 monolayer cell cultures showed typical colonies with an epithelial cobblestone-like phenotype, whereas SW620 cells display an ovoid morphology and form small aggregates (left panels). Interestingly, HCT116 and SW620 BMAL1-KD cells appeared to grow as more densely packed clones with significantly enhanced cell-to-cell contacts (right panels). The most prominent morphological changes were especially noted in SW480 BMAL1-KD cells. Indeed, BMAL1-KD induced an obvious phenotypic change in SW480 cells which lost their mesenchymal-like fusiform morphology, characteristic of a loosely connected mesenchymal phenotype (control cells) to a more tightly associated rounded epithelial phenotype (BMAL1-KD).

<figure><img src="image_9.png" /><figcaption>**Figure 4.** BMAL1-KD induced morphological changes on CRC cell lines. (A) Cells were viewed using phase contrast microscopy at the objective 20×. Compared to the control (Ctr) cells, CRC cells formed densely packed clones in the three BMAL1-KD CRC cell lines (white arrows). Morphological changes were most prominent in SW480 BMAL1-KD cells, with loss of an elongated and spindle-shaped morphology and acquisition of a typical epithelial shape with a cobblestone-like morphology. Scale bars, 50 µm. (B) F-actin distribution profile was analyzed by fluorescence microscopy with phalloidin (green). Nuclei staining with DAPI (blue) was inserted at the top-left of each image. Scale bars, 15 µm. The F-actin distribution profile was similar to those of E-cadherin and β-catenin in BMAL1-KD CRC cells, revealing a specific honeycomb-like epithelial organization of the adhesion belts delineated by E-cadherin, β-catenin, and F-actin.</figcaption></figure>
Furthermore, quantitative RT-PCR revealed a significant decrease in the expression of the mesenchymal marker N-cadherin in SW480 and SW620 BMAL1-KD cells, this gene being expressed at very low levels in HCT116 cell lines (Figure 3C). Flow cytometry analysis confirmed a decreased surface expression of N-cadherin on the three BMAL1-KD cell lines (Figure S1). Finally, quantitative RT-PCR results showed that the expression of the EMT transcription factor (EMT-TF) Twist decreased in the three BMAL1-KD CRC cell lines compared to their respective control, whereas no significant changes in the mRNA levels of the two other key EMT-TFs, Slug and Snail, were observed (Figure 3D).

## 2.4. BMAL1-KD Induces Morphological Changes in CRC Cell Lines

Cells undergoing EMT or its reverse process MET, exhibited profound morphological changes. As BMAL1-KD leans the epithelial–mesenchymal balance of the colon carcinoma cell toward epithelial properties, we investigated the cell morphology of the three CRC cell lines by phase contrast microscopy (Figure 4A). HCT116 monolayer cell cultures showed typical colonies with an epithelial cobblestone-like phenotype, whereas SW620 cells display an ovoid morphology and form small aggregates (left panels). Interestingly, HCT116 and SW620 BMAL1-KD cells appeared to grow as more densely packed clones with significantly enhanced cell-to-cell contacts (right panels). The most prominent morphological changes were especially noted in SW480 BMAL1-KD cells. Indeed, BMAL1-KD induced an obvious phenotypic change in SW480 cells which lost their mesenchymal-like fusiform morphology, characteristic of a loosely connected mesenchymal phenotype (control cells) to a more tightly associated rounded epithelial phenotype (BMAL1-KD).

<figure><img src="image_10.png" /><figcaption>**Figure 4.** BMAL1-KD induced morphological changes on CRC cell lines. (A) Cells were viewed using phase contrast microscopy at the objective 20×. Compared to the control (Ctr) cells, CRC cells formed densely packed clones in the three BMAL1-KD CRC cell lines (white arrows). Morphological changes were most prominent in SW480 BMAL1-KD cells, with loss of an elongated and spindle-shaped morphology and acquisition of a typical epithelial shape with a cobblestone-like morphology. Scale bars, 50 µm. (B) F-actin distribution profile was analyzed by fluorescence microscopy with phalloidin (green). Nuclei staining with DAPI (blue) was inserted at the top-left of each image. Scale bars, 15 µm. The F-actin distribution profile was similar to those of E-cadherin and β-catenin in BMAL1-KD CRC cells, revealing a specific honeycomb-like epithelial organization of the adhesion belts delineated by E-cadherin, β-catenin, and F-actin.</figcaption></figure>
Besides modified cell-cell junctions, epithelialization is associated also with actin cytoskeletal re-organization from stress fibers toward a cortical actin ring [28,29]. We then examined whether BMAL1 could remodel the actin cytoskeleton and affect F-actin distribution in these cell models. An increased formation of cortical actin rings was observed in all BMAL1-KD CRC cell lines but mostly evident in SW480 BMAL1-KD cells. Notably, F-actin was relatively evenly distributed in stress fibers in control of SW480 cells, whereas BMAL1-KD led to a decrease in cell size with an increased cortical distribution of F-actin (Figure 4B). All these results indicated that BMAL1-KD induces loss of mesenchymal traits and reinforces the epithelial phenotype in CRC cells.

## 2.5. BMAL1-KD Reduces Cell Migration and the Invasive Potential of CRC Cell Lines

As EMT promotes invasiveness of cancer cells [30], a scratch-wound healing assay and an in vitro Matrigel invasion assay were performed to check if the leaned epithelial-mesenchymal balance induced by BMAL1-KD could affect the migration and invasion properties of CRC cells. For the scratch-wound healing assay, a wound was generated in a confluent monolayer of cells. The cell's ability to migrate into the wound area was evaluated 24 h to 72 h after the wound by phase contrast microscopy. The HCT116 and SW480 BMAL1-KD cells displayed a significantly lower capacity of migration at 48 h and 72 h than their control cell lines (Figure 5A). On the opposite, wound healing was not significantly different between the control and BMAL1-KD SW620 cells ($p > 0.05$).

For the in vitro Matrigel invasion assay, cells were added to Matrigel-coated Boyden chambers and allowed to migrate across the matrix along a serum gradient (Figure 5B). In agreement with the wound healing assay, BMAL1-KD resulted in significantly lower levels of invasion in HCT116 and SW480 cell lines. However, the SW620 BMAL1-KD cell line showed no significant differences in its invasive properties compared to its control. Taken altogether, these results suggest that the BMAL1-KD primary CRC cell lines presented decreased migratory ability and invasiveness properties compared to the controls. Consistent with reinforced cell-cell junction integrity, these results provided more evidence that BMAL1-KD leaned epithelial-mesenchymal the balance of the primary colon carcinoma cell toward epithelial properties.

## 2.6. BMAL1-KD Increases the Expression of Cell Adhesion Molecules in CRC Cell Lines

Next, we further investigated whether BMAL1-KD could influence the expression of cell adhesion molecules, such as integrin β1 (CD29), integrin alpha5 (CD49e), and CD44, that were involved in the matrix-cell adhesion [31,32] and the apico-basolateral polarity [33] of CRC cells (Figure 6). Flow cytometry analysis showed that the downregulation of BMAL1 was associated with higher membranous levels of integrin β1, integrin alpha5, and CD44 in the two primary CRC cell lines (i.e., HCT116 and SW480). However, similar membranous levels of CD44, integrin β1, and integrin alpha5 were observed in SW620 cells before and after BMAL1-KD.

## 2.7. Response to Standard Therapy Is Increased in BMAL1-KD CRC Cell Lines

In addition to the migratory and invasive abilities, recent studies have shown that the EMT phenotype is closely related to chemoresistance of cancer cells [3,34–37]. We therefore studied whether BMAL1-KD affects the resistance of CRC cells to oxaliplatin (LOH), the third-generation platinum drug widely used as a first-line chemotherapeutic agent in patients with metastatic CRC. BMAL1-KD or the control CRC cell lines were treated with different concentrations of LOH (12.5–75 µM) for 48 h (Figure 7). The MTT assay was used to evaluate cell viability after LOH treatment. The results showed that BMAL1-KD led to a distinct decreased cell viability of HCT116 and SW480 BMAL1-KD cell lines compared to their proper control. However, no differences were noted between the control and the SW620 BMAL1-KD cells.
Besides modified cell-cell junctions, epithelialization is associated also with actin cytoskeletal re-organization from stress fibers toward a cortical actin ring [28,29]. We then examined whether BMAL1 could remodel the actin cytoskeleton and affect F-actin distribution in these cell models. An increased formation of cortical actin rings was observed in all BMAL1-KD CRC cell lines but mostly evident in SW480 BMAL1-KD cells. Notably, F-actin was relatively evenly distributed in stress fibers in control of SW480 cells, whereas BMAL1-KD led to a decrease in cell size with an increased cortical distribution of F-actin (Figure 4B). All these results indicated that BMAL1-KD induces loss of mesenchymal traits and reinforces the epithelial phenotype in CRC cells.

## 2.5. BMAL1-KD Reduces Cell Migration and the Invasive Potential of CRC Cell Lines

As EMT promotes invasiveness of cancer cells [30], a scratch-wound healing assay and an in vitro Matrigel invasion assay were performed to check if the leaned epithelial-mesenchymal balance induced by BMAL1-KD could affect the migration and invasion properties of CRC cells. For the scratch-wound healing assay, a wound was generated in a confluent monolayer of cells. The cell's ability to migrate into the wound area was evaluated 24 h to 72 h after the wound by phase contrast microscopy. The HCT116 and SW480 BMAL1-KD cells displayed a significantly lower capacity of migration at 48 h and 72 h than their control cell lines (Figure 5A). On the opposite, wound healing was not significantly different between the control and BMAL1-KD SW620 cells ($p > 0.05$).

For the in vitro Matrigel invasion assay, cells were added to Matrigel-coated Boyden chambers and allowed to migrate across the matrix along a serum gradient (Figure 5B). In agreement with the wound healing assay, BMAL1-KD resulted in significantly lower levels of invasion in HCT116 and SW480 cell lines. However, the SW620 BMAL1-KD cell line showed no significant differences in its invasive properties compared to its control. Taken altogether, these results suggest that the BMAL1-KD primary CRC cell lines presented decreased migratory ability and invasiveness properties compared to the controls. Consistent with reinforced cell-cell junction integrity, these results provided more evidence that BMAL1-KD leaned epithelial-mesenchymal the balance of the primary colon carcinoma cell toward epithelial properties.

## 2.6. BMAL1-KD Increases the Expression of Cell Adhesion Molecules in CRC Cell Lines

Next, we further investigated whether BMAL1-KD could influence the expression of cell adhesion molecules, such as integrin β1 (CD29), integrin alpha5 (CD49e), and CD44, that were involved in the matrix-cell adhesion [31,32] and the apico-basolateral polarity [33] of CRC cells (Figure 6). Flow cytometry analysis showed that the downregulation of BMAL1 was associated with higher membranous levels of integrin β1, integrin alpha5, and CD44 in the two primary CRC cell lines (i.e., HCT116 and SW480). However, similar membranous levels of CD44, integrin β1, and integrin alpha5 were observed in SW620 cells before and after BMAL1-KD.

## 2.7. Response to Standard Therapy Is Increased in BMAL1-KD CRC Cell Lines

In addition to the migratory and invasive abilities, recent studies have shown that the EMT phenotype is closely related to chemoresistance of cancer cells [3,34–37]. We therefore studied whether BMAL1-KD affects the resistance of CRC cells to oxaliplatin (LOH), the third-generation platinum drug widely used as a first-line chemotherapeutic agent in patients with metastatic CRC. BMAL1-KD or the control CRC cell lines were treated with different concentrations of LOH (12.5–75 µM) for 48 h (Figure 7). The MTT assay was used to evaluate cell viability after LOH treatment. The results showed that BMAL1-KD led to a distinct decreased cell viability of HCT116 and SW480 BMAL1-KD cell lines compared to their proper control. However, no differences were noted between the control and the SW620 BMAL1-KD cells.
<figure><img src="image_10.png" /><figcaption>**Figure 5.** BMAL1-KD inhibits cell migration and invasion of CRC cell lines. (A) A scratch-wound healing assay was applied for cell migration assay. Top, Artificial wounds were made in confluent monolayers of control and BMAL1-KD CRC cell lines. Migration of CRC cells towards the wound, photographed from 24 h to 72 h. The top panel shows one of three independent experiments. Bottom, HCT116 BMAL1-KD and SW480 BMAL1-KD cells displayed lower levels of migratory activity in comparison to control cell lines (n = 3, *p* < 0.05, **p* < 0.01, ***p* < 0.001). No significant differences in migratory activity were observed between BMAL1-KD and control SW620 cells. Data are shown as means ± SEM. (B) Invasion assay. Top, The invasive potential of BMAL1-KD CRC cell lines and their control were analyzed in Matrigel Transwells after 96 h by fluorescence microscopy. DAPI was used to stain the nuclei and to determine the number of invasive cells. A representative image of random fields is shown for each cell type. Scale bar, 50 µm. Bottom, Graphs represent the mean of the number of invaded cells per field in three independent experiments. Only the primary BMAL1-KD CRC cells had a lower invasion capacity than control cells (n = 3, ***p* < 0.001). Data are shown as means ± SEM.</figcaption></figure>
<figure><img src="image_11.png" /><figcaption>**Figure 5.** BMAL1-KD inhibits cell migration and invasion of CRC cell lines. (A) A scratch-wound healing assay was applied for cell migration assay. Top, Artificial wounds were made in confluent monolayers of control and BMAL1-KD CRC cell lines. Migration of CRC cells towards the wound, photographed from 24 h to 72 h. The top panel shows one of three independent experiments. Bottom, HCT116 BMAL1-KD and SW480 BMAL1-KD cells displayed lower levels of migratory activity in comparison to control cell lines (n = 3, *p* < 0.05, **p* < 0.01, ***p* < 0.001). No significant differences in migratory activity were observed between BMAL1-KD and control SW620 cells. Data are shown as means ± SEM. (B) Invasion assay. Top, The invasive potential of BMAL1-KD CRC cell lines and their control were analyzed in Matrigel Transwells after 96 h by fluorescence microscopy. DAPI was used to stain the nuclei and to determine the number of invasive cells. A representative image of random fields is shown for each cell type. Scale bar, 50 µm. Bottom, Graphs represent the mean of the number of invaded cells per field in three independent experiments. Only the primary BMAL1-KD CRC cells had a lower invasion capacity than control cells (n = 3, ***p* < 0.001). Data are shown as means ± SEM.</figcaption></figure>
<figure><img src="image_11.png" /><figcaption>**Figure 6.** BMAL1-KD increased the expression of cell adhesion molecules on CRC cell lines. Flow cytometry analysis of surface expression of CD29 (integrin β1), CD49e (integrin α5), and CD44 adhesion molecules. Left, a representative direct staining of four independent experiments is shown. The dark-colored bars correspond to cells incubated with phycoerythrin (PE)-conjugated antibodies against CD29, CD49e, and CD44, and grey bars correspond to cells incubated with the isotype-matched control antibody. Mean fluorescence intensity values for each marker are shown at the right of each histogram. Right, BMAL1-KD led to increased expression of membranous CD29, CD49e, and CD44, mostly in primary BMAL1-KD CRC cell lines (n = 3, *p* < 0.05, **p** < 0.01, ***p** < 0.001). Data are shown as means ± SEM.</figcaption></figure>

<figure><img src="image_12.png" /><figcaption>**Figure 7.** BMAL1-KD sensitizes CRC cell lines to oxaliplatin treatment. The sensitivity of CRC cells to oxaliplatin (LOH) drug was evaluated by the MTT cell viability assay. BMAL1-KD and their control cells were treated for 24 h to 48 h with increasing concentrations of LOH (from 12.5 to 75 µM). Cell viability after drug treatment is presented as a percentage relative to untreated cells. The mean values for three experiments are shown; error bars correspond to 95% confidence intervals. BMAL1-KD led to a distinct decrease in cell viability of HCT116 and SW480 cell lines under LOH treatment compared to their control cells, whereas no significant differences were observed with the metastasis SW620 cell line (n = 3, *p* < 0.05, **p** < 0.01, ***p** < 0.001). Data are shown as means ± SEM.</figcaption></figure>
<figure><img src="image_12.png" /><figcaption>**Figure 6.** BMAL1-KD increased the expression of cell adhesion molecules on CRC cell lines. Flow cytometry analysis of surface expression of CD29 (integrin β1), CD49e (integrin α5), and CD44 adhesion molecules. Left, a representative direct staining of four independent experiments is shown. The dark-colored bars correspond to cells incubated with phycoerythrin (PE)-conjugated antibodies against CD29, CD49e, and CD44, and grey bars correspond to cells incubated with the isotype-matched control antibody. Mean fluorescence intensity values for each marker are shown at the right of each histogram. Right, BMAL1-KD led to increased expression of membranous CD29, CD49e, and CD44, mostly in primary BMAL1-KD CRC cell lines (n = 3, *p* < 0.05, **p* < 0.01, ***p* < 0.001). Data are shown as means ± SEM.</figcaption></figure>

<figure><img src="image_13.png" /><figcaption>**Figure 7.** BMAL1-KD sensitizes CRC cell lines to oxaliplatin treatment. The sensitivity of CRC cells to oxaliplatin (LOH) drug was evaluated by the MTT cell viability assay. BMAL1-KD and their control cells were treated for 24 h to 48 h with increasing concentrations of LOH (from 12.5 to 75 µM). Cell viability after drug treatment is presented as a percentage relative to untreated cells. The mean values for three experiments are shown; error bars correspond to 95% confidence intervals. BMAL1-KD led to a distinct decrease in cell viability of HCT116 and SW480 cell lines under LOH treatment compared to their control cells, whereas no significant differences were observed with the metastasis SW620 cell line (n = 3, *p* < 0.05, **p* < 0.01, ***p* < 0.001). Data are shown as means ± SEM.</figcaption></figure>
### 3. Discussion

Epithelial-to-mesenchymal transition (EMT) is a key process driving CRC metastasis and chemoresistance. In this report, we provided evidence, for the first time, that BMAL1-KD promotes mesenchymal-to-epithelial transition (MET)-like changes of colorectal carcinoma cell lines. Consistently, several typical mesenchymal markers, such as N-cadherin and vimentin, as well as the EMT transcription factor Twist decreased, while several classic epithelial markers, such as E-cadherin, CK-20, and EpCAM, increased. Among all the molecular markers implicated in EMT, E-cadherin is referred to as the “caretaker” of the epithelial phenotype, and its loss is considered a key feature of the EMT [22–24]. First, E-cadherin associated with β-catenin plays a critical role in adherent junctions’ formation and epithelial integrity maintenance. In addition, the E-cadherin/β-catenin complex also plays a crucial role in modulating Wnt signaling since, by maintaining the integrity of epithelial cell-cell contacts, it allows to keep Wnt/β-catenin signals in check. Thus, disruption of the E-cadherin/β-catenin complex leads to β-catenin release and subsequent translocation to the nucleus, where it promotes the transcription of EMT-associated genes upon activation of Wnt signaling [26,27]. The nuclear accumulation of β-catenin, occurring at the invasive front of CRC, is an important hallmark of EMT [38]. Conversely, E-cadherin overexpression in cancer cells impedes their growth and metastasis [5,22]. Our study revealed that BMAL1-KD reinforced epithelial properties of CRC cells with elevated plasma membrane co-localization of E-cadherin and β-catenin as well as a decreased β-catenin nuclear location, suggesting a decreased Wnt/β-catenin pathway activity in BMAL1-KD CRC cell lines. This modification of β-catenin location is mostly evident in SW480 BMAL1-KD cell line, which could be due to two reasons: firstly, BMAL1-KD increased E-cadherin RNA and protein levels to a greater extent in SW480 when compared to HCT116 or SW620 cells. Moreover, SW480 cells express a truncated adenomatous polyposis coli (APC) protein that is deficient in its ability to promote β-catenin degradation [39,40]. Without efficient degradation induced by APC, the cytoplasmic (including cytosolic and membranous) accumulation of β-catenin is more influenced by increased E-cadherin expression and their membrane co-localization. Collectively, these elements may explain why SW480 cells presented a more evident modification of β-catenin distribution in cytoplasm and nuclear after BMAL1-KD.

Acquisition of mesenchymal phenotypes destabilizes cell-cell adhesive junctions and affects cell migratory and invasive properties [5]. Thus, the expression of cell adhesion molecules, such as E-cadherin, EpCAM, CD44, and CD166 (ALCAM), are lost at the invasive front of CRC. Our results revealed that BMAL1-KD, in addition to E-cadherin and EpCAM, increased membranous expression of several cell surface proteins implicated in the matrix-cell adhesion [31,32] and the apico-basolateral polarity [36] of CRC cells such as integrin β1 (CD29), integrin α5 (CD49e), and CD44. Thus, the reduction in migratory and invasive capacity of BMAL1-KD CRC cells should be due not only to the higher expression of epithelial key denominators (E-cadherin and EpCAM) but also because of adhesion molecules such as CD29, CD49e, and CD44. In support of this notion, it was reported that induction of CD29/CD49e clustering reduces invasion and restores apicobasolateral polarity in invasive CRC cells [33]. Moreover, CD29 plays a crucial role in supporting tumor cell attachment to the extracellular matrix. A recent study demonstrated that miR-10a silencing suppresses migration and invasion in vitro, promoting MET-like changes in SW480 cells that lead to enhanced cell-matrix adhesion induced by CD29 [31]. Moreover, inducing MET-like changes in CRC cells, BMAL1-KD reduced both migration and invasion capacities. These results are consistent with other studies showing that knocking down CLOCK, the important partner of BMAL1, led to the impairment of CRC cell migration, while CLOCK overexpression promoted their migration [41]. In addition, our immunofluorescence results in BMAL1-KD CRC cells showed a specific honeycomb-like epithelial organization of the adhesion belts delineated by E-cadherin, β-catenin, and F-actin, a phenomenon associated with the epithelialization. It was shown that BMAL1 and CLOCK could stabilize and activate RhoA [41]. Constitutive activated RhoA could enhance cortical disassembly in
### 3. Discussion

Epithelial-to-mesenchymal transition (EMT) is a key process driving CRC metastasis and chemoresistance. In this report, we provided evidence, for the first time, that BMAL1-KD promotes mesenchymal-to-epithelial transition (MET)-like changes of colorectal carcinoma cell lines. Consistently, several typical mesenchymal markers, such as N-cadherin and vimentin, as well as the EMT transcription factor Twist decreased, while several classic epithelial markers, such as E-cadherin, CK-20, and EpCAM, increased. Among all the molecular markers implicated in EMT, E-cadherin is referred to as the “caretaker” of the epithelial phenotype, and its loss is considered a key feature of the EMT [22–24]. First, E-cadherin associated with β-catenin plays a critical role in adherent junctions’ formation and epithelial integrity maintenance. In addition, the E-cadherin/β-catenin complex also plays a crucial role in modulating Wnt signaling since, by maintaining the integrity of epithelial cell-cell contacts, it allows to keep Wnt/β-catenin signals in check. Thus, disruption of the E-cadherin/β-catenin complex leads to β-catenin release and subsequent translocation to the nucleus, where it promotes the transcription of EMT-associated genes upon activation of Wnt signaling [26,27]. The nuclear accumulation of β-catenin, occurring at the invasive front of CRC, is an important hallmark of EMT [38]. Conversely, E-cadherin overexpression in cancer cells impedes their growth and metastasis [5,22]. Our study revealed that BMAL1-KD reinforced epithelial properties of CRC cells with elevated plasma membrane co-localization of E-cadherin and β-catenin as well as a decreased β-catenin nuclear location, suggesting a decreased Wnt/β-catenin pathway activity in BMAL1-KD CRC cell lines. This modification of β-catenin location is mostly evident in SW480 BMAL1-KD cell line, which could be due to two reasons: firstly, BMAL1-KD increased E-cadherin RNA and protein levels to a greater extent in SW480 when compared to HCT116 or SW620 cells. Moreover, SW480 cells express a truncated adenomatous polyposis coli (APC) protein that is deficient in its ability to promote β-catenin degradation [39,40]. Without efficient degradation induced by APC, the cytoplasmic (including cytosolic and membranous) accumulation of β-catenin is more influenced by increased E-cadherin expression and their membrane co-localization. Collectively, these elements may explain why SW480 cells presented a more evident modification of β-catenin distribution in cytoplasm and nuclear after BMAL1-KD.

Acquisition of mesenchymal phenotypes destabilizes cell-cell adhesive junctions and affects cell migratory and invasive properties [5]. Thus, the expression of cell adhesion molecules, such as E-cadherin, EpCAM, CD44, and CD166 (ALCAM), are lost at the invasive front of CRC. Our results revealed that BMAL1-KD, in addition to E-cadherin and EpCAM, increased membranous expression of several cell surface proteins implicated in the matrix-cell adhesion [31,32] and the apico-basolateral polarity [36] of CRC cells such as integrin β1 (CD29), integrin α5 (CD49e), and CD44. Thus, the reduction in migratory and invasive capacity of BMAL1-KD CRC cells should be due not only to the higher expression of epithelial key denominators (E-cadherin and EpCAM) but also because of adhesion molecules such as CD29, CD49e, and CD44. In support of this notion, it was reported that induction of CD29/CD49e clustering reduces invasion and restores apicobasolateral polarity in invasive CRC cells [33]. Moreover, CD29 plays a crucial role in supporting tumor cell attachment to the extracellular matrix. A recent study demonstrated that miR-10a silencing suppresses migration and invasion in vitro, promoting MET-like changes in SW480 cells that lead to enhanced cell-matrix adhesion induced by CD29 [31]. Moreover, inducing MET-like changes in CRC cells, BMAL1-KD reduced both migration and invasion capacities. These results are consistent with other studies showing that knocking down CLOCK, the important partner of BMAL1, led to the impairment of CRC cell migration, while CLOCK overexpression promoted their migration [41]. In addition, our immunofluorescence results in BMAL1-KD CRC cells showed a specific honeycomb-like epithelial organization of the adhesion belts delineated by E-cadherin, β-catenin, and F-actin, a phenomenon associated with the epithelialization. It was shown that BMAL1 and CLOCK could stabilize and activate RhoA [41]. Constitutive activated RhoA could enhance cortical disassembly in
addition to induced actin polymerization in the cell interior [42]. Thus, BMAL1-KD alters the dynamics of F-actin/G-actin turnover and induced an increased cortical distribution of F-actin, which finally reduced cancer cell migration, and invasion. Globally, BMAL1-KD induced fewer phenotypic changes in the metastatic cell line SW620 cells and does not affect their migration and invasion capabilities compared to HCT116 and SW480 cells. Compared to the two primary CRC cell lines, metastatic SW620 cells showed very limited MET-like changes after BMAL1-KD. A likely explanation for this muted response is that SW620 cells present severely diminished core-clock gene oscillations as well as a low baseline expression of BMAL1 [21,43–45]. Consistent with this, our results showed also that BMAL1-KD in metastatic SW620 cells did not alter NR1D1 and CLOCK expression in contrast to our results in both primary CRC cell lines [21]. All these accumulated findings indicate that the transcription–translation feedback loops involving different circadian genes are likely to be dysfunctional in SW620 cells. Consequently, even though they were derived from the same patient, SW480 and SW620 cell lines presented only 5.5% overlap of genes with oscillating expression profiles [45]. The loss of the intact circadian clock during tumor progression probably diminished the influence of BMAL1-KD on SW620 cells compared to SW480 cells. However, only a few studies have investigated the links between the circadian timing system and the metastatic process, and there are no direct clinic proofs to show a more obvious dysfunction of transcription–translation feedback loops among circadian proteins in metastatic CRC than in primary CRC. In addition, it has been shown in vitro that a rich variety of circadian phenotypes (strong, weak, and no-oscillation of BMAL1) exist in the different CRC cell lines [43]. Moreover, it will be important to take into account the gene mutations (KRAS, BRAF, PIK3CA, CTNNB1, BRCA2, CDKN2A, P53, APC) present in the different CRC cell lines. All these distinct mutations, as well as the different molecular circadian clock statuses, provide unique pathological backgrounds in CRC cell lines, which could produce various cell fates after BMAL1-KD. In our study, HCT116, SW480, and SW620 cells were all mtKRAS cell lines, suggesting that our results may be valid in a context of mutated KRAS, which accounts for 35-45% of all metastatic CRC [46]. Further investigations using metastatic CRC cell lines with other mutations are required to evaluate if a link between the circadian clock and the meta-static process, in particular EMT, may exist.

By repressing the expression of the EMT-TF Twist as well as increasing EpCAM and E-cadherin at the plasma membrane of CRC cells, BMAL1-KD promoted MET-like changes, supporting the conclusion that the circadian transcription factor BMAL1 plays a pivotal role in CRC metastasis progression, regulating the balance between EMT and MET. In a previous study, we observed that BMAL1-KD enhanced PI3K/AKT/mTOR signaling that triggered different CRC cell fates (proliferation, apoptosis) [21]. While numerous studies have presented the activation of the PI3K/AKT axis as a key signaling pathway controlling EMT induction and tumor progression [47,48], BMAL1-KD favors, in the current study, the MET process of CRC cell lines. This discrepancy may be explained by the fact that EMT can be influenced either directly or indirectly by multiple signaling pathways, which cooperate with the PI3K/AKT axis, such as TGF-β, NF-κB, Ras, and Wnt/β-catenin. Although the activation of the PI3K/AKT axis could promote EMT, inhibition of Wnt/β-catenin signaling as well as a decrease in Twist expression would, on the other hand, engage BMAL1-KD cells in MET. Further studies are required to check whether other BMAL1-dependent signaling pathways are involved in the EMT process.

In recent years, circadian biology is becoming a critical factor for improving drug efficacy and diminishing drug toxicity in cancer treatment (chronotherapy) [49]. The third-generation platinum drug oxaliplatin (LOH) is widely used, in combination with 5-fluorouracil and leucovorin, as a first-line chemotherapeutic agent in patients with metastatic CRC. Although chemotherapy using LOH improves overall survival for advanced or metastatic CRC patients, intrinsic or acquired LOH-resistance is still a major factor in the poor prognosis associated with advanced CRC [50]. However, the link between the circadian cycle and drug resistance mechanisms has been the subject of only
addition to induced actin polymerization in the cell interior [42]. Thus, BMAL1-KD alters the dynamics of F-actin/G-actin turnover and induced an increased cortical distribution of F-actin, which finally reduced cancer cell migration, and invasion. Globally, BMAL1-KD induced fewer phenotypic changes in the metastatic cell line SW620 cells and does not affect their migration and invasion capabilities compared to HCT116 and SW480 cells. Compared to the two primary CRC cell lines, metastatic SW620 cells showed very limited MET-like changes after BMAL1-KD. A likely explanation for this muted response is that SW620 cells present severely diminished core-clock gene oscillations as well as a low baseline expression of BMAL1 [21,43–45]. Consistent with this, our results showed also that BMAL1-KD in metastatic SW620 cells did not alter NR1D1 and CLOCK expression in contrast to our results in both primary CRC cell lines [21]. All these accumulated findings indicate that the transcription–translation feedback loops involving different circadian genes are likely to be dysfunctional in SW620 cells. Consequently, even though they were derived from the same patient, SW480 and SW620 cell lines presented only 5.5% overlap of genes with oscillating expression profiles [45]. The loss of the intact circadian clock during tumor progression probably diminished the influence of BMAL1-KD on SW620 cells compared to SW480 cells. However, only a few studies have investigated the links between the circadian timing system and the metastatic process, and there are no direct clinic proofs to show a more obvious dysfunction of transcription–translation feedback loops among circadian proteins in metastatic CRC than in primary CRC. In addition, it has been shown in vitro that a rich variety of circadian phenotypes (strong, weak, and no-oscillation of BMAL1) exist in the different CRC cell lines [43]. Moreover, it will be important to take into account the gene mutations (KRAS, BRAF, PIK3CA, CTNNB1, BRCA2, CDKN2A, P53, APC) present in the different CRC cell lines. All these distinct mutations, as well as the different molecular circadian clock statuses, provide unique pathological backgrounds in CRC cell lines, which could produce various cell fates after BMAL1-KD. In our study, HCT116, SW480, and SW620 cells were all mtKRAS cell lines, suggesting that our results may be valid in a context of mutated KRAS, which accounts for 35-45% of all metastatic CRC [46]. Further investigations using metastatic CRC cell lines with other mutations are required to evaluate if a link between the circadian clock and the meta-static process, in particular EMT, may exist.

By repressing the expression of the EMT-TF Twist as well as increasing EpCAM and E-cadherin at the plasma membrane of CRC cells, BMAL1-KD promoted MET-like changes, supporting the conclusion that the circadian transcription factor BMAL1 plays a pivotal role in CRC metastasis progression, regulating the balance between EMT and MET. In a previous study, we observed that BMAL1-KD enhanced PI3K/AKT/mTOR signaling that triggered different CRC cell fates (proliferation, apoptosis) [21]. While numerous studies have presented the activation of the PI3K/AKT axis as a key signaling pathway controlling EMT induction and tumor progression [47,48], BMAL1-KD favors, in the current study, the MET process of CRC cell lines. This discrepancy may be explained by the fact that EMT can be influenced either directly or indirectly by multiple signaling pathways, which cooperate with the PI3K/AKT axis, such as TGF-β, NF-κB, Ras, and Wnt/β-catenin. Although the activation of the PI3K/AKT axis could promote EMT, inhibition of Wnt/β-catenin signaling as well as a decrease in Twist expression would, on the other hand, engage BMAL1-KD cells in MET. Further studies are required to check whether other BMAL1-dependent signaling pathways are involved in the EMT process.

In recent years, circadian biology is becoming a critical factor for improving drug efficacy and diminishing drug toxicity in cancer treatment (chronotherapy) [49]. The third-generation platinum drug oxaliplatin (LOH) is widely used, in combination with 5-fluorouracil and leucovorin, as a first-line chemotherapeutic agent in patients with metastatic CRC. Although chemotherapy using LOH improves overall survival for advanced or metastatic CRC patients, intrinsic or acquired LOH-resistance is still a major factor in the poor prognosis associated with advanced CRC [50]. However, the link between the circadian cycle and drug resistance mechanisms has been the subject of only
few analyses [16,51–53]. Among the multiple mechanisms, EMT is closely associated with LOH-resistance in CRC [54]. Indeed, CRC cells undergoing EMT induced by TGF-β treatment [36] or EMT-TF Twist overexpression [55] become resistant to LOH chemotherapy. Acquired LOH resistance can also be induced by downregulation of the EMT-related miR-200c and miR-141 [56]. On the other side, selected LOH chemoresistant CRC cells acquire an EMT phenotype [34]. However, few reports have so far linked the circadian core genes to the EMT process. For instance, PER2 may suppress EMT and tumor malignancy in breast cancer cell lines [57] and in glioma cells [58,59]. In our work, consistent with the MET induction, BMAL1-KD increased the LOH sensitivity of CRC cell lines, demonstrating that BMAL1 is implicated in EMT-mediated chemoresistance in CRC. Curiously, it was reported that BMAL1 overexpression also sensitized CRC cell lines to LOH [53]. In fact, we suppose that overexpression or knockdown of BMAL1 could disrupt the circadian oscillation pattern of BMAL1 and both influence the expression of core circadian genes and BMAL1-dependent target genes [21,60], leading, ultimately, all to an increased sensitivity to LOH treatment. All these data suggested that an appropriate BMAL1 expression level is probably a key element in affecting CRC chemoresistance to LOH.

In summary, our study showed that the circadian transcription factor BMAL1 plays
crucial functions in CRC metastasis progression, promoting EMT, which is a key step in the
tumors' invasiveness and chemoresistance. Thus, therapeutic targeting of BMAL1 could be
considered as a potential strategy to enhance the efficacy of chemotherapy and improve
outcomes for patients with metastatic CRC.

4. Materials and Methods

4.1. RNAseq Analysis

The TCGA consortium RNAsequencing data matrix of Z-scores from the CRC co-
hort [61] were downloaded on the cBioPortal website [62], with access that provided
associated and clinical data for each patient of the cohort. CRC patients were stratified on
Z-scores of BMAL1 found in RNAseq as extremes of Gaussian distribution with 6 patients
having Z-scores less than minus 1.96 (BMAL1_low) and with 5 patients having Z-scores
over plus 1.96 (BMAL1_high). Gene set enrichment analysis was performed with Gene
Set Enrichment Analysis (GSEA) software version 4.0.3 [63] on MsigDB version 6.2 [64].
Bioinformatics analyses were performed in R software environment version 3.5.3: boxplot
with ggplot2 graph definition [65], principal component analysis with FactoMineR R pack-
age [66] and expression heatmap with pheatmap R package. Statistical significance was
attested with two-sided Student's t-test between groups of patients with an alpha error
below 0.05.

4.2. *Cell Culture*

Human colon cancer cell lines HCT116, SW480, and SW620 were grown in Dulbecco's modified Eagle's medium (DMEM) supplemented with GlutaMAX (Gibco-BRL, Life Technology, Saint Aubin, France) and 10% fetal bovine serum (FBS, Hyclone, Logan, UT, USA).

4.3. *The shRNA Cloning in Lentiviral Vector*

One short hairpin RNAs (shRNAs) against human BMAL1 and a control scrambled sequence (Scr) were cloned separately into a lentiviral vector (pLKO-shBMAL1-GFP-puro and pLKO-Scr-GFP-puro), as described previously [67].

4.4. Lentivirus Production and Cell Transduction

Lentivirus particles were prepared and used to transduce human colon cancer cell lines HCT116, SW480, and SW620 as previously described [21,67].

4.5. Cytoplasmic and Nuclear Extracts Preparation

CRC cell cultures were scrapped in RIPA buffer containing protease inhibitors. Then, cytoplasmic and nuclear fractions of different cell lines were extracted by using NE-PER™
few analyses [16,51–53]. Among the multiple mechanisms, EMT is closely associated with LOH-resistance in CRC [54]. Indeed, CRC cells undergoing EMT induced by TGF-β treatment [36] or EMT-TF Twist overexpression [55] become resistant to LOH chemotherapy. Acquired LOH resistance can also be induced by downregulation of the EMT-related miR-200c and miR-141 [56]. On the other side, selected LOH chemoresistant CRC cells acquire an EMT phenotype [34]. However, few reports have so far linked the circadian core genes to the EMT process. For instance, PER2 may suppress EMT and tumor malignancy in breast cancer cell lines [57] and in glioma cells [58,59]. In our work, consistent with the MET induction, BMAL1-KD increased the LOH sensitivity of CRC cell lines, demonstrating that BMAL1 is implicated in EMT-mediated chemoresistance in CRC. Curiously, it was reported that BMAL1 overexpression also sensitized CRC cell lines to LOH [53]. In fact, we suppose that overexpression or knockdown of BMAL1 could disrupt the circadian oscillation pattern of BMAL1 and both influence the expression of core circadian genes and BMAL1-dependent target genes [21,60], leading, ultimately, all to an increased sensitivity to LOH treatment. All these data suggested that an appropriate BMAL1 expression level is probably a key element in affecting CRC chemoresistance to LOH.

In summary, our study showed that the circadian transcription factor BMAL1 plays crucial functions in CRC metastasis progression, promoting EMT, which is a key step in the tumors' invasiveness and chemoresistance. Thus, therapeutic targeting of BMAL1 could be considered as a potential strategy to enhance the efficacy of chemotherapy and improve outcomes for patients with metastatic CRC.

# 4. Materials and Methods

## 4.1. RNAseq Analysis

The TCGA consortium RNAsequencing data matrix of Z-scores from the CRC cohort [61] were downloaded on the cBioPortal website [62], with access that provided associated and clinical data for each patient of the cohort. CRC patients were stratified on Z-scores of BMAL1 found in RNAseq as extremes of Gaussian distribution with 6 patients having Z-scores less than minus 1.96 (BMAL1_low) and with 5 patients having Z-scores over plus 1.96 (BMAL1_high). Gene set enrichment analysis was performed with Gene Set Enrichment Analysis (GSEA) software version 4.0.3 [63] on MsigDB version 6.2 [64]. Bioinformatics analyses were performed in R software environment version 3.5.3: boxplot with ggplot2 graph definition [65], principal component analysis with FactoMineR package [66] and expression heatmap with pheatmap R package. Statistical significance was attested with two-sided Student's t-test between groups of patients with an alpha error below 0.05.

## 4.2. Cell Culture

Human colon cancer cell lines HCT116, SW480, and SW620 were grown in Dulbecco's modified Eagle's medium (DMEM) supplemented with GlutaMAX (Gibco-BRL, Life Technology, Saint Aubin, France) and 10% fetal bovine serum (FBS, Hyclone, Logan, UT, USA).

## 4.3. The shRNA Cloning in Lentiviral Vector

One short hairpin RNAs (shRNAs) against human BMAL1 and a control scrambled sequence (Scr) were cloned separately into a lentiviral vector (pLKO-shBMAL1-GFP-puro and pLKO-Scr-GFP-puro), as described previously [67].

## 4.4. Lentivirus Production and Cell Transduction

Lentivirus particles were prepared and used to transduce human colon cancer cell lines HCT116, SW480, and SW620 as previously described [21,67].

## 4.5. Cytoplasmic and Nuclear Extracts Preparation

CRC cell cultures were scrapped in RIPA buffer containing protease inhibitors. Then, cytoplasmic and nuclear fractions of different cell lines were extracted by using NE-PER™
Nuclear and Cytoplasmic Extraction Reagents (Thermo Fisher Scientific, Illkirch, France) according to the manufacturer's manual. The different extracts were stored at -80° C until Western blot analysis.

## 4.6. Western-Blot Analysis

Western blot analysis was performed with sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE) as previously described [21]. Antibodies were used against BMAL1 (Abcam, Cambridge, UK, ab3350), E-cadherin (Cell signaling Technology, Inc., Boston, MA, USA, 14472), vimentin (Cell signaling Technology, 73260), β-catenin (Cell signaling Technology, 8480). HSC-70 (Stressgen, CA, USA, SPA-815), β-actin (Cell signaling Technology, 47778) and Lamin β1 (Abcam, ab16048) were used as control for protein loading. Image J 1.4.3 was used to quantification the western results. The *t*-test was used for statistical analysis.

## 4.7. Quantitative Real-Time PCR

Quantitative real-time PCR was performed as previous described by using LightCycler 480 SYBR Green I master kit (Roche Applied Science, Basel, Switzerland, 04707516001) [68]. Primers are listed in Table S1. Hybridization temperature for all primers was 60 °C. The relative quantification of target RNA by using 36B4 as a reference was computed with the Relquant software (Roche, Bâle, Switzerland) with the "delta delta Ct" (ΔΔCt) method. The *t*-test was used for statistical analysis.

## 4.8. Flow Cytometry Analysis

Cells were detached with Accutase (Sigma-Aldrich, St. Quentin Fallavier, France, A6964) and cell surface expression of molecules was analyzed by flow cytometry. Briefly, cells were permeabilized or not with BD Cytofix/Cytoperm reagent (BD Pharmingen, Le Pont De Claix, France), and 105 cells were suspended in DMEM medium supplemented with 1% FCS. Indirect stainings were performed for cell surface expression of EpCAM (Bio-Techne Ltd., Lille, France, AF640), E-cadherin (Bio-Techne Ltd., AF648), N-cadherin (Abcam, ab1221), and intracellular expression of Cytokeratin-20 (CK-20, Santa Cruz, Dallas, TX, USA, sc-17113) and vimentin (Santa Cruz, sc-7557), using species-specific secondary antibodies conjugated to phycoerythrin. Phycoerythrin (PE)-conjugated antibodies against CD29 (Immunotools, Friesoythe, Germany, 21270294), CD49e (Immunotools, 21336491), and CD44 (Immunotools, 21270444) were used to perform direct stainings of cell surface expression markers. After 3 washes, 10,000 cells were analyzed on a LSR FortessaTM flow cytometer (BecTon Dickinson, NJ, USA), and the data were analyzed using FlowJo software (Tree Star Inc., Ashland, OR, USA). The experiment was repeated at least three times. The *t*-test was used for statistical analysis.

## 4.9. Immunofluorescence and Confocal Microscopy

Immunofluorescence stainings were performed as previously described [21] using anti-β-catenin (Cell signaling Technology, 8480) and anti-E-cadherin (Cell signaling Technology, 14472) antibodies and Hoechst 33342 (Cell signaling Technology, 4082) to stain for nuclear. The confocal images were captured by a confocal LEICA SP5-AOBS microscope with a 63X/1.4 NA oil-immersion objective.

## 4.10. Migration and Invasion Assays

Migration was assessed using a wound healing assay (scratch assay). Briefly, CRC cells were cultured in 6-well plates for 24 h to achieve 100% confluence followed by starvation in serum-free DMEM medium containing 10 µg/mL mitomycin C (Thermo Fisher Scientific, BP25312) for 2 h to completely inhibit cell proliferation. After generating a wound in the cell monolayer with a sterile 200 µL pipette tip, detached cells and cellular debris were removed three times by washing with culture medium. The cells were fed with fresh growth medium and the width of the gap between the invasion fronts was measured using
Nuclear and Cytoplasmic Extraction Reagents (Thermo Fisher Scientific, Illkirch, France) according to the manufacturer's manual. The different extracts were stored at -80° C until Western blot analysis.

### 4.6. Western-Blot Analysis

Western blot analysis was performed with sodium dodecyl sulfate polyacrylamide gel electrophoresis (SDS-PAGE) as previously described [21]. Antibodies were used against BMAL1 (Abcam, Cambridge, UK, ab3350), E-cadherin (Cell signaling Technology, Inc., Boston, MA, USA, 14472), vimentin (Cell signaling Technology, 73260), β-catenin (Cell signaling Technology, 8480). HSC-70 (Stressgen, CA, USA, SPA-815), β-actin (Cell signaling Technology, 47778) and Lamin β1 (Abcam, ab16048) were used as control for protein loading. Image J 1.4.3 was used to quantification the western results. The *t*-test was used for statistical analysis.

### 4.7. Quantitative Real-Time PCR

Quantitative real-time PCR was performed as previous described by using LightCycler 480 SYBR Green I master kit (Roche Applied Science, Basel, Switzerland, 04707516001) [68]. Primers are listed in Table S1. Hybridization temperature for all primers was 60 °C. The relative quantification of target RNA by using 36B4 as a reference was computed with the Relquant software (Roche, Bâle, Switzerland) with the "delta delta Ct" (ΔΔCt) method. The *t*-test was used for statistical analysis.

### 4.8. Flow Cytometry Analysis

Cells were detached with Accutase (Sigma-Aldrich, St. Quentin Fallavier, France, A6964) and cell surface expression of molecules was analyzed by flow cytometry. Briefly, cells were permeabilized or not with BD Cytofix/Cytoperm reagent (BD Pharmingen, Le Pont De Claix, France), and 105 cells were suspended in DMEM medium supplemented with 1% FCS. Indirect stainings were performed for cell surface expression of EpCAM (Bio-Techne Ltd., Lille, France, AF640), E-cadherin (Bio-Techne Ltd., AF648), N-cadherin (Abcam, ab1221), and intracellular expression of Cytokeratin-20 (CK-20, Santa Cruz, Dallas, TX, USA, sc-17113) and vimentin (Santa Cruz, sc-7557), using species-specific secondary antibodies conjugated to phycoerythrin. Phycoerythrin (PE)-conjugated antibodies against CD29 (Immunotools, Friesoythe, Germany, 21270294), CD49e (Immunotools, 21336491), and CD44 (Immunotools, 21270444) were used to perform direct stainings of cell surface expression markers. After 3 washes, 10,000 cells were analyzed on a LSR FortessaTM flow cytometer (Becton Dickinson, NJ, USA), and the data were analyzed using FlowJo software (Tree Star Inc., Ashland, OR, USA). The experiment was repeated at least three times. The *t*-test was used for statistical analysis.

### 4.9. Immunofluorescence and Confocal Microscopy

Immunofluorescence stainings were performed as previously described [21] using anti-β-catenin (Cell signaling Technology, 8480) and anti-E-cadherin (Cell signaling Technology, 14472) antibodies and Hoechst 33342 (Cell signaling Technology, 4082) to stain for nuclear. The confocal images were captured by a confocal LEICA SP5-AOBS microscope with a 63X/1.4 NA oil-immersion objective.

### 4.10. Migration and Invasion Assays

Migration was assessed using a wound healing assay (scratch assay). Briefly, CRC cells were cultured in 6-well plates for 24 h to achieve 100% confluence followed by starvation in serum-free DMEM medium containing 10 µg/mL mitomycin C (Thermo Fisher Scientific, BP25312) for 2 h to completely inhibit cell proliferation. After generating a wound in the cell monolayer with a sterile 200 µL pipette tip, detached cells and cellular debris were removed three times by washing with culture medium. The cells were fed with fresh growth medium and the width of the gap between the invasion fronts was measured using
Image-J software (NIH, Bethesda, MD, USA) at time points from 0 to 72 h after scratching to calculate the rate of wound closure.

Invasion assays were performed in 24-well Transwell chambers with 8 µm pore polycarbonate filter inserts (Corning, Becton Dickinson Labware, Franklin Lakes, NJ, USA). CRC cells were serum starved and exposed to mitomycin C before seeding (5 × 10⁴) on Matrigel-coated inserts (BD Biosciences, Bedford, MA, USA) in 100 µL of serum-free medium inserts. The lower chambers were filled with 0.5 mL of complete media with 10% FBS. After 24 h and/or 48 h, assays were terminated by scraping the top of the membrane to remove non-migratory cells. The cells on the lower surface of the insert were fixed in 4% paraformaldehyde and stained with 4',6-Diamidino-2-Phenylindole, Dihydrochloride (DAPI, Thermo Fisher Scientific, D1306) and counted. Quantification of cells was carried out by counting 10 random fields using a 10X objective. Assays were performed in triplicates.

### 4.11. Drug Sensitivity Assay

The viability of CRC cells to Oxaloplatin (LOH) was measured using the MTT assay. Cells were seeded in quadruplicate at 3 × 10³ per well in 96-well plates and incubated overnight at 37 °C in a humidified environment containing 5% CO₂. Oxaliplatin was added to the cell culture medium at various concentrations (ranging from 12.5 to 75 µM). Untreated cells were used as control. After 48 h of treatment, 20 µL of 5 mg/mL MTT (HiMedia, Levallois-Perret, France, RM1131-1G) solution was added to each well and incubated at 37 °C for 3 h. Then, the supernatant was removed and 100 µL of 1:1 solution of DMSO and methanol were added to dissolve the formazan precipitate. Absorbance of the solution was read using a microplate reader (Bio-Rad Laboratories, Hercules, CA, USA) at 540 nm, with a background correction at 655 nm. The results were analyzed using GraphPad PRISM software (GraphPad, San Diego, CA, USA) and were plotted as the mean ± SEM of the absorbance for each dose of drug tested from three independent experiments.

### 4.12. Statistical Analysis

The data were analyzed by GraphPad Prism 8.0.1. Statistically significant differences were determined by *t*-test or two-way ANOVA, and *p* < 0.05 was considered to be a statistically significant result. Error bars represent standard error of the mean (±SEM) of independent experiments.

**Supplementary Materials:** The following are available online at https://www.mdpi.com/article/10.3390/ijms22105247/s1, Figure S1: Cell surface expression of N-cadherin on the three BMAL1-KD cell lines, Table S1: Sequences of primers used for real-time quantitative PCR amplification.

**Author Contributions:** Y.C. and J.G.-M. designed the research. Y.Z., A.D., C.D., L.E.B.d.S., Y.X., É.H., Y.C. and J.G.-M. performed the research. Y.Z., A.D., C.D., L.E.B.d.S., Y.X., É.H., Y.C. and J.G.-M. analyzed the data. Y.C. and J.G.-M. wrote the original draft paper. H.A., A.F., A.B., Y.C. and J.G.-M. wrote—reviewed and edited the paper. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research received no external funding.

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** The data presented in this study are available upon request from the corresponding author.

**Acknowledgments:** This work was supported by INSERM and the University Paris-Saclay. Y.Z. and Y.X. were supported by the China Scholarship Council. A.B. was funded by the Plan Cancer of the French National Cancer Institute (INCa) through the ATIP-Avenir program.

**Conflicts of Interest:** The authors declare that they have no conflict of interest.
Image-J software (NIH, Bethesda, MD, USA) at time points from 0 to 72 h after scratching to calculate the rate of wound closure.

Invasion assays were performed in 24-well Transwell chambers with 8 µm pore polycarbonate filter inserts (Corning, Becton Dickinson Labware, Franklin Lakes, NJ, USA). CRC cells were serum starved and exposed to mitomycin C before seeding (5 × 10⁴) on Matrigel-coated inserts (BD Biosciences, Bedford, MA, USA) in 100 µL of serum-free medium inserts. The lower chambers were filled with 0.5 mL of complete media with 10% FBS. After 24 h and/or 48 h, assays were terminated by scraping the top of the membrane to remove non-migratory cells. The cells on the lower surface of the insert were fixed in 4% paraformaldehyde and stained with 4',6-Diamidino-2-Phenylindole, Dihydrochloride (DAPI, Thermo Fisher Scientific, D1306) and counted. Quantification of cells was carried out by counting 10 random fields using a 10X objective. Assays were performed in triplicates.

### 4.11. Drug Sensitivity Assay

The viability of CRC cells to Oxaloplatin (LOH) was measured using the MTT assay. Cells were seeded in quadruplicate at 3 × 10³ per well in 96-well plates and incubated overnight at 37 °C in a humidified environment containing 5% CO₂. Oxaliplatin was added to the cell culture medium at various concentrations (ranging from 12.5 to 75 µM). Untreated cells were used as control. After 48 h of treatment, 20 µL of 5 mg/mL MTT (HiMedia, Levallois-Perret, France, RM1131-1G) solution was added to each well and incubated at 37 °C for 3 h. Then, the supernatant was removed and 100 µL of 1:1 solution of DMSO and methanol were added to dissolve the formazan precipitate. Absorbance of the solution was read using a microplate reader (Bio-Rad Laboratories, Hercules, CA, USA) at 540 nm, with a background correction at 655 nm. The results were analyzed using GraphPad PRISM software (GraphPad, San Diego, CA, USA) and were plotted as the mean ± SEM of the absorbance for each dose of drug tested from three independent experiments.

### 4.12. Statistical Analysis

The data were analyzed by GraphPad Prism 8.0.1. Statistically significant differences were determined by *t*-test or two-way ANOVA, and *p* < 0.05 was considered to be a statistically significant result. Error bars represent standard error of the mean (±SEM) of independent experiments.

**Supplementary Materials:** The following are available online at https://www.mdpi.com/article/10.3390/ijms22105247/s1, Figure S1: Cell surface expression of N-cadherin on the three BMAL1-KD cell lines, Table S1: Sequences of primers used for real-time quantitative PCR amplification.

**Author Contributions:** Y.C. and J.G.-M. designed the research. Y.Z., A.D., C.D., L.E.B.d.S., Y.X., É.H., Y.C. and J.G.-M. performed the research. Y.Z., A.D., C.D., L.E.B.d.S., Y.X., É.H., Y.C. and J.G.-M. analyzed the data. Y.C. and J.G.-M. wrote the original draft paper. H.A., A.F., A.B., Y.C. and J.G.-M. wrote—reviewed and edited the paper. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research received no external funding.

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** The data presented in this study are available upon request from the corresponding author.

**Acknowledgments:** This work was supported by INSERM and the University Paris-Saclay. Y.Z. and Y.X. were supported by the China Scholarship Council. A.B. was funded by the Plan Cancer of the French National Cancer Institute (INCa) through the ATIP-Avenir program.

**Conflicts of Interest:** The authors declare that they have no conflict of interest.
References

1. Siegel, R.L.; Medhanie, G.A.; Fedewa, S.A.; Jemal, A. State Variation in Early-Onset Colorectal Cancer in the United States, 1995–2015. *J. Natl. Cancer Inst.* **2019**, *111*, 1104–1106. [CrossRef]

2. Zhang, Y.; Ma, J.; Zhang, S.; Deng, G.; Wu, X.; He, J.; Pei, H.; Shen, H.; Zeng, S. A prognostic analysis of 895 cases of stage III colon cancer in different colon subsites. *Int. J. Color. Dis.* **2015**, *30*, 1173–1183. [CrossRef] [PubMed]

3. Thiery, J.P.; Acloque, H.; Huang, R.Y.J.; Nieto, M.A. Epithelial-Mesenchymal Transitions in Development and Disease. *Cell* **2009**, *139*, 871–890. [CrossRef]

4. Zhang, J.; Tian, X.-J.; Xing, J. Signal Transduction Pathways of EMT Induced by TGF-β, SHH, and WNT and Their Crosstalks. *J. Clin. Med.* **2016**, *5*, 41. [CrossRef] [PubMed]

5. Jeanes, A.; Gottardi, C.J.; Yap, A.S. Cadherins and cancer: How does cadherin dysfunction promote tumor progression? *Oncogene* **2008**, *27*, 6920–6929. [CrossRef] [PubMed]

6. Takahashi, J.S. Transcriptional architecture of the mammalian circadian clock. *Nat. Rev. Genet.* **2017**, *18*, 164–179. [CrossRef]

7. Fu, L.; Kettner, N.M. The Circadian Clock in Cancer Development and Therapy. *Prog. Mol. Biol. Transl. Sci.* **2013**, *119*, 221–282. [CrossRef] [PubMed]

8. Dallmann, R.; Okyar, A.; Lévi, F. Dosing-Time Makes the Poison: Circadian Regulation and Pharmacotherapy. *Trends Mol. Med.* **2016**, *22*, 430–445. [CrossRef] [PubMed]

9. Lévi, F.; Okyar, A.; Dulong, S.; Innominato, P.F.; Clairambault, J. Circadian Timing in Cancer Treatments. *Annu. Rev. Pharmacol. Toxicol.* **2010**, *50*, 377–421. [CrossRef]

10. Sulli, G.; Lam, M.T.Y.; Panda, S. Interplay between Circadian Clock and Cancer: New Frontiers for Cancer Treatment. *Trends Cancer* **2019**, *5*, 475–494. [CrossRef]

11. Straif, K.; Baan, R.; Grosse, Y.; Secretan, B.; El Ghissassi, F.; Bouvard, V.; Altieri, A.; Benbrahim-Tallaa, L.; Cogliano, V. Carcinogenicity of shift-work, painting, and fire-fighting. *Lancet Oncol.* **2007**, *8*, 1065–1066. [CrossRef]

12. Giacchetti, S. Chronotherapy of colorectal cancer. *Chronobiol. Int.* **2002**, *19*, 207–219. [CrossRef] [PubMed]

13. Focan, C.; Kreutz, F.; Graas, M.-P.; Longrée, L.; Focan-Henrard, D.; Demolin, G.; Moeneclaey, N. Phase I–II study to assess the feasibility and activity of the triple combination of 5-fluorouracil/folinic acid, carboplatin and irinotecan (CPT-11) administered by chronomodulated infusion for the treatment of advanced colorectal cancer. Final report of the BE-1603 study. *Pathol. Biol.* **2013**, *61*, e27–e31. [CrossRef] [PubMed]

14. Bunger, M.K.; Wilsbacher, L.D.; Moran, S.M.; Clendenin, C.; Radcliffe, L.A.; HogenEsch, J.B.; Simon, M.C.; Takahashi, J.S.; Bradfield, C.A. Mop3 is an essential component of the master circadian pacemaker in mammals. *Cell* **2000**, *103*, 1009–1017. [CrossRef]

15. Jiang, W.; Zhao, S.; Jiang, X.; Zhang, E.; Hu, G.; Hu, B.; Zheng, P.; Xiao, J.; Lu, Z.; Lu, Y.; et al. The circadian clock gene Bmal1 acts as a potential anti-oncogene in pancreatic cancer by activating the p53 tumor suppressor pathway. *Cancer Lett.* **2016**, *371*, 314–325. [CrossRef]

16. Tang, Q.; Cheng, B.; Xie, M.; Chen, Y.; Zhao, J.; Zhou, X.; Chen, L. Circadian Clock Gene Bmal1 Inhibits Tumorigenesis and Increases Paclitaxel Sensitivity in Tongue Squamous Cell Carcinoma. *Cancer Res.* **2017**, *77*, 532–544. [CrossRef]

17. Fu, L.; Lee, C.C. The circadian clock: Pacemaker and tumour suppressor. *Nat. Rev. Cancer* **2003**, *3*, 350–361. [CrossRef]

18. Ramos, C.A.; Ouyang, C.; Qi, Y.; Chung, Y.; Cheng, C.-T.; LaBarge, M.A.; Seewaldt, V.L.; Ann, D.K. A Non-canonical Function of BMAL1 Metabolically Limits Obesity-Promoted Triple-Negative Breast Cancer. *iScience* **2020**, *23*, 100839. [CrossRef]

19. Oshima, T.; Takenoshita, S.; Akaike, M.; Kunisaki, C.; Fujii, S.; Nozaki, A.; Numata, K.; Shiozawa, M.; Rino, Y.; Tanaka, K.; et al. Expression of circadian genes correlates with liver metastasis and outcomes in colorectal cancer. *Oncol. Rep.* **2011**, *25*, 1439–1446. [CrossRef] [PubMed]

20. Burgermeister, E.; Battaglin, F.; Eladly, F.; Wu, W.; Herweck, F.; Schulte, N.; Betge, J.; Härtel, N.; Kather, J.N.; Weis, C.-A.; et al. Aryl hydrocarbon receptor nuclear translocator-like (ARNTL/BMAL1) is associated with bevacizumab resistance in colorectal cancer via regulation of vascular endothelial growth factor A. *EBioMedicine* **2019**, *45*, 139–154. [CrossRef] [PubMed]

21. Zhang, Y.; Devocelle, A.; Souza, L.; Foudi, A.; Bento, S.T.; Desterke, C.; Sherrard, R.; Ballesta, A.; Adam, R.; Giron-Michel, J.; et al. BMAL1 knockdown triggers different colon carcinoma cell fates by altering the delicate equilibrium between AKT/mTOR and P53/P21 pathways. *Aging* **2020**, *10*, 8067–8083. [CrossRef]

22. Frixen, U.H.; Behrens, J.; Sachs, M.; Eberle, G.; Voss, B.; Warda, A.; Löchner, D.; Birchmeier, W. E-cadherin-mediated cell-cell adhesion prevents invasiveness of human carcinoma cells. *J. Cell Biol.* **1991**, *113*, 173–185. [CrossRef]

23. Kinsella, A.R.; Lepts, G.C.; Hill, C.L.; Jones, M. Reduced E-cadherin expression correlates with increased invasiveness in colorectal carcinoma cell lines. *Clin. Exp. Metastasis* **1994**, *12*, 335–342. [CrossRef]

24. Huang, J.; Ye, J.; Wu, D.; Shen, J.; Wu, P.; Ni, C.; Chen, J.; Zhao, J.; Zhang, T.; Wang, X. Enrichment of colorectal cancer stem cells through epithelial-mesenchymal transition via CDH1 knockdown. *Mol. Med. Rep.* **2012**, *6*, 507–512. [CrossRef]

25. Tian, X.; Liu, Z.; Niu, B.; Zhang, J.; Tan, T.K.; Lee, S.R.; Zhao, Y.; Harris, D.C.H.; Zheng, G. E-Cadherin/β-Catenin Complex and the Epithelial Barrier. *J. Biomed. Biotechnol.* **2011**, *2011*, 567305. [CrossRef]

26. Shang, S.; Hua, F.; Hu, Z.-W. The regulation of β-catenin activity and function in cancer: Therapeutic opportunities. *Oncotarget* **2017**, *8*, 33972–33989. [CrossRef] [PubMed]

27. Gilles, C.; Polette, M.; Mestdagt, M.; Nawrocki-Raby, B.; Ruggeri, P.; Birembaut, P.; Foidart, J.-M. Transactivation of vimentin by beta-catenin in human breast cancer cells. *Cancer Res.* **2003**, *63*, 2658–2664. [CrossRef] [PubMed]
## References

1. Siegel, R.L.; Medhanie, G.A.; Fedewa, S.A.; Jemal, A. State Variation in Early-Onset Colorectal Cancer in the United States, 1995–2015. *J. Natl. Cancer Inst.* **2019**, *111*, 1104–1106. [CrossRef]

2. Zhang, Y.; Ma, J.; Zhang, S.; Deng, G.; Wu, X.; He, J.; Pei, H.; Shen, H.; Zeng, S. A prognostic analysis of 895 cases of stage III colon cancer in different colon subsites. *Int. J. Color. Dis.* **2015**, *30*, 1173–1183. [CrossRef] [PubMed]

3. Thiery, J.P.; Acloque, H.; Huang, R.Y.J.; Nieto, M.A. Epithelial-Mesenchymal Transitions in Development and Disease. *Cell* **2009**, *139*, 871–890. [CrossRef]

4. Zhang, J.; Tian, X.-J.; Xing, J. Signal Transduction Pathways of EMT Induced by TGF-β, SHH, and WNT and Their Crosstalks. *J. Clin. Med.* **2016**, *5*, 41. [CrossRef] [PubMed]

5. Jeanes, A.; Gottardi, C.J.; Yap, A.S. Cadherins and cancer: How does cadherin dysfunction promote tumor progression? *Oncogene* **2008**, *27*, 6920–6929. [CrossRef] [PubMed]

6. Takahashi, J.S. Transcriptional architecture of the mammalian circadian clock. *Nat. Rev. Genet.* **2017**, *18*, 164–179. [CrossRef]

7. Fu, L.; Kettner, N.M. The Circadian Clock in Cancer Development and Therapy. *Prog. Mol. Biol. Transl. Sci.* **2013**, *119*, 221–282. [CrossRef] [PubMed]

8. Dallmann, R.; Okyar, A.; Lévi, F. Dosing-Time Makes the Poison: Circadian Regulation and Pharmacotherapy. *Trends Mol. Med.* **2016**, *22*, 430–445. [CrossRef] [PubMed]

9. Lévi, F.; Okyar, A.; Dulong, S.; Innominato, P.F.; Clairambault, J. Circadian Timing in Cancer Treatments. *Annu. Rev. Pharmacol. Toxicol.* **2010**, *50*, 377–421. [CrossRef]

10. Sulli, G.; Lam, M.T.Y.; Panda, S. Interplay between Circadian Clock and Cancer: New Frontiers for Cancer Treatment. *Trends Cancer* **2019**, *5*, 475–494. [CrossRef]

11. Straif, K.; Baan, R.; Grosse, Y.; Secretan, B.; El Ghissassi, F.; Bouvard, V.; Altieri, A.; Benbrahim-Tallaa, L.; Cogliano, V. Carcinogenicity of shift-work, painting, and fire-fighting. *Lancet Oncol.* **2007**, *8*, 1065–1066. [CrossRef]

12. Giacchetti, S. Chronotherapy of colorectal cancer. *Chronobiol. Int.* **2002**, *19*, 207–219. [CrossRef] [PubMed]

13. Focan, C.; Kreutz, F.; Graas, M.-P.; Longrée, L.; Focan-Henrard, D.; Demolin, G.; Moeneclaey, N. Phase I–II study to assess the feasibility and activity of the triple combination of 5-fluorouracil/folinic acid, carboplatin and irinotecan (CPT-11) administered by chronomodulated infusion for the treatment of advanced colorectal cancer. Final report of the BE-1603 study. *Pathol. Biol.* **2013**, *61*, e27–e31. [CrossRef] [PubMed]

14. Bunger, M.K.; Wilsbacher, L.D.; Moran, S.M.; Clendenin, C.; Radcliffe, L.A.; HogenEsch, J.B.; Simon, M.C.; Takahashi, J.S.; Bradfield, C.A. Mop3 is an essential component of the master circadian pacemaker in mammals. *Cell* **2000**, *103*, 1009–1017. [CrossRef]

15. Jiang, W.; Zhao, S.; Jiang, X.; Zhang, E.; Hu, G.; Hu, B.; Zheng, P.; Xiao, J.; Lu, Z.; Lu, Y.; et al. The circadian clock gene Bmal1 acts as a potential anti-oncogene in pancreatic cancer by activating the p53 tumor suppressor pathway. *Cancer Lett.* **2016**, *371*, 314–325. [CrossRef]

16. Tang, Q.; Cheng, B.; Xie, M.; Chen, Y.; Zhao, J.; Zhou, X.; Chen, L. Circadian Clock Gene Bmal1 Inhibits Tumorigenesis and Increases Paclitaxel Sensitivity in Tongue Squamous Cell Carcinoma. *Cancer Res.* **2017**, *77*, 532–544. [CrossRef]

17. Fu, L.; Lee, C.C. The circadian clock: Pacemaker and tumour suppressor. *Nat. Rev. Cancer* **2003**, *3*, 350–361. [CrossRef]

18. Ramos, C.A.; Ouyang, C.; Qi, Y.; Chung, Y.; Cheng, C.-T.; LaBarge, M.A.; Seewaldt, V.L.; Ann, D.K. A Non-canonical Function of BMAL1 Metabolically Limits Obesity-Promoted Triple-Negative Breast Cancer. *iScience* **2020**, *23*, 100839. [CrossRef]

19. Oshima, T.; Takenoshita, S.; Akaike, M.; Kunisaki, C.; Fujii, S.; Nozaki, A.; Numata, K.; Shiozawa, M.; Rino, Y.; Tanaka, K.; et al. Expression of circadian genes correlates with liver metastasis and outcomes in colorectal cancer. *Oncol. Rep.* **2011**, *25*, 1439–1446. [CrossRef] [PubMed]

20. Burgermeister, E.; Battaglin, F.; Eladly, F.; Wu, W.; Herweck, F.; Schulte, N.; Betge, J.; Härtel, N.; Kather, J.N.; Weis, C.-A.; et al. Aryl hydrocarbon receptor nuclear translocator-like (ARNTL/BMAL1) is associated with bevacizumab resistance in colorectal cancer via regulation of vascular endothelial growth factor A. *EBioMedicine* **2019**, *45*, 139–154. [CrossRef] [PubMed]

21. Zhang, Y.; Devocelle, A.; Souza, L.; Foudi, A.; Bento, S.T.; Desterke, C.; Sherrard, R.; Ballesta, A.; Adam, R.; Giron-Michel, J.; et al. BMAL1 knockdown triggers different colon carcinoma cell fates by altering the delicate equilibrium between AKT/mTOR and P53/P21 pathways. *Aging* **2020**, *10*, 8067–8083. [CrossRef]

22. Frixen, U.H.; Behrens, J.; Sachs, M.; Eberle, G.; Voss, B.; Warda, A.; Löchner, D.; Birchmeier, W. E-cadherin-mediated cell-cell adhesion prevents invasiveness of human carcinoma cells. *J. Cell Biol.* **1991**, *113*, 173–185. [CrossRef]

23. Kinsella, A.R.; Lepts, G.C.; Hill, C.L.; Jones, M. Reduced E-cadherin expression correlates with increased invasiveness in colorectal carcinoma cell lines. *Clin. Exp. Metastasis* **1994**, *12*, 335–342. [CrossRef]

24. Huang, J.; Ye, J.; Wu, D.; Shen, J.; Wu, P.; Ni, C.; Chen, J.; Zhao, J.; Zhang, T.; Wang, X. Enrichment of colorectal cancer stem cells through epithelial-mesenchymal transition via CDH1 knockdown. *Mol. Med. Rep.* **2012**, *6*, 507–512. [CrossRef]

25. Tian, X.; Liu, Z.; Niu, B.; Zhang, J.; Tan, T.K.; Lee, S.R.; Zhao, Y.; Harris, D.C.H.; Zheng, G. E-Cadherin/β-Catenin Complex and the Epithelial Barrier. *J. Biomed. Biotechnol.* **2011**, *2011*, 567305. [CrossRef]

26. Shang, S.; Hua, F.; Hu, Z.-W. The regulation of β-catenin activity and function in cancer: Therapeutic opportunities. *Oncotarget* **2017**, *8*, 33972–33989. [CrossRef] [PubMed]

27. Gilles, C.; Polette, M.; Mestdagt, M.; Nawrocki-Raby, B.; Ruggeri, P.; Birembaut, P.; Foidart, J.-M. Transactivation of vimentin by beta-catenin in human breast cancer cells. *Cancer Res.* **2003**, *63*, 2658–2664. [CrossRef] [PubMed]
28. Das, S.; Becker, B.N.; Hoffmann, F.M.; Mertz, J.E. Complete reversal of epithelial to mesenchymal transition requires inhibition of both ZEB expression and the Rho pathway. *BMC Cell Biol.* **2009**, *10*, 94. [CrossRef] [PubMed]

29. Zavadil, J.; Bitzer, M.; Liang, D.; Yang, Y.-C.; Massimi, A.; Kneitz, S.; Piek, E.; Böttinger, E.P. Genetic programs of epithelial cell plasticity directed by transforming growth factor-beta. *Proc. Natl. Acad. Sci. USA* **2001**, *98*, 6686–6691. [CrossRef] [PubMed]

30. Sánchez-Tilló, E.; Liu, Y.; De Barrios, O.; Siles, L.; Fanlo, L.; Cuatrecasas, M.; Darling, D.S.; Dean, D.C.; Castells, A.; Postigo, A. EMT-activating transcription factors in cancer: Beyond EMT and tumor invasiveness. *Cell. Mol. Life Sci.* **2012**, *69*, 3429–3456. [CrossRef]

31. Liu, Y.; Zhang, Y.; Wu, H.; Li, Y.; Zhang, Y.; Liu, M.; Li, X.; Tang, H. miR-10a suppresses colorectal cancer metastasis by modulating the epithelial-to-mesenchymal transition and anoikis. *Cell Death Dis.* **2017**, *8*, e2739. [CrossRef]

32. Zlobec, I.; Lugli, A.; Baker, K.; Roth, S.; Minoo, P.; Hayashi, S.; Terracciano, L.; Jass, J.R. Role of APAF-1, E-cadherin and peritumoural lymphocytic infiltration in tumour budding in colorectal cancer. *J. Pathol.* **2007**, *212*, 260–268. [CrossRef] [PubMed]

33. Starchenko, A.; Graves-Deal, R.; Yang, Y.-P.; Li, C.; Zent, R.; Singh, B.; Coffey, R.J. Clustering of integrin α5 at the lateral membrane restores epithelial polarity in invasive colorectal cancer cells. *Mol. Biol. Cell* **2017**, *28*, 1288–1300. [CrossRef] [PubMed]

34. Yuang, A.D.; Fan, F.; Camp, E.R.; Van Buren, G.; Liu, W.; Somcio, R.; Gray, M.J.; Cheng, H.; Hoff, P.M.; Ellis, L.M. Chronic Oxaliplatin Resistance Induces Epithelial-to-Mesenchymal Transition in Colorectal Cancer Cell Lines. *Clin. Cancer Res. Off. J. Am. Assoc. Cancer Res.* **2006**, *12*, 4147–4153. [CrossRef] [PubMed]

35. Cortez, M.A.; Valdecanas, D.; Zhang, X.; Zhan, Y.; Bhardwaj, V.; Calin, G.A.; Komaki, R.; Giri, D.K.; Quini, C.C.; Wolfe, T.; et al. Therapeutic Delivery of miR-200c Enhances Radiosensitivity in Lung Cancer. *Mol. Ther. J. Am. Soc. Gene Ther.* **2014**, *22*, 1494–1503. [CrossRef]

36. Mao, L.; Li, Y.; Zhao, J.; Li, Q.; Yang, B.; Wang, Y.; Zhu, Z.; Sun, H.; Zhai, Z. Transforming growth factor-β1 contributes to oxaliplatin resistance in colorectal cancer via epithelial to mesenchymal transition. *Oncol. Lett.* **2017**, *14*, 647–654. [CrossRef]

37. Mårtensson, A.; Oberg, A.; Jung, A.; Cederquist, K.; Stenling, R.; Palmqvist, R. Beta-catenin expression in relation to genetic instability and prognosis in colorectal cancer. *Oncol. Rep.* **2007**, *17*, 447–452.

38. Lugli, A.; Zlobec, I.; Minoo, P.; Baker, K.; Tornillo, L.; Terracciano, L.; Jass, J.R. Prognostic significance of the wnt signalling pathway molecules APC, beta-catenin and E-cadherin in colorectal cancer: A tissue microarray-based analysis. *Histopathology* **2007**, *50*, 453–464. [CrossRef]

39. Munemitsu, S.; Albert, I.; Souza, B.; Rubinfeld, B.; Polakis, P. Regulation of intracellular beta-catenin levels by the adenomatous polyposis coli (APC) tumor-suppressor protein. *Proc. Natl. Acad. Sci. USA* **1995**, *92*, 3046–3050. [CrossRef]

40. Henderson, B.R. Nuclear-cytoplasmic shuttling of APC regulates β-catenin subcellular localization and turnover. *Nat. Cell Biol.* **2000**, *2*, 653–660. [CrossRef] [PubMed]

41. Wang, Y.; Sun, N.; Lu, C.; Bei, Y.; Qian, R.; Hua, L. Upregulation of circadian gene 'hClock' contribution to metastasis of colorectal cancer. *Int. J. Oncol.* **2017**, *50*, 2191–2199. [CrossRef]

42. Sullivan, R.; Price, L.S.; Koffer, A. Rho Controls Cortical F-actin Disassembly in Addition to, but Independently of, Secretion in Mast Cells. *J. Biol. Chem.* **1999**, *274*, 38140–38146. [CrossRef]

43. Relógio, A.; Thomas, P.; Medina-Pérez, P.; Reischl, S.; Bervoets, S.; Gloc, E.; Riemer, P.; Mang-Fatehi, S.; Maier, B.; Schäfer, R.; et al. Ras-Mediated Deregulation of the Circadian Clock in Cancer. *PLoS Genet.* **2014**, *10*, e1004338. [CrossRef]

44. El-Athman, R.; Fuhr, L.; Relógio, A. A Systems-Level Analysis Reveals Circadian Regulation of Splicing in Colorectal Cancer. *EBioMedicine* **2018**, *33*, 68–81. [CrossRef]

45. Fuhr, L.; El-Athman, R.; Scrima, R.; Cela, O.; Carbone, A.; Knoop, H.; Li, Y.; Hoffmann, K.; Laukkanen, M.O.; Corcione, F.; et al. The Circadian Clock Regulates Metabolic Phenotype Rewiring Via HKDC1 and Modulates Tumor Progression and Drug Response in Colorectal Cancer. *EBioMedicine* **2018**, *33*, 105–121. [CrossRef] [PubMed]

46. Tan, C.; Du, X. KRAS mutation testing in metastatic colorectal cancer. *World J. Gastroenterol.* **2012**, *18*, 5171–5180.

47. Xu, W.; Yang, Z.; Lu, N. A new role for the PI3K/Akt signaling pathway in the epithelial-mesenchymal transition. *Cell Adhes. Migr.* **2015**, *9*, 317–324. [CrossRef]

48. Wei, R.; Xiao, Y.; Song, Y.; Yuan, H.; Luo, J.; Xu, W. FAT4 regulates the EMT and autophagy in colorectal cancer cells in part via the PI3K-AKT signaling axis. *J. Exp. Clin. Cancer Res.* **2019**, *38*, 1–14. [CrossRef] [PubMed]

49. Van Der Jeught, K.; Xu, H.-C.; Li, Y.-J.; Lu, X.-B.; Ji, G. Drug resistance and new therapies in colorectal cancer. *World J. Gastroenterol.* **2018**, *24*, 3834–3848. [CrossRef] [PubMed]

50. Cederroth, C.R.; Albrecht, U.; Bass, J.; Brown, S.A.; Dyhrfjeld-Johnsen, J.; Gachon, F.; Green, C.B.; Hastings, M.H.; Helfrich-Förster, C.; Hogenesch, J.B.; et al. Medicine in the Fourth Dimension. *Cell Metab.* **2019**, *30*, 238–250. [CrossRef]

51. Korkmaz, T.; Aygenli, F.; Emisoglu, H.; Ozcelik, G.; Canturk, A.; Yilmaz, S.; Ozturk, N. Opposite Carcinogenic Effects of Circadian Clock Gene BMAL1. *Sci. Rep.* **2018**, *8*, 16023. [CrossRef]

52. Sakamoto, W.; Takenoshita, S. Overexpression of both CLOCK and BMAL1 inhibits entry to S phase in human colon cancer cells. *Fukushima J. Med. Sci.* **2015**, *61*, 111–124. [CrossRef] [PubMed]

53. Zeng, Z.-L.; Luo, H.-Y.; Yang, J.; Wu, W.-J.; Chen, D.-L.; Huang, P.; Xu, R.-H. Overexpression of the Circadian Clock Gene Bmal1 Increases Sensitivity to Oxaliplatin in Colorectal Cancer. *Clin. Cancer Res.* **2014**, *20*, 1042–1052. [CrossRef]

54. Martinez-Balibrea, E.; Martínez-Cardús, A.; Ginés, A.; De Porras, V.R.; Moutinho, C.; Layos, L.; Manzano, J.L.; Bugés, C.; Bystrup, S.; Esteller, M.; et al. Tumor-Related Molecular Mechanisms of Oxaliplatin Resistance. *Mol. Cancer Ther.* **2015**, *14*, 1767–1776. [CrossRef] [PubMed]
28. Das, S.; Becker, B.N.; Hoffmann, F.M.; Mertz, J.E. Complete reversal of epithelial to mesenchymal transition requires inhibition of both ZEB expression and the Rho pathway. *BMC Cell Biol.* **2009**, *10*, 94. [CrossRef] [PubMed]

29. Zavadil, J.; Bitzer, M.; Liang, D.; Yang, Y.-C.; Massimi, A.; Kneitz, S.; Piek, E.; Böttinger, E.P. Genetic programs of epithelial cell plasticity directed by transforming growth factor-beta. *Proc. Natl. Acad. Sci. USA* **2001**, *98*, 6686–6691. [CrossRef] [PubMed]

30. Sánchez-Tilló, E.; Liu, Y.; De Barrios, O.; Siles, L.; Fanlo, L.; Cuatrecasas, M.; Darling, D.S.; Dean, D.C.; Castells, A.; Postigo, A. EMT-activating transcription factors in cancer: Beyond EMT and tumor invasiveness. *Cell. Mol. Life Sci.* **2012**, *69*, 3429–3456. [CrossRef]

31. Liu, Y.; Zhang, Y.; Wu, H.; Li, Y.; Zhang, Y.; Liu, M.; Li, X.; Tang, H. miR-10a suppresses colorectal cancer metastasis by modulating the epithelial-to-mesenchymal transition and anoikis. *Cell Death Dis.* **2017**, *8*, e2739. [CrossRef]

32. Zlobec, I.; Lugli, A.; Baker, K.; Roth, S.; Minoo, P.; Hayashi, S.; Terracciano, L.; Jass, J.R. Role of APAF-1, E-cadherin and peritumoural lymphocytic infiltration in tumour budding in colorectal cancer. *J. Pathol.* **2007**, *212*, 260–268. [CrossRef] [PubMed]

33. Starchenko, A.; Graves-Deal, R.; Yang, Y.-P.; Li, C.; Zent, R.; Singh, B.; Coffey, R.J. Clustering of integrin α5 at the lateral membrane restores epithelial polarity in invasive colorectal cancer cells. *Mol. Biol. Cell* **2017**, *28*, 1288–1300. [CrossRef] [PubMed]

34. Yuang, A.D.; Fan, F.; Camp, E.R.; Van Buren, G.; Liu, W.; Somcio, R.; Gray, M.J.; Cheng, H.; Hoff, P.M.; Ellis, L.M. Chronic Oxaliplatin Resistance Induces Epithelial-to-Mesenchymal Transition in Colorectal Cancer Cell Lines. *Clin. Cancer Res. Off. J. Am. Assoc. Cancer Res.* **2006**, *12*, 4147–4153. [CrossRef] [PubMed]

35. Cortez, M.A.; Valdecanas, D.; Zhang, X.; Zhan, Y.; Bhardwaj, V.; Calin, G.A.; Komaki, R.; Giri, D.K.; Quini, C.C.; Wolfe, T.; et al. Therapeutic Delivery of miR-200c Enhances Radiosensitivity in Lung Cancer. *Mol. Ther. J. Am. Soc. Gene Ther.* **2014**, *22*, 1494–1503. [CrossRef]

36. Mao, L.; Li, Y.; Zhao, J.; Li, Q.; Yang, B.; Wang, Y.; Zhu, Z.; Sun, H.; Zhai, Z. Transforming growth factor-β1 contributes to oxaliplatin resistance in colorectal cancer via epithelial to mesenchymal transition. *Oncol. Lett.* **2017**, *14*, 647–654. [CrossRef]

37. Mårtensson, A.; Oberg, A.; Jung, A.; Cederquist, K.; Stenling, R.; Palmqvist, R. Beta-catenin expression in relation to genetic instability and prognosis in colorectal cancer. *Oncol. Rep.* **2007**, *17*, 447–452.

38. Lugli, A.; Zlobec, I.; Minoo, P.; Baker, K.; Tornillo, L.; Terracciano, L.; Jass, J.R. Prognostic significance of the wnt signalling pathway molecules APC, beta-catenin and E-cadherin in colorectal cancer: A tissue microarray-based analysis. *Histopathology* **2007**, *50*, 453–464. [CrossRef]

39. Munemitsu, S.; Albert, I.; Souza, B.; Rubinfeld, B.; Polakis, P. Regulation of intracellular beta-catenin levels by the adenomatous polyposis coli (APC) tumor-suppressor protein. *Proc. Natl. Acad. Sci. USA* **1995**, *92*, 3046–3050. [CrossRef]

40. Henderson, B.R. Nuclear-cytoplasmic shuttling of APC regulates β-catenin subcellular localization and turnover. *Nat. Cell Biol.* **2000**, *2*, 653–660. [CrossRef] [PubMed]

41. Wang, Y.; Sun, N.; Lu, C.; Bei, Y.; Qian, R.; Hua, L. Upregulation of circadian gene 'hClock' contribution to metastasis of colorectal cancer. *Int. J. Oncol.* **2017**, *50*, 2191–2199. [CrossRef]

42. Sullivan, R.; Price, L.S.; Koffer, A. Rho Controls Cortical F-actin Disassembly in Addition to, but Independently of, Secretion in Mast Cells. *J. Biol. Chem.* **1999**, *274*, 38140–38146. [CrossRef]

43. Relógio, A.; Thomas, P.; Medina-Pérez, P.; Reischl, S.; Bervoets, S.; Gloc, E.; Riemer, P.; Mang-Fatehi, S.; Maier, B.; Schäfer, R.; et al. Ras-Mediated Deregulation of the Circadian Clock in Cancer. *PLoS Genet.* **2014**, *10*, e1004338. [CrossRef]

44. El-Athman, R.; Fuhr, L.; Relógio, A. A Systems-Level Analysis Reveals Circadian Regulation of Splicing in Colorectal Cancer. *EBioMedicine* **2018**, *33*, 68–81. [CrossRef]

45. Fuhr, L.; El-Athman, R.; Scrima, R.; Cela, O.; Carbone, A.; Knoop, H.; Li, Y.; Hoffmann, K.; Laukkanen, M.O.; Corcione, F.; et al. The Circadian Clock Regulates Metabolic Phenotype Rewiring Via HKDC1 and Modulates Tumor Progression and Drug Response in Colorectal Cancer. *EBioMedicine* **2018**, *33*, 105–121. [CrossRef] [PubMed]

46. Tan, C.; Du, X. KRAS mutation testing in metastatic colorectal cancer. *World J. Gastroenterol.* **2012**, *18*, 5171–5180.

47. Xu, W.; Yang, Z.; Lu, N. A new role for the PI3K/Akt signaling pathway in the epithelial-mesenchymal transition. *Cell Adhes. Migr.* **2015**, *9*, 317–324. [CrossRef]

48. Wei, R.; Xiao, Y.; Song, Y.; Yuan, H.; Luo, J.; Xu, W. FAT4 regulates the EMT and autophagy in colorectal cancer cells in part via the PI3K-AKT signaling axis. *J. Exp. Clin. Cancer Res.* **2019**, *38*, 1–14. [CrossRef] [PubMed]

49. Van Der Jeught, K.; Xu, H.-C.; Li, Y.-J.; Lu, X.-B.; Ji, G. Drug resistance and new therapies in colorectal cancer. *World J. Gastroenterol.* **2018**, *24*, 3834–3848. [CrossRef] [PubMed]

50. Cederroth, C.R.; Albrecht, U.; Bass, J.; Brown, S.A.; Dyhrfjeld-Johnsen, J.; Gachon, F.; Green, C.B.; Hastings, M.H.; Helfrich-Förster, C.; Hogenesch, J.B.; et al. Medicine in the Fourth Dimension. *Cell Metab.* **2019**, *30*, 238–250. [CrossRef]

51. Korkmaz, T.; Aygenli, F.; Emisoglu, H.; Ozcelik, G.; Canturk, A.; Yilmaz, S.; Ozturk, N. Opposite Carcinogenic Effects of Circadian Clock Gene BMAL1. *Sci. Rep.* **2018**, *8*, 16023. [CrossRef]

52. Sakamoto, W.; Takenoshita, S. Overexpression of both CLOCK and BMAL1 inhibits entry to S phase in human colon cancer cells. *Fukushima J. Med. Sci.* **2015**, *61*, 111–124. [CrossRef] [PubMed]

53. Zeng, Z.-L.; Luo, H.-Y.; Yang, J.; Wu, W.-J.; Chen, D.-L.; Huang, P.; Xu, R.-H. Overexpression of the Circadian Clock Gene Bmal1 Increases Sensitivity to Oxaliplatin in Colorectal Cancer. *Clin. Cancer Res.* **2014**, *20*, 1042–1052. [CrossRef]

54. Martinez-Balibrea, E.; Martínez-Cardús, A.; Ginés, A.; De Porras, V.R.; Moutinho, C.; Layos, L.; Manzano, J.L.; Bugés, C.; Bystrup, S.; Esteller, M.; et al. Tumor-Related Molecular Mechanisms of Oxaliplatin Resistance. *Mol. Cancer Ther.* **2015**, *14*, 1767–1776. [CrossRef] [PubMed]
55. Deng, J.-J.; Zhang, W.; Xu, X.-M.; Zhang, F.; Tao, W.-P.; Ye, J.-J.; Ge, W. Twist mediates an aggressive phenotype in human colorectal cancer cells. *Int. J. Oncol.* **2016**, *48*, 1117–1124. [CrossRef] [PubMed]

56. Zhao, X.; Fan, J.; Zhi, F.; Li, A.; Li, C.; Berger, A.E.; Boorgula, M.P.; Barkataki, S.; Journeya, J.-P.; Chen, Y.; et al. Mobilization of epithelial mesenchymal transition genes distinguishes active from inactive lesional tissue in patients with ulcerative colitis. *Hum. Mol. Genet.* **2015**, *24*, 4615–4624. [CrossRef] [PubMed]

57. Hwang-Verslues, W.W.; Chang, P.-H.; Jeng, Y.-M.; Kuo, W.-H.; Chiang, P.-H.; Chang, Y.-C.; Hsieh, T.-H.; Su, F.-Y.; Lin, L.-C.; Abbondante, S.; et al. Loss of corepressor PER2 under hypoxia up-regulates OCT1-mediated EMT gene expression and enhances tumor malignancy. *Proc. Natl. Acad. Sci. USA* **2013**, *110*, 12331–12336. [CrossRef] [PubMed]

58. De, A.; Beligala, D.H.; Sharma, V.P.; Burgos, C.A.; Lee, A.M.; Geusz, M.E. Cancer stem cell generation during epithelial-mesenchymal transition is temporally gated by intrinsic circadian clocks. *Clin. Exp. Metastasis* **2020**, *37*, 617–635. [CrossRef]

59. De, A.; Beligala, D.H.; Sharma, V.P.; Fry, B.R.; Geusz, M.E. Abstract 858: The circadian clock of glioma cells undergoing epithelial-mesenchymal transition. *Cancer Res.* **2017**, *77*, 858. [CrossRef]

60. Ma, S.; Wang, Z.; Cao, J.; Dong, Y.; Chen, Y. BMAL1 but not CLOCK is associated with monochromatic green light-induced circadian rhythm of melatonin in chick pinealocytes. *Endocr. Connect.* **2019**, *8*, 57–68. [CrossRef]

61. Cancer Genome Atlas Network. Comprehensive molecular characterization of human colon and rectal cancer. *Nature* **2012**, *487*, 330–337. [CrossRef] [PubMed]

62. Gao, J.; Aksoy, B.A.; Dogrusoz, U.; Dresdner, G.; Gross, B.; Sumer, S.O.; Sun, Y.; Jacobsen, A.; Sinha, R.; Larsson, E.; et al. Integrative Analysis of Complex Cancer Genomics and Clinical Profiles Using the cBioPortal. *Sci. Signal.* **2013**, *6*, pl1. [CrossRef] [PubMed]

63. Subramanian, A.; Tamayo, P.; Mootha, V.K.; Mukherjee, S.; Ebert, B.L.; Gillette, M.A.; Paulovich, A.; Pomeroy, S.L.; Golub, T.R.; Lander, E.S.; et al. Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles. *Proc. Natl. Acad. Sci. USA* **2005**, *102*, 15545–15550. [CrossRef]

64. Liberzon, A.; Birger, C.; Thorvaldsdóttir, H.; Ghandi, M.; Mesirov, J.P.; Tamayo, P. The Molecular Signatures Database (MSigDB) Hallmark Gene Set Collection. *Cell Syst.* **2015**, *1*, 417–425. [CrossRef]

65. Wickham, H. *ggplot2: Elegant Graphics for Data Analysis*; Springer: New York, NY, USA, 2009.

66. Lê, S.; Josse, J.; Husson, F. FactoMineR: An R Package for Multivariate Analysis. *J. Stat. Softw.* **2008**, *25*, 1–18. [CrossRef]

67. Lordier, L.; Bluteau, D.; Jalil, A.; Legrand, C.; Pan, J.; Rameau, P.; Jouni, D.; Bluteau, O.; Mercher, T.; Leon, C.; et al. RUNX1-induced silencing of non-muscle myosin heavy chain IIB contributes to megakaryocyte polyploidization. *Nat. Commun.* **2012**, *3*, 717. [CrossRef] [PubMed]

68. Zhang, Y.; Giacchetti, S.; Parouchev, A.; Hadadi, E.; Li, X.; Dallmann, R.; Xandri-Monje, H.; Portier, L.; Adam, R.; Lévi, F.; et al. Dosing time dependent in vitro pharmacodynamics of Everolimus despite a defective circadian clock. *Cell Cycle* **2018**, *17*, 33–42. [CrossRef]
55. Deng, J.-J.; Zhang, W.; Xu, X.-M.; Zhang, F.; Tao, W.-P.; Ye, J.-J.; Ge, W. Twist mediates an aggressive phenotype in human colorectal cancer cells. *Int. J. Oncol.* **2016**, *48*, 1117–1124. [CrossRef] [PubMed]

56. Zhao, X.; Fan, J.; Zhi, F.; Li, A.; Li, C.; Berger, A.E.; Boorgula, M.P.; Barkataki, S.; Journeya, J.-P.; Chen, Y.; et al. Mobilization of epithelial mesenchymal transition genes distinguishes active from inactive lesional tissue in patients with ulcerative colitis. *Hum. Mol. Genet.* **2015**, *24*, 4615–4624. [CrossRef] [PubMed]

57. Hwang-Verslues, W.W.; Chang, P.-H.; Jeng, Y.-M.; Kuo, W.-H.; Chiang, P.-H.; Chang, Y.-C.; Hsieh, T.-H.; Su, F.-Y.; Lin, L.-C.; Abbondante, S.; et al. Loss of corepressor PER2 under hypoxia up-regulates OCT1-mediated EMT gene expression and enhances tumor malignancy. *Proc. Natl. Acad. Sci. USA* **2013**, *110*, 12331–12336. [CrossRef] [PubMed]

58. De, A.; Beligala, D.H.; Sharma, V.P.; Burgos, C.A.; Lee, A.M.; Geusz, M.E. Cancer stem cell generation during epithelial-mesenchymal transition is temporally gated by intrinsic circadian clocks. *Clin. Exp. Metastasis* **2020**, *37*, 617–635. [CrossRef]

59. De, A.; Beligala, D.H.; Sharma, V.P.; Fry, B.R.; Geusz, M.E. Abstract 858: The circadian clock of glioma cells undergoing epithelial-mesenchymal transition. *Cancer Res.* **2017**, *77*, 858. [CrossRef]

60. Ma, S.; Wang, Z.; Cao, J.; Dong, Y.; Chen, Y. BMAL1 but not CLOCK is associated with monochromatic green light-induced circadian rhythm of melatonin in chick pinealocytes. *Endocr. Connect.* **2019**, *8*, 57–68. [CrossRef]

61. Cancer Genome Atlas Network. Comprehensive molecular characterization of human colon and rectal cancer. *Nature* **2012**, *487*, 330–337. [CrossRef] [PubMed]

62. Gao, J.; Aksoy, B.A.; Dogrusoz, U.; Dresdner, G.; Gross, B.; Sumer, S.O.; Sun, Y.; Jacobsen, A.; Sinha, R.; Larsson, E.; et al. Integrative Analysis of Complex Cancer Genomics and Clinical Profiles Using the cBioPortal. *Sci. Signal.* **2013**, *6*, pl1. [CrossRef] [PubMed]

63. Subramanian, A.; Tamayo, P.; Mootha, V.K.; Mukherjee, S.; Ebert, B.L.; Gillette, M.A.; Paulovich, A.; Pomeroy, S.L.; Golub, T.R.; Lander, E.S.; et al. Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles. *Proc. Natl. Acad. Sci. USA* **2005**, *102*, 15545–15550. [CrossRef]

64. Liberzon, A.; Birger, C.; Thorvaldsdóttir, H.; Ghandi, M.; Mesirov, J.P.; Tamayo, P. The Molecular Signatures Database (MSigDB) Hallmark Gene Set Collection. *Cell Syst.* **2015**, *1*, 417–425. [CrossRef]

65. Wickham, H. *ggplot2: Elegant Graphics for Data Analysis*; Springer: New York, NY, USA, 2009.

66. Lê, S.; Josse, J.; Husson, F. FactoMineR: An R Package for Multivariate Analysis. *J. Stat. Softw.* **2008**, *25*, 1–18. [CrossRef]

67. Lordier, L.; Bluteau, D.; Jalil, A.; Legrand, C.; Pan, J.; Rameau, P.; Jouni, D.; Bluteau, O.; Mercher, T.; Leon, C.; et al. RUNX1-induced silencing of non-muscle myosin heavy chain IIB contributes to megakaryocyte polyploidization. *Nat. Commun.* **2012**, *3*, 717. [CrossRef] [PubMed]

68. Zhang, Y.; Giacchetti, S.; Parouchev, A.; Hadadi, E.; Li, X.; Dallmann, R.; Xandri-Monje, H.; Portier, L.; Adam, R.; Lévi, F.; et al. Dosing time dependent in vitro pharmacodynamics of Everolimus despite a defective circadian clock. *Cell Cycle* **2018**, *17*, 33–42. [CrossRef]