Genome-wide analysis of differentially expressed genes and splicing isoforms in clear cell renal cell carcinoma

Identification of ATP6V0A4 as a potential biomarker in renal cell carcinoma using integrated bioinformatics analysis

Clear cell renal cell carcinoma (ccRCC) is the most common pathological type of renal cancer, and is associated with a high mortality rate, which is related to high rates of tumor recurrence and metastasis. The aim of the present study was to identify reliable molecular biomarkers with high specificity and sensitivity for ccRCC. A total of eight ccRCC-related expression profiles were downloaded from Gene Expression Omnibus for integrated bioinformatics analysis to screen for significantly differentially expressed genes (DEGs). Reverse transcription-quantitative (RT-q)PCR, western blotting and immunohistochemistry staining assays were performed to evaluate the expression levels of candidate biomarkers in ccRCC tissues and cell lines. In total, 255 ccRCC specimens and 165 adjacent normal kidney specimens were analyzed, and 344 significant DEGs, consisting of 115 upregulated DEGs and 229 downregulated DEGs, were identified. The results of Gene Ontology analysis suggested a significant enrichment of DEGs in 'organic anion transport' and 'small molecule catabolic process' in biological processes, in 'apical plasma membrane' and 'apical part of the cell' in cell components, and in 'anion transmembrane transporter activity' and 'active transmembrane transporter activity' in molecular functions. Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis indicated that the DEGs were significantly enriched in the 'phagosome', the 'PPAR signaling pathway', 'complement and coagulation cascades', the 'HIF-1 signaling pathway' and 'carbon metabolism'. Next, 7 hub genes (SUCNR1, CXCR4, VCAN, CASR, ATP6V0A4, VEGFA and SERPINE1) were identified and validated using The Cancer Genome Atlas database. Survival analysis showed that low expression of ATP6V0A4 was associated with a poor prognosis in patients with ccRCC. Additionally, received operating characteristic curves indicated that ATP6V0A4 could distinguish ccRCC samples from normal kidney samples. Furthermore, RT-qPCR, western blotting and immunohistochemistry staining results showed that ATP6V0A4 was significantly downregulated in ccRCC tissues and cell lines. In conclusion, ATP6V0A4 may be involved in tumor progression and regarded as a potential therapeutic target for the recurrence and metastasis of ccRCC.

PBX1, EMCN and ERG are associated with the sub-clusters and the prognosis of VHL mutant clear cell renal cell carcinoma

The molecular heterogeneity of primary clear cell renal cell carcinoma (ccRCC) has been reported. However, the classifications of Von Hippel-Lindau (VHL) mutant ccRCC are unclear. Here, VHL mutant ccRCC from The Cancer Genome Atlas and E-MTAB-1980 datasets were divided into two sub-clusters through non-negative matrix factorization algorithm. Most VHL mutant ccRCC patients in sub-cluster2 were with pathological T1 stage and VHL mutant ccRCC patients in sub-cluster1 were with decreased overall survival. DNA replication and homologous recombination scores were higher, while, WNT signaling pathway and regulation of autophagy scores were lower in sub-cluster1 VHL mutant ccRCC. Moreover, PBX1 transcriptional scores and mRNA expressions were lower in sub-cluster1 VHL mutant ccRCC patients and were associated with the overall survival of VHL mutant ccRCC. Furthermore, PBX1 associated genes EMCN and ERG were down-regulated in sub-cluster1 VHL mutant ccRCC and overall survival was decreased in EMCN or ERG lowly expressed VHL mutant ccRCC patients. Also, PBX1 and EMCN were down-regulated in ccRCC tissues, compared with normal kidney tissues. At last, we constructed risk models based on PBX1, EMCN and EGR expression features. With the increase of the risk score, the number of death of VHL mutant ccRCC patients was increased.

RHBDF2 gene functions are correlated to facilitated renal clear cell carcinoma progression

Background

The rhomboids are a family of multi-transmembrane proteins, many of which have been implicated in facilitating tumor progression. Little is yet known, however, about rhomboid-associated biomarkers in cancers. An analysis of such biomarkers could yield important insights into the role of the rhomboids in cancer pathology.

Methods

