The clinical significance of epigenetic and RNAPII variabilities occurring in clear cell renal cell carcinoma as a potential prognostic marker

(1)

Translational Oncology 20 (2022) 101420

Available online 10 April 2022

Original Research

The clinical significance of epigenetic and RNAPII variabilities occurring in clear cell renal cell carcinoma as a potential prognostic marker

N ´ ora Ord ¨ ¨ og, Barbara N Borsos, Hajnalka Majoros, Zsuzsanna Ujfaludi, Gabriella Pankotai-Bod ´ o, Sarolta Bank ´ o, Farkas Sük ¨ osd, Levente Kuthi

^**

, Tibor Pankotai

^*

Institute of Pathology, Albert Szent-Gy¨orgyi Medical School, University of Szeged, 1 Allom´ ´as Street, Szeged H-6725, Hungary

A R T I C L E I N F O Keywords:

ccRCC epigenetics

H3K4me³and H3K9me³ RNAPII

Renal cell carcinoma γH2A.X

A B S T R A C T

Patients diagnosed with clear cell renal cell carcinoma (ccRCC) have poor prognosis for recurrence and approximately 30–40% of them will later develop metastases. For this reason, the appropriate diagnosis and the more detailed molecular characterisation of the primary tumour, including its susceptibility to metastasis, are crucial to select the proper adjuvant therapy by which the most prosperous outcome can be achieved. Nowadays, clinicopathological variables are used for classification of the tumours. Apart from these, molecular biomarkers are also necessary to improve risk classification, which would be the most beneficial amongst modern adjuvant therapies. As a potential molecular biomarker, to follow the transcriptional kinetics in ccRCC patients (n=30), we analysed epigenetic changes (γH2A.X, H3K4me³, and H3K9me³) and the alterations in the level of RNA polymerase II (RNAPII) by immunohistochemical staining on dissected tissue sections. The variabilities between the tumorous and non-tumorous parts of the tissue were detected using quantitative image analysis by monitoring 30 cells from different positions of either the tumorous or the non-tumorous part of the tissue sections. Data obtained from the analyses were used to identify potential prognostic features and to associate them with the progression. These markers might have a value to predict patient outcomes based on their individual cellular background. These results also support that detection of any alteration in the level of H3K4me³, H3K9me³, and γH2A.X can account for valuable information for presuming the progression of ccRCC and the clinical benefits to select the most efficient personalised therapy.

Introduction

Among the various types of renal cell carcinoma (RCC), clear cell renal cell carcinoma (ccRCC) is the most frequently occurred subtype, accounting for approximately 70–80% of all cases, and it has the sixth highest mortality rate in Europe and US [1]. According to the WHO database, RCC is considered one of the most common epithelial cancers worldwide, with 430,000 newly diagnosed and 180,000 death cases in 2020. The mortality rate and severeness of the disease are highly related to late-stage diagnoses as only 10% of the patients developed typical symptoms, like haematuria, backache, or abdominal mass [2]. The current WHO classification distinguishes 14 RCC subsets, including clear cell RCC, papillary RCC, and chromophobe RCC, as the most frequent subtypes based on morphological and histopathological characteristics [3–5]. Due to the lack of molecular biomarkers, RCC classification is

difficult, and the adequate identification of stages requires reliable tools.

Additionally, it is also important to develop proper tools for predicting clinical outcomes to improve the standard treatment guidelines for RCC patients. In modern diagnostics, novel methods include prognostic modelling which relies on DNA-, RNA-, and protein-based biomarkers, to predict the possible behaviour of the tumour [6,7]. For patients diagnosed with high-risk tumour, adjuvant therapeutic inventions are crucial, while for patients diagnosed with low-risk cancer, overtreatment can be dangerous [8–10]. Applying large-scale genomic and proteomic analyses, the establishment of biomarkers with diagnostic quality is feasible, and alterations associated with RCC can be identified [11]. These techniques also have the potential to reveal novel therapeutic targets which can be implicated into unprecedented diagnostic and therapeutic strategies. However, to describe the overall perfor- mance of the tumour progression, an individual biomarker does not

* Corresponding author at: Institute of Pathology, Szegedi Tudomanyegyetem Szent-Gyorgyi Albert Orvostudomanyi Kar, Allom´ ´as street 1, 6725 Szeged, Hungary.

** co-correspondance.

E-mail addresses: kuthi.levente@med.u-szeged.hu (L. Kuthi), pankotai.tibor@szte.hu (T. Pankotai).

Contents lists available at ScienceDirect

Translational Oncology

journal homepage: www.elsevier.com/locate/tranon

https://doi.org/10.1016/j.tranon.2022.101420

Received 9 November 2021; Received in revised form 2 March 2022; Accepted 5 April 2022

(2)

Translational Oncology 20 (2022) 101420

2 provide valuable data in terms of sensitivity and specificity [12–14].

In the progression of ccRCC, mutations of the Von Hippel-Lindau Tumor Supressor (VHL) tumour suppressor gene are one of the main drivers [15,16]. Additionally, mutations of the Polybromo 1 (PBRM1), BRCA1 Associated Protein 1 (BAP1), and SET domain containing 2 (SETD2) are also found to be associated with ccRCC but with a lower frequency [17,18]. This genetic heterogeneity observed in ccRCC leads to diffi- culties in achieving a prosperous treatment and promotes recurrence.

