The Cannabinoid Receptor 2 Q63R Variant Modulates the Relationship between Childhood Obesity and Age at Menarche

Abstract

Background

The ovary is an important site where gene variants modulate pubertal timing. The cannabinoid receptor 2 (CB2) is expressed in the ovary, plays a role in folliculogenesis and ovulation, and can be modulated by estrogens. Obesity is strictly associated with early menarche and is characterized by sex hormone and endocannabinoid derangement.

Aim

In this study, we investigated the role of the CB2 receptor in determining the age at menarche in obese girls.

Methods

We studied a cohort of 240 obese girls (age 11.9±3 years; BMI z-score 2.8±0.8). The age at menarche (if it had already occurred) was recorded at the time of the visit or via phonecall. The CNR2 rs35761398 polymorphism, which leads to the CB2 Q63R variant, was detected by the TaqMan assay.

Results

In total, 105 patients were homozygous for the R63-coding allele (RR), 113 were QR and 22 were QQ. Variance analysis revealed a significantly earlier age of menarche in subjects carrying the Q63 allele, which was also found after adjusting for BMI z-score (11±1.2 vs. 11.6±1.2 years, p = 0.0003). Logistic regression analysis demonstrated that patients homozygous for the Q allele had a 2.2-fold higher risk (odds ratio = 2.2; CI1.1–3.4; p = 0.02) of presenting with an early menarche (age at menarche <12 years).

Conclusion

We demonstrated for the first time the association between the CB2 Q63R functional variant and the age at menarche in a cohort of Italian obese girls.

About the Authors:

Giulia Bellini

Contributed equally to this work with: Giulia Bellini, Anna Grandone

Affiliation: Department of Experimental Medicine, Pharmacology Division, Second University of Naples, Naples, Italy

Anna Grandone

Contributed equally to this work with: Giulia Bellini, Anna Grandone

Affiliation: Department of Woman, Child and of General and Specialist Surgery, Second University of Naples, Naples, Italy

Marco Torella

Affiliation: Department of Woman, Child and of General and Specialist Surgery, Second University of Naples, Naples, Italy

Emanuele Miraglia del Giudice

Affiliation: Department of Woman, Child and of General and Specialist Surgery, Second University of Naples, Naples, Italy

Bruno Nobili

Affiliation: Department of Woman, Child and of General and Specialist Surgery, Second University of Naples, Naples, Italy

Laura Perrone

Affiliation: Department of Woman, Child and of General and Specialist Surgery, Second University of Naples, Naples, Italy

Sabatino Maione

Affiliation: Department of Experimental Medicine, Pharmacology Division, Second University of Naples, Naples, Italy

Francesca Rossi

* E-mail: [email protected]

Affiliation: Department of Woman, Child and of General and Specialist Surgery, Second University of Naples, Naples, Italy

Introduction

The age at menarche is a marker of pubertal timing in females. Pubertal timing is widely variable, a result of the interactions of both environmental and genetic determinants. The age at menarche is associated with obesity, type 2 diabetes, cardiovascular disease, breast cancer and all-cause mortality and is characterized by a complex genetic architecture [1–3]. The mechanisms that determine pubertal timing and underlie its links to disease risks remain unclear.

In the ovary, some gene variants such as the lin-28 homolog B (LIN28B) gene polymorphism have been strongly associated with the age at menarche, thus modulating pubertal timing [4,5]. Additionally, early menarche has been associated with childhood obesity [6].

Multiple endocrine pathways link childhood obesity to early pubertal onset [7–10].

The endocannabinoid (EC) system includes the EC endogenous lipid transmitters, their G-protein-coupled cannabinoid receptors type 1 and type 2 (CB1 and CB2), and the enzymes for EC biosynthesis and degradation. ECs are synthesized and immediately released to activate CB1 and CB2 “on demand”, although they can also accumulate in intracellular adiposomes [11].

Recent data suggest that the EC system plays a role in folliculogenesis and ovulation and is widely expressed in the ovarian medulla and cortex. The CB2 receptor is preferentially expressed with respect to CB1 in the granulosa cells (GCs) of follicles at all stages of development as well as in the corpus luteum and corpus albicans. CB2 stimulation induces follicle growth and oocyte maturation [12, 13]. Moreover, the EC system is disturbed in obesity, resulting in an enhanced tone that is associated with increased EC plasma levels [14].

Based on this evidence, we are tempted to speculate that a different CB2 receptor response could contribute to early menarche in obesity. Thus, we investigated a functional variant of the CB2 receptor, Q63R, for its potential influence on the age at menarche in a cohort of 240 Italian obese girls.

