WO2021215532A1 - 月経前症候群の重症度検出方法 - Google Patents

月経前症候群の重症度検出方法 Download PDF

Info

Publication number
WO2021215532A1
WO2021215532A1 PCT/JP2021/016470 JP2021016470W WO2021215532A1 WO 2021215532 A1 WO2021215532 A1 WO 2021215532A1 JP 2021016470 W JP2021016470 W JP 2021016470W WO 2021215532 A1 WO2021215532 A1 WO 2021215532A1
Authority
WO
WIPO (PCT)
Prior art keywords
gene
pms
severity
genes
expression level
Prior art date
Application number
PCT/JP2021/016470
Other languages
English (en)
French (fr)
Inventor
夏海 長森
高良 井上
Original Assignee
花王株式会社
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 花王株式会社 filed Critical 花王株式会社
Priority to CN202180030487.6A priority Critical patent/CN115443337A/zh
Priority to US17/920,652 priority patent/US20230144067A1/en
Priority to EP21792472.9A priority patent/EP4141114A1/en
Publication of WO2021215532A1 publication Critical patent/WO2021215532A1/ja

Links

Classifications

    • CCHEMISTRY; METALLURGY
    • C12BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12QMEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q1/00Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
    • C12Q1/68Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
    • C12Q1/6876Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes
    • C12Q1/6883Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes for diseases caused by alterations of genetic material
    • CCHEMISTRY; METALLURGY
    • C12BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12QMEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q1/00Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
    • C12Q1/68Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
    • C12Q1/6869Methods for sequencing
    • CCHEMISTRY; METALLURGY
    • C12BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12QMEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q2600/00Oligonucleotides characterized by their use
    • C12Q2600/112Disease subtyping, staging or classification
    • CCHEMISTRY; METALLURGY
    • C12BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12QMEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q2600/00Oligonucleotides characterized by their use
    • C12Q2600/158Expression markers
    • CCHEMISTRY; METALLURGY
    • C12BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12QMEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
    • C12Q2600/00Oligonucleotides characterized by their use
    • C12Q2600/16Primer sets for multiplex assays

