哈林甜橙染色体水平基因组组装及MADS-box基因调控次生果发育的转录组学解析

《Frontiers in Plant Science》：Chromosome-level genome assembly of Hamlin sweet orange and transcriptomic insights into MADS-box gene regulation of secondary fruit development

　　

引言：MADS-box基因家族是调控植物生长发育尤其是花与果形态发生的关键转录因子（transcription factors, TFs）。然而该家族在甜橙（Citrus sinensis）中的进化特征及其在次生果发育中的具体调控作用仍不清楚。本研究通过比较不同果实表型的栽培种揭示上述机制。 方法：以具有典型次生果表型的新哈尔脐橙（C. sinensis ‘Newhall’，CSN）和无该性状的哈林甜橙（C. sinensis ‘Hamlin’，CSH）为材料，构建CSH高质量染色体水平de novo基因组，开展比较基因组分析，并在芽期、盛花期和落瓣期进行转录组测序及加权基因共表达网络分析（weighted gene co-expression network analysis, WGCNA）。 结果：CSH基因组组装N50为32.89 Mb。比较基因组显示串联重复（tandem duplication, TD）是两栽培种分化的主要驱动力。基因家族鉴定表明CSN含139个、CSH含119个MADS-box基因。转录组与WGCNA鉴定出与次生果发育显著相关的模块，并确定CSN SEP3与CSN AG6为关键候选调控基因。 讨论：该研究为CSH提供高质量参考基因组，并从转录组层面阐释MADS-box基因在柑橘分化与果实形态发生中的作用。CSN SEP3和CSN AG6作为次生果发育关键调控因子的鉴定，为柑橘果实形态分子育种奠定理论基础。

研究背景：柑橘为全球重要经济作物，新哈尔脐橙因优质、香气浓、市场接受度高而成为重要鲜食品种。次生果是脐橙最显著形态标志，由额外心皮轮发育而成，本质上源于心皮数增加，是脐橙特有性状，在拟南芥（Arabidopsis thaliana）等模式植物中无对应表型，导致相关研究相对匮乏。MADS-box基因家族是高度保守的转录因子（transcription factors, TFs）家族，其Type I仅含MADS结构域，参与配子体发育、胚乳形成与胁迫响应；Type II含MADS、K-box与C端结构域，其中MIKC^C型主导花器官发育，MIKC*型主要调控雄配子体发育。尽管MADS-box研究已从模式植物扩展至多种经济作物，但甜橙中该家族进化特征及在次生果发育中的调控网络仍不明确，因此研究人员以同属华盛顿脐橙 lineage 的CSN与CSH为对照系统开展研究。

主要技术方法：研究人员材料采自江西赣州国家脐橙工程技术研究中心种质资源圃与信丰县栽培园，设CSN与CSH两个栽培种，在花芽期、盛花期、落瓣期取样，每栽培种每时期3个生物学重复。技术方法包括HiFi读段经Hifiasm组装contig、Hi-C数据经ALLHIC锚定至近染色体水平并用Juicebox校正；重复注释采用Repbase与de novo转座元件库结合RepeatMasker；结构注释整合同源预测、ab initio预测与RNA-Seq辅助并由Evidence Modeler合并；MADS-box基因用Pfam中K-box（PF01486）与SRF-TF（PF00319）隐马尔可夫模型经HMMER识别并以NCBI CDD校验；系统发育用MAFFT比对、FastTree建最大似然树；转录组用HISAT2比对、featureCounts定量FPKM、DESeq2做差异表达；WGCNA以方差稳定变换（variance stabilizing transformation, VST）归一，按中位绝对偏差（median absolute deviation, MAD）筛高变基因，软阈值β=12建网，取模块隶属度（module membership, MM）与基因显著性（gene significance, GS）大于0.6筛枢纽基因；内源激素用UPLC-MS/MS定量GA、IAA、细胞分裂素、JA、SA与ABA等。

研究结果：

3.1 哈林甜橙基因组的de novo组装与注释：研究人员完成首个球平高连续度CSH染色体水平基因组，总长673.06 Mb，仅312个scaffold，N50 32.89 Mb，未知碱基数1,900 bp，gap数19，BUSCO完整度96.95%。与CSN相比，CSH scaffold数减少约3.25倍、gap减少约32倍。全基因组共线性比对识别131个结构变异（structural variations, SVs），其中串联重复扩张42个、插入36个、缺失26个，染色体CSN_Chr5_1富集24个SV。结论：串联重复扩张是两甜橙栽培种基因组进化与表型分化的主要动力。

