Truth Seeker
New Member
Abstract
Antipsychotic-induced weight gain has emerged as a serious complication in the treatment of patients with atypical antipsychotic drugs. The cannabinoid receptor 1 (CNR1) is expressed centrally in the hypothalamic region and associated with appetite and satiety, as well as peripherally. An antagonist of CNR1 (rimonabant) has been effective in causing weight loss in obese patients indicating that CNR1 might be important in antipsychotic-induced weight gain. Twenty tag SNPs were analyzed in 183 patients who underwent treatment (with either clozapine, olanzapine, haloperidol, or risperidone) for chronic schizophrenia were evaluated for antipsychotic-induced weight gain for up to 14 weeks. The polymorphism rs806378 was nominally associated with weight gain in patients of European ancestry treated with clozapine or olanzapine. 'T' allele carriers (CT+TT) gained more weight (5.96%), than the CC carriers (2.76%, p=0.008, FDR q-value=0.12). This translated into approximately 2.2kg more weight gain in patients carrying the T allele than the patients homozygous for the CC genotype (CC vs CT+TT, 2.21±4.51 vs 4.33±3.89kg; p=0.022). This was reflected in the allelic analysis (C vs T allele, 3.84 vs 5.83%, p=0.035). We conducted electrophoretic mobility shift assays which showed that the presence of the T allele created a binding site for arylhydrocarbon receptor translocator (ARNT), a member of the basic helix—loop—helix/Per—Arnt—Sim protein family. In this study, we provide evidence that the CNR1 gene may be associated with antipsychotic-induced weight gain in chronic schizophrenia patients. However, these observations were made in a relatively small patient population; therefore these results need to be replicated in larger sample sets.
INTRODUCTION
The use of atypical antipsychotic drugs (AADs) for treatment of schizophrenia has increased and they now capture over 90% of the antipsychotic drug market. Unlike the typical antipsychotic drugs, AADs have the primary advantage of being associated with very low or negligible risk of inducing extrapyramidal symptoms. However, one of the major concerns with AADs is the development of weight gain leading to full metabolic syndrome and early cardiac death (reviewed by Muller and Kennedy, 2006). Obesity in the general population has become a major public health challenge and association of AADs with weight gain adds another dimension to this problem.
The development of antipsychotic induced weight gain in a subset of schizophrenia patients (∼30%), concordance among twins, and high heritability of weight gain suggest a strong genetic component. Two polymorphisms, −759C/T in the serotonin receptor 5-HT2C, and the −2548 G>A promoter polymorphism in the leptin gene, are most consistently associated with antipsychotic-induced weight gain (eg, De Luca et al, 2007; Ellingrod et al, 2007; Templeman et al, 2005; Zhang et al, 2003). Other gene findings have been reviewed elsewhere (Muller and Kennedy, 2006). However, it is more than plausible to assume that, as for other complex inherited phenotypes, numerous and yet undiscovered genetic markers are involved in antipsychotic-induced weight gain.
The cannabinoid receptor 1 (CNR1 or CB1) is a seven transmembrane G-coupled receptor (Gi/o, Gs or Gq/11), expressed widely in several brain regions associated with reward, hunger, appetite and satiety and in peripheral tissues including the gut, liver (hepatocytes), adipocytes, skeletal muscle and pancreas (Ameri et al, 1999; Matias and Di Marzo, 2007; Woods, 2007). CNR1 knockout mice present a lean phenotype (Cota et al, 2003; Di Marzo et al, 2001), reduced food intake and resistance to diet-induced obesity (Osei-Hyiaman et al, 2005; Ravinet Trillou et al, 2003, 2004). Similarly, administration of rimonabant, a CNR1 receptor antagonist, leads to reduced food intake, decreased bodyweight, normalized plasma leptin, insulin, free fatty acids and corrected insulin resistance in wild-type animals (Jamshidi and Taylor, 2001; Ravinet Trillou et al, 2003). Conversely, injection of endogenous ligands of CNR1, such as anandamide or 2-arachidonoylglycerol, into the hypothalamus of pre-satiated rats increases dietary food intake. Furthermore, this increase in food intake could be prevented by pretreatment of the rats with rimonabant (Jamshidi and Taylor, 2001; Kirkham et al, 2002). Weston-Green et al (2008) showed that chronic administration of olanzapine decreases the binding of the agonist, [3H] CP-55940, to CNR1 in the dorsal vagal complex of the brain stem and causes increased weight gain in rats, whereas haloperidol failed to elicit similar effects. In humans, treatment of AIDS and cancer-related anorexia with dronabinol, a CNR1 agonist, leads to significant increase in appetite and body weight compared with placebo (Beal et al, 1995; Nelson et al, 1994). Furthermore, increased CNR1 mRNA levels in visceral adipose tissue of overweight and metabolic syndrome patients have been reported (Sarzani et al, 2009). In addition, hepatic stimulation of the CNR1 receptor leads to increased lipogenesis through augmented expression of the lipogenic transcription factor SREBP-1c (sterol regulatory element-binding protein-1c), an effect inhibited by treatment with rimonabant (Osei-Hyiaman et al, 2005).
At the molecular level, CNR1 interacts with several molecules in the hypothalamus. For example CNR1 downregulates anorexigenic molecules such as cocaine and amphetamine-regulated transcript in the arcuate nucleus and corticotrophin-releasing hormone (CRH) in the paraventricular nucleus. CNR1 also upregulates orexigenic neuropeptide Y and melanin-concentrating hormone (MCH) in the lateral hypothalamus (Cota et al, 2003; Gamber et al, 2005; Osei-Hyiaman et al, 2005). Some of the effects of CNR1 are under regulation of leptin, which prevents the release of orexigenic and anorexigenic molecules (MCH and CRH, respectively; for detailed review see (Kunos, 2007; Matias and Di Marzo, 2007; Woods, 2007). Similarly, ghrelin is known to exert its orexigenic effects at least in part through CNR1 (Kola et al, 2008).
The CNR1 gene is located on 6q14—15, and a common polymorphism in intron 2, rs20232239, has been associated with obesity in the general population (Benzinou et al, 2008) as well as with substance dependence and CNR1 mRNA levels (Zhang et al, 2004). Furthermore, the C allele of this polymorphism has been associated with greater CB1 binding in the prefrontal cortex and alcohol cue-elicited brain activation (Hutchison et al, 2008).
Based on the above studies, the CNR1 gene represents a compelling candidate for genetic studies of antipsychotic-induced weight gain. In this study, we analyzed 20 tag polymorphisms in the CNR1 gene, including the synonymous polymorphism rs1049353 (1359 G>A, Thr453Thr) and obesity-associated rs20232239 for involvement in antipsychotic-induced weight gain.
MATERIALS AND METHODS
Subjects
A total of 183 patients with chronic schizophrenia or schizoaffective disorder were included in this study. Approval from the institutional ethics committee and informed consent were obtained for all patients. The details of the sample are provided in Table 1. Patients were recruited at three different sites: Sample A: Charite University Medicine, Berlin, Germany (DJ Müller and I Puls, n=47). Patients 18—60 years old diagnosed with schizophrenia or schizoaffective disorder according to DSM-IV criteria were included. Other inclusion criteria were a minimum positive and negative symptom scale (PANSS) score of 60 and admission either due to first manifestation or relapse of pre-existent condition with significant deterioration. Patients were given antipsychotic medication and assessed up to 6 weeks (Table 1). Exclusion criteria included pregnancy or breast feeding, organic brain disorder or severe head injuries, severe or unstable comorbid medical conditions (Hepatitis C, HIV, Thyroid disorder, diabetes mellitus), and not stable, substance dependence (polytoxicomania), clinically relevant mental retardation, and severe personality disorder.
Sample B: Case Western Reserve University in Cleveland, Ohio (HY Meltzer, JA Lieberman, n=79). Patients 18—60 years old diagnosed with schizophrenia according to DSMIIR criteria who were either treatment refractory or intolerant to treatment with typical antipsychotics were included. Before starting treatment with clozapine there was a washout period of 7—14 days during which patients received no medications unless it was clinically essential. Clozapine dosage was not fixed and serum levels were monitored to ascertain compliance. Patients were treated with clozapine for 6 weeks and percent change in weight was calculated from baseline weight (for details see (Masellis et al, 1998)). Exclusion criteria were pregnancy, organic brain disorder or severe head injuries, previous medical conditions (Hepatitis C, HIV, Thyroid disorder, diabetes mellitus), which require treatment and are not stable, substance dependence (polytoxicomania), clinically relevant mental retardation and severe personality disorder.
Sample C: Hillside Hospital in Glen Oaks, New York (JA Lieberman, n=57). Patients 18—60 years old, with chronic schizophrenia or schizoaffective disorder according to DSM-IV were included. These patients had shown suboptimal response to previous treatment with typical antipsychotic drugs, mainly defined by persistent positive symptoms and a poor level of functioning over the past 2 years. Patients were randomly assigned to receive either clozapine (500mg/day), olanzapine (20mg/day), risperidone (8mg/day) or haloperidol in a 14-week double-blind study. Patients with a history of non-response to clozapine, risperidone, or olanzapine, defined as an unambiguous lack of improvement despite a contiguous adequate trial of risperidone or olanzapine for at least 6 weeks, or clozapine for at least 14 weeks were excluded. Similarly, patients with a history of clozapine, olanzapine, risperidone, or haloperidol intolerance as well as those who received a depot antipsychotic within 30 days before randomization were also excluded (for details see (Volavka et al, 2002)). Sample B and C have been described in more detail in previous studies (Basile et al, 2002; De Luca et al, 2007; Muller et al, 2005a,2005b).
Ethnicity was determined using ancestry information provided by the participants. Patients were classified for ethnicity by assessing ancestry over last three generations.
Genetic Analysis
10ml of venous blood was obtained from study participants and DNA was extracted at the Centre for Addiction and Mental Health using the high salt method (Lahiri and Nurnberger, 1991). All the genotyping was carried out using TaqMan SNP genotyping assays (Applied Biosystems, Foster City, CA). Genotype calling was confirmed by two independent researchers and 5% of the total sample was re-genotyped to assure concordance. The concordance rate for this study was 99.86% and discordant genotypes were treated as missing in the statistical analysis. A total of 20 tagSNPs covering ∼32,500bp and representing 100% of the common variation in the CNR1 gene region (minor allele frequency >0.05) were included in the study (Figure 1). Among the 20 SNPs rs806381, rs2023239, rs9450902 and rs10485170 were force-included as they were recently reported to be associated with weight gain in the general population (Benzinou et al, 2008).
Data Analysis
All the statistical analyses were performed using SPSS 15.0. The categorical variables were tested using a χ2-test and the continuous variables were tested using ANOVA, ANCOVA, or linear regression. Pearson's correlation was calculated between age, sex, duration of treatment, baseline Brief Psychiatric Rating Scale (BPRS). We decided a priori that the comparison of mean weight change (%) across genotypic categories would be the association-determining test. Only SNPs found to show a trend or significance were used for post hoc analyses. Haplotypes were constructed and percent change in weight was compared using QTPHASE routine in UNPHASED 2.0 (Dudbridge, 2003). Only haplotypes with frequencies greater than 5% were included in the analyses. Linkage disequilibrium (LD) and tagSNPs (r2=0.8, minor allele frequency >0.05) were determined using Haploview 4.1 (Barrett et al, 2005). LD blocks were constructed using solid spine of LD. Power calculations were carried out using Quanto 1.2.4 (Gauderman and Morrison, 2006). Assuming a minor allele frequency of 0.25 in a sample size of n=183, and a mean difference of 2.5% between carriers of the risk genotype, we had more than 80% power in a dominant and additive effect model. False discovery rate (FDR) q-values were calculated using the program QVALITY. Similar to other publications, a q-value threshold 0.20 was used to indicate SNPs with most significant associations (van den Berg et al, 2009).
Electrophoretic Mobility Shift Assay
Biotin-labelled and unlabelled double-stranded oligonucleotides for each allele of rs806378 (T and C) and unlabelled oligonucleotides with the specific consensus-binding site for arylhydrocarbon receptor translocator (ARNT) and XBP1 (Clauss et al, 1996) were synthesized (Table 2, Eurofin MWG Operon, Huntsville, AL). Nuclear extract was prepared from SH-SY5Y neuroblastoma cells (ATCC: CRL2266) using the NE-PER Nuclear and cytoplasmic extraction kit (Pierce, IL). The complementary pairs of nucleotides suspended in STE Buffer (10mM Tris pH 8.0, 50mM NaCl, 1mM EDTA) were annealed using a thermocycler (Step 1: 95°C for 5min, Step 2: 95°C (−1°C/cycle, 70 times), Step 3: HOLD at 4°C). The reannealed double-stranded oligonucleotides (20fmol) were incubated with 2μl of the nuclear extract, 1 × binding buffer, 2.5% glycerol, 5mM MgCl2, 50ng/μl poly (dI dC), 0.05% NP-40 for 20min at room temperature. For the competition experiments, unlabelled oligonucleotides at 200-fold molar excess (4pmol) were preincubated with the nuclear extract before adding the biotin-labelled oligobucleotides. DNA protein complexes were fractioned on a 5% polyacrylamide gel in 0.5 × TBE for 35min at 4°C. Complexes were electro-transferred to nylon membrane at 380mA for 50min, and cross-linked at 120mJ/cm2 for 45s. The biotinylated oligonucleotides were detected by chemiluminescence using the Chemiluminescent nucleic acid detection module of the Light Shift Chemiluminescent EMSA kit (Pierce, IL). The processed nylon membrane was exposed to an X-ray film for 2min.
RESULTS
Genetic Analysis
All the SNPs were in Hardy—Weinberg equilibrium (p>0.05, Table 3, Supplementary Table 1). Several clinical variables were different across the three sites, however, similar weight gain was observed across the clinical sites (Table 1). In the total sample set significant correlation was observed between weight change (%) and duration of treatment (p=0.045) but not with age (p=0.747) and baseline BPRS score (p=0.986). In the subgroup of patients of European ancestry on clozapine or olanzapine no significant correlation of duration of treatment (p=0.532), age (p=0.703), and baseline BPRS score (p=0.878) with weight change (%) was observed. Similarly, weight change (%) was not different between the male and female patients in the total (p=0.776) as well as the sub-group of patients of European ancestry on clozapine or olanzapine (p=0.822).
Association Study in the Total Sample
No association of any of the polymorphisms with weight change (%) was observed when we considered the total sample set of patients on various antipsychotics (p>0.05). Consistent with the literature (eg, reviewed by (Nasrallah, 2008)), we observed highest weight gain rates with clozapine (5.37%±6.27) and olanzapine (7.29%±8.17) in our sample compared with haloperidol (2.76%±4.29) and risperidone (3.54%±6.32). Clozapine and olanzapine have relatively similar receptor-binding profile (Nasrallah, 2008). Therefore, we assessed the subgroup of patients who received either clozapine or olanzapine for association with polymorphisms in CNR1 gene. No association of any of the SNPs with weight change was observed (p>0.05, Table 3, Supplementary Table 1).
Association Study in Patients of European Ancestry
The subjects of the study included patients of both European and African origin and the allele frequency of the SNPs analyzed in this study differed between these groups (data not shown). Therefore, a sub-analysis was carried out, stratified for ethnicity. Of the 20 tagSNPs, five SNPs showed a trend or significant association with antipsychotic-induced weight gain in the European subgroup (Table 3 and Supplementary Table 1). To determine the contribution of these SNPs we performed backward linear regression. Duration of treatment was entered as an independent variable, and change in weight was the dependent variable. The polymorphism rs806378 emerged as the only significant SNP in genotypic comparison (B=2.14, β=0.234, p=0.035) as well as when the dominant effect of the minor allele T allele (B=3.33, β=0.286, p=0.010) was considered. In the genotypic model rs806378 explained 4.3% of the variance in the weight change, whereas it explained 7% in the dominant effect model. Carriers of the 'T' allele (CT+TT) gained more weight (5.96% ±5.54) than the CC genotype carriers (2.76% ±5.51, p=0.008, FDR q-value 0.12). This translated into patients with the T allele gaining approximately 2.2kg more than the patients homozygous for the CC genotype (CC vs CT+TT, 2.21±4.51 vs 4.33±3.89kg; p=0.022). This was also reflected in the allelic analysis (C vs T allele, 3.84 vs 5.83%, p=0.035; Figure 1, Supplementary Table 2). In literature, age, sex, duration of treatment and severity of schizophrenia (baseline BPRS score) have been associated with the propensity of the patients to gain weight. Using ANCOVA, we corrected for the confounding effect of these variables, the association of rs806378 with weight change (%) was significant (CC vs CT+TT; p=0.017).
Haplotypes were constructed using a 3-marker sliding-window across the gene region. Although none of the haplotypes reached significance at the global level, haplotypes including rs806378 or rs9450902 and rs10485170 showed a trend for association (Figure 1). LD was calculated for the European sub-sample and several blocks of LD were observed similar to previous studies (Figure 1). A high correlation was observed between rs806378 and rs806380 (r2=0.76); and rs9450902 and rs10485170 (r2=0.84). Thus, based on the cumulative evidence from genotypic and allelic comparisons, rs806378 appears to be the SNP-driving association with antipsychotic-induced weight gain in this region.
Association Study in Patients of African Ancestry
The African-American patients showed a higher percent weight gain than patients with European ancestry when the total sample was considered (8.01±7.57 (n=55) vs 4.17±5.88 (n=117); p=0.004) or when only patients on the higher risk medication, that is, clozapine and olanzapine were considered (9.33±7.99 (n=41) vs 4.41±5.66 (n=88), p=0.01). In this small sample a significant association was observed only for the synonymous polymorphism rs1049353 (CC vs CT; 8.4±7.59 (n=36) vs 16.05±8.29 (n=5), p=0.043).
Electrophoretic Mobility Shift Assay (EMSA)
The SNP rs806378 is present in either exon 3 (5′untranslated region) of the alternate transcript (CB1-E) or is in intron 1 of the transcript (CB1-A, Figure 1). In other genes, regions within the 5′untranslated region and in introns have been reported to affect gene expression (Makkonen et al, 2009; Ott et al, 2009).
We carried out an in silico prediction (MatInspector, Genomatix) of transcription factor binding using a 30-bp oligonucleotide surrounding SNP rs806378. The presence of the ancestral C-allele was predicted to bind X-binding protein-1 (XBP-1, ttatacgACGTgatgag, the SNP is in bold and core binding site in capital letters), a member of the cAMP-response element-binding protein (Core similarity: 1.00, matrix similarity: 0.88). The T-allele was predicted to bind to arylhydrocarbon receptor translocator (ARNT, ttatacaaCGTGatgag), a member of the basic helix—loop—helix/Per—Arnt—Sim protein family (bHLH/PAS; core similarity: 1.00; matrix similarity: 0.897) but not to XBP-1.
An electrophoretic mobility shift assay was carried out and putative binding of ARNT to the T-allele was observed (Figure 2). Furthermore, in the competition assay, unlabelled oligonucleotide containing the consensus ARNT binding site was able to compete out the biotin-labeled T-allele carrying oligonucleotide, suggesting that the nuclear protein binding to this sequence is ARNT. The ancestral C-allele did not show binding to either ARNT or XBP-1 or any other nuclear protein.
Source, Graphs and Figures: A Common Polymorphism in the Cannabinoid Receptor 1 (CNR1) Gene is Associated with Antipsychotic-Induced Weight Gain in Schizophrenia
Antipsychotic-induced weight gain has emerged as a serious complication in the treatment of patients with atypical antipsychotic drugs. The cannabinoid receptor 1 (CNR1) is expressed centrally in the hypothalamic region and associated with appetite and satiety, as well as peripherally. An antagonist of CNR1 (rimonabant) has been effective in causing weight loss in obese patients indicating that CNR1 might be important in antipsychotic-induced weight gain. Twenty tag SNPs were analyzed in 183 patients who underwent treatment (with either clozapine, olanzapine, haloperidol, or risperidone) for chronic schizophrenia were evaluated for antipsychotic-induced weight gain for up to 14 weeks. The polymorphism rs806378 was nominally associated with weight gain in patients of European ancestry treated with clozapine or olanzapine. 'T' allele carriers (CT+TT) gained more weight (5.96%), than the CC carriers (2.76%, p=0.008, FDR q-value=0.12). This translated into approximately 2.2kg more weight gain in patients carrying the T allele than the patients homozygous for the CC genotype (CC vs CT+TT, 2.21±4.51 vs 4.33±3.89kg; p=0.022). This was reflected in the allelic analysis (C vs T allele, 3.84 vs 5.83%, p=0.035). We conducted electrophoretic mobility shift assays which showed that the presence of the T allele created a binding site for arylhydrocarbon receptor translocator (ARNT), a member of the basic helix—loop—helix/Per—Arnt—Sim protein family. In this study, we provide evidence that the CNR1 gene may be associated with antipsychotic-induced weight gain in chronic schizophrenia patients. However, these observations were made in a relatively small patient population; therefore these results need to be replicated in larger sample sets.
INTRODUCTION
The use of atypical antipsychotic drugs (AADs) for treatment of schizophrenia has increased and they now capture over 90% of the antipsychotic drug market. Unlike the typical antipsychotic drugs, AADs have the primary advantage of being associated with very low or negligible risk of inducing extrapyramidal symptoms. However, one of the major concerns with AADs is the development of weight gain leading to full metabolic syndrome and early cardiac death (reviewed by Muller and Kennedy, 2006). Obesity in the general population has become a major public health challenge and association of AADs with weight gain adds another dimension to this problem.
The development of antipsychotic induced weight gain in a subset of schizophrenia patients (∼30%), concordance among twins, and high heritability of weight gain suggest a strong genetic component. Two polymorphisms, −759C/T in the serotonin receptor 5-HT2C, and the −2548 G>A promoter polymorphism in the leptin gene, are most consistently associated with antipsychotic-induced weight gain (eg, De Luca et al, 2007; Ellingrod et al, 2007; Templeman et al, 2005; Zhang et al, 2003). Other gene findings have been reviewed elsewhere (Muller and Kennedy, 2006). However, it is more than plausible to assume that, as for other complex inherited phenotypes, numerous and yet undiscovered genetic markers are involved in antipsychotic-induced weight gain.
The cannabinoid receptor 1 (CNR1 or CB1) is a seven transmembrane G-coupled receptor (Gi/o, Gs or Gq/11), expressed widely in several brain regions associated with reward, hunger, appetite and satiety and in peripheral tissues including the gut, liver (hepatocytes), adipocytes, skeletal muscle and pancreas (Ameri et al, 1999; Matias and Di Marzo, 2007; Woods, 2007). CNR1 knockout mice present a lean phenotype (Cota et al, 2003; Di Marzo et al, 2001), reduced food intake and resistance to diet-induced obesity (Osei-Hyiaman et al, 2005; Ravinet Trillou et al, 2003, 2004). Similarly, administration of rimonabant, a CNR1 receptor antagonist, leads to reduced food intake, decreased bodyweight, normalized plasma leptin, insulin, free fatty acids and corrected insulin resistance in wild-type animals (Jamshidi and Taylor, 2001; Ravinet Trillou et al, 2003). Conversely, injection of endogenous ligands of CNR1, such as anandamide or 2-arachidonoylglycerol, into the hypothalamus of pre-satiated rats increases dietary food intake. Furthermore, this increase in food intake could be prevented by pretreatment of the rats with rimonabant (Jamshidi and Taylor, 2001; Kirkham et al, 2002). Weston-Green et al (2008) showed that chronic administration of olanzapine decreases the binding of the agonist, [3H] CP-55940, to CNR1 in the dorsal vagal complex of the brain stem and causes increased weight gain in rats, whereas haloperidol failed to elicit similar effects. In humans, treatment of AIDS and cancer-related anorexia with dronabinol, a CNR1 agonist, leads to significant increase in appetite and body weight compared with placebo (Beal et al, 1995; Nelson et al, 1994). Furthermore, increased CNR1 mRNA levels in visceral adipose tissue of overweight and metabolic syndrome patients have been reported (Sarzani et al, 2009). In addition, hepatic stimulation of the CNR1 receptor leads to increased lipogenesis through augmented expression of the lipogenic transcription factor SREBP-1c (sterol regulatory element-binding protein-1c), an effect inhibited by treatment with rimonabant (Osei-Hyiaman et al, 2005).
At the molecular level, CNR1 interacts with several molecules in the hypothalamus. For example CNR1 downregulates anorexigenic molecules such as cocaine and amphetamine-regulated transcript in the arcuate nucleus and corticotrophin-releasing hormone (CRH) in the paraventricular nucleus. CNR1 also upregulates orexigenic neuropeptide Y and melanin-concentrating hormone (MCH) in the lateral hypothalamus (Cota et al, 2003; Gamber et al, 2005; Osei-Hyiaman et al, 2005). Some of the effects of CNR1 are under regulation of leptin, which prevents the release of orexigenic and anorexigenic molecules (MCH and CRH, respectively; for detailed review see (Kunos, 2007; Matias and Di Marzo, 2007; Woods, 2007). Similarly, ghrelin is known to exert its orexigenic effects at least in part through CNR1 (Kola et al, 2008).
The CNR1 gene is located on 6q14—15, and a common polymorphism in intron 2, rs20232239, has been associated with obesity in the general population (Benzinou et al, 2008) as well as with substance dependence and CNR1 mRNA levels (Zhang et al, 2004). Furthermore, the C allele of this polymorphism has been associated with greater CB1 binding in the prefrontal cortex and alcohol cue-elicited brain activation (Hutchison et al, 2008).
Based on the above studies, the CNR1 gene represents a compelling candidate for genetic studies of antipsychotic-induced weight gain. In this study, we analyzed 20 tag polymorphisms in the CNR1 gene, including the synonymous polymorphism rs1049353 (1359 G>A, Thr453Thr) and obesity-associated rs20232239 for involvement in antipsychotic-induced weight gain.
MATERIALS AND METHODS
Subjects
A total of 183 patients with chronic schizophrenia or schizoaffective disorder were included in this study. Approval from the institutional ethics committee and informed consent were obtained for all patients. The details of the sample are provided in Table 1. Patients were recruited at three different sites: Sample A: Charite University Medicine, Berlin, Germany (DJ Müller and I Puls, n=47). Patients 18—60 years old diagnosed with schizophrenia or schizoaffective disorder according to DSM-IV criteria were included. Other inclusion criteria were a minimum positive and negative symptom scale (PANSS) score of 60 and admission either due to first manifestation or relapse of pre-existent condition with significant deterioration. Patients were given antipsychotic medication and assessed up to 6 weeks (Table 1). Exclusion criteria included pregnancy or breast feeding, organic brain disorder or severe head injuries, severe or unstable comorbid medical conditions (Hepatitis C, HIV, Thyroid disorder, diabetes mellitus), and not stable, substance dependence (polytoxicomania), clinically relevant mental retardation, and severe personality disorder.
Sample B: Case Western Reserve University in Cleveland, Ohio (HY Meltzer, JA Lieberman, n=79). Patients 18—60 years old diagnosed with schizophrenia according to DSMIIR criteria who were either treatment refractory or intolerant to treatment with typical antipsychotics were included. Before starting treatment with clozapine there was a washout period of 7—14 days during which patients received no medications unless it was clinically essential. Clozapine dosage was not fixed and serum levels were monitored to ascertain compliance. Patients were treated with clozapine for 6 weeks and percent change in weight was calculated from baseline weight (for details see (Masellis et al, 1998)). Exclusion criteria were pregnancy, organic brain disorder or severe head injuries, previous medical conditions (Hepatitis C, HIV, Thyroid disorder, diabetes mellitus), which require treatment and are not stable, substance dependence (polytoxicomania), clinically relevant mental retardation and severe personality disorder.
Sample C: Hillside Hospital in Glen Oaks, New York (JA Lieberman, n=57). Patients 18—60 years old, with chronic schizophrenia or schizoaffective disorder according to DSM-IV were included. These patients had shown suboptimal response to previous treatment with typical antipsychotic drugs, mainly defined by persistent positive symptoms and a poor level of functioning over the past 2 years. Patients were randomly assigned to receive either clozapine (500mg/day), olanzapine (20mg/day), risperidone (8mg/day) or haloperidol in a 14-week double-blind study. Patients with a history of non-response to clozapine, risperidone, or olanzapine, defined as an unambiguous lack of improvement despite a contiguous adequate trial of risperidone or olanzapine for at least 6 weeks, or clozapine for at least 14 weeks were excluded. Similarly, patients with a history of clozapine, olanzapine, risperidone, or haloperidol intolerance as well as those who received a depot antipsychotic within 30 days before randomization were also excluded (for details see (Volavka et al, 2002)). Sample B and C have been described in more detail in previous studies (Basile et al, 2002; De Luca et al, 2007; Muller et al, 2005a,2005b).
Ethnicity was determined using ancestry information provided by the participants. Patients were classified for ethnicity by assessing ancestry over last three generations.
Genetic Analysis
10ml of venous blood was obtained from study participants and DNA was extracted at the Centre for Addiction and Mental Health using the high salt method (Lahiri and Nurnberger, 1991). All the genotyping was carried out using TaqMan SNP genotyping assays (Applied Biosystems, Foster City, CA). Genotype calling was confirmed by two independent researchers and 5% of the total sample was re-genotyped to assure concordance. The concordance rate for this study was 99.86% and discordant genotypes were treated as missing in the statistical analysis. A total of 20 tagSNPs covering ∼32,500bp and representing 100% of the common variation in the CNR1 gene region (minor allele frequency >0.05) were included in the study (Figure 1). Among the 20 SNPs rs806381, rs2023239, rs9450902 and rs10485170 were force-included as they were recently reported to be associated with weight gain in the general population (Benzinou et al, 2008).
Data Analysis
All the statistical analyses were performed using SPSS 15.0. The categorical variables were tested using a χ2-test and the continuous variables were tested using ANOVA, ANCOVA, or linear regression. Pearson's correlation was calculated between age, sex, duration of treatment, baseline Brief Psychiatric Rating Scale (BPRS). We decided a priori that the comparison of mean weight change (%) across genotypic categories would be the association-determining test. Only SNPs found to show a trend or significance were used for post hoc analyses. Haplotypes were constructed and percent change in weight was compared using QTPHASE routine in UNPHASED 2.0 (Dudbridge, 2003). Only haplotypes with frequencies greater than 5% were included in the analyses. Linkage disequilibrium (LD) and tagSNPs (r2=0.8, minor allele frequency >0.05) were determined using Haploview 4.1 (Barrett et al, 2005). LD blocks were constructed using solid spine of LD. Power calculations were carried out using Quanto 1.2.4 (Gauderman and Morrison, 2006). Assuming a minor allele frequency of 0.25 in a sample size of n=183, and a mean difference of 2.5% between carriers of the risk genotype, we had more than 80% power in a dominant and additive effect model. False discovery rate (FDR) q-values were calculated using the program QVALITY. Similar to other publications, a q-value threshold 0.20 was used to indicate SNPs with most significant associations (van den Berg et al, 2009).
Electrophoretic Mobility Shift Assay
Biotin-labelled and unlabelled double-stranded oligonucleotides for each allele of rs806378 (T and C) and unlabelled oligonucleotides with the specific consensus-binding site for arylhydrocarbon receptor translocator (ARNT) and XBP1 (Clauss et al, 1996) were synthesized (Table 2, Eurofin MWG Operon, Huntsville, AL). Nuclear extract was prepared from SH-SY5Y neuroblastoma cells (ATCC: CRL2266) using the NE-PER Nuclear and cytoplasmic extraction kit (Pierce, IL). The complementary pairs of nucleotides suspended in STE Buffer (10mM Tris pH 8.0, 50mM NaCl, 1mM EDTA) were annealed using a thermocycler (Step 1: 95°C for 5min, Step 2: 95°C (−1°C/cycle, 70 times), Step 3: HOLD at 4°C). The reannealed double-stranded oligonucleotides (20fmol) were incubated with 2μl of the nuclear extract, 1 × binding buffer, 2.5% glycerol, 5mM MgCl2, 50ng/μl poly (dI dC), 0.05% NP-40 for 20min at room temperature. For the competition experiments, unlabelled oligonucleotides at 200-fold molar excess (4pmol) were preincubated with the nuclear extract before adding the biotin-labelled oligobucleotides. DNA protein complexes were fractioned on a 5% polyacrylamide gel in 0.5 × TBE for 35min at 4°C. Complexes were electro-transferred to nylon membrane at 380mA for 50min, and cross-linked at 120mJ/cm2 for 45s. The biotinylated oligonucleotides were detected by chemiluminescence using the Chemiluminescent nucleic acid detection module of the Light Shift Chemiluminescent EMSA kit (Pierce, IL). The processed nylon membrane was exposed to an X-ray film for 2min.
RESULTS
Genetic Analysis
All the SNPs were in Hardy—Weinberg equilibrium (p>0.05, Table 3, Supplementary Table 1). Several clinical variables were different across the three sites, however, similar weight gain was observed across the clinical sites (Table 1). In the total sample set significant correlation was observed between weight change (%) and duration of treatment (p=0.045) but not with age (p=0.747) and baseline BPRS score (p=0.986). In the subgroup of patients of European ancestry on clozapine or olanzapine no significant correlation of duration of treatment (p=0.532), age (p=0.703), and baseline BPRS score (p=0.878) with weight change (%) was observed. Similarly, weight change (%) was not different between the male and female patients in the total (p=0.776) as well as the sub-group of patients of European ancestry on clozapine or olanzapine (p=0.822).
Association Study in the Total Sample
No association of any of the polymorphisms with weight change (%) was observed when we considered the total sample set of patients on various antipsychotics (p>0.05). Consistent with the literature (eg, reviewed by (Nasrallah, 2008)), we observed highest weight gain rates with clozapine (5.37%±6.27) and olanzapine (7.29%±8.17) in our sample compared with haloperidol (2.76%±4.29) and risperidone (3.54%±6.32). Clozapine and olanzapine have relatively similar receptor-binding profile (Nasrallah, 2008). Therefore, we assessed the subgroup of patients who received either clozapine or olanzapine for association with polymorphisms in CNR1 gene. No association of any of the SNPs with weight change was observed (p>0.05, Table 3, Supplementary Table 1).
Association Study in Patients of European Ancestry
The subjects of the study included patients of both European and African origin and the allele frequency of the SNPs analyzed in this study differed between these groups (data not shown). Therefore, a sub-analysis was carried out, stratified for ethnicity. Of the 20 tagSNPs, five SNPs showed a trend or significant association with antipsychotic-induced weight gain in the European subgroup (Table 3 and Supplementary Table 1). To determine the contribution of these SNPs we performed backward linear regression. Duration of treatment was entered as an independent variable, and change in weight was the dependent variable. The polymorphism rs806378 emerged as the only significant SNP in genotypic comparison (B=2.14, β=0.234, p=0.035) as well as when the dominant effect of the minor allele T allele (B=3.33, β=0.286, p=0.010) was considered. In the genotypic model rs806378 explained 4.3% of the variance in the weight change, whereas it explained 7% in the dominant effect model. Carriers of the 'T' allele (CT+TT) gained more weight (5.96% ±5.54) than the CC genotype carriers (2.76% ±5.51, p=0.008, FDR q-value 0.12). This translated into patients with the T allele gaining approximately 2.2kg more than the patients homozygous for the CC genotype (CC vs CT+TT, 2.21±4.51 vs 4.33±3.89kg; p=0.022). This was also reflected in the allelic analysis (C vs T allele, 3.84 vs 5.83%, p=0.035; Figure 1, Supplementary Table 2). In literature, age, sex, duration of treatment and severity of schizophrenia (baseline BPRS score) have been associated with the propensity of the patients to gain weight. Using ANCOVA, we corrected for the confounding effect of these variables, the association of rs806378 with weight change (%) was significant (CC vs CT+TT; p=0.017).
Haplotypes were constructed using a 3-marker sliding-window across the gene region. Although none of the haplotypes reached significance at the global level, haplotypes including rs806378 or rs9450902 and rs10485170 showed a trend for association (Figure 1). LD was calculated for the European sub-sample and several blocks of LD were observed similar to previous studies (Figure 1). A high correlation was observed between rs806378 and rs806380 (r2=0.76); and rs9450902 and rs10485170 (r2=0.84). Thus, based on the cumulative evidence from genotypic and allelic comparisons, rs806378 appears to be the SNP-driving association with antipsychotic-induced weight gain in this region.
Association Study in Patients of African Ancestry
The African-American patients showed a higher percent weight gain than patients with European ancestry when the total sample was considered (8.01±7.57 (n=55) vs 4.17±5.88 (n=117); p=0.004) or when only patients on the higher risk medication, that is, clozapine and olanzapine were considered (9.33±7.99 (n=41) vs 4.41±5.66 (n=88), p=0.01). In this small sample a significant association was observed only for the synonymous polymorphism rs1049353 (CC vs CT; 8.4±7.59 (n=36) vs 16.05±8.29 (n=5), p=0.043).
Electrophoretic Mobility Shift Assay (EMSA)
The SNP rs806378 is present in either exon 3 (5′untranslated region) of the alternate transcript (CB1-E) or is in intron 1 of the transcript (CB1-A, Figure 1). In other genes, regions within the 5′untranslated region and in introns have been reported to affect gene expression (Makkonen et al, 2009; Ott et al, 2009).
We carried out an in silico prediction (MatInspector, Genomatix) of transcription factor binding using a 30-bp oligonucleotide surrounding SNP rs806378. The presence of the ancestral C-allele was predicted to bind X-binding protein-1 (XBP-1, ttatacgACGTgatgag, the SNP is in bold and core binding site in capital letters), a member of the cAMP-response element-binding protein (Core similarity: 1.00, matrix similarity: 0.88). The T-allele was predicted to bind to arylhydrocarbon receptor translocator (ARNT, ttatacaaCGTGatgag), a member of the basic helix—loop—helix/Per—Arnt—Sim protein family (bHLH/PAS; core similarity: 1.00; matrix similarity: 0.897) but not to XBP-1.
An electrophoretic mobility shift assay was carried out and putative binding of ARNT to the T-allele was observed (Figure 2). Furthermore, in the competition assay, unlabelled oligonucleotide containing the consensus ARNT binding site was able to compete out the biotin-labeled T-allele carrying oligonucleotide, suggesting that the nuclear protein binding to this sequence is ARNT. The ancestral C-allele did not show binding to either ARNT or XBP-1 or any other nuclear protein.
Source, Graphs and Figures: A Common Polymorphism in the Cannabinoid Receptor 1 (CNR1) Gene is Associated with Antipsychotic-Induced Weight Gain in Schizophrenia
