Skip to main content

Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

Schizophrenia risk from complex variation of complement component 4

An Author Correction to this article was published on 15 December 2021

This article has been updated

Abstract

Schizophrenia is a heritable brain illness with unknown pathogenic mechanisms. Schizophrenia’s strongest genetic association at a population level involves variation in the major histocompatibility complex (MHC) locus, but the genes and molecular mechanisms accounting for this have been challenging to identify. Here we show that this association arises in part from many structurally diverse alleles of the complement component 4 (C4) genes. We found that these alleles generated widely varying levels of C4A and C4B expression in the brain, with each common C4 allele associating with schizophrenia in proportion to its tendency to generate greater expression of C4A. Human C4 protein localized to neuronal synapses, dendrites, axons, and cell bodies. In mice, C4 mediated synapse elimination during postnatal development. These results implicate excessive complement activity in the development of schizophrenia and may help explain the reduced numbers of synapses in the brains of individuals with schizophrenia.

This is a preview of subscription content, access via your institution

Access options

Rent or buy this article

Prices vary by article type

from$1.95

to$39.95

Prices may be subject to local taxes which are calculated during checkout

Figure 1: Structural variation of the complement component 4 (C4) gene.
Figure 2: Haplotypes formed by C4 structures and SNPs.
Figure 3: Brain RNA expression of C4A and C4B in relation to copy numbers of C4A, C4B, and the C4–HERV.
Figure 4: Association of schizophrenia to C4 and the extended MHC locus.
Figure 5: C4 structures, C4A expression, and schizophrenia risk.
Figure 6: C4 protein at neuronal cell bodies, processes and synapses.
Figure 7: C4 in retinogeniculate synaptic refinement.

Similar content being viewed by others

Change history

