miRNA-mRNA 双重表达谱在前列腺癌分子调控机制中的作用

前列腺癌（prostate cancer，PCa）是男性最常见的癌症之一，也是西方男性癌症死亡的第2 大原因。有资料显示

，我国70 岁以上人群前列腺癌的发生率高达25%。早期前列腺癌不易被发现，大多数前列腺癌患者初诊时已经发生癌细胞转移，失去治愈的机会。在过去的几十年中，虽然前列腺癌的诊断技术已经得到了显著提高，而且前列腺特异性抗原（prostate-specific antigen，PSA）也被认为有助于提高前列腺癌的诊断率，但当前也有很多研究认为应用PSA 诊断前列腺癌并不完全成功。有研究表明，相比PSA，碱性磷酸酶（ALP）在评估治疗反应上的准确性更高

。miRNA（microRNA）是由18～25 个核苷酸组成的非编码RNA，参与基因表达的转录后调控

。研究认为

，体液中miRNA 表达差异与前列腺癌患者的预后相关，可能成为评估前列腺癌患者预后的分子标记物。本文通过对前列腺癌的mRNA 和miRNA的双重表达谱的生物信息学分析，探讨miRNA 和mRNA 在前列腺癌中的分子调控机制并构建调控网络，报道如下。

1 资料与方法

1.1 前列腺癌的miRNA 表达谱数据及差异表达分析通过GEO 数据库（https://www.ncbi.nlm.nih.gov/geo/）获取前列腺癌的miRNA 表达谱数据，数据平台为GSE8126，该数据中包括60 个前列腺癌组织样本和15 个正常前列腺组织样本。通过对数据的预处理（中心化，对数变换和标准化），获得326 个miRNA的表达值。采用matlab 生物信息工具箱中的

检验筛选差异表达的miRNA，定义差异表达的miRNA为：

＜0.05 且FDR＜0.01，其中Fold_change＞1.5 认为该miRNA 是上调的，Fold_change ＜1/1.5 认为该miRNA 是下调的。

上一秒还好得跟一个人似的，下一秒就掐起来，是田朵和小宁的常态。田朵和小宁都是独生子，在家里被宠惯了，虽然两个人结婚快两年了，可是他们的脾气一点也没收敛，经常为鸡毛蒜皮的小事吵架，吵急了的时候，还会把离婚挂在嘴边。

1.2 前列腺癌的mRNA 表达谱数据及差异表达分析通过GEO 数据库（https://www.ncbi.nlm.nih.gov/geo/）获取前列腺癌的mRNA 表达谱数据，数据平台为GSE6956，该数据中包括60 个前列腺癌组织样本和15 个正常前列腺组织样本。通过对数据的预处理（中心化，对数变换和标准化），获得13 787 个mRNA的表达值。差异表达mRNA的筛选标准为

＜0.05且FDR＜0.01；其中Fold_change＞1.5 认为该mRNA是上调的，Fold_change＜1/1.5 认为该mRNA 是下调的。

1.3.2 基于广义典则相关分析的miRNA-mRNA 调控 采用改进的广义典则相关分析筛选出所有上调和下调的miRNA 和mRNA，探讨miRNA 和mRNA的调控关系并构建调控网络。提取前10 个miRNA和前20 个mRNA 作为前列腺癌的特征生物标记，采用R 软件的mixOmics 软件包（http://r-project.org）进行分析。

2.3 miRNA 和mRNA的调控关系

2.1 miRNA 表达谱数据的差异表达分析 共发现117 个差异表达的miRNA，其中有15 个差异表达的miRNA 是下调的，如miR-521 和miR-527 等；有20个差异表达的miRNA 是上调的，如miR-511 和miR-526a 等。将上调和下调的差异表达的miRNA进行双向层次聚类，见图1。

1.3.1 基于Pearson 相关的miRNA-mRNA 反向调控从MicroCosm Targets 数据库（http://www.ebi.ac.uk/enright-srv/microcosm/htdocs/targets/v5/）下载筛选出的 Fold_change 最大的上调的 miRNA 和Fold_change 最小的下调的miRNA的靶基因，计算它们与靶基因的Person 相关系数并分析这两个miRNA 和它们的靶基因的反向调控关系。

1.3.3 miRNA-mRNA 网络构建 基于广义典则相关分析，确定典则相关系数的cutoff 值并根据cutoff值，保留获得的miRNA-mRNA 调控关系，根据该调控关系构建相应的调控网络。网络中的miRNA 如果调控的mRNA的数目较多，称为该miRNA 具有较高的度，为网络中的hub miRNA，认为其是与前列腺癌相关的重要生物标记。

