Skip to main content

ORIGINAL RESEARCH article

Front. Genet., 14 June 2018
Sec. Behavioral and Psychiatric Genetics

Detection of Significant Association Between Variants in Cannabinoid Receptor 1 Gene (CNR1) and Personality in African–American Population

\r\nYinghao YaoYinghao Yao1Yi XuYi Xu1Junsheng ZhaoJunsheng Zhao1Yunlong MaYunlong Ma1Kunkai SuKunkai Su1Wenji YuanWenji Yuan1Jennie Z. MaJennie Z. Ma2Thomas J. PayneThomas J. Payne3Ming D. Li,,*Ming D. Li1,4,5*
  • 1State Key Laboratory for Diagnosis and Treatment of Infectious Diseases, The First Affiliated Hospital, Collaborative Innovation Center for Diagnosis and Treatment of Infectious Diseases, Zhejiang University School of Medicine, Hangzhou, China
  • 2Department of Public Health Sciences, University of Virginia, Charlottesville, VA, United States
  • 3Department of Otolaryngology and Communicative Sciences, University of Mississippi Medical Center, Jackson, MS, United States
  • 4Research Center for Air Pollution and Health, Zhejiang University, Hangzhou, China
  • 5Institute of Neuroimmune Pharmacology, Seton Hall University, South Orange, NJ, United States

Background: Several studies have revealed significant associations between single nucleotide polymorphisms (SNPs) in the cannabinoid receptor 1 (CNR1) gene and a broad spectrum of psychiatric disorders such as major depressive disorder (MDD), attention deficit hyperactivity disorder (ADHD), and schizophrenia. Personality traits that are highly related to susceptibility to these conditions have been associated with the CNR1 variants in subjects of Caucasian origin. However, there are no reported studies regarding the effects of CNR1 polymorphisms on personality traits in the African-American (AA) population.

Methods: We performed an imputation-based association analysis for 26 CNR1 variants with five dimensions of personality in 3,046 AAs.

Results: SNPs rs806372 and rs2180619 showed a significant association with extraversion after Bonferroni correction for multiple testing (p < 0.0019). Further, several extraversion-associated SNPs were significantly associated with conscientiousness, agreeableness, and openness. SNP priority score analysis indicated that SNPs rs806368, rs806371, and rs2180619 play a role in the modulation of personality and psychiatric conditions.

Conclusion: CNR1 is important in determining personality traits in the AA population.

Introduction

Personality traits can affect human behavior and predict health outcomes. They commonly are measured by the revised NEO Personality Inventory (NEO-PI-R) (Costa et al., 1992). The NEO-PI-R examines the following five factors: agreeableness, conscientiousness, extraversion, neuroticism, and openness (McCrae et al., 2005). These five dimensions, which account for nearly all of the differences in personality between individuals, have been linked to emotional stability, active motivation, and cognition (Terracciano et al., 2008b; DeYoung et al., 2010). Furthermore, the score for each personality trait has been used to predict many psychiatric disorders. For example, high scores of neuroticism, extraversion, and openness are considered indicators of bipolar disorder (Barnett et al., 2011), a high neuroticism score is associated with major depression disorder (MDD) and anxiety (Hettema et al., 2006), and a low conscientiousness score renders one prone to MDD (Kendler and Myers, 2010). In addition, substance abuse behaviors are related to personality. For example, tobacco smokers and cocaine and heroin users score high on neuroticism and low on conscientiousness (Terracciano and Costa, 2004; Terracciano et al., 2008a; Choi et al., 2017).

Family, twin, and adoption studies have demonstrated that personality traits are highly heritable, with an estimated heritability ranging from 40 to 60% (Horn et al., 1976; Floderus-Myrhed et al., 1980; Jang et al., 1996). Recent advances in genome-wide association studies (GWAS) have provided a genetic map for various personality traits (Amin et al., 2012, 2013; de Moor et al., 2012; Kim et al., 2013, 2015; Genetics of Personality Consortium et al., 2015; Okbay et al., 2016; Smith et al., 2016; van den Berg et al., 2016; Lo et al., 2017; Luciano et al., 2018). Of them, neuroticism and extraversion are the two best-researched traits, uncovering a number of related loci in independent samples (Calboli et al., 2010; Genetics of Personality Consortium et al., 2015; Okbay et al., 2016; Smith et al., 2016; van den Berg et al., 2016; Luciano et al., 2018). However, although the most significant loci for several personality traits have been identified, only a small fraction of heritability could be explained by these top hits, indicating that single nucleotide polymorphisms (SNPs) with effect sizes well below genome-wide significance likely account for the “missing heritability” (Manolio et al., 2009; Yang et al., 2010). Thus, the genetic architecture of personality requires further investigation, especially for those loci below the commonly accepted GWAS threshold (Boyle et al., 2017).

Although it has not been revealed to be significant at the genome-wide significant level by any GWAS, cannabinoid receptor 1 (CNR1) is a plausible candidate gene for certain personality traits. Juhasz et al. (2009a) first reported a significant association of CNR1 with personality traits in a Caucasian population. By genotyping seven tag SNPs plus three predictive expression quantitative loci, they detected several haplotypes that were significantly associated with neuroticism and agreeableness, and substantial phenotypic variances were explained by the CNR1 gene variants. Furthermore, a recent study showed an association between CNR1 rs7766029 and neuroticism (Aleksandrova et al., 2012). In addition, a genome-wide linkage study based on the Erasmus Rucphen Family sample (Pardo et al., 2005) detected a significant linkage signal for extraversion at a genomic region near CNR1 (Amin et al., 2013). The protein encoded by this gene belongs to the G-protein receptor family and modulates neurotransmitter release by coupling with a decreased intracellular cAMP concentration (Elphick and Egertova, 2001). CB1 receptors are densely expressed in the central nervous system and act as neuromodulators to inhibit the release of glutamate and GABA (Elphick and Egertova, 2001).

