Next Article in Journal
Role of the NUDT Enzymes in Breast Cancer
Next Article in Special Issue
Comparative Transcriptomes of the Body Wall of Wild and Farmed Sea Cucumber Isostichopus badionotus
Previous Article in Journal
Identification of Enzymatic Bottlenecks for the Aerobic Production of Malate from Glycerol by the Systematic Gene Overexpression of Anaplerotic Enzymes in Escherichia coli
Previous Article in Special Issue
Analysis of Stop Codons within Prokaryotic Protein-Coding Genes Suggests Frequent Readthrough Events
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 22
Issue 5
10.3390/ijms22052265
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Somatic Copy Number Alterations and Associated Genes in Clear-Cell Renal-Cell Carcinoma in Brazilian Patients

by
Flávia Gonçalves Fernandes
1,
Henrique Cesar Santejo Silveira
1,
João Neif Antonio Júnior
2,
Rosana Antunes da Silveira
1,
Luis Eduardo Zucca
2,
Flavio Mavignier Cárcano
2,3,
André Octavio Nicolau Sanches
2,
Luciano Neder
4,
Cristovam Scapulatempo-Neto
4,
Sergio Vicente Serrano
2,3,
Eric Jonasch
5,
Rui Manuel Reis
1,6,7,* and
Adriane Feijó Evangelista
1,*
1
Molecular Oncology Research Center, Barretos Cancer Hospital, Barretos 14784-400, Brazil
2
Department of Medical Oncology, Barretos Cancer Hospital, Barretos 14784-400, Brazil
3
Barretos School of Health Sciences Dr Paulo Prata-FACISB, Barretos 14785-002, Brazil
4
Department of Pathology, Barretos Cancer Hospital, Barretos 14784-400, Brazil
5
Department of Genitourinary Medical Oncology, Division of Cancer Medicine, The University of Texas MD Anderson Cancer Center, Houston, TX 77030, USA
6
Life and Health Sci Research Institute (ICVS), Medical School, University of Minho, 4710-057 Braga, Portugal
7
ICVS/3B’s-PT Government Associate Laboratory, 4710-057 Braga/Guimarães, Portugal
*
Authors to whom correspondence should be addressed.
Int. J. Mol. Sci. 2021, 22(5), 2265; https://doi.org/10.3390/ijms22052265
Submission received: 23 October 2020 / Revised: 13 January 2021 / Accepted: 22 January 2021 / Published: 25 February 2021
(This article belongs to the Special Issue Molecular Mechanisms of Gene Expression: “Bioinformatics of Gene Regulations and Structure”)
Browse Figures
Versions Notes

Abstract

:
Somatic copy number aberrations (CNAs) have been associated with clear-cell renal carcinoma (ccRCC) pathogenesis and are a potential source of new diagnostic, prognostic and therapeutic biomarkers. Recurrent CNAs include loss of chromosome arms 3p, 14q, 9p, and gains of 5q and 8q. Some of these regional CNAs are suspected of altering gene expression and could influence clinical outcomes. Despite many studies of CNAs in RCC, there are currently no descriptions of genomic copy number alterations in a Brazilian ccRCC cohort. This study was designed to evaluate the chromosomal profile of CNAs in Brazilian ccRCC tumors and explore clinical associations. A total of 92 ccRCC Brazilian patients that underwent nephrectomy at Barretos Cancer Hospital were analyzed for CNAs by array comparative genomic hybridization. Most patients in the cohort had early-stage localized disease. The most significant alterations were loss of 3p (87.3%), 14q (35.8%), 6q (29.3%), 9p (28.6%) and 10q (25.0%), and gains of 5q (59.7%), 7p (29.3%) and 16q (20.6%). Bioinformatics analysis revealed 19 genes mapping to CNA significant regions, including SETD2, BAP1, FLT4, PTEN, FGFR4 and NSD1. Moreover, gain of 5q34-q35.3 (FLT4 and NSD1) and loss of 6q23.2-q23.3 (MYB) and 9p21.3 (MLLT3) had gene expression levels that correlated with TCGA data and was also associated with advanced disease features, such as larger tumors, Fuhrman 3, metastasis at diagnosis and death. The loss of region 14q22.1 which encompasses the NIN gene was associated with poor overall survival. Overall, this study provides the first CNA landscape of Brazilian patients and pinpoints genomic regions and specific genes worthy of more detailed investigations. Our results highlight important genes that are associated with copy number changes involving large chromosomal regions that are potentially related to ccRCC tumorigenesis and disease biology for future clinical investigations.

1. Introduction