2 结果

经济增长收敛性是指在封闭的经济条件下，对于一个有效经济范围的不同经济单位，初期的静态指标和其经济增长速度之间存在负相关关系，即落后地区比发达地区有更高的经济增长率 （刘强，2001）［21］。将其应用于能源领域，β收敛是指某地区能源强度的增长速度与其能源强度的初始水平呈负相关关系，分为绝对收敛和条件收敛两种。绝对β收敛是假设不同地区具有相同的外部条件，如经济发展程度、技术水平、产业结构等，随着时间的推移，所有地区的能源强度都将收敛于相同水平；条件β收敛是将外界因素纳入到模型中，探讨外界因素对能源强度收敛的影响。因而，本文通过建立条件β收敛的模型来探讨产业转移对能源强度收敛的影响。

2.3.2 基于广义典则相关的miRNA-mRNA 调控 对于提取出的全部差异表达的上调和下调的miRNAmRNA 双重表达谱数据集采用改进的广义典则相关分析，结果显示，该法可以有效区分前列腺癌样本和正常人群样本，见图3。根据第一主成分的权系数，提取出与前列腺癌相关度最高的前10 个miRNA 和前20 个mRNA，见图4A、图4B。在这些提取的miRNA 和mRNA 中，有些miRNA 和mRNA 被已有报道证实是与前列腺癌相关的，如发现MED13L的第一主成分权系数在前列腺癌患者中最高。有研究发现MED 家族的基因，如MED12的体细胞突变会导致前列腺癌的发生

。对于miRNA，其中miR-29b和miR-29b-1*的第一主成分权系数在正常人群中最高。此外，miR-488*的第一主成分权系数在前列腺癌患者中最高。按照典则相关系数的cutoff 为0.9的标准，保留了所有典则相关系数大于0.9的miRNA-mRNA 调控关系，绘制miRNA 和mRNA的调控网络图，其中miR-29b 和miR-29b-1*是网络中度最高的两个结点（度=21），这两个miRNAs 在网络中调控的基因数最多，见图5。

1.3 miRNA 和mRNA的调控关系

纳入标准：既往未行脊柱融合术；全脊柱X线片冠状面上Cobb角＞ 80°且bending位上柔韧度＜ 30%，或矢状面上Cobb角＞ 70°；自愿接受本研究治疗方案且签署知情同意书。排除标准：伴有颈椎不稳、寰枢椎脱位等不能耐受牵引的患者；严重药物依赖性疾病、心理疾病、严重脊柱外伤或脊髓损伤、心肺功能严重障碍等不能耐受治疗的患者；其他不适合因素。

2.3.1 基于Pearson 相关的miRNA-mRNA 反向调控Fold_change 最小的下调miRNA 为miR-137，在MicroCosm Targets 数据库中发现了其对应的1991个靶向基因，其中上调的基因为HOXC4、NPIPL3 和WSB1。miR-137 和HOXC4的相关系数为

=-0.5566，

=2.15E-07；miR-137 和NPIPL3的相关系数为

=-0.7366，

=4.98E-04；miR-137 和WSB1的相关系数为

=-0.7959，

=1.44E-17；miR-137 和这3 个基因呈反向调控关系。Fold_change 最大的上调miRNA为miR-367，在MicroCosm Targets 数据库中发现其对应的2045 个靶向基因，其中下调的基因分别为MEFA 和TNFAIP6。miR-367 和MEFA的相关系数为

=-0.7840，P=8.99E-17，miR-367 和TNFAIP6的相关系数为

=-0.7627，

=1.85E-15；miR-367 和MEFA 及TNFAIP6 呈反向调控关系。

2.2 mRNA 表达谱数据的差异表达分析 共识别出9054 个差异表达的mRNA，其中有35 个差异表达的mRNA 是下调的，有36 个差异表达的mRNA 是上调的。双向层次聚类结果见图2，mRNA 表达谱能够很好的区分不同组织起源的样本。

3 讨论

前列腺癌是男性生殖系最常见的恶性肿瘤，研究表明分子调控机制在前列腺癌的预后和诊断中有着非常重要的作用。本研究对前列腺癌的mRNA 和miRNA的双重表达谱进行了生物信息学分析，采用了差异表达分析提取了与前列腺癌相关的上调和下调的miRNA 和mRNA，并通过Pearson 相关分析和广义典则相关分析探讨了miRNA 和mRNA的分子调控机制。广义典则相关分析可用于特征降维，特征融合和相关分析等。将广义典则相关分析应用于特征miRNA 和mRNA的提取，其中特别结合了miRNA 和mRNA的靶向失调关系作为预设矩阵，并根据miRNA 或mRNA 在第一个成分的权系数提取与前列腺癌相关的特征miRNA 和特征mRNA。

