结节性硬化症儿童脑微结构改变：基于神经突方向离散度与密度成像的研究

《Brain Sciences》：Brain Microstructural Alterations in Children with Tuberous Sclerosis Complex: A Neurite Orientation Dispersion and Density Imaging Study

　　

目的：应用神经突方向离散度与密度成像（NODDI）评估结节性硬化症（TSC）儿童的脑微结构改变，并探讨NODDI参数与临床表型（包括癫痫、早发性癫痫、婴儿痉挛性癫痫、癫痫持续状态及药物难治性癫痫）之间的关联。方法：本前瞻性横断面研究纳入68例TSC儿童和39名健康对照（性别匹配；两组年龄差异显著，所有分析均对年龄进行调整，并采用样条建模及儿科范围敏感性分析），所有受试者均接受多壳扩散磁共振成像（MRI）检查。计算NODDI参数：神经突密度指数（NDI）、方向离散度指数（ODI）及各向同性体积分数（FISO）。采用Mann–Whitney U检验评估组间差异，并进一步行年龄调整分析。采用逻辑回归分析NODDI参数与临床变量之间的关联。采用受试者工作特征（ROC）分析评估推导队列内的判别性能，并进行bootstrap及交叉验证的内部验证。结果：经年龄调整后，TSC组全脑NDI、ODI及FISO较对照组显著降低（均p ≤ 0.001）。合并癫痫的TSC患者较对照组NDI（p < 0.001）、ODI（p = 0.001）及FISO（p = 0.003）显著降低，而无癫痫患者与对照组之间无显著差异。早发性癫痫及婴儿痉挛性癫痫与NDI显著降低相关（分别降低16.7%和19.0%；p < 0.001）。在简约模型中，NDI与癫痫相关（p = 0.008；bootstrap校正AUC 0.712）；在探索性多变量模型中，年龄调整后NDI（p = 0.008）和ODI（p = 0.035）与癫痫相关。NDI的AUC为0.739（95% CI 0.617–0.861），截断值为0.095时敏感度为65.5%，特异度为78.8%，反映推导队列内的判别能力。结论：NODDI可灵敏检测TSC儿童广泛的脑微结构改变。NDI以及调整后程度较弱的ODI与癫痫存在显著关联，可作为癫痫相关微结构负荷的定量影像标志物；然而，横断面设计无法确定这些改变是由癫痫所致还是反映TSC underlying疾病过程，其是否可预测未来癫痫发作或结局尚需纵向、外部验证研究。两组间显著的年龄差异为统计分析无法完全解决的关键局限性。NODDI宜视为多模态生物标志物框架的组成部分之一，该框架还应整合分子、免疫及电生理指标。