References

  1. Cannon, T. D. et al. Cortex mapping reveals regionally specific patterns of genetic and disease-specific gray-matter deficits in twins discordant for schizophrenia. Proc. Natl Acad. Sci. USA 99, 3228–3233 (2002)

    Article  ADS  CAS  PubMed  PubMed Central  Google Scholar 

  2. Cannon, T. D. et al. Progressive reduction in cortical thickness as psychosis develops: a multisite longitudinal neuroimaging study of youth at elevated clinical risk. Biol. Psychiatry 77, 147–157 (2015)

    Article  PubMed  Google Scholar 

  3. Garey, L. J. et al. Reduced dendritic spine density on cerebral cortical pyramidal neurons in schizophrenia. J. Neurol. Neurosurg. Psychiatry 65, 446–453 (1998)

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Glantz, L. A. & Lewis, D. A. Decreased dendritic spine density on prefrontal cortical pyramidal neurons in schizophrenia. Arch. Gen. Psychiatry 57, 65–73 (2000)

    Article  CAS  PubMed  Google Scholar 

  5. Glausier, J. R. & Lewis, D. A. Dendritic spine pathology in schizophrenia. Neuroscience 251, 90–107 (2013)

    Article  CAS  PubMed  Google Scholar 

  6. Schizophrenia Working Group of the Psychiatric Genomics Consortium. Biological insights from 108 schizophrenia-associated genetic loci. Nature 511, 421–427 (2014)

  7. Shi, J. et al. Common variants on chromosome 6p22.1 are associated with schizophrenia. Nature 460, 753–757 (2009)

    Article  ADS  CAS  PubMed  PubMed Central  Google Scholar 

  8. Stefansson, H. et al. Common variants conferring risk of schizophrenia. Nature 460, 744–747 (2009)

    Article  ADS  CAS  PubMed  PubMed Central  Google Scholar 

  9. International Schizophrenia Consortium et al. Common polygenic variation contributes to risk of schizophrenia and bipolar disorder. Nature 460, 748–752 (2009)

    Article  CAS  Google Scholar 

  10. Schizophrenia Psychiatric Genome-Wide Association Study Consortium. Genome-wide association study identifies five new schizophrenia loci. Nature Genet . 43, 969–976 (2011)

  11. Howson, J. M., Walker, N. M., Clayton, D. & Todd, J. A. Confirmation of HLA class II independent type 1 diabetes associations in the major histocompatibility complex including HLA-B and HLA-A. Diabetes Obes. Metab. 11 (Suppl 1), 31–45 (2009)

    Article  PubMed  PubMed Central  Google Scholar 

  12. Raychaudhuri, S. et al. Five amino acids in three HLA proteins explain most of the association between MHC and seropositive rheumatoid arthritis. Nature Genet. 44, 291–296 (2012)

    Article  CAS  PubMed  Google Scholar 

  13. Escudero-Esparza, A., Kalchishkova, N., Kurbasic, E., Jiang, W. G. & Blom, A. M. The novel complement inhibitor human CUB and Sushi multiple domains 1 (CSMD1) protein promotes factor I-mediated degradation of C4b and C3b and inhibits the membrane attack complex assembly. FASEB J . 27, 5083–5093 (2013)

    Article  CAS  PubMed  Google Scholar 

  14. Carroll, M. C., Campbell, R. D., Bentley, D. R. & Porter, R. R. A molecular map of the human major histocompatibility complex class III region linking complement genes C4, C2 and factor B. Nature 307, 237–241 (1984)

    Article  ADS  CAS  PubMed  Google Scholar 

  15. Carroll, M. C., Belt, T., Palsdottir, A. & Porter, R. R. Structure and organization of the C4 genes. Phil. Trans. R. Soc. Lond. B 306, 379–388 (1984)

    Article  ADS  CAS  Google Scholar 

  16. Dangel, A. W. et al. The dichotomous size variation of human complement C4 genes is mediated by a novel family of endogenous retroviruses, which also establishes species-specific genomic patterns among Old World primates. Immunogenetics 40, 425–436 (1994)

    Article  CAS  PubMed  Google Scholar 

  17. Horton, R. et al. Variation analysis and gene annotation of eight MHC haplotypes: the MHC Haplotype Project. Immunogenetics 60, 1–18 (2008)

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Bánlaki, Z., Doleschall, M., Rajczy, K., Fust, G. & Szilagyi, A. Fine-tuned characterization of RCCX copy number variants and their relationship with extended MHC haplotypes. Genes Immun. 13, 530–535 (2012)

    Article  PubMed  CAS  Google Scholar 

  19. Law, S. K., Dodds, A. W. & Porter, R. R. A comparison of the properties of two classes, C4A and C4B, of the human complement component C4. EMBO J. 3, 1819–1823 (1984)

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Isenman, D. E. & Young, J. R. The molecular basis for the difference in immune hemolysis activity of the Chido and Rodgers isotypes of human complement component C4. J. Immunol. 132, 3019–3027 (1984)

    CAS  PubMed  Google Scholar 

  21. Illarionova, A. E., Vinogradova, T. V. & Sverdlov, E. D. Only those genes of the KIAA1245 gene subfamily that contain HERV(K) LTRs in their introns are transcriptionally active. Virology 358, 39–47 (2007)

    Article  CAS  PubMed  Google Scholar 

  22. Nakamura, A., Okazaki, Y., Sugimoto, J., Oda, T. & Jinno, Y. Human endogenous retroviruses with transcriptional potential in the brain. J. Hum. Genet. 48, 575–581 (2003)

    Article  CAS  PubMed  Google Scholar 

  23. Suntsova, M. et al. Human-specific endogenous retroviral insert serves as an enhancer for the schizophrenia-linked gene PRODH. Proc. Natl Acad. Sci. USA 110, 19472–19477 (2013)

    Article  ADS  CAS  PubMed  PubMed Central  Google Scholar 

  24. Yang, Y. et al. Diversity in intrinsic strengths of the human complement system: serum C4 protein concentrations correlate with C4 gene size and polygenic variations, hemolytic activities, and body mass index. J. Immunol. 171, 2734–2745 (2003)

    Article  CAS  PubMed  Google Scholar 

  25. Browning, S. R. & Browning, B. L. Rapid and accurate haplotype phasing and missing-data inference for whole-genome association studies by use of localized haplotype clustering. Am. J. Hum. Genet. 81, 1084–1097 (2007)

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  26. Iossifov, I. et al. The contribution of de novo coding mutations to autism spectrum disorder. Nature 515, 216–221 (2014)

    Article  ADS  CAS  PubMed  PubMed Central  Google Scholar 

  27. Mayilyan, K. R., Arnold, J. N., Presanis, J. S., Soghoyan, A. F. & Sim, R. B. Increased complement classical and mannan-binding lectin pathway activities in schizophrenia. Neurosci. Lett. 404, 336–341 (2006)

    Article  CAS  PubMed  Google Scholar 

  28. Hakobyan, S., Boyajyan, A. & Sim, R. B. Classical pathway complement activity in schizophrenia. Neurosci. Lett. 374, 35–37 (2005)

    Article  CAS  PubMed  Google Scholar 

  29. Stevens, B. et al. The classical complement cascade mediates CNS synapse elimination. Cell 131, 1164–1178 (2007)

    Article  CAS  PubMed  Google Scholar 

  30. Schafer, D. P. et al. Microglia sculpt postnatal neural circuits in an activity and complement-dependent manner. Neuron 74, 691–705 (2012)

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Bialas, A. R. & Stevens, B. TGF-β signaling regulates neuronal C1q expression and developmental synaptic refinement. Nature Neurosci. 16, 1773–1782 (2013)

    Article  CAS  PubMed  Google Scholar 

  32. Kaiser, T. & Feng, G. Modeling psychiatric disorders for developing effective treatments. Nature Med. 21, 979–988 (2015)

    Article  CAS  PubMed  Google Scholar 

  33. Shatz, C. J. & Kirkwood, P. A. Prenatal development of functional connections in the cat’s retinogeniculate pathway. J. Neurosci. 4, 1378–1397 (1984)

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Sretavan, D. W. & Shatz, C. J. Prenatal development of retinal ganglion cell axons: segregation into eye-specific layers within the cat’s lateral geniculate nucleus. J. Neurosci. 6, 234–251 (1986)

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Chen, C. & Regehr, W. G. Developmental remodeling of the retinogeniculate synapse. Neuron 28, 955–966 (2000)

    Article  CAS  PubMed  Google Scholar 

  36. Fischer, M. B. et al. Regulation of the B cell response to T-dependent antigens by classical pathway complement. J. Immunol. 157, 549–556 (1996)

    CAS  PubMed  Google Scholar 

  37. Huttenlocher, P. R. & Dabholkar, A. S. Regional differences in synaptogenesis in human cerebral cortex. J. Comp. Neurol. 387, 167–178 (1997)

    Article  CAS  PubMed  Google Scholar 

  38. Huttenlocher, P. R. Synaptic density in human frontal cortex—developmental changes and effects of aging. Brain Res. 163, 195–205 (1979)

    Article  CAS  PubMed  Google Scholar 

  39. Petanjek, Z. et al. Extraordinary neoteny of synaptic spines in the human prefrontal cortex. Proc. Natl Acad. Sci. USA 108, 13281–13286 (2011)

    Article  ADS  CAS  PubMed  PubMed Central  Google Scholar 

  40. Buckner, R. L. & Krienen, F. M. The evolution of distributed association networks in the human brain. Trends Cogn. Sci. 17, 648–665 (2013)

    Article  PubMed  Google Scholar 

  41. Feinberg, I. Schizophrenia: caused by a fault in programmed synaptic elimination during adolescence? J. Psychiatr. Res. 17, 319–334 (1982–1983)

    Article  PubMed  Google Scholar 

  42. Kirov, G. et al. De novo CNV analysis implicates specific abnormalities of postsynaptic signalling complexes in the pathogenesis of schizophrenia. Mol. Psychiatry 17, 142–153 (2012)

    Article  CAS  PubMed  Google Scholar 

  43. Fromer, M. et al. De novo mutations in schizophrenia implicate synaptic networks. Nature 506, 179–184 (2014)

    Article  ADS  CAS  PubMed  PubMed Central  Google Scholar 

  44. Purcell, S. M. et al. A polygenic burden of rare disruptive mutations in schizophrenia. Nature 506, 185–190 (2014)

    Article  ADS  CAS  PubMed  PubMed Central  Google Scholar 

  45. Datwani, A. et al. Classical MHCI molecules regulate retinogeniculate refinement and limit ocular dominance plasticity. Neuron 64, 463–470 (2009)

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. Lee, H. et al. Synapse elimination and learning rules co-regulated by MHC class I H2-Db. Nature 509, 195–200 (2014)

    Article  ADS  CAS  PubMed  PubMed Central  Google Scholar 

  47. van den Elsen, J. M. et al. X-ray crystal structure of the C4d fragment of human complement component C4. J. Mol. Biol. 322, 1103–1115 (2002)

    Article  CAS  PubMed  Google Scholar 

  48. Dodds, A. W., Ren, X. D., Willis, A. C. & Law, S. K. The reaction mechanism of the internal thioester in the human complement component C4. Nature 379, 177–179 (1996)

    Article  ADS  CAS  PubMed  Google Scholar 

  49. Handsaker, R. E. et al. Large multiallelic copy number variations in humans. Nature Genet. 47, 296–303 (2015)

    Article  CAS  PubMed  Google Scholar 

  50. Torborg, C. L. & Feller, M. B. Unbiased analysis of bulk axonal segregation patterns. J. Neurosci. Methods 135, 17–26 (2004)

    Article  CAS  PubMed  Google Scholar 

  51. Fernando, M. M. et al. Assessment of complement C4 gene copy number using the paralog ratio test. Hum. Mutat. 31, 866–874 (2010)

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Rudduck, C., Beckman, L., Franzen, G., Jacobsson, L. & Lindstrom, L. Complement factor C4 in schizophrenia. Hum. Hered. 35, 223–226 (1985)

    Article  CAS  PubMed  Google Scholar 

  53. Schroers, R. et al. Investigation of complement C4B deficiency in schizophrenia. Hum. Hered. 47, 279–282 (1997)

    Article  CAS  PubMed  Google Scholar 

  54. Mayilyan, K. R., Dodds, A. W., Boyajyan, A. S., Soghoyan, A. F. & Sim, R. B. Complement C4B protein in schizophrenia. World J. Biol. Psychiatry 9, 225–230 (2008)

    Article  PubMed  Google Scholar 

  55. Jia, X. et al. Imputing amino acid polymorphisms in human leukocyte antigens. PLoS ONE 8, e64683 (2013)

    Article  ADS  CAS  PubMed  PubMed Central  Google Scholar 

  56. Nonaka, M., Nakayama, K., Yeul, Y. D. & Takahashi, M. Complete nucleotide and derived amino acid sequences of sex-limited protein (Slp), nonfunctional isotype of the fourth component of mouse complement (C4). J. Immunol. 136, 2989–2993 (1986)

    CAS  PubMed  Google Scholar 

