CN111666895B - 基于深度学习的神经干细胞分化方向预测系统及方法 - Google Patents
基于深度学习的神经干细胞分化方向预测系统及方法 Download PDFInfo
- Publication number
- CN111666895B CN111666895B CN202010514541.1A CN202010514541A CN111666895B CN 111666895 B CN111666895 B CN 111666895B CN 202010514541 A CN202010514541 A CN 202010514541A CN 111666895 B CN111666895 B CN 111666895B
- Authority
- CN
- China
- Prior art keywords
- module
- neural stem
- relu
- model
- dsconv
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 210000001178 neural stem cell Anatomy 0.000 title claims abstract description 126
- 230000024245 cell differentiation Effects 0.000 title claims abstract description 71
- 238000000034 method Methods 0.000 title claims abstract description 43
- 238000013135 deep learning Methods 0.000 title claims abstract description 28
- 210000004027 cell Anatomy 0.000 claims abstract description 100
- 230000004069 differentiation Effects 0.000 claims abstract description 79
- 238000012549 training Methods 0.000 claims abstract description 35
- 238000013527 convolutional neural network Methods 0.000 claims abstract description 19
- 230000014509 gene expression Effects 0.000 claims description 57
- 238000007781 pre-processing Methods 0.000 claims description 30
- 238000013528 artificial neural network Methods 0.000 claims description 23
- 210000002569 neuron Anatomy 0.000 claims description 17
- 210000001130 astrocyte Anatomy 0.000 claims description 16
- 210000004248 oligodendroglia Anatomy 0.000 claims description 16
- 210000000130 stem cell Anatomy 0.000 claims description 14
- 238000004422 calculation algorithm Methods 0.000 claims description 10
- 238000010606 normalization Methods 0.000 claims description 10
- 238000011176 pooling Methods 0.000 claims description 9
- 230000006870 function Effects 0.000 claims description 8
- 238000012258 culturing Methods 0.000 claims description 6
- 238000009472 formulation Methods 0.000 claims description 5
- 239000000203 mixture Substances 0.000 claims description 5
- 210000005036 nerve Anatomy 0.000 claims description 5
- 230000004913 activation Effects 0.000 claims 4
- 238000005516 engineering process Methods 0.000 abstract description 5
- 238000000684 flow cytometry Methods 0.000 abstract description 3
- 238000012545 processing Methods 0.000 abstract description 3
- 238000011005 laboratory method Methods 0.000 abstract 1
- 238000005457 optimization Methods 0.000 abstract 1
- 238000012360 testing method Methods 0.000 description 23
- 108090000742 Neurotrophin 3 Proteins 0.000 description 7
- 239000002609 medium Substances 0.000 description 6
- 238000011160 research Methods 0.000 description 6
- 238000002474 experimental method Methods 0.000 description 5
- 230000001939 inductive effect Effects 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 238000003125 immunofluorescent labeling Methods 0.000 description 4
- 239000000411 inducer Substances 0.000 description 4
- 108020005198 Long Noncoding RNA Proteins 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 230000007547 defect Effects 0.000 description 3
- 239000007850 fluorescent dye Substances 0.000 description 3
- 210000003958 hematopoietic stem cell Anatomy 0.000 description 3
- 239000000463 material Substances 0.000 description 3
- 150000003384 small molecules Chemical class 0.000 description 3
- 238000010186 staining Methods 0.000 description 3
- 238000007792 addition Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000018109 developmental process Effects 0.000 description 2
- 238000002372 labelling Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 239000003147 molecular marker Substances 0.000 description 2
- 230000001537 neural effect Effects 0.000 description 2
- 238000003062 neural network model Methods 0.000 description 2
- ORILYTVJVMAKLC-UHFFFAOYSA-N Adamantane Natural products C1C(C2)CC3CC1CC2C3 ORILYTVJVMAKLC-UHFFFAOYSA-N 0.000 description 1
- 108090000715 Brain-derived neurotrophic factor Proteins 0.000 description 1
- 102000004219 Brain-derived neurotrophic factor Human genes 0.000 description 1
- 206010008342 Cervix carcinoma Diseases 0.000 description 1
- 108010005939 Ciliary Neurotrophic Factor Proteins 0.000 description 1
- 102100031614 Ciliary neurotrophic factor Human genes 0.000 description 1
- 239000006144 Dulbecco’s modified Eagle's medium Substances 0.000 description 1
- 108010025020 Nerve Growth Factor Proteins 0.000 description 1
- 102000015336 Nerve Growth Factor Human genes 0.000 description 1
- 208000006105 Uterine Cervical Neoplasms Diseases 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 210000004504 adult stem cell Anatomy 0.000 description 1
- 238000013473 artificial intelligence Methods 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 210000000988 bone and bone Anatomy 0.000 description 1
- 230000009087 cell motility Effects 0.000 description 1
- 239000006285 cell suspension Substances 0.000 description 1
- 210000003169 central nervous system Anatomy 0.000 description 1
- 208000015114 central nervous system disease Diseases 0.000 description 1
- 201000010881 cervical cancer Diseases 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 239000004205 dimethyl polysiloxane Substances 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 238000000799 fluorescence microscopy Methods 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 239000001963 growth medium Substances 0.000 description 1
- 229940088597 hormone Drugs 0.000 description 1
- 239000005556 hormone Substances 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000011534 incubation Methods 0.000 description 1
- 230000006698 induction Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 230000003902 lesion Effects 0.000 description 1
- 230000033001 locomotion Effects 0.000 description 1
- 238000010801 machine learning Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000011089 mechanical engineering Methods 0.000 description 1
- 230000009456 molecular mechanism Effects 0.000 description 1
- 230000004660 morphological change Effects 0.000 description 1
- 210000000944 nerve tissue Anatomy 0.000 description 1
- 230000007170 pathology Effects 0.000 description 1
- 229920000435 poly(dimethylsiloxane) Polymers 0.000 description 1
- -1 polydimethylsiloxane Polymers 0.000 description 1
- 229920000642 polymer Polymers 0.000 description 1
- 239000011148 porous material Substances 0.000 description 1
- 238000004393 prognosis Methods 0.000 description 1
- 238000007637 random forest analysis Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 239000000758 substrate Substances 0.000 description 1
- 238000012706 support-vector machine Methods 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
- 238000002560 therapeutic procedure Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V20/00—Scenes; Scene-specific elements
- G06V20/60—Type of objects
- G06V20/69—Microscopic objects, e.g. biological cells or cellular parts
- G06V20/698—Matching; Classification
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12M—APPARATUS FOR ENZYMOLOGY OR MICROBIOLOGY; APPARATUS FOR CULTURING MICROORGANISMS FOR PRODUCING BIOMASS, FOR GROWING CELLS OR FOR OBTAINING FERMENTATION OR METABOLIC PRODUCTS, i.e. BIOREACTORS OR FERMENTERS
- C12M41/00—Means for regulation, monitoring, measurement or control, e.g. flow regulation
- C12M41/30—Means for regulation, monitoring, measurement or control, e.g. flow regulation of concentration
- C12M41/36—Means for regulation, monitoring, measurement or control, e.g. flow regulation of concentration of biomass, e.g. colony counters or by turbidity measurements
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/21—Design or setup of recognition systems or techniques; Extraction of features in feature space; Blind source separation
- G06F18/214—Generating training patterns; Bootstrap methods, e.g. bagging or boosting
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/04—Architecture, e.g. interconnection topology
- G06N3/045—Combinations of networks
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Biomedical Technology (AREA)
- Chemical & Material Sciences (AREA)
- Data Mining & Analysis (AREA)
- General Health & Medical Sciences (AREA)
- Zoology (AREA)
- Wood Science & Technology (AREA)
- Evolutionary Computation (AREA)
- Organic Chemistry (AREA)
- Molecular Biology (AREA)
- Artificial Intelligence (AREA)
- Computational Linguistics (AREA)
- Computing Systems (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Multimedia (AREA)
- Biophysics (AREA)
- Analytical Chemistry (AREA)
- Biotechnology (AREA)
- Microbiology (AREA)
- Sustainable Development (AREA)
- Evolutionary Biology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Biochemistry (AREA)
- Genetics & Genomics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Image Analysis (AREA)
- Investigating Or Analysing Biological Materials (AREA)
Abstract
本发明涉及基于深度学习的神经干细胞分化方向预测系统及方法。本发明结合流式细胞技术对于细胞的高通量处理能力,创造性地利用实验手段,收集全景流式细胞仪采集的神经干细胞分化培养的细胞图像,建立起神经干细胞三系分化的细胞训练数据集,再利用卷积神经网络进行模型训练优化,建立了一套精准高效、操作简易、省时低消耗的预测体系,能够克服现有实验室方法评判神经干细胞分化的各种技术的限制,方便快捷、高效精准的预测神经干细胞的分化方向。
Description
技术领域
本发明涉及生物医学和人工智能技术领域,具体地说,涉及基于深度学习的神经干细胞分化方向预测系统及方法。
背景技术
神经干细胞(Neural stem cells,NSCs)是一类存在于中枢神经系统中的干细胞,具有以下特征:1)能形成神经组织;2)具有自我繁殖和自我更新能力; 3)具有分化为神经元、星形胶质细胞和少突胶质细胞的潜能。近年来,神经干细胞在中枢神经系统疾病的基础生物学探究和基于细胞的治疗方面都具有巨大的潜力。目前,神经干细胞的治疗前景主要受限于不能精准控制培养过程中的干细胞行为,如何精准调控神经干细胞的定向分化成为一个至关重要的研究热点。而在常规的实验室基础研究中,判断干预后神经干细胞的分化方向往往需要复杂的实验验证过程,要了解干细胞分化的结果需要花费1个月甚至更长的时间,而在诱导NSCs因子的筛选应用中,需要预先知道在测试因子作用下的NSCs将会分化成哪种类型的细胞。
常规实验手段包括免疫荧光染色,流式细胞染色等手段,效率低,耗时长,耗费高。创建一种集快速、精准、高效、简易为一体的预测判断神经干细胞分化方向的方法,可以有效提高相关科学研究的进度,为NSCs的治疗应用前景创造有力的支持。有一些相关研究也一直在致力于缩短判断干细胞分化方向的时间,据Nature Methods杂志(doi:10.1038/nmeth.1487)报道,密歇根大学生物工程系与机械工程系和我国台湾成功大学医学院骨关节研究中心等处的研究人员开发了一种新型培养基,将成体干细胞在其中培养24小时,就能预测其分化方向,这种新型培养基是一种主要由聚二甲基硅氧烷弹性聚合物制成的硬度可调整的基质支架。干细胞在分化过程中的牵引力会发生改变,因此研究者认为,支架牵引力的变化能作为细胞分化的早期预测因素,可以根据底物的柔软度观察细胞分化,并使用荧光显微技术检测支架的弯曲度以确定牵引力的大小。此方法有效缩短了判断干细胞分化的时间,但受制于实验设备、人工操作,其精准性仍然具有很大的局限性,且荧光显微镜技术、染色技术,包括细胞基质的合成都存在着耗时耗费的缺陷。
近年来,伴随着计算机技术的迅速发展,其在生物医学研究中也开始受到广泛关注。例如专利文献CN107577924A,公开日2018.01.12,公开了一种基于深度学习的lncRNA亚细胞位置预测算法,包括以下步骤:S1、组建lncRNA 序列亚细胞位置预测训练基准数据集;S2、从lncRNA序列中提取原始的k-mer 频率特征;S3、采用上采样算法对训练数据进行预处理;S4、采用三层堆叠的栈式自编码器对提取的频率信息进行特征提取,将得到的频率特征和利用模型抽取的特征分别接入随机森林分类器和支持向量机模型训练;S5、将四个分类器的输出概率经过基于神经网络的整合算法进行融合训练输出预测结果。又如专利文献CN110334565A,公开日2019.10.15,公开了一种显微镜病理照片的宫颈癌病变细胞分类系统,包括:检测网络模型和分类网络模型,所述检测网络模型利用不完全分类的细胞标注数据进行训练,训练完成后用于对显微镜照片中的细胞进行检测;所述分类网络模型根据不完全的细胞分类标注数据进行训练,训练完成后用于对检测到的细胞进行分类,并对细胞的具体类别进行预测。
近年来已经有一些将机器学习应用于干细胞分化的研究,据Nature Methods杂志(doi:10.1038/nmeth.4182)报道分化会改变干细胞和祖细胞的分子特性,导致形状和运动特性的改变,由此利用光斑显微镜和细胞运动的图像块,建立了一个深层神经网络预测在分化主要造血祖细胞时的谱系选择。但此研究主要集中在造血干细胞的谱系选择,基于造血干细胞的变化图像进行预测,无法普遍性地应用到神经干细胞上,且此项研究基于显微成像技术,在数据量方面有一定的缺陷,无法达到高通量数据训练达到的更高精度和效率的判断体系。
综合上述,现有的神经干细胞分化方向的判断和预测系统不能够满足日益增长的科研需求,因此,建立一种更加简便、高效、精准、低成本的高通量预测体系十分关键。
发明内容
本发明的目的是针对现有技术中的不足,提供一种基于深度学习的神经干细胞分化方向预测方法和系统。
第一方面,本发明提供了一种基于深度学习的神经干细胞分化方向预测方法,包括以下步骤:
S100,利用全景流式细胞仪收集神经干细胞分化培养的细胞图像,包括分化为神经元、星形胶质细胞和少突胶质细胞三个方向的细胞图像;
S200,将全景流式细胞仪收集的神经干细胞分化培养的细胞图像,输入卷积神经网络模型,由卷积神经网络模型自动提取有效的特征,训练得到神经干细胞分化方向预测模型;
S300,将待分析分化方向的神经干细胞,利用全景流式细胞仪收集其细胞图像并输入所述神经干细胞分化方向预测模型,得到预测结果。
作为一个优选例,所述卷积神经网络模型表达式为:
模型子网络模块1(m1)的表达式为:
z1,1=Conv1,1(x);
a1,1=Relu(z1,1);
z1,2=Conv1,2(a1,1);
a1,2=Relu(z1,2);
u=a1,2;
模型子网络模块2(m2)的表达式为:
z2,1;s=Conv2,1;sx;
z2,1=DSConv2,1(x);
a2,1=Relu(z2,1);
z2,2=DSConv2,2(a2,1);
p2,1=Pool(z2,2);
u=p2,1+z2,1;s;
模型子网络模块3(m3)的表达式为:
z3,1;s=Conv3,1;sx;
z3,1=DSConv3,1(x);
a3,1=Relu(z3,1);
z3,2=DSConv3,2(a3,1);
p3,1=Pool(z3,2);
u=p3,1+z3,1;s;
模型子网络模块4(m4)的表达式为:
a4,1=Relu(x);
z4,1=DSConv4,1(a4,1);
a4,2=Relu(z4,1);
z4,2=DSConv4,2(a4,2);
a4,3=Relu(z4,2);
z4,3=DSConv4,3(a4,3);
u=z4,3+x;
模型子网络模块5(m5)的表达式为:
z5,1;s=Conv5,1;sx;
a5,1=Relu(x);
z5,1=DSConv5,1(a5,1);
a5,2=Relu(z5,1);
z5,2=DSConv5,2(a5,1);
p5,1=Pool(z5,2);
u=p5,1+z5,1;s;
模型子网络模块6(m6)的表达式为:
z6,1=Conv6,1(x);
a6,1=Relu(z6,1);
z6,2=Conv6,2(a6,1);
a6,2=Relu(z6,2);
u=a62;
模型整体表达式为:
u1=m1(x);
u2=m2(u1);
u3=m3(u2);
u4,1=m4,1(u3);
u4,2=m4,2(u4,1);
u4,3=m4,3(u4,2);
u4=m4,3(u4,3);
u5=m5(u4);
u6=m6(u5);
p=GPool(u6);
zw=Fc(p);
以上符号的含义为:
作为另一优选例,步骤S100中,神经干细胞分化培养的培养基配方和采集细胞图像的时间点如下:
作为另一优选例,步骤S100中,训练集中的分化为星形胶质方向、神经元细胞方向、少突胶质细胞方向和NT3处理组的细胞图像的比例为 (1.9~2.1):(1.25~1.45):1:(1~1.1)。
作为另一优选例,所述神经干细胞分化方向预测模型包括图片预处理模块和神经网络模块;所述图片预处理模块用于接收全景流式细胞仪收集的神经干细胞分化培养的细胞图像,预处理得到标准化的图片数据;所述神经网络模块用于接收标准化的图片数据,得到最终的神经干细胞分化方向的预测结果。
更优选地,所述图片预处理模块包括通道合并模块和图片标准化模块;所述通道合并模块用于输入具有相同的高度H和宽度W的单颜色通道细胞图像,将单通道图片沿通道合并为多通道“图片”表示,合并后的图片张量表示为[C, H,W],其中C表示输入的颜色通道数量;所述图片标准化模块用于承接通道合并模块,输入数据为合并后的多通道图片张量,符号表示为[C,H,W],将输入数据标准化为[C,45,30]的张量表示。
更优选地,所述图片标准化模块使用双立方插值算法将[C,H,W]的图像张量转化为[C,45,30],将经过插值操作的图像张量作正则化操作。
作为另一优选例,所述神经干细胞分化方向的预测结果还包括神经干细胞的分化比例。
第二方面,本发明提供了一种基于深度学习的神经干细胞分化方向预测系统,包括:
全景流式细胞仪,用于收集神经干细胞分化培养的细胞图像;
神经干细胞分化方向预测模型,用于接收全景流式细胞仪收集的待预测分化方向的神经干细胞分化培养的细胞图像,得到预测结果;所述神经干细胞分化方向预测模型是按照以下方法构建的:
S100,利用全景流式细胞仪收集神经干细胞分化培养的细胞图像,包括分化为神经元、星形胶质细胞和少突胶质细胞三个方向的细胞图像;
S200,将全景流式细胞仪收集的神经干细胞分化培养的细胞图像,输入卷积神经网络模型,由卷积神经网络模型自动提取有效的特征,训练得到神经干细胞分化方向预测模型;
所述卷积神经网络模型表达式为:
模型子网络模块1(m1)的表达式为:
z1,1=Conv1,1(x);
a1,1=Relu(z1,1);
z1,2=Conv1,2(a1,1);
a1,2=Relu(z1,2);
u=a1,2;
模型子网络模块2(m2)的表达式为:
z2,1;s=Conv2,1;sx;
z2,1=DSConv2,1(x);
a2,1=Relu(z2,1);
z2,2=DSConv2,2(a2,1);
p2,1=Pool(z2,2);
u=p2,1+z2,1;s;
模型子网络模块3(m3)的表达式为:
z3,1;s=Conv3,1;sx;
z3,1=DSConv3,1(x);
a3,1=Relu(z3,1);
z3,2=DSConv3,2(a3,1);
p3,1=Pool(z3,2);
u=p3,1+z3,1;s;
模型子网络模块4(m4)的表达式为:
a4,1=Relu(x);
z4,1=DSConv4,1(a4,1);
a4,2=Relu(z4,1);
z4,2=DSConv4,2(a4,2);
a4,3=Relu(z4,2);
z4,3=DSConv4,3(a4,3);
u=z4,3+x;
模型子网络模块5(m5)的表达式为:
z5,1;s=Conv5,1;sx;
a5,1=Relu(x);
z5,1=DSConv5,1(a5,1);
a5,2=Relu(z5,1);
z5,2=DSConv5,2(a5,1);
p5,1=Pool(z5,2);
u=p5,1+z5,1;s;
模型子网络模块6(m6)的表达式为:
z6,1=Conv6,1(x);
a6,1=Relu(z6,1);
z6,2=Conv6,2(a6,1);
a6,2=Relu(z6,2);
u=a62;
模型整体表达式为:
u1=m1(x);
u2=m2(u1);
u3=m3(u2);
u4,1=m4,1(u3);
u4,2=m4,2(u4,1);
u4,3=m4,3(u4,2);
u4=m4,3(u4,3);
u5=m5(u4);
u6=m6(u5);
p=GPool(u6);
zw=Fc(p);
以上符号的含义为:
作为一个优选例,所述神经干细胞分化方向预测模型包括图片预处理模块和神经网络模块;所述图片预处理模块用于接收全景流式细胞仪收集的神经干细胞分化培养的细胞图像,预处理得到标准化的图片数据;所述图片预处理模块包括通道合并模块和图片标准化模块;所述通道合并模块用于输入具有相同的高度H和宽度W的单颜色通道细胞图像,将单通道图片沿通道合并为多通道“图片”表示,合并后的图片张量表示为[C,H,W],其中C表示输入的颜色通道数量;所述图片标准化模块用于承接通道合并模块,输入数据为合并后的多通道图片张量,符号表示为[C,H,W],将输入数据标准化为[C,45,30]的张量表示;所述神经网络模块用于接收标准化的图片数据,得到最终的神经干细胞分化方向的预测结果。
本发明优点在于:
1、本发明将深度学习方法应用到神经干细胞分化方向的预测上,建立了基于深度学习的神经干细胞(NSCs)分化方向预测系统及方法。现有的NSCs 分化方向的判断手段都是基于常规实验手段,受到多种设备条件和人为操作方面的干扰,同时其鉴定的速度和准确度也受制于分子标记在分化过程中的表达程度。本发明创新性地构建了从实验所得的数据集,合理设计训练数据集,构建起从NSCs有效分化为神经元、星形胶质细胞和少突胶质细胞的细胞图片数据集。结合深度学习方法提取分化早期的细胞形态学变化,建立起对于NSCs 分化的预测体系,前瞻性地判断不同方式处理后NSCs分化方向。
2、除了混合测试集外,本发明还在独立测试数据中引入不同神经因子/小分子/纳米材料体系等处理后的NSCs向神经元方向分化,这些诱导物分别作用于不同的受体和信号通路,结果表明本发明的系统可以利用无荧光染色的流式单细胞图像在测试因子/小分子/纳米材料短时间处理细胞后精准判断出该细胞的未来分化方向,其中,神经元、星形胶质细胞和少突胶质细胞方向最短分别需耗时1天、12小时和1天,预测得到分化比例分别为0.827、0.959和0.806,在5天、2天和3天分别为0.830、0.975和0.920,此预测的细胞分化比例与免疫荧光染色实验获得的因子/小分子/纳米材料处理后的神经干细胞的在相应分化方向的比例高度一致。从而证实了本发明模型的普适性、精准性和抗干扰性。
3、本发明所用数据为未经荧光染色标记的全景流式细胞仪所收集的细胞图像,设备要求基础,人为操作简单易行,构建该系统的成本和测试成本很低,且处理时间较短,不受限于已知的分子标记手段,在细胞显示出分子标记可识别的变化之前,在短至1天的处理时间点精确预测出神经干细胞的分化方向,能在保证精准预测的基础上有效缩短判断时间。从而可以大幅度简化和加快神经干细胞分化相关新型诱导方式的鉴定,在干细胞治疗神经中枢系统疾病的发展方面有很大的应用价值。
4、本发明在神经网络构建方面采用了多层深度可分卷积层,并加入残差捷径,使得模型的训练更加迅速且稳健,能够有效提高模型的泛化能力。
5、本发明在数据预处理方面使用双立方插值算法将收集得到的细胞图像统一为45x30的大小,简便快捷。
附图说明
附图1:本发明神经干细胞分化方向预测模型训练测试流程示意图。
附图2:本发明基于深度学习的神经干细胞分化方向预测方法的流程示意图。
附图3:本发明基于深度学习的神经干细胞分化方向预测系统的结构框图。
具体实施方式
为使得本发明的发明目的、特征、优点能够更加的明显和易懂,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而非全部实施例。实施例1本发明基于深度学习的神经干细胞分化方向预测模型
以下为模型构建过程,包括两个主要部分,分别是训练数据的构建和神经网络模块的建立。
一、数据的构建过程
训练数据来自于神经干细胞三系分化体系构建:
1.将第三代到第五代状态良好的神经干细胞用无因子培养基 (DMEM/F12+1*N2+1*B27+1%双抗)制备成单细胞悬液,以5万/cm2的密度接种于PLO包被过的孔板中,置于37℃,5%CO2恒温培养箱培养过夜使细胞贴壁。
2.过夜后将无因子培养液替换为分化培养基,各组分化培养基配方如下:
3.同时,构建了一个非独立测试组NT3组和一个独立测试组LDH-NT3 (层状双氢氧化物搭载NT3,Layered double hydroxide-NT3)组,方法如下:
4.在设置好的时间点收集细胞,利用全景流式细胞仪收集细胞图像,其中神经干细胞三系分化的和NT3处理组数据中分别取80%作为训练数据 (119,533个细胞图像),20%作为非独立测试数据(29,895个细胞图像), LDH-NT3处理组作为不参与模型训练的独立测试组(20,801个细胞图像)。
分化组 | 细胞图像数目 |
诱导分化为星形胶质细胞 | 55466 |
诱导分化为神经元 | 37771 |
诱导分化为少突胶质细胞 | 27687 |
NT3诱导分化组 | 28504 |
一共149428,其中119533个图像作为训练数据,29895作为非独立测试数据。
二、模型的建立过程
模型主要包括两个主要部分:
(一)图片预处理模块
利用全景流式细胞仪收集处理后的细胞图像,数据预处理方面使用双立方插值算法将收集得到的细胞图像统一为45x30的大小,作为模型的训练数据。
所述图片预处理模块包括两个子模块:
(1)通道合并模块
本模块的输入数据为细胞的单颜色通道图片,每个颜色通道来源于相应的细胞着色通道。这些单颜色通道图片必须具有相同的高度H和宽度W。通道合并模块将这些单通道图片沿通道合并为多通道“图片”表示。若一次输入的颜色通道数量为C,则合并后的图片张量表示为[C,H,W]。
(2)图片标准化模块
本模块承接通道合并模块,即输入数据为合并后的多通道图片张量,符号表示为[C,H,W]。由于不同批次的输入数据有可能具有不同的高度H和宽度 W,本模块的作用就是将输入数据标准化为[C,45,30]的张量表示。具体方法为:
①使用双立方插值算法将[C,H,W]的图像张量转化为[C,45,30];
②将经过插值操作的图像张量作正则化操作:
x=x/255.0
其中x为图片张量。
(二)神经网络模块
神经网络模块用于承接图片标准化模块,输入数据为经过标准化的图片张量,具体表示为[C,45,30],经过神经网络得到最终的预测分类。
本发明基于深度学习的神经干细胞分化方向预测模型的运行流程如下:
①输入细胞图片进入图片预处理模块得到标准化的图片数据;
②标准化的图片数据直接进入神经网络模块,由神经网络根据数据自动提取有效的特征,得到最终的分类判断。
其中,分类判断包括:
标签 | 描述 |
0 | 诱导分化为星形胶质细胞 |
1 | 诱导分化为少突胶质细胞 |
2 | 诱导分化为神经元 |
基于以上构建的训练数据,已经训练完成的神经网络结构如下:
类型 | 卷积核(数量)尺寸/步长(或注释) |
模组1 | 子网络模块1 |
模组2 | 子网络模块2 |
模组3 | 子网络模块3 |
模组4 | 4x子网络模块4 |
模组5 | 子网络模块5 |
模组6 | 子网络模块6 |
池化 | (-)6x4/1 |
Softmax | 分类输出 |
子网络模块1
类型 | 卷积核(数量)尺寸/步长(或注释) |
卷积 | (32)3x3/1 |
卷积 | (64)3x3/1 |
子网络模块2
子网络模块3
子网络模块4
子网络模块5
子网络模块6
类型 | 卷积核(数量)尺寸/步长(或注释) |
卷积 | (1024)3x3/1 |
卷积 | (2048)3x3/1 |
接下来将介绍神经网络模型的数学表述,所采取的符号规则为:
模型子网络模块1(m1)的表达式为:
z1,1=Conv1,1(x);
a1,1=Relu(z1,1);
z1,2=Conv1,2(a1,1);
a1,2=Relu(z1,2);
u=a1,2。
模型子网络模块2(m2)的表达式为:
z2,1;s=Conv2,1;sx;
z2,1=DSConv2,1(x);
a2,1=Relu(z2,1);
z2,2=DSConv2,2(a2,1);
p2,1=Pool(z2,2);
u=p2,1+z2,1;s。
模型子网络模块3(m3)的表达式为:
z3,1;s=Conv3,1;sx;
z3,1=DSConv3,1(x);
a3,1=Relu(z3,1);
z3,2=DSConv3,2(a3,1);
p3,1=Pool(z3,2);
u=p3,1+z3,1;s。
模型子网络模块4(m4)的表达式为:
a4,1=Relu(x)
z4,1=DSConv4,1(a4,1);
a4,2=Relu(z4,1);
z4,2=DSConv4,2(a4,2);
a4,3=Relu(z4,2);
z4,3=DSConv4,3(a4,3);
u=z4,3+x。
模型子网络模块5(m5)的表达式为:
z5,1;s=Conv5,1;sx;
a5,1=Relu(x);
z5,1=DSConv5,1(a5,1);
a5,2=Relu(z5,1);
z5,2=DSConv5,2(a5,1);
p5,1=Pool(z5,2);
u=p5,1+z5,1;s。
模型子网络模块6(m6)的表达式为:
z6,1=Conv6,1(x);
a6,1=Relu(z6,1);
z6,2=Conv6,2(a6,1);
a6,2=Relu(z6,2);
u=a62。
模型整体表达式为:
u1=m1(x);
u2=m2(u1);
u3=m3(u2);
u4,1=m4,1(u3);
u4,2=m4,2(u4,1);
u4,3=m4,3(u4,2);
u4=m4,3(u4,3);
u5=m5(u4);
u6=m6(u5);
p=GPool(u6);
zw=Fc(p);
模型训练方法如下:
使用Pytorch框架在2块NVIDIA GTX 1080Ti显卡上训练神经网络。训练优化器为Adam优化器,相应的训练参数:学习率为0.001,beta1为0.9,beta2 为0.999,epsilon为1e-8。
模型训练测试流程示意图如图1所示。
下面是非独立测试集模型测试结果:
下面是非独立测试集免疫荧光染色实验统计结果和模型测试结果的对比:
下面是独立测试集模型测试结果:
下面给出另一个实例来说明本发明模型的性能。
另设计实验,用不同的诱导物诱导神经干细胞分化,包括神经因子(NT4, NGF,CNTF以及BDNF)和激素(MT)来处理细胞,并用全景流式细胞仪得到细胞图像,得到共计38,486个细胞图像,测试该模型的性能。
独立测试集诱导神经干细胞分化的方法如下:
下面是独立测试集免疫荧光染色实验统计结果和模型测试结果的对比:
结果显示我们训练得到的神经网络模型可以精准、快速、便捷地预测未标记的神经干细胞分化方向,并且可以在诱导物仅处理一天即预测出准确分化方向和分化比例,无需做荧光染色等试验操作。在传统依赖于分子标记手段无法识别时候,精准判断后期神经元分化比例,并在多种诱导物作用下都能取得很高的精确度,不受诱导物类型,以及诱导分化的信号通路和分子机制影响,具有优秀的应用前景。
实施例2本发明基于深度学习的神经干细胞分化方向预测方法
请参见图2,本实施例提供一种基于深度学习的神经干细胞分化方向预测方法,包括以下步骤:
S100,利用全景流式细胞仪收集神经干细胞分化培养的细胞图像,包括分化为神经元、星形胶质细胞和少突胶质细胞三个方向的细胞图像;
S200,将全景流式细胞仪收集的神经干细胞分化培养的细胞图像,输入卷积神经网络模型,由卷积神经网络模型自动提取有效的特征,训练得到神经干细胞分化方向预测模型;
S300,将待分析分化方向的神经干细胞,利用全景流式细胞仪收集其细胞图像并输入所述神经干细胞分化方向预测模型,得到预测结果。
所述卷积神经网络模型表达式为:
模型子网络模块1(m1)的表达式为:
z1,1=Conv1,1(x);
a1,1=Relu(z1,1);
z1,2=Conv1,2(a1,1);
a1,2=Relu(z1,2);
u=a1,2;
模型子网络模块2(m2)的表达式为:
z2,1;s=Conv2,1;sx;
z2,1=DSConv2,1(x);
a2,1=Relu(z2,1);
z2,2=DSConv2,2(a2,1);
p2,1=Pool(z2,2);
u=p2,1+z2,1;s;
模型子网络模块3(m3)的表达式为:
z3,1;s=Conv3,1;sx;
z3,1=DSConv3,1(x);
a3,1=Relu(z3,1);
z3,2=DSConv3,2(a3,1);
p3,1=Pool(z3,2);
u=P3,1+z3,1;s;
模型子网络模块4(m4)的表达式为:
a4,1=Relu(x);
z4,1=DSConv4,1(a4,1);
a4,2=Relu(z4,1);
z4,2=DSConv4,2(a4,2);
a4,3=Relu(z4,2);
z4,3=DSConv4,3(a4,3);
u=z4,3+x;
模型子网络模块5(m5)的表达式为:
z5,1;s=Conv5,1;sx;
a5,1=Relu(x);
z5,1=DSConv5,1(a5,1);
a5,2=Relu(z5,1);
z5,2=DSConv5,2(a5,1);
p5,1=Pool(z5,2);
u=p5,1+z5,1;s;
模型子网络模块6(m6)的表达式为:
z6,1=Conv6,1(x);
a6,1=Relu(z6,1);
z6,2=Conv6,2(a6,1);
a6,2=Relu(z6,2);
u=a62;
模型整体表达式为:
u1=m1(x);
u2=m2(u1);
u3=m3(u2);
u4,1=m4,1(u3);
u4,2=m4,2(u4,1);
u4,3=m4,3(u4,2);
u4=m4,3(u4,3);
u5=m5(u4);
u6=m6(u5);
p=GPool(u6);
zw=Fc(p);
以上符号的含义为:
具体网络结构为:
类型 | 卷积核(数量)尺寸/步长(或注释) |
模组1 | 子网络模块1 |
模组2 | 子网络模块2 |
模组3 | 子网络模块3 |
模组4 | 4x子网络模块4 |
模组5 | 子网络模块5 |
模组6 | 子网络模块6 |
池化 | (-)6x4/1 |
Softmax | 分类输出 |
子网络模块1
类型 | 卷积核(数量)尺寸/步长(或注释) |
卷积 | (32)3x3/1 |
卷积 | (64)3x3/1 |
子网络模块2
子网络模块3
子网络模块4
子网络模块5
子网络模块6
类型 | 卷积核(数量)尺寸/步长(或注释) |
卷积 | (1024)3x3/1 |
卷积 | (2048)3x3/1 |
步骤S100中,神经干细胞分化培养的培养基配方和采集细胞图像的时间点如下:
步骤S100中,训练集中的分化为星形胶质方向、神经元细胞方向、少突胶质细胞方向和NT3处理组的细胞图像的比例为 (1.9~2.1):(1.25~1.45):1:(1~1.1)。
所述神经干细胞分化方向预测模型包括图片预处理模块和神经网络模块;所述图片预处理模块用于接收全景流式细胞仪收集的神经干细胞分化培养的细胞图像,预处理得到标准化的图片数据;所述神经网络模块用于接收标准化的图片数据,得到最终的神经干细胞分化方向的预测结果。
所述图片预处理模块包括通道合并模块和图片标准化模块;所述通道合并模块用于输入具有相同的高度H和宽度W的单颜色通道细胞图像,将单通道图片沿通道合并为多通道“图片”表示,合并后的图片张量表示为[C,H,W],其中C表示输入的颜色通道数量;所述图片标准化模块用于承接通道合并模块,输入数据为合并后的多通道图片张量,符号表示为[C,H,W],将输入数据标准化为[C,45,30]的张量表示。
所述图片标准化模块使用双立方插值算法将[C,H,W]的图像张量转化为 [C,45,30],将经过插值操作的图像张量作正则化操作。
所述神经干细胞分化方向的预测结果还包括神经干细胞的分化比例。
实施例3本发明基于深度学习的神经干细胞分化方向预测系统
请参见图3,本实施例提供一种基于深度学习的神经干细胞分化方向预测系统,包括:
全景流式细胞仪100,用于获取神经干细胞分化培养的细胞图像。全景流式细胞仪是本领域常规实验设备之一,它不仅能够获得大量细胞的群体分析数据,而且还能够实时看到细胞图像,让每一步的分析结果都可以通过图像进行确认。较主流的全景流式细胞仪设备有德国Merck Millipore公司的FlowSight 多维全景流式细胞仪等。
神经干细胞分化方向预测模型200,用于接收全景流式细胞仪收集的待预测分化方向的神经干细胞分化培养的细胞图像,得到预测结果;所述神经干细胞分化方向预测模型是按照以下方法构建的:
S100,利用全景流式细胞仪收集神经干细胞分化培养的细胞图像,包括分化为神经元、星形胶质细胞和少突胶质细胞三个方向的细胞图像;
S200,将全景流式细胞仪收集的神经干细胞分化培养的细胞图像,输入卷积神经网络模型,由卷积神经网络模型自动提取有效的特征,训练得到神经干细胞分化方向预测模型;
所述卷积神经网络模型表达式为:
模型子网络模块1(m1)的表达式为:
z1,1=Conv1,1(x);
a1,1=Relu(z1,1);
z1,2=Conv1,2(a1,1);
a1,2=Relu(z1,2);
u=a1,2;
模型子网络模块2(m2)的表达式为:
z2,1;s=Conv2,1;sx;
z2,1=DSConv2,1(x);
a2,1=Relu(z2,1);
z2,2=DSConv2,2(a2,1);
p2,1=Pool(z2,2);
u=p2,1+z2,1;s;
模型子网络模块3(m3)的表达式为:
z3,1;s=Conv3,1;sx;
z3,1=DSConv3,1(x);
a3,1=Relu(z3,1);
z3,2=DSConv3,2(a3,1);
p3,1=Pool(z3,2);
u=p3,1+z3,1;s;
模型子网络模块4(m4)的表达式为:
a4,1=Relu(x);
z4,1=DSConv4,1(a4,1);
a4,2=Relu(z4,1);
z4,2=DSConv4,2(a4,2);
a4,3=Relu(z4,2);
z4,3=DSConv4,3(a4,3);
u=z4,3+x;
模型子网络模块5(m5)的表达式为:
z5,1;s=Conv5,1;sx;
a5,1=Relu(x);
z5,1=DSConv5,1(a5,1);
a5,2=Relu(z5,1);
z5,2=DSConv5,2(a5,1);
p5,1=Pool(z5,2);
u=p5,1+z5,1;s;
模型子网络模块6(m6)的表达式为:
z6,1=Conv6,1(x);
a6,1=Relu(z6,1);
z6,2=Conv6,2(a6,1);
a6,2=Relu(z6,2);
u=a62;
模型整体表达式为:
u1=m1(x);
u2=m2(u1);
u3=m3(u2);
u4,1=m4,1(u3);
u4,2=m4,2(u4,1);
u4,3=m4,3(u4,2);
u4=m4,3(u4,3);
u5=m5(u4);
u6=m6(u5);
p=GPool(u6);
zw=Fc(p);
以上符号的含义为:
具体网络结构为:
类型 | 卷积核(数量)尺寸/步长(或注释) |
模组1 | 子网络模块1 |
模组2 | 子网络模块2 |
模组3 | 子网络模块3 |
模组4 | 4x子网络模块4 |
模组5 | 子网络模块5 |
模组6 | 子网络模块6 |
池化 | (-)6x4/1 |
Softmax | 分类输出 |
子网络模块1
类型 | 卷积核(数量)尺寸/步长(或注释) |
卷积 | (32)3x3/1 |
卷积 | (64)3x3/1 |
子网络模块2
子网络模块3
子网络模块4
子网络模块5
子网络模块6
类型 | 卷积核(数量)尺寸/步长(或注释) |
卷积 | (1024)3x3/1 |
卷积 | (2048)3x3/1 |
步骤S100中,神经干细胞分化培养的培养基配方和采集细胞图像的时间点如下:
步骤S100中,训练集中的分化为星形胶质方向、神经元细胞方向、少突胶质细胞方向和NT3处理组的细胞图像的比例为 (1.9~2.1):(1.25~1.45):1:(1~1.1)。
所述神经干细胞分化方向预测模型包括图片预处理模块201和神经网络模块202。
所述图片预处理模块201用于接收全景流式细胞仪100获取的神经干细胞分化培养的细胞图像,预处理得到标准化的图片数据。具体地,所述图片预处理模块201包括通道合并模块2011和图片标准化模块2012。所述通道合并模块2011用于输入单颜色通道细胞图像,这些单颜色通道细胞图像必须具有相同的高度H和宽度W;通道合并模块2011将这些单通道图片沿通道合并为多通道“图片”表示;若一次输入的颜色通道数量为C,则合并后的图片张量表示为[C,H,W]。所述图片标准化模块2012用于承接通道合并模块2011,输入数据为合并后的多通道图片张量,符号表示为[C,H,W]。由于不同批次的输入数据有可能具有不同的高度H和宽度W,本模块的作用就是将输入数据标准化为[C,45,30]的张量表示。具体方法为:①使用双立方插值算法将[C,H,W]的图像张量转化为[C,45,30];②将经过插值操作的图像张量作正则化操作。
所述神经网络模块202用于接收标准化的图片数据,得到最终的神经干细胞分化方向预测结果。
作为一个优选例,神经干细胞分化方向的预测结果还包括神经干细胞的分化比例。
作为另一优选例,待预测分化方向的神经干细胞分化培养的细胞图像获取时间点为神经干细胞分化培养0.5-1d。
作为另一优选例,本发明基于深度学习的神经干细胞分化方向预测系统还包括显示器300,用于显示神经干细胞分化方向判断和分化比例判断结果。
作为另一优选例,本发明基于深度学习的神经干细胞分化方向预测系统还包括通信模块400,用于与全景流式细胞仪100连接,从而获取神经干细胞分化培养的细胞图像。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员,在不脱离本发明方法的前提下,还可以做出若干改进和补充,这些改进和补充也应视为本发明的保护范围。
Claims (7)
1.一种基于深度学习的神经干细胞分化方向预测方法,其特征在于,包括以下步骤:
S100,利用全景流式细胞仪收集神经干细胞分化培养的细胞图像,包括分化为神经元、星形胶质细胞和少突胶质细胞三个方向的细胞图像;
S200,将全景流式细胞仪收集的神经干细胞分化培养的细胞图像,输入卷积神经网络模型,由卷积神经网络模型自动提取有效的特征,训练得到神经干细胞分化方向预测模型;
S300,将待分析分化方向的神经干细胞,利用全景流式细胞仪收集其细胞图像并输入所述神经干细胞分化方向预测模型,得到预测结果,
所述卷积神经网络模型表达式为:
模型子网络模块1,m1的表达式为:
z1,1=Conv1,1(x);
a1,1=Relu(z1,1);
z1,2=Conv1,2(a1,1);
a1,2=Relu(z1,2);
u=a1,2;
模型子网络模块2,m2的表达式为:
z2,1;s=Conv2,1;sx;
z2,1=DSConv2,1(x);
a2,1=Relu(z2,1);
z2,2=DSConv2,2(a2,1);
p2,1=Pool(z2,2);
u=p2,1+z2,1;s;
模型子网络模块3,m3的表达式为:
z3,1;s=Conv3,1;sx;
z3,1=DSConv3,1(x);
a3,1=Relu(z3,1);
z3,2=DSConv3,2(a3,1):
p3,1=Pool(z3,2);
u=p3,1+z3,1;s;
模型子网络模块4,m4的表达式为:
a4,1=Relu(x);
z4,1=DSConv4,1(a4,1);
a4,2=Relu(z4,1);
z4,2=DSConv4,2(a4,2);
a4,3=Relu(z4,2);
z4,3=DSConv4,3(a4,3);
u=z4,3+x;
模型子网络模块5,m5的表达式为:
z5,1;s=Conv5,1;sx;
a5,1=Relu(x);
z5,1=DSConv5,1(a5,1);
a5,2=Relu(z5,1);
z5,2=DSConv5,2(a5,1);
p5,1=Pool(z5,2);
u=p5,1+z5,1;s;
模型子网络模块6,m6的表达式为:
z6,1=Conv6,1(x);
a6,1=Relu(z6,1);
z6,2=Conv6,2(a6,1);
a6,2=Relu(z6,2);
u=a62;
模型整体表达式为:
u1=m1(x);
u2=m2(u1);
u3=m3(u2);
u4,1=m4,1(u3);
u4,2=m4,2(u4,1);
u4,3=m4,3(u4,2);
u4=m4,3(u4,3);
u5=m5(u4);
u6=m6(u5);
p=GPool(u6);
zw=Fc(p);
以上符号的含义如下:
m表示子网络模块函数;x表示输入张量;u表示子网络模块输出向量;表示输出张量;z表示卷积运算结果;zw表示全连接层输出结果;p表示池化运算结果;a表示激活值;Fc表示全连接层;Conv表示卷积层;Pool表示池化层;Relu表示Relu激活函数;DSConv表示深度可分卷积层;GPool表示全局池化层;Softmax表示Softmax函数,
步骤S100中,神经干细胞分化培养的培养基配方和采集细胞图像的时间点如下:
步骤S100中,训练集中的分化为星形胶质方向、神经元细胞方向、少突胶质细胞方向和NT3处理组的细胞图像的比例为(1.9~2.1)∶(1.25~1.45)∶1∶(1~1.1)。
2.根据权利要求1所述的基于深度学习的神经干细胞分化方向预测方法,其特征在于,所述神经干细胞分化方向预测模型包括图片预处理模块和神经网络模块;所述图片预处理模块用于接收全景流式细胞仪收集的神经干细胞分化培养的细胞图像,预处理得到标准化的图片数据;所述神经网络模块用于接收标准化的图片数据,得到最终的神经干细胞分化方向的预测结果。
3.根据权利要求2所述的基于深度学习的神经干细胞分化方向预测方法,其特征在于,所述图片预处理模块包括通道合并模块和图片标准化模块;所述通道合并模块用于输入具有相同的高度H和宽度W的单颜色通道细胞图像,将单通道图片沿通道合并为多通道“图片”表示,合并后的图片张量表示为[C,H,W],其中C表示输入的颜色通道数量;所述图片标准化模块用于承接通道合并模块,输入数据为合并后的多通道图片张量,符号表示为[C,H,W],将输入数据标准化为[C,45,30]的张量表示。
4.根据权利要求3所述的基于深度学习的神经干细胞分化方向预测方法,其特征在于,所述图片标准化模块使用双立方插值算法将[C,H,W]的图像张量转化为[C,45,30],将经过插值操作的图像张量作正则化操作。
5.根据权利要求1所述的基于深度学习的神经干细胞分化方向预测方法,其特征在于,所述神经干细胞分化方向的预测结果还包括神经干细胞的分化比例。
6.一种基于深度学习的神经干细胞分化方向预测系统,其特征在于,包括:
全景流式细胞仪,用于收集神经干细胞分化培养的细胞图像;
神经干细胞分化方向预测模型,用于接收全景流式细胞仪收集的待预测分化方向的神经干细胞分化培养的细胞图像,得到预测结果;所述神经干细胞分化方向预测模型是按照以下方法构建的:
S100,利用全景流式细胞仪收集神经干细胞分化培养的细胞图像,包括分化为神经元、星形胶质细胞和少突胶质细胞三个方向的细胞图像;
S200,将全景流式细胞仪收集的神经干细胞分化培养的细胞图像,输入卷积神经网络模型,由卷积神经网络模型自动提取有效的特征,训练得到神经干细胞分化方向预测模型;
所述卷积神经网络模型表达式为:
模型子网络模块1,m1的表达式为:
z1,1=Conv1,1(x);
a1,1=Relu(z1,1);
z1,2=Conv1,2(a1,1);
a1,2=Relu(z1,2);
u=a1,2;
模型子网络模块2,m2的表达式为:
z2,1;s=Conv2,1;sx;
z2,1=DSConv2,1(x);
a2,1=Relu(z2,1);
z2,2=DSConv2,2(a2,1);
p2,1=Pool(z2,2);
u=p2,1+z2,1;s;
模型子网络模块3,m3的表达式为:
z3,1;s=Conv3,1;sx;
z3,1=DSConv3,1(x);
33,1=Relu(z3,1);
z3,2=DSConv3,2(a3,1);
p3,1=Pool(z3,2);
u=p3,1+z3,1;s;
模型子网络模块4,m4的表达式为:
a4,1=Relu(x);
z4,1=DSConv4,1(a4,1):
a4,2=Relu(z4,1);
z4,2=DSConv4,2(a4,2);
34,3=Relu(z4,2);
z4,3=DSConv4,3(a4,3);
u=z4,3+x;
模型子网络模块5,m5的表达式为:
z5,1;s=Conv5,1;sx;
a5,1=Relu(x);
z5,1=DSConv5,1(a51);
a5,2=Relu(z5,1);
z5,2=DSConv5,2(a5,1);
p5,1=Pool(z5,2);
u=p5,1+z5,1;s;
模型子网络模块6,m6的表达式为:
z6,1=Conv6,1(x);
a6,1=Relu(z6,1);
z6,2=Conv6,2(a6,1);
a6,2=Relu(z6,2);
u=a62;
模型整体表达式为:
u1=m1(x);
u2=m2(u1);
u3=m3(u2);
u4,1=m4,1(u3);
u4,2=m4,2(u4,1);
u4,3=m4,3(u4,2);
u4=m4,3(u4,3);
u5=m5(u4);
u6=m6(u5);
p=GPool(u6);
zw=Fc(p);
以上符号的含义如下:
m表示子网络模块函数;x表示输入张量;u表示子网络模块输出向量;表示输出张量;z表示卷积运算结果;表示全连接层输出结果;p表示池化运算结果;a表示激活值;Fc表示全连接层;Conv表示卷积层;Pool表示池化层;Relu表示Relu激活函数;DSConv表示深度可分卷积层;GPool表示全局池化层;Softmax表示Softmax函数,
所述步骤S100中,神经干细胞分化培养的培养基配方和采集细胞图像的时间点如下:
所述,步骤S100中,训练集中的分化为星形胶质方向、神经元细胞方向、少突胶质细胞方向和NT3处理组的细胞图像的比例为(1.9~2.1):(1.25~1.45):1:(1~1.1)。
7.根据权利要求6所述的基于深度学习的神经干细胞分化方向预测系统,其特征在于,所述神经干细胞分化方向预测模型包括图片预处理模块和神经网络模块;所述图片预处理模块用于接收全景流式细胞仪收集的神经干细胞分化培养的细胞图像,预处理得到标准化的图片数据;所述图片预处理模块包括通道合并模块和图片标准化模块;所述通道合并模块用于输入具有相同的高度H和宽度W的单颜色通道细胞图像,将单通道图片沿通道合并为多通道“图片”表示,合并后的图片张量表示为[C,H,W],其中C表示输入的颜色通道数量;所述图片标准化模块用于承接通道合并模块,输入数据为合并后的多通道图片张量,符号表示为[C,H,W],将输入数据标准化为[C,45,30]的张量表示;所述神经网络模块用于接收标准化的图片数据,得到最终的神经干细胞分化方向的预测结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010514541.1A CN111666895B (zh) | 2020-06-08 | 2020-06-08 | 基于深度学习的神经干细胞分化方向预测系统及方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010514541.1A CN111666895B (zh) | 2020-06-08 | 2020-06-08 | 基于深度学习的神经干细胞分化方向预测系统及方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111666895A CN111666895A (zh) | 2020-09-15 |
CN111666895B true CN111666895B (zh) | 2023-05-26 |
Family
ID=72385801
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010514541.1A Active CN111666895B (zh) | 2020-06-08 | 2020-06-08 | 基于深度学习的神经干细胞分化方向预测系统及方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111666895B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114121161B (zh) * | 2021-06-04 | 2022-08-05 | 深圳太力生物技术有限责任公司 | 一种基于迁移学习的培养基配方开发方法及系统 |
CN113850767A (zh) * | 2021-09-14 | 2021-12-28 | 北京大学 | 基于深度学习的细胞牵引力测量方法、装置、设备及介质 |
CN115880689A (zh) * | 2021-09-26 | 2023-03-31 | 瑞新(福州)科技有限公司 | 细胞识别的方法、装置和系统 |
WO2023221951A2 (zh) * | 2022-05-14 | 2023-11-23 | 北京大学 | 基于细胞动态图像机器学习的细胞分化 |
Family Cites Families (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2008100168A1 (en) * | 2007-02-16 | 2008-08-21 | Universidade De Coimbra | Method for the functional identification of new neurons, neural progenitors, astrocytes and immature cells from stem cell cultures and uses thereof |
CA2926575C (en) * | 2013-10-14 | 2023-01-31 | Hadasit Medical Research Services & Development Limited | Method of obtaining terminally differentiated neuronal lineages and uses thereof |
SG11201603414YA (en) * | 2013-11-11 | 2016-05-30 | C O Agency For Science Technology And Res | Method for differentiating induced pluripotent stem cells into renal proximal tubular cell-like cells |
CN108138130A (zh) * | 2015-08-31 | 2018-06-08 | 爱平世股份有限公司 | 多能干细胞制造系统和生产诱导多能干细胞的方法 |
US10282588B2 (en) * | 2016-06-09 | 2019-05-07 | Siemens Healthcare Gmbh | Image-based tumor phenotyping with machine learning from synthetic data |
CN106202997B (zh) * | 2016-06-29 | 2018-10-30 | 四川大学 | 一种基于深度学习的细胞分裂检测方法 |
US10282589B2 (en) * | 2017-08-29 | 2019-05-07 | Konica Minolta Laboratory U.S.A., Inc. | Method and system for detection and classification of cells using convolutional neural networks |
CN108280320B (zh) * | 2018-01-23 | 2020-12-29 | 上海市同济医院 | 一种基于深度学习的快速精准高通量药物筛选系统 |
US20210254049A1 (en) * | 2018-04-20 | 2021-08-19 | Cellino Biotech, Inc. | Directed cell fate specification and targeted maturation |
CN109360193A (zh) * | 2018-09-27 | 2019-02-19 | 北京基石生命科技有限公司 | 一种基于深度学习的原代肿瘤细胞分割识别方法及系统 |
CN110334565A (zh) * | 2019-03-21 | 2019-10-15 | 江苏迪赛特医疗科技有限公司 | 一种显微镜病理照片的宫颈癌病变细胞分类系统 |
CN110136775A (zh) * | 2019-05-08 | 2019-08-16 | 赵壮志 | 一种细胞分裂及抗干扰检测系统及方法 |
-
2020
- 2020-06-08 CN CN202010514541.1A patent/CN111666895B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN111666895A (zh) | 2020-09-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111666895B (zh) | 基于深度学习的神经干细胞分化方向预测系统及方法 | |
CN108280320B (zh) | 一种基于深度学习的快速精准高通量药物筛选系统 | |
Pepperkok et al. | High-throughput fluorescence microscopy for systems biology | |
Doan et al. | Leveraging machine vision in cell-based diagnostics to do more with less | |
Yao et al. | Cell type classification and unsupervised morphological phenotyping from low-resolution images using deep learning | |
Guo et al. | Breast cancer histology image classification based on deep neural networks | |
CN112784767A (zh) | 基于白细胞显微图像的细胞实例分割算法 | |
Zhou et al. | Computer vision meets microfluidics: a label-free method for high-throughput cell analysis | |
CN110807426B (zh) | 一种基于深度学习的寄生虫检测系统及方法 | |
CN113241154A (zh) | 一种人工智能血液涂片细胞标注系统及方法 | |
Majidpour et al. | NSGA-II-DL: Metaheuristic optimal feature selection with Deep Learning Framework for HER2 classification in Breast Cancer | |
Kabeya et al. | Deep convolutional neural network-based algorithm for muscle biopsy diagnosis | |
CN108319815A (zh) | 一种用于细胞虚拟染色的方法及其系统 | |
CN114864075A (zh) | 一种基于病理图像的胶质瘤级别分析方法及装置 | |
Cho et al. | Numerical learning of deep features from drug-exposed cell images to calculate IC50 without staining | |
Li et al. | CellBin enables highly accurate single-cell segmentation for spatial transcriptomics | |
CN116757998A (zh) | 基于ai的ctc细胞和ctc样细胞的筛查方法及装置 | |
Calvert et al. | Yeast cell cycle analysis: combining DNA staining with cell and nuclear morphology | |
Chu et al. | High resolution U-Net for quantitatively analyzing early spatial patterning of human induced pluripotent stem cells on micropatterns | |
Smith et al. | Machine vision and automation in secondary metabolite bioprocess control | |
CN111896456A (zh) | 一种基于微流控和高光谱成像的单细胞分析方法 | |
Fernandez et al. | Analytical validation of the PreciseDx digital prognostic breast cancer test in early-stage breast cancer | |
Rappez et al. | Label-free phenotyping of human microvessel networks | |
Wang et al. | Induced pluripotent stem cells detection via ensemble Yolo network | |
Mukhopadhyay et al. | Machine learning aided single cell image analysis improves understanding of morphometric heterogeneity of human mesenchymal stem cells |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |