基于脑白质完整性与DNA甲基化的遗传影像多模态融合分析方法技术

本发明专利技术公开了一种基于脑白质完整性与DNA甲基化的遗传影像多模态融合分析方法。本发明专利技术通过分析受试者全脑白质纤维完整性FA值，将患有某神经精神疾病患者与正常人进行组间比较，在Allen脑库中键入差异区域名称，获得在这些脑区高表达的基因名称，再从基因组甲基化样本中提取这些基因对应CpG位点的甲基化水平，采用相关分析统计这些CpG位点甲基化水平与FA值之间的关联。本发明专利技术对神经精神类疾病脑结构发育模式及其背后的表观遗传基础进行多模态融合分析，这种分析方法巧妙地将FA与DNA甲基化两模态数据结合起来分析，对现有的影像加表观遗传多模态数据信息的深入挖掘带来重大的应用及指导价值。用及指导价值。用及指导价值。

【技术实现步骤摘要】
基于脑白质完整性与DNA甲基化的遗传影像多模态融合分析方法


[0001]本专利技术同时涉及人脑结构连接与生物
，公开了一种基于脑白质完整性与DNA甲基化的表观遗传/影像多模态融合分析方法，将影像与表观遗传相融合的多模态分析方法用于疾病的早发现早诊断早治疗中，对揭示疾病进展的临床医学研究过程具有重要的意义。

技术介绍

[0002]近年来快速发展的磁共振弥散张量成像(Diffusion Tensor Imaging,DTI)技术，再结合对纤维走形和排列敏感性较高的部分各向异性(fractional anisotropy,FA)值，可对全脑白质纤维连接微结构进行定量刻画。近年来借助DTI和FA，人们可以从全脑水平上对疾病神经病理学方面的发病机制进行一系列详尽研究。尤其是对于当前的神经精神疾病，多数学者认为其病理机制可能并非与个别脑区的结构完整性异常有关，而与多个脑区的异常有关，因此从全脑水平上结合DTI和FA就能很好地阐明这类疾病的潜在神经病理机制。
[0003]上述脑白质纤维连接是从宏观角度刻画疾病的神经病理学机制，然而，微观上，脑白质纤维连接可能受基因表达及神经可塑性变化的影响，而基因表达又在很大程度上受表观遗传的调控。表观遗传是指DNA序列无改变，但基因表达发生变化，进而改变基因功能并引起表型变化。在表观遗传修饰中，DNA甲基化是研究得较深入且较重要的调节基因表达的方式之一。近年来尸检研究在神经精神疾病患者脑内发现一系列DNA甲基化水平改变。然而，脑组织既昂贵又较难获取，且不能用于纵向研究，为突破这一局限，一些研究试图探寻患者外周血细胞内的DNA甲基化标记，这将有助于对该病的早期检测，而且可用于预测患者疗效反应，以期指导个体化精准治疗。此外，有研究表明大脑和外周血样本两者之间特定基因组区域的CpG位点甲基化水平显著相关，这就提示：脑组织中一些DNA甲基化标记可通过外周血样本的相应位点反映出来；外周血可以作为脑组织的替代物，满足大规模或纵向研究的需要
[10]。
[0004]现有技术多局限于单独影像或遗传模态，不利于多模态数据的深入挖掘。因此，将脑白质纤维连接完整性FA值与外周血DNA甲基化水平相融合的表观遗传/影像多模态分析方法，可巧妙地将脑白质完整性与DNA甲基化两模态数据结合起来进行分析，将对现有的影像加表观遗传多模态数据信息的深入挖掘、学习带来重大的应用及指导价值。

技术实现思路

