《Journal of Food Measurement and Characterization》：Contrasting oxygen and carbon dioxide contributions to respiratory quotient variability in fresh and fresh-cut produce under modified atmospheres
编辑推荐：
呼吸商（RQ，定义为二氧化碳产生速率RCO?与氧气消耗速率RO?之比）整合了二氧化碳产生与氧气消耗，但其比值结构在氧气消耗较低时会放大变异性。本研究量化了处理依赖的RQ动态，并比较了在气调条件下新鲜及鲜切农产品中氧气与二氧化碳对RQ变异性的贡献。研究人员采用广
呼吸商（RQ，定义为二氧化碳产生速率RCO?与氧气消耗速率RO?之比）整合了二氧化碳产生与氧气消耗，但其比值结构在氧气消耗较低时会放大变异性。本研究量化了处理依赖的RQ动态，并比较了在气调条件下新鲜及鲜切农产品中氧气与二氧化碳对RQ变异性的贡献。研究人员采用广义估计方程（GEE）分析了生菜、菠菜、樱桃番茄和混合沙拉在四种气体处理和三种切割等级下的时间序列数据，并以经验二阶响应面支撑蒙特卡洛模拟、龙卷风分析、偏秩相关系数（PRCC）及分块半偏R2分析。GEE识别出空气处理、菠菜、樱桃番茄、选定气体×切割交互项及选定时间依赖项的显著系数，其中Mix C×细切交互项呈现最大正向交互系数（β=0.4378，p=0.042）。在基线过滤分析中，CO?的龙卷风跨度超过O?跨度（3.258对0.733），缩尾处理后结果相似（3.367对0.743）；但在限定高RO?子集时排名反转（O?：1.136；CO?：0.797）。缩尾PRCC显示CO?呈名义显著但极弱的负关联（?0.0708，p=0.0148），O?无实质关联（0.0006，p=0.9824）。移除O?相关预测变量块比移除CO?相关块导致更大的模型拟合下降（0.0677对0.0377）。这些对比排名表明大气对RQ的影响既依赖度量指标也依赖响应域。因此RQ最好作为随处理变化的轨迹与RO?和RCO?并列解释，而非作为发酵的独立生化标记。
研究背景方面，呼吸作用是决定新鲜及鲜切农产品采后行为的核心因素，它调控底物消耗、气体交换、代谢稳定性及胁迫响应发育。氧气可用性、二氧化碳积累、商品生理、组织损伤与贮藏时长在气调条件下共同决定呼吸行为。传统上呼吸商（RQ，Respiratory Quotient，即RCO?/RO?）被视为呼吸平衡的综合指标，接近1代表以有氧糖呼吸为主，升高可能提示氧限制、底物转换、代谢胁迫或发酵贡献增加，但RQ本身不能单独确认发酵，需乙醇或乙醛等代谢物佐证。其统计解释困难在于RQ由两项生理通量作比，当分母RO?趋近低位时，微小变化或测量误差会被比值结构放大，产生异方差与重尾响应，致使跨处理均值比较不足以刻画时间行为与大气敏感性。已有改良大气研究多建模RO?与RCO?，少有将RQ本身作为首要响应变量，亦缺乏将重复测量推断与扰动、条件关联、分块拟合贡献三类敏感性测度整合的框架。研究人员因此在《Journal of Food Measurement and Characterization》发表该研究，旨在把RQ作为派生生理响应，区分实验处理效应、大气扰动下模型响应幅值、大气变量条件贡献三个维度。
关键技术方法上，样本来自本地商业供应商的罗马生菜、菠菜、樱桃番茄（Solanum lycopersicum var. cerasiforme）及混合沙拉（40%生菜、40%菠菜、20%樱桃番茄），每样约30 g，设三种切割等级（整叶/粗切3–6 cm/细切0.5–1.0 cm，番茄对应整果/半切/薄片），四种气体（空气、Mix A：3% O?+5% CO?、Mix B：4% O?+10% CO?、Mix C：5% O?+3% CO?，余为N?），4℃下按0–168 h多点测量PO?、PCO?、RO?、RCO?。过滤RO?<0.01 mL kg?1 h?1观测后剩1185条。方法链含：广义估计方程（GEE，Gaussian族恒等链接、独立工作相关、稳健三明治标准误）估处理效应；普通最小二乘拟合二阶响应面（RQ~PO?+PCO?+PO?2+PCO?2+PO?×PCO?）作替代模型；对PO?/PCO?拟合计Johnson SB边际分布并抽10000次做蒙特卡洛；5–95百分位单变量扰动得龙卷风跨度；对观测数据做偏秩相关（PRCC）；分块删项算半偏R2。辅以缩尾与高RO?子集（>0.1 mL kg?1 h?1）做稳健性检验。
研究结果部分， Temporal evolution of RQ 小节中，GEE显示空气处理相对Mix A正向（β=0.4098），菠菜相对生菜正向（β=0.3237），樱桃番茄负向（β=?0.3094）；显著气体×切割交互含Air×粗切（β=0.2494）与Mix C×细切（β=0.4378）；时间交互项表明RQ轨迹非单一时间主效应，而随气氛×切割组合变化，细切组织在低氧高CO?下RQ早期突升且与RO?并行下降，符合分母放大与可能氧限制。Monte Carlo and Tornado sensitivity analysis 小节中，基线过滤下中位预测RQ=0.800，O?龙卷风跨度0.733、CO?跨度3.258；缩尾后0.743对3.367，CO?主导排名在轻度缩尾下稳定；蒙特卡洛有13.43%输出为负、1.16%超RQ=5，仅作替代模型不确定性披露。Effect of higher-RO? subset selection on surrogate-model performance and sensitivity ranking 小节中，高RO?子集调整R2由0.0999升至0.4168、SEE由0.536降至0.248，O?跨度反转为1.136、CO?为0.797，证明敏感性排名依赖保留的响应域。PRCC and conditional monotonic association 小节中，缩尾PRCC为CO? ?0.0708（p=0.0148）名义显著但极弱、O? 0.0006（p=0.9824）无实质关联，与龙卷风不矛盾因二者问答不同。Integrated sensitivity interpretation 小节中，删O?块致R2降0.0677、删CO?块降0.0377，即CO?在百分位扰动幅度更大、O?块对替代模型拟合增量贡献更大，多指标排名互补而非冲突。Implications for respiratory quotient characterization and uncertainty assessment 小节指出RQ应作处理依赖轨迹与RO?、RCO?并列，RQ升+RO?降可作实验条件下呼吸失衡的操作预警，但不能代发酵生化证据。Relationship with previous characterization and modelling studies 小节将其与Templalexis与Xanthopoulos基于Michaelis–Menten的CO?抑制O?消耗敏感性框架对照，说明从单通量拓展到比值变量后敏感性维度更复杂。Methodological implications and future work 小节承认替代模型非机制、未测发酵代谢物、蒙特卡洛用独立边际分布未保联合依赖，未来可联动力学质量平衡与温度动态。
讨论与结论翻译部分，研究人员在Conclusions中总结：GEE证实RQ受气体处理、商品类型、选定气体×切割交互及处理特异时间效应影响，Mix C×细切正向交互表明大气组成与组织破损共同左右呼吸平衡。假设“RQ提供超越任一单通量信息”获部分支持——它整合RCO?与RO?变化并凸显呼吸失衡组合，但低RO?时比值结构放大变异性，故RQ应补充而非替代RO?与RCO?单独报告；RQ升伴RO?降在细切组织受限大气中可作操作预警但非发酵直接证据。通过整合重复测量推断与互补敏感性测度，研究证明RQ不能归为单一域不变的大气驱动：基线及缩尾下CO?产生更大百分位替代响应跨度，但删O?块致更大拟合下降，PRCC实用幅值可忽略，高RO?子集中龙卷风排名反转。因而大气重要性取决于分析目标是处理响应、扰动幅值、条件关联还是分块拟合贡献，以及保留的分母域。方法论进展在于提供度量与域感知的RQ解释框架，RQ宜作与时间分辨的筛查信号与RO?、RCO?并列；基线替代模型解释力有限，结论支持比较性、数据集条件性解释而非高精度域不变预测。
