Genetic variation in the 3'‑untranslated region of PAK1 influences schizophrenia susceptibility

  • Authors:
    • Juan Jiang
    • Jianxiong Long
    • Weijun Ling
    • Guifeng Huang
    • Li Su
  • View Affiliations

  • Published online on: January 12, 2017     https://doi.org/10.3892/etm.2017.4039
  • Pages: 1101-1108
Metrics: Total Views: 0 (Spandidos Publications: | PMC Statistics: )
Total PDF Downloads: 0 (Spandidos Publications: | PMC Statistics: )


Abstract

The present study aimed to investigate the association of two polymorphisms (rs2844337 and rs11237200) in the P21 protein (cell division control protein 42/Rac)‑activated kinase 1 gene with susceptibility to schizophrenia (SCZ) in Chinese Zhuang and Chinese Han populations. A total of 700 patients with SCZ and 700 healthy controls were recruited. Rs2844337 and rs11237200 polymorphisms were genotyped using Sequenom technology. A total of 591 patients completed the Positive and Negative Syndrome Scale (PANSS) assessment. Data were statistically analyzed using PLINK version 1.07 and SPSS version 17.0. In the Chinese Han population, the genotypic (P=0.038) and allelic (P=0.033) frequencies of the 3'‑untranslated region (UTR) genetic variation of rs2844337 in patients were significantly decreased compared to that in controls; these frequencies were significantly associated with SCZ susceptibility in the additive model (Padj=0.032) and in the recessive model (Padj=0.031). Moreover, the TG haplotype constructed by rs2844337 and rs11237200 polymorphisms remained significantly associated with SCZ risk following adjustment for gender and age and applying a Bonferroni correction in the Chinese Han population (Padj=0.003, PBC=0.009). The adjacent 5'‑UTR genetic variation of rs11237200 was significantly associated with the total score (Padj=0.006), positive scale score (Padj=0.014) and general psychopathology scale scores (Padj=0.009) in the recessive model of the Chinese Han population. However, these polymorphisms were not significantly associated with SCZ susceptibility or the PANSS scores in the Chinese Zhuang population. In conclusion, variations in the PAK1 gene influenced the susceptibility and severity of the clinical symptoms of SCZ in the Chinese Han population investigated in the present study.

Introduction

