Varying combinations of genetic and environmental risk factors are thought to underpin phenotypic heterogeneity between individuals in psychiatric conditions such as schizophrenia. While epigenome-wide association studies in schizophrenia have identified extensive alteration of mean DNA methylation levels, less is known about the location and impact of DNA methylation variance, which could contribute to phenotypic and treatment response heterogeneity. To explore this question, we conducted the largest meta-analysis of blood DNA methylation variance in schizophrenia to date, leveraging three cohorts comprising 1036 individuals with schizophrenia and 954 non-psychiatric controls. Surprisingly, only a small proportion (0.1%) of the 213 variably methylated positions (VMPs) associated with schizophrenia (Benjamini-Hochberg FDR <
0.05) were shared with differentially methylated positions (DMPs
sites with mean changes between cases and controls). These blood-derived VMPs were found to be overrepresented in genes previously associated with schizophrenia and amongst brain-enriched genes, with evidence of concordant changes at VMPs in the cerebellum, hippocampus, prefrontal cortex, or striatum. Epigenetic covariance was also observed with respect to clinically significant metrics including age of onset, cognitive deficits, and symptom severity. We also uncovered a significant VMP in individuals with first-episode psychosis (n = 644) from additional cohorts and a non-psychiatric comparison group (n = 633). Collectively, these findings suggest schizophrenia is associated with significant changes in DNA methylation variance, which may contribute to individual-to-individual heterogeneity.