Download references

Acknowledgements

The authors would like to remember the late T. Stanley with appreciation and express their gratitude for his support. We thank S. Hyman, E. Lander, C. Bargmann, and C. Patil for conversations about the project and comments on drafts of the manuscript; M. Webster for expert advice on immunohistochemistry; B. Browning for expert advice on imputation; the Stanley Medical Research Institute Brain Collection and the NHGRI Gene and Tissue Expression (GTEx) Project for access to RNA and tissue samples; C. Emba for assistance with experiments; and C. Usher for contributions to manuscript figures. This work was supported by R01 HG 006855 (to S.A.M.), by the Stanley Center for Psychiatric Research (to S.A.M. and B.S.), by U01 MH105641 (to S.A.M.), by R01 MH077139 (to the PGC), and by T32 GM007753 (to A.S. and M.B.).

Author information

Authors and Affiliations

Authors

Consortia

Contributions

S.A.M. and A.S. conceived the genetic studies. A.S. performed the laboratory experiments and computational analyses to understand the molecular and population genetics of the C4 locus (Figs 1 and 2). A.S., K.T., N.K., and V.V.D. analysed C4 expression variation in human brain (Figs 3 and 5b, d). G.G., R.E.H., and S.A.R. contributed to genetic analyses. A.S. and A.D. did the imputation and association analysis (Figs 4 and 5a, c). M.J.D. provided advice on the association analyses. Investigators in the Schizophrenia Working Group of the Psychiatric Genomics Consortium collected and phenotyped cohorts and contributed genotype data for analysis. B.S. and M.C.C. contributed expertise and reagents for experiments described in Fig. 6 and 7. H.d.R and T.R.H. performed the C4 immunocytochemistry and immunohistochemistry experiments respectively, with advice from A.R.B. (Fig. 6). A.R.B. and J.P. analysed the role of C4 in synaptic refinement in the mouse visual system (Fig. 7). M.B. analysed C4 expression in mice. S.A.M and A.S. wrote the manuscript with contributions from all authors.