Schizophrenia (SCZ) is characterized by cognitive impairment and disorganized behavior, which may result in debilitation and require costly maintenance. The global prevalence of SCZ varies; the median incidence and the median lifetime morbid risk of SCZ are 15.2/100,000 and 7.2/1,000 individuals, respectively (1,2). In China, a systematic review has demonstrated that the lifetime prevalence of SCZ is 5.44/1,000 individuals (3). In addition, SCZ induces a financial and social burden worldwide (4). Furthermore, the genetic contribution to SCZ is ~81% (5). Genetic studies on SCZ have predominantly focused on the identification of candidate genes and the corresponding genetic loci (69); >1,000 genetic association studies on SCZ risk have also been conducted (http://www.schizophreniaforum.org/). However, studies on candidate genes have yet to clearly elucidate the genetic basis of SCZ (10).

The P21 protein [cell division control protein (Cdc) 42/Rac]-activated kinase 1 (PAK1) gene is located in the 11q13-q14 chromosome, and encodes PAK proteins. These proteins are critical effectors that link the Rho family of GTPases to cytoskeleton reorganization and nuclear signaling, and also function as targets of small GTP binding proteins (Cdc42 and Rac). Rubio et al (11) previously revealed that damage on the Cdc42/PAK1 and Duo/Rac-1/PAK1 pathways in the anterior cingulate cortex may enhance myosin light chain phosphorylation, which is responsible for the cytoskeletal dysfunction and dendritic spine loss associated with SCZ. A study on Caenorhabditis elegans has demonstrated that disrupted-in-SCZ 1 (DISC1) leads to the abnormal branching formation of schizophrenic motor neurons by activating the UNC-73-RAC-PAK signaling pathway (12). Datta et al (13) found that mRNA levels of PAK1 are significantly downregulated in the schizophrenic deep layer of three pyramidal cells. This evidence indicates that the PAK1 gene may be involved in the molecular pathological mechanism of SCZ. However, studies have yet to describe the genetic association of the PAK1 genetic variation and SCZ.

A total of 179 single nucleotide polymorphisms (SNPs) are present in the PAK1 gene. In the rs2844337 polymorphism, a G allele is inserted in the 3′-untranslated region (UTR) in the 77322663 position of the PAK1 gene. Genetic variation in 3′-UTRs of eukaryotic genes may affect gene expression through the binding of these genes to specific microRNAs (miRNAs or miRs), potentially resulting in serious pathology (14). Caputo et al (15) found that two 3′-UTR functional variants (rs11030100 and rs11030099) of the brain-derived neurotrophic factor gene (BDNF) influence SCZ susceptibility by altering miRNA-BDNF binding. In human-induced pluripotent stem cells, a rare functional SNP (rs1130354) in the 3′-UTR of the dopamine receptor D2 gene (DRD2) disrupts the targeting of miR-326 to the mRNA of DRD2; as a result, dopamine signaling in patients with SCZ is dysregulated (16). The G allele of the rs11122396 polymorphism, which is located in the 3′-UTR of the DISC1 gene, may be functionally related to the schizophrenic phenotypes of miR-135b-5p; however, this G allele may no longer regulate the mRNA level of DISC1 (17). Furthermore, rs11237200 is located near the 5′-UTR in the 77474957 position of the PAK1 gene. The ribosomal binding site near the 5′-UTR may be implicated in the inhibition of translation and may accelerate degradation by binding to small RNAs (18). In a meta-analysis, schizophrenic expression quantitative trait loci were suggested to be located near the 5′-UTR (19). Thus, the genetic variation in rs2844337 (3′-UTR) and rs11237200 (5′-UTR) of the PAK1 gene may be involved in SCZ.

Fifty-six ethnic groups are found in China. Among these groups, Han is the largest, followed by Zhuang which ranks second, with a population of 52.82 million (20). More than 90% of the Zhuang population inhabits the Guangxi Zhuang Autonomous Region in South China, which is located near the border of Vietnam. At present, studies on the genetic association of SCZ have predominantly focused on the Han population in China, and few studies have investigated SCZ affecting individuals in other ethnic groups. The authors of the present study have recently published two studies that focus on the association between genetic variations and SCZ susceptibility in the Chinese Zhuang population (21,22). These studies revealed that the rs12807809 polymorphism of the neurogranin gene influences the severity of the clinical symptoms of SCZ in the Chinese Han population, but not in the Chinese Zhuang population (21). Considering ethnic differences, a case-control analysis was performed in the present study to investigate the potential genetic association between SCZ and two PAK1 polymorphisms (rs2844337 and rs11237200) in Chinese Han and Chinese Zhuang populations. Moreover, SCZ predominantly manifests heterogeneous clinical symptoms; genetic variations are possibly related to different clinical symptoms. Therefore, the association between the PAK1 polymorphisms and the clinical symptoms of SCZ were assessed, on the basis of the Positive and Negative Syndrome Scale (PANSS) (23).

Materials and methods

Inclusion and exclusion criteria

In the present study, the Chinese Han or Chinese Zhuang patients were of the same ethnic origin for at least three generations. Patients with SCZ were diagnosed independently by two senior psychiatrists on the basis of the Diagnostic and Statistical Manual of Mental Disorders, 4th edition, criteria (24); the final diagnoses made by these psychiatrists were consistent. Patients who satisfied the following conditions were excluded: Mental disorders caused by organic brain syndrome; mental retardation; nervous system disease; history of severe brain injury; uncooperative patients with severe excitement or impulse; and pregnant or breast-feeding women. Healthy controls were age- and gender-matched with the patients with SCZ. Healthy controls who met the following criteria were excluded: Family history of mental disease; severe head trauma or birth injury; or febrile convulsion in childhood or infant stage.

Subjects

A total of 1,400 unrelated individuals were recruited from the Guangxi Zhuang Autonomous Region between April 2010 and May 2012. The case group (400 Chinese Han and 300 Chinese Zhuang) was recruited from Guangxi Encephalopathy Hospital (Liuzhou, China) and Community Management Patients of Guangxi Liujiang County (Zhuang-inhabited region). The healthy controls (400 Chinese Han and 300 Chinese Zhuang) included individuals who visited Guangxi Liuzhou People's Hospital (Liuzhou, China) and The First Affiliated Hospital of Guangxi University of Chinese Medicine (Nanning, China) for a regular checkup and healthy volunteers living in Liujiang County. All participants provided informed written consent after the protocols were clearly understood or explained by legal guardians. All protocols in the present study were approved by the Ethics Committee of Guangxi Medical University (Nanning, China).

PANSS symptom assessment

PANSS (23) was used to assess the symptoms of patients with SCZ. The following scale scores were calculated: PANSS total score and three subscale scores (positive, negative and general psychopathology scale scores). Symptoms were assessed by two senior psychiatrists on the day that the patients were admitted to hospital.

Genomic DNA extraction

A total of 5 ml peripheral blood was drawn from each subject and collected in vacutainer tubes containing the anticoagulant ethylenediaminetetraacetic acid; the blood specimen was shaken and stored at 4°C. Genomic DNA was isolated from leukocytes in seven days using a TIANamp blood DNA kit (Tiangen Biotech, Beijing, China). Extracted DNA was divided into three tubes and stored at −80°C.

Genotyping

Genotyping was performed using a MassARRAY microarray platform on an iPLEX GOLD technology system (Sequenom, Inc., San Diego, CA, USA). The primers of rs2844337 and rs11237200 polymorphisms were designed using Sequenom Assay Designer 3.1 (Sequenom, Inc.). Primer sequences are listed in Table I. In the present study, primer design, synthesis and genotypic tests were performed by Bomiao Biological Co., Ltd., (Beijing, China). Approximately 5% of the DNA samples were selected randomly and retested to assess the quality of genotypic data; the concordance rate was 100%.

Table I.

Primer sequences of PAK1 gene polymorphisms rs2844337and rs11237200.

Table I.

Primer sequences of PAK1 gene polymorphisms rs2844337and rs11237200.

SNPsForward primer (5′-3′)Reverse primer (5′-3′)
rs2844337 ACGTTGGATGTGTGCTTCTAGCAGCTGTTC ACGTTGGATGCATAGCCGAGAGCATTTCAC
rs11237200 ACGTTGGATGAAGAGTCCGGAATCAGGTTG ACGTTGGATGGGCTTAGAAATGTCAAGGTC

[i] PAK1, P21 protein (Cdc42/Rac)-activated kinase 1; SNP, single nucleotide polymorphism.

Statistical analysis

SPSS version 17.0 (SPSS, Inc., Chicago, IL, USA) was used for evaluating the frequencies of allele and genotype, and the comparison of age and gender. Linkage disequilibrium (LD) analyses were performed with Haploview version 4.2 software (Whitehead Institute for Biomedical Research, Cambridge, MA, USA). The following genetic analyses were conducted using PLINK version 1.07 (The Center for Human Genetic Research, Massachusetts General Hospital and the Broad Institute of Harvard & MIT, Cambridge, MA, USA): Hardy-Weinberg equilibrium (HWE) was detected via a chi-square goodness of fit test; the genetic association of rs2844337 and rs11237200 polymorphisms with SCZ susceptibility was evaluated in terms of the four genetic models (additive, dominant, recessive and allelic models), through unconditional logistic regression analysis; and the correlation of rs2844337 and rs11237200 polymorphisms with PANSS scores was estimated through linear regression analysis. P-values were two tailed, and P<0.05 was considered to indicate a statistically significant difference. In addition, Bonferroni correction was used for multiple testing. Furthermore, log-odds (LOD) scores were also calculated using Haploview. LOD >1 was considered to indicate a linkage disequilibrium, while LOD <-2 was considered to indicate that a linkage disequilibrium did not exist.

Results

Subject characteristics

In the Chinese Han population, the SCZ patient group consisted of 269 males and 131 females with an average age of 32.29±11.56 years, whereas the healthy control group comprised 248 males and 152 females with an average age of 33.09±11.17 years. In the Chinese Han population, the patients with SCZ did not significantly differ from the healthy controls in terms of gender (χ2=2.411; P=0.120) and age (t=0.992; P=0.321). In the Chinese Zhuang population, the SCZ patient group contained 207 males and 93 females with an average age of 33.68±11.99 years; the healthy control group included 198 males and 102 females with an average age of 32.37±12.27 years. In the Chinese Zhuang population, the patients with SCZ did not significantly differ from the healthy controls in terms of gender (χ2=0.615; P=0.894) and age (t=1.319; P=0.188). In both populations, the patients with SCZ did not significantly differ from the healthy controls in terms of gender (χ2=2.859; P=0.091) and age (t=−0.167; P=0.868).

A total of 591 patients with SCZ completed the PANSS assessment. Chinese Han patients included 385 individuals (216 males and 124 females) with an average age of 32.30±11.66 years; the Chinese Zhuang patients consisted of 206 individuals (151 males and 55 females) with an average age of 31.85±11.28 years. The Chinese Zhuang patients did not significantly differ from the Chinese Han patients in terms of the distribution of gender (χ2=1.929; P=0.165) or age (t=0.449; P=0.654).

HWE evaluation

In the healthy control group, the genotypic distributions of the rs2844337 and rs11237200 polymorphisms did not significantly deviate from the HWE (PHWE>0.05; Table II).

Table II.

Distribution of genotype and allele frequency and HWE test for rs2844337 and rs11237200.

Table II.

Distribution of genotype and allele frequency and HWE test for rs2844337 and rs11237200.

Genotype (A1A1/A1A2/A2A2) Allele (A1/A2)


SNPsEthnicityCaseControlχ2P-value PHWECaseControlχ2P-value
rs2844337Han2/64/23510/74/2206.5380.0380.26968/53494/5144.5250.033
(G/T)Zhuang6/61/1617/61/1590.0870.9570.63573/38375/3790.0440.835
Total8/125/39617/135/3793.9940.1360.256141/917169/8932.8400.092
rs11237200Han21/142/23720/161/2172.0920.3510.184184/616201/5951.1050.293
(A/G)Zhuang17/96/18316/107/1760.7480.6881.000130/462139/4590.2810.596
Total38/238/42036/268/3932.7290.2560.304314/1078340/10541.3020.254

[i] SNP, single nucleotide polymorphisms; A1, minor allele; A2, major allele; Case, subjects diagnosed with schizophrenia; PHWE: P-value for Hardy-Weinberg equilibrium test; G, guanine; T, thymine; A, adenine. Bold text indicates a statistically significant difference (P<0.05) between the case and control groups.

Distribution of genotype and allele frequencies

Genotypic (P=0.038) and allelic (P=0.033) frequencies of the 3′-UTR of the rs2844337 polymorphism in patients significantly decreased compared to that in controls in the Chinese Han population but not in the Chinese Zhuang population (Table II). By contrast, the genotypic and allelic frequencies of the rs11237200 polymorphism did not significantly differ in either population.

Genetic association between polymorphism and SCZ susceptibility

In the Chinese Han population, the rs2844337 polymorphism was significantly associated with SCZ susceptibility in additive (Padj=0.032), recessive (Padj=0.031) and allelic (P=0.033) models; however, these differences were not statistically significant following Bonferroni correction (additive, PBC=0.064; recessive, PBC=0.062; allelic, PBC=0.066). Conversely, the rs11237200 polymorphism was not significantly associated with SCZ susceptibility in any of the three genetic models. Likewise, these two polymorphisms were not significantly associated with SCZ susceptibility in the Chinese Zhuang population (Table III).

Table III.

Association of SNP rs2844337and rs11237200 with schizophrenia risk in different ethnic group.

Table III.

Association of SNP rs2844337and rs11237200 with schizophrenia risk in different ethnic group.

CrudeAdjusted


EthnicitySNPModelOR (95% CI)P-valueOR (95% CI) Padj
Hanrs2844337Additive0.699 (0.501–0.977)0.0360.694 (0.497–0.970)0.032
Dominant0.736 (0.508–1.066)0.1050.733 (0.505–1.062)0.101
Recessive0.197 (0.043–0.905)0.0370.186 (0.040–0.860)0.031
Allelic0.696 (0.498–0.973)0.033
rs11237200Additive0.880 (0.696–1.112)0.2840.879 (0.695–1.111)0.279
Dominant0.824 (0.623–1.092)0.1780.822 (0.621–1.089)0.172
Recessive1.047 (0.559–1.964)0.8861.050 (0.560–1.970)0.879
Allelic0.884 (0.703–1.112)0.293
Zhuangrs2844337Additive0.964 (0.680–1.367)0.8360.985 (0.693–1.400)0.933
Dominant0.973 (0.651–1.455)0.8941.005 (0.670–1.508)0.980
Recessive0.849 (0.281–2.568)0.7730.828 (0.273–2.516)0.740
Allelic0.963 (0.677–1.370)0.835
rs11237200Additive0.931 (0.712–1.218)0.6010.926 (0.708–1.212)0.576
Dominant0.884 (0.636–1.227)0.4600.878 (0.631–1.220)0.437
Recessive1.078 (0.534–2.176)0.8351.073 (0.531–2.169)0.843
Allelic0.929 (0.708–1.219)0.596

[i] SNP, single nucleotide polymorphism; OR, odds ratio; CI, confidence interval; A1, minor allele; A2, major allele; Additive, A1A1/A1A2/A2A2; Dominant, (A1A1+A1A2)/A2A2; Recessive, A1A1/(A1A2+A2A2); Allelic, A1/A2; Padj, P-value adjusted for age and gender. Bold text indicates a statistically significant association (P<0.05) between SNP and schizophrenia susceptibility.

LD and haplotype analysis

LD analyses were performed with D' values of 0.881 and 0.802 in the Chinese Han and Chinese Zhuang samples, respectively, suggesting that LD may exist between rs2844337 and rs11237200. Detailed results are presented in Table IV. In addition, three haplotypes were constructed using the rs2844337 and rs11237200 polymorphisms. In the Chinese Han population, the TG haplotype posed a significantly higher SCZ risk compared with other haplotypes (ORadj=1.52; Padj=0.003). This risk remained statistically significant following Bonferroni correction (PBC=0.009). By contrast, these haplotypes were not significantly associated with SCZ risk in the Chinese Zhuang population (Table V).

Table IV.

Linkage disequilibrium for rs2844337 and rs11237200.

Table IV.

Linkage disequilibrium for rs2844337 and rs11237200.

SamplesD'D' conf. boundsLODr2
Han0.8810.65–0.966.440.026
Zhuang0.8020.53–0.924.510.025

[i] LOD, log-odds.

Table V.

Association of haplotype rs2844337/rs11237200 with schizophrenia risk.

Table V.

Association of haplotype rs2844337/rs11237200 with schizophrenia risk.

CrudeAdjusted


EthnicityHaplotypeORP-value ORadj Padj
HanT-A0.7910.1590.7750.129
G-G0.7090.0480.720.062
T-G1.510.0031.520.003
ZhuangT-A0.7440.1320.7370.122
G-G0.9920.9631.010.952
T-G1.250.1571.240.178

[i] Adjusted, adjusted for gender and age; OR, odds ratio; Padj, P-value adjusted for gender and age; T, thymine; A, adenine; G, guanine. Bold text indicates a statistically significant association (P<0.05) between haplotype and schizophrenia susceptibility.

Association of the PAK1 polymorphisms with PANSS scores

In the Chinese Han population, rs11237200 was significantly associated with the total score (Padj=0.006), positive scale score (Padj=0.014) and general psychopathology scale score (Padj=0.009) in the recessive model. By contrast, rs2844337 was not significantly associated with the PANSS score in any of the three genetic models (Table VI). In the Chinese Zhuang population, the PAK1 polymorphisms were not significantly associated with the PANSS scores in any of the three genetic models (Table VII).

Table VI.

Linear regression results for genotypic association between rs2844337/rs11237200 and PANSS scores in the Han population (adjusted for age and gender).

Table VI.

Linear regression results for genotypic association between rs2844337/rs11237200 and PANSS scores in the Han population (adjusted for age and gender).

Additive modelDominant modelRecessive model



SNPVariablesβ (95% CI) Padjβ (95% CI) Padjβ (95% CI) Padj
rs2844337Total score1.904 (−4.035–7.844)0.5302.072 (−4.159–8.304)0.5150.608 (−30.810–32.020)0.970
Positive0.882 (−1.404–3.168)0.4500.875 (−1.524–3.274)0.4752.435 (−9.665–14.530)0.694
Negative0.145 (−1.824–2.113)0.8860.076 (−1.989–2.142)0.9422.102 (−8.305–12.510)0.693
General psychopathology0.852 (−2.014–3.719)0.5611.094 (−1.912–4.101)0.476−3.963 (−19.120–11.190)0.609
rs11237200Total score0.818 (−2.936–4.573)0.670−1.733 (−6.305–2.838)0.45813.770 (3.977–23.560)0.006
Positive0.583 (−0.828–1.994)0.418−0.130 (−1.849–1.590)0.8824.653 (0.963–8.342)0.014
Negative−0.301 (−1.561–0.960)0.640−1.047 (−2.579–0.484)0.1812.813 (−0.496–6.122)0.097
General psychopathology0.516 (−1.292–2.324)0.576−0.579 (−2.782–1.623)0.6066.283 (1.563–11.000)0.009

[i] SNP, single nucleotide polymorphisms; Additive model, A1A1/A1A2/A2A2; Dominant model, (A1A1+A1A2)/A2A2; Recessive model, A1A1/(A1A2+A2A2); A1, minor allele; A2, major allele; β, partial regression coefficent of the linear regression results; CI, confidence interval of partial regression coefficient; Padj, P-value adjusted for age and gender; PANSS, Positive and Negative Syndrome Scale. Bold text indicates a statistically significant genotypic association (P<0.05) between SNP and PANSS for schizophrenia.

Table VII.

Linear regression results for genotypic association between rs2844337/rs11237200 and PANSS scores in the Zhuang population (adjusted for age and gender).

Table VII.

Linear regression results for genotypic association between rs2844337/rs11237200 and PANSS scores in the Zhuang population (adjusted for age and gender).

Additive modelDominant modelRecessive model



SNPVariablesβ (95% CI) Padjβ (95% CI) Padjβ (95% CI) Padj
rs2844337Total score−0.024 (−8.092–8.045)0.9950.744 (−7.720–9.207)0.864−23.010 (−69.140–23.130)0.330
Positive0.673 (−2.409–3.756)0.6690.837 (−2.395–4.068)0.613−2.931 (−20.750–14.890)0.748
Negative−0.397 (−2.983–2.188)0.764−0.177 (−2.889–2.535)0.898−7.896 (−22.790–7.000)0.301
General psychopathology−0.395 (−4.432–3.642)0.848−0.026 (−4.262–4.209)0.990−12.210 (−35.280–10.870)0.302
rs11237200Total score0.775 (−4.406–5.955)0.770−1.211 (−7.609–5.187)0.71110.350 (−2.881–23.580)0.127
Positive−0.707 (−2.692–1.278)0.486−1.816 (−4.257–0.626)0.1473.196 (−1.896–8.288)0.220
Negative1.695 (−0.028–3.418)0.0551.661 (−0.470–3.793)0.1283.998 (−0.449–8.446)0.080
General psychopathology−0.228 (−2.789–2.332)0.861−1.082 (−4.236–2.072)0.5023.202 (−3.374–9.779)0.341

[i] SNP, single nucleotide polymorphisms; Additive model, A1A1/A1A2/A2A2; Dominant model, (A1A1+A1A2)/A2A2; Recessive model, A1A1/(A1A2+A2A2); A1, minor allele; A2, major allele; β, partial regression coefficent of the linear regression results; CI, confidence interval of partial regression coefficient; Padj, P-value adjusted for age and gender; PANSS, Positive and Negative Syndrome Scale. Bold text indicates a statistically significant genotypic association (P<0.05) between SNP and PANSS for schizophrenia.

Discussion

To the best of our knowledge, this research is the first to investigate the association of rs2844337 with susceptibility to SCZ, and the association of rs11237200 with the clinical symptoms of SCZ. The present results indicated that the rs2844337 polymorphism of PAK1 was a risk factor of the susceptibility to SCZ in the Chinese Han, but not in the Chinese Zhuang, population. Moreover, the results of the quantitative trait locus analysis suggested that rs11237200 influenced the SCZ symptoms of the Chinese Han but not of the Chinese Zhuang population.

Notably, rs2844337 and rs11237200 are located in the PAK1 gene, which encodes a family member of serine/threonine p21 activating kinases, namely PAK proteins (25). The dysregulation of PAK proteins contributes to the regulation of actin stabilization and the reduction of dendritic spines in patients with SCZ (26). PAK proteins can phosphorylate and inactivate LIM-kinase during actin depolymerization, and alter F-actin stabilization (27). Dominant negative PAK1 can reduce the number of basal dendrites and the number of the primary branches on apical dendrites in immature cortical neurons, and can alter the morphological characteristics of hippocampal neurons in the spine (28). Huang et al (29) demonstrated that dendritic axons are markedly simplified, and synapse density is reduced in PAK1/PAK3 double-knockout mice. Based on these results, the present hypothesis is that the PAK1 gene may be involved in the molecular pathological mechanism of SCZ.

3′-UTR has an important role in the post-transcriptional regulation of gene expression; in particular, 3′-UTR binds to miRNA. The genetic variation in this region induces various pathophysiological changes (30). Rs10759, a genetic variation in the 3′-UTR of the regulator of G protein signaling 4 (RGS4), is associated with SCZ (31). In a previous in vitro luciferase reporter assay, rs10759 was demonstrated to influence SCZ susceptibility as miR-124 repressed the binding of the 3′-UTR of the RGS4 gene to its target mRNA (32). Saetre et al (33) found that the polymorphism of DISC1 in the 3′-UTR enhances susceptibility to SCZ. Rossi et al (17) conducted a functional study and demonstrated that genetic variation in the 3′-UTR of DISC1 increases the gene expression levels and affects SCZ-associated phenotypes by disrupting miRNA-mRNA target recognition. In the present study, the 3′-UTR of the rs2844337 polymorphism of the PAK1 gene was significantly associated with the susceptibility of the Chinese Han population to SCZ. However, the biological function of rs2844337 polymorphism is yet to be verified. Furthermore, studies have yet to determine whether rs2844337 participates in the SCZ pathomechanism by influencing the binding of miRNA to the PAK1 gene, thereby affecting PAK1 expression levels and stability. Therefore, further study is required to verify this mechanism.

Considering the rs11237200 polymorphism near the 5′-UTR, it was indicated that this polymorphism did not influence SCZ susceptibility; alternatively, this polymorphism affected the quantitative traits of SCZ susceptibility. Nadalin et al (34) demonstrated that the peroxisome proliferator-activated receptor α-L162V polymorphism is not associated with SCZ; however, it influenced the overall PANSS psychopathology score. Wang et al (35) reported that the allelic variant of rs66642155 is not significantly associated with SCZ; nevertheless, this variant affects the positive symptoms of SCZ. SCZ is considered as a clinically heterogeneous disease, and the symptom heterogeneity of SCZ is defined as a symptom dimension (36). Symptom dimensions represent multiple sets of symptoms, which usually occur simultaneously and account for the majority of the phenotypic variations detected among patients (37). Different symptom dimensions are potentially under the control of a set of genes; these dimensions then facilitate the genetic distribution of SCZ (3739). Therefore, the association between the genetic variations and the symptoms assessed by PANSS may indicate the genetic mechanism of SCZ. At present, the diagnosis of SCZ is primarily based on descriptive clinical criteria (40). Therefore, it can be suggested that traditional diagnostic criteria should be replaced with a novel classification and diagnostic criteria based on symptom dimensions. The genetic variations associated with the clinical symptoms of SCZ should also provide the genetic basis of novel criteria.

The present results showed that the rs2844337 and rs11237200 polymorphisms investigated herein influenced SCZ susceptibility and/or quantitative traits in the Chinese Han population, but not in the Chinese Zhuang population. This disparity is attributed to differing genetic effects on the Chinese Han and Chinese Zhuang populations. Genetic effects on complex diseases are heterogeneous among different ethnic groups (41). For instance, the Chinese Zhuang population has higher lipid levels, lower hyperlipidemia prevalence and higher blood pressure compared with the Chinese Han population (4244). Studies on the genetic association of these two populations have revealed that multiple genetic variants, including polymorphisms in the lipoprotein lipase gene, hepatic lipase gene and sterol regulatory element-binding protein-2 gene, of the Chinese Zhuang population are distinct from those of the Chinese Han population (4547). However, it is possible that in the present study this ethnic difference was confounded by the smaller sample size of the Chinese Zhuang subjects, which may result in a relatively weaker power of the Chinese Zhuang samples.

In the present study, several limitations should be considered. Firstly, rs2844337 is a functional site located in the 3′-UTR of the PAK1 gene, and was significantly associated with SCZ susceptibility. However, the biological function of rs2844337 remains unclear. As such, further studies should explore the effect of rs2844337 on the binding force of miRNA to the PAK1 gene, and on the expression and translation of the PAK1 gene. Secondly, the sample size of the Chinese Zhuang was not equal to that of the Chinese Han. Therefore, the results should be interpreted cautiously; future studies should be conducted using a larger sample size of Chinese Zhuang individuals. Thirdly, only two SNPs on the PAK1 gene were investigated in the analysis of the association of rs2844337 and rs11237200 with SCZ, which cannot represent all genetic variations in the 3′-UTR and/or 5′-UTR of the PAK1 gene. Therefore, it is necessary for further studies to focus on whether other PAK1 gene SNPs affect the susceptibility and severity of the clinical symptoms of SCZ in the Chinese Han population. Fourthly, genetic effects should be determined to explain the cause of SCZ. However, the role of environmental factors, including maternal exposure to viruses during pregnancy (48), social disorganization of neighborhoods (49) and prenatal stress (50) should be considered.

In conclusion, the present study supports the hypothesis that genetic variations in the PAK1 gene influence the susceptibility and severity of the clinical symptoms of SCZ in the Chinese Han population. Further studies should be conducted to investigate the biological function of the rs2844337 polymorphism and the pathomechanism involved in SCZ.

Acknowledgements

The present study was supported by the Guangxi Science Foundation (grant no. 2013GXNSFAA019352); the National Natural Science Foundation of China (grant no. 81460518); the Science and Technology Research Projects of Guangxi University (grant no. 2013YB043); and the Youth Science Foundation of Guangxi Medical University (grant no. GXMUYSF201322).

References

1 

McGrath J, Saha S, Chant D and Welham J: Schizophrenia: A concise overview of incidence, prevalence, and mortality. Epidemiol Rev. 30:67–76. 2008. View Article : Google Scholar : PubMed/NCBI

2 

van Os J and Kapur S: Schizophrenia. Lancet. 374:635–645. 2009. View Article : Google Scholar : PubMed/NCBI

3 

Long J, Huang G, Liang W, Liang B, Chen Q, Xie J, Jiang J and Su L: The prevalence of schizophrenia in mainland China: Evidence from epidemiological surveys. Acta Psychiatr Scand. 130:244–256. 2014. View Article : Google Scholar : PubMed/NCBI

4 

Mueser KT and McGurk SR: Schizophrenia. Lancet. 363:2063–2072. 2004. View Article : Google Scholar : PubMed/NCBI

5 

Sullivan PF: The genetics of schizophrenia. PLoS Med. 2:e2122005. View Article : Google Scholar : PubMed/NCBI

6 

O'Donovan MC, Craddock N, Norton N, Williams H, Peirce T, Moskvina V, Nikolov I, Hamshere M, Carroll L, Georgieva L, et al: Identification of loci associated with schizophrenia by genome-wide association and follow-up. Nat Genet. 40:1053–1055. 2008. View Article : Google Scholar : PubMed/NCBI

7 

Harrison PJ and Owen MJ: Genes for schizophrenia? Recent findings and their pathophysiological implications. Lancet. 361:417–419. 2003. View Article : Google Scholar : PubMed/NCBI

8 

Stefansson H, Ophoff RA, Steinberg S, Andreassen OA, Cichon S, Rujescu D, Werge T, Pietiläinen OP, Mors O, Mortensen PB, et al: Common variants conferring risk of schizophrenia. Nature. 460:744–747. 2009.PubMed/NCBI

9 

Schizophrenia Working Group of the Psychiatric Genomics Consortium, . Biological insights from 108 schizophrenia-associated genetic loci. Nature. 511:421–427. 2014. View Article : Google Scholar : PubMed/NCBI

10 

Farrell MS, Werge T, Sklar P, Owen MJ, Ophoff RA, O'Donovan MC, Corvin A, Cichon S and Sullivan PF: Evaluating historical candidate genes for schizophrenia. Mol Psychiatry. 20:555–562. 2015. View Article : Google Scholar : PubMed/NCBI

11 

Rubio MD, Haroutunian V and Meador-Woodruff JH: Abnormalities of the Duo/Ras-related C3 botulinum toxin substrate 1/p21-activated kinase 1 pathway drive myosin light chain phosphorylation in frontal cortex in schizophrenia. Biol Psychiatry. 71:906–914. 2012. View Article : Google Scholar : PubMed/NCBI

12 

Chen SY, Huang PH and Cheng HJ: Disrupted-in-Schizophrenia 1-mediated axon guidance involves TRIO-RAC-PAK small GTPase pathway signaling. Proc Natl Acad Sci USA. 108:5861–5866. 2011. View Article : Google Scholar : PubMed/NCBI

13 

Datta D, Arion D, Corradi JP and Lewis DA: Altered expression of CDC42 signaling pathway components in cortical layer 3 pyramidal cells in schizophrenia. Biol Psychiatry. 78:775–785. 2015. View Article : Google Scholar : PubMed/NCBI

14 

Lu J and Clark AG: Impact of microRNA regulation on variation in human gene expression. Genome Res. 22:1243–1254. 2012. View Article : Google Scholar : PubMed/NCBI

15 

Caputo V, Sinibaldi L, Fiorentino A, Parisi C, Catalanotto C, Pasini A, Cogoni C and Pizzuti A: Brain derived neurotrophic factor (BDNF) expression is regulated by microRNAs miR-26a and miR-26b allele-specific binding. PLoS One. 6:e286562011. View Article : Google Scholar : PubMed/NCBI

16 

Shi S, Leites C, He D, Schwartz D, Moy W, Shi J and Duan J: MicroRNA-9 and microRNA-326 regulate human dopamine D2 receptor expression, and the microRNA-mediated expression regulation is altered by a genetic variant. J Biol Chem. 289:13434–13444. 2014. View Article : Google Scholar : PubMed/NCBI

17 

Rossi M, Kilpinen H, Muona M, Surakka I, Ingle C, Lahtinen J, Hennah W, Ripatti S and Hovatta I: Allele-specific regulation of DISC1 expression by miR-135b-5p. Eur J Hum Genet. 22:840–843. 2014. View Article : Google Scholar : PubMed/NCBI

18 

Lavi-Itzkovitz A, Peterman N, Jost D and Levine E: Quantitative effect of target translation on small RNA efficacy reveals a novel mode of interaction. Nucleic Acids Res. 42:12200–12211. 2014. View Article : Google Scholar : PubMed/NCBI

19 

Kim Y, Xia K, Tao R, Giusti-Rodriguez P, Vladimirov V, van den Oord E and Sullivan PF: A meta-analysis of gene expression quantitative trait loci in brain. Transl Psychiatry. 4:e4592014. View Article : Google Scholar : PubMed/NCBI

20 

China Bureau of Statistics, . China statistical yearbookPopulation. China Statistics press; Nanning: 2014

21 

Su L, Long J, Liang B, Gu L, Pan R, Li L, Zhou Y, Tang X, Jiang J, Chen Q, et al: Association between the NOTCH4 gene rs3131296 polymorphism with schizophrenia risk in the Chinese Zhuang population and Chinese Han population. Acta Neuropsychiatr. 26:240–245. 2014. View Article : Google Scholar : PubMed/NCBI

22 

Su L, Long J, Pan R, Xie X, Mao X, Zhou Y, Chen Q and Wei B: Influence of NRGN rs12807809 polymorphism on symptom severity in individuals with schizophrenia in the Han population but not the Zhuang population of south China. Acta Neuropsychiatr. 27:221–227. 2015. View Article : Google Scholar : PubMed/NCBI

23 

Kay SR, Fiszbein A and Opler LA: The positive and negative syndrome scale (PANSS) for schizophrenia. Schizophr Bull. 13:261–276. 1987. View Article : Google Scholar : PubMed/NCBI

24 

Allen F, Avram HM, Ruth R and Michael BF: The DSM-IV Classification and PsychopharmacologyThe Fourth Generation of Progress. Raven Press; New York, NY: 1995

25 

Zheng M, Liu J, Zhu M, Yin R, Dai J, Sun J, Shen W, Ji Y, Jin G, Ma H, et al: Potentially functional polymorphisms in PAK1 are associated with risk of lung cancer in a Chinese population. Cancer Med. 4:1781–1787. 2015. View Article : Google Scholar : PubMed/NCBI

26 

Konopaske GT, Lange N, Coyle JT and Benes FM: Prefrontal cortical dendritic spine pathology in schizophrenia and bipolar disorder. JAMA Psychiatry. 71:1323–1331. 2014. View Article : Google Scholar : PubMed/NCBI

27 

Edwards DC, Sanders LC, Bokoch GM and Gill GN: Activation of LIM-kinase by Pak1 couples Rac/Cdc42 GTPase signalling to actin cytoskeletal dynamics. Nat Cell Biol. 1:253–259. 1999. View Article : Google Scholar : PubMed/NCBI

28 

Hayashi K, Ohshima T, Hashimoto M and Mikoshiba K: Pak1 regulates dendritic branching and spine formation. Dev Neurobiol. 67:655–669. 2007. View Article : Google Scholar : PubMed/NCBI

29 

Huang W, Zhou Z, Asrar S, Henkelman M, Xie W and Jia Z: p21-Activated kinases 1 and 3 control brain size through coordinating neuronal complexity and synaptic properties. Mol Cell Biol. 31:388–403. 2011. View Article : Google Scholar : PubMed/NCBI

30 

Conne B, Stutz A and Vassalli JD: The 3′ untranslated region of messenger RNA: A molecular ‘hotspot’ for pathology? Nat Med. 6:637–641. 2000. View Article : Google Scholar : PubMed/NCBI

31 

Campbell DB, Ebert PJ, Skelly T, Stroup TS, Lieberman J, Levitt P and Sullivan PF: Ethnic stratification of the association of RGS4 variants with antipsychotic treatment response in schizophrenia. Biol Psychiatry. 63:32–41. 2008. View Article : Google Scholar : PubMed/NCBI

32 

Gong Y, Wu CN, Xu J, Feng G, Xing QH, Fu W, Li C, He L and Zhao XZ: Polymorphisms in microRNA target sites influence susceptibility to schizophrenia by altering the binding of miRNAs to their targets. Eur Neuropsychopharmacol. 23:1182–1189. 2013. View Article : Google Scholar : PubMed/NCBI

33 

Saetre P, Agartz I, De Franciscis A, Lundmark P, Djurovic S, Kähler A, Andreassen OA, Jakobsen KD, Rasmussen HB, Werge T, et al: Association between a disrupted-in-schizophrenia1 (DISC1) single nucleotide polymorphism and schizophrenia in a combined Scandinavian case-control sample. Schizophr Res. 106:237–241. 2008. View Article : Google Scholar : PubMed/NCBI

34 

Nadalin S, Giacometti J and Buretić-Tomljanović A: PPARα-L162V polymorphism is not associated with schizophrenia risk in a Croatian population. Prostaglandins Leukot Essent Fatty Acids. 91:221–225. 2014. View Article : Google Scholar : PubMed/NCBI

35 

Wang S, Li W, Zhang H, Wang X, Yang G, Zhao J, Yang Y and Lv L: Association of microRNA137 gene polymorphisms with age at onset and positive symptoms of schizophrenia in a Han Chinese population. Int J Psychiatry Med. 47:153–168. 2014. View Article : Google Scholar : PubMed/NCBI

36 

Andreasen NC and Carpenter WT Jr: Diagnosis and classification of schizophrenia. Schizophr Bull. 19:199–214. 1993. View Article : Google Scholar : PubMed/NCBI

37 

Rietkerk T, Boks MP, Sommer IE, Liddle PF, Ophoff RA and Kahn RS: The genetics of symptom dimensions of schizophrenia: Review and meta-analysis. Schizophr Res. 102:197–205. 2008. View Article : Google Scholar : PubMed/NCBI

38 

Boks MP, Leask S, Vermunt JK and Kahn RS: The structure of psychosis revisited: The role of mood symptoms. Schizophr Res. 93:178–185. 2007. View Article : Google Scholar : PubMed/NCBI

39 

Labbe A, Bureau A, Moreau I, Roy MA, Chagnon Y, Maziade M and Merette C: Symptom dimensions as alternative phenotypes to address genetic heterogeneity in schizophrenia and bipolar disorder. Eur J Hum Genet. 20:1182–1188. 2012. View Article : Google Scholar : PubMed/NCBI

40 

World Health Organization, . The ICD-10 classification of mental and behavioural disorders: Clinical descriptions and diagnostic guidelines. Geneva: World Health Organization. 1993.

41 

Ioannidis JP, Ntzani EE and Trikalinos TA: ‘Racial’ differences in genetic effects for complex diseases. Nat Genet. 36:1312–1318. 2004. View Article : Google Scholar : PubMed/NCBI

42 

Yin R, Chen Y, Pan S, He F, Liu T, Yang D, Wu J, Yao L, Lin W, Li R and Huang J: Comparison of lipid levels, hyperlipidemia prevalence and its risk factors between Guangxi Hei Yi Zhuang and Han populations. Arch Med Res. 37:787–793. 2006. View Article : Google Scholar : PubMed/NCBI

43 

Ruixing Y, Hui L, Jinzhen W, Weixiong L, Dezhai Y, Shangling P, Jiandong H and Xiuyan L: Association of diet and lifestyle with blood pressure in the Guangxi Hei Yi Zhuang and Han populations. Public Health Nutr. 12:553–561. 2009. View Article : Google Scholar : PubMed/NCBI

44 

Yin R, Li H, Wu J, Lin W, Yang D, Pan S, Huang J and Long X: Effects of alcohol consumption and other lifestyle behaviors on blood pressure for the middle-aged and elderly in the Guangxi Hei Yi Zhuang and Han populations. Alcohol. 41:541–550. 2007. View Article : Google Scholar : PubMed/NCBI

45 

Yin R, Wang Y, Chen G, Lin W, Yang D and Pan S: Lipoprotein lipase gene polymorphism at the PvuII locus and serum lipid levels in Guangxi Hei Yi Zhuang and Han populations. Clin Chem Lab Med. 44:1416–1421. 2006. View Article : Google Scholar : PubMed/NCBI

46 

Wu J, Yin R, Lin W, Pan S and Yang D: Hepatic lipase gene −514C/T polymorphism in the Guangxi Hei Yi Zhuang and Han populations. Lipids. 43:733–740. 2008. View Article : Google Scholar : PubMed/NCBI

47 

Yaju D, Ruixing Y, Yiyang L, Yijiang Z, Weixiong L, Shangling P and Dezhai Y: Polymorphism of the sterol regulatory element-binding protein-2 gene and its association with serum lipid levels in the Guangxi Hei Yi Zhuang and Han populations. Am J Med Sci. 337:14–22. 2009. View Article : Google Scholar : PubMed/NCBI

48 

Blomström A, Karlsson H, Wicks S, Yang S, Yolken RH and Dalman C: Maternal antibodies to infectious agents and risk for non-affective psychoses in the offspring-a matched case-control study. Schizophr Res. 140:25–30. 2012. View Article : Google Scholar : PubMed/NCBI

49 

Veling W, Susser E, Selten JP and Hoek HW: Social disorganization of neighborhoods and incidence of psychotic disorders: A 7-year first-contact incidence study. Psychol Med. 45:1789–1798. 2015. View Article : Google Scholar : PubMed/NCBI

50 

Markham JA and Koenig JI: Prenatal stress: Role in psychotic and depressive diseases. Psychopharmacology (Berl). 214:89–106. 2011. View Article : Google Scholar : PubMed/NCBI

Related Articles

Journal Cover

March-2017
Volume 13 Issue 3

Print ISSN: 1792-0981
Online ISSN:1792-1015

Sign up for eToc alerts

Recommend to Library

Copy and paste a formatted citation
x
Spandidos Publications style
Jiang J, Long J, Ling W, Huang G and Su L: Genetic variation in the 3'‑untranslated region of PAK1 influences schizophrenia susceptibility. Exp Ther Med 13: 1101-1108, 2017
APA
Jiang, J., Long, J., Ling, W., Huang, G., & Su, L. (2017). Genetic variation in the 3'‑untranslated region of PAK1 influences schizophrenia susceptibility. Experimental and Therapeutic Medicine, 13, 1101-1108. https://doi.org/10.3892/etm.2017.4039
MLA
Jiang, J., Long, J., Ling, W., Huang, G., Su, L."Genetic variation in the 3'‑untranslated region of PAK1 influences schizophrenia susceptibility". Experimental and Therapeutic Medicine 13.3 (2017): 1101-1108.
Chicago
Jiang, J., Long, J., Ling, W., Huang, G., Su, L."Genetic variation in the 3'‑untranslated region of PAK1 influences schizophrenia susceptibility". Experimental and Therapeutic Medicine 13, no. 3 (2017): 1101-1108. https://doi.org/10.3892/etm.2017.4039