冰山之下的暗流：血浆蛋白质组的全面解码
长期以来，关于生理周期的分子研究大多依赖于靶向的激素检测、小规模的免疫介质分析，或是对经血和子宫内膜组织的局部取样。这些方法虽然有其价值，但在样本量和时间分辨率上往往捉襟见肘。血浆蛋白质作为连接各个组织和生物学通路的“信使”，提供了一个独特且动态的窗口，能够以宏观的视角读取局部与全身的生理状态。
为了绘制这张宏观地图，研究人员调用了英国生物样本库（UK Biobank）这一庞大资源，对2760名健康女性的血浆样本进行了无靶向的全面蛋白质组学分析。借助于先进的邻近延伸分析技术，研究共对2917种蛋白质进行了精准定量。在这片浩瀚的蛋白质海洋中，研究人员通过严格的统计学分析（假发现率控制在0.05以下），最终锁定了198种与生理周期天数显著相关的蛋白质。
在这198种蛋白质中，关联性最强的前五名分别是PROK1、MMP10、RLN1、CXCL13和CGA。其中，CGA是卵泡刺激素、黄体生成素和人绒毛膜促性腺激素共有的α亚基，而FSHB则是卵泡刺激素的β亚基。这两种经典生殖内分泌标志物的出现，宛如航海图上的核心坐标，不仅验证了该分析系统的极高可靠性，也为探索其他未知蛋白质的节律奠定了坚实的基础。
更有意思的是，当我们退一步审视整个血浆蛋白质组时，会发现生理周期带来的改变呈现出一种高度的“克制”与“精准”。数据表明，生理周期天数能够解释其变异程度超过5%的蛋白质，仅仅只有9种（占比仅为0.31%）。作为对比，年龄和身体质量指数（BMI）则分别对391种（占比13.4%）和332种（占比11.3%）蛋白质的变异有着超过5%的解释力。
这一强烈的数据对比揭示了一个深刻的生物学事实：生理周期对人体的调控并非一场狂风暴雨式的全方位颠覆，而是类似于一场极其精准的外科手术。它仅对特定的生理网络施加深远影响，而年龄和体型等因素则在更为广阔的蛋白质谱系中留下岁月的痕迹。
时间的具象化：四重乐章般的蛋白质表达轨迹
为了探寻这198种蛋白质在28天周期内的动态表达规律，研究人员运用了K-均值聚类算法（K-means clustering），将这些复杂的升降起伏梳理成了四个特征鲜明的时间表达簇。这就如同将杂乱无章的音符，谱写成了四个相互交织、此起彼伏的乐章，完美地呼应了子宫内膜和卵巢在不同阶段的生理状态。
第一个乐章在月经的第一天奏响，达到了表达的巅峰。这个由53种蛋白质组成的簇群，包含了MMP10、PAEP和LEFTY2等关键分子。通路富集分析显示，这个时期的蛋白质高度集中在“蛋白酶抑制剂活性”通路上。我们可以想象，当月经来潮时，子宫内膜面临着大规模的脱落与组织重构，这绝不是简单的剥落，而是一个需要被极其严格控制的组织降解过程。蛋白酶及其抑制剂在这里扮演了清道夫与交通警察的双重角色，确保旧组织的清除既高效又不会对深层基底造成破坏。
第二个乐章随着周期的推进，在卵泡期早期缓缓展开。这个簇群包含了69种蛋白质，以经典的激素标志物CGA和FSHB，以及TNFSF10为代表。此时，机体正处于为卵泡发育和内膜再生蓄力的阶段。数据指出，这些蛋白质在“细胞因子活性”、“激素活性”和“胰岛素样生长因子结合”通路中显著富集。免疫系统和内分泌系统在此刻紧密合作，通过炎症信号的微调和生长因子的释放，悄然拉开了新一轮生命周期的帷幕。
第三个乐章当时间轴推进到卵泡期晚期和排卵期前后时迎来了高潮。这20种蛋白质（如STC2、SFRP4和CRLF1）呈现出明显的峰值。令人惊叹的是，在这个被生殖生长通路主导的阶段，除了GH1、PRL等激素外，催产素（OXT）也展现出了显著的周期性波动，并在排卵期前后达到顶峰。这种设计的巧妙之处在于，催产素的飙升不仅可能增强女性在排卵期的性唤起，更有助于促进生殖道的收缩与蠕动，从而为迎接可能的受精过程做好充分的物理和生理准备。
最后的第四个乐章在黄体期中晚期深沉作响。包含56种蛋白质的这个簇群，深深根植于“细胞因子信号传导”、“模式识别受体活性”以及MAPK通路的调节中。此时，在孕酮的驱动下，子宫内膜正在经历复杂的基质重塑和成熟过程。免疫细胞在内膜中穿梭巡逻，构建起一道既能容受胚胎着床，又能防御外界病原体的微环境屏障。值得注意的是，作为肾素-血管紧张素-醛固酮系统（RAAS）的限速酶，肾素（REN）也在黄体晚期出现峰值。这暗示了局部产生的肾素在调节子宫内膜血流、乃至触发月经启动中不可或缺的角色。
寻找源头：从微观单细胞到宏观全身循环
既然我们在外周血中捕捉到了这些随周期起伏的蛋白质，那么它们的源头究竟在哪里？为了回答这个问题，研究人员首先利用人类蛋白质图谱（Human Protein Atlas）中40种不同组织的表达数据进行了富集分析。结果毫无悬念，这些周期相关蛋白质在子宫内膜和胎盘组织中表现出了极其显著的特异性富集。
但这还不够深入。考虑到子宫内膜复杂的细胞构成，研究人员引入了来自10名健康供体的单细胞RNA测序（scRNA-seq）数据。通过对多达71032个单细胞的转录组进行剥茧抽丝般的解析，他们在细胞微观层面上为这些蛋白质找到了确切的“生产车间”。在112种能被定位的蛋白质中，研究人员锁定了184个显著的基因表达关联。这些蛋白质主要由三种细胞大量“制造”：基质成纤维细胞（37种）、无纤毛上皮细胞（36种）以及有纤毛上皮细胞（23种）。
这种微观定位为我们理解宏观生理现象提供了直接的证据支撑。以PAEP（妊娠相关子宫内膜蛋白）为例，单细胞数据清晰地显示，它主要由无纤毛上皮细胞高度特异性地表达。PAEP是一种受孕酮严格调控的分泌蛋白，对于构建子宫内膜容受性（即内膜接纳受精卵着床的能力）至关重要。另一个极具代表性的分子是CXCL8（白细胞介素-8），作为一种关键的趋化因子，它的主要职责是招募中性粒细胞，推动细胞外基质的重塑和血管生成的信号传导。在月经来潮前后，CXCL8水平的显著攀升，正是子宫内膜进行自我更新和血管网重建的微观缩影。
通过这样的跨尺度分析，我们清晰地看到了一条从局部到全身的信息传递链条。子宫内膜的各种特定细胞在激素的指令下，有条不紊地合成并分泌特定的蛋白质。这些蛋白质透过血管壁进入全身血液循环，不仅在局部执行着组织重塑和免疫调节的功能，更将生殖系统的状态如实地广播给全身各个脏器。
疾病的伏笔：当生理节律偏离正常轨道
生理状态下的波动，往往是理解病理状态的一面镜子。这些在健康女性体内按部就班波动的蛋白质，是否也参与了那些困扰无数女性的生殖系统疾病的发生发展？
为了探究这一问题，研究人员将视野扩大到了拥有完整Olink蛋白质组学数据的23674名女性队列中，并对42种常见女性生殖系统疾病进行了全表型关联研究（PheWAS）。数据呈现出的结果令人十分瞩目。
在疾病的患病率（即基线已有的疾病）分析中，研究发现了60对显著的蛋白质-疾病关联：
• 蛋白质 CGA 与月经过多、频繁和不规则出血的关联最为紧密，其蛋白质水平每增加一个标准差，相关疾病的几率比达到 OR = 1.35
• 对于卵巢癌，蛋白质 ITIH4 水平的升高带来了 OR = 1.50，而 IGDCC4 水平的升高则显示出保护性趋势，其 OR = 0.51
在发病率（即随访期间新发的疾病）分析中，线索变得更加引人深思：
• 蛋白质 CHRDL2 水平的升高与子宫平滑肌瘤的新发风险密切相关，高达 OR = 2.05
• 由无纤毛上皮细胞分泌的 PAEP，其水平的增加更是与未来的月经过多及不规则出血呈现出极强的相关性，达到了惊人的 OR = 2.44
这些强有力的数据表明，生理周期相关蛋白质并非仅仅是生理起伏的被动旁观者或标记物。虽然子宫内膜异位症、子宫平滑肌瘤和异常子宫出血在临床表现上大相径庭，但它们在深层生物学机制上却共享着炎症、血管生成以及类固醇激素响应的复杂通路。当这些本该精准起伏的蛋白质偏离了它们应有的节律轨道，或许正是机体在为未来的组织病变埋下伏笔。
探寻因果：是谁在暗中扣动了疾病的扳机？
为了搞清楚这些关联究竟是疾病导致的结果，还是引发疾病的真正元凶，研究人员动用了孟德尔随机化（Mendelian randomization, MR）。
孟德尔随机化的巧妙在于，它利用了我们在受精那一刻就随机分配并终生不变的遗传变异（这里使用的是蛋白质定量性状位点，cis-pQTLs）作为工具变量。因为基因变异不受后天环境和疾病反向因果的影响，这就相当于在人群中进行了一场天然的随机对照试验。
将这些遗传工具变量与涵盖42种女性生殖疾病的大规模全基因组关联分析（GWAS）数据进行碰撞后，研究人员最终锁定了14对具有潜在因果关系的蛋白质-疾病组合。在这其中，FSHB（卵泡刺激素β亚基）的表现尤为引人侧目。
在传统的观念中，卵泡刺激素（FSH）主要被视为作用于卵巢的激素。然而，孟德尔随机化分析提供了一组颠覆性的数据：FSHB水平的遗传性升高，能够直接增加子宫内膜异位症、异常子宫出血以及子宫平滑肌瘤的患病风险。事实上，组织学研究已经发现，子宫内膜的上皮细胞和基质细胞上同样存在FSH的受体。这就意味着，内膜组织能够直接“听懂”FSH的指令。
基于这一发现，我们不禁要重新思考疾病的发生机制。子宫内膜异位症是一种高度依赖雌激素的疾病，较高水平的FSHB可能会驱动更强的雌激素合成，从而为异位内膜组织的存活和增殖提供温床。但更深层的可能在于，异常放大的FSH信号或许直接篡改了内膜基质与上皮细胞之间的互动模式，引发了异常的炎症反应或蜕膜化过程。
除此之外，分析还发现，较高水平的SEZ6L2也与子宫内膜异位症的风险增加存在因果关联，而高水平的PAEP则是导致月经过多的直接风险因素。通过反向孟德尔随机化分析的重重考验，并未发现任何显著的反向因果证据，这使得这些蛋白质在发病机制中的因果地位更加稳固。
重塑时间标尺：75个蛋白质构建的预测模型
既然这些外周血中的蛋白质能够如此精准地刻画生理周期，并预示疾病风险，那么我们能否利用它们反向推断身体当前所处的时间坐标？
为了验证这一设想，研究人员将核心队列随机划分为70%的训练集（1931人）和30%的验证集（829人），并引入了LASSO回归模型。经过十折交叉验证的严格筛选，模型最终在接近3000种蛋白质中，锁定了一个由75种蛋白质组成的精干特征集合。
基于这75种蛋白质的权重，研究人员构建出了一个专属的“周期蛋白质评分”（cycleDayProtS）。当把这个评分应用到完全独立的验证集人群中时：该评分与女性自我报告的周期天数展现出了高度的线性相关性。为了更直观地评估这个新型“生物钟”的准确度，研究将其与传统的重量级参照物——血清雌二醇进行了对决。
单次抽血测定的解释力（对生理周期变异的预测度）：
血清雌二醇R² = 0.002
VS
周期蛋白质评分R² = 0.418
差异令人咋舌。当试图将雌二醇加入到蛋白质模型中共同预测时，发现联合模型的预测能力（R²=0.415）并没有得到任何提升。这个数据客观地宣告了一个事实：在判断女性处于生理周期哪一阶段时，由数十种血浆蛋白协同构成的系统性特征，远远比单一的经典激素标志物更具代表性和鲁棒性，它为我们提供了一个极其强大且稳定的分子读数器。
倾听身体的内部对话
这项针对近3000种血浆蛋白的大规模图谱研究，不仅仅是在学术界更新了一个数据库，它更像是一场对女性生理认知视角的重塑。长久以来，医学界习惯于将目光聚焦在骨盆深处的那个狭小腔体内，将子宫和卵巢视作孤立的内分泌器官。而如今，通过高通量蛋白质组学、单细胞测序以及孟德尔随机化等一整套先进分析工具的抽丝剥茧，我们终于能够倾听到从子宫深处传递向全身血液的内部对话。
这些随着日升月落般起伏的蛋白质，证明了生理周期是一个动员了全身免疫系统、血管重塑系统和代谢网络的宏大工程。它们在维持正常生育能力的同时，也可能因为微小的基因扰动而在暗中推动着子宫平滑肌瘤或子宫内膜异位症的发展。
更为重要的是，研究建立的75蛋白预测模型，展现出了在单次抽血中精准捕捉时间刻度的巨大潜力。在未来的女性健康管理、药物临床试验甚至是个性化医疗设计中，这种高精度的分子时钟将成为不可或缺的基础设施。我们不再需要依赖主观且容易存在偏差的记忆，血液中流淌的蛋白质密码，已经为我们记录下了生命节律中最真实、最客观的时间。
参考文献