Diagnostic support device, machine learning device, diagnostic support method, machine learning method, machine learning program, and Alzheimer's prediction program

Abstract

The possibility that an ADNC patient will develop Alzheimer's disease is predicted with high precision.　A diagnosis assistance device 2 that predicts the possibility that a test subject having ADNC will develop Alzheimer's disease within a prescribed period, said diagnosis assistance device 2 comprising a prediction unit 23 that predicts, according to a machine-learned prediction algorithm D4, the possibility that the test subject will develop Alzheimer's disease within the prescribed period .

Description

[Technical Field]

[0001] The present invention relates to a technology for predicting whether a patient with ADNC (Alzheimer's disease neuropathologic change) will develop Alzheimer's disease within a predetermined period of time, and in particular to a technology for making such predictions using artificial intelligence. [Background technology]

[0002] To treat Alzheimer's disease (AD), it is necessary to detect the onset of the disease early. It is desirable to develop a technology for diagnosing the disease before it develops. To address this issue, VSRAD (registered trademark) (Voxel-Based Specific Regional Analysis System for Alzheimer's Disease) has been developed as a diagnostic support system for early AD (Alzheimer's Disease) (Patent Reference 1). VSRAD (registered trademark) is image processing and statistical analysis software for reading from MRI images the degree of atrophy around the parahippocampal gyrus, which is characteristic of early AD, including the prodromal stage. Furthermore, the present inventors have investigated whether patients with mild cognitive impairment (MCI), which is considered a precursor to Alzheimer's disease, can develop progressive pMCI, which progresses to Alzheimer's disease. We have developed a technology to predict whether a patient has progressive MCI (progressive MCI) or non-progressive sMCI (stable MCI) that does not progress to Alzheimer's disease (Patent Document 2).

[0003] MCI is a concept that lies on the border between normal (NL) and Alzheimer's disease, but recent guidelines state that setting a cutoff value is not recommended for any test. In other words, it is difficult to clearly distinguish between normal and MCI, and between MCI and Alzheimer's disease.

[0004] Meanwhile, in recent years, in vivo pathological diagnosis has been used to detect ADNC (Alzheimer's disease neuropathologic changes) in order to predict the onset of Alzheimer's disease. ADNC is identified by the presence of both senile plaque (amyloid beta) deposition and neurofibrillary tangles (tau degeneration). Currently, amyloid beta deposition and tau degeneration can be detected by cerebrospinal fluid testing. In the future, it is expected that amyloid beta accumulation and tau degeneration will also be detectable by blood testing, and amyloid PET and tau PET have also been developed. As such, it is now possible to clearly determine whether or not a patient has ADNC. From the perspective of early treatment, it is believed that ADNC diagnosis will become more important than clinical diagnosis in the future.

[0005] Development of disease-modifying therapies (DMTs) that reduce the accumulation of amyloid beta is also progressing. For example, in the fall of 2019, it was reported that aducanumab, developed by Biogen and Eisai, slowed the progression of MCI and early Alzheimer's disease. [Prior art documents] [Patent documents]

[0006] [Patent Document 1] International Publication No. 2013 / 047278 [Patent Document 2] Patent No. 6483890 Summary of the Invention [Problem to be solved by the invention]

[0007] It is known that some patients diagnosed with ADNC do not necessarily develop Alzheimer's disease. Therefore, it is not appropriate to target all ADNC patients for DMT, and it is currently unclear which patients should receive DMT and when.

[0008] Therefore, an object of the present invention is to predict with high accuracy the possibility that an ADNC patient will develop Alzheimer's disease within a predetermined period of time. [Means for solving the problem]

[0009] The present invention includes the following aspects. Section 1. A diagnostic support device that predicts the likelihood that a subject with Alzheimer's neuropathological changes will develop Alzheimer's disease within a predetermined period, comprising: A diagnostic support device comprising a prediction unit that makes the prediction according to a machine-learned prediction algorithm. Section 2. Item 1, a diagnosis support device, a region segmentation unit that segments the brain image acquired from the subject into gray matter, white matter, and cerebrospinal fluid portions and separates the lateral ventricles from the cerebrospinal fluid portions; a region of interest setting unit that sets a plurality of regions of interest in each of the gray matter, the white matter, and the lateral ventricle; a t-value and p-value calculation unit that calculates t-values ​​and p-values ​​in each region of interest for the volume of each region of interest; a z-value calculation unit that calculates a z-value for each region of interest based on the t-value and p-value; Furthermore, The prediction unit performs the prediction based on the z value. Section 3. Item 2. The diagnosis support device according to item 2, The region segmentation unit separates the corpus callosum from the surrounding white matter by determining the boundary between the corpus callosum and the surrounding white matter using parameters of surface tension and viscosity of a fluid. Section 4. Item 2 or 3, the diagnosis support device, When a white matter lesion is present in the white matter, the region of interest setting unit extracts the white matter lesion and replaces it with an average signal value of the white matter of the subject, and then sets the region of interest in the white matter. Section 5. Item 1, a diagnosis support device, a region separation unit that separates gray matter from the brain image acquired from the subject; a region of interest setting unit that sets a plurality of regions of interest in the gray matter; a volume calculation unit that calculates the volume of each region of interest; a z-value calculation unit that calculates a z-value of each region of interest based on the volume; Furthermore, The prediction unit performs the prediction based on the z value. Section 6. Item 5. The diagnosis support device according to item 5, The diagnosis support device further comprises a covariate correction unit that performs covariate correction on the volume. Section 7. The diagnosis support device according to any one of items 2 to 6, The prediction unit performs the prediction as a posterior probability from a distance to a hyperplane using a sigmoid function. Section 8. A machine learning device that learns the prediction algorithm according to any one of items 1 to 7, A machine learning device comprising a learning unit that learns the prediction algorithm based on training data created from brain images of multiple people and diagnostic results indicating whether each person developed Alzheimer's disease before the specified period has elapsed since the acquisition of the brain images. Section 9. Item 9. The machine learning device according to item 8, The learning unit is a machine learning device configured with a support vector machine. Section 10. Item 8 or 9, the machine learning device according to item 8 or 9, The machine learning device, wherein the brain image is an MRI image. Section 11. The machine learning device according to any one of items 8 to 10, The machine learning device further includes a teacher data creation unit that creates the teacher data based on brain images of the plurality of people and diagnostic results indicating whether each person developed Alzheimer's disease before the specified period of time has elapsed since the brain images were acquired. Section 12. Item 12. The machine learning device according to item 11, The teacher data creation unit a region segmentation unit that segments a brain image acquired from the person into gray matter, white matter, and cerebrospinal fluid portions and separates a lateral ventricle from the cerebrospinal fluid portion; a region of interest setting unit that sets a plurality of regions of interest in each of the gray matter, the white matter, and the lateral ventricle; a t-value and p-value calculation unit that calculates t-values ​​and p-values ​​in each region of interest for the volume of each region of interest; a z-value calculation unit that calculates a z-value for each region of interest based on the t-value and p-value; Equipped with A machine learning device, wherein the training data includes the diagnosis result and the z-value. Section 13. Item 12. The machine learning device according to item 11, The teacher data creation unit a region separation unit that separates gray matter from a brain image acquired from the person; a region of interest setting unit that sets a plurality of regions of interest in the gray matter; a volume calculation unit that calculates the volume of each region of interest; a z-value calculation unit that calculates a z-value of each region of interest based on the volume; Equipped with A machine learning device, wherein the training data includes the diagnosis result and the z-value. Section 14. Item 14. The machine learning device according to item 13, The machine learning device further comprises a covariate correction unit that performs covariate correction on the volume. Section 15. A diagnostic support method for predicting the likelihood that a subject with Alzheimer's disease neuropathological changes will develop Alzheimer's disease within a predetermined period, comprising: A diagnostic support method comprising a prediction step of making the prediction according to a machine-learned prediction algorithm. Section 16. A machine learning method for learning the prediction algorithm according to item 12, A machine learning method comprising a learning step of learning the prediction algorithm based on training data created from brain images of multiple people and diagnostic results indicating whether each person developed Alzheimer's disease before the specified period of time has passed since the acquisition of the brain images. Section 17. Item 16. A machine learning program that causes a computer to learn the prediction algorithm according to item 15, A machine learning program that causes the computer to execute a learning step of learning the prediction algorithm based on training data created from brain images of multiple people and diagnostic results indicating whether each person developed Alzheimer's disease before the specified period of time has passed since the acquisition of the brain images. Section 18. a training data creation step of creating training data from brain images of a plurality of people and diagnostic results indicating whether or not each person has developed Alzheimer's disease before the predetermined period has elapsed since the acquisition of the brain images; a learning step of learning a prediction algorithm based on the training data; a prediction step of predicting the likelihood that a subject with Alzheimer's neuropathological changes will develop Alzheimer's disease within a predetermined period of time according to the prediction algorithm; An Alzheimer's disease prediction program that causes a computer to execute the following: The teacher data creation step includes: isolating grey matter from brain images obtained from said person; establishing a plurality of regions of interest in the gray matter; calculating a volume for each region of interest; calculating a z-score for each region of interest based on the volume; creating the training data by associating the diagnosis results with the z values; Equipped with The prediction step isolating gray matter from brain images acquired from the subject; establishing a plurality of regions of interest in the gray matter; calculating a volume for each region of interest; calculating a z-score for each region of interest based on the volume; making the prediction based on the z-value; An Alzheimer's prediction program. [Effects of the Invention]

[0010] According to the present invention, it is possible to predict with high accuracy the possibility that an ADNC patient will develop Alzheimer's disease. [Brief explanation of the drawings]

[0011] [Figure 1] 1 is a block diagram showing a schematic configuration of a prediction system according to a first embodiment of the present invention. [Figure 2] FIG. 1 is a block diagram showing functions of a machine learning device according to a first embodiment of the present invention. [Figure 3] 1 is a flowchart showing the overall procedure of a machine learning method according to a first embodiment of the present invention. [Figure 4] 1 is a flowchart showing the procedure of a teacher data creation step in a machine learning method according to a first embodiment of the present invention. [Figure 5] 1 is a specific example of a brain image segmentation method. [Figure 6] FIG. 10 is an explanatory diagram of the effect of brain image segmentation. [Figure 7] FIG. 10 is an explanatory diagram of the effect of brain image segmentation. [Figure 8] FIG. 10 is an explanatory diagram of the effect of brain image segmentation. [Figure 9] FIG. 10 is an explanatory diagram of the effect of brain image segmentation. [Figure 10] FIG. 1 is an explanatory diagram of the effect of lateral ventricle isolation. [Figure 11] This is an explanatory diagram of an example in which the three-dimensional structure of the corpus callosum is determined and its boundaries are clarified. [Figure 12] 1 is a block diagram showing functions of a diagnosis support device according to a first embodiment of the present invention. [Figure 13]1 is a graph showing the time to onset of AD in pMCI patients. [Figure 14] FIG. 10 is a block diagram showing a schematic configuration of a prediction system according to a second embodiment of the present invention. [Figure 15] FIG. 10 is a block diagram showing the functions of a machine learning device according to a second embodiment of the present invention. [Figure 16] 10 is a flowchart showing the overall procedure of a machine learning method according to a second embodiment of the present invention. [Figure 17] 10 is a flowchart showing the procedure of a teacher data creation step in a machine learning method according to a second embodiment of the present invention. [Figure 18] 10A is a flowchart showing the procedure of the region separation step, and FIG. 10B is a flowchart showing the procedure of the image correction step. [Figure 19] FIG. 10 is a block diagram showing functions of a diagnosis support device according to a second embodiment of the present invention. [Figure 20] This is a graph showing the relationship between the number of years elapsed and the rate of developing AD for each group classified by cerebrospinal fluid tests and prediction results from a prediction algorithm. DETAILED DESCRIPTION OF THE INVENTION

[0012] [First embodiment] A first embodiment of the present invention will be described below with reference to the accompanying drawings, but the present invention is not limited to the following embodiment.

[0013] (Overall composition) FIG. 1 is a block diagram showing a schematic configuration of a prediction system 100 according to this embodiment. The prediction system 100 includes a machine learning device 1 and a diagnostic support device 2. The machine learning device 1 learns a prediction algorithm for predicting the likelihood that a subject with Alzheimer's disease neuropathologic change (ADNC) (hereinafter referred to as an ADNC subject) will develop Alzheimer's disease within a predetermined period of time. The diagnostic support device 2 predicts the likelihood that an ADNC subject will develop Alzheimer's disease within a predetermined period of time according to the prediction algorithm learned by the machine learning device 1. The machine learning device 1 and the diagnostic support device 2 may be implemented as separate devices, or the machine learning device 1 and the diagnostic support device 2 may be implemented as a single device.

[0014] Below, configuration examples of the machine learning device 1 and the diagnosis support device 2 will be described.

[0015] (machine learning device) FIG. 2 is a block diagram showing the functions of the machine learning device 1 according to this embodiment. The machine learning device 1 can be configured, for example, as a general-purpose personal computer, and includes, as its hardware configuration, a CPU (not shown), a main memory device (not shown), and an auxiliary memory device 11. In the machine learning device 1, the CPU reads various programs stored in the auxiliary memory device 11 into the main memory device and executes them, thereby performing various arithmetic processing. The auxiliary memory device 11 can be configured, for example, as a hard disk drive (HDD) or a solid-state drive (SSD). The auxiliary memory device 11 may be built into the machine learning device 1, or may be provided as an external memory device separate from the machine learning device 1.

[0016] The machine learning device 1 has a function of learning a prediction algorithm D4 for predicting the possibility that an ADNC subject will develop Alzheimer's disease within a predetermined period (for example, within 5 years). Patients who have already developed Alzheimer's disease (AD) Patients with progressive mild dementia (pMCI) who will develop Alzheimer's disease within a specified period Patients with mild, non-progressive dementia (sMCI) who will not develop Alzheimer's disease in the future In this embodiment, among ADNC patients, AD and pMCI are referred to as the "ADNC spectrum." That is, the machine learning device 1 has a function of learning prediction algorithms D4 and D5 for predicting the possibility that an ADNC subject is in the ADNC spectrum.

[0017] To achieve this function, the machine learning device 1 includes functional blocks: a teacher data creation unit 12 and a learning unit 13. The teacher data creation unit 12 is a functional block that creates teacher data D3 from brain images D1 and diagnosis results D2 of multiple people. The multiple people are preferably patients diagnosed with ADNC, but are not limited to this and may also include people diagnosed with mild cognitive impairment and healthy individuals. The learning unit 13 is a functional block that learns prediction algorithms D4 and D5 based on the teacher data D3. The teacher data creation unit 12 and the learning unit 13 are realized by executing a machine learning program stored in the auxiliary storage device 11.

[0018] The machine learning device 1 is capable of accessing a diagnostic information database DB. The diagnostic information database DB stores brain images D1 of multiple people and diagnostic results D2 indicating whether or not each person is in the ADNC spectrum. The diagnostic results D2 are diagnostic results indicating whether or not Alzheimer's disease has developed before a predetermined period of time has elapsed since the acquisition of the brain image D1. "Before a predetermined period of time has elapsed since the acquisition of the brain image D1" includes not only the period from the acquisition of the brain image D1 until the predetermined period of time has elapsed, but also the period before the acquisition of the brain image D1. nothing That is, the diagnosis result D2 is a diagnosis result indicating whether each person has developed AD within a predetermined period of time since the acquisition of the brain image D1. Not only Furthermore, the results of the diagnosis indicating whether each person had AD at the time of acquisition of the brain image D1 are included. nothing .

[0019] In this embodiment, the brain image D1 is a three-dimensional MRI image. It is desirable to prepare a certain number of brain images D1 and diagnostic results D2 for each of the Alzheimer's disease patient group, the Alzheimer's disease patient group, and the non-Alzheimer's disease patient group, such that a statistically significant difference can be obtained. The diagnostic information database DB may be owned by a single medical institution or may be shared by multiple medical institutions.

[0020] The teacher data creation unit 12 includes functional blocks for creating teacher data D3, such as a brain image acquisition unit 121, a region division unit 122, an image correction unit 123, a region of interest setting unit 124, a volume calculation unit 125, a t-value and p-value calculation unit 126, a z-value calculation unit 127, and a diagnosis result acquisition unit 128.

[0021] The brain image acquisition unit 121 acquires a brain image D1 from the diagnosis information database DB. The region division unit 122 to the z-value calculation unit 127 perform calculation processes on the acquired brain image D1, such as setting multiple regions of interest (ROI) in the brain region and calculating the z-value of each ROI. The specific calculation process contents of each unit of the region division unit 122 to the z-value calculation unit 127 will be described later.

[0022] The diagnostic result acquisition unit 128 acquires the diagnostic results D2 for each person whose brain image D1 was acquired from the diagnostic information database DB. The teacher data creation unit 12 creates teacher data D3 for each person by associating the z-value of each region of interest with the diagnostic results D2, and stores the teacher data D3 in the auxiliary storage device 11.

[0023] The learning unit 13 includes a first learning unit 131 and a second learning unit 132. The first learning unit 131 learns a prediction algorithm D4 based on teacher data D3 and stores the learned prediction algorithm D4 in the auxiliary storage device 11. The second learning unit 132 further learns the prediction algorithm D4 and stores the learned prediction algorithm D5 in the auxiliary storage device 11. The machine learning method is not particularly limited, but in this embodiment, the first learning unit 131 and the second learning unit 132 are configured as support vector machines.

[0024] (machine learning methods) The machine learning method according to this embodiment is implemented using the machine learning device 1 shown in Fig. 2. Fig. 3 is a flowchart showing the overall procedure of the machine learning method according to this embodiment. Fig. 4 is a flowchart showing the procedure of the teacher data creation step in the machine learning method according to this embodiment.

[0025] 3, the brain image acquiring unit 121 and the diagnosis result acquiring unit 128 acquire brain images D1 and diagnosis results D2 of multiple people, respectively, from the diagnostic information database DB. Note that the brain images D1 and diagnosis results D2 of one person may be acquired, or the brain images D1 and diagnosis results D2 of multiple people may be acquired at once.

[0026] In step S2, the teacher data creation unit 12 creates teacher data D3 from the acquired brain image D1 and the diagnosis result D2.

[0027] 4 is a flowchart showing a specific processing procedure of step S2 for creating training data. Step S2 mainly includes steps S21 to S27.

[0028] In step S21, the region segmentation unit 122 separates and removes non-brain tissue from the acquired brain image D1. The brain image from which non-brain tissue has been separated and removed is then divided into gray matter, white matter, and cerebrospinal fluid regions, and the lateral ventricles are separated from the cerebrospinal fluid region. In this embodiment, the region segmentation unit 122 uses signal intensity-dependent maximum likelihood and posterior probability methods to segment the brain image, preventing brain lesions from being ignored by conventional standardization methods such as SPM. To prevent white matter lesions from being mixed into the gray matter, a multi-channel segmentation technique is implemented that incorporates FLAIR images into the segmentation.

[0029] Specifically, as shown in Figure 5, FLAIR images with low spatial information are supplemented with 3D brain image data to provide high spatial information, and then only white matter lesions are extracted and filled in (substituted) with the average signal value of the subject's white matter. As a result, as shown in Figure 6, this embodiment enables separation with unprecedented precision.

[0030] 7 to 9 show other examples demonstrating this effect. In Fig. 7 to Fig. 9, with the conventional method, white areas are mixed in at the top two locations of the gray matter, and white matter is also missing at the top two locations. By using the method of this embodiment, separation with unprecedented accuracy is possible.

[0031] Thereafter, if necessary, the image quality of the brain image may be evaluated, and if the image quality is below a certain level, processing such as displaying a warning may be performed.

[0032] In step S22, the image correcting unit 123 nonlinearly converts the brain divided in step S21 into coordinates in the MNI space. In the conversion, the image correcting unit 123 converts the tensor quantity for each voxel into a signal value.

[0033] In step S23, the region-of-interest setting unit 124 sets multiple regions of interest in brain regions included in the brain image, i.e., in each region of the gray matter, white matter, and lateral ventricles. In this embodiment, the region-of-interest setting unit 124 divides the gray matter into 290 regions, the white matter into 33 regions, and the lateral ventricles separated from the other ventricles into two regions (the left ventricle and the right ventricle), and sets each divided region as a region of interest.

[0034] As described above, the region segmentation unit 122 separates the lateral ventricles from the cerebrospinal fluid portion. Normal cerebral atrophy causes the brain surface to shrink toward the center (creating a gap between the skull and the brain surface), but when a white matter lesion is present, the lateral ventricles expand compensatory and shrink from the inside to the outside. Because of this, in this embodiment, the region segmentation unit 122 separates the lateral ventricles. This process allows the boundaries between the lateral ventricles and the gray matter and white matter to be determined with high accuracy, thereby improving the accuracy of discrimination.

[0035] The effect of lateral ventricle separation in this embodiment is shown in Figure 10. Comparing the three examples obtained using the conventional method with this embodiment shown in the lower right, the boundary between the lateral ventricle and the white matter can be obtained with high accuracy, as indicated by the broken line in the figure. This allows for more accurate identification than conventional methods.

[0036] The gray matter was analyzed using Automated Anatomical Labeling (AAL) at 108 locations, which was created by the inventors. The brain can be divided into eight regions, including the entorhinal cortex, which is related to Alzheimer's disease, 118 regions from Brodmann, and 56 regions from the Loni Probabilistic Brain Atlas 40 (LPBA40). The quality can be divided by custom-created regions of interest in MNI space.

[0037] In conventional methods, the size of the corpus callosum could only be evaluated by its cross-sectional area in the sagittal section. However, in this embodiment, in order to evaluate it in three dimensions, the boundary between the corpus callosum and the surrounding white matter is set and divided by adjusting the parameters of the surface tension and viscosity of the fluid.

[0038] Because the corpus callosum is seamlessly connected to the subcortical white matter, a special technique is required to create a region of interest. More specifically, in a 3D brain image, a virtual fluid is placed at the frontal and occipital locations of the corpus callosum, and its boundary is determined by simulating the three-dimensional spreading of the virtual fluid within the brain. A typical example is to imagine a droplet of water equivalent to cerebrospinal fluid (CSF), and determine the shape of the frontal and occipital sides of the corpus callosum from the shape of the droplet as it spreads freely based on its surface tension and viscosity. This then determines the shape of the gray and white matter in contact with the corpus callosum. This allows for a simple yet highly accurate method to clarify the boundary surface, which actually contains complex, intricate three-dimensional shapes.

[0039] Figure 11 shows an example in which the three-dimensional structure of the corpus callosum is determined and its boundaries are clearly defined using the method of this embodiment. This clarifies the boundaries between gray matter and white matter, improving discrimination accuracy. Furthermore, by excluding the corpus callosum that falls within the region of interest for the lateral ventricles, highly accurate measurement of the lateral ventricle volume is possible.

[0040] In conventional methods such as SPM, this process is performed using Bayesian estimation to calculate posterior probabilities, with intermediate values ​​between 0 and 1 considered as partial volumes. However, this approach ignores outliers due to lesions, which is not in line with the objectives of the present invention. In this embodiment, Bayesian estimation is used only in the early stages, such as affine transformation and skull stripping, and segmentation is performed using maximum likelihood estimation based on image signal values. White matter lesions are a problem in this process. However, as shown in Figure 5, this embodiment solves this problem by using a 3D brain image to complement the spatial information of the FLAIR image, extracting white matter lesions from the FLAIR image, which has excellent contrast, and pasting them onto the 3D image.

[0041] The lateral ventricles can be divided using a template that the inventors have prepared in advance in the MNI space.

[0042] In step S24, the volume calculation unit 125 calculates the volume of each region of interest. In this embodiment, the volume calculation unit 125 calculates the volume using the Jacobian matrix obtained during tensor transformation. The reason for calculating volume rather than concentration is that volume values ​​can be universally used. For example, even if the z-value after statistical processing is the same, the volume value may differ. In conventional methods, the volume is calculated for each voxel and then calculated as the sum of the voxels within the region of interest. In this embodiment, the volume is calculated for each region of interest as a single unit. Theoretically, both methods produce the same results, but in practice, the volume value for each voxel is susceptible to noise, so calculating the volume for each region of interest provides higher accuracy.

[0043] In step S25, the t-value and p-value calculation unit 126 calculates the t-value by replacing the t-distribution with a normal distribution. For this purpose, the IXI database, which is used as a standard control group, is used. The IXI database contains approximately 100 normal brains for each age group, so there is no problem with replacing the t-distribution with a normal distribution.

[0044] Specifically, if the values ​​to be examined in the population (volumes of each region of interest) are normally distributed (or can be assumed to be normally distributed if the number of subjects in the study is large), the t-value can be calculated using the following formula to determine whether there is a statistically significant difference in the mean values ​​between two groups (healthy individuals and Alzheimer's patients).

[0045]

number

[0046] The degree of freedom is n-1.

[0047] The p-value indicates the t-value obtained from the above formula, at what t-value boundary can a result be considered statistically significant. The z-value can be calculated by replacing the p-value on the T-distribution with the p-value on the Z-distribution.

[0048] In step S26, the z-value calculation unit 127 calculates the z-value for each region of interest based on the t-value and p-value for that region, thereby calculating the z-values ​​for multiple regions of interest from the brain image D1.

[0049] The z-value is a value for statistical verification calculated from the t-value and p-value. Specifically, the z-value is a value representing the standard deviation indicating where the volume of the region of interest in a patient's region of interest in a certain part of a healthy subject corresponds in this normal distribution after the distribution is calculated and fitted to a normal distribution. In the case of a normal distribution (mean = 0, standard deviation = 1), the standard deviation value is calculated as the z-value, but in this embodiment, since a t-test is performed, the value obtained is the t-value. If the population is normally distributed, this value will be approximately the same as the z-value. In this case, the z-value is the z-value obtained in the z-test, meaning the z-value representing the standard deviation.

[0050] Because the z-score is a standardized value, it is suitable as an input value for subsequent AI. This is because the weighting is not biased in the early stages when AI learning extracts the desired feature from the input value.

[0051] In step S27, training data D3 is created by associating data consisting of a plurality of regions of interest and z values ​​with the diagnosis results D2.

[0052] Steps S21 to S27 described above complete step S2 shown in Fig. 3. The created teacher data D3 is stored in the auxiliary storage device 11, and steps S1 and S2 are repeated until a sufficient amount of teacher data D3 is accumulated in the auxiliary storage device 11 (YES in step S3).

[0053] Next, in step S4, the first learning unit 131 learns a prediction algorithm D4 (SVMst) based on the teacher data D3 stored in the auxiliary storage device 11. In this embodiment, the learning unit 13 performs learning by a support vector machine (SVM) using an RBF (radial basis function) kernel. At this time, the leave-one-out method is used for cross-validation, and the hardware The optimal value of the hyperparameter is calculated, and a highly versatile discriminant boundary is calculated between a subject group of Alzheimer's disease patients, a subject group of subjects who have developed Alzheimer's disease, and a subject group of subjects who have not developed Alzheimer's disease. The trained prediction algorithm D4 is stored in the auxiliary storage device 11.

[0054] Next, in step S5, the second learning unit 132 further inputs the Mini-Mental State Examination (MMSE; a set of questions developed by Folstein et al. in the United States in 1975 for diagnosing dementia) score from the diagnostic information database DB into the prediction algorithm D4 and performs additional learning to generate a prediction algorithm D5 (SVMcog). The prediction algorithm D5 is stored in the auxiliary storage device 11.

[0055] Note that learning by the first learning unit 131 and learning by the second learning unit 132 may be performed in parallel. That is, the first learning unit 131 may learn the diagnosis result D2 as a correct answer label and the z-value of each region of interest as a diagnostic input variable, and the second learning data 132 may learn the diagnosis result D2 as a correct answer label and the z-value of each region of interest and the MMSE score as input variables.

[0056] (diagnostic support device) Below, a description will be given of an embodiment in which a disease is determined using the trained prediction algorithm D4.

[0057] FIG. 12 is a block diagram showing the functions of a diagnostic support device 2 according to this embodiment. Similar to the machine learning device 1 shown in FIG. 2, the diagnostic support device 2 can be configured, for example, as a general-purpose personal computer. That is, the diagnostic support device 2 includes, as its hardware configuration, a CPU (not shown), a main memory device (not shown), and an auxiliary memory device 21. In the diagnostic support device 2, the CPU reads various programs stored in the auxiliary memory device 21 into the main memory device and executes them, thereby performing various arithmetic processing. The auxiliary memory device 21 can be configured, for example, as a hard disk drive (HDD) or a solid-state drive (SSD), and stores trained prediction algorithms D4 and D5. The auxiliary memory device 21 may be built into the diagnostic support device 2, or may be provided as an external memory device separate from the diagnostic support device 2.

[0058] The diagnostic support device 2 is connected to the MRI device 3, and brain images of the subject acquired by the MRI device 3 are transmitted to the diagnostic support device 2. Note that the brain images of the subject acquired by the MRI device 3 may be temporarily stored in a recording medium, and the brain images may be input to the diagnostic support device 2 via the recording medium.

[0059] The diagnostic support device 2 has a function of predicting the possibility that the subject will develop Alzheimer's disease within a predetermined period (e.g., within five years) (i.e., the possibility that the subject is in the ADNC spectrum) based on brain images of the subject. To realize this function, the diagnostic support device 2 includes an image processing unit 22 and a prediction unit 23 as functional blocks.

[0060] FIG. 13 shows the period from diagnosis to onset for 284 patients (pMCI) who were diagnosed with mild cognitive impairment in the ADNI database and subsequently developed Alzheimer's disease. This data reveals that 87.3% of pMCI patients developed the disease within three years, 95.8% within four years, and 97.5% within five years. Therefore, the predetermined period is not particularly limited, but is preferably 3 to 5 years.

[0061] The image processing unit 22 sets multiple regions of interest in the brain region for brain images acquired from an external source, performs arithmetic processing such as calculating z-values ​​for each region of interest, and outputs the z-values ​​for each region of interest to the prediction unit 23. To generate z-values ​​for each region of interest, the image processing unit 22 includes a brain image acquisition unit 221, a region division unit 222, an image correction unit 223, a region of interest setting unit 224, a volume calculation unit 225, a t-value and p-value calculation unit 226, and a z-value calculation unit 227. These functional blocks have the same functions as the brain image acquisition unit 121, the region division unit 122, the image correction unit 123, the region of interest setting unit 124, the volume calculation unit 125, the t-value and p-value calculation unit 126, and the z-value calculation unit 127 of the teacher data creation unit 12 shown in FIG. 2, respectively.

[0062] A brain image of the subject is acquired by brain image acquisition unit 221. Thereafter, each unit from region division unit 222 to z-value calculation unit 227 performs the processes of steps S21 to S27 shown in Fig. 4 to generate z-values ​​for each region of interest.

[0063] The prediction unit 23 predicts the possibility that the subject has the ADNC spectrum according to the prediction algorithm D4. In this embodiment, the prediction unit 23 predicts the possibility that the subject has the ADNC spectrum based on the z-value of each region of interest generated by the image processing unit 22. The prediction result is displayed, for example, on a display 4 connected to the diagnosis support device 2. The possibility of an ADNC spectrum can be calculated as a posterior probability (0 to 1) from the distance to a hyperplane (a hyperplane in elementary geometry, which generalizes a two-dimensional plane to other dimensions) using a sigmoid function. Alternatively, the diagnosis support device 2 may simply predict whether or not an object has the ADNC spectrum.

[0064] As described above, the diagnostic support device 2 uses the prediction algorithm D4 to predict the possibility that the subject is in the ADNC spectrum. Here, the prediction algorithm D4 is obtained by machine learning in the machine learning device 1, and by performing machine learning using a sufficient amount of training data D3, it is possible to improve the prediction accuracy of the diagnostic support device 2. In this way, in this embodiment, by using artificial intelligence, it is possible to predict with high accuracy the possibility that the subject is in the ADNC spectrum.

[0065] [Additional Notes] The present invention is not limited to the above-described embodiments, and various modifications are possible within the scope of the claims. Forms obtained by appropriately combining the technical means disclosed in the embodiments are also included in the technical scope of the present invention.

[0066] For example, in the above embodiment, MRI images are used as brain images, but X-ray CT images, SPECT images, PET images, etc. may also be used. Furthermore, changes over time in MRI images using tensor-based morphometry may also be used.

[0067] Furthermore, in the above embodiment, the machine learning device 1 includes both the teacher data creation unit 12 and the learning unit 13, but the teacher data creation unit 12 and the learning unit 13 may be configured to be implemented in separate devices. That is, teacher data D3 created in a device other than the machine learning device 1 may be input to the machine learning device 1, and the machine learning device 1 may only learn the prediction algorithms D4 and / or D5.

[0068] Similarly, the image processing unit 22 and the prediction unit 23 of the diagnostic support device 2 may be configured to be implemented by separate devices. In this case, the z-values ​​of each region of interest created in a device other than the diagnostic support device 2 may be input to the diagnostic support device 2, and the diagnostic support device 2 may perform only prediction based on the prediction algorithms D4 and / or D5.

[0069] In the above embodiment, the training data D3 for training the prediction algorithms D4 and D5 was generated from brain images of multiple ADNC patients and diagnostic results indicating whether each patient is on the ADNC spectrum. However, the present invention is not limited to this. For example, training data may be generated from brain images of AD patients, MCI patients, and healthy individuals without Alzheimer's disease. Even when a diagnostic support device is constructed using a prediction algorithm trained based on such training data, it is possible to predict the likelihood of a subject developing Alzheimer's disease with higher accuracy than conventional technology, as shown in the examples described below.

[0070] The technology described in Patent Document 2 differs from the present invention in that it predicts whether a patient with mild cognitive impairment (MCI) will develop Alzheimer's disease within a specified period. However, since there are no clear criteria for determining whether or not a patient has MCI, it is difficult to accurately select a prediction target. In contrast, the present invention does not clearly distinguish between healthy individuals and MCI, or between MCI and Alzheimer's disease, but instead targets the ADNC spectrum as a prediction target. In other words, as long as the patient is in the ADNC spectrum, it can be either an existing case of AD or a case that will develop AD in the future. This allows accurate assignment of training labels during learning, preventing a decrease in the accuracy of prediction results. Furthermore, the present invention predicts that a subject will progress to Alzheimer's disease in the future, even if they have not yet developed it, and therefore makes it possible to appropriately select ADNC patients who should be candidates for future disease-modifying therapy (DMT) using expensive drugs such as aducanumab.

[0071] In the above-described embodiment, the region dividing unit 122, 222 divides the brain image into gray matter, white matter, and cerebrospinal fluid regions and separates the lateral ventricles from the cerebrospinal fluid regions. Alternatively, the region of interest setting unit 124, 224 may set multiple regions of interest in the gray matter. The t-value and p-value calculation unit 126, 226 calculates the t-value and p-value for each region of interest for its volume. The z-value calculation unit 127, 227 calculates the z-value for each region of interest based on the t-value and p-value.

[0072] [Second embodiment] A second embodiment of the present invention will be described below with reference to the accompanying drawings. Note that the present invention is not limited to the following embodiment. Furthermore, components having the same functions as those in the first embodiment are given the same reference numerals, and their description will be omitted.

[0073] (Overall composition) FIG. 14 is a block diagram showing a schematic configuration of a prediction system 100′ according to this embodiment. The prediction system 100′ includes a machine learning device 1′ and a diagnostic support device 2′. The machine learning device 1′ learns a prediction algorithm for predicting the likelihood that an ADNC subject will develop Alzheimer's disease within a predetermined period of time. The diagnostic support device 2′ predicts the likelihood that an ADNC subject will develop Alzheimer's disease within a predetermined period of time according to the prediction algorithm learned by the machine learning device 1′. The machine learning device 1′ and the diagnostic support device 2′ may be implemented as separate devices, or the machine learning device 1′ and the diagnostic support device 2′ may be configured as a single device.

[0074] Below, configuration examples of the machine learning device 1 and the diagnosis support device 2 will be described.

[0075] (machine learning device) 15 is a block diagram showing the functions of a machine learning device 1' according to this embodiment. The hardware configuration of the machine learning device 1' may be the same as that of the machine learning device 1 shown in FIG.

[0076] The machine learning device 1' has a function of learning prediction algorithms D4' and D5' for predicting the likelihood that an ADNC subject will develop Alzheimer's disease within a predetermined period (for example, within five years).

[0077] To achieve this function, the machine learning device 1' includes functional blocks of a teacher data creation unit 12' and a learning unit 13. The teacher data creation unit 12' is a functional block that creates teacher data D3' from brain images D1' and diagnosis results D2' of multiple people. The multiple people are patients diagnosed with ADNC and healthy individuals.

[0078] The learning unit 13 is a functional block that learns prediction algorithms D4' and D5' based on the training data D3'. The training data creation unit 12' and the learning unit 13 are realized by executing a machine learning program stored in the auxiliary storage device 11.

[0079] The machine learning device 1' is capable of accessing a diagnostic information database DB. The diagnostic information database DB stores brain images D1' of multiple people and diagnostic results D2' indicating whether each person is in the ADNC spectrum and whether they are healthy. In this embodiment, the brain images D1' are three-dimensional MRI images.

[0080] The teacher data creation unit 12' includes functional blocks for creating teacher data D3', such as a brain image acquisition unit 121, a region separation unit 122', an image correction unit 123', a region of interest setting unit 124', a volume calculation unit 125', a covariate correction unit 126', a z-value calculation unit 127', and a diagnostic result acquisition unit 128.

[0081] The brain image acquisition unit 121 acquires a brain image D1 from the diagnosis information database DB. The region separation unit 122' to the z-value calculation unit 127' perform calculation processes on the acquired brain image D1, such as setting multiple regions of interest (ROI) in the brain region and calculating the z-value of each ROI. The specific calculation process contents of each unit of the region separation unit 122' to the z-value calculation unit 127' will be described later.

[0082] The diagnostic result acquisition unit 128 acquires the diagnostic results D2' for each person for whom the brain image D1' was acquired from the diagnostic information database DB. The teacher data creation unit 12 creates teacher data D3' for each person by associating the z-value of each region of interest with the diagnostic results D2, and stores the teacher data D3' in the auxiliary storage device 11.

[0083] The learning unit 13 includes a first learning unit 131 and a second learning unit 132. The first learning unit 131 learns a prediction algorithm D4' based on the teacher data D3' and stores the learned prediction algorithm D4' in the auxiliary storage device 11. The second learning unit 132 further learns the prediction algorithm D4 and stores the learned prediction algorithm D5' in the auxiliary storage device 11. The machine learning method is not particularly limited, but in this embodiment, as in the first embodiment, the first learning unit 131 and the second learning unit 132 are configured as support vector machines.

[0084] (machine learning methods) The machine learning method according to this embodiment is implemented using a machine learning device 1' shown in Fig. 15. Fig. 16 is a flowchart showing the overall procedure of the machine learning method according to this embodiment. Fig. 17 is a flowchart showing the procedure of the teacher data creation step in the machine learning method according to this embodiment.

[0085] 16, the brain image acquiring unit 121 and the diagnosis result acquiring unit 128 acquire brain images D1' and diagnosis results D2' of multiple people, respectively, from the diagnostic information database DB. Note that the brain images D1' and diagnosis results D2' for one person may be acquired, or the brain images D1' and diagnosis results D2' for multiple people may be acquired at once.

[0086] In step S2', the teacher data creation unit 12' creates teacher data D3' from the acquired brain image D1' and diagnosis result D2'.

[0087] 17 is a flowchart showing a specific processing procedure of step S2' for creating training data. Step S2' mainly includes steps S21' to S26'.

[0088] In step S21', the region separating unit 122' separates gray matter tissue from the brain image D1'. Specifically, the region separating unit 122' performs the processes of steps S211 to S213 shown in FIG.

[0089] In step S211, the brain image is divided into voxels and then registered. Specifically, in step S212 (described later), the brain image is subjected to four types of linear transformation (affine transformation): translation, rotation, scaling, and shear, to correct the spatial position and angle of the brain image in order to match its shape with a standard brain image for accurate gray matter separation. Specifically, 12 transformation parameters (translation, rotation, scaling, and shear) are calculated for each of the x, y, and z directions in three-dimensional space so as to minimize the sum of squares of the error between the brain image and the standard brain image template. Next, the brain image is affine-transformed using the calculated parameters, thereby achieving spatial registration of the brain image with the standard brain image, whose position, size, etc. are preset. Note that in this registration process, it is effective to use not only linear transformation but also nonlinear transformation to further approximate the shape to the standard brain image. This alignment transforms the brain into a standard brain image, which results in deformation of the cubic voxels that were previously divided. Therefore, after the alignment process, the brain image is again divided into voxels.

[0090] In step S212, gray matter extraction processing is performed using the newly divided voxels. The T1-weighted brain image contains three types of tissue: gray matter, which corresponds to neurons; white matter, which corresponds to lighter nerve fibers; and cerebrospinal fluid, which is nearly black. Therefore, the gray matter extraction processing focuses on gray matter tissue and extracts it voxel by voxel. Specifically, the brain image, which contains gray matter, white matter, and cerebrospinal fluid, is clustered into these three clusters to separate the gray matter. For this clustering processing, a density model and a model of the existence probability of the three tissues at different spatial locations can be used.

[0091] The density model models the fact that the distribution of voxel density values ​​differs depending on the tissue. If the tissues are arranged in descending order of density value (closest to white), the order is white matter, gray matter, and cerebrospinal fluid. Here, we assume that the density histograms after separating each tissue will be normally distributed.

[0092] The model of the probability of the existence of three tissues at a spatial location is a model that expresses the difference in spatial distribution of tissues due to individual differences in brain images in terms of probability. Here, we assume that each voxel belongs to one of the tissues and that the probability of the existence of each tissue at a spatial location is known in advance.

[0093] The optimal tissue distribution is estimated so that both of the above assumptions hold. The existence probability calculated for each voxel of gray matter, white matter, and cerebrospinal fluid tissue from brain images of many healthy individuals is used as a template to obtain brain images in which gray matter tissue is extracted in three dimensions. Separating gray matter tissue voxel-by-voxel based on existence probability results in small irregularities at the boundary surfaces, resulting in unnatural shapes.

[0094] Therefore, in step S213, image smoothing (first image smoothing) is performed on the gray matter brain image from which the gray matter tissue has been extracted. Here, the image is smoothed using a three-dimensional Gaussian kernel to improve the image's S / N ratio and to equalize its smoothness with that of the template image used in the subsequent second normalization. The FWHM (full width at half maximum) of the filter used for this smoothing is, for example, approximately 8 mm. Specifically, the three-dimensional brain image is three-dimensionally convoluted with a three-dimensional Gaussian function. This can be achieved by sequentially performing one-dimensional convolution in each of the x, y, and z directions.

[0095] Through the above processes of S211 to S213, the region separating unit 122' separates gray matter tissue from the brain image D1'.

[0096] Next, in step S22', the image corrector 123' transforms the gray matter brain image separated and smoothed in step S21' to a standardized template (hereinafter referred to as the "template") so that the image can be accurately segmented into regions of interest in the subsequent processing steps. Specifically, the image corrector 123' performs steps S221 to S223 shown in FIG. 18(b).

[0097] In step S221, a transformation called anatomical standardization is applied to the smoothed gray matter brain image, which involves global correction of the overall brain size and local correction of the local size to absorb individual differences in the anatomical shape and size of the gray matter brain image.

[0098] Specifically, image processing is performed using linear and nonlinear transformations to minimize the sum of squares of the error between the smoothed gray matter brain image and the template. The gray matter brain image template used here is an average image obtained from brain images of gray matter tissue extracted from many healthy individuals. This anatomical standardization process first performs global correction of position, size, and angle using linear transformation, and then performs local shape correction such as unevenness using nonlinear transformation. The linear transformation performed here is an affine transformation similar to that performed in step S211. Furthermore, the nonlinear transformation transforms the original image using nonlinear transformations using DARTEL in both the x and y directions.

[0099] In step S222, the gray matter brain image (hereinafter also referred to as the standardized brain image) deformed by anatomical standardization is divided into voxels again, and then image smoothing is performed again (second image smoothing). This process improves the S / N ratio of the standardized brain image and equalizes the image smoothness with that of a group of images of healthy subjects used as a standard for later comparison. The FWHM (full width at half maximum) of the filter is, for example, approximately 12 to 15 mm. Specifically, this can be achieved by performing the same process as the first image smoothing process in step S213, except for the FWHM value. By performing image smoothing again in this way, it is possible to reduce individual differences that do not result in complete agreement with the anatomical standardization process.

[0100] In step S223, density correction is performed on the standardized brain image after the second image smoothing. This involves correcting the voxel density values, which correspond to pixel values ​​in voxels. This is done by adding or subtracting a fixed value to the voxel values ​​of the standardized brain image so that the average voxel value of the standardized brain image matches the average voxel value of the gray matter brain image template.

[0101] Through the above processes of S221 to S223, the image corrector 123' corrects the gray matter brain image in accordance with the template.

[0102] In step S23', the region of interest setting unit 124' sets N regions of interest (ROIs) in the gray matter separated by the region separation unit 123'. In this embodiment, 290 regions of interest (N=290) are set by dividing the gray matter based on four types of atlases. The four types of atlases are Automated Anatomical Labeling (AAL) 108 locations, These included 8 areas such as the entorhinal cortex, which is related to Alzheimer's disease, 118 Brodmann areas, and 56 Loni Probabilistic Brain Atlas 40 (LPBA40). These 290 areas were selected as areas of interest. It is defined as an area.

[0103] In step S24', the volume calculation unit 125' calculates the volume X for each region of interest. In this embodiment, the gray matter separated from the brain image is partially or entirely compressed or expanded by anatomical standardization, with the compressed portions of the image appearing white and the expanded portions appearing black. The volume calculation unit 125' corrects the volume of each region of interest in the anatomically standardized gray matter based on the density of the image, thereby calculating the original volume in the space corresponding to each region of interest before standardization.

[0104] However, it has been confirmed that the volume of each region of interest calculated by the volume calculation unit 125' is biased due to age and intracranial volume. Therefore, in step S25', the covariate correction unit 126' performs a correction calculation to eliminate the influence of these biases, and the value obtained by this covariate correction is used as the X value, allowing the atrophy state of gray matter tissue in each region of interest to be evaluated under the same conditions.

[0105] Naturally, the gray matter tissue in the training data exhibits a bias not found in the gray matter tissue of healthy individuals, and this bias is a feature of Alzheimer's disease. Therefore, in this embodiment, the X values ​​of the region of interest in the gray matter tissue of healthy individuals are used as reference data, and the X values ​​of the region of interest in the gray matter tissue of the training data are used as comparison data, thereby enabling the feature of the training data to be obtained. It is known that the distribution of the X values ​​of the region of interest in the gray matter tissue of healthy individuals will be a normal distribution if the sample volume is sufficiently large, and this normal distribution can be determined by the mean value μ and standard deviation σ.

[0106] Therefore, in this embodiment, before performing machine learning using training data, the mean value μ and standard deviation σ are calculated for each region of interest using brain image data of healthy individuals obtained in advance from the IXI database, using the processing steps described above, in order to identify the normal distribution state of the X-values ​​for each region of interest. By identifying these 290 pairs of mean values ​​μ and standard deviations σ, it becomes possible to convert X-values ​​to z-values ​​in machine learning and diagnostic support processing. The calculated mean values ​​μ and standard deviation values ​​σ may be stored in the auxiliary storage device 11, etc.

[0107] Furthermore, the z-value calculation unit 127' calculates a z-value from the X-value for each region of interest in the brain image of the ADNC patient based on the average value μ and standard deviation σ for each region of interest. Specifically, the z-value is calculated by substituting the X-value, average value μ, and standard deviation σ into the following formula: z=(X-μ) / σ

[0108] In step S26', data consisting of the region of interest and z-values ​​in the brain image of the ADNC patient is associated with the diagnosis result D2 to create training data D3'.

[0109] The above steps S21' to S26' complete S2' shown in Figure 16. The created teacher data D3' is stored in the auxiliary storage device 11, and steps S1' and S2' are repeated until a sufficient amount of teacher data D3' is accumulated in the auxiliary storage device 11 (YES in step S3).

[0110] Next, in step S4, the first learning unit 131 learns a prediction algorithm D4' (SVMst) based on the training data D3' stored in the auxiliary storage device 11, and in step S5, the second learning unit 132 inputs the Mini-Mental State Test score into the prediction algorithm D4' and performs additional learning, thereby generating a prediction algorithm D5' (SVMcog).

[0111] (diagnostic support device) Below, a description will be given of an embodiment in which disease is determined using the trained prediction algorithm D4'.

[0112] 19 is a block diagram showing the functions of a diagnosis support device 2' according to this embodiment. The hardware configuration of the diagnosis support device 2' may be the same as that of the diagnosis support device 2 shown in FIG.

[0113] The diagnostic support device 2' has a function of predicting the possibility that the subject will develop Alzheimer's disease within a predetermined period (e.g., within five years) based on the brain image of the subject. To realize this function, the diagnostic support device 2' includes an image processing unit 22' and a prediction unit 23 as functional blocks.

[0114] The image processing unit 22′ separates gray matter from an externally acquired brain image, sets multiple regions of interest in the gray matter, performs computational processing such as calculating z-values ​​for each region of interest, and outputs the z-values ​​for each region of interest to the prediction unit 23. To generate z-values ​​for each region of interest, the image processing unit 22′ includes a brain image acquisition unit 221, a region separation unit 222′, an image correction unit 223′, a region of interest setting unit 224′, a volume calculation unit 225′, a covariate correction unit 226′, and a z-value calculation unit 227′. These functional blocks have the same functions as the brain image acquisition unit 121, the region separation unit 123′, the region of interest setting unit 124′, the volume calculation unit 125′, the covariate correction unit 126′, and the z-value calculation unit 127′ of the training data creation unit 12′ shown in FIG. 15 .

[0115] A brain image of the subject is acquired by the brain image acquisition unit 221. Thereafter, each unit from the region separation unit 223′ to the z-value calculation unit 227′ performs the processes of steps S21′ to S26′ shown in FIG. 17 to generate z-values ​​for each region of interest in the gray matter.

[0116] The prediction unit 23 predicts the possibility that the subject has the ADNC spectrum according to the prediction algorithm D4'. In this embodiment, the prediction unit 23 predicts the possibility that the subject has the ADNC spectrum based on the z-value of each region of interest generated by the image processing unit 22'. The prediction result is displayed, for example, on a display 4 connected to the diagnosis support device 2'.

[0117] After a diagnosis result is obtained for a subject, training data may be created by associating the subject's brain image data with the diagnosis result, and the prediction algorithm may be retrained using the training data. This allows the prediction accuracy of the prediction algorithm to improve over time. [Example]

[0118] Examples of the present invention will be described below, but the present invention is not limited to the following examples.

[0119] Example 1 In Example 1, the North American ADNI database (NA-ADNI) was used as the diagnostic information database DB shown in Figure 2. The inventors extracted 1,314 cases for which MRI brain image data existed from NA-ADNI. The breakdown was 359 AD patients, 412 MCI patients, and 543 healthy subjects (NL). Of the MCI patients, 284 were pMCI patients who progressed to AD during follow-up, and 128 were sMCI patients who had been followed for four years or more without progressing to AD.

[0120] 645 cases were randomly extracted from the above 1,314 cases and used as training data, and prediction algorithm D4 (SVMst) in Figure 2 was trained based on training data D3. Similarly, MMSE scores were input into prediction algorithm D4 to generate prediction algorithm D5 (SVMcog) in Figure 2.

[0121] Using SVMst and SVMcog, the possibility that a subject is on the ADNC spectrum was predicted using the above-mentioned training data as evaluation data. Specifically, the accuracy of prediction was measured using the accuracy rate (Accuracy), sensitivity, specificity, positive predictive value (PPV), negative predictive value (NPV), F1 value, Matthews correlation coefficient (MCC), relative risk, post-diagnosis odds (Odds), and AUC (area under the ROC curve; AUC ranges from 0 to 1, and a value of 0 is 0). The closer to 1, the higher the discrimination ability. When the discrimination ability is random, AUC=0.5.) was calculated. The results are shown in Table 1.

[0122] [Table 1]

[0123] From the above results, it was found that by using the prediction algorithms D4 and D5 in Example 1, it is possible to predict with high accuracy the possibility of an ADNC spectrum.

[0124] Example 2 In Example 2, the prediction algorithm D4 (SVMst) and prediction algorithm D5 (SVMcog) generated in Example 1 were evaluated for overlearning. Specifically, 669 cases, excluding the 645 cases extracted in Example 1 from the above 1,314 NA-ADNI cases, were used as evaluation data to calculate the prediction accuracy of the possibility that a subject was in the ADNC spectrum. As a comparative example, the prediction algorithm (VSRAD) disclosed in Patent Document 1 was also prepared, and the prediction accuracy of VSRAD was calculated in the same manner. The results are shown in Table 2.

[0125] [Table 2]

[0126] These results prove that the trained prediction algorithms D4 and D5 are not overtrained and have higher prediction accuracy than VSRAD.

[0127] Example 3 In Example 3, the prediction accuracy of the possibility of AD was evaluated for the prediction algorithm D4 (SVMst) generated in Example 1. Specifically, data extracted from the Japanese ADNI database (JADNI) and the Australian ADNI database (AIBL) in addition to NA-ADNI were used as evaluation data. As a comparative example, VSRAD disclosed in Patent Document 1 was also prepared, and the same data as the training data in Example 1 was analyzed by VSRAD to calculate the prediction accuracy of the trained VSRAD. The results, including a breakdown of AD and NL in NA-ADNI, JADNI, and AIBL, are shown in Table 3.

[0128] [Table 3]

[0129] These results show that the trained prediction algorithm D4 showed higher prediction accuracy than conventional prediction algorithms, even when predicting the possibility of AD. Furthermore, prediction algorithm D4 showed similar prediction accuracy when evaluated using data extracted from multiple databases other than the one used as training data, proving its high versatility.

[0130] Example 4 In Example 4, the prediction accuracy of the possibility of developing AD within a specified period at the pre-symptomatic stage of the ADNC spectrum was examined for the prediction algorithm D4 (SVMst) and prediction algorithm D5 (SVMcog) generated in Example 1, and the prediction algorithm (VSRAD) disclosed in Patent Document 1. The results are shown in Table 4. As a result, if the VSRAD test was positive, the relative risk of developing AD in the future was 1.9 times higher than if the VSRAD test was negative, whereas the relative risks for SVMst and SVMcog were 3.5 times and 3.6 times higher, respectively. This indicates that if the ADNC spectrum was predicted by the prediction algorithms D4 and D5, the risk of developing AD was 3.5 to 3.6 times higher than if the ADNC spectrum was not predicted.

[0131] [Table 4]

[0132] From the above results, it was found that by using the diagnostic support device (diagnostic support method) in Example 2, it is possible to predict with high accuracy the possibility that an ADNC patient will develop Alzheimer's disease within a predetermined period of time.

[0133] Example 5 In Example 5, the prediction accuracy of the possibility of AD by the prediction algorithm D4 (SVMst) generated in Example 1 was compared with that of two radiologists with over 20 years of experience. Specifically, a total of 200 cases (100 AD cases and 100 NL cases) were randomly selected from the NA-ADNI database. Furthermore, 10 MRI brain images of AD and 10 MRI brain images of NL were presented to the radiologists to learn how to diagnose AD and NL. After several days, the radiologists were asked to first diagnose whether the 200 cases, including the 20 cases already presented, as AD or NL, and then to re-diagnose them using the VSRAD results. Furthermore, using the 200 cases as evaluation data, prediction of AD or NL was performed using the prediction algorithm D4, and the prediction accuracy was calculated. Table 5 shows a comparison of the diagnostic accuracy by the two radiologists and the prediction accuracy by the prediction algorithm D4.

[0134] [Table 5]

[0135] From the above results, for example, the diagnostic accuracy rates of radiologists with VSRAD support were 70% and 73%, respectively, while the prediction accuracy rate of SVMst was 90.5%, meaning that the prediction accuracy of SVMst was clearly higher than the diagnostic accuracy of radiologists.

[0136] Example 6 In Example 6, we investigated whether the diagnosis support device (diagnosis support method) of the present invention can predict cerebral amyloid-β deposition. Regarding the criteria for the presence or absence of cerebral amyloid-β deposition, a cerebrospinal fluid amyloid-β level of 192 pg / ml or less was considered positive (present) in the NA-ADNI database. In the NA-ADNI database, 415 cases were diagnosed as ADNC spectrum using the prediction algorithm D4 (SVMst) in Example 1, and 90.6% (376 cases) of these were positive for cerebral amyloid-β deposition. From this, it can be inferred that the diagnosis support device (diagnosis support method) of the present invention accurately identifies the pathology of AD.

[0137] Example 7 In Example 7, the degree to which the onset of AD can be predicted within a certain period of time by the prediction algorithm D4 (SVMst) according to Example 1 is evaluated by plotting a progression-free survival curve. The study was conducted using a NA-ADNI database of pMCI and sMCI cases, including those in which biomarkers were measured in cerebrospinal fluid tests and those who underwent AV-45 amyloid PET testing. Figure 20 shows the relationship between the number of people (n) in each group, the number of years since the onset of AD, and the rate of developing AD. A(+) and pT(+) indicate that the cerebrospinal fluid biomarkers were positive for amyloid beta and phosphorylated protein, respectively.

[0138] Table 6 shows the hazard ratio and its confidence interval for each biomarker. tT(+) indicates that tau protein is a positive biomarker in the cerebrospinal fluid. ADNC refers to the A(+) and pT(+) group, and is presumed to have the pathological condition of AD. The hazard ratio for ADNC with A(+) and pT(+) is 2.18, whereas the hazard ratio when predicted as positive by the prediction algorithm D4 (SVMst) according to Example 1 is 3.59, indicating a higher risk of developing AD than ADNC patients.

[0139] [Table 6] [Explanation of symbols]

[0140] 1. Machine learning device 1' Machine learning device 11 Auxiliary storage 12 Training Data Creation Department 12' Training Data Creation Department 121 Brain Image Acquisition Unit 122 Area division part 122' Area separation part 123 Image Correction Unit 123' Image correction section 124 Region of interest setting section 124' Region of interest setting section 125 Volume calculation unit 125' Volume calculation section 126 t-value and p-value calculation section 126' Covariate correction section 127 z-value calculation unit 127' z-value calculation section 128 Diagnostic result acquisition unit 13 Learning Department 131 First Learning Section 132 Second Learning Section 2. Diagnostic support device 2' Diagnostic support device 21 Auxiliary storage 22 Image processing section 22' Image processing section 221 Brain Imaging Unit 222 Area division part 222' Area separation section 223 Image Correction Unit 223' Image correction unit 224 Region of interest setting section 224' Region of interest setting section 225 Volume calculation unit 225' Volume calculation section 226 t-value and p-value calculation section 226' Covariate correction section 227 z-value calculation unit 227' z-value calculation section 23 Prediction Department 3 MRI machine 4. Display D1 Brain imaging D1' brain image D2 Diagnosis result D2' diagnosis result D3 Training data D3' training data D4 Prediction Algorithm D4' prediction algorithm D5 Prediction Algorithm D5' prediction algorithm DB Diagnostic information database

Claims

1. A diagnostic support device that predicts the possibility that a subject having Alzheimer's disease neuropathological changes is either a patient who has already developed Alzheimer's disease or a patient with progressive mild dementia who will develop Alzheimer's disease within a predetermined period, comprising: a region segmentation unit that segments the brain image acquired from the subject into gray matter, white matter, and cerebrospinal fluid portions and separates the lateral ventricles from the cerebrospinal fluid portions; a region of interest setting unit that sets a plurality of regions of interest in each of the gray matter, the white matter, and the lateral ventricle; a t-value and p-value calculation unit that calculates t-values ​​and p-values ​​in each region of interest for the volume of each region of interest; a z-value calculation unit that calculates a z-value for each region of interest based on the t-value and the p-value; a prediction unit that performs the prediction according to a machine-learned prediction algorithm; Equipped with The prediction unit performs the prediction based on the z-value.

2. The diagnosis support device according to claim 1, The region segmentation unit separates the corpus callosum from the surrounding white matter by determining the boundary between the corpus callosum and the surrounding white matter using parameters of surface tension and viscosity of a fluid.

3. 3. The diagnosis support device according to claim 2, When a white matter lesion is present in the white matter, the region of interest setting unit extracts the white matter lesion and replaces it with an average signal value of the white matter of the subject, and then sets the region of interest in the white matter.

4. A diagnostic support device that predicts the likelihood that a subject with Alzheimer's disease neuropathological changes is either a patient who has already developed Alzheimer's disease or a patient with progressive mild dementia who will develop Alzheimer's disease within a predetermined period, comprising: a region separation unit that separates gray matter from the brain image acquired from the subject; a region of interest setting unit that sets a plurality of regions of interest in the gray matter; a volume calculation unit that calculates the volume of each region of interest; a z-value calculation unit that calculates a z-value of each region of interest based on the volume; a prediction unit that performs the prediction according to a machine-learned prediction algorithm; Equipped with The prediction unit performs the prediction based on the z-value.

5. The diagnosis support device according to claim 4, The diagnosis support device further comprises a covariate correction unit that performs covariate correction on the volume.

6. The diagnosis support device according to any one of claims 1 to 5, The prediction unit performs the prediction as a posterior probability from a distance to a hyperplane using a sigmoid function.

7. A machine learning device that learns a prediction algorithm for predicting the likelihood that a subject with Alzheimer's disease neuropathological changes is either a patient who has already developed Alzheimer's disease or a patient with progressive mild dementia who will develop Alzheimer's disease within a predetermined period, comprising: a training data creation unit that creates training data based on brain images of a plurality of people and a diagnosis result indicating whether each person has developed Alzheimer's disease before the predetermined period has elapsed since the acquisition of the brain images; and a learning unit that learns the prediction algorithm based on the training data; Equipped with The teacher data creation unit a region segmentation unit that segments a brain image acquired from the person into gray matter, white matter, and cerebrospinal fluid portions and separates a lateral ventricle from the cerebrospinal fluid portion; a region of interest setting unit that sets a plurality of regions of interest in each of the gray matter, the white matter, and the lateral ventricle; a t-value and p-value calculation unit that calculates t-values ​​and p-values ​​in each region of interest for the volume of each region of interest; a z-value calculation unit that calculates a z-value for each region of interest based on the t-value and the p-value; Equipped with A machine learning device, wherein the training data includes the diagnosis result and the z-score.

8. The machine learning device according to claim 7, The learning unit is a machine learning device configured with a support vector machine.

9. The machine learning device according to claim 7 or 8, The machine learning device, wherein the brain image is an MRI image.

10. A machine learning device that learns a prediction algorithm to predict the likelihood that a subject with Alzheimer's disease neuropathological changes is either a patient who has already developed Alzheimer's disease or a patient with progressive mild dementia who will develop Alzheimer's disease within a predetermined period of time, comprising: a training data creation unit that creates training data based on brain images of a plurality of people and a diagnosis result indicating whether each person has developed Alzheimer's disease before the predetermined period has elapsed since the acquisition of the brain images; and a learning unit that learns the prediction algorithm based on the training data; Equipped with The teacher data creation unit a region separation unit that separates gray matter from a brain image acquired from the person; a region of interest setting unit that sets a plurality of regions of interest in the gray matter; a volume calculation unit that calculates the volume of each region of interest; a z-value calculation unit that calculates a z-value of each region of interest based on the volume; Equipped with A machine learning device, wherein the training data includes the diagnosis result and the z-score.

11. The machine learning device according to claim 10, The machine learning device further comprises a covariate correction unit that performs covariate correction on the volume.

12. A diagnostic support method in which a computer predicts the likelihood that a subject having Alzheimer's disease neuropathological changes is either a patient who has already developed Alzheimer's disease or a patient with progressive mild dementia who will develop Alzheimer's disease within a predetermined period, the method comprising: a region segmentation step of segmenting the brain image acquired from the subject into gray matter, white matter, and cerebrospinal fluid portions and isolating the lateral ventricles from the cerebrospinal fluid portion; a region of interest setting step of setting a plurality of regions of interest in each of the gray matter, the white matter, and the lateral ventricle; a t-value and p-value calculation step of calculating t-values ​​and p-values ​​in each region of interest for the volume of each region of interest; a z-value calculation step of calculating a z-value for each region of interest based on the t-value and the p-value; a prediction step of making the prediction according to a machine-learned prediction algorithm; Equipped with A diagnostic support method, wherein in the predicting step, the prediction is made based on the z-value.

13. A machine learning method for learning a prediction algorithm for predicting the likelihood that a subject with Alzheimer's disease neuropathological changes is either a patient who has already developed Alzheimer's disease or a patient with progressive mild dementia who will develop Alzheimer's disease within a predetermined period, comprising: a training data creation step of creating training data based on brain images of a plurality of people and a diagnosis result indicating whether or not each person has developed Alzheimer's disease before the predetermined period has elapsed since the acquisition of the brain images; a learning step of learning the prediction algorithm based on the training data; Equipped with The teacher data creation step includes: a region segmentation step of segmenting the brain image acquired from the person into gray matter, white matter, and cerebrospinal fluid portions and isolating the lateral ventricles from the cerebrospinal fluid portions; a region of interest setting step of setting a plurality of regions of interest in each of the gray matter, the white matter, and the lateral ventricle; a t-value and p-value calculation step of calculating t-values ​​and p-values ​​in each region of interest for the volume of each region of interest; a z-value calculation step of calculating a z-value for each region of interest based on the t-value and the p-value; Equipped with A machine learning method, wherein the training data includes the diagnosis result and the z-score.

14. A machine learning program that trains a computer to learn a prediction algorithm for predicting the likelihood that a subject with Alzheimer's disease neuropathological changes is either a patient who has already developed Alzheimer's disease or a patient with progressive mild dementia who will develop Alzheimer's disease within a predetermined period of time, a training data creation step of creating training data based on brain images of a plurality of people and a diagnosis result indicating whether or not each person has developed Alzheimer's disease before the predetermined period has elapsed since the acquisition of the brain images; a learning step of learning the prediction algorithm based on the training data; causing the computer to execute The teacher data creation step includes: a region segmentation step of segmenting the brain image acquired from the person into gray matter, white matter, and cerebrospinal fluid portions and isolating the lateral ventricles from the cerebrospinal fluid portions; a region of interest setting step of setting a plurality of regions of interest in each of the gray matter, the white matter, and the lateral ventricle; a t-value and p-value calculation step of calculating t-values ​​and p-values ​​in each region of interest for the volume of each region of interest; a z-value calculation step of calculating a z-value for each region of interest based on the t-value and the p-value; Equipped with A machine learning program, wherein the training data includes the diagnostic results and the z-scores.

15. a training data creation step of creating training data from brain images of a plurality of people and diagnostic results indicating whether or not each person developed Alzheimer's disease before a predetermined period of time has elapsed since the acquisition of the brain images; a learning step of learning a prediction algorithm based on the training data; a prediction step of predicting, according to the prediction algorithm, the possibility that a subject having Alzheimer's disease neuropathological changes is either a patient who has already developed Alzheimer's disease or a patient with progressive mild dementia who will develop Alzheimer's disease within the predetermined period; An Alzheimer's disease prediction program that causes a computer to execute The teacher data creation step includes: isolating grey matter from brain images obtained from said person; establishing a plurality of regions of interest in the gray matter; calculating a volume for each region of interest; calculating a z-score for each region of interest based on the volume; creating the training data by associating the diagnosis results with the z-values; Equipped with The prediction step separating gray matter from brain images acquired from the subject; establishing a plurality of regions of interest in the gray matter; calculating a volume for each region of interest; calculating a z-score for each region of interest based on the volume; making the prediction based on the z-value; An Alzheimer's prediction program.

Images