These mutations in PBRM1, SETD2, and BAP1 genes are strongly associated with epigenetic changes since the proteins produced from these genes are responsible for modulating the chromatin structural changes [18,19]. The PBRM1 gene encodes the bromodomain-containing protein BAF180, which is a subunit of the Polybromo- and BRG1-associated factors containing complex (PBAF) SWItching defective/Sucrose Non-Fermenting (SWI/SNF) chromatin remodelling complex. BAP1 encodes the histone deubiquitinating enzyme BRCA1-Associated Protein 1, while SETD2 encodes a methyltransferase that is specifically trimethy- lates at lysine 36 of histone H3 (H3K36me³) [20]. The loss of these genes leads to abnormal epigenetic alterations, providing new opportunities for the diagnostic application of these proteins and for developing novel therapeutic strategies in ccRCC. Some drugs, targeting the epigenetic system, are currently under investigation; however, strategies that combine therapies targeting the epigenetic machinery with conven- tional therapies are still in progress [21–23].

Here, we provide an overview of the epigenetic landscape of 30 patients diagnosed with ccRCC by using immunohistochemical profiling and automated quantification of the γH2A.X, H3K4me³, H3K9me³, and RNAPII level in tumorous and non-tumorous parts of the tissue. In this report, we highlight that the pattern of these marks is dependent on the tumour stages and grades. Therefore, our results establish the prognostic value of γH2A.X, H3K4me³, H3K9me³, and RNAPII patterns in ccRCC.

Materials and methods Selection of the cohort

A total of 30 patients, who underwent radical or partial nephrectomy at the University of Szeged, were selected for further investigation. This study was conducted with the permission of the Regional and Institu- tional Human Medical Biological Research Ethics Committee, University of Szeged (No. 188/2019-SZTE). All patients gave his/her written informed consent for the publication of all data and images. A histopa- thologist evaluated all the cases for the histological subtype, tumour stage, and ISUP grade. Haematoxylin and eosin staining was performed on each section to confirm the tumour content of the renal tissues. After surgery, tissue samples were dissected from both the tumorous and normal parts of the tissues, then they were flash frozen, and deposited in the Biobank of Pathology, University of Szeged. All the samples were stored at − 80^◦C till further processing.

Preparations of the normal and tumorous tissues

For immunofluorescence staining, the tissues were embedded in Shandon Cryomatrix gel (Thermo Fisher Scientific), and 5 µm sections were prepared on Superfrost Ultra Plus slides by cryostat (Cryostar NX50).

Immunostaining of frozen tissues

The tissues were fixed for 10 minutes using acetone, and then washed 3 times with PBS solution. Subsequently, the tissues were permeabilized with 0.3% Triton-X-100/PBS for 20 minutes at 25^◦C. Alternatively, for permeabilization 0.3% Na-DOC/0.3% Triton-X-100/PBS was used for 20 minutes at 25^◦C. Then sections were blocked with 5% BSA/0.3% Triton X-100/PBS for 1 hour. Samples were incubated with primary antibodies diluted in 1% BSA/PBST: anti-γH2A.X (Millipore 05-636) in 1:100, anti- H3K9me³(Abcam ab8898) in 1:100, anti-H3K4me³(Abcam ab8580) in 1:300, and RNAPII (1BP7G5 IGBMC) in 1:250 dilution. After washing Table 1

The components of the cohort

Code Sample ID number Gender Age Tumour size (mm) Stadium Metastasis Grade PBRM1 (CNV) BAP1 (CNV) SETD2 (CNV) VHL (CNV)

#1 T1207 female 62 75 pT3a adrenal gland High N N N N

#2 T1257 male 50 52 pT1b not detected Low L N N N

#3 T1337 female 61 30 pT1a not detected Low N N N N

#4 T1434 female 60 125 pT3a not detected High N N N L

#5 T1441 male 66 62 pT3a not detected High H N H L

#6 T1447 female 59 45 pT1b not detected Low N N N H

#7 T1453 male 78 60 pT3a not detected Low N N N N

#8 T1526 male 73 44 pT1b not detected High N N N N

#11 T1533 male 67 49 pT3a not detected High N N N N

#12 T1538 female 40 42 pT3a not detected Low H N N N

#13 T1622 female 62 112 pT3a not detected High N H N N

#15 T1652 female 71 125 pT3a adrenal gland High N H H H

#16 T1725 female 40 51 pT1b not detected Low N N N N

#18 T1742 female 55 32 pT1a not detected Low N H N N

#20 T747 male 68 90 pT3a not detected High L N N N

#21 T788 female 68 72 pT3a bone marrow High N H H H

#23 T795 male 57 50 pT3a not detected High L N N N

#24 T803 female 51 35 pT3a not detected Low L N N N

#25 T808 female 63 33 pT1a not detected High N H N H

#28 T827 male 78 56 pT1b not detected High N N N N

#29 T888 male 39 72 pT2a not detected High N H N N

N. Ord¨ og et al. ¨

(3)

steps, the following secondary antibodies were used: GAR Alexa 555 (Invitrogen, A21429) in 1:100, and GAM Alexa 488 (Molecular Probes, A11029) in 1:500 dilution. Finally, cells were mounted with DAPI- containing ProLong Gold antifade reagent (Life technologies). Samples were visualised with Olympus FluoView FV10i confocal microscopy.

The same exposition time was used for every image capturing. Images were quantified with ImageJ software.

Image analysis of the confocal pictures

For ImageJ analysis, 10-10 cells from 3 different parts of each tissue section were selected. The selection was established on the haematoxylin and eosin staining-based separation of the tumorous and non- tumorous tissue section which was defined by expert pathologists.