Definitions

  • the present invention relates to a method for detecting the severity of PMS.
  • PMS premenstrual syndrome
  • nucleic acids such as DNA and RNA in a biological sample
  • a comprehensive analysis method has been established and abundant information can be obtained by one analysis, and the function of the analysis result based on many research reports on monobasic polymorphism and RNA function. It has the advantage of being easy to link.
  • Nucleic acid derived from a living body can be extracted from tissues such as blood, body fluids, secretions and the like.
  • Patent Document 1 states that RNA contained in skin surface lipids (SSL) is used as a sample for analysis of a living body, and marker genes for epidermis, sweat glands, hair follicles, and sebaceous glands are detected in SL. Is described.
  • SSL skin surface lipids
  • Non-Patent Documents 1 to 3, Patent Document 2 There are genes whose expression changes with the production of the female hormone estrogen or progesterone.
  • AQP3 aquaporin3
  • HIF-1 ⁇ hyperoxia inducible factor 1 subunit ⁇
  • MFN2 mitofusin2
  • DEFB4A defensin beta 4A, or hBD-2
  • hBD-3 etc. have been reported so far.
  • Patent Document 1 International Publication No. 2018/008319
  • Patent Document 2 Japanese Patent Application Laid-Open No. 2015-228829
  • Hum Reprod, 2018, 33 (11): 2060-2073 Non-Patent Document 2
  • Biol Reprod, 2018, 99 (2): 308-318 Non-Patent Document 3
  • the present invention relates to the genes shown below: SH3BGRL3, NDUFA11, PSME2, SNORA24, SSR4, ZNF91, SNORA70, HSPA6, ALDH2, BRK1, C20orf111, CSNK1A1, STK10, GHITM, RPL21, MED13, RPL32, UBE2R2, RPS3, , ALDH1A3, TSPYL1, NUB1, KRT6A, CCL22, ELOVL3, TMEM66, GDI2, MYL12B, CCNI, PSMC6, RABGEF1, ARF1, ACTN1, CHMP2A, GAK, STX3, SERP1 and GAS7, and a group consisting of translations of the gene.
  • the present invention also presents nucleic acids derived from the following genes: SH3BGRL3, NDUFA11, PSME2, SNORA24, SSR4, ZNF91, SNORA70, HSPA6, ALDH2, BRK1, C20orf111, CSNK1A1, STK10, GHITM, RPL21, MED13, RPL32, UBE2R2, RPS PLEKHM2, DDX6, ALDH1A3, TSPYL1, NUB1, KRT6A, CCL22, ELOVL3, TMEM66, GDI2, MYL12B, CCNI, PSMC6, RABGEF1, ARF1, ACTN1, CHMP2A, GAK, STX3, SERP1 and GAS7, and translations of the genes.
  • markers of PMS severity including at least one selected from the group.
  • PMS Premenstrual Syndrome
  • luteal phase menstruation phase
  • Frank Archives of Neurology and Psychiatry, 1931, 26: 1053
  • PMS is defined as "physical and psychological symptoms that begin 3 to 10 days before the onset of menstruation and diminish or disappear with the onset of menstruation.”
  • Typical symptoms of PMS include lower abdominal pain, lower back pain, lower abdominal tension, headache, stiff shoulders, cold hands and feet, increased appetite, diarrhea, constipation, swelling, breast pain, and breast tension.
  • Prone to acne, rough skin, tiredness, sleepiness, etc., and mental symptoms include frustration, anger, aggressiveness, depression, fragile tears, increased anxiety, etc. ..
  • severe PMS means a state in which the physical or mental symptoms of PMS are severe
  • mild PMS means that the physical or mental symptoms of PMS are severe.
  • the severity of menstrual-related symptoms such as PMS has been evaluated based on, for example, the score of a questionnaire of Menstrual Distress Questionnaire (MDQ).
  • MDQ Menstrual Distress Questionnaire
  • PMS can be evaluated as mild.
  • the present invention relates to a PMS severity marker that enables detection of the severity of PMS, and a method for detecting the severity of PMS using the marker.
  • the present invention provides a method for detecting the severity of PMS using a PMS severity marker including a novel marker. Since the PMS severity marker used in the present invention can be collected from lipids on the surface of the skin (SSL), it can be obtained easily and non-invasively. Therefore, according to the present invention, it is possible to detect the severity of PMS without imposing a burden on the subject.
  • the present invention more simply and objectively detects the severity of PMS, which previously had to be evaluated based on the subject's interview during the luteal phase as described above, based on the expression of the gene or its translation product. Make it possible.
  • the present invention also makes it possible to pre-detect the severity of PMS based on the expression of a gene or its translation product in the follicular phase prior to the luteal phase in which PMS occurs.
  • the invention provides a PMS severity marker.
  • PMS severity marker As shown in the examples below, the present inventor has found that the expression of 51 genes shown in Table 1 below is related to the severity of PMS. Furthermore, 40 of these 51 genes shown in Table 2 were genes for which no association with PMS has been previously reported. Therefore, the nucleic acids derived from the genes shown in Table 1 and the translation products of the genes can be used as markers for the severity of PMS, respectively. In addition, nucleic acids derived from the genes shown in Table 2 and translation products of the genes are novel PMS severity markers.
  • nucleic acid derived from a gene includes mRNA transcribed from the gene, cDNA prepared from the mRNA, reaction product from the cDNA (for example, PCR product, cloned DNA, etc.), and the like. ..
  • the severity of PMS of a subject can be detected based on the expression level of each gene shown in Table 1 or its translation product.
  • the expression level of the gene shown in Table 1 or its translation product can be measured by quantifying the nucleic acid derived from the gene or the translation product of the gene. For example, by quantifying the mRNA transcribed from the gene, the cDNA prepared from the mRNA, the reaction product from the cDNA, or the translation product (protein, etc.) obtained by translating the mRNA, the gene or The expression level of the translation product can be measured.
  • Nucleic acid derived from the gene shown in Table 1 or a translation product of the gene can be prepared according to a conventional method from a biological sample collected from a subject, for example, cells, body fluids, secretions and the like.
  • the nucleic acid is an mRNA prepared from a subject's skin surface lipids (SSL).
  • the translation product is a translation product prepared from the subject's SSL.
  • “detection” of the severity of PMS can be paraphrased in terms such as prediction, examination, measurement, judgment or evaluation support.
  • the terms “detection”, “prediction”, “examination”, “measurement”, “judgment” or “evaluation” of the severity of PMS in the present invention refer to the diagnosis of the severity of PMS by a doctor. Does not include.
  • skin surface lipid refers to a fat-soluble fraction existing on the surface of the skin, and is sometimes called sebum.
  • SSL mainly contains secretions secreted from exocrine glands such as sebaceous glands in the skin, and is present on the surface of the skin in the form of a thin layer covering the skin surface.
  • SSL contains RNA expressed in skin cells (see Patent Document 1).
  • skin is a general term for regions including the epidermis, dermis, hair follicles, and tissues such as sweat glands, sebaceous glands, and other glands on the body surface, unless otherwise specified.
  • any means used to recover or remove SSL from the skin can be adopted.
  • an SSL-absorbing material, an SSL-adhesive material, or an instrument that scrapes the SSL from the skin, which will be described later, can be used.
  • the SSL-absorbing material or the SSL-adhesive material is not particularly limited as long as it is a material having an affinity for SSL, and examples thereof include polypropylene and pulp. More detailed examples of the procedure for collecting SSL from the skin include a method of absorbing SSL into a sheet-like material such as oil blotting paper and oil blotting film, a method of adhering SSL to a glass plate, tape, etc., a spatula, a scraper, etc.
  • the SSL-absorbent material is preferably used in a dry state.
  • the part of the skin from which the SSL is collected is not particularly limited, and examples thereof include the skin of any part of the body such as the head, face, neck, trunk, and limbs, and the part where sebum is secreted, for example, the skin of the face. Is preferable.
  • the SSL collected from the subject may be stored for a certain period of time.
  • the collected SSL is preferably stored under low temperature conditions as soon as possible after collection in order to suppress the decomposition of the contained RNA as much as possible.
  • the storage temperature condition of the RNA-containing SSL in the present invention may be 0 ° C. or lower, preferably ⁇ 20 ⁇ 20 ° C. to ⁇ 80 ⁇ 20 ° C., more preferably ⁇ 20 ⁇ 10 ° C. to ⁇ 80 ⁇ 10 ° C. , More preferably ⁇ 20 ⁇ 20 ° C. to ⁇ 40 ⁇ 20 ° C., further preferably ⁇ 20 ⁇ 10 ° C.
  • the storage period of the RNA-containing SSL under the low temperature condition is not particularly limited, but is preferably 12 months or less, for example, 6 hours or more and 12 months or less, more preferably 6 months or less, for example, 1 day or more and 6 months or less. More preferably, it is 3 months or less, for example, 3 days or more and 3 months or less.
  • RNA extraction or purification of RNA from a biological sample for example, a phenol / chloroform method, an AGPC (acid guanidinium thiocyanate-phenol-chromoform extension) method, or TRIzol ( Method using columns such as registered trademark), RNA®, QIAzol (registered trademark), method using special magnetic material particles coated with silica, method using Solid Phase Reversible Imaging magnetic material particles, ISOGEN, etc. Extraction with a commercially available RNA extraction reagent or the like can be used.
  • a method usually used for extracting or purifying a protein from a biological sample for example, a commercially available protein extraction reagent such as QIAzol Lysis Reagent (Qiagen) can be used.
  • a commercially available protein extraction reagent such as QIAzol Lysis Reagent (Qiagen) can be used.
  • the severity of PMS in a subject can be detected by examining the expression of one or more genes or their translation products shown in Table 1 using the nucleic acid or translation product prepared by the above procedure.
  • the invention comprises measuring the expression of at least one selected from the group consisting of the genes shown in Table 1 or their translations, the method of detecting the severity of PMS in a subject. I will provide a.
  • the subject in the method for detecting the severity of PMS according to the present invention is, for example, a female of a human or non-human mammal who needs to detect the severity of PMS. Just do it.
  • the subject is an animal having SSL on the skin, more preferably a human.
  • the method of the present invention comprises measuring the expression of at least one gene (hereinafter, also referred to as a target gene (group)) selected from the group consisting of the genes shown in Table 1 in a subject.
  • the method of the invention measures the expression of at least one (hereinafter, also referred to as a target product (group)) selected from the group consisting of translation products of the genes shown in Table 1 in a subject. Including that.
  • both the expression of the target gene (group) and the expression of the target product (group) may be measured.
  • the expression of the target gene (group) is measured.
  • the target gene (group) is at least one selected from the group consisting of the genes shown in Table 2, and the target product (group) is shown in Table 2. At least one selected from the group consisting of gene translations. Further, in the method of the present invention, at least one selected from the group consisting of the genes shown in Table 2 or translation products thereof and other genes (that is, CTNNB1, AQP3, ANXA1, ARHGDIB, DEFB4A, shown in Table 1).
  • the expression of the target gene (group) or target product (group) in the method of the present invention is measured using RNA or a translation product contained in the SSL of the subject.
  • the method of the invention may further comprise collecting the SSL of the subject.
  • the method of the invention may further comprise extracting RNA or translation product from SSL taken from the subject.
  • the expression level of the target gene can be measured by quantifying the amount of RNA derived from the target gene contained in the sample derived from the subject, for example, RNA extracted from SSL.
  • the amount of subject-derived RNA may be measured according to the procedure for gene expression analysis using RNA commonly used in the art.
  • RNA commonly used in the art.
  • the method of gene expression analysis using RNA after converting RNA into cDNA by reverse transcription, the cDNA or its amplification product is quantified by real-time PCR, multiplex PCR, microarray, sequencing, chromatography or the like. The method can be mentioned.
  • the expression level of the translation product can be measured by using protein quantification methods commonly used in the art, such as ELISA, immunostaining, fluorescence, electrophoresis, chromatography, mass analysis and the like.
  • the measured expression level may be the absolute amount of expression of the target gene (group) or target product (group) in the biological sample, or relative to the total expression level of other standard substances, genes or translation products. It may be the expression level.
  • the subject-derived RNA is converted to cDNA by reverse transcription, and then the cDNA is subjected to PCR to purify the resulting reaction product.
  • primers targeting the specific RNA to be analyzed may be used, but it is preferable to use random primers for more comprehensive nucleic acid storage and analysis.
  • a general reverse transcriptase or reverse transcriptase kit can be used for the reverse transcription.
  • a highly accurate and efficient reverse transcriptase or reverse transcriptase kit is used, and examples thereof include M-MLV Reverse Transcriptase and its variants, or a commercially available reverse transcriptase or reverse transcriptase kit.
  • PrimeScript registered trademark
  • Reverse Transcriptase series Takara Bio Co., Ltd.
  • SuperScript registered trademark
  • Reverse Transcriptase series Thermo Scientific
  • SuperScript registered trademark
  • III Reverse Transcriptase SuperScript (registered trademark) III Reverse Transcriptase
  • VILO cDNA Synthesis kit all of which are Thermo Scientific
  • PCR of the obtained cDNA only the specific DNA may be amplified using a primer pair targeting the specific DNA to be analyzed, or a plurality of DNAs may be amplified using a plurality of primer pairs. good.
  • the PCR is multiplex PCR.
  • Multiplex PCR is a method of simultaneously amplifying a plurality of gene regions by simultaneously using a plurality of primer pairs in a PCR reaction system. Multiplex PCR can be carried out using a commercially available kit (for example, Ion AmpliSeqTranscriptome Human Gene Expression Kit; Life Technologies Japan Co., Ltd., etc.).
  • the temperature is adjusted to preferably 42 ° C. ⁇ 1 ° C., more preferably 42 ° C. ⁇ 0.5 ° C., still more preferably 42 ° C. ⁇ 0.25 ° C.
  • the reaction time is preferably adjusted to 60 minutes or longer, more preferably 80 to 120 minutes.
  • the temperature of the annealing and extension reactions in the PCR can be appropriately adjusted depending on the primers used, but for example, when the above-mentioned multiplex PCR kit is used, it is preferably 62 ° C.
  • the annealing and extension reactions are preferably carried out in one step.
  • the time of the annealing and extension reaction steps can be adjusted depending on the size of the DNA to be amplified and the like, but is preferably 14 to 18 minutes.
  • the conditions of the denaturation reaction in the PCR can be adjusted depending on the DNA to be amplified, preferably 95 to 99 ° C. for 10 to 60 seconds. Reverse transcription and PCR at temperature and time as described above can be performed using a thermal cycler commonly used for PCR.
  • the purification of the reaction product obtained by the PCR is preferably performed by size separation of the reaction product.
  • size separation By size separation, the PCR reaction product of interest can be separated from the primers and other impurities contained in the PCR reaction solution.
  • DNA size separation can be performed by, for example, a size separation column, a size separation chip, magnetic beads that can be used for size separation, or the like.
  • Preferred examples of magnetic beads that can be used for size separation include Solid Phase Reversible Imaging (SPRI) magnetic beads such as Aple XP.
  • SPRI Solid Phase Reversible Imaging
  • the purified PCR reaction product may be subjected to further treatment necessary for subsequent quantitative analysis.
  • the purified PCR reaction product can be prepared into a suitable buffer solution, the PCR primer region contained in the PCR-amplified DNA can be cleaved, or the adapter sequence can be added to the amplified DNA. May be further added.
  • the purified PCR reaction product is prepared into a buffer solution, the PCR primer sequence is removed and adapter ligation is performed on the amplified DNA, and the obtained reaction product is amplified as necessary for quantitative analysis. Library can be prepared.
  • RNA expression can be quantified based on the number of reads created by sequencing (read count).
  • the method of the present invention detects the severity of PMS in a subject based on the expression level of the target gene (group) or the target product (group).
  • the genes shown in Table 1 can be divided into two groups shown in Tables 3 and 4 according to the pattern of fluctuations in their expression levels with the severity of PMS. These patterns of variation can be used to detect the severity of PMS in a subject.
  • the genes shown in Table 3 are genes whose expression increases when PMS is severe and decreases when PMS is mild.
  • the genes shown in Table 4 are genes whose expression decreases when PMS is severe and whose expression increases when PMS is mild.
  • the severity of PMS of a subject can be detected.
  • the underlined 22 and 18 genes, respectively are genes for which no association with PMS has been reported (genes shown in Table 2).
  • the expression level of each target gene (group) or target product (group) is directly used as a parameter for detecting the severity of PMS.
  • the expression level of a target gene or target product of a subject at a specific time of the menstrual cycle is measured regularly (for example, monthly), and fluctuations in the expression level are tracked. This makes it possible to detect the severity of PMS in the subject.
  • the subject by comparing the expression level of the target gene or target product measured from the subject at a specific time of the menstrual cycle with the reference expression level of the target gene or target product at the same time of the menstrual cycle, the subject.
  • the severity of PMS can be detected.
  • a predetermined value for example, a statistical value of the expression level of the target gene or target product at a specific time of the menstrual cycle (for example, average expression level) can be used. ..
  • the reference expression level may be set for each subject. For example, it may be calculated from a statistical value (for example, an average value) of the expression level of a target gene or a target product at a specific time of the menstrual cycle measured multiple times from the same subject.
  • the reference expression level may be a statistical value (for example, an average value) of the expression level of the target gene or target product at a specific time of the menstrual cycle measured from the population.
  • a plurality of genes or translation products are used as the target gene group or target product group, it is preferable to determine the reference expression level for each of the genes or translation products.
  • the target gene or target product is the gene shown in Table 3 or its translation product
  • the subject's PMS can be detected as severe, and the expression level is the reference expression level. If lower, the subject's PMS can be detected as non-severe (or mild).
  • the target gene or target product is the gene shown in Table 4 or its translation product
  • the expression level is lower than the reference expression level
  • the subject's PMS can be detected as severe, and the expression level is the reference expression level. If higher, the subject's PMS can be detected as non-severe (or mild).
  • two or more of the genes shown in Tables 3 to 4 or their translation products may be combined and used as a target gene or target product.
  • a certain percentage for example, 50% or more, preferably 70% or more, more preferably 90% or more, still more preferably 100% of the genes or translation products are the criteria for the above-mentioned expression level.
  • the severity of PMS in a subject can be detected based on whether or not the condition is met. For example, if the expression level of 50% or more of the target gene group or target product group selected from Table 3 is higher than the reference expression level, the PMS of the subject can be detected as severe.
  • the expression level of 50% or more of the target gene group or target product group selected from Table 4 is lower than the reference expression level, the PMS of the subject can be detected as severe.
  • the expression level of 50% or more of the target gene group or target product group selected from Table 3 is higher than the reference expression level, and the expression level of 50% or more of the target gene group or target product group selected from Table 4 is higher. If the expression level is lower than the reference expression level, the subject's PMS can be detected as severe.
  • the combination of the target gene group or the target product group in the method of the present invention is not limited to the above example.
  • the severity of PMS in a subject is detected based on a predictive model constructed using expression data of a target gene (group) or target product (group).
  • expression data data on the expression level of the target gene (group) or target product (group) at a specific time of the menstrual cycle (for example, follicular phase or luteal phase) can be used.
  • PMS is performed by machine learning using expression data for each of one or more target genes (groups) or target products (groups) as an explanatory variable and the severity of PMS (for example, severe or mild) as an objective variable. It is possible to build an optimal predictive model for detecting the severity of the disease. For example, from the genes shown in Table 1 or Table 2 or their translation products, a plurality of genes having a larger expression difference between the severe group and the mild group in PMS severity are selected, and each of the expression data is explained. It can be used as a variable.
  • the prediction model may be constructed after compressing the data by dimension reduction, if necessary. For example, a plurality of genes are selected from the genes shown in Table 1 or Table 2, and those genes or their translation products are extracted as a target gene group or a target product group. Next, a principal component analysis is performed on the expression level of the extracted target gene group or target product group. Optimal for detecting the severity of PMS by machine learning with one or more principal components calculated by the principal component analysis as explanatory variables and the severity of PMS (for example, severe or mild) as the objective variable. Predictive model can be constructed.
  • the target gene (group) or target product (group) used to construct the predictive model has a higher expression difference between the severe group and the mild group in the severity of PMS among the genes shown in Table 1 or their translation products. It is preferable to select one or more from the larger ones. For example, one or more genes or translations thereof selected from the 40 genes shown in Table 2, for which no association with PMS has been previously reported, are preferred target genes (groups) or targets used in the construction of predictive models. Listed as a product (group). The same applies to the target gene (group) or target product (group) used for principal component analysis.
  • genes having a larger expression difference between the severe PMS group and the mild PMS group in the luteal phase are selected.
  • About 5 to 20 genes, preferably about 5 to 15 genes are selected. More specifically, for example, from Table 2, No. 5 genes 1 to 5, No. 10 genes from 1 to 10 or No. 15 genes from 1 to 15, preferably No. 10 genes from 1 to 10 or No. Select 15 genes from 1 to 15.
  • a predictive model can be constructed using data on the expression level of the selected gene in the luteal phase as an explanatory variable.
  • the genes having a larger expression difference between the severe PMS group and the mild PMS group in the luteal phase are selected. More specifically, for example, from Table 2, No. 5 genes 1 to 5, No. 10 genes from 1 to 10 or No. 15 genes from 1 to 15, preferably No. Select 5 genes from 1 to 5.
  • a predictive model can be constructed using data on the expression level of the selected gene in the follicular phase as an explanatory variable.
  • the dimensions and number of principal components used to build the prediction model are not particularly limited, but components with a large contribution rate, for example, the first principal component is included and the total contribution rate is high.
  • the explanatory variable may be one or more principal components such that the content is preferably 20% or more, more preferably 50% or more.
  • principal component analysis is performed on the expression data of the target gene group in the luteal phase, and the first principal component is contained and the total contribution rate is 20% or more, preferably 50% or more, for example, preferably 22.
  • a prediction model can be constructed by selecting principal components so as to be ⁇ 55% or more, more preferably 50 to 70%, and even more preferably 55 to 65%, and using the selected principal components as explanatory variables.
  • principal component analysis is performed on the expression data of the target gene group in the follicular phase, and the first principal component is contained and the total contribution rate is 80% or more, preferably 85% or more, more preferably.
  • a known algorithm such as an algorithm used for machine learning can be used.
  • machine learning algorithms include algorithms such as Random Forest, Neural Net, Linear Kernel Support Vector Machine (SVM (linear)), and rbf Kernel Support Vector Machine (SVM (rbf)).
  • SVM Linear Kernel Support Vector Machine
  • rbf Kernel Support Vector Machine
  • the severity of PMS of the subject can be determined. Can be detected.
  • it comprises detecting the severity of PMS based on the expression level of the gene or translation product. More preferably, it comprises detecting the severity of PMS based on the expression level of the gene.
  • it comprises detecting the severity of PMS using a predictive model based on the expression level of the gene or translation product.
  • the prediction model uses the expression level of the gene or translation product in the severe and mild groups of PMS, or at least one principal component calculated by principal component analysis of the expression level as an explanatory variable, and premenstrual syndrome. Constructed with the severity of the syndrome as the objective variable, The method according to any one of [1] to [5].
  • Nucleic acids derived from the following genes are preferably: SH3BGRL3, NDUFA11, PSME2, SNORA24, SSR4, ZNF91, SNORA70, HSPA6, ALDH2, BRK1, C20orf111, CSNK1A1, STK10, GHITM, RPL21, MED13, RPL32, UBE2R2.
  • RPS3, PLEKHM2, DDX6, ALDH1A3, TSPYL1, NUB1, KRT6A, CCL22, ELOVL3, TMEM66, GDI2, MYL12B, CCNI, PSMC6, RABGEF1, ARF1, ACTN1, CHMP2A, GAK, STX3, SERP1 and GAS7 The use according to [8], wherein at least one selected from the group consisting of is used.
  • nucleic acids derived from the following genes selected from the group consisting of CTNNB1, AQP3, ANXA1, ARHGDIB, DEFB4A, BIRC3, S100A7, SFN, HMGN2, HIF1A, and MFN2, and translation products of the genes.
  • the use according to any one of [8] to [10], wherein the nucleic acid derived from the gene is the mRNA of the gene contained in the lipid on the surface of the skin.
  • a marker of the severity of PMS including at least one selected from.
  • nucleic acids derived from the following genes are further selected from the group consisting of CTNNB1, AQP3, ANXA1, ARHGDIB, DEFB4A, BIRC3, S100A7, SFN, HMGN2, HIF1A, and MFN2, and translation products of the genes.
  • the marker according to [12] which comprises at least one of the markers.
  • the marker according to [12] or [13], wherein the nucleic acid derived from the gene is the mRNA of the gene contained in the lipid on the surface of the skin.
  • Example 1 Identification of PMS severity marker 1 SSL collection 38 healthy subjects (female 20-45 years old, BMI 18.5 or more and less than 25.0) were used as subjects. In healthy subjects, it was confirmed in advance that the menstrual cycle was stable every 25 to 30 days. Each period of the menstrual cycle (menstrual phase: any day 2-4 days after the start of menstruation, follicular phase: any day 10-12 days after the start of menstruation, luteal phase: 21-23 days after the start of menstruation On that day), sebum was collected from the entire face of each subject using a degreasing film (5.0 cm ⁇ 8.0 cm, 3M company). The oil removal film was transferred to a glass vial and stored at ⁇ 80 ° C.
  • RNA preparation and sequencing The oil-removing film from 1) above was cut to an appropriate size, and RNA was extracted using QIAzol Lysis Reagent (Qiagen) according to the attached protocol. The extracted RNA was reverse transcribed under the conditions of 42 ° C. for 105 minutes, and the obtained cDNA was amplified by multiplex PCR. Multiplex PCR was performed at an annealing and elongation temperature of 62 ° C. The obtained PCR product was purified by Aple XP (Beckman Coulter, Inc.).
  • the solution of the obtained purified product is attached to the 5 ⁇ VILO RT Reaction Mix attached to the SuperScript VILO cDNA Synthesis kit and the Ion AmpliSeqTranscriptome Human Gene Expression Prime 5 ⁇ Hixnipion Primer (Life Technologies Japan Co., Ltd.). Mix and Ion Ampliseq Transcriptome Human Gene Expression Core Panel were mixed, and then the primer sequences were digested, adapter-ligated and purified, and amplified according to the protocol attached to each kit to prepare a library. The prepared library was loaded into an Ion 540 Chip and sequenced using an Ion S5 / XL system (Life Technologies Japan KK). The gene from which each read sequence was derived was determined by gene mapping each read sequence obtained by sequencing to hg19 AmpliSeq Transcriptome ERCC v1, which is a reference sequence of the human genome.
  • the gene that serves as is extracted.
  • the missing values were supplemented by Singular Value Decomposition (SVD) imputation.
  • SVD Singular Value Decomposition
  • Fifty-one genes whose expression fluctuated between groups were obtained (Table 5), of which 40 genes shown in Table 6 were genes for which no association with PMS was previously reported.
  • Nucleic acids derived from the genes in Table 5 and translation products of the genes can be used as markers of PMS severity, respectively.
  • nucleic acids derived from the genes shown in Table 6 and translation products of the genes are novel PMS severity markers.
  • Example 2 Construction of PMS severity prediction model based on luteal phase RNA expression level data
  • the RNA expression level data (read count RPM value) obtained in Example 1 is converted into the logarithmic value of the base 2 (Log 2 RPM). Value).
  • Method 2 Selection of features From the 40 genes shown in Table 6 for which no association with PMS has been reported among the genes extracted in Example 1, the features used for building a predictive model are selected by the following two methods. bottom.
  • Method 1 Among the genes shown in Table 6, 5 to 15 kinds of genes from those having a small p-value of the expression difference between the severe PMS group and the mild PMS group, specifically, No. 1 in Table 6. 5 types from 1 to 5, No. in Table 6 10 types from 1 to 10, or No. 1 in Table 6.
  • Expression level data (Log 2 RPM value) of 15 kinds of genes from 1 to 15 were selected as feature amounts for PMS severity detection.
  • Method 2 The gene expression level data shown in Table 6 was subjected to principal component analysis using the tidyverse package of the statistical analysis environment R, and variable compression was performed. Data with a new dimension (1st to 10th principal components, Table 7) obtained by principal component analysis were selected as features for PMS severity detection.
  • Model construction was performed using the carnet package of the statistical analysis environment R.
  • a prediction model was constructed using the features selected by methods 1 and 2 in 2) as explanatory variables and the PMS severity (severe or mild) as the objective variable.
  • Four algorithms for each feature Random Forest, Neural network, Linear kernel support vector machine (SVM (linear)), and rbf kernel support vector machine (SVM (rbf)).
  • SVM Linear kernel support vector machine
  • the prediction model was trained by performing 10-fold cross-validation using.
  • RNA expression level data Log 2 RPM value
  • principal component data of test data was input to the model after training to calculate the predicted value of PMS severity.
  • the accuracy rate which is a value indicating the degree of agreement between the predicted value and the measured value, was calculated, and the model with the largest value was selected as the optimum predicted model.
  • Results Table 8 shows the features used to detect the severity of PMS, the optimal algorithm with the maximum accuracy, and its accuracy. Based on the expression level data of genes related to PMS severity in the luteal phase or their main components, PMS severity could be detected with high accuracy. Therefore, it was shown that the severity of PMS can be detected using the data on the expression level of SSL-derived RNA in the luteal phase. Further, in the case of the prediction model based on the main component, the accuracy of the prediction model based on the first to fifth main components is higher than that of the prediction model based on the main component 1, while it is equivalent to the prediction model based on the first to sixth main components. there were.
  • one or more principal components such that the total contribution rate is 20% or more, for example, 22 to 55% or more, preferably 50 to 70%, and more preferably 55 to 65% are used as explanatory variables.
  • a more accurate PMS severity prediction model can be constructed.
  • the accuracy of the prediction model using 10 genes was higher than that of the prediction model using 5 genes, while it was equivalent to the prediction model using 15 genes.
  • a more accurate PMS severity prediction model can be constructed by using the expression data of about 10 to 15 genes or more genes from the one with the largest expression difference depending on the PMS severity. ..
  • Example 3 Construction of PMS Severity Prediction Model Based on Follicular Phase RNA Expression Data 1) Dataset Division and Feature Selection Using the subject's follicular phase SSL-derived RNA expression dataset for the genes shown in Table 6. Except for the above, training data and test data were created in the same procedure as in 1) of Example 2. Next, in the same procedure as in 2) of Example 2, the 1st to 15th principal components (Table 9) obtained by the principal component analysis of the Log 2 RPM value of the RNA expression level data and the expression level data were subjected to PMS severity. Selected as a feature for detection.
  • Results Table 10 shows the features used to detect the severity of PMS, the optimal algorithm with the maximum accuracy, and its accuracy. Based on the expression level data of genes related to PMS severity in the follicular phase or their main components, PMS severity could be detected with high accuracy. Therefore, it was shown that the severity of PMS associated with the next menstruation can be detected by using the data on the expression level of SSL-derived RNA in the follicular phase. Further, in the case of the prediction model based on the main components, the accuracy of the prediction model based on the first to fifteenth main components is the highest, and the total contribution rate is 80% or more, preferably 85% or more, more preferably 90% or more.
  • PMS severity prediction model can be constructed by using one or more principal components as explanatory variables.
  • a prediction model constructed directly from the gene expression level it is more accurate to use expression data of about 5 to 20, preferably about 5 to 15, genes having a high correlation with PMS severity. It was shown that a high PMS severity prediction model can be constructed.

