近日，中国农业科学院茶叶研究所等合作单位在期刊Horticulture Research上发表了一篇题为“Integrated single-nucleus transcriptomic and metabolomic insights into bud-to-leaf development and metabolite synthesis in tea plant”的研究成果。
该研究通过整合单核转录组学（snRNA-Seq）、批量转录组学（RNA-Seq）和代谢组学技术，系统揭示了茶树从芽到叶发育过程中的细胞分化轨迹和代谢物合成调控机制。研究旨在解决芽叶发育与次生代谢物积累的协同调控问题，为茶树品质改良提供理论依据。
百迈客生物为该研究提供了植物广靶代谢组检测服务！
茶树是一种多年生常绿木本植物，其芽叶经过加工可制成茶饮。茶的品质与风味取决于新芽中次生代谢物的合成与积累，包括黄酮类、氨基酸、生物碱和萜类化合物。尽管学界已对次生代谢物的生物合成与调控机制展开深入研究，但协调芽叶发育与代谢的整合机制仍不甚明了。
该研究通过单核RNA测序（snRNA-Seq）、批量RNA测序（RNA-Seq）和代谢组学技术，全面解析了茶树芽叶的发育轨迹与代谢特征。
研究人员对茶树新梢（包括顶芽、第1叶、第2叶和第3叶）进行了snRNA-Seq分析（图1a）。组织切片显示，叶片中存在上表皮（UE）、下表皮（LE）、栅栏组织（PM）、增殖木质部细胞（PC-XY）、韧皮部细胞（PH）和海绵组织（SM）等多种细胞类型（图1b，图S1a-d）。经过数据质控，共获得48,047个高质量细胞核，每个样本保留8,855-12,179个细胞核。降维聚类分析鉴定出17个不同的细胞簇，并根据标记基因将其归类为8种细胞类型：PM、UE、LE、PC-XY、PH、SM、增殖细胞（PC）和保卫细胞（GC）（图1c-d，图S2a-b）。值得注意的是，PM细胞的比例从芽到叶发育过程中逐渐增加，而PC细胞的比例则下降（图1e）。特别关注的第14细胞簇主要分布在芽中（图1f），可进一步分为6个亚簇，其基因富集于蛋白质同源二聚化、胞质过程、多细胞生物发育等关键代谢过程，表明该簇在促进芽细胞发育中起关键作用。
图1 茶树不同发育组织细胞类型的图谱
通过Monocle 2软件包分析，发现了四个关键的分支点，对应不同的细胞分化状态（图2a-b）。在早期发育阶段，芽显示出初始的分化模式，PM和PC细胞表现出典型的变化（图2c）。通过分支特异性分析，鉴定出与PM和PC分化相关的Top 50基因，这些基因富集于苯丙烷生物合成、单萜类合成、光合作用、α-亚麻酸代谢和脂肪酸延伸等通路（图2d）。对分支节点4的进一步分析揭示了显著改变的基因，这些基因富集于叶绿体组织和光合作用等过程。这些结果表明，在叶片发育过程中，PM细胞逐渐适应光合功能，而PC细胞随着组织成熟而减少，反映了叶片发育过程中组织生存和功能所需的适应性变化。
图2 通过伪时间分析分析的分化轨迹和细胞命运决定
植物激素是植物生长发育的关键调节因子。为了阐明茶树芽叶发育过程中植物激素的合成和积累模式，研究人员采用了多组学方法。可视化分析揭示了6种激素相关合成基因的表达模式（图3a-f）。积累分析显示，六种激素在芽叶不同发育阶段的积累量存在差异：IAA、ABA、JA和CK随着组织发育而减少，而GAs则增加（图3g）。这些动态变化表明，不同的激素在协调芽叶发育过程中扮演着不同的角色，例如JA和ABA可能参与早期防御反应，而GAs则促进细胞伸长和组织成熟。
图3 茶芽和茶叶发育过程中激素积累的模式
通过分析snRNA-seq和RNA-seq数据，研究人员可视化了参与细胞分裂和叶片发育的关键基因（如TCPs、GRFs、GIFs、CUCs）的动态表达模式（图4a-b）。转录组测序鉴定出几个在芽叶发育过程中具有动态表达模式的miRNAs。有趣的是，CsmiRNA156s、CsmiRNA164s和CsmiRNA396s的表达水平随着组织发育而增加，而CsmiRNA319s则下降（图4c）。进一步分析确定了关键的miRNA-靶基因对（图4d）。例如，转基因研究表明，CsmiRNA396b通过靶向抑制AtGRF3表达，从而抑制拟南芥的叶片大小和重量（图4e-g）。这一结果归因于miR396通过影响GRF表达作为有丝分裂细胞分裂的负调控因子。这些发现表明，miRNAs通过靶向特定的调控基因，在茶树芽叶发育轨迹中发挥着关键作用。
图4 控制茶叶发育的调节机制
茶树新梢富含次生代谢物，尤其是类黄酮（图5a）。代谢组学分析揭示了类黄酮（包括黄酮醇、儿茶素和花青素）在芽叶发育过程中的动态积累模式（图5b）。大多数类黄酮的含量随着叶片成熟而增加。PM细胞类型的KEGG富集分析显示，其显著富集于类黄酮生物合成途径（图5c）。基因表达可视化证实，主要的类黄酮生物合成基因随着类黄酮的积累持续上调。相比之下，氨基酸（特别是L-茶氨酸）则呈现相反的趋势，在芽中积累较高，在叶片发育过程中下降。这种差异反映了代谢物的功能特化：类黄酮提供防御，而氨基酸有助于早期生长。
图5 黄酮类代谢流可视化
系统发育树分析将茶树TCP转录因子分为两类：I类TCPs（如CsTCP14）促进组织发育，II类TCPs（如CsTCP3）抑制细胞分裂（图6a）。snRNA-Seq和RNA-Seq联合分析显示，CsTCP3和CsTCP14在组织发育过程中呈现相反的表达式模式（图6b-c）。CsTCP3的表达随着组织成熟而增加，主要定位于PM细胞，与类黄酮积累呈正相关；而CsTCP14的表达则在发育过程中下降。转基因植物实验进一步验证了这些结果：CsTCP3过表达（CsTCP3-OE）增强了类黄酮积累，而CsTCP14过表达则抑制了类黄酮积累（图6d-e）。此外，在转基因植物中观察到显著的叶片表型差异：与野生型相比，CsTCP3过表达抑制了叶片大小并减少了叶缘锯齿，而CsTCP14过表达则促进了叶片增大并抑制了锯齿。功能分析表明，抑制CsTCP3表达会导致CsUGT94P1下调并相应减少黄酮醇苷含量。总之，CsTCP3和CsTCP14在调控类黄酮生物合成和叶片发育中扮演拮抗角色。
图6 类黄酮的积累模式及形成机制
该研究在单细胞分辨率下解析了茶树芽叶发育的细胞图谱和代谢网络，明确了CsmiRNA396b、CsTCP3/CsTCP14等核心元件的功能。成果为茶树分子育种和品质调控提供了新靶点，并为木本植物发育生物学研究提供了方法论参考。未来可基于此框架探索环境因子对茶叶品质的影响。
关于百迈客生物：
北京百迈客生物科技有限公司(简称：百迈客生物)成立于2009年，是一家提供基因多组学测序服务和单细胞组学&空间组学仪器设备的高新技术企业。业务主要包括科技服务、智能制造两大业务板块。全球化的业务布局，在中国、欧洲等地区拥有以博士和硕士为主体的研发及服务团队，先后在Cell、Nature、Science等国际期刊上发表文章数千篇，累计影响因子超万分。拥有MGI、GeneMind、 Illumina、 PacBio、Nanopore、AB SCIEX、Waters、BMKMANU、10X等二代测序、三代测序和质谱检测平台，自主创新的百灵实验室全自动生产线、BMKCloud多组学大数据智能交付平台及亚细胞级S系列空间组学产品，为全球科研单位、育种机构、医药公司等提供高品质基因多组学服务和产品。