In this study, we carried out the univariate Cox regression analysis and compared gene expression patterns of several rhomboid genes in 30 types of cancers by using The Cancer Genome Atlas (TCGA) database and the methods delineated in Gene Expression Profiling Interactive Analysis (GEPIA). We then used datasets GSE47032, GSE126964, GSE68417 and 75 paired pathological specimens to verify the influences of the rhomboid genes in cancer progression. Moreover, we carried out Weighted Gene Correlation Network Analysis (WGCNA) to investigate gene-related functions and we exploited potential correlations between rhomboid genes expression and immune cell infiltration in cancer tissues. Furthermore, we constructed gene-knockdown cancer cell lines to investigate rhomboid gene functions.

Results

We find that kidney renal clear cell carcinoma (KIRC) disease progression is affected by fluctuations in the expression of a number of the rhomboid family of genes and, more specifically, high levels of RHBDF2 gene expression are a good indicator of poor prognosis of the disease, as patients with high RHBDF2 expression levels exhibit less favorable survival rates compared to those with low RHBDF2 levels. Silencing of the RHBDF2 gene in KIRC cell lines leads to significantly diminished cell proliferation and migration; this is in good agreement with the identification of an enhanced presence of a number of cell growth and migration promoting signaling molecules in KIRC tumors. We found that, although high level of RHBDF2 correlated with increased infiltration of lymphocytes in cancer tissues, artificially overexpressed RHBDF2 led to an inhibition of the activity of the infiltrated immune cells through sustaining PD-L1 protein level. Furthermore, we show that RHBDF2 related cell migration and PD-L1 regulation were potentially mediated by EGFR signaling pathway.

Conclusions

RHBDF2 gene functions are correlated to facilitated renal clear cell carcinoma progression and may serve as a critical prognostic biomarker for the disease.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12935-021-02277-0.

Integrated bioinformatic analysis identifies COL4A3, COL4A4, and KCNJ1 as key biomarkers in Wilms tumor

Wilms tumor (WT) is one of the most common pediatric solid tumors, affecting 1 in 10,000 children, worldwide. A subset of WT patients has poor prognosis, which is associated with a high risk of advanced and/or recurrent disease. Therefore, candidate markers are urgently needed for the diagnosis and effective treatment of WT. We evaluated three mRNA microarray datasets to identify the differences between normal kidney tissue and WT tissue. Gene expression profiling revealed 130 differentially expressed genes (DEGs). Enrichment analysis and gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were performed for the DEGs. Subsequently, we established a protein-protein interaction (PPI) network to reveal the associations among the DEGs and selected 10 hub genes, all of which were downregulated in WT. The expression of COL4A3, COL4A4, KCNJ1, MME, and SLC12A1 in WT tissues was significantly lower than that in normal renal tissues. Survival analyses using the Kaplan-Meier method showed that patients with WT and low expression of COL4A3, COL4A4, and KCNJ1 exhibited remarkably poor overall survival. The correlations among COL4A3, COL4A4, and KCNJ1 in WT were analyzed using cBioPortal; COL4A3, COL4A4, and KCNJ1 were positively correlated with each other. Thus, these genes were considered clinically significant and might therefore play important roles in carcinogenesis and the development of WT.

HnRNP A1 - mediated alternative splicing of CCDC50 contributes to cancer progression of clear cell renal cell carcinoma via ZNF395

BACKGROUND

Aberrant alternative splicing events play critical roles in carcinogenesis and progression of many cancers, while sparse studies regarding to alternative splicing are available for clear cell renal cell carcinoma (ccRCC). We identified that alternative splicing of coiled-coil domain containing 50 (CCDC50) was dysregulated in ccRCC, whereas the clinical significance of this splicing event and its splicing regulation mechanisms were still elusive.

METHODS

Bioinformatic algorithm was utilized to identify significant exon skipping events in ccRCC via exon sequencing data from The Cancer Genome Atlas. Semi-quantitative real-time polymerase chain reaction and western blot were used to validate the aberrant expression of different transcripts in renal cancer tissues, cell lines and corresponding noncancerous controls. Short hairpin RNA targeting CCDC50 and overexpressing plasmids for each transcript were introduced into ccRCC cell lines, followed by a series of in vitro and in vivo functional experiments. Moreover, a panel of splicing factors were identified and their roles on splicing regulation of CCDC50 precursor mRNA (pre-mRNA) were studied. Furthermore, RNAseq data were analyzed to elucidate downstream molecules of CCDC50. Two-way analysis of variance and unpaired Student t test were used in statistical analysis.

RESULTS