Landscapes

  • Chemical & Material Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Proteomics, Peptides & Aminoacids (AREA)
  • Health & Medical Sciences (AREA)
  • Organic Chemistry (AREA)
  • Wood Science & Technology (AREA)
  • Analytical Chemistry (AREA)
  • Zoology (AREA)
  • Genetics & Genomics (AREA)
  • Engineering & Computer Science (AREA)
  • Pathology (AREA)
  • Immunology (AREA)
  • Microbiology (AREA)
  • Molecular Biology (AREA)
  • Biotechnology (AREA)
  • Biophysics (AREA)
  • Physics & Mathematics (AREA)
  • Biochemistry (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • General Engineering & Computer Science (AREA)
  • General Health & Medical Sciences (AREA)
  • Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
  • Investigating Or Analysing Biological Materials (AREA)

Abstract

簡便且つ被験者への負担の少ない月経前症候群の重症度の検出方法の提供。月経前症候群の重症度の検出を可能にするマーカー。該マーカーを用いた被験体の月経前症候群の重症度の検出方法。

Description

月経前症候群の重症度検出方法
 本発明は、月経前症候群の重症度を検出する方法に関する。
 女性の月経前及び/又は月経中の身体的及び精神的不快は、月経随伴症状と呼ばれ、月経前症候群(Premenstrual Syndrome;以下、PMSという)、及び月経困難症に代表される。PMSは女性の労働生産性にも関係することが指摘されていることから、PMSや月経困難症に代表される月経随伴症状を緩和又は軽減し、女性のQOL(Quality of Life)を向上させることは、保健的だけでなく、社会的にも重要である。
 生体試料中のDNAやRNA等の核酸の解析によりヒトの生体内の現在、さらには将来の生理状態を調べる技術が開発されている。核酸を用いた解析は、網羅的な解析方法が確立されており一度の解析で豊富な情報を得られること、及び一塩基多形やRNA機能などに関する多くの研究報告に基づいて解析結果の機能的な紐付けが容易であることといった利点を有する。生体由来の核酸は、血液等の組織、体液、分泌物などから抽出することができる。特許文献1には、皮膚表上脂質(skin surface lipids;SSL)に含まれるRNAを生体の解析用の試料として用いること、SSLから表皮、汗腺、毛包及び皮脂腺のマーカー遺伝子が検出されたことが記載されている。
 女性ホルモンであるエストロゲン又はプロゲステロンの産生に伴って発現が変化する遺伝子が存在する。例えばこれまでに、AQP3(aquaporin 3)、HIF-1α(hypoxia inducible factor 1 subunit α)、MFN2(mitofusin 2)、DEFB4A(defensin beta 4A、又はhBD-2)、hBD-3などが報告されている(非特許文献1~3、特許文献2)。
(特許文献1)国際公開公報第2018/008319号
(特許文献2)特開2015-228829号公報
(非特許文献1)Hum Reprod, 2018, 33(11):2060-2073
(非特許文献2)Biol Reprod, 2018, 99(2):308-318
(非特許文献3)Mol Cell Endocrinol, 2013, 371:79-86
 本発明は、下記に示す遺伝子:SH3BGRL3、NDUFA11、PSME2、SNORA24、SSR4、ZNF91、SNORA70、HSPA6、ALDH2、BRK1、C20orf111、CSNK1A1、STK10、GHITM、RPL21、MED13、RPL32、UBE2R2、RPS3、PLEKHM2、DDX6、ALDH1A3、TSPYL1、NUB1、KRT6A、CCL22、ELOVL3、TMEM66、GDI2、MYL12B、CCNI、PSMC6、RABGEF1、ARF1、ACTN1、CHMP2A、GAK、STX3、SERP1及びGAS7、ならびに該遺伝子の翻訳産物からなる群より選択される少なくとも1つの発現を測定することを含む、被験体のPMSの重症度を検出する方法を提供する。
 また本発明は、下記遺伝子に由来する核酸:SH3BGRL3、NDUFA11、PSME2、SNORA24、SSR4、ZNF91、SNORA70、HSPA6、ALDH2、BRK1、C20orf111、CSNK1A1、STK10、GHITM、RPL21、MED13、RPL32、UBE2R2、RPS3、PLEKHM2、DDX6、ALDH1A3、TSPYL1、NUB1、KRT6A、CCL22、ELOVL3、TMEM66、GDI2、MYL12B、CCNI、PSMC6、RABGEF1、ARF1、ACTN1、CHMP2A、GAK、STX3、SERP1及びGAS7、ならびに該遺伝子の翻訳産物からなる群より選択される少なくとも1つを含む、PMSの重症度のマーカーを提供する。
発明の詳細な説明
 本明細書中で引用された全ての特許文献、非特許文献、及びその他の刊行物は、その全体が本明細書中において参考として援用される。
 月経前症候群(Premenstrual Syndrome;PMS)は、1931年にFrank(Archives of Neurology and Psychiatry, 1931, 26:1053)により初めて報告された、月経の1週間程前(黄体期)に発症する身体的症状群又は精神的症状群である。医学的には、PMSは、「月経開始の3~10日前から始まる身体的、精神的症状で月経開始とともに減退ないし消失するもの」と定義される。PMSの代表的な症状は、身体的症状群としては、下腹痛、腰痛、下腹部が張る、頭痛、肩こり、手足の冷え、食欲が増す、下痢、便秘、むくみ、乳房が痛い、乳房が張る、ニキビができやすい、肌荒れ、疲れやすい、眠くなる、等が挙げられ、精神的症状群としては、イライラする、怒りやすい、攻撃的になる、憂うつ、涙もろい、不安が高まる、等が挙げられる。
 本明細書において、PMSが重症であるとは、上記PMSの身体的症状又は精神的症状が重度である状態をいい、PMSが軽症であるとは、上記PMSの身体的症状又は精神的症状が軽度である状態をいう。従来、PMS等の月経随伴症状の重症度は、例えばMenstrual Distress Questionnaire(MDQ)のアンケートのスコアに基づいて評価されていた。黄体期におけるMDQスコアが高いほど、PMSの重症度は高いと評価される。例えば、心理測定尺度集III(サイエンス社、2001発行)に記載されたMDQに従って、黄体期のMDQスコアが66点以上である場合をPMSが重症であると評価し、48点以下である場合をPMSが軽症であると評価することができる。
 本発明は、PMSの重症度の検出を可能にするPMS重症度マーカー、及び該マーカー用いたPMSの重症度の検出方法に関する。
 本発明は、新規なマーカーを含むPMS重症度マーカーを用いてPMSの重症度を検出する方法を提供する。本発明で使用するPMS重症度マーカーは、皮膚表上脂質(SSL)から収集可能であるため、簡便且つ非侵襲的に取得することができる。したがって本発明によれば、被験者に負担をかけずにPMSの重症度を検出することが可能である。
 本発明は、従来は上記のような黄体期における被験者の問診に基づいて評価しなければならなかったPMSの重症度を、遺伝子又はその翻訳産物の発現に基づいてより簡易かつ客観的に検出することを可能にする。また本発明は、PMSが起こる黄体期に先立つ卵胞期における遺伝子又はその翻訳産物の発現に基づいてPMSの重症度を事前に検出することを可能にする。
 一態様において、本発明は、PMS重症度マーカーを提供する。後述の実施例に示されるとおり、本発明者は、下記表1に示す51種の遺伝子の発現が、PMSの重症度との間に関連性を有することを見出した。さらに、これら51種のうち表2に示した40種の遺伝子は、従来PMSとの関連性が報告されていない遺伝子であった。したがって、表1に示す遺伝子に由来する核酸、及び該遺伝子の翻訳産物は、各々、PMSの重症度のマーカーとして使用することができる。また、表2に示す遺伝子に由来する核酸、及び該遺伝子の翻訳産物は、新規のPMS重症度マーカーである。
Figure JPOXMLDOC01-appb-T000001
Figure JPOXMLDOC01-appb-T000002
 本明細書中に開示される遺伝子の名称は、NCBI([www.ncbi.nlm.nih.gov/])に記載のあるOfficial Symbolに従う。本明細書において、「遺伝子に由来する核酸」は、該遺伝子から転写されたmRNA、該mRNAから調製されたcDNA、該cDNAからの反応産物(例えばPCR産物、クローンDNA等)、などを包含する。
 本発明においては、表1に示す各遺伝子又はその翻訳産物の発現量に基づいて被験体のPMSの重症度の検出が可能である。本発明において、表1に示す遺伝子又はその翻訳産物の発現量は、該遺伝子に由来する核酸又は該遺伝子の翻訳産物を定量することで測定することができる。例えば、該遺伝子から転写されたmRNA、該mRNAから調製されたcDNA、該cDNAからの反応産物、又は該mRNAが翻訳されて合成される翻訳産物(タンパク質等)を定量することで、該遺伝子又はその翻訳産物の発現量を測定することができる。
 表1に示す遺伝子に由来する核酸又は該遺伝子の翻訳産物は、被験体から採取した生体試料、例えば細胞、体液、分泌物等から常法に従って調製することができる。好ましくは、該核酸は、被験体の皮膚表上脂質(skin surface lipids;SSL)から調製されるmRNAである。また好ましくは、該翻訳産物は、被験体のSSLから調製される翻訳産物である。
 本発明において、月経前症候群の重症度の「検出」とは、予測、検査、測定、判定又は評価支援などの用語で言い換えることもできる。なお、本発明における月経前症候群の重症度の「検出」、「予測」、「検査」、「測定」、「判定」又は「評価」という用語は、医師による月経前症候群の重症度の診断を含むものではない。
 本明細書において、「皮膚表上脂質(SSL)」とは、皮膚の表上に存在する脂溶性画分をいい、皮脂と呼ばれることもある。一般に、SSLは、皮膚にある皮脂腺等の外分泌腺から分泌された分泌物を主に含み、皮膚表面を覆う薄い層の形で皮膚表上に存在している。SSLは、皮膚細胞で発現したRNAを含む(特許文献1参照)。また本明細書において、「皮膚」とは、特に限定しない限り、体表の表皮、真皮、毛包、ならびに汗腺、皮脂腺及びその他の腺などの組織を含む領域の総称である。
 被験体の皮膚からのSSLの採取には、皮膚からのSSLの回収又は除去に用いられているあらゆる手段を採用することができる。好ましくは、後述するSSL吸収性素材、SSL接着性素材、又は皮膚からSSLをこすり落とす器具を使用することができる。SSL吸収性素材又はSSL接着性素材としては、SSLに親和性を有する素材であれば特に限定されず、例えばポリプロピレン、パルプ等が挙げられる。皮膚からのSSLの採取手順のより詳細な例としては、あぶら取り紙、あぶら取りフィルム等のシート状素材へSSLを吸収させる方法、ガラス板、テープ等へSSLを接着させる方法、スパーテル、スクレイパー等によりSSLをこすり落として回収する方法、などが挙げられる。SSLの吸着性を向上させるため、脂溶性の高い溶媒を予め含ませたSSL吸収性素材を用いてもよい。一方、SSL吸収性素材は、水溶性の高い溶媒や水分を含んでいるとSSLの吸着が阻害されるため、水溶性の高い溶媒や水分の含有量が少ないことが好ましい。SSL吸収性素材は、乾燥した状態で用いることが好ましい。SSLが採取される皮膚の部位としては、特に限定されず、頭、顔、首、体幹、手足等の身体の任意の部位の皮膚が挙げられ、皮脂の分泌が多い部位、例えば顔の皮膚が好ましい。
 被験体から採取されたSSLは一定期間保存されてもよい。採取されたSSLは、含有するRNAの分解を極力抑えるために、採取後できるだけ速やかに低温条件で保存することが好ましい。本発明における該RNA含有SSLの保存の温度条件は、0℃以下であればよく、好ましくは-20±20℃~-80±20℃、より好ましくは-20±10℃~-80±10℃、さらに好ましくは-20±20℃~-40±20℃、さらに好ましくは-20±10℃~-40±10℃、さらに好ましくは-20±10℃、さらに好ましくは-20±5℃である。該RNA含有SSLの該低温条件での保存の期間は、特に限定されないが、好ましくは12か月以下、例えば6時間以上12ヶ月以下、より好ましくは6ヶ月以下、例えば1日間以上6ヶ月以下、さらに好ましくは3ヶ月以下、例えば3日間以上3ヶ月以下である。
 採取したSSLからのRNAの抽出には、生体試料からのRNAの抽出又は精製に通常使用される方法、例えば、フェノール/クロロホルム法、AGPC(acid guanidinium thiocyanate-phenol-chloroform extraction)法、又はTRIzol(登録商標)、RNeasy(登録商標)、QIAzol(登録商標)などのカラムを用いた方法、シリカをコーティングした特殊な磁性体粒子を用いる方法、Solid Phase Reversible Immobilization磁性体粒子を用いる方法、ISOGEN等の市販のRNA抽出試薬による抽出などを用いることができる。また、SSLからの翻訳産物の抽出には、生体試料からのタンパク質の抽出又は精製に通常使用される方法、例えば、QIAzol Lysis Reagent(Qiagen)等の市販のタンパク質抽出試薬を用いることができる。
 以上の手順で調製された核酸又は翻訳産物を用いて表1に示す1つ以上の遺伝子又はその翻訳産物の発現を調べることで、被験体のPMSの重症度を検出することができる。したがって、別の一態様において、本発明は、表1に示す遺伝子又はその翻訳産物からなる群より選択される少なくとも1つの発現を測定することを含む、被験体のPMSの重症度を検出する方法を提供する。
 本発明によるPMSの重症度を検出する方法(以下、本発明の方法)における被験体は、例えばヒト又は非ヒト哺乳動物の雌性であって、PMSの重症度の検出を必要とする者であればよい。好ましくは、該被験体は、皮膚上にSSLを有する動物であり、より好ましくはヒトである。
 一実施形態において、本発明の方法は、被験体における表1に示す遺伝子からなる群より選択される少なくとも1つの遺伝子(以下、標的遺伝子(群)ともいう)の発現を測定することを含む。別の一実施形態において、本発明の方法は、被験体における表1に示す遺伝子の翻訳産物からなる群より選択される少なくとも1つ(以下、標的産物(群)ともいう)の発現を測定することを含む。あるいは、標的遺伝子(群)の発現と標的産物(群)の発現の両方が測定されてもよい。好ましくは、標的遺伝子(群)の発現が測定される。
 表1に示す遺伝子のなかでも、表2に示すものは、従来PMSとの関連性が報告されていない遺伝子である。したがって、本発明の方法の好ましい実施形態において、該標的遺伝子(群)は、表2に示す遺伝子からなる群より選択される少なくとも1つであり、該標的産物(群)は、表2に示す遺伝子の翻訳産物からなる群より選択される少なくとも1つである。さらに、本発明の方法においては、表2に示す遺伝子又はその翻訳産物からなる群より選択される少なくとも1つと、他の遺伝子(すなわち、表1に記載のCTNNB1、AQP3、ANXA1、ARHGDIB、DEFB4A、BIRC3、S100A7、SFN、HMGN2、HIF1A、及びMFN2)又はそれらの翻訳産物からなる群より選択される少なくとも1つとの組み合わせが標的遺伝子群又は標的産物群として用いられ、それらの発現が測定されてもよい。
 好ましくは、本発明の方法における該標的遺伝子(群)又は標的産物(群)の発現の測定は、被験体のSSLに含まれるRNA又は翻訳産物を用いて行われる。一実施形態において、本発明の方法は、被験体のSSLを採取することをさらに含んでいてもよい。一実施形態において、本発明の方法は、該被験体から採取したSSLからRNA又は翻訳産物を抽出することをさらに含んでいてもよい。
 被験体由来の試料、例えばSSLから抽出されたRNAに含まれる標的遺伝子由来のRNAの量を定量することで、標的遺伝子の発現量を測定することができる。被験体由来RNAの量の測定は、当該分野で通常用いられるRNAを用いた遺伝子発現解析の手順に従って実施すればよい。RNAを用いた遺伝子発現解析の手法の例としては、RNAを逆転写によりcDNAに変換した後、リアルタイムPCR、マルチプレックスPCR、マイクロアレイ、シーケンシング、クロマトグラフィーなどにより該cDNA又はその増幅産物を定量する方法が挙げられる。翻訳産物の発現量は、当該分野で通常用いられるタンパク質定量法、例えばELISA、免疫染色、蛍光法、電気泳動、クロマトグラフィー、質量分析などを用いて測定することができる。測定される発現量は、生体試料における標的遺伝子(群)又は標的産物(群)の発現の絶対量であってもよく、又は、他の標準物質や、遺伝子又は翻訳産物の総発現量に対する相対発現量であってもよい。
 好ましい実施形態において、本発明の方法では、該被験体由来RNAを逆転写によりcDNAに変換し、次いで該cDNAをPCRにかけ、得られた反応産物を精製する。該逆転写には、解析したい特定のRNAを標的としたプライマーを用いてもよいが、より包括的な核酸の保存及び解析のためにはランダムプライマーを用いることが好ましい。該逆転写には、一般的な逆転写酵素又は逆転写試薬キットを使用することができる。好適には、正確性及び効率性の高い逆転写酵素又は逆転写試薬キットが用いられ、その例としては、M-MLV Reverse Transcriptase及びその改変体、あるいは市販の逆転写酵素又は逆転写試薬キット、例えばPrimeScript(登録商標)Reverse Transcriptaseシリーズ(タカラバイオ社)、SuperScript(登録商標)Reverse Transcriptaseシリーズ(Thermo Scientific社)等が挙げられる。SuperScript(登録商標)III Reverse Transcriptase、SuperScript(登録商標)VILO cDNA Synthesis kit(いずれもThermo Scientific社)などが好ましく用いられる。得られたcDNAのPCRでは、解析したい特定のDNAを標的としたプライマーペアを用いて該特定のDNAのみを増幅してもよいが、複数のプライマーペアを用いて複数のDNAを増幅してもよい。好ましくは、該PCRはマルチプレックスPCRである。マルチプレックスPCRは、PCR反応系に複数のプライマー対を同時に使用することで、複数の遺伝子領域を同時に増幅する方法である。マルチプレックスPCRは、市販のキット(例えば、Ion AmpliSeqTranscriptome Human Gene Expression Kit;ライフテクノロジーズジャパン株式会社等)を用いて実施することができる。
 当該逆転写及びPCRの反応条件を調整することで、PCR反応産物の収率がより向上し、ひいてはそれを用いた解析の精度がより向上する。好適には、該逆転写での伸長反応において、温度を好ましくは42℃±1℃、より好ましくは42℃±0.5℃、さらに好ましくは42℃±0.25℃に調整し、一方、反応時間を好ましくは60分間以上、より好ましくは80~120分間に調整する。また好適には、該PCRにおけるアニーリング及び伸長反応の温度は、使用するプライマーに依存して適宜調整され得るが、例えば上記のマルチプレックスPCRキットを用いる場合、好ましくは62℃±1℃、より好ましくは62℃±0.5℃、さらに好ましくは62℃±0.25℃である。したがって、該PCRでは、好ましくはアニーリング及び伸長反応が1ステップで行われる。該アニーリング及び伸長反応のステップの時間は、増幅すべきDNAのサイズ等に依存して調整され得るが、好ましくは14~18分間である。該PCRにおける変性反応の条件は、増幅すべきDNAに依存して調整され得るが、好ましくは95~99℃で10~60秒間である。上記のような温度及び時間での逆転写及びPCRは、一般的にPCRに使用されるサーマルサイクラーを用いて実行することができる。
 当該PCRで得られた反応産物の精製は、反応産物のサイズ分離によって行われることが好ましい。サイズ分離により、目的のPCR反応産物を、PCR反応液中に含まれるプライマーやその他の不純物から分離することができる。DNAのサイズ分離は、例えば、サイズ分離カラムや、サイズ分離チップ、サイズ分離に利用可能な磁気ビーズなどによって行うことができる。サイズ分離に利用可能な磁気ビーズの好ましい例としては、Ampure XP等のSolid Phase Reversible Immobilization(SPRI)磁性ビーズが挙げられる。
 精製したPCR反応産物に対して、その後の定量解析を行うために必要なさらなる処理を施してもよい。例えば、DNAのシーケンシングのために、精製したPCR反応産物を、適切なバッファー溶液へと調製したり、PCR増幅されたDNAに含まれるPCRプライマー領域を切断したり、増幅されたDNAにアダプター配列をさらに付加したりしてもよい。例えば、精製したPCR反応産物をバッファー溶液へと調製し、増幅DNAに対してPCRプライマー配列の除去及びアダプターライゲーションを行い、得られた反応産物を、必要に応じて増幅して、定量解析のためのライブラリーを調製することができる。これらの操作は、例えば、SuperScript(登録商標)VILO cDNA Synthesis kit(ライフテクノロジーズジャパン株式会社)に付属している5×VILO RT Reaction Mix、及びIon AmpliSeq Transcriptome Human Gene Expression Kit(ライフテクノロジーズジャパン株式会社)に付属している5×Ion AmpliSeq HiFi Mix、及びIon AmpliSeq Transcriptome Human Gene Expression Core Panelを用いて、各キット付属のプロトコルに従って行うことができる。シーケンシングには、好ましくは次世代シーケンサー(例えばIon S5/XLシステム、ライフテクノロジーズジャパン株式会社)が用いられる。シーケンシングで作成されたリードの数(リードカウント)に基づいて、RNA発現を定量することができる。
 好ましくは、本発明の方法では、該標的遺伝子(群)又は該標的産物(群)の発現量に基づいて、被験体のPMSの重症度を検出する。例えば、表1に示す遺伝子は、PMSの重症度に伴うその発現量の変動のパターンによって表3及び表4に示す2つの群に分けることができる。これらの変動のパターンを利用して被験体のPMSの重症度を検出することができる。
 表3に示す遺伝子は、PMSが重症のときに発現が増加し、PMSが軽症のときに発現が低下する遺伝子である。一方、表4に示す遺伝子は、PMSが重症のときに発現が低下し、PMSが軽症のときに発現が増加する遺伝子である。表3及び表4に示す遺伝子又はその翻訳産物のいずれか1つ以上についての発現量を基準にすることで、被験体のPMSの重症度(例えば重症か否か)を検出することができる。表3及び表4に示す遺伝子のうち、下線で示した、それぞれ22種及び18種の遺伝子は、従来PMSとの関連性が報告されていない遺伝子(表2に示す遺伝子)である。
Figure JPOXMLDOC01-appb-T000003
Figure JPOXMLDOC01-appb-T000004
 本発明の方法の一実施形態においては、個々の標的遺伝子(群)又は標的産物(群)の発現量を直接、PMSの重症度を検出するためのパラメーターとして用いる。例えば、月経周期の特定の時期(例えば卵胞期又は黄体期、以下同じ)における被験体の標的遺伝子又は標的産物の発現量を定期的に(例えば毎月)測定し、その発現量の変動を追跡することで、該被験体のPMSの重症度を検出することができる。あるいは、月経周期の特定の時期において被験体から測定された標的遺伝子又は標的産物の発現量を、月経周期の同じ時期についての該標的遺伝子又は標的産物の基準発現量と比べることで、該被験体のPMSの重症度を検出することができる。標的遺伝子又は標的産物の基準発現量には、予め決定した値、例えば、該標的遺伝子又は標的産物の月経周期の特定の時期での発現量の統計値(例えば平均発現量)を用いることができる。該基準発現量は、被験体ごとに設定されてもよい。例えば、同じ被験体から複数回測定した月経周期の特定の時期における標的遺伝子又は標的産物の発現量の統計値(例えば平均値など)から算出してもよい。あるいは、該基準発現量は、集団から測定した月経周期の特定の時期における標的遺伝子又は標的産物の発現量の統計値(例えば平均値など)であってもよい。標的遺伝子群又は標的産物群として複数の遺伝子又は翻訳産物を用いる場合は、該遺伝子又は翻訳産物の各々について基準発現量を求めることが好ましい。
 例えば、標的遺伝子又は標的産物が表3に示す遺伝子又はその翻訳産物の場合、その発現量が基準発現量よりも高い場合、被験体のPMSは重症と検出され得、その発現量が基準発現量よりも低い場合、被験体のPMSは重症ではない(又は軽症)と検出され得る。一方、標的遺伝子又は標的産物が表4に示す遺伝子又はその翻訳産物の場合、その発現量が基準発現量よりも低い場合、被験体のPMSは重症と検出され得、その発現量が基準発現量よりも高い場合、被験体のPMSは重症ではない(又は軽症)と検出され得る。
 あるいは、表3~4に示す遺伝子又はその翻訳産物を2つ以上組み合わせて標的遺伝子又は標的産物として用いてもよい。複数の遺伝子又は翻訳産物を用いる場合には、一定割合、例えば50%以上、好ましくは70%以上、より好ましくは90%以上、さらに好ましくは100%の遺伝子又は翻訳産物が上述した発現量の基準を満たすか否かに基づいて、被験体のPMSの重症度を検出することができる。例えば、表3から選択した標的遺伝子群又は標的産物群の50%以上の発現量が基準発現量よりも高い場合、被験体のPMSは重症と検出され得る。あるいは、表4から選択した標的遺伝子群又は標的産物群の50%以上の発現量が基準発現量よりも低い場合、被験体のPMSは重症と検出され得る。あるいは、表3から選択した標的遺伝子群又は標的産物群の50%以上の発現量が基準発現量よりも高く、かつ表4から選択した標的遺伝子群又は標的産物群の50%以上の発現量が基準発現量よりも低い場合、被験体のPMSは重症と検出され得る。但し、本発明の方法における標的遺伝子群又は標的産物群の組み合わせは、上記の例に限定されない。
 別の一実施形態においては、標的遺伝子(群)又は標的産物(群)の発現データを用いて構築した予測モデルに基づいて、被験体のPMSの重症度を検出する。発現データとしては、月経周期の特定の時期(例えば卵胞期又は黄体期)における標的遺伝子(群)又は標的産物(群)の発現量についてのデータを用いることができる。例えば、1つ以上の標的遺伝子(群)又は標的産物(群)の各々についての発現データを説明変数とし、PMSの重症度(例えば、重症か軽症か)を目的変数とした機械学習により、PMSの重症度を検出するための最適な予測モデルを構築することができる。例えば、表1又は表2に示す遺伝子又はその翻訳産物の中からPMS重症度における重症群と軽症群のとの間で発現差のより大きいものを複数選択し、それらの発現データの各々を説明変数として用いることができる。
 あるいは、予測モデルの構築に複数の遺伝子又は翻訳産物の発現データを使用する場合、必要に応じて、次元削減によってデータを圧縮してから予測モデルの構築を行ってもよい。例えば、表1又は表2に示す遺伝子から複数の遺伝子を選択し、それらの遺伝子又はその翻訳産物を標的遺伝子群又は標的産物群として抽出する。次いで、該抽出した標的遺伝子群又は標的産物群の発現量について主成分分析を実施する。該主成分分析で算出された1つ以上の主成分を説明変数とし、PMSの重症度(例えば、重症か軽症か)を目的変数とした機械学習により、PMSの重症度を検出するための最適な予測モデルを構築することができる。
 予測モデルの構築に用いる標的遺伝子(群)又は標的産物(群)は、表1に示す遺伝子又はその翻訳産物のなかでPMSの重症度における重症群と軽症群との間での発現差がより大きいものから1種以上を選択するのが好ましい。例えば、従来PMSとの関連性が報告されていない、表2に示す40遺伝子から選択される1つ以上の遺伝子又はその翻訳産物が、予測モデルの構築に用いられる好ましい標的遺伝子(群)又は標的産物(群)として挙げられる。主成分分析に用いる標的遺伝子(群)又は標的産物(群)も同様である。
 一実施形態においては、表2に示す遺伝子の中から、黄体期においてPMS重症群と軽症群との間での発現差がより大きい遺伝子(例えば、表2のNo.1~20の遺伝子)から5~20個程度、好ましくは5~15個程度の遺伝子を選択する。より具体的には、例えば表2からNo.1~5の5遺伝子、No.1~10の10遺伝子又はNo.1~15の15遺伝子、好ましくはNo.1~10の10遺伝子又はNo.1~15の15遺伝子を選択する。選択した遺伝子の黄体期における発現量についてのデータを説明変数に用いて予測モデルを構築することができる。別の一実施形態においては、表2に示す遺伝子の中から、黄体期においてPMS重症群と軽症群との間での発現差がより大きい遺伝子(例えば、表2のNo.1~20の遺伝子)から5~20個程度、好ましくは5~15個程度の遺伝子を選択する。より具体的には、例えば表2からNo.1~5の5遺伝子、No.1~10の10遺伝子又はNo.1~15の15遺伝子、好ましくはNo.1~5の5遺伝子を選択する。選択した遺伝子の卵胞期における発現量についてのデータを説明変数に用いて予測モデルを構築することができる。
 主成分分析でデータ圧縮を行った場合、予測モデルの構築に用いる主成分の次元や数は、特に限定されないが、寄与率の大きい成分、例えば、第1主成分を含み、かつ合計寄与率が、好ましくは20%以上、より好ましくは50%以上になるような1つ以上の主成分を説明変数とするとよい。一実施形態においては、黄体期における標的遺伝子群の発現データに対して主成分分析を行い、第1主成分を含みかつ合計寄与率が20%以上、好ましくは50%以上、例えば、好ましくは22~55%又はそれ以上、より好ましくは50~70%、さらに好ましくは55~65%になるよう主成分を選択し、選択した主成分を説明変数に用いて予測モデルを構築することができる。別の一実施形態においては、卵胞期における標的遺伝子群の発現データに対して主成分分析を行い、第1主成分を含みかつ合計寄与率が、80%以上、好ましくは85%以上、より好ましくは90%以上になるよう主成分を選択し、選択した主成分を説明変数に用いて予測モデルを構築することができる。
 予測モデルの構築におけるアルゴリズムは、機械学習に用いるアルゴリズムなどの公知のものを利用することができる。機械学習アルゴリズムの例としては、ランダムフォレスト(Random Forest)、ニューラルネットワーク(Neural net)、線形カーネルのサポートベクターマシン(SVM (linear))、rbfカーネルのサポートベクターマシン(SVM (rbf))などのアルゴリズムが挙げられる。構築した予測モデルに検証用のデータを入力して予測値を算出し、該予測値が実測値と最も適合するモデル、例えば、予測値と実測値の一致度合いを示す値である正解率(Accuracy)が高いモデルを、最適な予測モデルとして選択することができる。
 構築された予測モデルに、被験体における該標的遺伝子(群)もしくは標的産物(群)の発現データ(又はそれらから算出された主成分)を入力することで、該被験体のPMSの重症度を検出することができる。
 本発明の例示的実施形態として、以下の物質、製造方法、用途、方法等をさらに本明細書に開示する。但し、本発明はこれらの実施形態に限定されない。
〔1〕下記に示す遺伝子:SH3BGRL3、NDUFA11、CTNNB1、AQP3、PSME2、SNORA24、SSR4、ZNF91、ANXA1、SNORA70、HSPA6、ALDH2、ARHGDIB、DEFB4A、BRK1、C20orf111、CSNK1A1、BIRC3、S100A7、STK10、GHITM、SFN、RPL21、MED13、RPL32、UBE2R2、RPS3、PLEKHM2、DDX6、ALDH1A3、TSPYL1、NUB1、KRT6A、CCL22、ELOVL3、TMEM66、GDI2、MYL12B、CCNI、PSMC6、HMGN2、HIF1A、RABGEF1、ARF1、ACTN1、MFN2、CHMP2A、GAK、STX3、SERP1及びGAS7、ならびに該遺伝子の翻訳産物からなる群より選択される少なくとも1つの発現を測定することを含む、被験体の月経前症候群の重症度を検出する方法。
〔2〕好ましくは、下記に示す遺伝子:SH3BGRL3、NDUFA11、PSME2、SNORA24、SSR4、ZNF91、SNORA70、HSPA6、ALDH2、BRK1、C20orf111、CSNK1A1、STK10、GHITM、RPL21、MED13、RPL32、UBE2R2、RPS3、PLEKHM2、DDX6、ALDH1A3、TSPYL1、NUB1、KRT6A、CCL22、ELOVL3、TMEM66、GDI2、MYL12B、CCNI、PSMC6、RABGEF1、ARF1、ACTN1、CHMP2A、GAK、STX3、SERP1及びGAS7、ならびに該遺伝子の翻訳産物からなる群より選択される少なくとも1つの発現を測定することを含む、〔1〕記載の方法。
〔3〕好ましくは、さらに、下記に示す遺伝子:CTNNB1、AQP3、ANXA1、ARHGDIB、DEFB4A、BIRC3、S100A7、SFN、HMGN2、HIF1A、及びMFN2、ならびに該遺伝子の翻訳産物からなる群より選択される少なくとも1つの発現を測定することを含む、〔2〕記載の方法。
〔4〕好ましくは、前記遺伝子又は翻訳産物の発現量に基づいて月経前症候群の重症度を検出することを含み、
 より好ましくは、前記遺伝子の発現量に基づいて月経前症候群の重症度を検出することを含む、
〔1〕~〔3〕のいずれか1項記載の方法。
〔5〕好ましくは、前記遺伝子又は翻訳産物の黄体期又は卵胞期での発現量に基づいて月経前症候群の重症度を検出することを含む、〔4〕記載の方法。
〔6〕好ましくは、前記遺伝子又は翻訳産物の発現量に基づく予測モデルを用いて月経前症候群の重症度を検出することを含み、
 該予測モデルが、月経前症候群の重症群及び軽症群における該遺伝子又は翻訳産物の発現量、又は該発現量の主成分分析で算出された少なくとも1種の主成分を、説明変数とし、月経前症候群の重症度を目的変数として構築される、
〔1〕~〔5〕のいずれか1項記載の方法。
〔7〕好ましくは、前記遺伝子の発現量が、前記被験体の皮膚表上脂質に含まれる該遺伝子のmRNAを定量することにより測定される、〔1〕~〔6〕のいずれか1項記載の方法。
〔8〕下記の遺伝子に由来する核酸:SH3BGRL3、NDUFA11、CTNNB1、AQP3、PSME2、SNORA24、SSR4、ZNF91、ANXA1、SNORA70、HSPA6、ALDH2、ARHGDIB、DEFB4A、BRK1、C20orf111、CSNK1A1、BIRC3、S100A7、STK10、GHITM、SFN、RPL21、MED13、RPL32、UBE2R2、RPS3、PLEKHM2、DDX6、ALDH1A3、TSPYL1、NUB1、KRT6A、CCL22、ELOVL3、TMEM66、GDI2、MYL12B、CCNI、PSMC6、HMGN2、HIF1A、RABGEF1、ARF1、ACTN1、MFN2、CHMP2A、GAK、STX3、SERP1及びGAS7、ならびに該遺伝子の翻訳産物からなる群より選択される少なくとも1つの、月経前症候群の重症度マーカーとしての使用。
〔9〕好ましくは、下記の遺伝子に由来する核酸:SH3BGRL3、NDUFA11、PSME2、SNORA24、SSR4、ZNF91、SNORA70、HSPA6、ALDH2、BRK1、C20orf111、CSNK1A1、STK10、GHITM、RPL21、MED13、RPL32、UBE2R2、RPS3、PLEKHM2、DDX6、ALDH1A3、TSPYL1、NUB1、KRT6A、CCL22、ELOVL3、TMEM66、GDI2、MYL12B、CCNI、PSMC6、RABGEF1、ARF1、ACTN1、CHMP2A、GAK、STX3、SERP1及びGAS7、ならびに該遺伝子の翻訳産物からなる群より選択される少なくとも1つが使用される、〔8〕記載の使用。
〔10〕好ましくは、さらに、下記の遺伝子に由来する核酸:CTNNB1、AQP3、ANXA1、ARHGDIB、DEFB4A、BIRC3、S100A7、SFN、HMGN2、HIF1A、及びMFN2、ならびに該遺伝子の翻訳産物からなる群より選択される少なくとも1つが使用される、〔9〕記載の使用。
〔11〕好ましくは、前記遺伝子に由来する核酸が、皮膚表上脂質に含まれる該遺伝子のmRNAである、〔8〕~〔10〕のいずれか1項記載の使用。
〔12〕下記の遺伝子に由来する核酸:SH3BGRL3、NDUFA11、PSME2、SNORA24、SSR4、ZNF91、SNORA70、HSPA6、ALDH2、BRK1、C20orf111、CSNK1A1、STK10、GHITM、RPL21、MED13、RPL32、UBE2R2、RPS3、PLEKHM2、DDX6、ALDH1A3、TSPYL1、NUB1、KRT6A、CCL22、ELOVL3、TMEM66、GDI2、MYL12B、CCNI、PSMC6、RABGEF1、ARF1、ACTN1、CHMP2A、GAK、STX3、SERP1及びGAS7、ならびに該遺伝子の翻訳産物からなる群より選択される少なくとも1つを含む、月経前症候群の重症度のマーカー。
〔13〕好ましくは、さらに、下記遺伝子に由来する核酸:CTNNB1、AQP3、ANXA1、ARHGDIB、DEFB4A、BIRC3、S100A7、SFN、HMGN2、HIF1A、及びMFN2、ならびに該遺伝子の翻訳産物からなる群より選択される少なくとも1つを含む、〔12〕記載のマーカー。
〔14〕好ましくは、前記遺伝子に由来する核酸が、皮膚表上脂質に含まれる該遺伝子のmRNAである、〔12〕又は〔13〕記載のマーカー。
 以下、実施例に基づき本発明をさらに詳細に説明するが、本発明はこれに限定されるものではない。
実施例1 PMS重症度マーカーの同定
1)SSL採取
 健常者(20~45歳女性、BMI18.5以上25.0未満)38名を被験者とした。健常者は、予め月経周期が約25~30日で毎回安定していることを確認されていた。月経周期の各時期(月経期:月経開始後2~4日のいずれかの日、卵胞期:月経開始後10~12日のいずれかの日、黄体期:月経開始後21~23日のいずれかの日)に、あぶら取りフィルム(5.0cm×8.0cm、3M社)を用いて各被験者の全顔から皮脂を回収した。該あぶら取りフィルムはガラスバイアルに移し、RNA抽出に使用するまで-80℃で、約1ヶ月間保存した。また、皮脂採取後に、月経関連症状を評価する質問紙である、心理測定尺度集III(サイエンス社、2001発行)に記載のMenstrual Distress Questionnaire(MDQ)を用いて、被験者の月経周期の各時期(月経期、卵胞期及び黄体期)の月経随伴症状(月経困難症、月経前症候群等)の重症度を評価した。
