# Loss of chromosome Y leads to down regulation of KDM5D and KDM6C epigenetic modifiers in clear cell renal cell carcinoma

Madeleine Arseneault, Jean Monlong, Naveen S Vasudev, Ruhina S Laskar, Maryam Safisamghabadi, Patricia Harnden, Lars Egevad, Nazanin Nourbehesht, Pudchalaluck Panichnantakul, Ivana Holcatova, et al.

► To cite this version:

Madeleine Arseneault, Jean Monlong, Naveen S Vasudev, Ruhina S Laskar, Maryam Safisamghabadi, et al.. Loss of chromosome Y leads to down regulation of KDM5D and KDM6C epigenetic modifiers in clear cell renal cell carcinoma. Scientific Reports, 2017, 7 (1), pp.44876. 10.1038/srep44876. cea-04505246

HAL Id: cea-04505246

https://cea.hal.science/cea-04505246

Submitted on 14 Mar 2024

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

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

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

Distributed under a Creative Commons Attribution 4.0 International License
SCIENTIFIC REPORTS

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

OPEN

# Loss of chromosome Y leads to down regulation of KDM5D and KDM6C epigenetic modifiers in clear cell renal cell carcinoma

Received: 27 September 2016

Accepted: 15 February 2017

Published: 23 March 2017

Madeleine Arseneault<sup>1,2,*</sup>, Jean Monlong<sup>1,2,*</sup>, Naveen S. Vasudev<sup>3</sup>, Ruhina S. Laskar<sup>4</sup>, Maryam Safisamghabadi<sup>1,2</sup>, Patricia Harnden<sup>3</sup>, Lars Egevad<sup>5</sup>, Nazanin Nourbehesht<sup>1,2</sup>, Pudchalaluck Panichnantakul<sup>1,2</sup>, Ivana Holcatova<sup>6</sup>, Antonin Brisuda<sup>7</sup>, Vladimir Janout<sup>8</sup>, Helena Kollarova<sup>8</sup>, Lenka Foretova<sup>9</sup>, Marie Navratilova<sup>9</sup>, Dana Mates<sup>10</sup>, Viorel Jinga<sup>11</sup>, David Zaridze<sup>12</sup>, Anush Mukeria<sup>12</sup>, Pouria Jandaghi<sup>1,2</sup>, Paul Brennan<sup>4</sup>, Alvis Brazma<sup>13</sup>, Jorg Tost<sup>14</sup>, Ghislaine Scelo<sup>4</sup>, Rosamonde E. Banks<sup>3</sup>, Mark Lathrop<sup>1,2</sup>, Guillaume Bourque<sup>1,2</sup> & Yasser Riazalhosseini<sup>1,2</sup>

Recent genomic studies of sporadic clear cell renal cell carcinoma (ccRCC) have uncovered novel driver genes and pathways. Given the unequal incidence rates among men and women (male:female incidence ratio approaches 2:1), we compared the genome-wide distribution of the chromosomal abnormalities in both sexes. We observed a higher frequency for the somatic recurrent chromosomal copy number variations (CNVs) of autosomes in male subjects, whereas somatic loss of chromosome X was detected exclusively in female patients (17.1%). Furthermore, somatic loss of chromosome Y (LOY) was detected in about 40% of male subjects, while mosaic LOY was detected in DNA isolated from peripheral blood in 9.6% of them, and was the only recurrent CNV in constitutional DNA samples. LOY in constitutional DNA, but not in tumor DNA was associated with older age. Amongst Y-linked genes that were downregulated due to LOY, *KDM5D* and *KDM6C* epigenetic modifiers have functionally-similar X-linked homologs whose deficiency is involved in ccRCC progression. Our findings establish somatic LOY as a highly recurrent genetic defect in ccRCC that leads to downregulation of hitherto unsuspected epigenetic factors, and suggest that different mechanisms may underlie the somatic and mosaic LOY observed in tumors and peripheral blood, respectively.

<footnote n="1"><sup>1</sup>Department of Human Genetics, McGill University, 1205 Dr Penfield Avenue, Montreal, QC, H3A 0B1, Canada.</footnote>

<footnote n="2"><sup>2</sup>McGill University and Genome Quebec Innovation Centre, 740 Doctor Penfield Avenue, Montreal, QC, H3A 0G1, Canada.</footnote>

