《Life》：Genomic Insights into Carotenoid Biosynthesis and Other Putative Functional Traits of Hot Spring-Derived Paracoccus bogoriensis PH1
编辑推荐：
本研究对一株分离自泰国北部温泉的产类胡萝卜素Paracoccus菌株PH1进行了基因组序列表征。功能分析聚焦于应激耐受、类胡萝卜素生产、微量营养素生物合成、益生特性、生物被膜形成及安全性特征。平均核苷酸一致性（Average Nucleotide Identi
本研究对一株分离自泰国北部温泉的产类胡萝卜素Paracoccus菌株PH1进行了基因组序列表征。功能分析聚焦于应激耐受、类胡萝卜素生产、微量营养素生物合成、益生特性、生物被膜形成及安全性特征。平均核苷酸一致性（Average Nucleotide Identity, ANI）与基于16S rRNA基因的系统发育分析证实了其分类学地位为Paracoccus bogoriensis。基因组注释鉴定出四氢嘧啶（ectoine）和类胡萝卜素的生物合成基因簇，以及应激耐受基因，与此前报道的表型特征相吻合。此外，基因组挖掘揭示了B族维生素、短链脂肪酸（short-chain fatty acids, SCFAs）和吡咯喹啉醌（pyrroloquinoline quinone, PQQ）生物合成的推定功能，以及与益生潜力相关的决定因子，包括耐酸耐胆盐、生物被膜形成和细胞黏附。未检测到抗生素耐药决定因子或毒力相关基因。这些基因组学发现凸显了该温泉来源Paracoccus菌株的多方面生物技术潜力，有必要在相关应用情境中进行进一步的功能评价。
一、研究背景与研究目的
类胡萝卜素是一类具有功能性和商业价值的生物活性色素，呈现黄、橙、深红等光谱色彩，兼具抗氧化、抗癌、抗炎、免疫调节以及保护皮肤和心血管系统等生物活性。全球类胡萝卜素市场预计从2024年的21亿美元增长至2029年的29亿美元，年复合增长率（compound annual growth rate, CAGR）为6.7%。细菌因生长迅速、培养不受季节限制，被视为天然类胡萝卜素最有前景的来源之一，且不与人争粮，契合可持续发展目标中关于粮食安全和可持续农业的要求。温泉来源的耐热菌株因其独特的生长条件在生物技术应用中备受关注。前期研究已报道温泉来源的Paracoccus菌株PH1是多种类胡萝卜素的天然生产者，主要产角黄素（canthaxanthin）[26]，但该产色素表型及其他潜在功能性状的基因组基础尚未被表征。为此，研究人员通过全基因组测序与挖掘，系统探究了该菌株的应激耐受能力、类胡萝卜素生产能力、微量营养素合成潜力、益生特性以及基于毒力因子与抗生素耐药性的安全性特征，以期为该菌株在生物技术和水产养殖等领域的应用开发提供基因组学依据。
二、主要技术方法
研究人员采用Illumina NextSeq 550平台对菌株PH1进行双端测序，利用Unicycler（v. 0.5.0）进行de novo组装，CheckM（v. 1.1.2）评估基因组完整度与污染。分类学鉴定采用OrthoANI分析（OAT软件v. 0.93.1）及基于16S rRNA基因的邻接法（neighbour-joining, NJ）系统发育分析。基因组结构注释使用NCBI Prokaryotic Genome Annotation Pipeline（PGAP v. 6.10），次级代谢基因簇（biosynthetic gene cluster, BGC）分析使用antiSMASH（v. 8.0.1）网络服务器。功能注释通过BLASTp（v. 2.17.0）、基于KEGG数据库的BlastKOALA（v. 3.1）和Prokka（v. 1.14.5）完成。抗生素耐药基因鉴定使用ResFinder（v. 4.7.2），致病性评估使用PathogenFinder（v. 1.1）。
三、研究结果
1. 基因组组装指标与结构注释
菌株PH1测序产生2,531,752条双端原始读段，质控后保留2,393,938条读段用于组装，平均测序深度为238×。草图基因组包含87个contigs，N50为95,067 bp，基因组大小为2,872,165 bp，GC含量为66.63%。PGAP结构注释鉴定出编码序列（CDS）、RNA基因和假基因（表1），基因组完整度为99.70%，为该菌株后续功能分析提供了坚实基础。
2. 菌株分类学鉴定
OrthoANI分析显示菌株PH1与模式菌株Paracoccus bogoriensis BOG6（DSM 16578）的平均核苷酸一致性为96.33%，高于公认的物种划分阈值（95–96%）。基于16S rRNA基因的系统发育分析中，菌株PH1与P. bogoriensis BOG6以100%自展支持率聚类为独立种系分支，且与P. homiensis和P. fistulariae所在谱系（96%）及P. rhizosphaerae、P. beibuensis、P. hibisci和P. tibetensis亚分支（70%）存在明确区隔。该结果与ANI数据和前期基于16S rRNA基因的鉴定一致，确认菌株PH1的分类学地位为Paracoccus bogoriensis。
3. 四氢嘧啶生物合成基因簇与应激耐受功能
基因组中含四氢嘧啶（ectoine）生物合成基因簇，由ectABC基因组成，编码EctA（L-2,4-二氨基丁酸乙酰转移酶）、EctB（二氨基丁酸-2-酮戊二酸转氨酶）和EctC（L-四氢嘧啶合酶）。ect基因邻位发现天冬氨酸激酶编码基因ask和四氢嘧啶合成特异性转录调控因子编码基因ectR，构成ectABC-ask操纵子，并鉴定出四氢嘧啶分解代谢途径酶DoeABCD。序列比较显示菌株PH1的Ect相关酶及DoeABCD与P. bogoriensis BOG6相似性最高。基因组中还鉴定出热休克反应和热保护相关蛋白，包括DnaK和DnaJ分子伴侣、GroEL/GroES伴侣蛋白、小热休克蛋白Hsp20、Hsp33家族伴侣蛋白HslO及热诱导转录抑制因子HrcA。除Hsp20（相似度45.76–48.21%）在物种间差异较大外，其余蛋白在Paracoccus各物种中高度保守（80.77–96.22%）。这一结果提示菌株PH1的耐热性难以仅用经典热休克蛋白的序列分歧解释，可能涉及基因调控或表达水平的差异。
4. 类胡萝卜素生物合成基因簇及途径
类胡萝卜素BGC由核心基因crtE、crtB、crtI和crtY，以及叶黄素合成必需基因crtW、crtX和crtZ组成。基因组中还鉴定出idi基因，编码2型异戊烯基二磷酸δ-异构酶，对类胡萝卜素前体调控至关重要。BLASTp分析显示，菌株PH1的idi和crt基因编码蛋白与P. bogoriensis BOG6的相似度为96.28–99.40%。KEGG途径分析表明，甲基赤藓醇4-磷酸（methylerythritol 4-phosphate, MEP）途径是菌株PH1合成类胡萝卜素前体的主要代谢途径，该菌具备合成虾青素、叶黄素、角黄素、海胆酮（echinenone）和β-胡萝卜素等多种高价值类胡萝卜素的基因组能力，与HPLC检测到的多种类胡萝卜素产物相印证。
5. B族维生素、短链脂肪酸和吡咯喹啉醌生物合成潜力
Prokka注释与KEGG途径分析揭示菌株PH1基因组中存在与B族维生素（B1、B2、B6、B7、B9和B12）生物合成相关的推定基因。此外还注释到与短链脂肪酸（SCFAs）和辅酶吡咯喹啉醌（PQQ）合成相关的基因。SCFAs可提供能量、降低肠腔pH，并作为信号分子增强肠道屏障和调节宿主免疫应答。PQQ具有氧化还原循环能力、自由基清除活性及胺氧化相关活性，并具有改善记忆、保护心脏和缓解空肠黏膜炎症损伤等健康效应潜力。但研究人员明确指出，基因存在不能确证实际的生物合成或生物学相关水平的积累，仍需通过LC-MS/MS等实验验证。
6. 推定益生特性与初步基因组安全性
基因组中鉴定出与益生特性相关的推定基因：包括胆酸盐：钠共转运体（bile acid:sodium symporter, BASS）家族蛋白编码基因（可能通过钠转运输出胞质胆酸盐以保护细胞）、F0F1-ATP合酶复合物基因（atpA-atpH，参与酸性环境胞内pH稳态维持）、分子伴侣基因（groL、dnaJ、dnaK、clpB）、硫氧还蛋白系统基因（trxA、trxB）、膜完整性相关脂肪酸生物合成基因（fabA、fabB/fabF、fabD、fabG、fabI、fabZ）和细胞被膜生物合成基因（rfbB、rfbC、glmU、dapA）。鞭毛生物合成基因（flgA、flgC、flgG、flgH、flgK）与黏附功能相关，胞外多糖/荚膜生物合成及输出基因（galE、eps、cps、wza）可能参与生物被膜形成。安全性方面，ResFinder分析未检出氨基糖苷类、喹诺酮类、β-内酰胺类、四环素类、大环内酯类等多种抗生素类别的耐药基因。PathogenFinder分析显示病原性可能性为0.278%，属于低可能性。基因组中仅检出单一vapE基因，未发现其他毒力簇组分，通常不足以赋予致病性。
四、讨论与结论
讨论部分指出，菌株PH1与P. bogoriensis BOG6共享四氢嘧啶和类胡萝卜素生物合成等关键代谢特征。四氢嘧啶作为相容性溶质可保护细胞抵御渗透压、高温和紫外（UV）胁迫，与该菌株分离自约40 ℃温泉水并能于40 ℃生长的耐热特性相符。值得注意的是，BOG6以产虾青素著称，而PH1以产角黄素为主要类胡萝卜素（未优化条件下相对含量53.5%），这种差异可能归因于代谢流向的不同：菌株PH1的代谢流更倾向经CrtW介导的酮基化途径（β-胡萝卜素→海胆酮→角黄素），而非BOG6中CrtZ介导的羟基化途径（β-胡萝卜素→β-隐黄素→虾青素）。海胆酮是PH1中第二丰富的类胡萝卜素（相对含量4.9%），为首次在P. bogoriensis菌株中报道。但这一解释仍属假设，尚未获得CrtW与CrtZ相对活性的直接证据。B族维生素合成基因与MEP途径存在功能关联，如维生素B1（硫胺素焦磷酸, TPP）是MEP途径中1-脱氧-D-木酮糖5-磷酸合酶（DXS）的辅因子，维生素B2是黄素氧还蛋白/黄素氧还蛋白还原酶系统为MEP途径最终氧化还原步骤提供还原当量所必需。基因组安全评估虽显示PH1的病原性可能性低、耐药决定因子缺失，但受99.70%基因组完整度限制，未组装部分可能存在未检出的耐药或毒力决定因子，安全预测须经实验验证。
结论部分总结指出，OrthoANI值与16S rRNA系统发育分析支持菌株PH1为Paracoccus bogoriensis。基因组注释与特定BGC分析揭示了该菌在四氢嘧啶和类胡萝卜素生物合成方面的显著代谢能力，表明其具有应激耐受和多种类胡萝卜素生产相关的基因组潜力。基因组序列中还注释到B族维生素、短链脂肪酸（SCFAs）和吡咯喹啉醌（PQQ）生物合成功能及益生特性相关功能。未检出已知抗生素耐药决定因子为PH1提供了初步的基因组安全基线，表明其可传播耐药性的潜力较低。虽然这一基因组特征支持其在生物技术和环境应用中的潜在用途，但全面表型验证仍须在使用前开展。