Pre-mRNA of CCDC50 generated two transcripts, full-length transcript (CCDC50-FL) and truncated transcript (CCDC50-S) with exon 6 skipped. CCDC50-S was overexpressed in ccRCC tissues and cell lines compared to noncancerous counterparts, but CCDC50-FL was only detected in noncancerous tissues and normal renal epithelial cells. Higher percent spliced-in index was associated with better survival in ccRCC patients. In vitro and in vivo functional experiments indicated that CCDC50-S transcript promoted the proliferation, migration, invasion and tumorigenesis of ccRCC, while CCDC50-FL exerted opposite tumor suppressive functions. Besides, we identified that heterogeneous nuclear ribonucleoprotein A1 (HnRNP A1) could promote the skipping of exon 6, which resulted in higher portion of CCDC50-S and oncogenic transformation. Moreover, zinc finger protein 395 (ZNF395) was identified as a downstream protein of CCDC50-S, and the interaction initiated oncogenic pathways which were involved in ccRCC progression.

CONCLUSIONS

Aberrant alternative splicing of CCDC50 is regulated by HnRNP A1 in ccRCC. This splicing event contributes to cancer progression through the downstream pathway involving ZNF395.

Identification of prognosis-related alternative splicing events in kidney renal clear cell carcinoma

Alternative splicing (AS) contributes to protein diversity by modifying most gene transcriptions. Cancer generation and progression are associated with specific splicing events. However, AS signature in kidney renal clear cell carcinoma (KIRC) remains unknown. In this study, genome‐wide AS profiles were generated in 537 patients with KIRC in the cancer genome atlas. With a total of 42 522 mRNA AS events in 10 600 genes acquired, 8164 AS events were significantly associated with the survival of patients with KIRC. Logistic regression analysis of the least absolute shrinkage and selection operator was conducted to identify an optimized multivariate prognostic predicting mode containing four predictors. In this model, the receptor‐operator characteristic curves of the training set were built, and the areas under the curves (AUCs) at different times were >0.88, thus indicating a stable and powerful ability in distinguishing patients' outcome. Similarly, the AUCs of the test set at different times were >0.73, verifying the results of the training set. Correlation and gene ontology analyses revealed some potential functions of prognostic AS events. This study provided an optimized survival‐predicting model and promising data resources for future in‐depth studies on AS mechanisms in KIRC.

Involvement of the CB2 cannabinoid receptor in cell growth inhibition and G0/G1 cell cycle arrest via the cannabinoid agonist WIN 55,212-2 in renal cell carcinoma

BACKGROUND

The anti-tumor properties of cannabinoids have been investigated in many in vitro and in vivo studies. Many of these anti-tumor effects are mediated via cannabinoid receptor types 1 and 2 (CB₁ and CB₂), comprising the endocannabinoid system (ECS). In this study, we investigated the ECS based on CB ₁ and CB ₂ receptor gene and protein expression in renal cell carcinoma (RCC) cell lines. In view of their further use for potential treatments, we thus investigated the roles of CB₁ and CB₂ receptors in the anti-proliferative action and signal transduction triggered by synthetic cannabinoid agonists [such as JWH-133 and WIN 55,212-2 (WIN-55)] in RCC cell lines.

METHODS

Human RCC cell lines were used for this study. The CB ₁ and CB ₂ gene expression levels were analyzed using real-time PCR. Flow cytometric, immunocytochemical and western blot analyses were performed to confirm CB₁ and CB₂ receptor protein expression. The anti-proliferative effects of synthetic cannabinoids were investigated on cell viability assay. The CB₁ and CB₂ receptors were blocked pharmacologically with the antagonists SR141716A and AM-630, respectively, to investigate the effects of the agonists JWH-133 and WIN-55. Cell cycle, apoptosis and LDH-based cytotoxicity were analyzed on cannabinoid-treated RCC cells.

RESULTS

The CB1 and CB2 genes expression was shown by real-time PCR and flow cytometric and western blot analysis indicating a higher level of CB₂ receptor as compared to CB₁ in RCC cells. Immunocytochemical staining also confirmed the expression of the CB₁ and CB₂ proteins. We also found that the synthetic cannabinoid agonist WIN-55 exerted anti-proliferative and cytotoxic effects by inhibiting the growth of RCC cell lines, while the CB₂ agonist JWH-133 did not. Pharmacologically blocking the CB1 and CB2 receptors with their respective antagonists SR141716A and AM-630, followed by the WIN-55 treatment of RCC cells allowed uncovering the involvement of CB2, which led to an arrest in the G0/G1 phase of the cell cycle and apoptosis.