2)RNA調製及びシーケンシング
 上記1)のあぶら取りフィルムを適当な大きさに切断し、QIAzol Lysis Reagent(Qiagen)を用いて、付属のプロトコルに準じてRNAを抽出した。抽出したRNAは、42℃、105分間の条件で逆転写し、得られたcDNAをマルチプレックスPCRで増幅した。マルチプレックスPCRは62℃のアニーリング及び伸長温度で実施した。得られたPCR産物は、Ampure XP(ベックマン・コールター株式会社)で精製した。得られた精製産物の溶液を、SuperScript VILO cDNA Synthesis kitに付属している5×VILO RT Reaction Mix、及びIon AmpliSeqTranscriptome Human Gene Expression Kit(ライフテクノロジーズジャパン株式会社)に付属している5×Ion Ampliseq HiFi Mix、及びIon AmpliSeq Transcriptome Human Gene Expression Core Panelを混合し、次いで各キット付属のプロトコルに従ってプライマー配列の消化、アダプターライゲーションと精製、及び増幅を行い、ライブラリーを調製した。調製したライブラリーをIon 540 Chipにローディングし、Ion S5/XLシステム(ライフテクノロジーズジャパン株式会社)を用いてシーケンシングした。シーケンシングで得られた各リード配列をヒトゲノムのリファレンス配列であるhg19 AmpliSeq Transcriptome ERCC v1に対して遺伝子マッピングすることで、各リード配列の由来する遺伝子を決定した。
