GhFIR_A08：一个未表征的调控子通过整合到GhWRKY16-GhRL6转录层级中促进棉花纤维起始

《Journal of Advanced Research》：GhFIR_A08: An uncharacterized regulator promotes cotton fiber initiation by integrating into a GhWRKY16-GhRL6 transcriptional hierarchy

　　

棉花（Gossypium spp.）是一种重要的全球经济作物，提高产量仍是主要育种目标。纤维发育受复杂调控网络控制，其中起始阶段（fiber initiation）至关重要，因为它决定纤维数量并最终影响产量。然而，大多数研究集中于上游转录因子（transcription factors, TFs），下游功能基因仍大部分未表征。为探索未注释基因GhFIR_A08（Fiber Initiation Regulator on chromosome A08）的基因功能和调控网络（该基因基于先前的结构变异（structural variant, SV）分析鉴定，并在纤维起始阶段高表达），研究人员比较了当前高质量棉属（Gossypium）基因组组装，以研究FIRs基因家族的起源及转座子介导的拷贝数变异（copy number variation, CNV）机制。通过生成稳定纯合转基因系（敲除和过表达）在受体背景中评估基因功能。利用MBP pull-down结合质谱（mass spectrometry, MS）、酵母双杂交（yeast two-hybrid, Y2H）、双分子荧光互补（bimolecular fluorescence complementation, BiFC）和pull-down实验筛选并验证相互作用蛋白。通过酵母单杂交（yeast one-hybrid, Y1H）、双荧光素酶报告基因（dual-luciferase reporter, Dual-LUC）、电泳迁移率变动分析（electrophoretic mobility shift assay, EMSA）、RNA-seq及相关代谢物测量探索调控网络。研究揭示了该基因家族在棉属中的进化起源和多倍化驱动的扩增，以及一个转座子介导的单倍型特异性GhFIR_A08缺失事件。过表达GhFIR_A08增加了胚珠突起数量并提高了衣分（lint percentage, LP），而敲除则产生相反表型。进一步发现了一个上游调控层级：GhWRKY16诱导另一个激活子GhRL6；GhWRKY16和GhRL6均结合GhFIR_A08启动子的相应位点，形成一个连贯的GhWRKY16-GhRL6-GhFIR_A08模块来激活GhFIR_A08。功能上，GhFIR_A08与Gh14-3-3、GhHSP70、GhSPL13和GhPPR协同作用，通过调节激素信号、蛋白质稳态和氧化还原过程来调控棉花纤维起始。结论：未表征基因GhFIR_A08被揭示为正调控纤维起始以提高LP，并为旨在提高棉花产量的分子育种提供了有价值的遗传资源。

**论文解读：GhFIR_A08——整合进GhWRKY16-GhRL6转录层级调控棉花纤维起始的未表征因子**

**研究背景、存在问题与研究目的**

棉花（Gossypium spp.）是全球重要的天然纺织纤维来源，其产量和品质直接影响纺织产业链的稳定与升级。衣分（lint percentage, LP）作为纤维重量占籽棉比例的关键产量性状，在棉花驯化和高产育种中备受关注。LP主要由胚珠表皮中成功起始的纤维细胞数量决定，这一细胞命运决定过程受转录调控、激素信号和代谢途径构成的多层复杂网络精密调控。尽管已有研究揭示了MYB-bHLH-WDR（MBW）复合物、激素（如生长素、油菜素内酯、乙烯）信号通路以及多组学整合网络在纤维起始中的作用，但当前研究仍存在关键空白：大部分研究聚焦于上游转录因子和信号调控子，而纤维起始是一个涉及细胞快速突起和扩展的生理过程，必然伴随“发育应激”信号，如活性氧（reactive oxygen species, ROS）积累、氧化还原稳态调节以及高效蛋白质合成与质量控制。然而，将上游调控信号和发育应激信号直接转化为细胞物理突起的“分子执行器”仍未得到表征。此外，棉花异源四倍体基因组的多倍化和亚基因组分化导致基因冗余、表达偏差和功能多样化，基于单核苷酸多态性（single nucleotide polymorphism, SNP）的变异解释力有限。相比之下，结构变异（structural variant, SV）可改变基因剂量、破坏编码区或重构调控架构，产生大效应变异。因此，系统鉴定群体水平的结构变异为阐明纤维起始的遗传基础、优先筛选因果变异和关键功能效应子提供了有效策略。

