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

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.

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 r2≥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 χ2 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 (padj) in all cases in which an initial significant relationship was found. Squared correlation coefficients (r2) 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 (r2>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/m2)	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	Powerb	HWE
				Case (n=234)	Control (n=399)
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). HaplotypeACGGTGT 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 haplotypeACGATAC 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)
	AA	AB	BB	AA	AB	BB	pa	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).

Haplotypea	Case	Control	p b	padjc	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)
		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 haplotypeACGGTGT 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 haplotypeACGGTGT 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.

Haplotypea	Subgroups	Case	Control	p b	padjc	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.

Click here for additional data file.

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