Retinoid metabolites of vitamin A are intrinsically linked to neural development, connectivity and plasticity, and have been implicated in the pathophysiology of schizophrenia. We hypothesised that a greater burden of common and rare genomic variation in genes involved with retinoid biogenesis and signalling could be associated with schizophrenia and its cognitive symptoms. Common variants associated with schizophrenia in the largest genome-wide association study were aggregated in retinoid genes and used to formulate a polygenic risk score (PRS Ret ) for each participant in the Australian Schizophrenia Research Bank. In support of our hypothesis, we found PRS Ret to be significantly associated with the disorder. Cases with severe cognitive deficits, while not further differentiated by PRS Ret , were enriched with rare variation in the retinoic acid receptor beta gene RARB, detected through whole-genome sequencing. RARB rare variant burden was also associated with reduced cerebellar volume in the cases with marked cognitive deficit, and with covariation in grey matter throughout the brain. An excess of rare variation was further observed in schizophrenia in retinoic acid response elements proximal to target genes, which we show are differentially expressed in the disorder in two RNA sequencing datasets. Our results suggest that genomic variation may disrupt retinoid signalling in schizophrenia, with particular significance for cases with severe cognitive impairment.

There is a long-standing interest in exploring the relationship between blood-based biomarkers and psychiatric disorders, despite their causal role being difficult to resolve in observational studies. In this study, we leverage genome-wide association study data for a large panel of heritable serum biochemical traits to refine our understanding of causal effect in biochemical-psychiatric trait pairings. We observed widespread positive and negative genetic correlation between psychiatric disorders and biochemical traits. Causal inference was then implemented to distinguish causation from correlation, with strong evidence that C-reactive protein (CRP) exerts a causal effect on psychiatric disorders. Notably, CRP demonstrated both protective and risk-increasing effects on different disorders. Multivariable models that conditioned CRP effects on interleukin-6 signaling and body mass index supported that the CRP-schizophrenia relationship was not driven by these factors. Collectively, these data suggest that there are shared pathways that influence both biochemical traits and psychiatric illness.

Brain morphology differs markedly between individuals with schizophrenia, but the cellular and genetic basis of this heterogeneity is poorly understood. Here, we sought to determine whether cortical thickness (CTh) heterogeneity in schizophrenia relates to interregional variation in distinct neural cell types, as inferred from established gene expression data and person-specific genomic variation. This study comprised 1849 participants in total, including a discovery (140 cases and 1267 controls) and a validation cohort (335 cases and 185 controls). To characterize CTh heterogeneity, normative ranges were established for 34 cortical regions and the extent of deviation from these ranges was measured for each individual with schizophrenia. CTh deviations were explained by interregional gene expression levels of five out of seven neural cell types examined: (1) astrocytes; (2) endothelial cells; (3) oligodendrocyte progenitor cells (OPCs); (4) excitatory neurons; and (5) inhibitory neurons. Regional alignment between CTh alterations with cell type transcriptional maps distinguished broad patient subtypes, which were validated against genomic data drawn from the same individuals. In a predominantly neuronal/endothelial subtype (22% of patients), CTh deviations covaried with polygenic risk for schizophrenia (sczPRS) calculated specifically from genes marking neuronal and endothelial cells (r = −0.40, p = 0.010). Whereas, in a predominantly glia/OPC subtype (43% of patients), CTh deviations covaried with sczPRS calculated from glia and OPC-linked genes (r = −0.30, p = 0.028). This multi-scale analysis of genomic, transcriptomic, and brain phenotypic data may indicate that CTh heterogeneity in schizophrenia relates to inter-individual variation in cell-type specific functions. Decomposing heterogeneity in relation to cortical cell types enables prioritization of schizophrenia subsets for future disease modeling efforts.