3.2 甜橙MADS-box基因全基因组鉴定：研究人员在CSN与CSH中分别鉴定139和119个MADS-box蛋白，系统发育分为Type I（Mα、Mβ、Mγ）与Type II（MIKC*、MIKC^C）；MIKC^C进一步含SEP、AGL6、AP1、AG、SOC1、FLC、XALI、AP3、PI、TT16、AGL15、AGL17、SVP等亚族。CSN在SEP、AP1-FUL、AGL17、SOC1、TT16亚族基因数略多于CSH。理化性质显示多数蛋白理论等电点偏碱、稳定性高，亚细胞定位多为细胞核。结论：甜橙MADS-box家族具进化保守性与功能分化潜力。

3.3 甜橙MADS-box基因特征：染色体定位显示Chr1与Chr2富集MADS-box基因，CSN与CSH整体共线性高，但Chr1上段与Chr2下段有删除或易位。启动子cis-acting元件中CSN的AGL15、SVP光响应元件及Type I家系胁迫元件（MBS、LTR）更丰富。基序与结构分析见Motif 1、Motif 2对应核心MADS域，MIKC型具MADS-K-C架构，部分Mα/Mβ同源基因缺典型K-box。结论：Type I亚族靠串联重复扩张，Type II受更强纯化选择。

3.4 新哈尔脐橙次生果发育的转录组与激素谱：PCA显示时期分离大于栽培种分离，qRT-PCR与RNA-seq的R²=0.844。MADS-box差异基因中CSN_AG_1、CSN_TT16_3、CSN_SEP_3、CSN_AGL_3在CSN的T1/T2特异高表达。激素方面，盛花期CSH的GA（12.74 ng/g）高于CSN（3.55 ng/g），CSN的转运型细胞分裂素tZR（20.97 ng/g）低于CSH（49.93 ng/g）；CSN的tZR/ABA比0.17高于CSH的0.057，CSH的GA/ABA比0.035高于CSN的0.012。结论：CSN盛花期高tZR驱动花托分生组织过度增殖，为次生果形成提供生理基础。

3.5 WGCNA识别新哈尔脐橙次生果发育关键调控基因：β=12下得15个模块，紫色模块与目标性状正相（r=0.52, p=0.027），T2期CSN显著上调；以MM>0.6且GS>0.6筛得10个hub基因，其中MADS-box基因为CSN_SEP_3（evm.model.chr3.1.2090）与CSN_AG_6（evm.model.chr8.1.2342）。131个SV均不位于这些hub基因上下游5 kb。结论：CSN_SEP_3与CSN_AG_6是次生果发育关键候选基因，其差异表达更可能来自细微序列变异而非大尺度SV。

讨论部分总结：研究人员指出脐橙次生果发育可由心皮身份基因过表达驱动；CSN盛花期高tZR使WUSCHEL（WUS）未被AG完全沉默，AG与WUS共表达窗口延长分生活性，而CSH无高细胞分裂素故AG正常抑制WUS。蓝色模块中SOC1/SVP/AGL17-like下调可解除对AG/SEP通路抑制，构成“抑制释放”模型。CSH基因组N50 32.89 Mb、仅19 gap为柑橘最高连续度组装；131个SV中串联重复扩张占32.1%，Chr5_1的SV热点与果实品质及NLR抗病基因簇QTL重叠。CSN较CSH多20个MADS-box基因，支持无性繁殖芽变中Type I经串联重复扩张、Type II受纯化选择的进化模式。

结论部分翻译：本研究对哈林甜橙进行染色体水平基因组de novo组装与注释，并系统表征哈林甜橙与新哈尔脐橙中MADS-box基因的鉴定与基因结构特征。通过整合转录组与激素谱，研究人员鉴定出新哈尔脐橙次生果发育相关的关键MADS-box基因。然而，这些MADS-box基因的潜在调控机制仍待未来研究阐明。论文发表于《Frontiers in Plant Science》。