Corresponding author

Correspondence to Steven A. McCarroll.

Ethics declarations

Competing interests

The authors declare no competing financial interests.

Additional information

Lists of participants and their affiliations appear in the Supplementary Information.

Extended data figures and tables

Extended Data Figure 1 Association of schizophrenia to common variants in the MHC locus in individual case-control cohorts, and schematic of the repeat module containing C4.

af, Data for several schizophrenia case-control cohorts that were genome-scanned before we began this work (ad) exhibits peaks of association near chr6: 32 Mb (blue vertical line) on the human genome reference sequence (GRCh37/hg19). Note that association patterns vary from cohort to cohort, reflecting statistical sampling fluctuations and potentially fluctuations in allele frequencies of the (unknown) causal variants in different cohorts. Cohorts such as in b, e and f suggest the existence of effects at multiple loci within the MHC region. Even in the cohorts with simpler peaks (a, c, d), the pattern of association across the individual SNPs at chr6: 32 Mb does not correspond to the LD around any known variant. This motivated the focus in the current work on cryptic genetic influences in this region that could cause unconventional association signals that do not resemble the LD patterns of individual variants. g, A complex form of genome structural variation resides near chr6: 32 Mb. Shown here are three of the known alternative structural forms of this genomic region. The most prominent feature of this structural variation is the tandem duplication of a genomic segment that contains a C4 gene, 3′ fragments of the STK19 and TNXB genes, and a pseudogenized copy of the CYP21A2 gene. This cassette is present in 1–3 copies on the three alleles depicted above; the boundaries below each haplotype demarcate the sequence that is duplicated. Haplotypes with multiple copies of this module (middle and bottom) contain multiple functional copies of C4, whereas the additional gene fragments or copies denoted STK19P, CYP21A2P, and TNXA are typically pseudogenized. Rare haplotypes with a gain or loss of intact CYP21A2 have also been observed18. Although C4A and C4B contain multiple sequence variants, they are defined based on the differences encoded by exon 26, which determine the relative affinities of C4A and C4B for distinct molecular targets19,20 (Fig. 1). Many additional forms of this locus appear to have arisen by non-allelic homologous recombination and gene conversion (ref. 18 and Fig. 1).

