# A dual model of normal vs isogenic Nrf2-depleted murine epithelial cells to explore oxidative stress involvement

Jacques Dupuy, Edwin Fouché, Céline Noirot, Pierre Martin, Charline Buisson, Françoise Guéraud, Fabrice Pierre, Cécile Héliès-Toussaint

► To cite this version:

Jacques Dupuy, Edwin Fouché, Céline Noirot, Pierre Martin, Charline Buisson, et al.. A dual model of normal vs isogenic Nrf2-depleted murine epithelial cells to explore oxidative stress involvement.
Scientific Reports, 2024, 14, 10.1038/s41598-024-60938-2. hal-04677587

HAL Id: hal-04677587

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

Submitted on 26 Aug 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_1.png" /></figure>

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

Check for updates

OPEN

# A dual model of normal vs isogenic Nrf2-depleted murine epithelial cells to explore oxidative stress involvement

Jacques Dupuy<fnref n="1" />, Edwin Fouché<fnref n="1" />, Céline Noirot², Pierre Martin², Charline Buisson<fnref n="1" />,
Françoise Guéraud<fnref n="1" />, Fabrice Pierre<fnref n="1" /> & Cécile Héliès-Toussaint<fnref n="1" />,✉

Cancer-derived cell lines are useful tools for studying cellular metabolism and xenobiotic toxicity, but they are not suitable for modeling the biological effects of food contaminants or natural biomolecules on healthy colonic epithelial cells in a normal genetic context. The toxicological properties of such compounds may rely on their oxidative properties. Therefore, it appears to be necessary to develop a dual-cell model in a normal genetic context that allows to define the importance of oxidative stress in the observed toxicity. Given that the transcription factor nuclear factor erythroid 2-related factor 2 (Nrf2) is considered to be the master regulator of antioxidant defenses, our aim was to develop a cellular model comparing normal and Nrf2-depleted isogenic cells to qualify oxidative stress-related toxicity. We generated these cells by using the CRISPR/Cas9 technique. Whole-genome sequencing enabled us to confirm that our cell lines were free of cancer-related mutations. We used 4-hydroxy-2-nonenal (HNE), a lipid peroxidation product closely related to oxidative stress, as a model molecule. Here we report significant differences between the two cell lines in glutathione levels, gene regulation, and cell viability after HNE treatment. The results support the ability of our dual-cell model to study the role of oxidative stress in xenobiotic toxicity.

**Keywords** Nrf2, Oxidative stress, Dual-cell model, HNE

Intestinal cells are confronted to a large number of xenobiotics, including environmental and dietary contaminants, neo-formed products of the diet, and biomolecules that are naturally produced during digestion. These compounds can exhibit cytotoxic and genotoxic properties in the colon, mainly due to the generation of reactive oxygen species (ROS) and lipid peroxidation products. Under normal conditions, ROS are essential for normal life as they are responsible for maintaining cellular redox homeostasis<fnref n="1" /> and have been considered major signaling molecules²,³. However, high ROS concentrations lead to an imbalance between ROS production and the cellular antioxidant systems, resulting in oxidative stress. When cellular repair capacities are exceeded, excessive levels of ROS, electrophiles, or prooxidant molecules can lead to protein or DNA adducts that could contribute to increased oncogenic risk⁴. Given that the transcription factor nuclear factor erythroid 2-related factor 2 (Nrf2) is a master regulator of ROS metabolism and the antioxidant defense, and that its induction leads to protection against oxidative damage⁵⁻⁷, we thought it would be of great interest to assess the relative weight of oxidative stress in xenobiotic management.

Due to their availability and ease of use, cell lines represent a tool of choice in xenobiotic toxicity studies. Although cell lines are very useful for modeling cancer promotion and progression, their carcinogenic phenotype and genotype mean that they are poorly suited for studying the biological effects of food contaminants on healthy colonic mucosa. Therefore, we believed it would be appropriate to develop a cellular model with a controlled and normal genetic background to avoid mutations due to carcinogenic transformation, which could distort mechanical conclusions on contaminant toxicity. We used a normal murine cell model to assess contaminant toxicity in a normal, controlled genetic context. We developed this cellular model from conditionally immortalized colonic

<footnote n="1">¹National Research Institute for Agriculture and Environment (INRAE), Toxalim (Research Centre in Food Toxicology), INRAE, ENVT, INP-Purpan, UPS, Université de Toulouse, 180 Chemin de Tournefeuille, BP93173, 31027 Toulouse Cedex 3, France. ²National Research Institute for Agriculture and Environment (INRAE), Université Fédérale de Toulouse, INRAE, BioinfOmics, GenoToul Bioinformatics Facility, 31326 Castanet-Tolosan, France. ✉email: cecile.toussaint@inrae.fr</footnote>
epithelial cell lines derived from C57BL/6J mice⁸, and generated an isogenic Nrf2-knockdown cell line by using the CRISPR/Cas9 technique. The normal genetic background of these cell lines was confirmed by whole-genome sequencing (WGS). After generating this dual-cell model, we used 4-hydroxy-2-nonenal (HNE), a molecule formed in the intestine and colonic lumen in a pro-peroxidant context—for example, after a diet high in heme iron and low in antioxidants⁹—as a model of reactive product of lipid peroxidation that induces oxidative stress. The effects of HNE on normal colon cells and their isogenic counterparts lacking Nrf2 enabled us to identify the burden of oxidative stress in cellular responses to xenobiotic attack in our murine epithelial cells.

# Results

## Identification of a normal colonic epithelial cell line

The murine colonic epithelial cell line used in these experiments (named Co) is derived from the Apc$^{+/+}$ clones described by Forest et al.⁸. In these experiments, we used WGS to determine the non-precancerous and non-cancerous status of the Co cell line, with a particular focus on pro-oncogenes, tumor suppressor genes, and high-impact mutations (see Table 1a for sample information). More specifically, we studied mutations in a panel of genes from the main pathways linked to colorectal cancer: the WNT and NOTCH pathways, the c-myc gene, the PI3K and P53 genes and other cancer-related genes, Nrf2 and the associated Nrf1 and Nrf3 genes, the Kelch-like ECH-associated protein 1 (Keap1) gene, and genes involved in cell cycle regulation. Only mutations

