The invention relates to methods for diagnosing and treating psychiatric disorders in a patient. In particular, the intention relates to methods for identifying genetic markers that are associated with treatment resistant psychiatric disorders and administering a treatment for the treatment resistant psychiatric disorder to the patient based on identifying the genetic markers. The genetic markers may include polymorphisms such as single nucleotide polymorphisms (SNPs) and the alternate treatment may include foregoing administering typical or atypical antipsychotic drugs (APDs), where the patient has been diagnosed with a treatment resistant psychiatric disorder, and instead administering a different treatment regimen or a new APD as disclosed herein.
Psychiatric disorders such as schizophrenia, bipolar disorder, and depression with psychotic features are estimated to affect millions or Americans. Many patients with psychiatric disorders are successfully treated with typical or atypical APDs. For example, following treatment with typical APDs (i.e., first generation APDs) or atypical APDs (i.e., second generation APDs), about 70% of people with schizophrenia have significant remission of positive symptoms.69 However, many patients with psychiatric disorders are resistant to treatment with typical or atypical APDs. For example, about 30% of people with schizophrenia are resistant to treatment with APDs, where treatment resistant schizophrenia (TRS) may be defined as presence of moderate to severe positive symptoms despite treatment with two or more trials of APDs other than clozapine at usually adequate doses for 4-6 weeks.6 Some population studies have suggested that an even higher proportion of treatment resistant (TR) patients exists.1-4
Furthermore, treatment with certain APDs may have undesirable side effects. For example, olanzapine, clozapine, and quetiapine may cause an increased risk of developing metabolic abnormalities such as the impaired glucose metabolism/diabetes and weight gain/obesity.48 As such, if a patient has a treatment resistant psychiatric disorder, administering some APDs is not beneficial and may be detrimental. However, there currently are no available diagnostic tests for predicting whether a patient having a psychiatric disorder will be resistant to treatment with APDs. Here, genetic markers that are associated with treatment resistant psychiatric disorders are identified. By using these genetic markers, clinicians can stratify psychiatric patients based on the markers and provide a more personalized treatment for patients with treatment resistant psychiatric disorders and avoid administering typical or atypical APDs and instead administer a different treatment regimen or a new APD as disclosed herein.
Disclosed are methods, compositions, and kits for diagnosing, prognosing, and treating psychiatric disorders. The disclosed methods may include identifying genetic markers that may be utilized to diagnose or prognose a patient. Based on identifying the genetic markers, the disclosed methods further may include treating the diagnosed or prognosed patient by not administering a typical or atypical APD and instead administering a different treatment regimen, a APD approved for treating a treatment resistant psychiatric disorder, such as clozapine, or a new APD such as a trace amine-associated receptor 1 (TAAR1) agonist as disclosed herein. The disclosed compositions and kits may include reagents or drugs for diagnosing or prognosing a patient or reagents or drugs for treating a patient.
The disclosed methods and kits may be utilized for diagnosing and treating patients having a treatment resistant psychiatric disorder. The methods may include assessing whether a patient has a treatment resistant psychiatric disorder or assessing the likelihood of the patient developing a treatment resistant psychiatric disorder. The methods may include detecting genetic markers in a genomic nucleic acid sample from the patient and/or receiving (e.g., as a caregiver) results of a test indicating whether a genetic marker has been detected in a genomic nucleic acid sample from the patient. The methods may include detecting, either directly or indirectly, the presence or absence of one or more polymorphisms (e.g., specific single nucleotide polymorphisms (SNPs)) in the genomic nucleic acid sample and/or receiving (e.g., as a caregiver) results of a test indicating whether a polymorphism (e.g., a SNP) is present or absent in the genomic nucleic acid sample. The kits may include reagents for practicing the claimed methods, such as reagents for detecting the specific polymorphisms (e.g., a SNPs), and the reagents of the kits may be configured as an array of reagents.
The polymorphisms detected in the methods or by the kits may include one or more of the following SNPs: rs2237457, rs10825879, rs1421321, rs6763175, rs6500291, rs6924295, rs9675544, rs7085850, rs7933576, rs10047718, rs6043303, rs9313540, rs6432686, rs976760, rs7023648, rs388033, rs4724374, rs12510912, rs6480553, rs12369177, rs7775338, rs6421889, rs3926649, rs12584486, rs17776723, rs1563817, rs9817986, rs9818122, rs13100616, rs7135617, rs11609649, and rs9854706. Specifically, the SNPs detected in the methods or by the kits may include one or more of the following SNPs: rs2237457, rs10825879, rs1421321, rs6763175, rs6500291, rs6924295, and more specifically rs2237457.
In the methods and kits, the minor allele and/or the major allele associated with a SNP may be detected. The methods may include and the kits may be used for determining whether a patient is homozygous or heterozygous for a minor allele and/or major allele associated with a SNP.
In the methods and kits, the SNPs may be detected by any suitable method including biological methods and assays. These may include but are not limited to, nucleotide sequencing, probe hybridization, and primer specific PCR. A reagent for detecting a SNP may include an oligonucleotide that hybridizes specifically to the SNP to form a complex. The oligonucleotide probe further may include a modification or label to detect the hybridization complex. Additionally, a reagent for detecting a SNP may include an oligonucleotide that hybridizes within 100, 50, 20, 10, or 5 nucleotides 5′ or 3′ of the SNP, for example, where the oligonucleotide is utilized as a primer and detecting comprises sequencing the SNP. The oligonucleotide primer further may include a modification or label to detect the hybridization complex.
In the disclosed methods and kits, the patient may be resistant to treatment with typical APDs (i.e., “first generation” APDs) or atypical APDs (i.e., “second generation” APDs), which may include but are not limited to Benperidol (Anguil®, Benguil®, Frenactil®, Glianimon®), Bromperidol (Bromodol®, Impromen®), Droperidol (Droleptan®, Inapsine®), Haloperidol (Haldol®, Serenace®), Moperone (Luvatren®), Pipamperone (Dipiperon®, Piperonil®), Timiperone (Celmanil®, Tolopelon®), Diphenylbutylpiperidine, Fluspirilene (Imap®), Penfluridol (Semap®), Pimozide (Orap®), Acepromazine (Plegicil®), Chlorpromazine (Largactil®, Thorazine®), Cyamemazine (Tercian®), Dixyrazine (Esucos®), Fluphenazine (Modecate®, Permitil®, Prolixin®), Levomepromazine (Levinan®, Levoprome®, Nozinan®), Mesoridazine (Lidanil®, Serentil®), Perazine (Peragal®, Perazin®, Pernazinum®, Taxilan®), Pericyazine (Neulactil®, Neuleptil®), Perphenazine (Trilafon®), Pipotiazine (Lonseren®, Piportil®), Prochlorperazine (Compazine®), Promazine (Prozine®, Sparine®), Promethazine (Avomine®, Phenergan®), Prothipendyl (Dominal®), Thioproperazine (Majeptil®), Thioridazine (Aldazine®, Mellaril®, Melleril®), Trifluoperazine (Stelazine®), Triflupromazine (Vesprin®), Chlorprothixene (Cloxan®, Taractan®, Truxal®), Clopenthixol (Sordinol®), Flupentixol (Depixol®, Fluanxol®), Tiotixene (Navane®, Thixit®), Zuclopenthixol (Acuphase®, Cisordinol®, Clopixol®), Clotiapine (Entumine®, Etomine®, Etumine®), Loxapine (Adasuve®, Loxitane®), Prothipendyl (Dominal®), Carpipramine (Defekton®, Prazinil®), Clocapramine (Clofekton®, Padrasen®), Molindone (Moban®), Mosapramine (Cremin®), Sulpiride (Meresa®), Sultopride (Barnetil®, Topral®), Veralipride (Agreal®), Amisulpride (Solian®), Amoxapine (Asendin®), Aripiprazole (Abilify®), Asenapine (Saphris®, Sycrest®), Clozapine (Clozaril®), Blonanserin (Lonasen®), Iloperidone (Fanapt®, Fanapta®, Zomaril®), Lurasidone (Latuda®), Melperone (Buronil®, Buronon®, Eunerpan®, Melpax®, Neuril®), Olanzapine (Zyprexa®), Paliperidone (Invega®), Perospirone (Lullan®), Quetiapine (Seroquel®), Remoxipride (Roxiam®), Risperidone (Risperdal®), Sertindole (Serdolect®, Serlect®), Trimipramine (Surmontil®), Ziprasidone (Geodon®, Zeldox®), and Zotepine (Lodopin®, Losizopilon®, Nipolept®, Setous®).
In the disclosed methods, the patient may be administered a treatment for a treatment resistant psychiatric disorder, which may include foregoing administering the typical or atypical ASD based on detecting the SNP and/or receiving results of tests indicating that the SNP is present. Treatments and therapies other than administering typical or atypical ASDs may include but are not limited to administering psychosocial treatment and electroconvulsive therapy, administering a APD approved for treating treatment resistant psychiatric disorders (e.g., clozapine), and/or administering new ASDs as disclosed herein (e.g., trace amine-associated receptor 1 (TAAR1) agonists). Alternatively in the disclosed methods, where the SNP is not detected and/or results are receiving indicating that the SNP is not present, the patient may be administered the typical or atypical ASD.
The results disclosed herein implicate L-Dopa Decarboxylase in psychiatric disorders, in particular, treatment resistant psychiatric disorders such as treatment resistant schizophrenia. As such, the presently disclosed methods also include methods of treating a patient having or at risk for developing a psychiatric disorder (e.g., a treatment resistant psychiatric disorders such as such as treatment resistant schizophrenia) by administering to the patient an agent that modulates L-Dopa Decarboxylase expression or activity (e.g., an agent that increases L-Dopa Decarboxylase expression or activity). Because DCC is an enzyme that is rate-limiting in the synthesis of trace amines and neurotransmitters implicated in psychiatric disorders, the present methods also contemplate treating a patient having or at risk for developing a psychiatric disorder (e.g., a treatment resistant psychiatric disorders such as such as treatment resistant schizophrenia) by administering to the patient an agent that is an agonist of the trace amine-associated receptor 1 (TAAR1) for which trace amines are a ligand.
The present invention is described herein using several definitions, as set forth below and throughout the application.
As used in this specification and the claims, the singular forms “a,” “an,” and “the” include plural forms unless the content clearly dictates otherwise. For example, “an SNP,” should be interpreted to mean “one or more SNPs.”
As used herein, “about”, “approximately,” “substantially,” and “significantly” will be understood by persons of ordinary skill in the art and will vary to some extent on the context in which they are used. If there are uses of the term which are not clear to persons of ordinary skill in the art given the context in which it is used, “about” and “approximately” will mean up to plus or minus 10% of the particular term and “substantially” and “significantly” will mean more than plus or minus 10% of the particular term.
As used herein, the terms “include” and “including” have the same meaning as the terms “comprise” and “comprising.” The terms “comprise” and “comprising” should be interpreted as being “open” transitional terms that permit the inclusion of additional components further to those components recited in the claims. The terms “consist” and “consisting of” should be interpreted as being “closed” transitional terms that do not permit the inclusion additional components other than the components recited in the claims. The term “consisting essentially of” should be interpreted to be partially closed and allowing the inclusion only of additional components that do not fundamentally alter the nature of the claimed subject matter.
As used herein, the term “plurality” means “two or more.” For example, a “plurality of SNPs” should be interpreted to mean “2, 3, 4, 5 or more SNPs.”
As used herein, the term “patient,” which may be used interchangeably with the terms “subject” or “individual,” refers to one who receives medical care, attention or treatment and may encompass a human patient. As used herein, the term “patient” is meant to encompass a person who has a psychiatric disorder or is at risk for developing a psychiatric disorder, which includes but is not limited to schizophrenia, bipolar disorder, and depression (e.g., depression with psychotic features). For example, the term “patient” is meant to encompass a person at risk for developing schizophrenia or a person diagnosed with schizophrenia (e.g., a person who may be symptomatic for schizophrenia but who has not yet been diagnosed). The term “patient” also is meant to encompass a person at risk for developing bipolar disorder or a person diagnosed with bipolar disorder (e.g., a person who may be symptomatic for bipolar disorder but who has not yet been diagnosed). The term “patient” further is meant to encompass a person at risk for developing depression (e.g., depression with psychotic features). As such, the term “patient” further is meant to encompass a person at risk for developing depression with psychotic features or a person diagnosed with depression with psychotic features (e.g., a person who may be symptomatic for depression with psychotic features but who has not yet been diagnosed).
The term “patient” further is meant to encompass a person at risk for developing a treatment resistant psychiatric disorder, which may include but is not limited to schizophrenia, bipolar disorder, and depression with psychotic features. For example, the term “patient” is meant to encompass a person at risk for developing treatment resistant schizophrenia (TRS) or a person diagnosed with TRS (e.g., a person who may be symptomatic for TRS but who has not yet been diagnosed). The term “patient” also is meant to encompass a person at risk for developing treatment resistant bipolar disorder or a person diagnosed with treatment resistant bipolar disorder (e.g., a person who may be symptomatic for treatment resistant bipolar disorder but who has not yet been diagnosed). The term “patient” further is meant to encompass a person at risk for developing treatment resistant depression (e.g., treatment resistant depression with psychotic features). As such, the term “patient” further is meant to encompass a person at risk for developing treatment resistant depression with psychotic features or a person diagnosed with treatment resistant depression with psychotic features (e.g., a person who may be symptomatic for treatment resistant depression with psychotic features but who has not yet been diagnosed).
As used herein, the term “patient” also is meant to encompass a person who has a psychiatric disorder or is at risk for developing a psychiatric disorder that is associated with increased or decreased biological activity of L-Dopa Decarboxylase (DCC), which is an enzyme that is rate limiting in the synthesis of trace amines and neurotransmitters. For example, a patient having a psychiatric disorder that is associated with decreased biological activity of DCC may be treated by administering an agent that increases biological activity of DCC, and/or an agent that increases expression of DCC. Because DCC is an enzyme that is rate-limiting in the synthesis of trace amines and neurotransmitters implicated in psychiatric disorders, a patient having a psychiatric disorder that is associated with decreased biological activity of DCC may be treated by administering an agent that is an agonist of the trace amine-associated receptor 1 (TAAR1) for which trace amines are a ligand.
A patient may include a patient having or at risk for developing a psychiatric disease or disorder that is resistant to treatment with typical or atypical antipsychotic drugs (APDs). Typical and atypical APDs for treating psychiatric diseases or disorders may include, but are not limited to, drugs such as Benperidol (Anguil®, Benguil®, Frenactil®, Glianimon®), Bromperidol (Bromodol®, Impromen®), Droperidol (Droleptan®, Inapsine®), Haloperidol (Haldol®, Serenace®), Moperone (Luvatren®), Pipamperone (Dipiperon®, Piperonil®), Timiperone (Celmanil®, Tolopelon®), Diphenylbutylpiperidine, Fluspirilene (Imap®), Penfluridol (Semap®), Pimozide (Orap®), Acepromazine (Plegicil®), Chlorpromazine (Largactil®, Thorazine®), Cyamemazine (Tercian®), Dixyrazine (Esucos®), Fluphenazine (Modecate®, Permitil®, Prolixin®), Levomepromazine (Levinan®, Levoprome®, Nozinan®), Mesoridazine (Lidanil®, Serentil®), Perazine (Peragal®, Perazin®, Pernazinum®, Taxilan®), Pericyazine (Neulactil®, Neuleptil®), Perphenazine (Trilafon®), Pipotiazine (Lonseren®, Piportil®), Prochlorperazine (Compazine®), Promazine (Prozine®, Sparine®), Promethazine (Avomine®, Phenergan®), Prothipendyl (Dominal®), Thioproperazine (Majeptil®), Thioridazine (Aldazine®, Mellaril®, Melleril®), Trifluoperazine (Stelazine®), Triflupromazine (Vesprin®), Chlorprothixene (Cloxan®, Taractan®, Truxal®), Clopenthixol (Sordinol®), Flupentixol (Depixol®, Fluanxol®), Tiotixene (Navane®, Thixit®), Zuclopenthixol (Acuphase®, Cisordinol®, Clopixol®), Clotiapine (Entumine®, Etomine®, Etumine®), Loxapine (Adasuve®, Loxitane®), Prothipendyl (Dominal®), Carpipramine (Defekton®, Prazinil®), Clocapramine (Clofekton®, Padrasen®), Molindone (Moban®), Mosapramine (Cremin®), Sulpiride (Meresa®), Sultopride (Barnetil®, Topral®), Veralipride (Agreal®), Amisulpride (Solian®), Amoxapine (Asendin®), Aripiprazole (Abilify®), Asenapine (Saphris®, Sycrest®), Clozapine (Clozaril®), Blonanserin (Lonasen®), Iloperidone (Fanapt®, Fanapta®, Zomaril®), Lurasidone (Latuda®), Melperone (Buronil®, Buronon®, Eunerpan®, Melpax®, Neuril®), Olanzapine (Zyprexa®), Paliperidone (Invega®), Perospirone (Lullan®), Quetiapine (Seroquel®), Remoxipride (Roxiam®), Risperidone (Risperdal®), Sertindole (Serdolect®, Serlect®), Trimipramine (Surmontil®), Ziprasidone (Geodon®, Zeldox®), and Zotepine (Lodopin®, Losizopilon®, Nipolept®, Setous®). A patient having a treatment resistant psychiatric disease or disorder accordingly may include a patient that is resistant to treatment with, for example, any of these foregoing APDs.
In some embodiments of the disclosed methods, the patient may have undergone treatment prior to the disclosed methods being performed and the patient may have been diagnosed with a treatment resistant psychiatric disorder prior to the method being performed. For example, prior to disclosed methods being performed the patient may have been diagnosed and treated for schizophrenia using an antipsychotic drug (APD) selected from the group consisting of Benperidol (Anguil®, Benguil®, Frenactil®, Glianimon®), Bromperidol (Bromodol®, Impromen®), Droperidol (Droleptan®, Inapsine®), Haloperidol (Haldol®, Serenace®), Moperone (Luvatren®), Pipamperone (Dipiperon®, Piperonil®), Timiperone (Celmanil®, Tolopelon®), Diphenylbutylpiperidine, Fluspirilene (Imap®), Penfluridol (Semap®), Pimozide (Orap®), Acepromazine (Plegicil®), Chlorpromazine (Largactil®, Thorazine®), Cyamemazine (Tercian®), Dixyrazine (Esucos®), Fluphenazine (Modecate®, Permitil®, Prolixin®), Levomepromazine (Levinan®, Levoprome®, Nozinan®), Mesoridazine (Lidanil®, Serentil®), Perazine (Peragal®, Perazin®, Pernazinum®, Taxilan®), Pericyazine (Neulactil®, Neuleptil®), Perphenazine (Trilafon®), Pipotiazine (Lonseren®, Piportil®), Prochlorperazine (Compazine®), Promazine (Prozine®, Sparine®), Promethazine (Avomine®, Phenergan®), Prothipendyl (Dominal®), Thioproperazine (Majeptil®), Thioridazine (Aldazine®, Mellaril®, Melleril®), Trifluoperazine (Stelazine®), Triflupromazine (Vesprin®), Chlorprothixene (Cloxan®, Taractan®, Truxal®), Clopenthixol (Sordinol®), Flupentixol (Depixol®, Fluanxol®), Tiotixene (Navane®, Thixit®), Zuclopenthixol (Acuphase®, Cisordinol®, Clopixol®), Clotiapine (Entumine®, Etomine®, Etumine®), Loxapine (Adasuve®, Loxitane®), Prothipendyl (Dominal®), Carpipramine (Defekton®, Prazinil®), Clocapramine (Clofekton®, Padrasen®), Molindone (Moban®), Mosapramine (Cremin®), Sulpiride (Meresa®), Sultopride (Barnetil®, Topral®), Veralipride (Agreal®), Amisulpride (Solian®), Amoxapine (Asendin®), Aripiprazole (Abilify®), Asenapine (Saphris®, Sycrest®), Clozapine (Clozaril®), Blonanserin (Lonasen®), Iloperidone (Fanapt®, Fanapta®, Zomaril®), Lurasidone (Latuda®), Melperone (Buronil®, Buronon®, Eunerpan®, Melpax®, Neuril®), Olanzapine (Zyprexa®), Paliperidone (Invega®), Perospirone (Lullan®), Quetiapine (Seroquel®), Remoxipride (Roxiam®), Risperidone (Risperdal®), Sertindole (Serdolect®, Serlect®), Trimipramine (Surmontil®), Ziprasidone (Geodon®, Zeldox®), and Zotepine (Lodopin®, Losizopilon®, Nipolept®, Setous®). Accordingly, the methods contemplated herein include methods for diagnosing, prognosing, and treating patients with treatment resistant schizophrenia (TRS) where the patient has been observed to be resistant to treatment with any of the aforementioned APDs.
Patients encompassed by the disclosed methods of diagnosis, prognosis, and treatment may include patients having treatment resistant schizophrenia (TRS), where TRS may be defined using generally accepted criteria in the field, for example, as developed by Kane et al. (1988) in the Multicenter Clozapine Trial (MCT).8 These criteria included: 1) persistent positive psychotic symptoms: item score≧4 (moderate) on at least 2 of 4 positive symptom items on the Brief Psychiatric Rating Scale (BPRS) (rated on a 1-7 scale)—hallucinatory behavior, suspiciousness, unusual thought content, and conceptual disorganization; 2) current presence of at least moderately severe illness as rated by the total BPRS score (score≧45 on the 18-item scale) and a score of ≧4 (moderate) on the Clinical Global Impression (CGI) scale (Guy 1976)9; 3) persistence of illness: no stable period of good social and/or occupational functioning within the last 5 years (inability to maintain work and relationships); 4) drug-refractory condition: at least 3 periods of treatment in the preceding 5 years with conventional APDs (from at least 2 chemical classes) at doses≧1000 mg per day of chlorpromazine for 6 weeks, each without significant symptom relief, and failure to improve by at least 20% in total BPRS score or intolerance to a 6-week prospective trial of haloperidol at a dose of 10-60 mg per day. From an operational point of view, only the positive symptom severity, number and duration of clinical trials (two) are utilized for clinical or research purposes. Treatment resistant schizophrenia usually does not respond to increased dosages of APDs or adding adjunctive agents such as benzodiazepines, antidepressants, anticonvulsants or lithium carbonate.
As used herein, a “caregiver” means a person who administers care to a patient. A “caregiver” may include, among others, a physician, physician's assistant, and a nurse. In the disclosed methods, a caregiver may receive the results of a genetic test prior to administering treatment to a patient. For example, a caregiver may receive the results of a test indicating whether the patient possesses a polymorphism, and the caregiver may administer treatment to the patient based on the results of the test. The test indicating whether the patient possesses a polymorphism may be performed by the caregiver and/or the test may be performed another person or testing facility.
As used herein the terms “diagnose” or “diagnosis” or “diagnosing” refer to distinguishing or identifying a disease, syndrome or condition or distinguishing or identifying a person having or at risk for developing a particular disease, syndrome or condition. As used herein the terms “prognose” or “prognosis” or “prognosing” refer to predicting an outcome of a disease, syndrome or condition. The methods contemplated herein include diagnosing a treatment resistant psychiatric disorder in a patient (e.g., diagnosing TRS in a patient). The methods contemplated herein also include determining a prognosis for a patient having a psychiatric disorder (e.g., determining a prognosis for a patient having TRS).
As used herein, the terms “treating” or “to treat” each mean to alleviate symptoms, eliminate the causation of resultant symptoms either on a temporary or permanent basis, and/or to prevent or slow the appearance or to reverse the progression or severity of resultant symptoms of the named disease or disorder. As such, the methods disclosed herein encompass both therapeutic and prophylactic administration. In particular, the methods contemplated herein include treating a patient having or at risk for developing a treatment resistant psychiatric disorder (e.g. treating a patient having TRS).
The present methods may include detecting a SNP in a patient sample. The term “sample” or “patient sample” is meant to include biological samples such as tissues and bodily fluids. “Bodily fluids” may include, but are not limited to, blood, serum, plasma, saliva, cerebral spinal fluid, pleural fluid, tears, lactal duct fluid, lymph, sputum, and semen. A sample may include nucleic acid, protein, or both.
The detected SNP is present in nucleic acid. The term “nucleic acid” or “nucleic acid sequence” refers to an oligonucleotide, nucleotide or polynucleotide, and fragments or portions thereof, which may be single or double stranded, and represents the sense or antisense strand. A nucleic acid may include DNA or RNA, and may be of natural or synthetic origin. For example, a nucleic acid may include mRNA or cDNA. Nucleic acid may include nucleic acid that has been amplified (e.g., using polymerase chain reaction). Nucleic acid may include genomic nucleic acid.
As used herein, the term “assay” or “assaying” means qualitative or quantitative analysis or testing.
As used herein the term “sequencing,” as in determining the sequence of a polynucleotide, refers to methods that determine the base identity at multiple base positions or that determine the base identity at a single position.
The term “amplification” or “amplifying” refers to the production of additional copies of a nucleic acid sequence. Amplification is generally carried out using polymerase chain reaction (PCR) technologies known in the art.
The term “oligonucleotide” is understood to be a molecule that has a sequence of bases on a backbone comprised mainly of identical monomer units at defined intervals. The bases are arranged on the backbone in such a way that they can enter into a bond with a nucleic acid having a sequence of bases that are complementary to the bases of the oligonucleotide. The most common oligonucleotides have a backbone of sugar phosphate units. Oligonucleotides of the method which function as primers or probes are generally at least about 10-15 nucleotides long and more preferably at least about 15 to 25 nucleotides long, although shorter or longer oligonucleotides may be used in the method. The exact size will depend on many factors, which in turn depend on the ultimate function or use of the oligonucleotide. An oligonucleotide (e.g., a probe or a primer) that is specific for a target nucleic acid will “hybridize” to the target nucleic acid under suitable conditions. As used herein, “hybridization” or “hybridizing” refers to the process by which an oligonucleotide single strand anneals with a complementary strand through base pairing under defined hybridization conditions. Oligonucleotides used as primers or probes for specifically amplifying (i.e., amplifying a particular target nucleic acid sequence) or specifically detecting (i.e., detecting a particular target nucleic acid sequence) a target nucleic acid generally are capable of specifically hybridizing to the target nucleic acid.
The present methods and kits may utilize or contain primers, probes, or both. The term “primer” refers to an oligonucleotide that hybridizes to a target nucleic acid and is capable of acting as a point of initiation of synthesis when placed under conditions in which primer extension is initiated (e.g., primer extension associated with an application such as PCR). For example, primers contemplated herein may hybridize to one or more polynucleotide sequences comprising the SNPs disclosed herein. A “probe” refers to an oligonucleotide that interacts with a target nucleic acid via hybridization. A primer or probe may be fully complementary to a target nucleic acid sequence or partially complementary. The level of complementarity will depend on many factors based, in general, on the function of the primer or probe. For example, probes contemplated herein may hybridize to one or more polynucleotide sequences comprising the SNPs disclosed herein. A primer or probe may specifically hybridize to a target nucleic acid (e.g., hybridize under stringent conditions as discussed herein). In particular, primers and probes contemplated herein may hybridize specifically to one or more polynucleotide sequences that comprise the SNPs disclosed herein and may be utilized to distinguish a polynucleotide sequence comprising a minor allele from a polynucleotide sequence comprising the major allele.
An “oligonucleotide array” refers to a substrate comprising a plurality of oligonucleotide primers or probes. The arrays contemplated herein may be used to detect the SNPs disclosed herein.
As used herein, the term “specific hybridization” indicates that two nucleic acid sequences share a high degree of complementarity. Specific hybridization complexes form under stringent annealing conditions and remain hybridized after any subsequent washing steps. Stringent conditions for annealing of nucleic acid sequences are routinely determinable by one of ordinary skill in the art and may occur, for example, at 65° C. in the presence of about 6×SSC. Stringency of hybridization may be expressed, in part, with reference to the temperature under which the wash steps are carried out. Such temperatures are typically selected to be about 5° C. to 20° C. lower than the thermal melting point (Tm) for the specific sequence at a defined ionic strength and pH. The Tm is the temperature (under defined ionic strength and pH) at which 50% of the target sequence hybridizes to a perfectly matched probe. Equations for calculating Tm and conditions for nucleic acid hybridization are known in the art.
As used herein, a “target nucleic acid” refers to a nucleic acid molecule containing a sequence that has at least partial complementarity with a probe oligonucleotide, a primer oligonucleotide, or both. A primer or probe may specifically hybridize to a target nucleic acid.
The present methods may be performed to detect the presence or absence of the disclosed SNPs. Methods of determining the presence or absence of a SNP may include a variety of steps known in the art, including one or more of the following steps: reverse transcribing mRNA that comprises the SNP to cDNA, amplifying nucleic acid that comprises the SNP (e.g., amplifying genomic DNA that comprises the SNP), hybridizing a probe or a primer to nucleic acid that comprises the SNP (e.g., hybridizing a probe to mRNA, cDNA, or amplified genomic DNA that comprises the SNP), and sequencing nucleic acid that comprises the SNP (e.g., sequencing cDNA, genomic DNA, or amplified DNA that comprises the SNP).
A “polymorphism” refers to the occurrence of two or more alternative genomic sequences or alleles between or among different genomes or individuals. “Polymorphic” refers to the condition in which two or more variants of a specific genomic sequence can be found in a population. A “polymorphic site” is the locus at which the variation occurs. A single nucleotide polymorphism (SNP) is the replacement of one nucleotide by another nucleotide at the polymorphic site. Deletion of a single nucleotide or insertion of a single nucleotide also gives rise to single nucleotide polymorphisms. “Single nucleotide polymorphism” preferably refers to a single nucleotide substitution. Typically, between different individuals, the polymorphic site can be occupied by two different nucleotides which results in two different alleles with the most common allele in the population (i.e., the ancestral allele) being referred to as the “major allele” and the less common allele in the population being referred to as the “minor allele.” An individual may be homozygous or heterozygous for the major allele or the minor allele of the polymorphism. “Mutation” as utilized herein, is intended to encompass a single nucleotide substitution, which may be recognized as a single nucleotide polymorphism.
The present methods may be performed to determine whether a patient is heterozygous or homozygous for a major or minor allele of a polymorphism (e.g., a SNP). The term “heterozygous” refers to having different alleles at one or more genetic loci in homologous chromosome segments. As used herein “heterozygous” may also refer to a sample, a cell, a cell population or a patient in which different alleles (e.g., major or minor alleles of SNPs) at one or more genetic loci may be detected. Heterozygous samples may also be determined via methods known in the art such as, for example, nucleic acid sequencing. For example, if a sequencing electropherogram shows two peaks at a single locus and both peaks are roughly the same size, the sample may be characterized as heterozygous. Or, if one peak is smaller than another, but is at least about 25% the size of the larger peak, the sample may be characterized as heterozygous. In some embodiments, the smaller peak is at least about 15% of the larger peak. In other embodiments, the smaller peak is at least about 10% of the larger peak. In other embodiments, the smaller peak is at least about 5% of the larger peak. In other embodiments, a minimal amount of the smaller peak is detected.
As used herein, the term “homozygous” refers to having identical alleles (e.g., major or minor alleles of SNPs) at one or more genetic loci in homologous chromosome segments. “Homozygous” may also refer to a sample, a cell, a cell population, or a patient in which the same alleles at one or more genetic loci may be detected. Homozygous samples may be determined via methods known in the art, such as, for example, nucleic acid sequencing. For example, if a sequencing electropherogram shows a single peak at a particular locus, the sample may be termed “homozygous” with respect to that locus.
Suitable SNPs for the presently disclosed methods, kits, and arrays are disclosed, for example in Table 2. The chromosome location and position of the SNP are listed in Table 2. For example, where the SNP referred to by dbSNP reference ID No. rs2237457 is located on chromosome 7, position 50693638, which may be a thymidine (T) in the minor allele, or a cytosine (C) in the major allele.
The present methods contemplate detecting a single nucleotide polymorphism (SNP). For example, the present methods may detect rs2237457 in either one or both alleles of the patient. (See rs2237457 SNP entry at the National Center for Biotechnology Information, which entry is incorporated herein by reference and refers to a C←→T transition at the reference nucleotide position, where the C-allele is the major allele and the T-allele is the minor allele). This SNP previously was assessed in regard to its association with type 2 diabetes.31 The present methods may detect a T-allele or a C-allele corresponding to the polymorphism (i.e., a T-nucleotide or a C-nucleotide at the position associated with the rs2237457). The present methods may detect whether a patient is homozygous or heterozygous for a T-allele or C-allele (i.e., whether the patient is TT, TC, or CC at the reference nucleotide position for rs2237457).
Specific genetic markers identified in the methods may include the polymorphism rs2237457. Based on the polymorphism being identified in the patient, for example where the patient is homozygous for the T-allele of rs2237457, the patient may be diagnosed as having a treatment resistant psychiatric disorder. As such, the patient may be administered a treatment for a treatment resistant psychiatric disorder, which may include but is not limited to one or more of the following treatments: foregoing administering to the patient a typical or atypical APD, administering to the patient an APD that is approved for treating a treatment resistant psychiatric disorder such as clozapine (e.g., Clozaril® and FazaClo®), and/or administering a new APD having an active ingredient that is a trace amine-associated receptor 1 (TAAR1) agonist.
The present methods may detect the polymorphism directly by analyzing chromosomal nucleic acid having the polymorphic variant sequence. Alternatively, the present method may detect the polymorphism indirectly by analyzing an isoform nucleic acid expressed from the polymorphic variant sequence, an isoform polypeptide expressed from the polymorphic variant sequence, or the expression of another nucleic acid or protein whose expression is regulated by the polymorphic sequence. For example, the present methods contemplate detecting the single nucleotide polymorphism (SNP) rs2237457 (TT, TC, or CC), where the ancestral allele is the C-allele. This reference SNP is located in the gene for growth factor receptor-bound protein 10 (GRB10). As disclosed herein, a 30 kb block flanking rs2237457 in GRB10, which is 70 kb upstream of L-Dopa Decarboxylase (DCC), was identified as an exclusive cis-acting eQTL for DCC in human dorsolateral prefrontal cortex by BrainCloud. Accordingly, the present methods contemplate detecting expression levels of DCC in a patient, and based on detecting reduced expression levels of DCC, administering an agent to the patient that increases expression of DCC. The present methods also contemplate detecting biological activity levels of DCC in a patient, and based on detecting reduced biological activity levels of DCC, administering an agent to the patient that increases biological activity of DCC. Because DCC is an enzyme that is rate-limiting in the synthesis of trace amines and neurotransmitters implicated in psychiatric disorders, the present methods also contemplate detecting biological activity levels of DCC in a patient, and based on detecting reduced biological activity levels of DCC, administering an agent to the patient that is an agonist of the trace amine-associated receptor 1 (TAAR1) for which trace amines are a ligand.
The results disclosed herein implicate L-Dopa Decarboxylase in psychiatric disorders, in particular, treatment resistant psychiatric disorders such as treatment resistant schizophrenia. As such, the presently disclosed methods include methods of treating a patient having or at risk for developing a psychiatric disorder (e.g., treatment resistant psychiatric disorders such as such as treatment resistant schizophrenia) by administering to the patient an agent that modulates L-Dopa Decarboxylase expression or activity (e.g., an agent that increases L-Dopa Decarboxylase expression or activity). Because DCC is an enzyme that is rate-limiting in the synthesis of trace amines and neurotransmitters implicated in psychiatric disorders, the present methods also contemplate treating a patient having or at risk for developing a psychiatric disorder (e.g., a treatment resistant psychiatric disorders such as such as treatment resistant schizophrenia) by administering to the patient an agent that is an agonist of the trace amine-associated receptor 1 (TAAR1) for which trace amines are a ligand.
Regarding TAAR1 agonists, the TAAR1 is a G protein-coupled receptor that in humans is encoded by the TAAR1 gene and is characterized by its ability to bind low-concentration, endogenous monoamines called trace amines and effect increases in cyclic AMP concentrations in the cell. (See Miller, J. Neurochem. 116(2):164-176 (2010); and Lindemann et al. Genomics 85(3):372-85 (2004)) Trace amines typically are found at low concentrations (e.g., 0.1-10 nM concentrations) and constitute less than 1% of total biogenic amines in the mammalian nervous system. (See Zucchi et al., Br. J. Pharmacol. 149(8):967-78 (2006)). These endogenous trace amines include para/meta-tyramine, tryptamine, phenylethylamine, and para/meta-octopamine, which share structural similarities with the three common monoamines: serotonin, dopamine, and norepinephrine. Each ligand has a different potency for TAAR1 (i.e. EC50), which may be measured by the trace amine's ability to effect increases in cyclic AMP (cAMP) concentration after the trace amine binds the receptor. In the presently disclosed methods, the agonists for TAAR1 preferably have an EC50 for the TAAR1 of less than about 1000 nM, 500 nM, 400 nM, 300 nM, 200 nM, or 100 nM. Suitable TAAR1 agonists for use in the disclosed methods may include, but are not limited to, m-tyramine, p-tyramine, phenethylamine, m-octopamine, p-octopamine, N-methylphenethylamine, N-methyltyramine, phenylethanolamine, synephrine, 3-methoxytyramine, 3-iodothyronamine, tryptamine, N-methyltryptamine, dopamine, norepinephrine, serotonin, histamine, amphetamine, methamphetamine, 3,4-methylenedioxymethamphetamine, 2,5-dimethoxy-4-iodoamphetamine, phenylpropanolamine, alpha-methylphenethylamine, beta-methylphenethylamine, N-methylphenethylamine, 2-methylphenethylamine, 3-methylphenethylamine, 4-methylphenethylamine, (S)-4-[(ethylphenylamino)methyl]-4,5-dihydrooxazol-2-ylamine (RO5166017), (4S)-4-[(2S)-2-phenylbutyl]-4,5-dihydro-1,3-oxazol-2-amine (RO5256390), (4S)-4-[3-fluoro-2methylphenyl]-4,5-dihydro-1,3-oxazol-2-amine (RO5263397), and (4S)-4-[3,4-dichlorophenyl]-4,5-dihydro-1,3-oxazol-2-amine (RO5203648).
Also disclosed herein are kits for performing the disclosed methods. For example, the disclosed kits may include reagents for diagnosing, prognosing, and/or treating a psychiatric disease or disorder in a patient. The presently disclosed kits may include specific reagents such as one or more of: (a) reagents for detecting the genotype of a patient at the polymorphism rs2237457 (e.g., reagents for sequencing nucleic acid comprising rs2237457 and/or reagents for probing nucleic acid comprising rs2237457); and (b) a pharmaceutical agent for treating a psychiatric disease or disorder. The reagents in the kit may include nucleic acid reagents (e.g., primers and/or probes) and non-nucleic acid reagents (e.g., polymerases and buffers). The pharmaceutical agent of the kits may include an atypical drug for treating the psychiatric disease or disorder such as a TAAR1 agonist formulated for administration to the patient or clozapine formulated for administration to the patient.
The following examples are illustrative and are not intended to limit the claimed subject matter.
A Genetic Locus in 7p12.2 Associated with Treatment Resistant Schizophrenia
Abstract
Identifying biomarkers for treatment resistant schizophrenia (TRS) could assist decision making regarding antipsychotic drug (APD) treatment and clarification of the homogeneity of TRS. We conducted a genome-wide association study (GWAS) of Caucasian patients with schizophrenia (79 TR and 95 non-TR). Rs2237457 had the lowest p value in an allelic association test (p=5.53×10-6); it is located in 7p12.2 with reported imprinted inheritance. Haploview disclosed a 30 kb block flanking this SNP in GRB100, 70 kb upstream of L-Dopa Decarboxylase (DDC), an enzyme which is rate limiting in the synthesis of trace amines and neurotransmitters implicated in schizophrenia. The association with TRS improved to p<10−11 after excluding heterozygotes, justified by unclear inheritance of the risk allele. This SNP or haplotype was identified as an exclusive cis-acting eQTL for DDC in human dorsolateral prefrontal cortex by BrainCloud. A replication sample produced a weaker result, but in the same direction. Rs2237457 may provide an initial screening test to identify Caucasian TRS and implicates trace amines in the pathophysiology of TRS.
1. Introduction
It is estimated that 1% of the adult population meets current criteria for schizophrenia. Following treatment with typical or atypical antipsychotic drugs (APDs), about 70% of people with schizophrenia have significant reductions in delusions and hallucinations, core psychotic symptoms. The other 30% of schizophrenia patients are referred to as treatment resistant (TRS) or refractory schizophrenia.1-4 TRS is operationally defined as persistence of moderate to severe positive symptoms despite two or more trials of 4-6 weeks duration with typical or atypical APDs other than clozapine, the only drug approved for TRS.5,6 Although most TRS patients are cognitively impaired, and have mild to severe negative symptoms, poor outcome, or unresponsive negative symptoms and cognitive impairment, these features are not the basis for TRS as defined here. The frequency and severity of these dimensions of schizophrenia are similar in TRS and non-TRS schizophrenia.7
Excluding inadequate APD trials, TRS may result from pharmacodynamic or pharmacokinetic reasons. The latter are rare and will not be considered here. Since schizophrenia is a heterogeneous, complex disease, differences in the cause of psychosis could be the basis for difference in drug response, suggesting TRS is also heterogeneous, specifically due to various genetic or epigenetic factors. During the past decade, there has been a dramatic increase in GWAS studies in schizophrenia to identify risk genes, including genes associated with specific phenotypes and endophenotypes. This approach permits interrogation of the entire human genome at high resolution and high throughput in small to large number of unrelated individuals, unconstrained by prior hypotheses, as is the case with a candidate gene approach. The latter may lead to an incomplete (biased) identification of the biologic pathways relating genes to the trait.
Previous GWAS or candidate gene studies have sought to identify genetic marker(s) of treatment resistance. SNPs, structural variants, or haplotypes in DISC1,8 UCP4(SLC25A27),9 MCP-1(SCYA2),10 CYP3A5,11 DRD3,11 HTR2A,12 HTR3A,13 HTR3B,14 SLAMF1,15 NFKB1,15 FAS,16 BDNF,17 and others have been reported to be associated with TRS. The reported effect sizes are generally small and have yet to be replicated. Patients with TRS have been reported to have increased cortical atrophy on magnetic resonance imaging (MRI) compared with patients who respond to treatment.18,19 Cortical atrophy may be due to deficits in neurogenesis, synaptic plasticity, and neurodegeneration. Genes which encode proteins that affect neurogenesis, synaptogenesis, axonal guidance, dendritic spine formation, and synaptic plasticity have been shown to be risk genes for schizophrenia20,21 and could also be candidate genes for TRS.
In this study, GWAS was initially used to identify genetic markers for TRS in Caucasians, based upon persistence of positive symptoms despite two or more adequate trials. A specific genomic locus or haplotype in intron 4 of growth factor receptor-bound protein 10 (GRB10) was found to be associated with TRS, stimulating an expanded follow-up analysis. BrainCloud, a well-established gene expression database based upon dorsal lateral prefrontal cortex (DLPFC) from the general population, was used to test whether identified genetic markers were related to gene expression changes.22 More information is available at the Lieber Institute for Brain Development's website.
2. Methods and Materials
2.1 Subjects and Phenotypes
Written informed consent was obtained from all patients after a full description of the risks and benefits of participation. The GWAS included 174 self-described Caucasian patients diagnosed with schizophrenia or schizoaffective disorder by DSM-IV criteria. They were recruited between 1999 and 2010 at clinical facilities associated with Departments of Psychiatry at Case Western Reserve University and Vanderbilt University. The subjects selected for GWAS had participated in prospective clinical trial or longitudinal study of the effect of clozapine. An additional 195 subjects with these diagnoses who had been prospectively classified as TR or NTR were identified from our sample bank for genotyping to obtain additional data. Nearly all of these patients had participated in prospective clinical trials. All study procedures were approved by the Institute Review Boards of the two universities. The data used for classification as TR or NTR, as well as diagnosis, were based upon research interviews supervised and reviewed by one of us (HYM), as well as all available medical records and interviews with significant others when possible.5,6
2.2 Genotyping
Genome-wide SNP genotyping was performed using Illumina 610K quad BeadChip® (Illumina, San Diego, Calif., USA) by the Center for Human Genome Variation at Duke University (Durham, N.C., USA). Genotype calling was determined by Illumina GenomeStudio® software (Illumina) using default parameters. Taqman® assay for the top hits was performed on an ABI 7900HT Sequence Detection System (Applied Biosystems, Foster City, Calif., USA) at Northwestern University Genomic Core Facility (Chicago, Ill., USA). Taqman Genotyper Software (Applied Biosystems) was used to do allelic discrimination end-point analysis. Concordance rates between Taqman assay and Illumina GWAS data for the selected top five SNPs was 100%.
2.3 GWAS Data Analysis
All GWAS quality control and association testing was conducted with PLINK 1.0.7 software23 (more information available at the PLINK website). In addition to the individual SNP, haplotype-based association testing was also conducted. Other statistical analysis was performed by SPSS. Genotype-phenotype association study was carried out by comparing the frequencies of alleles or genotypes between the TR and NTR groups using the following methods: allelic and genotypic association by Pearson's Chi-square, Cochran-Armitage Trend test, or Logistic regression with multiple covariates. Computer corrected significance values (FDR-BH) were calculated. Linkage disequilibrium (LD) was determined with Haploview (more information available at the Broad Institute's website. SNP imputation was achieved with default settings using PUNK binary file sets of the Hapmap genotype data (release 22) available for download at the PLINK website. The same HapMap CEU reference data (2,543,887 SNPs) were used for all imputations. SNPs with A/T or CJG transitions which have a potential to go unflagged by flip command are eliminated. A combined imputation and association analysis of case/control study was performed. SNP Annotation and Proxy Search (SNAP) (more information available at the Broad Institute's website) was also used to search SNPs in LD from 1000 Genome Database or Hapmap database. QUANTO 1.2 was used to calculate the power of the test more information available from the University of Southern California's website).
2.4 Mapping Cis eQTL and Cis Regulatory Elements
BrainCloud provides genome-wide gene expression data and genetic control based upon a study of human postmortem DLPFC from 109 postnatal Caucasian subjects from prenatal period and throughout the life span.22 We included all SNP loci within 200 kb interval (100 kb up- and downstream) of transcripts/genes which were found to be associated with TR. ENCODE (more information available at the University of California-Santa Cruz's website and UCSF Brain Methylation Database24 were applied for prediction of identified genetic variants function as cis regulatory elements.
3. Results
3.1 Clinical Characteristics and Demographic Feature of TR Vs NTR.
Clinical characteristics and demographic feature for the GWAS sample: TR (n=79) and NTR (n=95) patients, and for the additional subjects who were genotyped (TR: n=70 and NTR: n=125) are reported in Table 1.
(GWAS/Additional)
#t-test;
$ANCOVA test adjusted by Gender.
There was no significant difference between the four groups in gender. Age at onset was nearly significantly earlier in the TR patients, consistent with previous studies.25 There were significantly more NTR than TR patients who were schizoaffective (chi-square=9.554, p=0.002 in the GWAS sample and the replication sample; chi-square=5.595, p=0.018).
3.2 Association Study
SNPs which differentiated TR and NTR schizophrenia by allelic association test were identified by p values and were readily distinguished by Manhattan plot and Q-Q plot (see
The Cochran-Armitage trend test was also performed to test for association. Rs2237457, a GRB10 SNP, one of the top bits associated with TR, resides in a well-known imprinted genomic region, 7p12.2 (GRB10-DDC), in human and mouse (more information available at GeneImprint's website). There is extensive evidence that brain-specific paternal expression of GRB10 regulates both fetal brain development26 and subsequent ‘social dominance’ behavior and self-grooming.27 Menheniott also reported genomic imprinting of DDC in heart and reciprocal allelic expression with neighboring GRB10 in mouse.28 Therefore, due to unclear parental origin of risk allele, we reanalyze the GWAS dataset by excluding the 71 cases with heterozygous genotype for rs2237457, leaving 53 NTR cases and 50 TR cases. Interestingly, the association with TR was significantly enhanced, from p 5.53×10−6 with the heterozygotes to p=4.99×10−12, by allelic association test, and from p=1.66×105 to p=1.04×10−6 by the Cochran-Armitage Trend test. With the allelic association test, rs2237457 survived correction for multiple testing with adjusted p=2.74×10−6 by Bonferroni correction and adjusted p=1.37×10−6 by FDR-BH correction. There were no other SNPs across the whole genome close to the same level of significance as rs2237457, or SNPs in LD with rs2237457, in a 30 kb region of GRB10. The rs2237457 genotype was used to predict TR status of 103 schizophrenia patients with homozygous genotypes. The sensitivity was 48%; the specificity was 94.34%; the positive predictive value was 88.89%; the negative predictive value was 65.79%; the positive likelihood ratio was 8.48; the negative likelihood ratio was 0.55. Thus, this method identified 48% patients with TRS, while 52% diagnosed as TRS were not identified by this allele, suggesting that at least half of TRS cases are due to other genetic variants or non-genetic influences. As will be discussed, the TRS cases identified by rs2237457 could be a subtype of TRS.
Logistic regression was performed for the association of TRS with rs2237457, after controlling for age of onset, gender, and PCA1/2 of population stratification (p=7.64×10−5 for all cases, p=6.13×10−5 for homozygous cases). None of the covariates had a significant impact on this association (p>0.05). QUANTO 1.2 was used to calculate the power of the test. We used a disease model with unmatched case-control ratio of 1:1.2. The mode of inheritance was log-additive. Kp=0.30 (population risk) since the prevalence of TRS in the schizophrenia population is 30%. The T allele of rs2237457 was found to have a frequency of 0.36. If each additional T allele increases the odds of having TRS by 2.60, and 79 TR cases and 95 NTR controls were genotyped, the power (chance) to detect an association with significance p<0.01 was 95%.
3.3 Functional Analysis of the Top Hits
Regional association plots of p values for SNPs flanking the six genomic loci highly associated with TR were performed. A highly LD haplotype region of rs10825879 was identified as associated with TR. This 390.5 kb region is between ZWINT and MIR3924. ZWINT interacts with SNAP25 and rab3c,29 both of which have been reported to be involved in the regulation of release of neurotransmitters from presynaptic sites by synaptic vesicle exocytosis.
Another top SNP was rs1421321, located in the intron region of Cnpy1. Cnpy1 is expressed in the midbrain-hindbrain boundary in zebrafish, binds fibroblast growth factor receptor 1 (FGFR1), and plays a role in fibroblast growth factor (FGF) signaling.30 FGF controls the growth and patterning of several brain structures, while later in life, they have been reported to regulate synaptogenesis, axonal growth, neuroprotection, learning, memory, and the maintenance and repair of neuronal tissue. Converging evidence supports a possible role in the pathophysiology of schizophrenia.31
EDEM1 is the closest gene to rs6763175. It protects neural cells from mitochondrial damage by paraquat32 and has also been linked to schizophrenia in the GAIN study.33 Another top SNP, rs6924295, was discovered in an intronic region of CYP39A1, which belongs to the category of drug metabolism. This SNP is ˜15 kb upstream of SLC25A27, also called neuronal uncoupling protein 4 (UCP4). Previous studies linked a haplotype (rs3757241-rs10807344-rs9395206-rs2270450) in SLC25A27 to ultra-resistant schizophrenia (refractory to clozapine) in Caucasians, suggesting the treatment resistant may be related to mitochondrial dysfunction and oxidative stress.9 Although this haplotype or individual SNP is not in LD with rs6924295, our data confirmed this haplotype associated with TR (p=0.036) and ultra-resistant TR (p=0.03) with borderline significance (Table 3A and Table 3B).
The strongest TRS association signal was observed at chromosome 7 in a functionally relevant type 2 diabetes candidate gene, GRB10. GRB10 has been shown to bind to the activated insulin receptor and acts as a negative regulator of insulin action and glucose uptake.34 Overexpression of GRB10 in mice causes postnatal growth retardation and insulin resistance. Insulin resistance has been reported to be present in never-medicated, first episode schizophrenia patients in some35 but not all studies.36 This Illumina 610k SNP array contained a total of 15 SNPs in intron 4 of GRB10, 6 of which were significantly associated with TRS (P<10−5). SNP imputation data showed 20 more SNPs in this region of intron 4 associated with TRS at a similar level of significance. All these SNPs were in high LD (r square>0.8), which lead to a conservative 30 kb region in intron 4 of GRB10. Rs2237457 provided the lowest p value for association with TRS. Previous studies showed that this SNP was strongly associated with a Type 2 diabetes trait, OGTT GAUC, from a GWAS study in an old Amish population.37 T is the risk allele for Type 2 diabetes and is also the risk allele for TRS. A GWAS from the Framingham Heart Study 100k project demonstrated that this SNP was associated with BMI in a total of 1341 participants from 310 families having adiposity traits measured over 30 years of follow up.38 That study suggest that rs2237457 is one of the top 500 SNPs associated with adiposity traits (−log10 p value=4.3). We performed a test of SNP (rs2237457) by covariate interaction through introducing an interaction term in the logistic regression model as follows: environmental factor included gender, age of onset, PCA1/2, and BMI, Y=β0+β1.ADD+P2.COV1÷β3.ADD×COV1+e. For example, COV1 was the BMI of the patients (TR and NTR). The p value for β3 is 0.98 for all cases and 0.56 for cases with homozygous genotypes, suggesting that there was no significant interaction between rs2237457 and BMI. Other environmental factors alone or in combination demonstrated no interaction with rs2237457 genotype (p>0.01) in association with TR. A summary of the metabolic data is presented in Table 4.
One third of the patients had fasting blood glucose and HgA1C level data. No significant interactions between rs2237457 and fasting glucose or HgA1C was found. ANCOVA was also used to evaluate the genotypic effect of rs2237457 on baseline metabolic measures, including BMI, HgA1C, fasting blood glucose, high-density lipoprotein (HDL), low-density lipoprotein (LDL), triglyceride (TG), TG/HDL, and total cholesterol levels, using gender as a covariate. ANCOVA showed no significant interaction between gender and metabolic measures. The main effect of rs2237457 genotype (TT, TC, CC) on all metabolic measures was not statistically significant (p>0.05). Hence, rs2237457's contribution to the TRS trait was independent of adiposity or diabetes traits.
We also noted that this 30 kb haplotype resides in 70 kb upstream of the closest gene, dopa decarboxylase (DDC). DDC, also known as aromatic L-amino acid decarboxylase (AADC), catalyzes the decarboxylation of L-3,4-dihydroxyphenylalanine (DOPA) to dopamine (DA), L-5-hydroxytryptophan to serotonin (5-HT), and L-tryptophan to tryptamine. There is extensive information related DDC to schizophrenia.39-44 This will be discussed subsequently.
3.4 cis eQTL Analysis with Braincloud
The cis eQTL analysis included all SNP loci available within a 200 kb interval (100 kb up and downstream of rs2237457) with transcripts in the BrainCloud database. Single expression probes selected for DDC or GRB10 were located in the 3′ end of mRNA. For the cis association analysis, there were 105 and 130 SNPs within the 200 kb interval for DDC and GRB10, respectively. As shown in
3.5 Joint Analysis with Additional Samples
Additional Caucasians samples (70 TR and 125 NTR) were available from two research clinics of HYM. DNA from these patients was genotyped for rs2237457 by Taqman assay. No cases with genomic DNA were excluded. There was a trend for rs2237457 to be more common in TRS in this sample, but less significant than the GWAS sample (Table 5).
#Test includes Cochran-Armitage trend test (TREND), Genotypic (2 df) test (GENO), Dominant gene action (1 df) test (DOM), Recessive gene action (1 df) test (REC) in addition to the basic allelic test (ALLELIC).
Nevertheless, after the two sample were combined (TR=149; NTR=220), as shown in Table 5, rs2237457 remained significantly associated with TR. The % TT in the entire TRS sample (66.1%) was significantly higher than that in the NTR sample (15.0%). We calculated that each additional T allele increased the odds of being TRS by 1.74 (odds ratio in additive mode). With genotype data from 149 TRS cases and 220 NTR controls, with significance p<0.01, the power to detect an association was 79%.
4. Discussion
4.1 Subtyping TR Schizophrenia, First-Episode and Late-Onset
The main finding of this study was identification of a possible biomarker for TRS based upon a GWAS: rs2237457 and its flanking 30 kb haplotype inside GRB10, which is upstream of DDC, a gene linked to the synthesis of the neurotransmitters DA, norepinephrine (NE), 5-HT, trace amines, and the action of antipsychotic drugs. Thus, there is a plausible biological connection between TRS and DDC. The association between this marker and IRS status was tested in a second sample which produced a weaker result, but in the same direction. The allelic test of the combined sample of patients homozygous for the risk allele, justified because parental imprinting invalidated the classification of heterozygotes, met genome wide significance, p=9.73×1011. These results also reflect genetic heterogeneity, such that some TRS patients have an etiologically distinct form of TRS. When a common genetic feature such as the locus in GRB10-DDC identified in this study exerts more influence than individual heterogeneity, this shared genetic liability opens up the possibility of genetically-based prediction and refinement in diagnosis. We identified one subgroup (TT) of TRS. The remainder, i.e. all CC and some TC individuals, represent one or more additional types of TRS. Thus, TRS is a complex trait. The definitional lack of response of positive symptoms to APDs in TRS patients is not explained by APD dosage, drug metabolism, pharmacokinetic differences, environmental, or other non-genetic factors, alone or in combination. Therefore, it is reasonable to hypothesize that some TRS patients represent a biologically, and, possibly, genetically distinct group, since there is evidence that positive symptoms in patients with schizophrenia,45 and response of positive symptoms to APD treatment46,47 are likely to be heritable. This is most likely to be the case for first-episode TR patients who never respond to APDs at adequate dosages, as opposed to the subset of treatment resistant patients who initially experience improvement in positive symptoms during treatment with APDs, but subsequently become refractory, so-called late-onset TR. These two types of TRS patients have been previously discussed.43 In this study, we did not stratify the TR cases based on this dimension due to lack of data.
Several lines of evidence suggest a complicated link between type 2 diabetes or impaired glucose metabolism and schizophrenia35,36,48-52 However, in this study, we did not observe a significant impact of rs2237457 genotype on metabolic measures. However, metabolic data was available for only one third of the cases, giving limited power to test the possibility of comorbidity modifying treatment resistance.
4.2 Bioinformatics Analysis on GRB10-DDC Region
Further analysis revealed that this 30 kb intronic region associated with TRS in GRB10 spans a variety of DNase I hypersensitivity clusters, suggesting a strong regulatory activity in this region, and enrichment of the modified histone mark, H3K4Me1, suggesting enhancers. According to UCSF brain methylation database,24 the TRS associated intronic haplotype we identified harbors an H3K4me3 site (promoter). There is a strong overlap of H3K4me3 with an unmethylated intragenic CGI, suggesting this is an alternate promoter site. This promoter site is uniquely neuronal or tissue-specific since it does not overlap with the histone methylation tracks (H3K4Me3) derived from ENCODE which has signal overlaid from nine non-neuronal cell lines from a variety of human tissues. Therefore, we propose that the sequence variation in this genetic region of GRB10 will have an impact on the neuronal specific promoter, resulting in a change of gene expression, or generate isoforms of GRB10. Many common variants associated with disease are highly likely to be expression quantitative trait loci (eQTLs). These considerations suggest that this haplotype could be a cis eQTL for GRB10, or nearby genes such as DDC.
We successfully mapped the TR-associated SNP or haplotype in this region as an exclusive cis eQTL for the target gene, DDC. Although an early genetic study did not suggest that genetic polymorphism (including ˜1 kb upstream promoter region) in DDC gene is associated with liability for developing schizophrenia or bipolar disorder,53 this 30 kb region, ˜70 kb upstream of the transcriptional starting site of DDC, has never been investigated.
4.3 DDC and Trace Amines
DDC is not considered to be the rate-limiting enzyme for catecholamine or indoleamine synthesis, but is rate-limiting for the synthesis of trace amines (TAs) such as 2-phenylethylamine (2-PE), p-tyramine, and tryptamine, which are thought to present at least two orders of magnitude below the level of neurotransmitters such as DA, NE and 5-HT. TAs, through the action of trace amine-associated receptor 1 (TAAR1), are believed to function as neuromodulators to maintain neuronal activity of monoamine neurotransmitters, possibly also glutamate, within defined physiological limit.54-57 The TAAR1 agonists have been shown to have antipsychotic activity in animal models.57 DDC expression and activity is altered by APD administration, including clozapine, in rodents, suggesting it plays a role in response to APD treatment.58 DDC was showed to be regulated by immediate and delayed onset mechanisms involving enzyme activation and gene expression induction.59 Abi-Dargham summarized eight PET studies on striatal DDC activity in patients with schizophrenia using [18F] DOPA or [11C] DOPA.42,60,67 Six of them reported increased DDC activity.42,62-66 As a rate-limiting enzyme for all TAs synthesis, DDC activity could play a critical role in determining response to APDs.
A deficit in DDC activity in a subgroup of TR schizophrenia leading to inadequate stimulation of TAAR1 could contribute to inadequate response to APDs. If so, TAAR1 agonists may be expected to be particularly useful in the treatment of TRS patients. The limitations in this study include a relatively small number of subjects, and the absence of data in ethnic groups other than Caucasians. In summary, this is the first GWAS study of TRS in Caucasians. We have identified a novel genetic signature for some TRS patients and a genetic marker in 7p12.2 region, which may provide an initial screening test to identify some TR Caucasian schizophrenia.
5. Supplemental Material:
5.1 Quality Control of Genotyping Data
Quality control was done by examining a summary statistics of PLINK. Estimation of IBD sharing (relatedness) was done in order to avoid the cases from most recent common ancestor in homogeneous random mating population. All 174 cases were considered as unrelated Individuals based on the pairwise IBD values greater than 0.025 as a cutoff. Only individuals and SNPs meeting the following criteria were included in further analyses. There was no flagged individual with discrepancy on gender by sex check. Call rate for individual had to be higher than 90%. Genotyping rate for each SNP had to be higher than 95%. Minor allele frequency (MAF) for each SNP had to be higher than 1%. SNPs with significant deviation from Hardy-Weinberg Equilibrium (P<10-4) were excluded. After frequency and genotyping pruning, there were 552,419 SNPs left for further analysis. Total genotyping rate in the remaining individuals was 99.9571%.
All patients recruited in this study were self-described Caucasians. Detection of possible population stratification that might influence association analysis was carried out using principal component analysis (PCA) with genotype data for all autosome SNPs passing above QC. The PCA plot in
To search for evidence of systematic bias from unrecognized population structure, analytical approach, and genotyping artifacts, Quantile-quantile (Q-Q) plots were used to examine P-value distributions. The QQ-plot showed no significant deviation from base line, suggesting that population stratification was not significant in this dataset. Genomic inflation factor based on median was 1.097 for self-described Caucasians. There was no significant improvement in this value after excluding the cases (from 1.097 to 1.090), leaving only 166 population-stratified Caucasians in the box of
5.2 SNPs Interaction
Logistic regression conditioned by the target SNP was performed on all SNPs with P<10−5 but adding the allelic dosage of one of those SNPs, for example rs2237457, as a covariate. Top SNP rs2237457 appeared to explain the multiple associated SNPs (e.g. rs2282927, rs6593077, rs2237455, rs2237468, and rs2329487) in that region of GRB10 since no significant association was left after controlling the allelic dosage of rs2237457 (Table 6A).
In addition, we did conditional haplotype testing. Given these six SNPs were in high LD, we explicitly modeled the haplotypic association with TR. Likelihood ratio test suggested that the overall variation at this locus of GRB10 appeared to significantly influence the TR trait., with p=9.10×10−4 (all cases) and p=2.42×10−5 (homozygous cases only) (Table 6B).
There was no independent effect from the six SNPs after testing each individual SNP against haplotypic background (p>0.1 or identical model), suggesting that the weaker association at nearby rs2237457 did not represent an independent signal. When a specific haplotype, CCACCC, with the higher frequency (0.494/0.573), was compared to the rest of other haplotypes, the likelihood ratio test was p=3.64×10−4 (OR=2.09) for all cases and p=3.88×10−6 (OR=3.021) for homozygous cases only. When TACTTT was singled out, the likelihood ratio test gave p=5.82×10−6 (OR=2.83) for all cases and p=4.91×10−6 (OR=9.425) for homozygous cases.
In addition, we noted that other top hits with p<10−5 such as rs10825879, rs6500291, rs6763175, still maintained the most significant associations with TR after adjustment for the allelic dosage of rs2237457, suggesting that rs2237457 alone only contributed part of the genetic association with TR trait (Table 6C).
In other words, combined the genetic risk from all top hits would provide a better marker for genetic prediction of TR. Only the top 32 SNPs (LD prune) with p<10−4 was included in Epistasis analysis with the default parameter. There was no significant epistasis event identified with default p<10−4 (Table 6D).
It will be readily apparent to one skilled in the art that varying substitutions and modifications may be made to the invention disclosed herein without departing from the scope and spirit of the invention. The invention illustratively described herein suitably may be practiced in the absence of any element or elements, limitation or limitations which is not specifically disclosed herein. The terms and expressions which have been employed are used as terms of description and not of limitation, and there is no intention in the use of such terms and expressions of excluding any equivalents of the features shown and described or portions thereof, but it is recognized that various modifications are possible within the scope of the invention. Thus, it should be understood that although the present invention has been illustrated by specific embodiments and optional features, modification and/or variation of the concepts herein disclosed may be resorted to by those skilled in the art, and that such modifications and variations are considered to be within the scope of this invention.
Citations to a number of patent and non-patent references are made herein. The cited references are incorporated by reference herein in their entireties. In the event that there is an inconsistency between a definition of a term in the specification as compared to a definition of the term in a cited reference, the term should be interpreted based on the definition in the specification.
The present application claims the benefit of priority under 35 U.S.C. §119(e) to U.S. Provisional Patent Application No. 61/886,288, filed on Oct. 3, 2013, the content of which is incorporate herein by reference in its entirety.
Number | Date | Country | |
---|---|---|---|
61886288 | Oct 2013 | US |