甲藻染色体真的是胆甾相液晶吗？

胆甾相模型对Hi-C的预测

胆甾相液晶（CLC）模型提出于五十多年前，它把甲藻染色体描述为一叠扁平的盘层。每个盘层内部，染色质纤维平行于一个指向矢排列；相邻盘层之间，这个指向矢按固定角度旋转，由此定义胆甾相螺距。盘层按与DNA一级序列相同的顺序依次排列，染色体外环从盘层向外伸出，连接同一个盘层或相邻盘层。由于这些环被认为比凝聚的核心更易接近，它们也被推测是转录活跃基因所在的位置。

论文的第一步，是把这一几何描述转化为可检验的预测。作者生成了10^6个模拟CLC结构，数量大致相当于一次典型bulk Hi-C实验所测序的细胞数，并随机加入长度服从指数分布的染色体外环。随后他们计算这群结构会产生怎样的Hi-C接触图。结果是一种网格状图案：高接触概率的条纹彼此平行，并垂直于中心对角线；条纹的数量和间距由盘层的数量和直径决定。改变环所占比例、螺距以及盘层直径的渐变，这一网格依然存在，说明它是稳健的特征，而不是参数调出来的假象。

模拟网格与真实Hi-C不符

来自Fugacium kawagutii和Symbiodinium microadriaticum的实验Hi-C图在性质上明显不同。它们的主要特征是拓扑关联结构域（TAD）：这些一级序列区块彼此接触的频率，远高于它们与染色体其余部分的接触频率。CLC模拟没有产生这样的区块。定量的TAD指标也反映了这一点——对CLC而言，绝缘分数和方向性指数以固定周期平滑振荡，而在真实染色体中则不规则、突变式地变化。

这种不匹配也体现在定量层面。CLC模型给出的接触图对实验数据的预测能力有限，Pearson相关系数在F. kawagutii中为r = 0.62，在S. microadriaticum中为r = 0.59。这是论文的第一条证据，它比较的是模型生成的图与实验测得的图，而不是对分子结构的直接观察。

两种重建方法都没有得到盘层

为了弄清真实Hi-C数据究竟支持什么样的构象，作者没有只用一种建模方法，而是用了两种，以明确降低算法偏差。用CSynth做共识建模时，Hi-C接触被当作单体之间的弹簧力，输出单一的最优拟合构象；该方法以每个单体5 kb的分辨率，应用于全部F. kawagutii和S. microadriaticum染色体支架。用Integrative Genome Modeling（IGM）做群体建模时，则通过最大似然估计生成一群在统计上与bulk Hi-C数据一致的构象；它为每个物种最大的染色体支架生成了100个构象，分辨率同样是5 kb。

两种方法都没有恢复出胆甾相盘层。CSynth构象紧凑、呈椭球状；IGM构象则可能拉长或呈球形，包含多个紧凑结构域和较长的染色体外环。两种方法之间也存在有启发性的分歧：CSynth的椭球形状与透射电镜中看到的形状相似，而IGM有时会产生哑铃形和不规则形状，这些在显微学文献中较少见。但IGM构象对实验Hi-C图和接触概率曲线的复现要好得多，F. kawagutii的相关系数为r = 0.91，S. microadriaticum为r = 0.84，其模拟图沿对角线呈现方块，而不是条纹。值得注意的是，CSynth可以被调成重现类似CLC的构象，而IGM不能。

第三个物种支持该模式但有保留

Breviolum minutum提供了部分重复验证。其可用的Hi-C数据估计只覆盖基因组的39%，而F. kawagutii为70%，S. microadriaticum为72%，因此它的支架可能并不对应完整染色体；作者据此只重建了最大的10个Hi-C支架，而不是整条染色体。共识建模给出的是带大量染色体外环的球形构象（每个支架5.3个环），与另外两个物种中环很少的椭球不同；群体建模给出类似的球形，每个构象3.7个环。作者把这些差异归因于支架不完整，而不是B. minutum真实的结构差异。即便如此，IGM群体仍以r = 0.93复现了观测到的Hi-C图，两者都显示类似TAD的方块，而不是类似CLC的条纹。

TAD可以嫁接到胆甾相堆叠上

论文并没有声称CLC几何在原则上与TAD不相容。作者按观测到的dinoTAD尺寸把胆甾相盘层分组，并在组与组之间引入位点特异性边界，只在组内打乱盘层顺序、组间不打乱，结果模拟接触图沿对角线出现了高接触概率的方块——这正是TAD的特征。此时绝缘分数和方向性指数也变得不规则，与真实数据一致。

