短时渗透、温度和化学胁迫下甲藻Prorocentrum cordatum的qPCR内参基因验证

《Microorganisms》：Validation of Reference Genes for qPCR in the Dinoflagellate Prorocentrum cordatum Under Short-Term Osmotic, Thermal, and Chemical Stressors

可靠的qPCR归一化需要在特定实验条件下验证内参基因。研究人员在甲藻Prorocentrum cordatum中评估了8个候选内参基因，该藻分别暴露于短期温度、渗透和铜胁迫下，跨两个数据集。表达稳定性采用geNorm、NormFinder和BestKeeper进行排序，随后进行共识排名聚合。依数据集不同，最稳定的内参基因对为SAM和RPS9（分别编码S-腺苷甲硫氨酸合成酶和核糖体蛋白S9）或SAM和Crn（编码coronin）。相反，编码苹果酸脱氢酶（MDH）、延伸因子2（EF2）和亲环蛋白（Cyp）的基因排名最低。研究人员使用以实验组和培养系列为固定效应的独立线性模型，评估了归一化对三个目标基因——Mei2样基因Mei2_s3929、丝氨酸/苏氨酸蛋白激酶基因Phkg2和α/β水解酶基因Hydr——的影响。在数据集1中，低渗处理仅在最优归一化（RPS9 + SAM）下显著上调Mei2_s3929；在次优归一化（EF2 + MDH）下该效应无统计学显著性。相比之下，数据集2中Mei2_s3929在光照和黑暗对照之间的差异在两种归一化策略下均保持显著。这些结果表明，内参基因的选择可改变微藻qPCR中的统计效力，并影响目标基因表达变化的统计学检测，即使候选基因满足常规稳定性标准亦然。

研究背景方面，定量PCR（qPCR）是靶向基因表达分析的灵敏方法，尤其适用于在限定实验条件下检测有限数量的基因。然而，qPCR测量的可靠性不仅取决于引物特异性和扩增效率，还取决于用于校正RNA数量、RNA质量、逆转录效率和模板输入变异的归一化策略。归一化通常使用一个或多个假定在比较样本间表达相对稳定的内参基因，但所得表达估计的准确性在很大程度上取决于这些基因的选择。常规管家基因未必组成型表达，其转录水平可能随发育状态、组织或细胞类型、环境条件或实验处理而变化，不应未经实验验证即予使用；MIQE（定量实时PCR实验发表最低信息）指南建议报告内参基因选择的依据并证明其适用于研究设计。甲藻是生态上重要的单细胞藻类，其对变化环境条件和各种胁迫的响应影响其生长、生理和发育策略选择。已有若干研究评估了甲藻qPCR归一化的内参基因，候选集包括编码微管蛋白、肌动蛋白、甘油醛-3-磷酸脱氢酶（GAPDH）、18S核糖体RNA、翻译相关蛋白、代谢酶及其他保守细胞蛋白的基因。尽管某些候选基因在多项研究中表现良好，但最优内参基因因物种、处理、暴露持续时间和实验设计而异，这些观察结果反对将任何常规管家基因视为甲藻qPCR的通用内参。本研究在海洋浮游甲藻Prorocentrum cordatum中评估了8个候选内参基因，该藻暴露于短期低渗和高渗条件、低温和高温以及铜处理，特别关注鉴定适用于短期胁迫（即暴露不超过4 h）期间基因表达归一化的内参。短期胁迫实验在内参基因选择上需特别谨慎：一方面，对某些过程而言，这样的暴露时间尚不足以使变化在表达水平上显现；另一方面，急性暴露可触发快速而短暂的转录变化，即使细胞死亡率仍较低。因此，在长期驯化、跨生长阶段或不同环境处理下表现稳定的基因，可能不适用于急性胁迫响应期间的转录水平归一化。从理解自然环境中发生过程的角度看，此类研究也很重要，因为它们为预测甲藻对突变（例如压载水释放到新栖息地或极端天气事件）的细胞响应提供了数据。P. cordatum此前已用于基因表达研究，但现有的内参基因验证数据未能建立适用于所有实验条件的通用归一化策略。可靠的归一化需要结合互补的统计方法评估候选基因稳定性，并验证内参基因选择对目标基因表达的后果。研究人员选择的候选基因代表多样的细胞功能，包括中心代谢、蛋白质合成、核糖体结构和细胞骨架调控。表达稳定性采用三种统计算法评估，并对所得排名进行整合以获得共识评估。最后，通过比较三个目标基因在最稳定和最不稳定候选组合归一化后的表达模式，检验所选内参的生物学相关性。该设计使研究人员不仅能够确定哪些候选基因在测试条件下表现出最大的统计稳定性，还能确定内参基因选择如何影响胁迫相关目标基因表达变化的检测和解释。