Renal-cell carcinoma (RCC) comprises a complex set of diseases with different histopathological subtypes, clinical courses, and therapeutic responses [1,2,3,4]. Kidney cancer represented 2–3% of all cancers, affecting over 403,000 new individuals and accounted for around 175,000 cancer-related deaths worldwide in 2018 [5]. The 5-year survival rate for localized kidney cancer is up to 90%, but this rate drops dramatically for patients with distant metastasis (11.9%), accounting for 30% of the advanced disease diagnosed at presentation [5,6]. The clear-cell subtype is the most common histological subtype of renal-cell carcinoma (70–85%) and accounts for the most RCC-specific deaths [7]. In Brazil, the estimated number for 2018 was 10,688 new cases and 4084 deaths [5]. At Barretos Cancer Hospital, which receives patients from all regions in Brazil, urologic tumors accounted for 12.7% of all cancer cases in 2017, and over 150 new cases of kidney cancer were diagnosed in the same year (1.4% of all cases). More than half of the patients were diagnosed at stages II and IV, with a 5-year survival rate of 56% and 10%, respectively (https://infogram.com/rhc_hcb).
Recurrent chromosomal abnormalities and copy number imbalances are related to clear-cell renal-cell carcinoma (ccRCC) tumorigenesis [8]. Most of the ccRCC copy number alterations (CNAs) consist of extensive chromosomal losses or gains that can involve multiple genes, in which changes of gene dosage may alter the levels of gene’s expression [9,10,11,12,13]. The impact of CNAs has been recently associated with changes in the metabolic activity since copy number losses of KEAP1 are involved in the transcriptional control of glutathione metabolism genes [14]. The comparison of the CNAs profile with the germline alterations in single cells from VHL/PBRM1-negative ccRCCs showed higher amplification of these alterations, suggesting an important role of CNAs in the progression of these tumors [15]. Indeed, recurrent CNAs have been shown to interfere with some of the affected gene’s functions [11,12,13], but there is presently no established association with specific genes affected by CNAs and the ccRCC tumor biology or disease progression. Importantly, CNAs may constitute useful diagnostic, therapeutic and prognostic biomarkers [16,17].
Loss of chromosome 3p is the major cytogenetic hallmark of ccRCC [18,19]. The loss or inactivation of the von Hippel–Lindau (VHL) tumor-suppressor gene at 3p25.3 is considered the crucial step in ccRCC carcinogenesis; however, other alterations are also required [4,10]. The Cancer Genome Atlas (TCGA) consortium has characterized and made publicly available over 20,000 primary cancer and matched normal samples from 33 cancer types, including ccRCC and other kidney cancer types [10]. The TCGA datasets include genomic, transcriptomic, and protein level biological information that has validated the most frequent alterations in ccRCC, providing new insights about potential therapeutic targets, such as PI(3)K/AKT pathway and SWI/SNF chromatin remodeling complex [10]. Losses of chromosomes 14q, 9p and deletions on 6q and 8p are also reported and have been described to be correlated with advanced tumors and worse patient survival [17,20,21,22]. Gain of chromosome 5q is the only CNA associated with better patient survival [17]. Understanding the role of discrete CNAs and its impacts on gene dosage are essential to better characterize the ccRCC pathogenesis, and may be helpful in future clinical settings.
Despite multiples studies on RCC CNAs, there are no current reports of chromosomal copy number alterations in a Brazilian ccRCC cohort and their possible clinicopathological correlations. Also, data regarding the relationship between CNA and molecular aspects of RCC in Brazil are still scarce [23]. Therefore, this study aimed to evaluate the CNA profile in a series of Brazilian ccRCC patients and explore putative clinicopathological associations related to CNA-driven genes’ expression as well as describe new potential molecular features and possible molecular targets that could be worthy for further investigation.

2. Results

2.1. Study Population—Clinicalpathological Features

The main clinicopathological features of the present ccRCC series are summarized in Table 1. There were slightly more male patients than females, with a mean age of 58 years old (ranging from 20 to 81 years). More than half of the patients were white and presented localized disease (stages I and II). Approximately 20% of the ccRCC patients had metastasis at diagnosis and received systemic therapy. The mean follow-up time was four years and almost 70% of the patients were alive without disease. Disease progression and cancer-specific death was seen in 19.5% of the cases (Table 1). The major risk factors for ccRCC seen in this cohort were hypertension (54.3%), smoking (33.6%) and diabetes (22.8%) (Table S1). The 5-year overall survival rate for ccRCC cancer patients was 77.3%, with an estimated median survival time of 53.5 months. Kaplan–Meier plots showed that OS was significantly lower for patients in advanced stages (p = 0.01) and higher histopathological grade (p < 0.001), metastasis at diagnosis (p = 0.002), and tumors larger than 7 centimeters (p = 0.003) (Table S1, Figure S1 and Figure S2).

2.2. Analysis of Somatic Copy Number Alterations and Integration with In Silico Transcriptomic Profile

The somatic CNA landscape is shown in Figure 1. A median of 67.5 alterations per patient tumor was found. An average of 10–12% of each patient’s tumor genome had copy number imbalances. STAC analysis identified 158 regions of interest (Table S2). Most significant regions included those related to loss of 3p (87.3%), 14q (35.8%), 9p (28.6%) and 10q (25.0%), and gains of 5q (59.7%), 6q (29.3%), 7p (29.3%) and 16q (20.6%) (Figure 1).
Bioinformatics analysis using Nexus Copy Number software disclosed 19 genes mapped to CNA significant regions (Table 2). Compared to TCGA data, there was a modest difference in overall CNA frequencies.
According to Cancer Genome Interpreter (CGI) analysis, alterations in SETD2, BAP1, FLT4 and PTEN genes were reported as known drivers in several tumor types, and events in FGFR4 and NSD1 genes were classified as predicted drivers. Next, to look for putative correlations between CNAs and gene’s expression profile, the 19 genes were investigated for gains or losses of expression in TCGA dataset through the Oncomine platform (504 ccRCC samples) (504 ccRCC samples). Losses of BAP1, MYB, MLLT3 and TFG showed concurrent underexpression. Likewise, gains of CDH11 and FLT4 genes showed overexpression in TCGA dataset (Figure 2).

2.3. CNAs and Gene Expression Clinical Impact and Survival Analysis

The significant regions mentioned above were further investigated for their clinical associations (Table 3). Except for 3p loss, all other alterations, i.e., deletions on 9p (9p-) and 14q (14q-), and gain on 7p (7p+) were significantly associated with features of advanced disease, such as tumor metastasis, T3/T4 stages, higher Fuhrman grade (3–4), larger tumors and dismal outcome (death).
The multivariate analysis showed that the gain in the 5q has been shown to increase the risk of large tumors (over 7 centimeters) (p = 0.03; 95% CI = 1.08–8.14). Patients with a gain of 7p21.2 are more likely to be diagnosed with tumors that affect adjacent structures (T3 and T4) (p = 0.02; 95% CI = 1.11–8.04) and in advanced histological degrees (p = 0.01; 95% CI = 1.28–9.10). Patients with loss of 9p21.3 are more likely to have advanced staging (p = 0.005; 95% CI = 1.51–10.20) and with metastatic disease (p = 0.003; 95% IC = 1.72–16.02). Finally, patients with a loss of 14q22.1 had an increased risk of death approximately 3 times compared to patients without this event (p = 0.02; 95% CI = 1.12–7.84) (Table 4).
In addition, loss of chromosome 14q, which encompasses the NIN gene was associated with shorter survival (p = 0.03) (Figure 3). Interestingly, 5q gain (5q+) was also associated with poor outcomes (Table 3).
We further investigated the putative associations between CNA burden (total number of CNAs and PGA) and clinicopathological features (Table 5). Patients harboring tumors with more than 50 CN alterations (25th percentile) presented larger tumor size (p = 0.002), poor outcome (death) (p = 0.04) and higher Fuhrman grades (p = 0.01) in the multivariate analysis (Table 5). However, when using PGA for multivariate analysis, it was not significative, i.e., tumor size (≥7 cm), outcome (death) and Fuhrman grade showed p = 0.08, p = 0.33 and p = 0.13, respectively.
We next investigated whether some genes of interest that were disclosed in our study could be also altered in the TCGA dataset through Oncomine platform (Table 6). According to TCGA data, loss of BAP1 in 3p, an oncogenic aberration known to be a biomarker of RCC was only associated with Fuhrman grade 1 (p = 0.02). Losses of MYB and MLLT3, despite being predicted as passengers alterations were significantly associated with advanced disease: Loss of MLLT3 gene was associated with metastasis at diagnosis (p = 0.005), stage IV disease (p = 0.02), larger tumors (p = 0.01) and poor outcomes (p = 0.02), while loss of MYB was correlated with Fuhrman grade 3 (p = 0.02). Additionally, gains on 5q genes showed significant association with increased cancer aggressiveness. Loss of FTL4 was also associated with larger tumors (p = 0.007) and loss of NSD1, considered a predicted driver aberration was associated with larger tumors (p = 0.007) and advanced disease (stage IV; p = 0.01). Aberrations of TFG and CDH11, despite showing changes on gene expression levels are still reported as predicted passenger alterations and did not show any clinicopathological associations.

3. Discussion

Copy number alterations (CNA) are important drivers of ccRCC tumorigenesis [4]. This study provides the first comprehensive report of CNA profiles in Brazilian ccRCC patients. We observed significant numbers of CNA, including loss of 3p (87.3%), 14q (35.8%), 6q (29.3%), 9p (28.6%) and 10q (25.0%), and gains of 5q (59.7%), 7p (29.3%) and 16q (20.6%). Some regions, including the gain of 5q34–q35.3 (FLT4 and NSD1) and loss of 6q23.2–q23.3 (MYB), 9p21.3 (MLLT3), and 14q22.1 (NIN) had gene expression levels validated by in silico re-analysis of TCGA data and were associated with patients clinicopathologic features.
Overall, the frequency of CNA in chromosomal arms and whole chromosomes reveal similar profiles to those demonstrated by others, especially the TCGA [10]. We then re-analyzed TCGA aCGH array data with the same analysis parameter settings used for Brazilian patients. The CNA profile in Brazilian patients differs from TCGA only for chromosomes 6, 9, and 10. A possible explanation is that tumor events take place over time, and some important drivers, in general, occur at the beginning of the carcinogenesis process [24]. It has been reported that 3p loss is a truncal event in RCC [25,26]. Furthermore, the gain of 5q and loss of 14q are usually shared by different patient cohorts and are likely clonal events that occur early in the ccRCC disease process [25,26]. In our data, both changes are recurrent and similar to TCGA. Despite few differences, variations across populations have been reported in some cancer types. Recently, variations in VHL and PBRM1 mutation rates were observed according to ancestry, with fewer mutations in African kidney cancer patients [27]. Moreover, African American kidney cancer types from TCGA data harbor a lower level of chromosomal instability [28]. Further studies on ccRCC mutations in the Brazilian population and local ancestry inference are necessary to estimates the number of copies at genomic sites to unravel such differences. Additionally, loss of 3p in our patient cohort was similar between metastatic and non-metastatic tumors (88% vs. 84%, respectively). Although it is associated with initial histological grade by univariate analysis, there was an absence of multivariate significance and no correlation with survival. Since most of the samples in our study are primary tumors of clinical stages I/II, lower frequencies of chromosome imbalance would be expected. One limitation of our cohort is that it may not be possible to identify relevant aberrations specifically related to metastatic disease or advanced disease.
Gunawan and colleagues suggest that a gain of 5q is a clinically favorable alteration in ccRCC with an increase in overall survival rate [17]. Moreover, gains of 5q were further reported with a frequency twice as high in non-metastatic tumors than in metastatic ones [29]. However, 5q gain was not associated with significant differences in overall survival in our cohort.
Furthermore, the 14q loss, associated with loss of HIF1A, is another recurrent alteration in ccRCC and has been predicted to drive more aggressive disease [10]. Several studies showed a significant correlation of alterations in this region with a worse prognosis, distant metastasis and decreased overall survival in accordance with our findings [16,21,30,31]. It was also demonstrated that the loss of chromosome 14 leads to a decrease in HIF1-a protein levels since HIF1A gene maps at 14q23.2 [30]. The 14q22.1 region which encompasses the NIN gene was associated with decreased overall survival and was also identified in our analysis of the TCGA data. In our findings, 77% of patients with loss of NIN also showed concomitant loss of HIF1A. The NIN gene has been initially described as playing a role to centrosomal function [32]. The involvement of this gene in cancer is unclear. Future studies are warranted to clarify the role of NIN gene in the ccRCC tumorigenesis.
Concerning additional cancer genes identified in the significantly altered regions, some were characterized as known drivers (SETD2, BAP1, FLT4 and PTEN) and others as predicted drivers (FGFR4 and NSD1), according to the CGI. The CGI tool identifies as driver genes according to their mode of action in tumorigenesis (i.e., whether they function as oncogenes or tumor suppressors), based on either experimentally verified sources or in silico prediction [33]. However, despite the losses of MYB and MLLT3 were considered to be predicted passenger by CGI analysis, our findings show concordance with gene expression and association with advanced disease. Another important known driver is the loss PBRM1, found in 83% of our cohort despite not being identified in significant regions by the analysis of STAC. Deletions of PTEN, SETD2 and BAP1 genes, are described as biomarkers for RCC and are responsive to drugs of preclinical and clinical trials [34,35,36,37,38,39]. Also, deletions of JAK2, CD274 and MYD88 genes, and amplification of FGFR4 and NPM1 genes, are also involved in clinical trials for other tumor types [40,41,42,43]. Interestingly, alterations in JAK2, MYD88 and NPM1 genes, already have drugs FDA approved for other cancers [44,45,46]. Additionally, the loss of NSD1 gene, a possible driver alteration associated with tumors ≥7 cm and advanced stages. NSD1 gene located at chromosome 5, which encodes for the nuclear receptor-binding SET-domain containing protein 1 (NSD1). Little is known regarding the function of NSD1, but its genomic location is near FGFR4. Since the NSD1 gene alteration was not validated in silico in our gene expression analysis, it is possible that other genes in the region could play a driver role in such region in ccRCC.
Altogether, it is essential to consider that changes in CNAs could not affect gene expression levels due to transcriptional adaptative mechanisms. A recent study used a new method, known as transcriptional adaptation to CNA profiling (TACNA), and found a strong correlation of expression levels changes by CNAs with mean methylation from TCGA datasets [47]. These findings suggest that dysregulation of the epigenome could be one of the underlying mechanisms driving the translational response to CNAs, and further studies could help to understand the impact of the results regarding tumor suppressors involved in the chromatin modification mechanisms. Herein, we also addressed the poorly explored relationship between the overall CNA load and clinical data in ccRCC. Arai and colleagues in an unsupervised hierarchical cluster analysis based on aCGH results were also able to separate patients into 2 groups according to the CNA pattern [48]. In our current study, the group of patients with less than 50 CNAs was significantly associated with evidence of localized diseases such as small tumors, early stages and histological grade, absence of metastases, and better outcomes. Increased CNAs has been reported as indicative of chromosomal instability and have been associated with poor cancer outcome [49]. Despite a high number of studies associating tumor mutation burden with immune cell infiltrates, studies reporting the impact of CNA load in predict the response to targeted therapies and the outcome are scarce. Identifying biomarkers could better select treatments and improve survival in patients with metastatic ccRCC and CNAs could provide new approaches for identifying potential subtypes and pathways suitable for targeted therapies and improved outcome.

4. Materials and Methods

4.1. Sample Collection

This study comprises a retrospective cohort of 92 treatment-naïve patients diagnosed with clear-cell renal-cell carcinoma (ccRCC) that underwent nephrectomy between 2008 and 2014 at Barretos Cancer Hospital. All clinical data were obtained from medical records and managed through RedCap Research Electronic Data Capture (Vanderbilt University, Nashville, TN, USA) (https://hcbredcap.com.br). The study protocol was approved by the Internal Review Board (IRB No. #8042014). Tumor samples were dissected from surgical specimens, snap-frozen in liquid nitrogen and stored at −80 ∘ C until DNA preparation. Histopathological evaluation was performed independently by two pathologists (CSN and LN) according to the last WHO classification system [50]. Nuclear grade was scored according to the Fuhrman classification system [51]. Patients without histologically confirmation or non-classic ccRCC histopathology were excluded.

4.2. DNA Isolation

Tumor DNA was isolated from frozen tissue with >60% tumor cellularity and necrosis was <20%. DNA isolation was performed using QIAsymphony DNA mini Kit following the manufacturer’s instructions. DNA quality and integrity were assessed by NanoDrop and quantified by Qubit® Fluorometric Quantitation (Life Technologies). Samples with insufficient material for research or of DNA poor quality were removed from the analysis.

4.3. DNA Labeling and Hybridization

CGH-array was performed according to the standard protocol from Agilent SurePrint G3 CGH+SNP 4 × 180 K (Agilent Technologies, Santa Clara, CA, USA) using a diploid DNA reference as control (Agilent Technologies). Normal DNA (reference) was labeled with Cyanine-3 and tumor DNA (test) with Cyanine-5 (Agilent Technologies). Labeled DNAs were purified, and the quality of labeling was assessed by NanoDrop (Thermo Fisher Scientific, Waltham, MA, USA). Equal quantities of DNA (test and reference) from each patient were hybridized into the slides for 24 h and cleaned with washed and stabilization buffers, plus acetonitrile. The slides were scanned using a SureScan Microarray Scanner (Agilent Technologies, Santa Clara, CA, USA) following the manufacturer’s recommendations and as previously reported by our group [52]. The recommended cut-off Derivative of the Log Ratio Spread (DLRS) considered for the aCGH experiments was 0.30.

4.4. Bioinformatics and Statistical Analysis

CNA analysis was performed using Nexus Copy Number software (BioDiscovery Inc., El Segundo, CA, USA) (http://www.biodiscovery.com/nexus-copy-number/). In this study, we used the definition that CNA (Copy Number Alteration) refers to somatic events. BioDiscovery’s SNP-FASST2 Segmentation Algorithm, is an extension of the FASST2 Segmentation Algorithm (a Hidden Markov Model [HMM]-based approach that unlike other common HMM methods does not aim to estimate the copy number state at each probe but uses many states to cover more possibilities, such as mosaic events, and then make calls based on a second level threshold). With the SNP-FASST2 algorithm, B-allele frequency probes are assigned to a range of possible states and a combination of the BAF and Log-R states are used to make the final copy number and allelic event calls. The recommended parameters for Agilent CGH+SNP arrays were: selection of ProcessedSignal from Agilent Feature Extraction (FE) software (Agilent Technologies), exclusion of saturated (IsStaturated) and non-uniform (IsFeatNonUnifOl) spots as flagged from FE, the significance threshold for segmentation was set at 5.0 × 10 − 6 also requiring a minimum of 3 probes per segment. The log ratio thresholds for single copy gain and single copy loss were set at 0.18 and −0.18, respectively. The log ratio thresholds for two or more copy gain and homozygous loss were set at 0.7 and −1.1 respectively. The Homozygous Frequency Threshold was set to 0.85. The Homozygous Value Threshold was set to 0.75. The Heterozygous Imbalance Threshold was set to 0.3. The minimum length was set at 1000 base pairs (bp). A 3:1 sex chromosome gain threshold was set to 1.2 and a 4:1 sex chromosome gain threshold was set to 1.7. Male sex chromosome big loss threshold was set to −1.1. The Significance Testing for Aberrant Copy number (STAC) [53], which is based on global frequency statistic approach, was applied to identify regions of the genome with a statistically high frequency of aberration. It was considered a p-value ≤ 0.05 with 1000 permutations as input parameter, as well as the aggregate cut-off filter of ≥ 20%. An additional filter was applied to regions showing 100% overlap with previously identified CNVs from the Database of Genomic Variants (DGV), inside the Nexus Copy Number platform. The exception were regions with somatic cancer genes identified by Cancer Gene Census. Significant regions were queried for oncogenic alterations through Cancer Genome Interpreter platform (CGI) (http://www.cancergenomeinterpreter.org). The threshold value used to stratify the CNA burden was within the 25th percentile, accounting for approximately 50 alterations, similarly to previous studies [54,55]. We also used the percentage of genome alteration (PGA) in the analysis. The Oncomine™ transcriptomic cancer profiling platform (Compendia Bioscience, Ann Arbor, MI, USA) (https://www.oncomine.org/) was used to explore changes in gene expression of the regions found. The premium version of Oncomine™ has currently 715 independent datasets and 86,733 samples normalized by the same parameters, including the TCGA data. Genes were investigated for gains and losses of expression in TCGA dataset (totaling 579 ccRCC clinical samples with DNA copy number data) [10,56]. The R Package for statistical analysis (https://www.R-project.org) was locally implemented in Galaxy platform (https://usegalaxy.org) for the associations. Correlation between frequencies of molecular alterations and clinical, histological and pathological features were performed using chi-square or Fisher’s exact test, both one-sided and two-sided tests, with a significance set at p≤ 0.05. The p-values were adjusted according to the Benjamini–Hochberg false-discovery (FDR) method using the p.adjust function in R (https://www.R-project.org). For this analysis, it was used the 12 regions identified by STAC and the following clinical features (Fuhrman grade, metastasis, tumor size, tumors T3–T4, clinical stage and death). The multivariate analysis by the Binary Logistic Regression method was performed using the Backward model with variables with significance equal to or below 0.2 after FDR correction in the univariate analysis. The IBM SPSS Statistics software version 23 (https://www.ibm.com/analytics/spss-statistics-software) was used to perform the multivariate and survival analysis. Overall survival (OS) was defined as the time from nephrectomy until the last contact or death. OS was assessed using the Kaplan–Meier method and log-rank comparisons, considering p≤ 0.05 as the significance value.

5. Conclusions

In summary, we report the first CNA landscape of Brazilian ccRCC. We found the most frequent aberrant regions were loss of 3p (87.3%), 14q (35.8%), and gains of 5q (59.7%) in accordance with previous studies. We identified that CNAs in 5q34–q35.3 (FLT4 and NSD1), 6q23.2–q23.3 (MYB), 9p21.3 (MLLT3) and 14q22.1 (NIN) associated with patient’s aggressive disease and identified putative CNA-associated genes. This study can pave the way to the understand ccRCC biology and the CNA biomarkers in Brazilian patients.

Supplementary Materials

The following are available at https://www.mdpi.com/1422-0067/22/5/2265/s1. Table S1: Summary of the main risk factors associated with ccRCC patients; Table S2: All 158 regions of interest identified by STAC analysis; Figure S1: Overall survival of the 92 ccRCC patients; Figure S2: Kaplan Meier plots according to stage (A), histological grade (B), metastasis (C) and tumor size (D).

Author Contributions

Conceptualization: R.M.R., A.F.E. and E.J.; methodology: H.C.S.S., R.A.d.S., L.N. and C.S.-N.; data curation: E.J., J.N.A.J., L.E.Z., F.M.C. and S.V.S.; writing—original draft preparation: F.G.F.; writing—review and final editing: F.G.F., H.C.S.S., J.N.A.J., R.A.d.S., L.E.Z., F.M.C., A.O.N.S., L.N., C.S.-N., S.V.S., E.J., R.M.R. and A.F.E.; visualization: F.G.F. and A.F.E.; supervision: R.M.R. and A.F.E.; project administration: R.M.R. and A.F.E.; funding acquisition: R.M.R. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by Sister Institution Network Fund (SINF) of Barretos Cancer Hospital and Public Ministry of Labor Campinas (Research, Prevention and Education of Occupational Cancer, Brazil). Fellowships from Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) and Barretos Cancer Hospital.

Institutional Review Board Statement

The study protocol was approved by the Internal Review Board (IRB No. #8042014).

Informed Consent Statement

Due the retrospective nature of the study, the Local Ethical Review Committee waived the need for patients written informed consent.

Data Availability Statement

Data is available from the corresponding authors on reasonable request.

Acknowledgments

We thank the Barretos Cancer Hospital biobank for processing the samples and the OncoGU (Grupo de Pesquisa em Oncologia Genitourinária) for the clinical oncology advisory.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Linehan, W.M.; Srinivasan, R.; Schmidt, L.S. The genetic basis of kidney cancer: A metabolic disease. Nat. Rev. Urol. 2010, 7, 277–285. [Google Scholar] [CrossRef] [PubMed]
  2. Linehan, W.M. Genetic basis of kidney cancer: Role of genomics for the development of disease-based therapeutics. Genome Res. 2012, 22, 2089–2100. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Jonasch, E.; Gao, J.; Rathmell, W.K. State of the Art Review: Renal cell carcinoma. BMJ 2014. [Google Scholar] [CrossRef] [PubMed]
  4. Ricketts, C.J.; Cubas, A.A.D.; Fan, H.; Smith, C.C.; Lang, M.; Reznik, E.; Bowlby, R.; Gibb, E.A.; Akbani, R.; Beroukhim, R.; et al. The Cancer Genome Atlas Comprehensive Molecular Characterization of Renal Cell Carcinoma. Cell Rep. 2018, 23, 3698. [Google Scholar] [CrossRef] [PubMed]
  5. Bray, F.; Ferlay, J.; Soerjomataram, I.; Siegel, R.L.; Torre, L.A.; Jemal, A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries: Global Cancer Statistics 2018. CA Cancer J. Clin. 2018, 68, 394–424. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Middleton, R.G. Surgery for metastatic renal cell carcinoma. J. Urol. 1967, 97, 973–977. [Google Scholar] [CrossRef]
  7. Srigley, J.R.; Delahunt, B.; Eble, J.N.; Egevad, L.; Epstein, J.I.; Grignon, D.; Hes, O.; Moch, H.; Montironi, R.; Tickoo, S.K.; et al. The International Society of Urological Pathology (ISUP) Vancouver Classification of Renal Neoplasia. Am. J. Surg. Pathol. 2013, 37, 1469–1489. [Google Scholar] [CrossRef]
  8. Quddus, M.B.; Pratt, N.; Nabi, G. Chromosomal aberrations in renal cell carcinoma: An overview with implications for clinical practice. Urol. Ann. 2019, 11, 6–14. [Google Scholar] [CrossRef]
  9. Higgins, J.P. Gene Array Studies in Renal Neoplasia. Sci. World J. 2006, 6, 502–511. [Google Scholar] [CrossRef]
  10. The Cancer Genome Atlas Research Network. Comprehensive molecular characterization of clear cell renal cell carcinoma. Nature 2013, 499, 43–49. [Google Scholar] [CrossRef] [Green Version]
  11. Shen, C.; Beroukhim, R.; Schumacher, S.E.; Zhou, J.; Chang, M.; Signoretti, S.; Kaelin, W.G. Genetic and Functional Studies Implicate HIF1A as a 14q Kidney Cancer Suppressor Gene. Cancer Discov. 2011, 1, 222–235. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Klatte, T.; Kroeger, N.; Rampersaud, E.N.; Birkhäuser, F.D.; Logan, J.E.; Sonn, G.; Riss, J.; Rao, P.N.; Kabbinavar, F.F.; Belldegrun, A.S.; et al. Gain of chromosome 8q is associated with metastases and poor survival of patients with clear cell renal cell carcinoma. Cancer 2012, 118, 5777–5782. [Google Scholar] [CrossRef] [PubMed]
  13. Huang, Y.S.; Liu, W.B.; Han, F.; Yang, J.T.; Hao, X.L.; Chen, H.Q.; Jiang, X.; Yin, L.; Ao, L.; Cui, Z.H.; et al. Copy number variations and expression of MPDZ are prognostic biomarkers for clear cell renal cell carcinoma. Oncotarget 2017, 8, 78713–78725. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Hakimi, A.A.; Reznik, E.; Lee, C.H.; Creighton, C.J.; Brannon, A.R.; Luna, A.; Aksoy, B.A.; Liu, E.M.; Shen, R.; Lee, W.; et al. An Integrated Metabolic Atlas of Clear Cell Renal Cell Carcinoma. Cancer Cell 2016, 29, 104–116. [Google Scholar] [CrossRef] [Green Version]
  15. Zhou, W.; Yang, F.; Xu, Z.; Luo, M.; Wang, P.; Guo, Y.; Nie, H.; Yao, L.; Jiang, Q. Comprehensive Analysis of Copy Number Variations in Kidney Cancer by Single-Cell Exome Sequencing. Front. Genet. 2019, 10, 1379. [Google Scholar] [CrossRef] [Green Version]
  16. Klatte, T.; Rao, P.N.; de Martino, M.; LaRochelle, J.; Shuch, B.; Zomorodian, N.; Said, J.; Kabbinavar, F.F.; Belldegrun, A.S.; Pantuck, A.J. Cytogenetic profile predicts prognosis of patients with clear cell renal cell carcinoma. J. Clin. Oncol. Off. J. Am. Soc. Clin. Oncol. 2009, 27, 746–753. [Google Scholar] [CrossRef]
  17. Gunawan, B.; Huber, W.; Holtrup, M.; Heydebreck, A.V.; Efferth, T.; Poustka, A.; Ringert, R.H.; Jakse, G.; Füzesi, L. Prognostic Impacts of Cytogenetic Findings in Clear Cell Renal Cell Carcinoma: Gain of 5q31–qter Predicts a Distinct Clinical Phenotype with Favorable Prognosis. Cancer Res. 2001, 61, 7731–7738. [Google Scholar]
  18. Kovacs, G.; Frisch, S. Clonal chromosome abnormalities in tumor cells from patients with sporadic renal cell carcinomas. Cancer Res. 1989, 49, 651–659. [Google Scholar]
  19. Wei, E.Y.; Hsieh, J.J. A river model to map convergent cancer evolution and guide therapy in RCC. Nat. Rev. Urol. 2015, 12, 706–712. [Google Scholar] [CrossRef] [Green Version]
  20. Moch, H.; Presti, J.C.; Sauter, G.; Buchholz, N.; Jordan, P.; Mihatsch, M.J.; Waldman, F.M. Genetic aberrations detected by comparative genomic hybridization are associated with clinical outcome in renal cell carcinoma. Cancer Res. 1996, 56, 27–30. [Google Scholar]
  21. Brunelli, M.; Eccher, A.; Gobbo, S.; Ficarra, V.; Novara, G.; Cossu-Rocca, P.; Bonetti, F.; Menestrina, F.; Cheng, L.; Eble, J.N.; et al. Loss of chromosome 9p is an independent prognostic factor in patients with clear cell renal cell carcinoma. Mod. Pathol. 2008, 21, 1–6. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Elfving, P.; Mandahl, N.; Lundgren, R.; Limon, J.; Bak-Jensen, E.; Fernö, M.; Olsson, H.; Mitelman, F. Prognostic implications of cytogenetic findings in kidney cancer. Br. J. Urol. 1997, 80, 698–706. [Google Scholar] [CrossRef] [PubMed]
  23. De Oliveira, D.; Dall’Oglio, M.F.; Reis, S.T.; Zerati, M.; Souza, I.C.; Leite, K.R.; Srougi, M. Chromosome 9p deletions are an independent predictor of tumor progression following nephrectomy in patients with localized clear cell renal cell carcinoma. Urol. Oncol. 2014, 32. [Google Scholar] [CrossRef]
  24. Watkins, T.B.K.; Lim, E.L.; Petkovic, M.; Elizalde, S.; Birkbak, N.J.; Wilson, G.A.; Moore, D.A.; Grönroos, E.; Rowan, A.; Dewhurst, S.M.; et al. Pervasive chromosomal instability and karyotype order in tumour evolution. Nature 2020, 587, 126–132. [Google Scholar] [CrossRef] [PubMed]
  25. Huang, Y.; Wang, J.; Jia, P.; Li, X.; Pei, G.; Wang, C.; Fang, X.; Zhao, Z.; Cai, Z.; Yi, X.; et al. Clonal architectures predict clinical outcome in clear cell renal cell carcinoma. Nat. Commun. 2019, 10. [Google Scholar] [CrossRef] [PubMed]
  26. Mitchell, T.J.; Turajlic, S.; Rowan, A.; Nicol, D.; Farmery, J.H.R.; O’Brien, T.; Martincorena, I.; Tarpey, P.; Angelopoulos, N.; Yates, L.R.; et al. Timing the Landmark Events in the Evolution of Clear Cell Renal Cell Cancer: TRACERx Renal. Cell 2018, 173, 611–623. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Carrot-Zhang, J.; Chambwe, N.; Damrauer, J.S.; Knijnenburg, T.A.; Robertson, A.G.; Yau, C.; Zhou, W.; Berger, A.C.; Huang, K.L.; Newberg, J.Y.; et al. Comprehensive Analysis of Genetic Ancestry and Its Molecular Correlates in Cancer. Cancer Cell 2020, 37, 639–654. [Google Scholar] [CrossRef] [PubMed]
  28. Yuan, J.; Hu, Z.; Mahal, B.A.; Zhao, S.D.; Kensler, K.H.; Pi, J.; Hu, X.; Zhang, Y.; Wang, Y.; Jiang, J.; et al. Integrated Analysis of Genetic Ancestry and Genomic Alterations across Cancers. Cancer Cell 2018, 34, 549–560. [Google Scholar] [CrossRef] [Green Version]
  29. Kardas, I.; Mrózek, K.; Babinska, M.; Krajka, K.; Hadaczek, P.; Lubinski, J.; Roszkiewicz, A.; Kuziemska, E.; Limon, J. Cytogenetic and molecular findings in 75 clear cell renal cell carcinomas. Oncol. Rep. 2005, 13, 949–956. [Google Scholar] [CrossRef]
  30. Monzon, F.A.; Alvarez, K.; Peterson, L.; Truong, L.; Amato, R.J.; Hernandez-McClain, J.; Tannir, N.; Parwani, A.V.; Jonasch, E. Chromosome 14q loss defines a molecular subtype of clear-cell renal cell carcinoma associated with poor prognosis. Mod. Pathol. 2011, 24, 1470–1479. [Google Scholar] [CrossRef] [Green Version]
  31. La Rochelle, J.; Klatte, T.; Dastane, A.; Rao, N.; Seligson, D.; Said, J.; Shuch, B.; Zomorodian, N.; Kabbinavar, F.; Belldegrun, A.; et al. Chromosome 9p deletions identify an aggressive phenotype of clear cell renal cell carcinoma. Cancer 2010, 116, 4696–4702. [Google Scholar] [CrossRef]
  32. Mogensen, M.M.; Malik, A.; Piel, M.; Bouckson-Castaing, V.; Bornens, M. Microtubule minus-end anchorage at centrosomal and non-centrosomal sites: The role of ninein. J. Cell Sci. 2000, 113, 3013–3023. [Google Scholar]
  33. Tamborero, D.; Rubio-Perez, C.; Deu-Pons, J.; Schroeder, M.P.; Vivancos, A.; Rovira, A.; Tusquets, I.; Albanell, J.; Rodon, J.; Tabernero, J.; et al. Cancer Genome Interpreter annotates the biological and clinical relevance of tumor alterations. Genome Med. 2018, 10, 25. [Google Scholar] [CrossRef]
  34. Fan, Q.W.; Cheng, C.K.; Nicolaides, T.P.; Hackett, C.S.; Knight, Z.A.; Shokat, K.M.; Weiss, W.A. A dual phosphoinositide-3-kinase alpha/mTOR inhibitor cooperates with blockade of epidermal growth factor receptor in PTEN-mutant glioma. Cancer Res. 2007, 67, 7960–7965. [Google Scholar] [CrossRef] [Green Version]
  35. Iorio, F.; Knijnenburg, T.A.; Vis, D.J.; Bignell, G.R.; Menden, M.P.; Schubert, M.; Aben, N.; Gonçalves, E.; Barthorpe, S.; Lightfoot, H.; et al. A Landscape of Pharmacogenomic Interactions in Cancer. Cell 2016, 166, 740–754. [Google Scholar] [CrossRef] [Green Version]
  36. Pfister, S.; Markkanen, E.; Jiang, Y.; Sarkar, S.; D’Angiolella, V.; Dianov, G.; Ryan, A. A novel synthetic lethal interaction between the histone mark H3K36me3 and checkpoint kinases. Eur. J. Cancer 2014, 50, 71. [Google Scholar] [CrossRef]
  37. Peña-Llopis, S.; Vega-Rubín-de Celis, S.; Liao, A.; Leng, N.; Pavía-Jiménez, A.; Wang, S.; Yamasaki, T.; Zhrebker, L.; Sivanand, S.; Spence, P.; et al. BAP1 loss defines a new class of renal cell carcinoma. Nat. Genet. 2012, 44, 751–759. [Google Scholar] [CrossRef]
  38. Yap, T.A.; Walton, M.I.; Hunter, L.J.K.; Valenti, M.; de Haven Brandon, A.; Eve, P.D.; Ruddle, R.; Heaton, S.P.; Henley, A.; Pickard, L.; et al. Preclinical pharmacology, antitumor activity and development of pharmacodynamic markers for the novel, potent AKT inhibitor CCT128930. Mol. Cancer Ther. 2011, 10, 360–371. [Google Scholar] [CrossRef] [Green Version]
  39. Feng, C.; Sun, Y.; Ding, G.; Wu, Z.; Jiang, H.; Wang, L.; Ding, Q.; Wen, H. PI3KB inhibitor TGX221 selectively inhibits renal cell carcinoma cells with both VHL and SETD2 mutations and links multiple pathways. Sci. Rep. 2015, 5, 9465. [Google Scholar] [CrossRef] [Green Version]
  40. Eghtedar, A.; Verstovsek, S.; Estrov, Z.; Burger, J.; Cortes, J.; Bivins, C.; Faderl, S.; Ferrajoli, A.; Borthakur, G.; George, S.; et al. Phase 2 study of the JAK kinase inhibitor ruxolitinib in patients with refractory leukemias, including postmyeloproliferative neoplasm acute myeloid leukemia. Blood 2012, 119, 4614–4618. [Google Scholar] [CrossRef]
  41. Chase, A.; Bryant, C.; Score, J.; Haferlach, C.; Grossmann, V.; Schwaab, J.; Hofmann, W.K.; Reiter, A.; Cross, N.C.P. Ruxolitinib as potential targeted therapy for patients with JAK2 rearrangements. Haematologica 2013, 98, 404–408. [Google Scholar] [CrossRef]
  42. Motzer, R.J.; Rini, B.I.; McDermott, D.F.; Redman, B.G.; Kuzel, T.M.; Harrison, M.R.; Vaishampayan, U.N.; Drabkin, H.A.; George, S.; Logan, T.F.; et al. Nivolumab for Metastatic Renal Cell Carcinoma: Results of a Randomized Phase II Trial. J. Clin. Oncol. 2015, 33, 1430–1437. [Google Scholar] [CrossRef]
  43. Li, S.Q.; Cheuk, A.T.; Shern, J.F.; Song, Y.K.; Hurd, L.; Liao, H.; Wei, J.S.; Khan, J. Targeting wild-type and mutationally activated FGFR4 in rhabdomyosarcoma with the inhibitor ponatinib (AP24534). PLoS ONE 2013, 8, e76551. [Google Scholar] [CrossRef] [Green Version]
  44. Patel, J.P.; Gönen, M.; Figueroa, M.E.; Fernandez, H.; Sun, Z.; Racevskis, J.; Van Vlierberghe, P.; Dolgalev, I.; Thomas, S.; Aminova, O.; et al. Prognostic relevance of integrated genetic profiling in acute myeloid leukemia. N. Engl. J. Med. 2012, 366, 1079–1089. [Google Scholar] [CrossRef] [Green Version]
  45. Treon, S.P.; Tripsas, C.K.; Meid, K.; Warren, D.; Varma, G.; Green, R.; Argyropoulos, K.V.; Yang, G.; Cao, Y.; Xu, L.; et al. Ibrutinib in previously treated Waldenström’s macroglobulinemia. N. Engl. J. Med. 2015, 372, 1430–1440. [Google Scholar] [CrossRef] [Green Version]
  46. Zhao, S.; Adebiyi, M.G.; Zhang, Y.; Couturier, J.P.; Fan, X.; Zhang, H.; Kellems, R.E.; Lewis, D.E.; Xia, Y. Sphingosine-1-phosphate receptor 1 mediates elevated IL-6 signaling to promote chronic inflammation and multitissue damage in sickle cell disease. FASEB J. 2018, 32, 2855–2865. [Google Scholar] [CrossRef] [Green Version]
  47. Bhattacharya, A.; Bense, R.D.; Urzúa-Traslaviña, C.G.; de Vries, E.G.E.; van Vugt, M.A.T.M.; Fehrmann, R.S.N. Transcriptional effects of copy number alterations in a large set of human cancers. Nat. Commun. 2020, 11, 715. [Google Scholar] [CrossRef]
  48. Arai, E.; Ushijima, S.; Tsuda, H.; Fujimoto, H.; Hosoda, F.; Shibata, T.; Kondo, T.; Imoto, I.; Inazawa, J.; Hirohashi, S.; et al. Genetic clustering of clear cell renal cell carcinoma based on array-comparative genomic hybridization: Its association with DNA methylation alteration and patient outcome. Clin. Cancer Res. 2008, 14, 5531–5539. [Google Scholar] [CrossRef] [Green Version]
  49. Burrell, R.A.; McClelland, S.E.; Endesfelder, D.; Groth, P.; Weller, M.C.; Shaikh, N.; Domingo, E.; Kanu, N.; Dewhurst, S.M.; Gronroos, E.; et al. Replication stress links structural and numerical cancer chromosomal instability. Nature 2013, 494, 492–496. [Google Scholar] [CrossRef] [Green Version]
  50. Eble, J.N.; Weltgesundheitsorganisation. Pathology and Genetics of Tumours of the Urinary System and Male Genital Organs, Proceedings of the Editorial and Consensus Conference in Lyon, France, 14–18 December 2002, reprint ed.; Cancer, I.A., Ed.; Number 6 in World Health Organization Classification of Tumours; IARC Press: Lyon, France, 2006. [Google Scholar]
  51. Fuhrman, S.A.; Lasky, L.C.; Limas, C. Prognostic significance of morphologic parameters in renal cell carcinoma. Am. J. Surg. Pathol. 1982, 6, 655–663. [Google Scholar] [CrossRef]
  52. Bidinotto, L.T.; Torrieri, R.; Mackay, A.; Almeida, G.C.; Viana-Pereira, M.; Cruvinel-Carloni, A.; Spina, M.L.; Campanella, N.C.; Menezes, W.P.D.; Clara, C.A.; et al. Copy Number Profiling of Brazilian Astrocytomas. G3 Genes Genomes Genet. 2016, 6, 1867–1878. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Diskin, S.J.; Eck, T.; Greshock, J.; Mosse, Y.P.; Naylor, T.; Stoeckert, C.J.; Weber, B.L.; Maris, J.M.; Grant, G.R. STAC: A method for testing the significance of DNA copy number aberrations across multiple array-CGH experiments. Genome Res. 2006, 16, 1149–1158. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Hieronymus, H.; Murali, R.; Tin, A.; Yadav, K.; Abida, W.; Moller, H.; Berney, D.; Scher, H.; Carver, B.; Scardino, P.; et al. Tumor copy number alteration burden is a pan-cancer prognostic factor associated with recurrence and death. eLife 2018, 7, e37294. [Google Scholar] [CrossRef] [PubMed]
  55. Lu, Z.; Chen, H.; Li, S.; Gong, J.; Li, J.; Zou, J.; Wu, L.; Yu, J.; Han, W.; Sun, H.; et al. Tumor copy-number alterations predict response to immune-checkpoint-blockade in gastrointestinal cancer. J. Immunother. Cancer 2020, 8, e000374. [Google Scholar] [CrossRef] [PubMed]
  56. Beroukhim, R.; Brunet, J.P.; Di Napoli, A.; Mertz, K.D.; Seeley, A.; Pires, M.M.; Linhart, D.; Worrell, R.A.; Moch, H.; Rubin, M.A.; et al. Patterns of gene expression and copy-number alterations in VHL disease-associated and sporadic clear cell carcinoma of the kidney. Cancer Res. 2009, 69, 4674–4681. [Google Scholar] [CrossRef] [Green Version]
Ijms 22 02265 g001 550
Figure 1. Somatic DNA copy number aberrations (CNAs) of 92 patients of ccRCC. In (A) The linear distribution of DNA gains (red) and losses (blue) and in (B) Chromosomal representation of all the aberrations.
Figure 1. Somatic DNA copy number aberrations (CNAs) of 92 patients of ccRCC. In (A) The linear distribution of DNA gains (red) and losses (blue) and in (B) Chromosomal representation of all the aberrations.
Ijms 22 02265 g001
Ijms 22 02265 g002 550
Figure 2. Circos Plot showing the correlation of CNAs and gene expression from Oncomine data. The periphery circle represents regions of Copy Number Aberrations along the genome per chromosome. Red represents copy number gains and blue copy number losses. The inner circle contains genes with concomitant alteration in gene expression. Green dots represents down-regulated genes and red up-regulated genes.
Figure 2. Circos Plot showing the correlation of CNAs and gene expression from Oncomine data. The periphery circle represents regions of Copy Number Aberrations along the genome per chromosome. Red represents copy number gains and blue copy number losses. The inner circle contains genes with concomitant alteration in gene expression. Green dots represents down-regulated genes and red up-regulated genes.
Ijms 22 02265 g002
Ijms 22 02265 g003 550
Figure 3. Overall survival of chromosome 14 q 22.1 which encompasses the NIN gene.
Figure 3. Overall survival of chromosome 14 q 22.1 which encompasses the NIN gene.
Ijms 22 02265 g003
Table
Table 1. Patient clinicopathological features.
Table 1. Patient clinicopathological features.
Clinical DataCategoryPatients (Frequency)
GenderMale50 (54.3%)
Female42 (45.7%)
Alive69 (75.0%)
Status follow-upDeceased23 (25.0%)
AgeAverage (min-max)58 (20–81)
I/II60 (65.2%)
Clinical stageIII/IV32 (44.8%)
TT1a-b47 (51.1%)
T2a-b23 (25.0%)
T3a-c17 (18.5%)
T45 (5.4%)
N077 (83.7%)
N19 (9.8%)
NNx6 (6.5%)
MM075 (81.5%)
M117 (18.5%)
<7 cm58 (63.0%)
Tumor size≥7 cm34 (37.0%)
Furhman Grade120 (21.7%)
249 (53.3%)
316 (17.4%)
47 (7.6%)
White64 (69.6%)
Black3 (3.3%)
Yellow1 (1.0%)
Brown *20 (21.7%)
Self-Declared EthnicityUnknown4 (4.4%)
Systemic TherapySunitinib12 (13.0%)
Interferon4 (4.3%)
Pazopanib1 (1.1%)
None74 (80.4%)
Other1 (1.1%)
Alive without disease64 (69.6%)
Alive with disease progression2 (2.1%)
Progression and death from the disease18 (19.6%)
Disease status (last follow-up)Died without disease recurrence8 (8.7%)
* Brown = Brazilian mixed; T = size or direct extent of the primary tumor; N = degree of spread to regional. lymph nodes; M = presence of distant metastasis.
Table
Table 2. Frequency of significant regions according to STAC with frequency ≥20% in the present casuistic and in TCGA.
Table 2. Frequency of significant regions according to STAC with frequency ≥20% in the present casuistic and in TCGA.
ChromosomeCytobandEventPatients
Frequency
(%)		TCGA
Frequency
(%)		Cancer
Genes
Chr3p21.2–p21.3DEL83.772.2BAP1, SETD2
Chr3p22.2DEL83.772.2MYD88
Chr3p25.2DEL83.772.2PPARG
Chr3q11.2-q13.11DEL20.625.7TFG, CBLB
Chr5q34-q35.3AMP59.747.6RANBP17, NPM1, FGFR4, NSD1, FLT4
Chr6q23.2–q23.3DEL29.319.0MYB
Chr7p21.1AMP29.322.9ETV1
Chr9p21.3DEL28.217.9MLLT3
Chr9p24.1DEL27.120.7JAK2, CD274
Chr10q23.31DEL25.013.6PTEN
Chr14q22.1DEL35.929.0NIN
Chr16q21AMP20.611.9CDH11
AMP = Amplification; DEL = Deletion.
Table
Table 3. Association of clinical features with significant regions identified by STAC.
Table 3. Association of clinical features with significant regions identified by STAC.
ChromosomeCytobandEventClinical
Features		p ValueFDR
Chr3p21.2DELFuhrman Grade (1-2)0.020.15
Chr3p25.2DELFuhrman Grade (1–2)0.010.14
Chr5q34–q35.1AMPMetastasis0.030.18
Tumor size (≥7 cm)0.010.14
Fuhrman Grade (3–4)0.010.14
Chr7p21.2AMPTumors T3/T40.020.15
Chr9p21.3DELMetastasis0.010.14
Tumor size (≥7 cm)0.020.15
Tumors T3–T40.040.19
Clinical stage (3–4)0.010.14
Tumor size (≥7 cm)0.040.19
Chr9p24.1DELTumors T3-T40.040.19
Chr14q22.1DELClinical stage (3–4)0.040.19
Outcome (death)0.020.15
Metastasis0.030.19
Fuhrman Grade (3–4)0.010.14
AMP = Amplification; DEL = Deletion.
Table
Table 4. Multivariate analysis of clinical features associated with STAC significant regions.
Table 4. Multivariate analysis of clinical features associated with STAC significant regions.
Clinical FeatureRegionRR95% CIp Value
Tumors T3/T47p21.2 (AMP)3.01.11–8.040.02
Tumor size (>7 cm)5q35.2–q35.3 (AMP)2.91.08–8.140.03
9p24.1 (DEL)3.31.25–9.040.01
Clinical Stage (3–4)9p21.3 (DEL)3.91.51–10.200.005
Outcome (Death)14q22.1 (DEL)2.91.12–7.840.02
Fuhrman Grade (3–4)7p21.2 (AMP)3.41.28–9.100.01
Metastasis9p21.3 (DEL)5.21.72–16.020.003
AMP = Amplification; DEL = Deletion; RR = Relative Risk; CI = Confidence interval.
Table
Table 5. Multivariate analysis of CNA burden.
Table 5. Multivariate analysis of CNA burden.
Clinical FeatureCategoryRR95% CIp Value
Tumor size (≥7 cm)> 50 CNA7.82.14–28.550.002
Outcome (death)> 50 CNA3.71.02–14.020.04
Fuhrman Grade (3–4)> 50 CNA14.11.80–111.130.01
RR = Relative Risk; CI = Confidence interval.
Table
Table 6. Enrichment analysis for genes associated with gene expression through Oncomine platform.
Table 6. Enrichment analysis for genes associated with gene expression through Oncomine platform.
ChrCytobandGeneCN
Event		ExpressionClinical
Features		p-Value
3p21.2–p21.3BAP1LossDownregulationFuhrman Grade (1–2)0.02
5q34–q35.3FLT4GainUpregulationTumor size (>7 cm)0.007
5q34–q35.3NSD1GainUpregulationTumor size (>7 cm)0.007
Clinical Stage (3-4)0.01
6q23.2–q23.3MYBLossDownregulationFuhrman Grade (3-4)0.02
9p21.3MLLT3LossDownregulationMetastasis0.005
Clinical Stage (3–4)0.02
Tumor size (>7 cm)0.01
Outcome (death)0.02
Fuhrman Grade (3–4)0.01
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Fernandes, F.G.; Silveira, H.C.S.; Júnior, J.N.A.; da Silveira, R.A.; Zucca, L.E.; Cárcano, F.M.; Sanches, A.O.N.; Neder, L.; Scapulatempo-Neto, C.; Serrano, S.V.; et al. Somatic Copy Number Alterations and Associated Genes in Clear-Cell Renal-Cell Carcinoma in Brazilian Patients. Int. J. Mol. Sci. 2021, 22, 2265. https://doi.org/10.3390/ijms22052265

AMA Style

Fernandes FG, Silveira HCS, Júnior JNA, da Silveira RA, Zucca LE, Cárcano FM, Sanches AON, Neder L, Scapulatempo-Neto C, Serrano SV, et al. Somatic Copy Number Alterations and Associated Genes in Clear-Cell Renal-Cell Carcinoma in Brazilian Patients. International Journal of Molecular Sciences. 2021; 22(5):2265. https://doi.org/10.3390/ijms22052265

Chicago/Turabian Style

Fernandes, Flávia Gonçalves, Henrique Cesar Santejo Silveira, João Neif Antonio Júnior, Rosana Antunes da Silveira, Luis Eduardo Zucca, Flavio Mavignier Cárcano, André Octavio Nicolau Sanches, Luciano Neder, Cristovam Scapulatempo-Neto, Sergio Vicente Serrano, and et al. 2021. "Somatic Copy Number Alterations and Associated Genes in Clear-Cell Renal-Cell Carcinoma in Brazilian Patients" International Journal of Molecular Sciences 22, no. 5: 2265. https://doi.org/10.3390/ijms22052265

APA Style

Fernandes, F. G., Silveira, H. C. S., Júnior, J. N. A., da Silveira, R. A., Zucca, L. E., Cárcano, F. M., Sanches, A. O. N., Neder, L., Scapulatempo-Neto, C., Serrano, S. V., Jonasch, E., Reis, R. M., & Evangelista, A. F. (2021). Somatic Copy Number Alterations and Associated Genes in Clear-Cell Renal-Cell Carcinoma in Brazilian Patients. International Journal of Molecular Sciences, 22(5), 2265. https://doi.org/10.3390/ijms22052265

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop