干旱、盐碱、高温、低温、病害等逆境胁迫严重制约着全球农作物的产量与品质。据联合国粮食及农业组织统计，全球约90%的耕地受到一种或多种环境胁迫的影响，主要作物产量损失平均超过50%。传统育种方法依赖表型选择和经验积累，周期长、效率低。高通量测序（NGS）技术的成熟，正在彻底改变植物抗逆研究的范式，使研究者能够从基因组结构、转录调控、表观修饰到微生物组互作，全方位解码植物适应逆境的分子密码。本文将系统梳理NGS在植物抗逆研究中的技术原理、应用场景、标准流程及未来趋势。

图1 小麦干旱胁迫响应模式 （图片来源：网，侵删）

一、背景：当"单一组学"遭遇"复杂性状"

植物对逆境的适应是一个涉及多基因、多通路、多层调控的复杂生物学过程。传统研究依赖QTL定位和候选基因克隆，虽成功解析了部分关键基因（如DREB、NAC、WRKY等转录因子家族），但面临三重局限：

分辨率不足：单一环境胁迫往往激活数百至数千个基因，传统方法难以捕捉完整的调控网络；静态快照：基因表达具有时空特异性，单一时间点的数据无法反映动态响应过程；调控盲区：DNA甲基化、组蛋白修饰、非编码RNA等表观遗传机制在胁迫记忆中扮演关键角色，却长期被忽视。

NGS技术的突破彻底改变了这一格局。通过整合基因组学（Genomics）、转录组学（Transcriptomics）、表观基因组学（Epigenomics）和代谢组学（Metabolomics）等，研究者能够构建从"基因型"到"表型"的完整链条，系统揭示植物感知胁迫、传递信号、启动应答、恢复稳态的全过程。

二、技术应用：

1. 抗旱基因挖掘与分子机制解析

干旱是制约农业生产的首要胁迫因子。多组学整合研究已成为解析作物抗旱机制的标准范式。在番茄抗旱研究中，研究者整合转录组与代谢组数据，鉴定出332个核心差异表达基因和553种代谢物，发现氨基酸代谢、类黄酮合成和ABA信号通路是干旱响应的关键节点。表观遗传机制在干旱记忆中同样关键。苹果干旱响应的多组学研究首次绘制了全基因组表观遗传图谱，整合转录组、全基因组亚硫酸氢盐测序（WGBS）和染色质免疫共沉淀测序（ChIP-seq），揭示了组蛋白H3K9me2修饰与DNA甲基化的时空动态规律，为表观遗传育种提供了分子靶点。

2. 抗盐耐碱的离子稳态调控研究

土壤盐渍化影响全球超过9.5亿公顷耕地。植物通过"拒盐"（限制Na⁺吸收）和"耐盐"（区隔化Na⁺、合成渗透调节物质）两种策略应对盐胁迫。转录组测序揭示了盐胁迫下离子转运蛋白（如HKT1;1、SOS1、NHX家族）的表达调控模式。在葡萄中，VaNAC17转录因子通过调节内源茉莉酸（JA）生物合成和清除ROS，显著增强抗旱性和耐盐性。宏基因组研究则发现，根际耐盐微生物可通过产生胞外多糖（EPS）改善土壤结构，协助宿主维持离子稳态。

3. 抗寒防冻的冷驯化机制研究

低温胁迫影响作物分布和产量形成。冷驯化（cold acclimation）是一个涉及膜脂重构、渗透调节、抗氧化防御的复杂过程。转录组研究鉴定了CBF（C-repeat binding factor）冷响应通路的核心基因CBF1/2/3，它们通过激活下游COR（cold-regulated）基因表达，增强细胞膜稳定性。表观基因组分析发现，冷胁迫诱导组蛋白H3K4me3修饰在CBF位点的富集，促进染色质开放和基因转录。在月季中，RcMYB8通过上调ROS清除相关基因和脯氨酸积累，增强对盐和干旱胁迫的耐受性。

4. 重金属胁迫的解毒与修复研究

工业污染导致的土壤重金属积累威胁食品安全。植物通过螯合、区隔化、外排等机制应对重金属胁迫。转录组测序揭示了重金属胁迫下植物螯合肽（PCs）、金属硫蛋白（MTs）及ABC转运蛋白的表达模式。表观遗传研究发现，DNA甲基化水平的动态变化可调节重金属转运蛋白基因的表达活性，影响植物的富集与耐受能力。宏基因组研究则发现，某些根际内生菌可通过生物吸附和生物转化降低重金属生物有效性，为植物-微生物联合修复提供了新思路。

5. 多胁迫交叉耐受与系统获得性抗性