<table><caption>Table 1. Sequence comparison analyses.</caption><thead><tr><th colspan="10">(a) Raw data</th></tr><tr><th>Sample name</th><th>Raw reads</th><th>% Duplicates</th><th>Reads aligned</th><th>≥ 10X</th><th>Mean cov</th><th>Vars</th><th>SNP</th><th>Indel</th></tr></thead><tbody><tr><td>Co</td><td>481,788,378</td><td>29.6%</td><td>479,112,014</td><td>45.03</td><td>17.39</td><td>1,354,919</td><td>993,532</td><td>362,827</td></tr></tbody></table><table><thead><tr><th colspan="10">(b) Variant detection in a panel of genes from the main colorectal cancer-related pathways</th></tr><tr><th>Pathways</th><th>Genes tested</th><th>Transcript length</th><th>Average depth (Mosdeth)</th><th>Feature count on exons</th><th>Gene filter</th><th>Moderate</th><th>High</th><th></th></tr></thead><tbody><tr><td rowspan="5">WNT</td><td>APC</td><td>15,753</td><td>14.27</td><td>1710</td><td>9</td><td></td><td></td><td></td></tr><tr><td>CTNNb1/Beta catenin</td><td>5099</td><td>50.23</td><td>736</td><td>5</td><td>1</td><td></td><td></td></tr><tr><td>GSK3 beta</td><td>9821</td><td>11.49</td><td>3248</td><td>10</td><td></td><td></td><td></td></tr><tr><td>Axin 1</td><td>3915</td><td>28.44</td><td>1662</td><td>14</td><td></td><td></td><td></td></tr><tr><td>Axin 2</td><td>5506</td><td>36.50</td><td>2768</td><td>7</td><td></td><td></td><td></td></tr><tr><td rowspan="4">NOTCH</td><td>NOTCH 1</td><td>14,232</td><td>38.73</td><td>7990</td><td>12</td><td>1</td><td>1</td><td></td></tr><tr><td>NOTCH 2</td><td>10,777</td><td>24.76</td><td>4501</td><td>20</td><td></td><td></td><td></td></tr><tr><td>NOTCH 3</td><td>8044</td><td>64.99</td><td>6399</td><td>10</td><td>1</td><td></td><td></td></tr><tr><td>NOTCH 4</td><td>8080</td><td>35.11</td><td>5541</td><td>43</td><td></td><td></td><td></td></tr><tr><td>myc</td><td>c-myc</td><td>2767</td><td>35.26</td><td>1155</td><td>1</td><td></td><td></td><td></td></tr><tr><td rowspan="4">TGF beta</td><td>SMAD 2</td><td>2037</td><td>11.51</td><td>104</td><td>4</td><td></td><td></td><td></td></tr><tr><td>SMAD 3</td><td>5542</td><td>64.97</td><td>883</td><td>11</td><td></td><td></td><td></td></tr><tr><td>SMAD 4</td><td>3741</td><td>21.34</td><td>631</td><td>12</td><td>1</td><td></td><td></td></tr><tr><td>BMP4</td><td>2159</td><td>75.70</td><td>1411</td><td>3</td><td></td><td></td><td></td></tr><tr><td>Hippo</td><td>Yap1</td><td>9257</td><td>13.06</td><td>1978</td><td>38</td><td></td><td></td><td></td></tr><tr><td>PI3K</td><td>Akt1</td><td>3948</td><td>14.27</td><td>2955</td><td></td><td></td><td></td><td></td></tr><tr><td>P53</td><td>BRCA1</td><td>8533</td><td>13.15</td><td>786</td><td>9</td><td></td><td></td><td></td></tr><tr><td rowspan="6">Cell cycle</td><td>TRP 53</td><td>4056</td><td>29.77</td><td>1069</td><td>6</td><td></td><td></td><td></td></tr><tr><td>Cdk2</td><td>2432</td><td>71.45</td><td>1385</td><td>2</td><td></td><td></td><td></td></tr><tr><td>Cdk4</td><td>3249</td><td>13.52</td><td>1978</td><td>1</td><td></td><td></td><td></td></tr><tr><td>Cdk6</td><td>3151</td><td>33.03</td><td>1132</td><td>2678</td><td>12</td><td>3</td><td></td></tr><tr><td>E2f1</td><td>3150</td><td>87.34</td><td>2670</td><td>13</td><td></td><td></td><td></td></tr><tr><td>E2f3</td><td>4794</td><td>38.12</td><td>790</td><td>8</td><td></td><td></td><td></td></tr><tr><td rowspan="6">Others</td><td>cdkn2a</td><td>1172</td><td>48.88</td><td>893</td><td>9</td><td></td><td></td><td></td></tr><tr><td>B raf</td><td>13,863</td><td>12.48</td><td>1755</td><td>20</td><td></td><td></td><td></td></tr><tr><td>K ras</td><td>5347</td><td>25.54</td><td>1035</td><td>1</td><td></td><td></td><td></td></tr><tr><td>PTEN</td><td>10,632</td><td>16.14</td><td>1060</td><td>7</td><td></td><td></td><td></td></tr><tr><td>COX-2</td><td>4580</td><td>12.54</td><td>672</td><td></td><td></td><td></td><td></td></tr><tr><td>Fas</td><td>1801</td><td>14.84</td><td>169</td><td>10</td><td></td><td></td><td></td></tr><tr><td rowspan="4">NRF2</td><td>Bcl2</td><td>9129</td><td>7.14</td><td>2195</td><td>13</td><td></td><td></td><td></td></tr><tr><td>Nrf1</td><td>5350</td><td>83.79</td><td>3783</td><td>2</td><td>1</td><td></td><td></td></tr><tr><td>Nrf2</td><td>3796</td><td>20.67</td><td>702</td><td>4</td><td></td><td></td><td></td></tr><tr><td>Nrf3</td><td>2973</td><td>43.92</td><td>1663</td><td>8</td><td></td><td></td><td></td></tr><tr><td></td><td>Keap1</td><td>4664</td><td>24.99</td><td>2978</td><td>7</td><td>1</td><td></td><td></td></tr></tbody></table>
with a “high” impact, defined by McLaren et al.¹⁰ as “the variant is assumed to have a disruptive impact in the protein, probably causing protein truncation, loss of function or triggering nonsense-mediated decay,” and a moderate impact, “a non-disruptive variant that might change protein effectiveness,” were considered. With the exception of Notch1 and Cdk6, there are no “high” mutations in the Co cell line (Table 1b; Table S4), confirming its non-cancerous status.

## Creation of a Nrf2-knockdown clonal cell line using a CRISPR-directed gene editing approach

### Clone selection

We used the CRISPR/Cas9 technique to disable the Nrf2 gene in the Co cell line, rendering it incapable of producing a functional protein. We detected Nrf2 gene mutations with polymerase chain reaction (PCR), using primers specific to each guide RNA (Supplementary Fig. S1a). After prescreening to identify mismatches, we selected four clones and analyzed Nrf2 messenger RNA (mRNA) expression with reverse transcription-quantitative polymerase chain reaction (RT-qPCR) to select the clones with a substantial level of Nrf2 invalidation. We used hypoxanthine phosphoribosyl transferase 1 (Hprt1) as a control gene. We selected one clone, named Co-dN (for colon cell depleted for Nrf2) in which the Nrf2 mRNA level was 60% lower than in the parental Co cell line (Fig. 1a).

### NRF2 protein expression

We assessed Nrf2 protein expression with a western blot of total protein extracts from the Co and Co-dN cell lines, treated or not with HNE to induce Nrf2 expression (Fig. 1b; S2a). Vinculin served as an internal control (Fig. 1b; S2b). There were two bands in the Co cell line (around 75 and 90 kDa). The 90 kDa band was more

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

**Figure 1.** Generation and characterization of a normal mouse colonic epithelial cell line in which the Nrf2 gene is disabled. (a) Relative mRNA expression level of *Nrf2*, normalized to the *Hprt1* gene in Co and Co-dN cells. (b) Western blot analysis of Nrf2 protein expression in Co and Co-dN cells without or with HNE treatment (40 µM for 6 h); vinculin (124 kDa) was used as an internal control. (b1) Molecular weight marker; Nrf2 protein in Co (+/- HNE) and Co-dN (+/- HNE) cell extracts. (b2) Molecular weight marker; vinculin protein in Co (+/- HNE) and Co-dN (+/- HNE) cell extracts. (c). Reported slope before HNE treatment of both cell types based on xCELLigence data. (d) Relative mRNA expression level of *Hmox1*, a target gene of the transcription factor Nrf2. The detoxification pathway was induced by treating cells with HNE (40 µM for 6 h). Gene expression was normalized to the control gene *Hprt1*. For each figure, black boxes are Co cells, open boxes are Co-dN cells, and error bars indicate the mean ± standard error of the mean (n = 3 for **a** and **d**, and n = 5 for **c**). Student t-test (**a**, **c**) or two-way ANOVA followed by Tukey's multiple comparison test (**d**) were used for statistical analysis (ns = non-significant, \*P < 0.05, \*\*\*P < 0.0001).
intense in HNE-treated Co cells (about 200%, Supplementary Table S2). There was no band in Co-dN cells, without or with HNE treatment. The 90 kDa band is specific to the Nrf2 protein, as reported by Kopacz et al.¹¹, while the 75 kDa band corresponds to a non-specific band. Taken together, we detected the Nrf2 protein in Co cells, induced after HNE treatment, a finding that correlates with heme oxygenase 1 (*Hmox1*) mRNA expression (Fig. 1d). On the other hand, there was no Nrf2-specific band in the Co-dN cells even after HNE stimulation, confirming that the active form of Nrf2 had been genetically disabled.

### Cellular proliferation

The xCELLigence technique enables continuous monitoring of morphological changes in adherent cells¹². We used this technique to analyze cell proliferation in real time. The initial mean slope of the proliferation curve was 0.297 ± 0.014 h⁻¹ for Co cells and 0.128 ± 0.006 h⁻¹ for Co-dN cells. Thus, there was a significant (P < 0.0001) 57% reduction in the proliferation rate of the Nrf2-depleted Co-dN cell line (Fig. 1c). The difference between the two groups was analyzed with Student's t-test.

### Nrf2 target gene induction

We analyzed the modulation of the Nrf2 pathway via a target gene: *Hmox1*. We treated Co and Co-dN cells with HNE (40 µM for 6 h) to activate the Nrf2-mediated HNE detoxification pathway. *Hmox1* was strongly induced by HNE in Co cells (43-fold, P < 0.0001), whereas induction was 10 times weaker in Co-dN cells (4.3-fold, non-significant; Fig. 1d). The differences were assessed with two-way analysis of variance (ANOVA) followed by Tukey's multiple comparison test.

### Sequence analysis