Extended Data Figure 2 Schematic of strategy for identifying the segregating structural forms of the C4 locus.

a, Molecular assays for measuring copy number of the key, variable C4 structural features—the length polymorphism (HERV insertion) that distinguishes the long (L) from the short (S) genomic form of C4, and the C4A/C4B isotypic difference. Each primer–probe–primer assay is represented with the combination of arrows (primers) and asterisk (probe) in its approximate genomic location (though not to scale). b, Measurement of copy number of C4 gene types in the genomes of 162 individuals (from HapMap CEU sample). The absolute, integer copy number of each C4 gene type in each genome is precisely inferred from the resulting data. To ensure high accuracy, the data are further evaluated for a checksum relationship (A + B = L + S) and for concordance with earlier data from Southern blotting of 89 of the same HapMap individuals51. c, To measure the copy number of compound structural forms of C4 (involving combinations of L/S and A/B), we perform long-range PCR followed by quantitative measurement of the A/B isotype-distinguishing sequences in droplets. d, Analysis of transmissions in father–mother–offspring trios enables inference of the C4 gene contents of individual copies (alleles) of chromosome 6. Three example trios are shown in this schematic. e, Examples of the inferred structural forms of the C4 locus (more shown in Fig. 1c). For the common C4 structures (AL–BL, AL–BS, AL–AL, and BS), genomic order of the C4 gene copies is known from earlier assemblies of sequence contigs in individuals homozygous for MHC haplotypes due to consanguinity17 and other molecular analyses of the C4 locus18. For the rarer C4 structures, the genomic order of C4 gene copies is hypothesized or provisional.

Extended Data Figure 3 Linkage disequilibrium relationships (r2) of MHC SNPs to forms of C4 structural variation.

a, b, Correlations of SNPs in the MHC locus with copy number of C4 gene types (a) and larger-scale structural forms (haplotypes) (b) of the C4 locus. Dashed, vertical lines indicate the genomic location of the C4 locus. C4 structural forms show only partial correlation (r2) to the allelic states of nearby SNPs, reflecting the relationship shown in Fig. 2, in which a structural form of the C4 locus often segregates on multiple different SNP haplotypes.

Extended Data Figure 4 RNA expression of C4A and C4B in relation to copy number of C4A, C4B, and the C4–HERV (long form of C4), in eight panels of post-mortem brain tissue.