自然界中植物常面临多重胁迫叠加。多组学研究发现，干旱、盐、冷胁迫存在交叉对话（crosstalk）：激素信号整合：ABA是干旱和盐胁迫的核心信号，而JA/ET通路主导生物胁迫响应，两者通过关键转录因子（如MYC2）实现交叉调控； ROS信号枢纽：活性氧（ROS）既是胁迫损伤因子，也是信号分子，通过氧化还原修饰调控下游基因表达；表观遗传记忆：前期胁迫经历可通过DNA甲基化和组蛋白修饰"印记"在基因组中，使后代获得更强的抗逆能力。

三、常见实验步骤：

1. 实验设计与样本采集

时间序列设计：胁迫处理需设置多个时间点，捕捉早期响应（信号感知）和后期适应（稳态重建）基因； 组织特异性：根（吸收器官）、叶（光合作用）、茎（物质运输）分别取样，揭示器官特异性响应； 对照设置：每个时间点设置未处理对照，采用生物学重复（n≥3），确保统计效力； 样本保存：液氮速冻后-80℃保存，避免RNA降解和代谢物转化。

2. 核酸提取与质检

DNA提取：采用CTAB法或提取试剂盒提取高质量基因组DNA。

RNA提取：采用Trizol法或商业化试剂盒提取总RNA。NanoDrop和Bioanalyzer检测完整性和质量，RIN值作为质控指标。

小RNA提取：需富集18-30 nt的小RNA片段，可采用PAGE胶切胶回收或商业化小RNA提取试剂盒。

关键质控：

DNA：OD260/280在1.8-2.0，OD260/230大于2.0，无降解

RNA：OD260/280在1.9-2.0，RIN值≥7.0，28S:18S rRNA条带比例≥1.5

3. 文库构建与测序

全基因组重测序（GWAS）：

DNA片段化至350-550bp；测序深度：自然群体≥10×，亲本及极端池≥30×；

平台选择：Illumina NovaSeq/MGI DNBSEQ，平衡通量与成本。

转录组测序（RNA-seq）：

rRNA去除：植物RNA中rRNA占比>90%，需使用试剂盒进行去除；文库构建：逆转录一链合成、用dUTP二链合成，片段化、接头连接、文库扩增；测序模式：PE150，数据量6-10Gb/样本。

小RNA测序：

Illumina HiSeq/MiSeq，SE50，每个样本至少10-15 million reads。

表观基因组测序：

全基因组甲基化测序（WGBS）：重亚硫酸盐处理将未甲基化C转为U，配对末端测序，一般覆盖度>30×；

CUT&Tag：ConA磁珠捕获细胞/细胞盒，抗体靶向富集目标蛋白，Tn5转座酶切割目标区域，富集靶蛋白结合DNA片段；

ATAC-seq：Tn5转座酶切割开放染色质区域，快速绘制全基因组可及性图谱。

4. 生物信息分析

基因组分析流程：

质控：FastQC评估数据质量，Trimmomatic去除接头和低质量序列；

比对：BWA-mem或Bowtie2比对参考基因组；

变异检测：GATK流程识别SNP、InDel，筛选非同义突变和胁迫相关基因变异；

GWAS：GEMMA、GAPIT等软件进行关联分析，结合表型数据定位抗逆QTL。

转录组分析流程：

比对/定量：STAR/HISAT2比对，featureCounts/StringTie定量；

差异表达：DESeq2/edgeR筛选差异基因（|log2FC|>1, FDR<0.05）；

功能注释：GO富集（生物过程、分子功能）、KEGG通路分析；

调控网络：WGCNA构建共表达网络，识别关键模块和hub基因。

小RNA数据分析：

质控：去除接头、低质量序列、短于18 nt或长于30 nt的序列。

miRNA鉴定：比对miRBase数据库（v22.1）鉴定已知miRNA；预测发夹结构鉴定新miRNA。

差异表达：Student's t-test，筛选阈值p<0.05且|log2FC|≥1。

靶基因预测：psRNATarget/TargetFinder预测靶基因。

表观基因组分析：

甲基化分析：Bismark比对，识别DMR（差异甲基化区域），关联基因表达；

CUT&Tag分析：PEAK峰 calling，Motif分析识别转录因子结合位点，关联靶基因；

多组学整合：整合甲基化、ChIP-seq、RNA-seq数据，构建"表观-转录"调控网络。

高级分析：

机器学习：随机森林、SVM预测抗逆性状，深度学习识别调控元件；

代谢组关联：Spearman相关分析基因表达与代谢物积累，构建"基因-代谢"关联网络。

四、总结：从"基因清单"到"智能设计"的技术跃迁

NGS技术已在植物抗逆研究领域完成从"单一组学"到"多组学整合"、从"静态描述"到"动态预测"的深刻转变。其核心价值在于：全景式解析：一次性捕获基因组变异、转录动态、表观修饰、代谢重塑，构建胁迫响应的完整图景；精准定位：从QTL区间到候选基因，从转录因子到靶标基因，从表观标记到染色质状态，实现多层次精准定位；预测能力：机器学习模型整合多组学数据，可预测基因型-表型关联，指导分子设计育种。

【翌圣生物·植物抗逆研究NGS完整解决方案】