基于生物信息学的乳腺癌患者预后相关miRNA预测模型构建与深度剖析

基于生物信息学的乳腺癌患者预后相关miRNA预测模型构建与深度剖析一、引言1.1研究背景与意义1.1.1乳腺癌现状乳腺癌是全球女性最常见的恶性肿瘤之一，严重威胁着女性的生命健康和生活质量。近年来，其发病率呈逐年上升趋势。根据世界卫生组织国际癌症研究机构（IARC）发布的2020年全球最新癌症负担数据，2020年乳腺癌取代肺癌成为全球发病率第一的癌症，全球新发乳腺癌病例达226万例，占全球女性新发癌症总数的24.5%，死亡病例68万例，位居世界第五。在中国，由于人口基数较大，乳腺癌的发病和死亡人数也不容小觑，2020年国内新发乳腺癌病例约为41.6万，癌症死亡人数约为11.7万。乳腺癌的发病与多种因素相关，包括遗传因素、生活方式、激素水平等。例如，具有乳腺癌家族遗传史的女性，其发病风险显著高于普通人群；现代女性生活节奏加快，长期熬夜、精神持续紧张、作息不规律、饮食不健康等问题突出，导致自身抵抗力下降，也增加了患癌风险；另外，晚婚晚育、生育次数减少、哺乳时期过短，以及超重和肥胖、缺乏运动等，都是乳腺癌的诱发因素。乳腺癌分为腔面A、腔面B、HER-2阳性和三阴性四种亚型，其中三阴性乳腺癌因缺乏雌激素受体（ER）、孕激素受体（PR）和人表皮生长因子受体2（HER2）表达，治疗手段相对有限，预后较差，被认为是最“毒”的亚型。乳腺癌早期症状主要表现为乳腺肿块、乳头溢液、乳头内陷等，但这些症状往往不明显，容易被忽视，一旦发现时可能已处于中晚期，错过了最佳治疗时机。目前，乳腺癌的检查手段主要有钼靶、彩超等，通常将钼靶和B超相结合，以提高确诊率。1.1.2miRNA与癌症的关系miRNA（microRNA）是一类长度约20-25个核苷酸的非编码小RNA分子，其在基因表达调控中发挥着至关重要的作用。2024年诺贝尔生理学或医学奖授予了美国科学家维克托・安布罗斯（VictorAmbros）和加里・鲁弗肯（GaryRuvkun），以表彰他们在microRNA发现及其在转录后基因调控中作用研究方面的贡献。miRNA主要通过与靶mRNA的互补配对结合，抑制mRNA的翻译过程，或者促使mRNA降解，从而实现对基因表达的负向调控。每个miRNA可以有多个靶基因，而几个miRNA也可以调节同一个基因，这种复杂的调控网络使得miRNA能够精细地调控细胞的各种生理过程。越来越多的研究表明，miRNA与癌症的发生、发展密切相关。在癌症的发生过程中，miRNA的表达水平常常发生异常改变，一些miRNA可作为癌基因促进肿瘤的生长和转移，而另一些则可作为抑癌基因抑制肿瘤的发展。例如，在乳腺癌中，miR-21表达上调，它可以通过抑制其靶基因的表达，促进乳腺癌细胞的增殖、侵袭和转移；而miR-34家族成员（miR-34a、miR-34b和miR-34c）作为肿瘤抑制miRNA，在大多数癌症类型中，由于肿瘤抑制因子p53经常发生突变、缺失或下调，miR-34a也相应下调，从而失去对肿瘤细胞的抑制作用，导致癌症的发生发展。此外，miRNA还参与了癌症的耐药过程，一些miRNA的异常表达与乳腺癌细胞对化疗药物、靶向药物的耐药性相关，影响着癌症的治疗效果和患者的预后。1.1.3构建miRNA预测模型的意义准确评估乳腺癌患者的预后情况，对于制定个性化的治疗方案、提高患者的生存率和生活质量具有重要意义。传统的乳腺癌预后评估主要依赖于临床病理特征，如肿瘤大小、淋巴结转移情况、病理分期等，但这些指标存在一定的局限性，无法全面准确地预测患者的预后。近年来，随着对miRNA与癌症关系研究的深入，发现miRNA在乳腺癌的预后评估中具有巨大的潜力。构建基于miRNA的预测模型，可以综合考虑多个miRNA的表达水平及其相互作用，更全面、准确地评估乳腺癌患者的预后风险。通过分析大量乳腺癌患者的miRNA表达谱与临床预后数据之间的关联，筛选出与预后密切相关的miRNA作为生物标志物，进而建立预测模型，能够为临床医生提供更有价值的预后信息。例如，通过该模型可以预测患者的复发风险、生存时间等，对于高风险患者，可以加强术后的辅助治疗和随访监测，采取更积极的治疗策略，以降低复发风险，提高生存率；对于低风险患者，则可以避免过度治疗，减少不必要的医疗负担和毒副作用，提高患者的生活质量。此外，miRNA预测模型还有助于深入了解乳腺癌的发病机制和生物学行为，为开发新的治疗靶点和治疗方法提供理论依据，推动乳腺癌的精准医疗发展。1.2国内外研究现状1.2.1miRNA在乳腺癌预后中的研究进展在乳腺癌预后研究领域，miRNA作为极具潜力的生物标志物，近年来吸引了众多国内外学者的广泛关注，相关研究成果丰硕。国外方面，诸多研究深入挖掘了特定miRNA与乳腺癌预后的紧密联系。例如，有研究表明miR-21在乳腺癌组织中呈现高表达状态，并且这种高表达与患者的不良预后显著相关。深入探究其机制发现，miR-21能够通过靶向作用于关键的抑癌基因，如PTEN等，抑制其表达，进而激活PI3K/Akt信号通路，最终促进乳腺癌细胞的增殖、侵袭与转移，严重影响患者的预后。还有研究聚焦于miR-155，发现其在乳腺癌中表达上调，可通过调控多个靶基因，如SOCS1等，参与免疫调节和肿瘤细胞的增殖、转移过程，与乳腺癌的不良预后密切相关。在国内，科研人员也在积极探索miRNA在乳腺癌预后中的价值。有研究从中国人群的乳腺癌样本出发，筛选出了一系列与预后相关的miRNA。如发现miR-125b在乳腺癌组织中低表达，且低表达水平与患者的远处转移和不良预后相关。进一步研究证实，miR-125b可通过靶向调控相关基因，抑制乳腺癌细胞的迁移和侵袭能力，发挥抑癌作用，其表达缺失可能导致乳腺癌细胞的恶性行为增强，从而影响患者的预后。另有研究对中国乳腺癌患者的miRNA表达谱进行分析，发现miR-34a的表达水平与乳腺癌的病理分期、淋巴结转移等临床病理特征相关，低表达的miR-34a预示着患者的不良预后，这可能是因为miR-34a作为p53的下游靶点，在肿瘤抑制过程中发挥重要作用，其表达下调会削弱对肿瘤细胞的抑制作用，进而影响患者的生存预后。此外，部分研究还关注到外泌体miRNA在乳腺癌预后中的作用。外泌体是一种由细胞分泌的纳米级囊泡，可携带多种生物活性分子，包括miRNA。研究发现，乳腺癌患者血清外泌体中的miR-221和miR-222表达水平升高，与肿瘤的复发和转移密切相关，有望作为乳腺癌预后监测的新型标志物。这是由于外泌体miRNA能够在细胞间传递信息，影响受体细胞的生物学行为，促进肿瘤的进展，因此检测外泌体miRNA的表达水平可以为乳腺癌患者的预后评估提供新的视角和依据。1.2.2乳腺癌预后模型的研究现状目前，乳腺癌预后模型的构建方法丰富多样，每种方法都各有优劣。传统的预后模型主要基于临床病理因素构建，如TNM分期系统，它依据肿瘤大小（T）、淋巴结转移情况（N）和远处转移（M）对乳腺癌进行分期，为临床医生提供了初步的预后评估依据。TNM分期系统在临床实践中应用广泛，具有直观、简单的优点，能够帮助医生快速了解患者的病情严重程度。然而，该系统也存在一定局限性，它无法全面考虑肿瘤的生物学特性，对于一些分期相同但预后差异较大的患者，难以准确预测其预后情况。随着分子生物学技术的飞速发展，基于基因表达谱的预后模型逐渐兴起。这类模型通过分析乳腺癌组织中的基因表达数据，筛选出与预后相关的基因，进而构建模型来预测患者的预后。例如，OncotypeDX基因检测通过检测21个基因的表达水平，计算复发评分，用于评估早期乳腺癌患者的复发风险和化疗获益。这种模型能够从分子层面揭示肿瘤的生物学行为，提供更精准的预后信息，有助于医生制定个性化的治疗方案。但它也存在一些不足之处，如检测成本较高，对检测技术和设备要求严格，限制了其在一些资源有限地区的广泛应用。近年来，机器学习算法在乳腺癌预后模型构建中得到了广泛应用。支持向量机（SVM）、随机森林（RF）、人工神经网络（ANN）等算法被用于分析大量的临床和分子数据，构建预后模型。SVM算法通过寻找一个最优的分类超平面，将不同预后的患者进行区分，具有良好的泛化能力和较高的预测准确性。随机森林算法则通过构建多个决策树，并对其结果进行综合，能够有效处理高维数据和特征选择问题，提高模型的稳定性和预测性能。人工神经网络算法具有强大的非线性拟合能力，能够学习复杂的数据模式，但存在训练时间长、容易过拟合等问题。这些机器学习算法在乳腺癌预后模型构建中展现出了巨大的潜力，能够综合考虑多种因素，提高预后预测的准确性，但它们也面临着数据质量、模型可解释性等挑战。例如，机器学习模型通常需要大量高质量的数据进行训练，数据的缺失、噪声等问题可能会影响模型的性能；同时，由于机器学习模型的复杂性，其决策过程往往难以解释，这在一定程度上限制了其在临床实践中的应用。1.3研究目的与内容1.3.1研究目的本研究旨在通过对乳腺癌患者的miRNA表达谱进行深入分析，构建基于miRNA的乳腺癌患者预后预测模型，并对该模型进行验证和评估，以提高乳腺癌患者预后预测的准确性。具体目标如下：从大量的miRNA中筛选出与乳腺癌患者预后密切相关的miRNA，确定其作为生物标志物的潜力，为乳腺癌预后评估提供新的指标。运用合适的数据分析方法和机器学习算法，构建基于miRNA的乳腺癌患者预后预测模型，实现对患者预后风险的精准预测。对构建的预测模型进行内部和外部验证，评估其预测性能，包括准确性、敏感性、特异性等指标，确保模型的可靠性和稳定性。分析预测模型中miRNA与乳腺癌临床病理特征之间的关系，探讨miRNA在乳腺癌发生、发展过程中的作用机制，为乳腺癌的个性化治疗和精准医疗提供理论依据。1.3.2研究内容本研究将围绕乳腺癌患者预后相关miRNA预测模型的构建及分析展开，具体研究内容包括以下几个方面：数据收集与预处理：收集乳腺癌患者的临床病理资料，包括年龄、肿瘤大小、淋巴结转移情况、病理分期、分子分型等信息，同时收集患者的肿瘤组织样本，提取RNA并进行质量检测，确保数据的可靠性和完整性。运用生物信息学方法对收集到的miRNA表达数据进行预处理，包括数据标准化、归一化、缺失值处理等，消除数据中的噪声和误差，为后续分析提供高质量的数据基础。差异表达miRNA筛选：利用生物信息学工具和统计学方法，比较乳腺癌患者与正常对照组之间的miRNA表达谱差异，筛选出在乳腺癌组织中显著差异表达的miRNA。进一步分析差异表达miRNA与乳腺癌患者预后指标（如总生存期、无病生存期等）之间的相关性，确定与预后密切相关的miRNA，为模型构建提供候选生物标志物。预测模型构建：运用机器学习算法，如支持向量机（SVM）、随机森林（RF）、逻辑回归（LR）等，以筛选出的预后相关miRNA作为特征变量，构建乳腺癌患者预后预测模型。通过交叉验证等方法对模型进行优化和调参，提高模型的预测性能和泛化能力。模型验证与评估：采用内部验证和外部验证相结合的方式对构建的预测模型进行验证。内部验证使用同一数据集进行多次交叉验证，评估模型在训练数据上的性能表现；外部验证则使用独立的数据集对模型进行测试，检验模型在不同数据上的泛化能力。通过计算准确率、召回率、F1值、受试者工作特征曲线（ROC）下面积等指标，全面评估模型的预测性能，判断模型的优劣。模型分析与机制探讨：分析预测模型中miRNA与乳腺癌临床病理特征之间的关联，探讨miRNA在乳腺癌发生、发展及预后中的作用机制。例如，研究miRNA对乳腺癌细胞增殖、凋亡、侵袭和转移等生物学行为的影响，以及miRNA与乳腺癌相关信号通路的交互作用，为乳腺癌的治疗提供新的靶点和思路。1.4研究方法与技术路线1.4.1研究方法生物信息学分析：运用生物信息学工具，如R语言、Python、GraphPadPrism9、GEPIA2、Ualcan、Kaplan-MeierPlotter、miRNet、TargetScan、miRTarBase等，对乳腺癌患者的miRNA表达数据进行分析。利用R语言和Python进行数据预处理，包括数据标准化、归一化、缺失值处理等；使用GraphPadPrism9对数据进行可视化展示，绘制柱状图、箱线图、散点图等，直观呈现数据特征；通过GEPIA2和Ualcan数据库分析miRNA与乳腺癌患者临床病理特征之间的相关性；借助Kaplan-MeierPlotter数据库进行生存分析，评估miRNA对乳腺癌患者总生存期和无病生存期的影响；运用miRNet、TargetScan、miRTarBase等数据库预测miRNA的靶基因，并进行基因功能富集分析，探讨miRNA的作用机制。统计分析：采用SPSS26.0、R语言等统计软件进行数据分析。对于计量资料，如miRNA表达水平，采用t检验或方差分析比较不同组之间的差异；对于计数资料，如患者的生存状态，采用卡方检验分析miRNA与生存情况之间的关联。通过单因素和多因素Cox回归分析筛选出与乳腺癌患者预后独立相关的miRNA，并计算风险比（HR）和95%置信区间（CI），评估miRNA对预后的影响程度。运用受试者工作特征曲线（ROC）评估模型的预测性能，计算曲线下面积（AUC），AUC越接近1，表明模型的预测准确性越高。机器学习算法：运用支持向量机（SVM）、随机森林（RF）、逻辑回归（LR）等机器学习算法构建乳腺癌患者预后预测模型。在Python的scikit-learn库中实现这些算法，通过交叉验证等方法对模型进行优化和调参，选择最优的模型参数，以提高模型的预测性能和泛化能力。例如，对于SVM算法，通过调整核函数类型（如线性核、高斯核等）、惩罚参数C等，寻找最佳的分类超平面；对于RF算法，调整决策树的数量、最大深度等参数，提高模型的稳定性和准确性。利用网格搜索、随机搜索等方法进行参数搜索，以找到最优的模型参数组合。1.4.2技术路线本研究的技术路线如图1-1所示：数据获取：从TCGA数据库、GEO数据库等公共数据库下载乳腺癌患者的miRNA表达谱数据及相应的临床病理资料，包括患者的年龄、性别、肿瘤大小、淋巴结转移情况、病理分期、分子分型、生存时间等信息。同时，收集本中心的乳腺癌患者样本，进行RNA提取和测序，获取miRNA表达数据，以增加数据的多样性和可靠性。数据预处理：运用生物信息学工具对下载和测序得到的miRNA表达数据进行预处理。首先，对数据进行标准化和归一化处理，消除不同实验条件和技术平台带来的差异，使数据具有可比性；然后，进行缺失值处理，对于缺失值较少的样本，采用均值、中位数等方法进行填充，对于缺失值较多的样本，考虑将其剔除。此外，对临床病理资料进行整理和清洗，确保数据的准确性和完整性。差异表达miRNA筛选：利用R语言的edgeR、limma等包，对预处理后的miRNA表达数据进行差异分析，筛选出在乳腺癌组织与正常组织中差异表达的miRNA。设定差异表达的阈值，如|logFC|>1且adj.P.Val<0.05，以确保筛选出的miRNA具有统计学意义。进一步对差异表达miRNA进行火山图、热图等可视化分析，直观展示差异表达情况。预后相关miRNA筛选：通过单因素Cox回归分析，评估差异表达miRNA与乳腺癌患者预后指标（如总生存期、无病生存期）之间的相关性，筛选出与预后相关的miRNA。为了避免过拟合，采用LASSO回归分析对与预后相关的miRNA进行进一步筛选和降维，确定最终的预后相关miRNA。预测模型构建：运用支持向量机（SVM）、随机森林（RF）、逻辑回归（LR）等机器学习算法，以筛选出的预后相关miRNA作为特征变量，构建乳腺癌患者预后预测模型。在Python的scikit-learn库中实现这些算法，通过交叉验证等方法对模型进行优化和调参，选择最优的模型参数。例如，对于SVM算法，调整核函数类型、惩罚参数C等；对于RF算法，调整决策树的数量、最大深度等参数。利用网格搜索、随机搜索等方法进行参数搜索，以找到最优的模型参数组合。模型验证与评估：采用内部验证和外部验证相结合的方式对构建的预测模型进行验证。内部验证使用同一数据集进行多次交叉验证，如10折交叉验证，评估模型在训练数据上的性能表现；外部验证则使用独立的数据集对模型进行测试，检验模型在不同数据上的泛化能力。通过计算准确率、召回率、F1值、受试者工作特征曲线（ROC）下面积等指标，全面评估模型的预测性能，判断模型的优劣。模型分析与应用：分析预测模型中miRNA与乳腺癌临床病理特征之间的关联，探讨miRNA在乳腺癌发生、发展及预后中的作用机制。将构建的预测模型应用于临床实践，为乳腺癌患者的预后评估提供参考，辅助临床医生制定个性化的治疗方案。[此处插入技术路线图]图1-1技术路线图二、相关理论基础2.1miRNA概述2.1.1miRNA的结构与功能miRNA是一类由内源基因编码的长度约为22个核苷酸的非编码单链RNA分子。其基因常以单拷贝、多拷贝或基因簇等多种形式存在于基因组中，具有固定的基因座位，部分位于蛋白基因的基因间隔区，其余在内含子，个别甚至在编码区的互补链。这表明它们的转录独立于其他基因，拥有自身独特的转录调控机制。miRNA的生成是一个复杂且精细的过程。首先，在细胞核中，miRNA基因在RNA聚合酶II的作用下转录形成初级转录本（pri-miRNA）。pri-miRNA长度可达几百到几千个核苷酸，具有复杂的二级结构，包含多个茎环结构。随后，pri-miRNA在一种被称为“微处理器复合体”（由Drosha酶和DGCR8蛋白组成）的作用下，被切割成约70-90个核苷酸的发夹状前体miRNA（pre-miRNA）。pre-miRNA在转运蛋白Exportin-5的协助下，通过消耗GTP，从细胞核转运至细胞质。在细胞质中，Dicer酶进一步识别并切割pre-miRNA，生成长度约为22个核苷酸的双链miRNA。双链miRNA中的一条链会与RNA诱导沉默复合体（RISC）结合，这条链被称为“引导链”，而另一条链（“乘客链”）通常会被降解，最终形成具有活性的成熟miRNA。成熟的miRNA主要通过与靶mRNA的3'非翻译区（3'UTR）进行碱基互补配对，从而对基因表达进行调控。当miRNA与靶mRNA的互补程度较高时，会促使靶mRNA降解；而当互补程度较低时，则主要抑制靶mRNA的翻译过程。例如，在线虫中，lin-4miRNA通过与lin-14mRNA的3'UTR不完全互补结合，抑制lin-14蛋白的翻译，从而调控线虫的发育进程。这种调控方式使得miRNA能够在转录后水平上精细地调节基因表达，参与细胞的增殖、分化、凋亡、代谢等多种生理过程。此外，一个miRNA可以调控多个靶基因，一个靶基因也可以受到多个miRNA的调控，这种复杂的调控网络使得miRNA在细胞内发挥着至关重要的作用。例如，在哺乳动物细胞中，miR-17-92基因簇编码的多个miRNA，它们协同作用，参与调控细胞的增殖、分化和凋亡等过程。miR-17-92基因簇中的miR-17和miR-20a可以通过靶向调控E2F1等基因，抑制细胞的增殖；而miR-92a则可以通过抑制PTEN基因的表达，促进细胞的存活和增殖。2.1.2miRNA与基因调控miRNA在基因转录后调控中扮演着核心角色，其调控机制主要通过与靶mRNA的相互作用来实现。当miRNA与靶mRNA的3'UTR区域碱基互补配对时，会形成miRNA-mRNA双链结构。这一结构的形成会招募相关的蛋白因子，从而引发不同的调控效应。在某些情况下，miRNA与靶mRNA的互补程度较高，几乎完全匹配，此时AGO蛋白（是RISC复合物的核心组成部分）会招募核酸内切酶，对靶mRNA进行切割，导致靶mRNA降解。例如，在植物中，miR-165/166可以与PHB等基因的mRNA完全互补配对，通过切割靶mRNA，调控植物的生长发育过程。研究表明，当miR-165/166表达上调时，PHB基因的mRNA水平显著下降，从而影响植物细胞的分化和组织形态建成。而当miRNA与靶mRNA的互补程度较低，存在错配碱基时，miRNA主要通过抑制靶mRNA的翻译过程来调控基因表达。具体来说，miRNA-mRNA双链结构会阻碍核糖体在mRNA上的移动，抑制蛋白质合成的起始或延伸步骤。例如，在动物细胞中，miR-122是肝脏中高度表达的一种miRNA，它可以与多个参与肝脏代谢和脂质合成的基因的mRNA3'UTR区域不完全互补结合，抑制这些基因的翻译，从而调控肝脏的生理功能。研究发现，当miR-122表达缺失时，相关基因的蛋白质表达水平升高，导致肝脏脂质代谢紊乱。除了上述经典的调控方式外，miRNA还可以通过其他机制影响基因表达。有研究发现，miRNA可以调控mRNA的稳定性，延长或缩短mRNA的半衰期。miRNA还可能参与染色质修饰等表观遗传调控过程，影响基因的转录活性。例如，在某些肿瘤细胞中，miRNA可以通过与DNA甲基转移酶等表观遗传调控因子相互作用，影响基因启动子区域的甲基化状态，进而调控基因的表达。此外，miRNA与靶基因之间的调控关系并非一成不变，而是受到多种因素的影响。细胞的生理状态、环境信号、其他非编码RNA等都可能干扰miRNA与靶mRNA的相互作用，从而改变基因的表达调控模式。例如，在细胞受到应激刺激时，一些miRNA的表达水平会发生变化，它们与靶mRNA的结合能力也会受到影响，进而引发细胞内基因表达谱的重塑，以适应应激环境。2.1.3miRNA在癌症中的作用机制miRNA在癌症的发生、发展过程中扮演着极为关键的角色，其作用机制复杂多样，既可以作为癌基因促进肿瘤的进展，也能够作为抑癌基因抑制肿瘤的生长。在癌基因方面，部分miRNA的异常高表达会对肿瘤的发生发展起到推动作用。以miR-21为例，大量研究表明，在乳腺癌、肺癌、胃癌等多种癌症中，miR-21均呈现高表达状态。其主要通过靶向作用于PTEN基因来发挥致癌作用。PTEN是一种重要的抑癌基因，它能够通过抑制PI3K/Akt信号通路，阻止细胞的异常增殖和存活。而miR-21与PTENmRNA的3'UTR区域互补配对，抑制PTEN的翻译过程，导致PTEN蛋白表达水平下降。PTEN的缺失使得PI3K/Akt信号通路被过度激活，促进细胞的增殖、迁移和侵袭，同时抑制细胞凋亡，最终导致肿瘤的发生和发展。有研究通过在乳腺癌细胞中抑制miR-21的表达，发现PTEN蛋白水平显著回升，PI3K/Akt信号通路活性受到抑制，癌细胞的增殖和侵袭能力明显减弱。另一个典型的致癌性miRNA是miR-155。在淋巴瘤、白血病以及多种实体瘤中，miR-155的表达显著上调。miR-155可以靶向调控多个基因，其中包括SOCS1等抑癌基因。SOCS1是细胞因子信号传导的负调节因子，能够抑制JAK/STAT信号通路的过度激活。miR-155通过抑制SOCS1的表达，导致JAK/STAT信号通路持续激活，促进肿瘤细胞的增殖、存活和免疫逃逸。研究还发现，miR-155在肿瘤微环境中也发挥着重要作用，它可以调节肿瘤相关巨噬细胞的功能，促进肿瘤的生长和转移。相反，一些miRNA作为抑癌基因，在癌症中表达下调，失去了对肿瘤细胞的抑制作用。miR-34家族是研究较为深入的抑癌性miRNA。miR-34a、miR-34b和miR-34c在多种癌症中表达缺失或低表达。它们主要通过靶向调控多个与细胞增殖、凋亡和周期调控相关的基因来发挥抑癌作用。miR-34家族可以靶向抑制SIRT1、CDK4、E2F3等基因的表达。SIRT1是一种去乙酰化酶，能够促进细胞的存活和增殖；CDK4和E2F3则参与细胞周期的调控。miR-34家族通过抑制这些基因的表达，诱导肿瘤细胞凋亡，阻滞细胞周期，从而抑制肿瘤的生长。在肺癌细胞中，过表达miR-34a可以显著抑制癌细胞的增殖能力，诱导细胞凋亡，并使细胞周期停滞在G1期。miR-125b也是一种重要的抑癌性miRNA。在乳腺癌、卵巢癌等多种癌症中，miR-125b的表达水平明显降低。miR-125b可以靶向作用于多个致癌基因，如ERBB2、MCL1等。ERBB2是一种受体酪氨酸激酶，其过表达与多种癌症的发生发展密切相关；MCL1是一种抗凋亡蛋白，能够抑制细胞凋亡。miR-125b通过抑制这些致癌基因的表达，抑制肿瘤细胞的增殖、迁移和侵袭，促进细胞凋亡。有研究在乳腺癌细胞中恢复miR-125b的表达，发现癌细胞的恶性行为受到明显抑制，肿瘤生长速度减缓。2.2乳腺癌相关知识2.2.1乳腺癌的发病机制乳腺癌的发病机制是一个多因素、多步骤的复杂过程，涉及遗传因素、激素水平失衡、生活方式及环境因素等多个方面，这些因素相互作用，共同影响着乳腺癌的发生发展。遗传因素在乳腺癌的发病中起着关键作用，约5%-10%的乳腺癌与遗传相关。其中，乳腺癌易感基因1（BRCA1）和乳腺癌易感基因2（BRCA2）是最为重要的遗传易感基因。BRCA1和BRCA2属于抑癌基因，其编码的蛋白质参与DNA损伤修复、细胞周期调控等重要生物学过程。当BRCA1或BRCA2基因发生突变时，会导致DNA损伤修复功能缺陷，使得细胞基因组的不稳定性增加，从而更容易积累致癌突变，进而增加乳腺癌的发病风险。研究表明，携带BRCA1或BRCA2基因突变的女性，其一生中患乳腺癌的风险可高达40%-80%。除BRCA1和BRCA2外，还有其他一些基因的突变也与乳腺癌的发病相关，如p53、PTEN、ATM等基因。p53基因是一种重要的肿瘤抑制基因，其突变会导致细胞凋亡受阻，细胞增殖失控，从而促进肿瘤的发生；PTEN基因编码的蛋白质具有磷酸酶活性，能够抑制PI3K/Akt信号通路，当PTEN基因发生突变或缺失时，PI3K/Akt信号通路被过度激活，促进乳腺癌细胞的生长和存活；ATM基因参与DNA损伤应答和细胞周期调控，ATM基因突变会影响细胞对DNA损伤的修复能力，增加乳腺癌的发病风险。激素水平失衡是乳腺癌发生的重要危险因素之一，雌激素和孕激素在乳腺癌的发病中起着关键作用。雌激素可以通过与雌激素受体（ER）结合，激活一系列信号通路，促进乳腺细胞的增殖和分化。长期高水平的雌激素刺激会导致乳腺细胞过度增殖，增加基因突变的概率，从而增加乳腺癌的发病风险。如初潮年龄早、绝经年龄晚、月经周期短等因素，都与乳腺癌的发病风险增加相关，这是因为这些因素会导致女性一生中暴露于雌激素的时间延长。此外，雌激素还可以通过诱导氧化应激、DNA损伤等机制，促进乳腺癌的发生发展。孕激素也与乳腺癌的发病有关，虽然孕激素本身对乳腺细胞具有一定的保护作用，但在某些情况下，孕激素可以协同雌激素促进乳腺细胞的增殖，尤其是在孕激素受体（PR）阳性的乳腺癌中，孕激素的作用更为明显。生活方式和环境因素对乳腺癌的发病也有着重要影响。长期高脂肪、高热量饮食会导致体重增加和肥胖，肥胖是乳腺癌的一个重要危险因素。肥胖会导致体内雌激素水平升高，脂肪组织可以将雄激素转化为雌激素，使得体内雌激素的合成增加；肥胖还会引起慢性炎症反应，炎症细胞分泌的细胞因子可以促进肿瘤细胞的生长和转移。缺乏运动也是乳腺癌的危险因素之一，适当的运动可以降低乳腺癌的发病风险。运动可以促进新陈代谢，调节激素水平，增强免疫力，减少炎症反应，从而对乳腺癌起到预防作用。此外，长期饮酒、吸烟、暴露于电离辐射等环境因素也与乳腺癌的发病相关。酒精可以干扰雌激素的代谢，增加乳腺癌的发病风险；吸烟会导致体内自由基增多，引起DNA损伤，促进肿瘤的发生；电离辐射可以直接损伤DNA，导致基因突变，增加乳腺癌的发病风险。2.2.2乳腺癌的临床病理特征乳腺癌的临床病理特征对于疾病的诊断、治疗和预后评估具有重要意义，主要包括分期、分型和常见病理特征等方面。乳腺癌的分期主要依据TNM分期系统，该系统通过评估原发肿瘤（T）、区域淋巴结（N）和远处转移（M）的情况来确定肿瘤的分期。原发肿瘤的大小和侵犯范围是分期的重要依据之一，T0表示原发癌瘤未查出，Tis指原位癌，包括导管原位癌和小叶原位癌，此时癌细胞尚未突破基底膜，属于早期阶段；T1肿瘤直径小于2cm，T2肿瘤直径为2-5cm，T3肿瘤直径大于5cm，T4则表示肿瘤侵及胸壁或皮肤，如炎性乳腺癌就属于T4分期，其病情较为严重。区域淋巴结的累及情况也对分期有重要影响，N0指同侧腋窝无淋巴结肿大，N1表示腋窝肿大淋巴结可推动，N2提示腋窝肿大淋巴结融合，N3表示有同侧胸骨旁和锁骨上淋巴结转移。远处转移情况是分期的关键指标，M0指无远处转移，M1为有远处转移。根据T、N、M的不同组合，乳腺癌可分为0期（TisN0M0）、I期（T1N0M0）、II期（T0-3N0-1M0）、III期（T0-3N1-2M0，或任何T4、任何N3）和IV期（有远处转移M1的任何情况）。分期越晚，患者的预后往往越差。乳腺癌的分型多样，根据病理组织学特点，可分为非浸润性癌和浸润性癌两大类。非浸润性癌包括导管原位癌、小叶原位癌和乳头湿疹样乳腺癌。导管原位癌是指癌细胞局限于乳腺导管内，未突破导管壁基底膜，其癌细胞沿着乳腺导管生长，可表现为粉刺型、非粉刺型等不同亚型，粉刺型导管原位癌的癌细胞核分级较高，坏死明显，预后相对较差；小叶原位癌是癌细胞未突破末梢乳管或腺泡基底膜，常表现为多中心性，双侧乳腺受累的可能性较大；乳头湿疹样乳腺癌较为少见，主要表现为乳头乳晕区的湿疹样改变，病理上可见Paget细胞，常伴有导管内癌。浸润性癌又可分为浸润性非特殊癌和浸润性特殊癌。浸润性非特殊癌是最常见的类型，占乳腺癌的70%-80%，包括浸润性导管癌、浸润性小叶癌、髓样癌、硬癌、单纯癌等。浸润性导管癌是临床上最常见的类型，癌细胞突破导管壁基底膜，向间质浸润生长，其组织学形态多样，预后相对较差；浸润性小叶癌的癌细胞呈单行串珠状或细条索状浸润于纤维间质中，常累及双侧乳腺，预后与浸润性导管癌相似。浸润性特殊癌相对少见，包括黏液腺癌、乳头状癌、腺样囊性癌等。黏液腺癌的癌细胞分泌大量黏液，形成黏液湖，癌细胞漂浮其中，预后相对较好；乳头状癌的癌细胞呈乳头状生长，预后也较好；腺样囊性癌的癌细胞呈筛状、腺样排列，恶性程度较低。乳腺癌的常见病理特征还包括肿瘤的组织学分级、分子分型等。组织学分级主要依据腺管形成的程度、细胞核的多形性和核分裂计数来评估。腺管形成程度高、细胞核多形性小、核分裂计数少的肿瘤，分级较低，提示肿瘤的分化程度高，预后相对较好；反之，分级较高，预后较差。分子分型是基于基因表达谱和免疫组织化学检测结果进行分类的，主要包括LuminalA型（ER+/PR+，HER-2-）、LuminalB型（ER+/PR+，HER-2+）、HER-2过表达型（ER-/PR-/HER-2+）和三阴性乳腺癌（ER-/PR-/HER-2-）。LuminalA型对内分泌治疗敏感，预后较好；LuminalB型对内分泌治疗也有一定效果，但HER-2的过表达提示预后相对较差；HER-2过表达型对HER-2靶向治疗敏感，但预后仍不理想；三阴性乳腺癌缺乏有效的治疗靶点，预后最差。2.2.3乳腺癌的治疗方法乳腺癌的治疗是一个综合性的过程，需要根据患者的病情、身体状况和分子分型等因素，制定个性化的治疗方案，以达到最佳的治疗效果。目前，主要的治疗方法包括手术治疗、化学治疗、放射治疗、靶向治疗和内分泌治疗等。手术治疗是乳腺癌的主要治疗手段之一，适用于早期和部分中期乳腺癌患者。手术方式主要包括乳腺癌根治术、改良根治术、保乳手术和乳房重建术。乳腺癌根治术切除范围广泛，包括整个乳房、胸大肌、胸小肌以及腋窝淋巴结，这种手术方式虽然能彻底切除肿瘤，但对患者的身体损伤较大，术后可能会出现上肢水肿、肩关节活动受限等并发症。改良根治术在保留胸大肌或胸大、小肌的基础上，切除乳房和腋窝淋巴结，与根治术相比，其对患者身体的损伤相对较小，术后上肢功能恢复较好，目前在临床上应用较为广泛。保乳手术则是在切除肿瘤的同时，尽可能保留乳房的外形和功能，适用于肿瘤较小、位置合适且患者有保乳意愿的情况。保乳手术需要严格掌握适应证，术后通常需要配合放射治疗，以降低局部复发的风险。乳房重建术是在乳腺癌手术切除乳房后，通过自体组织移植或假体植入等方法，重建乳房的外形，提高患者的生活质量。乳房重建术可以在乳腺癌手术的同时进行，也可以在术后一段时间进行，根据患者的具体情况和需求选择合适的时机。化学治疗是利用化学药物杀死癌细胞的治疗方法，适用于各期乳腺癌患者，尤其是晚期和复发转移的患者。化疗药物可以通过静脉注射、口服或局部注射等方式进入体内，作用于全身的癌细胞。化疗的目的主要包括术前新辅助化疗、术后辅助化疗和晚期姑息化疗。术前新辅助化疗是在手术前进行化疗，其作用是使肿瘤缩小，降低肿瘤分期，提高手术切除率，同时还可以观察肿瘤对化疗药物的敏感性，为术后化疗方案的选择提供依据。常用的新辅助化疗方案有蒽环类联合紫杉类，如表柔比星联合紫杉醇等。术后辅助化疗是在手术后进行化疗，旨在杀死残留的癌细胞，降低复发和转移的风险。术后辅助化疗的方案根据患者的病情和分子分型等因素选择，一般需要进行6-8个周期。晚期姑息化疗则是针对晚期或复发转移的乳腺癌患者，通过化疗来缓解症状，延长生存期。化疗药物在杀死癌细胞的同时，也会对正常细胞产生一定的损伤，导致一系列不良反应，如恶心、呕吐、脱发、骨髓抑制等。为了减轻这些不良反应，临床上通常会采取相应的对症治疗措施，如使用止吐药物缓解恶心呕吐，使用升白细胞药物治疗骨髓抑制等。放射治疗是利用放射线杀死癌细胞的局部治疗方法，主要用于乳腺癌术后辅助放疗和局部晚期乳腺癌的治疗。对于保乳手术的患者，术后放疗是必不可少的，它可以降低局部复发的风险，提高患者的生存率。放疗的范围通常包括患侧乳房、腋窝淋巴结和胸壁等部位。放疗的剂量和疗程根据患者的具体情况而定，一般需要进行5-6周的放疗。对于乳腺癌根治术或改良根治术的患者，如果存在高危因素，如腋窝淋巴结转移数较多、肿瘤较大、切缘阳性等，也需要进行术后放疗。此外，对于局部晚期乳腺癌患者，放疗可以与化疗、靶向治疗等联合应用，提高治疗效果。放射治疗也会产生一些不良反应，如皮肤损伤、放射性肺炎、上肢水肿等。在放疗过程中，医生会采取相应的措施来减轻这些不良反应，如使用皮肤保护剂预防皮肤损伤，控制放疗剂量和范围减少放射性肺炎的发生等。靶向治疗是针对肿瘤细胞的特定分子靶点进行治疗的方法，具有特异性强、疗效好、不良反应相对较小的特点。在乳腺癌中，HER-2是一个重要的治疗靶点。HER-2过表达的乳腺癌患者可以使用抗HER-2靶向药物进行治疗，如曲妥珠单抗、帕妥珠单抗等。曲妥珠单抗是第一个被批准用于治疗HER-2阳性乳腺癌的靶向药物，它可以与HER-2受体结合，抑制肿瘤细胞的增殖和存活，同时还可以激活机体的免疫细胞，发挥免疫杀伤作用。帕妥珠单抗则是与曲妥珠单抗联合使用，通过不同的作用机制，进一步增强对HER-2阳性乳腺癌的治疗效果。除了抗HER-2靶向药物外，近年来还出现了一些其他的靶向药物，如PI3K抑制剂、mTOR抑制剂等，这些药物针对乳腺癌细胞的不同信号通路进行靶向治疗，为乳腺癌患者提供了更多的治疗选择。靶向治疗的不良反应相对化疗较轻，但也可能会出现一些特殊的不良反应，如心脏毒性、腹泻等。在使用靶向药物治疗过程中，需要密切监测患者的不良反应，并及时进行处理。内分泌治疗是通过调节体内激素水平来抑制乳腺癌细胞生长的治疗方法，主要适用于激素受体阳性（ER+和/或PR+）的乳腺癌患者。内分泌治疗的药物主要包括选择性雌激素受体调节剂（SERM）、芳香化酶抑制剂（AI）和促性腺激素释放激素类似物（GnRHa）等。SERM如他莫昔芬，它可以与雌激素受体结合，阻断雌激素对乳腺癌细胞的刺激作用，从而抑制肿瘤细胞的生长。他莫昔芬适用于绝经前和绝经后的激素受体阳性乳腺癌患者，一般需要服用5-10年。AI如阿那曲唑、来曲唑和依西美坦等，它们通过抑制芳香化酶的活性，减少体内雌激素的合成，从而达到治疗乳腺癌的目的。AI主要适用于绝经后的激素受体阳性乳腺癌患者，其疗效优于他莫昔芬，不良反应相对较少。GnRHa如戈舍瑞林、亮丙瑞林等，主要用于绝经前激素受体阳性乳腺癌患者，通过抑制垂体分泌促性腺激素，降低体内雌激素水平，达到治疗的效果。内分泌治疗的不良反应相对较轻，常见的有潮热、盗汗、骨质疏松、阴道干涩等。为了预防和减轻这些不良反应，患者在治疗过程中可以采取相应的措施，如补充钙剂和维生素D预防骨质疏松，使用润滑剂缓解阴道干涩等。二、相关理论基础2.3数据分析与模型构建方法2.3.1数据处理工具与软件在本研究中，使用了多种专业的数据处理工具与软件，以确保数据的高效处理和分析的准确性。RStudio作为一款集成开发环境，为R语言编程提供了便利的操作界面。R语言拥有丰富的包资源，其中edgeR包主要用于RNA-seq数据的差异表达分析。它基于负二项分布模型，能够有效地处理计数数据，准确地识别出在不同样本组之间差异表达的miRNA。通过edgeR包，可以对乳腺癌患者和正常对照组的miRNA表达数据进行分析，筛选出具有显著表达差异的miRNA，为后续研究提供关键的候选分子。caret包在机器学习模型的构建和评估中发挥了重要作用。该包提供了一系列函数和工具，用于数据预处理、模型训练、参数调优以及模型评估等环节。在构建乳腺癌患者预后预测模型时，利用caret包可以方便地对数据进行划分，将数据集分为训练集和测试集，确保模型的泛化能力。它还支持多种机器学习算法的实现，如支持向量机、随机森林等，通过统一的接口进行模型训练和评估，大大提高了研究效率。在使用支持向量机算法构建模型时，caret包可以帮助我们快速地调整模型参数，如核函数类型、惩罚参数等，以找到最优的模型配置，提高模型的预测性能。GraphPadPrism9是一款功能强大的数据分析和绘图软件，常用于生物医学研究领域。在本研究中，利用GraphPadPrism9对数据进行可视化展示，将复杂的数据以直观的图表形式呈现出来。通过绘制柱状图，可以清晰地比较不同组之间miRNA表达水平的差异；箱线图则能够展示数据的分布情况，帮助我们快速了解数据的集中趋势和离散程度；散点图可以用于分析miRNA表达水平与其他临床病理指标之间的相关性。这些可视化图表不仅有助于我们更好地理解数据，还为研究结果的展示和论文撰写提供了有力支持。在分析miRNA表达水平与乳腺癌患者年龄的关系时，使用散点图可以直观地观察到两者之间是否存在线性或非线性关系，为进一步的统计分析提供了线索。2.3.2差异表达分析差异表达分析是筛选与乳腺癌相关miRNA的关键步骤，通过该分析可以确定在乳腺癌组织和正常组织中表达水平存在显著差异的miRNA，这些差异表达的miRNA可能在乳腺癌的发生、发展过程中发挥重要作用。本研究采用R语言的edgeR包进行miRNA的差异表达分析。首先，对原始的miRNA表达数据进行标准化处理，消除不同样本之间由于测序深度、实验批次等因素造成的差异，使数据具有可比性。在RNA-seq实验中，不同样本的测序深度可能不同，这会导致基因表达量的估计出现偏差。通过标准化处理，可以将不同样本的表达数据调整到同一尺度上，从而更准确地进行差异表达分析。在标准化处理后，使用edgeR包中的精确检验方法对两组样本（乳腺癌组织和正常组织）的miRNA表达数据进行统计检验。该方法基于负二项分布模型，考虑了基因表达的生物学变异性和技术重复性，能够更准确地检测出差异表达的miRNA。在计算过程中，会得到每个miRNA的差异表达倍数（logFC）和相应的P值。差异表达倍数表示乳腺癌组织中miRNA的表达水平相对于正常组织的变化倍数，P值则用于衡量这种差异的统计学显著性。为了控制假阳性率，采用Benjamini-Hochberg方法对P值进行校正，得到校正后的P值（adj.P.Val）。设定差异表达的阈值为|logFC|>1且adj.P.Val<0.05。|logFC|>1表示miRNA在两组样本中的表达差异达到2倍及以上，具有较为明显的表达变化；adj.P.Val<0.05则保证了这种差异具有统计学意义，即不是由于随机因素造成的。通过这两个阈值的筛选，可以得到在乳腺癌组织和正常组织中显著差异表达的miRNA列表。进一步对筛选出的差异表达miRNA进行火山图和热图可视化分析。火山图以logFC为横坐标，-log10(adj.P.Val)为纵坐标，将每个miRNA在图中用一个点表示。图中的点根据其logFC和adj.P.Val的值分布在不同位置，差异表达显著的miRNA（|logFC|>1且adj.P.Val<0.05）会分布在图的两侧，远离中心线，从而直观地展示出差异表达miRNA的分布情况。热图则通过颜色的深浅来表示miRNA在不同样本中的表达水平，将差异表达miRNA在乳腺癌组织和正常组织中的表达情况以矩阵的形式呈现出来。通过热图，可以清晰地看到不同miRNA在两组样本中的表达模式，以及样本之间的表达差异，有助于发现潜在的miRNA表达特征和规律。2.3.3回归分析回归分析是本研究中用于筛选与乳腺癌患者预后相关miRNA，并构建预后预测模型的重要统计方法，主要包括单因素Cox回归分析、LASSO回归分析和多因素Cox回归分析。单因素Cox回归分析用于评估每个miRNA与乳腺癌患者预后（如总生存期、无病生存期）之间的相关性。Cox回归模型是一种半参数模型，它不需要对生存时间的分布做出假设，能够同时考虑多个因素对生存时间的影响。在单因素Cox回归分析中，将每个miRNA的表达水平作为自变量，患者的生存时间和生存状态作为因变量，通过计算风险比（HR）和95%置信区间（CI）来评估miRNA对预后的影响程度。如果某个miRNA的HR>1且95%CI不包含1，说明该miRNA的高表达与患者的不良预后相关，即高表达该miRNA的患者死亡风险更高；反之，如果HR<1且95%CI不包含1，则说明该miRNA的高表达与患者的良好预后相关，高表达该miRNA的患者死亡风险较低。通过单因素Cox回归分析，可以初步筛选出与乳腺癌患者预后相关的miRNA，为后续的分析提供候选变量。然而，单因素Cox回归分析可能会受到多重共线性的影响，导致结果不准确。为了克服这一问题，采用LASSO（LeastAbsoluteShrinkageandSelectionOperator）回归分析对与预后相关的miRNA进行进一步筛选和降维。LASSO回归是一种压缩估计方法，它在回归系数的估计过程中加入了L1正则化项，能够使一些不重要的变量的系数变为0，从而实现变量选择和降维的目的。在LASSO回归分析中，通过交叉验证的方法选择最优的正则化参数lambda，以确保模型的准确性和稳定性。经过LASSO回归分析，可以从单因素Cox回归筛选出的miRNA中进一步筛选出对预后影响最为显著的miRNA，减少模型中的变量数量，提高模型的可解释性和预测性能。多因素Cox回归分析则是将LASSO回归筛选出的miRNA纳入模型中，同时考虑其他临床病理因素（如年龄、肿瘤大小、淋巴结转移情况等），构建多因素Cox回归模型。通过多因素Cox回归分析，可以确定这些miRNA是否为乳腺癌患者预后的独立危险因素，并计算出每个因素的HR和95%CI。多因素Cox回归模型能够综合考虑多个因素对预后的影响，更准确地评估患者的预后风险。在多因素Cox回归模型中，如果某个miRNA的HR和95%CI在调整其他因素后仍然具有统计学意义，说明该miRNA是乳腺癌患者预后的独立危险因素，对患者的预后具有重要的预测价值。通过多因素Cox回归分析构建的模型，可以为临床医生提供更全面、准确的预后评估工具，帮助他们制定个性化的治疗方案。2.3.4模型评估方法模型评估是构建乳腺癌患者预后预测模型过程中的关键环节，通过多种评估方法可以全面、客观地评价模型的性能，确保模型的可靠性和有效性，为临床应用提供有力支持。受试者工作特征曲线（ROC）是评估模型预测性能的常用方法之一。ROC曲线以真阳性率（sensitivity）为纵坐标，假阳性率（1-specificity）为横坐标，通过绘制不同阈值下模型预测结果的真阳性率和假阳性率，展示模型在不同分类阈值下的性能表现。在乳腺癌患者预后预测模型中，将模型预测的患者生存状态（生存或死亡）与实际生存状态进行比较，计算不同阈值下的真阳性率和假阳性率，从而绘制出ROC曲线。曲线下面积（AUC）是衡量ROC曲线性能的重要指标，AUC取值范围在0.5到1之间，AUC越接近1，表明模型的预测准确性越高；当AUC=0.5时，说明模型的预测效果与随机猜测无异。在本研究中，通过计算AUC来评估构建的乳腺癌患者预后预测模型的准确性，AUC值越高，说明模型能够更好地区分高风险和低风险患者，为临床决策提供更有价值的信息。除了ROC曲线，还使用准确率（Accuracy）、召回率（Recall）和F1值等指标来评估模型的性能。准确率是指模型预测正确的样本数占总样本数的比例，反映了模型整体的预测准确性。在乳腺癌患者预后预测模型中，准确率=（真阳性样本数+真阴性样本数）/总样本数。召回率，也称为灵敏度或真阳性率，是指实际为阳性的样本中被正确预测为阳性的比例，体现了模型对正样本的识别能力。在本研究中，召回率=真阳性样本数/（真阳性样本数+假阴性样本数）。F1值是综合考虑准确率和召回率的指标，它是准确率和召回率的调和平均数，能够更全面地反映模型的性能。F1值的计算公式为：F1=2×（准确率×召回率）/（准确率+召回率）。通过这些指标的计算，可以从不同角度评估模型的性能，了解模型在预测乳腺癌患者预后时的优势和不足。此外，采用交叉验证的方法对模型进行内部验证，以提高模型评估的可靠性。常见的交叉验证方法有K折交叉验证，如10折交叉验证。在10折交叉验证中，将数据集随机分为10个大小相近的子集，每次选取其中9个子集作为训练集，用于训练模型，剩下的1个子集作为测试集，用于评估模型的性能。重复这个过程10次，每次使用不同的子集作为测试集，最后将10次的评估结果进行平均，得到模型的性能指标。通过交叉验证，可以充分利用数据集的信息，减少因数据集划分方式不同而导致的模型性能波动，更准确地评估模型的泛化能力。在本研究中，通过10折交叉验证对构建的乳腺癌患者预后预测模型进行内部验证，确保模型在不同的数据子集上都具有较好的性能表现，提高了模型的可靠性和稳定性。三、乳腺癌患者数据收集与预处理3.1数据来源3.1.1TCGA数据库本研究主要从癌症基因组图谱（TheCancerGenomeAtlas，TCGA）数据库中获取乳腺癌患者的miRNA测序数据和临床资料。TCGA是一个大规模的癌症基因组学研究项目，它整合了来自全球多个研究机构的癌症数据，包括肿瘤组织和正常组织的基因组、转录组、蛋白质组等多组学数据，以及详细的临床病理信息。该数据库涵盖了多种癌症类型，其中乳腺癌数据丰富，包含了大量的样本和全面的信息，为研究乳腺癌的发病机制、诊断和预后提供了宝贵的资源。在数据获取过程中，首先通过R语言的TCGAbiolinks包与TCGA数据库进行交互。该包提供了一系列函数，方便用户从TCGA数据库中查询、下载和处理数据。使用getGDCprojects函数获取TCGA中所有项目的信息，从中筛选出乳腺癌项目（项目ID为“TCGA-BRCA”）。然后，利用GDCquery函数构建查询对象，设置参数data.category为“TranscriptomeProfiling”，表示获取转录组数据；data.type为“miRNAExpressionQuantification”，明确获取的是miRNA表达定量数据；workflow.type为“BCGSCmiRNAProfiling”，指定数据处理流程类型。通过这样的设置，确保获取到的miRNA测序数据是经过特定标准化处理流程的，以保证数据的质量和一致性。执行GDCquery函数后，得到符合条件的数据集信息。接着，使用GDCdownload函数下载数据，将数据保存到本地指定目录。下载完成后，利用GDCprepare函数对下载的数据进行预处理，将原始数据转换为R语言中便于分析的数据格式，并保存为Rdata文件。通过这一系列操作，成功从TCGA数据库中获取了乳腺癌患者的miRNA测序数据。在获取临床资料方面，同样使用TCGAbiolinks包。通过设置相应的参数，查询并下载乳腺癌患者的临床信息，包括患者的年龄、性别、肿瘤大小、淋巴结转移情况、病理分期、分子分型、生存时间等。这些临床资料与miRNA测序数据相对应，为后续分析miRNA表达与临床病理特征之间的关系，以及构建预后预测模型提供了重要的基础数据。3.1.2GEO数据库除了TCGA数据库，还从基因表达综合数据库（GeneExpressionOmnibus，GEO）中下载乳腺癌患者的mRNA测序数据。GEO是一个公共的功能基因组学数据存储库，由美国国立生物技术信息中心（NCBI）维护。它收录了来自世界各地的高通量基因表达实验数据，包括芯片数据、测序数据等，涵盖了多种生物样本和研究领域。GEO数据库的数据来源广泛，数据类型丰富，为生物医学研究提供了大量的公开数据资源。从GEO数据库下载mRNA测序数据的过程如下：首先，在GEO数据库的官方网站（/geo/）上进行检索。在搜索框中输入关键词“breastcancer”和“mRNAsequencing”，并设置其他筛选条件，如物种为“human”，样本类型为“tumor”和“normal”，以缩小搜索范围，获取与乳腺癌相关的mRNA测序数据集。经过筛选，选择了符合研究需求的数据集，如GSE122768等。这些数据集包含了乳腺癌患者肿瘤组织和正常组织的mRNA测序数据，以及相应的样本信息。对于选定的数据集，点击进入数据集详情页面。在该页面中，可以查看数据集的详细描述、样本信息、实验方法等内容。如果数据集提供了原始测序数据的下载链接，可直接点击链接进行下载；若数据集只提供了处理后的表达矩阵文件，则下载该文件。在下载过程中，根据数据的大小和网络状况，可能需要花费一定的时间。下载完成后，将数据保存到本地指定目录。对于下载的mRNA测序数据，若为原始测序数据（通常为fastq格式文件），需要进行一系列的预处理步骤，包括质量控制、序列比对、基因表达定量等。使用FastQC软件对原始测序数据进行质量评估，检查数据的质量分布、碱基组成、GC含量等指标，判断数据是否存在质量问题。如果数据质量不佳，可使用Trimmomatic等软件进行数据清洗，去除低质量的碱基和接头序列。然后，利用Hisat2等比对软件将清洗后的数据比对到人类参考基因组上，生成比对结果文件（如bam格式文件）。最后，使用StringTie等软件对bam文件进行基因表达定量分析，计算每个基因的表达量，得到mRNA表达矩阵。若下载的是已经处理好的表达矩阵文件，则可直接进行后续的数据分析。三、乳腺癌患者数据收集与预处理3.2数据筛选与清洗3.2.1数据质量控制在从TCGA和GEO数据库获取乳腺癌患者的miRNA测序数据、mRNA测序数据以及临床资料后，首要任务是进行严格的数据质量控制，以确保后续分析结果的可靠性和准确性。数据质量控制主要聚焦于去除低质量数据和异常值，采用了一系列严谨的方法和步骤。对于测序数据，首先运用FastQC软件进行全面的质量评估。FastQC能够生成详细的质量报告，从多个维度展示数据质量信息。在碱基质量分布方面，它会分析每个位置碱基的质量得分情况，若发现某些位置的碱基质量得分普遍较低，可能是测序过程中出现了问题，这些数据需要进一步处理或舍弃。数据的GC含量也是重要的评估指标，正常情况下，测序数据的GC含量应在一定范围内波动，如果GC含量过高或过低，都可能暗示数据存在异常，如样本污染、文库构建问题等。通过FastQC的评估，能够快速发现数据中潜在的质量问题，为后续的数据处理提供依据。在处理低质量数据时，若发现某样本的低质量碱基比例超过一定阈值（如20%），则考虑对该样本进行数据清洗。使用Trimmomatic软件进行清洗操作，它可以根据设定的参数，去除低质量的碱基和接头序列。在设置参数时，通常将碱基质量阈值设为20，即当碱基质量得分低于20时，认为该碱基质量较低，需要被去除；对于接头序列，Trimmomatic能够自动识别并去除常见的接头，确保数据的纯净度。如果经过清洗后，样本的数据量仍然过少，无法满足后续分析的要求，如有效reads数低于100万，则将该样本剔除。异常值的检测和处理同样至关重要。采用箱线图（BoxPlot）方法来识别异常值。箱线图能够直观地展示数据的分布情况，通过绘制miRNA或mRNA表达量的箱线图，可以清晰地看到数据的四分位数、中位数以及异常值的分布。在箱线图中，若某个数据点距离上四分位数或下四分位数超过1.5倍的四分位间距（IQR），则将其判定为异常值。以miRNA表达量数据为例，假设某样本的miRNA表达量在箱线图中明显偏离其他样本，超出了1.5倍IQR的范围，那么这个样本就可能是异常值。对于判定为异常值的样本，进一步检查其数据来源、实验操作记录等信息，以确定异常的原因。如果是由于实验误差、样本标记错误等可纠正的原因导致的异常值，对其进行修正；若无法确定原因或认为是真实的异常情况，但该异常值对整体分析结果可能产生较大影响时，将其从数据集中剔除。除了上述基于统计方法的异常值检测，还结合领域知识进行判断。在乳腺癌研究领域，已知某些临床特征与miRNA或mRNA表达量之间存在一定的关联。如果某个样本的miRNA表达量与该样本的临床特征严重不符，如根据临床诊断为早期乳腺癌的患者，其某些与肿瘤进展相关的miRNA表达量却异常高，且与大量其他早期乳腺癌患者的数据表现差异显著，那么这个样本也可能被视为异常值，需要进一步分析和处理。3.2.2数据标准化数据标准化是数据预处理过程中的关键环节，它能够消除不同样本之间由于实验条件、测序深度等因素造成的差异，使数据具有可比性，为后续的数据分析和模型构建提供可靠的基础。本研究采用了多种数据标准化方法，对miRNA测序数据、mRNA测序数据以及临床资料进行了细致的标准化处理。对于miRNA测序数据，由于其数据类型为计数数据，采用了TPM（TranscriptsPerMillion）标准化方法。TPM标准化考虑了测序深度和基因长度对表达量的影响，能够更准确地反映miRNA的真实表达水平。具体计算步骤如下：首先，对于每个miRNA，将其原始计数（count）除以该样本中所有miRNA的总计数，得到每个miRNA在该样本中的相对表达比例；然后，将这个相对表达比例乘以1000000，得到以每百万转录本为单位的表达量，即TPM值。假设某样本中miRNAA的原始计数为100，该样本中所有miRNA的总计数为10000000，则miRNAA的相对表达比例为100/10000000=0.00001，其TPM值为0.00001×1000000=10。通过TPM标准化，不同样本之间的miRNA表达量可以在同一尺度上进行比较，避免了由于测序深度不同而导致的表达量偏差。对于mRNA测序数据，使用了DESeq2包中的标准化方法。DESeq2基于负二项分布模型，能够有效处理RNA-seq数据中的技术变异和生物学变异。在标准化过程中，DESeq2首先对原始计数数据进行归一化处理，通过估计每个样本的大小因子（sizefactor）来校正测序深度的差异。大小因子反映了每个样本相对于一个参考样本的测序深度倍数，通过将每个样本的原始计数除以其对应的大小因子，实现了对测序深度的标准化。DESeq2还考虑了基因长度对表达量的影响，通过对基因长度进行校正，使不同长度的基因表达量具有可比性。在处理过程中，DESeq2会自动识别并处理数据中的异常值和离群点，进一步提高了数据的质量和可靠性。临床资料中的数据类型较为复杂，包括数值型数据（如年龄、肿瘤大小等）和分类数据（如病理分期、分子分型等），因此采用了不同的标准化方法。对于数值型数据，使用Z-score标准化方法。Z-score标准化是基于原始数据的均值（mean）和标准差（standarddeviation）进行的，其计算公式为：新数据=（原数据-均值）/标准差。以年龄数据为例，假设所有患者年龄的均值为50岁，标准差为10岁，某患者的年龄为60岁，则经过Z-score标准化后，该患者的年龄数据变为（60-50）/10=1。经过Z-score标准化后，数据的均值变为0，标准差变为1，不同数值型变量之间的量纲差异被消除，便于后续的数据分析和模型构建。对于分类数据，采用独热编码（One-HotEncoding）方法进行处理。独热编码是将每个类别变量转换为一个二进制向量，向量中只有一个元素为1，其余元素为0，以此来表示不同的类别。在处理病理分期数据时，假设病理分期分为I期、II期、III期和IV期。使用独热编码后，I期可以表示为[1,0,0,0]，II期表示为[0,1,0,0]，III期表示为[0,0,1,0]，IV期表示为[0,0,0,1]。通过独热编码，分类数据被转换为数值型向量，能够更好地与其他数值型数据进行整合和分析，同时也便于机器学习算法对分类数据的处理。3.3数据分组3.3.1训练组与验证组划分在构建乳腺癌患者预后预测模型时，为了确保模型具有良好的泛化能力，需要对数据进行合理分组，划分为训练组和验证组。本研究使用R语言中的caret包来完成这一任务。caret包提供了createDataPartition函数，能够根据指定的比例将数据集随机划分为训练组和验证组。首先，加载caret包。使用library(caret)命令，确保R语言环境中已经安装了caret包，若未安装，则需要先使用install.packages("caret")命令进行安装。在加载包之后，读取经过预处理的乳腺癌患者数据集，假设数据集存储在名为data的对象中。接下来，使用createDataPartition函数进行数据分组。在函数中，第一个参数为要划分的目标变量，通常选择与患者预后密切相关的变量，如生存状态（survival_status）。本研究将生存状态作为目标变量，即data$survival_status。第二个参数p表示训练组所占的比例，经过综合考虑和多次试验，本研究设定p=0.8，即训练组占数据集的80%，验证组占20%。list参数设置为FALSE，表示返回的结果是一个向量，而不是列表。times参数表示重复划分的次数，这里设置为1，即只进行一次划分。具体代码如下：set.seed(123)#设置随机数种子，确保结果可重复trainIndex<-createDataPartition(data$survival_status,p=0.8,list=FALSE,times=1)trainData<-data[trainIndex,]#训练组数据testData<-data[-trainIndex,]#验证组数据通过上述代码，首先使用set.seed(123)设置随机数种子为123，这一步非常关键，它可以确保每次运行代码时，数据划分的结果是一致的，便于研究结果的复现。然后，createDataPartition函数根据生存状态变量将数据集划分为训练组和验证组，并返回训练组的索引。利用这些索引，从原始数据集中提取出训练组数据trainData和验证组数据testData。3.3.2分组合理性验证数据分组的合理性对模型的构建和验证具有重要影响。合理的分组能够使训练组充分代表总体数据的特征，同时验证组能够有效地评估模型在未知数据上的性能。若分组不合理，可能导致训练组和验证组的数据分布差异较大，从而使模型在训练组上表现良好，但在验证组上性能急剧下降，即出现过拟合现象。为了验证分组的合理性，本研究从多个方面进行了分析。首先，对训练组和验证组的样本数量进行了比较。通过计算可知，训练组样本数量为[X1]，验证组样本数量为[X2]，两者比例接近设定的80%和20%，表明样本数量的划分较为合理。这种合理的样本数量分配，能够保证训练组有足够的数据用于模型训练，使模型能够学习到数据中的各种特征和规律；同时，验证组也有足够的样本用于评估模型的性能，避免因样本数量过少而导致评估结果不准确。接着，分析了两组样本中乳腺癌患者的临床病理特征分布情况。以肿瘤大小为例，在训练组中，肿瘤大小的均值为[mean1]，标准差为[sd1]；在验证组中，肿瘤大小的均值为[mean2]，标准差为[sd2]。通