PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of jmedgeneJournal of Medical GeneticsVisit this articleSubmit a manuscriptReceive email alertsContact usBMJ
 
J Med Genet. 2007 August; 44(8): 532–536.
Published online 2007 April 5. doi:  10.1136/jmg.2006.047944
PMCID: PMC2597930

Connexin 50 gene on human chromosome 1q21 is associated with schizophrenia in matched case–control and family‐based studies

Abstract

Background

The gap junction subunit connexin permits direct intercellular exchange of ions and molecules including glutamate, and plays an important role in the central nervous system. The connexin 40 (Cx40) and connexin 50 (Cx50) genes are located on chromosome 1q21.1, a region strongly linked with schizophrenia. These lines of evidence suggest that Cx40 and Cx50 may play a role in schizophrenia.

Methods

Using an allele‐specific PCR assay, four polymorphisms each were genotyped for Cx40 and Cx50 in 190 Caucasian patients with schizophrenia and 190 controls matched for sex, age and ethnicity. Following up, Cx50 rs989192 and rs4950495 were investigated in 99 Canadian and 163 Portuguese trios and nuclear families with schizophrenia probands. Hardy–Weinberg equilibrium and linkage disequilibrium (LD) block identification was carried out with HaploView, and association analysis for alleles and haplotypes with a permutation test of 10 000 simulations was carried out using the UNPHASED software program.

Results

Distributions of genotype frequencies of all markers were in Hardy–Weinberg equilibrium in Caucasian patients, controls and families. One rs989192‐rs4950495 LD block was found in patients but not in controls. We found a significant association between the Cx50 rs989192‐rs4950495 haplotype and schizophreniay (χ2 = 29.55, p<0.01). The A‐C haplotype had a higher frequency in patients (χ2 = 7.153, p<0.01). Family studies also showed that the A‐C haplotype was transmitted more often to patients with schizophrenia (χ2 = 8.43, p<0.01). No association of Cx40 with schizophrenia was found for allele, genotype or haplotype analyses.

Conclusions

Our matched case–control and family study indicate that Cx50, but not Cx40, may play a role in the genetic susceptibility to schizophrenia.

Keywords: connexin 50, connexin 40, schizophrenia, matched case–control study, family study

Gap junctions (GJs) are channel‐forming structures in contacting plasma membranes. Two hemichannels, each built of six connexin subunits, can dock to each other to form conduits between cells. GJs allow direct metabolic and electrical communication between almost all cell types in the mammalian brain.1 These small molecules involve second messengers (inositol triphosphate (IP3), cyclic nucleotides and Ca2+), neurotransmitters (glutamate, adenosine triphosphate and nitric oxide) and metabolites.2,3 During development, there is extensive regulation of the spatial and temporal expression of connexins, this regulation being implicated in the control of differentiation and growth.4 Injection of antibodies to connexins into blastomeres of amphibian or rodent embryos, or transgenic mice with altered connexin expression results in severe developmental anomalies.5 Differentiation of neurons and astrocytes from embryonal carcinoma cells is accompanied by decreased connexin 43 expression and loss of coupling.6 This differentiation is significantly attenuated if gap junction intercellular communication (GJIC) is blocked.7 Therefore, gap junctions play very important roles in mammalian neural development.1,8 Schizophrenia has generally been considered to be a neurodevelopmental disorder.9 Moreover, an increasing number of reports have demonstrated that gamma‐aminobutyric acid (GABA)ergic interneurons themselves are interconnected through electrical synapses.2 These interconnections, if impaired, may be involved in the aetiology of schizophrenia. Furthermore, changes in GJIC between neurons have been observed in vivo after applying treatment with antipsychotic drugs or withdrawal of amphetamine.2,10 Taking these data together, connexin genes can be considered as potential candidates in the development of schizophrenia. In fact, Cx36 was tested in a large catatonic schizophrenia pedigree, but no association was found.11

