★ 汇报要点（口头汇报提纲）
供汇报时照讲；下方各节为支撑证据与图表。一句话：已用真实多组学数据 + 样本级严谨统计，确认了两条 PE 上调机制轴，给出候选靶点与文献佐证；两条"创新主线"如实判为未支持；并明确需客户补三项数据。
① 开场与定位
- 这是 PE 胎盘多组学项目的当前进展汇报——真实数据、样本级统计、诚实结论。
- 数据：8 例配套样本（对照 Z2/Z3/Z4/Z5 + PE C1/C2/C3/C5）、5 大组学、约 1.3 TB，已落盘分析服务器。
- 交付物：网站 pe.sinogenomics.com（结果 + 执行报告 + 靶点优先级 + 文献佐证）。
② 核心结论（可讲得有底气）
- 找到两条统计确认的 PE 上调程序：
- ▸ 分泌/抗血管轴（sFlt-1/FLT1、ENG、INHBA、FSTL3、CRH、LEP、PAPPA2）——PE 公认核心机制。
- ▸ 干扰素/炎症/抗原呈递轴（IFN-α/γ、CXCL9/10/11、HLA-G、IDO1）。
- 两轴 GSEA 通路水平 FDR<0.05–0.001、跨多种细胞一致、且有已核实文献支持；各轴均给出候选靶点优先级。
③ 方法严谨性（主动说，建立信任）
- 用 样本级 pseudobulk DESeq2 + GSEA，而非逐细胞检验——后者把 6.5 万细胞当独立样本会使"显著基因数"虚高（已纠正）。
- 修正了一处分组标签反置的历史错误（已用 PE 标志物验证：对照=Z、PE=C）。
- 所有数字可溯源；文献 PMID 经 NCBI 独立核实，无编造。
④ 诚实阴性 + 关键请求
- 设计里两条"创新主线"（铁死亡/脂质盾、Retromer/sFlt-1 错误分泌）当前数据不支持——如实说，避免误导后续。
- 受 8 例样本功效限制，仅最强信号过严格 FDR。
- 请客户尽快补三项数据（否则招牌交付物做不了）：① 药物外植体处理转录组（→ 药物逆转打分）；② 空转切片图像（→ Spot 级原位制图）；③ 血清代谢组定量表。
💬 预期问答（被问到时这样答）
| 显著基因为何只有几十个？ | 因为用了正确的样本级统计；逐细胞会虚高到数千，那是伪重复假象。8 例样本下只有最强效应过 FDR，属正常、可信。 |
| 铁死亡 / 脂质那条线呢？ | 当前数据不支持（GSEA 不显著、效应量极小），如实标注，不强行下结论；可作为后续更大样本的探索方向。 |
| 药物靶点 / 逆转分析做了吗？ | 需要外植体药物处理数据，客户尚未交付；数据到位即可做 CMap 逆转打分。当前先给出基于确认轴的候选靶点优先级。 |
| 结论可靠吗？ | 两条确认轴有多组学 + 空间 + 文献三重佐证；靶点优先级为假设生成式（已注明），供验证选点。 |
| 空间分析为何有限？ | 缺切片明场/H&E 图像，无法做 Spot 级原位制图 / RCTD；现用厂商区域级差异结果作支撑。 |
01 执行摘要
一句话结论：在样本级严谨统计下，PE 胎盘的「分泌 / 抗血管生成」程序（sFlt-1/FLT1、ENG、INHBA、FSTL3、CRH、LEP、PAPPA2 等）被一致、显著上调——这是 PE 公认核心机制，在本数据中获跨细胞类型确证。设计中两条「创新主线」（脂质/铁死亡、Retromer）在当前数据与功效下未获支持，如实报告。
① 已确认 FDR<0.05
分泌/抗血管轴 GSEA 在多种细胞一致富集上调：EVT NES 1.63 (q=0.009)、VCT 1.61 (q=0.004)、Hofbauer 1.58 (q=0.01)、全体 1.53 (q=0.04)。pseudobulk 中 CRH、LEP 显著上调（PE 经典标志物）。
② 方法严谨性 自算
以 pseudobulk DESeq2（样本级）为主统计；逐细胞 Wilcoxon 把 6.5 万细胞当独立样本→显著数虚高（如 EVT“8000+ 显著”），本页降级为探索。低功效下用 GSEA 检验通路协同变化。
③ 诚实阴性 未支持
主线1（GPX4/AGPS/GNPAT 铁死亡）与主线2（VPS35/SNX1 Retromer）：样本级 per-gene 效应量极小（GPX4 +0.10、VPS35 +0.08）、GSEA 不显著（q>0.14）。当前数据不支持这两条假设，避免据此误导后续。
⚠️ 功效与缺失数据：仅 8 例样本 → 仅最强效应（log2FC>3 的 CRH/LEP）过 FDR；更细机制需更多样本或更深设计。招牌交付物仍被缺数据阻断：药物外植体转录组（从未交付，无法做药物逆转打分）、空转切片图像（缺失，无法 Spot 级原位制图/RCTD）。
02 项目执行报告
需求来源：企业微信沟通（2026-05-13 起）、腾讯会议转写（2026-06-09）、需求分析报告（2026-06-12）。
🎯 目标与需求
- 跨组学整合解析 PE 胎盘机制，筛选转化靶点。
- 双轨：假设驱动（脂质/铁死亡、Retromer/sFlt-1）+ 无偏数据驱动；收敛门禁 Jaccard>0.15 & Fisher p<0.01。
- 细胞亚群精细化（EVT/SCT/VCT + 免疫/基质/内皮）。
- 关联药物干预（外植体处理 → CMap 逆转打分）。
- 导出可发表级矢量图（SVG/PDF，≥300 DPI）。
🧪 方法与环境
- server199：40 核 / 125 GB / 2× RTX 5000。
- 单细胞：Scanpy（QC/Leiden/UMAP/注释）；差异用 pydeseq2（pseudobulk DESeq2，样本级）为主，逐细胞 Wilcoxon 作探索。
- 通路：gseapy GSEA pre-rank（机制基因集 + Hallmark/KEGG）。
- 蛋白质组：DIA-NN 基因矩阵（n=4/组）Welch t。
- 空转/空代：整合厂商已交付差异/富集结果。
🗓️ 执行时间线
05-13 ~ 06-02 需求对接 → 硬盘交付 1.3 TB，MD5 校验（血清代谢组 003/004 百度网盘补传）。
06-09 / 06-12 首次会议锁定方案；客户百度网盘补充单细胞矩阵、空转报告。
06-16 ~ 06-17 单细胞再分析 + 分组校正 + pseudobulk DESeq2 + GSEA + 蛋白/空转/空代整合。
进行中 完整富集/通路图、更多细胞类型 pseudobulk、出图统一矢量化。
待数据 药物 CMap 逆转、空转 Spot 级制图、原始血清代谢组处理。
03 客户提供资料清单
原始数据经硬盘邮寄（2026-06-02 收讫），落盘 server199:/disk1/BIO/PE/胎盘组学数据/，约 1.3 TB。
| 组学 | 内容 | 规模 | 厂商批次 | 状态 |
|---|---|---|---|---|
| 临床元数据 | 标本采集记录表（8 例） | ~46 KB | — | 已用 |
| 单细胞转录组 | FASTQ + Cell Ranger 矩阵 | ≈613 GB | DZOE2023121956 | 已分析 |
| 空间转录组 | 10x Visium FASTQ + 报告 | ≈479 GB | DZOE2023121143 | 差异已整合 图像缺失 |
| 空间代谢组 | MALDI-MSI 原始+差异+报告 | ≈107 GB | DZOE2024041979 | 差异已整合 |
| Bulk 蛋白质组 | DIA-NN 定量（8 样本） | ≈4.0 GB | DZOE2024081145 | 已分析 |
| 血清代谢组 | 母血/脐血 非靶向 LC/GC-MS（.D） | ≈51 GB | DZOE2024081148 | 待 MS 处理 |
| 外植体药物转录组 | 对照胎盘外植体药物处理 | 待定 | — | 待客户提供 |
04 分析状态看板
✅ 已完成
- 单细胞 QC / 聚类 / 12 类注释
- 分组标签校正（C=PE, Z=对照）
- pseudobulk DESeq2 差异（全体 + 各细胞类型）
- GSEA 通路/机制检验
- Bulk 蛋白质组差异
- 空转 & 空间代谢组差异整合
🔄 进行中
- 完整 GO/KEGG/Reactome 富集制图
- 更多细胞类型 pseudobulk 深化
- 空间代谢×转录 通路级共定位
- 统一矢量出图
⏳ 待客户数据 / 阻断
- 药物外植体转录组 → CMap 逆转打分
- 空转切片图像 → Spot 级制图 / RCTD
- 原始血清代谢组 .D → 专用 MS 处理
- 专家方案签字确认
05 单细胞转录组结果 自算
—
🔧 已修正数据问题：历史中间结果（6-16）分组标签反了（C 被当成对照）。已按样本表 + PE 标志物校正为 对照=Z2/Z3/Z4/Z5、PE=C1/C2/C3/C5，并验证：
图1. 65,828 细胞 UMAP，12 类细胞（含 EVT/SCT/VCT 滋养层）。
图2. 细胞组成（组内占比）：PE 中 SCT 升高（0.9%→14.4%）。注：组成差异部分可能受解离/取材影响。
差异表达（主统计：pseudobulk DESeq2，样本级 n=4 vs 4）
图3. 全体细胞 pseudobulk 火山图。显著上调以 CRH/LEP/PLAC4/KRT5 等为主；FLT1/FSTL3/PAPPA2 虽上调但未过 FDR（功效所限）。
各细胞类型显著基因数 padj<0.05 & |log2FC|>1
| 细胞类型 | 显著 | ↑PE | ↓PE |
|---|
数量普遍较少，反映 8 例样本的真实统计功效；并非分析缺陷。对比逐细胞 Wilcoxon（全体 2783“显著”）即知后者虚高。
探索性（仅作可视化）：逐细胞 Wilcoxon DE、score_genes 评分图等因 pseudoreplication 不作显著性依据；分组校正验证：
06 机制假设评估（GSEA + 多组学） 自算
对设计的机制主线，用 pseudobulk 排序做 GSEA pre-rank（低功效下检验通路协同变化的稳健方法），并逐基因汇总三层证据。
图4. 机制基因集 GSEA NES（PE vs 对照，各细胞类型）；* = FDR<0.05。仅「分泌/抗血管」一致显著上调。
分泌/抗血管 确认
EVT/VCT/Hofbauer/全体 NES 1.5–1.6，FDR<0.05。PE 公认核心机制，本数据稳健支持。
主线1 脂质/铁死亡 未支持
GSEA 不显著（全体 NES 1.03, q=0.62）；GPX4/AGPS/GNPAT 样本级效应量极小。
主线2 Retromer 未支持
NES 方向不一致、均不显著；VPS35/SNX1 几无变化。原“Retromer 下调”假设未获证。
机制基因 · 多组学证据表
PE 上调PE 下调未检出/不显著
| 基因 | 单细胞 pseudobulk log2FC | pseudobulk padj | 蛋白 log2FC | 空转 | 判读 |
|---|
📌 读法：分泌轴基因（FLT1/ENG/INHBA/FSTL3/PAPPA2/CRH/LEP）方向一致上调（GSEA 集体显著）；脂质/Retromer 基因样本级 padj 多不显著、效应量小——单个基因看“都在动”，但严格统计下不成立。蛋白层 n=4，按方向+nominal 解读。
07 通路富集与候选靶点 自算
在已确认方向的基础上给出：（A）数据驱动发现的另一条显著通路；（B）聚焦分泌/抗血管轴的候选靶点优先级（透明打分）。
第二条确认轴（GSEA 发现）：PE 胎盘 干扰素 / 炎症 / 抗原呈递程序显著上调——全体 IFN-α 反应 NES 1.88、IFN-γ 1.86（FDR 0.001）；EVT 抗原加工呈递 1.91、TLR 信号 1.78、胞质 DNA 感知 1.79。提示无菌性炎症 / 天然免疫激活，亦为 PE 公认特征。
GO / KEGG / Reactome 通路富集（全体细胞，PE 上调，FDR<0.25；点大小=前沿基因数，颜色=−log10 FDR）——三大数据库一致指向 I 型干扰素 / 抗病毒信号（IFN-α/β、ISG15、JAK-STAT），独立印证炎症轴。
发现通路（GSEA，FDR 最低）
| 细胞 | 通路 | NES | FDR |
|---|
Hallmark + KEGG，pseudobulk 排序 GSEA；正 NES = PE 上调。
图6. 分泌/抗血管轴候选靶点优先级。标注：LE=GSEA 前沿基因、Prot↑=蛋白同向、ST↑=空转同向、FDR✓=样本级显著。
分泌 / 抗血管轴 — 候选靶点优先级
| # | 基因 | pseudobulk log2FC (最强细胞类型) | 全体 padj | 蛋白 | 空转 | 前沿 | 已知 PE 角色 | 评分 |
|---|
干扰素 / 炎症轴 — 候选靶点优先级
图7. 干扰素/炎症轴候选靶点优先级；前沿基因取自已确认的 IFN / 抗原呈递 GSEA 通路（180 个前沿基因）。
| # | 基因 | log2FC | padj | 蛋白 | 空转 | 免疫角色 | 评分 |
|---|
⚠️ 定位：两轴均为假设生成式优先级——通路集水平已由 GSEA 确认（前沿 q=0.001），但多数单基因受 n=8 功效限制未达 per-gene FDR；评分综合效应量 + 多组学一致性 + 已知生物学，用于指导后续验证选点。药物逆转（CMap）打分仍因外植体数据未交付而无法进行。
08 跨组学整合
—
蛋白质组定量基因数（DIA-NN, 8 样本） · 自算
—
空间转录组 PE vs 对照 区域比较 · 厂商
—
空间代谢组 差异代谢物（PE vs 对照，早+晚） · 厂商
Bulk 蛋白质组 自算
—
图5. 蛋白质组 PE vs 对照火山图（DIA-NN）。ENG/INHBA/FSTL3 等分泌轴蛋白上调（nominal）。
空间转录组（厂商差异）厂商
—
| 区域比较 (PE vs 对照) | ↑ | ↓ |
|---|
机制基因空间上调：CRH/ENG/FLT1/FSTL3/HTRA4/INHBA/PAPPA2（分泌轴，与上方一致）。
空间代谢组差异（厂商）厂商
—
PE vs 对照（早期切片）火山图
PE vs 对照（晚期切片）火山图
🧬 脂质相关：— 提示 PE 胎盘存在代谢物层面的脂质重塑；但需注意：转录层面的脂质酶/铁死亡假设（主线1）并未获支持——两者层级不同，勿混为一谈。真正的 Spot 级「代谢×转录」原位共定位需空转切片图像（缺失），列为后续。
09 局限与下一步
诚实的局限
- 样本量小（n=4/组）：FDR 功效有限，仅最强信号显著；细机制需更多样本。
- 逐细胞统计已纠偏：改用 pseudobulk/GSEA，避免伪重复夸大。
- 组成差异（SCT 升高）可能受解离/取材影响，非纯生物学结论。
- 蛋白/空代为 nominal/厂商结果，跨层一致性为方向性证据。
- 缺数据：药物外植体、空转图像、原始血清代谢组未到位。
下一步
- 本团队：完整富集与通路图；空间代谢×转录通路级整合；候选靶点（聚焦已确认的分泌/抗血管轴）优先级排序；矢量出图。
- 需客户/上游：① 药物外植体转录组（解锁 CMap 逆转）；② 空转切片图像（解锁 Spot 级制图/RCTD）；③ 血清代谢组定量表或授权处理 .D；④ 专家方案确认。
10 候选靶点文献佐证 已核实
为候选靶点附真实文献依据；每条引用均经 NCBI E-utilities 独立核实（PMID ↔ 标题/年份一致），点击跳转 PubMed，无编造。共 — 条。注：HLA-G / IDO1 / CD74 等在 PE 中多为下调或方向有争议，按「失调」理解（见各条）。
分泌 / 抗血管轴
干扰素 / 炎症轴
数据溯源与方法： 自算 单细胞（Scanpy；pseudobulk DESeq2 via pydeseq2；GSEA via gseapy）与蛋白质组差异由本团队在 server199 计算；厂商 空转、空间代谢组差异为测序厂商交付。主统计为样本级；逐细胞结果仅作探索。本页不含任何占位/推测数字；未交付数据对应模块均明确标注。
项目 2605136 · 子痫前期胎盘多组学整合 · 报告生成 2026-06-17 · pe.sinogenomics.com
项目 2605136 · 子痫前期胎盘多组学整合 · 报告生成 2026-06-17 · pe.sinogenomics.com