结节性硬化症（TSC）是一种罕见的常染色体显性遗传性神经皮肤综合征，发病率约为1/6000至1/10000活产儿，由TSC1或TSC2基因致病性变异引起，导致哺乳动物雷帕霉素靶蛋白（mTOR）信号通路组成性激活，进而破坏细胞增殖、分化与迁移，在多器官系统形成错构瘤。中枢神经系统受累是TSC发病和死亡的主要来源，癫痫影响高达85%的患者，三分之二受累儿童发展为药物难治性癫痫（DRE），约90%的TSC儿童存在不同程度的TSC相关神经精神障碍（TAND）。常规结构MRI（T1加权、T2加权及液体衰减反转恢复序列）虽能可靠显示皮质结节、放射状迁移线、室管膜下结节及室管膜下巨细胞星形细胞瘤等特征性病变，并有研究建立了结节大小、钙化/囊性外观、结节-脑比值、局灶性皮质发育不良样特征及放射状迁移线与致痫性的关联，还提出了基于信号强度和形态特征的结节分型系统，但结构MRI仅提供宏观形态学信息，无法捕捉可能同等或更为关键的微结构改变。扩散张量成像（DTI）虽已用于研究TSC白质微结构完整性，发现持续癫痫儿童存在广泛的各向异性分数（FA）降低和平均扩散率（MD）改变，且微结构异常可延伸至可见结节之外，但DTI受限于单室模型及解析复杂纤维构型的能力有限，且FA在某些病理条件下可能因特定纤维群选择性变性而 paradoxically 不变甚至升高。这些局限促使研究者开发能提供更具特异性和生物学可解释性的组织微结构测量技术的先进扩散MRI方法。NODDI采用多室生物物理模型，区分神经突内受限扩散、神经突外阻碍扩散及脑脊液自由扩散三种微结构环境，生成细胞内体积分数（ICVF，亦称NDI，反映神经突密度）、ODI（量化神经突角度离散度）及FISO（代表自由水比例）三个参数图，具有对微结构改变特异性更高、在复杂纤维构型区域敏感性更强、能区分神经突密度与方向离散度相对贡献等优势，已应用于神经发育、阿尔茨海默病、帕金森病、多发性硬化、癫痫及局灶性皮质发育不良等多种神经系统疾病。NODDI在TSC中的早期应用已取得初步结果：Taoka等首次证明TSC患者白质NDI降低，尤其在含皮质结节的联络纤维束中，且NDI对TSC相关白质改变的敏感性高于常规DTI参数；Shao等进一步显示NODDI参数（尤其NDI）在勾画致痫结节方面较常规MRI和DTI提供更优对比度；Anaby等报告全脑白质NDI与临床严重程度评分呈显著负相关。组织病理学证据显示，巨细胞、巨神经元及星形胶质细胞增生不仅存在于结节内，也存在于周围MRI表现正常的组织中，为观察到的微结构改变提供了生物学基础。然而，现有研究存在若干关键知识空白：样本量相对较小，无法进行基于临床表型的稳健亚组分析；NODDI参数与特定临床表现（尤其癫痫特征，包括发作类型、发病年龄、药物难治性及癫痫持续状态）之间的关系尚未系统研究；NODDI参数对识别TSC患者癫痫及不良神经结局风险的预测价值尚未确立；遗传病理与继发性癫痫活动对观察到的微结构改变的相对贡献尚不完全清楚。本研究旨在应用NODDI全面评估TSC儿童的脑微结构改变，具体目标为：（1）量化TSC患者与健康对照之间全脑NDI、ODI及FISO差异；（2）探讨NODDI参数与不同临床表型（包括癫痫、早发性癫痫、婴儿痉挛性癫痫、癫痫持续状态及药物难治性癫痫）之间的横断面关联；（3）评估NODDI参数识别合并癫痫的TSC患者的判别性能。研究人员假设TSC患者存在广泛的微结构改变，且NODDI参数与癫痫表型相关并可区分合并癫痫的患者。本研究为前瞻性横断面研究，于2025年3月至12月在深圳市儿童医院（中国广东省罕见神经系统疾病三级转诊中心）开展，研究方案经深圳理工大学伦理审查委员会批准，所有程序严格遵循《赫尔辛基宣言》原则，并获得所有参与者父母或法定监护人的书面知情同意，≥7岁儿童在适当时获取知情同意，报告遵循STROBE指南。筛查116名儿童，排除4名临床数据不完整者和5名图像质量不佳者，最终分析样本为68例TSC儿童和39名健康对照；质量控制后，67例TSC患者和38名对照可用于TBSS，66例TSC患者和37名对照可用于GBSS。TSC患者队列共68例（男34例，女34例；中位年龄8.5岁，IQR 7–16岁；范围0.58–16岁），连续纳入自儿童神经内科门诊和住院部，由儿童神经科医生、遗传学家和神经放射学家组成的多学科团队根据2021年国际结节性硬化症共识会议诊断标准确立诊断，所有患者接受全面临床评估，包括详细神经系统检查、基因检测（如有）、脑电图监测及系统器官评估。纳入标准：（1）根据2021年国际标准确诊TSC；（2）发病年龄<18岁；（3）完整临床数据和脑MRI DICOM数据；（4）法定监护人书面知情同意。排除标准：（1）既往癫痫手术（包括迷走神经刺激或胼胝体切开术等姑息性手术）；（2）脑炎、颅内感染或其他混淆影像学解释的神经系统疾病史；（3）MRI序列不完整（T2加权、FLAIR或多壳扩散MRI）；（4）图像质量差或运动伪影过大（平移>3 mm或旋转>3°）无法可靠拟合NODDI。健康对照队列共39名（男21名，女18名；中位年龄11岁，IQR 10–14.58岁；范围6–18岁），招募自当地社区，仅按性别与TSC患者组匹配；由于招募限制未实现年龄匹配，由此产生的显著年龄差异在所有分析中进行统计学处理。对照纳入标准：（1）无个人或家族神经或精神疾病史；（2）无物质滥用或伴意识丧失的头部外伤史；（3）神经系统检查正常；（4）常规脑MRI（T1加权、T2加权和FLAIR）由两名神经放射学家独立确认无异常发现；（5）法定监护人书面知情同意。临床数据通过电子病历系统提取并通过结构化访谈补充，使用标准化病例报告表记录性别、MRI采集时年龄、TSC诊断时年龄、TSC1/TSC2基因突变状态、TSC/癫痫/发育障碍家族史、癫痫发病年龄、发作类型（按2017年ILAE操作分类）、抗癫痫药物（ASM）方案、既往和当前非药物治疗（包括生酮饮食、mTOR抑制剂或癫痫手术）、发作频率和预后、神经精神共病及根据2021年监测建议的其他器官系统受累情况。癫痫相关变量严格按既定国际标准定义：癫痫为至少两次间隔>24小时的非诱发性发作，或一次发作且未来10年复发风险高（≥60%）；早发性癫痫为首发年龄<1岁；婴儿癫痫性痉挛综合征（IESS）按2022年ILAE分类标准诊断，要求：（1）发病年龄1–24个月；（2）典型屈曲性、伸展性或混合性癫痫性痉挛成簇出现；（3）发作间期EEG显示高峰节律紊乱或改良高峰节律紊乱；（4）可能的发育倒退（第1–2项为强制）；癫痫持续状态（SE）为持续发作活动≥30分钟，或反复发作且发作间期意识未恢复至基线；药物难治性癫痫（DRE）为按2010年ILAE定义，在最大耐受剂量下充分试用≥2种适当ASM后仍未达到持续无发作（≥12个月或≥3倍最长治疗前无发作间期）。所有参与者接受3.0 Tesla Siemens Prisma系统多壳扩散MRI检查，使用64通道头颈线圈，泡沫垫和胶带固定以减少头动，无法配合的儿童按机构镇静方案口服水合氯醛（50 mg/kg，最大1 g）或静脉右美托咪定。采用自旋回波单次激发平面回波成像（EPI）序列，参数：TR = 12,500 ms；TE = 81 ms；FOV = 224 × 224 mm2；采集矩阵 = 114 × 114；层厚 = 2.0 mm；80层连续轴位；各向同性体素 = 2.0 × 2.0 × 2.0 mm3。多壳扩散编码方案包括83个非共线扩散加权方向，分布于四个b值：b = 0 s/mm2（2卷）、b = 1000 s/mm2（27方向）、b = 2000 s/mm2（27方向）、b = 3000 s/mm2（27方向）。并行成像采用blipped CAIPIRINHA层面加速与平面内GRAPPA（AF = 4）组合。采集一个反向相位编码方向的b0卷用于畸变校正。总采集时间约17分55秒。扩散MRI数据使用FSL（6.0.7.4版）、MRtrix3（3.0.4版）和ANTs（2.3.5版）实现的综合流程预处理，所有处理步骤在原生扩散空间进行，各阶段严格质量控制，完整预处理流程见图2，描绘从原始DICOM转换经伪影校正到NODDI参数估计的序列级联。原始DICOM文件使用dcm2niix转换为NIfTI格式，生成扩散加权图像（DWI）、b向量和b值。使用b0参考卷通过FSL脑提取工具（BET）以0.3分数强度阈值进行脑提取生成脑掩膜，掩膜经视觉检查并在必要时手动校正。应用综合伪影校正级联：噪声校正使用MRtrix3的dwidenoise以5体素滑动窗口实现MP-PCA去噪，去除热噪声而不引入空间模糊；Gibbs振铃去除使用MRtrix3的mrdegibbs，采用外推方法抑制组织边界截断伪影；磁化率畸变校正使用FSL的topup，利用相反相位编码方向采集的b0图像校正B0场不均匀引起的几何畸变；运动与涡流校正使用FSL的eddy同时校正受试者运动、涡流诱导畸变和信号丢失，启用层到体积配准以考虑体积内运动，最终输出包括去扭曲DWI卷和重定向b向量；偏置场校正使用ANTs的N4BiasFieldCorrection校正射频线圈不均匀性引起的强度不均匀。预处理后使用NODDI MATLAB工具箱（1.05版）估计NODDI参数，该工具箱实现多室生物物理模型。NODDI模型将扩散加权信号分解为三个室：细胞内室（神经突内受限扩散）、细胞外室（阻碍扩散）及各向同性室（脑脊液自由扩散），总信号模型结合这些室。模型参数使用非线性最小二乘优化逐体素估计，目标函数最小化残差平方和，优化从10个随机起始点初始化以确保全局收敛。生成三个参数图：NDI（神经突体积分数，范围0–1）、ODI（神经突角度离散度，范围0–1）、FISO（自由水分数，范围0–1）。全脑感兴趣区分析计算脑实质内所有体素的平均NDI、ODI和FISO值，将BET脑掩膜应用于每个参数图，计算所有幕上和幕下区域的体素级均值。由于BET提供的是实质掩膜而非组织特异性分割，全脑均值代表灰质、白质和残余部分体积/脑脊液体素的复合，并非纯组织估计；FISO尤其易受脑脊液部分体积效应影响。为评估这一点，敏感性分析在排除FISO > 0.5的体素后重新计算全脑NDI、ODI和FISO，从而去除自由水主导的体素。组织特异性推断由下述GBSS（灰质）和TBSS（白质）分析提供。白质微结构差异使用FSL中的基于纤维束的空间统计（TBSS）评估，每个受试者的FA图像使用FNIRT非线性对齐至FMRIB58_FA标准空间，使用FA阈值0.2创建平均FA骨架；关键的是，全程使用FA而非NDI进行配准和骨架构建。然后将NDI、ODI和FISO值通过垂直于局部极大值搜索投影到该共同FA衍生骨架上，对所有三个指标使用相同的变形场和投影向量，确保所有NODDI参数在相同解剖位置采样。骨架化图使用FSL的Randomise（5000次置换）分析，采用TFCE和FWE校正，p < 0.05。灰质微结构改变使用定制GBSS流程评估：第一阶段组织分割与模板创建，组织分数估计中自由水分数由NODDI衍生的FISO图近似，白质分数通过ANTs的Atropos分割FA图获得，灰质分数计算为余量；研究人员注意到FISO是模型衍生的自由水参数而非常规脑脊液组织概率图，因此灰质骨架锚定于基于模型而非结构组织的估计；为评估这一依赖性，验证分析使用ANTs的Atropos（三分类分割）从T1加权图像独立衍生脑脊液、灰质和白质分数，并重复模板、骨架和投影步骤，灰质骨架及NDI和ODI降低的空间分布和方向基本不变（显著体素空间重叠91.4%；峰值位移0.8 mm），因此主要GBSS结果报告使用基于FISO/FA的分数以与NODDI框架一致，结构MRI实现作为确认。伪T1图像生成，使用ANTs的antsMultivariateTemplateConstruction2.sh（ANTs 2.3.3版）以对称微分同胚配准（SyN）创建人群模板；非线性变形，反转变形场并应用于将灰质分数、NDI和ODI图变形至模板空间。第二阶段灰质骨架化：平均灰质分数图像，所有模板空间灰质分数图像平均化；骨架化，平均灰质图像骨架化以代表皮质带中厚度；骨架阈值化，保留灰质分数>0.65且存在于>75%参与者的体素；为限制部分体积和误分类体素的影响，推断限于灰质分数>0.65且存在于>75%参与者的体素，但这些标准并未完全消除对基于模型的组织估计的依赖；参数投影，将NDI和ODI值从灰质分数最大的体素投影到骨架。第三阶段骨架填充与统计分析：骨架填充，未通过灰质分数阈值的体素使用高斯加权平均（σ = 2 mm）填充；统计分析，使用FSL的Randomise（FSL 6.0.4版；5000次置换）进行体素级比较，采用TFCE和FWE校正，p < 0.05。所有统计分析使用SPSS（24.0版）和R（4.2.1版）进行，双侧p < 0.05认为具有统计学显著性。人口学和临床特征中，连续变量正态分布数据以均值±SD表示，非正态分布以中位数（IQR）表示；分类变量以频率（百分比）表示；组间比较酌情使用独立样本t检验、Mann–Whitney U检验或卡方检验。NODDI参数组间比较中，TSC vs. 对照使用Mann–Whitney U检验比较全脑平均NDI、ODI和FISO。亚组分析为避免表型特异性差异与癫痫vs.无癫痫的更广泛区分相混淆，所有表型比较限于55例合并癫痫的患者；在该癫痫限制队列内，将各表型患者与无该表型的癫痫患者进行比较：早发性癫痫（28 vs. 27）、IESS（13 vs. 42）、SE（13 vs. 42）和DRE（24 vs. 31）；13例无发作患者从这些对比中排除。原始全队列对比中比较组包括13例无发作患者，仅作为假设生成分析。由于这些比较为探索性，在每个NODDI参数内跨四个癫痫限制表型对比分别应用Benjamini–Hochberg FDR校正（每参数4次检验；共12次检验），报告未校正p值和FDR校正q值。主要TSC vs. 对照比较预先指定为每参数单一假设，未纳入FDR校正；描述性三组比较（对照 vs. 无癫痫TSC vs. 合并癫痫TSC）同样未校正，仅作描述性呈现。无论显著性如何，发现均解释为假设生成。鉴于组间年龄差异显著，NODDI参数的组间比较额外使用ANCOVA对年龄进行调整；结果报告包含和不包含年龄调整。ANCOVA前验证模型假设，包括残差正态性（Shapiro–Wilk检验，p > 0.05）和回归斜率同质性（年龄×组交互项，均p > 0.05）。为考虑潜在非线性发育效应，敏感性分析中额外使用限制性立方样条（3节点）建模年龄。还进行了限于重叠儿科年龄范围（15–18岁）参与者的敏感性分析以最小化残余年龄混杂。由于统计调整无法完全解释组间发育差异，使用三种互补方法：第一，使用限制性立方样条（3节点）非线性建模年龄并与线性模型比较；第二，进行限于重叠儿科年龄范围的敏感性分析，由于对照组从6岁起招募而TSC组低至0.58岁，两组均代表的年龄范围为6–18岁，因此分析限于6–16岁参与者（TSC n = 68；对照 n = 31），排除4名>16岁对照和最年轻TSC患者，同时保留完整TSC队列；还进行了仅排除>18岁对照（n = 4）的单独分析；在15–18岁窄带（两组均代表，TSC n = 4；对照 n = 4）进行更窄的描述性分析，结果中仅作描述性报告；第三，描述性年龄分层分析在两组均代表的年龄带（<6、6–10、10–15、>15岁）内比较TSC患者和对照的NODDI参数。这些分析旨在评估组间差异对残余发育混杂的稳健性，而非消除它。多变量逻辑回归识别癫痫的独立预测因子，结果以比值比（OR）及95%置信区间（CI）表示。鉴于无癫痫参与者有限（n = 13），主要模型刻意简约：（i）单变量NDI，（ii）NDI + 年龄（两个预测因子，EPV ≈ 6.5）。通过bootstrap重采样（1000次重复）评估稳定性，报告乐观校正AUC和95% CI。四预测因子模型（NDI、ODI、FISO、年龄）仅作为探索性分析拟合，因每变量事件比≈3.25低于常规稳定性阈值；其估计值附带此警示报告。使用Spearman相关矩阵和方差膨胀因子（VIF > 5视为有问题）评估多重共线性。ROC分析评估推导队列内的判别性能，不构成验证性预测；AUC报告95%置信区间（DeLong法）。内部验证使用bootstrap乐观校正（1000次重采样）和5折交叉验证，报告乐观校正AUC。最佳截断值由Youden指数确定，承认为样本特异性。无外部验证队列可用。各阶段实施严格质量控制：图像质量由两名独立神经放射学家检查所有数据以发现运动伪影和畸变，平移≥3 mm或旋转≥3°的数据集排除；预处理QC系统审查脑体积、配准质量、运动参数和离群值计数；NODDI拟合QC计算测量与模型信号之间的RMSE；配准QC视觉检查变形图像，Jacobian行列式>0.3确认有效变形；数据排除最终样本量：全脑分析N = 68（TSC）和N = 39（对照）；TBSS N = 67和N = 38；GBSS N = 66和N = 37。样本量使用G*Power 3.1针对双尾独立t检验确定（α = 0.05，power = 0.80，效应量d = 0.6），每组需45名参与者；考虑脱落（10%）和技术失败（15%），目标≥70例TSC患者和≥40名对照。最终68例TSC患者和39名对照的样本提供80% power检测0.5 SD的差异。对于逻辑回归，55个事件（80%癫痫率）满足每预测变量≥10个事件的标准。结果方面，共纳入68例TSC儿童（男34例，女34例；各50%）和39名健康对照（男21名，女18名；分别为53.8%和46.2%）。TSC组平均年龄8.1 ± 3.93岁（中位8.5岁，IQR 7–16岁；范围0.58–16岁），对照组平均年龄13.09 ± 5.72岁（中位11岁，IQR 10–14.58岁；范围6–18岁）。组间年龄差异显著（p < 0.001），对照组年龄较大。该年龄差异为潜在混杂因素，所有后续分析酌情对年龄进行调整。任何NODDI参数均未观察到显著的年龄×组交互（均p > 0.10），限制性立方样条建模确认主要组效应独立于非线性年龄效应。排除>18岁对照（n = 4）的敏感性分析中，结果基本不变（NDI：p < 0.001；ODI：p = 0.001；FISO：p = 0.002）。在限于6–16岁参与者（TSC n = 68；对照 n = 31）的进一步敏感性分析中，TSC患者仍显示NDI（p < 0.001）、ODI（p = 0.002）和FISO（p = 0.004）较对照显著降低，且跨癫痫亚组的梯度模式得以保留。在15–18岁较窄重叠范围（两组均代表，TSC n = 4；对照 n = 4）中，组间差异方向与主要分析一致，但样本量太小无法进行有意义的统计推断，因此仅作描述性报告。年龄分层描述性分析中，TSC患者全脑NDI和ODI在每个两组均代表的年龄带均低于对照（6–10岁：NDI p = 0.012，ODI p = 0.021；10–15岁：NDI p = 0.008，ODI p = 0.019），表明组间差异并非仅归因于不同的年龄分布。然而，由于两组未进行年龄匹配且NODDI参数在整个儿童期和青春期持续变化，残余发育混杂无法排除。基因筛查结果显示，在54例接受基因检测的患者（54/68，79.4%）中，43例（占检测者79.6%）携带TSC2变异，8例（14.8%）携带TSC1变异，3例（5.6%）未检测到变异。队列组成显示55例（80.9%）合并癫痫，13例（19.1%）无发作。55例癫痫患者的发作类型分布：局灶性发作48例（87.3%）、早发性癫痫28例（50.9%）、婴儿癫痫性痉挛综合征（IESS）13例（23.6%）、癫痫持续状态（SE）13例（23.6%）、药物难治性癫痫（DRE）24例（43.6%）。全脑NODDI参数比较中，TSC组全脑NDI、ODI和FISO值较健康对照显著降低：NDI（p < 0.001）、ODI（p < 0.001）、FISO（p = 0.001）。NDI降低表明神经突密度降低，提示轴突和树突丢失；ODI降低反映方向离散度降低，意味着纤维复杂性降低；FISO降低提示自由水分数降低，反映组织组成改变。这些发现表明TSC存在广泛的全脑微结构改变；这些改变是否延伸至可见皮质结节之外的组织，本分析未直接检验。这些全脑差异并非归因于脑脊液部分体积污染：排除FISO > 0.5的体素后，TSC患者仍显示NDI（p < 0.001）、ODI（p = 0.001）和FISO（p = 0.004）降低，效应方向和幅度基本不变。GBSS在灰质中观察到一致的NDI和ODI组织特异性降低，而TBSS在FWE校正后未显示显著白质差异。TSC合并癫痫（n = 55）与无癫痫（n = 13）患者之间任何NODDI参数均无显著差异（NDI：0.088 ± 0.017 vs. 0.096 ± 0.014，t = 1.657，p = 0.102；ODI：0.102 ± 0.013 vs. 0.104 ± 0.012，t = 0.533，p = 0.595；FISO：0.032 [0.028–0.035] vs. 0.032 [0.028–0.036]，Z = ?0.700，p = 0.944）。然而，三组间呈现清晰的梯度模式：TSC合并癫痫参数值最低，其次为无癫痫TSC，健康对照最高。与健康对照相比，TSC合并癫痫患者NDI（p < 0.001）、ODI（p = 0.001）和FISO（p = 0.003）显著降低，而无癫痫TSC患者无显著差异。该梯度提示underlying疾病过程是微结构改变的主要驱动因素，癫痫代表更严重病理的标志。TSC合并与不合并癫痫患者之间无统计学显著差异应结合无癫痫亚组样本小（n = 13）、该比较未调整及亚组间年龄差异大来解释；下述回归分析将NDI和ODI作为连续变量并调整年龄建模，因此检测关联的敏感性更高。由于仅13例患者无癫痫，该比较把握度不足，无显著差异应视为不确定而非等效证据。三组间梯度模式（对照 > 无癫痫TSC > 合并癫痫TSC）因此作为描述性和假设生成呈现，并不确立癫痫导致微结构改变。在55例癫痫患者中，早发性癫痫患者（n = 28）较无早发性癫痫的癫痫患者（n = 27）NDI（p < 0.001；降低16.7%）和ODI（p = 0.007；降低8.5%）显著降低。IESS患者（n = 13）显示最显著降低：较无IESS的癫痫患者（n = 42）NDI（p < 0.001；降低19.0%）和ODI（p = 0.006；降低9.4%）。这些发现表明在关键发育窗口出现的严重癫痫表型与神经突密度和组织结构更大程度降低相关；由于横断面设计，无法确定该关联的方向。在次要全队列对比中，比较组包括13例无发作患者，效应量相似但略大，符合纳入相对保留亚组于比较组时的预期；这些次要分析不应解释为表型特异性效应。在55例癫痫患者中，SE患者（n = 13）与无SE患者（n = 42）之间任何NODDI参数均无显著差异；观察到NDI（p = 0.068，q = 0.204）和ODI（p = 0.137，q = 0.274）降低趋势，提示SE可能贡献额外神经突丢失，但未通过校正。DRE患者（n = 24）较无DRE的癫痫患者（n = 31）FISO显著升高：0.034（IQR 0.030–0.039）vs. 0.031（IQR 0.027–0.033）；Z = ?2.800，p = 0.005，q = 0.020。NDI（p = 0.799，q = 0.799）或ODI（p = 0.558，q = 0.744）无显著差异。FISO升高可能反映自由水含量增加，包括药物难治性病例中胶质增生相关组织重塑、细胞外间隙改变或血脑屏障通透性改变。在次要全队列分析中DRE–FISO关联未通过FDR校正（q = 0.084）；因此此处报告的癫痫限制结果为初步发现，尽管鉴于亚组样本量仍属探索性，需在独立队列中验证。TBSS分析在FWE校正后未显示TSC患者与健康对照之间统计学显著的白质差异。相反，GBSS分析显示显著的灰质改变。GBSS分析显示TSC组与健康对照之间NDI组间差异的轴位切片，显著区域（黄色，FWE校正，p < 0.05）主要分布于皮质区域，包括额叶（前额叶皮质、辅助运动区）、顶叶（顶上小叶、中央后回）、颞叶（颞上回、颞极）和枕叶（视觉皮质、楔叶）区域。ODI组间差异的轴位切片显示相似皮质区域显著ODI降低（FWE校正，p < 0.05），表明涉及神经突密度降低和方向离散度降低的广泛灰质微结构改变。皮质优势与TSC已知组织病理学特征一致，包括主要影响灰质的皮质结节形成和异常神经元迁移。所有三个NODDI指标均投影到共同FA衍生白质骨架；NDI未用于配准或骨架构建。基于T1的GBSS流程重新实现产生基本一致的灰质骨架及NDI和ODI降低的空间模式（显著体素空间重叠91.4%；峰值位移0.8 mm），表明发现并非基于FISO组织分数估计的伪影。逻辑回归分析中，在探索性四预测因子模型（EPV ≈ 3.25）中，调整后NDI（B = ?0.154，SE 0.058，p = 0.008，OR 0.858，95% CI 0.765–0.961）和ODI（B = +0.166，SE 0.079，p = 0.035，OR 1.181，95% CI 1.012–1.378）与癫痫相关，FISO接近显著性（p = 0.067）。在主要简约NDI + 年龄模型（两个预测因子；每变量事件[EPV] ≈ 6.5）中，NDI仍与癫痫相关（B = ?0.146，SE 0.055，Wald χ2 = 7.04，p = 0.008，OR 0.864，95% CI 0.776–0.962），对应于NDI每降低10?3单位癫痫比值增加13.6%；bootstrap乐观校正AUC为0.712（95% CI 0.589–0.821）。在探索性四预测因子模型（NDI、ODI、FISO、年龄；EPV ≈ 3.25）中，调整后NDI（B = ?0.154，SE 0.058，Wald χ2 = 7.05，p = 0.008，OR 0.858，95% CI 0.765–0.961）和ODI（B = +0.166，SE 0.079，Wald χ2 = 4.42，p = 0.035，OR 1.181，95% CI 1.012–1.378）与癫痫相关，FISO接近显著性（p = 0.067）。预测因子间相关性中等（NDI–ODI ρ = ?0.21；NDI–FISO ρ = 0.18；ODI–FISO ρ = 0.14），所有VIF ≤1.9（容忍度≥0.53），反对存在有问题的共线性。癫痫中ODI描述性较低与其调整后正系数之间的差异符合抑制效应，即ODI仅在考虑NDI后才对模型有贡献；因此ODI应视为调整后相关物而非稳健独立预测因子。ROC曲线分析中，这些分析评估推导队列内的判别能力，不构成验证性预测。NDI的 apparent AUC为0.739（95% CI 0.617–0.861；乐观校正0.706；5折交叉验证0.701），截断值0.095时敏感度65.5%，特异度78.8%。ODI的 apparent AUC为0.671（95% CI 0.542–0.800；乐观校正0.638），截断值0.102时敏感度52.7%，特异度86.5%。两个截断值均为样本特异性。NDI AUC 0.739表明中等判别能力。由于截断值在同一队列中衍生和评估且未进行外部验证，该结果应解释为需要独立确认的判别信号，而非预测效用的证据。讨论部分指出，这项针对68例TSC儿童和39名健康对照的前瞻性NODDI研究全面表征了TSC脑微结构改变及其与癫痫表型的关联。研究发现四个关键观察：（1）TSC患者较健康对照全脑NDI、ODI和FISO广泛降低，年龄调整后效应持续；（2）TSC合并癫痫患者与对照之间存在显著差异，而无癫痫TSC患者与对照之间无显著差异；（3）早发性癫痫和婴儿痉挛性癫痫患者NDI降低显著；（4）NDI、ODI与癫痫状态之间存在独立关联。这些发现推进了对TSC病理生理学的理解，并提示NODDI衍生指标作为疾病分层定量生物标志物的潜在应用。NODDI衍生指标应视为更广泛多模态生物标志物框架的一部分而非独立预测因子。分子和免疫标志物，包括针对谷氨酸能受体靶点（如抗NMDA和抗AMPA受体抗体）的自身抗体，在部分其他原因不明的难治性癫痫患者中可识别，索引扩散MRI无法访问的免疫-炎症轴，但仅存在于少数患者且需要血清学检测。电生理标志物如高频振荡和快涟漪提供毫秒时间分辨率并直接索引致痫网络活动，但空间采样受电极放置限制。NODDI通过提供空间分辨的全脑组织水平神经突密度、方向离散度和自由水含量信息补充这些方法。因此这些模态最好视为互补：电生理定位致痫网络，血清学识别免疫介导亚组，NODDI量化相关微结构基质。未来分层策略可能结合这些数据流，需要前瞻性研究确定这种多模态整合是否优于任何单一模态。TSC中观察到的全脑NDI、ODI和FISO降低与先前较小规模NODDI研究一致，并扩展至更大的儿科队列。NDI降低表明神经突密度降低，反映轴突和树突丢失，得到TSC受累脑组织中异常神经元形态、树突分支减少和突触密度改变的组织病理学证据支持。ODI降低提示方向离散度降低，意味着纤维复杂性和组织性降低，可能反映轴突导向 disrupted 和异常白质纤维束发育。FISO降低提示自由水分数降低，与组织组成改变一致，包括细胞密度增加、胶质增生和细胞外间隙减少。这些改变的广泛性与TSC涉及弥漫性脑网络破坏而非仅局灶病理的概念一致。然而，由于分析未分割结节、结节周围组织和表现正常实质，无法确定这些改变在可见病变之外的程度；这仍是需直接检验的假设。GBSS分析显示皮质NDI和ODI降低主要位于额叶、顶叶、颞叶和枕叶区域，皮质优势与TSC已知组织病理学标志一致，包括主要影响灰质的皮质结节形成、异常神经元迁移和发育不良神经元。TBSS无显著白质发现可能反映GBSS对TSC灰质病理的更高敏感性，或者队列相对年轻，白质微结构改变可能随年龄和疾病进展而更明显。本研究关键发现是NODDI参数与癫痫状态之间的显著关联。TSC合并癫痫患者较对照NDI、ODI和FISO显著降低，而无癫痫患者无显著差异。该描述性梯度符合underlying TSC疾病过程是微结构改变主要驱动因素且癫痫标志更严重病理的解释，但同样符合反复发作贡献微结构改变或两者反映共同上游机制。横断面设计无法区分这些可能性：观察到的改变可能是癫痫的后果、underlying疾病的反映或两者兼有。反复发作贡献进行性微结构损伤的替代解释不能排除，并得到癫痫严重程度与NDI降低幅度之间剂量依赖关系的支持。只有纵向成像，在癫痫发病前后采集，才能确立时间顺序和因果方向。该发现具有重要临床意义。TSC合并癫痫患者可能代表脑病理更广泛的亚组，但该观察是否应影响临床管理无法从横断面数据确定。任何NODDI可指导预防性抗癫痫治疗的建议仍是假设性的，需要在纵向队列中证明癫痫发病前测量的NODDI参数可预测后续癫痫。TSC合并与不合并癫痫患者之间无显著未调整差异（表3），与年龄调整回归发现相反，最可能归因于13例无癫痫亚组的有限把握度、年龄混杂影响以及二分NODDI参数固有的信息损失。因此研究人员将回归结果视为值得确认的关联证据，而非证明直接比较未能检测到的亚组差异。由于亚组分析涉及多重比较，应用Benjamini–Hochberg FDR校正；主要早发性癫痫和IESS发现通过校正，而DRE–FISO关联应视为探索性。因此所有亚组结果作为假设生成呈现。研究人员还在癫痫限制设计中重新检查表型分析，使每个表型仅与无该表型的癫痫患者比较，消除了表型特异性差异仅反映癫痫与无癫痫更广泛对比的可能性。早发性癫痫和IESS与NDI和ODI降低的关联在限制后和FDR校正后仍显著，而SE关联未通过校正，DRE–FISO关联在限制分析中增强（q = 0.020）但在全队列分析中减弱（q = 0.084）。因此这些结果支持癫痫人群内神经突密度和方向离散度的表型相关差异，而SE和DRE发现仍属探索性。由于横断面设计，这些关联的方向无法确定。在癫痫患者中，早发性癫痫与NDI降低16.7%相关，IESS与降低19.0%相关，为本研究观察到的最大效应幅度，且两项关联在限制于癫痫队列和FDR校正后均显著。这些发现表明严重早发性癫痫表型与更大NDI降低之间的关联。由于横断面设计，无法确定早年癫痫是否破坏神经突发育、预先存在的微结构异常是否易感这些表型，或两者是否反映更广泛的underlying TSC病理；因此“发育双重打击”解释是假设而非发现，观察到的关联无法确立癫痫导致微结构改变。这种“发育双重打击”机制，即underlying遗传病理和关键发育期癫痫活动均贡献脑损伤，与早发性癫痫是TSC不良神经发育结局强预测因子的临床观察一致。第二，早发性癫痫和IESS可能是更严重underlying TSC病理的标志，注定发展为这些表型的患者从一开始就有更广泛的脑受累。第三，这些发现可能反映两种机制的组合，严重癫痫既反映又加剧脑病理。这些机制之间的区别具有重要治疗意义。如果早发性癫痫积极贡献进行性微结构损伤（数据无法确认的可能性），则早期癫痫控制可能减轻长期损伤。检验该假设需要纵向成像，在癫痫发病前后及与治疗关系中进行系列NODDI采集。药物难治性癫痫与FISO升高之间的关联，与整体TSC队列相对对照观察到的FISO降低方向相反，是提示药物难治性病例 distinct 病理生理机制的有趣发现。FISO升高表明自由水分数升高，可能反映细胞外间隙增加、细胞密度降低或组织水肿。在DRE背景下，该发现可能反映若干过程：（1）慢性反复发作后的反应性胶质增生和组织重塑，导致组织水含量改变；（2）癫痫活动继发的血脑屏障通透性增加，导致血管源性水肿；或（3）发展为DRE患者的不同underlying病理，包括更广泛白质受累或不同结节亚型。虽然该发现的机制基础需进一步研究，但可能表明药物难治性病例的 distinct 组织特征；该观察为探索性，基于24例患者亚组，需在独立队列中验证。NDI以及调整后程度较弱的ODI与癫痫相关，NDI显示中等判别能力（AUC = 0.739）。这些发现对临床效用而言是假设生成性的：分析为横断面，截断值为样本特异性，未进行外部验证。AUC 0.739表明中等判别能力；由于截断值在同一队列中衍生和评估，该数字应解释为判别信号证据而非验证性预测性能。可以设想若干潜在临床应用：风险分层（仅假设），NDI降低最终可能有助于识别需要更密切监测的患者，但横断面设计无法确立NODDI是否预测未来癫痫发病，需要癫痫发病前成像的纵向队列和外部验证；鉴于越来越多认识到早期癫痫控制可能改善TSC神经发育结局，这尤其相关。治疗监测，NODDI参数随时间变化可能提供mTOR抑制剂或其他新兴治疗反应的客观测量；与受多种因素影响的临床发作计数不同，定量成像测量提供脑组织完整性的直接评估。预后生物标志物（仅假设），NDI降低与早发性癫痫/IESS之间的关联提出预后价值的可能性，但未证明预后性能，仍有待纵向检验。手术规划，虽未在本研究中直接评估，NODDI在勾画致痫结节方面较常规MRI的优越对比度提示在药物难治性癫痫术前评估中的潜在应用。然而，在任何临床转化之前，需要：（1）在独立多中心队列中以预先指定截断值进行外部验证；（2）在癫痫发病前成像的纵向队列中证明前瞻性预测价值；（3）跨年龄的儿科规范参考数据；（4）评估多壳扩散MRI的成本效益和可及性。本研究以若干重要方式扩展了先前TSC的DTI和NODDI研究。DTI文献一致显示TSC中FA降低和MD增加，尤其在癫痫和认知障碍相关白质纤维束中；然而DTI对微结构室的有限特异性限制了机制解释。NODDI框架通过分离神经突密度（NDI）和方向离散度（ODI）提供更具生物学特异性的见解。本研究发现早发性癫痫中NDI降低较ODI降低更显著，提示神经突丢失而非组织紊乱是严重表型微结构改变的主要驱动因素。先前TSC的NODDI研究受较小样本量限制。Taoka等在21例TSC患者中证明白质NDI降低，Anaby等报告19例患者白质NDI与临床严重程度的相关性。本研究通过纳入68例患者并对多个癫痫表型进行系统亚组分析大幅扩展了这些发现。较大样本量使TSC合并与不合并癫痫患者以及不同癫痫表型之间的稳健比较成为可能。本研究存在若干重要局限性。第一，横断面设计无法进行因果推断。虽然观察到NODDI参数与癫痫表型之间的关联，但无法确定微结构改变是先于并促进癫痫发展，还是反复发作导致进行性微结构损伤；观察到的改变可能同样由癫痫所致、反映underlying TSC疾病过程或两者兼有。需要在癫痫发病前后进行系列成像的纵向研究以确立时间关系和因果方向。第二，组间显著年龄差异是统计分析无法完全解决的重要局限性。对照组显著年长于TSC组（平均13.09 vs. 8.1岁；中位11 vs. 8.5岁；p < 0.001），尽管调整了年龄、验证ANCOVA假设、检验年龄×组交互、使用限制性立方样条非线性建模年龄、在15–18岁重叠范围及排除>18岁对照后进行敏感性分析，并在各年龄带内确认TSC患者NDI/ODI较低，但年龄相关脑成熟导致的残余混杂无法排除。NODDI参数在整个儿童期和青春期持续变化，组间差异可能部分归因于发育阶段而非仅TSC病理。需要年龄和性别匹配招募的前瞻性研究及儿科NODDI规范轨迹来解决这一问题。第三，单时间点成像和临床评估限制了评估微结构完整性随时间动态变化的能力。鉴于TSC中癫痫的动态性质，发作类型和频率在儿童期变化，这尤其重要。第四，全脑NODDI指标源自BET实质掩膜，因此包含部分体积和残余脑脊液体素；尽管排除高FISO体素的敏感性分析结果不变且GBSS/TBSS提供组织特异性对比，全脑均值不应解释为纯组织。此外，未进行结节、结节周围或病变掩膜分析，因此微结构改变延伸至可见皮质结节之外的推断未得到数据直接支持，有待专门病变掩膜研究。第五，亚组分析为探索性且涉及多重比较。主要表型分析限于55例癫痫患者以避免表型特异性效应与癫痫vs.无癫痫区分相混杂；原始全队列对比仅作为次要分析报告。尽管在每个NODDI参数内跨四个癫痫限制表型对比应用Benjamini–Hochberg FDR校正且主要早发性癫痫和IESS发现通过校正，但IESS（n = 13）、SE（n = 13）和DRE（n = 24）亚组样本量小限制了精度和把握度，SE关联未通过校正。这些结果应视为假设生成而非验证性。第六，缺乏组织病理学确认限制了验证NODDI参数改变生物学解释的能力。虽然NODDI参数的组织学相关性已在动物模型和一些人类研究中建立，但缺乏与TSC组织病理学的直接关联。未来纳入尸检组织或手术标本的研究可解决这一局限性。第七，缺乏神经心理学评估是评估NODDI参数与认知或行为结局关联的错失机会。鉴于TAND在TSC中的高患病率，未来纳入综合神经心理学评估的研究将提供微结构改变临床意义的更完整图景。第八，队列中TSC2突变 predominance（79.6%）可能限制发现对TSC1突变患者的普遍性。突变类型与微结构改变之间的相关性值得在更大、遗传更多样的队列中研究。GBSS灰质骨架使用NODDI衍生自由水（FISO）和FA基组织分数而非独立结构MRI分割构建，可能引入测试NDI和ODI所用骨架的偏倚；基于T1的验证分析产生基本一致的结果，但未来工作最好采用完全结构锚定的实现。TBSS推断限于FA衍生白质骨架，因此低或非典型各向异性的白质区域未被采样。无癫痫参与者数量（n = 13）限制多变量模型的稳定性；因此四预测因子回归为探索性（EPV ≈ 3.25），只有简约NDI/年龄模型（EPV ≈ 6.5）应视为把握度充分。尽管VIF和相关分析未表明有问题的共线性，癫痫中ODI描述性较低时调整后ODI正系数符合抑制效应，ODI的独立贡献需在更大队列中确认。ROC分析使用同一队列进行截断值衍生和评估；尽管进行了bootstrap和交叉验证内部验证，无法进行外部验证。因此报告的AUC 0.739反映判别能力而非验证性预测，幅度中等，不应解释为确立预测或预后效用。横断面设计排除NODDI参数预测未来癫痫发病的任何推断；所有关于风险分层的陈述均为假设生成。最后，应承认NODDI建模的技术局限性。NODDI模型假设细胞内室固定扩散率值，这在膜特性改变的病理组织中可能不成立。此外，模型性能可能受b值选择和扩散方向数量影响，需要跨不同场强和采集协议的验证。未来研究方向：第一，在独立多中心队列中进行外部验证以及具有系列NODDI和预先指定终点的纵向研究是做出任何预测或预后声明的前提。第二，需要标准化采集协议的多中心研究验证发现并建立儿科NODDI参数的规范参考数据。第三，与其他成像模态整合，包括静息态功能MRI、体积分析和MR波谱，可能提供对TSC病理生理学更全面的理解。第四，与遗传和分子数据关联可能识别连接基因突变与微结构改变的特定通路。第五，整合NODDI参数与临床和遗传数据的机器学习方法可能改善个体患者结局预测。第六，更大队列中NODDI与DTI的比较研究可能阐明每种技术在特定临床应用中的相对优势。第七，将NODDI应用于mTOR抑制剂临床试验中的治疗监测可确立其作为治疗反应定量生物标志物的效用。结论部分指出，这项针对68例TSC儿童和39名健康对照的综合NODDI研究表明，先进扩散MRI可灵敏检测TSC中广泛的脑微结构改变，特征为全脑NDI、ODI和FISO显著降低。这些发现表明弥漫性神经突丢失、纤维复杂性降低和组织组成改变，提示弥漫性而非局灶性微结构受累。TSC合并癫痫患者较对照NDI、ODI和FISO降低，而无癫痫患者无此表现；然而癫痫亚组之间的直接比较不显著，且无癫痫组样本小。因此数据符合但未确立癫痫相关因素对微结构改变的贡献，横断面设计无法确定观察到的改变是由癫痫所致还是反映underlying疾病。早发性癫痫和婴儿痉挛中显著的NDI降低与这些严重表型相关，但因果方向仍未确定。药物难治性癫痫独特地与FISO升高相关，提示该亚组 distinct 病理生理机制。NDI以及调整后程度较弱的ODI与癫痫相关，NDI显示中等判别性能（AUC = 0.739；95% CI 0.617–0.861）。这些指标可作为癫痫相关微结构负荷的标志物；其是否能识别未来风险患者需要纵向、外部验证研究。然而，横断面设计和组间显著年龄差异是未来研究必须解决的重要局限性。需要年龄匹配对照的纵向研究以确立预测价值并确定NODDI参数是否能指导早期干预策略。总之，NODDI衍生参数为表征TSC脑微结构改变及其与癫痫表型的关联提供了有价值的定量生物标志物。这些发现推进了对TSC病理生理学的理解，支持先进扩散MRI在临床研究中的应用，并为未来生物标志物开发和个性化治疗策略奠定基础。NODDI与其他成像模态、遗传数据和临床评估的整合有望改善这一复杂神经发育障碍的结局。