A genome‐wide scan for schizophrenia susceptibility loci provided considerable evidence of linkage to chromosome 1q21‐q22.12 The linkage of chromosome 1q with schizophrenia was replicated in several other reports.13,14,15,16 Two connexin genes, Cx50 (GJA8) and Cx40 (GJA5), are located at 1q21. To our knowledge, studies on the genetic association of Cx50 and Cx40 in schizophrenia are lacking. Therefore, we genotyped four polymorphisms for each of these two connexin genes in a matched case‐control sample from Toronto, Canada, and found a significant association between the Cx50 rs989192‐rs4950495 haplotypes and schizophrenia. This association was confirmed in an independent study of Portuguese families with probands with schizophrenia.

SUBJECTS AND METHODS

Subjects

Informed consent was obtained from each participant. For the case–control study, 190 Caucasian patients and 190 controls matched for age and sex (115 men, 75 women; mean (SD) age 43.5 (10.3) years) were analysed. The Structured Clinical Interview for DSM‐IV Axis I Disorders (SCID‐I) was administered by trained research assistants to each patient. The controls were screened for major psychiatric disorders or substance misuse, and excluded if either was detected, either current or in the past. We also investigated 99 Canadian Caucasian families from Toronto (40 triads and 59 diads) with probands diagnosed by the Diagnostic and statistical manual of mental disorders, 4th edition (DSM‐IV) criteria (72 male and 27 female probands), among whom 33 patients overlapped with the case–control study. Another 163 Portuguese trios and nuclear families with probands with DSM‐IV schizophrenia (105 male and 58 female probands, all Caucasians) from both mainland Portugal (Coimbra) and the Azores islands were investigated. The collection of patient samples from Portugal occurred within the framework of a NIMH‐sponsored grant to four of the authors (CP, MP, JLK, MA). The diagnosis was confirmed using the Diagnostic Interview for Genetic Studies. In all samples from Canada and Portugal, a consensus‐based best estimate procedure provided the final decision for the diagnostic classification of the patients. Schizoaffective disorder was included.

Genotyping

Genotypes for the single nucleotide polymorphisms (SNPs) rs989192, rs4950495, rs3766503 and rs1532399 in Cx50, and rs1043806, rs1692141, rs791277 and rs1495956 in Cx40 were assessed by the 5'‐exonuclease activity (TaqMan) assay17 using an automated sequence detection system (ABI Prism 7000; Applied Biosystems, Foster City, California, USA). Figure 11 shows the locations of these polymorphisms.

figure mg47944.f1
Figure 1 Polymorphisms and linkage disequilibrium (LD) in Cx40 and Cx50. Only Caucasian patients with schizophrenia and their controls are shown in the LD diagram. The numbers in the squares are D′ values (×100).

Statistical analyses

Statistical power of the sample size was calculated using the program QUANTO 1.1 for matched case–control and family‐based studies.18 The program HaploView V3.2 was used to test for concordance with Hardy–Weinberg equilibrium, calculate D′ between SNP pairs, and identify linkage disequilibrium (LD) blocks, which were defined with confidence intervals.19 The program UNPHASED V3.0.6 was used for an association analysis for alleles and haplotypes, with the rare haplotype being set to 0.05.20 The p values for each marker and haplotype were adjusted from permutation test with 10 000 simulations. Comparison of genotype frequencies between patients and controls was performed using the χ2 test. Odds ratios (OR) with 95% CI were estimated for the effects of high‐risk alleles.21

RESULTS

Matched case–control study

Our 190 patients and 190 matched controls had 80% power to detect a relative risk as low as 1.89 if the significance was set at 0.05. For Caucasian patients and controls, distributions of genotype frequencies of the eight polymorphisms in Cx40 and Cx50 were all in Hardy–Weinberg equilibrium (p>0.10). The polymorphisms rs1043806 and rs1692141 in Cx40 were tightly linked in one LD block in both patients and controls. However, rs989192 and rs4950495 in Cx50 formed another LD block only in patients with schizophrenia, but not in controls (fig 11).