During the analysis, the size of the nuclei of each cell and their mean

grey/intensity values were recorded and added to the database. The CTCF values (Corrected Cell Fluorescence) were calculated accordingly:

integrated density–(nuclear area * average background intensity). The calculated values were visualized by box plots in Sigma Plot 12.5, and the significance of each sample was calculated by independent samples t-test in IBM SPSS Statistics 27.0.

Haematoxylin and eosin staining

Tissues were stained with haematoxylin dye, supplemented with aluminium potassium sulphate dodecahydrate (Sigma-Aldrich).

Following washing with distilled water, sections were stained with 1 g/l eosin-B solution (Sigma-Aldrich). Then samples were washed twice with distilled water. Next, the tissues were washed with isopropyl-alcohol and distilled water, then the sections were mounted. Finally, the slides Fig. 1. Adjustment of fixation and permeabilization steps of immunohistochemistry method on normal and tumorous tissue section parts from patients with clear cell renal cell carcinoma (ccRCC). A) Fixation with acetone or formaldehyde, B) Permeabilization with solely Triton X-100 or the combinatorial usage of Triton X-100 and Na-DOC. DAPI (blue) was used to visualize the nuclei. Scale bars represent 30 µm.

(4)

4 were scanned with 3DHISTECH Pannoramic Midi digital slide scanner for monitoring histoarchitectural changes.

CNV (copy number variation) measurement

CNV values were determined by RT-qPCR reactions performed on Thermo Scientific PIKO 96-well Thermal Cycler (Thermo Fischer Sci- entific) platform using SYBR Green chemistry (Promega, GoTaq qPCR Master Mix) as it is described in Ujfaludi et al. 2022 (manuscript accepted in POR).

Statistical analysis for correlation study

Pearson correlation was executed to verify the strength of a linear association between the examined markers. For this, first mean value was calculated from the CTCF values of each specimen, then the mean value of the tumorous part was divided with the mean value of the non- tumorous part of each tissue section in case of each marker. Subse- quently, using IBM SPSS Statistics 27.0, these values were applied for determination of the Pearson correlation coefficient (referred to as ‘r’) between the examined markers.

Results

Selection and classification of the cohort

The age distribution of the 30 candidates selected for the current study was between 39 to 78 years, with an average age of 62. 14 tissue sections were derived from females, and 16 were dissected from male patients, representing an almost 1:1 gender ratio. According to the haematoxylin-eosin staining of the dissected tissue, patients were cate- gorized based on either tumour stage (stadium pT1-3) or grade (low

[ISUP grade 1+2] or high [ISUP grade 3+4]) by an expert pathologist. In 3 patients, metastases were also detected in the bone marrow or in the adrenal gland. The patients selected for this study were also characterised in order to test their von Hippel-Lindau Tumor Supressor (VHL), Polybromo 1 (PBRM1), BRCA1-Associated Protein 1 (BAP1), and SET domain containing 2 (SETD2) status. We detected pathological mutation of the VHL gene in approximately one third of the samples. Although in most cases, the copy number variation of these genes was close to the control, in some cases they were found to be haploid: (I) PBRM1 in three, (II) SETD2 in four, and (III-IV) BAP1 and VHL in five patients. The detailed clinical and genetic parameters are represented in Table 1. All the 30 patients were further analysed for H3K4me³, H3K9me³, H2AXS139P (referred to as γH2A.X), and RNA polymerase II (RNAPII) levels.

First, the reliability and specificity of the antibodies against H3K4me³, H3K9me³, γH2A.X, and RNAPII on frozen material were optimised by setting the immunostaining procedure on slices derived from ccRCC and the surrounding non-tumorous parts of the tissue (Fig. 1). As expected, strong nuclear expression of the RNAPII protein was observed in all the conditions we tested, underlying the specificity of the antibody used for this study [24]. Intriguingly, the widely used formaldehyde fixation resulted in both nuclear and cytoplasmic RNAPII positivity in the non-tumorous parts and loss of signal in the ccRCC regions (Fig. 1A). To achieve comparable staining efficiency from the different tissue parts, we tried acetone fixation which appeared to be applicable for further investigations (Fig. 1A). Moreover, we also tested different permeabilization methods, but it did not affect the staining efficiency (Fig. 1B).

Fig. 2. Immunofluorescent spectra of H3K4me³detected in normal and ccRCC tumorous tissue parts. A) N represents non-tumorous tissue parts, while T corresponds to ccRCC sections. Representative images of H3K4me³staining in peritumoral (left) and tumour tissues (right) from ccRCC patients. For the better visualization of the changes in H3K4me³level, one cell is represented in a magnified resolution in each specimen. Scale bars represent 30 µm and 1.5 µm, respectively. B) The observed alterations in the quantification of relative fluorescence intensity (CTCF) of H3K4me³level between the tumorous and non-tumorous tissue parts obtained from automated image acquisition (n=30) were illustrated in box plots and evaluated based on independent samples t-test, ***P <0.001.

(5)

Global levels of H3K4me³, RNAPII, and H3K9me³show altered patterns in ccRCC

To reveal whether H3K4me³level, a widely applied euchromatic mark, correlates with the progression of ccRCC, immunohistochemistry staining was performed on the dissected tissues from all the 30 patients who were classified according to stages (1–4) and grades (low and high).

