《Plants》：CRISPR/Cas9-Mediated pds Knockout in Potato Reveals Network-Level Transcriptomic Reorganization Beyond Pigment Loss
编辑推荐：
本研究通过CRISPR/Cas9技术敲除马铃薯八氢番茄红素脱氢酶（pds）基因，结合RNA测序（RNA-seq）和加权基因共表达网络分析（WGCNA），系统解析了白化表型背后超越色素代谢的全局转录重编程。研究发现白化组织存在9700余个差异表达基因（DEGs），涉及光合作用抑制、碳代谢紊乱、胁迫响应激活及激素信号通路重塑，揭示白化实为质体功能失调触发的多层级调控网络失衡。该研究为pds作为基因编辑可视化标记的合理应用提供了新见解，并为作物基因编辑的转录组效应评估建立了网络分析框架。
引言
马铃薯（Solanum tuberosum L.）是全球第四大粮食作物，其二倍体杂交自交系和真实种子繁殖技术为遗传改良提供了新途径。然而，基因型依赖的转化和再生效率仍是CRISPR技术应用的瓶颈。八氢番茄红素脱氢酶（pds）基因因其功能缺失导致白化表型，常被用作基因编辑的可视化标记。尽管色素缺失的生化机制明确，但pds敲除是否引发类胡萝卜素合成途径之外的转录组变化尚不清楚。本研究旨在建立马铃薯pds靶向的CRISPR/Cas9编辑平台，通过表型分类和转录组分析，揭示白化表型的分子基础及其超越色素代谢的调控网络。
结果
2.1. pds靶向CRISPR/Cas9编辑系统的建立
研究团队构建了二元表达载体pCAMBIA2300-CAS9-8964，其sgRNA靶向pds基因第四外显子，Cas9切割位点位于PAM序列上游3 bp。载体包含AtU6启动子驱动的sgRNA支架和拟南芥泛素（UBQ）启动子驱动的FLAG标签Cas9编码序列。通过农杆菌介导转化二倍体马铃薯基因型DM，获得21个再生芽，经T-DNA诊断PCR验证17个为转基因阳性，转化效率达80.95%。其中15个株系成功生根，用于后续表型和分子验证。
2.2. 靶点基因分型与白化表型
对15个阳性转基因株系进行pds靶点PCR扩增和Sanger测序，结合TIDE分析发现11个株系存在编辑事件。编辑类型包括插入/缺失（indel）和碱基替换，其中移码突变占44.4%，可能导致蛋白质功能丧失。表型评估采用白化指数（AlbinoLevel 0-2）分级：完全绿色为0级，部分白化分枝为1级，完全白化植株为2级。两个株系（DM1-1和DM1-18）出现白化侧枝（AlbinoLevel 1），主茎仍保持绿色。白化分枝表现为矮化、叶片极小、叶缘粉红色色素沉积，继代培养后仍维持白化表型。移码突变和较大缺失与白化表型相关趋势明显，但基因型-表型关系非绝对。
2.3. RNA-seq数据质量评估与全局转录组概览
对DM1-18株系的白化幼苗（StAlbino）、非白化绿色组织（StUnchanged）和野生型对照（StCK）进行RNA-seq分析，每组三个生物学重复。测序数据质量高，与参考基因组比对效率达85.09%-88.70%。主成分分析（PCA）和Pearson相关性显示组间分离良好。转录组数据还用于优化DM v8.1参考基因组的基因结构注释，更新5876个基因的UTR区域，并识别2262个新基因。pds基因座注释存在不一致，读长覆盖分析支持其真实外显子区域在StAlbino中表达显著降低。差异表达基因（DEGs）分析显示，StAlbino与StCK和StUnchanged比较均出现超过9000个DEGs，而StCK与StUnchanged间仅4319个DEGs。通过交集分析定义"白化核心集"（Albino Core Set）和"非白化特异集"（Non-Albino Specific Set），用于区分表型相关转录特征。
2.4. 功能富集揭示部分与完全白化间的差异转录重编程
基因本体（GO）和KEGG富集分析显示，StCK与StUnchanged比较中，光合作用相关过程（如光捕获、光系统I）显著富集，且呈负向富集（NES < 0），反映部分白化时光合模块减弱。同时，解毒和激素信号相关基因集正向富集（NES > 0），包括III类过氧化物酶、谷胱甘肽S-转移酶和ERF/AP2、WRKY、bZIP转录因子激活。StCK与StAlbino比较则呈现更广泛的转录变化，涉及细胞壁生物合成、激素信号、次级代谢、多糖代谢和质体胁迫/伴侣通路。KEGG分析显示植物激素信号、淀粉和蔗糖代谢、碳固定等通路显著富集，表明完全白化伴随叶绿体功能破坏、氧化胁迫响应增强及碳分配向防御相关过程重组。
2.5. WGCNA模块识别与表型关联
基于方差稳定计数（DESeq2-vst）构建加权基因共表达网络，软阈值功率β=74时满足无标度准则（R2=0.83）。动态树切割和模块特征向量（ME）合并后获得18个模块。将样本按白化严重程度（StCK=0, StUnchanged=1, StAlbino=2）排序，相关性分析发现9个模块与表型显著相关。MEorangered4（r = -0.998）和MEblack（r = -0.996）等模块随白化加剧表达下调，而MEcoral2（r = 0.962）和MEdarkgreen（r = 0.855）等模块表达上调。有序组趋势检验（Jonckheere-Terpstra）证实多数模块在白化梯度中呈现单调变化。
2.6. 性状相关模块的功能富集揭示协同代谢与信号转变
KEGG富集分析显示，MEdarkgreen模块基因随白化加剧上调，富集于DNA复制、核糖体、蛋白酶体及芪类/姜辣素生物合成通路，提示白化进程中次级代谢和物理屏障相关通路激活。MEorangered4模块显著富集淀粉和蔗糖代谢（q = 3.38×10-4），其基因表达随白化单调下调，包括碳水化合物代谢酶和抗性相关基因。MEblack模块持续下调，富集MAPK信号和蛋白质稳态通路（q = 0.0379），包含内质网蛋白折叠、泛素连接酶及WRKY/ERF/bZIP转录因子。MEgrey60模块富集植物激素信号转导（q = 0.0198），涵盖生长素、油菜素甾醇、赤霉素和乙烯信号核心组分，表明白化组织中激素感知和转录调控广泛衰减。模块与DEG集交叉验证显示，MEdarkgreen与白化核心集重叠显著（OR = 23.65），MEorangered4富集于全局响应集（OR = 5.54），而MEblack、MEgrey60和MEbrown4主要与转基因背景集（AC-only）相关。
讨论
3.1. 白化表型伴随广泛转录与代谢改变
pds敲除导致的马铃薯白化表型不仅是色素缺乏综合征，还涉及多途径转录和代谢改变。白化株系初级碳代谢（如MEorangered4）下调，而次级代谢和细胞壁重塑（如MEdarkgreen）激活，类似质体缺陷突变体中的逆行信号响应。pds缺失可能通过消耗类胡萝卜素衍生的脱落酸前体（apocarotenoids）和积累活性氧（ROS）破坏质体-细胞核通信，导致核基因表达重编程。白化表型可视为质体功能失调下调控缓冲失效和代谢失衡的结果。
3.2. 白化植物的被动损伤或程序化响应
模块行为更支持程序化调控响应而非被动损伤。胁迫富集模块（MEgrey60/MEblack/MEdarkgreen）中过氧化氢酶（CAT）、超氧化物歧化酶（SOD）、抗坏血酸过氧化物酶（APX）和热激蛋白（HSP）家族基因诱导，与StAlbino中氧化负荷和蛋白质稳态需求升高一致。光合亚基（如RBCS、PsaE）在MEorangered4中协调下调。激素信号模块（如MEgrey60）中ABA/JA响应转录因子（ERF/WRKY/bZIP）诱导，符合质体胁迫下激素介导的 rewiring。StUnchanged样本虽无可见白化，但呈现模块漂移，表明表型缺失不等于分子扰动缺失，可能代表接近转录不稳定性阈值的中间状态。
3.3. 重新审视基于pds的可视化筛选工具
pds敲除的白化表型虽便于筛选，但转录组异质性挑战了仅凭表型的筛选策略。"正常"编辑植株（StUnchanged）可能携带显著转录变化，导致假阴性；而可见白化可能源于载体整合效应或组培胁迫，造成假阳性。建议整合表型与分子标记（如模块水平读数、早期标记基因、光系统亚基）以提高筛选特异性。未来可开发荧光报告系统或代谢物比色法辅助编辑效率评估。
材料与方法
4.1. 组织培养与CRISPR/Cas9敲除载体构建
二倍体马铃薯（Solanum phureja DM1）于MS培养基保存。通过CRISPOR设计sgRNA（靶点5′-TCACAAACCGATACTACTGG-3′），构建pCAMBIA2300-CAS9-8964载体，经同源重组克隆和测序验证。
4.2. 农杆菌介导的马铃薯转化
农杆菌GV3101携带载体侵染马铃薯茎段，共培养后转移至含卡那霉素和Timentin的选择培养基诱导抗性芽，生根后获得转基因植株。
4.3. 阳性转化子鉴定与表型评估
基因组DNA提取后PCR检测pds靶点，Sanger测序结合TIDE分析编辑效率。表型按白化指数分级记录。
4.4. TIDE分析编辑效率
Sanger测序色谱图经TIDE平台解卷积，评估插入/缺失频率和移码突变比例。
4.5. 转录组测序与差异表达分析
TRIzol提取总RNA，建库后Illumina PE150测序。质控后与参考基因组比对，DESeq2进行差异表达分析（|log2FC| ≥ 1, FDR < 0.01）。功能富集采用GO/KEGG过表征分析（ORA）和基因集富集分析（GSEA）。
4.6. 共表达网络分析
基于VST矩阵构建WGCNA网络，软阈值功率β=74。模块特征向量与白化表型关联分析，功能富集评估模块生物学意义。
结论
本研究建立了马铃薯pds靶向CRISPR/Cas9平台，并通过网络分析阐释了白化表型的多层级调控基础。白化组织呈现超越色素代谢的转录重编程，涉及碳代谢抑制、胁迫信号衰减和激素整合调控等功能层级。这些发现将pds从简单可视化标记提升为量化系统响应指标，为作物基因编辑的转录组效应评估提供了新框架。