We did not find any significant differences between patients with schizophrenia and controls (p>0.05) for the allele and genotype frequencies for each individual marker in Cx40 and Cx50 (table 11).). Performing haplotype analyses using two‐marker pairs, we found that the rs989192‐rs4950495 haplotypes were significantly associated with schizophrenia (χ2 = 29.55, p = 0.001; table 22).). Compared with controls, patients had a higher frequency of an A‐C haplotype (χ2 = 7.153, p<0.01, OR = 1.91, 95% CI = 1.17 to 3.13) and low frequencies of a G‐C haplotype (χ2 = 33.44, p<0.001, OR = 5.51, 95% CI = 2.54 to 11.96). Cx50 haplotypes from all four markers showed a significant association with schizophrenia (χ2 = 28.02, p<0.01), but Cx40 did not (χ2 = 0.87, p>0.58).

Table thumbnail
Table 1 Comparison of allele and genotype frequencies of eight markers in Cx40 and Cx50 in 190 Canadian Caucasian patients with schizophrenia and age and sex matched healthy controls
Table thumbnail
Table 2 Analysis of Cx40 and Cx50 haplotypes in 190 Caucasian patients with schizophrenia and 190 matched healthy controls

Family study

The sample size for the family study had 80% power to detect a relative risk as low as 1.67 if significance was set at 0.05. We selected two polymorphisms in Cx50, rs989192 and rs4950495, based on the significant results from our matched case–control study. The distributions of genotype frequencies of these two polymorphisms in both the Canadian and Portuguese families were in Hardy–Weinberg equilibrium (p>0.05). These two markers also formed one LD block in the families. In the Portuguese families, the rs989192A allele was transmitted more often to patients with schizophrenia (transmitted:untransmitted (T:U)  = 62:36.54, χ2 = 5.248, p = 0.022, OR = 1.79, 95% CI = 1.17 to 2.75). There was also a difference between the transmitted and untransmitted ratio of the rs4950495C allele (T:U = 76:40.39, χ2 = 8.823, p = 0.003, OR = 2.03, 95% CI = 1.36 to 3.04) (fig 22).). Haplotype analysis for these two markers showed a significant association between Cx50 and schizophrenia (χ2 = 8.707, p = 0.037). The A‐C haplotype was transmitted more often (χ2 = 6.096, p = 0.014, OR = 1.9, 95% CI = 1.23 to 2.95) and the G‐T haplotype less often (χ2 = 6.218, p = 0.013, OR = 1.77, 95% CI = 1.2 to 2.6) to patients with schizophrenia (table 33).). In the Canadian families, we found a trend towards significance for rs989192 (A allele T:U = 30:16.92, χ2 = 3.236, p = 0.072) (fig 22),), but no association between the Cx50 haplotype and schizophrenia (χ2 = 2.671, p = 0.263) (table 33).). Combining the data from the Canadian and Portuguese families together, we found that the rs989192A and rs4950495C alleles were transmitted significantly more often to patients with schizophrenia (T:U = 92:53.73, χ2 = 8.271, p = 0.004, OR = 1.82, 95% CI = 1.28 to 2.59 and T:U = 114:69.1, χ2 = 9.316, p = 0.002, OR = 1.77, 95% CI = 1.29 to 2.44, respectively) (fig 22).). We also found an association between the Cx50 rs989192‐rs4950495 haplotype and schizophrenia (χ2 = 10.315, p = 0.012). The A‐C haplotype was transmitted 90.98 times and not transmitted 52.61 times (χ2 = 8.433, p = 0.004, OR = 1.83, 95% CI = 1.28 to 2.62) and the G‐T haplotype was transmitted less often to patients (T:U = 608:649.7, χ2 = 7.876, p = 0.005, OR = 1.67, 95% CI = 1.22 to 2.28) (table 33).). No significant differences in maternal and paternal transmission were detected in our family studies (data not shown).