CONCLUSIONS

This study elucidated the involvement of CB₂ in the in vitro inhibition of RCC cells, and future applications of CB₂ agonists in the prevention and management of RCC are discussed.

Global isoform-specific transcript alterations and deregulated networks in clear cell renal cell carcinoma

Extensive genome-wide analyses of deregulated gene expression have now been performed for many types of cancer. However, most studies have focused on deregulation at the gene-level, which may overlook the alterations of specific transcripts for a given gene. Clear cell renal cell carcinoma (ccRCC) is one of the best-characterized and most pervasive renal cancers, and ccRCCs are well-documented to have aberrant RNA processing. In the present study, we examine the extent of aberrant isoform-specific RNA expression by reporting a comprehensive transcript-level analysis, using the new kallisto-sleuth-RATs pipeline, investigating coding and non-coding differential transcript expression in ccRCC. We analyzed 50 ccRCC tumors and their matched normal samples from The Cancer Genome Altas datasets. We identified 7,339 differentially expressed transcripts and 94 genes exhibiting differential transcript isoform usage in ccRCC. Additionally, transcript-level coexpression network analyses identified vasculature development and the tricarboxylic acid cycle as the most significantly deregulated networks correlating with ccRCC progression. These analyses uncovered several uncharacterized transcripts, including lncRNAs FGD5-AS1 and AL035661.1, as potential regulators of the tricarboxylic acid cycle associated with ccRCC progression. As ccRCC still presents treatment challenges, our results provide a new resource of potential therapeutics targets and highlight the importance of exploring alternative methodologies in transcriptome-wide studies.

Identify clear cell renal cell carcinoma related genes by gene network

Clear cell renal cell carcinoma (ccRCC) is the most prominent type of kidney cancer in adults. The patients within metastatic ccRCC have a poor 5-year survival rate that is less than 10%. It is essential to identify ccRCC -related genes to help with the understanding of molecular mechanism of ccRCC. In this literature, we aim to identify genes related to ccRCC based on a gene network. We collected gene expression level data of ccRCC from the Cancer Genome Atlas (TCGA) for our analysis. We constructed a co-expression gene network as the first step of our study. Then, the network sparse boosting approach was performed to select the genes which are relevant to ccRCC. Results of our study show there are 15 genes selected from the all genes we collected. Among these genes, 7 of them have been demonstrated to play a key role in development and progression or in drug response of ccRCC. This finding offers clues of gene markers for the treatment of ccRCC.

TCF21 functions as a tumor suppressor in colorectal cancer through inactivation of PI3K/AKT signaling

Colorectal cancer (CRC) has become a major public health problem, ranking as the third most common type of cancer. Our previous study has revealed that TCF21 is frequently silenced by promoter hypermethylation in both CRC cell lines and primary CRC, with TCF21 methylation being significantly correlated with lymph node invasion. In this study, we further analyze the expression of TCF21 in CRC tissues and investigate the role of TCF21 in CRC in vitro and in vivo. We also explore the possible pathway regulated by TCF21. We thus demonstrate that decreased levels of TCF21 are associated with the pathological stage, clinical stage and lymph node metastasis, indicating a poor prognosis in CRC patients; overexpression of TCF21 inhibits cell proliferation, migration and invasion in the colorectal cell lines HCT116 and HT29. Furthermore, TCF21 functions as a tumor suppressor probably through inactivation of PI3K/AKT signaling and inhibition of MMPs. Our results suggest that enhancement of TCF21 levels may be a potential strategy to facilitate the prevention and treatment of CRC in the clinic.

Identification of prognostic genes in kidney renal clear cell carcinoma by RNA‑seq data analysis

The present study aimed to analyze RNA-seq data of kidney renal clear cell carcinoma (KIRC) to identify prognostic genes. RNA‑seq data were downloaded from The Cancer Genome Atlas. Feature genes with a coefficient of variation (CV) >0.5 were selected using the genefilter package in R. Gene co‑expression networks were constructed with the WGCNA package. Cox regression analysis was performed using the survive package. Furthermore, a functional enrichment analysis was conducted using Database for Annotation, Visualization and Integrated Discovery tools. A total of 533 KIRC samples were collected, from which 6,758 feature genes with a CV >0.5 were obtained for further analysis. The KIRC samples were divided into two sets: The training set (n=319 samples) and the validation set (n=214 samples). Subsequently, gene co‑expression networks were constructed for the two sets. A total of 12 modules were identified, and the green module was significantly associated with survival time. Genes from the green module were revealed to be implicated in the cell cycle and p53 signaling pathway. In addition, a total of 11 hub genes were revealed, and 10 of them (CCNA2, CDC20, CDCA8, GTSE1, KIF23, KIF2C, KIF4A, MELK, TOP2A and TPX2) were validated as possessing prognostic value, as determined by conducting a survival analysis on another gene expression dataset. In conclusion, a total of 10 prognostic genes were identified in KIRC. These findings may help to advance the understanding of this disease, and may also provide potential biomarkers for therapeutic development.