We performed WGS on the Co-dN cell line. The IGV software allowed us to visualize insertion/deletion (indel) in the Nrf2 gene. There is a deletion (Fig. 2a) in exon 2 (Fig. 2b). The search for variants between the initial mouse genotype and the exomes of the two cell lines revealed few differences (Fig. 2c, d; Supplementary Table S3). The Venn diagram (Fig. 2d) shows the total number of variants, which are referenced in supplementary Table S4. A large number of variants: 12 and 22 variants for Co and Co-dN respectively are associated to a predicted gene 16505 (Gm 16505, MGP_C57BL6NJ_G0029685) without any identified function (Table S4). In the Co cell line, 3 genes presented mutations, Notch 1 and 3 and Smad4 but only Notch1 was associated with high impact mutation. In the Co-dN cell line, the *Cdk6* mutation already identified in Co cells is still detectable, confirming that the two cell lines are isogenic. The Notch4 and Trp53 mutations were considered as moderate. Only Nrf2 was associated with high impact mutation (protein modification). There are no other “high” mutations in Co-dN cells, indicating that in our experiment, the CRISPR/Cas9 technique did not induce off-target effects.

### Cellular response to HNE

#### Cellular detoxification

To investigate the role of Nrf2 in the cellular response to toxic lipid oxidation products, we used RT-qPCR to analyze the expression of three Nrf2-target genes involved in HNE detoxification: glutathione-S-transferase A4 (*Gsta4*), aldo-keto reductase 1b8 (*Akr1b8*), and aldehyde dehydrogenase 2 (*Aldh2*)¹³ in both cell lines. We used HNE treatment (40 µM for 6 h) to induce its own detoxification pathway. The results are presented in Fig. 3a-c. Under control conditions, there was no difference between the cell lines in *Gsta4* and *Akr1b8* expression, but Nrf2 depletion resulted in significantly lower *Aldh2* expression (P < 0.05). Treating the cells with HNE resulted in overexpression of *Gsta4* (eightfold, P < 0.0001) and *Akr1b8* (2.5-fold, P < 0.01) in Co cells, whereas it did not affect *Aldh2* expression. On the other hand, HNE did not affect the expression of these genes in Co-dN cells. Two-way ANOVA followed by Tukey's multiple comparison test was used for statistical analysis.

#### HPLC profiling and metabolites identification

We treated both cell lines for 30 min with 40 µM HNE containing a radioactive tracer ([³H]-HNE) and analyzed the supernatant with radio-high-performance liquid chromatography (HPLC). There are peaks corresponding to HNE and its metabolites, including isomers of glutathione conjugates (HNE-GSH), a mixture of oxidized metabolites (Ox), the saturated derivative of HNE (sat-HNE), and the original HNE molecule, as reported by Baradat et al.¹⁴. The results, presented as the percentages of observed peaks, showed that HNE metabolites were more abundant in Co cells (P < 0.05), while the remaining original HNE molecule (+32%) and its saturated derivatives (+14%) were more abundant in Co-dN cells (Fig. 3d). Two-way ANOVA followed by Tukey's multiple comparison test was used for statistical analysis.

#### Oxidative stress response associated with GSH

Based on its antioxidant properties, we pretreated cells with curcumin (CCM) for 24 h to induce the cellular antioxidant capacity¹⁵. Furthermore, CCM is a known inducer of Nrf2¹⁶⁻¹⁸. The cellular antioxidant capacity is linked to the intracellular GSH content. We determined the GSH/GSSG ratio in both cell lines and after treatments (40 µM HNE and/or 10 µM CCM pretreatment). There was an overall lower GSH/GSSG ratio in Co-dN cells compared with Co cells (Fig. 4a, P < 0.01, -27% for CTRL, -50% for 40 µM HNE, -34% for 10 µM CCM, and -17% for CCM pretreatment and HNE treatment). These differences reflected greater oxidative stress in Co-dN cells. HNE treatment resulted in a significant decrease in the GSH/GSSG ratio in both cell lines, revealing its oxidative stress effect (-46% in Co, P < 0.01; -63% in Co-dN, P < 0.01). Pretreatment with CCM 24 h before HNE treatment did not protect Co or Co-dN cells from HNE-induced oxidative stress (no significant difference between HNE and CCM).
<figure><img src="image_4.png" /><figcaption>b</figcaption></figure>

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

<table><caption>d</caption><thead><tr><th></th><th></th><th>Co</th><th>Co-dN</th></tr></thead><tbody><tr><td>Out of nextflow</td><td></td><td>1,354,919</td><td>1,562,679</td></tr><tr><td rowspan="3">All variant<br/>qual&gt; 30 && dp&gt;5</td><td>All</td><td>1,246,185</td><td>1,546,056</td></tr><tr><td>High</td><td>1,12</td><td>1,434</td></tr><tr><td>Moderate</td><td>6,783</td><td>9,117</td></tr><tr><td rowspan="3">Gene of interest</td><td>Variant</td><td>3,317</td><td>3,649</td></tr><tr><td>Effect High</td><td>9</td><td>12</td></tr><tr><td>Moderate effect</td><td>25</td><td>36</td></tr></tbody></table>

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

**Figure 2.** Identification of the deletion of the Nrf2 gene and exploration of variants. (a) Identification of insertion/deletion in the Nrf2 gene, visualized using the IGV software and gRNA localization. (b) Variants predictor results (extracted from e!Ensembl) showing the “high”-impact mutation in exon 2 of the *Nfe2l2* cDNA. (c) The number of variants per gene of interest and per effect. (d) Venn diagram representing common variants with high and moderate effects between the two cell lines.

GSH synthesis depends on the regulation of glutamate-cysteine ligase (GCL). Therefore, we analyzed the regulation of the genes that encode the two subunits of this enzyme (Gclc and Gclm) with RT-qPCR (Fig. 4 b, c). There were no differences between Co and Co-dN cells in the control condition. HNE treatment resulted in a strong induction of both genes in Co cells (P < 0.001). This induction was completely lost in Co-dN cells (P < 0.001).

The balance between cystine uptake and glutamate release, via the glutamate-cystine antiporter xCT, encoded by the *Slc7a11* gene, also regulates cellular GSH levels¹⁹. We analyzed the regulation of the *Slc7a11* gene by using RT-qPCR (Fig. 4d). There were no differences between Co and Co-dN cells in the control condition, but there was a strong induction by HNE treatment in Co cells (P < 0.0001) but not Co-dN cells.

We then analyzed the cystine uptake using a radioactive tracer ([¹⁴C]-cystine), treating cells with increasing concentrations of HNE (Fig. 4e) or pretreating cells with CCM (10 µM for 24 h) and then treating or not treating
<figure><img src="image_7.png" /><figcaption>**Figure 3.** Exploration of the cellular detoxification capacities of the two cell lines. (a)–(c) Relative expression of Nrf2 target genes (*Gsta4*, *Akr1b8*, and *Aldh2*) in Co and Co-dN cells. The detoxification pathway was induced by treating cells with HNE (40 µM for 6 h). Gene expression was normalized to *Hprt1* expression. (d) Identification of HNE metabolites in both cell lines. [³H]-HNE (40 µM for 30 min) was used as a radioactive tracer. The cell supernatant was analyzed by radio-HPLC. The results are expressed as the percentage of HPLC peaks. Black boxes correspond to Co cells, open boxes to Co-dN cells, and error bars indicate the mean ± standard error of the mean from three independent experiments. Two-way ANOVA followed by Tukey's multiple comparison test was used for statistical analysis (ns = non-significant, \*P < 0.05, \*\*P < 0.01, \*\*\*P < 0.0001).</figcaption></figure>

with HNE (40 µM; Fig. 4f). The results, presented as radioactive cystine incorporated into the cells relative to the total cellular protein content, showed that HNE treatment increased cystine uptake in Co and Co-dN cells, but it was significantly lower in Co-dN cells (P < 0.0001). The decreased cystine uptake after treatment with 80 µM HNE could be due to the cytotoxic effect of the molecule (Fig. 4e). Pretreatment with CCM or combination of CCM pretreatment and subsequent HNE treatment induced cystine uptake in Co cells, with higher uptake in the CCM + HNE condition (P < 0.001; Fig. 4f). The effects were virtually abolished in Co-dN cells. Two-way ANOVA followed by Tukey's multiple comparison test was used for statistical analysis.

## Cell survival under oxidative stress

### Cytotoxicity

To assess the involvement of Nrf2 in cell viability in an oxidative context, we treated both cell lines with increasing concentrations of HNE (from 20 to 80 µM) during proliferation to represent events that take place deep in the colon crypts. There was a dose-dependent decrease in the survival of both cell lines, with Co-dN cells showing lower viability (20, 40 and 80 µM HNE: 78%, 42%, and 8%, respectively, for Co cells, and 48%, 21%, and 1.5%, respectively, for Co-dN cells; Fig. 5a). Going a step further, we pretreated the cells with CCM (10 µM for 24 h) to induce the Nrf2 transcription pathway. Induction of Nrf2 in Co cells resulted in improved survival (20, 40 and 80 µM HNE: 95%, 82%, and 38%, respectively, for Co cells, and 62%, 32%, and 6.7% for Co-dN cells, respectively). We calculated values equivalent to the $IC_{50}$—that is, the HNE concentration that induced 50% of cell death—to represent the toxic effect of HNE under the different conditions. These HNE concentrations were 35 and 20 µM for Co and Co-dN cells, respectively. Induction of Nrf2 by CCM pretreatment protected the cells, with $IC_{50}$-equivalent concentrations of 70 µM for Co cells and 27 µM for Co-dN cells.