<footnote n="3"><sup>3</sup>Leeds Institute of Cancer and Pathology, University of Leeds, Cancer Research Building, St James's University Hospital, Leeds, LS9 7TF, UK.</footnote>

<footnote n="4"><sup>4</sup>International Agency for Research on Cancer (IARC), 150 cours Albert Thomas, 69008 Lyon, France.</footnote>

<footnote n="5"><sup>5</sup>Karolinska Institutet, Department of Pathology, SE-171 77 Stockholm, Sweden.</footnote>

<footnote n="6"><sup>6</sup>First Faculty of Medicine, Institute of Hygiene and Epidemiology, Charles University in Prague, Studničkova 7, Praha 2, 128 00 Prague, Czech Republic.</footnote>

<footnote n="7"><sup>7</sup>University Hospital Motol, V Úvalu 84, 150 06 Prague, Czech Republic.</footnote>

<footnote n="8"><sup>8</sup>Department of Preventive Medicine, Faculty of Medicine, Palacky University, Hnevotinska 3, 775 15 Olomouc, Czech Republic.</footnote>

<footnote n="9"><sup>9</sup>Department of Cancer Epidemiology and Genetics, Masaryk Memorial Cancer Institute and MF MU, Zluty Kopec 7, 656 53 Brno, Czech Republic.</footnote>

<footnote n="10"><sup>10</sup>National Institute of Public Health, Dr Leonte Anastasievici 1-3, sector 5, Bucuresti 050463, Romania.</footnote>

<footnote n="11"><sup>11</sup>Carol Davila University of Medicine and Pharmacy, Th. Burghele Hospital, 20 Panduri Street, 050659 Bucharest, Romania.</footnote>

<footnote n="12"><sup>12</sup>Russian N.N. Blokhin Cancer Research Centre, Kashirskoye shosse 24, Moscow 115478, Russian Federation.</footnote>

<footnote n="13"><sup>13</sup>European Molecular Biology Laboratory, European Bioinformatics Institute, EMBL-EBI, Wellcome Trust Genome Campus, Hinxton, CB10 1SD, UK.</footnote>

<footnote n="14"><sup>14</sup>Laboratory for Epigenetics & Environment, Centre National de Génotypage, CEA-Institut de Génomique, 2 rue Gaston Crémieux, 91000 Evry, France.</footnote>

<footnote>*These authors contributed equally to this work. Correspondence and requests for materials should be addressed to Y.R. (email: yasser.riazalhosseini@mcgill.ca)</footnote>
Chromosomal aneuploidy is a common phenomenon in many cancers, and the analysis of copy number varia-
tions (CNVs) across multiple samples has helped identify relevant driver genes for human cancers. For example,
several oncogenes including *MYC*, *EGFR*, *ERBB2* and *CCND1* are recurrently amplified through chromosomal
or focal gains, while multiple tumor suppressors such as *ATM*, *PTEN* and *CDKN2A* are commonly deleted in
different cancers¹.

Clear cell renal cell carcinoma (ccRCC), which accounts for 75–80% of all renal cell carcinomas, is charac-
terized by loss of chromosome 3p in about 90% of the sporadic cases². Remarkably, 3p harbors the four most
commonly mutated genes in ccRCC whose cancer-driving activities have been established in the disease; *VHL*³,
*PBRM1*⁴, *SETD2*⁵, and *BAP1*⁶, which are mutated in 80%, 40%, 19% and 12% of cases, respectively⁷⁻⁹. Inactivation
of VHL leads to constitutive stabilization of the hypoxia inducible transcription factors (HIF), and abnormal
activation of their downstream genes, which contribute to cancer development¹⁰. The remaining three genes
encode proteins involved in chromatin remodeling and histone modifications, highlighting the important role of
epigenome aberration in the disease². While the incidence of ccRCC is increasing worldwide, the male-to-female
incidence ratios are typically within the range of 1.5–2:1.0¹¹, arguing for a sex-specific analysis of the genomic
abnormalities. Here, we set out to investigate the occurrence and the extent of germline and somatic CNVs in spo-
radic ccRCC in male and female patients separately, and to further characterize those affecting sex chromosomes.

