《Materials Today Advances》：Statistical analyses of feature orientation in 3D reticulate biological assemblies
编辑推荐：
生物结构与功能密切相关，功能线索常隐匿于结构复杂的三维细节之中。分支状、纤维状和多孔状生物 assembly 的这些细节可通过体容积成像（volumetric imaging）在多重尺度上进行研究。然而，三维可视化本身无法对样本结构进行定量评估。结构属性的方向
生物结构与功能密切相关，功能线索常隐匿于结构复杂的三维细节之中。分支状、纤维状和多孔状生物 assembly 的这些细节可通过体容积成像（volumetric imaging）在多重尺度上进行研究。然而，三维可视化本身无法对样本结构进行定量评估。结构属性的方向性量化——如反映结构元素取向与共 alignment（各向异性，anisotropy）的纹理（texture）——对于解读结构-功能关系、预测力学行为以及优化设计具有重要意义。研究人员开发了 VectoRose——一种用于球形直方图可视化与三维向量数据定量分析的新型开源软件包，并利用方向统计学方法，对三组生物 assembly（分支状、纤维状和多孔状）进行了案例分析。VectoRose 产生的输出结果为全局与局部结构模式中的结构-功能关系提供了定量认识。对于分支结构，通过比较月桂叶与驯鹿地衣的运输网络， researchers 量化了平面型与体积填充型树状模式之间的差异。对鸟类和蜥蜴蛋壳膜中三维纤维 assembly 的统计分析可能反映了蛋孵化策略的进化分歧。最后，对绵羊股骨多孔小梁骨模式的研究为探讨骨生物力学适应性的结构-功能关系提供了一种定量方法。借助适用于同源或异源结构的稳健三维各向异性统计描述，研究人员预见定量纹理分析将为材料研究其他领域的结构-功能研究赋予可重复性、细微差别与精妙之处。
**研究背景与问题**
生物实体的结构组织反映其用途与性能，这一"形式追随功能"（form follows function）范式已存在数百年。然而，结构可能承担多种功能，且某些结构特征在进化中长期保留，即使其原有功能已不再相关，这使得建立结构与功能之间的联系变得非易事。生物结构从分子尺度到生物体尺度呈层级排列，包括细胞状、镶嵌状、纤维状、网状或分支状等常见模式。当代三维成像方法可从密度、孔隙率、局部有序与无序、形态计量属性、纹理或网络连通性等方面提供定量认识，但若结构-功能联系难以捉摸，不同结构之间的比较便可能被轻易否定。
结构各向异性（structural anisotropy）——即主导全局与局部取向及局部邻域内组分共 alignment 的程度——是解答此类问题的关键。小梁骨（trabecular bone）长期以来一直是结构各向异性分析的热门对象，在力学加载背景下，结构各向异性直观上是功能历史的替代指标。然而，经典各向异性分析方法存在局限：感兴趣区域必须足够大、几何形状不能干扰计算（如应为球形），且常仅报告各向异性程度而忽略取向信息。Reznikov 等人近期提出的方法虽能计算贯穿整个三维图像的密集向量场以保留解剖学背景并解析局部各向异性变异，但这些密集三维场无法以有意义的数值方式直接报告，亟需统计汇总工具。
**研究开展与结论**
为此，研究人员开发了 VectoRose 软件包，用于三维各向异性向量场的可视化与定量分析。该软件通过生成三维球形直方图（spherical histograms）并计算方向统计量，实现对所有取向的汇总，且不受样本空间位置限制。研究聚焦于 Woodcock 形状参数（γ ?，γ>1 为 cluster 分布，γ<1 为 girdle 分布）和强度参数（ζ?，ζ?0 表示紧凑分布，ζ→0 表示弥散分布），对分支结构、纤维结构和多孔结构三类案例进行了分析，论文发表于《Materials Today Advances》。
该研究的重要意义在于：VectoRose 将经典各向异性方法的易报告性与基于 mapping 的方法所保留的丰富空间背景相结合，为三维结构-功能分析开辟了新途径，适用于生命科学和材料科学等领域。
**主要技术方法**
研究人员采用的核心技术流程包括：体容积图像获取（X 射线微计算机断层扫描（μCT）或聚焦离子束-扫描电子显微镜（FIB-SEM）断层成像）、图像分割（强度阈值分割或深度学习辅助分割）、基于表面网格的三维各向异性 mapping（采用 Reznikov 等人提出的方法，在 Dragonfly 3D World 中实现，通过用户指定影响半径的球形采样体积计算局部主导取向和各向异性程度），以及 VectoRose 分析（导入 NumPy 数组格式向量场，排除各向异性程度为零的向量，将负 z 分量向量反演以确保取向限于上半球，采用 Tregenza 球面离散化（5806 个近似等面积面元）构建球形直方图，计算 Woodcock 参数、平均合成长度、球面中位方向等方向统计量，并生成边缘直方图、联合直方图和条件直方图）。样本来源包括：教学收集的月桂叶、驯鹿地衣，先前研究项目的家鸡和普通豹纹守宫蛋壳膜，以及先前研究项目的绵羊股骨。
**研究结果**
**案例研究 1：分支结构——运输网络的维度差异**
研究人员比较了月桂叶（Laurus sp.）的叶脉网络与驯鹿地衣（Cladonia sp.）的菌丝簇三维 arrangement。
月桂叶呈现网状叶脉，整体取向局限于叶平面，低阶叶脉在局部呈三维形态并产生面外向量分量。各向异性程度分布普遍较低（模态值 0.1–0.2），球形直方图上形成弱 girdle 分布（γ=0.396，ζ=1.499），极直方图进一步证实其近似平面性质。
驯鹿地衣则为体积填充的分支菌丝网络，菌丝向各方向分支，共 alignment 可忽略。整体各向异性程度直方图显示低值（模态值 0.2–0.3），球形直方图弥散且斑驳，强度参数极低（ζ=0.267）确认了地衣中缺乏优选取向，极直方图近似理想均匀分布。
该案例量化了直观差异：叶运输网络大致跨越平面，而地衣体跨越体积。
**案例研究 2：纤维结构——蛋壳膜纤维 assembly 的功能分歧**
研究人员比较了家鸡（Gallus gallus domesticus）刚性蛋壳与豹纹守宫（Eublepharis macularius）柔性蛋壳的 eggshell membrane 纤维 arrangement。
家鸡蛋壳膜为平行于蛋表面的直纤维层，纤维在平面内呈多样取向。各向异性程度集中于 0.4–0.6，球形取向直方图显示强 girdle（γ=0.042，ζ=3.017），对应于蛋壳表面切平面。
豹纹守宫蛋壳膜纤维 arrangement 约束较少，细纤维呈波浪形或卷曲状，可见旋涡（swirl），集体取向渐变，部分纤维显示面外取向。各向异性程度分布窄，中心约 0.4–0.5。球形取向直方图显示大斑点，看似由两个弥散 cluster 叠加于极 faint girdle 之上，呈现 girdle 与 cluster 之间的过渡分布（γ=0.820，ζ=1.833）。极直方图揭示了特定的主导方位角。约 20°的取向 cluster 偏移表明纤维以 ±10° 相对于蛋表面波动，可能使膜能像针织纺织品般拉伸收缩，±10° 波动可产生约 1.5% 的额外长度（1/cos(10°)≈1.015），对应体积增加约 5%，与半柔性龟蛋在高湿与低湿条件下孵化结束时的质量差异观察相符。而鸟类蛋壳膜因被外层方解石层锚定而无法收缩或扩张。
**案例研究 3：多孔结构——小梁骨中的应力轨迹**
研究人员分析了绵羊（Ovis sp.）股骨两个部位：股骨头（与髋关节的球窝关节相关）和外侧髁后部（与膝关节的铰链样关节相关）。
股骨头部观察到大量垂直于软骨下骨表面的细薄小梁，各向异性程度低，向股骨颈逐渐增高。主导取向向量图显示大多数局部取向在矢状 yz 平面内呈放射状。球形直方图汇总为弥散的部分 girdle 模式（γ=0.794，ζ=1.609），近似与矢状平面对齐，而非占据上半球的弥散 cluster。极直方图进一步表明这是主要在球面一侧可见的部分 girdle。
外侧髁中，大量细薄小梁亦垂直于软骨下骨表面取向。各向异性程度直方图显示中高度值，模态值较高（0.4–0.5），但与股骨头相比左尾更长。取向图揭示FC stark 差异：大多数向量呈现单向（z 对齐）取向。球形直方图揭示 cluster 模式（γ=2.104，ζ=1.847），垂直取向，而非矢状 girdle。该 cluster 对应于羊自然站立位时股骨远侧端与胫骨近侧端的大致接触点。
这些意外结果——股骨头为 girdle 而非半球形 cluster，外侧髁为集中 cluster 而非 girdle——揭示了骨的结构-功能关系。研究人员初始假设基于解剖学上可能的关节运动范围，但发现强调了更窄的生理性习惯运动范围：绵羊髋关节最大运动范围为屈曲-伸展，内收-外展最小；且家养绵羊多数时间站立而非移动，解释了解剖范围向习惯生理性运动范围的限制。小梁各向异性反映的是重复发生的运动，而非所有可能的运动。
**讨论总结**
VectoRose 将先前不相容的两个维度——易报告性与丰富空间背景的 inclusion——加以调和。与经典方法类似，其输出为紧凑的清晰图件和数值结果，便于报告和比较；又通过 leverage 基于 mapping 的各向异性计算方法，捕捉整个体积内三维形态的丰富性和复杂性。三个案例分别被选为直观差异（2D 与 3D）、感知相似性（蛋壳膜纤维）和预期符合关节生物力学原理（骨）的代表，但有趣的是，三维各向异性量化在案例 2 和 3 中反驳了初始假设。
VectoRose 的进一步功能包括：从取向直方图表面选取特定取向以提取含特定取向结构的空间位置（类似于暗场选区电子衍射，DF-SAED）；绘制条件直方图以可视化特定取向制度内的幅度分布或特定幅度水平的不同取向分布；构建联合直方图显示所有幅度与取向组合的向量分布。该软件的通用性通过拓扑 distinct 的架构（分支状、纤维状、多孔状）得以说明。
**研究结论**
研究人员新开发的分析方法为密集各向异性场提供了汇总信息，从各向异性程度、主导结构取向以及这两个变量之间的相互作用中衍生出定量认识。这种方法为生命科学和材料科学等领域的结构-功能研究开辟了丰富机遇。