Materials and Methods

Subjects

We studied a cohort of 240 obese girls. Obesity was defined as a BMI exceeding the 95th percentile for age and sex. Patients with autoimmune disorders were excluded due to the previously described association between the CB2 Q63R variant and autoimmune disorders [15–18].

Assent from the girls and written informed consent from their parents was obtained. The study was conducted with the approval of the Ethic Committee of the Second University of Naples in compliance with the national legislation and the Declaration of Helsinki.

Assessment of menarchal age was obtained by the status quo method when recorded at the time of the visit if it had already occurred and by the recall or retrospective method via phone call if it had not yet occurred.

Genotyping

Genomic DNA was extracted from peripheral whole blood with a DNA extraction kit (Roche Diagnostics, Branchburg, NJ, USA). Molecular screening for the CNR2 rs35761398 polymorphism (CAA/CGG) underlying the CB2 Q63R substitution was performed using a TaqMan assay (Real Master Mix Probe, 5 PRIME, Germany). The primers and probes used were the following: sense 5'-GTGCTCTATCTGATCCTGTC-3' and anti-sense 5'-TAGTCACGCTGCCAATC-3'; AA-probe 5'-CCCACCAACTCCGC-3' and GG-probe 5'-CCCACCGGCTCCG-3' (PRIMM, Milan, Italy). Both PCR and post-PCR allelic discrimination were performed on a7900 HT Fast System Thermal Cycler (Applied Biosystems, Foster City, CA, USA).

