菜豆热响应开花相关编码与非编码RNA调控网络（Regulatory Networks of coding and non-coding RNAs involved in Heat-Responsive Flowering in Common Beans）的鉴定

《Journal of Plant Growth Regulation》：Identification of the Regulatory Networks of coding and non-coding RNAs involved in Heat-Responsive Flowering in Common Beans

　　

高温胁迫对豆类生产构成关键威胁，主要通过诱导生殖失败产生影响。为阐明与耐热性相关的复杂调控网络，研究人员对耐热菜豆（Phaseolus vulgaris）品种Perola及热敏感品种Karaca?ehir-90从营养生长至结荚期进行高温处理（32°C/27°C），开展了整合多组学分析（mRNA、miRNA、lncRNA、circRNA）。研究结果表明，Perola的耐热性并非由胁迫响应的强度驱动，而是由独特的"转录韧性（transcriptional resilience）"机制介导——该策略优先维持基因表达的精确稳态，有效缓冲生殖器官免受敏感基因型中观察到的系统性失调影响。相反，Karaca?ehir-90表现出"失调的转录重编程（dysregulated transcriptional reprogramming）"，其特征为花芽中超过2600个基因显著失调；此种广泛重编程与关键开花调节因子SPL14被基因型特异性"miRNA爆发"（如novel_42）抑制及lncRNA TCONS_00080848过度诱导呈负相关，而后者与抗氧化性查尔酮合酶（Chalcone Synthase, CHS）基因抑制呈负相关。Perola规避了这些适应不良的调控级联反应，转而部署"代谢缓冲（metabolic buffering）"策略，其标志为Albumin-2和LTP1提前＞400倍诱导以保护生殖器官。此外，整合竞争性内源RNA（competing endogenous RNA, ceRNA）网络分析显示，敏感品种退化为指示膜损伤的脂质分解代谢，而耐受品种主动维持糖酵解通量及生物能稳定性。这些结果重新定义了耐热性为预防复杂调控改变并优先进行代谢保护的能力，并鉴定出特定的ncRNA-mRNA互作（尤其是novel_42/SPL14及Albumin-2）可作为培育气候适应性作物的潜在靶标。

论文解读：菜豆热响应开花相关编码与非编码RNA调控网络的整合多组学解析

该论文发表于《Journal of Plant Growth Regulation》。普通菜豆（Phaseolus vulgaris）是全球重要的蛋白质与微量营养素来源，高温胁迫（昼温＞30°C、夜温＞22°C）通过降低花粉活力、抑制花粉管生长及减少结荚导致严重减产，热带适耕区预计2050年因变暖失去栽培适宜性。已有转录组研究聚焦蛋白编码基因，而非编码RNA（non-coding RNAs, ncRNAs）——包括microRNA（miRNA）、长链非编码RNA（long non-coding RNA, lncRNA）及环状RNA（circular RNA, circRNA）——在菜豆热胁迫开花调控中作用不明；植物中竞争性内源RNA（competing endogenous RNA, ceRNA）假说框架下lncRNA/circRNA作为miRNA海绵（sponge）调控靶mRNA的机制在菜豆耐热开花中尚待揭示。为此，研究人员选用耐热品种Perola与热敏品种Karaca?ehir-90，设置对照（25°C/20°C）与高温（32°C昼/27°C夜）处理，取茎尖分生组织（shoot apical meristem, SAM）、花芽（flower bud）及开放花（open flower）进行整合转录组（mRNA、lncRNA、circRNA、miRNA）分析、ceRNA网络构建及qRT-PCR验证，结合生理生化指标（相对含水量RWC、丙二醛MDA、脯氨酸Pro、开花数等），系统比较两基因型热响应差异，重新定义菜豆耐热性的分子内涵并鉴定关键ncRNA-mRNA调控模块。

主要关键技术方法

研究人员选用菜豆耐热品种Perola与热敏品种Karaca?ehir-90（土耳其种子基因库来源），设置对照（25°C/20°C）与高温胁迫（32°C昼/27°C夜，渐变升温持续14天覆盖V4至R2期）处理，采集茎尖分生组织（BBCH 15）、花芽（BBCH 55）及初开花朵（BBCH 64）三种组织（n=3生物学重复）。测定生理指标（RWC、SPAD、脯氨酸、MDA、开花数等）。总 RNA经rRNA去除建库，进行链特异性mRNA/lncRNA/circRNA文库（Illumina NovaSeq PE150）及小RNA文库（SE50）测序；lncRNA通过CPC2、CNCI及Pfam三联算法筛选（≥200 nt、≥2外显子）；circRNA用CIRI2与find_circ双流程鉴定（≥2独特背剪接连接reads）；miRNA用miRDeep2与miREvo预测。差异表达分析采用DESeq2（|log₂FC|≥1.3，FDR＜0.05；ceRNA构建用|log₂FC|≥2.0），miRNA靶基因用psRNATarget与TargetFinder预测并经Pearson反相关过滤（|r|＞0.9，p＜0.01）。GO与KEGG富集用clusterProfiler，ceRNA网络用Cytoscape与NetworkX可视化。qRT-PCR用Actin-11（mRNA/lncRNA/circRNA）与U6（miRNA）归一化，2^?ΔΔCt法验证。

