JP7097854B2 - How to determine the risk of uterine fibroids - Google Patents

Description

The present invention relates to a method for determining the risk of uterine fibroids.

In order to use it for determining the risk of a disease, the association between a single nucleotide polymorphism (hereinafter, also referred to as "SNP") and a disease is being identified. NCBI SNP Database is a database that summarizes human SNPs, and manages SNPs with rs numbers. The rs number in the present specification shall also mean the registration number in this NCBI SNP Database.

The relationship between the SNP specified by the rs number in the present specification and those disclosed in non-patent documents and the like as diseases, pathological conditions, conditions, etc. related to the SNP is as follows.
rs5945326: SNP for type 2 diabetes (Non-Patent Document 1)
rs4432842: SNP regarding birth weight (Non-Patent Document 2)
rs4503880: SNP on liver enzyme level (γGTP) (Non-Patent Document 3)
rs2062377: SNP on bone density (spine) (Non-Patent Document 4)
rs2062377: SNP on bone density (base of femur) (Non-Patent Document 5)
rs12688220: SNP related to pancreatitis (Non-Patent Document 6)
rs651007: SNP regarding bad cholesterol level (Non-Patent Document 7)
rs5955543: SNP on nephroblastoma (Wilms tumor) (Non-Patent Document 8)
rs7896783: SNP regarding fibrinogen concentration in blood (Non-Patent Document 9)
rs1934179: SNP related to hypospadias (Non-Patent Document 10)

Hara K, Fujita H, Johnson TA, Yamauchi T, Yasuda K, Horikoshi M, et al. Genome-wide association study identifies three novel loci for type 2 diabetes. Hum. Mol. Genet. 2014; 23: 239-46. Horikoshi M, Yaghootkar H, Mook-Kanamori DO, Sovio U, Taal HR, Hennig BJ, et al. New loci associated with birth weight identify genetic links between intrauterine growth and adult height and metabolism. Nat. Genet. 2013; 45: 76 -82. Chambers JC, Zhang W, Sehmi J, Li X, Wass MN, Van der Harst P, et al. Genome-wide association study identifies loci influencing concentrations of liver enzymes in plasma. Nat. Genet. 2011; 43: 1131-8. Estrada K, Styrkarsdottir U, Evangelou E, Hsu YH, Duncan EL, Ntzani EE, et al. Genome-wide meta-analysis identifies 56 bone mineral density loci and reveals 14 loci associated with risk of fracture. Nat. Genet. 2012; 44 : 491-501. Estrada K, Styrkarsdottir U, Evangelou E, Hsu YH, Duncan EL, Ntzani EE, et al. Genome-wide meta-analysis identifies 56 bone mineral density loci and reveals 14 loci associated with risk of fracture. Nat. Genet. 2012; 44 : 491-501. Whitcomb DC, LaRusch J, Krasinskas AM, Klei L, Smith JP, Brand RE, et al. Common genetic variants in the CLDN2 and PRSS1-PRSS2 loci alter risk for alcohol-related and sporadic pancreatitis. Nat. Genet. 2012; 44: 1349-54. Kim YJ, Go MJ, Hu C, Hong CB, Kim YK, Lee JY, et al. Large-scale genome-wide association studies in East Asians identify new genetic loci influencing metabolic traits. Nat. Genet. 2011; 43: 990- Five. Turnbull C, Perdeaux ER, Pernet D, Naranjo A, Renwick A, Seal S, et al. A genome-wide association study identifies susceptibility loci for Wilms tumor. Nat. Genet. 2012; 44: 681-4. Sabater-Lleal M, Huang J, Chasman D, Naitza S, Dehghan A, Johnson AD, et al. Multiethnic meta-analysis of genome-wide association studies in> 100 000 subjects identifies 23 fibrinogen-associated Loci but no strong evidence of a causal association between circulating fibrinogen and cardiovascular disease. Circulation 2013; 128: 1310-24. van der Zanden LF, van Rooij IA, Feitz WF, Knight J, Donders AR, Renkema KY, et al. Common variants in DGKK are strongly associated with risk of hypospadias. Nat. Genet. 2011; 43: 48-50.

An object of the present invention is to provide a method for determining a risk of uterine fibroids (hereinafter, also referred to as "this disease").