研究人员开展了以下研究：使用Prorocentrum cordatum克隆培养物（菌株CCAP1136/16），在f/2培养基（盐度25）中于12 h光:12 h暗周期、100 μmol photons m^?2 s^?1和18 °C条件下培养至指数生长期。实验培养物（250 mL，初始细胞浓度4 × 10⁴ cells mL^?1）生长6–7天。为评估胁迫下的基因表达稳定性和目标基因响应，进行了两个独立实验集（数据集1和数据集2）。数据集1中，样本在标准光照（100 μmol photons m^?2 s^?1）下接受4 h各种胁迫条件：低温处理（low_T，2 °C）、盐度增加10（hyper_OS，添加0.94 M NaCl）、盐度降低10（hypo_OS，添加蒸馏水）和铜毒性（Cu，最终浓度3.13 μM CuSO₄，对应铜含量200 μg）。数据集2用于评估短期热胁迫：高温处理（high_T）在37 °C下孵育，暴露时间缩短至1 h以防止细胞死亡率急剧增加；由于无法在37 °C培养箱内放置光源，high_T处理在黑暗中进行，因此设立了两种1 h暴露的对照条件：光照对照（C_light，18 °C、100 μmol photons m^?2 s^?1）和黑暗对照（C_dark，18 °C黑暗）。孵化后样本离心、上清去除、细胞沉淀液氮速冻并储存于?80 °C。研究结论为：内参基因选择可改变微藻qPCR中的统计效力并影响目标基因表达变化的统计学检测，即使候选基因满足常规稳定性标准；SAM在短期胁迫条件下始终排名最稳定之列，RPS9和SAM（数据集1）以及Crn和SAM（数据集2）为最可靠的内参对，而MDH以及依数据集而定的EF2或Cyp最不稳定；在数据集1中，低渗响应仅在最优归一化下显著，而在数据集2中光照与黑暗对照的差异在两种归一化策略下均显著，表明内参基因选择对统计检测的影响取决于目标基因和实验对比，而非次优归一化下统一的灵敏度丧失。重要意义在于强调甲藻qPCR内参基因应针对每种物种和实验设计进行经验性验证，而非仅基于其作为内参的常规用途进行选择；并首次提出coronin（Crn）作为潜在新型内参基因，但其适用性需在其他物种和根本不同的实验设计中验证。

关键技术方法方面（不超过250字）：研究采用两个独立实验数据集（数据集1含对照及低温、高渗、低渗、铜处理，数据集2含高温及光/暗对照），以P. cordatum CCAP1136/16菌株为材料；总RNA经裂解缓冲液和珠磨法提取、DNase I处理，逆转录合成cDNA；针对8个候选内参基因和3个目标基因设计特异性引物（MDH引物沿用Guo和Ki），经扩增效率标准曲线（E = (10^?1/slope ? 1) × 100%）和熔解曲线验证；qPCR采用SYBR mix在CFX96系统上运行，样本最大化策略、技术三重复并设无模板对照（NTC）和无逆转录酶对照（no-RT）；表达稳定性以geNorm、NormFinder和BestKeeper评估，并以蒙特卡罗交叉熵（CE）排名聚合获得共识排名，辅以留二法刀切分析验证Vn/n+1；目标基因表达以最优和次优内参对归一化，采用含Group和Series固定效应的线性模型及emmeans估计边际均值（EMMs）进行比较（数据集1用Dunnett校正，数据集2用Tukey校正），并以偏F检验评估Series贡献。