figure mg47944.f2
Figure 2 Transmission disequilibrium test (TDT) with Cx50 rs989192 and rs4950495 in Canadian and Portuguese families.
Table thumbnail
Table 3 Transmission disequilibrium test with Cx50 rs989192‐rs4950495 haplotypes in 99 Canadian Caucasian and 163 Portuguese families with schizophrenia.

DISCUSSION

For the first time, to our knowledge, an association between the Cx genes on chromosome 1q21 and schizophrenia has been investigated, and we found a significant association between the Cx50 rs989192‐rs4950495 haplotype and schizophrenia. Compared with controls, patients with schizophrenia had a higher frequency of an A‐C haplotype (p = 0.007). In the Canadian and Portuguese families, the A‐C haplotype was transmitted more often to patients with schizophrenia (p = 0.004). No association of Cx40 with schizophrenia was found for allele, genotype or haplotype analyses. Our matched case–control and family studies indicate that Cx50, but not Cx40, may play a role in the genetic susceptibility to schizophrenia.

The rs989192A‐rs4950495C haplotype shows a significant association with schizophrenia. Whether this intronic haplotype can affect Cx50 gene expression or link to other functional markers is unknown. Although other rs989192‐rs4950495 haplotypes show different results between our matched case–control and family studies, we found that the haplotype frequencies were similar in Canadian patients and in the frequency of transmission patients in families, particularly compared with the frequencies in controls ((tablestables 2 and 33).). Population stratification in this study can be excluded to a certain extent because the patients and controls were matched for sex, age and ethnicity, and this association is replicated in the family study. However, we should consider effects from genotyping error and rare haplotypes,22 although the rare haplotype was set to 0.05. The results from the Canadian Caucasian families showed no association with schizophrenia although there was a trend towards significance with rs989192. The result may mainly be due to the small sample size of only 40 triads and 59 diads.

KEY POINTS

  • The connexin genes Cx40 and Cx50 are considered candidates for schizophrenia because of their functions in the brain and their chromosomal location at 1q21, a region strongly linked with schizophrenia.
  • We found that the Cx50 rs989192‐rs4950495 haplotype was associated with schizophrenia in matched case–control samples from Toronto, Canada, and this result was replicated in a Portuguese family study.
  • No association between Cx40 and schizophrenia was found for allele, genotype or haplotype analyses.

Compared with other connexins, Cx50 channels are among the most CO2‐sensitive connexin channels.23,24 At negative voltages, Cx50 currents do not deactivate completely, in contrast, to Cx46 currents, which do.25 These characteristics of Cx50 may indicate its physiological functions and pathology in brain. It has been reported that Cx50 is not present in rat brain tissues,26 other than as some unexplained non‐punctate labelling in the choroid plexus.27 However, from the Allen brain atlas (www.brain‐map.org) we can see that Cx50 is expressed in the brain, including the hippocampus, lateral ventricle, prefrontal cortex, amygdala and cerebellum. In addition, using real‐time PCR, we have detected Cx50 expression in the brain of rats undergoing a conditioned avoidance response test, and in yoked and untrained controls (data not shown). Cx40 is expressed in the conduction system and atria,28 and in neurons,29 astrocytes30 and neuronal precursors.31 The Cx40 gene is close to Cx50, at a distance of 129.5 kb. Although no significant association of Cx40 with schizophrenia was found in the present study, Cx40 markers can be considered agood null controls.32