Copy number of C4 structural features was measured by ddPCR; RNA expression levels were measured by RT-ddPCR. ae, Data for tissues from the Stanley Medical Research Institute (SMRI) Array Consortium consisting of anterior cingulate cortex (a), cerebellum (b), corpus callosum (c), orbital frontal cortex (d), and parietal cortex (e). f, Data for the frontal cortex samples from the NHGRI Genes and Tissues Expression (GTEx) Project. g, h, Data for tissues from the SMRI Neuropathology Consortium (anterior cingulate cortex and cerebellum, respectively). These data were then used to inform (by linear regression) the derivation of a linear model for predicting each individual’s RNA expression of C4A and C4B as a function of the numbers of copies of AL, BL, AS, and BS. The derivation of this model, and the regression coefficients induced, are described in Supplementary Methods. In the rightmost plot of each panel, expression of C4A (per genomic copy) is normalized to expression of C4B (per genomic copy) to more specifically visualize the effect of the C4–HERV by controlling for genomic copy number and for any trans-acting influences shared by C4A and C4B; the inferred regression coefficients (Supplementary Methods) suggest that the observed effect is mostly due to increased expression of C4A.

Extended Data Figure 5 Detailed analysis of the association of schizophrenia to genetic variation at and around C4, in data from 28,799 schizophrenia cases and 35,986 controls.

(Psychiatric Genomics Consortium, ref. 6.) SCZ, schizophrenia; β, estimated effect size per copy of the genomic feature or allele indicated; SE, standard error. Detailed association analyses of HLA alleles are in Extended Data Figs 6 and 7. The single asterisk (*) indicates that we specifically tested C4B-null status because a 1985 study52 reported an analysis of 165 schizophrenia patients and 330 controls in which rare C4B-null status associated with elevated risk of schizophrenia, though two subsequent studies53,54 found no association of schizophrenia to C4B-null genotype. We sought to evaluate this using the large data set in this study, finding no association to C4B-null status. The double asterisk (**) indicates total copy number of C4 is also strongly correlated to copy number of the CYP21A2P pseudogene, which is present on duplicated copies of the sequence shown in Extended Data Fig. 1g.

Extended Data Figure 6 Evaluation of the association of schizophrenia with HLA alleles and coding-sequence polymorphisms.

ae, Associations to HLA alleles and coding-sequence polymorphisms are shown in black; to provide the context of levels of association to nearby SNPs, associations to other SNPs are shown in grey. The series of conditional analyses shown in be parallels the analyses in Fig. 4. Further detail on the most strongly associating HLA alleles (including conditional association analysis) is provided in Extended Data Fig. 7.

Extended Data Figure 7 Detailed association analysis for the most strongly associating classical HLA alleles.

The most strongly associating HLA loci were HLA-B (in primary analyses, Fig. 4a and Extended Data Fig. 6a) and HLA-DRB1 and HLA-DQB1 (in analyses controlling for the signal defined by rs13194504, Fig. 4c and Extended Data Fig. 6b). At these loci, the most strongly associating classical HLA alleles were HLA-B*0801, HLA-DRB1*0301, and HLA-DQB*02, respectively. These HLA alleles are all in strong but partial LD with C4 BS, the most protective of the C4 alleles; they are also in partial LD with the low-risk allele at rs13194505, representing the distinct signal several megabases to the left (Fig. 4). In joint analyses with each of these HLA alleles, genetically predicted C4A expression and rs13194505 continued to associate strongly with schizophrenia, while the HLA alleles did not. In further joint analyses with rs13194504 and genetically predicted C4A expression, 0 of 2,514 tested HLA SNP, amino acid and classical-allele polymorphisms (from ref. 55, including all variants with minor allele frequency (MAF) >0.005) associated with schizophrenia as strongly as rs13194504 or predicted C4A expression did.

Extended Data Figure 8 Expression of C4A RNA in brain tissue (five brain regions) from 35 schizophrenia cases and 70 non-schizophrenia controls, from the Stanley Medical Research Institute Array Consortium.

C4A RNA expression levels were measured by ddPCR. P values are derived from Mann–Whitney U-test.

Extended Data Figure 9 Secretion of C4, and specificity of the monoclonal anti-C4 antibody for C4 protein in human brain tissue and cultured primary cortical neurons.

a, Brain tissue (from an individual affected with schizophrenia) was stained with a fluorescent secondary antibody, C4 antibody, or C4 antibody that was pre-adsorbed with purified C4 protein. Confocal images demonstrate the loss of immunoreactivity in the secondary-only and pre-adsorbed conditions. b, Primary human neurons were stained with a fluorescent secondary antibody, C4 antibody or C4 antibody that was pre-adsorbed with purified C4 protein. Confocal images demonstrate the loss of immunoreactivity in the secondary-only and pre-adsorbed conditions. Scale bars, 25 μm. c, Secretion of C4 protein by cultured primary neurons. Western blot for C4 protein analysis. (+) Purified human C4 protein. (–) Unconditioned medium, a negative control. HN-conditioned shows the same medium after conditioning by cultured human neurons at days 7 (d7) and 30 (d30). Details of western blot protocol, antibody catalogue numbers and concentrations used are in Supplementary Methods. C4 molecular weight, ~210 kDa.