We next used the xCELLigence technique. By this approach, we can monitor cell proliferation, real-time cytotoxicity, and survival of Co and Co-dN cells after HNE treatment and/or CCM pretreatment. This tool is complementary to and more accurate than the cell proliferation test. Both cell lines were treated after reaching confluence with an increasing concentration of HNE (from 20 to 80 µM) (Fig. 5b). Positives values of the slope represent growing cells, while a negative slope value reflects a decrease in the number of cells. The decrease in the number of cells was proportional to the HNE concentration in Co cells (P < 0.02 for 40 µM HNE and P < 0.0001
<figure><img src="image_8.png" /><figcaption>**Figure 4.** Assessment of GSH-associated oxidative stress response. (a) The cellular GSH content is represented by the GSH/GSSG ratio. Where indicated, cells were pretreated with CCM (10 µM for 24 h) and/or treated with HNE (40 µM for 24 h) prior to GSH measurement. (b)–(d) Relative expression of Nrf2 target genes (*Gclm*, *Gclc*, and *Slc7a11*) in the two cell lines. The detoxification pathway was induced by treating cells with HNE (40 µM for 6 h). Gene expression was normalized to *Hprt1* expression. (e) and (f) Cystine uptake was assessed after administration of radiolabeled cystine, relative to the total protein content. (g) Cells were treated with HNE (40 or 80 µM for 6 h); (h) where indicated (CCM and CCM+HNE), cells were pretreated with CCM (10 µM for 24 h) and/or treated with HNE (40 µM for 6 h) prior to radioactive measurement. Black boxes represent Co cells, open boxes Co-dN cells, and error bars indicate the mean ± standard error of the mean (n = 3). Two-way ANOVA followed by Tukey's multiple comparison test was used for statistical analysis (ns = non-significant, \*P < 0.05, \*\*P < 0.01, \*\*\*P < 0.001, \*\*\*\*P < 0.0001).</figcaption></figure>

for 80 µM HNE; Fig. 5a). Co-dN cells were more sensitive to the toxic effects of HNE (P < 0.0001 for 40 and 80 µM HNE), with no significant difference between the two concentrations. The difference between the two cell types was highly significant (P < 0.0001) (Fig. 5c). Taken together, CCM alone had no significant effect on cell growth, but it did reduce cellular toxicity induced by 40 µM HNE in Co cells, with no significant difference between CCM+HNE and CCM alone, but a significant difference compared with HNE 40 µM (P < 0.006) (Fig. 5d). CCM did not prevent toxicity in Co-dN cells. Two-way ANOVA followed by Tukey's multiple comparison test was used for statistical analysis.

### Genotoxicity

We determined the genotoxic effect of HNE on both cell lines by measuring the induction of γH2AX. Cells were treated with increasing concentrations of HNE, without reaching cytotoxic doses (i.e., 5 and 10 µM), and the percentage of relative cell count (%RCC) was determined. The genotoxic effect of HNE was greater in Co-dN cells. The induction of γH2AX was 1.5-fold in Co cells (both HNE concentrations) and 1.2-fold at 5 µM HNE and twofold at 10 µM HNE in Co-dN cells (Fig. 6).

### Discussion

The main objective of this work was to generate a suitable cellular model, based on a non-tumoral murine cell line, to explore the role of Nrf2 in the cellular toxicity induced by food contaminants or neo-formed compounds (such as lipoperoxidation products) from digestion. The development of a cellular model free of cancer-related mutations is essential for studying the responses of normal cells to xenobiotic stimuli, including genotoxic agents, and for identifying potential cancer-initiating molecules. As Nrf2 is the main regulator of a panel of antioxidant and cell protection genes, mitigating cellular damage due to oxidative stress, a Nrf2-depleted cell
<figure><img src="image_9.png" /><figcaption>**Figure 5.** Cell survival in an oxidative context. (a) The Co and Co-dN cell lines were treated with increasing concentrations of HNE (0–80 µM), with or without CCM pretreatment (10 µM for 24 h). The results are presented as a percentage of the viability of untreated cells. (b) Cell proliferation measured by the xCELLigence technique. Real-time cellular impedance was measured in each well and the signal was analyzed through integrated software (RTCA Analyzer). Average from 5 independent experiments are represented, blue lines are Co cells, red/pink lines are Co-dN cells, in brackets are HNE concentrations. The green arrow represents the initiation time of treatment by HNE. (c) and (d) Cell proliferation and detachment analyzed using the xCELLigence technique. Increasing concentrations of HNE (0–80 µM) were applied to confluent cells (c). Where indicated, cells were pretreated with CCM (10 µM for 24 h) prior to HNE treatment (40 µM for 6 h) (d). Black boxes or dots are Co cells, open boxes or dots are Co-dN cells, and error bars indicate the mean ± standard error of the mean (n = 6). Two-way ANOVA followed by Tukey's multiple comparison test was used for statistical analysis (circles represent nonspecific differences, ns = non-significant, \*P < 0.05, \*\*P < 0.01, \*\*\*P < 0.0001).</figcaption></figure>

line could allow researchers to assess the implication of Nrf2 in the cellular response to oxidative stress. To our knowledge, we are the first group to generate Nrf2-deficient normal murine epithelial colonic cells. Following CRISPR/Cas9-mediated invalidation of Nrf2, we used HNE as a representative toxic lipid peroxidation product, reflecting the heme-induced cytotoxic and genotoxic alkenals produced in the colon during digestion, to highlight the importance of Nrf2 in cellular protection against toxic injuries⁹. In an attempt to identify the mechanisms involved in these processes, we focused on the GSH pathway, HNE metabolism, target gene modulation, and the cell proliferation capacity. We also used CCM, a well-known inducer of Nrf2¹⁶,¹⁸, to confirm our results.

The Co cell line used in these experiments was derived from the Apc⁺/⁺ clones described by Forest et al.⁸ and characterized by Plaisancié et al.²⁰. To ensure the non-transformed status of our cell line, we performed WGS, focusing specifically on cancer-related genes. Despite a relatively low coverage rate for our sequences, we carried out an in-depth analysis on cancer-related genes to ensure that coverage was sufficient for our analyses. We used the snpEff classification of variants to identify mutations in the tested genes. We compared the Co cell sequences to the C57Bl/6J mouse genome, focusing on genes from different pathways involved in cancer hallmarks. None of these genes were mutated in our cell line. We paid particular attention to the WNT pathway as *Apc* mutations are involved in a majority of human colorectal cancer cases²¹. Moreover, the WNT pathway has been shown to activate Nrf2¹³. We also tested tumor suppressor genes such as Tp53 and proto-oncogenes such as those in the phosphoinositide 3-kinase (PI₃K) pathway. Among all of the tested genes, listed in the Tables S3 and S4, only 3
<figure><img src="image_10.png" /><figcaption>**Figure 6.** The genotoxic effect of HNE was determined by measuring γH2AX induction. γH2AX induction is represented by bars (left Y axis). Cell viability was assessed in parallel (right Y axis), represented by lines. Cells were treated with increasing concentrations of HNE (0–10 µM). Black boxes represent Co cells, open boxes represent Co-dN cells, and error bars indicate the mean ± standard error of the mean (n = 3). Two-way ANOVA followed by Tukey's multiple comparison test was used for statistical analysis (ns = non-significant, \*\*P < 0.01, \*\*\*\*P < 0.0001).</figcaption></figure>

differed from the wild-type. The identified variants associated with Notch4 and Tp53 are non-disruptive variants and the variants with disruptive effect on protein (high impact) were associated with Nrf2, confirming that our cell line is cancer-free²².