Several genetic association studies have implicated the CNR1 in the risk of several psychiatric disorders (Juhasz et al., 2009a,b, 2017; Lazary et al., 2009). It has been reported that genetic variants in the promoter region of CNR1 and serotonin transporter gene (SLC6A4) interactively increase the risk of high anxiety scores (Lazary et al., 2009), and interaction of CNR1 variants with recent negative life events is considered to be an important risk factor for development of depression symptoms (Juhasz et al., 2009a) and migraine (Juhasz et al., 2017). Furthermore, additional data indicate a nominal association between CNR1 variants and metabolic syndrome in patients with schizophrenia (Ujike et al., 2002; Martinez-Gras et al., 2006; Yu et al., 2013) and attention deficit hyperactivity disorder (ADHD) (Lu et al., 2008).

Given that the CNR1 gene has been implicated in the pathology of a wide range of psychiatric disorders, examining its effects on personality traits is of great interest, as they are regarded as risk factors for these disorders. Although a significant association of the gene with personality has been uncovered in Caucasian populations, little is known about non-Caucasian populations such as persons of African origin. Thus, in the present study, we focused on determining whether there exists any significant association between the CNR1 variants and five dimensions of personality in the African-American (AA) population.

Materials and Methods

Subjects

The 3,046 AA subjects used in this study were recruited during 2005–2011 as part of the Mid-South Tobacco Case-Control (MSTCC) study (Yang et al., 2015). Of note, all AA samples used in this study were recruited from the Jackson area of Mississippi, where residents do not relocate or migrate as often as residents of other large cities within the United States, implying that they are relatively more homogeneous. These subjects aged 18 or older were all biologically unrelated. The NEO Five-Factor Inventory (NEO-FFI) questionnaire (McCrae et al., 2005) was used to assess the five dimensions of personality. The characteristics of the MSTCC AA sample are shown in Table 1. After a detailed explanation of the project was provided to potential participants, informed written consent was obtained. The Institutional Review Boards of the University of Virginia and University of Mississippi Medical Center approved this study.

TABLE 1
www.frontiersin.org

TABLE 1. Characteristics of the five dimensions of personality traits.

Questionnaires

The background information collected from participants included age, sex, ethnicity, education, and smoking history. Sixty items of NEO-FFI (Costa and McCrae, 1997) were used to define the five personality dimensions of agreeableness, extraversion, openness to experience, conscientiousness, and neuroticism (12 items per factor). Items were answered on a 5-point Likert-type scale ranging from strongly disagree (score = 0) to strongly agree (score = 4).

Genotyping and Imputation

For all participants, peripheral blood DNA was extracted using Qiagen DNA isolation kits. Genotyping was conducted by HumanExome BeadChip (v1.1) according to the manufacturer’s instruction (Illumina Inc., San Diego, CA, United States). Genotypes were called with Illumina GenomeStudio software. Following genotyping, we conducted whole-genome imputation using IMPUTE2 (Howie et al., 2012) in 5-Mb chunks after prephasing with SHAPEIT2 (Delaneau et al., 2013). The haplotype panels released by the 1000 Genome Project (1000 Genomes Project Consortium et al., 2010) Phase31 (October 2014) were used as a reference.

For this study, SNPs located in the start and end positions of CNR1 plus 20-Kb regions on both sides were extracted (Chr6: 88829583-88896063, NCBI build 37). As shown in Supplementary Figure S1, most SNPs have a low minor allele frequency (MAF) (<0.01), and any SNP with a MAF of less than 0.01 were removed from the following analysis. The recommended cutoff of 0.3 for the “info” metric was used to filter poorly imputed variations (Marchini and Howie, 2010). After quality control, 26 imputed SNPs were included in the association analysis (Supplementary Table S1). Finally, we compared the MAF between imputed SNPs in this study and corresponding variants in 1000 Genome Project African population (Figure 1) and found them to be quite consistency of each other. A detailed genotyping information for all samples used in the study is provided in Supplementary Data Sheet S1.

FIGURE 1
www.frontiersin.org

FIGURE 1. Minor allele frequency (MAF) comparison between imputed SNPs in current study and SNPs in 1000 Genome Project African population.

Individual SNP-Based Association Analysis

We used the PLINK program (v. 1.07) (Purcell et al., 2007) to perform linear regression for each personality trait under an additive genetic model with age and sex included as covariates. Briefly, fcGENE (Roshyara and Scholz, 2014) was fitted into imputation workflow to convert genotype probability data to the PLINK file format with a tolerant cutoff value of 0.5. Standard quality controls for converted genotype data were implemented by PLINK.

To account for the probabilistic nature of imputed genotypes, frequentist tests, implemented by SNPTEST (v2.5.2) (Marchini et al., 2007), were used to evaluate the effects of CNR1 variants on personality traits. We applied two types of frequentist tests (Marchini and Howie, 2010). The Score test uses genotype likelihood to perform association tests in the way that each possible genotype is weighted by its imputation probability. The Expected test treats the expected allele count (also called allele dosage) as an independent variable and phenotype as dependent variable and then relates them by regression analysis. Sex and age were included as covariates in the association analysis. Based on the design for a candidate gene study, we adopted a relatively conservative Bonferroni correction for 26 SNPs, with a p-value less than 0.0019 (0.05/26) being considered statistically significant. Quanto v1.2.42 was employed to estimate the power of this study.

SNP Priority Score

For the identified SNPs, we calculated SNP priority scores by integrating evidence from single marker statistics, functional annotation, and previous studies. First, SNPs shown to be associated with any personality trait in this study were each assigned one point, which is accumulated if more than one trait is marked. Second, each SNP annotated by HaploReg v4.1 (Ward and Kellis, 2012) as exonic, splicing, 5′-UTR, 3′-UTR, promoter or enhancer was assigned one point. Third, SNPs reported to be associated with personality traits, psychiatric disorders, and emotional behaviors in other studies from the literature were highlighted again. Finally, we added all the scores as the SNP priority score as described previously (Imamura et al., 2016).

Results

