基于您提供的论文内容，我将按照要求依次回答四个问题。 --- **中文标题** 柚（Citrus grandis‘Wanbai’）GATA基因家族全基因组鉴定与表达分析揭示GA3诱导枝条伸长候选基因

《International Journal of Molecular Sciences》：Genome-Wide Identification and Expression Profiling of the GATA Gene Family Reveals Candidate Genes Associated with GA3-Induced Shoot Elongation in Citrus grandis ‘Wanbai’

　　

---

**摘要翻译**

植物GATA转录因子是参与信号转导、生长发育的锌指蛋白。然而，其在柚（pummelo）株高和赤霉素相关枝条生长中的潜在作用尚未得到充分表征。本研究通过理化性质、染色体分布、系统发育关系、基因结构、保守基序、共线性及假定上游调控区预测顺式作用元件等分析，对CgGATA基因家族成员进行了全基因组鉴定与表征。研究人员采用逆转录定量PCR（RT-qPCR）检测了其在不同组织中的表达谱、长期赤霉素（GA3）处理下‘万白柚’（‘Wanbai Pummelo’）梢尖的表达变化，以及高秆与矮化柚基因型中的表达模式。共鉴定出23个CgGATA基因，并将其划分为四个系统发育亚家族。同一亚家族成员通常具有相似的基因结构和保守基序组成。比较共线性分析发现，柚与甜橙之间的共线GATA基因对多于柚与拟南芥（Arabidopsis thaliana）。假定上游调控区含有与光、植物激素、生长及胁迫响应相关的预测元件。这些基因在不同组织中表现出各异的表达谱，其中多个基因在茎和叶中的相对表达值较高。长期GA3处理促进了幼苗枝条伸长，并在28 d时伴随多个CgGATA基因转录丰度的变化。综合上述表达谱、高秆与矮化对比及预测的赤霉素响应元件，CgGATA11和CgGATA18被优先选为研究柚株高变异和枝条伸长的候选基因。这些发现为功能研究提供了基础，但并未证明任一基因直接调控赤霉素信号传导或株高。

---

**论文解读**

**一、研究背景与科学问题**

柚（Citrus grandis，异名Citrus maxima）是芸香科（Rutaceae）重要的经济果树。柑橘树体结构直接影响栽植密度、冠层光分布以及修剪、喷药和采收效率，因此紧凑型或矮化表型对集约化果园体系具有重要价值。近年来，已鉴定出多个调控柑橘树形的基因，如影响分枝角度和内源激素谱的CsTAC1、与株高相关的BBX22-HY5转录级联，以及调控赤霉素生物合成和株高的miR156-SPL模块。株高受遗传、发育和环境信号的共同调控，其中赤霉素（gibberellins，GAs）是细胞分裂、细胞伸长和茎生长的关键调节因子。植物GATA转录因子是IV型锌指DNA结合蛋白，其成员参与光形态建成、叶绿素生物合成、叶绿体发育、根和叶发育、开花、萌发及环境信号响应。在拟南芥中，GATA因子GNC和GNL/CGA1作用于DELLA蛋白和光敏色素互作因子下游，抑制赤霉素促进的生长。然而，柚GATA家族在长期GA3处理、枝条生长和株高表型对比中的特征尚未被充分研究。因此，开展柚GATA基因家族的全基因组鉴定与表达分析，对于解析其与GA3诱导的枝条伸长及株高变异的关联具有重要意义。

**二、研究内容与主要结论**

研究人员在柚参考基因组中鉴定出23个CgGATA基因（CgGATA1–CgGATA23），并对其进行了系统分析。这些基因编码的蛋白质长度从103个氨基酸（CgGATA5）到542个氨基酸（CgGATA13）不等，分子量介于11.55至60.30 kDa之间。预测的等电点范围为4.68至10.30，其中11个蛋白的等电点大于7.0。不稳定指数分析显示CgGATA14为稳定蛋白，其余22个为不稳定蛋白。所有蛋白均具有负的总平均亲水性（GRAVY）值，且均被预测定位于细胞核。染色体定位表明，22个基因分布于9条已装配染色体上，而CgGATA23位于未锚定的ChrUn序列上；5号染色体含最多基因（5个），1号和4号染色体各4个。系统发育分析显示，23个CgGATA蛋白可分为四个亚家族（I–IV），每个亚家族均包含柚、甜橙和拟南芥的代表，表明主要GATA谱系在物种分化前已确立。基因结构和保守基序分析表明，同一亚家族成员通常具有相似的基序组成和外显子-内含子结构，其中基序1对应保守的GATA锌指结构域，存在于所有CgGATA蛋白中。共线性分析检测到柚基因组内多个染色体间的共线关系，种间比较发现柚与甜橙之间（28对）的共线GATA基因对多于柚与拟南芥之间（19对），这与物种间的亲缘关系一致。顺式作用元件预测表明，22个基因的假定上游调控区（起始密码子上游2000 bp）含有多种光响应、植物激素响应、生长发育及胁迫响应相关的元件，其中赤霉素响应元件存在于13个基因的启动子区域。

RT-qPCR分析揭示了CgGATA基因在不同组织（根、茎、叶、花、果、种子）中的差异化表达模式，多个基因在茎和叶中相对表达较高，部分基因在果实或种子中表达较高。在长期GA3处理实验中，50 mg L?1 GA3处理28 d显著促进了‘万白柚’幼苗的枝条伸长，处理组平均株高为10.04 ± 1.71 cm，对照组为6.65 ± 1.18 cm。在28 d时，15个CgGATA基因的表达发生变化：11个基因（CgGATA1、CgGATA3、CgGATA4、CgGATA7、CgGATA11、CgGATA12、CgGATA14、CgGATA15、CgGATA17、CgGATA18和CgGATA21）表达显著上调，4个基因（CgGATA2、CgGATA6、CgGATA9和CgGATA20）表达显著下调。在高秆‘沙田柚’与矮化‘矮晚柚’的春季梢尖对比中，双因素方差分析显示基因型与发育阶段的交互作用在所有基因区块中均极显著。综合GA3处理表达变化、高矮基因型对比和顺式元件预测，CgGATA11和CgGATA18被优先选为候选基因，两者在GA3处理后表达上调、在高秆基因型中相对表达值较高，且启动子区含有预测的赤霉素响应元件（GARE）。但需强调的是，这些结果仅为功能研究提供假设基础，尚未证明其直接调控赤霉素信号或株高。

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

