肾透明细胞癌预后相关miRNAs的生物信息学分析

目的寻找可作为肾透明细胞癌（ccRCC）生物标志物的miRNA,以及ccRCC与正常组织间miRNA差异表达情况。 方法利用TCGA数据库下载ccRCC中miRNA表达数据,分析肿瘤与正常组织间差异表达miRNA。使用Kaplan-Meier曲线对患者进行生存分析,筛选出表达情况与临床预后相关的miRNA。通过生物信息学对miRNA的靶基因进行预测,然后运用FunRich软件和ClueGO对靶基因进行GO和KEGG富集分析。 结果通过TCGA数据库分析发现,ccRCC较正常组织差异表达miRNA共54个,其中上调33个,下调21个。通过生存分析发现hsa-miR-21和hsa-miR-155与患者预后相关,P≤0.05。进一步通过Perl软件在Targetscan、miRDB、miRTarBase、miRPath这四个数据库中预测miRNA靶基因并将结果取交集,共发现129个靶基因。GO和KEGG分析结果表明,这些靶基因主要与转录因子活性、信号转导以及FoxO、TNF等信号通路密切相关。 结论通过生物信息学分析发现了ccRCC与正常组织的差异表达miRNA;其中hsa-miR-21和hsa-miR-155与患者总体生存率相关,并通过调控靶基因参与相关的信号通路进而影响ccRCC的发生发展进程,提示hsa-miR-21和hsa-miR-155可能是ccRCC潜在的生物标志物。     

肾透明细胞癌Caki-1细胞系差异表达基因的生物信息学分析

目的比较肾透明细胞癌Caki-1细胞系与正常肾上皮细胞系ASE-5063中的差异表达基因（DEGs）,寻找潜在的肾透明细胞癌特异性分子标志物。 方法利用GEO数据库自带的GEO2R在线分析工具分析基因芯片GSE78179,将筛选出的DEGs分别导入Metascape、STRING以及Cytoscape进行综合分析并筛选出核心基因。最后使用FunRich等软件对筛选出的核心基因进行GO和KEGG富集分析。 结果共筛选出562个DEGs,其中上调基因345个,下调基因217个。进一步使用MCODE筛选出36个关键基因,GO功能分析发现这些基因与细胞粘附分子活性、趋化因子活性、细胞通讯和信号转导等密切相关;KEGG通路富集结果则表明差异基因主要集中在趋化因子信号通路、TNF信号通路以及NF-κB信号通路等多种与肿瘤相关的通路上。 结论运用生物信息学方法筛选出肾透明细胞癌Caki-1细胞系中DEGs,其中数个核心基因广泛参与多种肿瘤的病理进程,但尚未在肾透明细胞癌有相关研究报道,提示其可能是治疗肾透明细胞癌的潜在靶点。     

The expression of ketohexokinase is diminished in human clear cell type of renal cell carcinoma

For identification and targeting of tumor-associated marker proteins, the proteome of clear cell type of renal cell carcinoma (RCC) and normal kidney tissues was analyzed by 2-DE. Ketohexokinase (also called fructokinase), which catalyzes the phosphorylation of fructose to fructose 1-phosphate, was identified by MALDI-TOF MS and found to be expressed at low rates in the renal tumor tissues. We found a decreased amount of ketohexokinase mRNA in RCC compared to that observed in the normal kidney tissues by Northern blot. The activity of ketohexokinase in 20 clear cell RCC specimens and the 20 corresponding normal kidneys was investigated, and its activity was shown to be approximately 1.4-fold lower in the RCC specimens than in the normal kidney. Ketohexokinase activity in tumor stage pT3 RCC was 1.5-fold lower than in pT1 RCC. The level of ketohexokinase activity in histological grade 3 RCC was 1.8-fold lower than that in grade 1 cancer. In addition, using in situ hybridization, it was revealed that ketohexokinase in the normal kidney tissue was confined to the proximal tubular epithelial cells, while the expression of ketohexokinase in RCC tissues was extremely low. Our research results show that the expression of human ketohexokinase was diminished in clear cell RCC.     

Integrated glycoproteomic characterization of clear cell renal cell carcinoma

1. Download : Download high-res image (269KB)
2. Download : Download full-size image

     

Characterization of an in vitro system of human renal papillary collecting duct cells

