《Frontiers in Plant Science》：Transcriptomic and biochemical analyses revealed the dynamic regulatory mechanisms underlying leaf morphogenesis in Camellia sinensis
编辑推荐：
茶树叶片发育是一个以连续形态发生为特征的生物学过程，直接影响茶叶产量和品质。尽管这些阶段至关重要，但调控这些发育阶段的动态机制仍知之甚少。在本研究中，研究人员利用茶树的芽、第一叶、第二叶和第三叶作为实验材料。通过整合转录组测序与生化分析，研究人员系统阐明了叶片
茶树叶片发育是一个以连续形态发生为特征的生物学过程，直接影响茶叶产量和品质。尽管这些阶段至关重要，但调控这些发育阶段的动态机制仍知之甚少。在本研究中，研究人员利用茶树的芽、第一叶、第二叶和第三叶作为实验材料。通过整合转录组测序与生化分析，研究人员系统阐明了叶片形态发生背后的转录-代谢共调控网络。研究结果揭示，不同发育阶段叶片转录组存在显著差异，芽与成熟叶之间的分化最为明显。芽与第三叶的比较鉴定出最多数量的差异表达基因（6,854个），这些基因主要集中在KEGG代谢通路中。核心代谢物，如儿茶素(catechins)、游离氨基酸(free amino acids)和生物碱(alkaloids)，主要在芽和幼叶中富集，且大多数随叶片成熟而显著下降。然而，O-甲基化儿茶素(O-methylated catechins)在成熟叶中显著富集。这些核心代谢物的表达受到苯丙烷和类黄酮通路中的差异表达基因（CHS、ANS、SCPL、OMTs等）、茶氨酸生物合成基因（ALaDC、ALaAT、TS等）以及嘌呤生物碱生物合成基因（IMPDH、SAMS、TCS等）的调控。相关性分析表明，大多数游离氨基酸（L-The、Glu、Gln、Asp、Ala和GABA）与生物碱和儿茶素（C、GC、CG、EC和ECG）呈显著正相关，而与O-甲基化儿茶素成分的相关性较弱。实时定量PCR分析显示，108个基因显著差异表达，其功能主要涉及细胞壁合成与重塑（CER3、PRP4、EXPA1等）、光合作用（ndhL）、植物激素信号调控（SAUR23、CYP85A1、IAA14等）、细胞周期（CYCD1-1、CYCA2-4）、信号转导（MKK6、MLO6）、防御与次生代谢（STR1）以及转录因子（GRF1、MYB17、bHLH94等）。此外，与叶片发育相关的功能基因启动子区域显著富集了MYB和MYC转录因子，以及对激素（ABRE、ERE、TCA-element）、光照（GT1-motif、GATA-motif、TCT-motif）和胁迫信号（ARE、STRE、W box）的顺式调控元件。本研究鉴定了茶树叶片形态发生过程中的核心基因和代谢机制，为茶叶优质栽培提供了理论基础。
茶树叶片的发育是一个依次经历形态建成、功能成熟的生物学过程，直接决定茶叶的产量和品质。尽管已有研究初步揭示了叶片发育过程中基因表达与代谢物积累的关系，但整个发育进程（从芽到成熟叶）的动态调控机制，特别是转录调控与生化代谢之间的协同作用，仍未得到系统阐明。为此，研究人员以茶树品种‘Lingtou Dancong’（LTDC）的芽（B）、第一叶（L1）、第二叶（L2）和第三叶（L3）为材料，通过整合转录组测序（RNA-seq）与生化分析，系统解析了叶片形态发生过程中的转录-代谢共调控网络。该研究揭示了不同发育阶段叶片转录组和代谢组的显著差异，鉴定了核心代谢物（儿茶素、游离氨基酸、生物碱）的时空积累规律及其关键调控基因，并分析了顺式调控元件的富集特征。研究结果深化了对茶树叶片发育与品质形成协同关系的理解，为优质茶叶栽培提供了理论依据。该论文发表在《Frontiers in Plant Science》。
关键技术方法：本研究采用广东省英德县种植的茶树品种LTDC的芽、第一叶、第二叶和第三叶（2024年4月采集）。主要技术方法包括：（1）转录组测序（RNA-seq）：由广州吉迪奥生物技术有限公司完成，使用LTDC基因组进行基因注释，以DESeq2鉴定差异表达基因（DEGs），阈值设为校正p<0.05且|log2倍变化|>2；（2）高效液相色谱（HPLC）检测儿茶素和生物碱：参照GB/T 8313-2018，使用ZORBAX Eclipse C18柱，紫外检测波长278 nm；（3）游离氨基酸检测：采用生理阳离子交换柱（Sykam S433D系统），以茚三酮柱后衍生，紫外检测570和440 nm；（4）实时定量PCR（qRT-PCR）验证：使用LightCycler480系统，以2
-ΔΔCt法分析相对表达量；（5）顺式调控元件分析：利用PlantCARE数据库对功能基因启动子进行顺式元件预测。
研究结果：
3.1 叶片发育表型与RNA-seq主成分分析：通过表型观察和RNA-seq的主成分分析（PCA），发现四个发育阶段（B、L1、L2、L3）的转录组显著分离，PC1和PC2解释总变异的96.8%，芽与成熟叶（L3）差异最大，表明叶片发育伴随广泛的转录组重塑。
3.2 差异表达基因分析：在B vs L1、B vs L2、B vs L3比较中，DEG总数随叶片成熟显著增加，其中B vs L3鉴定出最多DEGs（6,854个）。KEGG富集分析显示，DEGs主要富集在代谢通路，尤其是全球代谢、碳水化合物代谢、脂质代谢和次生代谢物生物合成，且富集基因数随成熟度增加。
3.3 叶片形态发生过程中儿茶素及相关基因的变化：儿茶素总量呈现先升后降趋势，B和L3最低，L1最高。大多数儿茶素单体（C、EC、GC、ECG、CG）在芽中积累最高，随叶片成熟显著下降；而O-甲基化儿茶素（Me-ECG和Me-EGCG）在L3中显著上调。基因表达分析表明，苯丙烷和类黄酮通路的关键基因（如CHS、LAR、ANS、SCPL）在芽中高表达，而DFR、LAR、ANR及OMTs在L3中上调，与甲基化儿茶素积累一致。
3.4 叶片形态发生过程中游离氨基酸及相关基因的变化：37种游离氨基酸中，大多数鲜味相关氨基酸（L-茶氨酸、谷氨酸、谷氨酰胺、丙氨酸等）在芽中显著富集，随叶片成熟下降。茶氨酸生物合成通路基因（GDH、AlaAT、TS、AlaDC）在芽中高表达，而GOGAT、GS、GDH在成熟叶中上调，可能与光呼吸和氮转运有关。
3.5 叶片形态发生过程中生物碱及相关基因的变化：咖啡因、可可碱、茶碱含量在芽中最高，随叶片成熟显著下降，在L2和L3中最低。茶叶碱（theacrine）未检测到。咖啡因生物合成关键基因（ADK、IMPDH、SAMS、TCS、XMT、MXMT）在芽中高表达，与代谢物积累模式一致。
3.6 不同发育阶段茶叶主要代谢物的相关性特征：Pearson相关性分析表明，大多数游离氨基酸（L-The、Glu、Gln、Asp、Ala、GABA）与生物碱（咖啡因、可可碱、茶碱）和儿茶素（C、GC、CG、EC、ECG）呈显著正相关，而与O-甲基化儿茶素相关性弱，表明芽和幼叶中多种代谢通路协同增强，而甲基化儿茶素具有相对独立的调控途径。
3.7 叶片发育阶段关键差异表达基因的表达分析：从转录组中筛选108个与叶片发育相关的DEGs，qRT-PCR验证30个候选基因。结果显示，细胞壁合成与重塑基因（EXPA1、PRP4等）、光合作用基因（ndhL）和生长素信号基因（SAUR23、IAA14）在成熟叶中表达升高；细胞周期基因（CYCD1-1、CYCA2-4）在芽中高表达后下调；转录因子（GRF1、MYB17、bHLH94）主要在芽中高表达，而bHLH72和MYB60在叶片发育中上调。
3.8 不同发育阶段差异功能基因的顺式作用元件分析：对77个功能基因启动子进行顺式元件分析，发现MYB和MYC结合位点高度富集（分别526和314个），光响应元件（Box4、GATA-motif等）和激素/胁迫响应元件（ERE、ARE、ABRE等）也普遍存在。不同基因簇（如细胞壁、光合、激素信号等）具有特异性元件富集模式，表明叶片发育受激素、光照和胁迫信号的协同调控。
总结讨论：本研究通过转录组和生化分析揭示了茶树叶片形态发生的动态调控网络，阐明了关键品质代谢物的时空积累规律及其遗传调控机制。讨论部分指出，总儿茶素含量在芽和幼叶中高，成熟叶中甲基化儿茶素增加，这与前人研究一致；氨基酸和生物碱在芽中高积累，与基因表达模式吻合；相关性分析揭示代谢物间的协同与特异性关系；顺式元件分析表明叶片发育受多种信号整合调控。研究局限性包括：仅使用单一品种LTDC、单季节采样、仅基于相关性分析的多组学整合，缺乏因果验证。未来需结合多品种、多季节、环境处理及转基因功能验证。研究结论：本研究系统阐明了茶树叶片形态发生的转录-代谢协同调控机制，明确了与关键品质化合物相关的发育分配模式及核心基因，丰富了茶树发育生物学和品质形成的理论框架，并对优质高效栽培、品种改良和原料采收期优化具有重要实践价值。
