WO2024025536A1 - Precision medicine for anxiety disorders: objective assessment, risk prediction, pharmacogenomics, and repurposed drugs - Google Patents
Precision medicine for anxiety disorders: objective assessment, risk prediction, pharmacogenomics, and repurposed drugs Download PDFInfo
- Publication number
- WO2024025536A1 WO2024025536A1 PCT/US2022/038673 US2022038673W WO2024025536A1 WO 2024025536 A1 WO2024025536 A1 WO 2024025536A1 US 2022038673 W US2022038673 W US 2022038673W WO 2024025536 A1 WO2024025536 A1 WO 2024025536A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- anxiety
- biomarkers
- subject
- score
- biomarker
- Prior art date
Links
- 208000019901 Anxiety disease Diseases 0.000 title claims abstract description 180
- 239000003814 drug Substances 0.000 title claims abstract description 41
- 230000002974 pharmacogenomic effect Effects 0.000 title abstract description 9
- 229940040939 repurposed drug Drugs 0.000 title description 2
- 239000000090 biomarker Substances 0.000 claims abstract description 121
- 230000036506 anxiety Effects 0.000 claims abstract description 120
- 230000014509 gene expression Effects 0.000 claims abstract description 62
- 229940079593 drug Drugs 0.000 claims abstract description 36
- 238000011282 treatment Methods 0.000 claims abstract description 22
- 239000008280 blood Substances 0.000 claims abstract description 18
- 210000004369 blood Anatomy 0.000 claims abstract description 17
- 230000008859 change Effects 0.000 claims abstract description 11
- 238000000034 method Methods 0.000 claims description 35
- 108090000623 proteins and genes Proteins 0.000 claims description 33
- 150000001875 compounds Chemical class 0.000 claims description 25
- 239000000523 sample Substances 0.000 claims description 17
- 238000002560 therapeutic procedure Methods 0.000 claims description 13
- 230000001225 therapeutic effect Effects 0.000 claims description 12
- 230000000116 mitigating effect Effects 0.000 claims description 9
- 239000013074 reference sample Substances 0.000 claims description 8
- 108010049586 Norepinephrine Plasma Membrane Transport Proteins Proteins 0.000 claims description 5
- 102000005030 SLC6A2 Human genes 0.000 claims description 5
- 210000001175 cerebrospinal fluid Anatomy 0.000 claims description 5
- 230000002093 peripheral effect Effects 0.000 claims description 5
- 102000004169 proteins and genes Human genes 0.000 claims description 5
- 210000001519 tissue Anatomy 0.000 claims description 5
- 102100021603 DNA excision repair protein ERCC-6-like 2 Human genes 0.000 claims description 4
- 101000898683 Homo sapiens DNA excision repair protein ERCC-6-like 2 Proteins 0.000 claims description 4
- YBHILYKTIRIUTE-UHFFFAOYSA-N berberine Chemical compound C1=C2CC[N+]3=CC4=C(OC)C(OC)=CC=C4C=C3C2=CC2=C1OCO2 YBHILYKTIRIUTE-UHFFFAOYSA-N 0.000 claims description 4
- 229940093265 berberine Drugs 0.000 claims description 4
- QISXPYZVZJBNDM-UHFFFAOYSA-N berberine Natural products COc1ccc2C=C3N(Cc2c1OC)C=Cc4cc5OCOc5cc34 QISXPYZVZJBNDM-UHFFFAOYSA-N 0.000 claims description 4
- 229960004170 clozapine Drugs 0.000 claims description 4
- QZUDBNBUXVUHMW-UHFFFAOYSA-N clozapine Chemical compound C1CN(C)CCN1C1=NC2=CC(Cl)=CC=C2NC2=CC=CC=C12 QZUDBNBUXVUHMW-UHFFFAOYSA-N 0.000 claims description 4
- LNEPOXFFQSENCJ-UHFFFAOYSA-N haloperidol Chemical compound C1CC(O)(C=2C=CC(Cl)=CC=2)CCN1CCCC(=O)C1=CC=C(F)C=C1 LNEPOXFFQSENCJ-UHFFFAOYSA-N 0.000 claims description 4
- 229930014626 natural product Natural products 0.000 claims description 4
- -1 240253_at Proteins 0.000 claims description 3
- WSEQXVZVJXJVFP-HXUWFJFHSA-N (R)-citalopram Chemical compound C1([C@@]2(C3=CC=C(C=C3CO2)C#N)CCCN(C)C)=CC=C(F)C=C1 WSEQXVZVJXJVFP-HXUWFJFHSA-N 0.000 claims description 2
- RTHCYVBBDHJXIQ-MRXNPFEDSA-N (R)-fluoxetine Chemical compound O([C@H](CCNC)C=1C=CC=CC=1)C1=CC=C(C(F)(F)F)C=C1 RTHCYVBBDHJXIQ-MRXNPFEDSA-N 0.000 claims description 2
- VOXZDWNPVJITMN-ZBRFXRBCSA-N 17β-estradiol Chemical compound OC1=CC=C2[C@H]3CC[C@](C)([C@H](CC4)O)[C@@H]4[C@@H]3CCC2=C1 VOXZDWNPVJITMN-ZBRFXRBCSA-N 0.000 claims description 2
- SVUOLADPCWQTTE-UHFFFAOYSA-N 1h-1,2-benzodiazepine Chemical compound N1N=CC=CC2=CC=CC=C12 SVUOLADPCWQTTE-UHFFFAOYSA-N 0.000 claims description 2
- 102100031854 60S ribosomal protein L14 Human genes 0.000 claims description 2
- CEUORZQYGODEFX-UHFFFAOYSA-N Aripirazole Chemical compound ClC1=CC=CC(N2CCN(CCCCOC=3C=C4NC(=O)CCC4=CC=3)CC2)=C1Cl CEUORZQYGODEFX-UHFFFAOYSA-N 0.000 claims description 2
- 230000007067 DNA methylation Effects 0.000 claims description 2
- 102100029674 E3 ubiquitin-protein ligase TRIM9 Human genes 0.000 claims description 2
- 102100021261 Frizzled-10 Human genes 0.000 claims description 2
- 101000704267 Homo sapiens 60S ribosomal protein L14 Proteins 0.000 claims description 2
- 101000795280 Homo sapiens E3 ubiquitin-protein ligase TRIM9 Proteins 0.000 claims description 2
- 101000819451 Homo sapiens Frizzled-10 Proteins 0.000 claims description 2
- 101000579354 Homo sapiens PHD finger protein 21A Proteins 0.000 claims description 2
- 101000780111 Homo sapiens Serine/threonine-protein phosphatase 6 regulatory ankyrin repeat subunit A Proteins 0.000 claims description 2
- YQEZLKZALYSWHR-UHFFFAOYSA-N Ketamine Chemical compound C=1C=CC=C(Cl)C=1C1(NC)CCCCC1=O YQEZLKZALYSWHR-UHFFFAOYSA-N 0.000 claims description 2
- WHXSMMKQMYFTQS-UHFFFAOYSA-N Lithium Chemical compound [Li] WHXSMMKQMYFTQS-UHFFFAOYSA-N 0.000 claims description 2
- DUGOZIWVEXMGBE-UHFFFAOYSA-N Methylphenidate Chemical compound C=1C=CC=CC=1C(C(=O)OC)C1CCCCN1 DUGOZIWVEXMGBE-UHFFFAOYSA-N 0.000 claims description 2
- UEQUQVLFIPOEMF-UHFFFAOYSA-N Mianserin Chemical compound C1C2=CC=CC=C2N2CCN(C)CC2C2=CC=CC=C21 UEQUQVLFIPOEMF-UHFFFAOYSA-N 0.000 claims description 2
- 102100028222 PHD finger protein 21A Human genes 0.000 claims description 2
- MEFKEPWMEQBLKI-AIRLBKTGSA-N S-adenosyl-L-methioninate Chemical compound O[C@@H]1[C@H](O)[C@@H](C[S+](CC[C@H](N)C([O-])=O)C)O[C@H]1N1C2=NC=NC(N)=C2N=C1 MEFKEPWMEQBLKI-AIRLBKTGSA-N 0.000 claims description 2
- 102100034285 Serine/threonine-protein phosphatase 6 regulatory ankyrin repeat subunit A Human genes 0.000 claims description 2
- YJYPHIXNFHFHND-UHFFFAOYSA-N agomelatine Chemical compound C1=CC=C(CCNC(C)=O)C2=CC(OC)=CC=C21 YJYPHIXNFHFHND-UHFFFAOYSA-N 0.000 claims description 2
- 229960002629 agomelatine Drugs 0.000 claims description 2
- 229960004372 aripiprazole Drugs 0.000 claims description 2
- 229940049706 benzodiazepine Drugs 0.000 claims description 2
- FFGPTBGBLSHEPO-UHFFFAOYSA-N carbamazepine Chemical compound C1=CC2=CC=CC=C2N(C(=O)N)C2=CC=CC=C21 FFGPTBGBLSHEPO-UHFFFAOYSA-N 0.000 claims description 2
- 229960000623 carbamazepine Drugs 0.000 claims description 2
- 229960001653 citalopram Drugs 0.000 claims description 2
- 235000014113 dietary fatty acids Nutrition 0.000 claims description 2
- 229960005309 estradiol Drugs 0.000 claims description 2
- 229930182833 estradiol Natural products 0.000 claims description 2
- 229930195729 fatty acid Natural products 0.000 claims description 2
- 239000000194 fatty acid Substances 0.000 claims description 2
- 150000004665 fatty acids Chemical class 0.000 claims description 2
- 229960002464 fluoxetine Drugs 0.000 claims description 2
- 229960003878 haloperidol Drugs 0.000 claims description 2
- BCGWQEUPMDMJNV-UHFFFAOYSA-N imipramine Chemical compound C1CC2=CC=CC=C2N(CCCN(C)C)C2=CC=CC=C21 BCGWQEUPMDMJNV-UHFFFAOYSA-N 0.000 claims description 2
- 229960004801 imipramine Drugs 0.000 claims description 2
- 229960003299 ketamine Drugs 0.000 claims description 2
- 229910052744 lithium Inorganic materials 0.000 claims description 2
- RDOIQAHITMMDAJ-UHFFFAOYSA-N loperamide Chemical compound C=1C=CC=CC=1C(C=1C=CC=CC=1)(C(=O)N(C)C)CCN(CC1)CCC1(O)C1=CC=C(Cl)C=C1 RDOIQAHITMMDAJ-UHFFFAOYSA-N 0.000 claims description 2
- 229960001571 loperamide Drugs 0.000 claims description 2
- 229960001344 methylphenidate Drugs 0.000 claims description 2
- 229960003955 mianserin Drugs 0.000 claims description 2
- 239000000203 mixture Substances 0.000 claims description 2
- 238000012986 modification Methods 0.000 claims description 2
- 230000004048 modification Effects 0.000 claims description 2
- WIQRCHMSJFFONW-UHFFFAOYSA-N norfluoxetine Chemical compound C=1C=CC=CC=1C(CCN)OC1=CC=C(C(F)(F)F)C=C1 WIQRCHMSJFFONW-UHFFFAOYSA-N 0.000 claims description 2
- 239000002773 nucleotide Substances 0.000 claims description 2
- 125000003729 nucleotide group Chemical group 0.000 claims description 2
- RZJQGNCSTQAWON-UHFFFAOYSA-N rofecoxib Chemical compound C1=CC(S(=O)(=O)C)=CC=C1C1=C(C=2C=CC=CC=2)C(=O)OC1 RZJQGNCSTQAWON-UHFFFAOYSA-N 0.000 claims description 2
- 229960000371 rofecoxib Drugs 0.000 claims description 2
- 210000003296 saliva Anatomy 0.000 claims description 2
- VGKDLMBJGBXTGI-SJCJKPOMSA-N sertraline Chemical compound C1([C@@H]2CC[C@@H](C3=CC=CC=C32)NC)=CC=C(Cl)C(Cl)=C1 VGKDLMBJGBXTGI-SJCJKPOMSA-N 0.000 claims description 2
- 229960002073 sertraline Drugs 0.000 claims description 2
- 210000002966 serum Anatomy 0.000 claims description 2
- AEQFSUDEHCCHBT-UHFFFAOYSA-M sodium valproate Chemical compound [Na+].CCCC(C([O-])=O)CCC AEQFSUDEHCCHBT-UHFFFAOYSA-M 0.000 claims description 2
- 210000002700 urine Anatomy 0.000 claims description 2
- 229940102566 valproate Drugs 0.000 claims description 2
- YQNWZWMKLDQSAC-UHFFFAOYSA-N vortioxetine Chemical compound CC1=CC(C)=CC=C1SC1=CC=CC=C1N1CCNCC1 YQNWZWMKLDQSAC-UHFFFAOYSA-N 0.000 claims description 2
- 229960002263 vortioxetine Drugs 0.000 claims description 2
- 239000000091 biomarker candidate Substances 0.000 abstract description 19
- 238000012913 prioritisation Methods 0.000 abstract description 17
- 238000003745 diagnosis Methods 0.000 abstract description 16
- 238000013459 approach Methods 0.000 abstract description 13
- 208000020016 psychiatric disease Diseases 0.000 abstract description 11
- 238000013461 design Methods 0.000 abstract description 6
- 230000004044 response Effects 0.000 abstract description 4
- 238000013517 stratification Methods 0.000 abstract description 3
- 238000012360 testing method Methods 0.000 description 40
- 238000004458 analytical method Methods 0.000 description 33
- 238000010200 validation analysis Methods 0.000 description 25
- 241000282414 Homo sapiens Species 0.000 description 14
- 238000002483 medication Methods 0.000 description 13
- CZNNQWMLAHSKRA-NVAFIHLTSA-N dTDP 1-ester with 2,6-dideoxy-L-erythro-hexopyranos-3-ulose Chemical compound C1C(=O)[C@@H](O)[C@H](C)O[C@@H]1OP(O)(=O)OP(O)(=O)OC[C@@H]1[C@@H](O)C[C@H](N2C(NC(=O)C(C)=C2)=O)O1 CZNNQWMLAHSKRA-NVAFIHLTSA-N 0.000 description 11
- 102100023441 Centromere protein J Human genes 0.000 description 10
- 101100148755 Saccharomyces cerevisiae (strain ATCC 204508 / S288c) SAS4 gene Proteins 0.000 description 9
- 230000003247 decreasing effect Effects 0.000 description 9
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 9
- 206010010144 Completed suicide Diseases 0.000 description 6
- 230000000694 effects Effects 0.000 description 6
- 230000002068 genetic effect Effects 0.000 description 6
- 238000002493 microarray Methods 0.000 description 6
- 201000010099 disease Diseases 0.000 description 5
- 238000009511 drug repositioning Methods 0.000 description 5
- 238000000692 Student's t-test Methods 0.000 description 4
- 238000010171 animal model Methods 0.000 description 4
- 208000035475 disorder Diseases 0.000 description 4
- 230000037361 pathway Effects 0.000 description 4
- 208000028173 post-traumatic stress disease Diseases 0.000 description 4
- 238000012353 t test Methods 0.000 description 4
- 230000008236 biological pathway Effects 0.000 description 3
- 210000004556 brain Anatomy 0.000 description 3
- 238000007876 drug discovery Methods 0.000 description 3
- 239000000463 material Substances 0.000 description 3
- 230000036651 mood Effects 0.000 description 3
- 239000002547 new drug Substances 0.000 description 3
- 239000002417 nutraceutical Substances 0.000 description 3
- 235000021436 nutraceutical agent Nutrition 0.000 description 3
- 230000002441 reversible effect Effects 0.000 description 3
- 201000000980 schizophrenia Diseases 0.000 description 3
- 208000024891 symptom Diseases 0.000 description 3
- 230000000007 visual effect Effects 0.000 description 3
- 102100035902 Glutamate decarboxylase 1 Human genes 0.000 description 2
- 101000873546 Homo sapiens Glutamate decarboxylase 1 Proteins 0.000 description 2
- 208000028017 Psychotic disease Diseases 0.000 description 2
- 238000000540 analysis of variance Methods 0.000 description 2
- 230000006399 behavior Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 230000003993 interaction Effects 0.000 description 2
- 239000003550 marker Substances 0.000 description 2
- 230000002265 prevention Effects 0.000 description 2
- 230000003449 preventive effect Effects 0.000 description 2
- 239000000092 prognostic biomarker Substances 0.000 description 2
- 229940124811 psychiatric drug Drugs 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 102100036321 5-hydroxytryptamine receptor 2A Human genes 0.000 description 1
- 108091005932 CCKBR Proteins 0.000 description 1
- 102000038594 Cdh1/Fizzy-related Human genes 0.000 description 1
- 108091007854 Cdh1/Fizzy-related Proteins 0.000 description 1
- 208000035473 Communicable disease Diseases 0.000 description 1
- 102100031250 Disks large-associated protein 1 Human genes 0.000 description 1
- 102100024741 Dynein light chain 2, cytoplasmic Human genes 0.000 description 1
- 102100023734 G protein-coupled receptor kinase 4 Human genes 0.000 description 1
- 102100035212 Gamma-aminobutyric acid type B receptor subunit 1 Human genes 0.000 description 1
- 102100036016 Gastrin/cholecystokinin type B receptor Human genes 0.000 description 1
- 208000019779 Gender disease Diseases 0.000 description 1
- 102100036755 Glutathione peroxidase 7 Human genes 0.000 description 1
- 102100026122 High affinity immunoglobulin gamma Fc receptor I Human genes 0.000 description 1
- 101000783617 Homo sapiens 5-hydroxytryptamine receptor 2A Proteins 0.000 description 1
- 101000756842 Homo sapiens Alpha-2A adrenergic receptor Proteins 0.000 description 1
- 101000907924 Homo sapiens Centromere protein J Proteins 0.000 description 1
- 101000844784 Homo sapiens Disks large-associated protein 1 Proteins 0.000 description 1
- 101000908706 Homo sapiens Dynein light chain 2, cytoplasmic Proteins 0.000 description 1
- 101100066427 Homo sapiens FCGR1A gene Proteins 0.000 description 1
- 101000829481 Homo sapiens G protein-coupled receptor kinase 4 Proteins 0.000 description 1
- 101001071391 Homo sapiens Glutathione peroxidase 7 Proteins 0.000 description 1
- 101001120056 Homo sapiens Phosphatidylinositol 3-kinase regulatory subunit alpha Proteins 0.000 description 1
- 101000606537 Homo sapiens Receptor-type tyrosine-protein phosphatase delta Proteins 0.000 description 1
- 101000974846 Homo sapiens Sodium/potassium-transporting ATPase subunit beta-2 Proteins 0.000 description 1
- 101000835093 Homo sapiens Transferrin receptor protein 1 Proteins 0.000 description 1
- 208000019022 Mood disease Diseases 0.000 description 1
- 102100029166 NT-3 growth factor receptor Human genes 0.000 description 1
- 208000012902 Nervous system disease Diseases 0.000 description 1
- 102000048238 Neuregulin-1 Human genes 0.000 description 1
- 108090000556 Neuregulin-1 Proteins 0.000 description 1
- 208000025966 Neurological disease Diseases 0.000 description 1
- 102100026169 Phosphatidylinositol 3-kinase regulatory subunit alpha Human genes 0.000 description 1
- 238000002123 RNA extraction Methods 0.000 description 1
- 102100039666 Receptor-type tyrosine-protein phosphatase delta Human genes 0.000 description 1
- 102000005038 SLC6A4 Human genes 0.000 description 1
- 108010012996 Serotonin Plasma Membrane Transport Proteins Proteins 0.000 description 1
- 102100022791 Sodium/potassium-transporting ATPase subunit beta-2 Human genes 0.000 description 1
- 102100026144 Transferrin receptor protein 1 Human genes 0.000 description 1
- 230000001154 acute effect Effects 0.000 description 1
- 239000000654 additive Substances 0.000 description 1
- 230000000996 additive effect Effects 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 230000007815 allergy Effects 0.000 description 1
- 238000009640 blood culture Methods 0.000 description 1
- 230000011496 cAMP-mediated signaling Effects 0.000 description 1
- 230000008777 canonical pathway Effects 0.000 description 1
- 210000004027 cell Anatomy 0.000 description 1
- 238000004113 cell culture Methods 0.000 description 1
- 239000003153 chemical reaction reagent Substances 0.000 description 1
- 238000009223 counseling Methods 0.000 description 1
- 238000002790 cross-validation Methods 0.000 description 1
- 238000013499 data model Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 229940011871 estrogen Drugs 0.000 description 1
- 239000000262 estrogen Substances 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000013401 experimental design Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000010234 longitudinal analysis Methods 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000006996 mental state Effects 0.000 description 1
- 238000010208 microarray analysis Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000003012 network analysis Methods 0.000 description 1
- 238000010855 neuropsychological testing Methods 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 229940012843 omega-3 fatty acid Drugs 0.000 description 1
- 235000020660 omega-3 fatty acid Nutrition 0.000 description 1
- 239000006014 omega-3 oil Substances 0.000 description 1
- 238000001543 one-way ANOVA Methods 0.000 description 1
- 238000003068 pathway analysis Methods 0.000 description 1
- 230000003234 polygenic effect Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000010076 replication Effects 0.000 description 1
- 238000012552 review Methods 0.000 description 1
- 238000013515 script Methods 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 108010090953 subunit 1 GABA type B receptor Proteins 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- 108010064892 trkC Receptor Proteins 0.000 description 1
- 238000009424 underpinning Methods 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61K—PREPARATIONS FOR MEDICAL, DENTAL OR TOILETRY PURPOSES
- A61K31/00—Medicinal preparations containing organic active ingredients
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61P—SPECIFIC THERAPEUTIC ACTIVITY OF CHEMICAL COMPOUNDS OR MEDICINAL PREPARATIONS
- A61P25/00—Drugs for disorders of the nervous system
- A61P25/22—Anxiolytics
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Q—MEASURING 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/00—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
- C12Q1/68—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H20/00—ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance
- G16H20/10—ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance relating to drugs or medications, e.g. for ensuring correct administration to patients
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H20/00—ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance
- G16H20/70—ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance relating to mental therapies, e.g. psychological therapy or autogenous training
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
- G16H50/20—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for computer-aided diagnosis, e.g. based on medical expert systems
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H50/00—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
- G16H50/30—ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for calculating health indices; for individual health risk assessment
Definitions
- Anxiety disorders are highly prevalent and insidious in their effects on one’s ability to do things and quality of life.
- Psychiatric patients may have increased anxiety, as well as increased reasons for anxiety, due to their often, adverse life trajectory. Improvements are needed to adequately diagnose and treat all individuals suffering with anxiety disorders.
- biomarkers for anxiety provide a means of assessing state severity, short-term risk, and long-term risk.
- the biomarkers can also be used for drug repurposing.
- Treating involves administered at least one course of treatment, treatment may include psychiatric counseling, administering certain physical intervention, and/or prescribing and/or administering at least one therapeutic compound. Treating may include at least one of the following outcomes, curing, mitigating, managing or otherwise recuing the severity and/or the number of frequency anxious behaviors, especially behaviors that negatively impact a individual’s quality of life.
- a first embodiment is a therapeutic compound selected from one or more compounds in Tables 4 and Table 6, for use in a method of mitigating and anxiety disorder in a subject in need thereof for example wherein the therapy is a drug, a natural compound, and combinations thereof.
- a second embodiment is the therapeutic compound for use according to the first embodiment, wherein the drug is selected from the group consisting of: estradiol, loperamide, omage-3 fatty acids, lithium, fluoxetine, norfluoxetine, imipramine, citalopram, carbamazepine,
- a third embodiment is the therapeutic compound for use according to the first embodiment, wherein the subject has an anxiety disorder.
- a fourth embodiment is a therapeutic compound for use according to the first embodiment, wherein the subject is a male subject or a female subject.
- a fifth embodiment is a computer-implemented method for assessing the level of anxiety in a subject, and for assessing risk of developing anxiety disorder in the future in a subject, the method comprising: computing a score based on RNA level, protein level, DNA methylation, a single nucleotide polymorphism, for a panel of biomarkers chosen from Table 2, Table 5, and Figure 2, in a sample obtained from a subject; computing a score based on a reference expression level of the panel of biomarkers; and identifying a difference between the score in the sample obtained from the subject and the score in the reference sample, wherein the difference in the score in the sample obtained from the subject and the score in the reference sample indicates an anxiety disorder, and a risk of future anxiety disorders.
- a sixth embodiment is the method according to the fifth embodiment, wherein the subject is: a male or a female and the biomarker is at least one biomarker selected from the group consisting of: ERCC6L2, SLC6A2, FZD10, and RPL14; or a female and the biomarker is at least one biomarker selected from the group consisting of: PHF21A, SLC6A2, 240253_at, and TRIM9; or a male and the biomarker is at least one biomarker selected from the group consisting of: ERCC6L2, ANKRD28, ACPI, and RP
- a seventh embodiment is the method according to the fifth or the sixth embodiments, further including the steps of: identifying a difference between the score in the sample obtained from the subject and the score in the reference sample, administering treatment to the subject, wherein the treatment reduces the difference between the score in the sample from the subject and the score in the reference sample, and mitigates anxiety in the subject, or the risk for developing anxiety disorder in the subject.
- An eighth embodiment is the method according to the seventh embodiment, further including, the step of: measuring a change in the score after treatment.
- a ninth embodiment is the method according to the eighth embodiment, wherein the treatment is at least one treatment selected from the group consisting of: lifestyle modification,
- SUBSTITUTE SHEET (RULE 26) electro-magnetic intervention, and administering a therapeutic compound.
- a tenth embodiment is the method according to the ninth embodiment, wherein the therapy is selected by a computer- implemented method selected from the group consisting of one or more psychiatric compounds from Table 6, and wherein each therapy selection is based on a panel of one or more individual biomarkers.
- An eleventh embodiment is the method according to the tenth embodiment, wherein the therapy is selected based on a panel of individual biomarkers changed in a subject, by a computer-implemented method for therapy selection, and comprises, of one or more compounds in Table 4.
- An twelfth embodiment is a method for assessing and mitigating anxiety disorder, comprising determining an expression level of a panel of biomarkers listed in Table 2, Table 5 and Figure 2, in a sample for example wherein the sample comprises a peripheral tissue, blood, saliva, cerebrospinal fluid (CSF), serum, urine, or stool; wherein the expression level of the biomarkers in the sample is different relative to a reference expression level; and includes the steps of: identifying the subject currently having or at risk of having an anxiety disorder based on a biomarker panel score relative to a biomarker panel score of a reference; selecting a therapy based on the score from the biomarkers, wherein the therapy includes; and administering one or more compounds selected from the group consisting of compounds selected from the groups consisting of compounds selected from Tables 4 and Table 6.
- CSF cerebrospinal fluid
- a thirteenth embodiment is a composition comprising one or more compounds selected from Tables 4 and Table 6, for use in a method for treating anxiety disorder.
- FIG. 1 A Discovery cohort, within subject changes in anxiety state (N of 58).
- FIG. 1 B Discovery cohort 141 male and 17 female psychiatric participants have at least one switch between a low anxiety state visit and a high anxiety state visit.
- FIG. 1 C Convergent functional genomics. Multiple independent lines of evidence for identification and prioritization of anxiety biomarkers.
- FIG. 1 E Graphic representation of Discovery, Prioritization and Validation of biomarkers.
- FIG. 2A Graph of predictions for state- high anxiety state (SAS4>60). See also TABLE 6.
- FIG. 2B Graph of predictions for state- high STAI (STA>55). Light colored bars, cross sectional; dark colored bars, longitudinal.
- FIG. 2C Graph of anxiety trait predictions first year. Light colored bars, cross sectional; dark colored bars, longitudinal.
- FIG. 2D Graph of anxiety trait predictions-all future year. Light colored bars, cross sectional; dark colored bars, longitudinal.
- FIG. 5 SASA4 scaler for measuring anxiety State (Niculescu, et al., 006, 2015). The score is an average of 4 items, Ranges from 0-100.
- FIG. 6 Graphical depiction of the correlation between SAS4 and STAI State.
- Our SAS4 scale shows a moderate to strong correlation with a well-validated clinical measure of anxiety, the STAI State, while being more temporally related to that particular moment in time, quantitative, and simple.
- SUBSTITUTE SHEET (RULE 26) worsening (hospitalizations with anxiety as the primary cause), in another independent cohort of subjects.
- the present disclosure is generally directed at methods for assessing anxiety disorders and early identification of risk for future episodes of acute or chronic anxiety, as well as methods for matching patients and drugs for prevention and mitigation of anxiety and for monitoring response to treatment.
- the methods may further include the generation of a report providing a risk score and/or personalized treatment options.
- the present disclosure generally is directed to drugs for mitigating anxiety in subjects.
- Particular drugs have been found that can mitigate anxiety in subjects universally; that is, drugs that can be used for mitigating anxiety across genders an in the presence of confounding factors such as other psychiatric disorders.
- Some drugs, however, have been found that can be used more effectively for mitigating anxiety dependent on gender, confounding factors, and combinations thereof.
- the present disclosure is directed to blood gene expression biomarkers that are more universal in nature; that is, blood biomarkers that can be used for predicting anxiety disorders in one of more genders. Accordingly, a longitudinal within- participant design and large cohorts were used. Additionally, subtypes of anxiety were identified based on mental state at the time that anxiety was manifest.
- FIG. 1 Steps 1-3: Discovery, Prioritization and Validation of Biomarkers for Anxiety (A.) Cohorts used in study, depicting flow of discovery, prioritization, validation of biomarkers from each step and independent testing cohorts. (B.) Discovery cohort longitudinal within- subject analysis. Phchp### is study ID for each subject. V# denotes visit number.
- C. Convergent Functional Genomics lines of evidence from human to non-human.
- D In the validation step biomarkers are assessed for stepwise change from the validation group with severe Anxiety, to the discovery groups of subjects with high Anxiety, low Anxiety, to the validation group with severe Anxiety, using ANOVA.
- N number of testing visits.
- the histograms depict a top increase and a top decrease
- Subjects were recruited from the patient population at the Indianapolis VA Medical Center. All subjects understood and signed informed consent forms detailing the research goals, procedure, caveats and safeguards, per IRB approved protocol. Subjects completed diagnostic assessments and extensive structured neuropsychological testing at each testing visit, 3-6 months apart or whenever a new psychiatric
- SUBSTITUTE SHEET (RULE 26) hospitalization occurred. At each testing visit, they received a series of rating scales, including a self-report visual analog scale (1-100) for quantitatively assessing state anxiety at that particular moment in time (Simplified Anxiety Scale- SAS-4). This 4- item scale looks at anxiety overall, as well as fear, anger, and uncertainty. Each of the items are a VAS of 0 to 100, related to that moment in time. As such, it generates temporal, quantitative, and targeted data.
- the prioritization step that occurs after discovery is based on a field-wide convergence with literature that includes genetic data and animal model data, that are unrelated to medication effects.
- the discovery, validation and replication by testing in independent cohorts of the biomarkers, with our design occurs despite the subjects having different genders, diagnoses, being on various different medications, and other lifestyle variables.
- RNA extraction Whole blood (2.5 ml) was collected into each PaxGene tube by routine venipuncture. PaxGene tubes contain proprietary reagents for the stabilization of RNA. Total RNA was extracted and processed as previously described.
- Microarrays Microarray work was carried out using previously described methodology 1 .
- probesets that have a 1.2-fold change are then assigned either a 1 (increased in high anxiety) or a -1 (decreased in high anxiety) in each comparison. Fold changes between 1.1 and 1.2 are given 0.5, and fold changes less than 1.1 are given 0. These values were then summed for each probeset across all the comparisons and subjects, yielding a range of raw scores.
- the probesets above the 33.3% of raw scores were carried forward in analyses (Figure 1), and received an internal score of 2 points; those above 50% 4 points, and those above 80% 6 points 2 3 4 . We have developed in our labs R scripts to automate and conduct all these large dataset analyses in bulk, checked against human manual scoring 4 .
- Step 2 Prioritization using Convergent Functional Genomics (CFG) Databases.
- the AP derived and DE derived lists of genes were combined, and the gene expression data corresponding to them was used for the validation analysis.
- the 3 groups (low anxiety, high anxiety, clinical severe anxiety) were assembled out of Affymetrix. cell data that was RMA normalized by gender and diagnosis. We transferred the log transformed expression data to an Excel sheet, and non-log transformed the data by taking 2 to the power of the
- IPA Ingenuity Pathway Analysis, version 24390178, Qiagen
- David Functional Annotation Bioinformatics Microarray Analysis National Institute of Allergy and Infectious Diseases
- KEGG Kyoto Encyclopedia of Genes and Genomes
- CFG beyond Anxiety evidence for involvement in other psychiatric and related disorders.
- SUBSTITUTE SHEET (RULE 26) or STAI), and predict future risk of anxiety (trait -future hospitalizations with anxiety).
- biomarkers were combined by simple summation of the increased risk biomarkers minus the decreased risk biomarkers. Predictions were performed using R-studio. For cross-sectional analyses, we used biomarker expression levels, z-scored by gender and diagnosis. For longitudinal analyses, we combined four measures: biomarker expression levels, slope (defined as ratio of levels at current testing visit vs. previous visit, divided by time between visits), maximum levels (at any of the current or past visits), and maximum slope (between any adjacent current or past visits). For decreased biomarkers, we used the minimum rather than the maximum for level calculations. All four measures were Z-scored, then combined in an additive fashion into a single measure. The longitudinal analysis was carried out in a sub-cohort of the testing cohort consisting of subjects that had at least two visits (timepoints).
- Receiver-operating characteristic (ROC) analyses between marker levels and anxiety state were performed by assigning subjects visits with a anxiety SAS-4 score of >60 into the high anxiety category, and subjects with STAI scores >55 in the high anxiety category.
- ROC Receiver-operating characteristic
- SUBSTITUTE SHEET (RULE 26) (2 D) are at least nominally significant. Some gender and diagnosis group are missing from the graph as they did not have any significant biomarkers or that the cohort was too small with limited data for the z-scoring by gender-dx.
- Cross-sectional is based on levels at one visit. Longitudinal is based on levels at multiple visits (integrates levels at most recent visit, maximum levels, slope into most recent visit, and maximum slope). Dividing lines represent the cutoffs for a test performing at chance levels (white), and at the same level as the best biomarkers for all subjects in cross-sectional (gray) and longitudinal (black) based predictions. Biomarkers perform better than chance. Biomarkers performed better when personalized by gender and diagnosis. * survived Bonferroni correction for the number of candidate biomarkers tested. Predicting Trait- Future Psychiatric Hospitalization with Anxiety as a Symptom/Reason for Admission.
- a Cox regression was performed using the time in days from the testing visit date to first hospitalization date in the case of patients who had been hospitalized, or 365 days for those who did not.
- the odds ratio was calculated such that a value greater than 1 always indicates increased risk for hospitalization, regardless if the biomarker is increased or decreased in expression.
- SUBSTITUTE SHEET (RULE 26) geographically and/or be lost to follow-up.
- the Cox regression was performed using the time in days from visit date to first hospitalization date in the case of patients who had hospitalizations with anxiety, or from visit date to last note date in the electronic medical records for those who did not.
- the first average should be higher than the second average in increased biomarkers, and the reverse is true for decreased biomarkers. 15 out of 19 biomarkers were thus concordant.
- the digitized biomarkers were also used for matching with existing psychiatric medications and alternative treatments (nutraceuticals and others). We used our large datasets and literature databases to match biomarkers to medications that had effects on gene expression opposite to their expression in high anxiety. Each medication matched to a biomarker got the biomarker score of 1, 0.5 or 0. The scores for the medications were added, normalized for the number of biomarkers that were 1 or 0.5 in that patient, resulting in a percentile match. Thus, psychiatric medications matched to the patient and ranked in order of impact on the panel.
- Step 1 Discovery we identified candidate blood gene expression biomarkers that: 1. change in expression in blood between self-reported low anxiety and high anxiety states, 2. track the anxiety state across visits in a subject, and 3. track the anxiety state in multiple subjects.
- SAS-4 visual analog measure for anxiety state
- the SAS-4 quantitates anxiety state at a particular moment in time, and normalizes anxiety
- Step 2 Prioritization we used a Convergent Functional Genomics (CFG) approach to prioritize the candidate biomarkers identified in the discovery step (33% cutoff, internal score of >2pt.) by using published literature evidence (genetic, gene expression and proteomic), from human and animal model studies, for involvement in anxiety disorders (Figure 1 and Table 2). For Axniety there were 284 probesets (corresponding to 238 unique genes) that had a total score (combined discovery score and prioritization CFG score) of 6 and above. These were carried forward to the validation step. This represents approximately a 10-fold enrichment of the probesets on the Affymetrix array.
- CFG Convergent Functional Genomics
- HTR2A is at the overlap of a network containing GABBR1, GAD1 and SLC6A4, and one centered on CDH1 that also contains PIK3R1 and IGFR1.
- a third network includes DLGAP1, PTPRD and DYNLL2. These networks may have biological significance and could be targeted therapeutically.
- the gene expression data in the test cohorts was normalized (Z-scored) across genders and various psychiatric diagnoses, before those different demographic groups were combined.
- CFE4 convergent functional evidence
- SUBSTITUTE SHEET (RULE 26) first year hospitalization with anxiety, trait all future hospitalizations with anxiety- up to 3 points each if it significantly predicts in all subjects, 2 points if in gender, 1 points if in gender/diagnosis).
- the total score can be up to 36 points: 24 from our empirical data, and 12 from literature data used for CFG. We weigh our new empirical data more than the literature data, as it is functionally related to mood in 3 independent cohorts (discovery, validation, testing).
- the goal is to highlight, based on the totality of our data and of the evidence in the field to date, biomarkers that have all around evidence: track anxiety, have convergent evidence for involvement in anxiety disorders, and predict anxiety state and future clinical events (Table 2). Such biomarkers merit priority evaluation in future clinical trials.
- Bioinformatic analyses using the gene expression signature of the panel of top biomarkers for anxiety identified new potential therapeutics for anxiety, such as estrogen.
- SUBSTITUTE SHEET (RULE 26) Table 1. Aggregate Demographics, BP-bipolar, MDD- depression, MOOD- mood nos, SZ- schizophrenia, SZA- schizoaffective, PSYCH- psychosis nos, PTSD-post-traumatic stress disorder.
- CMAP from genes ADRA2A, ATP1B2, CCKBR, FZD1O,GAD1,GPX7,GRK4 ,NRG1,NTRK3,SLC6A2,TFRC). Drugs that have opposite gene expression effects to the gene expression signature.
- SUBSTITUTE SHEET (RULE 26) Table 6. Pharmacogenomics. Top biomarkers (from Table 2) that are targets of existing drugs and are changed in expression in opposite direction to high anxiety. (I)- increased in expression, (D)- decreased in expression. DE- differential expression, AP-Absent/Present.
- Biomarkers may be useful for matching patients to medications and measuring response to treatment (pharmacogenomics) ( Figure 4, Tables 2 and 6), as well as new drug discovery and repurposing (Table 4).
- SUBSTITUTE SHEET (RULE 26) anxiety disorders episode in their lifetime, that anxiety disorders can severely affect quality of life, sometimes leading to suicides, and that not all patients respond to current treatments or the treatments are addictive, the need for and importance of efforts such as ours cannot be overstated.
Landscapes
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Public Health (AREA)
- Chemical & Material Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Epidemiology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Primary Health Care (AREA)
- Organic Chemistry (AREA)
- Biomedical Technology (AREA)
- Medicinal Chemistry (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Pharmacology & Pharmacy (AREA)
- Proteomics, Peptides & Aminoacids (AREA)
- Databases & Information Systems (AREA)
- Data Mining & Analysis (AREA)
- Veterinary Medicine (AREA)
- Pathology (AREA)
- Animal Behavior & Ethology (AREA)
- Wood Science & Technology (AREA)
- Zoology (AREA)
- Biochemistry (AREA)
- Molecular Biology (AREA)
- General Engineering & Computer Science (AREA)
- Genetics & Genomics (AREA)
- Microbiology (AREA)
- Immunology (AREA)
- Biotechnology (AREA)
- Biophysics (AREA)
- Analytical Chemistry (AREA)
- Physics & Mathematics (AREA)
- Neurology (AREA)
- Neurosurgery (AREA)
- Chemical Kinetics & Catalysis (AREA)
- General Chemical & Material Sciences (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Child & Adolescent Psychology (AREA)
- Developmental Disabilities (AREA)
- Hospice & Palliative Care (AREA)
Abstract
We used a powerful longitudinal within-subject design in individuals with psychiatric disorders to discover blood gene expression change between self-reported low anxiety and high anxiety states. We prioritized the list of candidate biomarkers with a Convergent Functional Genomics approach using other evidence in the field and validated our top biomarkers from discovery and prioritization in an independe cohort of psychiatric subjects with clinically severe anxiety. We then tested if these candidate biomarkers are able to predict anxiety severity state, and future clinical worsening (hospitalizations with anxiety as a contributory cause), in another independent cohort of psychiatric subjects. We showed increased accuracy of individual biomarkers with a personalized approach, by gender and diagnosis, particularly in women. Finally, we identified which of our biomarkers are targets of existing drugs and thus can be used for pharmacogenomic population stratification and measuring of response to treatment, as well as used biomarker gene expression signatures to identify drugs that can be repurposed for treating anxiety.
Description
Precision Medicine for Anxiety Disorders: Objective Assessment, Risk Prediction, Pharmacogenomics, and Repurposed Drugs
STATEMENT OF GOVERNMENT SUPPORT
[0001] This invention was made with government support under MH117431 awarded by the National Institutes of Health and CX000139 merit award by the Veterans Administration. The government has certain rights in the invention.
BACKGROUND
[0002] Anxiety disorders are highly prevalent and insidious in their effects on one’s ability to do things and quality of life. Psychiatric patients may have increased anxiety, as well as increased reasons for anxiety, due to their often, adverse life trajectory. Improvements are needed to adequately diagnose and treat all individuals suffering with anxiety disorders.
SUMMARY
[0003] Provided here are newly identified blood gene expression biomarkers for anxiety. The biomarkers provide a means of assessing state severity, short-term risk, and long-term risk. The biomarkers can also be used for drug repurposing.
[0004] Some aspects of the invention include methods for treating an individual experiencing or at a heighted risk for developing symptoms is anxiety disorders. Treating involves administered at least one course of treatment, treatment may include psychiatric counseling, administering certain physical intervention, and/or prescribing and/or administering at least one therapeutic compound. Treating may include at least one of the following outcomes, curing, mitigating, managing or otherwise recuing the severity and/or the number of frequency anxious behaviors, especially behaviors that negatively impact a individual’s quality of life.
[0005] A first embodiment is a therapeutic compound selected from one or more compounds in Tables 4 and Table 6, for use in a method of mitigating and anxiety disorder in a subject in need thereof for example wherein the therapy is a drug, a natural compound, and combinations thereof.
[0006] A second embodiment is the therapeutic compound for use according to the first embodiment, wherein the drug is selected from the group consisting of: estradiol, loperamide, omage-3 fatty acids, lithium, fluoxetine, norfluoxetine, imipramine, citalopram, carbamazepine,
1
SUBSTITUTE SHEET (RULE 26)
aripiprazole, haloperidol, mianserin, SAM, vortioxetine, agomelatine, sertraline, benzodiazepine, clozapine, ketamine, valproate, clozapine, rofecoxib, berberine, methylphenidate, or berberine. [0007] A third embodiment is the therapeutic compound for use according to the first embodiment, wherein the subject has an anxiety disorder.
[0008] A fourth embodiment is a therapeutic compound for use according to the first embodiment, wherein the subject is a male subject or a female subject.
[0009] A fifth embodiment is a computer-implemented method for assessing the level of anxiety in a subject, and for assessing risk of developing anxiety disorder in the future in a subject, the method comprising: computing a score based on RNA level, protein level, DNA methylation, a single nucleotide polymorphism, for a panel of biomarkers chosen from Table 2, Table 5, and Figure 2, in a sample obtained from a subject; computing a score based on a reference expression level of the panel of biomarkers; and identifying a difference between the score in the sample obtained from the subject and the score in the reference sample, wherein the difference in the score in the sample obtained from the subject and the score in the reference sample indicates an anxiety disorder, and a risk of future anxiety disorders.
[0010] A sixth embodiment is the method according to the fifth embodiment, wherein the subject is: a male or a female and the biomarker is at least one biomarker selected from the group consisting of: ERCC6L2, SLC6A2, FZD10, and RPL14; or a female and the biomarker is at least one biomarker selected from the group consisting of: PHF21A, SLC6A2, 240253_at, and TRIM9; or a male and the biomarker is at least one biomarker selected from the group consisting of: ERCC6L2, ANKRD28, ACPI, and RP
[0011] A seventh embodiment is the method according to the fifth or the sixth embodiments, further including the steps of: identifying a difference between the score in the sample obtained from the subject and the score in the reference sample, administering treatment to the subject, wherein the treatment reduces the difference between the score in the sample from the subject and the score in the reference sample, and mitigates anxiety in the subject, or the risk for developing anxiety disorder in the subject.
[0012] An eighth embodiment is the method according to the seventh embodiment, further including, the step of: measuring a change in the score after treatment.
[0013] A ninth embodiment is the method according to the eighth embodiment, wherein the treatment is at least one treatment selected from the group consisting of: lifestyle modification,
2
SUBSTITUTE SHEET (RULE 26)
electro-magnetic intervention, and administering a therapeutic compound.
[0014] A tenth embodiment is the method according to the ninth embodiment, wherein the therapy is selected by a computer- implemented method selected from the group consisting of one or more psychiatric compounds from Table 6, and wherein each therapy selection is based on a panel of one or more individual biomarkers.
[0015] An eleventh embodiment is the method according to the tenth embodiment, wherein the therapy is selected based on a panel of individual biomarkers changed in a subject, by a computer-implemented method for therapy selection, and comprises, of one or more compounds in Table 4.
[0016] An twelfth embodiment is a method for assessing and mitigating anxiety disorder, comprising determining an expression level of a panel of biomarkers listed in Table 2, Table 5 and Figure 2, in a sample for example wherein the sample comprises a peripheral tissue, blood, saliva, cerebrospinal fluid (CSF), serum, urine, or stool; wherein the expression level of the biomarkers in the sample is different relative to a reference expression level; and includes the steps of: identifying the subject currently having or at risk of having an anxiety disorder based on a biomarker panel score relative to a biomarker panel score of a reference; selecting a therapy based on the score from the biomarkers, wherein the therapy includes; and administering one or more compounds selected from the group consisting of compounds selected from the groups consisting of compounds selected from Tables 4 and Table 6.
[0017] A thirteenth embodiment is a composition comprising one or more compounds selected from Tables 4 and Table 6, for use in a method for treating anxiety disorder.
BRIEF DESCRIPTION OF THE FIGURES
[0018] FIG. 1 A. Discovery cohort, within subject changes in anxiety state (N of 58).
[0019] FIG. 1 B. Discovery cohort 141 male and 17 female psychiatric participants have at least one switch between a low anxiety state visit and a high anxiety state visit.
[0020] FIG. 1 C. Convergent functional genomics. Multiple independent lines of evidence for identification and prioritization of anxiety biomarkers.
[0021] FIG. 1 D. Graphic presentation of data. Left panel, EFNAS, P=2.02 E-04.; right panel TG51, P=2.3E-03.
[0022] FIG. 1 E. Graphic representation of Discovery, Prioritization and Validation of biomarkers.
3
SUBSTITUTE SHEET (RULE 26)
[0023] FIG. 2A. Graph of predictions for state- high anxiety state (SAS4>60). See also TABLE 6.
[0024] FIG. 2B. Graph of predictions for state- high STAI (STA>55). Light colored bars, cross sectional; dark colored bars, longitudinal.
[0025] FIG. 2C. Graph of anxiety trait predictions first year. Light colored bars, cross sectional; dark colored bars, longitudinal.
[0026] FIG. 2D. Graph of anxiety trait predictions-all future year. Light colored bars, cross sectional; dark colored bars, longitudinal.
[0027] FIG. 3. String network analysis. For top candidate biomarkers (n=82 genes, 95 probesets).
[0028] FIG. 4. Summary of an analysis of Phchp328, a specific Female subject, MDD, 37 years old, Caucasian; SAS4=92.5; STAI State=68/80. Anxiety Score 86.67% (13/15) (High), Chronic Anxiety Risk: 88.89% (High).
[0029] FIG. 5. SASA4 scaler for measuring anxiety State (Niculescu, et al., 006, 2015). The score is an average of 4 items, Ranges from 0-100.
[0030] FIG. 6. Graphical depiction of the correlation between SAS4 and STAI State. Our SAS4 scale shows a moderate to strong correlation with a well-validated clinical measure of anxiety, the STAI State, while being more temporally related to that particular moment in time, quantitative, and simple. Pearson Correlation Coefficient 0.667503766; sample size 784; test statistic 25.06861378; p-value 1.4989E-102.
DETAILED DESCRIPTION
[0031] We endeavored to find objective blood biomarkers for anxiety disorders. First, we used a powerful longitudinal within-subject design in individuals with psychiatric disorders to discover blood gene expression changes between self-reported low anxiety and high anxiety states, as measured by a visual analog scale we developed, the Simplified Anxiety Scale (SAS- 4). Second, we prioritized the list of candidate biomarkers with a Bayesian-like Convergent Functional Genomics approach, comprehensively integrating previous human and animal model evidence in the field. Third, we validated our top biomarkers from discovery and prioritization in an independent cohort of psychiatric subjects with clinically severe anxiety. We prioritized a list of 95 candidate biomarkers that had the most evidence from the first three steps. Fourth, we tested if these candidate biomarkers are able to predict anxiety severity state, and future clinical
4
SUBSTITUTE SHEET (RULE 26)
worsening (hospitalizations with anxiety as the primary cause), in another independent cohort of subjects. We tested the biomarkers in all subjects in the test cohort, as well as in a more personalized fashion by gender and anxiety disorder diagnosis, showing increased accuracy with the personalized approach, particularly in women. Fifth, we analyzed the biological pathways and networks the biomarkers are involved in, as well as which of our top biomarkers have evidence for involvement in other psychiatric and related disorders. Sixth, we identified which of our biomarkers are targets of existing drugs and thus can be used for pharmacogenomic population stratification and measuring of response to treatment, as well as used the top biomarkers gene expression signature to identify drugs and natural compounds that can be repurposed for treating anxiety.
[0032] Unless defined otherwise, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which the disclosure belongs. Although any methods and materials similar to or equivalent to those described herein can be used in the practice or testing of the present disclosure, the preferred methods and materials are described below.
[0033] The present disclosure is generally directed at methods for assessing anxiety disorders and early identification of risk for future episodes of acute or chronic anxiety, as well as methods for matching patients and drugs for prevention and mitigation of anxiety and for monitoring response to treatment. The methods may further include the generation of a report providing a risk score and/or personalized treatment options. Further, the present disclosure generally is directed to drugs for mitigating anxiety in subjects. Particular drugs have been found that can mitigate anxiety in subjects universally; that is, drugs that can be used for mitigating anxiety across genders an in the presence of confounding factors such as other psychiatric disorders. Some drugs, however, have been found that can be used more effectively for mitigating anxiety dependent on gender, confounding factors, and combinations thereof.
[0034] In additional embodiments, the present disclosure is directed to blood gene expression biomarkers that are more universal in nature; that is, blood biomarkers that can be used for predicting anxiety disorders in one of more genders. Accordingly, a longitudinal within- participant design and large cohorts were used. Additionally, subtypes of anxiety were identified based on mental state at the time that anxiety was manifest.
5
SUBSTITUTE SHEET (RULE 26)
Materials and Methods
Cohorts
[0035] We used three independent cohorts: discovery (major psychiatric disorders with changes in state anxiety), validation (major psychiatric disorders with clinically severe anxiety), and testing (an independent major psychiatric disorders cohort for predicting state anxiety, and for predicting trait anxiety (future hospitalization with anxiety as the primary reason) (Figure 1). [0036] Referring now to FIG, 1. Figure 1. Steps 1-3: Discovery, Prioritization and Validation of Biomarkers for Anxiety (A.) Cohorts used in study, depicting flow of discovery, prioritization, validation of biomarkers from each step and independent testing cohorts. (B.) Discovery cohort longitudinal within- subject analysis. Phchp### is study ID for each subject. V# denotes visit number. C.) Convergent Functional Genomics lines of evidence from human to non-human. (D) In the validation step biomarkers are assessed for stepwise change from the validation group with severe Anxiety, to the discovery groups of subjects with high Anxiety, low Anxiety, to the validation group with severe Anxiety, using ANOVA. N=number of testing visits. The histograms depict a top increase and a top decrease Validator biomarkers (E) Pyramid schemes and scoring at each of the steps. Discovery probesets are identified based on their score for tracking Mood with a maximum of internal points of 6 (33% (2pt), 50% (4pt) and 80% (6pt)). Prioritization with CFG for prior evidence of involvement in mood disorders. In the prioritization step probesets are converted to their associated genes using Affymetrix annotation and GeneCards. Genes are prioritized and scored using CFG for Anxiety evidence with a maximum of 12 external points. Genes scoring at least 6 points out of a maximum possible of 18 total internal and external scores points are carried to the validation step. Validation in an independent cohort of psychiatric patients with clinically severe Anxiety (STAI State >55 and SAS4 >=60). 4 biomarkers were nominally significant, and 57 biomarkers were stepwise changed. We selected the top CFE3 score >8 (n= 95) for further testing in independent cohorts. [0037] Similar to our previous studies, the psychiatric subjects are part of a larger longitudinal cohort of adults that we are continuously collecting. Subjects were recruited from the patient population at the Indianapolis VA Medical Center. All subjects understood and signed informed consent forms detailing the research goals, procedure, caveats and safeguards, per IRB approved protocol. Subjects completed diagnostic assessments and extensive structured neuropsychological testing at each testing visit, 3-6 months apart or whenever a new psychiatric
6
SUBSTITUTE SHEET (RULE 26)
hospitalization occurred. At each testing visit, they received a series of rating scales, including a self-report visual analog scale (1-100) for quantitatively assessing state anxiety at that particular moment in time (Simplified Anxiety Scale- SAS-4). This 4- item scale looks at anxiety overall, as well as fear, anger, and uncertainty. Each of the items are a VAS of 0 to 100, related to that moment in time. As such, it generates temporal, quantitative, and targeted data.
[0038] At each testing visit we collected whole blood (10 ml) in two RNA-stabilizing PAXgene tubes, labeled with an anonymized ID number, and stored at -80 degrees C in a locked freezer until the time of future processing. Whole-blood RNA was extracted for microarray gene expression studies from the PAXgene tubes, as detailed below.
[0039] For this study, our within-subject discovery cohort, from which the biomarker data were derived, consisted of 58 subjects (41 males, 17 females) with multiple testing visits, who each had at least one diametric change in anxiety state from low anxiety state (SAS-4 score of < 40/100) to a high anxiety state (SAS-4 score of >60/100), or vice versa, from one testing visit to another (Figures 1). There were 2 subjects with 5 visits each, 3 subjects with 4 visits each, 21 subjects with 3 visits each, and 32 subjects with 2 visits each resulting in a total of 149 blood samples for subsequent gene expression microarray studies (Figure 1, Table 1).
[0040] Our independent validation cohort, in which the top biomarker findings were validated for being even more changed in expression, consisted of 40 subjects (32 male and 8 female) with clinically severe anxiety (SAS-4 scores >60, and concordant high anxiety STAI State scores > 55 (Table 1).
[0041] We use independent test cohorts for predicting high anxiety state (SAS-4 <60) consisted of 161 male and 36 female subjects and another test cohorts, using (STAI > 55) consisted of 159 male and 36 female subjects with psychiatric disorders, demographically matched with the discovery cohort, with one or multiple testing visit in our study, with either low anxiety, intermediate anxiety, or high anxiety states based on the two state scales (Figure 1 and Table 1).
[0042] For trait predictions of future hospitalizations with anxiety as the primary reason (Figure 1 and Table 1), are a subset of the independent test cohort for which we had longitudinal follow-up with electronic medical records. The subjects’ subsequent number of hospitalizations with anxiety was tabulated from electronic medical records.
Medications.
7
SUBSTITUTE SHEET (RULE 26)
[0043] The subjects in the discovery cohort were all diagnosed with various psychiatric disorders (Table 1) and had various medical co-morbidities. Their medications were listed in their electronic medical records and documented by us at the time of each testing visit. Medications can have a strong influence on gene expression. However, there was no consistent pattern of any particular type of medication, as our subjects were on a wide variety of different medications, psychiatric and non-psychiatric. Furthermore, the independent validation and testing cohort’s gene expression data was Z-scored by gender and by diagnosis before being combined, to normalize for any such effects. Some subjects may be non-compliant with their treatment and may thus have changes in medications or drug of abuse not reflected in their medical records. That being said, our goal is to find biomarkers that track anxiety, regardless if the reason for it is endogenous biology or it is driven by medications or drugs. In fact, one would expect some of these biomarkers to be targets of medications, as we show in this paper.
Moreover, the prioritization step that occurs after discovery is based on a field-wide convergence with literature that includes genetic data and animal model data, that are unrelated to medication effects. Overall, the discovery, validation and replication by testing in independent cohorts of the biomarkers, with our design, occurs despite the subjects having different genders, diagnoses, being on various different medications, and other lifestyle variables.
Blood gene expression experiments
[0044] RNA extraction. Whole blood (2.5 ml) was collected into each PaxGene tube by routine venipuncture. PaxGene tubes contain proprietary reagents for the stabilization of RNA. Total RNA was extracted and processed as previously described.
Microarrays. Microarray work was carried out using previously described methodology1.
[0045] Of note, all genomic data was normalized (RMA for technical variability, then z- scoring for biological variability), by gender and psychiatric diagnosis, before being combined and analyzed.
Biomarkers
Step 1: Discovery
[0046] We have used the subject’s score from a visual-analog scale (SAS-4) scale, assessed at the time of blood collection (Figure 1). We analyzed gene expression differences between visits with low anxiety (defined as a score of 0-40) and visits with high anxiety (defined as a score of 60 -100), using a powerful wi thin-subject design, then an across-subjects
8
SUBSTITUTE SHEET (RULE 26)
summation (Figure 1).
[0047] We analyzed the data in two ways: an Absent-Present (AP) approach, and a differential expression (DE) approach, as in previous work by us on suicide biomarkers. The AP approach may capture turning on and off of genes, and the DE approach may capture gradual changes in expression. Analyses were performed as previously described2'4. In brief, we imported all Affymetrix microarray data as CEL. files into Partek Genomic Suites 6.6 software package (Partek Incorporated, St Louis, MI, USA). Using only the perfect match values, we ran a robust multi-array analysis (RMA) by gender and diagnosis, background corrected with quantile normalization and a median polish probeset summarization of all chips, to obtain the normalized expression levels of all probesets for each chip. Then, to establish a list of differentially expressed probesets we conducted a within-subject analysis, using a fold change in expression of at least 1.2 between consecutive high- and low anxiety visits within each subject. Probesets that have a 1.2-fold change are then assigned either a 1 (increased in high anxiety) or a -1 (decreased in high anxiety) in each comparison. Fold changes between 1.1 and 1.2 are given 0.5, and fold changes less than 1.1 are given 0. These values were then summed for each probeset across all the comparisons and subjects, yielding a range of raw scores. The probesets above the 33.3% of raw scores were carried forward in analyses (Figure 1), and received an internal score of 2 points; those above 50% 4 points, and those above 80% 6 points2 3 4. We have developed in our labs R scripts to automate and conduct all these large dataset analyses in bulk, checked against human manual scoring4.
[0048] Gene Symbol for the probesets were identified using NetAffyx (Affymetrix) for Affymetrix HG-U133 Plus 2.0 GeneChips, followed by GeneCards to confirm the primary gene symbol. In addition, for those probesets that were not assigned a gene symbol by NetAffyx, we used GeneAnnot, or if need be UCSC to obtain gene symbol for these uncharacterized probesets, followed by GeneCard. Genes were then scored using our manually curated CFG databases as described below (Figure 1C).
Step 2: Prioritization using Convergent Functional Genomics (CFG) Databases.
[0049] We have established in our laboratory (Laboratory of Neurophenomics) manually curated databases of the human gene expression/protein expression studies (postmortem brain, peripheral tissue/fluids : CSF, blood and cell cultures), human genetic studies (association, copy
9
SUBSTITUTE SHEET (RULE 26)
number variations and linkage), and animal model gene expression and genetic studies, published to date on psychiatric disorders. Only findings deemed significant in the primary publication, by the study authors, using their particular experimental design and thresholds, are included in our databases. Our databases include only primary literature data and do not include review papers or other secondary data integration analyses to avoid redundancy and circularity. We also favored unbiased discovery studies over candidate genes hypothesis-driven studies. These large and constantly updated databases have been used in our CFG cross validation and prioritization platform (Figure IE). For this study, data from 354 papers on anxiety were present in the databases at the time of the CFG analyses (July 2019) (human genetic studies-93, human brain studies-10, human peripheral tissue/fluids- 96, non-human genetic studies-17, non-human brain studies-123, non-human peripheral tissue/fluids- 17). Analyses were performed as previously described. We have developed in our lab a computerized CFG Wizard to automate and score in bulk large lists of genes by integrating evidence from these large databases, checked against manual scoring. Analyses were performed as previously described2, 3.
Step 3: Validation analyses
[0050] We examined which of the top candidate genes (score of 6 or above after the first two steps), were changed in expression even more in an independent validation cohort (n=40) of clinically severe anxiety as measured by a STAI State >55, as well as a SAS-4 > 60. A total score of 6 or above after the first two steps permits the inclusion of potentially novel genes with maximal internal score of 6 from Discovery but no external evidence CFG score from Prioritization.
[0051] Subjects with low anxiety as well as subjects with high anxiety from the discovery cohort who did not have clinically severe anxiety were used, along with the independent validation cohort. We looked for stepwise change from the discovery cohort low anxiety group to the discovery cohort high anxiety group to the validation cohort clinically severe group.
[0052] The AP derived and DE derived lists of genes were combined, and the gene expression data corresponding to them was used for the validation analysis. The 3 groups (low anxiety, high anxiety, clinical severe anxiety) were assembled out of Affymetrix. cell data that was RMA normalized by gender and diagnosis. We transferred the log transformed expression data to an Excel sheet, and non-log transformed the data by taking 2 to the power of the
10
SUBSTITUTE SHEET (RULE 26)
transformed expression value. We then Z-scored the values by gender and diagnosis. We then imported the Excel sheets with the Z-scored by gender and diagnosis expression data into Partek, and statistical analyses were performed using a one-way ANOVA for the stepwise changed probesets, and also did a stringent Bonferroni correction for all the probesets tested in ANOVA (Figure 1).
Top Candidate Biomarkers (after the first 3 Steps)
[0053] Adding the scores from the first three steps into an overal convergent functional evidence (CFE) score (Figure 1), we ended up with a list of 95 top candidate biomarkers (95 probesets in n=82 genes), that had a CFE score of 8 and above (out of a maximum possible of 24 ; >33%) (see also Supplementary Information- Pathways, Predictions and Reproducibility).
These 95 top candidate biomarkers were carried forward into additional analyses for biological understanding, and for testing for clinical utility (Step 4).
Biological Understanding
Pathway Analyses
[0054] IPA (Ingenuity Pathway Analysis, version 24390178, Qiagen), David Functional Annotation Bioinformatics Microarray Analysis (National Institute of Allergy and Infectious Diseases) version 6.7 (August 2016), and Kyoto Encyclopedia of Genes and Genomes (KEGG) (through DAVID) were used to analyze the biological roles, including top canonical pathways and diseases (Table 3).
Networks
[0055] For network analyses we performed STRING Interaction network (https ://string- db.org) by inputting the genes into the search window and performed Multiple Proteins Homo sapiens analysis (Figure 3).
CFG beyond Anxiety : evidence for involvement in other psychiatric and related disorders.
[0056] We also used a CFG approach to examine evidence from other psychiatric and related disorders, as exemplified for the list of top biomarkers after Step 4 testing (Table 2). This was not used to prioritize genes, but rather to understand the molecular basis of clinical comorbidities.
Testing for Clinical Utility in Independent Cohorts
[0057] We tested in independent cohorts of psychiatric patients the ability of each of the top candidate biomarkers (n=95) to assess current severity of anxiety (state -measured by SAS-4
11
SUBSTITUTE SHEET (RULE 26)
or STAI), and predict future risk of anxiety (trait -future hospitalizations with anxiety). We conducted our analyses across all patients, as well as personalized by gender and diagnosis. [0058] The test cohort for predicting high anxiety (state), and the test cohort for predicting future hospitalizations with anxiety (trait), were assembled out of data that was RMA normalized by gender and diagnosis. The cohort was completely independent from the discovery and validation cohorts, there was no subject overlap with them. Individual markers used for predictions were Z scored by gender and diagnosis, to be able to combine different biomarkers into panels and to avoid potential artefacts due to different ranges of expression in different gender and diagnoses. For panels, biomarkers were combined by simple summation of the increased risk biomarkers minus the decreased risk biomarkers. Predictions were performed using R-studio. For cross-sectional analyses, we used biomarker expression levels, z-scored by gender and diagnosis. For longitudinal analyses, we combined four measures: biomarker expression levels, slope (defined as ratio of levels at current testing visit vs. previous visit, divided by time between visits), maximum levels (at any of the current or past visits), and maximum slope (between any adjacent current or past visits). For decreased biomarkers, we used the minimum rather than the maximum for level calculations. All four measures were Z-scored, then combined in an additive fashion into a single measure. The longitudinal analysis was carried out in a sub-cohort of the testing cohort consisting of subjects that had at least two visits (timepoints).
Predicting State- High Anxiety.
[0059] Receiver-operating characteristic (ROC) analyses between marker levels and anxiety state were performed by assigning subjects visits with a anxiety SAS-4 score of >60 into the high anxiety category, and subjects with STAI scores >55 in the high anxiety category. We used the pROC package of R (Xavier Robin et al. BMC Bioinformatics 2011). (Table 3, Figure 2). Additionally, a one-tailed t-test was performed between high anxiety group vs. the rest, and Pearson R (one-tail) was calculated between anxiety scores and biomarker levels.
[0060] Referring now to Fig. 2. Figure 2, Best Single Biomarkers Predictors for Anxiety State and Trait. From top candidate biomarkers after Steps 1-3 (Discovery, Prioritization, Validation-Bold) (n=95). Bar graph shows best predictive biomarkers in each group. All markers with * are nominally significant p<0.05. Table underneath the figures displays the actual number of biomarkers for each group whose ROC AUC p-values (2 A-C) and Cox Odds Ratio p-values
12
SUBSTITUTE SHEET (RULE 26)
(2 D) are at least nominally significant. Some gender and diagnosis group are missing from the graph as they did not have any significant biomarkers or that the cohort was too small with limited data for the z-scoring by gender-dx. Cross-sectional is based on levels at one visit. Longitudinal is based on levels at multiple visits (integrates levels at most recent visit, maximum levels, slope into most recent visit, and maximum slope). Dividing lines represent the cutoffs for a test performing at chance levels (white), and at the same level as the best biomarkers for all subjects in cross-sectional (gray) and longitudinal (black) based predictions. Biomarkers perform better than chance. Biomarkers performed better when personalized by gender and diagnosis. * survived Bonferroni correction for the number of candidate biomarkers tested. Predicting Trait- Future Psychiatric Hospitalization with Anxiety as a Symptom/Reason for Admission.
[0061] We conducted analyses for predicting future psychiatric hospitalizations with anxiety as a symptom/reason for the visit or admission in the first year following each testing visit, in subjects that had at least one year of follow-up in the VA system, in which we have access to complete electronic medical records. ROC analyses between biomarkers measures (cross- sectional, longitudinal) at a specific testing visit and future hospitalization admission were performed as described above, based on assigning if subjects had been admitted to the hospital with depression or not. Additionally, a one tailed t-test with unequal variance was performed between groups of subject visits with and without future hospitalization with anxiety. Pearson R (one-tail) correlation was performed between Hospitalization frequency (number of Hospitalization with anxiety divided by duration of follow-up) and marker levels. A Cox regression was performed using the time in days from the testing visit date to first hospitalization date in the case of patients who had been hospitalized, or 365 days for those who did not. The odds ratio was calculated such that a value greater than 1 always indicates increased risk for hospitalization, regardless if the biomarker is increased or decreased in expression.
[0062] We also conducted Cox regression and Pearson R analyses for all future hospitalizations with anxiety, including those occurring beyond one year of follow-up, in the years following testing (on average 7.35 years per subject, range 0.07 to 14.74 years), as these calculations, unlike the ROC and t-test, account for the actual length of follow-up, which varied from subject to subject. The ROC and t-test might in fact, if used, under-represent the power of the markers to predict, as the more severe psychiatric patients are more likely to move
13
SUBSTITUTE SHEET (RULE 26)
geographically and/or be lost to follow-up. The Cox regression was performed using the time in days from visit date to first hospitalization date in the case of patients who had hospitalizations with anxiety, or from visit date to last note date in the electronic medical records for those who did not.
Therapeutics
Pharmacogenomics.
[0063] We analyzed which of the top biomarkers for anxiety after Steps 1- 4 are known to be changed in expression by existing drugs in a direction opposite to the one in disease, using our CFG databases.
New drug discovery/repurposing
[0064] We also analyzed which drugs and natural compounds are an opposite match for the gene expression signatures of our top biomarkers, using the Connectivity Map (Broad Institute, MIT) (Table 4). Of note, not all the probesets from the HG-U133 Plus 2.0 array we used were present in the HGU-133A array used for the Connectivity Map. We stayed with exact probeset level matches, not gene level imputation. We also used the NIH LINCS tools database L1000CDS2 https://maayanlab.cloud/L1000CDS2/#/index to conduct similar analyses, at a gene level.
Report generation
[0065] We present an example of how a report to doctors might look, using the above insights. We chose as a case study a visit from a female subject ( phchp328vl) with anxiety and depression who had died by suicide, a case previously discussed in a suicide biomarker paper of ours (Levey et al. 20163) (Figure 4). We used the panel of the top biomarkers for anxiety from Table 2 (n= 19).
Example of report for physicians
[0066] Referring now to Figure 4: Example of report for physicians. Using the panel of the top biomarkers for anxiety from Table 2 (n=19) . This subject (Phchp328) was previously described by us in a suicidality biomarker study (Levey et al. 2016), as high risk for suicide, and died by suicide a year after completing our study. No information was provided to the patient’s clinicians by us at that time due to anonymity and privacy rules in research studies.
[0067] The raw expression values of the 19 biomarkers for 794 microarrays gene expression were Z-scored by gender and diagnosis. We calculated as thresholds the average expression
14
SUBSTITUTE SHEET (RULE 26)
value for a biomarker in the high anxiety group SAS4 >60, and in the low anxiety group SAS4 < 40. The first average should be higher than the second average in increased biomarkers, and the reverse is true for decreased biomarkers. 15 out of 19 biomarkers were thus concordant.
[0068] We also calculated as thresholds the average expression value for a biomarker in the first-year hospitalizations group, and in the not hospitalized in the firs- year group. We did the same thing for all future hospitalizations, and no future hospitalizations. The first average should be higher than the second average in increased biomarkers, and the reverse is true for decreased biomarkers. 18 out of 19 biomarkers were thus concordant for first year, and for all future.
[0069] The Z-scored expression value of each increased in expression biomarker was compared to the average value for the biomarker in the high anxiety group SAS4 >60, and the average value of the low anxiety group SAS4 < 40, resulting in scores of 1 if above high anxiety, 0 if below low anxiety, and 0.5 if it was in between. The reverse was done for decreased in expression biomarkers. This digitalization of the scores was done to avoid overfitting to our particular cohort and provide an easily understandable and interpretable readout for clinicians.
[0070] The digitized biomarkers were then added into a polygenic risk score and normalized for the number of biomarkers in the panel, resulting in a percentile score for anxiety. We did the same thing for first year hospitalizations, and all future hospitalizations, generating a combined score for chronic anxiety risk.
[0071] The digitized biomarkers were also used for matching with existing psychiatric medications and alternative treatments (nutraceuticals and others). We used our large datasets and literature databases to match biomarkers to medications that had effects on gene expression opposite to their expression in high anxiety. Each medication matched to a biomarker got the biomarker score of 1, 0.5 or 0. The scores for the medications were added, normalized for the number of biomarkers that were 1 or 0.5 in that patient, resulting in a percentile match. Thus, psychiatric medications matched to the patient and ranked in order of impact on the panel.
Results
[0072] In Step 1 Discovery, we identified candidate blood gene expression biomarkers that: 1. change in expression in blood between self-reported low anxiety and high anxiety states, 2. track the anxiety state across visits in a subject, and 3. track the anxiety state in multiple subjects. We used a visual analog measure for anxiety state (SAS-4). At a phenotypic level, the SAS-4 quantitates anxiety state at a particular moment in time, and normalizes anxiety
15
SUBSTITUTE SHEET (RULE 26)
measurements in each subject, comparing them to the lowest and highest anxiety that subject ever experienced (Figure SI). It has a moderate to strong correlation with current clinical scales for anxiety state (the STAI State, Figure S2).
[0073] We used a powerful within -subject and then across-subject design in a longitudinally followed cohort of subjects (n= 58 subjects, with 149 visits) who displayed at least a 50% change in the anxiety measure (from below 40/100 to above 60/100) between at least two consecutive testing visits, to identify differentially expressed genes that track anxiety state. Using our 33% of maximum raw score threshold (internal score of 2pt) 2’3, we identified 10,573 unique probesets (corresponding to 7195 unique genes) from Affymetrix Absent/Present (AP) analyses and Differential Expression (DE) analyses (Figure 1). These were carried forward to the prioritization step. This represents approximately a 5-fold enrichment of the 54,625 probesets on the Affymetrix array.
[0074] In Step 2 Prioritization, we used a Convergent Functional Genomics (CFG) approach to prioritize the candidate biomarkers identified in the discovery step (33% cutoff, internal score of >2pt.) by using published literature evidence (genetic, gene expression and proteomic), from human and animal model studies, for involvement in anxiety disorders (Figure 1 and Table 2). For Axniety there were 284 probesets (corresponding to 238 unique genes) that had a total score (combined discovery score and prioritization CFG score) of 6 and above. These were carried forward to the validation step. This represents approximately a 10-fold enrichment of the probesets on the Affymetrix array.
[0075] In Step 3 Validation, we validated for change in clinically severe anxiety these prioritized biomarkers, in a demographically matched cohort of (n=40 clinically severe anxiety) by assessing which markers were stepwise changed in expression: from clinically severe anxiety in validation cohort, to high anixety in discovery cohort, to low anxiety in discovery cohort in the validation cohort (Figure 1). 224 probesets were not stepwise changed, and 57 were stepwise changed. Of these, 4 probesets (corresponding to 4 unique genes) were nominally significant. [0076] Adding the scores from the first three steps into an overall convergent functional evidence (CFE) score (Figure 1), we ended up with a list of 95 top candidate biomarkers (n=82 genes, 95 probesets), that had a CFE3 score >= 8, better than than 33% of the maximum possible score of 24 after the first three steps, which we decided to use as an empirical cutoff. This represents approximately an over 500-fold enrichment of the probesets on the Affymetrix array.
16
SUBSTITUTE SHEET (RULE 26)
These 95 top candidate biomarkers were carried forward into analyses for understanding biological underpinnings. Last but not least, they were tested for predictive ability and clinical utility in additional independent cohorts.
Biological Understanding
Biological Pathways.
[0077] We conducted biological pathway analyses using the two lists of top candidate biomarkers for anxiety (n=82 genes, 95 probesets), which suggest that cAMP signaling is involved (Table 3). Depression was a top disease identified by the pathway analyses using DAVID, pointing out to the issue of co-morbidity, and Ingenuity identified neurological disorders as the top diseases.
Networks and Interactions.
[0078] We conducted a STRING analysis (Figure 3) of the top candidate biomarkers that revealed groups of interacting proteins. In particular, HTR2A is at the overlap of a network containing GABBR1, GAD1 and SLC6A4, and one centered on CDH1 that also contains PIK3R1 and IGFR1. A third network includes DLGAP1, PTPRD and DYNLL2. These networks may have biological significance and could be targeted therapeutically.
[0079] Testing for Clinical Utility In Step 4 Testing, we examined in completely independent cohorts from the ones used for discovery or validation whether the 26 top candidate biomarkers can assess high anxiety state (n=190 subjects), as well as predict of future psychiatric hospitalizations with anxiety (n=170 subjects) (Figure 2, and Table 2), using electronic medical records follow-up data of our study subjects (over a decade from initial visit at the time of the analyses) (Figure 1, Table 1). The gene expression data in the test cohorts was normalized (Z-scored) across genders and various psychiatric diagnoses, before those different demographic groups were combined. We used biomarker levels information cross-sectionally, as well as expanded longitudinal information about biomarker levels at multiple visits, as predictors. We tested the biomarkers in all subjects in the independent test cohort, as well as in a more personalized fashion by gender and psychiatric diagnosis.
Convergent Functional Evidence (CFE) (Table 2 )
[0080] For the top candidate biomarkers (n=95), we tabulated into a convergent functional evidence (CFE4) score all the evidence from discovery (up to 6 points), CFG prioritization (up to 12 points), validation (up to 6 points), and testing (state high anxiety, state clinical anxiety, trait
17
SUBSTITUTE SHEET (RULE 26)
first year hospitalization with anxiety, trait all future hospitalizations with anxiety- up to 3 points each if it significantly predicts in all subjects, 2 points if in gender, 1 points if in gender/diagnosis). The total score can be up to 36 points: 24 from our empirical data, and 12 from literature data used for CFG. We weigh our new empirical data more than the literature data, as it is functionally related to mood in 3 independent cohorts (discovery, validation, testing). The goal is to highlight, based on the totality of our data and of the evidence in the field to date, biomarkers that have all around evidence: track anxiety, have convergent evidence for involvement in anxiety disorders, and predict anxiety state and future clinical events (Table 2). Such biomarkers merit priority evaluation in future clinical trials. The top biomarkers (n=l 9) are depicted in Table 2.
Targeted Therapeutics
Pharmacogenomics.
[0081] A number of individual top biomarkers are known to be modulated by medications in current clinical use and by nutraceuticals (Table 6, Figure 4). This is of potential utility in patient stratification and pharmacogenomics approaches. Omega- 3 fatty acids may be a widely deployable preventive treatment, with minimal side-effects, including in women who are or may become pregnant.
New drug discovery/repurposing.
[0082] Bioinformatic analyses using the gene expression signature of the panel of top biomarkers for anxiety (Table 4) identified new potential therapeutics for anxiety, such as estrogen.
18
SUBSTITUTE SHEET (RULE 26)
Table 1. Aggregate Demographics, BP-bipolar, MDD- depression, MOOD- mood nos, SZ- schizophrenia, SZA- schizoaffective, PSYCH- psychosis nos, PTSD-post-traumatic stress disorder.
SUBSTITUTE SHEET (RULE 26)
Table 2. Top Anxiety Biomarkers: Convergent Functional Evidence (CFE). After Step 4 Testing in independent cohorts for state and trait predictive ability. For Step 4 Predictions, C: -cross- sectional (using levels from one visit), L: -longitudinal (using levels and slopes from multiple visits). In ALL, by Gender, and personalized by Gender and Diagnosis (Gender/Dx).M-Males, F- Females. MDD-depression, BP-bipolar, SZ-schizophrenia, SZA-schizoaffective, PSYCHOSIS- schizophrenia and schizoaffective combined, PTSD-post-traumatic stress disorder.
20
24
25
26
SUBSTITUTE SHEET (RULE 26)
Table 3 A, and B. Biology of Anxiety Biomarkers. Top CFE3 =>8 (n= 95 probesets, 82 genes).
SUBSTITUTE SHEET (RULE 26)
Table 4: Therapeutics: Drug repurposing for Anxiety using Connectivity Map5 (CMAP) For Top Biomarkers from Table 2 (n=19). Direction of expression in high anxiety. 2 out of 4 Decreased and 10 out of 15 Increased probesets were present in HG-U133A, the array used for
CMAP (from genes ADRA2A, ATP1B2, CCKBR, FZD1O,GAD1,GPX7,GRK4 ,NRG1,NTRK3,SLC6A2,TFRC). Drugs that have opposite gene expression effects to the gene expression signature.
Table 5. Some predictive biomarkers for anxiety determined from cross-sectional and longitudinal studies.
36
SUBSTITUTE SHEET (RULE 26)
Table 6. Pharmacogenomics. Top biomarkers (from Table 2) that are targets of existing drugs and are changed in expression in opposite direction to high anxiety. (I)- increased in expression, (D)- decreased in expression. DE- differential expression, AP-Absent/Present.
39
40
Discussion
[0083] We describe a novel and comprehensive effort to discover and validate blood biomarkers of relevance to anxiety, including testing them in independent cohorts to evaluate predictive ability and clinical utility. These biomarkers also open a window into understanding the biology of anxiety disorders, as well as indicate new and more precise therapeutic approaches. Biomarkers may be useful for matching patients to medications and measuring response to treatment (pharmacogenomics) (Figure 4, Tables 2 and 6), as well as new drug discovery and repurposing (Table 4).
[0084] Overall, this work is a major step forward towards understanding, diagnosing and treating anxiety disorders. We hope that our trait biomarkers for future risk may be useful in preventive approaches before the full-blown disorder manifests itself (or re-occurs). Prevention could be accomplished with social, psychological, or biological interventions (i.e. early targeted use of medications or nutraceuticals). Given the fact that 1 in 5 people will have a clinical
42
SUBSTITUTE SHEET (RULE 26)
anxiety disorders episode in their lifetime, that anxiety disorders can severely affect quality of life, sometimes leading to suicides, and that not all patients respond to current treatments or the treatments are addictive, the need for and importance of efforts such as ours cannot be overstated.
43
SUBSTITUTE SHEET (RULE 26)
Claims
1. A therapeutic compound, selected from one or more compounds from Tables 4 and Table 6, in a method of mitigating and anxiety disorder in a subject in need thereof for example wherein the therapy is a drug, a natural compound, and combinations thereof.
2. The therapeutic compound, for use according to claim 1, wherein the drug is estradiol, loperamide, omage-3 fatty acids, lithium, fluoxetine, norfluoxetine, imipramine, citalopram, carbamazepine, aripiprazole, haloperidol, mianserin, SAM, vortioxetine, agomelatine, sertraline, benzodiazepine, clozapine, ketamine, valproate, clozapine, rofecoxib, berberine, methylphenidate, or berberine.
3. The therapeutic compound, for use according to claim 1, wherein the subject has anxiety disorder.
4. The therapeutic compound, for use according to claim 1, wherein the subject is a male subject or a female subject.
5. A computer-implemented method for assessing the level of anxiety in a subject, and for assessing risk of developing anxiety disorder in the future in a subject, the method comprising: computing a score based on RNA level, protein level, DNA methylation, a single nucleotide polymorphism, for a panel of biomarkers chosen from Table 2, Table 5, and Figure 2, in a sample obtained from a subject; computing a score based on a reference expression level of the panel of biomarkers; and identifying a difference between the score in the sample obtained from the subject and the score in the reference sample, wherein the difference in the score in the sample obtained from the subject and the score in the reference sample indicates an anxiety disorder, and a risk of future anxiety disorders.
44
SUBSTITUTE SHEET (RULE 26)
6. The computer implemented method according to claim 5, wherein the subject is: a male or a female and the biomarker is at least one biomarker selected from the group consisting of: ERCC6L2, SLC6A2, FZD10, and RPL14; or a female and the biomarker is at least one biomarker selected from the group consisting of: PHF21A, SLC6A2, 240253_at, and TRIM9; or a male and the biomarker is at least one biomarker selected from the group consisting of: ERCC6L2, ANKRD28, ACPI, and RP
7. The methods of claims 5 and 6, further including the steps of: identifying a difference between the score in the sample obtained from the subject and the score in the reference sample, administering treatment to the subject, wherein the treatment reduces the difference between the score in the sample from the subject and the score in the reference sample, and mitigates anxiety in the subject, or the risk for developing anxiety disorder in the subject.
8. The method of claim 7, further including, the step of: measuring a change in the score after treatment.
9. The method of claim 8, wherein the treatment is at least one treatment selected from the group consisting of: lifestyle modification, electro-magnetic intervention, and administering a therapeutic compound.
10. The method of claim 9, wherein the therapy is selected by a computer- implemented method selected from the group consisting of one or more psychiatric compounds from Table 6, and wherein each therapy selection is based on a panel of one or more individual biomarkers.
11. The method of claim 7, wherein the therapy is selected based on a panel of individual biomarkers changed in a subject, by a computer-implemented method for therapy selection, and comprises, of one or more compounds in Table 4.
12. A method for assessing and mitigating anxiety disorder, comprising determining an expression level of a panel of biomarkers listed in Table 2, Table 5 and Figure 2, in a sample for example wherein the sample comprises a peripheral tissue, blood, saliva, cerebrospinal fluid
45
SUBSTITUTE SHEET (RULE 26)
(CSF), serum, urine, or stool; wherein the expression level of the biomarkers in the sample is different relative to a reference expression level; and includes the steps of: identifying the subject currently having or at risk of having an anxiety disorder based on a biomarker panel score relative to a biomarker panel score of a reference panel; selecting a therapy based on the score from the biomarkers; and administering one or more compounds selected from the group consisting of compounds selected from the groups consisting of compounds selected from Tables 4 and Table 6.
13. A composition comprising one or more compounds from Tables 4 and Table 6, for use in a method for treating anxiety disorder.
46
SUBSTITUTE SHEET (RULE 26)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/US2022/038673 WO2024025536A1 (en) | 2022-07-28 | 2022-07-28 | Precision medicine for anxiety disorders: objective assessment, risk prediction, pharmacogenomics, and repurposed drugs |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/US2022/038673 WO2024025536A1 (en) | 2022-07-28 | 2022-07-28 | Precision medicine for anxiety disorders: objective assessment, risk prediction, pharmacogenomics, and repurposed drugs |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2024025536A1 true WO2024025536A1 (en) | 2024-02-01 |
Family
ID=89707011
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/US2022/038673 WO2024025536A1 (en) | 2022-07-28 | 2022-07-28 | Precision medicine for anxiety disorders: objective assessment, risk prediction, pharmacogenomics, and repurposed drugs |
Country Status (1)
Country | Link |
---|---|
WO (1) | WO2024025536A1 (en) |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20140274764A1 (en) * | 2013-03-15 | 2014-09-18 | Pathway Genomics Corporation | Method and system to predict response to treatments for mental disorders |
WO2018209341A1 (en) * | 2017-05-12 | 2018-11-15 | Indiana University Research And Technology Corporation | Precision medicine for treating and preventing suicidality |
-
2022
- 2022-07-28 WO PCT/US2022/038673 patent/WO2024025536A1/en unknown
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20140274764A1 (en) * | 2013-03-15 | 2014-09-18 | Pathway Genomics Corporation | Method and system to predict response to treatments for mental disorders |
WO2018209341A1 (en) * | 2017-05-12 | 2018-11-15 | Indiana University Research And Technology Corporation | Precision medicine for treating and preventing suicidality |
Non-Patent Citations (2)
Title |
---|
HOMMERS L. G., RICHTER J., YANG Y., RAAB A., BAUMANN C., LANG K., SCHIELE M. A., WEBER H., WITTMANN A., WOLF C., ALPERS G. W., ARO: "A functional genetic variation of SLC6A2 repressor hsa-miR-579-3p upregulates sympathetic noradrenergic processes of fear and anxiety", TRANSLATIONAL PSYCHIATRY, GB, vol. 8, no. 1, GB , XP093136733, ISSN: 2158-3188, DOI: 10.1038/s41398-018-0278-4 * |
ZHANG SHENGLI, CHEN ZHIHAN, WU YANFANG, LIN DIANTIAN, HE JUANJUAN, LIU JIANWEN, GAO FEI, LIN HE: "Efficacy and safety of thalidomide on psychological symptoms and sleep disturbances in the patient with refractory ankylosing spondylitis", ANNALS OF PALLIATIVE MEDICINE, vol. 10, no. 3, 1 March 2021 (2021-03-01), pages 2512 - 2519, XP093136730, ISSN: 2224-5820, DOI: 10.21037/apm-19-409 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Le-Niculescu et al. | Precision medicine for mood disorders: objective assessment, risk prediction, pharmacogenomics, and repurposed drugs | |
Yang et al. | Genomic atlas of the proteome from brain, CSF and plasma prioritizes proteins implicated in neurological disorders | |
Smith et al. | Classification of genes: standardized clinical validity assessment of gene–disease associations aids diagnostic exome analysis and reclassifications | |
Saw et al. | Chromosome 1q21. 2 and additional loci influence risk of spontaneous coronary artery dissection and myocardial infarction | |
Zhao et al. | Short leukocyte telomere length predicts risk of diabetes in american indians: the strong heart family study | |
Endam et al. | Association of IL1A, IL1B, and TNF gene polymorphisms with chronic rhinosinusitis with and without nasal polyposis: a replication study | |
Xie et al. | Interactive effect of stressful life events and the serotonin transporter 5-HTTLPR genotype on posttraumatic stress disorder diagnosis in 2 independent populations | |
Dupuis et al. | Prevalence of meconium ileus marks the severity of mutations of the Cystic Fibrosis Transmembrane Conductance Regulator (CFTR) gene | |
Karason et al. | A strong heritability of psoriatic arthritis over four generations—the Reykjavik Psoriatic Arthritis Study | |
US10719579B2 (en) | Validating biomarker measurement | |
Allen-Brady et al. | Significant linkage evidence for a predisposition gene for pelvic floor disorders on chromosome 9q21 | |
WO2018160548A1 (en) | Markers for coronary artery disease and uses thereof | |
Liebers et al. | Discriminating bipolar depression from major depressive disorder with polygenic risk scores | |
C. Silva et al. | Cross-tissue analysis of blood and brain epigenome-wide association studies in Alzheimer’s disease | |
Yang et al. | Genomic and multi-tissue proteomic integration for understanding the biology of disease and other complex traits | |
McGrath et al. | Insights from Mendelian randomization and genetic correlation analyses into the relationship between endometriosis and its comorbidities | |
Akbarzadeh et al. | Evaluating machine learning-powered classification algorithms which utilize variants in the GCKR gene to predict metabolic syndrome: Tehran Cardio-metabolic Genetics Study | |
Yan et al. | ICH-LR2S2: a new risk score for predicting stroke-associated pneumonia from spontaneous intracerebral hemorrhage | |
Chen et al. | Identification and validation of BCL6 and VEGFA as biomarkers and ageing patterns correlating with immune infiltrates in OA progression | |
Edmiston et al. | Gene expression profiling of peripheral blood leukocytes identifies potential novel biomarkers of chronic obstructive pulmonary disease in current and former smokers | |
US20230265517A1 (en) | Novel dna methylation markers associated with renal function and method for predictiing renal function | |
WO2020237203A1 (en) | Methods for objective assessment of memory, early detection of risk for alzheimer's disease, matching individuals with treatments, monitoring response to treatment, and new methods of use for drugs | |
WO2019240997A1 (en) | Methods for assessment and early detection of stress, selecting and monitoring treatment, and new use for drugs | |
WO2024025536A1 (en) | Precision medicine for anxiety disorders: objective assessment, risk prediction, pharmacogenomics, and repurposed drugs | |
US20230220472A1 (en) | Deterimining risk of spontaneous coronary artery dissection and myocardial infarction and sysems and methods of use thereof |
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: 22953318 Country of ref document: EP Kind code of ref document: A1 |