The CRISPR/Cas9 technique has rapidly become a reference method for generating mutations in specific genes, leading to gene silencing²³. Comparison in isogenic cell lines of the disruption of a specific gene is particularly interesting for studying specific pathways. After selection, we obtained a Nrf2-depleted clone, which we called Co-dN. Our data confirmed that Nrf2 protein knockdown was successful in these cells. Although Nrf2 mRNA was still partially detectable (Fig. 1a), the protein was not detectable under our conditions, even after HNE stimulation. The Co and Co-dN cell lines are isogenic, and we confirmed the absence of genetic differences between them by performing WGS. Comparing the difference in high impact variants in Co-dN and Co cells, only two genes were impacted, namely Nrf2 and Cdk6. Thus, we confirmed that the CRISPR/Cas9 technique was able to generate a Nrf2-depleted cell line, with no off-target mutations (Fig. 2; Table 1). Only a few cases of functional knockout of the Nrf2 gene by CRISPR/Cas9 have been reported—in human lung cancer cells (A549), human embryonic stem cells²⁴,²⁵, and, recently, in the human colorectal cancer cell line DLD1²⁶, where isogenic wild-type cells, CRISPR/Cas9 generated Nrf2-KO cells and gain-of-function Nrf2 cells were compared. HTC 116 cells, a human colon cancer cell line, have also been used to identify p53 protein activity in drug-induced colon cancer cell death²⁷. Rayner and colleagues²⁸ reported that the CRISPR/Cas9 technique causes chromosomal instabilities and rearrangements in cancer cell lines.

We had previously shown that heme iron, derived from red meat consumption, can promote colorectal cancer, probably due to its strong ability to catalyze the peroxidation of dietary polyunsaturated fatty acids (PUFAs), generating cytotoxic and genotoxic alkenals in the colon²⁹. These molecules have a longer half-life than ROS and are considered second messengers of oxidative stress²,³. Among the reactive aldehydes produced by lipid peroxidation, malondialdehyde (MDA) and HNE are the major products. HNE, the end-product of n-6 PUFA peroxidation, is considered one of the most bioactive aldehydes²,³⁰ and the most cytotoxic to colon cells¹⁴,³¹,³². High levels of HNE in the cytosol can activate Nrf2 via its adduction onto Keap1 cysteines, inducing conformational changes, followed by Nrf2 release and translocation to the nucleus, where its binds to the antioxidant-responsive element (ARE) present in the promoter of several genes³,³¹. Nrf2 has been shown to regulate the expression of antioxidant enzymes, resulting in an adaptive response to alkenal-related injury³,³³,³⁴. Furthermore, HNE has been shown to activate Nrf2 in a variety of cells, including colon cancer cells²,²⁶,³⁵⁻³⁸. In the present study, we used
HNE to decipher the molecular mechanisms involving Nrf2 in cellular protection against cytotoxic molecules resulting from heme-induced lipoperoxidation.

We assessed three of the main HNE detoxification pathways in epithelial cells by analyzing the enzymes Gsta4, Aldh2 and Akr1b8³¹. Gsta4 is the glutathione-S-transferase specific for HNE conjugation to GSH and the main cellular HNE detoxification pathway. It is highly expressed in the neoplastic epithelia of invasive carcinomas and is activated by HNE and Nrf2³⁹. Akr1b8 is responsible for the reduction of reactive aldehydes to their less toxic alcoholic forms⁴⁰. Aldh2 participates in the oxidation of the carbonyl group of HNE leading to the formation of 4-hydroxy-2-nonenoic acid (HNA)⁴¹. The results showed the involvement of Nrf2 in the HNE detoxification process involving Gsta4 and Akr1b8, but not Aldh2. To further investigate the role of Nrf2 in detoxification processes, we treated cells with radiolabeled HNE. In Co-dN cells, the original HNE molecule and its saturated derivatives tended to be more abundant, while HNE-GSH and oxidized metabolites were less present, demonstrating that Nrf2 is the cornerstone of the HNE detoxification process.

GSH is one of the most important components of the cellular antioxidant defense against many reactive compounds. It involves the cystine/glutamate antiporter xCT encoded by *Slc7a11*. Imported cystine is reduced by GCL, via the catalytic and modulatory subunits (Gclc/Gclm), and glutathione synthetase (Gss) to generate GSH³⁴. In addition, HNE is detoxified via conjugation with cysteine to form a non-toxic HNE-cysteine conjugate¹⁴,⁴¹. Tonelli et al.⁴² reported that Nrf2 tightly regulates GSH levels and directly controls the expression of the GCL complex (Gclc/Gclm), playing a fundamental role in maintaining cellular redox homeostasis, in addition to regulating GSH synthesis. We measured the GSH/GSSG ratio; *Gclc*, *Gclc*, and *Slc7a11* gene expression; and cystine uptake in Co and Co-dN cells after HNE treatment and with or without CCM pretreatment. Taken together, the results presented in Figs. 3 and 4 demonstrate that the detoxification process involved in the response of mouse colonic cells to HNE-induced injury involves the GSH pathway, and depletion of Nrf2 results in reduced cell protection against toxic injuries. HNE can induce genes involved in its detoxification process (*Slc7a11* and *Gsta4*)²⁵,³⁴,⁴¹,⁴², and Nrf2 activation may be linked to this process, as it has been suggested that HNE modifies the cysteine residues of Keap1, the negative regulator of Nrf2²⁶,³⁶,⁴¹. Okazaki et al.⁴³ reported that in cancer cells, activation of Nrf2 induces transcriptional activation of *Slc7a11*, *Gclc*, and *Gclm* in response to oxidative and electrophilic stresses, such as HNE. For these authors, cancer cells are addicted to Nrf2 and adopt a highly specialized metabolism favoring GSH synthesis, and therefore have powerful antioxidant detoxification capabilities. These changes in cellular metabolism, combined with activation of the Nrf2 pathway to promote a highly efficient detoxification process in cancer cells, could be linked to a greater resistance to anticancer drugs.

We demonstrated the importance of Nrf2 in cell survival after HNE injury, as increasing concentrations of HNE resulted in the decreased survival of Co and Co-dN cells, with Co-dN cells showing greater sensitivity. Induction of Nrf2 by CCM protected Co cells, as the IC₅₀-equivalent concentration doubled (from 35 to 70 µM). However, CCM had virtually no effect on Co-dN cells. It is important to note that we cannot rule out the possibility that CCM may act on other cellular pathways.

Nrf2 overexpression is known to induce cell proliferation²⁴,⁴⁴,⁴⁵. We used the xCELLigence technique to study the rate of cell proliferation as well as cell viability. Nrf2 depletion led to slower cell growth, confirming the importance of Nrf2 for the proliferation rate in our model. Furthermore, this technique confirmed the toxic effects of HNE (Fig. 5b, c). Pretreatment of cells with CCM prior to HNE treatment protected cells from HNE-induced damage, and had no effect on Co-dN cells (Fig. 5c). These results confirm those obtained in the viability tests (Fig. 5a).

We also showed that Nrf2 could modulate the genotoxic activities of HNE. The genotoxic effects of HNE have already been reported³⁰,³¹,⁴⁶. In our cells subjected to non-cytotoxic concentrations of HNE, there was greater genotoxic damage in Co-dN cells. It could be of great interest to investigate whether HNE, through its Nrf2-activating effects, could be considered a cancer-initiating agent.

We used CCM, a known inducer of Nrf2¹⁶,¹⁸, to confirm the involvement of Nrf2 in cell protection against oxidative stress. This phytochemical has been shown to stabilize Nrf2 protein in mouse epidermal cells¹⁸. This results in activation of Nrf2 and induction of cytoprotective targets genes. Thus, CCM could act as an antioxidant and anti-inflammatory compound, countering ROS formation¹⁵. Considering the “double-edged sword” action of Nrf2⁴⁷, it is interesting to evaluate the effects of Nrf2 stimulation in non-cancerous cells. In Co cells, CCM exerted protective effects, notably by protecting cells from HNE-induced damage leading to cell death. Because of its antitumor and Nrf2-inducing activities, this natural compound has the potential to be used to prevent colorectal cancer¹⁵.

Certain limitations need to be considered when interpreting our results. The strategy for creating a conditionally immortalized murine colon cell line involved introducing the SV40 large T antigen as the immortalization gene, which carries a temperature-sensitive mutation. At the permissive temperature of 33 °C, the transformation of the SV40 large T protein enables it to bind to the p53 protein, inhibiting its role in senescence and thus facilitating cell immortalization. Therefore, cells need to be cultured at 33 °C, which is not the physiological temperature. When cells are cultured at 37 °C (the non-permissive temperature), the SV40 large T protein no longer binds to p53, and the cells can partially differentiate into colonocytes. Furthermore, the HNE concentrations we used are not physiological concentrations found in the colon lumen. However, we used this molecule as a representative of all the toxic alkenals that might be present in the colon lumen. The CRISPR/Cas9 technique led to an incomplete deletion of Nrf2 due to the random insertion/deletion in the genomic sequence. Nevertheless, the Co-dN cells with partial deletion of Nrf2 could be considered a more realistic model given the evidence for the Nrf2 decline with aging.