研究人员采用多学科技术手段开展研究。首先，通过BLAST和HMM搜索（PF00320）从柚参考基因组HWB.v1.0中鉴定GATA家族成员，并经NCBI CDD验证保守结构域。其次，利用ExPASy ProtParam预测蛋白质理化性质，Cell-PLoc 2.0预测亚细胞定位。系统发育分析采用ClustalW v2.1多序列比对和MEGA11邻接法构建进化树（1000次bootstrap）。基因结构和保守基序分别通过TBtools和MEME Suite v5.5.9（最多14个基序）分析。共线性分析使用MCScanX模块。顺式元件预测采用PlantCARE分析起始密码子上游2000 bp序列。表达分析采用RT-qPCR，以β-Actin为内参基因，使用2?ΔΔCt法计算相对表达量。植物材料包括‘万白柚’幼苗（不同组织及GA3处理）、高秆‘沙田柚’和矮化‘矮晚柚’的春季梢尖（S1快速伸长期和S2成熟期），样本均设三次生物学重复。

**四、研究结果**

**2.1 CgGATA基因的鉴定、理化性质与染色体分布**

共鉴定出23个CgGATA基因，其编码蛋白长度、分子量、等电点、不稳定指数、脂肪族指数和亲水性等理化性质存在差异，所有蛋白均预测定位于细胞核。染色体定位显示基因不均匀分布于9条染色体和1条未锚定序列上。

**2.2 GATA蛋白的系统发育关系**

基于柚、甜橙和拟南芥GATA蛋白构建的邻接树将其分为四个亚家族（I–IV），其中亚家族I含10个CgGATA蛋白，亚家族II含8个，亚家族III含4个，亚家族IV仅含CgGATA13。每个亚家族均包含三个物种的代表，且柚与甜橙的同源蛋白聚在一起，支持GATA谱系分化早于物种分歧的观点。

**2.3 CgGATA家族基因结构与保守基序组成**

MEME分析鉴定出14个保守基序，基序1为GATA锌指结构域，存在于所有23个蛋白中。同一亚家族成员通常具有相似的基序组成和基因结构；亚家族III成员具有较复杂的外显子-内含子结构。

**2.4 CgGATA基因的种内及种间共线性**

种内共线性分析检测到1、4、5、7、8、9号染色体上的共线关系。种间分析分别鉴定出19对（柚-拟南芥）和28对（柚-甜橙）共线GATA基因对，与物种亲缘关系远近一致。

**2.5 CgGATA基因假定上游调控区的顺式作用元件预测**

在2000 bp上游区域内预测到28种顺式元件，包括11种光响应相关、9种激素响应相关、3种生长发育相关和5种胁迫响应相关元件。赤霉素响应元件存在于CgGATA4、CgGATA5、CgGATA8、CgGATA9、CgGATA10、CgGATA11、CgGATA12、CgGATA13、CgGATA15、CgGATA17、CgGATA18、CgGATA21和CgGATA22的启动子区。

**2.6 CgGATA基因在不同组织及高矮基因型中的表达谱**

组织表达分析显示CgGATA基因具有不同的组织表达模式，多个基因在茎和叶中相对表达较高。高矮基因型对比中，双因素方差分析显示基因型×阶段交互作用在所有基因区块中均极显著（p < 0.0001），说明高矮基因型间CgGATA表达差异依赖于发育阶段。

**2.7 长期GA3处理对幼苗株高和CgGATA表达的影响**

线性混合效应模型分析表明，GA3处理与时间存在显著交互作用，GA3处理组在28 d时株高显著高于对照组。同一时间点，11个基因表达上调，4个基因下调，8个基因无显著变化。

**2.8 枝条伸长和株高变异研究的CgGATA候选基因优选**

综合表达组织偏好、GA3响应、高矮基因型差异和启动子GARE预测，CgGATA11和CgGATA18被选为优先候选基因。但该优选仅提出功能验证的假说，并非证明其因果作用。

**五、讨论与结论**

讨论部分指出，GATA转录因子在植物发育和信号转导中具有保守性和多样性。本研究的系统发育和共线性分析支持GATA家族在芸香科中的进化保守性。顺式元件预测暗示CgGATA基因可能受光、激素和胁迫信号的调控，但短序列基序的存在并不能证明其调控活性。组织表达谱和小种质对比受采样设计限制，无法完全分离组织效应与植株年龄、发育阶段和采样季节的影响。GA3处理实验仅测定了28 d单一时间点，无法区分直接的转录响应与发育或反馈调节引起的继发变化。讨论引用了拟南芥GNC/GNL和杨树PdeGATA3的功能研究作为CgGATA基因潜在功能的参考，但强调本研究属于候选基因筛选，而非功能验证。结论部分重申：共鉴定出23个CgGATA基因，分为4个亚家族；长期GA3处理促进了幼苗枝条伸长，并在28 d时与15个基因的表达变化相关；综合多维度信息，将CgGATA11和CgGATA18优先列为研究柚枝条伸长和株高变异的候选基因，但其确切作用仍需遗传学、生物化学和生理学实验加以验证。