The present inventors have diligently studied to solve the above problems. As a result, it was found that individual single nucleotide polymorphisms, which are seemingly not related to this disease at first glance, are related to this disease when viewed as a single set. Then, by using the relationship, the present invention for determining the risk of the present disease has been completed.

That is, in the method of the present invention, a single nucleotide polymorphism set containing at least rs5945326, rs4432842, rs4503880, rs2062377, rs12688220, rs651007, rs5955543, rs7896783, and rs1934179, which were found to be related to the present disease (hereinafter, "this SNP"). The risk of this disease is determined based on the genotype information of "set").

In the method of the present invention, the "single nucleotide polymorphism set" means one set of a plurality of single nucleotide polymorphisms, and the association with this disease has been found by this one set. ..

Further, the "genotype information" in the method of the present invention is a single nucleotide polymorphism, which is classified into two homozygous types (AA, BB) and a heterozygous type (AB) in the single nucleotide polymorphism. It means genotype information, and "genotype information of this SNP set" means a set of genotype information of each single nucleotide polymorphism specified in this SNP set, in other words. For example, it is a set of information about the genotypic base of each SNP in the base sequence indicated by each rs number. The genotype information of this SNP set is as shown in FIG.

According to the present invention, the risk of this disease can be determined.

The genotype information of this SNP set is shown. An example of a conversion table showing the relationship between the genotype information of this SNP set and the value associated with the junction type for each SNP is shown. The ROC curve and AUC of the model using this SNP set are shown. Further, an SNP set containing N-1 SNPs obtained by arbitrarily removing one SNP from this SNP set containing N SNPs is also referred to as a "comparison SNP set", and when representing each comparison SNP set, it is referred to as a "comparison SNP set". It is described as a comparative SNP set 1 and a comparative SNP set 2. The ROC curve and AUC of the model using the comparative SNP set 1 are shown. The ROC curve and AUC of the model using the comparative SNP set 2 are shown. The ROC curve and AUC of the model using the comparative SNP set 3 are shown. The ROC curve and AUC of the model using the comparative SNP set 4 are shown. The ROC curve and AUC of the model using the comparative SNP set 5 are shown. The ROC curve and AUC of the model using the comparative SNP set 6 are shown. The ROC curve and AUC of the model using the comparative SNP set 7 are shown. The ROC curve and AUC of the model using the comparative SNP set 8 are shown. The ROC curve and AUC of the model using the comparative SNP set 9 are shown. The ROC curve and AUC of the model using the comparative SNP set 10 are shown.

An embodiment of the present invention will be described. The following embodiments are examples for explaining the present invention, and the present invention is not intended to be limited to this embodiment. The present invention can be carried out in various forms as long as it does not deviate from the gist thereof.

In this embodiment, uterine fibroids are benign tumors that develop from the muscles of the uterus.

Further, in the present embodiment, the present disease is generally a disease diagnosed in accordance with the guidelines published by the Medical Society regarding the present disease, a disease described in the column of efficacy / effect in the attached document of a medical drug, and the like. Alternatively, it can be understood to mean at least one of the diseases understood as a general term in the pharmaceutical / medical industry.

In the method of the present embodiment, the risk of the disease is determined using a predetermined number of single nucleotide polymorphism sets that are seemingly unrelated to the disease.

The risk of this disease refers to the possibility of contracting this disease such as susceptibility to susceptibility to this disease and susceptibility to susceptibility to this disease. "Determining the risk" includes, for example, outputting the possibility of contracting the disease at present or in the future by dividing it into several levels, or outputting it numerically. Determining the risk of the disease includes assessing genetic factors or susceptibility to the disease, such as whether it is prone to disease or not.

In determining the risk of this disease, whether or not the subject who receives the determination of the risk of this disease is actually suffering (or developing) this disease at the time of determining the risk of this disease. Does not matter.

The method of the present embodiment is a set of genotypes in which the genotypes of each SNP specified in the SNP set are classified into two homozygous types (AA, BB) and heterozygous types (AB). Use the genotype information of the set. Then, based on the genotype information of this SNP set, the risk of this disease of the subject is determined.

