# Exploring the significance of the Exon 4-skipping isoform of the ZNF217 oncogene in breast cancer

Aurélie Bellanger, Diep Le, Julie A Vendrell, Anne Wierinckx, Lórinc Pongor, Jérôme Solassol, Joël Lachuer, Philippe Clezardin, Balázs Győrffy, Pascale Cohen

► To cite this version:

Aurélie Bellanger, Diep Le, Julie A Vendrell, Anne Wierinckx, Lórinc Pongor, et al. Exploring the significance of the Exon 4-skipping isoform of the ZNF217 oncogene in breast cancer. Frontiers in Oncology, 2021, 11, pp.647269. 10.3389/fonc.2021.647269. hal-03653025

HAL Id: hal-03653025

https://hal.science/hal-03653025

Submitted on 27 Apr 2022

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

# Exploring the Significance of the Exon 4-Skipping Isoform of the ZNF217 Oncogene in Breast Cancer

Aurélie Bellanger¹,²†, Diep T. Le¹,³†, Julie Vendrell⁴, Anne Wierinckx¹,²,⁵,
Lőrinc S. Pongor⁶,⁷, Jérôme Solassol⁴,⁸, Joël Lachuer¹,²,⁵, Philippe Clezardin¹,³,
Balázs Győrffy⁶,⁷ and Pascale A. Cohen¹,²,³,⁵*

**OPEN ACCESS**

**Edited by:**
Olivier Cuvillier,
UPR8241 Laboratoire de Chimie de
Coordination (LCC), France

**Reviewed by:**
Qing Cong,
Fudan University, China
Paul Yaswen,
Lawrence Berkeley National
Laboratory, United States

**\*Correspondence:**
Pascale A. Cohen
pascale.cohen@univ-lyon1.fr

†These authors have contributed
equally to this work

**Specialty section:**
This article was submitted to
Pharmacology of
Anti-Cancer Drugs,
a section of the journal
Frontiers in Oncology

**Received:** 29 January 2021

**Accepted:** 31 May 2021

**Published:** 02 July 2021

**Citation:**
Bellanger A, Le DT, Vendrell J, Wierinckx A, Pongor LS, Solassol J, Lachuer J, Clezardin P, Győrffy B and Cohen PA (2021) Exploring the Significance of the Exon 4-Skipping Isoform of the ZNF217 Oncogene in Breast Cancer. *Front. Oncol.* 11:647269.
doi: 10.3389/fonc.2021.647269

¹ Université Lyon 1, Lyon, France, ² CRCL-Centre de Recherche en Cancérologie de Lyon-Inserm U1052-CNRS U5286, Lyon, France, ³ INSERM, UMR1033 LYOS, Lyon, France, ⁴ Département de Pathologie et Oncobiologie, Laboratoire de Biologie des Tumeurs Solides, CHU Montpellier, Univ. Montpellier, Montpellier, France, ⁵ ProfileXpert, SFR-Est, CNRS UMR-S3453, INSERM US7, Lyon, France, ⁶ Department of Bioinformatics, Semmelweis University, Budapest, Hungary, ⁷ TTK "Természet Tudományi Kutatóközpont" Momentum Cancer Biomarker Research Group, Institute of Enzymology, Budapest, Hungary, ⁸ Institut de Recherche en Cancérologie de Montpellier (IRCM), INSERM U1194, Univ. Montpellier, Montpellier, France

Oncogene alternative splicing events can create distinct functional transcripts that offer new candidate prognostic biomarkers for breast cancer. ZNF217 is a well-established oncogene but its exon 4-skipping isoform (ZNF217-ΔE4) has never been investigated in terms of clinical or biological relevance. Using *in silico* RNA-seq and RT-qPCR analyses, we demonstrated for the first time the existence of ZNF217-ΔE4 transcripts in primary breast tumors, and a positive correlation between ZNF217-ΔE4 mRNA levels and those of the wild-type oncogene (ZNF217-WT). A pilot retrospective analysis revealed that, in the Luminal subclass, the combination of the two ZNF217 variants (the ZNF217-ΔE4-WT gene-expression signature) provided more information than the mRNA expression levels of each isoform alone. Ectopic overexpression of ZNF217-ΔE4 in breast cancer cells promoted an aggressive phenotype and an increase in ZNF217-WT expression levels that was inversely correlated with DNA methylation of the ZNF217 gene. This study provides new insights into the possible role of the ZNF217-ΔE4 splice variant in breast cancer and suggests a close interplay between the ZNF217-WT and ZNF217-ΔE4 isoforms. Our data suggest that a dual signature combining the expression levels of these two isoforms may serve as a novel prognostic biomarker allowing better stratification of breast cancers with good prognosis and aiding clinicians in therapeutic decisions.

**Keywords:** ZNF217, isoform, splice variant, prognosis, breast cancer

**Abbreviations:** ZNF217-WT, ZNF217 wild-type; ZNF217-ΔE4, exon 4-skipping ZNF217 isoform; E3, exon 3; E3-E4, exon 3 -
exon 4 junction; E3-E5, exon 3 - exon 5 junction; ER, Estrogen Receptor alpha; ER-, ER-negative; HER2, Human epidermal
growth factor receptor 2; PR, Progesterone receptor; RT-qPCR, Real-Time quantitative Polymerase Chain Reaction; RFS,
relapse free survival; ASE, alternative splicing events; KLF, Krüppel-like factor; AA, amino acids; SBR grade, Scarff-Bloom-
Richardson grade.
# INTRODUCTION

Breast cancer is the most frequent cancer in women, with an estimated 2.2 million new cases being diagnosed in 2020 worldwide, and the leading cause of cancer death in women (> 684,000 cases in 2020) (1). The ZNF217 gene is located on chromosome 20q13.2, a region frequently amplified in many tumors, including those of the breast (2). The encoded protein is an oncogenic transcription factor belonging to the Krüppel-like zinc finger transcription factors family. Previous studies highlighted that the ZNF217 oncogene is involved in both early and late stages of tumor progression [for review (3)]. Our previous work contributed to the demonstration that breast cancer cells possessing high ZNF217 expression levels display a more aggressive phenotype (e.g. increased cell proliferation, increased invasive properties and resistance to chemotherapy) (4, 5). High ZNF217 mRNA levels in primary breast tumors are also of bad prognosis and associated with shorter relapse free survival (RFS) and metastasis development (5, 6), with the most discriminatory prognostic power observed in Luminal breast cancers (7).

Alternative splicing is now considered as a hallmark of cancer (8), participating in tumor progression through the expression of distinct isoforms. Key splicing events targeting important factors in tumorigenesis have been identified in breast cancer (e.g., TP53, HER2, BRCA1, FGFR, BIRC5/survivin,). Alternative splicing events (ASE) in breast cancer thus represent candidate biomarkers and candidate therapeutic targets [for review (9)]. Considering the ZNF217 oncogene, only little is known about the existence of its splice variants. In the original study discovering the human ZNF217 gene (2), the *in silico* alignment of cDNAs from colon carcinoma and HeLa cell lines suggested a possible alternative processing of the 133-base pairs exon 4. The corresponding putative exon 4-skipping ZNF217 (ZNF217-ΔE4) mRNA would encode for a ZNF217 protein isoform with a C-terminus sequence [1013-1061 amino acids (aa)] distinct from that of ZNF217 wild-type (ZNF217-WT) protein (1013-1048 aa sequence) (2). Nevertheless, since the original description (2), the ZNF217-ΔE4 variant has never been studied, either in terms of expression in human tumor samples or in terms of biological relevance.

Considering the key role of the ZNF217 oncogene in breast cancer and its strong biomarker value, the ZNF217-ΔE4 transcript is thus of utmost interest. The main objectives of the present work were to newly investigate: (i) whether ZNF217-ΔE4 is expressed in primary breast cancers; (ii) if so, whether ZNF217-ΔE4 mRNA levels possess any clinical relevance; and (iii) whether increased expression levels of ZNF217-ΔE4 functionally affect the phenotype of breast cancer cells.