在前期比较基因组分析中，A08染色体倒位区域被解析为三个主要单倍型区块，其中A08_HB2与LP显著关联。携带g0单倍型的品种LP显著降低，并含有一个约60 kb的缺失，导致未注释基因Ghi_A08G06761（GhFIR_A08）丢失。GhFIR_A08在陆地棉TM-1纤维起始阶段具有高度特异的时空表达，但其生物学功能和调控机制未知。本研究旨在确定GhFIR_A08是否真正调控LP，并阐明其潜在机制，从而扩展纤维起始的调控框架，并为基于分子设计的棉花高产育种提供可操作靶点。

**主要关键技术方法**

1. **系统发育与进化分析**：利用BLASTP搜索棉属及外群物种的蛋白质数据库，通过MAFFT多序列比对和IQ-TREE最大似然法构建系统发育树，分析FIRs基因家族进化及多倍化驱动的扩增。
2. **结构变异（SV）鉴定与转座子分析**：收集159个陆地棉、17个海岛棉等种质资源基因组，通过MUMmer比对和SyRI检测缺失事件；利用EDTA注释转座元件（transposable element, TE），并分析LTR反转座子插入时间。
3. **转基因功能验证**：在Jin668受体背景中，通过农杆菌介导法获得GhFIR_A08过表达（OE）和CRISPR/Cas9敲除（KO）稳定纯合转基因系，结合扫描电子显微镜（SEM）观察纤维起始表型，并测定LP。
4. **蛋白互作与转录调控分析**：通过MBP pull-down偶联质谱（MBP pull-down/MS）筛选候选互作蛋白，结合Y2H、BiFC和pull-down验证；利用Y1H、EMSA和Dual-LUC（在棉花原生质体和本氏烟草中）鉴定直接结合并激活GhFIR_A08启动子的上游转录因子，并构建转录层级。
5. **转录组与代谢物分析**：对0 DPA（开花当天）胚珠进行RNA-seq，鉴定差异表达基因（DEGs）并进行GO/KEGG富集分析；定量测定H₂O₂和茉莉酸（jasmonic acid, JA）相关代谢物水平。

**研究结果**

**1. Gypsy/Tekay LTR反转座子导致A08位点GhFIR_A08缺失**
通过系统发育和染色体定位分析，发现FIRs基因家族在二倍体棉花中有2个基因，在异源四倍体陆地棉中有4个基因，且染色体8衍生的拷贝经历更强的纯化选择。对159个陆地棉种质等群体的基因组比对显示，A08缺失事件在栽培棉中频率显著高于半野生棉，缺失断点位于完整的Gypsy/Tekay LTR反转座子内，推测通过非同源末端连接（non-homologous end joining, NHEJ）产生，且该缺失可能在驯化育种中受到人工选择。

**2. GhFIR_A08动态调控纤维起始**
RT-qPCR显示GhFIR_A08在纤维起始关键期（-1至1 DPA）胚珠中特异性高表达。过表达GhFIR_A08显著增加0 DPA胚珠表面纤维初始突起密度并提高LP，而敲除则产生相反表型，表明GhFIR_A08是纤维起始的正调控因子。

**3. GhWRKY16和GhRL6直接转录激活GhFIR_A08**
通过Y1H筛选，发现GhWRKY16（已报道的正调控因子）和GhRL6（MYB相关类型转录因子）可结合GhFIR_A08启动子。EMSA证实两者直接结合其各自顺式元件（W-box和MYB结合位点）。Dual-LUC实验在棉花原生质体和本氏烟草中均证实二者显著激活GhFIR_A08启动子活性。VIGS沉默GhRL6导致GhFIR_A08表达下降和纤维初始数量减少。