研究结果方面：3.1 细胞培养生长与死亡率测试——指数期细胞（平均生长率0.114 ± 0.006 day^?1）暴露于各种短期胁迫，死亡率范围为“Cu”处理的3.0% ± 0.02%至“hypo_OS”处理的6.6% ± 0.03%；37 °C孵育4 h导致绝大多数细胞死亡，故数据集2处理时间缩短至1 h（死亡率5.8% ± 0.04%）。3.2 候选内参基因的表达谱——数据集1中候选基因中度至高度表达，平均Cq值20.27–24.79，SAM表达最高（平均20.27 ± 0.79），MDH最低（平均24.79 ± 0.88）；RPS9的Cq范围最窄（2.28个循环，SD = 0.65），GAPDH变异最宽（3.63个循环，SD = 0.87）。数据集2中平均Cq值20.14–24.70，SAM最高（平均20.14 ± 0.64），Cyp最低（平均24.70 ± 0.47）；RPS9范围最窄（1.76个循环，SD = 0.51），MDH变异最宽（3.80个循环，SD = 1.14）。总体Cq变异处于qPCR归一化的合适范围。3.3 候选内参基因的表达稳定性——geNorm分析中，数据集1的RPS9和SAM的M值最低（0.211），EF2最高（0.376）；数据集2的GAPDH和SAM的M值最低（0.191），MDH最高（0.417）；所有候选基因M值均低于常用阈值1.5。BestKeeper分析中，数据集1的RPS9的BestKeeper SD最低（0.516）、CV最低（2.097%），GAPDH最高（0.752和3.459%）；数据集2的Cyp最低（0.373和1.507%），MDH最高（0.958和3.954%）；两个数据集中所有8个候选基因均满足BestKeeper SD ≤ 1.0的标准。NormFinder分析中，数据集1的SAM稳定性值最低（SV = 0.04），EF2和MDH最高（各0.19）；数据集2的Crn最稳定（SV = 0.05），MDH最高（0.33）。共识排名（蒙特卡罗交叉熵法）显示，数据集1中RPS9居首，依次为SAM、Crn、Cyp、RPL7、GAPDH、MDH、EF2；数据集2中Crn居首，依次为SAM、RPS9、GAPDH、EF2、RPL7、Cyp、MDH；表明单个内参基因的相对适用性取决于各数据集所代表的实验条件。geNorm配对变异分析中，完整数据集的V2/3值数据集1为0.073、数据集2为0.082，均低于阈值0.15；留二法刀切分析中数据集1的120个14/16子样本平均V2/3为0.073（95% CI：0.061–0.079，最大0.079），数据集2的66个10/12子样本平均V2/3为0.081（95% CI：0.049–0.090，最大0.090），所有Vn/n+1值均低于0.15，支持每个数据集使用两个内参基因归一化：数据集1为RPS9和SAM，数据集2为Crn和SAM。3.4 通过目标基因表达谱验证内参基因——数据集1中Mei2_s3929对低渗胁迫的响应取决于归一化策略：最优归一化下hypo_OS与对照差异显著（对比估计 = ?0.689 ΔCq单位，Dunnett校正p = 0.035，约1.61倍增加）；次优归一化下估计值亦为负（?0.574 ΔCq单位）但无统计学显著性（校正p = 0.250）；铜、低温暴露或高渗胁迫下无论何种归一化策略均无显著变化；Phkg2和Hydr在任一条件下两种内参组合均无显著差异。数据集2中Mei2_s3929在光照对照（C_light）与黑暗对照（C_dark）之间在两种归一化策略下均显著不同：最优归一化下对比估计为?0.581 ΔCq单位（Tukey校正p = 0.011，C_light相对表达约高1.50倍）；次优归一化下差异仍显著（估计 = ?0.404 ΔCq单位，校正p = 0.019，约高1.32倍）；high_T条件在两种归一化策略下与C_light均无显著差异（最优：校正p = 0.170；次优：校正p = 0.294），C_dark与high_T的直接比较亦无统计学显著性（最优：校正p = 0.128；次优：校正p = 0.141）。偏F检验显示，数据集1最优归一化下Mei2_s3929的Series总体效应显著（F_4,7 = 9.12，p = 0.0066），数据集2在最优和次优归一化下均显著（分别为F_3,6 = 17.66，p = 0.0022和F_3,6 = 19.06，p = 0.0018），经Holm校正后仍显著；Phkg2和Hydr在数据集1中未检测到显著的Series效应。