Extended Data Figure 10 Mouse C4 genes and additional analyses of the dLGN eye segregation phenotype in C4 mutant mice and wild-type and heterozygous littermate controls.

a, The functional specialization of C4 into C4A and C4B in humans does not have an analogy in mice. Although the mouse genome contains both a C4 gene and a C4-like gene (classically called Slp), and these genes are also present as a tandem duplication within the mouse MHC locus, analysis of the encoded protein sequences indicates a distinct specialization, as illustrated by the protein phylogenetic tree. Top, mouse Slp is indicated in grey to reflect its potential pseudogenization: Slp is already known to have mutations at a C1s cleavage site, which are thought to abrogate activation of the protein through the classical complement pathway56; and the M. musculus reference genome sequence (mm10) at Slp shows a 1-bp deletion (relative to C4) within the coding region at chr17:34815158, which would be predicted to cause a premature termination of the encoded protein. In some genome data resources, mouse Slp and C4 have been annotated respectively as ‘C4a’ (for example, NM_011413.2) and ‘C4b’ (for example, NM_009780.2) based on synteny with the human C4A and C4B genes, but the above sequence analysis indicates that they are not paralogous to C4A and C4B. b, Sequence differences between C4A and C4B—which are otherwise 99.5% identical at an amino acid level—are concentrated at the ‘isotypic site’ where they shape each isotype’s relative affinity for different molecular targets19,20. At the isotypic site, mouse C4 contains a combination of the residues present in human C4A and C4B. c, Expression of mouse C4 mRNA in whole retina and lateral geniculate nucleus (LGN) from P5 animals and in purified retinal ganglion cells (RGCs) from P5 and P15 animals. These time points were chosen as P5 is a time of more robust synaptic refinement in the retinogeniculate system compared to P15. The same assays detected no C4 RNA in control RNA isolated from C4−/− mice (not shown). n = 3 samples for p5 retina, LGN, and P15 RGCs, n = 4 samples for P5 RGCs; *P < 0.05 by ANOVA with post-hoc Tukey–Kramer multiple-comparisons test. d, Representative images of dLGN innervation by contralateral projections (red in bottom image), ipsilateral projections (green in bottom image), and their overlap (yellow in bottom image). Scale bar, 100 μm. e, Quantification of the percentage of total dLGN area receiving both contralateral and ipsilateral projections shows a significant increase in C4−/− compared to wild-type littermates (ANOVA, n = 5 mice per group, P < 0.01). These data are consistent with results using R value analysis as shown in Fig. 7. f, Quantification of total dLGN area showed no significant difference between wild-type and C4−/− mice (ANOVA, n = 5 per group, P > 0.05). g, Quantification of dLGN area receiving ipsilateral innervation showed a significant increase in ipsilateral territory in the C4−/− mice compared to wild-type littermates (AVOVA, n = 5 mice per group, P > 0.01). This result is consistent with defects in eye specific segregation. Scale bar, 100 μm. h, The number of RGCs in the retina was estimated by counting the number of Brn3a+ cells in wild-type and C4−/− mice. No differences were observed between wild-type and C4−/− mice (t-test, n = 4 mice per group, P > 0.05). Scale bar, 100 μm.

Supplementary information

Supplementary Information

This file contains Supplementary Methods, Supplementary Tables 1-3, a full list of the collaborators from the PGC Schizophrenia Working Group and additional references. This file was replaced on 11 April 2016 to update affiliation 210. (PDF 668 kb)

PowerPoint slides

Rights and permissions

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Sekar, A., Bialas, A., de Rivera, H. et al. Schizophrenia risk from complex variation of complement component 4. Nature 530, 177–183 (2016). https://doi.org/10.1038/nature16549

Download citation

  • Received:

  • Accepted:

  • Published:

  • Issue Date:

  • DOI: https://doi.org/10.1038/nature16549

This article is cited by

Comments

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.

Search

Quick links

Nature Briefing

Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.

Get the most important science stories of the day, free in your inbox. Sign up for Nature Briefing