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.
DISCUSSION
The cannabinoid receptor 1 as part of the endocannabinoid system has emerged as an important receptor modulating diverse effects through its expression in multiple parts of the brain and peripheral tissues. CNR1 is considered to be important for modulating weight gain as evidenced by the effect of CNR1-antagonist rimonabant, which has been shown to foster weight loss both in animal studies as well as in clinical trials on humans. Genetic polymorphisms in the CNR1 gene have been associated with basal metabolic index, obesity and various metabolic parameters (Aberle et al, 2007; Baye et al, 2008; Benzinou et al, 2008; Gazzerro et al, 2007; Jaeger et al, 2008; Peeters et al, 2007; Russo et al, 2007), however some negative studies have also been reported (Aberle et al, 2007; Lieb et al, 2009; Muller et al, 2007). Significant association of polymorphisms in CNR1 with anorexia nervosa and bulimia nervosa have also been reported (Monteleone et al, 2009; Siegfried et al, 2004), despite one negative report (Muller et al, 2008). Based on these many significant findings, we hypothesized that CNR1 might be associated with antipsychotic-induced weight gain.
In this study, we have tested the association of SNPs in the CNR1 gene with antipsychotic-induced weight gain in a group of well-characterized schizophrenia patients. The CNR1 gene was extensively covered with markers; however, of the 20 tagSNPs investigated, we did not observe any significant association in the total sample set. Reasons for this negative finding in the whole sample could be due to the inclusion of patients who received drugs, which have lesser propensity to cause weight gain (such as haloperidol and risperidone) and different ancestries. Similarly, the lack of association observed in the patients on higher risk medication for weight gain (clozapine and olanzapine) could be potentially due to inclusion of patients of different ancestries (European and African ancestry), which may have different susceptibility loci. We therefore stratified the samples according to their self-reported ethnicities and medication and observed significant association in patients of both African and European ancestry. In the patients with African ancestry, we observed a significant association with the synonymous polymorphism rs1049353 (C>T). However, the sample set is small and the observation became non-significant after correction for multiple comparisons (Bonferroni corrected p [0.043 × 20]=0.86).
We observed a significant effect of the polymorphism rs806378 (C>T) with the propensity to gain weight in the larger sample of patients of European ancestry who had received either clozapine or olanzapine. This polymorphism is present in the 5′UTR of the CB1-E transcript or intron2 of CB1-A (Figure 1). Our analyses indicate that carriers of the T-allele have a higher risk for weight gain. Although no functional effect has yet been described, the derived 'T' allele of this polymorphism leads to generation of a transcription factor-binding site for ARNT, a member of the bHLH/PAS protein family. Using EMSA we provided evidence that ARNT is indeed the nuclear protein binding to this sequence and that the ancestral C-allele does not bind to either ARNT or any other nuclear protein. ARNT is known to heterodimerize with other bHLH-PAS proteins, primarily with single minded-1 (SIM1) and arylhydrocarbon receptor AHR. The AHR—ARNT heterodimer preferentially binds to TNGCGTG, the dioxin responsive element, whereas the SIM1—ARNT heterodimer binds to GT(G/A)CGTG (Yang et al, 2004). The latter appears to be more similar to the site containing the rs806378 polymorphism (CAACGTG). Interestingly SIM1 is essential for differentiation of neurons in the paraventricular nucleus, a region of hypothalamus critical for regulating physiological processes such as energy homeostasis and blood pressure. Furthermore, haploinsufficiency and disruption of SIM1 are associated with hyperphagic obesity in humans (Faivre et al, 2002; Holder et al, 2000). Similarly, mice with heterozygous deletion of Sim1 develop hyperphagia and obesity (Michaud et al, 2001). Thus the presence of this transcription factor-binding site may lead to a novel regulatory mechanism involving ARNT—SIM1-mediated transcriptional alteration.
The limitations of the significant findings in our stratified sub-samples include the relatively limited sample size, inclusion of self-reported ancestry, and lack of information on drug dose, cannabis use and diet. Our observation of different percent weight gain among CC vs CT+TT genotype (p=0.008) of rs806378 in patients of European ancestry on high risk medication did not withstand conservative corrections for multiple testing (Bonferroni corrected p=0.15, Table 3). However, the FDR q-value of 0.12 suggests that this observation is at least nominally significant. Furthermore, we had a strong biological rationale to test this gene for association with antipsychotic-induced weight gain (Benzinou et al, 2008; Hutchison et al, 2008; Weston-Green et al, 2008). The recent report of an association between the same T-allele of rs806378 with general obesity in a sample of French individuals further lends support to this observation (Benzinou et al, 2008). As mentioned earlier, SREBP expression is activated by CNR1 and recently Le Hellard et al (2009) reported a significant association of insulin-induced gene 2 (INSIG2) in a German sample. INSIG2 is expressed in the endoplasmic reticulum and sequesters SREBP by forming a trimeric complex along with SREBP cleavage-activating protein. This provides further evidence that genes involved in this pathway might have an important role in antipsychotic-induced weight gain.
In conclusion, the observations made in this study indicate a possible role of polymorphisms in CNR1 gene in antipsychotic-induced weight gain. However, to completely elucidate the role of CNR1these observations need to be further investigated and replicated in larger independent sample sets. In addition, it may be interesting to examine CNR1 for association with weight gain in other psychiatric disorders treated with atypical antipsychotic drugs.
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.
DISCUSSION
The cannabinoid receptor 1 as part of the endocannabinoid system has emerged as an important receptor modulating diverse effects through its expression in multiple parts of the brain and peripheral tissues. CNR1 is considered to be important for modulating weight gain as evidenced by the effect of CNR1-antagonist rimonabant, which has been shown to foster weight loss both in animal studies as well as in clinical trials on humans. Genetic polymorphisms in the CNR1 gene have been associated with basal metabolic index, obesity and various metabolic parameters (Aberle et al, 2007; Baye et al, 2008; Benzinou et al, 2008; Gazzerro et al, 2007; Jaeger et al, 2008; Peeters et al, 2007; Russo et al, 2007), however some negative studies have also been reported (Aberle et al, 2007; Lieb et al, 2009; Muller et al, 2007). Significant association of polymorphisms in CNR1 with anorexia nervosa and bulimia nervosa have also been reported (Monteleone et al, 2009; Siegfried et al, 2004), despite one negative report (Muller et al, 2008). Based on these many significant findings, we hypothesized that CNR1 might be associated with antipsychotic-induced weight gain.
In this study, we have tested the association of SNPs in the CNR1 gene with antipsychotic-induced weight gain in a group of well-characterized schizophrenia patients. The CNR1 gene was extensively covered with markers; however, of the 20 tagSNPs investigated, we did not observe any significant association in the total sample set. Reasons for this negative finding in the whole sample could be due to the inclusion of patients who received drugs, which have lesser propensity to cause weight gain (such as haloperidol and risperidone) and different ancestries. Similarly, the lack of association observed in the patients on higher risk medication for weight gain (clozapine and olanzapine) could be potentially due to inclusion of patients of different ancestries (European and African ancestry), which may have different susceptibility loci. We therefore stratified the samples according to their self-reported ethnicities and medication and observed significant association in patients of both African and European ancestry. In the patients with African ancestry, we observed a significant association with the synonymous polymorphism rs1049353 (C>T). However, the sample set is small and the observation became non-significant after correction for multiple comparisons (Bonferroni corrected p [0.043 × 20]=0.86).
We observed a significant effect of the polymorphism rs806378 (C>T) with the propensity to gain weight in the larger sample of patients of European ancestry who had received either clozapine or olanzapine. This polymorphism is present in the 5′UTR of the CB1-E transcript or intron2 of CB1-A (Figure 1). Our analyses indicate that carriers of the T-allele have a higher risk for weight gain. Although no functional effect has yet been described, the derived 'T' allele of this polymorphism leads to generation of a transcription factor-binding site for ARNT, a member of the bHLH/PAS protein family. Using EMSA we provided evidence that ARNT is indeed the nuclear protein binding to this sequence and that the ancestral C-allele does not bind to either ARNT or any other nuclear protein. ARNT is known to heterodimerize with other bHLH-PAS proteins, primarily with single minded-1 (SIM1) and arylhydrocarbon receptor AHR. The AHR—ARNT heterodimer preferentially binds to TNGCGTG, the dioxin responsive element, whereas the SIM1—ARNT heterodimer binds to GT(G/A)CGTG (Yang et al, 2004). The latter appears to be more similar to the site containing the rs806378 polymorphism (CAACGTG). Interestingly SIM1 is essential for differentiation of neurons in the paraventricular nucleus, a region of hypothalamus critical for regulating physiological processes such as energy homeostasis and blood pressure. Furthermore, haploinsufficiency and disruption of SIM1 are associated with hyperphagic obesity in humans (Faivre et al, 2002; Holder et al, 2000). Similarly, mice with heterozygous deletion of Sim1 develop hyperphagia and obesity (Michaud et al, 2001). Thus the presence of this transcription factor-binding site may lead to a novel regulatory mechanism involving ARNT—SIM1-mediated transcriptional alteration.
The limitations of the significant findings in our stratified sub-samples include the relatively limited sample size, inclusion of self-reported ancestry, and lack of information on drug dose, cannabis use and diet. Our observation of different percent weight gain among CC vs CT+TT genotype (p=0.008) of rs806378 in patients of European ancestry on high risk medication did not withstand conservative corrections for multiple testing (Bonferroni corrected p=0.15, Table 3). However, the FDR q-value of 0.12 suggests that this observation is at least nominally significant. Furthermore, we had a strong biological rationale to test this gene for association with antipsychotic-induced weight gain (Benzinou et al, 2008; Hutchison et al, 2008; Weston-Green et al, 2008). The recent report of an association between the same T-allele of rs806378 with general obesity in a sample of French individuals further lends support to this observation (Benzinou et al, 2008). As mentioned earlier, SREBP expression is activated by CNR1 and recently Le Hellard et al (2009) reported a significant association of insulin-induced gene 2 (INSIG2) in a German sample. INSIG2 is expressed in the endoplasmic reticulum and sequesters SREBP by forming a trimeric complex along with SREBP cleavage-activating protein. This provides further evidence that genes involved in this pathway might have an important role in antipsychotic-induced weight gain.
In conclusion, the observations made in this study indicate a possible role of polymorphisms in CNR1 gene in antipsychotic-induced weight gain. However, to completely elucidate the role of CNR1these observations need to be further investigated and replicated in larger independent sample sets. In addition, it may be interesting to examine CNR1 for association with weight gain in other psychiatric disorders treated with atypical antipsychotic drugs.
Source, Graphs and Figures: A Common Polymorphism in the Cannabinoid Receptor 1 (CNR1) Gene is Associated with Antipsychotic-Induced Weight Gain in Schizophrenia