Table 1 and Supplementary Figure S2 present the score statistics and distributions of five personality dimensions for all participants. All five dimensions showed a pattern of normal distribution. SNPs included in this report are shown in Supplementary Table S1, which spans from the 5′ promoter to 3′ downstream of the CNR1 gene region (Figure 2). We observed a very high consistency of allele frequency distribution between our imputed SNPs and those deposited in the dbSNP database (1000 Genomes Project Consortium et al., 2012) (R2 = 0.79, p = 9.13 × 10-10; Figure 1).

FIGURE 2
www.frontiersin.org

FIGURE 2. Schematic diagram of the human CNR1 and the positions of the SNPs investigated. Horizontal black arrows indicate the direction of transcription. The white bar represents the exonic region according to the National Center for Biotechnology Information (NCBI) and the University of California at Santa Cruz Browser (UCSC). The gray region within the white bar corresponds to the alternative splicing coding region, and the untranslated region is shown in black.

Association of Individual SNPs With Personality

Detailed association analysis of individual SNPs with five dimensions of personality is shown in Supplementary Tables S2–S6. Overall, we found that the p values estimated by the Score and Expected tests were highly concordant (R2 > 0.99), which is in line with a previous report that these two approaches show a great approximation when the effect size of risk alleles is small (Marchini and Howie, 2010). On the other hand, the frequentist test appeared to detect more significant results than PLINK did, indicating that genotype uncertainty represents a major factor with the imputed data.

For the case-control study on quantitative traits under the hypothesis of a gene-only model that directly relates phenotype to genotype, we have 94.97% power to detect a polymorphism (with an allele frequency > 1%) that could explain 1% of the variance of the traits at the 0.0001 two-tailed significance level with a sample size of 3,046, as used in this study.

For extraversion, 13 SNPs showed nominal significance (P < 0.05; see Supplementary Table S2). Of these significant SNPs, the association of rs806368 (p = 1.79 × 10-4), rs806370 (p = 2.68 × 10-4), rs806372 (p = 2.66 × 10-4), and rs2180619 (p = 1.53 × 10-3) remained significant after correction for multiple testing. To identify the potential significance-driven signal(s), we controlled the most significant SNP (rs806368) for the other three SNPs under investigation. A residual association was observed at SNP rs2180619 (p = 0.017), indicating that there might exist independent signal(s) within this region. Notably, all the significant SNPs showed negative associations with extraversion, indicating the CNR1 gene has a regulatory effect on the extraversion score.

For conscientiousness, the strongest association was detected on rs2180619 (p = 3.01 × 10-4) in the 5′ promoter region of CNR1 (Table 2 and Supplementary Table S3), an important region likely involved in regulating mRNA expression of the gene (Zhang et al., 2004). The other two SNPs significantly associated with conscientiousness were rs9353526 (p = 9.97 × 10-4) and rs9362466 (p = 1.36 × 10-3), both of which are located in the upstream region of CNR1. After conditioning on rs2180619, we found no evidence of any association signal for these two SNPs.

TABLE 2
www.frontiersin.org

TABLE 2. Identified SNPs associated with dimensions of personality.

For openness, five SNPs showed nominal significance (p < 0.05, Supplementary Table S4). Of them, SNPs rs806368, rs806370, and rs806372 also showed significantly association with extraversion (Table 2). It is worthwhile to note that SNP rs806368 has been reported to be significantly associated with openness in a sample with Caucasian origins (Juhasz et al., 2009a).

For agreeableness, we found four marginally associated loci, with the SNP rs2180619 being the strongest (p = 0.007, Supplementary Table S5). Importantly, these SNPs were detected only with PLINK, but not with SNPTEST, indicating a potentially less reliable result generated from missing genotypes. However, the agreeableness-associated SNPs were linked to extraversion (Table 2). For neuroticism, we did not observe any significance for the SNPs investigated (Supplementary Table S6).

Priority of Important SNPs Implicated in Personality Traits and Neuropsychiatry Disorders

Through integrating multiple lines of evidence for associated SNPs and calculating SNP priority score, we revealed that variants in CNR1 contributed to the genetic component of personality and their related psychiatric disorders, particularly for rs806368, rs806371, and rs2180619, whose priority scores equaled or were greater than 4 (range 1–9; Figure 3 and Supplementary Table S7).

FIGURE 3
www.frontiersin.org

FIGURE 3. Prioritized SNPs with current and previous evidence. Filled orange boxes indicate evidence from association analysis with nominal significance (P < 0.05). Filled red boxes indicate evidence from association analysis with multiple testing correction significance (P < 0.0019). Filled green boxes indicate evidence from functional annotation. Filled blue boxes indicate evidence from the published literature. SNP priority scores were counted by adding each piece of evidence.

Discussion

In the present study, we identified several variants in CNR1 showing significant associations with various dimensions of personality in an AA sample, especially for extraversion, for which four SNPs were identified. Importantly, SNPs associated with the other four dimensions of personality were all associated with extraversion. Considering that variants in CNR1 have been reported to be associated with an increased risk of mental health conditions, these data indicate the presence of genetic influences shared by personality traits and psychiatric disorders (Hettema et al., 2006; Kendler and Myers, 2010).

Our study might benefit greatly from the AA samples used in this study compared with those reported GWAS studies on various personality traits (de Moor et al., 2012; Smith et al., 2016), where samples with multiple ethnic origins, locations, or studies designed for different purposes were included. By using two approaches to analyze the same imputed dataset, we found that the frequentist test performed better than the regression model. When compared with PLINK in the context of genotype uncertainty, the frequentist tests, including the Score and Expected tests, improved statistical significance. One of the possible explanations for this result is that SNPTEST fully accounts for the uncertainty in imputed genotypes by weighting the imputation probability (Score test) or by calculating the allele dosage (Expected test). Conversely, the PLINK-based regression analysis uses only imputed genotypes that have a posterior probability above some threshold, which might lead to both false-positive findings and loss of power (Marchini and Howie, 2010). By implementing a gene-only model-based power analysis, we believe that the current genetic study is able to detect significant effects of CNR1 variants on personality traits in this valuable AA sample.

