A variant in the 3ʹ-untranslated region of the MC2R gene decreases the risk of schizophrenia in a female Han Chinese population

Abstract

Objective

Schizophrenia is a complex mental disorder with high heritability. The hypothalamic–pituitary–adrenal (HPA) axis, which is the stress system of the neuroendocrine system, is considered to impact psychotic disorders. We hypothesized that polymorphisms of HPA axis genes might be involved in the development of schizophrenia.

Methods

A case–control study comprising 234 patients with schizophrenia and 399 matched healthy controls was conducted to investigate the association between the human melanocortin 2 receptor (MC2R) gene and schizophrenia risk. Seven tag single nucleotide polymorphisms (SNPs) (rs16941303, rs16941314, rs2186944, rs28926188, rs7230126, rs948322, and rs948331) of MC2R were genotyped by direct sequencing.

Results

No significant associations were observed between any of the alleles, genotypes, or haplotypes examined within the MC2R gene and the risk of schizophrenia in the total group or in subgroups stratified by smoking or alcoholism. However, a subgroup analysis stratified by sex revealed that under the additive model, the C allele of the MC2R rs948331 SNP significantly decreased the risk of schizophrenia in females (odds ratio=0.18).

Conclusion

The C allele of the MC2R rs948331 locus may be a protective factor, reducing the risk of schizophrenia in the female Han Chinese population.

Keywords

Melanocortin 2 receptor schizophrenia polymorphism Chinese Han female hypothalamic–pituitary–adrenal axis

Introduction

Schizophrenia, a complex disorder typified by delusions, hallucinations and thought disorders, has both genetic and environmental triggers.¹ Stress sensitivity and activity of the hypothalamic–pituitary–adrenocortical (HPA) axis have been implicated in the development of schizophrenia.² When exposed to stress, the HPA axis is activated by binding of corticotropin-releasing hormone (CRH) receptors, which stimulates adrenocorticotrophic hormone (ACTH) release³ and aids in inhibiting CRH expression in amygdala neurons via negative feedback. ACTH can induce the formation of myelin sheaths and dendrites by increasing the circulating sugar content, thus promoting brain maturation and shortening the vulnerable period. It can also act as a neuroregulator and neurotransmitter that can affect nervous system development, whereby it plays roles in preventing convulsion, improving learning and memory functions, and reducing attention deficits.^4–6 ACTH can have pharmacological effects by binding to the melanocortin 2 receptor (MC2R, also known as the ACTH receptor).⁷ Recently, five MCR genes (MC1R, MC2R, MC3R, MC4R, and MC5R) were identified. MC2R is an important member of this family that is mainly distributed in the adrenal cortex and binds ACTH as a ligand.⁸

Thus, it is plausible that variants of this receptor gene may be involved in disease pathogenesis and drug resistance. Single nucleotide polymorphisms (SNPs) in the coding region of the MC2R gene could influence the molecular conformation of the ATCH receptor. Such structural differences could influence the efficiency of ligand binding and the strength of the transduction of signals sent through intracellular pathways, leading to diseases such as infantile spasm,⁹ alopecia areata,¹⁰ heroin addiction,¹¹ and familial glucocorticoid deficiency.^12,13

Few studies have examined the link between MC2R gene polymorphisms and schizophrenia susceptibility. Using a case–control method, we investigated MC2R gene polymorphisms to determine whether they contributed to the development of schizophrenia in a Chinese population.

Materials and methods

Patients

Patients with schizophrenia were recruited from the First Affiliated Hospital of Changsha Medical College; their diagnoses were confirmed by at least two independent experienced psychiatrists according to the criteria in the Diagnostic and Statistical Manual of Mental Disorders, 4th Edition, Text Revision. Individuals without schizophrenia were recruited for the control group. This study followed the ethical standards of the 1964 Helsinki Declaration and its later amendments, as well as national and international guidelines. The research procedures were approved by the Ethics Committee of the Changsha Medical College, China (Ethical Committee of Changsha Medical College (CMCEC-1801022)). Written informed consent for genetic analyses was obtained from all subjects.

Genotyping