Further and more detailed studies are necessary to elucidate the role of Cx50 in schizophrenia. We also have checked a marker in the promoter region of Cx50, rs7536277; unfortunately, no variant was found in our samples. Cx expression and GJIC in astrocytes and in neurons are targets for a number of neurotransmitters, growth factors, peptides, cytokines and endogenous bioactive lipids.2 Although few reports have been published, Cx50 may interact with GABA receptors (GABRA5, GABRB3 and GABRG2; http://string.embl.de/). Insights into the interactions between connexins and neurotransmitters will advance our knowledge about connexins in the aetiology of schizophrenia.

ACKNOWLEDGEMENTS

This study was supported by the National Alliance for Research on Schizophrenia and Depression (NARSAD) (X. N.).

Abbreviations

Cx40 - connexin 40

Cx50 - connexin 50

GABA - gamma‐aminobutyric acid

GJ - gap junction

GJIC - inositol triphosphate

LD - linkage disequilibrium

SCID‐I - Structured Clinical Interview for DSM‐IV Axis I Disorders

SNP - single nucleotide polymorphism

T:U - transmitted:untransmitted ratio

Footnotes

Competing interests: None declared.

References

1. Sohl G, Maxeiner S, Willecke K. Expression and functions of neuronal gap junctions. Nat Rev Neurosci 2005. 6191–200.200 [PubMed]
2. Rouach N, Avignone E, Meme W, Koulakoff A, Venance L, Blomstrand F, Giaume C. Gap junctions and connexin expression in the normal and pathological central nervous system. Biol Cell 2002. 94457–475.475 [PubMed]
3. Kirchhoff F, Dringen R, Giaume C. Pathways of neuron‐astrocyte interactions and their possible role in neuroprotection. Eur Arch Psychiatry Clin Neurosci 2001. 251159–69 [PubMed]
4. Naus C C, Bani‐Yaghoub M. Gap junctional communication in the developing central nervous system. Cell Biol Int 1998. 22751–763.763 [PubMed]
5. Naus C C. Gap junctions and tumour progression. Can J Physiol Pharmacol 2002. 80136–141.141 [PubMed]
6. Belliveau D J, Naus C C. Cellular localization of gap junction mRNAs in developing rat brain. Dev Neurosci 1995. 1781–96.96 [PubMed]
7. Bani‐Yaghoub M, Bechberger J F, Underhill T M, Naus C C. The effects of gap junction blockage on neuronal differentiation of human NTera2/clone D1 cells. Exp Neurol 1999. 15616–32.32 [PubMed]
8. Wei C J, Xu X, Lo C W. Connexins and cell signaling in development and disease. Annu Rev Cell Dev Biol 2004. 20811–838.838 [PubMed]
9. Waddington J L, Lane A, Scully P J, Larkin C, O'Callaghan E. Neurodevelopmental and neuroprogressive processes in schizophrenia. Antithetical or complementary, over a lifetime trajectory of disease? Psychiatr Clin North Am 1998. 21123–149.149 [PubMed]
10. O'Donnell P, Grace A A. Different effects of subchronic clozapine and haloperidol on dye‐coupling between neurons in the rat striatal complex. Neuroscience 1995. 6763–767.767 [PubMed]
11. Meyer J, Mai M, Ortega G, Mossner R, Lesch K P. Mutational analysis of the connexin 36 gene (CX36) and exclusion of the coding sequence as a candidate region for catatonic schizophrenia in a large pedigree. Schizophr Res 2002. 5887–91.91 [PubMed]
12. Brzustowicz L M, Hodgkinson K A, Chow E W, Honer W G, Bassett A S. Location of a major susceptibility locus for familial schizophrenia on chromosome 1q21‐q22. Science 2000. 288678–682.682 [PubMed]
13. Hwu H G, Liu C M, Fann C S, Ou‐Yang W C, Lee S F. Linkage of schizophrenia with chromosome 1q loci in Taiwanese families. Mol Psychiatry 2003. 8445–452.452 [PubMed]
14. Gasperoni T L, Ekelund J, Huttunen M, Palmer C G, Tuulio‐Henriksson A, Lonnqvist J, Kaprio J, Peltonen L, Cannon T D. Genetic linkage and association between chromosome 1q and working memory function in schizophrenia. Am J Med Genet B Neuropsychiatr Genet 2003. 1168–16.16 [PubMed]
15. Ekelund J, Hovatta I, Parker A, Paunio T, Varilo T, Martin R, Suhonen J, Ellonen P, Chan G, Sinsheimer J S, Sobel E, Juvonen H, Arajarvi R, Partonen T, Suvisaari J, Lonnqvist J, Meyer J, Peltonen L. Chromosome 1 loci in Finnish schizophrenia families. Hum Mol Genet 2001. 101611–1617.1617 [PubMed]
16. Cai G, Wu X, Li T, Collier D A, Liu X, Feng B, Deng H, Tong D, Li J, Ou J. Linkage analysis of susceptibility genes for familial schizophrenia on chromosome 1 in Chinese population. (In Chinese. ) Zhonghua Yi Xue Yi Chuan Xue Za Zhi2002. 19491–494.494 [PubMed]
17. Heid C A, Stevens J, Livak K J, Williams P M. Real time quantitative PCR. Genome Res 1996. 6986–994.994 [PubMed]
18. Gauderman W J, Morrison J M. QUANTO 1. 1: A computer program for power and sample size calculations for genetic‐epidemiology studies, http://hydra.usc.edu/gxe 2006
19. Barrett J C, Fry B, Maller J, Daly M J. Haploview: analysis and visualization of LD and haplotype maps. Bioinformatics 2005. 21263–265.265 [PubMed]
20. Dudbridge F. Pedigree disequilibrium tests for multilocus haplotypes. Genet Epidemiol 2003. 25115–121.121 [PubMed]
21. Bland J M, Altman D G. Statistics notes. The odds ratio. BMJ 2000. 3201468 [PMC free article] [PubMed]
22. Moskvina V, Craddock N, Holmans P, Owen M J, O'Donovan M C. Effects of differential genotyping error rate on the type I error probability of case‐control studies. Hum Hered 2006. 6155–64.64 [PubMed]
23. Xu X, Berthoud V M, Beyer E C, Ebihara L. Functional role of the carboxyl terminal domain of human connexin 50 in gap junctional channels. J Membr Biol 2002. 186101–112.112 [PMC free article] [PubMed]
24. Beahm D L, Hall J E. Hemichannel and junctional properties of connexin 50. Biophys J 2002. 822016–2031.2031 [PubMed]
25. Srinivas M, Kronengold J, Bukauskas F F, Bargiello T A, Verselis V K. Correlative studies of gating in Cx46 and Cx50 hemichannels and gap junction channels. Biophys J 2005. 881725–1739.1739 [PubMed]
26. Vis J C, Nicholson L F, Faull R L, Evans W H, Severs N J, Green C R. Connexin expression in Huntington's diseased human brain. Cell Biol Int 1998. 22837–847.847 [PubMed]
27. Mack A, Neuhaus J, Wolburg H. Relationship between orthogonal arrays of particles and tight junctions as demonstrated in cells of the ventricular wall of the rat brain. Cell Tissue Res 1987. 248619–625.625 [PubMed]
28. Severs N J. Gap junction remodeling in heart failure. J Card Fail 2002. 8(6 Suppl)S293–S299.S299 [PubMed]
29. Chang Q, Gonzalez M, Pinter M J, Balice‐Gordon R J. Gap junctional coupling and patterns of connexin expression among neonatal rat lumbar spinal motor neurons. J Neurosci 1999. 1910813–10828.10828 [PubMed]
30. Dermietzel R, Gao Y, Scemes E, Vieira D, Urban M, Kremer M, Bennett M V, Spray D C. Connexin43 null mice reveal that astrocytes express multiple connexins. Brain Res Brain Res Rev 2000. 3245–56.56 [PubMed]
31. Rozental R, Morales M, Mehler M F, Urban M, Kremer M, Dermietzel R, Kessler J A, Spray D C. Changes in the properties of gap junctions during neuronal differentiation of hippocampal progenitor cells. J Neurosci 1998. 181753–1762.1762 [PubMed]
32. Balding D J. A tutorial on statistical methods for population association studies. Nat Rev Genet 2006. 7781–791.791 [PubMed]

Articles from Journal of Medical Genetics are provided here courtesy of BMJ Group