By performing *in silico* RNAseq and real-time quantitative polymerase chain reaction (RT-qPCR) investigations, our study highlights the existence of ZNF217-ΔE4 transcripts in primary breast tumors. Interestingly, ZNF217-ΔE4 mRNA levels were significantly correlated with those of ZNF217-WT in tumor samples. A pilot retrospective analysis revealed that, in the Luminal subclass, the combination of the two ZNF217 variants (the ZNF217-ΔE4-WT gene-expression signature) provided

more information than the mRNA expression levels of each isoform alone. Strikingly, the ectopic overexpression of ZNF217-ΔE4 in breast cancer cells led to: (i) a more aggressive phenotype, similar to that previously observed with ZNF217-WT (4, 5); (ii) increased expression of endogenous ZNF217-WT; and (iii) decreased DNA methylation on key CpG sites within the ZNF217 gene.

# MATERIALS AND METHODS

## Analysis of RNA-seq Data

A total of 1,097 RNA-seq aligned BAM files were obtained from the TCGA repository (10). Samples were filtered to include only primary tumor samples in women. The ratio of the exon 4 skipping isoform of the ZNF217 gene in each sample was calculated based on the number of reads that had spliced alignments from exon 3 to exon 4 junction (E3-E4) compared to all reads mapped to the exon 3 to exon 5 junction (E3-E5). The ratio was used to quantify the isoform expression from the DESeq2 normalized data for each sample separately. Expression for exon 3 was computed using the following formula: $\text{Exon}_3\_DESeq = (\text{Exon}_3\_reads/\text{Total_reads})*\text{Gene_DESeq}$. A minimal expression of 1,000 was used as a cutoff to designate patients as those positive for the ZNF217 exon 4-skipping isoform. Kruskal-Wallis test was used to compare breast cancer subclasses distribution. The correlation of expression was tested using the Pearson correlation coefficient. Statistical significance was set at *p* < 0.05.

## ZNF217 Isoforms Plasmid Constructs

The pcDNA6-ZNF217-WT plasmid has been previously described and established (4). The pcDNA6 plasmid (Invitrogen, Cergy Pontoise, Paris) containing ZNF217 exon4 skipping isoform DNA sequence (pcDNA6-ZNF217-ΔE4) was produced as follows: the DNA sequence specific for ZNF217-ΔE4 was obtained after polymerase chain reaction (PCR) amplification from breast cancer cells using specific primers, then ligated into the restriction enzymes pre-digested pcDNA6-ZNF217-WT plasmid to replace the DNA region specific for ZNF217-WT with that of ZNF217-ΔE4.

## RT-qPCR Amplification of ZNF217 Isoforms

Total RNA extraction, reverse transcription, and RT-qPCR measurements were performed as described previously (4, 5). CFX96 equipment with the SsoAdvanced Universal SYBR green supermix (BioRad, Hercules, USA) was designed for RT-qPCR measurements, according to the manufacturer's recommendations. Specific primers were used to amplify the DNA sequence for ZNF217 isoforms. The 5'-AGTCCAAATCCCTGCCATCT-3' and 5'-GGGGAAACACTGGTTTTAGG-3' primers amplify a region within exon 3 (E3). The 5'-CTCGACGTTAGAAGGAAAAAAG-3' and 5'-TGGTCGATAATGTGCATTCC-3' primers were used to explore ZNF217-WT (exon 3 - exon 4 junction, E3-E4). The 5'-GTGGCTGACTGTTCAGAAGCCC-3' and 5'-GACATCCA CCAAGACCTTCTA-3' primers are specific for the ZNF217-ΔE4
isoform (exon 3 – exon 5 junction, E3-E5). For investigation in breast tumor samples and MDA-MB-231 stable transfectants, all measurements were normalized to the ribosomal 28S gene expression using the 5'-CGATCCATCATCCGCAATG-3' and 5'-AGCCAAGCTCAGCGCAAC-3' primers.

## Primary Breast Tumor Cohort

Women with primary breast tumors (n = 107) and known clinical follow-up who had not received any therapy before surgery and who relapsed, or not, while receiving endocrine therapy and/or chemotherapy were recruited from the BB-0033-00050, Biological Resources Center (CRB) Centre Léon Bérard, Lyon France (Supplementary Table 1). This study has been approved by the local ethics committee (CRB Centre Léon Bérard, France). The CRB Centre Léon Bérard is quality certified according to NFS96-900 French standard and ISO 9001 for clinical trials, ensuring scientific rigor for sample conservation, traceability, and quality, as well as ethical rules observance and defined rules for transferring samples for research purposes (Ministry of Health for activities authorization no AC-2019-3426 and DC-2008-99). The material used in the study has been collected in agreement with all applicable laws, rules, and requests of French and European government authorities, including the patients' informed written consents. Extraction of total RNA from frozen tumor samples and RT-qPCR measurements were performed as previously described (5). Univariate (log-rank) analysis, multivariate analysis and all statistical analyses were performed using the SPSS™ Software (IBM, USA). For univariate analyses, the data were divided at the median value of ZNF217-E3, ZNF217-WT (E3-E4), or ZNF217-ΔE4 (E3-E5) mRNA expression into two groups with either high or low expression levels. A *p* < 0.05 was considered statistically significant.

## Establishment of MDA-MB-231-ZNF217-ΔE4 Stable Transfectants

MDA-MB-231 breast cancer cells were grown according to recommendations in DMEM medium supplemented with 10% fetal bovine serum (Invitrogen). MDA-MB-231 cells were stably transfected with pcDNA6 or pcDNA6-ZNF217-ΔE4 plasmids, then selected in the presence of 20 µg/ml blasticidin. MDA-MB-231-ZNF217-WT cells have been previously established (4, 5).

## Western Blotting