The SNP set used in the method of the present embodiment is a set including SNPs that have not been previously recognized to be related to the disease. That is, usually, even if the SNPs included in this SNP set are analyzed individually, the risk of this disease cannot be determined. However, in the method of the present embodiment, the risk of the present disease can be determined by analyzing the genotype information of each SNP contained in the present SNP set as a set. Further, when the case where the present SNP set is analyzed and the case where the comparative SNP set is analyzed are compared, the statistically significant result is obtained in the case where the present SNP set is analyzed. That is, in the method of the present embodiment, by analyzing the SNP set to determine the risk of the disease, it is possible to provide a risk determination method having high accuracy or high predictive ability.

Hereinafter, in relation to each SNP included in this SNP set, the rs number, the chromosome number in which each SNP exists (in the case of a sex chromosome, indicated by X or Y), and the position on the chromosome of each SNP. , And the base sequence corresponding to the rs number are listed. In the base sequence indicated by each rs number, SNPs are enclosed in []. For example, when it is expressed as [A / G], it indicates that there is a single nucleotide polymorphism of A or G at the position of the base sequence. In addition, information such as a base sequence and a disease related to each SNP can be obtained, for example, by searching the NCBI SNP Database based on the rs number. Such information can be referred to by the Database and is incorporated herein by reference. The positions on the chromosome described below correspond to the version GRCh37 of the assembly genome.

rs5945326
Chromosome number X
Position on the chromosome 15289922
Nucleotide sequence ATTCCTTCTCTCCTGGGGGTTACAC [A / G] CATTGCTACGGGCCCATCCCAGAGC (SEQ ID NO: 1)

rs4432842
Chromosome number 5
Position on the chromosome 57172078
Nucleotide sequence CTCCCTCACTTAAAAAACTTACTT [T / C] TTCATGCCATTCAAACAGTCAACAT (SEQ ID NO: 2)

rs45038880
Chromosome number 18
Chromosome position 56084054
Nucleotide sequence GAGTGTGTCCAGTTCTTGTAAAATT [T / C] CAATCTTGTTATATGTTCAAAAGCCT (SEQ ID NO: 3)

rs2062377
Chromosome number 8
Position on the chromosome 120007420
Nucleotide sequence GGAAAGGACTACACTTTTGCTGCAA [A / T] TCAATAGTGAAGAAAATACATTTTG (SEQ ID NO: 4)

rs126888220
Chromosome number X
Chromosome position 106244767
Nucleotide sequence ATGTCCTTTGAGCATCATTTTTTAC [T / C] CCCATTGGGTGCTTTTACATTTTGTCT (SEQ ID NO: 5)

rs651007
Chromosome number 9
Position on the chromosome 136153875
Nucleotide sequence TGGCTGCTTTCAGGGAAGAAGTTCA [A / G] GGCAGTGGGCAGGGAGCTCTTGCACAC (SEQ ID NO: 6)

rs5955543
Chromosome number X
Position on the chromosome 176988397
Nucleotide sequence GGCTGATTACTCCCCATGGGAGGAAG [A / G] GGCTGGTTAGGAAGTGCATGGGCC (SEQ ID NO: 7)

rs7896783
Chromosome number 10
Position on the chromosome 65162153
Nucleotide sequence AAAATTGTTAGTTCCTTTGTGTTGA [A / G] TGAATAAAATTGAAAATTGTGACTGTGT (SEQ ID NO: 8)

rs1934179
Chromosome number X
Position on the chromosome 50182184
Nucleotide sequence TGTATTTCTTCCAATAGTGACTGGC [T / C] TTTAGGAGCCAATTGATAGAAAAAT (SEQ ID NO: 9)

In the method of the present embodiment, each SNP constituting the SNP set can be specified by referring to the base sequence specified by the rs number, but the rs number described in the present specification is merged with another rs number. And when a new rs number is assigned, the corresponding rs number in the present specification also means the rs number after merging and other rs numbers to be merged. Further, when the rs number described in the present specification is a number assigned by merging a plurality of rs numbers, the corresponding rs number in the present specification also means another original rs number.

Further, the above-mentioned base sequence indicated by each rs number regarding SNP is shown as a specific base sequence, but even if the base sequence in a portion other than the corresponding SNP in the base sequence is changed due to a difference in race or the like. good.

The method of the present embodiment can be used for subjects of any race, but can be particularly preferably used for Asians. Among Asians, it can be more preferably used by East Asian subjects such as Japanese. In addition, the method of this embodiment may be used for a subject of any gender.

Hereinafter, one aspect of a method for determining the risk of this disease by analyzing the genotype information of this SNP set will be described. However, the determination method is not limited to the following.