To perform comparable immunohistochemical detection on non- tumorous and ccRCC parts, following microtome dissection of the tissues, the segments obtained from either the tumour mass or the surrounding normal part of the tissue were placed on the same microscopic slide as it is shown in Supplementary Fig. 1. Subsequently, in case of each patient, the staining protocol was performed simultaneously on the non-tumorous and tumorous tissue parts. Regardless of tumour grades, the H3K4me³level significantly increased in 50% of the cohort, while we detected less H3K4me³in 5 tumorous specimens (Fig. 2 and Sup- plementary Fig. 2). Although the H3K4me³level generally altered in 63% of the samples, we did not find any correlation with cancer progression considering the patients’ tumour stages and grades (Table 1).

We also monitored the H3K4me³status and the cancer survival rate, and we found that the alteration of this PTM resulted in better survivability (Supplementary Fig. 3). Based on these, we believe this marker presumably has a prognostic value which has to be further validated on an extended clinical cohort.

To study whether the alteration in the euchromatic H3K4me³level correlates with the transcriptional profile changes of these patients, we performed co-immunofluorescence staining with H3K4me³and RNAPII antibodies. Henceforth, the RNAPII levels of the tumorous and adjacent normal tissue parts were analysed in all the examined ccRCC patients simultaneously with H3K4me³. These results revealed a significantly strong correlation between the alteration of H3K4me³and RNAPII levels, which presumes a transcriptomic profile change in these tumours

(Fig. 3 and Supplementary Fig. 4). In 14 samples increased levels, while in 3 specimens a decreased level of RNAPII was detected in the tumorous part compared to the corresponding adjacent normal tissue part. The elevated level of RNAPII mainly concerns high grade ccRCC in these patients (10 versus 4 samples).

We hypothesized that the increased RNAPII binding and the elevated level of H3K4me³might contribute to an imbalance in the euchromatic–heterochromatic milieu of these cells. To decipher this issue, H3K9me³heterochromatic mark was selected and detected by the same immunostaining procedure described before. We observed altered H3K9me³level in 20 patients: decreased in 11 cases, while increased in 9 samples (Fig. 4 and Supplementary Fig. 5).

Higher rate of DNA damage is detected in ccRCC

As mentioned previously, ccRCC is characterised as an aggressive tumour type, in which high genomic mutation frequency is peculiar [25]. It also affects certain tumour suppressor genes, including VHL and TP53, which play a significant role in DNA repair [26]. However, the severity of DNA damage and the maintenance of genome stability in ccRCC have yet to be characterised. Since the metastatic capability of the tumour is notably influenced by the type of DNA damage and mutation frequency, we explored the possible alterations in the ccRCC-related DNA repair process. To follow the DNA repair foci formation in primary ccRCC tumours, we applied immunostaining with an antibody against one of the initial DDR factors, γH2A.X (Fig. 5 and Supplementary Fig. 6). We observed increased γH2A.X levels in 50% of the specimens irrespectively of tumour grading, which underlined our hypothesis that persistent DNA damage occurs in these cells, which is presumably one of the driving forces of the metastatic progression in ccRCC cells. Additionally, in ccRCC patients, we revealed significant correlations between γH2A.X and H3K4me³or H3K9me³, which support Fig. 3. Immunofluorescent spectra of RNAPII detected in normal and ccRCC tumorous tissue parts. A) N represents non-tumorous tissue parts, while T corresponds to ccRCC sections. Representative images of RNAPII staining in peritumoral (left) and tumour tissues (right) from ccRCC patients. For the better visualization of the changes in RNAPII level, one cell is represented in a magnified resolution in each specimen. Scale bars represent 30 µm and 1.5 µm, respectively. B) The observed alterations in the quantification of relative fluorescence intensity (CTCF) of RNAPII level between the tumorous and non-tumorous tissue parts obtained from automated image acquisition (n=30) were illustrated in box plots and evaluated based on independent samples t-test, ***P <0.001.

(6)

6

Fig. 4. Immunofluorescent spectra of H3K9me³detected in normal and ccRCC tumorous tissue parts. A) N represents non-tumorous tissue parts, while T corresponds to ccRCC sections. Representative images of H3K9me³staining in peritumoral (left) and tumour tissues (right) from ccRCC patients. For the better visualization of the changes in H3K9me³level, one cell is represented in a magnified resolution in each specimen. Scale bars represent 30 µm and 1.5 µm, respectively. B) The observed alterations in the quantification of relative fluorescence intensity (CTCF) of H3K9me³level between the tumorous and non-tumorous tissue parts obtained from automated image acquisition (n=30) were illustrated in box plots and evaluated based on independent samples t-test, ***P <0.001.

Fig. 5. Immunofluorescent spectrum of γH2A.X (H2AXS139P) detected in normal and ccRCC tumorous tissue parts. A) N represents non-tumorous tissue parts, while T corresponds to ccRCC sections. Representative images of γH2A.X staining in peritumoral (left) and tumour tissues (right) from ccRCC patients. For the better visualization of the changes in γH2A.X level, one cell is represented in a magnified resolution in each specimen. Scale bars represent 30 µm and 1.5 µm, respectively.

B) The observed alterations in the quantification of relative fluorescence intensity (CTCF) of γH2A.X level between the tumorous and non-tumorous tissue parts obtained from automated image acquisition (n=30) were illustrated in box plots and evaluated based on independent samples t-test, ***P <0.001.

(7)

the assumption that chromatin structural changes, occurred in the tumorous regions, associate with persistent DNA damage which is more prevalent in tumorous regions than in the adjacent non-tumorous parts (Fig. 6).