The relations between CNR1 and personality traits as well as other psychiatric disorders, such as depression and anxiety, have been established by Juhasz et al. (2009a), where 10 SNPs in CNR1 were genotyped, but they failed to detect any significant association with any personality trait at the individual SNP level in 1,269 participants. Interestingly, the haplotype trend regression analysis revealed significant associations with neuroticism and agreeableness after correction for multiple testing. Furthermore, those investigators detected SNP-by-SNP interactions for rs806379 with several other SNPs on both neuroticism and agreeableness. Another genetic association study revealed the association of rs7766029 with neuroticism (Aleksandrova et al., 2012).

Our individual SNP-based association analysis indicated significant association of CNR1 variants with extraversion in AAs. Rather than selecting haplotype-tagged SNPs, we imputed SNPs within the region of interest according to a high-density reference panel. An established statistical method was used in this study, which takes the genotype uncertainty into account and appeared to be more sensitive than PLINK in terms of our imputed data. Although almost all the investigated SNPs in the study reported by Juhasz et al. (2009a) were included in our analysis (except for rs7766029 because of the lack of reliable imputation quality), comparable results were not obtained in the two studies. We detected a positive association only between rs806368 and openness regardless of the genetic model. Such divergent results might be explained by the following factors: (Costa et al., 1992) different genetic models were used in these two studies. We detected single marker associations with the additive model, whereas Juhasz et al. (2009a) reported their results from three models (additive, dominant, and recessive); (McCrae et al., 2005) different questionnaires were used in these two studies. Juhasz et al. (2009a) used the Big Five Inventory (BFI-44) to assess personality, whereas we used the NEO Personality Inventory. Although there exist excellent correlations between the BFI-44 Inventory and NEO, a slight discrepancy might arise in phenotypic assessments; and (DeYoung et al., 2010) population stratification: we observed differential allele frequencies between AAs and European-American (EA) populations (Supplementary Table S8). For example, the allele frequency of rs1049359 is rarer in AAs (MAF = 0.029) than in EAs (MAF = 0.258).

The extraversion-associated SNP rs806368 remained significant after Bonferroni correction. In addition, this SNP showed a marginal association with both agreeableness and openness (Table 2). Rs806368 is located in the 3′-UTR region of CNR1. A strong association between rs306368 and impulsivity has been reported for the Southwest California Indians (also known as Mission Indians) (Ehlers et al., 2007). Impulsivity is considered as a facet of personality and acts as a major component of various psychiatric diseases such as ADHD (Nigg, 2001), bipolar disorder (Henry et al., 2001), and antisocial personality disorder (Horn et al., 2003). The study reported by Eysenck (1967) indicates that extraversion is a combination of impulsiveness and sociability, which is consistent with the genetic findings that rs806368 is associated with extraversion and impulsivity.

A pharmacogenomics study highlighted the significant role of rs806371 in both MDD etiology and the clinical response to citalopram treatment (Mitjans et al., 2013). Mitjans et al. (2013) reported a higher frequency of rs806371 G carriers in MDD patients with melancholia and psychotic symptoms than in controls. Moreover, the C allele carriers of rs806371 demonstrate increased remission status of MD symptoms after 12 weeks of treatment with citalopram (Mitjans et al., 2013).

Another interesting SNP revealed by this study is rs2180619, located in the promoter region of CNR1, which is not only significantly associated with both extraversion and conscientiousness, but also nominally with agreeableness (Table 2). By genotyping 706 individuals for the 5′-HTTLPR in the SLC6A4 and four SNPs in CNR1, Lazary et al. (2009) showed that the risk of anxiety was 4.6-fold greater in individuals with genotypes G/G for rs2180619 and S/S for 5-HTTLPR compared with other genotypes. Furthermore, Heitland et al. (2012) found that individuals with genotypes G/G or G/A for rs2180619 showed robust extinction of fear compared with individuals with the A/A genotype. Interestingly, Mroczkowski et al. (2016) found that a high neuroticism score significantly increases the risk of three anxiety disorders (agoraphobia, panic disorder, and social phobia), and extraversion and conscientiousness were negatively associated with separation anxiety disorder and social phobia. Although these relationships were constructed in individuals with obsessive compulsive disorder, the significances survived after adjusting for age at onset of this disorder (Mroczkowski et al., 2016).

Relationships among personality traits are worthy of consideration. For example, extraversion and openness shared genetic overlap in the context of rs806368, rs806370, rs806371, rs806372, and rs9444584; extraversion, conscientiousness, and agreeableness were highly related for rs9353526, rs9362466, and rs2180619 (Table 2). We further found the direction of effects for associated SNPs were opposite between neuroticism and the other four personality traits. It was reported recently that neuroticism shows a negative genetic correlation with the other four personality traits (Lo et al., 2017). Taken together, the multiple dimensions of personality-associated CNR1 variants support the concept that the presence of shared genetic factors contributes to the correlation among personalities (Lo et al., 2017), suggesting a fundamental role of CNR1 in personality formation.

Dysfunction of CNR1 may result in neuropsychological disturbances, leading to psychiatric disorders. Martin et al. (2002) first used CB1 knockout mice to evaluate the potential role of CNR1 in emotional behavior (Ledent et al., 1999). By using a light/dark box, they observed a significant decrease in the number of entries and in the time spent in the lit compartment for CB1 knockout mice, indicating that these mice present an anxiogenic-like response. By using a chronic unpredictable mild stress model (CMS) to simulate anhedonia, the CB1 knockout mice exhibited greater anhedonia, a core symptom of depression and defining feature of melancholia (Klonsky, 2008; Battle, 2013). More recently, Bowers and Ressler (2016) administered Cnr1 antagonist SR141716A to C57BL/6J mice, after which both male and female adult mice were sensitive to Cnr1 antagonist-mediated increases in anxiety-like behavior, thus supporting the protective role of CNR1 against the development of anxiety disorder. These studies confirmed that dysfunction of the CB1 receptor changes several behavioral responses and increases susceptibility to emotional disorders, such as anxiety, depression, and aggressiveness.