讨论部分总结：准确的qPCR目标基因表达评估根本上取决于数据归一化所用稳定表达内参基因的选择。研究人员在暴露于短期渗透、温度和化学胁迫的P. cordatum中评估了8个候选内参基因的表现。对甲藻qPCR内参基因验证的既往研究综述显示，已发表工作中涉及的基因集部分重叠，最常评估的候选包括经典真核管家基因（微管蛋白、肌动蛋白、GAPDH和18S核糖体RNA），以及其他参与翻译和基础细胞代谢的蛋白编码基因（eIF4E、SAM、MDH和Cyp）；这些基因的适用性似乎强烈依赖背景。因此，不同物种和处理采用了不同的内参组合，意味着内参基因稳定性不能可靠地跨实验条件和分类群外推。研究人员选择的候选池包括GAPDH、MDH、Cyp、SAM、EF2以及编码大核糖体亚基蛋白L7（RPL7）、小核糖体亚基蛋白S9（RPS9）和coronin（Crn）的基因；与18S rRNA基因不同，核糖体蛋白基因在甲藻类似研究中很少使用，但已在其他藻类（包括绿藻和多种异鞭毛藻）中验证为内参，因此扩展候选集以增加所测内参的功能多样性。使用参与不同细胞内过程和代谢途径的蛋白编码基因可降低共调控的可能性，从而可能提高表达分析的准确性；还选择了编码coronin（一种参与肌动蛋白动力学、囊泡运输和胞质分裂调控的肌动蛋白结合蛋白）的基因，据研究人员所知，尚无其他研究将coronin定义为藻类或其他生物qPCR归一化的标准或经验证内参基因。三种广泛接受的方法（geNorm、NormFinder和BestKeeper）评估显示，数据集1中geNorm和NormFinder排名高度一致而BestKeeper结果明显不同；数据集2中geNorm和NormFinder均将SAM列为两个最稳定候选之一，但另一高排名基因在算法间不同，BestKeeper则将Cyp和RPS9列为最稳定。此类排名差异在包括甲藻在内的多种生物中已有报道，并不意外，因为三种算法使用不同的数学框架；然而SAM在两个数据集中始终位列最稳定候选之列，支持其作为P. cordatum短期胁迫条件下内参基因的适用性。共识排名确定数据集1的RPS9和SAM、数据集2的Crn和SAM为最可靠内参对，而MDH以及依数据集而定的EF2或Cyp最不稳定。值得注意的是，SAM、RPS9和Crn表现出不同的表达水平（SAM平均Cq在数据集1和2中分别为20.27和20.13，而RPS9为24.58、Crn为24.11），使用这些基因对归一化可减少丰度特异性技术变异的影响。与P. cordatum唯一既往研究（Guo和Ki，检测温度和化学胁迫，暴露12、24或48 h，确定TUA和MDH最稳定、其次GAPDH）比较凸显了准确选择和验证候选内参基因的重要性：即使在单一微藻物种内，潜在内参基因的选择也非通用，支持甲藻qPCR内参基因应针对每种物种和实验设计进行经验性验证的观点。除实验设计不同外，方法与Guo和Ki亦有差异：后者仅依赖geNorm且未通过目标基因表达分析验证所选基因，而研究人员采用三种算法并整合结果，还通过以最高或最低排名候选归一化目标基因评估内参基因选择的实际后果。本实验限于指数期收获细胞的短期处理，覆盖的生理范围较窄；所有测试内参基因均表现出相当高的表达稳定性，均满足所应用的稳定性标准。尽管如此，数据集1中Mei2_s3929的低渗响应在最优内参对（RPS9 + SAM）归一化时具有统计学显著性，而次优对（EF2 + MDH）则不然；相比之下，数据集2中C_light–C_dark差异在两种归一化策略下均保持显著。因此，内参基因选择对统计检测的影响取决于目标基因和实验对比，而非次优归一化下统一的灵敏度丧失，这与显示次优内参基因可改变组间差异是否达到统计学显著性的qPCR研究一致。在所有目标基因中，仅Mei2_s3929对所检测的实验条件表现出敏感性；在甲藻中，MEI2样蛋白似乎与生活史有性阶段的转变以及胁迫诱导的囊胞形成（包括冷诱导临时囊胞形成）相关，但本研究数据未显示低温条件下Mei2_s3929表达的变化，相反其在数据集1低渗条件下上调、在数据集2黑暗对照（C_dark）下下调，而在结合高温（37 °C）和黑暗的high_T处理中未检测到显著变化。值得注意的是，S. acuminata的 Pellicle 囊胞形成在冷和黑暗联合条件下被诱导，由于这些胁迫同时施加，其各自对基因表达变化的贡献无法区分，因此S. acuminata中某些Mei2样基因表达的改变可能反映对黑暗而非低温的响应；这些观察凸显了环境条件与Mei2样基因表达关系的复杂性，并提示MEI2样蛋白可能具有超出生活史转变调控的更广泛功能。Phkg2和Hydr在所检测实验条件下均未显示表达的显著变化；研究人员未基于本实验设计得出关于其表达稳定性或细胞功能的结论。本研究主要目的是证明使用最优经验证内参基因进行qPCR归一化的重要性，并更广泛地强调在候选内参基因将要使用的特定实验条件下对其进行彻底验证的关键需求。研究结论（翻译）：这些结果表明，内参基因的选择可改变微藻qPCR中的统计效力，并影响目标基因表达变化的统计学检测，即使候选基因满足常规稳定性标准亦然。