2号染色体位点调控玫瑰花瓣农杆菌介导的瞬时GFP表达

《Plant Cell Reports》：A chromosome 2 locus controls transient Agrobacterium-mediated GFP expression in rose petals

　　

瞬时表达分析在玫瑰基因功能验证中具有重要价值，而玫瑰的稳定转化仍繁琐且依赖基因型。研究人员在96个玫瑰基因型的花瓣中评估了农杆菌(Agrobacterium)介导的瞬时绿色荧光蛋白(GFP)表达，并利用全基因组关联分析鉴定相关的宿主基因组区域。GFP表达在浸润后3天(dpi)和5天(dpi)按0至4的半定量标准评分。各基因型在两个时间点表达差异显著，且3 dpi与5 dpi评分强相关(r = 0.90, p < 0.001)。5 dpi时平均评分介于0至3.98之间，其中Sebastian Kneipp、Friesia和Comtessa AL等品种为最允许型基因型。基于37,161个高质量单核苷酸多态性(SNP)的GWAS在2号染色体69–73 Mbp处揭示主要关联，并在未锚定0号染色体contig上发现可能对应同一区域的次峰，顶端SNP达-log10(P) = 7.85。该区间作为高优先级候选区域包含多个潜在宿主因子，注释涉及生长素运输、泛素介导的蛋白周转、核糖核酸(RNA)代谢、膜相关防御反应及核运输等基因。结果为花瓣瞬时分析提供了实用玫瑰基因型，并为解析宿主对农杆菌介导瞬时表达的调控提供了遗传切入点。

研究背景与立项依据

玫瑰是重要经济观赏作物，其基因组资源日益丰富，但稳定转化与再生仍受限于基因型，导致花色、香气、衰老、胁迫响应及株型等性状的基因功能验证困难。农杆菌(Agrobacterium)介导基因转移效率取决于细菌附着、转移DNA(T-DNA)转运、核导入、表达及整合等一系列宿主与细菌过程，宿主基因型可显著影响这些环节。瞬时表达分析能规避稳定转基因耗时问题，在本氏烟(Nicotiana benthamiana)等体系中广泛用于启动子与候选基因快速测试；玫瑰花瓣与叶片的农杆菌浸润瞬时体系亦已建立，并服务于离层、脱水耐性与色素合成分析。然而玫瑰花瓣瞬时表达基因型依赖性的遗传基础尚未在基因组水平解析。全基因组关联分析(GWAS)适合多样种质的表型-标记关联，配合蔷薇(Rosa chinensis)参考基因组与WagRhSNP 68K Axiom阵列可实现四倍体玫瑰面板关联作图。基于此，研究人员在96基因型玫瑰关联面板中量化花瓣农杆菌浸润后瞬时GFP表达，判断变异是否具遗传基础，并筛选适用基因型、定位主效区域及提名候选基因。

论文发表于《Plant Cell Reports》，研究通过对96个以四倍体为主的玫瑰基因型（87四倍体、8三倍体、1二倍体）花瓣在初开至全开期(S4–S5)进行农杆菌GV3101携带UBQ10启动子驱动intron-GFP载体浸润，3与5 dpi蓝光下按0–4级半定量评分，表型经Box–Cox变换后采用GWASpoly（适配自多倍体、检验加性与单式显性模型）基于37,161个SNP结合亲缘矩阵与两主成分混合模型作图，显著性阈值取Meff校正-log₁₀(P)=5.54，关联区间借Rosa chinensis基因组坐标注释，显著SNP侧翼基因经BLASTx比对GDR与NCBI蛋白库提名。该文核心结论为：玫瑰花瓣瞬时GFP表达呈强基因型依赖；5 dpi表型稳健且与3 dpi高度相关；GW题定位到2号染色体69–73 Mbp主效区间（顶端SNP RhK5_12481_722P于70.13 Mbp，-log₁₀(P)=7.85，单式显性模型），0号染色体峰经BLAST确认为同一物理区域未锚定contig；候选基因涵盖ABCB19、Cullin3、DExH-box RNA解旋酶、RNA结合蛋白、过敏反应诱导膜蛋白、单脱氢抗坏血酸还原酶4及核孔GP210样等，涉及生长素运输、泛素-蛋白酶体、RNA代谢、核运输与防御；等位剂量分析显示部分SNP低剂量或特定纯合/四式剂量类群GFP评分显著升高；该2号区域与同面板香气相关标记簇重叠但GFP与香气仅弱相关(r=0.27)，暗示同一广区间不同连锁位点分别调控。研究意义在于提供Sebastian Kneipp、Friesia、Comtessa AL等高允许型基因型及染色体2标记，加速玫瑰功能基因组学早期验证，并为宿主调控农杆菌瞬时表达提供可验证遗传入口。