[0005]本专利技术的目的在于克服现有技术的不足之处而提供一种基于脑白质完整性与DNA甲基化的表观遗传/影像多模态融合分析方法，用于人类疾病尤其是神经精神疾病的病因及进展的多模态全面研究。
[0006]为实现上述目的，本专利技术提供一种基于脑白质完整性与DNA甲基化的表观遗传/影像多模态融合分析方法，包括以下步骤：
[0007]步骤S1：获取多个被试对象(包括某神经精神疾病患者和健康对照)的DTI数据，并进行预处理，剔除不符合规定的被试数据；
[0008]步骤S2：计算DTI数据中全脑各体素的FA值；
[0009]步骤S3：一般线性回归分析统计某神经精神疾病患者与健康对照FA值的组间差异；
[0010]步骤S4：在Allen脑库中键入步骤3所得到这些差异区域的脑区名称，获得在这些脑区高表达的基因名称；
[0011]步骤S5：获取受试者(患者)的外周血全基因组DNA甲基化数据。从患者外周血全基因组甲基化样本中提取步骤4中所获得的基因对应的CpG位点的甲基化水平；
[0012]步骤S6：采用Pearson或Spearman相关分析统计步骤5所得到的CpG位点甲基化水平与步骤3中得到的有组间差异的脑白质连接完整性FA值之间的关联；实现疾病相关的脑结构发育模式及其背后的表观遗传基础的多模态融合分析。
[0013]作为优选的方案，所述步骤S1中，预处理包括以下步骤：
[0014]采用FSL 5.0工具包(http://www.fmrib.ox.ac.uk/fsl)对DTI影像数据进行头动和涡流扭曲矫正；用弥散工具包软件Diffusion Toolkit software(TrackVis.org)通过线性最小二乘法拟合方法来估量每个体素的弥散张量模型，得到各自的b0图和FA图。
[0015]进一步地，所述步骤S2中，FA计算包括以下步骤：
[0016]应用SPM8工具包(Statistical parametric Mapping，http://www.fil.ion.ucl.ac.uk/spm)对FA图进行处理。通过在本地弥散空间的线性变化将每一个体的T1加权像配准到b0图像；并将变换参数应用于对应的FA，得到配准后的FA图。再利用仿射变换和一系列非线性弯曲将配准的结构影像绘制到蒙特利尔神经研究所(Montreal Neurological Institute，MNI)T1模板；将这一步得到的空间变形参数应用于各样本配准后的FA图，将其转换到标准MNI空间(3
×3×
3mm)。使用半高全宽(full width at half maximum,FWHM)为8mm的高斯核对标准化后的FA图进行平滑。
[0017]更进一步地，所述步骤S3中，某神经精神疾病患者与健康对照FA值的组间差异的一般线性回归分析包括以下步骤：
[0018]组间比较通过SPM8进行一般线性回归分析，将年龄，性别，教育年限和头动参数作为协变量进行回归分析，采用高斯随机场方法(Gaussian Random Field,GRF)进行多重比较校正(体素水平P<0.01，簇水平P<0.05)。
[0019]更进一步地，所述步骤S4中Allen脑库的使用包括以下步骤：
[0020]进入网站http://human.brain
‑
map.org/；点击左侧Differential Search；目标组织(Target Structures)选择上述步骤S3中所得到的FA值有组间差异的脑区的名称；对比结构(Contrast Structures)选择全脑(Br)；捐献者(Donors)选择全部捐献者(All Donors)；点击Search；点击左下角下载(Download this data)；筛选Fold Change值＞2的基因。
[0021]更进一步地，所述步骤S5中获取受试者外周血全基因组DNA甲基化数据，包括以下步骤：
[0022]采用Illumina450K甲基化芯片检测受试者(患者)的全基因组的CpG位点的甲基化水平，去掉所检测到的CpG位点中存多态性的位点，并获得剩余CpG位点的甲基化Beta值。从
患者剩余CpG位点中提取步骤S4中所获得的基因对应的CpG位点的甲基化水平。
[0023]更进一步地，所述步骤S6中采用Pearson或Spearman相关分析统计步骤S5所得到的CpG位点甲基化水平与步骤S3中得到的有组间差异的脑白质连接完整性FA值之间的本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于脑白质完整性与DNA甲基化的表观遗传/影像多模态融合分析方法，其特征在于：包括以下步骤：S1：获取多个被试对象的DTI数据，并进行预处理，剔除不符合规定的被试数据；所述被试对象包括包括受试者和健康对照；S2：计算DTI数据中全脑各体素的FA值；S3：一般线性回归分析统计受试者与健康对照FA值的组间差异；S4：在Allen脑库中键入步骤S3所得到这些差异区域的脑区名称，获得在这些脑区高表达的基因名称；S5：获取受试者的外周血全基因组DNA甲基化数据；从受试者外周血全基因组甲基化样本中提取步骤S4中所获得的基因对应的CpG位点的甲基化水平；S6：采用Pearson或Spearman相关分析统计步骤S5所得到的CpG位点甲基化水平与步骤S3中得到的有组间差异的脑白质连接完整性FA值之间的关联；实现疾病相关的脑结构发育模式及其背后的表观遗传基础的多模态融合分析。2.根据权利要求1所述的基于脑白质完整性与DNA甲基化的表观遗传/影像多模态融合分析方法，其特征在于：所述步骤S1中，预处理包括以下步骤：采用FSL 5.0工具包对DTI影像数据进行头动和涡流扭曲矫正；用弥散工具包软件Diffusion Toolkit software通过线性最小二乘法拟合方法来估量每个体素的弥散张量模型，得到各自的b0图和FA图。3.根据权利要求1或2所述的基于脑白质完整性与DNA甲基化的表观遗传/影像多模态融合分析方法，其特征在于：所述步骤S2中，FA计算包括以下步骤：应用SPM8工具包对FA图进行处理；通过在本地弥散空间的线性变化将每一个体的T1加权像配准到b0图像；并将变换参数应用于对应的FA，得到配准后的FA图；再利用仿射变换和一系列非线性弯曲将配准的结构影像绘制到蒙特利尔神经研究所T1模板；将这一步得到的空间变形参数应用于各样本配准后的FA图，将其转换到标准MNI空间；使用半高全宽为8mm的高斯核对标准化后的FA图进行平滑。4.根据权利要求1或2所述的基于脑白质完整性与DNA甲基化的表观遗传/影像多模态融合分析方法，其特征在于：所述步骤S3中，受试者与健康对照FA值的组间差异的一般线性回归分析包括以下步骤：组间比较通过SPM8进行一般线性回归分析，将年龄，性别，教育年限和头动参数作为协变量进行回归分析，采用高斯随机场方法进行多重比较校正，体素水平P<0.01，簇水平P<0.05。5.根据权利要求3所述的基于脑白质完整性与DNA甲基化的表观遗传/影像多模态融合分析方法，其特征在于：所述步骤S3中，受试者与...

【专利技术属性】
技术研发人员：宗小芬，胡茂林，刘忠纯，徐顺生，翁深宏，
申请(专利权)人：武汉大学，
类型：发明
国别省市：