指纹引导的光声/超声显微成像实现整合汗腺导管-微血管单元的三维定量 mapping

《Photoacoustics》：Quantitative 3D mapping of the integrated sweat duct-microvasculature unit via fingerprint-guided photoacoustic/ultrasound microscopy

　　

皮肤体温调节依赖于汗腺导管与真皮微血管之间的协同相互作用。然而，单一模态成像无法同时可视化并表征这一功能耦合。本研究提出了一种指纹引导的双模态光声/超声（PA/US）显微成像框架，用于整合汗腺导管-微血管单元（SDMU）的三维 mapping。该策略将高频超声与深度学习（U-Net）相结合以提取汗腺导管和指纹特征；同时，通过深度依赖性黑色素的结构相似性指数（SSIM）分析，在PA成像中精确定位微血管分布。关键在于，指纹引导的非刚性配准消除了运动伪影，使得在生理过程中能够动态追踪同一SDMU。定量分析表明，微血管的功能耦合高度集中在汗腺导管周围75微米半径范围内，且超声可见汗腺导管结构的扩张与微血管招募呈强正相关。此外，一种生理学信息驱动的操作性双响应标准进一步揭示了功能激活的显著空间异质性，从而有助于稳健区分激活与未激活的SDMU。本研究不仅将SDMU确立为研究皮肤体温调节的可量化功能单元，还强调了其作为评估糖尿病相关神经血管功能障碍候选生物标志物的潜力。

研究背景：皮肤体温调节关键依赖于外泌汗腺导管与邻近皮肤微血管的紧密协同活动。汗腺导管通过出汗调节体温，这一过程需要充足的血液供应，因此体温升高常伴随血管舒张；同时汗腺导管形态也随温度变化而发生微妙改变。汗腺导管与微血管均为受交感神经系统调控的功能性微器官，在热应激下交感神经驱动血管舒张与汗腺导管激活，二者并非独立运作，而是形成由交感神经系统协调的功能耦合单元，这一"汗腺导管-微血管"功能耦合单元的协同机制对维持人体稳态至关重要。然而目前对该协同机制的理解仍不充分，同步成像与定量表征这些功能耦合单元是阐明皮肤体温调节协同机制不可或缺的方法。当前主要无创检测与成像方法分别聚焦于汗腺导管或血管：通风胶囊法、红外热成像和光学相干断层扫描（OCT）可测量汗腺导管功能或结构，OCT血管成像（OCTA）和光声显微成像（PAM）主要提供血管形态信息。通风胶囊法虽为量化出汗的既定方法，但无法评估皮肤表面及单个汗腺导管的确切功能；中波红外（MWIR）热成像可提供出汗量与汗孔激活信息，但宏透镜景深有限，难以清晰可视化汗孔的细微变化；OCT虽能清晰识别和提取表皮汗腺导管的螺旋结构，但难以同时捕获血管与汗腺导管的高性能图像；OCTA利用固有运动对比区分功能性血管与静态组织背景，但无红细胞血流的皮肤血管无法可视化，且图像质量高度依赖血管阈值选择；PAM结合激光诱导光学照射与超声检测以映射生物组织的光吸收差异，可克服OCTA在血管成像中的问题，具有无创、高空间分辨率及通过血红蛋白和黑色素吸收提供优异内源性对比的优势，能详细可视化体内血管网络与黑色素分布，并依靠超声检测避免严重光学散射，实现微观水平微血管结构的灵敏检测与可视化，但PAM无法提供汗腺导管结构。因此亟需开发能同时成像并量化汗腺导管与血管的新方法。本研究采用新型双模态PA/US显微成像系统，实现指尖"汗腺导管-微血管单元"的同步 mapping 与相关性分析，而非孤立分析这些结构。研究中"超声衍生汗腺结构"特指超声可见的汗腺导管/管状组件，整个解剖学外泌汗腺（包括更深部的盘曲分泌部）未被当前超声成像方法直接分辨。该系统提供SDMU微环境的全面三维体内可视化；为克服生理监测中固有的运动伪影，研究实施了指纹引导的非刚性配准框架，实现对同一单个SDMU从静息状态到运动诱导热应激状态的精确动态追踪；利用该平台首次定量阐明了SDMU的时空耦合动力学，分析揭示功能性微血管耦合高度集中在以汗腺导管为中心的75 μm核心半径内，并发现自主神经调节的显著空间异质性：交感神经系统并非均匀全局激活，而是采用高效的"选择性招募"策略，仅激活约70%的SDMU以满足出汗的局部代谢需求。通过将SDMU确立为功能耦合单元，本研究不仅为皮肤温度微调节机制提供了新见解，还提示SDMU有潜力作为糖尿病周围神经病变早期检测与评估的功能性生物标志物。