Discussion

ccRCC is considered a highly invasive tumour type, associated with poor prognosis and high chance of recurrence [27]. To diagnose these patients, a classical pathological scoring system is still widely used to predict the prognosis of ccRCC; however, it mainly relies on anatomical information and has no biological characteristics [28]. Additionally, mutations have been recently used to characterise the status of VHL, PBRM1, SETD2, and BAP1 genes [17,19]. On the other hand, due to the alteration in histone post-translational modifications (PTMs), chromatin structural changes have been associated with carcinogenesis and turned to be potential key regulators of cancer-related pathways [29,30].

Dysregulation in histone PTMs can be originated from the malfunction of either PTM writers or erasers and also from the genome-wide transcriptional reprogramming occurring during tumorigenesis. Addition- ally, DNA methylation can be used to characterise epigenetic changes in ccRCC, and the methylome dataset can be applied to predispose prognostic signatures [31]. This emphasizes the biological and diagnostic relevance of the epigenetic marks used to characterise ccRCC tumours regarding their invasiveness and response to the applied therapy. In this study, we characterised tumorous and adjacent normal tissues dissected from 30 patients diagnosed with ccRCC, and we validated the signatures of H3K4me³, and RNAPII associated with actively transcribed genomic regions. Although previous studies described an increased level of H3K4me³during ccRCC progression, they found a decrease in the level of this PTM at later stages of ccRCC [32]. In contrast to this finding, we detected a synergistic elevation both in H3K4me³and RNAPII level which confirms the reliability of our data. Our findings are further supported by the fact that extensive transcriptional reprogramming and overactivation occur during tumour progression in various types of cancer. The present study also establishes a strong correlation between H3K4me³and RNAPII marks. Finally, the sample processing and automated analysis we applied in our experimental setup also minimise the appearance of experimental bias. Our data is in accordance with a

previous study that established a correlation between the level of H3K4me³and later tumour stages, indicating the prognostic implication of this marker [33]. In accordance with this, our study also supports the idea that global histone modifications, including H3K4me³, are also altered during ccRCC progression.

As it was mentioned before, the altered DNA methylation pattern is also a hallmark of cancerous malformations, i.e. DNA hypermethylation contributes to heterochromatinization which leads to the appearance of H3K9me³mark at the given genomic locus, which then recruits Het- erochromatic Protein 1 (HP1), resulting in constitutive heterochromatin formation [34]. We also found that the alteration in the global level of H3K9me³corresponding with changes in the level of H3K4me³and RNAPII was correlated with the presence of ccRCC. This finding is in accordance with already published data which demonstrated that patients diagnosed with ccRCC exhibited low level of heterochromatin-related H3K9 mono- and di-methylations as well as increased level of euchromatin-related H3K9ac and H3K18ac [35–37].

Taken together, these data reveal that the global level of different histone modifications is altered in unfavourably localized ccRCC. The correlation between histone methylation and acetylation presumes epigenetic confusion in ccRCC cells due to aberrant activity and/or expression of histone-modifying enzymes. These findings are further supported by our previous study, in which we demonstrated that the PBRM1, BAP1, and SETD2 expression levels were completely altered in patients developing ccRCC (accepted manuscript at POR), and it was also published that Jumonji histone demethylases are also dysregulated in ccRCC patients [38]. Despite the correlation between histone methylation and acetylation, these modifications may have different relevance to ccRCC stages [39]. Thus, epigenetic alterations may represent a novel therapeutic potential: Histone Deacetylase (HDAC) inhibitors show an antineoplastic effect on ccRCC patients as valproic acid inhibits the growth of RCC tumours and vorinostat increases the anticancer activity of the novel therapeutic drugs [40]. Additionally, the globally altered histone methylation levels detected in ccRCC can be restored with histone demethylase inhibitors, such as Lysine-specific histone demethylase 1 (LSD1) [41]. We believe that our study has the potential in combining clinicopathological variables with methylation signatures, by which the current risk classification of non-metastatic ccRCC can be improved.

Fig. 6. Bivariate Pearson correlation calculated between H3K4me³, RNAPII, H3K9me³and γH2A.X is represented, *P <0.05 and **P <0.01.

(8)

8 Finally, in ccRCC tumour-derived specimens, we observed increased γH2A.X level, which is the hallmark of persistence DNA damage [42,43].

In correlation with the perpetual presence of γH2A.X, in ccRCC patients considerable number of DNA damage or insufficient DNA repair takes place. In agreement with this, previous studies revealed that during ccRCC progression, the DNA repair mechanism is also malfunctional [44]. First, the application of PARP and ATR inhibitors on PBRM1-deficient cells, is known to be synthetically lethal [45]. In the absence of exogenous DNA damage, PBRM1-defective ccRCC-derived cells exhibited elevated levels of replication stress and micronuclei formation. It unveils that the loss of PBRM1 during ccRCC progression is strongly associated with genome instability [45]. Our findings also supported this idea as insufficient DNA repair occurs in these cells, which opens the possibility of using PARPi or other chemotherapeutic drugs targeting DNA repair in ccRCC-related tumour therapy. Second, a large number of differentially expressed DNA repair genes were in silico identified, which reflect on the patient outcomes in the development of ccRCC [46]. These data also support our findings that the improper function of DNA repair genes is probably the consequence of epigenetic changes which contribute to an aggravated mutagenic rate in ccRCC and largely affect the tumour progression. Hence, these markers can be used to predict the clinical benefits and efficacy of the applied therapy.