Summary We have developed an in vitro model of human papillary collecting duct cells isolated from cadaver kidneys using methods similar to those we previously reported for the isolation of human proximal tubule cells. To date we have isolated papillary collecting duct cells from 100 normal human kidneys. Papillae were dissected and digested in Cellgro containing 400 U/ml collagenase. Cells were plated on fibronectin-coated culture flasks at a density of 10⁴ live cells/ml in Cellgro supplemented with insulin and 10% fetal bovine serum. Confluent monolayers, which were able to withstand 600 mOsm for 8 h, were obtained within 10 to 15 d. Cells of primary isolates and first passages exhibited epithelial cell ultrastructure including cell junctions, microvilli, and cilia. A dark-brown reaction product was observed in these cells when stained by the immunoperoxidase method with peroxidase-labeled peanut lectin (Arachis hypogaea), which binds specifically to human distal tubule and collecting duct cells. These cells were negative for Factor-VIII (a marker for endothelial cells) and γ-glutamyltransferase (a marker for proximal tubule cells). High activities of the glycolytic enzyme pyruvate kinase and arginine vasopressin-stimulated cAMP production in these cells are consistent with a distal nephron origin. The results indicate that human collecting duct cells can be isolated and cultured to provide an in vitro system to probe pathogenetic mechanisms of potential nephrotoxins. Part of this work was presented at a Symposium of the Center for Alternatives to Animal Testing, April 4–5, 1989, Johns Hopkins Medical Institutions, Baltimore, MD 21205. This work was supported in part by grants R01-AI24179, PO1-A804393 for the Public Health Service, U.S. Department of Health and Human Services, and by a grant from the National Kidney Foundation, Baltimore, MD affiliate.     

Prognostic value of a gene signature in clear cell renal cell carcinoma

Renal cancer is a common urogenital system malignance. Novel biomarkers could provide more and more critical information on tumor features and patients’ prognosis. Here, we performed an integrated analysis on the discovery set and established a three-gene signature to predict the prognosis for clear cell renal cell carcinoma (ccRCC). By constructing a LASSO Cox regression model, a 3-messenger RNA (3-mRNA) signature was identified. Based on the 3-mRNA signature, we divided patients into high- and low-risk groups, and validated this by using three other data sets. In the discovery set, this signature could successfully distinguish between the high- and low-risk patients (hazard ratio (HR), 2.152; 95% confidence interval (CI),1.509–3.069; p < 0.0001). Analysis of internal and two external validation sets yielded consistent results (internal: HR, 2.824; 95% CI, 1.601–4.98; p < 0.001; GSE29609: HR, 3.002; 95% CI, 1.113–8.094; p = 0.031; E-MTAB-3267: HR, 2.357; 95% CI, 1.243–4.468; p = 0.006). Time-dependent receiver operating characteristic (ROC) analysis indicated that the area under the ROC curve at 5 years was 0.66 both in the discovery and internal validation set, while the two external validation sets also suggested good performance of the 3-mRNA signature. Besides that, a nomogram was built and the calibration plots and decision curve analysis indicated the good performance and clinical utility of the nomogram. In conclusion, this 3-mRNA classifier proved to be a useful tool for prognostic evaluation and could facilitate personalized management of ccRCC patients.     

Resistance to tyrosine kinase inhibitors in clear cell renal cell carcinoma: From the patient's bed to molecular mechanisms

The introduction of anti-angiogenic drugs especially tyrosine kinase inhibitors (TKIs) was a breakthrough in the treatment of renal cell carcinoma (RCC). Although TKIs have significantly improved outcome in patients with metastatic disease, the majority still develop resistance over time. Because different combinations and sequences of TKIs are tested in clinical trials, resistance patterns and mechanisms underlying this phenomenon should be thoroughly investigated. From a clinical point of view, resistance occurs either as a primary phenomenon (intrinsic) or as a secondary phenomenon related to various escape/evasive mechanisms that the tumor develops in response to vascular endothelial growth factor (VEGF) inhibition. Intrinsic resistance is less common, and related to the primary redundancy of available angiogenic signals from the tumor, causing unresponsiveness to VEGF-targeted therapies. Acquired resistance in tumors is associated with activation of an angiogenic switch which leads to either upregulation of the existing VEGF pathway or recruitment of alternative factors responsible for tumor revascularization. Multiple mechanisms can be involved in different tumor settings that contribute both to evasive and intrinsic resistance, and current endeavor aims to identify these processes and assess their importance in clinical settings and design of pharmacological strategies that lead to enduring anti-angiogenic therapies.     

肾透明细胞癌患者铁死亡相关基因的预后作用

