《International Journal of Molecular Sciences》：Adaptive Structural and Transcriptional Responses Contribute to Cold Tolerance Variation in Xinluzhong 61 and Tahe 2 Cotton Cultivars
编辑推荐：
本研究阐明了耐冷棉花品种新陆中61（C61）和冷敏感棉花品种塔河2（C2）之间耐冷性差异的分子、形态学和解剖学机制。研究人员将幼苗置于0 °C冷胁迫下处理12和24小时，随后对子叶、真叶和茎进行比较转录组学（RNA-Seq）和全面的解剖学分析。转录组分析发现，
本研究阐明了耐冷棉花品种新陆中61（C61）和冷敏感棉花品种塔河2（C2）之间耐冷性差异的分子、形态学和解剖学机制。研究人员将幼苗置于0 °C冷胁迫下处理12和24小时，随后对子叶、真叶和茎进行比较转录组学（RNA-Seq）和全面的解剖学分析。转录组分析发现，在12小时时，C2和C61中分别鉴定出8834个和14,664个差异表达基因（differentially expressed genes, DEGs），在24小时时分别鉴定出14,399个和16,791个DEGs。KEGG富集分析显示，12小时时主要涉及苯丙烷生物合成、光合作用以及角质/木栓质/蜡质生物合成，24小时时则转向吞噬体、半胱氨酸和甲硫氨酸代谢。表型和解剖学观察证实，C61发育出密集的茎腺毛（C2缺乏），在24小时胁迫后真叶中维持显著更大的叶片厚度、栅栏组织厚度和栅栏-海绵比，并表现出更厚的茎木质部。这些发现突出了核心适应性性状，并为改善棉花耐冷性提供了有价值的遗传靶点。
**论文解读文章**
**研究背景与意义**
棉花（Gossypium hirsutum L.）是全球重要的经济作物，为纺织工业提供天然纤维，同时在食用油和动物饲料生产中具有重要贡献。棉花起源于热带和亚热带地区，对低温胁迫尤为敏感，尤其在种子萌发和幼苗建立阶段。低温胁迫会引发一系列有害的生理生化反应，导致生长受阻、生物量积累减少、生殖发育受损，最终造成纤维产量和品质显著下降。在中国新疆等主要棉区，早春寒潮（晚霜）频繁干扰出苗，低温（<15 °C）直接损伤光合机构，诱导过量活性氧（reactive oxygen species, ROS）积累，造成脂质、蛋白质和核酸的氧化损伤，进而表现为叶片栅栏/海绵组织厚度减小、维管结构受损等解剖学改变。
尽管已有研究分别探讨了转录组或解剖学响应，但缺乏系统性地将时间序列转录组分析与跨关键幼苗器官（子叶、真叶和茎）的解剖学观察相结合，以建立分子重编程与结构适应之间直接联系的研究。例如，参与次生代谢或结构发育的基因转录激活与表型性状（如细胞壁成分改变或腺毛密度等保护结构）之间的关联尚需明确。为了填补这一空白，研究人员选取了广泛种植于新疆地区、耐冷性不同的两个陆地棉品种——耐冷型‘新陆中61’（C61）和冷敏感型‘塔河2’（C2），通过时间序列比较转录组学（0 °C严酷冷胁迫）并结合子叶、真叶和茎的解剖学分析，旨在揭示C61耐冷性优越的关键基因、生物学过程和结构修饰。研究假设C61的优越耐冷性源于更快速、更强烈地激活特定防御通路，包括抗氧化系统、次生代谢物生物合成以及物理屏障强化相关通路，这些变化应同时反映在其转录组特征和解剖学特征中。该研究综合多层次分析，为棉花耐冷性适应机制提供了全面理解，并为未来改良该重要作物提供了有价值的遗传靶点。该论文发表在《International Journal of Molecular Sciences》。
**关键技术与方法**
研究人员以中国新疆地区广泛种植的陆地棉品种新陆中61（耐冷型，C61）和塔河2（冷敏感型，C2）为材料。主要技术方法包括：① 时间序列冷胁迫处理：将三叶期幼苗置于0 °C恒温生长箱中处理12小时和24小时，同时设置28 °C正常生长对照（仅24小时取样）；② 转录组测序（RNA-Seq）：提取叶片总RNA，构建文库后在Illumina NovaSeq 6000平台进行150 bp双端测序，使用HISAT2比对到陆地棉TM-1参考基因组（v2.1），DESeq2筛选差异表达基因（DEGs，|log2(fold change)| ≥ 1且调整后p值 < 0.05）；③ 功能富集分析：利用clusterProfiler进行Gene Ontology（GO）和KEGG通路富集分析；④ 解剖学分析：对子叶、真叶和茎进行石蜡切片（番红-固绿染色），使用光学显微镜观察并测量叶片厚度、栅栏组织厚度、海绵组织厚度、茎皮层厚度、木质部厚度等参数；⑤ 抗氧化酶活性测定：检测超氧化物歧化酶（SOD）、过氧化物酶（POD）和过氧化氢酶（CAT）活性；⑥ qRT-PCR验证：对10个随机选择的DEGs进行定量验证。
**研究结果**
**2.1. 低温胁迫对两个棉花品种幼苗的影响**
通过表型观察，0 °C处理后C61幼苗保持直立，仅轻微叶片萎蔫，茎结构完整；C2则出现严重倒伏、叶片广泛萎蔫、水渍状病斑和茎韧性下降。直接比较证实C61在冷胁迫下维持显著更好的植株形态和叶片完整性，验证其更强的耐冷性。
**2.2. 冷胁迫下棉花品种C61和C2的转录组分析**
高质量RNA-seq数据（每个样本42.5–58.4百万条clean reads，Q30 > 93%，GC含量约44%）表明数据可靠。C61与C2之间DEG数量随胁迫时间增加：12小时时鉴定出273个DEGs，24小时时增至637个，表明长时间冷胁迫逐渐放大耐冷与敏感品种之间的转录差异。主成分分析（PCA）显示生物学重复紧密聚类，处理时长是主要变异来源。
**2.3. 低温胁迫下两个棉花品种DEGs的GO富集分析**
12小时冷胁迫时，C61中DEGs显著富集于苯丙烷代谢、ROS响应和细胞壁生物合成等生物学过程（BP），细胞组分（CC）包括过氧化物酶体、微体和质体小球，分子功能（MF）以氧化还原酶、超氧化物歧化酶和抗氧化活性为主，反映C61增强的ROS解毒和次生代谢。24小时时，BP主要富集于腺毛分化、韧皮部/木质部组织发生、肌球蛋白丝动力学、硫转运和昆虫防御反应；CC集中于肌动蛋白细胞骨架和收缩装置；MF涉及运动活性、跨膜转运和催化活性，而受体激酶活性和ATP结合呈负富集，表明长时间胁迫下这些功能下调。
**2.4. 低温胁迫下两个棉花品种DEGs的KEGG富集分析**
12小时时，DEGs最显著富集的通路包括苯丙烷生物合成、光合作用、角质/木栓质/蜡质生物合成，以及脂肪酸代谢、萜类骨架生物合成和亚油酸代谢，提示C61通过强化细胞壁屏障、优化能量代谢和激活次生代谢通路增强耐冷性。24小时时，最显著富集的是吞噬体通路，其次是半胱氨酸和甲硫氨酸代谢、萜类骨架生物合成和甘油酯代谢，从12小时的防御和光合转向免疫激活和抗氧化代谢，体现时间依赖性适应策略。
**2.5. 冷胁迫下两个棉花品种的抗氧化酶活性**
SOD活性：C61在12小时显著上升（p < 0.01）并持续至24小时，表现出快速诱导型防御；C2在24小时较对照显著下降，显示SOD介导的抗氧化能力受损。POD活性：两个品种均同步上调，12小时显著增加，24小时达峰值，表明POD是冷胁迫响应的共同关键组分。CAT活性：C61在12小时短暂升高后24小时急剧下降；C2则维持较高本底CAT活性，但12和24小时均进行性下降（p < 0.001）。结果表明C61依赖诱导型抗氧化系统（持续SOD和瞬时CAT），C2更依赖本底CAT活性，揭示不同的抗氧化策略。
**2.6. 低温处理对两个棉花品种腺毛的影响**
GO分析中腺毛分化为24小时最显著富集的BP条目。茎组织表型观察显示，C61茎表面密布显著腺毛，C2茎相对光滑且腺毛显著较少。定量结果证实C61茎腺毛密度显著高于C2，支持腺毛分化有助于棉花耐冷性的结论。
**2.7. 冷胁迫下子叶解剖结构响应**
与对照相比，0 °C处理24小时后，两个品种的子叶厚度、栅栏组织厚度和海绵组织厚度均显著降低，叶脉维管束出现破裂。C61的子叶栅栏组织厚度、海绵组织厚度、叶片厚度、栅栏-海绵比和叶片结构紧密度均高于C2，表明C61具有更高的潜在光合能力和更慢的水分散失速率。
**2.8. 冷胁迫下真叶解剖结构响应**
低温处理后，C61真叶的栅栏组织厚度、海绵组织厚度和叶片厚度分别下降20.11%、42.82%和34.47%，但栅栏-海绵比增加40%；C2相应下降14.73%、20.58%和16.86%，栅栏-海绵比仅增加7.02%。正常条件下C61的栅栏组织厚度、海绵组织厚度和叶片厚度均高于C2，低温处理后C61的叶片厚度和栅栏-海绵比显著高于C2，表明C61具有更强的冷抗性。
**2.9. 茎解剖结构响应**
24小时冷胁迫后，C61茎皮层厚度轻微下降4.23%，而表皮、韧皮部、木质部厚度和髓径分别增加14.97%、21.14%、35.98%和8.96%；C2皮层厚度下降47.17%，韧皮部厚度基本稳定，木质部厚度和髓径分别增加19.77%和12.56%。C61的韧皮部和木质部厚度在冷处理与对照间差异显著，且C61木质部增厚幅度大于C2，表明其更强的耐冷性。
**讨论与结论**
讨论部分指出，低温胁迫严重制约棉花生长与地理分布。C61的优越耐冷性源于协调的时间依赖性代谢重编程和持续的结构适应。12小时时快速激活苯丙烷生物合成、ROS响应和细胞壁修饰通路，形成物理屏障并增强抗氧化能力；24小时时腺毛分化显著富集，C61茎上发育出密集腺毛，形成物理和化学屏障。解剖学上，C61维持更大的叶片厚度和栅栏组织完整性，茎木质部增厚，保障水分和养分运输。这些发现整合了分子、形态和解剖学证据，揭示了C61耐冷性的多层级机制。研究还确认ICE-CBF-COR通路在冷胁迫下被激活。
**结论部分翻译**：本研究表明，棉花耐冷性由整合的分子、形态学和解剖学适应所调控。耐冷型品种新陆中61表现出时间依赖性转录重编程：早期激活苯丙烷生物合成、光合作用和角质/木栓质通路以快速防御，随后诱导吞噬体和氨基酸代谢以长期驯化。与塔河2相比，C61表现出更广泛的DEG调控（12小时：14,664 vs. 8,834；24小时：16,791 vs. 14,399）、密集的茎腺毛、优越的真叶解剖结构和更厚的茎木质部。这些发现阐明了棉花耐冷性的协同机制，并为耐冷棉花品种的分子育种提供了关键遗传靶点。