这是一个相容性结果，而不是对模型的挽救。作者要说明的是：位点特异性边界加上构象异质性（哺乳动物细胞中环挤压的组成要素），才是产生类似TAD图案的原因；单靠胆甾相堆叠做不到。该实验指出了哪个额外要素是必需的，也说明TAD的存在本身不能作为胆甾相盘层的证据。

有序性是局部的而非胆甾相的

作者用取向序参量（OOP）量化取向有序性，该参量由不同一级序列间隔处单体之间的切向量计算，并对其做傅里叶变换。简单螺旋和连续的CLC盘层会给出振荡特征；平衡球状聚合物是没有取向有序的受限随机游走，其OOP单调衰减到零，傅里叶变换没有峰。F. kawagutii、S. microadriaticum和B. minutum的共识构象两者都不符合。它们在约10^5 bp处出现负峰，表明该尺度上存在部分反平行取向，而这一尺度对应基因阵列的平均长度。在dinoTAD尺度（约10^6 bp）上未检测到取向有序，可能是因为这些结构域在结构上高度异质。

第二种度量是局部向列序参量，在75个结构或构象的小体素和大体素内分别计算，用来捕捉依赖三维空间邻近而非序列间隔的有序性。两项度量共同支持论文的概括：甲藻染色体具有中等程度的取向有序和向列有序，高于随机预期，但低于CLC。作者还指出，OOP在基因阵列长度处的峰提示基因阵列可能以反平行方式折叠，而发散基因阵列中同向基因的转录可能产生这种排列——这是一个作为可能性提出的解释，而不是已被证明的机制。

转录并未集中在染色体边缘

CLC模型除了结构预测，还有功能预测：如果染色体表面的染色体外环是易接近的位点，那么转录活跃基因就应在那里富集。作者把RNA-seq读段映射到模拟构象上来检验这一点。每个物种最大的75个染色体支架被等比缩放到相同大小，长轴水平放置；读段在柱坐标系的体素中计数，沿方位角积分，再除以柱坐标雅可比因子，以校正每个体素中DNA量的差异。

染色体表面的转录水平与内部相当，并没有升高。作者据此认为，甲藻染色体虽然高度凝聚，但仍可被RNA聚合酶接近。这一结论针对的是模拟构象中转录的空间分布，它直接否定了CLC模型中染色体外环部分所依据的可及性理由。

重建方法无法确定的部分

这些结论建立在由Hi-C数据推断出的聚合物构象之上，作者也明确说明了这种推断的局限。不同计算方法可能对同一输入给出不同输出，这正是使用两个软件包的原因；两者在染色体形状上的分歧被如实报告，而不是被掩盖。B. minutum的分析受基因组覆盖不完整的限制。此外，围绕结构的长期争论很大程度上由二维透射电镜横截面塑造，而纤维取向从未在三维中被成像；对条带图案的解释一直依赖把三维CLC模型投影到二维平面。作者指出，产生此类条带的染色体在整体制片后呈螺旋状，而近期一项研究在横切面中发现纤维并不平行。

论文自己提出的后续优先方向正来自这一空白：在甲藻染色体内部进行三维成像和染色质纤维追踪，并对新生RNA和转录进行三维成像。论文还提出，应把FISH、DamID和SPRITE等正交数据类型整合进群体建模——IGM已经能做到这一点，但这些数据在甲藻中尚不可得。作者进一步指出，他们的构象与其他已提出的模型也不一致——包括相互缠绕的双螺旋、扭曲的环面、多线带、层级螺旋，以及带大量环的螺旋核心——因为Hi-C图缺少环状染色体特有的X形图案，沿对角线没有可变的接触概率，也缺少多染色体缠绕所应有的显著跨染色体接触。

这篇论文支持的判断

最站得住脚的读法是一个分开的裁决，而不是一概否定。整体性主张——甲藻染色体是胆甾相盘层的堆叠，且转录局限于外周环——在四项独立比较中都不成立：CLC模拟图缺少TAD、接触概率标度与CLC不符、构象比CLC更松散且有序性更低、转录未在表面富集。较弱的结论则保留下来：这些染色体是凝聚的聚合物，在基因阵列尺度上有局部有序，但没有全局的胆甾相结构，而这种凝聚本身并不阻断转录。

对于评估甲藻染色体结构主张的读者，实际后果是：TAD的存在不能作为胆甾相组织的诊断依据，因为作者只需在胆甾相堆叠中加入位点特异性边界，就能生成类似TAD的图。真正有区分力的证据是网格状条纹图案和接触概率标度，而真实数据在这两项上都不符合。作者把结论限定在所研究的三个物种内；B. minutum支架不完整，是把该物种的形状结果视为暂定、而非物种差异的最明确理由。

本文引用的文献