SNP selection was performed using Haploview software (https://sourceforge.net/projects/haploview/) choosing only SNPs with a minor allele frequency higher than 0.05 and r²≥0.8 based on the HapMap database (CHB, Chinese Han population) (http://coriell.org/1/NHGRI/Collections/HapMap-Collections/HapMap-Project). After screening, seven tag SNPs (rs16941303, rs16941314, rs2186944, rs28926188, rs7230126, rs948322, and rs948331) in the MC2R gene were selected. Genomic DNA was extracted from peripheral blood leukocytes using the standard phenol-chloroform method. Multiplex polymerase chain reaction (PCR) assays were performed on a GeneAmp 9700 PCR thermocycler (Applied Biosystems, Foster City, California, USA). The reactions were performed in a total volume of 10 µL, including 1 µL of genomic DNA, 1 µL of primer mix (1 µM of each primer), and 1 U of Hotstar Taq polymerase (Qiagen Inc., Valencia, CA, USA). The cycling parameters were as follows: 95°C for 2 minutes; 10 cycles at 95°C for 20 s, 65°C for 40 s, and 72° for 30 s; 25 cycles at 95°C for 20 s, 55°C for 30 s, and 72°C for 1 minute; then the samples were held at 4°C. The genotypes of the MC2R polymorphisms were determined by Sanger sequencing. A 3730XL genetic analyzer (Applied Biosystems) was used for sequencing. GeneScan ™3.7 software (Applied Biosystems) was used for data analysis.

Statistical analyses

The Hardy–Weinberg equilibrium (HWE) of all polymorphisms under study was assessed using a χ² test. The distributions of the genetic models (allelic, dominant, additive, and recessive models) were compared between the control and case groups using a t test. All calculations were performed with Plink 1.0 software (http://softwaretopic.informer.com/plink-software/). Corresponding odds ratios (ORs) and 95% confidence intervals (CIs) were evaluated by logistic regression analysis. A stringent Bonferroni correction was performed to adjust the p-value (p_adj) in all cases in which an initial significant relationship was found. Squared correlation coefficients (r²) for each pair of polymorphic sites were calculated with Haploview software. A value of p<0.05 indicated statistical significance. The allele data from patients stratified into subgroups according to alcoholism and smoking were also analyzed. A power analysis was conducted using G*Power software (www.gpower.hhu.de).

Results

A total of 234 patients with schizophrenia and 399 healthy controls were included in this study. The case group consisted of 234 unrelated patients including 130 males and 104 females (age 33.1±7.2 years, mean±standard deviation (SD)). Within the case group, 7.4% of the patients had a history of alcoholism, and 35.0% had a history of smoking. In addition, 399 unrelated healthy individuals (302 males and 97 females; age 27.25±8.4 years (mean±SD)) who did not have schizophrenia according to clinical examinations performed at The First Affiliated Hospital of Changsha Medical College were recruited for the control group. Within the control group, 62.2% of the participants had a history of alcoholism, and 46.8% had a history of smoking (Table 1). No significant deviations from HWE were observed for any of the seven tag SNPs (Table 2). Linkage disequilibrium analysis indicated that rs16941303 and rs2186944, rs16941303 and rs28926188, and rs2186944 and rs28926188 exhibited strong linkage (r²>0.85) (Figure 1). Moreover, our study obtained >71.4% power at the 5% significance level (two-tailed), which indicates that the sample size was sufficiently large to detect an association with an OR of 1.5 (Table 2).

Table 1.

Characteristics of the case and control groups.

	Case (n=234)	Control (n=399)
Age, mean±SD (years)	33.1±7.2	27.25±8.4
Sex, male:female (%)	55.6:44.4	75.7:24.3
Smoking, n	82	187
Alcoholism, n	24	248
BMI, mean±SD (kg/m²)	22.48±6.7	21.87±5.4
Age of onset, mean±SD (years)	21.13±3.3	21.47±6.4

SD, standard deviation; BMI, body mass index.

Table 2.

Allele distributions of MC2R gene polymorphisms in schizophrenia patients (case) and healthy controls (control).

SNP ID	Position	SNP property	Minor allele (A<B)	Minor allele		p	OR [95% CI]^a	Power^b	HWE
SNP ID	Position	SNP property	Minor allele (A<B)	Case (n=234)	Control (n=399)	p	OR [95% CI]^a	Power^b	HWE
rs16941303	18:13912243	intron 2	G<A	85 (0.182)	135 (0.169)	0.572	0.917 [0.681–1.237]	71.4	0.881
rs16941314	18:13914355	intron 2	T<C	60 (0.128)	109 (0.137)	0.672	1.075 [0.767–1.508]		0.850
rs2186944	18:13915542	5ʹ-flanking	A<G	85 (0.182)	136 (0.170)	0.614	1.080 [0.801–1.456]		0.834
rs28926188	18:13882728	3ʹ-UTR_exon2	A<G	84 (0.179)	148 (0.185)	0.791	0.961 [0.714–1.291]		0.672
rs7230126	18:13890199	intron 2	C<T	82 (0.175)	146 (0.183)	0.729	0.948 [0.704–1.278]		0.905
rs948322	18:13886846	intron 2	A<G	91 (0.194)	180 (0.226)	0.192	0.828 [0.624–1.099]		0.841
rs948331	18:13881279	3ʹ-flanking	C<T	178 (0.380)	322 (0.404)	0.416	0.907 [0.717–1.146]		0.528

^aThe OR and 95% CI are calculated for the minor allele of each polymorphism.

^bPower calculations assume α=0.05 and a small effect size (0.1) or OR=1.5.

A, major allele; B, minor allele; SNP, single nucleotide polymorphism; OR, odds ratio; 95% CI, 95% confidence interval; HWE, Hardy–Weinberg equilibrium.

Figure 1.

Linkage disequilibrium analysis of seven single nucleotide polymorphisms of MC2R in the case and control groups.

No association was detected between the allelic, additive, dominant, and recessive models of any of the seven SNPs in the MC2R gene and schizophrenia in the Han Chinese population before or after Bonferroni correction (Table 2 and 3). Seven major haplotypes were identified among the rs16941303, rs16941314, rs2186944, rs28926188, rs7230126, rs948322, and rs948331 loci (lowest frequency threshold >0.01). Haplotype_ACGGTGT was the most prevalent in both the case and control groups, and no significant difference in its frequency was observed between patients with schizophrenia and controls. The frequency of haplotype_ACGATAC was significantly increased in the control group versus the case group (p=0.039, OR [95% CI]=0.153 [0.020–1.191]), but significance was not obtained after correction (p=0.273) (Table 4).

Table 3.

Genetic model distributions of MC2R gene polymorphisms in the case and control groups.

SNP ID (A<B)	Case			Control			Additive model (AA vs. BB)		Dominant model (AA+AB vs. BB)		Recessive model (AA vs. AB+BB)
SNP ID (A<B)	AA	AB	BB	AA	AB	BB	p^a	OR [95% CI]^b	p^a	OR [95% CI]^b	p^a	OR [95% CI]^b
rs16941303 (G<A)	5	75	154	11	113	275	0.69	0.81 [0.27–2.36]	0.42	1.15 [0.82–1.62]	0.63	0.77 [0.26–2.24]
rs16941314 (T<C)	4	52	178	7	95	297	0.94	0.95 [0.27–3.30]	0.65	0.92 [0.63–1.33]	0.97	0.97 [0.28–3.36]
rs2186944 (A<G)	5	75	154	11	114	274	0.69	0.80 [0.27–2.35]	0.46	1.14 [0.81–1.60]	0.63	0.77 [0.26–2.24]
rs28926188 (A<G)	5	74	155	15	118	266	0.27	0.56 [0.20–1.57]	0.91	1.02 [0.72–1.43]	0.27	0.56 [0.20–1.56]
rs7230126 (C<T)	6	70	158	13	120	266	0.60	0.77 [0.29–2.06]	0.83	0.96 [0.68–1.36]	0.78	0.78 [0.29–2.08]
rs948322 (A<G)	8	75	151	21	138	240	0.21	0.58 [0.25–1.35]	0.27	0.83 [0.59–1.16]	0.29	0.64 [0.28–1.46]
rs948331 (C<T)	35	108	91	68	186	145	0.20	0.71 [0.42–1.21]	0.52	0.90 [0.64–1.25]	0.49	0.86 [0.55–1.33]

^ap values were calculated using Fisher’s exact test.

^bThe OR [95% CI] was estimated by logistic regression analysis.

A, major allele; B, minor allele; SNP, single nucleotide polymorphism; OR, odds ratio; 95% CI, 95% confidence interval.

Table 4.

Frequencies of haplotypes (>1%) of the seven single nucleotide polymorphisms of MC2R in patients with schizophrenia (case) and healthy controls (control).

Haplotype^a	Case	Control	p^b	p_adj^c	OR [95% CI]
A C G A T A C	1 (0.002)	11 (0.014)	0.039	0.273	0.153 [0.020–1.191]
A C G G C A C	82 (0.175)	144 (0.180)	0.813	–	0.965 [0.715–1.301]
A C G G T G C	12 (0.025)	30 (0.038)	0.224	–	0.656 [0.331–1.301]
A C G G T G T	231 (0.493)	362 (0.454)	0.176	–	1.172 [0.931–1.475]
A T G G T G T	57 (0.122)	107 (0.134)	0.545	–	0.899 [0.638–1.268]
G C A A T A C	5 (0.011)	18 (0.023)	0.115	–	0.458 [0.169–1.240]
G C A A T G C	75 (0.161)	117 (0.146)	0.485	–	1.119 [0.816–1.534]

^aThe Plink program was used to estimate common (frequency >0.01) haplotypes of the seven SNPs of MC2R (rs16941303, rs16941314, rs2186944, rs28926188, rs7230126, rs948322, and rs948331).

^bEach haplotype was compared with all other haplotypes combined.

^cThe Bonferroni correction was applied to correct the p value.

SNP, single nucleotide polymorphism; OR, odds ratio; 95% CI, 95% confidence interval.

Analyses were performed on subgroups stratified by sex, smoking, and alcoholism. As shown in Table 5, the allele distributions under the allelic, additive, dominant, and recessive models for the MC2R rs16941303, rs16941314, rs2186944, rs28926188, rs7230126, and rs948322 loci did not significantly differ between the case and control groups. Interestingly, significant associations were observed under the allelic (p=0.046, OR [95% CI]=0.665 [0.445–0.992]) and additive models (p=0.0008, OR [95% CI]=0.18 [0.06–0.49]) for rs948331 and schizophrenia in female participants but not in male participants or in the smoking and alcoholism groups. The association under the allelic model disappeared after Bonferroni correction (p=0.362). A significant association remained under the additive model even after Bonferroni correction (p=0.0056).

Table 5.

The association between genetic models of MC2R single nucleotide polymorphisms and schizophrenia stratified by sex, smoking, and alcoholism.

SNPs	Subgroup (case/control)	Allelic model (A vs. B)		Additive model (AA vs .BB)		Dominant model (AA+AB vs. BB)		Recessive model (AA vs. AB+BB)
SNPs	Subgroup (case/control)	p^a	OR [95% CI]	p^a	OR [95% CI]	p^a	OR [95% CI]	p^a	OR [95% CI]
rs16941303	Female (104/97)	0.649	1.126 [0.674–1.881]	0.76	1.33 [0.21–8.20]	0.52	0.82 [0.46–1.48]	0.71	1.41 [0.23–8.63]
	Male (130/302)	0.311	0.823 [0.565–1.199]	0.47	0.56 [0.12–2.67]	0.15	1.38 [0.89–2.13]	0.39	0.51 [0.11–2.39]
	Smoking (82/151)	0.843	1.051 [0.641–1.724]	0.22	0.16 [0.01–2.96]	0.81	1.07 [0.61–1.89]	0.22	0.16 [0.01–2.96]
	Alcoholism (24/210)	0.911	0.957 [0.445–2.059]	0.87	0.79 [0.04–14.96]	0.72	1.17 [0.49–2.82]	0.86	0.76 [0.04–14.21]
rs16941314	Female (104/97)	0.968	0.988 [0.546–1.787]	0.38	2.81 [0.29–27.60]	0.79	0.91 [0.47–1.76]	0.37	2.85 [0.29–27.89]
	Male (130/302)	0.696	1.108 [0.710–1.668]	0.37	0.38 [0.04–3.17]	0.87	0.96 [0.60–1.54]	0.38	0.38 [0.05–3.21]
	Smoking (82/151)	0.567	0.849 [0.485–1.485]	0.71	0.65 [0.07–6.38]	0.46	1.26 [0.68–2.35]	0.67	0.61 [0.06–5.95]
	Alcoholism (24/210)	0.599	1.295 [0.491–3.413]	0.94	0.89 [0.05–17.09]	0.71	0.82 [0.29–2.31]	0.97	0.94 [0.05–17.92]
rs2186944	Female (104/97)	0.649	0.887 [0.531–1.482]	0.76	1.33 [0.21–8.20]	0.52	0.82 [0.46–1.48]	0.71	1.41 [0.23–8.63]
	Male (130/302)	0.341	1.200 [0.824–1.747]	0.47	0.56 [0.12–2.65]	0.17	1.36 [0.88–2.09]	0.39	0.51 [0.11–2.39]
	Smoking (82/151)	0.843	0.951 [0.579–1.560]	0.22	0.16 [0.01–2.96]	0.81	1.07 [0.61–1.89]	0.22	0.16 [0.01–2.96]
	Alcoholism (24/210)	0.911	1.044 [0.485–2.247]	0.87	0.79 [0.04–14.96]	0.72	1.17 [0.49–2.82]	0.86	0.76 [0.04–14.21]
rs28926188	Female (104/97)	0.475	0.827 [0.492–1.391]	0.78	1.29 [0.21–7.97]	0.35	0.75 [0.41–1.36]	0.71	1.41 [0.23–8.63]
	Male (130/302)	0.712	1.072 [0.741–1.549]	0.20	0.37 [0.08–1.68]	0.32	1.24 [0.81–1.90]	0.17	0.35 [0.08–1.56]
	Smoking (82/151)	0.912	0.972 [0.592–1.597]	0.23	0.17 [0.01–3.05]	0.73	1.11 [0.63–1.95]	0.22	0.16 [0.01–2.96]
	Alcoholism (24/210)	0.878	1.061 [0.493–2.284]	0.88	0.79 [0.04–15.07]	0.68	1.20 [0.50–2.88]	0.86	0.76 [0.04–14.21]
rs7230126	Female (104/97)	0.397	0.807 [0.492–1.324]	0.15	0.36 [0.09–1.45]	0.72	0.90 [0.50–1.61]	0.17	0.38 [0.10–1.52]
	Male (130/302)	0.978	0.995 [0.677–1.460]	0.84	1.16 [0.28–4.74]	0.92	0.98 [0.63–1.52]	0.83	1.17 [0.29–4.73]
	Smoking (82/151)	0.110	1.483 [0.912–2.412]	0.34	2.23 [0.43–11.48]	0.13	1.56 [0.88–2.74]	0.45	1.87 [0.37–9.50]
	Alcoholism (24/210)	0.524	1.272 [0.606–2.669]	0.84	0.73 [0.04–13.72]	0.31	1.56 [0.66–3.69]	0.77	0.64 [0.04–11.75]
rs948322	Female (104/97)	0.812	0.395 [0.501–1.313]	0.17	0.38 [0.09–1.54]	0.71	0.90 [0.50–1.59]	0.17	0.38 [0.10–1.52]
	Male (130/302)	0.342	0.839 [0.585–1.204]	0.59	0.75 [0.26–2.16]	0.33	0.81 [0.53–1.24]	0.71	0.82 [0.29–2.33]
	Smoking (82/151	0.314	1.273 [0.794–2.039]	0.31	2.11 [0.50–8.81]	0.42	1.25 [0.72–2.19]	0.38	1.88 [0.46–7.74]
	Alcoholism (24/210)	0.798	1.263 [0.617–2.584]	0.75	1.43 [0.16–12.53]	0.50	1.34 [0.57–3.17]	0.83	1.26 [0.15–10.71]
rs948331	Female (104/97)	0.046	0.665 [0.445–0.992]	0.0008	0.18 [0.06–0.49]	0.17	0.67 [0.37–1.19]	0.06	0.49 [0.23–1.03]
	Male (130/302)	0.728	1.054 [0.783–1.419]	0.52	1.25 [0.63–2.47]	0.87	1.04 [0.68–1.59]	0.66	1.13 [0.65–1.98]
	Smoking (82/151	0.128	1.251 [0.916–1.993]	0.03	2.75 [1.09–6.94]	0.19	1.46 [0.83–2.56]	0.30	1.47 [0.71–3.05]
	Alcoholism (24/210)	0.389	1.303 [0.712∼2.382]	0.16	2.94 [0.66–13.19]	0.56	1.31 [0.54–3.19]	0.40	1.58 [0.55–4.55]

^ap values were calculated using Fisher’s exact test.

^bThe OR [95% CI] was estimated by logistic regression analysis.

A, major allele; B, minor allele; SNP, single nucleotide polymorphism; OR, odds ratio; 95% CI, 95% confidence interval.

The haplotypes were also analyzed based on sex, smoking, and alcoholism (Table 6), and the frequency of haplotype_ACGGTGT in patients with schizophrenia was significantly lower than that in controls among the female subgroup (p=0.041, OR [95% CI]=1.508 [1.015–2.240]). The frequency of haplotype_ACGGTGT in patients with schizophrenia was also significantly lower than that in the controls in the smoking subgroup (p=0.031, OR [95% CI]=0.659 [0.448–0.968]). However, these significant associations disappeared after Bonferroni correction (females: p=0.287; smoking: 0.217).

Table 6.

The association between haplotypes (>1%) of the seven single nucleotide polymorphisms of MC2R and schizophrenia risk stratified by sex, smoking, and alcoholism.

Haplotype^a	Subgroups	Case	Control	p^b	p_adj^c	OR [95% CI]
A C G A T A C	Female (104/97)	–	–	–	–	–
	Male (130/302)	1.00 (0.004)	11.00 (0.018)	0.136		0.206 [0.037–1.135]
	Smoking (82/151)	–	–	–	–	–
	Alcoholism (24/210)	–	–	–	–	–
A C G G C A C	Female (104/97)	36.99 (0.178)	41.00 (0.211)	0.393	–	0.806 [0.491–1.323]
	Male (130/302)	45.00 (0.173)	102.96 (0.170)	0.981	–	1.004 [0.683–1.476]
	Smoking (82/151	35.00 (0.213)	46.95 (0.154)	0.127	–	1.458 [0.896–2.373]
	Alcoholism (24/210)	10.00 (0.208)	71.97 (0.171)	0.557	–	1.248 [0.595–2.620]
A C G G T G C	Female (104/97)	4.00 (0.019)	10.01 (0.052)	0.077	–	0.360 [0.111–1.168]
	Male (130/302)	7.78 (0.030)	19.21 (0.032)	0.863	–	0.929 [0.399–2.164]
	Smoking (82/151)	6.00 (0.037)	5.33 (0.018)	0.211	–	2.094 [0.641–6.833]
	Alcoholism (24/210)	1.23 (0.026)	9.46 (0.023)	0.906	–	1.121 [0.167–7.526]
A C G G T G T	Female (104/97)	109.02 (0.524)	82.00 (0.423)	0.041	0.287	1.508 [1.015–2.240]
	Male (130/302)	122.01 (0.469)	281.48 (0.466)	0.954	–	0.991 [0.740–1.328]
	Smoking (82/151)	71.01 (0.433)	160.43 (0.528)	0.031	0.217	0.659 [0.448–0.968]
	Alcoholism (24/210)	22.78 (0.475)	209.22 (0.498)	0.677	–	0.881 [0.484–1.602]
A T G G T G T	Female (104/97)	22.98 (0.110)	23.00 (0.119)	0.798	–	0.923 [0.499–1.706]
	Male (130/302)	33.99 (0.131)	82.97 (0.137)	0.747	–	0.932 [0.607–1.431]
	Smoking (82/151)	22.99 (0.140)	33.56 (0.110)	0.378	–	1.291 [0.731–2.280]
	Alcoholism (24/210)	4.22 (0.088)	51.78 (0.123)	0.454	–	0.673 [0.238–1.906]
G C A A T A C	Female (104/97)	–	–	–	–	–
	Male (130/302)	4.78 (0.018)	16.66 (0.028)	0.412	–	0.652 [0.233–1.823]
	Smoking (82/151)	–	–	–	–	–
	Alcoholism (24/210)	1.00 (0.021)	3.26 (0.008)	0.178	–	2.673 [0.393–18.168]
G C A A T G C	Female (104/97)	29.99 (0.144)	34.98 (0.180)	0.323	–	0.765 [0.449–1.303]
	Male (130/302)	45.22 (0.174)	82.34 (0.136)	0.173	–	1.316 [0.886–1.956]
	Smoking (82/151)	28.00 (0.171)	47.67 (0.157)	0.751	–	1.086 [0.651–1.811]
	Alcoholism (24/210)	7.99 (0.167)	67.54 (0.161)	0.955	–	1.023 [0.458–2.283]

^aThe Plink program was used to estimate common (frequency >0.01) haplotypes of the seven SNPs of MC2R (rs16941303, rs16941314, rs2186944, rs28926188, rs7230126, rs948322, and rs948331).

^bEach haplotype was compared with all other haplotypes combined.

^cThe Bonferroni correction was applied to correct the p value.

SNP, single nucleotide polymorphism; OR, odds ratio; 95% CI, 95% confidence interval.

Discussion

The HPA axis is an important aspect of the neuroendocrine system.¹⁴ Under normal physiological conditions, the final products released by the HPA axis, corticosteroids (such as cortisol), regulate physiological and psychological stress through complex feedback mechanisms, helping the body maintain homeostasis.^15,16 Continuous stress induces hyperactivity of the HPA axis, increasing corticosteroid secretion, and eventually disrupting the HPA axis function.¹⁷ In the HPA axis, ACTH can promote adrenal glucocorticoid secretion by binding to ACTH receptor 1, which then increases plasma cortisol levels.¹⁸ In contrast, there are inconsistencies among the reported results with regard to plasma ACTH levels in patients with schizophrenia. Duval et al. did not find any differences between their patient and control groups,¹⁹ while Jiang et al.²⁰ reported that the plasma ACTH levels of patients with schizophrenia were higher than the normal level. These contradictory results may be caused by the different methods of ACTH determination used in these respective studies.

The MC2R gene maps to chromosome l8p l1.2.²¹ Missense mutations and variants in the MC2R gene may result in the loss of signal transduction following ligand binding or may affect the affinity of ligand binding, but both types of defects would result in impaired receptor function.²² Further studies are needed to investigate functional variants that are either regulated naturally or alternatively spliced upon ACTH release. Interestingly, one report described an association between polymorphisms of the MC2R gene in relation to other psychotic disorders including infantile spasm 9 and heroin addiction.¹¹ In this case, the heterogeneity of different diseases might explain the differences between their data and those of the present study. Follow-up studies with a larger population and other MC2R polymorphisms are required to determine the complex relationship of genetic and environmental factors in the development of schizophrenia.

Interestingly, we detected a significant association between the additive model of MC2R rs948331 and schizophrenia in a female Han Chinese population. To our knowledge, this is the first time that a significant relationship has been identified between any variant in the MC2R gene and schizophrenia in the Han Chinese population. This SNP is located in the 3ʹ-untranslated region of MC2R, but its physiological effect is unknown. We hypothesize that this polymorphism may decrease MC2R expression, resulting in a reduction of plasma cortisol, which would then increase the synthesis and secretion of CRH via a negative feedback loop to protect against schizophrenia. Based on the knowledge that autoimmune disorders can be triggered in women following menarche or childbirth, resulting in increased risk of a myriad of autoimmune diseases, it is likely that other factors that affect the autoimmune axis and the HPA axis may indirectly affect the cognitive functions affected by these axes.²³ Most neuroendocrine studies have focused on pharmacology, disease course, and symptoms without considering sex differences.²⁴ Our results suggest that sex might be an important factor in the genetic association between the MC2R rs948331 locus alleles and the pathogenesis of schizophrenia.

Several limitations of this study must be considered. First, the positive association found in this study may be a result of type I error; however, this limitation was addressed using the Bonferroni correction. Second, the sample size of this study is relatively small; however, a power analysis showed a 71.4% probability of detecting an effect size of 0.5 at the 0.05 significance level. Therefore, based on these caveats, we suggest that the association between the MC2R rs948331 locus and susceptibility to schizophrenia should be replicated in larger studies and among additional ethnic populations. Third, the definitions of types of alcoholism and smoking (chronic and acute) among patients in the study were not precise, which could have influenced the consistency of our results.

Conclusions

Our study showed for the first time that under the additive model, the C allele of the MC2R rs948331 locus might be a protective factor against schizophrenia among the female Chinese Han population.

Supplemental Material

sj-pdf-1-imr-10.1177_03000605211029504 - Supplemental material for A variant in the 3ʹ-untranslated region of the MC2R gene decreases the risk of schizophrenia in a female Han Chinese population

Supplemental material, sj-pdf-1-imr-10.1177_03000605211029504 for A variant in the 3ʹ-untranslated region of the MC2R gene decreases the risk of schizophrenia in a female Han Chinese population by Liang Tang, Qin Xiang, Ju Xiang, Jianming Li and Danna Chen in Journal of International Medical Research

Footnotes

Acknowledgment

We are grateful to all of the individuals who participated in this study.

Declaration of conflicting interest

The authors declare that there is no conflict of interest.

Funding

The present study was funded by the National Natural Science Foundation of China (Grant No. 81873780, 61702054); Hunan Natural Science Foundation Youth Program (2019JJ50697, 2018JJ3568, 2018JJ3573); Changsha Outstanding Innovative Young People Training Scheme (kq2009095, kq2009093, kq2004077, kq2009092); Foundation of the Education Department of Hunan Province (19A058, 19B072, 18A493); Foundation of the Health and Family Planning Commission of Hunan Province (20201918, 20201910); Application Characteristic Discipline of Hunan Province; Hunan Key Laboratory Cultivation Base of the Research and Development of Novel Pharmaceutical Preparations (No. 2016TP1029); and Project of Changsha Science and Technology (kq2004077).

ORCID iD

Liang Tang

References

Chowdari

Bamne

Nimgaonkar

VL.

Genetic association studies of antioxidant pathway genes and schizophrenia.

Antioxid Redox Signal 2011; 15: 2037–2045.

De Luca

Tharmalingam

Zai

, et al. Association of HPA axis genes with suicidal behaviour in schizophrenia. J Psychopharmacol 2010; 24: 677–682.

Keck

Kern

Erhardt

, et al. Combined effects of exonic polymorphisms in CRHR1 and AVPR1B genes in a case/control study for panic disorder. Am J Med Genet B Neuropsychiatr Genet 2008; 147B: 1196–1204.

Wang

Murphy

Dow

, et al. Systemic adrenocorticotropic hormone administration down-regulates the expression of corticotropin-releasing hormone (CRH) and CRH-binding protein in infant rat hippocampus. Pediatr Res 2004; 55: 604–610.

Scantlebury

Chun

, et al. Adrenocorticotropic hormone protects learning and memory function in epileptic Kcna1-null mice. Neurosci Lett 2017; 645: 14–18.

Hernan

Alexander

Lenck-Santini

, et al. Attention deficit associated with early life interictal spikes in a rat model is improved with ACTH. PLoS One 2014; 9: e89812.

Yang

Harmon

CM.

Molecular signatures of human melanocortin receptors for ligand binding and signaling.

Biochim Biophys Acta Mol Basis Dis 2017; 1863: 2436–2447.

Fragoso

MCBV

Alencar

Lerario

, et al. Genetics of primary macronodular adrenal hyperplasia. J Endocrinol 2015; 224: R31–R43.

Liu

Fang

, et al . Genetic polymorphisms of mc2r gene associated with responsiveness to adrenocorticotropic hormone therapy in infantile spasms. Chin Med J (Engl) 2008; 121: 1627–1632.

10.

Guo

Deng

Yang

, et al. Melanocortin receptor type 2 (mc2r, acth receptor) expression in patients with alopecia areata. Exp Dermatol 2010; 19: 1020–1022.

11.

Proudnikov

Hamon

Ott

, et al. Association of polymorphisms in the melanocortin receptor type 2 (MC2R, ACTH receptor) gene with heroin addiction. Neurosci Lett 2008; 435: 234–239.

12.

Akin

Coban

, et al. A novel mutation in the MC2R gene causing familial glucocorticoid deficiency type 1. Neonatology 2011; 100: 277–281.

13.

Switonski

Mankowska

Salamon

Family of melanocortin receptor (MCR) genes in mammals-mutations, polymorphisms and phenotypic effects.

J Appl Genet 2013; 54: 461–472.

14.

Richards

JR.

Cannabinoid hyperemesis syndrome: a disorder of the HPA axis and sympathetic nervous system?

Med Hypotheses 2017; 103: 90–95.

15.

Ludwig

Roy

Dwivedi

Role of HPA and the HPG axis interaction in testosterone-mediated learned helpless behavior.

Mol Neurobiol 2018; 56: 394–405.

16.

Kellner

Herzog

Holsboer

, et al. Circadian changes in the sensitivity of the corticotropin-releasing hormones-stimulated HPA system after arginine vasopressin and atrial natriuretic hormone in human male controls. Psychoneuroendocrinology 2018; 20: 515.

17.

Eddy

Wertheim

Hale

, et al. A systematic review and meta-analysis of the effort-reward imbalance model of workplace stress and HPA axis measures of stress. Psychosom Med 2017; 80: 1.

18.

Yang

Chen

Kesterson

Jr , et al. Structural insights into the role of the ACTH receptor cysteine residues on receptor function. Am J Physiol Regul Integr Comp Physiol 2007; 293: R1120–R1126.

19.

Duval

Mokrani

Monreal

, et al. Dopamine and serotonin function in untreated schizophrenia: clinical correlates of the apomorphine and d-fenfluramine tests. Psychoneuroendocrinology 2003; 28: 627–642.

20.

Jiang

Chen

, et al. A study of hypothalamus pituitary adrenal axis function in schizophrenia. Clin Psychiatry 2008; 18: 292–294. http://www/cnki.com.cn/Article/CJFDTotal-LCJS200805004.htm.

21.

Mountjoy

Robbins

Mortrud

, et al. The cloning of a family of genes that encode the melanocortin receptors. Science 1992; 257: 1248–1251.

22.

Metherell

Chan

Clark

AJL.

The genetics of ACTH resistance syndromes.

Best Pract Res Clin Endocrinol Metab 2006; 20: 547–560.

23.

Quintero

Amador-Patarroyo

Montoya-Ortiz

, et al. Autoimmune disease and gender: plausible mechanisms for the female predominance of autoimmunity. J Autoimmun 2012; 38: J109–J119.

24.

Karanikas

Garyfallos

Role of cortisol in patients at risk for psychosis mental state and psychopathological correlates: a systematic review.

Psychiatry Clin Neurosci 2015; 69: 268–282.

Supplementary Material

Please find the following supplemental material available below.

For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.

For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.

0.00 MB

4.74 MB