3)MDQスコアによる被験者の群分け
 MDQではスコアが高いほど月経随伴症状の重症度が高いことを示す。そのため、被験者38名を、黄体期に得られたMDQのスコアに従って、スコアが66点以上である15名(平均スコア87.4)をPMS重症群、スコアが48点以下である15名(平均スコア40.3)をPMS軽症群に群分けした。中間層である8名に関しては、以下のデータ解析から除外した。
4)データ解析
 上記2)で取得した卵胞期及び黄体期における被験者のSSL由来RNAのシーケンシングで得られた各リードのリードカウントを、各RNAの発現量のデータとし、リードカウントが10未満のリードは欠損値として扱った。サンプル間の総リードカウントの違いを補正するため、各リードカウントをRPM(Reads per million mapped reads)値に変換した。80%以上の被験者で欠損値ではない発現量データが得られた1884種の遺伝子を選択した。これら1884種の遺伝子について、PMSが発症する黄体期におけるそれらの発現量データ(RPM値)に基づいて群間(PMS重症群及び軽症群)でStudent’s t testを行い、p<0.05となる遺伝子を抽出した。なお、欠損値についてはSingular Value Decomposition(SVD)imputationにて補填した。群間で発現が変動する遺伝子が51個得られ(表5)、そのうち表6に示す40種の遺伝子が従来PMSとの関連が報告されていない遺伝子であった。これら表5中の遺伝子に由来する核酸、及び該遺伝子の翻訳産物は、各々、PMS重症度のマーカーとして使用できる。さらに、これら表5中の遺伝子のうち、表6に示す遺伝子に由来する核酸、及び該遺伝子の翻訳産物は、新規のPMS重症度マーカーである。