Funding

This research was funded by the National Research, Development and Innovation Office grant GINOP-2.2.1-15-2017-00052 and NKFI-FK 132080. T.P. was funded by National Research, Development and Innovation Office grant GINOP-2.2.1-15-2017-00052, the János Bolyai Research Scholarship of the Hungarian Academy of Sciences BO/27/20, UNKP-20-5-SZTE-265, by NKFI-FK 132080, UNKP-21-5-SZTE-563 and EMBO short-term fellowship 8513. L.K. was funded by the University of Szeged, Faculty of Medicine Research Fund-Hetényi Géza Grant (Grant No. 5S 340 A202) and the New National Excellence Programme (Grant No. UNKP-21-4-SZTE-131).

Author contribution

Study design, N.O., B.N.B., H.M., L.K., and T.P.; investigation, N.¨ O., ¨ Z.U.; data collections, N.O., B.N.B., H.M., Z.U., G.P.B., S.B., F.S., L.K., T. ¨ P.; data analysis, N.O., B.N.B., U.Z., T.P.; writing—original draft prep-¨ aration, N.O., T.P.; writing¨ —review and editing, N.O., B.N.B., H.M., Z. ¨ U., G.P.B., S.B., F.S., L.K., T.P.; visualization, N.O., B.N.B., Z.U.; funding ¨ acquisition, B.N.B., L.K., T.P. All authors have read and agreed to the published version of the manuscript.

Declaration of competing interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

Acknowledgment

We thank Ivett Csik´os for her expert technical help.

Supplementary materials

Supplementary material associated with this article can be found, in the online version, at doi:10.1016/j.tranon.2022.101420.

References

[1] O. Bratu, D. Mischianu, D. Marcu, D. Spinu, L. Iorga, A. Cherciu, I. Balescu, N. Bacalbasa, C. Diaconu, C. Savu, C. Savu, R. Anghel, Renal tumor biomarkers (Review), Exp. Ther. Med. 22 (2021) 1297.

[2] A. Thorstenson, M. Bergman, A.H. Scherman-Plogell, S. Hosseinnia, B. Ljungberg, J. Adolfsson, S. Lundstam, Tumour characteristics and surgical treatment of renal cell carcinoma in Sweden 2005-2010: a population-based study from the national Swedish kidney cancer register, Scand. J. Urol. 48 (2014) 231–238.

[3] S.R. Prasad, P.A. Humphrey, J.R. Catena, V.R. Narra, J.R. Srigley, A.D. Cortez, N.

C. Dalrymple, K.N. Chintapalli, Common and uncommon histologic subtypes of renal cell carcinoma: imaging spectrum with pathologic correlation, Radiographics 26 (2006) 1795–1806, discussion 1806-1710.

[4] J.R. Srigley, B. Delahunt, Uncommon and recently described renal carcinomas, Mod. Pathol. 22 (2) (2009) S2–S23. Suppl.

[5] J.R. Srigley, B. Delahunt, J.N. Eble, L. Egevad, J.I. Epstein, D. Grignon, O. Hes, H. Moch, R. Montironi, S.K. Tickoo, M. Zhou, P. Argani, I.R.T. Panel, The International Society of Urological Pathology (ISUP) Vancouver Classification of Renal Neoplasia, Am. J. Surg. Pathol. 37 (2013) 1469–1489.

[6] N. Goossens, S. Nakagawa, X. Sun, Y. Hoshida, Cancer biomarker discovery and validation, Transl. Cancer Res 4 (2015) 256–269.

[7] N.J. Farber, C.J. Kim, P.K. Modi, J.D. Hon, E.T. Sadimin, E.A. Singer, Renal cell carcinoma: the search for a reliable biomarker, Transl. Cancer Res. 6 (2017) 620–632.

[8] G.S. Bhattacharyya, D.C. Doval, C.J. Desai, H. Chaturvedi, S. Sharma, S.

P. Somashekhar, Overview of breast cancer and implications of overtreatment of early-stage breast cancer: an indian perspective, JCO Glob. Oncol. 6 (2020) 789–798.

[9] E.M.C. Tacconi, M. Tuthill, A. Protheroe, Review of adjuvant therapies in renal cell carcinoma: evidence to date, Onco. Targets Ther. 13 (2020) 12301–12316.

[10] T. Chandrasekar, Z. Klaassen, H. Goldberg, R.K. Sayyid, G.S. Kulkarni, N.

E. Fleshner, High competing risks minimize real-world utility of adjuvant targeted therapy in renal cell carcinoma: a population-based analysis, Oncotarget 9 (2018) 16731–16743.

[11] Q.K. Li, C.P. Pavlovich, H. Zhang, C.R. Kinsinger, D.W. Chan, Challenges and opportunities in the proteomic characterization of clear cell renal cell carcinoma (ccRCC): A critical step towards the personalized care of renal cancers, Semin.

Cancer Biol. 55 (2019) 8–15.

[12] R. Mayeux, Biomarkers: potential uses and limitations, NeuroRx 1 (2004) 182–188.

[13] S.A. Byrnes, B.H. Weigl, Selecting analytical biomarkers for diagnostic applications: a first principles approach, Expert Rev. Mol. Diagn. 18 (2018) 19–26.

[14] H. Moch, J. Srigley, B. Delahunt, R. Montironi, L. Egevad, P.H. Tan, Biomarkers in renal cancer, Virchows. Arch. 464 (2014) 359–365.