**4. 层级调控级联GhWRKY16-GhRL6-GhFIR_A08**
进一步研究表明，GhWRKY16可直接结合并激活GhRL6启动子（Y1H、EMSA、Dual-LUC验证）。共表达实验显示GhWRKY16不仅增强GhRL6启动子活性，还能进一步增强GhRL6对GhFIR_A08启动子的反式激活。在RNAi-GhWRKY16转基因系中，GhRL6和GhFIR_A08表达均显著下调，确立了一个连贯的GhWRKY16-GhRL6-GhFIR_A08调控模块。

**5. 筛选和鉴定GhFIR_A08的潜在相互作用蛋白**
通过MBP pull-down/MS从0 DPA胚珠中鉴定出4个候选互作蛋白：GhHSP70、GhPPR、Gh14-3-3和GhSPL13。亚细胞定位显示GhFIR_A08主要定位于质膜和内质网，互作蛋白也分布在不同亚细胞区域。Y2H、BiFC和pull-down实验验证了这些相互作用。在GhFIR_A08过表达背景中沉默GhHSP70或GhPPR显著降低了纤维初始突起数量，表明这些互作蛋白在GhFIR_A08介导的纤维起始中具有功能必要性。

**6. GhFIR_A08的转录组分析揭示纤维发育中的潜在通路**
RNA-seq分析显示，过表达和敲除系分别有10,371和差异基因，594个核心DEGs共同响应GhFIR_A08扰动。GO/KEGG富集分析显示这些基因显著富集于JA响应、氧化应激、脂肪酸代谢以及蛋白质加工和转运通路。在敲除系中JA相关基因上调，而在过表达系中氧化应激和蛋白质折叠相关基因激活。代谢物测定证实敲除系中H₂O₂水平最高，JA稳态被扰乱（前体OPDA积累而活性衍生物JA-Ile、MeJA减少），表明GhFIR_A08与JA稳态及下游信号输出密切相关。

**讨论与结论**

讨论部分总结了GhFIR_A08蛋白的结构特征：预测为单次跨膜蛋白，C端含疏水跨膜结构域，N端大结构域面向胞外，与内质网/质膜定位一致；其三维结构与拟南芥AtHUP40高度同源，为后续生化功能研究提供线索。进化分析表明，FIRs基因家族在四倍体棉花中经历不对称选择压力，染色体8拷贝受强纯化选择，染色体12拷贝受松弛选择，这种不对称性可能为功能特化提供遗传基础。群体PAV单倍型分析显示，基因缺失单倍型Hap2在北方早熟棉区频率较高，因其与缩短生育期相关，尽管对LP有不利影响，反映了多性状权衡下的育种选择。调控框架方面，研究提出了GhWRKY16-GhRL6-GhFIR_A08多层转录架构，将GhWRKY16延伸至下游执行基因，并将RADIALIS-LIKE因子引入纤维起始网络。互作蛋白分析表明，GhFIR_A08可能作为质膜/内质网关联的“执行节点”，通过招募Gh14-3-3（与油菜素内酯信号相连）、GhHSP70（分子伴侣）、GhPPR（细胞器RNA代谢）和GhSPL13（发育调控）整合激素信号、蛋白质稳态和氧化还原调控，从而促进纤维起始。转录组和代谢数据进一步支持这一模型。

**结论部分翻译**：总之，本研究鉴定GhFIR_A08为一个先前未表征的下游效应子，作为棉花纤维起始的候选正调控因子。通过整合进化、遗传和多组学分析，揭示了一个转座子介导的单倍型特异性缺失事件，强调了结构变异作为多倍体棉花功能多样化关键驱动力的作用。机制上，提出了一个潜在的GhWRKY16-GhRL6-GhFIR_A08调控模块，将上游转录调控与下游蛋白复合体组装及与胚珠突起相关的细胞过程联系起来。这些发现为纤维起始中转录因子中心调控研究与功能执行之间提供了有价值的概念桥梁。尽管如此，其精确的机制定位和下游生化途径仍需进一步的遗传上位性分析、天然背景互作验证和详细的生化表征。除了机制洞察，GhFIR_A08的鉴定为旨在提高衣分的精准育种提供了一个有前景的分子靶点。总之，本工作推进了我们对纤维发育调控的理解，并阐明了结构变异如何塑造多倍体作物农艺性状的进化。