Figure JPOXMLDOC01-appb-T000005
Figure JPOXMLDOC01-appb-T000006
実施例2 黄体期RNA発現量データに基づくPMS重症度予測モデルの構築
1)データセット分割
 実施例1で得られた被験者の黄体期のSSL由来RNA発現量データセットのうち、全解析対象者数の67%に当たる各群(PMS重症群及び軽症群)10名分ずつについてのデータをPMS重症度予測モデルの訓練データとし、残り33%に当たる各群5名分ずつのデータをモデル精度の評価に使用するテストデータとした。負の二項分布に従うRPM値を正規分布に近似するため、実施例1で求めたRNAの発現量データ(リードカウントのRPM値)を、底2の対数値に変換したRPM値(Log2RPM値)に変換した。
2)特徴量の選択
 実施例1で抽出した遺伝子のうち、PMSとの関連が報告されていない表6に示す40個の遺伝子から、以下の2つの方法により予測モデル構築に用いる特徴量を選択した。
(方法1)表6に示す遺伝子のうち、PMS重症群と軽症群との発現差のp値が小さいものから5~15種の遺伝子、具体的には表6のNo.1~5の5種、表6のNo.1~10の10種、又は表6のNo.1~15の15種の遺伝子の発現量データ(Log2RPM値)を、PMS重症度検出の特徴量として選択した。
(方法2)表6に示す遺伝子の発現量データについて、統計解析環境Rのtidyverseパッケージを使用して主成分分析を行い、変数圧縮を実施した。主成分分析によって得られた、新しい次元を持つデータ(第1~第10主成分、表7)をPMS重症度検出の特徴量として選択した。
Figure JPOXMLDOC01-appb-T000007
3)モデル構築
 モデル構築は統計解析環境Rのcarnetパッケージを使用して行った。前記2)で方法1及び2により選択した特徴量を説明変数として用い、PMS重症度(重症又は軽症)を目的変数として用いて、予測モデルを構築した。特徴量毎に、ランダムフォレスト(Random Forest)、ニューラルネットワーク(Neural net)、線形カーネルのサポートベクターマシン(SVM (linear))、及びrbfカーネルのサポートベクターマシン(SVM (rbf))の4種のアルゴリズムを用いて、10倍交差検証を行って予測モデルを学習させた。各アルゴリズムについて、学習後のモデルにテストデータのRNAの発現量データ(Log2RPM値)又は主成分データを入力してPMS重症度の予測値を計算した。予測値と実測値の一致度合いを示す値である正解率(Accuracy)を計算し、その値が最も大きいモデルを最適な予測モデルとして選抜した。