There are limitations associated with this study. First, the SNPs investigated were generated primarily by imputation. Although we have great confidence in the accuracy of our imputed genotypes, as demonstrated in the literature, including our own work (Marchini and Howie, 2010; Howie et al., 2011), individual genotyping of variants in CNR1 is desirable in future studies to support the conclusions about the association of CNR1 variants with personality in AAs. Second, our imputation-based SNP selections provide a fine-mapping of CNR1, but the SNPs examined could not accurately tag the causal genetic variants. Last but not least, the detected associations of variants in CNR1 with various personality traits require to be replicated with independent AA samples. However, such type of AA samples with appropriate phenotypes is unavailable to us.

Conclusion

This is the first study to demonstrate the significant effects of CNR1 polymorphisms on personality traits in AAs. The SNPs identified in this study revealed significant associations with multiple personality traits, indicating the fundamental role of CNR1 variants in the expression of personality. Furthermore, previous publications have related our identified personality-associated SNPs with psychiatric disorders, which provides additional evidence supporting the presence of genetic correlation between personalities and mental health conditions. Understanding these relationships may reveal opportunities to prevent and treat neuropsychiatric disorders.

Author Contributions

JM, TP, and ML participated in clinical data collection. YY, JZ, WY, KS, and YM performed the laboratory experiments. YY, YX, KS, and YM conducted the data analysis. YY, TP, and ML participated in writing of the manuscript. ML conceived the study and involved in every step of the work. All authors approved the final draft of the paper.

Funding

This study was supported in part by the China Precision Medicine Initiative (2016YFC0906300), Research Center for Air Pollution and Health of Zhejiang University, the State Key Laboratory for Diagnosis and Treatment of Infectious Diseases of the First Affiliated Hospital of Zhejiang University, and U.S. National Institutes of Health grant DA012844.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

We thank Dr. David L. Bronson for excellent editing of this manuscript.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2018.00199/full#supplementary-material

Footnotes

  1. ^http://mathgen.stats.ox.ac.uk/impute/1000GP_Phase3.html
  2. ^http://biostats.usc.edu/Quanto.html

References

1000 Genomes Project Consortium, Abecasis, G. R., Altshuler, D., Auton, A., Brooks, L. D., Durbin, R. M., et al. (2010). A map of human genome variation from population-scale sequencing. Nature 467, 1061–1073. doi: 10.1038/nature09534

PubMed Abstract | CrossRef Full Text | Google Scholar

1000 Genomes Project Consortium, Abecasis, G. R., Auton, A., Brooks, L. D., DePristo, M. A., Durbin, R. M., et al. (2012). An integrated map of genetic variation from 1,092 human genomes. Nature 491, 56–65. doi: 10.1038/nature11632

PubMed Abstract | CrossRef Full Text | Google Scholar

Aleksandrova, L., Souza, R., Bagby, M., Casey, D., Hodgins, D., Smith, G., et al. (2012). Genetic underpinnings of neuroticism: a replication study. J. Addict. Res. Ther. 3, 1–4. doi: 10.4172/2155-6105.1000119

CrossRef Full Text | Google Scholar

Amin, N., Hottenga, J. J., Hansell, N. K., Janssens, A. C., de Moor, M. H., Madden, P. A., et al. (2013). Refining genome-wide linkage intervals using a meta-analysis of genome-wide association studies identifies loci influencing personality dimensions. Eur. J. Hum. Genet. 21, 876–882. doi: 10.1038/ejhg.2012.263

PubMed Abstract | CrossRef Full Text | Google Scholar

Amin, N., Schuur, M., Gusareva, E. S., Isaacs, A., Aulchenko, Y. S., Kirichenko, A. V., et al. (2012). Agenome-wide linkage study of individuals with high scores on NEO personality traits. Mol. Psychiatry 17, 1031–1041. doi: 10.1038/mp.2011.97

PubMed Abstract | CrossRef Full Text | Google Scholar

Barnett, J. H., Huang, J., Perlis, R. H., Young, M. M., Rosenbaum, J. F., Nierenberg, A. A., et al. (2011). Personality and bipolar disorder: dissecting state and trait associations between mood and personality. Psychol. Med. 41, 1593–1604. doi: 10.1017/S0033291710002333

PubMed Abstract | CrossRef Full Text | Google Scholar

Battle, D. E. (2013). Diagnostic and statistical manual of mental disorders (DSM). Codas 25, 191–192.

Google Scholar

Bowers, M. E., and Ressler, K. J. (2016). Sex-dependence of anxiety-like behavior in cannabinoid receptor 1 (Cnr1) knockout mice. Behav. Brain Res. 300, 65–69. doi: 10.1016/j.bbr.2015.12.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Boyle, E. A., Li, Y. I., and Pritchard, J. K. (2017). An expanded view of complex traits: from polygenic to omnigenic. Cell 169, 1177–1186. doi: 10.1016/j.cell.2017.05.038

PubMed Abstract | CrossRef Full Text | Google Scholar

Calboli, F. C., Tozzi, F., Galwey, N. W., Antoniades, A., Mooser, V., Preisig, M., et al. (2010). A genome-wide association study of neuroticism in a population-based sample. PLoS One 5:e11504. doi: 10.1371/journal.pone.0011504

PubMed Abstract | CrossRef Full Text | Google Scholar

Choi, J. S., Payne, T. J., Ma, J. Z., and Li, M. D. (2017). Relationship between personality traits and nicotine dependence in male and female smokers of African-American and European-American samples. Front. Psychiatry 8:122. doi: 10.3389/fpsyt.2017.00122

PubMed Abstract | CrossRef Full Text | Google Scholar

Costa, P. T. Jr., and McCrae, R. R. (1997). Stability and change in personality assessment: the revised NEO Personality Inventory in the year 2000. J. Pers. Assess. 68, 86–94. doi: 10.1207/s15327752jpa6801_7

PubMed Abstract | CrossRef Full Text | Google Scholar

Costa, P. T., McCrea, R. R., and Psychological Assessment Resources Inc. (1992). Revised NEO Personality Inventory (NEO PI-R) and NEO Five-Factor Inventory (NEO-FFI), Vol. VI. Odessa, FL: Psychological Assessment Resources, 101.

Google Scholar