探讨铁死亡相关基因在肾透明细胞癌患者中的表达及其预后价值。通过TCGA数据库下载KIRC的相关测序数据与检索到的铁死亡相关基因取交集,进行铁死亡相关基因的差异分析。之后利用单变量和多变量Cox回归分析,筛选具有预后价值的基因,构建预测患者生存情况的风险评分模型,并对模型进行验证。对高低风险组进行GO与KEGG通路富集,探讨风险差异的可能原因;通过ssGSEA分析,评估高低风险组间的免疫浸润情况。在KIRC患者的肿瘤组织和正常组织中,共得到21个差异的铁死亡相关基因;通过单因素Cox回归分析,获得 28 个与KIRC预后相关的基因;之后进行Lasso回归与多因素Cox回归分析,结果显示有10个基因被纳入模型,计算公式为:风险值(Risk score)=(0.024 5)×ALOX5表达值+(0.126 0)×CBS表达值+(0.199 5)×CD44表达值+(0.218 3)×CHAC1表达值+(-0.295 9)×HMGCR表达值+(0.036 7)×MT1G表达值+(0.061 4)×SLC7A11表达值+(-0.080 7)×FDFT1表达值+(0.160 3)×PEBP1表达值+(-0.220 5)×GOT1表达值。生存状态图表明,高风险组死亡病例数多于低风险组;ROC曲线表明风险评分模型具备一定预测能力;K-M生存分析显示,高风险组总体生存率低于低风险组(P=5.73×10^-13)。GO与KEGG富集分析提示,高低风险组间免疫情况及IL-17信号通路存在显著差异;进一步的ssGSEA富集显示,高低风险组间大部分免疫细胞的评分存在显著差异。基于铁死亡相关基因的预后风险评分模型可用于KIRC的预后预测,针对铁死亡相关基因设计靶点可能是治疗KIRC的一种新选择。     

Nephrocystin-1 interacts directly with Ack1 and is expressed in human collecting duct

Nephronophthisis is characterised by renal fibrosis, tubular basement membrane disruption and corticomedullary cyst formation leading to end stage renal failure. Mutations in NPHP1 account for the underlying genetic defect in 25% of patients with nephronophthisis. Loss of urine concentration ability may be an early feature of nephronophthisis. Using yeast-2-library screening with the SH3 domain of nephrocystin-1 as bait, we identify Ack1 as a novel interaction partner. This interaction is confirmed using exogenous over-expression followed by co-immunoprecipitation. Ack1 is an activated Cdc42-associated kinase, and like nephrocystin-1, is a known interactor of p130Cas. Nephrocystin-1 partially colocalises with Ack1 at cell-cell contacts in IMCD3 cells. In human kidney, nephrocystin-1 expression is limited to cell-cell junctions in renal collecting duct cells. These data define Ack1 as a novel interaction partner of nephrocystin-1 and implicate cell-cell junctions and the renal collecting duct in the pathology of nephronophthisis.     

MicroRNA-187, down-regulated in clear cell renal cell carcinoma and associated with lower survival,inhibits cell growth and migration though targeting B7-H3

Aberrantly expressed microRNAs (miRNAs) are frequently associated with the aggressive malignant behavior of human cancers, including clear cell renal cell carcinoma (ccRCC). Based on the preliminary deep sequencing data, we hypothesized that miR-187 may play an important role in ccRCC development. In this study, we found that miR-187 was down-regulated in both tumor tissue and plasma of ccRCC patients. Lower miR-187 expression levels were associated with higher tumor grade and stage. All patients with high miR-187 expression survived 5 years, while with low miR-187 expression, only 42% survived. Suppressed in vitro proliferation, inhibited in vivo tumor growth, and decreased motility were observed in cells treated with the miR-187 expression vector. Further studies showed that B7 homolog 3 (B7-H3) is a direct target of miR-187. Over-expression of miR-187 decreased B7-H3 mRNA level and repressed B7-H3-3′-UTR reporter activity. Knockdown of B7-H3 using siRNA resulted in similar phenotype changes as that observed for overexpression of miR-187. Our data suggest that miR-187 is emerging as a novel player in the disease state of ccRCC. miR-187 plays a tumor suppressor role in ccRCC.     

Orai1 and STIM1 are critical for cell migration and proliferation of clear cell renal cell carcinoma