**Results and Discussion**

**Loss of chromosome Y is common in ccRCC.** Using whole-genome sequencing (WGS) data of ccRCC and matched constitutional DNA sample pairs, which we have reported recently¹², we interrogated CNVs in DNA from 52 male and 41 female patients (discovery set; Supplementary Table 1) by analyzing coverage of sequencing reads mapped to each chromosome (see Methods). In line with previous literature, the most frequent somatic CNV was the loss of 3p detected in 91% of samples, followed by recurrent gains of chromosomes 5q (32%), 7 (23.6%), 12 (13%), and losses of chromosomes 14q (30%), 8p (29%) and 9 (16%). Overall, tumors from male patients exhibited higher prevalence for the recurrent chromosomal aberrations, in particular for gain of 7q (28% in males vs. 17% in females) and deletion of 9p (25% in males vs. 10% in females) (Fig. 1a). In contrast, we observed that loss of chromosome X (LOX) exclusively happens in female patients (17.1% of female cases). Given that several X-linked genes escape X-inactivation, and have therefore two functional copies in females but one in males, this observation suggests that presence of a copy of chromosome X may potentially be essential for the survival of cancer cells. Curiously, whereas no tumors from male patients displayed LOX, loss of chromosome Y (LOY) was the second most frequent somatic chromosome aneuploidy in these tumors (36.5% of male subjects, N = 19; Fig. 1a). The fraction of cells estimated to be affected by somatic LOY in these patients ranged from 11% to 75%, and in 14 patients somatic LOY was detected in at least 20% of the cells (Fig. 1b). Next, we examined the presence of CNVs in constitutional DNA isolated from peripheral samples collected from the same patients. Of significance, LOY was the only recurrent aneuploidy in constitutional DNA of our samples that was detected in 5 male patients (9.6%; Supplementary Figure 1), of which 4 showed the deletion in more than 20% of cells (Fig. 1b). Corroborating previous studies¹³,¹⁴, the observed LOY was associated with older age in patients (P = 0.04); the average age of patients with LOY in the peripheral blood was 68.9 year in comparison to 58.8 year in those without this abnormality. Notably, we did not observe any association between age of patients and extent of somatic LOY in tumors of the affected patients.

**LOY is a whole-chromosome event.** Given the high prevalence of LOY in tumors and peripheral DNA of male patients, we further analyzed LOY in our sample series, particularly whether the observed LOY spans the whole chromosome or is focal. Analysis of sequencing read coverage along chromosome Y showed that the loss is observed throughout the chromosome in samples affected by LOY (Fig. 2a), suggesting that the deletion affects the whole chromosome. Based on availability of DNA, we subjected samples from seven of the patients affected with somatic LOY to verification by an orthogonal Y-Chromosome deletion detection assay surveying the presence of twenty specific regions of the Y chromosome by polymerase chain reaction (PCR) (see Methods). Somatic LOY at the chromosomal level was confirmed in all examined tumors, evident from an attenuated amplification of Y-chromosome-specific loci in DNA isolated from tumor samples compared to that of the matched constitutional DNA. This pattern was not observed in samples of other male patients who had not been identified as being affected by somatic LOY based on the analysis of their WGS data (Fig. 2b).

To confirm these findings, we screened tumor and matched control DNA sample pairs of an additional
48 male ccRCC patients (validation set) for LOY using the above PCR-based assay. This analysis revealed
somatic LOY in 20 (42.7%) of the validation sample set, demonstrating that this is a common genomic aber-
ration in ccRCC, detected in 39.6% overall (discovery and validation sets; n = 100) of male ccRCC patients
(Supplementary Figure 2). Analysis of association between somatic LOY and clinical annotations including
tumor stage or grade did not show any significant relationships.