de Moor, M. H., Costa, P. T., Terracciano, A., Krueger, R. F., de Geus, E. J., Toshiko, T., et al. (2012). Meta-analysis of genome-wide association studies for personality. Mol. Psychiatry 17, 337–349. doi: 10.1038/mp.2010.128

PubMed Abstract | CrossRef Full Text | Google Scholar

Delaneau, O., Zagury, J. F., and Marchini, J. (2013). Improved whole-chromosome phasing for disease and population genetic studies. Nat. Methods 10, 5–6. doi: 10.1038/nmeth.2307

PubMed Abstract | CrossRef Full Text | Google Scholar

DeYoung, C. G., Hirsh, J. B., Shane, M. S., Papademetris, X., Rajeevan, N., and Gray, J. R. (2010). Testing predictions from personality neuroscience. Brain structure and the big five. Psychol. Sci. 21, 820–828. doi: 10.1177/0956797610370159

PubMed Abstract | CrossRef Full Text | Google Scholar

Ehlers, C. L., Slutske, W. S., Lind, P. A., and Wilhelmsen, K. C. (2007). Association between single nucleotide polymorphisms in the cannabinoid receptor gene (CNR1) and impulsivity in southwest California Indians. Twin Res. Hum. Genet. 10, 805–811. doi: 10.1375/twin.10.6.805

PubMed Abstract | CrossRef Full Text | Google Scholar

Elphick, M. R., and Egertova, M. (2001). The neurobiology and evolution of cannabinoid signalling. Philos. Trans. R. Soc. Lond. Ser. B Biol. Sci. 356, 381–408. doi: 10.1098/rstb.2000.0787

PubMed Abstract | CrossRef Full Text | Google Scholar

Eysenck, H. J. (1967). The Biological Basis of Personality, Vol. XVII. Springfield, IL: Thomas, 399.

Google Scholar

Floderus-Myrhed, B., Pedersen, N., and Rasmuson, I. (1980). Assessment of heritability for personality, based on a short-form of the Eysenck Personality Inventory: a study of 12,898 twin pairs. Behav. Genet. 10, 153–162. doi: 10.1007/BF01066265

PubMed Abstract | CrossRef Full Text | Google Scholar

Genetics of Personality Consortium, de Moor, M. H., van den Berg, S. M., Verweij, K. J., Krueger, R. F., Luciano, M., et al. (2015). Meta-analysis of genome-wide association studies for neuroticism, and the polygenic association with major depressive disorder. JAMA Psychiatry 72, 642–650. doi: 10.1001/jamapsychiatry.2015.0554

PubMed Abstract | CrossRef Full Text | Google Scholar

Heitland, I., Klumpers, F., Oosting, R. S., Evers, D. J., Leon Kenemans, J., and Baas, J. M. (2012). Failure to extinguish fear and genetic variability in the human cannabinoid receptor 1. Transl. Psychiatry 2:e162. doi: 10.1038/tp.2012.90

PubMed Abstract | CrossRef Full Text | Google Scholar

Henry, C., Mitropoulou, V., New, A. S., Koenigsberg, H. W., Silverman, J., and Siever, L. J. (2001). Affective instability and impulsivity in borderline personality and bipolar II disorders: similarities and differences. J. Psychiatr. Res. 35, 307–312. doi: 10.1016/S0022-3956(01)00038-3

PubMed Abstract | CrossRef Full Text | Google Scholar

Hettema, J. M., Neale, M. C., Myers, J. M., Prescott, C. A., and Kendler, K. S. (2006). A population-based twin study of the relationship between neuroticism and internalizing disorders. Am. J. Psychiatry 163, 857–864. doi: 10.1176/ajp.2006.163.5.857

PubMed Abstract | CrossRef Full Text | Google Scholar

Horn, J. M., Plomin, R., and Rosenman, R. (1976). Heritability of personality traits in adult male twins. Behav. Genet. 6, 17–30. doi: 10.1007/BF01065675

CrossRef Full Text | Google Scholar

Horn, N. R., Dolan, M., Elliott, R., Deakin, J. F., and Woodruff, P. W. (2003). Response inhibition and impulsivity: an fMRI study. Neuropsychologia 41, 1959–1966. doi: 10.1016/S0028-3932(03)00077-0

CrossRef Full Text | Google Scholar

Howie, B., Fuchsberger, C., Stephens, M., Marchini, J., and Abecasis, G. R. (2012). Fast and accurate genotype imputation in genome-wide association studies through pre-phasing. Nat. Genet. 44, 955–959. doi: 10.1038/ng.2354

PubMed Abstract | CrossRef Full Text | Google Scholar

Howie, B., Marchini, J., and Stephens, M. (2011). Genotype imputation with thousands of genomes. G3 1, 457–470. doi: 10.1534/g3.111.001198

PubMed Abstract | CrossRef Full Text | Google Scholar

Imamura, M., Takahashi, A., Yamauchi, T., Hara, K., Yasuda, K., Grarup, N., et al. (2016). Genome-wide association studies in the Japanese population identify seven novel loci for type 2 diabetes. Nat. Commun. 7:10531. doi: 10.1038/ncomms10531

PubMed Abstract | CrossRef Full Text | Google Scholar

Jang, K. L., Livesley, W. J., and Vernon, P. A. (1996). Heritability of the big five personality dimensions and their facets: a twin study. J. Pers. 64, 577–591. doi: 10.1111/j.1467-6494.1996.tb00522.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Juhasz, G., Chase, D., Pegg, E., Downey, D., Toth, Z. G., Stones, K., et al. (2009a). CNR1 gene is associated with high neuroticism and low agreeableness and interacts with recent negative life events to predict current depressive symptoms. Neuropsychopharmacology 34, 2019–2027. doi: 10.1038/npp.2009.19

PubMed Abstract | CrossRef Full Text | Google Scholar

Juhasz, G., Csepany, E., Magyar, M., Edes, A. E., Eszlari, N., Hullam, G., et al. (2017). Variants in the CNR1 gene predispose to headache with nausea in the presence of life stress. Genes Brain Behav. 16, 384–393. doi: 10.1111/gbb.12352