The intracellular Ca²⁺ regulation has been implicated in tumorigenesis and tumor progression. Notably, store-operated Ca²⁺ entry (SOCE) is a major Ca²⁺ entry mechanism in non-excitable cells, being involved in cell proliferation and migration in several types of cancer. However, the expression and biological role of SOCE have not been investigated in clear cell renal cell carcinoma (ccRCC). Here, we demonstrate that Orai1 and STIM1, not Orai3, are crucial components of SOCE in the progression of ccRCC. The expression levels of Orai1 in tumor tissues were significantly higher than those in the adjacent normal parenchymal tissues. In addition, native SOCE was blunted by inhibiting SOCE or by silencing Orai1 and STIM1. Pharmacological blockade or knockdown of Orai1 or STIM1 also significantly inhibited RCC cell migration and proliferative capability. Taken together, Orai1 is highly expressed in ccRCC tissues illuminating that Orai1-mediated SOCE may play an important role in ccRCC development. Indeed, Orai1 and STIM1 constitute a native SOCE pathway in ccRCC by promoting cell proliferation and migration.     

Systematic analysis of alternative splicing signature unveils prognostic predictor for kidney renal clear cell carcinoma

There is growing evidence that alternative splicing (AS) plays an important role in cancer development. However, a comprehensive analysis of AS signatures in kidney renal clear cell carcinoma (KIRC) is lacking and urgently needed. It remains unclear whether AS acts as diagnostic biomarkers in predicting the prognosis of KIRC patients. In the work, gene expression and clinical data of KIRC were obtained from The Cancer Genome Atlas (TCGA), and profiles of AS events were downloaded from the SpliceSeq database. The RNA sequence/AS data and clinical information were integrated, and we conducted the Cox regression analysis to screen survival-related AS events and messenger RNAs (mRNAs). Correlation between prognostic AS events and gene expression were analyzed using the Pearson correlation coefficient. Protein-protein interaction analysis was conducted for the prognostic AS-related genes, and a potential regulatory network was built using Cytoscape (version 3.6.1). Meanwhile, functional enrichment analysis was conducted. A prognostic risk score model is then established based on seven hub genes (KRT222, LENG8, APOB, SLC3A1, SCD5, AQP1, and ADRA1A) that have high performance in the risk classification of KIRC patients. A total 46,415 AS events including 10,601 genes in 537 patients with KIRC were identified. In univariate Cox regression analysis, 13,362 survival associated AS events and 8,694 survival-specific mRNAs were detected. Common 3,105 genes were screen by overlapping 13,362 survival associated AS events and 8,694 survival-specific mRNAs. The Pearson correlation analysis suggested that 13 genes were significantly correlated with AS events (Pearson correlation coefficient >0.8 or <−0.8). Then, We conducted multivariate Cox regression analyses to select the potential prognostic AS genes. Seven genes were identified to be significantly related to OS. A prognostic model based on seven genes was constructed. The area under the ROC curve was 0.767. In the current study, a robust prognostic prediction model was constructed for KIRC patients, and the findings revealed that the AS events could act as potential prognostic biomarkers for KIRC.     

Identification of a three-miRNA signature as a novel potential prognostic biomarker in patients with clear cell renal cell carcinoma

Current studies suggest that some microRNAs (miRNAs) are associated with prognosis in clear cell renal cell carcinoma (ccRCC). In this paper, we aimed to identify a miRNAs signature to improve prognostic prediction for ccRCC patients. Using ccRCC RNA-Seq data of The Cancer Genome Atlas (TCGA) database, we identified 177 differentially expressed miRNAs between ccRCC and paracancerous tissue. Then all the ccRCC tumor samples were divided into training set and validation set randomly. Three-miRNA signature including miR130b, miR-18a, and miR-223 were constructed by the least absolute shrinkage and selection operator (LASSO) Cox regression model in training set. According to optimal cut-off value of three-miRNA signature risk score, all the patients could be classified into high-risk group and low-risk group significantly. Survival of patients was significantly different between two groups (hazard ratio, 5.58, 95% confidence interval, 3.17-9.80; P < 0.0001), and three-miRNA signature performed favorably prognostic and predictive accuracy. The results were further validated in the validation set and total set. Multivariate Cox regression analyses and subgroup analyses showed that three-miRNA signature was an independent prognostic factor. Two nomograms that integrated three-miRNA signature and three clinicopathological risk factors were constructed to predict overall survival and disease-free survival after surgery for ccRCC patients. Functional enrichment analysis showed the possible roles of three-miRNA signature in some cancer-associated biological processes and pathways. In conclusion, we developed a novel three-miRNA signature that performed reliable prognostic for patient survival with ccRCC, it might facilitate ccRCC patients counseling and individualize management.     

Candidate biomarkers in renal cell carcinoma