4)結果
 表8に、PMS重症度の検出に用いた特徴量、最大のAccuracyを与えた最適アルゴリズム、及びそのAccuracyを示す。PMS重症度に関連する遺伝子の黄体期における発現量データ、又はそれらの主成分に基づいて、高いAccuracyでPMS重症度を検出することができた。したがって、SSL由来RNAの黄体期における発現量のデータを用いてPMS重症度の検出が可能であることが示された。また、主成分に基づく予測モデルの場合、第1~第5主成分による予測モデルのAccuracyは、主成分1による予測モデルよりも高く、一方、第1~第6主成分による予測モデルと同等であった。これらの結果から、合計寄与率が20%以上、例えば22~55%又はそれ以上、好ましくは50~70%、より好ましくは55~65%になるような1つ以上の主成分を説明変数とすることで、より精度の高いPMS重症度予測モデルが構築可能であることが示された。また、遺伝子発現量から直接構築した予測モデルの場合、10遺伝子による予測モデルのAccuracyは、5遺伝子による予測モデルよりも高く、一方、15遺伝子による予測モデルと同等であった。PMS重症度による発現差の大きいものから10~15個程度の遺伝子又はそれより多くの遺伝子の発現データを用いることで、より精度の高いPMS重症度予測モデルが構築可能であることが示された。
Figure JPOXMLDOC01-appb-T000008
実施例3 卵胞期RNA発現量データに基づくPMS重症度予測モデルの構築
1)データセット分割及び特徴量の選択
 表6に示す遺伝子についての被験者の卵胞期のSSL由来RNA発現量データセットを用いた以外は、実施例2の1)と同様の手順で訓練データとテストデータを作成した。次いで実施例2の2)と同様の手順で、RNAの発現量データのLog2RPM値及び発現量データの主成分分析によって得られた第1~第15主成分(表9)をPMS重症度検出の特徴量として選択した。
Figure JPOXMLDOC01-appb-T000009
2)モデル構築
 1)で選択した特徴量を説明変数として用いて、実施例2の3)と同様の手順で予測モデルを学習させ、正解率(Accuracy)が最も大きいモデルを最適な予測モデルとして選抜した。
