Critical role of long non-coding RNA-microRNA-mRNA regulatory network in the development of advanced bladder cancer
杜洋,庞明睿,袁静萍,刘雯,喻希,刘修恒,王磊
DOI: https://doi.org/10.3760/cma.j.cn421213-20200718-01248
2021-01-01
Abstract:Objective:To explore the RNA associated with poor prognosis in the development of advanced bladder cancer (BCa), we constructed and validated an long non-coding RNA-microRNA-mRNA(lncRNA-miRNA-mRNA) network, so as to investigate its molecular regulatory role in the development of BCa.Methods:Using transcriptome profiles from The Cancer Genome Atlas and Gene Expression Omnibus databases, 431 and 492 cases of BCa tissue samples were collected. Clinical stage Ⅲ and Ⅳ, high grade, high grade with stage ≥pathological stage 2(pT2) were classified as advanced clinical stage BCa (ACS-BCa). Clinical stage Ⅱ andⅠ, low grade were classified as low clinical stage BCa (LCS-BCa). We performed differential expression (DE) analysis, weighted gene co-expression network analysis, functional enrichment analysis, survival analysis, prediction of miRNA targeting, and Pearson correlation analysis. Statistically, t′ test, Fisher′s exact test, Pearson correlation analysis, Log-rank test and Flase discovery rate (FDR) calculated with benjamini hochberg (BH) method were used. FDR < 0.05 and P<0.05 were used as sifting thresholds. Subsequently, 22 cases of BCa tissue samples were collected. High grade and muscle-invasion with stage ≥T2, high grade or invasion with stage TaT1 were classified as ACS-BCa, and low grade with stage TaT1 was classified as LCS-BCa. The real time quantitative polymerane chain reaction (qPCR), Western blotting, and immunohistochemistry were carried out. The t′ test was used to compare the statistically significant differences between experimental groups. Results:Totally, 8 lncRNAs-28 mRNAs and 8 miRNAs-28 mRNAs pairs shared similar expression patterns, constituting a core competing endogenous RNA(ceRNA) regulatory network related to the invasion, progression, metastasis and poor prognosis of ACS-BCa. Respectively, the RNA relative expression fold changes of collagen type Ⅵ α1(COL6A1), cadherin type Ⅱ (CDH11), a disintegrin and metalloproteinase with thrombaspondin motif type I(ADAMTS12), long non-coding RNA 01705、01929(LINC01705, LINC01929) in ACS-BCa tissues were significantly greater than those in LCS-BCa tissues (-1.213 vs. -4.212, t=4.170, P<0.01; -10.75 vs. -14.64, t=4.602, P<0.01; -7.594 vs. -11.93, t=5.152, P<0.01; -10.18 vs. 16.55, t=4.493, P<0.01; -11.62 vs. -15.08, t=2.869, P<0.05). The mean density of COL6A1, CDH11 and ADAMTS12 proteins in ACS-BCa tissues was significantly higher than that in early bladder cancer tissues (0.213 8 vs. 0.090 9, t=2.302, P<0.05; 0.183 9 vs. 0.078 9, t=5.992, P<0.01; 0.279 1 vs. 0.147 8, t=4.527, P<0.01). The RNA relative expression fold changes of homo sapiens microRNA 6784-5p、6875-5p、128-2(hsa-miR-6875-5p, hsa-miR-6784-5p and hsa-miR-128-2) in BCa tissues were significantly lower than those in normal tissues (-15.78 vs. -12.41, t=4.995, P<0.01; -15.53 vs. -12.68, t=3.226, P<0.05; -13.79 vs. -11.43, t=3.400, P<0.01). The relative gray values of COL6A1 protein in BCa cell line T24 were significantly lower than those in normal bladder epithelial cell line SV-HUV-1 (0.677 vs. 1.000, t=7.584, P<0.01). The relative gray values of CDH11 protein in BCa cell lines T24 and EJ were significantly lower than those in normal bladder epithelial cell line SV-HUV-1 (0.656 vs. 1.000, t=6.056, P<0.01). Conclusion:lncRNAs, miRNAs and mRNAs in the ceRNA network tend to play a critical role in the development of BCa.