PubMed Abstract | CrossRef Full Text | Google Scholar

Juhasz, G., Lazary, J., Chase, D., Pegg, E., Downey, D., Toth, Z. G., et al. (2009b). Variations in the cannabinoid receptor 1 gene predispose to migraine. Neurosci. Lett. 461, 116–120. doi: 10.1016/j.neulet.2009.06.021

PubMed Abstract | CrossRef Full Text | Google Scholar

Kendler, K. S., and Myers, J. (2010). The genetic and environmental relationship between major depression and the five-factor model of personality. Psychol. Med. 40, 801–806. doi: 10.1017/S0033291709991140

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, B. H., Kim, H. N., Roh, S. J., Lee, M. K., Yang, S., Lee, S. K., et al. (2015). GWA meta-analysis of personality in Korean cohorts. J. Hum. Genet. 60, 455–460. doi: 10.1038/jhg.2015.52

PubMed Abstract | CrossRef Full Text | Google Scholar

Kim, H. N., Roh, S. J., Sung, Y. A., Chung, H. W., Lee, J. Y., Cho, J., et al. (2013). Genome-wide association study of the five-factor model of personality in young Korean women. J. Hum. Genet. 58, 667–674. doi: 10.1038/jhg.2013.75

PubMed Abstract | CrossRef Full Text | Google Scholar

Klonsky, E. D. (2008). What is emptiness? Clarifying the 7th criterion for borderline personality disorder. J. Pers. Disord. 22, 418–426. doi: 10.1521/pedi.2008.22.4.418

PubMed Abstract | CrossRef Full Text | Google Scholar

Lazary, J., Lazary, A., Gonda, X., Benko, A., Molnar, E., Hunyady, L., et al. (2009). Promoter variants of the cannabinoid receptor 1 gene (CNR1) in interaction with 5-HTTLPR affect the anxious phenotype. Am. J. Med. Genet. B Neuropsychiatr. Genet. 150B, 1118–1127. doi: 10.1002/ajmg.b.31024

PubMed Abstract | CrossRef Full Text | Google Scholar

Ledent, C., Valverde, O., Cossu, G., Petitet, F., Aubert, J. F., Beslot, F., et al. (1999). Unresponsiveness to cannabinoids and reduced addictive effects of opiates in CB1 receptor knockout mice. Science 283, 401–404. doi: 10.1126/science.283.5400.401

PubMed Abstract | CrossRef Full Text | Google Scholar

Lo, M. T., Hinds, D. A., Tung, J. Y., Franz, C., Fan, C. C., Wang, Y., et al. (2017). Genome-wide analyses for personality traits identify six genomic loci and show correlations with psychiatric disorders. Nat. Genet. 49, 152–156. doi: 10.1038/ng.3736

PubMed Abstract | CrossRef Full Text | Google Scholar

Lu, A. T., Ogdie, M. N., Jarvelin, M. R., Moilanen, I. K., Loo, S. K., McCracken, J. T., et al. (2008). Association of the cannabinoid receptor gene (CNR1) with ADHD and post-traumatic stress disorder. Am. J. Med. Genet. B Neuropsychiatr. Genet. 147B, 1488–1494. doi: 10.1002/ajmg.b.30693

PubMed Abstract | CrossRef Full Text | Google Scholar

Luciano, M., Hagenaars, S. P., Davies, G., Hill, W. D., Clarke, T. K., Shirali, M., et al. (2018). Association analysis in over 329,000 individuals identifies 116 independent variants influencing neuroticism. Nat. Genet. 50, 6–11. doi: 10.1038/s41588-017-0013-8

PubMed Abstract | CrossRef Full Text | Google Scholar

Manolio, T. A., Collins, F. S., Cox, N. J., Goldstein, D. B., Hindorff, L. A., Hunter, D. J., et al. (2009). Finding the missing heritability of complex diseases. Nature 461, 747–753. doi: 10.1038/nature08494

PubMed Abstract | CrossRef Full Text | Google Scholar

Marchini, J., and Howie, B. (2010). Genotype imputation for genome-wide association studies. Nat. Rev. Genet. 11, 499–511. doi: 10.1038/nrg2796

PubMed Abstract | CrossRef Full Text | Google Scholar

Marchini, J., Howie, B., Myers, S., McVean, G., and Donnelly, P. (2007). A new multipoint method for genome-wide association studies by imputation of genotypes. Nat. Genet. 39, 906–913. doi: 10.1038/ng2088

PubMed Abstract | CrossRef Full Text | Google Scholar

Martin, M., Ledent, C., Parmentier, M., Maldonado, R., and Valverde, O. (2002). Involvement of CB1 cannabinoid receptors in emotional behaviour. Psychopharmacology 159, 379–387. doi: 10.1007/s00213-001-0946-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Martinez-Gras, I., Hoenicka, J., Ponce, G., Rodriguez-Jimenez, R., Jimenez-Arriero, M. A., Perez-Hernandez, E., et al. (2006). (AAT)n repeat in the cannabinoid receptor gene, CNR1: association with schizophrenia in a Spanish population. Eur. Arch. Psychiatry Clin. Neurosci. 256, 437–441. doi: 10.1007/s00406-006-0665-3

PubMed Abstract | CrossRef Full Text | Google Scholar

McCrae, R. R., Costa, P. T. Jr., and Martin, T. A. (2005). The NEO-PI-3: a more readable revised NEO Personality Inventory. J. Pers. Assess. 84, 261–270. doi: 10.1207/s15327752jpa8403_05

PubMed Abstract | CrossRef Full Text | Google Scholar

Mitjans, M., Serretti, A., Fabbri, C., Gasto, C., Catalan, R., Fananas, L., et al. (2013). Screening genetic variability at the CNR1 gene in both major depression etiology and clinical response to citalopram treatment. Psychopharmacology 227, 509–519. doi: 10.1007/s00213-013-2995-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Mroczkowski, M. M., Goes, F. S., Riddle, M. A., Grados, M. A., Bienvenu, O. J., Greenberg, B. D., et al. (2016). Dependent personality, separation anxiety disorder and other anxiety disorders in OCD. Pers. Ment. Health 10, 22–28. doi: 10.1002/pmh.1321