本病案对患者进行了综合性治疗，尤其对基础病进行有效干预，患者加用胰岛素治疗，使血糖维持在相对正常的水平，提高了机体抗感染能力，为我们后续治疗提供了保障。

除了对病毒分离培养与鉴定方法之外，还可以采用免疫电镜法，这种检测方法拥有简便、直观、快速、准确的优点，但是由于检测过程中需要使用价格昂贵的电镜设备，所以诊断费用比较高。此种方法不适合进行大规模的诊断。免疫电镜法就是可以通过电镜能够清晰的观察到病毒结构、形态以及立体结构等特征，然后可以对疾病作出判断。将患病仔猪的小肠或者病毒细胞培养物用固定液进行固定，然后用机器切成超薄的切片，在电镜下进行观察病毒的结构、形态。这种方法能够有效区别猪流行性腹泻病毒和猪传染性胃肠炎病毒，因为免疫电镜法不仅可以用已知的猪流行性腹泻病毒高免血清检测未知抗体，而且还能够借助已知的猪流行性腹泻病毒抗原检测未知抗体[3]。

本研究发现，MED 家族基因MED13L的第一主成分权系数在前列腺癌患者中最高。有研究发现，在局限性前列腺癌患者的异质性肿瘤病灶中检测到体细胞肿瘤变异，发现cfDNA 突变负荷与疾病严重程度相关。在cfDNA 中发现了体细胞组织的改变，其中就包括MED 家族中MED12的非同义变异

。此外，TRRAP 是本次识别的与前列腺癌高度相关的基因。TRRAP 曾在阴茎鳞状细胞癌

和卵巢癌中检测到突变

。在通过外显子测序及Sanger 测序探索66 个家族中140 例患有前列腺癌的日本患者的生殖系变异研究中，除了已知的易感基因BRCA2 和HOXB13 外，还发现了TRRAP 变体以及BRCA2、HOXB13 和TRRAP的共同变体

。

了解癌症中miRNA的调控机制对于理解其生理作用以及癌症相关功能障碍和失调的影响至关重要

。近年来的研究表明，miRNAs的异常表达与前列腺癌的发生、侵袭、转移及预后密切相关。虽然已经报道了许多miRNAs 可作为前列腺癌治疗的新生物标记物，但报告的数据不一致，仍需进一步实验

。因此，研究前列腺癌相关的miRNA 及其调控的靶向基因对于识别前列腺癌相关的生物标记及这些生物标记在前列腺癌的临床应用有着十分重要的价值。

本研究发现，miR-29b 与前列腺癌高度关联，并在miRNA-mRNA 分子调控网络中调控的基因数最多。目前，已有文献证实miR-29b 可以阻止前列腺癌的发生和发展。Zhu C等

发现前列腺癌患者肿瘤组织中miR-29b的表达明显低于癌旁组织，有望成为前列腺癌诊断的新生物标志物。Ivanovic RF等

研究认为，恢复miR-29b 水平的治疗可能是控制前列腺癌转移的有效方案。Sur S等

在瘤内注射模拟miR-29b，显著抑制了裸鼠前列腺癌异种移植瘤的生长，表明miR-29b 可以作为前列腺癌治疗的潜在分子，为未来开发前列腺癌的药物靶向提供帮助。Vladimir AV等

在年轻男性前列腺癌患者的分子特征研究中发现，miR-29b的表达与老年男性前列腺癌患者显著不同，这项研究为前列腺癌亚型的分子机制提供了参考。Kato M等

的研究表明，miR-29b 与其他五种miRNA 协同调节促进转移的基因（LOXL2），揭示了前列腺癌发病机制的新分子机制。此外，miR-488 是本研究中发现的与前列腺癌相关度最高的前10 个miRNA。miR-488 在多种肿瘤中作为抑癌基因发挥作用，如乳腺癌、卵巢癌和结直肠癌等

。Wang J等

研究发现，miR-488 在前列腺癌患者肿瘤组织中的表达低于正常组织，其通过靶向PFKFB3 可以抑制前列腺癌细胞的增殖和糖酵解，因此可能会作为潜在的前列腺癌新的治疗靶点。本研究也存在一定的局限性。如数据库中的前列腺癌相关的miRNA-mRNA 双重表达谱数据集相对很少，因此不能保证分析结果的稳定性；此外，本研究结果是基于miRNA-mRNA 双重表达谱数据的生物信息学分析，识别出的特征miRNA 和mRNA 还需要应用分子生物实验进行验证。