[15] E. Kim, S. Zschiedrich, Renal Cell Carcinoma in von Hippel-Lindau Disease-From Tumor Genetics to Novel Therapeutic Strategies, Front. Pediatr. 6 (2018) 16.

[16] H. Kim, B.Y. Shim, S.J. Lee, J.Y. Lee, H.J. Lee, I.H. Kim, Loss of Von Hippel-Lindau (VHL) tumor suppressor gene function: VHL-HIF pathway and advances in treatments for metastatic Renal Cell Carcinoma (RCC), Int. J. Mol. Sci. 22 (2021).

[17] S. Bihr, R. Ohashi, A.L. Moore, J.H. Ruschoff, C. Beisel, T. Hermanns, A. Mischo, C. Corro, J. Beyer, N. Beerenwinkel, H. Moch, P. Schraml, Expression and mutation patterns of PBRM1, BAP1 and SETD2 mirror specific evolutionary subtypes in clear cell renal cell carcinoma, Neoplasia 21 (2019) 247–256.

[18] F. Piva, M. Santoni, M.R. Matrana, S. Satti, M. Giulietti, G. Occhipinti, F. Massari, L. Cheng, A. Lopez-Beltran, M. Scarpelli, G. Principato, S. Cascinu, R. Montironi, BAP1, PBRM1 and SETD2 in clear-cell renal cell carcinoma: molecular diagnostics and possible targets for personalized therapies, Expert Rev. Mol. Diagn. 15 (2015) 1201–1210.

[19] J. Brugarolas, PBRM1 and BAP1 as novel targets for renal cell carcinoma, Cancer J.

19 (2013) 324–332.

[20] A.A. de Cubas, W.K. Rathmell, Epigenetic modifiers: activities in renal cell carcinoma, Nat. Rev. Urol. 15 (2018) 599–614.

[21] Y. Cheng, C. He, M. Wang, X. Ma, F. Mo, S. Yang, J. Han, X. Wei, Targeting epigenetic regulators for cancer therapy: mechanisms and advances in clinical trials, Signal Transduct. Target Ther. 4 (2019) 62.

[22] T.K. Kelly, D.D. De Carvalho, P.A. Jones, Epigenetic modifications as therapeutic targets, Nat. Biotechnol. 28 (2010) 1069–1078.

[23] J.C. Angulo, C. Manini, J.I. Lopez, A. Pueyo, B. Colas, S. Ropero, The role of epigenetics in the progression of clear cell renal cell carcinoma and the basis for future epigenetic treatments, Cancers (Basel) 13 (2021).

[24] B.N. Borsos, I. Huliak, H. Majoros, Z. Ujfaludi, A. Gyenis, P. Pukler, I.M. Boros, T. Pankotai, Human p53 interacts with the elongating RNAPII complex and is required for the release of actinomycin D induced transcription blockage, Sci. Rep.

7 (2017) 40960.

[25] S. Wu, X. Li, A genomic instability-derived risk index predicts clinical outcome and immunotherapy response for clear cell renal cell carcinoma, Bioengineered 12 (2021) 1642–1662.

[26] C.M. Robinson, M. Ohh, The multifaceted von Hippel-Lindau tumour suppressor protein, FEBS Lett. 588 (2014) 2704–2711.

[27] J.J. Hsieh, M.P. Purdue, S. Signoretti, C. Swanton, L. Albiges, M. Schmidinger, D.

Y. Heng, J. Larkin, V. Ficarra, Renal cell carcinoma, Nat. Rev. Dis. Primers 3 (2017) 17009.

[28] B. Delahunt, J.K. McKenney, C.M. Lohse, B.C. Leibovich, R.H. Thompson, S.

A. Boorjian, J.C. Cheville, A novel grading system for clear cell renal cell carcinoma incorporating tumor necrosis, Am. J. Surg. Pathol. 37 (2013) 311–322.

[29] J.E. Audia, R.M. Campbell, Histone Modifications and Cancer, Cold Spring Harb.

Perspect. Biol. 8 (2016), a019521.

[30] T. Pankotai, O. Komonyi, L. Bodai, Z. Ujfaludi, S. Muratoglu, A. Ciurciu, L. Tora, J. Szabad, I. Boros, The homologous Drosophila transcriptional adaptors ADA2a and ADA2b are both required for normal development but have different functions, Mol. Cell. Biol. 25 (2005) 8215–8227.

(9)

[31] B.N. Lasseigne, J.D. Brooks, The role of DNA methylation in renal cell carcinoma, Mol Diagn Ther 22 (2018) 431–442.

[32] A. Kumar, N. Kumari, U. Sharma, S. Ram, S.K. Singh, N. Kakkar, K. Kaushal, R. Prasad, Reduction in H3K4me patterns due to aberrant expression of methyltransferases and demethylases in renal cell carcinoma: prognostic and therapeutic implications, Sci. Rep. 9 (2019) 8189.

[33] S. Li, L. Shen, K.N. Chen, Association between H3K4 methylation and cancer prognosis: a meta-analysis, Thorac. Cancer 9 (2018) 794–799.

[34] L. Monaghan, M.E. Massett, R.P. Bunschoten, A. Hoose, P.A. Pirvan, R.M.

J. Liskamp, H.G. Jorgensen, X. Huang, The emerging role of H3K9me3 as a potential therapeutic target in acute myeloid leukemia, Front. Oncol. 9 (2019) 705.

