在藤上和藤下成熟的番茄果实中，基因表达的遗传变异性：与最差异显著品种中不同转录模式相关的顺式调控元件

《Plants》：Genetic Variability of Gene Expression in Tomato Fruits Ripened on and off the Vine: Cis-Regulatory Elements Associated with Differential Transcription Patterns in the Most Discrepant Variety

番茄果实成熟过程中，基因表达差异及调控元件分析。采用cDNA-AFLP结合生物信息学工具，比较了栽培品种CAI和突变体NOR与野生种LA0722、LA1385在植物上成熟与离体成熟时的基因表达模式，发现显著差异并识别最不一致的品种CAI。通过顺式调控元件分析，揭示植物上成熟基因的调控元件多位于远端（-1500至-1001bp），而离体成熟基因的近端调控元件（-501至+1bp）更丰富，且胁迫响应元件占主导。研究结果为番茄耐储性遗传改良提供新依据。

本研究聚焦于番茄果实成熟过程中基因表达调控机制的探索，通过整合传统分子技术与现代生物信息学方法，揭示了不同基因型在树上成熟与离树成熟条件下的转录组差异及调控元件特征。研究采用cDNA-AFLP技术结合广义 Procrustes分析（GPA）和 cis-调控元件预测工具，系统比较了栽培品种CAI、突变体NOR以及野生近缘种LA0722和LA1385在两种成熟模式下的基因表达谱，并深入解析了关键基因的调控序列。

### 研究背景与核心问题
番茄作为重要的经济作物，其果实品质受多因素调控，其中成熟环境（树上或离树）对品质的影响尤为显著。树上成熟的果实因持续激素信号刺激和水分供应，通常保留更优的感官品质（颜色、风味、质地等）和更长的货架期。然而，传统研究多关注单一成熟环境下的基因表达，缺乏对品种间遗传差异与调控网络整合性研究的深入探讨。本研究的核心目标在于：
1. 解析不同基因型在两种成熟模式下的差异表达基因
2. 识别具有显著遗传变异的品种
3. 揭示差异表达基因的调控元件特征及其功能关联

### 关键技术路径
研究采用多维度技术整合策略：
- **cDNA-AFLP转录组分析**：通过限制性内切酶酶切结合选择性扩增，生成覆盖全基因组的表达谱，特别适用于比较多个基因型的转录组差异。
- **广义 Procrustes分析（GPA）**：构建包含表型数据（货架期、成熟阶段天数）与分子数据的整合分析框架，实现表型-转录组-基因型的多维关联。
- ** cis-调控元件预测系统**：结合PlantCARE和PLACE数据库，对差异基因的启动子区域进行精细解析，定位关键调控元件。
- **RT-qPCR验证体系**：通过设计特异性引物对关键基因进行定量验证，确保转录组分析的可靠性。

### 主要研究发现
#### 一、转录组特征与品种差异
1. **基因表达总量差异**：
- 树上成熟果实（如CAI品种）表现出更高的转录活性，平均检测到132.4个转录相关片段（TDFs），其中NOR突变体因转录因子干扰，实际检测到882个TDFs，显著高于其他品种。
- 离树成熟果实中TDFs数量平均减少7.6%，可能与乙烯信号减弱相关（离树后乙烯浓度从5.4ppm降至1.4ppm）。

2. **品种特异性表达特征**：
- **CAI栽培品种**：在6种引物组合中均呈现最低TDFs数量（694个），但通过GPA分析被识别为最不一致的品种，其转录组与表型数据（货架期、成熟阶段）的吻合度最低（共识值最高达0.83）。
- **NOR突变体**：虽然属于栽培品种改良目标，但其转录活性最高（882个TDFs），且包含独特的胁迫响应基因。
- **野生品种LA0722/LA1385**：在特定引物组合（如A/B组合）下表现出高可变表达特征，但整体基因表达量低于突变体NOR。

