《Frontiers in Plant Science》：Transcriptional signatures associated with environmental and developmental variation in tobacco aroma-related pathways across cultivars and regions
编辑推荐：
背景：烟草香气是一种复杂品质性状，受基因型、发育进程和区域环境之间相互作用的塑造，然而将生态环境背景与香气相关代谢联系起来的转录程序仍未被充分理解。方法：研究人员分析了种植于三个代表性区域的十个烟草（Nicotiana tabacum L.）品种的转录组，并在
背景：烟草香气是一种复杂品质性状，受基因型、发育进程和区域环境之间相互作用的塑造，然而将生态环境背景与香气相关代谢联系起来的转录程序仍未被充分理解。方法：研究人员分析了种植于三个代表性区域的十个烟草（Nicotiana tabacum L.）品种的转录组，并在三个发育阶段进行取样。对标准化后的RNA-seq表达量进行了主成分分析、差异表达分析、软聚类、通路富集分析、基因表达变异分析以及蛋白质-蛋白质相互作用网络分析。结果：种植区域是转录组变异的主要来源，其次是发育阶段，而品种效应相对较小。叶片从S1到S3的成熟过程与光合作用及碳同化途径的抑制、胁迫响应、氧化还原调控、细胞壁和膜重塑、氨基酸代谢及次生代谢的诱导相关。香气相关过程，包括类胡萝卜素、苯丙烷类及氧化还原相关途径，表现出依赖于区域的表达模式。强香型品种在区域间表现出比清香型品种更高的转录变异性，提示其调控网络具有更强的环境响应性。网络分析进一步突出了温度与氧化胁迫相关模块中的候选枢纽基因，包括热激相关蛋白A0A1S4A548和过氧化物氧还蛋白Q A0A1S4B900。结论：这些发现为理解环境背景和发育成熟如何与烟草香气相关代谢重编程相联系提供了转录组层面的框架，并鉴定了可用于未来功能验证的候选基因。
烟草香气是一种由基因型、发育进程、环境条件和采后加工共同作用的复杂品质性状。尽管已知香气形成涉及多类代谢物，如碳水化合物、类胡萝卜素衍生物、苯丙烷类化合物和萜类等，但协调环境信号与香气相关代谢的转录调控程序尚不明确。此前研究多集中于香气化合物鉴定或采后处理，而跨品种、跨区域转录组层面的系统研究缺乏。为进一步阐明生态环境与发育成熟如何影响香气相关代谢通路，研究人员选取十个烟草品种，在湖南郴州（CZ）、湖南龙山（LS）、云南玉溪（YX）三个代表性产区，于中心花期（S1）、果柄成熟期（S2）和叶中部成熟期（S3）三个发育阶段采集叶片样品，共构建了380个RNA-seq文库，开展大规模比较转录组分析。研究发现，种植区域是转录组变异的最主要来源，其次为发育阶段，而品种效应相对较弱。叶片从S1到S3的成熟伴随着光合作用、碳同化相关通路的大范围下调，以及胁迫响应、氧化还原调控、细胞壁和膜重塑、氨基酸代谢和次生代谢途径的激活。强香型品种如K326和粤烟96表现出更高的基因表达变异系数（CV），而清香型品种红大则较低，提示强香型品种可能具有更强的转录可塑性和环境响应性。通过网络分析，研究人员在温度响应、氧化应激和激素信号模块中鉴定了多个候选枢纽基因，例如热激相关蛋白A0A1S4A548和过氧化物氧还蛋白Q A0A1S4B900。该研究为理解环境与发育如何协同塑造香气相关代谢重编程提供了转录组框架，并为后续功能验证和品质改良提供了候选靶点。
研究人员采用的主要技术方法包括：从三个产区（CZ、LS、YX）采集十个品种、三个发育阶段的叶片中（其中YX区域S3阶段额外增加了生物学重复），提取RNA后构建Illumina测序文库，利用HISAT2比对到烟草参考基因组，featureCounts进行基因计数；使用DESeq2进行差异表达分析（FDR<0.01，|log2(Fold Change)|>1）；通过Mfuzz软聚类分析动态表达模式；利用clusterProfiler进行KEGG通路富集；以CV衡量基因表达变异；基于STRING数据库构建蛋白质-蛋白质相互作用网络，并用Cytoscape进行可视化。
研究结果部分按原文小标题归纳如下。在“Profiling scheme in tobacco strains, stages, and regions”中，研究人员分析了10个品种、3个地区、3个发育阶段，构建了380个RNA-seq文库，并利用PCA、差异表达分析、Mfuzz聚类、CV分析和PPI网络等策略系统表征转录组变异。在“Growing region is the dominant source of transcriptomic variation”中，主成分分析显示三个地区的样本形成显著分离的簇，表明地区环境是整体转录组变异的主要来源；发育阶段S1与S2相近而S3明显分离，品种间则重叠较多，说明地区效应大于发育效应，品种效应最弱。在“Developmental-stage effects on aroma-related transcriptional programs”中，Mfuzz聚类鉴定出6类发育表达模式，共享发育基因集的功能富集显示下调基因主要富集于光合作用、叶绿体、类囊体膜相关过程，上调基因在LS中富集于核糖体和次生代谢，而在YX中主要与线粒体功能相关，说明发育重编程具有区域依赖性。在“Differential expression analysis reveals maturation-associated metabolic remodeling”中，S3与S1比较共鉴定877至4870个DEGs，上调基因富集于糖脂代谢、鞘脂生物合成、氨基酸代谢、谷胱甘肽代谢及次生代谢物合成等通路，下调基因富集于光合作用、天线蛋白、碳固定、类胡萝卜素生物合成和碳代谢通路，表明成熟时从光合生长转向结构重塑与胁迫适应。在“Transcriptomic variability differs among aroma-type cultivars and regions”中，强香品种K326和粤烟96的CV中位数较高，清香品种红大最低；三个代表性品种的比较显示，K326在郴州S3 vs S1中上调基因富集于蛋白折叠、热响应、渗透胁迫、氧化应激和活性氧等术语，而红大在玉溪主要富集于昼夜节律和氧化还原稳态，说明不同香气类型品种的转录可塑性和响应途径存在差异。在“Candidate hub genes in aroma-related environmental response modules”中，针对温度响应模块，PPI网络显示A0A1S4A548（类叶绿体基质70 kDa热激相关蛋白）和A0A1S3ZX95（热激因子结合蛋白类似蛋白）具有较高连接度；氧化应激模块中，A0A1S4B900（叶绿体过氧化物氧还蛋白Q）和A0A1S4A3V7（硫氧还蛋白依赖性过氧化物氧还蛋白）是候选枢纽基因；激素信号模块中，多个生长素输出载体组分和生长素响应蛋白表达下降，提示生长素信号可能与发育及环境响应相关。
在讨论部分，研究人员指出区域效应可能反映了气候、土壤、栽培条件和微生物群落的综合影响，环境信号通过光合作用、碳分配、激素信号、氧化还原平衡和次生代谢整合，进而影响香气相关代谢潜力。发育阶段方面，S1与S2相近而S3发生重大转变，光合相关基因广泛抑制，代谢、脂质、细胞壁和次生代谢相关基因上调，表明成熟过程中资源从光合作用向结构和防御代谢重分配。强香品种转录变异性较高，可能赋予其更强的表型可塑性和环境响应能力，从而在不同条件下重塑胁迫、氧化还原和次生代谢通路，但这一关系因未直接测定香气化合物而需谨慎解释。PPI分析所鉴定的候选枢纽基因涉及蛋白稳态、氧化还原和激素信号，为后续功能验证提供了起点。研究结论部分总结如下：本研究提供了转录组层面的框架，表明种植区域是烟草叶片转录组变异的主要来源，其次是发育阶段，品种效应相对较弱；叶片成熟（特别是S1到S3）伴随光合、叶绿体和碳固定相关通路广泛下调，同时激活代谢、胁迫响应、氧化还原和次生代谢过程；强香品种表现出更高的转录变异性，提示其具有更强的转录可塑性和环境响应性；功能富集和PPI分析进一步鉴定了与温度响应、氧化应激、氧化还原稳态和激素介导信号相关的候选模块和枢纽基因。尽管未直接定量香气化合物，这些发现突出了参与环境响应型香气相关代谢重塑的候选通路和调控基因，为未来功能验证、代谢组整合以及旨在提高烟叶品质的栽培和育种策略提供了基础。