主要关键技术方法概括

研究人员采用德国联邦植物品种局栽培的96基因型玫瑰关联面板（三重复随机区组、取中部花瓣S4–S5期），以农杆菌GV3101浸润后3/5 dpi半定量GFP评分；基因分型使用WagRhSNP 68K Axiom SNP阵列，滤后37,161 SNP；表型Box–Cox变换，GWASpoly自多倍体混合模型（加性/单式显性、留染色体外亲缘矩阵+2 PC、Meff阈值-log₁₀(P)=5.54）定位；显著区间依Rosa chinensis 2018组装坐标注释，转录本BLASTx比对GDR/NCBI提名候选；等位剂量效应借多倍体剂量类群箱线分析；表型间相关用Pearson/Spearman检验。

研究结果

花瓣瞬时GFP表达基因型差异显著：96基因型3 dpi均分0–3.81、5 dpi 0–3.98；33基因型3 dpi为零信号，13基因型5 dpi仍零；5 dpi强表达类增多，Sebastian Kneipp达3.98；Kruskal–Wallis检验3/5 dpi基因型效应均极显著(p≈10^-15)；3与5 dpi评分r=0.90，故取5 dpi作GWAS端点。

主要关联映射至2号染色体：曼哈顿图显示2号69–73 Mbp主峰（19显著SNP）及0号contig次峰；顶端RhK5_12481_722P 70.13 Mbp，-log₁₀(P)=7.85（单式显性）；12条0号显著SNP转录本BLAST回贴2号同源位置，确认物理同区。

关联区间候选基因：21显著SNP落于ABCB19同源、Cullin3同源、DExH-box解旋酶、RNA结合蛋白、HR响应膜蛋白、MDAR4、核孔GP210样及未知蛋白；因面板连锁衰减达兆碱基层级，区间而非单基因解释信号。

代表SNP等位剂量分层：Rh12GR_23266_2817P与RhK5_12481_722P低剂量类中位GFP更高；RhK5_3122_150P/RhK5_2002_1055P特定高剂量类最高；高表达剂量类中54.3%栽培种评分≥2。

染色体2区间与同面板香气标记簇重叠但GFP-香气仅r=0.27（方差解释<8%），分布异质，提示广区间不同连锁位点分别控制。

讨论部分总结

基因型依赖瞬时表达是实用筛选性状：高允许型（Sebastian Kneipp、Friesia、Comtessa AL）可直接用于花瓣功能分析，规避稳定转化瓶颈；5 dpi动态范围优于3 dpi。2号69–73 Mbp为主效区间，0号峰为同区未锚定contig；96面板分辨率有限，区间含多候选而非单因果基因。候选具生物学合理性：Cullin3/泛素系统参与T复合体蛋白清除；RNA解旋酶/结合蛋白影响转录本加工；核孔GP210样关涉T-DNA核导入；ABCB19关涉生长素运输与组织生理；HR样蛋白关涉免疫抑制转化。表型仅为可见瞬时GFP，非细菌迁移/稳定整合，低表达可能源于进入、T-DNA转运、沉默或防御激活；评分系统足以揭主效QTL但需影像定量、组织学及T-DNA负荷测量提升分辨率；扩大面板与更新组装可缩小区间。

结论译文：多样玫瑰面板揭示农杆菌浸润花瓣后瞬时GFP表达广泛遗传变异；GWAS检出2号染色体主效区及可能同区0号峰；关联标记与高表达基因型即时服务玫瑰功能基因组学，区间候选基因为宿主调控农杆菌瞬时表达提供可检验假说。