Cellular responses to HSV-1 infection are linked to specific types of alterations in the host transcriptome

Pathogen invasion triggers a number of cellular responses and alters the host transcriptome. Here we report that the type of changes to cellular transcriptome is related to the type of cellular functions affected by lytic infection of Herpes Simplex Virus type I in Human primary fibroblasts. Specifically, genes involved in stress responses and nuclear transport exhibited mostly changes in alternative polyadenylation (APA), cell cycle genes showed mostly alternative splicing (AS) changes, while genes in neurogenesis, rarely underwent these changes. Transcriptome wide, the infection resulted in 1,032 cases of AS, 161 incidences of APA, 1,827 events of isoform changes, and up regulation of 596 genes and down regulations of 61 genes compared to uninfected cells. Thus, these findings provided important and specific links between cellular responses to HSV-1 infection and the type of alterations to the host transcriptome, highlighting important roles of RNA processing in virus-host interactions.

Soluble Serum αKlotho Is a Potential Predictive Marker of Disease Progression in Clear Cell Renal Cell Carcinoma

Renal cell carcinoma (RCC) accounts for approximately 3% of adult malignancies, and clear cell RCC (ccRCC), that has a high metastatic index and high relapse rate, is the most common histological subtype. The identification of new biomarkers in ccRCC is fundamental for stratifying patients into prognostic risk groups and to guide therapy. The renoprotective antiaging gene, αKlotho, has recently been found to work as a tumor suppressor in different human cancers. Here, we evaluated αKlotho expression in tissue and serum of ccRCC patients and correlated it with disease progression. Tissue αKlotho expression was studied by quantitative RT-PCR and immunohistochemistry. In addition, soluble serum αKlotho levels were preoperatively measured in 160 patients who underwent nephrectomy for RCC with ELISA. Estimates of cancer-specific (CSS) and progression-free survival (PFS) were calculated according to the Kaplan-Meier method. Multivariate analysis was performed to identify the most significant variables for predicting CSS and PFS. αKlotho protein levels were significantly decreased in RCC tissues compared with normal tissues (P < 0.01) and the more advanced the disease, the more evident the down-regulation. This trend was also observed in serum samples. Statistically significant differences resulted between serum αKlotho levels and tumor size (P = 0.003), Fuhrman grade (P = 0.007), and clinical stage (P = 0.0004). CSS and PFS were significantly shorter in patients with lower levels of αKlotho (P < 0.0001 and P = 0.0004, respectively). At multivariate analysis low serum levels of αKlotho were independent adverse prognostic factors for CSS (HR = 2.11; P = 0.03) and PFS (HR = 2.18; P = 0.03).These results indicate that a decreased αKlotho expression is correlated with RCC progression, and suggest a key role of declining αKlotho in the onset of cancer metastasis.

Variation in genomic landscape of clear cell renal cell carcinoma across Europe

The incidence of renal cell carcinoma (RCC) is increasing worldwide, and its prevalence is particularly high in some parts of Central Europe. Here we undertake whole-genome and transcriptome sequencing of clear cell RCC (ccRCC), the most common form of the disease, in patients from four different European countries with contrasting disease incidence to explore the underlying genomic architecture of RCC. Our findings support previous reports on frequent aberrations in the epigenetic machinery and PI3K/mTOR signalling, and uncover novel pathways and genes affected by recurrent mutations and abnormal transcriptome patterns including focal adhesion, components of extracellular matrix (ECM) and genes encoding FAT cadherins. Furthermore, a large majority of patients from Romania have an unexpected high frequency of A:T>T:A transversions, consistent with exposure to aristolochic acid (AA). These results show that the processes underlying ccRCC tumorigenesis may vary in different populations and suggest that AA may be an important ccRCC carcinogen in Romania, a finding with major public health implications.