First, the subject's sample is used to identify the genotype of each SNP contained in this SNP set in the sample. The sample used for detecting SNP is not particularly limited as long as it is a sample containing chromosomal DNA. Examples of such a sample include body fluid samples such as saliva, blood, and urine; cell samples such as oral mucosa; and body hair such as hair. For the detection of SNP, the chromosomal DNA isolated by a conventional method from these samples may be directly used, or the isolated chromosomal DNA may be amplified and the amplified chromosomal DNA may be used.

Detection of SNPs can be performed by conventional gene polymorphism analysis methods. For example, DNA chip method (DNA microarray), sequence analysis using a conventional sequencer using the Sanger method, next-generation sequencer (NGS; Next Generation Sequencer), PCR (Polymerase Chain Reaction), hybridization, invader method, etc. However, it is not limited to these.

In the DNA chip method, a DNA chip in which a large number of DNA fragments (probes) including an SNP site are arranged on a substrate is used, chromosomal DNA is hybridized with the probe on the chip, and the binding site is detected by fluorescence or current. The sequence of chromosomal DNA is analyzed by. Examples of the DNA chip used for SNP analysis include a chip on which an oligonucleotide probe capable of detecting a base sequence containing an SNP site is arranged.

In addition, sequence analysis can be performed by a normal Sanger method. For example, a sequence reaction is performed using a primer set at a position of several tens of bases on the 5'side of a base showing a polymorphism, and the analysis result is used to determine what kind of base the corresponding position is. Can be done. Before the sequence reaction, it is preferable to amplify the fragment containing the SNP site by PCR or the like in advance. From the viewpoint of efficiency, NGS technology may be used.

Further, the detection of SNP can be performed, for example, by examining the presence or absence of amplification by conventional PCR. For example, a primer having a sequence corresponding to a region containing a base indicating a polymorphism and having a 3'end corresponding to each polymorphism is prepared. PCR can be performed using each primer to determine which type of polymorphism is present depending on the presence or absence of amplification products. In addition, the LAMP method (Loop-Mediated Isothermal Amplification; Patent No. 33133358), the NASBA method (Nucleic Acid Sequence-Based Amplification; Patent No. 2843586; Patent No. 2843586; Patent No. 2843586; Patent No. 2843586; Patent No. 2843586; It is also possible to check the presence or absence of amplification by referring to Japanese Patent No. 343929). In addition, a single-chain amplification method or an analysis method using NGS may be used.

It is also possible to amplify a DNA fragment containing an SNP site and determine which type of polymorphism is due to the difference in mobility of the amplified product in electrophoresis. Examples of such a method include the PCR-SCSP (single-strand conformation polymorphism) method (Genomics. 1992 Jan 1; 12 (1): 139-146.). Specifically, first, the DNA containing the target SNP is amplified, and the amplified DNA is dissociated into a single-stranded DNA. The dissociated single-stranded DNA can then be separated on a non-denatured gel and the type of polymorphism can be determined by the difference in mobility of the separated single-stranded DNA on the gel.

Furthermore, when a base showing a polymorphism is included in the restriction enzyme recognition sequence, it can be analyzed by the presence or absence of cleavage by the restriction enzyme (RFLP (Restriction Fragment Length Polymorphism) method). In this case, first, the DNA sample is cleaved with a restriction enzyme. The DNA fragment can then be separated and the size of the detected DNA fragment can be used to determine what type of polymorphism it is.

It is also possible to analyze the types of polymorphisms by examining the presence or absence of hybridization. That is, it is also possible to check which base the SNP is by preparing a probe corresponding to each base and checking which probe it hybridizes to.

In this way, the genotype data of the subject can be determined for each SNP of this SNP set. Here, the "genotype data of the subject" means the genotype information possessed by the subject.

Next, the risk of this disease is determined based on the genotype information of this SNP set. Any model can be used to determine the risk. The model is not particularly limited, but for example, a logistic regression model in which the feature amount calculated from the genotype data of the subject is input using the genotype information of this SNP set and the risk of this disease is output is used. Can be used. In the logistic regression model, parameters are machine-learned in advance using genotype data of humans suffering from this disease and genotype data of humans not suffering from this disease as training data.