The cellular model of the non-cancerous Co cell line was enhanced by the generation of an isogenic cell line with Nrf2 depletion, enabling us to highlight the importance of this transcription factor in the cellular response to oxidative stress-induced molecules, and in particular the involvement of GSH in the detoxification of oxidative stress. Oxidative stress is a major feature of aging, in which the ability to respond to oxidative stress is reduced.
Interestingly, during aging, basal Nrf2 protein levels and activity are reduced, as is the expression of genes such as those of the GCL family and Hmox-1<sup>48</sup>. Decreased Nrf2 expression and transcriptional activity may contribute to the loss of homeostasis and affect function that correlate with aging<sup>49-51</sup>. Thus, this dual-cell model represents an original tool for studying the role of Nrf2 in the response to oxidative stress induced by food-derived xenobiotics, luminal peroxidation, and ROS<sup>52</sup>, and could be used to explore metabolic perturbations that emerge during aging.

# Materials and methods

## HNE synthesis

HNE was synthesized as previously described<sup>53</sup>. Briefly, the Grignard reaction between fumaraldehyde mono-acetal and 1-pentylmagnesium bromide afforded HNE-dimethylacetal. HNE was obtained by acid hydrolysis.

## Cell culture conditions

Murine colon epithelial cells were used as described in Surya et al.<sup>32</sup>. The Co cell line used in these experiments was issued from the Apc<sup>+/+</sup> clones described in Forest et al.<sup>8</sup>. This transformed cell line firstly generated by Pierre<sup>8</sup> is present in the lab since then. Cell lines were cultured at permissive temperature of 33 °C in Dulbecco-modified essential medium (DMEM) supplemented with 10% fetal calf serum, 1% penicillin/streptomycin, 2% glutamine, 10 U/ml interferon γ and 10 U/ml epidermal growth factor (EGF) (Merck, France). Excepted when indicated, most experiments were performed at non-permissive temperature of 37 °C without INFγ and EGF to inhibit the SV40 transgene and to reach a differentiated state. HNE treatments were performed in serum-free DMEM containing 2% glutamine for 24 or 6 h as indicated in the experiments. Pretreatment with curcumin (CCM, 10 µM) was performed the day before HNE treatment at 37 °C for 24 h.

## Whole genome sequencing

Genomic DNA for whole-genome sequencing (WGS) was extracted from murine cells using the GenElute Mammalian Genomic DNA Miniprep Kit as described by the manufacturers (Sigma, G1N10). DNAseq was performed at the GeT-PlaGe core facility, INRAe Toulouse. DNA-seq libraries have been prepared according to Illumina's protocols using the Illumina TruSeq Nano DNA HT Library Prep Kit. Briefly, DNA was fragmented by sonication, size selection was performed using SPB beads, and adaptors were ligated to be sequenced. Library quality was assessed using an Advanced Analytical Fragment Analyzer and libraries were quantified by qPCR using the Kapa Library Quantification Kit. DNA-seq experiments have been performed on an Illumina NovaSeq6000 using a paired-end read length of 2×150 pb with the Illumina NovaSeq6000 Reagent Kits. The raw sequences were deposited at SRA under accession PRJNA800346.

## Bioinformatics analysis

FastQC<sup>54</sup> was used to check the quality and no relevant contamination hit was found after the alignment against E. coli, Yeast and PhiX. Reads mapping and variants calling were performed using Nextflow v20.10.0 and Sarek pipeline v2.6.1<sup>55</sup>. Genome Analysis Toolkit (GATK) best practices were followed as implemented in the Sarek pipeline. The following steps were performed: aligning the reads with BWA v0.7.17-r1188<sup>56</sup> against the mouse genome assembly C57BL_6NJ_v1 (GCA_001632555.1), marking duplicate reads (MarkDuplicates), base quality recalibration (BQSR) and calling germline small variants (HaplotypeCaller in GVCF mode) with GATK v4.1.7.0<sup>57</sup> on merged bam by sample. Finally, quality control is performed with MultiQC v1.8<sup>58</sup>. Out of the nextflow pipeline the variants were annotated with VEP(5) version 103.1 with assembly C57BL_6NJ_v1. Variant with at list quality of 30 and depth of 5 were kept to further analyses.

## Nrf2 invalidation

The CRISPR/Cas9 technique used to generate the Nrf2-depleted cell line used the 3 "all-in-one plasmids" containing 3 pre-designed gRNAs (Sigma, Merck Company) that target 3 different sequences at exon 2 of the murine Nrf2 mRNA gene, transcript variant 1, NM_010902.4 (Supplementary Fig S1a). These gRNAs coupled with promoters and the GFP (Green Fluorescent Protein) were inserted into vector plasmids (pCMV-cas9-GFP). These plasmids where amplified and purified (Column Plasmid DNA Mini-preps Kit, BioBasic). Cells were transfected 48 h post plating (30 000 cells/cm2) using X-treme GENE HP DNA Transfection Reagent (Roche) with a ratio of 2 µl XtremeHP/2 µg plasmid DNA. Cells were treated with the three different gRNAs. Transformed cells were analyzed and separated into three distinct GFP positive populations (weak, medium and high) using cell sorting flow cytometry (FACSARIA-II, BD Biosciences, CPTP Toulouse). Only the populations that integrated the largest amount of plasmid (High GFP) were plated at one cell per well of 96-well culture plates rate. DNA mutations in Nrf2 target gene were sought by PCR using primers specific to each gRNA (Supplementary Fig. S1).

## Western blot

Nrf2 expression in whole cell protein extracts was determined by Western blot. Two sets of cells (Co and Co-dN) were treated, or not, by HNE (40 µM for 6 h) before protein extraction. Protein concentration was determined using BCA protein Assay kit (Pierce, Rockford, IL, USA). Proteins were separated using an SDS-PAGE gradient pre-cast gels 4-12%, and transferred onto a nitrocellulose membrane 0.2 µM (Trans-blot turbo mini 0.2 µM nitrocellulose Transfer pack) using Trans-Blot® Turbo™ (7 min, 25v, 1.3 A) (all from Bio-Rad). Detection of the primary antibody (Nrf2 rabbit polyclonal IgG, Genetex Biotech, 1:500) and vinculin (rabbit IgG, Proteintech Europe, 1:5000) was performed with the WesternBreeze system (WesternBreeze™ Chemiluminescent Kit, anti-rabbit for Nrf2 and vinculin, Thermofisher). Protein bands were visualized using Chemidoc with the associated
software (Bio-Rad, ImageLab 5.2.1 software) and the intensity of each bands were evaluated. Full size, images and band intensity are shown in supplementary Fig. S2 and Table S2.

## Gene expression

Total cellular RNA was extracted with Tri reagent (Molecular Research Center). Total RNA samples (1 µg) were then reverse-transcribed with the iScript™ Reverse Transcription Supermix (Bio-Rad) for real-time quantitative polymerase chain reaction (qPCR) analyses. The primers for Sybr Green assays are presented in Supplementary Table S1. Amplifications were performed on a ViiA 7 Real-Time PCR System (Applied Biosystems). The qPCR data were normalized to the level of Hypoxanthine Phosphoribosyl transferase - 1 (Hprt1) messenger RNA and analyzed by LinRegPCR v.11 software.

## GSH/GSSG assay

The GSH/GSSG ratio was evaluated with the GSH/GSSG™ Assay (Promega, France) based in a luminescence detection of GSH-dependent conversion of a GSH probe, Luciferin-NT, to luciferin by a glutathione S-transferase enzyme coupled to a firefly luciferase reaction. Cells were seeded into a 96-well plate. After treatment, total glutathione and oxidized glutathione were measured according to the manufacturer's instructions. Results are expressed as GSH/GSSG ratio.

## Cystine uptake assay

L-[¹⁴C]-Cystine (PerkinElmer, 100 mCi/mmol) uptake was measured according to Giraudi et al. with some modifications⁵⁹. Cells were seeded in a 12-well plate pretreated with CCM and incubated for 6 h with HNE as previously described. Cells were then rinsed with warmed uptake buffer. Cystine uptake was started by incubating the cells in uptake buffer containing L-[¹⁴C]-Cystine 0.1 µCi at room temperature for 10 min, and stopped by rinsing with ice-cold unlabeled uptake buffer. Cells were then lysed by adding NaOH 0.1 N. Lysate was mixed with a scintillation cocktail (UltimaGold, PerkinElmer), protein quantification was done with the BCA protein Assay kit (Pierce, Rockford, IL, USA) and radioactivity was determined using a scintillation counter (Hewlett Packard). Results were expressed as nmol of incorporated cystine reported to total protein content.

## Cellular viability

Cell viability assay was performed by using cell proliferation reagent WST-1 (Roche Life Science) according to manufacturer's protocol. Briefly, cells were seeded and grown in 96-well plates at 33 °C until they reach 50% confluence, then transferred at 37 °C in the non-permissive medium with 10 µM curcumin when indicated (CCM, HC) for 24 h. HNE was then added (0–20–40 and 80 µM) for the next 24 h in a serum-free medium at 37 °C. Cells were incubated in WST-1 reagent for 1 h in the dark. The absorbance was evaluated at 440 nm using a Tecan apparatus. The results from 6 independent experiments were expressed as percentage of non-treated control cells.