[35] D. Minardi, G. Lucarini, A. Filosa, G. Milanese, A. Zizzi, R.D. Primio, R. Montironi, G. Muzzonigro, Prognostic role of global DNA-methylation and histone acetylation in pT1a clear cell renal carcinoma in partial nephrectomy specimens, J. Cell. Mol.

Med. 13 (2009) 2115–2121.

[36] D. Mosashvilli, P. Kahl, C. Mertens, S. Holzapfel, S. Rogenhofer, S. Hauser, R. Buttner, A. Von Ruecker, S.C. Muller, J. Ellinger, Global histone acetylation levels: prognostic relevance in patients with renal cell carcinoma, Cancer Sci. 101 (2010) 2664–2669.

[37] E. Cocco, M. Leo, C. Canzonetta, S. Di Vito, A. Mai, D. Rotili, A. Di Napoli, A. Vecchione, C. De Nunzio, P. Filetici, A. Stoppacciaro, KAT3B-p300 and H3AcK18/H3AcK14 levels are prognostic markers for kidney ccRCC tumor aggressiveness and target of KAT inhibitor CPTH2, Clin Epigenetics 10 (2018) 44.

[38] S. Ramadoss, X. Chen, C.Y. Wang, Histone demethylase KDM6B promotes epithelial-mesenchymal transition, J. Biol. Chem. 287 (2012) 44508–44517.

[39] S. Ramakrishnan, L. Ellis, R. Pili, Histone modifications: implications in renal cell carcinoma, Epigenomics 5 (2013) 453–462.

[40] S. Ramakrishnan, R. Pili, Histone deacetylase inhibitors and epigenetic modifications as a novel strategy in renal cell carcinoma, Cancer J. 19 (2013) 333–340.

[41] L. Zhu, J. Wang, W. Kong, J. Huang, B. Dong, Y. Huang, W. Xue, J. Zhang, LSD1 inhibition suppresses the growth of clear cell renal cell carcinoma via upregulating P21 signaling, Acta Pharm Sin B 9 (2019) 324–334.

[42] J. Ji, Y. Zhang, C.E. Redon, W.C. Reinhold, A.P. Chen, L.K. Fogli, S.L. Holbeck, R.

E. Parchment, M. Hollingshead, J.E. Tomaszewski, Q. Dudon, Y. Pommier, J.

H. Doroshow, W.M. Bonner, Phosphorylated fraction of H2AX as a measurement for DNA damage in cancer cells and potential applications of a novel assay, PLoS One 12 (2017), e0171582.

[43] T. Pankotai, C. Bonhomme, D. Chen, E. Soutoglou, DNAPKcs-dependent arrest of RNA polymerase II transcription in the presence of DNA breaks, Nat. Struct. Mol.

Biol. 19 (2012) 276–282.

[44] E. Guo, C. Wu, J. Ming, W. Zhang, L. Zhang, G. Hu, The clinical significance of DNA damage repair signatures in clear cell renal cell carcinoma, Front Genet 11 (2020), 593039.

[45] R.M. Chabanon, D. Morel, T. Eychenne, L. Colmet-Daage, I. Bajrami, N. Dorvault, M. Garrido, C. Meisenberg, A. Lamb, C. Ngo, S.R. Hopkins, T.I. Roumeliotis, S. Jouny, C. Henon, A. Kawai-Kawachi, C. Astier, A. Konde, E. Del Nery, C. Massard, S.J. Pettitt, R. Margueron, J.S. Choudhary, G. Almouzni, J.C. Soria, E. Deutsch, J.A. Downs, C.J. Lord, S. Postel-Vinay, PBRM1 deficiency confers synthetic lethality to DNA repair inhibitors in cancer, Cancer Res. 81 (2021) 2888–2902.

[46] D.J. Clark, S.M. Dhanasekaran, F. Petralia, J. Pan, X. Song, Y. Hu, F. da Veiga Leprevost, B. Reva, T.M. Lih, H.Y. Chang, W. Ma, C. Huang, C.J. Ricketts, L. Chen, A. Krek, Y. Li, D. Rykunov, Q.K. Li, L.S. Chen, U. Ozbek, S. Vasaikar, Y. Wu, S. Yoo, S. Chowdhury, M.A. Wyczalkowski, J. Ji, M. Schnaubelt, A. Kong, S. Sethuraman, D.M. Avtonomov, M. Ao, A. Colaprico, S. Cao, K.C. Cho, S. Kalayci, S. Ma, W. Liu, K. Ruggles, A. Calinawan, Z.H. Gumus, D. Geiszler, E. Kawaler, G.C. Teo, B. Wen, Y. Zhang, S. Keegan, K. Li, F. Chen, N. Edwards, P.M. Pierorazio, X.S. Chen, C.

P. Pavlovich, A.A. Hakimi, G. Brominski, J.J. Hsieh, A. Antczak, T. Omelchenko, J. Lubinski, M. Wiznerowicz, W.M. Linehan, C.R. Kinsinger, M. Thiagarajan, E.

S. Boja, M. Mesri, T. Hiltke, A.I. Robles, H. Rodriguez, J. Qian, D. Fenyo, B. Zhang, L. Ding, E. Schadt, A.M. Chinnaiyan, Z. Zhang, G.S. Omenn, M. Cieslik, D.W. Chan, A.I. Nesvizhskii, P. Wang, H. Zhang, C. Clinical proteomic tumor analysis, integrated proteogenomic characterization of clear cell renal cell carcinoma, Cell, 179 (2019) 964–983, e931.