《Evolution & Development》：Modeling the Histogenetic Switch Between Enameloid and Enamel
编辑推荐：
釉质（enamel）和釉质样组织（enameloid）是高度矿化的组织，存在于脊椎动物牙齿结构的表面。共同的发育特征表明，微小的变化可能驱动了牙本质-釉质（dentine–enamel）和牙本质-釉质样（dentine–enameloid）结构之间的多次转换。
釉质（enamel）和釉质样组织（enameloid）是高度矿化的组织，存在于脊椎动物牙齿结构的表面。共同的发育特征表明，微小的变化可能驱动了牙本质-釉质（dentine–enamel）和牙本质-釉质样（dentine–enameloid）结构之间的多次转换。已有假说认为，上皮细胞活性的延迟（H1）或延长（H2），以及上皮细胞产生速率的改变（H3），可以解释这些转换。为了检验这些演化发育（evo-devo）假说，研究人员利用现生脊椎动物共有的简化特性，构建了一个基于细胞的组织发生模型（histogenetic model）。通过改变分泌细胞参数，模拟能够再现牙本质、釉质样组织和釉质的多种比例，从而探索它们在釉质样-釉质转换中的作用。研究人员对12个细胞参数的探索表明，H1和H3在一定程度上可以解释这种转换，而H2可能需要与其他修饰结合。结果还表明，矿化前沿形成的位置和时间的变化是限制釉质-釉质样-牙本质发育的关键。因此，除了基因功能的获得或丧失之外，微小的发育修饰也可能在牙齿组织的出现和演化中发挥了关键作用。
论文解读：模拟釉质样组织与釉质之间的组织发生转换
研究背景：脊椎动物的牙齿结构由牙本质（dentine）、釉质样组织（enameloid）和釉质（enamel）等矿化组织构成，这些组织在矿化程度、基质组成和发育机制上存在差异。釉质和釉质样组织均为高度矿化的覆盖组织，但釉质由上皮细胞来源的釉质基质蛋白（enamel matrix proteins, EMPs）构成，而釉质样组织则具有间充质和上皮的双重来源。化石和现生脊椎动物的分布模式提示，在演化过程中可能发生了多次牙本质-釉质与牙本质-釉质样结构之间的转换，尤其是早期硬骨鱼类中釉质样向釉质的转变，这标志着釉质的演化起源。现有假说包括：延迟上皮细胞活动（H1）、延长上皮细胞活动（H2）以及改变上皮细胞产物比例（H3，即从胶原为主的分泌转向EMPs为主）。然而，这些假说缺乏定量验证，且可能低估了发育修饰的多样性。因此，本研究旨在通过计算模型系统性地探索细胞参数变化对组织形成的影响，以检验这些演化发育（evo-devo）假说。
研究内容：研究人员开发了一个基于细胞的组织发生模型，模拟牙本质、釉质样组织和釉质的形成过程。模型采用简化的现生脊椎动物特征，并利用EmbryoMaker软件进行三维细胞尺度模拟。模型包含三个发育阶段：基质形成、矿化和成熟。研究人员识别了12个关键参数，包括胶原、囊泡和EMPs的分泌起始时间、持续时间和速率，并通过网格搜索方式进行了551,124次模拟。模拟结果根据七项生物学标准进行筛选，最终保留42,571个符合条件的案例，其中33,435个为釉质样案例（釉质比例<5%），1,912个为釉质案例（釉质比例>95%）。
技术方法：主要技术方法包括：(1) 使用EmbryoMaker软件构建细胞尺度模型，模拟细胞分泌和扩散过程；(2) 进行主成分分析（PCA）和线性判别分析（LDA）以识别影响组织类型的关键参数；(3) 通过成对比较分析，枚举所有可能的参数变化组合，识别能够独立驱动釉质样-釉质转换的“主要转换路径”（primary shifts）；(4) 计算每种路径的“参与率”（implication rate）以评估其重要性。计算在法国国家自然历史博物馆的PCIA平台上完成。
研究结果：
1. 主成分分析（PCA）：前两个主成分轴仅解释了25.2%的方差，无法清晰区分釉质和釉质样案例，表明多种参数配置可能产生相同结果。主要贡献参数包括DEE（EMPs分泌持续时间）、DCE（胶原分泌持续时间）和SVO（囊泡分泌起始时间）。
2. 线性判别分析（LDA）：交叉验证成功率为99.5%，SVO（囊泡分泌起始时间）和SEE（EMPs分泌起始时间）是区分两种组织的关键参数。SVO在釉质样案例中较高，而SEE在釉质案例中较高，支持H1假说。
3. 比较分析：通过63,927,720对成对比较，识别出186,355种独特的参数变化组合，其中38种组合构成独立的主要转换路径，每种路径包含1至6个参数变化。最高参与率的路径是SVO-（囊泡分泌提前，58.6%），其次是H1（SEE+，53.0%）。H2（DEE+）必须结合其他变化（如DVO+）才能形成有效路径，最高参与率为14.6%；H3（RCE-和REE+）也需结合其他变化，单独不能驱动转换。补充分析表明，当RCE、SCE和DCE同时降至零时，H31/2（RCE-）可以独立触发转换。此外，研究还发现矿化前沿的形成时间和位置是决定组织类型的关键因素，SVO-通过影响矿化前沿的时空位置来促进釉质形成。
讨论：研究结果支持H1作为独立驱动转换的假说，而H2和H3作为单一因素无效，需结合其他修饰。研究还提出了5种新的独立转换路径，其中SVO-的参与率最高。矿化前沿的形成时机和位置被强调为核心机制：提前形成靠近上皮细胞的矿化前沿可隔离胶原和EMPs，从而促进釉质生成。在体内研究中，这些发现提示部分发育差异可能只是伴随性变化，而非关键驱动因素。例如，H31/2和H32/2在体内常同时出现，但模拟显示它们不共存于同一主要路径中，表明其中一项变化可能发生在转换之后。
结论：研究成功开发了细胞尺度组织发生模型，模拟了釉质、釉质样和牙本质的形成过程。通过系统探索12个参数，识别出38种独立的主要转换路径，其中SVO-和H1（SEE+）是最重要的单一驱动因素。研究强调，微小的发育修饰，如分泌时间或矿化前沿位置的变化，可能在牙齿组织的演化中发挥关键作用，远非简单的基因功能得失所能解释。