#### 二、成熟模式与基因表达调控
1. **时空表达模式分化**：
- 树上成熟相关基因在启动子近端（-501~+1bp）富集42%的胁迫响应元件（如HSE、ABA响应元件），而离树成熟基因在远端调控区（-1500~ -1001bp）分布更密集（48.8%的激素相关元件）。
- 关键调控元件分布差异导致两种成熟模式下的基因表达时序不同：树上成熟基因受近端强调控元件驱动，离树成熟基因则依赖远端组合调控。

2. **核心差异基因群**：
- **胁迫响应基因**：包括热休克蛋白（HSP）基因Solyc03g115230、过氧化物酶（GPX）基因Solyc08g080940，这些基因在离树成熟中表达量显著提升（qPCR验证显示上调2-3倍）。
- **激素信号通路基因**：Solyc11g020040（DNA结合蛋白）和Solyc03g083910（运输蛋白）在离树成熟中特异性表达，其启动子区域含有高浓度茉莉酸（JA）响应元件（CGTCA/TGACG）。

#### 三、调控元件功能解析
1. **元件分布规律**：
- 树上成熟基因：启动子区（-500~+1bp）以短串联元件为主（6bp长度占比43.8%），包含高密度ABA响应元件（ABRE）和茉莉酸响应元件（茉莉酸甲酯响应元件）。
- 离树成熟基因：远端调控区（-1500~ -1001bp）富含重复序列，其中35%为未知功能元件（如AAGAA、AP-1等），可能与发育阶段特异性调控相关。

2. **功能关联网络**：
- **胁迫响应模块**：包含3类核心元件（HSE、ABA响应、JA响应），通过形成转录因子结合位点网络，协同调控HSP和GPX的表达。
- **代谢调控冗余性**：研究发现，同一基因家族（如ClpB热休克蛋白）在野生品种中可能存在多态性调控元件，形成功能冗余的调控网络。

### 创新性突破与产业应用
1. **GPA分析法的优化应用**：
- 首次将表型数据（货架期、成熟阶段天数）与转录组数据（6种引物组合的TDFs）整合进行三维GPA分析，成功筛选出CAI品种作为最不一致对象（GPA共识值达0.82，残差值最高）。
- 通过残差值排序发现，LA0722在三个维度（表型、树上/离树转录组）的变异度最高（总残差值3.24），但未达到最不一致品种的筛选阈值。

2. **调控元件的定向改造策略**：
- 在CAI品种中，离树成熟相关基因（如Solyc11g020040）的启动子区存在独特的TAAT框元件，与野生品种的TATA框形成互补调控结构。
- 通过设计TAAT/TATA双功能启动子，成功在NOR突变体中实现ABA响应基因的上调表达（实验数据显示提高47%的转录效率）。

### 研究局限与未来方向
1. **技术局限性**：
- cDNA-AFLP对低丰度基因检测灵敏度有限（<5%表达量基因难以捕获）。
- 元件预测工具（PlantCARE）对短串联重复元件（<8bp）识别率不足（约32%漏检率）。

2. **后续研究方向**：
- 建立离体成熟基因的时空表达图谱（0-7天关键节点监测）
- 开发基于TAAT/TATA双框的基因编辑载体（CRISPR-Cas9靶向元件改造）
- 构建货架期相关的多组学整合数据库（转录组+表观组+代谢组）

### 结论
本研究首次系统揭示了栽培品种与野生种在离体成熟调控网络中的本质差异，发现CAI品种的转录组调控网络存在关键断裂点，其启动子区TAAT框元件与野生品种的TATA框形成互补调控结构。该发现为：
1. 开发基于胁迫响应基因的货架期改良技术（如过氧化物酶基因过表达工程）
2. 设计双功能启动子系统（同时响应ABA和茉莉酸信号）
3. 建立多维度GPA分析模型（整合转录组、代谢组、表型数据）

研究结果为番茄品质改良提供了新的分子靶点，特别是在离体成熟调控网络中发现的冗余性元件（如未知功能AAGAA元件），可能成为突破现有遗传改良瓶颈的关键突破口。