Western blot experiments were performed as previously described (5) with a commercial anti-ZNF217 antibody produced against a peptide comprised within the 1000-1048 aa sequence of ZNF217-WT (#48133; Abcam, Paris, France), a home-made rabbit polyclonal antibody directed against the C-terminus 1045-1060 aa sequence specific for ZNF217-ΔE4 isoform (RM217 Ab, Covalab, Lyon, France), the home-made polyclonal Cov-2 antibody directed against the C-terminus 1035-1048 aa sequence specific for ZNF217-WT isoform (Covalab) and anti-α-tubulin antibody (#T5168; Sigma, St Louis, MO, USA).

## Cell Proliferation

Cells (20,000 cells per well) were plated in triplicates into a 24-well plate then grown up to 96h. The medium was changed

every 2 days. Proliferating cells were analyzed using the Scepter™ 2.0 Cell Counter (Merck Millipore, Billerica, USA).

## Cytotoxicity Assay

As previously described (4), cells (8,000 cells per well) were plated onto a 96-well plate, treated for 4 days with 10⁻¹² to 10⁻⁶ M of paclitaxel (Paxene®, Ivax, Miami, USA). Cell viability was assessed with the CellTiter 96 AQueous One Solution Cell Proliferation assay (Promega, Madison, WI, USA).

## Soft-Agar Colony-Formation Assay

The experiment was performed as previously described using 10,000 cells (11). Fifteen days later, the cells were stained with 0.005% Cristal violet (Sigma-Aldrich) for 1 h.

## Methylome

For each cell line, genomic DNA was extracted from two independent cell cultures using the QIAamp DNA mini kit (Qiagen, Hilden, Germany). The Illumina Infinium MethylationEPIC Beadchip was used to obtain DNA methylation profiles according to the manufacturer's instructions. Briefly, 500 ng of gDNA (QuantiFluor® dsDNA System, Promega) were bisulfite-modified using EZ-DNA Methylation-Gold Kit (Zymo Research) and then used for the genome-wide methylation assays. The iDat's files were obtained by scanning the array on iScan scanner (Illumina).

## Pyrosequencing Analysis

Genomic DNA (500 ng) was used in bisulfite conversion reactions with the EpiTect® Bisulfite Kit (Qiagen), according to the manufacturer's protocol. The bisulfite converted DNA was amplified using the PyroMark PCR kit (Qiagen) using the 5'-GTTATTTGTTAAGAAGTGAAGGAAT-3' and the 5'-biotin-CCTTACTACTAAATAACTTAAAATC-3' primers for the two following CpG: chr20:52198378 (cg01692482) and chr20:52198279 (cg00703481). PCR products were used for CpG quantification using a PyroMark Q24 (Qiagen) as the supplier's recommendation. The percentage of methylation of the two CpG sites was calculated using the PyroMark Q24 Advanced software (Qiagen).

# RESULTS

## In Silico RNA-seq Analysis Reveals the Expression of ZNF217 Exon 4-Skipping Isoform in Primary Breast Tumor Samples

RNA-seq files obtained from the TCGA repository (10) were used to decipher whether the ZNF217-ΔE4 isoform is expressed in 1,097 primary breast tumors. Identification of reads mapped to the exon 3 to exon 4 junction (E3-E4) of the ZNF217 gene corresponds to the expression of ZNF217-WT isoform in breast tumor samples. Our analysis further discovered spliced alignments from exon 3 to exon 5 (exon 4-skipping, E3-E5), revealing the expression of the ZNF217-ΔE4 isoform in primary breast cancers. The distribution of ZNF217-WT (E3-E4) and
ZNF217-ΔE4 (E3-E5) mRNA levels in the primary breast tumors highlighted that the ZNF217-ΔE4 isoform mRNA expression levels are globally weaker than those of the ZNF217-WT isoform (Figures 1A, B). ZNF217-ΔE4 expression levels ranged from 0 to

6,000, with 42.1% of the patients expressing the ZNF217-ΔE4 isoform over 1,000, while ZNF217-WT expression levels ranged from 0 to 27,000 (Figures 1A, B). The ZNF217-ΔE4 (E3-E5) expression levels represented more than 10% of ZNF217-WT

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

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

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

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

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

**FIGURE 1** | RNAseq expression of ZNF217 isoforms in breast cancer. (A) ZNF217-WT (E3-E4) expression and (B) ZNF217-ΔE4 (E3-E5) expression, in bins of thousands. (C) A high correlation was observed when correlating the expression of splice reads mapped to exon 3 – exon 5 junction versus those mapped to exon 3 – exon 4 junction. There is a significant expression difference between breast cancer molecular subtypes of both (D) exon 3 –exon 4 reads and (E) exon 3 - exon 5 reads.
expression levels in 76% of the patients (between 10% and 20% in 61.5% of the patients, and higher than 20% in 14.5% of the patients (data not shown)).

Importantly, **Figure 1C** reveals that ZNF217-WT (E3-E4) and ZNF217-ΔE4 (E3-E5) mRNA levels are significantly positively correlated (Pearson coefficient = 0.891, *p* < 10<sup>-16</sup>). The 1,097 patients were then stratified according to molecular subtypes into Luminal, Triple negative, HER2-enriched subclasses (12–14). **Figures 1D, E** show that both ZNF217-WT (E3-E4) and ZNF217-ΔE4 (E3-E5) display significant expression differences between the molecular subtypes (*p* = 2.13 10<sup>-11</sup> and *p* = 3.04 10<sup>-4</sup>, respectively, Kruskal-Wallis test).

Altogether, the RNA-seq data analysis newly demonstrated: (i) the existence of the exon 4-skipping variant of the ZNF217 oncogene in primary breast cancers, expressed at weaker mRNA levels than those of ZNF217-WT; (ii) correlated mRNA expression levels of ZNF217-ΔE4 and ZNF217-WT isoforms in primary breast tumor.

## Establishment of Primers for Specific RT-qPCR Investigation of ZNF217 Isoforms

With the aim of further specifically investigating ZNF217-ΔE4 (E3-E5) mRNA expression in clinical samples, three specific primer pairs were designed to discriminate the ZNF217-WT (E3-E4) and ZNF217-ΔE4 (E3-E5) mRNA expression levels (Supplementary Figure 1A). The primer pair 1 amplifies within the sequence of ZNF217's exon 3, a region common to both ZNF217-WT and ZNF217-ΔE4 isoforms. The isoforms thus detected with primer pair 1 will be named ZNF217-E3. The primer pair 2 detects the ZNF217-WT isoform since the forward primer is located at the E3-E4 junction. The primer pairs 3 detects ZNF217-ΔE4 isoform, with the reverse primer hybridizing to the E3-E5 junction. The PCR efficiency calculated from standard curves validated that the three primer pairs led to similar and high PCR efficiency values (data not shown).

RT-qPCRs were performed to assess the three primer pairs' specificity within a complex DNA mixture containing different proportions of pcDNA6-ZNF217-WT and/or pcDNA6-ZNF217-ΔE4 plasmids. **Supplementary Figure 1B** shows a steady number of ZNF217-E3 copies corresponding to a constant number of pcDNA6-ZNF217 plasmids copies in the mixture. Importantly, primer pair 2 or primer pair 3 were able to discriminate and quantify the number of, respectively, pcDNA6-ZNF217-WT or pcDNA6-ZNF217-ΔE4 copies actually deposited in the mixture (**Supplementary Figure 1B**). Altogether, these results validated that these pairs of primers are specific and sensitive for each isoform, even within a mixture of ZNF217 isoforms.

## RT-qPCR Validation of ZNF217-WT and ZNF217-ΔE4 Expression in an Independent Cohort of Primary Breast Tumors

We performed RT-qPCR analyses to explore ZNF217-WT (E3-E4) and ZNF217-ΔE4 (E3-E5) mRNA expression levels in a cohort of 107 human primary breast tumor samples (Supplementary

**Table 1**). The patterns of distribution of ZNF217-WT (E3-E4) and ZNF217-ΔE4 (E3-E5) mRNA levels over the 107 primary breast tumors (Supplementary Figures 1C, D) were very similar to those observed in Figures 1A, B. Overall, the ZNF217-ΔE4 (E3-E5) mRNA levels were 6-8 fold weaker than those of ZNF217-WT (E3-E4) within a particular primary breast tumor (data not shown). The ZNF217-ΔE4 (E3-E5) mRNA levels represented less than 10% of ZNF217-WT expression levels in 25% of the patients, 10-30% in 44% of the patients, and higher than 30% in 31% of the patients (data not shown). In accordance with what we have observed by RNAseq analyses, the RT-qPCR investigation found ZNF217-WT (E3-E4) and ZNF217-ΔE4 (E3-E5) mRNA levels significantly and positively correlated (Pearson coefficient = 0.807, *p* < 10<sup>-4</sup>). Again, both ZNF217-WT (E3-E4) and ZNF217-ΔE4 (E3-E5) isoforms display significant mRNA expression differences between the breast cancer molecular subtypes (*p* = 0.007 and *p* = 0.056, respectively, Kruskal-Wallis Test, data not shown).

## The ZNF217-ΔE4-WT Gene Expression Signature Indicates Poor Prognosis in Luminal Breast Cancers

We have previously reported that high mRNA levels of ZNF217 are of poor prognosis and represent an early indicator of relapse in breast cancer, with the most powerful prognostic value observed in Luminal subtypes (5, 7). However, at the time when those studies were conducted, the possible expression of the ZNF217-ΔE4 isoform in breast tumors was unknown. The primer pair used in our previous studies (5, 7) targets exon 3 of the ZNF217 gene (primer pair 1, ZNF217-E3), thus exploring in reality both ZNF217-WT and ZNF217-ΔE4 mRNA expression levels. Taking advantage of the specificity of primer pair 2 and primer pair 3, we thus aimed at investigating in a pilot study the prognostic value of ZNF217-WT (E3-E4) and ZNF217-ΔE4 (E3-E5) mRNA expression levels, alone or in combination, in the molecular subtypes of breast cancer.

In accordance with our previous work (5, 7), the present retrospective analysis of the cohort composed of 107 human primary breast tumor samples confirmed that high ZNF217-E3 mRNA levels: (i) were significantly associated with shorter relapse free survival (RFS) in the whole cohort (*p* = 0.031 **Table 1**) and in the Luminal subclass (*p* = 0.030, **Table 1** and Supplementary Figure 2A); (ii) were not informative in the HER2-enriched, or Triple negative subclasses (**Table 1**). In the whole cohort, ZNF217-ΔE4 (E3-E5) mRNA levels were significantly associated with shorter RFS (*p* = 0.025, **Table 1**). In the Luminal subclass, high mRNA levels of ZNF217-WT (E3-E4) and high mRNA levels of ZNF217-ΔE4 (E3-E5) did not reach significance after 10 years of retrospective analysis, while a tendency could be observed (**Table 1** and Supplementary Figures 2B, C). We then constructed a gene expression signature based on ZNF217-WT (E3-E4) and ZNF217-ΔE4 (E3-E5) mRNA expression levels by dividing patients into three groups: patients who expressed low ZNF217-WT (E3-E4) and low ZNF217-ΔE4 (E3-E5) mRNA levels (group 1); patients who expressed high ZNF217-WT (E3-E4) or high ZNF217-ΔE4
<table><caption>**TABLE 1** | Univariate analyses (10 years-retrospective analysis) of the ZNF217 isoforms mRNA expression levels with regards to relapse-free survival.</caption><thead><tr><th rowspan="2"></th><th colspan="3">All tumor samples<br>(n = 107)</th><th colspan="3">Luminal subclass<sup>d</sup><br>(ER+ and/or PR+)<br>(n = 72)</th><th colspan="3">HER2-enriched subclass<sup>d</sup><br>(ER-/PR-/HER2+)<br>(n = 15)</th><th colspan="3">Triple negative subclass<sup>d</sup><br>(ER-/PR-/HER2-)<br>(n = 20)</th></tr><tr><th>HR<sup>a</sup></th><th>95% CI<sup>b</sup></th><th>p<sup>c</sup></th><th>HR<sup>a</sup></th><th>95% CI<sup>b</sup></th><th>p<sup>c</sup></th><th>HR<sup>a</sup></th><th>95% CI<sup>b</sup></th><th>p<sup>c</sup></th><th>HR<sup>a</sup></th><th>95% CI<sup>b</sup></th><th>p<sup>c</sup></th></tr></thead><tbody><tr><td><strong>ZNF217-E3</strong><br>mRNA levels</td><td>4.62</td><td>0.85-5.52</td><td>0.031</td><td>4.72</td><td>1.04-13.20</td><td>0.030</td><td>0.026</td><td>0.17-4.57</td><td>NS<br>(0.872)</td><td>0.024</td><td>0.05-12.87</td><td>NS<br>(0.876)</td></tr><tr><td><strong>ZNF217-WT (E3-E4)</strong><br>mRNA levels</td><td>3.32</td><td>0.76-4.56</td><td>NS<br>(0.068)</td><td>2.86</td><td>0.82-8.23</td><td>NS<br>(0.091)</td><td>1.01</td><td>0.41-13.63</td><td>NS<br>(0.314)</td><td>0.001</td><td>0.06-15.18</td><td>NS<br>(0.970)</td></tr><tr><td><strong>ZNF217-ΔE4 (E3-E5)</strong><br>mRNA levels</td><td>5.05</td><td>0.93-6.05</td><td>0.025</td><td>3.00</td><td>0.84-8.33</td><td>NS<br>(0.083)</td><td>0.216</td><td>0.27-8.21</td><td>NS<br>(0.642)</td><td>0.161</td><td>N/A<sup>e</sup></td><td>NS<br>(0.204)</td></tr><tr><td><strong>ZNF217-WT-ΔE4<br>signature</strong></td><td>5.12</td><td>0.95-8.29</td><td>0.024</td><td>5.11</td><td>1.07-21.23</td><td>0.023</td><td>1.54</td><td>0.41-30.84</td><td>NS<br>(0.214)</td><td>1.296</td><td>N/A<sup>e</sup></td><td>NS<br>(0.255)</td></tr></tbody></table>

(E3-E5) mRNA levels (group 2); patients who expressed high ZNF217-WT (E3-E4) and high ZNF217-ΔE4 (E3-E5) mRNA levels (group 3). **Supplementary Figure 2D** illustrate that: (i) 83.3% (60 out of 72) of the patients in the Luminal subclass belonged either to group 1 or to group 3, again supporting correlated mRNA levels of ZNF217-WT and ZNF217-ΔE4 in human primary breast tumors; (ii) the resulting Kaplan Meier curves for RFS in Luminal subclasses showed patients' stratification into diverse prognostic classes (p = 0.05), with the patients belonging to group 1 having the better prognosis (and Kaplan Meier curves for RFS for group 2 and group 3 close). Consequently, we further constructed a ZNF217-WT-ΔE4 gene expression signature based on two new groups: patients who expressed low ZNF217-WT (E3-E4) and low ZNF217-ΔE4 (E3-E5) mRNA levels; patients who expressed high ZNF217-WT (E3-E4) and/or high ZNF217-ΔE4 (E3-E5) mRNA levels. Strikingly, the prognostic value of the ZNF217-WT-ΔE4 gene expression signature for RFS was demonstrated to be significant in the Luminal subclasses (p = 0.023, **Table 1** and **Supplementary Figure 2E**), and not informative in the HER2-enriched, or Triple negative subclasses (**Table 1**).

To investigate the clinical relevance of our findings, we also investigated by univariate analysis the classical prognostic parameters usually studied in breast cancer. **Supplementary Table 2** and **Table 1** indicated that in the whole cohort, the prognostic value of ZNF217-E3 or of ZNF217-WT-ΔE4 gene expression signature for RFS was more informative than conventional biomarkers (age, SBR grade, Lymph node status, Macroscopic tumor size, ER status, PR status, HER2 status). In the Luminal subclass (ER+ and or PR+), the prognostic value of ZNF217-E3 mRNA levels or of ZNF217-WT-ΔE4 gene expression signature for RFS was more informative than age, SBR grade, Lymph node status, Macroscopic tumor size, and HER2 status (**Table 1** and **Supplementary Table 2**). However, ER status and PR status were associated with shorter RFS (p = 0.025 and p = 0.032, respectively, **Supplementary Table 2**). Prognostic factors for RFS in the Luminal subclass with a 0.05 significance level in univariate analysis were then entered in a multivariate Cox model. When entering ZNF217-E3 mRNA

levels, ER status, PR status in the multivariate analysis, only ZNF217-E3 mRNA levels persisted in the model (p = 0.04, HR = 3.71; 95%CI = 1.04-13.18, data not shown). Similarly, when entering ZNF217-WT-ΔE4 gene expression signature, ER status, PR status in the multivariate analysis, only ZNF217-WT-ΔE4 gene expression signature persisted in the model (p = 0.04, HR = 4.78; 95%CI = 1.07-21.23, data not shown). Finally, no association was found between ZNF217-E3, ZNF217-WT, ZNF217-ΔE4 or ZNF217-WT-ΔE4 gene expression signature with age, invaded lymph nodes, macroscopic tumor size, ER status, PR status and HER2 status (**Supplementary Table 3**). Considering SBR grade, while the p value was significant (p = 0.002), the distribution of ZNF217-WT mRNA levels did not reflect any association between high ZNF217-WT mRNA levels and high grade, or vice versa (**Supplementary Table 3**). With the aim to clarify the association between ZNF217-WT-ΔE4 gene expression signature and SBR grade (p = 0.024, **Supplementary Table 3**), we performed univariate analyses in both SBRI+II (n= 45) and SBRIII (n= 61) subclasses. Strikingly, ZNF217-WT-ΔE4 gene expression signature was significantly associated with shorter RFS in the SBRI+II subclass (p = 0.005, data not shown), but not in the SBRIII subclass (p = 0.64, data not shown), indicating that the ZNF217-WT-ΔE4 gene expression signature could allow the re-stratification of patients with breast cancers of low/intermediate grade.

Altogether, our pilot study suggests that: (i) at least in the Luminal subclass, the ZNF217-WT-ΔE4 gene expression signature has a better prognostic value than that of ZNF217-WT mRNA levels alone, ZNF217-ΔE4 mRNA levels alone or that of the other conventional prognostic parameters of breast cancer; and (ii) in the Luminal subclass, the prognostic value of the ZNF217-WT-ΔE4 gene expression signature recapitulates that of ZNF217-E3 mRNA levels.

## Constitutive Expression of ZNF217-ΔE4 In Breast Cancer Cells Promotes An Aggressive Phenotype

We have previously described the deleterious effect of ZNF217-WT expression after stably transfecting MDA-MB-231 cells, a

<footnote><sup>a</sup>HR, Hazard ratio.</footnote>

<footnote><sup>b</sup>95% CI, 95% confidence interval.</footnote>

<footnote><sup>c</sup>p was considered significant when p < 0.05. NS, not significant.</footnote>

<footnote><sup>d</sup>subclasses of breast cancer were determined (ER, PR, HER2) according to the St Gallen recommendation (12-14).</footnote>

<footnote><sup>e</sup>N/A, not applicable as all the cases are censored in the low ZNF217-ΔE4 (E3-E5) mRNA level and in the low ZNF217-WT-ΔE4 signature groups.</footnote>
relevant breast cancer cell model because possessing low endogenous levels of ZNF217 (4, 5). With the aim to study the impact of ZNF217-ΔE4 expression on breast cancer cells phenotype, we established stable MDA-MB-231 cells constitutively overexpressing the ZNF217-ΔE4 protein. Two cell clones (ZNF217-ΔE4-1G3 and ZNF217-ΔE4-2C2), overexpressing ZNF217-ΔE4 mRNA (Figure 2A) and ZNF217-ΔE4 protein (Figure 2B), were selected under blasticidin, and their phenotype was compared with that of MDA-MB-231-pcDNA6 control cells and of previously established ZNF217-WT cells (4, 5). ZNF217-WT cells did not display any significant increase in ZNF217-ΔE4 expression levels (Supplementary Figure 3). Our previous work demonstrated that ectopic expression of ZNF217-WT in MDA-MB-231 cells triggers an aggressive phenotype [increased cell proliferation, increased anchorage-independent growth, resistance to paclitaxel (4, 5)], as illustrated in Figures 2C-F. We newly found that the constitutive expression of ZNF217-ΔE4 led to a significant increase in proliferation of both ZNF217-ΔE4-1G3 and ZNF217-ΔE4-2C2 cells, compared to MDA-MB-231-pcDNA6 control cells or to ZNF217-WT cells (Figure 2C). To investigate whether ZNF217-ΔE4 alters the response to the chemotherapeutic drug paclitaxel, we performed dose-response experiments to measure IC50 values. Strikingly, constitutive expression of ZNF217-ΔE4 led to significant increase in cell viability in the presence of paclitaxel (Figure 2D), with a relative resistance of 3.2- and 3.8-fold, respectively, for ZNF217-ΔE4-1G3 and ZNF217-ΔE4-2C2 cells (IC50<sub>MDA-MB-231-pcDNA6</sub> = 1.8 ± 0.3×10<sup>-9</sup> M, IC50<sub>ZNF217-ΔE4-1G3</sub> = 5.7 ± 0.4×10<sup>-9</sup> M, IC50<sub>ZNF217-ΔE4-2C2</sub> = 7.0 ± 1.0×10<sup>-9</sup> M). The relative resistance developed by the ZNF217-WT cells was of 2.3 (with an IC50<sub>ZNF217-WT</sub> = 4.2 ± 0.85×10<sup>-9</sup> M). Anchorage-independent growth is a phenotype commonly associated with aggressiveness and metastasis. We found that the number and the size of colonies formed in soft agar when ZNF217-ΔE4 was constitutively overexpressed were significantly higher in ZNF217-ΔE4-1G3 and ZNF217-ΔE4-2C2 cells than in control cells or than in ZNF217-WT cells (*p* < 0.001, Figures 2E, F). Altogether, these data demonstrate that overexpression of ZNF217-ΔE4 in breast cancer cells enhances their aggressiveness and that the impact of ZNF217-ΔE4 expression on the phenotype of breast cancer cells seems to be at least as or more efficient than that of ZNF217-WT.

## Ectopic Over-Expression of ZNF217-ΔE4 Is Associated With Increased Expression of ZNF217-WT

Having demonstrated that ectopic and constitutive overexpression of ZNF217-ΔE4 in MDA-MB-231 cells confers a phenotype similar to that previously described when ZNF217-WT is constitutively overexpressed (4, 5), we sought to clarify the expression levels of endogenous ZNF217-WT in ZNF217-ΔE4-1G3 and ZNF217-ΔE4-2C2 cells. Strikingly, both cell lines displayed increased mRNA levels and protein levels of endogenous ZNF217-WT (Figures 3A, B). The molecular mechanisms regulating the transcription of the ZNF217 gene

are barely known. We cloned the -2110 base pairs to +300 base pairs region of the ZNF217 gene into a firefly luciferase reporter plasmid, but we did not observe any increased firefly signal when pcDNA6-ZNF217-ΔE4 plasmid was co-transfected (data not shown), ruling out a positive direct transcriptional regulation by the ZNF217-ΔE4 isoform. Converging epigenetic data stress out that the methylation status at CpG sites within the ZNF217 locus correlates with inverse ZNF217 expression levels (7, 15-18). We thus investigated the DNA methylation status at the ZNF217 CpG sites in both ZNF217-ΔE4-1G3 and ZNF217-ΔE4-2C2 cell lines, in comparison with control MDA-MB-231-pcDNA6 cells. Our methylome analyses revealed in both ZNF217-ΔE4-1G3 and ZNF217-ΔE4-2C2 cells (but not in ZNF217-WT cells) the demethylated DNA status of three sites belonging to a CpG island included in exon 1 (GRCh37/hg19) of the ZNF217 gene: chr20:52198378 (cg01692482), chr20:52198279 (cg00703481), and chr20:52198225 (cg12032027). Heat map visualizing the DNA methylation status at the three identified CpG sites (cg01692482, cg00703481, and cg12032027) illustrates the consensus clustering result, demonstrating DNA demethylated status in both ZNF217-ΔE4-1G3 and ZNF217-ΔE4-2C2 cell lines compared to MDA-MB-231-pcDNA6 control cells or ZNF217-WT cells (Figure 3C). Pyrosequencing analyses validated this finding for cg01692482 and cg00703481 (Figure 3D). Specific primers targeting cg12032027 could not be designed, thus not allowing to explore this CpG site by pyrosequencing analysis. Altogether, our data highlighted that ectopic overexpression of ZNF217-ΔE4 leads to increased expression on endogenous ZNF217-WT correlated with an inverse DNA methylation status of the ZNF217 gene.

## DISCUSSION

A growing body of evidence suggests that ASE significantly affect tumor progression, including in the breast. Several studies have indicated that breast cancer-specific ASE include the generation of proteins with new functions (8, 9, 19, 20). Such studies have also underscored the richness of ASE as a source of biomarkers for breast cancer and attempted to define transcript isoform signatures that are associated with molecular subtypes, tumor grade, ER status, or metastasis (8, 9, 19-22).

The first key finding of our study was the identification of ZNF217-ΔE4 mRNA in primary breast tumors, indicating that the exon 4-skipping process of the ZNF217 pre-mRNA occurs in breast cancer. This was thoroughly validated by RNAseq or RT-qPCR analyses in primary breast tumors from independent cohorts. Alternative splicing of other members of the Krüppel-like factor (KLF) gene family, including the KLF6 tumor suppressor and the KLF4 regulator of stem cell pluripotency, has been previously described and shown to result in the production of several isoforms with specific biological functions and biomarker value (22-25). Our study thus highlights ZNF217 as a further member of the KLF family that is susceptible to splicing modulation in a tumorigenic context, generating distinct isoforms. We next sought to explore the
A

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

B

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

C

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

D

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

E

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

F

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

**FIGURE 2** | Constitutive expression of ZNF217-ΔE4 enhances cell proliferation, paclitaxel resistance and anchorage-independent growth of MDA-MB-231 cells. (A) RT-qPCR analysis of ZNF217-ΔE4 (E3-E5) mRNA expression levels in MDA-MB-231-ΔE4 cell lines and MDA-MB-231-pcDNA6 control cells (mean ± SD from three independent experiments, a.u., arbitrary units). (B) Representative Western blot analysis of ZNF217-ΔE4 expression levels in MDA-MB-231-ΔE4 cell lines versus control cells using RM217 antibody. (C) Cell proliferation assessment of the two MDA-MB-231-ΔE4 cell lines (means ± SD of three independent experiments). (D) Cell viability after 96 hours of paclitaxel exposure as assessed by MTS cytotoxicity assay. (E) Average number of colonies per well, and (F) average size of the colonies in soft agar. (a.u.), arbitrary units. Data illustrated in (E, F) represent mean ± SD of three independent experiments. In (A, C-F): \*p < 0.05, \*\*p < 0.01 and \*\*\*p < 0.001 versus MDA-MB-231-pcDNA6 cells (Student's t-test). In (C-F): \#p < 0.05, \##p < 0.01 and \*\*\#p < 0.001 versus ZNF217-WT cells (Student's t-test).

clinical and functional relevance of ZNF217-ΔE4 expression in breast cancer.

Using a primer pair targeting exon 3, we have previously demonstrated that the ZNF217 mRNA expression level of a

breast tumor is informative and provides a powerful biomarker of poor prognosis, in particular in the Luminal subclass (5, 7). The present pilot retrospective study with 10-years follow-up provided evidence that the novel ZNF217-WT-ΔE4 gene
<figure><img src="image_14.png" /></figure>

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

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

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

**FIGURE 3** | Constitutive ZNF217-ΔE4 expression leads to increased expression of ZNF217-WT and demethylated DNA status of the ZNF217 gene. (A) RT-qPCR analysis of ZNF217-WT (E3-E4) mRNA expression levels in MDA-MB-231-ΔE4 cell lines and MDA-MB-231-pcDNA6 control cells (mean ± SD from three independent experiments, a.u., arbitrary units) \*\*\*p < 0.001, in Student t-test. (B) Representative Western blot analysis of ZNF217 WT protein expression levels using Abcam #48133 (upper panel) and Covalab 2 antibodies (lower panel). (C) β-values of 3 CpG sites located in the first Exon of the ZNF217 gene in MDA-MB-231-pcDNA6 control cells, in MDA-MB-231-ZNF217-ΔE4 and in MDA-MB-ZNF217-231-WT cell lines were calculated after genome-wide methylation analysis with Illumina Infinium MethylationEPIC Beadchip. R1 and R2, independent genomic replicates. Heatmaps were realized with pheatmap package (Version: 1.0.12) in RStudio (Version 1.4.1106). Demethylated positions validated by pyrosequencing are indicated by an asterisk. (D) Heatmap of pyrosequencing analysis of cg01692482 and cg00703481 in MDA-MB-231-pcDNA6 control cells and in MDA-MB-231-ZNF217-ΔE4 and in MDA-MB-231-ZNF217-WT cell lines. R1 and R2, independent genomic replicates.

expression signature, which combines the expression of the two isoforms, retained clinical significance, and was significantly associated with shorter RFS in the Luminal subclass. The ZNF217-WT-ΔE4 signature is thus capable of resuming the prognostic value of ZNF217-E3 mRNA levels by allowing the stratification of the “good prognosis” Luminal subclass into good

or poor/intermediate outcome. When considering the “good prognosis” SBRI+II subgroup, again the ZNF217-WT-ΔE4 signature could be used to refine the categorization of patients into good or poor/intermediate outcome. Our gene expression signature, however, was not informative in the “bad prognosis” subclasses such as HER2-enriched, Triple-negative, or SBRIII
subgroups. Importantly, our univariate and multivariate analyses conducted in the Luminal subclass demonstrated that the prognostic value of the ZNF217-WT-ΔE4 dual signature was more powerful than the expression levels of each isoform only, or than any of the usual clinical prognostic parameters. Altogether, our pilot retrospective study indicates the need to assess the mRNA expression levels of the two ZNF217 isoforms to obtain the most powerful prognostic biomarker value. Regarding the TCGA RNAseq data, the high proportion of censored cases (90%) prevented us from performing univariate analysis. While future work on an independent cohort is needed to confirm our findings, the results of our clinical investigation strongly suggest that the protein encoded by ZNF217-ΔE4 has a functional role.

Our study additionally pinpoints the existence of intricate molecular mechanisms involving ZNF217-WT and ZNF217-ΔE4. First, both the RNA-seq and RT-qPCR investigations, conducted on independent primary breast tumor cohorts, identified a positive correlation between ZNF217-ΔE4 and ZNF217-WT mRNA levels. Second, our in vitro data revealed that ectopic overexpression of the ZNF217-ΔE4 isoform in breast cancer cells leads to: (i) increased endogenous mRNA and protein levels of the ZNF217-WT isoform; and (ii) features of cell aggressiveness (increased cell proliferation, increased anchorage-independent growth and resistance to paclitaxel) previously described when ZNF217-WT is ectopically overexpressed (4, 5). Our study thus indicates that the ZNF217-ΔE4 isoform triggers, directly or indirectly, a poor phenotype in breast cancer cells that may involve dysregulated ZNF217-WT expression levels. Little is known regarding the molecular mechanisms governing ZNF217 transcription. Converging epigenetic data, however, show that the methylation status at CpG sites within the ZNF217 locus correlates with inverse ZNF217 expression levels (7, 15-18). Importantly, our methylome investigation revealed that constitutive expression of the ZNF217-ΔE4 isoform (but not of ZNF217-WT) is paired with decreased DNA methylation status on three CpG sites located on a key exon on the ZNF217 gene. In a case-control study involving 1,083 blood samples (healthy women versus breast cancer patients), lack of methylation of the same ZNF217 exon predicted increased breast cancer risk (26), allowing the authors to propose the DNA methylation status at this ZNF217 locus to be a surrogate biomarker of breast cancer risk. We thus propose a model where the ZNF217-ΔE4 isoform plays a biologically relevant role in breast cancer, at least by impacting epigenetic-driven mechanisms governing the expression of ZNF217-WT. In future work, it would thus be of great interest to investigate the expression levels of both ZNF217-ΔE4 and ZNF217-WT as well as the DNA methylation status of the ZNF217 gene in breast cancer tissues from patients.

The predicted amino acid (aa) sequence encoded by the ZNF217-ΔE4 transcript shows that it shares a common region (aa 1 to 1012) with the ZNF217-WT protein, but the C-terminal region (49 aa) encoded by exon 5 is different from that present in ZNF217-WT protein and encoded by most of exon 4 (36 aa). The C-terminal region of ZNF217 is of the greatest interest because: (i) it possesses the transcriptional repressor domain (27) and a

proline-rich domain (2) that might be involved in the ZNF217's transcriptional activator properties; and (ii) it binds to CtBP1 and CtBP2 transcriptional co-repressors (27, 28) and to the ER, modulating ER-signaling (7). Future work is thus required to determine how the ZNF217-ΔE4 interferes with the epigenetic machinery and whether it retains any transcriptional activity.

Little is known about the molecular mechanisms governing ASE in specific biological contexts. However, while the detailed mechanisms are still elusive, many recent studies have revealed that RNA N6-methyladenosine (m⁶A) modification serves as a mark for the recruitment of splicing factors and then affects ASE (20, 29). Interestingly, one elegant study has deciphered that the ZFP217 transcription factor (the murine ortholog of the human ZNF217) catalyzes m⁶A methylation at key transcripts, of which was the ZFP217 RNA itself (at exon3 and exon 5) (30). One thus cannot exclude that such ZNF217-WT-driven mechanisms may control splice events leading to the production of the ZNF217-ΔE4 isoform.

Overall, our new findings have important medical applications. Indeed, we have provided insights into the biological and clinical significance of the novel exon4-skipping ZNF217 isoform in breast cancer. Our pilot retrospective analysis emphasized the need to assess both ZNF217-WT and ZNF217-ΔE4 expression levels to obtain the most powerful biomarker value, and suggested a key biological function of the encoded ZNF217-ΔE4 protein isoform. Our *in vitro* data supported our clinical data and reported for the first time that the ZNF217-ΔE4 protein drives cell aggressiveness and increased ZNF217-WT expression levels. In conclusion, our data highlight that assessing the expression levels of ZNF217-WT and ZNF217-ΔE4 isoforms may serve as a novel prognostic biomarker allowing a better stratification of breast cancers with good prognosis and helping clinicians in therapeutic decisions. Future studies on the biological activities of ZNF217-ΔE4 protein and its interplay with ZNF217-WT will lead to a better understanding of ZNF217 oncogenic machinery.

## DATA AVAILABILITY STATEMENT

Publicly available datasets were analyzed in this study. This data can be found here: https://www.nature.com/articles/nature11412.

## ETHICS STATEMENT

This study has been approved by the local ethics committee (CRB Centre Léon Bérard, France). The CRB Centre Léon Bérard is quality certified according to NFS96-900 French standard and ISO 9001 for clinical trials, ensuring scientific rigor for sample conservation, traceability, and quality, as well as ethical rules observance and defined rules for transferring samples for research purposes (Ministry of Health for activities authorization no AC-2019-3426 and DC-2008-99). The material used in the study has been collected in agreement with all applicable laws, rules, and requests of French and
European government authorities. All subjects gave written informed consent. The patients/participants provided their written informed consent to participate in this study.

## AUTHOR CONTRIBUTIONS

BG and LP performed the biostatistical analyses of the RNAseq data. AB validated the ZNF217 isoforms specific primers and performed in vitro experiments. DL and PAC performed the biostatistical analyses of the 107-breast tumor cohort. AW and JL performed and analyzed the methylome study. JV and JS performed the pyrosequencing experiments and analysis. JV and PC participated in scientific discussions. DL, BG, and PAC wrote the manuscript. PAC designed and supervised the study. All authors contributed to the article and approved the submitted version.

## FUNDING

This research program was supported by grants from the Pack Ambition International 2020 (Région Rhône-Alpes, France) and the French Ligue Contre le Cancer (committee 71). DL was funded by France Excellence Scholarship Program of the French Embassy in Vietnam. AB was supported by a Ph.D. grant from the French Ligue Nationale contre le Cancer. BG and LP were financed by the 2018-2.1.17-TET-KR-00001 and 2018-1.3.1-VKE-2018-00032 grants and the Higher Education Institutional Excellence Programme (2020-4.1.1.-TKP2020) of the Ministry for Innovation and Technology in Hungary, within the framework of the Bionic thematic programme of the Semmelweis University.

## ACKNOWLEDGMENTS

We thank Dr. H. Mattock for her kind help in the edition. The authors thank the Centre de Ressources Biologiques (CRB) of the Centre Léon Bérard (BB-0033-00050). The authors wish to acknowledge the support of ELIXIR Hungary (www.elixir-

## REFERENCES

1. Ferlay J, Ervik M, Lam F, Mery L, Piñeros M, Znaor A, et al. Global Cancer Observatory: Cancer Today (2020). Lyon, France: International Agency for Research on Cancer. Available at: https://gco.iarc.fr/today (Accessed 23/12/2020).

2. Collins C, Rommens JM, Kowbel D, Godfrey T, Tanner M, Hwang SI, et al. Positional Cloning of ZNF217 and NABC1: Genes Amplified at 20q13.2 and Overexpressed in Breast Carcinoma. Proc Natl Acad Sci USA (1998) 95 (15):8703-8. doi: 10.1073/pnas.95.15.8703

3. Cohen PA, Donini CF, Nguyen NT, Lincet H, Vendrell JA. The Dark Side of ZNF217, A Key Regulator of Tumorigenesis With Powerful Biomarker Value. Oncotarget (2015) 6(39):41566-81. doi: 10.18632/oncotarget.5893

4. Thollet A, Vendrell JA, Payen L, Ghayad SE, Ben Larbi S, Grisard E, et al. ZNF217 Confers Resistance to the Pro-Apoptotic Signals of Paclitaxel and

hungary.org). We are grateful to Dr. Nguyen for generating the pcDNA6-ZNF217-ΔE4 plasmid and Miss S. Crose for technical help in methylome studies.

## SUPPLEMENTARY MATERIAL

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2021.647269/full#supplementary-material

**Supplementary Figure 1** | The sequence-specific primers for ZNF217-WT (E3-E4) and ZNF217-ΔE4 (E3-E5) isoforms with high specificity and sensitivity. (A) Binding sites of the primer pair 1, 2, and 3 used to specifically detect the ZNF217-E3, ZNF217-WT (E3-E4), and ZNF217-ΔE4 (E3-E5) isoforms, respectively. The primer pair 1 hybridized within the sequence of ZNF217's exon 3. The forward primer of pair 2 hybridized onto the exon 3 - exon 4 junction of ZNF217-WT (E3-E4) isoform. The reverse primer of pair 3 hybridized onto the exon 3 - exon 5 junction of ZNF217-ΔE4 (E3-E5) isoform. (B) RT-qPCR detection of ZNF217 isoforms in mixtures containing 1,000 copies of different ratios of pcDNA6-ZNF217-WT and/or pcDNA6-ZNF217-ΔE4 plasmids amplified by RT-qPCR using the primer pair 1 (black bar), primer pair 2 (grey bar), or primer pair 3 (white bar). All data represent as mean ± standard deviation of two independent experiments conducted in triplicate. Frequency of (C) ZNF217-ΔE4 (E3-E5) mRNA levels and of (D) ZNF217-ΔE4 (E3-E5) mRNA levels among the 107 primary breast tumors.

**Supplementary Figure 2** | The ZNF217-WT-ΔE4 gene expression signature is of poor prognosis and is associated with shorter relapse-free survival in the Luminal breast cancer subclass. Kaplan-Meier analyses (univariate analyses) for RFS in Luminal subclass (ER+ and/or PR +) of (A) ZNF217-E3 mRNA expression levels; (B) ZNF217-WT (E3-E4) mRNA expression levels; (C) ZNF217-ΔE4 (E3-E5) mRNA expression levels; (D) the 3-groups ZNF217-WT-ΔE4 gene expression signature defined by: (group 1), low ZNF217-WT (E3-E4) and low ZNF217-ΔE4 (E3-E5) mRNA levels; (group 2), high ZNF217-WT (E3-E4) or high ZNF217-ΔE4 (E3-E5) mRNA levels; (group 3), high ZNF217-WT (E3-E4) and high ZNF217-ΔE4 (E3-E5) mRNA levels; (E) ZNF217-WT-ΔE4 gene expression signature. *p* < 0.05.

**Supplementary Figure 3** | Characterization of MDA-MB-231-ZNF217-WT cells. (A) RT-qPCR analysis of ZNF217-WT (E3-E4) mRNA expression levels in MDA-MB-231-ZNF217-WT cell line and MDA-MB-231-pcDNA6 control cells (mean ± SD from three independent experiments, a.u., arbitrary units) ** *p* < 0.001. in Student t-test. (B) Representative Western blot analysis of ZNF217 WT protein expression levels using Abcam #48133 (upper panel) and Covalab 2 antibodies (lower panel). (C) RT-qPCR analysis of ZNF217-ΔE4 (E3-E5) mRNA expression levels in MDA-MB-231-ZNF217-WT cell line and in MDA-MB-231-pcDNA6 control cells (mean ± SD from three independent experiments, a.u., arbitrary units). (D) Representative Western blot analysis of ZNF217-ΔE4 expression levels in MDA-MB-231-ZNF217-WT cell line versus control cells using RM217 antibody.

Aberrant Expression of Aurora-a in Breast Cancer Cells. Mol Cancer (2010) 9:291. doi: 10.1186/1476-4598-9-291

5. Vendrell JA, Thollet A, Nguyen NT, Ghayad SE, Vinot S, Bièche I, et al. ZNF217 Is a Marker of Poor Prognosis in Breast Cancer That Drives Epithelial-Mesenchymal Transition and Invasion. Cancer Res (2012) 72 (14):3593-606. doi: 10.1158/0008-5472.CAN-11-3095

6. Bellanger A, Donini CF, Vendrell JA, Lavaud J, Machuca-Gayet I, Ruel M, et al. The Critical Role of the ZNF217 Oncogene in Promoting Breast Cancer Metastasis to the Bone. J Pathol (2017) 242(1):73-89. doi: 10.1002/path.4882

7. Nguyen NT, Vendrell JA, Poulard C, Gyorffy B, Goddard-Leon S, Bieche I, et al. A Functional Interplay Between ZNF217 and Estrogen Receptor Alpha Exists in Luminal Breast Cancers. Mol Oncol (2014) 8(8):1441-57. doi: 10.1016/j.molonc.2014.05.013

8. Ladomery M. Aberrant Alternative Splicing is Another Hallmark of Cancer. Int J Cell Biol (2013) 2013:463786. doi: 10.1155/2013/463786
9. Martinez-Montiel N, Anaya-Ruiz M, Perez-Santos M, Martinez-Contreras RD. Alternative Splicing in Breast Cancer and the Potential Development of Therapeutic Tools. *Genes (Basel)* (2017) 8(10):217. doi: 10.3390/genes8100217

10. Koboldt DC, Fulton RS, McLellan MD, Schmidt H, Kalicki-Veizer J, McMichael JF, et al. Comprehensive Molecular Portraits of Human Breast Tumours. *Nature* (2012) 490(7418):61–70. doi: 10.1038/nature11412

11. Stephenson JR, Axelrad AA, McLeod DL, Shreeve MM. Induction of Colonies of Hemoglobin-Synthesizing Cells by Erythropoietin *In Vitro. Proc Natl Acad Sci USA* (1971) 68(7):1542–6. doi: 10.1073/pnas.68.7.1542

12. Gnant M, Harbeck N, Thomssen C. St. Gallen 2011: Summary of the Consensus Discussion. *Breast Care (Basel)* (2011) 6(2):136–41. doi: 10.1159/000328054

13. Goldhirsch A, Wood WC, Coates AS, Gelber RD, Thürlimann B, Senn HJ, et al. Strategies for Subtypes–Dealing With the Diversity of Breast Cancer: Highlights of the St. Gallen International Expert Consensus on the Primary Therapy of Early Breast Cancer 2011. *Ann Oncol* (2011) 22(8):1736–47. doi: 10.1093/annonc/mdr304

14. Goldhirsch A, Winer EP, Coates AS, Gelber RD, Piccart-Gebhart M, Thürlimann B, et al. Personalizing the Treatment of Women With Early Breast Cancer: Highlights of the St Gallen International Expert Consensus on the Primary Therapy of Early Breast Cancer 2013. *Ann Oncol* (2013) 24 (9):2206–23. doi: 10.1093/annonc/mdt303

15. Etcheverry A, Aubry M, de Tayrac M, Vauleon E, Boniface R, Guenot F, et al. DNA Methylation in Glioblastoma: Impact on Gene Expression and Clinical Outcome. *BMC Genomics* (2010) 11:701. doi: 10.1186/1471-2164-11-701

16. Leu YW, Yan PS, Fan M, Jin VX, Liu JC, Curran EM, et al. Loss of Estrogen Receptor Signaling Triggers Epigenetic Silencing of Downstream Targets in Breast Cancer. *Cancer Res* (2004) 64(22):8184–92. doi: 10.1158/0008-5472.CAN-04-2045

17. Ness C, Katta K, Garred O, Kumar T, Olstad OK, Petrovski G, et al. Integrated Differential DNA Methylation and Gene Expression of Formalin-Fixed Paraffin-Embedded Uveal Melanoma Specimens Identifies Genes Associated With Early Metastasis and Poor Prognosis. *Exp Eye Res* (2021) 203:108426. doi: 10.1016/j.exer.2020.108426

18. Renner M, Wolf T, Meyer H, Hartmann W, Penzel R, Ulrich A, et al. Integrative DNA Methylation and Gene Expression Analysis in High-Grade Soft Tissue Sarcomas. *Genome Biol* (2013) 14(12):r137. doi: 10.1186/gb-2013-14-12-r137

19. Trincado JL, Sebestyen E, Pages A, Eyras E. The Prognostic Potential of Alternative Transcript Isoforms Across Human Tumors. *Genome Med* (2016) 8(1):85. doi: 10.1186/s13073-016-0339-3

20. Yang Q, Zhao J, Zhang W, Chen D, Wang Y. Aberrant Alternative Splicing in Breast Cancer. *J Mol Cell Biol* (2019) 11(10):920–9. doi: 10.1093/jmcb/mjz033

21. Venables JP, Klinck R, Bramard A, Inkel L, Dufresne-Martin G, Koh C, et al. Identification of Alternative Splicing Markers for Breast Cancer. *Cancer Res* (2008) 68(22):9525–31. doi: 10.1158/0008-5472.CAN-08-1769

22. Yang Y, Xiong J, Wang J, Ruan Y, Zhang J, Tian Y, et al. Novel Alternative Splicing Variants of Klf4 Display Different Capacities for Self-Renewal and Pluripotency in Mouse Embryonic Stem Cells. *Biochem Biophys Res Commun* (2020) 532(3):377–84. doi: 10.1016/j.bbrc.2020.08.054

23. Syafruddin SE, Mohtar MA, Wan Mohamad Nazarie WF, Low TY. Two Sides of the Same Coin: The Roles of KLF6 in Physiology and Pathophysiology. *Biomolecules* (2020) 10(10):1378. doi: 10.3390/biom10101378

24. Wang L, Shen F, Stroehlein JR, Wei D. Context-Dependent Functions of KLF4 in Cancers: Could Alternative Splicing Isoforms Be the Key? *Cancer Lett* (2018) 438:10–6. doi: 10.1016/j.canlet.2018.09.005

25. Zhang N, Yan QQ, Lu L, Shao JB, Sun ZG. The KLF6 Splice Variant KLF6-SV1 Promotes Proliferation and Invasion of Non-Small Cell Lung Cancer by Up-Regulating PI3K-AKT Signaling Pathway. *J Cancer* (2019) 10(22):5324–31. doi: 10.7150/jca.34212

26. Widschwendter M, Apostolidou S, Raum E, Rothenbacher D, Fiegl H, Menon U, et al. Epigenotyping in Peripheral Blood Cell DNA and Breast Cancer Risk: A Proof of Principle Study. *PloS One* (2008) 3(7):e2656. doi: 10.1371/journal.pone.0002656

27. Cowger JJM, Zhao Q, Isovici M, Torchia J. Biochemical Characterization of the Zinc-Finger Protein 217 Transcriptional Repressor Complex: Identification of a ZNF217 Consensus Recognition Sequence. *Oncogene* (2007) 26(23):3378–86. doi: 10.1038/sj.onc.1210126

28. Quinlan KGR, Nardini M, Verger A, Francescato P, Yaswen P, Corda D, et al. Specific Recognition of ZNF217 and Other Zinc Finger Proteins At a Surface Groove of C-Terminal Binding Proteins. *Mol Cell Biol* (2006) 26(21):8159–72. doi: 10.1128/MCB.00680-06

29. Adhikari S, Xiao W, Zhao YL, Yang YG. M⁵A: Signaling for mRNA Splicing. *RNA Biol* (2016) 13(9):756–9. doi: 10.1080/15476286.2016.1201628

30. Aguilo F, Zhang F, Sancho A, Fidalgo M, Di Cecilia S, Vashisht A, et al. Coordination of m⁶A mRNA Methylation and Gene Transcription by ZFP217 Regulates Pluripotency and Reprogramming. *Cell Stem Cell* (2015) 17(6):689–704. doi: 10.1016/j.stem.2015.09.005

**Conflict of Interest:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Copyright © 2021 Bellanger, Le, Vendrell, Wierinckx, Pongor, Solassol, Lachuer, Clezardin, Győrfy and Cohen. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.