A Regional Burden of Sequence-Level Variation in the 22q11.2 Region Influences Schizophrenia Risk and Educational Attainment

Background: Genomic loci where recurrent pathogenic copy number variants are associated with psychiatric phenotypes in the population may also be sensitive to the collective impact of multiple functional low-frequency single nucleotide variants (SNVs). Methods: We examined the cumulative impact of low-frequency, functional SNVs within the 22q11.2 region on schizophrenia risk in a discovery cohort and an independent replication cohort (N = 1933 and N = 11,128, respectively), as well as the impact on educational attainment (EA) in a third, independent, general population cohort (N = 2081). In the discovery and EA cohorts, SNVs were identified using genotyping arrays; in the replication cohort, whole-exome sequencing was available. For verification, we compared the regional SNV count for schizophrenia cases in the discovery cohort with a normative count distribution derived from a large population dataset (N = 26,500) using bootstrap procedures. Results: In both schizophrenia cohorts, an increased regional SNV burden (≥4 low-frequency SNVs) in the 22q11.2 region was associated with schizophrenia (discovery cohort: odds ratio = 7.48, p = .039; replication cohort: odds ratio = 1.92, p = .004). In the EA cohort, an increased regional SNV burden at 22q11.2 was associated with decreased EA (odds ratio = 4.65, p = .049). Comparing the SNV count for schizophrenia cases with a normative distribution confirmed the unique nature of the distribution for schizophrenia cases (p = .002). Conclusions: In the general population, an increased burden of low-frequency, functional SNVs in the 22q11.2 region is associated with schizophrenia risk and a decrease in EA. These findings suggest that in addition to structural variation, a cumulative regional burden of low-frequency, functional SNVs in the 22q11.2 region can also have a relevant phenotypic impact.