**LOY results in downregulation of epigenetic modifier genes.** We further examined the possible
effect of somatic LOY at the RNA level by interrogating a RNA-Seq dataset on gene expression in normal and
tumor samples from male patients within the discovery set¹². We found that 11 genes had significantly differ-
ent patterns of expression in tumors of the patients with and without somatic LOY (false-discovery rate (FDR)
<0.01; Supplementary Table 2). These 11 genes were located on chromosome Y, and while expressed in nor-
mal kidney tissue, exhibited lower expression in tumors of patients harboring somatic LOY, indicating that this
aberration may have functional consequences through deregulation of the affected genes. Moreover, the level of
expression of each gene was found to be inversely correlated to the proportion of cells affected by LOY (Fig. 3).
This observation was confirmed using gene expression data generated by microarrays, which was available for
29 tumors of the validation set¹⁵ (Supplementary Figure 3). We surveyed the list of genes affected for potential
<figure><img src="image_3.png" /><figcaption>**Figure 1.** Copy number analysis in ccRCC. (a) Bar graphs show the frequency of copy number variations across the genome in ccRCC tumors. Frequencies are presented in samples from female and male cases separately. (b) Status of chromosome Y in DNA isolated from tumors (Y-axis) and patient-matched peripheral blood (X-axis) is shown for individual male subjects. In samples affected by LOY, the normalized coverage of chromosome Y, shown on Y and X axes for tumor and normal samples, respectively, is lower than the expected value of 0.5. The color codes define patient groups with different states for LOY.</figcaption></figure>