3)結果
 表10に、PMS重症度の検出に用いた特徴量、最大のAccuracyを与えた最適アルゴリズム、及びそのAccuracyを示す。PMS重症度に関連する遺伝子の卵胞期における発現量データ、又はそれらの主成分に基づいて、高いAccuracyでPMS重症度を検出することができた。したがって、SSL由来RNAの卵胞期における発現量のデータを用いて、次回月経に伴うPMSの重症度の検出が可能であることが示された。また、主成分に基づく予測モデルの場合、第1~第15主成分による予測モデルのAccuracyが最も高く、合計寄与率が80%以上、好ましくは85%以上、より好ましくは90%以上になるような1つ以上の主成分を説明変数とすることで、より精度の高いPMS重症度予測モデルが構築可能であることが示された。一方、遺伝子発現量から直接構築した予測モデルの場合、PMS重症度との相関が高いものから5~20個程度、好ましくは5~15個程度の遺伝子の発現データを用いることで、より精度の高いPMS重症度予測モデルが構築可能であることが示された。
Figure JPOXMLDOC01-appb-T000010
 

Claims (10)

  1.  下記に示す遺伝子:SH3BGRL3、NDUFA11、PSME2、SNORA24、SSR4、ZNF91、SNORA70、HSPA6、ALDH2、BRK1、C20orf111、CSNK1A1、STK10、GHITM、RPL21、MED13、RPL32、UBE2R2、RPS3、PLEKHM2、DDX6、ALDH1A3、TSPYL1、NUB1、KRT6A、CCL22、ELOVL3、TMEM66、GDI2、MYL12B、CCNI、PSMC6、RABGEF1、ARF1、ACTN1、CHMP2A、GAK、STX3、SERP1及びGAS7、ならびに該遺伝子の翻訳産物からなる群より選択される少なくとも1つの発現を測定することを含む、被験体の月経前症候群の重症度を検出する方法。
  2.  さらに、下記遺伝子:CTNNB1、AQP3、ANXA1、ARHGDIB、DEFB4A、BIRC3、S100A7、SFN、HMGN2、HIF1A、及びMFN2、ならびに該遺伝子の翻訳産物からなる群より選択される少なくとも1つの発現を測定することを含む、請求項1記載の方法。
  3.  前記遺伝子又は翻訳産物の発現量に基づいて月経前症候群の重症度を検出することを含む、請求項1又は2記載の方法。
  4.  前記遺伝子又は翻訳産物の黄体期又は卵胞期での発現量に基づいて月経前症候群の重症度を検出することを含む、請求項3記載の方法。
  5.  前記遺伝子又は翻訳産物の発現量に基づく予測モデルを用いて月経前症候群の重症度を検出することを含み、
     該予測モデルが、月経前症候群の重症群及び軽症群における該遺伝子又は翻訳産物の発現量、又は該発現量の主成分分析で算出された少なくとも1種の主成分を、説明変数とし、月経前症候群の重症度を目的変数として構築される、
    請求項1~4のいずれか1項記載の方法。
  6.  前記遺伝子の発現量が、前記被験体の皮膚表上脂質に含まれる該遺伝子のmRNAを定量することにより測定される、請求項1~5のいずれか1項記載の方法。
  7.  下記の遺伝子に由来する核酸:SH3BGRL3、NDUFA11、PSME2、SNORA24、SSR4、ZNF91、SNORA70、HSPA6、ALDH2、BRK1、C20orf111、CSNK1A1、STK10、GHITM、RPL21、MED13、RPL32、UBE2R2、RPS3、PLEKHM2、DDX6、ALDH1A3、TSPYL1、NUB1、KRT6A、CCL22、ELOVL3、TMEM66、GDI2、MYL12B、CCNI、PSMC6、RABGEF1、ARF1、ACTN1、CHMP2A、GAK、STX3、SERP1及びGAS7、ならびに該遺伝子の翻訳産物からなる群より選択される少なくとも1つの、月経前症候群の重症度のマーカーとしての使用。
  8.  前記遺伝子に由来する核酸が、皮膚表上脂質に含まれる該遺伝子のmRNAである、請求項7記載の使用。
  9.  下記の遺伝子に由来する核酸:SH3BGRL3、NDUFA11、PSME2、SNORA24、SSR4、ZNF91、SNORA70、HSPA6、ALDH2、BRK1、C20orf111、CSNK1A1、STK10、GHITM、RPL21、MED13、RPL32、UBE2R2、RPS3、PLEKHM2、DDX6、ALDH1A3、TSPYL1、NUB1、KRT6A、CCL22、ELOVL3、TMEM66、GDI2、MYL12B、CCNI、PSMC6、RABGEF1、ARF1、ACTN1、CHMP2A、GAK、STX3、SERP1及びGAS7、ならびに該遺伝子の翻訳産物からなる群より選択される少なくとも1つを含む、月経前症候群の重症度のマーカー。
  10.  前記遺伝子に由来する核酸が、皮膚表上脂質に含まれる該遺伝子のmRNAである、請求項9記載のマーカー。 
PCT/JP2021/016470 2020-04-24 2021-04-23 月経前症候群の重症度検出方法 WO2021215532A1 (ja)

Priority Applications (3)

Application Number Priority Date Filing Date Title
CN202180030487.6A CN115443337A (zh) 2020-04-24 2021-04-23 经前综合征的重症度检测方法
US17/920,652 US20230144067A1 (en) 2020-04-24 2021-04-23 Method for detecting severity of premenstrual syndrome
EP21792472.9A EP4141114A1 (en) 2020-04-24 2021-04-23 Method for detecting severity of premenstrual syndrome

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2020-077036 2020-04-24
JP2020077036 2020-04-24

Publications (1)

Publication Number Publication Date
WO2021215532A1 true WO2021215532A1 (ja) 2021-10-28

Family

ID=78269345

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2021/016470 WO2021215532A1 (ja) 2020-04-24 2021-04-23 月経前症候群の重症度検出方法

Country Status (5)

Country Link
US (1) US20230144067A1 (ja)
EP (1) EP4141114A1 (ja)
JP (1) JP2021175395A (ja)
CN (1) CN115443337A (ja)
WO (1) WO2021215532A1 (ja)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2015228829A (ja) 2014-06-05 2015-12-21 ポーラ化成工業株式会社 抗菌剤の選択方法
WO2018008319A1 (ja) 2016-07-08 2018-01-11 花王株式会社 核酸試料の調製方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2015228829A (ja) 2014-06-05 2015-12-21 ポーラ化成工業株式会社 抗菌剤の選択方法
WO2018008319A1 (ja) 2016-07-08 2018-01-11 花王株式会社 核酸試料の調製方法
JP2018000206A (ja) * 2016-07-08 2018-01-11 花王株式会社 核酸試料の調製方法

Non-Patent Citations (10)

* Cited by examiner, † Cited by third party
Title
BAHRAMI, A ET AL.: "Menstrual problems in adolescence: relationship to serum vitamins A and E, and systemic inflammation", ARCH. GYNECOL. OBSTET., vol. 301, 16 November 2019 (2019-11-16), pages 189 - 197, XP037036182, DOI: 10.1007/s00404-019-05343-1 *
BERTONE-JOHNSON E. R., RONNENBERG A. G., HOUGHTON S. C., NOBLES C., ZAGARINS S. E., TAKASHIMA-UEBELHOER B. B., FARAJ J. L., WHITCO: "Association of inflammation markers with menstrual symptom severity and premenstrual syndrome in young women", HUMAN REPRODUCTION, OXFORD JOURNALS, GB, vol. 29, no. 9, 1 September 2014 (2014-09-01), GB , pages 1987 - 1994, XP055869562, ISSN: 0268-1161, DOI: 10.1093/humrep/deu170 *
BIOL. REPROD., vol. 99, no. 2, 2018, pages 308 - 318
DATABASE Nucleotide [online] 15 July 2006 (2006-07-15), "Homo sapiens SH3 domain binding glutamic acid- rich protein like 3, mRNA (cDNA clone MGC: 19619 IMAGE:4040656", XP055868305, retrieved from NCBI Database accession no. BC030135.2 *
DATABASE Nucleotide [online] 24 June 2018 (2018-06-24), "Homo sapiens C-C motif chemokine ligand 22 (CCL22", XP055869560, retrieved from NCBI Database accession no. NM_002990.4 *
DATABASE Nucleotide [online] 7 October 2008 (2008-10-07), "Homo sapiens mRNA for NADH-ubiquinone oxidoreductase subunit B14.7 (NDUFA11 gene", XP055869545, retrieved from NCBI Database accession no. AJ539081.1 *
DATABASE Nucleotide 16 October 2008 (2008-10-16), "Homo sapiens full open reading frame cDNA clone RZPDo834D0927D for gene PSME2, proteasome (prosome, macropain) activator subunit 2 (PA28 beta); complete eds, without stopcodon", XP055869554, retrieved from NCBI Database accession no. CR541657.1 *
HUM. REPROD., vol. 33, no. 11, 2018, pages 2060 - 2073
MOL. CELL. ENDOCRINOL., vol. 371, 2013, pages 79 - 86
See also references of EP4141114A1

Also Published As

Publication number Publication date
EP4141114A1 (en) 2023-03-01
US20230144067A1 (en) 2023-05-11
CN115443337A (zh) 2022-12-06
JP2021175395A (ja) 2021-11-04

Similar Documents

Publication Publication Date Title
Wong et al. Use of RT-PCR and DNA microarrays to characterize RNA recovered by non-invasive tape harvesting of normal and inflamed skin
JP2020074769A (ja) 被験体の皮膚細胞に由来する核酸の調製方法
Benson et al. An analysis of select pathogenic messages in lesional and non-lesional psoriatic skin using non-invasive tape harvesting
WO2021230379A1 (ja) パーキンソン病の検出方法
WO2021215532A1 (ja) 月経前症候群の重症度検出方法
Mosca et al. Cell-free DNA in the plasma of patients with systemic sclerosis
WO2021215531A1 (ja) 月経周期マーカー
JP2022026407A (ja) プラセボ効果を評価する方法
JP2008054590A (ja) 運動ストレス状態の評価方法
JP2023069499A (ja) 乾燥による皮膚表面形状悪化の検出方法
WO2022220299A1 (ja) 乳幼児顔湿疹の症度の検出方法
JP2022097302A (ja) 睡眠状態マーカー及びそれを用いた睡眠状態の検出方法
JP2022097303A (ja) 疲労マーカー及びそれを用いた疲労の検出方法
WO2022220298A1 (ja) 乳幼児おむつ皮膚炎の症度の検出方法
WO2022114200A1 (ja) アトピー性皮膚炎の症度の検出方法
JP2010115178A (ja) 皮膚の炎症の予測方法及びその用途
JP2024045048A (ja) 肥満改善剤の有効性予測方法
JP2023002299A (ja) 内臓脂肪面積の検出方法
JP2022174645A (ja) 生物学的年齢の予測方法
JP2023002298A (ja) 糖尿病指標値の検出方法
JP2023073135A (ja) 更年期障害の重症度の検出方法
JP2024045047A (ja) 肥満改善剤の有効性予測方法
WO2022009988A1 (ja) アトピー性皮膚炎の重症度の検出方法
WO2022114201A1 (ja) アトピー性皮膚炎の症度変化の検出方法
JP2022026406A (ja) 睡眠改善剤の有効性を評価する方法

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 21792472

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 2021792472

Country of ref document: EP

Effective date: 20221124

NENP Non-entry into the national phase

Ref country code: DE