In addition, as a model for determining the risk of a disease, instead of the logistic regression model, an arbitrary kernel function such as a multi-layer perceptron, a neural network such as CNN (Statistical Neural Network) and RNN (Recurrent Natural Network), or a Gaussian kernel is used. Various other models such as support vector machines, random forests modeled as regression trees, multiple regression analysis, models using hidden Markov models, statistical models and stochastic models can also be adopted. It is also possible to adopt a model that makes a comprehensive judgment by combining various models.

Next, an example of risk determination of this disease using a model will be described. First, the genotype data of the subject who determines the risk of this disease is converted into a feature amount that can be input to the model. The feature amount in the method of the present embodiment is, for example, whether the genotype data of the subject is homozygous (AA), homozygous (BB), or heterozygous (AB) for each SNP of this SNP set. It is a parameter indicating whether or not. The genotype is represented by nucleotides such as "GG" indicating that both SNPs of homologous chromosomes are G (guanine) and "AG" indicating that one is G (guanine) and the other is A (adenine). Therefore, the genotypic data of the subject is converted into parameters that can be input to the model using the genotypic information of this SNP set. However, if the model does not need to be converted to such parameters, then the above conversion is not required.

The conversion of the genotype data of the subject into the feature amount can be performed, for example, by adding a value to the genotype data of the subject for each SNP included in this SNP set. For example, for each SNP, the value (eg, for example) of the SNP depends on whether the subject's genotype data corresponds to homozygous (AA), homozygous (BB), or heterozygous (AB). , 0 or 1). This makes it possible to convert the genotype data of the subject into a feature amount. In the following, the case where the value corresponding to each SNP is 0 or 1 will be described as an example, but the value corresponding to the SNP is not limited to two values of 0 or 1.

The value associated with the junction type can be determined for each SNP. For example, one SNP associates a value 1 when the subject's genetic data is homozygous (AA) and associates a value 0 when it is homozygous (BB) and heterozygous (AB). In this way, the other SNPs associate a value of 1 when the subject's genotype data is heterozygous (AB), and a value of 0 when they are homozygous (AA) and homozygous (BB). May be associated with each other. In addition, if the genotype data of the subject is heterozygous (AB) and homozygous (BB), the value 1 is associated, and if it is homozygous (AA), the value 0 is associated. good.

As described above, the genotype data of the subject can be converted into a feature amount. The value used for the correspondence in the conversion to the feature amount can be arbitrarily determined. For example, based on the above non-patent document, a value 1 is associated with a genotype that is originally related to a disease to which each SNP is related, and a genotype that is not related to a disease to which each SNP is related is associated with a genotype. The value 0 can be associated.

Such a relationship between the junction type for each SNP and the value associated with the junction type can be expressed as, for example, a conversion table as shown in FIG. 2 based on the genotype information of the present SNP set as shown in FIG. can. In the conversion table of FIG. 2, when the genotype is shaded, the value associated with the SNP is set to 1, and when the genotype is not matched, the value associated with the SNP is set to 0. In the notation of specific genotypes in FIGS. 1 and 2, A indicates adenine, G indicates guanine, C indicates cytosine, and T indicates thymine. However, the format of the feature amount conversion table is not limited to FIG.

Finally, the risk of the disease in the subject is determined based on the genotype information of the SNP set. More specifically, using the conversion table based on the genotype information of this SNP set, the genotype data of the subject is calculated as a feature amount converted so that it can be input to the model, and the feature amount is calculated as a predetermined determination model. Can be entered in to determine the subject's risk of this disease.

In the judgment model, the feature amount is weighted to indicate that there is a positive correlation with the risk of this disease or a weight indicating that there is a negative correlation with the risk of this disease for each SNP of this SNP set. can do. For example, the values (features) associated with rs5945326, rs12688220, rs1934179, and rs5955543 are weighted to indicate that they have a positive correlation with the risk of this disease, and rs4432842, rs651007, rs7896783, rs4503880, and rs2062377. The associated value (feature amount) can be weighted to indicate that it has a negative correlation with the risk of this disease.

For example, when weighting the feature amount, this disease is when the genotype of rs5945326 is AG, the genotype of rs12688220 is TC, the genotype of rs1934179 is TC, and the genotype of rs5955543 is AG. The genotype of rs4432842 is TC, the genotype of rs651007 is AG, the genotype of rs7896783 is GG, the genotype of rs4503880 is TC or CC, and the genotype of rs2062377. When the genotype is AT, weighting can be given to indicate that there is a negative correlation with the risk of the disease. In addition, in the case of the genotype of each SNP associated with the value 0 as the feature amount, it can be evaluated that there is no correlation with the risk of this disease or it is low enough to be ignored.