综上所述，MED 家族基因MED13L的第一主成分权系数在前列腺癌患者中最高，TRRAP 是与前列腺癌高度相关的基因；miR-29b 与前列腺癌高度关联，并在miRNA-mRNA 分子调控网络中调控的基因数最多，而miR-488 是与前列腺癌相关度最高的前10 个miRNA。

[1] 罗于杰,蓝建华,黄国华.miR-29a 过表达对人前列腺癌细胞凋亡的调控机制研究[J].中国煤炭工业医学杂志,2021,24(1):15-20.

[2] Wenter V,Herlemann A,Fendler WP,et al.Radium-223 for primary bone metastases in patients with hormone-sensitive prostate cancer after radical prostatectomy[J].Oncotarget,2017,8(27):2021-2026.

[3] Khan S,Ayub H,Khan T,et al.MicroRNA biogenesis,gene silencing mechanisms and role in breast,ovarian and prostate cancer[J].Biochimie,2019,167:12-24.

[4] McDonald AC,Vira M,Walter V,et al.Circulating microRNAs in plasma among men with low-grade and high-grade prostate cancer at prostate biopsy [J].Prostate,2019,79 (9):961-968.

[5] Kämpjärvi K,Kim NH,Keskitalo S,et al.Somatic MED12 mutations in prostate cancer and uterine leiomyomas promote tumorigenesis through distinct mechanisms[J].The Prostate,2016,76(1):22-31.

[6] Chen E,Cario CL,Leong L,et al.Cell-Free DNA Detection of Tumor Mutations in Heterogeneous,Localized Prostate Cancer Via Targeted,Multiregion Sequencing [J].JCO precision oncology,2021,5:PO.20.00428.

[7] Wang Y,Wang K,Chen Y,et al.Mutational landscape of penile squamous cell carcinoma in a Chinese population[J].International Journal of Cancer,2019,145(5):1280-1289.

[8] Kang KT,Kwon YW,Kim DK,et al.TRRAP stimulates the tumorigenic potential of ovarian cancer stem cells [J].BMB Reports,2018,51(10):514-519.

[9] Takahide H,Hiroshi M,Hirofumi N,et al.Germline Variants of Prostate Cancer in Japanese Families [J].PLoS One,2017,11(10):e0164233.

[10] Zainab AS,Siu SS,Choijamts M,et al.Regulatory Mechanism of MicroRNA Expression in Cancer [J].International Journal of Molecular Sciences,2020,21(5):1273.

[11] Sharma N,Baruah MM.The microRNA signatures: aberrantly expressed miRNAs in prostate cancer [J].Clin Transl Oncol,2019,21(2):126-144.

[12] Zhu C,Hou X,Zhu J,et al.Expression of miR -30c and miR-29b in prostate cancer and its diagnostic significance[J].Oncology Letters,2018,16(3):3140-3144.

[13] Ivanovic RF,Viana NI,Morais DR,et al.miR-29b enhances prostate cancer cell invasion independently of MMP-2 expression[J].Cancer Cell International,2018,18(1):18.

[14] Sur S,Steele R,Shi X,et al.miRNA-29b Inhibits Prostate Tumor Growth and Induces Apoptosis by Increasing Bim Expression[J].Cells,2019,8(11):455.

[15] Vladimir AV,Rafael P,Beatriz AW,et al.microRNA Expression Profiling in Young Prostate Cancer Patients [J].Journal of Cancer,2020,11(14):4106-4114

[16] Kato M,Kurozumi A,Goto Y,et al.Regulation of metastasispromoting LOXL2 gene expression by antitumor microRNAs in prostate cancer [J].Journal of Human Genetics,2017,62(1):123-132.

[17] Masuda T,Shinden Y,Noda M,et al.Circulating Pre-microRNA-488 in Peripheral Blood Is a Potential Biomarker for Predicting Recurrence in Breast Cancer [J].Anticancer Res,2018,38(8):4515-4523.

[18] Guo JY,Wang XQ,Sun LF.MicroRNA-488 inhibits ovarian cancer cell metastasis through regulating CCNG1 and p53 expression[J].Eur Rev Med Pharmacol Sci,2020,24(6):2902-2910.

[19] Deng X,Li D,Ke X,et al.Mir-488 alleviates chemoresistance and glycolysis of colorectal cancer by targeting PFKFB3[J].J Clin Lab Anal,2021,35(1):e23578.

[20] Wang J,Li X,Xiao Z,et al.MicroRNA-488 inhibits proliferation and glycolysis in human prostate cancer cells by regulating PFKFB3[J].FEBS Open Bio,2019,9(10):1798-1807.