Physiological and Biochemical Responses to Heat Stress

Perola在高温下保持较高RWC（花组织81.2% vs 敏感型72.4%）、株高无明显受抑、花组织温度较低（28.3°C vs 31.5°C）；脯氨酸积累更显著（SAM中近3倍），MDA水平仅为敏感型约一半，开花率维持约80%（敏感型下降~60%）。结论：Perola耐热性体现为优越的水分与氧化稳态及生殖发育维持。

Global Transcriptomic Landscape of Coding and Non-Coding RNAs Under Heat Stress

36个文库经质控比对到菜豆参考基因组（v2.1），生物学重复相关性R＞0.96。鉴定得37316个mRNA、3660条高置信lncRNA（多为正义重叠、较短、少外显子）、1144个新circRNA（主要外显子来源，染色体3/9/7富集）及189个miRNA（10保守+179新）。结论：数据集质量可靠，提供菜豆热胁迫下编码与非编码RNA动态全景。

Identification and Functional Annotation of DE-mRNAs, DE-lncRNAs, DE-miRNAs, and DE-circRNAs

Perola花芽仅884个差异表达基因（differentially expressed genes, DEGs），Karaca?ehir-90花芽达2624个（大规模失调转录重编程）；Perola花芽中Albumin-2（Phvul.007G276300）与LTP1较敏感型分别＞400倍与显著上调，属代谢缓冲策略。lncRNA在Perola各组织波动极小（花芽83个DE-lncRNA），Karaca?ehir-90花芽达237个。miRNA在Karaca?ehir-90花芽爆发性差异表达（149个DE-miRNA），Perola仅29个；novel_42（靶向SPL14）在敏感型维持/升高抑制开花，Perola下调释放SPL14抑制；TCONS_00080848（lncRNA）在敏感型花芽过表达且与CHS抑制负相关。KEGG显示Perola花芽富集糖酵解/戊糖代谢维持生物能，Karaca?ehir-90倾向脂肪酸代谢（膜损伤迹象）。结论：耐热性源于转录稳态维持与代谢缓冲，非广谱强诱导胁迫响应。

Construction and Analysis of ceRNA Regulatory Network

整合得908个miRNA-lncRNA-mRNA、591个lncRNA-mRNA共表达对及38个miRNA-circRNA互作。严格过滤后构建组织特异性ceRNA网络：Perola生殖组织聚为一类，Karaca?ehir-90无协调模式；关键模块含novel_42靶向SPL14、TCONS_00080848与CHS互作、pvu-miR2119靶向丙酮酸激酶（Phvul.001G064000）。结论：ceRNA网络揭示ncRNA通过海绵吸附或协同调控决定开花相关基因在热胁迫下的命运，基因型间网络拓扑截然不同。

Enriched Pathway-Specific Marker Gene Expression and ncRNA Regulatory Dynamics

途径标记基因热图显示Karaca?ehir-90花芽CHS、DFR、ANS显著抑制（黄酮类抗氧化能力丧失），FT被抑制、ELF4上调（开花受阻）；Perola维持CHS稳定、诱导特定FT同源物。lncRNA TCONS_00080848与CHS强烈相关。结论：ncRNA（尤其lncRNA）参与微调黄酮类合成与开花信号，敏感型丢失此保护。

qRT-PCR Validation of DEmRNAs, DEmiRNAs, DElncRNAs, and DEcircRNAs

21个候选转录本（8 mRNA、7 lncRNA、4 miRNA、2 circRNA）qRT-PCR与RNA-seq趋势吻合度＞80%，Albumin-2在Perola花芽大幅诱导及novel_42/SPL14模块差异获确认。结论：转录组与ceRNA网络可靠。

讨论与结论翻译

讨论指出Perola通过维持转录稳定性与代谢稳态（增强糖酵解通量、黄酮类合成、渗透保护）而非放大转录扰动实现耐热与 timely flowering；Karaca?ehir-90广泛转录重编程伴能量代谢转向脂质分解、CHS抑制及miRNA爆发（novel_42压制SPL14），导致生殖失败。ceRNA网络阐明lncRNA（如TCONS_00080848）、circRNA及miRNA（novel_42、pvu-miR2119）协同调控开花与代谢途径，为气候适应性菜豆分子育种提供靶标（如编辑novel_42或过表达Albumin-2）。

结论：本研究通过对耐热与热敏菜豆品种整合编码与非编码RNA多组学分析表明，Perola在高温胁迫下通过协调能量代谢与抗氧化途径维持开花；Karaca?ehir-90表现为代谢失衡、miRNA活性增强（如pvu-miR2119）及偏向脂质分解代谢，与生殖性能下降相关。整合ceRNA网络提示miRNA与lncRNA参与调控上述响应，所鉴定的调控模块（TCONS_00080848/CHS、novel_42/SPL14、pvu-miR2119/丙酮酸激酶）为菜豆耐热机制研究与后续功能验证提供基础。