★ FEATURED ARTICLE 摘要蓝细菌Cyanobacteria是唯一演化出产氧光合作用oxygenic photosynthesis的原核生物曾从根本上改变地球大气与海洋化学组成。该类群呈现从单细胞球状、多重裂殖、简单丝状到具异形胞heterocyst和真分枝的显著形态多样性但其在单细胞unicellularity与多细胞multicellularity之间反复转换的基因组机制仍不明确。作者整合192个蓝细菌菌株与8个外类群共200个完整基因组的系统基因组学phylogenomics、贝叶斯分子钟Bayesian molecular clock与祖先蛋白编码基因含量重建对主要形态转换节点的基因获得与丢失进行量化。分析显示多细胞相关性状存在独立的获得与丢失功能汇聚functional convergence在信号转导、细胞分裂与转录调控通路中尤为突出。丝状谱系持续保留细胞协调与结构整合基因而逆转为单细胞常丢失运动性、调控复杂性与细胞骨架相关基因许多丢失基因可追溯到深层祖先创新说明复杂性降低往往逆转此前的基因组获得。上述结果整体支持蓝细菌形态演化具有模块化与可逆性并揭示功能汇聚如何重塑表型转换中的基因组架构。keyword蓝细菌、基因组演化、功能汇聚、多细胞性、单细胞性、系统基因组学、贝叶斯分子钟、祖先基因重建Yang L., Bianchini G., Paps J., Sanchez-Baracaldo P. (2026). The genomic origins of cyanobacterial morphological diversity. Proceedings of the National Academy of Sciences of the United States of America, 123(39), e2608186123. https://doi.org/10.1073/pnas.2608186123期刊Proceedings of the National Academy of Sciences of the United States of AmericaPNASIF9.5发表时间2026年9月21日在线收稿2026年3月14日接受2026年7月14日研究总结蓝细菌形态基因组起源研究背景与目标产氧光合作用与大氧化事件单细胞与丝状多细胞反复转换200个完整基因组系统框架主要结果祖先单细胞与早期丝状节点新生与丢失同源组不对称信号分裂转录功能汇聚多细胞基因深层起源与逆转丢失方法体系BUSCO单拷贝与OrthoFinderPhyloBayes分子钟与sMap祖态eggNOGCOG与Count扩张收缩应用意义模块化可逆形态模型早期地球碳氮氧反馈背景介绍传统蓝细菌分类按形态分为五类即单细胞型、多重裂殖型、简单丝状型、异形胞型与真分枝型但系统发育研究已证明该形态划分并非单系。已有分子生物学研究识别若干多细胞关键基因例如隔片连接蛋白基因sepJ与fraE、异形胞分化膜重组基因fraH、隔定位分裂蛋白基因cyDiv以及功能未完全明确的all2344但全基因组尺度上从单细胞到丝状、再逆转单细胞的基因创新与丢失路径仍不系统。作者以200个高质量完整基因组为基础联合单拷贝直系同源基因与核糖体RNA构建系统树用化石与地球化学标定进行深时定年并重建祖先形态与同源组含量以阐明产氧光合生物在地球早期如何通过基因获得、丢失与功能汇聚实现形态复杂性的动态演化。重要结果结果1 系统发育时序树与祖先形态重建显示丝状多细胞的多次独立起源及逆转作者以765个单拷贝BUSCO直系同源基因联合16S与23S rRNA构建最大似然树结合贝叶斯分子钟标定得到高支持度的蓝细菌时序系统框架。祖先形态二态重建显示全群最近共同祖先节点1为单细胞的概率达82.2%根龄约2.9 Ga95%最高后验密度2.637至3.055 Ga最早强支持的丝状节点A约2.3 Ga后验80.3%95%最高后验密度2.097至2.497 Ga与大氧化事件时间窗口重叠。微蓝细菌Micro-cyanobacteria在节点C约1.4 Ga发生一次逆转单细胞后验77.4%巨蓝细菌Macro-cyanobacteria在节点D、G、B分别发生独立逆转后验92.6%、99.9%、71.2%其中节点B之后又在节点F与节点H二次恢复丝状。假鱼腥藻目Pseudanabaenales节点E呈丝状强支持后验97.1%但定年不确定性大约0.8 Ga。上述结果说明丝状多细胞并非单次起源后线性保留而是多次获得、丢失与再获得。Figure1蓝细菌最大似然时序系统树与祖先形态重建。基于192个蓝细菌与8个外类群、765个单拷贝BUSCO同源基因及16S或23S rRNA在最佳分区模型下以千次超快自举构建时间轴由PhyloBayes联合rRNA与化石或地球化学标定估计。节点1为全群最近共同祖先单细胞后验82.2%根龄约2.9 Ga节点A为宏微蓝细菌共同丝状祖先后验80.3%约2.3 Ga节点E为假鱼腥藻目丝状祖先后验97.1%约0.8 Ga节点C为微蓝细菌单系逆转单细胞后验77.4%约1.4 Ga节点D、G、B为巨蓝细菌独立逆转单细胞后验分别为92.6%、99.9%、71.2%节点F与H为逆转后再次获得丝状。枝色区分Macro与Micro支系图右标注形态状态与主要目级类群。结果2 核心同源组功能组成与新生同源组演化热图揭示巨蓝细菌功能创新更强OrthoFinder将200个蛋白质组聚类为53,396个同源组homology groupHG其中7,968个可获得COGClusters of Orthologous Groups注释占核心可注释信息比例有限整体核心获得功能注释者约9%。在已注释部分代谢类基因约24%细胞过程与信号类约23%信息存储与处理类约15%。按从根节点到各现存谱系累计“新生同源组”群内出现、外类群缺失的热图显示巨蓝细菌各支系的功能HG计数普遍高于微蓝细菌该差异与巨蓝细菌更大的细胞直径、异形胞分化及底栖或水华生态策略相一致说明形态复杂性提升伴随更广的特有功能基因创新。Figure2蓝细菌核心同源组COG功能构成与新生同源组计数热图。左图按COG大类统计可注释核心同源组比例信息存储与处理含翻译、转录、复制修复等细胞过程与信号含细胞分裂、信号转导、细胞壁膜、运动、细胞骨架等代谢含能量、碳水、氨基酸、核苷酸、辅酶、脂质、无机离子与次级代谢右图热图按系统树各分支展示自根至现存物种的新生同源组逐功能类别累计数量行按COG类别、列按代表性蓝细菌 lineage 排列Macro-cyanobacteria整体计数高于Micro-cyanobacteria反映不同形态与细胞大小类群的功能创新强度。结果3 节点水平新生与丢失同源组动态及两次初级丝状转换的相反基因组轨迹沿系统树逐节点比较新生与丢失同源组深层祖先节点呈净基因获得浅层形态转换节点多呈净丢失。首次初级单细胞向丝状转换节点A获得1,792个新生同源组伴随复制修复、膜生物合成与信号转导类扩张第二次独立丝状转换假鱼腥藻目节点E则以大量丢失为主涉及2,823个同源组丧失。该对比表明相同丝状表型可通过“基因组加法”或“基因组减法”分别实现前者依赖新调控与结构基因创新后者可能通过去除自主代谢与冗余调控、强化细胞间依赖来获得协同形态。Figure3系统发育各节点新生与丢失同源组数量。柱状或折线按节点展示新生同源组数在蓝细菌群内出现且外类群缺失与丢失同源组数在群内缺失且于姊妹类或外类群存在深层节点如A1至A3与A以净获得为主浅层形态转换节点以净丢失为主。节点A初次丝状转换获得1,792个新生同源组节点E假鱼腥藻目丝状转换伴随2,823个同源组丢失图中同时标注Macro与Micro主要转换支系的净变化方向。结果4 形态转换COG功能类别汇聚以信号转导、细胞分裂与细胞壁丢失为主针对丝状转换节点A、E、F与单细胞逆转节点B、C、D、G作者将同源组按新生、丢失、扩张与收缩注释到COG类别。丝状转换中新生HG无一致共享类别扩张HG仅在信号转导T、翻译后修饰或伴侣O、氨基酸代谢E有限汇聚收缩无一致模式但丢失HG在防御V、伴侣或蛋白周转O、碳水G、核苷酸F等类别出现汇聚性减少。单细胞逆转中新生HG汇聚仅见于节点B与G之间丢失HG在所有四个节点共同影响细胞壁或膜或包膜生物合成M、信号转导T、复制重组修复L扩张与收缩汇聚较弱。结果说明形态去复杂化比形态复杂化更具可重复的功能丢失模式。Figure4形态转换节点的COG功能汇聚分析。A为向丝状形态转换汇总节点A、E、F的各类同源组演化模式新生HG未显示共享COG类别扩张HG在T信号转导、O翻译后修饰或伴侣、E氨基酸转运代谢有限汇聚收缩HG无一致模式丢失HG在V防御、O伴侣或蛋白周转、G碳水化合物转运代谢、F核苷酸转运代谢呈汇聚减少。B为向单细胞逆转汇总节点B、C、D、G新生HG仅B与G间有限汇聚丢失HG跨全部节点共同涉及M细胞壁或膜或包膜、T信号转导、L复制或重组或修复扩张与收缩HG汇聚较弱整体反映运动、调控与维持功能被系统精简。结果5 多细胞性相关同源组深层起源及其在单细胞逆转中的优先丢失作者以Boden等定义的多细胞基因集做BLASTp回溯定位其所属同源组在系统树中的起源与丢失。fraE与sepJ已存在于全蓝细菌最近共同祖先提示细胞间连接元件早于丝状形态all2344、cyDiv、sepI在节点A3假鱼腥藻目与宏微共同祖先之前出现可能支持不完全分离导致的物理关联节点A进一步获得fraC、fraD、hetR、hetZ、bglK、patU3构成隔片连接、发育调控、细胞间通信与异形胞分化核心。比较转换节点丢失HG与深层新生HG来源发现单细胞与部分丝状逆转更优先丢失蓝细菌特异性创新同源组而非前蓝细菌深层保守同源组其中导向聚球藻类的节点C丢失所考察全部多细胞基因拟色球藻目等逆转亦系统缺失fraC与fraD。该结果说明复杂性降低常按“逆向撤销”祖先多细胞创新的方式进行。Figure5多细胞性相关同源组的系统起源及其在形态逆转中的丢失来源。上图按系统树展示fraE、sepJ于全群最近共同祖先出现all2344、cyDiv、sepI于节点A3出现fraC、fraD、hetR、hetZ、bglK、patU3于节点A集中获得下图比较各形态转换节点丢失同源组与深层祖先新生同源组的重叠显示逆转节点优先丢失源自蓝细菌根至节点B区间的特异性创新同源组而非更早期外类群保守同源组。插入表列明各多细胞基因对应HG、所属COG功能组及在单细胞谱系中的保留或缺失状态。方法学参考200个完整基因组从NCBI下载并BUSCO评估完整性CD-HIT构建非冗余蛋白与OrthoFinder聚类同源组765单拷贝BUSCO加16S23S用MAFFT比对trimAl修剪IQ-TREE最佳分区最大似然建树与千次自举PhyloBayes固定拓扑贝叶斯分子钟并加入化石或地球化学标定sMap步进石模型重建二态祖先形态与五态比较系统发育解析脚本提取新生与丢失同源组Count基于Wagner简约重建扩张与收缩基因家族eggNOG-mapper注释COG并统计功能类别BLASTp追踪多细胞基因所属同源组起源与丢失上述流程中可复用要点包括以BUSCO cyanobacteria_odb10评估基因组完整性对高缺基因但为共生或极端简化的菌株单独标注用CD-HIT以完全一致阈值生成每基因组非冗余代表序列再以OrthoFinder进行全群同源组聚类系统树使用MAFFT自动比对、trimAl自动化修剪IQ-TREE ModelFinder选最佳分区千次超快自举评估支持度分子钟采用PhyloBayes的GTR加类别模型以rRNA与若干化石或地质标定点约束深时节点祖先形态用sMap比较非对称速率模型与均等速率模型、以边际似然混合结果基因家族动态用Count做Wagner简约祖先拷贝数推演功能注释用eggNOG-mapper v2获得COG类别多细胞基因追踪以BLASTp取每查询最优命中并映射回OrthoFinder同源组从而避免仅看单基因而忽略旁系同源扩张。总结作者通过对200个蓝细菌及相关外类群基因组的系统发育、分子钟与同源组动态分析证实产氧光合蓝细菌的形态演化不是从单细胞到多细胞的单一线性过程而是以模块化、可逆转的方式反复进行。全群最近共同祖先更接近单细胞约2.9 Ga丝状多细胞至少在宏微共同祖先节点A约2.3 Ga、假鱼腥藻目节点E更晚时期独立起源两次前者以大量新基因获得为特征后者以大量基因丢失为特征。微蓝细菌与巨蓝细菌均发生多次逆转单细胞部分逆转支系又可再次获得丝状形态。功能层面丝状创新依赖信号转导、细胞分裂、转录调控、细胞壁与细胞骨架整合单细胞逆转则趋同丢失运动、调控、细胞壁生物合成与复制修复类基因且优先去除蓝细菌深层特有创新而非更古老保守基因。该框架为早期地球大氧化、新元古代氧化及氮碳循环模型中纳入“形态可塑性—基因组流”耦合提供了机制化依据。受深时定年不确定性与五态形态模型可识别性限制后续应结合更多化石标定、基因水平转移筛查与实验验证进一步细化多细胞工具箱的因果功能。参考文献Yang L., Bianchini G., Paps J., Sanchez-Baracaldo P. (2026). The genomic origins of cyanobacterial morphological diversity. Proceedings of the National Academy of Sciences of the United States of America, 123(39), e2608186123. https://doi.org/10.1073/pnas.2608186123分析代码与数据集GitHubjiushixiaoyang或Cyanobacteria_morphology_diversityhttps://github.com/jiushixiaoyang/Cyanobacteria_morphology_diversity基因组数据源NCBI GenBank与RefSeqhttps://www.ncbi.nlm.nih.govBUSCO基因组完整性评估https://busco.ezlab.orgOrthoFinder同源组聚类https://github.com/davidemms/OrthoFinderIQ-TREE系统发育建模https://iqtree.github.ioPhyloBayes贝叶斯分子钟http://www.atgc-montpellier.fr/phylobayessMap祖先状态重建https://github.com/gbrianchini/sMapeggNOG-mapper功能注释与COGhttps://eggnog-mapper.embl.deCOG数据库https://www.ncbi.nlm.nih.gov/research/cogRNAmmer rRNA预测MAFFT比对与trimAl修剪相关文档见各软件官方站点。
阅读完成 · 觉得有帮助？