研究人员开展的研究：研究招募10名无疾病史的健康参与者，所有程序均详细说明并获得书面知情同意，参与者佩戴防护眼镜，手指涂覆耦合剂置于充水槽的膜密封底部，以一致且最小压力成像，指尖接受575 nm脉冲光照射，扫描步长15 μm，扫描区域6×6 mm2。研究获福建师范大学伦理委员会批准。研究开发了三阶段计算工作流：阶段1为二维基线数据准备，对原始双模态数据进行预处理以提取高保真结构标志，超声图像数据应用拉普拉斯加权多尺度增强算法生成高对比度正面指纹图像，同时利用U-Net模型获得汗腺导管横截面数据及正面形态的最大强度投影（MIP）重建；阶段2为汗孔特征配准，针对静息与运动后生理状态间的空间错位，开发指纹引导的非刚性配准框架，先应用尺度不变特征变换（SIFT）提取并匹配增强指纹图像关键点，再用随机抽样一致性（RANSAC）算法滤除失败匹配的离群点，接着应用薄板样条（TPS）变换计算非刚性形变场，将运动后数据扭曲以实现与运动前状态的完美对齐，经交互式人工验证重叠区域后导出精确配准的汗孔坐标对；阶段3为耦合单元定量分析，利用配准的汗孔坐标将原始三维体积血管数据进行空间映射并加载至基于PyVista的三维交互界面，围绕每个选定汗腺导管构建多半径同心圆柱感兴趣区（ROI），应用局部动态Z轴配准补偿血管舒张引起的轴向位移，最终提取精确配准的微血管体素以定量评估SDMU的空间耦合强度与功能招募。研究采用基于U-Net架构的深度学习网络进行指尖边界与汗腺导管双分割任务，使用2000张人工标注横截面超声图像，随机分为训练集（70%）、验证集（20%）和测试集（10%），采用交叉验证优化超参数，批大小8，训练200轮，使用RMSprop优化器，初始学习率0.0001，衰减率0.995，模型在独立验证集上Dice系数达约98.5%。超声信号主要用于捕获指尖结构信息，在三维重建前应用引导滤波进行图像增强以优化指尖边界与汗腺导管结构可视化。为精确定位指尖微血管网络，需先 delineate 真皮-表皮交界（DEJ），因黑色素主要局限于表皮层，研究利用SSIM定量评估650 nm采集的相邻正面PA切片中黑色素的深度依赖性分布，SSIM能有效捕获强度变化与几何结构相似性，对追踪生物纹理沿深度轴的连续性高度稳健。研究采用多阶段图像处理工作流解决低对比度、信号伪影和背景噪声问题，包括预处理（对比度增强如CLAHE）、血管特征滤波（Frangi滤波器突出管状微血管结构）、分割（Otsu法自动阈值提取二值血管掩膜）、血管骨架化（提取单像素拓扑中心线）和血管描述符（基于骨架化网络进行几何与拓扑分析，用精确距离变换法计算血管直径并以伪彩色编码直径热图呈现）。利用575 nm（P1）和650 nm（P2）采集的原始PA信号可区分血红蛋白与黑色素吸收，计算两波长差值获得主要血管信号，有效抑制表皮纹理中黑色素的干扰。研究基于HSV颜色空间建立多模态图像融合框架，整合超声图像呈现的表皮纹理与汗腺导管形态信息与血管指纹图像中高对比度皮下微血管分布。研究实施超声指纹引导的非刚性配准与验证以实现汗腺导管动态追踪，配准前特征点平均欧氏距离（MED）约18.25像素，TPS非刚性配准后误差显著降至1.38±0.3像素，证明算法在消除运动诱导空间畸变方面具有像素级精度。