## Cellular genotoxicity

Genotoxic effects of increasing concentrations of HNE was measured by γH2AX in-cell Western assay according to Khoury et al.⁶⁰. Cells were seeded into 96-well plates at 15 × 10³ cells per well in permissive conditions. After reached subconfluence, cells were transferred to 37 °C and treated with HNE (0, 5, 10 µM) for 24 h. The results are expressed as γH2AX fold induction. To consider the cytotoxicity of HNE, the percentage of relative cell counts, measured as DNA content had been evaluated and the viability curves were added to the graphs.

## Real-time cellular impedance

The real-time cell impedance analyzer (xCELLigence system) was used according to the manufacturers' instructions (Roche Applied Science, Mannheim, Germany, and ACEA Biosciences, San Diego, CA, USA). Each cell line was added to the 16-well E-plates at exactly 10 000 cells per well. After 96 h at permissive conditions, when cells had reached confluence, the medium was changed to non-permissive medium and temperature, without serum and different HNE concentrations (20; 40; 80 µM). When indicated, cells were pretreated with curcumin (CCM, 10 µM, complete medium, permissive conditions) 24 h before HNE treatments. Real-time cellular impedance was measured in each well (cell index values) and the signal was analyzed through integrated software (RTCA Analyzer). The results are indicated as the slope (h - 1), with positive values representing proliferating cells (increasing cell number), and negatives values represent detachment of cells (decreasing cell number). Results are from 5 independent experiments.

HNE metabolism was studied using radiolabeled HNE synthesized as previously described⁶¹ (specific activity 222 GBq/mol). Cells were treated by 40 µM [³H]-HNE for 30 min, then cellular supernatants were analyzed by radio-HPLC as described before¹⁴, except for the HPLC gradient used that has been modified to improve peak separation: both mobile phases contained water/acetonitrile/acetic acid with phase A (97.5/2.5/0.1) and B (40/60/0.1); the elution gradient was 0 min 18% B, 2–10 min 26% B, 23–37 min 41% B, 42 min 65% B, 47 min 100% B, 60 min 100% B, 65 min 18% B.

## Statistical analysis

All data were expressed as mean ± standard error of the mean (SEM) of three (or more) independent experiments. Statistical significance was determined by two-ways analysis of variance (ANOVA), followed by Tukey's multiple comparison test when needed, except for Fig. 1a and 1c, where a t-test had been performed, using GraphPad Prism 9 software. Statistical significance was indicated by \*P < 0.05, \*\*P < 0.01, \*\*\*P < 0.005 and \*\*\*\*P < 0.001.
## Data availability

All data associated with this study are present in the paper or the supplementary data, and are available from the corresponding author upon reasonable request. The datasets generated (raw sequences) and analyzed during the current study are available in the NCBI/SRA repository under accession PRJNA800346.

Received: 16 December 2022; Accepted: 29 April 2024

Published online: 13 May 2024

## References

1. Harris, I. S. & DeNicola, G. M. The complex Interplay between antioxidants and ROS in cancer. *Trends Cell Biol.* **30**, 440–451 (2020).

2. Gall Troselj, K. et al. Oxidative stress and cancer heterogeneity orchestrate NRF2 roles relevant for therapy response. *Molecules* **27**, 66 (2022).

3. Lee, D.-Y., Song, M.-Y. & Kim, E.-H. Role of oxidative stress and Nrf2/KEAP1 signaling in colorectal cancer: Mechanisms and therapeutic perspectives with phytochemicals. *Antioxidants* **10**, 743 (2021).

4. Mena, S., Ortega, A. & Estrela, J. M. Oxidative stress in environmental-induced carcinogenesis. *Mutat. Res.* **674**, 36–44 (2009).

5. Hayes, J. D. & Dinkova-Kostova, A. T. The Nrf2 regulatory network provides an interface between redox and intermediary metabolism. *Trends Biochem. Sci.* **39**, 199–218 (2014).

6. Yamamoto, M., Kensler, T. W. & Motohashi, H. The KEAP1-NRF2 system: A Thiol-based sensor-effector apparatus for maintaining redox homeostasis. *Physiol. Rev.* **98**, 1169–1203 (2018).

7. Purohit, V., Simeone, D. M. & Lyssiotis, C. A. Metabolic regulation of redox balance in cancer. *Cancers* **11**, 66 (2019).

8. Forest, V., Clement, M., Pierre, F., Meflah, K. & Menanteau, J. Butyrate restores motile function and actin cytoskeletal network integrity in apc mutated mouse colon epithelial cells. *Nutr. Cancer* **45**, 84–92 (2003).

9. Keller, J. et al. Heme-iron-induced production of 4-hydroxynonenal in intestinal lumen may have extra-intestinal consequences through protein-adduct formation. *Antioxidants* **9**, 1293 (2020).

10. McLaren, W. et al. The ensembl variant effect predictor. *Genome Biol.* **17**, 122 (2016).

11. Kopacz, A. et al. Overlooked and valuable facts to know in the NRF2/KEAP1 field. *Free Radic. Biol. Med.* **192**, 37–49 (2022).

12. Peyre, L. et al. Comparative study of bisphenol A and its analogue bisphenol S on human hepatic cells: A focus on their potential involvement in nonalcoholic fatty liver disease. *Food Chem. Toxicol.* **70**, 9–18 (2014).

13. Rada, P. et al. WNT-3A regulates an Axin1/NRF2 complex that regulates antioxidant metabolism in hepatocytes. *Antioxid. Redox Signal.* **22**, 555–571 (2015).

14. Baradat, M. et al. 4-Hydroxy-2(E)-nonenal metabolism differs in Apc+/+ cells and in ApcMin/+ cells: It may explain colon cancer promotion by heme iron. *Chem. Res. Toxicol.* **24**, 1984–1993 (2011).

15. Weng, W. & Goel, A. Curcumin and colorectal cancer: An update and current perspective on this natural medicine. *Semin. Cancer Biol.* https://doi.org/10.1016/j.semcancer.2020.02.011 (2020).

16. Zhao, S.-G. et al. Curcumin attenuates insulin resistance in hepatocytes by inducing Nrf2 nuclear translocation. *Hepatogastroenterology* **58**, 2106–2111 (2011).

17. Patel, S. S. et al. Cellular and molecular mechanisms of curcumin in prevention and treatment of disease. *Crit. Rev. Food Sci. Nutr.* **60**, 887–939 (2020).

18. Shin, J. W. et al. Curcumin induces stabilization of Nrf2 protein through Keap1 cysteine modification. *Biochem. Pharmacol.* **173**, 113820 (2020).

19. Koppula, P., Zhang, Y., Zhuang, L. & Gan, B. Amino acid transporter SLC7A11/xCT at the crossroads of regulating redox homeostasis and nutrient dependency of cancer. *Cancer Commun.* **38**, 12 (2018).

20. Plaisancié, P. et al. Study of the colonic epithelial-mesenchymal dialogue through establishment of two activated or not mesenchymal cell lines: Activated and resting ones differentially modulate colonocytes in co-culture. *PLoS ONE* **17**, e0273858 (2022).

21. Fearon, E. R. Molecular genetics of colorectal cancer. *Annu. Rev. Pathol.* **6**, 479–507 (2011).

22. Hanahan, D. & Weinberg, R. A. The hallmarks of cancer. *Cell* **100**, 57–70 (2000).

23. Giuliano, C. J., Lin, A., Girish, V. & Sheltzer, J. M. Generating single cell-derived knockout clones in mammalian cells with CRISPR/Cas9. *Curr. Protoc. Mol. Biol.* **128**, e100 (2019).

24. Bialk, P., Wang, Y., Banas, K. & Kmiec, E. B. Functional gene knockout of NRF2 increases chemosensitivity of human lung cancer A549 cells in vitro and in a xenograft mouse model. *Mol. Ther. Oncolytics* **11**, 75–89 (2018).

25. Kim, S.-J. et al. Generation of a Nrf2 homozygous knockout human embryonic stem cell line using CRISPR/Cas9. *Stem Cell Res.* **19**, 46–48 (2017).

26. Zhang, Y. et al. Pirin, an Nrf2-regulated protein, is overexpressed in human colorectal tumors. *Antioxidants* **11**, 262 (2022).

27. Garufi, A., Pistritto, G., D'Orazi, V., Cirone, M. & D'Orazi, G. The impact of NRF2 inhibition on drug-induced colon cancer cell death and p53 activity: A pilot study. *Biomolecules* **12**, 461 (2022).