Although the human genome has been decoded, the knowledge about the pathogenesis of diseases including cancer is still limited. By focusing on renal cell carcinoma (RCC) we here summarize the data of various research groups analyzing the protein/peptide expression profiles of tumor lesions/cell lines or serum obtained from patients and respective controls. Different powerful approaches such as 2-DE, PROTEOMEX/SERPA/SPEARS, and T cell epitope discovery upon elution of MHC class I-bound peptides in combination with MS/LC-MS/MS revealed 500 differentially expressed proteins. The overlap in target recognition limits the pool to 299 unique protein identities, but only few thereof (12%) have been validated. The management, analysis, and interpretation of the distinct data sets derived from 27 publications required bioinformatic restructuring of the results. However, the comprehensive analysis of the results expands the knowledge about the pathophysiology of RCC in particular of the most prominent clear cell subtype by providing information on the differentially expressed proteins, their regulation status in RCC compared to normal kidney epithelium next to additional information on MHC-presented T cell epitopes and on serological targets. Despite the low number of validated differentially expressed proteins some of them might serve as candidate biomarkers for the diagnosis and/or as therapeutic targets.     

Sialylated glycoconjugates in chromophobe cell renal carcinoma compared with other renal cell tumors. Indication of its development from the collecting duct epithelium.

The present study was designed to shed light on the extraordinary histochemical properties of the chromophobe cell renal carcinoma detected by Hale's colloidal iron reaction. Special emphasis was laid on the lectin histochemical analysis of cytoplasmic glycoconjugates. Binding of peanut agglutinin (PNA) and Erythrina cristagalli agglutinin (ECA) after enzymatic release of sialic acid and direct binding of Dolichos biflorus agglutinin (DBA) correlates well with the expression of binding sites for Sambucus nigra agglutinin (SNA) and Maackia amurensis agglutinin (MAA) revealing abundant sialylated carbohydrate moieties within the cytoplasm. This characteristic binding pattern differs considerably from the faint staining observed in the majority of other renal carcinomas, thus confirming that the chromophobe cell renal carcinoma is a distinct entity. However, the lectin binding pattern of renal oncocytoma obviously resembles that of chromophobe carcinoma indicating a close relationship between these renal tumors. Detailed analysis of adjacent renal parenchyma revealed a lectin binding pattern quite similar to that described in the chromophobe carcinomas exclusively in the intercalated cells lining the collecting duct. This finding suggests that the chromophobe cell renal carcinoma originates from the collecting duct epithelium. The detection of small complexes consisting of altered epithelia which display the morphological characteristics of chromophobe carcinoma and the histochemical properties of intercalated cells probably indicates the emergence of preneoplastic lesions preceding the development of chromophobe carcinoma. Even though further studies are clearly needed to elucidate the physiological role of the cellular glycoconjugates detected, the present results already provide valuable insight into the histogenesis and pathogenesis of the chromophobe cell renal carcinoma.     

Single-cell analysis reveals metastatic cell heterogeneity in clear cell renal cell carcinoma

Renal cell carcinoma (RCC) is one of the leading causes of cancer-related death worldwide. Tumour metastasis and heterogeneity lead to poor survival outcomes and drug resistance in patients with metastatic RCC (mRCC). In this study, we aimed to assess intratumoural heterogeneity (ITH) in mRCC cells by performing a combined analysis of bulk data and single-cell RNA-sequencing data, and develop novel biomarkers for prognosis prediction on the basis of the potential molecular mechanisms underlying tumorigenesis. Eligible single-cell cohorts related to mRCC were acquired using the Gene Expression Omnibus (GEO) dataset to identify potential mRCC subpopulations. We then performed gene set variation analysis to understand the differential function in primary RCC and mRCC samples. Subsequently, we applied weighted correlation network analysis to identify coexpressing gene modules that were related to the external trait of metastasis. Protein-protein interactions were used to screen hub subpopulation-difference (sub-dif) markers (ACTG1, IL6, CASP3, ACTB and RAP1B) that might be involved in the regulation of RCC metastasis and progression. Cox regression analysis revealed that ACTG1 was a protective factor (HR < 1), whereas the other four genes (IL6, CASP3, ACTB and RAP1B) were risk factors (HR > 1). Kaplan-Meier survival analysis suggested the potential prognostic value of these sub-dif markers. The expression of sub-dif markers in mRCC was further evaluated in clinical samples by immunohistochemistry (IHC). Additionally, the genetic features of sub-dif marker expression patterns, such as genetic variation profiles, correlations with tumour-infiltrating lymphocytes (TILs), and targeted signalling pathway activities, were assessed in bulk RNA-seq datasets. In conclusion, we established novel subpopulation markers as key prognostic factors affecting EMT-related signalling pathway activation in mRCC, which could facilitate the implementation of a treatment for mRCC patients.