研究结果：3.1 运动诱导的汗腺导管与微血管变化——研究首先评估不同生理条件下汗腺导管形态与血管分布变化，运动后由于活跃出汗，分割的汗腺导管呈现显著形态变化，体积、直径和长度均显著增加（p<0.0001），提示超声可见导管组件的运动相关形态响应，可能与活跃出汗期间经外泌汗腺导管通路的瞬时管腔体积变化和汗液运输有关；需注意本研究超声分割主要捕获超声可见汗腺导管结构而非整个解剖学外泌汗腺，尤其是更深部盘曲分泌部。定量血管分析显示运动期间总体血管密度较静息显著增加（P<0.001），而平均血管直径和总血管计数无统计学显著变化（ns），表明运动诱导的指尖微血管反应主要反映为PA可检测的灌注微血管结构表征增加，而非单个血管的均匀扩大。虽然总体定量结果证实热应激下汗腺导管及其周围微血管被激活，但仅评估总体平均值无法确定每个局部汗腺导管结构是否与匹配的微血管响应存在空间耦合，因此研究进一步将单个汗腺导管及其相关血管区域作为配对SDMU进行评估。3.2 SDMU的空间与功能表征——为探究单个汗腺导管与下方微血管的耦合机制，须首先确保这些动态结构在不同生理状态间绝对空间对齐，在超声指纹和汗孔引导的全局配准后，研究引入局部引导的光声血管模板匹配算法，以静息时微血管二值化拓扑特征为参考模板识别运动后匹配血管区域，两个配准步骤有效消除手指颤抖和运动后血管舒张的影响，确保定量分析严格局限于同一SDMU。研究基于积分强度（定义为血管密度×平均信号强度）量化耦合强度并计算运动诱导差异以反映血管响应程度，评估距汗腺导管下方血管区域多个径向距离（R2至R15）的耦合强度差异，分析表明耦合强度差异与半径的关系不呈线性衰减，而是最显著且一致的微血管招募峰值出现在R5半径（约75 μm）附近；进一步评估该最优R5半径下不同汗腺导管的耦合深度显示最强血管响应出现在约1100 μm物理深度。研究在优化空间窗口内对超声可见导管结构直径和体积与相关血管积分强度进行定量相关分析，观察到超声可见导管/管状结构形态变化与优化R5耦合区域内相应微血管响应呈稳健正相关：归一化最大直径变化（r=0.660，p=0.00047）和平均直径变化（r=0.621，p=0.0012）均与功能性微血管招募呈显著正相关，最显著的是超声可见导管结构三维体积变化呈最强相关性（r=0.701，p=0.00013），这种相关强度的逐步增加表明三维体积提供了优化SDMU窗口内导管/管状形态变化的全面超声衍生描述符，R5半径内观察到的显著相关性支持所选空间边界的可靠性。这些发现表明运动后导管形态变化与局部血管招募在优化SDMU内空间相关，支持将这些共配准结构视为成像定义的功能单元。3.3 SDMU的空间异质性与选择性激活——该分类原理源于人体体温调节的既定生理学原则，热应激或运动期间外泌汗腺分泌主要由交感胆碱能激活介导，是散热的主要效应机制之一；同时皮肤血流通过体温调节血管运动响应调节，包括交感控制的活动性血管舒张，增加局部皮肤灌注以促进散热。因此出汗与皮肤血管招募代表体温调节系统两个协调但不同的生理输出；此外出汗响应不仅取决于单个腺体的分泌能力，还取决于皮肤表面功能性活跃汗腺的招募，导致热挑战期间腺体激活的显著空间异质性。基于外泌汗腺激活与皮肤血管调节的生理耦合，研究将SDMU定义为仅当两个成像衍生标准同时满足时才表现协调激活响应：第一，超声可见汗腺导管形态归一化变化大于零（Δ>0），表明与汗腺导管活动增加相关的局部结构扩张；第二，优化75 μm径向边界内PA衍生微血管耦合强度变化大于零（Δ>0），表明局部血管招募增强；未满足任一标准的SDMU被分类为不符合双响应标准。基于此严格定义，研究对单个指尖29个配对SDMU进行纵向分析，观察到激活过程的显著空间异质性：约70%（20/29）SDMU为协调激活，其余30%（9/29）保持休眠；激活单元表现同步血管舒张和稳健局部微血管灌注，而未激活SDMU几乎无形态或血管响应。代表性三维体积可视化显示，激活SDMU特征为汗腺导管（蓝色）增大和集中于汗腺导管周围密集高强度微血管网络（红黄色），而未激活SDMU汗腺导管体积、血管密度或光声强度无显著变化。这些结果证实热应激期间SDMU并非被交感神经系统均匀激活，而是采用选择性招募的高效策略，仅激活空间特定子集SDMU以满足局部体温调节需求，同时节约全身血容量和代谢能量，这种异质性调节机制是维持皮肤热稳态的基础。