Genotypes of random samples were confirmed by PCR followed by direct sequencing. The PCR program consisted of 94°C for 4 min followed by 31cycles of 94°C for 30 s, 60°C for 30 s and 72°C for 30 s. Primers were the following: forward 5'-GAGTGGTCCCCAGAAGACAG-3' and reverse 5'-CACAGAGGCTGTGAAGGTCA-3'. PCR products were sequenced using an ABI PRISM 9600 automated sequencer (Applied Biosystems, Foster City, CA, USA) and the Big Dye Terminator reaction kit (Applera, Foster City, CA, USA) according to the manufacturer’s instructions. All of the primers were chosen using Primer3 software (http://primer3.sourceforge.net/).

Statistics

Differences between categorical variables were analyzed using a chi-squared test. A linear logistic regression was performed to analyze clinical data with respect to the CB2 Q63R variant. A p value < 0.05 was considered to be statistically significant. All of the analyses were performed using StatGraphics CENTURION XV.II (Adalta, Arezzo, Italy; STATPOINT TECHNOLOGIES INC., Virginia, USA).

Results

We analyzed a cohort of 240 obese girls(age 12.9±2.6 years; BMI z-score 2.8±0.8) for the CNR2 rs35761398 polymorphism (CB2 Q63R variant) by the TaqMan assay. Clinical data are summarized in Table 1.

[Figure omitted. See PDF.]

Table 1. Clinical data of 240 obese Italian girls according to CB2 Q63R variant.

https://doi.org/10.1371/journal.pone.0140142.t001

In total, 105 patients were homozygous for the R-encoding allele (R63), 113 were heterozygous (Q63R) and 22 were homozygous for the Q-encoding allele (Q63).

Allele frequencies were in Hardy Weinberg equilibrium (p = 0.28).

Variance analysis demonstrated a significantly earlier age of menarche in subjects carrying the Q allele also after adjusting for the BMI z-score (p = 0.0016). No differences in metabolic parameters were observed (Table 1). Logistic regression analysis demonstrated that Q63 subjects had a 2.05-fold higher risk of presenting menarche before 12 years of age (Odds Ratio (OR) = 2.05; CI 1.21–3.45; p = 0.0068).

After stratifying by BMI z-score, the differences in the age at menarche among genotypes were not significant among patients in the lowest BMI z-score tertile (Fig 1).

[Figure omitted. See PDF.]

Fig 1. Differences in age at menarche according to BMI z-score tertiles.

*p<0.001.

https://doi.org/10.1371/journal.pone.0140142.g001

Discussion

In this study, we demonstrated, for the first time, an association between the rs35761398 polymorphism in the CNR2 gene, which underlies the Q63R functional variant, and the age at menarche in a cohort of 240 Italian obese girls. CB2 differentially modulates its effector cells and downstream pathways depending on the presence of glutamine or arginine at codon 63 of the N-terminal domain; thus, the Q63 variant is more functional than the R63 mutant. Indeed, lymphocytes derived from RR-subjects display minor proliferation inhibition compared with those derived from QQ subjects [15, 19].

Our data reveal that the CB2 Q63 variant is associated with an earlier age of menarche in obese girls, and there was a significantly higher probability of menarche before 12 years of age.

Interestingly, the CNR2 gene, which encodes for CB2, maps to 1p36, which is a region that displays multipoint LOD scores greater than 1.00 in the Weight-Adjusted Genome Scan Analysis for Mapping Quantitative Trait Loci for Menarchal Age, although genome wide scan association studies have never identified the specific contribution of CNR2 [20]. Thus, the possibility exists that the effect exerted by the Q63R variant on the age at menarche is not detectable in women of normal weight, while it could be unmasked by facilitating environmental factors such as obesity. Consistently, by dividing patients into three different subgroups according to three different BMI ranges, we observed that the association persists only in the subgroup with the highest BMI range.

Common features of childhood obesity are elevated leptin levels and insulin resistance. Leptin, which is a hormone released by adipocytes, plays a major role in puberty initiation and progression by directly and positively acting on the hypothalamic Gonadotropin Releasing Hormone (GnRH) secretory system. Moreover, leptin receptors have been identified in both GCs and theca cells (TCs) of follicles [21]. Although we have not measured leptin levels in our cohort, leptin levels are high in obese subjects. Leptin also downregulates orexigenic factors, which include ECs [22]. Therefore, whereas the leptin signaling is uniformly high in obesity, it is possible that the CB2 receptor will better respond to the decreased EC levels when the functional variant Q63 is expressed.

Alternatively, it is plausible that the leptin resistance that accompanies obesity could mask the effect of leptin in terms of EC tone.

Insulin acts on various organs including the ovary and adipocytes to increase sex steroid bioavailability and to decrease sex hormone binding levels [23]. Thus, elevated sex steroid levels in obese prepubertal children could early activate hypothalamic-pituitary puberty. Interestingly, CB2 expression can be modulated by estrogens. In particular, 17β-estradiol increases CB2 expression, possibly by the recruitment of a putative estrogen responsive element (ERE) in the CNR2 gene [24]. Thus, the increased levels of estrogens in obese girls could enhance CB2 expression also in the follicles and could be a more efficient trigger of their maturation if the Q63 is present.

Finally, much evidence demonstrates that the EC tone is higher in obese subjects [14, 22]. Thus, the increased EC levels as well as the accompanying obesity could further stimulate the ovary through the CB2 receptor and predispose females to early menarche, whereby a more functional variant is present and more responsive.

All of these mechanisms might be enough to induce an early response of the follicles as a peripheral target, and in turn, early menarche onset.

Although CB2 receptor stimulation induces follicle growth and oocyte maturation[12], very little is known about the mechanism underlying this process.

It is noteworthy that CB2 receptor stimulation by selective agonists activates the serine/threonine protein kinase Akt signaling pathway, which is involved in cell growth, survival, proliferation and metabolism. This CB2-mediated activation prevents apoptosis of cardiac cells during ischemia/reperfusion and of central neurons after axonal damage as well as accelerates in vitro oligodendrocyte differentiation [25–28].

Several lines of evidence indicate that the Akt signaling pathway is a critical regulator of follicle growth, differentiation and survival, which determines the pool of primordial follicles and the transition from the quiescent to the growing phase, modulating the granulosa cell apoptosis throughout folliculogenesis and regulating the spindle organization in oocytes during meiosis resumption [29–32]. Moreover, suppression of Akt activity determines a delayed resumption of meiosis in mouse oocytes, and female mice deficient in the Akt1 isoform demonstrate reduced fertility and dysregulation of follicular and oocyte development, including a delay in estrus onset [33].

In conclusion, we hypothesize that the CB2 receptor may be part of the peripheral pathway that mediates the ovarian response to peripheral estrogens, possibly through the Akt pathway. In obese girls, all the permissive signals for puberty are dysregulated (such as estrogens, ECs and leptin), and the more functional CB2 receptor variant (Q63) could enhance ovarian stimulation, thus contributing to early menarche by affecting either puberty initiation or the rate of puberty progression. These data may add a new piece to the complex puzzle of female puberty regulation. Our findings could moreover help to identify obese patients who have an increased risk of developing early menarche as well as to individualize treatment and follow-up.

Acknowledgments

The authors would like to thank all of the girls and their parents for their participation in the study.

Author Contributions

Conceived and designed the experiments: GB AG FR. Performed the experiments: GB AG MT. Analyzed the data: GB AG EMdG FR. Contributed reagents/materials/analysis tools: FR BN. Wrote the paper: GB AG FR. Supervised manuscript writing: EMdG BN SM FR. Supervised experiments: FR. Supervised clinical data: LP. Supervised molecular analysis: SM. Supervised statistic analysis: EMdG.

Citation: Bellini G, Grandone A, Torella M, Miraglia del Giudice E, Nobili B, Perrone L, et al. (2015) The Cannabinoid Receptor 2 Q63R Variant Modulates the Relationship between Childhood Obesity and Age at Menarche. PLoS ONE 10(10): e0140142. https://doi.org/10.1371/journal.pone.0140142

References

1. Banerjee I, Clayton P. The genetic basis for the timing of human puberty. J Neuroendocrinol. 2007;19: 831–838. pmid:17927661

2. Dreyfus J, Jacobs DR Jr, Mueller N, Schreiner PJ, Moran A, Carnethon MR, et al. Age at Menarche and Cardiometabolic Risk in Adulthood: The Coronary Artery Risk Development in Young Adults Study. J Pediatr. 2015. pii:S0022-3476(15)00393-5.

3. Perry JR, Day F, Elks CE, Sulem P, Thompson DJ, Ferreira T, et al. Parent-of-origin-specific allelic associations among 106 genomic loci for age at menarche. Nature. 2014;514: 92–97. pmid:25231870

4. Ong KK, Elks CE, Li S, Zhao JH, Luan J, Andersen LB, et al. Genetic variation in LIN28B is associated with the timing of puberty. Nat Genet. 2009;41: 729–733. pmid:19448623

5. Grandone A, Cirillo G, Messa F, Coppola R, Festa A, Perrone L, et al. LIN28B polymorphism could modulate the relationship between childhood obesity and age at menarche. J Adolesc Health. 2013;52: 375. pmid:23427786

6. Ahmed ML, Ong KK, Dunger DB. Childhood obesity and the timing of puberty. Trends Endocrinol Metab. 2009;20: 237–242. pmid:19541497

7. Chou SH, Mantzoros C. 20 years of leptin: role of leptin in human reproductive disorders. J Endocrinol. 2014;223: T49–T62. pmid:25056118

8. Condorelli RA, Calogero AE, Vicari E, Mongioi' L, Favilla V, Morgia G, et al. The gonadal function in obese adolescents: review. J Endocrinol Invest. 2014 [Epub ahead of print].

9. Wagner IV, Sabin MA, Pfäffle RW, Hiemisch A, Sergeyev E, Körner A, et al. Effects of obesity on human sexual development. Nat Rev Endocrinol. 2012;8: 246–254. pmid:22290357

10. Wan W, Deng X, Archer KJ, Sun SS. Pubertal pathways and the relationship to anthropometric changes in childhood: The Fels longitudinal study. Open J Pediatr. 2012;2. pmid:24307981

11. Maccarrone M, Bab I, Bíró T, Cabral GA, Dey SK, Di Marzo V, et al. Endocannabinoid signaling at the periphery: 50 years after THC. Trends Pharmacol Sci. 2015;36: 277–296. pmid:25796370

12. El-Talatini MR, Taylor AH, Elson JC, Brown L, Davidson AC, Konje JC. Localisation and function of the endocannabinoid system in the human ovary. PLoSOne. 2009;4: e4579.

13. Cecconi S, Rossi G, Castellucci A, D'Andrea G, Maccarrone M. Endocannabinoid signaling in mammalian ovary. Eur J Obstet Gynecol Reprod Biol. 2014;178: 6–11. pmid:24948047

14. Matias I, Petrosino S, Racioppi A, Capasso R, Izzo AA, Di Marzo V. Dysregulation of peripheral endocannabinoid levels in hyperglycemia and obesity: Effect of high fat diets. Mol Cell Endocrinol. 2008;286: S66–78. pmid:18343566

15. Sipe JC, Arbour N, Gerber A, Beutler E. Reduced endocannabinoid immunemodulation by a common cannabinoid 2 (CB2) receptor gene polymorphism: possible risk for autoimmune disorders. J Leukoc Biol. 2005;78: 231–238. pmid:15845647

16. Rossi F, Mancusi S, Bellini G, Roberti D, Punzo F, Vetrella S, et al. CNR2 functional variant (Q63R) influences childhood immune thrombocytopenic purpura. Haematologica. 2011;96: 1883–1885. pmid:21828121

17. Rossi F, Bellini G, Tolone C, Luongo L, Mancusi S, Papparella A, et al. The cannabinoid receptor type 2 Q63R variant increases the risk of celiac disease: implication for a novel molecular biomarker and future therapeutic intervention. Pharmacol Res. 2012;66: 88–94. pmid:22465144

18. Bellini G, Olivieri AN, Grandone A, Alessio M, Gicchino MF, Nobili B, et al. Association between cannabinoid receptor type 2 Q63R variant and oligo/polyarticular juvenile idiopathic arthritis. Scand J Rheumatol. 2015;14:1–4.

19. Carrasquer A, Nebane NM, Williams WM, Song ZH. Functional consequences of nonsynonymous single nucleotide polymorphisms in the CB2 cannabinoid receptor. Pharmacogenet Genomics. 2010;20: 157–166. pmid:20124950

20. Rothenbuhler A, Fradin D, Heath S, Lefevre H, Bouvattier C, Lathrop M, et al. Weight-adjusted genome scan analysis for mapping quantitative trait Loci for menarchal age. J ClinEndocrinolMetab. 2006;91: 3534–3537.

21. Abir R, Ao A, Jin S, Barnett M, Raanani H, Ben-Haroush A, et al. Leptin and its receptors in human fetal and adult ovaries. Fertil Steril. 2005; 84: 1779–1782. pmid:16359991

22. Di Marzo V, Goparaju SK, Wang L, Liu J, Bátkai S, Járai Z, et al. Leptin-regulated endocannabinoids are involved in maintaining food intake. Nature.2001;410: 822–825. pmid:11298451

23. Pinkney J, Streeter A, Hosking J, Mostazir M, Jeffery A, Wilkin T. Adiposity, chronic inflammation, and the prepubertal decline of sex hormone binding globulin in children: evidence for associations with the timing of puberty. J Clin Endocrinol Metab. 2014;99: 3224–3232. Erratum in: J Clin Endocrinol Metab. 2015;100: 763. pmid:24926948

24. Rossi F, Bellini G, Luongo L, Mancusi S, Torella M, Tortora C, et al. The 17-β-oestradiol inhibits osteoclast activity by increasing the cannabinoid CB2 receptor expression. Pharmacol Res. 2013;68: 7–15. pmid:23142558

25. Li Q, Wang F, Zhang YM, Zhou JJ, Zhang Y. Activation of cannabinoid type 2 receptor by JWH133 protects heart against ischemia/reperfusion-induced apoptosis. Cell Physiol Biochem. 2013;31: 693–702. pmid:23711495

26. Viscomi MT, Oddi S, Latini L, Pasquariello N, Florenzano F, Bernardi G, et al. Selective CB2 receptor agonism protects central neurons from remote axotomy-induced apoptosis through the PI3K/Akt pathway. J Neurosci. 2009;29: 4564–4570. pmid:19357281

27. Gomez O, Sanchez-Rodriguez A, Le M, Sanchez-Caro C, Molina-Holgado F, Molina-Holgado E. Cannabinoid receptor agonists modulate oligodendrocyte differentiation by activating PI3K/Akt and the mammalian target of rapamycin (mTOR) pathways. Br J Pharmacol. 2011;163: 1520–1532. pmid:21480865

28. Gomez O, Sanchez-Rodriguez MA, Ortega-Gutierrez S, Vazquez-Villa H, Guaza C, Molina-Holgado F, et al. A Basal Tone of 2-Arachidonoylglycerol Contributes to Early Oligodendrocyte Progenitor Proliferation by Activating Phosphatidylinositol 3-Kinase (PI3K)/AKT and the Mammalian Target of Rapamycin (MTOR) Pathways. J Neuroimmune Pharmacol. 2015;10: 309–317. pmid:25900077

29. Cecconi S, Mauro A, Cellini V, Patacchiola F. The role of Akt signalling in the mammalian ovary. Int J Dev Biol. 2012;56: 809–817. pmid:23417403

30. Hoshino Y, Sato E. Protein kinase B (PKB/Akt) is required for the completion of meiosis in mouse oocytes. Dev Biol. 2008;314: 215–223. pmid:18177853

31. Kimura T, Tomooka M, Yamano N, Murayama K, Matoba S, Umehara H, et al. AKT signaling promotes derivation of embryonic germ cells from primordial germ cells. Development. 2008;135: 869–879. pmid:18216173

32. Tomek W, Smiljakovic T. Activation of Akt (protein kinase B) stimulates metaphase I to metaphase II transition in bovine oocytes. Reproduction. 2005;130: 423–430. pmid:16183860

33. Kalous J, Kubelka M, Solc P, Susor A, Motlik J. AKT (protein kinase B) is implicated in meiotic maturation of porcine oocytes. Reproduction. 2009;138: 645–654. pmid:19633130