Such weighting representing the correlation with the risk of this disease is a parameter using genotype data of humans suffering from this disease and genotype data of humans not suffering from this disease as learning data. Is identified by machine learning. At this time, if a certain model is weighted to indicate that a certain SNP has a positive correlation with the risk of this disease, the SNP has a positive correlation with the risk of this disease in other models as well. It is usual that weighting is given to indicate that there is. That is, it is difficult to imagine a situation in which the correlation with the risk of this disease is reversed in a certain SNP depending on the type of model or the like. The specific value of the weighting differs depending on the model and is not particularly limited.

Here, a group of SNPs that are weighted to indicate that there is a positive correlation with the risk of this disease in this SNP set is called a "positive correlation SNP set" and has a negative correlation with the risk of this disease. A group of SNPs that are weighted to represent is called a "negatively correlated SNP set". This SNP set includes a positively correlated SNP set and a negatively correlated SNP set, and based on such genotype information of this SNP set, the risk of this disease of the subject is increased. It is possible to comprehensively judge both the factors that reduce the risk and the factors that reduce the risk.

The determination result obtained as described above is also used as an aid when a specialist of this disease diagnoses this disease. In addition, the determination result of the risk of the present disease may be corrected based on the risk of the present disease determined as described above and the result of the questionnaire from the subject. In addition, advice on life improvement may be output to the subject based on the risk of this disease and the result of the questionnaire from the subject.

The present invention can also provide test reagents such as primers and probes. Examples of such a probe include a probe that includes the SNP site and can determine the type of base of the SNP site depending on the presence or absence of hybridization. Examples of the primer include a primer that can be used for PCR for amplifying the SNP site, and a primer that can be used for sequence analysis of the SNP site. In addition to these primers and probes, the test reagent of the present embodiment may contain a polymerase or buffer for PCR, a reagent for hybridization, and the like.

Hereinafter, this embodiment will be described in more detail with reference to Examples. However, this embodiment is not limited to these examples.

The relationship between this SNP set and this disease was verified as follows.

With the consent of more than 73,000 users of the gene analysis service, saliva samples and information on the morbidity of various diseases were collected. The morbidity information is, for example, a numerical value that becomes 1 when suffering from this disease and 0 when not suffering from this disease. The genotype data for each user was identified from the saliva sample, and a database was constructed in which the genotype data of the users and various morbidity information were associated with each other. From this database, a case-control set of 1500 subjects suffering from this disease and 1500 controls not suffering from this disease was constructed.

α＝０．１

The genotypes of each SNP in this SNP set of subjects and controls were then classified into two homozygous (AA, BB) and heterozygous (AB). Then, when the genotype matches the genotype of the shaded conversion table shown in FIG. 2, the value of x _i is set to 1, and when the genotype does not match, it is set to 0, and x ₁ to x _N are set to the following mathematical formulas. It was used as an explanatory variable of the logistic regression model represented by (1). For example, in the case of rs5945326, the value of x ₁ was set to 1 when the genotype was "AG", and the value of x ₁ was set to 0 when the genotype was "AA" or "GG". In this embodiment, N = 10. In addition, the objective variable of the logistic regression model represented by the following mathematical formula was a value p (morbidity information) between 0 and 1 representing the probability of suffering from this disease.

α = 0.1

1. 1. Model verification by AUC The accuracy of the judgment method using this SNP set will be described. From the above database, a data set in which the genotype information of the user and the morbidity information were associated with each other was created for the test. The genotypes of each SNP in this SNP set of each user in the data set are classified into homozygous (AA, BB) and heterozygous (AB), and the classified genotypes are shaded as shown in FIG. When the genotype was matched, the value of x _i was evaluated as 1, and when the genotype was not matched, the value was evaluated as 0, and x ₁ to x _N were calculated as feature quantities.