讨论部分总结：本研究采用指纹引导双模态PA/US显微成像系统系统研究运动诱导热应激下人指尖SDMU的空间 mapping 与功能相关性，发现不仅证实热应激下汗腺导管与周围微血管的协同激活，还将SDMU确立为统一的空间耦合功能单元，该耦合单元有潜力作为糖尿病神经病变早期检测与评估的新型定量生物标志物。4.1 运动诱导热应激：汗腺导管与微血管的总体响应与机制——与既往皮肤体温调节和汗腺导管形态研究一致，结果表明分割的超声可见汗腺导管/管状结构运动后发生显著形态变化，特征为等效直径、长度和体积增加（P<0.0001），支持超声衍生导管/管状指标对运动诱导生理应激的敏感性，但不应解释为全腺体增大或分泌盘曲扩张的直接证据；更重要的是血管网络密度显著增加（P<0.001），而平均微血管直径或计数无明显差异，提示运动诱导指尖充血主要反映为灌注微血管结构招募或表征增加而非现有血管均匀扩张。解释超声衍生汗腺结构指标时需注意解剖学区分：外泌汗腺导管包括更深部盘曲分泌部和向皮肤表面延伸的导管组件，分泌盘曲主要负责汗液产生，而导管通路将分泌液传导至皮肤表面并通过导管修饰和重吸收参与最终汗液组成；因此图10量化的体积、等效直径和长度应解释为超声可见导管/管状结构指标而非整个外泌汗腺或分泌盘曲的直接测量，观察到的运动后变化可能反映导管/管状变形、管腔填充或活跃出汗伴随的局部组织几何变化，但当前超声数据集无法区分这些机制或独立量化分泌盘曲扩张；未来工作将进一步优化PA/US平台以实现更深部外泌汗结构（特别是位于真皮深部的分泌盘曲）的深度分辨可视化及共配准微血管评估，以更好确定运动诱导超声可见导管变化与活跃出汗期间全汗腺响应和局部血管招募的关系。4.2 SDMU：具有确定空间边界和强相关激活的统一功能单元——本研究关键创新在于利用双模态成像系统和超声指纹引导的非刚性配准框架将SDMU确立为独立功能单元，局部光声血管模板匹配算法有效消除手指颤抖和运动后血管舒张引起的配准误差，确保定量分析仅在静息和热应激条件下对同一SDMU进行，这是研究其空间与功能耦合的前提。空间边界分析确定微血管耦合最优径向距离约为距汗腺导管/管状结构75 μm（R5半径），此处微血管招募最显著且一致；此外最强血管响应出现在约1100 μm物理深度，界定SDMU分析所用垂直成像范围。这些发现应解释为超声可见导管/管状变化与PA可检测微血管招募关联最强的优化成像窗口，而非整个外泌汗腺或分泌盘曲的直接解剖边界；1100 μm深度与解剖学和三维重建研究报告的外泌汗腺结构真皮深部位置相符，但当前PA/US测量未独立分割分泌盘曲。定量相关分析揭示超声可见导管/管状形态参数（最大直径、平均直径和三维体积）与相应微血管区域积分强度呈显著正相关，体积与积分强度相关性最强（r=0.701，p=0.00013），表明三维体积是导管/管状形态响应的敏感超声衍生描述符，因其整合直径和长度变化；更重要的是这些相关性表明导管/管状形态变化与局部微血管招募在优化SDMU窗口内空间相关，支持SDMU作为共配准成像定义的功能单元，同时避免整个分泌盘曲被直接可视化或出汗量被定量测量的更强解释。4.3 空间异质性与选择性招募：概念验证评估——基于已建立的空间边界，研究引入定量双阈值评估方法系统评估单个SDMU功能状态；作为概念验证，应用该评估框架纵向追踪代表性受试者29个SDMU，揭示热应激条件下显著空间异质性，发现约70%（20/29）追踪SDMU被选择性激活，表现超声可见导管/管状形态和微血管灌注协调增加，其余30%（9/29）保持休眠，形态和血管变化极小；这一局部观察提示运动诱导热应激期间单个SDMU并非均匀响应，而是呈现"按需"招募模式，仅部分成像定义SDMU显示协调导管/管状和血管响应；这种异质性可能归因于导管/管状敏感性、微血管密度或交感神经支配的局部差异，与既往汗腺导管招募和出汗响应在皮肤区域和个体间存在差异的观察一致；从功能角度看，SDMU选择性激活可能实现散热的空间分辨调节。4.4 局限性、临床意义与未来方向——尽管本研究建立了稳健方法学框架，仍有若干局限性需考虑：第一，虽然总体生理响应在参与者队列中评估，但详细单单元异质性评估集中于单个代表性受试者的29个SDMU，因此观察到的70/30激活比例应解释为探索性单参与者观察而非SDMU激活的群体水平估计；该局限性尤为显著，因既往研究报告促汗响应和皮肤血管响应均存在显著个体间差异；未来研究需更大更多样化的独立参与者队列、重复测量和受试者水平统计模型以量化个体内和个体间变异性、验证生理学信息驱动的操作性双响应标准并确定SDMU衍生指标的可推广性。从临床转化角度看，所提PA/US框架应视为对评估出汗和皮肤微结构的既定方法的补充：OCT可实现汗腺导管无创高分辨率可视化并已用于外泌汗腺出汗的结构和动态评估，但传统结构OCT主要依赖光学散射或折射率对比，汗腺导管形态与富血红蛋白真皮微血管的同步共配准评估仍具挑战；通风胶囊技术是热或运动应激期间连续测量局部出汗率的既定方法，但提供面积平均蒸发输出，无法评估单个汗腺导管形态变化，也无法确定这些单元如何与周围微血管区域空间耦合。与这些方法相比，当前PA/US框架的主要优势在于能在同一三维坐标系统内整合超声衍生汗腺导管形态与PA衍生微血管信息，这种共配准双对比策略允许将SDMU作为局部功能单元评估而非孤立的导管形态或整体出汗输出，该能力可能与神经病变相关评估特别相关，因为糖尿病和自主神经病变可同时涉及促汗功能障碍和微血管损伤；通过同时量化汗腺导管/管状形态和局部微血管招募，PA/US框架可能提供比单独OCT结构成像或通风胶囊出汗率测量更具空间分辨的神经血管和促汗耦合读数；然而由于本研究在健康参与者中进行，SDMU衍生指标的临床诊断价值仍有待专门患者队列验证，临床转化需要疾病特异性队列的前瞻性验证、重测信度和操作者间稳健性评估以及临床可解释SDMU响应标准的建立。原发性局灶性多汗症特征为超出体温调节需求的过度出汗，通常影响手掌、足底、腋窝或颅面部区域，可显著损害患者生活质量，但其诊断和严重程度评估仍具挑战，因当前评估主要依赖病史和区域出汗输出测试，尽管这些方法提供有用临床证据，但出汗发作可能间歇性和情绪触发，而Minor淀粉-碘测试主要提供出汗区域的定性或半定量空间 mapping 而非客观单导管功能信息。所提指纹引导双模态PA/US显微成像框架未来可扩展至多汗症队列，该扩展在生理学上合理，因原发性多汗症与交感胆碱能促汗活动过度相关，交感胆碱能神经释放的乙酰胆碱是外泌汗腺分泌的主要刺激；此外活跃出汗伴随皮肤血流和局部微血管灌注增加，可能为分泌活动和散热提供血管支持；在此背景下共配准PA/US成像可实现手掌或指尖区域单个SDMU变化和整体SDMU激活模式的监测与量化，从而帮助确定过度出汗是否与显著汗腺导管响应、局部微血管耦合改变或SDMU激活异常空间异质性相关，这些特征可能为未来多汗症评估提供潜在定量指标；此类应用仍属探索性，需针对既定临床评估、出汗输出测量和Minor淀粉-碘测试进行验证。

研究结论翻译：总之，研究人员开发了指纹引导双模态PA/US显微成像平台，能够在体内精确 mapping 和量化SDMU的空间-功能动力学。结果表明运动诱导热应激引起同步神经血管响应，揭示皮肤充血主要由靶向微血管招募而非全局血管舒张驱动。通过在高定位75 μm径向边界和1100 μm轴向深度内定量 delineate SDMU，研究验证了超声可见导管/管状形态扩张与局部血管激活之间的稳健正相关。最终，本研究将SDMU从概念模型转变为严格定义的成像功能单元，加深了对皮肤体温调节的基本理解；此外该定量范式可能为未来评估微血管和自主神经异常（特别是糖尿病周围神经病变）提供无创框架。