PubMed Abstract | CrossRef Full Text | Google Scholar

Nigg, J. T. (2001). Is ADHD a disinhibitory disorder? Psychol. Bull. 127, 571–598.

Google Scholar

Okbay, A., Baselmans, B. M., De Neve, J. E., Turley, P., Nivard, M. G., Fontana, M. A., et al. (2016). Genetic variants associated with subjective well-being, depressive symptoms, and neuroticism identified through genome-wide analyses. Nat. Genet. 48, 624–633. doi: 10.1038/ng.3552

PubMed Abstract | CrossRef Full Text | Google Scholar

Pardo, L. M., MacKay, I., Oostra, B., van Duijn, C. M., and Aulchenko, Y. S. (2005). The effect of genetic drift in a young genetically isolated population. Ann. Hum. Genet. 69, 288–295. doi: 10.1046/j.1529-8817.2005.00162.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Purcell, S., Neale, B., Todd-Brown, K., Thomas, L., Ferreira, M. A., Bender, D., et al. (2007). PLINK: a tool set for whole-genome association and population-based linkage analyses. Am. J. Hum. Genet. 81, 559–575.

PubMed Abstract | Google Scholar

Roshyara, N. R., and Scholz, M. (2014). fcGENE: a versatile tool for processing and transforming SNP datasets. PLoS One 9:e97589. doi: 10.1371/journal.pone.0097589

PubMed Abstract | CrossRef Full Text | Google Scholar

Smith, D. J., Escott-Price, V., Davies, G., Bailey, M. E., Colodro-Conde, L., Ward, J., et al. (2016). Genome-wide analysis of over 106 000 individuals identifies 9 neuroticism-associated loci. Mol. Psychiatry 21, 749–757. doi: 10.1038/mp.2016.49

PubMed Abstract | CrossRef Full Text | Google Scholar

Terracciano, A., and Costa, P. T. Jr. (2004). Smoking and the Five-Factor Model of personality. Addiction 99, 472–481. doi: 10.1111/j.1360-0443.2004.00687.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Terracciano, A., Lockenhoff, C. E., Crum, R. M., Bienvenu, O. J., and Costa, P. T. Jr. (2008a). Five-Factor Model personality profiles of drug users. BMC Psychiatry 8:22. doi: 10.1186/1471-244X-8-22

PubMed Abstract | CrossRef Full Text | Google Scholar

Terracciano, A., Lockenhoff, C. E., Zonderman, A. B., Ferrucci, L., and Costa, P. T. Jr. (2008b). Personality predictors of longevity: activity, emotional stability, and conscientiousness. Psychosom. Med. 70, 621–627. doi: 10.1097/PSY.0b013e31817b9371

PubMed Abstract | CrossRef Full Text | Google Scholar

Ujike, H., Takaki, M., Nakata, K., Tanaka, Y., Takeda, T., Kodama, M., et al. (2002). CNR1, central cannabinoid receptor gene, associated with susceptibility to hebephrenic schizophrenia. Mol. Psychiatry 7, 515–518. doi: 10.1038/sj.mp.4001029

PubMed Abstract | CrossRef Full Text | Google Scholar

van den Berg, S. M., de Moor, M. H., Verweij, K. J., Krueger, R. F., Luciano, M., Arias Vasquez, A., et al. (2016). Meta-analysis of genome-wide association studies for extraversion: findings from the genetics of personality consortium. Behav. Genet. 46, 170–182. doi: 10.1007/s10519-015-9735-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Ward, L. D., and Kellis, M. (2012). HaploReg: a resource for exploring chromatin states, conservation, and regulatory motif alterations within sets of genetically linked variants. Nucleic Acids Res. 40, D930–D934. doi: 10.1093/nar/gkr917

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, J., Benyamin, B., McEvoy, B. P., Gordon, S., Henders, A. K., Nyholt, D. R., et al. (2010). Common SNPs explain a large proportion of the heritability for human height. Nat. Genet. 42, 565–569. doi: 10.1038/ng.608

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, J., Wang, S., Yang, Z., Hodgkinson, C. A., Iarikova, P., Ma, J. Z., et al. (2015). The contribution of rare and common variants in 30 genes to risk nicotine dependence. Mol. Psychiatry 20, 1467–1478. doi: 10.1038/mp.2014.156

PubMed Abstract | CrossRef Full Text | Google Scholar

Yu, W., De Hert, M., Moons, T., Claes, S. J., Correll, C. U., and van Winkel, R. (2013). CNR1 gene and risk of the metabolic syndrome in patients with schizophrenia. J. Clin. Psychopharmacol. 33, 186–192. doi: 10.1097/JCP.0b013e318283925e

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, P. W., Ishiguro, H., Ohtsuki, T., Hess, J., Carillo, F., Walther, D., et al. (2004). Human cannabinoid receptor 1: 5′ exons, candidate regulatory regions, polymorphisms, haplotypes and association with polysubstance abuse. Mol. Psychiatry 9, 916–931. doi: 10.1038/sj.mp.4001560

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: CNR1, SNP, African-Americans, NEO-FFI, personality traits

Citation: Yao Y, Xu Y, Zhao J, Ma Y, Su K, Yuan W, Ma JZ, Payne TJ and Li MD (2018) Detection of Significant Association Between Variants in Cannabinoid Receptor 1 Gene (CNR1) and Personality in African–American Population. Front. Genet. 9:199. doi: 10.3389/fgene.2018.00199

Received: 13 September 2017; Accepted: 17 May 2018;
Published: 14 June 2018.

Edited by:

Francis J. McMahon, Intramural Research Program (NIH), United States

Reviewed by:

Radka Kaneva, Medical University, Sofia, Bulgaria
Sarah E. Bergen, Karolinska Institutet, Sweden

Copyright © 2018 Yao, Xu, Zhao, Ma, Su, Yuan, Ma, Payne and Li. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Ming D. Li, ml2km@zju.edu.cn; ming.li@shu.edu; limd586@outlook.com

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.