functionally-relevant candidates. Among these genes, *TMSB4Y* has recently been identified as a tumor suppressor gene downregulated in male breast cancers¹⁶, but not connected to ccRCC. Likewise, deletion of *KDM5D* has been detected in 52% of prostate cancers¹⁷. *KDM5D* encodes a lysine-specific histone H3 demethylase, which plays an important role in epigenetic regulation¹⁸. Furthermore, it has been shown that knockdown of KDM5D through RNA-interference (RNAi) increases cell proliferation and reduces apoptosis in prostate cancer¹⁹, suggesting a tumor suppressor function for this gene. Intriguingly, *KDM5C*, the X-linked homologue of *KDM5D* is recurrently mutated in ccRCC⁸,⁹,¹², and its inactivation leads to genomic instability in ccRCC through deregulation of H3K4 methylation²⁰. *KDM5D* shows 85% sequence identity to *KDM5C*, and the products of these two genes possess a similar function in demethylating tri-methyl H3K4¹⁸,²⁰. Given this functional similarity, we surveyed the mutational status of *KDM5C* in our discovery set, and investigated possible relationships between mutational status of *KDM5C* and *KDM5D* in tumors of male patients. In female patients, *KDM5C* was deleted in tumors of 7 cases through somatic LOX, and was affected by focal somatic deletions in two additional patients. Furthermore, somatic mutations of *KDM5C* were present in tumors of 3 patients who were also affected with LOX (P = 0.003, Fisher's exact test, Supplementary Figure 4). Overall, *KDM5C* was affected with somatic genomic aberrations in 9 out of 41 (22%) female cases. As *KDM5C* escapes the X-inactivation²¹, the concomitant mutations of *KDM5C* and LOX in the same tumors may suggest that this gene is a classical tumor suppressor affected with bi-allelic inactivation in ccRCC. In male cases, we identified *KDM5C* mutations in tumors of 3 patients (5.8%), of which one was also affected by somatic LOY (Supplementary Figure 4). We did not detect any mutation or a focal CNV affecting *KDM5D* in tumors of the male patients who did not exhibit LOY.

Our list of LOY-associated down-regulated genes (Supplementary Table 2) includes another epigenome modifier with an X-linked homologue that is also recurrently mutated in ccRCC; *UTY/KDM6C*. KDM6C demethylates
<figure><img src="image_4.png" /><figcaption>**Figure 2.** LOY affects whole chromosome. (a) Sequencing coverage across chromosome Y is shown in constitutional DNA samples without (top) and with LOY (middle), and in a tumor sample with LOY (bottom). (b) The cartoon on top depicts the location of the loci examined by PCR on Y chromosome. The dot graph on bottom shows average of relative amplification values (Tumor/normal samples of the same patient) for each locus in patients with (red) and without (blue) somatic LOY (SLOY). Error bars show the range across patients of each group.</figcaption></figure>

H3K27, a function similar to that of KDM6A²². These genes also share over 83% in sequence similarity, resulting in highly conserved active sites in their products. Mutations of KDM6A leading to its inactivation have been recurrently observed in ccRCC⁸,²³, highlighting this gene as a potential key tumor suppressor in renal cancer. In addition to being affected by somatic LOX in 7 female patients, KDM6A was also affected by focal deletion in a female patient in our cohort.

**KDM5D expression reduces viability of renal cancer cells.** Given the reported tumor-suppressive function of KDM5D in prostate cancer¹⁹, and of its X-link homolog KDM5C in renal cancer²⁰, we set out to examine whether KDM5D expression has an anti-tumor activity in renal cancer. We first evaluated KDM5D expression levels in several renal cancer cell lines, which have been derived from tumors resected from male patients. Amongst cell lines examined, ACHN cell line did not show any expression for KDM5D (Fig. 4a). This observation was in line with a previous study reporting the loss of chromosome Y in ACHN cell line²⁴. We therefore selected this cell line for functional analysis of KDM5D expression. Ectopic expression of KDM5D cells reduced cell viability to 65% as compared to control transfection (Fig. 4b,c), suggesting the potential involvement of KDM5D depletion in renal cancer pathology.

**Conclusions**

Emerging data emphasizes an association between LOY in peripheral blood and higher risk of cancer²⁵. Likewise focal or chromosome-level somatic LOY occurs recurrently in different malignancies; however, current knowledge of mechanisms by which LOY may contribute to cancer is limited. Recent genomic studies of ccRCC have highlighted the importance of molecular aberrations that impair the function of chromatin remodeling and epigenetic modifiers in ccRCC development⁵,²⁰,²⁶⁻²⁹. Our study expands these findings by highlighting the prevalence of somatic LOY among men affected by ccRCC, and suggesting a functional relevance for this aberration through down-regulation of previously unrecognized epigenetic modifiers *KDM5D* and *KDM6C*. Given the functional
<figure><img src="image_5.png" /><figcaption>**Figure 3.** Somatic LOY leads to downregulation of Y-linked genes. Expression of Y chromosome genes downregulated in patients affected by somatic LOY is compared to the proportion of cells estimated to harbor somatic LOY in individual tumor samples.</figcaption></figure>

similarities between these genes and their X-linked homologs, it is plausible that down-regulation of *KDM5D* and *KDM6C*, through somatic LOY, may contribute to ccRCC development or progression. Our *in vitro* data shows that over expression of KDM5D in cancer cells that are affected by LOY reduces cell viability. These findings indicate that down-regulation of KDM5D through LOY may contribute to the pathogenesis of renal cancer. However, further detailed analysis through future functional studies is warranted to understand the exact function and pathway context of KDM5D in renal cancer.

## Methods

**Patient samples and DNA isolation.** Clinical information for patients included in this study is presented in Supplementary Table 1. Patients undergoing nephrectomy for suspected renal cancer during the period December 2008 to March 2011 at St James's University Hospital in Leeds, UK; University Hospital Motol, Prague, Czech Republic; Masaryk Memorial Cancer Institute, Brno, Czech Republic; Th. Burghele Hospital, Bucharest, Romania; and N. N. Blokhin Cancer Research Centre, Moscow, Russia, were recruited to the study after informed consent was obtained. Recruitment in Central and Eastern Europe was coordinated by the International Agency for Research on Cancer (IARC). All experiments and methods were performed in accordance to the ethics guidelines from the International Cancer Genome Consortium (ICGC) and to the relevant national regulations and with sampling and clinical data collection being undertaken according to predefined standard operating procedures (SOPs) based on guidelines from ICGC. Ethical approvals were obtained from the Leeds (East) Local Research Ethics Committee, the IARC Ethics Committee, as well as from local ethics committee for recruiting centers in Czech Republic, Romania, and Russia. DNA from fresh-frozen tumor tissue samples and buffy coat was isolated using Autopure (Qiagen) as described previously¹², and were quantified by Quant-iT PicoGreen dsDNA Assay Kit (Invitrogen, ON, CAN).

**Inference of LOY from WGS data.** WGS data of tumor and blood DNA samples studied here were reported previously¹². To detect aneuploidy and LOY from WGS data, we first measured read coverage across the genome in 5 Kbp bins. In each sample, the coverage was normalized by the median coverage across the autosomes. We then estimated, for each sample, the median normalized coverage in each chromosome arm. The only exception was chromosome Y which was considered as a whole. In order to avoid noise due to mapability issues, we used only the top 1000 bins with the lowest median divergence from the expected baseline in the normal samples. We used this normalized median coverage per chromosome arm to test aneuploidy in each sample. For each chromosome arm (or chromosome Y), a mixture of two Gaussian distributions was fitted to the empirical distribution of the median normalized coverage across samples. The main Gaussian was used as the
<figure><img src="image_6.png" /></figure>

<figure><img src="image_7.png" /><figcaption>**Figure 4.** Effect of KDM5D on viability of renal cancer cells. (a) Expression levels of KDM5D mRNA in renal cancer cell lines derived from tumors procured from male patients, as measured by qRT-PCR. GAPDH served as a housekeeping gene for measurement of relative gene expression. (b) Over expression of KDM5D in ACHN cell line reduces cell viability. Values are the mean ± SD of six independent experiments. \*\*P < 0.01 when compared to the corresponding results from control (ctrl) (Mann-Whitney U test). (c) Over expression of KDM5D following transfection was confirmed using qRT-PCR.</figcaption></figure>

null distribution (Supplementary Figure 5) to derive P-values. A chromosome arm was flagged as aneuploid if the Bonferroni-adjusted P-value was smaller than 0.01 and at least 10% of cells were affected. The proportion of cell with aneuploidy was estimated as the proportion of missing/excess coverage. For LOY, we expect a normalized coverage of 0.5 and the proportion of cells with LOY was (0.5-coverage)/0.5.

We used a logistic regression to test the association of LOY with age. Finally, the CNVs used for *KDM5C* or *KDM6A* deletion investigation were detected by PopSV³⁰ using the normal samples as reference and 5 Kbp bins.

**PCR-based detection of LOY.** To examine the status of LOY in DNA of tumor and blood samples, Y Chromosome Deletion Detection System assay, Version 2 (Promega, WI, USA) was used as instructed by the manufacturer. Briefly, 20 specific regions of the Y chromosome were amplified by PCR using 5 multiplex master mixes, and PCR products were loaded on a QIAxcel instrument (Qiagen, ON, CAN). Densities of PCR products were estimated by BioCalculator software (v.3.2) and a normalization was performed by the control primer pair included in each multiplex master mix to control the amplification efficacy. We also included samples from three male subjects without LOY and one female sample to control the performance of the assay. Similar to the analysis on WGS data, the probes were first normalized by the median probe amplification value across the normal samples. Then the median of the normalized amplification was computed for each sample. It summarized the overall amplification of chromosome Y in each sample. These values were used to produce Supplementary Figure 2 and to identify LOY. Following the same analysis as for the WGS data, the mixture of Gaussian distributions was fitted on the normalized amplification of the normal samples. Samples which deviated significantly (*P* < 0.01) from the expected amplification and with an estimated proportion of affected of cells >10% were flagged as being affected by LOY.

**Gene expression analysis.** Transcriptome profiles of the tumor samples included in this study (previously reported in our earlier publication¹²), were used to examine differential gene expression between male subjects affected with somatic LOY and those without this abnormality. RNA-seq data was available for tumors of 34 patients, of which 21 had RNA-seq for matched normal kidney samples. Differentially expressed genes between tumors affected with somatic LOY and those without this abnormality were identified using Student's T-Test on log2-transformed RPKM data, and the Benjamini-Hochberg method was used to correct multiple testing. Genes with a FDR< 0.01 were considered differentially expressed. A linear regression was used to test the association between the proportion of cells with somatic LOY and gene expression (RPKM).
Gene expression microarray data for 29 tumors of validation samples had previously been reported¹⁵, and were used to confirm the anti-correlation between the proportion of cells with somatic LOY and gene expression levels (log2 intensity).

**Cell viability assay.** Renal cancer cell lines 786-O, A704, Caki-2, ACHN were obtained from ATCC (Rockville, USA) and cultured in RPMI, EMEM and McCoy medium supplemented with 10% (v/v) fetal bovine serum (FBS), 100 U/ml penicillin and 100 lg/ml streptomycin. Cells were incubated at 37°C and 5% (v/v) CO2. For viability assays, 5000 cells were transfected with 100 ng of either KDM5D cDNA-expressing (courtesy of Dr. Stephane Richard) or control empty vector (Sigma, Oakville, Canada) in 96-well plates using Lipofectamine 3000 (Invitrogen) according to the manufacturer's instructions. CellTiter-Glo assay (Promega, WI, USA) was used to assess cell viability after 72 hours post-transfection.

**Quantitative real-time PCR (qRT-PCR).** Total RNA was extracted from cells using miRNeasy kit (Qiagen, Toronto, Canada) according to the supplier protocols. 1 µg RNA was reverse transcribed into complementary DNA (cDNA) using Transcription First Strand cDNA Synthesis Kit (Roche, Laval, Canada) following instructions provided by the manufacturer. Real-time PCR reactions were prepared using LightCycler 480 SYBR green I master kit (Roche), and were run on a LightCycler 480 instrument (Roche) according to the manufacturer's recommendations. Triplicate PCR reactions were performed for each sample to ensure reliability. Expression of *KDM5D* mRNA was normalized to the expression of the housekeeping gene *GAPDH*, and was reported as $2^{-\Delta Ct}$. All the primers were purchased from IDT (Coralville, IA, US). The sequences of primers were CGTGGAAGGACTCATGACCA (GAPDH forward), GCCATCACGCCACAGTTTC (GAPDH reverse), CGCAGCTTTGAAGAGCTAAG (KDM5D forward) and CAGCTGTGGAGTGTCCATCC (KDM5D reverse).

**References**

1. Zack, T. I. et al. Pan-cancer patterns of somatic copy number alteration. *Nat Genet* **45**, 1134–40 (2013).

2. Frew, I. J. & Moch, H. A clearer view of the molecular complexity of clear cell renal cell carcinoma. *Annu Rev Pathol* **10**, 263–89 (2015).

3. Latif, F. et al. Identification of the von Hippel-Lindau disease tumor suppressor gene. *Science* **260**, 1317–20 (1993).

4. Benusiglio, P. R. et al. A germline mutation in PBRM1 predisposes to renal cell carcinoma. *J Med Genet* **52**, 426–30 (2015).

5. Carvalho, S. et al. SETD2 is required for DNA double-strand break repair and activation of the p53-mediated checkpoint. *Elife* **3**, e02482 (2014).

6. Popova, T. et al. Germline BAP1 mutations predispose to renal cell carcinomas. *Am J Hum Genet* **92**, 974–80 (2013).

7. Varela, I. et al. Exome sequencing identifies frequent mutation of the SWI/SNF complex gene PBRM1 in renal carcinoma. *Nature* **469**, 539–42 (2011).

8. Dalgliesh, G. L. et al. Systematic sequencing of renal carcinoma reveals inactivation of histone modifying genes. *Nature* **463**, 360–3 (2010).

9. Comprehensive molecular characterization of clear cell renal cell carcinoma. *Nature* **499**, 43–9 (2013).

10. Harris, A. L. Hypoxia-a key regulatory factor in tumour growth. *Nat Rev Cancer* **2**, 38–47 (2002).

11. Capitanio, U. & Montorsi, F. Renal cancer. *Lancet* **387**, 894–906 (2016).

12. Scelo, G. et al. Variation in genomic landscape of clear cell renal cell carcinoma across Europe. *Nat Commun* **5**, 5135 (2014).

13. Pierre, R. V. & Hoagland, H. C. Age-associated aneuploidy: loss of Y chromosome from human bone marrow cells with aging. *Cancer* **30**, 889–94 (1972).

14. Zhou, W. et al. Mosaic loss of chromosome Y is associated with common variation near TCL1A. *Nat Genet* **48**, 563–8 (2016).

15. Wozniak, M. B. et al. Integrative genome-wide gene expression profiling of clear cell renal cell carcinoma in Czech Republic and in the United States. *PLoS One* **8**, e57886 (2013).

16. Wong, H. Y. et al. TMSB4Y is a candidate tumor suppressor on the Y chromosome and is deleted in male breast cancer. *Oncotarget* **6**, 44927–40 (2015).

17. Perinchery, G. et al. Deletion of Y-chromosome specific genes in human prostate cancer. *J Urol* **163**, 1339–42 (2000).

18. Lee, M. G., Norman, J., Shilatifard, A. & Shiekhattar, R. Physical and functional association of a trimethyl H3K4 demethylase and Ring6a/MBLR, a polycomb-like protein. *Cell* **128**, 877–87 (2007).

19. Jangravi, Z. et al. Two Splice Variants of Y Chromosome-Located Lysine-Specific Demethylase 5D Have Distinct Function in Prostate Cancer Cell Line (DU-145). *J Proteome Res* **14**, 3492–502 (2015).

20. Rondinelli, B. et al. Histone demethylase JARID1C inactivation triggers genomic instability in sporadic renal cancer. *J Clin Invest* **125**, 4625–37 (2015).

21. Agulnik, A. I. et al. A novel X gene with a widely transcribed Y-linked homologue escapes X-inactivation in mouse and human. *Hum Mol Genet* **3**, 879–84 (1994).

22. Walport, L. J. et al. Human UTY(KDM6C) is a male-specific N-methyl lysyl demethylase. *J Biol Chem* **289**, 18302–13 (2014).

23. van Haaften, G. et al. Somatic mutations of the histone H3K27 demethylase gene UTX in human cancer. *Nat Genet* **41**, 521–3 (2009).

24. du Manoir, S. et al. Detection of complete and partial chromosome gains and losses by comparative genomic in situ hybridization. *Hum Genet* **90**, 590–610 (1993).

25. Forsberg, L. A. et al. Mosaic loss of chromosome Y in peripheral blood is associated with shorter survival and higher risk of cancer. *Nat Genet* **46**, 624–8 (2014).

26. Kakarougkas, A. et al. Requirement for PBAF in transcriptional repression and repair at DNA breaks in actively transcribed regions of chromatin. *Mol Cell* **55**, 723–32 (2014).

27. Simon, J. M. et al. Variation in chromatin accessibility in human kidney cancer links H3K36 methyltransferase loss with widespread RNA processing defects. *Genome Res* **24**, 241–50 (2014).

28. Pfister, S. X. et al. SETD2-dependent histone H3K36 trimethylation is required for homologous recombination repair and genome stability. *Cell Rep* **7**, 2006–18 (2014).

29. Riazalhosseini, Y. & Lathrop, M. Precision medicine from the renal cancer genome. *Nat Rev Nephrol* **12**, 655–666 (2016).

30. Monlong, J. et al. Human copy number variants are enriched in regions of low-mappability. *bioRxiv* (2016).

**Acknowledgements**

This study was supported by funding from EU FP7 under grant agreement number 241669 (the CAGEKID project, http://www.cng.fr/cagekid), and funding from Génome Québec as well as Ministère de l'Enseignement supérieur, de la Recherche, de la Science et de la Technologie (MESRST). We also acknowledge the support of Cancer Research UK Centre and ECMC infrastructure funding in Leeds for contributions to sample collection.
The Czech-Brno group was supported by MH CZ - DRO (MMCI, 00209805). JM is supported by funds from National Sciences and Engineering Research Council (NSERC-448167-2013).

**Author Contributions**

Y.R. conceived the study and designed the experiments with contribution from G.B. N.V., P.H., L.E., I.H., A.B., V.J., H.K., L.F., M.N., D.M., Vi.J., D.Z., A.M., P.B., G.S. and R.E.B. were responsible for patient selection, sample collection, sample preparation and pathological reviews. M.A. and N.N. prepared DNA, performed experiments to detect LOY and analyzed the data with contribution from P.J. J.M. developed the computational methods for WGS analysis, analyzed WGS data and generated the figures. P.P., R.S.L., P.J. and Al.B. performed the gene expression analysis. M.L. provided critical advice on data analysis and statistical approaches. M.S. and P.J. performed functional analysis with cell line models. M.A. and Y.R. wrote the manuscript, with assistance from J.T., G.S. and R.E.B.

**Additional Information**

**Supplementary information** accompanies this paper at http://www.nature.com/srep

**Competing Interests:** The authors declare no competing financial interests.

**How to cite this article:** Arseneault, M. et al. Loss of chromosome Y leads to down regulation of KDM5D and KDM6C epigenetic modifiers in clear cell renal cell carcinoma. *Sci. Rep.* 7, 44876; doi: 10.1038/srep44876 (2017).

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

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

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/

© The Author(s) 2017