28. Rayner, E. et al. CRISPR-Cas9 causes chromosomal instability and rearrangements in cancer cell lines, detectable by cytogenetic methods. *CRISPR J.* **2**, 406–416 (2019).

29. Guéraud, F. et al. Dietary polyunsaturated fatty acids and heme iron induce oxidative stress biomarkers and a cancer promoting environment in the colon of rats. *Free Radic. Biol. Med.* **83**, 192–200 (2015).

30. Gentile, F. et al. DNA damage by lipid peroxidation products: Implications in cancer, inflammation and autoimmunity. *AIMS Genet.* **04**, 103–137 (2017).

31. Guéraud, F. 4-Hydroxynonenal metabolites and adducts in pre-carcinogenic conditions and cancer. *Free Radic. Biol. Med.* **111**, 196–208 (2017).

32. Surya, R. et al. Red meat and colorectal cancer: Nrf2-dependent antioxidant response contributes to the resistance of preneoplastic colon cells to fecal water of hemoglobin- and beef-fed rats. *Carcinogenesis* **37**, 635–645 (2016).

33. Choi, E.-J. & Jeon, S.-M. NRF2-driven redox metabolism takes center stage in cancer metabolism from an outside-in perspective. *Arch. Pharm. Res.* **43**, 321–336 (2020).

34. Dodson, M., Castro-Portuguez, R. & Zhang, D. D. NRF2 plays a critical role in mitigating lipid peroxidation and ferroptosis. *Redox Biol.* **6**, 101–107. https://doi.org/10.1016/j.redox.2019.101107 (2019).

35. Jung, K.-A. & Kwak, M.-K. Enhanced 4-hydroxynonenal resistance in KEAP1 silenced human colon cancer cells. *Oxid. Med. Cell Longev* **2013**, 423965 (2013).

36. Siow, R. C. M., Ishii, T. & Mann, G. E. Modulation of antioxidant gene expression by 4-hydroxynonenal: Atheroprotective role of the Nrf2/ARE transcription pathway. *Redox Rep.* **12**, 11–15 (2007).

37. Huang, Y., Li, W. & Kong, A.-N.T. Anti-oxidative stress regulator NF-E2-related factor 2 mediates the adaptive induction of antioxidant and detoxifying enzymes by lipid peroxidation metabolite 4-hydroxynonenal. *Cell Biosci.* **2**, 40 (2012).

38. Balogun, E. et al. Curcumin activates the haem oxygenase-1 gene via regulation of Nrf2 and the antioxidant-responsive element. *Biochem. J.* **371**, 887–895 (2003).
39. Yang, Y., Huycke, M. M., Herman, T. S. & Wang, X. Glutathione S-transferase alpha 4 induction by activator protein 1 in colorectal cancer. *Oncogene* **35**, 5795–5806 (2016).

40. Singh, M., Kapoor, A. & Bhatnagar, A. Oxidative and reductive metabolism of lipid-peroxidation derived carbonyls. *Chem. Biol. Interact.* **234**, 261–273 (2015).

41. Dalleau, S., Baradat, M., Guéraud, F. & Huc, L. Cell death and diseases related to oxidative stress: 4-Hydroxynonenal (HNE) in the balance. *Cell Death Differ.* **20**, 1615–1630 (2013).

42. Tonelli, C., Chio, I. I. C. & Tuveson, D. A. Transcriptional regulation by Nrf2. *Antioxid. Redox Signal.* https://doi.org/10.1089/ars.2017.7342 (2017).

43. Okazaki, K., Papagiannakopoulos, T. & Motohashi, H. Metabolic features of cancer cells in NRF2 addiction status. *Biophys. Rev.* https://doi.org/10.1007/s12551-020-00659-8 (2020).

44. Sadeghi, M. R., Jeddi, F., Soozangar, N., Somi, M. H. & Samadi, N. The role of Nrf2-Keap1 axis in colorectal cancer, progression, and chemoresistance. *Tumour Biol.* **39**, 1010428317705510 (2017).

45. Zhang, H.-S. et al. Nrf2 promotes breast cancer cell migration via up-regulation of G6PD/HIF-1α/Notch1 axis. *J. Cell. Mol. Med.* **23**, 3451–3463 (2019).

46. Jaganjac, M. et al. The relevance of pathophysiological alterations in redox signaling of 4-hydroxynonenal for pharmacological therapies of major stress-associated diseases. *Free Radic. Biol. Med.* https://doi.org/10.1016/j.freeradbiomed.2019.11.023 (2019).

47. Gonzalez-Donquiles, C. et al. The NRF2 transcription factor plays a dual role in colorectal cancer: A systematic review. *PLoS ONE* **12**, 66 (2017).

48. Zhang, H., Davies, K. J. A. & Forman, H. J. Oxidative stress response and Nrf2 signaling in aging. *Free Radic. Biol. Med.* **88**, 314–336 (2015).

49. George, M., Tharakan, M., Culberson, J., Reddy, A. P. & Reddy, P. H. Role of Nrf2 in aging, Alzheimer's and other neurodegenerative diseases. *Ageing Res. Rev.* **82**, 101756 (2022).

50. Zinovkin, R. A., Kondratenko, N. D. & Zinovkina, L. A. Does Nrf2 play a role of a master regulator of mammalian aging?. *Biochem. Moscow* **87**, 1465–1476 (2022).

51. Miller, C. J. et al. Disruption of Nrf2/ARE signaling impairs antioxidant mechanisms and promotes cell degradation pathways in aged skeletal muscle. *Biochim Biophys. Acta* **1822**, 1038–1050 (2012).

52. Dickman, C. T. D., Towle, R., Saini, R. & Garnis, C. Molecular characterization of immortalized normal and dysplastic oral cell lines. *J. Oral. Pathol. Med.* **44**, 329–336 (2015).

53. Chandra, A. & Srivastava, S. K. A synthesis of 4-hydroxy-2-trans-nonenal and 4-(3H) 4-hydroxy-2-trans-nonenal. *Lipids* **32**, 779–782 (1997).

54. Babraham Bioinformatics—*FastQC A Quality Control tool for High Throughput Sequence Data*. https://www.bioinformatics.babraham.ac.uk/projects/fastqc/.

55. Garcia, M. et al. Sarek: A portable workflow for whole-genome sequencing analysis of germline and somatic variants. *F1000Res* **9**, 63 (2020).

56. Li, H. & Durbin, R. Fast and accurate long-read alignment with Burrows-Wheeler transform. *Bioinformatics* **26**, 589–595 (2010).

57. McKenna, A. et al. The Genome Analysis Toolkit: A MapReduce framework for analyzing next-generation DNA sequencing data. *Genome Res.* **20**, 1297–1303 (2010).

58. Ewels, P., Magnusson, M., Lundin, S. & Käller, M. MultiQC: Summarize analysis results for multiple tools and samples in a single report. *Bioinformatics* **32**, 3047–3048 (2016).

59. Giraudi, P. J., Bellarosa, C., Coda-Zabetta, C. D., Peruzzo, P. & Tiribelli, C. Functional Induction of the cystine-glutamate exchanger system Xc-activity in SH-SY5Y cells by unconjugated bilirubin. *PLoS ONE* **6**, e29078 (2011).

60. Khoury, L., Zalko, D. & Audebert, M. Validation of high-throughput genotoxicity assay screening using γH2AX in-cell western assay on HepG2 cells. *Environ. Mol. Mutag.* **54**, 737–746 (2013).

61. Alary, J. et al. Mercapturic acid conjugates as urinary end metabolites of the lipid peroxidation product 4-Hydroxy-2-nonenal in the rat. *Chem. Res. Toxicol.* **8**, 34–39 (1995).

## Acknowledgements

This project has received funding from La Ligue Contre le Cancer, Midi-Pyrénées.

## Author contributions

C.H.-T. and J.D. created and characterized the reported cell lines. C.B. performed the western blot analyses, E.F. performed the RTqPCR experiments, C.N. and P.M. performed the sequences analyses and comparisons, C.H.-T. wrote the manuscript. All authors and specially F.G. and F.P. reviewed the manuscript.

## Competing interests

The authors declare no competing interests.

## Additional information

**Supplementary Information** The online version contains supplementary material available at https://doi.org/10.1038/s41598-024-60938-2.

**Correspondence** and requests for materials should be addressed to C.H.-T.

**Reprints and permissions information** is available at www.nature.com/reprints.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
<figure><img src="image_11.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International

License, which permits use, sharing, adaptation, distribution and reproduction in any medium or
format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the
Creative Commons licence, and indicate if changes were made. The images or other third party material in this
article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the
material. If material is not included in the article's Creative Commons licence and your intended use is not
permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from
the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2024