The feature amount of this SNP set for each user is input to the above logistic regression model (hereinafter, also referred to as "judgment model"), it is predicted whether or not each user has this disease, and the false positive is obtained. The rate and the true positive rate were calculated, and the ROC (Receiver Operating Characteristic) curve and the AUC (Area Under the Curve) were obtained, respectively. More specifically, a 5-fold cross-validation is performed on the judgment model, five ROC curves (ROC fold 1 to ROC fold 5) are obtained, and the average (Mean ROC) and standard deviation (± 1st d. Dev.) Are obtained. I asked. The broken line (Luck) in FIG. 3 is a case where whether or not the patient is suffering from this disease is randomly output, and corresponds to the ROC curve of the model having no predictive ability.

Further, in the same manner, a logistic regression model (hereinafter, also referred to as “comparison determination model”) was created for each comparative SNP set excluding one SNP from the present SNP set in the same manner as described above. Then, the feature amount related to each comparative SNP is input to each comparative judgment model, it is predicted whether or not each user has this disease, the false positive rate and the true positive rate are calculated, and the ROC curve and AUC are obtained. I asked for each. The results are shown in FIGS. 4 and later.

When the disease was determined using this SNP set, the AUC was 0.76 ± 0.02, which was significantly higher than that in the case of random output (AUC = 0.5), and the determination using this SNP set. It can be confirmed that the predictive power of the model is high.

On the other hand, in the case of the comparison determination model using each comparative SNP set, the AUC is lower than that in the case of using the present SNP set. Therefore, although the AUC of the comparison judgment model using each comparison SNP set is higher than that in the case of random output (AUC = 0.5), the AUC (0.76 ± 0.02) of the judgment model using this SNP set. It can be confirmed that it is generally lower than.

Therefore, by making a judgment using all the SNPs included in this SNP set, it is possible to determine whether or not the patient is suffering from this disease with higher accuracy than when using each comparative SNP set excluding one SNP from this SNP set. It turned out that it can be predicted with.

2. 2. Verification by Wilcoxon's Rank Sum Test To confirm that the judgment model using this SNP set is significantly superior to the comparative judgment model using each comparative SNP set, Wilcoxon's rank, which is a type of nonparametric test, is used. A sum test was performed. Specifically, we set a null hypothesis that there is no difference between the AUC of the judgment model using this SNP set and the AUC of the comparison judgment model using each comparative SNP set, and the significance level is 0.01, and Wilcoxon's ranking is set. A sum test was performed.

As a result, the p-values were all 3.96 × 10 ^-18 , and it was confirmed that the null hypothesis was rejected. That is, there is a statistically significant difference between the AUC of the judgment model using this SNP set and the AUC of the comparison judgment model using each comparison SNP set, and the judgment model using this SNP set uses each comparison SNP set. It can be said that the model is superior to the comparative judgment model used.

As described above, the method of the present embodiment has the effect that the accuracy of predicting whether or not the patient has the disease is significantly higher than the accuracy of the random prediction. Further, in the method of the present embodiment, there is a significant difference between the result of determination of this disease based on the genotype information of this SNP set and the result of determination of this disease based on the genotype information of the comparative SNP set. Has the effect of being. This effect is believed to be based on the potential correlation previously unseen between the genotype information of the SNP set and the disease. The logistic regression model and other models exemplified above are based on the genotype information of this SNP set, and learn data and disease information on the genotypes of humans suffering from this disease and those not suffering from this disease. It is obtained by machine learning the parameters by using as. That is, each model is only one phenotype representing the above-mentioned potential correlation, and the type of model used in the implementation of the method of the present embodiment is not particularly limited.

The method of the present invention determines the risk of the disease and contributes to its prevention and / or treatment in the fields related to medical care and health care.

Claims (2)

Using a machine-learned model using genotype data of humans with uterine fibroids and genotype data of humans without uterine fibroids as training data, there is a positive correlation with uterine fibroids rs5945326, Genotypic information for a single-base polymorphic set containing at least rs12688220, rs1934179, and rs5955543 and rs4432842, rs651007, rs7896783, rs4503880, and rs2062377 that are negatively correlated with uterine fibroids.
rs5945326 genotype is AG, rs12688220 genotype is TC, rs1934179 genotype is TC, and rs5955543 genotype is AG, rs4432842 genotype is TC, rs651007 genotype is AG, rs7896783 genotype is GG, A method for determining the risk of uterine myoma based on genotype information regarding whether the genotype of rs4503880 is TC or CC and the genotype of rs2062377 is AT. The method according to claim 1, wherein a body fluid sample, a cell sample, or hair of a subject who receives a risk determination is used.

Images